Soumi Srimani1,2
Cosima Xenia Schmidt1
Maria Pilar Gómez-Serranillos2
Henrik Oster1
Pradeep K. Divakar2*- 1Institute of Neurobiology, Center of Brain, Behavior & Metabolism (CBBM), University of Lübeck, Lübeck, Germany
- 2Department of Pharmacology, Pharmacognosy and Botany, Faculty of Pharmacy, Complutense University of Madrid, Madrid, Spain
Background: Most mammalian cells harbor molecular circadian clocks that synchronize physiological functions with the 24-h day-night cycle. Disruption of circadian rhythms, through genetic or environmental changes, promotes the development of disorders like obesity, cardiovascular diseases, and cancer. At the cellular level, circadian, mitotic, and redox cycles are functionally coupled. Evernic (EA) and usnic acid (UA), two lichen secondary metabolites, show various pharmacological activities including anti-oxidative, anti-inflammatory, and neuroprotective action. All these effects have likewise been associated with a functional circadian clock.
Hypothesis/Purpose: To test, if the lichen compounds EA and UA modulate circadian clock function at the cellular level.
Methods: We used three different cell lines and two circadian luminescence reporter systems for evaluating dose- and time-dependent effects of EA/UA treatment on cellular clock regulation at high temporal resolution. Output parameters studied were circadian luminescence rhythm period, amplitude, phase, and dampening rate.
Results: Both compounds had marked effects on clock rhythm amplitudes and dampening independent of cell type, with UA generally showing a higher efficiency than EA. Only in fibroblast cells, significant effects on clock period were observed for UA treated cells showing shorter and EA treated cells showing longer period lengths. Transient treatment of mouse embryonic fibroblasts at different phases had only minor clock resetting effects for both compounds.
Conclusion: Secondary metabolites of lichen alter cellular circadian clocks through amplitude reduction and increased rhythm dampening.
Introduction
The circadian clock system helps to maintain the adaptation of physiological and psychological functions to the changing environmental conditions induced by the Earth’s rotation around its axis (Neumann et al., 2019). Studies reported that most cells and tissues of our body harbor molecular clocks which are in synchronization with day-night rhythms and coordinated by the suprachiasmatic nucleus (SCN) of the hypothalamus (Stephan and Zucker, 1972; Balsalobre et al., 1998; Yamazaki et al., 2000). At the cellular level, clock genes Clock (Circadian Locomotor Output Cycles Kaput), Bmal1 (Brain and Muscle ARNT-Like 1), Per1 (Period 1), Per2 (Period 2), Cry1 (Cryptochrome 1), Cry2 (Cryptochrome 2), and Rev-Erbα regulate these circadian rhythms (Hastings and Herzog, 2004). Interlocked transcriptional-translational feedback loops control these clock genes (Reppert and Weaver, 2002; Bell-Pedersen et al., 2005; Rosbash et al., 2007). Factors like sedentary lifestyle, prolonged stress or consumption of processed foods or high-sugar or high-fat diets alter circadian rhythms of cells and tissues and promote the development of a range of cardio-metabolic disorders and cancers (Tahara et al., 2015; Panda, 2016; Husse et al., 2017; Kiehn et al., 2017; Neumann et al., 2019). Many of these factors also result in the cellular accumulation of pro-oxidants and induce oxidative stress (OS), an imbalance between intracellular production of reactive oxygen species (ROS) and antioxidant defense mechanism (Charradi et al., 2013; Sharifi-Rad et al., 2020). OS activates a series of transcription factors including NF-κβ (nuclear factor kappa-light-chain-enhancer of activated B cells), AP-1 (Activator protein 1), p53 (Tumor Protein p53), HIF-1α (Hypoxia-inducible factor 1-alpha), PPAR-γ (Peroxisome proliferator-activated receptor gamma), β-catenin/Wnt, and Nrf2 (nuclear factor erythroid 2-related factor 2; Reuter et al., 2010). Activation of these transcription factors results in the expression of more than 500 different genes encoding for growth factors, inflammatory cytokines, chemokines, cell cycle regulators, and anti-inflammatory molecules, and hence OS results in cytotoxic effects in many mammalian organs and tissues (Behl et al., 1994; Uttara et al., 2009; Reuter et al., 2010). Increased ROS generation and decreased antioxidative enzyme activity (like CAT, SOD, GPx, and GST) have been reported in animals deficient in clock proteins, supporting the evidence that the activity and expression of these genes in different brain regions follow the diurnal rhythm and thus they are under the control of the endogenous circadian system (Pablos et al., 1998; Baydas et al., 2002; Kondratov et al., 2006; Krishnan et al., 2008; Yuan et al., 2017). Hence, along with altered circadian rhythms, prolonged OS also leads to sustained inflammation and chronic diseases including obesity, cancer, diabetes, cardiovascular, neurodegenerative, and pulmonary disorders (Lu and Zee, 2006; Perwez Hussain and Harris, 2007; Reuter et al., 2010; Tahara et al., 2015; Panda, 2016; Sharifi-Rad et al., 2020).
Lichens are symbiotic organisms between a fungal (the mycobiont) and an algal and/or cyanobacterial (the photobiont) partner (Calcott et al., 2018). In addition to the predominant myco- and photobionts, additional fungi, and non-photosynthetic bacteria are often associated with a lichen thallus (Grube et al., 2015; Spribille et al., 2016; Smith et al., 2020). Primary metabolites (such as simple sugars, ribitol, or glucose) produced by the photobiont partner of lichens are used for the mycobiont’s nutrition (Brodo et al., 2001; Calcott et al., 2018; Spribille et al., 2022). Lichens are well known to produce a large number of secondary metabolites with almost 1,000 known substances, the large majority of which are exclusively found in lichen forming fungi (Gómez-Serranillos et al., 2014; Calcott et al., 2018). These are deposited extracellularly, mainly in the medullary layer of lichen thallus or in the cortical layer. Phenolic compounds that originate from polyketide pathways, such as depsides, depsidones, and usnic acids are found almost exclusively in lichens. Atranorin and usnic acid are the most common and abundant cortical substances in lichens. They help in protecting the organism from UV radiation, environmental toxicity, pathogens, herbivores, and from other physical hazards (Ranković and Kosanić, 2015). Evernic acid (EA; molecular weight: 332.308 g), a depside, is produced in the medullary layer mainly by a common lichen species, Evernia prunastri, usually found in oak trees (Ter Heide et al., 1975; Joulain and Tabacchi, 2009). Other than Evernia prunastri, lichen genera such as Ramalina and Hypogymnia also produce EA (Olivier-Jimenez et al., 2019). Usnic acid (UA; molecular weight: 344.315 g), a dibenzofuran derivative, occurs in two enantiomers and is most commonly found in the cortex of Usnea species (Shukla et al., 2010). Other than in Usnea, it is also found in several lichen genera as Cladonia, Lecanora, Ramalina, Xanthoparmelia, Flavoparmelia and Alectoria (Cocchietto et al., 2002).
The production of primary metabolites by the photobionts is highly regulated by external cues like light, temperature, moisture, and gaseous concentrations in the atmosphere (Kallio and Heinonen, 1971; Kershaw and Smith, 1978; Okada et al., 1978; Korhonen and Kallio, 1987; Carré and Edmunds, 1993; Ott and Schieleit, 1994; Mittag, 2001; Segovia et al., 2003; Dodd et al., 2005; Eymann et al., 2017; Cano-Ramirez et al., 2018). kaiA, kaiB, and kaiC are the core clock genes in cyanobacteria as the circadian rhythms of it have been studied vigorously (Kondo and Ishiura, 2000; Dvornyk et al., 2003). Whereas the circadian clock genes have not been largely studied in fungi, this has only been well documented in a model organism like Neurospora crassa (Collett et al., 2002; Froehlich et al., 2002; Dunlap and Loros, 2005, 2006). Highly conserved circadian clock homologs are present in most plants and algae, while white collar-1 (WC-1), the circadian core clock component in fungi shows greater similarity in the sequences with the animal core clock component BMAL1/CLOCK (Tauber et al., 2004; Linde et al., 2017; Brody, 2019). More recently, the presence of core circadian clock genes frequency (frq), wc-1, white collar-2 (wc-2), and Frequency Interacting RNA Helicase (frh) have been reported in a lichenized fungus (Valim et al., 2022). They have also reported that the frq gene was activated in a light-dependent manner, similar to Neurospora crassa. Though it remains unclear whether lichens as a holobiont have any circadian rhythms, and thus the biosynthesis of the secondary metabolites is probably controlled by the circadian rhythms of the lichens.
Lichenized secondary metabolites impart a wide range of pharmacological activities like antioxidant, neuroprotective, cytotoxic, antimicrobial, anti-inflammatory, analgesic, and enzyme inhibitory action (Gómez-Serranillos et al., 2014; Fernández-Moriano et al., 2015a; Ranković and Kosanić, 2015; Yamamoto et al., 2015; Solárová et al., 2020). Anti-proliferative and anti-tumoral effect of UA (10 μg/ml and 50 μg/ml concentrations) on human lung cancer tumoral cells by stimulating APOPT1, CYCS, APAF1, CASP3, and CASP9 genes expression, has also been reported (Çoban et al., 2017). Fernández-Moriano et al. (2017) reported the neuroprotective potential of UA and EA with strong radical scavenging properties (ORAC and DPPH tests) and improved cellular redox status (by inhibiting H2O2 induced intracellular ROS overproduction, by restoring GSH/GSSG ratio and by lowering lipid peroxidation level through MDA) for the first time on 2017 in two models of central nervous system-like cells (U373-MG and SH-SY5Y cell lines). They have also reported that EA pre-treatment also improved the enzymatic and non-enzymatic antioxidant defense through the activation of the Nrf2 signaling pathway. UA exerts significant but more moderate effects (Fernández-Moriano et al., 2017). Further, they could be used in the therapy of oxidative stress–related diseases (Sahin et al., 2019).
Oxidative stress, mitochondrial impairment, neuroinflammation, and impaired protein degradation are the most important pathophysiologies behind the onset of neurodegenerative diseases like Parkinson’s disease, Alzheimer’s disease, etc. (Mandel et al., 2003). The brain has higher uptake of glucose and oxygen, which causes neurodegeneration via OS-related dysregulation of Nrf2-ARE defense system (Cui et al., 2016). Nrf2 is activated in response to OS which also regulates the gene expression of antioxidant enzymes, GSH-related genes, and antioxidant proteins via the antioxidant responsive element (ARE) in a cell-type dependent manner (Itoh et al., 1997; Ishii et al., 2000; Mann and Forman, 2015). Previous studies reported that astrocytes play a major role in GSH synthesis and thus protect neurons from oxidative stress induced degeneration (Sagara et al., 1993; Fernandez-Checa et al., 1997; Ishii and Mann, 2014). Moreover, circadian rhythms of antioxidative gene expression were altered by 6-hydroxydopamine, a neurotoxin via downregulating the expression of Bmal1, Per2, and other clock genes in both the in vitro and in vivo PD animal models (Wang et al., 2018). As disruption of the circadian rhythm may be associated with the inadequate inactivation of Nrf2 and dysregulation of astrocyte-neuron interactions, control of the circadian clock is particularly important for the normalization of brain functions and neuronal protections. Previous studies show that EA and UA are highly lipophilic in nature and, thus, cross lipid bilayers to reach the cytosol and mitochondria (Joseph et al., 2009; Shcherbakova et al., 2021). Hence, we can speculate that lichenized substances impart radical scavenging activities by altering the core clock gene expressions which also activates Nrf2 and prevents neurodegenerations. However, until now it remains unclear whether these varieties of pharmacological activities imparted by lichenized secondary metabolites take place by altering the circadian clocks of the organism. Thus, in this present study, we aimed to find out whether two lichen metabolites, EA and UA, modulate circadian clock function at the cellular level.
Materials and Methods
Reagents
DMEM (1×) Glutamax-I, DPBS, penicillin-streptomycin (P/S), fetal bovine serum (FBS), 1× 0.05% trypsin-EDTA (T/E), B-27, HEPES (1 M), and sodium bicarbonate 7.5% solution were purchased from Gibco, Thermo Fisher Scientific (Waltham, USA). EA and D-luciferin were purchased from Cayman Chemical, Ann Arbor, USA, and PanReac AppliChem, Darmstadt, Germany, respectively. DMEM low-glucose powder, D-(+)-glucose, UA, dexamethasone (Dex), and dimethyl-sulfoxide (DMSO) were purchased from Sigma Aldrich, St. Louis, USA. Acridine Orange Propidium Iodide (AO-PI) was purchased from Logos Biosystems, Gyeonggi-do, South Korea.
Cell Culture
Human bone osteosarcoma epithelial cells (U2OS), originally known as 2T cells, were derived from the bone tissue of a female teenager suffering from osteosarcoma. U2OS cells exhibit epithelial adherent morphology (Ponten and Saksela, 1967). An embryonic mouse hypothalamic cell line (N44) was isolated and immortalized from mouse embryonic (days 15–18 post coitum) hypothalamic primary cultures by retroviral transfer of SV40-T antigen (Cedarlane, Burlington, Canada). An embryonic mouse fibroblast (MEF) cell line was created and immortalized from PER2::LUCIFERASE reporter mice (Yoo et al., 2005) following standard protocols. N44-BL (Tsang et al., 2020), U2OS-BL (Maier et al., 2009), and MEF-P2L circadian reporter cells were authenticated by short tandem repeat (STR) profiling and tested negative for mycoplasma by PCR. Cells were maintained in DMEM with 4.5 g/L of D-glucose, supplemented with 10% (FBS) and 10,000 U P/S at 37°C with 5% CO2. Cells stably expressing Bmal1:luc reporter via lentiviral transduction were generated polyclonally by puromycin selection (Brown et al., 2005; Lin et al., 2019; Tsang et al., 2020).
Bioluminescence Assay
The setup of the bioluminescence experiment is illustrated in Figure 1A. In brief, 2*105 cells/ml of mouse hypothalamic (N44-BL), human osteosarcoma (U2OS-BL), and mouse embryonic fibroblast (MEF-P2L) cells stably expressing circadian clock reporter constructs (Bmal1::luciferase—BL or Per2::luciferase—P2L) were seeded in either 96-well plates (200 μl cell suspension/well) or 35-mm Petri dishes (2 ml cell suspension per dish) and grown to about 90% confluency (ca. 24 h at 37°C with 5% CO2). On the next day, cells were synchronized by adding 100 nM Dex for 2 h at 37°C with 5% CO2. To remove the phenol red from the DMEM, cells were washed with DPBS after the synchronization. After that, the medium was replaced with the same volume of recording medium (1 g/L DMEM low-glucose powder, 10 mM D-glucose, 3 mM 7.5% sodium bicarbonate, 10 mM HEPES, 1% 10,000 U penicillin/streptomycin, 2% B-27 supplement and 0.1 mM D-luciferin; Pilorz et al., 2020) in the dark. Plates were sealed with transparent film and luminescence was recorded for 3–5 days at 34°C using either Berthold TriStar LB 941 (Berthold Technologies, Wildbach, DE; N-44 and U2OS cells) or SpectraMax L 1-channel luminescence plate readers (Molecular Devices, San Jose, USA; MEF cells; Tsang et al., 2012). For the phase response experiment, cells were seeded in 35-mm Petri dishes. After 24 h incubation at 37°C with 5% CO2, cells were synchronized for 2 h with Dex followed by washing with DPBS and a change to the recording medium. Dishes were covered with glass coverslips and sealed with vacuum grease. Bioluminescence was measured at 32.5°C in LumiCycle 32 luminometer (Actimetrics, Evanston, USA; Tsang et al., 2012).
Figure 1. Characterization of circadian luminescence rhythms in N44, U2OS, and MEF cells. (A) Experimental setup; for more details, see text. (B) Representative raw and (C) representative normalized data and sine curve fittings of Bmal1:luc (BL) and Per2:luc (P2L) luminescence recordings of hypothalamus-derived mouse N44 (N44-BL), human bone osteosarcoma epithelial (U2OS-BL), and mouse embryonic fibroblast (MEF-P2L) cells after synchronization with dexamethasone (Dex). (D–F) Luminescence period lengths (D), rhythm phases (E), and dampening rates (F) of the same data set. Single measures are shown; error bars indicate means ± SEM (*p < 0.05, **p < 0.01, ***p < 0.001; 1-way ANOVA with Dunnett’s post-test; n = 6). ns, non-significant.
Cell Treatment and Viability Assay
Usnic acid (UA; 2,6-Diacetyl-7,9-dihydroxy-8,9b-dimethyldibenzofuran-1,3(2H,9bH)-dione) and Evernic acid (EA; 2-hydroxy-4-[(2-hydroxy-4-methoxy-6-methylbenzoyl)oxy]-6-methyl-benzoic acid) were dissolved in dimethyl-sulfoxide (DMSO) followed by dilution in phosphate-buffered saline (PBS) to the desired concentrations (final concentration of DMSO lower than 0.1%). Cell viability assay was conducted after 24 h pre-treatment with different doses of secondary metabolites of lichens using AO-PI staining. Cells were seeded in 35-mm Petri dishes (2 ml cell suspension per dish) and grown to about 90% confluency (ca. 24 h at 37°C with 5% CO2). On the next day, cell suspension was discarded, followed by a wash with DPBS. Cells were pre-treated with different concentrations of secondary metabolites and control (PBS/DMSO), added with 2 ml of cell culture medium, and incubated for approximately 24 h at 37°C + 5% CO2. On the following day cells were washed with 1 ml of DPBS/ 35 mm dish. After that, 500 μl of T/E per 35 mm dish were added and incubated for approximately 2–3 min at 37°C + 5% CO2. 2 ml cell culture medium was added to every 35 mm Petri dish and the cell suspension was centrifuged at 200× g for 5 min followed by the removal of the supernatant and resuspension in 1 ml of cell culture medium per sample. Eighteen microliters of this cell suspension was mixed with 2 μl of AO-PI stain. Ten microliters of this mixture was transferred to a channel in the Ultra-low Fluorescence counting slide for the estimation of viable cells (in the percentage of the total cells) were recorded in the LUNA-FL Dual Fluorescence Cell Counter by Logos Biosystems, Gyeonggi-do, South Korea. 94.35%, 96.35%, 96.2%, 96.3%, and 71.8% (n = 2) N44-BL cells were viable upon 24 h pre-treatment with control along with 10 μg/ml, 25 μg/ml, 50 μg/ml, and 100 μg/ml of UA, respectively. Whereas, 98.9%, 97.8%, and 96.5% (n = 2) U2OS-BL cells were viable after 24 h pre-treatment with control along with 10 μg/ml and 50 μg/ml of UA, respectively. Moreover, significant differences in the rhythm parameters were absent in cells treated with <10 μg/ml and concentrations between 10 μg/ml to 50 μg/ml of EA or UA against control. We then proceeded with our final data analysis and representation with a lower, i.e., 10 μg/ml and a higher, i.e., 50 μg/ml concentrations treatment of the studied secondary metabolites of lichens as these concentrations were not lethal to the cells. For the dose-response experiments, cells were kept in a recording medium containing different concentrations (10 μg/ml and 50 μg/ml) of UE or EA. Bioluminescence was recorded for 3–5 days. Control cells were exposed to recording medium containing the same volume of PBS/DMSO. Same doses of UA, EA or Control were added to the cells after 5 days of recording for resynchronization and bioluminescence was recorded for another 2–3 days. For phase response experiments, bioluminescence recordings were started after the addition of recording medium to the cells. At the indicated time points, equal volumes of PBS/DMSO, UA, or EA were added individually to the designated dishes without wash-out.
Data Transformation and Statistical Analysis
Circadian parameters of luminescence recordings were determined by individually adjusting for long-term trends in raw luminescence readouts by 24-h moving average baseline subtraction as described (Tsang et al., 2020). A damped sine wave was fitted to baseline-adjusted data after testing for rhythmicity using the Cosinor algorithm (Nelson et al., 1979; Refinetti, 2004). Phase shifts were determined by comparing the fit peak times of substance and PBS/DMSO treated cells (Tsang et al., 2012). Amplitude effects were calculated after normalizing to starting values to account for differences in cell density or reporter signal. Treatment times for phase response profiles were determined post hoc using the intersection of the ascending cross-section of the sine fit with the x-axis as a reference and converted into degrees. Dampening rate constants of peak magnitudes were calculated from damped sine wave regressions using GraphPad Prism 9.1.2 software (GraphPad, La Jolla, USA). All statistical analyses were performed using GraphPad Prism. Two-tailed unpaired Student’s t-tests were used for pairwise comparisons. Multi-group analyses were performed by 1-way ANOVA and Dunnett’s post-hoc test. To compare changes in rhythm parameters in dose response experiments, 2-way ANOVAs with Tukey’s post-hoc test was used. A p-value of < 0.05 was used as cut-off for significance (Tsang et al., 2012, 2020).
Results
Quantitative Characterization of Circadian Luminescence Rhythms in Three Different Cell Lines
Figure 1B depicts the representative raw bioluminescence recordings of all three cell lines. Figure 1C represents the curve-fitted and detrended rhythm of the same data set. For the determination of circadian rhythm parameters (period, phase, and dampening rate) damped sine curves were fitted to the normalized data. The mean period lengths of MEF-P2L, N44-BL, and U2OS-BL cells were 22.2 ± 0.3 h, 23.9 ± 0.1 h, and 24.8 ± 0.06 h, respectively (Figure 1D). Significant differences in period length between all three cell lines were found. Phase shifts (Figure 1E) were determined by comparing the fit peak times in these cells. The mean phases of MEF-P2L, N44-BL, and U2OS-BL cells were 29.74 ± 0.21 h, 17.72 ± 0.07 h, and 18.34 ± 0.14 h, respectively. Dampening rate constants of peak magnitudes (Figure 1F) were calculated from damped sine wave regressions. Mean dampening rates of MEF-P2L, N44-BL, and U2OS-BL cells were 0.027 ± 0.003 h−1, 0.020 ± 0.001 h−1, and 0.028 ± 0.001 h−1, respectively. Significant phase differences were observed between MEF-P2L cells and the other two cell lines. Significant differences in dampening rates were observed between N44-BL cells and the other two cell lines. In summary, irrespective of the type of cell line and the luciferase reporter among untreated cells, significantly different period length and dampening rates were observed. As expected from the feedback nature of the circadian molecular clockwork (Figure 1E), phasing was significantly different between P2L and BL rhythms.
Lichen Metabolite-Induced Dose- and Cell Type-Dependent Changes in Circadian Luminescence Rhythms
We conducted dose-response bioluminescence experiments in all the three cell lines to compare the extents to which different concentrations of UA or EA affect cellular P2L and BL rhythms. Supplementary Figures S1A,C depict the representative raw bioluminescence recording in N44-BL cells treated with different doses of UA (2.5 μg/ml, 5 μg/ml, 10 μg/ml and 50 μg/ml; Supplementary Figure S1A) and EA (5 μg/ml, 10 μg/ml, 25 μg/ml, 50 μg/ml and 100 μg/ml; Supplementary Figure S1C). Curve-fitted and detrended rhythms of the same data sets are presented under Supplementary Figure S1B (UA) and Supplementary Figure S1D (EA). For resynchronization, after 5 days of recording, the same doses of UA, EA, or Control were added to the cells with the pre-treatment and bioluminescence was recorded for another 2–3 days. In our exploratory analyses, significant resynchronization effect by the lichenized substances was not seen in any of the experiment (data are not shown here). Supplementary Figures S1A,B shows the representative raw (Supplementary Figure S1A) and curve-fitted detrended (Supplementary Figure S1B) bioluminescence Bmal1:luc rhythm of the N44 cells treated with different doses of UA against PBS/DMSO control. A significant difference in the rhythm parameters was absent in cells treated with <10 μg/ml and concentrations between 10 μg/ml to 50 μg/ml of EA or UA against control (data are not shown here). Whereas, cells became arrhythmic upon treatment with 100 μg/ml of EA (Supplementary Figures S1C,D). We speculated that the dampening in the rhythm at 50 μg/ml and arrhythmicity upon 100 μg/ml of UA or EA could be due to the cell death in such higher concentrations. We then conducted a cell viability assay using AO-PI stain in the N44 and U2OS cells after 24 h pre-treatment with different concentrations of UA against PBS/ DMSO control. 94.35%, 96.35%, 96.2%, 96.3%, and 71.8% (n = 2) N44-BL cells were viable upon 24 h pre-treatment with control along with 10 μg/ml, 25 μg/ml, 50 μg/ml and 100 μg/ml of UA, respectively. Whereas, 98.9%, 97.8%, and 96.5% (n = 2) U2OS-BL cells were viable after 24 h pre-treatment with control along with 10 μg/ml and 50 μg/ml of UA, respectively. We then proceeded with our final data analysis and representation with a lower, i.e., 10 μg/ml and a higher, i.e., 50 μg/ml concentrations treatment of the studied secondary metabolites of lichens as these concentrations were not lethal to the cells. Figures 2A, 3A, 4A show representative BL and P2L rhythms of the respective N44-BL (Figure 2A), U2OS-BL (Figure 3A), and MEF-P2L cells (Figure 4A) treated with 10 and 50 μg/ml of UA against solvent control. On the other hand, Figures 2B, 3B, 4B depict representative BL and P2L curve-fitted normalized luminescence rhythms upon treatment with 10 and 50 μg/ml of EA against solvent control in N44-BL (Figure 2B), U2OS-BL (Figure 3B), and MEF-P2L cells (Figure 4B), respectively. Irrespective of reporter and cell type, 50 μg/ml of UA and EA yielded robust dampening of the luminescence rhythm. UA at 50 μg/ml dampened the BL rhythm more strongly than that of P2L. In general, BL rhythms in U2OS cells showed greater sensitivity towards any concentrations of both secondary lichen metabolites compared to solvent controls.
Figure 2. Effects of lichenized secondary metabolites on circadian luminescence rhythms of N44-BL cells. Representative curve fits of Bmal1:luc luminescence rhythms in N44 cells after treatment with different doses of usnic acid (A, UA/Blue) or evernic acid (B, EA/Red). Amplitude (C), Period (D), Phase (E), and Dampening effects (F) on bioluminescence rhythms in N44 cells after treatment with EA (10 or 50 μg/ml; red) or UA (10 or 50 μg/ml; blue) vs. solvent control (PBS/ DMSO; black). Single measures are shown; error bars indicate means ± SEM (***p < 0.001; 2-way ANOVA with Tukey’s post-test). ###p < 0.001 between UA 10 vs. EA 10 or UA 50 vs. EA 50.
Figure 3. Effects of lichenized secondary metabolites on circadian luminescence rhythms of U2OS-BL cells. Representative curve fits of Bmal1:luc luminescence rhythms in U2OS cells after treatment with different doses of usnic acid (A, UA/Blue) or evernic acid (B, EA/Red). Amplitude (C), Period (D), Phase (E), and Dampening effects (F) on bioluminescence rhythms in U2OS cells after treatment with EA (10 or 50 μg/ml; red) or UA (10 or 50 μg/ml; blue) vs. solvent control (PBS/ DMSO; black). Single measures are shown; error bars indicate means ± SEM (*p < 0.05, ***p < 0.001; 2-way ANOVA with Tukey’s post-test).
Figure 4. Effects of lichenized secondary metabolites on circadian luminescence rhythms of MEF-P2L cells. Representative curve fits of Per2:luc luminescence rhythms in MEF cells after treatment with different doses of usnic acid (A, UA/Blue) or evernic acid (B, EA/Red). Amplitude (C), Period (D), Phase (E), and Dampening effects (F) on bioluminescence rhythms in MEF cells after treatment with EA (10 or 50 μg/ml; red) or UA (10 or 50 μg/ml; blue) vs. solvent control (PBS/ DMSO; black). Single measures are shown; error bars indicate means ± SEM (**p < 0.01, ***p < 0.001; 2-way ANOVA with Tukey’s post-test). ###p < 0.001 between UA 10 vs. EA 10 or UA 50 vs. EA 50.
Comparison of Rhythm Parameter Changes After Treatment With Lichen Secondary Metabolites
Normalized amplitude, period, phase, and dampening rate in different reporter cell models are depicted in Figures 2–4 in panels C,D,E,F, respectively. 2C–F shows the rhythm parameters of N44-BL cells whereas Figures 3C–F and Figures 4C–F represent the rhythm parameters of U2OS-BL and MEF-P2L cells, respectively. Irrespective of the reporter cell line, a dampening of the amplitude was observed in all cells treated with either UA or EA (Figures 2C, 3C, 4C). In all the three cell lines, significantly lowered amplitudes were observed upon treatment with 50 μg/ml of UA or EA compared to solvent controls. U2OS-BL (Figure 3C) and MEF-P2L (Figure 4C) cells treated with 10 μg/ml of UA had significantly dampened amplitudes. UA treated cells showed lower amplitudes than solvent controls and EA treated cells. An exception was observed only in the case of N44-BL cells treated with 10 μg/ml of UA and EA (Figure 2C). At 10 μg/ml EA treated cells showed lower amplitudes than UA cells treated at the same concentration. Subtle differences in period length were observed only in MEF-P2L cells treated with 10 μg/ml (shortened) of UA and 50 μg/ml of EA (lengthened; Figure 4D). EA 10 μg/ml treated MEF-P2L cells had a significantly longer period than 10 μg/ml of UA treated cells. Irrespective of the reporter cell line, in cells treated with 50 μg/ml of UA and EA significant phase differences were observed, though directional effects on phase differed between cell types (Figures 2E; Figure 3E; 4E). 50 μg/ml of UA or EA treated MEF-P2L (4E) and N44-BL (Figure 2E) cells showed significantly advanced phase compared to solvent controls whereas U2OS-BL (Figure 3E) showed significantly delayed phases upon treatment. N44-BL cells treated with 50 μg/ml of UA showed stronger phase delays than cells treated with EA at 50 μg/ml (Figure 2E). In contrast, among cells treated with 10 μg/ml of secondary metabolites, MEF-P2L cells showed significant phase advance upon treatment with 10 μg/ml of UA against solvent control (Figure 4E). Increased dampening was observed in BL reporter cells upon treatment with 50 μg/ml of UA in N44 cells and with both UA and EA in U2OS cells (Figures 2F, 3F, 4F). N44-BL cells treated with 50 μg/ml of UA had significantly higher dampening rates than EA treated cells (Figure 2F). In summary, we have found that higher concentrations of secondary metabolites of lichens had marked effects on amplitude and dampening, while period and phase showed only modest responses.
Phase Responses of MEF Per2:luc Rhythms After Acute Treatment With Lichen Secondary Metabolites
A phase response curve (PRC) is the graphical illustration of phase shifts as a function of the circadian treatment phase induced by stimuli like light, food, temperature, or chemicals. As 50 μg/ml of UA and EA had profound dampening effects in all the cell models and significant differences in amplitude, period, and phase were observed from dose-dependent experiments on MEF-P2L cells treated with 10 μg/ml of UA or EA (Figure 4E), we wanted to investigate if these changes of circadian parameters are dependent on the phase of the treatment. We treated MEF-P2L cells with 10 μg/ml of UA or EA (or solvent) at different phases of the Per2:luc luminescence rhythm. Figure 5A shows representative data of treatments before the peak of luminescence (at ca. 70°) whereas treatment after the peak of luminescence (i.e., at ca. 160°), before the trough (at ca. 255°), and after the trough (at ca. 355°) are shown in Figures 5B–D, respectively. Panels on the left show raw data, middle panels the respective curve-fitted detrended data. Pink arrows indicate the approximate treatment times. Except the treatment at ca. 355°, PRCs reveal that UA treated cells showed advanced phase compared to EA treated cells (Figure 5E), though overall phase effects were rather small and differences in phase shifts were not significant at any time point. While PRCs on Bmal1:luc rhythm report that EA treated cells showed an advanced phase compared UA treated cells, the overall phase effects were inconclusive (Supplementary Figure S2). Supplementary Figure S2 depicts the representative PRC of N44-BL (Supplementary Figure S2A) and U2OS-BL (Supplementary Figure S2B) cells.
Figure 5. Phase response curves (PRCs) of MEF Per2:luc rhythms after timed treatment with secondary lichen metabolites. (A–D) Raw bioluminescence data (left) and normalized curve fits (right) of luminescence in MEF-P2L cells treated with 10 μg/ml of UA (blue) or EA (red) vs. solvent control (PBS/DMSO; black) at 70° (A), 160° (B), 255° (C), or 355° (D) of the pre-treatment Per2:luc rhythm. Pink arrows depict treatment time points. (E) PRC of the same data set (n = 3 per time point; data are averages ± SEM).
Discussion
To test our hypothesis, we chose three different commonly used cellular circadian models, i.e., mouse embryonic fibroblasts (MEF) and hypothalamic neurons (N44) along with human osteosarcoma (U2OS) cell lines with two different circadian luciferase reporters (Bmal1 and Per2). In this way, changes in circadian clock rhythms upon treatment consistently observed in all models would likely also be applicable to other tissues and cell types. Along this line, previous studies reported a broad range of pharmacological activities of UA and EA like neuroprotective, cytoprotective, and antioxidant actions in healthy neuronal, cardiac, gastric cell lines, and cytotoxic, antiproliferative and anticarcinogenic effects in different cancer cell lines (MCF-7, HeLa, HCT-116, FemX, U937, LS174, SH-SY5Y, HL-60, A2780, SK-BR-3, HT29, and U373 MG; Marante et al., 2003; de Paz et al., 2010; Bačkorová et al., 2011; Bessadottir et al., 2012; Rabelo et al., 2012; Ranković et al., 2012; Brisdelli et al., 2013; Kosanić et al., 2013; Geng et al., 2018). In a previous study, Fernández-Moriano et al. (2017) had reported the neuroprotective potential of UA and EA with strong radical scavenging properties (ORAC and DPPH tests) for the first time in 2017 in two models of central nervous system-like cells (U373-MG and SH-SY5Y cell lines). They first measured approximately 55%–60% cell viability against control after inducing oxidative stress by the exogenous H2O2. Twenty-four hours of pre-treatment with 5 μg/ml EA showed the most promising and effective cytoprotection against the oxidative damage in both of the cell models tested. Whereas, that of UA in astrocytes and neurons were 2.5 μg/ml and 1 μg/ml, respectively (Fernández-Moriano et al., 2017). In our cellular models, except U2OS cells other two cells were from mouse embryonic hypothalamic neurons and fibroblasts. Moreover, exogenous oxidative stress was not induced in any of these cells.
We observed that irrespective of the cell model and luciferase reporter type, 50 μg/ml of UA or EA significantly lowered amplitudes and accelerated the dampening of cellular circadian rhythms (Figures 2C, 3C, 4C). The lowered dampening could result from cell death by the action of higher concentrations of secondary metabolites of lichens. To find out whether lichenized substances cause cell death, we conducted a cell viability assay using AO-PI stain in the N44 and U2OS cells after 24 h pre-treatment with different concentrations of UA against PBS/ DMSO control. Our results showed that the 94.35%, 96.35%, 96.2%, 96.3%, and 71.8% (n = 2) N44-BL cells were viable upon 24 h pre-treatment with control along with 10 μg/ml, 25 μg/ml, 50 μg/ml, and 100 μg/ml of UA, respectively. Whereas, after 24 h pre-treatment with control along with 10 μg/ml and 50 μg/ml of UA the 98.9%, 97.8%, and 96.5% (n = 2) U2OS-BL cells were viable, respectively. We have observed that UA had stronger dampening effects on BL rhythms than on P2L (Figures 2A, 3A; 4A). Dose response experiments showed that BL rhythms in U2OS cells were more sensitive upon treatment with any concentration of both compounds (Figures 3A,B). At this point, the mechanism behind these effects remains unclear. Reduced amplitudes—but not stronger dampening of rhythms could result from lowered cellular levels of adenosine tri-phosphate (ATP) as ATP is required for the bioluminescence signal (Gould and Subramani, 1988). Usnic acid (10 μg/ml) treatment of the breast cancer cell line T47D and EA pre-treatment of primary neurons have been shown to reduce cellular levels of ATP by inducing the phosphorylation of adenosine monophosphate kinase (AMPK; Bessadottir et al., 2012) and proton leakage through the mitochondrial membrane (Lee et al., 2021).
Circadian rhythm alterations like disrupted sleep and dysregulated circadian gene expression patterns are early markers of neurodegenerative disease progression (Hatfield et al., 2004; Breen et al., 2014; Song et al., 2015; Musiek and Holtzman, 2016; Cronin et al., 2017; Musiek et al., 2018). Fernández-Moriano et al. (2017) reported the neuroprotective potential of UA and EA with strong radical scavenging properties in U373-MG astrocytoma and SH-SY5Y neuroblastoma cell lines. Dysfunctional astrocytes cause neurodegeneration by inducing OS, followed by downregulating Nrf2 gene expression (Sagara et al., 1993; Fernandez-Checa et al., 1997; Itoh et al., 1997; Ishii et al., 2000; Ishii and Mann, 2014; Mann and Forman, 2015). Thus, it plays a critical role in brain health (Yamanaka et al., 2008; Macauley et al., 2011; Lian et al., 2015). Lananna et al. (2018) reported that BMAL1 protein regulates astrocyte activation via a cell-autonomous mechanism.
Pre-treatment of primary astrocytes with EA counteracts MPP+ induced COX2 (Cyclooxygenase 2) upregulation and glial activation by blocking the NF-κβ signaling pathway (Lee et al., 2021). It has been suggested that the regulation of NF-κβ activity in macrophages is negatively associated with Bmal1 expression (Oishi et al., 2017). Besides, in U2OS cells, like CRY1, a core clock repressor, NF-κβ activation shortens the period length and represses BMAL1-mediated E-box transcription as its subunit RELA repressed the transcriptional activity of the BMAL1/CLOCK at the circadian E-box cis-element and it also altered the diurnal locomotor behavior by modifying the rhythms in the SCN (longer period; Shen et al., 2021). They have also reported that gene expression of BMAL1, IL-6, and TNF-α were upregulated by RELA. EA also imparts neuroprotective effects as its pre-treatment significantly reduces MPP+ induced neurite shortening in the primary neurons through reduction of ROS and reduced ATP production (Lee et al., 2021). In line with this, motor function recovery is significantly accelerated in murine Parkinson’s disease models when mice are treated with EA (Lee et al., 2021).
Usnic and evernic acids prevent ROS-induced oxidative damage in various cell types via activating the NRF2 signaling pathway, inhibiting caspase-3 activity, upregulating protein levels of Bcl-2, and downregulating protein level of BAX (Fernández-Moriano et al., 2015b, 2017; Krajka-Kuźniak and Baer-Dubowska, 2021). In line with their role in clock regulation, the loss of NRF2 function in mouse fibroblasts and hepatocytes alters circadian rhythms through increased Cry2 expression and repressed CLOCK/BMAL1 regulated E-box transcription in a time-dependent manner (Wible et al., 2018). Evernic acid pre-treatment restores MPP+ induced cell viability by inducing BCL-2 and suppressing BAX protein levels in primary murine neurons (Lee et al., 2021). In the mouse central clock, p53 acts as a transcription factor that blocks BMAL1/CLOCK binding to the Per2 promoter (Miki et al., 2013), and other mechanisms of p53-clock interaction have been described (Gotoh et al., 2014; Jiang et al., 2016). Many of these pathways offer plausible mechanisms mediating the phase alterations and changes in amplitude and dampening of circadian clock rhythms after EA/UA treatment in the tested cell models. Further experiments are necessary to test this hypothesis.
In conclusion, cellular circadian clock functions in human osteosarcoma and embryonic mouse hypothalamic neurons and fibroblasts are altered by treatment with UA and EA, suggesting that lichen secondary metabolites may provide interesting novel therapeutic options for the treatment of chronodisruption-associated diseases. Our study sets a baseline for further exploration of potential natural products for therapeutic applications of chronodisruption-associated diseases.
Limitations
Serum shock or medium change synchronization effect against synchronization with Dex was not studied with any of the cellular models. Though we cannot exclude the interaction of the synchronization method and EA/UA treatment because of the simultaneous application.
Bioluminescence of all the dose-response and resynchronization experiments were recorded at 34°C in the plate readers, whereas bioluminescence for all the phase response experiments was recorded at 32.5°C due to the existing setup of the Lumicycle at this temperature.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
SS: conceptualization, methodology, formal analysis, investigation, writing—original draft. CS: supervision, validation, writing—review and editing. HO: conceptualization, methodology, formal analysis, supervision, resources, software, funding acquisition, project administration, validation, visualization, writing—review and editing. MG-S and PKD: conceptualization, supervision, resources, funding acquisition, project administration, writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by an Erasmus+ fellowship (SS). HO was supported by the German Research Foundation (DFG; RTG-1957). PKD was supported by the Spanish Ministerio de Ciencia e Innovación (grant number PID2019-105312GB-I00).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Prof. Dr. Achim Kramer for donating U2OS cells, Dr. Brid Bode for the MEF cells preparation from PER2::LUCIFERASE mice, and Dr. A. M. Neumann for help with cell culture experiments.
Abbreviations
EA, evernic acid; UA, usnic acid; BL, Bmal1:luciferase; P2L, Per2:luciferase; N44, hypothalamus-derived mouse N44; U2OS, human bone osteosarcoma epithelial; MEF, mouse embryonic fibroblast; Dex, dexamethasone; PRC, phase response curve.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2022.907308/full#supplementary-material.
Supplementary Figure 1 | Lichenized secondary metabolite’s dose effects and resynchronization changes in N44-BL cells. (A,B) Representative raw (A) and representative curve fits (B) of Bmal1:luc luminescence rhythms and resynchronization effects of different concentrations of UA treatment on N44 cells. (C,D) Representative raw (C) and representative curve fits (D) of Bmal1:luc luminescence rhythms upon treatment with different concentrations of EA on N44 cells.
Supplementary Figure 2 | Phase response curves (PRCs) of N44 and U2OS Bmal1:luc rhythms after timed treatment with secondary lichen metabolites. PRC of the N44 cells (A) and the U2OS cells (B) treated with UA (blue) or EA (red) against solvent control (PBS/ DMSO; black) at different degrees of the Bmal1:luc rhythm (n = 3 per time point; data are averages ± SEM).
References
Bačkorová, M., Bačkor, M., Mikeš, J., Jendželovský, R., and Fedoročko, P. (2011). Variable responses of different human cancer cells to the lichen compounds parietin, atranorin, usnic acid and gyrophoric acid. Toxicol. In vitro 25, 37–44. doi: 10.1016/j.tiv.2010.09.004
Balsalobre, A., Damiola, F., and Schibler, U. (1998). A serum shock induces circadian gene expression in mammalian tissue culture cells. Cell 93, 929–937. doi: 10.1016/s0092-8674(00)81199-x
Baydas, G., Gursu, M. F., Yilmaz, S., Canpolat, S., Yasar, A., Cikim, G., et al. (2002). Daily rhythm of glutathione peroxidase activity, lipid peroxidation and glutathione levels in tissues of pinealectomized rats. Neurosci. Lett. 323, 195–198. doi: 10.1016/s0304-3940(02)00144-1
Behl, C., Davis, J. B., Lesley, R., and Schubert, D. (1994). Hydrogen peroxide mediates amyloid β protein toxicity. Cell 77, 817–827. doi: 10.1016/0092-8674(94)90131-7
Bell-Pedersen, D., Cassone, V. M., Earnest, D. J., Golden, S. S., Hardin, P. E., Thomas, T. L., et al. (2005). Circadian rhythms from multiple oscillators: lessons from diverse organisms. Nat. Rev. Genet. 6, 544–556. doi: 10.1038/nrg1633
Bessadottir, M., Egilsson, M., Einarsdottir, E., Magnusdottir, I. H., Ogmundsdottir, M. H., Omarsdottir, S., et al. (2012). Proton-shuttling lichen compound usnic acid affects mitochondrial and lysosomal function in cancer cells. PLoS One 7:e51296. doi: 10.1371/journal.pone.0051296
Breen, D. P., Vuono, R., Nawarathna, U., Fisher, K., Shneerson, J. M., Reddy, A. B., et al. (2014). Sleep and circadian rhythm regulation in early Parkinson disease. JAMA Neurol. 71, 589–595. doi: 10.1001/jamaneurol.2014.65
Brisdelli, F., Perilli, M., Sellitri, D., Piovano, M., Garbarino, J. A., Nicoletti, M., et al. (2013). Cytotoxic activity and antioxidant capacity of purified lichen metabolites: an in vitro study. Phytother. Res. 27, 431–437. doi: 10.1002/ptr.4739
Brodo, I. M., Sharnoff, S. D., and Sharnoff, S. (2001). Lichens of North America. New Haven, CT: Yale University Press.
Brody, S. (2019). Circadian rhythms in fungi: structure/function/evolution of some clock components. J. Biol. Rhythms 34, 364–379. doi: 10.1177/0748730419852832
Brown, S. A., Fleury-Olela, F., Nagoshi, E., Hauser, C., Juge, C., Meier, C. A., et al. (2005). The period length of fibroblast circadian gene expression varies widely among human individuals. PLoS Biol. 3:e338. doi: 10.1371/journal.pbio.0030338
Calcott, M. J., Ackerley, D. F., Knight, A., Keyzers, R. A., and Owen, J. G. (2018). Secondary metabolism in the lichen symbiosis. Chem. Soc. Rev. 47, 1730–1760. doi: 10.1039/c7cs00431a
Cano-Ramirez, D. L., de Fraine, T. S., Griffiths, O. G., and Dodd, A. N. (2018). Photosynthesis and circadian rhythms regulate the buoyancy of marimo lake balls. Curr. Biol. 28, R869–R870. doi: 10.1016/j.cub.2018.07.027
Carré, I. A., and Edmunds, L. N. (1993). Oscillator control of cell division in Euglena: cyclic AMP oscillations mediate the phasing of the cell division cycle by the circadian clock. J. Cell Sci. 104, 1163–1173. doi: 10.1242/jcs.104.4.1163
Charradi, K., Elkahoui, S., Limam, F., and Aouani, E. (2013). High-fat diet induced an oxidative stress in white adipose tissue and disturbed plasma transition metals in rat: prevention by grape seed and skin extract. J. Physiol. Sci. 63, 445–455. doi: 10.1007/s12576-013-0283-6
Çoban, Z. D., Karaer, T., Atmaca, B., Demir, H. K., and Güran, ş. (2017). Usnic acid uses mitochondrial apopitotic pathway in it’s antitumoral role. CMJ 39, 539–545. doi: 10.7197/223.v39i31705.347452
Cocchietto, M., Skert, N., Nimis, P., and Sava, G. (2002). A review on usnic acid, an interesting natural compound. Naturwissenschaften 89, 137–146. doi: 10.1007/s00114-002-0305-3
Collett, M. A., Garceau, N., Dunlap, J. C., and Loros, J. J. (2002). Light and clock expression of the Neurospora clock gene frequency is differentially driven by but dependent on white collar-2. Genetics 160, 149–158. doi: 10.1093/genetics/160.1.149
Cronin, P., McCarthy, M. J., Lim, A. S., Salmon, D. P., Galasko, D., Masliah, E., et al. (2017). Circadian alterations during early stages of Alzheimer’s disease are associated with aberrant cycles of DNA methylation in BMAL1. Alzheimers Dement. 13, 689–700. doi: 10.1016/j.jalz.2016.10.003
Cui, Z., Zhong, Z., Yang, Y., Wang, B., Sun, Y., Sun, Q., et al. (2016). Ferrous iron induces Nrf2 expression in mouse brain astrocytes to prevent neurotoxicity. J. Biochem. Mol. Toxicol. 30, 396–403. doi: 10.1002/jbt.21803
de Paz, G. A., Raggio, J., Gómez-Serranillos, M. P., Palomino, O. M., González-Burgos, E., Carretero, M. E., et al. (2010). HPLC isolation of antioxidant constituents from Xanthoparmelia spp. J. Pharm. Biomed. Anal. 53, 165–171. doi: 10.1016/j.jpba.2010.04.013
Dodd, A. N., Salathia, N., Hall, A., Kévei, E., Tóth, R., Nagy, F., et al. (2005). Plant circadian clocks increase photosynthesis, growth, survival and competitive advantage. Science 309, 630–633. doi: 10.1126/science.1115581
Dunlap, J. C., and Loros, J. J. (2005). Analysis of circadian rhythms in Neurospora: overview of assays and genetic and molecular biological manipulation. Methods Enzymol. 393, 3–22. doi: 10.1016/S0076-6879(05)93001-2
Dunlap, J. C., and Loros, J. J. (2006). How fungi keep time: circadian system in Neurospora and other fungi. Curr. Opin. Microbiol. 9, 579–587. doi: 10.1016/j.mib.2006.10.008
Dvornyk, V., Vinogradova, O., and Nevo, E. (2003). Origin and evolution of circadian clock genes in prokaryotes. Proc. Natl. Acad. Sci. U S A 100, 2495–2500. doi: 10.1073/pnas.0130099100
Eymann, C., Lassek, C., Wegner, U., Bernhardt, J., Fritsch, O. A., Fuchs, S., et al. (2017). Symbiotic interplay of fungi, algae and bacteria within the lung lichen Lobaria pulmonaria L. Hoffm. as assessed by state-of-the-art metaproteomics. J. Proteome Res. 16, 2160–2173. doi: 10.1021/acs.jproteome.6b00974
Fernandez-Checa, J. C., Kaplowitz, N., Garcia-Ruiz, C., Colell, A., Miranda, M., Mari, M., et al. (1997). GSH transport in mitochondria: defense against TNF-induced oxidative stress and alcohol-induced defect. Am. J. Physiol. - Gastrointest. Liver Physiol. 273, G7–G17. doi: 10.1152/ajpgi.1997.273.1.G7
Fernández-Moriano, C., Divakar, P. K., Crespo, A., and Gómez-Serranillos, M. P. (2015a). Neuroprotective activity and cytotoxic potential of two Parmeliaceae lichens: identification of active compounds. Phytomedicine 22, 847–855. doi: 10.1016/j.phymed.2015.06.005
Fernández-Moriano, C., González-Burgos, E., and Gómez-Serranillos, M. P. (2015b). Mitochondria-targeted protective compounds in Parkinson’s and Alzheimer’s diseases. Oxidative Med. Cell. Longev. 2015:408927. doi: 10.1155/2015/408927
Fernández-Moriano, C., Divakar, P. K., Crespo, A., and Gómez-Serranillos, M. P. (2017). Protective effects of lichen metabolites evernic and usnic acids against redox impairment-mediated cytotoxicity in central nervous system-like cells. Food Chem. Toxicol. 105, 262–277. doi: 10.1016/j.fct.2017.04.030
Froehlich, A. C., Liu, Y., Loros, J. J., and Dunlap, J. C. (2002). White Collar-1, a circadian blue light photoreceptor, binding to the frequency promoter. Science 297, 815–819.
Geng, X., Zhang, X., Zhou, B., Zhang, C., Tu, J., Chen, X., et al. (2018). Usnic acid induces cycle arrest, apoptosis and autophagy in gastric cancer cells in vitro and in vivo. Med. Sci. Monit. 24:556. doi: 10.12659/msm.908568
Gómez-Serranillos, M. P., Fernández-Moriano, C., González-Burgos, E., Divakar, P. K., and Crespo, A. (2014). Parmeliaceae family: phytochemistry, pharmacological potential and phylogenetic features. RSC Adv. 4, 59017–59047. doi: 10.1039/C4RA09104C
Gotoh, T., Vila-Caballer, M., Santos, C. S., Liu, J., Yang, J., and Finkielstein, C. V. (2014). The circadian factor period 2 modulates p53 stability and transcriptional activity in unstressed cells. J. Mol. Cell Biol. 25, 3081–3093. doi: 10.1091/mbc.E14-05-0993
Gould, S. J., and Subramani, S. (1988). Firefly luciferase as a tool in molecular and cell biology. Anal. Biochem. 175, 5–13. doi: 10.1016/0003-2697(88)90353-3
Grube, M., Cernava, T., Soh, J., Fuchs, S., Aschenbrenner, I., Lassek, C., et al. (2015). Exploring functional contexts of symbiotic sustain within lichen-associated bacteria by comparative omics. ISME J. 9, 412–424. doi: 10.1038/ismej.2014.138
Hastings, M. H., and Herzog, E. D. (2004). Clock genes, oscillators and cellular networks in the suprachiasmatic nuclei. J. Biol. Rhythms 19, 400–413. doi: 10.1177/0748730404268786
Hatfield, C. F., Herbert, J., Van Someren, E. J., Hodges, J. R., and Hastings, M. H. (2004). Disrupted daily activity/rest cycles in relation to daily cortisol rhythms of home dwelling patients with early Alzheimer’s dementia. Brain 127, 1061–1074. doi: 10.1093/brain/awh129
Husse, J., Kiehn, J. T., Barclay, J. L., Naujokat, N., Meyer-Kovac, J., Lehnert, H., et al. (2017). Tissue-specific dissociation of diurnal transcriptome rhythms during sleep restriction in mice. Sleep 40:zsx068. doi: 10.1093/sleep/zsx068
Ishii, T., and Mann, G. E. (2014). Redox status in mammalian cells and stem cells during culture in vitro: critical roles of Nrf2 and cystine transporter activity in the maintenance of redox balance. Redox Biol. 2, 786–794. doi: 10.1016/j.redox.2014.04.008
Ishii, T., Itoh, K., Takahashi, S., Sato, H., Yanagawa, T., Katoh, Y., et al. (2000). Transcription factor Nrf2 coordinately regulates a group of oxidative stress-inducible genes in macrophages. J. Biol. Chem. 275, 16023–16029. doi: 10.1074/jbc.275.21.16023
Itoh, K., Chiba, T., Takahashi, S., Ishii, T., Igarashi, K., Katoh, Y., et al. (1997). An Nrf2/small Maf heterodimer mediates the induction of phase II detoxifying enzyme genes through antioxidant response elements. Biochem. Biophys. Res. Commun. 236, 313–322. doi: 10.1006/bbrc.1997.6943
Jiang, W., Zhao, S., Jiang, X., Zhang, E., Hu, G., Hu, B., et al. (2016). The circadian clock gene Bmal1 acts as a potential anti-oncogene in pancreatic cancer by activating the p53 tumor suppressor pathway. Cancer Lett. 371, 314–325. doi: 10.1016/j.canlet.2015.12.002
Joseph, A., Lee, T., Moland, C. L., Branham, W. S., Fuscoe, J. C., Leakey, J. E., et al. (2009). Effect of (+)-usnic acid on mitochondrial functions as measured by mitochondria-specific oligonucleotide microarray in liver of B6C3F1 mice. Mitochondrion 9, 149–158. doi: 10.1016/j.mito.2009.02.002
Joulain, D., and Tabacchi, R. (2009). Lichen extracts as raw materials in perfumery. Part 1: oakmoss. Flavour Fragrance J. 24, 49–61. doi: 10.1002/ffj.1916
Kallio, P., and Heinonen, S. (1971). Influence of short-term low temperature on net photosynthesis in some subarctic lichens. Rep. Kevo. Subarctic Res. Stat. 8, 63–72.
Kershaw, K. A., and Smith, M. M. (1978). Studies on lichen-dominated systems. XXI. The control of seasonal rates of net photosynthesis by moisture, light and temperature in Stereocaulon paschale. Can. J. Bot. 56, 2825–2830.
Kiehn, J. T., Koch, C. E., Walter, M., Brod, A., and Oster, H. (2017). Circadian rhythms and clocks in adipose tissues: current insights. Chronophysiol. Ther. 2017, 7–17. doi: 10.2147/CPT.S116242
Kondo, T., and Ishiura, M. (2000). The circadian clock of cyanobacteria. BioEssays 22, 10–15. doi: 10.1002/(SICI)1521-1878(200001)22:1<10::AID-BIES4>3.0.CO;2-A
Kondratov, R. V., Kondratova, A. A., Gorbacheva, V. Y., Vykhovanets, O. V., and Antoch, M. P. (2006). Early aging and age-related pathologies in mice deficient in BMAL1, the core component of the circadian clock. Genes Dev. 20, 1868–1873. doi: 10.1101/gad.1432206
Korhonen, P., and Kallio, P. (1987). The effect of different night conditions on the CO2 fixation in a lichen Xanthoria parietina. Photosynth. Res. 12, 3–11. doi: 10.1007/BF00019146
Kosanić, M., Manojlović, N., Janković, S., Stanojković, T., and Ranković, B. (2013). Evernia prunastri and Pseudoevernia furfuraceae lichens and their major metabolites as antioxidant, antimicrobial and anticancer agents. Food Chem. Toxicol. 53, 112–118. doi: 10.1016/j.fct.2012.11.034
Krajka-Kuźniak, V., and Baer-Dubowska, W. (2021). Modulation of Nrf2 and NF-κB signaling pathways by naturally occurring compounds in relation to cancer prevention and therapy. Are combinations better than single compounds? Int. J. Mol. Sci. 22:8223. doi: 10.3390/ijms22158223
Krishnan, N., Davis, A. J., and Giebultowicz, J. M. (2008). Circadian regulation of response to oxidative stress in Drosophila melanogaster. Biochem. Biophys. Res. Commun. 374, 299–303. doi: 10.1016/j.bbrc.2008.07.011
Lananna, B. V., Nadarajah, C. J., Izumo, M., Cedeño, M. R., Xiong, D. D., Dimitry, J., et al. (2018). Cell-autonomous regulation of astrocyte activation by the circadian clock protein BMAL1. Cell Rep. 25, 1–9. doi: 10.1016/j.celrep.2018.09.015
Lee, S., Suh, Y. J., Yang, S., Hong, D. G., Ishigami, A., Kim, H., et al. (2021). Neuroprotective and anti-inflammatory effects of evernic acid in an MPTP-induced Parkinson’s disease model. Int. J. Mol. Sci. 22:2098. doi: 10.3390/ijms22042098
Lian, H., Yang, L., Cole, A., Sun, L., Chiang, A. C. A., Fowler, S. W., et al. (2015). NFκB-activated astroglial release of complement C3 compromises neuronal morphology and function associated with Alzheimer’s disease. Neuron 85, 101–115. doi: 10.1016/j.neuron.2014.11.018
Lin, H. H., Qraitem, M., Lian, Y., Taylor, S. R., and Farkas, M. E. (2019). Analyses of BMAL1 and PER2 oscillations in a model of breast cancer progression reveal changes with malignancy. Integr. Cancer Ther. 18:1534735419836494. doi: 10.1177/1534735419836494
Linde, A. M., Eklund, D. M., Kubota, A., Pederson, E. R., Holm, K., Gyllenstrand, N., et al. (2017). Early evolution of the land plant circadian clock. New Phytol. 216, 576–590. doi: 10.1111/nph.14487
Lu, B. S., and Zee, P. C. (2006). Circadian rhythm sleep disorders. Chest 130, 1915–1923. doi: 10.1378/chest.130.6.1915
Macauley, S. L., Pekny, M., and Sands, M. S. (2011). The role of attenuated astrocyte activation in infantile neuronal ceroid lipofuscinosis. J. Neurosci. 31, 15575–15585. doi: 10.1523/JNEUROSCI.3579-11.2011
Maier, B., Wendt, S., Vanselow, J. T., Wallach, T., Reischl, S., Oehmke, S., et al. (2009). A large-scale functional RNAi screen reveals a role for CK2 in the mammalian circadian clock. Genes Dev. 23, 708–718. doi: 10.1101/gad.512209
Mandel, S., Grünblatt, E., Riederer, P., Gerlach, M., Levites, Y., Youdim, M. B., et al. (2003). Neuroprotective strategies in Parkinson’s disease. CNS Drugs 17, 729–762. doi: 10.2165/00023210-200317100-00004
Mann, G. E., and Forman, H. J. (2015). Introduction to special issue on ‘Nrf2 regulated redox signaling and metabolism in physiology and medicine. Free Radic. Biol. Med. 88, 91–92. doi: 10.1016/j.freeradbiomed.2015.08.002
Marante, F. J., Castellano, A. G., Rosas, F. E., Aguiar, J. Q., and Barrera, J. B. (2003). Identification and quantitation of allelochemicals from the lichen Lethariella canariensis: phytotoxicity and antioxidative activity. J. Chem. Ecol. 29, 2049–2071. doi: 10.1023/a:1025682318001
Miki, T., Matsumoto, T., Zhao, Z., and Lee, C. C. (2013). p53 regulates Period2 expression and the circadian clock. Nat. Commun. 4, 1–11. doi: 10.1038/ncomms3444
Mittag, M. (2001). Circadian rhythms in microalgae. Int. Rev. Cytol. 206, 213–247. doi: 10.1016/s0074-7696(01)06023-5
Musiek, E. S., and Holtzman, D. M. (2016). Mechanisms linking circadian clocks, sleep and neurodegeneration. Science 354, 1004–1008. doi: 10.1126/science.aah4968
Musiek, E. S., Bhimasani, M., Zangrilli, M. A., Morris, J. C., Holtzman, D. M., Ju, Y. E. S., et al. (2018). Circadian rest-activity pattern changes in aging and preclinical Alzheimer disease. JAMA Neurol. 75, 582–590. doi: 10.1001/jamaneurol.2017.4719
Nelson, W., Tong, Y. L., Lee, J. K., and Halberg, F. (1979). Methods for cosinor-rhythmometry. Chronobiologia 6, 305–323.
Neumann, A. M., Schmidt, C. X., Brockmann, R. M., and Oster, H. (2019). Circadian regulation of endocrine systems. Auton. Neurosci. 216, 1–8. doi: 10.1016/j.autneu.2018.10.001
Oishi, Y., Hayashi, S., Isagawa, T., Oshima, M., Iwama, A., Shimba, S., et al. (2017). Bmal1 regulates inflammatory responses in macrophages by modulating enhancer RNA transcription. Sci. Rep. 7, 1–14. doi: 10.1038/s41598-017-07100-3
Okada, M., Inoue, M., and Ikeda, T. (1978). Circadian rhythm in photosynthesis of the green alga Bryopsis maxima. Plant Cell Physiol. 19, 197–202.
Olivier-Jimenez, D., Chollet-Krugler, M., Rondeau, D., Beniddir, M. A., Ferron, S., Delhaye, T., et al. (2019). A database of high-resolution MS/MS spectra for lichen metabolites. Sci. Data 6, 1–11. doi: 10.1038/s41597-019-0305-1
Ott, S., and Schieleit, P. (1994). Influence of exogenous factors on the ethylene production by lichens. I. Influence of water content and water status conditions on ethylene production. Symbiosis 16, 187–201.
Pablos, M. I., Reiter, R. J., Ortiz, G. G., Guerrero, J. M., Agapito, M. T., Chuang, J. I., et al. (1998). Rhythms of glutathione peroxidase and glutathione reductase in brain of chick and their inhibition by light. Neurochem. Int. 32, 69–75. doi: 10.1016/s0197-0186(97)00043-0
Panda, S. (2016). Circadian physiology of metabolism. Science 354, 1008–1015. doi: 10.1126/science.aah4967
Perwez Hussain, S., and Harris, C. C. (2007). Inflammation and cancer: an ancient link with novel potentials. Int. J. Cancer 121, 2373–2380. doi: 10.1002/ijc.23173
Pilorz, V., Kolms, B., and Oster, H. (2020). Rapid jetlag resetting of behavioral, physiological, and molecular rhythms in proestrous female mice. J. Biol. Rhythms 35, 612–627. doi: 10.1177/0748730420965291
Ponten, J., and Saksela, E. (1967). Two established in vitro cell lines from human mesenchymal tumours. Int. J. Cancer 2, 434–447. doi: 10.1002/ijc.2910020505
Rabelo, T. K., Zeidán-Chuliá, F., Vasques, L. M., dos Santos, J. P. A., da Rocha, R. F., de Bittencourt Pasquali, M. A., et al. (2012). Redox characterization of usnic acid and its cytotoxic effect on human neuron-like cells (SH-SY5Y). Toxicol. In vitro 26, 304–314. doi: 10.1016/j.tiv.2011.12.003
Ranković, B., and Kosanić, M. (2015). Lichen Secondary Metabolites. New York; Cham: Springer International Publishing.
Ranković, B., Kosanić, M., Stanojković, T., Vasiljević, P., and Manojlović, N. (2012). Biological activities of Toninia candida and Usnea barbata together with their norstictic acid and usnic acid constituents. Int. J. Mol. Sci. 13, 14707–14722. doi: 10.3390/ijms131114707
Refinetti, R. (2004). Non-stationary time series and the robustness of circadian rhythms. J. Theor. Biol. 227, 571–581. doi: 10.1016/j.jtbi.2003.11.032
Reppert, S. M., and Weaver, D. R. (2002). Coordination of circadian timing in mammals. Nature 418, 935–941. doi: 10.1038/nature00965
Reuter, S., Gupta, S. C., Chaturvedi, M. M., and Aggarwal, B. B. (2010). Oxidative stress, inflammation and cancer: how are they linked? Free Radic. Biol. Med. 49, 1603–1616. doi: 10.1016/j.freeradbiomed.2010.09.006
Rosbash, M., Bradley, S., Kadener, S., Li, Y., Luo, W., Menet, J. S., et al. (2007). Transcriptional feedback and definition of the circadian pacemaker in Drosophila and animals. Cold Spring Harb. Symp. Quant. Biol. 72, 75–83. doi: 10.1101/sqb.2007.72.062
Sagara, J. I., Miura, K., and Bannai, S. (1993). Maintenance of neuronal glutathione by glial cells. J. Neurochem. 61, 1672–1676. doi: 10.1111/j.1471-4159.1993.tb09802.x
Sahin, E., Dabagoglu Psav, S., Avan, I., Candan, M., Sahinturk, V., Koparal, A. T., et al. (2019). Vulpinic acid, a lichen metabolite, emerges as a potential drug candidate in the therapy of oxidative stress-related diseases, such as atherosclerosis. Hum. Exp. Toxicol. 38, 675–684. doi: 10.1177/0960327119833745
Segovia, M., Gordillo, F. J., and Figueroa, F. L. (2003). Cyclic-AMP levels in the lichen evernia prunastri are modulated by light quantity and quality. J. Photochem. Photobiol. B: Biol. 70, 145–151. doi: 10.1016/s1011-1344(03)00074-5
Sharifi-Rad, M., Anil Kumar, N. V., Zucca, P., Varoni, E. M., Dini, L., Panzarini, E., et al. (2020). Lifestyle, oxidative stress and antioxidants: back and forth in the pathophysiology of chronic diseases. Front. Physiol. 11:694. doi: 10.3389/fphys.2020.00694
Shcherbakova, A., Strömstedt, A. A., Göransson, U., Gnezdilov, O., Turanov, A., Boldbaatar, D., et al. (2021). Antimicrobial and antioxidant activity of Evernia prunastri extracts and their isolates. World J. Microbiol. Biotechnol. 37, 1–14. doi: 10.1007/s11274-021-03099-y
Shen, Y., Endale, M., Wang, W., Morris, A. R., Francey, L. J., Harold, R. L., et al. (2021). NF-κB modifies the mammalian circadian clock through interaction with the core clock protein BMAL1. PLoS Genet. 17:e1009933. doi: 10.1371/journal.pgen.1009933
Shukla, V., Joshi, G. P., and Rawat, M. S. M. (2010). Lichens as a potential natural source of bioactive compounds: a review. Phytochem. Rev. 9, 303–314. doi: 10.1007/s11101-010-9189-6
Smith, H. B., Dal Grande, F., Muggia, L., Keuler, R., Divakar, P. K., Grewe, F., et al. (2020). Metagenomic data reveal diverse fungal and algal communities associated with the lichen symbiosis. Symbiosis 82, 133–147. doi: 10.1007/s13199-020-00699-4
Solárová, Z., Liskova, A., Samec, M., Kubatka, P., Büsselberg, D., and Solár, P. (2020). Anticancer potential of lichens’ secondary metabolites. Biomolecules 10:87. doi: 10.3390/biom10010087
Song, H., Moon, M., Choe, H. K., Han, D. H., Jang, C., Kim, A., et al. (2015). Aβ-induced degradation of BMAL1 and CBP leads to circadian rhythm disruption in Alzheimer’s disease. Mol. Neurodegener. 10, 1–15. doi: 10.1186/s13024-015-0007-x
Spribille, T., Resl, P., Stanton, D. E., and Tagirdzhanova, G. (2022). Evolutionary biology of lichen symbioses. New Phytol. 234, 1566–1582. doi: 10.1111/nph.18048
Spribille, T., Tuovinen, V., Resl, P., Vanderpool, D., Wolinski, H., Aime, M. C., et al. (2016). Basidiomycete yeasts in the cortex of ascomycete macrolichens. Science 353, 488–492. doi: 10.1126/science.aaf8287
Stephan, F. K., and Zucker, I. (1972). Circadian rhythms in drinking behavior and locomotor activity of rats are eliminated by hypothalamic lesions. Proc. Natl. Acad. Sci. USA 69, 1583–1586. doi: 10.1073/pnas.69.6.1583
Tahara, Y., Shiraishi, T., Kikuchi, Y., Haraguchi, A., Kuriki, D., Sasaki, H., et al. (2015). Entrainment of the mouse circadian clock by sub-acute physical and psychological stress. Sci. Rep. 5, 1–11. doi: 10.1038/srep11417
Tauber, E., Last, K. S., Olive, P. J., and Kyriacou, C. P. (2004). Clock gene evolution and functional divergence. J. Biol. Rhythms 19, 445–458. doi: 10.1177/0748730404268775
Ter Heide, R., Provatoroff, N., Traas, P. C., De Valois, P. J., Van der Plasse, N., Wobben, H. J., et al. (1975). Qualitative analysis of the odoriferous fraction of oakmoss (Evernia prunastri). J. Agric. Food Chem. 23, 950–957.
Tsang, A. H., Koch, C. E., Kiehn, J. T., Schmidt, C. X., and Oster, H. (2020). An adipokine feedback regulating diurnal food intake rhythms in mice. eLife 9:e55388. doi: 10.7554/eLife.55388
Tsang, A. H., Sánchez-Moreno, C., Bode, B., Rossner, M. J., Garaulet, M., Oster, H., et al. (2012). Tissue-specific interaction of Per1/2 and Dec2 in the regulation of fibroblast circadian rhythms. J. Biol. Rhythms 27, 478–489. doi: 10.1177/0748730412462838
Uttara, B., Singh, A. V., Zamboni, P., and Mahajan, R. (2009). Oxidative stress and neurodegenerative diseases: a review of upstream and downstream antioxidant therapeutic options. Curr. Neuropharmacol. 7, 65–74. doi: 10.2174/157015909787602823
Valim, F. H., Dal Grande, F., Otte, J., Singh, G., Merges, D., Schmitt, I., et al. (2022). Identification and expression of functionally conserved circadian clock genes in lichen-forming fungi. Res. Square [Preprint]. doi: 10.21203/rs.3.rs-1499077/v1
Wang, Y., Lv, D., Liu, W., Li, S., Chen, J., Shen, Y., et al. (2018). Disruption of the circadian clock alters antioxidative defense via the SIRT1-BMAL1 pathway in 6-OHDA-induced models of Parkinson’s disease. Oxidative Med. Cell. Longev. 2018:4854732. doi: 10.1155/2018/4854732
Wible, R. S., Ramanathan, C., Sutter, C. H., Olesen, K. M., Kensler, T. W., Liu, A. C., et al. (2018). NRF2 regulates core and stabilizing circadian clock loops, coupling redox and timekeeping in Mus musculus. eLife 7:e31656. doi: 10.7554/eLife.31656
Yamamoto, Y., Hara, K., Kawakami, H., and Komine, M. (2015). “Lichen substances and their biological activities,” in Recent Advances in Lichenology, eds D. Upreti, P. Divakar, V. Shukla, and R. Bajpai (New Delhi: Springer), 181–199. doi: 10.1007/978-81-322-2235-4_10
Yamanaka, K., Chun, S. J., Boillee, S., Fujimori-Tonou, N., Yamashita, H., Gutmann, D. H., et al. (2008). Astrocytes as determinants of disease progression in inherited amyotrophic lateral sclerosis. Nat. Neurosci. 11, 251–253. doi: 10.1038/nn2047
Yamazaki, S., Numano, R., Abe, M., Hida, A., Takahashi, R. I., Ueda, M., et al. (2000). Resetting central and peripheral circadian oscillators in transgenic rats. Science 288, 682–685. doi: 10.1126/science.288.5466.682
Yoo, S. H., Ko, C. H., Lowrey, P. L., Buhr, E. D., Song, E. J., Chang, S., et al. (2005). A noncanonical E-box enhancer drives mouse Period2 circadian oscillations in vivo. Proc. Natl. Acad. Sci. U S A 102, 2608–2613. doi: 10.1073/pnas.0409763102
Keywords: evernic acid, usnic acid, circadian clocks, amplitude, dampening, in vitro models
Citation: Srimani S, Schmidt CX, Gómez-Serranillos MP, Oster H and Divakar PK (2022) Modulation of Cellular Circadian Rhythms by Secondary Metabolites of Lichens. Front. Cell. Neurosci. 16:907308. doi: 10.3389/fncel.2022.907308
Received: 08 April 2022; Accepted: 20 May 2022;
Published: 23 June 2022.
Edited by:
Frédéric Gachon, The University of Queensland, AustraliaReviewed by:
Stephen Beesley, Florida State University, United StatesMarian Henryk Lewandowski, Jagiellonian University, Poland
Copyright © 2022 Srimani, Schmidt, Gómez-Serranillos, Oster and Divakar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pradeep K. Divakar, pdivakar@farm.ucm.es