Research note | Open | Published:
Zero-valent iron sand filtration reduces concentrations of virus-like particles and modifies virome community composition in reclaimed water used for agricultural irrigation
BMC Research Notesvolume 12, Article number: 223 (2019)
Zero-valent iron sand filtration can remove multiple contaminants, including some types of pathogenic bacteria, from contaminated water. However, its efficacy at removing complex viral populations, such as those found in reclaimed water used for agricultural irrigation, has not been fully evaluated. Therefore, this study utilized metagenomic sequencing and epifluorescent microscopy to enumerate and characterize viral populations found in reclaimed water and zero-valent iron-sand filtered reclaimed water sampled three times during a larger greenhouse study.
Zero-valent iron-sand filtered reclaimed water samples had significantly less virus-like particles than reclaimed water samples at all collection dates, with the reclaimed water averaging between 108 and 109 and the zero-valent iron-sand filtered reclaimed water averaging between 106 and 107 virus-like particles per mL. In addition, for both sample types, viral metagenomes (viromes) were dominated by bacteriophages of the order Caudovirales, largely Siphoviridae, and genes related to DNA metabolism. However, the proportion of sequences homologous to bacteria, as well as the abundance of genes possibly originating from a bacterial host, was higher in the viromes of zero-valent iron-sand filtered reclaimed water samples. Overall, zero-valent iron-sand filtered reclaimed water had a lower total concentration of virus-like particles and a different virome community composition compared to unfiltered reclaimed water.
The intense use of groundwater for agricultural and other activities has led to substantial aquifer depletions globally [1, 2]. Consequently, demand has grown for technologies, such as zero-valent iron (ZVI) sand filters, to treat alternative irrigation water sources and allow for their use. These filters—initially designed to remove chlorinated compounds from groundwater—are composed of mixtures of sand and zero-valent iron . Currently, they are being developed to remove multiple contaminants, including microorganisms, from diverse water sources [4,5,6,7,8,9,10]. Specifically, ZVI reduces Escherichia coli levels in water [4, 5], and can reduce titers of multiple viruses, including Murine norovirus, and several bacteriophage species [6, 7, 11]. However, ZVI has not been evaluated on its ability to remove complex viral populations, such as those within reclaimed water used for agricultural irrigation.
In reclaimed water, virus-like particles (VLPs) are estimated to be 1000-fold higher than in potable water, with the majority of viruses homologous to bacteriophages . Bacteriophages are among the most abundant biological entities on earth and are critical components in food web-dynamics and nutrient cycling . Moreover, phages can influence its bacterial host’s phenotype through the horizontal transfer of genes [e.g. antibiotic resistance genes (ARGs)] [14,15,16,17]. This is potentially of concern for wastewater treatment plants, which are reported to be reservoirs for ARGs . To address the data gaps noted above, this study evaluated DNA viruses in reclaimed water and ZVI-sand filtered reclaimed water collected during a larger study .
Samples were collected during a larger greenhouse study that evaluated impacts of ZVI-sand filtration on antimicrobials, E. coli and total bacterial communities in reclaimed water used to irrigate lettuce . Briefly, in Summer 2016, 240 L of chlorinated effluent was collected biweekly from a tertiary wastewater treatment plant in the Mid-Atlantic U.S. that employed grinding, activated sludge processing, and secondary clarification, and then stored treated effluent in an open-air lagoon where it was chlorinated prior to land application. After collection, water was delivered to the greenhouse and stored in 189 L storage barrels (Algreen Products Inc., Ontario, Canada).
ZVI-sand filter and filtration process
A commercially-available biosand filter (HydrAid®BioSand Water Filter, NativeEnergy, Burlington, VT, USA) was modified for use in this experiment . Briefly, the filter vessel is made of opaque plastic and has a total volume of ~ 55.5 L. Fine filtration sand, provided with the filter [20, 21] and ZVI particles (Peerless Metal Powders and Abrasives Company, Detroit, MI) were sieved (particle size range of 400 µm to 625 µm) and mixed, generating a 50:50 volume/volume mixture. The porosity of the filter was approximately 0.52 , the average volumetric flow rate was ~ 5.6 L/min, the filtration rate was 18 L/min/m2 and the approximate ZVI contact time was 2.58 min .
To mimic the use of sand filters in agricultural settings where irrigation water is not needed every day due to periodic rain events, reclaimed water was filtered every five days. During each filtration event, a 20 L composite of the stored water was gravity filtered to accommodate the irrigation needs of the greenhouse study. To maintain the filter between filtration events, it was kept submerged in reclaimed water, and right before filtration, the five-day old water within the filter was completely flushed out by pouring 20 L reclaimed water through the filter and discarding it. A new 20 L composite reclaimed water sample was then filtered. From the total ~ 20 L ZVI filtrate, a 1 L subsample of filtrate (ZW sample) was collected, along with 1 L of unfiltered reclaimed water (RW sample), once a month for the present study on 6/21/2016, 7/30/2016, and 8/9/2016.
Viral enumeration was performed using the “filter mount method” . Briefly, aliquots (1 μL for RW and 100 μL for ZW) of formalin fixed samples were suspended in sterile deionized water (total volume of 1000 μL), vacuum filtered onto a 13 mm 0.02-µm Anodisc filter (Whatman, USA), and stained with SYBR Gold (Thermo Fisher Scientific, USA). Triplicate slides for each sample were made and counted within 24 h with an Olympus BX61 microscope (20 fields, 1000× magnification). VLPs were quantified with iVision software and differences in VLP counts between samples types were tested using a paired t-test with Bonferroni correction.
Virome preparation and sequencing
Each sample (1 L) was vacuum filtered through a 0.2 μm membrane filter (Pall Corporation, Port Washington, NY) to remove the cellular fraction and collected in sterile receiving flasks. Viral particles were then concentrated using 100 μL of a Fe solution (4.83 g FeCl3 into 100 mL H2O) and 1 mL of a resuspension buffer (0.1 M EDTA-0.2 M MgCl2-0.2 M ascorbate buffer) . After concentration, DNA was extracted from each concentrate (500 µL) via the AllPrep DNA/RNA Mini Kit (Qiagen, CA, USA). DNA libraries were prepared using the modified Nextera XT protocol. Briefly, each of the viral DNA extracts were used in a tagmentation reaction, followed by 13 cycles of PCR with the Nextera i7 and i5 index primers and 2× Kapa master mix, and then sequenced on the Illumina HiSeq 4000 platform (Illumina, San Diego, CA, USA).
Metagenome assembly and analysis
Viromes were assembled and assigned taxonomy/function as described previously . Briefly, paired-end reads were trimmed, merged, and de novo assembled using Trimmomatic ver. 0.36 (slidingwindow:4:30 minlen:60) , FLASh ver. 1.2.11 , and metaSPAdes ver. 3.10.1 (without read error correction), respectively . MetaGene was then used to predict open reading frames (ORFs) from each assembly . Peptide sequences encoded by the predicted ORFs were queried against the peptide SEED and Phage SEED databases (retrieved 11/17/2017) using protein–protein BLAST (BLASTp ver. 2.6.0+) (E value ≤ 1e−3) [30, 31].
Coverage was calculated for each contig by: (i) recruiting quality-controlled reads to contigs via Bowtie2 ver. 2.3.3 (very sensitive local mode) (ii) processing the BAM file for artificial duplicates using Picard (https://broadinstitute.github.io/picard/), and then (iii) using the “depth” function of Samtools ver. 1.4.1 to compute the per-contig coverage [32, 33]. ORF coverage was denoted by start and stop coordinates within each contig. To normalize the abundances, contig and ORF coverages were divided by the sum of coverage per million, similar to the transcripts per million (TPM) metric used in RNA-Seq . Assignment data were visualized using ggplot2 and heatplus [35, 36].
At each sampling date, the VLPs were significantly (p < 0.05) less abundant in the ZW samples compared to the RW samples (Fig. 1). RW samples contained an average of 1.6 × 109, 6.7 × 108, and 7.0 × 108 VLPs/mL in June, July, and August, respectively. The ZW samples contained an average of 8.6 × 106, 2.8 × 107, and 4.2 × 107 VLPs/mL in June, July, and August, respectively.
Sequencing effort and assembly
Viral DNA was extracted from the six samples; however, it was not possible to obtain enough DNA from the June ZW sample for shotgun sequencing. From the samples that were sequenced, there were a total of 136,267,357 read pairs (Additional file 1: Table S1), with an average of 27,253,471 read pairs per virome (± 3,234,104 SD). Metagenomic assembly produced a total of 825,658 contigs, with an average of 165,132 contigs (± 30,305 SD) and an average of 278,196 ORFs (± 63,500 SD) per virome.
To assess the percentage of functional similarity between RW and ZW viromes, peptide ORFs originating from the same sampling dates (July and August) were clustered using CD-HIT (60% peptide similarity) . In July, 42% of the RW peptide ORFs clustered with 68% of the ZW peptide ORFs. For August, the percentage of shared function increased for the reclaimed sample; 60% of the RW peptide ORFs clustered with 61% of the ZW peptide ORFs (Additional file 2: Fig. S1).
Similar to other virome studies , between 32–38% of contigs could be assigned a taxonomic representative (Additional file 1: Table S2). For both the RW and ZW viromes, viral phyla (51–67%) accounted for the majority of the normalized abundance of the taxonomically assigned contigs, followed by bacteria (11–29%) and unknown (17–23%). However, the proportion of bacteria-assigned contigs was greater in the ZW viromes ( ~ 29%) than the RW (11–17%) (Fig. 2a).
The most abundant viral taxonomic classifications for each virome ( ~ 98% of all viral classified taxa) belonged to the dsDNA phages of the order Caudovirales (Fig. 2b), largely Siphoviridae (51–54%), followed by Myoviridae (28–31%), and Podoviridae (13–16%). The remaining ~ 2% of viral sequences were assigned as unclassified phage and viruses with putative hosts archaea, amoeba, plants, or vertebrates.
For both sample types the SEED Subsystem, DNA metabolism, was the most abundant, accounting for 20–30% of the total normalized abundance assigned to functionally classified peptide ORFs. This was followed by phage elements (11–17%), and protein metabolism (8–10%) (Fig. 3a). Annotated SEED Subsystem assignments were parsed for those assigned as resistant to antibiotics and toxic compounds, which were only between 1 and 2%. Among the putative antibiotic and toxic compound annotations, genes for Beta-lactamase were dominant in both of the August viromes. Additionally, the ZW viromes had a greater normalized abundance than all of the RW viromes for: cobalt-zinc-cadmium resistance, copper homeostasis, multidrug resistance efflux pumps, fluoroquinolone resistance, methicillin resistance in staphylococci, zinc resistance, mercuric reductase, arsenic resistance, and mercury resistance operon (Fig. 3b).
Reclaimed water is an important emerging resource that can help alleviate stress on freshwater systems [38, 39]. However, concerns remain regarding microbiological and chemical constituents persisting in reclaimed water, and whether treatment technologies can be used for further remediation. Here, we found that reclaimed water harbors 108–109 VLPs/mL, comparable to the abundances published in previous similar studies [12, 40]. After ZVI-sand filtration, the number of VLPs was significantly lower at all dates, ranging between 106 and 107 VLPs/mL, a roughly 1–2 log reduction. Previous studies have suggested that virus removal from water during ZVI-sand filtration is likely attributed to adsorption and inactivation via iron oxides within the iron [6, 41]. You et al. suggested that as water flows through a ZVI-sand filter, new iron oxides are formed continuously, generating additional adsorption sites that can extend the life of the filter .
Our findings are similar to recent results on the reduction efficiency of sand filtration alone for viruses ϕX174, MS2, and AiV ( < 1–2 log) and lower than previous studies on ZVI-sand filtration, which reported that φX174 and MS-2 were reduced by 4–6 logs [6, 7]. However, these studies focused largely on the removal of a few specific viral species and, even so, have found that removal efficiencies vary among species . Here, we used epifluorescence microscopy to look at the entire viral population. This includes hundreds to thousands of different species, with a range of capsid sizes and isoelectric points, which may explain the smaller removal efficiency . Moreover, while the log-reduction is lower than expected for the overall population, the total VLP concentration post ZVI-sand filtration is still comparable to those described in well and potable water [12, 43].
In terms of viral taxonomic composition, both RW and ZW viromes were dominated by Siphoviridae, which are abundant in human waste and reclaimed water [12, 44]. These viruses present a unique risk, as the majority of cultured representatives are capable of lysogeny and, thus, may facilitate horizontal gene transfer among bacteria . Additionally, in both sample types the functional profiles were largely composed of genes related to DNA metabolism (Fig. 3a). A previous study of wastewater also found these genes to be highly abundant and attributed this to the elevated metabolic activity within wastewater treatment plants .
In addition, ZW viromes had a greater relative abundance of bacterial assigned contigs. In virome studies, sequences with significant homology to bacteria are sometimes unknown prophage embedded in a bacterial genome present in the database, or phages carrying host genes [12, 46]. During use, ZVI produces ROS, which can promote prophage induction [47, 48]. Thus, during filtration, ROS may stimulate the induction of integrated prophages. However, while the abundance of some genes is higher in the ZW viromes, the overall number of gene copies is likely still higher in the RW sample due to the increased number of total VLPs. Therefore, additional work is necessary to determine whether the bacterial assigned contigs are indeed prophage and whether this impacts the dissemination of bacterial genes in water reuse systems.
Our study was limited by sample size (n = 6), which prevented a rigorous statistical analysis. In addition, we could not obtain enough DNA from the June ZW water sample for shotgun sequencing and, therefore, a comparison between RW and ZW water for June was not possible. Finally, because we were not able to include a second sand-only filter control—due to the set-up of the larger study—we were unable to tease out the effects of ZVI versus sand in terms of virus removal.
zero-valent iron filtered reclaimed water
open reading frame
virus like particle
Konikow LF. Long-term groundwater depletion in the United States. Groundwater. 2015;53(1):2–9.
Konikow LF, Kendy E. Groundwater depletion: a global problem. Hydrogeol J. 2005;13(1):317–20.
Chiu P. Applications of zero-valent iron (ZVI) and nanoscale ZVI to municipal and decentralized drinking water systems—a review. In: Ahuja, S, Hristovski, K, editors. Novel solutions to water pollution. 2013; p. 237–49.
Lee C, Kim JY, Lee WI, Nelson KL, Yoon J, Sedlak DL. Bactericidal effect of zero-valent iron nanoparticles on Escherichia coli. Environ Sci Technol. 2008;42(13):4927–33.
Ingram D, Callahan MT, Ferguson S, Hoover D, Shelton D, Millner P, Camp M, Patel J, Kniel K, Sharma M. Use of zero-valent iron biosand filters to reduce Escherichia coli O157: H12 in irrigation water applied to spinach plants in a field setting. J Appl Microbiol. 2012;112(3):551–60.
You Y, Han J, Chiu PC, Jin Y. Removal and inactivation of waterborne viruses using zerovalent iron. Environ Sci Technol. 2005;39(23):9263–9.
Shi C, Wei J, Jin Y, Kniel KE, Chiu PC. Removal of viruses and bacteriophages from drinking water using zero-valent iron. Separat Purif Technol. 2012;84:72–8.
Fu F, Dionysiou DD, Liu H. The use of zero-valent iron for groundwater remediation and wastewater treatment: a review. J Hazard Mater. 2014;267:194–205.
Tratnyek PG, Scherer MM, Johnson TL, Matheson LJ. Permeable reactive barriers of iron and other zero-valent metals. Environmental Science and Pollution Control Series. 2003; p. 371–422.
Simon F-G, Meggyes T. Removal of organic and inorganic pollutants from groundwater using permeable reactive barriers: part 1. Treatment processes for pollutants. Land Contam Reclamation. 2000;8(2):103–16.
Shearer AE, Kniel KE. Enhanced removal of norovirus surrogates, Murine norovirus and Tulane virus, from aqueous systems by zero-valent iron. J Food Protect. 2018;81:1432–8.
Rosario K, Nilsson C, Lim YW, Ruan Y, Breitbart M. Metagenomic analysis of viruses in reclaimed water. Environ Microbiol. 2009;11(11):2806–20.
Wommack KE, Colwell RR. Virioplankton: viruses in aquatic ecosystems. Microbiol Mol Biol Rev. 2000;64(1):69–114.
Reyes A, Haynes M, Hanson N, Angly FE, Heath AC, Rohwer F, Gordon JI. Viruses in the faecal microbiota of monozygotic twins and their mothers. Nature. 2010;466(7304):334–8.
Willner D, Furlan M, Haynes M, Schmieder R, Angly FE, Silva J, Tammadoni S, Nosrat B, Conrad D, Rohwer F. Metagenomic analysis of respiratory tract DNA viral communities in cystic fibrosis and non-cystic fibrosis individuals. PloS ONE. 2009;4(10):e7370.
Wang X, Kim Y, Ma Q, Hong SH, Pokusaeva K, Sturino JM, Wood TK. Cryptic prophages help bacteria cope with adverse environments. Nat Commun. 2010;1:147.
Brüssow H, Canchaya C, Hardt W-D. Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev. 2004;68(3):560–602.
Rizzo L, Manaia C, Merlin C, Schwartz T, Dagot C, Ploy M, Michael I, Fatta-Kassinos D. Urban wastewater treatment plants as hotspots for antibiotic resistant bacteria and genes spread into the environment: a review. Sci Total Environ. 2013;447:345–60.
Kulkarni P, Raspanti GA, Bui AQ, Bradshaw RN, Kniel KE, Chiu PC, Sharma M, Sapkota A, Sapkota AR. Zerovalent iron-sand filtration can reduce the concentration of multiple antimicrobials in conventionally treated reclaimed water. Environ Res. 2019;172:301–9.
Manz DH. Slow sand filter for use with intermittently flowing water supply and method of use thereof. In: Google Patents. 2000.
Quest T LLC. Hydraid (R) BioSand water filter installation manual. In: C Engineering editor. Grand rapids 2010; p. 1–20.
Center for Affordable Water and Sanitation Technology (CAWAST). How was the biosand filter pore volume determined? 2015. https://biosandfilters.info/faq/how-was-biosand-filter-pore-volume-determined. Accessed 11 Jan 2017.
Williamson KE, Wommack KE, Radosevich M. Sampling natural viral communities from soil for culture-independent analyses. Appl Environ Microbiol. 2003;69(11):6628–33.
John SG, Mendez CB, Deng L, Poulos B, Kauffman AKM, Kern S, Brum J, Polz MF, Boyle EA, Sullivan MB. A simple and efficient method for concentration of ocean viruses by chemical flocculation. Environ Microbiol Rep. 2011;3(2):195–202.
Chopyk J, Allard S, Nasko DJ, Bui A, Mongodin EF, Sapkota AR. Agricultural freshwater pond supports diverse and dynamic bacterial and viral populations. Front Microbiol. 2018;9:792.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–200.
Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27(21):2957–63.
Nurk S, Meleshko D, Korobeynikov A, Pevzner P. metaSPAdes: a new versatile de novo metagenomics assembler. arXiv preprint arXiv:160403071. 2016.
Noguchi H, Park J, Takagi T. MetaGene: prokaryotic gene finding from environmental genome shotgun sequences. Nucleic Acids Res. 2006;34(19):5623–30.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Overbeek R, Fonstein M, D’souza M, Pusch GD, Maltsev N. The use of gene clusters to infer functional coupling. Proc Natl Acad Sci. 1999;96(6):2896–901.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Szcześniak MW, Gaffney DJ, Elo LL, Zhang X. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17(1):13.
Wickham H. ggplot2: elegant graphics for data analysis. Berlin: Springer Science & Business Media; 2009.
Ploner A. Heatplus: heatmaps with row and/or column covariates and colored clusters. R package version. 2012;2(0).
Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.
Parsons LR, Sheikh B, Holden R, York DW. Reclaimed water as an alternative water source for crop irrigation. HortScience. 2010;45(11):1626–9.
Agency USEP. Potable Reuse Compendium, ed. 2017.
Tamaki H, Zhang R, Angly FE, Nakamura S, Hong PY, Yasunaga T, Kamagata Y, Liu WT. Metagenomic analysis of DNA viruses in a wastewater treatment plant in tropical climate. Environ Microbiol. 2012;14(2):441–52.
Ryan JN, Harvey RW, Metge D, Elimelech M, Navigato T, Pieper AP. Field and laboratory investigations of inactivation of viruses (PRD1 and MS2) attached to iron oxide-coated quartz sand. Environ Sci Technol. 2002;36(11):2403–13.
Michen B, Graule T. Isoelectric points of viruses. J Appl Microbiol. 2010;109(2):388–97.
Rinta-Kanto JM, Lehtola MJ, Vartiainen T, Martikainen PJ. Rapid enumeration of virus-like particles in drinking water samples using SYBR green I-staining. Water Res. 2004;38(10):2614–8.
Breitbart M, Hewson I, Felts B, Mahaffy JM, Nulton J, Salamon P, Rohwer F. Metagenomic analyses of an uncultured viral community from human feces. J Bacteriol. 2003;185(20):6220–3.
Suttle CA. Viruses in the sea. Nature. 2005;437(7057):356–61.
Park E-J, Kim K-H, Abell GC, Kim M-S, Roh SW, Bae J-W. Metagenomic analysis of the viral communities in fermented foods. Appl Environ Microbiol. 2011;77(4):1284–91.
Nanda AM, Thormann K, Frunzke J. Impact of spontaneous prophage induction on the fitness of bacterial populations and host-microbe interactions. J Bacteriol. 2015;197(3):410–9.
DeMarini DM, Lawrence BK. Prophage induction by DNA topoisomerase II poisons and reactive-oxygen species: role of DNA breaks. Mutat Res Fundam Mol Mech Mutagen. 1992;267(1):1–17.
ARS and PK designed the experiment. PK and RB collected samples and filtered the RW. JC processed both sample types for viral concentrates, conducted the direct counts, analyzed the metagenomic data, and wrote the manuscript. DN contributed to data analysis. KEK, PC and MS designed the ZVI sand filter. All authors read and approved the final manuscript.
A special thanks to Russell Reynells who contributed to the filter design, Cary Coppock and Eric Handy who helped construct the filter, and Shirley Micallef for use of her greenhouse space. We would also like to thank Eric Wommack for use of his lab space and epifluorescence microscope to conduct the viral direct counts.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the NCBI’s Sequence Read Archive under the accession numbers SAMN09231814-SAMN09231818, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA472060.
Consent for publication
Ethics approval and consent to participate
This work was supported by the United States Department of Agriculture-National Institute of Food and Agriculture, Grant Number 2016-68007-25064, awarded to the University of Maryland School of Public Health, that established CONSERVE: A Center of Excellence at the Nexus of Sustainable Water Reuse, Food and Health. Funding was also provided by the Computation and Mathematics for Biological Networks (COMBINE) NSF Research Traineeship program (DGE-1632976) and USDA Agricultural Research Project “Characterization and Mitigation of Bacterial Pathogens in the Fresh Produce Production and Processing Continuum” (Project Number 8042-32420-006-00-D).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.