Impact of UV irradiation at full scale on bacterial communities in drinking water

Water in a full-scale drinking water treatment plant was irradiated with ultraviolet (UV) doses of 250, 400, and 600 J/m2, and the effect on bacterial communities investigated using 16s rRNA gene amplicon sequencing, heterotrophic plate counts (HPCs), coliform, and Escherichia coli counts. The bacteria in the irradiated water were also analyzed following storage for 6 days at 7 °C, to approximate the conditions in the distribution system. The log10 reduction of HPCs at 400 J/m2 was 0.43 ± 0.12. Phylogenetic examination, including DESeq2 analysis, showed that Actinobacteria was more resistant to UV irradiation, whereas Bacteroidetes was sensitive to UV. Phylum Proteobacteria contained monophyletic groups that were either sensitive or resistant to UV exposure. The amplicon sequence variants (ASVs) resistant to UV irradiation had a greater average GC content than the ASVs sensitive to UV, at 55% ± 1.7 (n = 19) and 49% ± 2.5 (n = 16), respectively. Families Chitinophagaceae, Pelagibacteraceae, Holophagaceae, Methylophilaceae, and Cytophagaceae decreased linearly in relative abundance, with increasing UV dose (P < 0.05, Pearson’s correlation). When irradiated water was stored, Chitinophagaceae, Comamonadaceae, and Flavobacteriaceae families decreased in relative abundance, whereas ACK-M1, Mycobacteriaceae, and Nitrosomonadaceae were increasing in relative abundance. This suggests that the impact of UV irradiation cannot only be considered directly after application but that this treatment step likely continues to influence microbial dynamics throughout the distribution system.


INTRODUCTION
Ultraviolet (UV) irradiation is widely used as a disinfection method for drinking water treatment. The technique became increasingly popular in the 1990s when its ability to disinfect water containing Cryptosporidium and Giardia was recognized 1 . Unlike other disinfection methods such as chlorination or ozonation, UV irradiation requires no addition of chemicals and low-pressure UV produces insignificant amounts of disinfection byproducts 1,2 . Various bacteria have different UV susceptibility: a 4 − log 10 reduction of a lab-grown environmental isolate of Mycobacterium avium requires a dose of 128 J/m 2 UV 254 nm 3 , whereas for the same reduction, cultivated environmental isolate of Escherichia coli requires a dose of 81 J/m 2 4 . The dose required for disinfection can also be affected by suspended particles in the water, which can absorb and scatter UV light and affect UV efficiency 5,6 .
The disinfection mechanism resulting from exposure to UV is mainly damage to nucleic acids by irradiation 7 . Nucleotides absorb UV light with wavelengths of between 200 and 300 nm with a peak absorption between 260 and 265 nm 8 . The absorption of light triggers the formation of mutagenic DNA lesions, such as cyclobutane pyrimidine dimers and 6-4 photoproducts 9 . Both pyrimidines and purines can absorb UV light, although pyrimidines are considered to be more photoreactive 10,11 . When nucleotides are damaged by UV light, the DNA replication is blocked, resulting in cell inactivation 12,13 . Some microorganisms are able to repair UV damage either by photoreactivation or dark repair 9,14 .
The impact of UV on target microorganisms has largely been studied using cultivation-based techniques of monocultures at laboratory scale 1 . At full scale, a biodosimetry test is used to calibrate the irradiation dose for UV reactors by spiking a known concentration of a specific cultured microorganism and calculating the log-reduction. This is compared with results from a calibrated laboratory UV reactor to calculate the final UV dose of the full-scale UV reactor at a specific UV transmission and flow 15 . The validity of these tests to assess disinfection of drinking water, however, is debatable, as the majority of microorganisms in drinking water cannot currently be cultivated 16 , bacteria in drinking water are diverse, and bacteria in the environment have an increased UV resistance compared with laboratory-cultivated strains 1 . Exposure to UV can also cause some bacteria to enter a viable but nonculturable state as a response to environmental stress 17,18 .
Molecular DNA-based methods analyze the microbial community without the need for cultivation and, as UV irradiation causes DNA lesions and reduces the number of amplifiable target templates in the PCR reaction 13,19 , amplicon-based methods can describe which types of bacteria and genes are affected by UV [20][21][22] . Microbial inactivation of Pseudomonas aeruginosa and Enterococcus faecium by UV was assessed by cultivation and quantitative PCR (qPCR) 19 and impact of UV on adenovirus concentrations were measured by cell culture infectivity and long-range PCR with subsequent qPCR 23 . The impact of UV irradiation on the number of bacteria in drinking water has been quantified with 16s rRNA gene amplification 24 .
Although the value of UV irradiation for reduction of microbial pathogens in drinking water is not disputed, it is not known how UV irradiation impacts other bacteria that are undesirable. Some bacteria can cause water quality issues, by producing unpleasant odors and tastes 25 , or exacerbate corrosion of infrastructure 26 . Descriptions of the bacterial community in water using 16s rRNA gene amplicon sequencing have described changes in drinking water during distribution 27 ; however, few studies exist, examining the impact of UV exposure on the bacterial community in drinking water. An initial study using 16s rRNA gene amplicon sequencing conducted by Nocker et al. 21 showed how some phyla of bacteria were affected by a single dose of UV in a full-scale drinking water treatment plant. The approach of using 16S rRNA gene amplicon sequencing has been used to observe changes in the bacterial community for UV wastewater disinfection 22 and UV irradiation of marine water 28 .
The current study is a detailed examination of the impact of three different UV irradiation doses (250, 400, and 600 J/m 2 ) in drinking water at full scale. The contents and diversity of the bacterial community was assessed using 16S rRNA gene amplicon sequencing. The results of this molecular analysis were compared with traditional cultivation-based methods determining heterotrophic plate counts (HPCs), coliforms, and E. coli. To understand how changes might occur during distribution of water irradiated with UV, bacterial communities of irradiated water, stored at temperatures simulating those of a distribution system, were also investigated.

RESULTS
Effect of UV irradiation on bacterial community structure The bacterial community in water irradiated with different UV doses (250, 400, and 600 J/m 2 ) and following 6 days water storage at 7°C was investigated using 16S rRNA gene amplicon sequencing (Fig. 1). Canonical correspondance analysis (CCA) of the relative abundances of amplicon sequence variants (ASVs), which takes into account the parameters of UV dose and storage, showed increasingly dissimilar bacterial communities with higher UV dose (P < 0.001, CCA followed by analysis of variance (ANOVA)). This was also observed with ordination analysis using the Bray-Curtis distance between samples ( Supplementary Fig. 1). Although storage of the water affected the community composition of the UV-irradiated samples (Supplementary Fig. 2A, P < 0.001, R 2 = 0.62, Analysis of Dissimilarities (ADONIS), storage had no impact on the bacterial community that had not been irradiated by UV according to Bray-Curtis analysis ( Supplementary  Fig. 2B, P = 0.36, R 2 = 0.10, ADONIS).
To further investigate how irradiation shaped the bacterial community structures, alpha diversity was assessed by comparing evenness (Pielou's measure), species richness (number of observed ASVs), and diversity (Shannon index) ( Supplementary  Fig. 3). When all water samples exposed to UV irradiation (250, 400, and 600 J/m 2 ) were grouped, the evenness of the community increased compared with the non-irradiated water from 0.49 ± 0.016 to 0.54 ± 0.025 (P < 0.001, one-way ANOVA). Species richness in water samples ranged from 109 to 140, but no changes in species richness were observed by UV irradiation (P > 0.05, oneway ANOVA). Changes in Shannon diversity reflected those in evenness, with UV exposure resulting in increased Shannon diversity in the bacterial communities, from 2.33 ± 0.094 for nonirradiated communities to 2.61 ± 0.12 when all water samples exposed to UV were grouped for analysis (P < 0.001, one-way ANOVA).
The effects of UV irradiation were observed as changes in the relative abundance of specific bacterial taxa as the UV dose increased (Fig. 2). Following agglomeration of ASVs into phylum level, Bacteroidetes decreased in relative abundance (P < 0.05, Pearson's correlation), whereas Actinobacteria increased in relative abundance (P < 0.05, Pearson's correlation). Five families showed linear decreases in relative abundance (P < 0.05, Pearson's correlation), including Pelagibacteraceae (R = −0. Phylogenetic relationships of bacteria impacted by UV ASVs from bacteria affected by UV irradiation were identified by differential abundance analysis using DESeq2 where negative and positive log 2 fold changes were defined as sensitive-and resistant to UV, respectively. Thirty-five out of 164 ASVs with a significant change in differential abundance (P adjusted < 0.05) were identified (Supplementary Figs 4 and 5). Specific clades were identified by their sensitivity to UV (Fig. 3), phylum Bacteroidetes was sensitive to UV, with ASVs identified as Sediminibacterium, Sphingobacteriaceae, and Cytophagaceae having a negative log 2 fold change. Phylum Actinobacteria was resistant to UV with ASVs classified as families ACK-M1 and C111 with a positive log 2 fold change. The Proteobacteria phylum included ASVs that were both sensitive and resistant to UV: ASVs classified as Pelagibacteraceae (order Ricketsiales, Alphaproteobacteria), Methylophilaceae (Betaproteobacteria), and Limnohabitans (Betaproteobacteria) showed a negative log 2 fold change, whereas others including Rhodospirillaceae (Alphaproteobacteria), Ralstonia (Betaproteobacteria), Polynucleobacter (Betaproteobacteria), and Rhodoferax (Betaproteobacetria) demonstrated positive log 2 fold changes. Phylogeny was determined to be a cause for UV sensitivity when the community sensitive to UV, which consisted of 16 ASVs showed a difference in composition compared with all identified ASVs in the samples at phylum level (P = 0.020, Fisher's exact test). The GC content of the 16S rRNA gene region of ASVs identified by DESeq2 analysis was calculated (Supplementary Table 2) where ASVs with positive log 2 fold change (n = 19) had a significantly greater average GC content compared with ASVs with negative log 2 fold change (n = 16), 55% ± 1.7 and 49% ± 2.5, respectively (P < 0.001, one-way ANOVA).

Fig. 1
Canonical correspondence analysis (CCA) of bacterial communities and UV dose. UV dose is used as a numeric variable and sample type as a factor. Bacterial communities were analyzed from water exposed to UV (0, 250, 400, and 600 J/m 2 ; blue, green, orange, and red circles, respectively) or after storage (6 day storage at 7°C, stars); n = 5 for UV dose of 0 J/m 2 and n = 6 for storage not exposed to UV; n = 3 for all UV doses >0 and sample type.
Impact of water storage on bacterial communities exposed to UV irradiation As water exposed by UV irradiation is often distributed to the consumer over a period of days and at low temperatures, changes that could occur in the bacterial community following irradiation and distribution were assessed (Fig. 1). To replicate these conditions, water samples were stored for 6 days at 7°C following UV exposure. The alpha diversity measure evenness recovered after storage to similar levels as before UV irradiation when the UV-irradiated and -stored samples were grouped together and compared with the non-irradiated samples, 0.49 ± 0.026 (mean ± SD) and 0.48 ± 0.034, respectively ( Supplementary Fig. 3). The storage did not affect the species richness (P > 0.05, one-way ANOVA). As with evenness, the Shannon diversity also recovered after storage for the grouped irradiated and stored samples to the same levels as the non-irradiated samples, Shannon index 2.36 ± 0.14 and 2.29 ± 0.16, respectively.
To assess which bacteria were able to repair the UV damage or regrow faster than other members in the community, the relative abundance of ASVs representing bacterial taxa in the UVirradiated water were compared before and after storage (Fig. 4). Following UV irradiation, the Pelagibacteraceae family comprised 37-50% of the community ( Fig. 2) but after storage this increased to 46-59% ( Fig. 4, P < 0.05, one-way ANOVA). The families Flavobacteriaceae, Trebouxiophyceae, C111, and Hyphomicrobiaceae were not initially affected or increased in relative abundance by UV irradiation (Fig. 2 and Supplementary Fig. 9), but decreased in relative abundance after storage ( Fig. 4, P < 0.05, one-way ANOVA). A relative decrease in the Trebouxiophyceae family and changes in relative abundance for 6 (out of 34) additional families (P < 0.05, one-way ANOVA) were observed when the water that had not been UV-irradiated was stored ( Supplementary Fig. 16).
To assess the impact of UV exposure on the bacterial community in distributed water, the relative abundance of bacterial taxa in stored water that had, or had not been irradiated with UV were compared (Fig. 5). Five families of bacteria decreased in relative abundance when water exposed to UV was stored, including Flavobacteriaceae, Chitinophagaceae, Cytophagaceae, Comamonadaceae, and Sphingobacteriaceae (one-way ANOVA, P < 0.05), and 11 families increased in relative abundance following UV exposure and storage, including ACK-M1, Oxalobacteraceae, Nitrosomonadaceae, Rhodospirillaceae, mb2424, Sinobacteraceae, LD19 Trebouxiphyceae, Gemmatimonadaceae, Mycobacteriaceae, and Acetobacteraceae (P < 0.05, one-way ANOVA).

Impact of UV on conventional microbial indicators
HPCs were assessed for water irradiated with 250, 400, and 600 J/ m 2 , and after storage. HPCs in the UV-exposed water decreased to an average of 1.2 ± 0.83 CFU/mL from 4.5 ± 1.4 CFU/mL (Fig. 6, P < 0.001, one-way ANOVA). There was no significant change in HPCs between the three different doses (P > 0.05, one-way ANOVA, followed by a Tukey's test): log 10 reductions in HPCs were 0.53 ± 0.12 and 0.43 ± 0.12 for 250 J/m 2 and 400 J/m 2 , respectively. Although the number of heterotrophic bacteria increased during storage in the non-irradiated samples, from 4.5 ± 1.4 to 8.7 ± 6.2 CFU/mL, this was not statistically significant (P > 0.05, one-way ANOVA). HPCs in water that had been UV-irradiated did not change during storage, regardless of dose (P > 0.05, one-way ANOVA, followed by a Tukey's test). No differences in fast growing HPCs after 3 days were seen for any of the samples (P > 0.05, oneway ANOVA, followed by a Tukey's test; Supplementary Fig. 17). To assess coliform counts, 60 L of water was concentrated using dead-end hollow-fiber ultrafiltration 29 following exposure of the water to 0, 400, and 600 J/m 2 (n = 3 for each dose). Coliforms (1.7 ± 0.57 coliforms/10 L) and E. coli (0.97 ± 0.93 E. coli/10 L) were detected in the concentrates of the untreated water but not in any sample of water exposed to UV (Supplementary Table 1).

DISCUSSION
As bacteria in drinking water are diverse and likely have different susceptibility to UV 1 , this study examined bacterial communities in water irradiated with different UV doses, using 16S rRNA gene amplicon sequencing. The impact of storage on the relative abundance of different taxa was also determined to describe how exposing water to UV irradiation could alter the microbiology of the distributed water. The community descriptions were compared with conventional methods using indicator bacteria to assess and compare how these approaches reflect the impact of UV irradiation on microbial water quality.
Examining the relationship between phylogeny and UV exposure revealed bacterial clades that were both sensitive and resistant to UV. Actinobacteria phylum was resistant, Bacteroidetes phylum showed sensitivity to UV, whereas some monophyletic groups in the Proteobacteria phylum showed sensitivity and others resistance to UV (Fig. 3). Low and high GC content calculated from the different ASV DNA sequences correlated with how the taxa were categorized as sensitive or resistant to UV, respectively (Supplementary Table 2). It has been proposed that the high genomic GC content in Actinobacteria contributes to their UV resistance 30 , perhaps as TT and TC nucleotide sites are more photoreactive 31 . Reichenberger et al. 32 proposed that phyla could be classified into GC-rich (Actinobacteria), GC-intermediate (Proteobacteria), and GC-poor (Bacteroidetes), linking UV sensitivity to a genomic, and thus heritable, trait. Nocker et al. 21 showed resistance of Actinobacteria to UV in drinking water and suggested that the Proteobacteria phylum is sensitive to UV. In addition to GC content, Actinobacteria, Bacteroidetes, and Proteobacteria have distinct membrane structures that can result in these bacteria having different degrees of resistance to UV 33 . A study from McKinney and Pruden 20 observed that two Gram-positive organisms (Staphylococcus aureus and E. faecium, both Firmicutes) were more resistant to UV irradiation than two Gram-negative bacteria (E. coli and P. aeruginosa, phylum Proteobacteria and Firmicutes, respectively), which could be due to the thicker peptidoglycan layer in Gram-positive organisms protecting the cells from UV damage. Actinobacteria are Gram-positive and have been isolated in environments exposed to strong solar UV radiation 30,34 . In addition, Actinobacteria are spore forming 35 and this can also contribute to their resistance to UV 36,37 . In contrast the phyla Bacteroidetes are Gram-negative rods that do not form endospores 38 and this combination could be a reason for the sensitivity to UV demonstrated by the members of the Bacteroidetes phylum observed in the current study. Strains classified in the Bacteroidetes phylum are abundant in high-altitude lakes receiving strong UV radiation from the sun 39 , suggesting tolerance to UV exposure; however, the Bacteroidetes phylum is diverse and the members of the Bacteroidetes in this study may be UV-sensitive, as the source water used in the treatment plant in this study is from a low-altitude lake. Proteobacteria are Gram-negative and the division is the largest and most diverse among prokaryotes 40 . In this study within the Proteobacteria, the response to UV seemed largely related to GC content. In Alphaproteobacteria, the Rickettsiales order was sensitive to UV (low GC content, ranging from~30% to 40% 41 ), whereas the Rhodospirillaceae family (high GC content of~65% 41 ) was resistant to UV, and in Betaproteobacteria, the Methylophilaceae family showed sensitivity to UV (low GC content of~35-40% 42 ), whereas Ralstonia and Rhodoferax were resistant to UV (high GC content of~67% 43 and~60% 44 , respectively). The exceptions were the sensitivity of Polynucleobacter, with a GC content of~45% 45 to UV, and that members of Limnohabitans (GC content of 59% 44 ) were both resistant and sensitive. These observations could be due, in part, to the few isolates of these bacteria that have been sequenced, to determine GC content, which may differ from those present in this study.
Depending on the time for drinking water to reach consumers, some bacteria in the water could repair UV damage 46 received at the treatment plant. In this study, the evenness and diversity of the bacterial community were affected by UV, although both metrics indicated that the community rebounded following storage. This could largely be attributed to changes in the relative abundance of Pelagibacteraceae, which was highly abundant in the community. The relative abundance of this family recovered after UV exposure and storage (Fig. 4), possibly due to the ability of these bacteria to repair UV damage and/or by surviving the UV exposure, and then regrowing faster than the other surviving taxa. As Pelagibacteraceae, or the SAR11 clade, are abundant in drinking water 47 and marine water 48,49 ; they may have acquired UV resistance to survive the UV exposure from the sun in their preferred niche and, with a large surface-to-volume ratio, efficient nutrient uptake in oligotrophic environments 50 allows them to grow faster than other members of the bacterial community 51 .
The relative abundance of families Flavobacteriaceae, Trebouxiophyceae, C111, and Hyphomicrobiaceae was initially not affected or increased by UV irradiation (Supplementary Fig. 14) but during storage decreased in relative abundance (Fig. 4). The initial relative UV resistance of these taxa could be due to the 16S rRNA gene not Fig. 3 Phylogenetic tree of ASVs identified as sensitive and resistant to UV. The inferred phylogenetic tree with branches colored by sensitivity and resistant to UV was created with ASVs recognized as having a signficant change with differential abundance analysis using DESeq2 (P adjusted < 0.05). Samples not exposed to UV (n = 5) were compared to UV irradiated samples (n = 9), with the dose as a numeric parameter. A negative log2fold change was considered a UV-sensitive taxa (red) and positive log2fold change as a UV-resistant taxa (blue). The scale bar shows an estimate of the substitutions per nucleotide position. See Supplementary Fig. 4 and Supplementary Table 2 for output data from DESeq2. See Supplementary Fig. 5 depicting the ASVs classified as sensitive and resistant to UV among all identified ASVs. receiving any DNA lesions, but damage by UV irradiation can also cause cell death via by, e.g., creating protein-DNA crosslinks 52 , which would not be detected with 16S rRNA gene amplicon sequencing. In addition, although valuable in many applications, the use of viability dyes (e.g., propidium monoazide (PMA) and propidium iodide (PI)) together with qPCR or flow cytometry (Supplementary Table 5) to assess UV disinfection is difficult, as UV irradiation damages DNA and not cell membranes.
UV irradiation is increasingly considered as a disinfection process for drinking water, and with specific observations of how the bacterial community would change during distribution from a treatment plant applying UV, it was also important to understand how this community compared with distributed water that had not been UV-irradiated water. After storage, some families, including Chitinophagaceae, Comamonadaceae, and Flavobacteriaceae decreased in relative abundance in water that had been UV-irradiated compared with untreated, stored water, whereas other families increased in relative abundance, including ACK-M1, Mycobacteriaceae, and Nitrosomonadaceae (Fig. 5). The ACK-M1 family is included in the Actinomycetales order and the presence of Actinomycetes has been associated with taste and odor problems, due to their production of odorants, including 2methylisoborneol (2-MIB) 25 . 2-MIB is an odorant that is difficult to remove with conventional water treatment processes; however, Chitinophagaceae bacterium, among other bacteria isolated from an activated carbon filter, are able to biodegrade 2-MIB 53 . This suggests that UV treatment may both increase the abundance of bacteria associated with production of odorants and decrease those with the potential to biodegrade the problematic compounds. This may be a concern for drinking water treatment plants, as preventing taste and odor problems in drinking water is of great importance 25 .
In addition to aesthetic compounds, the presence of opportunistic pathogens such as Legionellae, Mycobacteria, and P. aeruginosa are a concern in distributed water and can also persist in biofilms 54 . In this study, family Mycobacteriaceae (classified as Mycobacterium; Supplementary Fig. 15) was detected in higher relative abundance in water after UV exposure and storage, than in untreated, stored water (Fig. 5). Shin et al. 55 have described M. avium as more resistant to UV irradiation than other waterborne pathogens, and as this species has also shown resistance to other disinfection methods such as chlorine, monochloramine, chlorine dioxide, and ozone 56 , it could be of concern for drinking water producers.
Not much is known about how a UV installation will affect the biofilm in the drinking water distribution system (DWDS) in the long term 57 . The biofilm community in DWDS is complex and can differ depending on location and material 58,59 . The genus Nitrosomonas, included in the Nitrosomonadaceae family, which increased in relative abundance by UV irradiation and storage (Fig. 5), has been found in DWDS biofilms 60,61 . Comamonadaceae has been observed in DWDS biofilms 62 and Flavobacterium, Fig. 4 Comparison of bacterial taxa at family level, in water samples exposed to UV, with and without storage. Relative abundance of bacterial taxa in water samples exposed to UV, and then stored for 6 days at 7°C (red), or not stored (gold), were compared. Only families with significant difference in relative abundance between groups (P < 0.05, one-way ANOVA) are shown, see Supplementary Figs 11-15 for full data set. Families are ordered with the most relative change in the top left with subsequent decrease. The boxes show the interquartile range, the line inside each box represents the median, whiskers show maximum and minimum values, and the crosses show outliers. n = 9 for water samples exposed to UV and n = 9 for water exposed to UV and incubated. Samples irradiated at different UV doses were grouped for this analysis.
included in the Flavobacteriaceae family, has been observed in drinking water wells 63 , both decreasing in relative abundance by UV irradiation and storage in this study. If a treatment plant using UV disinfection selects for certain bacteria, this will inevitably shape the DWDS biofilm, which will affect the water quality in the end.
Despite the observation in this, and other studies that bacteria are affected differently by UV irradiation, biodosimetry tests are based on single target organisms. This makes it difficult to assess UV irradiation in full scale 21,24,64 . However, in this study families Chitinophagaceae, Pelagibacteraceae, Holophagaceae, Methylophilaceae, and Cytophagaceae showed sensitivity to UV in a linear correlation with dose (Fig. 2). These bacterial families have been widely observed in drinking water, freshwater, and river water [65][66][67][68] , and the 16S rRNA gene from these families could be a feasible target biomarker for evaluating UV irradiation dose at full scale. As water treatment plants have their own unique microbiome 69 , 16S rRNA gene amplicon sequencing to screen for these families, and potentially other biomarker target organisms, could then support targeted gene analysis, e.g., qPCR, to be applied routinely to assess UV disinfection. Longer fragment length of the qPCR target (compared with that used in this sequencing-based study) could increase the resolution of the methodology, as more UV damage can be observed with qPCR of longer amplicons 19 .
The need for molecular biology techniques was supported by observations in this study and others 70 of HPCs, where log reductions attributed to UV irradiation were lower than previously determined with lab-cultivated organisms 1 . Several taxa that have been identified from drinking water by growth on HPCs 71-73 were among the taxa identified by sequencing in this study as UVsensitive, including Methylophilaceae, Sphingobacteriales, and Comamonadaceae correlating with the reduction in HPCs by UV irradiation (Fig. 6); however, no change in HPCs was observed when comparing the UV-irradiated water before and after storage, which was seen by molecular techniques (Fig. 1). Mofidi and Linden 74 showed regrowth of HPCs following UV disinfection at 200, 600, and 1400 J/m 2 , and 7 days of storage at 20°C to the same counts as the non-irradiated control, with regrowth attributed to the robustness of the natural bacterial community in drinking water, and likely temperatures, as this was not observed in the current study, where a lower temperature was used for storage. Thus, HPCs are particularly unsuitable for monitoring quality of UV-irradiated water in systems where water temperatures are low, with overall low numbers making logreduction comparison difficult or meaningless. This was also the case for coliforms and E. coli, which were completely removed by UV irradiation, as expected 4,75 .
Although UV treatment is widely used in full-scale drinking water treatment plants, it is difficult to apply traditional Fig. 5 Comparison of bacterial taxa at family level following storage of water. Relative abundance of bacterial taxa in water samples that were stored for 6 days at 7°C, and were either not exposed (blue) or exposed (red) to UV. Only families with significant difference in relative abundance between groups (P < 0.05, one-way ANOVA) are shown; see Supplementary Figs 11-15 for full data set. Families are ordered with the most relative change in the top left with subsequent decrease. The boxes show the interquartile range, the line inside each box represents the median, whiskers show maximum and minimum values, and the crosses show outliers. n = 6 for no UV samples and n = 9 for UV-irradiated samples. Samples irradiated at different UV doses were grouped for this analysis.
microbiological methods to assess dose and quality of the distributed water. A variety of treatment processes produce water containing microorganisms with different UV susceptibility and different wavelengths of UV can be more efficient on certain organisms 76 . In addition to treatment plants, additional applications of UV disinfection include point-of-use and within the pipe network to prevent regrowth of opportunistic pathogens, including UV LED devices 57 . The need to assess the impact of these technologies, both for process monitoring purposes and to investigate the bacterial community that is exposed to consumers through water, could use a molecular biology approach, based on biomarkers for specific taxa. In addition, as the 16S rRNA gene is a fraction of the complete genome in bacteria, shotgun metagenomics examining the genome of UV-irradiated bacteria could investigate the role of GC content on survival, the sensitivity of specific genes involved in biofilm formation, DNA repair, antibiotic resistance, and virulence to understand how UV irradiation affects the microbiota and quality of drinking water. As the DWDS biofilm also interacts with drinking water quality, the impact of UVexposed water on biofilm formation and composition, in lab-or full scale, would also be of interest.

Sampling of water
Drinking water was sampled at Görvälnverket waterworks in Stockholm, Sweden. The waterworks operated by Norrvatten AB produces 1.6 m 3 /s drinking water from the surface water taken from Lake Mälaren. Water is produced using flocculation with aluminum sulfate, followed by sedimentation, rapid sand filtration, granulated activated carbon, UV irradiation, addition of monochloramine, and pH adjustment. The UV aggregate is a low-pressure UV (Trojan UV) emitting light at a wavelength of 254 nm with the USEPA standard 8 , routinely at a UV dose of 400 J/m 2 . The aggregate is 1000 mm in diameter consisting of ten rows of four UV lamps diagonally placed through the pipe. The UV intensities are measured with five light intensity meters evenly distributed throughout the aggregate measuring UV intensity at 200-300 nm. When sampling was conducted, the flow rate through the UV aggregate was~500 L/s and UV transmittance~84%. Water exposed to three different doses of UV irradiation (250, 400, and 600 J/m 2 ) and water that had not been irradiated (0 J/m 2 ) were sampled both directly after exposure and following storage for 6 days at 7°C. The conditions for storage were chosen to resemble the water temperature (~1-8°C) in spring (February-April) and residence time (maximum 14 days) in this distribution system. General water quality data can be found in Supplementary Table 4. Sampling was done in March on two successive days with sampling for 0 J/m 2 (n = 3 for all analyses except n = 2 for DNA analysis) and 250 J/m 2 (n = 3) on day 1 and 0, 400, and 600 J/m 2 (n = 3, for each dose) on day 2. Each sample was collected in 1 L sterile borosilicate bottles for DNA analysis, 500 mL sterile borosilicate bottles for HPCs, and 6 × 10 L sterile plastic cans for dead-end ultrafiltration and Colilert analysis. For the storage experiment, water was collected in 2 L sterile borosilicate bottles for DNA and HPC analyses, wrapped in aluminum foil for 6 days at 7°C (0 J/m 2 n = 6 and each UV dose n = 3). For DNA analysis, 1 L water samples were filtered onto 0.22 µm filters (Merck, Germany) and filters were stored at −20°C until DNA extraction. Dead-end ultrafiltration 29 was done using 60 L of water exposed to UV irradiation doses of 0, 400, and 600 J/m 2 . Rexeed 25AX filters (Asahi Kasei Corporation, Tokyo, Japan) were pretreated with fetal calf serum (PAA Laboratories, Waltham, MA, USA) to prevent adhesion of microorganisms. After filtration, the concentrate was eluted by back-flushing with 500 mL elution buffer (phosphate-buffered saline containing 1% Tween 80 and 0.01% Antifoam A (both from Sigma-Aldrich, St. Louis, MO, USA). The final volume typically consisted of 600-700 mL.

Conventional water quality parameters
HPCs were determined by the SS-EN ISO 6222:1999 standard with the pour plate method using 1 mL of water and addition of melted yeast peptone agar, followed by incubation at 22°C for 3 (72 h) and 7 days (168 h) in triplicate. Concentrations of coliforms and E. coli were assessed using 100 mL of the dead-end ultrafiltered eluate water with the Colilert®−18 (IDEXX Laboratories Sverige AB) according to the manufacturer's instructions.
Bacterial community analysis DNA was extracted from filter papers using the FastDNA Spin Kit for Soil (MP Biomedicals). Filters were cut into strips and added to Lysing Matrix E tubes with pre-added sodium phosphate, then extracted following the manufacturer's instructions. Empty filters were extracted for use as negative controls. Extracted DNA was stored at −20°C. The V3-V4 region of the 16S rRNA gene was amplified using primers 341 F (5′-CCTACGGGNGGCWGCAG-3′) and 785 R (5′-GACTACHVGGGTATCTAATCC-3′) 77 . The PCR reaction (25 µL) contained 12.6 µL MilliQ-water, 10 µL 5 Prime Hot MasterMix (Quantabio, USA), 0.4 µL (20 mg mL −1 ) bovine serum albumin, 0.5 µL (10 µM) forward and reverse primers, and 1 µL template DNA. The PCR cycling settings were 94°C for 3 min and 35 cycles of 94°C for 45 s, 50°C for 1 min, 72°C for 1.5 min, and a final step of 72°C for 10 min. Samples were analyzed in triplicate PCR reactions and pooled together. Pooled amplicons were visualized by agarose gel and DNA concentration quantified using a Qubit 2.0 dsDNA BR Assay Kit (Thermo Fisher Scientific, USA). No amplicons were visible on the agarose gel for the negative controls, which were then excluded from further analyses. Fifty nanograms of each pooled amplicon were pooled together, purified using Select-Asize DNA clean and concentrator (Zymo Research, catalog #4080), and quantified using Qubit. The concentration was adjusted to 2 nM and the library was denatured and diluted according to the manufacturer's instruction (Illumina, USA), with 10% PhiX added to the sequencing run and sequenced using the MiSeq Reagent Kit v3 (600-cycles) (Illumina), according to the manufacturer's instructions.

Bioinformatics and statistics
Raw data from sequencing was demultiplexed with deML 78 . Demultiplexed reads were imported and processed in QIIME 2, version 2018.8.0 79 ; forward and reverse reads were truncated using DADA2 80 at 250 bp, classified using the Greengenes database 81 , and phylogenetic tree created with FastTree 82 . The feature table and phylogenetic tree from QIIME 2 were imported to the Phyloseq package 83 in R 84 . Singletons and ASVs at a frequency <0.005% of the total number of reads 85 and ASVs with <15 reads in 4 samples were removed. Reads were normalized with the transform_sample_counts function to relative abundances. CCA was done using the vegan package 86 with the ASV table, UV as a numeric variable, and sample type (after UV exposure or stored for 6 days at 7°C) as a factor, Fig. 6 Heterotrophic plate counts of water exposed to different UV doses, and following storage. Water exposed to different UV doses (dark gray), and following storage for 6 days at 7°C (light gray). n = 6 for UV dose of 0 J/m 2 and n = 6 for storage not exposed to UV. n = 3 for all UV doses >0 and sample type. The boxes show the interquartile range, the line inside each box represents the median, whiskers show maximum and minimum values, and the crosses show outliers. the CCA output was followed by an ANOVA permutation test for CCA. Principal Coordinates Analysis plot using Bray-Curtis dissimilarity was done and visualized using Phyloseq and the ggplot2 package 87 , and ADONIS statistics performed using the vegan package. Data were agglomerated at different taxonomic ranks in Phyloseq and wrangled using Tidyverse R package 88 . Linear regression plots were plotted in R with ggplot2 87 , ggpubr 89 , and ggpmisc 90 packages, and the linear regression model calculated with the lm function, Pearson's correlation, and correlation coefficient calculated using the cor.test function in R. All other statistical analyses were done using base R. Non-normalized reads were used for differential abundance analysis using DESeq2 91 within Phyloseq (P adjusted < 0.05), excluding samples which were stored and accounting for the UV dose as a numeric parameter and visualized in ggplot2. Output ASVs from DESeq2 were visualized in a phylogenetic tree using ggtree 92 . Species diversity indices were calculated with non-normalized reads with the Phyloseq package (Shannon Index and observed ASVs) and custom scripts (Pielou's measure).
Log reductions of HPCs (H) were calculated using Eq. 1: where N 0 is the HPCs of the non-irradiated water and N the HPCs of the irradiated water at a specific dose. The SDs of log reductions (σ) were calculated with Eq. (2): using log 10 values of HPCs where σ N0 is the SD of N 0 and σ N the SD of N, n 0 is the sample size of the non-irradiated water, and n the sample size of the irradiated water at a specific dose.

DATA AVAILABILITY
DNA sequences are available at the NCBI Sequence Read Archive, accession number: PRJNA579449.