An Asian Origin for Subtype IV BK Virus Based on Phylogenetic Analysis

Similarly to other members of the Polyomaviridae family, BK virus (BKV) is thought to have co-evolved with its human host. BKV has four subtypes that are distinguishable by immunological reactivity, with two (subtypes I and IV) being most prevalent in human populations. Subtype I is the major subtype worldwide, whereas subtype IV is prevalent in East Asia and Europe but rare in Africa. The geographic distribution pattern of subtype IV BKV is in apparent disagreement with the hypothesis that BKV co-evolved with humans, since subtype IV rarely occurs in Africa. To elucidate the origin of subtype IV, 53 complete subtype IV sequences derived from East Asians and Europeans were subjected to a detailed phylogenetic analysis using the maximum-likelihood and neighbor-joining methods. We identified six subgroups (a1, a2, b1, b2, c1, and c2) that formed a tree represented by the formula: “(a1, a2), ((b1, b2), (c1, c2)).” Interestingly, we found a close correlation between subtype IV subgroups and population geography; thus, all subgroups except c2 were prevalent in particular East Asian populations, with c2 occurring in both Europe and Northeast Asia. From these findings, we conclude that subtype IV of BKV now prevalent in modern humans is derived from a virus that infected ancestral Asians. We introduce two hypotheses to explain how ancestral Asians became infected with subtype IV BKV. [Reviewing Editor: Dr. Joshua Plotkin] Yuriko Nishimoto and Huai-Ying Zheng are two authors contributed equally to this article.     

Relationships between BK virus lineages and human populations

BK polyomavirus (BKV) is ubiquitous in human populations, infecting children asymptomatically and then persisting in the kidney, in which it can cause nephropathy in renal transplant patients. BKV isolates are classified into four subtypes (I-IV) using serological or genotyping methods, and subtype I is further divided into four subgroups, Ia, Ib-1, Ib-2, and Ic, based on DNA sequence variations. To clarify whether there is an association between BK virus lineages and human populations, we examined BKV-positive urine samples collected from immunocompetent individuals at various locations in Europe, Africa, and Asia. Partial BKV DNA sequences (n=299) in these samples were determined and subjected to phylogenetic and single nucleotide polymorphism analysis to classify BKV isolates around the world. The validity of the classification was confirmed by analyses based on complete BKV DNA sequences. Subtype I was the major subtype throughout the studied regions, and subtype IV was prevalent only in Asia and Europe. Subtype-I subgroups showed close relationships to major geographical areas. It has recently been shown that JC virus (a human polyomavirus closely related to BKV) co-evolved with human populations, and the present study thus suggests that host-linked evolution is the general mode of polyomavirus evolution. Additionally, our results indicate certain unique aspects of the relationship between BKV and humans.     

Evolution of BK Virus Based on Complete Genome Data

The human polyomavirus BK virus (BKV) is ubiquitous in humans, infecting children asymptomatically. BKV is the only primate polyomavirus that has subtypes (I–IV) distinguishable by immunological reactivity. Nucleotide (nt) variations in a major capsid protein (VP1) gene region (designated the epitope region), probably responsible for antigenic diversity, have been used to classify BKV isolates into subtypes. Here, with all the protein-encoding gene sequences, we attempted to elucidate the evolutionary relationships among 28 BKV isolates belonging to subtypes I, III, and IV (no isolate belonging to subtype II, a minor one, was included). First, using the GTR + Γ + I model, maximum likelihood trees were reconstructed for individual viral genes as well as for concatenated viral genes. On the resultant trees, the 28 BKV isolates were consistently divided into three clades corresponding to subtypes I, III, and IV, although bootstrap probabilities are not always high. Then we used more sophisticated likelihood models, one of which takes account of codon structure, to elucidate the phylogenetic relationships among BKV subtypes, but the phylogeny of the deep branchings remained ambiguous. Furthermore, the possibility of positive selection in the evolution of BKV was examined using the nonsynonymous/synonymous rate ratio as a measure of selection. An analysis based on entire genes could not detect any strong evidence for positive selection, but that based on the epitope region identified a few sites potentially under positive selection (these sites were among those showing subtype linked polymorphisms). These author Yuriko Nishimoto and Tomokazu Takasaka contributed equally to this article. [Reviewing Editor: Dr. Rasmus Nielsen]     

Two hidden taxa in the Japanese encephalitis vector mosquito,Culex tritaeniorhynchus,and the potential for long-distance migration from overseas to Japan

The Culex vishnui subgroups, particularly Culex tritaeniorhynchus, are considered the primary vectors of the Japanese encephalitis virus (JEV) in Asia. Recent molecular phylogenetic analyses of JEV isolates from Asian countries have shown that JEVs with diverse genetic variants are present in Asia. Furthermore, some JEV strains have been found to have crossed the East China Sea and been introduced into Japan. In this study, the possibility of overseas migration of the JE vector mosquito, Cx. tritaeniorhynchus was examined from the genetic, physical, and meteorological perspectives. Molecular phylogenetic analysis was performed based on both whole coding sequences and on the barcoding region of the mitochondrial cytochrome c oxidase subunit I (COI) gene of Cx. vishnui subgroups collected from Asian countries. Culex tritaeniorhymchus was classified into two genetically independent taxa by COI sequences: the Japanese type (Ct-J), which inhabits Japan except for the Amami Islands of southern Japan, and the continental type (Ct-C), which inhabits the Asian region except for Japan. It was confirmed that approximately 10% of Cx. tritaeniorhynchus trapped during the summer in western Kyushu were Ct-C, and that they could fly for up to 38 h continuously. The meteorological analysis also confirmed that the atmospheric flow occurring over the continent coincided with the date of Ct-C capture. This is the first report showing the existence of two taxa in Cx. tritaeniorhynchus. Their physical and physiological characteristics suggest the possibility of long-distance migration from overseas regions to Japan across the East China Sea. Future efforts are expected to provide evidence to support the occurrence of long-distance migration of Cx. tritaeniorhynchus with JEV.     

JC virus strains indigenous to northeastern Siberians and Canadian Inuits are unique but evolutionally related to those distributed throughout Europe and Mediterranean areas

Human polyomavirus JC virus (JCV) isolates around the world are classified into more than 10 geographically distinct genotypes (designated as subtypes). Evolutionary relationships among JCV subtypes were recently examined, and the following pattern of JCV evolution was indicated. The ancestral JCV first divided into three superclusters, designated Types A, B, and C. A split in Type A generated two subtypes, EU-a and -b, containing mainly European and Mediterranean isolates. The split in Type B generated Af 2 (the major African subtype), Bl-c (a minor European subtype), and various Asian subtypes. Type C generated a single subtype (Afl), consisting of isolates derived from western Africa. In this study, JCV isolates prevalent among northeastern Siberians and Canadian Inuits were evaluated in the context of the above-described pattern of JCV evolution. The Siberian/Arctic JCV isolates were classified as belonging mainly to Type A, based on the result of a preliminary phylogenetic analysis. We then examined, using the whole-genome approach, the phylogenetic relationships among worldwide Type A isolates. In neighbor-joining and maximum-likelihood analyses, Type A JCVs worldwide consistently diverged into three subtypes, EU-a, -b, and -c, with high bootstrap probabilities. EU-c was constructed only by northeastern Siberian isolates, derived mainly from Nanais living in the lower Amur River region, and was shown to have been generated by the first split in Type A. Most Siberian/Arctic isolates derived from Chukchis, Koryaks, and Canadian Inuits formed a distinct cluster within the EU-a subtype, with a high bootstrap probability. Based on the present findings, we discuss ancient human migrations, accompanied by Type A JCVs, across Asia and to Arctic areas of North America.     

Phylogenetic analysis of the rDNA intergenic spacer subrepeats and its implication for the domestication history of foxtail millet, Setaria italica

We sequenced ribosomal DNA intergenic spacer subrepeats and their flanking regions of foxtail millet landraces from various regions in Europe and Asia and its wild ancestor to elucidate phylogenetic differentiation within each of types I–III found in our previous work and to elucidate relationships among these three types. Type I was classified into seven subtypes designated as Ia–Ig based on subrepeat sequences; C repeats downstream of those subrepeats are also polymorphic. Of these, subtypes Ia–Id and Ig were found in foxtail millet landraces. Subtypes Ia and Ib were distributed broadly throughout Asia and Europe. Subtype Ic was distributed in China, Korea and Japan. Subtype Id has a 20-bp deletion in subrepeat 3 and has a unique C repeat sequence. This subtype was found in a morphologically primitive landrace group from Afghanistan and northwestern Pakistan and differed greatly from other type I subtypes, implying that these landraces were domesticated independently. Subtypes Ig was found in a landrace from Pakistan and Ia and Ie–Ig were in six wild ancestor accessions. Type II was also highly polymorphic and four subtypes were found and designated as subtypes IIa–IId, but sequence analyses indicated type III as monomorphic. The present work indicates that type III should be classified as a subtype of type II (subtype IIe). Sequence polymorphism of subrepeats of types I–III indicated that subrepeats of subtype IIa are greatly divergent from others. Relationships among types I–III are much more complicated than anticipated based on previous RFLP work.     

Occurrence of the African subgroup (Ia) of BK polyomavirus in younger Japanese children

BK polyomavirus (BKV) is ubiquitous among humans, usually infecting them asymptomatically during childhood. BKV persists in renal tissue of individuals and their progeny are excreted in urine, particularly in immunocompromised patients. JC virus, another human polyomavirus, has been considered to be transmitted from parents to children during prolonged cohabitation. However, BKV has been supposed to be transmitted not only within but also outside the family. In the present study, to clarify this possibility, we analyzed phylogenetically 35 BKV which were excreted in the urine by Japanese children and adults undergoing stem cell transplantation. Subtypes I, III and IV were detected in 15, two and one children and in 15, one and one adults, respectively. Among 15 subtype I isolates from children, three, four and eight belonged to subgroups Ia, Ib-1 and Ic, respectively. All the three children from whom Ia was detected were less than 9 years old. In contrast in the adults, three subtype I belonged to Ib-1 and the other 12 to Ic. These findings may reflect the recent transmission of BKV Ia strains to Japanese children.     

Two unusual human immunoglobulin V kappa genes

The V kappa genes A10 and A14 which have been previously localized within the human kappa locus were analysed now. A10 hybridizes under stringent conditions only weakly or not at all to probes characteristic for the four V kappa subgroups. According to their DNA sequences and the derived amino-acid sequences A10 and A14 do not fit well into the subgroup classification. They seem to be about as closely related to the subgroup I and III genes and less related to those of subgroups II and IV. Hybridization experiments indicate that A10 and A14 belong to a small V kappa gene family. After discussing the various features of the sequences we suggest neither to assign A10 and A14 to one of the existing subgroups nor to establish a new one but to apply to them the subgroup designation N which may be changed when all V kappa genes are known and can be classified together.     

A meta-analysis of interleukin-10-1082 promoter polymorphism associated with gastric cancer risk

We aimed to explore the role of allele A/G single nucleotide polymorphism (SNP) of gene Interleukin 10 (IL-10) promoter-1082 in the susceptibility to gastric cancer through a systematic review and meta-analysis. Each initially included article was scored for quality appraisal. Desirable data were extracted and registered into databases. Twenty studies were ultimately eligible for the meta-analysis of IL-10-1082 A/G SNP. We adopted the most probably appropriate genetic model (dominant model), with the combined group of GG-plus-GA genotypes compared with the AA genotype. Potential sources of heterogeneity were sought out via subgroup analyses and sensitivity analyses, and publication biases were estimated. Between IL-10-1082 GG-plus-GA genotypes with the risk of developing gastric cancer, statistically significant association could be noted with overall gastric cancer, being mainly in Asian subgroup, large sample subgroup, high quality subgroup, intestinal-type subgroup, cardia-type subgroup, and some genotyping method subgroups. Our meta-analysis indicates that IL-10-1082 GG-plus-GA genotypes are associated with the overall risk of developing gastric cancer and seem to be more susceptible to overall gastric cancer in Asian populations. IL-10-1082 GG-plus-GA genotypes are more associated with the pathologically intestinal-type gastric cancer or anatomically cardia-type gastric cancer.     

Genetic differentiation of the gobies Gymnogobius castaneus and G. taranetzi in the region surrounding the sea of Japan as inferred from a mitochondrial gene genealogy

The phylogenetic relationships between gobies of the genus Gymnogobius were analyzed using mitochondrial cytochrome b gene sequences, focusing on the species currently classified as G. taranetzi and G. castaneus that occur in Japan, South Korea, and Russia. Gobies of the two species collected at 12 localities in Japan, South Korea, and Russia formed a monophyletic clade (called the "castaneus species complex" here) with G. breunigii as the sister clade. Within the species complex, six lineages were recognized: (L1) G. castaneus from the Akigawa River, Tokyo, Japan; (L2) G. castaneus from Yuza, Yamagata, Japan; (L3) G. taranetzi from Russia and South Korea; (L4) G. castaneus from the Tonegawa River, Chiba, Japan; (L5a) G. taranetzi from Shimane, Japan; and (L5b) G. castaneus + G. taranetzi from the Japan Sea coast of northern Japan. The two local lineages of G. castaneus (L1 and L2) are highly divergent from the others. The Japanese populations of G. taranetzi have diverged from the continental G. taranetzi populations, while one mitochondrial lineage (L5b) is shared with G. castaneus of northeast Japan. Therefore, the current species G. taranetzi and G. castaneus as defined morphologically are polyphyletic, necessitating a taxonomic revision. The genetic differentiation of isolated local lineages and the evolution of taranetzi-and castaneus-type gobies have likely occurred repeatedly in brackish/freshwater habitats around the Sea of Japan. We discussed the time of divergence for these gobies based on a tree with the molecular clock assumption.     

Characterization and phylogenetic analysis of Varicella-zoster virus strains isolated from Korean patients

Varicella-zoster virus (VZV) is a causative agent of chickenpox in primary infection and shingles after its reactivation from latency. Complete or almost-complete genomic DNA sequences for various VZV strains have been reported. Recently, clinical VZV strains were isolated from Korean patients whose genome was sequenced using high-throughput sequencing technology. In this study, we analyzed single nucleotide polymorphism (SNP) of VZV strains to genetically characterize Korean clinical isolates. Phylogenetic analyses revealed that three Korean strains, YC01, YC02, and YC03, were linked to clade 2. Comprehensive SNP analysis identified 86 sites specific for the 5 VZV clades. VZV strains isolated from Korea did not form a phylogenetic cluster. Rather, YC02 and YC03 clustered strongly with Chinese strain 84-7 within clade 2, more specifically cluster 2a. Signature sequences for the cluster 2a were identified and found to play an important role in the separation of cluster 2a strains from other clade 2 strains, as shown in substitution studies. Further genetic analysis with additional strains isolated from Japan, China, and other Asian countries would provide a novel insight into the significance of two distinct subclades within clade 2.     

The variable region of human immunoglobulins. I. Serologic and structural correlations of antigenic markers common to V lambda I and V lambda IV proteins (isotypic cross-reactivities)

Sixty-eight human lambda-chain sources including representatives of the five different subgroups were studied for the distribution of the serologic markers ST, 111, and VOR The variable region character of these antigenic determinants has been demonstrated for those proteins whose primary structures had been partially or completely determined previously. The lambda-chain sub-groups I and IV are readily distinguished antigenically but the specificity designated VOR is shared between them and is absent in the II, III, and V subtypes. The results of these studies are discussed with respect to the possible relationships of the structural genes controlling the synthesis of the lambda-chain subgroups, and some potential phylogenetic and functional meaning of the heterogeneities is revealed.     

Y-chromosomal DNA haplogroups and their implications for the dual origins of the Koreans

We have analyzed eight Y-chromosomal binary markers (YAP, RPS4Y₇₁₁, M9, M175, LINE1, SRY₊₄₆₅, 47z, and M95) and three Y-STR markers (DYS390, DYS391, and DYS393) in 738 males from 11 ethnic groups in east Asia in order to study the male lineage history of Korea. Haplogroup DE-YAP was found at a high frequency only in Japan but was also present at low frequencies in northeast Asia, including 2.5% in Korea, suggesting a northern origin for these chromosomes. Haplogroup C-RPS4Y₇₁₁ was present in Korea and Manchuria at moderate frequencies: higher than in populations from southeast Asia, but lower than those in the northeast, which may imply a northern Asian expansion of these lineages, perhaps from Mongolia or Siberia. The major Y-chromosomal expansions in east Asia were those of haplogroup O-M175 (and its sublineages). This haplogroup is likely to have originated in southern east Asia and subsequently expanded to all of east Asia. The moderate frequency of one sublineage in the Koreans, haplogroup O-LINE1 (12.5%), could be a result of interaction with Chinese populations. The age of another sublineage, haplogroup O-SRY₊₄₆₅, and Y-STR haplotype diversity provide evidence for relatively recent male migration, originally from China, through Korea into Japan. In conclusion, the distribution pattern of Y-chromosomal haplogroups reveals the complex origin of the Koreans, resulting from genetic contributions involving the northern Asian settlement and range expansions mostly from southern-to-northern China.     

PanSNPdb: the Pan-Asian SNP genotyping database

The HUGO Pan-Asian SNP consortium conducted the largest survey to date of human genetic diversity among Asians by sampling 1,719 unrelated individuals among 71 populations from China, India, Indonesia, Japan, Malaysia, the Philippines, Singapore, South Korea, Taiwan, and Thailand. We have constructed a database (PanSNPdb), which contains these data and various new analyses of them. PanSNPdb is a research resource in the analysis of the population structure of Asian peoples, including linkage disequilibrium patterns, haplotype distributions, and copy number variations. Furthermore, PanSNPdb provides an interactive comparison with other SNP and CNV databases, including HapMap3, JSNP, dbSNP and DGV and thus provides a comprehensive resource of human genetic diversity. The information is accessible via a widely accepted graphical interface used in many genetic variation databases. Unrestricted access to PanSNPdb and any associated files is available at: http://www4a.biotec.or.th/PASNP.     

The robust phylogeny of Korean wild boar (<Emphasis Type="Italic">Sus scrofa coreanus</Emphasis>) using partial D-loop sequence of mtDNA

In order to elucidate the precise phylogenetic relationships of Korean wild boar (Sus scrofa coreanus), a partial mtDNA D-loop region (1,274 bp, NC_000845 nucleotide positions 16576-1236) was sequenced among 56 Korean wild boars. In total, 25 haplotypes were identified and classified into four distinct subgroups (K1 to K4) based on Bayesian phylogenetic analysis using Markov chain Monte Carlo methods. An extended analysis, adding 139 wild boars sampled worldwide, confirmed that Korean wild boars clearly belong to the Asian wild boar cluster. Unexpectedly, the Myanmarese/Thai wild boar population was detected on the same branch as Korean wild boar subgroups K3 and K4. A parsimonious median-joining network analysis including all Asian wild boar haplotypes again revealed four maternal lineages of Korean wild boars, which corresponded to the four Korean wild boar subgroups identified previously. In an additional analysis, we supplemented the Asian wild boar network with 34 Korean and Chinese domestic pig haplotypes. We found only one haplotype, C31, that was shared by Chinese wild, Chinese domestic and Korean domestic pigs. In contrast to our expectation that Korean wild boars contributed to the gene pool of Korean native pigs, these data clearly suggest that Korean native pigs would be introduced from China after domestication from Chinese wild boars.     

The genetic diversity of the Vigna angularis complex in Asia.

A selected set of accessions of components of the azuki bean (Vigna angularis) complex comprising 123 cultivated accessions and 23 wild or weedy accessions from Bhutan, China (including Taiwan), India, Japan, Korea, and Nepal was analyzed using amplified fragment length polymorphism (AFLP) methodology. Using 12 AFLP primer pairs, 580 unambiguous bands were generated, 313 (53.9%) of which were polymorphic among azuki bean accessions. All 580 bands were used to assess phenotypic (band) and genetic (nucleotide) diversity among the 146 azuki bean accessions. The results indicate five major groups of azuki bean germplasm primarily associated with geographic origin of accessions and their status: wild, weedy, or cultivated. These five groups are (i) Himalayan wild, (ii) Nepal-Bhutan cultivated, (iii) Chinese wild, (iv) Taiwan wild - Bhutan cultivated, and (v) northeast Asian accessions. Within the northeast Asian accessions, three subgroups are present. These consist of (v1) Japanese complex - Korean cultivated, (v2) Japanese cultivated, and (v3) Chinese cultivated accessions. The results suggest domestication of azuki bean occurred at least twice, once in the Himalayan region of southern Asia and once in northeast Asia. The remarkable diversity of azuki bean germplasm in the Himalayan region compared with other regions suggests this is a rich source of germplasm for plant breeding. The results suggest there are important gaps in the germplasm collections of azuki bean and its close relatives from various parts of Asia and that specific collecting missions for Vigna germplasm related to azuki bean in the highlands of subtropical Asia are needed.     

东北地区中鼩鼱遗传多样性与分子系统地理学

本研究对我国东北地区3个种群（大兴安岭、小兴安岭和长白山山脉）的中鼩鼱（99个样本）Cyt b基因全序列进行分析,共获得81个单倍型。整体单倍型多态性为0.9940,显示东北地区中鼩鼱整体遗传多样性较高,3个种群的遗传多样性无明显差异。分子变异方差分析显示3个种群之间的分化占总分化的2.48%,种群内采样地的分化占总分化的0.34%,采样地内部遗传分化占总分化的97.19%。遗传分化结果显示中鼩鼱遗传分化与地理距离存在非正相关关系。结合GenBank下载欧亚其他地区中鼩鼱Cyt b基因全序列构建系统发生树,分化为3个分支：日本北海道分布的中鼩鼱构成1个小分支;韩国济州岛的中鼩鼱构成1个小分支;我国东北地区和其他欧亚大陆的中鼩鼱构成1个大分支,在这个大支上未显现进一步的遗传分化结构。中介网络分析也暗示除日本北海道和韩国济州岛外,欧亚大陆分布的中鼩鼱没有明显的地理系统分化,表现出遗传分化与地理关联性较弱。种群历史分析显示东北地区中鼩鼱经历过种群扩张。     

Phylogenetic analysis of Orientia tsutsugamushi strains based on the sequence homologies of 56-kDa type-specific antigen genes

Close and distant relationship among 31 strains of Orientia tsutsugamushi (20, two, one and eight strains were isolated in Japan, Korea, China and southeast Asia, respectively) were clarified using phylogenetic analyses based on homologies of 56-kDa type-specific antigen genes. Isolates in Japan, Korea and China were located in eight separate clusters in the phylogenetic tree, and each was designated as JG (Japanese Gilliam type), JP-1 and JP-2 (Japanese Karp 1 and 2 types), Kato, Kawasaki, Kuroki, Shimokoshi and LX-1 types. All isolates originated in southeast Asia, including the prototype Gilliam and Karp strains isolated in Burma and New Guinea, respectively, were distantly located in the phylogenetic tree from those isolates in Japan, Korea and China, indicating that strains of O. tsutsugamushi distributed in northeastern and southeastern Asia are different types.     

Gene Flow between the Korean Peninsula and Its Neighboring Countries

SNP markers provide the primary data for population structure analysis. In this study, we employed whole-genome autosomal SNPs as a marker set (54,836 SNP markers) and tested their possible effects on genetic ancestry using 320 subjects covering 24 regional groups including Northern ( = 16) and Southern ( = 3) Asians, Amerindians ( = 1), and four HapMap populations (YRI, CEU, JPT, and CHB). Additionally, we evaluated the effectiveness and robustness of 50K autosomal SNPs with various clustering methods, along with their dependencies on recombination hotspots (RH), linkage disequilibrium (LD), missing calls and regional specific markers. The RH- and LD-free multi-dimensional scaling (MDS) method showed a broad picture of human migration from Africa to North-East Asia on our genome map, supporting results from previous haploid DNA studies. Of the Asian groups, the East Asian group showed greater differentiation than the Northern and Southern Asian groups with respect to Fst statistics. By extension, the analysis of monomorphic markers implied that nine out of ten historical regions in South Korea, and Tokyo in Japan, showed signs of genetic drift caused by the later settlement of East Asia (South Korea, Japan and China), while Gyeongju in South East Korea showed signs of the earliest settlement in East Asia. In the genome map, the gene flow to the Korean Peninsula from its neighboring countries indicated that some genetic signals from Northern populations such as the Siberians and Mongolians still remain in the South East and West regions, while few signals remain from the early Southern lineages.