PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2017; 12(9): e0184067.
Published online 2017 September 13. doi:  10.1371/journal.pone.0184067
PMCID: PMC5597190

Identification of vaccine-derived rotavirus strains in children with acute gastroenteritis in Japan, 2012-2015

Mei Kaneko, Conceptualization, Data curation, Formal analysis, Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing,1 Sayaka Takanashi, Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing,1,* Aksara Thongprachum, Data curation, Investigation,2 Nozomu Hanaoka, Formal analysis, Methodology, Writing – review & editing,3 Tsuguto Fujimoto, Formal analysis, Methodology, Writing – review & editing,3 Koo Nagasawa, Data curation, Formal analysis, Investigation, Methodology, Writing – review & editing,3 Hirokazu Kimura, Formal analysis, Methodology, Resources, Writing – review & editing,3 Shoko Okitsu, Formal analysis, Methodology, Supervision, Writing – review & editing,1,2 Masashi Mizuguchi, Supervision, Writing – review & editing,1 and Hiroshi Ushijima, Conceptualization, Formal analysis, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing1,2
Dong-Yan Jin, Editor

Abstract

Two live attenuated oral rotavirus vaccines, Rotarix and RotaTeq, have been introduced as voluntary vaccination in Japan since 2011 and 2012, respectively. Effectiveness of the vaccines has been confirmed, whereas concerns such as shedding of the vaccine strains and gastroenteritis cases caused by vaccine strains are not well assessed. We aimed to identify the vaccine strains in children with acute gastroenteritis (AGE) to investigate the prevalence of AGE caused by vaccination or horizontal transmission of vaccine strains. A total of 1,824 stool samples were collected from children with AGE at six outpatient clinics in 2012–2015. Among all, 372 group A rotavirus (RVA) positive samples were screened for vaccine components by real-time RT-PCR which were designed to differentiate vaccine strains from rotavirus wild-type strains with high specificity. For samples possessing both vaccine and wild-type strains, analyses by next-generation sequencing (NGS) were conducted to characterize viruses existed in the intestine. As a result, Rotarix-derived strains were identified in 6 of 372 (1.6%) RVA positive samples whereas no RotaTeq strain was detected. Among six samples, four possessed Rotarix-derived strains while two possessed both Rotarix-derived strains and wild-type strains. In addition, other pathogens such as norovirus, enterovirus and E.coli were detected in four samples. The contribution of these vaccine strains to each patient’s symptoms was unclear as all of the cases were vaccinated 2–14 days before sample collection. Proportion of average coverage for each segmented gene by NGS strongly suggested the concurrent infection of the vaccine-derived strain and the wild-type strain rather than reassortment of these two strains in one sample. This is the first study to report the prevalence of vaccine-derived strains in patients with RVA AGE in Japan as 1.6% without evidence of horizontal transmission. The results emphasized the importance of continuous monitoring on vaccine strains and their clinical impacts on children.

Introduction

Rotavirus infection is a leading cause of acute gastroenteritis (AGE) among children under five years of age causing approximately 215,000 deaths in a year worldwide [1]. Rotaviruses are nonenveloped triple layered viruses which contain 11 segments of double-stranded RNA coding for six structural proteins (VP1-VP4, VP6 and VP7) and six nonstructural proteins (NSP1-NSP5/6) (2). Group A rotaviruses (RVA) are a leading cause of AGE in children and are further classified based on genotypes of two outer capsid proteins, VP7 (G-type) and VP4 (P-type) [2]. In addition to G/P-genotyping, a complete genome classification system was developed by rotavirus classification working group (RCWG) [3]. The VP7-VP4-VP6-VP1-VP2-VP3-NSP1-NSP2-NSP3-NSP4-NSP5/6 genes of each strain were described as Gx-Px-Ix-Rx-Cx-Mx-Ax-Nx-Tx-Ex-Hx (x = Arabic numbers starting from one, representing each genotype) to represent a genotype constellation. Since RNA segments from different rotaviruses within the same group would reassort with high frequency during dual infections of cells, a comprehensive classification on 11 genome segments allows us to identify reassortments to confirm interspecies transmission [2,4,5].

Two live attenuated oral rotavirus vaccines, Rotarix (GlaxoSmithKine, Rixensart, Belgium) and RotaTeq (Merck & Co., Whitehouse Station, NJ, USA) were first licensed in 2004 in Mexico and the Dominican Republic, and in 2006 in the United States (US), respectively. Rotarix is a two-dose monovalent human rotavirus vaccine while RotaTeq is a three-dose pentavalent bovine-human reassortant vaccine [6,7]. Eighty six countries have introduced them into their national immunization programs since WHO recommended routine immunization of rotavirus vaccines to all infants in 2009 [8,9]. Both vaccines were found to be highly effective against severe rotavirus gastroenteritis and rotavirus-associated hospitalization in early vaccine-introducing countries such as the US, Mexico, Belgium and Australia [1014]. Following vaccine introduction in 2006 (RotaTeq) and 2008 (Rotarix) in the US, the effectiveness for ≥1 dose of any rotavirus vaccine was 80% (95% CI, 74%–85%) against rotavirus hospitalizations and emergency department visits in 2009–2011 [10].

Despite the benefits of rotavirus vaccines, several concerns need to be mentioned. The first is shedding of vaccine strains. After the first dose of administration, the peak time of shedding occurred between day 4 and day 7, with positive detection rate of 80–90% by real-time RT-PCR [15]. For immunocompetent infants, the shedding was detected as early as one day and as late as 28 days after the administration. The duration of shedding is prolonged in immunocompromised infants, for instance, for more than six months in infants with severe combined immune deficiency (SCID) [16,17]. The second is AGE caused by vaccine strains or vaccine-derived reassortants. Some cases resulted from administration of the vaccines while others from horizontal transmissions of vaccine strains. According to a previous study in the US, vaccine or vaccine-reassortant rotavirus strains were detected in stool samples from five of 106 (4.7%) immunocompetent children who required treatment for rotavirus gastroenteritis [18]. To provide the evidence of horizontal transmission of the vaccine strains, a study documented the occurrence of RotaTeq vaccine-derived rotavirus transmission from a vaccinated infant to his older unvaccinated sibling, resulting in symptomatic AGE which required emergency department care [19]. In addition, several studies have reported severe and chronic gastroenteritis caused by vaccine strains in children with SCID [17,20].

Japan is still in the initial stage of using rotavirus vaccines to control rotavirus infections among children. Rotarix and RotaTeq were introduced as voluntary vaccination in November 2011 and July 2012, respectively. The effectiveness of the vaccines has been confirmed since two vaccines were introduced with the average vaccine coverage being 47.6% between July 2012 and November 2014 [2123]. However, the concerns as mentioned above remain to be clarified. Risk-benefit analyses on rotavirus vaccines are essential in the introduction of these two vaccines into the national immunization program. Therefore, the objective of this study is to investigate the prevalence of rotavirus vaccine derived strain in Japanese children with AGE to assess the risks of rotavirus vaccines by elucidating horizontal transmission of the vaccine strains.

Materials and methods

Study samples

A total of 1,824 stool samples from patients with AGE at six outpatient pediatric clinics in Hokkaido, Tokyo, Shizuoka, Kyoto, Osaka and Saga were collected from July 2012 to June 2015. The demographic characteristics and clinical symptoms were collected from medical records based on the factors described in a previous study [18]. Age of the patients ranged from neonate to 15 years (0–178 months; median age of 17 months). All stool samples were initially screened by a method of multiplex RT-PCR which was designed for detection of 11 diarrhea causing viruses including rotavirus (group A, B, and C), norovirus (GI and GII), adenovirus, astrovirus, sapovirus, aichivirus, parechovirus and enterovirus [24]. Additionally, stool samples were cultured and tested for bacteria when patient exhibited suspicious symptoms such as bloody stool.

Viral RNA extraction

Viral RNA was extracted from 10% (w/v) stool suspension with carrier RNA using QIAamp Viral RNA mini kit (QIAGEN, Hilden, Germany) according to the manufacturer`s instruction. For next-generation sequencing (NGS), RNA was extracted using the same method but without adding carrier RNA.

Real-time RT-PCR for detection of vaccine components

Rotarix NSP2 assay and RotaTeq VP6 assay are quantitative real-time RT-PCR assays for detection of Rotarix and RotaTeq vaccine components in stool samples [25]. The assays were designed for vaccine specific targets in genomes of Rotarix NSP2 gene (281 bp) and RotaTeq VP6 gene (100 bp), which were the appropriate targets to differentiate the vaccine strains from wild-type RVA strains. Rotavirus NSP3 assay was established for detection of wide range of clinically important rotaviruses [26] and was used to titrate wild-type strains in this study. All the assays were conducted using Applied Biosystems StepOnePlus Real-Time PCR System (Applied Biosystems, Foster city, CA, USA) in accordance with the procedures described in previous studies [25,26]. Primers and probes used in the assays are described in S1 Table.

Quantitative real-time RT-PCR for other viral pathogens

Quantitative real-time RT-PCR for norovirus and enterovirus were carried out as previously described [27,28].

Conventional RT-PCR and sanger sequencing

Reverse transcription was performed with denatured RNA template using SuperScript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) and random primer (Takara, Shiga, Japan) according to the manufacturer's instructions. For those samples that had difficulty in amplifying target genes in the downstream PCR, gene specific primers were used instead of random primer to increase the specificity and sensitivity in PCR (S2 Table).

PCR was performed to amplify whole or part of four rotavirus segmented genes: VP7, VP4, VP6 and NSP2. The primers used for each gene are shown in S2 Table. The reaction was carried out using GoTaq DNA Polymerase (Promega, Madison, WI, USA) in accordance with manufacturer's instructions and PCR conditions in the previous studies shown in S2 Table. For those samples that had difficulties in amplification, PCR was then performed using PrimeSTAR GXL DNA polymerase (Takara, Shiga, Japan) which had higher efficiency [29].

Sequencing was conducted for PCR amplicons using ABI Prism BigDye termination cycle sequencing reaction kit (Applied Biosystems, Foster, CA, USA) according to the manufacturer`s instruction. The sequence data were collected by an ABI Prism 310 Genetic Analyzer (Applied Biosystems, Foster city, CA, USA). Nucleotide alignments were edited and analyzed using MEGA 6 [30], BioEdit [31] and Basic Local Alignment Search Tool (BLAST) [32]. The obtained sequences were aligned and compared with the corresponding genes of Rotarix original strain (accession No.JX943614, JX943612, JX943613 and JX943605 for VP7, VP4, VP6 and NSP2 gene, respectively) [25].

Analysis by NGS

RNA sequencing library was prepared using NEBNext Ultra RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA). Nucleotide sequencing was carried out on Miseq II (Illumina, San Diego, CA, USA) using MiSeq Reagent Kit v3 to generate 300-cycle paired-end reads. To construct contigs from the obtained short reads, de novo assembling were conducted using CLC genomics workbench (QIAGEN, Hilden, Germany) and VirusTap [33]. Contigs derived from either rotavirus wild-type strains or Rotarix strain were identified by BLAST search. Genotypes of each contigs were determined by an automated genotyping tool RotaC 2.0 (http://Rotac.regatools.be/) according to the guideline proposed by RCWG.

GenBank accession numbers

The representative gene sequences of RVA strains obtained in this study were deposited into GenBank under the accession numbers KY616885-KY616907 (S3 Table).

Ethical clearance

Written informed consents were obtained from the guardians on behalf of the minors/children enrolled in our study. The study was approved by the ethical committees of The University of Tokyo and Nihon University under the registration number of 1139 and 25-13-0, respectively.

Results

Demographic characteristics and clinical symptoms of positive cases in Rotarix NSP2 assay

In total, 372 stool samples were positive for RVA and analyzed further in this study. Among them, partial NSP2 gene of Rotarix strain was detected in seven samples (1.9%), whereas no products of RotaTeq VP6 gene strain were detected. The demographic characteristics and clinical symptoms of seven positive cases in Rotarix NSP2 assay are summarized in Table 1. Except for patient No.3, all patients had been administered either 1st dose or 2nd dose of Rotarix and the samples were collected one to 14 days after the vaccinations. All cases exhibited mild and various combinations of AGE symptoms.

Table 1
Demographic characteristics and clinical symptoms of seven positive cases in Rotarix NSP2 assay.

Viral titration and sanger sequencing of positive samples in Rotarix NSP2 assay

The results of viral titration and sequencing of these seven positive samples are shown in Table 2. Sequencing of four genes (VP7, VP4, VP6, and NSP2) by conventional RT-PCR identified Rotarix-derived strains in four samples (sample No.2, 4, 5 and 7), and both Rotarix-derived strains and rotavirus wild-type strains in two (sample No.1 and 6).

Table 2
Results of viral titration and sanger sequencing of four genes of seven positive samples in Rotarix NSP2 assay.

For sample No.1, VP7 (1051 bp), VP4 (803 bp), VP6 (1237 bp) and NSP2 (1030 bp) genes were highly similar (99.9–100% identity) with those of DS-1-like G1P[8] wild-type strain (HC12016) detected in Osaka in 2012 [34]. Along with those wild-type genes, sequences of VP6 (341 bp) and NSP2 (281 bp) which were 100% identical with those of Rotarix original strain were also obtained by conventional RT-PCR. The viral titer of rotavirus wild-type strain was 13.08 log 10 copies/g of stool.

For sample No.6, VP7 (1014 bp), VP4 (811 bp), VP6 (339 bp) and NSP2 (953 bp) genes were highly similar (95.9–99.9% identity) with Wa-like G1P[8] wild-type strain (CK00110) detected in 2011 in Victoria, Australia [35]. Along with those wild-type genes, sequence of NSP2 (281 bp) which was 100% identical with that of Rotarix original strain was also obtained. The viral titer of rotavirus wild-type stain could not be calculated probably due to its low viral titer.

For sample No.3, VP7 (1018 bp), VP4 (810 bp), VP6 (1305 bp) and NSP2 (950 bp) genes were highly similar (99.5–99.9% identity) with those of rotavirus wild-type strains and they had genome constellation of G9-P[8]-I1-N1. Based on the results of sequencing on four genes, the results of Rotarix NSP2 assay was considered false positive.

For two samples (No.1 and No.6) that possessed both Rotarix-derived strains and rotavirus wild-type strains, analyses by NGS were carried out to characterize the viruses that concurrently existed in the intestine. Four out of seven patients (No.1, No.5, No.6 and No.7) had concurrent infection with other pathogens such as norovirus, E.coli and enterovirus based on the results of multiplex RT-PCR and stool culture. The viral titers of these viruses are also described in Table 2.

Analysis by NGS

In sample No.1, 8,786 contigs were constructed by CLC genomics workbench and 24 contigs by VirusTap. All the 11 genes of both rotavirus wild-type strain and Rotarix strain were finally identified by BLAST search. Fig 1 shows the maps of short reads to cover each contig of VP6 gene, as a representative example. The number of reads to construct one contig is apparently much larger in wild-type strain than in Rotarix strain. The number of average coverage for each contig is an indicator to show the average times that each nucleotide is sequenced [36]. Accordingly, we calculated the percentage of wild-type strain and Rotarix strain for each gene based on the number of average coverage. As a result, rotavirus wild-type strain existed with higher percentage (92.9–99.3%) compared to Rotarix strain (0.7–7.1%) in all the 11 genes as shown in Fig 2.

Fig 1
Maps of short reads to cover each contig of VP6 gene in sample No.1.
Fig 2
Percentage of wild-type strain and Rotarix strain in sample No.1.

Nearly whole nucleotide sequences of all the 11 genes of wild-type strain were obtained in sample No.1. According to the results of genotyping and BLAST search, the genome constellation of the wild-type strain was G1-P[8]-I2-R2-C2-M2-A2-N2-T2-E2-H2, which represented a DS-1-like constellation. In addition, all the genes were highly similar (99.8–100% identity) with those of DS-1-like G1P[8] strain detected in Osaka in 2012 (HC12016, accession nos. AB848004-AB848014) [34].

In sample No.6, 2,500 contigs were constructed by CLC genomics workbench and 10 contigs by VirusTap. Based on the result of BLAST search, contigs derived from rotaviruses could not be identified. Coxsackievirus B4 (CB4), one of the members of genus Enterovirus within the family Picornaviridae, was detected among contigs constructed by both CLC genomic workbench and Virus Tap.

Discussion

This study was conducted right after the two rotavirus vaccines had been introduced as voluntary vaccination in Japan. With the vaccine coverage across the country being approximately 50% during the study period, we detected six out of 372 (1.6%) RVA cases of AGE which possessed rotavirus vaccine-derived strains. This is relatively low prevalence compared to the one in a previous study in the US where they reported 4.7% of such cases with the vaccine coverage being around 60% [18,37].

In patient No.2 and No.4, the obtained sequences of VP7, VP4, VP6 and NSP2 gene were 100% identical with those of Rotarix original strain. No other viral or bacterial pathogen was detected from these patients. They had developed mild symptoms of AGE within six days after being vaccinated with 1st dose of Rotarix. As neither of them had siblings nor had attended a daycare, we considered that Rotarix-derived strains in these two patients were vaccines shed after the vaccination, rather than horizontal transmission event. Since we were not able to rule out non-infectious diarrhea, such as intolerances, food allergies, and secondary to the use of unrecorded medicines, the contribution of the vaccine strains to these AGE patients’ symptoms remained unclear.

For patient No.5, we suspected that the main cause of AGE was E.coli O6, based on the clinical symptom of bloody stool, a finding suggestive of bacterial rather than viral AGE. In a similar manner, the AGE case of patient No.7 probably had been caused by norovirus GII.3 based on the results of real-time RT-PCR which showed viral titer of norovirus was high while that of Rotarix-derived strain was low [38].

In patient No.1 and No.6 from whom both the Rotarix-derived strains and rotavirus wild-type strains were detected, we firstly performed electropherotype profiling by PAGE to confirm the concurrent infection of the two strains and to capture any signs of rearrangement or reassortment in vain (S1 Text). Next, we carried out NGS analyses aiming to characterize the viruses that concomitantly existed in the intestine. In sample No.1, rotavirus wild-type strain accounted for more than 92.9% in all 11 genes. We interpreted NGS data as strongly suggestive of concurrent infection of vaccine-derived strain and the wild-type strains without clear evidence of reassortant strains. Attenuated rotaviruses (such as Rotarix strain) are shed at much lower levels than wild-type strains, and so would the re-assorted strain. Therefore, the higher prevalence of the wild-type rotavirus genomic segments versus vaccine or potential reassortant strain should be expected. Only virus isolation via liquid limited dilutions followed by complete genomic sequencing could distinguish concurrent infection versus reassortment event. In our case, however, this strategy was hindered due to low amount of remaining sample volume. Though norovirus was detected as well, the viral titer was quite low while those of rotavirus wild-type strain (DS-1-like G1P[8]) and Rotarix strain were high. Based on these analyses, we inferred that rotavirus wild-type strain was the probable pathogen of AGE in sample No.1. The infant might have got infected with the wild-type strain soon after he was vaccinated with 1st dose of Rotarix and right before the immune response against rotavirus was established. The Rotarix strain detected in this sample was likely the shed vaccine strain rather than the strain from any horizontal transmission event.

In sample No.6, contigs of rotavirus wild-type strain and Rotarix-derived strain were not successfully constructed. A previous study has demonstrated that the sensitivity and quantitative performance of the NGS-based approach for detection of DNA viral sequences was equivalent to that of real-time PCR in cell-free clinical samples with a sequence depth of 5,000,000 reads [39]. Although Rotarix strain in sample No.6 was quantified as 5.89 log 10 copies/g of stool by real-time RT-PCR, instability of RNA virus and impurities contained in clinical stool sample might have affected the analytical sensitivity of NGS. On the other hand, the result of viral titration by real-time RT-PCR for enterovirus and analyses by NGS both demonstrated the predominance of CB4. Therefore, we inferred that CB4 was the probable pathogen of AGE in this sample.

Of note, among all six AGE cases which possessed Rotarix-derived strains, four (sample No.1, 5, 6 and 7) were suspected to be caused by other pathogens. Most likely, the infants were infected with other pathogens during the shedding period of Rotarix strain. This result exhibited a complex relationship between vaccine strains and clinical illnesses. To properly assess the potential risks of the vaccines, clinical illnesses associated with vaccine or vaccine-derived strains should be evaluated carefully by considering vaccination history and testing other viral or bacterial pathogens.

One false-positive sample (sample No.3) was detected in real-time RT-PCR. One of the reasons to explain this error would be the similarity of nucleotide sequences in both primer and probe binding sites. In order to improve the specificity of detecting vaccine strains, secondary screening assay of real-time RT-PCR, such as Rotarix VP4 assay, or genome sequencing is needed [25].

Although the generalizability of the findings in the current study is limited due to the confined surveillance sites, to the best of our knowledge, this is the first study to assess the prevalence of vaccine-derived strains in patients with AGE in Japan. RotaTeq strains were not detected in this study partly because the vaccine coverage of RotaTeq was lower than that of Rotarix during the study period. However, the AGE cases caused by direct vaccination or horizontal transmission of RotaTeq strains or its reassortants have been reported in other countries [18,19,40]. As the coverage of rotavirus vaccines has increased in Japan, investigation on both Rotarix and RotaTeq strains among AGE patients should be continued to monitor any horizontal transmission case or reassortment event associated with vaccine strains. By monitoring the prevalence of vaccine strains or vaccine-derived strains, we could assess the potential for sustained circulation of these strains.

Conclusion

The prevalence of vaccine-derived strains in patients with RVA AGE was 1.6% in Japan in 2012-2015. Neither horizontal transmissions nor reassortment event of vaccine strains were observed. The results emphasized the importance of continuous monitoring on vaccine strains and their clinical impacts on children.

Supporting information

S1 Table

Primers and probes used in real-time RT-PCR assays.

(PDF)

S2 Table

Primers used in conventional RT-PCR and sanger sequencing.

(PDF)

S3 Table

GenBank accession numbers for representative gene sequences of seven positive samples in Rotarix NSP2 assay.

(PDF)

S1 Text

Supplementary material for PAGE.

(PDF)

Acknowledgments

We sincerely thank Hideaki Kikuta, Atsuko Yamamoto, Masaaki Kobayashi, Kumiko Sugita, Shuichi Nishimura and Tsuneyoshi Baba for collection of samples and epidemiological data. This work was supported by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research (Grant number: 24390266, 16H05360) and Public Foundation of the Vaccination Research Center (2015–34).

Funding Statement

This work was supported by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research (Grant number: 24390266, 16H05360) and Public Foundation of the Vaccination Research Center (2015-34). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files.

References

1. Tate JE, Burton AH, Boschi-Pinto C, Parashar UD, Agocs M, Serhan F, et al. Global, Regional, and National Estimates of Rotavirus Mortality in Children <5 Years of Age, 2000–2013. Clin Infect Dis. 2016;62(Suppl 2):S96–105. [PubMed]
2. Estes MK, Greenberg HB. Rotaviruses In: Kinpe DM, Howley PM, editors. Fields virology sixth edition Philadelphia: Lippincott; 2013. 1347–401.
3. Matthijnssens J, Ciarlet M, Heiman E, Arijs I, Delbeke T, McDonald SM, et al. Full genome-based classification of rotaviruses reveals a common origin between human Wa-Like and porcine rotavirus strains and human DS-1-like and bovine rotavirus strains. J Virol. 2008;82(7):3204–19. doi: 10.1128/JVI.02257-07 [PMC free article] [PubMed]
4. Kirkwood CD. Genetic and Antigenic Diversity of Human Rotaviruses: Potential Impact on Vaccination Programs. J Infect Dis. 2010;202(S1):S43–8. [PubMed]
5. Matthijnssens J, Van Ranst M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr Opin Virol. 2012;2(4):426–33. doi: 10.1016/j.coviro.2012.04.007 [PubMed]
6. Ward RL, Bernstein DI. Rotarix: A Rotavirus Vaccine for the World. Clin Infect Dis. 2009;48(2):222–8. doi: 10.1086/595702 [PubMed]
7. Matthijnssens J, Joelsson DB, Warakomski DJ, Zhou T, Mathis PK, van Maanen M-H, et al. Molecular and biological characterization of the 5 human-bovine rotavirus (WC3)-based reassortant strains of the pentavalent rotavirus vaccine, RotaTeq. Virology. 2010;403(2):111–27. doi: 10.1016/j.virol.2010.04.004 [PubMed]
8. World Health Organization. Rotavirus Vaccination. Meeting of the immunization Strategic Advisory Group of Experts, April 2009. Wkly Epidemiol Rec. 2009;84(23):220–36. [PubMed]
9. Loharikar A, Dumolard L, Chu S, Hyde T, Goodman T, Mantel C. Status of New Vaccine Introduction—Worldwide, September 2016. Morb Mortal Wkly Reports. 2016;65(41):1136–40. [PubMed]
10. Payne DC, Boom JA, Staat MA, Edwards KM, Szilagyi PG, Klein EJ, et al. Effectiveness of pentavalent and monovalent rotavirus vaccines in concurrent use among US children <5 years of age, 2009–2011. Clin Infect Dis. 2013;57(1):13–20. doi: 10.1093/cid/cit164 [PMC free article] [PubMed]
11. Leshem E, Moritz RE, Curns AT, Zhou F, Tate JE, Lopman B a., et al. Rotavirus vaccines and health care utilization for diarrhea in the United States (2007–2011). Pediatrics. 2014;134(1):15–23. doi: 10.1542/peds.2013-3849 [PubMed]
12. Velázquez RF, Linhares AC, Muñoz S, Seron P, Lorca P, DeAntonio R, et al. Efficacy, safety and effectiveness of licensed rotavirus vaccines: a systematic review and meta-analysis for Latin America and the Caribbean. Pediatrics. 2017;17(1):14 doi: 10.1186/s12887-016-0771-y [PMC free article] [PubMed]
13. Sabbe M, Berger N, Blommaert A, Ogunjimi B, Grammens T, Callens M, et al. Sustained low rotavirus activity and hospitalisation rates in the post-vaccination era in Belgium, 2007 to 2014. Eurosurveillance. 2016;21(27). [PubMed]
14. Cortese MM, Immergluck LC, Held M, Jain S, Chan T, Grizas AP, et al. Cost effectiveness of both (Monovalent and Pentavalent) rotavirus vaccines. Pediatrics. 2013;7(1):13–20. [PubMed]
15. Hsieh YC, Wu FT, Hsiung C a., Wu HS, Chang KY, Huang YC. Comparison of virus shedding after lived attenuated and pentavalent reassortant rotavirus vaccine. Vaccine. 2014;32(10):1199–204. doi: 10.1016/j.vaccine.2013.08.041 [PubMed]
16. Werther RL, Crawford NW, Boniface K, Kirkwood CD, Smart JM. Rotavirus vaccine induced diarrhea in a child with severe combined immune deficiency. J Allergy Clin Immunol. 2009;124(3):600. [PubMed]
17. Kaplon J, Cros G, Ambert-Balay K, Leruez-Ville M, Chomton M, Fremy C, et al. Rotavirus vaccine virus shedding, viremia and clearance in infants with severe combined immune deficiency. Pediatr Infect Dis J. 2015;34(3):326–8. doi: 10.1097/INF.0000000000000560 [PubMed]
18. Boom JA, Sahni LC, Payne DC, Gautam R, Lyde F, Mijatovic-Rustempasic S, et al. Symptomatic infection and detection of vaccine and vaccine-reassortant rotavirus strains in 5 children: a case series. J Infect Dis. 2012;206(8):1275–9. doi: 10.1093/infdis/jis490 [PMC free article] [PubMed]
19. Payne DC, Edwards KM, Bowen MD, Keckley E, Peters J, Esona MD, et al. Sibling transmission of vaccine-derived rotavirus (RotaTeq) associated with rotavirus gastroenteritis. Pediatrics. 2010;125(2):e438–41. doi: 10.1542/peds.2009-1901 [PubMed]
20. Klinkenberg D, Blohm M, Hoehne M, Mas Marques A, Malecki M, Schildgen V, et al. Risk of Rotavirus Vaccination for Children with SCID. Pediatr Infect Dis J. 2015;34(1):114–5. [PubMed]
21. Committee for Immunization and vaccines. Evaluations and analyses on rotavirus vaccines (Report 1–3). Ministry of Health, Labour and Welfare (Japan). 2015.
22. Hashizume M, Nakagomi T, Nakagomi O. An Early Detection of Decline in Rotavirus Cases during the 2013/2014 Season in Japan as Revealed by Time-series Analysis of National Surveillance Data. Trop Med Health. 2015;43(3):177–81. doi: 10.2149/tmh.2015-23 [PMC free article] [PubMed]
23. Araki K, Hara M, Sakanishi Y, Shimanoe C, Nishida Y, Matsuo M, et al. Estimating rotavirus vaccine effectiveness in Japan using a screening method. Hum Vaccin Immunother. 2015;12(5):1–6. [PMC free article] [PubMed]
24. Thongprachum A, Khamrin P, Pham NTK, Takanashi S, Okitsu S, Shimizu H, et al. Multiplex RT-PCR for Rapid Detection of Viruses Commonly Causing Diarrhea in Pediatric Patients. J Med Virol. 2017;89:818–24. doi: 10.1002/jmv.24711 [PubMed]
25. Gautam R, Esona MD, Mijatovic-Rustempasic S, Ian Tam K, Gentsch JR, Bowen MD. Real-time RT-PCR assays to differentiate wild-type group A rotavirus strains from Rotarix(®) and RotaTeq(®) vaccine strains in stool samples. Hum Vaccin Immunother. 2014;10(3):767–77. doi: 10.4161/hv.27388 [PMC free article] [PubMed]
26. Jothikumar N, Kang G, Hill VR. Broadly reactive TaqMan assay for real-time RT-PCR detection of rotavirus in clinical and environmental samples. J Virol Methods. 2009;155(2):126–31. doi: 10.1016/j.jviromet.2008.09.025 [PubMed]
27. Kageyama T, Kojima S, Shinohara M, Uchida K, Fukushi S, Hoshino FB, et al. Broadly reactive and highly sensitive assay for Norwalk-like viruses based on real-time quantitative reverse transcription-PCR. J Clin Microbiol. 2003;41(4):1548–57. doi: 10.1128/JCM.41.4.1548-1557.2003 [PMC free article] [PubMed]
28. Nijhuis M, Van Maarseveen N, Schuurman R, Verkuijlen S, De Vos M, Hendriksen K, et al. Rapid and Sensitive Routine Detection of All Members of the Genus Enterovirus in Different Clinical Specimens by Real-Time PCR. J Clin Microbiol. 2002;40(10):3666–70. doi: 10.1128/JCM.40.10.3666-3670.2002 [PMC free article] [PubMed]
29. Fujii Y, Shimoike T, Takagi H, Murakami K, Todaka-Takai R, Park Y, et al. Amplification of all 11 RNA segments of group A rotaviruses based on reverse transcription polymerase chain reaction. Microbiol Immunol. 2012;56(9):630–8. doi: 10.1111/j.1348-0421.2012.00479.x [PubMed]
30. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. doi: 10.1093/molbev/mst197 [PMC free article] [PubMed]
31. Hall T, Biosciences I, Carlsbad C. BioEdit: An important software for molecular biology. GERF Bull Biosci. 2011;2(June):60–1.
32. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402. [PMC free article] [PubMed]
33. Yamashita A, Sekizuka T, Kuroda M. VirusTAP: Viral Genome-Targeted Assembly Pipeline. Front Microbiol. 2016;7(32):1–5. [PMC free article] [PubMed]
34. Yamamoto SP, Kaida A, Kubo H, Iritani N. Gastroenteritis outbreaks caused by a DS-1-like G1P[8] Rotavirus strain, Japan, 2012–2013. Emerg Infect Dis. 2014;20(6):1030–3. doi: 10.3201/eid2006.131326 [PMC free article] [PubMed]
35. Zeller M, Donato C, Trovão NS, Cowley D, Heylen E, Donker NC, et al. Genome-wide evolutionary analyses of G1P[8] strains isolated before and after rotavirus vaccine introduction. Genome Biol Evol. 2015;7(9):2473–83. doi: 10.1093/gbe/evv157 [PMC free article] [PubMed]
36. Sims D, Sudbery I, Ilott NE, Heger A, Ponting CP. Sequencing depth and coverage: key considerations in genomic analyses. Nat Rev Genet. 2014;15(2):121–32. doi: 10.1038/nrg3642 [PubMed]
37. Hill HA, Elam-Evans LD, Yankey D, Singleton JA, Kolasa M. National, State, and Selected Local Area Vaccination Coverage Among Children Aged 19–35 Months—United States, 2014. Morb Mortal Wkly Reports. 2015;64(33):889–96. [PubMed]
38. Elfving K, Andersson M, Msellem MI, Welinder-Olsson C, Petzold M, Björkman A, et al. Real-time PCR threshold cycle cutoffs help to identify agents causing acute childhood diarrhea in Zanzibar. J Clin Microbiol. 2014;52(3):916–23. doi: 10.1128/JCM.02697-13 [PMC free article] [PubMed]
39. Kawada J-I, Okuno Y, Torii Y, Okada R, Hayano S, Ando S, et al. Identification of Viruses in Cases of Pediatric Acute Encephalitis and Encephalopathy Using Next-Generation Sequencing. Sci Rep. 2016;6:33452 doi: 10.1038/srep33452 [PMC free article] [PubMed]
40. Donato CM, Ch’ng LS, Boniface KF, Crawford NW, Buttery JP, Lyon M, et al. Identification of strains of RotaTeq rotavirus vaccine in infants with gastroenteritis following routine vaccination. J Infect Dis. 2012;206(3):377–83. doi: 10.1093/infdis/jis361 [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science