Comparison of metagenomic and targeted methods for sequencing human pathogenic viruses from wastewater.

Wastewater-based epidemiology is a powerful tool for monitoring the emergence and spread of viral pathogens at the population scale. Typical polymerase chain reaction (PCR)-based methods of quantitative and genomic monitoring of viruses in wastewater provide high sensitivity and specificity. However, these methods are limited to the surveillance of target viruses in a single assay and require prior knowledge of the target genome(s). Metagenomic sequencing methods may represent a target-agnostic approach to viral wastewater monitoring, allowing for the detection of a broad range of target viruses, including potentially novel and emerging pathogens. In this study, targeted and untargeted metagenomic sequencing methods were compared with tiled-PCR sequencing for the detection and genotyping of viral pathogens in wastewater samples. Deep shotgun metagenomic sequencing was unable to generate sufficient genome coverage of human pathogenic viruses for robust genomic epidemiology, with samples dominated by bacteria. Hybrid-capture enrichment of shotgun libraries for respiratory viruses led to significant increases in genome coverage for a range of targets. Tiled-PCR sequencing led to further improvements in genome coverage compared to hybrid capture for severe acute respiratory syndrome coronavirus 2, enterovirus D68, norovirus GII, and human adenovirus F41 in wastewater samples. In conclusion, untargeted shotgun sequencing was unsuitable for genomic monitoring of the low virus concentrations in wastewater samples analyzed in this study. Hybrid-capture enrichment represented a viable method for simultaneous genomic epidemiology of a range of viral pathogens, while tiled-PCR sequencing provided the optimal genome coverage for individual viruses with the minimum sequencing depth. IMPORTANCE Most public health initiatives that monitor viruses in wastewater have utilized quantitative polymerase chain reaction (PCR) and whole genome PCR sequencing, mirroring techniques used for viral epidemiology in individuals. These techniques require prior knowledge of the target viral genome and are limited to monitoring individual or small groups of viruses. Metagenomic sequencing may offer an alternative strategy for monitoring a broad spectrum of viruses in wastewater, including novel and emerging pathogens. In this study, while amplicon sequencing gave high viral genome coverage, untargeted shotgun sequencing of total nucleic acid samples was unable to detect human pathogenic viruses with enough sensitivity for use in genomic epidemiology. Enrichment of shotgun libraries for respiratory viruses using hybrid-capture technology provided genotypic information on a range of viruses simultaneously, indicating strong potential for wastewater surveillance. This type of targeted metagenomics could be used for monitoring diverse targets, such as pathogens or antimicrobial resistance genes, in environmental samples.

Postpandemic rebound of adeno-associated virus type 2 (AAV2) infections temporally associated with an outbreak of unexplained severe acute hepatitis in children in the United Kingdom.

Over 1000 cases of unexplained severe acute hepatitis in children have been reported to date worldwide. An association with adeno-associated virus type 2 (AAV2) infection, a human parvovirus, prompted us to investigate the epidemiology of AAV in the United Kingdom. Three hundred pediatric respiratory samples collected before (April 03, 2009-April 03, 2013) and during (April 03, 2022) the COVID-19 pandemic were obtained. Wastewater samples were collected from 50 locations in London (August 2021-March 2022). Samples were tested for AAV using real-time polymerase chain reaction followed by sequencing. Selected adenovirus (AdV)-positive samples were also sequenced. The detection frequency of AAV2 was a sevenfold higher in 2022 samples compared with 2009-2013 samples (10% vs. 1.4%) and highest in AdV-positive samples compared with negatives (10/37, 27% vs. 5/94, 5.3%, respectively). AAV2-positive samples displayed high genetic diversity. AAV2 sequences were either very low or absent in wastewater collected in 2021 but increased in January 2022 and peaked in March 2022. AAV2 was detected in children in association with AdV of species C, with a highest frequency in 2022. Our findings are consistent with the expansion of the population of children unexposed to AAV2, leading to greater spread of the virus once distancing restrictions were lifted.

Wastewater monitoring for detection of public health markers during the COVID-19 pandemic: Near-source monitoring of schools in England over an academic year.

BACKGROUND: Schools are high-risk settings for infectious disease transmission. Wastewater monitoring for infectious diseases has been used to identify and mitigate outbreaks in many near-source settings during the COVID-19 pandemic, including universities and hospitals but less is known about the technology when applied for school health protection. This study aimed to implement a wastewater surveillance system to detect SARS-CoV-2 and other public health markers from wastewater in schools in England. METHODS: A total of 855 wastewater samples were collected from 16 schools (10 primary, 5 secondary and 1 post-16 and further education) over 10 months of school term time. Wastewater was analysed for SARS-CoV-2 genomic copies of N1 and E genes by RT-qPCR. A subset of wastewater samples was sent for genomic sequencing, enabling determination of the presence of SARS-CoV-2 and emergence of variant(s) contributing to COVID-19 infections within schools. In total, >280 microbial pathogens and >1200 AMR genes were screened using RT-qPCR and metagenomics to consider the utility of these additional targets to further inform on health threats within the schools. RESULTS: We report on wastewater-based surveillance for COVID-19 within English primary, secondary and further education schools over a full academic year (October 2020 to July 2021). The highest positivity rate (80.4%) was observed in the week commencing 30th November 2020 during the emergence of the Alpha variant, indicating most schools contained people who were shedding the virus. There was high SARS-CoV-2 amplicon concentration (up to 9.2x106 GC/L) detected over the summer term (8th June - 6th July 2021) during Delta variant prevalence. The summer increase of SARS-CoV-2 in school wastewater was reflected in age-specific clinical COVID-19 cases. Alpha variant and Delta variant were identified in the wastewater by sequencing of samples collected from December to March and June to July, respectively. Lead/lag analysis between SARS-CoV-2 concentrations in school and WWTP data sets show a maximum correlation between the two-time series when school data are lagged by two weeks. Furthermore, wastewater sample enrichment coupled with metagenomic sequencing and rapid informatics enabled the detection of other clinically relevant viral and bacterial pathogens and AMR. CONCLUSIONS: Passive wastewater monitoring surveillance in schools can identify cases of COVID-19. Samples can be sequenced to monitor for emerging and current variants of concern at the resolution of school catchments. Wastewater based monitoring for SARS-CoV-2 is a useful tool for SARS-CoV-2 passive surveillance and could be applied for case identification and containment, and mitigation in schools and other congregate settings with high risks of transmission. Wastewater monitoring enables public health authorities to develop targeted prevention and education programmes for hygiene measures within undertested communities across a broad range of use cases.

Evaluation of variant calling algorithms for wastewater-based epidemiology using mixed populations of SARS-CoV-2 variants in synthetic and wastewater samples.

Wastewater-based epidemiology has been used extensively throughout the COVID-19 (coronavirus disease 19) pandemic to detect and monitor the spread and prevalence of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and its variants. It has proven an excellent, complementary tool to clinical sequencing, supporting the insights gained and helping to make informed public-health decisions. Consequently, many groups globally have developed bioinformatics pipelines to analyse sequencing data from wastewater. Accurate calling of mutations is critical in this process and in the assignment of circulating variants; yet, to date, the performance of variant-calling algorithms in wastewater samples has not been investigated. To address this, we compared the performance of six variant callers (VarScan, iVar, GATK, FreeBayes, LoFreq and BCFtools), used widely in bioinformatics pipelines, on 19 synthetic samples with known ratios of three different SARS-CoV-2 variants of concern (VOCs) (Alpha, Beta and Delta), as well as 13 wastewater samples collected in London between the 15th and 18th December 2021. We used the fundamental parameters of recall (sensitivity) and precision (specificity) to confirm the presence of mutational profiles defining specific variants across the six variant callers. Our results show that BCFtools, FreeBayes and VarScan found the expected variants with higher precision and recall than GATK or iVar, although the latter identified more expected defining mutations than other callers. LoFreq gave the least reliable results due to the high number of false-positive mutations detected, resulting in lower precision. Similar results were obtained for both the synthetic and wastewater samples.

Rapid Assessment of SARS-CoV-2 Variant-Associated Mutations in Wastewater Using Real-Time RT-PCR.

Within months of the COVID-19 pandemic being declared on March 20, 2020, novel, more infectious variants of SARS-CoV-2 began to be detected in geospatially distinct regions of the world. With international travel being a lead cause of spread of the disease, the importance of rapidly identifying variants entering a country is critical. In this study, we utilized wastewater-based epidemiology (WBE) to monitor the presence of variants in wastewater generated in managed COVID-19 quarantine facilities for international air passengers entering the United Kingdom. Specifically, we developed multiplex reverse transcription quantitative PCR (RT-qPCR) assays for the identification of defining mutations associated with Beta (K417N), Gamma (K417T), Delta (156/157DEL), and Kappa (E154K) variants which were globally prevalent at the time of sampling (April to July 2021). The assays sporadically detected mutations associated with the Beta, Gamma, and Kappa variants in 0.7%, 2.3%, and 0.4% of all samples, respectively. The Delta variant was identified in 13.3% of samples, with peak detection rates and concentrations observed in May 2021 (24%), concurrent with its emergence in the United Kingdom. The RT-qPCR results correlated well with those from sequencing, suggesting that PCR-based detection is a good predictor for variant presence; although, inadequate probe binding may lead to false positive or negative results. Our findings suggest that WBE coupled with RT-qPCR may be used as a rapid, initial assessment to identify emerging variants at international borders and mass quarantining facilities. IMPORTANCE With the global spread of COVID-19, it is essential to identify emerging variants which may be more harmful or able to escape vaccines rapidly. To date, the gold standard to assess variants circulating in communities has been the sequencing of the S gene or the whole genome of SARS-CoV-2; however, that approach is time-consuming and expensive. In this study, we developed two duplex RT-qPCR assays to detect and quantify defining mutations associated with the Beta, Gamma, Delta, and Kappa variants. The assays were validated using RNA extracts derived from wastewater samples taken at quarantine facilities. The results showed good correlation with the results of sequencing and demonstrated the emergence of the Delta variant in the United Kingdom in May 2021. The assays developed here enable the assessment of variant-specific mutations within 2 h after the RNA extract was generated which is essential for outbreak rapid response.

Epstein-Barr Virus Genome Deletions in Epstein-Barr Virus-Positive T/NK Cell Lymphoproliferative Diseases.

The main target cells for Epstein-Barr virus (EBV) infection and persistence are B lymphocytes, although T and NK cells can also become infected. In this paper, we characterize the EBV present in 21 pediatric and adult patients who were treated in France for a range of diseases that involve infection of T or NK cells. Of these 21 cases, 5 pediatric patients (21%) and 11 adult patients (52%) were of Caucasian origin. In about 30% of the cases, some of the EBV genomes contain a large deletion. The deletions are different in every patient but tend to cluster near the BART region of the viral genome. Detailed investigation of a family in which several members have persistent T or NK cell infection by EBV indicates that the virus genome deletions arise or are selected independently in each individual patient. Genome sequence polymorphisms in the EBV in these T or NK cell diseases reflect the geographic origin of the patient and not a distinct type of EBV (the 21 cases studied included examples of both type 1 and type 2 EBV infection). Using virus produced from type 1 or type 2 EBV genomes cloned in bacterial artificial chromosome (BAC) vectors, we demonstrate infection of T cells in cord blood from healthy donors. Our results are consistent with transient infection of some T cells being part of normal asymptomatic infection by EBV in young children. IMPORTANCE EBV contributes to several types of human cancer. Some cancers and nonmalignant lymphoproliferative diseases involving T or NK cells contain EBV. These diseases are relatively frequent in Japan and China and have been shown sometimes to have deletions in the EBV genome in the disease cells. We identify further examples of deletions within the EBV genome associated with T or NK cell diseases, and we provide evidence that the virus genomes with these deletions are most likely selected in the individual cases, rather than being transmitted between people during infection. We demonstrate EBV infection of cord blood T cells by highly characterized, cloned EBV genomes and suggest that transient infection of T cells may be part of normal asymptomatic infection by EBV in young children.

A comparison of precipitation and filtration-based SARS-CoV-2 recovery methods and the influence of temperature, turbidity, and surfactant load in urban wastewater.

Wastewater-based epidemiology (WBE) has become a complimentary surveillance tool during the SARS-CoV-2 pandemic. Viral concentration methods from wastewater are still being optimised and compared, whilst viral recovery under different wastewater characteristics and storage temperatures remains poorly understood. Using urban wastewater samples, we tested three viral concentration methods; polyethylene glycol precipitation (PEG), ammonium sulphate precipitation (AS), and CP select™ InnovaPrep® (IP) ultrafiltration. We found no major difference in SARS-CoV-2 and faecal indicator virus (crAssphage) recovery from wastewater samples (n = 46) using these methods, PEG slightly (albeit non-significantly), outperformed AS and IP for SARS-CoV-2 detection, as a higher genome copies per litre (gc/l) was recorded for a larger proportion of samples. Next generation sequencing of 8 paired samples revealed non-significant differences in the quality of data between AS and IP, though IP data quality was slightly better and less variable. A controlled experiment assessed the impact of wastewater suspended solids (turbidity; 0-400 NTU), surfactant load (0-200 mg/l), and storage temperature (5-20 °C) on viral recovery using the AS and IP methods. SARS-CoV-2 recoveries were >20% with AS and <10% with IP in turbid samples, whilst viral recoveries for samples with additional surfactant were between 0-18% for AS and 0-5% for IP. Turbidity and sample storage temperature combined had no significant effect on SARS-CoV-2 recovery (p > 0.05), whilst surfactant and storage temperature combined were significant negative correlates (p < 0.001 and p < 0.05, respectively). In conclusion, our results show that choice of methodology had small effect on viral recovery of SARS-CoV-2 and crAssphage in wastewater samples within this study. In contrast, sample turbidity, storage temperature, and surfactant load did affect viral recovery, highlighting the need for careful consideration of the viral concentration methodology used when working with wastewater samples.

Comparative analysis of the chicken IFITM locus by targeted genome sequencing reveals evolution of the locus and positive selection in IFITM1 and IFITM3.

BACKGROUND: The interferon-induced transmembrane (IFITM) protein family comprises a class of restriction factors widely characterised in humans for their potent antiviral activity. Their biological activity is well documented in several animal species, but their genetic variation and biological mechanism is less well understood, particularly in avian species. RESULTS: Here we report the complete sequence of the domestic chicken Gallus gallus IFITM locus from a wide variety of chicken breeds to examine the detailed pattern of genetic variation of the locus on chromosome 5, including the flanking genes ATHL1 and B4GALNT4. We have generated chIFITM sequences from commercial breeds (supermarket-derived chicken breasts), indigenous chickens from Nigeria (Nsukka) and Ethiopia, European breeds and inbred chicken lines from the Pirbright Institute, totalling of 206 chickens. Through mapping of genetic variants to the latest chIFITM consensus sequence our data reveal that the chIFITM locus does not show structural variation in the locus across the populations analysed, despite spanning diverse breeds from different geographic locations. However, single nucleotide variants (SNVs) in functionally important regions of the proteins within certain groups of chickens were detected, in particular the European breeds and indigenous birds from Ethiopia and Nigeria. In addition, we also found that two out of four SNVs located in the chIFITM1 (Ser36 and Arg77) and chIFITM3 (Val103) proteins were simultaneously under positive selection. CONCLUSIONS: Together these data suggest that IFITM genetic variation may contribute to the capacities of different chicken populations to resist virus infection.

Sequence Variation of Epstein-Barr Virus: Viral Types, Geography, Codon Usage, and Diseases.

One hundred thirty-eight new Epstein-Barr virus (EBV) genome sequences have been determined. One hundred twenty-five of these and 116 from previous reports were combined to produce a multiple-sequence alignment of 241 EBV genomes, which we have used to analyze variation within the viral genome. The type 1/type 2 classification of EBV remains the major form of variation and is defined mostly by EBNA2 and EBNA3, but the type 2 single-nucleotide polymorphisms (SNPs) at the EBNA3 locus extend into the adjacent gp350 and gp42 genes, whose products mediate infection of B cells by EBV. A small insertion within the BART microRNA region of the genome was present in 21 EBV strains. EBV from saliva of U.S. patients with chronic active EBV infection aligned with the wild-type EBV genome with no evidence of WZhet rearrangements. The V3 polymorphism in the Zp promoter for BZLF1 was found to be frequent in nasopharyngeal carcinoma cases from both Hong Kong and Indonesia. Codon usage was found to differ between latent and lytic cycle EBV genes, and the main forms of variation of the EBNA1 protein have been identified.IMPORTANCE Epstein-Barr virus causes most cases of infectious mononucleosis and posttransplant lymphoproliferative disease. It contributes to several types of cancer, including Hodgkin's lymphoma, Burkitt's lymphoma, diffuse large B cell lymphoma, nasopharyngeal carcinoma, and gastric carcinoma. EBV genome variation is important because some of the diseases associated with EBV have very different incidences in different populations and geographic regions, and differences in the EBV genome might contribute to these diseases. Some specific EBV genome alterations that appear to be significant in EBV-associated cancers are already known, and current efforts to make an EBV vaccine and antiviral drugs should also take account of sequence differences in the proteins used as targets.

Accurate characterization of the IFITM locus using MiSeq and PacBio sequencing shows genetic variation in Galliformes.

BACKGROUND: Interferon inducible transmembrane (IFITM) proteins are effectors of the immune system widely characterized for their role in restricting infection by diverse enveloped and non-enveloped viruses. The chicken IFITM (chIFITM) genes are clustered on chromosome 5 and to date four genes have been annotated, namely chIFITM1, chIFITM3, chIFITM5 and chIFITM10. However, due to poor assembly of this locus in the Gallus Gallus v4 genome, accurate characterization has so far proven problematic. Recently, a new chicken reference genome assembly Gallus Gallus v5 was generated using Sanger, 454, Illumina and PacBio sequencing technologies identifying considerable differences in the chIFITM locus over the previous genome releases. METHODS: We re-sequenced the locus using both Illumina MiSeq and PacBio RS II sequencing technologies and we mapped RNA-seq data from the European Nucleotide Archive (ENA) to this finalized chIFITM locus. Using SureSelect probes capture probes designed to the finalized chIFITM locus, we sequenced the locus of a different chicken breed, namely a White Leghorn, and a turkey. RESULTS: We confirmed the Gallus Gallus v5 consensus except for two insertions of 5 and 1 base pair within the chIFITM3 and B4GALNT4 genes, respectively, and a single base pair deletion within the B4GALNT4 gene. The pull down revealed a single amino acid substitution of A63V in the CIL domain of IFITM2 compared to Red Jungle fowl and 13, 13 and 11 differences between IFITM1, 2 and 3 of chickens and turkeys, respectively. RNA-seq shows chIFITM2 and chIFITM3 expression in numerous tissue types of different chicken breeds and avian cell lines, while the expression of the putative chIFITM1 is limited to the testis, caecum and ileum tissues. CONCLUSIONS: Locus resequencing using these capture probes and RNA-seq based expression analysis will allow the further characterization of genetic diversity within Galliformes.

Evaluation of the genotoxicity of 10 selected dietary/environmental compounds with the in vitro micronucleus cytokinesis-block assay in an interlaboratory comparison.

Complex exposure to xenobiotics is one of the reasons for the reported increase of respiratory diseases, cancer and immunological disturbances. Among such xenobiotics there are food mutagens whose effects on human health in the low level and/or chronic exposure still remains unknown. In the present manuscript, the compounds ethanol (EtOH), 4-hydroxynonenal (4-HNE), malondialdehyde (MDA), 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 3,3',4,4'-tetrachlorobiphenyl (PCB 153), benzo[a]pyrene (BaP), 2-amino-3-methylimidazol[4,5-f]quinoline (IQ), 2-amino-1-methyl-6-phenylimidazol[4,5-b]pyridine (PhIP), N-Nitrosodimethylamine (NDMA) and acrylamide (AA) were evaluated in an interlaboratory comparison in the in vitro cytokinesis-block micronucleus assay (CBMN) with objective of assessing the induction of micronuclei, buds and nucleoplasmic bridges in dose responses. Statistically significant increase in MNBN frequency in binucleated cells was recorded by both laboratories for the compound PhIP (2.5µM). The compounds PCB (250 microM) and AA (500 microM) induced statistically significant increase of MNBN although it was recorded by one of the two laboratories. Induction of buds and nucleoplasmic bridges was only observed for BaP (100 microM) and AA (500 microM) by one of the laboratories. Data generated in this study may assist in the interpretation of the mother/newborn biomonitoring study being carried out within project NewGeneris and will contribute to overall knowledge on the genotoxic potential of dietary/environmental toxicants.