Genus-specific substitution rate variability among picornaviruses

Directed Evolution of Seneca Valley Virus in Tumorsphere and Monolayer Cell Cultures of a Small-Cell Lung Cancer Model

The Seneca Valley virus (SVV) is an oncolytic virus from the picornavirus family, characterized by a 7.3-kilobase RNA genome encoding for all the structural and functional viral proteins. Directed evolution by serial passaging has been employed for oncolytic virus adaptation to increase the killing efficacy towards certain types of tumors. We propagated the SVV in a small-cell lung cancer model under two culture conditions: conventional cell monolayer and tumorspheres, with the latter resembling more closely the cellular structure of the tumor of origin. We observed an increase of the virus-killing efficacy after ten passages in the tumorspheres. Deep sequencing analyses showed genomic changes in two SVV populations comprising 150 single nucleotides variants and 72 amino acid substitutions. Major differences observed in the tumorsphere-passaged virus population, compared to the cell monolayer, were identified in the conserved structural protein VP2 and in the highly variable P2 region, suggesting that the increase in the ability of the SVV to kill cells over time in the tumorspheres is acquired by capsid conservation and positively selecting mutations to counter the host innate immune responses.

Molecular evolutionary characteristics of SARS-CoV-2 emerging in the United States

SARS‐CoV‐2 is a newly discovered beta coronavirus at the end of 2019, which is highly pathogenic and poses a serious threat to human health. In this paper, 1875 SARS‐CoV‐2 whole genome sequences and the sequence coding spike protein (S gene) sampled from the United States were used for bioinformatics analysis to study the molecular evolutionary characteristics of its genome and spike protein. The MCMC method was used to calculate the evolution rate of the whole genome sequence and the nucleotide mutation rate of the S gene. The results showed that the nucleotide mutation rate of the whole genome was 6.677 × 10⁻⁴ substitution per site per year, and the nucleotide mutation rate of the S gene was 8.066 × 10⁻⁴ substitution per site per year, which was at a medium level compared with other RNA viruses. Our findings confirmed the scientific hypothesis that the rate of evolution of the virus gradually decreases over time. We also found 13 statistically significant positive selection sites in the SARS‐CoV‐2 genome. In addition, the results showed that there were 101 nonsynonymous mutation sites in the amino acid sequence of S protein, including seven putative harmful mutation sites. This paper has preliminarily clarified the evolutionary characteristics of SARS‐CoV‐2 in the United States, providing a scientific basis for future surveillance and prevention of virus variants.

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.

The Effect of Sample Bias and Experimental Artefacts on the Statistical Phylogenetic Analysis of Picornaviruses

Statistical phylogenetic methods are a powerful tool for inferring the evolutionary history of viruses through time and space. The selection of mathematical models and analysis parameters has a major impact on the outcome, and has been relatively well-described in the literature. The preparation of a sequence dataset is less formalized, but its impact can be even more profound. This article used simulated datasets of enterovirus sequences to evaluate the effect of sample bias on picornavirus phylogenetic studies. Possible approaches to the reduction of large datasets and their potential for introducing additional artefacts were demonstrated. The most consistent results were obtained using “smart sampling”, which reduced sequence subsets from large studies more than those from smaller ones in order to preserve the rare sequences in a dataset. The effect of sequences with technical or annotation errors in the Bayesian framework was also analyzed. Sequences with about 0.5% sequencing errors or incorrect isolation dates altered by just 5 years could be detected by various approaches, but the efficiency of identification depended upon sequence position in a phylogenetic tree. Even a single erroneous sequence could profoundly destabilize the whole analysis by increasing the variance of the inferred evolutionary parameters.

Epidemiology and molecular characterization of Porcine teschovirus in Hunan, China

Porcine teschoviruses (PTVs) have been shown to be widely distributed in pig populations. In this study, 261 faecal and 91 intestinal content samples collected from pigs at 29 farms in Hunan, China, were tested for the presence of PTV by reverse transcription-polymerase chain reaction (RT-PCR). An overall PTV-positivity rate of 19.03% was detected by RT-PCR, and a high PTV infection rate was circulating in asymptomatic fattening and nursery pigs. In total, 40 PTV isolates (PTV-HuNs) were obtained. Alignment of their coding sequences with those of other known PTVs revealed that the genomic sequence of the polyprotein contains 6,606-6,621 nucleotides, encoding a 2,202-2,207-amino acid sequence. Phylogenetic analyses based on the VP1 gene and capsid protein gene exhibited 13 main lineages corresponding to PTV serotypes 1-13, and seven PTV serotypes (PTV 2-6, 9, and 11) were identified in the isolates obtained in our study; this is the first report of PTV 5, 9 and 11 in China. Recombination analysis among the PTV-HuNs indicated that nine recombination events have occurred, including both inter- and intraserotype events. In addition, results demonstrated that only limited positive selection is acting on the global population of PTV isolates, and purifying selection is predominant. In conclusion, this study revealed a high infection rate of PTVs circulating in asymptomatic fattening and nursery pigs. The 40 PTV-HuNs showed high genetic diversity, and genetic analysis of all available PTV sequences revealed that strong purifying selection and recombination play important roles in the genetic diversity and evolution of the virus.

Genome wide analysis of the evolution of Senecavirus A from swine clinical material and assembly yard environmental samples

Senecavirus A (SVA), previously known as Seneca Valley virus, was first isolated in the United States in 2002. SVA was associated with porcine idiopathic vesicular disease in Canada and the USA in 2007 and 2012, respectively. Recent increase in SVA outbreaks resulting in neonatal mortality of piglets and/or vesicular lesions in sows in Brazil, the USA and Canada point to the necessity to study the pathogenicity and molecular epidemiology of the virus. Here, we report the analysis of the complete coding sequences of SVA from 2 clinical cases and 9 assembly yard environmental samples collected in 2015 in Canada, along with 22 previously released complete genomes in the GenBank. With this combined data set, the evolution of the SVA over a 12-month period in 2015/2016 was evaluated. These SVA isolates were characterized by a rapid accumulation of genetic variations driven mainly by a high nucleotide substitution rate and purifying selection. The SVA sequences clustered in clearly defined geographical areas with reported cases of SVA infection. No transmission links were identified between assembly yards, suggesting that point source introductions may have occurred. In addition, 25 fixed non-synonymous mutations were identified across all analyzed strains when compared to the prototype SVA strain (SVV-001). This study highlights the importance of monitoring SVA mutations for their role in increased virulence and impact on SVA diagnostics.

Impact of genetic changes, pathogenicity and antigenicity on Enterovirus- A71 vaccine development

Enterovirus-A71 (EV-A71) is an etiological agent of the hand, foot and mouth disease (HFMD). EV-A71 infection produces high fever and ulcers in children. Some EV-A71 strains produce severe infections leading to pulmonary edema and death. Although the protective efficacy of the inactivated vaccine (IV) was ≥90% against mild HFMD, there was approximately 80% protection against severe HFMD. The monovalent EV-A71 IV elicits humoral immunity but lacks long-term immunogenicity. Spontaneous mutations of the EV-A71 genome could lead to antigenicity changes and the virus may not be neutralized by antibodies elicited by the IV. A better alternative would be the live attenuated vaccine (LAV) that elicits cellular and humoral immunity. The LAV induces excellent antigenicity and chances of reversion is reduced by presence of multiple mutations which could reduce pathogenicity. Besides CV-A16, outbreaks have been caused by CV-A6 and CV-A10, hence the development of bivalent and trivalent vaccines is required.

Molecular epidemiology of senecavirus A associated with vesicular disease in pigs in Brazil

Senecavirus A (SV-A) may cause vesicular disease and neonatal mortality in pigs, and was first detected in Brazil in 2015. Samples including tissues and serum from pigs with suspected vesicular diseases were collected from January to August in 2015 from farms in the states of Minas Gerais, Santa Catarina, Goiás and Rio Grande do Sul, Brazil, and tested for the presence of SV-A by reverse transcriptase PCR. All samples were negative for foot and mouth disease virus, as well as 13 other infectious agents associated with vesicular diseases in pigs. SV-A was detected by PCR in 65/265 (24.5%) specimens. A 530 base pair fragment sequenced from the VP1 protein coding region indicated a high genetic distance from SV-A in other countries, but a common origin among the Brazilian isolates.

Phylodynamic Characterization of an Ocular-Tropism Coxsackievirus A24 Variant

Recent phylodynamic studies have focused on using tree topology patterns to elucidate interactions among the epidemiological, evolutionary, and demographic characteristics of infectious agents. However, because studies of viral phylodynamics tend to focus on epidemic outbreaks, tree topology signatures of tissue-tropism pathogens might not be clearly identified. Therefore, this study used a novel Bayesian evolutionary approach to analyze the A24 variant of coxsackievirus (CV-A24v), an ocular-tropism agent of acute hemorrhagic conjunctivitis. Analyses of the 915-nucleotide VP1 and 690-nt 3Dpol regions of 21 strains isolated in Taiwan and worldwide during 1985-2010 revealed a clear chronological trend in both the VP1 and 3Dpol phylogenetic trees: the emergence of a single dominant cluster in each outbreak. The VP1 sequences included three genotypes: GI (prototype), GIII (isolated 1985-1999), and GIV (isolated after 2000); no VP1 sequences from GII strains have been deposited in GenBank. Another five genotypes identified in the 3Dpol region had support values >0.9. Geographic and demographic transitions among CV-A24v clusters were clearly identified by Bayes algorithm. The transmission route was mapped from India to China and then to Taiwan, and each prevalent viral population declined before new clusters emerged. Notably, the VP1 and 3Dpol genes had high nucleotide sequence similarities (94.1% and 95.2%, respectively). The lack of co-circulating lineages and narrow tissue tropism affected the CV-A24v gene pool.

Evolution and phylogeographic dissemination of endemic porcine picornaviruses in Vietnam

Members of the Picornaviridae are important and often zoonotic viruses responsible for a variety of human and animal diseases. However, the evolution and spatial dissemination of different picornaviruses circulating in domestic animals are not well studied. We examined the rate of evolution and time of origin of porcine enterovirus G (EV-G) and porcine kobuvirus species C lineages (PKV-C) circulating in pig farms in Vietnam and from other countries. We further explored the spatiotemporal spread of EV-G and PKV-C in Southwest Vietnam using phylogeographic models. Multiple types of EV-G are co-circulating in Vietnam. The two dominant EV-G types among isolates from Vietnam (G1 and G6) showed strong phylogenetic clustering. Three clades of PKV-C (PKV-C1-3) represent more recent introductions into Vietnam; PKV-C2 is closely related to PKV-C from Southwest China, indicating possible cross-border dissemination. In addition, high virus lineage migration rates were estimated within four districts in Dong Thap province in Vietnam for both EV-G types (G1, G6) and all PKV-C (C1-3) clades. We found that Chau Thanh district is a primary source of both EV-G and PKV-C clades, consistent with extensive pig trading in and out of the district. Understanding the evolution and spatial dissemination of endemic picornaviruses in pigs may inform future strategies for the surveillance and control of picornaviruses.

Distribution and genetic characterization of Enterovirus G and Sapelovirus A in six Spanish swine herds

The prevalence of Enterovirus G (EV-G) and Sapelovirus A (PSV-1) was investigated in Spanish swine herds by means of cross-sectional studies. Faecal samples from clinically healthy pigs were collected from six farms, and analysed by RT-PCR. The results indicated a high prevalence of EV-G detected in nearly all the animals older than 3 weeks of age. Otherwise, PSV-1 was only detected in 3-week-old piglets from one of the farms. Genetic analyses performed in the VP1 region of the EV-G indicated circulation of diverse strains in the same farm, related to genotypes G1, G2, G3, G4, G6, G9, G12, G13 and G14. Moreover, co-infection of several PSV-1 variants in the same animal was evident, typical of viral quasispecies. Evolutionary pressure analysis indicated that microevolution of PSV-1 seems to be driven by negative selection. This study gives further insights in the epidemiology of EV-G and PSV-1.

Phylodynamic reconstruction of the spatiotemporal transmission and demographic history of coxsackievirus B2

BACKGROUND

Studies regarding coxsackievirus (CV) tend to focus on epidemic outbreaks, an imbalanced topology is considered to be an indication of acute infection with partial cross-immunity. In enteroviruses, a clear understanding of the characteristics of tree topology, transmission, and its demographic dynamics in viral succession and circulation are essential for identifying prevalence trends in endemic pathogens such as coxsackievirus B2 (CV-B2). This study applied a novel Bayesian evolutionary approach to elucidate the phylodynamic characteristics of CV-B2. A dataset containing 51 VP1 sequences and a dataset containing 34 partial 3D(pol) sequencing were analyzed, where each dataset included Taiwan sequences isolated during 1988-2013.

RESULTS

Four and five genotypes were determined based on the 846-nucleotide VP1 and 441-nucleotide 3D(pol) (6641-7087) regions, respectively, with spatiotemporally structured topologies in both trees. Some strains with tree discordance indicated the occurrence of recombination in the region between the VP1 and 3D(pol) genes. The similarities of VP1 and 3D(pol) gene were 80.0%-96.8% and 74.7%-91.9%, respectively. Analyses of population dynamics using VP1 dataset indicated that the endemic CV-B2 has a small effective population size. The balance indices, high similarity, and low evolutionary rate in the VP1 region indicated mild herd immunity selection in the major capsid region.

CONCLUSIONS

Phylodynamic analysis can reveal demographic trends and herd immunity in endemic pathogens.

Gene Diversification of an Emerging Pathogen: A Decade of Mutation in a Novel Fish Viral Hemorrhagic Septicemia (VHS) Substrain since Its First Appearance in the Laurentian Great Lakes

Viral Hemorrhagic Septicemia virus (VHSv) is an RNA rhabdovirus, which causes one of the world's most serious fish diseases, infecting >80 freshwater and marine species across the Northern Hemisphere. A new, novel, and especially virulent substrain—VHSv-IVb—first appeared in the Laurentian Great Lakes about a decade ago, resulting in massive fish kills. It rapidly spread and has genetically diversified. This study analyzes temporal and spatial mutational patterns of VHSv-IVb across the Great Lakes for the novel non-virion (Nv) gene that is unique to this group of novirhabdoviruses, in relation to its glycoprotein (G), phosphoprotein (P), and matrix (M) genes. Results show that the Nv-gene has been evolving the fastest (k = 2.0x10^-3 substitutions/site/year), with the G-gene at ~1/7 that rate (k = 2.8x10^-4). Most (all but one) of the 12 unique Nv- haplotypes identified encode different amino acids, totaling 26 changes. Among the 12 corresponding G-gene haplotypes, seven vary in amino acids with eight total changes. The P- and M- genes are more evolutionarily conserved, evolving at just ~1/15 (k = 1.2x10^-4) of the Nv-gene’s rate. The 12 isolates contained four P-gene haplotypes with two amino acid changes, and six M-gene haplotypes with three amino acid differences. Patterns of evolutionary changes coincided among the genes for some of the isolates, but appeared independent in others. New viral variants were discovered following the large 2006 outbreak; such differentiation may have been in response to fish populations developing resistance, meriting further investigation. Two 2012 variants were isolated by us from central Lake Erie fish that lacked classic VHSv symptoms, having genetically distinctive Nv-, G-, and M-gene sequences (with one of them also differing in its P-gene); they differ from each other by a G-gene amino acid change and also differ from all other isolates by a shared Nv-gene amino acid change. Such rapid evolutionary differentiation may allow new viral variants to evade fish host recognition and immune responses, facilitating long-time persistence along with expansion to new geographic areas.

Species-specific transmission of novel picornaviruses in lemurs

The roles of host genetics versus exposure and contact frequency in driving cross-species transmission remain the subject of debate. Here, we used a multitaxon lemur collection at the Saint Louis Zoo in the United States as a model to gain insight into viral transmission in a setting of high interspecies contact. Lemurs are a diverse and understudied group of primates that are highly endangered. The speciation of lemurs, which are endemic to the island of Madagascar, occurred in geographic isolation apart from that of continental African primates. Although evidence of endogenized viruses in lemur genomes exists, no exogenous viruses of lemurs have been described to date. Here we identified two novel picornaviruses in fecal specimens of ring-tailed lemurs (Lemur catta) and black-and-white ruffed lemurs (Varecia variegata). We found that the viruses were transmitted in a species-specific manner (lesavirus 1 was detected only in ring-tailed lemurs, while lesavirus 2 was detected only in black-and-white ruffed lemurs). Longitudinal sampling over a 1-year interval demonstrated ongoing infection in the collection. This was supported by evidence of viral clearance in some animals and new infections in previously uninfected animals, including a set of newly born triplets that acquired the infection. While the two virus strains were found to be cocirculating in a mixed-species exhibit of ring-tailed lemurs, black-and-white ruffed lemurs, and black lemurs, there was no evidence of cross-species transmission. This suggests that despite high-intensity contact, host species barriers can prevent cross-species transmissions of these viruses.

IMPORTANCE

Up to 75% of emerging infectious diseases in humans today are the result of zoonotic transmission. However, a challenge in understanding transmission dynamics has been the limited models of cross-species transmission. Zoos provide a unique opportunity to explore parameters defining viral transmission. We demonstrated that ongoing virus transmission in a mixed lemur species exhibit was species specific. This suggests that despite high contact intensity, host species barriers contribute to protection from cross-species transmission of these viruses. While the combinations of species might differ, most zoological parks worldwide commonly feature mixed-species exhibits. Collectively, this report demonstrates a widely applicable approach toward understanding infectious disease transmission.

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.

An evolutionary analysis of the Secoviridae family of viruses

The plant-infecting Secoviridae family of viruses forms part of the Picornavirales order, an important group of non-enveloped viruses that infect vertebrates, arthropods, plants and algae. The impact of the secovirids on cultivated crops is significant, infecting a wide range of plants from grapevine to rice. The overwhelming majority are transmitted by ecdysozoan vectors such as nematodes, beetles and aphids. In this study, we have applied a variety of computational methods to examine the evolutionary traits of these viruses. Strong purifying selection pressures were calculated for the coat protein (CP) sequences of nine species, although for two species evidence of both codon specific and episodic diversifying selection were found. By using Bayesian phylogenetic reconstruction methods CP nucleotide substitution rates for four species were estimated to range from between 9.29×10⁻³ to 2.74×10⁻³ (subs/site/year), values which are comparable with the short-term estimates of other related plant- and animal-infecting virus species. From these data, we were able to construct a time-measured phylogeny of the subfamily Comovirinae that estimated divergence of ninety-four extant sequences occurred less than 1,000 years ago with present virus species diversifying between 50 and 250 years ago; a period coinciding with the intensification of agricultural practices in industrial societies. Although recombination (modularity) was limited to closely related taxa, significant and often unique similarities in the protein domains between secovirid and animal infecting picorna-like viruses, especially for the protease and coat protein, suggested a shared ancestry. We discuss our results in a wider context and find tentative evidence to indicate that some members of the Secoviridae might have their origins in insects, possibly colonizing plants in a number of founding events that have led to speciation. Such a scenario; virus infection between species of different taxonomic kingdoms, has significant implications for virus emergence.

Cell tropism predicts long-term nucleotide substitution rates of mammalian RNA viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.

RNA viruses are the fastest evolving human pathogens, making their treatment and control difficult. Compared to DNA viruses, RNA viruses replicate with much lower fidelity, which can explain why RNA viruses evolve significantly faster than most DNA viruses. However, there is tremendous variation among the evolutionary rates of different RNA viruses, which is not explained by variation in mutation rates. Here we present a survey of mammalian RNA virus rates of evolution, and a comprehensive comparison of these rates to different properties of virus genomic architecture and ecology. We found that cell tropism is the most significant predictor of long-term rates of mammalian RNA virus evolution. For instance, viruses targeting epithelial cells evolve significantly faster than viruses that target neurons. Our results provide mechanistic insight into why viruses that infect respiratory and gastrointestinal epithelia have been difficult to control.

Viral substitution rate variation can arise from the interplay between within-host and epidemiological dynamics

The evolutionary rates of RNA viruses can differ from one another by several orders of magnitude. Much of this variation has been explained by differences in viral mutation rates and selective environments. However, substitution rates also vary considerably across viral populations belonging to the same species. In particular, viral lineages from epidemic regions tend to have higher substitution rates than those from endemic regions, and lineages from populations with higher contact rates tend to have higher substitution rates than those from populations with lower contact rates. We address the mechanism behind these patterns by using a nested modeling approach, whereby we integrate within-host viral replication dynamics with a population-level epidemiological model. Through numerical simulations and analytical approximations, we show that variation in viral substitution rates over the course of an infection, coupled with differences in age of infection of transmitting hosts under different epidemiological scenarios, can explain these evolutionary patterns. We further derive analytical estimates of expected substitution rate differences under epidemic versus endemic epidemiological conditions. By comparing these estimates to empirical data for four viral species, we show that these factors are sufficient to explain observed variation in substitution rates in three of four cases. This work shows that even in neutrally evolving viral populations, epidemiological dynamics can alter substitution rates via the interplay between within-host replication dynamics and population-level disease dynamics.

Recombination in the evolution of human rhinovirus genomes

Human rhinoviruses (HRV) are highly prevalent human respiratory pathogens that belong to the genus Enterovirus. Although recombination within the coding region is frequent in other picornavirus groups, most evidence of recombination in HRV has been restricted to the 5' untranslated region. We analysed the occurrence of recombination within published complete genome sequences of members of all three HRV species and additionally compared sequences from HRV strains spanning 14 years. HRV-B and HRV-C showed very little evidence of recombination within the coding region. In contrast, HRV-A sequences appeared to have undergone a large number of recombination events, typically involving whole type groups. This suggests that HRV-A may have been subject to extensive recombination during the period of diversification into types. This study demonstrates the rare and sporadic nature of contemporary recombination of HRV strains and contrasts with evidence of extensive recombination within HRV-A and between members of different species during earlier stages in its evolutionary diversification.

One misdated sequence of rabbit hemorrhagic disease virus prevents accurate estimation of its nucleotide substitution rate

BACKGROUND

The literature is ripe with phylogenetic estimates of nucleotide substitution rates, especially of measurably evolving species such as RNA viruses. However, it is not known how robust these rate estimates are to inaccuracies in the data, particularly in sampling dates that are used for molecular clock calibration. Here we report on the rate of evolution of the emerging pathogen Rabbit hemorrhagic disease virus (RHDV), which has significantly different rates of evolution for the same outer capsid (VP60) gene published in the literature. In an attempt to reconcile the conflicting data and further elucidate details of RHDV 's evolutionary history, we undertook fresh Bayesian analyses and employed jackknife control methods to produce robust substitution rate and time to most recent common ancestor (TMRCA) estimates for RHDV based on the VP60 and RNA-dependent RNA polymerase genes.

RESULTS

Through these control methods, we were able to identify a single misdated taxon, a passaged lab strain used for vaccine production, which was responsible for depressing the RHDV capsid gene's rate of evolution by 65%. Without this isolate, the polymerase and the capsid protein genes had nearly identical rates of evolution: 1.90x10-3 nucleotide substitutions/site/year, ns/s/y, (95% highest probability density (HPD) 1.25x10-3-2.55x10-3) and 1.91x10-3 ns/s/y (95% HPD 1.50x10-3-2.34x10-3), respectively.

CONCLUSIONS

After excluding the misdated taxon, both genes support a significantly higher substitution rate as well as a relatively recent emergence of RHDV, and obviate the need for previously hypothesized decades of unobserved diversification of the virus. The control methods show that using even one misdated taxon in a large dataset can significantly skew estimates of evolutionary parameters and suggest that it is better practice to use smaller datasets composed of taxa with unequivocal isolation dates. These jackknife controls would be useful for future tip-calibrated rate analyses that include taxa with ambiguous dates of isolation.

Molecular epidemiology and evolution of human enterovirus serotype 68 in Thailand, 2006-2011

Background

Publications worldwide have reported on the re-occurrence of human enterovirus 68 (EV68), a rarely detected pathogen usually causing respiratory illness. However, epidemiological data regarding this virus in particular on the Asian continent has so far been limited.

Methodology/Findings

We investigated the epidemiology and genetic variability of EV68 infection among Thai children with respiratory illnesses from 2006–2011 (n = 1810). Semi-nested PCR using primer sets for amplification of the 5′-untranslated region through VP2 was performed for rhino-enterovirus detection. Altogether, 25 cases were confirmed as EV68 infection indicating a prevalence of 1.4% in the entire study population. Interestingly, the majority of samples were children aged >5 years (64%). Also, co-infection with other viruses was found in 28%, while pandemic H1N1 influenza/2009 virus was the most common co-infection. Of EV68-positive patients, 36% required hospitalizations with the common clinical presentations of fever, cough, dyspnea, and wheezing. The present study has shown that EV68 was extremely rare until 2009 (0.9%). An increasing annual prevalence was found in 2010 (1.6%) with the highest detection frequency in 2011 (4.3%). Based on analysis of the VP1 gene, the evolutionary rate of EV68 was estimated at 4.93×10⁻³ substitutions/site/year. Major bifurcation of the currently circulating EV68 strains occurred 66 years ago (1945.31 with (1925.95–1960.46)95% HPD). Among the current lineages, 3 clusters of EV68 were categorized based on the different molecular signatures in the BC and DE loops of VP1 combined with high posterior probability values. Each cluster has branched off from their common ancestor at least 36 years ago (1975.78 with (1946.13–1984.97)95% HPD).

Conclusion

Differences in epidemiological characteristic and seasonal profile of EV68 have been found in this study. Results from Bayesian phylogenetic investigations also revealed that EV68 should be recognized as a genetically diverse virus with a substitution rate identical to that of enterovirus 71 genotype B (4.2×10⁻³ s/s/y).

Partitioning the genetic diversity of a virus family: approach and evaluation through a case study of picornaviruses

The recent advent of genome sequences as the only source available to classify many newly discovered viruses challenges the development of virus taxonomy by expert virologists who traditionally rely on extensive virus characterization. In this proof-of-principle study, we address this issue by presenting a computational approach (DEmARC) to classify viruses of a family into groups at hierarchical levels using a sole criterion-intervirus genetic divergence. To quantify genetic divergence, we used pairwise evolutionary distances (PEDs) estimated by maximum likelihood inference on a multiple alignment of family-wide conserved proteins. PEDs were calculated for all virus pairs, and the resulting distribution was modeled via a mixture of probability density functions. The model enables the quantitative inference of regions of distance discontinuity in the family-wide PED distribution, which define the levels of hierarchy. For each level, a limit on genetic divergence, below which two viruses join the same group, was objectively selected among a set of candidates by minimizing violations of intragroup PEDs to the limit. In a case study, we applied the procedure to hundreds of genome sequences of picornaviruses and extensively evaluated it by modulating four key parameters. It was found that the genetics-based classification largely tolerates variations in virus sampling and multiple alignment construction but is affected by the choice of protein and the measure of genetic divergence. In an accompanying paper (C. Lauber and A. E. Gorbalenya, J. Virol. 86:3905-3915, 2012), we analyze the substantial insight gained with the genetics-based classification approach by comparing it with the expert-based picornavirus taxonomy.

Characterization of a canine homolog of human Aichivirus

Many of our fatal "civilization" infectious diseases have arisen from domesticated animals. Although picornaviruses infect most mammals, infection of a companion animal is not known. Here we describe the identification and genomic characterization of the first canine picornavirus. Canine kobuvirus (CKoV), identified in stool samples from dogs with diarrhea, has a genomic organization typical of a picornavirus and encodes a 2,469-amino-acid polyprotein flanked by 5' and 3' untranslated regions. Comparative phylogenetic analysis using various structural and nonstructural proteins of CKoV confirmed it as the animal virus homolog most closely related to human Aichivirus (AiV). Bayesian Markov chain Monte Carlo analysis suggests a mean recent divergence time of CKoV and AiV within the past 20 to 50 years, well after the domestication of canines. The discovery of CKoV provides new insights into the origin and evolution of AiV and the species specificity and pathogenesis of kobuviruses.