Herbicides are a class of chemical pesticides that are widely used to remove or control nuisance plants in agriculture and horticulture . Residues amounting to 99.7% of the applied load are dispersed and enter the aquatic environment through run-off and leaching with damaging consequences for ecosystem health .
Paraquat, N,N-dimethyl-4,4-bipyridinium dichloride, is one of the most widely used ionic herbicides for controlling the growth and spread of broadleaved weeds, grasses and aquatic weeds. Contamination of water courses and lakes leads to its rapid accumulation by aquatic organisms, and especially fish . It is a systemic herbicide that is highly toxic to plants, as it diverts electrons away from the donor side of the photosystem I (PSI) complex by accepting electrons from iron-sulphur (Fe-S) centres and/or ferredoxin, thus preventing electron transfer to NADP . Oxidation of paraquat generates highly phytotoxic reactive oxygen species (ROS) such as superoxide anion (O− 2), hydroxyl radicals (·OH), and hydrogen peroxide (H2O2) . Paraquat is also potentially lethal to a wide variety of non-target organisms, particularly primary producers  . There is growing concern over the exposure to small quantities of paraquat in the aquatic environment, as it may have negative impacts on the endocrine systems of humans and wildlife. Therefore, continued monitoring of the levels and persistence of paraquat in the environment is urgently needed . Despite its use being prohibited in several countries (e.g. throughout Europe) it continues to be manufactured and exported to various nations, globally . As weeds become resistant to the world’s most popular weed killer, Roundup, paraquat has been marketed as an alternative, with seven million pounds of paraquat being used in the United States on nearly 15 million hectares (https://www.nytimes.com/2016/12/20/business/paraquat-weed-killer-pesticide.html). A recent report suggests that paraquat be closely related to a less immediately apparent effect—Parkinson’s disease .
Therefore, determination of the presence of paraquat is still important for conducting a risk assessment of aquatic ecosystems. Chemical analysis of the environmental matrix is the most direct approach to reveal the status of contamination, of which reversed-phase high performance liquid chromatography (HPLC) with ultraviolet (UV) detection is considered a particularly sensitive technique for the determination of paraquat levels. However, conventional chemical analyses have several drawbacks, such as complex procedures of sample preparation and the need for expensive chemicals and equipment. In addition, most chemical analyses do not provide ecologically relevant information on the temporal changes or interactive effects of pollutants . Moreover, paraquat is usually used in small quantities (1.0 - 5.0 ppm) and, therefore, the concentrations of paraquat residues in water are extremely low .
Aquatic bioassays are an important means of assessing the quality of water containing pollutants (both mixtures of pollutants and unknown pollutants) and are useful for providing safety standards for water management in an ecological context. Many different species of plants and protists have been used as test organisms, but primary producers are of paramount importance for monitoring the functioning and health of ecosystems. Since the 1940s, duckweed (Lemnaceae) has been extensively used as a model organism for conducting fundamental and applied research in environmental sciences, particularly phytotoxicity testing and bioremediation, and has also been used as a bioresource   . Duckweeds are commonly found in freshwater and brackish ecosystems in temperate climates and serve not only as an important food source for various water birds and fish, but also as habitats for small invertebrates. Most ecotoxicological studies conducted to date focus on duckweeds belonging to the genus Lemna, and in particular L. gibba and L. minor .
Recently, a toxicity test based on growth (change in area) inhibition of 3 day-old Spirodela fronds germinated from turions has recently been proposed for the International Organization for Standardization (ISO 20227) . During its life cycle, Spirodela produces turions (meaning “shoot”), a type of buds that come from modified shoot apices and are often rich in starch and sugar, so that they can act as storage organs (Figure 1). Turions are known to be produced in response to adverse conditions such as decreasing day length or falling temperature . Turions sink to the bottom of a pond or lake and hibernate when the water freezes, but rise again to germinate and grow into complete plants in spring.
The Spirodela test had similar levels of sensitivity to the Lemna test (ISO 20079) for nine herbicides, four inorganic and organic compounds and nine metals . Results from a study by Oláh et al.  indicated differential sensitivity of turions and mature fronds to three metals (cadmium (Cd), chromium (Cr) and nickel (Ni)), with higher tolerance levels for turions than normal fronds. Thus, further investigation of the relative sensitivity of plants of different ages and at different life stages is required.
Figure 1. The lifecycle of Spirodela polyrhiza. This plant forms a dormant stage, turion, during the dormant phase (a); after artificial germination of the turion sank to the bottom (b); germinated turoin forms young fronds (c) and these are cultivated under optimal conditions (d); then they grow back in the mature stage (e) and repeat their lifecycle.
The appropriate selection of the assessed endpoints in (phyto) toxicity tests is the key factor for a successful and valid test result. Conventional toxicity testing methods (ISO20079) have employed numerous endpoints including: frond number, plant number, root number, dry or fresh biomass, root length, frond diameter, carbon uptake and chlorophyll (Chl) content . Chlorophyll a fluorescence of photosystem II (PS II) reaction centres, a rapid and sensitive tool for evaluating toxicity in algae and higher plants  , has also been successfully used in Lemna tests .
Recently, Park et al.  reported a well-defined toxicant concentration-dependent inhibition of root regrowth for three Lemna species. This method has several operational advantages over other conventional techniques (ISO20079); for example, it requires a shorter duration (72 h), smaller test solution volume (approximately 3.0 mL) and use of non-axenic plant material.
Depending on the circumstances, fronds germinated from turions may be stored in the laboratory for a long time before testing. It is therefore important to ascertain whether frond samples retrieved from long-term storage are suitable material for tests. The prospective ERA for herbicides has for many years been based exclusively on Lemna sp. as the only aquatic macrophyte, although it has only relatively recently been extended to dicotyledonous submerged macrophytes such as Myriophyllum sp. The inclusion of new monocotyls such as Spirodela as a test battery for herbicide toxicity may therefore not be an overpaid attempt. In this study, we aimed to: 1) compare the intrinsic physiological traits and paraquat sensitivity of 3-day-(young) and 8-week-old (mature) S. polyrhiza fronds germinated from turions, based on measurements of frond area and fresh biomass weight, root regrowth, pigment (Chl a and Chl b) content and Chl a fluorescence, 2) verify the inter-relationship between the tested endpoints in response to paraquat exposure; 3) evaluate whether extremely low concentrations of paraquat in aquatic ecosystems exert toxicological action on S. polyrhiza; and 4) evaluate whether the endpoints employed were sensitive enough to detect the toxic action of paraquat.
2. Methods and Materials
2.1. Sample Collection and Maintenance
Mature plants of Spirodela polyrhiza (L.) Schleiden were collected from the Upo wetland located at Changnyeong-gun, Gyeongsangnam-do, Korea (35˚34'05''60N, 128˚24'06''90E). Fronds were maintained in a 1.5 L glass tank containing Steinberg medium  at 25˚C under continuous white light (30 - 40 µmol photons m−2∙s−1) supplied by cool daylight fluorescence tubes (FL 20 SS/18D, Philips). The medium was replaced regularly at 7 day intervals, and its pH was adjusted to 6.9 ± 0.2 using either 1 M hydrochloric acid (HCl) or 1 M sodium hydroxide (NaOH).
To induce the production of turions, 1.0 L culture flasks were stored at 15˚C for 4 - 6 weeks. Upon induction, each turion was separated from its mother frond . All turions were then collected from the bottom of the flasks and transferred to 50 mL Falcon tubes containing Steinberg medium. The Falcon tubes were covered with foil, and the stock seed banks were stored at 4˚C until needed for tests.
2.2. Toxicity Test
Prior to the tests, turions were germinated in Petri dishes (85.6 mm × 12.6 mm) containing Steinberg medium at 25˚C under continuous light (photon irradiance of 90 - 100 µmol photons m−2∙s−1). Fronds germinated from turions were maintained under the same environmental conditions for either 3 days or 8 weeks. Just prior to experimentation, roots were cut from all fronds using stainless steel scissors and an individual rootless plant was added to each well of a 24-well plastic plate (85.4 × 127.6 mm; well diameter, 15.6 mm; SPL, Seoul, Korea). There were 4 plants per paraquat concentration, 6 concentrations per plate and three replicate plates (n = 3). Different nominal concentrations of paraquat (Table 1) were prepared by diluting the original stock solutions (CAS No. 1910-42-5) with Steinberg medium. To perform the toxicity tests, a 3.0 mL test solution was added to each well. Tests were run for 72 h in an environmentally controlled chamber maintained at 25 ± 1˚C and continuous light (100 ± 10 µmol photons m−2∙s−1). Static tests were performed, i.e. the test solutions were not renewed during the period of exposure. A fully randomized design was used to account for any variability in environmental conditions within the culture chamber.
2.3. Physiological Traits of Intrinsic Fronds and Responses to Paraquat Exposure
2.3.1. Frond Area, Weight and Regrown Root Length
Fronds maintained for 3 days or 8 weeks with no other experimental treatments were randomly selected for comparison of their intrinsic physiological traits. Following the 72 h exposure to different concentrations of paraquat, S. polyrhiza fronds with regrown roots were picked using tweezers and placed upside-down on a glass slide. Because the fronds are wet, roots could be easily straightened with a light touch. Photographs of the fronds with regrown roots were captured using image analyser (e.g. Image J); the length of the longest root was measured to evaluate the root regeneration ability of fronds following paraquat exposure. The justification for why only the longest root to measure was to see the response of the fronds that will regenerate their first root becoming longest.
Table 1. Final concentration range and mode of action used for testing toxicity of paraquat with 2 ages of Spirodela polyrhiza.
Fronds were then dried on a paper towel and weighed to determine their fresh biomass.
2.3.2. Chl Contents
Chl a and b contents (μg∙mL−1) were calculated using the following equation :
where A666 and A653 represent the absorbance at 666 and 653 nm. Based on the Chl content of the extract, the Chl content per gram of frond fresh biomass was calculated.
2.3.3. Chl a Fluorescence
Chl a fluorescence was measured using an Imaging PAM (Walz, Germany), as a proxy for photosynthetic performance. To measure Fv/Fm and rETRmax, samples were incubated in the dark for 10 - 15 min. Pulses (0.15 µmol photons m−2∙s−1) from light emitting diodes (LED) were used to determine the initial fluorescence yield (Fo), which denotes the fluorescence yield when all photosystem II (PSII) reaction centres are open, with fully oxidized plastoquinone A (QA). Then, a saturation pulse of approximately 5000 µmol photons m−2∙s−1 emitted by an LED lamp was applied to produce the maximum fluorescence yield (Fm), which is induced by a short saturating pulse of actinic light that reduces all QA molecules. The value of Fv/Fm was then calculated using the following equation:
where Fv is the variable fluorescence.
Rapid light curves were produced using 10s pulses of actinic light increased stepwise from 0 to 1517 μmol photons m−2∙s−1 . Effective quantum yield (ΦPSII) was calculated using the following equation:
where is the maximum light-acclimated fluorescence yield, and F is the light-acclimated fluorescence yield.
The maximum electron transport rate (ETRmax) was calculated using the hyperbolic tangent equation adapted from Jassby and Platt :
where α and I indicate ETR and PFD, respectively, under light limiting conditions.
2.4. Statistical Analysis
Data were analyzed using one-way analysis of variance (ANOVA) at a significance level of P < 0.05, after homogeneity test. To determine differences among treatments, post-hoc comparisons were performed using the least significant difference (LSD) test. Toxicity test results were reported as EC50 values (effective concentration at which 50% inhibition occurs), with 95% confidence intervals estimated using the linear interpolation method (ToxCalc 5.0; Tidepool Science, California, USA).
3. Results and Discussion
3.1. Turion Germination
Determination of the optimal conditions for turion germination is important since several environmental factors can simultaneously influence germination of the turions of S. polyrhiza , but standardization of optimal germination for toxicity testing has not yet been established . The germination of S. polyrhiza under different photon irradiances, pH and temperatures is shown in Figure 1. Turions germinated in all light levels from 5 to 100 μmol photons m−2∙s−1 but not in the dark (Figure 2(a)); turion germination is light dependent normally . The maximal germination was recorded at pH 5 and 7 but at pH 3 no germination was observed (Figure 2(b)). There was 100% germination of turions at 25˚C and 30˚C but there was no evidence of germination at either 5˚C or 15˚C (Figure 2(c)).
In nature, germination of Spirodela turions is controlled by temperature and light , although exposure to combinations of red and far-red light, long photoperiods, gibberellin and kinetin (plant hormones) can accelerate the germination process   . In laboratory studies, turion germination was found to be controlled by exposure to temperature and light of 22˚C - 25˚C and 60 to 120 μmol photons m−2∙s−1 for 3 - 7 days  .
It has been observed that turions float to the water surface by formation of internal oxygen bubbles prior to their germination  . Photosynthetic activities thus appear to be involved in the germination process since turion surfacing may be geared by photosynthetic O2 production, which is also affected by bicarbonate concentrations as a major aquatic carbon source . In this respect, temperature and pH would be influential factors for initiation of germination since the amount of carbon (as CO2 or H2CO3) in the water is dependent on these
Figure 2. Effects of photon flux density (PFD) (a); pH (b) and temperature (c) on the germination of turion in Spirodela polyrhiza. The bar denotes mean values and error bar indicates the standard deviation. LSD shows the least significant difference at the 5% level. No bar indicates no germination.
parameters. It is therefore interesting to note that germination reaches its maximum at pH 7 where the availability of H2CO3 is at its greatest.
3.2. Standard Conditions
Environmental factors influence growth and photosynthetic activity of aquatic plants, with pH, light and temperature being probably the most important parameters. It has been shown that these parameters also modify the toxic effects of chemical pollutants in aquatic plants . Thus, it is important that the culture conditions under which toxicity tests are carried out are fully reported/disclosed and that the tests are, preferably, carried out under standardized conditions for a particular species and endpoints to be measured. For this reason, we first determined the optimal conditions for growth of fronds and roots of S. polyrhiza (Figure 3 and Figure 4). The maximal growth of both frond and root was recorded at 100 μmol photons m−2∙s−1 and 25˚C while no significant difference was observed in the range of pH between 4 and 10. It is interesting to note that there was some growth at low temperatures despite the lack of germination at these temperatures. Also, no growth was detected at temperatures lower than 10˚C. All these data appear to imply the existence of different control mechanisms for turion germination, frond growth and root regrowth.
Figure 3. Effects of photon flux density (PFD) (a); pH (b) and temperature (c) on the relative growth rate (%, d−1) of germinated fronds in Spirodela polyrhiza. The bar denotes mean values and error var indicates the standard deviation. LSD shows the least significant difference at the 5% level. No bar indicates that all samples are dead.
Figure 4. Effects of photon flux density (PFD) (a); pH (b) and temperature (c) on the root regrowth length (mm) of germinated fronds in Spirodela polyrhiza. The bar denotes mean values and error var indicates the standard deviation. LSD shows the least significant difference at the 5% level. No bar indicates that all samples are dead.
3.3. Intrinsic Physiological Traits of Young and Mature Fronds in S. polyrhiza
Different aged fronds (3 days and 8 weeks old) showed intrinsic differences in physiological traits (Table 2). The area and fresh weight of young fronds were approximately 19.06% and 29.35% that of mature fronds. However, roots of young fronds (17.58 ± 2.98 mm) were longer than those of mature fronds (11.98 ± 3.23 mm). No statistically significant differences were detected in Chl a and Chl b contents between young and mature fronds. When Spirodela fronds become mature, there are enlargements in the mesophyll cells with prominent intercellular spaces and air chambers in the abaxial surface of maturing fronds . Morphological changes in the internal structure of fronds during development may explain the significantly larger frond area and weight in mature fronds than in newly germinated fronds.
Spirodela plants have simple, short and thin roots that lack branches and hairs , although the root systems of young and mature fronds have not been compared. Roots appear to serve an important role in balancing and maintaining the floating body in a stable and upright position . Differences in root length between young and mature fronds might be explained by the higher photosynthetic capacity of young fronds. Thus, in young plants photosynthetic products are re-directed in favor of root production that provides the smaller and lighter floating fronds stable positioning.
On the other hand, we found significantly higher Chl a and Chl b content in young than in mature fronds, which was similar to the report made by Kim . As the chlorophyll contents were calculated on a fresh weight basis, lower chlorophyll in mature fronds may be ascribed to higher proportion of non-photosynthetic structures such as intracellular spaces and air chambers.
Values of Fv/Fm, rETRmax and IK of young fronds (0.752 ± 0.034, 20.58 ± 8.28 and 122.34 ± 47.64 µmol photons m−2∙s−1, respectively) were significantly higher than those of mature fronds (0.491 ± 0.160, 12.18 ± 6.355 and 71.22 ± 26.28 µmol
Table 2. Physiological traits of young and mature fronds of S. polyrhiza.
FA; relative growth rate of frond area, FW; fresh weight, RL; root regrowth length, Chl a; chlorophyll a content, Chl b; chlorophyll b content, Fv/Fm; Optimal quantum yield, rETRmax; the maximum electron transport rate, SD; standard deviation. **P < 0.001; *P < 0.05.
photons m−2∙s−1, respectively). However, alpha values were similar between the two frond types (0.168 for young and 0.160 for mature fronds).
The maximum quantum yield of PSII chemistry (Fv/Fm) is an estimate of the photochemical conversion efficiency of PSII in the dark. Higher values of Fv/Fm indicate higher light utilization efficiency and greater ability to adapt to low light conditions . Lemna spp. typically have Fv/Fm maxima of about 0.68 - 0.73  . Compared with Lemna spp., young fronds of S. polyrhiza were found to have similar Fv/Fm, but mature fronds showed lower Fv/Fm. Fv/Fm values appear to be in line with chlorophyll contents in that Fv/Fm was greater in young fronds with higher Chl a and Chl b than in mature fronds.
rETR is an empirical estimate of the rate of flow of electrons through the electron transport chain and rETRmax is defined as the maximum rETR at saturating irradiance . Higher rETRmax in young fronds means that there is higher probability of a photochemical event resulting in more active and effective electron transport upon absorption of a photon by the antennae of PS II. The initial slope of a P vs I curve represents the light-limited phase of photosynthesis, and its gradient is affected by the efficiency with which the plant can absorb the limited light available . Saturation irradiance (Ik) is defined as the point at which the extrapolated initial slope crosses rETRmax. Comparison of rETRmax, alpha and Ik showed that the number of photosynthetic units (PSUs) but not the size of PSUs differed between young and mature fronds of S. polyrhiza since ETRmax and Ik values were higher while alpha values were the same . Young fronds seem to have more PSUs than mature ones, and hence the higher rETRmax and Ik values.
3.4. Comparative Sensitivity of Young and Mature Fronds to Paraquat
Paraquat sensitivity of S. polyrhiza fronds at both stages was dependent on the endpoints. After a 72 h exposure to paraquat, the most sensitive paraquat toxicity endpoints were root length and rETRmax in both frond types, with EC50 values of 0.66 ± 0.03 and 0.76 ± 0.09 µg∙L−1, respectively, in young fronds and 5.53 ± 1.5 and 2.28 ± 0.39 µg∙L−1, respectively, in mature fronds (Table 3). Coefficients of variation for root re-growth and rETRmax were 4.55% and 11.84%, respectively, for young fronds and 27.12% and 17.11%, respectively, for mature fronds (Table 3). Paraquat was reported to be one of the most toxic herbicides for Lemna spp. with EC50 values of 7.1 - 10.6 µg∙L−1 for root regrowth and 6.6 - 8.0 µg∙L−1 for rETRmax . The paraquat sensitivity of Spirodela was higher in young fronds and similar in mature fronds to that of Lemna spp.
It is notable that there was statistically significant difference in the sensitivity to paraquat toxicity between young and mature fronds, which would bring to question as to which stages, young or mature fronds would then be more suitable for toxicity testing. Sensitivity of bioassay methods is important for determination of whether to use them or not for water quality risk assessment. Effective bioassays should produce results within the relevant environmental ranges.
Table 3. EC50 and CV values for inhibition of various parameters in two different ages (young and mature) of S. polyrhiza exposed to different concentrations of paraquat. Mean ± standard deviation are shown (n = 3 plates, 24 plants per plate with 4 plants per each concentration).
FA; relative growth rate of frond area, FW; fresh weight, RL; root regrowth length, Chl a; chlorophyll a content, Chl b; chlorophyll b content, Fv/Fm; Optimal quantum yield, rETRmax; the maximum electron transport rate, EC50; the effective concentration at which 50% inhibition occurs, CV; coefficient of values, NA; not applicable.
A current guideline for the allowable concentrations in drinking waters set by WHO is 10 µg∙L−1 for paraquat. This study shows that endpoints of root length and rETRmax in both young and mature frond of S. polyrhiza are sensitive enough to detect toxic impacts of water samples containing paraquat in excess of allowable guidelines. However, the sensitivity of young fronds is almost over 10 times the permissible levels of paraquat in drinking water, and unless even a slightest presence of paraquat should be monitored and prevented, mature fronds would more successfully be employed for management decisions. Spirodela methods also show a high level of precision and reproducibility which are essential for adoption of toxicity testing methods. A desirable level of repeatability expressed by CVs is 30% or less according to Environment Canada (2007) . For young and mature fronds with the endpoints of root length and rETRmax CVs for EC50 values were found to lie within this acceptable range (Table 3).
The ultimate goal of bioassay testing is to provide representative and inclusive criteria of exposure conditions, thereby improving risk assessment and water quality management. In this respect, multiple rather than single endpoint assays may have greater potential for a more comprehensive risk assessment of toxicants. Such an approach allows to gain important insights into the mechanisms of toxicity and to obtain information on the relative sensitivity of measured endpoints to toxicity concentration and/or exposure duration, thereby identifying specific endpoints that can effectively detect interferences caused by certain phytotoxics . In the past, little attention has been paid to the roots in Lemna since it was generally considered that root fragility made their handling for measurements difficult and that it was impractical to obtain sufficient numbers of individual plants with identical root lengths to initiate tests. However, more recently the ecotoxicological significance of the root endpoint has been re-evaluated and root length is now considered to be a sensitive, precise and ecologically significant endpoint in comparison with more traditional frond growth and biomass endpoints  .
One of the most frequently used methods for monitoring the status of the photosynthetic apparatus in plants is in vivo chlorophyll a fluorescence, a non-destructive, straightforward and rapid technique that is applicable in both laboratory and field studies. It is used as a potential indicator of exposure to environmental and chemical stresses, including herbicides. The impact of certain herbicides, such as commonly used ones like diuron, atrazine and simazine, on the photochemical activity of PSII has long been recognized    . Such PSII inhibitors restrict photosynthetic activity through their binding to the D1 protein in thylakoids and blocking electron transport through the PSII reaction center with effects being manifested by changes in various chlorophyll fluorescence parameters . Frond area, fresh weight, root length, chlorophylls and chlorophyll a fluorescence have never been measured simultaneously as endpoints. FigureS1 shows correlative relationships between two different endpoints. There was a strong relationship between the frond area and weight and between Chl a and Chl b to a lesser degree in both young and mature fronds. It was notable that root length rather than frond growth was highly correlated with photosynthetic performance. Photosynthetic electron transport events support the biochemical reactions needed for plant growth since the electron transport rate is closely related to the photosynthetic activity including oxygen evolution or CO2 uptake . Therefore, a direct or an indirect effect of a pollutant on photosynthetic processes is observed prior to an effect on the growth process . Little relationship between PSII inhibition and frond growth may indicate possible detoxification of cells and recovery of growth, thereby the inhibition of PSII inhibition no longer reflecting the inhibition of growth . In contrast, the close relationship between PSII inhibition and root re-growth may imply a direct or an indirect effect of a pollutant on photosynthetic processes observed prior to an effect on the root growth process .
We found intrinsic physiological traits difference between young and mature frond of S. polyrhiza: young fronds invest more photosynthetic products to develop root elongation and chlorophylls which might be an adaptive strategy for future development via active photosynthesis and stabilization of the body positioning using longer roots, whereas mature fronds tilt to increase frond area and weight via structural development. For organic xenobiotics in particular, such as paraquat and other artificial toxicants, there is evidence that some standardised endpoints, such as total growth, fresh or dry weight or number of individuals, may under- or overestimate the actual risk level, especially if the doses of the toxic compounds are sublethal. In many cases there are other valuable endpoints such as more detailed analyses of stress responses at the subcellular and biochemical levels, but these assessments are still quite expensive because of the probes and detection systems or laboratory equipment required and are not included in standardised toxicity tests. In this respect, our Spirodela root test with either turion germinated or mature fronds would be a valuable asset for the risk assessment of paraquat, as it is a technique that allows toxicity to be assessed more quickly and easily, but without loss of sensitivity.
This research was supported by a grant from Incheon National University Research (2015-1444).
Figure S1. Matrix view of Pearson’s correlations for different parameters in young (a) and mature (b) Spirodela polyrhiza. The color and size of the circles indicate the strength of the correlation. The color bar on the right provides a scale of the correlation coefficients. Light blue, blue and navy colors mean the correlation coefficients are over 0.5, 0.7, and 0.9, respectively, while the light pink, orange and red color indicates correlation coefficients < 0.5. The lable within the circle refers to the correlation coefficients values. An asterisk (*) refers to the statistical significance at the 5% level. Note that for clear presentation, all matrix views show only the upper portion of the matrix to avoid duplication. FA; relative growth rate of frond area, FW; fresh weight, RL; root regrowth length, Chl a; chlorophyll a content, Chl b; chlorophyll b content, Fv/Fm; Optimal quantum yield, rETRmax; the maximum electron transport rate.
 Fatima, M., Mandiki, S., Douxfils, J., Silvestre, F., Coppe, P. and Kestemont, P. (2007) Combined Effects of Herbicides on Biomarkers Reflecting Immune-Endocrine Interactions in Goldfish: Immune and Antioxidant Effects. Aquatic Toxicology, 81, 159-167.
 Prado, R., García, R., Rioboo, C., Herrero, C., Abalde, J. and Cid, A. (2009) Comparison of the Sensitivity of Different Toxicity Test Endpoints in a Microalga Exposed to the Herbicide Paraquat. Environment International, 35, 240-247.
 Gabryelak, T. and Klekot, J. (1985) The Effect of Paraquat on the Peroxide Metabolism Enzymes in Erythrocytes of Freshwater Fish Species. Comparative Biochemistry and Physiology—Part C: Toxicology & Pharmacology, 81, 415-418.
 Moustaka, J. and Moustakas, M. (2014) Photoprotective Mechanism of the Non-Target Organism Arabidopsis thaliana to Paraquat Exposure. Pesticide Biochemistry and Physiology, 111, 1-6.
 Franqueira, D., Cid, A., Torres, E., Orosa, M. and Herrero, C. (1999) A Comparison of the Relative Sensitivity of Structural and Functional Cellular Responses in the Alga Chlamydomonas eugametos Exposed to the Herbicide Paraquat, Archives of Environmental Contamination and Toxicology, 36, 264-269.
 Mensah, P.K., Palmer, C.G. and Muller, W.J. (2014) Lethal and Sublethal Effects of Pesticides on Aquatic Organisms: The Case of a Freshwater Shrimp Exposure to Roundup®. In: Pesticides-Toxic Aspects, IntechOpen, London, 163-185.
 Zaranyika, M. and Nyoni, S. (2013) Degradation of Paraquat in the Aquatic Environment: A Proposed Enzymatic Kinetic Model That Takes into Account Adsorption/Desorption of the Herbicide by Colloidal and Sediment Particles. International Journal of Research in Chemistry and Environment, 3, 26.
 Vaccari, C., El Dib, R. and de Camargo, J.L.V. (2017) Paraquat and Parkinson’s Disease: A Systematic Review Protocol According to the OHAT Approach for Hazard Identification. Systematic Reviews, 6, 98-98.
 Park, J., Brown, M.T., Depuydt, S., Kim, J.K., Won, D.-S. and Han, T. (2017) Comparing the Acute Sensitivity of Growth and Photosynthetic Endpoints in Three Lemna Species Exposed to Four Herbicides. Environmental Pollution, 220, 818-827.
 Iqbal, J., Javed, A. and Baig, M.A. (2019) Growth and Nutrient Removal Efficiency of Duckweed (Lemna minor) from Synthetic and Dumpsite Leachate under Artificial and Natural Conditions, PLoS ONE, 14, e0221755.
 Sońta, M., Rekiel, A. and Batorska, M. (2019) Use of Duckweed (Lemna L.) in Sustainable Livestock Production and Aquaculture—A Review. Annals of Animal Science, 19, 257-271.
 Drost, W., Matzke, M. and Backhaus, T. (2007) Heavy Metal Toxicity to Lemna minor: Studies on the Time Dependence of Growth Inhibition and the Recovery after Exposure. Chemosphere, 67, 36-43.
 Baudo, R., Foudoulakis, M., Arapis, G., Perdaen, K., Lanneau, W., Paxinou, A.-C., Kouvdou, S. and Persoone, G. (2015) History and Sensitivity Comparison of the Spirodela polyrhiza Microbiotest and Lemna Toxicity Tests. Knowledge and Management of Aquatic Ecosystems, 416, 23.
 Appenroth, K.J. and Nickel, G. (2010) Turion Formation in Spirodela polyrhiza: The Environmental Signals That Induce the Developmental Process in Nature. Physiologia Plantarum, 138, 312-320.
 Oláh, V., Hepp, A. and Mészáros, I. (2015) Comparative Study on the Sensitivity of Turions and Active Fronds of Giant Duckweed (Spirodela polyrhiza (L.) Schleiden) to Heavy Metal Treatments. Chemosphere, 132, 40-46.
 Schreiber, U., Quayle, P., Schmidt, S., Escher, B.I. and Mueller, J.F. (2007) Methodology and Evaluation of a Highly Sensitive Algae Toxicity Test Based on Multiwell Chlorophyll Fluorescence Imaging. Biosensors and Bioelectronics, 22, 2554-2563.
 Muller, R., Schreiber, U., Escher, B.I., Quayle, P., Nash, S.M.B. and Mueller, J.F. (2008) Rapid Exposure Assessment of PSII Herbicides in Surface Water Using a Novel Chlorophyll a Fluorescence Imaging Assay. Science of the Total Environment, 401, 51-59.
 Hulsen, K., Minne, V., Lootens, P., Vandecasteele, P. and Hofte, M. (2002) A Chlorophyll a Fluorescence-Based Lemna minor Bioassay to Monitor Microbial Degradation of Nanomolar to Micromolar Concentrations of Linuron. Environmental Microbiology, 4, 327-337.
 Appenroth, K.J., Palharini, L. and Ziegler, P. (2013) Low-Molecular Weight Carbohydrates Modulate Dormancy and Are Required for Post-Germination Growth in Turions of Spirodela polyrhiza. Plant Biology, 15, 284-291.
 Lichtenthaler, H.K. and Wellmurn, A.R. (1983) Determinations of Total Carotenoids and Chlorophylls a and b of Leaf Extracts in Different Solvents. Biochemical Society Transactions, 11, 591-592.
 Schreiber, U. (2004) Pulse-Amplitude-Modulation (PAM) Fluorometry and Saturation Pulse Method: An Overview. In: Chlorophyll a Fluorescence, Springer, Berlin, 279-319.
 Jassby, A.D. and Platt, T. (1976) Mathematical Formulation of the Relationship between Photosynthesis and Light for Phytoplankton. Limnology and Oceanography, 21, 540-547.
 Appenroth, K.-J., Opfermann, J., Hertel, W. and Augsten, H. (1989) Photophysiology of Turion Germination in Spirodela polyrhiza (L.) SCHLEIDEN. II. Influence of After-Ripening on Germination Kinetics. Journal of Plant Physiology, 135, 274-279.
 Appenroth, K.-J., Dathe, W., Hertel, W. and Augsten, H. (1991) Photophysiology of Turion Germination in Spirodela polyrhiza (L.) SCHLEIDEN. VII. Action of Jasmonic Acid. Journal of Plant Physiology, 138, 345-349.
 Oláh, V., Hepp, A. and Mészáros, I. (2016) Assessment of Giant Duckweed (Spirodela polyrhiza L. Schleiden) Turions as Model Objects in Ecotoxicological Applications. Bulletin of Environmental Contamination and Toxicology, 96, 596-601.
 Park, J. and Dinh, T.B. (2019) Contrasting Effects of Monochromatic LED Lighting on Growth, Pigments and Photosynthesis in the Commercially Important Cyanobacterium Arthrospira maxima. Bioresource Technology, 291, Article ID: 121846.
 Schmid, R. and Dring, M.J. (1992) Circadian Rhythm and Fast Responses to Blue Light of Photosynthesis in Ectocarpus (Phaeophyta, Ectocarpales). Planta, 187, 53-59.
 Nestler, H., Groh, K.J., Schonenberger, R., Eggen, R.I. and Suter, M.J.-F. (2012) Linking Proteome Responses with Physiological and Biochemical Effects in Herbicide-Exposed Chlamydomonas reinhardtii. Journal of Proteomics, 75, 5370-5385.
 Gopalapillai, Y., Vigneault, B. and Hale, B.A. (2014) Root Length of Aquatic Plant, Lemna minor L., as an Optimal Toxicity Endpoint for Biomonitoring of Mining Effluents. Integrated Environmental Assessment and Management, 10, 493-497.
 Beaumont, G., Bastin, R. and Therrien, H. (1976) Physiological Effects of Sublethal Doses of Atrazine on Lemna minor L. II. Effect on Photosynthesis and on Respiration. Naturaliste Canadien, 103, 535-541.
 Merlin, G., Eulaffroy, P. and Blake, G. (1993) Use of Fluorescence Induction Kinetics of Lemna minor as a Tool for Chemical Stress Evaluation. Science of the Total Environment, 134, 761-772.
 Küster, A. and Altenburger, R. (2007) Development and Validation of a New Fluorescence-Based Bioassay for Aquatic Macrophyte Species. Chemosphere, 67, 194-201.
 Kumar, K.S. and Han, T. (2011) Toxicity of Single and Combined Herbicides on PSII Maximum Efficiency of an Aquatic Higher Plant, Lemna sp. Toxicology and Environmental Health Sciences, 3, 97-105.
 Murata, N., Takahashi, S., Nishiyama, Y. and Allakhverdiev, S.I. (2007) Photoinhibition of Photosystem II under Environmental Stress. Biochimica et Biophysica Acta (BBA)—Bioenergetics, 1767, 414-421.
 Holbrook, G.P., Beer, S., Spencer, W.E., Reiskind, J.B., Davis, J.S. and Bowes, G. (1988) Photosynthesis in Marine Macroalgae: Evidence for Carbon Limitation. Canadian Journal of Botany, 66, 577-582.
 Juneau, P., Sumitomo, H., Matsui, S., Itoh, S., Kim, S.-G. and Popovic, R. (2003) Use of Chlorophyll Fluorescence of Closterium ehrenbergii and Lemna gibba for Toxic Effect Evaluation of Sewage Treatment Plant Effluent and Its Hydrophobic Components. Ecotoxicology and Environmental Safety, 55, 1-8.
 Hayat, Q., Hayat, S., Alyemeni, M. and Ahmad, A. (2012) Salicylic Acid Mediated Changes in Growth, Photosynthesis, Nitrogen Metabolism and Antioxidant Defense System in Cicer arietinum L. Plant, Soil and Environment, 58, 417-423.