Accumulation, morpho-physiological and oxidative stress induction by single and binary treatments of fluoride and low molecular weight phthalates in Spirodela polyrhiza L. Schleiden

The present study examined the interactive effects of fluoride and phthalates on their uptake, generation of reactive oxygen species and activation of antioxidative defence responses in Spirodela polyrhiza L. Schleiden. A hydroponic study was conducted in which S. polyrhiza cultured in Hoagland’s nutrient medium, was exposed to fluoride (50 ppm) and different concentrations viz., 75, 150 300 ppm of diethyl phthalate (DEP) and diallyl phthalate (DAP) individually as well as in combination for the time period of 24, 72, 120 and 168 h respectively. A significant decline in fresh weight, dry to fresh weight ratio, total chlorophyll, carotenoid content and increased anthocyanin content was observed. Fluoride and phthalates was found to be readily accumulated by S. polyrhiza in all the exposure periods. Interestingly, when binary treatments were given in nutrient medium, uptake of both fluoride and phthalate was found to be influenced by each other. In combined treatments, DEP stimulated fluoride uptake, while its own uptake was restricted by fluoride. In contrary to this, fluoride stimulated DAP uptake. Moreover, combined stress further caused significant decrement in carbohydrate, protein content and increment in MDA levels, phenolic content and electrolyte leakage. Nevertheless, phthalates showed more pronounced oxidative stress and growth inhibition compared to fluoride. To cope up with the oxidative damage, enhanced level of antioxidant enzymatic activities was observed in S. polyrhiza under both fluoride and phthalate stress as compared to control. Scanning electron microscope imaging of leaf stomata revealed that combined stress of fluoride with phthalates caused distortion in the shape of guard cells. Confocal micrographs confirmed the generation of reactive oxygen species, cell damage, disruption in membrane integrity, and enhanced levels of glutathione in plant cells. This study focussed on ecotoxicological and interactive significance of fluoride led phthalate uptake or vice versa which was also assumed to confer tolerance attributes.

. Accumulation of DEP, DAP and fluoride content by S. polyrhiza at different exposure periods.ou. Data is shown as (mean ± SE, two-way ANOVA, Tukey's HSD), * indicates significant at (p ≤ 0.05). ND = not detected. Mean values with same letters signifies no significant difference between two values of different concentrations for the same exposure period (using one-way ANOVA).
www.nature.com/scientificreports www.nature.com/scientificreports/ Effect on photosynthetic pigments. In the present investigation, single and binary treatments of DEP and fluoride induced significant decrease in total chlorophyll content as compared to the respective controls during all the durations, but binary treatments of fluoride with 300 ppm DEP showed minimum chlorophyll content. Also, individual as well as combined treatment of fluoride with DAP showed decreasing trend in total chlorophyll content at all the durations and the maximum content was found after 168 h of treatment period (Fig. 2). Negative β-coefficients indicated that treatment of fluoride and phthalate decreased the chlorophyll content to some extent, while negative βcoefficients for interaction of F x phthalate depicted further decrement in the content ( Table 2).
On the other hand, carotenoid content in plant showed decreasing trend at all the applied concentrations of DEP, DAP and fluoride and the minimum carotenoid content was observed at combined concentration of fluoride with 300 ppm DEP as well as DAP after the exposure period of 72 hours (Fig. 2).
All individual treatments with either fluoride, DEP and DAP exhibited significant increase in anthocyanin pigent when compared to control plants (Fig. 2). However, almost all binary treatments led to significantly higher anthocyanin content in comparison to respective single treatment of fluoride, DEP and DAP concentrations at all the exposure periods. Notably, combined effect of fluoride with DEP showed more anthocyanin accumulation as compared to DAP. www.nature.com/scientificreports www.nature.com/scientificreports/ Effect on soluble protein, carbohydrate and MDA content. All the individual treatments either with fluoride, DEP and DAP showed significant reduction (p ≤ 0.05) in protein content at all the exposure periods (Fig. 3). However, binary treatment of 50 ppm of fluoride and 75 ppm of DEP exhibited high value of protein content (1.20 mg/gfw) at 24 h time period, while at other exposure periods of 72, 120 and 168 h, binary treatments of fluoride and DEP showed slight increased protein content as compared to maximum individual treatment (300 ppm) of DEP. Notably, maximum reduced protein content (0.47 mg/gfw and 0.22 mg/gfw) was found when plants were exposed to binary treatments of 50 ppm of fluoride and 300 ppm of DEP and DAP for 168 h respectively. Interestingly, combined treatment of fluoride and DAP induced marked reduction in protein content in all the combinations (Fig. 3). Negative β-coefficients for fluoride and phthalate revealed reduction in protein content in plant, while interaction between fluoride and phthalate was positively regressed for protein content at the exposure periods of 72, 120 and 168 h ( Table 2).    www.nature.com/scientificreports www.nature.com/scientificreports/ 103.99 mg/gfw, 101.19 mg/gfw at 24 h time period as compared to control (T2). At other exposure periods of 72, 120 and 168 h, the significant reduction of 35, 17.8, 18.7% respectively was found at combined treatments of fluoride with 300 ppm DEP. However, exposure of plant to individual DAP and combined treatment of fluoride and DAP induced significant reduction and showed dose dependent relationship at all the exposure periods. Maximum decrease of 49.59, 66.2, 64.16 and 85.8% in carbohydrate content was found at 24, 72, 120 and 168 h respectively at highest combined treatment of fluoride and DAP. Multiple linear regression also revealed that F-, DEP and DAP treatment significantly reduced carbohydrate content, whereas binary treatment showed slight alterations (Table 2).
Single and joint effects of fluoride, DEP and DAP were also measured in terms of significant enhancement in MDA content shown in Fig. 3. Exposure of plant to individual and binary treatments provoked dose-responsive enhancement in MDA content. Maximum significant increase (p ≤ 0.05) in MDA content was recorded as 146.5, 374.1, 303.1 and 203.1% at binary treatments of fluoride with 300 ppm DEP whereas, combination of fluoride with highest concentration of DAP exhibited maximum enhancement of 131.7, 447.7, 453.6 and 197.9% after exposure period of 24, 72, 120 and 168 h respectively. Positive β regressions implied increment in MDA content. Interaction of fluoride with phthalate was positively regressed implying further increase in MDA content during combined stress ( Table 2).
Effect on total phenolic, proline content and electrolyte leakage. Present investigation revealed that individual treatment of fluoride, DEP and DAP led to significant increase in phenolic content, when plants of S. polyrhiza were exposed for 24, 72, 120 and 168 h time period (Fig. 4). However, content was maximum in combination of fluoride with highest concentration of either DEP or DAP as compared to control. Maximum percent increase was recorded as 160. 6 (Table 2).
Effect on enzymatic activities. All applied treatments to experimental plants exhibited significantly higher SOD activity when compared to control, which was mainly pronounced in combined treatments of fluoride with DEP and DAP during all the durations (Fig. 5). Highest SOD activity was recorded at the binary treatment of fluoride with 300 ppm DEP after 24 h time period in comparison to individual treatments. With increasing time interval, reduced activity was observed in binary treatments. However, the treatment of fluoride with highest concentration of DAP showed maximum SOD activity after the time interval of 120 and 168 hours as compared to single treatments.
All individual treatments significantly elevated CAT activity in exposed fronds as compared to control during all the exposure periods, although all binary treatments triggered even higher CAT activity. Maximum CAT activity was observed at the highest combination of fluoride with DEP as well as DAP during the time interval of 24, 72 and 120 h, while after the exposure period of 168 h, reduction in CAT activity was found.
In the present investigation, APX activity showed significant increase with increasing dose and duration. Maximum activity was observed when binary treatments of fluoride with 300 ppm DEP were given to S. polyrhiza after the time period of 168 h, while combined exposure of fluoride with 300 ppm DAP showed highest APX activity after 120 h time period. However, binary treatments induced reduced CAT activity after 24 and 72 h time period when compared to control and individual treatments of DAP.
In S. polyrhiza, enhanced GPOX activity was measured upon exposure to individual fluoride, DEP and DAP concentrations and their combinations comparison to control. Binary treatment of fluoride with DEP showed reduced GPX activity as compared to individual treatment after 24 hours and then become maximum at the highest combination of fluoride with 300 ppm DEP. While, reduced GPX activity was observed at same combined treatment after 72, 120 and 168 h. However, individual treatment of DAP as well as combined treatment of DAP with fluoride showed increasing trend at all the durations and maximum activity was found after the time period of 120 h.
In this observation, all individual treatments induced significant increase in GR activity, while binary treatment of fluoride with 75 ppm DEP or DAP showed decreased GR activity when compared to control during all the durations. However, maximum increase in the activity was found at highest combined concentration of fluoride with 300 ppm DEP after 168 h, while those treated with combination of fluoride with 300 ppm DAP showed higher GR activity after the exposure period of 24 h.

Discussion
Plants play a pivotal role as a bioindicator to give us insights into cumulative and interactive effects of environmental contaminants. Deleterious effects on them give us alarming signals for protecting our environment and human beings. In our study, we reported the effects of binary concentrations of fluoride and phthalates on accumulation, growth and biochemical parameters of S. polyrhiza. The potential of aquatic plants for uptake of organic compounds and removal of nutrients from wastewater has been vastly explored 16,17 . The present investigation on S. polyrhiza showed that DEP, DAP and fluoride tend to accumulate in fronds of an aquatic macrophyte, which is in accordance with the data reported for other aquatic as well as terrestrial plants 18,19 . As illustrated in Table 1, accumulation of DEP and DAP in the plants decreased with the increase in the exposure periods. Higher DEP concentrations in growth medium led to enhanced DEP uptake which then gradually decreased with increasing exposure time, while uptake of DAP decreased with increasing concentrations though the accumulated content was much higher than DEP. Although majority of studies reported the individual effects of fluoride and www.nature.com/scientificreports www.nature.com/scientificreports/ phthalates on various plant species, but there is no report about the fluoride-phthalate interaction mechanism in any plant species. Since, the laboratory experiments are performed under controlled conditions, so the results obtained partly give us an idea about the consequences under environmental conditions, which are much more complex and unstable. From these results, it was revealed that fluoride and phthalates considerably influenced each other uptake, when present together in the nutrient medium. In fronds exposed to binary treatments of fluoride with DEP, higher fluoride content was observed as compared to accumulated DEP content, while binary treatment with DAP, higher DAP uptake was observed, which indicated that fluoride stimulated DAP uptake. This suggested synergistic relationship existing between fluoride and phthalates. However, the effect of fluoride on alleviation of DEP toxicity might be attributable partly to the fluoride induced decreased accumulation of DEP in fronds. Different responses in accumulation and interaction of fluoride with DEP and DAP attributed to genotype of plant, contaminant, pH of the medium, temperature and humidity 21,22 . Interestingly, uptake of organic compounds by this free floating macrophyte is driven by simple process of diffusion and then penetration into a leaf as solutions 23,24 . Since, plants do not have any specific transporters for these organic compounds, so their movements into the plants mainly depends on their physicochemical properties such as hydrophobicity, polarity, aqueous solubility and molecular weight of the organic contaminant 24,25 . Our results revealed that longer exposure periods resulted in less accumulation of phthalates in plants, which might be due to metabolic pathways for the biotransformation of organic contaminants adopted by the plant. Exposure of aquatic plants to organic compounds resulted in i) fast uptake or sequestration of the compound into vacuoles ii) Transformation or degradation of the compound via volatilization, lignification or metabolization to carbon dioxide and water iii) Assimilation into plant tissues as non-toxic compounds 24,26 .
In our experiment, we found that accumulation of F, DEP and DAP in fronds of an aquatic macrophyte was accompanied by considerable cellular and metabolic changes, ultimately enhancing tolerance ability of plant. One of the visible change was growth inhibition of plant due to toxic effects of added F, DEP and DAP in the medium after all the exposure periods. Present study demonstrated duration-dependent response of plant against binary www.nature.com/scientificreports www.nature.com/scientificreports/ as well as individual treatments stress on growth parameters of plant. Noticeably, maximum effect was observed after longer duration combined treatments of fluoride with DEP or DAP. Reduction in growth and dry to fresh weight ratio under toxic conditions might be due to disruption in metabolic pathways and adverse effect on cell division, cellular growth or enhanced accumulation rate of both fluoride and phthalate in combined treatments 27 . Also, there may be decreased transport of nutrients from the nutrient medium to the plant which caused growth retardation. Results are in coherence with studies that reported the effects of fluoride and phthalates on biomass reduction in S. polyrhiza 28,29 . Combined effects of fluoride and phthalates further inhibited plant growth compared to single treatment of F − , DEP and DAP.
A reduction in total chlorophyll and carotenoid content was observed in both individual and combined treatments of F − , DEP and DAP during all the exposure periods. Influence of applied treatments of DEP and DAP on photosynthetic pigments is a consequence of nutrient deficiency such as iron or altered activity of δ-ALA-D www.nature.com/scientificreports www.nature.com/scientificreports/ (delta-aminolevulunic acid dehydratase) that induced inhibition of chlorophyll biosynthesis 30 . Moreover, another reason behind the reduction of photosynthetic pigments may be due to lipid peroxidation in the membrane of chloroplasts 31 . Perusal of literature reported reduction in photosynthetic content in various aquatic plants exposed to stress such as Lemna gibba and Lemna minor 32,33 . However, fronds exposed to DAP showed highest amount of total chlorophyll content during 168 h time period as compared to other exposure periods, suggesting that individual as well combined stress induced disruption or disorganization of thylakoid membranes and activation of chlrophyllase enzyme led to enhanced pigment content 34 . Moreover, we observed a significant decrease in carotenoids content at combined stress of F with DEP as well as DAP, implying their role in inducing inhibition of biosynthesis of carotenoids. As imposed treatments with both DEP and DAP did not induce much reduction in carotenoid content as compared to control, suggesting minimal effect on carotenoids. Similar reports are there which revealed the toxic effects of anthracene (PAH) on photosynthetic pigments of Lemna gibba 35 . In addition, anthocyanin pigment content significantly enhanced in this investigation, implying stress activated the gene which stimulated anthocyanin biosynthesis. Though, the effect of DEP on anthocyanin pigment was negligible as compared to DAP in this observation. Increased content is a part of strategy adopted by the plant to minimize the deleterious impacts of stress.
Similarly, evaluation of total soluble protein content revealed significant toxic effects of F − , DEP and DAP treatments on aquatic plant, which is an important parameter for determining phytotoxic effects of environmental contaminants. Highest reduction in protein content was recorded in fronds when exposed to combined concentration of fluoride and phthalates. Overall, DEP and DAP are toxic and significantly decreased protein content of S. polyrhiza, probably due to increased protease activity or other catabolic enzymes. Also, decreased level of protein content in response to fluoride and phthalates could be explained by the fact that toxicant formed during the combination of fluoride and phthalate might altered cytochrome oxidase activity and various respiratory pathways 36 . These results are consistent with earlier findings that both fluoride and phthalate provoked oxidative stress and lead to protein degradation in plant species 28 . Slight increase of protein content at some durations attributed to the alteration in protease activity to cope up with the stress. Moreover, in response to the applied treatments, carbohydrate showed decreasing trend which might help in the adjustment of osmotic imbalance by the conversion of starch to sugar in the presence of hydrolytic enzymes or to supply more ATP to the stressed cells via increasing respiration rate [37][38][39] . Although the individual as well as combined toxic effect of DAP on carbohydrate content was much more pronounced as compared to DEP at all the concentrations and exposure periods, implying toxic behaviour of DAP on plant.
On the other hand, single and combined effects of fluoride, DEP and DAP, measured in terms of formation of MDA content, byproduct of lipid peroxidation of plasma membrane, found to be increased in the present study, since, it is well documented that both fluoride and phthalate triggered oxidative stress, that damages cell membrane and consequently enhanced MDA content. Furthermore, this result is supported by the findings of Ting-Ting et al. 40 and Cai et al. 20 that reported increased MDA content in mung beans under phthalates (DBP and DEHP) stress and tea plant under fluoride stress respectively. Higher phthalate concentrations can induce oxidative stress, which has already been reported for several plant species.The mechanism behind the formation of reactive oxygen species (ROS) is a consequence of interaction of the contaminants with the lipid rich membrane, subsequently led to conformational changes in the membrane by the activation of NADPH-oxidase located there, hence generating ROS 41 .
To gain further insights into the osmotic adjustment potential of the plant, we evaluated proline content of the plant. Our data revealed that the accumulation of proline enhanced significantly during all the treatments, whilst plant tends to accumulate maximum proline contents at individual fluoride concentrations. Although, combined concentrations of fluoride with DAP and DEP showed more proline accumulation as compared to individual concentrations in our study. Being an osmoprotectant, proline contributes in the detoxification of ROS and protects integrity of membrane 42 . Enhanced proline accumulation is an indicator of stress in S. polyrhiza. This also confirmed a dynamic bond between increased lipid peroxidation and proline content in this observation. Our results are in coherence with earlier findings of impact of fluoride and phthalates on proline content in plants 40,43 . Additionally, we detected disruption in the membrane integrity in the form of ion leakage and found that to be maximum at combined concentrations of fluoride and phthalates. These results correlated well with the increased lipid peroxidation and leakage of ions during combined stress, suggesting additional evidence that expression of proline and MDA is important during severity of combined stress, as they play a crucial role in protecting cellular membrane from disruption, consequently counteract the effects of stress in S. polyrhiza. Apart from this, increased accumulation of phenolic content in the fronds attributed to the activation of hexose-monophosphate pathway along with the release of bound phenols in stressed fronds 44 .
Generally, reactive oxygen species (ROS) are inescapable byproducts of metabolic processes such as photosynthesis and respiration in plants 45,46 . However, environmental contaminants such as heavy metals, herbicides, pesticides and other xenobiotics pose stress in plants and overwhelmingly generate ROS and free radicals, consequently compromising various metabolic and cellular activities. Meanwhile, plants possess an armory of antioxidant enzymes to scavenge ROS and free radicals and to counteract the deleterious effects. While, activity of antioxidant enzymes varies, since it depends on plant species, concentration, and exposure duration of contaminants affecting the plant 42,47 . Activation of enzymatic activities (SOD, CAT, APX, GPOX and GR) play a crucial role in detoxification and thus, represent altered redox status of the cell under stress 8 .
Among them, SOD is an essential enzyme of antioxidant defense system, catalyses the diminution of superoxide radicals to molecular oxygen and hydrogen peroxide at a very high pace 48 . In our present investigation of evaluation of enzymes activity, we observed that the SOD activity in fronds under F − , DEP and DAP enhanced, with the maximum elevation found at combined stress. Surprisingly, SOD activity was recorded in plants under DAP stress enhanced significantly with increasing concentration and exposure period, while SOD activity under DEP stress reduces initially and then increases after longer exposure periods, though showed increasing trend with the proceeding concentrations. Probable reason behind increased SOD activity is attributable to the adverse effects of combined concentrations mediated formation of enormous superoxide radicals. Ting-Ting et al. 40 stated that, exposure to higher concentrations of DBP and DEHP also resulted in significant rise in SOD activity in mung bean seedlings. Karmakar et al. 49 also reported increased SOD activity in aquatic macrophytes under fluoride stress. Subsequent increase in the activity reflected that SOD enzyme was stimulated by the formation of superoxide radicals, to protect the plant from the toxicity of fluoride and phthalates.
Further, removal of hydrogen peroxide involves the activation of catalase enzyme (CAT) from antioxidant machinery. CAT activity showed increased activity with the highest activity was recorded at maximum combined concentrations of fluoride with DEP as well as DAP. Although the present investigation revealed increased CAT activity with progressive increase in dose and duration, but the activity reduced at the highest exposure period as compared to the CAT activity found after other initial exposure periods, suggesting that scavenging potential of this enzyme declined in plants exposed to stress for longer duration. Thus, it is obvious that combined stressinduced damage might have increased at high concentrations and durations. Also imbalance in dynamic equilibrium existing between ROS population and pace of elimination or scavenging via enzymes resulted in deleterious effects on plants. Results are in agreement with the previous reports of the authors who reported decreased CAT activity in plants under phthalate and fluoride stres 49,50 .
Despite of the distinct antioxidative enzymes involved in scavenging of ROS, POD also plays a chief role in the breakdown of hydrogen peroxide. APX is one of the types of POD that uses ascorbate as an electron donor to eliminate hydrogen peroxide, which is the first step in ascorbate-glutathione cycle 48,51 . Importantly, to cope up with oxidative stress, activities of SOD and APX need not to be required sufficiently high, but balance in the equilibrium existing between SOD and APX is needed 45 . Similar trend was displayed by APX as by SOD in this observation. Compared to control, a dose-duration dependent maximum increase by 88.2% at 168 h under maximum concentrations of fluoride with DEP, while combined treatment with DAP showed maximum activity at 120 h. A possible explanation underlying enhanced APX activity was probably increased level of peroxides in the chloroplast of plant and stimulation of APX activity provides resistance to stress. High response of APX to phthalate in plant was also well demonstrated by Liao et al. 51 . Authors also suggested the defensive role of APX exposed to stress by its striking role to remove peroxides 48 .
Meanwhile, GPOX another type of POD, exhibited increased activity initially under combined exposure of fluoride with DEP exposure after 24 h and then reduced drastically at other exposure periods. Contrary to this, combined treatments of fluoride with DAP showed enhanced GPX activity in concentration-duration dependent manner. Interestingly, plant showed its efficiency to up-regulate both types of peroxidases to counteract the effects of peroxides generated under stressful environments.
Likewise APX, GR is also one of the main components of ascorbate-glutathione cycle which is generally located in the chloroplast, but also found in mitochondria and cytosol 48,52 . GR plays an important role in maintaining reduced/oxidized glutathione (GSH/GSSG) ratio by catalysing the reduction of GSSG to GSH and thus assured efficient recovery of glutathione during oxidative damage 1,53 . Present observation showed that GR activity has increased with enhanced dosage of fluoride, DEP and DAP individually as well as combined, implying activation of ascorbate-glutathione cycle at a very high pace to scavenge ROS generated in chloroplast of stressed fronds. Since, maintenance of glutathione pool is essential so that it does not intervene in phytochelatin biosynthesis. Basically, phytochelatins play a critical role in sequesteration of various environmental contaminants and inactivate them via conjugates formation 54,55 . Hence, a crucial role is played by GR in the protection of plant chloroplasts. Our results are consistent with the earlier reports on the increased GR activity in plants under phthalate exposure 56,57 .
Oxidative stress faced by S. polyrhiza under fluoride and phthalates was confirmed by viewing electron micrographs of adaxial surface of stressed fronds which showed many closed stomata as compared to unstressed fronds www.nature.com/scientificreports www.nature.com/scientificreports/ (Fig. 6). S. polyrhiza attained DEP stress tolerance attributes by showing alterations in stomatal movements. Furthermore, our results are supported by confocal micrographs which depicted no or less fluorescence in the roots of control plants, while high fluorescence was observed in treated roots implying that S. polyrhiza accumulated reactive oxygen species and fluorescence helps in determining apoptosis (Fig. 7). A hypothetical schematic representation depicting detoxification strategy adopted by S. polyrhiza is presented in Fig. 8. conclusions and future prospectives. In this investigation, fluoride and phthalates influenced the uptake of each other and revealed synergistic and/or antagonistic effects during combined treatments. The greatest impact of fluoride-phthalates induced stress was recorded at the highest concentration of combined treatments. Analysis showed that individual treatment of fluoride and phthalates induced oxidative damage to some extent but combined treatment further triggered oxidative burst and was more prominent in initial exposure periods, which then alleviated with increase in exposure periods. This indicates that this aquatic plant may adapt to the combined stress with due course of time by activating antioxidant machinery, thus maintaining redox imbalance. It may be possible that these contaminants (Fluoride-phthalates) bind to the plant cell wall and are metabolized or sequestered by plant in the vacuoles and converted into non-toxic form. Higher DAP uptake mediated by fluoride in combined treatments led to more ecotoxicological effects on plant as compared to DEP. Delineation of fate of fluoride-phthalates within plant tissues is required for better understanding of the mechanisms of toxicity and its sequesteration within plant cells. Diethyl phthalate (99.0% purity, CAS: 84-66-2) and diallyl phthalate (97.0% purity, CAS: 131-17-9) were purchased from Hi-Media, Mumbai (India). Other chemicals used as ingredients of Hoagland's nutrient medium (3%) were of analytical grade.

Materials and
Harvesting of plant material was done after 24, 72, 120 and 168 hours. Plant material was thoroughly washed, blotted to remove the excess water and then stored in deep freezer (−80 °C). Plant material was then assayed for various growth and biochemical parameters.
Analysis of accumulated phthalate content in S. polyrhiza using reverse phase high performance liquid chromatography (RP-HPLC). From plant samples, phthalates were extracted by using the ultrasonication method of Ma et al. 59 with slight modifications. 500 mg of plant sample was extracted with 20 ml acetone followed by ultrasonication of 30 min. Hexane (20 ml) was added for further extraction and again sonication was done for 30-35 min. The extract obtained was filtered and evaporated to 1 ml at 40 °C in rotary evaporator. Acetonitrile (4 ml) was added to remaining extract to make up 5 ml volume, filtered using 0.22 µm syringe filter and then analysed using reverse phase high performance liquid chromatography (RP-HPLC). chromatographic conditions. Phthalates in our present study were analysed using Shimadzu reverse phase high performance liquid chromatography coupled with PDA detector (RP-HPLC-PDA) Nexera system (Shimadzu, Japan). The system was equipped with degasser, autosampler, solvent mixing unit, column oven, PDA detector and CBM-20A (system controller). Chromatographic column used for separation of phthalates was C 18 having dimensions of 150 × 4.6 mm with pore size of 5 µm. The mobile phase used for separation were acetonitrile (Solvent A) and HPLC grade water (Solvent B). The gradient elution system was initiated with 50% A and 50% B, followed by 40% B for 2.00 min., 30% B held for 4.00 min, then 35% B for 5.00 min, followed by 35% B for 8.00 min and ended at 20.01 min. Flow rate was kept at 0.85 ml/min and injection volume was 10 μL. The maximum pressure limit of the instrument was 15000 psi. Column temperature was 38 °C and wavelength was set at 232 nm. Phthalates in the samples were detected on the basis of comparing retention time of standard peaks at absorbance of 232 nm (Supplementary File S1).

Method validation.
A method was developed and validated for the estimation of phthalates using reverse phase-high performance liquid chromatography (RP-HPLC). The stock solution (500 mg/l) of DAP and DEP was prepared in HPLC grade acetonitrile and working solutions (5-400 mg/l) were prepared through serial dilution of stock. (2019) 9:20006 | https://doi.org/10.1038/s41598-019-56110-w www.nature.com/scientificreports www.nature.com/scientificreports/ The method validation was carried out using normalized guidelines of International Conference on Harmonisation of Technical Requirement for Registration of Pharmaceutical for Human Use 60 . Firstly, the calibration curve was obtained using peak area versus concentration of each phthalate. The analyzed validation parameters are precision, accuracy, limit of detection (LOD) and limit of quantification (LOQ). The precision is repeatability of a system which is expressed as percentage of relative standard deviation (% RSD) and was calculated from six replicates of phthalate. Whereas, accuracy is the proximity between the conventional value and the observed value and expressed as percentage recovery. The determination of LOD and LOQ was done on the basis of signal-to-noise ratio (S/N) which is 3.3 and 10 respectively 60 .
Analysis of accumulated fluoride content in S. polyrhiza using ion chromatography. Fluoride content accumulated by plant was measured by using the protocol of Zhou et al. 61 with slight modifications. Plant samples were oven dried at 105 °C for first 25 min and then dried at 80 °C for next 48 hours. About 1 g of oven dried plant sample was taken in crucible and pulverised properly. The plant material was further soaked in 3 ml of 10% Mg(NO 3 ) 2 and 1 ml of 10% NaOH for 30 min. The samples were then evaporated and incubated at 200 °C in an oven. After carbonisation, samples were ashed in muffle furnace for 6 hours at 550 °C. Crucibles containing ashed material was finally cooled and rinsed with millipore water. Rinses of these samples were combined and total volume was made upto 30 ml. Fluoride content in plant samples was then analysed using ion chromatography (Metrohm Ion Chromatography, Orion-960) with known concentration of fluoride solution for calibration 28 . Determination of growth parameters. For determining fresh weight, plant material was washed several times, dried on filter paper to remove excess water and then weighed. Plants were oven dried at 105 °C for first 20 min and then dried at 80 °C for two days to obtain constant dry weight. Percentage change in fresh weight and dry weight to fresh weight ratio was calculated.
Determination of photosynthetic pigments. Chlorophyll content was extracted according to protocol given by Arnon 62 . Plant material was crushed in acetone (80%) and the extract was centrifuged for 10 min. Absorbance of chl a and chl b was taken at 663 and 645 nm respectively using ELISA plate reader and concentrations were expressed in µg/ml. Carotenoid content was measured by the method of Lichtenthaler and Wellburn 63 . Absorbance of the extract was observed at 470 nm and the concentration was expressed in µg/ml. Anthocyanin content was assayed by the protocol given by Mancinelli 64 . 1 g of fresh plant material was crushed in 3 ml of acidified methanol (99 ml methanol: 1 ml HCl). Extract was centrifuged for 20 min at 12000 rpm at 4 °C. Supernatant was collected and incubated at 4 °C for 24 h. Absorbance was taken at 530 and 657 nm. Results were expressed as OD/gfw. Determination of soluble protein, carbohydrate and MDA content. The measurement of soluble protein content was done by homogenising the plant material in potassium phosphate buffer (pH = 7). Extract was centrifuged and protein content was determined according Bradford method 65 , where bovine serum albumin was used as a standard. Lipid peroxidation was determined by measuring malondialdehyde (MDA) content 66 . Homogenisation of plant material was done in 0.1% trichloroacetic acid (TCA) and centrifuged for 10 min. Reaction mixture containing supernatant, 0.5% thiobarbituric acid (TBA) and 20% TCA was heated at 95 °C in a boiling water bath for 30 min and then rapidly cooled on ice. After centrifugation at 10,000 rpm for 10 min, the absorbance of the supernatant was read at 532 nm and non-specific absorption was corrected by subtracting the absorbance value observed at 600 nm. The concentration was expressed as µmol/g fresh weight. For carbohydrate determination, anthrone method was used 67 . Glucose was used as a standard.
Determination of total phenolic, proline content and electrolyte leakage. Total phenolic content was determined according to a procedure given by Singleton and Rossi 72 . Results were expressed as mg/gfw. Proline content in plant was determined by the method given by Bates et al. 73 . Results were expressed as µmol/g tissue.
Electrolyte leakage was determined as per the method given by Dionisio-sese and Tobita 74 . Fronds of S. polyrhiza (0.2 g) were cut into small pieces and added in the test tubes filled with distilled water (10 ml). Tubes were incubated at 40 °C for 2 h. After cooling, initial electrical conductivity (EC1) was measured using conductivity meter. Then samples were autoclaved at 121 °C for 20 min to release all the electrolytes. Samples were then cooled to room temperature and final electrical conductivity was recorded (EC2). Percent electrolyte leakage was calculated using following equation: www.nature.com/scientificreports www.nature.com/scientificreports/ Scanning electron microscopy (SEM). Scanning electron microscopy was performed to detect the stomatal response under stressful conditions. The stomatal response of control and treated fronds of S. polyrhiza was observed by using scanning electron microscope (Carl Zeiss-EvoLS 10). For sample preparation, leaf samples were fixed overnight in 2.5% glutaraldehyde prepared in 0.1 M potassium phosphate buffer (PPB). The leaf samples were washed with distilled water and then further dehydrated in different ethanol series (30%, 50%, 70% and 90%) for 15 min individually. Adaxial surface of leaves were then placed on metal stubs The surface features of fronds were observed under the scanning electron microscope at a voltage of 15 KV and stomata were viewed under the resolution of 500-4000 nm.
Confocal laser scanning microscopic (CLSM) studies. For confocal microscopy, roots of control and treated plant samples were washed with distilled water and then treated with different fluorescent dyes i.e. Propidium iodide (PI), dichlorofluorescein diacetate (H2DCFDA) and monochlorobimane (MCB) respectively in dark for 10 min separately for the study of cell viability, detection of ROS and GSH levels occurred during stress in plants. A 50 μM solution of PI was prepared according to the protocol of Gutierrez-Alcala, et al. 75 Detection of ROS levels in the roots of S. polyrhiza was done by following the method given by Ortega-Villasante et al. 76 Localization of glutathione (GSH) in roots of S. polyrhiza was done by the method of Fricker and Meyer 77 . For sample preparation, roots were washed and placed on glass slide over a drop of water to prevent dehydration and then covered with cover slip to visualise under confocal laser microscope. For propidium iodide (PI), He-Ne gas laser was used to excite the electrons at wavelength of 535 nm, multiline argon gas laser was used for 2, 7 dichloroflurescein (DCF) to excite the electrons at the wavelength of 488 nm and for monochlorobimane (MCB) excitation wavelength of 380 nm was used.
Statistical analysis. All the experimentation was performed in triplicates and expressed as the mean ± standard error. The data was subjected to two-way analysis of variance (ANOVA) for assessing the effect of fluoride and phthalates on S. polyrhiza. Analysis of all the experiments were done using self coded software MS-Excel 2010. The Tukey's post hoc multiple comparison test was done at 0.05 level of significance for the comparisons against control values. Multiple linear regression was performed to determine the relationship between dependent and independent variables and was calculated using following equation: where y = dependent variable, a = y-intercept, b 1 = partial regression coefficient for x 1 on y eliminating effect of x 2 , b 2 = partial regression coefficient for x 2 on y eliminating effect of x 1, x 1 = independent variable (Fluoride), x 2 = independent variable (Phthalate). In order to explain the relative effects of independent variables (fluoride and phthalate) on dependent variable, unitless β-regression coefficients (β1, β2 and β3) were used. ethical approval. This article does not include any studies related with human subjects or animal models performed by any of the authors.