Recombination strategies and evolutionary dynamics of the Human enterovirus A global gene pool

A newly developed whole genome sequencing protocol enables early tracking of enterovirus D68 molecular evolution

BACKGROUND

Human enterovirus D68 (EV-D68) has been associated with an increase in mild-to-severe pediatric respiratory diseases worldwide. The rate of circulation of this virus is largely underestimated in the population and genetic evolutionary data are usually available only for partial sequences. To achieve a timely genomic surveillance, a reliable, high-throughput EV-68 sequencing assay is required. Here we report an improved high-throughput EV-D68 whole-genome sequencing assay performed directly on clinical samples that is suitable for short-read sequencing platforms. Between June and December 2022, a total 37 (1.9 %) respiratory samples were EV-D68 positive and together with 52 additional samples with a median cycle of quantification (Cq) of 28.3, ranging from 18 to 36.8 Cq were included in the validation analyses. Overall, all the primers had good performance and no mismatches were detected in more than 85 % of sequences (932 whole-genome dataset). Using a cut-off of Cq < 32 in at least 85.5 % of samples a whole-genome or partial genome was obtained, confirming an acceptable positive sequencing rate for the designed method. A total of 65 whole-genome sequences were obtained and have a mean coverage of 98.4 % across the genome, with a median depth of 6158x (range 2815x-7560x). Based on the obtained data, this method is cost effective resulting in an easy-to-perform protocol helpful for tracing the evolution of EV-D68 in protein different from VP1. EV-D68 could become a significant pathogen for public health in the next future, and thus this protocol for whole genome sequencing could help clinical and molecular virologists to be ready for molecular epidemiology surveillance.

Mamastrovirus species are shaped by recombination and can be reliably distinguished in ORF1b genome region

Astroviruses are a diverse group of small non-enveloped positive sense single-stranded RNA viruses that infect animals and birds. More than half of all known genome sequences of mammalian astroviruses are not assigned to provisional species, and the biological mechanisms that could support segregation of astroviruses into species are not well understood. The systematic analysis of recombination in Mamastrovirus genomes available in GenBank was done to identify mechanisms providing genetic distinction between astroviruses. Recombination breakpoints were present in all Mamastrovirus genome regions, but occurred most commonly at the ORF1b/ORF2 junction. Recombination was ubiquitous within, but never between established and putative new species, and may be suggested as an additional species criterion. The current species criterion for the genus Mamastrovirus based on ORF2 amino acid sequence p-distances did not reliably distinguish several established species and was of limited use to identify distinct groups among unclassified astroviruses that were isolated recently, predominantly from cattle and pigs. A 17% nucleotide sequence distance cut-off in ORF1b fairly distinguished the established species and several groups among the unclassified viruses, providing better correspondence between phylogenetic grouping, reproductive isolation and the virus hosts. Sequence distance criteria (17% in nucleotide sequence of ORF1b and 25% in amino acid sequence of ORF2) and the recombination pattern corresponded fairly well as species criteria, but all had minor exclusions among mammalian astroviruses. A combination of these taxonomic criteria supported the established Mamastrovirus species and suggested redefining a few provisional species that were proposed earlier and introducing at least six novel species among recently submitted rat and bovine astroviruses.

Genetic diversity and spread of recombinant coxsackievirus A4 in hand, foot, and mouth disease cases in Bangkok, Thailand: 2017-2023

Coxsackievirus A4 (CVA4) has recently become one of the most common causative agents of hand, foot, and mouth disease. The current study investigated the genetic diversity and spread of recombinant CVA4 by analyzing circulating genotypes and recombinant strains in Bangkok, Thailand, from 2017 to 2023. Partial VP1, 3Dpol, and whole genome sequencing of CVA4 samples collected from collaborating hospitals were conducted. Phylogenetic analysis of CVA4 VP1 and 3Dpol genome regions revealed discordance, indicating recombination. The predominant CVA4 genotype was C3, primarily observed in 2019. The predominant genotype in 2017 was C1. D2, commonly found in China, was occasionally observed. In nucleotide similarity analysis, intertypic recombination between CVA4 and EV-A during the evolutionary history of the virus was evident, particularly in the nonstructural region. The estimated emergence of genotypes C1 and C3 in Thailand occurred around 2014, with an evolutionary rate of 5.8 × 10- 3 nucleotide substitutions per site per year. Genotype D2 exhibited notable variability across both the entire genome and the structural protein region compared to genotype C. Monitoring the genetic diversity and circulation of recombinant CVA4 is crucial for identifying newly emerging virus strains, enabling prompt public health responses and containment efforts, and enhancing surveillance in Thailand.

Genotyping and phylogeographic dynamics of coxsackievirus A16

Coxsackievirus A16 (CV-A16) is one of the major pathogens of Hand, Foot and Mouth disease. Here, we analyzed 287 full-length genome sequences of CV-A16 found worldwide from 1994 to 2019 to see the genomic evolution characteristics. Full-length genome-based phylogenetic tree divided the viruses into five different genotypes, G-a to G-e. The CV-A16 strains circulating in China dominate G-a and G-c, but can also be found in other genotypes including G-b and G-e. Phylogeographic analysis showed a high diversity of CV-A16 distribution. In addition, recombination was shown to drive the genomic evolution of CV-A16 during past decades. However, the structural proteins still remain relative conserved while there is extensive genomic recombination. This study updates the phylogenetic and phylogeographic information of CV-A16 and provides insights into the genetic characteristics of CV-A16.

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.

Epidemiology of hand, foot and mouth disease and genomic surveillance of coxsackievirus A10 circulating in Zhejiang Province, China during 2017 to 2022

BACKGROUND

Coxsackievirus A10 (CA10) is one of the etiological agents associated with hand, foot and mouth disease (HFMD).

OBJECTIVES

We aimed to perform a retrospective analysis of the molecular epidemiological characteristics and genetic features of HFMD associated with CA10 infections in Zhejiang Province from 2017 to 2022.

STUDY DESIGN

Epidemiologic features were summarized. Throat swab specimens were collected and tested. The VP1 regions were sequenced for genotyping. CA10 positive samples were isolated. Whole genomes of CA10 isolations were sequenced. Nucleotide and amino acid changes were characterized. Phylogenetic trees were constructed.

RESULTS

The number of HFMD cases fluctuated from 2017 to 2022. Children aged below 3 years accounted for the majority (66.29%) and boys were more frequently affected than girls. Cases peaked in June. The positivity rate of HEV was 62.69%. A total of 90 strains of CA10 were isolated and 53 genomes were obtained. All CA10 in this study could be assigned to two genogroups, C (C2) and F (F1 and F3).

CONCLUSION

The clinical manifestations of HFMD associated with HEV are complex and diverse. CA10 infection may be emerging as a new and major cause of HFMD because an upward trend was observed in the proportion of CA10 cases after the use of EV71 vaccines. Different genogroups of CA10 had different geographic distribution patterns. Surveillance should be strengthened and further comprehensive studies should be continued to provide a scientific basis for HFMD prevention and control.

Molecular Characteristics of Enterovirus B83 Strain Isolated from a Patient with Acute Viral Myocarditis and Global Transmission Dynamics

This study determined the global genetic diversity and transmission dynamics of enterovirus B83 (EV-B83) and proposed future disease surveillance directions. Blood samples were collected from a patient with viral myocarditis, and viral isolation was performed. The complete genome sequence of the viral isolate was obtained using Sanger sequencing. A dataset of 15 sequences (from three continents) that had sufficient time signals for Bayesian phylogenetic analysis was set up, and the genetic diversity and transmission dynamics of global EV-B83 were analyzed using bioinformatics methods, including evolutionary dynamics, recombination event analysis, and phylogeographic analysis. Here, we report the complete genome sequence of an EV-B83 strain (S17/YN/CHN/2004) isolated from a patient with acute viral myocarditis in Yunnan Province, China. All 15 EV-B83 strains clustered together in a phylogenetic tree, confirming the classification of these isolates as a single EV type, and the predicted time for the most recent common ancestor appeared in 1998. Recombinant signals were detected in the 5'-untranslated region and 2A-3D coding regions of the S17 genome. The phylogeographic analysis revealed multiple intercontinental transmission routes of EV-B83. This study indicates that EV-B83 is globally distributed. Our findings add to the publicly available EV-B83 genomic sequence data and deepen our understanding of EV-B83 epidemiology.

Laboratory diagnosis of nonpolio enteroviruses: A review of the current literature

Infections by nonpolio enteroviruses (EVs) are highly prevalent, particularly among children and neonates, where they may cause substantial morbidity and mortality. Laboratory diagnosis of these viral infections is important in patient prognosis and guidance of clinical management. Although the laboratory diagnosis of nonpolio EVs is mainly based on molecular techniques, classical virus-isolation techniques are still used in reference laboratories. Other techniques, such as antigen detection and serology, are becoming obsolete and rarely used in diagnosis. An important part of diagnosis and surveillance of EV infections is viral typing by VP1 gene sequencing using conventional Sanger technique and more recently, full-genome next-generation sequencing. The latter allows the typing of all EVs, better investigation of EV outbreaks, detection of coinfection, and identification of severity markers in the EV genome.

Current status of hand-foot-and-mouth disease

Hand-foot-and-mouth disease (HFMD) is a viral illness commonly seen in young children under 5 years of age, characterized by typical manifestations such as oral herpes and rashes on the hands and feet. These symptoms typically resolve spontaneously within a few days without complications. Over the past two decades, our understanding of HFMD has greatly improved and it has received significant attention. A variety of research studies, including epidemiological, animal, and in vitro studies, suggest that the disease may be associated with potentially fatal neurological complications. These findings reveal clinical, epidemiological, pathological, and etiological characteristics that are quite different from initial understandings of the illness. It is important to note that HFMD has been linked to severe cardiopulmonary complications, as well as severe neurological sequelae that can be observed during follow-up. At present, there is no specific pharmaceutical intervention for HFMD. An inactivated Enterovirus A71 (EV-A71) vaccine that has been approved by the China Food and Drug Administration (CFDA) has been shown to provide a high level of protection against EV-A71-related HFMD. However, the simultaneous circulation of multiple pathogens and the evolution of the molecular epidemiology of infectious agents make interventions based solely on a single agent comparatively inadequate. Enteroviruses are highly contagious and have a predilection for the nervous system, particularly in child populations, which contributes to the ongoing outbreak. Given the substantial impact of HFMD around the world, this Review synthesizes the current knowledge of the virology, epidemiology, pathogenesis, therapy, sequelae, and vaccine development of HFMD to improve clinical practices and public health efforts.

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.

Enterovirus characterized from cerebrospinal fluid in a cohort from the Eastern United States

BACKGROUND

Enteroviruses (EVs) are predominant causes of a spectrum of neurological diseases. To better understand the origins of the outbreaks of disease associated with EV, it is essential to develop an efficient surveillance system that identifies the circulating EVs and correlate their genomic evolution with the disease presentations.

METHODS

The clinical presentations of patients with positive EV from cerebrospinal fluid (CSF) between 2014 and 2022, diagnosed at the Johns Hopkins Medical Microbiology Laboratory, were compared from year to year. EV typing and whole genome sequencing were performed and correlated to the spectrum of disease.

RESULTS

A total of 95 CSF specimens were positive for EV between 2014 and 2022. The percentage positivity ranged from the lowest of 1.1% in 2020 to the highest of 3.2% in 2015. The median ages declined from 22 years in 2014 to less than one year starting in 2016 to 34 in 2022. Typing using VP1 sequencing revealed that E30 and E6 were associated with meningitis in adults but coxsackieviruses (CVs-B3 and B5) were detected from pediatric patients with fever. Whole genome sequencing revealed multiple recombination events. In 2020, a recombinant CV-A9 was detected in a CSF sample associated with unusual presentation of sepsis, profound acute bilateral sensory neural hearing loss, and myofasciitis.

CONCLUSIONS

EV genomic surveillance is needed for a better understanding of the genetic determinants of neurovirulence. Whole genome sequencing can reveal recombination events missed by traditional molecular surveillance methods.

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.

Case report: Clinical and virological characteristics of aseptic meningitis caused by a recombinant echovirus 18 in an immunocompetent adult

Echovirus 18 has been recognized as an important causative pathogen of aseptic meningitis in young children worldwide, and echovirus 18-induced meningitis is rarely found in adults with immunocompetence. In this case study, we report the clinical and virological characteristics of aseptic meningitis caused by recombinant echovirus 18 in an adult with immunocompetence. A 31-year-old woman with immunocompetence was admitted to our hospital with fever, dizziness, severe headache, nausea, and vomiting for the past 1 day and was diagnosed with viral meningitis based on the clinical manifestations and laboratory results from cerebrospinal fluid (CSF). The patient received antiviral treatment with ribavirin and interferon as soon as the enterovirus infection was identified using qRT-PCR and was cured after 4 days. From the oropharyngeal swab and CSF samples, two echovirus 18 strains were isolated with a single nucleotide difference located at the 5' UTR. Phylogenetic analyses based on the VP1 gene showed that the two strains belonged to the subgenotype C2 and were clustered with sequences obtained from China after 2015, while the results from the 3D polymerase region showed that the two strains were closely related to the E30 strains. Bootscanning results using the 5' UTR to 2A region and the 2B to 3' UTR region showed that potential intertypic recombination had occurred in the 2B gene. Recombination analyses further confirmed that the two strains (echovirus 18) presented genome recombination with echovirus 30 in the nucleotide regions of the 2B gene. To the best of our knowledge, this is the first report of echovirus 18-induced meningitis in an adult with immunocompetence from mainland China, highlighting the need for close surveillance of echovirus 18 both in children and adults in the future.

Enterovirus A Shows Unique Patterns of Codon Usage Bias in Conventional Versus Unconventional Clade

Enterovirus A (EV-A) species cause hand, foot and mouth disease (HFMD), threatening the health of young children. Understanding the mutual codon usage pattern of the virus and its host(s) has fundamental and applied values. Here, through examining multiple codon usage parameters, we found that the codon usage bias among EV-A strains varies and is clade-specific. EVA76, EVA89, EVA90, EVA91 and EVA92, the unconventional clade of EV-A strains, show unique codon usage pattern relative to the two conventional clades, including EVA71, CVA16, CVA6 and CVA10, etc. Analyses of Effective Number of Codon (ENC), Correspondence Analysis (COA) and Parity Rule 2 (PR2), etc., revealed that the codon usage patterns of EV-A strains are shaped by mutation pressure and natural selection. Based on the neutrality analysis, we determined the dominant role of natural selection in the formation of the codon usage bias of EV-A. In addition, we have determined the codon usage compatibility of potential hosts for EV-A strains using codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses, and found that EV-A showed host-specific codon adaptation patterns in different clades. Finally, we confirmed that the unique codon usage pattern of the unconventional clade affected protein expression level in human cell lines. In conclusion, we identified novel characteristics of codon usage bias in distinct EV-A clades associated with their host range, transmission and pathogenicity.

Global emergence of Enterovirus 71: a systematic review

Background

Hand, foot, and mouth disease (HFMD) is a viral infection caused by a virus from the enterovirus genus of picornavirus family that majorly affects children. Though most cases of HFMD do not cause major problems, the outbreaks of Enterovirus 71 (EV71) can produce a high risk of neurological sequelae, including meningoencephalitis, lung difficulties, and mortality. In Asia, HFMD caused by EV71 has emerged as an acutely infectious disease of highly pathogenic potential, which demands the attention of the international medical community.

Main body of the abstract

Some online databases including NCBI, PubMed, Google Scholar, ProQuest, Scopus, and EBSCO were also accessed using keywords relating to the topic for data mining. The paid articles were accessed through the Centre Library facility of Siksha O Anusandhan University. This work describes the structure, outbreak, molecular epidemiology of Enterovirus 71 along with different EV71 vaccines. Many vaccines have been developed such as inactivated whole-virus live attenuated, subviral particles, and DNA vaccines to cure the patients. In Asia–Pacific nations, inactivated EV71 vaccination still confronts considerable obstacles in terms of vaccine standardization, registration, price, and harmonization of pathogen surveillance and measurements.

Short conclusion

HFMD has emerged as a severe health hazard in Asia–Pacific countries in recent decades. In Mainland China and other countries with high HFMD prevalence, the inactivated EV71 vaccination will be a vital tool in safeguarding children's health. When creating inactivated EV71 vaccines, Mainland China ensured maintaining high standards of vaccine quality. The Phase III clinical studies were used to confirm the safety and effectiveness of vaccinations.

Graphical Abstract

Coxsackievirus A6 Recombinant Subclades D3/A and D3/H Were Predominant in Hand-Foot-And-Mouth Disease Outbreaks in the Paediatric Population, France, 2010–2018

Coxsackievirus A6 (CVA6) emerged as the most common enterovirus of seasonal outbreaks of hand-foot-and-mouth disease (HFMD). We investigated CVA6 genetic diversity among the clinical phenotypes reported in the paediatric population during sentinel surveillance in France between 2010 and 2018. CVA6 infection was confirmed in 981 children (mean age 1.52 years [IQR 1.17–2.72]) of whom 564 (58%) were males. Atypical HFMD was reported in 705 (72%) children, followed by typical HFMD in 214 (22%) and herpangina in 57 (6%) children. Throat specimens of 245 children were processed with a target-enrichment new-generation sequencing approach, which generated 213 complete CVA6 genomes. The genomes grouped within the D1 and D3 clades (phylogeny inferred with the P1 genomic region). In total, 201 genomes were classified among the recombinant forms (RFs) A, B, F, G, H, and N, and 12 genomes were assigned to 5 previously unreported RFs (R–V). The most frequent RFs were A (58%), H (19%), G (6.1%), and F (5.2%). The yearly number of RFs ranged between 1 (in 2012 and 2013) and 6 (2018). The worldwide CVA6 epidemic transmission began between 2005 and 2007, which coincided with the global spread of the recombinant subclade D3/RF-A.

Robust AAV Genotyping Based on Genetic Distances in Rep Gene That Are Maintained by Ubiquitous Recombination

Adeno-associated viruses (AAVs) are a convenient tool for gene therapy delivery. According to the current classification, they are divided into the species AAV A and AAV B within the genus Dependoparvovirus. Historically AAVs were also subdivided on the intraspecies level into 13 serotypes, which differ in tissue tropism and targeted gene delivery capacity. Serotype, however, is not a universal taxonomic category, and their assignment is not always robust. Cross-reactivity has been shown, indicating that classification could not rely on the results of serological tests alone. Moreover, since the isolation of AAV4, all subsequent AAVs were subdivided into serotypes based primarily on genetic differences and phylogenetic reconstructions. An increased interest in the use of AAV as a gene delivery tool justifies the need to improve the existing classification. Here, we suggest genotype-based AAV classification below the species level based on the rep gene. A robust threshold was established as 10% nt differences within the 1248 nt genome fragment, with 4 distinct AAV genotypes identified. This distinct sub-species structure is maintained by ubiquitous recombination within, but not between, rep genes of the suggested genotypes.

Genomic surveillance of coxsackievirus A10 reveals genetic features and recent appearance of genogroup D in Shanghai, China, 2016–2020

Coxsackievirus A10 (CVA10) is one of the major causative agents of hand, foot and mouth disease (HFMD). To investigate the epidemiological characteristics as well as genetic features of CVA10 currently circulating in Shanghai, China, we collected a total of 9,952 sporadic HFMD cases from January 2016 to December 2020. In the past five years, CVA10 was the fourth prevalent causatives associated with HFMD in Shanghai and the overall positive rate was 2.78%. The annual distribution experienced significant fluctuations over the past five years. In addition to entire VP1 sequencing, complete genome sequencing and recombination analysis of CVA10 isolates in Shanghai were further performed. A total of 64 near complete genomes and 11 entire VP1 sequences in this study combined with reference sequences publicly available were integrated into phylogenetic analysis. The CVA10 sequences in this study mainly belonged to genogroup C and presented 91%–100% nucleotide identity with other Chinese isolates based on VP1 region. For the first time, our study reported the appearance of CVA10 genogroup D in Chinese mainland, which had led to large-scale outbreaks in Europe previously. The recombination analysis showed the recombination break point located between 5,100 nt and 6,700 nt, which suggesting intertypic recombination with CVA16 genogroup D. To conclusion, CVA10 genogroup C was the predominant genogroup in Shanghai during 2016–2020. CVA10 recombinant genogroup D was firstly reported in circulating in Chinese mainland. Continuous surveillance is needed to better understand the evolution relationships and transmission pathways of CVA10 to help to guide disease control and prevention.

•

Systematic profiles of genetic features of CVA10 near complete genome.

•

First report of the appearance of CVA10 genogroup D in Chinese mainland.

•

Genomic comparisons indicate the potential recombinant origin of CVA10 genogroup D.

Global spatiotemporal transmission patterns of human enterovirus 71 from 1963 to 2019

Enterovirus 71 (EV71) can cause large outbreaks of hand, foot, and mouth disease (HFMD) and severe neurological diseases, which is regarded as a major threat to public health, especially in Asia-Pacific regions. However, the global spatiotemporal spread of this virus has not been identified. In this study, we used large sequence datasets and a Bayesian phylogenetic approach to compare the molecular epidemiology and geographical spread patterns of different EV71 subgroups globally. The study found that subgroups of HFMD presented global spatiotemporal variation, subgroups B0, B1, and B2 have caused early infections in Europe and America, and then subgroups C1, C2, C3, and C4 replaced B0-B2 as the predominant genotypes, especially in Asia-Pacific countries. The dispersal patterns of genotype B and subgroup C4 showed the complicated routes in Asia and the source might in some Asian countries, while subgroups C1 and C2 displayed more strongly supported pathways globally, especially in Europe. This study found the predominant subgroup of EV71 and its global spatiotemporal transmission patterns, which may be beneficial to reveal the long-term global spatiotemporal transmission patterns of human EV71 and carry out the HFMD 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.

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.

Genetic analysis of Enterovirus D68 associated with pneumonia in children from South India

EV-D68 is an emerging enterovirus infection associated with severe acute respiratory illness (SARI), acute flaccid myelitis (AFM) and acute flaccid paralysis (AFP). While EV-D68 outbreaks and sporadic cases are reported globally, a single case has been reported from India. The present study aims to investigate the molecular epidemiology and clinical characteristics of EV-D68-associated SARI cases from South India. We screened influenza-negative archived throat swab specimens from Influenza-Like Illness (ILI) and SARI cases (n=959; 2016 to 2018 period) for enteroviruses by pan-enterovirus real-time RT-PCR. Thirteen samples positive for enteroviruses were typed by PCR and sequencing based on VPI, VP2 and/or 5'NCR regions. One EV-D68 RNA sample was subjected to next-generation sequencing for whole genome characterisation. Among 13 enterovirus cases, four were ECHO-11, three EV-D68, two CV-A16 and one each EV-71, CV-B1, CV-B2 and CV-A9. All three cases of EV-D68 infection were reported in children below 2 years of age from Kerala state of South India during June and July 2017. The patients developed pneumonia without any neurological complications. Sequencing based on VPI and 5'NCR regions showed that EV-D68 strains belong to the novel subclade B3. The EV-D68 complete genome identified with two unique amino acid substitutions in VP1 (T-246-I) and 3D (K-344-R) regions. This study reiterates the EV-D68 novel subclade B3 circulation in India and indicates the urgent need for structured EV-D68 surveillance in the country to describe the epidemiology.

Rules governing genetic exchanges among viral types from different Enterovirus A clusters

The species Enterovirus A (EV-A) consists of two conventional clusters and one unconventional cluster. At present, sequence analysis shows no evidence of recombination between conventional and unconventional EV-A types. However, the factors underlying this genetic barrier are unclear. Here, we systematically dissected the genome components linked to these peculiar phenomena, using the viral reverse genetic tools. We reported that viral capsids of the unconventional EV-A types expressed poorly in human cells. The trans-encapsidation outputs across conventional and unconventional EV-A types were also with low efficiency. However, replicons of conventional types bearing exchanged 5'-untranslated region (UTR) or non-structural regions from the unconventional types were replication-competent. Furthermore, we created a viable recombinant EVA71 (conventional type) with its P3 region replaced by that from EVA89 (unconventional type). Thus, our data for the first time reveal the potential for fertile genetic exchanges between conventional and unconventional EV-A types. It also discloses that the mysterious recombination barriers may lie in uncoordinated capsid expression and particle assembly by different EV-A clusters.

Epidemiology and Sequence-Based Evolutionary Analysis of Circulating Non-Polio Enteroviruses

Enteroviruses (EVs) are positive-sense RNA viruses, with over 50,000 nucleotide sequences publicly available. While most human infections are typically associated with mild respiratory symptoms, several different EV types have also been associated with severe human disease, especially acute flaccid paralysis (AFP), particularly with endemic members of the EV-B species and two pandemic types—EV-A71 and EV-D68—that appear to be responsible for recent widespread outbreaks. Here we review the recent literature on the prevalence, characteristics, and circulation dynamics of different enterovirus types and combine this with an analysis of the sequence coverage of different EV types in public databases (e.g., the Virus Pathogen Resource). This evaluation reveals temporal and geographic differences in EV circulation and sequence distribution, highlighting recent EV outbreaks and revealing gaps in sequence coverage. Phylogenetic analysis of the EV genus shows the relatedness of different EV types. Recombination analysis of the EV-A species provides evidence for recombination as a mechanism of genomic diversification. The absence of broadly protective vaccines and effective antivirals makes human enteroviruses important pathogens of public health concern.

Molecular typing and characterization of a novel genotype of EV-B93 isolated from Tibet, China

EV-B93 is a novel serotype within the Enterovirus B species and is uncommon worldwide. Currently, only one full-length genomic sequence (the prototype strain) has been deposited in the GenBank database. In this study, three EV-B93 were identified, including one from an acute flaccid paralysis (AFP) patient (named 99052/XZ/CHN/1999, hereafter XZ99052) and two from healthy children (named 99096/XZ/CHN/1999 and 99167/XZ/CHN/1999, hereafter XZ99096 and XZ99167, respectively) from Tibet in 1999 during the polio eradication program. The identity between the nucleotide and amino acid sequences of the Tibet EV-B93 strain and the EV-B93 prototype strain is 83.2%–83.4% and 96.8%–96.9%, respectively. The Tibet EV-B93 strain was found to have greater nucleotide sequence identity in the P3 region to another enterovirus EV-B107 as per a phylogenetic tree analysis, which revealed that recombination occurred. Seroepidemiology data showed that EV-B93 has not produced an epidemic in Tibet and there may be susceptible individuals. The three Tibet EV-B93 strains are temperature-resistant with prognosticative virulence, suggesting the possibility of a potential large-scale outbreak of EV-B93. The analyzed EV-B93 strains enrich our knowledge about this serotype and provide valuable information on global EV-B93 molecular epidemiology. What is more, they permit the appraisal of the serotype's potential public health impact and aid in understanding the role of recombination events in the evolution of enteroviruses.

Genetic recombination in fast-spreading coxsackievirus A6 variants: a potential role in evolution and pathogenicity

Hand, foot, and mouth disease (HFMD) is a common global epidemic. From 2008 onwards, many HFMD outbreaks caused by coxsackievirus A6 (CV-A6) have been reported worldwide. Since 2013, with a dramatically increasing number of CV-A6-related HFMD cases, CV-A6 has become the predominant HFMD pathogen in mainland China. Phylogenetic analysis based on the VP1 capsid gene revealed that subtype D3 dominated the CV-A6 outbreaks. Here, we performed a large-scale (near) full-length genetic analysis of global and Chinese CV-A6 variants, including 158 newly sequenced samples collected extensively in mainland China between 2010 and 2018. During the global transmission of subtype D3 of CV-A6, the noncapsid gene continued recombining, giving rise to a series of viable recombinant hybrids designated evolutionary lineages, and each lineage displayed internal consistency in both genetic and epidemiological features. The emergence of lineage –A since 2005 has triggered CV-A6 outbreaks worldwide, with a rate of evolution estimated at 4.17 × 10⁻³ substitutions site^-1 year⁻¹ based on a large number of monophyletic open reading frame sequences, and created a series of lineages chronologically through varied noncapsid recombination events. In mainland China, lineage –A has generated another two novel widespread lineages (–J and –L) through recombination within the enterovirus A gene pool, with robust estimates of occurrence time. Lineage –A, –J, and –L infections presented dissimilar clinical manifestations, indicating that the conservation of the CV-A6 capsid gene resulted in high transmissibility, but the lineage-specific noncapsid gene might influence pathogenicity. Potentially important amino acid substitutions were further predicted among CV-A6 variants. The evolutionary phenomenon of noncapsid polymorphism within the same subtype observed in CV-A6 was uncommon in other leading HFMD pathogens; such frequent recombination happened in fast-spreading CV-A6, indicating that the recovery of deleterious genomes may still be ongoing within CV-A6 quasispecies. CV-A6-related HFMD outbreaks have caused a significant public health burden and pose a great threat to children’s health; therefore, further surveillance is greatly needed to understand the full genetic diversity of CV-A6 in mainland China.

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.

Assessing the application of a pseudovirus system for emerging SARS-CoV-2 and re-emerging avian influenza virus H5 subtypes in vaccine development

Background

Highly pathogenic emerging and re-emerging viruses continuously threaten lives worldwide. In order to provide prophylactic prevention from the emerging and re-emerging viruses, vaccine is suggested as the most efficient way to prevent individuals from the threat of viral infection. Nonetheless, the highly pathogenic viruses need to be handled in a high level of biosafety containment, which hinders vaccine development. To shorten the timeframe of vaccine development, the pseudovirus system has been widely applied to examine vaccine efficacy or immunogenicity in the emerging and re-emerging viruses.

Methods

We developed pseudovirus systems for emerging SARS coronavirus 2 (SARS-CoV-2) and re-emerging avian influenza virus H5 subtypes which can be handled in the biosafety level 2 facility. Through the generated pseudovirus of SARS-CoV-2 and avian influenza virus H5 subtypes, we successfully established a neutralization assay to quantify the neutralizing activity of antisera against the viruses.

Results

The result of re-emerging avian influenza virus H5Nx pseudoviruses provided valuable information for antigenic evolution and immunogenicity analysis in vaccine candidate selection. Together, our study assessed the potency of pseudovirus systems in vaccine efficacy, antigenic analysis, and immunogenicity in the vaccine development of emerging and re-emerging viruses.

Conclusion

Instead of handling live highly pathogenic viruses in a high biosafety level facility, using pseudovirus systems would speed up the process of vaccine development to provide community protection against emerging and re-emerging viral diseases with high pathogenicity.

Genetic characterization of VP1 of coxsackieviruses A2, A4, and A10 associated with hand, foot, and mouth disease in Vietnam in 2012-2017: endemic circulation and emergence of new HFMD-causing lineages

While conducting sentinel surveillance of hand, foot, and mouth disease (HFMD) in Vietnam, we found a sudden increase in the prevalence of coxsackievirus A10 (CV-A10) in 2016 and CV-A2 and CV-A4 in 2017, the emergence of which has been reported recently to be associated with various clinical manifestations in other countries. However, there have been only a limited number of molecular studies on those serotypes, with none being conducted in Vietnam. Therefore, we sequenced the entire VP1 genes of CV-A10, CV-A4, and CV-A2 strains associated with HFMD in Vietnam between 2012 and 2017. Phylogenetic analysis revealed a trend of endemic circulation of Vietnamese CV-A10, CV-A4, and CV-A2 strains and the emergence of thus-far undescribed HFMD-causing lineages of CV-A4 and CV-A2. The Vietnamese CV-A10 strains belonged to a genotype comprising isolates from patients with HFMD from several other countries; however, most of the Vietnamese strains were grouped into a local lineage. Recently, emerging CV-A4 strains in Vietnam were grouped into a unique lineage within a genotype comprising strains isolated from patients with acute flaccid paralysis from various countries. New substitutions were detected in the putative BC and HI loops in the Vietnamese CV-A4 strains. Except for one strain, Vietnamese CV-A2 isolates were grouped into a unique lineage of a genotype that includes strains from various countries that are associated with other clinical manifestations. Enhanced surveillance is required to monitor their spread and to specify their roles as etiological agents of HFMD or "HFMD-like" diseases, especially for CV-A4 and CV-A2. Further studies including whole-genome sequencing should be conducted to fully understand the evolutionary changes occurring in these newly emerging strains.

Multirecombinant Enterovirus A71 Subgenogroup C1 Isolates Associated with Neurologic Disease, France, 2016-2017

In 2016, an upsurge of neurologic disease associated with infection with multirecombinant enterovirus A71 subgenogroup C1 lineage viruses was reported in France. These viruses emerged in the 2000s; 1 recombinant is widespread. This virus lineage has the potential to be associated with a long-term risk for severe disease among children.

Using Statistical Phylogenetics for Investigation of Enterovirus 71 Genotype A Reintroduction into Circulation

Neurovirulent enterovirus 71 (EV-A71) caused a massive epidemic in China in 2008-2011. While subgenotype C4 was the major causative agent, a few isolates were almost identical to the prototype EV-A71 strain and belonged to genotype A. This variant was allegedly extinct since 1970, and its identification in this epidemic suggests reintroduction of the archive virus. Regression analysis of genetic distances (TempEst software) was of moderate utility due to the low resolution of classical phylogenetic methods. Bayesian phylogenetic analysis (BEAST software) suggested artificial introduction event based on highly aberrant phylogenetic tree branch rates that differed by over three standard deviations from the mean substitution rate for EV71. Manual nucleotide-level analysis was used to further explore the virus spread pattern after introduction into circulation. Upon reintroduction, the virus accumulated up to seven substitutions in VP1, most of them non-synonymous and located within the capsid's canyon or at its rims, compatible with readaptation of a lab strain to natural circulation.

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.

Genomic analysis of serologically untypable human enteroviruses in Taiwan

BACKGROUND

Human enteroviruses contain over 100 serotypes. We have routinely conducted enterovirus surveillance in northern Taiwan; but about 10% of isolates could not be serotyped using traditional assays. Next-generation sequencing (NGS) is a powerful tool for genome sequencing.

METHODS

In this study, we established an NGS platform to conduct genome sequencing for the serologically untypable enterovirus isolates.

RESULTS

Among 130 serologically untypable isolates, 121 (93%) of them were classified into 29 serotypes using CODEHOP (COnsensus-DEgenerate Hybrid Oligonucleotide Primer)-based RT-PCR to amplify VP1 genes (VP1-CODEHOP). We further selected 52 samples for NGS and identified 59 genome sequences from 51 samples, including 8 samples containing two virus genomes. We also detected 23 genome variants (nucleotide identity < 90% compared with genome sequences in the public domain) which were potential genetic recombination, including 9 inter-serotype recombinants and 14 strains with unknown sources of recombination.

CONCLUSIONS

We successfully integrated VP1-CODEHOP and NGS techniques to conduct genomic analysis of serologically untypable enteroviruses.

Recombinant CV-A6 strains related to hand-foot-mouth disease and herpangina at primary care centers (Barcelona, Spain)

Aim: To describe the genetic diversity of enteroviruses (EV) causing hand, foot and mouth disease (HFMD) and herpangina, especially of coxsackievirus (CV)-A6, from patients attended at pediatric primary care centers during the 2017-2018 season. Methods: Phylogenetic analysis of partial VP1 region was performed for genetic characterization. The complete VP1 and 3Dpol proteins were sequenced for lineage determination and detection of recombination events. Results: An 80% of samples were EV laboratory-confirmed. CV-A6 was the most detected (70%) and associated with atypical HFMD (78%). The comparison of VP1 and 3Dpol phylogenies showed evidence of recombination in three strains, in which two shifted to CV-A16 3Dpol. Conclusion: The study provides recent information regarding the nonrecombinant and recombinant EVs related to HFMD at primary care centers.

Genetic landscape and macro-evolution of co-circulating Coxsackieviruses A and Vaccine-derived Polioviruses in the Democratic Republic of Congo, 2008-2013

Enteroviruses (EVs) are among the most common viruses infecting humans worldwide but only a few Non-Polio Enterovirus (NPEV) isolates have been characterized in the Democratic Republic of Congo (DR Congo). Moreover, circulating vaccine-derived polioviruses (PVs) [cVDPVs] isolated during multiple outbreaks in DR Congo from 2004 to 2018 have been characterized so far only by the sequences of their VP1 capsid coding gene. This study was carried to i) investigate the circulation and genetic diversity of NPEV and polio vaccine isolates recovered from healthy children and Acute Flaccid Paralysis (AFP) patients, ii) evaluate the occurrence of genetic recombination among EVs belonging to the Enterovirus C species (including PVs) and iii) identify the virological factors favoring multiple emergences of cVDPVs in DR Congo. The biological material considered in this study included i) a collection of 91 Sabin-like PVs, 54 cVDPVs and 150 NPEVs isolated from AFP patients between 2008 and 2012 in DR Congo and iii) a collection of 330 stool specimens collected from healthy children in 2013 in the Kasai Oriental and Maniema provinces of DR Congo. Studied virus isolates were sequenced in four distinct sub-genomic regions 5’-UTR, VP1, 2C^ATPase and 3D^pol. Resulting sequences were compared through comparative phylogenetic analyses. Virus isolation showed that 19.1% (63/330) healthy children were infected by EVs including 17.9% (59/330) of NPEVs and 1.2% (4/330) of type 3 Sabin-like PVs. Only one EV-C type, EV-C99 was identified among the NPEV collection from AFP patients whereas 27.5% of the 69 NPEV isolates typed in healthy children belonged to the EV-C species: CV-A13 (13/69), A20 (5/69) and A17 (1/69). Interestingly, 50 of the 54 cVDPVs featured recombinant genomes containing exogenous sequences in at least one of the targeted non-structural regions of their genomes: 5’UTR, 2C^ATPase and 3D^pol. Some of these non-vaccine sequences of the recombinant cVDPVs were strikingly related to homologous sequences from co-circulating CV-A17 and A20 in the 2C^ATPase region as well as to those from co-circulating CV-A13, A17 and A20 in the 3D^pol region. This study provided the first evidence uncovering CV-A20 strains as major recombination partners of PVs. High quality AFP surveillance, sensitive environmental surveillance and efficient vaccination activities remain essential to ensure timely detection and efficient response to recombinant cVDPVs outbreaks in DR Congo. Such needs are valid for any epidemiological setting where high frequency and genetic diversity of Coxsackieviruses A13, A17 and A20 provide a conducive viral ecosystem for the emergence of virulent recombinant cVDPVs.

The strategy of the Global Polio Eradication Initiative is based on the surveillance of patients suffering from Acute Flaccid Paralysis (AFP) and mass vaccination with live-attenuated vaccine strains of polioviruses (PVs) in endemic areas. However, vaccine strains of PVs can circulate and replicate for a long time when the vaccine coverage of the population is low. Such prolonged circulation and replication of vaccine strains of PVs can result to the emergence of circulating vaccine-derived polioviruses [cVDPVs] that are as virulent as wild PVs. In this study, we performed the molecular characterization of a large collection of 377 virus isolates recovered from paralyzed patients between 2008 and 2012 in DR Congo and healthy children in 2013 in the Kasai Oriental and Maniema provinces of DR Congo. We found that the genetic diversity of enteroviruses of the species Enterovirus C is more important than previously reported. Interestingly, 50 of the 54 cVDPVs featured recombinant genomes containing exogenous sequences of the 2C ATPase and/or 3D polymerase coding genes acquired from co-circulating Coxsackieviruses A13, A17 and A20. Coxsackieviruses A20 strains were identified for the first time as major partners of genetic recombination with co-circulating live-attenuated polio vaccine strains.

Our findings highlight the need to reinforce and maintain high quality surveillance of PVs and efficient immunization activities in order to ensure early detection and control of emerging cVDPVs in all settings where high frequency and diversity of Coxsackieviruses A13, A17 and A20 have been documented.

Hand, foot, and mouth disease associated with coxsackievirus A10: more serious than it seems

INTRODUCTION

Hand, foot, and mouth disease (HFMD) is a common viral childhood illness, that has been a severe public health concern worldwide, particularly in the Asia-Pacific region. According to epidemiological data of HFMD during the past decade, the most prevalent causal viruses were enterovirus (EV)-A71, coxsackievirus (CV)-A16, CV-A6, and CV-A10. The public health burden of CV-A10-related diseases has been underestimated as their incidence was lower than that of EV-A71 and CV-A16 in most HFMD outbreaks. However, cases of CV-A10 infection are more severe, and its genome is more variable, which has alerted the research community worldwide. Areas covered: In this paper, studies on the epidemiology, laboratory diagnosis, clinical manifestations, molecular epidemiology, seroepidemiology, animal models of CV-A10, and vaccines and antiviral strategies against this genotype are reviewed. In addition, the genetic evolution of circulating strains was analyzed. Expert opinion: Multivalent vaccines against EV-A71, CV-A16, CV-A6, and CV-A10 should be a next-step HFMD vaccine strategy.

High Permissiveness for Genetic Exchanges between Enteroviruses of Species A, including Enterovirus 71, Favors Evolution through Intertypic Recombination in Madagascar

Human enteroviruses of species A (EV-A) are the leading cause of hand-foot-and-mouth disease (HFMD). EV-A71 is frequently implicated in HFMD outbreaks and can also cause severe neurological manifestations. We investigated the molecular epidemiological processes at work and the contribution of genetic recombination to the evolutionary history of EV-A in Madagascar, focusing on the recently described EV-A71 genogroup F in particular. Twenty-three EV-A isolates, collected mostly in 2011 from healthy children living in various districts of Madagascar, were characterized by whole-genome sequencing. Eight different types were identified, highlighting the local circulation and diversity of EV-A. Comparative genome analysis revealed evidence of frequent recent intra- and intertypic genetic exchanges between the noncapsid sequences of Madagascan EV-A isolates. The three EV-A71 isolates had different evolutionary histories in terms of recombination, with one isolate displaying a mosaic genome resulting from recent genetic exchanges with Madagascan coxsackieviruses A7 and possibly A5 and A10 or common ancestors. The engineering and characterization of recombinants generated from progenitors belonging to different EV-A types or EV-A71 genogroups with distantly related nonstructural sequences indicated a high level of permissiveness for intertypic genetic exchange in EV-A. This permissiveness suggests that the primary viral functions associated with the nonstructural sequences have been highly conserved through the diversification and evolution of the EV-A species. No outbreak of disease due to EV-A has yet been reported in Madagascar, but the diversity, circulation, and evolution of these viruses justify surveillance of EV-A circulation and HFMD cases to prevent possible outbreaks due to emerging strains.IMPORTANCE Human enteroviruses of species A (EV-A), including EV-A71, are the leading cause of hand-foot-and-mouth disease (HFMD) and may also cause severe neurological manifestations. We investigated the circulation and molecular evolution of EV-A in Madagascar, focusing particularly on the recently described EV-A71 genogroup F. Eight different types, collected mostly in 2011, were identified, highlighting the local circulation and diversity of EV-A. Comparative genome analysis revealed evidence of frequent genetic exchanges between the different types of isolates. The three EV-A71 isolates had different evolutionary histories in terms of recombination. The engineering and characterization of recombinants involving progenitors belonging to different EV-A types indicated a high degree of permissiveness for genetic exchange in EV-A. No outbreak of disease due to EV-A has yet been reported in Madagascar, but the diversity, circulation, and evolution of these viruses justify the surveillance of EV-A circulation to prevent possible HFMD outbreaks due to emerging strains.

Molecular Epidemiology of Echovirus 18 Circulating in Mainland China from 2015 to 2016

Echovirus 18 (E18), a serotype of Enterovirus B (EV-B) species, is an important pathogen in aseptic meningitis. E18 had rarely been detected in mainland China, but became the predominant pathogen associated with viral encephalitis (VE) and meningitis in Hebei province for the first time in 2015. To investigate the molecular epidemiology and genetic characteristics of E18 in mainland China, sixteen E18 strains from patient throat swabs with hand, foot, and mouth disease (HFMD) in six provinces in China collected between 2015 and 2016, and four E18 strains isolated from 18 patient cerebrospinal fluid specimens with VE in Hebei Province in 2015 were obtained and sequenced. Combined with the sequences from the GenBank database, we performed an extensive genetic analysis. Phylogenetic analysis of VP1 gene sequences revealed that all E18 strains from mainland China after 2015 belonged to subgenotype C2. There were no obvious specific differences in phylogenetic and variation analyses of E18 genome sequences between HFMD and VE/meningitis strains. Potential multiple recombination may have occurred in the 5'-untranslated region and in the P2 and P3 nonstructural protein-encoding regions of E18 strains from China. The current E18 strains were potential multiple-recombinant viruses. Overall, these findings supported that E18 caused HFMD, VE, and meningitis, although there were no significant associations between clinical features and viral genomic characteristics.

Molecular surveillance of coxsackievirus A16 reveals the emergence of a new clade in mainland China

Coxsackievirus A16 (CV-A16) of the genotypes B1a and B1b have co-circulated in mainland China in the past decades. From 2013 to 2017, a total of 3,008 specimens from 3,008 patients with mild hand, foot, and mouth disease were collected in the present study. Viral RNA was tested for CV-A16 by a real-time RT-PCR method, and complete VP1 sequences and full-length genome sequences of CV-A16 strains from this study were determined by RT-PCR and sequencing. Sequences were analyzed using a series of bioinformatics programs. The detection rate for CV-A16 was 4.1%, 25.9%, 10.6%, 28.1% and 12.9% in 2013, 2014, 2015, 2016 and 2017, respectively. Overall, the detection rate for CV-A16 was 16.5% (497/3008) in this 5-year period in Shenzhen, China. One hundred forty-two (142/155, 91.6%) of the 155 genotype B1 strains in the study belonged to subgenotype B1b, and 13 (13/155, 8.4%) strains belonged to subgenotype B1a. Two strains (CVA16/Shenzhen174/CHN/2017 and CVA16/Shenzhen189/CHN/2017) could not be assigned to a known genotype. Phylogenetic analysis of these two strains and other Chinese CV-A16 strains indicated that these two CV-A16 strains clustered independently in a novel clade whose members differed by 8.4%-11.8%, 8.4%-12.1%, and 14.6%-14.8% in their nucleotide sequences from those of Chinese B1a, B1b, and genotype D strains, respectively. Phylogenetic analysis of global CV-A16 strains further indicated that the two novel CV-A16 strains from this study grouped in a previously uncharacterized clade, which was designated as the subgenogroup B3 in present study. Meanwhile, phylogenetic reconstruction revealed two other new genotypes, B1d and B4, which included a Malaysian strain and two American strains, respectively. The complete genome sequences of the two novel CV-A16 strains showed the highest nucleotide sequence identity of 92.3% to the Malaysian strain PM-15765-00 from 2000. Comparative analysis of amino acid sequences of the two novel CV-A16 strains and their relatives suggested that variations in the nonstructural proteins may play an important role in the evolution of modern CV-A16.

Molecular epidemiology and phylogenetics of human enteroviruses: Is there a forest behind the trees?

Enteroviruses are among the best studied small non-enveloped enteric RNA viruses. Most enteroviruses are easy to isolate in cell culture, and many non-polio enterovirus strains were archived worldwide as a byproduct of the WHO poliovirus surveillance system. Common outbreaks and epidemics, most prominently the epidemic of hand-foot-and-mouth disease with severe neurological complications in East and South-East Asia, justify practical interest of non-polio enteroviruses. As a result, there are over 50 000 enterovirus nucleotide sequences available in GenBank. Technical possibilities have been also improving, as Bayesian phylogenetic methods with an integrated molecular clock were introduced a decade ago and provided unprecedented opportunities for phylogenetic analysis. As a result, hundreds of papers were published on the molecular epidemiology of enteroviruses. This review covers the modern methodology, structure, and biases of the sequence dataset available in GenBank. The relevance of the subtype classification, findings of co-circulation of multiple genetic variants, previously unappreciated complexity of viral populations, and global evolutionary patterns are addressed. The most relevant conclusions and prospects for further studies on outbreak emergence mechanisms are discussed.

A complex mosaic of enteroviruses shapes community-acquired hand, foot and mouth disease transmission and evolution within a single hospital

Human enteroviruses (EV) pose a major risk to public health. This is especially so in the Asia-Pacific region where increasing numbers of hand, foot and mouth disease (HFMD) cases and large outbreaks of severe neurological disease associated with EV-A71 have occurred. Despite their importance, key aspects of the emergence, epidemiology and evolution of EVs remain unclear, and most studies of EV evolution have focused on a limited number of genes. Here, we describe the genomic-scale evolution of EV-A viruses sampled from pediatric patients with mild disease attending a single hospital in western Sydney, Australia, over an 18-month period. This analysis revealed the presence of eight viral serotypes-Coxsackievirus (CV) A2, A4, A5, A6, A8, A10, A16 and EV-A71-with up to four different serotypes circulating in any 1 month. Despite an absence of large-scale outbreaks, high levels of geographical and temporal mixing of serotypes were identified. Phylogenetic analysis revealed that multiple strains of the same serotype were present in the community, and that this diversity was shaped by multiple introductions into the Sydney population, with only a single lineage of CV-A6 exhibiting in situ transmission over the entire study period. Genomic-scale analyses also revealed the presence of novel and historical EV recombinants. Notably, our analysis revealed no association between viral phylogeny, including serotype, and patient age, sex, nor disease severity (for uncomplicated disease). This study emphasizes the contribution of EV-A viruses other than EV-A71 to mild EV disease including HFMD in Australia and highlights the need for greater surveillance of these viruses to improve strategies for outbreak preparedness and vaccine design.

Antigenic characteristics and genomic analysis of novel EV-A90 enteroviruses isolated in Xinjiang, China

Enterovirus A90 (EV-A90) is a novel serotype of enterovirus A species that is rarely reported. Here, we isolated five enteroviruses from patients with acute flaccid paralysis in Hotan and Kashgar cities in Xinjiang, China that were identified as EV-A90 by molecular typing. The VP1 sequences of these Xinjiang EV-A90 strains showed 88.4–89% nucleotide sequence identity to the prototype EV-A90 strain; however, genome analysis indicated complex recombination events in P2 and P3 regions. Next, the seroprevalence of EV-A90 was examined in 49 serum specimens collected in Hotan and Kashgar, and 37.5% were EV-A90 antibody positive (>1:8), with a geometric mean titre (GMT) of 1:10.47. The low positive rate and GMT suggest a low-level EV-A90 epidemic in Xinjiang. Two of the five Xinjiang EV-A90 strains were temperature sensitive, and three were temperature resistant, and a comparative genomics analysis suggested that an amino acid substitution (H1799Y) in the 3D^pol region was related to temperature sensitivity. Although the epidemic strength is low, some EV-A90 strains were temperature resistant, which is suggestive of strong virulence and transmission capacity. This study expanded the number of EV-A90 in GenBank and provided basic data that may be useful for studying the molecular epidemiology of EV-A90.

Recombination in the rabies virus and other lyssaviruses

Recombination is a common event in RNA viruses; however, in the rabies virus there have been only a few reports of isolated recombination events. Comprehensive analysis found traces of recent recombination events within Arctic, Arctic-like and Africa 1b rabies virus groups, as well as recombination between distinct lyssaviruses. Recombination breakpoints were not linked to gene boundaries and could be detected all over the genome. However, there was no evidence that recombination is an important factor in the genetic variability of the rabies virus. It is therefore likely that recombination in the rabies virus is limited by ecological factors (e.g., rare co-circulation of distinguishable lineages and a narrow window for productive coinfection in most carnivore hosts), rather than molecular barriers (e.g., incompatibility of genome fragments).

Genetic diversity and recombination of enterovirus G strains in Japanese pigs: High prevalence of strains carrying a papain-like cysteine protease sequence in the enterovirus G population

To study the genetic diversity of enterovirus G (EV-G) among Japanese pigs, metagenomics sequencing was performed on fecal samples from pigs with or without diarrhea, collected between 2014 and 2016. Fifty-nine EV-G sequences, which were >5,000 nucleotides long, were obtained. By complete VP1 sequence analysis, Japanese EV-G isolates were classified into G1 (17 strains), G2 (four strains), G3 (22 strains), G4 (two strains), G6 (two strains), G9 (six strains), G10 (five strains), and a new genotype (one strain). Remarkably, 16 G1 and one G2 strain identified in diarrheic (23.5%; four strains) or normal (76.5%; 13 strains) fecal samples possessed a papain-like cysteine protease (PL-CP) sequence, which was recently found in the USA and Belgium in the EV-G genome, at the 2C–3A junction site. This paper presents the first report of the high prevalence of viruses carrying PL-CP in the EV-G population. Furthermore, possible inter- and intragenotype recombination events were found among EV-G strains, including G1-PL-CP strains. Our findings may advance the understanding of the molecular epidemiology and genetic evolution of EV-Gs.

Molecular evolution of types in non-polio enteroviruses

Non-polio enteroviruses are a ubiquitous and divergent group of non-enveloped RNA viruses. Novel types are reported regularly in addition to over 100 known types; however, mechanisms of emergence of novel types remain obscure. Here, the 33 most common types represented by 35-629 non-redundant partial VP1 sequences in GenBank were studied in parallel using Bayesian coalescent molecular clock analysis to investigate common evolutionary trends among enterovirus types. Inferred substitution rates were in the range of 0.41×10^-2 to 3.07×10^-2 substitutions per site per year. The most recent common ancestors of known isolates of each type presumably existed between 55 and 200 years ago. Phylogenetic analysis results suggested that global type populations underwent bottlenecks that could repeatedly reset the common ancestor dates. Nevertheless, species-level analysis suggested that the contemporary enterovirus types emerged within the last millennium. Analysis of 2657 complete VP1 sequences of the 24 most common types indicated that the type criterion based upon 75 % nucleotide sequence identity remains generally valid, despite exponential growth of the number of known sequences and a high rate of mutation fixation. However, in few types there was evidence that enteroviruses can drift slightly beyond the type threshold, up to 73 % identity, and both amino acid and nucleotide sequences should be considered for type identification. Analysis of sequence distances within types implied that sequence-identity-based identification of genotypes is rational within some, but not all, types and distinct genotype cut-offs (9-20 %) may be useful for different types.

Independent evolution of tetraloop in enterovirus oriL replicative element and its putative binding partners in virus protein 3C

Background

Enteroviruses are small non-enveloped viruses with a (+) ssRNA genome with one open reading frame. Enterovirus protein 3C (or 3CD for some species) binds the replicative element oriL to initiate replication. The replication of enteroviruses features a low-fidelity process, which allows the virus to adapt to the changing environment on the one hand, and requires additional mechanisms to maintain the genome stability on the other. Structural disturbances in the apical region of oriL domain d can be compensated by amino acid substitutions in positions 154 or 156 of 3C (amino acid numeration corresponds to poliovirus 3C), thus suggesting the co-evolution of these interacting sequences in nature. The aim of this work was to understand co-evolution patterns of two interacting replication machinery elements in enteroviruses, the apical region of oriL domain d and its putative binding partners in the 3C protein.

Methods

To evaluate the variability of the domain d loop sequence we retrieved all available full enterovirus sequences (>6, 400 nucleotides), which were present in the NCBI database on February 2017 and analysed the variety and abundance of sequences in domain d of the replicative element oriL and in the protein 3C.

Results

A total of 2,842 full genome sequences was analysed. The majority of domain d apical loops were tetraloops, which belonged to consensus YNHG (Y = U/C, N = any nucleotide, H = A/C/U). The putative RNA-binding tripeptide 154–156 (Enterovirus C 3C protein numeration) was less diverse than the apical domain d loop region and, in contrast to it, was species-specific.

Discussion

Despite the suggestion that the RNA-binding tripeptide interacts with the apical region of domain d, they evolve independently in nature. Together, our data indicate the plastic evolution of both interplayers of 3C-oriL recognition.

Molecular characterization of two novel echovirus 18 recombinants associated with hand-foot-mouth disease

Human echovirus 18 (E-18) is a member of the enterovirus B species. To date, sixteen full-length genome sequences of E-18 are available in the GenBank database. In this study, we describe the complete genomic characterization of two E-18 strains isolated in Yunnan, China. Pairwise comparisons of the nucleotide sequences and the deduced amino acid sequences revealed that the two Yunnan E-18 strains had 87.5% nucleotide identity and 96.3–96.5% amino acid identity with the Chinese strain. Phylogenetic and bootscanning analyses revealed the two E-18 strains had the highest identity with other several EV-B serotypes than the other E-18 strains in the P3 coding region, especially, 3B region of the Swine Vesicular disease virus (SVDV) strain HK70, indicated that frequent intertypic recombination might have occurred in the two Yunnan strains. This study contributes the complete genome sequences of E-18 to the GenBank database and provides valuable information on the molecular epidemiology of E-18 in China.