Enumerating viable phytoplankton using a culture-based Most Probable Number assay following ultraviolet-C treatment

Journal of Applied Phycology, Sep 2017

Ballast water management systems (BWMS) must be tested to assess their compliance with standards for the discharge of organisms, for example in the ≥ 10- and < 50-μm size category, which is dominated by phytoplankton. Assessment of BWMS performance with the vital stains fluorescein diacetate + 5-chlorofluorescein diacetate, required by regulations in the USA, is problematic in the case of ultraviolet-C (UVC) radiation. This is because UVC targets nucleotides—and thus reproduction, hence viability—rather than membrane integrity, which is assayed by the stains. The Serial Dilution Culture-Most Probable Number (SDC-MPN) method, long used to enumerate fragile phytoplankton from natural communities, is appropriate for counting viable phytoplankton. We developed QA/QC “best practice” criteria for its application as a robust and repeatable assay of viable cells in cultures of phytoplankton before and after experimental treatment, then constructed dose-response curves for UVC-induced loss of viable cells in 12 species of phytoplankton from seven divisions. Sensitivity to UVC, expressed as the dose required to reduce viability by 99%—the criterion for type approval of treatment systems—varied more than 10-fold and was not correlated with cell size. The form of the dose-response curves varied between taxa, with most having a threshold dose below which there was no reduction in viability. Analysis of the patterns of growth indicates that if recovery from treatment occurred, it was complete in 1 or 2 days in > 80% of cases, long before the assays were terminated. We conclude that the SDC-MPN assay as described is robust and adaptable for use on natural phytoplankton.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:


Enumerating viable phytoplankton using a culture-based Most Probable Number assay following ultraviolet-C treatment

Enumerating viable phytoplankton using a culture-based Most Probable Number assay following ultraviolet-C treatment Hugh L. MacIntyre 0 1 2 3 John J. Cullen 0 1 2 3 Trina J. Whitsitt 0 1 2 3 Brian Petri 0 1 2 3 0 Trojan Technologies , 3020 Gore Rd, London, ON N5V 4T7 , Canada 1 Department of Oceanography, Dalhousie University , PO Box 15000, Halifax, NS B3H 4R2 , Canada 2 Hugh L. MacIntyre 3 Abbreviations IMO International Maritime Organization SDC-MPN Serial-Dilution Culture-Most Probable Number USCG United States Coast Guard UVC Ultraviolet-C radiation Ballast water management systems (BWMS) must be tested to assess their compliance with standards for the discharge of organisms, for example in the ≥ 10- and < 50-μm size category, which is dominated by phytoplankton. Assessment of BWMS performance with the vital stains fluorescein diacetate + 5-chlorofluorescein diacetate, required by regulations in the USA, is problematic in the case of ultraviolet-C (UVC) radiation. This is because UVC targets nucleotides-and thus reproduction, hence viability-rather than membrane integrity, which is assayed by the stains. The Serial Dilution Culture-Most Probable Number (SDC-MPN) method, long used to enumerate fragile phytoplankton from natural communities, is appropriate for counting viable phytoplankton. We developed QA/QC Bbest practice^ criteria for its application as a robust and repeatable assay of viable cells in cultures of phytoplankton before and after experimental treatment, then constructed dose-response curves for UVCinduced loss of viable cells in 12 species of phytoplankton from seven divisions. Sensitivity to UVC, expressed as the dose required to reduce viability by 99%-the criterion for type approval of treatment systems-varied more than 10fold and was not correlated with cell size. The form of the dose-response curves varied between taxa, with most having a threshold dose below which there was no reduction in viability. Analysis of the patterns of growth indicates that if recovery from treatment occurred, it was complete in 1 or 2 days in > 80% of cases, long before the assays were terminated. We conclude that the SDC-MPN assay as described is robust and adaptable for use on natural phytoplankton. UVC treatment; Phytoplankton; Ballast water; Stress; Viability; Delayed repair Introduction Harmful algal blooms (HABs) can be costly in both economic and ecological terms. There is much debate over the conditions that give rise to them (e.g., Davidson et al. 2012; Glibert et al. 2012) , but one contributory factor is fundamental: viable cells of the bloom-forming organism must be present for a bloom to develop. There is a significant potential for transport and introduction of HAB taxa through ships’ ballast water (e.g., Doblin et al. 2004; Bolch and de Salas 2007; Roy et al. 2012) , which was one stimulus for the International Maritime Organization (IMO) adopting the International Convention for the Control and Management of Ships’ Ballast Water and Sediments at a Diplomatic Conference in London in 2004 (IMO 2004). The treaty was ratified by the requisite number of signatory nations in September 2016, meaning that ballast water treatment will be subject to international law one year later. National regulation is already in effect in the USA, governed by the Nonindigenous Aquatic Nuisance Prevention and Control Act (1990) and the National Invasive Species Act (1996). Regulatory authority resides with the United States Coast Guard (USCG), which entered into a formal Memorandum of Understanding with the US Environmental Protection Agency to develop procedures for evaluating ballast water management systems (BWMS). The regulations governing ballast water management address many categories of potentially invasive organisms. Two size classes of plankton are regulated: organisms that are ≥ 10 to < 50 μm in minimum dimension and those that are ≥ 50 μm in minimum dimension. Potential treatment methods include filtration, ultraviolet irradiation, heating, ultrasound, magnetic separation, electrolysis, chlorination, ozonation, and use of other chemica l biocides (reviewed by Ts olaki a nd Diamadopoulos 2010) . Germicidal UV (ultraviolet-C, UVC) radiation is a mature technology that is used extensively in disinfection of wastewater, drinking water, and recycled water and which was demonstrated to be an effective disinfectant for marine microbes more than a century ago (Henri et al. 1910) and to be effective with an alga more than 60 years ago (Redford and Myers 1951) . To be granted type approval, a BWMS must meet performance standards as demonstrated by a third-party testing facility to the satisfaction of the responsible administration. The standard for the 10–50-μm size range is fewer than 10 viable (IMO 2008) or living (USCG 2012) organisms mL−1, reduced from untreated concentrations of at least 1000 organisms mL−1, corresponding to a 99% reduction in viable or living organisms. The path to type approval therefore requires not only the development of effective treatment technologies but also reliable assays that can verify performance metrics. In this study, we use a 99% reduction as a metric to evaluate UVC treatment, recognizing that it is relevant for type approval of BWMS. It is important to recognize that this is different from doses that would be applied to meet regulatory compliance. Because a BWMS must treat to the regulatory standard of 10 viable/living cells mL−1, the applied doses might be higher or lower depending on the initial concentrations of cells in the 10–50-μm size range. The effectiveness of UVC is described using dose-response relationships (e.g., Hijnen et al. 2006; Haji Malayeri et al. 2016) . The primary targets of UVC radiation are nucleotides (Sinha and Hader 2002; Cadet et al. 2014) , and in the case of the phytoplankton, their dimerization can impair a cell’s ability to reproduce at doses that have comparatively minor effects on secondary targets such as photosynthesis (e.g., Redford and Myers 1951). Consequently, a definitive assay would be a test of the cells’ ability to reproduce (Breeuwer and Abee 2000) , measured in the relative reduction in numbers of viable cells as a function of UVC exposure. However, there are very few published dose-response relationships for phytoplankton that could be used to develop appropriate BWMS, although manufacturers undoubtedly have conducted their own tests. Calculations on data for the chlorophytes Tetraselmis suecica and Chlamydomonas reinhardtii (Chaudhari et al. 2015; Olsen et al. 2015) suggest that a dose of 300–490 mJ cm−2 would be required to effect the minimum 99% reduction of viable cell numbers required to meet the test standard (for phytoplankton in this context, Bcell^ is equivalent to a regulated Borganism^). The data are not truly comparable as the former used medium-pressure lamps (broad-spectrum UV) while the latter used low-pressure lamps (narrow-spectrum UVC). Even so, this is very high compared to the range of doses that give 99% disinfection in viruses, bacteria, and freshwater protozoa such as Cryptosporidium and Giardia and that have been used to design systems for treatment of drinking water (4–111 mJ cm−2 in data compiled by Hijnen et al. 2006; Haji Malayeri et al. 2016) . In contrast, field experiments with mixed assemblages suggest that phytoplankton might be less resistant than bacteria (First and Drake 2014) . There is, therefore, a need to quantify dose-response curves for viable phytoplankton cells to assess how effective UVC is as a treatment technology. Reliable enumerations of viable organisms are also needed for assessing the effectiveness of BWMS. Although testing based on SDC-MPN assays for viable cells (see Wright and Welschmeyer 2015) has been explicitly rejected by the US Coast Guard for type approval of UVC BWMS, the primary reason was because tests of viability do not use the criterion for counting living organisms, as required by US regulations (https://www.uscg.mil/foia/docs/FINAL-AG.pdf). Notably, BWMS type approvals have been granted by other administrations based on SDC-MPN assays (IMO 2017), and SDC-MPN is explicitly recognized as the appropriate assay for testing UVC-treated drinking water in the USA (USEPA 2006) . Discharge standards expressed in living or viable organisms mL−1 imply that assessment of BWMS performance requires classification and enumeration of individual organisms (for phytoplankton, cells) rather than a bulk measurement of vitality or viability (Zetsche and Meysman 2012) . Even so, there is significant value in Bindicative^ bulk parameters, such as variable fluorescence (Wright et al. 2015; Bradie et al. 2017; First et al. 2017) , for the rapid assessment of BWMS performance, with the recognition that such assays must ultimately be related to corresponding organism counts and status (i.e., viable or living). The need for enumeration of viable phytoplankton cells stands. While Bregrowth^ experiments—incubation of treated samples under conditions favorable for growth (e.g., Liebich et al. 2012; Stehouwer et al. 2015; Wright et al. 2015) —can demonstrate unequivocally the persistence of viable cells, they do not translate to the concentrations of viable cells remaining after treatment. However, quantitative estimates can be made by adaptation of a venerable technique used in microbiology, SDC-MPN (McCrady 1915; Cochran 1950) . The assay is based on testing for the presence of viable cells in serial dilutions of a sample by monitoring for growth. As dilution increases toward infinity and the probability of a single viable cell remaining in a diluted sample declines to 0, so does the probability of detecting growth in the sample. The most probable number of organisms in the undiluted sample is determined statistically from the number of tubes in a replicated dilution series in which there is evidence of growth (Cochran 1950; Throndsen 1978) . The SDC-MPN method has been used in several contexts in studies of microalgae and cyanobacteria. It has been used to enumerate cyanobacteria and diatom, raphidophyte, and dinoflagellate propagules from soils and sediments (Imai et al. 1986; Othman and Wollum 1987; An et al. 1992; Itakura et al. 1997; Harris et al. 1998; McQuoid 2002) and to estimate the abundances of difficult-to-preserve subpopulations within natural communities of phytoplankton (Knight-Jones 1951; Throndsen 1978; Furuya and Marumo 1983; Jochem 1990; Throndsen and Kristiansen 1991; Backe-Hansen and Throndsen 2002; Andersen and Throndsen 2003; Cerino and Zingone 2006; reviewed recently by Cullen and MacIntyre 2016) . More recently, it has been used to estimate viability after UVC treatment (Oemcke et al. 2004; First and Drake 2014; Casas-Monroy et al. 2016; Chaudhari et al. 2015; Wright and Welschmeyer 2015) . Although the SDC-MPN method is a benchmark test in microbiology, against which other methods of enumerating bacteria can be assessed (e.g., Russell 2000) , and although it has been used to enumerate phytoplankton for more than six decades, its applicability for enumerating total viable phytoplankton in a regulatory framework will require further validation, including Bbest practice^ procedures that have yet to be fully developed. In particular, as addressed by Cullen and MacIntyre (2016) , it is critical to assess whether the lack of observed growth in an assay tube is a confirmed absence of viable cells, or a false-negative error arising because the assay was not conducted for long enough. This primary uncertainty with the use of SDC-MPN on natural phytoplankton—the assumption that a species will grow under the assay conditions—is obviated when using lab cultures that are diluted and then incubated in growth conditions that are known to support vigorous growth. Experiments on cultures thereby facilitate optimization, replication, and interpretation of results. In this study, we have developed rigorous Bbest practice^ criteria for laboratory experiments using the SDC-MPN method and have quantified repeatable UVC dose-response relationships on cultures of phytoplankton from seven divisions. Although the qualitative and quantitative responses differed between taxa, replicated experiments assessed with SDCMPN showed that UVC was an effective disinfectant in all cases. We have also tested for the likelihood that the assay is biased by undetected delayed recovery and conclude that it is not an issue. Because the experiments were conducted on cells in a reproducible growth condition, the dose-response curves could be used to ground-truth other viability assays without the need to repeat the experiments, which involved monitoring > 3000 individual cultures over 2 ½ years. Materials and methods Culturing and chlorophyll a fluorescence Cultures of phy toplankton were obtained from the National Center for Marine Algae and Microbiota (NCMA, East Boothbay, ME, USA). The isolates were grown at 18 ± 0.3 °C on a 12:12 L:D photocycle at 130 ± 5 μmol photons m−2 s−1 of PAR, provided by Osram FL40SS-W/37 fluorescent bulbs. Cultures were maintained at low optical density in 40-mL volumes in balanced, semicontinuous, nutrient-replete growth (MacIntyre and Cullen 2005; Wood et al. 2005) by inoculation into fresh medium approximately once per generation (1–4 days, Table 1). The cultures were grown on f/2 (Guillard 1975) or L1 (Guillard and Hargraves 1993) seawater media. The growth media were prepared from tangential-flow-filtered coastal seawater (salinity c. 30) that was enriched with nutrients and delivered into autoclaved glassware through a 0.2-μm Whatman Polycap AS disposable filter capsule. Culture transfers were performed in a laminar flow hood that was wiped down with methanol prior to use. All tubes were flamed on opening and on closure to prevent contamination. All vessels used in culturing were cleaned prior to use by soaking in 5% Micro-90 cleaning solution (Z281506, Sigma Aldrich) and then in 2% HCl, rinsing copiously with E-Pure water (Barnstead Nanopure/APS Water Services Corporation, Lake Balboa, CA, USA) after each cleaning step. They were then sterilized in an autoclave. The cultures were grown in 50-mL borosilicate glass tubes with HDPE screw-caps that could be inserted while still sealed into the interrogation volumes of two fluorometers, so that fluorescence could be measured without compromising sterility (Brand et al. 1981) . Specific fluorescence-based growth rates (μF, day−1; Table 2) were calculated daily (Wood et al. 2005) from the change in dark-acclimated and dilutioncorrected Chla fluorescence (Eq. 1, Table 3). Fluorescence was measured daily with a Turner 10AU fluorometer (Turner Designs, USA) 3–5 h after the start of the night phase of the L:D cycle. Variable Chla fluorescence was also measured daily with a FIRe fluorometer (Satlantic, Canada). Both fluorometers were blanked with filtered seawater and the FIRe was standardized daily with a 100 μmol L−1 solution of rhodamine b (R6625, Sigma Aldrich) in E-Pure water. Estimates of the minimum, maximum, and variable fluorescence (F0, Fm, Fv; Table 2) were extracted from the single-turnover flash induction curves using Fireworx software (Audrey Barnett, http://sourceforge.net/projects/fireworx/). Curve fits were performed with the curvature parameter in the fit, p, set to 0. An essential element of the experimental approach was maintenance of the cultures in balanced growth using ACD (μm) tgen (days) semicontinuous culture (see MacIntyre and Cullen 2005) , thereby avoiding the variations of physiological status and life stage that are inherent during unbalanced growth (Harrison et al. 1990; Han et al. 2002; Guan and Lu 2010) , e.g., during typical batch-culture experiments. The requirement for balanced growth ensures that all cells are actively growing and fully acclimated to defined growth conditions. Under these conditions, the cells are demonstrably all viable (MacIntyre and Cullen 2016) , which is the appropriate starting point for tests of a stressor. Stock cultures were judged to be in balanced growth when the coefficients of variation (CV) for both daily growth rate and Fv/Fm were < 10% over 10 generations (Eq. 2, Table 3) for cultures for which μF was ≥ 0.5 day−1. Cultures in which μF was < 0.5 day−1 were judged to be in balanced growth when the coefficient of determination of a linear regression of ln-transformed, dilution-corrected fluorescence on time was > 0.995 and the CV for Fv/Fm was < 10% over 10 generations (Eq. 2, Table 3). The second set of criteria were adopted because stochastic error in estimating fluorescence effectively precluded meeting the criterion for < 10% variation in μF at low growth rates. When a culture had achieved balanced growth, an inoculum of 15–25 mL was transferred into a sterile 750-mL polystyrene tissue-culture flask. Its total volume was increased by adding fresh medium daily to maintain constant fluorescence at that time in the photoperiod until the volume required for experimentation was reached. The experimental cultures were monitored daily using the same protocols as for the stock lines to confirm that the test material remained in balanced growth. The only difference was that subsamples were decanted into sterile 50-mL culture tubes in the laminar flow hood and darkacclimated prior to measuring fluorescence. After measurement, the subsample was returned to the culture flask under sterile conditions in the hood. Limit of detection and average cell diameter The fluorometer’s lower limit of detection (LLDCell; Eq. 3, Table 3), the lowest cell concentration that could be distinguished from a blank, was measured for each species by regression (Anderson 1989; Miller and Miller 2005) . Dark-acclimated fluorescence was measured on aliquots of steady-state cultures that had been diluted to 20, 40, 60, 80, and 100% of the initial concentration with fresh medium. Cell concentration in each steady-state culture (N0, cells mL−1; Table 2) was determined with an Accuri C6 flow cytometer (Becton Dickinson, Canada). Counts were gated on Chla fluorescence (emission > 670 nm) at a level above the signals generated from 0.2-μmfiltered E-Pure water, to eliminate the contribution of noise and particles other than the phytoplankton. Fluorescence in the diluted samples was then regressed on cell concentration, obtained from N0 and % dilution, and LLDCell calculated from the standard error of the fluorescence estimates and the slope of the regression (Eq. 3, Table 3). The corresponding fluorescence intensity, LLDF (Arb.), was determined from LLDCell using the coefficients of the regression. Average cell diameter (ACD; Tables 1 and 2) was calculated from cross-sectional area, which was measured on a minimum of 50 cells per culture, using an AX10 microscope with AxioCam ERc5s camera (Zeiss, Germany). Cell areas were analyzed in ZEN software, calibrated using the lines on a Reichert Bright-Line hemocytometer (Hausser Scientific, USA). UV treatment Cultures were treated with defined doses of UVC radiation using a custom-built, low-pressure collimated beam (Bolton and Linden 2003) . The output of the collimated beam (mW cm−2) was measured before each use with a NIST-traceable ILT1700 radiometer (International Light Technologies, USA), which was calibrated annually. The dose (mJ cm−2) was applied by setting an appropriate time of exposure for 50 mL of culture stirred in a glass reaction vessel 5 cm in diameter and 2.5 cm in depth. The doses applied were all within the range of doses applied in municipal water treatment (USEPA 2006; NWRI 2012) . The appropriate time of exposure for each dose was calculated from the average intensity (mW cm−2) in the reaction vessel, calculated from the measured intensity at the sample surface and corrected for the reflection, Petri, water, and divergence factors (Bolton and Linden 2003). Transmittance was measured in a 1cm quartz cuvette with a P100 UV254 meter (Real Tech Inc., Canada). Transmittance at 254 nm in the cultures was 81 ± 7% (mean ± standard deviation; n = 136) for independent samples. The experimental regimen was designed to mimic that of a UV-based treatment system undergoing land-based testing for IMO approval: the IMO G8 protocol (IMO 2008) specified that such testing must include a minimum dark hold of 5 days between uptake and discharge. Accordingly, dose-response curves were based on division of a culture into 50-mL aliquots, irradiation of each, incubation in darkness for 5 days, and a second irradiation at the same dose. Where reported, the UV dose refers to the dose applied in each treatment, so is equivalent to half of the total dosage received by each sample. The culture tubes were maintained in darkness by shrouding them in four thicknesses of aluminum foil in the incubator. The following terms are used to distinguish between treatments: T0 samples are samples in balanced growth that have received no further manipulation; treated samples were either irradiated and maintained in darkness for 5 days as described above (UV-treated), or were held in darkness for 5 days but without UV irradiation (dark controls). SDC-MPN growth assays Cell viability was determined using the serial dilution culture method with MPN (McCrady Specification 1915; Cochran 1950) . Its application for phytoplankton (Throndsen 1978) was modified by using fluorescence rather than microscopy to monitor growth (cf. First and Drake 2014) . The SDC-MPN assay has a quantitative range of detection of 1.8–1600 cells when conducted with 5 replicate tubes at dilutions of 10−1, 10−2, and 10−3 (FDA 2010) . The combination of the dynamic range of the assay and the assay volumes gives the detectable range of viable cells (e.g., 0.36–320 cells mL−1 for 5-mL assay volumes). If the concentration of viable cells in the parent sample, Nviable, exceeds the detectable range, more dilution is required. So, to maximize compatibility between the dilutions and estimated Nviable, dilution ranges were set as 10X, 10X−1, and 10X−2, in which X is an integer based on the estimated number of viable cells and the median tabulated value of MPN (Eq. 4, Table 3). The estimate of viable cells came from either (i) total cell concentration, N0 (cells mL−1, determined by flow cytometry), in untreated T0 samples (demonstrably equivalent to Nviable, see BResults^ and MacIntyre and Cullen 2016) ; or (ii) from estimated survival of viable cells after treatment, obtained from N0 and preliminary range-finding tests on UVC-treated samples. For the range of Nviable assayed (3.6 × 10−1 to 1.8 × 105 cells mL−1), X ranged from 0 to − 6. Test samples were serially diluted with fresh, sterile growth media to 10X, 10X − 1, and 10X − 2, calculated according to Eq. 4 (Table 3). Five 5-mL replicates of each dilution were inoculated Eq. number 1 2a 2b 3 4 5 6a 6b 7a 7b into autoclaved 8-mL borosilicate tubes with HDPE caps. These grow-out tubes were incubated at a temperature and irradiance tested to provide rapid growth (see BResults^), and Chla fluorescence was monitored with the Turner 10AU fluorometer every 2–7 days, depending on growth rate. During the incubations, Chla fluorescence in the culture tubes could increase due to an increase of cell numbers and/ or due to recovery from a treatment-induced decrease in percell fluorescence. To account for treatment-induced changes in per-cell fluorescence that might be confounded with changes in cell number, per-cell Chla fluorescence was measured with the flow cytometer before and after 5 days in the dark with or without UVC treatment. It declined strongly with increasing UVC dose in Isochrysis, Phaeocystis, and Heterosigma and declined below the T0 value without a dependence on UVC dose in the other taxa. The 1st percentile in a frequency distribution of per-cell fluorescence, measured immediately after the 5-day treatment regime and normalized to the T0 value, was 0.125, an 8-fold reduction. Recovery of the T0 fluorescence intensity in the light to its initial value without cell division would therefore account for less than an 8× increase for 99% of the population. Consequently, we set the criterion for scoring positive growth as an increase of 8× in F (but see BResults^ and BDiscussion^). Tubes were scored as positive for growth if fluorescence increased 8× above either LLDF or the fluorescence measured immediately after inoculation, whichever was higher. Tubes were scored as negative for growth if the criterion for a positive score was not met during the period in which a detectable increase in fluorescence would be predicted from the growth rate of cultures that did grow, with an added margin of 5 generations (see below). When all grow-out tubes had been scored as either positive or negative, the most probable number of viable cells in the undiluted sample (± the 95% confidence intervals) was determined from lookup tables (FDA 2010) . The concentration of viable cells in the undiluted sample, Nviable, was calculated (Eq. 5, Table 3) from MPN by normalizing to the incubation volume and reconciling the dilution range used (10X, 10X − 1, and 10X − 2) to the range of 10−1, 10−2, and 10−3 used in the lookup tables. Viability was expressed relative to the untreated sample as relative viability (RV, dimensionless; Table 2), the ratio Nviable/N0. Results BBest practice^ procedures for the SDC-MPN growth as says We tested cell viability by monitoring growth in SDCMPN assays using Chla fluorescence. This required developing a set of best practice protocols that were constrained by rigorous and conservative quality control criteria designed to minimize experimental errors and maximize repeatability between experiments on the same species; these are summarized in Table 3. To optimize growth conditions in the assays while ensuring that they did not confound the outcome of tests, growth of three taxa (Thalassiosira, Amphidinium, and Heterosigma, Table 1) was tested in a matrix of irradiance x temperature for cultures that had been subject to a low dose of UVC (> 10−1 RV), a high dose (< 10−6 RV), or no dose (dark controls). Of the treatments tested, the highest apparent growth rates for each of the three taxa were consistently obtained at 22 °C and 260 μmol photons m−2 s−1 (Fig. 1). Critically, growth responses for each treatment were consistent over all temperatures and light intensities tested, i.e., controls and cultures subjected to the low dose of UVC always grew, and cultures subjected to the high dose of UVC always failed to grow and bleached. Consequently, the choice of growth conditions would affect only the time required to detect growth and score the test—because slower growth requires a longer assay—but would have no effect on the outcome of the test. For all other taxa tested, the growth rate at 22 °C and 260 μmol photons m−2 s−1 was tested prior to any SDC-MPN assays. The grow-outs were performed under these conditions except where μF was lower than at the conditions under which the cultures were acclimated (18 °C and 130 μmol photons m−2 s−1). This was the case for Lingulodinium, which was assayed under the acclimated conditions. Patterns of growth after inoculation demonstrated that it was necessary to monitor the grow-out tubes repeatedly to obtain accurate scores, rather than measuring fluorescence on inoculation and at an arbitrarily chosen end date. This is illustrated in Fig. 2, which shows the time-courses of fluorescence in replicate tubes of two fast- and one slow-growing cultures. Of the former, growth of Thalassiosira would be assessed accurately as 5 positive scores if a single test were done 7 or 14 days after inoculation. However, a decline in fluorescence after 7–11 days in Heterosigma would result in underestimates of the true score if single tests were done at 14 or 21 days (4 and 2 positives, respectively, rather than 5). In Lingulodinium, the relatively slow growth and the loss of fluorescence in two of the tubes during stationary phase would lead to positive scores of 0, 3, 5, and 3 if the test were only scored once at 7, 14, 21, or 28 days after inoculation, respectively. If a scoring threshold lower than the 8× increase in fluorescence were chosen, some of the tabulated positive scores would increase, but the general pattern and conclusions would be the same. It was equally important to define criteria for negative scoring to prevent a sample that had been growing slowly below the LLD from being misclassified because it was not monitored for long enough to show growth. We developed an objective method to do so, based on linear regression of the growth curve to a day-0 value. Other researchers (Redford and Myers 1951; Romero-Martínez et al. 2016) have used the intercept from the regression of a growth curve as a direct Fig. 1 Fluorescence-based specific growth rates, μF, of the diatom Thalassiosira weissflogii at different combinations of growth irradiance (Eµ) and temperature following treatment. Cultures were acclimated to 130 μmol photons m−2 s−1 at 18 °C and then subjected to the experimental treatment: 5 days in darkness at 18 °C ± UVC exposure before and after the dark holding period. Data are for a control cells, which were not irradiated, and b for cells exposed to a low or high dose of UVC. Note that the growth rates for the last are negative, indicating bleaching and/or loss of cells following the return to light. Data are shown as the mean ± standard deviation of 3 replicate cultures estimate of the number of viable cells—this treatment differs in using the regression as a quality-control criterion on the MPN estimate rather than as an estimate in its own right. An increase in fluorescence in the least-diluted samples can be fitted by linear regression of ln-transformed values greater than LLDF on time (Fig. 3a). The slope of the regression is the specific growth rate, μF, and the intercept, FInit, is an estimate of the fluorescence of viable cells on inoculation, assuming no lag in the growth response (see Cullen and MacIntyre 2016) . Because they are all treatments of the same sample, differing only in the degree of dilution, the increase in fluorescence in any tubes containing viable cells in the second and third tier of dilutions can be predicted on two assumptions: that the growth rate is the same and that the fluorescence of viable cells on inoculation is reduced by the ratio of the relative dilutions, 10−1 and 10−2. These are the Bpredicted^ lines in Fig. 3b, c. At sufficiently high dilution, the predicted line will correspond to the growth of a hypothetical fractional cell at inoculation, so the tubes that contain one or more cells will display growth curves above the line, consistent with the results in Fig. 3. Monitoring tubes that show no growth for the equivalent of a further 5 generations after the predicted time to reach LLDF (tend; Eq. 6a, Table 3) provides increased certainty in concluding that there is—and will be—no growth in a tube (Fig. 3b, c). Applying the findings of the analyses summarized here, the following procedure was established: any tube in which fluorescence has exceeded 8× LLDF or 8× the initial fluorescence by tend is scored as positive (Eq. 7a, Table 3), and any tube that has not shown an 8× increase in fluorescence by this point is scored as negative (Eq. 7b, Table 3). The duration of tend is determined by the 3-generation period used to define growth in the least-diluted tubes and by the 5-generation margin of observation beyond the predicted time of appearance above LLDF in the second and third tiers of dilutions (Fig. 4a; Eq. 6). The rationale for the latter was empirical: a frequency distribution of the observed minus expected time to LLDF in the second and third tiers of dilutions showed that 99% of tubes exhibiting growth would grow to LLDF within 5 generations (Fig. 4b; Y = 5 in Eq. 6, Table 3). If the 95th percentile of the distribution were used instead, the margin would be 3 generations (Fig. 4b; Y = 3 in Eq. 6, Table 3). Under these conditions, tend would be determined by the 3-generation duration of the window for any tube exhibiting growth before the predicted date (67% of the cases). For any tube exhibiting growth later than the prediction, tend would be determined by the time required to satisfy the criterion for confirming growth. Analysis of the data collected to define LLDF shows that the minimum increase in F required to demonstrate growth is much lower than the 8× increase in F corresponding to 3 generations. The minimal increase in fluorescence required to demonstrate growth is defined as the value at which the lower 95% confidence interval (Sokal and Rohlf 2011) is equal to LLDF (Fig. 4c). The calculated values for the different taxa (all species except Lingulodinium, Table 1) were 0.09–0.45 (median, 0.29) of LLDF. That is, an increase in fluorescence of less than 50% beyond that corresponding to LLDF could be reliably detected. The period of growth required to confirm a positive score could therefore be reduced significantly below 3 generations (8× increase in F) without compromising the accuracy of the test. Assessment of viability following UVC treatment by SDC MPN assay The protocols described above were used to estimate the relative viability of 12 species at different doses of Fig. 3 Time-courses of in vivo fluorescence in three serial dilutions of a culture of Eutreptiella gymnastica, with 5 replicates per dilution. Note log scales. Tubes in which fluorescence increases above 8× the lower limit of detection (LLDF) are scored as positive in the SDC-MPN assay. a The dashed line (Bfitted^) is an extrapolation of the exponential growth phase, determined by regression of ln-transformed fluorescence on time for points above LLDF (see text for details). The intercept is FInit and the slope is μF. b, c The dashed lines (Bpredicted^) are the expected timecourses of fluorescence for tubes that are 10 and 100 times more dilute that in a, respectively, assuming that a hypothetical fractional cell at the highest dilution can grow exponentially after inoculation; in practice, cultures with one or more viable cells at inoculation will show growth curves above LLDF and tubes with no cells will not. The shaded bars indicate a 5-generation period from the time at which fluorescence is predicted to exceed LLDF. Tubes without an increase in fluorescence at the end of this period (tend, right margin of the gray bar) are scored as negatives. The final MPN score in this example is 5–3–1 UVC. Because their sensitivity to UVC was unknown in advance, this necessitated a range-finding experiment for each to define the appropriate range of UVC doses and to determine the appropriate dilutions for the SDC-MPN assays (Eq. 4, Table 3). Based on the results of the range-finding experiment, each species was then tested over a range of UVC doses that lead to a maximum reduction in relative viability, RV, to 10−4–10−6. An SDC-MPN assay was also performed on the untreated parent culture (T0). As is expected if the SDC-MPN assay accurately enumerates actively growing cells, total measured cell concentrations in the T0 samples, measured with flow cytometry (N0), fell within the 95% confidence intervals of the estimates of viable cells in the T0 SDC-MPN assays in all cases. The limit of detectable relative viability, defined as the point at which the calculated concentration would be reduced to one viable cell in a grow-out tube, was 10–6.3 for Fig. 4 Determination of tend in SDC-MPN assays. a Observed time dependence of bulk fluorescence in the most-dilute tubes in an SDCMPN assay (10−6 dilution from Fig. 3). The expected growth (Bpredicted^; dashed line) was based on the FInit method (Fig. 3). Time is expressed in generation equivalents, tgen (Table 2). The duration of the 5-generation margin of error and the 3-generation (8×) increase in fluorescence is shown as shaded boxes. Note that observed F in the sole tube that exhibited growth exceeded LLDF 3 generations ahead of the expected time. b Frequency distributions of observed-expected differences in the time to reach LLDF in the second and third tiers of dilutions in all species (n = 910). The expected times were calculated from FInit (Fig. 3). Both observed and expected times to LLD were calculated by regression of ln(F) on time and expressed in terms of the generation time in the least-diluted tubes. The 95th and 99th percentiles of the distribution (dashed lines) were 3.14 and 5.06 generations, respectively. c Variation in F with cell concentration in dilutions of a culture of Eutreptiella gymnastica used to estimate LLDF (see text for details). Estimates of F from dilution by linear regression and the 95% confidence intervals of the estimates are shown as light and heavy lines, respectively. The minimal increase in fluorescence (ΔF) required to demonstrate growth is defined as the value at which the lower 95% confidence interval is equal to LLDF (inset). In this case, LLDF + ΔF was 1.35× LLDF (0.43 generations). The calculation requires that there be no variation in per-cell fluorescence Synechococcus, 10–5.8–10−5.9 for Isochrysis and Phaeocystis, and 10–4.0–10−5.1 for the other, larger eukaryotic taxa. Because dose-response curves for relative viability showed both Bshouldering^ and Btailing^ (Weavers and Wickramanayake 2001; Hijnen et al. 2006) , data were fitted to a biphasic model (Ren et al. 2016) modified to allow a shoulder (i.e., the lack of effect below a threshold dose). Data were fitted to the following equations simultaneously: RVD ¼ RVD¼0 for D ≤ DTh and for D > DTh RVD ¼ RVD¼0⋅½ð1−αÞ⋅expð−k1½D−DTh Þ þ α⋅expð−k2½D−DTh Þ where RVD is relative viability (Nviable/N0, dimensionless) for UVC dose D (mJ cm−2) applied before and after the 5-day dark hold; DTh is a threshold dose (mJ cm−2) below which viability is unaffected by dose; RVD = 0 is the relative viability after treatments at or below the threshold dose—this serves as the statistical estimate for the dark-control samples that are held in darkness but without UVC exposure, as distinct from the observed value. The response to UVC doses above the threshold is biphasic and dependent on dose in exceedance of DTh. The biphasic response is partitioned by the dimensionless coefficient α, which varies between 0 and 1. The sensitivity coefficients k1 and k2 (mJ−1 cm2) correspond to the two inferred targets with relative influences of 1 – α and α, respectively. In the case where k1 = k2 = k, the decline in viability follows first-order kinetics and Eq. 8b reduces to: RVD ¼ RVD¼0⋅expð−k½D−DTh Þ for D > DTh ð8cÞ Because of the nonlinear distribution of the data, these equations were fitted after log-transformation of RVD. The dose required to reduce relative viability to 10−2—the requirement for the IMO’s approval during land-based testing rather than the dose that would be required to meet the regulatory standard of 10 viable/living cells mL−1—is D[0.01] (mJ cm−2), which was obtained by solving Eq. 8 for the case where RVD = 0.01. Equation 8 describes a response that is a function of dose in exceedance of a threshold that is ≥0 mJ cm−2. Equations 8a, 8b (biphasic), or 8a, 8c (first-order) were fit simultaneously with a Levenberg-Marquardt method in Kaleidagraph 4.5 (Synergy Software) using the Boolean expression D ≥ DTh = 1 if true, otherwise 0, with DTh constrained to be ≥ 0. Fits for the biphasic model are reported where the RMSE for Eqs. 8a, 8b was lower than for Eqs. 8a, 8c and where the estimates of k1 and k2 were both positive. (In two cases, estimates of k2 were negative, which results in a local minimum in the response curve, with apparent increased viability at higher doses.) ð8aÞ ð8bÞ Coefficients of UVC sensitivity are presented as D′[0.01], D0Th, and k′, after scaling to the values found for Phaeocystis, the taxon most sensitive to UVC (Table 4). There were widespread differences between species in the relationship between relative viability and UVC dose (Table 4, Fig. 5). These could be grouped into three categories. In the first, exemplified by Synechococcus (Fig. 5a) and also observed in Thalassiosira and Isochrysis (not shown), the dose-response curve was biphasic with a threshold dose (a Bshoulder^ as described by Hijnen et al. 2006) . In the second form of response, exemplified by Scrippsiella (Fig. 5b) and also observed in Heterosigma (not shown), the response was biphasic but without a shoulder. The third form of response was a first-order (i.e., log-linear) response but only above a threshold dose, DTh. This is exemplified by Phaeocystis (Fig. 5c). The sensitivity of an organism to UVC treatment, e.g., as represented by D′[0.01], the dose required to reduce RV by a factor of 100, scaled to the minimum value found among the taxa (Table 4), varied by 32× (Table 4). This measure is determined by the threshold (D0Th ), the sensitivity to UVC doses above the threshold (k′), and effects of the 5-day dark treatment (RVD = 0). With the exception of Heterosigma (not shown) and Scrippsiella, all taxa had scaled threshold doses, D0Th, that were significantly higher than 0 (Fig. 5; Table 4), varying by 37× in magnitude (Fig. 5, Table 4). The lowest value was observed in Phaeocystis and the highest value was observed in Pyramimonas. The sensitivity coefficients, k′, which are also scaled by the dose required to effect a 100× reduction in viability in Phaeocystis, varied by more than an order of magnitude, too. The range was 34× (Table 4) for the first-order fits. These were not significantly correlated with the threshold dose (p = 0.10; n = 7 for species with first-order responses—data are not included for species with biphasic responses because of differences in the model structures). Changes associated with 5 days in the dark without UVC treatment were indicated by significant non-0 intercepts for 5 species (3 shown in Table 4): RVD = 0 was < 1 in Pyramimonas, Scrippsiella, and Synechococcus (Table 4) and also in Heterosigma and Isochrysis (not shown). None of these indices of UVC sensitivity was significantly (p < 0.05) correlated with cell size, as ACD, in the eukaryotic taxa—ACD could not be measured accurately in Synechococcus, which is, nonetheless, much smaller than the eukaryotes. The D′[0.01] in Synechococcus was in the mid-range of values among species. The wide differences in both the form of response and sensitivity to UVC were repeatable for independent tests of the same species. Replicate data for each taxon (different symbols in Fig. 5) were collected on independent cultures at intervals of weeks to months. As the cultures were all in balanced growth—i.e., fully acclimated to the growth conditions—there should have been no consistent differences between the replicate curves. This was confirmed by comparison Response Table 4 Parameters of the dose-response curve for UVC radiation and relative viability in 9 species of phytoplankton (Fig. 5). Absolute values of the fits are proprietary, so UVC doses and sensitivity coefficients are reported relative to the dose required to effect a 10−2 reduction in relative viability in Phaeocystis globosa, the most sensitive taxon studied. Relative viability was fitted to UVC dose using Eq. 8 (see text for details). The intercept corresponds to the effects of 5 days in darkness with no UVC treatment. D′[0.01] is the UVC dose required to reduce relative viability by a factor of 100. Taxa are sorted by type of response (see text for details). Dose-response curves for Isochrysis, Heterosigma, and Thalassiosira (not reported) were all biphasic. Coefficients are reported ± the standard error of the estimate, with the exception of the intercept and threshold dose, which are reported with the 95% confidence interval. If the 95% confidence interval of the intercept or the threshold dose encompassed 0, the curve was refitted without that term. These are indicated by dashes Intercept (RVD = 0, dimensionless) Threshold dose (D0Th, rel.) Sensitivity coefficient (k′, rel.) Partitioning coefficient (α, dimensionless) of the normalized residual errors in Eq. 8 at each dose between different test dates: in all cases, there was no significant difference (ANOVA, p > 0.05). We were unable to measure a dose-response curve for Lingulodinium (Table 1). The cultures were handled with extreme care while mixing and pipetting, and the number of viable cells in the T0 sample was in good agreement with the measured cell concentration (Nviable = 3400 cells mL−1, 95% CI 1160–8000 cells mL−1, vs N0 = 3543 cells mL−1). Even so, growth was detected in only 4 of the 12 samples subjected to the experimental treatments. Relative viability in these was low (2.3 × 10−4–2.0 × 10−3) and was not significantly correlated with UVC dose (p = 0.51). Although viable cells were detected in some of the UV-treated samples, no viable cells were detected in the dark controls. Delayed repair of UVC damage The descriptions of doseresponse above are predicated on the assumption that there was no delay in repair and recovery beyond the duration of the incubations—i.e., that the observed estimates of viability are an accurate reflection of the rates of disinfection. Delayed repair could be expressed as incomplete recovery of UVCimpaired growth rate during the incubation, thereby compromising the determination of tend (Eq. 6, Fig. 3), or a delay in the commencement of recovery—i.e., a lag, sensu Cullen and MacIntyre (2016) —long enough to prevent accurate detection of growth by the end of the assay. Incomplete recovery of growth rate is addressed through an analysis of the growth curves in the SDC-MPN dilutions. The specific growth rate (μF, Fig. 6a) reflects integrated cellular function and persistent debilitation (i.e., incomplete recovery) would result in significantly reduced growth rates as a function of UVC dose. In all species except Scrippsiella, there was no significant (p < 0.05) dose dependence of growth rates in the SDC-MPN assays (e.g., Fig. 6b). Even in Scrippsiella, there was only a significant decline in growth rate at the 2 highest doses, where relative viability was reduced below 0.075% (Fig. 5). The lack of dose dependency indicates that when growth became observable, i.e., when fluorescence rose above LLDF (see Fig. 3), the surviving cells had recovered fully from any damage imposed by the treatment. The stability of observed growth rates as a function of UVC dose indicates that if for any species there was a period of repair and recovery before growth resumed, i.e., a lag, it must have occurred before growth was detectable. The upper limit of the recovery period is therefore the time taken for fluorescence to reach LLDF in the least-diluted tubes of each SDCMPN assay for a given UVC dose. The duration of the lag can be estimated by comparison of FInit in treated samples with FInit in the corresponding untreated (T0) samples, with both normalized to the concentration of viable cells in the sample: F′ = F/Nviable and F0Init = FInit/Nviable (Fig. 7a). In untreated cultures, the cells are demonstrably all viable (MacIntyre and Cullen 2016) , so F0Init should accurately reflect the contribution of viable cells at time 0. If there is a dose-dependent delay in growth, effectively shifting the growth curve to the right, the value of F0Init in a treated culture will decline below F0Init in the untreated culture. This is illustrated in Fig. 7b. In 8 of the 12 species tested, there was a significant (p < 0.05) decline with UVC dose in log-transformed F0Init in treated cells (e.g., Fig. 7c), consistent with a progressive lag before recovery in growth rate. The exceptions were Alexandrium, Isochrysis, Synechococcus, and Thalassiosira. If repair and subsequent growth are assumed to be sequential, the duration of a lag period in which repair occurs can be calculated as the time that it takes for F0Init in treated cells to increase to the value in untreated (T0) cultures (Fig. 7b). Assuming exponential growth at the same rate observed when fluorescence is above LLDF: 0 C⋅μ F ½Treated −1 A where Lag is the duration of the lag phase (days); F0Init [Treated] and F0Init [Untreated] are the mean value of F0Init ð9Þ different scales on inset graphs. T0 values are for cultures in balanced growth. Replicates in the response curves, shown by different symbols, are from independent cultures and were assayed up to 3 months apart from each other. Relative viability is defined as the concentration of viable cells, Nviable, determined by SDC-MPN assay, normalized to the T0 cell concentration, N0, determined by flow cytometry. Error bars are the 95% confidence intervals of each MPN estimate. The lines are fits to Eq. 8. Parameter values are reported in Table 4 from the upper tier of dilutions of untreated samples and the value from a single tube of a treated sample, respectively; and μF[Treated] is the growth rate in the treated sample. Estimates of lag time are negative if the value of F0Init in a treated sample is higher than the mean value in the untreated samples. The relationship between lag phase and UVC dose for Pyramimonas (Figs. 6b and 7c) is shown in Fig. 8a. Values ranged from −0.1 to 2.0 days. For comparison, it took 4–12 days for the same tubes to reach LLDF. Only 2 species, Pyramimonas and Eutreptiella, had lag phases that were significantly (p < 0.05) correlated with UVC dose or with log-transformed relative viability. Even so, the variability—as the range of estimates from a second-order polynomial fit to the data—was low compared to the time required for the cultures to reach LLDF: 1.2 vs 7.2 days in Fig. 6 a Variation in fluorescence, F, in the upper tier of SDC-MPN assays of untreated (T0) and UVC-treated samples of Eutreptiella gymnastica (5 replicates for both assays). Only values above the limit of detection, LLDF (see Fig. 3), are shown. Regression of lntransformed data gives the specific growth rate (μF, day−1) and intercept (FInit, dimensionless). b Variation in specific growth rates with normalized UVC dose, D′ (cf. Fig. 5), in Pyramimonas parkeae. Data are from positively scoring tubes from the least-dilute tier of MPN assays. Data for the corresponding (untreated) T0 samples are shown for comparison. Open and closed symbols are for independent experiments Pyramimonas and 2.9 vs 10.6 days in Eutreptiella. There was no significant correlation between the duration of the lag phase and log-transformed relative viability for all species combined (Fig. 8b). Note that in 22% of the cases, the estimated lag time is negative, not unexpected for a calculation that is sensitive to uncertainty in the estimate of μF and where the estimates of μF in the treated and untreated cultures were very close (cf. Fig. 6a, b). The highest estimate of the lag time was 5.4 days but 58% of the estimates were < 1 day and 81% were < 2 days. For comparison, the time required to reach LLDF in the same samples was 3–26 days (median 9 and mean 10 days). When the lag time was expressed in generation-time equivalents by dividing by ln(2) ⋅ μ−1, the median value was 0.7 generations (Fig. 8c). The range was − 2.0 to 3.0 generations (mean 0.8): 59% of the values were < 1 generation, and 88% of the values were < 2 generations. Discussion Our application of the SDC-MPN method to assess loss of viability in a wide diversity of taxa (12 species from 7 divisions) confirms that UVC is an effective biocide for phytoplankton within the range of doses applied in municipal water treatment Fig. 7 a Variation in dilution-corrected fluorescence normalized to the estimated number of viable cells, F′, in the upper tier of MPN growth assays of untreated (T0) and UVC-treated samples of Eutreptiella gymnastica (same data as Fig. 5a; 5 replicates for each treatment). Only values above the limit of detection, LLDF (see Fig. 3), are shown. Regression of ln-transformed data gives the intercept ( F0Init, dimensionless). Close agreement in the two regressions is consistent with no lag in the resumption of growth of viable cells after treatment. The subscripts [Treated] and [Untreated] discriminate between estimates from these two treatments. b Variation in F′ with time illustrating the lag phase in recovery from a hypothetical experimental treatment. Regression of ln-transformed data gives the specific growth rate (μF, day−1) and intercept ( F0Init, dimensionless). The lag is the period (double-headed arrow) required for recovery of F′ in the treated sample to the T0 value for the untreated culture. This is calculated from the difference between the two values of F0Init and μF for the treated sample: in practice, it encompasses any trajectory through the stippled area between the two values. c Variation in F0Init with UVC dose in Pyramimonas parkeae. Data are from positively scoring tubes from the least-dilute tier of MPN assays. Data for the corresponding (untreated) T0 samples are shown for comparison. UVC dose is given in relative units (cf. Fig. 5). Open and closed symbols are for independent samples Fig. 8 a Variation in calculated lag periods in positively scored tubes from the least-dilute tier of SDC-MPN assays of Pyramimonas parkeae as a function of UVC exposure (same data as in Figs. 6b and 7c ). Data are from positively scoring tubes from the least-dilute tier of the assays. Data for the corresponding (untreated) T0 samples are shown for comparison. UVC dose is given in relative units (cf. Fig. 4). Open and closed symbols are for independent samples. Lag was correlated (R = 0.65, p < 0.05) with UVC dose. b Variation in calculated lag periods (mean of 1–5 tubes per (USEPA 2006; NWRI 2012) . Quantitative parameterization of the dose-response curves showed that there were wide differences in the pattern and magnitude of response to UVC (Fig. 5, Table 4). The finding depends critically on the development of defensible Bbest practice^ criteria for the SDC-MPN assay. Best practice criteria for SDC-MPN assays By definition, the SDC-MPN assay is a quantitative test of the ability of cells to proliferate by cell division. As such, it is the appropriate assay for a stressor such as UVC, for which the primary mode of action is damage to nucleotides (Sinha and Hader 2002) and the subsequent reduction in growth potential (Breeuwer and Abee 2000) . Because the assay is based on the ability to reproduce, rather than a proxy of a vital function, it is equally appropriate for assessing the effect of any treatment that is intended to render microorganisms harmless. The SDCMPN method has been used to estimate the loss of viability due to UVC treatment (Oemcke et al. 2004; First and Drake 2014; Chaudhari et al. 2015; Wright and Welschmeyer 2015; Casas-Monroy et al. 2016) , among other applications. However, to our knowledge, none of these studies has developed a formal best practice protocol that includes, critically, the definition of an appropriate duration for testing: reported tests were terminated at an arbitrary end point that ranged from 7 days (determination by variable fluorescence, First and Drake 2014) to B1 to 2 months in summer and 2 to 4 months in winter^ (determination by microscopy, Knight-Jones 1951) . Only by defining an appropriate end point for the assay can the absence of growth be ascribed with confidence to the absence of viable cells in a tube rather than an inability to detect slowly growing cells. We have described the formalism of defining the duration of the test in Table 3 (Eqs. 3–7) after verifying that the growth conditions were optimal for each culture that was tested. experimental exposure) with relative viability, RV, in the 12 species for which dose-response curves were measured (n = 133). There was no correlation (p = 0.52) between the variables. The dashed line is the median of the estimated lags and the dotted lines are the first and third quartiles of the data distribution. c Frequency distribution of estimated lag periods, expressed in generation-time equivalents (see text), for the same data The criteria for defining positive and negative scores (Eq. 7, Table 3) are objectively derived. While these are specific to detection by fluorescence, the same principles can be applied to other means of detecting growth. By definition, growth cannot be detected if it is invisible to the detector being used (i.e., is below the LLD), which sets the lower limit of a change in sensor output that can be used to demonstrate positive growth. The criterion that we adopted for detection by fluorescence—an increase of 8× above LLDF or the initial fluorescence, whichever was higher—is an aggressively stringent criterion. It was based on changes in per-cell fluorescence that might bias interpretation of fluorescence as cell number if they occurred when fluorescence readings were higher than LLDF. Per-cell fluorescence (F cell−1) is likely to vary systematically throughout growth in a batch culture with up- and downregulation of the absorption cross section and quantum yield of fluorescence in response to changes in light and nutrient availability (Falkowski et al. 1981; Kolber et al. 1988; Geider et al. 1993) . In practice, though, changes during the period of exponential growth are likely to be minimal compared to recovery from stresses imposed by the prolonged darkness and/or UVC. The time required to affirm positive scores could be reduced, with the attendant savings in personnel time and effort, by accepting a less demanding threshold, and there are good grounds for doing so. As the cultures all had consistent growth rates, based on measurements of F above the LLD—with the exception of samples of Scrippsiella subjected to the two highest UVC doses—it is logical to assume that any repair and recovery from the imposed stresses was complete before growth was detectable. The estimated period for repair and recovery was 1–2 days in the vast majority (81%) of samples, consistent with the rapid rate of repair of cyclobutane pyrimidine dimers induced by UVC treatment in the chlorophyte T. suecica (Hull et al. 2017) , and significantly shorter than the average of 10 days required to reach the LLD. It is reasonable to assume that under these conditions, F cell−1 is constant and so increases in cell concentration could be assessed based on the observations of F. A pragmatic alternative to the 8×-LLDF criterion would be to accept an increase in F of 50% (based on the 95% confidence intervals of mean F cell−1 in independent cultures of 0.07–0.20× the mean) above LLDF as evidence of a positive score. Adoption of the more relaxed criterion would reduce monitoring time of tubes that score positive by several days and would eliminate redundant monitoring after growth was reliably detected. Optimizing the precision of SDC-MPN assays The under lying principle of the SDC-MPN assay is that the concentration of viable cells is diluted to the point at which a statistically interpretable number of tubes should contain no viable cells. This condition is satisfied when the dilution series encompasses a concentration of 1 viable cell per tube (Cochran 1950) , a condition that can be satisfied if the likely range of Nviable can be constrained with confidence. Because the concentration of viable cells is generally not known, the range of the dilution series—(DR)q − 1, where DR is the dimensionless dilution ratio (fixed at 10 in this study; Eq. 4, Table 3) and q is the number of tiers—should be large enough to ensure that no experiment yields all sterile or all fertile tubes. To accommodate uncertainty in the estimate of Nviable, the range can be expanded by increasing the number of tiers, at the cost of requiring additional resources. The dilution ratio can be increased, but at the expense of the method’s precision, which can be improved with compensatory increases in the number of replicate tubes per tier (Cochran 1950) . Given that increasing the number of dilutions, tiers or tubes per tier can require a significant investment of time and resources, it is useful to examine the quantitative relationships between the precision of SDC-MPN and the design of the test. Hurley and Roscoe (1983) presented a general equation to explore the influences of experimental design on the precision of the estimate of Nviable: σlnðNviableÞ ¼ ½V ⋅CFðiÞ 2⋅r ðN viableÞ2⋅∑iq¼1 exp½V ⋅CFðiÞ⋅N viable −1 #−1=2 ð10Þ where σlnðNViableÞ is the standard error of ln(Nviable), V (mL) is the volume incubated (Table 2), q is the number of tiers with r replicates per tier, and CF(i) (dimensionless) is the concentration factor in tier i. Each value of CF(i) is the concentration factor in tier i, determined from CF in the least-diluted tier and the dilution ratio between tiers, DR (dimensionless). In this study, the former is 10X and the latter is 10 (Eq. 4, Table 3). Different dilution ratios can be accommodated in Eq. 4 by using different bases for the power function that defines DR. A sensitivity analysis was used to examine the effect of varying q and r on the expected precision of the estimates of Nviable. It included the less important influences of V and Nviable. The analysis used 70,400 combinations of input variables, including those used by testing organizations charged with type approval of BWMS (e.g., IMO 2016, see Table 5) . The dilution ratio and number of tubes explained 96% of the variability in a log-transformed regression (70,397 df). Expressing the independent variables in linear form: σlnðNviableÞ ¼ 0:65⋅ DR0:266 pffirffi ð11Þ The result is generally consistent with the simpler relationship presented by Cochran (1950) , with a mean ratio (Eq. 11/ Cochran) of 0.93, s.d. = 0.03. The lower and upper 95% confidence limits on an MPN estimate of Nviable are N viable=exp 2σlnðNviableÞ and N viable⋅exp 2σlnðNviableÞ , respectively. The factor exp 2σlnðNviableÞ serves as an appropriate substitute for a standard error (cf. Cochran 1950) . For the design used in this study, DR = 10 and r = 5; the factor for estimated 95% confidence intervals is 2.9, generally consistent with the tabulated SDC-MPN confidence limits for individual experiments (Fig. 5). Precision is maximized by increasing the number of replicate tubes (r) and minimizing the dilution ratio (DR), as illustrated in Fig. 9. Doubling the number of tubes from 5 (this study) to 10 and reducing the dilution ratio from 10 (this study) to 2 would reduce the mean uncertainty in the estimate of viable cells (i.e., σlnðNViableÞ ) by 56% (Fig. 9). This would reduce the ratio of the 95% confidence interval from 0.26– 3.86× the estimate to 0.63–1.59× the estimate (medians in Fig. 9). A reduction in DR would necessitate increasing the number of tiers in the assay (q) if Nviable were not well constrained. Consequently, the improvement in precision would come at the expense of personnel time in monitoring the increased number of tubes, q × r. Minimizing the duration of SDC-MPN assays For logistical reasons, it is advantageous to minimize tend. With the current criteria, tend is determined by the 5-generation safety margin after the most dilute sample is predicted to reach the LLD rather than the 3-generation (8× increase) period taken to verify growth when fluorescence is higher than the LLD (Eqs. 6 and 7, Table 3). Only reducing the duration of the safety margin would reduce the duration of the test, tend. This is predicated on three factors, the specific growth rate, the number of viable cells in the sample, and the sensitivity of the sensor used to monitor growth (Eq. 6, Table 3, see also Cullen and MacIntyre 2016) . The duration of tests in this study was Concentration (Nviable, cells mL−1) Culture volume (V, mL) Number of tiers (q) Replicate tubes per tier (r) Dilution ratio (DR) Initial dilution (CF(1)) a Number of tiers is restricted to a total dilution range of 100 ≤ DRq − 1 ≤ 10,000 constrained by the availability of a sensitive instrument and ranged from 15 to 90 days where Nviable was 10 cells mL−1, and considerably longer (19–115 days) where the predicted concentration of viable cells was reduced to the limit of the test for 5-mL volumes, 0.36 cells mL−1. Reducing the length of the safety margin from 5 to 3 generations (corresponding to 99 to 95% confidence in detecting late-growing tubes in the most-diluted tier of the assay; Fig. 4b), enhancing growth rate (Fig. 1), and using a more sensitive detector (i.e., with a lower Fig. 9 Relative uncertainty in estimates of Nviable from SDC-MPN, σlnðNviableÞ, as a function of the dilution ratio (DR) and number of tubes per dilution tier (r), calculated according to Eq. 10 (Hurley and Roscoe 1983). The second y-axis shows the ratio of the upper 95% confidence limit to the estimate of Nviable, exp 2σlnðNviableÞ : Most of the variability in 70,400 solutions was explained by DR and r, shown here. The median value for σlnðNviableÞ is plotted for each combination; vertical bars show 5th and 95th percentiles for DR = 10 (1200 solutions) and 2 (400 solutions). The ranges for DR = 4 and 6 are not shown; they are intermediate. A general solution is presented in Eq. 11 Table 5 Input variables for an analysis of the influence of experimental configuration on the precision of SDC-MPN (Fig. 9). A range of nominal concentrations of viable cells, Nviable, was used in Eq. 10, along with all combinations of the input variables listed below. Examples of configurations used by three testing organizations (DHI, Denmark and Singapore; NIVA, Norway; and Golden Bear, USA) come from submissions to IMO (IMO 2016) . The concentration factor for the least-dilute first tier, CF(1), was calculated from Nviable and the dilution ratio DR: the original sample is diluted by DR successively until the number of detectable organisms per tube, V ⋅ CF(1) ⋅ Nviable, is ≤DR. Consequently, the number of viable organisms per tube in tier 1 is ≥ 1.0 and in tier 2 it is ≤ 1.0 Input variable Input values (total number) Nominal configurations DHI 6 3 3 replicates 5 each 10 NIVA Golden Bear 6 3–5 3–5 10 3 7 5 6 LLDCell) would all shorten the duration of the test, with the attendant savings in effort and time. We consider these in turn. We calculate values of tend for the third tier of an MPN assay for all eukaryotic species on which dose-response curves were measured (i.e., all except Lingulodinium, Table 1) using Eq. 6b. The starting concentration of viable cells, N0, was set at 10 cells mL−1; μF was the acclimated growth rate; and LLDCell was measured (see Fig. 4c). The calculation uses a 5-generation safety margin (Table 3). The calculated values of tend range from 13 to 40 days (median 21 days), as shown in Fig. 10a. Reduction of the safety margin to 3 generations reduces the corresponding values of tend to 11–34 (median 18) days (Fig. 10b). The reduction is disproportionately weighted to slow-growing cells because it scales with μF−1 (Table 3, Eq. 6). Because tend is inversely proportional to growth rate, the duration can be shortened by increasing the growth rate. Doing so by optimizing the growth conditions of temperature and irradiance (cf. Fig. 1) reduces tend by a factor proportional to [change in μF]−1 (Eq. 6, Table 3). This is illustrated in Fig. 10b for Thalassiosira and Alexandrium, the fastest- and slowest-growing cultures used in this analysis, respectively. Optimization of growth rate reduced tend from 11 to 7 days in the former and from 34 to 23 days in the latter, both calculated with a 3-generation safety margin. Any attempt to increase growth rate by increasing temperature and/or irradiance should be tested before conducting the assays, though (e.g., Fig. 1). Temperature optima vary between taxa, with significant genotypic variation by latitude (Thomas et al. 2012) , but growth rates frequently decline very rapidly at superoptimal temperatures: the median difference between the temperature optimum and the temperature at which growth was completely Fig. 10 Variation in tend, the time required to complete a 3-tiered assay, as a function of growth rate (μF) and the starting concentration of viable cells (Eq. 6, Table 3). The latter is expressed relative to the lower limit of detection for the strain (Nviable/LLDCell). Note the log scale for time. The duration of the assay is shortest under conditions of high initial concentration of viable cells (Nviable) and/or high sensitivity of detection (low LLDCell) and high growth rate. The points overlaid on the plot are the values for the eukaryotic taxa in Table 1, except Lingulodinium, with a starting concentration of 10 cells mL−1 and the measured LLDCell. An increase in growth rate and/or a higher sensitivity results in shorter assay duration. This is illustrated for the fastest- and slowest-growing taxa, Thalassiosira weissflogii and Alexandrium andersoni. a Values of tend abolished in 163 cultures of phytoplankton was 7.1 °C (Cullen and MacIntyre 2016) . In a similar vein, changes in growth irradiance should be tested prior to being imposed during grow-out to avoid exceeding genotypic and phenotypic limits on the irradiance at which photosynthesis and growth become light-saturated (MacIntyre et al. 2002). Imposition of significantly higher irradiances will result in photoinhibition of growth rather than more rapid growth (e.g., Kana and Glibert 1987) . Finally, the use of a more sensitive fluorometer affords a further reduction in tend by a factor proportional to [ln(change in LLDF−1)] (Eq. 6, Table 3). This is illustrated in Fig. 10c, again for Thalassiosira and Alexandrium. The use of a customized in situ FIRe fluorometer (Satlantic, Halifax, NS, Canada) reduced LLDCell by approximately an order of magnitude, from 117 to 13 cells mL−1 for Thalassiosira and from 91 to 14 cells mL−1 for Alexandrium. This would result in further reductions in tend of 1 and 4 days, respectively (Fig. 10c). As these factors are additive, the cumulative reductions for the two species were from 13 to 6 days for Thalassiosira and from 40 to 19 days for Alexandrium. The bulk of the reductions, 57 and 52%, respectively, can be attributed to optimization of growth rates, with an additional 29% reduction by using a 3-generation safety margin and 14–19% from using a more sensitive fluorometer. The duration of the assay could also be reduced by increasing the number of viable cells per tube through preconcentration of the parent sample. While this does have the advantage of reducing the inherent sampling bias in the case of very low cell concentrations (Frazier et al. 2013) , there is a risk of damaging cells are based on the growth rate in balanced growth (Table 1), LLDCell determined for a Turner 10AU fluorometer, and a 5-generation margin for detection (Eq. 6b). b Estimates as in a but for a 3-generation margin for detection. The white arrows and points show the reductions in tend achieved for Thalassiosira (1) and Alexandrium (2) by optimizing the irradiance and temperature conditions in the assay to maximize growth rates. c Estimates as in a but using optimized growth rates (Fig. 1). The white arrows and points show the reductions in tend for Thalassiosira (1) and Alexandrium (2) achieved by using a customized Satlantic FIRe fluorometer with limits of detection c. one order of magnitude below those of the 10AU during concentration, resulting in an underestimate of the number of viable cells. Again, pretesting is warranted before implementing the approach. Dose-response of UVC treatment Adoption of the criterion of balanced growth prior to experimentation allowed definitive attribution of losses in viability to the experimental treatments (e.g., darkness ± UVC) rather than the experimental treatment plus confounding factors associated with growth conditions. The test of this assertion is repeatability, defined here as the agreement between independent assays conducted days or months apart. For our experiments, repeatability is verified by the observation that in all cases (repeated tests on each of 12 taxa), the results for replicate assays conducted days to months apart (different symbols in Fig. 5) were not statistically distinguishable from each other. The repeatability of responses from cells in balanced growth provides a significant advantage beyond confidence in the results: the UVC dose-response curves could be used as benchmarks for future tests of viability assays without repeating the SDC-MPN analyses, if the same cultures were grown under the same conditions. This would be of considerable value when developing new assays for viability, as the SDC-MPN assays described here required repeated monitoring of more than 3000 individual cultures over a period of 2 ½ years. Accepting that the dose-response curves presented here are accurate and that the criteria for test duration are appropriate requires an assumption that results are not significantly influenced by delayed repair of damage caused by UVC. Any recovery that occurred after the grow-outs were completed and had been scored would cause the SDC-MPN assay to underestimate the viability of the test populations. The first line of evidence that this does not occur is the observation that there was almost no effect of treatment on the growth rates of the treated samples by the time that growth could be observed in the grow-outs (i.e., when fluorescence had risen above LLDF). Because the primary target of UVC damage is the nucleotides (Sinha and Hader 2002; Cadet et al. 2014) , persistent, nonlethal damage would be manifested as impaired protein synthesis and a subsequent reduction in growth rate. By inference, recovery of growth rate after treatment demonstrates that any nonlethal damage arising due to treatment must have been repaired before growth was observed. Analysis of the ratio of the fluorescence attributable to viable cells relative to viable cell concentration, F0Init, indicates that the repair period was 1–2 days in the overwhelming majority of cases (Fig. 8b). This corresponds to a median lag of less than 1 generation, which is consistent with the rapid rate of repair of UVC damage observed in the chlorophyte T. suecica (Hull et al. 2017) . It was significantly shorter than the median time for the growing populations to reach the LLD, so it would have been completed before growth was observable. Consecutive repair and then growth is consistent with prior observations of recovery from UV radiation, in which the time-scale of repair of cyclobutane pyrimidine dimers is hours (van de Poll et al. 2002; Rastogi et al. 2011) , resulting in resumption of growth in 2 days or less (Takahashi et al. 2006). There are very few published UVC dose-response curves for phytoplankton, and responses published for other taxa (Hijnen et al. 2006; Haji Malayeri et al. 2016) are not directly comparable with our data for two reasons. First, the dose dependencies presented here are proprietary so cannot be compared in absolute units. Second, most published curves are based on single-dose UVC exposure, in contrast to the regimen used here of UVC exposure/dark hold/UVC exposure. This was intended to simulate treatment of ballast water at intake and discharge of a 5-day voyage and by its nature is a dualstressor treatment for a photoautotrophic organism. If representative, the dose-response curves for Microcystis aeruginosa (Sakai et al. 2011), C. reinhardtii (Chaudhari et al. 2015) , and T. suecica (Olsen et al. 2015) indicate that phytoplankton are more resistant to UVC than the viruses, bacteria, and heterotrophic protists reviewed by Hijnen et al. (2006) and Haji Malayeri et al. (2016). However, the doses required to treat the species used in this study are within the ranges used to treat municipal wastewater (USEPA 2006; NWRI 2012) . A striking feature of the majority of dose-response curves presented here is the threshold below which UVC had no discernible effect on viability, often referred to as Bshouldering^ in the literature on disinfection (Weavers and Wickramanayake 2001; Hijnen et al. 2006) . A minority of samples also had biphasic responses, often referred to as Btailing.^ Our empirical model (Eq. 8) reproduces tailing as the combined effects of two targets with differing influences, represented by a partitioning coefficient, α, and respective sensitivities to UVC, k1 and k2 (mJ−1 cm). We used Eq. 8 to describe patterns in dose-response curves and to group them by species (Fig. 5, Table 4) without considering their mechanistic foundations. A general model of disinfection was proposed by Jensen (2010) , in which relative viability is described from the cumulative distribution function of a normal distribution. The model captures much of the variability observed in real-life responses, including both shouldering and tailing, but does so from a single damage point, rather than the two that are specified in Eq. 8. Although elegantly simple, the structure is almost certainly not representative of metabolic control, which is likely to be distributed through multiple components of a metabolic pathway (Kacser and Burns 1973) . Application of Jensen’s model to our results is complicated by the fact that our study was based on application of two stresses, prolonged darkness without and with UVC. This makes it more unlikely that damage would correspond to a single structural target. Even so, the model presents a useful conceptual framework for considering the mechanisms that might underlie the observed responses. A period of darkness is a predictable part of the day/night cycle and may be important for maintaining optimal growth rates (Brand and Guillard 1981) . Nevertheless, extended darkness results in debilitation of cellular pathways as the demand for maintenance respiration is not balanced by autotrophic carbon fixation. While some mixotrophic phytoplankton are capable of heterotrophy in darkness (Martínez and Orús 1991; Manoharan et al. 1999; Lee 2001) , in many others, mixotrophic carbon acquisition is light-dependent (Li et al. 1999; Walsby and Juttner 2006; Brutemark and Granéli 2011) . Respiratory demands are therefore met by catabolic consumption of energy-storage compounds and structural components of the cells, resulting in degradation of photosynthetic competence (Yentsch and Reichert 1963; Popels et al. 2007; McKie-Krisberg et al. 2015) . Ultimately, prolonged darkness can trigger programmed cell death (Franklin and Berges 2004; Segovia and Berges 2005) . Statistically significant reductions in viability in the dark controls were observed in only 3 of the 9 species described here, accounting for mean reductions of 55–72%. This suggests that the stress imposed by darkness alone had not triggered mass mortality in the cultures. Light is critical to recovery from UVC irradiation. The primary damage is from DNA lesions, notably cyclobutane pyrimidine dimers and pyrimidine 6-4 pyrimidone photoproducts, CPDs and 6-4PPs (reviewed by Sinha and Hader 2002; Rastogi et al. 2010a) . Nucleotide repair is effected by multiple pathways and illumination is a direct component of one. Photorepair is catalyzed by light-sensitive flavoproteins, the cryptochrome/photolyase family, which are involved in lightsensing and signaling processes that both regulate protein synthesis and directly bind and repair DNA (reviewed by Oliveri et al. 2014; Lepetit and Dietzel 2015) . There is also excision repair of damaged nucleotides (e.g., Rastogi et al. 2010a; Chaudhari et al. 2015) . In contrast to photorepair, this is not light-dependent, so is often referred to as dark repair. There have been very few studies of photolyases in the phytoplankton; even so, it is clear that they differ in structure and location across taxa (Brunelle et al. 2007; Coesel et al. 2009; Heijde et al. 2010) . Variability in rates of repair and the relative magnitudes of photorepair vs dark repair may underlie the wide ranges in DTh and k observed here. Irradiation with UVC results secondarily in the formation of reactive oxygen species that cause further damage, critically to the psbA (D1) protein in the PSII reaction center (Rastogi et al. 2010b; Wang et al. 2015) . Photosynthesis provides the energy required for catabolic metabolism—including both repair and cell division—beyond what can be supported by internal energy-storage reserves. Sakai et al. (2011) measured viability following UVC irradiation in cultures of Microcystis with and without photorepair. A pronounced shoulder in the dose-response curve for cultures permitted photorepair and its absence in cultures without it, and the convergence of the curves at high doses suggests that loss of viability at low dose was primarily from nucleotide damage but the loss at high doses was primarily from another source. A dose-dependent transition between targets is implicit in Eq. 8b and may underlie the biphasic response observed in 5 of the 12 taxa studied. Although ascribing mechanism to the wide differences in dose-response curves observed here can only be conjectural, it is clear that there are differences. The Bbest practice^ approach used in this study, when coupled with maintenance of the cultures in balanced growth prior to experimentation, allowed demonstration that the measurements were consistent over time and between independent cultures. Beyond their utility in defining design specifications for shipboard treatment systems, the response curves are true indicators of viability—i.e., growth and reproduction (Breeuwer and Abee 2000) —following treatment. They could and should be used as benchmarks for testing other assays of viability, many of which are (strictly) assays of vitality (e.g., intact membranes, nucleotide functionality, and/or metabolic competence, Breeuwer and Abee 2000) . The tests can be conducted by reproducing the growth conditions and UVC exposures for the tested strains— relative viability can be specified for each exposure without the need to repeat the SDC-MPN assay. Given the recent ratification of the IMO’s Ballast Water Convention, and the need to assess treatment efficacy of ballast water treatment systems for certification and for shipboard compliance testing, there is a pressing need for rapid, reliable, and accurate assays. Candidates that have been proposed include the vital stains fluorescein diacetate and chloromethylfluorescein diacetate (FDA and CMFDA) required by US regulations (USEPA 2010) and variable fluorescence (e.g., First and Drake 2013, 2014) . However, these assay membrane integrity and metabolic competence, respectively, rather than nucleotide functionality, which is the target for UVC damage. We have tested both against the SDC-MPN method for the 12 taxa described here and found that neither provides an accurate estimate of UVC-induced loss of viability (MacIntyre et al., unpublished). Conclusion We argue that the appropriate metric for measuring invasive potential of organisms in ballast water is viability, the ability to reproduce, and not vitality, the manifestation of functional components of the cellular apparatus. Adaption of the Most Probable Number assay—a benchmark test used in microbiology—and development of appropriate Bbest practice^ criteria for its application with phytoplankton allow the loss in viability following application of a biocide to be assessed repeatedly and consistently in laboratory experiments. When used to assay the responses of cells maintained in balanced growth under well-defined conditions, the SDC-MPN assay provides repeatable estimates of the dose-response function for UVC-based loss of viability. This is suitable for comparison with proxies of treatment effectiveness, even if experiments are conducted at different times. Acknowledgements We thank Shannah Rastin, Jessica Miller, and the members of Dalhousie’s Science Coop program who worked hard to collect the data: Gordon Byron, Sarah De Mendonça, Katerina Fupšová, Adam Harman, Anna Haverstock, Jessica Hurtubise, Isabelle Jubinville, Benjamin King, John MacDonald, Emma Marotte, Rhea Newman, Elayne Tan, and Jude van der Meer. We also thank an anonymous reviewer for a careful and constructive evaluation that significantly improved the manuscript. Funding information This work was supported by Collaborative Research and Development Grant—Project 445451-12 from the Natural Sciences and Engineering Research Council of Canada and Trojan Technologies, with supplementary funding from the Strategic Cooperative Education Incentive Program of Nova Scotia and Trojan Technologies. Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / / creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. An KH , Lassus P , Maggi P , Bardouil M , Truquet P ( 1992 ) Dinoflagellate cyst changes and winter environmental conditions in Vilaine Bay, southern Brittany (France) . Bot Mar 35 : 61 - 67 Andersen P , Throndsen J ( 2003 ) Estimating cell numbers . In: Hallegraeff GM , Anderson DM , Cembella AD (eds) Manual on harmful marine microalgae . UNESCO , Paris, pp 99 - 129 Anderson DJ ( 1989 ) Determination of the lower limit of detection . Clin Chem 35 : 2152 - 2153 Backe-Hansen P , Throndsen J ( 2002 ) Occurrence of pico- and smaller nanoplanktonic flagellates in the inner Oslofjord, eastern Norway, during the breeding season of the blue mussel (Mytilus edulis L .) Sarsia 87 : 65 - 74 Bolch CJS , de Salas MF ( 2007 ) A review of the molecular evidence for ballast water introduction of the toxic dinoflagellates Gymnodinium catenatum and the Alexandrium Btamarensis complex^ to Australasia . Harmful Algae 6 : 465 - 485 Bolton JR , Linden KG ( 2003 ) Standardization of methods for fluence (UV dose) determination in bench-scale UV experiments . J Environ Eng ASCE 129 : 209 - 215 Bradie J , Broeg K , Gianoli C , He J , Heitmüller S , Lo Curto A , Nakata A , Rolke M , Schillak L , Stehouwer P , Vanden Byllaardt J , Veldhuis M , Welschmeyer N , Younan L , Zaake A , Bailey SA ( 2017 ) A shipboard comparison of analytic methods for ballastwater compliance monitoring . J Sea Res . https://doi.org/10.1016/j.seares2017. 01 .006 Brand LE , Guillard RRL ( 1981 ) The effects of continuous light and light intensity on the reproduction rates of twenty-two species of marine phytoplankton . J Exp Mar Biol Ecol 50 : 119 - 132 Brand LE , Guillard RRL , Murphy LS ( 1981 ) A method for the rapid and precise determination of acclimated phytoplankton reproduction rates . J Plankton Res 3 : 193 - 201 Breeuwer P , Abee T ( 2000 ) Assessment of viability of microorganisms employing fluorescence techniques . Int J Food Microbiol 55 : 193 - 200 Brunelle SA , Hazard ES , Sotka EE , Van Dol ah FM ( 2007 ) Characterization of a dinoflagellate cryptochrome blue-light receptor with a possible role in circadian control of the cell cycle . J Phycol 43 : 509 - 518 Brutemark A , Granéli E ( 2011 ) Role of mixotrophy and light for growth and survival of the toxic haptophyte Prymnesium parvum . Harmful Algae 10 : 388 - 394 Cadet J , Wagner JR , Shafirovich V , Geacintov NE ( 2014 ) One-electron oxidation reactions of purine and pyrimidine bases in cellular DNA . Int J Radiat Biol 90 : 423 - 432 Casas-Monroy O , Chan P-S , Linley RD , Vanden Byllaardt J , Kydd J , Bailey SA ( 2016 ) Comparison of three techniques to evaluate the number of viable phytoplankton cells in ballast water after ultraviolet irradiation treatment . J Appl Phycol 28 : 2821 - 2830 Cerino F , Zingone A ( 2006 ) A survey of cryptomonad diversity and seasonality at a coastal Mediterranean site . Eur J Phycol 41 : 363 - 378 Chaudhari VR , Vyawahare A , Bhattacharjee SK , Rao BJ ( 2015 ) Enhanced excision repair and lack of PSII activity contribute to higher UV survival of Chlamydomonas reinhardtii cells in dark . Plant Physiol Biochem 88 : 60 - 69 Cochran WG ( 1950 ) Estimation of bacterial densities by means of the most probable number . Biometrics 6 : 105 - 116 Coesel S , Mangogna M , Ishikawa T , Heijde M , Rogato A , Finazzi G , Todo T , Bowler C , Falciatore A ( 2009 ) Diatom PtCPF1 is a new cryptochrome/photolyase family member with DNA repair and transcription regulation activity . EMBO Rep 10 : 655 - 661 Cullen JJ , MacIntyre HL ( 2016 ) On the use of the serial dilution culture method to enumerate viable phytoplankton in natural communities of plankton subjected to ballast water treatment . J Appl Phycol 28 : 279 - 298 Davidson K , Gowen RJ , Tett P , Bresnan E , Harrison PJ , McKinney A , Milligan S , Mills DK , Silke J , Crooks AM ( 2012 ) Harmful algal blooms: how strong is the evidence that nutrient ratios and forms influence their occurrence? Estuar Coast Shelf Sci 115 : 399 - 413 Doblin MA , Popels LC , Coyne KJ , Hutchins DA , Cary SC , Dobbs FC ( 2004 ) Transport of the harmful bloom alga Aureococcus anophagefferens by oceangoing ships and coastal boats . Appl Environ Microbiol 70 : 6495 - 6500 Falkowski PG , Owens TG , Ley AC , Mauzerall DC ( 1981 ) Effect of growth irradiance levels on the ratio of reaction centers in two species of marine phytoplankton . Plant Physiol 68 : 969 - 973 FDA ( 2010 ) Bacteriological analytical manual. Appendix 2: most probable number from serial dilutions . FDA . http://www.fda.gov/food/ foodscienceresearch/laboratorymethods/ucm109656.htm First MR , Drake LA ( 2013 ) Approaches for determining the effects of UV radiation on microorganisms in ballast water . Manag Biol Invasion 4 : 87 - 99 First MR , Drake LA ( 2014 ) Life after treatment: detecting living microorganisms following exposure to UV light and chlorine dioxide . J Appl Phycol 26 : 227 - 235 First MR , Robbins-Wamsley SH , Riley SC , Drake LA ( 2017 ) Assessment of variable fluorescence fluorometry as an approach for rapidly detecting living photoautotrophs in ballast water . J Sea Res . https://doi.org/10.1016/j.seares. 2017 . 02 .012 Franklin DJ , Berges JA ( 2004 ) Mortality in cultures of the dinoflagellate Amphidinium carterae during culture senescence and darkness . Proc R Soc B 271 : 2099 - 2107 Frazier M , Miller AW , Lee H II , Reusser DA ( 2013 ) Counting at low concentrations: the statistical challenges of verifying ballast water discharge standards . Ecol Appl 23 : 339 - 351 Furuya K , Marumo R ( 1983 ) The structure of the phytoplankton community in the subsurface chlorophyll maxima in the western North Pacific Ocean . J Plankton Res 5 : 393 - 406 Geider RJ , La Roche J , Greene RM , Olaizola M ( 1993 ) Response of the photosynthetic apparatus of Phaeodactylum tricornutum (Bacillariophyceae) to nitrate, phosphate or iron starvation . J Phycol 29 : 755 - 766 Glibert PM , Burkholder JM , Kana TM ( 2012 ) Recent insights about relationships between nutrient availability, forms, and stoichiometry, and the distribution, ecophysiology, and food web effects of pelagic and benthic Prorocentrum species . Harmful Algae 14 : 231 - 259 Guan WC , Lu SH ( 2010 ) The short- and long-term response of Scrippsiella trochoidea (Pyrrophyta) to solar ultraviolet radiation . Photosynthetica 48 : 287 - 293 Guillard RRL ( 1975 ) Culture of phytoplankton for feeding marine invertebrates . In: Smith WL , Chanley MH (eds) Culture of marine invertebrate animals . Plenum , New York, pp 108 - 132 Guillard RRL , Hargraves PE ( 1993 ) Stichochrysis immobilis is a diatom, not a chrysophyte . Phycologia 32 : 234 - 236 Haji Malayeri A , Mohseni M , Cairns W , Bolton JR ( 2016 ) Fluence (UV dose) required to achieve incremental log inactivation of bacteria, protozoa, viruses and algae . International Ultraviolet Association Inc News 18 ( 3 ): 4 − 6 and supplemental tables Han MS , Kim YP , Cattolico RA ( 2002 ) Heterosigma akashiwo (Raphidophyceae) resting cell formation in batch culture: strain identity versus physiological response . J Phycol 38 : 304 - 317 Harris ASD , Jones KJ , Lewis J ( 1998 ) An assessment of the accuracy and reproducibility of the most probable number (MPN) technique in estimating numbers of nutrient stressed diatoms in sediment samples . J Exp Mar Biol Ecol 231 : 21 - 30 Harrison PJ , Thompson PA , Calderwood GS ( 1990 ) Effects of nutrient and light limitation on the biochemical composition of phytoplankton . J Appl Phycol 2 : 45 - 56 Heijde M , Zabulon G , Corellou F , Ishikawa T , Brazard J , Usman A , Sanchez F , Plaza P , Martin M , Falciatore A , Todo T , Bouget F-Y , Bowler C ( 2010 ) Characterization of two members of the cryptochrome/photolyase family from Ostreococcus tauri provides insights into the origin and evolution of cryptochromes . Plant Cell Environ 33 : 1614 - 1626 Henri V , Helbronner A , de Recklinghausen M ( 1910 ) Nouvelles recherches sur la stérilization de grandes quantités d'eau par les rayons ultraviolets . Comp Rend Hebd Seances l'Acad Sci 151 : 677 - 680 Hijnen WAM , Beerendonk EF , Medema GJ ( 2006 ) Inactivation credit of UV radiation for viruses, bacteria and protozoan (oo)cysts in water: a review . Water Res 40 : 3 - 22 Hull NM , Isola MR , Petri B , Chan P-S , Linden KG ( 2017 ) Algal DNA repair kinetics support culture-based enumeration for validation of ultraviolet disinfection ballast water treatment systems . Environ Sci Technol Lett . https://doi.org/10.1021/acs.estlett.7b00076 Hurley MA , Roscoe ME ( 1983 ) Automated statistical-analysis of microbial enumeration by dilution series . J Appl Bacteriol 55 : 159 - 164 Imai I , Itoh K , Terada K , Kamizono M ( 1986 ) Distribution of dormant cells of Chattonella (Raphidophyceae) and occurrence of summer red tide in Suo-Nada, western Seto Inland Sea . Bull Jpn Soc Sci Fish 52 : 1665 - 1671 IMO ( 2004 ) International convention for the control and management of ballast water and sediments . International Maritime Organization. http://www.imo.org/home.asp IMO ( 2008 ) Guidelines for approval of ballast water management systems (G8) . International Maritime Organization. http://www.imo. org/home.asp IMO ( 2016 ) Review of the guidelines for approval of ballast water management systems (G8): analysis methods for determining the viability of organisms in the 10 to 50 μm size class, submitted by Denmark and Norway . PPR 4/7 , 12 October 2016 . International Maritime Organization Sub-Committee on Pollution Prevention and Response IMO ( 2017 ) Report to the Marine Environment Protection Committee PPR\04\PPR 4-21 . International Maritime Organization SubCommittee on Pollution Prevention and Response 4th Session Itakura S , Imai I , Itoh K ( 1997 ) BSeed bank^ of coastal planktonic diatoms in bottom sediments of Hiroshima Bay, Seto Inland Sea, Japan . Mar Biol 128 : 497 - 508 Jensen JN ( 2010 ) Disinfection model based on excess inactivation sites: implications for linear disinfection curves and the Chick-Watson dilution coefficient . Environ Sci Technol 44 : 8162 - 8168 Jochem FJ ( 1990 ) On the seasonal occurrence of autotrophic naked nanoflagellates in Kiel Bight, western Baltic . Estuar Coast Shelf Sci 31 : 189 - 202 Kacser H , Burns JA ( 1973 ) The control of flux . Symp Soc Exp Biol 27 : 65 - 104 Kana TM , Glibert PM ( 1987 ) Effect of irradiances up to 2000 μE m−2 s−1 on marine Synechococcus WH7803-I. Growth, pigmentation, and cell composition . Deep-Sea Res 34 : 479 - 495 Knight-Jones EW ( 1951 ) Preliminary studies of nanoplankton and ultraplankton systematics and abundance by a quantitative culture method . J Cons Int Explor Mer 17 : 140 - 155 Kolber Z , Zehr J , Falkowski P ( 1988 ) Effects of growth irradiance and nitrogen limitation on photosynthetic energy conversion in photosystem II . Plant Physiol 88 : 923 - 929 Lee YK ( 2001 ) Microalgal mass culture systems and methods: their limitation and potential . J Appl Phycol 13 : 307 - 315 Lepetit B , Dietzel L ( 2015 ) Light signaling in photosynthetic eukaryotes with 'green' and 'red' chloroplasts . Environ Exp Bot 114 : 30 - 47 Li AS , Stoecker DK , Adolf JE ( 1999 ) Feeding, pigmentation, photosynthesis and growth of the mixotrophic dinoflagellate Gyrodinium galatheanum . Aquat Microb Ecol 19 : 163 - 176 Liebich V , Stehouwer PP , Veldhuis M ( 2012 ) Re-growth of potential invasive phytoplankton following UV-based ballast water treatment . Aquat Invasions 7 : 29 - 36 MacIntyre HL , Cullen JJ ( 2005 ) Using cultures to investigate the physiological ecology of microalgae . In: Andersen RA (ed) Algal culture techniques . Academic, New York, pp 287 - 326 MacIntyre HL , Cullen JJ ( 2016 ) Classification of phytoplankton cells as live or dead using the vital stains fluorescein diacetate and 5- chloromethylfluorescein diacetate . J Phycol 52 : 572 - 589 MacIntyre HL , Kana TM , Anning T , Geider RJ ( 2002 ) Photoacclimation of photosynthesis irradiance response curves and photosynthetic pigments in microalgae and cyanobacteria . J Phycol 38 : 17 - 38 Manoharan K , Lee T , Cha J , Kim J , Lee W , Chang M , Park C , Cho J ( 1999 ) Acclimation of Prorocentrum minimum (Dinophyceae) to prolonged darkness by use of an alternative carbon source from triacylglycerides and galactolipids . J Phycol 35 : 287 - 292 Martínez F , Orús MI ( 1991 ) Interactions between glucose and inorganic carbon metabolism in Chlorella vulgaris UAM 101 . Plant Physiol 95 : 1150 - 1155 McCrady MH ( 1915 ) The numerical interpretation of fermentation-tube results . J Infect Dis 17 : 183 - 212 McKie-Krisberg ZM , Gast RJ , Sanders RW ( 2015 ) Physiological responses of three species of Antarctic mixotrophic phytoflagellates to changes in light and dissolved nutrients . Microb Ecol 70 : 21 - 29 McQuoid MR ( 2002 ) Pelagic and benthic environmental controls on the spatial distribution of a viable diatom propagule bank on the Swedish west coast . J Phycol 38 : 881 - 893 Miller JN , Miller JC ( 2005 ) Statistics and chemometrics for analytical chemistry, 5th edn . Pearson Education Limited , Harlow NWRI ( 2012 ) Ultraviolet disinfection guidelines for drinking water and water reuse, Third edn . National Water Research Institute, Fountain Valley Oemcke D , Parker N , Mountfort D ( 2004 ) Effect of UV irradiation on viability of micro scale and resistant forms of marine organisms: implications for the treatment of ships' ballast water . J Mar Environ Eng 7 : 153 - 172 Oliveri P , Fortunato AE , Petrone L , Ishikawa-Fujiwara T , Kobayashi Y , Todo T , Antonova O , Arboleda E , Zantke J , Tessmar-Raible K , Falciatore A ( 2014 ) The cryptochrome/photolyase family in aquatic organisms . Mar Genomics 14 : 23 - 37 Olsen RO , Hess-Erga O-K , Larsen A , Thuestad G , Tobiesen A , Hoell IA ( 2015 ) Flow cytometric applicability to evaluate UV inactivation of phytoplankton in marine water samples . Mar Pollut Bull 96 : 279 - 285 Othman R , Wollum AG II ( 1987 ) Evaluation of some methods for enumeration of nitrogen-fixing blue-green algae in soil . J Elisha Mitchell Sci Soc 103 ( 3 ): 94 - 100 Popels LC , MacIntyre HL , Warner ME , Zhang YH , Hutchins DA ( 2007 ) Physiological responses during dark survival and recovery in Aureococcus anophagefferens (Pelagophyceae) . J Phycol 43 : 32 - 42 Rastogi R , Richa KA , Tyagi M , Sinha R ( 2010a ) Molecular mechanisms of ultraviolet radiation-induced DNA damage and repair . J Nucleic Acids 2010 : 592980 Rastogi RP , Singh SP , Haeder D-P , Sinha RP ( 2010b ) Detection of reactive oxygen species (ROS) by the oxidant-sensing probe 2′,7′- dichlorodihydrofluorescein diacetate in the cyanobacterium Anabaena variabilis PCC 7937 . Biochem Biophys Res Commun 397 : 603 - 607 Rastogi RP , Singh SP , Hader DP , Sinha RP ( 2011 ) Ultraviolet-B-induced DNA damage and photorepair in the cyanobacterium Anabaena variabilis PCC 7937 . Environ Exp Bot 74 : 280 - 288 Redford EL , Myers J ( 1951 ) Some effects of ultraviolet radiations on the metabolism of Chlorella . J Cell Comp Physiol 38 : 217 - 243 Ren ZJ , Zhang L , Shi Y ( 2016 ) UV disinfection of Staphylococcus aureus in ballast water: effect of growth phase on the disinfection kinetics and the mechanization at molecular level . Sci China Technol Sci 59 ( 2 ): 330 - 336 Romero-Martínez L , Moreno-Andrés J , Acevedo-Merino A , Nebot E ( 2016 ) Evaluation of ultraviolet disinfection of microalgae by growth modeling: application to ballast water treatment . J Appl Phycol 28 : 2831 - 2842 Roy S , Parenteau M , Casas-Monroy O , Rochon A ( 2012 ) Coastal ship traffic: a significant introduction vector for potentially harmful dinoflagellates in eastern Canada . Can J Fish Aquat Sci 69 : 627 - 644 Russell SM ( 2000 ) Comparison of the traditional three-tube Most Probable Number method with the Petrifilm, SimPlate, BioSys optical, and Bactometer conductance methods for enumerating Escherichia coli from chicken carcasses and ground beef . J Food Prot 63 : 1179 - 1183 Sakai H , Katayama H , Oguma K , Ohgaki S ( 2011 ) Effect of photoreactivation on ultraviolet inactivation of Microcystis aeruginosa . Water Sci Technol 63 : 1224 - 1229 Segovia M , Berges JA ( 2005 ) Effect of inhibitors of protein synthesis and DNA replication on the induction of proteolytic activities, caspaselike activities and cell death in the unicellular chlorophyte Dunaliella tertiolecta . Eur J Phycol 40 : 21 - 30 Sinha RP , Hader DP ( 2002 ) UV-induced DNA damage and repair: a review . Photochem Photobiol Sci 1 : 225 - 236 . Sokal RR , Rohlf FJ ( 2011 ) Biometry, 4th edn . W.H Freeman and Company, New York Stehouwer PP , Buma A , Peperzak L ( 2015 ) A comparison of six different ballast water treatment systems based on UV radiation, electrochlorination and chlorine dioxide . Environ Technol 36 : 2094 - 2104 Takahashi A , Shibata N , Nishikawa S , Ohnishi K , Ishioka N , Ohnishi T ( 2006 ) UV-B light induces an adaptive response to UV-C exposure via photoreactivation activity in Euglena gracilis . Photochem Photobiol Sci 5 : 467 - 471 Thomas MK , Kremer CT , Klausmeier CA , Litchman E ( 2012 ) A global pattern of thermal adaptation in marine phytoplankton . Science 338 : 1085 - 1088 Throndsen J ( 1978 ) The dilution-culture method . In: Sournia A (ed) Phytoplankton manual . UNESCO , Paris, pp 218 - 224 Throndsen J , Kristiansen S ( 1991 ) Micromonas pusilla (Prasinophyceae) as part of picoplankton and nanoplankton communities of the Barents Sea . Polar Res 10 : 201 - 207 Tsolaki E , Diamadopoulos E ( 2010 ) Technologies for ballast water treatment: a review . J Chem Technol Biotechnol 85 : 19 - 32 USCG ( 2012 ) Standards for living organisms in ships' ballast water discharged in US waters . Fed Regist 77 : 17254 - 17320 . United States Coast Guard USEPA ( 2006 ) Ultraviolet disinfection guidance manual for the final long term. 2. Enhanced surface water treatment rule (EPA 815 -R- 06 - 007) . Vol EPA 815 -R- 06 -007. United States Environmental Protection Agency Office of Water (4601) , Washington DC USEPA ( 2010 ) Generic protocol for the verification of ballast water treatment technology . vol Version 5 .1. EPA/600/R-10/146. United States Environmental Protection Agency Environmental Technology Verification Program , Washington van de Poll WH , Hanelt D , Hoyer K , Buma AGJ , Breeman AM ( 2002 ) Ultraviolet-B-induced cyclobutane-pyrimidine dimer formation and repair in arctic marine macrophytes . Photochem Photobiol 76 : 493 - 500 Walsby AE , Juttner F ( 2006 ) The uptake of amino acids by the cyanobacterium Planktothrix rubescens is stimulated by light at low irradiances . FEMS Microbiol Ecol 58 : 14 - 22 Wang BL , Wang X , Hu YW , Chang MX , Bi YH , Hu ZY ( 2015 ) The combined effects of UV-C radiation and H2O2 on Microcystis aeruginosa, a bloom-forming cyanobacterium . Chemosphere 141 : 34 - 43 Weavers LK , Wickramanayake GB ( 2001 ) Kinetics of the inactivation of microorganisms . In: Block SS ( ed) Disinfection, sterilization, and preservation, 5th edn . Lippincott Williams & Wilkins , Philadelphia, pp 65 - 78 Wood MA , Everroad RC , Wingard LM ( 2005 ) Measuring growth rates in microalgal cultures . In: Andersen RA (ed) Algal culture techniques . Academic, New York, pp 269 - 285 Wright DA , Welschmeyer NA ( 2015 ) Establishing benchmarks in compliance assessment for the ballast water management convention by port state control . J Mar Eng Technol 14 : 9 - 18 Wright DA , Welschmeyer NA , Peperzak L ( 2015 ) Alternative, indirect measures of ballast water treatment efficacy during a shipboard trial: a case study . J Mar Eng Technol 14 : 1 - 8 Yentsch CS , Reichert CA ( 1963 ) The effects of prolonged darkness on photosynthesis, respiration and chlorophyll in the marine flagellate Dunaliella euchlora . Limnol Oceanogr 8 : 338 - 342 Zetsche EM , Meysman FJR ( 2012 ) Dead or alive? Viability assessment of micro- and mesoplankton . J Plankton Res 34 : 493 - 509

This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs10811-017-1254-8.pdf

Hugh L. MacIntyre, John J. Cullen, Trina J. Whitsitt, Brian Petri. Enumerating viable phytoplankton using a culture-based Most Probable Number assay following ultraviolet-C treatment, Journal of Applied Phycology, 2017, 1-22, DOI: 10.1007/s10811-017-1254-8