Long-read sequencing has revolutionised metagenomics through improved metagenome assembly, taxonomic classification and functional characterisation. Automation can enhance the throughput, reproducibility, and accuracy of library preparation. However, the validation of automated library preparation protocols remains undetermined for metagenomic workflows, which are particularly sensitive to methodological perturbation. Here, we compare long-read metagenomic sequencing of environmental samples through parallel manual and automated protocols. Although automated library preparation led to minor reduction in read and contig lengths, taxonomic classification rate and alpha diversity was slightly higher than manual libraries, including the detection of more rare taxa. Despite this, no significant difference in microbial community structure was identified between manual and automated libraries. Despite minor differences in sequencing and classification metrics, automated and manual library preparation resulted in comparable characterization of environmental community metagenomes. These findings demonstrate the suitability of automation for high-throughput long-read metagenomics, with broad applicability to automated long-read sequencing for improved efficiency and reproducibility.
The reviewers found this study to be well designed and conducted, with reliable results relevant to researchers working in the field of wastewater surveillance. Very minor comments were made regarding the choice of the four methods.
AbstractMetagenomic analysis has enabled insights into soil community structure and dynamics. Long-read sequencing for metagenomics can enhance microbial ecology by improving taxonomic classification, genome assembly, and functional annotation. However, protocols for purifying high-molecular weight DNA from soil are not yet optimised. We used a statistical design of experiments approach to enhance mechanical lysis of soil samples, increasing the length of purified DNA fragments. Low energy input into mechanical lysis improved DNA integrity, resulting in longer sequenced reads. Our optimized settings of 4 m s−1 for 10 s increased fragment length by 70% compared to the manufacturer’s recommendations. Longer reads from low intensity lysis produced longer contiguous sequences after assembly, potentially improving a range of down-stream analyses. Importantly, there was minimal bias exhibited in the microbial community composition due to lysis efficiency variations. We therefore propose a framework for improving the fragment lengths of DNA purified from diverse soil types, improving soil science research with long-read sequencing.
Metagenomics has been transformative in our understanding of the diversity and function of soil microbial communities. Applying long read sequencing to whole genome shotgun metagenomics has the potential to revolutionise soil microbial ecology through improved taxonomic classification, functional characterisation and metagenome assembly. However, optimisation of robust methods for long read metagenomics of environmental samples remains undeveloped. In this study, Oxford Nanopore sequencing using samples from five commercially available soil DNA extraction kits was compared across four soil types, in order to optimise read length and reproducibility for comparative long read soil metagenomics. Average extracted DNA lengths varied considerably between kits, but longer DNA fragments did not translate consistently into read lengths. Highly variable decreases in the length of resulting reads from some kits were associated with poor classification rate and low reproducibility in microbial communities identified between technical repeats. Replicate samples from other kits showed more consistent conversion of extracted DNA fragment size into read length and resulted in more congruous microbial community representation. Furthermore, extraction kits showed significant differences in the community representation and structure they identified across all soil types. Overall, the QIAGEN DNeasy PowerSoil Pro Kit displayed the best suitability for reproducible long-read WGS metagenomic sequencing, although further optimisation of DNA purification and library preparation may enable translation of higher molecular weight DNA from other kits into longer read lengths. These findings provide a novel insight into the importance of optimising DNA extraction for achieving replicable results from long read metagenomic sequencing of environmental samples.
ABSTRACT 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.
Abstract The lifecycle of Zymoseptoria tritici requires a carefully regulated asymptomatic phase within the wheat leaf following penetration of the mesophyll via stomata. Here we compare the roles in this process of two key fungal signalling pathways, mutants of which were identified through forward genetics due to their avirulence on wheat. Whole‐genome resequencing of avirulent Z. tritici T‐DNA transformants identified disruptive mutations in ZtBCK1 from the kinase cascade of the cell wall integrity (CWI) pathway, and the adenylate cyclase gene ZtCYR1. Targeted deletion of these genes abolished the pathogenicity of the fungus and led to similar in vitro phenotypes to those associated with disruption of putative downstream kinases, both supporting previous studies and confirming the importance of these pathways in virulence. RNA sequencing was used to investigate the effect of ZtBCK1 and ZtCYR1 deletion on gene expression in both the pathogen and host during infection. ZtBCK1 was found to be required for the adaptation to the host environment, controlling expression of infection‐associated secreted proteins, including known virulence factors. Meanwhile, ZtCYR1 is implicated in controlling the switch to necrotrophy, regulating expression of effectors associated with this transition. This represents the first study to compare the influence of CWI and cAMP signalling on in planta transcription of a fungal plant pathogen, providing insights into their differential regulation of candidate effectors during invasive growth.
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.
The world has moved into a new stage of managing the SARS-CoV-2 pandemic with minimal restrictions and reduced testing in the population, leading to reduced genomic surveillance of virus variants in individuals. Wastewater-based epidemiology (WBE) can provide an alternative means of tracking virus variants in the population but decision-makers require confidence that it can be applied to a national scale and is comparable to individual testing data. We analysed 19,911 samples from 524 wastewater sites across England at least twice a week between November 2021 and February 2022, capturing sewage from >70% of the English population. We used amplicon-based sequencing and the phylogeny based de-mixing tool Freyja to estimate SARS-CoV-2 variant frequencies and compared these to the variant dynamics observed in individual testing data from clinical and community settings. We show that wastewater data can reconstruct the spread of the Omicron variant across England since November 2021 in close detail and aligns closely with epidemiological estimates from individual testing data. We also show the temporal and spatial spread of Omicron within London. Our wastewater data further reliably track the transition between Omicron subvariants BA1 and BA2 in February 2022 at regional and national levels. Our demonstration that WBE can track the fast-paced dynamics of SARS-CoV-2 variant frequencies at a national scale and closely match individual testing data in time shows that WBE can reliably fill the monitoring gap left by reduced individual testing in a more affordable way.
Monitoring the spread of viral pathogens in the population during epidemics is crucial for mounting an effective public health response. Understanding the viral lineages that constitute the infections in a population can uncover the origins and transmission patterns of outbreaks and detect the emergence of novel variants that may impact the course of an epidemic. Population-level surveillance of viruses through genomic sequencing of wastewater captures unbiased lineage data, including cryptic asymptomatic and undiagnosed infections, and has been shown to detect infection outbreaks and novel variant emergence before detection in clinical samples. Here, we present an optimised protocol for quantification and sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in influent wastewater, used for high-throughput genomic surveillance in England during the COVID-19 pandemic. This protocol utilises reverse compliment PCR for library preparation, enabling tiled amplification across the whole viral genome and sequencing adapter addition in a single step to enhance efficiency. Sequencing of synthetic SARS-CoV-2 RNA provided evidence validating the efficacy of this protocol, while data from high-throughput sequencing of wastewater samples demonstrated the sensitivity of this method. We also provided guidance on the quality control steps required during library preparation and data analysis. Overall, this represents an effective method for high-throughput sequencing of SARS-CoV-2 in wastewater which can be applied to other viruses and pathogens of humans and animals.
Background Many countries have moved into a new stage of managing the SARS-CoV-2 pandemic with minimal restrictions and reduced testing in the population, leading to reduced genomic surveillance of virus variants in individuals. Wastewater-based epidemiology (WBE) can provide an alternative means of tracking virus variants in the population but is lacking verifications of its comparability to individual testing data. Methods We analysed more than 19,000 samples from 524 wastewater sites across England at least twice a week between November 2021 and February 2022, capturing sewage from >70% of the English population. We used amplicon-based sequencing and the phylogeny based de-mixing tool Freyja to estimate SARS-CoV-2 variant frequencies and compared these to the variant dynamics observed in individual testing data from clinical and community settings. Findings We show that wastewater data can reconstruct the spread of the Omicron variant across England since November 2021 in close detail and aligns closely with epidemiological estimates from individual testing data. We also show the temporal and spatial spread of Omicron within London. Our wastewater data further reliably track the transition between Omicron subvariants BA1 and BA2 in February 2022 at regional and national levels. Interpretation Our demonstration that WBE can track the fast-paced dynamics of SARS-CoV-2 variant frequencies at a national scale and closely match individual testing data in time shows that WBE can reliably fill the monitoring gap left by reduced individual testing in a more affordable way. Funding Department of Health and Social Care, UK, Natural Environmental Research Council, UK, COG-UK Evidence before this study Genomic monitoring of wastewater for SARS-CoV-2 variants has been introduced in several countries and shown to effectively detect the spread of known variants in multiple studies. However, verification of its alignment with individual testing data at a national scale has so far been reported only for Austria, where sampling covered around 5.4million people. Further and larger scale verifications of the reliability of wastewater-based epidemiology (WBE) are needed to increase confidence in its use for public health monitoring. Added value of this study We provide evidence that WBE was able to closely track the spread of the emerging SARS-CoV-2 variant Omicron, as well as its sub lineage dynamics, at a regional and national scale across England. Our sampling covered >70% of the English population, equivalent to 39.4 million people. We thereby demonstrate the scalability of our approach to national levels. We also show how WBE is able to track dynamics in different regions of the UK and at a finer scale within London. Its close alignment, in estimated epidemiological timings, with results from intensive individual testing in the same timeframe provides evidence that wastewater-based monitoring can be a reliable alternative when large scale data from individual testing is not available. Implications of all the available evidence Altogether, evidence is accumulating that WBE is a reliable approach for monitoring SARS-CoV-2 variant dynamics and informing public health measures across spatial scales. ### Competing Interest Statement The authors have declared no competing interest. ### Clinical Protocols ### Funding Statement Funding was provided by DHSC UK (2020\_097) and NERC (NE/V003860/1). This report is independent research funded by the Department of Health and Social Care. COG-UK is supported by funding from the Medical Research Council (MRC) part of UK Research & Innovation (UKRI), the National Institute of Health Research (NIHR) [grant code: MC\_PC_19027], and Genome Research Limited, operating as the Wellcome Sanger Institute. The authors acknowledge use of data generated through the COVID-19 Genomics Programme funded by the Department of Health and Social Care. The views expressed are those of the authors and not necessarily those of the Department of Health and Social Care or UKHSA. ### Author Declarations I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained. Yes The details of the IRB/oversight body that provided approval or exemption for the research described are given below: Use of surplus nucleic acid derived from routine diagnostics and associated patient data was approved through the COG-UK consortium by the Public Health England Research Ethics and Governance Group (R&D NR0195). I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals. Yes I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance). Yes I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable. Yes Wastewater sequencing data are publicly available on the European Nucleotide Archive under Study ID PRJEB55313. The clinical case data used in this study are visualised at . A filtered, privacy conserving version of the lineage-LTLA-week dataset is publicly available online () and gives access to almost all used data, despite a small number of cells having been suppressed to conserve patient privacy.
The fungal wheat pathogen Zymoseptoria tritici causes major crop losses as the causal agent of the disease Septoria tritici blotch. The infection cycle of Z. tritici displays two distinct phases, beginning with an extended symp-tomless phase of 1-2 weeks, before the fungus induces host cell death and tissue collapse in the leaf. Recent evidence suggests that the fungus uses little host-derived nutrition during asymptomatic colonisation, raising questions as to the sources of energy required for this initial growth phase. Autophagy is crucial for the path-ogenicity of other fungal plant pathogens through its roles in supporting cellular differentiation and growth under starvation. Here we characterised the contributions of the autophagy genes ZtATG1 and ZtATG8 to the development and virulence of Z. tritici. Deletion of ZtATG1 led to inhibition of autophagy but had no impact on starvation-induced hyphal differentiation or virulence, suggesting that autophagy is not required for Z. tritici pathogenicity. Contrastingly, ZtATG8 deletion delayed the transition to necrotrophic growth, despite having no influence on filamentous growth under starvation, pointing to an autophagy-independent role of ZtATG8 during Z. tritici infection. To our knowledge, this study represents the first to find autophagy not to contribute to the virulence of a fungal plant pathogen, and reveals novel roles for different autophagy-associated proteins in Z. tritici.
This protocol describes the procedure for transforming plasmids into the plant pathogenic fungus Zymoseptoria tritici by Agrobacterium tumefaciens-mediated transformation. It is based on the protocol first described by Zwiers and De Waard (2001; doi: 10.1007/s002940100216) and utilises A.tumefaciens ternary vectors developed by Sidhu et al. (2015a; doi: 10.1016/j.fgb.2015.04.015), Sidhu et al. (2015b; doi: 10.1016/j.fgb.2015.03.021) and Kilaru et al. (2015; doi: 10.1016/j.fgb.2015.03.018). Section 1 details the recipes for media used in the protocol Section 2 describes how to make competent A. tumefaciens cells Section 3 details the process of transforming A. tumefaciens competent cells with plasmid vectors Section 4 describes the process of A. tumefaciens-mediated Z. tritici transformation
This SOP describes the procedure for generating cDNA libraries from SARS-CoV-2 viral nucleic acid extracts from wastewater samples using an adaptation of the Nimagen EasySeq™ SARS-CoV-2 Protocol v4.01 (https://www.nimagen.com/gfx/Covid19/protcol_NimaGen_covid_wgs_v401.pdf). Version 2.0 of this protocol is an update to accommodate the change to using lyophilized indexes in the EasySeq™ RC-PCR SARS-CoV-2 (novel coronavirus) Whole Genome Sequencing Kit (NimaGen, SKU: RC-COV096). Updates to the SARS-CoV-2 WGS Panel probe mixes have also been made (currently v4.02) to improve coverage of new dominant variants of concern (VoCs), although these can be used with no change to this protocol. For the latest BED files containing primer details for analysis, check the download section of https://www.nimagen.com/covid19. Version 3.0 includes a minor update to the options for post PCR cleanup, following optimization of the composition of the AmpliCleanTM beads provided with the kit by NimaGen. Additionally, advice is provided for sequencing poor quality libraries resulting from a high proportion of SARS-CoV-2 negative samples. The first step of this protocol is an addition to the original method, which aims to remove any PCR inhibitors from extracted RNA samples and potentially concentrate the RNA further, before carrying out reverse transcription to cDNA. During the second step, which follows the original EasySeq™ SARS-CoV-2 Protocol v4.01, with some minor modifications, reverse complement PCR is performed using a panel of probes that will generate tiled amplicons covering the whole SARS-CoV-2 genome. The final step involves an improved library clean up that will remove unused/unincorporated PCR probes, which will be more abundant in samples that contain low viral RNA concentrations such as those from wastewater.