Frequency and dynamics of recombination within different species of human enteroviruses

A review of the recombination events, mechanisms and consequences of Coxsackievirus A6

Hand, foot, and mouth disease (HFMD) is one of the most common class C infectious diseases, posing a serious threat to public health worldwide. Enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) have been regarded as the major pathogenic agents of HFMD; however, since an outbreak caused by coxsackievirus A6 (CV-A6) in France in 2008, CV-A6 has gradually become the predominant pathogen in many regions. CV-A6 infects not only children but also adults, and causes atypical clinical symptoms such as a more generalized rash, eczema herpeticum, high fever, and onychomadesis, which are different from the symptoms associated with EV-A71 and CV-A16. Importantly, the rate of genetic recombination of CV-A6 is high, which can lead to changes in virulence and the rapid evolution of other characteristics, thus posing a serious threat to public health. To date, no specific vaccines or therapeutics have been approved for CV-A6 prevention or treatment, hence it is essential to fully understand the relationship between recombination and evolution of this virus. Here, we systematically review the genetic recombination events of CV-A6 that have occurred worldwide and explore how these events have promoted virus evolution, thus providing important information regarding future HFMD surveillance and prevention.

Genetic diversity and recombination of bovine enterovirus strains in China

Bovine enterovirus (BEV) consisting of enterovirus species E (EV-E) and F (EV-F) is the causative agent associated with respiratory and gastrointestinal diseases in cattle. Here, we reported the characterization, genetic diversity, and recombination of novel BEV strains isolated from the major cattle-raising regions in China during 2012-2018. Twenty-seven BEV strains were successfully isolated and characterized. Molecular characterization demonstrated that the majority of these novel BEV strains (24/27) were EV-E, while only few strains (3/27) were EV-F. Sequence analysis revealed the diversity of the circulating BEV strains such as species and subtypes where different species or subtype coinfections were detected in the same regions and even in the same cattle herds. For the EV-E, two novel subtypes, designated as EV-E6 and EV-E7, were revealed in addition to the currently reported EV-E1-EV-E5. Comparative genomic analysis revealed the intraspecies and interspecies genetic exchanges among BEV isolates. The representative strain HeN-B62 was probably from AN12 (EV-F7) and PS-87-Belfast (EV-F3) strains. The interspecies recombination between EV-E and EV-F was also discovered, where the EV-F7-AN12 might be from EV-E5 and EV-F1, and EV-E5-MexKSU/5 may be recombined from EV-F7 and EV-E1. The aforementioned results revealed the genetic diversity and recombination of novel BEV strains and unveiled the different BEV species or subtype infections in the same cattle herd, which will broaden the understanding of enterovirus genetic diversity, recombination, pathogenesis, and prevention of disease outbreaks.

IMPORTANCE

Bovine enterovirus (BEV) infection is an emerging disease in China that is characterized by digestive, respiratory, and reproductive disorders. In this study, we first reported two novel EV-E subtypes detected in cattle herds in China, unveiled the coinfection of two enterovirus species (EV-E/EV-F) and different subtypes (EV-E2/EV-E7, EV-E1/EV-E7, and EV-E3/EV-E6) in the same cattle herds, and revealed the enterovirus genetic exchange in intraspecies and interspecies recombination. These results provide an important update of enterovirus prevalence and epidemiological aspects and contribute to a better understanding of enterovirus genetic diversity, evolution, and pathogenesis.

The evolutionary phylodynamics of human parechovirus A type 3 reveal multiple recombination events in South Korea

Human parechovirus A (HPeV-A) is a causative agent of respiratory and gastrointestinal illnesses, acute flaccid paralysis encephalitis, meningitis, and neonatal sepsis. To clarify the characteristics of HPeV-A infection in children, 391 fecal specimens were collected from January 2014 to October 2015 from patients with acute gastroenteritis in Seoul, South Korea. Of these, 221/391 (56.5%) HPeV-A positive samples were found in children less than 2 years old. Three HPeV-A genotypes HPeV-A1 (117/221; 52.94%), HPeV-A3 (100/221; 45.25%), and HPeV-A6 (4/221; 1.81%) were detected, among which HPeV-A3 was predominant with the highest recorded value of 58.6% in 2015. Moreover, recombination events in the Korean HPeV-A3 strains were detected. Phylogenetic analysis revealed that the capsid-encoding regions and noncapsid gene 2A of the four Korean HPeV-A3 strains are closely related to the HPeV-A3 strains isolated in Canada in 2007 (Can82853-01), Japan in 2008 (A308/99), and Taiwan in 2011 (TW-03067-2011) while noncapsid genes P2 (2B-2C) and P3 (3A-3D) are closely related to those of HPeV-A1 strains BNI-788St (Germany in 2008) and TW-71594-2010 (Taiwan in 2010). This first report on the whole-genome analysis of HPeV-A3 in Korea provides insight into the evolving status and pathogenesis of HPeVs in children.

Investigating the virulence of coxsackievirus B6 strains and antiviral treatments in a neonatal murine model

Coxsackievirus B6 (CVB6), a member of the human enterovirus family, is associated with severe diseases such as myocarditis in children. However, to date, only a limited number of CVB6 strains have been identified, and their characterization in animal models has been lacking. To address this gap, in this study, a neonatal murine model of CVB6 infection was established to compare the replication and virulence of three infectious-clone-derived CVB6 strains in vivo. The results showed that following challenge with a lethal dose of CVB6 strains, the neonatal mice rapidly exhibited a series of clinical signs, such as weight loss, limb paralysis, and death. For the two high-virulence CVB6 strains, histological examination revealed myocyte necrosis in skeletal and cardiac muscle, and immunohistochemistry confirmed the expression of CVB6 viral protein in these tissues. Real-time PCR assay also revealed higher viral loads in the skeletal and cardiac muscle than in other tissues at different time points post infection. Furthermore, the protective effect of passive immunization with antisera and a neutralizing monoclonal antibody against CVB6 infection was evaluated in the neonatal mouse model. This study should provide insights into the pathogenesis of CVB6 and facilitate further research in the development of vaccines and antivirals against CVBs.

Molecular characterization and pathogenicity evaluation of enterovirus G isolated from diarrheic piglets

IMPORTANCE

Enterovirus G is a species of positive-sense single-stranded RNA viruses associated with several mammalian diseases. The porcine enterovirus strains isolated here were chimeric viruses with the PLCP gene of porcine torovirus, which grouped together with global EV-G1 strains. The isolated EV-G strain could infect various cell types from different species, suggesting its potential cross-species infection risk. Animal experiment showed the pathogenic ability of the isolated EV-G to piglets. Additionally, the EV-Gs were widely distributed in the swine herds. Our findings suggest that EV-G may have evolved a novel mechanism for broad tropism, which has important implications for disease control and prevention.

Genetic diversity and characterization of rhinoviruses from Chinese clinical samples with a global perspective

IMPORTANCE

Based on clinical samples collected in China, we detected and reported 22 types for the first time in China, as well as three types for the first time in Asia, and reported their genetic characteristics and diversity. We identified a novel type of Rhinovirus (RV), A110, highlighting its unique genetic features. We annotated the genomic structure and serotype of all the existing RV sequences in the database, and four novel RV types were identified and their genetic diversity reported. Combined with the sequence annotation, we constructed a complete VP1 data set of RV and conducted the first large-scale evolutionary dynamics analysis of RV. Based on a high-quality data set, we conducted a comprehensive analysis of the guanine-cytosine (GC) content variations among serotypes of RVs. This study provides crucial theoretical support and valuable data for understanding RV's genetic diversity and developing antiviral strategies.

Molecular Epidemiology and Evolution of Coxsackievirus A14

As the proportion of non-enterovirus 71 and non-coxsackievirus A16 which proportion of composition in the hand, foot, and mouth pathogenic spectrum gradually increases worldwide, the attention paid to other enteroviruses has increased. As a member of the species enterovirus A, coxsackievirus A14 (CVA14) has been epidemic around the world until now since it has been isolated. However, studies on CVA14 are poor and the effective population size, evolutionary dynamics, and recombination patterns of CVA14 are not well understood. In this study, 15 CVA14 strains were isolated from HFMD patients in mainland China from 2009 to 2019, and the complete sequences of CVA14 in GenBank as research objects were analyzed. CVA14 was divided into seven genotypes A-G based on an average nucleotide difference of the full-length VP1 coding region of more than 15%. Compared with the CVA14 prototype strain, the 15 CVA14 strains showed 84.0-84.7% nucleotide identity in the complete genome and 96.9-97.6% amino acid identity in the encoding region. Phylodynamic analysis based on 15 CVA14 strains and 22 full-length VP1 sequences in GenBank showed a mean substitution rate of 5.35 × 10-3 substitutions/site/year (95% HPD: 4.03-6.89 × 10-3) and the most recent common ancestor (tMRCA) of CVA14 dates back to 1942 (95% HPD: 1930-1950). The Bayesian skyline showed that the effective population size had experienced a decrease-increase-decrease fluctuation since 2004. The phylogeographic analysis indicated two and three possible migration paths in the world and mainland China, respectively. Four recombination patterns with others of species enterovirus A were observed in 15 CVA14 strains, among which coxsackievirus A2 (CVA2), coxsackievirus A4 (CVA4), coxsackievirus A6 (CVA6), coxsackievirus A8 (CVA8), and coxsackievirus A12 (CVA12) may act as recombinant donors in multiple regions. This study has filled the gap in the molecular epidemiological characteristics of CVA14, enriched the global CVA14 sequence database, and laid the epidemiological foundation for the future study of CVA14 worldwide.

First Identification and Pathogenicity Evaluation of an EV-G17 Strain Carrying a Torovirus Papain-like Cysteine Protease (PLCP) Gene in China

Enterovirus G (EV-G) is prevalent in pig populations worldwide, and a total of 20 genotypes (G1 to G20) have been confirmed. Recently, recombinant EV-Gs carrying the papain-like cysteine protease (PLCP) gene of porcine torovirus have been isolated or detected, while their pathogenicity is poorly understood. In this study, an EV-G17-PLCP strain, 'EV-G/YN23/2022', was isolated from the feces of pigs with diarrhea, and the virus replicated robustly in numerous cell lines. The isolate showed the highest complete genome nucleotide (87.5%) and polyprotein amino acid (96.6%) identity in relation to the G17 strain 'IShi-Ya4' (LC549655), and a possible recombination event was detected at the 708 and 3383 positions in the EV-G/YN23/2022 genome. EV-G/YN23/2022 was nonlethal to piglets, but mild diarrhea, transient fever, typical skin lesions, and weight gain deceleration were observed. The virus replicated efficiently in multiple organs, and the pathological lesions were mainly located in the small intestine. All the challenged piglets showed seroconversion for EV-G/YN23/2022 at 6 to 9 days post-inoculation (dpi), and the neutralization antibody peaked at 15 dpi. The mRNA expression levels of IL-6, IL-18, IFN-α, IFN-β, and ISG-15 in the peripheral blood mononuclear cells (PBMCs) were significantly up-regulated during viral infection. This is the first documentation of the isolation and pathogenicity evaluation of the EV-G17-PLCP strain in China. The results may advance our understanding of the evolution characteristics and pathogenesis of EV-G-PLCP.

Novel recombinant porcine enterovirus G viruses lacking structural proteins are maintained in pig farms in Japan

Type 1 recombinant enterovirus G (EV-G), which carries the papain-like cysteine protease (PLCP) gene of torovirus between its 2C/3A regions, and type 2 recombinant EV-G, which carries the torovirus PLCP gene with its flanking regions having non-EV-G sequences in place of the viral structural genes, have been detected in pig farms in several countries. In a previous study, we collected 222 fecal samples from 77 pig farms from 2104 to 2016 and detected one type 2 recombinant EV-G genome by metagenomics sequencing. In this study, we reanalyzed the metagenomic data and detected 11 type 2 recombinant EV-G genomes. In addition, we discovered new type 2 recombinant EV-G genomes of the two strains from two pig farms samples in 2018 and 2019. Thus, we identified the genomes of 13 novel type 2 recombinant EV-Gs isolated from several pig farms in Japan. Type 2 recombinant EV-G has previously been detected only in neonatal piglets. The present findings suggest that type 2 recombinant EV-G replicates in weaning piglets and sows. The detection of type 1 recombinant EV-Gs and type 2 recombinant EV-Gs at 3-year and 2-year intervals, respectively, from the same pig farm suggests that the viruses were persistently infecting or circulating in these farms.

Prevalence and genetic diversity in bovine parechovirus infecting Japanese cattle

The first bovine parechovirus (Bo_ParV) was reported in 2021, and currently, only two nearly complete genome sequences of Bo_ParV are available. In this study, we detected Bo_ParVs in 10 out of 158 bovine fecal samples tested using real-time RT-PCR, and Bo_ParVs were isolated from three of these samples using MA104 cells. Analysis of the P1 region revealed that Bo_ParVs shared high pairwise amino acid sequence similarity (≥ 95.7% identity), suggesting antigenic similarity among Bo_ParVs, whereas nucleotide sequence identity values (≥ 84.8%) indicated more variability. A recombination breakpoint was identified in the 2B region, which may influence the evolution of this virus.

Patterns and Temporal Dynamics of Natural Recombination in Noroviruses

Noroviruses infect a wide range of mammals and are the major cause of gastroenteritis in humans. Recombination at the junction of ORF1 encoding nonstructural proteins and ORF2 encoding major capsid protein VP1 is a well-known feature of noroviruses. Using all available complete norovirus sequences, we systematically analyzed patterns of natural recombination in the genus Norovirus both throughout the genome and across the genogroups. Recombination events between nonstructural (ORF1) and structural genomic regions (ORF2 and ORF3) were found in all analyzed genogroups of noroviruses, although recombination was most prominent between members of GII, the most common genogroup that infects humans. The half-life times of recombinant forms (clades without evidence of recombination) of human GI and GII noroviruses were 10.4 and 8.4-11.3 years, respectively. There was evidence of many recent recombination events, and most noroviruses that differed by more than 18% of nucleotide sequence were recombinant relative to each other. However, there were no distinct recombination events between viruses that differed by over 42% in ORF2/3, consistent with the absence of systematic recombination between different genogroups. The few inter-genogroup recombination events most likely occurred between ancient viruses before they diverged into contemporary genogroups. The recombination events within ORF1 or between ORF2/3 were generally rare. Thus, noroviruses routinely exchange full structural and nonstructural blocks of the genome, providing a modular evolution.

Enterovirus B types cause severe infection in infants aged 0-3 months

BACKGROUND

Enterovirus (EV) infections are being increasingly seen in younger infants, often being more severe than in older children. The risk factors of EV infection in infants have been inadequately investigated till date.

METHODS

We conducted a retrospective study on hospitalized children with laboratory-confirmed EV infection (50 infants aged 0-3 months and 65 older than 3 months) at a tertiary care center in China. Prevalence, clinical characteristics, and genetic features of the virus were analyzed, and independent predictors for severe infection were assessed.

RESULTS

Clinical findings showed that severe infection was more common in infants aged 0-3 months than in older children (78.0% vs. 35.4%, p < 0.001), with higher morbidity of pneumonia, meningitis, and sepsis (p < 0.01). EV-B types were detected more frequently in infants aged 0-3 months than in older children (88.0% vs. 7.7%, p < 0.001). Echovirus 11 was the most identified EV-B, and it recombined with E6 in P2 and P3 regions. Risk factors for severe EV infection included EV-B types infection, age less than 3 months, elevated alanine aminotransferase level, abnormal platelet count, and abnormal cerebrospinal fluid characteristics.

CONCLUSIONS

Our data indicated that EV-B types mainly cause severe infection in infants aged 0-3 months. Therefore, knowledge about EV-B types could have implications in designing effective intervention and prevention strategies for young infants with severe EV infection.

First detection and molecular characteristics of bopivirus from goats in China

A metavirome analysis was performed and detected bopivirus in the diarrhoeal fecal samples of goats in China. A total of 136 fecal samples were collected from yeanlings between the dates of June 2021 and January 2022 in Sichuan province, China. Moreover, "Bopivirus B" strains were detected by a specific RT-PCR targeting the 3D gene of the virus. The results showed that the overall detection rate of "Bopivirus B" was 19.12% (26/136). Additionally, there was a higher detection rate (24.05%, 19/79) in the fecal samples collected from yeanlings with diarrhea compared to those from asymptomatic animals (12.28%, 7/57). In these samples, no other common diarrhea-causing pathogens were detected except for three enteric viruses, namely caprine enterovirus, caprine kobuvirus and caprine hunnivirus (with detection rates of 13.97, 13.97, and 8.82%, respectively). Subsequently, full-length VP4, VP2, VP3, and VP1 genes from "Bopivirus B"-positive samples were amplified, cloned, sequenced, and analyzed. The phylogenetic analysis performed on the VP1 genes revealed that the identified bopivirus belonged to genotype B1 (seven strains) and B2 (three strains) and presented a high genetic diversity. Furthermore, a complete genome sequence of a "Bopivirus B" strain (SWUN/B1/2022) was obtained using PCR from fecal sample of a diarrhoeal yeanling. The complete genome was 7,309 nucleotides in length with a standard picornavirus genome organization, and shares 93.10% and 91.10% nucleotide similarity with bopivirus B1 genotype strain ovine/TB14/2010-HUN and bopivirus B2 genotype strain goat/AGK16/2020-HUN, respectively. According to the species classification criteria put forward by the International Committee on Taxonomy of Viruses and VP1 genotype, the strain SWUN/B1/2022 belongs to the bopivirus B1. This strain has unique amino acid substitutions in the VP4, VP2, VP3, and VP1 genes. Moreover, genomic recombination analysis revealed that this strain may be a minor parental strain of bopivirus B1 ovine/TB14/2010-HUN. Evolutionary analysis based on the 2C and 3CD genes revealed that the new bopivirus B1 strain SWUN/B1/2022 presents a unique evolutionary pattern. This study provided evidence to suggest that "Bopivirus B" is circulating with substantial genetic diversity in goats in China at present, and the mixed infection of "Bopivirus B" with other enteric viruses should be considered to be a composite factor in the occurrence of viral diarrhea in goats.

Genetic characteristics of Coxsackievirus A6 from children with hand, foot and mouth disease in Beijing, China, 2017-2019

OBJECT

To investigate the evolution and genetic characteristics of Coxsackievirus A6 (CVA6) which acted as the predominant pathogen of hand, foot and mouth disease (HFMD) in children in Beijing, China, 2017-2019.

METHODS

Throat swab specimens were collected for general Enterovirus (EV), enterovirus A71 (EV-A71) and CVA16 detection by Real-time PCR. These general EV-positive samples were identified by semi-nested RT-PCR method and sequencing. The CVA6 VP1 gene and genome sequences were amplified and sequenced. The phylogenetic, variation and recombination analyses were performed.

RESULTS

A total of 1721 HFMD patients were enrolled in this study, with the male to female ratio of 1.62:1. The majority of cases were less than five years, which accounted for 73.50%. The overall detection rate of EV was 88.32% (1520/1721). A total of 8 EV types were identified, including CVA6 (55.86%), CVA16 (26.32%), EV-A71 (2.24%), CVA10 (2.04%), CVA4 (1.05%), CVA5 (0.59%), CVA2 (0.33%), and CVA8 (0.07%), while 175 (11.51%) EV were untyped. The main pathogen of HFMD was CVA6 from 2017 to 2018, while CVA6 and CVA16 were the main causative pathogens in 2019. The nucleotide and amino acid sequence identities of the 120 CVA6 complete VP1 gene sequences in this study were 91.2%-100.0% and 97.7%-100.0%, respectively. Compared with the prototype strain (Gdula) of CVA6, the nucleotide and amino acid sequence identities were 81.7%-84% and 94.7%-96.3%, respectively. The phylogenetic tree indicated that all 120 CVA6 sequences belonged to sub-genotype D3, while 119 CVA6 sequences belonged to evolutionary branch D3a, except one from 2017 belonged to D3b. Recombination analysis based on the complete genome sequences showed that potential multiple recombination may have occurred in 2B and 3D protein coding regions with EV-A114.

CONCLUSIONS

The main pathogens of HFMD were CVA6 and CVA16 in Beijing, China, 2017-2019; while these CVA6, as recombination strains, belonged to the D3a evolutionary branch.

Molecular Epidemiology Reveals the Co-Circulation of Two Genotypes of Coxsackievirus B5 in China

Coxsackievirus B5 (CVB5) is an important enterovirus B species (EV-Bs) type. We used the full-length genomic sequences of 53 viral sequences from the national hand, foot, and mouth disease surveillance network in the Chinese mainland (2001-2021). Among them, 69 entire VP1 coding region nucleotide sequences were used for CVB5 genotyping and genetic evolution analysis. Phylogenetic analysis based on a data set of 448 complete VP1 sequences showed that CVB5 could be divided into four genotypes (A-D) worldwide. Sequences from this study belonged to genotypes B and D, which dominated transmission in the Chinese mainland. Two transmission lineages of CVB5 have been discovered in the Chinese mainland, lineage 2 was predominant. Markov chain Monte Carlo analysis indicated that the tMRCA of CVB5 in the Chinese mainland could be traced to 1955, while the global trend could be traced to 1862, 93 years earlier than China. The evolution rate of CVB5 was higher in the Chinese mainland than worldwide. The spatiotemporal dynamics analysis of CVB5 assessed that virus transportation events were relatively active in Central, Northeast, North and Northwest China. Recombination analysis revealed frequent intertypic recombination in the non-structural region of CVB5 genotypes B and D with the other EV-Bs, revealing eight recombination lineages. Our study showed the molecular evolution and phylogeography of CVB5 that could provide valuable information for disease prevention.

Contribution of phylogenetics to understanding the evolution and epidemiology of dengue virus

Dengue virus (DENV) is one of the most important arboviral pathogens in the tropics and subtropics, and nearly one‐third of the world's population is at risk of infection. The transmission of DENV involves a sylvatic cycle between nonhuman primates (NHP) and Aedes genus mosquitoes, and an endemic cycle between human hosts and predominantly Aedes aegypti. DENV belongs to the genus Flavivirus of the family Flaviviridae and consists of four antigenically distinct serotypes (DENV‐1‐4). Phylogenetic analyses of DENV have revealed its origin, epidemiology, and the drivers that determine its molecular evolution in nature. This review discusses how phylogenetic research has improved our understanding of DENV evolution and how it affects viral ecology and improved our ability to analyze and predict future DENV emergence.

Molecular characteristics of a coxsackievirus A12 strain in Zhejiang of China, 2019

Background

Enterovirus A (EV-A), such as enterovirus A71 (EV-A71), generally causes hand, foot, and mouth disease (HFMD). However, limited studies focused on uncommon enterovirus serotypes such as coxsackievirus A12 (CV-A12). This study aimed to provide evidence to determine the molecular characteristics of a CV-A12 strain isolated in Zhejiang province, China.

Methods

In routine surveillance of HFMD, we identified a child case with CV-A12 infection in 2019 in Zhejiang province, China. Enterovirus was examined by using real-time reverse transcription-PCR (qRT-PCR). A partial VP1 sequence was amplified to determine the serotype, and then a full-length CV-A12 genome was sequenced. Nucleotide and amino acid similarity was calculated with those CV-A12 strains available in GenBank. Recombination was detected using RDP 4 and SimPlot. Furthermore, phylogenetic analysis was conducted by using BEAST 1.10, and protein modeling was performed with I-TASSER webserver.

Results

A full-length CV-A12 genome PJ201984 was isolated in a Chinese child with HFMD. The similarities with complete coding sequences of the CV-A12 strains in GenBank ranged between 79.3–100% (nucleotide) and 94.4–100% (amino acid), whereas it was 88.7–100.0% (nucleotide) and 97.2–100% (amino acid) when excluding the CV-A12 prototype strain Texas-12. In PJ201984, amino acid variations were more divergent in P2 and P3 regions than those in P1; the majority of those variations in VP1 (13/15) and VP4 (7/8) were similar to those documented in recently isolated CV-A12 strains in China. Furthermore, recombination was identified in P2 region, which involved a CV-A5 strain collected in China. Phylogenetic analysis revealed that PJ201984 clustered together with multiple CV-A12 strains isolated in China and the Netherlands during 2013–2018, as compared to another cluster consisting of CV-A12 strains in China and France during 2009–2015. Additionally, protein models of VP1 and VP4 in PJ201984 were well predicted to be similar to VP1 protein of EV-A71 and VP4 protein of coxsackievirus A21, respectively.

Conclusions

The full-length CV-A12 genome was characterized to have common recombination in P2 region and be phylogenetically related to those CV-A12 strains isolated in recent years, suggesting a continual spread in China. It warrants strengthening the routine surveillance for uncommon enterovirus serotypes, particularly on possible recombination and variation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-022-01892-1.

Genetic and immunological characterization of G9 group A porcine rotaviruses in China

G9 group A rotaviruses (RVAs) are considered emerging pathogens in pigs and humans, and pigs are considered a potential host reservoir for human G9 RVAs. In this study, RVAs of two genotypes, G9P[23] and G9P[13], were successfully isolated and the genomic sequences were obtained, the genome constellation is G9-P[23]-I5-R1-C1-M1-A8-N1-T1-E1-H1 and G9-P[13]-I5-R1-C1-M1-A8-N1-T7-E1-H1 respectively. One strain which amplified from clinic faecal sample had an unique genome constellation G9-P[23]-I1-R1-C1-M1-A8-N1-T1-E1-H1. All the genomic segments of three porcine G9 RVAs were closely related to those of porcine and/or porcine-like human RVAs, demonstrating that the three viruses were porcine-human reassortant strains. To study the immunogenicity of the porcine G9 RVAs, 6-week-old female BALB/c mice were immunized with inactivated vaccines derived from porcine RVAs and then mated. The highest titres of neutralizing antibodies against G9P[23] and G9P[13] porcine RVAs (1,291 ± 35.22 and 1:232 ± 39.28 respectively) were produced in mice 7 days after the second immunization. Suckling mice born to the vaccinated dams were protected by maternal antibodies against challenge with homologous strains. Overall, our data demonstrate the occurrence of porcine-human reassortants of G9 RVAs, and extend our understanding of the immunogenicity of porcine G9 rotaviruses. They also provide a basis for the development of a porcine G9 RVA vaccine.

Recombination in Positive-Strand RNA Viruses

RNA recombination is a major driver of genetic shifts tightly linked to the evolution of RNA viruses. Genomic recombination contributes substantially to the emergence of new viral lineages, expansion in host tropism, adaptations to new environments, and virulence and pathogenesis. Here, we review some of the recent progress that has advanced our understanding of recombination in positive-strand RNA viruses, including recombination triggers and the mechanisms behind them. The study of RNA recombination aids in predicting the probability and outcome of viral recombination events, and in the design of viruses with reduced recombination frequency as candidates for the development of live attenuated vaccines. Surveillance of viral recombination should remain a priority in the detection of emergent viral strains, a goal that can only be accomplished by expanding our understanding of how these events are triggered and regulated.

Enterovirus D68 epidemic, UK, 2018, was caused by subclades B3 and D1, predominantly in children and adults, respectively, with both subclades exhibiting extensive genetic diversity

Enterovirus D68 (EV-D68) has recently been identified in biennial epidemics coinciding with diagnoses of non-polio acute flaccid paralysis/myelitis (AFP/AFM). We investigated the prevalence, genetic relatedness and associated clinical features of EV-D68 in 193 EV-positive samples from 193 patients in late 2018, UK. EV-D68 was detected in 83 (58 %) of 143 confirmed EV-positive samples. Sequencing and phylogenetic analysis revealed extensive genetic diversity, split between subclades B3 (n=50) and D1 (n=33), suggesting epidemiologically unrelated infections. B3 predominated in children and younger adults, and D1 in older adults and the elderly (P=0.0009). Clinical presentation indicated causation or exacerbation of respiratory distress in 91.4 % of EV-D68-positive individuals, principally cough (75.3 %), shortness of breath (56.8 %), coryza (48.1 %), wheeze (46.9 %), supplemental oxygen required (46.9 %) and fever (38.9 %). Two cases of AFM were observed, one with EV-D68 detectable in the cerebrospinal fluid, but otherwise neurological symptoms were rarely reported (n=4). Both AFM cases and all additional instances of intensive care unit (ICU) admission (n=5) were seen in patients infected with EV-D68 subclade B3. However, due to the infrequency of severe infection in our cohort, statistical significance could not be assessed.

Viral Population Diversity during Co-Infection of Foot-And-Mouth Disease Virus Serotypes SAT1 and SAT2 in African Buffalo in Kenya

African buffalo are the natural reservoirs of the SAT serotypes of foot-and-mouth disease virus (FMDV) in sub-Saharan Africa. Most buffalo are exposed to multiple FMDV serotypes early in life, and a proportion of them become persistently infected carriers. Understanding the genetic diversity and evolution of FMDV in carrier animals is critical to elucidate how FMDV persists in buffalo populations. In this study, we obtained oropharyngeal (OPF) fluid from naturally infected African buffalo, and characterized the genetic diversity of FMDV. Out of 54 FMDV-positive OPF, 5 were co-infected with SAT1 and SAT2 serotypes. From the five co-infected buffalo, we obtained eighty-nine plaque-purified isolates. Isolates obtained directly from OPF and plaque purification were sequenced using next-generation sequencing (NGS). Phylogenetic analyses of the sequences obtained from recombination-free protein-coding regions revealed a discrepancy in the topology of capsid proteins and non-structural proteins. Despite the high divergence in the capsid phylogeny between SAT1 and SAT2 serotypes, viruses from different serotypes that were collected from the same host had a high genetic similarity in non-structural protein-coding regions P2 and P3, suggesting interserotypic recombination. In two of the SAT1 and SAT2 co-infected buffalo identified at the first passage of viral isolation, the plaque-derived SAT2 genomes were distinctly grouped in two different genotypes. These genotypes were not initially detected with the NGS from the first passage (non-purified) virus isolation sample. In one animal with two SAT2 haplotypes, one plaque-derived chimeric sequence was found. These findings demonstrate within-host evolution through recombination and point mutation contributing to broad viral diversity in the wildlife reservoir. These mechanisms may be critical to FMDV persistence at the individual animal and population levels, and may contribute to the emergence of new viruses that have the ability to spill-over to livestock and other wildlife species.

OUP accepted manuscript

Chronic enterovirus infections can cause significant morbidity, particularly in immunocompromised patients. This study describes a fatal case associated with a chronic untypeable enterovirus infection in an immunocompromised patient admitted to a Dutch university hospital over nine months. We aimed to identify the enterovirus genotype responsible for the infection and to determine potential evolutionary changes. Long-read sequencing was performed using viral targeted sequence capture on four respiratory and one faecal sample. Phylogenetic analysis was performed using a maximum likelihood method, along with a root-to-tip regression and time-scaled phylogenetic analysis to estimate evolutionary changes between sample dates. Intra-host variant detection, using a Fixed Ploidy algorithm, and selection pressure, using a Fixed Effect Likelihood and a Mixed Effects Model of Evolution, were also used to explore the patient samples. Near-complete genomes of enterovirus C104 (EV-C104) were recovered in all respiratory samples but not in the faecal sample. The recovered genomes clustered with a recently reported EV-C104 from Belgium in August 2018. Phylodynamic analysis including ten available EV-C104 genomes, along with the patient sequences, estimated the most recent common ancestor to occur in the middle of 2005 with an overall estimated evolution rate of 2.97 × 10⁻³ substitutions per year. Although positive selection pressure was identified in the EV-C104 reference sequences, the genomes recovered from the patient samples alone showed an overall negative selection pressure in multiple codon sites along the genome. A chronic infection resulting in respiratory failure from a relatively rare enterovirus was observed in a transplant recipient. We observed an increase in single-nucleotide variations between sample dates from a rapidly declining patient, suggesting mutations are weakly deleterious and have not been purged during selection. This is further supported by the persistence of EV-C104 in the patient, despite the clearance of other viral infections. Next-generation sequencing with viral enrichment could be used to detect and characterise challenging samples when conventional workflows are insufficient.

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

Coxsackievirus B4 sewage-isolate induces pancreatitis after oral infection of mice

Numerous serotypes which belong to the genus Enterovirus (EV) show variability in their virulence and clinical manifestations. They are also known to undergo changes caused by mutations and recombination during their circulation in the environment and the population. Various EV serotypes are prevalent in groundwater, wastewater and surface waters. Our previous studies showed that oral infection induces pancreatitis depending on specific conditions, such as gravidity, in an outbred murine model. Our aim in the present study was to further explore the pancreatic histopathology in an outbred mouse model following oral infection with clinical isolates from a patient who had aseptic meningitis and an isolate from a treated-sewage sample recovered from the residential area of the patient. The isolates were identified as coxsackievirus B4 (CVB4) in tissue culture. The CVB4 sewage-isolate induced pancreatitis after oral infection. In contrast, pancreatitis was absent following infection with the clinical isolates. Comparison of polyprotein sequences showed that the treated-sewage strains differed from the patient's isolates by 9 and 11 amino acids. We conclude that the isolates of clinical and environmental origin differed in their pathogenic properties and showed genetic variation.

The importance of enterovirus surveillance in a post-polio world

Poliovirus is known to most people in the world as the cause of polio, a devastating paralytic disease from the past. Success in polio eradication has understandably translated into stricter containment plans for poliovirus, coordinated by WHO. In this Personal View, we discuss the impact of recent biosafety level 3+ guidelines for handling potential poliovirus-containing diagnostic specimens, which has resulted in closure of many national WHO poliovirus reference laboratories. This reduction in laboratory capacity has a knock-on effect of capability to detect and characterise non-polio enteroviruses in samples obtained from patients with neurological symptoms. The development is of concern given the widespread circulation of non-polio enteroviruses, their role as the most common cause of meningitis worldwide, and their involvement in other severe neurological conditions, such as acute flaccid myelitis and encephalitis. These disease presentations have increased substantially in the past decade, and have been associated with major outbreaks of enterovirus D68 and enterovirus A71, leaving many who survived with lasting paralysis and disabilities. To address this growing gap in diagnostic and surveillance capability, we have established the European Non-Poliovirus Enterovirus Network (also known as ENPEN) as a supra-national, non-commercial, core reference consortium. Our consortium will develop, test, and implement generic surveillance platforms for non-polio enteroviruses and other emerging viral diseases.

Modular Evolution of Coronavirus Genomes

The viral family Coronaviridae comprises four genera, termed Alpha-, Beta-, Gamma-, and Deltacoronavirus. Recombination events have been described in many coronaviruses infecting humans and other animals. However, formal analysis of the recombination patterns, both in terms of the involved genome regions and the extent of genetic divergence between partners, are scarce. Common methods of recombination detection based on phylogenetic incongruences (e.g., a phylogenetic compatibility matrix) may fail in cases where too many events diminish the phylogenetic signal. Thus, an approach comparing genetic distances in distinct genome regions (pairwise distance deviation matrix) was set up. In alpha, beta, and delta-coronaviruses, a low incidence of recombination between closely related viruses was evident in all genome regions, but it was more extensive between the spike gene and other genome regions. In contrast, avian gammacoronaviruses recombined extensively and exist as a global cloud of genes with poorly corresponding genetic distances in different parts of the genome. Spike, but not other structural proteins, was most commonly exchanged between coronaviruses. Recombination patterns differed between coronavirus genera and corresponded to the modular structure of the spike: recombination traces were more pronounced between spike domains (N-terminal and C-terminal parts of S1 and S2) than within domains. The variability of possible recombination events and their uneven distribution over the genome suggest that compatibility of genes, rather than mechanistic or ecological limitations, shapes recombination patterns in coronaviruses.

Coxsackievirus A16 in Southern Vietnam

Background: Hand, Foot and Mouth Disease (HFMD) is a major public health concern in the Asia-Pacific region. Most recent HFMD outbreaks have been caused by enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), CVA10, and CVA6. There has been no report regarding the epidemiology and genetic diversity of CVA16 in Vietnam. Such knowledge is critical to inform the development of intervention strategies.

Materials and Methods: From 2011 to 2017, clinical samples were collected from in- and outpatients enrolled in a HFMD research program conducted at three referral hospitals in Ho Chi Minh City (HCMC), Vietnam. Throat or rectal swabs positive for CVA16 with sufficient viral load were selected for whole genome sequencing and evolutionary analysis.

Results: Throughout the study period, 320 CVA16 positive samples were collected from 2808 HFMD patients (11.4%). 59.4% of patients were male. The median age was 20.8 months (IQR, 14.96–31.41). Patients resided in HCMC (55.3%), Mekong Delta (22.2%), and South East Vietnam (22.5%). 10% of CVA16 infected patients had moderately severe or severe HFMD. CVA16 positive samples from 153 patients were selected for whole genome sequencing, and 66 complete genomes were obtained. Phylogenetic analysis demonstrated that Vietnamese CVA16 strains belong to a single genogroup B1a that clusters together with isolates from China, Japan, Thailand, Malaysia, France and Australia. The CVA16 strains of the present study were circulating in Vietnam some 4 years prior to its detection in HFMD cases.

Conclusion: We report for the first time on the molecular epidemiology of CVA16 in Vietnam. Unlike EV-A71, which showed frequent replacement between subgenogroups B5 and C4 every 2–3 years in Vietnam, CVA16 displays a less pronounced genetic alternation with only subgenogroup B1a circulating in Vietnam since 2011. Our collective findings emphasize the importance of active surveillance for viral circulation in HFMD endemic countries, critical to informing outbreak response and vaccine development.

Genetic diversity and evolution of enterovirus A71 subgenogroup C1 from children with hand, foot, and mouth disease in Thailand

Enterovirus A71 (EV-A71) can cause hand, foot, and mouth disease (HFMD) in children and may be associated with severe neurological complications. There have been numerous reports of increased incidence of EV-A71 subgenogroup C1 (EV-A71 C1) infections associated with neurological diseases since the first occurrence in Germany in 2015. Here, we describe 11 full-length genome sequences of 2019 EV-A71 C1 strains isolated from HFMD patients in Thailand from 2019 to early 2020. The genetic evolution of 2019 EV-A71 C1 was traced in the outbreaks, and the emergence of multiple lineages was detected. Our results demonstrated that 2019 EV-A71 C1 from Thailand emerged through recombination between its nonstructural protein gene and those of other EV-A genotypes. Bayesian-based phylogenetic analysis showed that the 2019 EV-A71 C1 Thai strains share a common ancestor with variants in Europe (Denmark and France). The substitution rate for the 2019 EV-A71 C1 genome was estimated to be 4.38 × 10^-3 substitutions/(site∙year^-1) (95% highest posterior density interval: 3.84-4.94 × 10^-3 substitutions/[site∙year^-1]), approximating that observed between previous EV-A71 C1 outbreaks. These data are essential for understanding the evolution of EV-A C1 during the ongoing HFMD outbreak and may be relevant to disease outcomes in children worldwide.

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

A novel interspecies recombinant enterovirus (Enterovirus A120) isolated from a case of acute flaccid paralysis in China

EV-A120 is a recently identified serotype of the enterovirus A species. Only one full-length genomic sequence is currently available in GenBank, and very few studies have been conducted on EV-A120 globally. Thus, additional information and research on EV-A120 are needed to explore its genetic characteristics, phylogeny, and relationship with enteroviral disease. In this study, we report the phylogenetic characteristics of a EV-A120 strain (Q0082/XZ/CHN/2000) from Tibet, China. The amino acid sequence similarity and nucleotide sequence similarity of the full-length genomic sequence of this EV-A120 strain and the EV-A120 prototype strain were 96.3% and 79.9%, respectively, showing an evolutionary trend. Recombination analysis found intraspecies recombination in the 5′ -UTR, 2B, 2C, and 3D regions. Serum neutralization testing of the EV-A120 (Q0082) strain was also carried out. Low serum-positive rates and geometric mean titres (GMTs) indicated that the extent of EV-A120 transmission and exposure in the population was very limited compared with that in the outbreaks of EV-A71 and CV-A16 in China since 2008. The EV-A120 strain (Q0082) is non-temperature sensitive, indicating its potential to spread in the population. In summary, this study reports the full-length genomic sequence of EV-A120 and provides important information for its global molecular epidemiology.

The First Whole Genome Sequence and Characterisation of Avian Nephritis Virus Genotype 3

Avian nephritis virus (ANV) is classified in the Avastroviridae family with disease associations with nephritis, uneven flock growth and runting stunting syndrome (RSS) in chicken and turkey flocks, and other avian species. The whole genome of ANV genotype 3 (ANV-3) of 6959 nucleotides including the untranslated 5′ and 3′ regions and polyadenylated tail was detected in a metagenomic virome investigation of RSS-affected chicken broiler flocks. This report characterises the ANV-3 genome, identifying partially overlapping open reading frames (ORFs), ORF1a and ORF1b, and an opposing secondary pseudoknot prior to a ribosomal frameshift stemloop structure, with a separate ORF2, whilst observing conserved astrovirus motifs. Phylogenetic analysis of the Avastroviridae whole genome and ORF2 capsid polyprotein classified the first complete whole genome of ANV-3 within Avastroviridae genogroup 2.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.

Enterovirus subtyping in a routine UK laboratory setting between 2013 and 2017

BACKGROUND

Human enteroviruses (EV) are the leading cause of viral meningitis. EV genotyping is predominantly performed through amplification and sequencing of viral capsid protein-1 (VP1), frequently by national reference laboratories (NRLs).

OBJECTIVE

To determine the frequency of genotyping failure in our NRL-submitted samples and apply a superior alternative assay to resolve untyped specimens.

STUDY DESIGN

We initially audited genotyping data received for a cohort of patients in the East Midlands, UK by the NRL between 2013 and 2017, then identified an alternative RT-PCR typing method by literature review and evaluated primers from both assays in silico against comprehensive publicly available genomic data. The alternative assay was further optimised and applied to archived nucleic acids from previously untypable samples.

RESULTS

Genotyping data showed a significant increase in untypable EV strains through the study period (p = 0.0073). Typing failure appeared unrelated to sample type or viral load. In silico analyses of 2,201 EV genomes showed high levels of mismatch between reference assay primers and clinically significant EV-species, in contrast to a selected alternative semi-nested RT-PCR VP1-typing assay. This alternative assay, with minor modifications, successfully genotyped 23 of 24 previously untypable yet viable archived specimens (EV-A, n = 4; EV-B, n = 19). Phylogenetic analyses identified no predominant strain within NRL untypable isolates, suggesting sub-optimal reference assay sensitivity across EV species, in agreement with in silico analyses.

CONCLUSION

This modified highly sensitive RT-PCR assay presents a suitable alternative to the current English national reference VP1-typing assay and is recommended in other settings experiencing typing failure.

Multiple genotypes of enterovirus G carrying a papain-like cysteine protease (PL-CP) sequence circulating on two pig farms in Japan: first identification of enterovirus G10 carrying a PL-CP sequence

Two and three genotypes of enterovirus G (EV-G) carrying a papain-like cysteine protease (PL-CP) sequence were detected on two pig farms and classified into genotypes G1 and G10, and G1, G8, and G17, respectively, based on VP1 sequences. A G10 EV-G virus bearing a PL-CP sequence was detected for the first time. Phylogenetic analysis of the P2 and P3 regions grouped the viruses by farm with high sequence similarity. Furthermore, clear recombination break points were detected in the 2A region, suggesting that PL-CP EV-G-containing strains gained sequence diversity through recombination events among the multiple circulating EV-G genotypes on the farms.

An Extended Primer Grip of Picornavirus Polymerase Facilitates Sexual RNA Replication Mechanisms

Picornaviruses have both asexual and sexual RNA replication mechanisms. Asexual RNA replication mechanisms involve one parental template, whereas sexual RNA replication mechanisms involve two or more parental templates. Because sexual RNA replication mechanisms counteract ribavirin-induced error catastrophe, we selected for ribavirin-resistant poliovirus to identify polymerase residues that facilitate sexual RNA replication mechanisms. We used serial passage in ribavirin, beginning with a variety of ribavirin-sensitive and ribavirin-resistant parental viruses. Ribavirin-sensitive virus contained an L420A polymerase mutation, while ribavirin-resistant virus contained a G64S polymerase mutation. A G64 codon mutation (G64Fix) was used to inhibit emergence of G64S-mediated ribavirin resistance. Revertants (L420) or pseudorevertants (L420V and L420I) were selected from all independent lineages of L420A, G64Fix L420A, and G64S L420A parental viruses. Ribavirin resistance G64S mutations were selected in two independent lineages, and novel ribavirin resistance mutations were selected in the polymerase in other lineages (M299I, M323I, M392V, and T353I). The structural orientation of M392, immediately adjacent to L420 and the polymerase primer grip region, led us to engineer additional polymerase mutations into poliovirus (M392A, M392L, M392V, K375R, and R376K). L420A revertants and pseudorevertants (L420V and L420I) restored efficient viral RNA recombination, confirming that ribavirin-induced error catastrophe coincides with defects in sexual RNA replication mechanisms. Viruses containing M392 mutations (M392A, M392L, and M392V) and primer grip mutations (K375R and R376K) exhibited divergent RNA recombination, ribavirin sensitivity, and biochemical phenotypes, consistent with changes in the fidelity of RNA synthesis. We conclude that an extended primer grip of the polymerase, including L420, M392, K375, and R376, contributes to the fidelity of RNA synthesis and to efficient sexual RNA replication mechanisms.IMPORTANCE Picornaviruses have both asexual and sexual RNA replication mechanisms. Sexual RNA replication shapes picornavirus species groups, contributes to the emergence of vaccine-derived polioviruses, and counteracts error catastrophe. Can viruses distinguish between homologous and nonhomologous partners during sexual RNA replication? We implicate an extended primer grip of the viral polymerase in sexual RNA replication mechanisms. By sensing RNA sequence complementarity near the active site, the extended primer grip of the polymerase has the potential to distinguish between homologous and nonhomologous RNA templates during sexual RNA replication.

Recombination Analysis of Non-Poliovirus Members of the Enterovirus C Species; Restriction of Recombination Events to Members of the Same 3DPol Cluster

Enteroviruses (EVs) are highly prevalent viruses worldwide. Recombination is known to occur frequently in EVs belonging to species Enterovirus A, Enterovirus B, and Enterovirus C. Although many recombinant vaccine-derived poliovirus (VDPV) strains have been reported, our knowledge on recombination in non-polio EVs in the species Enterovirus C is limited. Here, we combined a dataset consisting of 11 newly generated full-length Enterovirus C sequences and 180 publicly available sequences to study recombination dynamics in non-polio EVs. To identify recombination patterns, maximum likelihood phylogenetic trees of different genomic regions were constructed, and segregation analyses were performed. Recombination was observed between members of the same 3DPol cluster, but was rarely observed between members of different clusters. We hypothesize that this restriction may have arisen through their different compartmentalization in respiratory and enteric tracts related to differences in cellular tropisms so that the opportunity to recombine may not be available.

Genetic diversity of enterovirus G detected in faecal samples of wild boars in Japan: identification of novel genotypes carrying a papain-like cysteine protease sequence

The genetic diversity of enterovirus G (EV-G) was investigated in the wild-boar population in Japan. EV-G-specific reverse transcription PCR demonstrated 30 (37.5 %) positives out of 80 faecal samples. Of these, viral protein 1 (VP1) fragments of 20 samples were classified into G1 (3 samples), G4 (1 sample), G6 (2 samples), G8 (4 samples), G11 (1 sample), G12 (7 samples), G14 (1 sample) and G17 (1 sample), among which 11 samples had a papain-like cysteine protease (PL-CP) sequence, believed to be the first discoveries in G1 (2 samples) or G17 (1 sample) wild-boar EV-Gs, and in G8 (2 samples) or G12 (6 samples) EV-Gs from any animals. Sequences of the non-structural protein regions were similar among EV-Gs possessing the PL-CP sequence (PL-CP EV-Gs) regardless of genotype or origin, suggesting the existence of a common ancestor for these strains. Interestingly, for the two G8 and two G12 samples, the genome sequences contained two versions, with or without the PL-CP sequence, together with the homologous 2C/PL-CP and PL-CP/3A junction sequences, which may explain how the recombination and deletion of the PL-CP sequences occured in the PL-CP EV-G genomes. These findings shed light on the genetic plasticity and evolution of EV-G.

Phylogenetic characteristics and molecular epidemiological analysis of novel enterovirus EV-B83 isolated from Tibet, China

Enterovirus B83 (EV-B83) is a new member of the enterovirus B group. Currently, there are only two full-length genomic sequences of EV-B83 in the GenBank database and few VP1 region sequences. The aetiology and epidemiology of EV-B83 is unclear. 24 stool specimens were collected from twelve AFP patients and 298 stool specimens were collected from 298 healthy children in support of polio eradication activities in Tibet in 1999. Two polioviruses (isolated by L20B cell) and one non-polio enterovirus (isolated by RD cell) were isolated from AFP patients and nine polioviruses (isolated by L20B cell) and 90 non-polio enteroviruses (isolated by RD cell) were isolated from health children. Through molecular typing, we confirmed that the six of non-polio enteroviruses belong to EV-B83. The sequence similarity between the VP1 region of the Tibet isolates and that of the EV-B83 prototype strain was 80%. The maximum-likelihood phylogenetic tree of the partial VP1 region in EV-B83 demonstrated that EV-B83 formed four genotypes globally during the evolution process. The six Tibet EV-B83 strains formed the D genotype alone. Recombination analysis of Tibet EV-B83 showed that CV-B4, CV-A9, EV-B80, and EV-B106 may act as recombinant donors in multiple regions. The serum neutralization test showed that the antibody-positive rate was 58.8% and GMT was 1:19.70, which was higher than the previously reported results of EV-B106 and EV-B80. Temperature sensitivity test results showed that the six Tibet EV-B83 strains were temperature-insensitive with stronger virulence and potential infectivity, which was consistent with the results of the serum neutralization test. This study enriched the genome-wide sequence, epidemiological characteristics, and provided basic data for the follow-up study of EV-B83.

Genomic and serologic characterization of enterovirus A71 brainstem encephalitis

Objective

In 2016, Catalonia experienced a pediatric brainstem encephalitis outbreak caused by enterovirus A71 (EV-A71). Conventional testing identified EV in the periphery but rarely in CSF. Metagenomic next-generation sequencing (mNGS) and CSF pan-viral serology (VirScan) were deployed to enhance viral detection and characterization.

Methods

RNA was extracted from the CSF (n = 20), plasma (n = 9), stool (n = 15), and nasopharyngeal samples (n = 16) from 10 children with brainstem encephalitis and 10 children with meningitis or encephalitis. Pathogens were identified using mNGS. Available CSF from cases (n = 12) and pediatric other neurologic disease controls (n = 54) were analyzed with VirScan with a subset (n = 9 and n = 50) validated by ELISA.

Results

mNGS detected EV in all samples positive by quantitative reverse transcription polymerase chain reaction (qRT-PCR) (n = 25). In qRT-PCR-negative samples (n = 35), mNGS found virus in 23% (n = 8, 3 CSF samples). Overall, mNGS enhanced EV detection from 42% (25/60) to 57% (33/60) (p-value = 0.013). VirScan and ELISA increased detection to 92% (11/12) compared with 46% (4/12) for CSF mNGS and qRT-PCR (p-value = 0.023). Phylogenetic analysis confirmed the EV-A71 strain clustered with a neurovirulent German EV-A71. A single amino acid substitution (S241P) in the EVA71 VP1 protein was exclusive to the CNS in one subject.

Conclusion

mNGS with VirScan significantly increased the CNS detection of EVs relative to qRT-PCR, and the latter generated an antigenic profile of the acute EV-A71 immune response. Genomic analysis confirmed the close relation of the outbreak EV-A71 and neuroinvasive German EV-A71. A S241P substitution in VP1 was found exclusively in the CSF.

Evidence for multiple recombination events within foot-and-mouth disease viruses circulating in West Eurasia

Phylogenetic studies on foot-and-mouth disease viruses (FMDVs) circulating in the West Eurasian region have largely focused on the genomic sequences encoding the structural proteins that determine the serotype. The present study has compared near-complete genome sequences of FMDVs representative of the viruses that circulate in this region. The near-complete genome sequences (ca. 7,600 nt) were generated from multiple overlapping RT-PCR products. These amplicons were from FMDVs belonging to serotypes O, A and Asia-1, including members of the O-PanAsia-II and the A-Iran05 lineages, and of Group-II and Group-VII (Sindh-08) within serotype Asia-1, which are currently predominant and widespread in West Eurasia. These new sequences were analysed together with other sequences obtained from GenBank. Comparison of different regions of the FMDVs genomes revealed evidence for multiple, inter-serotypic, recombination events between FMDVs belonging to the serotypes O, A and Asia-1. It is concluded from the present study that dramatic changes in virus sequences can occur in the field through recombination between different FMDV genomes. These analyses provide information about the ancestry of the serotype O, A and Asia-1 FMDVs that are currently circulating within the West Eurasian region.

Parechovirus A Pathogenesis and the Enigma of Genotype A-3

Parechovirus A is a species in the Parechovirus genus within the Picornaviridae family that can cause severe disease in children. Relatively little is known on Parechovirus A epidemiology and pathogenesis. This review aims to explore the Parechovirus A literature and highlight the differences between Parechovirus A genotypes from a pathogenesis standpoint. In particular, the curious case of Parechovirus-A3 and the genotype-specific disease association will be discussed. Finally, a brief outlook on Parechovirus A research is provided.

National Epidemiology and Evolutionary History of Four Hand, Foot and Mouth Disease-Related Enteroviruses in China from 2008 to 2016

Hand, foot and mouth disease (HFMD) is a major public health concern in China. The most predominant enteroviruses that cause HFMD have traditionally been attributed to enterovirus A71 (EVA71) and coxsackievirus A16 (CVA16). Since its first large outbreak in 2008, the dominant HFMD pathogens are constantly changing. In 2013 and 2015, CVA6 exceeded both EVA71 and CVA16 to become the leading cause of HFMD in some provinces. However, there still lacks a comprehensive overview on the molecular epidemiology and evolution of HFMD-related enteroviruses at the national level. In this study, we performed systematic epidemiological analyses of HFMD-related enteroviruses using the data of 64 published papers that met the inclusion criteria, and conducted phylogenetic analyses based on 12,080 partial VP1 sequences identified in China before 31st June 2018. We found that EVA71 prevalence has decreased sharply but other enteroviruses have increased rapidly from 2008 to 2016 and that one subtype of each enterovirus is represented during the epidemic. In addition, four genotypes EVA71_C4, CVA16_B1, CVA6_D and CVA10_C are the most predominant enterovirus strains and collectively they cause over 90% of all HFMD cases in China according to the phylogenetic trees using representative partial VP1 sequences. These four major enterovirus genotypes have different geographical distributions, and they may co-circulate with other genotypes and serotypes. These results suggest that more molecular epidemiological studies should be performed on several enteroviruses simultaneously, and such information should have implications for virological surveillance, disease management, vaccine development and policy-making on the prevention and control of HFMD.

Recombination in Enteroviruses, a Multi-Step Modular Evolutionary Process

RNA recombination is a major driving force in the evolution and genetic architecture shaping of enteroviruses. In particular, intertypic recombination is implicated in the emergence of most pathogenic circulating vaccine-derived polioviruses, which have caused numerous outbreaks of paralytic poliomyelitis worldwide. Recent experimental studies that relied on recombination cellular systems mimicking natural genetic exchanges between enteroviruses provided new insights into the molecular mechanisms of enterovirus recombination and enabled to define a new model of genetic plasticity for enteroviruses. Homologous intertypic recombinant enteroviruses that were observed in nature would be the final products of a multi-step process, during which precursor nonhomologous recombinant genomes are generated through an initial inter-genomic RNA recombination event and can then evolve into a diversity of fitter homologous recombinant genomes over subsequent intra-genomic rearrangements. Moreover, these experimental studies demonstrated that the enterovirus genome could be defined as a combination of genomic modules that can be preferentially exchanged through recombination, and enabled defining the boundaries of these recombination modules. These results provided the first experimental evidence supporting the theoretical model of enterovirus modular evolution previously elaborated from phylogenetic studies of circulating enterovirus strains. This review summarizes our current knowledge regarding the mechanisms of recombination in enteroviruses and presents a new evolutionary process that may apply to other RNA viruses.

The high genetic similarity between rhinoviruses and enteroviruses remains as a pitfall for molecular diagnostic tools: A three-year overview

BACKGROUND

Enteroviruses (EVs) and rhinoviruses (RVs) belong to the Enterovirus genus within the Picornaviridae family, and show genetic similarities. These viruses are related to mild diseases, but EVs infections can sometimes lead to more severe complications. Current diagnostic molecular techniques should discriminate between the four EV and the three RV species that infect humans. The aim was to revise the EV and RV PCR-confirmed specimens by sequencing for genetic characterisation.

MATERIAL AND METHODS

Respiratory tract specimens were collected from patients with suspicion of respiratory infection. Respiratory viruses' laboratory-confirmation was performed by commercial multiplex real-time RT-PCR assays. Genetic characterisation of all EV and in a selection of RV was performed based on the phylogenetic analyses of partial VP1 and VP4/2 sequences, respectively.

RESULTS

From 19,957 tested specimens, 309 (1.5%) were EV-positive, 2546 (12%) were RV-positive, and 233 (1%) were EV/RV co-detections. The phylogenetic analyses revealed that: among single EV detections, 177/309 (57%) were characterised as EV, 2/309 (1%) as RV, and 130/309 (42%) could not be typed; among single 1771 RV detections (Ct < 35), 1651/1771 (93%) were characterised as RV, 3/1771 (0.3%) as EV and 117/1771 (6.7%) could not be typed. Among EV/RV co-detections, 62/233 (27%) were characterised as EV, 130/233 (56%) as RV and 41/233 (18%) could not be typed.

CONCLUSIONS

A diagnostic method well considered for routine laboratory-confirmation of respiratory viruses should discriminate EV and RV targets. RVs are usually associated with mild respiratory disease, but the potential relatedness of EVs to neurological complications makes their monitoring mandatory. Therefore, an accurate detection and differentiation should be required in commercial diagnostic solutions.

A novel defective recombinant porcine enterovirus G virus carrying a porcine torovirus papain-like cysteine protease gene and a putative anti-apoptosis gene in place of viral structural protein genes

Enterovirus G (EV-G) belongs to the family of Picornaviridae. Two types of recombinant porcine EV-Gs carrying papain-like cysteine protease (PLCP) gene of porcine torovirus, a virus in Coronaviridae, are reported. Type 1 recombinant EV-Gs are detected in pig feces in Japan, USA, and Belgium and carry the PLPC gene at the junction site of 2C/3A genes, while PLPC gene replaces the viral structural genes in type 2 recombinant EV-G detected in pig feces in a Chinese farm. We identified a novel type 2 recombinant EV-G carrying the PLCP gene with flanking sequences in place of the viral structural genes in pig feces in Japan. The ~0.3 kb-long upstream flanking sequence had no sequence homology with any proteins deposited in GenBank, while the downstream ~0.9 kb-long flanking sequence included a domain having high amino acid sequence homology with a baculoviral inhibitor of apoptosis repeat superfamily. The pig feces, where the novel type 2 recombinant EV-G was detected, also carried type 1 recombinant EV-G. The amount of type 1 and type 2 recombinant EV-G genomes was almost same in the pig feces. Although the phylogenetic analysis suggested that these two recombinant EV-Gs have independently evolved, type 1 recombinant EV-G might have served as a helper virus by providing viral structural proteins for dissemination of the type 2 recombinant EV-G.

•

A novel type 2 recombinant EV-G was discovered in pig feces in Japan.

•

Type 2 recombinant EV-G carried the PLCP torovirus gene with unknown flanking genes, in place of the viral structural proteins.

•

Amount of type 2 recombinant EV-G in the pig feces was almost same with type 1.

•

Type 2 recombinant EV-G belonged to be a different cluster from the cluster of type 1.

Identification of three linear B cell epitopes using monoclonal antibodies against bovine enterovirus VP2 protein

Bovine enterovirus (BEV) VP2 protein is a structural protein that plays an important role in inducing protective immunity in the host. The function of VP2 has been characterized, but there is little information on its B cell epitopes. Three monoclonal antibodies (mAbs) directed against BEV VP2 were generated and characterized from mice immunized with the recombinant VP2 protein. Three minimal linear epitopes ¹⁵²FQEAFWLEDG¹⁶¹, ¹⁶⁸LIYPHQ¹⁷³, and ⁴⁶DATSVD⁵¹ reactive to the three mAbs were identified using western blotting analysis. Three-dimensional model of the BEV-E virion and the VP2 monomer showed that epitope ¹⁵²FQEAFWLEDG¹⁶¹ is exposed on surface of the virion and epitopes ⁴⁶DATSVD⁵¹ and ¹⁶⁸LIYPHQ¹⁷³ are located inside the virion. Alignment of the amino acid sequences corresponding to the regions containing the three minimal linear epitopes in the VP2 proteins and their cross-reactivity with the three mAbs showed that epitope ¹⁶⁸LIYPHQ¹⁷³ is completely conserved in all BEV strains. Epitope ⁴⁶DATSVD⁵¹ is highly conserved among BEV-E strains and partly conserved among BEV-F strains. However, epitope ¹⁵²FQEAFWLEDG¹⁶¹ is not conserved among BEV-F strains. Using the mAbs of 3H4 and 1E10, we found that VP2 localized in the cytoplasm during viral replication and could be used to monitor the viral antigen in infected tissues using immunohistochemistry. A preliminary 3H4-epitope-based indirect ELISA allowed us to detect anti-BEV-strain-HY12 antibodies in mice. This study indicates that the three mAbs could be useful tools for investigating the structure and function of the viral VP2 protein and the development of serological diagnostic techniques for BEV infection.

Strength in numbers: Mechanisms of viral co-infection

RNA virus populations are diverse due to a variety of factors, including lack of proofreading of the viral RNA-dependent RNA polymerase. These diverse viral populations include defective viruses incapable of productive infection. Recent studies have determined the existence of several modes of viral transmission outside of canonical pathways, including en bloc transmission of multiple viruses into a single host cell via membrane vesicles. Additionally, it has recently been determined that viral aggregation and bacteria can facilitate the delivery of multiple viruses to a single cell. Co-infection of RNA viruses is important since it has the potential to enhance viral fitness. Furthermore, through complementation and recombination, co-infection could potentially promote "resurrection" of otherwise defective viral genomes and has the potential to expand viral diversity.