Recombination among human non-polio enteroviruses: implications for epidemiology and evolution

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.

Rapid detection of four major HFMD-associated enteroviruses by multiplex HiFi-LAMP assays

Hand, foot, and mouth disease (HFMD) caused by various enteroviruses is a major public health concern globally. Human enterovirus 71(EVA71), coxsackievirus A16 (CVA16), coxsackievirus A6 (CVA6), and coxsackievirus A10 (CVA10) are four major enteroviruses responsible for HFMD. Rapid, accurate, and specific point-of-care (POC) detection of the four enteroviruses is crucial for the prevention and control of HFMD. Here, we developed two multiplex high-fidelity DNA polymerase loop-mediated isothermal amplification (mHiFi-LAMP) assays for simultaneous detection of EVA71, CVA16, CVA6, and CVA10. The assays have good specificity and exhibit high sensitivity, with limits of detection (LOD) of 11.2, 49.6, 11.4, and 20.5 copies per 25 μL reaction for EVA71, CVA16, CVA6, and CVA10, respectively. The mHiFi-LAMP assays showed an excellent clinical performance (sensitivity 100.0%, specificity 83.3%, n = 47) when compared with four singleplex RT-qPCR assays (sensitivity 93.1%, specificity 100%). In particular, the HiFi-LAMP assays exhibited better performance (sensitivity 100.0%, specificity 100%) for CVA16 and CVA6 than the RT-qPCR assays (sensitivity 75.0-92.3%, specificity 100%). Furthermore, the mHiFi-LAMP assays detected all clinical samples positive for the four enteroviruses within 30 min, obviously shorter than about 1-1.5 h by the RT-qPCR assays. The new mHiFi-LAMP assays can be used as a robust point-of-care testing (POCT) tool to facilitate surveillance of HFMD at rural and remote communities and resource-limited settings.

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.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.

Phylogenetic characteristics and recombination analysis of echovirus 5 associated with severe acute respiratory infection in China

IMPORTANCE

This study is the first report of echovirus 5 (E5) associated with severe acute respiratory infection and obtained the first E5 whole-genome sequence in China. Combined with the sequences available in the GenBank database, the first genotyping, phylogenetic characteristics, recombination, and genetic evolutionary analysis of E5 was performed in this study. Our findings providing valuable information on global E5 molecular epidemiology.

Genetic Characteristics and Phylogeographic Dynamics of Echovirus

Echoviruses belong to the genus Enterovirus in the Picornaviridae family, forming a large group of Enterovirus B (EV-B) within the Enteroviruses. Previously, Echoviruses were classified based on the coding sequence of VP1. In this study, we performed a reliable phylogenetic classification of 277 sequences isolated from 1992 to 2019 based on the full-length genomes of Echovirus. In this report, phylogenetic, phylogeographic, recombination, and amino acid variability landscape analyses were performed to reveal the evolutional characteristics of Echovirus worldwide. Echoviruses were clustered into nine major clades, e.g., G1-G9. Phylogeographic analysis showed that branches G2-G9 were linked to common strains, while the branch G1 was only linked to G5. In contrast, strains E12, E14, and E16 clustered separately from their G3 and G7 clades respectively, and became a separate branch. In addition, we identified a total of 93 recombination events, where most of the events occurred within the VP1-VP4 coding regions. Analysis of amino acid variation showed high variability in the a positions of VP2, VP1, and VP3. This study updates the phylogenetic and phylogeographic information of Echovirus and indicates that extensive recombination and significant amino acid variation in the capsid proteins drove the emergence of new strains.

Enteroviral 2B Interacts with VDAC3 to Regulate Reactive Oxygen Species Generation That Is Essential to Viral Replication

Enterovirus (EV) 71 caused episodes of outbreaks in China and Southeast Asia during the last few decades. We have previously reported that EV71 induces reactive oxygen species (ROS). However, the underlying mechanism remains elusive. Co-immunoprecipitation-proteomic analysis revealed that enteroviral 2B protein interacted with mitochondrial voltage-dependent anion channel 3 (VDAC3). Knockdown (KD) of VDAC3 expression specifically inhibited enteroviral replication. Single-round viral replication was also inhibited in KD cells, suggesting that VDAC3 plays an essential role in replication. Consistent with this, VDAC3 gene KD significantly reduced the EV71-induced mitochondrial ROS generation. Exogenous 2B expression could induce the mitochondrial ROS generation that was significantly reduced in VDAC3-KD cells or in the Mito-TEMPO-treated cells. Moreover, VDAC3 appears to be necessary for regulation of antioxidant metabolism. VDAC3 gene KD led to the enhancement of such pathways as hypotaurine/taurine synthesis in the infected cells. Taken together, these findings suggest that 2B and VDAC3 interact to enhance mitochondrial ROS generation, which promotes viral replication.

Human pathogenic RNA viruses establish noncompeting lineages by occupying independent niches

Numerous pathogenic viruses are endemic in humans and cause a broad variety of diseases, but what is their potential for causing new pandemics? We show that most human pathogenic RNA viruses form multiple, cocirculating lineages with low turnover rates. These lineages appear to be largely noncompeting and occupy distinct epidemiological niches that are not regionally or seasonally defined, and their persistence appears to stem from limited outbreaks in small communities so that only a small fraction of the global susceptible population is infected at any time. However, due to globalization, interaction and competition between lineages might increase, potentially leading to increased diversification and pathogenicity. Thus, endemic viruses appear to merit global attention with respect to the prevention of future pandemics.

Many pathogenic viruses are endemic among human populations and can cause a broad variety of diseases, some potentially leading to devastating pandemics. How virus populations maintain diversity and what selective pressures drive population turnover is not thoroughly understood. We conducted a large-scale phylodynamic analysis of 27 human pathogenic RNA viruses spanning diverse life history traits, in search of unifying trends that shape virus evolution. For most virus species, we identify multiple, cocirculating lineages with low turnover rates. These lineages appear to be largely noncompeting and likely occupy semiindependent epidemiological niches that are not regionally or seasonally defined. Typically, intralineage mutational signatures are similar to interlineage signatures. The principal exception are members of the family Picornaviridae, for which mutations in capsid protein genes are primarily lineage defining. Interlineage turnover is slower than expected under a neutral model, whereas intralineage turnover is faster than the neutral expectation, further supporting the existence of independent niches. The persistence of virus lineages appears to stem from limited outbreaks within small communities, so that only a small fraction of the global susceptible population is infected at any time. As disparate communities become increasingly connected through globalization, interaction and competition between lineages might increase as well, which could result in changing selective pressures and increased diversification and/or pathogenicity. Thus, in addition to zoonotic events, ongoing surveillance of familiar, endemic viruses appears to merit global attention with respect to the prevention or mitigation of future pandemics.

Detection of recombinant breakpoint in the genome of human enterovirus E11 strain associated with a fatal nosocomial outbreak

BACKGROUND

The aim of this study was to characterize the genome of a recombinant Enterovirus associated with severe and fatal nosocomial infection; it was typed as Echovirus 11 (E-11) according to the VP1 gene. Enterovirus infection is generally asymptomatic and self-limited, but occasionally it may progress to a more severe clinical manifestation, as in the case described here. Recombination plays a crucial role in the evolution of Enteroviruses (EVs) and has been recognized as the main driving force behind the emergence of epidemic strains associated with severe infection. Therefore, it is of utmost importance to monitor the circulation of recombinant strains for surveillance purposes.

METHODS

Enterovirus-RNA was detected in the serum and liver biopsy of patients involved in the nosocomial cluster by commercial One-Step qRT-PCR method and the Enterovirus strains were isolated in vitro. The EVs typing was determined by analyzing the partial-length of the 5'UTR and VP1 sequences with the web-based open-access Enterovirus Genotyping Tool Version 0.1. The amplicons targeting 5'UTR, VP1 and overlapping fragments of the entire genome were sequenced with the Sanger method. Phylogenetic analysis was performed comparing the VP1 and the full-genome sequences of our strains against an appropriate reference set of Enterovirus prototypes of the Picornaviridae genera and species retrieved from the Enterovirus Genotyping Tool. Recombination analysis was performed using RDP4 software.

RESULTS

The Neighbor-Joining tree of the VP1 gene revealed that the 4 patients were infected with an identical molecular variant of Echovirus 11 (E-11). While the phylogenetic and the RDP4 analysis of the full-genome sequences provided evidence that it was a chimeric strain between an E-11 and a Coxsackievirus B (CV-B).

CONCLUSIONS

The chimeric structure of the E-11 genome might have contributed to the severe infection and epidemic feature of the strain, but further biological characterizations are needed. The evidence reported in this study, highlights the limit of typing techniques based on the VP1 gene, as they fail to identify the emergence of recombinant strains with potentially more pathogenic or epidemic properties, thus providing only partial information on the epidemiology and pathogenesis of Enteroviruses.

Origin and evolution analysis and genetic characteristics of echovirus 9 in China

Background

Echovirus 9 (E9) is associated with a wide variety of diseases and medical conditions, and the clinical symptoms of sporadic cases caused by E9 often are severe. With a high global prevalence, E9 has caused multiple outbreaks worldwide. However, little is known about the genetic and geographic population dynamics of E9.

Method

A total of 131 VP1 gene sequences, including15 generated in this study and 116 obtained from GenBank, were used to coestimate time-resolved phylogenies to infer viral evolution and transmission in worldwide. Overlapping fragments representing whole genomes were amplified by reverse transcription polymerase chain reaction (RT-PCR) using specific primers. Then, we reported the genetic characteristics of fifteen E9 strains in the Chinese Mainland. Similarity plots and bootscanning analysis were used to determine recombination patterns of E9.

Results

The estimated mean evolutionary rate of global E9 VP1 gene was 4.278 × 10⁻³ substitutions per site per year (95% confidence interval [CI], 3.822 × 10⁻³/site/year to 4.710 × 10⁻³/site/year), and the common ancestor of E9 likely emerged around 1868 (95% CI, 1840 to 1892). The full-length genomic sequences of the fifteen E9 strains showed 76.9–79.6% nucleotide identity and 95.3–95.9% amino acid identity with E9 Barty strain. 11 of 15 E9 whole genome sequence present four recombination patterns, and E9 recombinants have extensive genetic exchanges in the 2C and P3 regions with other Enterovirus B (EV-B) circulated in China. Four of six E9 strains were temperature sensitive, and two were temperature resistant, and a comparative genomics analysis suggested that 411, 865 and 867 amino acid substitution in the P1 region was related to temperature sensitivity.

Conclusion

This study highlights a persistent transmission network of E9 in worldwide, provides valuable information regarding the molecular epidemiology of E9.

Whole-genome analysis of coxsackievirus B3 reflects its genetic diversity in China and worldwide

Background

Coxsackievirus B3 (CVB3) has emerged as an active pathogen in myocarditis, aseptic meningitis, hand, foot, and mouth disease (HFMD), and pancreatitis, and is a heavy burden on public health. However, CVB3 has not been systematically analyzed with regard to whole-genome diversity and recombination. Therefore, this study was undertaken to systematically examine the genetic characteristics of CVB3 based on its whole genome.

Methods

We combined CVB3 isolates from our national HFMD surveillance and global sequences retrieved from GenBank. Phylogenetic analysis was performed to examine the whole genome variety and recombination forms of CVB3 in China and worldwide.

Results

Phylogenetic analysis showed that CVB3 strains isolated worldwide could be classified into clusters A–E based on the sequence of the entire VP1 region. The predominant CVB3 strains in China belonged to cluster D, whereas cluster E CVB3 might be circulated globally compared to other clusters. The average nucleotide substitution rate in the P1 region of CVB3 was 4.82 × 10–3 substitutions/site/year. Myocarditis was more common with cluster A. Clusters C and D presented more cases of acute flaccid paralysis, and cluster D may be more likely to cause HFMD. Multiple recombination events were detected among CVB3 variants, and there were twenty-three recombinant lineages of CVB3 circulating worldwide.

Conclusions

Overall, this study provides full-length genomic sequences of CVB3 isolates with a wide geographic distribution over a long-term time scale in China, which will be helpful for understanding the evolution of this pathogen. Simultaneously, continuous surveillance of CVB3 is indispensable to determine its genetic diversity in China as well as worldwide.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

Enterovirus D: A Small but Versatile Species

Enteroviruses (EVs) from the D species are the causative agents of a diverse range of infectious diseases in spite of comprising only five known members. This small clade has a diverse host range and tissue tropism. It contains types infecting non-human primates and/or humans, and for the latter, they preferentially infect the eye, respiratory tract, gastrointestinal tract, and nervous system. Although several Enterovirus D members, in particular EV-D68, have been associated with neurological complications, including acute myelitis, there is currently no effective treatment or vaccine against any of them. This review highlights the peculiarities of this viral species, focusing on genome organization, functional elements, receptor usage, and pathogenesis.

Molecular characterization of non-polio enteroviruses isolated from acute flaccid paralysis patients in Uganda

Enteroviruses (EVs) are RNA viruses that can cause many clinical syndromes including acute flaccid paralysis (AFP). Within the global polio laboratory network, EVs are categorized either as polioviruses or non-polio enteroviruses (NPEVs). Specific NPEVs have been described in polio-like residual paralytic events in AFP patients. Retrospective analysis of 112 NPEV isolates from AFP patients was performed and thirty one NPEV types were identified of which 91% were Enterovirus B and 9% were Enterovirus A species. The NPEVs were distributed across the country with most patients in the eastern region (41/89; 46.1%). The highest proportion of patients were children less than 5 years (77/89; 86.5%) and male patients were more common (54/89; 60.7%). Echovirus 11 (11/89; 12.4%) was frequently observed and phylogenetic analysis of these sequences revealed high diversity. Coxsackievirus B5 (CV-B5), CV-B6, E21, and EV-B69 were only seen in patients with residual paralysis. Analyses of the EV-A71 sequence indicated a unique genogroup.

Temporal phylogeny and molecular characterization of echovirus 30 associated with aseptic meningitis outbreaks in China

BACKGROUND

An outbreak of aseptic meningitis occurred from June to August 2016, in Inner Mongolia Autonomous Region, China.

METHODS

To determine its epidemiological characteristics, etiologic agent, and possible origin, specimens were collected for virus isolation and identification, followed by molecular epidemiological analysis.

RESULTS

A total of 363 patients were clinically diagnosed from June 1st to August 31st 2016, and most cases (63.1%, n = 229) were identified between June 22nd and July 17th, with children aged 6 to 12 years constituting the highest percentage (68.9%, n = 250). All viral isolates from this study belonged to genotype C of echovirus 30 (E30), which dominated transmission in China. To date, two E30 transmission lineages have been identified in China, of which Lineage 2 was predominant. We observed fluctuant progress of E30 genetic diversity, with Lineage 2 contributing to increased genetic diversity after 2002, whereas Lineage 1 was significant for the genetic diversity of E30 before 2002.

CONCLUSIONS

We identified the epidemiological and etiological causes of an aseptic meningitis outbreak in Inner Mongolia in 2016, and found that Lineage 2 played an important role in recent outbreaks. Moreover, we found that Gansu province could play an important role in E30 spread and might be a possible origin site. Furthermore, Fujian, Shandong, Taiwan, and Zhejiang provinces also demonstrated significant involvement in E30 evolution and persistence over time in China.

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.

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.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

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.

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.

Reinforced poliovirus and enterovirus surveillance in Romania, 2015-2016

Due to the risk of poliovirus importation from Ukraine in 2015, a combined surveillance program monitoring the circulation of enteroviruses (EVs) in healthy children from at-risk areas and in the environment was conducted in Romania. Virological testing of stool samples collected from 155 healthy children aged from two months to six years and of 186 sewage water samples collected from different areas was performed. A total of 58 (37.42%) stool samples and 50 (26.88%) sewage water samples were positive for non-polio EVs, but no poliovirus was detected. A high level of circulation of echovirus (E) types 6 and 7 and coxsackievirus (CV) type B5 was observed.

Neonatal hand, foot, and mouth disease due to coxsackievirus A6 in Shanghai

BACKGROUND

Evidence of hand, foot, and mouth disease (HFMD) in neonates is limited. The aim of this study was to evaluate the clinical symptoms, pathogens, possible transmission routes, and prognosis of neonatal HFMD in Shanghai.

METHODS

This was a case-control study based on the HFMD registry surveillance system. All neonates and infected family members were enrolled between 2016 and 2017 in Shanghai. Neonates with HFMD were followed for at least half a year. Detailed questionnaires, medical history, and physical examination were recorded. Routine blood examination, liver and renal function, immunophenotypes of peripheral blood lymphocytes (CD3, CD4, and CD8 T-cells; NK cells), immunoglobulin (Ig) M, IgG, and IgA, and cytokine interleukin (IL-1β, IL-2R, IL-6, IL-8, IL-10, and TNF-α) levels were measured. All rectal swab specimens were collected and genotyped for enterovirus, and phylogenetic analysis based on the VP1 sequences of coxsackievirus A6 (CV-A6) was performed to investigate molecular and evolutionary characteristics. T-test or nonparametric test was used to evaluate the differences. Logistic analysis was applied to calculate the risk of clinical manifestations in the group of HFMD neonates and their paired siblings.

RESULTS

There were 16 neonates among the 12,608 diagnosed patients with HFMD, accounting for 0.13%. All neonatal infections were transmitted by other members of the family, mainly the elder siblings, and were caused by CV-A6. CV-A6 was the emerging and predominant causative agent of HFMD in Shanghai. None of the neonates with HFMD experienced fever, onychomadesis, or severe complications. However, two elder sibling patients showed lethargy, and one developed hypoperfusion. In the elder siblings with HFMD, the proportion of white blood cells was generally higher than in neonates with HFMD. The immunologic function of the neonates with HFMD was basically normal. The levels of inflammatory markers were higher in both neonates and elder siblings with HFMD compared to age-matched controls. The clinical symptoms receded about 1 week after onset. None of the neonates had sequelae.

CONCLUSIONS

In our study, CV-A6 infection in neonates was benign, but had the character of family clustering. Due to the two-child policy in China, elder siblings may be the main route of HFMD transmission.

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.

Investigating the etiologic agents of aseptic meningitis outbreak in Iranian children

Introduction:

This study aimed to determine the viral agent (s) and their genome burden associated with an aseptic meningitis outbreak that occurred in Yasuj, Iran.

Materials and Methods:

During April to August 2015, 104 CSF samples from 104 patients under 14 years old admitted to the hospital of Yasuj, Iran, with aseptic meningitis associated clinical signs were collected. 200 μl CSF specimens was prepared for DNA and RNA viral genome extraction each and then subjected to diagnostic Taq-man real time PCR assays for the present of Enteroviruses, HSV, VZV, mumps, measles and rubella in the samples.

Results:

The majority of them had experienced clinical meningitis sign. Primary laboratory differentiated tests were in favor of viral meningitis. Among a total of 104 patients diagnosed with clinically aseptic meningitis, enterovirus as the most significant viral agent was detected in 53 subjects. However, mumps, HSV and VZV, as the endemic causes of viral meningitis, were detected in 6, 6 and 2 of the affected patients. It was revealed that two HSV and one VZV affected patients were co-infected with enteroviruses. All affected children with relatively variable viral load recovered without any sequels.

Conclusion:

The present study revealed enterovirus as the main predominant cause of pediatric aseptic meningitis that broke out in Yasuj-Iran. Also, the co-circulation of mumps, HSV and VZV, as the endemic cause during the same aseptic meningitis outbreak, was demonstrated in some cases.

Dissemination of Internal Ribosomal Entry Sites (IRES) Between Viruses by Horizontal Gene Transfer

Members of Picornaviridae and of the Hepacivirus, Pegivirus and Pestivirus genera of Flaviviridae all contain an internal ribosomal entry site (IRES) in the 5'-untranslated region (5'UTR) of their genomes. Each class of IRES has a conserved structure and promotes 5'-end-independent initiation of translation by a different mechanism. Picornavirus 5'UTRs, including the IRES, evolve independently of other parts of the genome and can move between genomes, most commonly by intratypic recombination. We review accumulating evidence that IRESs are genetic entities that can also move between members of different genera and even between families. Type IV IRESs, first identified in the Hepacivirus genus, have subsequently been identified in over 25 genera of Picornaviridae, juxtaposed against diverse coding sequences. In several genera, members have either type IV IRES or an IRES of type I, II or III. Similarly, in the genus Pegivirus, members contain either a type IV IRES or an unrelated type; both classes of IRES also occur in members of the genus Hepacivirus. IRESs utilize different mechanisms, have different factor requirements and contain determinants of viral growth, pathogenesis and cell type specificity. Their dissemination between viruses by horizontal gene transfer has unexpectedly emerged as an important facet of viral evolution.

Mosaic genome of Human Coxsackievirus A4 associated with herpangina and HFMD in Yancheng, China, 2016 and 2018

OBJECTIVES

To better understand the spectrums of pathogens causing herpangina and circulation of Coxsackievirus A4 in Yancheng, China.

METHODS

Stool samples from herpangina and HFMD cases were collected. Real Time PCR Kits was used to identify Enterovirus 71, CV-A16 and CV-A6, and nested reverse transcription PCR (nRT-PCR) to detect the other enterovirus types. Complete VP1 and genome sequence of CV-A4 were amplified by using nRT-PCR. Genetic, phylogenetic and recombination analysis were performed.

RESULTS

Co-circulation of three recombinant CV-A4 groups, including one novel (C2 lineage), was identified in Yancheng, China, 2016 and 2018. One was the major causative agent of herpangina, and another two were responsible for HFMD. Phylogenetic and recombination analysis indicated that the non-structural region of their genome originated from the same ancestry and subsequently adaptation. C2 lineage of CV-A4 group may be introduced from countries outside China and its genome occurred recombination in China.

CONCLUSION

Novel recombinant CV-A4 was mainly associated with herpanginain in Yancheng, 2018, China. C2 lineage of CV-A4 group with recombinant non-structural region was also identified in HFMD patients.

Parechovirus A prevalence in adults in The Netherlands

Human parechoviruses (HPeV) of the species Parechovirus A are highly prevalent disease-causing pathogens in children worldwide. HPeVs are capable of causing severe disease in adults as well, but the prevalence in adults may be much lower. The aim of our present study was to determine the prevalence of HPeV in clinical samples from adults sent in for diagnostic procedures in a tertiary hospital in the Netherlands. From a total of 10,645 samples obtained from 6175 patients, 20 samples from 11 patients (0.18%) tested positive for HPeV by RT-PCR. Two patients were positive for HPeV-1, two for HPeV-3, and one for HPeV-6. Six HPeVs could not be typed. Eight of the 11 HPeV-positive patients were immunocompromised. Due to comorbidity, we were unable to attribute the patients’ clinical symptoms to the HPeV infection. The HPeV prevalence in adults found in this study is low compared to HPeV prevalence in children. This may be largely explained by the high seropositivity rates in adults, although there could be other mechanisms involved.

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.

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.

Discovery of a virus of the species Enterovirus F in goats

Here, we report two novel enteroviruses, designated as SD-S67 and SD-S68, isolated from a goat farm. Their complete genome sequences were determined and found to be 7455 and 7465 nucleotides in length, respectively. Molecular characterization revealed that SD-S67 is closely related to bovine enterovirus strain 261 and that SD-S68 to caprine enterovirus strain CEV-JL14. Phylogenetic analysis showed that SD-S67 clustered with members of the species Enterovirus F, and that SD-S68 clustered with enteroviruses of goats and sheep. Recombination analysis showed that SD-S67 is likely to have undergone several recombination events in the process of its evolution. To the best of our knowledge, this is the first report of an enterovirus F isolate from a goat and of a coinfection with enteroviruses of different species in the same goat herd.

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.

Electronic supplementary material

The online version of this article (10.1186/s12929-019-0541-x) contains supplementary material, which is available to authorized users.

Current Understanding of Human Enterovirus D68

Human enterovirus D68 (EV-D68), a member of the species Enterovirus D of the Picornaviridae family, was first isolated in 1962 in the United States. EV-D68 infection was only infrequently reported until an outbreak occurred in 2014 in the US; since then, it has continued to increase worldwide. EV-D68 infection leads to severe respiratory illness and has recently been reported to be linked to the development of the neurogenic disease known as acute flaccid myelitis (AFM), mostly in children, seriously endangering public health. Hitherto, treatment options for EV-D68 infections were limited to supportive care, and as yet there are no approved, specific antiviral drugs or vaccines. Research on EV-D68 has mainly focused on its epidemiology, and its virologic characteristics and pathogenesis still need to be further explored. Here, we provide an overview of current research on EV-D68, including the genotypes and genetic characteristics of recent epidemics, the mechanism of infection and virus-host interactions, and its relationship to acute flaccid myelitis (AFM), in order to broaden our understanding of the biological features of EV-D68 and provide a basis for the development of effective antiviral agents.

Phylogenetic analysis and phenotypic characterisatics of two Tibet EV-C96 strains

BACKGROUND

Enterovirus C96 (EV-C96) is a newly named type of enterovirus belonging to species C, and the prototype strain (BAN00-10488) was firstly isolated in 2000 from a stool specimen of a patient with acute flaccid paralysis in Bangladesh. In this study, we report the genomic and phenotypic characteristics of two EV-C96 strains isolated from individuals from the Tibet Autonomous Region of China.

METHODS

Human rhabdomyosarcoma (RD), human laryngeal epidermoid carcinoma (HEp-2), and human cervical cancer (Hela) cells were infected with the Tibet EV-C96 strains, and enterovirus RNA in the cell culture was detected with a real time RT-PCR-based enterovirus screening method. The temperature sensitivity of Tibet EV-C96 strains were assayed on a monolayer of RD cells in 24-well plates. Full-length genome sequencing was performed by a 'primer-walking' strategy, and the evolutionary history of EV-C96 was studied by maximum likelihood analysis.

RESULTS

Strain 2005-T49 grew in all three kinds of cells, and it was not temperature sensitive. In contrast, none of the three cells produced CPE for strain 2012-94H. Phylogenetic analysis of the two Tibetan viruses, other EV-C96 strains, and EV-C prototypes showed that EV-C96 strains were grouped into three clusters (Cluster1-3) based on their VP1 sequences, which may represent three genotypes. Phylogenetic trees based on the P2 and P3 regions highlighted the difference between Chinese EV-C96 strains and the EV-C96 prototype strain BAN-10488. All Chinese strains formed a cluster separate from BAN-10488, which clustered with CV-A1/CV-A22/CV-A19.

CONCLUSIONS

There is genetic variability between EV-C96 strains which suggest that at least few genetic lineages co-exist and there has been some degree of circulation in different geographical regions for some time. Some recombination events must have occurred during EV-C96 evolution as EV-C96 isolates cluster with different EV-C prototype strains in phylogenetic trees in different genomic regions. However, recombination does not seem to have occurred frequently as EV-C96 isolates from different years and locations appear to cluster together in all genomic regions analysed. These findings expand the understanding of the characterization of EV-C96 and are meaningful for the surveillance of the virus.

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.

Discovery of Enterovirus A71-like nonstructural genomes in recent circulating viruses of the Enterovirus A species

Enterovirus A71 (EV-A71) is an important nonpolio enterovirus that causes severe neurological complications. In 1998, Taiwan experienced an EV-A71 outbreak that caused 78 deaths. Since then, periodic epidemics of EV-A71 associated with newly emerging strains have occurred. Several of these strains are known to be recombinant; however, how these strains arose within such a short period of time remains unknown. Here, we sequenced 64 full-length genomes from clinical isolates collected from 2005 to 2016 and incorporated all 91 Taiwanese genomes downloaded from the Virus Pathogen Resource to extensively analyze EV-A71 recombination in Taiwan. We found that the B3 subgenotype was a potential recombinant parent of the EV-A71 C2-like and C4 strains by intratypic recombination. Such B3-similar regions were also found in many cocirculating coxsackieviruses belonging to Enterovirus A species (EV-A) through a series of intertypic recombinations. Therefore, locally enriched outbreaks of cocirculating viruses from different genotypes/serotypes may facilitate recombination. Most recombination breakpoints we found had nonrandom distributions and were located within the region spanning from the boundary of P1 (structural gene) and P2 (nonstructural) to the cis-acting replication element at P2, indicating that specific genome reassembly of structural and nonstructural genes may be subject to natural selection. Through intensive recombination, 11 EV-A71-like signatures (including one in 3A, two in 3C, and eight in 3D) were found to be present in a variety of recently cocirculating EV-A viruses worldwide, suggesting that these viruses may be targets for wide-spectrum antiviral development.

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.

Phylogenetic Characterizations of Highly Mutated EV-B106 Recombinants Showing Extensive Genetic Exchanges with Other EV-B in Xinjiang, China

Human enterovirus B106 (EV-B106) is a new member of the enterovirus B species. To date, only three nucleotide sequences of EV-B106 have been published, and only one full-length genome sequence (the Yunnan strain 148/YN/CHN/12) is available in the GenBank database. In this study, we conducted phylogenetic characterisation of four EV-B106 strains isolated in Xinjiang, China. Pairwise comparisons of the nucleotide sequences and the deduced amino acid sequences revealed that the four Xinjiang EV-B106 strains had only 80.5–80.8% nucleotide identity and 95.4–97.3% amino acid identity with the Yunnan EV-B106 strain, indicating high mutagenicity. Similarity plots and bootscanning analyses revealed that frequent intertypic recombination occurred in all four Xinjiang EV-B106 strains in the non-structural region. These four strains may share a donor sequence with the EV-B85 strain, which circulated in Xinjiang in 2011, indicating extensive genetic exchanges between these strains. All Xinjiang EV-B106 strains were temperature-sensitive. An antibody seroprevalence study against EV-B106 in two Xinjiang prefectures also showed low titres of neutralizing antibodies, suggesting limited exposure and transmission in the population. This study contributes the whole genome sequences of EV-B106 to the GenBank database and provides valuable information regarding the molecular epidemiology of EV-B106 in China.

T-RECs: rapid and large-scale detection of recombination events among different evolutionary lineages of viral genomes

Background

Many computational tools that detect recombination in viruses are not adapted for the ongoing genomic revolution. A computational tool is needed, that will rapidly scan hundreds/thousands of genomes or sequence fragments and detect candidate recombination events that may later be further analyzed with more sensitive and specialized methods.

Results

T-RECs, a Windows based graphical tool, employs pairwise alignment of sliding windows and can perform (i) genotyping, (ii) clustering of new genomes, (iii) detect recent recombination events among different evolutionary lineages, (iv) manual inspection of detected recombination events by similarity plots and (v) annotation of genomic regions.

Conclusions

T-RECs is very effective, as demonstrated by an analysis of 555 Norovirus complete genomes and 2500 sequence fragments, where a recombination hotspot was identified at the ORF1-ORF2 junction.

How much reduction of virus is needed for recycled water: A continuous changing need for assessment?

To ensure the safety of wastewater reuse for irrigation of food crops and drinking water pathogenic viruses must be reduced to levels that pose no significant risk. To achieve this goal minimum reduction of viruses by treatment trains have been suggested. For use of edible crops a 6-log reduction and for production of potable drinking water a 12-log reduction has been suggested. These reductions were based on assuming infective virus concentrations of 10⁵ to 10⁶ per liter. Recent application of molecular methods suggests that some pathogenic viruses may be occurring in concentrations of 10⁷ to 10⁹ per liter. Factors influencing these levels include the development of molecular methods for virus detection, emergence of newly recognized viruses, decrease in per capita water use due to conservation measures, and outbreaks. Since neither cell culture nor molecular methods can assess all the potentially infectious virus in wastewater conservative estimates should be used to assess the virus load in untreated wastewater. This review indicates that an additional 2- to 3-log reduction of viruses above current recommendations may be needed to ensure the safety of recycled water. Information is needed on peak loading of viruses. In addition, more virus groups need to be quantified using better methods of virus quantification, including more accurate methods for measuring viral infectivity in order to better quantify risks from viruses in recycled water.

Poliovirus Polymerase Leu420 Facilitates RNA Recombination and Ribavirin Resistance

RNA recombination is important in the formation of picornavirus species groups and the ongoing evolution of viruses within species groups. In this study, we examined the structure and function of poliovirus polymerase, 3D(pol), as it relates to RNA recombination. Recombination occurs when nascent RNA products exchange one viral RNA template for another during RNA replication. Because recombination is a natural aspect of picornavirus replication, we hypothesized that some features of 3D(pol) may exist, in part, to facilitate RNA recombination. Furthermore, we reasoned that alanine substitution mutations that disrupt 3D(pol)-RNA interactions within the polymerase elongation complex might increase and/or decrease the magnitudes of recombination. We found that an L420A mutation in 3D(pol) decreased the frequency of RNA recombination, whereas alanine substitutions at other sites in 3D(pol) increased the frequency of recombination. The 3D(pol) Leu420 side chain interacts with a ribose in the nascent RNA product 3 nucleotides from the active site of the polymerase. Notably, the L420A mutation that reduced recombination also rendered the virus more susceptible to inhibition by ribavirin, coincident with the accumulation of ribavirin-induced G→A and C→U mutations in viral RNA. We conclude that 3D(pol) Leu420 is critically important for RNA recombination and that RNA recombination contributes to ribavirin resistance.

IMPORTANCE

Recombination contributes to the formation of picornavirus species groups and the emergence of circulating vaccine-derived polioviruses (cVDPVs). The recombinant viruses that arise in nature are occasionally more fit than either parental strain, especially when the two partners in recombination are closely related, i.e., members of characteristic species groups, such as enterovirus species groups A to H or rhinovirus species groups A to C. Our study shows that RNA recombination requires conserved features of the viral polymerase. Furthermore, a polymerase mutation that disables recombination renders the virus more susceptible to the antiviral drug ribavirin, suggesting that recombination contributes to ribavirin resistance. Elucidating the molecular mechanisms of RNA replication and recombination may help mankind achieve and maintain poliovirus eradication.

Whole Genome Sequencing of Enterovirus species C Isolates by High-Throughput Sequencing: Development of Generic Primers

Enteroviruses are among the most common viruses infecting humans and can cause diverse clinical syndromes ranging from minor febrile illness to severe and potentially fatal diseases. Enterovirus species C (EV-C) consists of more than 20 types, among which the three serotypes of polioviruses, the etiological agents of poliomyelitis, are included. Biodiversity and evolution of EV-C genomes are shaped by frequent recombination events. Therefore, identification and characterization of circulating EV-C strains require the sequencing of different genomic regions. A simple method was developed to quickly sequence the entire genome of EV-C isolates. Four overlapping fragments were produced separately by RT-PCR performed with generic primers. The four amplicons were then pooled and purified prior to being sequenced by a high-throughput technique. The method was assessed on a panel of EV-Cs belonging to a wide-range of types. It can be used to determine full-length genome sequences through de novo assembly of thousands of reads. It was also able to discriminate reads from closely related viruses in mixtures. By decreasing the workload compared to classical Sanger-based techniques, this method will serve as a precious tool for sequencing large panels of EV-Cs isolated in cell cultures during environmental surveillance or from patients, including vaccine-derived polioviruses.

A Novel Recombinant Enterovirus Type EV-A89 with Low Epidemic Strength in Xinjiang, China

Enterovirus A89 (EV-A89) is a novel member of the EV-A species. To date, only one full-length genome sequence (the prototype strain) has been published. Here, we report the molecular identification and genomic characterization of a Chinese EV-A89 strain, KSYPH-TRMH22F/XJ/CHN/2011, isolated in 2011 from a contact of an acute flaccid paralysis (AFP) patient during AFP case surveillance in Xinjiang China. This was the first report of EV-A89 in China. The VP1 coding sequence of this strain demonstrated 93.2% nucleotide and 99.3% amino acid identity with the EV-A89 prototype strain. In the P2 and P3 regions, the Chinese EV-A89 strain demonstrated markedly higher identity than the prototype strains of EV-A76, EV-A90, and EV-A91, indicating that one or more recombination events between EV-A89 and these EV-A types might have occurred. Long-term evolution of these EV types originated from the same ancestor provides the spatial and temporal circumstances for recombination to occur. An antibody sero-prevalence survey against EV-A89 in two Xinjiang prefectures demonstrated low positive rates and low titres of EV-A89 neutralization antibody, suggesting limited range of transmission and exposure to the population. This study provides a solid foundation for further studies on the biological and pathogenic properties of EV-A89.

Prevalence and molecular characterizations of enterovirus D68 among children with acute respiratory infection in China between 2012 and 2014

EV-D68 is associated with respiratory tract infections (RTIs). Since its first isolation, EV-D68 has been detected sporadically. However, the US and Canada have experienced outbreaks of EV-D68 infections between August and December 2014. This study aimed to investigate the molecular epidemiology and clinical characteristics of EV-D68 in Chongqing, Southwestern China. From January 2012 to November 2014, 1876 nasopharyngeal aspirate specimens (NPAs) were collected from hospitalized children with RTIs. Among the 1876 NPAs, EV-D68 was detected in 19 samples (1.0%, 19/1876). Of these, 13 samples were detected in September and October 2014 (9.8%, 13/132). Phylogenetic analysis showed that all 13 strains detected in the 2014 Chongqing had high homology with the main strains of the 2014 US outbreak. Among the children with EV-D68 infection, 13 (68%) had a history of recurrent wheezing. A total of 13 children had a discharge diagnosis of asthma. Of these, 11 children were diagnosed with acute asthma exacerbation. EV-D68 was the predominant pathogen that evoked asthma exacerbation in September and October 2014. In conclusion, our results found that a history of recurrent wheezing may be a risk factor for the detection of EV-D68 and viral-induced asthma exacerbation may be a clinical feature of EV-D68 infection.