Recombination in circulating enteroviruses

Evolutionary and Genetic Recombination Analyses of Coxsackievirus A6 Variants Associated with Hand, Foot, and Mouth Disease Outbreaks in Thailand between 2019 and 2022

Coxsackievirus (CV)-A6 infections cause hand, foot, and mouth disease (HFMD) in children and adults. Despite the serious public health threat presented by CV-A6 infections, our understanding of the mechanisms by which new CV-A6 strains emerge remains limited. This study investigated the molecular epidemiological trends, evolutionary dynamics, and recombination characteristics of CV-A6-associated HFMD in Thailand between 2019 and 2022. In the HFMD patient samples collected during the 4-year study period, we identified enterovirus (EV) RNA in 368 samples (48.7%), of which CV-A6 (23.7%) was the predominant genotype, followed by CV-A4 (6%), EV-A71 (3.7%), and CV-A16 (3.4%). According to the partial viral protein (VP) 1 sequences, all these CV-A6 strains belonged to the D3 clade. Based on the viral-RNA-dependent RNA polymerase (RdRp) gene, four recombinant forms (RFs), RF-A (147, 84.5%), RF-N (11, 6.3%), RF-H (1, 0.6%), and newly RF-Y (15, 8.6%), were identified throughout the study period. Results from the similarity plot and bootscan analyses revealed that the 3D polymerase (3Dpol) region of the D3/RF-Y subclade consists of sequences highly similar to CV-A10. We envisage that the epidemiological and evolutionarily insights presented in this manuscript will contribute to the development of vaccines to prevent the spread of CV-A6 infection.

Genetic and Cross Neutralization Analyses of Coxsackievirus A16 Circulating in Taiwan from 1998 to 2021 Suggest Dominant Genotype B1 can Serve as Vaccine Candidate

Coxsackievirus A16 (CVA16) is well known for causing hand-foot-and-mouth disease (HFMD) and outbreaks were frequently reported in Taiwan in the past twenty years. The epidemiology and genetic variations of CVA16 in Taiwan from 1998 to 2021 were analyzed in this study. CVA16 infections usually occurred in early summer and early winter, and showed increased incidence in 1998, 2000-2003, 2005, 2007-2008, and 2010 in Taiwan. Little or no CVA16 was detected from 2017 to 2021. CVA16 infection was prevalent in patients between 1 to 3 years old. A total of 69 isolates were sequenced. Phylogenetic analysis based on the VP1 region showed that CVA16 subgenotype B1 was dominantly isolated in Taiwan from 1998 to 2019, and B2 was identified only from isolates collected in 1999 and 2000. There was a high frequency of synonymous mutations in the amino acid sequences of the VP1 region among CVA16 isolates, with the exception of position 145 which showed positive selection. The recombination analysis of the whole genome of CVA16 isolates indicated that the 5'-untranslated region and the non-structural protein region of CVA16 subgenotype B1 were recombined with Coxsackievirus A4 (CVA4) and enterovirus A71 (EVA71) genotype A, respectively. The recombination pattern of subgenotype B2 was similar to B1, however, the 3D region was similar to EVA71 genotype B. Cross-neutralization among CVA16 showed that mouse antisera from various subgenotypes viruses can cross-neutralize different genotype with high neutralizing antibody titers. These results suggest that the dominant CVA16 genotype B1 can serve as a vaccine candidate for CVA16.

Molecular characterization of a porcine sapelovirus strain isolated in China

Porcine sapelovirus (PSV) infections have been associated with a wide spectrum of symptoms, ranging from asymptomatic infection to clinical signs including diarrhoea, pneumonia, reproductive disorders, and polioencephalomyelitis. Although it has a global distribution, there have been relatively few studies on PSV in domestic animals. We isolated a PSV strain, SHCM2019, from faecal specimens from swine, using PK-15 cells. To investigate its molecular characteristics and pathogenicity, the genomic sequence of strain SHCM2019 was analysed, and clinical manifestations and pathological changes occurring after inoculation of neonatal piglets were observed. The virus isolated using PK-15 cells was identified as PSV using RT-PCR, transmission electron microscopy (TEM), and immunofluorescence assay (IFA). Sequencing results showed that the full-length genome of the SHCM2019 strain was 7,567 nucleotides (nt) in length, including a 27-nucleotide poly(A) tail. Phylogenetic analysis demonstrated that this virus was a PSV isolate belonging to the Chinese strain cluster. Recombination analysis indicated that there might be a recombination breakpoint upstream of the 3D region of the genome. Pathogenicity experiments demonstrated that the virus isolate could cause diarrhoea and pneumonia in piglets. In breif, a recombinant PSV strain, SHCM2019, was isolated and shown to be pathogenic. Our results may provide a reference for future research on the pathogenic mechanism and evolutionary characteristics of PSV.

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.

Deformed wing virus: using reverse genetics to tackle unanswered questions about the most important viral pathogen of honey bees

Deformed wing virus (DWV) is the most important viral pathogen of honey bees. It usually causes asymptomatic infections but, when vectored by the ectoparasitic mite Varroa destructor, it is responsible for the majority of overwintering colony losses globally. Although DWV was discovered four decades ago, research has been hampered by the absence of an in vitro cell culture system or the ability to culture pure stocks of the virus. The recent developments of reverse genetic systems for DWV go some way to addressing these limitations. They will allow the investigation of specific questions about strain variation, host tropism and pathogenesis to be answered, and are already being exploited to study tissue tropism and replication in Varroa and non-Apis pollinators. Three areas neatly illustrate the advances possible with reverse genetic approaches: (i) strain variation and recombination, in which reverse genetics has highlighted similarities rather than differences between virus strains; (ii) analysis of replication kinetics in both honey bees and Varroa, in studies that likely explain the near clonality of virus populations often reported; and (iii) pathogen spillover to non-Apis pollinators, using genetically tagged viruses to accurately monitor replication and infection.

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.

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.

The heterogeneous pathogenesis of type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) results from the destruction of pancreatic β-cells that is mediated by the immune system. Multiple genetic and environmental factors found in variable combinations in individual patients are involved in the development of T1DM. Genetic risk is defined by the presence of particular allele combinations, which in the major susceptibility locus (the HLA region) affect T cell recognition and tolerance to foreign and autologous molecules. Multiple other loci also regulate and affect features of specific immune responses and modify the vulnerability of β-cells to inflammatory mediators. Compared with the genetic factors, environmental factors that affect the development of T1DM are less well characterized but contact with particular microorganisms is emerging as an important factor. Certain infections might affect immune regulation, and the role of commensal microorganisms, such as the gut microbiota, are important in the education of the developing immune system. Some evidence also suggests that nutritional factors are important. Multiple islet-specific autoantibodies are found in the circulation from a few weeks to up to 20 years before the onset of clinical disease and this prediabetic phase provides a potential opportunity to manipulate the islet-specific immune response to prevent or postpone β-cell loss. The latest developments in understanding the heterogeneity of T1DM and characterization of major disease subtypes might help in the development of preventive treatments.

Yellow Fever: Integrating Current Knowledge with Technological Innovations to Identify Strategies for Controlling a Re-Emerging Virus

Yellow fever virus (YFV) represents a re-emerging zoonotic pathogen, transmitted by mosquito vectors to humans from primate reservoirs. Sporadic outbreaks of YFV occur in endemic tropical regions, causing a viral hemorrhagic fever (VHF) associated with high mortality rates. Despite a highly effective vaccine, no antiviral treatments currently exist. Therefore, YFV represents a neglected tropical disease and is chronically understudied, with many aspects of YFV biology incompletely defined including host range, host–virus interactions and correlates of host immunity and pathogenicity. In this article, we review the current state of YFV research, focusing on the viral lifecycle, host responses to infection, species tropism and the success and associated limitations of the YFV-17D vaccine. In addition, we highlight the current lack of available treatments and use publicly available sequence and structural data to assess global patterns of YFV sequence diversity and identify potential drug targets. Finally, we discuss how technological advances, including real-time epidemiological monitoring of outbreaks using next-generation sequencing and CRISPR/Cas9 modification of vector species, could be utilized in future battles against this re-emerging pathogen which continues to cause devastating disease.

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.

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.

Higher abundance of enterovirus A species in the gut of children with islet autoimmunity

Enteroviruses (EVs) are prime candidate environmental triggers of islet autoimmunity (IA), with potential as vaccine targets for type 1 diabetes prevention. However, the use of targeted virus detection methods and the selective focus on EVs by most studies increases the risk for substantial investigation bias and an overestimated association between EV and type 1 diabetes. Here we performed comprehensive virome-capture sequencing to examine all known vertebrate-infecting viruses without bias in 182 specimens (faeces and plasma) collected before or at seroconversion from 45 case children with IA and 48 matched controls. From >2.6 billion reads, 28 genera of viruses were detected and 62% of children (58/93) were positive for ≥1 vertebrate-infecting virus. We identified 129 viruses as differentially abundant between the gut of cases and controls, including 5 EV-A types significantly more abundant in the cases. Our findings further support EV’s hypothesised contribution to IA and corroborate the proposal that viral load may be an important parameter in disease pathogenesis. Furthermore, our data indicate a previously unrecognised association of IA with higher EV-A abundance in the gut of children and provide a catalog of viruses to be interrogated further to determine a causal link between virus infection and type 1 diabetes.

Impact of Enterovirus Molecular Assay Turnaround Time on Hospitalization Length During an Echovirus 30 Meningitis Outbreak, France, Fall 2014

Background:

The impact of Enterovirus Real Time-Polymerase Chain Reaction assay (EV RT-PCR) on hospitalization lengths of patients with aseptic meningitis has been investigated but the impact of early EV RT-PCR results released on time before patient discharge remains unclear during Echovirus meningitis outbreaks.

Objective:

To assess a potential correlation between EV RT-PCR turn-around time and hospitalization lengths during an Echovirus meningitis outbreak.

Method:

Eighteen patients demonstrating a positive EV RT-PCR assay performed on Cerebrospinal Fluid (CSF) samples collected between October 1st 2014 and December 31st 2014 were retrospectively included. Viral protein 1 (VP1) gene region was amplified and sequenced using a classical Sanger sequencing reaction. Clinical data were retrospectively collected from patient’s records. Quantitative variables expressed as median values and ranges were compared using Mann Whitney U test. Correlations were performed using simple regression analysis.

Results:

Phylogenetic VP1 sequence analyses identified that the outbreak was related to an Echovirus 30 strain in 7 out of the 10 cases with available sequencing data. The three remaining sequences analyses evidenced Echovirus 14, 9 and 7 strains. Hospitalization length was statistically shorter in children without comorbidity (n=5) than in adult patients (n=10) or neonates and children with comorbidity (n=3) (p=0.003 and 0.01 respectively), whereas EV RT-PCR turnaround time was not statistically different between these groups. Correlation between hospitalization length and EV RT-PCR turnaround time was poor (R2=0.06), especially in adults (R2=0.01)

Conclusion:

Our data indicated that EV RT-PCR turnaround time was not correlated to hospitalization length during a short Echovirus meningitis outbreak.

Complete genome analysis demonstrates multiple introductions of enterovirus 71 and coxsackievirus A16 recombinant strains into Thailand during the past decade

Hand, foot, and mouth disease (HFMD) caused by enteroviruses remains a public health threat, particularly in the Asia-Pacific region during the past two decades. Moreover, the introduction of multiple subgenotypes and the emergence of recombinant viruses is of epidemiological importance. Based on either the full genome or VP1 sequences, 32 enteroviruses (30 from HFMD patients, 1 from an encephalitic patient, and 1 from an asymptomatic contact case) isolated in Thailand between 2006 and 2014 were identified as 25 enterovirus 71 (EV71) isolates (comprising 20 B5, 1 C2, 2 C4a, and 2 C4b subgenotypes) and 7 coxsackievirus A16 (CA16) isolates (comprising 6 B1a and 1 B1b subgenotypes). The EV71 subgenotype C4b was introduced into Thailand for the first time in 2006 and was replaced by subgenotype C4a strains in 2009. Phylogenetic, similarity plot and bootscan analyses of the complete viral genomes identified 12 recombinant viruses among the 32 viral isolates. Only one EV71-B5 isolate out of 20 was a recombinant virus with one region of intratypic or intertypic recombination, while all four EV71-C4 isolates were recombinant viruses having undergone double recombination, and all seven CA16 isolates were recombinant viruses. The recombination breakpoints of these recombinants are located solely within the P2 and P3 regions. Surveillance for circulating strains and subgenotype replacement are important with respect to molecular epidemiology and the selection of the upcoming EV71 vaccine. In addition, the clinical importance of recombinant viruses needs to be further explored.

Contemporary Circulating Enterovirus D68 Strains Have Acquired the Capacity for Viral Entry and Replication in Human Neuronal Cells

Since the EV-D68 outbreak during the summer of 2014, evidence of a causal link to a type of limb paralysis (AFM) has been mounting. In this article, we describe a neuronal cell culture model (SH-SY5Y cells) in which a subset of contemporary 2014 outbreak strains of EV-D68 show infectivity in neuronal cells, or neurotropism. We confirmed the difference in neurotropism in vitro using primary human neuron cell cultures and in vivo with a mouse paralysis model. Using the SH-SY5Y cell model, we determined that a barrier to viral entry is at least partly responsible for neurotropism. SH-SY5Y cells may be useful in determining if specific EV-D68 genetic determinants are associated with neuropathogenesis, and replication in this cell line could be used as rapid screening tool for identification of neurotropic EV-D68 strains. This may assist with better understanding of pathogenesis and epidemiology and with the development of potential therapies.

Enterovirus D68 (EV-D68) has historically been associated with respiratory illnesses. However, in the summers of 2014 and 2016, EV-D68 outbreaks coincided with a spike in polio-like acute flaccid myelitis/paralysis (AFM/AFP) cases. This raised concerns that EV-D68 could be the causative agent of AFM during these recent outbreaks. To assess the potential neurotropism of EV-D68, we utilized the neuroblastoma-derived neuronal cell line SH-SY5Y as a cell culture model to determine if differential infection is observed for different EV-D68 strains. In contrast to HeLa and A549 cells, which support viral infection of all EV-D68 strains tested, SH-SY5Y cells only supported infection by a subset of contemporary EV-D68 strains, including isolates from the 2014 outbreak. Viral replication and infectivity in SH-SY5Y were assessed using multiple assays: virus production, cytopathic effects, cellular ATP release, and VP1 capsid protein production. Similar differential neurotropism was also observed in differentiated SH-SY5Y cells, primary human neuron cultures, and a mouse paralysis model. Using the SH-SY5Y cell culture model, we determined that barriers to viral binding and entry were at least partly responsible for the differential infectivity phenotype. Transfection of genomic RNA into SH-SY5Y generated virions for all EV-D68 isolates, but only a single round of replication was observed from strains that could not directly infect SH-SY5Y. In addition to supporting virus replication and other functional studies, this cell culture model may help identify the signatures of virulence to confirm epidemiological associations between EV-D68 strains and AFM and allow for the rapid identification and characterization of emerging neurotropic strains.

Whole Genome Sequencing of Enteroviruses Species A to D by High-Throughput Sequencing: Application for Viral Mixtures

Human enteroviruses (EV) consist of more than 100 serotypes classified within four species for enteroviruses (EV-A to -D) and three species for rhinoviruses, which have been implicated in a variety of human illnesses. Being able to simultaneously amplify the whole genome and identify enteroviruses in samples is important for studying the viral diversity in different geographical regions and populations. It also provides knowledge about the evolution of these viruses. Therefore, we developed a rapid, sensitive method to detect and genetically classify all human enteroviruses in mixtures. Strains of EV-A (15), EV-B (40), EV-C (20), and EV-D (2) viruses were used in addition to 20 supernatants from RD cells infected with stool extracts or sewage concentrates. Two overlapping fragments were produced using a newly designed degenerated primer targeting the conserved CRE region for enteroviruses A-D and one degenerated primer set designed to specifically target the conserved region for each enterovirus species (EV-A to -D). This method was capable of sequencing the full genome for all viruses except two, for which nearly 90% of the genome was sequenced. This method also demonstrated the ability to discriminate, in both spiked and unspiked mixtures, the different enterovirus types present.

Mechanisms and consequences of positive-strand RNA virus recombination

Genetic recombination in positive-strand RNA viruses is a significant evolutionary mechanism that drives the creation of viral diversity by the formation of novel chimaeric genomes. The process and its consequences, for example the generation of viruses with novel phenotypes, has historically been studied by analysis of the end products. More recently, with an appreciation that there are both replicative and non-replicative mechanisms at work, and with new approaches and techniques to analyse intermediate products, the viral and cellular factors that influence the process are becoming understood. The major influence on replicative recombination is the fidelity of viral polymerase, although RNA structures and sequences may also have an impact. In replicative recombination the viral polymerase is necessary and sufficient, although roles for other viral or cellular proteins may exist. In contrast, non-replicative recombination appears to be mediated solely by cellular components. Despite these insights, the relative importance of replicative and non-replicative mechanisms is not clear. Using single-stranded positive-sense RNA viruses as exemplars, we review the current state of understanding of the processes and consequences of recombination.

Exploratory re-encoding of yellow fever virus genome: new insights for the design of live-attenuated viruses

Virus attenuation by genome re-encoding is a pioneering approach for generating effective live-attenuated vaccine candidates. Its core principle is to introduce a large number of synonymous substitutions into the viral genome to produce stable attenuation of the targeted virus. Introduction of large numbers of mutations has also been shown to maintain stability of the attenuated phenotype by lowering the risk of reversion and recombination of re-encoded genomes. Identifying mutations with low fitness cost is pivotal as this increases the number that can be introduced and generates more stable and attenuated viruses. Here, we sought to identify mutations with low deleterious impact on the in vivo replication and virulence of yellow fever virus (YFV). Following comparative bioinformatic analyses of flaviviral genomes, we categorised synonymous transition mutations according to their impact on CpG/UpA composition and secondary RNA structures. We then designed seventeen re-encoded viruses with 100–400 synonymous mutations in the NS2A-to-NS4B coding region of YFV Asibi and Ap7M (hamster-adapted) genomes. Each virus contained a panel of synonymous mutations designed according to the above categorisation criteria. The replication and fitness characteristics of parent and re-encoded viruses were compared in vitro using cell culture competition experiments. In vivo laboratory hamster models were also used to compare relative virulence and immunogenicity characteristics. Most of the re-encoded strains showed no decrease in replicative fitness in vitro. However, they showed reduced virulence and, in some instances, decreased replicative fitness in vivo. Importantly, the most attenuated of the re-encoded strains induced robust, protective immunity in hamsters following challenge with Ap7M, a virulent virus. Overall, the introduction of transitions with no or a marginal increase in the number of CpG/UpA dinucleotides had the mildest impact on YFV replication and virulence in vivo. Thus, this strategy can be incorporated in procedures for the finely tuned creation of substantially re-encoded viral genomes.

Metagenomic analysis of Sichuan takin fecal sample viromes reveals novel enterovirus and astrovirus

The Sichuan takin inhabits the bamboo forests in the Eastern Himalayas and is considered as a national treasure of China with the highest legal protection and conservation status considered as vulnerable according to The IUCN Red List of Threatened Species. In this study, fecal samples of 71 Sichuan takins were pooled and deep sequenced. Among the 103,553 viral sequences, 21,961 were assigned to mammalian viruses. De novo assembly revealed genomes of an enterovirus and an astrovirus and contigs of circoviruses and genogroup I picobirnaviruses. Complete genome sequencing and phylogenetic analysis showed that Sichuan takin enterovirus is a novel serotype/genotype of the species Enterovirus G, with evidence of recombination. Sichuan takin astrovirus is a new subtype of bovine astrovirus, probably belonging to a new genogroup in the genus Mamastrovirus. Further studies will reveal whether these viruses can also be found in Mishmi takin and Shaanxi takin and their pathogenic potentials.

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.

Comparative Genomic Analysis for Genetic Variation in Sacbrood Virus of Apis cerana and Apis mellifera Honeybees From Different Regions of Vietnam

Sacbrood virus (SBV) is one of the most common viral infections of honeybees. The entire genome sequence for nine SBV infecting honeybees, Apis cerana and Apis mellifera, in Vietnam, namely AcSBV-Viet1, AcSBV-Viet2, AcSBV-Viet3, AmSBV-Viet4, AcSBV-Viet5, AmSBV-Viet6, AcSBV-Viet7, AcSBV-Viet8, and AcSBV-Viet9, was determined. These sequences were aligned with seven previously reported complete genome sequences of SBV from other countries, and various genomic regions were compared. The Vietnamese SBVs (VN-SBVs) shared 91-99% identity with each other, and shared 89-94% identity with strains from other countries. The open reading frames (ORFs) of the VN-SBV genomes differed greatly from those of SBVs from other countries, especially in their VP1 sequences. The AmSBV-Viet6 and AcSBV-Viet9 genome encodes 17 more amino acids within this region than the other VN-SBVs. In a phylogenetic analysis, the strains AmSBV-Viet4, AcSBV-Viet2, and AcSBV-Viet3 were clustered in group with AmSBV-UK, AmSBV-Kor21, and AmSBV-Kor19 strains. Whereas, the strains AmSBV-Viet6 and AcSBV-Viet7 clustered separately with the AcSBV strains from Korea and AcSBV-VietSBM2. And the strains AcSBV-Viet8, AcSBV-Viet1, AcSBV-Viet5, and AcSBV-Viet9 clustered with the AcSBV-India, AcSBV-Kor and AcSBV-VietSBM2. In a Simplot graph, the VN-SBVs diverged stronger in their ORF regions than in their 5' or 3' untranslated regions. The VN-SBVs possess genetic characteristics which are more similar to the Asian AcSBV strains than to AmSBV-UK strain. Taken together, our data indicate that host specificity, geographic distance, and viral cross-infections between different bee species may explain the genetic diversity among the VN-SBVs in A. cerana and A. mellifera and other SBV strains.

First Report of a Fatal Case Associated with EV-D68 Infection in Hong Kong and Emergence of an Interclade Recombinant in China Revealed by Genome Analysis

A fatal case associated with enterovirus D68 (EV-D68) infection affecting a 10-year-old boy was reported in Hong Kong in 2014. To examine if a new strain has emerged in Hong Kong, we sequenced the partial genome of the EV-D68 strain identified from the fatal case and the complete VP1, and partial 5′UTR and 2C sequences of nine additional EV-D68 strains isolated from patients in Hong Kong. Sequence analysis indicated that a cluster of strains including the previously recognized A2 strains should belong to a separate clade, clade D, which is further divided into subclades D1 and D2. Among the 10 EV-D68 strains, 7 (including the fatal case) belonged to the previously described, newly emerged subclade B3, 2 belonged to subclade B1, and 1 belonged to subclade D1. Three EV-D68 strains, each from subclades B1, B3, and D1, were selected for complete genome sequencing and recombination analysis. While no evidence of recombination was noted among local strains, interclade recombination was identified in subclade D2 strains detected in mainland China in 2008 with VP2 acquired from clade A. This study supports the reclassification of subclade A2 into clade D1, and demonstrates interclade recombination between clades A and D2 in EV-D68 strains from China.

VIRO-TypeNed, systematic molecular surveillance of enteroviruses in the Netherlands between 2010 and 2014

VIRO-TypeNed is a collaborative molecular surveillance platform facilitated through a web-based database. Genetic data in combination with epidemiological, clinical and patient data are shared between clinical and public health laboratories, as part of the surveillance underpinning poliovirus eradication. We analysed the combination of data submitted from 2010 to 2014 to understand circulation patterns of non-polio enteroviruses (NPEV) of public health relevance. Two epidemiological patterns were observed based on VIRO-TypeNed data and classical surveillance data dating back to 1996: (i) endemic cyclic, characterised by predictable upsurges/outbreaks every two to four years, and (ii) epidemic, where rare virus types caused upsurges/outbreaks. Genetic analysis suggests continuous temporal displacement of virus lineages due to the accumulation of (silent) genetic changes. Non-synonymous changes in the antigenic B/C loop suggest antigenic diversification, which may affect population susceptibility. Infections were frequently detected at an age under three months and at an older, parenting age (25–49 years) pointing to a distinct role of immunity in the circulation patterns. Upsurges were detected in the summer and winter which can promote increased transmissibility underlying new (cyclic) upsurges and requires close monitoring. The combination of data provide a better understanding of NPEV circulation required to control and curtail upsurges and outbreaks.

Characterization of enterovirus 71 infection and associated outbreak of Hand, Foot, and Mouth Disease in Shawo of China in 2012

Infection of enterovirus 71 (EV71) and associated hand, foot, and mouth disease (HFMD) are recognized as emerging public health issues worldwide. Hundreds of thousands of children are annually infected with EV71 and develop HFMD in China alone. Studies of EV71 infection are critical to the treatment and prevention of the associated HFMD outbreaks. In this report, we studied an outbreak of 105 HFMD cases in Shawo Township of China between September to October 2012. More than 90% of cases were children younger than 9 years old, with over 50% of cases aged 3–6 years old. Laboratory studies detected a high prevalence of EV71 and suggested EV71 as the most common enterovirus causing HFMD in Shawo. Sequencing analysis showed that the EV71 strains from Shawo belong to the C4 subgenotype, and are phylogenetically more related to those from the distant city of Nanchang than those from the nearby city of Wuhan with distinct variations. More girls were found to be associated with EV71 in Shawo whereas more boys were associated with EV71 in Wuhan and Nanchang. Our studies further the understanding of the molecular epidemiological features of HFMD and infection by enteroviruses in China.

Complete Genome Analysis of an Enterovirus EV-B83 Isolated in China

Enterovirus B83 (EV-B83) is a recently identified member of enterovirus species B. It is a rarely reported serotype and up to date, only the complete genome sequence of the prototype strain from the United States is available. In this study, we describe the complete genomic characterization of an EV-B83 strain 246/YN/CHN/08HC isolated from a healthy child living in border region of Yunnan Province, China in 2008. Compared with the prototype strain, it had 79.6% similarity in the complete genome and 78.9% similarity in the VP1 coding region, reflecting the great genetic divergence among them. VP1-coding region alignment revealed it had 77.2–91.3% with other EV-B83 sequences available in GenBank. Similarity plot analysis revealed it had higher identity with several other EV-B serotypes than the EV-B83 prototype strain in the P2 and P3 coding region, suggesting multiple recombination events might have occurred. The great genetic divergence with previously isolated strains and the extremely rare isolation suggest this serotype has circulated at a low epidemic strength for many years. This is the first report of complete genome of EV-B83 in China.

Homology differences between complete Sacbrood virus genomes from infected Apis mellifera and Apis cerana honeybees in Korea

Sacbrood virus (SBV) represents a serious threat to the health of managed honeybees. We determined four complete SBV genomic sequences (AmSBV-Kor1, AmSBV-Kor2, AcSBV-Kor3, and AcSBV-Kor4) isolated from Apis mellifera and Apis cerana in various regions of South Korea. A phylogenetic tree was constructed from the complete genomic sequences of these Korean SBVs (KSBVs) and 21 previously reported SBV sequences from other countries. Three KSBVs (not AmSBV-Kor1) clustered with previously reported Korean genomes, but separately from SBV genomes from other countries. The KSBVs shared 90-98 % identity, and 89-97 % identity with the genomes from other countries. AmSBV-Kor1 was least similar (~90 % identity) to the other KSBVs, and was most similar to previously reported strains AmSBV-Kor21 (97 %) and AmSBV-UK (93 %). Phylogenetic analysis of the partial VP1 region sequences indicated that SBVs clustered by host species and country of origin. The KSBVs were aligned with nine previously reported complete SBV genomes and compared. The KSBVs were most different from the other genomes at the end of the 5' untranslated region and in the entire open reading frame. A SimPlot graph of the VP1 region confirmed its high variability, especially between the SBVs infecting A. mellifera and A. cerana. In this genomic region, SBVs from A. mellifera species contain an extra continuous 51-nucleotide sequence relative to the SBVs from A. cerana. This genomic diversity may reflect the adaptation of SBV to specific hosts, viral cross-infections, and the spatial distances separating the KSBVs from other SBVs.

Precise genotyping and recombination detection of Enterovirus

Enteroviruses (EV) with different genotypes cause diverse infectious diseases in humans and mammals. A correct EV typing result is crucial for effective medical treatment and disease control; however, the emergence of novel viral strains has impaired the performance of available diagnostic tools. Here, we present a web-based tool, named EVIDENCE (EnteroVirus In DEep conception, http://symbiont.iis.sinica.edu.tw/evidence), for EV genotyping and recombination detection. We introduce the idea of using mixed-ranking scores to evaluate the fitness of prototypes based on relatedness and on the genome regions of interest. Using phylogenetic methods, the most possible genotype is determined based on the closest neighbor among the selected references. To detect possible recombination events, EVIDENCE calculates the sequence distance and phylogenetic relationship among sequences of all sliding windows scanning over the whole genome. Detected recombination events are plotted in an interactive figure for viewing of fine details. In addition, all EV sequences available in GenBank were collected and revised using the latest classification and nomenclature of EV in EVIDENCE. These sequences are built into the database and are retrieved in an indexed catalog, or can be searched for by keywords or by sequence similarity. EVIDENCE is the first web-based tool containing pipelines for genotyping and recombination detection, with updated, built-in, and complete reference sequences to improve sensitivity and specificity. The use of EVIDENCE can accelerate genotype identification, aiding clinical diagnosis and enhancing our understanding of EV evolution.

Transmission and Demographic Dynamics of Coxsackievirus B1

The infectious activity of coxsackievirus B1 (CV-B1) in Taiwan was high from 2008 to 2010, following an alarming increase in severe neonate disease in the United States (US). To examine the relationship between CV-B1 strains isolated in Taiwan and those from other parts of the world, we performed a phylodynamic study using VP1 and partial 3Dpol (414 nt) sequences from 22 strains of CV-B1 isolated in Taiwan (1989-2010) and compared them to sequences from strains isolated worldwide. Phylogenetic trees were constructed by neighbor-joining, maximum likelihood, and Bayesian Monte Carlo Markov Chain methods. Four genotypes (GI-IV) in the VP1 region of CV-B1 and three genotypes (GA-C) in the 3Dpol region of enterovirus B were identified and had high support values. The phylogenetic analysis indicates that the GI and GIII strains in VP1 were geographically distributed in Taiwan (1993-1994) and in India (2007-2009). On the other hand, the GII and GIV strains appear to have a wider spatiotemporal distribution and ladder-like topology A stair-like phylogeny was observed in the VP1 region indicating that the phylogeny of the virus may be affected by different selection pressures in the specified regions. Further, most of the GI and GII strains in the VP1 tree were clustered together in GA in the 3D tree, while the GIV strains diverged into GB and GC. Taken together, these data provide important insights into the population dynamics of CV-B1 and indicate that incongruencies in specific gene regions may contribute to spatiotemporal patterns of epidemicity for this virus.

Phylogenetic analysis of the major causative agents of hand, foot and mouth disease in Suzhou City, Jiangsu province, China, in 2012-2013

Hand, foot and mouth disease (HFMD) is a serious public health problem that has emerged over the past several decades. Pathogen detection by the Chinese national HFMD surveillance system has focused mainly on enterovirus 71 (EV71) and coxsackievirus A16 (CA16). Therefore, epidemiological information regarding the other causative enteroviruses is limited. To identify the pandemic enterovirus in Suzhou, Jiangsu province, China, clinical samples from patients with HFMD were collected from 2012 to 2013 and analyzed. The results revealed that CA16 was the most dominant HFMD pathogen in 2012, whereas CA6 and CA10 were the dominant pathogens in 2013. Phylogenetic analysis revealed that the C4a sub-genogroup of EV71 and the B1a and B1b sub-genogroups of CA16 continued to evolve and circulate in Suzhou. The CA6 strains were assigned to six genotypes (A–F) and the CA10 strains were assigned to seven genotypes (A–G), with clear geographical and temporal distributions. All of the CA6 strains in Suzhou belonged to genogroup F, and there were several lineages circulating in Suzhou. All of the CA10 strains in Suzhou belonged to genogroup G, and they had the same genetic origin. Co-infections of EV71/CA16 and CA6/CA10 were found in the samples, and bootscan analysis of 5′-untranslated regions (UTRs) revealed that some CA16 strains in Suzhou had genetic recombination with EV71. This property might allow CA16 to alter its evolvability and circulating ability. This study underscores the need for surveillance of CA6 and CA10 in the Yangtze River Delta and East China.

Efficient replication of recombinant Enterovirus B types, carrying different P1 genes in the coxsackievirus B5 replicative backbone

Recombination is an important feature in the evolution of the Enterovirus genus. Phylogenetic studies of enteroviruses have revealed that the capsid genomic region (P1) is type specific, while the parts of the genome coding for the non-structural proteins (P2-P3) are species specific. Hence, the genome may be regarded as consisting of two modules that evolve independently. In this study, it was investigated whether the non-structural coding part of the genome in one type could support replication of a virus with a P1 region from another type of the same species. A cassette vector (pCas) containing a full-length cDNA copy of coxsackievirus B5 (CVB5) was used as a replicative backbone. The P1 region of pCas was replaced with the corresponding part from coxsackievirus B3 Nancy (CVB3N), coxsackievirus B6 Schmitt (CVB6S), and echovirus 7 Wallace (E7W), all members of the Enterovirus B species. The replication efficiency after transfection with clone-derived in vitro transcribed RNA was studied and compared with that of pCas. All the recombinant viruses replicated with similar efficiencies and showed threshold cycle (Ct) values, tissue culture infectivity dose 50 %, and plaque-forming unit titers comparable to viruses generated from the pCas construct. In addition to this, a clone without the P1 region was also constructed, and Western Blot and immunofluorescence staining analysis showed that the viral genome could be translated and replicated despite the lack of the structural protein-coding region. To conclude, the replicative backbone of the CVB5 cassette vector supports replication of intraspecies constructs with P1 regions derived from other members of the Enterovirus B species. In addition to this, the replicative backbone can be both translated and replicated without the presence of a P1 region.

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

Human enteroviruses (EV) belong to the Picornaviridae family and are among the most common viruses infecting humans. They consist of up to 100 immunologically and genetically distinct types: polioviruses, coxsackieviruses A and B, echoviruses, and the more recently characterized 43 EV types. Frequent recombinations and mutations in enteroviruses have been recognized as the main mechanisms for the observed high rate of evolution, thus enabling them to rapidly respond and adapt to new environmental challenges. The first signs of genetic exchanges between enteroviruses came from polioviruses many years ago, and since then recombination has been recognized, along with mutations, as the main cause for reversion of vaccine strains to neurovirulence. More recently, non-polio enteroviruses became the focus of many studies, where recombination was recognized as a frequent event and was correlated with the appearance of new enterovirus lineages and types. The accumulation of multiple inter- and intra-typic recombination events could also explain the series of successive emergences and disappearances of specific enterovirus types that could in turn explain the epidemic profile of circulation of several types. This review focuses on recombination among human non-polio enteroviruses from all four species (EV-A, EV-B, EV-C, and EV-D) and discusses the recombination effects on enterovirus epidemiology and evolution.

Molecular evolution of VP3, VP1, 3C(pro) and 3D(pol) coding regions in coxsackievirus group A type 24 variant isolates from acute hemorrhagic conjunctivitis in 2011 in Okinawa, Japan

A large acute hemorrhagic conjunctivitis (AHC) outbreak occurred in 2011 in Okinawa Prefecture in Japan. Ten strains of coxsackievirus group A type 24 variant (CA24v) were isolated from patients with AHC and full sequence analysis of the VP3, VP1, 3C(pro) and 3D(pol) coding regions performed. To assess time-scale evolution, phylogenetic analysis was performed using the Bayesian Markov chain Monte Carlo method. In addition, similarity plots were constructed and pairwise distance (p-distance) and positive pressure analyses performed. A phylogenetic tree based on the VP1 coding region showed that the present strains belong to genotype 4 (G4). In addition, the present strains could have divided in about 2010 from the same lineages detected in other countries such as China, India and Australia. The mean rates of molecular evolution of four coding regions were estimated at about 6.15 to 7.86 × 10(-3) substitutions/site/year. Similarity plot analyses suggested that nucleotide similarities between the present strains and a prototype strain (EH24/70 strain) were 0.77-0.94. The p-distance of the present strains was relatively short (<0.01). Only one positive selected site (L25H) was identified in the VP1 protein. These findings suggest that the present CA24v strains causing AHC are genetically related to other AHC strains with rapid evolution and emerged in around 2010.

Long-term immunogenicity studies of formalin-inactivated enterovirus 71 whole-virion vaccine in macaques

Enterovirus 71 (EV71) has caused epidemics of hand, foot and mouth diseases in Asia during the past decades and no vaccine is available. A formalin-inactivated EV71 candidate vaccine (EV71vac) based on B4 subgenotype has previously been developed and found to elicit strong neutralizing antibody responses in mice and humans. In this study, we evaluated the long-term immunogenicity and safety of this EV71vac in a non-human primate model. Juvenile macaques were immunized at 0, 3 and 6 weeks either with 10 or 5 µg doses of EV71vac formulated with AlPO4 adjuvant, or PBS as control. During the 56 weeks of studies, no fever nor local redness and swelling at sites of injections was observed in the immunized macaques. After single immunization, 100% seroconversion based on 4-fold increased in neutralization titer (Nt) was detected in EV71vac immunized monkeys but not PBS controls. A dose-dependent IgG antibody response was observed in monkeys receiving EV71vac immunization. The Nt of EV71vac immunized macaques had reached the peak after 3 vaccinations, then decreased gradually; however, the GMT of neutralizing antibody in the EV71vac immunized macaques were still above 100 at the end of the study. Correspondingly, both dose- and time-dependent interferon-γ and CD4+ T cell responses were detected in monkeys receiving EV71vac. Interestingly, similar to human responses, the dominant T cell epitopes of macaques were identified mainly in VP2 and VP3 regions. In addition, strong cross-neutralizing antibodies against most EV71 subgenotypes except some C2 and C4b strains, and Coxsackievirus A16 were observed. In summary, our results indicate that EV71vac elicits dose-dependent T-cell and antibody responses in macaques that could be a good animal model for evaluating the long-term immune responses elicited by EV71 vaccines.

Use of the 5' untranslated region and VP1 region to examine the molecular diversity in enterovirus B species

Human enteroviruses evolve quickly. The 5' untranslated region (UTR) is fundamentally important for efficient viral replication and for virulence; the VP1 region correlates well with antigenic typing by neutralization, and can be used for virus identification and evolutionary studies. In order to investigate the molecular diversity in EV-B species, the 5' UTR and VP1 regions were analysed for 208 clinical isolates from a single public-health laboratory (serving New South Wales, Australia), representing 28 EV-B types. Sequences were compared with the 5' UTR and VP1 regions of 98 strains available in GenBank, representing the same 28 types. The genetic relationships were analysed using two types of software (mega and BioNumerics). The sequence analyses of the 5' UTR and VP1 regions of 306 EV-B strains demonstrated that: (i) comparing the two regions gives strong evidence of epidemiological linkage of strains in some serotypes; (ii) the intraserotypic genetic variation within each gene reveals that they evolve distinctly largely due to their different functions; and (iii) mutation and possible recombination in the two regions play significant roles in the molecular diversity of EV-B. Understanding the tempo and pattern of molecular diversity and evolution is of great importance in the pathogenesis of EV-B enteroviruses, information which will assist in disease prevention and control.

Whole genomic characterization of a Korean human parechovirus type 1 (HPeV1) identifies recombination events

A human parechovirus (HPeV) strain CAU10-NN was detected from a stool sample of a 2-year-old healthy female infant in South Korea using a metagenomic approach. The CAU10-NN virus was isolated using a cell culture system and its whole genome was analyzed. The RNA genome of the CAU10-NN strain consists of 7,348 nucleotides (nt), excluding a poly(A) tail. A large open reading frame of 6,540 nt that encodes a putative polyprotein precursor of 2,180 amino acids is flanked by a 5'-untranslated region (UTR) of 708 nt and 3'-UTR of 88 nt followed by a poly(A) tail. Phylogenetic analysis shows that the CAU10-NN strain belongs to HPeV1. SimPlot and Bootscan analyses reveal that the virus genome is composed of regions related to corresponding genomic regions of other HPeVs. Recombination analysis indicates that the CAU10-NN strain might be a product of more than one genomic recombination event that occurred among HPeV1, HPeV3, and HPeV4 strains.

Complete genome characterization of a novel enterovirus type EV-B106 isolated in China, 2012

Human enterovirus B106 (EV-B106) is a recently identified member of enterovirus species B. In this study, we report the complete genomic characterization of an EV-B106 strain (148/YN/CHN/12) isolated from an acute flaccid paralysis patient in Yunnan Province, China. The new strain had 79.2–81.3% nucleotide and 89.1–94.8% amino acid similarity in the VP1 region with the other two EV-B106 strains from Bolivia and Pakistan. When compared with other EV serotypes, it had the highest (73.3%) VP1 nucleotide similarity with the EV-B77 prototype strain CF496-99. However, when aligned with all EV-B106 and EV-B77 sequences available from the GenBank database, two major frame shifts were observed in the VP1 coding region, which resulted in substantial (20.5%) VP1 amino acid divergence between the two serotypes. Phylogenetic analysis and similarity plot analysis revealed multiple recombination events in the genome of this strain. This is the first report of the complete genome of EV-B106.

Whole genomic sequence and replication kinetics of a new enterovirus C96 isolated from Guangdong, China with a different cell tropism

Enterovirus 96 (EV-C96) is a newly described serotype within the enterovirus C (EV-C) species, and its biological and pathological characters are largely unknown. In this study, we sequenced the whole genome of a novel EV-C96 strain that was isolated in 2011 from a patient with acute flaccid paralysis (AFP) in Guangdong province, China and characterized the properties of its infection. Sequence analysis revealed the close relationship between the EV-C96 strains isolated from the Guangdong and Shandong provinces of China, and suggested that recombination events occurred both between these EV-C96 strains and with other EV-C viruses. Moreover, the virus replication kinetics showed EV-C96 Guangdong strain replicated at a high rate in RD cells and presented a different cell tropism to other strains isolated from Shandong recently. These findings gave further insight into the evolutionary processes and extensive biodiversity of EV-C96.

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

We analysed natural recombination in 79 Human enterovirus A strains representing 13 serotypes by sequencing of VP1, 2C and 3D genome regions. The half-life of a non-recombinant tree node in coxsackieviruses 2, 4 and 10 was only 3.5 years, and never more than 9 years. All coxsackieviruses that differed by more than 7 % of the nucleotide sequence in any genome region were recombinants relative to each other. Enterovirus 71 (EV71), on the contrary, displayed remarkable genetic stability. Three major EV71 clades were stable for 19-29 years, with a half-life of non-recombinant viruses between 13 and 18.5 years in different clades. Only five EV71 strains out of over 150 recently acquired non-structural genome regions from coxsackieviruses, while none of 80 contemporary coxsackieviruses had non-structural genes transferred from the three EV71 clades. In contrast to earlier observations, recombination between VP1 and 2C genome regions was not more frequent than between 2C and 3D regions.

Enterovirus 74 infection in children

Enterovirus 74 (EV74) is a rarely detected viral infection of children. In 2010, EV74 was identified in New Zealand in a 2 year old child with acute flaccid paralysis (AFP) through routine polio AFP surveillance. A further three cases of EV74 were identified in children within six months. These cases are the first report of EV74 in New Zealand. In this study we describe the near complete genome sequence of four EV74 isolates from New Zealand, which shows only limited sequence identity in the non-structural proteins when compared to the other two known EV74 sequences. As is typical of enteroviruses multiple recombination events were evident, particularly in the P2 region and P3 regions. This is the first complete EV74 genome sequenced from a patient with acute flaccid paralysis.

Human enterovirus 71 epidemics: what's next?

Human enterovirus 71 (EV71) epidemics have affected various countries in the past 40 years. EV71 commonly causes hand, foot and mouth disease (HFMD) in children, but can result in neurological and cardiorespiratory complications in severe cases. Genotypic changes of EV71 have been observed in different places over time, with the emergence of novel genotypes or subgenotypes giving rise to serious outbreaks. Since the late 1990s, intra- and inter-typic recombination events in EV71 have been increasingly reported in the Asia-Pacific region. In particular, 'double-recombinant' EV71 strains belonging to a novel genotype D have been predominant in mainland China and Hong Kong over the last decade, though co-circulating with a minority of other EV71 subgenotypes and coxsackie A viruses. Continuous surveillance and genome studies are important to detect potential novel mutants or recombinants in the near future. Rapid and sensitive molecular detection of EV71 is of paramount importance in anticipating and combating EV71 outbreaks.

Analysis of the nonstructural and structural polyprotein regions, and complete genome sequences of Israel acute paralysis viruses identified from honeybees (Apis mellifera) in Korea

Phylogenetic trees were constructed for 24 partial nucleotide sequences of the nonstructural polyprotein (ORF1) and structural polyprotein regions (ORF2) of Korean IAPV genotypes, as well as eight previously reported IAPV sequences from various countries. Most of the Korean genotypes formed a distinct cluster, separate from other country genotypes. To investigate this phenomenon in more detail, three complete IAPV genome sequences were identified from different regions in Korea, i.e., Korea1, Korea2, and Korea3. These sequences were aligned with eight previously reported complete genome sequences and various genome regions were compared. The Korean IAPVs were very similar to those from China and Israel, but highly diverged from USA and Australian genotypes. Interestingly, they showed greater variability than the USA and Australian genotypes in ORF1, but highly similar to the Australian genotype in the ORF2 region. Thus, genetic recombination may account for the spatial distance between the Korean IAPV genotypes and those from other countries.

Genetic characterization of EV71 isolates from 2004 to 2010 reveals predominance and persistent circulation of the newly proposed genotype D and recent emergence of a distinct lineage of subgenotype C2 in Hong Kong

Background

Enterovirus 71 (EV71) is a common etiological agent of hand, foot and mouth disease (HFMD) in children. EV71 epidemics have been reported in Hong Kong in recent years, and yet the genetic information of EV71 strains circulating in our locality is limited. The objective of this study was to investigate the genetic evolution of these EV71 isolates in Hong Kong over a 7-year period.

Methods

Twenty-two EV71 isolates from Hong Kong during 2004–2010 were included for phylogenetic analysis using partial VP2-VP3, 2C and 3D regions. Eight EV71 strains were selected for complete genome sequencing and recombination analysis.

Results

Among the 22 EV71 isolates, 20 belonged to subgenotype C4 and 2 belonged to subgenotype C2 based on the phylogenetic analysis of partial VP2-VP3, 2C and 3D gene regions. Phylogenetic, similarity plot and bootscan analyses using complete genome sequences of seven EV71 isolates of subgenotype C4 supported that the “double-recombinant” strains of subgenotype C4 persistently circulating in Hong Kong should belong to a newly proposed genotype D. Further analysis revealed two clusters, subgenotypes C4b and C4a (proposed genotypes D1a and D1b respectively), with “genotype D1b” strains being predominant in recent years in Hong Kong. A distinct lineage of EV71 subgenotype C2 has emerged in Hong Kong in 2008. The evolutionary rate of EV71 was 3.1 × 10^-3 nucleotide substitutions per site per year similar to that of other enterovirus, such as EV68, but was relatively lower than those of echovirus 30 and poliovirus. Molecular clock analysis using VP1 gene dated the time to the most recent common ancestor of all EV71 genotypes to 1900s, while the EV71 “double-recombinant” strains of “genotype D” were detected as early as 1998.

Conclusions

This study provides the molecular basis for proposing a new “genotype D” of EV71 and assigning a discrete lineage of subgenotype C2. EV71 strains of “genotype D” have been circulating in Hong Kong for over 7 years, with “genotype D1b” being predominant.

Epidemics and Frequent Recombination within Species in Outbreaks of Human Enterovirus B-Associated Hand, Foot and Mouth Disease in Shandong China in 2010 and 2011

The epidemiology and molecular characteristics of human enterovirus B (HEV-B) associated with hand, foot and mouth disease (HFMD) outbreaks in China are not well known. In the present study, we tested 201 HEV isolates from 233 clinical specimens from patients with severe HFMD during 2010-2011 in Linyi, Shandong, China. Of the 201 isolates, 189 were fully typed and 18 corresponded to HEV-B species (six serotypes CVA9, CVB1, CVB4, Echo 6, Echo 25 and Echo 30) using sensitive semi-nested polymerase chain reaction analysis of VP1 gene sequences. Phylogenetic analysis based on the VP1 region showed that eight E30SD belonged to a novel sub-genogroup D2; E25SD belonged to a novel sub-genogroup D6; E6SD belonged to sub-lineage C6 and five CVB1SD belonged to subgroup 4C; and B4SD belonged sub-lineage D2. The full viral genomes of the CVB1SD, E6SD, E25SD and E30SD isolates were sequenced. Analysis of phylogenetic and similarity plots indicated that E25SD recombined with E25-HN-2, E30FDJS03 and E4AUS250 at noncontiguous P2A-P3D regions, while E30SD, E30FDJ03, E25-HN-2 and E9 DM had shared sequences in discrete regions of P2 and P3. Both E6SD and B1SD shared sequences with E1-HN, B4/GX/10, B5-HN, and A9-Alberta in contiguous regions of most of P2 and P3. Genetic algorithm recombination detection analysis further confirmed the existence of multiple potential recombination points. In conclusion, analysis of the complete genomes of E25SD, E30SD, CVB1SD and E6SD isolated from HFMD patients revealed that they formed novel subgenogroup. Given the prevalence and recombination of these viruses in outbreaks of HFMD, persistent surveillance of HFMD-associated HEV-B pathogens is required to predict potential emerging viruses and related disease outbreaks.

Specific primer amplification of the VP1 region directed by 5' UTR sequence analysis: enterovirus testing and identification in clinical samples from hand-foot-and-mouth disease patients

Many genotypes of the enterovirus (EV) pathogens can cause clinical hand-foot-and-mouth disease (HFMD). Therefore, rapid identification and monitoring of HFMD pathogens can be difficult, especially from the original clinical specimens. In this study, both universal pan-enterovirus and EV71/CA16 VP1-specific primer sets were designed and used to examine clinical specimens from HFMD patients. Based on the initial sequence analysis of the 5'-untanslated region (5'-UTR) and VP1 amplification products, additional primers for the VP1 region were redesigned for further genotyping of the remaining small portion non-EV71/non-CA16 specimens. With a known panel, it was possible to identify 15 out of 16 members using 5'-UTR sequence typing and VP1 typing, suggesting good detectability and genotyping of this method. One strain that was not typed by 5'-UTR was shown to be a recombinant virus. When this method was applied to examine clinical specimens from 44 suspected HFMD patients, 41 were detected as EV positive. In only one case, the VP1 sequence could not be identified. Four types of EVs, including CA16 (26/41, 63.4%), EV71-C4 (6/41, 14.6%), CA6 (5/41, 12.2%) and CA10 (3/41, 7.3%), were detected. In conclusion, 5' UTR amplification sequencing and subsequent VP1 specific primer amplification ensures a high detection rate and good genotyping accuracy in the examination of clinical samples. This detection strategy can be used for routine evaluation and monitoring of HFMD to follow local trends of EV infection.

Genetic characterization of a Coxsackie A9 virus associated with aseptic meningitis in Alberta, Canada in 2010

Background

An unusually high incidence of aseptic meningitis caused by enteroviruses was noted in Alberta, Canada between March and October 2010. Sequence based typing was performed on the enterovirus positive samples to gain a better understanding of the molecular characteristics of the Coxsackie A9 (CVA-9) strain responsible for most cases in this outbreak.

Methods

Molecular typing was performed by amplification and sequencing of the VP2 region. The genomic sequence of one of the 2010 outbreak isolates was compared to a CVA-9 isolate from 2003 and the prototype sequence to study genetic drift and recombination.

Results

Of the 4323 samples tested, 213 were positive for enteroviruses (4.93%). The majority of the positives were detected in CSF samples (n = 157, 73.71%) and 81.94% of the sequenced isolates were typed as CVA-9. The sequenced CVA-9 positives were predominantly (94.16%) detected in patients ranging in age from 15 to 29 years and the peak months for detection were between March and October. Full genome sequence comparisons revealed that the CVA-9 viruses isolated in Alberta in 2003 and 2010 were highly homologous to the prototype CVA-9 in the structural VP1, VP2 and VP3 regions but divergent in the VP4, non-structural and non-coding regions.

Conclusion

The increase in cases of aseptic meningitis was associated with enterovirus CVA-9. Sequence divergence between the prototype strain of CVA-9 and the Alberta isolates suggests genetic drifting and/or recombination events, however the sequence was conserved in the antigenic regions determined by the VP1, VP2 and VP3 genes. These results suggest that the increase in CVA-9 cases likely did not result from the emergence of a radically different immune escape mutant.