Previous Article | Next Article 
Journal of Virology, September 2007, p. 10188-10194, Vol. 81, No. 18
0022-538X/07/$08.00+0 doi:10.1128/JVI.00337-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
Evidence of Intragenic Recombination in G1 Rotavirus VP7 Genes
Tung Gia Phan,1
Shoko Okitsu,1
Niwat Maneekarn,2 and
Hiroshi Ushijima1*
Department of Developmental Medical Sciences, Institute of International Health, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan,1
Department of Microbiology, Faculty of Medicine, Chiang Mai University, Chiang Mai, Thailand2
Received 15 February 2007/
Accepted 25 June 2007

ABSTRACT
The G1 rotavirus is the most widespread genotype causing acute
gastroenteritis in children. In an attempt to investigate the
occurrence of intragenic recombination, 131 complete coding
region sequences of VP7 genes of the G1 rotaviruses in GenBank
were examined. Three hitherto-unreported intragenic recombinant
rotaviruses were identified. It was noteworthy that two different
types (interlineage and intersublineage) of intragenic recombination
in rotaviruses were also found. This is the first report to
demonstrate the existence of intragenic recombinations between
interlineage and intersublineage in G1 rotaviruses.

TEXT
Acute gastroenteritis is a significant cause of morbidity and
mortality of children in both developed and developing countries.
Despite much progress in the understanding of the pathogenesis
of the disease and its management with the widespread use of
oral rehydration therapies, acute gastroenteritis consistently
ranks as one of the top causes of death worldwide (
9,
10,
20).
The rotaviruses, which comprise a genus in the family
Reoviridae,
are icosahedral in appearance. With VP4 spikes, the rotavirus
is about 110 nm in diameter. The virion of this virus is a nonenveloped,
triple-layered capsid containing 11 segments of double-stranded
RNA genome. The rotavirus genome encodes six structural and
six nonstructural proteins (
3). This virus is estimated to be
responsible for 111 million episodes of diarrhea requiring only
home care, 25 million clinic visits, 2 million hospitalizations,
and 325,000 to 592,000 deaths every year in children under five
years old (
12). Over past decades, G1 rotaviruses have been
the most widespread genotype causing acute gastroenteritis in
children from many countries covering all continents of the
world (
18). Nucleotide substitution and genomic reassortment
have been proposed to be the most important mechanisms of rotavirus
evolution in nature (
3,
4,
18). The rapidly increasing detection
of G1 rotavirus, in association with the genetic heterogeneity,
raises intriguing questions such as whether rotavirus evolution
is driven by intragenic recombination. Thus, the objective of
this study was to assess the occurrence of intragenic recombination
in the VP7 genes of the G1 rotaviruses.
A total of 131 sequences of the G1 rotaviruses, including our 36 sequence data from China, Japan, and Vietnam (15, 22, 26), which did not include any gaps in the alignment for the entire coding region of VP7 genes, were collected from GenBank. Sequence alignment was performed using CLUSTAL X (21). Phylogenetic trees with 100 bootstrap replicates of the nucleotide alignment datasets were generated using the neighbor-joining method (17). Genetic distance was calculated using Kimura's two-parameter method (PHYLIP) (6). SimPlot was used to detect recombinant rotavirus sequences as well as the breakpoints (7).
All 131 sequences of the G1 rotaviruses in this study were classified into different lineages and sublineages (Fig. 1) according to the recent G1 rotavirus classification scheme in which the nucleotide homology of rotavirus strains within each sublineage ranged from 98% to 100%, indicating a genetic difference of only less than 2% among them; the nucleotide sequence divergence between sublineages within the same lineage was from 3% to 4%; and sequence variation among strains between lineages was considerably higher, ranging from 5% to 16% (15). Interestingly, by using SimPlot, we found the VP7 genes of three G1 rotavirus strains to have intragenic recombinations between interlineage and intersublineage.
Ban-59.
The G1 rotavirus strain Ban-59 (U26366) was isolated from an
infant with acute gastroenteritis in Bangladesh during 1988,
and only this strain was assigned into lineage VII (
15). Figure
2A shows evidence of the novel recombinant G1 rotavirus bearing
different lineage sequence when the nucleotide sequence of strain
Ban-59 was compared with that of strain JP421 belonging to lineage
IV using SimPlot. The recombination breakpoints were estimated
at positions 469 and 649. Before position 469 and after position
649, the identities of Ban-59 and JP421 were distinctly different,
ranging from 90% to 93%. From positions 469 to 649, their identities
were extremely high (100%). In contrast, the examination of
the sequences for nucleotides 469 to 649 among G1 rotavirus
lineages revealed the low identities, ranging from 82% to 96%.
Of note, this region contained antigenic regions B (nucleotides
474 to 506) and E (nucleotides 615 to 620). Figure
2B also revealed
that Ban-59 clustered into different lineages when the different
part of VP7 gene-based grouping was performed. From nucleotides
49 to 468 and from nucleotides 650 to 1029, the lineage VII
of Ban-59 remained. However, Ban-59 was classified into lineage
IV. This kind of phenomenon was recognized as intragenic recombination
between interlineage.
Strain D.
More than 10 years before the isolation of strain Ban-59, strain
D (AB118022) was identified in the United States in 1974. Strain
D, together with three Brazilian strains isolated in 1991 and
1992 (Brazil-4, Brazil-5, and Brazil-6), formed sublineage IIIa
(
15). Using SimPlot, we determined the recombination event and
the breakpoint at position 449 in the VP7 gene when strain D
was compared with the reference strain Brazil-5, belonging to
sublineage IIIa, and with the reference strain Egypt-8, belonging
to sublineage IIIb (Fig.
3A). In comparison with Egypt-8, strain
D shared a high level of nucleotide identity (98%) in the region
from nucleotide 49 to nucleotide 448 and a lower level of the
nucleotide identity (97%) in the region from nucleotides 449
to 1029. In contrast, strain D shared a low level of nucleotide
identity (96%) in the region from nucleotides 49 to 448 and
a high level of the nucleotide identity (98%) in the region
from nucleotides 449 to 1029 with Brazil-5. Figure
3B demonstrated
that strain D clustered into different sublineages when the
different part of VP7 gene-based grouping was performed. Obviously,
strain D was recognized as the intragenic recombinant between
intersublineage IIIa and intersublineage IIb. It was found that
the recombinant region (nucleotides 49 to 448) contained antigenic
region A (nucleotides 309 to 353).
Russia-1407.
The G1 rotavirus strain Russia-1407 (S83903) was originally
detected in Russia, and it was assigned to sublineage IIe (
15).
When the nucleotide sequence of Russia-1407 was compared with
that of Oh-64 by using SimPlot, an apparent site of genetic
recombination was found at position 549 in the VP7 gene (Fig.
4A). Before this site, the identities of Russia-1407 and Oh-64
were rather low (only 96%). After this site, however, they were
highly similar (99%). The results demonstrated that the nucleotide
sequences of regions from positions 49 to 548 in these two strains
were rather different, but their sequences from positions 549
to 1029 were identical. The phylogenetic trees also supported
the SimPlot results. Figure
4B showed that Russia-1407 was grouped
into two distinct sublineages, IIa and IIe, according to its
low and high identities of nucleotide sequences to Oh-64. Taken
together, the results clearly indicated that Russia-1407 was
the novel intragenic recombinant of intersublineage. It was
found that the recombinant region (nucleotides 549 to 1029)
contained antigenic regions C (nucleotides 672 to 713), D (nucleotides
921 to 925), E (nucleotides 615 to 920), and F (nucleotides
747 to 776).
G1 rotaviruses have been reported as the most common genotype
in many countries (
8,
11,
18). Even though many studies have
conducted surveillance on rotavirus infection in various countries,
reports documenting intragenic recombination in rotaviruses
are limited. To date, only three intragenic recombinant rotaviruses,
CHW17 and CH55 from China (recombination between G1 and G3)
and ArgRes1723 from Argentina (recombination between sublineages
Ib and Ic within G4), were detected (
14,
19). Remarkably, we
found three novel intragenic recombinant rotaviruses, accounting
for 2.3% (3 of 131). This seems to be a high frequency for rotaviruses,
which are not generally thought to undergo intragenic recombination.
Three novel intragenic recombinants were found in three different
continents, including Asia (Ban-59), Europe (Russia-1407), and
North America (D), where the prevalence of G1 rotavirus infection
was very high, ranging from 60.3% to 73.7% (
18). However, among
three previously reported intragenic recombinants, two were
found in Asia (CHW17 and CH55) and one was found in South America,
where the prevalence of G3 and G4 rotavirus infection was low
(only 6.5% and 8.8.3%, respectively) (
18). Obviously, the prevalence
of intragenic recombinants in G3 and G4 rotaviruses is higher
than that in G1 rotavirus.
Phylogenetic analysis showed that these rotavirus strains clustered into different lineages and sublineages within the G1 genotype when gene-based grouping was performed for a different region of VP7. These observations suggested that there are two distinct kinds of intragenic recombination, interlineages and intersublineages. These results are noteworthy because this is the first report, to the best of our knowledge, showing intragenic recombination in interlineages and intersublineages of G1 rotaviruses. These phenomenons improve our current knowledge of the G1 rotavirus evolution as well as of the origin of genetic heterogeneity. Altogether, the recombination in rotaviruses could greatly affect phylogenetic groupings and confuse molecular epidemiological studies. Among three previously described intragenic recombinants, strains CHW17 and CH55 appeared to result from a double-crossover event (19) and strain ArgRes1723 appeared to result from a single-crossover event (14). Consistent with these findings, the single and double crossovers were also found in two novel intragenic recombinants (strain D and strain Russia-1407) and in one novel intragenic recombinant (strain Ban-59) detected in the study, respectively.
To date, two rotavirus vaccines, Rotarix and RotaTeq, have recently been released onto the market and licensed in more than 30 countries (1, 13, 16, 23, 24). Rotarix and RotaTeq have been proven to have significant clinical efficacies against G1 rotavirus gastroenteritis, with efficacy values of 96% and 95%, respectively (1, 23, 24). However, the occurrence of the G1 rotavirus gastroenteritis after immunization in children is possible. And these G1 rotavirus strains were different from the vaccine viruses (16, 23, 24). Therefore, the evidence of an intragenic recombination event in these cases should be investigated. Obviously, the intragenic recombination in the G1 rotaviruses in this study is important to aid the explanation of vaccine failure because it could produce immunity escape through exchanging antigenic regions between different G1 rotavirus strains with different antigenicities. Even antigenic regions were located in the recombinant regions of Ban-59, D, and Russia-1407; however, direct inspection of the alignment of the deduced sequences of antigenic regions of VP7 revealed that these antigenic regions are conserved among them and among the reference rotavirus strains used in the present study (15). It is likely that structural variation of the VP7 protein within a genotype develops by successive mutations in the gene (2). Therefore, the intragenic recombination in VP7 genes of three strains might induce the conformational changes of G1 rotavirus VP7 protein, which probably led to the changing antigenicity of these intragenic recombinant strains. To date, a study on localization of amino acids involved in conformational changes of rotavirus VP7 structure is not available (15). The identification of these amino acids is of significance and should be investigated by further studies.
Although the crucial contributions of genetic reassortment and nucleotide substitution in producing antigenic variants for rotaviruses (3, 4) have been recognized, the intragenic recombination events in a few rotavirus strains were reported. Several studies showed the mixed infections with two different serotypes or with different variants belonging to the same serotype in one individual, therefore increasing the probability of reassortment or intragenic recombination in the evolutionary pathway of rotaviruses (3-5). RNA intragenic recombination is a mechanism for virus evolution (25). The primary mechanism involved in intragenic recombination in RNA viruses is the copy-choice model in which intragenic recombination is known to depend on various immunological and intracellular constraints, such as (i) successful coinfection of the host and in a single cell by two parental strains, (ii) efficient replication of parental viral genomes with template switching, and (iii) adaptation to different environments to be transmitted (25). The observation of the rotavirus intragenic recombinant strains detected in the study probably underscored that these strains theoretically might fulfill all prerequisites for their intragenic recombination. However, it is unclear that a copy-choice mechanism could work in rotaviruses where all replication occurs in particles and the replication complexes are anchored and separated from each other. The predicted amino acid sequences, including regions adjacent to the suggested breakpoints, were further analyzed. No same amino acid breakpoints were identified since the nucleotide breakpoints at positions 469 and 649 in Ban-59 are involved in coding for the amino acids leucine (position 141) and glutamine (position 201), respectively, and the nucleotide breakpoints at positions 449 and 549 in D and Russia-1407 are involved in coding for the amino acids tyrosine (position 134) and proline (position 167), respectively. Interestingly, the alignment of complete VP7 protein using CLUSTAL X showed the highly conserved region from amino acid 124 to amino acid 210 in G1 rotaviruses (data not shown). More interestingly, all of four amino acid breakpoints of three recombinant G1 rotaviruses were located in this region. Altogether, this highly conserved region in the VP7 gene was a favorable condition for recombination in G1 rotaviruses. However, further research should be conducted to elucidate the molecular mechanism of intragenic recombination in rotaviruses.
In conclusion, our results have described the genetic characterization of novel intragenic recombinant rotaviruses and have increased the evidence of worldwide distribution of intragenic recombinant rotaviruses. Intragenic recombination events in rotaviruses are not completely understood; further research should be conducted to investigate whether intragenic recombination can be potentially dangerous for host species, and it likely limits the rotavirus control programs and has major implications for rotavirus vaccine design.

ACKNOWLEDGMENTS
This study was supported by grants-in-aid from the Ministry
of Education, Culture, Sports, Sciences, and Technology, Japan,
and the Ministry of Health, Labor, and Welfare, Japan.

FOOTNOTES
* Corresponding author. Mailing address: Department of Developmental Medical Sciences, Institute of International Health, Graduate School of Medicine, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033. Phone: 81-3-5841-3590. Fax: 81-3-5841-3629. E-mail:
ushijima{at}m.u-tokyo.ac.jp 
Published ahead of print on 3 July 2007. 

REFERENCES
1 - Bernstein, D. I., and R. L. Ward. 2006. Rotarix: development of a live attenuated monovalent human rotavirus vaccine. Pediatr. Ann. 35:38-43.[Medline]
2 - Blackhall, J., R. Bellinzoni, N. Mattion, M. K. Estes, J. I. La Torre, and G. Magnusson. 1992. A bovine rotavirus serotype 1: serologic characterization of the virus and nucleotide sequence determination of the structural glycoprotein VP7 gene. Virology 189:833-837.[CrossRef][Medline]
3 - Estes, M. K. 2001. Rotaviruses and their replication, p. 1747-1785. In D. M. Knipe, P. M. Howley, et al. (ed.), Fields virology, 4th ed. Lippincott Williams & Wilkins, Philadelphia, PA.
4 - Gouvea, V., and M. Brantly. 1995. Is rotavirus a population of reassortants? Trends Microbiol. 3:159-162.[CrossRef][Medline]
5 - Jain, V., B. K. Das, M. K. Bhan, R. I. Glass, and J. R. Gentsch. 2001. Great diversity of group A rotavirus strains and high prevalence of mixed rotavirus infections in India. J. Clin. Microbiol. 39:3524-3529.[Abstract/Free Full Text]
6 - Kimura, M. 1980. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16:111-120.[CrossRef][Medline]
7 - Lole, K. S., R. C. Bollinger, R. S. Paranjape, D. Gadkari, S. S. Kulkarni, N. G. Novak, R. Ingersoll, H. W. Sheppard, and S. C. Ray. 1999. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 73:152-160.[Abstract/Free Full Text]
8 - Maneekarn, N., and H. Ushijima. 2000. Epidemiology of rotavirus infection in Thailand. Pediatr. Int. 42:415-421.[CrossRef][Medline]
9 - Mulholland, E. K. 2004. Global control of rotavirus disease. Adv. Exp. Med. Biol. 549:161-168.[Medline]
10 - Murray, C. J., and A. D. Lopez. 1997. Mortality by cause for eight regions of the world: global burden of disease study. Lancet 349:1269-1276.[CrossRef][Medline]
11 - Okitsu-Negishi, S., T. A. Nguyen, T. G. Phan, and H. Ushijima. 2004. Molecular epidemiology of viral gastroenteritis in Asia. Pediatr. Int. 46:245-252.[CrossRef][Medline]
12 - Parashar, U. D., E. G. Hummelman, M. A. Miller, and R. I. Glass. 2003. Global illness and deaths caused by rotavirus disease in children. Emerg. Infect. Dis. 9:565-572.[Medline]
13 - Parashar, U. D., and R. I. Glass. 2006. Public health. Progress toward rotavirus vaccines. Science 12:851-852.
14 - Parra, G. I., K. Bok, M. Martinez, and J. A. Gomez. 2004. Evidence of rotavirus intragenic recombination between two sublineages of the same genotype. J. Gen. Virol. 85:1713-1716.[Abstract/Free Full Text]
15 - Phan, T. G., P. Khamrin, T. D. Quang, S. K. Dey, S. Takanashi, S. Okitsu, N. Maneekarn, and H. Ushijima. 2007. Detection and genetic characterization of group A rotavirus strains circulating among children with acute gastroenteritis in Japan. J. Virol. 81:4645-4653.[Abstract/Free Full Text]
16 - Ruiz-Palacios, G. M., I. Perez-Schael, F. R. Velazquez, H. Abate, T. Breuer, S. C. Clemens, B. Cheuvart, F. Espinoza, P. Gillard, B. L. Innis, Y. Cervantes, A. C. Linhares, P. Lopez, M. Macias-Parra, E. Ortega-Barria, V. Richardson, D. M. Rivera-Medina, L. Rivera, B. Salinas, N. Pavia-Ruz, J. Salmeron, R. Ruttimann, J. C. Tinoco, P. Rubio, E. Nunez, M. L. Guerrero, J. P. Yarzabal, S. Damaso, N. Tornieporth, X. Saez-Llorens, R. F. Vergara, T. Vesikari, A. Bouckenooghe, R. Clemens, B. De Vos, M. O'Ryan, and the Human Rotavirus Vaccine Study Group. 2006. Safety and efficacy of an attenuated vaccine against severe rotavirus gastroenteritis. N. Engl. J. Med. 354:11-22.[Abstract/Free Full Text]
17 - Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425.[Abstract]
18 - Santos, N., and Y. Hoshino. 2005. Global distribution of rotavirus serotypes/genotypes and its implication for the development and implementation of an effective rotavirus vaccine. Rev. Med. Virol. 15:29-56.[CrossRef][Medline]
19 - Suzuki, Y., T. Gojobori, and O. Nakagomi. 1998. Intragenic recombinations in rotaviruses. FEBS Lett. 8:183-187.
20 - Thapar, N., and I. R. Sanderson. 2004. Diarrhoea in children: an interface between developing and developed countries. Lancet 363:641-653.[CrossRef][Medline]
21 - Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 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]
22 - Trinh, Q. D., T. A. Nguyen, T. G. Phan, P. Khamrin, H. Yan, L. P. Hoang, N. Maneekarn, Y. Li, F. Yagyu, S. Okitsu, and H. Ushijima. 2007. Sequence analysis of the VP7 gene of human rotavirus G1 isolated in Japan, China, Thailand, and Vietnam in the context of changing distribution of rotavirus G-types. J. Med. Virol. 79:1009-1016.[CrossRef][Medline]
23 - Vesikari, T., A. Karvonen, L. Puustinen, S. Q. Zeng, E. D. Szakal, A. Delem, and B. De Vos. 2004. Efficacy of RIX4414 live attenuated human rotavirus vaccine in Finnish infants. Pediatr. Infect. Dis. J. 23:937-943.[CrossRef][Medline]
24 - Vesikari, T., D. O. Matson, P. Dennehy, P. Van Damme, M. Santosham, Z. Rodriguez, M. J. Dallas, J. F. Heyse, M. G. Goveia, S. B. Black, H. R. Shinefield, C. D. Christie, S. Ylitalo, R. F. Itzler, M. L. Coia, M. T. Onorato, B. A. Adeyi, G. S. Marshall, L. Gothefors, D. Campens, A. Karvonen, J. P. Watt, K. L. O'Brien, M. L. DiNubilem, H. F. Clark, J. W. Boslego, P. A. Offit, P. M. Heaton, and the Rotavirus Efficacy and Safety Trial (REST) Study Team. 2006. Safety and efficacy of a pentavalent human-bovine (WC3) reassortant rotavirus vaccine. N. Engl. J. Med. 354:23-33.[Abstract/Free Full Text]
25 - Worobey, M., and E. C. Holmes. 1999. Evolutionary aspects of recombination in RNA viruses. J. Gen. Virol. 80:2535-2543.[Free Full Text]
26 - Xin, K. Q., S. Morikawa, Z. Y. Fang, A. Mukoyama, K. Okuda, and H. Ushijima. 1993. Genetic variation in VP7 gene of human rotavirus serotype 1 (G1 type) isolated in Japan and China. Virology 197:813-816.[CrossRef][Medline]
Journal of Virology, September 2007, p. 10188-10194, Vol. 81, No. 18
0022-538X/07/$08.00+0 doi:10.1128/JVI.00337-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Martinez-Laso, J., Roman, A., Rodriguez, M., Cervera, I., Head, J., Rodriguez-Avial, I., Picazo, J. J.
(2009). Diversity of the G3 genes of human rotaviruses in isolates from Spain from 2004 to 2006: cross-species transmission and inter-genotype recombination generates alleles. J. Gen. Virol.
90: 935-943
[Abstract]
[Full Text]
-
Hon, C.-C., Lam, T. T.-Y., Yip, C.-W., Wong, R. T.-Y., Shi, M., Jiang, J., Zeng, F., Leung, F. C.-C.
(2008). Phylogenetic evidence for homologous recombination within the family Birnaviridae. J. Gen. Virol.
89: 3156-3164
[Abstract]
[Full Text]