Detection of Common, Emerging and Uncommon VP4, and VP7 Human Group A Rotavirus Genotypes from Urban Sewage Samples in Uruguay

Rotavirus in Water Environments: A Systematic Review and Meta-Analysis

Background

Despite the adoption of rotavirus vaccines, sporadic outbreaks of the virus have been reported in many parts of the world. These outbreaks are facilitated by several factors including the ease of transmission of rotavirus through water environments. This systematic review aimed to determine the global prevalence of rotavirus in water environments.

Methodology

Comprehensive Boolean searches were conducted in PubMed, SCOPUS, and Web of Science. A total of 75 eligible studies were included in the study, from which data was extracted for both systematic review and meta-analysis. Extracted prevalence data was grouped according to six water categories: drinking water, untreated sewage, treated sewage, surface water, groundwater, and others. A single-group prevalence meta-analysis was conducted in RStudio version 4.3.3 subjecting the data to the random-effects model.

Results

The included studies were conducted in 32 countries that span 5 continents: Africa, Asia, Europe, North America, and South America. The pooled prevalence of rotavirus in water environments was 40.86%. Among the individual water environments, untreated sewage had the highest prevalence (68.27%), followed by treated sewage (53.07%), surface water (33.40%), groundwater (25.64%) and drinking water (9.46%). Continental stratification of the prevalence data was as follows: Africa (51.75%), Asia (32.48%), Europe (55.90%), North America (41.80%), and South America (28.51%).

Conclusion

There is a high prevalence of rotavirus in water environments, especially in untreated sewage, and in Europe. Further research is needed to find more efficient methods that can effectively eliminate rotavirus to insignificant levels in water environments.

Bioinformatics-based screening tool identifies a wide variety of human and zoonotic viruses in Trujillo-Peru wastewater

Peru was one of the most affected countries during the COVID-19 pandemic. Moreover, multiple other viral diseases (enteric, respiratory, bloodborne, and vector-borne) are endemic and rising. According to Peru's Ministry of Health, various health facilities in the country were reallocated for the COVID-19 pandemic, thereby leading to reduced action to curb other diseases. Many viral diseases in the area are under-reported and not recognized. The One Health approach, in addition to clinical testing, incorporates environmental surveillance for detection of infectious disease outbreaks. The purpose of this work is to use a screening tool that is based on molecular methods, high throughput sequencing and bioinformatics analysis of wastewater samples to identify virus-related diseases circulating in Trujillo-Peru. To demonstrate the effectiveness of the tool, we collected nine untreated wastewater samples from the Covicorti wastewater utility in Trujillo-Peru on October 22, 2022. High throughput metagenomic sequencing followed by bioinformatic analysis was used to assess the viral diversity of the samples. Our results revealed the presence of sequences associated with multiple human and zoonotic viruses including Orthopoxvirus, Hepatovirus, Rhadinovirus, Parechovirus, Mamastrovirus, Enterovirus, Varicellovirus, Norovirus, Kobuvirus, Bocaparvovirus, Simplexvirus, Spumavirus, Orthohepevirus, Cardiovirus, Molliscipoxvirus, Salivirus, Parapoxvirus, Gammaretrovirus, Alphavirus, Lymphocryptovirus, Erythroparvovirus, Sapovirus, Cosavirus, Deltaretrovirus, Roseolovirus, Flavivirus, Betacoronavirus, Rubivirus, Lentivirus, Betapolyomavirus, Rotavirus, Hepacivirus, Alphacoronavirus, Mastadenovirus, Cytomegalovirus and Alphapapillomavirus. For confirmation purposes, we tested the samples for the presence of selective viruses belonging to the genera detected above. PCR based molecular methods confirmed the presence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), monkeypox virus (MPXV), noroviruses GI and GII (NoVGI and NoVGII), and rotavirus A (RoA) in our samples. Furthermore, publicly available clinical data for selected viruses confirm our findings. Wastewater or other environmental media surveillance, combined with bioinformatics methods, has the potential to serve as a systematic screening tool for the identification of human or zoonotic viruses that may cause disease. The results of this method can guide further clinical surveillance efforts and allocation of resources. Incorporation of this bioinformatic-based screening tool by public health officials in Peru and other Latin American countries will help manage endemic and emerging diseases that could save human lives and resources.

Molecular Characterization of Gastroenteric Viruses in Wastewater from Cities in Uruguay

Group A Rotavirus, Human Astrovirus, and Norovirus (RVA, HAstV, and NoV) are recognized as the major causative agents of acute gastroenteritis in children and adults worldwide. The aim of this study was to determine the prevalence and molecular epidemiology of RVA, HAstV, and NoV in wastewater from three cities in Uruguay. Thirty-six samples from Bella Unión, Salto, and Fray Bentos cities were analyzed using quantitative and qualitative PCR. RVA was the most frequently detected virus (50%), followed by HAstV (39%), NoV GII (36%), and NoV GI (25%). RVA strains were characterized as P[8] and G3 based on the VP4 and VP7 genes, respectively. Among NoV-positive samples, genotypes GI.2, GI.3, GI.5, GI.6, GI.7, GII.2, GII.6, and GII.4 were detected, and only one HAstV genotype (MLB1) was found. Our wastewater-based epidemiological approach provides a snapshot of the overall genetic diversity of these viruses in three cities of the Uruguay River basin during 2017-2018. These findings reinforce the importance of this environmental surveillance tool for monitoring epidemiological trends of enteric viruses circulating in the population, which can be used to guide public health intervention.

Full genotype characterization of Brazilian canine G3P[3] strains during a 10-year survey (2012-2021) of rotavirus infection in domestic dogs and cats

There is a dearth of information on the molecular epidemiology of rotaviruses in pets in Brazil. The aim of this study was to monitor rotavirus infections in household dogs and cats, determine full-genotype constellations, and obtain data on evolutionary relationships. Between 2012 and 2021, 600 fecal samples from dogs and cats (516 and 84, respectively) were collected at small animal clinics in São Paulo state, Brazil. Rotavirus screening was conducted using ELISA, PAGE, RT-PCR, sequencing, and phylogenetic analysis. Rotavirus type A (RVA) was detected in 0.5% (3/600) of the animals. No non-RVA types were detected. The three canine RVA strains were found to have a novel genetic constellation, G3-P[3] -I2-R3-C2-M3-A9-N2-T3-E3-H6, which has never been reported in dogs. As expected, all of the viral genes, except those encoding NSP2 and VP7, were closely related to the corresponding genes from canine, feline, and canine-like-human RVA strains. A novel N2 (NSP2) lineage was identified, grouping together Brazilian canine, human, rat and bovine strains, suggesting that genetic reassortment had occurred. Uruguayan G3 strains obtained from sewage contained VP7 genes that were phylogenetically close to those of the Brazilian canine strains, which suggests that these strains are widely distributed in pet populations in South American countries. For the NSP2 (I2), NSP3 (T3), NSP4 (E3), NSP5 (H6), VP1 (R3), VP3 (M3), and VP6 (I2) segments, phylogenetic analysis revealed possibly new lineages. The epidemiological and genetic data presented here point out the necessity for collaborative efforts to implement the One Health strategy in the field of RVA research and to provide an updated understanding of RVA strains circulating canines in Brazil.

Temporal variations of human and animal Rotavirus A genotypes in surface water used for drinking water production

Rotavirus is a major cause of gastroenteritis among infants and children. In this study, nested PCR assays were developed to amplify partial regions of the VP7, VP4, and VP6 genes of Rotavirus A (RVA) for amplicon-based Illumina MiSeq sequencing to investigate RVA genotypes in environmental water samples. Eight sets of inner primers were first designed and screened for use in the nested PCR assays, and four sets of them could produce amplicons. Six sets of outer primers were then designed and combined with the four sets of inner primers that worked. The assays were evaluated for sensitivity using raw water samples collected from one drinking water treatment plant between April 2019 and March 2020 (Sample Set 1; N = 12) and seven DWTPs between 2018 and 2020 (Sample Set 2; N = 18). In total, 43 amplicons from Set 1 were sequenced and diverse sequences from human, porcine, bovine, equine, and feline RVA were observed. Human G8, G3, and G2 genotypes were obtained, with G8 predominant (relative abundance, 36–87%) in samples taken during the rotavirus epidemic season between April and June. Porcine G5, G11, and G4, and bovine G10 and G6 genotypes were also detected. VP4 sequence analysis revealed that the human P[8] genotype was present throughout the year, whereas P[4] and P[9] were present only in the epidemic season. The vaccine strains P[5] and P[8] (RotaTeq®) were also detected. Our approach enables the identification of prevalent human and animal RVA genotypes and their host species that potentially caused fecal contamination in water sources.

Epidemiological significance of the occurrence and persistence of rotaviruses in water and sewage: a critical review and proposal for routine microbiological monitoring

Globally, waterborne gastroenteritis attributable to rotaviruses is on the increase due to the rapid increase in population growth, poor socioeconomic conditions, and drastic changes in climatic conditions. The burden of diarrhea is quite alarming in developing nations where the majority of the populations still rely on untreated surface water that is usually polluted for their immediate water needs. Humans and animals of all ages are affected by rotaviruses. In humans, the preponderance of cases occurs in children under 5 years. Global efforts in advancing water/wastewater treatment technologies have not yet realized the objective of complete viral removal from wastewater. Most times, surface waters are impacted heavily by inadequately treated wastewater run-offs thereby exposing people or animals to preventable health risks. The relative stability of rotaviruses in aquatic matrices during wastewater treatment, poor correlation of bacteriological indicators with the presence of rotaviruses, and their infectiousness at a low dose informed the proposal for inclusion in the routine microbiological water screening panel. Environmental monitoring data have been shown to provide early warnings that can complement clinical data used to monitor the impact of current rotavirus vaccination in a community. This review was therefore undertaken to critically appraise rotavirus excretion and emission pathways, and the existence, viability and persistence in the receiving aquatic milieu. The efficiency of the current wastewater treatment modality for rotavirus removal, correlation of the current bacteriological water quality assessment strategy, public health risks and current laboratory methods for an epidemiological study were also discussed.

Detection and evaluation of rotavirus surveillance methods as viral indicator in the aquatic environments

Group A rotaviruses (RVAs) have been introduced as the most important causative agents of acute gastroenteritis in the young children. One of every 260 children born globally will die due to rotavirus (RV) before 5 years old. The RV is widely known as a viral indicator for health (fecal contamination) because this pathogen has a high treatment resistance nature, which has been listed as a relevant waterborne pathogen by the World Health Organization (WHO). Therefore, monitoring of environmental is important, and RV is one of the best-known indicators for monitoring. It has been proved that common standards for microbiological water quality do not guarantee the absence of viruses. On the other hand, in order to recover and determine RV quantity within water, standard methods are scarce. Therefore, dependable prediction of RV quantities in water sample is crucial to be able to improve supervision efficiency of the treatment procedure, precise quantitative evaluation of the microbial risks as well as microbiological water safety. Hence, this study aimed to introduce approaches to detecting and controlling RV in environmental waters, and discussed the challenges faced to enable a clear perception on the ubiquity of the RV within different types of water across the world.

Quantification and Trends of Rotavirus and Enterovirus in Untreated Sewage Using Reverse Transcription Droplet Digital PCR

The quantification and trends in concentrations for naturally occurring rotaviruses (RV) and enteroviruses (EV) in untreated sewage in various wastewater systems have not often been compared. There is now greater interest in monitoring the infections in the community including live vaccine efficacy by evaluating untreated sewage. The goals of this study were to 1) survey the concentrations of naturally occurring RV and EV in untreated sewage using a reverse transcription-droplet digital polymerase chain reaction (RT-ddPCR) and 2) investigate the use of a new adsorption elution (bag-mediated filtration system (BMFS) using ViroCap filters) against more traditional polyethylene glycol (PEG) precipitation for virus concentration. Sewage samples were collected from lagoons in Kenya and Michigan (MI), the United States (USA) and from wastewater treatment plants (WWTPs) in the USA. RVs were detected at geometric mean concentrations in various locations, California (CA) 1.31 × 10⁵ genome copies/L (gc/L), Kenya (KE) 2.71 × 10⁴ gc/L and Virginia (VA) 1.48 × 10⁵ gc/L, and EVs geometric means were 3.72 × 10⁶ gc/L (CA), 1.18 × 10⁴ gc/L (Kenya), and 6.18 × 10³ gc/L (VA). The mean RV concentrations using BMFS-ViroCap in split samples compared to PEG precipitation methods demonstrated that the levels were only 9% (#s BMFS/PEG) in the Michigan lagoons which was significantly different (p < 0.01). This suggests that RV concentrations in Kenya are around 1.69 × 10⁶ gc/L. Overall, there was no difference in concentrations for the other sampling locations across the methods of virus recovery (i.e., PEG precipitation and HA filters) using one-way ANOVA (F = 1.7, p = 0.2739) or Tukey-Kramer pairwise comparisons (p > 0.05). This study provides useful data on RV and EV concentrations in untreated sewage in Kenya and the USA. It also highlights on the usefulness of the RT-ddPCR for absolute quantification of RV and EV in sewage samples. The BMFS using ViroCap filters while less efficient compared to the more traditional PEG precipitation method was able to recover RVs and EVs in untreated sewage and may be useful in poor resource settings while underestimating viruses by 1 to 1.5 logs.

A Real-Time RT-PCR Assay for Genotyping of Rotavirus

Background

Human rotavirus (HRV) is the causative agent of severe gastroenteritis in children and responsible for two million hospitalizations and more than a half-million deaths annually. Sequence characteristics of the gene segments encoding the VP7 and VP4 proteins are used for the genotype classification of rotavirus. A wide variety of molecular methods are available, mainly based on reverse transcription PCR for rapid, specific and sensitive genotyping of rotaviruses. This study describes an alternative real-time PCR assay for genotyping of rotavirus.

Methods

The samples of stools studied in this research have been collected from patients referred to Children's Medical Centers, Tehran, Iran. Rotavirus detection and genotyping were performed using the RT-PCR and semi-nested RT-PCR, respectively. Samples were then genotyped with a new real-time PCR.

Results

The real-time PCR was able to genotype all positive samples with a mean Ct of 28.2. Besides, a concordance rate of 100% was detected between real-time PCR and semi-nested RT-PCR.

Conclusion

In this study, the genotyping of rotavirus with real-time PCR showed that this method can provide several favorable features, including high sensitivity and specificity, and a wide dynamic range for rotavirus genotyping.

Detection, Quantification, and Microbial Risk Assessment of Group A Rotavirus in Rivers from Uruguay

The aim of this study was to detect, quantify, and assess the risk of infection and illness for Group A Rotavirus (RVA) in the watersheds of the Santa Lucia and Uruguay rivers in Uruguay. Monthly sampling was carried out for one year in six sites in the watershed of the Santa Lucía River and four in the Uruguay River. All the collection sites are used for recreational activities. Viral concentration was performed with the adsorption-elution method, and detection and quantification of RVA was carried out by TaqMan quantitative PCR (qPCR). Quantitative microbial risk assessment was applied to estimate the daily and annual risk of RVA infection, as well as the daily risk of illness considering direct exposure through recreational activity. RVA was detected in 42% (20/48) of the analyzed samples in the Uruguay River and 40% (29/72) in the Santa Lucía River. The virus was present in all the analyzed points in both watersheds. A pattern of seasonality, characterized by a higher detection frequency of the virus during coldest month of the year, was observed in both basins. The mean concentration for RVA was 1.3 × 10⁵ genomic copies/L. The microbiological risk assessment shows that Santa Lucía watershed presented the highest daily risk of infection (6.41E-01) and illness (3.20E-01) estimated for the point downstream of Florida City; meanwhile for Uruguay River, the highest probabilities of infection (6.82E-01) and illness (3.41E-01) were estimated for the collection site for drinking water intake in Salto city. These results suggest that RVA contamination of these important rivers negatively impact on their microbiological quality since they are used for recreation and drinking water intake, demonstrating that the disposal of waste from cities located in their riverside confers a constant threat of infection for the general population, especially for children.

Molecular detection and genotyping of group A rotavirus in two wastewater treatment plants, Iran

In different countries especially developing countries, treatment of urban wastewater might be ineffective removal of pathogens such as group A rotavirus. The objective of this study is to evaluate the efficiency of rotavirus removal in two wastewater treatment plants (WWTPs) in Isfahan, Iran. To meet the study objectives, 96 sewage samples from influent (n = 48) and final effluents (n = 48) were collected by grab sampling. Two different concentration methods included pellet and two-phase used for concentrating sewage samples. The presence of rotavirus antigens in all concentrated sewage samples was screened by enzyme-linked immunosorbent method. To analyze the study samples, real-time PCR technique with SYBR Green I fluorescent dye and nested multiplex PCR for rotavirus genotyping were utilized respectively. Result indicated positive rotavirus percentage in two methods of ELISA and real-time PCR was equal to 61.45% (59 cases) and 44.79% (43 cases). In addition, analyzing seasonal distribution of rotavirus shows different distributions as below: in spring (18.64%), summer (20.33%), autumn (35.60%), and winter (25.42%). Finally, rotaviruses illustrate significantly higher frequency in cold seasons. G10 and G1 types are considered the most, among common genotypes which identified in 11 (25.58%) and 5 (11.62%), out of the 43 positive samples in WWTPs, followed by non-typeable genotypes (13.95%) and mix genotypes (11.62%); and different genotypes including G2, G3, G4, G8, G9, and G12 were equal to 2.33, 9.30, 9.30, 2.33, 7, and 7% in the WWTPs, respectively. Such high prevalence underlines the significance of environmental surveillance. Also, to eliminate potential pathogens especially enteric viruses from sewage, the improvement of treatment systems is essential.

Detection of Norovirus and Rotavirus Present in Suspended and Dissolved Forms in Drinking Water Sources

We investigated the present forms of genogroup II norovirus and group A rotavirus in surface water used for drinking water production. River water samples (N = 15) collected at a drinking water treatment plant (DWTP) monthly from June 2017 to August 2018 were fractioned by filtration through 10- and 0.45-μm-pore-size membranes, and viruses present in suspended and dissolved forms were quantitatively detected. Norovirus GII was present in > 10-μm- and 0.45-10-μm-suspended and dissolved forms with detection rates of 33%, 60%, and 87%, respectively. Rotavirus A was detected more frequently than norovirus GII in each form (> 10 μm suspended, 73%; 0.45-10 μm suspended, 93%; dissolved, 100%). We also analyzed surface water samples from 21 DWTPs all over Japan in non-epidemic and epidemic seasons of gastroenteritis. Norovirus GII was detected in 48% and 81% of samples with the concentrations of up to 4.1 and 5.3 log₁₀ copies/L in dissolved form in non-epidemic and epidemic seasons, respectively, and GII.4 Sydney 2012 was predominant genotype followed by GII.2. Rotavirus A was detected in 95% and 86% of samples with the maximum concentrations of 5.5 and 6.3 log₁₀ copies/L in dissolved form in respective seasons. Concentration of norovirus GII was similar in 0.45-10-μm suspended and dissolved forms, while there was a significant difference for rotavirus A (P < 0.01, pared t test), indicating that rotavirus A was less associated with suspended solids in the surface water samples compared to norovirus GII. Our observations provide important implications for understanding of viral behavior in environmental waters.

Hepatitis A Virus, Hepatitis E Virus, and Rotavirus in Foods of Animal Origin Traded at the Borders of Brazil, Argentina, and Uruguay

The aim of this study was to investigate hepatitis A virus (HAV), hepatitis E (HEV), and rotavirus (RV) in fresh and processed meat traded on the border of Brazil with Argentina and Uruguay. In total, 159 samples of raw and processed foods of animal origin were collected in Paso de los Libres, Argentina (n = 53 raw meat, n = 24 processed meat) and Rivera, Uruguay (n = 55 raw meat, n = 18 processed meat), or were seized by the Brazilian International Agricultural Surveillance System-VIGIAGRO (Brazil-Argentina border) (n = 8 raw meat, n = 1 bush meat). All samples were tested for the presence of HAV, HEV, and RV genomes. HAV genes were detected in 18.23% of samples and RV genes in 23.89%. No HEV-positive samples were detected. HAV was also detected in two of the VIGIAGRO samples. Processed meats from Argentina and Uruguay had a higher rate of HAV and RV than raw meat (P > 0.05). The median HAV in the Argentinian and Uruguayan samples was 6.9 × 10⁴ and 3.5 × 10³ copies/g, respectively. The presence of RV viral genes in raw meats from Argentina was significant, and this was not observed in processed meats. The presence of HAV and RV genes in a significant portion of products from Argentina and Uruguay is a potential source of human infection. This also indicates precarious conditions of acquisition, processing, and manipulation, which could be improved by improved regulation of food across borders.

An Environmental Surveillance in Uruguay Reveals the Presence of Highly Divergent Types of Human Enterovirus Species C and a High Frequency of Species A and B Types

Information about Human Enterovirus circulation in Uruguay is scarce. The aim of this study was to generate the first description about their circulation in the country through the study of sewage samples collected before and after the switch from Oral Poliovirus Vaccine to Inactivated Poliovirus Vaccine. Viruses were concentrated by an adsorption-elution to a negatively charged membrane, and real-time quantitative PCR and qualitative PCR methods were used to detect, quantify, and characterize enteroviruses. Positive samples were inoculated in RD cells and two passages were performed. Additionally, RD+ samples were subsequently passed onto L20B cells. Human Enteroviruses were detected in 67.6% of the samples, with concentrations between 4.9 and 6.6 Log₁₀ genomic copies per liter. 10% of positive samples replicated in RD cells, of which none in L20B cells. Molecular characterization of Human Enterovirus strains directly detected from sewage sample concentrates allowed the identification of highly divergent members of species C such as Enterovirus C99 and Coxsackievirus A13, as well as the frequent detection of species A and B members (particularly Coxsackievirus A16 and Echovirus 6, respectively). Other detected types were Coxsackievirus A2, A22, B1, B5, Echovirus 5, and 9. The characterization of viruses isolated in cell culture revealed the presence of Echovirus 6 and Coxsackievirus B3. Despite the absence of poliovirus, a wide circulation of different enterovirus types was evidenced in Uruguayan sewage samples, highlighting that the local populations are exposed to different kinds of diseases originated by several human enterovirus.

Species A Rotavirus (RVA) Isolated from Sewage in Nigeria, 2014: Close Genetic Relatedness of Partial G, P, and NSP4 Gene Sequences Encoding G1 with Cogent Genes of Other Asian and African Rotaviruses

Rotavirus has been identified as a major cause of gastroenteritis in Nigeria. There is limited information on the intragenotype diversity of Nigerian rotavirus isolates. We therefore investigated the molecular characteristics of some rotavirus gene sequences detected in sewage from Nigeria. Seven sewage samples, out of a total of 68, tested positive for rotavirus RNA (10.3%). Genotype G1P[4] was the most common genotype (5 isolates) and one isolate for genotypes G1P[8] and G3P[6]. Phylogenetic analysis of the partial VP7 gene of 3 G1P[4] isolates analyzed identified them as genotype G1 Lineage 2 along with Chinese strains with 99.1% to 100% amino acid similarity. Amino acid substitutions D-97→E and S-147→D/N were observed within the 7-1a and 7-2 domains of VP7 gene among the study G1P4 isolates in reference to vaccine strain RotaTeq®. Phylogenetic analysis of the G3P[6] study isolate identified it as genotype G3 Lineage 3, forming a monophyletic cluster with 100% bootstrap value with other West African strains G3 isolates. Phylogenetic analysis of GIP[4] VP4 genes identified them as P4 Lineage 5, while 3 NSP4 gene sequences belonged to genotype E1, while 1 belonged to E2. The results from this study represent phylogenetic analysis of partial gene sequences of environmental group A rotavirus (RVA) isolates from Nigeria.

Human Bocavirus: Detection, Quantification and Molecular Characterization in Sewage and Surface Waters in Uruguay

Human bocavirus (HBoV) infections are related to respiratory and gastroenteric diseases. The aim of this study was to investigate the presence of HBoV in both sewage and surface waters in Uruguay. Sixty-eight sewage samples from the cities of Salto, Paysandú, Bella Unión, Fray Bentos, Treinta y Tres and Melo and 36 surface water samples from the cities of Salto, Florida and Santa Lucía were studied. HBoV was screened by multiplex qPCR for the detection of the four subtypes, followed by monoplex qPCRs for the independent quantification of each subtype. A qualitative PCR followed by DNA sequencing and phylogenetic analysis was carried out for molecular characterization of HBoV strains. HBoV was present in a high frequency (69%) in sewage and only one positive sample (3%) was found in surface water. Concerning sewage samples, HBoV1 was detected in 11 (23%) out of the 47 positives samples, with a mean concentration of 8.2 × 10⁴ genomic copies/Liter (gc/L), HBoV3 was detected in 35 (74%) of the positive samples with a mean concentration of 4.1 × 10⁶ gc/L and subtypes 2 and/or 4 were detected in 39 (83%) of the positive samples with a mean concentration of 7.8 × 10⁶ gc/L. After the phylogenetic analysis performed by a Bayesian approach, the four HBoV subtypes were confirmed. This is the first study determining a high frequency of HBoV and the presence of the four HBoV subtypes in aquatic matrices in Latin America, mainly in sewage. Although HBoV was scarcely detected in surface water, a waterborne transmission is likely to occur if people enter in contact with polluted surface waters for recreational activities such as fishing or swimming since an elevated frequency of HBoV was detected in raw sewage which is usually directly discharged into surface waters.