J Gen Virol Tips for Better Browsing
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 89 (2008), 1446-1456; DOI 10.1099/vir.0.83607-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table and Figures
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.
Agricola
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.

Complete sequence characterization of isolates of Getah virus (genus Alphavirus, family Togaviridae) from China

You-gang Zhai1, Huan-Yu Wang1, Xiao-hong Sun1, Shi-hong Fu1, Huan-qin Wang1, Houssam Attoui2, Qing Tang1 and Guo-dong Liang1

1 State Key Laboratory for Infectious Disease Prevention and Control, Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, 100 Ying Xin St, Beijing 100052, PR China
2 Department of Arbovirology, Institute for Animal Health, Pirbright, Woking, Surrey GU24 0NF, UK

Correspondence
Guo-dong Liang
gdliang{at}hotmail.com


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Ten virus isolates belonging to species Getah virus (GETV) have been obtained during surveys for arboviruses in China since 1964. Seven of these isolates (YN0540, YN0542, SH05-6, SH05-15, SH05-16, SH05-17 and GS10-2) were obtained during the current study. The full-length sequences of three Chinese isolates (M1, isolated in 1964; HB0234, isolated in 2002; YN0540, isolated in 2005) were determined. The full-length sequences of these isolates were respectively 11 696, 11 686 and 11 690 nt, and showed more than 97 % intraspecies identity. Deletions were found in the capsid protein of strain M1 and non-structural protein nsP3 of strain HB0234. The E2 gene and 3' UTR of all ten isolates were also characterized. The E2 gene of the Chinese GETV isolates showed nucleotide sequence identities of 98–100 % when compared with other GETV isolates. In the 3' UTR of the Chinese isolates, an insertion of 10 consecutive adenine residues (nt 189–198) appeared in strain M1, and 9 or 3 consecutive adenines were found towards the 3' end of the third RES in strains SH05-6 and SH05-15, respectively. The 3' UTRs of the Chinese isolates showed a deletion between positions 45 and 54 and nucleotide transitions at positions 43, 64 and 148. Sequence and phylogenetic analyses showed that there was a relatively high degree of conservation among GETV isolates. The isolation of GETV from various provinces in China and also in Russia and Mongolia (including regions of the northern tundra) are an indication of changes in the world distribution of this re-emerging virus.

The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are EU015061–EU015076.

Details of primers used for RT-PCR and sequencing and alignments of GETV and SAGV E2 and 3' UTR sequences are available as supplementary material with the online version of this paper.


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The genus Alphavirus belongs to the family Togaviridae and includes 29 virus species that are recognized by the International Committee for the Taxonomy of Viruses (Weaver et al., 2005Down). Alphaviruses are important pathogens of livestock animals and humans worldwide. The alphavirus particle contains a single copy of the virus genome, which is a linear positive-sense single-stranded RNA. The genome length of the alphaviruses that have been sequenced ranges from 11 442 nt for Semliki Forest virus to 11 851 nt for Ross River virus (Faragher et al., 1988Down; Strauss & Strauss, 1994Down; Pfeffer et al., 1998Down). The genome has a 5' terminal cap and the 3' end has a poly(A) tract. The 5' two-thirds of the genome encodes the non-structural proteins nsP1 to nsP4, followed by the 26S junction region, which promotes transcription of the intracellular subgenomic 26S RNA. This latter mRNA contains the structural gene that encodes the C (capsid), E3, E2, 6K and E1 proteins and represents the remaining 3' third of the genome (Pfeffer et al., 1998Down).

Most alphaviruses are transmitted by haematophagous arthropods and can be grouped into eight antigenic complexes based on their serological cross-reactions. These are the Barmah Forest, Eastern equine encephalitis, Middleburg, Ndumu, Semliki forest, Trocara, Venezuelan equine encephalitis and Western equine encephalitis complexes (Griffin, 2001Down; Weaver et al., 2005Down). Antigenic cross-reactivities between alphaviruses reflect sequence conservation in the C protein and E1 glycoprotein, while antibodies directed against the E2 protein are usually virus-specific (Griffin, 2001Down).

The species Getah virus (GETV) is a member of the Semliki Forest (SFV) complex (Weaver et al., 2005Down), based on genetic and antigenic properties. GETV was first isolated from Culex gelidus in Malaysia in 1955 (Berge, 1975Down); the virus was isolated in rubber plantations, hence the name Getah, which means rubber in Malay.

To date seven Alphavirus species, namely GETV, Sindbis virus (SINV), Ross River virus (RRV), Chikungunya virus (CHIKV), Mayaro virus (MAYV) Eastern equine encephalitis virus (EEEV) and Western equine encephalitis virus (WEEV), have been isolated from various geographical locations in China (Zhang et al., 1989Down; Li et al., 1992Down; Zhao et al., 1997Down; Liang et al., 2000Down; Huang et al., 2001Down; Lv et al., 2001Down; Wang & Liang, 2003Down; Wang et al., 2006Down). Chinese alphaviruses may show significant sequence divergence from their relatives isolated elsewhere; for example, isolate XJ-160 of SINV from China exhibited 18 % overall nucleotide (8.6 % amino acid) variation from the prototype SINV (Liang et al., 2000Down).

The closest phylogenetic relatives of GETV within the SFV complex are RRV and Sagiyama virus (SAGV) (Calisher et al., 1980Down; Powers et al., 2001Down; Attoui et al., 2007Down). RRV is a known human pathogen and SAGV was isolated from various mammals (Attoui et al., 2007Down). Previously published data show that the taxonomic status of SAGV should be reassessed (Shirako & Yamaguchi, 2000Down; Attoui et al., 2007Down). Sequence analysis suggested that SAGV should be reconsidered as a member of the species GETV.

GETV is widely distributed in the countries of South-east Asia and in northern Australia along the Pacific Ocean (Norder et al., 1996Down; Turell et al., 2003Down; Bryant et al., 2005Down; Wang et al., 2006Down; Chang et al., 2006Down). The virus is a known pathogen of horses and pigs. In horses, the disease is characterized by pyrexia, rash, oedema of the hind legs and intumescence of lymph nodes. The infection in pigs results in miscarriage (Kamada et al., 1980Down; Kono et al., 1980Down; Kumanomido et al., 1988aDown, bDown; Shibata et al., 1991Down; Brown & Timoney, 1998Down). Neutralizing antibodies to GETV have been identified in human sera and in birds in Malaysia, northern Australia and Hainan province in China (Kanamitsu et al., 1979Down; Marchette et al., 1980Down; Li et al., 1992Down). Despite these findings, the virus has not been linked to illness in humans.

In 1964, an alphavirus designated M1 was isolated in Hainan province (Li et al., 1992Down) and was subsequently identified as an isolate of GETV by serological analyses. During investigations of arboviruses in China in recent years, several new GETV isolates have been obtained from mosquitoes. In 2002, 38 years after the first isolation, two more GETV strains were isolated, from wild-caught mosquitoes collected in Hebei province in northern China (Wang et al., 2006Down). This isolation reflects the first known incursion of GETV into the northern regions of China.

In total, 10 GETV strains have been isolated in China, including the M1 strain from Hainan, two strains from Hebei and seven strains (obtained during this study) from Yunnan, Ganshu and Shanghai (see Fig. 1Down). The geographical distribution of these isolates ranged between latitudes 1 ° and 4 ° N and longitudes 9 ° and 12 ° E, reflecting the wide distribution of GETV in China. The sequence of the M1 strain has been reported previously (Zhao et al., 2000Down; Wen et al., 2007Down). In this study, we report the full-length sequence characterization of various GETV isolates from China, with particular emphasis on the genetic characterization of the envelope genes and 3' untranslated regions (UTRs). The sequence characterization suggests that GETV has spread recently to China and that the virus has been evolving slowly. Arbovirus surveillance studies carried out before 2000 in China failed to identify the circulation of GETV in this country.


Figure 1
View larger version (30K):
[in this window]
[in a new window]

 
Fig. 1. Geographical locations from which the Chinese GETV isolates were obtained. Shaded areas indicate provinces from which isolates were obtained as follows: 1964, one strain from Hainan; 2002, two strains from Hebei; 2005 six strains (two strains from Yunnan, four strains from Shanghai); 2006, one strain from Gansu.

 

   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Mosquito collection, treatment and virus isolation.
Mosquitoes were collected from villages where livestock were bred in Shanghai and Yunnan province in 2005 and in Gansu province in 2006. Collection locations were within 10 m of henhouses, hog pens and sheep pens. Collected mosquitoes were frozen for 30 min at –20 °C; they were then placed on an ice plate to determine mosquito species and to exclude blood-fed or male mosquitoes. Fifty to 100 mosquitoes were sorted into a collecting tube and stored in liquid nitrogen. Pooled mosquitoes were added to 2 ml minimal essential medium (HyClone) supplemented with 2 mM glutamine, 0.12 % NaHCO3, 100 U penicillin ml–1 and 100 U streptomycin ml–1, followed by grinding in a pre-cooled sterile plastic grinding tube. The ground samples were centrifuged at 12 000 r.p.m. in a microcentrifuge for 20 min at 4 °C and the clarified supernatant was assayed for virus by inoculation onto C6/36 cells and BHK-21 cells. The development of cytopathic effects (CPE) was followed daily. Identification of virus isolates was accomplished by antigen ELISA, using mouse-derived specific antibodies to alphavirus [including Mayaro virus (MAYV), SAGV, RRV, CHIKV, EEEV, WEEV, GETV, SINV and SFV] and by RT-PCR using generic alphavirus primers M2W (5'-YAgAgCDTTTTCgCAYSTRgCHW-3', W=A or T; Y=C or T; R=A or G; H=A, C or T; D=A, G or T; S=C or G), M2W2 (5'-TGYCCNVTGMDNWSYVCNGARGAYCC-3', V=A, C or G; M=A or C, N=A, T, C or G) and cM3W (5'-ACATDAANKgNgTNgTRTCRAANCCDAYCC-3', K=G or T) (Pfeffer et al., 1997Down), followed by sequencing.

Virus propagation.
Ten viral isolates designated M1 (isolated in 1964 in Hainan province; Li et al., 1992Down), HB0215-3, HB0234 (isolated in 2002 in Hebei province; Wang et al., 2006Down), YN0540, YN0542, SH05-6, SH05-15, SH05-16, SH05-17 and GS10-2 (isolated during this study) were propagated in BHK-21 cells at 37 °C or Aedes albopictus C6/36 mosquito cells at 28 °C in minimal essential medium (HyClone) supplemented with 2 % fetal bovine serum (FBS), 2 mM glutamine, 0.12 % NaHCO3, 100 U penicillin ml–1 and 100 U streptomycin ml–1.

Preparation of viral RNA, reverse transcription and PCR amplification of the GETV genome.
Viral RNA was extracted from 140 µl supernatant of virus-infected BHK-21 cell cultures using a viral RNA kit (Qiagen) according to the manufacturer's instruction. cDNA was synthesized using You-Prime-Ready-To-Go beads (Amersham Bioscience) and random hexanucleotide primers. Forty-three primers were designed from the complete sequences of GETV (GenBank accession number AY702913 [GenBank] ; South Korea) and SAGV (AB032553 [GenBank] ; Japan) genomes. These primers were used for full-length amplification and sequencing of Chinese GETV isolates. The list of primers is provided in Supplementary Table S1, available in JGV Online.

The 5' UTR of the viral RNA was amplified using the 5'-Full RACE Core Set (Takara) and the 3' UTR was amplified using primers Ma05F4 (5'-CGGCAATGACATGGGTGCAGC-3', positions 11185–11205), Get3'UR (5'-CTGTAAAATATTAAAAAAACAAATTAGACGC-3', 11670–11698) and M06R1 [5'-CTGTCAGCGAATTCGGTACC(T)18G-3', 11698]. PCR was performed using Ex Taq DNA polymerase (Takara). Amplified products were examined by agarose gel electrophoresis (1 % gel) and purified by QIA quick Gel Extraction kit (Qiagen) and then sequenced directly. At the same time, the amplified products of 5'-RACE and the 3' UTR were cloned into pGEM-T easy vector (Promega) and sequenced using M13 universal primers. Each nucleotide position was sequenced three times. The contig sequences were assembled using ContigExpress Project software (Vector NTI Advance 1001 package; Invitrogen).

Sequence analysis and phylogenetic comparisons.
Nucleotide and amino acid sequence alignments were generated by CLUSTAL_X version 1.8 (Thompson et al., 1997Down). Analysis of nucleotide and deduced amino acid sequence identities was performed using GeneDoc and DNASTAR (Lasergene). Neighbour-joining phylogenetic trees based on nucleotide sequences were constructed using MEGA version 3.1 (Kumar et al., 2004Down). The robustness of phylogenetic constructions was evaluated by bootstrapping using 500 replications. Previously published GETV and other alphavirus sequences used in this study include sequences from O'Nyong-nyong virus (ONNV; GenBank accession number M20303 [GenBank] ), Barmah Forest virus (BFV; U73745 [GenBank] ), MAYV (AF237947 [GenBank] ), SFV (X04129 [GenBank] ), RRV (K00046), CHIKV (NC_004162 [GenBank] ), GETV (AF268025 [GenBank] , AF339484 [GenBank] , AY702913 [GenBank] and EF011023 [GenBank] ) and SAGV (AB032553 [GenBank] and AF242890 [GenBank] ).


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Virus isolation and identification
Isolates that produced CPE in cell culture were targeted for serological and molecular identification. Antigen ELISA using standard anti-alphavirus antibodies identified isolates YN0540, YN0542, SH05-6, SH05-15, SH05-16, SH05-17 and GS10-2 as alphaviruses. These virus isolates caused CPE in BHK-21 cells, which was characterized by aggregation and fusion by day 3 post-infection. RT-PCR amplification using alphavirus primers (particularly GETV-derived primers) followed by sequence analysis of PCR products (and the amino acid sequences that they encode) showed that they exhibited high amino acid sequence identity with GETV. Isolates from Yunnan province were designated YN0540 and YN0542, the single isolate from Gansu province was designated GS10-2 and the Shanghai isolates were designated SH05-6, SH05-15, SH05-16 and SH05-17. Previously reported isolates HB0215-3 and HB0234 (Wang et al., 2006Down) were also included in this study.

Sequence analysis
Full-length genome sequences were determined for three isolates, M1, HB0234 and YN0540. Details are provided in Table 1Down. It is noteworthy that the non-structural ORF of strain HB0234 was 3 nucleotides shorter than those of the two other strains, as a result of a deletion of a codon at positions 5473–5475 when compared with SAGV (GenBank accession number AB032533 [GenBank] ). A deletion of 3 nucleotides (positions 99–102) was identified in the capsid gene of strain M1 when compared with SAGV (AB032533 [GenBank] ). Sequence analysis of 3' UTRs of isolates HB0215-3, YN0542, SH05-6 and SH05-17 indicated that their UTRs were 401 nt long, while those of isolates SH05-15, SH05-16 and GS10-2 were respectively 404, 410 and 402 nt long.


View this table:
[in this window]
[in a new window]

 
Table 1. Genome lengths and other details for the fully sequenced strains M1, HB0234 and YN0540

 
Comparison of the full-length sequences of GETV isolates M1, HB0234 and YN0540
The full-length genome sequences of GETV strains M1, HB0234 and YN0540, as well as three GETV strains from South Korea (GenBank accession number AY702913 [GenBank] ), Russia (EF631998 [GenBank] ) and Mongolia (EF631999 [GenBank] ) and one isolate of SAGV (AB032553 [GenBank] ) obtained from databases, were compared.

The nucleotide sequence identity between the three Chinese isolates (M1, HB0234 and YN0540) ranged from 98.1 to 99.1 %, while amino acid sequence identity ranged from 98.7 to 99.8 % (Tables 2Down and 3Down). The nucleotide sequence identities between the Chinese and other isolates ranged from 96.9 to 98.8 %, while the amino acid sequence identity ranged from 98.7 to 99.8 % (98.9–99.6 % in the non-structural polyprotein and 98.0–99.8 % in the structural polyprotein).


View this table:
[in this window]
[in a new window]

 
Table 2. Sequence divergence of M2021, M1, HB0234, YN0540 and South Korea GETV strains from SAGV

Numbers and percentages (in parentheses) of divergent nucleotides (including insertions and deletions) and amino acids are shown. There is only an incomplete sequence for the original Malaysian GETV strain (M2021). Junction refers to the non-coding region between the non-structural and structural ORFs.

 

View this table:
[in this window]
[in a new window]

 
Table 3. Amino acid deletions or changes found in virus isolates analysed in this study

Positions are given relative to the prototype Malaysian strain M2021 of GETV. Deletions are indicated by Del.

 
Analysis of all sequences included in this study (GETV AY702913 [GenBank] , YN0540, HB0234, M1, LEIV 16275 Mag, LEIV 17741 MPR and SAGV) showed that the lowest nucleotide sequence identity was that between the South Korean GETV isolate and SAGV (96.9 %). The lowest amino acid sequence identity was found between the Chinese isolate GETV-HB0234 or GETV-M1 and SAGV (98.7 %).

Amino acid changes were identified in various genes, as shown in Table 3Up. Deletion of a single amino acid was detected in the nsP3 of isolate HB0234. The nsP3–nsP4 junction of the same isolate is not interrupted by the usual leaky opal stop codon (TGA) that is found at position 518 of GETV; rather, this is replaced by a TGG (coding for Trp).

In the capsid protein, seven positions were found to be variable among all analysed GETV isolates. Deletion of a single amino acid was found in the capsid protein of strain M1. Again, the M1 strain showed a single amino acid change in the E3 protein. Seventeen sites showed changes in the E2 protein in the various Chinese isolates. In the 6K protein, only SAGV presented a single change. The results are detailed in Tables 3Up and 4Down.


View this table:
[in this window]
[in a new window]

 
Table 4. Identity analysis of the sequences of GETV and SAGV

Percentage identity was determined from pairwise comparisons of full-length nucleotide sequences (above the diagonal) and amino acid sequences (below the diagonal).

 
Analysis of the E2 sequence
The length of the E2 gene of each of the seven strains (HB0215-3, YN0542, SH05-6, SH05-15, SH05-16, SH05-17 and GS10-2) was found to be 1266 nt, encoding a 422 aa glycoprotein. The E2 sequences of all ten Chinese isolates of GETV were analysed and detailed identity values are reported in Table 5Down. These values ranged between 97.7 and 99.9 % for nucleic acids and 98.3 and 100 % for amino acids among the Chinese isolates and between 94.5 and 95.3 % for nucleic acids and 96.9 to 97.6 % for amino acids when the Chinese isolates were compared with the original Malaysian GETV isolate M2021 (GenBank accession number AF339484 [GenBank] ). Interestingly, the Mongolian isolate (LEIV 17741 MPR; GenBank accession number EF631999 [GenBank] ) showed 100 % amino acid identity with four Chinese isolates (SH05-6, SH05-15, YN0540 and YN0542) isolated in the south of China.


View this table:
[in this window]
[in a new window]

 
Table 5. Sequence identities for the E2 gene of GETV and SAGV

Percentage identity was determined from pairwise comparisons of nucleotide sequences (above the diagonal) and amino acid sequences (below the diagonal) of the E2 gene.

 
As shown in Supplementary Fig. S1, 17 positions were found to be variable in the sequence of the E2 protein of GETV isolates (compared with the original Malaysian isolate). This reflects a relatively good conservation in the E2 amino acid sequence, even though they have been isolated from different geographical locations and distinct vectors.

Analysis of the 3' UTR of the Chinese isolates of GETV
Alphaviruses possess a highly conserved 3' sequence element (3' CSE; approximately 19 nt long) which immediately precedes the poly(A) tail (Pfeffer et al., 1998Down). Both the poly(A) tail and the 3' CSE are required for virus replication and, more specifically, for efficient minus-strand RNA synthesis (Kuhn et al., 1990Down, 1992Down; Raju et al., 1999Down; Hardy & Rice, 2005Down). The terminal 19 nt conserved sequence was identical in all GETV isolates, including the M1 isolate that was previously reported to have lost this conserved sequence (Zhao et al., 2000Down).

Within the 3' UTR, repeated sequence elements (RSE) of 40–60 nt have been described for all alphaviruses investigated so far (Strauss & Strauss, 1994Down). These were identified as being three repeats of about 54 nt in GETV (Supplementary Fig. S2; underlined) (Pfeffer et al., 1998Down). An analysis of the 3' UTR of GETV isolates identified considerable variations, as shown in Supplementary Fig. S2. Insertions and/or deletions were found at several positions in the 3' UTR. A deletion was identified in the 3' UTR of a previously published GETV sequence (GenBank accession number AY702913 [GenBank] ) between nucleotides 50 and 153 (positions are relative to the SAGV sequence).

Sequence alignment showed that all isolates including the previously published Taiwan district (GenBank accession number AF242890 [GenBank] ), Russian (EF631998 [GenBank] ) and Mongolian (EF631999 [GenBank] ) isolates show the same deletion (10 nt) at positions 45–54 in the 3' UTR (except YN0540, which showed a deletion of only 9 nt, at positions 46–54).

The Chinese isolates, including the Taiwan district isolate, all showed a T to C transition at position 64 and a C to T transition at positions 43 and 148 (relative to the SAGV sequence, GenBank accession number AB032533 [GenBank] ). These strains were isolated up to 40 years apart (the M1 strain was isolated in 1964). Insertions (or deletions) of adenine residues were identified at several positions (Supplementary Fig. S2; positions 173, 299, 250 and 372). These included an insertion of 10 consecutive adenine residues in the 3' UTR of strain M1 (positions 189–198), which is located towards the end of the second RSE. Insertions of 9 and 3 consecutive adenine residues, respectively, were also identified towards the end of the third RSE in isolates SH05-6 and SH05-15 (Supplementary Fig. S2).

Phylogenetic analysis
The P-distance model and Kimura's two-parameter model were used for tree building. The alpha shape parameter used for the gamma distribution analysis (for the Kimura two-parameter) was calculated using the PAML package (Yang, 1997Down). The shape parameter measures how variable the rates are among sites. When the value of alpha is >1, most sites have rates around 1 (similar rates). When alpha is ≤1 (there is a relatively large amount of rate variation), most sites have very low rates (<1), but there are evolutionary hot spots with higher rates.

A phylogenetic tree (Fig. 2aDown) was constructed using the complete genome sequences of several alphaviruses including the GETV isolates M1, HB0234 and YN0540 from China and other GETV isolates from South Korea, Russia and Mongolia. Because the sequence of the original Malaysian strain of GETV is incomplete, only the sequence encoding nsp4 and the structural polyprotein (4789 nt) of this particular strain was used in the alignment. The sequence encoding nsP1, nsP2 and nsP3 of this isolate was accounted for as missing information in the MEGA 3.1 package options.


Figure 2
View larger version (38K):
[in this window]
[in a new window]

 
Fig. 2. Phylogenetic relationship between Chinese isolates of GETV and other alphaviruses (particularly those of the SFV and Barmah forest serocomplexes) based on complete nucleotide or capsid protein or E2 nucleotide sequences. (a) Phylogenetic tree constructed using full-length genome sequences of Chinese isolates M1, HB0234 and YN0540. (b) Phylogenetic tree constructed using the nucleotide sequences of the capsid gene of Chinese isolates M1, HB0234 and YN0540. (c) Phylogenetic tree constructed using the nucleotide sequences of the E2 protein of GETV isolates. Bootstrap values for 500 replications are indicated at each node. Bars, 0.05 (a, b) or 0.005 (c) substitutions per site.

 
Another phylogenetic tree (Fig. 2bUp) was constructed using nucleotide sequences encoding the capsid protein and included, alongside several other alphaviruses, the GETV isolates M1, HB0234 and YN0540 from China, GETV from Malaysia, SAGV from Japan and other GETV isolates from South Korea, Russia and Mongolia.

A third phylogenetic tree (Fig. 2cUp) was constructed from the E2 protein sequences of all Chinese isolates of GETV alongside the original GETV from Malaysia, SAGV from Japan and other GETV isolates from South Korea, Russia and Mongolia.

Trees constructed for given genes using either the P-distance or Kimura's two-parameter model were found to have identical topologies.

The trees showed that the GETV and SAGV isolates clustered onto a distinct evolutionary branch, rooted by the original GETV from Malaysia. The various trees showed that the most recent isolates of GETV are located at the end branches, rooted by the earlier isolates (the Chinese isolates represent the latest evolutionary events pertaining to GETV). The only noticeable exception is the case of the recent Mongolian and Russian isolates (made in 2000). These two strains do not cluster among those that have been made more recently around the same time. Rather, the Russian isolate made in 2000 roots the M1 Chinese strain (isolated in 1964) which, in turn, roots the Mongolian isolate, also made in 2000. These data suggest that GETV has been evolving slowly (possibly through slow circulation of the virus) and recent years account for the faster evolution of the more recent Chinese isolates. We have calculated the various shape parameters for the available sets of GETV E1, E2 and capsid sequences. Shape parameter values were between 0.21 and 0.23. These low values of the shape parameter are an indication that most of the sites are evolving slowly, but that there are evolutionary hot-spots in these sequences that are driving the evolution of these sequences of the genome.


   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
GETV was originally isolated in Malaysia in 1955. The following year, a related virus was isolated from mosquitoes in Japan and designated Sagiyama virus (Scherer et al., 1962Down). Differences in antigenicity between GETV and SAGV were only identified by a complement fixation (CF) test, while virus neutralization and haemagglutination assays could not distinguish GETV from SAGV (Kono et al., 1980Down). Previous sequence analyses showed that GETV and SAGV should be classified as members of the same species within the genus Alphavirus (Shirako & Yamaguchi, 2000Down; Attoui et al., 2007Down).

In this study, we analysed the sequences of ten GETV strains that were isolated in China during investigations for arboviruses between 1964 and 2006. Seven strains were isolated during the current study in regions within the range of latitude 1 ° (Hainan) and 4 °(Hebei and Ganshu) N and longitude 9 ° (Ganshu and Yunnan) and 12 ° (Shanghai) E. These isolations showed that GETV is distributed widely in China.

Non-plaque-purified virus material was used for the sequencing strategies, which is based on that published by Shirako & Yamaguchi (2000)Down for the sequence determination of SAGV. The PCR products were not cloned but rather sequenced directly, in order to avoid sequence variations that might result from the presence of quasispecies.

Our analysis of the full-length sequences of the Chinese GETV isolates, South Korean GETV isolate and SAGV (Japanese isolate) showed that these viruses exhibit high identity, therefore confirming previous findings about the relationship of GETV and SAGV. Up to 5.5 % amino acid sequence divergence was calculated between the original Malaysian isolate of GETV and the Chinese isolates.

The E2 protein is responsible for eliciting virus-specific antibodies (Griffin, 2001Down). We found 17 variable sites among the 14 GETV strains from China, Malaysia, Russia, Mongolia and South Korea (Supplementary Fig. S2). Our results showed that some amino acid changes in E2 are specific to isolates from a particular geographical location. By contrast, the number of variable sites between the E1 proteins of the various Chinese GETV isolates is quite low (only four sites among M1, HB0234 and YN0540), which reflects sequence conservation in the E1 of GETV.

In some arthropod-borne viruses, the nucleocapsid is the antigen responsible for reaction with CF antibodies (Casals & Whitman, 1961Down; Bose & Sagik, 1970Down). A CF assay distinguished between SAGV and GETV, while haemagglutination and seroneutralization did not (Kono et al., 1980Down). Wekesa et al. (2001)Down analysed the capsid gene of GETV, SAGV isolates in Japan and the Malaysian GETV (M2021) and found that the difference between GETV and SAGV was located at amino acid position 100 (Pro to Leu change). We found that this particular change was also present between the Chinese isolates and SAGV. All of the Chinese isolates have the same amino acid difference (Pro to Leu) from SAGV.

These results show the particularity of the Malaysian isolate and the close relationship between the Chinese GETV isolates and SAGV. In agreement with these results and previous findings (Powers et al., 2001Down; Attoui et al., 2007Down); the latter virus should therefore be identified as a member of the species Getah virus and not as a member of the species Ross River virus (according to the current classification; Weaver et al., 2005Down).

Read-through of leaky opal (TGA) stop codons found upstream of the nsP3/nsP4 cleavage site in most alphaviruses generates the non-structural polyprotein nsP1234 (Kim et al., 2004Down; Lavergne et al., 2006Down). The read-through efficiency in alphaviruses was estimated to be around 5–20 % (Strauss & Strauss, 1994Down), and a TGA to TGG mutation is supposed to increase the expression of nsP1234, which might relate to the efficiency of virus replication. However, in both SFV and SINV (Shirako & Strauss, 1994Down; Kim et al., 2004Down), this has not been the case, because free nsP4 was found to be degraded rapidly by the N-end-rule pathway (Takkinen et al., 1991Down; Ferreira et al., 2003Down).

Like a few other alphaviruses, strain HB0234 did not possess a stop codon between nsp3 and nsp4, but rather this was replaced by a sense codon. This suggests that strain HB0234 would express the non-structural polyprotein nsP1234 without being dependent on read-through. Our results are in agreement with the previous finding concerning SFV and SINV, as we did not observe any significant differences in proliferation between HB0234 and those Chinese GETV strains that possess the leaky stop codon upstream of the nsP3/nsP4 cleavage site (data not shown).

The 19 nt CSE of the 3' UTR is identical in the various GETV isolates, which suggests an important role of the 3' UTR in viral replication in the host (Hajjou et al., 1996Down). Mutations were previously reported in the 3' UTR among different alphaviruses, including high rates of mutations in RRV and Venezuelan equine encephalitis virus (Pfeffer et al., 1998Down). We found a poly(A) insertion in the 3' UTR which is located at the end of the second and/or the third RSE in GETV. The number of poly(A) residues in the insertion that we identified in strain M1 was different from that reported by Zhao et al. (2000)Down (four adenines; GenBank accession number EF011023 [GenBank] ), which suggest that this might be related to the passage history of the virus.

Several arbovirus investigations were conducted between the 1970s and 1990s in provinces of China (Zhang et al., 1989Down; Chen & Tao, 1996Down; Liang et al., 2000Down; Huang et al., 2001Down; Lv et al., 2001Down; Wang & Liang., 2003Down). None of these investigations reported any isolation of GETV. It has only been since 2000 that GETV strains have been isolated from wild-caught mosquitoes.

The area of GETV dissemination in Russia and Mongolia includes steppe, mixed forest, northern taiga and forest-tundra zones, reaching the tundra zone in the north (L'vov et al., 2000). GETV is the only alphavirus that occurs under such severe climatic conditions in Russia (L'vov et al., 2000). A plausible explanation for the position of the Russian and Mongolian isolates in the various trees would be that these particular strains have not been circulating extensively (and hence have not accumulated as many mutations as the more recent Chinese isolates), given the extreme climatic conditions in the regions where they have been isolated. Various mosquitoes and birds have been found to be infected by GETV (Takashima & Hashimoto, 1985Down; Kono, 1988Down) and the virus is supposed not to induce persistent infection in mammals (Wekesa et al., 2001Down; Chang et al., 2006Down). Persistent infection in naturally infected mammals as a possible mechanism of overwintering would be ideal to explain the clustering of the Russian and Mongolian isolates of GETV (made in 2000) nearer to the roots than to the exterior branches of the various phylogenetic trees (as is the case of the recent Chinese isolates, made around similar dates).

Finally, the expansion of GETV in China is a recent phenomenon. Beside the original M1 strain that was isolated in 1964 (almost 40 years before the most recent isolations), previous investigation in Gansu, Xinjiang, Yunnan and Hebei provinces over 10 years ago failed to identify any GETV in China (Zhang et al., 1989Down; Chen & Tao, 1996Down; Liang et al., 2000Down; Huang et al., 2001Down; Lv et al., 2001Down; Wang & Liang, 2003Down). This should prompt further investigations as a follow up for the dissemination of this virus in China. The recent isolations of 15 strains of GETV in Russia and Mongolia (L'vov et al., 2000) and the widespread distribution of this virus in China are indications of a change in the world distribution of the re-emerging GETV. This change and the isolation of GETV from various mosquito species of the genera Culex, Aedes and Armigeres identify this virus as a pathogen that is an increasing threat to animal health, particularly horses and pigs.


   ACKNOWLEDGEMENTS
 
We are grateful to Dr Wu-yang Zhu and Dr Xiao-yan Gao in our department for helpful comments on the manuscript. This work was supported by a grant (to G.-d. L.) from the Ministry of Science and Technology of China (no. 2003BA712A08-01) and by the Japan Sciences Health Foundation.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Attoui, H., Sailleau, C., Mohd Jaafar, F., Belhouchet, M., Biagini, P., Cantaloube, J. F., de Micco, P., Mertens, P. & Zientara, S. (2007). Complete nucleotide sequence of Middelburg virus, isolated from the spleen of a horse with severe clinical disease in Zimbabwe. J Gen Virol 88, 3078–3088.[Abstract/Free Full Text]

Berge, T. O. (1975). Getah. In International Catalog of Arboviruses, 2nd edn, pp. 278–279. Washington, DC: US Department of Health, Education, and Welfare.

Bose, H. R. & Sagik, B. P. (1970). Immunological activity associated with the nucleocapsid and envelope components of an arbovirus. J Virol 5, 410–412.[Abstract/Free Full Text]

Brown, C. M. & Timoney, P. J. (1998). Getah virus infection of Indian horses. Trop Anim Health Prod 30, 241–252.[CrossRef][Medline]

Bryant, J. E., Crabtree, M. B., Nam, V. S., Yen, N. T., Duc, H. M. & Miller, B. R. (2005). Isolation of arboviruses from mosquitoes collected in northern Vietnam. Am J Trop Med Hyg 73, 470–473.[Abstract/Free Full Text]

Calisher, C. H., Shope, R. E., Brandt, W., Casals, J., Karabatsos, N., Murphy, F. A., Tesh, R. B. & Wiebe, M. E. (1980). Proposed antigenic classification of registered arboviruses I. Togaviridae, Alphavirus. Intervirology 14, 229–232.[Medline]

Casals, J. & Whitman, L. (1961). Group C, a new serological group of hitherto undescribed arthropod-borne viruses. Immunological studies. Am J Trop Med Hyg 10, 250–258.[Abstract/Free Full Text]

Chang, C. Y., Huang, C. C., Huang, T. S., Deng, M. C., Jong, M. H. & Wang, F. I. (2006). Isolation and characterization of a Sagiyama virus from domestic pigs. J Vet Diagn Invest 18, 156–161.[Abstract/Free Full Text]

Chen, B. Q. & Tao, S. J. (1996). Arbovirus survey in China in recent ten years. Chin Med J (Engl) 109, 13–15.[Medline]

Faragher, S. G., Meek, A. D. J., Rice, C. M. & Dalgarno, L. (1988). Genome sequences of a mouse-avirulent and a mouse-virulent strain of Ross River virus. Virology 163, 509–526.[CrossRef][Medline]

Ferreira, D., Hernandez, R., Horton, M. & Brown, D. T. (2003). Morphological variants of Sindbis virus produced by a mutation in the capsid protein. Virology 307, 54–66.[CrossRef][Medline]

Griffin, D. E. (2001). Alphaviruses. In Fields Virology, 4th edn, pp. 917–962. Edited by D. M. Knipe & P. M. Howley. Philadelphia: Lippincott Williams and Wilkins.

Hajjou, M., Hill, K. R., Subramaniam, S. V., Hu, J. Y. & Raju, R. (1996). Nonhomologous RNA-RNA recombination events at the 3' nontranslated region of the Sindbis virus genome: hot spots and utilization of nonviral sequences. J Virol 70, 5153–5164.[Abstract/Free Full Text]

Hardy, R. W. & Rice, C. M. (2005). Requirements at the 3' end of the sindbis virus genome for efficient synthesis of minus-strand RNA. J Virol 79, 4630–4639.[Abstract/Free Full Text]

Huang, W.-L., Zhang, H.-L., Wang, G.-L., Zhang, Y.-Z., Zhou, G.-L., Mi, Z.-Q., Shi, H.-F., Zi, D.-Y. & Liang, G.-D. (2001). Study on biological characters of Sindbis virus Yunnan strain. Virol Sin 16, 97–100 (in Chinese with English abstract).

Kamada, M., Ando, Y., Fukunaga, Y., Kumanomido, T., Imagawa, H., Wada, R. & Akiyama, Y. (1980). Equine Getah virus infection: isolation of the virus from racehorses during an enzootic in Japan. Am J Trop Med Hyg 29, 984–988.[Abstract/Free Full Text]

Kanamitsu, M., Taniguchi, K., Urasawa, S., Ogata, T., Wada, Y., Wada, Y. & Saroso, J. S. (1979). Geographic distribution of arbovirus antibodies in indigenous human populations in the Indo-Australian archipelago. Am J Trop Med Hyg 28, 351–363.[Abstract/Free Full Text]

Kim, K. H., Rumenapf, T., Strauss, E. G. & Strauss, J. H. (2004). Regulation of Semliki Forest virus RNA replication: a model for the control of alphavirus pathogenesis in invertebrate hosts. Virology 323, 153–163.[CrossRef][Medline]

Kono, Y. (1988). Getah virus disease. In: The Arboviruses. Epidemiology and Ecology, pp. 21–36. Edited by T. P. Monath. Boca Raton, FL: CRC Press.

Kono, Y., Sentsui, H. & Ito, Y. (1980). An epidemic of Getah virus infection among racehorses: properties of the virus. Res Vet Sci 29, 162–167.[Medline]

Kuhn, R. J., Hong, Z. & Strauss, J. H. (1990). Mutagenesis of the 3' nontranslated region of Sindbis virus RNA. J Virol 64, 1465–1476.[Abstract/Free Full Text]

Kuhn, R. J., Griffin, D. E., Zhang, H., Niesters, H. G. & Strauss, J. H. (1992). Attenuation of Sindbis virus neurovirulence by using defined mutations in nontranslated regions of the genome RNA. J Virol 66, 7121–7127.[Abstract/Free Full Text]

Kumanomido, T., Wada, R., Kanemaru, T., Kamada, M., Hirasawa, K. & Akiyama, Y. (1988a). Clinical and virological observations on swine experimentally infected with Getah virus. Vet Microbiol 16, 295–301.[CrossRef][Medline]

Kumanomido, T., Kamada, M., Wada, R., Kenemaru, T., Sugiura, T. & Akiyama, Y. (1988b). Pathogenicity for horses of original Sagiyama virus, a member of the Getah virus group. Vet Microbiol 17, 367–373.[CrossRef][Medline]

Kumar, S., Tamura, K. & Nei, M. (2004). MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 5, 150–163.[Abstract/Free Full Text]

Lavergne, A., de Thoisy, B., Lacoste, V., Pascalis, H., Pouliquen, J. F., Mercier, V., Tolou, H., Dussart, P., Morvan, J. & other authors (2006). Mayaro virus: complete nucleotide sequence and phylogenetic relationships with other alphaviruses. Virus Res 117, 283–290.[CrossRef][Medline]

Li, X. D., Qiu, F. X., Yang, H., Rao, Y. N. & Calisher, C. H. (1992). Isolation of Getah virus from mosquitoes collected on Hainan Island, China, and results of a serosurvey. Southeast Asian J Trop Med Public Health 23, 730–734.[Medline]

Liang, G.-D., Li, L., Zhou, G.-L., Fu, S.-H., Li, Q.-P., Li, F.-S., He, H.-H., Jin, Q., He, Y. & other authors (2000). Isolation and complete nucleotide sequence of a Chinese Sindbis-like virus. J Gen Virol 81, 1347–1351.[Abstract/Free Full Text]

Lv, X.-J., Fu, S.-H., Yang, Y.-L., He, H.-H., Zhang, G.-Y., Chen, X.-W., Liang, G.-D., Jin, Q. & Hou, Y.-D. (2001). XJ-90260 virus isolated in China was identified as Western equine encephalitis virus. Chin J Virol 17, 112–115 (in Chinese).

L'vov, S. D., Gromashevski'i, V. L., Aristova, V. A, Morozova, T. N., Skvortsova, T. M., Gushchina, E. A., Petrova, E. S. & L'vov, D. K. (2000). Isolation of Getah virus (Togaviridae, Alfavirus) strains in North-Eastern Asia. Vopr Virusol 45, 14–18 (in Russian).[Medline]

Marchette, N. J., Rudnick, A. & Garcia, R. (1980). Alphaviruses in Peninsular Malaysia. II. Serological evidence of human infection. Southeast Asian J Trop Med Public Health 11, 14–23.[Medline]

Norder, H., Lundstrom, J. O., Kozuch, O. & Magnius, L. O. (1996). Genetic relatedness of Sindbis virus strains from Europe, Middle East, and Africa. Virology 222, 440–445.[CrossRef][Medline]

Pfeffer, M., Proebster, B., Kinney, R. M. & Kaaden, O. R. (1997). Genus-specific detection of alphaviruses by a semi-nested reverse transcription-polymerase chain reaction. Am J Trop Med Hyg 57, 709–718.[Abstract/Free Full Text]

Pfeffer, M., Kinney, R. M. & Kaaden, O. R. (1998). The alphavirus 3'-nontranslated region: size heterogeneity and arrangement of repeated sequence elements. Virology 240, 100–108.[CrossRef][Medline]

Powers, A. M., Brault, A. C., Shirako, Y., Strauss, E. G., Kang, W., Strauss, J. H. & Weaver, S. C. (2001). Evolutionary relationships and systematics of the alphaviruses. J Virol 75, 10118–10131.[Abstract/Free Full Text]

Raju, R., Hajjou, M., Hill, K. R., Botta, V. & Botta, S. (1999). In vivo addition of poly(A) tail and AU-rich sequences to the 3' terminus of the Sindbis virus RNA genome: a novel 3'-end repair pathway. J Virol 73, 2410–2419.[Abstract/Free Full Text]

Scherer, W. F., Funkenbusch, M., Buescher, E. L. & Izumit, T. (1962). Sagiyama virus, a new group A arthropod-borne virus from Japan. I. Isolation, immunologic classification, and ecologic observations. Am J Trop Med Hyg 11, 255–268.[Abstract/Free Full Text]

Shibata, I., Hatano, Y., Nishimura, M., Suzuki, G. & Inaba, Y. (1991). Isolation of Getah virus from dead fetuses extracted from a naturally infected sow in Japan. Vet Microbiol 27, 385–391.[CrossRef][Medline]

Shirako, Y. & Strauss, J. H. (1994). Regulation of Sindbis virus RNA replication: uncleaved P123 and nsP4 function in minus strand RNA synthesis whereas cleaved products from P123 are required for efficient plus strand RNA synthesis. J Virol 68, 1874–1885.[Abstract/Free Full Text]

Shirako, Y. & Yamaguchi, Y. (2000). Genome structure of Sagiyama virus and its relatedness to other alphaviruses. J Gen Virol 81, 1353–1360.[Abstract/Free Full Text]

Strauss, J. H. & Strauss, E. G. (1994). The alphaviruses: gene expression, replication, and evolution. Microbiol Rev 58, 491–562.[Abstract/Free Full Text]

Takashima, I. & Hashimoto, N. (1985). Getah virus in several species of mosquitoes. Trans R Soc Trop Med Hyg 79, 546–550.[CrossRef][Medline]

Takkinen, K., Peranen, J. & Kaariainen, L. (1991). Proteolytic processing of Semliki Forest virus-specific nonstructural polyprotein. J Gen Virol 72, 1627–1633.[Abstract/Free Full Text]

Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F. & Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25, 4876–4882.[Abstract/Free Full Text]

Turell, M. J., O'Guinn, M. L., Wasieloski, L. P., Jr, Dohm, D. J., Lee, W. J., Cho, H. W., Kim, H. C., Burkett, D. A., Mores, C. N. & other authors (2003). Isolation of Japanese encephalitis and Getah viruses from mosquitoes (Diptera: Culicidae) collected near Camp Greaves, Gyonggi Province, Republic of Korea, 2000. J Med Entomol 40, 580–584.[Medline]

Wang, H. Y. & Liang, G. D. (2003). Arbovirus in China in recent ten years. Chin J Public Health Rev 19, 473–476 (in Chinese).

Wang, H. Q., Liu, W. B., Yang, D. R., Liang, Y., Wang, J. W., Zhang, L. S., Liu, J. W., Tao, S. J., Lv, X. J. & other authors (2006). Isolation and identification of arboviruses in Hebei Province. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 20, 52–55 (in Chinese).[Medline]

Weaver, S. C., Frey, T. K., Huang, H. V., Kinney, R. M., Rice, C. M., Roehrig, J. T., Shope, R. E. & Strauss, E. G. (2005). Togaviridae. In Virus Taxonomy. Eighth Report of the International Committee on Taxonomy of Viruses, pp. 999–1008. Edited by C. M. Fauquet, M. A. Mayo, J. Maniloff, U. Desselberger & L. A. Ball. London: Elsevier/Academic Press.

Wekesa, S. N., Inoshima, Y., Murakami, K. & Sentsui, H. (2001). Genomic analysis of some Japanese isolates of Getah virus. Vet Microbiol 83, 137–146.[CrossRef][Medline]

Wen, J. S., Zhao, W. Z., Liu, J. W., Zhou, H., Tao, J. P., Yan, H. J., Liang, Y., Zhou, J. J. & Jiang, L. F. (2007). Genomic analysis of a Chinese isolate of Getah-like virus and its phylogenetic relationship with other alphaviruses. Virus Genes 35, 597–603.[CrossRef][Medline]

Yang, Z. (1997). PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 13, 555–556.[Free Full Text]

Zhang, H. L., Shi, H. F., Liu, L. H., Yu, Y. X., Zi, D. Y., Li, Z. X., Zhang, T. S., Cui, W. Q., Wang, Z. W. & other authors (1989). Isolation of chikungunya virus from bat in Yunnan Province and serological investigations. Bing Du Xue Bao 5, 31–36 (in Chinese).

Zhao, C. S., Jiang, L. H., Yu, X. L., Chen, W. J., Rao, Y. N. & Lu, Z. Z. (1997). Isolation of Ross River virus and its antibody prevalence in Hainan province. Chin J Vet Sci 17, 241–243 (in Chinese).

Zhao, W., Zhou, G. & He, H. (2000). Cloning and primary analysis of 3' end genome of two alphaviruses isolated from Hainan Province of China. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 14, 213–217 (in Chinese).[Medline]

Received 21 November 2007; accepted 21 February 2008.



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table and Figures
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via CrossRef
Google Scholar
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.
Agricola
Right arrow Articles by Zhai, Y.-g.
Right arrow Articles by Liang, G.-d.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS