Target Screening of Plant Secondary Metabolites in River Watersby Liquid Chromatography Coupled to High Resolution Mass Spectrometry (Lc-hrms)

DOI: https://doi.org/10.21203/rs.3.rs-36114/v1

Abstract

Background:Substantial efforts have been made to monitor potentially hazardous anthropogenic contaminants in surface waters while for plant secondary metabolites(PSMs) almost no data on occurrence in the water cycle are available.These metabolitesenter river waters through various pathways such as leaching, surface run-off and rain sewers or input of litter from vegetation andmight add to the biological activity of the chemical mixture.To reduce this data gap we conducted aLC-HRMS target screening in river waters from two different catchments for 150plant metabolites which were selected from a larger database considering their expected abundance in the vegetation, their potential mobility, persistence and toxicity in the water cycle and commercial availability of standards.

Results: The  screening revealed the presence of 12out of 150toxic PSMs including coumarins (bergapten, scopoletin, fraxidin, esculetin and psoralen), a flavonoid (formononetin) and alkaloids (lycorine and narciclasine). The compounds narciclasine and lycorine were detected at concentrations up to 3 µg/L while esculetin and fraxidin occurred at concentrations above 1 µg/L.Nine compounds occurred at concentrations above 0.1 µg/L, the Threshold for Toxicological Concern (TTC) for non-genotoxic and non-endocrine disrupting chemicals in drinking water.

Conclusions: Our study provides a first insight into potentially biologically active PSMs in surface waters and recommends their consideration in monitoring and risk assessment of water resources. This is currently hampered by a lack of effect data including toxicity to aquatic organisms, endocrine disruption and genotoxicity and demands for involvement of these compounds in biotesting.

1. Background

Plants produce a large variety of chemical compounds, which may be categorized as primary and secondary metabolites. Primary metabolites are necessary for growth and maintenance of cellular functions of the plant while secondary metabolites play an important role for example as defence (against herbivores, microbes, viruses or competing plants) and signal compounds to attract pollinating or seed dispersing animals (Lavecchia et al. 2013; Pichersky and Gang 2000; Wink 2003). Many PSMs can be seen as nature’s own pesticides and have the potential to contribute to adverse effects of chemical mixtures in aquatic ecosystems together with anthropogenic chemicals (Singh et al. 2003). Often, the production of PSMs is specific for taxonomic groups, species, genera or families. The amounts produced are typically lower than those of primary metabolites (Pichersky and Gang 2000). Secondary metabolites are diverse in their chemical nature. Most of them belong to four major classes of compounds, namely terpenes, phenols, glucosides and nitrogen-containing alkaloids (Edreva et al. 2007; Isah 2019). Due to their biological activity, PSMs have been used for drug development including anti-inflammatory, antioxidant and antiviral agents. However, many medicinally active PSMs also show toxic side effects (Lavecchia et al. 2013; Edreva et al. 2007; Puri et al. 1998; Furst 2016). For instance, bergapten, a furanocoumarin, has shown antitumor effects in a variety of cell types, but is also a potential photosensitizer that can cause severe phytophotodermatitis after either skin contact or ingestion followed by sun UV exposure (Dugrand-Judek et al. 2015; Schlatter et al. 1991; Walter et al. 1982). Some PSMs detected in the water cycle have been shown to cause severe impacts on human health such as aristolochic acids from Aristolochia clematitis causing Balkan endemic nephropathy (Tung et al. 2020) and the potent carcinogen ptaquiloside. The latter is produced by bracken fern and emitted into the water cycle particularly during rain events at toxicologically relevant concentrations (Clauson-Kaas et al. 2016; Ojika et al. 1987). Recently, natural carboline alkaloids have been demonstrated to exhibit synergistic mutagenic effects with anthropogenic aromatic amines (Muz et al. 2017).

Plants produce toxic PSMs particularly under environmental stress and release these compounds to the environment through various means such as root exudates, volatilization and animal contact as part of their defence mechanism (Aulakh et al. 2001; Al-Shatti et al. 2014; Lavecchia et al. 2013; Edreva et al. 2007). Previous research demonstrated their bioactivity and toxicity from plant extracts (Altemimi et al. 2017; Sasidharan et al. 2011) contributing to mixture toxic risk in river water (Bucheli 2014; O'Driscoll et al. 2016). In silico predictions suggest that many PSMs are persistent and mobile in the environment (Gunthardt et al. 2018). The authors identified priority phytotoxins characterized by a half-life longer than 20 days, a log DOC (organic carbon-water partition coefficient) below 4.5, rodent or aquatic toxicity and high abundance of the producing plant in Switzerland (Gunthardt et al. 2018). Assuming similar vegetation in Germany, these priority phytotoxins were used as a basis for target selection in the present study. The PSMs may be transported to river water through leaching, rain sewers and surface run-off potentially posing a risk to aquatic organisms but also to human health in case of exposure, if the water is used for human consumption. Recently, target and suspect screening of PSMs identified 12 compounds in Swiss small creeks from three compound classes including formononetin, an estrogenic isoflavone, the indole alkaloid gramine and several pyrrolizidine alkaloids (Gunthardt et al. 2020). Formononetin in concert with other isoflavones has been detected in Swiss and USA surface waters already earlier (Hoerger et al. 2009a; Kolpin et al. 2010). Thus, PSMs may add to the complex mixtures of anthropogenic organic micropollutants in water resources. Therefore, there are indications that PSMs may jeopardize water quality and affect aquatic ecosystems and human health in concert with anthropogenic compounds.

The objective of the present study was to perform a first river water target screening of PSM selected from a larger database for their expected mobility, persistence and toxicity, their expected abundance and their commercial availability as standards in two selected catchments with primarily natural vegetation and agricultural land use, respectively. We focused particularly on sampling during or after rain events in the vegetation season to enhance the probability of detection of PSMs leaching to the river water. Since hardly any quantitative toxicity data for PSMs is available, we compared water concentrations with TTC suggested for drinking water contaminants for which no toxicity data exist for a preliminary estimate of risks (Mons et al. 2013).

2. Experimental Section

2.1 Study site and water sampling

The investigated catchments are located in the north-west part of the federal state of Saxony (close to the city of Leipzig) and in Saxony-Anhalt (Bode catchment), Germany. Both catchments were selected due to their land coverage with natural and/or agricultural vegetation along the river banks. The 50 km2 large floodplain forest in Leipzig reaches along the rivers Elster, Pleiße and Luppe (EPL catchment) together with some smaller tributaries. The natural old-growth deciduous forest was historically used for the extraction of coppice and clay. It is mainly characterized by Quercus robur, Fraxinus excelsior, Acer pseudoplatanus, Ulmus minor, Alnus glutinosa, Tilia cordata, Carpinus betlulus, Acer plantanoides and Acer campestre. In spring, the forest scrub is dominated by monocotylendonous and tuberous plants such as Allium ursanium, Anemone nemorosa, Galanthus nivalis and Leucojum vernum (Klimo and Hager 2001). Thus, during their periodic growth and decay, they might leave their secondary metabolites footprints in the environmental components (e.g., soil, river). The Bode catchment is characterized by large diversity of natural and agricultural vegetation. Land use is dominated by forest (such as broad-leaved forest, coniferous forest and mixed forest) in the mountain areas and agriculture in the lowland areas. Main crops include winter wheat, triticale, winter barley, rye, rape, sugar beet and corn.

Grab water samples were collected in the vegetation growing season of the years 2018 and 2019 in spring and summer during and after rain events when leaching of PSMs was expected. A total of 38 rain event river samples were collected from both locations – 23 samples from 6 streams situated in Leipzig and 15 samples from 15 streams in the Bode catchment (Fig. 1). In the ELP catchment, samples were taken repeatedly from the same spot at different rain event days. These samples were complemented with 18 dry weather samples (8 and 10 samples from ELP and Bode catchments, respectively) from different seasons for comparison (for more information on samples see Table S1 in SI). Samples were collected in 125 mL glass bottles which had been cleaned and oven dried prior to sampling. To minimize the interferences, all sampling bottles and laboratory vessels were washed and rinsed with ethyl acetate, acetone and methanol before use. Field, trip and laboratory blanks were also included to control interferences during the sampling campaign and transportation. Samples were chilled with ice packs during transportation to the UFZ laboratory, and then stored at -24 °C until analysis.

 

2.2 Target secondary metabolite selection

Due to limited information on PSM, in silico evaluations were performed to assess their likelihood to occur in water. Prioritization for target screening was built on previous work by Gunthardt et al., identifying plant toxins with a high probability to reach the aqueous environment due to mobility and persistence (Gunthardt et al. 2018). The selected metabolites represent structurally diverse natural compounds from plant species. Thus, taking these PSM as the basic population of candidate compounds, we produced a shortlist of 150 metabolites also considering commercial availability and the probability of occurrence due to the abundance of the plants identified as the origin of these metabolites (for more information see Table S2 in SI). Furthermore, only metabolites containing the elements nitrogen, oxygen and sulfur, in addition to carbon and hydrogen were selected to allow for a likely ionization by an electrospray ion source (Moschet et al. 2013).

2.3 Chemical analysis

Water samples were prepared for direct injection by adding 25 µL of an internal standard mixture (40 ng/L), 25 µL of methanol (LC-MS grade) and 10 µL of ammonium formate buffer (2 M, pH = 3.5) to each 1-mL sample aliquot. For the chemical analysis, 100 µL of the sample were injected into a Thermo Ultimate 3000 LC system (consisting of a ternary pump, autosampler and column oven) coupled to a quadrupole-orbitrap instrument (Thermo QExactive Plus) equipped with a heated electrospray ionisation (ESI) source.

Liquid Chromatography. LC separation was performed on a Kinetex C18 EVO column (50 × 2.1 mm, 2.6 µm particle size) using a gradient elution with 0.1% of formic acid (eluent A) and methanol containing 0.1% of formic acid (eluent B) at a flow rate of 300 µL/min. After 1 min elution with 5% B, the fraction of B was linearly increased to 100% within 12 min and 100% B were kept for 11 min. Subsequently, the column was rinsed with a mixture of isopropanol + acetone 50:50 / eluent B / eluent A (85% / 10% / 5%) to remove hydrophobic matrix constituents from the column. Finally, the column was re-equilibrated to initial conditions for 5.7 min. The column was operated at 40 °C.

Mass spectrometry. The heated ESI source and the transfer capillary were both operated at 300 °C, with a spray voltage of 3.8 kV (pos. mode), a sheath gas flow rate of 45 a.u. and an auxiliary gas flow rate of 1 a.u. The full scan MS1 was recorded in a m/z range from 100–1500 with a nominal resolving power of 140,000 (referenced to m/z 200). For metabolite confirmation, data dependent MS/MS acquisition was performed at a resolving power of 70,000 in additional runs. The MS was calibrated externally every two days using the calibration mixtures of the vendor, the mass accuracy was always below 5 ppm for all analyses. All MS analyses were performed in ESI positive mode (ESIpos) since we expected a better ionizability of SPMs than in/with ESI negative mode.

2.4 Target screening

2.4.1 Qualitative target screening

The LC-HRMS raw data were converted to mzML format using ProteoWizard (version 2.1.0) (Holman et al. 2014). The centroid data were subjected to MZmine (version 2.38) for peak detection followed by peak alignment and identification (target compound annotation) (Katajamaa et al. 2006; Katajamaa and Oresic 2005; Pluskal et al. 2010). Settings for each step of the data processing are given in SI (Table S3). Further evaluation and visualization were performed using Excel 2013 (Microsoft office) and R (version 3.4.3).

Targets were identified by matching m/z and retention time between water samples and standard compounds with a mass and retention time tolerance of 5 ppm and +/- 0.1 min, respectively. Prior to clearing of false positives from the annotated list, the cut-off intensity was set to 104 to exclude signals due to noise and background. For blank correction, seven blanks were analysed together with the samples to remove noise and background contaminants. Duplicates resulting from multiple annotation were removed manually using peak resolution and intensity (for detailed steps on workflow see Figure S1 in SI). For the tentatively identified target compounds, an inclusion list was developed for data dependent acquisition (MS/MS). MS/MS experiments were conducted on authentic standard compounds and the samples to confirm the chemical structure. Diagnostic MS/MS fragments were matched with the MS/MS of reference standards. For the target compounds with low intensity in unresolved chromatograms, Parallel reaction monitoring analysis was conducted for better chromatographic peaks visualization. The XCalibur v4.0.27.10 (Thermo Fisher Scientific) software was used for analysis of extracted ion chromatograms (EICs) and mass spectra (MS1 and MS2).

2.4.2 Quantification of detected metabolites

TraceFinder 3.2 (ThermoFisher Scientific) was used for the quantification of the 12 confirmed target PSMs using extracted ion chromatograms of the full scan data. In TraceFinder, the use of only one identifier mass (precursor ion) bares the risk of false positive identification and quantification of contaminants. Thus, additional fragment ions were used to confirm the presence of target compounds and to eliminate errors in identification. For some metabolites, ions used for confirmation were not clearly detectable due to low intensity. In such cases, confirmation was complemented using Xcalibur. A series of calibration standards ranging from 1 to 5000 ng/L were used. All the calibration standards were treated exactly the same way as river water samples. The target compounds were quantified using the internal standards with the nearest retention time. The method detection limit (MDLs) (Table 1) for the detected PSMs were determined following US-EPA procedure (US-EPA 2011). The calculated concentrations below the MDLs were excluded.

2.4.3 Risk estimates

Due to a lack of toxicity data for our target compounds, we based a tentative risk estimate on TTC for non-genotoxic and non-endocrine disrupting compounds of 0.1 µg/L. We defined the ratio between measured concentration of the compounds i (ci) and TTC as risk quotient (RQ), and calculated mixture RQs as the sum of individual RQs (equation. 1) assuming a mixture RQ below one as safe for exposed humans and aquatic organisms.

(∑RQ = ∑(ci/TTC) Eq. 1

3. Result And Discussions

3.1 Metabolites detected in river waters

Peak picking followed by annotation (based on MS1 fullscan) resulted in 145 target peaks with m/z and retention time matching to the target metabolite with a tolerance of 5 ppm and +/- 0.1 min, respectively. Some target metabolites were annotated multiple times due to picking multiple peaks at a single precursor ion mass with given retention time tolerance. Removal of false positives and peak filtering using intensity and resolution reduced the target list to 106 peaks. Based on additional MS/MS fragment comparison with reference standards, we confirmed the presence of 12 target metabolites in the river waters. The identified compounds belong to different classes of natural compounds including coumarins, alkaloids, isoflavone and others. In general, the identified metabolites contain one or more phenolic groups representing a class of compounds found most abundantly in vegetation (Puri et al. 1998; Afendi et al. 2011). The names and chemical structures of the identified metabolites are given in Fig. 2.

3.2 Distribution of measured metabolites in river waters

In 18 out of 38 rain event samples PSMs were detected (9 samples from each location – Leipzig and Bode catchment), while in none of the dry weather control samples any of the target PSMs were found. In the Bode catchment, larger numbers of metabolites were detected in rivers impacted by agricultural than natural vegetation (Fig. 1 and Table S4 in SI).

Among the identified 12 metabolites, 11 were detected in samples from the ELP catchment and 10 were found in the Bode catchment. In both catchments, the coumarin derivatives scopoletin and fraxidin were detected with the highest frequency with 9 samples from each catchments (Figs. 3 and 4). Esculetin, another coumarin derivative was the third most frequently detected PSM in the Bode catchment with 20% while it was detected in 9% of the samples from the ELP catchment. The high detection frequency of fraxidin and esculetin is in good agreement with its formation by Fraxinus excelsior, a frequent tree in central European floodplains including the ones under investigation here. Scopoletin is produced by Scopolia species, but also the very frequently occurring stinging nettle Urtica dioica. However, all three compounds are present in a wide range of plants, which might contribute to emissions (Whang et al. 2005). The isoflavone formononetin is the third most frequently occurring PSM in the ELP catchment with 22%, while it was found in 7% of the Bode catchment samples. Other compounds were detected only in specific samples from a specific season such as lycorine and narciclasine occurring in Amaryllidacea, which show a high abundance within short growth periods such as Galanthus species in early spring. Although coniferyl aldehyde is a lignin component of many plants, it could be detected only in the EPL catchment and the coumarin psoralen only in Bode catchment. Natural compounds stemming likely from human consumption such as the piperidin alkaloid piperine as a component of pepper and nicotine from tobacco could be detected in both catchments.

3.3 Quantification of PSMs in river water

The target PSMs were detected in a concentration range of 1–3400 ng/L (Table 1). The concentrations of identified metabolites in individual samples is given in SI (Table S4). The highest concentrations were detected for lycorine and narciclasine with maximum concentrations of 2 and 3 µg/L during the times of high abundance of Galanthus sp. (snowdrop) and Leucojum vernum (spring snowflake). The concentrations of these phytotoxins strongly exceed the TTC of 0.1 µg/L for non-genotoxic and non-endocrine disrupting compounds and would be of concern, if these water resources would be used for drinking water production. Both compounds are highly bioactive and toxic causing among others nausea and emesis in human and animals (Kornienko and Evidente 2008; Kretzing et al. 2011). Acetylcholinesterase inhibition effects of lycorine and narciclasine were reported (Hostettmann et al. 2006; Cahlikova et al. 2013). So, apart from its toxicity, lycorine also has more positive effects, as many SPMs, such as antibacterial, anti-viral, anti-malarial, anti-allergy effects, inhibits protein and DNA synthesis and has cardiovascular protection and antitumor effects (Jahn et al. 2012).

Table 1

The concentration range (min – max, ng/L) of identified plant metabolites in river water.

Plant secondary metabolite

Chemical Formula

CAS No

m/z (M + H+)

Retention time (min)

MDL (ng/L)

Concentration range (min – max, ng/L)

Frequency of detection

Leipzig

Bode

Leipzig

Bode

Coumarin

C9H6O2

91-64-5

147.0441

7.3

11

12

43

1

1

Esculetin

C9H6O4

305-01-1

179.0336

4.2

50

116–1658

104–157

2

3

Fraxidin

C11H10O5

525-21-3

223.0600

7.8

4

56–1145

19–155

9

8

Scopoletin

C10H8O4

92-61-5

193.0496

7.1

2

9 – 47

7–49

9

9

Bergapten

C12H8O4

484-20-8

217.0495

10.1

4

510

541

1

1

Psoralen

C11H6O3

66-97-7

187.0388

9.1

3

ND

141–224

0

2

Lycorine

C16H17NO4

476-28-8

288.1225

1.0

3

1015–2331

11

2

1

Narciclasine

C14H13NO7

29477-83-6

308.0765

5.7

150

507–3353

ND

2

-

Nicotine

C10H14N2

54-11-5

163.1228

0.9

1.6

2–6

4 – 35

2

2

Piperine

C17H19NO3

94-62-2

286.1434

11.9

0.9

1–338

4 – 294

4

2

Formononetin

C16H12O4

485-72-3

269.0804

10.8

3

8–35

123

5

1

Coniferyl aldehyde

C10H10O3

458-36-6

179.0701

7.6

8

13–46

ND

2

-

ND – Not Detected

The second group of PSMs exceeding the TTC of 0.1 µg/L is the coumarin derivatives with fraxidin and esculetin concentrations of 19 to 1145 ng/L and 116 to 1658 ng/L, respectively, while coumarin and scopoletin remained below 50 ng/L (Table 1). In general, samples from the ELP catchment showed higher concentrations of coumarins than those from the Bode catchment. All four compounds have been isolated from Fraxinus excelsior (Qazi et al. 2018; Venugopala et al. 2013; Klimo and Hager 2001; Whang et al. 2005), a characteristic tree along the rivers in both catchments. Coumarins comprise a very large class of substances, found in several higher plants and constitute fused benzene and pyrone rings (Nakamura et al. 2013; Shinbo et al. 2006; Whang et al. 2005). Simple coumarins have been found to be biologically active with anti-stress, anti-fatigue, anti-gastric ulcer, anti-depressive, immuno-enhancing and anti-inflammatory effects (Whang et al. 2005; Witaicenis et al. 2010). Scopoletin, isolated from Scopolia carniolica (Solanaceae), was shown to inhibit acetylcholinesterase (Hostettmann et al. 2006).

Two furanocoumarins, bergapten and psoralen, have been detected in only two samples, each, from both catchments but in all cases above the TTC with 510 and 541 ng/L for bergapten and 141 and 224 ng/L for psoralen. They are synthesized by several plants, especially by those of the Apiaceace family (Schlatter et al. 1991). They are generally known for their strong photosensitizing activity when applied topically or accidentally get in contact to the skin. The exposure to furanocoumarins combined with long UV radiation causes cytotoxic reactions (e.g. erythema) and genotoxic responses by binding to nucleobases in DNA (Schlatter et al. 1991; Walter et al. 1982).

The flavonoid, formononetin was detected in five samples from the ELP catchment at a concentration range of 8–35 ng/L and in one sample from Bode catchment with 123 ng/L again exceeding the TTC. The latter was taken from an agricultural area (Kolpin et al. 2010; Hoerger et al. 2009b). The PSM occurs in many leguminous plants such as clover (Trifilium), an abundant species in fertile meadows and pastures but also beans such as green beans, lima beans and soy (Afendi et al. 2011, Nakamura et al. 2013, Puri et al. 1998). Formononetin has been shown to display estrogenic properties and induce angiogenesis activities (Li et al. 2015).

In both catchments, also the PSMs nicotine and piperine have been found in concentrations of 2 to 35 ng/L and 1–338 ng/L again with two samples exceeding the TTC. The input of both metabolites to the river water is very likely due to human activities, while no plants containing these compounds in the catchments are known. Both PSMs are widely consumed by humans and related to tobacco smoking and food flavouring, respectively. Nicotine is highly addictive and acts as a receptor agonist at most nicotinic acetylcholine receptors (nAChRs) (Godin et al. 2019). Piperine is a major component of Piper species (e.g Piper nigrum, Piper longum, Piper officinarum and Piper retrofractum), which are globally marketed as flavoring agent and cooking spice with a long history of human health benefits and a wide consumption (Shoba et al. 1998; Schnabel et al. 2020). Piperine has been found to have numerous medicinal applications such as antioxidant, antiplatelet, anti-inflammatory, antihypertensive, hepatoprotective, antithyroid, antitumor, antiasthmatic activity and has also been used as fertility enhance (Derosa et al. 2016). Apart from its numerous benefits, it may also have adverse effects including hemorrhagic necrosis and edema in gastrointestinal tract, urinary bladder and adrenal glands observed in animal tests with rats (Piyachaturawat et al. 1983). Zwart et al detected piperine in waste water treatment plant effluent and classified it as one of the most potent nonsteroidal estrogens (Zwart et al. 2018).

3.4 Co-occurrence of PSMs

Similar to anthropogenic compounds, also PSMs occur in mixtures. In all of the samples, where we detected our target PSMs, we found at least two of them, at two sites (11% of all samples), we detected even six co-occurring PSMs (Fig. 4(a)). The compounds fraxidin and scopoletin were common to all samples, with only one exception in the Bode catchment. Based on TTC of 0.1 µg/L, mixture risks exceeded a RQ of 1 at 16 out of 18 sites, at 5 sites mixture RQ was above 5. At one site each, even RQs of 10 and 50 were exceeded (Fig. 4(b)). This may indicate that toxic risks by frequently occurring PSMs may not be negligible and should be included in risk assessment of chemical mixtures in water resources.

4. Conclusion

By target screening of 150 prioritized PSMs in river water from two small catchments in Germany we were able to detect 12 compounds of different classes (e.g. coumarins, flavonoids, alkaloids and others), some of them occurring frequently. In many cases, concentrations of these compounds, which are known to exhibit substantial biological activity and toxicity, exceeded the TTC for drinking water individually and as mixture in almost all samples in few cases by more than one order of magnitude. This finding clearly indicates that PSMs and other natural compounds should be included into monitoring and risk assessment and should be considered in the context of drinking water abstraction. A potential risk may occur particularly during rain events promoting the leaching of PSMs to surface waters and massive occurrence of toxic plants in specific seasons. Large-scale seasonal target and suspect screening of PSMs together with toxicity testing of frequently occurring and high concentration compounds is required to estimate the contribution of PSMs to overall water pollution and to identify seasons and situations potentially posing a risk to drinking water production. Toxic risks to aquatic ecosystems might be relevant particularly in areas where vegetation undergoes drastic changes for example by massive occurrence of toxic invasive species or by substantial changes in land use.

Abbreviations

LCHRMS

Liquid chromatography coupled to high resolution mass spectrometry; PSMs:plant secondary metabolites; ELP:Elster, Luppe and Pleiße catchment; BD:Bode catchment; TTC:threshold for toxicological concern; RQ:risk quotient; MDLs:method detection limit.

Declarations

   

Acknowledgement

The authors would like to thank European Commission for funding the NaToxAq project (Grant agreement ID: 722493); analyses were run on the high-throughput platform CITEPro (Chemicals in the Environment Profiler) funded by the Helmholtz Association. A free academic licence of JChem for office (Excel) used for the structure based property and Marvin/JChem 17.21.0, 2019 ChemAxon (https:www.chemaxon.com). QGIS 3.4.15 was used to map the sampling sites.

Acknowledgement

The authors would like to thank European Commission for funding the NaToxAq project (Grant agreement ID: 722493); analyses were run on the high-throughput platform CITEPro (Chemicals in the Environment Profiler) funded by the Helmholtz Association. A free academic licence of JChem for office (Excel) used for the structure based property and Marvin/JChem 17.21.0, 2019 ChemAxon (https:www.chemaxon.com). QGIS 3.4.15 was used to map the sampling sites.

Author’s Contributions

MYN: Investigation, experimental analysis, target compound selection, data evaluation and visualization, writing (original draft). MK: Conceptualization, investigation, writing – review and editing; WB: Conceptualization, supervision, writing – review and editing; CDS, BFG and TDB: In silico prediction and selection of target compounds, review and editing. All authors read and approved the final manuscript. 

Funding

This study has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 722493.

Corresponding author

Correspondence to Werner Brack

Availability of data and materials

The datasets obtained and analyzed in the current study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interest

The author’s declare that they have no competing interest.

References

  1. Afendi FM, Okada T, Yamazaki M, Hirai-Morita A, Nakamura Y, Nakamura K, Ikeda S, Takahashi H, Altaf-Ul-Amin M, Darusman LK, Saito K, Kanaya S (2011) KNApSAcK Family Databases: Integrated Metabolite–Plant Species Databases for Multifaceted Plant Research. Plant Cell Physiol 53(2):e1(1–12). doi:10.1093/pcp/pcr165
  2. Al-Shatti AH, Redha A, Suleman P, Al-Hasan R (2014) The allelopathic potential of Conocarpus lancifolius (Engl.) leaves on dicot (Vigna sinensis L.), Monocot (Zea mays L.) and soil-Borne pathogenic fungi. American Journal of Plant Sciences 5(19):2889–2903. doi:10.4236/ajps.2014.519304
  3. Altemimi A, Lakhssassi N, Baharlouei A, Watson DG, Lightfoot DA (2017) Phytochemicals: Extraction, Isolation, and Identification of Bioactive Compounds from Plant Extracts. Plants (Basel) 6(4):1–23. doi:10.3390/plants6040042
  4. Aulakh MS, Wassmann R, Bueno C, Kreuzwieser J, Rennenberg H (2001) Characterization of Root Exudates at Different Growth Stages of Ten Rice (Oryza sativa L.) Cultivars. Plant Biol 3(2):139–148. doi:10.1055/s-2001-12905
  5. Bucheli TD (2014) Phytotoxins: Environmental Micropollutants of Concern? Environmental Science Technology 48(22):13027–13033. doi:10.1021/es504342w
  6. Cahlikova L, Hrabinova M, Kulhankova A, Benesova N, Chlebek J, Jun D, Novak Z, Macakova K, Kunes J, Kuca K, Opletal L (2013) Alkaloids from Chlidanthus fragrans and their acetylcholinesterase, butyrylcholinesterase and prolyl oligopeptidase activities. Nat Prod Commun 8(11):1541–1544
  7. Clauson-Kaas F, Ramwell C, Hansen HCB, Strobel BW (2016) Ptaquiloside from bracken in stream water at base flow and during storm events. Water Res 106:155–162. doi:10.1016/j.watres.2016.09.049
  8. Dugrand-Judek A, Olry A, Hehn A, Costantino G, Ollitrault P, Froelicher Y, Bourgaud F (2015) The Distribution of Coumarins and Furanocoumarins in Citrus Species Closely Matches Citrus Phylogeny and Reflects the Organization of Biosynthetic Pathways. PLOS ONE 10(11):e0142757. doi:10.1371/journal.pone.0142757
  9. Edreva A, Velikova V, Tsonev T, Dagnon S, Gürel AL, Aktaş L (2007) Stress-protective role of secondary metabolites: Diversity of functions and mechanisms. Gen Appl Plant Physiol 34(1–2):67–78
  10. Furst R (2016) Narciclasine - an Amaryllidaceae Alkaloid with Potent Antitumor and Anti-Inflammatory Properties. Planta Med 82(16):1389–1394. doi:10.1055/s-0042-115034
  11. Godin JR, Roy P, Quadri M, Bagdas D, Toma W, Narendrula-Kotha R, Kishta OA, Damaj MI, Horenstein NA, Papke RL, Simard AR (2019) A silent agonist of alpha7 nicotinic acetylcholine receptors modulates inflammation ex vivo and attenuates EAE. Brain Behav Immun. doi:10.1016/j.bbi.2019.12.014
  12. Gunthardt BF, Hollender J, Hungerbuhler K, Scheringer M, Bucheli TD (2018) Comprehensive Toxic Plants-Phytotoxins Database and Its Application in Assessing Aquatic Micropollution Potential. J Agric Food Chem 66(29):7577–7588. doi:10.1021/acs.jafc.8b01639
  13. Gunthardt BF, Schonsee CD, Hollender J, Hungerbuhler K, Scheringer M, Bucheli TD (2020) "Is there anybody else out there?" - First Insights from a Suspect Screening for Phytotoxins in Surface Water. Chimia 74(3):129–135. doi:10.2533/chimia.2020.129
  14. Hoerger CC, Wettstein FE, Hungerbuhler K, Bucheli TD (2009a) Occurrence and Origin of Estrogenic Isoflavones in Swiss River Waters. Environmental Science Technology 43(16):6151–6157. doi:10.1021/es901034u
  15. Hoerger CC, Wettstein FE, Hungerbühler K, Bucheli TD (2009b) Occurrence and Origin of Estrogenic Isoflavones in Swiss River Waters. Environmental Science Technology 43(16):6151–6157. doi:10.1021/es901034u
  16. Employing ProteoWizard to Convert Raw Mass Spectrometry Data
    10.1002/0471250953.bi1324s46
    Holman JD, Tabb DL, Mallick P (2014) Employing ProteoWizard to Convert Raw Mass Spectrometry Data. Current Protocols in Bioinformatics 46 (1):13.24.11–13.24.19. doi:10.1002/0471250953.bi1324s46
  17. Hostettmann K, Borloz A, Urbain A, Marston A (2006) Natural Product Inhibitors of Acetylcholinesterase. Curr Org Chem 10(8):825–847. doi:10.2174/138527206776894410
  18. Isah T (2019) Stress and defense responses in plant secondary metabolites production. Biol Res 52(1):39. doi:10.1186/s40659-019-0246-3
  19. Jahn S, Seiwert B, Kretzing S, Abraham G, Regenthal R, Karst U (2012) Metabolic studies of the Amaryllidaceous alkaloids galantamine and lycorine based on electrochemical simulation in addition to in vivo and in vitro models. Anal Chim Acta 756:60–72. doi:10.1016/j.aca.2012.10.042
  20. Katajamaa M, Miettinen J, Oresic M (2006) MZmine: toolbox for processing and visualization of mass spectrometry based molecular profile data. Bioinformatics 22(5):634–636. doi:10.1093/bioinformatics/btk039
  21. Katajamaa M, Oresic M (2005) Processing methods for differential analysis of LC/MS profile data. BMC Bioinformatics 6:179. doi:10.1186/1471-2105-6-179
  22. Klimo E, Hager H (2001) The Floodplain Forests in Europe: Current Situations and Perspectives. Brill
  23. Kolpin DW, Hoerger CC, Meyer MT, Wettstein FE, Hubbard LE, Bucheli TD (2010) Phytoestrogens and Mycotoxins in Iowa Streams: An Examination of Underinvestigated Compounds in Agricultural Basins. J Environ Qual 39(6):2089–2099. doi:10.2134/jeq2010.0121
  24. Kornienko A, Evidente A (2008) Chemistry, biology, and medicinal potential of narciclasine and its congeners. Chem Rev 108(6):1982–2014. doi:10.1021/cr078198u
  25. Kretzing S, Abraham G, Seiwert B, Ungemach FR, Krugel U, Teichert J, Regenthal R (2011) In vivo assessment of antiemetic drugs and mechanism of lycorine-induced nausea and emesis. Arch Toxicol 85(12):1565–1573. doi:10.1007/s00204-011-0719-9
  26. Lavecchia T, Rea G, Antonacci A, Giardi MT (2013) Healthy and adverse effects of plant-derived functional metabolites: the need of revealing their content and bioactivity in a complex food matrix. Crit Rev Food Sci Nutr 53(2):198–213. doi:10.1080/10408398.2010.520829
  27. Li S, Dang Y, Zhou X, Huang B, Huang X, Zhang Z, Kwan YW, Chan SW, Leung GPH, Lee SMY, Hoi MPM (2015) Formononetin promotes angiogenesis through the estrogen receptor alpha-enhanced ROCK pathway. Sci Rep 5:16815–16815. doi:10.1038/srep16815
  28. Mons MN, Heringa MB, van Genderen J, Puijker LM, Brand W, van Leeuwen CJ, Stoks P, van der Hoek JP, van der Kooij D (2013) Use of the Threshold of Toxicological Concern (TTC) approach for deriving target values for drinking water contaminants. Water Res 47(4):1666–1678. doi:10.1016/j.watres.2012.12.025
  29. Moschet C, Piazzoli A, Singer H, Hollender J (2013) Alleviating the Reference Standard Dilemma Using a Systematic Exact Mass Suspect Screening Approach with Liquid Chromatography-High Resolution Mass Spectrometry. Anal Chem 85(21):10312–10320. doi:10.1021/ac4021598
  30. Muz M, Krauss M, Kutsarova S, Schulze T, Brack W (2017) Mutagenicity in Surface Waters: Synergistic Effects of Carboline Alkaloids and Aromatic Amines. Environmental Science Technology 51(3):1830–1839. doi:10.1021/acs.est.6b05468
  31. Nakamura K, Shimura N, Otabe Y, Hirai-Morita A, Nakamura Y, Ono N, Ul-Amin MA, Kanaya S (2013) KNApSAcK-3D: a three-dimensional structure database of plant metabolites. Plant Cell Physiol 54(2):e4. doi:10.1093/pcp/pcs186
  32. O'Driscoll C, Ramwell C, Harhen B, Morrison L, Clauson-Kaas F, Hansen HC, Campbell G, Sheahan J, Misstear B, Xiao L (2016) Ptaquiloside in Irish Bracken Ferns and Receiving Waters, with Implications for Land Managers. Molecules 21(5):543. doi:10.3390/molecules21050543
  33. Ojika M, Wakamatsu K, Niwa H, Yamada K (1987) Ptaquiloside, a potent carcinogen isolated from bracken fern pteridiumaquilinum var. latiusculum: structure elucidation based on chemical and spectral evidence, and reactions with amino acids, nucleosides, and nucleotides. Tetrahedron 43(22):5261–5274. doi:10.1016/S0040-4020(01)87702-4
  34. Pichersky E, Gang DR (2000) Genetics and biochemistry of secondary metabolites in plants: an evolutionary perspective. Trends Plant Sci 5(10):439–445. doi:10.1016/s1360-1385(00)01741-6
  35. Piyachaturawat P, Glinsukon T, Toskulkao C (1983) Acute and subacute toxicity of piperine in mice, rats and hamsters. Toxicol Lett 16(3–4):351–359. doi:10.1016/0378-4274(83)90198-4
  36. Pluskal T, Castillo S, Villar-Briones A, Orešič M (2010) MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics 11(1):395. doi:10.1186/1471-2105-11-395
  37. Puri B, Hall A, Baxter H, Harborne JB, Moss GP (1998) Phytochemical Dictionary: A Handbook of Bioactive Compounds from Plants, Second Edition. CRC Press, London
  38. Qazi SS, Lombardo DA, Abou-Zaid MM (2018) A Metabolomic and HPLC-MS/MS Analysis of the Foliar Phenolics, Flavonoids and Coumarins of the Fraxinus Species Resistant and Susceptible to Emerald Ash Borer. Molecules 23(11):2734. doi:10.3390/molecules23112734
  39. Sasidharan S, Chen Y, Saravanan D, Sundram KM, Yoga Latha L (2011) Extraction, isolation and characterization of bioactive compounds from plants' extracts. Afr J Tradit Complement Altern Med 8(1):1–10
  40. Schlatter J, Zimmerli B, Dick R, Panizzon R, Schlatter C (1991) Dietary intake and risk assessment of phototoxic furocoumarins in humans. Food Chem Toxicol 29(8):523–530. doi:10.1016/0278-6915(91)90044-8
  41. Schnabel A, Cotinguiba F, Athmer B, Yang C, Westermann B, Schaks A, Porzel A, Brandt W, Schumacher F, Vogt T (2020) A piperic acid CoA ligase produces a putative precursor of piperine, the pungent principle from black pepper fruits. Plant J 102(3):569–581. doi:10.1111/tpj.14652
  42. Shinbo Y, Nakamura Y, Altaf-Ul-Amin M, Asahi H, Kurokawa K, Arita M, Saito K, Ohta D, Shibata D, Kanaya S (2006) KNApSAcK: A Comprehensive Species-Metabolite Relationship Database. In: Saito K, Dixon RA, Willmitzer L (eds) Plant Metabolomics. Springer Berlin Heidelberg, Berlin, pp 165–181. doi:10.1007/3-540-29782-0_13
  43. Shoba G, Joy D, Joseph T, Majeed M, Rajendran R, Srinivas PSSR (1998) Influence of Piperine on the Pharmacokinetics of Curcumin in Animals and Human Volunteers. Planta Med 64(04):353–356. doi:10.1055/s-2006-957450
  44. Singh HP, Batish DR, Kohli RK (2003) Allelopathic Interactions and Allelochemicals: New Possibilities for Sustainable Weed Management. Crit Rev Plant Sci 22(3–4):239–311. doi:10.1080/713610858
  45. Tung KK, Chan CK, Zhao Y, Chan KJ, Liu G, Pavlović NM, Chan W (2020) Occurrence and Environmental Stability of Aristolochic Acids in Groundwater Collected from Serbia: Links to Human Exposure and Balkan Endemic Nephropathy. Environ Sci Technol 54(3):1554–1561. doi:10.1021/acs.est.9b05337
  46. US-EPA (2011) 40 CFR Part 136, Appendix B to Part 136 - Definition and Procedure for the Determination of the Method Detection Limit-Revision 1.11. 40
  47. Venugopala KN, Rashmi V, Odhav B (2013) Review on natural coumarin lead compounds for their pharmacological activity. Biomed Res Int 2013:963248–963248. doi:10.1155/2013/963248
  48. Walter JF, Gange RW, Mendelson IR (1982) Psoralen-containing sunscreen induces phototoxicity and epidermal ornithine decarboxylase activity. J Am Acad Dermatol 6(6):1022–1027. doi:10.1016/s0190-9622(82)70086-6
  49. Whang WK, Park HS, Ham I, Oh M, Namkoong H, Kim HK, Hwang DW, Hur SY, Kim TE, Park YG, Kim JR, Kim JW (2005) Natural compounds, fraxin and chemicals structurally related to fraxin protect cells from oxidative stress. Exp Mol Med 37(5):436–446. doi:DOI 10.1038/emm.2005.54
  50. Wink M (2003) Evolution of secondary metabolites from an ecological and molecular phylogenetic perspective. Phytochemistry 64(1):3–19. doi:10.1016/s0031-9422(03)00300-5
  51. Witaicenis A, Seito LN, Di Stasi LC (2010) Intestinal anti-inflammatory activity of esculetin and 4-methylesculetin in the trinitrobenzenesulphonic acid model of rat colitis. Chem Biol Interact 186(2):211–218. doi:10.1016/j.cbi.2010.03.045
  52. Zwart N, Nio SL, Houtman CJ, de Boer J, Kool J, Hamers T, Lamoree MH (2018) High-Throughput Effect-Directed Analysis Using Downscaled in Vitro Reporter Gene Assays To Identify Endocrine Disruptors in Surface Water. Environ Sci Technol 52(7):4367–4377. doi:10.1021/acs.est.7b06604