Next Article in Journal
A Turkey-origin H9N2 Avian Influenza Virus Shows Low Pathogenicity but Different Within-host Diversity in Experimentally Infected Turkeys, Quail and Ducks
Previous Article in Journal
Construction and Immunogenicity of Virus-Like Particles of Feline Parvovirus from the Tiger
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 3
10.3390/v12030318
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Prevalence of Rotavirus A in Raw Sewage Samples from Northeast Spain

by
Marcelle Silva-Sales
*,†,
Sandra Martínez-Puchol
,
Eloy Gonzales-Gustavson
,
Ayalkibet Hundesa
and
Rosina Gironès
Microbiology, Virology and Biotechnology Section, Department of Genetics, Microbiology and Statistics, Faculty of Biology, University of Barcelona, Av. Diagonal 643, Barcelona 08028, Spain
*
Author to whom correspondence should be addressed.
†
Present address: Laboratory of Virology and Cell Culture, Department of Bioscience and Technology, Institute of Tropical Pathology and Public Health, Federal University of Goiás; Rua 235 s/n, Setor Leste Universitário, 74605-050 Goiânia, Goiás, Brazil.
Viruses 2020, 12(3), 318; https://doi.org/10.3390/v12030318
Submission received: 28 January 2020 / Revised: 8 March 2020 / Accepted: 14 March 2020 / Published: 16 March 2020
Browse Figures
Versions Notes

Abstract

:
Rotavirus A (RVA) is the most common virus associated with infantile gastroenteritis worldwide, being a public health threat, as it is excreted in large amounts in stool and can persist in the environment for extended periods. In this study, we performed the detection of RVA and human adenovirus (HAdV) by TaqMan qPCR and assessed the circulation of RVA genotypes in three wastewater treatment plants (WWTPs) between 2015 and 2016 in Catalonia, Spain. RVA was detected in 90% and HAdV in 100% of the WWTP samples, with viral loads ranging between 3.96 × 104 and 3.30 × 108 RT-PCR Units/L and 9.51 × 104 and 1.16 × 106 genomic copies/L, respectively. RVA VP7 and VP4 gene analysis revealed the circulation of G2, G3, G9, G12, P[4], P[8], P[9] and P[10]. Nucleotide sequencing (VP6 fragment) showed the circulation of I1 and I2 genotypes, commonly associated with human, bovine and porcine strains. It is important to mention that the RVA strains isolated from the WWTPs were different from those recovered from piglets and calves living in the same area of single sampling in 2016. These data highlight the importance of monitoring water matrices for RVA epidemiology and may be a useful tool to evaluate and predict possible emergence/reemergence of uncommon strains in a region.

Graphical Abstract

1. Introduction

Rotavirus A (RVA) is the most common pathogen associated with acute gastroenteritis in children worldwide. It accounts for approximately 215,000 deaths annually of children up to 5 years old [1]. In Europe, RVA is responsible for 75,000–150,000 infantile hospitalization due to acute gastroenteritis associated with RVA infections. In Spain, the annual incidence of acute gastroenteritis associated with RVA in primary care ranges between 15.4 and 19.5 cases per 1000 children up to 5 years and 20 cases per 1000 children up to 3 years [2].
RVA is one of the nine groups (A–I) of rotavirus established according to the International Committee on Taxonomy of Viruses (ICTV) (ICTV, updated in July 2018 and available at https://talk.ictvonline.org/taxonomy/). The virion is non-enveloped, icosahedral, with a triple-layered protein capsid and consists of 11 segments of double-stranded genomic RNA (dsRNA), which encode six structural (VP1–4, 6 and 7) and six non-structural (NSP1–6) proteins. The binary classification system for RVA is based on molecular analysis of the two outer capsid proteins, VP7 and VP4, which define the G (for glycoprotein) and P (for protease sensitive) genotypes, respectively. To date, there are 36 G genotypes and 51 P genotypes described in the literature (Rotavirus Classification Working Group, updated in May 2018 and available at https://rega.kuleuven.be/cev/viralmetagenomics/virus-classification/rcwg). Six G/P combinations: G1P[8], G2P[4], G3P[8], G4P[8], G9P[8] and G12P[8] are the most prevalent combinations detected in humans [3,4,5,6,7,8].
However, RVA genotypes derived from animals have also been identified in humans, demonstrating that animals can act as a source of infection for humans. It is believed that animal RVA that crosses the interspecies barrier is not able to infect or disseminate efficiently in the new host. However, if such strains acquired gene segments of human origin, this would increase the chance of efficiently infecting humans [9].
With an average population of 150 million pigs, the European Union is the world’s second largest producer of pig meat (with the countries that produce the most being France, Germany and Spain) after China; it is also the largest exporter of pig meat and pig meat products (accessed at https://ec.europa.eu/agriculture/pigmeat_en, on 14th October 2018). The Spanish region of Catalonia is responsible for 44% of all pig meat produced in the country; its production is greater than that of Italy, Belgium and the Netherlands, and it exports more than Brazil, France or Poland (accessed at http://generalitatgirona.gencat.cat/ca/detalls/Noticia/Nova-Noticia-04179, on 30th January 2017).
RVA is excreted at high concentration in stools from infected humans and animals (1 × 1011 virus particles per gram of stool), is very stable and resistant to different methods of wastewater treatments, thus it can persist in the environment for extended periods serving as a source of infection for susceptible humans and animals [10,11].
Vaccination is considered the most efficient strategy to reduce the impact of RVA infection. Since 2006, two anti-RVA vaccines have been licensed worldwide, Rotarix® and Rotateq®. Rotarix® (RV1) (GlaxoSmithKline, Rixensart, Belgium) is an attenuated monovalent vaccine, derived from a human G1P[8] strain, which is administered orally to infants at 2 and 4 months of age. Rotateq® (RV5) (Merck and Co., Whitehouse Station, NJ, USA; Sanofi Pasteur MSD, Lyon, France) is a pentavalent vaccine comprising five human-bovine reassortant strains: four expressing human G1, G2, G3, G4 and bovine P7[5], and one expressing P1A[8] (human) and G6 (bovine). This latter vaccine is administered orally to infants at 2, 4 and 6 months of age [12]. Both vaccines were licensed in Spain in late 2006, and although the Advisory Board on Vaccines of the Spanish Association of Pediatricians began to recommend vaccination against RVA in all eligible children in 2008, RV1 and RV5 are not funded by the public health system, being only available in the private market [13].
This study aims to increase the availability of information concerning the circulation of RVA, based on data obtained from three wastewater treatment plants (WWTPs) in Catalonia, Spain, over a year (April 2015–May 2016). We also report our analysis of RVA in animal samples (porcine and bovine) collected at three local farms, in the same area of the WWTP sample collection, to determinate whether these animals may contribute as a source of contamination for the local population. Due to the role of human adenovirus (HAdV) as an indicator of human fecal contamination in water matrices, HAdV was considered as an alternative internal control over the studied samples.

2. Material and Methods

2.1. Raw Sewage Samples

Raw sewage samples were collected over a 1-year period between April 2015 and May 2016, from three WWTPs located in Catalonia, Spain. Among them was EDAR Sant Adrià de Besòs, with a capacity of 525,000 m3/day and which receives sewage from more than half the population of the Barcelona metropolitan area in the municipalities of Badalona, Barcelona (three quarters), Montgat, Sant Adrià de Besòs, Santa Coloma de Gramenet, Tiana and part of Montcada and Reixac, comprising around 2,800,000 people. Another source was EDAR Baix Llobregat, with a capacity of 315,000 m3/day, which receives sewage from Barcelona (one quarter), Cornellà de Llobregat, El Prat de Llobregat, Esplugues de Llobregat, L’Hospitalet de Llobregat, Sant Joan Despi, Sant Boi de Llobregat (partially), Santa Coloma de Cervelló and Saint Just Desvern (partially), equivalent to more than 1,700,000 inhabitants.
The remaining plant was EDAR Granollers (Vallès Oriental), located in the town of Granollers, with a capacity of 30,000 m3/day and which receives sewage from the municipalities of Canovelles, Les Franqueses of the Vallès and Granollers, with a population of 112,154 inhabitants. A total of 20 raw sewage samples were collected over the four annual seasons. At each sampling point, a 42 mL raw sewage sample was collected in sterile plastic bottles, kept at 4 °C and transported to the laboratory for immediate analysis.

2.2. Animal Fecal Samples

In order to estimate the patterns of circulation of RVA in animals in the region, we collected 1 fecal sample from a 7-week-old diarrheic piglet and 8 samples pools (3 samples in each pool) collected from the floor of pens of diarrheic and non-diarrheic pig and cattle of different ages, in the city of Vic, Catalonia. To note, this city is known for its production of sausages and other pork derivatives. In all visited farms, diarrheal animals are kept separately from healthy animals. The samples were collected in 50 mL sterile plastic bottles, kept at 4 °C and transported to the laboratory for immediate analysis.

2.3. Virus Concentration

The samples were concentrated using the ultracentrifugation method as previously described [14]. To avoid false negative results and to evaluate the presence of inhibitors, each sample was inoculated with 250 µL of bacteriophage MS2 (1 × 106 gene copies—GC) as an internal control before the concentration assay and detected by Taqman® RT-qPCR for viral recovery analysis [15].

2.4. Nucleic Acid Extraction, Reverse Transcription (RT) and Quantitative PCR (qPCR)

Viral nucleic acid was extracted using the commercial QIAmp Viral RNA kit® (Qiagen Inc., Hilden, Germany). The RT-qPCR was performed for detection and quantification of RVA as previously described [16]. A standard curve (5 × 100–5 × 107 RT-PCR Units (RT-PCR U)/per reaction) was generated using 10-fold serial dilutions of a 473 bp qBlock® gene fragment specific for the NSP3 gene from RVA (IDT—Integrated DNA Technologies, Coralville, IA, USA). The complementary DNA (cDNA) was obtained for RVA genotyping analysis using pd(N)6 (Random Hexamer, Amersham Biosiences) and SuperScript® III as Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA). HAdV was analyzed by qPCR as described by [17] and was considered as a second internal process control in WWTPs samples, as is considered as an indicator of human fecal contamination in water matrices. The limits of detection (LOD) for RVA and HAdV were 5.29 × 104 RT-PCR U/L and 5.10 × 103 GC/L, respectively.

2.5. VP4, VP6 and VP7 Amplification

Complementary DNA (cDNA) product was used as the template for amplification of the VP6, VP7 and VP8* portion of the VP4 gene segments. RT-PCR was performed using the primers VP6F and VP6R [18], 9Con1/9Con2 [19] and 4Con3/4Con2 [20], to the amplification of genes VP6, VP7 and VP8*, respectively. G (VP7) and P (VP8*) genotyping methods were conducted according to protocols established by the World Health Organization (WHO), available at http://apps.who.int/iris/bitstream/10665/70122/1/WHO_IVB_08.17_eng.pdf.

2.6. Sequencing of the VP6 RT-PCR Products

The 379 bp VP6 amplicons (nt 747–1126) obtained from RVA-positive strains were purified using a DNA Cleaner and Concentrator kit (D4013, Zymo Research Corp., Irvine, CA, USA), following the manufacturer’s recommendations. The purified DNA amplicons were subjected to a sequencing reaction using the Big Dye Terminator Cycle Sequencing Kit version 3.1 (Applied Biosystems, Foster City, CA, USA) and the primers VP6F and VP6R at 3.2 pmol. The BioEdit Sequence Alignment Editor Program® was used to edit and align the sequences obtained [21]. Genetic analysis was performed using the neighbor-joining method with 2000 bootstrap replicates and the genetic distance was calculated using Kimura’s two-parameter matrix. A phylogenetic tree was constructed using MEGA v.7 software [22].

2.7. Statistical Analysis

All data were statistically analyzed using R software® 3.5.1 [23]. First, descriptive statistics of all the quantifications and recoveries were obtained. Additionally, the Kruskal–Wallis non-parametric test was used to evaluate the differences in the RVA and HAdV concentrations in raw sewage between seasons. Post-hoc pairwise comparisons after the Kruskal–Wallis test were performed using the Tukey and Kramer (Nemenyi) test with Tukey-Dist approximation for independent samples with the PMCMR package in R [24]. Finally, Spearman’s rank correlation coefficient was used to evaluate the quantifications observed for HAdV and RVA.

3. Results and Discussion

RVA and HAdV genome concentrations were assessed over a one-year study in raw sewage samples from three WWTPs from the region of Catalonia, Spain. RVA was detected in 18 out of 20 sewage samples (90%) by RT-qPCR, with a mean viral load ranging from 3.96 × 104 to 3.30 × 108 RT-PCR U/L over the studied period (Table 1). The presence of RVA in a large portion of samples indicates the high circulation of this virus in Catalonia. Both high- and low-income countries have indicated the presence of RVA in WWTP samples, thus demonstrating the great dispersion of this virus in the environment worldwide [25,26,27].
All samples were found positive for the bacteriophage MS2 (our internal control), with recovery efficiency ranging from 9.8% to 88.2%. This variation may be justified by the complex nature of the raw sewage samples. Previous studies have shown important variation rates in viral recovery between concentration methods [28,29,30] and once RVA and MS2 do not have the same characteristics (dsRNA vs. bacteriophage ssRNA), the efficiency of recovery for both viruses is not expected to be correlated.
The concentrations of HAdV and RVA varied significantly between different seasons (p-value < 0.000 in both viruses). Both viruses presented a similar pattern of circulation over the seasons. with significant p-values comparing viral detection in winter and spring with that in summer and autumn (Table 2); this suggests greater viral detection in the colder months in the northern hemisphere for both viruses, with a moderate correlation (R = 67.2%, p-value < 0.000) between viral detection and seasonality (Figure 1). In accordance with data on RVA seasonality in Europe, this virus was not detected in the months of June and July (summer) in Barcelona (Figure 2), corroborating previous studies that have shown a large circulation of RVA in the cooler months in temperate countries [31,32]. It is important to emphasize that lack of detection of the virus does not imply that there was no RVA circulation in this period; instead, viral detection was below the LOD, due to the low sensitivity of RT-qPCR for RVA detection.
The VP4 gene was successfully characterized in all samples collected in 2015, showing the circulation of the genotypes P[4], P[8], P[9] and P[10] in Barcelona. The VP7 gene could only be genotyped in samples collected in April 2016 from Besòs and Granollers, showing the circulation of genotypes G2, G3, G9 and G12 in raw sewage from Barcelona (Table 3). These results corroborate a previous study [33], which showed the circulation of the genotypes G2, G3, G9, P[4], P[8] and P[9] in Barcelona; more specifically, at the Sant Adrià de Besòs WWTP. A study conducted between 2010 and 2011 [26] also showed RVA circulation in urban wastewater samples from different Italian cities. Corroborating our present data, those authors detected a large variety of G and P genotypes, such as G1, G2, G3, G4, G6, G9, G26, P[4], P[6], P[8], P[9], P[14] and P[19], as well as combined genotypes: G1+G2-P[4]+P[8].
The G12 genotype was detected at the Granollers WWTP in April 2016. This genotype, previously considered unusual, has been reported in association with cases of hospitalization in several studies worldwide, as it is the 6th most prevalent human RVA genotype [34]. The G12 strain detected in this study was sequenced and the result showed a high degree of nucleotide similarity with other human wild-type G12 strains circulating recently in different continents: Africa, America and Europe.
RVA-positive strains, by the qPCR method, were subjected to VP6 amplification and Sanger sequencing. The results showed the circulation of the genotypes I1 and I2 in Catalonia (Figure 3). Samples collected in 2015 belonged to the I2 genotype and were grouped in a different cluster together with human and animal (antelope, cow and sheep) wild-type strains. On the other hand, samples collected in 2016 belonged to the I1 genotype, grouped together with human wild-type strains of different G/P genotypes circulating worldwide.
The analysis conducted on animal pool samples showed a high level of RVA circulation in different farms (see Supplementary Materials). Seven out of the nine samples (78%) were positive for RVA, either via RT-qPCR or RT-PCR methods. Samples recovered from animals under seven months of age showed positivity for RVA with viral loads ranging from 3.11 × 105 to 7.60 × 108. The two remaining pool samples were recovered from 7-month-old piglets and 8-month-old calves; thus, our results reinforce previous findings in the literature that describe RVA as one of the most common viruses detected in younger animals (1–8 weeks of age) [35,36]. Despite the small number of animals enrolled in our study, the differences concerning in RVA circulating in raw sewage and animals from the studied area, suggest that RVA detected in WWTPs may have a distinct origin than the surveyed animals.
It is important to highlight that the detection of high or low levels of virus in a specific geographic area depends not only on virus excretion levels, but also on the protocol used for detection [37]. We observed that there is no statistical significance between PCR and qPCR results for RVA detection, although previous studies showed greater sensitivity of qPCR in comparison to PCR [25,38]. The VP6 sequencing protocol was previously suggested to be more sensitive for RVA detection in environmental samples than amplification of the VP4 and VP7 genes [39].
The present study showed the high circulation of RVA in sewage samples from Catalonia, suggesting continuous RVA dissemination by individuals, allowing its detection despite the high dilution in sewage. Our results also highlight the importance of RVA surveillance studies in Spain. In our samples, RVA was detected in high concentrations in sewage samples, reflecting the high circulation of the virus in the population, although asymptomatically in most cases. Knowledge of RVA genotype circulation in the environment may help clinical monitoring of RVA infections, predicting which genotypes could be detected in clinical samples during RVA outbreaks.
Despite the relevance to environmental virology regarding RVA spreading in Catalonia in high concentrations, the small number of samples, the presence of possible inhibitors in sewage samples and the variable viral recovery values represent the main limitations of this study. The present data show RVA circulation in young animals (piglets and calves), though the sequences recovered from animal samples were not related to those found in sewage samples. Nevertheless, more studies are necessary to exclude the possibility of RVA infection through the consumption of raw or poorly cooked meat from infected animals.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1999-4915/12/3/318/s1, Table S1: Data on mean RVA load detected from fecal samples recovered from animals.

Author Contributions

Conceptualization, M.S.-S. and R.G.; methodology, M.S.-S. and A.H.; formal analysis, M.S.-S., S.M.-P. and E.G.-G.; investigation, M.S.-S. and A.H.; resources, R.G.; writing—original draft preparation, M.S.-S.; writing—review and editing, M.S.-S., S.M.-P. and E.G.-G.; visualization, M.S.-S.; supervision, R.G.; project administration, M.S.-S. and R.G.; funding acquisition, R.G. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Water Challenges for a Changing World Joint Programming Initiative (award: W2013-095-C03-01), the Spanish Ministry of Economy and Competitiveness (MINECO; award: AGL2014-55081-R) and the Generalitat de Catalunya via award 2014-SGR-914 to the Grup de Recerca Consolidat: Virus, bacteris i protozous d’interès en salut pública, aigua y aliments; VIRBAP.

Acknowledgments

Marcelle Silva-Sales would like to acknowledge the Brazilian National Council for Scientific and Technological Development (CNPq, in Portuguese: Conselho Nacional de Desenvolvimento Científico e Tecnológico) for support via a postdoctoral grant.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tate, J.E.; Burton, A.H.; Boschi-Pinto, C.; Parashar, U.D.; World Health Organization-Coordinated Global Rotavirus Surveillance Network. Global, Regional, and National Estimates of Rotavirus Mortality in Children. Clin. Infect. Dis. 2016, 62 (Suppl. S2), S96–S105. [Google Scholar] [CrossRef] [Green Version]
  2. Diez-Domingo, J.; Garces-Sanchez, M.; Gimenez-Sanchez, F.; Colomina-Rodriguez, J.; Martinon-Torres, F. [What have we learnt about rotavirus in Spain in the last 10 years?]. An. Pediatr. (Barc.) 2019. [Google Scholar] [CrossRef]
  3. Rahman, M.; Matthijnssens, J.; Yang, X.; Delbeke, T.; Arijs, I.; Taniguchi, K.; Iturriza-Gomara, M.; Iftekharuddin, N.; Azim, T.; Van Ranst, M. Evolutionary history and global spread of the emerging g12 human rotaviruses. J. Virol. 2007, 81, 2382–2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Matthijnssens, J.; Van Ranst, M. Genotype constellation and evolution of group A rotaviruses infecting humans. Curr. Opin. Virol. 2012, 2, 426–433. [Google Scholar] [CrossRef] [PubMed]
  5. Trojnar, E.; Sachsenroder, J.; Twardziok, S.; Reetz, J.; Otto, P.H.; Johne, R. Identification of an avian group A rotavirus containing a novel VP4 gene with a close relationship to those of mammalian rotaviruses. J. Gen. Virol. 2013, 94, 136–142. [Google Scholar] [CrossRef] [PubMed]
  6. Bucardo, F.; Rippinger, C.M.; Svensson, L.; Patton, J.T. Vaccine-derived NSP2 segment in rotaviruses from vaccinated children with gastroenteritis in Nicaragua. Infect. Genet. Evol. 2012, 12, 1282–1294. [Google Scholar] [CrossRef] [Green Version]
  7. Delogu, R.; Ianiro, G.; Camilloni, B.; Fiore, L.; Ruggeri, F.M. Unexpected spreading of G12P[8] rotavirus strains among young children in a small area of central Italy. J. Med. Virol. 2015, 87, 1292–1302. [Google Scholar] [CrossRef]
  8. da Silva, M.F.; Fumian, T.M.; de Assis, R.M.; Fialho, A.M.; Carvalho-Costa, F.A.; da Silva Ribeiro de Andrade, J.; Leite, J.P. VP7 and VP8* genetic characterization of group A rotavirus genotype G12P[8]: Emergence and spreading in the Eastern Brazilian coast in 2014. J. Med. Virol. 2017, 89, 64–70. [Google Scholar] [CrossRef]
  9. Martella, V.; Banyai, K.; Matthijnssens, J.; Buonavoglia, C.; Ciarlet, M. Zoonotic aspects of rotaviruses. Vet. Microbiol. 2010, 140, 246–255. [Google Scholar] [CrossRef] [Green Version]
  10. Carter, M.J. Enterically infecting viruses: Pathogenicity, transmission and significance for food and waterborne infection. J. Appl. Microbiol. 2005, 98, 1354–1380. [Google Scholar] [CrossRef]
  11. Bosch, A.; Guix, S.; Sano, D.; Pinto, R.M. New tools for the study and direct surveillance of viral pathogens in water. Curr. Opin. Biotechnol. 2008, 19, 295–301. [Google Scholar] [CrossRef] [PubMed]
  12. Heaton, P.M.; Ciarlet, M. Vaccines: The pentavalent rotavirus vaccine: Discovery to licensure and beyond. Clin. Infect. Dis. 2007, 45, 1618–1624. [Google Scholar] [CrossRef] [PubMed]
  13. Moreno-Perez, D.; Alvarez Garcia, F.J.; Alvarez Aldean, J.; Cilleruelo Ortega, M.J.; Garces Sanchez, M.; Garcia Sanchez, N.; Hernandez Merino, A.; Mendez Hernandez, M.; Merino Moina, M.; Montesdeoca Melian, A.; et al. En representación del Comité Asesor de Vacunas de la Asociación Española de Pediatría (CAV-AEP), Immunisation schedule of the Spanish Association of Paediatrics: 2018 recommendations. An. Pediatr. (Barc.) 2018, 88, 53. [Google Scholar] [CrossRef] [PubMed]
  14. Pina, S.; Jofre, J.; Emerson, S.U.; Purcell, R.H.; Girones, R. Characterization of a strain of infectious hepatitis E virus isolated from sewage in an area where hepatitis E is not endemic. Appl. Environ. Microbiol. 1998, 64, 4485–4488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Calgua, B.; Carratala, A.; Guerrero-Latorre, L.; de Abreu Correa, A.; Kohn, T.; Sommer, R.; Girones, R. UVC Inactivation of dsDNA and ssRNA Viruses in Water: UV Fluences and a qPCR-Based Approach to Evaluate Decay on Viral Infectivity. Food Environ. Virol. 2014, 6, 260–268. [Google Scholar] [CrossRef] [PubMed]
  16. Zeng, S.Q.; Halkosalo, A.; Salminen, M.; Szakal, E.D.; Puustinen, L.; Vesikari, T. One-step quantitative RT-PCR for the detection of rotavirus in acute gastroenteritis. J. Virol. Methods 2008, 153, 238–240. [Google Scholar] [CrossRef]
  17. Bofill-Mas, S.; Albinana-Gimenez, N.; Clemente-Casares, P.; Hundesa, A.; Rodriguez-Manzano, J.; Allard, A.; Calvo, M.; Girones, R. Quantification and stability of human adenoviruses and polyomavirus JCPyV in wastewater matrices. Appl. Environ. Microbiol. 2006, 72, 7894–7896. [Google Scholar] [CrossRef] [Green Version]
  18. Iturriza Gomara, M.; Wong, C.; Blome, S.; Desselberger, U.; Gray, J. Molecular characterization of VP6 genes of human rotavirus isolates: Correlation of genogroups with subgroups and evidence of independent segregation. J. Virol. 2002, 76, 6596–6601. [Google Scholar] [CrossRef] [Green Version]
  19. Gouvea, V.; Glass, R.I.; Woods, P.; Taniguchi, K.; Clark, H.F.; Forrester, B.; Fang, Z.Y. Polymerase chain reaction amplification and typing of rotavirus nucleic acid from stool specimens. J. Clin. Microbiol. 1990, 28, 276–282. [Google Scholar] [CrossRef] [Green Version]
  20. Gentsch, J.R.; Glass, R.I.; Woods, P.; Gouvea, V.; Gorziglia, M.; Flores, J.; Das, B.K.; Bhan, M.K. Identification of group A rotavirus gene 4 types by polymerase chain reaction. J. Clin. Microbiol. 1992, 30, 1365–1373. [Google Scholar] [CrossRef] [Green Version]
  21. Hall, T. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. In Nucleic Acids Symposium Series; Oxford University Press: Oxford, UK, 1999; pp. 95–98. [Google Scholar]
  22. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. R Development Core Team. R: A Language and Environment for Statistical Computing. Computing 2016. [Google Scholar] [CrossRef]
  24. Pohlert, T. The Pairwise Multiple Comparison of Mean Ranks Package (PMCMR). R Package 2014, 27, 9. [Google Scholar]
  25. Fumian, T.M.; Leite, J.P.; Rose, T.L.; Prado, T.; Miagostovich, M.P. One year environmental surveillance of rotavirus specie A (RVA) genotypes in circulation after the introduction of the Rotarix(R) vaccine in Rio de Janeiro, Brazil. Water Res. 2011, 45, 5755–5763. [Google Scholar] [CrossRef] [PubMed]
  26. Ruggeri, F.M.; Bonomo, P.; Ianiro, G.; Battistone, A.; Delogu, R.; Germinario, C.; Chironna, M.; Triassi, M.; Campagnuolo, R.; Cicala, A.; 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–249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Zhou, N.; Lv, D.; Wang, S.; Lin, X.; Bi, Z.; Wang, H.; Wang, P.; Zhang, H.; Tao, Z.; Hou, P.; et al. Continuous detection and genetic diversity of human rotavirus A in sewage in eastern China, 2013–2014. Virol. J. 2016, 13, 153. [Google Scholar] [CrossRef] [Green Version]
  28. Calgua, B.; Rodriguez-Manzano, J.; Hundesa, A.; Sunen, E.; Calvo, M.; Bofill-Mas, S.; Girones, R. New methods for the concentration of viruses from urban sewage using quantitative PCR. J. Virol. Methods 2013, 187, 215–221. [Google Scholar] [CrossRef] [Green Version]
  29. Gonzales-Gustavson, E.; Cardenas-Youngs, Y.; Calvo, M.; da Silva, M.F.; Hundesa, A.; Amoros, I.; Moreno, Y.; Moreno-Mesonero, L.; Rosell, R.; Ganges, L.; et al. Characterization of the efficiency and uncertainty of skimmed milk flocculation for the simultaneous concentration and quantification of water-borne viruses, bacteria and protozoa. J. Microbiol. Methods 2017, 134, 46–53. [Google Scholar] [CrossRef]
  30. Petterson, S.; Grondahl-Rosado, R.; Nilsen, V.; Myrmel, M.; Robertson, L.J. Variability in the recovery of a virus concentration procedure in water: Implications for QMRA. Water Res. 2015, 87, 79–86. [Google Scholar] [CrossRef]
  31. Levy, K.; Hubbard, A.E.; Eisenberg, J.N. Seasonality of rotavirus disease in the tropics: A systematic review and meta-analysis. Int. J. Epidemiol. 2009, 38, 1487–1496. [Google Scholar] [CrossRef] [Green Version]
  32. Patel, M.M.; Pitzer, V.E.; Alonso, W.J.; Vera, D.; Lopman, B.; Tate, J.; Viboud, C.; Parashar, U.D. Global seasonality of rotavirus disease. Pediatr. Infect. Dis. J. 2013, 32, e134–e147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Villena, C.; El-Senousy, W.M.; Abad, F.X.; Pinto, R.M.; Bosch, A. Group A rotavirus in sewage samples from Barcelona and Cairo: Emergence of unusual genotypes. Appl. Environ. Microbiol. 2003, 69, 3919–3923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Matthijnssens, J.; Bilcke, J.; Ciarlet, M.; Martella, V.; Banyai, K.; Rahman, M.; Zeller, M.; Beutels, P.; Van Damme, P.; Van Ranst, M. Rotavirus disease and vaccination: Impact on genotype diversity. Future Microbiol. 2009, 4, 1303–1316. [Google Scholar] [CrossRef] [PubMed]
  35. Miyazaki, A.; Kuga, K.; Suzuki, T.; Kohmoto, M.; Katsuda, K.; Tsunemitsu, H. Annual changes in predominant genotypes of rotavirus A detected in the feces of pigs in various developmental stages raised on a conventional farm. Vet. Microbiol. 2013, 163, 162–166. [Google Scholar] [CrossRef] [PubMed]
  36. Papp, H.; Borzak, R.; Farkas, S.; Kisfali, P.; Lengyel, G.; Molnar, P.; Melegh, B.; Matthijnssens, J.; Jakab, F.; Martella, V.; et al. Zoonotic transmission of reassortant porcine G4P[6] rotaviruses in Hungarian pediatric patients identified sporadically over a 15 year period. Infect. Genet. Evol. 2013, 19, 71–80. [Google Scholar] [CrossRef] [PubMed]
  37. Girones, R.; Ferrus, M.A.; Alonso, J.L.; Rodriguez-Manzano, J.; Calgua, B.; Correa Ade, A.; Hundesa, A.; Carratala, A.; Bofill-Mas, S. Molecular detection of pathogens in water--the pros and cons of molecular techniques. Water Res. 2010, 44, 4325–4339. [Google Scholar] [CrossRef]
  38. Vieira, C.B.; Mendes, A.C.; Guimaraes, F.R.; Fumian, T.M.; Leite, J.P.; Gaspar, A.M.; Miagostovich, M.P. Detection of enteric viruses in recreational waters of an urban lagoon in the city of Rio de Janeiro, Brazil. Mem. Inst. Oswaldo Cruz 2012, 107, 778–784. [Google Scholar] [CrossRef] [Green Version]
  39. Ferreira, F.F.; Guimaraes, F.R.; Fumian, T.M.; Victoria, M.; Vieira, C.B.; Luz, S.; Shubo, T.; Leite, J.P.; Miagostovich, M.P. Environmental dissemination of group A rotavirus: P-type, G-type and subgroup characterization. Water Sci. Technol. 2009, 60, 633–642. [Google Scholar] [CrossRef]
Viruses 12 00318 g001 550
Figure 1. Correlation analysis between human adenovirus (HAdV) and rotavirus A (RVA) circulation in Catalonia. HAdV and RVA were detected and viral loads were estimated through qPCR (HAdV) and RT-qPCR (RVA) analysis of 42 mL raw sewage samples. The correlation analysis was performed using R software package.
Figure 1. Correlation analysis between human adenovirus (HAdV) and rotavirus A (RVA) circulation in Catalonia. HAdV and RVA were detected and viral loads were estimated through qPCR (HAdV) and RT-qPCR (RVA) analysis of 42 mL raw sewage samples. The correlation analysis was performed using R software package.
Viruses 12 00318 g001
Viruses 12 00318 g002 550
Figure 2. Human adenovirus (HAdV) and rotavirus A (RVA) circulation in raw sewage samples by season. Seasonal distribution of HAdV and RVA detected by qPCR and RT-qPCR, respectively, in raw sewage samples collected between spring 2015 and spring 2016 in Catalonia.
Figure 2. Human adenovirus (HAdV) and rotavirus A (RVA) circulation in raw sewage samples by season. Seasonal distribution of HAdV and RVA detected by qPCR and RT-qPCR, respectively, in raw sewage samples collected between spring 2015 and spring 2016 in Catalonia.
Viruses 12 00318 g002
Viruses 12 00318 g003 550
Figure 3. Phylogenetic analysis of the rotavirus A (RVA) VP6 gene from raw sewage and animal samples. The raw sewage samples from the present study are indicated with a black dot followed by month, year and wastewater treatment plant of collection. Samples recovered from animals are indicated with white diamonds.
Figure 3. Phylogenetic analysis of the rotavirus A (RVA) VP6 gene from raw sewage and animal samples. The raw sewage samples from the present study are indicated with a black dot followed by month, year and wastewater treatment plant of collection. Samples recovered from animals are indicated with white diamonds.
Viruses 12 00318 g003
Table
Table 1. Human adenovirus (HAdV) and rotavirus A (RVA) mean viral load during the study period. Quantification analysis of viral load (mean) detected by qPCR (HAdV) and RT-qPCR (RVA) in raw sewage samples collected between spring 2015 and spring 2016 in Catalonia.
Table 1. Human adenovirus (HAdV) and rotavirus A (RVA) mean viral load during the study period. Quantification analysis of viral load (mean) detected by qPCR (HAdV) and RT-qPCR (RVA) in raw sewage samples collected between spring 2015 and spring 2016 in Catalonia.
PeriodHAdV (CG/L)RVA (RT-PCR U/L)
Spring 20159.51 × 1043.96 × 104
Summer 20151.71 × 105-
Fall 20151.98 × 1051.73 × 106
Winter 20156.50 × 1051.41 × 107
Spring 20161.16 × 1063.30 × 108
Table
Table 2. Pairwise comparison analysis (p-values) between viral load and season using the Tukey and Kramer (Nemenyi) test.
Table 2. Pairwise comparison analysis (p-values) between viral load and season using the Tukey and Kramer (Nemenyi) test.
HAdVRVA
SummerFallWinterSummerFallWinter
Fall0.97538 0.0743
Winter0.000772.0 × 10−8 3.6 × 10−60.0011
Spring0.009411.6 × 10−50.936181.3 × 10−82.8 × 10−60.4376
Table
Table 3. VP4, VP7 and VP6 rotavirus A (RVA) genotypes during the study period. RVA genotypes circulating in Catalonia detected by qualitative RT-PCR and/or Sanger sequencing. The empty spaces in the table indicate that it was not possible to genotype the RVA strains from the given season.
Table 3. VP4, VP7 and VP6 rotavirus A (RVA) genotypes during the study period. RVA genotypes circulating in Catalonia detected by qualitative RT-PCR and/or Sanger sequencing. The empty spaces in the table indicate that it was not possible to genotype the RVA strains from the given season.
SeasonRVA Genotypes
VP4VP7VP6
Spring 2015P[4] I2
Summer 2015
Fall 2015P[4], P[8], P[9], P[10] I2
Winter 2015P[4], P[8], P[9] I1
Spring 2016P[4], P[8]G2, G3, G9, G12I1

Share and Cite

MDPI and ACS Style

Silva-Sales, M.; Martínez-Puchol, S.; Gonzales-Gustavson, E.; Hundesa, A.; Gironès, R. High Prevalence of Rotavirus A in Raw Sewage Samples from Northeast Spain. Viruses 2020, 12, 318. https://doi.org/10.3390/v12030318

AMA Style

Silva-Sales M, Martínez-Puchol S, Gonzales-Gustavson E, Hundesa A, Gironès R. High Prevalence of Rotavirus A in Raw Sewage Samples from Northeast Spain. Viruses. 2020; 12(3):318. https://doi.org/10.3390/v12030318

Chicago/Turabian Style

Silva-Sales, Marcelle, Sandra Martínez-Puchol, Eloy Gonzales-Gustavson, Ayalkibet Hundesa, and Rosina Gironès. 2020. "High Prevalence of Rotavirus A in Raw Sewage Samples from Northeast Spain" Viruses 12, no. 3: 318. https://doi.org/10.3390/v12030318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop