Escherichia coli O157:H7 attachment and persistence within root biofilm of common treatment wetlands plants Rachel J. VanKempen-Fryling & Anne K. Camper NOTICE: this is the author’s version of a work that was accepted for publication in Ecological Engineering. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Ecological Engineering, [Vol. 98, January 2017], DOI#:10.1016/j.ecoleng.2016.10.018. VanKempen-Fryling RJ, Camper AK “Escherichia coli O157:H7 attachment and persistence within root biofilm of common treatment wetlands plants,” Ecological Engineering, 2017 Jan;98:64–69. Made available through Montana State University’s ScholarWorks scholarworks.montana.edu 6KRUWFRPPXQLFDWLRQ (VFKHULFKLDFROL2+DWWDFKPHQWDQGSHUVLVWHQFH ZLWKLQURRWELR¿OPRIFRPPRQWUHDWPHQWZHWODQGVSODQWV 5DFKHO-9DQ.HPSHQ)U\OLQJD$QQH.&DPSHUEѽ D&HQWHUIRU%LR¿OP(QJLQHHULQJ0RQWDQD6WDWH8QLYHUVLW\%R]HPDQ0786$ E'HSDUWPHQWRI&LYLO(QJLQHHULQJDQG&HQWHUIRU%LR¿OP(QJLQHHULQJ0RQWDQD6WDWH8QLYHUVLW\%R]HPDQ07 86$ .H\ZRUGV:HWODQGV(VFKHULFKLDFROL2+%LR¿OPV3ODQWV&RORQL]DWLRQ $%675$&7 3DWKRJHQ UHWHQWLRQDQG VXEVHTXHQW UHOHDVHZLWKLQ WKH UKL]RVSKHUHRIZDVWHZDWHU WUHDWPHQWZHWODQGVPD\EHD FRQFHUQ IRU KXPDQKHDOWK7RDGGUHVVWKLVFRQFHUQWKHHQWHULFSDWKRJHQ(VFKHULFKLDFROL2+ZLWKD'V5HGSODVPLGLQVHUWLRQZDV XVHG DV D PRGHO SDWKRJHQLF RUJDQLVP LQ DQ RSHQDLU FKHPRVWDW UHDFWRU ZLWK FRQVWDQW ÀRZ RI VLPXODWHG ZDVWHZDWHU &RORQL]DWLRQDQGSHUVLVWHQFHRIWKHRUJDQLVPZDVWUDFNHGRQURRWVRIWZRREOLJDWHZHWODQGSODQWVSHFLHV&DUH[XWULFXODWDDQG 6FKRHQRSOHFWXVDFXWXVRULJLQDOO\JURZQLQSLORWVFDOHZHWODQGUHDFWRUV7HÀRQQ\ORQVWULQJFOHDQDQGZLWKH[LVWLQJLQGLJHQRXV ELR¿OP ZDV XVHG DV DQ LQHUW VXUIDFH FRQWURO (SLÀXRUHVFHQFHPLFURVFRS\ DQG T3&5ZHUH XVHG WR YHULI\ ( FROL2+ DEXQGDQFHIRUXSWRZHHN,QLWLDODWWDFKPHQWZDVVHHQRQDOOVXUIDFHVZLWKFRORQL]DWLRQGHFUHDVLQJWKURXJKZHHNT3&5 VKRZHGSUHIHUHQWLDODVVRFLDWLRQRI WKHSDWKRJHQZLWKURRWVRYHU WKHQ\ORQ7KHUHZDVDVLJQL¿FDQWGLIIHUHQFHEHWZHHQSODQW W\SH6DFXWXVVKRZHGVLJQL¿FDQWO\KLJKHUQXPEHUVFRPSDUHGWR&XWULFXODWD(FROL2+ELQGLQJDQGSHUVLVWHQFHRQ URRWVXUIDFHVPD\EHDPHDQVRIVXUYLYDOLQWUHDWPHQWZHWODQGV ,1752'8&7,21 :DVWHZDWHU WUHDWPHQW LV LPSRUWDQW IRU WKH HQYLURQPHQW DQG KXPDQ KHDOWK 0DQ\ DOWHUQDWLYH PHWKRGV WR FRQYHQWLRQDO WUHDWPHQW KDYH EHHQ UHVHDUFKHG 7RQGHUVNL   LQFOXGLQJ XVLQJ WUHDWPHQWZHWODQGV 7: 7:DUH WKH DGRSWLRQ DQG DGDSWDWLRQRIDZHWODQG¶V DELOLW\ WR UHPRYH FRPPRQ FRQWDPLQDQWV LQ ZDWHU ,Q JHQHUDO 7: KDYH D QRQSHUPHDEOH OD\HU WR SUHYHQW OHDFKLQJ DQG DUH SODQWHG ZLWK VLQJOH RU PL[HG VSHFLHV SODQWV LQ D PHGLXP RI JUDYHO RU VRLO $NUDWRV DQG 7VLKULQW]LV  %DVNDU HW DO   7: FDQYDU\LQWHUPVRI ZDVWHZDWHU ORDGLQJ DQG UHWHQWLRQ WLPH .RWWL HW DO $XVWLQ  )RU ZDVWHZDWHU WUHDWPHQW 7:V SURYLGH ERWK DHURELF DQG DQDHURELF PLFURQLFKHV )DXOZHWWHU HW DO 3ODQWURRWVUHOHDVHR[\JHQDVZHOODVRWKHU H[XGDWHV LQWR WKH VXEVXUIDFH SURYLGLQJ DQ HQYLURQPHQWIRUDHURELFEDFWHULDWRVXUYLYHDWGHSWKV EHORZ WKH GLVVROYHG R[\JHQ OHYHO &DUEDOOHLUD HW DO   7KLV DOORZV DHURELF SURFHVVHV VXFK DV QLWUL¿FDWLRQDQGSKRVSKRUXVUHPRYDOOHDYLQJ WKHDQDHURELFPDWUL[WRFRQWLQXHRWKHUDQDHURELFUHDFWLRQV VXFKDVGHQLWUL¿FDWLRQ +HQFKHWDO7D\ORUHWDO   7KH URRW ELR¿OP RI 7: PD\ DOVR UHGXFH SDWKRJHQLF EDFWHULD SUHVHQW LQ WKH ZDVWHZDWHU /LFHQFH HW DO  %RXWLOLHUHWDO=KDQJHWDO ([LVWLQJELR¿OP FRPSHWLWLRQ DQG SUHGDWLRQ LQ WKH UKL]RVSKHUH OLPLWV DWWDFKPHQWDQGVHWWOLQJ %RXWLOLHUHWDO-DVSHUHWDO  +RZHYHU WKHUH LV H[WHQVLYHUHVHDUFKUHYHDOLQJ WKDW SDWKRJHQV PD\ VXUYLYH DQG EH UHOHDVHG LQ HIÀXHQW ZDWHU %DUFLQDDQG$UDQD/LXHWDO3DUXFK  (YHQDORZFRXQWUHOHDVHVXFKDVFHOOVIRU(VFKHULFKLD FROL 2+ LV HQRXJK WR FDXVH KXPDQ LQIHFWLRQ DQG GLVHDVH %LWWRQ86(3$  ( FROL 2+ LV D IRRG DQG ZDWHUERUQH SDWKRJHQ WUDQVPLWWHG YLD WKH IHFDORUDO URXWH XVXDOO\ WKURXJK FRQWDPLQDWHG ZDWHU RU YHJHWDEOHV /LFHQFH HW DO  'LQX DQG %DFK  2QJHQJ HW DO   3UHYLRXV UHVHDUFK KDV VKRZQ D GHFUHDVH RI ( FROL 2+ WKURXJK 7:V KRZHYHU WKHVH VWXGLHV DUH VROHO\ EDVHG RQ SODWH FRXQWV /LX HW DO  2QJHQJ HW DO  *DUFtD HW DO   'LUHFW PLFURVFRS\ RU PROHFXODU WHFKQLTXHV IRU HQXPHUDWLRQ ZRXOG OLNHO\ UHYHDO PXFK KLJKHU VXUYLYDO DV ZHOO DV WKH DELOLW\ IRU ODWHU UHOHDVH 7KLV LV GXH WR D SKHQRPHQRQ NQRZQ DV WKH YLDEOH EXW QRQFXOWXUDEOH 9%1&  VWDWH LQ ZKLFK EDFWHULD DUH QR ORQJHU FXOWXUDEOH RQDSODWHEXW FDQEH UHVXVFLWDWHG IURP WKH HQYLURQPHQW ZKHQ UHWXUQHG WR JURZWKSURPRWLQJ FRQGLWLRQV %DUFLQDDQG$UDQD'LQXDQG%DFK 3DUXFK  Fig. 1. Representative imagesof the fourexperimental surface types takenat2, 24and168h. Fluorescencewasobservedunder558nmwavelength.Due tonaturalfluorescence of photosystem II chlorophyll in the plant tissue, plant surface is shown in C. utriculata and S. acutus images. No stains or dyes were used in Teflon nylon imaging under the assumption stationary cells were attached to a surface. Scale bar represents 50m. Another serious limitation is the focus on influent to efflu- ent measurements, assuming that E. coli O157:H7 is not surviving within the TW. Long-term lab studies in river and reservoir water found E. coli O157:H7 can persist up to 6 months (Liu et al., 2009). E. coli O157:H7 may also be surviving within the biofilm matrix or incorporated into the plants, as seen in vegetable production (Fu et al., 2005; Dinu and Bach 2011) E. coli O157:H7 can invade plant tissuesof leafygreensandstalks (Deeringet al., 2012).Although this is not a health concern in TW as with edible plants, E. coli O157:H7 within vascular openings of plants could provide means of survival and potential release. The hypothesis was that indigenous biofilm and living plant tis- sue (with an established biofilm) impact the persistence of E. coli O157:H7 in a simulated TW. Two plant species, C. utriculata and S. acutus, were used in hydroponic reactors to determine interactions with the roots. In previous research these two obligate wetland species were the top performers for removal efficiency of carbon, nitrogen and sulfate in pilot-scale wetland systems (Taylor et al., 2011). Inert Teflon nylon material was used as a control both with and without a pre-established biofilm. Pathogen persistence on these surfaces was assessed using qPCR and microscopy. 2. Materials and methods 2.1. Growth of E. coli O157:H7 The E. coli O157:H7 was isolated from an outbreak in 1993 (Pyle et al., 1995). A constitutively expressed DsRed plasmid with a car- benicillin resistance selective gene (Clonetech) was incorporated (Klayman et al., 2009) for detection by fluorescent microscopy. Tryptic Soy Broth (TSB) was made using non-casein enzymatic digest of protein (Klayman et al., 2009) for plasmid retention and expression. E. coli O157:H7 was grown overnight to approximately 108 cellsmL−1 as determined by plate count. Survival and plas- mid stability in simulated wastewater (SW; Supplemental Table 1) was verified by plate counts, qPCR, and epifluorescencemicroscopy (VanKempen-Fryling et al., 2015). 2.2. Plant growth and strain collection Shoots with substantial root mass of C. utriculata and S. acutus grown in a controlled greenhouse environmentwereharvested and washed with nonsterile H2O to remove excess planting medium and loosely associated bacteria while keeping the attached biofilm. The rinsed plants (49.3mL±26.4mL, determined by the mass displacement) were transferred to 300mL hydroponic chemostat reactors under 24h light at 25 ◦C and grown for 3–7days in SW. Comparisons were made between living roots to an inert sur- face of Teflon nylon strings (D’Addario; GHS strings) similar in dimensions to the roots (0.79, 0.81 and 1.01mm diameter). Clean nylon strings were added to chemostats with SW immediately prior to E. coli O157:H7 inoculation. Alternatively, strings were added to effluent of the established plant mesocosms for 24h in batch, followed by the addition of SW (flow conditions) for 1 week to encourage biofilm growth prior to pathogen inoculation. Total weight of the nylon strings, cut into 17.5 cm lengths (half immersed), was 2.94g. 2.3. Determination of surface area To estimate the E. coli O157:H7 surface area attachment, stan- dard curves were generated relating wet mass to root surface area. Representative roots were excised, imaged, and analyzed with the Rootscan® program to determine the surface area. The wet mass was plotted against surface area to generate standard curves (Sup- plemental Fig. 1). C. utriculata’s two distinct root types (taproot and branched) formed unique standard curves; DNA extraction for C. utriculata samples were separated by root type. Teflon surface area was estimated by treating the strings as open ended cylinders. The 0% 20% 40% 60% 80% 100% 0 50 100 150 200 per cen t cl um ps per tot al s ize ra nge s 0% 20% 40% 60% 80% 100% 0 20 40 60 80 100 120 140 160 180 Per cen t cl um ps per tot al s ize ra nge s Single Cels 2-10 cels 11-100 cels >100 cels 0% 10% 20% 30% 40% 50% 60% 70% 80% 0 50 100 150 200 per cen t cl um ps per tot al s ize ra nge s 0% 20% 40% 60% 80% 100% 0 20 40 60 80 100 120 140 160 180Pe rce nt c lu mps pe r tot al s ize ra nge s Time (hours) A. D. C. B. Fig. 2.Distribution of cel clumping over total counts per time point for each surface: (A)C. utriculata, (B)S. acutus, (C) clean Teflon nylon, and (D) biofilm Teflon nylon. Root locations were combined for a total mean. Bars represent standard deviation from al experimental replicates. total surface area of the 24 sample strings, 8 of each diameter was 53.32 cm2. 2.4. Reactor series 20 mL of overnightE. coliO157:H7 culture was spun at 4816 x g in triplicate, washed in nanopure water, and concentrated into 1 mL. For each triplicate experimental condition (C. utriculata,S. acutus,clean nylon, and established biofilm nylon), 1 mL of inocu- lum was added and mixed via pipetting and rotating the plant or nylon. One additional non-inoculated reactor served as a control. Flow was at 3 mL min−1(100 min hydraulic residence times (HRT)) for one week. A time zero sample was taken from each reactor immediately after inoculation. Effluent samples were colected at 45-min intervals for 6 h. One replicate reactor was destructively sampled at 135 min, approximately 1 HRT (initial attachment). Another effluent sample was taken at 24 h and the second replicate reactor destructively sampled. The third replicate was destructively sampled at one week. The uninoculated control reactor effluent was sampled at 0 and 6 h and destructively sampled at one week. Al water samples were frozen for further analysis. Submerged roots from destructively sampled reactors were divided in half by length for differentiation of basal and root tip, rinsed with nanopure water to remove unattached cels, placed into 2 mL microcentrifuge tubes and frozen for DNA extraction. The nylon samples were taken in entirety, placed into 50 mL Falcon tubes, and frozen. A representative root or nylon sample was puled through a 2 mm square glass capilary approximately 90–110 mm in length (Wale Apparatus, PA) using a wire hook (Aircraft Spruce), the interior of the capilary rinsed with nanopure water to slough off unattached cels, and attachment observed by epifluorescence microscopy. 2.5. Epifluorescence microscopy The glass capilaries were viewed under 600×magnification (Nikon; Melvile, NY) and 558 nm wavelength for DsRed plasmid excitation on a Nikon Eclipse E-800 microscope with a cooled CCD camera (Photometrics). The root tip, mid-root, and top root were 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 2 24 168 2 24 168 2 24 168 2 24 168 2 24 168 2 24 168 Branched Thick Inactive Active Biofilm Clean C. utriculata S. acutus Nylon rel ativ e a bun da nce (g ene co pie s/c m2 ) Fig. 3.Relative abundance ofE. coliO157:H7 per surface area normalized to back- ground fluorescence. The data is differentiated over time and by root type or nylon condition for each respective surface. Table 1 DNA extraction quantification. Sample Type Concentration (ng/uL) A260/230 C. utriculatawater C. utriculataroot S. acutuswater S. acutusroot Nylon clean water Nylon clean surface Nylon biofilm water Nylon biofilm surface 80.08±40.56 96.38±30.02 57.34±51.16 72.09±37.17 59.81±29.91 48.37±36.63 72.22±45.80 40.03±24.86 2.01±0.21 2.02±0.15 1.49±0.61 1.88±0.25 2.04±0.21 1.94±0.70 1.95±0.48 1.99±0.34 scanned and 20 images taken in each area for samples colected at 2, 24 and 168 h. Natural photosystem I chlorophyl fluorescence alowed the plant surface to be viewed without stain. The images were semi quantitatively analyzed (MetaMorph®Image Analysis Software) to determine number of cels and colony clump size using the average pixel area forE. coliO157:H7 single cel fluorescence (Wilson et al., 2004; Behnke et al., 2011). Thresholding was used to subtract background autofluorescence and remaining objects were counted. 2.6. DNA extraction and qPCR analysis Thawed water samples were spun down at 4816×gfor 20 min at 4◦C. The peleted sample was resuspended and transferred to a PowerBead tube from the PowerSoil®DNA Isolation Kit (MoBio Laboratories; Carlsbad, CA, USA). Mass of the roots was determined prior to transfer to a PowerBead tube. DNA extraction folowed the kit protocol except a Fastprep®-24 bead beater (M.P. Biomedicals; Santa Ana, CA, USA) was used at 5.5 m sec−1for 45 s. Cel mate- rial on the nylon strings was removed by adding glass beads to the 5 mL mark of the 50 mL conical Falcon tube with 10 mL auto- claved nanopure H2O, and vortexed for 60 s. The liquid fraction was removed and centrifuged at 4816×gfor 5 min. The supernatant was removed and the pelet resuspended and added to a PowerBead tube for extraction as previously described. DNA was isolated in 100 L of DEPC water (Ambion®; Grand Island, NY, USA) and quan- tified using a Nanodrop® ND-1000 spectrophotometer (Thermo Scientific; Wilmington, DE, USA). Average quantification and purity per sample type is listed inTable 1. qPCR primers targeted the Shiga toxin 2 gene unique to the O157:H7 serotype (Jothikumar and Griffiths, 2002). Control exper- iments using plate counts verified that the cels had one copy of the fStx1andStx2genes (Shaikh and Tarr, 2003). A 25 L reaction mixture was prepared using SYBR Green Master Mix (Kapa Biosys- tems; Wilmington, MA, USA), 0.5 mM primer, and 3.1×10−3mM DNA diluted from the extraction. Programs for qPCR can be found VanKempen-Fryling et al., 2015. 2.7. Data analysis Multifactor analysis of variance (ANOVA) was performed to determine statistical significance between surface types and time. Comparisons were analyzed between plants and Teflon nylon strings, the root tip or basal root, established biofilm compared to clean nylon, and time. General linear models (GLM) were used to compare interactions using Minitab 17®. Statistical significance between specific factors for each variable were identified by pair- wise Tukey’spost hocmultiple comparison test using a 0.05p-value cut-off. The experimental replicates were input as random variables nested within time. 3. Results Initial washout ofE. coliO157:H7 during the first four HRT was similar for al systems (Supplemental Fig. 2). A twenty-four hour sample showed the continued release ofE. coliO157:H7 from the inoculated reactors, which would result from detachment of pre- viously attached cels. E. coliO157:H7 could be imaged on the surfaces throughout the time course (Fig. 1). Initial attachment was mainly single cels and there was no observable difference in abundance by location or sur- face type. At 24 h, microcolony formation was seen. The final time point at 168 h showed a lower, sustainedE. coliO157:H7 cel pres- ence and more variable microcolony formation. These observations were substantiated with the MetaMorph®Image Analysis software (Fig. 2) and reported as a percentage of the total events for al size ranges from smal and medium clusters to large, >100 cel forma- tions. The number of clumps greater than single cels increased on al surfaces except the biofilm-colonized nylon, where increases were limited to clumps >100 cels. This is more dramatic inS. acutus (23 to 36 events) and clean nylon experiments (96 to 109 events). Single cels were more varied between replicates, as seen through the large error bars. The total clump counts decreased from 2 to 168 h (22.6 average events to 13.3 average events), however the largest formations increased (average increase of 1 event for >100 cel clump size). qPCR determined relative target cel abundance on inoculated vs non-inoculated surfaces over time. For the root experiments, location (root tip or basal root) and condition (branched or thick forC. utriculata; red/live or black/dead coloration forS. acutus) of the root were analyzed. For nylon experiments, the biofilm presence/absence was compared. For al experimental surfaces (Fig. 3), the population decreased during the period between 2 and 168 h, averaging a reduction from 704 gene copies cm−2to 55 gene copies cm−2(P < 0.001). Attachment onS. acutuswas signifi- cantly higher (average 493 gene copies cm−2) thanC. utriculataby 209 gene copies cm=2(P = 0.002), as wel as both nylon replicates (P = 0.021) in al time points. Nylon replicates andC. utriculatawere statisticaly equivalent with the lowestE. coliO157:H7 attachment across al time points. ForC. utriculata,E. coliO157:H7 preferred branched roots to uniform taproots with an average difference of 82.1 gene copies cm−2to 281 gene copies cm−2(P < 0.001) across al time points. InS. acutus, no meaningful effect of live roots was observed. Localization was taken into account between root tip and the basal roots, but there was no preference. The two abiotic surfaces were also compared to identify any difference with respect to preexisting biofilm alone. There was no significant difference inE. coliO157:H7 attachment. A slight inter- action between surface and time was seen from initial attachment to 24h in clean nylon over biofilm-colonized nylon (Supplemental Fig. 3). 4. Discussion Epifluorescence microscopy and qPCR showed that E. coli O157:H7 attached and persisted on all surfaces under flow con- ditions for the duration of the experiment. Attachment was significantly lower on abiotic nylon controls compared to S. acu- tus roots. This supports the hypothesis that there is a benefit for survival or potential growth for bacterial pathogens within an established root biofilm (Tyler and Triplett 2008). Many plants have been shown to foster biofilm growth along the root struc- tures for reasons ranging from protection against plant pathogens to nitrogen fixation (Gruyer et al., 2013). Wetland plants provide labile organic carbon exudates during active periods and by slow release during decay (García et al., 2010; Iannelli et al., 2011), and also contribute ions and oxygen (Hawkes et al., 2007). Given previ- ous extensive research on interactions between plants and bacteria (DanhornandFuqua, 2007) it is possible thatE. coliO157:H7attach- ment and persistencemay benefit from these interactions (Deering et al., 2012). Less attachment to C. utriculata roots suggests that this plant may have a lower potential for pathogen retention in TW. Interestingly, there was no statistical difference in pathogen attachment and persistence between the clean nylon surface and nylon with an established biofilm. It was hypothesized that com- petition with the natural consortia of bacteria in an environment would protect against invading E. coli O157:H7 cells since this is commonwith other organisms (Hibbing et al., 2010) andpathogens in aquatic biofilms (Nocker et al., 2013). The mechanisms in this system remain unknown and will warrant further research. The fateof theundetectedE. coliO157:H7cells andwhether they were sloughed into the effluent feed and still viable or inactivated by predation or competition was not determined. Other studies have indicated these factors play a large role in lowering human pathogenretention in theenvironment (Boutilier et al., 2009; Jasper et al., 2013; Li et al., 2014), however for strains with a low infection rate such as E. coliO157:H7, thesemaynot be enough to completely eradicate subsequent release and the potential infection. If E. coli O157:H7 is able to attach to roots, it may have prolonged ability to survive in aTW. Further research is needed for optimizingpathogen removal in TW. Acknowledgements Betsey Pitts at the Center for BiofilmEngineering provided assis- tance with microscopy and the Land Resources and Environmental Sciences and Plant Sciences departments provided assistance with the plant work. This research was supported by the Molecular Bio- sciencesProgramFellowship, theMontana InstituteonEcosystems, and the Center for Biofilm Engineering. Appendix A. Supplementary data Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.ecoleng.2016.10. 018. References Akratos, C.S., Tsihrintzis, V.A., 2007. Effect of temperature, HRT, vegetation and porous media on removal efficiency of pilot-scale horizontal subsurface flow constructed wetlands. Ecol. Eng. 29, 173–191. Austin, G., 2014. Green infrastructure for landscape planning. In: Routledge (Ed.), Integrating Human and Natural Systems. Routledge New York, New York. Barcina, I., Arana, I.S., 2009. The viable but nonculturable phenotype: a crossroads in the life-cycle of non-differentiating bacteria? Rev. Environ. Sci. Biotechnol. 8, 245–255. Baskar, G., Deeptha, V.T., Annadurai, R., 2014. Comparison of treatment performance between constructed wetlands with different substrates. Int. J. Res. Eng. Adv. Technol. 2 (3), 1–5. Behnke, S., Parker, A.E., Woodall, D., Camper, A.K., 2011. Comparing the chlorine disinfection of detached biofilm clusters with those of sessile biofilms and planktonic cells in single- and dual-species cultures. Appl. Environ. Microbiol. 77 (20), 7176–7184. Bitton, G., 2005. Wastewater Microbiology. John Wiley & Songs, Inc., Hokboken. Boutilier, L., Jamieson, R., Gordon, R., Lake, C., Hart, W., 2009. Adsorption, sedimentation, and inactivation of E. coli within wastewater treatment wetlands. Water Res. 43, 4370–4380. Carballeira, T., Ruiz, I., Soto, M., 2016. Effect of plants and surface loading rate on the treatment efficiency of shallow subsurface constructed wetlands. Ecol. Eng. 90, 203–214. Danhorn, T., Fuqua, C., 2007. Biofilm formation by plant-associated bacteria. Annu. Rev. Microbiol. 61, 401–422. Deering, A.J., Mauer, L.J., Pruitt, R.E., 2012. Internalization of E. coli O157:H7 and Salmonella spp. in plants: a review. Food Res. Int. 45, 567–575. Dinu, L.-D., Bach, S., 2011. Induction of viable but nonculturable Escherichia coli O157:H7 in the phyllosphere of lettuce: a food safety risk factor. Appl. Environ. Microbiol. 77, 8295–8302. Faulwetter, J.L., Gagnon, V., Sundberg, C., Chazarenc, F., Burr, M.D., Brisson, J., Camper, A.K., Stein, O.R., 2009. Microbial processes influencing performance of treatment wetlands: a review. Ecol. Eng. 35 (6), 987–1004. Fu, Z., Rogelj, S., Kieft, T.L., 2005. Rapid detection of Escherichia coli O157:H7 by immunomagnetic separation and real-time PCR. Int. J. Food Microbiol. 99, 47–57. García, J.D., Rousseau, P.L., Morato, J., Lesage, E.L.S., Matamoros, V., Bayona, J.M., 2010. Contaminant removal processes in subsurface-flow constructed wetlands: a review. Crit. Rev. Environ. Sci. Technol. 40 (7), 561–661. García, J.A., Paredes, D., Cubillos, J.A., 2013. Effect of plants and the combination of wetland treatment type systems on pathogen removal in tropical climate conditions. Ecol. Eng. 58, 57–62. Gruyer, N., Dorais, M., Zagury, G.J., Alsanius, B.W., 2013. Removal of plant pathogens from recycled greenhouse wastewater using constructed wetlands. Agric. Water Manage. 117, 153–158. Hawkes, C.V., DeAngelis, K.M., Firestone, M.K., 2007. Chapter 1 – root interactions with soil microbial communities and processes. In: Whitbeck, Z.G.C.L. (Ed.), The Rhizosphere. Academic Press, Burlington, pp. 1–29. Hench, K.R., Bissonnette, G.K., Sexstone, A.J., Coleman, J.G., Garbutt, K., Skousen, J.G., 2003. Fate of physical, chemical, and microbial contaminants in domestic wastewater following treatment by small constructed wetlands. Water Res. 37 (4), 921–927. Hibbing, M.E., Fuqua, C., Parsek, M.R., Peterson, S.B., 2010. Nat. Rev. Microbiol. 8, 15–25. Iannelli, R., Bianchi, V., Salvato, M., Borin, M., 2011. Modelling assessment of carbon supply by different macrophytes for nitrogen removal in pilot vegetated mesocosms? Int. J. Environ. Anal. Chem. 91 (7–8), 708–726. Jasper, J.T., Nguyen, M.T., Jones, Z.L., Ismail, N.S., Sedlak, D.L., Sharp, J.O., Luthy, R.G., Horne, A.J., Nelson, K.L., 2013. Unit process wetlands for removal of trace organic contaminants and pathogens from municipal wastewater effluents. Environ. Eng. Sci. 30 (8), 421–436. Klayman, B.J., Volden, P.A., Stewart, P.S., Camper, A.K., 2009. Escherichia coli O157:H7 requires colonizing partner to adhere and persist in a capillary flow cell. Environ. Sci. Technol. 43, 2105–2111. Kotti, I.P., Gikas, G.D., Tsihrintzis, V.A., 2010. Effect of operational and design parameters on removal efficiency of pilot-scale FWS constructed wetlands and comparison with HSF systems. Ecol. Eng. 36, 862–875. Li, F., Shi, M., Zheng, X., Zhang, N., Zheng, H., Gao, C., 2014. A novel method of rural sewage disinfection via root extracts of hydrophytes. Ecol. Eng. 64, 344–349. Licence, K., Oates, K.R., Synge, B.A., Reid, T.M.S.1, 2001. An outbreak of E.coli 0157 infection with evidence of spread from animals to man through contamination of a private water supply. Epidemiol. Infect. 126, 135–138. Liu, Y., Wang, C., Tyrrell, G., Hrudey, S.E., Li, X.-F., 2009. Induction of Escherichia coli O157:H7 into the viable but non-culturable state by chloraminated water and river water, and subsequent resuscitation. Environ. Microbiol. Rep. 1, 155–161. Nocker, A., Burr, M., Camper, A.K., 2013. Chapter 1. Pathogens in Water and Biofilms. In: Microbiology of Waterborne Disease, 2nd Edition. S. Percival, M.V., Yates, D., Williams, N. Gray and R. Chalmers, eds. Ongeng, D., Muyanja, C., Ryckeboer, J., Geeraerd, A.H., Springael, D., 2011. Development and validation of a culture-based method suitable for monitoring environmental survival of Escherichia coli O157:H7 and Salmonella enterica serovar Typhimurium in developing countries. Ann. Microbiol 61 (4), 809–817. Paruch, A.M., 2011. Long-term survival of Escherichia coli in lightweight aggregate filter media of constructed wastewater treatment wetlands. Water Sci. Technol. 63 (3), 558–564. Pyle, B.H., Broadaway, S.C., McFeters, G.A., 1995. A rapid, direct method for enumerating respiring enterohemorrhagic Escherichia coli O157:H7 in Water. Appl. Environ. Microbiol. 61, 2614–2619. Shaikh, N., Tarr, P.I., 2003. Escherichia coli O157:H7 shiga toxin-encoding bacteriophages: integrations, excisions, truncations, and evolutionary implications. J. Bacteriol. 185 (12), 3596–3605. Taylor, C.R., Hook, P.B., Stein, O.R., Zabinski, C.A., 2011. Seasonal effects of 19 plant species on COD removal in subsurface treatment wetland microcosms. Ecol. Eng. 37 (5), 703–710. Tonderski, K.S., 2009. Molecular and microbial advances in wetland science. Ecol. Eng. 35 (6), 959–960. USEPA, 2013. Basic Information about E. coli O157:H7 in Drinking Water, United States Environmental Protection Agency, http://water.epa.gov/drink/ contaminants/basicinformation/ecoli.cfm. VanKempen-Fryling, R.J., Stein, O.R., Camper, A.K., 2015. Presence and persistence of wastewater pathogen Escherichia coli O157:H7 in hydroponic reactors of treatment wetland species. Water Sci. Technol. 72 (1), 135–140. Wilson, S., Hamilton, M.A., Hamilton, G.C., Schumann, M.R., Stoodley, P., 2004. Statistical quantification of detachment rates and size distributions of cell clumps from wild-type (PAO1) and cell signaling mutant (JP1) Pseudomonas aeruginosa biofilms. Appl. Environ. Microbiol. 70 (10), 5847–5852. Zhang, D., Luo, J., Lee, Z.M.P., Maspolim, Y., Gersberg, R.M., Liu, Y., Tan, S.K., Ng, W.J., 2016. Characterization of bacterial communities in wetland mesocosms receiving pharmaceutical-enriched wastewater. Ecol. Eng. 90, 215–224.