Evolution of enterovirus 70 in nature: all isolates were recently derived from a common ancestor

Genetic and phenotypic characterization of recently discovered enterovirus D type 111

Members of the species Enterovirus D (EV-D) remain poorly studied. The two first EV-D types (EV-D68 and EV-D70) have regularly caused outbreaks in humans since their discovery five decades ago but have been neglected until the recent occurrence of severe respiratory diseases due to EV-D68. The three other known EV-D types (EV-D94, EV-D111 and EV-D120) were discovered in the 2000s-2010s in Africa and have never been observed elsewhere. One strain of EV-D111 and all known EV-D120s were detected in stool samples of wild non-human primates, suggesting that these viruses could be zoonotic viruses. To date, EV-D111s are only known through partial genetic sequences of the few strains that have been identified so far. In an attempt to bring new pieces to the puzzle, we genetically characterized four EV-D111 strains (among the seven that have been reported until now). We observed that the EV-D111 strains from human samples and the unique simian EV-D111 strain were not phylogenetically distinct, thus suggesting a recent zoonotic transmission. We also discovered evidences of probable intertypic genetic recombination events between EV-D111s and EV-D94s. As recombination can only happen in co-infected cells, this suggests that EV-D94s and EV-D111s share common replication sites in the infected hosts. These sites could be located in the gut since the phenotypic analysis we performed showed that, contrary to EV-D68s and like EV-D94s, EV-D111s are resistant to acid pHs. We also found that EV-D111s induce strong cytopathic effects on L20B cells, a cell line routinely used to specifically detect polioviruses. An active circulation of EV-D111s among humans could then induce a high number of false-positive detection of polioviruses, which could be particularly problematic in Central Africa, where EV-D111 circulates and which is a key region for poliovirus eradication.

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.

Enterovirus 70 receptor utilization is controlled by capsid residues that also regulate host range and cytopathogenicity

Enterovirus type 70, an etiologic agent of acute hemorrhagic conjunctivitis, may bind different cellular receptors depending on cell type. To understand how EV70-receptor interaction is controlled, we studied two variants of the virus with distinct receptor utilization. EV70-Rmk, derived by passage in rhesus monkey kidney cells, replicates poorly in HeLa cells and does not cause cytopathic effects. Decay accelerating factor (DAF) is not a cell receptor for EV70-Rmk. Passage of EV70-Rmk in HeLa cells lead to isolation of EV70-Dne, which does not replicate in rhesus monkey kidney cells but grows to high titers in HeLa cells and causes cytopathic effects. DAF is sufficient for cell entry of EV70-Dne. EV70-Rmk replicates in human eye and brain-derived cell lines, whereas the Dne strain replicates only in HeLa cells and in conjunctiva-derived 15C4 cells. The two EV70 strains differ by five amino acid changes in the viral capsid. Single substitution of four of the five EV70-Rmk amino acids with the residue from EV70-Dne leads to lytic replication in HeLa cells. Conversely, substitution of any of the five EV70-Dne amino acids with the EV70-Rmk amino acid does not alter replication in HeLa cells. Three of these capsid amino acids are predicted to be located in the canyon encircling the fivefold axis of symmetry, one amino acid is found at the fivefold axis of symmetry, and one is located the interior of the capsid. The five EV70 residues define a region of the capsid that controls viral host range, DAF utilization, and cytopathogenicity.

Enterovirus 94, a proposed new serotype in human enterovirus species D

The genus Enterovirus (family Picornaviridae) contains five species with strains isolated from humans: Human enterovirus A (HEV-A), HEV-B, HEV-C, HEV-D and Poliovirus. In this study, a proposed new serotype of HEV-D was characterized. Four virus strains were isolated from sewage in Egypt and one strain from acute flaccid paralysis cases in the Democratic Republic of the Congo. The complete genome of one environmental isolate, the complete coding sequence of one clinical isolate and complete VP1 regions from the other isolates were sequenced. These isolates had 66.6–69.4 % nucleotide similarity and 74.7–76.6 % amino acid sequence similarity in the VP1 region with the closest enterovirus serotype, enterovirus 70 (EV70), suggesting that the isolates form a new enterovirus type, tentatively designated enterovirus 94 (EV94). Phylogenetic analyses including sequences of the 5′ UTR, VP1 and 3D regions demonstrated that EV94 isolates formed a monophyletic group within the species HEV-D. No evidence of recombination was found between EV94 and the other HEV-D serotypes, EV68 and EV70. Further biological characterization showed that EV94 was acid stable and had a wide cell tropism in vitro. Attempts to prevent replication with protective antibodies to known enterovirus receptors (poliovirus receptor, vitronectin α v β 3 receptor and decay accelerating factor) were not successful. Seroprevalence studies in the Finnish population revealed a high prevalence of this virus over the past two decades.

[Development and role of comparative sequence analysis in medical virology]

INTRODUCTION

Development of the polymerase chain reaction and deoxyribonucleic acid sequencing techniques has enabled precise identification, classification and taxonomy of viruses. COMPARATIVE SEQUEENCE ANALYSIS: (Comparative sequence analysis methods can be used in medical virology for many practical purposes. They may be classified into three broad categories: I - reconstruction of genealogical relationships between individual viral isolatesfor detection and monitoring of sources, reservoirs and modes of viral transmission; II - virus genotyping, that is determination of relationships between genetic types of viruses and their phenotypic properties, which has important implications for immunoprophylaxis, therapy and prognosis of viral diseases, and III investigation of functional properties of defined viral sequences, of special importance for explanation of viral pathogenesis and design of antiviral drugs.

FUTURE PROSPECTS

The combination of DNA sequencing with polymerase chain reaction following reverse transcription with the use of random primers offers a universal means for diagnosis of viral infections.

Molecular typing of enteroviruses: current status and future requirements. The European Union Concerted Action on Virus Meningitis and Encephalitis

Human enteroviruses have traditionally been typed according to neutralization serotype. This procedure is limited by the difficulty in culturing some enteroviruses, the availability of antisera for serotyping, and the cost and technical complexity of serotyping procedures. Furthermore, the impact of information derived from enterovirus serotyping is generally perceived to be low. Enteroviruses are now increasingly being detected by PCR rather than by culture. Classical typing methods will therefore no longer be possible in most instances. An alternative means of enterovirus typing, employing PCR in conjunction with molecular genetic techniques such as nucleotide sequencing or nucleic acid hybridization, would complement molecular diagnosis, may overcome some of the problems associated with serotyping, and would provide additional information regarding the epidemiology and biological properties of enteroviruses. We argue the case for developing a molecular typing system, discuss the genetic basis of such a system, review the literature describing attempts to identify or classify enteroviruses by molecular methods, and suggest ways in which the goal of molecular typing may be realized.

Evolutionary analysis of the picornavirus family

An exhaustive evolutionary analysis of the picornavirus family has been carried out using the amino acid sequences of several proteins of the viruses including: the capsid proteins (1D, 1B, and 1C) situated at the 5' end of the genome and responsible for the serotype of the viruses, and the viral polymerase (3D), located at the 3' end of the genome. The evolutionary relationships found among the viruses studied support the new classification, recently suggested, in contrast to the classical one, and the existence of a new genus for the picornavirus family. In the new taxonomic organization, five genera form the picornavirus family: (1) aphthoviruses, (2) cardioviruses, (3) hepatoviruses (previously classified as enteroviruses), (4) renteroviruses (which mainly constitute a combination of the previous genera rhinovirus and enterovirus), and (5) a new genus, with a new and unique representative: the echovirus 22. Our analysis also allowed us, for the first time, to propose the most probable sequence of speciation events to have given rise to the current picornavirus family. The bootstrap procedure was used to check the reliability of the phylogenetic trees obtained. The application of the method of the statistical geometry in distance space to internal branches of the tree revealed a high degree of evolutionary "noise," which makes the resolution of some internal branching points difficult.

The Soriano Award Lecture. Emerging infections of the nervous system

The epidemic of acquired immunodeficiency disease [AIDS] has focused interest on the origins of "new" infectious agents. Great plagues are well known from the distant past, but a number of novel diseases affecting the nervous system infections have emerged in recent years. The causes of such new disorders are diverse: whereas rapid mutations of microbes allow the evolution of truly novel agents, the appearance of new diseases is more often due to changes in human or vector populations or changes in societal mores that result in dissemination of preexistent microbes. Examples of recently emerging infections that involve the nervous system include the enterovirus 70 epidemics with poliomyelitis-like disease, the appearance of California virus encephalitis in the midwestern United States, the rapid spread of Lyme disease with its many neurological complications in the eastern United States, and the outbreak of bovine spongiform encephalopathy in the United Kingdom, in addition to the devastating epidemic of human immunodeficiency virus (HIV), which will cause nervous system disease in over half of those infected. As the world population increases and modern transportation brings us closer into a "global village" more new agents will emerge and more will be sustained. Knowledge of the molecular biology and ecology of the agents and awareness of how our actions can alter their behavior are our best defense.

Molecular evolution of the major capsid protein VP1 of enterovirus 70

Nucleotide sequences of the genome RNA encoding capsid protein VP1 (918 nucleotides) of 18 enterovirus 70 (EV70) isolates collected from various parts of the world in 1971 to 1981 were determined, and nucleotide substitutions among them were studied. The genetic distances between isolates were calculated by the pairwise comparison of nucleotide difference. Regression analysis of the genetic distances against time of isolation of the strains showed that the synonymous substitution rate was very high at 21.53 x 10(-3) substitution per nucleotide per year, while the nonsynonymous rate was extremely low at 0.32 x 10(-3) substitution per nucleotide per year. The rate estimated by the average value of synonymous and nonsynonymous substitutions (W.-H. Li, C.-C. Wu, and C.-C. Luo, Mol. Biol. Evol. 2:150-174, 1985) was 5.00 x 10(-3) substitution per nucleotide per year. Taking the average value of synonymous and nonsynonymous substitutions as genetic distances between isolates, the phylogenetic tree was inferred by the unweighted pairwise grouping method of arithmetic average and by the neighbor-joining method. The tree indicated that the virus had evolved from one focal place, and the time of emergence was estimated to be August 1967 +/- 15 months, 2 years before first recognition of the pandemic of acute hemorrhagic conjunctivitis. By superimposing every nucleotide substitution on the branches of the phylogenetic tree, we analyzed nucleotide substitution patterns of EV70 genome RNA. In synonymous substitutions, the proportion of transitions, i.e., C<==>U and G<==>A, was found to be extremely frequent in comparison with that reported on other viruses or pseudogenes. In addition, parallel substitutions (independent substitutions at the same nucleotide position on different branches, i.e., different isolates, of the tree) were frequently found in both synonymous and nonsynonymous substitutions. These frequent parallel substitutions and the low nonsynonymous substitution rate despite the very high synonymous substitution rate described above imply a strong restriction on nonsynonymous substitution sites of VP1, probably due to the requirement for maintaining the rigid icosahedral conformation of the virus.

Viral evolution and insects as a possible virologic turning table

Three lines of observation demonstrate the role of arthropods in transmission and evolution of viruses. a) Recent outbreaks of viruses from their niches took place and insects have played a major role in propagating the viruses. b) Examination of the list of viral families and their hosts shows that many infect invertebrates (I) and vertebrates (V) or (I) and plants (P) or all kingdoms (VIPs). This notion holds true irrespective of the genome type. At first glance the argument seems to be weak in the case of enveloped and non-enveloped RNA viruses with single-stranded (ss) segmented or non-segmented genomes of positive (+) or negative polarity. Here, there are several families infecting V or P only; no systematic relation to arthropods is found. c) In the non-enveloped plant viruses with ss RNA genomes there is a strong tendency for segmentation and individual packaging of the genome pieces. This is in contrast to ss+ RNA animal viruses and can only be explained by massive transmission by seed or insects or both, because individual packaging necessitates a multihit infection. Comparisons demonstrate relationships in the nonstructural proteins of double-stranded and ss+ RNA viruses irrespective of host range, segmentation, and envelope. Similar conclusions apply for the negative-stranded RNA viruses. Thus, viral supergroups can be created that infect V or P and exploit arthropods for infection or transmission or both. Examples of such relationships and explanations for viral evolution are reviewed and the arthropod orders important for cell culture are given.

The complete nucleotide sequence of a variant of Coxsackievirus A24, an agent causing acute hemorrhagic conjunctivitis

The complete nucleotide sequence was determined for the cDNAs that represent the RNA genome of the standard strain of a variant of coxsackievirus A24, the EH24/70, one of the agents causing acute hemorrhagic conjunctivitis. The genome is 7461 nucleotide long and is polyadenylated at the 3'-end terminus. Following a 750-nucleotide 5'-noncoding region, there was a long open reading frame of 6642 nucleotides, which serve to encode a viral polyprotein consisting of 2214 amino acids. Comparison of the deduced amino acid sequence of the polyprotein with those of known enteroviruses allowed us to predict the possible cleavage sites. The overall structure and the organization of the RNA genome is typical for an enterovirus. Based on the similarity of the nucleotide sequence of the 5' and 3' noncoding regions, together with the amino-acid sequence of the encoded proteins, EH24/70 appeared to be closely related to polioviruses and coxsackievirus A21.

Phylogenetic analysis of a coxsackievirus A24 variant: the most recent worldwide pandemic was caused by progenies of a virus prevalent around 1981

Nucleotide substitutions in the viral-encoded proteinase 3C (3Cpro) region (549 nucleotides) of the RNA genome of a coxsackievirus A24 variant (CA24v), one of the agents causing acute hemorrhagic conjunctivitis (AHC), were studied using 32 isolates collected from the Eastern hemisphere in 1970-1989. Based on regression analysis of nucleotide differences among isolates, the nucleotide substitution rate of CA24v 3Cpro was estimated to be 3.7 x 10(-3)/nucleotide/year. A phylogenetic tree constructed by the modified unweighted pair group method using arithmetic averages (UPGMA) indicated that CA24v had evolved from a common ancestor which appeared in one focal place in November 1963 +/- 21 months, about 7 years before the first isolation of CA24v in Singapore. The tree also revealed that all the recent epidemic isolates in 1985-1989 including Asian and Ghanian strains diverged from each other after 1981. This finding is consistent with the evidence that AHC due to CA24v had been confined to Southeast Asia and the Indian subcontinent until 1985, then suddenly and explosively spread to other areas where no CA24v isolations had been reported.

Phylogenetically different strains of a variant of coxsackievirus A 24 were repeatedly introduced but discontinued circulating in Japan

Variations in the nucleotide sequence of 3 C proteinase of coxsackievirus A 24 variant (CA 24 v) were analyzed to define the route of transmission and spread of the virus which was introduced to Japan on three separate occasions, 1985-86, 1988, and 1989. The nucleotide sequences of isolates from the same year's outbreak in Japan were identical or closely related, while the isolates from different outbreaks were less closely related to one another than to those from other countries in the same year. All Japanese isolates from Okinawa and other prefectures in 1985 and 1986 were closely related to the Taiwan strains in those same years, indicating common-source outbreaks. Two 1988 isolates from Chiba Prefecture, Japan, were closely related to those from Singapore in 1987, China in 1988 and Hong Kong in 1988. All seven Japanese isolates from Chiba Prefecture in 1989 comprised a group together with the Taiwan and Singapore strains in 1988. The results indicate that CA 24 v was introduced into Japan on each occasion from the outside. Furthermore, in contrast to the explosive epidemics in Okinawa Prefecture in 1985 and 1986, the virus which was repeatedly introduced to other areas in Japan did not circulate endemically, and disappeared within a short time.

Rapid antigenic-type replacement and DNA sequence evolution of canine parvovirus

Analysis of canine parvovirus (CPV) isolates with a panel of monoclonal antibodies showed that after 1986, most viruses isolated from dogs in many parts of the United States differed antigenically from the viruses isolated prior to that date. The new antigenic type (designated CPV type 2b) has largely replaced the previous antigenic type (CPV type 2a) among virus isolates from the United States. This represents the second occurrence of a new antigenic type of this DNA virus since its emergence in 1978, as the original CPV type (CPV type 2) had previously been replaced between 1979 and 1981 by the CPV type 2a strain. DNA sequence comparisons showed that CPV types 2b and 2a differed by as few as two nonsynonymous (amino acid-changing) nucleotide substitutions in the VP-1 and VP-2 capsid protein genes. One mutation, resulting in an Asn-Asp difference at residue 426 in the VP-2 sequence, was shown by comparison with a neutralization-escape mutant selected with a non-CPV type 2b-reactive monoclonal antibody to determine the antigenic change. The mutation selected by that monoclonal antibody, a His-Tyr difference in VP-2 amino acid 222, was immediately adjacent to residue 426 in the three-dimensional structure of the CPV capsid. The CPV type 2b isolates are phylogenetically closely related to the CPV type 2a isolates and are probably derived from a common ancestor. Phylogenetic analysis showed a progressive evolution away from the original CPV type. This pattern of viral evolution appears most similar to that seen in some influenza A viruses.

The nucleotide sequence of 3C proteinase region of the coxsackievirus A24 variant: comparison of the isolates in Taiwan in 1985-1988

Acute hemorrhagic conjunctivitis caused by coxsackievirus A24 variant (CA24v) first appeared in Taiwan in October 1985, followed by two other sequential epidemics in 1986 and 1988. In order to know the evolutionary relationship of the CA24v strains isolated in Taiwan, we first determined the nucleotide sequence of the 3C proteinase (3Cpro) region of the prototype strain (EH24/70), isolated in Singapore in 1970, by molecular cloning. The nucleotide sequence of the 3Cpro region thus sequenced showed striking homology with polioviruses and coxsackievirus A21. Viral RNA of eight isolates obtained from the three epidemics was reverse transcribed, amplified by the polymerase chain reaction, and cloned into M13 phage for the production of ssDNA for nucleotide sequencing by the dideoxy chain termination method. When the number of nucleotide difference was taken as a genetic distance between isolates, all isolates showed a very similar distance from the EH24/70, the earliest isolate of CA24v, indicating that they evolved at a constant evolutionary rate. Phylogenetic analysis by the unweighted pairwise grouping method of arithmetic average (UPGMA) indicated that the six isolates collected in 1985 and 1986 were closely related, while two 1988 isolates were more distant from them. The branching time between these two groups was estimated to be May 1984, 18 months before the first recognition of the CA24v epidemic in Taiwan. This is the first report of the nucleotide sequence of CA24v genome RNA and of an evolutionary analysis of the virus using the nucleotide sequence.

A study of the evolution of coxsackievirus A24 variant in Ghana by viral RNA fingerprinting analysis

An epidemic of acute haemorrhagic conjunctivitis (AHC) caused by coxsackievirus A24 variant (CA24v) was reported in Accra, Ghana in May 1987. We studied 7 of the viral strains collected from May to November, 1987, by RNA genome fingerprinting. Pairwise comparisons of the oligonucleotide maps showed that genetic similarity among them ranged from between 60.0 to 84.7%. Using base sequence variations deduced from genetic similarity among the isolates, isolation time of the strains and the rate of nucleotide substitution (estimated in a previous paper, Miyamura et al., 1990), we calculated divergence times and constructed a phylogenetic tree. This tree indicated that all the 7 strains had diverged from each other from 11 to 26 months before the AHC epidemic in Accra. CA24v may have been introduced into the country or the neighbouring area, at least, more than two years earlier, i.e. in the early half of 1985.

Evolutionary study on the Coxsackievirus A 24 variant causing acute hemorrhagic conjunctivitis by oligonucleotide mapping analysis of RNA genome

The evolution of the variant of Coxsackievirus A 24 (CA 24 v) which causes acute hemorrhagic conjunctivitis was explored. Using 15 isolates obtained from Southeast Asia during the period 1970-1986, the genetic distance between isolates was estimated from pairwise comparison of nucleotide changes deduced from common spots on oligonucleotide maps of the isolates. From regression analysis of the genetic distance and the time of isolation of the isolates, the evolutionary rate of CA 24 v was estimated to be 3.44 x 10(-4)/nucleotide/month. The phylogenetic relationship of these isolates was explored using the neighbor-joining method and the modified unweighted pair group method using arithmetic averages (UPGMA). The phylogenetic tree constructed indicates that CA 24 v appeared from one focal place in July 1968 +/- 25 months, very close to the time of the first world epidemic of, then newly recognized, acute hemorrhagic conjunctivitis.

Neutralizing monoclonal antibody against enterovirus-70 reacts with viral proteins 1C and 1D

A library of murine monoclonal antibodies against the prototype Enterovirus-70 (EV-70) strain, J670/71, was made for the purpose of studying the immunologically reactive determinants of the virus. Each of the monoclonal antibodies reacted with several other strains of Enterovirus-70 when tested by immunofluorescence. However, none of these monoclonal antibodies reacted with any other picornavirus tested. It was found that all of the monoclonal antibodies precipitated EV-70 viral proteins 1C and 1D in radio-immunoprecipitation assays. However, only one of these monoclonal antibodies, an IgG3 kappa, was capable of neutralizing the virus.

Antigenic relationships between avian paramyxoviruses. III. A mathematical model of antigenic drift and a computer-assisted approach for construction of a phylogenetic tree

The suggested model of antigenic kinship between related paramyxoviruses is based on another concept of antigenic determinant, as compared to the previously suggested combinatorial mathematical model by the authors. According to it, antigenic changes of any determinant do not proceed by "leaps" but can be changed gradually. Such changed determinant can induce a correspondingly changed type of antibodies which still preserve a certain kinship to the original type of the determinant (before its changing) revealed by cross reaction serological tests. Accordingly, there can be "families" of the determinants differing by degree of relatedness to (or, reversely, by antigenic distance from) the "original" ("ancestor") determinant. In addition to another interpretation of the antigenic kinship, the new mathematical model was used as an approach for revealing phylogenetic relationships between antigenically related viruses.

Correlation between plaque size and genetic variation of type 3 poliovirus from a vaccinate

Vaccine strain derivatives could be serially isolated from the feces of a healthy infant for about 30 days after the second administration of a trivalent poliovirus vaccine. These were all identified as poliovirus type 3, and no other types were detected. The titer in feces increased gradually until day 19, when it reached a maximum, after which it decreased. The plaque size altered roughly in parallel with the changes in the titer. On day 11 it began to increase, becoming the highest abruptly on day 15, and then fell again. Oligonucleotide mapping analysis of RNAs of isolated viruses showed that the samples from days 1 and 11 exhibited the same patterns as those shown by the vaccine strain RNA. However, on day 15 many of the preexisting spots had decreased, and new ones appeared. On day 27, the pattern was again similar to that of the original virus, rather than that of the day-15 virus. A nonmetric distance scaling and cluster analysis suggested that all the strains were derived from the same origin and that the polioviruses on days 1, 11, 21, and 27 had closer relationships, but the poliovirus on day 15 was different.