Continuous detection and genetic diversity of human rotavirus A in sewage in eastern China, 2013–2014
Virology Journal volume 13, Article number: 153 (2016)
Rotavirus is the leading viral agent for pediatric gastroenteritis. However, the case-based surveillance for rotavirus is limited in China, and its circulation in the environment is not well investigated.
From 2013 to 2014, rotavirus was detected in raw sewage samples of Jinan and Linyi by quantitative PCR (qPCR) and conventional reverse transcription PCR (RT-PCR). After sequenced and genotyped, sequences analysis was conducted.
A total of 46 sewage samples were collected monthly for the detection of rotavirus, and rotavirus was positive in 43 samples (93.5 %, 43/46). By quantitative assessment, the concentrations of rotavirus in raw sewage ranged from 4.1 × 103 to 1.3 × 106 genome copies (GC)/L in Jinan, and from 1.5 × 103 to 3.0 × 105 GC/L in Linyi. A total of 318 sequences of 5 G-genotypes and 318 sequences of 5 P-genotypes were obtained. G9 (91.8 %, 292/318) and P (56.0 %, 178/318) were the most common G- and P-genotype, respectively. Multiple transmission lineages were recognized in these genotypes. Interestingly, an intragenic recombination event between two G9 lineages was observed.
This study provided the first report of comprehensive environmental surveillance for rotavirus in China. The results suggest that the concentration of rotavirus in raw sewage was high, and multiple rotavirus transmission lineages continuously co-circulated in Shandong.
Rotavirus is the leading viral agent attributable to pediatric acute gastroenteritis, and belongs to genus Rotavirus, family Reoviridae. Its genome contains 11 double-stranded RNA (dsRNA) segments, encoding six structural proteins (VPs) and five (or six) nonstructural proteins (NSPs). Rotavirus is divided into seven groups, and group A, B and C can infect humans. Of these, rotavirus group A (RVA) is predominant, which contains at least 27 G- and 37 P- genotypes based on the VP7 and VP4 sequences . It was estimated that rotavirus infections led to about 0.45 million deaths in children younger than 5 years in 2008 globally, accounting for 37 % of deaths caused by diarrhea . Thus, long-term surveillance of rotavirus is important.
Currently, rotavirus surveillance was mostly carried out in clinical cases. However, asymptomatic rotavirus infection is common , and case-based surveillance can only focus on the symptomatic cases. The epidemiology of rotavirus in asymptomatic individuals was not well investigated. Because rotavirus particles can be excreted to urban sewage system from both symptomatic and asymptomatic individuals, we can obtain more comprehensive molecular characterizations of rotavirus by sewage analysis. Previous studies have proved the feasibility and effectiveness of this novel approach for rotavirus surveillance [4, 5]. In addition, the quantification of rotavirus based on quantitative PCR (qPCR) in sewage has been conducted to better evaluate the seasonal profile of rotavirus and perform risk assessment [4, 6].
In China, rotavirus was the most frequently detected pathogen among children with diarrhea . G3, G1, P and P were the most prevalent G- and P-genotypes . Though the Lanzhou Lamb rotavirus (LLR) vaccine consisting of the monovalent serotype G10P has been licensed since 2001, it was not included in China’s Expanded Program on Immunization. There were still more than 50,000 deaths of rotavirus-specific children under the age of 5 years from 2003 to 2012 in China . Moreover, in China, case-based rotavirus surveillance is limited, and environmental surveillance for rotavirus is rare. The concentration of rotavirus in raw sewage is still unknown. In view of these, based on the successful practices of sewage surveillance for other enteric viruses in China since 2008 [10–13], we detected rotavirus in raw sewage of Jinan and Linyi in Shandong Province from 2013 to 2014, to obtain the epidemiological characterizations of rotavirus circulating in the local population and quantitatively assess the rotavirus in raw sewage in eastern China.
In 2013 and 2014, raw sewage samples were collected monthly in the inlet collector canal of wastewater treatment plants (WWTPs) in Jinan and Linyi of Shandong Province, eastern China, with metropolitan populations of 2.6 million and 1.9 million, respectively. Each sample was transported, stored at a cold temperature (4 °C) and processed within 24 h.
An adsorption-elution method was used to concentrate viruses from sewage samples as described previously . Briefly, sewage sample (1 L) was centrifuged at 3000 × g for 30 min at 4 °C. MgCl2 (2.5 mol l−1) was added to the supernatant to a final concentration of 0.05 mol l−1, and the pH was adjusted to 3.5 by hydrochloric acid (0.5 mol l−1). After filtered through a mixed cellulose ester membrane filter (0.45 μm, ADVANTEC, Tokyo, Japan), the membranes were cut into pieces. The viruses were eluted from the membranes, using 10 ml of a 3 % beef extract solution (pH = 9), followed by ultrasonication for 5 min. The eluted solution was centrifuged at 3000 × g for 30 min, and the supernatant was filtered through a 0.2 μm filter (PALL, Ann Arbor, MI, USA) prior to RNA extraction.
Quantification of rotavirus
High Pure Viral Nucleic Acid Large Volume Kit (Roche, Indiana, USA) was used to extract viral RNA from 1 ml of concentrated sewage. Because the samples all tested negative for poliovirus type 1 Sabin strain, concentrated sewage samples were seeded with a known amount of Sabin 1 as a control to monitor the efficiency of RNA extraction and reverse transcription-qPCR (RT-qPCR). RT-qPCR for rotavirus was carried out by QuantiTect Probe RT-PCR Kit (Qiagen, Hilden, Germany). Briefly, the reaction mixture (25 μl) consisted of 12.5 μl 2× QuantiTect Probe RT-PCR Master Mix, 0.4 μmol l−1 final concentration of each primer, 0.2 μmol l−1 final concentration of each probe, 0.25 μl QuantiTect RT Mix and 10 μl template RNA. The primers and probe towards NSP3 gene were previously described by Freeman et al. . Each sample was measured in triplicate. All qPCR reactions were carried out in an ABI 7500 Real-Time System (Applied Biosystems, Foster City, California, USA). To determine the copy numbers of rotavirus genome, a standard curve was conducted by 10-fold serial dilution (101 ~ 109) of plasmid DNA containing the target gene. All of the concentrations of rotavirus presented were transformed to genome copies per litre (GC/L) according to the volumes used for each step of the procedure. To avoid cross contamination, molecular procedures were performed in different separated rooms and a negative control was included in every qPCR run.
In rotavirus-positive sewage samples confirmed by RT-qPCR, conventional RT-PCR was performed with two improved primer sets by using the Access RT-PCR System (Promega, Fitchburg, WI, USA). Primer set VP4F.142S (TAT GCI CCW GTI AMT TGG) and VP4R.805A (ATT GCA TTT CYT TCC AYA AYG) was used for amplification of a 664-nucleotide (nt) VP4 sequence according to nt position 142–805 of strain Human-tc/USA/Wa/1974 segment 4 (JX406750). Primer set VP7F.49S (ATG TAT GGT ATT GAA TAT ACC) and VP7R.932A (ACT TGC CAC CAT YTY TTC C) was used for amplification of a 884-nt VP7 sequence according to nt position 49–932 of strain Human-tc/USA/Wa/1974 segment 9 (JX406755).
Cloning and sequencing
RT-PCR products were analyzed by electrophoresis with 1.5 % agarose gels. All positive products were gel-purified by QIAquick gel extraction kit (Qiagen, Valencia, CA, USA) and ligated into the pGEM®-T Easy vector (Promega). By shock method, the ligation products were transformed into competent Escherichia coli JM109 cells. For each transformation reaction, ten positive clones were selected after blue and white screening. Then, the plasmid was extracted and sequenced with a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA), on an ABI 3130 genetic analyzer (Applied Biosystems). Molecular typing was carried out by using BLAST.
Nucleotide sequence alignments were carried out by BioEdit 18.104.22.168 . Recombination event was analyzed by RDP software package . Mega 5.0 was used to construct the phylogenetic tree by using the neighbor-joining method . Bootstrapping test was performed with 1000 duplicates. Nucleotide sequences used in phylogenetic analysis from this study were deposited in GenBank under accession numbers KU173953-KU174110.
Concentration of rotavirus
During the 2-year period, a total of 46 raw sewage samples (23 in Jinan and 23 in Linyi) were collected. After RT-qPCR detection, rotavirus was positive in 43 samples (93.5 %, 43/46), 24 (100 %, 24/24) in Jinan and 19 (86.4 %, 19/22) in Linyi, respectively. The RNA extraction and RT-qPCR efficiency of Sabin 1 was 74.5 ± 11.3 % (n = 46), and the efficiency for each sample was used to estimate the actual concentration of rotavirus in the sample. Figure 1 shows the monthly concentration of rotavirus in raw sewage. In Jinan, the concentrations of rotavirus ranged from 4.1 × 103 to 1.3 × 106 GC/L, and were relatively high in November 2013 and December 2014. In Linyi, the concentrations of rotavirus ranged from 1.5 × 103 to 3.0 × 105 GC/L, and were relatively high in October 2013.
Distribution of rotavirus genotypes
By conventional RT-PCR, rotavirus was detected in all rotavirus-positive samples confirmed by RT-qPCR. After cloning and sequencing, 318 G-genotypes sequences (200 in Jinan and 118 in Linyi) and 318 P-genotypes sequences (200 in Jinan and 118 in Linyi) were obtained, belonging to 5 G- (G9, G6, G8, G4 and G2) and 5 P-genotypes (P, P, P, P and P), respectively. G9 was the most common G-genotype, accounting for 91.8 % (292/318). Other G-genotypes were all detected at a low prevalence. P and P were two most common P-genotypes, accounting for 56.0 % (178/318) and 35.8 % (114/318). P and P were only found in Jinan, with the detection rate of 3.5 % (7/200) and 0.5 % (1/200). Figure 2 shows the month-by-month relative proportions of genotypes.
Homologous comparison was performed in all obtained sequences from Jinan and Linyi. To reduce workload, Shandong G9, P and P sequences displaying more than 99.5 % nucleotide identities within these genotypes have been deleted. Phylogenetic analysis was conducted between the remaining sequences and reference strains selected randomly from GenBank database.
Phylogenetic relationships of G-genotypes were illustrated in Fig. 3. Shandong G9 sequences can be divided into two clusters, and cluster II only contained Jinan sequences. Two Chinese reference strains from cases in Shanghai (KC200154) and Wuhan (EU708591) displayed close relationship with sequences in cluster I (98.8–99.5 %) and cluster II (98.6–99.2 %), respectively. Another two Chinese reference strains (DQ904517 and KF726066) showed high similarity with Shandong G2 (98.8–99.3 %) and G4 (96.9 %) sequences, respectively. Shandong G6 sequences were found to be more closely related to reference strains from animals (JX442777 and GU937882) than classical human G6 sequences (89.9–91.1 % vs. 81.1–86.4 %). A human G8 reference strain from Germany (GQ414545) and an animal G8 reference strain from Sudan (KC257096) presented similar nucleotide identity to Shandong G8 sequences (89.7–90.2 % vs. 90.6–91.1 %).
Phylogenetic relationships of P-genotypes were shown in Fig. 4. P sequences were located at one main phylogenetic branch except one sequence in Jinan (JN379-8), and great genetic difference was observed in this branch (91.5–99.5 % identity). Within selected P sequences, two clusters were recognized after excluding sequences JN378-5, JN369-7, JN372-9 and JN375-3. Additionally, Chinese reference strains from Wuhan showed close relationship with Shandong P sequences in the major cluster (94.7–99.5 % identity). Shandong P and P sequences can also be separated into two clusters. P sequence in this study displayed low similarity with human [P3] reference strain (76.1 %), but high similarity with rhesus rotavirus reference strains (87.2–87.8 %).
After RDP analysis and the alignment of 884-nt G9 sequences, an intragenic recombination event was observed in one G9 sequence (JN344-9) (Fig. 5). The crossover sites ranged from position 688 to 698 according to the sequence of strain Human-tc/USA/Wa/1974 segment 9 (JX406755).
Environmental sewage surveillance has been successfully conducted to obtain molecular characterizations of human enteric viruses in some regions [19–21]. In this study, this novel surveillance method was performed to monitor the circulation of rotavirus in Jinan and Linyi of eastern China. During the 2-year period of monitoring, a large number of rotavirus sequences were obtained and revealed the molecular characterizations of rotaviruses circulating in the local population.
At present, various detection rates of rotavirus in urban sewage have been reported [22, 23]. In this study, rotavirus was detected at a high rate (93.5 %) both by conventional RT-PCR and RT-qPCR, reflecting the high prevalence of rotavirus in the local population. Sewage surveillance has been discussed to be a useful tool to evaluate the circulation of rotavirus after rotavirus vaccination, and low detection rates of rotavirus in sewage were observed in some regions with universal rotavirus vaccination [24, 25]. In view of the high detection rate of rotavirus in this study, rotavirus vaccination should be included in the National Immunization Program in China. Rotavirus infection is seasonal, and shows a peak in cold seasons . In this study, RT-qPCR results showed that the concentrations of rotavirus in raw sewage tended to peak in autumn and winter months in Jinan, revealing that RT-qPCR could be an effective approach to exploring the seasonal profile of rotavirus in the population. In Linyi, no apparent seasonality was observed both by conventional RT-PCR and RT-qPCR. This phenomenon might indicate that asymptomatic rotavirus infection was common in this city [5, 23]. The infectious dose of rotavirus is low (<100 virus particles) . The relative abundance of rotavirus in sewage in this study shows that the public health will be threatened if drinking water or food is contaminated by sewage, although the actual level of rotavirus with infectivity may be overestimated.
Currently, G1-G4, G9, P, P and P are the prevalent genotypes worldwide . In this study, G9 and P were recognized predominantly, showing that they were the dominant genotypes in this region. This data also suggested that G9P might be the most prevalent G-P genotypes combination circulating in the local population, as previously described in other regions of China [29, 30]. As two most frequently detected G-genotypes of China in the past decades, G1 and G3 were not recognized in this study, and the detection rates of common G2 and G4 were also low, which may indicate that common G1-G4 were currently circulating at a very low prevalence in the two cities.
G9 infections are globally increasing at the exponential growth since 1993, and G9 have recently been the significant genotype in some regions. In this study, G9 was also observed as the predominant G-genotype and two main phylogenetic clusters were identified. The minor cluster was only recognized in Jinan in 2013, showing its local circulation in Jinan. The major cluster was continuously observed in Jinan and Linyi from 2013 to 2014, revealing that it was the main G9 transmission lineage in the two cities and the inter-city spread may exist. Interestingly, an intragenic recombination event between two G9 lineages was observed in this study. Intragenic recombination in rotavirus is rare until increasing reports recently [31, 32]. To our knowledge, this is the first report of intragenic recombination in G9 rotaviruses. Mixed infections between different variants belonging to the same genotype were discussed as the cause of intragenic recombination . Thus, co-circulation of two G9 transmission lineages in Jinan in 2013 might lead to the occurrence of this recombination event.
G6 and G8 strains were reported to be of animal origin [34, 35]. In this study, some G6 and G8 sequences were obtained. By BLAST, we found that these sequences showed the highest identities to animal sequences (although 10 % divergence is quite different), and no human rotavirus sequences showed higher identities to them. We cannot confirm that these sequences were from humans or animals. They also might be the novel variants circulating in the local population. All of these need to be confirmed by case-based surveillance in the future.
The genetic differences within P, P, P and P sequences were high, and more than one cluster was recognized, suggesting the co-circulation of multiple transmission lineages in these genotypes. P sequence in this study showed the highest identity to the rhesus rotavirus P strain (RRV strain, M18736). RRV is the recipient virus in human-rhesus reassortant vaccine , but the detection of RRV-like P-genotype is rare so far. This is the first recognition of it in sewage. Nevertheless, the meaning of this finding is unknown because there is no evidence to support the epidemiological relevance of RRV P in humans . Further research on rare RRV P is needed.
In conclusion, this study shows that the concentration of rotavirus in raw sewage was high, and multiple rotavirus transmission lineages continuously co-circulated in Jinan and Linyi, eastern China. Additionally, this study proves that environmental surveillance is an ideal measure to obtain epidemiological data on circulating viruses in a given population.
Lanzhou lamb rotavirus
Reverse transcription-quantitative PCR
Rotavirus group A
Wastewater treatment plants
Desselberger U. Rotaviruses. Virus Res. 2014;190:75–96.
Tate JE, Burton AH, Boschi-Pinto C, Steele AD, Duque J, Parashar UD, et al. 2008 estimate of worldwide rotavirus-associated mortality in children younger than 5 years before the introduction of universal rotavirus vaccination programmes: a systematic review and meta-analysis. Lancet Infect Dis. 2012;12:136–41.
Phillips G, Lopman B, Rodrigues LC, Tam CC. Asymptomatic rotavirus infections in England: prevalence, characteristics, and risk factors. Am J Epidemiol. 2010;171:1023–30.
Fumian TM, Leite JP, Rose TL, Prado T, Miagostovich MP. One year environmental surveillance of rotavirus specie A (RVA) genotypes in circulation after the introduction of the Rotarix® vaccine in Rio de Janeiro, Brazil. Water Res. 2011;45:5755–63.
Ruggeri FM, Bonomo P, Ianiro G, Battistone A, Delogu R, Germinario C, et al. Rotavirus genotypes in sewage treatment plants and in children hospitalized with acute diarrhea in Italy in 2010 and 2011. Appl Environ Microbiol. 2015;81:241–9.
Victoria M, Tort LF, García M, Lizasoain A, Maya L, Leite JP, et al. Assessment of gastroenteric viruses from wastewater directly discharged into Uruguay River, Uruguay. Food Environ Virol. 2014;6:116–24.
Yu J, Jing H, Lai S, Xu W, Li M, Wu J, et al. Etiology of diarrhea among children under the age five in China: results from a five-year surveillance. J Infect. 2015;71:19–27.
Liu N, Xu Z, Li D, Zhang Q, Wang H, Duan ZJ. Update on the disease burden and circulating strains of rotavirus in China: a systematic review and meta-analysis. Vaccine. 2014;32:4369–75.
Zhang J, Duan Z, Payne DC, Yen C, Pan X, Chang Z, et al. Rotavirus-specific and overall diarrhea mortality in Chinese children younger than 5 years: 2003 to 2012. Pediatr Infect Dis J. 2015;34:e233–7.
Tao Z, Wang H, Li Y, Xu A, Zhang Y, Song L, et al. Cocirculation of two transmission lineages of echovirus 6 in jinan, china, as revealed by environmental surveillance and sequence analysis. Appl Environ Microbiol. 2011;77:3786–92.
Tao Z, Song Y, Wang H, Zhang Y, Yoshida H, Ji S, et al. Intercity spread of echovirus 6 in Shandong Province, China: application of environmental surveillance in tracing circulating enteroviruses. Appl Environ Microbiol. 2012;78:6946–53.
Wang H, Tao Z, Li Y, Lin X, Yoshida H, Song L, et al. Environmental surveillance of human enteroviruses in Shandong Province, China, 2008 to 2012: serotypes, temporal fluctuation, and molecular epidemiology. Appl Environ Microbiol. 2014;80:4683–91.
Zhou N, Lin X, Wang S, Wang H, Li W, Tao Z, et al. Environmental surveillance for human astrovirus in Shandong Province, China in 2013. Sci Rep. 2014;4:7539.
Lodder WJ, de Roda Husman AM. Presence of noroviruses and other enteric viruses in sewage and surface waters in The Netherlands. Appl Environ Microbiol. 2005;71:1453–61.
Freeman MM, Kerin T, Hull J, McCaustland K, Gentsch J. Enhancement of detection and quantification of rotavirus in stool using a modified real-time RT-PCR assay. J Med Virol. 2008;153:238–40.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser (Oxf). 1999;41:95–8.
Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:vev003.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Kamel AH, Ali MA, El-Nady HG, Aho S, Pothier P, Belliot G. Evidence of the co-circulation of enteric viruses in sewage and in the population of Greater Cairo. J Appl Microbiol. 2010;108:1620–9.
Kokkinos PA, Ziros PG, Mpalasopoulou A, Galanis A, Vantarakis A. Molecular detection of multiple viral targets in untreated urban sewage from Greece. Virol J. 2011;8:195.
Hata A, Katayama H, Kitajima M, Furumai H. Wastewater analysis indicates that genetically diverse astroviruses, including strains belonging to novel clades MLB and VA, are circulating within Japanese Populations. Appl Environ Microbiol. 2015;81:4932–9.
Rodríguez-Díaz J, Querales L, Caraballo L, Vizzi E, Liprandi F, Takiff H, et al. Detection and characterization of waterborne gastroenteritis viruses in urban sewage and sewage-polluted river waters in Caracas, Venezuela. Appl Environ Microbiol. 2009;75:387–95.
Barril PA, Fumian TM, Prez VE, Gil PI, Martínez LC, Giordano MO, et al. Rotavirus seasonality in urban sewage from Argentina: effect of meteorological variables on the viral load and the genetic diversity. Environ Res. 2015;138:409–15.
Bucardo F, Lindgren PE, Svensson L, Nordgren J. Low prevalence of rotavirus and high prevalence of norovirus in hospital and community wastewater after introduction of rotavirus vaccine in Nicaragua. PLoS One. 2011;6:e25962.
O’Ryan ML, Lucero Y, Vidal R. Enteric viruses in wastewaters: an interesting approach to evaluate the potential impact of rotavirus vaccination on viral circulation. Expert Rev Vaccines. 2012;11:419–22.
Parashar UD, Bresee JS, Gentsch JR, Glass RI. Rotavirus. Emerg Infect Dis. 1998;4:561–70.
Glass RI, Parashar UD, Bresee JS, Turcios R, Fischer TK, Widdowson MA, et al. Rotavirus vaccines: current prospects and future challenges. Lancet. 2006;368:323–32.
Bányai K, László B, Duque J, Steele AD, Nelson EA, Gentsch JR, et al. Systematic review of regional and temporal trends in global rotavirus strain diversity in the pre rotavirus vaccine era: insights for understanding the impact of rotavirus vaccination programs. Vaccine. 2012;30S:A122–30.
Zhang S, Yin J, Yang J, Tian L, Li D, Zhang Q, et al. Epidemiology and genetic diversity of group A rotavirus in acute diarrhea patients in pre-vaccination era in southwest China. J Med Virol. 2016; doi: 10.1002/jmv.24606.
Zhang J, Liu H, Jia L, Payne DC, Hall AJ, Xu Z, et al. Active, population-based surveillance for rotavirus gastroenteritis in Chinese children: Beijing Municipality and Gansu Province, China. Pediatr Infect Dis J. 2015;34:40–6.
Parra GI, Bok K, Martínez M, Gomez JA. Evidence of rotavirus intragenic recombination between two sublineages of the same genotype. J Gen Virol. 2004;85:1713–6.
Phan TG, Okitsu S, Maneekarn N, Ushijima H. Evidence of intragenic recombination in G1 rotavirus VP7 genes. J Virol. 2007;81:10188–94.
Jain V, Das BK, Bhan MK, Glass RI, Gentsch JR, Indian Strain Surveillance Collaborating Laboratories. Great diversity of group A rotavirus strains and high prevalence of mixed rotavirus infections in India. J Clin Microbiol. 2001;39:3524–9.
Gerna G, Sarasini A, Parea M, Arista S, Miranda P, Brüssow H, et al. Isolation and characterization of two distinct human rotavirus strains with G6 specificity. J Clin Microbiol. 1992;30:9–16.
Martella V, Bányai K, Matthijnssens J, Buonavoglia C, Ciarlet M. Zoonotic aspects of rotaviruses. Vet Microbiol. 2010;140:246–55.
Midthun K, Greenberg HB, Hoshino Y, Kapikian AZ, Wyatt RG, Chanock RM. Reassortant rotaviruses as potential live rotavirus vaccine candidates. J Virol. 1985;53:949–54.
Martella V, Ciarlet M, Pratelli A, Arista S, Terio V, Elia G, et al. Molecular analysis of the VP7, VP4, VP6, NSP4, and NSP5/6 genes of a buffalo rotavirus strain: identification of the rare P rhesus rotavirus-like VP4 gene allele. J Clin Microbiol. 2003;41:5665–75.
We thank Wenfeng Li and Peng Chen for their contributions to the study.
This study was supported by a grant from the National Natural Science Foundation of China (grant number: 81573209), a grant from Taishan Scholar Program of Shandong Province (grant number: ts201511105), and a grant from Health and Family Planning Commission of Shandong Province (grant number: 2013WS0158).
Availability of data and material
All data generated or analyzed during this study are included in this published article.
AX designed the study. SW, XL, PW and HZ conducted sampling and experiments. ZB, HW and PH supervised the project. NZ and DL conducted sequences analysis and wrote this manuscript. ZT and YS corrected and edited the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
About this article
Cite this article
Zhou, N., Lv, D., Wang, S. et al. Continuous detection and genetic diversity of human rotavirus A in sewage in eastern China, 2013–2014. Virol J 13, 153 (2016). https://doi.org/10.1186/s12985-016-0609-0