1 full-genome characterization of g8p[8] rotavirus emerged among

19
1 Full-genome characterization of G8P[8] rotavirus emerged among children with diarrhea in 1 Croatia, in 2006 2 3 Roberto Delogu 1‡ , Alessandra Lo Presti 2‡ , Franco Maria Ruggeri 3# , Eleonora Cella 2 , Marta Giovanetti 2 , 4 Massimo Ciccozzi 2 , Suncanica Ljubin Sternak 4 , Suzana Bukovski-Simonoski 5 , Amarela Lukic-Grlic 6 , 5 Giovanni Ianiro 1 , Lucia Fiore 1 . 6 7 1 National Center for Immunobiologicals Research and Evaluation, 2 Deptartment of Infectious, 8 Parasitic and Immune-Mediated Diseases, and 3 Department of Veterinary Public Health and Food 9 Safety, Istituto Superiore di Sanità, Rome, Italy; 4 Croatian National Institute of Public Health, and 10 School of Medicine University of Zagreb, Zagreb, Croatia; 5 University Hospital for Infectious Diseases 11 “Dr Fran Mihaljevic”, Zagreb, Croatia; 6 Zagreb University Children’s Hospital, and School of 12 Medicine University of Zagreb, Zagreb, Croatia. 13 14 RDL and ALP have equally contributed to the work 15 16 Running title: Full-genome sequencing of G8P[8] rotavirus in Croatia 17 18 # Corresponding author: 19 Franco M Ruggeri, PhD 20 Dept. Veterinary Public Health & Food Safety 21 Istituto Superiore di Sanità, 22 Viale Regina Elena, 299, 00161 Rome, Italy 23 Phone: +39 06 4990 2980 24 Copyright © 2013, American Society for Microbiology. All Rights Reserved. J. Clin. Microbiol. doi:10.1128/JCM.00396-13 JCM Accepts, published online ahead of print on 20 February 2013 on April 3, 2018 by guest http://jcm.asm.org/ Downloaded from

Upload: dinhphuc

Post on 03-Feb-2017

215 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: 1 Full-genome characterization of G8P[8] rotavirus emerged among

1

Full-genome characterization of G8P[8] rotavirus emerged among children with diarrhea in 1

Croatia, in 2006 2

3

Roberto Delogu1‡, Alessandra Lo Presti2‡, Franco Maria Ruggeri3#, Eleonora Cella2, Marta Giovanetti2, 4

Massimo Ciccozzi2, Suncanica Ljubin Sternak4, Suzana Bukovski-Simonoski5, Amarela Lukic-Grlic6, 5

Giovanni Ianiro1, Lucia Fiore1. 6

7

1National Center for Immunobiologicals Research and Evaluation, 2Deptartment of Infectious, 8

Parasitic and Immune-Mediated Diseases, and 3Department of Veterinary Public Health and Food 9

Safety, Istituto Superiore di Sanità, Rome, Italy; 4Croatian National Institute of Public Health, and 10

School of Medicine University of Zagreb, Zagreb, Croatia; 5University Hospital for Infectious Diseases 11

“Dr Fran Mihaljevic”, Zagreb, Croatia; 6Zagreb University Children’s Hospital, and School of 12

Medicine University of Zagreb, Zagreb, Croatia. 13

14

‡RDL and ALP have equally contributed to the work 15

16

Running title: Full-genome sequencing of G8P[8] rotavirus in Croatia 17

18

#Corresponding author: 19

Franco M Ruggeri, PhD 20

Dept. Veterinary Public Health & Food Safety 21

Istituto Superiore di Sanità, 22

Viale Regina Elena, 299, 00161 Rome, Italy 23

Phone: +39 06 4990 2980 24

Copyright © 2013, American Society for Microbiology. All Rights Reserved.J. Clin. Microbiol. doi:10.1128/JCM.00396-13 JCM Accepts, published online ahead of print on 20 February 2013

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 2: 1 Full-genome characterization of G8P[8] rotavirus emerged among

2

Fax: +39 06 4938 7101 25

E-mail: [email protected] 26

27

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 3: 1 Full-genome characterization of G8P[8] rotavirus emerged among

3

Abstract 28

A 2006 epidemic G8P[8] rotavirus from Croatia was whole-genome sequenced, showing a Wa-like 29

genotype constellation. Its VP7 clustered with DS1-like G8 African rotaviruses, and a G8P[4] German 30

strain. Remaining genes clustered with contemporary Belgian G1P[8] rotaviruses, suggesting 31

reassortment between human G8 and G1P[8] rotaviruses, in Croatia or other European countries. 32

33

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 4: 1 Full-genome characterization of G8P[8] rotavirus emerged among

4

Rotavirus A (RVA) causes acute gastroenteritis in man and animals globally, killing 450,000 children 34

per year, mainly in developing countries (1). 35

RVA genome includes 11 double-stranded RNA segments, encoding six viral (VP1-4,6,7) and six non-36

structural proteins (NSP1-6) (2). Rotavirus outer layer protein VP7 and spike protein VP4 are 37

neutralization antigens, determining virus genotypes G (glycoprotein, VP7) and P (protease-sensitive, 38

VP4). At least 27 G and 37 P genotypes have been identified among human, animal and avian RVAs 39

(3-5). 40

Five genotypes (G1P[8], G2P[4], G3P[8], G4P[8], and G9P[8]) represent over 90% of global human 41

RVAs (6), that may explain the efficacy of both monovalent (P[8], G1) Rotarix and pentavalent (P[8], 42

G1-4) RotaTeq vaccines against rotavirus diarrhea worldwide (7). 43

Serotypes other than G1-4, and strains evolving by mutation or genomic reassortment between human 44

and/or animal RVAs represent an emerging threat for children (8-10). Millions of people are exposed to 45

animal rotaviruses worldwide, and several uncommon human genotypes (G6, G8, G12, P3[9], P5A[3]) 46

closely resemble animal RVA strains (11). 47

Molecular epidemiology and viral typing can help monitor emergence of novel RVAs, confirming 48

current vaccine efficacy against unconventional genotypes. 49

Reassortment events can involve any dsRNA segments, and G and P typing may be insufficient to 50

investigate RVA origin and inter-species transmission. Genotyping based on sequencing all 11 genome 51

segments (3, 4) is very helpful, distinguishing 9 VP1 (R), 9 VP2 (C), 8 VP3 (M), 35 VP4 (P), 16 VP6 52

(I), 27 VP7 (G), 16 NSP1 (A), 9 NSP2 (N), 12 NSP3 (T), 14 NSP4 (E), and 11 NSP5 (H) genotypes. 53

Bovine-like G8 rotavirus was first reported in an Indonesian child (12), and further G8 RVA human 54

cases occurred worldwide, including industrialized countries (13-21). Using full-genome sequencing, 55

G8P[8] African RVAs were shown to involve reassortment between at least four human, swine, and 56

bovine strains (22). 57

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 5: 1 Full-genome characterization of G8P[8] rotavirus emerged among

5

This paper reports the whole genome characterization of a G8P[8] human rotavirus strain emerged in 58

Croatia in 2006 (23). 59

Cases showed severe gastroenteritis, and included <5 year-old children from Croatian hospitals in 60

2005-2006. Stools were screened by Rotalex, and RVA characterized by RT-PCR (23). Agarose gel 61

purified amplicons were sequenced using RT-PCR primers (8, 18, 24, 25), assembling consensus 62

sequences for all gene segments (18, 22, 26). 63

Twelve datasets were built, including all NCBI sequences (http://www.ncbi.nlm.nih.gov/pubmed), 64

showing >89% similarity with query sequences (between 116 for VP2 and 127 for NSP4). 65

Sequences were aligned using ClustalX software (27), and edited by Bioedit software. The ModelTest 66

v3.0 (28) with hierarchical likelihood ratio test was used to select the best-fit models for sequence data 67

analysis. 68

Bayesian phylogenetic trees were constructed for VP7, VP4 and NSP4 datasets (29), using GTR+I+G 69

nucleotide substitution model for VP7, and HKY+G model for VP4/NSP4. 70

Markov Chain Monte Carlo searches were made (50x106 generations, tree-sampling every 5000 71

generation, and 10% burn-in fraction. Clade statistical support followed >0.90 posterior probability). 72

Maximum-likelihood (ML) phylogenetic trees of NSP1-3, NSP5-6, VP1-3 and VP6 were generated 73

with PAUP* 4.0 package (30), using GTR+I+G for VP6, NSP1, VP1; HKY+G model for NSP6/NSP5; 74

TrN+I+G for NSP2, NSP3, VP2 and VP3 datasets. Statistical robustness and reliability of branching 75

order within trees was confirmed by bootstrap analysis (1000 replicates), considering bootstrap >70% 76

as clades support. 77

Sequences were deposited in GenBank database (accession numbers JQ988894-JQ988904). 78

79

Initial sequencing of genes encoding, VP4 and NSP4 for 10/31 G8P[8] RVAs showed nearly complete 80

nucleotide identity suggesting circulation of a single G8P[8] strain in Croatia in 2005-6. The consensus 81

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 6: 1 Full-genome characterization of G8P[8] rotavirus emerged among

6

Croatia 2006 (CR2006) VP7 gene showed >98% nucleotide identity with other human G8 GenBank 82

strains, particularly from neighboring Slovenia, 2006, and Africa, 2000-2004. CR2006 VP4 and NSP4 83

recalled human G3P[8] strains reported in Italy, 2005 (31). To investigate Croatian G8P[8] origin, full-84

genome sequencing was performed, using RotaC (http://rotac.regatools.be) (32). CR2006 presented a 85

genotype G8-P[8]-I1-R1-C1-M1-A1-N1-T1-E1-H1. 86

Bayesian analysis identified CR2006 VP7 in a statistically supported cluster with African (Malawi, 87

Tunisia, Cameroon, Ethiopia), and European strains (Slovenia [Sl-885/06], Germany [GER1H-09]) 88

(Fig 1A). In VP4 and NSP4 trees (Fig 1B and 2), CR2006 also clustered with European and African 89

RVAs, including Hu/Slo/2006/G8P[8]/Sl-885/06, Hu/Italy/2005/G3P[8], 90

Hu/Ethiopia/2004/G8P[8]/ARN, Hu/Tunisia/2000-2004/G8P[8]/ARN. 91

CR2006 clustered with Wa-like genotype constellation Hu/Bel/2005-2009/G1P[8], 92

Hu/Ger/G3G12P[8]/GER126-08, and other global RVAs in NSP1-NSP6, VP1-VP3, and VP6 93

(Supplementary Fig S1-9). 94

G8P[8] rotavirus circulated abnormally in Croatia in 2006, reaching 15.8% of children hospitalized 95

with diarrhea (23). Full-genome sequencing of these viruses was performed to investigate links 96

between their emergence and genome characters. 97

The G8 rotavirus genotype, normally combined with P[1],[5],[11] or [21] types (11, 33, 34), is typical 98

of bovine strains, but has also been sporadically detected in humans since 1978 (35) associated with 99

different P-types (13, 16-18, 36-38). Human G8 strains normally exhibit P[4] or P[6] both in developed 100

countries and in Africa, where they are particularly prevalent (15-17, 26). Recently, full-genome 101

sequencing of a human G8P[8] rotavirus in Democratic Republic of Congo (DRC) in 2003 (18, 22) 102

suggested this strain to originate from reassortment between a bovine G8P[6] ancestor with three or 103

more different rotaviruses, including human Wa- and DS1-like strains. This strain exhibited a DS-1-104

like genotype constellation. 105

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 7: 1 Full-genome characterization of G8P[8] rotavirus emerged among

7

Differently, Croatian CR2006 G8P[8] presented the Wa-like constellation (3) G8-P[8]-I1-R1-C1-M1-106

A1-N1-T1-E1-H1, all genes showing close evolutionary relationship to prototype G1P[8] human strain 107

Wa, except the VP7 gene, related to African and Slovenian human G8 strains. Rotaviruses are usually 108

species-specific, and gene reassortment with human strains appears necessary for zoonotic adaptation 109

(26, 39-41). Characterization of epidemic CR2006 further supports this view. As VP7 plays a major 110

role in the immune response (42, 43), CR2006 emergence and epidemic spread may have involved 111

reassortment with common rotaviruses pre-existing in Croatia. 112

CR2006 was more frequent in younger children (23) than any other rotavirus genotypes co-circulating 113

in Croatia, 2006. This could be fortuitous, or suggest that maternal antibodies protected younger 114

children from CR2006 less than from other serotypes. Maternal immunity fades quickly (44, 45), and it 115

is conceivable that older Croatian children in 2006 were similarly unprotected hence susceptible to G8 116

or any other common rotavirus genotypes. 117

Phylogenetic analyses revealed that, except for the VP7 gene, CR2006 was more closely related to 118

humans G1P[8] strains identified between 2005-2009 in Belgium. Therefore, CR2006 origin likely 119

involved reassortment of a Wa-like rotavirus circulating in one or more countries of Europe, including 120

Belgium and possibly Croatia, with a co-circulating G8 virus. In 2006, a similar G8P[8] strain was 121

reported and partially sequenced in Slovenia (20), whereas no G8P[8] rotaviruses were detected in 122

other European countries. It is possible that the reassortment event leading to CR2006 may have 123

occurred in or near Croatia. Low VP7 correlation between CR2006 and G8 DS1-like rotaviruses 124

characterized in African DRC (18) indicates independent evolution of these G8 reassortants. 125

Conversely, CR2006 shared a similar gene 9 with human G8 strains from Northern/Central Africa in 126

the early 2000s, suggesting that the G8 gene followed immigration routes from Africa. A closely 127

similar VP7 was found in a G8P[4] human strain identified in Germany two years later (21), suggesting 128

that the CR2006 G8 gene participated to additional reassortments with European DS1-like rotaviruses. 129

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 8: 1 Full-genome characterization of G8P[8] rotavirus emerged among

8

The German strain was strictly correlated with the DRC G8 rotaviruses reported earlier, except for VP4 130

(18), and may represent a further evolution of G8P[8] rotaviruses imported from Africa. The VP7 of all 131

European G8 strains differed markedly from bovine G8 rotaviruses, contrary to a recent human-bovine 132

reassortment event within Europe. 133

Rotavirus genetic evolution is certainly influenced by immune pressure (46), that may favor G8 and 134

other uncommon rotavirus G-types spreading (19, 20, 47), as occurred for G9 strains originated from 135

swine (48, 49). International transmission of genotypes is facilitated by globalization and travelling, 136

urging worldwide rotavirus strain surveillance to anticipate emergence of novel strains and 137

animal/human reassortants, and monitor effectiveness of current human rotavirus vaccines. 138

139

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 9: 1 Full-genome characterization of G8P[8] rotavirus emerged among

9

Acknowledgments 140

Grant sponsor: Ministry of Health, Italy (to LF CCM ‘‘Epidemiologia molecolare di rotavirus in eta` 141

pediatrica in Italia. Creazione di una rete di sorveglianza per monitorare la diffusione e l’evoluzione di 142

genotipi virali’’ and to FMR Strategic Program: ‘‘Preparedness and response to emerging zoonoses and 143

exotic viral infection through an integrated medical and veterinary approach’’); Grant sponsor: 144

ISS/NIH (to FMR ‘‘Molecular and antigenic evolution of rotavirus strains of human and animal 145

origin’’); Grant sponsor: EuroRotaNet (http://www.eurorota.net). 146

147

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 10: 1 Full-genome characterization of G8P[8] rotavirus emerged among

10

References 148

1. Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Duque J, Parashar UD. 2012. 2008 149

estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the 150

introduction of universal rotavirus vaccination programmes: a systematic review and meta-151

analysis. The Lancet Infectious Diseases 12:136-141. 152

2. Estes MK, Kapikian AZ. 2007. Rotaviruses, p. 1917-1974. In Knipe DM, Howley PM, Griffin 153

DE, Lamb RA, Martin MA, Roizman B, Straus SE (ed.), Fields Virology, vol. 2. Kluwer / 154

Lippincott, Williams and Wilkins, Philadelphia, PA. 155

3. Matthijnssens J, Ciarlet M, Rahman M, Attoui H, Banyai K, Estes MK, Gentsch JR, 156

Iturriza-Gomara M, Kirkwood CD, Martella V, Mertens PP, Nakagomi O, Patton JT, 157

Ruggeri FM, Saif LJ, Santos N, Steyer A, Taniguchi K, Desselberger U, Van Ranst M. 158

2008. Recommendations for the classification of group A rotaviruses using all 11 genomic 159

RNA segments. Arch Virol 153:1621-1629. 160

4. Matthijnssens J, Ciarlet M, McDonald SM, Attoui H, Banyai K, Brister JR, Buesa J, 161

Esona MD, Estes MK, Gentsch JR, Iturriza-Gomara M, Johne R, Kirkwood CD, Martella 162

V, Mertens PP, Nakagomi O, Parreno V, Rahman M, Ruggeri FM, Saif LJ, Santos N, 163

Steyer A, Taniguchi K, Patton JT, Desselberger U, Van Ranst M. 2011. Uniformity of 164

rotavirus strain nomenclature proposed by the Rotavirus Classification Working Group 165

(RCWG). Arch Virol 156:1397-1413. 166

5. Trojnar E, Sachsenroder J, Twardziok S, Reetz J, Otto PH, Johne R. 2013. Identification 167

of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of 168

mammalian rotaviruses. The Journal of general virology 94:136-142. 169

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 11: 1 Full-genome characterization of G8P[8] rotavirus emerged among

11

6. Santos N, Hoshino Y. 2005. Global distribution of rotavirus serotypes/genotypes and its 170

implication for the development and implementation of an effective rotavirus vaccine. Rev Med 171

Virol 15:29-56. 172

7. Jiang V, Jiang B, Tate J, Parashar UD, Patel MM. 2010. Performance of rotavirus vaccines 173

in developed and developing countries. Human vaccines 6:532-542. 174

8. Iturriza-Gomara M, Isherwood B, Desselberger U, Gray J. 2001. Reassortment in vivo: 175

driving force for diversity of human rotavirus strains isolated in the United Kingdom between 176

1995 and 1999. J Virol 75:3696-3705. 177

9. Maunula L, Von Bonsdorff CH. 2002. Frequent reassortments may explain the genetic 178

heterogeneity of rotaviruses: analysis of Finnish rotavirus strains. J Virol 76:11793-11800. 179

10. Gentsch JR, Laird AR, Bielfelt B, Griffin DD, Banyai K, Ramachandran M, Jain V, 180

Cunliffe NA, Nakagomi O, Kirkwood CD, Fischer TK, Parashar UD, Bresee JS, Jiang B, 181

Glass RI. 2005. Serotype diversity and reassortment between human and animal rotavirus 182

strains: implications for rotavirus vaccine programs. J Infect Dis 192 Suppl 1:S146-159. 183

11. Martella V, Banyai K, Matthijnssens J, Buonavoglia C, Ciarlet M. 2010. Zoonotic aspects 184

of rotaviruses. Vet Microbiol 140:246-255. 185

12. Hasegawa A, Inouye S, Matsuno S, Yamaoka K, Eko R, Suharyono W. 1984. Isolation of 186

human rotaviruses with a distinct RNA electrophoretic pattern from Indonesia. Microbiol 187

Immunol 28:719-722. 188

13. Holmes JL, Kirkwood CD, Gerna G, Clemens JD, Rao MR, Naficy AB, Abu-Elyazeed R, 189

Savarino SJ, Glass RI, Gentsch JR. 1999. Characterization of unusual G8 rotavirus strains 190

isolated from Egyptian children. Arch Virol 144:1381-1396. 191

14. Santos N, Lima RC, Pereira CF, Gouvea V. 1998. Detection of rotavirus types G8 and G10 192

among Brazilian children with diarrhea. J Clin Microbiol 36:2727-2729. 193

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 12: 1 Full-genome characterization of G8P[8] rotavirus emerged among

12

15. Steele AD, Parker SP, Peenze I, Pager CT, Taylor MB, Cubitt WD. 1999. Comparative 194

studies of human rotavirus serotype G8 strains recovered in South Africa and the United 195

Kingdom. The Journal of general virology 80 ( Pt 11):3029-3034. 196

16. Cunliffe NA, Gentsch JR, Kirkwood CD, Gondwe JS, Dove W, Nakagomi O, Nakagomi T, 197

Hoshino Y, Bresee JS, Glass RI, Molyneux ME, Hart CA. 2000. Molecular and serologic 198

characterization of novel serotype G8 human rotavirus strains detected in Blantyre, Malawi. 199

Virology 274:309-320. 200

17. Adah MI, Nagashima S, Wakuda M, Taniguchi K. 2003. Close relationship between G8-201

serotype bovine and human rotaviruses isolated in Nigeria. J Clin Microbiol 41:3945-3950. 202

18. Matthijnssens J, Rahman M, Yang X, Delbeke T, Arijs I, Kabue JP, Muyembe JJ, Van 203

Ranst M. 2006. G8 rotavirus strains isolated in the Democratic Republic of Congo belong to 204

the DS-1-like genogroup. J Clin Microbiol 44:1801-1809. 205

19. Montenegro FM, Correia JB, Rodrigues Falbo A, Dove W, Nakagomi T, Nakagomi O, 206

Cuevas LE, Cunliffe NA, Hart CA. 2007. Anticipating rotavirus vaccines in Brazil: detection 207

and molecular characterization of emerging rotavirus serotypes G8 and G9 among children with 208

diarrhoea in Recife, Brazil. J Med Virol 79:335-340. 209

20. Steyer A, Poljsak-Prijatelj M, Bufon TL, Marcun-Varda N, Marin J. 2007. Rotavirus 210

genotypes in Slovenia: unexpected detection of G8P[8] and G12P[8] genotypes. J Med Virol 211

79:626-632. 212

21. Pietsch C, Petersen L, Patzer L, Liebert UG. 2009. Molecular characteristics of German 213

G8P[4] rotavirus strain GER1H-09 suggest that a genotyping and subclassification update is 214

required for G8. J Clin Microbiol 47:3569-3576. 215

22. Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, Palombo EA, 216

Iturriza-Gomara M, Maes P, Patton JT, Rahman M, Van Ranst M. 2008. Full genome-217

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 13: 1 Full-genome characterization of G8P[8] rotavirus emerged among

13

based classification of rotaviruses reveals a common origin between human Wa-Like and 218

porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol 82:3204-219

3219. 220

23. Tcheremenskaia O, Marucci G, De Petris S, Ruggeri FM, Dovecar D, Sternak SL, 221

Matyasova I, Dhimolea MK, Mladenova Z, Fiore L. 2007. Molecular epidemiology of 222

rotavirus in Central and Southeastern Europe. J Clin Microbiol 45:2197-2204. 223

24. Das BK, Gentsch JR, Cicirello HG, Woods PA, Gupta A, Ramachandran M, Kumar R, 224

Bhan MK, Glass RI. 1994. Characterization of rotavirus strains from newborns in New Delhi, 225

India. J Clin Microbiol 32:1820-1822. 226

25. Gentsch JR, Glass RI, Woods P, Gouvea V, Gorziglia M, Flores J, Das BK, Bhan MK. 227

1992. Identification of group A rotavirus gene 4 types by polymerase chain reaction. J Clin 228

Microbiol 30:1365-1373. 229

26. Esona MD, Geyer A, Page N, Trabelsi A, Fodha I, Aminu M, Agbaya VA, Tsion B, Kerin 230

TK, Armah GE, Steele AD, Glass RI, Gentsch JR. 2009. Genomic characterization of human 231

rotavirus G8 strains from the African rotavirus network: relationship to animal rotaviruses. J 232

Med Virol 81:937-951. 233

27. Thompson JD, Higgins DG, Gibson TJ. 1994. CLUSTAL W: improving the sensitivity of 234

progressive multiple sequence alignment through sequence weighting, position-specific gap 235

penalties and weight matrix choice. Nucleic Acids Res 22:4673-4680. 236

28. Posada D, Buckley TR. 2004. Model selection and model averaging in phylogenetics: 237

advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. 238

Syst Biol 53:793-808. 239

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 14: 1 Full-genome characterization of G8P[8] rotavirus emerged among

14

29. Zehender G, Ebranati E, Lai A, Santoro MM, Alteri C, Giuliani M, Palamara G, Perno 240

CF, Galli M, Lo Presti A, Ciccozzi M. 2010. Population dynamics of HIV-1 subtype B in a 241

cohort of men-having-sex-with-men in Rome, Italy. J Acquir Immune Defic Syndr 55:156-160. 242

30. Ciccozzi M, Babakir-Mina M, Lo Presti A, Marcuccilli F, Perno CF, Ciotti M. 2011. 243

Phylogenesis and Clinical Aspects of Pandemic 2009 Influenza A (H1N1) Virus Infection. The 244

open virology journal 5:22-26. 245

31. De Grazia S, Martella V, Colomba C, Cascio A, Arista S, Giammanco GM. 2009. Genetic 246

characterization of G3 rotaviruses detected in Italian children in the years 1993-2005. J Med 247

Virol 81:2089-2095. 248

32. Maes P, Matthijnssens J, Rahman M, Van Ranst M. 2009. RotaC: a web-based tool for the 249

complete genome classification of group A rotaviruses. BMC Microbiol 9:238. 250

33. Dhama K, Chauhan RS, Mahendran M, Malik SV. 2009. Rotavirus diarrhea in bovines and 251

other domestic animals. Vet Res Commun 33:1-23. 252

34. Monini M, Cappuccini F, Battista P, Falcone E, Lavazza A, Ruggeri FM. 2008. Molecular 253

characterization of bovine rotavirus strains circulating in northern Italy, 2003-2005. Vet 254

Microbiol 129:384-389. 255

35. Matsuno S, Hasegawa A, Mukoyama A, Inouye S. 1985. A candidate for a new serotype of 256

human rotavirus. J Virol 54:623-624. 257

36. Palombo EA, Clark R, Bishop RF. 2000. Characterisation of a "European-like" serotype G8 258

human rotavirus isolated in Australia. J Med Virol 60:56-62. 259

37. Armah GE, Pager CT, Asmah RH, Anto FR, Oduro AR, Binka F, Steele D. 2001. 260

Prevalence of unusual human rotavirus strains in Ghanaian children. J Med Virol 63:67-71. 261

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 15: 1 Full-genome characterization of G8P[8] rotavirus emerged among

15

38. Volotao EM, Soares CC, Maranhao AG, Rocha LN, Hoshino Y, Santos N. 2006. Rotavirus 262

surveillance in the city of Rio de Janeiro-Brazil during 2000-2004: detection of unusual strains 263

with G8P[4] or G10P[9] specificities. J Med Virol 78:263-272. 264

39. Rahman M, Matthijnssens J, Yang X, Delbeke T, Arijs I, Taniguchi K, Iturriza-Gomara 265

M, Iftekharuddin N, Azim T, Van Ranst M. 2007. Evolutionary history and global spread of 266

the emerging g12 human rotaviruses. J Virol 81:2382-2390. 267

40. Steyer A, Poljsak-Prijatelj M, Barlic-Maganja D, Marin J. 2008. Human, porcine and 268

bovine rotaviruses in Slovenia: evidence of interspecies transmission and genome reassortment. 269

The Journal of general virology 89:1690-1698. 270

41. Browning GF, Snodgrass DR, Nakagomi O, Kaga E, Sarasini A, Gerna G. 1992. Human 271

and bovine serotype G8 rotaviruses may be derived by reassortment. Arch Virol 125:121-128. 272

42. Ward R. 2009. Mechanisms of protection against rotavirus infection and disease. Pediatr Infect 273

Dis J 28:S57-59. 274

43. Offit PA, Shaw RD, Greenberg HB. 1986. Passive protection against rotavirus-induced 275

diarrhea by monoclonal antibodies to surface proteins vp3 and vp7. J Virol 58:700-703. 276

44. Dennehy PH, Cortese MM, Begue RE, Jaeger JL, Roberts NE, Zhang R, Rhodes P, 277

Gentsch J, Ward R, Bernstein DI, Vitek C, Bresee JS, Staat MA. 2006. A case-control study 278

to determine risk factors for hospitalization for rotavirus gastroenteritis in U.S. children. Pediatr 279

Infect Dis J 25:1123-1131. 280

45. Clemens J, Rao M, Ahmed F, Ward R, Huda S, Chakraborty J, Yunus M, Khan MR, Ali 281

M, Kay B, et al. 1993. Breast-feeding and the risk of life-threatening rotavirus diarrhea: 282

prevention or postponement? Pediatrics 92:680-685. 283

46. Kirkwood CD. 2010. Genetic and antigenic diversity of human rotaviruses: potential impact on 284

vaccination programs. J Infect Dis 202 Suppl:S43-48. 285

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 16: 1 Full-genome characterization of G8P[8] rotavirus emerged among

16

47. Esona MD, Geyer A, Banyai K, Page N, Aminu M, Armah GE, Hull J, Steele DA, Glass 286

RI, Gentsch JR. 2009. Novel human rotavirus genotype G5P[7] from child with diarrhea, 287

Cameroon. Emerg Infect Dis 15:83-86. 288

48. Ramachandran M, Gentsch JR, Parashar UD, Jin S, Woods PA, Holmes JL, Kirkwood 289

CD, Bishop RF, Greenberg HB, Urasawa S, Gerna G, Coulson BS, Taniguchi K, Bresee 290

JS, Glass RI. 1998. Detection and characterization of novel rotavirus strains in the United 291

States. J Clin Microbiol 36:3223-3229. 292

49. Santos N, Lima RC, Nozawa CM, Linhares RE, Gouvea V. 1999. Detection of porcine 293

rotavirus type G9 and of a mixture of types G1 and G5 associated with Wa-like VP4 specificity: 294

evidence for natural human-porcine genetic reassortment. J Clin Microbiol 37:2734-2736. 295

296

297

298

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 17: 1 Full-genome characterization of G8P[8] rotavirus emerged among

17

Legends to the Figures 299

300

FIG. 1. Bayesian phylogenetic analysis of Human Rotavirus VP7 and VP4 nucleotide sequences. The 301

trees were rooted by using the midpoint rooting method. Branch lengths were estimated with the best 302

fitting nucleotide substitution model according to a hierarchical likelihood ratio test. For the significant 303

statistical support, posterior probabilities > 0.90 are reported. Several clusters are replaced by a triangle 304

for simplicity, and its size is an indication of the number of sequences it represents. The number of the 305

strain for each place of isolation is indicated in parenthes. RV strain isolated in Croatia is in boldface 306

type. For each strain, the following data are given: species of origin/place of isolation/year/strain name. 307

The number of substitutions per site is indicated by the scale bar. 308

(A) Phylogenetic tree based on VP7 gene nucleotide sequence. (B) Detail of the phylogenetic 309

relationship between G8 genotype RV strains. 310

(C) Phylogenetic tree based on VP4 gene nucleotide sequence. (D) Detail of the phylogenetic 311

relationship between CR2006 and other RV strains. 312

313

FIG. 2. Bayesian phylogenetic analysis of Human Rotavirus NSP4 nucleotide sequences. The tree was 314

rooted by using the midpoint rooting method. Branch lengths were estimated with the best fitting 315

nucleotide substitution model according to a hierarchical likelihood ratio test. For the significant 316

statistical support, posterior probabilities > 0.90 are reported. Several clusters are replaced by a triangle 317

for simplicity, and its size is an indication of the number of sequences it represents. The number of the 318

strain for each place of isolation is indicated in parenthes. Phylogenetic clusters containing CR2006 are 319

boxed and the strain is in boldface type. For each strain, the following data are given: species of 320

origin/place of isolation/year/strain name. The number of substitutions per site is indicated by the scale 321

bar. 322

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 18: 1 Full-genome characterization of G8P[8] rotavirus emerged among

VP7

A.

G8

G1

G12G3

G10

G2

G4

0.2

Hu / Jap / 1982 / G3P9 / AU-1 La / Italy / 1996 / G3P14 / 30/96

Bo / India / 2007 / G8P14 (3)

Hu / Philip / G12P4 / L26

Hu / Africa / 2004 / G8P6 / SA4948JHB/04

Hu / Australia / 2002 / G8P14 (2)

Hu / G1P8 / KU

Bo / USA / 1996 / G8 / Cody I801

Hu / India / G9 / RMC321

Hu / China / 1996 / G2P4 / TB-ChenBo / G6 / UK - NCDV

Bo / Kor / G8 (28)

Bo / USA / 2003 / G8 (4)

Hu / Malawi / 2000 / G8P8 / MW4103

Hu / India / G8P1 / MP409

Bo / USA / 2007 / G8

Hu / Finland / HAL 1166Hu / Australia / G8 (3)

Hu / Congo / 2003 / G8P8 / DRC88

Po / Kor / 2006 (15)

Bo / Thailand / G10

Ov / Spain / G8P14 / OVR762

Hu / USA / 1976 / G2P4 / DS-1

Lama / Arg / 1998 / G8P14 (4)

Hu / Cameroon / 2000 / G8 / ARN (4)

Hu / Nigerian / 2003 / G8 (2)

Po / USA / 1976 / G4P6 / Gottfried

Po / G5P7 / OSU

Bo / Jap / G8 (5)

Hu / UK / 1985 / G8 (2)

Hu / USA / 1974 / G1P8 / Wa

Si / Kenia / 1999 / G8P6 / Si/KY1646

Hu / Egypt / G8P14 / EGY1850

Bo / Nigeria / 2003 / G8 / NGRBg8

Hu / Tunisia / 2002-04 / G8 (3)

Hu / Switz / G8 / U25.96

Hu / Croatia / 2006 / G8P8 / CR2006

Hu / Congo / 2003 / G8P6 / DRC86

Hu / Malawi / 2000-01 / G8 (7)

Si / Sa11 / G3Hu / Bangl / 2000-03 / G12P6 (2)

Hu / Africa / G8 / QEH14262

Bo / Thailand / G8 / A5Hu / Jap / 1994-05 / G8 (2)

Hu / Egypt / G8P14 / EGY2295

Hu / 69M / G8P10

Hu / India / G10P11 / I321

Hu / Ger / 2009 / G8P4 / GER1H-09

Hu / Ethiopia / 2004 / G8 / ARN6810

Si / USA / PTRV

Hu / Kenia / 2002 / G8 (3)

Hu / Brazil / G8P4 / R291

Hu / UK / 678

Hu / Slo / 2006 / G8P8 / SI-885/06

1

0,97

1

1

1

1

1

1

0,951

1

0,991

1

1

0,94

0,91

1

0,91

1

1

1

1

Hu / Cameroon / 2000 / G8 / ARN6780

Hu / Ethiopia / 2004 / G8 / ARN6810

Hu / Tunisia / 2000 / G8 / ARN6862

Hu / Malawi / 2000 / G8P8 / MW4103

Hu / Tunisia / 2004 / G8 / ARN6736

Hu / Ger / 2009 / G8P4 / GER1H-09

Hu / Cameroon / 2000 / G8 / ARN6809

Hu / Cameroon / 2000 / G8 / ARN6782

Hu / Croatia / 2006 / G8P8 / CR2006

Hu / Tunisia / 2002 / G8 / ARN6854

Hu / Slo / 2006 / G8P8 / SI-885/06

Hu / Cameroon / 2000 / G8 / ARN67871

0,94

B.

FIG. 1. Bayesian phylogenetic trees of Human Rotavirus VP7 and VP4 nucleotide sequences. The trees were rooted by using the midpoint rooting method. Branch lengths were estimated with the best fitting nucleotide substitution model according to a hierarchical likelihood ratio test. For the significant statistical support, posterior probabilities > 0.90 are reported. Several clusters are replaced by a triangle for simplicity, and its size is an indication of the number of sequences it represents. The number of the strain for each place of isolation is indicated in parenthes. RV strain isolated in Croatia is in boldface type. For each strain, the following data are given: species of origin/place of isolation/year/strain name. The number of substitutions per site is indicated by the scale bar. (A) Phylogenetic tree based on VP7 gene nucleotide sequence. (B) Detail of the phylogenetic relationship between G8 genotype RV strains. (C) Phylogenetic tree based on VP4 gene nucleotide sequence.

0.3

Hu / Ger / 2009 / G8P4 / GER1H-09

Hu / India / G9 / RMC321

Hu / Bel / 2003 / G12P8 / B4633-03

Hu / China / 2007 / G3P8 (4)

Hu / Croatia / 2006 / G8P8 / CR2006Hu / Slo / 2006 / G8P8 / SI-885/06

Hu / Bangl / / G12P8 / Dhaka25-02

Hu / India / 2006 / G1P8 (3)

Hu / Bel / 2009 / G1P8 (8)

Hu / USA / 1974 / G1P8 / Wa

Hu / Jap / P8 (2)

Hu / USA / 2010 / G9P11 / 116E/AG

Hu / Italy / 2005 / G9P8 / BIA2

Hu / Bel / 2008 / G1P8 (6)

Hu / Italy / 2006 / G9P8 (2)

Hu / Russia / 2004-09 / P8 (16)

Hu / Russia / 2004-07 / P8 (2)

Hu / Russia / 2008 / G4P8 (3)

Hu / Bangl / 2003 / G12P6 (2)

Hu / China / 2007 / G3P8 (4)

Hu / Congo / 2003 / G8P8 / DRC88

Hu / Africa / 2000-04 / G8 / ARN (5)

Hu / Russia / 2009 / G4P8 (4)

Hu / Thailand / 2009 / P8 (7)

La / Italy / 1996 / G3P14 / 30/96

Hu / China / 2007 / P8 (6)

Hu / USA / 2007 / G1P8

Hu / Bangl / 2005 / P8 (9)

Hu / China / 1996 / G2P4 / TB-ChenHu / USA / 1976 / G2P4 / DS-1

Hu / India / 2006 / P8 / mani-133/06

Hu / Congo / 2003 / G8P6 / DRC86

Hu / GBR / 1975 / G4P6 / ST3

Hu / Bel / 2006 / G1P8 (7)Hu / Russia / 2006 / G9P8 / RUS-Nov06-1289

Po / USA / 1976 / G4P6 / Gottfried

Hu / USA / G1P8 / US6668

Hu / Ger / 2008 / G3G12P8 / GER126-08

Hu / Malawi / 1997 / P8 (4)

Hu / Italy / 2005 / G9P8 (2)

0,93

0,99

1

0,95

1

1

0,93

1

1

0,98

1

1

1

0,99

VP4

C.

P[8]

P[4]

P[11]

P[6]

P[14]

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from

Page 19: 1 Full-genome characterization of G8P[8] rotavirus emerged among

0.3

Hu / Jap / G3 (2)

Hu / Tunisia / 2000-04 / G8 / ARN (3)

Hu / G3P2 / RV3

Hu / Jap / G3 / ITO

Hu / Bangl / 2000-03 / G12P6 (4)

Hu / Bel / 2009 / G1P8 (2)

Hu / G1P8 / KU

La / Italy / 1996 / G3P14 / 30/96

Hu / GBR / 1975 / G4 / ST3

Hu / Italy / 2006 / G3P6 / 128/07-34

Hu / South Africa / 1986 / G4P6 (4)

Hu / Australia / G4 (2)

Hu / Italy / 2005 / G3P8 / PA116

Hu / Russia / G1P8 / N.N.12813

Hu / USA / 2008 / G1P8

Hu / China / G3 (4)

Hu / Ger / 2008 / G3G12P8 / GER126-08

Hu / Russia / 2008-09 / G9P8 (2)

Hu / China / G1 (2)

Hu / Russia / 2004-09 / P8 (6)

Hu / USA / 2005-08 / G1P8 (3)

Hu / USA / G1P8 (3)

Hu / Slo / 2006 / G8P8 / SI-885/06

Hu / Ethiopia / 2004 / G8 / ARN6810

Hu / Bel / 2006 / G1P8 (2)

Hu / China / 1996 / G2P4 / TB-ChenHu / Congo / 2003 / G8 / DRC (2)

Hu / Russia / G1P8 / N.N.131

Hu / Russia / 2007 / G4P8 / Omsk07-78

Hu / USA / 1976 / G2P4 / DS-1

Hu / Italy / 2005 / G3P8 / PA106

Hu / Italy / 2005 / G3P8 / PA110

Hu / Jap / G3 (2)

Hu / USA / G1P8 (7)

Hu / China / 1994 / G4 (2)Hu / Jap / 1990 / G4 (2)

Hu / Jap / G3 (4)

Hu / Italy / 2005 / G3P8 / PA49

Hu / Italy / 2005 / G3P8 / PA57

Hu / Ger / 2009 / G8P4 / GER1H-09

Hu / India / G9/ RMC321

Hu / Bel / 2004-08 / G1P8 (4)

Po / OSU / G5P7

Hu / Bel / 2006-08 / G1P8 (3)

Hu / Russia / 2009 / G1P8 / D83Hu / India / 2006 / G1P8 / AM06-I

Hu / Russia / G4P8 / N.N.13518

Hu / Bel / 2007-09 / G1P8 (4)

Hu / Cameroon / 2000 / G8 / ARN6780

Hu / Jap /1982 / G3P9 / AU-1

Hu / Bel / 2005-06 / G1P8 (2)

Hu / Croatia / 2006 / G8P8 / CR2006

Hu / Italy / 2005 / G3P8 / PA105

Hu / Bel / 2007 / G1P8

Hu / India / G10 / I321Bo / Kor / CBN-1

Hu / USA / G1P8 / Wa (6)

Hu / Jap / G3-G4 (2)

1

1

11

1

1

1

1

1

0,96

1

1

1

1 1

1

1

1

1

1

1

1NSP4

E1

E2

E3E5

FIG. 2. Bayesian phylogenetic tree of Human Rotavirus NSP4 nucleotide sequence. The tree was rooted by using the midpoint rooting method. Branch lengths were estimated with the best fitting nucleotide substitution model according to a hierarchical likelihood ratio test. For the significant statistical support, posterior probabilities > 0.90 are reported. Several clusters are replaced by a triangle for simplicity, and its size is an indication of the number of sequences it represents. The number of the strain for each place of isolation is indicated in parenthes. Phylogenetic clusters containing CR2006 are boxed and the strain is in boldface type. For each strain, the following data are given: species of origin/place of isolation/year/strain name. The number of substitutions per site is indicated by the scale bar.

on April 3, 2018 by guest

http://jcm.asm

.org/D

ownloaded from