Global Distribution and Diversity of Marine Parmales.

Parmales (Bolidophyceae) is a minor eukaryotic phytoplankton group, sister to diatoms, which exists as two distinct forms of unicellular organisms: silicified cells and naked flagellates. Since their discovery, many field studies on Parmales have been performed; however, their global distribution has not yet been examined in detail. We herein compiled more than 3,000 marine DNA metabarcoding datasets targeting the V4 region of the 18S rRNA gene from the EukBank database. By linking this large dataset with the latest morphological and genetic information, we provide updated estimates on the diversity and distribution of Parmales in the global ocean at a fine taxonomic resolution. Parmalean amplicon sequence variants (ASVs) were detected in nearly 90% of the samples analyzed. However, the relative abundance of parmaleans in the eukaryotic community was less than 0.2% on average, and the estimated true richness of parmalean ASVs was approximately 316 ASVs, confirming their low abundance and diversity. A phylogenetic ana-lysis divided these algae into four clades, and three known morphotypes of silicified cells were classified into three different clades. The abundance of Parmales is generally high in the poles and decreases towards the tropics, and individual clades/subclades show further distinctions in their distribution. Collectively, the present results suggest clade/subclade-specific adaptation to different ecological niches.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs.

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Improving taxonomic classification of marine zooplankton by molecular approach: registration of taxonomically verified 18S and 28S rRNA gene sequences.

Background: Zooplankton plays an important role in the marine ecosystem. A high level of taxonomic expertise is necessary for accurate species identification based on morphological characteristics. As an alternative method to morphological classification, we focused on a molecular approach using 18S and 28S ribosomal RNA (rRNA) gene sequences. This study investigates how the accuracy of species identification by metabarcoding improves when taxonomically verified sequences of dominant zooplankton species are added to the public database. The improvement was tested by using natural zooplankton samples. Methods: rRNA gene sequences were obtained from dominant zooplankton species from six sea areas around Japan and registered in the public database for improving the accuracy of taxonomic classifications. Two reference databases with and without newly registered sequences were created. Comparison of detected OTUs associated with single species between the two references was done using field-collected zooplankton samples from the Sea of Okhotsk for metabarcoding analysis to verify whether or not the newly registered sequences improved the accuracy of taxonomic classifications. Results: A total of 166 sequences in 96 species based on the 18S marker and 165 sequences in 95 species based on the 28S marker belonging to Arthropoda (mostly Copepoda) and Chaetognatha were registered in the public database. The newly registered sequences were mainly composed of small non-calanoid copepods, such as species belonging to Oithona and Oncaea. Based on the metabarcoding analysis of field samples, a total of 18 out of 92 OTUs were identified at the species level based on newly registered sequences in the data obtained by the 18S marker. Based on the 28S marker, 42 out of 89 OTUs were classified at the species level based on taxonomically verified sequences. Thanks to the newly registered sequences, the number of OTUs associated with a single species based on the 18S marker increased by 16% in total and by 10% per sample. Based on the 28S marker, the number of OTUs associated with a single species increased by 39% in total and by 15% per sample. The improved accuracy of species identification was confirmed by comparing different sequences obtained from the same species. The newly registered sequences had higher similarity values (mean >0.003) than the pre-existing sequences based on both rRNA genes. These OTUs were identified at the species level based on sequences not only present in the Sea of Okhotsk but also in other areas. Discussion: The results of the registration of new taxonomically verified sequences and the subsequent comparison of databases based on metabarcoding data of natural zooplankton samples clearly showed an increase in accuracy in species identification. Continuous registration of sequence data covering various environmental conditions is necessary for further improvement of metabarcoding analysis of zooplankton for monitoring marine ecosystems.

Metagenomic analysis provides functional insights into seasonal change of a non-cyanobacterial prokaryotic community in temperate coastal waters.

The taxonomic compositions of marine prokaryotic communities are known to follow seasonal cycles, but functional metagenomic insights into this seasonality is still limited. We analyzed a total of 22 metagenomes collected at 11 time points over a 14-month period from two sites in Sendai Bay, Japan to obtain seasonal snapshots of predicted functional profiles of the non-cyanobacterial prokaryotic community. Along with taxonomic composition, functional gene composition varied seasonally and was related to chlorophyll a concentration, water temperature, and salinity. Spring phytoplankton bloom stimulated increased abundances of putative genes that encode enzymes in amino acid metabolism pathways. Several groups of functional genes, including those related to signal transduction and cellular communication, increased in abundance during the mid- to post-bloom period, which seemed to be associated with a particle-attached lifestyle. Alternatively, genes in carbon metabolism pathways were generally more abundant in the low chlorophyll a period than the bloom period. These results indicate that changes in trophic condition associated with seasonal phytoplankton succession altered the community function of prokaryotes. Our findings on seasonal changes of predicted function provide fundamental information for future research on the mechanisms that shape marine microbial communities.

Development of a time-series shotgun metagenomics database for monitoring microbial communities at the Pacific coast of Japan.

Although numerous metagenome, amplicon sequencing-based studies have been conducted to date to characterize marine microbial communities, relatively few have employed full metagenome shotgun sequencing to obtain a broader picture of the functional features of these marine microbial communities. Moreover, most of these studies only performed sporadic sampling, which is insufficient to understand an ecosystem comprehensively. In this study, we regularly conducted seawater sampling along the northeastern Pacific coast of Japan between March 2012 and May 2016. We collected 213 seawater samples and prepared size-based fractions to generate 454 subsets of samples for shotgun metagenome sequencing and analysis. We also determined the sequences of 16S rRNA (n = 111) and 18S rRNA (n = 47) gene amplicons from smaller sample subsets. We thereafter developed the Ocean Monitoring Database for time-series metagenomic data ( http://marine-meta.healthscience.sci.waseda.ac.jp/omd/ ), which provides a three-dimensional bird's-eye view of the data. This database includes results of digital DNA chip analysis, a novel method for estimating ocean characteristics such as water temperature from metagenomic data. Furthermore, we developed a novel classification method that includes more information about viruses than that acquired using BLAST. We further report the discovery of a large number of previously overlooked (TAG)n repeat sequences in the genomes of marine microbes. We predict that the availability of this time-series database will lead to major discoveries in marine microbiome research.

Micro-size plankton abundance and assemblages in the western North Pacific Subtropical Gyre under microscopic observation.

While primary productivity in the oligotrophic North Pacific Subtropical Gyre (NPSG) is changing, the micro-size plankton community has not been evaluated in the last 4 decades, prompting a re-evaluation. We collected samples over three years (2016-2018) from depths of 10 to 200 m (n = 127), and the micro-size plankton were identified and counted to understand the heterogeneity of micro-size plankton community structure. The assemblages were consistent to the those of 4 decades ago. Dinophyceae (dinoflagellates) were the most numerically abundant, followed by Cryptophyceae and Bacillariophyceae (diatoms). The other micro-size plankton classes (Cyanophyceae, Haptophyceae, Dictyochophyceae, Euglenophyceae, and Prasinophyceae) were not always detected, whereas only Trichodesmium spp. was counted in the Cyanophyceae. Other unidentified autotrophic and heterotrophic flagellates were also significantly present, and their numeric abundance was higher than or at the same level as was that of the Dinophyceae. In the Dinophyceae, Gymnodiniaceae and Peridiniales were abundant. The chlorophyll a concentration and these class-level assemblages suggested micro-size plankton is not a major primary producer in this area. We applied generalized additive models (GAMs) and principal coordination analyses (PCoAs) to evaluate the habitats of every plankton group and the heterogeneity of the assemblages. The GAMs suggested that every classified plankton abundance showed a similar response to salinity, and we observed differences in habitats in terms of temperature and nitrate concentrations. Based on the PCoAs, we observed unique communities at the 200 m depth layer compared with those at the other sampling layers. The site scores of PCoAs indicated that the micro-size plankton assemblages are most heterogeneous at the 10 m depth layer. At such depth, diazotrophic Cyanophyceae (Trichodesmium spp.) are abundant, particularly in less-saline water. Therefore, nitrogen fixation may contribute to the heterogeneity in the abundance and assemblages in the western NPSG.

A DNA metabarcoding approach for recovering plankton communities from archived samples fixed in formalin.

Plankton samples have been routinely collected and preserved in formalin in many laboratories and museums for more than 100 years. Recently, attention has turned to use DNA information from formalin-fixed samples to examine changes in plankton diversity over time. However, no molecular ecological studies have evaluated the impact of formalin fixation on the genetic composition of the plankton community structure. Here, we developed a method for extracting DNA from archived formalin-preserved plankton samples to determine their community structure by a DNA metabarcoding approach. We found that a lysis solution consisting of borate-NaOH buffer (pH 11) with SDS and proteinase K effectively cleaved the cross-link formed by formalin fixation. DNA was extracted from samples preserved for decades in formalin, and the diatom community of the extracted DNA was in good agreement with the microscopy analysis. Furthermore, we stored a plankton sample for 1.5 years and demonstrated that 18S rRNA gene community structures did not change significantly from non-formalin-fixed, time-zero samples. These results indicate that our method can be used to describe the original community structure of plankton archived in formalin for years. Our approach will be useful for examining the long-term variation of plankton diversity by metabarcoding analysis of 18S rRNA gene community structure.

Draft Whole-Genome Sequence of Triparma laevis f. inornata (Parmales, Bolidophyceae), Isolated from the Oyashio Region, Western North Pacific Ocean.

We present the first draft whole-genome sequence for the Parmales (Bolidophyceae, Heterokonta), a picoplanktonic sister group of diatoms, using a Triparma laevis f. inornata strain that was isolated from the Oyashio region in the western North Pacific Ocean.

Depth-dependent transcriptomic response of diatoms during spring bloom in the western subarctic Pacific Ocean.

Diatoms play important roles in primary production and carbon transportation in various environments. Large-scale diatom bloom occurs worldwide; however, metabolic responses of diatoms to environmental conditions have been little studied. Here, we targeted the Oyashio region of the western subarctic Pacific where diatoms bloom every spring and investigated metabolic response of major diatoms to bloom formation by comparing metatranscriptomes between two depths corresponding to different bloom phases. Thalassiosira nordenskioeldii and Chaetoceros debilis are two commonly occurring species at the study site. The gene expression profile was drastically different between the surface (late decline phase of the bloom; 10 m depth) and the subsurface chlorophyll maximum (SCM, initial decline phase of the bloom; 30 m depth); in particular, both species had high expression of genes for nitrate uptake at the surface, but for ammonia uptake at the SCM. Our culture experiments using T. nordenskioeldii imitating the environmental conditions showed that gene expression for nitrate and ammonia transporters was induced by nitrate addition and active cell division, respectively. These results indicate that the requirement for different nitrogen compounds is a major determinant of diatom species responses during bloom maturing.

A Preliminary Metagenome Analysis Based on a Combination of Protein Domains.

Metagenomic data have mainly been addressed by showing the composition of organisms based on a small part of a well-examined genomic sequence, such as ribosomal RNA genes and mitochondrial DNAs. On the contrary, whole metagenomic data obtained by the shotgun sequence method have not often been fully analyzed through a homology search because the genomic data in databases for living organisms on earth are insufficient. In order to complement the results obtained through homology-search-based methods with shotgun metagenomes data, we focused on the composition of protein domains deduced from the sequences of genomes and metagenomes, and we utilized them in characterizing genomes and metagenomes, respectively. First, we compared the relationships based on similarities in the protein domain composition with the relationships based on sequence similarities. We searched for protein domains of 325 bacterial species produced using the Pfam database. Next, the correlation coefficients of protein domain compositions between every pair of bacteria were examined. Every pairwise genetic distance was also calculated from 16S rRNA or DNA gyrase subunit B. We compared the results of these methods and found a moderate correlation between them. Essentially, the same results were obtained when we used partial random 100 bp DNA sequences of the bacterial genomes, which simulated raw sequence data obtained from short-read next-generation sequences. Then, we applied the method for analyzing the actual environmental data obtained by shotgun sequencing. We found that the transition of the microbial phase occurred because the seasonal change in water temperature was shown by the method. These results showed the usability of the method in characterizing metagenomic data based on protein domain compositions.

Ontogenetic analysis of siliceous cell wall formation in Triparma laevis f. inornata (Parmales, Stramenopiles).

Triparma laevis f. inornata is a unicellular alga belonging to the Bolidophyceae, which is most closely related to diatoms. Like diatoms, T. laevis f. inornata has a siliceous cell wall. The cell wall of T. laevis f. inornata consists of four round plates (three shields and one ventral plate) and one dorsal and three girdle plates. But, unlike diatoms, T. laevis f. inornata cells can grow when concentrations of silica are depleted. We took advantage of this ability, using TEM to study the ontogeny of the siliceous plate, pattern center formation, and development. Two types of pattern centers (annulus and sternum) were observed in the early and middle stage of plate formation. During their formation, the annuli were initially crescent-shaped but eventually their ends fused to make a ring. Only outward silica deposition of the branching ribs occurred on the growing annulus until it became a ring, resulting in an unfilled circle inside the annulus. The pattern center of the shield plate was always an annulus, but in ventral plates both annulus and sternum were observed. The annuli and sterna in T. laevis f. inornata round plates were very similar to the annuli and sterna in diatom valves. These results suggested that the round plates of Parmales are homologous to diatom valves. This information on the plate ontogeny of T. laevis f. inornata provides new insights into the evolution of the siliceous cell wall in the Parmales and diatoms.

Mitotic spindle formation in Triparma laevis NIES-2565(Parmales, Heterokontophyta).

The parmalean algae possess a siliceous wall and represent the sister lineage of diatoms; they are thought to be a key group for understanding the evolution of diatoms. Diatoms possess well-characterized and unique mitotic structures, but the mitotic apparatus of Parmales is still unknown. We observed the microtubule (MT) array during interphase and mitosis in Triparma laevis using TEM. The interphase cells had four or five centrioles (∼80 nm in length), from which MTs emanated toward the cytoplasm. In prophase, the bundle of MTs arose at an extranuclear site. The position of centrioles with respect to an MT bundle changed during its elongation. Centrioles were observed on the lateral side of a shorter MT bundle (∼590 nm) and on either side of an extended MT bundle (∼700 nm). In metaphase, the spindle consisted of two types of MTs-MT bundle that passed through a cytoplasmic tunnel in the center of the nucleus and single MTs (possibly kinetochore MTs) that extended from the poles into the nucleus. The nuclear envelope disappeared only at the regions where the kinetochore MTs penetrated. In telophase, daughter chromosomes migrated toward opposite poles, and the MT bundle was observed between segregating chromosomes. These observations showed that MT nucleation does not always occur at the periphery of centrioles through cell cycle and that the spindle of T. laevis has a similar configuration to that of diatoms.

Seasonal and geographical distribution of near-surface small photosynthetic eukaryotes in the western North Pacific determined by pyrosequencing of 18S rDNA.

In this study, we investigated the distribution of small photosynthetic eukaryotes in the near-surface layer of the western North Pacific at four stations, including two oceanic stations where the subarctic Oyashio and subtropical Kuroshio currents influence a transition region and the bay mouth and head of the Sendai Bay, from April 2012 to May 2013. Flow cytometry was applied to sort small photosynthetic eukaryotes (<5 µm), and high-throughput sequencing of 18S rDNA was performed. Our taxonomic analysis showed that 19/195 operational taxonomic units (OTUs) were frequently distributed among all sites. Composition analysis showed that the OTUs had characteristic patterns and were divided into four main groups. Two groups reflected the low-saline water and winter season, with the characteristic OTUs belonging to diatoms; Chaetoceros and Leptocylindrus were characteristic of low saline water, and two diatom genera (Minidiscus and Minutocellus) and Cryptomonadales-related OTUs were prevalent in the winter. Our results indicate that the community composition of small photosynthetic eukaryotes seasonally changes in a dynamic manner according to variations in water properties.

Sequencing and analysis of the complete organellar genomes of Parmales, a closely related group to Bacillariophyta (diatoms).

We sequenced the complete plastid and mitochondrial genomes of the unicellular marine phytoplankton Triparma laevis, belonging to the order Parmales (Heterokonta). The cells of Parmales are surrounded by silicified cell walls, similar to Bacillariophyta (diatoms). T. laevis was recognized as a sister group of Bacillariophyta using a molecular phylogenetic analysis based on SSU rDNA and rbcL sequences. Bacillariophyta are the most successful group of phytoplankton in the modern ocean, but the origin and early evolution of them have not been clearly established. Detailed molecular analyses of T. laevis may increase our understanding of the evolutionary relationships among Parmales and Bacillariophyta. The gene contents of the plastid and mitochondrial genomes are similar between T. laevis and Bacillariophyta. The gene order of the plastid genome is also similar to Bacillariophyta, whereas the gene order of the mitochondrial genome is not conserved in Bacillariophyta, but the structure is more compact than Bacillariophyta. Phylogenetic analyses, using plastid-encoded concatenated amino acid datasets and mitochondria-encoded concatenated amino acid datasets suggest that T. laevis is a sister group of Bacillariophyta. These results suggest that the characteristics of the organellar genomes of T. laevis are similar and conserve ancestral characteristics more than Bacillariophyta.

Diversity and oceanic distribution of the Parmales (Bolidophyceae), a picoplanktonic group closely related to diatoms.

Bolidomonas is a genus of picoplanktonic flagellated algae that is closely related to diatoms. Triparma laevis, a species belonging to the Parmales, which are small cells with a siliceous covering, has been shown to form a monophyletic group with Bolidomonas. We isolated several novel strains of Bolidophyceae that have permitted further exploration of the diversity of this group using nuclear, plastidial and mitochondrial genes. The resulting phylogenetic data led us to formally emend the taxonomy of this group to include the Parmales within the Bolidophyceae, to combine Bolidomonas within Triparma and to define a novel species, Triparma eleuthera sp. nov. The global distribution of Bolidophyceae was then assessed using environmental sequences available in public databases, as well as a large 18S rRNA V9 metabarcode data set from the Tara Oceans expedition. Bolidophyceans appear ubiquitous throughout the sampled oceans but always constitute a minor component of the phytoplankton community, corresponding to at most ~4% of the metabarcodes from photosynthetic groups (excluding dinoflagellates). They are ~10 times more abundant in the small size fraction (0.8-5 µm) than in larger size fractions. T. eleuthera sp. nov. constitutes the most abundant and most widespread operational taxonomic unit (OTU) followed by T. pacifica, T. mediterranea and the T. laevis clade. The T. mediterranea OTU is characteristic of Mediterranean Sea surface waters and the T. laevis clade OTU is most prevalent in colder waters, in particular off Antarctica.

Spatial variation of bacterial community composition at the expiry of spring phytoplankton bloom in Sendai Bay, Japan.

In order to characterize how bacterial communities are propagated over spatial scales in a coastal area, the bacterial community composition was examined along with a transect line set in a bay at an expiry of spring phytoplankton bloom. Four distinctive bacterial communities were found within the bay by a fingerprinting method of 16S rRNA gene amplicons. The most widely distributed one was distributed in the surface and middle layers at whole area of the bay. The water was characterized by low inorganic nutrients concentration and high bacterial abundance, suggesting that the bacterial community had been developed in the bloom. Pyrosequencing analyses of the gene amplicons indicated that Rhodobacteriaceae and Flavobacteriaceae were abundant in the bacterial community, though the most abundant bacterial taxon was SAR11. The second group was distributed in the bottom water at the coastal side of the bay where considerably high Chl. a concentration was observed, probably because of the sedimentation of phytoplankton bloom. The community diversity was high and Alteromonadaceae, Saprospiraceae, and some families of Actinobacter existed more in this community than the others. The third group was distributed in the deep water near the border with the outside of the bay. The ratio of SAR11 was the highest in this community; besides, Burkholderianceae and Rhodospilliraceae existed in relatively high abundances. Another bacterial community having intermediate characters was observed in the middle to bottom layers around a central part of the bay where vertical water mixing was observed. These findings suggest that spatially different bacterial communities were formed under the influences of phytoplankton bloom and/or hydrographic events such as oceanic seawater intrusion of the bay.

Application of cryopreservation to genetic analyses of a photosynthetic picoeukaryote community.

Cryopreservation is useful for long-term maintenance of living strains in microbial culture collections. We applied this technique to environmental specimens from two monitoring sites at Sendai Bay, Japan and compared the microbial diversity of photosynthetic picoeukaryotes in samples before and after cryopreservation. Flow cytometry (FCM) showed no considerable differences between specimens. We used 2500 cells sorted with FCM for next-generation sequencing of 18S rRNA gene amplicons and after removing low-quality sequences obtained 10,088-37,454 reads. Cluster analysis and comparative correlation analysis of observed high-level operational taxonomic units indicated similarity between specimens before and after cryopreservation. The effects of cryopreservation on cells were assessed with representative culture strains, including fragile cryptophyte cells. We confirmed the usefulness of cryopreservation for genetic studies on environmental specimens, and found that small changes in FCM cytograms after cryopreservation may affect biodiversity estimation.

Effects of silicon-limitation on growth and morphology of Triparma laevis NIES-2565 (Parmales, Heterokontophyta).

The order Parmales (Heterokontophyta) is a group of small-sized unicellular marine phytoplankton, which is distributed widely from tropical to polar waters. The cells of Parmales are surrounded by a distinctive cell wall, which consists of several siliceous plates fitting edge to edge. Phylogenetic and morphological analyses suggest that Parmales is one of the key organisms for elucidating the evolutionary origin of Bacillariophyceae (diatoms), the most successful heterokontophyta. The effects of silicon-limitation on growth and morphogenesis of plates were studied using a strain of Triparma laevis NIES-2565, which was cultured for the first time in artificial sea water. The cells of T. laevis were surrounded by eight plates when grown with sufficient silicon. However, plate formation became incomplete when cells were cultured in a medium containing low silicate (ca. <10 µM). Cells finally lost almost all plates in a medium containing silicate concentrations lower than ca. 1 µM. However, silicon-limitation did not affect growth rate; cells continued to divide without changing their growth rate, even after all plates were lost. Loss of plates was reversible; when cells without plates were transferred to a medium containing sufficient silicate, regeneration of shield and ventral plates was followed by the formation of girdle and triradiate plates. The results indicate that the response to silicon-limitation of T. laevis is different from that of diatoms, where cell division becomes inhibited under such conditions.

ISOLATION AND CHARACTERIZATION OF PARMALES (HETEROKONTA/HETEROKONTOPHYTA/STRAMENOPILES) FROM THE OYASHIO REGION, WESTERN NORTH PACIFIC(1).

A small siliceous species of marine phytoplankton, order Parmales (Heterokonta), was isolated and characterized for the first time with the aid of a fluorescent silicon tracer 2-(4-pyridyl)-5-([4-(2-dimethylaminoethylaminocarbamoyl)-methoxy]phenyl)oxazole (PDMPO). This dye was easily detected by clear fluorescence in newly produced silica cell plates. Our isolate was surrounded by eight smooth plates without any ornamentation, suggesting a similarity to Triparma laevis B. C. Booth. TEM observation showed the typical ultrastructure of photosynthetic heterokontophytes; with two chloroplast endoplasmic reticulate membranes, a girdle lamella, three thylakoid lamellae, and mitochondrion with tubular cristae. Molecular phylogenetic analyses of SSU rDNA and rbcL genes showed that the parmalean alga was within the bolidophycean clade of autotrophic naked flagellates and a sister group of diatoms. HPLC analysis detected chl a, c1 + c2 , and c3 ; fucoxanthin; and diadinoxanthin as major photosynthetic pigments, and a composition that is shared with Bolidophyceae and diatoms. Together, these data indicate a close evolutionary relationship between Parmales, Bolidophyceae, and diatoms. The PDMPO-staining procedure should accelerate isolation of other Parmales species, helping to establish their diversity and aiding quantitative study of their role in oceanic processes.