Comparison of Different Reverse Transcriptase-Polymerase Chain Reaction-Based Methods for Wastewater Surveillance of SARS-CoV-2: Exploratory Study.

BACKGROUND: Many countries have applied the wastewater surveillance of the COVID-19 pandemic to their national public health monitoring measures. The most used methods for detecting SARS-CoV-2 in wastewater are quantitative reverse transcriptase-polymerase chain reaction (RT-qPCR) and reverse transcriptase-droplet digital polymerase chain reaction (RT-ddPCR). Previous comparison studies have produced conflicting results, thus more research on the subject is required. OBJECTIVE: This study aims to compare RT-qPCR and RT-ddPCR for detecting SARS-CoV-2 in wastewater. It also aimed to investigate the effect of changes in the analytical pipeline, including the RNA extraction kit, RT-PCR kit, and target gene assay, on the results. Another aim was to find a detection method for low-resource settings. METHODS: We compared 2 RT-qPCR kits, TaqMan RT-qPCR and QuantiTect RT-qPCR, and RT-ddPCR based on sensitivity, positivity rates, variability, and correlation of SARS-CoV-2 gene copy numbers in wastewater to the incidence of COVID-19. Furthermore, we compared 2 RNA extraction methods, column- and magnetic-bead-based. In addition, we assessed 2 target gene assays for RT-qPCR, N1 and N2, and 2 target gene assays for ddPCR N1 and E. Reverse transcription strand invasion-based amplification (RT-SIBA) was used to detect SARS-CoV-2 from wastewater qualitatively. RESULTS: Our results indicated that the most sensitive method to detect SARS-CoV-2 in wastewater was RT-ddPCR. It had the highest positivity rate (26/30), and its limit of detection was the lowest (0.06 gene copies/µL). However, we obtained the best correlation between COVID-19 incidence and SARS-CoV-2 gene copy number in wastewater using TaqMan RT-qPCR (correlation coefficient [CC]=0.697, P<.001). We found a significant difference in sensitivity between the TaqMan RT-qPCR kit and the QuantiTect RT-qPCR kit, the first having a significantly lower limit of detection and a higher positivity rate than the latter. Furthermore, the N1 target gene assay was the most sensitive for both RT-qPCR kits, while no significant difference was found between the gene targets using RT-ddPCR. In addition, the use of different RNA extraction kits affected the result when the TaqMan RT-qPCR kit was used. RT-SIBA was able to detect SARS-CoV-2 RNA in wastewater. CONCLUSIONS: As our study, as well as most of the previous studies, has shown RT-ddPCR to be more sensitive than RT-qPCR, its use in the wastewater surveillance of SARS-CoV-2 should be considered, especially if the amount of SARS-CoV-2 circulating in the population was low. All the analysis steps must be optimized for wastewater surveillance as our study showed that all the analysis steps including the compatibility of the RNA extraction, the RT-PCR kit, and the target gene assay influence the results. In addition, our study showed that RT-SIBA could be used to detect SARS-CoV-2 in wastewater if a qualitative result is sufficient.

Wastewater-based surveillance is an efficient monitoring tool for tracking influenza A in the community.

Around the world, influenza A virus has caused severe pandemics, and the risk of future pandemics remains high. Currently, influenza A virus surveillance is based on the clinical diagnosis and reporting of disease cases. In this study, we apply wastewater-based surveillance to monitor the amount of the influenza A virus RNA at the population level. We report the influenza A virus RNA levels in 10 wastewater treatment plant catchment areas covering 40 % of the Finnish population. Altogether, 251 monthly composite influent wastewater samples (collected between February 2021 and February 2023) were analysed from supernatant fraction using influenza A virus specific RT-qPCR method. During the study period, an influenza A virus epidemic occurred in three waves in Finland. This study shows that the influenza A virus RNA can be detected from the supernatant fraction of 24 h composite influent wastewater samples. The influenza A virus RNA gene copy number in wastewater correlated with the number of confirmed disease cases in the Finnish National Infectious Diseases Register. The median Kendall's τ correlation strength was 0.636 (min= 0.486 and max=0.804) and it was statistically significant in all 10 WTTPs. Wastewater-based surveillance of the influenza A virus RNA is an independent from individual testing method and cost-efficiently reflects the circulation of the virus in the entire population. Thus, wastewater monitoring complements the available, but often too sparse, information from individual testing and improves health care and public health preparedness for influenza A virus pandemics.

Developing wastewater-based surveillance schemes for multiple pathogens: The WastPan project in Finland.

Wastewater comprises multiple pathogens and offers a potential for wastewater-based surveillance (WBS) to track the prevalence of communicable diseases. The Finnish WastPan project aimed to establish wastewater-based pandemic preparedness for multiple pathogens (viruses, bacteria, parasites, fungi), including antimicrobial resistance (AMR). This article outlines WastPan's experiences in this project, including the criteria for target selection, sampling locations, frequency, analysis methods, and results communication. Target selection relied on epidemiological and microbiological evidence and practical feasibility. Within the WastPan framework, wastewater samples were collected between 2021 and 2023 from 10 wastewater treatment plants (WWTPs) covering 40 % of Finland's population. WWTP selection was validated for reported cases of Extended Spectrum Beta-lactamase-producing bacterial pathogens (Escherichia coli and Klebsiella pneumoniae) from the National Infectious Disease Register. The workflow included 24-h composite influent samples, with one fraction for culture-based analysis (bacteria and fungi) and the rest of the sample was reserved for molecular analysis (viruses, bacteria, antibiotic resistance genes, and parasites). The reproducibility of the monitoring workflow was assessed for SARS-CoV-2 through inter-laboratory comparisons using the N2 and N1 assays. Identical protocols were applied to same-day samples, yielding similar positivity trends in the two laboratories, but the N2 assay achieved a significantly higher detection rate (Laboratory 1: 91.5 %; Laboratory 2: 87.4 %) than the N1 assay (76.6 %) monitored only in Laboratory 2 (McNemar, p < 0.001 Lab 1, = 0.006 Lab 2). This result indicates that the selection of monitoring primers and assays may impact monitoring sensitivity in WBS. Overall, the current study recommends that the selection of sampling frequencies and population coverage of the monitoring should be based on pathogen-specific epidemiological characteristics. For example, pathogens that are stable over time may need less frequent annual sampling, while those that are occurring across regions may require reduced sample coverage. Here, WastPan successfully piloted WBS for monitoring multiple pathogens, highlighting the significance of one-litre community composite wastewater samples for assessing community health. The infrastructure established for COVID-19 WBS is valuable for monitoring various pathogens. The prioritization of the monitoring targets optimizes resource utilization. In the future legislative support in target selection, coverage determination, and sustained funding for WBS is recomended.

Phenotypic Recovery of a Heterobasidion Isolate Infected by a Debilitation-Associated Virus Is Related to Altered Host Gene Expression and Reduced Virus Titer.

The fungal genus Heterobasidion includes forest pathogenic species hosting a diverse group of partitiviruses. They include the host debilitating Heterobasidion partitivirus 13 strain an1 (HetPV13-an1), which was originally observed in a slowly growing H. annosum strain 94233. In this study, a relatively fast-growing sector strain 94233-RC3 was isolated from a highly debilitated mycelial culture of 94233, and its gene expression and virus transcript quantities as well as the genomic sequence of HetPV13-an1 were examined. The sequence of HetPV13-an1 genome in 94233-RC3 was identical to that in the original 94233, and thus not the reason for the partial phenotypic recovery. According to RNA-seq analysis, the HetPV13-an1 infected 94233-RC3 transcribed eight genes differently from the partitivirus-free 94233-32D. Three of these genes were downregulated and five upregulated. The number of differentially expressed genes was considerably lower and the changes in their expression were small compared to those of the highly debilitated original strain 94233 with the exception of the most highly upregulated ones, and therefore viral effects on the host transcriptome correlated with the degree of the virus-caused debilitation. The amounts of RdRp and CP transcripts of HetPV13-an1 were considerably lower in 94233-RC3 and also in 94233 strain infected by a closely related mildly debilitating virus HetPV13-an2, suggesting that the virus titer would have a role in determining the effect of HetPV13 viruses on their hosts.

Multiple virus infections on Heterobasidion sp.

Heterobasidion viruses have previously been shown to affect each other's transmission and phenotypic effects on their hosts in a complex way. In this work, Heterobasidion parviporum strains hosting five coinfecting viruses simultaneously were constructed and used as donors in transmission experiments. They showed that viruses move more frequently between the mycelia of the same species than between the mycelia of H. parviporum and Heterobasidionannosum. One of the strains was used to show that coinfection of five viruses is relatively unstable in a natural environment and analyses of the growth rate and competitive ability of Heterobasidion strains hosting various virus combinations revealed that viral effects are not additive. The results also supported the view that the transmission of the promising virocontrol agent HetPV13-an1 may be enhanced by coinfecting viruses in the donor mycelium. However, its detrimental effects may be blocked by the presence of other viruses in the same mycelium. REPOSITORIES: GenBank accession number MN058080.

Distribution of Viruses Inhabiting Heterobasidion annosum in a Pine-Dominated Forest Plot in Southern Finland.

We investigated the diversity and spatial distribution of viruses infecting strains of the root rot fungus Heterobasidion annosum collected from pine stumps at a heavily infected forest site. Four different partitiviruses were detected in 14 H. annosum isolates at the study site, constituting approximately 29% of all Heterobasidion isolates investigated (N = 48). Two of the viruses detected were new partitiviruses designated here as Heterobasidion partitivirus 16 (HetPV16) and HetPV20, and two were previously known partitiviruses: HetPV7 and HetPV13. The two new partitiviruses found, HetPV16-an1 and HetPV20-an1, shared ~70% RdRp nucleotide sequence identity with the alphapartitivirus Rosellinia necatrix partitivirus 2, and less than 40% identity with known viruses of Heterobasidion spp. HetPV7-an1 was closely similar to HetPV7-pa1 isolated earlier from Heterobasidion parviporum, supporting the view of conspecific virus pools in different Heterobasidion species. Three fungal isolates were found to be co-infected with two different partitivirus strains (HetPV7-an1 and HetPV13-an2 or HetPV16-an1 and HetPV20-an1). Different isolates representing each host clone had variable virus compositions, and virus strains occurring in more than one host clone showed minor sequence variations between clones.

Heterobasidion Partitivirus 13 Mediates Severe Growth Debilitation and Major Alterations in the Gene Expression of a Fungal Forest Pathogen.

The fungal genus Heterobasidion includes some of the most devastating conifer pathogens in the boreal forest region. In this study, we showed that the alphapartitivirus Heterobasidion partitivirus 13 from Heterobasidion annosum (HetPV13-an1) is the main causal agent of severe phenotypic debilitation in the host fungus. Based on RNA sequencing using isogenic virus-infected and cured fungal strains, HetPV13-an1 affected the transcription of 683 genes, of which 60% were downregulated and 40% upregulated. Alterations observed in carbohydrate and amino acid metabolism suggest that the virus causes a state of starvation, which is compensated for by alternative synthesis routes. We used dual cultures to transmit HetPV13-an1 into new strains of H. annosum and Heterobasidion parviporum The three strains of H. parviporum that acquired the virus showed noticeable growth reduction on rich culturing medium, while only two of six H. annosum isolates tested showed significant debilitation. Based on reverse transcription-quantitative PCR (RT-qPCR) analysis, the response toward HetPV13-an1 infection was somewhat different in H. annosum and H. parviporum We assessed the effects of HetPV13-an1 on the wood colonization efficacy of H. parviporum in a field experiment where 46 Norway spruce trees were inoculated with isogenic strains with or without the virus. The virus-infected H. parviporum strain showed considerably less growth within living trees than the isolate without HetPV13-an1, indicating that the virus also causes growth debilitation in natural substrates.IMPORTANCE A biocontrol method restricting the spread of Heterobasidion species would be highly beneficial to forestry, as these fungi are difficult to eradicate from diseased forest stands and cause approximate annual losses of €800 million in Europe. We used virus curing and reintroduction experiments and RNA sequencing to show that the alphapartitivirus HetPV13-an1 affects many basic cellular functions of the white rot wood decay fungus Heterobasidion annosum, which results in aberrant hyphal morphology and a low growth rate. Dual fungal cultures were used to introduce HetPV13-an1 into a new host species, Heterobasidion parviporum, and field experiments confirmed the capability of the virus to reduce the growth of H. parviporum in living spruce wood. Taken together, our results suggest that HetPV13-an1 shows potential for the development of a future biocontrol agent against Heterobasidion fungi.

Population structure of a novel putative mycovirus infecting the conifer root-rot fungus Heterobasidion annosum sensu lato.

We describe a novel putative mycovirus infecting the conifer root-rot fungus Heterobasidion annosum sensu lato. This virus, designated as Heterobasidion RNA virus 6 (HetRV6), is taxonomically distant from all previously known viruses of Heterobasidion species, but somewhat related to the Curvularia thermal tolerance virus and the Fusarium graminearum virus 4. Based on a population analysis including 35 virus strains from Heterobasidion abietinum, Heterobasidion parviporum, Heterobasidion annosum sensu stricto and Heterobasidion occidentale, HetRV6 showed a considerable degree of geographical and host-related differentiation. The North American and Eurasian virus populations were clearly separated. In Eurasia, we observed cases of discrepancy between virus and host taxonomy, suggesting interspecies virus transfer. HetRV6 was also successfully transmitted between the three European species H. abietinum, H. annosum and H. parviporum. Based on growth rate tests on agar plates and spruce stem pieces, HetRV6 seemed to be cryptic or slightly mutualistic to its host.

Effect of nutrient loading on bacterioplankton community composition in lake mesocosms.

Changes in bacterioplankton community composition were followed in mesocosms set up in the littoral of Lake Vesijärvi, southern Finland, over two summers. Increasing nitrogen and phosphorus concentrations in the mesocosms represented different trophic states, from mesotrophic to hypertrophic. In 1998, the mesocosms were in a turbid state with a high biomass of phytoplankton, whereas in 1999, macrophytes proliferated and a clear-water state prevailed. The bacterial communities in the mesocosms also developed differently, as shown by denaturing gradient gel electrophoresis profiling of partial 16S rRNA gene fragments and by nonmetric multidimensional scaling analysis. In 1998, nutrient treatments affected the diversity and clustering of bacterial communities strongly, but in 1999, the bacterial communities were less diversified and not clearly affected by treatments. Canonical correspondence analysis indicated that bacterioplankton communities in the mesocosms were influenced by environmental physicochemical variables linked to the increasing level of eutrophication. Nitrogen concentration correlated directly with the bacterioplankton composition. In addition, the high nutrient levels had indirect effects through changes in the biomass and composition of phyto- and zooplankton. Sequencing analysis showed that the dominant bacterial divisions remained the same, but the dominant phylotypes changed during the 2-year period. The occurrence of Verrucomicrobia correlated with more eutrophic conditions, whereas the occurrence of Actinobacteria correlated with less eutrophic conditions.