Previous Article | Next Article ![]()
Journal of Clinical Microbiology, October 2002, p. 3684-3688, Vol. 40, No. 10
0095-1137/02/$04.00+0 DOI: 10.1128/JCM.40.10.3684-3688.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Laboratory of Clinical Virology, Kyushu Research Station, National Institute of Animal Health, Kagoshima 891-0105, Japan
Received 26 December 2001/ Returned for modification 21 April 2002/ Accepted 22 July 2002
|
|
|---|
|
|
|---|
In the late summer to autumn of 1997, a large-scale epidemic of Ibaraki disease occurred in western Japan. During the same period of time, numerous miscarriages and stillbirths occurred in apparently healthy cattle in the same area. In the 1997 epidemic, a new strain of IBAV was isolated not only from cattle showing Ibaraki disease but also from aborted fetuses, their dams, and biting midges, as well as from the blood of both affected and apparently healthy cattle. The new virus was distinguished from the existing strains of IBAV by a cross-neutralization test and by PCR-restriction fragment length polymorphism analysis (21).
The IBAV genome is composed of 10 double-stranded (ds) RNA segments in a bilayered capsid. Segments 3 and 7 encode the inner capsid proteins VP3 and VP7, respectively. These capsid proteins play the role of serogroup-specific antigens. Segments 2 and 6 encode the outer layer proteins VP2 and VP5, respectively, of which VP2 is a major neutralizing antigen and serotype-specific determinant (14). The 10 genomic segments have the potential to be exchanged within the serogroup (reassortment). In another orbivirus, bluetongue virus (BTV), reassortment events in tissue culture, in the vector, and in nature have been described (7, 9, 26). The molecular characterization of IBAV would be valuable for investigating the viral evolution and epidemiology of the disease.
In this report, four RNA segments of IBAV, including the new strain, were sequenced and compared with the published sequence data in order to demonstrate the genetic heterogeneity within the EHDV serogroup, along with the consequences of virus evolution.
|
|
|---|
Extraction of viral dsRNAs. Viral dsRNAs were extracted from the infected BHK-21 cells as described previously (28), and the individual RNA segments were separated on a 0.8% agarose gel. The separated viral dsRNAs were excised from the agarose gel and purified with the RNaid kit (Bio 101, La Jolla, Calif.) in accordance with the manufacturer's instructions.
Genomic amplification and cloning. The primers corresponding to the 5' and 3' termini of segments 2, 3, 6, and 7 were synthesized on the basis of the sequences of the No.2 strain of IBAV (Table 1) (15, 19, 20). Reverse transcription-PCR (RT-PCR) was performed with the Titan One Tube RT-PCR kit (Roche Diagnostics, Indianapolis, Ind.) containing a proofreading polymerase. The RT-PCR mixture consisted of 0.2 mM (each) deoxynucleoside triphosphate, 5 mM dithiothreitol, 1x RT-PCR buffer, 1.5 mM MgCl2, and an enzyme mixture in a total reaction volume of 50 µl. The purified viral RNA segment and both strand primers were incubated at 94°C for 5 min and quenched on ice prior to the addition of the reaction mixture. The RT was conducted at 48°C for 60 min. This was followed by an initial 3-min denaturation at 95°C and 10 subsequent cycles of denaturation at 95°C for 30 s, annealing at 55 to 63°C (depending on the primer pair) for 30 s, and elongation at 68°C for 2 min, and then by 25 cycles of denaturation at 95°C for 30 s, annealing at 55 to 63°C for 30 s, and elongation at 68°C for 2 min at first, though with each cycle the elongation time was incremented by 5 s.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers used for the cDNA amplification and sequencing
|
Sequencing. The sequencing primers are shown in Table 1. Three individual clones containing the viral genome were sequenced in both directions by using a model 4200 automated DNA sequencer (Li-Cor Inc., Lincoln, Nebr.) with the Thermo Sequenase cycle sequencing kit (U.S. Biochemicals, Cleveland, Ohio).
Genetic analysis. The nucleotide sequences obtained in this study were edited with the GENETYX software. The multiple alignments of the deduced amino acid sequences encoded by the open reading frames (ORF) were done with the CLUSTAL W program (31). The aligned sequences were used to construct the phylogenetic trees by the neighbor-joining method (27), and the distances were corrected by using Kimura's two-parameter method (16). The reliability of the branching orders was estimated by bootstrapping (1,000 samples). The phylogenetic trees were drawn with TreeView software (24).
Sequence data. The sequence data determined in this study and previously reported nucleotide sequences of IBAV and Australian and North American EHDV segments 2, 3, 6, and 7 were used in the genetic analysis. The sequence data were taken from the nucleotide database (the accession numbers are shown in the legends for Fig. 1 and 2).
![]() View larger version (12K): [in a new window] |
FIG. 1. Phylogenetic tree based on VP3 and VP7 of IBAV and EHDV. Trees were constructed by the neighbor-joining method and bootstrap test (n = 1,000). EHDV-1 was used as an outgroup to root the tree. Accession numbers are as follows: EHDV-1 segment 3, X61589; EHDV-1 segment 7, D10766; SV 124 (EHDV-2) segment 3, L33819; SV 124 segment 7, AF188643; CSIRO439 (EHDV-2) segment 3, S68010; No.2 (IBAV) segment 3, AB04930; No.2 segment 7, AB041934.
|
![]() View larger version (11K): [in a new window] |
FIG. 2. Phylogenetic tree based on VP2 and VP5 of IBAV and EHDV. Trees were constructed by the neighbor-joining method and bootstrap test (n = 1,000). EHDV-1 was used as an outgroup to root the tree. Accession numbers are as follows: EHDV-1 segment 2, D10767; EHDV-1 segment 6, X55782; SV 124 (EHDV-2) segment 2, L33818; No.2 (IBAV) segment 2, AB030735; No.2 segment 6, AB030736.
|
|
|
|---|
The ORFs of segments 3 and 7 of all the strains of IBAV comprised 2,700 and 1,050 nucleotides, respectively. The variations of segments 3 and 7 within the strains of IBAV and CSIRO439 ranged from 92.8 to 93.7% for segment 3 and from 92.6 to 94.1% for segment 7, whereas the variations between segments 3 and 7 of IBAV together with CSIRO439 and those of SV 124 of EHDV-2 were both approximately 80% (Table 2). The deduced amino acid sequences of VP3 and VP7 were relatively conserved at more than 99% for the IBAV strains, except for the VP7 sequences of strain No.2. The amino acid sequence identities for VP7 between IBAVs and CSIRO439 ranged from 96.0 to 98.9%. The VP3 and VP7 amino acid sequence identities between IBAVs and SV 124 of EHDV-2 ranged from 95.3 to 95.9% for VP3 and from 94.6 to 97.1% for VP7. Similar lower identities between the VP3 and VP7 proteins of CSIRO439 and those of SV 124 were observed.
|
View this table: [in a new window] |
TABLE 2. Nucleotide and deduced amino acid sequence identities of segments 2, 6, 3, and 7 among IBAV and EHDV strains
|
A remarkable genomic variation was observed on the nucleotide sequence of segment 2. In contrast with the identities (between 94.4 and 98.5%) of segments 2 of strains No.2, Ibaraki-2, Y87061, and CSIRO439, the identities of segments 2 of KSB-14/E/97 and these strains were less than 70%. Poor identities of the deduced amino acid sequence of KSB-14/E/97 VP2 with those of the other strains of IBAV and EHDV-2 were observed, and the identities had lower values than those between the VP2 sequences of SV 124 and the other strains of IBAV.
Phylogenetic analysis. The phylogenetic trees of VP2, VP3, VP5, and VP7, constructed by using the deduced amino acid sequences encoded by the complete ORFs of segments 2, 3, 6 and 7, were individually compared. In the tree based on the VP3 sequence data (Fig. 1) there were two major clusters. SV 124 of EHDV-2 constituted an individual cluster, and the other cluster was composed of strains of IBAV and CSIRO439 of EHDV-2. Because strain SV 124 was isolated in North America and the other strains were isolated in Asia and Australia, the separation of these clusters seemed to correspond to the geographical distribution of the virus strains. The phylogenetic tree of VP7 resembled that of VP3 in its cluster construction. The clusters in the tree of VP7 were divided into three groups by geographical distribution and might define their geographical gene pool.
The phylogenetic trees of VP5 and VP2 were rather complicated. In the tree of VP5, strain KSB-14/E/97 was separated into a separate cluster (Fig. 2). However, strain CSIRO439 was in the same cluster as the existing strains of IBAV. The phylogenetic difference of strain KSB-14/E/97 was demonstrated more clearly in the tree of VP2. Strains No.2, Ibaraki-2, and Y87061 of IBAV were all in the same cluster, which was distinguishable both from the cluster of KSB-14/E/97 and from the cluster of SV 124 (Fig. 2).
|
|
|---|
The virus strains selected differed in geographical origin: the IBAV isolate was from Japan, and one EHDV-2 strain was from Australia, and the other was from North America. We compared sequences of segments 3 and 7, encoding serogroup-specific proteins VP3 and VP7, respectively, from the EHDV serogroup, to one another. The deduced amino acid sequences of the proteins were conserved, having 5% variability, between IBAV and EHDV. In the nucleotide sequences, there is a 20% variability between IBAV and North American EHDV, while there is a less than a 7% difference between IBAV and Australian EHDV. The phylogenetic analysis of VP3 showed two distinct clusters between the strains of IBAV and Australian EHDV-2. However, the KSB-14/E/97 strain of IBAV formed a single cluster with the CSIRO439 strain of EHDV-2. The result suggested that IBAV and Australian EHDV share a genetic pool. However, within the EHDV serogroup, phylogenetic analysis of segment 7 produced a more accurate geographical separation than analysis of segment 3, and so North American EHDV might be distinguishable from IBAV and Australian EHDV. Gould et al. demonstrated that a feature of the phylogenic trees based on VP3 within the EHDV and BTV serogroups could be clearly delineated on the basis of their geographical origins (11, 12, 25), and phylogenetic analyses revealed that the BTV segment 7 genes do not reflect the geographic origins of the virus strains (2, 32). Unlike the phylogenetic analyses of BTV, the phylogenetic analysis of IBAV and EHDV based on VP7 showed more clearly a geographical distribution. The VP7 has been identified as a protein that binds to vector cells (30, 33). The results of the phylogenetic analysis based on VP7 in this experiment, which divided the viruses into three groups (the North American and Australian EHDV-2 and the IBAV groups), might indicate viral adaptability to vector insects.
The nucleotide sequences of segments 6 from the strains of IBAV and CSIRO439 were relatively divergent, as were the sequences of segments 3 and 7. The identities of the deduced amino acid sequences of VP5 among these viruses were more than 98.1% and might bring the same antigenicity. Although the biological functions of VP5 are not well understood, VP5 was involved in the determination of virus serotype, possibly by imposing conformational constraints on VP2. The higher identities of VP5s from KSB-14/E/97 and the other strains of IBAV and Australian EHDV-2 might give a consistent explanation of the partial cross-reactivity of KSB-14/E/97 to other viruses in the cross-neutralization tests.
Segment 2 displayed remarkable genetic heterogeneity within the EHDV serogroup. The nucleotide sequence identities of segments 2 from KSB-14-E/97 and the existing strains of IBAV were less than 70%, and those of segments 2 from the existing strains of IBAV and the North American EHDV were lower than that. VP2 is the main neutralization antigen of the virus and is the determinant of virus serotype. The low identities of VP2 of KSB-14/E/97 might result in low reactivity in cross-neutralization tests with existing strains of IBAV. The phylogenetic analysis of VP2 and VP5 revealed that KSB-14/E/97 was segregated in a cluster distinct from that of IBAV and Australian EHDV-2. This phylogenetic distinction was remarkable in the tree of VP2. The greatest antigenic and sequence variations are also those exhibited by VP2s of BTV and EHDV (4, 8, 13). It is considered that such low nucleotide identity is the result of reassortment rather than the result of cumulative genetic mutation. The findings of the great variation of VP2, together with the findings of the variations of VP3, VP5, and VP7 of KSB-14/E/97, suggested that this new variant of IBAV was generated by a reassortment event in segment 2 and by mutations, including base substitution, deletion, and insertion, in segments 3, 6, and 7.
In natural conditions, reassortment occurs given the existence of at least two different strains present at the same time, in the same area, and in the same affected host. Although we could not estimate the other parental virus, the new variant of IBAV may have been generated in nature by genetic reassortment between IBAV and another orbivirus.
Our data analyzing four structural protein genes indicated that IBAV was clearly segregated from the North American EHDV and showed a slight difference from the Australian EHDV in nucleotides and amino acid sequences. It may be possible to regard IBAV and Australian EHDV-2 as nearly identical topotypes. Little is known about the mechanisms that induce such a variation between IBAV and Australian EHDV. When orbivirus infects host cells, VP2 is responsible for the attachment to mammalian cells, while VP7 mediates the attachment to insect cells. Therefore, vertebrate and invertebrate hosts might play important roles in the genetic variations of VP2 and VP7, respectively.
In recent years, epidemics caused by arboviruses have tended to spread over large areas due to the increasing activity of vector insects. Changes in the population of vector insects might be influencing the evolution of the viruses. Further study is needed to investigate the mutation of the virus.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»