Evolutionary and functional characterization of lagomorph guanylate-binding proteins: a story of gain and loss and shedding light on expression, localization and innate immunity-related functions.

Guanylate binding proteins (GBPs) are an evolutionarily ancient family of proteins that are widely distributed among eukaryotes. They belong to the dynamin superfamily of GTPases, and their expression can be partially induced by interferons (IFNs). GBPs are involved in the cell-autonomous innate immune response against bacterial, parasitic and viral infections. Evolutionary studies have shown that GBPs exhibit a pattern of gene gain and loss events, indicative for the birth-and-death model of evolution. Most species harbor large GBP gene clusters that encode multiple paralogs. Previous functional and in-depth evolutionary studies have mainly focused on murine and human GBPs. Since rabbits are another important model system for studying human diseases, we focus here on lagomorphs to broaden our understanding of the multifunctional GBP protein family by conducting evolutionary analyses and performing a molecular and functional characterization of rabbit GBPs. We observed that lagomorphs lack GBP3, 6 and 7. Furthermore, Leporidae experienced a loss of GBP2, a unique duplication of GBP5 and a massive expansion of GBP4. Gene expression analysis by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) and transcriptome data revealed that leporid GBP expression varied across tissues. Overexpressed rabbit GBPs localized either uniformly and/or discretely to the cytoplasm and/or to the nucleus. Oryctolagus cuniculus (oc)GBP5L1 and rarely ocGBP5L2 were an exception, colocalizing with the trans-Golgi network (TGN). In addition, four ocGBPs were IFN-inducible and only ocGBP5L2 inhibited furin activity. In conclusion, from an evolutionary perspective, lagomorph GBPs experienced multiple gain and loss events, and the molecular and functional characteristics of ocGBP suggest a role in innate immunity.

A new HaCV-EBHSV recombinant lagovirus circulating in European brown hares (Lepus europaeus) from Catalonia, Spain.

In 2020/2021, several European brown hare syndrome virus (EBHSV) outbreaks were recorded in European hares (Lepus europaeus) from Catalonia, Spain. Recombination analysis combined with phylogenetic reconstruction and estimation of genetic distances of the complete coding sequences revealed that 5 strains were recombinants. The recombination breakpoint is located within the non-structural protein 2C-like RNA helicase (nucleotide position ~ 1889). For the genomic fragment upstream of the breakpoint, a non-pathogenic EBHSV-related strain (hare calicivirus, HaCV; GII.2) was the most closely related sequence; for the rest of the genome, the most similar strains were the European brown hare syndrome virus (EBHSV) strains recovered from the same 2020/2021 outbreaks, suggesting a recent origin. While the functional impact of the atypical recombination breakpoint remains undetermined, the novel recombinant strain was detected in different European brown hare populations from Catalonia, located 20-100 km apart, and seems to have caused a fatal disease both in juvenile and adult animals, confirming its viability and ability to spread and establish infection. This is the first report of a recombination event involving HaCV and EBHSV and, despite the recombination with a non-pathogenic strain, it appears to be associated with mortality in European brown hares, which warrants close monitoring.

On the virome's hidden diversity: lessons from RHDV.

Emerging infectious diseases are a major challenge to human and animal health. While predicting the emergence of pathogens is complex, the advent of high-throughput sequencing technologies has allowed the rapid identification of unknown microbiology diversity within organisms. Here, we discuss an example of a metatranscriptomics output to decipher viral evolution.

First detection and molecular characterization of rabbit hemorrhagic disease virus (RHDV) in Algeria.

Since the first detection of rabbit hemorrhagic disease (RHD), the rabbit hemorrhagic disease virus (RHDV) has been responsible for high morbidity and mortality worldwide, both in domestic and in wild rabbits. Despite the apparent control of RHD in rabbitries through vaccination, several studies highlighted the rapid evolution of RHDV by recombination, which may facilitate the emergence of new pathogenic strains. The aim of this study was to confirm the presence and characterize RHDV in Algeria. For this, rabbit samples were collected in the north of Algeria, between 2018 and 2021, from small farms where the virus was suspected after the sudden death of a high number of rabbits, and from healthy hunted wild rabbits. The domestic rabbits revealed clinical signs and lesions that were suggestive of RHD. RT-PCR showed that 79.31% of the domestic rabbit samples were positive for RHDV, while in 20.69%, including the hunted rabbits, the virus was not detected. Phylogenetic analysis of the Algerian strains allowed the confirmation and identification as GI.2 (RHDV2), and showed a close relation to GI.3P-GI.2 recombinant strains, suggesting a potential introduction from other countries, with an older strain potentially originated from neighboring Tunisia, while more recent isolates grouped with strains from North America. Our study reports for the first time the presence of GI.2 (RHDV2) in Algeria with multiple routes of introduction. Consequently, we propose that RHDV control in Algeria should be based on epidemiological surveys in association with an adequate prophylactic program.

Evolution of primate interferon-induced transmembrane proteins (IFITMs): a story of gain and loss with a differentiation into a canonical cluster and IFITM retrogenes.

Interferon-inducible transmembrane proteins (IFITMs) are a family of transmembrane proteins. The subgroup of immunity-related (IR-)IFITMs is involved in adaptive and innate immune responses, being especially active against viruses. Here, we suggest that IFITMs should be classified as (1) a canonical IFITM gene cluster, which is located on the same chromosome, and (2) IFITM retrogenes, with a random and unique location at different positions within the genome. Phylogenetic analyses of the canonical cluster revealed the existence of three novel groups of primate IFITMs (pIFITM) in the IR-IFITM clade: the prosimian pIFITMs(pro), the new world monkey pIFITMs(nwm) and the old world monkey pIFITMs(owm). Therefore, we propose a new nomenclature: IR-pIFITM1, IR-pIFITM2, IR-pIFITM3, IR-pIFITMnwm, IR-pIFITMowm, and IR-pIFITMpro. We observed divergent evolution for pIFITM5 and pIFITM10, and evidence for concerted evolution and a mechanism of birth-and-death evolution model for the IR-pIFITMs. In contrast, the IFITMs scattered throughout the genomes possessed features of retrogenes retrotransposed by class 1 transposable elements. The origin of the IFITM retrogenes correspond to more recent events. We hypothesize that the transcript of a canonical IFITM3 has been constantly retrotransposed using class 1 transposable elements resulting in the IFITM retro(pseudo)genes. The unique pattern of each species has most likely been caused by constant pseudogenization and loss of the retro(pseudo)genes. This suggests a third mechanism of evolution for the IR-IFITMs in primates, similar to the birth-and-death model of evolution, but via a transposable element mechanism, which resulted in retro(pseudo)genes.

Identification and characterization of a polyomavirus in the thornback skate (Raja clavata).

Members of the family Polyomaviridae have a circular double-stranded DNA genome that have been identified in various hosts ranging from mammals to arachnids. Here we report the identification and analysis of a complete genome sequence of a novel polyomavirus, Raja clavata polyomavirus (RcPyV1), from a cartilaginous fish, the thornback skate (Raja clavata). The genome sequence was determined using a metagenomics approach with an aim to provide baseline viral data in cartilaginous fish in different ecosystems. The RcPyV1 genome (4,195 nucleotides) had typical organization of polyomavirus, including early antigens (small T; Large T) encoded on one strand and late viral proteins (VP1; VP2) on the complementary strand. Maximum-likelihood phylogenetic analysis of the large T-antigen revealed that RcPyV1 clusters with a polyomavirus obtained from another cartilaginous fish, the guitarfish polyomavirus 1 (GfPyV1). These two share ~ 56% pairwise identity in LT and VP1 protein sequences. These analyses support the hypothesis that cartilaginous fishes have a specific lineage of polyomaviruses.

The potential role of scavenging flies as mechanical vectors of Lagovirus europaeus/GI.2.

The European rabbit (Oryctolagus cuniculus) populations of the Iberian Peninsula have been severely affected by the emergence of the rabbit haemorrhagic disease virus (RHDV) Lagovirus europaeus/GI.2 (RHDV2/b). Bushflies and blowflies (Muscidae and Calliphoridae families, respectively) are important RHDV vectors in Oceania, but their epidemiological role is unknown in the native range of the European rabbit. In this study, scavenging flies were collected between June 2018 and February 2019 in baited traps at one site in southern Portugal, alongside a longitudinal capture-mark-recapture study of a wild European rabbit population, aiming to provide evidence of mechanical transmission of GI.2 by flies. Fly abundance, particularly from Calliphoridae and Muscidae families, peaked in October 2018 and in February 2019. By employing molecular tools, we were able to detect the presence of GI.2 in flies belonging to the families Calliphoridae, Muscidae, Fanniidae and Drosophilidae. The positive samples were detected during an RHD outbreak and absent in samples collected when no evidence of viral circulation in the local rabbit population was found. We were able to sequence a short viral genomic fragment, confirming its identity as RHDV GI.2. The results suggest that scavenging flies may act as mechanical vectors of GI.2 in the native range of the southwestern Iberian subspecies O. cuniculus algirus. Future studies should better assess their potential in the epidemiology of RHD and as a tool for monitoring viral circulation in the field.

Case series: Four fatal rabbit hemorrhagic disease virus infections in urban pet rabbits.

Four pet rabbits (Oryctolagus cuniculus cuniculus) diagnosed with a fatal infection by rabbit hemorrhagic disease virus (RHDV GI.2) were identified in the same week and further investigated. All animals lived in an urban environment (Lisbon, Portugal), were between 8 months and 2 years old and none had been vaccinated against RHDV2 (GI.2). Three animals arrived at the clinic and died shortly afterward and it was only possible to collect material for RT-qPCR (RHDV) test. These rabbits tested positive for RHDV2, with high viral loads. In the fourth case, additional clinical and post-mortem gross and histological evaluations were performed. This 8 month old intact female indoor pet rabbit was presented with apathy, tachypnea and tachycardia. Radiographic projections revealed no clinical revealed no clinical abnormalities. Serum biochemistry revealed a significant increase in AST and ALT with a small hypoglycemia. Abdominal ultrasound revealed an acute hepatitis. Despite hospitalization support, after 30 h of admission, the rabbit lost consciousness and developed anorexia and pyrexia in the last minutes before death. Post-mortem analysis and molecular testing by RT-qPCR, confirmed the diagnosis of RHDV2 (GI.2) infection also with high viral load. In conclusion, this paper reports a case series that demonstrates the severe infectious ability and the high mortality associated with RHDV even in rabbits from urban environments. Furthermore, it highlights the importance of always considering rabbit hemorrhagic disease (RHD) as a differential diagnosis in pet rabbits with non-specific clinical signs, and should warn veterinarians that pet rabbits living indoors can also be infected with a fatal outcome.

Functional cross-species conservation of guanylate-binding proteins in innate immunity.

Guanylate binding proteins (GBPs) represent an evolutionary ancient protein family widely distributed among eukaryotes. They are interferon (IFN)-inducible guanosine triphosphatases that belong to the dynamin superfamily. GBPs are known to have a major role in the cell-autonomous innate immune response against bacterial, parasitic and viral infections and are also involved in inflammasome activation. Evolutionary studies depicted that GBPs present a pattern of gain and loss of genes in each family with several genes pseudogenized and some genes more divergent, indicative for the birth-and-death evolution process. Most species harbor large GBP gene clusters encoding multiple paralogs. Previous functional studies mainly focused on mouse and human GBPs, but more data are becoming available, broadening the understanding of this multifunctional protein family. In this review, we will provide new insights and give a broad overview about GBP evolution, conservation and their roles in all studied species, including plants, invertebrates and vertebrates, revealing how far the described features of GBPs can be transferred to other species.

Multi-event capture-recapture models estimate the diagnostic performance of serological tests for myxoma and rabbit haemorrhagic disease viruses in the absence of reference samples.

Estimation of the diagnostic performance of serological tests often relies on another test assumed as a reference or on samples of known infection status, yet both are seldom available for emerging pathogens in wildlife. Longitudinal disease serological data can be analysed through multi-event capture-mark-recapture (MECMR) models accounting for the uncertainty in state assignment, allowing us to estimate epidemiological parameters such as incidence and mortality. We hypothesized that by estimating the uncertainty in state assignment, MECMR models estimate the diagnostic performance of serological tests for rabbit haemorrhagic disease virus (RHDV) and myxoma virus (MYXV). We evaluated this hypothesis on longitudinal serological data of three tests of RHDV and one test of MYXV in two populations of the European rabbit (Oryctolagus cuniculus algirus). First, we selected the optimal cut-off threshold for each test using finite mixture models, a reference method not relying on reference tests or samples. Second, we used MECMR models to compare the diagnostic sensitivity (Se) and specificity (Sp) of the three tests for RHDV. Third, we compared the estimates of diagnostic performance by MECMR and finite mixture models across a range of cut-off values. The MECMR models showed that the RHDV test employing GI.2 antigens (Se: 100%) outperformed two tests employing GI.1 antigens (Se: 21.7% ± 8.6% and 8.7% ± 5.9%). At their selected cut-offs (2.0 for RHDV GI.2 and 2.4 for MYXV), the estimates of Se and Sp were concordant between the MECMR and finite mixture models. Over the duration of the study (May 2018 to September 2020), the monthly survival of European rabbits seropositive for MYXV was significantly higher than that of seronegative rabbits (82.7% ± 4.9% versus 61.5% ± 12.7%) at the non-fenced site. We conclude that MECMR models can reliably estimate the diagnostic performance of serological tests for RHDV and MYXV in European rabbits. This conclusion could extend to other diagnostic tests and host-pathogen systems. Longitudinal disease surveillance data analysed through MECMR models allow the validation of diagnostic tests for emerging pathogens in novel host species while simultaneously estimating epidemiological parameters.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages.

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

Evolution of Guanylate Binding Protein (GBP) Genes in Muroid Rodents (Muridae and Cricetidae) Reveals an Outstanding Pattern of Gain and Loss.

Guanylate binding proteins (GBPs) are paramount in the host immunity by providing defense against invading pathogens. Multigene families related to the immune system usually show that the duplicated genes can either undergo deletion, gain new functions, or become non-functional. Here, we show that in muroids, the Gbp genes followed an unusual pattern of gain and loss of genes. Muroids present a high diversity and plasticity regarding Gbp synteny, with most species presenting two Gbp gene clusters. The phylogenetic analyses revealed seven different Gbps groups. Three of them clustered with GBP2, GBP5 and GBP6 of primates. Four new Gbp genes that appear to be exclusive to muroids were identified as Gbpa, b, c and d. A duplication event occurred in the Gbpa group in the common ancestor of Muridae and Cricetidae (~20 Mya), but both copies were deleted from the genome of Mus musculus, M. caroli and Cricetulus griseus. The Gbpb gene emerged in the ancestor of Muridae and Cricetidae and evolved independently originating Gbpb1 in Muridae, Gbpb2 and Gbpb3 in Cricetidae. Since Gbpc appears only in three species, we hypothesize that it was present in the common ancestor and deleted from most muroid genomes. The second Gbp gene cluster, Gbp6, is widespread across all muroids, indicating that this cluster emerged before the Muridae and Cricetidae radiation. An expansion of Gbp6 occurred in M. musculus and M. caroli probably to compensate the loss of Gbpa and b. Gbpd is divided in three groups and is present in most muroids suggesting that a duplication event occurred in the common ancestor of Muridae and Cricetidae. However, in Grammomys surdaster and Mus caroli, Gbpd2 is absent, and in Arvicanthis niloticus, Gbpd1 appears to have been deleted. Our results further demonstrated that primate GBP1, GBP3 and GBP7 are absent from the genome of muroids and showed that the Gbp gene annotations in muroids were incorrect. We propose a new classification based on the phylogenetic analyses and the divergence between the groups. Extrapolations to humans based on functional studies of muroid Gbps should be re-evaluated. The evolutionary analyses of muroid Gbp genes provided new insights about the evolution and function of these genes.

Sleep-wake patterns and disturbances in Portuguese primary school children: a comparison between 1995 and 2016.

Purpose: It is often assumed sleep duration has decreased and sleep schedules have delayed over the last decades, as society modernized. We aimed to investigate changes in the sleep patterns of school-age children over time. Methods: We compared the sleep timings, durations, and disturbances of primary school-age children in 1995 and roughly two decades later, in 2016. Data from 666 children attending the 3rd and 4th grades of basic education were combined from two different cross-sectional school-based studies conducted within the same educational region of mainland Portugal using the same parent-report questionnaire (Children's Sleep-wake Patterns Questionnaire). Results: Mean sleep duration did not differ significantly between the two time points (schooldays: t = .118, p = .906; free days: t = 1.310, p = .191), albeit the percentage of children sleeping the recommended number of hours decreased significantly in 2016 when compared to 1995 (schooldays: χ2 = 4.406, p = .036; free days: χ2 = 16.859, p < .001). Wake-times advanced on free days in 2016. Difficulties on settling to sleep alone and returning to sleep were more prevalent in 2016, as well as fearing the dark and needing lights on or parent's presence to fall asleep. Conclusions: Sleep onset-related disturbances appear to have increased from 1995 to 2016. One possible explanation for this increase might be the change in parental practices preventing children from learning to fall asleep autonomously.

Multiple Introductions of Rabbit Hemorrhagic Disease Virus Lagovirus europaeus/GI.2 in Africa.

Rabbit hemorrhagic disease (RHD) causes high mortality and morbidity in European rabbits (Oryctolagus cuniculus). In Africa, the presence of the causative agent, the rabbit hemorrhagic disease virus (RHDV), was first confirmed in 1992 (genotype Lagovirus europaeus/GI.1). In 2015, the new genotype Lagovirus europaeus/GI.2 (RHDV2/b) was detected in Tunisia. Currently, GI.2 strains are present in several North and Sub-Saharan African countries. Considerable economic losses have been observed in industrial and traditional African rabbitries due to RHDV. Like other RNA viruses, this virus presents high recombination rates, with the emergence of GI.2 being associated with a recombinant strain. Recombination events have been detected with both pathogenic (GI.1b and GII.1) and benign (GI.3 and GI.4) strains. We obtained complete genome sequences of Tunisian GI.2 strains collected between 2018 and 2020 and carried out phylogenetic analyses. The results revealed that Tunisian strains are GI.3P-GI.2 strains that were most likely introduced from Europe. In addition, the results support the occurrence of multiple introductions of GI.2 into Africa, stressing the need for characterizing complete genome sequences of the circulating lagoviruses to uncover their origin. Continued monitoring and control of rabbit trade will grant a better containment of the disease and reduce the disease-associated economic losses.

Spillover event of recombinant Lagovirus europaeus/GI.2 into the Iberian hare (Lepus granatensis) in Spain.

Viruses that affect lagomorphs have decades of reported history of spillover events. One of these viruses is the causative agent of the so-called rabbit or 'lagomorph' haemorrhagic disease (e.g. Lagovirus europaeus/GI.1 and L. europaeus/GI.2). In particular, L. europaeus/GI.2 has shown a great capacity to recombine with existing lagoviruses. In fact, it has replaced the former GI.1 genotype in the wild, and recently, an increase on spillover events has been detected among several lagomorph species including European and North American species of hares. In this study, we report for the first time the infection of a wild Iberian hare with GI.2 (RHDV2/b), potential shedding and associated histopathological alterations. We identify the recombinant GI.4P-GI.2 as causative of the infection and discuss plausible causes regarding the origin of the spillover event and its potential consequences for the Iberian hare wild populations, which is an endemic species of the Iberian Peninsula as well as an important game and prey species for many predators, including endangered species.

A Review on the Methods Used for the Detection and Diagnosis of Rabbit Hemorrhagic Disease Virus (RHDV).

Since the early 1980s, the European rabbit (Oryctolagus cuniculus) has been threatened by the rabbit hemorrhagic disease (RHD). The disease is caused by a lagovirus of the family Caliciviridae, the rabbit hemorrhagic disease virus (RHDV). The need for detection, identification and further characterization of RHDV led to the development of several diagnostic tests. Owing to the lack of an appropriate cell culture system for in vitro propagation of the virus, much of the methods involved in these tests contributed to our current knowledge on RHD and RHDV and to the development of vaccines to contain the disease. Here, we provide a comprehensive review of the RHDV diagnostic tests used since the first RHD outbreak and that include molecular, histological and serological techniques, ranging from simpler tests initially used, such as the hemagglutination test, to the more recent and sophisticated high-throughput sequencing, along with an overview of their potential and their limitations.

Evolution of the guanylate binding protein (GBP) genes: Emergence of GBP7 genes in primates and further acquisition of a unique GBP3 gene in simians.

Guanylate binding proteins (GBPs) are major players in the host immunity, providing defense against bacterial and viral invaders. Multigene families may suffer different processes of evolution. Gene families related to the immune system usually follow the birth-and-death evolution process, where duplicated genes can be deleted, gain new functions or become non-functional. We analyzed publicly available primate GBP sequences and their genomic organization and observed that GBP7 genes appear to have emerged from a duplication of GBP4 and seem to be only present in primates. Furthermore, GBP3 genes are only present in Simiiformes and probably originated from GBP1 genes. Finally, a duplication event occurred in the GBP6 in Tarsiiformes and became functional which might also explain the duplication of GBP6 in New World monkeys and Cercopithecidae. Taken together, this study provides new knowledge on the evolution of GBPs in primates and suggests that a revision of the GBPs nomenclature is necessary.

Real-world intensive safety monitoring of biosimilars rituximab and trastuzumab in a Portuguese oncology hospital.

PURPOSE: The aim of this study was to assess the safety profiles of two biosimilar medicines (rituximab and trastuzumab) in the treatment of cancer patients within a Portuguese oncology hospital. METHODS: This hospital-based prospective observational study followed a cohort event monitoring approach focused on signalling suspected adverse drug reactions (ADRs). Patients undergoing treatment with rituximab biosimilar CT-P10 (Truxima®) or trastuzumab biosimilar CT-P6 (Herzuma®) were recruited over an 11-month and a 6-month period, respectively. Clinicians identified eligible patients and used paper-based forms to report all ADRs associated with biosimilar medicines. ADR case reports were assessed for seriousness, expectedness and causality in the Pharmacovigilance Unit of Coimbra. RESULTS: Ninety-four patients received biosimilar medicines (rituximab, n = 35; trastuzumab, n = 59). Of those, 4 patients (11.4%) experienced 16 ADRs with rituximab and 1 patient (1.7%) experienced 5 ADRs with trastuzumab. All case reports contained serious and expected ADRs that were at least probably related with biosimilar medicines under study. Based on the MedDRA PT coding, the most reported ADR for rituximab CT-P10 was chest discomfort (n = 4; 19.1%), followed by odynophagia (n = 2; 9.5%). Trastuzumab CT-P6 was associated with back pain, headache, pain in extremity, tachypnoea and tremor (each, n = 1; 4.8%). CONCLUSION: The results of this study suggest that using biosimilar rituximab and biosimilar trastuzumab to treat cancer patients in the real-world clinical setting is associated with acceptable safety profiles. No new safety problems were identified.

Recombination at the emergence of the pathogenic rabbit haemorrhagic disease virus Lagovirus europaeus/GI.2.

Rabbit haemorrhagic disease is a viral disease that emerged in the 1980s and causes high mortality and morbidity in the European rabbit (Oryctolagus cuniculus). In 2010, a new genotype of the rabbit haemorrhagic disease virus emerged and replaced the former circulating Lagovirus europaeus/GI.1 strains. Several recombination events have been reported for the new genotype Lagovirus europaeus/GI.2, with pathogenic (variants GI.1a and GI.1b) and benign (genotype GI.4) strains that served as donors for the non-structural part while GI.2 composed the structural part; another recombination event has also been described at the p16/p23 junction involving GI.4 strains. In this study, we analysed new complete coding sequences of four benign GI.3 strains and four GI.2 strains. Phylogenetic and recombination detection analyses revealed that the first GI.2 strains, considered as non-recombinant, resulted from a recombination event between GI.3 and GI.2, with GI.3 as the major donor for the non-structural part and GI.2 for the structural part. Our results indicate that recombination contributed to the emergence, persistence and dissemination of GI.2 as a pathogenic form and that all described GI.2 strains so far are the product of recombination. This highlights the need to study full-genomic sequences of lagoviruses to understand their emergence and evolution.

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.

Recombination is one of the major sources of genetic variation in viruses. RNA viruses, such as rabbit hemorrhagic disease virus (RHDV), are among the viruses with the highest recombination rates. Several recombination events have been described for RHDV, mostly as a consequence of their genomic architecture. Here, we undertook phylogenetic and recombination analyses of French and Swedish RHDV strains from 1994 to 2016 and uncovered a new intergenotypic recombination event. This event occurred in the late 1990s/early 2000s and involved nonpathogenic GI.3 strains as donors for the nonstructural part of the genome of these recombinants, while pathogenic GI.1d strains contributed to the structural part. These GI.3P-GI.1d recombinant strains did not entirely replace GI.1d (nonrecombinant) strains, but became the dominant strains in France and Sweden, likely due to a fitness advantage associated with this genomic architecture. GI.3P-GI.1d (P stands for polymerase) strains persisted until 2013 and 2016 in Sweden and France, respectively, and cocirculated with the new genotype GI.2 in France. Since strains from the first GI.2 outbreaks were GI.3P-GI.2, we hypothesize that GI.3P-GI.1d could be the parental strain. Our results confirm the outstanding recombination ability of RHDV and its importance in the evolution of lagoviruses, which was only revealed by studying complete genomic sequences.