Ribosomal DNA variation in foxtail millet, Setaria italica (L.) P. Beauv., and a survey of variation from Europe and Asia

 Restriction fragment length polymorphism (RFLP) and the structure of ribosomal RNA genes (rDNA) were investigated in 117 landraces of foxtail millet, Setaria italica (L.) P. Beauv. Five RFLP phenotypes were found when the genomic DNA was digested with BamHI; these were named types I–V. Of these types I, II and III were the most frequent. Type I was mainly distributed in the temperature zone, type II in the Taiwan-Philippines Islands and type III in South Asia. Restriction mapping of the cloned rDNA and comparison with RFLP phenotypes showed that the different types originated from a polymorphism in the length within the intergenic spacer (IGS) and BamHI site changes within the IGS. Received: 28 August 1996 / Accepted: 28 February 1997     

Ribosomal DNA Intergenic Spacer of the Swimming Crab, Charybdis japonica

We have determined the full sequence of the ribosomal DNA intergenic spacer (IGS) of the swimming crab, Charybdis japonica, by long PCR for the first time in crustacean decapods. The IGS is 5376 bp long and contains two nonrepetitive regions separated by one long repetitive region, which is composed mainly of four subrepeats (subrepeats I, II, III, and IV). Subrepeat I contains nine copies of a 60-bp repeat unit, in which two similar repeat types (60 bp-a and 60 bp-b) occur alternatively. Subrepeat II consists of nine successive repeat units with a consensus sequence length of 142 bp. Subrepeat III consists of seven copies of another 60-bp repeat unit (60 bp-c) whose sequence is complementary to that of subrepeat I. Immediately downstream of subrepeat III is subrepeat IV, consisting of three copies of a 391-bp repeat unit. Based on comparative analysis among the subrepeats and repeat units, a possible evolutionary process responsible for the formation of the repetitive region is inferred, which involves the duplication of a 60-bp subrepeat unit (60 bp-c) as a prototype. Received: 13 April 1999 / Accepted: 2 August 1999     

Study of genetic diversity and relationships among accessions of foxtail millet [Setaria italica (L.) P. Beauv.] in Korea, China, and Pakistan using SSR markers

In this study, 28 simple sequence repeat (SSR) primer sets were used to analyze the genetic diversity, population structure, and genetic relationships among 37 accessions of foxtail millet from Korea, China and Pakistan. A total of 298 alleles were detected with an average allele number of 10.6 per locus among 37 foxtail millet accessions. The number of alleles per locus ranged from 2 (b226) to 20 (b236). Of the 298 alleles, 138 alleles (46.3%) were rare (frequency < 0.05), 152 alleles (51.0%) were detected at an intermediate frequency (range, 0.05?C0.50), and eight alleles (2.7%) were abundant (frequency > 0.50), respectively. The average gene diversity values were 0.652, 0.692, and 0.491 and polymorphic information content values were 0.621, 0.653, and 0.438, for accessions from Korea, China, and Pakistan, respectively. The accessions from China showed higher SSR diversity than those from Korea and Pakistan. A phylogenetic tree constructed using the un-weighted pair group methods with arithmetic mean algorithm revealed three major groups of accessions that were not congruent with geographical distribution patterns with a few exceptions. The lack of correlation between the accession clusters and their geographic location indicates that the diffusion of foxtail millet from China to Korea might have occurred through multiple routes. Our results provide support for the origin and diffusion route of foxtail millet in East Asia. This SSR-based assessment of genetic diversity, genetic relationships, and population structure among genetic resources of foxtail millet landraces will be valuable to foxtail millet breeding and genetic conservation programs in Korea.     

Structural variation in the Waxy gene and differentiation in foxtail millet [Setaria italica (L.) P. Beauv.]: implications for multiple origins of the waxy phenotype

The origin and evolution of the waxy type of foxtail millet [Setaria italica (L.) P. Beauv] were studied by analyzing structural variation in the Waxy gene. Initially, the Waxy gene was amplified by RT-PCR, RACE and genomic PCR from a non-waxy strain to determine the structure of the wild-type gene. Secondly, we screened by PCR for polymorphisms at the Waxy locus in 79 strains with various waxy phenotypes. We then carried out genomic Southern analysis on 67 strains and identified seven RFLP classes which were designated as types I-VII. RFLP type was correlated with phenotype, such that types I and II corresponded to non-waxy, types III and VI to low-amylose, and types IV, V and VII to waxy phenotypes. The differences between RFLP types could be attributed to insertions in the Waxy gene. Types II and VI were caused by the insertion of a Tourist element into intron 1 and a SINE-like sequence into intron 12, respectively. Types III, IV, V and VII were characterized by the insertion of large sequences into the Waxy gene that may alter the expression of the gene. Thus, multiple, independent insertions in the Waxy gene appear to have caused the loss-of-function waxy phenotypes. Furthermore, the geographical distributions of the three RFLP types associated with the waxy phenotype (types IV, V and VII) were distinct, with type IV being found mainly in Taiwan and Japan, type V in Korea, and type VII in Myanmar. These results indicate a polyphyletic origin for the waxy phenotype in landraces of foxtail millet.     

Mitochondrial DNA analysis of Sporothrix schenckii in North and South America

Mitochondrial DNA (mtDNA) types based on restriction fragment length polymorphism (RFLP) patterns with HaeIII were investigated in clinical isolates of Sporothrix schenckii in North and South America. In addition to 14 mtDNA types (Types 1–14) so far reported, six new mtDNA types, Types 15–20 were found in this study. Type 3 was divided into two subtypes, Subtype 3A and Subtype 3B based on RFLP with Msp1. Type 14 was also divided into three subtypes, Subtype 14A, Subtype 14B and Subtype 14C based on RFLP with Hha1. Nineteen isolates in the United States consisted of 1 isolate of Type 1, 12 of Type 2, 2 of Type 4, 3 of Type 14 (1 of Subtype 14B and 2 of Subtype 14C) and 1 of Type 15. Twenty nine isolates in Venezuela consisted of 13 of Type 3 (Subtype 3B), 6 of Type 4, 1 of Type 18, 3 of Type 19 and 6 of Type 20. Thirteen isolates in Argentina consisted of 2 of Type 3 (Subtype 3A), 4 of Type 4, 4 of Type 16 and 3 of Type 17. One isolate in Brazil was Type 3 (Subtype 3A). Based on the phylogeny of 20 mtDNA types (Types 1–20) constructed by estimating sequence divergences of mtDNA, mtDNA types were clustered into two groups: Group A (Types 1–3, Type 11 and Types 14–19) and Group B (Types 4–10, Types 12–13 and Type 20). These results suggest that S. schenckiiisolates in North and South America mainly belong to Group A. This revised version was published online in June 2006 with corrections to the Cover Date.     

Molecular mechanisms of transformation of SkQ mitotropic quinones and the search for new approaches to creation of selective free radical traps

Features of the mechanism of action of positively charged benzoquinone derivatives (SkQ), which are the analogs of coenzyme Q (I), plastoquinone (II), and tocopherol (III), are discussed. It is usually considered that the main target of these compounds is mitochondria, where they accumulate due to the positive charge of the molecule. In the present work, it is shown with model systems that the reduced forms of compounds (I–III) under certain conditions can transform into electrically neutral cyclic zwitterions, which theoretically can escape from the matrix of energized mitochondria against the concentration gradient. A weak uncoupling effect of molecules I–III has been found on mitochondria. Its existence is in agreement with the abovementioned transformation of positively charged hydroquinones of type Ia–IIIa into electrically neutral molecules. The data obtained with model systems suggest that the target of SkQ hydroquinones as free radical traps may be not only mitochondria but also biochemical systems of the cytoplasm. Due to the presence of a large number of reactive oxygen species (ROS)-dependent signal systems in a cell, the functioning of cytoplasmic systems might be disturbed under the action of antioxidants. The problem of selective effect of antioxidants is discussed in detail in the present work, and a functional diagram of selective decrease of the “background level” of ROS based on differences in the intensity of background and “signal” ROS fluxes is considered.     

Genetic analysis of <Emphasis Type="Italic">NEKODE1</Emphasis> gene involved in panicle branching of foxtail millet, <Emphasis Type="Italic">Setaria italica</Emphasis> (L.) P. Beauv., and mapping by using QTL-seq

We carried out genetic analysis and mapping of a gene for the tip-branched panicle (Nekode or Neko-ashi in Japanese) in foxtail millet. We revealed that this trait is controlled by a single dominant gene by using two F₂ populations and designated the gene as NEKODE1. By using an F₂ population between closely related Taiwanese landraces with a new method based on next-generation sequencing (NGS), QTL-seq, we successfully and rapidly mapped the responsible gene (NEKODE1) on chromosome 9. We also mapped the gene by using SSR markers to verify that this gene is located at the position on chromosome 9, suggested by QTL-seq, and we obtained SSR markers closely linked to the gene and found several candidate genes for this trait in a foxtail millet genome sequence database. The use of a foxtail millet genome sequence and NGS enables rapid mapping of a gene(s) by using a segregation population derived from a cross even between closely related foxtail millet landraces.     

Alteration of mannose transport in fibroblasts from type I carbohydrate deficient glycoprotein syndrome patients

The aim of the present study was to explore how mannose enters fibroblasts derived from a panel of children suffering from different subtypes of type I carbohydrate deficient glycoprotein syndrome: seven carbohydrate deficient glycoprotein syndrome subtype Ia (phosphomannomutase deficiency), two carbohydrate deficient glycoprotein syndrome subtype Ib (phosphomannose isomerase deficiency) and two carbohydrate deficient glycoprotein syndrome subtype Ix (not identified deficiency). We showed that a specific mannose transport system exists in all the cells tested but has different characteristics with respect to carbohydrate deficient glycoprotein syndrome subtypes. Subtype Ia fibroblasts presented a mannose uptake equivalent or higher (maximum 1.6-fold) than control cells with a D-[2-3H]-mannose incorporation in nascent N-glycoproteins decreased up to 7-fold. Compared to control cells, the mannose uptake was greatly stimulated in subtype Ib (4.0-fold), due to lower Kuptake and higher Vmax values. Subtype Ib cells showed an increased incorporation of D-[2-3H]-mannose into nascent N-glycoproteins. Subtype Ix fibroblasts presented an intermediary status with mannose uptake equivalent to the control but with an increased incorporation of D-[2-3H]-mannose in nascent N-glycoproteins. All together, our results demonstrate quantitative and/or qualitative modifications in mannose transport of all carbohydrate deficient glycoprotein syndrome fibroblasts in comparison to control cells, with a relative homogeneity within a considered subtype of carbohydrate deficient glycoprotein syndrome. These results are consistent with the possible use of mannose as a therapeutic agent in carbohydrate deficient glycoprotein syndrome Ib and Ix.     

Characterization of subrepeat regions within rDNA intergenic spacers of the edible basidiomycete Lentinula edodes

For 16 commercial cultivars of Lentinula edodes, DNA fragments for the nuclear rDNA intergenic spacers IGS1 and IGS2 were amplified and analyzed. IGS1 contained a subrepeat region, named SR1, and IGS2 contained a pair of direct repeats and a subrepeat region, named SR2. Three and five types of subrepeats were found in SR1 and SR2, respectively. Heterogeneity in the lengths of IGS1 and IGS2 arose mainly from the number of different kinds of subrepeats within SR1 and SR2. The DNA fingerprints from the PCR products targeting SR1 and SR2 were specific for each of the 16 cultivars, and had enough variation for discrimination among the cultivars. This result suggests that the DNA fingerprints targeting SR1 and SR2 are useful for investigations of L. edodes cultivars.     

Ribosomal DNA Evolution at the Population Level: Nucleotide Variation in Intergenic Spacer Arrays of Daphnia Pulex

Nucleotide variation was surveyed in 21 subrepeat arrays from the ribosomal DNA intergenic spacer of three Daphnia pulex populations. Eighteen of these arrays contained four subrepeats. Contrary to expectations, each of the four positions within the array had a different consensus sequence. However, gene conversion, involving sequences less than the length of a subrepeat, had occurred between subrepeats in different positions. Three arrays had more than four subrepeats and were undoubtedly generated by unequal crossing over between standard-length arrays. The data strongly suggested that most unequal exchanges between arrays are intrachromosomal and that they occur much less frequently than unequal exchanges at the level of the entire rDNA repeat. Strong associations among variants at different positions allowed the recognition of five groups of arrays, two of which were found in more than one population. Five of the seven individuals surveyed had arrays from more than one group. Analysis of the distribution of nucleotide variation suggested that the populations were quite divergent, a result that is concordant with previous surveys of allozyme and mitochondrial DNA variation. It was suggested that some of the subrepeat array types are quite old, at least predating the recolonization of pond habitats in the midwestern United States after the last glaciation.     

禽流感病毒A型和H5亚型RT-PCR检测试剂盒研究

目的　检测和鉴定A型、H5亚型禽流感病毒 (AIV) ,研发一种高效实用的检测手段。方法　根据Ming ShiuhLee报道的文献设计、合成引物 ,采用反转录和PCR一步法对A型、H5亚型禽流感病毒cDNA进行扩增和电泳鉴定 ,组装成禽流感病毒RT PCR试剂盒 ,对H1～ 15亚型AIV参考株、38份AIV国内分离株进行检测试验。结果　建立了A型、H5亚型禽流感病毒RT PCR检测方法 ,并在此基础上组装试剂盒 ,用A型试剂盒检测时 ,全部AIV毒株均为阳性 ,能检测 1 10 2 4血凝单位禽流感病毒 ;用H5亚型试剂盒检测时 ,仅有H5亚型AIV参考株和 19株H5亚型AIV分离株呈阳性 ,其余H1～H4、H6～H15参考株和H7、H9分离株以及 1株H5分株均为阴性 ,能检测1 6 4血凝单位禽流感病毒。 2种试剂盒对实验感染鸡病料检出率均为 10 0 %。结论　研制的AIVA型、H5亚型RT PCR试剂盒具有特异性强、敏感性高、稳定性和重复性好的特点。     

Structural evolution of nrDNA ITS in Pinaceae and its phylogenetic implications

Nuclear ribosomal DNA (nrDNA) has been considered as an important tool for inferring phylogenetic relationships at many taxonomic levels. In comparison with its fast concerted evolution in angiosperms, nrDNA is symbolized by slow concerted evolution and substantial ITS region length variation in gymnosperms, particularly in Pinaceae. Here we studied structure characteristics, including subrepeat composition, size, GC content and secondary structure, of nrDNA ITS regions of all Pinaceae genera. The results showed that the ITS regions of all taxa studied contained subrepeat units, ranging from 2 to 9 in number, and these units could be divided into two types, longer subrepeat (LSR) without the motif (5'-GGCCACCCTAGTC) and shorter subrepeat (SSR) with the motif. Phylogenetic analyses indicate that the homology of some SSRs still can be recognized, providing important informations for the evolutionary history of nrDNA ITS and phylogeny of Pinaceae. In particular, the adjacent tandem SSRs are not more closely related to one another than they are to remote SSRs in some genera, which may imply that multiple structure variations such as recombination have occurred in the ITS1 region of these groups. This study also found that GC content in the ITS1 region is relevant to its sequence length and subrepeat number, and could provide some phylogenetic information, especially supporting the close relationships among Picea, Pinus, and Cathaya. Moreover, several characteristics of the secondary structure of Pinaceae ITS1 were found as follows: (1) the structure is dominated by several extended hairpins; (2) the configuration complexity is positively correlated with subrepeat number; (3) paired subrepeats often partially overlap at the conserved motif (5'-GGCCACCCTAGTC), and form a long stem, while other subrepeats fold onto itself, leaving part of the conserved motif exposed in hairpin loops.     

Nucleotide sequence variation in salmonid alphaviruses from outbreaks of salmon pancreas disease and sleeping disease

We compared 18 salmonid alphaviruses (SAV) including the reference F93-125 salmon pancreas disease virus (SPDV) and S49p sleeping disease virus (SDV) isolates by nucleotide sequence analyses of regions within the E1, nsP4 and nsP3 genes, and found these to comprise 3 distinct groups, which we have designated Subtypes 1, 2 and 3: Subtype 1, which comprised SAVs with sequences closely similar to the reference SPDV isolate, included SAVs from pancreas disease (PD) outbreaks in farmed salmon in Ireland and Scotland over a 10 yr period; viruses from recent outbreaks of sleeping disease (SD) in freshwater-reared trout farmed in England, Scotland and France were closely similar to and were grouped with the reference SDV isolate in Subtype 2; 3 viruses isolated from PD-affected salmon in Norway were genetically different from viruses belonging to Subtypes 1 and 2 and have been assigned to Subtype 3; 1 virus isolated from PD-affected salmon in the Western Isles, Scotland, in 2003 showed consistent nucleotide sequence differences from SAV Subtypes 1, 2 and 3, but was more closely related to the Subtype 1 SAVs. The occurrence of the different subtype SAVs appeared to have a geographical basis, which may prove useful in future molecular epidemiology studies of SAV-induced disease outbreaks.     

Quantitative changes in the expression of histone H1 and H2B subtypes and their relationship to the differentiation of mouse embryonal carcinoma cells

The pattern of histones from several mouse embryonal carcinoma cell (ECC) lines, differentiated cell lines, and adult organs was analyzed using acid-urea gels containing Triton X-100 and long SDS-gel electrophoresis. All cell lines had comparable histone types except for a unique H2B-like component that was found only in the ECC line PCC4. The mouse histone H1 has four different subtypes (H1a, H1b, H1c, and H1d), as resolved in SDS-gel electrophoresis. The expression of the four subtypes was shown to be cell line specific. Subtypes H1a and H1d are present in approximately the same relative amounts in all cell lines investigated. Subtype H1b is found in higher relative amounts than subtype H1c in ECC lines and testis. The ratio of H1b and H1c is reversed in differentiated cell lines and in kidney, white blood cells, liver and spleen. All four subtypes of H1 are phosphorylated although to a different extent in different cell lines. In ECC lines, subtypes H1b and especially H1d incorporate most of a ³²P label, whereas H1c is predominately phosphorylated in differentiated parietal endoderm cell lines. These data indicate that H1 subtypes differ depending on the stage of cell differentiation. Difference in ratio between H1 subtypes and in phosphorylation might influence the chromatin configuration and thus gene expression in these cells.     

临床分离耐甲氧西林溶血性葡萄球菌的SCCmec分型及同源性分析

目的了解临床分离耐甲氧西林溶血性葡萄球菌（MRSH）的SCCmec基因型别及相同SCCmec型别菌株的同源性。方法多重PCR进行SCCmec分型,ERIC-PCR法对相同SCCmec型别菌株进行同源性分析。结果83株临床分离MRSH菌株中,SCCmecI型有23株（27．7％）,SCCmecⅡ型有10株（12．1％）,SCCmecm型有24株（28．9％）,SCCmecIV型有1株（1．2％）,I、Ⅱ混合型有8株（9．6％）,I、Ⅲ混合型有6株（7．2％）,Ⅱ、11混合型有5株（6．0％）,I、Ⅱ、Ⅲ混合型有3株（3．6％）,未分型3株（3．6％）。ERIC—PCR结果显示,23株SCCmecI型分为11型,其中A型5株,B型5株,C型3株,其余8株各为1型,2株未分型;10株SCCmecⅡ型分为6型,其中D型4株,E型2株,3株各为1型,1株未分型;24株SCCmecm型分为9型,其中F型11株,G型2株,H型2株,I型2株,5株各为1型,2株未分型。结论临床分离MRSH中,SCCmecI、Ⅲ型为多,部分菌株呈混合型别;相同SCCmec型别的部分菌株之间可能存在克隆传播。     

The occurrence and identification of <Emphasis Type="Italic">Setaria italica</Emphasis> (L.) P. Beauv. (foxtail millet) grains from the Chengtoushan site (ca. 5800 cal B.P.) in central China,with reference to the domestication centre in Asia

Archaeobotanical remains of Setaria grains and chaff were found at the Chengtoushan site in south-central China (ca. 5800 cal b.p.). Grain shape was determined, using length to breadth ratios, and morphological variation in the upper lemma of modern domesticated and wild Setaria species were examined using scanning electron microscopy as a basis for identifying archaeobotanical remains. Grains of S. viridis, S. yunnanensis, and S.×pycnocoma are slender, whereas S. italica, S. italica var. germinica, S. lutescenes, S. faberi, S. glauca, S. pallidefusca and S. intermedia are round in shape. The papillae distributed on the upper lemma of S. italica are small (8–15 μm) with a non-ridged base, while other Setaria species have large papillae (15–20 μm) with a widely ridged base. The remains of the Setaria from the Neolithic layers at Chengtoushan included S. italica, based on these identification characters. These new finds of foxtail millet are the earliest discoveries from the Yangtze River basin of southern China and are also the earliest evidence for co-cultivation of foxtail millet with rice. The implications of these findings for understanding foxtail millet domestication centres are discussed.     

Pollen morphology of the tribe Rhinantheae (Orobanchaceae) and its systematic significances

Pollen morphology of 36 species representing 14 genera within the tribe Rhinantheae in the family Orobanchaceae was studied and illustrated with light microscopy (LM) and scanning electron microscopy (SEM). Five major pollen types were recognized on the basis of exine ornamentation. Within these major types, minor types (subtypes) were distinguished based on exine surface pattern, size, shape, amb form, colpi and colpus membrane. These types and subtypes are as follows: type I. retipilate: subtype Ia. regular retipilate: (1) pollen size < 27 μm, (2) pollen size > 27 μm, subtype Ib. irregular retipilate; type II. verrucate: subtype IIa. macro-verrucate, subtype IIb. verrucate, subtype IIc. sparse verrucate; type III. retirugulate; type IV. granulate; type V. micro-reticulate. A key to pollen morphology of genera studied within the Rhinantheae was made based on pollen morphology from our study and earlier work. Combining with other sources of information on the Rhinantheae, the systematic relationships of this tribe are discussed. Rhinantheae pollen displays considerable variation between genera and species, with taxonomically significant characters at genus and species level. Palynological characteristics provide evidence for interpreting the conflicting views concerning the “Pterygiella Complex”. The evolutionary trend in exine sculpture of Rhinantheae could be proposed, namely that retipilate sculpturing which is the most widespread type is more primitive than the other types (such as foveolate, granulate, regulate, reticulate, retirugulate and verrucate). The pollen data in present study and the view of Hong (1986), as well as the molecular data from Bennett and Mathews (2006) indicated that Asia and related regions were likely to the origin centre of the tribe Rhinantheae.     

Cloning and analysis of a 6.8-kb rDNA intergenic spacer region of the European ash (Fraxinus excelsior L.)

The 6.8-kb rDNA intergenic spacer region of F. excelsior was isolated from a CsCl/actinomycin-D gradient and cloned into pUC18 for further characterization. We observed the presence of subrepeats delimited by HaeIII enzyme sites. These subrepeats were sub-cloned and 11 clones were sequenced. These corresponded to subrepeated elements of either 32 bp or 41 bp that shared a 23-bp common sequence in the 5 end. Within each family of subrepeats, the percentage of common nucleotides was 84.4% for the 5 32-bp subrepeats and 67.4% for the 640-bp subrepeats. Non-repeated HaeIII fragments of 450 bp and 650 bp were also sub-cloned. To compare homology at the IGS region between the rDNA spacers of F. excelsior and the three related species (F. oxyphylla, F. americana, F. ornus), we conducted Southern hybridization analyses using each member of the 32-bp and 40-bp subrepeat families and the unique 450-bp and 650-bp fragments as probes. These analyses indicated that (1) the American ash is more genetically distant from the other three species that the latter are from each other and (2) F. oxyphylla and F. excelsior are more closely related to each other than to F. ornus.     

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]