Genetic diversity of JC virus in the modern Filipino population: implications for the peopling of the Philippines

The Philippines is generally believed to have been established by various peoples who migrated from neighboring areas. To gain new insights into the peopling of the Philippines, we used the JC virus (JCV) genotyping approach. We collected about 50 urine samples on each of two representative islands of the Philippines, Luzon and Cebu. DNA was extracted from the urine samples and used to amplify the 610-bp region (IG region) of the viral genome. For each island, we determined about 20 IG sequences, from which a neighbor-joining phylogenetic tree was constructed to classify the JCV isolates detected into distinct genotypes. The predominant genotype detected was SC, the Southeast Asian genotype. Minor JCV genotypes were SC/Phi, B1-a, and B3. SC/Phi was a subcluster of SC and has not been detected in areas other than the Philippines. B1-a was detected previously in mainland China, Pamalican Island (Palawan, Philippines), and Taiwan (an aboriginal tribe). B3 was classified in this study into two subgroups, one (B3-a) containing three Luzon isolates and several Chinese, Thai, and Uzbek isolates, the other (B3-b) containing two Luzon, one Cebu, and one Indonesian isolate. These findings suggest that the modern Filipino population was formed not only by Southeast Asians carrying SC but also by a few distinct ethnic groups carrying SC/Phi, B1-a, and B3-a or -b.     

Genetic diversity of JC virus in the Saami and the Finns: implications for their population history

The JC virus (JCV) genotyping method was used to gain insights into the population history of the Saami and the Finns, both speaking Finno-Ugric languages and living in close geographic proximity. Urine samples from Saami and Finns, collected in northern and southern Finland, respectively, were used to amplify a 610-bp JCV-DNA region containing abundant type-specific mutations. Based on restriction site polymorphisms in the amplified fragments, we classified JCV isolates into one of the three superclusters of JCV, type A, B, or C. All 15 Saami isolates analyzed and 41 of 43 Finnish isolates analyzed were classified as type A, the European type, and two samples from Finns were classified as type B, the African/Asian type. We then amplified and sequenced a 583-bp JCV-DNA region from the type A isolates of Saami and Finns. According to type-determining nucleotides within the region, we classified type A isolates into EU-a1, -a2, or -b. Most type A isolates from Saami were classified as EU-a1, while type A isolates from Finns were distributed among EU-a1, EU-a2, and EU-b. This trend in the JCV-genotype distribution was statistically significant. On a phylogenetic tree based on complete sequences, most of the type A isolates from Saami were clustered in a single clade within EU-a1, while those from Finns were distributed throughout EU-a1, EU-a2, and EU-b. These findings are discussed in the context of the population history of the Saami and the Finns. This study provides new complete JCV DNA sequences derived from populations of anthropological interest.     

Evolution of Human Polyomavirus JC: Implications for the Population History of Humans

The polyomavirus JC virus (JCV), the etiological agent of progressive multifocal leukoencephalopathy, is ubiquitous in the human population, infecting children asymptomatically, then persisting in the kidney. The main mode of transmission of JCV is from parents to children through long-term cohabitation. Twelve JCV subtypes that occupy unique domains in Europe, Africa, and Asia have been identified. Here, we attempted to elucidate the evolutionary relationships among JCV strains worldwide using the whole-genome approach with which a highly reliable phylogeny of JCV strains can be reconstructed. Sixty-five complete JCV DNA sequences, derived from various geographical regions and belonging to 11 of the 12 known subtypes, were subjected to phylogenetic analysis using three independent methods: the neighbor-joining, maximum parsimony, and maximum likelihood methods. The trees obtained with these methods consistently indicated that ancestral JCVs were divided into three superclusters, designated as Types A, B, and C. A split in Type A generated two subtypes, EU-a and -b, mainly containing European and Mediterranean strains. The first split in Type B generated Af2 (the major African subtype). Subsequent splits in Type B generated B1-c (a minor European subtype) and all seven Asian subtypes (B1-a, -b, -d, B2, MY, CY, and SC). Type C generated a single subtype (Af1), consisting of strains derived from western Africa. While the present findings provided a basis on which to classify JCV into types or subtypes, they have several implications for the divergence and migration of human populations. Received: 4 April 2001 / Accepted: 31 July 2001     

Phylogenetic analysis of major African genotype (Af2) of JC virus: Implications for origin and dispersals of modern Africans

Both mtDNA and the Y chromosome have been used to investigate how modern humans dispersed within and out of Africa. This issue can also be studied using the JC virus (JCV) genotype, a novel marker with which to trace human migrations. Africa is mainly occupied by two genotypes of JCV, designated Af1 and Af2. Af1 is localized to central/western Africa, while Af2 is spread throughout Africa and in neighboring areas of Asia and Europe. It was recently suggested that Af1 represents the ancestral type of JCV, which agrees with the African origin of modern humans. To better understand the origin of modern Africans, we examined the phylogenetic relationships among Af2 isolates worldwide. A neighbor-joining phylogenetic tree was constructed based on the complete JCV DNA sequences of 51 Af2 isolates from Africa and neighboring areas. According to the resultant tree, Af2 isolates diverged into two major clusters, designated Af2-a and -b, with high bootstrap probabilities. Af2-a contained isolates mainly from South Africa, while Af2-b contained those from the other parts of Africa and neighboring regions of Asia and Europe. These findings suggest that Af2-carrying Africans diverged into two groups, one carrying Af2-a and the other carrying Af2-b; and that the former moved to southern Africa, while the latter dispersed throughout Africa and to neighboring regions of Asia and Europe. The present findings are discussed with reference to relevant findings in genetic and linguistic studies.     

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.     

Phylogenetic relationships among JC virus strains in Japanese/Koreans and Native Americans speaking Amerind or Na-Dene

Many genetic studies using human mtDNA or the Y chromosome have been conducted to elucidate the relationships among the three Native American groups speaking Amerind, Na-Dene, and Eskimo-Aleut. Human polyomavirus JC (JCV) may also help to gain insights into this issue. JCV isolates are classified into more than 10 geographically distinct genotypes (designated subtypes here), which were generated by splits in the three superclusters, Types A, B, and C. A particular subtype of JCV (named MY) belonging to Type B is spread in both Japanese/Koreans and Native Americans speaking Amerind or Na-Dene. In this study, we evaluated the phylogenetic relationships among MY isolates worldwide, using the whole-genome approach, with which a highly reliable phylogeny of JCV isolates can be reconstructed. Thirty-six complete sequences belonging to MY (10 from Japanese/Koreans, 24 from Native Americans, and 2 from others), together with 54 belonging to other subtypes around the world, were aligned and subjected to phylogenetic analysis using the neighbor-joining and maximum-likelihood methods. In the resultant phylogenetic trees, the MY sequences diverged into two Japanese/Korean and five Native American clades with high bootstrap probabilities. Two of the Native American clades contained isolates mainly from Na-Denes and the others contained isolates mainly from Amerinds. The Na-Dene clades were not clustered together, nor were the Amerind clades. In contrast, the two Japanese/Korean clades were clustered at a high bootstrap probability. We concluded that there is no distinction between Amerinds and Na-Denes in terms of indigenous JCVs, although they are linguistically distinguished from each other.     

Phylogenetic Relationships Among JC Virus Strains in Japanese/Koreans and Native Americans Speaking Amerind or Na-Dene

Abstract Many genetic studies using human mtDNA or the Y chromosome have been conducted to elucidate the relationships among the three Native American groups speaking Amerind, Na-Dene, and Eskimo-Aleut. Human polyomavirus JC (JCV) may also help to gain insights into this issue. JCV isolates are classified into more than 10 geographically distinct genotypes (designated subtypes here), which were generated by splits in the three superclusters, Types A, B, and C. A particular subtype of JCV (named MY) belonging to Type B is spread in both Japanese/Koreans and Native Americans speaking Amerind or Na-Dene. In this study, we evaluated the phylogenetic relationships among MY isolates worldwide, using the whole-genome approach, with which a highly reliable phylogeny of JCV isolates can be reconstructed. Thirty-six complete sequences belonging to MY (10 from Japanese/Koreans, 24 from Native Americans, and 2 from others), together with 54 belonging to other subtypes around the world, were aligned and subjected to phylogenetic analysis using the neighbor-joining and maximum-likelihood methods. In the resultant phylogenetic trees, the MY sequences diverged into two Japanese/Korean and five Native American clades with high bootstrap probabilities. Two of the Native American clades contained isolates mainly from Na-Denes and the others contained isolates mainly from Amerinds. The Na-Dene clades were not clustered together, nor were the Amerind clades. In contrast, the two Japanese/Korean clades were clustered at a high bootstrap probability. We concluded that there is no distinction between Amerinds and Na-Denes in terms of indigenous JCVs, although they are linguistically distinguished from each other.     

Molecular Evolution of the Small Subunit of Ribulose Bisphosphate Carboxylase: Nucleotide Substitution and Gene Conversion

The nucleotide sequences encoding the mature portion of 31 ribulose 1.5-bisphosphate carboxylase small subunit (SSU) genes from 17 genera of plants, green algae and cyanobacteria were examined. Among the 465 pairwise sequence comparisons, SSU multigene family members within the same species were more similar to each other in nonsynonymous or replacement nucleotide substitutions (RNS) than they were to SSU sequences in any other organism. The concerted evolution of independent SSU gene lineages within closely related plant species suggests that homogenization of RNS positions has occurred at least once in the life of each genus. The rate of expected RNS among mature SSU sequences was calculated to be 1.25 X 10(-9)/site/yr for the first 70 million years (MY) of divergence with a significant slowing to 0.13 X 10(-9)/site/yr for the next 1,400 MY. The data suggest that mature SSU sequences do not accumulate more than 20% differences in the RNS positions without compensatory changes in other components of this enzyme system. During the first 70 MY of divergence between species, the rate of expected synonymous or silent nucleotide substitutions (SNS) is approximately 6.6 X 10(-9)/site/yr. This is five times the RNS rate and is similar to the silent rate observed in animals. In striking contrast, SNS and RNS do not show this correlation among SSU gene family members within a species. A mechanism involving gene conversion within the exons followed by selection for biased gene conversion products with conservation of RNS positions and divergence of SNS positions is discussed. A SSU gene tree based on corrected RNS for 31 SSU sequences is presented and agrees well with a species tree based on morphological and cytogenetic traits for the 17 genera examined. SSU gene comparisons may be useful in predicting phylogenetic relationships and in some cases divergence times of various plant, algal and cyanobacterial species.     

东北地区棕背鼠平种群遗传结构及亚种分类

Using partial sequences of mtDNA control region,we sought genetic diversity,population genetic structure and subspecies classification in 5 populations of Myodes rufocanus from Northeast China and Hokkaido,Japan. After sequencing a 616 bp segment of control region,we detected 45 haplotypes based on 57 variable sites among 55 individuals. Our results revealed that high-level haplotype diversity and slightly less nucleotide diversity occurred in these populations. AMOVA analysis showed that high-level genetic differentiation happened between populations due to 64.70% genetic diversity and Fst ranging from 0. 240 to 0. 814. In a phylogenetic analysis,the haplotypes from Northeast China split into two groups:Daxing’anling-Xiaoxing’anling-Wandashan and Hailin-Dahailin. This suggested that we can separate the populations from Northeast China into two subspecies,M . r. irkutensis and M. r. Changbaishanensis respectively corresponding with the former and latter areas.     

中国东北地区远东鼩鼱遗传多样性与遗传结构

对中国东北地区3个种群(大兴安岭、小兴安岭和长白山山脉)的远东鼩鼱(77个样本)Cyt b基因全序列进行分析,共获得64个单倍型。整体单倍型多态性为0.9920,核苷酸多态性为0.0105,表明该地区远东鼩鼱具有较高遗传多样性,且长白山山脉远东鼩鼱种群遗传多样性明显高于大兴安岭和小兴安岭种群。F-统计量、遗传相似系数和遗传距离分析结果均显示,种群间和采样地间的遗传距离与地理距离基本相符。方差分析显示,种群间的变异占总变异的33.4%,种群内的采样地间变异占总变异的10.2%,采样点内部变异占总变异的56.4%。种群历史分析显示,东北地区远东鼩鼱未经历过数量扩张。从GenBank下载了欧亚其他地区远东鼩鼱序列进行遗传结构研究。远东鼩鼱系统发生树分化为2大分支:一大支主要由大兴安岭和小兴安岭种群构成,两种群具有一定分化;另一大支又分为两个分支。中介网络图显示,远东鼩鼱具有3个谱系:一个谱系主要由大兴安岭和小兴安岭的单倍型样本构成,还包括长白山山脉的4个单倍型样本;另一谱系包括来自于中国东北地区3个种群的个别单倍型,还包括俄罗斯贝加尔湖单倍型和芬兰单倍型;最后一个谱系完全是由长白山山脉单倍型构成。遗传多样性、系统发生树和中介网络图结果均表明,长白山山脉为远东鼩鼱末次冰期避难所。     

Pattern of synonymous and nonsynonymous substitutions: An indicator of mechanisms of molecular evolution

Comparison of numbers of synonymous and nonsynonymous substitutions is useful for understanding mechanisms of molecular evolution. In this paper, I examine the statistical properties of six methods of estimating numbers of synonymous and nonsynonymous substitutions. The six methods are Miyata and Yasunaga’s (MY) method; Nei and Gojobori’s (NG) method; Li, Wu and Luo’s (LWL) method; Pamilo, Bianchi and Li’s (PBL) method; and Ina’s (Ina) two methods. When the transition/transversion bias at the mutation level is strong, the numbers of synonymous and nonsynonymous substitutions are estimated more accurately by the PBL and Ina methods than by the NG, MY and LWL methods. When the nucleotide-frequency bias is strong and distantly related sequences are compared, all the six methods give underestimates of the number of synonymous substitutions. The concept of synonymous and nonsynonymous categories is also useful for analysis of DNA polymorphism data.     

New methods for estimating the numbers of synonymous and nonsynonymous substitutions

New methods for estimating the numbers of synonymous and nonsynonymous substitutions per site were developed. The methods are unweighted pathway methods based on Kimura's two-parameter model. Computer simulations were conducted to evaluate the accuracies of the new methods, Nei and Gojobori's (NG) method, Miyata and Yasunaga's (MY) method, Li, Wu, and Luo's (LWL) method, and Pamilo, Bianchi, and Li's (PBL) method. The following results were obtained: (1) The NG, MY, and LWL methods give overestimates of the number of synonymous substitutions and underestimates of the number of nonsynonymous substitutions. The major cause for the biased estimation is that these three methods underestimate the number of synonymous sites and overestimate the number of nonsynonymous sites. (2) The PBL method gives better estimates of the numbers of synonymous and nonsynonymous substitutions than those obtained by the NG, MY, and LWL methods. (3) The new methods also give better estimates of the numbers of synonymous and nonsynonymous substitutions than those obtained by the NG, MY, and LWL methods. In addition, estimates of the numbers of synonymous and nonsynonymous sites obtained by the new methods are reasonably accurate. (4) In some cases, the new methods and the PBL method give biased estimates of substitution numbers. However, from the number of nucleotide substitutions at the third position of codons, we can examine whether estimates obtained by the new methods are good or not, whereas we cannot make an examination of estimates obtained by the PBL method. (5) When there are strong transition/transversion and nucleotide-frequency biases like mitochondrial genes, all of the above methods give biased estimates of substitution numbers. In such cases, Kondo et al.'s method is recommended to be used for estimating the number of synonymous substitutions, although their method cannot estimate the number of nonsynonymous substitutions and is time-consuming. These results, particularly result (1), call for reexaminations of some genes. This is because evolutionary pictures of genes have often been discussed on the basis of results obtained by the NG, MY, and LWL methods, which are favorable for the neutral theory of molecular evolution.     

Evidence for natural selection on leukocyte immunoglobulin-like receptors for HLA class I in Northeast Asians

Human leukocyte antigen (HLA) plays a critical role in innate and adaptive immunity and is a well-known example of genes under natural selection. However, the genetic aspect of receptors recognizing HLA molecules has not yet been fully elucidated. Leukocyte immunoglobulin (Ig)-like receptors (LILRs) are a family of HLA class I-recognizing receptors comprising activating and inhibitory forms. We previously reported that the allele frequency of the 6.7 kb LILRA3 deletion is extremely high (71%) in the Japanese population, and we identified premature termination codon (PTC)-containing alleles. In this study, we observed a wide distribution of the high deletion frequency in Northeast Asians (84% in Korean Chinese, 79% in Man Chinese, 56% in Mongolian, and 76% in Buryat populations). Genotyping of the four HapMap populations revealed that LILRA3 alleles were in strong linkage disequilibrium with LILRB2 alleles in Northeast Asians. In addition, PTC-containing LILRA3 alleles were detected in Northeast Asians but not in non-Northeast Asians. Furthermore, flow-cytometric analysis revealed that the LILRB2 allele frequent in Northeast Asians was significantly associated with low levels of expression. F(ST) and extended-haplotype-homozygosity analysis for the HapMap populations provided evidence of positive selection acting on the LILRA3 and LILRB2 loci. Taken together, our results suggest that both the nonfunctional LILRA3 alleles and the low-expressing LILRB2 alleles identified in our study have increased in Northeast Asians because of natural selection. Our findings, therefore, lead us to speculate that not only HLA class I ligands but also their receptors might be sensitive to the local environment.     

Asian genotypes of JC virus in Japanese-Americans suggest familial transmission

To examine the mode of JC virus (JCV) transmission, we collected urine samples from second- and third-generation Japanese-Americans in Los Angeles, Calif., whose parents and grandparents were all Japanese. From the urine samples of these Japanese-Americans, we mainly detected two subtypes (CY and MY) of JCV that are predominantly found among native Japanese. This finding provides support for the hypothesis that JCV is transmitted mainly within the family through long-term cohabitation.     

Apparent Selection Intensity for the Cytochrome Oxidase Subunit I Gene Varies with Mode of Reproduction in Echinoderms

When most amino acid substitutions in protein-coding genes are slightly deleterious rather than selectively neutral, life history differences can potentially modify the effective population size or the selective regime, resulting in altered ratios of non-synonymous to synonymous substitutions among taxa. We studied substitution patterns for the mitochondrial cytochrome oxidase subunit I (COI) gene in a sea star genus (Leptasterias spp.) with an obligate brood-protecting mode of reproduction and small-scale population genetic subdivision, and compared the results to available COI sequences in nine other genera of echinoderms with pelagic larvae: three sea stars, five sea urchins and one brittle star. We predicted that this life history difference would be associated with differences in the ratio of non-synonymous (dN) to synonymous (dS) substitution rates. Leptasterias had a significantly greater dN/dS ratio (both between species and within species), a significantly smaller transition/transversion rate ratio, and a significantly lower average nucleotide diversity within species, than did the non-brooding genera. Other explanations for the results, such as altered mutation rates or selective sweeps, were not supported by the data analysis. These findings highlight the potential influence of reproductive traits and other life history factors on patterns of nucleotide substitution within and between species.     

Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations

We examined rates of DNA sequence evolution in 12 populations of Escherichia coli propagated in a glucose minimal medium for 20,000 generations. Previous work saw mutations mediated by mobile elements in these populations, but the extent of other genomic changes was not investigated. Four of the populations evolved defects in DNA repair and became mutators. Some 500 bp was sequenced in each of 36 genes for 50 clones, including 2 ancestral variants, 2 clones from each population at generation 10,000, and 2 from each at generation 20,000. Ten mutations were found in total, all point mutations including mostly synonymous substitutions and nonsynonymous polymorphisms; all 10 were found in mutator populations. We compared the observed sequence evolution to predictions based on different scenarios. The number of synonymous substitutions is lower than predicted from measured mutation rates in E. coli, but the number is higher than rates based on comparing E. coli and Salmonella genomes. Extrapolating to the entire genome, these data predict about 250 synonymous substitutions on average per mutator population, but only about 3 synonymous substitutions per nonmutator population, during 20,000 generations. These data illustrate the challenge of finding sequence variation among bacterial isolates that share such a recent ancestor. However, this limited variation also provides a useful baseline for research aimed at finding the beneficial substitutions in these populations.     

Distinctively variable sequence-based nuclear DNA markers for multilocus phylogeography of the soybean- and rice-infecting fungal pathogen Rhizoctonia solani AG-1 IA

A series of multilocus sequence-based nuclear DNA markers was developed to infer the phylogeographical history of the Basidiomycetous fungal pathogen Rhizoctonia solani AG-1 IA infecting rice and soybean worldwide. The strategy was based on sequencing of cloned genomic DNA fragments (previously used as RFLP probes) and subsequent screening of fungal isolates to detect single nucleotide polymorphisms (SNPs). Ten primer pairs were designed based on these sequences, which resulted in PCR amplification of 200-320 bp size products and polymorphic sequences in all markers analyzed. By direct sequencing we identified both homokaryon and heterokaryon (i.e. dikaryon) isolates at each marker. Cloning the PCR products effectively estimated the allelic phase from heterokaryotic isolates. Information content varied among markers from 0.5 to 5.9 mutations per 100 bp. Thus, the former RFLP codominant probes were successfully converted into six distinctively variable sequence-based nuclear DNA markers. Rather than discarding low polymorphism loci, the combination of these distinctively variable anonymous nuclear markers would constitute an asset for the unbiased estimate of the phylogeographical parameters such as population sizes and divergent times, providing a more reliable species history that shaped the current population structure of R. solani AG-1 IA.     

Adaptive evolution of G-protein coupled receptor genes

The phylogeny and patterns of nucleotide substitutions in the visual pigment genes, adrenergic receptor genes, muscarinic receptor genes, and in the human mas oncogene were studied by comparing their DNA sequences. The evolutionary tree obtained shows that the visual pigment genes and mas oncogene form one cluster and that the receptor genes form another. In the evolution of rhodopsin genes, synonymous substitutions outnumber nonsynonymous substitutions. This is consistent with the neutral theory of molecular evolution. However, the early evolutionary stages of alpha- and beta-adrenergic and muscarinic receptors are notable for significantly more nonsynonymous substitutions than synonymous substitutions, suggesting the acquisition of novel functional adaptations. Variable rates of nonsynonymous changes in different domains of these proteins reveal DNA segments that might have been important in their functional adaptations.      

Evolutionary trajectory of fish Piscine novirhabdovirus (=Viral Hemorrhagic Septicemia Virus) across its Laurentian Great Lakes history: Spatial and temporal diversification

Piscine novirhabdovirus = Viral Hemorrhagic Septicemia Virus (VHSV) first appeared in the Laurentian Great Lakes with large outbreaks from 2005 to 2006, as a new and novel RNA rhabdovirus subgenogroup (IVb) that killed >30 fish species. Interlude periods punctuated smaller more localized outbreaks in 2007, 2010, and 2017, although some fishes tested positive in the intervals. There have not been reports of outbreaks or positives from 2018, 2019, or 2020. Here, we employ a combined population genetics and phylogenetic approach to evaluate spatial and temporal evolutionary trajectory on its G‐gene sequence variation, in comparison with whole‐genome sequences (11,083 bp) from a subset of 44 individual isolates (including 40 newly sequenced ones). Our results show that IVb (N = 184 individual fish isolates) diversified into 36 G‐gene haplotypes from 2003 to 2017, stemming from two originals (“a” and “b”). G‐gene haplotypes “a” and “b” differed by just one synonymous single‐nucleotide polymorphism (SNP) substitution, remained the most abundant until 2011, then disappeared. Group “a” descendants (14 haplotypes) remained most prevalent in the Upper and Central Great Lakes, with eight (51%) having nonsynonymous substitutions. Group “b” descendants primarily have occurred in the Lower Great Lakes, including 22 haplotypes, of which 15 (68%) contained nonsynonymous changes. Evolutionary patterns of the whole‐genome sequences (which had 34 haplotypes among 44 isolates) appear congruent with those from the G‐gene. Virus populations significantly diverged among the Upper, Central, and Lower Great Lakes, diversifying over time. Spatial divergence was apparent in the overall patterns of nucleotide substitutions, while amino acid changes increased temporally. VHSV‐IVb thus significantly differentiated across its less than two decades in the Great Lakes, accompanied by declining outbreaks and virulence. Continuing diversification likely allowed the virus to persist at low levels in resident fish populations, and may facilitate its potential for further and future spread to new habitats and nonacclimated hosts.     

Molecular evolution of a type 1 wild-vaccine poliovirus recombinant during widespread circulation in China

Type 1 wild-vaccine recombinant polioviruses were isolated from poliomyelitis patients in China from 1991 to 1993. We compared the sequences of 34 recombinant isolates over the 1,353-nucleotide (nt) genomic interval (nt 2480 to 3832) encoding the major capsid protein, VP1, and the protease, 2A. All recombinants had a 367-nt block of sequence (nt 3271 to 3637) derived from the Sabin 1 oral poliovirus vaccine strain spanning the 3'-terminal sequences of VP1 (115 nt) and the 5' half of 2A (252 nt). The remaining VP1 sequences were closely (up to 99.5%) related to those of a major genotype of wild type 1 poliovirus endemic to China up to 1994. In contrast, the non-vaccine-derived sequences at the 3' half of 2A were more distantly related (<90% nucleotide sequence match) to those of other contemporary wild polioviruses from China. The vaccine-derived sequences of the earliest (April 1991) isolates completely matched those of Sabin 1. Later isolates diverged from the early isolates primarily by accumulation of synonymous base substitutions (at a rate of approximately 3.7 x 10(-2) substitutions per synonymous site per year) over the entire VP1-2A interval. Distinct evolutionary lineages were found in different Chinese provinces. From the combined epidemiologic and evolutionary analyses, we propose that the recombinant virus arose during mixed infection of a single individual in northern China in early 1991 and that its progeny spread by multiple independent chains of transmission into some of the most populous areas of China within a year of the initiating infection.