Avian metapneumovirus: A five-plex digital droplet RT-PCR method for identification of subgroups A, B, C, and D.

End-point and real-time avian metapneumovirus (AMPV) RT-PCRs have been developed to detect one or two of the four recognized subgroups (A,B,C, and D) simultaneously or for broad range AMPV detection. Current subgroup specific tests target variable areas of the genome which makes these PCRs sensitive to specificity defects as recently documented. In the current study, a single five-plex digital droplet RT-PCR targeting the conserved viral polymerase gene of AMPV, which is less prone to genetic drift, has been designed. This digital droplet RT-PCR was capable of identifying each of the four AMPV subgroups. Each subgroup was identified according to a specifically assigned fluorescent amplitude. Specificity, which was tested including 31 AMPV strains, non-AMPV avian viruses and closely related human respiratory viruses, was 100%. The specific limit of detection for extracted viral RNA was estimated between 1 and 3 copies/µl. This tool simplifies the number of tests required for AMPV genotype diagnostics and should be theoretically less effected by viral genome evolution due to its target region. Ultimately, application of this test will contribute to an improved understanding of the global geographic distribution and subgroup host range of field strains.

Rabbit haemorrhagic disease virus Lagovirus europaeus/GI.1d strain: genome sequencing, in vivo virus replication kinetics, and viral dose effect.

BACKGROUND: Rabbit haemorrhagic disease virus Lagovirus europaeus/GI.1d variant (GI.1d/RHDV) was identified in 1990 in France, and until the emergence of the new genotype GI.2, it was the main variant circulating in the country. The early stages of RHDV infection have been described in a few studies of rabbits experimentally infected with earlier strains, but no information was given on the minimum infective dose. We report the genomic and phenotypic characterisation of a GI.1d/RHDV strain collected in 2000 in France (GI.1d/00-21). RESULTS: We performed in vivo assays in rabbits to study virus replication kinetics in several tissues at the early stage of infection, and to estimate the minimum infective dose. Four tested doses, negligible (10- 1 viral genome copies), low (104), high (107) and very high (1011) were quantified using a method combining density gradient centrifugation of the viral particles and an RT-qPCR technique developed to quantify genomic RNA (gRNA). The GI.1d/00-21 genome showed the same genomic organisation as other lagoviruses; however, a substitution in the 5' untranslated region and a change in the potential p23/2C-like helicase cleavage site were observed. We showed that the liver of one of the two rabbits inoculated via the oral route was infected at 16 h post-infection and all tissues at 39 h post-infection. GI.1d/00-21 induced classical RHD signs (depression) and lesions (haemorrhage and splenomegaly). Although infective dose estimation should be interpreted with caution, the minimum infective dose that infected an inoculated rabbit was lower or equal to 104 gRNA copies, whereas between 104 and 107 gRNA copies were required to also induce mortality. CONCLUSIONS: These results provide a better understanding of GI.1d/RHDV infection in rabbits. The genome analysis showed a newly observed mutation in the 5' untranslated region of a lagovirus, whose role remains unknown. The phenotypic analysis showed that the pathogenicity of GI.1d/00-21 and the replication kinetics in infected organs were close to those reported for the original GI.1 strains, and could not alone explain the observed selective advantage of the GI.1d strains. Determining the minimum dose of viral particles required to cause mortality in rabbits is an important input for in vivo studies.

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.

Genetic diversity and evolution of Hare Calicivirus (HaCV), a recently identified lagovirus from Lepus europaeus.

First recognized as highly pathogenic viruses, hare lagoviruses belonging to genotype GII.1 (EBHSV) infect various Lepus species. Genetically distinct benign lagoviruses (Hare Calicivirus, HaCV) have recently been identified but few data have been available so far on these strains. The analysis of 199 samples from hunted hares collected throughout France allowed the detection of 20 HaCV and showed that they were widely distributed in this country. Ten HaCV capsid protein gene sequences were characterized. A first HaCV capsid protein structural model was proposed, revealing a global structure similar to that of a pathogenic GII.1 strain. The HaCV sequences showed an even higher genetic diversity than previously appreciated, with the characterization of two genotypes (GII.2, GII.3) and several additional putative genotypes. The most recent common ancestor for HaCV VP60 gene was estimated to be much older than that for GII.1 pathogenic strains. These results give new insights into the phylogenetic relationships of HaCV within the Lagovirus genus.

Host specificity of avian metapneumoviruses.

To date, four subgroups of avian metapneumoviruses have been defined (AMPV-A, B, C and D) based on genetic and antigenic differences. The extent of infection in the three principal species (turkeys, chickens and ducks) by these subgroups is, however, not well defined. Here, a series of controlled and ethically approved experimental infections were performed in specific pathogen-free turkeys, chickens and ducks with each of the four AMPV subgroups. For subgroup C, one strain isolated from turkeys in the USA (turkey AMPV-C) and one isolated from ducks in France (duck AMPV-C) were compared. Globally, these extensive experimental trials demonstrated that AMPV-A, B, turkey C and D were well adapted to Galliformes, especially turkeys; however, chickens showed limited clinical signs and differences in seroconversion and transmission. Notably, chickens did not transmit AMPV-A to contacts and were shown for the first time to be susceptible to AMPV-D. The duck AMPV-C was well adapted to ducks; however, chickens and turkeys seroconverted and were positive by virus isolation. In addition, seroconversion of contact turkeys to duck AMPV-C demonstrated horizontal transmission of this virus in a non-palmiped species under our experimental conditions. Interestingly, in chickens and turkeys, duck AMPV-C isolation was possible despite a lack of detection of viral RNA. Likewise, the turkey AMPV-C virus was well adapted to turkeys yet was also isolated from chickens despite a lack of detection of viral RNA. These results would suggest a selection for viral genetic sequences that differ from the original strain upon adaptation to a 'non-conventional host'.

First Complete Genome Sequence of a Hare Calicivirus Strain Isolated from Lepus europaeus.

The first full-genome sequence of a hare calicivirus (HaCV), recently characterized as a novel member of the Caliciviridae, is described. This presumed nonpathogenic lagovirus is 7,433 nucleotides long, shows the same genomic organization as that of other lagoviruses, and has the highest nucleotide identity (79%) with pathogenic European brown hare syndrome viruses.

First complete genome sequence of a European non-pathogenic rabbit calicivirus (lagovirus GI.3).

We report the full genome sequence of the non-pathogenic rabbit lagovirus Lagovirus europaeus/GI.3/O cun/FR/2006/06-11 (GI.3/06-11), collected from a healthy French domestic rabbit in 2006, and initially described as 06-11 strain. The sequence reveals a genomic organization similar to lagoviruses. It was 7,436 bases long and contained two open reading frames (ORF). A dipeptide variation at the potential p23/2C-like helicase cleavage site (EE instead of ED) was observed, a feature only shared with non-recombinant pathogenic lagoviruses in GI.2 and with two European brown hare syndrome viruses (EBHSV) collected in 1982 in Sweden. GI.3/06-11 has only one initiation codon at the beginning of the ORF2 like the avirulent Italian rabbit calicivirus (RCV) and EBHSV. Previous genetic analyses based on the capsid gene sequences showed that GI.3/06-11 was closer to the RCV and pathogenic lagoviruses GI.1 strains than other lagoviruses. This study, by revealing that GI.3/06-11 genome sequence significantly clustered with pathogenic GI.2 strains, gives prominence of new genetic relationship among lagoviruses and should contribute to understand the emergence of pathogenic strains.

Large-scale lagovirus disease outbreaks in European brown hares (Lepus europaeus) in France caused by RHDV2 strains spatially shared with rabbits (Oryctolagus cuniculus).

Rabbit haemorrhagic disease virus (RHDV) is a lagovirus that causes rabbit haemorrhagic disease (RHD) in European rabbits (Oryctolagus cuniculus). In 2010, a new genotype called RHDV2 emerged in France. It exhibits a larger host range than classical RHDV strains by sporadically infecting different hare species, including the European hare (Lepus europaeus). Phylogenetic analyses revealed that closely related RHDV2 strains circulate locally in both hares and rabbits, and therefore that RHDV2 strains infecting hares do not belong to a lineage that has evolved only in this species. We showed that RHDV2 is widely distributed in France and that it was responsible for more than a third of cases of lagovirus disease in European hare populations in 2015. The oldest RHDV2 positive hare was sampled in November 2013 and we reported two hares co-infected by EBHSV and RHDV2. All together, our results raise important epidemiological and evolutionary issues. In particular, along with the potential emergence of recombinant EBHSV/RHDV2 strains in hares, the enlargement of the host range changes the host population structure of RHDV2 and may alter the impact of the virus on rabbit and hare populations.

Emerging highly pathogenic H5 avian influenza viruses in France during 2015/16: phylogenetic analyses and markers for zoonotic potential.

Several new highly pathogenic (HP) H5 avian influenza virus (AIV) have been detected in poultry farms from south-western France since November 2015, among which an HP H5N1. The zoonotic potential and origin of these AIVs immediately became matters of concern. One virus of each subtype H5N1 (150169a), H5N2 (150233) and H5N9 (150236) was characterised. All proved highly pathogenic for poultry as demonstrated molecularly by the presence of a polybasic cleavage site in their HA protein - with a sequence (HQRRKR/GLF) previously unknown among avian H5 HPAI viruses - or experimentally by the in vivo demonstration of an intravenous pathogenicity index of 2.9 for the H5N1 HP isolate. Phylogenetic analyses based on the full genomes obtained by NGS confirmed that the eight viral segments of the three isolates were all part of avian Eurasian phylogenetic lineage but differed from the Gs/Gd/1/96-like lineage. The study of the genetic characteristics at specific amino acid positions relevant for modulating the adaptation to and the virulence for mammals showed that presently, these viruses possess most molecular features characteristic of AIV and lack some major characteristics required for efficient respiratory transmission to or between humans. The three isolates are therefore predicted to have no significant pandemic potential.

Molecular comparisons of full length metapneumovirus (MPV) genomes, including newly determined French AMPV-C and -D isolates, further supports possible subclassification within the MPV Genus.

Four avian metapneumovirus (AMPV) subgroups (A-D) have been reported previously based on genetic and antigenic differences. However, until now full length sequences of the only known isolates of European subgroup C and subgroup D viruses (duck and turkey origin, respectively) have been unavailable. These full length sequences were determined and compared with other full length AMPV and human metapneumoviruses (HMPV) sequences reported previously, using phylogenetics, comparisons of nucleic and amino acid sequences and study of codon usage bias. Results confirmed that subgroup C viruses were more closely related to HMPV than they were to the other AMPV subgroups in the study. This was consistent with previous findings using partial genome sequences. Closer relationships between AMPV-A, B and D were also evident throughout the majority of results. Three metapneumovirus "clusters" HMPV, AMPV-C and AMPV-A, B and D were further supported by codon bias and phylogenetics. The data presented here together with those of previous studies describing antigenic relationships also between AMPV-A, B and D and between AMPV-C and HMPV may call for a subclassification of metapneumoviruses similar to that used for avian paramyxoviruses, grouping AMPV-A, B and D as type I metapneumoviruses and AMPV-C and HMPV as type II.

Emergence of a new lagovirus related to Rabbit Haemorrhagic Disease Virus.

Since summer 2010, numerous cases of Rabbit Haemorrhagic Disease (RHD) have been reported in north-western France both in rabbitries, affecting RHD-vaccinated rabbits, and in wild populations. We demonstrate that the aetiological agent was a lagovirus phylogenetically distinct from other lagoviruses and which presents a unique antigenic profile. Experimental results show that the disease differs from RHD in terms of disease duration, mortality rates, higher occurrence of subacute/chronic forms and that partial cross-protection occurs between RHDV and the new RHDV variant, designated RHDV2. These data support the hypothesis that RHDV2 is a new member of the Lagovirus genus. A molecular epidemiology study detected RHDV2 in France a few months before the first recorded cases and revealed that one year after its discovery it had spread throughout the country and had almost replaced RHDV strains. RHDV2 was detected in continental Italy in June 2011, then four months later in Sardinia.

An alternative method to determine the 5' extremities of non-segmented, negative sense RNA viral genomes using positive replication intermediate 3' tailing: application to two members of the Paramyxoviridae family.

Determining the sequence of non-segmented, negative sense RNA viral genomes is far from routine and often requires the application of several techniques. In this study, an existing method used currently just for determination of the genomic 3' extremity was used to determine both the 3' and 5' sequence extremities of a Newcastle disease virus and an avian metapneumovirus. This was achieved with a single 3' nucleotide tailing reaction of both the genomic RNA and the full length, positive sense, antigenomic RNA, followed by a single reverse transcription reaction targeted to the common polynucleotide tails, and then individual PCRs specific for each extremity using PCR primers derived from the sequence of the RT primer or from neighbouring virus sequences known previously. For each virus the method was employed separately. Sequences from both viruses were in agreement with those reported previously for other paramyxoviruses, yet one extra base at the 3' and one extra base at the 5' were identified for the avian metapneumovirus. In this study, importantly, the newly determined extremities maintained the complementarity known to exist between the extremities of these viruses. The method was equally successful with both viruses and can be tailored easily to function with other non-segmented, negative sense viruses through minor modification of only the primer sequences.