Endocytosis of integrin-binding human picornaviruses

Molecular Epidemiology and Evolution of Coxsackievirus A9

Nineteen CVA9 isolates were obtained between 2010 and 2019 from six provinces of mainland China, using the HFMD surveillance network established in China. Nucleotide sequencing revealed that the full-length VP1 of 19 CVA9 isolates was 906 bases encoding 302 amino acids. The combination of the thresholds of the phylogenetic tree and nucleotide divergence of different genotypes within the same serotype led to a value of 15-25%, and enabled CVA9 worldwide to be categorized into ten genotypes: A-J. The phylogenetic tree showed that the prototype strain was included in genotype A, and that the B, C, D, E, H, and J genotypes disappeared during virus evolution, whereas the F, I, and G genotypes showed co-circulation. Lineage G was the dominant genotype of CVA9 and included most of the strains from nine countries in Asia, North America, Oceania, and Europe. Most Chinese strains belonged to the G genotype, suggesting that the molecular epidemiology of China is consistent with that observed worldwide. The 165 partial VP1 strains (723 nt) showed a mean substitution rate of 3.27 × 10-3 substitution/site/year (95% HPD range 2.93-3.6 × 10-3), dating the tMRCA of CVA9 back to approximately 1922 (1911-1932). The spatiotemporal dynamics of CVA9 showed the spread of CVA9 obviously increased in recent years. Most CVA9 isolates originated in USA, but the epidemic areas of CVA9 are now concentrated in the Asia-Pacific region, European countries, and North America. Recombination analysis within the enterovirus B specie (59 serotypes) revealed eight recombination patterns in China at present, CVB4, CVB5, E30, CVB2, E11, HEV106, HEV85, and HEV75. E14, and E6 may act as recombinant donors in multiple regions. Comparison of temperature sensitivity revealed that temperature-insensitive strains have more amino acid substitutions in the RGD motif of the VP1 region, and the sites T283S, V284M, and R288K in the VP1 region may be related to the temperature tolerance of CVA9.

A phytochemical-based medication search for the SARS-CoV-2 infection by molecular docking models towards spike glycoproteins and main proteases

Identifying best bioactive phytochemicals from different medicinal plants using molecular docking techniques demonstrates a potential pre-clinical compound discovery against SARS-CoV-2 viral infection. The in silico screening of bioactive phytochemicals with the two druggable targets of SARS-CoV-2 by simple precision/extra precision molecular docking methods was used to compute binding affinity at its active sites. phyllaemblicin and cinnamtannin class of phytocompounds showed a better binding affinity range (-9.0 to -8.0 kcal mol-1) towards both these SARS-CoV-2 targets; the corresponding active site residues in the spike protein were predicted as: Y453, Q496, Q498, N501, Y449, Q493, G496, T500, Y505, L455, Q493, and K417; and Mpro: Q189, H164, H163, P168, H41, L167, Q192, M165, C145, Y54, M49, and Q189. Molecular dynamics simulation further established the structural and energetic stability of protein-phytocompound complexes and their interactions with their key residues supporting the molecular docking analysis. Protein-protein docking using ZDOCK and Prodigy server predicted the binding pose and affinity (-13.8 kcal mol-1) of the spike glycoprotein towards the human ACE2 enzyme and also showed significant structural variations in the ACE2 recognition site upon the binding of phyllaemblicin C compound at their binding interface. The phyllaemblicin and cinnamtannin class of phytochemicals can be potential inhibitors of both the spike and Mpro proteins of SARS-CoV-2; furthermore, its pharmacology and clinical optimization would lead towards novel COVID-19 small-molecule therapy.

Isolation and complete genome analysis of a novel duck picornavirus in China

A novel duck picornavirus, designated as duck/AH15/CHN/2015, was isolated and identified from Cherry Valley ducks with short beak and dwarfism syndrome in 2015 in Anhui province of China. Duck/AH15/CHN/2015 has the highest degree of amino acid sequence identity (approximately 43 %) with duck hepatitis A viruses (DHAV) Complete genome analysis revealed that duck/AH15/CHN/2015 possesses a typical picornavirus-like genomic organization, 5' UTR-L-P1 (VP0-VP3-VP1)-P2 (2A1-2A2- 2B-2C)-P3 (3A-3B-3C-3D)-3'UTR-poly (A). The 5'UTR contains a potential type IV internal ribosome entry site, while a conserved "barbell"-like structure is found at the 3'UTR, which is similar to DHAV. Compared to the closest related DHAVs, two unrelated 2A proteins were predicted in duck/AH15/CHN/2015, while three unrelated 2A proteins were presented in DHAVs. Based on the amino acid identity comparison and phylogenetic analysis of P1, 2C, and 3CD (3C and 3D), duck/AH15/CHN/2015 was closely related to but distinct from DHAVs, and it was proposed to be a member of a novel species in the genus Avihepatovirus of the family Picornaviridae.

A new perspective in sepsis treatment: could RGD-dependent integrins be novel targets?

Sepsis is a life-threatening condition caused by the response of the body to an infection, and has recently been regarded as a global health priority because of the lack of effective treatments available. Vascular endothelial cells have a crucial role in sepsis and are believed to be a major target of pathogens during the early stages of infection. Accumulating evidence suggests that common sepsis pathogens, including bacteria, fungi, and viruses, all contain a critical integrin recognition motif, Arg-Gly-Asp (RGD), in their major cell wall-exposed proteins that might act as ligands to crosslink to vascular endothelial cells, triggering systemic dysregulation resulting in sepsis. In this review, we discuss the potential of anti-integrin therapy in the treatment of sepsis and septic shock.

Antifungal Triazole Posaconazole Targets an Early Stage of the Parechovirus A3 Life Cycle

Viruses in species Parechovirus A (Picornaviridae) are associated with a wide variety of clinical manifestations. Parechovirus A3 (PeV-A3) is known to cause sepsis-like illness, meningitis, and encephalitis in infants and young children. To date, no specific therapies are available to treat PeV-A3-infected children. We had previously identified two FDA-cleared antifungal drugs, itraconazole (ITC) and posaconazole (POS), with potent and specific antiviral activity against PeV-A3. Time-of-addition and synchronized infection assays revealed that POS targets an early stage of the PeV-A3 life cycle. POS exerts an antiviral effect, evidenced by a reduction in viral titer following the addition of POS to Vero-P cells before infection, coaddition of POS and PeV-A3 to Vero-P cells, incubation of POS and PeV-A3 prior to Vero-P infection, and at attachment. POS exerts less of an effect on virus entry. A PeV-A3 enzyme-linked immunosorbent assay inhibition experiment, using an anti-PeV-A3 monoclonal antibody, suggested that POS binds directly to the PeV-A3 capsid. POS-resistant PeV-A3 strains developed by serial passage in the presence of POS acquired substitutions in multiple regions of the genome, including the capsid. Reverse genetics confirmed substitutions in capsid proteins VP0, VP3, and VP1 and nonstructural proteins 2A and 3A. Single mutants VP0_K66R, VP0_A124T, VP3_N88S, VP1_Y224C, 2A_S78L, and 3A_T1I were 4-, 9-, 12-, 34-, 51-, and 119-fold more resistant to POS, respectively, than the susceptible prototype strain. Our studies demonstrate that POS may be a valuable tool in developing an antiviral therapy for PeV-A3.

Recombination Events and Conserved Nature of Receptor Binding Motifs in Coxsackievirus A9 Isolates

Coxsackievirus A9 (CVA9) is an enterically transmitted enterovirus and one of the most pathogenic type among human enteroviruses. CVA9 isolates use a distinctive RGD (Arg-Gly-Asp) motif within VP1 capsid protein that defines its ability to bind to integrin receptor(s) for cellular entry. To investigate CVA9 evolution and pathogenicity, genetic relationships and recombination events were analyzed between 54 novel clinical isolates of CVA9, as well as 21 previously published full length CVA9 sequences from GenBank. Samples were investigated by partial sequencing of the novel VP1 and 3Dpol genes, as well as including the corresponding areas from GenBank sequences. Phylogenetic analyses were combined with clinical data in a further attempt to analyze whether sequence evolution reflects CVA9 pathogenicity in the phylogenies. Furthermore, VP1 gene was also analyzed for receptor binding sites including the RGD motif and the putative heparan sulfate (HS) site. Analysis of the 559-nucleotide-long VP1 sequences identified six clades. Although most of the strains within each clade showed geographical clustering, the grouping pattern of the isolates in the analysis of the VP1 gene was strikingly different from grouping of 3Dpol, which suggests that recombination events may have occurred in the region encoding the nonstructural proteins. Inclusion of clinical data did not provide any evidence of symptom based phylogenetic clustering of CVA9 isolates. Amino acid sequence analysis of the VP1 polypeptide demonstrated that the RGD motif was fully conserved among the isolates while the putative HS binding site was only found in one isolate. These data suggest that integrin binding is essential for virus tropism, but do not explain the symptom repertoire.

Cellular entry and uncoating of naked and quasi-enveloped human hepatoviruses

Many ‘non-enveloped’ viruses, including hepatitis A virus (HAV), are released non-lytically from infected cells as infectious, quasi-enveloped virions cloaked in host membranes. Quasi-enveloped HAV (eHAV) mediates stealthy cell-to-cell spread within the liver, whereas stable naked virions shed in feces are optimized for environmental transmission. eHAV lacks virus-encoded surface proteins, and how it enters cells is unknown. We show both virion types enter by clathrin- and dynamin-dependent endocytosis, facilitated by integrin β₁, and traffic through early and late endosomes. Uncoating of naked virions occurs in late endosomes, whereas eHAV undergoes ALIX-dependent trafficking to lysosomes where the quasi-envelope is enzymatically degraded and uncoating ensues coincident with breaching of endolysosomal membranes. Neither virion requires PLA2G16, a phospholipase essential for entry of other picornaviruses. Thus naked and quasi-enveloped virions enter via similar endocytic pathways, but uncoat in different compartments and release their genomes to the cytosol in a manner mechanistically distinct from other Picornaviridae.

The Hepatitis A virus is a common cause of liver disease in humans. It is unable to multiply on its own so it needs to enter the cells of its host and hijack them to make new virus particles.

Infected human cells produce two different types of Hepatitis A particles. The first, known as ‘naked’ virus particles, consist of molecules of ribonucleic acid (or RNA for short) that are surrounded by a protein shell. Naked virus particles are shed in the feces of infected individuals and are very stable, allowing the virus to spread in the environment to find new hosts.

At the same time, a second type of particle, known as the ‘quasi-enveloped’ virus, circulates in the blood of the infected individual. In a quasi-enveloped particle, the RNA and protein shell are completely enclosed within a membrane that is released from the host cell. This membrane protects the protein shell from human immune responses, enabling quasi-enveloped virus particles to spread in a stealthy fashion within the liver.

It was not clear how these two different types of virus particle are both able to enter cells despite their surface being so different. To address this question, Rivera-Serrano et al. used a microscopy approach to observe Hepatitis A particles infecting human liver cells.

The experiments showed that both types of virus particle actually use similar routes. First, the external membrane of the cell folded around the particles, creating a vesicle that trapped the viruses and brought them within the cell. Inside these vesicles, the naked virus particles soon fell apart, and their RNA was released directly into the interior of the cell.

However, the vesicles that carried quasi-enveloped virus travelled further into the cell and eventually delivered their contents to a specialized compartment, the lysosome, where the virus membrane was degraded. This caused the quasi-enveloped viruses to fall apart and release their RNA into the cell more slowly than the naked particles.

Several viruses, such as the one that causes polio, also have quasi-enveloped forms. Studying how these particles are able to infect human cells while hiding behind membranes borrowed from the host may help us target these viruses better.

Progress in human picornavirus research: New findings from the AIROPico consortium

Several research groups in Europe are active on different aspects of human picornavirus research. The AIROPico (Academia-Industry R&D Opportunities for Picornaviruses) consortium combined the disciplines of pathogenesis, diagnostics and therapy development in order to fill the gaps in our understanding of how picornaviruses cause human disease and how to combat them. AIROPico was the first EU consortium dedicated to human picornavirus research and development, and has largely accelerated and improved R&D on picornavirus biology, diagnostics and therapy. In this article, we present the progress on pathogenesis, diagnostics and treatment strategy developments for human picornaviruses resulting from the structured, translational research approach of the AIROPico consortium. We here summarize new insights in protection against infection by maternal or cross-protective antibodies, the visualisation of interactions between virus and neutralizing antibodies by cryoEM structural imaging, and the outcomes from a picornavirus-infected human 3D organoid. Progress in molecular detection and a fast typing assay for rhinovirus species are presented, as well as the identification of new compounds potentially interesting as therapeutic compounds.

Integrin αvβ3 promotes infection by Japanese encephalitis virus

Japanese encephalitis virus (JEV) is a mosquito-borne flavivirus that is one of the major causes of viral encephalitis diseases worldwide. The JEV envelope protein facilitates viral entry, and its domain III contains an Arg-Gly-Asp (RGD) motif, that may modulate JEV entry through the RGD-binding integrin. In this study, the roles of integrin αv and β3 on the infection of JEV were evaluated. Reduced expression of integrin αv/β3 by special shRNA confers 2 to 4-fold inhibition of JEV replication in BHK-21 cells. Meanwhile, antibodies specific for integrin αv/β3 displayed ~58% and ~33% inhibition of JEV infectivity and RGD-specific peptides produced ~36% of inhibition. Expression of E protein and JEV RNA loads were clearly increased in CHO cells transfected with cDNA encoding human integrin β3. Moreover, integrin αv mediates JEV infection in viral binding stage of life cycle. Therefore, our study suggested that integrin αv and β3 serve as a host factor associated with JEV entry into the target cells.

Integrins are not essential for entry of coxsackievirus A9 into SW480 human colon adenocarcinoma cells

Background

Coxsackievirus A9 (CV-A9) is a pathogenic enterovirus type within the family Picornaviridae. CV-A9 infects A549 human epithelial lung carcinoma cells by attaching to the αVβ6 integrin receptor through a highly conserved Arg-Gly-Asp (RGD) motif, which is located at the exposed carboxy-terminus of the capsid protein VP1 detected in all studied clinical isolates. However, genetically-modified CV-A9 that lacks the RGD motif (CV-A9-RGDdel) has been shown to be infectious in some cell lines but not in A549, suggesting that RGD-mediated integrin binding is not always essential for efficient entry of CV-A9.

Methods

Two cell lines, A549 and SW480, were used in the study. SW480 was the study object for the integrin-independent entry and A549 was used as the control for integrin-dependent entry. Receptor levels were quantitated by cell sorting and quantitative PCR. Antibody blocking assay and siRNA silencing of receptor-encoding genes were used to block virus infection. Peptide phage display library was used to identify peptide binders to CV-A9. Immunofluorescence and confocal microscopy were used to visualize the virus infection in the cells.

Results

We investigated the receptor use and early stages of CV-A9 internalization to SW480 human epithelial colon adenocarcinoma cells. Contrary to A549 infection, we showed that both CV-A9 and CV-A9-RGDdel internalized into SW480 cells and that function-blocking anti-αV integrin antibodies had no effect on the binding and entry of CV-A9. Whereas siRNA silencing of β6 integrin subunit had no influence on virus infection in SW480, silencing of β2-microglobulin (β2M) inhibited the virus infection in both cell lines. By using a peptide phage display screening, the virus-binding peptide identical to the N-terminal sequence of HSPA5 protein was identified and shown to block the virus infection in both A549 and SW480 cell lines. HSPA5 was also found to co-localize with CV-A9 at the SW480 cell periphery during the early stages of infection by confocal microscopy.

Conclusions

The data suggest that while αVβ6 integrin is essential for CV-A9 in A549 cell line, it is not required in SW480 cell line in which β2M and HSPA5 alone are sufficient for CV-A9 infection. This suggests that the choice of CV-A9 receptor(s) is dependent on the tissue/cellular environment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0619-y) contains supplementary material, which is available to authorized users.

Human Parechovirus 1 Infection Occurs via αVβ1 Integrin

Human parechovirus 1 (HPeV-1) (family Picornaviridae) is a global cause of pediatric respiratory and CNS infections for which there is no treatment. Although biochemical and in vitro studies have suggested that HPeV-1 binds to αVβ1, αVβ3 and αVβ6 integrin receptor(s), the actual cellular receptors required for infectious entry of HPeV-1 remain unknown. In this paper we analyzed the expression profiles of αVβ1, αVβ3, αVβ6 and α5β1 in susceptible cell lines (A549, HeLa and SW480) to identify which integrin receptors support HPeV-1 internalization and/or replication cycle. We demonstrate by antibody blocking assay, immunofluorescence microscopy and RT-qPCR that HPeV-1 internalizes and replicates in cell lines that express αVβ1 integrin but not αVβ3 or αVβ6 integrins. To further study the role of β1 integrin, we used a mouse cell line, GE11-KO, which is deficient in β1 expression, and its derivate GE11-β1 in which human integrin β1 subunit is overexpressed. HPeV-1 (Harris strain) and three clinical HPeV-1 isolates did not internalize into GE11-KO whereas GE11-β1 supported the internalization process. An integrin β1-activating antibody, TS2/16, enhanced HPeV-1 infectivity, but infection occurred in the absence of visible receptor clustering. HPeV-1 also co-localized with β1 integrin on the cell surface, and HPeV-1 and β1 integrin co-endocytosed into the cells. In conclusion, our results demonstrate that in some cell lines the cellular entry of HPeV-1 is primarily mediated by the active form of αVβ1 integrin without visible receptor clustering.

Therapeutic Use of Native and Recombinant Enteroviruses

Research on human enteroviruses has resulted in the identification of more than 100 enterovirus types, which use more than 10 protein receptors and/or attachment factors required in cell binding and initiation of the replication cycle. Many of these “viral” receptors are overexpressed in cancer cells. Receptor binding and the ability to replicate in specific target cells define the tropism and pathogenesis of enterovirus types, because cellular infection often results in cytolytic response, i.e., disruption of the cells. Viral tropism and cytolytic properties thus make native enteroviruses prime candidates for oncolytic virotherapy. Copy DNA cloning and modification of enterovirus genomes have resulted in the generation of enterovirus vectors with properties that are useful in therapy or in vaccine trials where foreign antigenic epitopes are expressed from or on the surface of the vector virus. The small genome size and compact particle structure, however, set limits to enterovirus genome modifications. This review focuses on the therapeutic use of native and recombinant enteroviruses and the methods that have been applied to modify enterovirus genomes for therapy.

Role of Heparan Sulfate in Cellular Infection of Integrin-Binding Coxsackievirus A9 and Human Parechovirus 1 Isolates

Heparan sulfate/heparin class of proteoglycans (HSPG) have been shown to function in cellular attachment and infection of numerous viruses including picornaviruses. Coxsackievirus A9 (CV-A9) and human parechovirus 1 (HPeV-1) are integrin-binding members in the family Picornaviridae. CV-A9 Griggs and HPeV-1 Harris (prototype) strains have been reported not to bind to heparin, but it was recently shown that some CV-A9 isolates interact with heparin in vitro via VP1 protein with a specific T132R/K mutation. We found that the infectivity of both CV-A9 Griggs and HPeV-1 Harris was reduced by sodium chlorate and heparinase suggestive of HSPG interactions. We analyzed the T132 site in fifty-four (54) CV-A9 clinical isolates and found that only one of them possessed T132/R mutation while the other nine (9) had T132K. We then treated CV-A9 Griggs and HPeV-1 Harris and eight CV-A9 and six HPeV-1 clinical isolates with heparin and protamine. Although infectivity of Griggs strain was slightly reduced (by 25%), heparin treatment did not affect the infectivity of the CV-A9 isolates that do not possess the T132R/K mutation, which is in line with the previous findings. Some of the HPeV-1 isolates were also affected by heparin treatment, which suggested that there may be a specific heparin binding site in HPeV-1. In contrast, protamine (a specific inhibitor of heparin) completely inhibited the infection of both prototypes and clinical CV-A9 and HPeV-1 isolates. We conclude that T132R/K mutation has a role in heparin binding of CV-A9, but we also show data, which suggest that there are other HSPG binding sites in CV-A9. In all, we suggest that HSPGs play a general role in both CV-A9 and HPeV-1 infections.

Plaque2.0-A High-Throughput Analysis Framework to Score Virus-Cell Transmission and Clonal Cell Expansion

Classical plaque assay measures the propagation of infectious agents across a monolayer of cells. It is dependent on cell lysis, and limited by user-specific settings and low throughput. Here, we developed Plaque2.0, a broadly applicable, fluorescence microscopy-based high-throughput method to mine patho-biological clonal cell features. Plaque2.0 is an open source framework to extract information from chemically fixed cells by immuno-histochemistry or RNA in situ hybridization, or from live cells expressing GFP transgene. Multi-parametric measurements include infection density, intensity, area, shape or location information at single plaque or population levels. Plaque2.0 distinguishes lytic and non-lytic spread of a variety of DNA and RNA viruses, including vaccinia virus, adenovirus and rhinovirus, and can be used to visualize simultaneous plaque formation from co-infecting viruses. Plaque2.0 also analyzes clonal growth of cancer cells, which is relevant for cell migration and metastatic invasion studies. Plaque2.0 is suitable to quantitatively analyze virus infections, vector properties, or cancer cell phenotypes.

Differential permissivity of human cerebrovascular endothelial cells to enterovirus infection and specificities of serotype EV-A71 in crossing an in vitro model of the human blood-brain barrier

Human cerebral microvascular endothelial cells (hCMEC/D3 cell line) form a steady polarized barrier when cultured in vitro on a permeable membrane. Their susceptibility to enterovirus (EV) strains was analysed to investigate how these viruses may cross the blood-brain barrier. A sample of 88 virus strains was selected on phylogenetic features amongst 43 epidemiologically relevant types of the four EV species A-D. The EV-A71 genome was replicated at substantial rates, whilst the infectious virus was released at extremely low but sustained rates at both barrier sides for at least 4 days. EV-A71 antigens were detected in a limited number of cells. The properties of the endothelial barrier (structure and permeability) remained intact throughout infection. The chronic EV-A71 infection was in sharp contrast to the productive infection of cytolytic EVs (e.g. echoviruses E-6 and E-30). The hCMEC/D3 barriers infected with the latter EVs exhibited elevated proportions of apoptotic and necrotic cells, which resulted in major injuries to the endothelial barriers with a dramatic increase of paracellular permeability and virus crossing to the abluminal side. The following intracellular rearrangements were also seen: early destruction of the actin cytoskeleton, remodelling of intracellular membranes and reorganization of the mitochondrion network in a small cluster near the perinuclear space.

Cellular density effect on RGD ligand internalization in glioblastoma for MRI application

Cellular density is a parameter measured for glioma grade and invasiveness diagnosis. The characterization of the cellular density can be performed, non invasively, by magnetic resonance imaging (MRI), since, this technique displays a good resolution. Nevertheless MRI sensitivity is critical. Development of smart contrast agents appears useful to increase MRI signal to noise ratio (SNR). Tumor invasiveness is correlated with high expression of integrins that can be targeted by RGD motif. In this study, MRI contrast agents or fluorescent probes linked to RGD-peptides were used, in a glioma model, to assess the relation between RGD uptake/signal improvement/cell density and consequently tumor invasiveness.

Experiments were performed in vitro with U87-MG glioma cells. Flow cytometry and microscopy experiments with RGD and iRGD-alexa488 demonstrated that cell internalization was dependent on cell density. The internalization involved a clathrin-dependent endocytosis. Cytoskeleton and particularly the microtubules were concerned. Actin filaments played a minor role. The internalization was also dependent on the glycolysis and the oxidative phosphorylations. The cellular density modulated the importance of the endocytosis pathways and of the metabolism but not the cytoskeleton contribution. The internalization of the RGD-peptide associated to gadolinium chelate increased the SNR of U87 cells. Moreover, following the cell density augmentation, the SNR increased with a low amplitude but a trend was clearly determined. In conclusion, RGD-peptide internalization appeared, in vitro, as a marker of cellular density. In perspective, the combination of these peptides with contrast agents associated to more sensitive MRI techniques could improve the MRI signal allowing the characterization of cellular density for tumor diagnosis.

Genomic characterization of Sebokele virus 1 (SEBV1) reveals a new candidate species among the genus Parechovirus

We determined the genomic features and the taxonomic classification of Sebokele virus 1 (SEBV1), a previously unclassified arbovirus isolated in 1972 from rodents collected in Botambi, Central African Republic. The complete genome sequence was obtained using a deep sequencing approach (Illumina technology) and dedicated bioinformatics workflows for data analysis. Molecular analysis identified SEBV1 as a picornavirus, most closely related to Ljungan viruses of the genus Parechovirus. The genome has a typical Ljungan virus-like organization, including the presence of two unrelated 2A protein motifs. Phylogenetic analysis confirmed that SEBV1 belongs to the parechovirus phylogroup and was most closely related to the Ljungan virus species. However, it appeared clearly distinct from all members of this phylogroup, suggesting that it represents a novel species of the genus Parechovirus.