Characterization of a putative ancestor of coxsackievirus B5

Influence of Amino Acid Substitutions in Capsid Proteins of Coxsackievirus B5 on Free Chlorine and Thermal Inactivation

The sensitivity of enteroviruses to disinfectants varies among genetically similar variants and coincides with amino acid changes in capsid proteins, although the effect of individual substitutions remains unknown. Here, we employed reverse genetics to investigate how amino acid substitutions in coxsackievirus B5 (CVB5) capsid proteins affect the virus' sensitivity to free chlorine and heat treatment. Of ten amino acid changes observed in CVB5 variants with free chlorine resistance, none significantly reduced the chlorine sensitivity, indicating a minor role of the capsid composition in chlorine sensitivity of CVB5. Conversely, a subset of these amino acid changes located at the C-terminal region of viral protein 1 led to reduced heat sensitivity. Cryo-electron microscopy revealed that these changes affect the assembly of intermediate viral states (altered and empty particles), suggesting that the mechanism for reduced heat sensitivity could be related to improved molecular packing of CVB5, resulting in greater stability or altered dynamics of virus uncoating during infection.

Molecular characterization of coxsackievirus B5 from the sputum of pneumonia children patients of Kunming, Southwest China

BACKGROUND

CVB5 can cause respiratory infections. However, the molecular epidemiological information about CVB5 in respiratory tract samples is still limited. Here, we report five cases in which CVB5 was detected in sputum sample of pneumonia children patients from Kunming, Southwest China.

METHODS

CVB5 isolates were obtained from sputum samples of patients with pneumonia. Whole-genome sequencing of CVB5 isolates was performed using segmented PCR, and phylogenetic, mutation and recombination analysis. The effect of mutations in the VP1 protein on hydration were analyzed by Protscale. The tertiary models of VP1 proteins were established by Colabfold, and the effect of mutations in VP1 protein on volume modifications and binding affinity were analyzed by Pymol software and PROVEAN.

RESULTS

A total of five CVB5 complete genome sequences were obtained. No obvious homologous recombination signals comparing with other coxsackie B viruses were observed in the five isolates. Phylogenetic analysis showed that the five CVB5 sputum isolates were from an independent branch in genogroup E. Due to the mutation, the structure and spatial of the VP1 protein N-terminus have changed significantly. Comparing to the Faulkner (CVB5 prototype strain), PROVEAN revealed three deleterious substitutions: Y75F, N166T (KM35), T140I (KM41). The last two of the three deleterious substitutions significantly increased the hydrophobicity of the residues.

CONCLUSIONS

We unexpectedly found five cases of CVB5 infection instead of rhinoviruses infection during our routine surveillance of rhinoviruses in respiratory tract samples. All five patients were hospitalized with pneumonia symptoms and were not tested for enterovirus during their hospitalization. This report suggests that enterovirus surveillance in patients with respiratory symptoms should be strengthened.

Genotype-dependent kinetics of enterovirus inactivation by free chlorine and ultraviolet (UV) irradiation

Inactivation kinetics of enterovirus by disinfection is often studied using a single laboratory strain of a given genotype. Environmental variants of enterovirus are genetically distinct from the corresponding laboratory strain, yet it is poorly understood how these genetic differences affect inactivation. Here we evaluated the inactivation kinetics of nine coxsackievirus B3 (CVB3), ten coxsackievirus B4 (CVB4), and two echovirus 11 (E11) variants by free chlorine and ultraviolet irradiation (UV). The inactivation kinetics by free chlorine were genotype- (i.e., susceptibility: CVB5 < CVB3 ≈ CVB4 < E11) and genogroup-dependent and exhibited up to 15-fold difference among the tested viruses. In contrast, only minor (up to 1.3-fold) differences were observed in the UV inactivation kinetics. The differences in variability between the two disinfectants could be rationalized by their respective inactivation mechanisms: inactivation by UV mainly depends on the genomic size and composition, which was similar for all viruses tested, whereas free chlorine targets the viral capsid protein, which exhibited critical differences between genogroups and genotypes. Finally, we integrated the observed variability in inactivation rate constants into an expanded Chick-Watson model to estimate the overall inactivation of an enterovirus consortium. The results highlight that the distribution of inactivation rate constants and the abundance of each genotype are essential parameters to accurately predict the overall inactivation of an enterovirus population by free chlorine. We conclude that predictions based on inactivation data of a single variant or reference pathogen alone likely overestimate the true disinfection efficiency of free chlorine.

Atomic Structures of Coxsackievirus B5 Provide Key Information on Viral Evolution and Survival

Coxsackie virus B5 (CVB5), a main serotype in human Enterovirus B (EVB), can cause severe viral encephalitis and aseptic meningitis among infants and children. Currently, there is no approved vaccine or antiviral therapy available against CVB5 infection. Here, we determined the atomic structures of CVB5 in three forms: mature full (F) particle (2.73 Å), intermediate altered (A) particle (2.81 Å), and procapsid empty (E) particle (2.95 Å). Structural analysis of F particle of CVB5 unveiled similar structures of “canyon,” “puff,” and “knob” as those other EV-Bs. We observed structural rearrangements that are alike during the transition from F to A particle, indicative of similar antigenicity, cell entry, and uncoating mechanisms shared by all EV-Bs. Further comparison of structures and sequences among all structure-known EV-Bs revealed that while the residues targeted by neutralizing MAbs are diversified and drive the evolution of EV-Bs, the relative conserved residues recognized by uncoating receptors could serve as the basis for the development of antiviral vaccines and therapeutics.

IMPORTANCE As one of the main serotypes in Enterovirus B, CVB5 has been commonly reported in recent years. The atomic structures of CVB5 shown here revealed classical features found in EV-Bs and the structural rearrangement occurring during particle expansion and uncoating. Also, structure- and sequence-based comparison between CVB5 and other structure-known EV-Bs screened out key domains important for viral evolution and survival. All these provide insights into the development of vaccine and therapeutics for EV-Bs.

Immunological Study of Reconstructed Common Ancestral Sequence of Adenovirus Hexon Protein

Aim: To reconstruct the ancestral sequence of human adenoviral hexon protein by combining sequence variations and structural information. And to provide a candidate hexon protein for developing new adenoviral vector capable of escaping the pre-existing immunity in healthy populations.

Methods: The sequences of 74 adenovirus-type strains were used to predict the ancestral sequence of human adenovirus hexon protein using FastML and MEGA software. The three-dimensional structure model was built using homology modeling methods. The immunological features of ancestral loop 1 and loop 2 regions of sequences were tested using protein segments expressed in a prokaryotic expression system and polypeptides synthesized with human serum samples.

Results: The tower region of the hexon protein had the highest sequence variability, while the neck and base regions remained constant among different types. The modern strains successfully predicted the common ancestral sequence of the human adenovirus hexon. The positive sera against neutralizing epitopes on the common ancestor of adenoviral hexon were relatively rare among healthy adults.

Conclusion: The existing strains inferred the common ancestor of human adenoviruses, with epitopes never observed in the current human strains. The predicted common ancestor hexon is a good prospect in the improvement of adenovirus vectors.

Coxsackievirus B5 virus-like particle vaccine exhibits greater immunogenicity and immunoprotection than its inactivated counterpart in mice

Coxsackievirus B group 5 (CVB5) represents one of the major pathogens that cause diseases such as hand, foot and mouth disease (HFMD) and aseptic meningitis et al. Currently, no specific drugs and vaccines are available, and a safe and effective CVB5 vaccine is of great value for control of the diseases. In this study, CVB5 P1 precursor and 3CD protease were co-expressed in Sf9 cells by using a baculovirus expression system. The P1 was processed by 3CD and self-assembled into CVB5 virus-like particles (VLPs). VP1 and VP3 capsid proteins of CVB5 could be detected by SDS-PAGE and Western blotting. Transmission electron microscopy revealed that the CVB5 VLPs were spherical particles with a diameter of about 30 nm, mimicking wild-type CVB5 virus. Our study showed that the total IgG and neutralizing antibodies induced by CVB5 VLPs were higher than those induced by inactivated vaccine. More importantly, the CVB5 VLPs conferred full protection to the CVB5-challenged suckling mice via passive immunity while protection efficiency of the inactivated vaccine was only 80%. The CVB5 VLPs vaccine could protect the limb muscles, brain, and heart tissues of suckling mice from CVB5-induced damage. These results demonstrated that the CVB5 VLPs vaccine possessed stronger immunogenicity and provided more robust immunoprotection than the inactivated CVB5 vaccine, suggesting that the CVB5 VLPs promise to be a CVB5 vaccine candidate in future.

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.

Slow Infection due to Lowering the Amount of Intact versus Empty Particles Is a Characteristic Feature of Coxsackievirus B5 Dictated by the Structural Proteins

Enterovirus B species typically cause a rapid cytolytic infection leading to efficient release of progeny viruses. However, they are also capable of persistent infections in tissues, which are suggested to contribute to severe chronic states such as myocardial inflammation and type 1 diabetes. In order to understand the factors contributing to differential infection strategies, we constructed a chimera by combining the capsid proteins from fast-cytolysis-causing echovirus 1 (EV1) with nonstructural proteins from coxsackievirus B5 (CVB5), which shows persistent infection in RD cells. The results showed that the chimera behaved similarly to parental EV1, leading to efficient cytolysis in both permissive A549 and semipermissive RD cells. In contrast to EV1 and the chimera, CVB5 replicated slowly in permissive cells and showed persistent infection in semipermissive cells. However, there was no difference in the efficiency of uptake of CVB5 in A549 or RD cells in comparison to the chimera or EV1. CVB5 batches constantly contained significant amounts of empty capsids, also in comparison to CVB5's close relative CVB3. During successive passaging of batches containing only intact CVB5, increasing amounts of empty and decreasing amounts of infective capsids were produced. Our results demonstrate that the increase in the amount of empty particles and the lowering of the amount of infective particles are dictated by the CVB5 structural proteins, leading to slowing down of the infection between passages. Furthermore, the key factor for persistent infection is the small amount of infective particles produced, not the high number of empty particles that accumulate.IMPORTANCE Enteroviruses cause several severe diseases, with lytic infections that lead to rapid cell death but also persistent infections that are more silent and lead to chronic states of infection. Our study compared a cytolytic echovirus 1 infection to persistent coxsackievirus B5 infection by making a chimera with the structural proteins of echovirus 1 and the nonstructural proteins of coxsackievirus B5. Coxsackievirus B5 infection was found to lead to the production of a high number of empty viruses (empty capsids) that do not contain genetic material and are unable to continue the infection. Coinciding with the high number of empty capsids, the amount of infective virions decreased. This characteristic property was not observed in the constructed chimera virus, suggesting that structural proteins are in charge of these phenomena. These results shed light on the mechanisms that may cause persistent infections. Understanding events leading to efficient or inefficient infections is essential in understanding virus-caused pathologies.

Alignment Modulates Ancestral Sequence Reconstruction Accuracy

Accurate reconstruction of ancestral states is a critical evolutionary analysis when studying ancient proteins and comparing biochemical properties between parental or extinct species and their extant relatives. It relies on multiple sequence alignment (MSA) which may introduce biases, and it remains unknown how MSA methodological approaches impact ancestral sequence reconstruction (ASR). Here, we investigate how MSA methodology modulates ASR using a simulation study of various evolutionary scenarios. We evaluate the accuracy of ancestral protein sequence reconstruction for simulated data and compare reconstruction outcomes using different alignment methods. Our results reveal biases introduced not only by aligner algorithms and assumptions, but also tree topology and the rate of insertions and deletions. Under many conditions we find no substantial differences between the MSAs. However, increasing the difficulty for the aligners can significantly impact ASR. The MAFFT consistency aligners and PRANK variants exhibit the best performance, whereas FSA displays limited performance. We also discover a bias towards reconstructed sequences longer than the true ancestors, deriving from a preference for inferring insertions, in almost all MSA methodological approaches. In addition, we find measures of MSA quality generally correlate highly with reconstruction accuracy. Thus, we show MSA methodological differences can affect the quality of reconstructions and propose MSA methods should be selected with care to accurately determine ancestral states with confidence.

Genetic diversity of the enteroviruses detected from cerebrospinal fluid (CSF) samples of patients with suspected aseptic meningitis in northern West Bank, Palestine in 2017

BACKGROUND

Human enterovirus genus showed a wide range of genetic diversity.

OBJECTIVES

To investigate the genetic diversity of the enteroviruses isolated in 2017 in northern West Bank, Palestine.

STUDY DESIGN

249 CSF samples from aseptic meningitis cases were investigated for HEV using two RT-PCR protocols targeting the 5' NCR and the VP1 region of the HEV genome. The phylogenetic characterization of the sequenced VP1 region of Echovirus18 (E18) and Coxsackievirus B5 (CVB5) isolated in Palestine along with 27 E18 and 27 CVB5 sequences available from the Genbank were described.

RESULTS

E18 and CVB5 account for 50% and 35% of the successfully HEV types, respectively. Phylogenetic tree of E18 and CVB5 showed three main clusters, with all Palestinian isolates uniquely clustering together with those from China and from different countries, respectively. Cluster I of E18, with 13 Palestinian and 6 Chinese isolates, showed the lowest haplotype-to-sequence ratio (0.6:1), haplotype diversity (Hd), nucleotide diversity (π), and number of segregating sites (S) compared to clusters II and III. Furthermore, cluster I showed negative Tajima's D and Fu-Li'sF tests with statistically significant departure from neutrality (P<0.01). In both E18 and CVB5 populations, high haplotype diversity, but low genetic diversity was evident. Inter-population pairwise genetic distance (Fst) and gene flow (Nm) showed that the Palestinian E18 and CVB5 clusters were highly differentiated from the other clusters.

CONCLUSIONS

The study divulged close genetic relationship between Palestinian HEV strains as confirmed by population genetics and phylogenetic analyses.

A neonatal mouse model of central nervous system infections caused by Coxsackievirus B5

As one of the key members of the coxsackievirus B group, coxsackievirus B5 (CV-B5) can cause many central nervous system diseases, such as viral encephalitis, aseptic meningitis, and acute flaccid paralysis. Notably, epidemiological data indicate that outbreaks of CV-B5-associated central nervous system (CNS) diseases have been reported worldwide throughout history. In this study, which was conducted to promote CV-B5 vaccine and anti-virus drug research, a 3-day-old BALB/c mouse model was established using a CV-B5 clinical isolate (CV-B5/JS417) as the challenge strain. Mice challenged with CV-B5/JS417 exhibited a series of neural clinical symptoms and death with necrosis of neuronal cells in the cerebral cortex and the entire spinal cord, hindlimb muscles, and cardiomyocytes. The viral load of each tissue at various post-challenge time points suggested that CV-B5 replicated in the small intestine and was subsequently transmitted to various organs via viremia; the virus potentially entered the brain through the spinal axons, causing neuronal cell necrosis. In addition, this mouse model was used to evaluate the protective effect of a CV-B5 vaccine. The results indicated that both the inactivated CV-B5 vaccine and anti-CVB5 serum significantly protected mice from a lethal infection of CV-B5/JS417 by producing neutralizing antibodies. In summary, the first CV-B5 neonatal mouse model has been established and can sustain CNS infections in a manner similar to that observed in humans. This model will be a useful tool for studies on pathogenesis, vaccines, and anti-viral drug evaluations.

Evolutionary histories of coxsackievirus B5 and swine vesicular disease virus reconstructed by phylodynamic and sequence variation analyses

Coxsackievirus (CV)-B5 is a common human enterovirus reported worldwide; swine vesicular disease virus (SVDV) is a porcine variant of CV-B5. To clarify the transmission dynamics and molecular basis of host switching between CV-B5 and SVDV, we analysed and compared the VP1 and partial 3D^pol gene regions of these two viruses. Spatiotemporal dynamics of viral transmission were estimated using a Bayesian statistical inference framework. The detected selection events were used to analyse the key molecules associated with host switching. Analyses of VP1 sequences revealed six CV-B5 genotypes (A1–A4 and B1–B2) and three SVDV genotypes (I–III). Analyses of partial 3D^pol revealed five clusters (A–E). The genotypes evolved sequentially over different periods, albeit with some overlap. The major hub of CV-B5 transmission was in China whereas the major hubs of SVDV transmission were in Italy. Network analysis based on deduced amino acid sequences showed a diverse extension of the VP1 structural protein, whereas most sequences were clustered into two haplotypes in the partial 3D^pol region. Residue 178 of VP1 showed four epistatic interactions with residues known to play essential roles in viral host tropism, cell entry, and viral decoating.

Seroepidemiology of coxsackievirus B5 in infants and children in Jiangsu province, China

Coxsackievirus B5 (CV-B5) is associated with various human diseases such as viral encephalitis, aseptic meningitis, paralysis, herpangina, and hand, foot and mouth disease (HFMD). However, there is currently no effective vaccine against CV-B5.The seroepidemiologic characteristics of CV-B5 remained unknown. A cohort study was carried out in 176 participants aged 6-35 months from January 2012 to January 2014. The serum samples were collected and tested for CV-B5 neutralizing antibodies (NtAbs) four times during these two years. The confirmed enterovirus cases were recorded through the surveillance system, and their throat or rectal swabs were collected for pathogen detection. According to the changes of CV-B5 NtAbs, two CV-B5 epidemics were detected among these participants during the two-year follow-up. Sixty-seven cases out of all participants had seroconversion in CV-B5 NtAbs. During the first epidemic from March 2012 to September 2012, CV-B5 seropositivity rate increased significantly (6.8%, 12/176 vs. 21.6%, 38/176, P = 0.000). The seroconversion rate and geometric mean fold-increase (GMFI) were 18.2% (32/176) and 55.7, respectively; During the second epidemic from September 2012 to January 2014, CV-B5 seropositivity rate also increased (21.6%, 38/176 vs. 38.6%, 68/176, P = 0.000), and the seroconversion rate and GMFI were 19.9% (35/176) and 46.5, respectively. Only one case had CV-B5 associated HFMD during the two-year follow-up, and CV-B5 from the throat swab isolate was GI.D3 subtype, which belonged to the major pandemic strain in mainland China. CV-B5 infection was common in infants and children in Jiangsu province, China. Therefore, it's necessary to strengthen the surveillance on CV-B5 and to understand the epidemic characteristics of CV-B5 infection.

Analysis of an echovirus 18 outbreak in Thuringia, Germany: insights into the molecular epidemiology and evolution of several enterovirus species B members

In October and November 2010, six children and one woman were presented with symptoms of aseptic meningitis in Jena, Thuringia, Germany. Enterovirus RNA was detected in the cerebrospinal fluid of all patients by RT-PCR, and preliminary molecular typing revealed echovirus 18 (E-18) as causative agent. Virus isolates were obtained from stool samples of three patients and several contact persons. Again, most isolates were typed as E-18. In addition, coxsackievirus B5 (CV-B5) and echovirus 25 (E-25) were found to co-circulate. As only few complete E-18 sequences are available in GenBank, the entire genomes of these isolates were determined using direct RNA-sequencing technology. We did not find evidence for recombination between E-18, E-25 or CV-B5 during the outbreak. Viral protein 1 gene sequences and the cognate 3D polymerase gene sequences of each isolate and GenBank sequences were analysed in order to define type-specific recombination groups (recogroups).

AAV ancestral reconstruction library enables selection of broadly infectious viral variants

Adeno-associated virus (AAV) vectors have achieved clinical efficacy in treating several diseases. However, enhanced vectors are required to extend these landmark successes to other indications and protein engineering approaches may provide the necessary vector improvements to address such unmet medical needs. To generate new capsid variants with potentially enhanced infectious properties and to gain insights into AAV's evolutionary history, we computationally designed and experimentally constructed a putative ancestral AAV library. Combinatorial variations at 32 amino acid sites were introduced to account for uncertainty in their identities. We then analyzed the evolutionary flexibility of these residues, the majority of which have not been previously studied, by subjecting the library to iterative selection on a representative cell line panel. The resulting variants exhibited transduction efficiencies comparable to the most efficient extant serotypes and, in general, ancestral libraries were broadly infectious across the cell line panel, indicating that they favored promiscuity over specificity. Interestingly, putative ancestral AAVs were more thermostable than modern serotypes and did not use sialic acids, galactose or heparan sulfate proteoglycans for cellular entry. Finally, variants mediated 19- to 31-fold higher gene expression in the muscle compared with AAV1, a clinically used serotype for muscle delivery, highlighting their promise for gene therapy.

Viral meningitis epidemics and a single, recent, recombinant and anthroponotic origin of swine vesicular disease virus

Swine vesicular disease, an important viral disease affecting domestic pigs, is shown to have a single and recent origin in humans, leading us closer to a full understanding of the sudden emergence of this enigmatic veterinary disease, and exemplifying the sometimes overlooked risk of human to animal disease transfers.

Background and objectives: Swine vesicular disease virus (SVDV) is a close relative of the human Enterovirus B serotype, coxsackievirus B5. As the etiological agent of a significant emergent veterinary disease, several studies have attempted to explain its origin. However, several key questions remain, including the full biological ancestry of the virus, and its geographical and temporal origin.

Methodology: We sequenced near-complete genomes of 27 SVDV and 13 coxsackievirus B5 samples, all originally isolated between 1966 and 2006, and analysed these in conjunction with existing sequences and historical information.

Results: While analyses incorporating 24 additional near-complete SVDV genomic sequences indicate clear signs of within-SVDV recombination, all 51 SVDV isolates remain monophyletic. This supports a hypothesis of a single anthroponotic transfer origin. Analysis of individual coding and non-coding regions supports that SVDV has a recombinant origin between coxsackievirus B5 and another Enterovirus B serotype, most likely coxsackievirus A9. Extensive Bayesian sequence-based analysis of the time of the most recent common ancestor of all analysed sequences places this within a few years around 1961. Epidemiological evidence points to China as an origin, but there are no available samples to test this conclusively.

Conclusions and implications: Historical investigation and the clinical aspects of the involved Enterovirus B serotypes, makes the current results consistent with a hypothesis stating that SVDV originated through co-infection, recombination, and a single anthroponotic event, during large viral meningitis epidemics around 1960/1961 involving the ancestral serotypes. The exact geographical origin of SVDV may remain untestable due to historical aspects.

In Silico Reconstruction of the Viral Evolutionary Lineage Yields a Potent Gene Therapy Vector

Adeno-associated virus (AAV) vectors have emerged as a gene-delivery platform with demonstrated safety and efficacy in a handful of clinical trials for monogenic disorders. However, limitations of the current generation vectors often prevent broader application of AAV gene therapy. Efforts to engineer AAV vectors have been hampered by a limited understanding of the structure-function relationship of the complex multimeric icosahedral architecture of the particle. To develop additional reagents pertinent to further our insight into AAVs, we inferred evolutionary intermediates of the viral capsid using ancestral sequence reconstruction. In-silico-derived sequences were synthesized de novo and characterized for biological properties relevant to clinical applications. This effort led to the generation of nine functional putative ancestral AAVs and the identification of Anc80, the predicted ancestor of the widely studied AAV serotypes 1, 2, 8, and 9, as a highly potent in vivo gene therapy vector for targeting liver, muscle, and retina.

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.

The evolution of Vp1 gene in enterovirus C species sub-group that contains types CVA-21, CVA-24, EV-C95, EV-C96 and EV-C99

Genus Enterovirus (Family Picornaviridae,) consists of twelve species divided into genetically diverse types by their capsid protein VP1 coding sequences. Each enterovirus type can further be divided into intra-typic sub-clusters (genotypes). The aim of this study was to elucidate what leads to the emergence of novel enterovirus clades (types and genotypes). An evolutionary analysis was conducted for a sub-group of Enterovirus C species that contains types Coxsackievirus A21 (CVA-21), CVA-24, Enterovirus C95 (EV-C95), EV-C96 and EV-C99. VP1 gene datasets were collected and analysed to infer the phylogeny, rate of evolution, nucleotide and amino acid substitution patterns and signs of selection. In VP1 coding gene, high intra-typic sequence diversities and robust grouping into distinct genotypes within each type were detected. Within each type the majority of nucleotide substitutions were synonymous and the non-synonymous substitutions tended to cluster in distinct highly polymorphic sites. Signs of positive selection were detected in some of these highly polymorphic sites, while strong negative selection was indicated in most of the codons. Despite robust clustering to intra-typic genotypes, only few genotype-specific ‘signature’ amino acids were detected. In contrast, when different enterovirus types were compared, there was a clear tendency towards fixation of type-specific ‘signature’ amino acids. The results suggest that permanent fixation of type-specific amino acids is a hallmark associated with evolution of different enterovirus types, whereas neutral evolution and/or (frequency-dependent) positive selection in few highly polymorphic amino acid sites are the dominant forms of evolution when strains within an enterovirus type are compared.

In silico resurrection of the major vault protein suggests it is ancestral in modern eukaryotes

Vaults are very large oligomeric ribonucleoproteins conserved among a variety of species. The rat vault 3D structure shows an ovoid oligomeric particle, consisting of 78 major vault protein monomers, each of approximately 861 amino acids. Vaults are probably the largest ribonucleoprotein structures in eukaryote cells, being approximately 70 nm in length with a diameter of 40 nm—the size of three ribosomes and with a lumen capacity of 50 million ų. We use both protein sequences and inferred ancestral sequences for in silico virtual resurrection of tertiary and quaternary structures to search for vaults in a wide variety of eukaryotes. We find that the vault’s phylogenetic distribution is widespread in eukaryotes, but is apparently absent in some notable model organisms. Our conclusion from the distribution of vaults is that they were present in the last eukaryote common ancestor but they have apparently been lost from a number of groups including fungi, insects, and probably plants. Our approach of inferring ancestral 3D and quaternary structures is expected to be useful generally.

Phylogenetic patterns of human coxsackievirus B5 arise from population dynamics between two genogroups and reveal evolutionary factors of molecular adaptation and transmission

The aim of this study was to gain insights into the tempo and mode of the evolutionary processes that sustain genetic diversity in coxsackievirus B5 (CVB5) and into the interplay with virus transmission. We estimated phylodynamic patterns with a large sample of virus strains collected in Europe by Bayesian statistical methods, reconstructed the ancestral states of genealogical nodes, and tested for selection. The genealogies estimated with the structural one-dimensional gene encoding the VP1 protein and nonstructural 3CD locus allowed the precise description of lineages over time and cocirculating virus populations within the two CVB5 clades, genogroups A and B. Strong negative selection shaped the evolution of both loci, but compelling phylogenetic data suggested that immune selection pressure resulted in the emergence of the two genogroups with opposed evolutionary pathways. The genogroups also differed in the temporal occurrence of the amino acid changes. The virus strains of genogroup A were characterized by sequential acquisition of nonsynonymous changes in residues exposed at the virus 5-fold axis. The genogroup B viruses were marked by selection of three changes in a different domain (VP1 C terminus) during its early emergence. These external changes resulted in a selective sweep, which was followed by an evolutionary stasis that is still ongoing after 50 years. The inferred population history of CVB5 showed an alternation of the prevailing genogroup during meningitis epidemics across Europe and is interpreted to be a consequence of partial cross-immunity.

Emergence, circulation, and spatiotemporal phylogenetic analysis of coxsackievirus a6- and coxsackievirus a10-associated hand, foot, and mouth disease infections from 2008 to 2012 in Shenzhen, China

Sporadic hand, foot, and mouth disease (HFMD) outbreaks and other infectious diseases in recent years have frequently been associated with certain human enterovirus (HEV) serotypes. This study explored the prevalences and genetic characteristics of non-HEV71 and non-coxsackievirus A16 (CV-A16) human enterovirus-associated HFMD infections in Shenzhen, China. A total of 2,411 clinical stool specimens were collected from hospital-based surveillance for HFMD from 2008 to 2012. The detection of HEV was performed by real-time reverse transcription-PCR (RT-PCR) and RT-seminested PCR, and spatiotemporal phylogenetic analysis was performed based on the VP1 genes. A total of 1,803 (74.8%) strains comprising 28 different serotypes were detected. In the past 5 years, the predominant serotypes were HEV71 (60.0%), followed by CV-A16 (21.2%) and two uncommon serotypes, CV-A6 (13.0%) and CV-A10 (3.3%). However, CV-A6 replaced CV-A16 as the second most common serotype between 2010 and 2012. As an emerging pathogen, CV-A6 became as common a causative agent of HFMD as HEV71 in Shenzhen in 2012. Phylogenetic analysis revealed that little variation occurred in the Chinese HEV71 and CV-A16 strains. The genetic characteristics of the Chinese CV-A6 and CV-A10 strains displayed geographic differences. The CV-A6 and CV-A10 strains circulating in Shenzhen likely originated in Europe. It was found that human enteroviruses have a high mutation rate due to evolutionary pressure and frequent recombination (3.2 × 10(-3) to 6.4 ×10(-3) substitutions per site per year for HEV71, CV-A6, CV-A16, and CV-A10). Since certain serotypes are potential threats to the public health, this study provides further insights into the significance of the epidemiological surveillance of HFMD.

Full genome sequence of a novel coxsackievirus B5 strain isolated from neurological hand, foot, and mouth disease patients in China

Coxsackievirus B5 (CVB5) belongs to the human enterovirus B species within the family Picornaviridae. We report the complete genome sequence of a novel CVB5 strain, CVB5/SD/09, that is associated with neurological hand, foot, and mouth disease in China. The complete genome consists of 7,399 nucleotides, excluding the 3' poly(A) tail, and has an open reading frame that maps between nucleotide positions 744 and 7301 and encodes a 2,185-amino-acid polyprotein. Phylogenetic analysis based on different genome region regions reveals that CVB5/SD/09 belongs to a novel CVB5 lineage, and similarity plotting and bootscanning analysis based on the whole genome of CVB5 in the present study and those available in GenBank indicate that the genome of CVB5/SD/09 has a mosaic-like structure, suggesting that recombination between different CVB5 strains may occur.

Update on molecular characterization of coxsackievirus B5 strains

Among Coxsackie B viruses, Coxsckievirus B5 is one of the most predominant serotypes in human, it is frequently associated with cases of neurological diseases, epidemics of meningitis and is a common cause of cardiomyopathy and diabetes. In the present study 27 isolates of Coxsackievirus B5 from North Africa, obtained from cerebrospinal fluid and stool samples of healthy individuals, patients with acute flaccid paralysis or aseptic meningitis were investigated by partial sequencing in the 5' half of the VP1 region and compared to the up-to-date published Coxsackievirus B5 sequences in the same genomic region. Four distinct genomic groups and ten different clusters were individualized. Most of the isolates from Algeria and Tunisia belonged to two clusters. For both, the sequences from North Africa clustered mainly with sequences from European countries, the majority isolated recently during the 2000s. The analysis of the alignment of amino-acids sequences in the VP1 gene revealed four major substitutions in strains from different clusters, we also noticed changes in the BC-loop region; this region is associated with viral antigenicity. This study permit to better identify circulating Coxsackievirus B5 strains throughout the world and their genetic relationship. The protein analysis showed changes that could imply some antigenic significance. J. Med. Virol. 83:1247-1254, 2011. © 2011 Wiley-Liss, Inc.

Recombination in eukaryotic single stranded DNA viruses

Although single stranded (ss) DNA viruses that infect humans and their domesticated animals do not generally cause major diseases, the arthropod borne ssDNA viruses of plants do, and as a result seriously constrain food production in most temperate regions of the world. Besides the well known plant and animal-infecting ssDNA viruses, it has recently become apparent through metagenomic surveys of ssDNA molecules that there also exist large numbers of other diverse ssDNA viruses within almost all terrestrial and aquatic environments. The host ranges of these viruses probably span the tree of life and they are likely to be important components of global ecosystems. Various lines of evidence suggest that a pivotal evolutionary process during the generation of this global ssDNA virus diversity has probably been genetic recombination. High rates of homologous recombination, non-homologous recombination and genome component reassortment are known to occur within and between various different ssDNA virus species and we look here at the various roles that these different types of recombination may play, both in the day-to-day biology, and in the longer term evolution, of these viruses. We specifically focus on the ecological, biochemical and selective factors underlying patterns of genetic exchange detectable amongst the ssDNA viruses and discuss how these should all be considered when assessing the adaptive value of recombination during ssDNA virus evolution.

Resurrection of an ancestral 5S rRNA

Background

In addition to providing phylogenetic relationships, tree making procedures such as parsimony and maximum likelihood can make specific predictions of actual historical sequences. Resurrection of such sequences can be used to understand early events in evolution. In the case of RNA, the nature of parsimony is such that when applied to multiple RNA sequences it typically predicts ancestral sequences that satisfy the base pairing constraints associated with secondary structure. The case for such sequences being actual ancestors is greatly improved, if they can be shown to be biologically functional.

Results

A unique common ancestral sequence of 28 Vibrio 5S ribosomal RNA sequences predicted by parsimony was resurrected and found to be functional in the context of the E. coli cellular environment. The functionality of various point variants and intermediates that were constructed as part of the resurrection were examined in detail. When separately introduced the changes at single stranded positions and individual double variants at base-paired positions were also viable. An additional double variant was examined at a different base-paired position and it was also valid.

Conclusions

The results show that at least in the case of the 5S rRNAs considered here, ancestors predicted by parsimony are likely to be realistic when the prediction is not overly influenced by single outliers. It is especially noteworthy that the phenotype of the predicted ancestors could be anticipated as a cumulative consequence of the phenotypes of the individual variants that comprised them. Thus, point mutation data is potentially useful in evaluating the reasonableness of ancestral sequences predicted by parsimony or other methods. The results also suggest that in the absence of significant tertiary structure constraints double variants that preserve pairing in stem regions will typically be accepted. Overall, the results suggest that it will be feasible to resurrect additional meaningful 5S rRNA ancestors as well as ancestral sequences of many different types of RNA.

Genus-specific substitution rate variability among picornaviruses

Picornaviruses have some of the highest nucleotide substitution rates among viruses, but there have been no comparisons of evolutionary rates within this broad family. We combined our own Bayesian coalescent analyses of VP1 regions from four picornaviruses with 22 published VP1 rates to produce the first within-family meta-analysis of viral evolutionary rates. Similarly, we compared our rate estimates for the RNA polymerase 3D(pol) gene from five viruses to four published 3D(pol) rates. Both a structural and a nonstructural gene show that enteroviruses are evolving, on average, a half order of magnitude faster than members of other genera within the Picornaviridae family.

Evolution of newly described enteroviruses

Several new enterovirus serotypes and a new human rhinovirus species have been characterized in the Enterovirus genus recently, raising a question about the origin of the new viruses. In this article we attempt to outline the general patterns of enterovirus evolution, ultimately leading to the emergence of new serotypes or species. Different evolutionary and epidemiological patterns can be deduced between different enterovirus species, between entero- and rhino-viruses and between different serotypes within a species. This article presents a hypothesis that the divergent evolution leading to a new serotype is likely to involve adaptation to a new ecological niche either within a single host species or due to interspecies transmission. By contrast, evolution within a serotype appears to occur primarily by genetic drift.

A large Finnish echovirus 30 outbreak was preceded by silent circulation of the same genotype

An outbreak of echovirus 30 (E-30) in 2009 was confirmed by both frequent isolation of the virus from sewage as well as from patient samples in Finland. Over the last 10 years E-30 had only been isolated sporadically in Finland. We here study the phylogenetic relationships of the strains from the outbreak in the context of E-30 circulation over the last 20 years. The analyzed region comprised 276 nucleotides in the 5' end of VP1 (nucleotides 132-407 in the VP1 of the E-30 Bastianni strain). The Finnish strains were clustered into at least four distinct genogroups, with seven clusters exceeding the genotype demarcation of 12% and the 2009 epidemic strains forming the largest genogroup VII. Moreover, we detected largely divergent genotypes in 2007 and 2009. Interestingly, close genetic relatives of the epidemic strains had already been isolated a few years before the outbreak. Phylodynamic analysis estimated 8.9 years (95% highest posterior density intervals 7.0-11.0) as the age of genogroup VII, indicating a probable origin and evolutionary history prior to its introduction and epidemic expansion in Finland. Finally, the most recent common ancestor for the current E-30 diversity dates back to 1939 (95% highest posterior density intervals 1913-1956).