A Portable Cruising Speed Net: Expanding Global Collection of Sea Surface Plankton Data

Plankton are central to planetary ecology, generating 50% of Earth’s atmospheric oxygen and forming the largest system of interconnected life at the base of the marine food chain. Yet, current oceanographic models aimed at predicting global climate change lack high-resolution biological data, emphasizing the need for innovative approaches to collect plankton biodiversity and distribution data over larger spatial, temporal, and taxonomic scales. The significant number of boats, ranging from small sailing yachts to large commercial vessels, that ply the world’s oceans every day could help scientists collect thousands of valuable plankton samples. Traditional Plankton Nets (TPN) are not suited to the speed of a recreational craft cruising in the high seas (i.e., at speeds >2 knots). We developed and validated the efficiency of a lightweight, easily deployable Cruising Speed Net (CSN) that enables the collection of ocean surface micro- and mesoplankton at speeds up to 5 knots. Field testing was conducted during two distinct research cruises along coastal and oceanic latitudinal gradients (SSV Robert C. Seamans in New Zealand and RV Investigator in the south-east Indian Ocean). DNA metabarcoding performed on the collected plankton samples showed the TPN and CSN yielded identical sequence-based diversity at low speed, with the CSN also effective at higher speed for characterizing latitudinal distribution of plankton communities. The CSN represents a valuable new tool for expanding the global collection of plankton data.


INTRODUCTION
Marine plankton play critical roles in marine ecosystems, being at the basis of oceanic food webs and responsible for generating at least half of the Earth's oxygen (Sekerci and Petrovskii, 2015). Plankton assemblages in the ocean vary spatially among biogeographic zones in much the same way as terrestrial biomes where different combinations of plant species occupy distinct land regions (Tappan and Loeblich, 1973). Heavily influenced by abiotic and biotic parameters of the water masses, plankton assemblages are also very sensitive to rapid changes in the ecology of the marine ecosystem related to global warming or other anthropogenic influences (Landry et al., 2020). Therefore, comprehensive plankton datasets are critical as indicators of ecological dynamics brought on by stressors, such as pollution or climate change (Hallegraeff, 2010).
Currently, oceanographic models aimed at predicting global ecological changes (Follows et al., 2007;Follows and Dutkiewicz, 2011;Ward et al., 2014) lack good quality, high resolution data on the nature and dynamics of marine plankton biodiversity. However, modeling the complexity of plankton assemblages is challenging due to small-and largescale spatiotemporal variations that do not necessarily follow classic ecological laws, a paradigm known as "the paradox of the plankton" (Hutchinson, 1961). Fine resolution, large-scale plankton data is therefore required for extrapolation (D'Alelio et al., 2019). Despite the availability of >60 years of extremely valuable and relevant data derived from projects such as the modern Continuous Plankton Recorder series (Reid et al., 2003(Reid et al., , 2010Richardson et al., 2006), there exists no complete global distribution map of plankton assemblages. Furthermore, modern comprehensive sampling of open ocean life is hindered by extremely expensive operational ship-time (typically more than 30,000 USD per day), with more traditional oceanographic research vessels having limited capacity to accommodate a broad range of research projects and/or sampling requirements (Lauro et al., 2014). These factors often influence research cruise logistics and result in minimized collection time with sampling only at specific stations (Sameoto et al., 2000). As a result, some areas of the ocean are almost devoid of plankton composition observations (taxonomic assemblages), including much of the South and tropical Pacific Oceans.
In this context, the significant number of boats, ranging from small sailing yachts to large commercial vessels, that ply the world's oceans every day represents an extraordinary opportunity for marine, and especially plankton-focused, research. There are approximately 120,000 active commercial ships (e.g., trading cargos, tourism cruises, fishing vessels) in the world fleet (Moser et al., 2015), and an estimated 10,000 yachts sailing the world's oceans daily (YOREPS, 2020). Most vessels follow longestablished routes dictated by predominant weather and trade winds (Cornell, 2014), providing a unique chance for dense spatiotemporal sampling of plankton. However, this concept is not new, and many studies have already exploited its potential (Lauro et al., 2014;Duarte, 2015;Simoniello et al., 2019). Most significant is the Ship of Opportunity Program run by the World Maritime Organization, which has deployed Continuous Plankton Recorders from large commercial vessels for more than 60 years; however, this collection technique is limited to a narrow range of plankton which are susceptible to capture, preservation and visual taxonomic identification in the towed recorders (Reid et al., 2010). The Global Ocean Sampling (GOS) expedition (Rusch et al., 2007;Williamson et al., 2008) paved the way for a new avenue of oceanography based on the use of sailing boats and advanced DNA-based ecology. This concept was expanded with the Tara Oceans expeditions (TO, 2009(TO, -2014, during which molecular, taxonomy based on morphological characteristics, and environmental data were collated from plankton assemblages (viruses to animals) across the world's oceans (Karsenti et al., 2011). These programs resulted in a visionary update to our knowledge of the extraordinary biodiversity and global structure of the marine planktonic ecosystem. However, developing a working understanding of ecological and evolutionary dynamics requires a far greater sampling effort.
One recent initiative, Plankton Planet 1 , equipped 27 international sailing yachts with a sampling toolkit enabling 1 www.planktonplanet.org citizen sailors to collect and isolate sea surface micro- (20-200 µm), meso-(200-2000 µm), and macroplankton (>2000 µm) using a Traditional Plankton Net (TPN, see example in Figure 1) while sailing the high seas. In a single year, this effort revealed a global picture of ocean eukaryotic plankton by using metabarcoding analyses on samples collected from 258 sites . However, a critical limitation to the Plankton Planet initiative was the operational requirement to reduce vessel speed to <2 knots to allow the 20 µm mesh net to tow stably within the water column. Such low speed made it very difficult for sailors to maneuver their vessels in the open sea, even precluding sampling in difficult weather conditions; this was identified as the primary limiting factor for achieving more extensive sampling at a more meaningful scientific scale. Thus, a modified net is needed -a low-cost, small size and effective plankton sampling device capable of collecting undamaged and unbiased surface ocean plankton assemblages at about 5 knots, which would not impact cruising speed and could be easily deployed from any vessel type, small sailing yachts to large commercial vessels.
In this study, we designed and built an easily deployable Cruising Speed Net (CSN) enabling the collection of ocean surface marine micro-and mesoplankton by tow from a vessel at a minimum speed of 5 knots. Two case studies were conducted for field-based testing of the plankton sampling devices and comparison of plankton biodiversity derived from metabarcoding analyses of nuclear small-subunit ribosomal DNA sequences. The first case study occurred on-board a brigantine sailing vessel (SSV Robert C. Seamans) dedicated to educational oceanographic research while cruising the eastern coast of New Zealand, and tested two hypotheses: (i) there will be no significant differences in plankton community composition captured by the CSN and TPN devices when used side-by-side at approximately 2 knots, and (ii) when increasing the vessel speed to approximately 5 knots using the CSN, the captured plankton community composition will remain the same. The second case study was conducted on-board a large oceanographic research vessel (RV Investigator) and tested the CSN at 5 knots along a 30 • latitudinal transect in the Indian Ocean (off Western Australia), evaluating the hypothesis: (iii) the CSN at 5 knots combined with metabarcoding analysis will be effective in capturing subtle plankton community composition changes associated with the latitudinal gradient. This study also offers complete specifications and general assembly guidelines for the CSN, with the hope that other researchers will use it and contribute to rapidly increasing our knowledge of the diversity, distribution, and evolution of open ocean plankton.

Description of Traditional Plankton Net (TPN) and Cruising Speed Net (CSN)
The TPN consisted of a conical 1.5 m long (20 µm nylon mesh; NITEX TM , Sefar Ltd., Switzerland) tow net with a 500 ml codend (including a 0.5 kg lead weight at the bottom but no mesh window) and a 30 cm diameter metal ring at the mouth, towed by three equally spaced bridle lines ( Figure 1A). The bridle lines were connected to a stainless-steel ring (3 cm diameter) to which three additional lines were also attached. The first line was ∼50 cm long with a 1 kg lead weight attached at the bottom, with the purpose of sinking the TPN slightly below the surface whilst being towed. The second line was 2 m long, leading to a buoy at the water surface and ensuring that the TPN remained at <2 m within the water column. The third line was a >20 m long (2 cm diameter) Dynice Dyneema rope used to tow the TPN from the vessel's rail. The very small nylon mesh size of the TPN combined with the necessary attachments (e.g., buoy and weight) created considerable drag, forcing a maximum tow speed of 2 knots among other problems such as extrusion of biota through the mesh and clogging. Vessel speeds of >3 knots invariably pulled the TPN to the surface, precluding plankton sampling.
CGW Consulting Engineers Ltd. (Nelson, New Zealand) assisted in designing a new and low-cost plankton sampling device able to easily collect plankton assemblages for microscopic and molecular identification, from any vessel type traveling at >3 knots. Design specifications for the sampling device were to: (i) be easily deployable at up to 5 knots, (ii) be relatively light and compact to minimize storage space use on small vessels, (iii) have a small opening of 10-20 cm diameter able to capture 2-5 m 3 of seawater in ca. 5 min through a 20 µm mesh size, and (iv) have a solid front head, with or without wings, shaped hydrodynamically to sink and maintain device stability at 2-3 m depth with minimal weights and lines to reduce drag and to facilitate ease of handling during sampling.
Three different candidate models were designed, built and field tested (Supplementary Figure 1). Upon preliminary field testing the first two models (Supplementary Figures 1A,B) were highly unstable in the water column and thus not considered for further improvement. The third model (Supplementary Figure 1C) met most of the requirements and was selected for additional field testing and fine-tuning. This CSN consisted of a specifically shaped PVC head, with a 20 cm wide isometric front eye or mouth of the net that was positioned at a 20 • angle, and a slightly conical 2 m long (20 µm nylon mesh; NITEX, Sefar Ltd., Switzerland) tow net with a 500 ml cod-end independently supplied by Whitlock Engineering Ltd., Nelson NZ ( Figure 1B and Supplementary Figure 2). Early field testing showed that the CSN performed well at low speed but became unstable at >4 knots, often resulting in the device flipping over and going backward, followed by rapid surfacing. This limitation was resolved by (i) slightly curving the two sides of the front eye in a concave shape (acting as an "underwater kite"), (ii) reducing the CSN attachments to two front bridle lines only, and (iii) adding a swivel to the main towing line ( Figure 1B). Further field testing demonstrated that the new modifications allowed the CSN to remain stable within the water column at towed speeds of up to 6.5 knots and it very effectively collected a variety of plankton species, as evidenced by microscopic observations 2 .

New Zealand Case Study
Sampling was conducted on-board a 41 m steel brigantine-rigged sailing vessel (SSV Robert C. Seamans) operated by the Sea Education Association (SEA; Woods Hole, MA, United States) for oceanographic research and education. Field testing occurred between 16 February 2019 and 20 March 2019 along the east coast of New Zealand's North and South Islands (35.0 • S, 174.01 • E to 42.8289 • S, 173.9039 • E) during SEA cruise S284. Samples were collected using the TPN and/or CSN as frequently as possible (i.e., at least once a day -time and weather permitting) following two models: In the first approach, employed in 21 distinct stations, both plankton net devices (TPN and CSN) were towed sideby-side, at 1-2 m below the sea surface from a metal boom extending off the vessel's port side at approximately 2 knots ship speed for 5 min (Table 1 and Figure 2). The purpose of sideby-side sampling was to compare the TPN and CSN's plankton capture efficiency, assessed through consistency in the recovered eDNA plankton diversity and community composition. The second sampling scheme involved the deployment of the CSN only, for the same duration, and aimed to assess the CSN's ability to capture representative plankton assemblages at higher sailing speeds (3-5.1 knots). On eight occasions, this higherspeed sampling was conducted immediately following the sideby-side sampling, hence considered the same stations (Table 1 and Figure 2). Before each new station, sterile gloves were worn and all gear (i.e., sampling nets, cod-ends, collection bottles, and filtration units) were thoroughly cleaned in local seawater mixed with a 2% bleach solution for at least 10 min, then rinsed with seawater from each sampling site to avoid cross contamination from other stations. Directly following each net tow, the codend was carefully removed and approximately 500 ml of the collected seawater was carefully poured into a sterile plastic bottle. Geographic location and time were recorded for each station at the start of each tow ( Table 1).

Australian Case Study
Sampling was conducted in the Indian Ocean on-board a 94 m oceanographic research vessel (RV Investigator, CSIRO, Hobart, Tasmania) between 15 May 2019 and 10 June 2019, along an approximately 30 • latitudinal transect off Western Australia (39.4910 • S, 109.9890 • E to 11.5040 • S, 110.0100 • E; Table 1 and Figure 2). For this case study, only the CSN was deployed, once per day in the afternoon and towed for 5 min at 5 knots following the same procedure as described above.

Filtration
The collected plankton material from TPN and/or CSN tow was carefully resuspended by shaking the plastic sample bottles, of which 200 ml was used for filtration thereby retaining microand mesoplankton samples for further analyses. Filtration was done using a sterile Whatman filter (pore size 3 µm; Merck KGaA, Darmstadt, Germany) placed in a Sterifil filtration system (Merck, Darmstadt, Germany) which used a 12V Seaflo 21 Series Water Pressure Pump 3.8LPM (Marine Deals, Auckland, New Zealand). Each filter, and its retained plankton material, was cut in half and stored into two separate 2 ml sterile cryotubes (one of the tubes with the second half of the cut filter was kept as a backup) containing 1.5 ml of LifeGuard Soil Preservation Solution (QIAGEN, Hilden, Germany). Samples were frozen at −20 • C immediately after filtration on the research vessels, transported to the Cawthron Institute (Nelson, New Zealand) post cruise and stored at −80 • C until further processing.

DNA Extraction
In the laboratory, filters were directly transferred into ZR BashingBead Lysis Tubes (2.0 mm; Zymo Research, CA, United States) containing 1 ml Lysis Buffer from the ZR Duet TM DNA/RNA MiniPrep Kit Plus (Zymo Research). All samples were homogenized via bead beating for 2 min (1600 MiniG Spex SamplePrep, NJ, United States), and centrifuged (10,000 × g, 5 min, 20 • C; Eppendorf Centrifuge 5430R, Hamburg, Germany). DNA was extracted from each sample following the manufacturer's protocol and DNA extraction blanks (i.e., negative controls) were included for each extraction series. The quantity and quality of extracted DNA were measured using a NanoPhotometer (Implen, Munich, Germany).

High-Throughput Sequencing
For the characterization of eukaryotic communities, an approximately 400 base pairs (bp) fragment of the V4 region of the nuclear small subunit ribosomal DNA (18S rRNA gene) was amplified from each sample by a single Polymerase Chain Reaction (PCR). The eukaryotic-specific primers were Uni18SF: 5 AGG GCA AKY CTG GTG CCA GC3 and Uni18SR: 5 GRC GGT ATC TRA TCG YCT T3 (Zhan et al., 2013). Illumina TM overhang adaptors were attached to the primers to allow dualindexing as described in Kozich et al. (2013). Thermocycling PCR conditions were: 94 • C for 3 min followed by 35 cycles of 94 • C for 30 s, 52 • C for 30 s, 72 • C for 1 min, with a final extension step at 72 • C for 7 min. Negative controls were included in each PCR analysis.
Amplicons were purified by magnetic separation following the Agencourt AMPure XP protocol (Agencourt Bioscience Corporation, Beverly, MA, United States), quantified (Qubit 2.0 Fluorometer, Invitrogen, Carlsbad, CA, United States) and diluted to 3 ng/µl. Two samples of 20 µl of RNA/DNA free water (UltraPure) were added as negative sequencing controls. Amplicons were sent to Auckland Genomics (The University of Auckland, Auckland, New Zealand) for library preparation. Sequencing adapters and sample specific indices were added to each amplicon using the Nextera Index kit (Illumina). Amplicons were pooled into a single library and paired-end sequences (2 × 250 bp) generated on an Illumina MiSeq using the TruSeq SBS kit (Illumina).

Bioinformatics (Both Case Studies)
Cutadapt version 2.10 (Martin, 2011) was used to remove the primer sequences from the raw reads with a single mismatch being allowed. These trimmed sequences were subsequently processed using the DADA2 package [version 1.12.1; Callahan et al. (2016) within R version 3.6.1; R Core Team (2017)]. The reads were truncated to 250 and 240 bp (forward and reverse reads, respectively) and filtered with a maximum number of "expected errors" (maxEE) threshold of two forward reads and two reverse reads. If reads did not meet this threshold, they were discarded from further analysis. A parametric error matrix was constructed within the program to dereplicate sequence variants for the forward and reverse reads. Paired-end reads were merged with a maximum mismatch of 1 bp and a minimum overlap of 10 bp and singletons were discarded. Chimeric sequences were removed within DADA2 using the consensus option in the removeBimeraDenovo function. The Amplicon Sequence Variants (ASVs) for the 18S rDNA gene were taxonomically classified against the PR2 (Guillou et al., 2013) database. The DADA2 assignTaxonomy function, based on the rdp classifier (Wang et al., 2007), was run with a bootstrap cutoff of 0.9. DNA negative controls were assessed and read numbers for ASVs found in the negative blanks were removed by proportional subtraction. The raw FASTQ sequence data were deposited in the Short Read Archive, BioProject PRJNA657626.

Statistical Analyses
Plankton biodiversity estimates were derived from ASVs, which may include both intra-and inter-organismal genetic variability. Consequently, this approach cannot accurately measure absolute abundances of plankton organisms, but rather estimate plankton composition based on the relative proportions of annotated plankton ASVs recovered from each environmental sample. Therefore, "plankton biodiversity" is hereafter referred to as "ASV-diversity." The New Zealand and Australian datasets were analyzed separately. Rarefaction curves of sequence reads from all samples were performed using ggrare in the R package ggplot2 (version 3.3.0). For community comparisons, untransformed numbers of reads of both datasets were rarefied down to 23,440 reads per sample using the rarefy_to_even_depth function in the phyloseq R package (version 3.6.1). Samples were fourth root transformed and Bray Curtis dissimilarity matrices were computed for each dataset, as advised in Clarke et al. (2017).  Table 1). Black circles, side-by-side sampling using the Traditional Plankton Net (TPN) and Cruising Speed Net (CSN) at 2 knots; white circles, additional and/or unique sampling at 3-5 knots using the CSN.

New Zealand Case Study
The New Zealand dataset was statistically analyzed based on two main factors: (i) TPN versus CSN (hereafter referred to as factor "net") and (ii) Samples collected at <3 knots versus samples collected at >3 knots (hereafter referred to as factor "speed"). Alpha-diversity indices (at ASV level) were calculated for the factors "net" and "speed, " including: Observed ASVs (S), Margalef richness (d) and Shannon diversity (H'), and displayed in boxplots. The effect of "net" and "speed" on these metrics was then examined using a pairwise t-test with 999 permutations. For beta-diversity analysis, low reads (<10) were removed from the dataset to avoid stochastic effects on the data. Significant differences between the plankton assemblages for the "net" and "speed" factors were tested using distancebased permutational analysis [PERMANOVA; Anderson and Epifanio (2009) implemented in Primer7 PRIMERE Ltd., United Kingdom; Clarke and Gorley (2015)] with "net" and "speed" as fixed factors ("speed" factor nested within "net").
Exploring differences in ASV-diversity and which taxa accounted for the differences between factors, a two-way crossed species contribution to sample (dis)similarities (SIMPER) analysis was performed based on presence/absence data, as described in San Martín et al. (2005), and splitting the New Zealand dataset into the factors "net" and "speed" with a cutoff percentage set at 70%. As there were no significant taxa associated with each subgroup of the factors (Supplementary  Table 2), the 30 most abundant ASVs of the fourth root transformed dataset, assigned at lowest possible taxonomic level, were visualized in a shade plot. A metric Multi Dimensional Scaling (mMDS) plot, based on fourth root transformation and Bray Curtis similarity matrix, was used to visualize similarities of samples for factors "net" and "speed."

Australian Case Study
The Australian dataset was visualized using Principal Coordinates Analysis (PCO) implemented in Primer7 (Clarke and Gorley, 2015), including the six most abundant plankton species (i.e., ASVs), assigned at lowest possible taxonomic level and overlaid as proportional pie charts.

RESULTS
The 76 plankton samples collected in this study (51 samples from the New Zealand and 25 samples from the Australian case studies) generated a total of 12,088,295 raw reads, which decreased to 5,982,714 reads (2,839 ASVs) post filtering (i.e., denoising, merging and chimera removal). Following the rarefaction analysis (Supplementary Figures 3, 4), each sample was rarefied to 23,440 sequence reads, which resulted in the removal of five low read count samples from the downstream analysis (2 CSN samples at station 21, CSN and TPN samples at stations 22, and CSN at station 34; Table 1). Overall, the taxonomic assignment revealed 74 orders, 133 families, 302 genera, and 379 species (Supplementary Table 1). Residual contamination was detected across all control samples (2,929 sequence reads representing 12 distinct ASVs) and these were removed from the actual dataset following the same procedure as described in Bell et al. (2018).

New Zealand Case Study
The comparison of TPN and CSN data from New Zealand revealed no significant differences between either the sampling devices (factor "net") or the vessel speeds (factor "speed") for any of the calculated diversity indices: Observed ASVs, Margalef richness and Shannon diversity (Figure 3 and Table 2A). PERMANOVA analysis resulted in no significant p-values for either "net" or "speed" factors (Table 2B).

TABLE 2 | (A) Pairwise tests for each calculated diversity index (Observed ASVs, Margalef richness, and Shannon diversity) for factors "net" and "speed." (B)
Permutational multivariate analysis of variance (PERMANOVA) based on Bray Curtis similarities for the "net" and "speed" factors.  Table 2). Therefore, average dissimilarity values, even when based on presence/absence transformed data, were too low to be considered for further analysis. For example, the highest listed ASV, ASV_23, contributed with just 0.6% average dissimilarity to differentiation between CSN and TPN. On the other hand, ASV_23 had similar values of 0.45% and 0.59% in the TPN and the CSN, respectively (Supplementary Table 2). Twelve ASVs accounted for a 10% cumulative difference between nets, with high average abundance across both. The same pattern was detected between the speed factor groups (Supplementary Table 2). Overall, no clear indicator taxa were found to be characteristic of just one factor group. Any differences between factors seemed mainly due to random appearances of organisms that sporadically appeared in a given sample within a group.

(Supplementary
This hypothesis was also supported by the shade (Figure 4) and mMDS (Figure 5) plots, indicating higher affinity of plankton community by sampling station than either the net or speed factors. Plotting the 30 most abundant planktonic ASVs amongst all New Zealand samples revealed zooplankton taxa with highest abundances (top of shade plot) were largely mirrored between factor groups (Figure 4). Gelatinous taxa such as the salp Thalia sp. and the tunicate Oikopleura dioica, but also copepods such as Temora turbinata, as well as two unidentified Hydrozoan and Echinodermata ASVs, were observed with higher abundances in early stations (19) for both TPN and CSN devices. In later stations (10-22), a community shift was apparent in both nets as evidenced by the highly abundant copepod species, Oithona similis and Paracalanus parvus. While a few taxa were restricted to a single sample (e.g., Euterpina acutifrons at station 18 and the ciliophoran Dictyocystidae at station 19), the majority of dominant taxa were recovered in multiple samples, regardless of nets or speed. Nevertheless, single stations revealed rare taxa, including a radiolarian species (Collozoum pelagicum), a salp species (Iasis cylindrica) and an unidentified cnidarian taxon, exclusively associated with the TPN device. Similarly, a ctenophoran taxon as well as another Collozoum sp. were only captured by the CSN device.
The mMDS analysis ( Figure 5) yielded a relatively strong clustering of planktonic assemblages from each sampling station, regardless of the "net" and "speed" factors (e.g., 17,19,and 20), which was consistent with the PERMANOVA results ( Table 2). There was also some evidence of regional separation along the latitudes and longitudes of the cruise track. Both geographic location and proximity to shore were influential; for example, stations 1, 3, 4, 5, and 6 reflect nearshore sampling in the Northland region of New Zealand while 7 and 9 cluster further away, potentially because they were collected at the continental shelf break. Stations around or south of the Coromandel Peninsula clustered on the left side of the mMDS plot (samples 9, 10, and 11). However, the 0.25 stress values represented a complicated fit of the samples' similarities in the two-dimensional space. Station 2 appeared as an outlier from its "Northland neighbors" and contained high abundances of the copepods Oithona similis and Paracalanus parvus that characterized later sampling sites 10-22 (Figure 4). Station 3, 5, and 12 exhibited high variability between sample subgroups, attributable to the isolated presence of rare taxa. The cnidarian Bougainvillia sp. appeared only in station 3's CSN sample while Thalia sp. was detected in only one tow at station 5 (CSN2) as well as with only the TPN at station 12 (Figure 4).

Australian Case Study
A PCO (Figure 6) demonstrated that the metabarcoding based plankton assemblages captured by the CSN were heavily influenced by the latitudinal gradient, which broadly corresponded with increasing sea surface temperature from 12 to 28 • C. Plankton composition exhibited successive assemblage changes from the southernmost station 24 (latitude 39.491 • S) to station 33 (latitude 29.001 • S), followed by a clustering of assemblages at the northernmost sampling stations 35-43 (latitudes 24.558 • -11.504 • S). Interestingly, samples from multiple stations at similar latitudes captured almost identical plankton assemblages despite being collected several days and 100's of kilometers apart (e.g., station 23/27/28 and stations 35 -46/46/47 and stations 41/43/44).  Five of the six most dominant plankton species in the dataset were copepods, and a clear shift in species distribution was noted along the latitudinal gradient (Figure 6). For example, two copepod species, Oithona similis and Paracalanus parvus were prevalent at the southernmost sampled latitudes in the Indian Ocean; these were also the two most dominant species at similar latitudes in the New Zealand dataset (Figure 1 and  Figure 4). Oithona similis was abundant between stations 24 and 28, Paracalanus parvus between stations 25 and 32, Clausocalanus furcatus between stations 26 and 43. Lower latitudes (from 26.007 • S) were dominated by two other copepod species (Delibus nudus and Acrocalanus gracilis) and an unidentified Hydrozoa.

DISCUSSION
In the present study, we designed a novel plankton sampling device that operates effectively at the average cruising speed (∼ 5 knots) of a standard sailing vessel. This innovative CSN exploits hydrodynamics rather than weight to remain stabilized below the sea surface, eliminating the unnecessary drag encountered with traditional nets.
Many distinct plankton sampling techniques have been developed over time (Tranter and Fraser, 1968). Some, such as Niskin bottles (Pennington et al., 2016) or pumps (Bishop et al., 1992), can effectively capture plankton assemblages, but these methods require the vessel to stop thereby limiting the sampling area that can be covered. Plankton nets are one of the oldest, simplest, effective and least expensive methods of sampling plankton (Gutkowska et al., 2012) and can be used for both morphological and molecular taxonomic analyses (Harvey et al., 2017;Schabacker et al., 2020). Thus, the more effective the net is in performance (e.g., sample capture and towing), as well as ease of handling, the greater the sampling opportunity and capacity.
The main advantage of towing a plankton net, horizontally or vertically, is that it considerably increases sampling volume and spatiotemporal coverage, e.g., towing a net at about 5 knots represents about 10 km per hour (Sameoto et al., 2000). A major limitation for open ocean towing is the size of the net mesh, usually between 5 and 80 µm for phytoplankton, 125 and 300 µm for mesozooplankton, and 1000 µm for micronekton (John et al., 2001;Castellani and Edwards, 2017). These mesh sizes can create considerable drag and typically require the vessel to not exceed a speed of approximately 2 knots, a situation that is often impossible to meet in marine areas with large waves and swell. Moreover, nets designed for higher towing speeds have reduced sampling capacity, as mesh sizes are increased (ca. > 300 µm) to reduce drag and handling becomes more cumbersome as weight and size of the net tends to increase, too. These constraints mean that existing high-speed nets can only be operated from large vessels able to tow such heavy equipment because they generally have sufficient storage capacity and more crew members.
The first high-speed plankton samplers were developed in the early 1900s, including the Hardy Continuous Plankton Recorder (Wiebe and Benfield, 2003). Off-the-shelf devices, such as the high-speed plankton collector "nackthai" (Schnack, 1974) perform well at up to 8 knots of vessel speed 3 . Another high-speed net example is the AVANI trawl; developed for collecting sea surface microplastics, it performs well in conditions of moderate sea state under 8 knots of boat speed over long transects (Eriksen et al., 2018). More recently, the Tara Ocean expedition developed a high-speed net able to collect macrozooplankton and neuston at 5-9 knots from the ocean atmosphere interface (Gorsky et al., 2019). All of these "high-speed" devices, however, present two major limitations. First, they are all equipped with nylon mesh sizes from >300 to 1000 µm, missing the vast majority of microand mesoplankton (<200 µm) taxa that are critical for holistic community assessments. In fact, the predicted 0.5-2.2 million plankton species do not even take into account the enormous diversity of small (20-200 µm) protist groups mainly associated with unknown taxa (De Vargas et al., 2015). Second, these devices are heavy and bulky. For example, the "nackthai" has a 150 cm long metallic frame with a 22 kg V-fin depressor while the 3 https://www.hydrobios.de/highspeedplanktoncollectornackthai/ Tara Ocean high-speed net consists of an 80 kg stainless steel head mounted with a 200 cm long nytrel net and cod-end. Furthermore, the high cost of these devices tends to preclude citizen science engagement in plankton collection.
The CSN represents a lightweight, affordable (estimated at ∼500 USD), compact and easily deployable alternative for the collection of ocean surface plankton (micro-and mesoplankton) data from a wider range of vessels and over much broader spatiotemporal scales than currently possible. We conducted two case studies to systematically test the performance and ability of the CSN to capture ocean surface plankton assemblages at two distinct cruising speeds and across two geographical gradients. The New Zealand case study showed that the CSN and TPN performed equally well at low (∼2 knots) speed, and that the same results were obtained when increasing the CSN towing speed to approximately 5 knots. Derived plankton composition was consistent between the two nets and speeds, with some variability likely explained by sporadic detection of rare and sparsely distributed taxa. No consistent bias toward over or underrepresentation of certain groups of organisms was noted. Some taxonomic outliers, such as different Collozoum species (Radiolaria), were detected in individual TPN or CSN samples which shows that single organisms randomly caught in the net can lead to some variability between devices due to sampling stochasticity but not device performances. Enhanced replication could have helped confirm our assumption that these taxa were incidental and do not indicate systematic taxonomic biases of one specific sampling approach. Indeed, other radiolarian taxa previously considered as rare in that area (Taylor, 1978), such as Acanthometron sp., appeared across all stations in equal abundance from both nets. Despite their important contribution to marine ecosystems (Countway et al., 2007;Edgcomb et al., 2011), radiolarians (including Collozoum sp.) are known to be difficult to capture under conventional plankton sampling due to their small size and fragile silica skeleton (Biard et al., 2017). In general, the sampling approach using both nets revealed geographic patterns between the northern and southern New Zealand stations, with the majority of TPN and CSN samples displaying similar community clustering. An overall shift was visible for gelatinous taxa (e.g., Thalia sp. and Oikopleura sp.) known to prefer the warmer, nutrient rich, coastal areas sampled along the Northland areas of New Zealand (Wiebe et al., 1979;Franco et al., 2016). Nearshore nutrient rich sampling stations also captured a high abundance of meroplanktonic forms, such as Chordata, Hydrozoa, and Echinodermata, as well as neritic copepods, including Temora turbinata, all consistent with previous observations (Bradford, 2010). Holoplanktonic copepod species, such as Oithona similis and Paracalanus parvus, dominated at sampling stations further offshore, confirming prior studies that reported higher abundance of these pelagic copepods in global oceanic waters Jose et al., 2017).
The Australian case study demonstrated that the CSN can be successfully operated from a large oceanographic research vessel traveling at a constant cruising speed of 5 knots. The device was very effective in capturing subtle biogeographic changes in plankton species distribution along a 30 • latitudinal gradient in the Indian Ocean, confirming data collected in the historic voyage of the International Indian Ocean Expedition (IIOE) in the 1960s by Tranter (1977). The most remarkable result was the very strong and linear successional shift in community distribution from the highest latitude in the South (station 24) to the lowest latitude in the North (station 43) of the transect (Figure 6), consistent with previous studies in the region (Eriksen et al., 2019). Also supporting published observations was the predominance and differential distribution of copepods species across latitudinal zones (Tranter, 1977;Schnack-Schiel et al., 2010). Interestingly, two oceanic copepod species (Oithona similis and Paracalanus parvus) were prevalent at similar (30 • -40 • S) latitudes in both case studies, substantiating the CSN's ability to reveal consistent geographically driven plankton diversity patterns and to reach a reliable comparability between different field trips. Plankton assemblages at lower latitudes were most significantly correlated with higher sea surface temperatures that favored gelatinous Hydrozoa (Wiebe et al., 1979) and tropical copepod taxa (Clausocalanus furcatus and Delibus nudus) Peralba et al., 2017). Finally, from local to global geographies, plankton are not homogenous but highly scattered and variable throughout the water column (Lima-Mendez et al., 2015). Many environmental variables such as salinity, oxygen and nutrients play a significant role in plankton distribution (Lynch et al., 2014;Landry et al., 2020), and should be incorporated in future CSN studies. We intentionally do not further discuss the fine-scale ecological patterns of plankton diversity reported in the present study, as they are beyond the scope of this paper. The aim of this second case study was to prove that the CSN device can be effectively used, from any vessel type able to reduce its cruising speed to 5 knots, to detect subtle shifts in sea surface planktonic assemblages and will boost upcoming studies seeking to identify driving environmental factors associated with current and future biodiversity patterns in open ocean environments.
Combining the CSN device with DNA metabarcoding to analyze plankton assemblages offers greater efficiency for taxa identification compared with conventional methods (e.g., microscopy), which can be time-consuming and require considerable taxonomic expertise and in spite of that, miss cryptic diversity and species devoid of conspicuous morphology. However, metabarcoding also presents some well acknowledged pitfalls, including possible biases due to primer or marker choice and the incompleteness of sequence reference databases (Harvey et al., 2017). The present study targeted the universal V4 region of 18S rRNA gene followed by sequence assignment using the PR2 database, enabling analysis of a broad range of taxonomic assemblages and including zoo-and phytoplanktonic taxa (De Vargas et al., 2015;Zamora-Terol et al., 2020). However, our data showed zooplankton (especially copepods) to be highly abundant compared with phytoplankton (e.g., Dinophyceae and Spirotrichea, see Supplementary Table 1) which were relatively scarce and mostly restricted to coastal sampling stations. Future studies aiming to capture a wider representation of phytoplankton taxa should consider using serial filtration methods before DNA extraction ) and/or include copepods-blocking oligos during PCR (Zamora-Terol et al., 2020). Additionally, simple DNA isolation methods should be considered while sampling plankton from sailing vessels with limited freezer (−20 • C) capacity. These may include the desiccation of filters or the use of DNA/RNA isolation buffers that can be kept cold (4 • C) for an extended period of time (Thomas et al., 2019;De Vargas et al., 2020). Finally, comparative morphological observations will also be critical for relating the sequence data to the plankton cell sizes and absolute abundances. Beyond purely scientific endeavors, it is also important to provide the general public with the widest range of accessible analytical toolkits to better understand the different "windows of the planktonic world" and improve the future management of our marine ecosystems. The recent development of a low-cost yet high-throughput PlanktoScope for citizen science applications  represents an exciting opportunity for future surveys.
Citizen Science (i.e., the involvement of volunteers who collect and/or process data as part of a scientific inquiry) has expanded considerably over the last 20 years across societies and biomes (Fore et al., 2001;Cohn, 2008;Newman et al., 2012). Citizen scientists now participate in various projects targeting various areas such as climate change, invasive species, conservation biology, ecological restoration, water quality monitoring, and population ecology (Silvertown, 2009). In the marine realm, the democratization of ocean observation has the potential to add a vast amount of observations every day. The rapidly decreasing costs of environmental DNA sequencing, imaging, sensing, and other analytical diagnostic tools (Cybulski et al., 2014;McLeod et al., 2015;Contreras-Naranjo et al., 2016;Srivathsan et al., 2019) creates new opportunities for researchers and citizen scientists to collect and analyze oceanographic data at unprecedented scales. This revolution has already started, with existing initiatives aiming to capitalize on the tens of thousands of sailing boats plying the world's oceans every day (Rusch et al., 2007;Williamson et al., 2008;Lauro et al., 2014;Simoniello et al., 2019;De Vargas et al., 2020). It is our ambition that the prototype CSN developed in this study will serve as an incentive for larger scale monitoring plankton surveys, empowering citizen science programs such as Plankton Planet and ultimately improving our knowledge on the diversity, distribution and evolution of open ocean plankton assemblages and the associated ecological processes. Future testing and improvements of the CSN are already underway and include (i) trialing the behavior and flow of the device when confronted with different particle sizes using Particle Imaging Velocimetry (PIV; Keane et al., 1994), (ii) visualizing the diversity and integrity of plankton cells captured by the CSN using the PlanktonScope , and (iii) working toward an even smaller and lighter version of the CSN similar to the Plankton Indicator device produced in the 1950s (Glover, 1953;Hardy, 1956;Planque and Reid, 2002).

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

AUTHOR CONTRIBUTIONS
XP, AZ, UA, and EM conceived and designed the experiments. DG and LB facilitated sampling cruises. UA, AJ, and AR performed the experiments. UA analyzed the data. UA, AJ, AZ, AR, DG, LB, EM, and XP contributed reagents, material, and tools and reviewed the manuscript drafts. UA, AZ, and XP wrote the manuscript. All authors contributed to the article and approved the submitted version.