
==== Front
Sci Data
Sci Data
Scientific Data
2052-4463
Nature Publishing Group UK London

39256479
3787
10.1038/s41597-024-03787-y
Data Descriptor
The NEREA Augmented Observatory: an integrative approach to marine coastal ecology
http://orcid.org/0000-0002-3543-738X
Campese Lucia lucia.campese@szn.it

1
http://orcid.org/0000-0002-8550-6400
Russo Luca 1
http://orcid.org/0000-0002-8571-4816
Abagnale Maria 2
http://orcid.org/0000-0003-3372-9423
Alberti Adriana 319
http://orcid.org/0000-0002-8077-1512
Bachi Giancarlo 4
http://orcid.org/0000-0003-4433-4793
Balestra Cecilia 5
http://orcid.org/0000-0002-3319-7036
Bellardini Daniele 16
Buondonno Angela 2
http://orcid.org/0000-0002-0816-6158
Cardini Ulisse 7
http://orcid.org/0000-0002-7188-9279
Carotenuto Ylenia 1
http://orcid.org/0000-0002-7960-7168
Checcucci Giovanni 4
http://orcid.org/0000-0002-6296-7132
Chiusano Maria Luisa 28
http://orcid.org/0000-0001-5486-9535
D’Ambra Isabella 19
http://orcid.org/0000-0002-6821-6962
d’Ippolito Giuliana 10
http://orcid.org/0000-0003-2959-8977
Di Capua Iole 19
http://orcid.org/0000-0002-2372-7072
Donnarumma Vincenzo 111
http://orcid.org/0000-0002-5453-461X
Fontana Angelo 1012
Furia Marta 2
Galarza-Verkovitch Denisse 1
Gallia Roberto 13
http://orcid.org/0000-0001-7467-8509
Labadie Karine 14
http://orcid.org/0000-0003-4787-0882
Leone Serena 15
http://orcid.org/0000-0001-6963-0068
Licandro Priscilla 1
Longo Antonio 116
http://orcid.org/0000-0001-5729-9574
Maselli Maira 2
http://orcid.org/0000-0002-7032-1403
Merquiol Louise 117
http://orcid.org/0000-0003-4180-0410
Murano Carola 19
http://orcid.org/0000-0003-3161-8367
Oliveira Pedro H. 14
Passarelli Augusto 13
Percopo Isabella 2
Perdereau Aude 14
http://orcid.org/0000-0003-4426-991X
Piredda Roberta 1
http://orcid.org/0000-0002-7462-5922
Raffini Francesca 115
http://orcid.org/0000-0003-3024-7568
Roncalli Vittoria 1
http://orcid.org/0000-0001-7473-6086
Ruscheweyh Hans-Joachim 18
http://orcid.org/0000-0002-7329-3559
Russo Ennio 110
http://orcid.org/0000-0001-5794-0050
Saggiomo Maria 2
Santinelli Chiara 49
http://orcid.org/0000-0001-9697-5301
Sarno Diana 29
http://orcid.org/0000-0003-3065-0314
Sunagawa Shinichi 18
Tramontano Ferdinando 1
Trano Anna Chiara 1
http://orcid.org/0000-0001-7026-0156
Uttieri Marco 19
http://orcid.org/0000-0001-7562-3454
Wincker Patrick 14
http://orcid.org/0000-0003-3028-5812
Zampicinini Gianpaolo 15
http://orcid.org/0000-0002-9876-4601
Casotti Raffaella raffaella.casotti@szn.it

1
http://orcid.org/0000-0003-2037-2277
Conversano Fabio fabio.conversano@szn.it

2
http://orcid.org/0000-0002-2189-503X
D’Alelio Domenico domenico.dalelio@szn.it

19
http://orcid.org/0000-0002-7473-394X
Iudicone Daniele daniele.iudicone@szn.it

1
http://orcid.org/0000-0003-0757-5934
Margiotta Francesca francesca.margiotta@szn.it

2
http://orcid.org/0000-0002-2475-1787
Montresor Marina marina.montresor@szn.it

1
1 https://ror.org/03v5jj203 grid.6401.3 0000 0004 1758 0806 Department of Integrative Marine Ecology, Stazione Zoologica Anton Dohrn, Villa Comunale, 80121 Naples, Italy
2 https://ror.org/03v5jj203 grid.6401.3 0000 0004 1758 0806 Department of Research Infrastructures for marine biological resources, Stazione Zoologica Anton Dohrn, Villa Comunale, 80121 Naples, Italy
3 https://ror.org/03xjwb503 grid.460789.4 0000 0004 4910 6535 Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), 91198 Gif-sur-Yvette, France
4 grid.419463.d 0000 0004 1756 3731 Biophysics Institute, CNR, Via G. Moruzzi 1, 56124 Pisa, Italy
5 https://ror.org/04y4t7k95 grid.4336.2 0000 0001 2237 3826 National Institute of Oceanography and Applied Geophysics, OGS, via Auguste Piccard 54, 34151 Trieste, Italy
6 https://ror.org/0107c5v14 grid.5606.5 0000 0001 2151 3065 DiSTAV, Department of Earth, Environment and Life Sciences, University of Genoa, Corso Europa 26, 16132 Genoa, Italy
7 https://ror.org/03v5jj203 grid.6401.3 0000 0004 1758 0806 Department of Integrative Marine Ecology, Stazione Zoologica Anton Dohrn, Genoa Marine Centre, Villa del Principe, Piazza del Principe 4, 16126 Genoa, Italy
8 https://ror.org/05290cv24 grid.4691.a 0000 0001 0790 385X Department of Agricultural Sciences, University Federico II of Naples, Piazza Carlo di Borbone, 1, 80055 Portici, Napoli Italy
9 NBFC, National Biodiversity Future Center, Piazza Marina 61, 90133 Palermo, Italy
10 https://ror.org/04zaypm56 grid.5326.2 0000 0001 1940 4177 Bio-Organic Chemistry Unit, Institute of Bio-molecular Chemistry, National Research Council, Via Campi Flegrei 34, 80078 Pozzuoli, Italy
11 https://ror.org/003109y17 grid.7763.5 0000 0004 1755 3242 Department of Life and Environmental Sciences, University of Cagliari, Via Tommaso Fiorelli 1, 09126 Cagliari, Italy
12 https://ror.org/05290cv24 grid.4691.a 0000 0001 0790 385X Laboratory of Bio-Organic Chemistry and Chemical Biology, Biology department, University of Naples “Federico II”, Via Cupa Nuova Cinthia 21, 80126 Napoli, Italy
13 https://ror.org/03v5jj203 grid.6401.3 0000 0004 1758 0806 Department of Research Infrastructures for marine biological resources, Stazione Zoologica Anton Dohrn, Via Nuovo Macello 16, 80055 Portici Naples, Italy
14 grid.434728.e 0000 0004 0641 2997 Génomique Métabolique, Genoscope, Institut François Jacob, Commissariat à l’Energie Atomique (CEA), CNRS, Université Evry, Université Paris-Saclay, 2 Rue Gaston Crémieux, 91057 Evry, France
15 https://ror.org/03v5jj203 grid.6401.3 0000 0004 1758 0806 Department of Biology and Evolution of Marine Organisms, Stazione Zoologica Anton Dohrn, Villa Comunale, 80121 Naples, Italy
16 https://ror.org/05290cv24 grid.4691.a 0000 0001 0790 385X Department of Biology, University of Napoli “Federico II”, Via Cupa Nuova Cinthia 21, 80126 Napoli, Italy
17 https://ror.org/035xkbk20 grid.5399.6 0000 0001 2176 4817 TELEMMe - MMSH, Aix-Marseille Université, 5 rue du Château de l’Horloge, 13097 Aix-en-Provence, France
18 https://ror.org/05a28rw58 grid.5801.c 0000 0001 2156 2780 Department of Biology, Institute of Microbiology and Swiss Institute of Bioinformatics, ETH Zurich, Zurich, 8093 Switzerland
19 grid.434728.e 0000 0004 0641 2997 Present Address: Génomique Métabolique, Genoscope, Institut François Jacob, Commissariat à l’Energie Atomique (CEA), CNRS, Université Evry, Université Paris-Saclay, 2 Rue Gaston Crémieux, 91057 Evry, France
10 9 2024
10 9 2024
2024
11 98926 4 2024
13 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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-nc-nd/4.0/.
The NEREA (Naples Ecological REsearch for Augmented observatories) initiative aims to establish an augmented observatory in the Gulf of Naples (GoN), designed to advance the understanding of marine ecosystems through a holistic approach. Inspired by the Tara Oceans expedition and building on the scientific legacy of the MareChiara Long-Term Ecological Research (LTER-MC) site, NEREA integrates traditional physical, chemical, and biological measurements with state-of-the-art methodologies such as metabarcoding and metagenomics. Here we present the first 10 months of NEREA data, collected from April 2019 to January 2020, encompassing physico-chemical parameters, plankton biodiversity (e.g., microscopy and flow cytometry), prokaryotic and eukaryotic metabarcoding, a prokaryotic gene catalogue, and a collection of 3818 prokaryotic Metagenome-Assembled Genomes (MAGs). NEREA’s efforts produce a significant volume of multifaceted data, which enhances our understanding of marine ecosystems and promotes the development of scientific hypotheses and ideas.

Subject terms

Ecology
Ecology
TechOceanS (101000858) AtlantECO (862923) MARCO-BOLO (101082021)BIOcean5D (101059915)https://doi.org/10.13039/501100010685 Stazione Zoologica Anton Dohrn (SZN) https://doi.org/10.13039/100008509 Open University (OU) https://doi.org/10.13039/501100001711 Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (Swiss National Science Foundation) OBAMA-NEXT (101081642) TechOceanS (101000858) EuroSea (862626)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcBackground & Summary

Oceans harbour invaluable biodiversity that plays a pivotal role in driving biogeochemical cycles and enhancing ecosystem productivity. Coastal ecosystems provide humankind with a wide array of goods and services, including food, biochemical compounds, fossil fuels, and recreational spaces1.

Coastal systems face threats from a variety of intensive human activities including agriculture, fisheries, aquaculture, shipping, urbanisation, and tourism. Addressing such increasing threats requires integrated approaches within well-established monitoring frameworks2, such as the Long-Term Ecological Research (LTER) sites, crucial for identifying significant ecological events under a changing global climate, while providing valuable insights for society management and policymakers3,4. Among LTERs, the MareChiara LTER site (LTER-MC) in the Gulf of Naples (GoN; Central Tyrrhenian Sea5,6, Western Mediterranean Sea), established in 1984, represents one of the world’s longest ecological time series focused on the study of plankton, as well as on the physico-chemical state of the water column7.

The GoN, spanning roughly 195 km of coastline, faces significant pressures from a variety of human activities. The city of Naples and its metropolitan area, with about 4 million inhabitants, see an influx of tourists during summer, further impacting both the mainland and the highly popular islands of Capri, Ischia and Procida. The region is also a hub for industrial and commercial activities, including one of the largest ports in the Mediterranean Sea, which amplifies the anthropogenic impact. Associated anthropogenic stresses, such as illegal dumping and discharges, intense maritime traffic and industrial activities are only examples of the continuous human pressures on this area8. It is within this context that the LTER-MC was established and continues to operate.

Implementing existing observation sites with data generated through innovative and diversified approaches is crucial for the development of “augmented observatories”, as stated by the G7 Turin declaration9, which encourages multidisciplinary marine biodiversity observations. This effort is essential to pursue the ‘Ocean Health Index’ (OHI) assessment10 which aims to employ scientific methods to quantitatively characterise specific biological, physical, economic, and social aspects and to guide policymakers towards the sustainable management of ocean resources. In line with this, a multidisciplinary and integrative investigation of plankton, ranging from virus to jellyfish, was carried out by the Tara Oceans project, an oceanographic expedition aimed at using cutting-edge molecular, microscopical and optical techniques to acquire multifaceted data on the wide variety of species contributing to marine plankton11. Plankton encompasses autotrophic, heterotrophic, and mixotrophic species from all kingdoms of life12, spanning six orders of magnitude in size13. These organisms play a pivotal role in driving the Earth’s biogeochemical cycles14 and account for half of the planet’s oxygen production15.

Based upon the conceptual and methodological framework developed by the Tara Oceans project and building up on the scientific legacy of the LTER-MC site, we here present NEREA (Naples Ecological REsearch for Augmented observatories)16. This initiative aims to establish an augmented observatory in the GoN to study the dynamics of marine plankton communities. The project integrates expertise from various fields, including physics, oceanography, chemistry, physiology, molecular biology, bioinformatics, data science, modelling, and theoretical ecology. This integration facilitates a synergistic approach where laboratory experiments and computational analyses complement field observations. Such observations are refined by experimental and theoretical insights, creating a feedback loop that enhances ecological research methods. The collaborative nature of NEREA produces a comprehensive dataset that advances our understanding of marine ecosystems and addresses ecological questions.

Additionally, NEREA is committed to fostering international collaborations and it is actively involved in the global-scale research program BioGeoSCAPES, aimed to improve our understanding of the microbial biogeochemistry of the oceans. NEREA is also part of the UN Decade of Ocean Sciences (UNDOS) under the Ocean Biomolecular Observing Network (OBON) Program.

Here we present the first 10 months of NEREA data, collected from April 2019 to January 2020, encompassing physico-chemical parameters, plankton biodiversity (microscopy and flow cytometry), prokaryotic and eukaryotic metabarcoding, a prokaryotic gene catalogue, and a collection of 3818 prokaryotic Metagenome-Assembled Genomes (MAGs).

Methods

Study area

The GoN is one of the most representative Mediterranean coastal basins, with key interplay between physical drivers and ecosystem responses (e.g., 17–20). It is characterised by five persistent water masses17: Coastal Surface Water (CSW), Tyrrhenian Surface Water (TSW), Tyrrhenian Intermediate Water (TIW), Atlantic Water (formerly referred to as Modified Atlantic Water), and Levantine Intermediate Water (LIW)6. Local climatic conditions, including wind patterns and river discharges, significantly influence these waters, particularly near the shorelines (e.g., 18,19,21). NEREA covers the whole GoN area with monthly sampling at the LTER-MC site (Lat 40°48.5′ N, Long 14° 15′ E; bottom depth: ca. 75 m), located two nautical miles off downtown Naples; this site experiences alternating influences from nutrient-rich coastal waters and the nutrient-poor waters of the Mid Tyrrhenian Sea7. The LTER-MC site (DEIMS iD: https://deims.org/0b87459a-da3c-45af-a3e1-cb1508519411) belongs to the Long-Term Ecological Research national and international networks (LTER-Italy, LTER-Europe and ILTER).

In addition to sampling at LTER-MC site, NEREA includes seasonal sampling at the River Sarno mouth and plume (Lat 40°43′ N, Long 14°27′ E; bottom depth: ca. 15 m) and a biannual sampling at a deep station located at the Dohrn Canyon (Lat 40°36′ N, Long 14°08′ E; bottom depth ca. 1300 m)22 near the Island of Capri (Fig. 1). The Sarno river is highly polluted23, bringing intermittent loads of nutrients and chemical pollutants to the coast from intensive agricultural activities and leather factories24. On the other hand, the Dohrn Canyon, extending 25 km and located around 12 nautical miles off Naples, is a deep bifurcated submarine structure. It represents an active geological site contributing to the upwelling of deep waters and generating turbidity currents at the seafloor25. While primarily investigated for its geological attributes25, the canyon biodiversity remains largely unknown26. It represents a biodiversity hotspot, offering critical ecological benefits and serving as a breeding ground for various marine species along the coastline27.Fig. 1 Geographical map of the Gulf of Naples (GoN) with NEREA sampling sites. LTER-MC: Long-Term Ecological Research MareChiara; NRC: NEREA-Capri (Dohrn Canyon); NRS: NEREA-Sarno (Sarno river).

NEREA operational units

NEREA’s multidisciplinary approach leads to the production of a wide variety of data (Fig. 2), which can be listed as follows:(i) Physical parameters;

(ii) Chemical data, including nutrients, particulate and dissolved organic matter, particulate lipid-derived oxylipins and stable isotope ratios;

(iii) Plankton biodiversity data, including picoplankton quantification by flow cytometry; species identification and quantification by light microscopy of phyto- micro- and mesozooplankton;

(iv) Plankton metabarcoding and metagenomic data.

Fig. 2 Presence and absence matrix of analysed parameters. The matrix displays the data availability for each parameter across NEREA samples. Torquoise squares indicate the presence of data, while white squares represent the absence of data for the corresponding parameter in each sample. DOC: Dissolved Organic Carbon; POC: Particulate Organic Carbon; LOFAs: Linear Oxygenated Fatty Acids; FCM: Flow Cytometry; metaB: metabarcoding; MetaG: metagenomics.

All scientists operate in a coordinated way both in laboratory and during onboard activities, where simultaneous sampling and observations are run (Fig. 3).Fig. 3 NEREA sampling strategy. A Conductivity, Temperature, and Depth (CTD) probe attached to a rosette equipped with multiple Niskin bottles collects distinct water samples. Plankton samples are also collected by deploying nets of appropriate mesh size (20 µm for phytoplankton and 200 µm for zooplankton). Part of the collected water is subjected to size fractionation, involving a series of sieves and filters with decreasing mesh and pore sizes, for the separation of plankton into distinct size categories. Samples obtained from Niskin bottles and nets are subsequently chemically analysed and processed via flow cytometry, microscopy, and DNA/RNA extraction followed by sequencing.

Physical state of the system

This task is firstly pursued onboard by the real-time production and observation of vertical profiles for multiple parameters, such as temperature (°C), salinity, density (σ), dissolved oxygen (% and mL L−1), fluorescence (µg L−1), transmittance (%), turbidity (NTU), pH and photosynthetically active radiation (PAR). The above-mentioned observations are obtained using a SeaBird Electronic SBE 911 plus v2 multi-parameter profiler. This probe acquires 24 data variables per second with an accuracy of 0.002 °C for temperature and 0.001 S/m for conductivity, and it is equipped on a rosette-type sampler with twelve Niskin bottles of 10 L each, which are employed in the collection of water samples at different depths, based on the physical structure of the water column.

The procedures described are essential for analysing the water column characteristics based on its vertical physical structure, enabling real-time planning for further water sampling. In the GoN, particularly near the coast, the water column tends to be mixed and homogeneous during winter but becomes highly stratified from late spring to autumn28. This seasonal variability leads to different environmental conditions, which in turn influence the composition of planktonic communities at various depths. Consequently, sampling strategies need to be adjusted seasonally. For example, in the context of meta-omic studies, at the coastal LTER-MC site only surface water is collected during winter, while during periods of stratification an additional sample is taken at 40 m depth to account for the different conditions and planktonic communities. In addition, a thermosalinograph is used during the transit between sampling sites, collecting continuous data for the superficial temperature, conductivity, and salinity.

Chemical state of the system

The seawater used for the chemical analyses is collected from Niskin bottles at different depths. The analyses involve the measurement of nutrients, total nitrogen (TN), total phosphorus (TP), dissolved organic carbon (DOC), particulate organic carbon (POC), particulate nitrogen (PN), particulate lipid-derived oxylipins (Linear Oxygenated Fatty Acids, LOFAs), chlorophyll a and other photosynthetic pigments.

At the LTER-MC site, nutrients, TN, TP, and DOC are collected at 0, 2, 5, 10, 20, 30, 40, 50, 60, and 70 m depths; POC and HPLC at 0, 10, 20, 40, and 60 m depths; LOFAs at 0 m depth, chlorophyll a at 0, 2, 5, 10, 20, 40, and 60 m depths. At the Dohrn Canyon station, all the chemical analyses are carried out on samples collected at 0, 10, 20, 40, 60, 80, 100 m depths. Nutrients, TN, TP, DOC and POC are analysed also at the deeper depths of 150, 200, 250, 300, 400, 600, 800 m. LOFAs are analysed at 0 m depth. At Sarno, all the chemical analyses are carried out on samples collected at 0 m.

Nutrients, TN, and TP

Determination of inorganic nutrients (NO3−, NO2−; PO4−, SiO4−) is performed on samples directly collected from the Niskin bottles into high-density polypropylene vials and frozen at −20 °C. Analyses are performed using a Flow-Sys III Systea auto-analyzer using the methods reported by29. For the determination of TN and TP, samples are taken directly from the Niskin bottle in 100 ml HDPE containers and immediately frozen at −20 °C. Analyses are performed after digestion at 120 °C for 30 minutes, following the method proposed by30.

DOC and POC

Samples for DOC and POC analysis are collected directly from the Niskin bottles into acid-washed Nalgene polycarbonate bottles. Samples for DOC are stored at −20 °C until filtration and analysis. Filtration is carried out using 0.22 µm polyethersulfone (PES) filters (Minisart®), DOC analysis is carried out by high-temperature catalytic oxidation with a Shimadzu TOC analyzer (TOC-Vcsn) following the methods described in31. For POC determination, filtrations are conducted on board using a variable volume (1–3 L) of seawater filtered on pre-combusted 25 mm GF/F filters. After filtration, filters are rinsed with deionized water and stored at −20 °C. The analyses are carried out using a Thermo Scientific FlashEA 1112 elemental analyzer (Thermo Fisher Scientific) following32 and using cyclohexanone-2,4-dinitrophenylhydrazone as standard.

LOFAs

For the analysis of particulate oxylipins a variable volume (0.75–3 L) of seawater is collected using Niskin bottles, pre-filtered onto 200 μm mesh nylon net, and then filtered through polycarbonate filters (47 mm diameter, 2 μm mesh size); filters are kept in 2 ml Eppendorf tubes and promptly frozen in liquid nitrogen. After thawing, filters are sonicated in Milli-Q water, left at room temperature for 30 min and then oxylipins extracted and quantified through Liquid Chromatography-Mass Spectrometry (LC-MS) after addition of a known amount of 16-hydroxyhexadecanoic acid as internal standard. According to the targeted-metabolomic approach used, only LOFAs deriving from hexadecatrienoic (C16), eicosapentaenoic (C20), and docosahexaenoic (C22) fatty acid precursors are quantified. Amount of LOFAs is normalised by the volume of the seawater filtered and reported as ng-LOFAs/L33.

Total chlorophyll a

The seawater for the determination of the total chlorophyll a is filtered (200 to 1080 ml) on GF/F (25 mm diameter), which are quickly frozen in liquid nitrogen. Once thawed, chlorophyll a is extracted in acetone 90% and its concentration determined according to34 using a spectrofluorometer (Shimadzu RF-5301 PC), which is daily calibrated using a chlorophyll a standard solution (from Anacystis nidulans; Sigma).

Photosynthetic pigments

For the analysis of the pigment spectrum using High Performance Liquid Chromatography (HPLC), 2-3 L of seawater are collected at each depth and filtered through GF/F filters (47 mm in diameter), which are then stored in liquid nitrogen. The pigment separations are conducted using an Agilent 1100 HPLC (Agilent technologies, United States) according to the method outlined in35 and modified by36. The HPLC system, equipped with an HP 1050 photodiode array detector and a HP 1046 A fluorescence detector, is specifically used for the measurement of chlorophyll degradation products. Calibration of the instrument is performed using external standard pigments supplied by the International Agency for 14C determination-VKI Water Quality Institute.

Carbon and nitrogen stable isotope ratios of plankton

Water samples for pico- and nanoplankton (<20 µm) are collected at a depth of 1 m using a Niskin bottle and then transferred into plastic bins. These samples are transported to the laboratory on ice, where they are pre-filtered through a 20 µm mesh. The filtered water is subsequently concentrated onto pre-combusted (4 hours at 500 °C) GF/F filters. These filters are individually frozen and later freeze-dried. Water samples for microplankton (20–200 µm) and mesozooplankton (200–2,000 µm) are collected through vertical tows from approximately 70 m depth to the surface using 20 µm and 200 µm plankton nets, respectively. Samples are transferred into plastic jars and transported to the laboratory on ice. In the lab, microplankton samples are concentrated onto 20 µm filters, and mesozooplankton samples onto 200 µm filters. These are then transferred into vials, frozen, and freeze-dried. For the determination of carbon and nitrogen stable isotopes, 1.0 ± 0.1 mg of each sample is sent to the Stable Isotope Facility at the University of California, Davis, USA. Analysis is performed using a PDZ Europa ANCA-GSL elemental analyser interfaced with a PDZ Europa 20–20 isotope ratio mass spectrometer (Sercon Ltd., Cheshire, UK). The instrumental precision is ±0.02‰ for δ13C and ± 0.03‰ for δ15N. Final δ13C and δ15N values are expressed relative to the international standards Vienna Pee Dee Belemnite (VPDB) for carbon and nitrogen. To ensure reproducibility, duplicate aliquots of 26 samples are analysed, yielding a reproducibility of ± 0.04‰ for δ13C and ± 0.06‰ for δ15N.

Plankton biodiversity - Flow Cytometry

Discrete samples for picoplankton are collected from the Niskin bottles and kept in the dark at 4 °C until analysis (less than 2 h). Subsamples of 1 ml are either run unstained for the picophytoplankton counts or fixed with a mix of 1% paraformaldehyde and 0.5% glutaraldehyde, frozen in liquid nitrogen and stored at −80 °C37. Once thawed, samples are stained with SYBRGreen I for heterotrophic bacteria38 and processed using a Becton-Dickinson FACVerse flow cytometer (BD BioSciences) equipped with a 488 nm solid-state laser and a standard filter set. Picophytoplankton (the phototrophic part) and bacteria are classified based on their scatter versus fluorescence properties and counted by the volumetric device of the instrument, after appropriate adjustments due to dilution. Gating is on red fluorescence for the picophytoplankton and green fluorescence from the SYBRGreen stain for the heterotrophic prokaryotes. One µm beads (PolySciences) are used as internal standards and for normalisation of scatter and fluorescence values (arbitrary units, a.u.). Files generated are elaborated using the FCS Express software (DeNovo Software Inc).

Plankton biodiversity - Light Microscopy

Identification and quantification by light microscopy of unicellular eukaryotic organisms, i.e., protists, is carried out following two different analytical protocols, depending on the size and abundance of the organisms that roughly correspond to phytoplankton and microzooplankton. We grouped under the “Phytoplankton” label mainly autotrophic unicellular species but also small-sized mixotrophic and heterotrophic taxa. “Microzooplankton” includes ciliates and large heterotrophic and mixotrophic dinoflagellates. Phytoplankton and microzooplankton samples are collected from the same Niskin bottles used for the meta-omic analyses. 250 ml of water in a dark glass bottle in duplicate is fixed with Lugol’s solution (1%) and stored in the dark. Phytoplankton and microzooplankton analysis are performed according to the Utermöhl method after the settling of a variable volume of sample39. Settled volumes vary in relation to the concentration of target cells ranging from 3 to 100 ml for phytoplankton and from 100 to 250 ml for microzooplankton. Counting is done over varying fractions of the sedimentation chamber, generally 2–4 transects for phytoplankton, half, or entire chamber for microzooplankton. Identification and counting are done using an inverted light microscope Zeiss Axiovert 200 (Carl Zeiss, Germany) at 400x (phytoplankton) and 200x (microzooplankton) magnifications. Specimens are identified at the lowest possible taxonomic level based on classic taxonomy books40–45 supplemented by more recent specialised papers. Lists of taxa have been updated and checked for nomenclature and synonyms as reported in the taxonomic reference websites (Algaebase http://www.algaebase.org/ and WoRMS http://www.marinespecies.org/).

In the laboratory, mesozooplankton samples are analysed after being fixed in 95% ethanol. Initially, the sample is concentrated and re-suspended in 200 ml of ethanol. Depending on the sample density, aliquots are taken using the Huntsman beaker subsampling method at ratios ranging from 1/4 to 1/3246. These aliquots are thoroughly mixed and then subsampled using a large bore graduated pipette. The aliquots are transferred into a 10 ml mini-Bogorov chamber for examination under a Leica M165C stereomicroscope. In most cases, copepods are identified at the highest taxonomic levels, i.e. species, genus and family. Among other mesozooplankton, some groups are identified at the species level (e.g., Cladocerans, Chaetognaths, Siphonophores), while other holoplankton (e.g., Amphipods, Mysids, Ostracods, Euphausiids, Pteropods, Hydromedusae) and meroplanktonic taxa (e.g., Decapods, Cirripeds, Echinoderms, Bivalves, Molluscs, Gastropods) are identified at higher taxonomic levels. All the mesozooplanktonic species abundances are reported as individuals m−3.

Plankton metabarcoding and metagenomics – sampling

The sampling strategy for the biological samples used herein, from the sampling devices to the size fraction division of organisms, is based on that used in the Tara Oceans project47,48, integrating and adapting those procedures to the coastal system studied in NEREA. The size fractions used are efficient for effectively separating the principal components of the marine plankton: 0.2–3 µm (picoplankton: e.g., archaea, heterotrophic bacteria, cyanobacteria); 3–20 µm (nanoplankton: small eukaryotic protists, e.g., diatoms, flagellates, chrysophyta, chlorophyta); 20–200 µm (microplankton: larger eukaryotic protists, e.g., diatoms; protozoa; foraminifera; tintinnids; rotifers; juvenile metazoans); 200–20000 µm (mesozooplankton: metazoans, e.g., copepods and cladocera).

At each station, 10-L tanks are filled with seawater taken directly from the Niskin bottles. From each bottle, seawater is collected using a double sieve in series, the first with a mesh of 200 µm and the second one of 20 µm. After that, using a peristaltic pump and two “tripods” (e.g., stainless steel filters holders of 142 mm diameter), the seawater is filtered in series on two polycarbonate filters (diameter 142 mm, porosity 3 µm and 0.22 µm, respectively) to divide the size fractions 3–20 µm and 0.22–3 µm. The 10 L of water collected is used to obtain two replicates (5 L for replicate) of each size class.

The size fraction 20–200 µm is obtained by rinsing the 20-µm mesh previously used, with 10 L of water resulting from the serial filtration of the fractions 3–20 µm and 0.22–3 µm, recovering in this way the material attached to the mesh. The 10 L of water in the tank containing the 20–200 µm fraction is then filtered in two rounds on two polycarbonate filters of 3 µm (5 L per duplicate). The filters collected from each of the above-mentioned size classes are finally placed into 5 ml cryovials and frozen in liquid nitrogen.

At 0 and 40 m depths, we collect the 0.22–3 µm, 3–20 µm and 20–200 µm size fractions, while we have chosen to collect only the 0.22–3 µm and >3 µm sizes in the deeper station of the Dohrn Canyon using 10 L of water for each of the two replicates due to the limited concentration of material found at 800 m depth.

The 200–2000 µm size fraction is gathered with a double WP2 plankton net, featuring a mouth area of 0.25 m² and a mesh opening of 200 µm, provided with a 500 ml plexiglass filtering cod-end. This net is vertically deployed from a depth of 50 m to the surface at LTER-MC, and both from 200 m to the surface and from 50 m to the surface at Dohrn Canyon. To calculate the volume of water that has been filtered, a mechanical flowmeter from Hydrobios is utilised. After collection, the 500 ml of the mesozooplankton sample are divided into two equal parts of 250 ml each. Both aliquots are filtered through a 200 µm mesh and then fixed in ethanol 95%. One aliquot is used for counting, while the other is used for reference library. The other 500 ml are filtered through a 200 µm mesh, treated, and stored in liquid nitrogen for meta-omics.

Nucleic acid extraction

Nucleic acid extractions are performed according to the protocol used within the Tara Oceans project49. Briefly, it is carried out by cryogenically grinding cryopreserved filters and zooplankton pellets. This process is followed by simultaneous RNA/DNA co-extraction utilising the NucleoSpin RNA and DNA Elution buffer kits from Macherey-Nagel, Düren, Germany. To disintegrate cells, the cryopreserved filters are ground as follows: first, each filter is treated with 1 ml of RA1 lysis buffer and 1% β-mercaptoethanol and then subjected to a grinding program: pre-cooling (2 min), first cycle of 10 knocks/s grinding (1 min), cooling time (1 min) and a final cycle of 10 knocks/s grinding (1 min). The powder resulting from cryogrinding is then subjected to nucleic acid extraction with the NucleoSpin RNA kits (Macherey-Nagel) and DNA Elution buffer kit (Macherey-Nagel). Cryoground powder is resuspended in 2 ml of RA1 lysis buffer mixed with 1% β-mercaptoethanol and then transferred to a NucleoSpin filter from the RNA Midi kit. This solution is centrifuged for 10 minutes at 1,500 g. Following this, an additional 1 ml of RA1 lysis buffer with 1% β-mercaptoethanol is added, and the filter is again centrifuged for 3 minutes at 1,500 g. The resulting eluate is then transferred to a new tube containing an equal volume of 70% ethanol. This mixture is loaded onto a NucleoSpin RNA Mini spin column and washed twice with a DNA washing solution. The DNA is eluted with 100 µl of DNA elution buffer three times and subsequently stored at −20 °C. Finally, the DNA is quantified using a dsDNA-specific fluorometric quantitation method (Broad Range and High Sensitivity Assays) using a Qubit fluorometer. Under such conditions, median quantities of 0.97 and 0.81 μg of DNA and RNA were obtained.

Library preparation from amplicon PCRs

Two markers are targeted for metabarcoding analyses: (i) the hypervariable region V9 of the 18S rDNA for eukaryotes (515 F: 5′ GTG YCA GCM GCC GCG GTA A 3′; 926 R: 5′ CCG YCA ATT YMT TTR AGT TT 3′)50, and (ii) the hypervariable region V4-V5 of the 16S rDNA gene for prokaryotes (1389 F: 5′ TTG TAC ACA CCG CCC 3′; 1510 R: 5′ CCT TCY GCA GGT TCA CCT AC 3′)51. The libraries are constructed using a Barcode IDentifier (BID). In particular, 12 different BIDs are incorporated into the amplification primers, enabling pooling of 6 to 12 PCR products before the library preparation. Then, from pools of PCR products tagged by different BIDs, a single library is constructed and indexed by a NextFlex DNA barcode. The libraries are prepared using 100 ng of amplicon directly end-repaired, A-tailed at 3’ end and ligated with Illumina adaptors using the NEBNext DNA Module and NextFlex DNA barcodes. After two 1x AMPure XP clean ups (only 1 for the 18S), the ligated products are amplified with Kapa Hifi HotStart NGS library Amplification kit and then purified by 1x AMPure XP. Libraries are initially quantified by Quant-it dsDNA HS on a Fluoroskan Ascent instrument (Thermo scientific). Following this, they are quantified by qPCR with the KAPA Library Quantification Kit for Illumina Libraries (Kapa Biosystems) on an MXPro instrument (Agilent Technologies). Library profiles are evaluated using a high throughput microfluidic capillary electrophoresis system (LabChip GX, Perkin Elmer, Waltham, MA). Metabarcoding libraries exhibit low diversity sequences at the beginning of the reads due to the primer sequences used for tag amplification. This low diversity can interfere with the correct cluster identification, leading to a significant loss of data output. For this reason, loading concentrations of these libraries and PhiX DNA spike-in (20% instead of 1%) are adapted to minimise the impacts on the run quality. Metabarcoding libraries are sequenced on a NovaSeq 6000 instrument (Illumina, San Diego, CA, USA) using 2 × 150 or 250 base-length read chemistry for 18SV9 or 16SV4-V5 libraries, respectively.

Library preparation for metagenomics

DNA is first fragmented to a size range of 100–1,000 bp using the E220 Covaris Focused ultrasonicator (Covaris, Woburn, MA, USA). These fragments undergo end-repair and are then adenylated at the 3′ end. NEXTflex HT Barcodes (Bioo Scientific, Austin, TX, USA) are subsequently added using products from the NEBNext DNA modules (New England Biolabs, Ipswich, MA, USA). The library construction process includes two consecutive clean-up steps using 1 × AMPure XP (Beckman Coulter, Brea, CA, USA). After these clean-ups, the ligated product is amplified through 12 PCR cycles using the Kapa Hifi Hotstart NGS library amplification kit (Kapa Biosystems, Wilmington, MA, USA), followed by 0.6x AMPure XP purification. Metagenomics libraries, similarly to metabarcoding libraries, are quantified using Quant-iT dsDNA HS on a Fluoroskan Ascent instrument (Thermo scientific) and by qPCR with the KAPA Library Quantification Kit for Illumina Libraries (Kapa Biosystems) on an MXPro instrument (Agilent Technologies). Subsequently, library profiles are assessed with a high-throughput microfluidic capillary electrophoresis system (LabChip GX, Perkin Elmer, Waltham, MA). Metagenomic libraries are sequenced by a NovaSeq 6000 instrument (Illumina, San Diego, CA, USA) by paired ends of 150 bases in length.

Metabarcoding - 16S ASV generation and taxonomic assignment

Raw 16S paired-end sequences were subjected to a data quality control step49 and subsequently imported into the QIIME2 pipeline v.2022.2.052. Leftover primers and adapter sequences were removed through cutadapt53. The Amplicon Sequence Variants (ASVs) table, which represents true biological sequences within each sample, was generated using the denoised-paired method including truncation, denoising, dereplication, merging, and chimera filtering of the DADA2 (Divisive Amplicon Denoising Algorithm 2) plugin inside QIIME254. Default parameters were used with the exception of the forward and reverse sequence length (–p-trunc-len-f and–p-trunc-len-r), that were set to 220 and 180, respectively. Processed reads that passed all these filters were used for taxonomy classification (Fig. 4a). The V4-V5 regions were extracted from the pre-formatted reference sequences and taxonomy file built on the SILVA 138 99% OTUs database55–57 and the vsearch (v.2.6.2) global alignment implemented in QIIME2 was used58.Fig. 4 Class-level abundance of Amplicon Sequence Variants (ASVs) and Metagenome-Assembled Genomes (MAGs). The panels show the log-scaled read counts for 16S (a) and 18S (b) ASVs and the number of prokaryotic MAGs (c) ASVs and MAGs representing more than 0.5% of the total abundance and assigned to a known class are displayed.

Metabarcoding - 18S ASV generation and taxonomic assignment

Illumina paired-end V9-18S raw reads (FASTQ format, 2 × 150 PE) were pre-processed with cutadapt53 and vsearch to remove primer sequences, trim low quality bases and unify mixed orientation reads produced in the ligation-based library preparation; the procedure was implemented in a custom bash script. Processed reads were then used to generate ASVs using the DADA2 R library54; the pipeline was adapted from the one described on the program website (https://benjjneb.github.io/dada2/tutorial.html); no further quality filtering was implemented at this stage, except for discarding all reads with ambiguities (parameter maxN = 0 of function filterAndTrim). Filtered forward (F) and reverse (R) reads were used to train the error model and then denoised by applying the trained error model to generate ASVs. Finally, F and R reads were merged and checked for chimeras, allowing no mismatches in read merging (default parameter maxMismatch = 0 of function mergePairs). ASVs were then classified with BLAST against the PR2 (v.5.01) reference database59, integrated with 1,293 environmental sequences from GoN protist strains and fungi. Highest bit scores matching with the best taxonomic resolution were then selected among the returned results (Fig. 4b).

Metagenomic data processing

Metagenomic data processing was performed as described in60. Briefly, all sequencing reads were quality filtered using BBMap (v.38.79). This step involved removing sequencing adapters, eliminating reads that matched the PhiX control genome, and discarding low-quality reads based on specific criteria (trimq = 14, maq = 20, maxns = 1, minlength = 45). For downstream analyses, either quality-controlled reads were used, or in some cases, reads were merged using bbmerge.sh with a minimum overlap of 16. All metagenomes were assembled individually with metaSPAdes (v.3.15.2)61. The assembled scaffolds were filtered to retain only those longer than 0.5 kbp. Gene prediction on these scaffolds was performed using Prodigal (v.2.6.3)62 with the parameters -c -q -m -p meta. For MAGs reconstruction, quality-controlled metagenomic reads from all samples were mapped back to the scaffolds (≥1 kbp) of each sample. The mapping was executed using BWA (v.0.7.17-r1188)63, configured to allow reads to map at secondary sites with the use of the -a flag. This step ensured that alignments were only considered if they were at least 45 bases long, had an identity of at least 97%, and provided a coverage of 80% or more of the read sequence. The alignments were then compiled into BAM files, which were processed using the jgi_summarize_bam_contig_depths script from MetaBAT2 (v.2.12.1)64. This script facilitated the calculation of within-sample and between-sample coverages for each scaffold. Finally, the scaffolds were binned using MetaBAT2, applying parameters –minContig 2000 and–maxEdges 500 to enhance the sensitivity of the binning process. Metagenomic bin quality was assessed using two methods: the ‘lineage workflow’ of CheckM (v.1.1.3)65 and Anvi’o (v.7.1)66. Bins were retained for further analysis if they met the following criteria from either tool: a completeness of at least 50% and a contamination level of 10% or less. Following the quality assessment, 1070 prokaryotic MAGs were taxonomically classified using GTDB-Tk (v.2.1.0)67 utilising default parameters against the GTDB R207 release68 (Fig. 4c). Gene prediction on these MAGs was carried out using Prodigal (v.2.6.3)62 with parameters -c -q -m -p single. Gene sequences predicted from prokaryotic MAGs of all samples and gene sequences predicted from assemblies of prokaryote-enriched samples were clustered at 95% identity, retaining the longest sequence within each cluster as the representative using CD-HIT (v.4.8.1)69. The specific parameters used for this clustering were -c 0.95 -M 0 -G 0 -aS 0.9 -g 1 -r 1 -d 0. The representative gene sequences obtained from the clustering were then aligned against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, using the DIAMOND aligner (v.2.0.15)70. The filtering criteria for these alignments included a minimum query and subject coverage of 70%, and a requirement for the bitScore to be at least 50% of the maximum expected bitScore, which references the score against itself. The 59 metagenomes were mapped to the 5,515,266 million cluster representatives using BWA (v.0.7.17-r1188)63 with the -a parameter to allow alignment at secondary sites. The resulting BAM files from these alignments were then refined to retain only the ones exhibiting a minimum percentage identity of 95% and a minimum alignment length of 45 bases. For the calculation of length-normalised gene abundance, the process began by counting inserts from the best unique alignments. In cases of ambiguously mapped inserts, fractional counts were allocated to the corresponding target genes based on their unique insert abundances. Subsequently, the total insert counts were divided by the length of each respective gene to normalise the gene abundance data accurately. Gene-length normalised read abundances were further converted into per-cell gene copy numbers, dividing them by the median abundance of 10 single-copy marker gene copies71 in each sample. For taxonomic profiling, the mOTUs database (v.3.1)71 was extended with 1070 prokaryotic MAGs forming 433 new species level clusters. Quality controlled sequencing reads of 59 metagenomes were profiled using the extended mOTUs database with default parameters.

Data Records

We assign a specific label to each NEREA sampling event. We define the station code using the letters NR for all the sites followed by MC, S and C respectively for MareChiara, Sarno and Dohrn Canyon. The sampling code at the MareChiara site follows the specific LTER-MC code campaign to unify the two samplings when conducted simultaneously.

Data described here are available in the Sample Registry “NEREA- Naples Ecological REsearch for Augmented observatories” on Zenodo (https://zenodo.org/communities/nerea/records), which contains the following projects:Physical Oceanography: CTD and Thermosalinograph72 (https://zenodo.org/records/10986789);

Biogeochemical data73 (https://zenodo.org/records/11035857);

LOFAs data: particulate lipid-derived oxylipins74 (https://zenodo.org/records/11058615);

HPLC data75 (https://zenodo.org/records/11035381);

Carbon and nitrogen stable isotope ratios of plankton76 (https://zenodo.org/records/12743875);

Plankton biodiversity data: Flow cytometry data77 (https://zenodo.org/records/11035933);

Plankton biodiversity data: Phytoplankton78 (https://zenodo.org/records/10987241);

Plankton biodiversity data: Microzooplankton79 (https://zenodo.org/records/10987215);

Plankton biodiversity data: Mesozooplankton80 (https://zenodo.org/records/10987255).

The datasets included in the above-listed projects share a similar structure: they are all xlsx files including three specific sheets:a “Readme”, which contains the dataset title, a brief description of the dataset and the methods used to produce it, and any additional note;

a “Coordinates” sheet, including the coordinates of all the sampling stations included in the datasets;

a “Data” sheet with the data.

Prokaryotic and eukaryotic metabarcoding Zenodo projects include xlsx files structured as described above. Moreover, they contain additional files with bioinformatic scripts and intermediates, including tsv and fasta files:

Metabarcoding data: 16S ASVs abundance table and taxonomic assignment81 (https://zenodo.org/records/12801913);

Metabarcoding data: 18S ASVs abundance table and taxonomic assignment82 (https://zenodo.org/records/12801941).

Metagenomic data are grouped in two Zenodo projects. The first project contains:the full prokaryotic gene catalogue (tsv);

the Metagenome-Assembled Genomes (MAGs) sequences (fasta) and annotation (tsv);

the Taxonomic profiling using mOTUs (MG-based operational taxonomic units) database83 (https://zenodo.org/records/11035656).

The second project contains the single-assemblies and the genes (fasta and gff files)84 (https://zenodo.org/records/11046519).

Plankton metabarcoding and metagenomic raw reads were registered in the European Nucleotide Archive (ENA) at EMBL-EBI, under the umbrella project permanent identifier PRJEB7464985 (https://www.ebi.ac.uk/ena/browser/view/PRJEB74649). Raw sequencing data are recorded in three component projects all linked to the umbrella project:(i) Prokaryotic metabarcoding sequencing data: The 16S raw sequence data files (FastQ format) are available under the BioProject ID PRJEB7464186 (https://www.ebi.ac.uk/ena/browser/view/PRJEB74641);

(ii) Eukaryotic metabarcoding sequencing data: the 18S raw sequence data files (FastQ format) are available under the project ID PRJEB7465887 (https://www.ebi.ac.uk/ena/browser/view/PRJEB74658);

(iii) Metagenomics: Raw sequence data files (FastQ format) are available under the project ID PRJEB7465988 (https://www.ebi.ac.uk/ena/browser/view/PRJEB74659).

Technical Validation

The NEREA sampling methodology follows the one used during the Tara Oceans expedition32. This methodology is based on the study from the literature of size, abundance, and richness values of all the plankton organisms used to characterise each size fraction division. DNA obtained from extraction is quantified by a dsDNA-specific fluorometric method (Broad Range and High Sensitivity Assays) using a Qubit Fluorometer instrument. All metabarcoding libraries are quantified by measurement of the Qubit dsDNA HS Assay and then a size profile analysis is performed using an Agilent 2100 Bioanalyzer and through qPCR with the KAPA Library Quantification Kit for Illumina Libraries on an MXPro instrument. For metagenomics, the size profiles are visualised using an Agilent Bioanalyzer DNA High Sensitivity chip.

Acknowledgements

The NEREA activities have received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 101000858 (TechOceanS), No 862626 (EuroSea), No 101081642 (OBAMA-NEXT). NEREA is an UN Ocean Decade Action No 64.2 under the OBON Program. NEREA is also partly funded under the National Recovery and Resilience Plan (NRRP), Mission 4 Component 2 Investment 1.4 - Call for tender No. 3138 of 16 December 2021, rectified by Decree n.3175 of 18 December 2021 of Italian Ministry of University and Research funded by the European Union – NextGenerationEU. The Long-Term Ecological Research site MareChiara program is funded by Stazione Zoologica Anton Dohrn, Naples (Italy). L.C. acknowledges MARCO-BOLO (ID: 101082021) and AtlantECO (ID: 862923) projects. R.C. acknowledges the OBAMA-NEXT (ID: 101081642) project. L.R. has been supported by the BiOcean5D project (ID: 101059915). V. D., L.M., E.R., and D. GV. acknowledge a fellowship funded by the Stazione Zoologica Anton Dohrn (SZN) within the SZN-Open University Ph.D. program. E.R. was also supported by the project “Antitumor Drugs and Vaccines from the Sea (ADViSE)” (CUP B43D18000240007–SURF 17061BP000000011; PG/2018/0494374) funded by POR Campania FESR 2014–2020 “Technology Platform for Therapeutic Strategies against Cancer”- Action 1.1.2 and 1.2.2. The sequencing lab at Genoscope acknowledges financial support from France Génomique National infrastructure, funded as part of “Investissement d'Avenir” program managed by Agence Nationale pour la Recherche (contract ANR‐10‐INBS‐09). S.S. acknowledges funding from the Swiss National Science Foundation (205320_215395) and support from the ETH IT services for calculations that were carried out on the ETH Euler cluster. The authors thank Alexandre Epinoux, Carmen Minucci and Francesco Manfellotto for technical assistance.

Author contributions

Conceptualization: R.C., F.C., D.DA., D.I., S.L., F.M., M.Mo.; Funding: R.C., D.I.; Sampling design and data collection: U.C., Y.C., R.C., D.DA., I. DA., I.DC., V.D., R.G., D.I., P.L., F.M., L.M., M.Mo., C.M., A.Pa., I.P., E.R., L.R., D.S., F.T., AC.T., D.GV.; Data processing: M.A., A.A., G.B., C.B., A.B., Y.C., R.C., G.C., I. DA., I.DC., G.dI., R.G., M.F., K.L., F.M., M.Ma., PH.O., A.Pa., I.P., A.Pe., R.P., F.R., HJ.R., E.R., M.S., C.S., D.S., AC.T., G.Z.; Data curation: A.A., G.B., C.B., D.B., L.C., Y.C., R.C., G.C., I. DA., I.DC., G.dI., R.G., M.F., K.L., A.L., F.M., M.Ma., PH.O., A.Pa., I.P., A.Pe., R.P., F.R., HJ.R., L.R., M.S., D.S., AC.T., G.Z.; Data management: L.C., R.C., F.M., L.R., AC.T.; Visualisation: L.C., D.DA., A.L., L.R.; Writing - Original Draft: R.C., D.DA., M.Mo., L.R.; Writing - Review and Editing: G.B., D.B., L.C., Y.C., ML. C., R.C., D.DA., I. DA., I.DC., R.G., K.L., F.M., L.M., PH.O., A.Pa., I.P., F.R., V.R., HJ.R., L.R., M.S., D.S., S.S., G.Z.; Supervision: Y.C., R.C., ML.C., D.DA., A.F., D.I., F.M., M.Mo., R.P., D.S., S.S., P.W.

Code availability

Coding scripts for 16S and 18S ASVs generation and taxonomic assignment can be found at [https://zenodo.org/records/12801913] and [https://zenodo.org/records/12801941], respectively.

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.

These authors contributed equally: Lucia Campese, Luca Russo.

These authors jointly supervised this work: Raffaella Casotti, Fabio Conversano, Domenico D’Alelio, Daniele Iudicone, Francesca Margiotta, Marina Montresor.
==== Refs
References

1. Barbier EB Marine ecosystem services Curr. Biol. 2017 27 R507 R510 10.1016/j.cub.2017.03.020 28586688
Barbier, E. B. Marine ecosystem services. Curr. Biol. 27, R507–R510 (2017).28586688 10.1016/j.cub.2017.03.020
2. Mirtl M Genesis, goals and achievements of Long-Term Ecological Research at the global scale: A critical review of ILTER and future directions Sci. Total Environ. 2018 626 1439 1462 10.1016/j.scitotenv.2017.12.001 29898550
Mirtl, M. et al. Genesis, goals and achievements of Long-Term Ecological Research at the global scale: A critical review of ILTER and future directions. Sci. Total Environ. 626, 1439–1462 (2018).29898550 10.1016/j.scitotenv.2017.12.001
3. Hughes BB Long-Term Studies Contribute Disproportionately to Ecology and Policy BioScience 2017 67 271 281 10.1093/biosci/biw185
Hughes, B. B. et al. Long-Term Studies Contribute Disproportionately to Ecology and Policy. BioScience 67, 271–281 (2017).10.1093/biosci/biw185
4. Wells HBM Dougill AJ Stringer LC The importance of long‐term social‐ecological research for the future of restoration ecology Restor. Ecol. 2019 27 929 933 10.1111/rec.13000
Wells, H. B. M., Dougill, A. J. & Stringer, L. C. The importance of long‐term social‐ecological research for the future of restoration ecology. Restor. Ecol. 27, 929–933 (2019).10.1111/rec.13000
5. Occhipinti-Ambrogi A Alien species along the Italian coasts: an overview Biol. Invasions 2011 13 215 237 10.1007/s10530-010-9803-y
Occhipinti-Ambrogi, A. et al. Alien species along the Italian coasts: an overview. Biol. Invasions 13, 215–237 (2011).10.1007/s10530-010-9803-y
6. Iacono R Napolitano E Palma M Sannino G The Tyrrhenian Sea Circulation: A Review of Recent Work Sustainability 2021 13 6371 10.3390/su13116371
Iacono, R., Napolitano, E., Palma, M. & Sannino, G. The Tyrrhenian Sea Circulation: A Review of Recent Work. Sustainability 13, 6371 (2021).10.3390/su13116371
7. Zingone A Time series and beyond: multifaceted plankton research at a marine Mediterranean LTER site Nat. Conserv. 2019 34 273 310 10.3897/natureconservation.34.30789
Zingone, A. et al. Time series and beyond: multifaceted plankton research at a marine Mediterranean LTER site. Nat. Conserv. 34, 273–310 (2019).10.3897/natureconservation.34.30789
8. Tornero V Ribera d’Alcalà M Contamination by hazardous substances in the Gulf of Naples and nearby coastal areas: A review of sources, environmental levels and potential impacts in the MSFD perspective Sci. Total Environ. 2014 466–467 820 840 10.1016/j.scitotenv.2013.06.106 23994731
Tornero, V. & Ribera d’Alcalà, M. Contamination by hazardous substances in the Gulf of Naples and nearby coastal areas: A review of sources, environmental levels and potential impacts in the MSFD perspective. Sci. Total Environ. 466–467, 820–840 (2014).23994731 10.1016/j.scitotenv.2013.06.106
9. Summary, E.G7 Future Seas Oceans Working Group. Progress since may 2016. Turin 2017.
10. Halpern BS An index to assess the health and benefits of the global ocean Nature 2012 488 615 620 10.1038/nature11397 22895186
Halpern, B. S. et al. An index to assess the health and benefits of the global ocean. Nature 488, 615–620 (2012).22895186 10.1038/nature11397
11. Karsenti E A Holistic Approach to Marine Eco-Systems Biology PLOS Biol. 2011 9 e1001177 10.1371/journal.pbio.1001177 22028628
Karsenti, E. et al. A Holistic Approach to Marine Eco-Systems Biology. PLOS Biol. 9, e1001177 (2011).22028628 10.1371/journal.pbio.1001177
12. Ruggiero MA A Higher Level Classification of All Living Organisms PLOS ONE 2015 10 e0119248 10.1371/journal.pone.0119248 25923521
Ruggiero, M. A. et al. A Higher Level Classification of All Living Organisms. PLOS ONE 10, e0119248 (2015).25923521 10.1371/journal.pone.0119248
13. Boyce DG Frank KT Leggett WC From mice to elephants: overturning the ‘one size fits all’ paradigm in marine plankton food chains Ecol. Lett. 2015 18 504 515 10.1111/ele.12434 25919397
Boyce, D. G., Frank, K. T. & Leggett, W. C. From mice to elephants: overturning the ‘one size fits all’ paradigm in marine plankton food chains. Ecol. Lett. 18, 504–515 (2015).25919397 10.1111/ele.12434
14. Falkowski PG Fenchel T Delong EF The microbial engines that drive Earth’s biogeochemical cycles Science 2008 320 1034 1039 10.1126/science.1153213 18497287
Falkowski, P. G., Fenchel, T. & Delong, E. F. The microbial engines that drive Earth’s biogeochemical cycles. Science 320, 1034–1039 (2008).18497287 10.1126/science.1153213
15. Field CB Behrenfeld MJ Randerson JT Falkowski P Primary Production of the Biosphere: Integrating Terrestrial and Oceanic Components Science 1998 281 237 240 10.1126/science.281.5374.237 9657713
Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary Production of the Biosphere: Integrating Terrestrial and Oceanic Components. Science 281, 237–240 (1998).9657713 10.1126/science.281.5374.237
16. Fanelli E Towards Naples Ecological REsearch for Augmented Observatories (NEREA): The NEREA-Fix Module, a Stand-Alone Platform for Long-Term Deep-Sea Ecosystem Monitoring Sensors 2020 20 2911 10.3390/s20102911 32455611
Fanelli, E. et al. Towards Naples Ecological REsearch for Augmented Observatories (NEREA): The NEREA-Fix Module, a Stand-Alone Platform for Long-Term Deep-Sea Ecosystem Monitoring. Sensors 20, 2911 (2020).32455611 10.3390/s20102911
17. Cianelli, D. et al. Dynamics of a very special Mediterranean coastal area: the Gulf of Naples. Mediterr. Ecosyst. Dyn. Conserv. 129–150 (2012).
18. Cianelli D Inshore/offshore water exchange in the Gulf of Naples J. Mar. Syst. 2015 145 37 52 10.1016/j.jmarsys.2015.01.002
Cianelli, D. et al. Inshore/offshore water exchange in the Gulf of Naples. J. Mar. Syst. 145, 37–52 (2015).10.1016/j.jmarsys.2015.01.002
19. Cianelli D Disentangling physical and biological drivers of phytoplankton dynamics in a coastal system Sci. Rep. 2017 7 15868 10.1038/s41598-017-15880-x 29158517
Cianelli, D. et al. Disentangling physical and biological drivers of phytoplankton dynamics in a coastal system. Sci. Rep. 7, 15868 (2017).29158517 10.1038/s41598-017-15880-x
20. Di Capua, I. et al. From Phenotypes to Genotypes and Back: Toward an Integrated Evaluation of Biodiversity in Calanoid Copepods. Front. Mar. Sci. 9, 833089 (2022).
21. Carrada GC Variability in the Hydrographic and Biological Features of the Gulf of Naples Mar. Ecol. 1980 1 105 120 10.1111/j.1439-0485.1980.tb00213.x
Carrada, G. C. et al. Variability in the Hydrographic and Biological Features of the Gulf of Naples. Mar. Ecol. 1, 105–120 (1980).10.1111/j.1439-0485.1980.tb00213.x
22. Milia A The Dohrn canyon: a response to the eustatic fall and tectonic uplift of the outer shelf along the eastern Tyrrhenian Sea margin, Italy Geo-Mar. Lett. 2000 20 101 108 10.1007/s003670000044
Milia, A. The Dohrn canyon: a response to the eustatic fall and tectonic uplift of the outer shelf along the eastern Tyrrhenian Sea margin, Italy. Geo-Mar. Lett. 20, 101–108 (2000).10.1007/s003670000044
23. De Pippo T Donadio C Guida M Petrosino C The Case of Sarno River (Southern Italy). Effects of geomorphology on the environmental impacts (8 pp) Environ. Sci. Pollut. Res. 2006 13 184 191 10.1065/espr2005.08.287
De Pippo, T., Donadio, C., Guida, M. & Petrosino, C. The Case of Sarno River (Southern Italy). Effects of geomorphology on the environmental impacts (8 pp). Environ. Sci. Pollut. Res. 13, 184–191 (2006).10.1065/espr2005.08.287
24. Thiele S Richter M Balestra C Glöckner FO Casotti R Taxonomic and functional diversity of a coastal planktonic bacterial community in a river-influenced marine area Mar. Genomics 2017 32 61 69 10.1016/j.margen.2016.12.003 28063827
Thiele, S., Richter, M., Balestra, C., Glöckner, F. O. & Casotti, R. Taxonomic and functional diversity of a coastal planktonic bacterial community in a river-influenced marine area. Mar. Genomics 32, 61–69 (2017).28063827 10.1016/j.margen.2016.12.003
25. Aiello G Iorio M Molisso F Sacchi M Integrated Morpho-Bathymetric, Seismic-Stratigraphic, and Sedimentological Data on the Dohrn Canyon (Naples Bay, Southern Tyrrhenian Sea): Relationships with Volcanism and Tectonics Geosciences 2020 10 319 10.3390/geosciences10080319
Aiello, G., Iorio, M., Molisso, F. & Sacchi, M. Integrated Morpho-Bathymetric, Seismic-Stratigraphic, and Sedimentological Data on the Dohrn Canyon (Naples Bay, Southern Tyrrhenian Sea): Relationships with Volcanism and Tectonics. Geosciences 10, 319 (2020).10.3390/geosciences10080319
26. Taviani M Angeletti L Cardone F Montagna P Danovaro R A unique and threatened deep water coral-bivalve biotope new to the Mediterranean Sea offshore the Naples megalopolis Sci. Rep. 2019 9 3411 10.1038/s41598-019-39655-8 30833591
Taviani, M., Angeletti, L., Cardone, F., Montagna, P. & Danovaro, R. A unique and threatened deep water coral-bivalve biotope new to the Mediterranean Sea offshore the Naples megalopolis. Sci. Rep. 9, 3411 (2019).30833591 10.1038/s41598-019-39655-8
27. Ciannelli L Ichthyoplankton assemblages and physical characteristics of two submarine canyons in the south central Tyrrhenian Sea Fish. Oceanogr. 2022 31 480 496 10.1111/fog.12596
Ciannelli, L. et al. Ichthyoplankton assemblages and physical characteristics of two submarine canyons in the south central Tyrrhenian Sea. Fish. Oceanogr. 31, 480–496 (2022).10.1111/fog.12596
28. Ribera d’Alcalà M Seasonal patterns in plankton communities in a pluriannual time series at a coastal Mediterranean site (Gulf of Naples): an attempt to discern recurrences and trends Sci. Mar. 2004 68 65 83 10.3989/scimar.2004.68s165
Ribera d’Alcalà, M. et al. Seasonal patterns in plankton communities in a pluriannual time series at a coastal Mediterranean site (Gulf of Naples): an attempt to discern recurrences and trends. Sci. Mar. 68, 65–83 (2004).10.3989/scimar.2004.68s165
29. Grasshoff, K., Kremlingl, K. & Ehrhardt, M. Methods of seawater analysis. 347–395 (1983).
30. Valderrama JC The simultaneous analysis of total nitrogen and total phosphorus in natural waters Mar. Chem. 1981 10 109 122 10.1016/0304-4203(81)90027-X
Valderrama, J. C. The simultaneous analysis of total nitrogen and total phosphorus in natural waters. Mar. Chem. 10, 109–122 (1981).10.1016/0304-4203(81)90027-X
31. Sugimura, Y. & Suzuki, Y. A high-temperature catalytic oxidation method for the determination of non-volatile dissolved organic carbon in seawater by direct injection of a liquid sample. Mar. Chem. 24, 105–131 (1988).
32. Hedges JI Stern JH Carbon and nitrogen determinations of carbonate-containing solids1 Limnol. Oceanogr. 1984 29 657 663 10.4319/lo.1984.29.3.0657
Hedges, J. I. & Stern, J. H. Carbon and nitrogen determinations of carbonate-containing solids1. Limnol. Oceanogr. 29, 657–663 (1984).10.4319/lo.1984.29.3.0657
33. Russo E Density-dependent oxylipin production in natural diatom communities: possible implications for plankton dynamics ISME J. 2020 14 164 177 10.1038/s41396-019-0518-5 31611654
Russo, E. et al. Density-dependent oxylipin production in natural diatom communities: possible implications for plankton dynamics. ISME J. 14, 164–177 (2020).31611654 10.1038/s41396-019-0518-5
34. Holm-Hansen O Lorenzen CJ Holmes RW Strickland JDH Fluorometric Determination of Chlorophyll ICES J. Mar. Sci. 1965 30 3 15 10.1093/icesjms/30.1.3
Holm-Hansen, O., Lorenzen, C. J., Holmes, R. W. & Strickland, J. D. H. Fluorometric Determination of Chlorophyll. ICES J. Mar. Sci. 30, 3–15 (1965).10.1093/icesjms/30.1.3
35. Vidussi F Claustre H Bustillosguzman J Cailliau C Marty J Determination of chlorophylls and carotenoids of marine phytoplankton: Separation of chlorophyll a from divinyl-chlorophyll a and zeaxanthin from lutein J. Plankton Res. 1996 18 2377 2382 10.1093/plankt/18.12.2377
Vidussi, F., Claustre, H., Bustillosguzman, J., Cailliau, C. & Marty, J. Determination of chlorophylls and carotenoids of marine phytoplankton: Separation of chlorophyll a from divinyl-chlorophyll a and zeaxanthin from lutein. J. Plankton Res. 18, 2377–2382 (1996).10.1093/plankt/18.12.2377
36. Brunet, C. & Mangoni, O. Determinazione quali-quantitativa dei pigmenti fitoplanctonici mediante HPLC. vol. 56 (Roma, 2010).
37. Marie D Partensky F Vaulot D Brussaard C Enumeration of Phytoplankton, Bacteria, and Viruses in Marine Samples Curr. Protoc. Cytom. 1999 10 11.11.1 11.11.15
Marie, D., Partensky, F., Vaulot, D. & Brussaard, C. Enumeration of Phytoplankton, Bacteria, and Viruses in Marine Samples. Curr. Protoc. Cytom. 10, 11.11.1–11.11.15 (1999).
38. Balestra C Alonso-Sáez L Gasol JM Casotti R Group-specific effects on coastal bacterioplankton of polyunsaturated aldehydes produced by diatoms Aquat. Microb. Ecol. 2011 63 123 131 10.3354/ame01486
Balestra, C., Alonso-Sáez, L., Gasol, J. M. & Casotti, R. Group-specific effects on coastal bacterioplankton of polyunsaturated aldehydes produced by diatoms. Aquat. Microb. Ecol. 63, 123–131 (2011).10.3354/ame01486
39. Elder, L. & Elbrachter, M. The Utermöhl method for quantitative phytoplankton analysis. Microscopic and molecular methods for quantitative phytoplankton analysis. in Microscopic and Molecular Methods for Quantitative Phytoplankton Analysis. UNESCO 13–20 (2010).
40. Cupp, E. E. Marine Plankton Diatoms of the West Coast of North America. vol. 5 (University of California, Berkeley, California and London, England, 1943).
41. Monks Wood Experimental Station. Provisional Atlas of the Marine Dinoflagellates of the British Isles. (Monks Wood Experimental Station, Huntingdon, 1981).
42. Rampi, L. & Zattera, A. Chiave per la determinazione dei Tintinnidi mediterranei. (ENEA, Roma, 1982).
43. Tomas, C. R. Identifying Marine Phytoplankton. (San Diego, 1997).
44. Lee AJ Folkard AR Factors Affecting Turbidity in the Southern North Sea ICES J. Mar. Sci. 1969 32 291 302 10.1093/icesjms/32.3.291
Lee, A. J. & Folkard, A. R. Factors Affecting Turbidity in the Southern North Sea. ICES J. Mar. Sci. 32, 291–302 (1969).10.1093/icesjms/32.3.291
45. Strüder-Kypke, M., Dr. Agatha, S., Warwick, J. & Montagnes, D. The user-friendly key to coastal planktonic ciliates. (2002).
46. van Guelpen L Markle DF Duggan DJ An evaluation of accuracy, precision, and speed of several zooplankton subsampling techniques ICES J. Mar. Sci. 1982 40 226 236 10.1093/icesjms/40.3.226
van Guelpen, L., Markle, D. F. & Duggan, D. J. An evaluation of accuracy, precision, and speed of several zooplankton subsampling techniques. ICES J. Mar. Sci. 40, 226–236 (1982).10.1093/icesjms/40.3.226
47. Pesant S Open science resources for the discovery and analysis of Tara Oceans data Sci. Data 2015 2 150023 10.1038/sdata.2015.23 26029378
Pesant, S. et al. Open science resources for the discovery and analysis of Tara Oceans data. Sci. Data 2, 150023 (2015).26029378 10.1038/sdata.2015.23
48. Gorsky, G. et al. Expanding Tara Oceans Protocols for Underway, Ecosystemic Sampling of the Ocean-Atmosphere Interface During Tara Pacific Expedition (2016–2018). Front. Mar. Sci. 6, (2019).
49. Belser C Integrative omics framework for characterization of coral reef ecosystems from the Tara Pacific expedition Sci. Data 2023 10 326 10.1038/s41597-023-02204-0 37264047
Belser, C. et al. Integrative omics framework for characterization of coral reef ecosystems from the Tara Pacific expedition. Sci. Data 10, 326 (2023).37264047 10.1038/s41597-023-02204-0
50. Amaral-Zettler LA McCliment EA Ducklow HW Huse SM A Method for Studying Protistan Diversity Using Massively Parallel Sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA Genes PLOS ONE 2009 4 e6372 10.1371/journal.pone.0006372 19633714
Amaral-Zettler, L. A., McCliment, E. A., Ducklow, H. W. & Huse, S. M. A Method for Studying Protistan Diversity Using Massively Parallel Sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA Genes. PLOS ONE 4, e6372 (2009).19633714 10.1371/journal.pone.0006372
51. Parada AE Needham DM Fuhrman JA Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples Environ. Microbiol. 2016 18 1403 1414 10.1111/1462-2920.13023 26271760
Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).26271760 10.1111/1462-2920.13023
52. Bolyen E Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2 Nat. Biotechnol. 2019 37 852 857 10.1038/s41587-019-0209-9 31341288
Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).31341288 10.1038/s41587-019-0209-9
53. Martin M Cutadapt removes adapter sequences from high-throughput sequencing reads EMBnet.journal 2011 17 10 12 10.14806/ej.17.1.200
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17, 10–12 (2011).10.14806/ej.17.1.200
54. Callahan BJ DADA2: High-resolution sample inference from Illumina amplicon data Nat. Methods 2016 13 581 583 10.1038/nmeth.3869 27214047
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).27214047 10.1038/nmeth.3869
55. Quast C The SILVA ribosomal RNA gene database project: improved data processing and web-based tools Nucleic Acids Res. 2013 41 D590 D596 10.1093/nar/gks1219 23193283
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).23193283 10.1093/nar/gks1219
56. Yilmaz P The SILVA and ‘All-species Living Tree Project (LTP)’ taxonomic frameworks Nucleic Acids Res. 2014 42 D643 648 10.1093/nar/gkt1209 24293649
Yilmaz, P. et al. The SILVA and ‘All-species Living Tree Project (LTP)’ taxonomic frameworks. Nucleic Acids Res. 42, D643–648 (2014).24293649 10.1093/nar/gkt1209
57. Robeson MS RESCRIPt: Reproducible sequence taxonomy reference database management PLOS Comput. Biol. 2021 17 e1009581 10.1371/journal.pcbi.1009581 34748542
Robeson, M. S. et al. RESCRIPt: Reproducible sequence taxonomy reference database management. PLOS Comput. Biol. 17, e1009581 (2021).34748542 10.1371/journal.pcbi.1009581
58. Rognes T Flouri T Nichols B Quince C Mahé F VSEARCH: a versatile open source tool for metagenomics PeerJ 2016 4 e2584 10.7717/peerj.2584 27781170
Rognes, T., Flouri, T., Nichols, B., Quince, C. & Mahé, F. VSEARCH: a versatile open source tool for metagenomics. PeerJ 4, e2584 (2016).27781170 10.7717/peerj.2584
59. Guillou L The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy Nucleic Acids Res. 2013 41 D597 D604 10.1093/nar/gks1160 23193267
Guillou, L. et al. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote Small Sub-Unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D597–D604 (2013).23193267 10.1093/nar/gks1160
60. Paoli L Biosynthetic potential of the global ocean microbiome Nature 2022 607 111 118 10.1038/s41586-022-04862-3 35732736
Paoli, L. et al. Biosynthetic potential of the global ocean microbiome. Nature 607, 111–118 (2022).35732736 10.1038/s41586-022-04862-3
61. Nurk S Meleshko D Korobeynikov A Pevzner PA metaSPAdes: a new versatile metagenomic assembler Genome Res. 2017 27 824 834 10.1101/gr.213959.116 28298430
Nurk, S., Meleshko, D., Korobeynikov, A. & Pevzner, P. A. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834 (2017).28298430 10.1101/gr.213959.116
62. Hyatt D Prodigal: prokaryotic gene recognition and translation initiation site identification BMC Bioinformatics 2010 11 119 10.1186/1471-2105-11-119 20211023
Hyatt, D. et al. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11, 119 (2010).20211023 10.1186/1471-2105-11-119
63. Li H Durbin R Fast and accurate short read alignment with Burrows–Wheeler transform Bioinformatics 2009 25 1754 1760 10.1093/bioinformatics/btp324 19451168
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).19451168 10.1093/bioinformatics/btp324
64. Kang DD MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies PeerJ 2019 7 e7359 10.7717/peerj.7359 31388474
Kang, D. D. et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7, e7359 (2019).31388474 10.7717/peerj.7359
65. Parks DH Imelfort M Skennerton CT Hugenholtz P Tyson GW CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes Genome Res. 2015 25 1043 1055 10.1101/gr.186072.114 25977477
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055 (2015).25977477 10.1101/gr.186072.114
66. Eren AM Anvi’o: an advanced analysis and visualization platform for ‘omics data PeerJ 2015 3 e1319 10.7717/peerj.1319 26500826
Eren, A. M. et al. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ 3, e1319 (2015).26500826 10.7717/peerj.1319
67. Chaumeil P-A Mussig AJ Hugenholtz P Parks DH GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database Bioinformatics 2020 36 1925 1927 10.1093/bioinformatics/btz848
Chaumeil, P.-A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 36, 1925–1927 (2020).10.1093/bioinformatics/btz848
68. Parks DH A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life Nat. Biotechnol. 2018 36 996 1004 10.1038/nbt.4229 30148503
Parks, D. H. et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat. Biotechnol. 36, 996–1004 (2018).30148503 10.1038/nbt.4229
69. Fu L Niu B Zhu Z Wu S Li W CD-HIT: accelerated for clustering the next-generation sequencing data Bioinformatics 2012 28 3150 3152 10.1093/bioinformatics/bts565 23060610
Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152 (2012).23060610 10.1093/bioinformatics/bts565
70. Buchfink B Reuter K Drost H-G Sensitive protein alignments at tree-of-life scale using DIAMOND Nat. Methods 2021 18 366 368 10.1038/s41592-021-01101-x 33828273
Buchfink, B., Reuter, K. & Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 18, 366–368 (2021).33828273 10.1038/s41592-021-01101-x
71. Ruscheweyh H-J Cultivation-independent genomes greatly expand taxonomic-profiling capabilities of mOTUs across various environments Microbiome 2022 10 212 10.1186/s40168-022-01410-z 36464731
Ruscheweyh, H.-J. et al. Cultivation-independent genomes greatly expand taxonomic-profiling capabilities of mOTUs across various environments. Microbiome 10, 212 (2022).36464731 10.1186/s40168-022-01410-z
72. Gallia R Passarelli A Conversano F 2024 Physical-oceanographic data (CTD and thermosalinograph) of NEREA Augmented Observatory Zenodo 10.5281/zenodo.10986789
Gallia, R., Passarelli, A. & Conversano, F. Physical-oceanographic data (CTD and thermosalinograph) of NEREA Augmented Observatory. Zenodo10.5281/zenodo.10986789 (2024).10.5281/zenodo.10986789
73. Margiotta F 2024 Biogeochemical data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11035857
Margiotta, F. et al. Biogeochemical data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11035857 (2024).10.5281/zenodo.11035857
74. Carotenuto Y d’Ippolito G 2024 LOFAs data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11058615
Carotenuto, Y. & d’Ippolito, G. LOFAs data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11058615 (2024).10.5281/zenodo.11058615
75. Saggiomo M Buondonno A 2024 HPLC data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11035381
Saggiomo, M. & Buondonno, A. HPLC data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11035381 (2024).10.5281/zenodo.11035381
76. D’Ambra I 2024 Stable isotope ratios data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.12743875
D’Ambra, I. Stable isotope ratios data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.12743875 (2024).10.5281/zenodo.12743875
77. Casotti R Balestra C Trano AC 2024 Flow cytometry data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11035933
Casotti, R., Balestra, C. & Trano, A. C. Flow cytometry data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11035933 (2024).10.5281/zenodo.11035933
78. Sarno D Percopo I 2024 Phytoplankton data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.10987241
Sarno, D. & Percopo, I. Phytoplankton data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.10987241 (2024).10.5281/zenodo.10987241
79. Percopo I Maselli M Sarno D 2024 Microzooplankton data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.10987215
Percopo, I., Maselli, M. & Sarno, D. Microzooplankton data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.10987215 (2024).10.5281/zenodo.10987215
80. Di Capua I 2024 Mesozooplankton data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.10987255
Di Capua, I. Mesozooplankton data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.10987255 (2024).10.5281/zenodo.10987255
81. Trano AC Casotti R Raffini F Piredda R 2024 16S Amplicon sequence variants (ASVs) data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.12801913
Trano, A. C., Casotti, R., Raffini, F. & Piredda, R. 16S Amplicon sequence variants (ASVs) data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.12801913 (2024).10.5281/zenodo.12801913
82. Montresor M Zampicinini G Piredda R Di Capua I 2024 18S Amplicon sequence variants (ASVs) data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.12801941
Montresor, M., Zampicinini, G., Piredda, R. & Di Capua, I. 18S Amplicon sequence variants (ASVs) data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.12801941 (2024).10.5281/zenodo.12801941
83. Campese L Ruscheweyh H-J 2024 Prokaryotic gene catalog, prokaryotic Metagenome-Assembled Genomes (MAGs) and taxonomic profiling of metagenomic data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11035656
Campese, L. & Ruscheweyh, H.-J. Prokaryotic gene catalog, prokaryotic Metagenome-Assembled Genomes (MAGs) and taxonomic profiling of metagenomic data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11035656 (2024).10.5281/zenodo.11035656
84. Campese L Ruscheweyh H-J 2024 Assemblies of metagenomic data of NEREA Augmented Observatory Zenodo 10.5281/zenodo.11046519
Campese, L. & Ruscheweyh, H.-J. Assemblies of metagenomic data of NEREA Augmented Observatory. Zenodo10.5281/zenodo.11046519 (2024).10.5281/zenodo.11046519
85. 2024 ENA European Nucleotide Archive PRJEB74649
ENA European Nucleotide Archive. https://identifiers.org/ena.embl:PRJEB74649 (2024).
86. 2024 ENA European Nucleotide Archive PRJEB74641
ENA European Nucleotide Archive. https://identifiers.org/ena.embl:PRJEB74641 (2024).
87. 2024 ENA European Nucleotide Archive PRJEB74658
ENA European Nucleotide Archive. https://identifiers.org/ena.embl:PRJEB74658 (2024).
88. 2024 ENA European Nucleotide Archive PRJEB74659
ENA European Nucleotide Archive. https://identifiers.org/ena.embl:PRJEB74659 (2024).
