Detection of Human Sapoviruses in Sewage in China by Next Generation Sequencing

Propagating and banking genetically diverse human sapovirus strains using a human duodenal cell line: investigating antigenic differences between strains

There are four genogroups and 18 genotypes of human sapoviruses (HuSaVs) responsible for acute gastroenteritis. To comprehend their antigenic and virological differences, it is crucial to obtain viral stocks of the different strains. Previously, we utilized the human duodenum-derived cell line HuTu80, and glycocholate, a conjugated bile acid, to replicate and propagate GI.1, GI.2, and GII.3 HuSaVs (H. Takagi et al., Proc Natl Acad Sci U S A 117:32078-32085, 2020, https://10.1073/pnas.2007310117). First, we investigated the impact of HuTu80 passage number on HuSaV propagation. Second, we demonstrated that taurocholate improved the initial replication success rate and viral RNA levels in fecal specimens relative to glycocholate. By propagating 15 HuSaV genotypes (GI.1-7, GII.1-5, -8, and GV.1-2) and accomplishing preparation of viral stocks containing 1.0 × 109 to 3.4 × 1011 viral genomic copies/mL, we found that all strains required bile acids for replication, with GII.4 showing strict requirements for taurocholate. The deduced VP1 sequences of the viruses during the scale-up of serial passaged virus cultures were either identical or differed by only two amino acids from the original sequences in feces. In addition, we purified virions from nine strains of different genotypes and used them as immunogens for antiserum production. Enzyme-linked immunosorbent assays (ELISAs) using rabbit and guinea pig antisera for each of the 15 strains of different genotypes revealed distinct antigenicity among the propagating viruses across genogroups and differences between genotypes. Acquisition of biobanked viral resources and determination of key culture conditions will be valuable to gain insights into the common mechanisms of HuSaV infection.

IMPORTANCE

The control of human sapovirus, which causes acute gastroenteritis in individuals of all ages, is challenging because of its association with outbreaks similar to those caused by human norovirus. The establishment of conditions for efficient viral propagation of various viral strains is essential for understanding the infection mechanism and identifying potential control methods. In this study, two critical factors for human sapovirus propagation in a conventional human duodenal cell line were identified, and 15 strains of different genotypes that differed at the genetic and antigenic levels were isolated and used to prepare virus stocks. The preparation of virus stocks has not been successful for noroviruses, which belong to the same family as sapoviruses. Securing virus stocks of multiple human sapovirus strains represents a significant advance toward establishing a reliable experimental system that does not depend on limited virus-positive fecal material.

Molecular epidemiology of enteroviruses from Guatemalan wastewater isolated from human lung fibroblasts

The Global Specialized Polio Laboratory at CDC supports the Global Poliovirus Laboratory Network with environmental surveillance (ES) to detect the presence of vaccine strain polioviruses, vaccine-derived polioviruses, and wild polioviruses in high-risk countries. Environmental sampling provides valuable supplementary information, particularly in areas with gaps in surveillance of acute flaccid paralysis (AFP) mainly in children less than 15 years. In collaboration with Guatemala's National Health Laboratory (Laboratorio Nacional de Salud Guatemala), monthly sewage collections allowed screening enterovirus (EV) presence without incurring additional costs for sample collection, transport, or concentration. Murine recombinant fibroblast L-cells (L20B) and human rhabdomyosarcoma (RD) cells are used for the isolation of polioviruses following a standard detection algorithm. Though non-polio-Enteroviruses (NPEV) can be isolated, the algorithm is optimized for the detection of polioviruses. To explore if other EV's are present in sewage not found through standard methods, five additional cell lines were piloted in a small-scale experiment, and next-generation sequencing (NGS) was used for the identification of any EV types. Human lung fibroblast cells (HLF) were selected based on their ability to isolate EV-A genus. Sewage concentrates collected between 2020-2021 were isolated in HLF cells and any cytopathic effect positive isolates used for NGS. A large variety of EVs, including echoviruses 1, 3, 6, 7, 11, 13, 18, 19, 25, 29; coxsackievirus A13, B2, and B5, EV-C99, EVB, and polioviruses (Sabin 1 and 3) were identified through genomic typing in NGS. When the EV genotypes were compared by phylogenetic analysis, it showed many EV's were genomically like viruses previously isolated from ES collected in Haiti. Enterovirus occurrence did not follow a seasonality, but more diverse EV types were found in ES collection sites with lower populations. Using the additional cell line in the existing poliovirus ES algorithm may add value by providing data about EV circulation, without additional sample collection or processing. Next-generation sequencing closed gaps in knowledge providing molecular epidemiological information on multiple EV types and full genome sequences of EVs present in wastewater in Guatemala.

Strong evolutionary constraints against amino acid changes in the P2 subdomain of sapovirus GI.1 capsid protein VP1

Sapovirus (SaV) is a nonenveloped RNA virus that causes acute gastroenteritis in humans. Although SaV is a clinically important pathogen in children, an effective vaccine is currently unavailable. The capsid protein VP1 of SaVs forms the outer shell of the virion and is highly diverse, as often seen in the virion-surface proteins of RNA viruses, creating an obstacle for vaccine development. We here report a unique phenomenon pertaining to the variation of SaV VP1. Phylogenetic and information entropy analyses using full-length VP1 sequences from a public database consistently showed that the amino acid sequences of the VP1 protein have been highly conserved over more than 40 years in the major epidemic genotype GI.1 but not in GI.2. Structural modeling showed that even the VP1 P2 subdomain, which is arranged on the outermost shell of the virion and presumably exposed to anti-SaV antibodies, remained highly homogeneous in GI.1 but not in GI.2. These results suggest strong evolutionary constraints against amino acid changes in the P2 subdomain of the SaV GI.1 capsid and illustrate a hitherto unappreciated mechanism, i.e., preservation of the VP1 P2 subdomain, involved in SaV survival. Our findings could have important implications for the development of an anti-SaV vaccine.

Environmental Surveillance of Human Astroviruses in Jinan City of China, 2020-2021

Human astroviruses (HAstVs) are a significant etiological agent of acute gastroenteritis in children. In order to investigate the circulation of HAstVs during the COVID-19 pandemic, a 2-year environmental surveillance was conducted in Jinan between 2020 and 2021. A total of 24 sewage samples were collected and concentrated. Real-time PCR indicated a positive rate of 83.3%, 79.2% (19/24), and 62.5% for classic, MLB, and VA types of HAstV in sewage samples, respectively, with genomic copies ranging from 6.4 × 103 to 3.7 × 107, 3.2 × 104 to 2.2 × 106, and 1.2 × 104 to 1.6 × 107 l-1. Next-generation sequencing (NGS) analysis on complete ORF2 amplicons from each sewage concentrate revealed the presence of 11 HAstV types, including HAstV-1, -2, -4, -5, MLB1, and VA1 to VA6, as well as non-human animal astroviruses. The most abundant HAstV types were HAstV-1, -4, and -5, which accounted for 70.3%, 12.6%, and 9.1% of total HAstV reads, respectively. Phylogenetic analysis revealed that the sequences obtained in this study were segregated into multiple transmission lineages, yet exhibited less genetic divergence among themselves than with foreign strains. These findings provide insight into the genotype diversity and genetic characterization of HAstVs during the COVID-19 pandemic, and highlight the effectiveness of utilizing NGS approaches to investigate sewage HAstVs.

The analysis of the genotype of Sapovirus outbreaks in Zhejiang Province

BACKGROUND

Sapovirus (SaV) infection is increasing globally. Concurrently, several SaV-outbreaks were observed in children of Zhejiang province, China, in recent years, In this study, the genotypes of Sapovirus from seven outbreaks in the Zhejiang province were analysed.

METHODS

A total of 105 faecal samples were collected from children aged between 4 and 17 years from the Zhejiang Provincial Center for Disease Control and Prevention between October 2021 and February 2023. Genotypes were processed using reverse transcription polymerase chain reaction and Sanger sequencing, while next-generation sequencing was used to generate a complete viral genome. Deduced amino acid sequences were analysed to detect VP1 gene mutations.

RESULTS

In total, 60 SaV-positive patients were detected at a 57.14% (60/105) positivity rate. Positive rates in the seven outbreaks were: 22.22% (2/9), 15.00% (3/20), 93.10% (27/29), 84.21% (16/19), 28.57% (2/7), 53.33% (8/15) and 33.33% (2/6), respectively. Four genotypes were identified in the seven outbreaks, of which, GI.1 accounted for 14.29% (1/7), GI.2 accounted for 14.29% (1/7), GI.6 and GII.5 accounted for 14.29% (1/7), and GI.6 accounted for 57.14% (4/7). All patients were children and outbreaks predominantly occurred in primary schools and during cold seasons. Additionally, the complete sequence from the GI.6 outbreak strain showed high homology (identity: 99.99%) with few common substitutions (Y300S, N302S and L8M) in VP1 protein.

CONCLUSIONS

SaV genotype diversity was observed in the seven outbreaks, with GI.6 being the main SaV genotype in Zhejiang province. It demonstrated high homology and may provide a platform for SaV prevention and control measures.

Characterization of a Human Sapovirus Genotype GII.3 Strain Generated by a Reverse Genetics System: VP2 Is a Minor Structural Protein of the Virion

We devised a reverse genetics system to generate an infectious human sapovirus (HuSaV) GII.3 virus. Capped/uncapped full-length RNAs derived from HuSaV GII.3 AK11 strain generated by in vitro transcription were used to transfect HuTu80 human duodenum carcinoma cells; infectious viruses were recovered from the capped RNA-transfected cells and passaged in the cells. Genome-wide analyses indicated no nucleotide sequence change in the virus genomes in the cell-culture supernatants recovered from the transfection or those from the subsequent infection. No virus growth was detected in the uncapped RNA-transfected cells, suggesting that the 5′-cap structure is essential for the virus’ generation and replication. Two types of virus particles were purified from the cell-culture supernatant. The complete particles were 39.2-nm-dia., at 1.350 g/cm3 density; the empty particles were 42.2-nm-dia. at 1.286 g/cm3. Two proteins (58-kDa p58 and 17-kDa p17) were detected from the purified particles; their molecular weight were similar to those of VP1 (~60-kDa) and VP2 (~16-kDa) of AK11 strain deduced from their amino acids (aa) sequences. Protein p58 interacted with HuSaV GII.3-VP1-specific antiserum, suggesting that p58 is HuSaV VP1. A total of 94 (57%) aa of p17 were identified by mass spectrometry; the sequences were identical to those of VP2, indicating that the p17 is the VP2 of AK11. Our new method produced infectious HuSaVs and demonstrated that VP2 is the minor protein of the virion, suggested to be involved in the HuSaV assembly.

Quantification of human enteric viruses as alternative indicators of fecal pollution to evaluate wastewater treatment processes

We investigated the potential use and quantification of human enteric viruses in municipal wastewater samples of Winnipeg (Manitoba, Canada) as alternative indicators of contamination and evaluated the processing stages of the wastewater treatment plant. During the fall 2019 and winter 2020 seasons, samples of raw sewage, activated sludge, effluents, and biosolids (sludge cake) were collected from the North End Sewage Treatment Plant (NESTP), which is the largest wastewater treatment plant in the City of Winnipeg. DNA (Adenovirus and crAssphage) and RNA enteric viruses (Pepper mild mottle virus, Norovirus genogroups GI and GII, Rotavirus Astrovirus, and Sapovirus) as well as the uidA gene found in Escherichia coli were targeted in the samples collected from the NESTP. Total nucleic acids from each wastewater treatment sample were extracted using a commercial spin-column kit. Enteric viruses were quantified in the extracted samples via quantitative PCR using TaqMan assays. Overall, the average gene copies assessed in the raw sewage were not significantly different (p-values ranged between 0.1023 and 0.9921) than the average gene copies assessed in the effluents for DNA and RNA viruses and uidA in terms of both volume and biomass. A significant reduction (p-value ≤ 0.0438) of Adenovirus and Noroviruses genogroups GI and GII was observed in activated sludge samples compared with those for raw sewage per volume. Higher GCNs of enteric viruses were observed in dewatered sludge samples compared to liquid samples in terms of volume (g of sample) and biomass (ng of nucleic acids). Enteric viruses found in gene copy numbers were at least one order of magnitude higher than the E. coli marker uidA, indicating that enteric viruses may survive the wastewater treatment process and viral-like particles are being released into the aquatic environment. Viruses such as Noroviruses genogroups GI and GII, and Rotavirus were detected during colder months. Our results suggest that Adenovirus, crAssphage, and Pepper mild mottle virus can be used confidently as complementary viral indicators of human fecal pollution.