Subscriber access provided by University of Otago Library
Environmental Processes
Extracellular DNA in Monochloraminated Drinking Water and its Influence on DNA-based Profiling of Microbial Community Bairoliya Sakcham, Amit Kumar, and Bin Cao Environ. Sci. Technol. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.estlett.9b00185 • Publication Date (Web): 21 Apr 2019 Downloaded from http://pubs.acs.org on April 22, 2019
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 19
Environmental Science & Technology Letters
1
Extracellular DNA in Monochloraminated Drinking Water and its Influence on DNA-
2
based Profiling of Microbial Community
3
Bairoliya Sakcham1,2, Amit Kumar1,3, Bin Cao1,3*
4 5
1
6
University, Singapore 637551
7
2
Interdisciplinary Graduate School, Nanyang Technological University, Singapore 637553
8
3
School of Civil and Environmental Engineering, Nanyang Technological University, Singapore
9
639798
Singapore Centre on Environmental Life Sciences Engineering, Nanyang Technological
10 11
*Corresponding author: Dr. Bin Cao, School of Civil and Environmental Engineering, Nanyang
12
Technological University, 50 Nanyang Ave, N1-01C-69, Singapore 639798; e-mail:
13
[email protected]; Tel: (+65) 6790 5277; Fax: (+65) 6791 0676.
14
1 ACS Paragon Plus Environment
Environmental Science & Technology Letters
Page 2 of 19
15
Abstract
16
The interaction between biofilms and disinfectant in drinking water distribution systems
17
(DWDS) as well as the role of this interaction in the formation of disinfection byproducts
18
(DBPs) has been extensively studied in recent years. In contrast, lysis of cells and/or release of
19
intracellular biomolecules from inactivated/damaged cells and their fate and implications are an
20
overlooked aspect in DWDS. In particular, DNA, once released into DWDS, may persist in
21
water as extracellular DNA (eDNA). In this study, we report for the first time that the total DNA
22
extracted from monochloraminated drinking water contains a high fraction of eDNA. Drinking
23
water samples were obtained from Locations 1 (~20-year-old pipeline) and 2 (~7-year-old
24
pipeline). At Location 1, 85-386 ng of eDNA was found per litre of sampled water, which
25
accounted for 52 ±12 % of total DNA, while at Location 2, 33-58 ng of eDNA was found per
26
litre of sampled water, accounting for 42±8% of the total DNA. We further showed that the
27
removal of eDNA reduced α-diversity, increased community evenness, and changed relative
28
abundance of detected taxa. Our findings lead to future research questions about the source, fate
29
and implications eDNA in DWDS.
2 ACS Paragon Plus Environment
Page 3 of 19
30
Environmental Science & Technology Letters
INTRODUCTION
31
Drinking water distribution systems (DWDS) house a huge diversity of microorganisms in the
32
form of planktonic cells in the bulk water and sessile microbial communities, i.e., biofilms, on
33
the inner pipe surfaces. Majority of the microorganisms (~90-98%) are found attached to the
34
inner surfaces of the distribution pipelines in the form of biofilms1. Biofilms are believed to
35
provide more favorable conditions for the sustenance of the microorganisms as they can protect
36
the residing organisms from the harsh conditions in the bulk water (e.g., low nutrient content,
37
dynamic flow patterns and/or presence of disinfectant residuals)2. Cells along with the biofilm
38
matrix and loose mineral deposits can make their way into the bulk water stream due to natural
39
dispersion of biofilms or changes in operational parameters (e.g., flow velocity, nutrient loading,
40
disinfection residual) causing biofilm detachment3.
41
To manage biofilms in DWDS, a secondary disinfectant monochloramine (MCA) is often used3.
42
MCA can inactivate cells in bulk water4-5, reducing biofilm formation. It can also effectively
43
penetrate biofilms and inactivate cells within the biofilms6-8. Various studies have found that the
44
interaction between the disinfectant and the organic matter in the form of DWDS biofilms
45
contributes greatly to the formation of disinfection by-products (DBPs), a group of emerging
46
contaminants potentially genotoxic and carcinogenic9. An overlooked aspect of the cell-
47
disinfectant interaction in the DWDS is the possible lysis of cells and release of intracellular
48
biomolecules into the drinking water. In particular, DNA, once released from the cells and
49
biofilms in DWDS, may persist in water as extracellular DNA (eDNA) because DNA reacts
50
slowly and moderately with MCA (kMCA =10-1 to 101 M-1 s-1) and chlorine (kchlorine =101 to 103
51
M-1 s-1), respectively10. DNA has been shown to be stable in various environmental conditions
52
and has been observed to persist and accumulate for long periods in freshwaters 11, marine
53
environments 12 and soils 13, engineered water systems14-15, especially in the presence of natural
54
organic matters 16. Although the occurrence and fate of eDNA in many natural environments
55
including soils and surface waters have been investigated, there have been no studies examining
56
the presence of eDNA in DWDS17. In addition, many studies characterized DWDS microbial
57
communities using DNA-based profiling methods18, however, eDNA has not been taken into
58
consideration. Whether and how eDNA influences DNA-based profiling of DWDS microbial
59
communities still remains a question. 3 ACS Paragon Plus Environment
Environmental Science & Technology Letters
60
In this study, we report for the first time that the total DNA extracted from monochloraminated
61
drinking water contains a high fraction of eDNA and the presence of eDNA obscures the DNA-
62
based estimation of microbial diversity and community profiling. Our results suggest that the
63
occurrence, fate and implications of eDNA, an overlooked component in drinking water, should
64
be carefully examined in future studies.
65
Page 4 of 19
MATERIALS AND METHODS
66
Sampling of Drinking Water. Drinking water samples were obtained from two different
67
locations, Location 1 (~20-year-old pipeline) and Location 2 (~7-year-old pipeline), both
68
distributing water treated with MCA as a secondary disinfectant. Glass fiber membrane (25 mm
69
in diameter, 0.4 µm nominal pore size, Grade GB-140; Advantec) packed in polypropylene in-
70
line filter holders (Sterlitech Corporation) were used to filter the drinking water directly from the
71
taps. The glass filter membranes have been shown to be capable of effectively retaining DNA19
72
and have been applied to extract environmental DNA from water samples20. The initial flow rate
73
for filtration ranged between 200-215 mL/min and a total of ~15 L of water was filtered through
74
each membrane. For each location, 4 samples were collected in parallel using a flow splitter from
75
the same tap (Figure S1 in the Supporting Information (SI)). The filter samples were then
76
subjected to various treatments and analyses (Figure S2) which are described below in details.
77
The SEM imaging methods and representative SEM images (Figure S3) can be found in SI.
78
Water quality parameters have been provided in Table S1.
79
DNase Treatment. After filtration, the membranes were taken out and divided into two equal
80
parts and one part was randomly assigned to DNase treatment (refer to SI for details). DNase
81
treatment has been previously shown to remove DNA only from membrane-compromised cells
82
without affecting the intact cells21. The samples were incubated with DNase (DNase-
83
Treated)/MilliQ water (Control) at 37 °C for 1 h and then washed twice with 0.1 M Phosphate
84
Buffered Saline (PBS) before proceeding for the DNA extraction. The effectiveness of the
85
DNase treatment was evaluated using a cell suspension (OD600 ~0.1) of Pseudomonas
86
aeruginosa (a bacterial strain isolated from DWDS) spiked with known amounts (2, 3 or 4 µg) of
87
genomic DNA.
4 ACS Paragon Plus Environment
Page 5 of 19
Environmental Science & Technology Letters
88
DNA Staining and Confocal Microscopy Imaging. A small portion (~2.5 × 2.5 mm) from the
89
center of each half was randomly assigned as a control or for DNase treatment. After DNase
90
treatment, samples were rinsed with sterile MilliQ water and stained using TOTO-1 (2 µM in
91
DMSO, excitation/emission spectrum of 488/488-640 nm) and counter stained with SYTO-60
92
(10 µM in DMSO, excitation/emission spectrum of 633/644-664 nm). The effectiveness of the
93
combination of dyes was evaluated using a cell suspension of P. aeruginosa imaged before and
94
after treatment with DNase to remove eDNA (refer to SI for details).
95
DNA Extraction and Quantification. DNA extraction was carried out using a modified
96
cetyltrimethylammonium bromide (CTAB) –Proteinase K method (refer to SI for details). DNA
97
was quantified using Quant-iT ™ dsDNA High-Sensitivity Assay Kit from Invitrogen. The
98
percentage of eDNA in the extracted samples was estimated based on the difference between the
99
amount of DNA extracted from the filter samples with and without DNase treatment.
100
qPCR, PCR Amplification and Sequencing. The abundance of eDNA was estimated using
101
qPCR based on the difference between the 16S rRNA gene copy numbers from the filter samples
102
with and without DNase treatment. Fungal contribution to eDNA was also assessed by
103
quantifying the 18S rRNA gene. PCR amplification was carried out in duplicate in a 20 µL
104
reaction volume using universal 16S rRNA gene primers. The amplicons were purified using
105
AMPure XP beads (Beckman Coulter, Brea, 138 CA) and quantified using Quant-iT™ dsDNA
106
High-Sensitivity Assay Kit. Following library preparation and quality control, pair ended
107
sequencing (300 × 300 bp) was carried out on an Illumina MiSeq platform. Experimental details
108
can be found in SI.
109
Sequencing Data Processing and Analysis. The sequences were processed using the DADA2
110
pipeline v1.8 in R v3.5.122. The resulting ASV (Amplicon Sequence Variants) table was filtered
111
for reads which were unassigned at Kingdom level or assigned as the Eukaryotic and Chloroplast
112
sequences. Only the prokaryotic sequences were kept for further analysis. ASV table is just a
113
higher resolution version of a traditional OTU table22. The ASV table obtained from the DADA2
114
pipeline was further analyzed in R (v3.5.1) using the “phyloseq”23 and “vegan”24 packages. The
115
resulting plots were made using the package “ggplot2”25. All analysis was conducted on data-sets
116
rarefied to the minimum number of sequences obtained for a sample. More details can be found 5 ACS Paragon Plus Environment
Environmental Science & Technology Letters
117
in the SI. Sequencing results have been deposited to the Sequence Read Archive (SRA) under
118
Biosample number SAMN10667470- SAMN10667485 and Bioproject number PRJNA512437.
119
Page 6 of 19
RESULTS AND DISCUSSION
120
Significant Fraction of eDNA in Total DNA Extracted from Drinking Water. Representative
121
CLSM images of the filter samples with and without DNase treatment stained with TOTO-1
122
(green fluorescence) and SYTO 60 (red fluorescence) are shown in Figure S4. TOTO-1 is
123
impermeable to intact cells and selectively stains eDNA and DNA in membrane-compromised
124
cells, while SYTO 60 stains intracellular DNA in intact cells as well26. Using a cell suspension
125
spiked with DNA, we have shown that these stains work appropriately in terms of cell
126
permeability/impermeability (Figure S5) and the DNase treatment effectively removed eDNA,
127
while it did not substantially reduce cell viability (>95% after the DNase treatment for all tested
128
conditions) (Figure S6). Treatment with DNase greatly reduced the TOTO-1 signal in the
129
images, suggesting a high abundance of eDNA in the samples. Further quantitative analyses
130
showed a 2-4-fold reduction in the biovolume of eDNA in the treated samples while the
131
biovolume for the intact cells remains unchanged (Figure S7A). In addition, greater proportion of
132
the green fluorescence signal colocalizes with the red fluorescence signal (Figure S7B),
133
indicating the presence of partially permeabilized cells, which is not surprising given the
134
monochloramine treatment of the collected water27. The used disinfectant monochloramine
135
might be a possible driver of DNA release in the DWDS (Figure S8).
136
Pairwise comparison was made for the DNA extracted from the control and the DNase-treated
137
halves from each filter sample and the amount of eDNA was estimated assuming that all eDNA
138
was removed by DNase treatment. The proportion of eDNA varied in different samples. At
139
Location 1, 85-386 ng of eDNA was found per liter of sampled water, which accounted for 52
140
±12 % of total DNA, while at Location 2, 33-58 ng of eDNA was found per liter of sampled
141
water, accounting for 42±8% of the total DNA (Figure 1A). These proportions for eDNA are
142
comparable with what have been reported in surface and lake water samples28. The proportion of
143
eDNA was also estimated using qPCR based on the total copy number of the 16S rRNA gene.
144
The qPCR-based estimation of eDNA facilitates the comparison of our data with relevant
145
literature in other systems such as soils and sediments13, 28 because most eDNA studies used 6 ACS Paragon Plus Environment
Page 7 of 19
Environmental Science & Technology Letters
146
marker genes such as the 16S rRNA gene to estimate the proportion of eDNA in the samples. It
147
also helps to confirm that by analysing the bacterial community through the 16S rRNA gene
148
amplicon sequencing approach, we are covering the contribution of a major proportion of the
149
eDNA. The qPCR results showed an eDNA fraction of 48±12% and 29±11% at Locations 1 and
150
2, respectively (Figure 1B). Single factor Analysis of Variance (ANOVA) showed no significant
151
difference between the estimated proportion of eDNA obtained through extraction and that
152
through 16S rRNA gene-based qPCR (p = 0.66 and 0.13 for Location 1 and 2, respectively),
153
suggesting that majority of the eDNA in the system is of bacterial and archaeal origin and is
154
expected to contribute to the analysis of the microbial community through the 16S rRNA gene
155
amplicon sequencing approach. For verification, fungal contribution to eDNA was also assessed
156
by quantifying the 18S rRNA gene (Table S2). Our results show that in 1 L water sample, a copy
157
number of about 103 – 104 were detected for the 18S rRNA gene, which is substantially lower
158
than the 16S rRNA genes (~108). Thus, taking 18S rRNA gene copy numbers into account does
159
not change the calculated proportion of eDNA based on the quantitative results using qPCR.
160
7 ACS Paragon Plus Environment
Environmental Science & Technology Letters
Page 8 of 19
161 162
Figure 1. Percentage of eDNA found in the samples based on (A) DNA extraction (n = 4) and
163
(B) total 16S rRNA gene copies quantified in the control and treated halves of the filters (n = 4).
164
Locations 1 and 2 refer to systems which are ~20- and ~7-year-old, respectively. S1, S2, S3 and
165
S4 are samples from different filters (annotated by different colors) collected in parallel from one
166
tap at each location. 8 ACS Paragon Plus Environment
Page 9 of 19
Environmental Science & Technology Letters
167
Influence of eDNA on Estimation of α- and β-diversity. The α-diversity indices for samples
168
with and without DNase treatment are shown in Figure 2A1 and 2B1 for Locations 1 and 2,
169
respectively. The Chao1 index and Abundance-based Coverage Estimator (ACE) are estimators
170
of richness, taking into consideration rare organisms that may have been missed due to under
171
sampling29, while the observed richness is the actual number of ASVs detected in the samples.
172
Both Chao1 and ACE have values comparable with the observed richness (Location 1, p = 1.0,
173
0.99; Location 2, p = 0.99, 0.99 for Chao1 and ACE, respectively), confirming that the sampling
174
was adequate. The DNase treatment, i.e., removal of eDNA, significantly reduced the Chao1 and
175
ACE indices for samples from Location 2 (p 0” indicates more detection in the
240
control samples; “FC