Complete genome reconstruction of the global and European regional dispersal history of the lumpy skin disease virus.

IMPORTANCE: Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015-2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.

RNA sequencing of avian paramyxovirus (Paramyxoviridae, Avulavirinae) isolates from wild mallards in Belgium, 2021: complete genomes and coinfections.

We used untargeted RNA sequencing to characterize three Avulavirinae isolates from pooled samples obtained from wild mallards in Belgium in 2021. The complete genome sequences of two avian Orthoavulavirus-1 (AOAV-1) strains and one avian Paraavulavirus-4 (APMV-4) strain were determined confirming hemagglutination inhibition testing of the virus isolates. In addition, the applied sequencing strategy identified an avian influenza virus (AIV) coinfection in all three virus isolates, confirming weak-positive AIV realtime RT-PCR results from the original sample material. In one AOAV-1 isolate, partial sequences covering all genome segments of an AIV of subtype H11N9 could be de novo assembled from the sequencing data. Besides an AIV coinfection, RNA metagenomic data from the APMV-4 isolate also showed evidence of Alpharetrovirus and Megrivirus coinfection. In total, two AOAV-1 of Class II, genotype I.2 and one APMV-4 complete genome sequences were assembled and compared to publicly available sequences, highlighting the importance of surveillance for poultry pathogens in wild birds. Beyond the insights from full genome characterization of virus isolates, untargeted RNA sequencing strategies provide additional insights in the RNA virome of clinical samples as well as their derived virus isolates that are particularly useful when targeting wild avifauna reservoirs of poultry pathogens.

Combined Phylogeographic Analyses and Epidemiologic Contact Tracing to Characterize Atypically Pathogenic Avian Influenza (H3N1) Epidemic, Belgium, 2019.

The high economic impact and zoonotic potential of avian influenza call for detailed investigations of dispersal dynamics of epidemics. We integrated phylogeographic and epidemiologic analyses to investigate the dynamics of a low pathogenicity avian influenza (H3N1) epidemic that occurred in Belgium during 2019. Virus genomes from 104 clinical samples originating from 85% of affected farms were sequenced. A spatially explicit phylogeographic analysis confirmed a dominating northeast to southwest dispersal direction and a long-distance dispersal event linked to direct live animal transportation between farms. Spatiotemporal clustering, transport, and social contacts strongly correlated with the phylogeographic pattern of the epidemic. We detected only a limited association between wind direction and direction of viral lineage dispersal. Our results highlight the multifactorial nature of avian influenza epidemics and illustrate the use of genomic analyses of virus dispersal to complement epidemiologic and environmental data, improve knowledge of avian influenza epidemiologic dynamics, and enhance control strategies.

Has Epizootic Become Enzootic? Evidence for a Fundamental Change in the Infection Dynamics of Highly Pathogenic Avian Influenza in Europe, 2021.

Phylogenetic evidence from the recent resurgence of high-pathogenicity avian influenza (HPAI) virus subtype H5N1, clade 2.3.4.4b, observed in European wild birds and poultry since October 2021, suggests at least two different and distinct reservoirs. We propose contrasting hypotheses for this emergence: (i) resident viruses have been maintained, presumably in wild birds, in northern Europe throughout the summer of 2021 to cause some of the outbreaks that are part of the most recent autumn/winter 2021 epizootic, or (ii) further virus variants were reintroduced by migratory birds, and these two sources of reintroduction have driven the HPAI resurgence. Viruses from these two principal sources can be distinguished by their hemagglutinin genes, which segregate into two distinct sublineages (termed B1 and B2) within clade 2.3.4.4b, as well as their different internal gene compositions. The evidence of enzootic HPAI virus circulation during the summer of 2021 indicates a possible paradigm shift in the epidemiology of HPAI in Europe.

A robust, cost-effective and widely applicable whole-genome sequencing protocol for capripoxviruses.

The diseases caused by capripoxviruses (CaPVs) are of major economic concern in sheep, goat and cattle as they are inexorably spreading into non-endemic regions. As CaPV strains are serologically indistinguishable and genetically highly homologous, typing closely related strains can only be achieved by whole genome sequencing. Unfortunately the number of publicly available genomes remains low as most sequencing methods rely on virus isolation. Therefore, we developed a robust, cost-effective and widely applicable method that allows to generate (nearly) complete CaPV genomes directly from clinical samples or commercial vaccine batches. A set of pan-CaPVs long-range PCRs spanning the entire genome was designed to generate PCR amplicons that can be sequenced on commonly used high-throughput sequencing platforms: MiSeq (Illumina), RSII (PacBio) and MinION (Oxford Nanopore Technologies). The robustness of the LR-PCR strategy was evaluated for all 3 members of CaPV directly from a variety of samples, including clinical samples (N = 7), vaccine batches (N = 6), and virus isolates (N = 2). The sequencing method described here allows to reconstruct (nearly) complete CaPV genomes in less than a week and will aid researchers studying closely-related CaPV strains worldwide.

Coding-Complete Sequences of Recombinant Lumpy Skin Disease Viruses Collected in 2020 from Four Outbreaks in Northern Vietnam.

Lumpy skin disease virus (LSDV) causes a severe, systemic, and economically important disease in cattle. Here, we report coding-complete sequences of recombinant LSDVs from four outbreaks in October and November 2020 in northeastern Vietnam.

Metagenomic sequencing determines complete infectious bronchitis virus (avian Gammacoronavirus) vaccine strain genomes and associated viromes in chicken clinical samples.

Infectious bronchitis virus (IBV, genus Gammacoronavirus) causes an economically important and highly contagious disease in chicken. Random primed RNA sequencing was applied to two IBV positive clinical samples and one in ovo-passaged virus. The virome of a cloacal swab pool was dominated by IBV (82% of viral reads) allowing de novo assembly of a GI-13 lineage complete genome with 99.95% nucleotide identity to vaccine strain 793B. In addition, substantial read counts (16% of viral reads) allowed the assembly of a near-complete chicken astrovirus genome, while lower read counts identified the presence of chicken calicivirus and avian leucosis virus. Viral reads in a respiratory/intestinal tissue pool were distributed between IBV (22.53%), Sicinivirus (Picornaviridae, 24%), and avian leucosis virus (37.04%). A complete IBV genome with 99.95% nucleotide identity to vaccine strain H120 (lineage GI-1), as well as a near-complete avian leucosis virus genome and a partial Sicinivirus genome were assembled from the tissue sample data. Lower read counts identified chicken calicivirus, Avibirnavirus (infectious bursal disease virus, assembling to 98.85% of segment A and 69.66% of segment B closely related to D3976/1 from Germany, 2017) and avian orthoreovirus, while three avian orthoavulavirus 1 reads confirmed prior real-time RT-PCR result. IBV sequence variation analysis identified both fixed and minor frequency variations in the tissue sample compared to its in ovo-passaged virus. Metagenomic methods allow the determination of complete coronavirus genomes from clinical chicken samples while providing additional insights in RNA virus sequence diversity and coinfecting viruses potentially contributing to pathogenicity.

Nearly Complete Genome Sequences of Two Bluetongue Viruses Isolated during the 2020 Outbreak in the Grand Duchy of Luxembourg.

Bluetongue is one of the major diseases of ruminants listed by the World Organisation for Animal Health. Bluetongue virus serotype 8 (BTV-8) has been considered enzootic in France since 2018. Here, we report the nearly complete genome sequences of two BTV-8 isolates from the 2020 outbreak in the Grand Duchy of Luxembourg.

Molecular and virological characterization of the first poultry outbreaks of Genotype VII.2 velogenic avian orthoavulavirus type 1 (NDV) in North-West Europe, BeNeLux, 2018.

After two decades free of Newcastle disease, Belgium encountered a velogenic avian orthoavulavirus type 1 epizootic in 2018. In Belgium, 20 cases were diagnosed, of which 15 occurred in hobby flocks, 2 in professional poultry flocks and 3 in poultry retailers. The disease also disseminated from Belgium towards the Grand Duchy of Luxembourg by trade. Independently, the virus was detected once in the Netherlands, almost simultaneously to the first Belgian detection. As such Newcastle disease emerged in the entire BeNeLux region. Both the polybasic sequence of the fusion gene cleavage site and the intracerebral pathotyping assay demonstrated the high pathogenicity of the strain. This paper represents the first notification of this specific VII.2 subgenotype in the North-West of Europe. Time-calibrated full genome phylogenetic analysis indicated the silent or unreported circulation of the virus prior to the emergence of three genetic clusters in the BeNeLux region without clear geographical or other epidemiological correlation. The Dutch strain appeared as an outgroup to the Belgian and Luxembourgian strains in the time-correlated genetic analysis and no epidemiological link could be identified between the Belgian and Dutch outbreaks. In contrast, both genetic and epidemiological outbreak investigation data linked the G.D. Luxembourg case to the Belgian outbreak. The genetic links between Belgian viruses from retailers and hobby flocks only partially correlated with epidemiological data. Two independent introductions into the professional poultry sector were identified, although their origin could not be determined. Animal experiments using 6-week- old specific pathogen-free chickens indicated a systemic infection and efficient transmission of the virus. The implementation of re-vaccination in the professional sector, affected hobby and retailers, as well as the restriction on assembly and increased biosecurity measures, possibly limited the epizootic and resulted in the disappearance of the virus. These findings emphasize the constant need for awareness and monitoring of notifiable viruses in the field.

Complete Coding Sequence of a Lumpy Skin Disease Virus from an Outbreak in Bulgaria in 2016.

Lumpy skin disease (LSD) is an emerging cattle disease with serious economic consequences. We report the complete coding sequence of LSD virus 210LSD-249/BUL/16, detected in a blood sample from a diseased cow during an outbreak in Bulgaria (Kabile Village, Yambol Region) in June 2016.

Increased viral read counts and metagenomic full genome characterization of porcine astrovirus 4 and Posavirus 1 in sows in a swine farm with unexplained neonatal piglet diarrhea.

Neonatal diarrhea in piglets may cause major losses in affected pig herds. The present study used random high-throughput RNA sequencing (metagenomic next generation sequencing, mNGS) to investigate the virome of sows from a farm with persistent neonatal piglet diarrhea in comparison to two control farms without diarrhea problems. A variety of known swine gastrointestinal viruses was detected in the control farms as well as in the problem farm (Mamastrovirus, Enterovirus, Picobirnavirus, Posavirus 1, Kobuvirus, Proprismacovirus). A substantial increase in normalized viral read counts was observed in the affected farm compared to the control farms. The increase was attributable to a single viral species in each of the sampled sows (porcine astrovirus 4 and Posavirus 1). The complete genomes of a porcine astrovirus 4 and two co-infecting Posavirus 1 were de novo assembled and characterized. The 6734 nt single-stranded RNA genome of porcine astrovirus 4 (PoAstV-4) strain Belgium/2019 contains three overlapping open reading frames (nonstructural protein 1ab, nonstructural protein 1a, capsid protein). Posavirus 1 strains Belgium/01/2019 and Belgium/02/2019 have a 9814 nt single-stranded positive-sense RNA genome encoding a single open reading frame (polyprotein precursor) containing the five expected Picornavirales-conserved protein domains. The study highlights the potential of mNGS workflows to study unexplained neonatal diarrhea in piglets and contributes to the scarce availability of both PoAstV-4 and Posavirus-1 whole genome sequences from Western Europe.

Atypical Pathogenicity of Avian Influenza (H3N1) Virus Involved in Outbreak, Belgium, 2019.

In 2019, an outbreak of avian influenza (H3N1) virus infection occurred among commercial poultry in Belgium. Full-genome phylogenetic analysis indicated a wild bird origin rather than recent circulation among poultry. Although classified as a nonnotifiable avian influenza virus, it was associated with reproductive tropism and substantial mortality in the field.

Evaluation of a commercial exogenous internal process control for diagnostic RNA virus metagenomics from different animal clinical samples.

Metagenomic next generation sequencing (mNGS) is increasingly recognized as an important complementary tool to targeted human and animal infectious disease diagnostics. It is, however, sensitive to biases and errors that are currently not systematically evaluated by the implementation of quality controls (QC) for the diagnostic use of mNGS. We evaluated a commercial reagent (Mengovirus extraction control kit, CeraamTools, bioMérieux) as an exogenous internal control for mNGS. It validates the integrity of reagents and workflow, the efficient isolation of viral nucleic acids and the absence of inhibitors in individual samples (verified using a specific qRT-PCR). Moreover, it validates the efficient generation of viral sequence data in individual samples (verified by normalized mengoviral read counts in the metagenomic analysis). We show that when using a completely random metagenomics workflow: (1) Mengovirus RNA can be reproducibly detected in different animal sample types (swine feces and sera, wild bird cloacal swabs), except for tissue samples (swine lung); (2) the Mengovirus control kit does not contain other contaminating viruses that may affect metagenomic experiments (using a cutoff of minimum 1 Kraken classified read per million (RPM)); (3) the addition of 2.17 × 106Mengovirus copies/mL of sample does not affect the virome composition of pig fecal samples or wild bird cloacal swab samples; (4) Mengovirus Cq values (using as cutoff the upper limit of the 99 % confidence interval of Cq values for a given sample matrix) allow the identification of samples with poor viral RNA extraction or high inhibitor load; (5) Mengovirus normalized read counts (cutoff RPM > 1) allow the identification of samples where the viral sequences are outcompeted by host or bacterial target sequences in the random metagenomic workflow. The implementation of two QC testing points, a first one after RNA extraction (Mengoviral qRT-PCR) and a second one after metagenomic data analysis provide valuable information for the validation of individual samples and results. Their implementation in addition to external controls validating runs or experiments should be carefully considered for a given sample type and workflow.

Complete Coding Sequence of a Lumpy Skin Disease Virus Strain Isolated during the 2016 Outbreak in Kazakhstan.

Lumpy skin disease virus (LSDV) causes an economically important disease in cattle. Here, we report the complete coding sequence of the LSDV isolate Kubash/KAZ/16, detected in a clinical sample from an infected cow from the outbreak reported on 7 July 2016 in Kazakhstan (Atyrau Region).

Updated unified phylogenetic classification system and revised nomenclature for Newcastle disease virus.

Several Avian paramyxoviruses 1 (synonymous with Newcastle disease virus or NDV, used hereafter) classification systems have been proposed for strain identification and differentiation. These systems pioneered classification efforts; however, they were based on different approaches and lacked objective criteria for the differentiation of isolates. These differences have created discrepancies among systems, rendering discussions and comparisons across studies difficult. Although a system that used objective classification criteria was proposed by Diel and co-workers in 2012, the ample worldwide circulation and constant evolution of NDV, and utilization of only some of the criteria, led to identical naming and/or incorrect assigning of new sub/genotypes. To address these issues, an international consortium of experts was convened to undertake in-depth analyses of NDV genetic diversity. This consortium generated curated, up-to-date, complete fusion gene class I and class II datasets of all known NDV for public use, performed comprehensive phylogenetic neighbor-Joining, maximum-likelihood, Bayesian and nucleotide distance analyses, and compared these inference methods. An updated NDV classification and nomenclature system that incorporates phylogenetic topology, genetic distances, branch support, and epidemiological independence was developed. This new consensus system maintains two NDV classes and existing genotypes, identifies three new class II genotypes, and reduces the number of sub-genotypes. In order to track the ancestry of viruses, a dichotomous naming system for designating sub-genotypes was introduced. In addition, a pilot dataset and sub-trees rooting guidelines for rapid preliminary genotype identification of new isolates are provided. Guidelines for sequence dataset curation and phylogenetic inference, and a detailed comparison between the updated and previous systems are included. To increase the speed of phylogenetic inference and ensure consistency between laboratories, detailed guidelines for the use of a supercomputer are also provided. The proposed unified classification system will facilitate future studies of NDV evolution and epidemiology, and comparison of results obtained across the world.

Serological and molecular epidemiology of foot-and-mouth disease viruses in agro-pastoralist livestock herds in the kachia grazing reserve, Nigeria.

The Kachia Grazing Reserve (KGR) is located in Kaduna state in north-western Nigeria and consists of 6 contiguous blocks housing 744 defined households (HH), all engaged in livestock keeping. It is considered as a homogenous epidemiological unit and a defined study area. In 2012, all cattle and sheep of 40 selected HH were sampled to determine sero-prevalence of antibodies to foot-and-mouth disease virus (FMDV) and of FMDV. The overall sero-prevalence of antibodies to the non-structural 3ABC protein (NSP-3ABC ELISA) was 28.9% (380/1,315) (30.6% cattle; 16.3% sheep), and in 4.5% (62/1,380) (5% cattle; 0.6% sheep) of the examined sera FMD viral RNA could be detected by real-time RT-PCR (rRT-PCR). Additionally, in 2012 and 2014 serum, epithelium and probang samples were collected from cattle in reported FMD outbreaks and the causative FMDVs were molecularly characterized. Approximately half (28/59) of the outbreak sera reacted positive in NSP-3ABC ELISA, and 88% (52/59) of the outbreak sera contained detectable viral RNA. Overall, antibodies against five FMDV serotypes (O, A, SAT1, SAT2 and SAT3) were detected by solid phase competitive ELISA with combinations of two or more serotypes being common. Of the 21 FMDVs that could be isolated 19 were sequenced and 18 were confirmed as SAT2 (lineage VII) while one was characterized as serotype O (EA-3 topotype). Phylogenetic analysis revealed a close relationship between Nigerian FMDV strains and strains in this region and even with strains in North-Africa. Our findings indicate that FMD constitutes an endemic health problem to cattle rearing in the agro-pastoralist community in the KGR and that the KGR is not a closed epidemiological unit. Insight into the local FMDV epidemiology and in the circulating FMDV serotypes/strains is of support to the relevant authorities in Nigeria when considering the need for an FMD control policy to improve animal production in grazing reserves.

Complete coding sequence of a novel picorna-like virus in a blackbird infected with Usutu virus.

Using random high-throughput RNA sequencing, the complete coding sequence of a novel picorna-like virus (a 9,228-nt contig containing 212,202 reads) was determined from a blackbird (Turdus merula) infected with Usutu virus. This sequence shares only 36% amino acid sequence identity with its closest homolog, arivirus 1, (an unclassified member of the order Picornavirales), and shares its dicistronic genome arrangement. The new virus was therefore tentatively named "blackbird arilivirus" (ari-like virus). The nearly complete genome sequence consists of at least 9,228 nt and contains two open reading frames (ORFs) encoding the nonstructural polyprotein (2235 amino acids) and structural polyprotein (769 amino acids). Two TaqMan RT-qPCR assays specific for ORF1 confirmed the presence of high levels of this novel virus in the original sample. Nucleotide composition analysis suggests that blackbird arilivirus is of dietary (plant) origin.

Complete Genome Sequences of Five Foot-and-Mouth Disease Viruses of Serotype A Isolated from Cattle in Nigeria between 2013 and 2015.

The complete genome sequences of 5 foot-and-mouth disease viruses of serotype A are reported here. These viruses originate from outbreaks in northern Nigeria in 2013 to 2015 and belong to the A/AFRICA/G-IV lineage.

Complete Genome Sequences of Four Foot-and-Mouth Disease Viruses of Serotype South African Territories 1 (SAT 1), Topotype X, Isolated from Cattle in Nigeria in 2015.

The complete genome sequences of four foot-and-mouth disease viruses of South African territories 1 (SAT 1) serotype are reported. These viruses originate from an outbreak in Nigeria in 2015 and belong to the novel SAT 1 topotype X from the west and central African virus pool.

Complete Genome Sequence of the Lumpy Skin Disease Virus Isolated from the First Reported Case in Greece in 2015.

Lumpy skin disease virus (LSDV) causes an economically important disease in cattle. Here, we report the complete genome sequence of the first LSDV isolate identified in mainland Europe. LSDV isolate Evros/GR/15 was isolated from the first cases reported on 18 August 2015 in the Evros region, Greece.