Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

The rate of molecular evolution of alpha-fetoprotein approaches that of pseudogenes

We conducted the present study in an attempt to correlate function with the rate of molecular evolution for serum albumin and alpha-fetoprotein. We found a high rate of silent substitution (between 5 X 10(-9) and 7 X 10(-9)/site/year) for both the albumin and alpha-fetoprotein genes, perhaps the highest so far reported for an expressed nuclear gene. The rates of effective substitution and amino acid changes were also very high, but in contrast to silent substitutions, they are higher for alpha-fetoprotein than for albumin by approximately 70%. For alpha-fetoprotein, the rate of effective substitution (1.5 X 10(-9)/site/year) may be approaching that for nonfunctional pseudogenes (about 3 X 10(-9)/site/year). Evolutionary divergence was also estimated at the amino acid level. It was found that the rate of change of alpha-fetoprotein (55% amino acids replaced in 100 Myr) approaches that of the fastest-evolving fibrinopeptides (92% amino acids replaced in 100 Myr). This high rate may indicate that alpha-fetoprotein can tolerate a great deal of molecular variation without its function being impaired in the process. Albumin evolves at a slower rate (39% amino acids replaced in 100 Myr), although still faster than either hemoglobin (17% amino acids replaced in 100 Myr) or cytochrome c (5% amino acids replaced in 100 Myr). The slower evolutionary rate may indicate that albumin has more refined functional specifications and hence can tolerate fewer mutational changes. The latter conclusion remains, however, to be reconciled with the condition of inherited analbuminemia, where a virtually complete absence of albumin produces surprisingly few symptoms.     

浙江地区青霉素耐药肺炎链球菌PBPs基因及氨基酸序列研究

为阐明温州地区青霉素耐药肺炎链球菌（PRSP）的青霉素结合蛋白（PBPs）的基因和氨基酸序列的变异特点,对温州医学院自2000年11月～2004年1月收集的26份肺炎链球菌进行分离、鉴定及青霉素药敏实验,并对每株链球菌的PBP1A、PBP2B、PBP2X基因进行PCR扩增和直接测序,通过序列比对与生物信息学分析。结果表明,研究中的PBP1A的主要突变位点是保守基序KTG之后的4个连续氨基酸替换Thr574Ala、Ser575Thr、Gln576Gly、Phe577Tyr和保守序列STMK内的氨基酸替换Thr371Ser;PBP2B的主要突变位点是保守序列SSN之后的氨基酸替换Thr451Ala;PBP2X的主要突变位点是保守基序STMK 内的氨基酸替换Thr338Ala。以上突变类型以及菌株的青霉素耐药水平与文献报道相符。研究检测的PRSP的PBPs基因中暂未发现本地区特有的（新的）基因突变,也未检测出文献报道的某些与青霉素抗性相关的氨基酸替换。     

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.     

Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.     

Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp nasonia (hymenoptera: pteromalidae)

We sequenced the nearly complete mtDNA of 3 species of parasitic wasps, Nasonia vitripennis (2 strains), Nasonia giraulti, and Nasonia longicornis, including all 13 protein-coding genes and the 2 rRNAs, and found unusual patterns of mitochondrial evolution. The Nasonia mtDNA has a unique gene order compared with other insect mtDNAs due to multiple rearrangements. The mtDNAs of these wasps also show nucleotide substitution rates over 30 times faster than nuclear protein-coding genes, indicating among the highest substitution rates found in animal mitochondria (normally <10 times faster). A McDonald and Kreitman test shows that the between-species frequency of fixed replacement sites relative to silent sites is significantly higher compared with within-species polymorphisms in 2 mitochondrial genes of Nasonia, atp6 and atp8, indicating directional selection. Consistent with this interpretation, the Ka/Ks (nonsynonymous/synonymous substitution rates) ratios are higher between species than within species. In contrast, cox1 shows a signature of purifying selection for amino acid sequence conservation, although rates of amino acid substitutions are still higher than for comparable insects. The mitochondrial-encoded polypeptides atp6 and atp8 both occur in F0F1ATP synthase of the electron transport chain. Because malfunction in this fundamental protein severely affects fitness, we suggest that the accelerated accumulation of replacements is due to beneficial mutations necessary to compensate mild-deleterious mutations fixed by random genetic drift or Wolbachia sweeps in the fast evolving mitochondria of Nasonia. We further propose that relatively high rates of amino acid substitution in some mitochondrial genes can be driven by a "Compensation-Draft Feedback"; increased fixation of mildly deleterious mutations results in selection for compensatory mutations, which lead to fixation of additional deleterious mutations in nonrecombining mitochondrial genomes, thus accelerating the process of amino acid substitutions.     

Adaptive molecular evolution of virulence genes of avian influenza - A virus subtype H5N1: An analysis of host radiation

The phenomenon of host radiation is strongly influenced by the rates of mutation of their virulence genes. We have studied the molecular evolution of virulence genes (HA, NS, PB2) of the Avian Influenza Virus H5N1 from avian to human hosts. We used a site-specific comparison of synonymous (silent) and non-synonymous (amino acid altering) nucleotide substitutions for the three chosen genes in parasite populations from different hosts. Analyses were made using Maximum Likelihood (ML) genealogies for the null and alternate hypothesis based on differential gamma distribution rates. The null hypothesis had a higher rate of substitution and was found to be more suitable for all the studied genes by Likelihood Ratio Test (LRT). The study showed the NS gene to be having the fastest rate of evolution.     

Complete structure of A/duck/Ukraine/63 influenza hemagglutinin gene: animal virus as progenitor of human H3 Hong Kong 1968 influenza hemagglutinin

We have explored the possibility that an animal viral reservoir contained a direct ancestor gene for the H3 hemagglutinin type present in influenza A viruses in humans since 1968. For this purpose, the duck/Ukraine/1/63 hemagglutinin gene was cloned and sequenced. From the comparison of its complete primary structure with that of several human H3 hemagglutinins as well as those of an H2 and an H7 hemagglutinin, we conclude that the duck/Ukraine/63 hemagglutinin sequence fully corroborates its previous identification by immunological and other methods as belonging to the H3 subtype. Moreover, the duck/Ukraine/63 amino acid sequence is more closely related structurally and presumably antigenically to the human Aichi/68 hemagglutinin, which formed the beginning of the H3N2 pandemic in humans, than to that of Victoria/75, which has undergone an additional 7 year drift period in humans. This observation could best be explained by a common ancestor hemagglutinin gene for duck/Ukraine/63 and human Aichi/68. On the basis of silent, accumulated base changes, we estimate that the strain carrying this postulated common progenitor hemagglutinin gene was circulating in the period 1949–1953 in the animal reservoir. This relatively recent divergence, as well as the closer kinship between the duck/Ukraine/63 and the human Aichi/68 hemagglutinin, as compared with the later Victoria/75, strongly suggests that the influenza A virus of the H3N2 subtype circulating in the human population since 1968 has derived its hemagglutinin gene from a strain in the animal reservoir. Undoubtedly, this occurred by reassortment between previously present human H2N2 virus and this animal strain. These results provide support at the molecular level for the general idea that the wide variety of influenza viruses known to be present in animals can serve as a gene reservoir for human influenza A viruses.     

广东人禽流感H5N1毒株M基因特性、进化和变异

通过对人禽流感H5N1毒株M基因序列的变异分析,揭示毒株M基因特征与进化。检测广东地区人禽流感H5N1毒株M基因核苷酸序列,同时检索全球人禽流感H5N1毒株M基因序列,采用DNAStar5.0软件对检索的人禽流感H5N1毒株M基因核苷酸序列进行比对和分析;并结合流行病学资料对变异毒株进行进化速度分析。结果发现,1997~2006年53株毒株M1基因和51株毒株M2核苷酸序列同源性均分成两组,1997年毒株为第一组(GⅠ),2003~2006年香港、越南、泰国、印尼、中国大陆、土耳其、伊拉克、阿塞拜疆、埃及毒株为第二组(GⅡ)。M1基因20个氨基酸位点置换,占7.94%(20/252),其中2003~2006年毒株M1基因有9个氨基酸位点不同于1997年毒株;M2基因22个氨基酸置换,占22.7%,其中2003~2006年毒株M2基因有4个氨基酸不同于1997年毒株。M2基因Ks值为26.8×10-6~42.6×10-6Nt/d,Ka值为4.39×10-6~6.98×10-6Nt/d;而M1基因的同义突变速度均远高于错义突变速度,显示M1基因受到机体免疫压力较小;检验发现M1基因进化存在负选择性压力。2003~2006年毒株M1基因通过氨基酸S224N置换,增加一个糖基化位点NSS224-226;而来自印尼的8株毒株M2基因发生C50F置换,引起蛋白二级结构改变。1997年中国香港人禽流感毒株自当时出现后,便未在以后人禽流感疫情中出现。2003~2006年毒株M1基因增加糖基化位点NSS224-226,可能与毒株致病性有关。人禽流感H5N1毒株M基因在自然界变异频繁,可能影响H5N1毒株的人-人传播能力。     

H5N1亚型禽流感病毒A/duck/Shandong/093/2004株的全基因克隆及序列分析

应用流感病毒通用引物［4］和H5N1亚型禽流感（Avian influenza virus, AIV）的型特异性引物,成功的扩增出H5N1亚型禽流感病毒A/duck/Shandong/093/2004株（简称A/D/SD/04）的全基因序列（包括5′和3′端的非编码区序列）。A/D/SD/04的基因组核苷酸全序列与18株网上公布的禽流感基因序列进行比较和分析,结果与4株鸭源H5N1的5～7个基因具99%以上的同源性;与14株H5N1有至少一个以上内部基因同源性在95%以上。与H9亚型AIV代表株A/Quail/Hongkong/G1/97（简称G1株）和A/Chicken/Beijing/1/94（简称BJ94）比较,除了非结构基因（Nonstructural gene, NS）与G1株的同源性为95.3%外,其余基因均在36.6%～92.1%之间。说明A/D/SD/04没有H9N2基因的直接整合,是H5N1毒株在自然界的重组株。推导的HA氨基酸序列分析,A/D/SD/04 的血凝素（Heamgglutinin,HA）裂解位点与比较的16株AIV的序列一致,是高致病性禽流感的分子特征（PQRERRRKKR/G）,第226位氨基酸是对禽类和哺乳细胞均具有亲嗜性的蛋氨酸（Met）。神经氨酸酶（Neuraminidase, NA）在第48位氨基酸（颈部）后有20个氨基酸的缺失,但非结构蛋白（NS）没有在79～84氨基酸发生缺失。碱性聚合酶2（PB2）的627位氨基酸是亲禽类细胞的谷氨酸（Glu, E）。结合生物学特性和分子特征,A/D/SD/04对小鼠的致病力是由多种因素决定,其可能是一株对鸡高度致病,并逐渐获得对哺乳动物致病能力的中间重组病毒。     

The causes of synonymous rate variation in the rodent genome. Can substitution rates be used to estimate the sex bias in mutation rate?

Miyata et al. have suggested that the male-to-female mutation rate ratio (alpha) can be estimated by comparing the neutral substitution rates of X-linked (X), Y-linked (Y), and autosomal (A) genes. Rodent silent site X/A comparisons provide very different estimates from X/Y comparisons. We examine three explanations for this discrepancy: (1) statistical biases and artifacts, (2) nonneutral evolution, and (3) differences in mutation rate per germline replication. By estimating errors and using a variety of methodologies, we tentatively reject explanation 1. Our analyses of patterns of codon usage, synonymous rates, and nonsynonymous rates suggest that silent sites in rodents are evolving neutrally, and we can therefore reject explanation 2. We find both base composition and methylation differences between the different sets of chromosomes, a result consistent with explanation 3, but these differences do not appear to explain the observed discrepancies in estimates of alpha. Our finding of significantly low synonymous substitution rates in genomically imprinted genes suggests a link between hemizygous expression and an adaptive reduction in the mutation rate, which is consistent with explanation 3. Therefore our results provide circumstantial evidence in favor of the hypothesis that the discrepancies in estimates of alpha are due to differences in the mutation rate per germline replication between different parts of the genome. This explanation violates a critical assumption of the method of Miyata et al., and hence we suggest that estimates of alpha, obtained using this method, need to be treated with caution.     

A maximum likelihood method for analyzing pseudogene evolution: implications for silent site evolution in humans and rodents.

We present a new likelihood method for detecting constrained evolution at synonymous sites and other forms of nonneutral evolution in putative pseudogenes. The model is applicable whenever the DNA sequence is available from a protein-coding functional gene, a pseudogene derived from the protein-coding gene, and an orthologous functional copy of the gene. Two nested likelihood ratio tests are developed to test the hypotheses that (1) the putative pseudogene has equal rates of silent and replacement substitutions; and (2) the rate of synonymous substitution in the functional gene equals the rate of substitution in the pseudogene. The method is applied to a data set containing 74 human processed-pseudogene loci, 25 mouse processed-pseudogene loci, and 22 rat processed-pseudogene loci. Using the informatics resources of the Human Genome Project, we localized 67 of the human-pseudogene pairs in the genome and estimated the GC content of a large surrounding genomic region for each. We find that, for pseudogenes deposited in GC regions similar to those of their paralogs, the assumption of equal rates of silent and replacement site evolution in the pseudogene is upheld; in these cases, the rate of silent site evolution in the functional genes is approximately 70% the rate of evolution in the pseudogene. On the other hand, for pseudogenes located in genomic regions of much lower GC than their functional gene, we see a sharp increase in the rate of silent site substitutions, leading to a large rate of rejection for the pseudogene equality likelihood ratio test.     

Evolutional analysis of human influenza A virus N2 neuraminidase genes based on the transition of the low-pH stability of sialidase activity

The 1957 and 1968 human pandemic influenza A virus strains as well as duck viruses possess sialidase activity under low-pH conditions, but human H3N2 strains isolated after 1968 do not possess such activity. We investigated the transition of avian (duck)-like low-pH stability of sialidase activities with the evolution of N2 neuraminidase (NA) genes in human influenza A virus strains. We found that the NA genes of H3N2 viruses isolated from 1971 to 1982 had evolved from the side branches of NA genes of H2N2 epidemic strains isolated in 1968 that were characterized by the low-pH-unstable sialidase activities, though the NA genes of the 1968 pandemic strains preserved the low-pH-stable sialidase. These findings suggest that the prototype of the H3N2 epidemic influenza strains isolated after 1968 probably acquired the NA gene from the H2N2 low-pH-unstable sialidase strain by second genetic reassortment in humans.     

Mammalian genes as molecular clocks?

In analyzing the silent nucleotide substitutions in some mammalian mitochondrial mRNA coding genes, we had found that the frequency of each of the four nucleotides in rat, mouse, and cow, but not in humans, is the same in the silent third codon position (Lanave C, Preparata G, Saccone C, Serio G (1984) J Mol Evol 20:86-93). Because our findings for these three species were compatible with a stationary Markov process for the evolution of nucleotide sequences, we applied such a model to calculate the effective evolutionary silent substitution rate (vs) and the divergence times among the species. In this paper we have analyzed the first and second codon positions in the same mammalian mitochondrial genes. We found that in the first and second codon positions the human mitochondrial genes satisfy the stationarity conditions. This has allowed us to use the stochastic model mentioned above to calculate the divergence times among mouse, rat, cow, and human. Furthermore, we have analyzed the silent substitution rate in one nuclear gene for these four mammals. We found that in this gene the effective silent substitution rate is about 3 times lower than in mitochondrial genes, and that humans are in this case stationary with respect to the other three mammals in the third codon position as well. Application of our Markov model to this latter gene yields divergence times consistent with our previous determinations.     

Immunological and genomic analyses of two serotypes of avian paramyxovirus isolated from wild ducks in Japan

A total of 11 avian paramyxoviruses isolated from migrating feral ducks in Niigata, Japan, were characterized by serological and genomic analyses. Hemagglutination inhibition and immuno-double-diffusion tests with antisera specific for the isolated hemagglutinin-neuraminidase polypeptides of reference strains indicated that, of these, eight isolates possessed hemagglutinin-neuraminidase antigen closely related to that of duck/Hong Kong/D3/75, and the remaining three isolates possessed antigen closely related to that of duck/Hong Kong/199/77. RNA analysis of the eight isolates identified serologically as duck/Hong Kong/D3/75 by oligonucleotide mapping revealed that these isolates were genetically very similar to each other but different from the reference strain and isolates reported previously. The oligonucleotide maps of duck/Hong Kong/199/77-like isolates appeared to be very similar to each other, suggesting the same origin, but not to the duck/Hong Kong/199/77 virus.     

The core 2 beta-1,6-N-acetylglucosaminyltransferase-mucin encoded by bovine herpesvirus 4 was acquired from an ancestor of the African buffalo

The Bo17 gene of bovine herpesvirus 4 (BoHV-4) is the only viral gene known to date that encodes a homologue of the cellular core 2 beta-1,6-N-acetylglucosaminyltransferase-mucin type (C2GnT-M). To investigate the origin and evolution of the Bo17 gene, we analyzed its distribution among BoHV-4 strains and determined the sequences of Bo17 from nine representative strains and of the C2GnT-M gene from six species of ruminants expected to encompass the group within which the gene acquisition occurred. Of 34 strains of BoHV-4, isolated from four different continents, all were found to contain the Bo17 gene. Phylogenetic analyses indicated that Bo17 was acquired from a recent ancestor of the African buffalo, implying that cattle subsequently acquired BoHV-4 by cross-species transmission. The rate of synonymous nucleotide substitution in Bo17 was estimated at 5 x 10(-8) to 6 x 10(-8) substitutions/site/year, consistent with previous estimates made under the assumption that herpesviruses have cospeciated with their hosts. The Bo17 gene acquisition was dated to around 1.5 million years ago. Bo17 sequences from BoHV-4 strains from African buffalo and from cattle formed two separate clades, estimated to have split about 700,000 years ago. Analysis of the ratio of nonsynonymous to synonymous nucleotide substitutions revealed a burst of amino acid replacements subsequent to the transfer of the cellular gene to the viral genome, followed by a return to a strong constraint on nonsynonymous changes during the divergence of contemporary BoHV-4 strains. The Bo17 gene represents the most recent of the known herpesvirus gene acquisitions and provides the best opportunity for learning more about this important process of viral evolution.     

Molecular evolution of hemagglutinin genes of H1N1 swine and human influenza A viruses

Summary The hemagglutinin (HA) genes of influenza type A (H1N1) viruses isolated from swine were cloned into plasmid vectors and their nucleotide sequences were determined. A phylogenetic tree for the HA genes of swine and human influenza viruses was constructed by the neighbor-joining method. It showed that the divergence between swine and human HA genes might have occurred around 1905. The estimated rates of synonymous (silent) substitutions for swine and human influenza viruses were almost the same. For both viruses, the rate of synonymous substitution was much higher than that of nonsynonymous (amino acid altering) substitution. It is the case even for only the antigenic sites of the HA. This feature is consistent with the neutral theory of molecular evolution. The rate of nonsynonymous substitution for human influenza viruses was three times the rate for swine influenza viruses. In particular, nonsynonymous substitutions at antigenic sites occurred less frequently in swine than in humans. The difference in the rate of nonsynonymous substitution between swine and human influenza viruses can be explained by the different degrees of functional constraint operating on the amino acid sequence of the HA in both hosts.     

Polymorphism and evolution of influenza A virus genes

The nucleotide sequences of four genes of the influenza A virus (nonstructural protein, matrix protein, and a few subtypes of hemagglutinin and neuraminidase) are compiled for a large number of strains isolated from various locations and years, and the evolutionary relationship of the sequences is investigated. It is shown that all of these genes or subtypes are highly polymorphic and that the polymorphic sequences (alleles) are subject to rapid turnover in the population, their average age being much less than that of higher organisms. Phylogenetic analysis suggests that most polymorphic sequences within a subtype or a gene appeared during the last 80 years and that the divergence among the subtypes of hemagglutinin genes might have occurred during the last 300 years. The high degree of polymorphism in this RNA virus is caused by an extremely high rate of mutation, estimated to be 0.01/nucleotide site/year. Despite the high rate of mutation, most influenza virus genes are apparently subject to purifying selection, and the rate of nucleotide substitution is substantially lower than the mutation rate. There is considerable variation in the substitution rate among different genes, and the rate seems to be lower in nonhuman viral strains than in human strains. The difference might be responsible for the so-called freezing effect in some viral strains.      

Nucleotide Sequence Analysis of Recent Epidemic Strains of Enterovirus 70

Nucleotide sequences of the genome of enterovirus 70 (EV70) isolated in Osaka in 1993 were determined, and compared with those of the past epidemic strains. Nucleotide substitution rates in 332 bp, 174 bp, and 178 bp of the genes encoding viral protein (VP)1, VP2, or VP3 were 9.0, 7.5, and 5.6% between Kumoi-2/93 and J670/71, respectively. Likewise, the putative amino acid substitution rates were 1.8, 0, and 0%. It seems that the epidemic strains of EV70 in Japan have been evolving at a constant high nucleotide substitution rate but almost all the substitutions were synonymous.