Italian Botanist 9: |—19 (2020) doi: 10.3897 /italianbotanist.9.46426 RESEARCH ARTICLE YQ g Italian Botanist http://italianbotanist.pensoft.net Assessment of the ecological structure of Posidonia oceanica (L.) Delile on the northern coast of Lazio, Italy (central Tyrrhenian, Mediterranean) Valentina Gnisci', Selvaggia Cognetti de Martiis', Alessandro Belmonte’, Carla Micheli*, Viviana Piermattei!, Simone Bonamano!, Marco Marcelli! | Department of Ecology and Biology (DEB), Laboratory of Experimental Oceanology and Marine Ecology, University of Tuscia, Porto di Civitavecchia, I-00053 Civitavecchia (RM), Italy 2 ENEA Italian Natio- nal Agency for New Technologies, Energy and Sustainable Economic Development. Department of Energe- tic Technologies, Laboratory of Biomass and Energetic Biotechnologies. Research Center Casaccia Box 2400, I-00123 Rome, Italy Corresponding author: Carla Micheli (carla.micheli@enea.it) Academic editor: Giovanni Astuti | Received 12 September 2019 | Accepted 7 December 2019 | Published 9 January 2020 Citation: Gnisci V, Martiis SC, Belmonte A, Micheli C, Piermattei V, Bonamano S, Marcelli M (2020) Assessment of the ecological structure of Posidonia oceanica (L.) Delile on the northern coast of Lazio, Italy (central Tyrrhenian, Mediterranean). Italian Botanist 9: 1-19. https://doi.org/10.3897/italianbotanist.9.46426 Abstract The ecological structure of Posidonia oceanica (L.) Delile meadows was evaluated on the northern coast of Lazio, Italy (central Tyrrhenian, Mediterranean sea). This is an infra-littoral zone with a wide range of an- thropogenic activities and high geo-morphological variability, which reflects heterogeneity in shoot densi- ty, leaf morphology and biomass in fragmented patches. Genetic variability in populations corresponds to the formation of 3 sub-clusters, in the diverse impacted zones (north, centre and south), being correlated to the geographical distance between sites. AMOVA estimated a high genetic variation showing 43.05% individual differences within populations with a marked differentiation among the populations (56.9%) indicated by Fst value (0.57). These results revealed the role of the genetic structure of seagrasses for de- termining selectivity of fragmented habitat, in response to natural drivers. They showed that site-specific self-recruitment is related to biodiversity capacity and to the geo-morphological characteristic of the coast. Keywords AMOVA, Functional descriptors, Genetic variation, RAPD Copyright Valentina Gnisci et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 2 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) Introduction Coastal areas are characterized by environmental disturbances due to both natural process and anthropogenic activities, with consequent impacts on marine ecosystems (Benoit and Comeau 2005). Posidonia oceanica (L.) Delile meadows are considered as indicators of these impacts (Pergent et al. 1995; Pergent-Martini et al. 2005; Martinez- Grego et al. 2008; Lopez y Royo et al. 2011; Bennett et al. 2011; Personnic et al. 2014), which may influence their natural structural variability, causing a general de- cline of biomass, density and coverage (Duarte 2002; Orth et al. 2006; Boudouresque et al. 2009; Marba et al. 2014). In the Mediterranean Sea, P oceanica forms wide- spread ecosystems with high levels of biodiversity and productivity (Boudouresque et al. 1984; 2006; 2012) and, like other seagrasses, it is important from the ecological, geological and economic point of view (Costanza et al. 1997; Spalding et al. 2003; Vassallo et al. 2013). Moreover, this plant colonizes different types of substrates (e.g. sand, rocky, etc.) and generates meadows that occur in a wide range of morphotypes, ranging from continuous beds on the seafloor to fragmented patches of different shape and size (Molinier and Picard 1952; Den Hartog 1970; Semroud 1993; Boudouresque et al. 2012; Bonhomme et al. 2015). The monitoring of structural and functional descriptors of P oceanica meadow health status (Buia et al. 2004) provides important information about the vitality and dynamic of meadows (Pergent-Martini et al. 2005), as well as the human influence on the environment (Montefalcone 2009). Recently, a new suite of descriptors has been proposed to assess the ecological status of P oceanica meadows (Rotini et al. 2013), taking into account the physiological features of the plants by means of biochemical (Arnold et al. 2012) and genetic markers (Micheli et al. 2005, 2015; Procaccini et al. 2007; Macreadie et al. 2014) as standard indicators. The objective of this work was to assess the ecological structure and genetic patterns of P oceanica meadows along the northern coast of Lazio to contribute to the conser- vation of seagrass by future management activity. For this aim, we used functional descriptors (phenotypic variables, including biomass and density) and RAPD (Ran- dom Amplified Polymorphic DNA) molecular markers. The PCR-based RAPD marker technique was recently used in seagrasses for estimating the influence of environmen- tal disturbances on phenotypic variables and genetic diversity of populations (Kim et al. 2019), and for distinguishing genetic variation of populations at individual level (Dilipan et al. 2017). At multiple spatial scales, RAPD analyses was conducted to assess the genetic variability of seagrass, evidencing genetic differences closely correlated with coastal currents (Jones et al. 2008). In previous studies, RAPD allowed us to examine the genetic structure in meadows located in the Ligurian Gulf (North Tyrrhenian coast) that showed a decrease in genetic diversity along an anthropogenic disturbance gradient at small scale (Micheli et al. 2005), revealing the role of currents (Micheli et al. 2010b; Rotini et al. 2013). Furthermore, we were also able to evaluate changes in genetic struc- ture of the meadows and their resilience over a decade (Micheli et al. 2012). Here we used RAPD to detect the genetic structure of P oceanica meadows in the fragmented habitat of the northern coast of Lazio, in order to compare it to that of Posidonia oceanica (L) Delile 3 other Tyrrhenian populations previously studied. Such information was detected in an area where natural drivers occur at exceptional conditions and where no similar works have been carried out before. Materials and methods Study area The study area between Marina di Tarquinia and Santa Severa (Lazio, Italy, Mediterranean Sea; Fig. 1) is located within two physiographic units, one extending from Monte Argen- tario to Capo Linaro, the other from Capo Linaro to Capo d’ Anzio. Here the appearance of the coast is intimately related to the morphology of the seabed as the isobaths show a very uneven underwater sea-floor. In fact, according to Anselmi et al. (1978), different coastal geo-morphological types (coastal morpho-type) are present in the study area. The northern part hosts the Mignone river floodplain, which is characterized by small sandy beaches and a rocky coastal terrace. From Civitavecchia to Santa Marinel- la, the coastline is dominated by the “Tolfa Mountains”, which form a promontory characterizing the coastal morpho-type (terraces coasts). This promontory separates the southern physiographic unit from the northern one, and is crossed by small streams (e.g. Marangone stream) with local continental contributions. The southern part pre- sents a small portion of coastal terraces and beaches with a prevalently sandy coast. This coastal area is characterized by the presence of the littoral currents having a prevailing direction from south to north following a coastal local dynamic connected to high geo-morphological variability of the sea bottom, which can generate barriers to gene flow. In fact, in the study area the prevailing wind events come from to the southeast (data provided from the weather station of C-CEMS), inducing a sea cur- rent direction to north (Bonamano et al. 2015; 2016). However, P oceanica flowers appeared in the meadows, and seeds were found in the central zone, near the Maran- gone stream, during November 2013, as was observed in other sites of the area studied (Gnisci 2014; Cognetti de Martiis 2016). The area is also characterized by a very large port (Civitavecchia harbour), two im- portant power plants located in the northern part of Civitavecchia, and a dense urban environment constituted by the municipalities of Civitavecchia and Santa Marinella, which altogether form a single urban aggregate (Fig. 1). Field and sampling work Field work activities were carried out during the late spring 2013, along 40 km of coastline. By SCUBA diving, shoots were collected into 18 sampling areas (6 shoots per sites, 3 sites per each station) from May 3“ to June 19" 2013 (Fig. 1, Table 1). For field work, we considered sufficient this temporal range as meteorological conditions 4 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) Table |. Sampling sites, geographical coordinates, depth and three types of substrate of Posidonia oce- anica. | = Rock, II = Rock+Sand, III = Sand+Matte. Stations Geographical coordinates Depth (m) Substrate Type stl 42°13'04.2"N, 11°41'25.2"E 10.0 Rock, Sand II St2 42°12'19.8"N, 11°42'16.2"E O85 Sand, Matte III St3 42°11'04.2"N, 11°42'25.8"E 13.4 Rock, Sand II St4 42°08'49.2"N, 11°43'52.8"E 10.3 Rock I St 42°08'04.8"N, 11°44'25.8"E 5.5 Rock, Sand II St6 42°05'05.0"N, 11°47'37.2"E 4.9 Rock I St7 42°04'35.1"N, 11°47'57.2"E 75 Sand, Matte III St8 42°04'15.1"N, 11°48'06.4"E 12.0 Rock I St9 42°03'27.9"N, 11°48'35.9"E 10.0 Rock I St10 42°02'59.2"N, 11°48'57.8"E 10.3 Rock, Sand II Stll 42°02'23.7"N, 11°48'59.8"E 11.0 Rock I St12 42°01'58.8"N, 11°48'55.8"E 8.9 Rock I St13 42°01'16.8"N, 11°50'40.8"E 10.9 Rock I St14 42°01'43.2"N, 11°51'16.2"E 10.0 Rock, Sand II Stl5 42°01'58.8"N, 11°52'36.0"E 13.0 Sand, Matte III St16 42°01'52.8"N, 11°54'52.8"E 13.2 Sand, Matte III Stl7 42°00'58.2"N, 11°56'40.8"E 9.4 Sand, Matte III St18 42°00'24.2"N, 11°55'00.7"E 13.4 Rock, Sand II : = | | 42.24 7 ae arta river » %* . . nee te qMarina di Tarquinia : : we 42.2 si ie manta ae 42.16— 42.12 Lazio as - ne 42.08" %, ° , one N = mare” eS a | Tr Santa Marinella ° ‘| Power plant : » ai @ Sampling station > [sm h N Habitat identified as a priority by : Habitat Directive (92/43/EEC) 0 5km 47 11.65" 11.7° 11.75" 11.8° 11.85° 11.9 11.95° 1r Figure |. Map of the study area in northern coast of Lazio (Italy, central Tyrrhenian, Mediterranean) with the sampling sites of Posidonia oceanica. Posidonia oceanica (L) Delile 5 were constant in the study area, and because P oceanica meadows undergo changes in a larger time scales. In each station, three sites were selected to estimate density and coverage of the meadow. By hierarchical sampling design, shoot density was calculated as shoot number/m? counting the number of shoots in nine random squares (40 x 40 cm) (Buia et al. 2004; Panayotidis et al. 1981). On the seafloor, the percentage in coverage of P oceanica was estimated counting the number of sub-squares occupied by the plants in three gridded squares of 1 x 1 m (Boudouresque et al. 2007). At the same time, eighteen shoots were randomly collected in each station and stored in sea-water at 4 °C, for laboratory analyses. Morphological measures and biomass In the laboratory, leaves of the shoots were washed in distilled water and biometric variables, such as number, length and width of juvenile, intermediate and adult leaves per shoot were measured in each station, according to Giraud’s classification (Giraud 1977). Leaf Area Index (LAI) average was calculated according to Buia et al. (2004). Furthermore, the leaves were cleaned (by a razor blade to remove epiphytes) and dried at 80 °C for 72 hours, to obtain the dried weight, i.e. the dry biomass (g dm/shoot). RAPD genetic analysis Genetic analyses of P oceanica populations were performed on the same shoots col- lected in each of the 18 sampled stations. According to Micheli et al. (2005), the young leaves were washed in distilled water, frozen in liquid nitrogen at -180 °C and stored at -80 °C before the extraction of genomic DNA, as reported in Dellaporta et al. (1983). Subsequently, DNA was amplified by TAQ Polymerase (Applied Biosystems), using PCR (Polymerase Chain Reaction) conditions similar to those described by Echt et al. (1992) with some modification involving reaction buffer (100 mM Tris-HCL, 500 mM KCL and 25 mM MegC]I2). Amplification reactions were carried out in a ther- mal cycler (Perkin-Elmer) with the following temperature program: 5 min denatura- tion cycle at 94 °C and, subsequently, 40 cycles of 30 s at 37 °C annealing temperature and 6 min at 72 °C synthesis temperature. For each DNA (1 pg/L) amplification, the AmpliTagq DNA Polymerase (0.4 pg/uL) and 10 different easily repeated oligonucleo- tides (10 mM) were used for their capacity to discriminate bands and scoring them as present/absent. Amplification products were separated by agarose gel electrophoresis (1.4%), stained with ethidium bromide, visualized on UV light and photographed. The different DNA fingerprints obtained were analyzed for the presence (1) or absence (0) of the bands, in order to determine the percentage of polymorphism calculated as the number of polymorphic bands out of the total number of bands. 6 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) Statistical analyses Morphological and structural variables of P oceanica meadows were statistically ana- lyzed by means of PAST (Hammer et al. 2001) and R Development Core Team (2008) software. Three way Analysis of Variance (ANOVA) was employed to test the effect of three main factors, i.e. location, depth and substrate typology, on meadow structure parameters and leaf morphology. Furthermore the non-parametric Kruskal-Wallis test was performed to check the level of significance (p < 0.05, p < 0.01 or p < 0.001) of the main factors (location, depth and substrate typology) effect on number of juvenile, intermediate and adult leaves in each shoot, which data no present a normal distri- bution. The Mann-Whitney post-hoc pairwise comparison was performed to valuate significant differences (p < 0.05, p < 0.01, p < 0.001) among different sampling area (North, Centre and South). Intra- and inter-population genetic distance has been de- termined using the NTSYS 2.0 program (Rohlf 1993) generating a similarity matrix (Nei and Li 1979). Cluster analysis was performed on the similarity matrix with the SAHN program using UPGMA (Un-weighted Pair Group Method with Arithmetic Mean) generating the dendrogram. A cophenetic correlation has been used to measure goodness of fit for the cluster analysis (r) by comparing it with the original pairwise dis- tance matrix (Euclidean squared distance) using a Mantel test (MXCOMP program). Then, Principal Coordinate Analysis (PCoA, NT-SYS software; Rohlf 1993) was performed to elucidate the distribution of the samples by deriving genetic distances. Subsequently AMOVA (Analysis of Molecular Variance) was performed with AR- LEQUIN 3.5 program (Excofher and Lischer 2010; http://cmpg.unibe.ch/software/arle- quin3) to elucidate the genetic variation and relationships between individuals within and among populations. Using AMOVA, we estimated the number of genetically homogene- ous groups of individuals based on their genotypes, at multiple loci, using the Bayesian ap- proach. The method assumes a model in which there are K populations, each characterized by a set of allele frequencies at each locus. For each individual, the software computes the proportion of the genome that can be assigned to the inferred K populations with a high probability. Because RAPD markers are dominant, we assumed that each band represented the phenotype at a single biallelic locus (Williams et al. 1990). The heterogeneity of band frequencies across all populations was tested using the Fixation index (Fst) in ARLEQUIN software. Fst is a measure of population differentiation (Meirmans 2006), and represents the most frequently used method to analyze population structure concerning gene flow. Results Leaf morphology, biomass, density and coverage of the meadows Results of leaf morphology (number of leaves/shoot, width and length), Leaf Area Index (LAI), biomass, shoot density and coverage of the plants collected in the 18 sampling stations, are reported in Table 2. Posidonia oceanica (L) Delile VVCL LOUIE + OVEY CSTE SE V84+8°LY 1004+ c9 LLL + VSP Cl + 9°86 VC+IUII C9 + L81 66TH HY 8S + 9°97 CST FEED OL +9°C1 10°0 + €6 VEL + CCT G64 SL 8'ET+ TOE VY +9C EST+LIV v9 + 107 (%) a3vJaA07) vO'LY CEL + 9LT GSH + SBTC 9°69 + 6'S9E G°8S + COST CSS +76LI 9°SS+ S'8CYV 9°8Z + 6061 9°8L + 6 CLE 8°8Z + 8°S9T ['S9 + 1°9€7 EOS + CET 69 + LTP 1°86 + C6CE E84 LTE CSYIF S'ESO ESY + T18T CW + LCT Tes + 6CLI LI8 + 9°€ZT (:u1/s1004s) Ayisuaq oo0yYs TeOE 170 + 020 170 +920 €C'0 + C80 IT'0 + 070 0c'0 + 670 Se"0 + ZZ'0 €C'0 + €9'0 60°0 + 670 c£'0 + 68°0 0¢'0 + £€8°0 £70 + 99'0 LE'0 + $60 670 + £01 670 + $S'0 170 + 6€°0 67'0 + 69°0 940 + £60 LC'0 + £80 6¢°0 + £90 (Gooys/up 3) ssewolg SC8P LV + OS'Y VL+LY VC+HOL 60+ 81 TL+9¢ VE+C6 VL+67C CL + 8€ CTHIS LI+8€ SI FEV CI+0°€ CC + 89 CCT +OYV CVEGL SL +S 81+ 0S LI+¢V SIL+ TE (ct /, Ut) (IVT) Xpuy easy FeoT c8°0C O€'L + 90°SE LL + GLE LST + ESE Cll + C67 VEL + OCE SOL + BEV LCI + SOE VCI+OTE CLI + OSE CCL FOE C8I + SOF GVI + CEE LYL+6PVE CYL + CVC SCL + EV LYL+CLE C6L+ ZY ULI + 69 Sl + LOY (m9) sug] Feo] GL’S S0°0 + 06°0 €l'0 + 160 60°0 + 88°0 80°0 + 08°0 I1'0 + 06°0 80°0 + 76'0 60°0 + C6'0 €C'0 + 06'0 80°0 + 76'0 yV10 + 780 80°0 + 68°0 60°0 + 260 60°0 + Z6'0 O10 + 98°0 60°0 + 98°0 L0°0 + 760 cl'0 + 98°0 80°0 + 160 O10 + 06'0 (03) PEM Jeo] IyZl €O'l + C6'S 60 +69 Te+89 CL +9°S CSI+LY CT HOS SI + 6S LI+eS 60+79 60 + 6S CLH+ZS CL+Z9 ST + 19 CL + 6S VC+O9 CI + 6S LI +9°9 TL + 6S SLES (ot) OOUS TENT UONLITeA JO JUDIOYYI07) (CS) UonelAsq prepurig 7 ursyy 81S LVS 91S STS VIS €13$ CTS TTS OLS 63S 83S LS 9$ GS vs €1$ (aN) Ls suoners ‘ous Surjdures yoeva url poraqjoo syurld y21uvaI0 PIUOpISOg JO IBeIZAOD pue Aisuap IOOYs ‘sseusorg ‘([y"]) XOpuy vary Jeary ‘(yISus] pue ‘YIprM “aquinu) saavoy Jo somnseous [ed1sojoydioyy *7 aIqeL 8 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) The presence of highly fragmented meadows (evidenced in Suppl. material 1: Fig. S1), in a large range of depths (from 4.9 m to 13.4 m. Table 1), was confirmed by the high variability of coverage found among stations, ranging from 6 to 98 % (St 16 and St 14 respectively), and by the coefficient of variation (72.4%) (Table 2). Three-way ANOVA (Suppl. material 3: Table S1), was tested to detect significant factors in a multifactor model, in which there is one dependent variable (shoot density of P oceanica) and three independent variables (depth, location and substrate typology). The three way ANOVA showed differences in sea-bottom coverage (P < 0.05) due to different substrate types (Table 3). In all areas, the mean shoot density estimated was 276 = 133 SD (number of shoot/m? = Standard Deviation), ranging from 141 + 62 SD to 653 + 145 SD (St 8 and St 5, respectively). Biomass mean value 0.70 + 0.213 SD (g dm/shoot) showed the highest values in the stations at intermediate depths, ranging from 1.03 = 0.49 SD to 0.40 + 0.11 SD. All the shoots showed similar values in leaf number (5.92 + 1.03 SD mean value) and width (0.9 = 0.05 SD mean value). Due to the differences in depth and water transparency of the 18 sites, we found a highly variable leaf length with values ranging from 24.24 = 14.22 SD to 49.44 + 17.11 SD (mean values of 35.05 + 7.3 SD), showing statistically significant differences in the three populations (ANOVA, P < 0.05) (Table 2, Table 3, Suppl. material 3: Table S1). In the Table 3, the multi-factors three-way ANOVA performed on shoot density of P oceanica, showed no significant difference among the sites located in the north, central and south areas (ANOVA, p = 0.1868); on the contrary significant differences were found among stations located both on different depths (ANOVA, p < 0.001) and on different substrates (ANOVA, p < 0.05). The non-parametric Kruskal-Wallis test was performed to verify the presence of significant differences in number of juvenile and intermediate leaves among the three sampling location (North, Centre and South). The test highlighted significant dif- ferences in the number of juvenile (p < 0.001) and intermediate (p < 0.001) leaves. Following Kruskal-Wallis results, the Mann-Whitney post-hoc pairwise comparison test was performed; it highlighted the lowest values of the juvenile leaf number, and the highest values of the intermediate leaf number in the central area compared to the other two locations. ‘The adult leaves number showed no significant differences among the different locations (Kruskal-Wallis, p = 0.0525) (Suppl. material 2: Fig. S2). Genetic variability of the population By UPGMA dendrogram (Fig. 2), based on SAHN clustering of the data, the simi- larity coefficient (Nei’s index) among the P oceanica samples was 0.61 + 0.01 based on 140 RAPD total bands (mono- and polymorphic) with fragments ranging in size from 0.3 kb to 3.8 kb. The percentage of polymorphism obtained in the population was 90.14% (Table 4). Posidonia oceanica (L) Delile 9 Table 3. Three way Analysis of Variance (ANOVA) applied to three main factors (i.e. location, depth, and substrate type) on meadow structure parameters and leaf morphology. Sum Sq. d.f. Mean Sq. F Prob > F Depth 226471 2 1NH3239.5 W256 p < 9.001 Location 31885.8 2, 15942.9 1.77 0.1868 Substrates 76990.4 2 38495.2 4.27 p < 0.05 Error 288410 32 9012.8 Total 896383.4 50 Table 4. RAPD marker sequences, amplification range and total polymorphism detected in the studied Posidonia oceanica meadows. Primer Sequence Amplification range Fragment amplification Tot Polymorphism (%) polymorphic monomorphic BY11 5’-ATCCACTGCA-3’ 0.3-2.4 Kb 18 0 18 100.00 BY12 5’-GGTCGCAGGC-3’ 0.3-3.8 Kb 12 2 14 85.71 BY13. = 5’--CCTTGACGCA-3’ 0.5-3.8 Kb 13 2 15 86.67 BY15 = 5’-CT'CACCGTCC-3’ 0.73-2.8 Kb isl 1 12 91.67 DN4 5’-GTCGTGCTAT-3’ 0.-2.4 Kb 9 2 11 81.82 DN5 5’-CCGACGGCAA-3’ 0.3-1.75 Kb Die, 0 Doak 100.00 DN6 = 5’-TGGACCGGTG-3’ 0.5-3.0 Kb 8 1 9 88.89 UB24 5’-GGGTGAACCG-3’ 0.75-2.25 Kb 12 0 12 100.00 UB26 =5’-CGCCCCCAGT-3’ 0.73-2.75 Kb 8 4 12 66.67 UB28 5’-GCTGGGCCGA-3’ 0.5-1.85 Kb 15 0 15 100.00 Average pergentage of polymorphism (%) 90.14 Mantel test, which compares Nei’s distance and cophenetic matrices, returned sta- tistically significant values, with a matrix correlation of r = 0.95 (normalized Mantel Sta- tistic Z); t = 25 (Mantel T-test) and p = 1.00000 (Probability random Z < observed Z). PCoA analysis showed that the individuals sampled in the 3 different impacted zones fell into 3 distinct groups (north, centre and south populations) (Fig. 3). The first two axes accounted for 57.8% of variation, with axis 1 and 2 explaining 39.7 % and 18.1% of the variation, respectively. AMOVA (Suppl. material 4: Table S2) showed 43.05% individual differences within populations with a marked differentia- tion among the populations (56.9%) indicated by Fst value (0.57). Discussion In the meadows, the genetic structure of P oceanica populations was clustered into three main groups located into three different zones (Fig. 2). The first population has been identified in the northern area, corresponding to sampling sites between Civi- tavecchia harbour and Marina di Tarquinia, an area mainly dominate by the presence of the Mignone river floodplain. As observed in UPGMA analysis, meadows located in this northern area are divided in two sub-populations, which are crossed by Mignone 10 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) river. The second population (in the centre) corresponds to the sampling stations be- tween Civitavecchia harbour and Capo Linaro with the presence of low cliffs and a very heterogeneous sea floor, with a prevalently rocky bottom. ‘The third population, located in the southern part of the study area, from Capo Linaro to Santa Severa, is characterized by a prevalent sandy sea floor. By AMOVA, a high level of genetic diversity was recognized within (43.05%) and among the populations (56.95%) and was confirmed by PCoA (Fig. 3) when the vari- ation was applied to genetic distance among individuals (57.8%). In this highly heterogeneous coastal area, we found 90.14% of total polymor- phism ranging from 66.67 % to 100 % in the individuals analyzed (Table 4). This high percentage of polymorphism demonstrates the high variability displayed by each individual, which is also evidenced by the UPGMA cluster analysis (Figure 2) showing a low similarity value (0.52) within each population. In the UPGMA cluster analysis, the 0.52 similarity value is lower than the one found for Monterosso meadow (0.66), where natural and anthropogenic pressure were low, and lower than the one found at the Mediterranean basin scale (0.81) (Micheli et al. 2005). Moreover, in the same sub-cluster (Figure 2) are grouped the individuals of popu- lation with the same morphological features of the leaves linked to their genetic capac- ity to fit the same environment. In the Table 2, the structural and functional descriptors have confirmed a high variability in all sampling areas, due to the different depth, which may particularly influence the morphology of leaves. According to other authors (Di Maida et al. 2013; Balestri et al. 2004; Pergent et al. 1995), the meadows that we found are frag- mented, discontinuous and are characterized by a wide range of shoot density and substrate coverage which reflect a spatial heterogeneity of P oceanica meadows, prob- ably due to the complex nature of the seafloor (Tables 1, 2). Our results highlight the influence of natural drivers, as the geo-morphological features of the fragmented habitat, in the connectivity of seagrass populations. The heterogeneity patterns of the meadows that we found along the northern coast of Lazio, describe the ability of plants to respond to the environment. This depends on the availability of genetic variation improving plant growth, development and resilience (Hughes and Stacho- wicz 2004; Reusch et al. 2005). At this regard, Jahnke et al. (2015) suggested not only a greater resistance and resilience of individuals of higher genetic and genotypic diversity under disturbance, but also that the high diversity in high impacted sites reflects the mismatches with respect to pre-environmental impact conditions. espe- cially because flowering and sexual recruitment are seldom observed. In fact, by the genetic diversity metrics study, they observed that the absence of low and medium levels of genetic variation in impacted locations were probably due to local extinc- tions of individuals that already exceeded their resistance capacity. Therefore, they suggested to use this valuable tool for restoration and mitigation projects, because the exceptional longevity of individuals (i.e. P oceanica population) could create a temporal mismatch that may falsely suggest good meadow health status, while gradual deterioration of allelic diversity might go unnoticed (Jahnke et al. 2015). Posidonia oceanica (L) Delile 11 = North “i & Centre | == South Figure 2. UPGMA analysis of Posidonia oceanica meadows. South © 148 @_18A 1 © vat 1781% ? pore % 15A * 188 © 164 Coordinate 2 Centre @ ec 7B @ 1201® 10C1 7A 12810@ © 9c zc @ Se_24, 11B1 @ 2 @ 11A1 Coordinate 1 Figure 3. Principal Coordinate analyses (PCoA) based on 140 RAPD bands of XX individuals of Posidonia oceanica. According to the above studies, our genetic data have highlighted the influence of different habitat conditions on plants found in the three different areas, which parallel the high variability of all structural and functional parameters (Table 2). The incidence 12 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) of multiple environmental cues led plants to a more plastic response, differing among individuals, and the genetic distance between populations correlates to the distance between sites showing different environmental conditions (substrate, marine dynamic, nutrients, etc.) and subject to both natural drivers and different stressors. Conclusion The first investigation on the ecological/genetic structure of P oceanica growing along the coasts of the northern Lazio has called the attention on the higher genetic variability of these meadows than others previously studied (Micheli et al. 2005; 2010b; 2012; Rotini et al. 2011; 2013). The observed genetic diversity and differentiation (Fst = 0.57) found among the P oceanica populations suggest an adaptation capacity of the species to share different environments (rocky and sandy). This reflects in the formation of dif- ferent local meadows with a spatial heterogeneity distribution along the studied area. UPGMA cluster of genetic distance revealed a high variability of seagrass population in the relatively close meadows as the result of a low level of gene flow among them. By using a dispersal simulation analysis, Serra et al. (2015) suggested that such pattern of genetic structure probably reflects historical events of recolonization from relict glacial areas and of past vicariance, which seems to persist because of the low connectivity among populations via marine currents. Whilst Jahnke et al. (2015) observed that the absence of low and medium levels of genetic variation in impacted locations are proba- bly due to local extinctions of individuals that already exceeded their resistance capacity. In the Mediterranean, P oceanica plays a crucial role for the ecosystem productivity through favoring biodiversity, revealing a whole range of intra-specific levels of diver- sity from complete clonality to high variability (Ribeiro et al. 2016). In this littoral zone of the northern coast of Lazio, natural drivers like the prevailing current direction towards north (Bonamano et al. 2015; 2016) could affect seed transport favoring the colonization of new habitats whilst the incidence of multiple environmental impacts and the intra-meadow genotypic heterogeneity drive the plant to a more plastic re- sponse among individuals (Micheli et al. 2015; 2012). Posidonia oceanica is considered a vulnerable species, largely susceptible to pertur- bations, and is thus included within Habitats Directive (92/43/EEC), which supports the Natura 2000 network of European protected areas. From a conservation perspec- tive, and considering future increase in both anthropogenic and climatic pressures (Mannino 2018), the assessment and management of environmental quality require an adequate knowledge of the ecosystems that can be obtained only by multidiscipli- nary investigations (Cozza et al. 2019; Domina et al. 2018). Therefore further stud- ies focused on the intra-variability of the meadows, by analyzing the allelic richness (Jahnke et al. 2015), are needed. These studies should include the changes of private (neutral) allelic in order to control the population demography vs genetic relationship stability necessary to reinforce their resilience as we have observed in other protected areas in the Mediterranean (Micheli et al. 2015; 2012). Finally, since seagrasses are Posidonia oceanica (L) Delile 13 globally declining, and their loss is often due to synergies among stressors (Ceccherelli et al. 2018), the ecological implication of the results obtained here, such as the genetic features and the environmental conditions, could be further aspects to be considered in the monitoring surveys. Acknowledgments Grateful acknowledgments are made to the Port Authority of Civitavecchia, divers and colleagues of Experimental Oceanology and Marine Ecology laboratory for sampling activities. Many thanks to University of Tuscia that has funded this research and VG and SCDM PhDs. We thank the anonymous reviewers whose comments and sugges- tions helped improve this manuscript. References Anselmi B, Brondi A, Falchi G, Ferretti O (1978) Lineamenti granulometrici e mineralogici generali dei sedimenti fluviali e costieri del territorio italiano. Memorie della Societa Geo- logica Italiana 19: 307-314. Arnold T, Mealey C, Leahey H, Miller AW, Hall-Spencer JM, Milazzo M, Maers K (2012) Ocean acidification and the loss of phenolic substances in marine plants. PLOS One 7(4): e35107. https://doi.org/10.1371/journal.pone.0035 107 Balestri E, Benedetti Cecchi L, Lardicci C (2004) Variability in patterns of growth and mor- phology of Posidonia oceanica exposed to urban and industrial wastes: contrasts with two reference locations. Journal of Experimental Marine Biology and Ecology 308: 1-21. https://doi.org/10.1016/j.jembe.2004.01.015 Bennett S, Roca G, Romero J, Alcoverro T (2011) Ecological status of seagrass ecosystems: an uncertainty analysis of meadow classification based on the Posidonia oceanica multivari- ate index (POMI). Marine Pollution Bulletin 62: 1616-1621. https://doi.org/10.1016/j. marpolbul.2011.06.016 Benoit G, Comeau A (2005) A sustainable future for the Mediterranean. In: The Blue Plan’s Environment and Development Outlook. Earth scan, London, 450 pp. Bonamano S, de Mendoz FP, Piermattei V, Martellucci R, Madonia A, Gnisci V, Mancini E, Fersini G, Burgio C, Marcelli M, Zappala G (2015) Mathematical models supporting the monitoring of Civitavecchia harbour (Rome). WIT Transactions on Modelling and Simu- lation 59: 443-454. https://doi.org/10.2495/CMEM150401 Bonamano S, Piermattei V, Madonia A, de Mendoza FP, Pierattini A, Martellucci R, Stefani C, Zappala G, Caruso G, Marcelli M (2016) The Civitavecchia Coastal Environment Moni- toring System (C-CEMS): a new tool to analyse the conflicts between coastal pressures and sensitivity areas. Ocean Science 12(1): 87-100. https://doi.org/10.5194/os-12-87-2016 Bonhomme D, Boudouresque CF, Astruch P, Bonhomme J, Bonhomme P, Goujard A, Thibaut T (2015) Typology of the reef formations of the Mediterranan seagrass Posidonia oceanica, 14 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) and the discovery of extensive reefs in the Gulf of Hyéres (Provence, Mediterranean). Sci- entific Report Port-Cros National Park 29: 41-73. Boudouresque CF, Jeudy de Grissac A, Olivier J (1984) First International workshop on Posi- donia oceanica Beds. GIS Posidonie, Marseille, France, 1-454. Boudouresque CF, Bernard G, Bonhomme P, Charbonnel EG, Diviacco E, Meinesz A, Pergent G, Pergent-Martini C, Ruitton S, Tunesi L (2006) Préservation et Conservation des Her- biers a Posidonia oceanica. RaMoGe pub., Monaco, 200 pp. Boudouresque CF, Bernard G, Bonhomme P, Charbonnel E, Le Diréach L, Ruitton S (2007) Monitoring methods for Posidonia oceanica seagrass meadows in Provence and the French Riviera. Scientific Report Port-Cros National Park 22: 17-38. Boudouresque CF, Bernard G, Pergent G, Shili A, Verlaque M (2009) Regression of Mediter- ranean seagrasses caused by natural processes and anthropogenic disturbances and stress: a critical review. Botanica Marina 52: 395-418. https://doi.org/10.1515/BOT.2009.057 Boudouresque CF, Bernard G, Bonhomme P, Charbonnel E, Diviacco G, Meinesz A, Pergent G, Pergent-Martini C, Ruitton S, Tunesi L (2012) Protection and conservation of Posido- nia oceanica meadows. RAMOGE and RAC/SPA publ., Tunis, 1-202. Buia MC, Gambi MC, Dappiano M (2004) I sistemi a fanerogame marine. In: Gambi MC, Dappiano M (Eds) Manuale di metodologie di campionamento e studio del benthos ma- rino Mediterraneo. Biologia Marina Mediterranea 11 (suppl. 1): 133-183. Ceccherelli G, Oliva S$, Pinna S, Piazzi L, Procaccini G, Marin-Guirao L, Dattolo E, Gallia R, La Manna G, Gennaro P, Costa MM, Barrote I, Silva J, Bulleri F (2018) Seagrass collapse due to synergistic stressors is not anticipated by phenological changes. Oecologia. 186(4): 1137-1152. https://doi.org/10.1007/s00442-018-4075-9 Cognetti de Martiis S (2016) Studio della produzione primaria di Posidonia oceanica L. Delile, 1813 nel Lazio settentrionale PhD thesis. La Tuscia University http://hdl.handle.net/2067/2965. Costanza R, d’Arge R, de Groot R, Farber S, Grasso M, Hannon B, Limburg K, Naeem S, O’Neill RV, Paruelo J, Raskin RG, Sutton P, van den Belt M (1997) The value of the world’s ecosys- tem services and natural capital. Nature 387: 253-26. https://doi.org/10.1038/387253a0 Cozza R, Rende F, Ferrari M, Bruno L, Pacenza M, Dattola L, Bitonti MB (2019) Biomoni- toring of Posidonia oceanica beds by a multiscale approach. Aquatic Botany 156: 14-24. https://doi.org/10.1016/j.aquabot.2019.04.002 Dellaporta SL, Wood J, Hicks JB (1983) A plant DNA mini preparation: version I]. RM.B.R. 1(4): 19-21. https://doi.org/10.1007/BF02712670 Den Hartog C (1970) The Seagrasses of the World. North Holland Publ. Co., Amsterdam. Dilipan E, Papenbrock J, Thangaradjou T (2017) Random Amplified Polymorphic DNA (RAPD) fingerprints evidencing high genetic variability among marine angiosperms of In- dia. Journal of the Marine Biological Association of the United Kingdom 97: 1307-1315. https://doi.org/10.1017/S0025315416000631 Di Maida G, Tomasello A, Sciandra M, Pirrotta M, Milazzo M, Calvo S (2013) Effect of different substrata on rhizomes growth, leaf biometry and shoot density of Posidonia oceanica. Marine Environment Research 87—88: 96-102. https://doi.org/10.1016/j.marenvres.2013.04.001 Domina G, Campisi P, Mannino AM, Sparacio I, Raimondo FM (2018) Environmental qual- ity assessment of the Sicilian coast using a multi-disciplinary approach. Acta Zoologica Bulgarica, Suppl. 11: 11-18. Posidonia oceanica (L) Delile 15 Duarte CM (2002) The future of seagrass meadows. Environmental Conservation 29: 192- 206. https://doi.org/10.1017/S0376892902000127 Echt GS, Erdahl LA, McCoy TJ (1992) Genetic segregation of random amplified polymorphic DNA in diploid cultivated alfaalfa. Genome 35: 84-87. https://doi.org/10.1139/g92-014 Excofher L, Lischer HEL (2010) Arlequin suite version 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10: 564-567. https://doi.org/10.1111/j.1755-0998.2010.02847.x Giraud G (1977) Contribution 4 la description et a la phénologie quantitative des herbiers 4 Posidonia oceanica (L.) Delile. Thése Doctorat 3°" Cycle, Univ, Aix-Marseille II, 150. https://doi.org/10.1515/botm.1977.20.8.487 Gnisci V (2014) Spatial and temporal analysis of health status of meadows along Civitavecchia coast (Latium, Italy). PhD thesis. La Tuscia University. http://hdl-handle.net/2067/2645 Hammer @, Harper DAT, Ryan PD (2001) PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4(1): 1-9. Hughes AR, Stachowicz JJ (2004) Genetic diversity enhances the resistance of a seagrass ecosys- tem to disturbance. Proceeding National Academy Science 101: 8998-9002. https://doi. org/10.1073/pnas.0402642101 Kim YK, Kim SH, Yi JM, Park SR, Lee KS (2019) Influence of environmental disturbances and reproductive strategy on genetic diversity and differentiation of Zostera marina populations on the southern coast of Korea. Marine Ecology 40: e12537. https://doi.org/10.1111/ maec. 12537 Jahnke M, Olsen JL, Procaccini G (2015) A meta-analysis reveals a positive correlation between genetic diversity metrics and environmental status in the long-lived seagrass Posidonia oce- anica. Molecular Ecology 24(10): 2336-2348. https://doi.org/10.1111/mec.13174 Jones TC, Gemmill CEC, Pilditch CA (2008) Genetic variability of New Zealand seagrass (Zostera muelleri) assessed at multiple spatial scales. Aquatic Botany 88: 39-46. https://doi. org/10.1016/j.aquabot.2007.08.017 Lopez y Royo C, Pergent G, Alcoverro T, Buia MC, Casazza G, Martinez-Crego B, Pérez M, Silvestre F, Romero J (2011) The seagrass Posidonia oceanica as indicator of coastal water quality: experimental inter calibration of classification systems. Ecological Indicators 11: 557-563. https://doi.org/10.1016/j.ecolind.2010.07.012 Macreadie PI, Schliep MT, Rasheed MA, Chartrand KM, Ralph PJ (2014) Molecular indica- tors of chronic seagrass stress: A new era in the management of seagrass ecosystems? Eco- logical Indicators 38: 279-281. https://doi.org/10.1016/j.ecolind.2013.11.017 Mannino AM (2018) Human activities trigger change in marine landscape. Flora Mediterranea 28: 295-311 https://doi.org/10.7320/FlMedit28.295 Marba N, Diaz-Almela E, Duarte CM (2014) Mediterranean seagrass (Posidonia oceani- ca) loss between 1842 and 2009. Biological Conservation 176: 183-190. https://doi. org/10.1016/j.biocon.2014.05.024 Martinez-Crego B, Vergés A, Alcoverro T, Romero J (2008) Selection of multiple seagrass in- dicators for environmental bio monitoring. Marine Ecology Progress Series 361: 93-109. https://doi.org/10.3354/meps07358 Meirmans PG (2006) Using the AMOVA framework to estimate a standardized genetic differ- entiation measure. Evolution 60: 2399-2402. https://doi.org/10.1554/05-631.1 16 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) Micheli C, Paganin P, Peirano A, Caye C, Meinesz A, Bianchi CN (2005) Genetic variability of Posidonia oceanica (L.) Delile in relation to local factors and biogeographic patterns. Aquatic Botany 82: 210-221. https://doi.org/10.1016/j.aquabot.2005.03.002 Micheli C, Spinosa F, Aliani S, Gasparini GP, Molcard A, Peirano A (2010b) Genetic input by Posidonia oceanica (L) Delile fruits dispersed by currents in the Ligurian Sea. Plant Biosys- tems 144: 333-339. https://doi.org/10.1080/11263501003764798 Micheli C, Cupido R, Lombardi C, Belmonte A, Peirano A (2012) Changes in genetic struc- ture of Posidonia oceanica at Monterosso al Mare (Ligurian Sea) and its resilience over a decade (1998-2009). Environmental Management 50: 598-606. https://doi.org/10.1007/ s00267-012-9917-3 Micheli C, D’Esposito D, Belmonte A, Peirano A, Valiante LM, Procaccini G (2015) Genetic diversity and structure in two protected Posidonia oceanica meadows. Marine Environmen- tal Research 109: 124-131. https://doi.org/10.1016/j.marenvres.2015.06.016 Molinier R, Picard J (1952) Recherches sur les herbiers de phanérogames marines du littoral Méditerranéen francais. Annales Institut Océanographique 27(3): 157-234. Montefalcone M (2009) Ecosystem health assessment using the seagrass Posidonia oce- anica: a review. Ecological Indicators 9: 595-604. https://doi.org/10.1016/j.ec- olind.2008.09.013 Nei M, Li WH (1979) Mathematical model for studying genetic variation in terms of restric- tion endonucleases. Proceedings of the National Academy of Sciences, USA 76: 5269-— 5273. https://doi.org/10.1073/pnas.76.10.5269 Orth RJ, Carruthers TJB, Dennison WC, Duarte CM, Fourqurean JW, Heck KL, Hughes AR, Kendrick GA, Kenworthy WJ, Olyarnik S, Short F Waycott M, Williams SL (2006) A global crisis for seagrass ecosystems. BioScience 56(12): 987-996. https://doi. org/10.1641/0006-3568(2006)56[987:AGCFSE]2.0.CO;2 Panayotidis P, Boudouresque CE, Marcot-Coqueugniol J (1981) Microstructure de lherbier de Posidonia oceanica (L) Delile. Botanica Marina 24: 115-124. https://doi.org/10.1515/ botm.1981.24.3.115 Pergent G, Pergent-Martini C, Boudouresque CF (1995) Utilisation de Pherbier 4 Posidonia oceanica comme indicateur biologique de la qualité du milieu littoral en Méditerranée: état de connaissances. Mésogée 54: 3-27. Pergent-Martini C, Leoni V, Pasqualini V, Ardizzone GD, Balestri E, Bedini R, Belluscio A, Belsher T, Borg J, Boudouresque CF, Boumaza S, Bouquegneau JM, Buia MC, Calvo S, Cebrian J, Charbonnel E, Cinelli EK Cossu A, Di Maida G, Dural B, Francour P, Gobert S, Lepoint G, Meinesz A, Molenaar H, Mansour HM, Panayotidis P, Peirano A, Pergent G, Piazzi L, Pirrotta M, Relini G, Romero J, Sanchez-Lizaso JL, Semroud R, Shembri P, Shili A, Tomasello A, Velimirov B (2005) Descriptors of Posidonia oceanica meadows: use and application. Ecological Indicators 5: 213-230. https://doi-org/10.1016/j.ec- olind.2005.02.004 Personnic S$, Boudouresque CF, Astruch P, Ballesteros E, Blouet S, Bellan-Santini D, Bon- homme P, Thibault-Botha D, Feunteun E, Harmelin-Vivien M, Pergent G, Pergent-Mar- tini C, Pastor J, Poggiale JC, Renaud F, Thibaut T, Ruitton S (2014) An ecosystem-based Posidonia oceanica (L) Delile sles approach to assess the status of a Mediterranean ecosystem, the Posidonia oceanica seagrass meadow. PLOS One 9(6): 1-17 (e98994). https://doi.org/10.137 1/journal.pone.0098994 Procaccini G, Piazzi L (2001) Genetic polymorphism and transplantation success in the Mediterranean seagrass Posidonia oceanica. Restoration Ecology 9: 332-338. https://doi. org/10.1046/j.1526-100x.2001.009003332.x Procaccini G, Olsen JL, Reusch TBH (2007) Contribution of genetics and genomics to sea- grass biology and conservation. Journal of Experimental Marine Biology and Ecology 350: 234-359. https://doi.org/10.1016/j.jembe.2007.05.035 R Development Core Team (2008) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org Reusch TBH, Ehlers A, Hammerli A, Worm B (2005) Ecosystem recovery after climatic extremes enhanced by genotypic diversity. Proceedings of the National Academy of Sciences of the United States of Americal102(8): 2826-2831. https://doi.org/10.1073/pnas.0500008102 Ribeiro AM, Foote AD, Kupczok A, Frazao B, Limborg MT, Pifeiro R, Abalde S, Rocha S, da Fonseca RR (2016) Marine genomics: News and views. Marine Genomics 31: 1-8. https:// doi.org/10.1016/j.margen.2016.09.002 Rohlf FJ (1993) NTSYS- PC Numerical Taxonomy and Multivariate Analysis System. Exeter Software. Setauket, NY. www.exetersoftware.com/downloads/ntsysguide2 1 Rotini A, Micheli C, Valiante LM, Migliore L (2011) Assessment of Posidonia oceanica (L.) Delile conservation status by standard and putative approaches: the case study of Santa Marinella meadow (Italy, W Mediterranean). Open Journal of Ecology 1: 48-56. https:// doi.org/10.4236/oje.2011.12006 Rotini A, Belmonte A, Barrote I, Micheli C, Peirano A, Santos RO, Silva J, Migliore L (2013) Effectiveness and consistency of a suite of descriptors for assessing the ecological status of seagrass meadows (Posidonia oceanica L. Delile). Estuarine, Coastal and Shelf Science 130: 252-259. https://doi.org/10.1016/j.ecss.2013.06.015 Semroud R (1993) Contribution 4 la connaissance de l’ecosysteme a Posidonia oceanica (L.) Delile dans la région d’Alger (Algérie): étude de quelques compartiments. PhD thesis. Serra IA, Innocenti AM, Di Maida G, Calvo S, Migliaccio M, Zambianchi E, Pizzigalli C, Arnaud-Haond S, Duarte CM, Serrao A, Procaccini G (2015) A meta-analysis reveals a positive correlation between genetic diversity metrics and environmental status in the long-lived seagrass Posidonia oceanica. Molecular Ecology 24(10): 2336-2348. https://doi. org/10.1111/mec.13174 Spalding M, Taylor M, Ravilious C, Short K Green E (2003) The distribution and status of seagrasses. In: Green EP, Short FT (Eds) World Atlas of Seagrass. University California Press, Berkeley CA, 5-26. Vassallo P, Paoli C, Rovere A, Montefalcone M, Morri C, Bianchi CN (2013) The value of the seagrass Posidonia oceanica: a natural capital assessment. Marine Pollution Bulletin 75: 157-167. https://doi.org/10.1016/j.marpolbul.2013.07.044 Williams JGK, Kubelik AR, Livak J, Antoni Rafalski J, Tingey SV (1990) DNA polymor- phisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acid Research 18: 6531-6535. https://doi.org/10.1093/nar/18.22.6531 18 Valentina Gnisci et al. / Italian Botanist 9: 1-19 (2020) Supplementary material | Figure S1. Heterogeneous spatial distribution of seagrass population along the northern coast of Lazio (Italy), Central Tyrrhenian, Mediterranean Authors: Valentina Gnisci, Selvaggia Cognetti de Martiis, Alessandro Belmonte, Carla Micheli, Viviana Piermattei, Simone Bonamano, Marco Marcelli Data type: occurrence Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODDbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/italianbotanist.9.46426.suppl1 Supplementary material 2 Figure S2. The non-parametric Kruskal-Wallis test Authors: Valentina Gnisci, Selvaggia Cognetti de Martiis, Alessandro Belmonte, Carla Micheli, Viviana Piermattei, Simone Bonamano, Marco Marcelli Data type: statistical data Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODDbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: hettps://doi.org/10.3897/italianbotanist.9.46426.suppl2 Supplementary material 3 Table $1. ANOVA results Authors: Valentina Gnisci, Selvaggia Cognetti de Martiis, Alessandro Belmonte, Carla Micheli, Viviana Piermattei, Simone Bonamano, Marco Marcelli Data type: statistical data Explanation note: Multifactor three-way ANOVA results to test the presence of sig- nificant differences in shoot density, leaf number, leaf length, leaf surface and LAI among the samples collected at different depth, location and substrate type. Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/italianbotanist.9.46426.suppl3 Posidonia oceanica (L) Delile 19 Supplementary material 4 Table S2. AMOVA Authors: Valentina Gnisci, Selvaggia Cognetti de Martiis, Alessandro Belmonte, Carla Micheli, Viviana Piermattei, Simone Bonamano, Marco Marcelli Data type: statistical data Copyright notice: This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODDbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited. Link: https://doi.org/10.3897/italianbotanist.9.46426.suppl4