==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 37261 10.1038/s41598-023-37261-3 Article Dim artificial light at night alters gene expression rhythms and growth in a key seagrass species (Posidonia oceanica) Dalle Carbonare L. laura.dallecarbonare@biology.ox.ac.uk 13 Basile A. 1 Rindi L. 2 Bulleri F. 2 Hamedeh H. 1 Iacopino S. 1 Shukla V. 1 Weits D. A. 1 Lombardi L. 2 Sbrana A. 2 Benedetti-Cecchi L. 2 Giuntoli B. 12 Licausi F. 3 Maggi E. elena.maggi@unipi.it 2 1 grid.263145.7 0000 0004 1762 600X Institute of Life Sciences, Scuola Superiore Sant’Anna, Piazza Martiri Della Libertà, 56127 Pisa, Italy 2 grid.5395.a 0000 0004 1757 3729 Dipartimento di Biologia, Universita’ di Pisa, CoNISMa, Via Luca Ghini 13, 56126 Pisa, Italy 3 grid.4991.5 0000 0004 1936 8948 Department of Biology, University of Oxford, Oxford, OX1 3RB UK 30 6 2023 30 6 2023 2023 13 106208 11 2022 19 6 2023 © The Author(s) 2023 https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Artificial light at night (ALAN) is a globally spreading anthropogenic stressor, affecting more than 20% of coastal habitats. The alteration of the natural light/darkness cycle is expected to impact the physiology of organisms by acting on the complex circuits termed as circadian rhythms. Our understanding of the impact of ALAN on marine organisms is lagging behind that of terrestrial ones, and effects on marine primary producers are almost unexplored. Here, we investigated the molecular and physiological response of the Mediterranean seagrass, Posidonia oceanica (L.) Delile, as model to evaluate the effect of ALAN on seagrass populations established in shallow waters, by taking advantage of a decreasing gradient of dim nocturnal light intensity (from < 0.01 to 4 lx) along the NW Mediterranean coastline. We first monitored the fluctuations of putative circadian-clock genes over a period of 24 h along the ALAN gradient. We then investigated whether key physiological processes, known to be synchronized with day length by the circadian rhythm, were also affected by ALAN. ALAN influenced the light signalling at dusk/night in P. oceanica, including that of shorter blue wavelengths, through the ELF3–LUX1–ZTL regulatory network, and suggested that the daily perturbation of internal clock orthologs in seagrass might have caused the recruitment of PoSEND33 and PoPSBS genes to mitigate the repercussions of a nocturnal stress on photosynthesis during the day. A long-lasting impairment of gene fluctuations in sites characterised by ALAN could explain the reduced growth of the seagrass leaves when these were transferred into controlled conditions and without lighting during the night. Our results highlight the potential contribution of ALAN to the global loss of seagrass meadows, posing questions about key interactions with a variety of other human-related stressors in urban areas, in order to develop more efficient strategies to globally preserve these coastal foundation species. Subject terms Ecophysiology Light responses Plant ecology Plant evolution Plant molecular biology Plant physiology Plant stress responses Ph.D. Program in AgrobiosciencesPh.D. Program in AgrobiodiversityUniversity of PisaPRA_2017_19 Benedetti-Cecchi L. issue-copyright-statement© Springer Nature Limited 2023 ==== Body pmcIntroduction Artificial light at night (ALAN) has come to the fore as one of the most impactful anthropogenic factors that affect wildlife1,2. With a rate of 6% increase of sky brightness per year, ALAN is spreading globally, affecting over 23% of the world’s land surfaces between 75° N and 60° S3,4. In particular, European countries and United States are those most exposed to light-polluted nights, with 88% and almost 50% of the total surface, respectively4. An alteration of the night skyglow by artificial light is expected to impact on the physiology of a wide range of organisms, including animals and plants, as well as the structure of ecological communities1. Most life forms on Earth have evolved complex biochemical and molecular circuits that are entrained by light/dark cycles and allow memory and anticipative adaptation to environmental fluctuations5. Such regulatory processes display an oscillatory pattern of about 24 h, that is termed as circadian rhythms or the internal clock5. In addition to the joint role that light plays for all the organisms to entrain circadian rhythms and as a source of information from the surrounding environment, light radiation in plants is also the main driving force to fix carbon into organic compounds through photosynthesis6. While the threat of light pollution and its effect on terrestrial ecosystems are well documented, the implication that ALAN might have in marine habitats has become a relevant issue only over the last decade7. Indeed, several of the world’s largest metropolises are located along the coasts, where they contribute to illuminate both the shoreline and the adjacent shallow sea7. In this scenario, Europe holds the record of the continent with the highest percentage (54%) of light polluted coasts8. Along the coastline, ALAN is mostly identified as a localized disturbance related to urban areas and originating from fixed light sources, either municipal or private. The spatial arrangement and type of lamps, however, create a quite heterogeneous light environment at the scale of meters, with higher and lower light intensity areas of variable extent and spectrum9. As for the light spectrum, recent advances in lighting technologies have boosted the use of broader-spectrum artificial lights (e.g., Light Emitting Diodes/LEDs) in urban areas, due to their improved energy efficiency and enhanced visual performances10. Despite their benefits, broader-spectrum lights have created more complex patterns of lighting as compared to natural moonlight, especially along the water column. In fact, the scattered blue component of the spectrum penetrates deeper than the red one, mostly affecting those species that are sensitive to short wavelengths1,8. As an additional source of spatial variability, the artificial skyglow extends the ‘sphere of potential influence’ of ALAN far beyond the localization of the point sources1, by creating gradients of light intensity tens of kilometres long, while moving away from urban centres along the coastline. Recent advances in marine ALAN research focused mostly on the effect of direct sources of LED lights on the physiology and behaviour of animals, either in intertidal or shallow subtidal systems (invertebrates and vertebrates; e.g.11–14, corals; e.g.15,16). To the best of our knowledge, effects on primary producers have been investigated only on microbial species15,17–19. This is at odds with the key role played by macroalgae and plants in marine habitats, both as oxygen producers and habitat-forming species (e.g., canopy-forming macroalgae, kelps and seagrasses20,21). The need for studies on ALAN effects on marine plants is further stimulated by research on trees commonly planted along roads (Populus ssp., Platanus ssp. and Salix spp.), which revealed species-specific sensitivity to artificial night lighting22. For instance, the proximity to streetlamps can greatly impact plant fitness and phenology in case of Platanus spp., leading to the retention of leaves into winter and the anticipation of budburst in spring23. Studies on herbaceous plants reported temporally variable, species-specific and light source dependent outcomes, both on growth-related variables and phenological aspects24–26. Interestingly, experimental exposure of bean plants to a gradient of artificial night light intensity showed a positive correlation of ALAN with plant biomass, suggesting potential effects of intensities as low as those experienced far from direct light sources27. The seagrass Posidonia oceanica (L.) Delile is the most important endemic seagrass along Mediterranean coastlines, covering about 2% of the seafloor28. It is a marine monocot species, whose meadows are adapted to grow across a broad depth range, from shallow waters down to 45-m depths, which expose them to different light and temperature regimes29. P. oceanica supports marine biodiversity through the provision of habitat and food and provides a wealth of services, ranging from the exportation of organic matter, carbon sequestration, sediment trapping and reduction of coastal erosion20. Like most seagrasses worldwide, occurrence along coastal areas poses P. oceanica meadows under threats related to urban/port development and activities, resulting, in the worst-case scenario, in their regression30,31. For those reasons, we chose P. oceanica (hereafter, Posidonia) as model species to study the effects of ALAN on shallow-water seagrass populations next to urban areas, taking advantage of a decreasing gradient of nocturnal light intensity along the coastline south of the city of Leghorn (Italy). We hypothesized that continuous exposure to artificial night lighting could impact the entrainment of gene expression rhythms in Posidonia. Among several physiological processes, circadian rhythms regulate photoprotection and photosynthesis32 and induction of flowering33. We therefore investigated whether ALAN modified the oscillatory expression of putative circadian genes in Posidonia. We further extended our analysis to photosynthetic and photoprotection genes that could participate in the output of the circadian clock, and to the regulators of photoperiodic flowering in the land angiosperm Arabidopsis thaliana. To determine if ALAN effect on gene expression could affect performances of Posidonia plants, we monitored in situ the maximum photosystem II (PSII) efficiency on dark-adapted leaves along the gradient of nocturnal light intensity (as a reliable measurement of the plant sensitivity to environmental changes34). Finally, to assess whether these effects were long-lasting, we monitored plant fitness in terms of growth under controlled conditions and without lighting during the night. Material and methods Sites location and sampling This study was carried out at three locations located south of the city of Leghorn (Italy, Western Mediterranean Sea): Rotonda (43°30′44.2″N 10°19′01.6″E), Scalinata (43°30′28.1″N 10°19′08.1″E) and Maroccone (43°29′03.9″N 10°19′51.9″E). All locations were characterized by a gently sloping rocky seabed, with occasional presence of pebbles and small boulders. Within a distance of 10–20 m from the shoreline, the biological assemblage was dominated by turf-forming and encrusting algae, with patches of Posidonia; due to the shallow water depth (about 0.5–1 m), seagrass leaves could reach the surface under low tide conditions. Despite their geographical proximity, the three locations differ in (1) night sky brightness (measured through a Sky Quality Meter, SQM—Unihedron; the unit of measurement is magnitude per squared arcsecond; Maroccone: 19.945 mag/arcsec2, Scalinata: 18.4975 mag/arcsec2, Rotonda: 15.3675 mag/arcsec2), (2) intensity of artificial light at night (defined as illuminance; the unit of measurement is lux; Maroccone: 0 ± 0.1 lx [≈ 0  µmol/s/m2], Scalinata: 1 ± 0.1 lx [≈ 0.01 µmol/s/m2], Rotonda: 4.0 ± 0.3 lx [≈ 0.1 µmol/s/m2]) and (3) distance from the industrial and merchant port (Fig. 1). We measured consistent values (as indicated in (1) and (2) above) of night sky brightness and illuminance during new moon nights in mid-June in 2017 and 2019, supporting the evidence that urban nocturnal light remains constant throughout the years. Night sky brightness and illuminance levels at the three locations decreased when moving from the centre to the periphery of the city. Despite being approximately only 500 m apart from Rotonda, Scalinata was characterized by relatively low nocturnal illuminance (more similar to Maroccone than Rotonda). We also used the LANcube system to collect data on the red and blue bands of the spectrum at the two illuminated sites35. Data on LANcube's color blue and red bands ratios to the clear band (which are linked to the percentage of blue and red to the total visible light, respectively) showed no statistically significant differences between the two sites (mean ± 1 SE (standard error); Red/Clear: Scalinata 1.7824 (1.0060), Rotonda 2.5103 (0.0036). Blue/Clear: Scalinata 0.3810 (0.2152), Rotonda 1.2067 (0.0071) (tested through one-way ANOVA).Figure 1 Diffusion of light pollution at night along the Leghorn coastline. (A) Map showing the skyglow over Europe as it appeared from space; the city of Leghorn is delimited by a white circle. Photo was taken on 13th November 2018 (credit ESA/NASA; https://www.esa.int/ESA_Multimedia/Images/2018/11/Mapping_the_night). (B) Satellite image of a portion of Tuscany (Italy), the distribution of Posidonia meadows is shown in green36; the circulation of the surface water along the coast is indicated by white arrows37. The city and the port of Leghorn are labelled on map together with the three sites where Posidonia plants were collected, distances between the city and sites are also reported for comparison (modified from Google Earth, version 9.152.0.1: https://earth.google.com/web/search/43.512270,+10.317117/@43.52372365,10.31830617,6.92320943a,12856.72487566d,35y,28.45564798h,37.3416394t,0r/data=CigiJgokCTrXa2BEwUVAEVRf3M-6wEVAGZUGYPbVpiRAIbgMprlboCRA). The values of illuminance detected for each site are reported in lux in brackets on the map. Images in (A) and (B) have been modified with the help of ‘CorelDRAW Graphics Suite 2021’ software. A total of 70 Posidonia samples were collected at 4-h intervals over a period of 24 h at 19:00, 23:00, 3:00, 7:00, 11:00, 15:00 and 19:00 (GTM + 1), at 1-m depth, in mid-June 2017 (number of replicates at each time point: n = 10), characterised by an average of 15 h of light and 9 h of darkness (dawn and dusk approximately at 5:00 and 21:00 GTM + 1 according to https://www.timeanddate.com/). Each replicate consisted of one tip of 2 cm-length selected from the youngest, fully mature leaves collected from the same seagrass meadow, avoiding those covered by epiphytes, and subsequently stored in RNAlater tissue collection (Sigma-Aldrich). The study on plant complies with relevant institutional, national, and international guidelines and legislation. The collection of P. oceanica material was performed in line with the guidelines stated in the permit granted by the city of Leghorn to the University of Pisa in order to collect plant material for research purpose (Determinazione N.3030 del 29/05/2020). Phylogenetic analysis and choice of putative circadian genes To the best of our knowledge, putative circadian clock genes and the daily fluctuation of the respective mRNAs have been only recently studied in two seagrass species, Cymodocea nodosa and Zostera marina38, but never in P. oceanica. Thus, we conducted a detailed phylogenetic analysis to retrieve the closest homologs of Arabidopsis thaliana circadian genes in Posidonia and compare their relatedness among various angiosperms, both monocotyledons and eudicotyledons, which thrive in terrestrial or aquatic environments. We focused on genes whose oscillatory expression leads to the establishment of the main endogenous circadian rhythm in A. thaliana: CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), REVEILLE (RVE8), PSEUDO-RESPONSE REGULATOR 5 (the two orthologues PRR5-1 and PRR5-2) and GIGANTEA (GI), which constitute the morning loop, and LUX ARRHYTHMO (LUX1), EARLY FLOWERING 3 (ELF3) and ZEITLUPE (ZTL), which, instead, participate to the evening loop39. We did not find evidence of a TOC1 gene in Posidonia as the closest homologous of the Arabidopsis TOC1 was more similar to the PRR5 gene family. We also performed a similar phylogenetic analysis on the flowering regulator CONSTANS (CO), whose expression is tightly controlled by the endogenous clock in terrestrial angiosperms33. We selected Posidonia and Zostera marina as examples of aquatic monocots, Nelumbo nucifera, the so-called ‘sacred lotus’, for the aquatic eudicots, Arabidopsis thaliana as terrestrial eudicot and Oryza sativa for the group of terrestrial monocots. Because of their key phylogenetic position in the evolution of plants on land, we also included the lycophyte Selaginella moellendorffii and the liverwort Marchantia polymorpha, whose homologs of the core clock genes have been recently identified40. Protein sequences were retrieved applying the BLAST algorithm41 to search the Phytozome (http://www.phytozome.net), DDBJ/ENA/GenBank (accession GEMD01000000, BioProject ID: PRJNA315106) and NCBI (https://www.ncbi.nlm.nih.gov) databases (Table S1) and Arabidopsis protein sequences were used as queries. The resulting sequences were then aligned back against the Arabidopsis thaliana protein database to ensure that they represent its closest homologs (Supplementary Files S1–S8). CCA1, PRR5, GI, ZTL, ELF3, LUX, RVE protein sequences for Marchantia polymorpha, Selaginella moellendorffii, Oryza sativa and Arabidopsis thaliana were based on40 (Table S1). The phylogenetic analysis was performed using MEGAX42. The phylogenetic trees were generated firstly by aligning the retrieved sequences using the MUSCLE algorithm43 and secondly by applying the Maximum Likelihood method and a JTT matrix-based model44. The bootstrap analysis (500 repeats) returned the percentage of trees in which the associated taxa clustered together. Choice of genes involved in light-dependent processes in Posidonia Given the capacity of a continuous source of light at night to modify the oscillatory expression of circadian genes in Posidonia, we then selected five representative genes (Chlorophyll a/b binding protein 4, PoLHCA4; Chlorophyll a/b binding protein CP29.2, PoLHCB4.2; Chlorophyll a/b binding protein 6A, 151, PoCAB-6A, PoCAB-151; Photosystem I reaction centre subunit V, PSAG), from the original group already analysed in45, and which are known to be involved in different aspects of the photosynthetic process. To test the hypothesis that excessive night lighting could cause photo-oxidative damage and, consequently, activate photoprotective mechanisms in Posidonia tip leaves, we analysed the expression of two photo-protection marker genes, Zeaxanthin epoxidase (PoZEP) and Photosystem II 22 kDa protein (PoPSBS), and of one photosynthesis-related gene, Ferredoxin-1 (PoSEND33), during the usual time course. Finally, we analysed the expression of Ribulose-bisphosphate carboxylase small chain 5B, PoSSU5B (encoding for Ribulose-1,5-Bisphosphate Carboxylase/Oxygenase, RuBisCO, small subunit) and the expression of Alcohol dehydrogenase (PoADH), considered as a marker for various stresses in plants. RNA extraction and cDNA synthesis Total RNA was extracted from Posidonia leaf tips, 2 cm long, according to the procedure described in46. The addition of fifty milligrams of polyvinylpyrrolidone (PVP) to each sample during the grinding step was the only modification introduced to the original protocol. Half microgram of RNA was subjected to DNase treatment along with cDNA synthesis, both performed with Maxima First Strand cDNA Synthesis Kit (Thermo Scientific) for subsequent quantitative real-time PCR amplification (RT-qPCR). Design of primers and quantifications of gene expression via RT-qPCR RT-qPCR primers from47 were used for PoPSBS, PoZEP, PoSSU5B, PoLHCA4, PoPSAG genes; primers from48 were used for PoCAB-6A, PoCAB-151, PoLHCB4.2 and PoSEND33. RT-qPCR primers for PoPRR5-1, PoPRR5-2, PoELF3, PoLUX, PoREV1, PoCCA1, PoZTL, PoGI, PoCO and PoADH were designed on Posidonia transcriptome available at DDBJ/ENA/GenBank (accession GEMD01000000, BioProject ID: PRJNA315106) and their sequences are listed in Table S2. RT-qPCR was carried out with the CFX384 Touch sequence detection system (Bio-Rad), using a SYBR Green PCR Master Mix (Thermo Scientific). 12.5 ng cDNA template were used in each 10 µl amplification reaction. Ubiquitin (GeneBank: GO347619) and 60S ribosomal protein L23 (GeneBank: GO347779), Posidonia housekeeping genes described in49, were chosen as reference genes in our experimental set up as their raw Ct values remained constant during the time course sampling (Fig. S1). The average Ct values of Ubiquitin and L23 was used as reference value to calculate gene expression of target genes. Relative quantification of the expression of each gene was performed using the comparative threshold cycle (Ct) method50. Data were reported as log2 x−FC, where x is a value based on primers efficiency E (x = 1 + E/100) and FC stands for expression fold change of the tested gene, calculated as difference in ΔCt from the reference (average ΔCt for samples collected in Maroccone). Primers percent efficiencies (E) are listed in Table S2 and x−FC values are listed in Table S3. Two technical replicates were used for each of the ten biological samples. Maximum photosystem II efficiency Maximum efficiency of photosystem II (PSII) of plants from Maroccone, Scalinata and Rotonda was evaluated in the field by measuring chlorophyll a fluorescence with a submersible fluorimeter (Diving PAM, Walz), on fully mature leaves, avoiding those covered by epiphytes, at 12:00, in mid-June 2019. Maximum quantum efficiency of PSII was determined as Fv/Fm, where Fv indicates the variable fluorescence, calculated as the difference between the maximum fluorescence after the saturation of the light pulse (Fm) and the minimum fluorescence (F0) after 5 min of dark adaptation (Table S4). Measurements of ALAN and illuminance were repeated and confirmed to be consistent with those collected in 2017 (as reported above in the description of the sites location). Five biological replicates were used in the experimental setup. Growth analysis Posidonia plants were collected from the three sites of Rotonda, Scalinata and Maroccone, in mid-June 2019, transferred in glass tanks (aquaria) where they were acclimated for 2 weeks and sequentially grown under the same and controlled environmental conditions (22 °C/18 °C day/night with a 16 h photoperiod and 80 µmol photon m−2 s−1 light intensity, no illumination during the night cycle). The aquaria were filled with water collected in Maroccone at the same time when plants were also harvested. Growth rate was measured on youngest fully mature leaves, newly emerged in the aquaria system, using a modified Zieman's method51. Each plant was marked by punching all leaves together with a hypodermic needle, just above the ligula of the most external leaf. Growth was monitored over 14 days, by tracking the distance of the puncture mark from the ligula. Five biological replicates (shoots) were considered for Maroccone and six for Scalinata and Rotonda plants. Statistical analysis Generalized additive mixed effect models (GAMMs) were used to evaluate the effect of light intensity on time series of residuals of gene expression derived from the regression-adjusted analysis. GAMMs allow flexibly in inferring the relationship between a response variable and a number of predictors and to deal with auto-correlated observations and complex variance–covariance structures52. For each site, gene expression was modelled as a cubic regression spline in the fixed part of the model. To account for the relationship between gene-expression and time among locations (light intensity), Site and Time (hours) were also included in the random part as random intercept and slope, respectively. GAMMs were fitted using the gamm function of the package gmcv, in R 3.5.1. Model assumptions were assessed visually using plots of residuals vs. fitted values, box-plots of residuals vs. experimental conditions (Site and Time) and QQ-plots of standardized residuals vs. normal quantiles53. Because study sites may be characterized by additional confounding environmental factors aside from ALAN, we performed GAMM on residuals of gene expression. For this purpose, long-term means (2000–2014) of daily maxima of Surface Sea Temperature (SST), Chlorophyll a concentration, Photosynthetic Active Radiation (PAR), diffusion attenuation coefficient, nitrate concentrations, phosphorous concentrations, dissolved O2 and iron concentration were retrieved for each study site from Bio-ORACLE database54. We finally regressed the gene expression against the linear combination of all environmental variables and extracted the residuals that were used as response variables in the analysis with GAMM. A linear model was used to estimate differences in maximum photosynthetic efficiency and growth rate between sites characterized by different night light intensities. This analysis was performed using the lm function of the package stats, in R 3.5.1. Results Phylogeny of Posidonia putative circadian clock genes Putative circadian clock proteins as LUX, RVE and ELF3 showed two main clusters that matched the distinction between angiosperms and non-angiosperms species (Fig. 2A,D,F,G). On the other hand, CCA1, GI and ZTL showed an additional sub-cluster within the group of the angiosperms, separating monocotyledon from eudicotyledon species (Fig. 2A,B,E,H). A special case was instead represented by PRR5 proteins, where the sub-cluster containing the closest Arabidopsis homolog showed a separation of Posidonia and Z. marina PRR5s from those of rice, lotus and Arabidopsis (Fig. 2C). For both Posidonia and Z. marina, whose flowering cycles are unpredictable and influenced by various environmental conditions55,56, we observed that their best CO homologous hits clustered separately from those of rice, Arabidopsis and lotus (Fig. 2A,I). Moreover, they subclustered with CO-like proteins found in Marchantia and Selaginella which do not produce flower-like structures.Figure 2 Phylogenetic analysis of circadian clock proteins and PoCO (CONSTANS) in Posidonia. (A) Phylogeny of land plants. Simplified phylogenetic tree adapted from57, showing the evolutionary relationship among plant lineages. Species selected for further phylogenetic analysis, which are presented in this study, are in red. (B–I) Unrooted phylogenetic trees of the main proteins regulating the morning and evening phases of the circadian rhythm (PoCCA1, PoPRR5-1, PoPRR5-2, PoLUX, PoGI, PoRVE, PoELF3 and PoZTL) and CONSTANS (PoCO) considered in this study, in aquatic monocots (Posidonia oceanica and Zostera marina) and an eudicot (Nelumbo nucifera), terrestrial eudicot (Arabidopsis thaliana) and monocot (Oryza sativa), lycophytes (Selaginella moellendorffii) and liverwort (Marchantia polymorpha). Phylogenetic relatedness was inferred by using the Maximum Likelihood method and JTT matrix-based model. Results of the bootstrap analysis are shown next to each branch. Artificial light at night impacts the expression of putative circadian clock genes in Posidonia Field measurements presented in this study were performed in mid-June 2017, when in the West Mediterranean Sea days are typically characterized by an average of 15 h of light and 9 h of darkness, with dawn and dusk occurring approximately at 5:00 and 21:00 GTM + 1. Among the eight putative circadian clock genes evaluated (Fig. 3A), all but PoLUX1 showed a significant oscillatory pattern of expression across 24 h under naturally dark conditions at night (S (Time): Maroccone, Fig. 3B, Table 1). Unexpectedly, PoLUX1 showed a rhythmic expression in its transcriptional level only under dim and intense night light conditions (i.e., at Scalinata and Rotonda sites; Table 1), whereas PoGI expression pattern, although significantly described by a cubic spline in Maroccone, showed a different oscillatory behaviour compared to the other two sites (Table 1). The Posidonia homolog of Arabidopsis ZTL exhibited a peak of expression at the end of the night, which resembled the pattern of oscillation observed in Z. marina but not in C. nodosa38.Figure 3 Posidonia circadian genes are influenced by lit nights. (A) Simplified representation of the relationships between the core genes of the circadian clock in the land angiosperm Arabidopsis thaliana. Dashed lines and solid lines indicate transcriptional or post-transcriptional regulation, respectively. The scheme is original, inspired by the text in58,59 and created it with the help of ‘CorelDRAW Graphics Suite 2021’ software. (B) Relative expression levels (Log2 x−FC) of five morning-loop genes (PoCCA1, PoRVE, PoPRR5-1, PoPRR5-2 and PoGI), on the left, and three evening-loop genes (PoZTL, PoELF3 and PoLUX1), on the right, measured in Posidonia leaf tips collected from Maroccone (dark blue), Scalinata (purple) and Rotonda (orange). Samples were collected every 4 h, starting from seven o’clock in the evening for a total period of 24 h. Data are presented as mean ± SE (number of replicates: n = 10) and analysed by Generalized Additive Mixed Models (GAMM) for significance. Table 1 Results of Generalized Additive Mixed Models (GAMMs) on gene expression levels. ^0.05 < p < 0.1, *p < 0.05, **p < 0.01, ***p < 0.001. PoCCA1 PoRVE1 PoPRR5-2 PoGI PoLUX1 PoELF3 PoPRR5-1 PoZTL Parametric coefficient  Intercept 0.06 (0.12) − 0.04 (0.11) 0.06 (0.09) − 0.04 (0.15) 0.04 (0.11) 0.03 (0.08) 0.13 (0.15) 0.05 (0.10)  Lux − 0.03 (0.05) − 0.01 (0.05) − 0.02 (0.04) 0.03 (0.06) − 0.01 (0.04) − 0.01 (0.03) − 0.08 (0.06) − 0.01 (0.04) Approximate significance of smooth terms edf edf edf edf edf edf edf edf  S (Time): Maroccone 4.293*** 4.478*** 4.225*** 2.743** 1.000 3.358*** 4.244*** 2.252**  S (Time): Scalinata 4.586*** 4.706*** 4.294*** 4.229*** 3.975** 4.382* 4.500*** 1.000  S (Time): Rotonda 4.354*** 4.280*** 3.843** 4.079*** 2.682* 1.000 4.278*** 3.154  R2 adj 0.697 0.564 0.468 0.488 0.153 0.191 0.639 0.106 PoLHCA4 PoLHCB 4.2 PoCAB-6A PoCAB-151 PoPSBS PoSEND33 PoZEP PoPSAG Parametric coefficient  Intercept 0.003 (0.06) − 0.02 (0.09) 0.01 (0.08) 0.003 (0.09) − 0.03 (0.07) − 0.04 (0.08) − 0.004 (0.05) − 0.01 (0.08)  Lux 0.001 (0.03) 0.005 (0.04) − 0.01 (0.03) − 0.004 (0.04) − 0.01 (0.03) 0.02 (0.03) 0.003 (0.02) 0.01 (0.03) Approximate significance of smooth terms edf edf edf edf edf edf edf edf  S (Time): Maroccone 3.940*** 1.000^ 3.719** 2.448** 1.000^ 1.000^ 1.000* 2.622***  S (Time): Scalinata 4.352*** 3.687*** 3.783*** 2.844* 4.215*** 4.522*** 3.531*** 4.110***  S (Time): Rotonda 4.427*** 1.993^ 1.000 1.000 4.484*** 4.559*** 1.516 3.558***  R2 adj 0.538 0.237 0.22 0.133 0.526 0.522 0.154 0.443 PoCO PoSSU5B PoADH Parametric coefficient  Intercept 0.05 (0.18) 0.04 (0.08) 0.11 (0.13)  Lux − 0.02 (0.08) − 0.03 (0.03) − 0.05 (0.05) Approximate significance of smooth terms edf edf edf  S (Time): Maroccone 2.093^ 2.441** 3.342***  S (Time): Scalinata 3.026* 3.155*** 2.848***  S (Time): Rotonda 1.000 2.965** 1.000**  R2 adj 0.082 0.299 0.34 Among the putative evening-loop genes, the transcriptional levels of PoELF3 in plants exposed to dim night lighting (i.e., at Scalinata) appeared comparable in comparison to dark conditions (i.e., at Maroccone), while the rhythmic expression completely disappeared with intense ALAN (i.e., at Rotonda). A loss of rhythmic expression in presence of dim night light in Scalinata was observed in case of PoZTL (Fig. 3B). Lit night effects on gene expression of light-dependent processes in Posidonia We observed that the expression profile of light-dependent genes in a natural night setting (Maroccone) significantly fitted into a mathematical model describing an oscillatory curve only for PoLHCA4, PoCAB-6A and PoPSAG (Fig. 4A,B, Table 1). Rhythmic oscillation disappeared at the most lit site (Rotonda) in case of PoCAB-6A (Table 1). Unexpectedly, a significant oscillatory pattern in PoSEND33 and PoPSBS appeared in Scalinata and Rotonda, where the transcript levels in Maroccone were almost constant throughout the day. Throughout 24 h of monitoring, an oscillatory expression of PoCO was observed both in Maroccone and Scalinata but not at Rotonda, where there are intense night light conditions (but the R2 adjusted value for the model was very low in this case) (Fig. 4C, Table 1).Figure 4 Light-dependent genes in relationship with artificial light at night in Posidonia. Relative expression level (Log2 x−FC) of genes of the light-harvesting complex (PoLHCA4, PoLHCB4.2, PoCAB-6A, PoPSAG and PoCAB-151) (A), photosynthesis and photo-protection (PoSEND33, PoPSBS and PoZEP) (B) and flowering (PoCO) (C), measured in leaves tip, at Maroccone (dark blue), Scalinata (purple) and Rotonda (yellow). Samples were collected every 4 h, starting from seven o’clock in the evening for a total period of 24 h. Data are presented as mean ± SE (number of replicates: n = 10) and analysed by Generalized Additive Mixed Models (GAMM) for significance. We then asked whether the observed effects on light-dependent genes might have consequences on in situ photosystem II (PSII) maximum efficiency of Posidonia. Data collected in situ on dark-adapted leaves showed a significant decrease in Fv/Fm values measured on plants at Scalinata (artificial night light intensity: 1 lx) in comparison to the control dark conditions in Maroccone (naturally dark condition) (Fig. 5, Table 2). However, we could not detect differences between the darkest and most illuminated site (at night).Figure 5 Maximum photosystem II efficiency in Posidonia under different nocturnal light conditions. Box plot depicts photosystem II (PSII) efficiency on dark-adapted Posidonia leaves (Fv/Fm values), in Maroccone (dark blue), Scalinata (purple) and Rotonda (yellow), where different light intensities at night are indicated in x axis for each site. Data (number of replicates: n = 5) were analysed for significance by a linear model (lm function of the package stats in R). Table 2 Results of linear models on maximum photosystem II efficiency and growth rate. *p < 0.05, ***p < 0.001. Maximum photosynthetic efficiency Growth rate Intercept 0.82 (0.01)*** 22.06 (2.24)*** Scalinata vs. Maroccone − 0.03 (0.01)* − 3.41 (3.03) Rotonda vs. Maroccone 0.004 (0.01) − 6.17 (2.85)* R2 adj 0.308 0.13 Posidonia growth is decreased by light at night We decided to test whether ALAN impacts Posidonia physiology when growth conditions are restored to dark nights. Therefore, we measured growth rates of Posidonia plants collected from the same three sites where we collected samples for gene expression analyses and measured chlorophyll fluorescence. To this end, the plants were collected from the three sites making sure that they retained a sufficient portion of their rhizome to enable survival, and they were subsequently installed and acclimatised at the bottom of tanks (aquaria) containing sea water from the Maroccone site. We aimed to collect plants when temperature, photoperiod and moon phase were similar to those at the time of our first sampling, for gene expression analyses. This was only possible 2 years later, in summer 2019. We confirmed that values of ALAN and sky brightness during new moon nights were comparable between 2017 and 2019. Indeed, the source of nocturnal urban lighting remained constant during the 2 years, both in terms of spectrum and intensity. The plants were grown for 2 weeks to measure the elongation rate of leaves produced after the transplant. In a controlled environment, with absence of any source of ALAN, plants harvested from Rotonda grew significantly less (with about 35% reduction) than those from Maroccone (Fig. 6A,B, Table 2). In addition to this analysis, we measured the expression of a photosynthesis-related and a stress-related gene in our sample set previously analysed for circadian genes. The transcript abundance of PoSSU5B, encoding for RuBisCO small subunit, was not significantly affected at Rotonda (Fig. 6C, Table 1). Intense night lighting negatively affected the gene expression of PoADH, resulting in a loss of rhythmicity at Rotonda (Fig. 6D, Table 1).Figure 6 ALAN effect on Posidonia leaf growth. (A) Example of a young Posidonia leaf selected for measurements of elongation, over a period of 14 days. Position of first and second measurements are marked on leaf and labelled as t1 and t2, respectively. Bar = 2 cm. (B) Box plots representing the growth rate (elongation in millimetres over 14 days) in Posidonia leaves collected from Maroccone (dark blue), Scalinata (purple) and Rotonda (orange). Measurements were taken as described in (A) and data (n = 5 for Maroccone, n = 6 for Scalinata and Rotonda) were analysed for significance by linear regression. Relative expression levels (Log2 x−FC) of PoSSU5B (C) and PoADH (D) measured in Posidonia leaf tips. Samples were collected form Maroccone (dark blue), Scalinata (purple) and Rotonda (orange), every 4 h, starting from seven o’clock in the evening for a total period of 24 h. Data are presented as mean ± SE (number of replicates: n = 10) and analysed by Generalized Additive Mixed Models (GAMM) for significance. Discussion Our understanding of the impact of ALAN on marine organisms, and primary producers, in particular, is still in its infancy, which is at odds with the global spread of night light pollution along the coastline8. In this study we investigated the molecular and physiological response of a key seagrass species, Posidonia oceanica, along the NW Mediterranean coastline. ALAN influenced the expression of light-responsive and putative circadian cycling genes at dusk/night in Posidonia. Regretfully, the timeframe of 24 h in which it was possible to reliably collect samples did not allow us to validate their bone fide association with the plant circadian clock. Nonetheless, we could measure daily fluctuations in gene expression which match those observed for their closest homologs in terrestrial angiosperms58,59. The different expression patterns in sites that differ for nocturnal lighting suggested a perturbation of the plant’s clock output. ALAN conditions were associated to reduced growth of the seagrass leaves under controlled conditions. Light is one of the most important environmental cues to entrain the circadian clock, especially in a marine environment where temperature fluctuations are attenuated in comparison to terrestrial systems. Due to the entrainment of circadian rhythms with the photoperiod, we initially focussed on genes that potentially participate in the internal clock of Posidonia. Through a phylogenetic analysis, based on protein sequence similarity across public databases, we identified the closest homologs of Arabidopsis circadian proteins in Posidonia. We also evaluated the relatedness of these sequences to their closest homologs in species that represent key steps of plant evolution in aquatic and terrestrial habitats (Fig. 2). Sequence similarity among CCA1, GI and ZTL orthologs correlated with the 140–150 million years old separation of angiosperm between monocot and dicot species60,61, although this was less clear in the case of the other five genes considered. Instead, we could not find a correlation with the habitats, suggesting that their colonization occurred after the divergence of these sequences in the clades described above and that no differentiation in these coding sequences contributed to the adaptation to terrestrial or aquatic environments62. By monitoring the fluctuation of these putative circadian-clock genes in Posidonia for a period of 24 h, we observed that the genes PoELF3, PoLUX1 and PoZTL lost or altered their rhythmicity in the presence of dimly bright nights (characterized, in this study, by 1–4 lx) (Fig. 3B). In Arabidopsis, the transcription factors ELF3 and LUX1, together with ELF4, constitute a complex that acts during the evening phase of the circadian rhythm, known as evening complex63. In particular, ELF3 and ZTL participate to blue light entrainment of the circadian clock both in Arabidopsis and rice, where they control circadian outputs as well as feedback regulation of light perception64–67. The orthologous blue light receptor ZTL represents an interesting case in Posidonia, where its RNA levels fluctuate on a 24 h period, with a peak at the end of the night phase, as observed also in the seagrass Z. marina (Fig. 3B)38. This is different from what observed in Arabidopsis, where ZTL abundance is instead controlled at the post-translational level68, and suggests the existence of a feedback mechanism that sensitizes cells to blue light under dark conditions. Since blue light is expected to penetrate the water column the farthest8, we speculated that the effect of dim nocturnal light on PoELF3 and PoZTL expression supports their conserved function in blue light signalling. Conversely, ALAN did not affect the fluctuation of the putative morning-loop genes considered in this study. Prompted by the observed impact of ALAN on the daily rhythmicity of putative circadian clock genes expression, we investigated whether physiological processes known to be synchronized with day length by the circadian rhythm, were also affected under ALAN conditions. The rhythmic abundance of components of the photosynthetic apparatus has been attested from marine algae to flowering plants, and it includes proteins involved in different aspects of the process: oxygen evolution, light harvesting, electron transport and CO2 fixation, suggesting a selective pressure in order to preserve this advantageous trait during plant evolution32,69. Among the Posidonia photosynthetic genes tested in this study, only PoCAB-6A was influenced by night lightning, displaying a loss of rhythmicity at the most lit site (4 lx, Rotonda) (Fig. 4A). Previously, PoCAB-6A expression has been shown to increase at the end of the day in Posidonia plants growing in deep waters45, suggesting that the loss of rhythmicity of PoCAB-6A in Rotonda could be associated with the misperception of the night period at this site. Unexpectedly, in case of ALAN (i.e., at Scalinata and Rotonda sites), we attested the acquisition of an oscillatory pattern by PoSEND33 and PoPSBS expression (Fig. 4B). Since PSBS is known for its photoacclimation and photo-protective action on the photosystem II70, we speculate that the perturbation of the plant’s internal clock might have caused their recruitment to mitigate the repercussions of nocturnal stress on photosynthesis during the day. In land plants, synchronisation of the internal circadian rhythm with the external day–night cycle allows that both carbon fixation during the day, coupled with photosynthesis, and starch consumption during the night occur at the best conditions to optimize plant growth and development32. The asynchrony between external and internal clocks has been shown to impair growth and reduce fitness71–73. Although we did not observe an effect on maximum photosynthetic efficiency by ALAN at the most lit site, Posidonia plants exposed to nocturnal lighting of 4 lx exhibited reduced growth when transferred into controlled growth conditions (aquaria) in the absence of any source of nocturnal lighting, in our laboratory (Fig. 6B). We concluded that Posidonia plants adapted to grow under specific and local conditions that characterise their native sites. Since the three sites are sufficiently close to each other, and mainly distinct by ALAN, these results, although preliminary, suggest that high intensity ALAN can exert a long-lasting effect on plant physiology. The few studies that have investigated the influence of ALAN on terrestrial plants reported temporally variable, species-specific and light source dependent outcomes. In particular, observed effects on growth-related variables were either positive24,26 or negative25, reflecting the complexity of species-specific processes and mechanisms involved in plant growth. In the model plant Arabidopsis, sugar and starch consumption, as well as their allocation throughout the plant, are under the control of the circadian clock; in particular, sugar consumption is regulated antagonistically by the morning protein CCA1 and the evening complex ELF3–ELF4–LUX74. Due to the observed effect of ALAN on PoELF3–PoLUX pattern of expression, we hypothesise that, similarly to Arabidopsis, clock-controlled metabolic pathways might be responsible for the observed growth reduction in the seagrass Posidonia oceanica. PoADH could have a possible role in this, as we observed that intense nocturnal lighting negatively affected its expression (Fig. 6D). Alcohol dehydrogenase plays an important role in fermentative metabolism which is actuated in response to various abiotic and abiotic stresses, as low oxygen conditions, cold, drought, salt stress and pathogen infection75–79. The loss of PoADH periodicity observed at the most lit site (Rotonda), could indicate a metabolic impairment and a more susceptibility to other stressors (e.g. cold temperatures, chemical pollutants and pathogens), in response to ALAN. Flowering is one of the many physiological processes which are under control of the circadian clock33. Although the molecular regulation of flowering in Posidonia has not been characterised yet, and its occurrence in the Northern Mediterranean is far more variable and rarer than in southern populations, so far only heating has been correlated as an environmental stressor able to alter flowering time55,80. In our dataset, we observed that the best Posidonia homologous of the Arabidopsis CONSTANS was unlikely to be affected by ALAN at any of the study sites (see very low R2 value for this analysis in Table 1) (Fig. 4C). However, it should be noted that our experiments were not carried out within the flowering season, starting from September–October, therefore further analyses should be conducted specifically in this period to evaluate whether ALAN affects this aspect of Posidonia phenology. Studies investigating the potential role of ALAN on terrestrial plants highlighted effects on both growth-related and phenological variables in response to different intensities of light, including values of about 20 lx (such as those experienced relatively far from streetlamps), but not as dim as around 1–5 lx27. Organisms colonizing shallow coastal areas are usually exposed to a range of artificial light of smaller intensities in comparison to terrestrial species81,82. However, they are also adapted to low values of natural light1, which makes the potential effect of gradients of dim artificial light on marine autotrophs of scientific concern. Interestingly, a change in photosynthetic biomass have been observed in a planktonic cyanobacteria (Microcystis aeruginosa) under an irradiance of about 6.6 lx at night (80 nmol m2 s1)83. Our results showed a relation between an increase in intensity of dim artificial light and the pattern of expression of genes putatively related to the circadian clock in a marine angiosperm colonizing shallow coastal waters. In particular, they suggest an impairment in the mechanisms related to the control of gene expression during the day/night cycle (see PoELF3, PoZTL and PoCAB-6A, Figs. 3 and 4), including the perception of shorter blue wavelengths (see PoLUX1, Fig. 3), that are those penetrating deeper underwater and to which many marine species are particularly sensitive. The impairment of the expression pattern of genes regulated by the putative evening-loop genes could be associated with the reduced growth of the seagrass leaves. It is worth noting here that although the selected sites did not apparently differ in any other key biotic or abiotic factors (i.e., depth, wave-exposure, substratum type and inclination, grazing pressure) for the clock itself, other stressors potentially related to the gradient of proximity to the harbour city of Leghorn, might have influenced our results84,85, which cannot be, thus, considered as a definitive test of a cause-effect relationship between ALAN and the variables analysed in this study. Artificial light pollution is currently impinging on more than 20% of world’s coasts, and this percentage is set to increase with the growth rate of coastal human population. The emergent interest on potential impacts of ALAN on marine environments is showing the wide variety of effects related to the disruption of circadian rhythms on different organisms and habitats, from the molecular to the community level. Here we present evidence on the molecular and physiological effects of light pollution on a model seagrass species, Posidonia oceanica, which provides habitat and food for invertebrates, fish, mammals, birds and sea turtles, nutrient recycling and sediment stabilization20. Seagrasses are key ecosystem engineers in both subtidal and intertidal coastal areas worldwide. A variety of human-related stressors have been recognized as causes of the loss of coastal seagrass meadows and the services they provide, with a rate of disappearance comparable to those of mangrove forests and coral reefs86. As ALAN frequently co-occurs with a number of human-related disturbances in urban areas, our results pose key questions about the potential for light pollution to cumulate additively or to interact synergistically or antagonistically with these stressors87,88, in order to develop more efficient strategies to preserve seagrass ecosystems in coastal areas. Supplementary Information Supplementary Information 1. Supplementary Information 2. Supplementary Information 3. Supplementary Information 4. Supplementary Information 5. Supplementary Information 6. Supplementary Information 7. Supplementary Information 8. Supplementary Information 9. Supplementary Table S1. Supplementary Table S3. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-023-37261-3. Author contributions L.B.C. and F.B. conceived the ideas; L.B.C., F.B., E.M., F.L. and B.G. designed methodology; E.M. and F.B. identified the plant material in Maroccone, Scalinata and Rotonda sites; L.D.C., A.B., H.H., S.I., V.S., D.A.W., A.S., B.G., F.L., E.M. and L.B.C. collected the data; L.D.C., L.L., B.G., F.L. analysed the data; L.R. performed statistical analyses; L.D.C., F.L., E.M., F.B. and L.B.C. wrote the first draft of the paper. All authors contributed critically to the drafts and gave final approval for publication. Funding This research was supported by the Ph.D. Program in Agrobiosciences (to L.D.C., S.I. and V.S.) and Agrobiodiversity (H.H.) of the Scuola Superiore Sant’Anna, and by the University of Pisa under grant PRA_2017_19 to L.B.C.. F.L. and L.D.C. were also supported by the ERC Consolidator Grant number 101001320. Data availability All data generated or analysed in this study are included in this manuscript and its supplementary information files. Competing interests The authors declare no competing interests. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ==== Refs References 1. Davies TW Smyth T Why artificial light at night should be a focus for global change research in the 21st century Glob. Change Biol. 2018 24 3 872 882 10.1111/gcb.13927 2. Sanders D Frago E Kehoe R Patterson C Gaston KJ A meta-analysis of biological impacts of artificial light at night Nat. Ecol. Evol. 2020 10.1038/s41559-020-01322-x 33139919 3. Hölker F Moss T Griefahn B Kloas W Voigt CC Henckel D Hänel A Kappeler PM Völker S Schwope A Franke S Uhrlandt D Fischer J Klenke R Wolter C Tockner K The dark side of light: A transdisciplinary research agenda for light pollution policy Ecol. Soc. 2010 15 4 1 11 10.5751/ES-03685-150413 4. Falchi F Cinzano P Duriscoe D Kyba CCM Elvidge CD Baugh K Portnov BA Rybnikova NA Furgoni R The new world atlas of artificial night sky brightness Sci. Adv. 2016 2 6 e1600377 10.1126/sciadv.1600377 27386582 5. Brody, S. Circadian rhythms. In Brenner’s Encyclopedia of Genetics, 2nd edn. 10.1016/B978-0-12-374984-0.00251-5 (2013). 6. Bennie J Davies TW Cruse D Gaston KJ Ecological effects of artificial light at night on wild plants J. Ecol. 2016 104 611 620 10.1111/1365-2745.12551 7. Depledge MH Godard-Codding CAJ Bowen RE Light pollution in the sea Mar. Pollut. Bull. 2010 60 9 1383 1385 10.1016/j.marpolbul.2010.08.002 20736097 8. Davies TW Duffy JP Bennie J Gaston KJ The nature, extent, and ecological implications of marine light pollution Front. Ecol. Environ. 2014 12 6 347 355 10.1890/130281 9. Maggi E Serôdio J Artificial light at night: A new challenge in microphytobenthos research Front. Mar. Sci. 2020 10.3389/fmars.2020.00329 10. Schubert EF Kim JK Solid-state light sources getting smart Science 2005 308 5726 1274 1278 10.1126/science.1108712 15919985 11. Davies TW Coleman M Griffith KM Jenkins SR Night-time lighting alters the composition of marine epifaunal communities Biol. Lett. 2015 11 4 20150080 10.1098/rsbl.2015.0080 25926694 12. Navarro-Barranco C Hughes LE Effects of light pollution on the emergent fauna of shallow marine ecosystems: Amphipods as a case study Mar. Pollut. Bull. 2015 15 235 240 10.1016/j.marpolbul.2015.02.023 13. Bolton D Mayer-Pinto M Clark GF Dafforn KA Brassil WA Becker A Johnston EL Coastal urban lighting has ecological consequences for multiple trophic levels under the sea Sci. Total Environ. 2017 576 1 9 10.1016/j.scitotenv.2016.10.037 27780095 14. Underwood CN Davies TW Queirós AM Artificial light at night alters trophic interactions of intertidal invertebrates J. Anim. Ecol. 2017 86 4 781 789 10.1111/1365-2656.12670 28452048 15. Ayalon I de Barros Marangoni LFB Benichou JIC Avisar D Levy O Red Sea corals under Artificial Light Pollution at Night (ALAN) undergo oxidative stress and photosynthetic impairment Glob. Change Biol. 2019 25 12 4194 4207 10.1111/gcb.14795 16. Tamir R Eyal G Cohen I Loya Y Effects of light pollution on the early life stages of the most abundant northern red sea coral Microorganisms 2020 8 2 193 10.3390/microorganisms8020193 32023896 17. Maggi E Benedetti-Cecchi L Trophic compensation stabilizes marine primary producers exposed to artificial light at night Mar. Ecol. Prog. Ser. 2018 606 1 5 10.3354/meps12769 18. Maggi E Bongiorni L Fontanini D Capocchi A Dal Bello M Giacomelli A Benedetti-Cecchi L Artificial light at night erases positive interactions across trophic levels Funct. Ecol. 2020 34 3 694 706 10.1111/1365-2435.13485 19. Maggi E Bertocci I Benedetti-Cecchi L Light pollution enhances temporal variability of photosynthetic activity in mature and developing biofilm Hydrobiologia 2020 847 1793 1802 10.1007/s10750-019-04102-2 20. Duarte CM The future of seagrass meadows Environ. Conserv. 2002 29 2 192 206 10.1017/S0376892902000127 21. Steneck RS Graham MH Bourque BJ Corbett D Erlandson JM Estes JA Tegner MJ Kelp forest ecosystems: Biodiversity, stability, resilience and future Environ. Conserv. 2002 29 4 436 459 10.1017/S0376892902000322 22. Briggs W Rich C Longcore T Physiology of plant responses to artificial lighting Ecological Consequences of Artificial Night Lighting 2006 Island Press 389 412 23. Massetti L Assessing the impact of street lighting on Platanus x acerifolia phenology Urban For. Urban Green. 2018 34 71 77 10.1016/j.ufug.2018.05.015 24. Bennie J Davies TW Cruse D Bell F Gaston KJ Artificial light at night alters grassland vegetation species composition and phenology J. Appl. Ecol. 2017 55 1 442 450 10.1111/1365-2664.12927 25. Bennie J Davies TW Cruse D Inger R Gaston KJ Artificial light at night causes top-down and bottom-up trophic effects on invertebrate populations J. Appl. Ecol. 2018 55 6 2698 2706 10.1111/1365-2664.13240 26. Hey MH DiBiase E Roach DA Carr DE Haynes KJ Interactions between artificial light at night, soil moisture, and plant density affect the growth of a perennial wildflower Oecologia 2020 193 503 510 10.1007/s00442-020-04679-9 32533357 27. Sanders D Kehoe R Cruse D van Veen FJF Gaston KJ Low levels of artificial light at night strengthen top-down control in insect food web Curr. Biol. 2018 28 2474 2478 10.1016/j.cub.2018.05.078 30057304 28. Gobert S Cambridge MT Velimirov B Pergent G Lepoint G Bouquegneau JM Dauby P Pergent-Martini C Walker DI Biology of Posidonia Seagrasses: Biology, Ecology and Conservation 2006 Springer 387 408 29. Dattolo E Marín-Guirao L Ruiz JM Procaccini G Long-term acclimation to reciprocal light conditions suggests depth-related selection in the marine foundation species Posidonia oceanica Ecol. Evol. 2017 7 4 1148 1164 10.1002/ece3.2731 28303185 30. Boudouresque CF Bernard G Pergent G Shili A Verlaque M Regression of Mediterranean seagrasses caused by natural processes and anthropogenic disturbances and stress: A critical review Bot. Mar. 2009 52 5 395 418 10.1515/BOT.2009.057 31. Vassallo P Paoli C Rovere A Montefalcone M Morri C Bianchi CN The value of the seagrass Posidonia oceanica: A natural capital assessment Mar. Pollut. Bull. 2013 75 1–2 157 167 10.1016/J.MARPOLBUL.2013.07.044 23953894 32. Dodd AN Salathia N Hall A Kévei E Tòth R Nagy F Hibberd JM Millar AJ Webb AAR Plant Circadian clocks increase photosynthesis, growth, survival, and competitive advantage Science 2005 39 5734 630 633 10.1126/science.1115581 33. Shim JS Kubota A Imaizumi T Circadian clock and photoperiodic flowering in arabidopsis: CONSTANS is a Hub for Signal integration Plant Physiol. 2017 173 1 5 15 10.1104/pp.16.01327 27688622 34. Murchie EH Lawson T Chlorophyll fluorescence analysis: A guide to good practice and understanding some new applications J. Exp. Bot. 2013 64 13 3983 3998 10.1093/jxb/ert208 23913954 35. Aubé M Roby J Kocifaj M Evaluating potential spectral impacts of various artificial lights on melatonin suppression, photosynthesis, and star visibility PLoS One 2013 8 7 e67798 10.1371/journal.pone.0067798 23861808 36. Telesca L Belluscio A Criscoli A Ardizzone G Apostolaki ET Fraschetti S Gristina M Knittweis L Martin CS Pergent G Alagna A Badalamenti F Garofalo G Gerakaris V Louise Pace M Pergent-Martini C Salomidi M Seagrass meadows (Posidonia oceanica) distribution and trajectories of change Sci. Rep. 2015 5 1 14 10.1038/srep12505 37. El-Geziry TM Bryden IG The circulation pattern in the Mediterranean Sea: Issues for modeller consideration J. Oper. Oceanogr. 2010 3 2 39 46 10.1080/1755876X.2010.11020116 38. Ruocco M Barrote I Hofman JD Pes K Costa MM Procaccini G Silva J Dattolo E Daily regulation of key metabolic pathways in two seagrasses under natural light conditions Front. Ecol. Evol. 2021 10.3389/fevo.2021.757187 39. Oakenfull RJ Davis SJ Shining a light on the Arabidopsis circadian clock Plant Cell Environ. 2017 40 11 2571 2585 10.1111/pce.13033 28732105 40. Linde AM Eklund DM Kubota A Pederson ERA Holm K Gyllenstrand N Nishihama R Cronberg N Muranaka T Oyama T Kohchi T Lagercrantz U Early evolution of the land plant circadian clock New Phytol. 2017 216 2 576 590 10.1111/nph.14487 28244104 41. Altschul SF Gish W Miller W Myers EW Lipman DJ Basic local alignment search tool J. Mol. Biol. 1990 215 3 403 410 10.1016/S0022-2836(05)80360-2 2231712 42. Kumar S Stecher G Li M Knyaz C Tamura K MEGA X: Molecular evolutionary genetics analysis across computing platforms Mol. Biol. Evol. 2018 35 6 1547 1549 10.1093/molbev/msy096 29722887 43. Edgar RC MUSCLE: Multiple sequence alignment with high accuracy and high throughput Nucleic Acids Res. 2004 32 5 1792 1797 10.1093/nar/gkh340 15034147 44. Jones DT Taylor WR Thornton JM The rapid generation of mutation data matrices from protein sequences Bioinformatics 1992 8 275 282 10.1093/bioinformatics/8.3.275 45. Procaccini G Ruocco M Marín-Guirao L Dattolo E Brunet C D'Esposito D Lauritano C Mazzuca S Serra IA Bernardo L Piro A Beer S Björk M Gullström M Buapet P Rasmusson LM Felisberto P Gobert S Runcie JW Silva J Olivé I Costa MM Barrote Is Santos R Depth-specific fluctuations of gene expression and protein abundance modulate the photophysiology in the seagrass Posidonia oceanica Sci. Rep. 2017 7 42890 10.1038/srep42890 28211527 46. Loreti E Valeri MC Novi G Perata P Gene regulation and survival under hypoxia requires starch availability and metabolism Plant Physiol. 2018 176 2 1286 1298 10.1104/pp.17.01002 29084901 47. Ruocco M Brunet C Lorenti M Lauritano C D’Esposito D Riccio M Procaccini G Posidonia oceanica photoadaptation to the depth gradient Biol. Mar. Mediterr. 2012 19 63 64 48. Mazzuca S Björk M Beer S Felisberto P Gobert S Procaccini G Runcie J Silva J Borges AV Brunet C Buapet P Establishing research strategies, methodologies and technologies to link genomics and proteomics to seagrass productivity, community metabolism, and ecosystem carbon fluxes Front. Plant Sci. 2013 4 38 10.3389/fpls.2013.00038 23515425 49. Serra IA Lauritano C Dattolo E Puoti A Nicastro S Innocenti AM Procaccini G Reference genes assessment for the seagrass Posidonia oceanica in different salinity, pH and light conditions Mar. Biol. 2012 159 6 1269 1282 10.1007/s00227-012-1907-8 50. Livak KJ Schmittgen TD Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method Methods 2001 25 4 402 408 10.1006/meth.2001.1262 11846609 51. Zieman JC Methods for the study of the growth and production of turtle grass, Thalassia testudinum König Aquaculture 1974 4 139 143 10.1016/0044-8486(74)90029-5 52. Wood S Generalized Additive Models: An Introduction with R, Chapman & Hall/CRC Texts in Statistical Science 2006 Taylor & Francis 53. Zuur A Ieno EN Walker N Saveliev AA Smith GM Mixed Effects Models and Extensions in Ecology with R 2009 Springer 54. Assis J Tyberghein L Bosch S Verbruggen H Serrão EA De Clerck O Bio-ORACLE v2.0: Extending marine data layers for bioclimatic modelling Glob. Ecol. Biogeogr. 2018 27 3 277 284 10.1111/geb.12693 55. Balestri E Flowering of the seagrass Posidonia oceanica in a north-western Mediterranean coastal area: Temporal and spatial variations Mar. Biol. 2004 145 1 61 68 10.1007/s00227-004-1301-2 56. Orth RJ Harwell MC Inglis GJ Ecology of seagrass seeds and seagrass dispersal processes Seagrasses: Biology, Ecology and Conservation 2006 Springer 111 113 57. Zhang L Chen F Zhang X Li Z Zhao Y Lohaus R Chang X Dong W Ho SYW Liu X Song A Chen J Guo W Wang Z Zhuang Y Wang H Chen X Hu J Liu Y Qin Y Wang K Dong S Liu Y Zhang S Yu X Wu Q Wang L Yan X Jiao Y Kong H Zhou X Yu C Chen Y Li F Wang J Chen W Chen X Jia Q Zhang C Jiang Y Zhang W Liu G Fu J Chen F Ma H Van de Peer Y Tang H The water lily genome and the early evolution of flowering plants Nature 2020 577 79 84 10.1038/s41586-019-1852-5 31853069 58. Hsu PY Harmer SL Wheels within wheels: The plant circadian system Trends Plant Sci. 2014 19 4 240 249 10.1016/j.tplants.2013.11.007 24373845 59. Nohales MA Kay SA Molecular mechanisms at the core of the plant circadian oscillator Nat. Struct. Mol. Biol. 2016 23 12 1061 1069 10.1038/nsmb.3327 27922614 60. Wolfe KH Gouy M Yang YW Sharp PM Li WH Date of the monocot–dicot divergence estimated from chloroplast DNA sequence data Proc. Natl. Acad. Sci. U.S.A. 1989 86 16 6201 6205 10.1073/pnas.86.16.6201 2762323 61. Chaw SM Chang CC Chen HL Li WH Dating the monocot-dicot divergence and the origin of core eudicots using whole chloroplast genomes J. Mol. Evol. 2004 58 424 441 10.1007/s00239-003-2564-9 15114421 62. Olsen J Rouzé P Verhelst B The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea Nature 2016 530 331 335 10.1038/nature16548 26814964 63. Huang H Nusinow DA Into the evening: Complex interactions in the arabidopsis circadian clock Trends Genet. 2016 32 10 674 686 10.1016/j.tig.2016.08.002 27594171 64. Somers DE Schultz TF Milnamow M Kay SA ZEITLUPE encodes a novel clock-associated PAS protein from Arabidopsis Cell 2000 101 3 319 329 10.1016/S0092-8674(00)80841-7 10847686 65. Covington MF Panda S Liu XL Strayer CA Wagner DR Kay SA ELF3 modulates resetting of the circadian clock in Arabidopsis Plant Cell 2001 13 6 1305 1315 10.1105/tpc.13.6.1305 11402162 66. Zhao J Huang X Ouyang X Chen W Du A Zhu L Wang S Deng XW Li S OsELF3–1, an ortholog of arabidopsis EARLY FLOWERING 3, regulates rice circadian rhythm and photoperiodic flowering PLoS One 2012 7 8 e43705 10.1371/journal.pone.0043705 22912900 67. Sanchez SE Rugnone ML Kay SA Light perception: A matter of time Mol. Plant 2020 13 3 363 385 10.1016/j.molp.2020.02.006 32068156 68. Kim WY Fujiwara S Suh SS Kim J Han L David K Putterill J Nam HG Somers DE ZEITLUPE is a circadian photoreceptor stabilized by GIGANTEA in blue light Nature 2007 449 356 360 10.1038/nature06132 17704763 69. Dodd AN Kusakina J Hall A Gould PD Hanaoka M The circadian regulation of photosynthesis Photosynth. Res. 2014 119 1–2 181 190 10.1007/s11120-013-9811-8 23529849 70. Takahashi S Badger MR Photoprotection in plants: A new light on photosystem II damage Trends Plant Sci. 2011 16 1 53 56 10.1016/j.tplants.2010.10.001 21050798 71. Green RM Tingay S Wang ZY Tobin EM Circadian rhythms confer a higher level of fitness to Arabidopsis plants Plant Physiol. 2002 129 576 584 10.1104/pp.004374 12068102 72. Michael TP Salomé PA Yu HJ Spencer TR Sharp EL McPeek MA Alonso JM Ecker JR McClung CR Enhanced fitness conferred by naturally occurring variation in the Circadian clock Science 2003 302 5647 1049 1053 10.1126/science.1082971 14605371 73. Yerushalmi S Yakir E Green RM Circadian clocks and adaptation in Arabidopsis Mol. Ecol. 2011 20 6 1155 1165 10.1111/j.1365-294X.2010.04962.x 21232073 74. Kim JA Kim HS Choi SH Jang JY Jeong MJ Lee SI The importance of the circadian clock in regulating plant metabolism Int. J. Mol. Sci. 2017 18 12 2680 10.3390/ijms18122680 29232921 75. Fukao T Bailey-Serres J Plant responses to hypoxia—Is survival a balancing act? Trends Plant Sci. 2004 9 9 449 456 10.1016/j.tplants.2004.07.005 15337495 76. Geng Y Wu R Wee CW Xie F Wei X Chan PMY Tham C Duan L Dinneny JR A spatio-temporal understanding of growth regulation during the salt stress response in Arabidopsis Plant Cell 2013 25 6 2132 2154 10.1105/tpc.113.112896 23898029 77. Song Y Liu L Wei Y Li G Yue X An L Metabolite profiling of adh1 mutant response to cold stress in Arabidopsis Front. Plant Sci. 2017 7 2072 10.3389/fpls.2016.02072 28123394 78. Kim JM To TK Matsui A Tanoi K Kobayashi NI Matsuda F Habu Y Ogawa D Sakamoto T Matsunaga S Bashir K Acetate-mediated novel survival strategy against drought in plants Nat. Plants 2017 3 7 1 7 10.1038/nplants.2017.97 79. Copley TR Aliferis KA Kliebenstein DJ Jabaji SH An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease BMC Plant Biol. 2017 17 1 1 18 10.1186/s12870-017-1020-8 28049439 80. Ruiz JM Marín-Guirao L García-Muñoz R Ramos-Segura A Bernardeau-Esteller J Pérez M Sanmartí N Ontoria Y Romero J Arthur R Alcoverro T Procaccini G Experimental evidence of warming-induced flowering in the Mediterranean seagrass Posidonia oceanica Mar. Pollut. Bull. 2018 134 49 54 10.1016/j.marpolbul.2017.10.037 29102072 81. Tamir R Lerner A Haspel C Dubinsky Z Iluz D The spectral and spatial distribution of light pollution in the waters of the northern Gulf of Aqaba (Eilat) Sci. Rep. 2017 7 42329 10.1038/srep42329 28186138 82. Davies TW McKee D Fishwick J Tidau S Smyth T Biologically important artificial light at night on the seafloor Sci. Rep. 2020 10 12545 10.1038/s41598-020-69461-6 32719492 83. Poulin C Bruyant F Laprise MH Cockshutt AM Marie-Rose Vandenhecke J Huot Y The impact of light pollution on diel changes in the photophysiology of Microcystis aeruginosa J. Plankton Res. 2014 36 286 291 10.1093/plankt/fbt088 84. Salomé PA McClung CR Pseudo-response Regulator 7 and 9 are partially redundant genes essential for the temperature responsiveness of the Arabidopsis circadian clock Plant Cell 2005 17 3 791 803 10.1105/tpc.104.029504 15705949 85. Harmer SL The Circadian system in higher plants Annu. Rev. Plant Biol. 2009 60 357 377 10.1146/annurev.arplant.043008.092054 19575587 86. Waycott M Duarte CM Carruthers TJB Orth RJ Dennison WC Olyarnik S Calladine A Fourqurean JW Heck KL Hughes AR Kendrick GA Kenworthy WJ Short FT Williams SL Accelerating loss of seagrasses across the globe threatens coastal ecosystems Proc. Natl. Acad. Sci. U.S.A. 2009 106 30 12377 12381 10.1073/pnas.0905620106 19587236 87. Halpern BS Walbridge S Selkoe KA Kappel CV Micheli F D'Agrosa C Bruno JF Casey KS Ebert C Fox HE Fujita R Heinemann D Lenihan HS Madin EMP Perry MT Selig ER Spalding M Steneck R Watson R A global map of human impact on marine ecosystems Science 2008 319 5865 948 952 10.1126/science.1149345 18276889 88. Claudet J Fraschetti S Human-driven impacts on marine habitats: A regional meta-analysis in the Mediterranean Sea Biol. Conserv. 2010 143 9 2195 2206 10.1016/j.biocon.2010.06.004