Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Assessment of Zooplankton Community Composition along a Depth Profile in the Central Red Sea

Abstract

The composition of zooplankton in the water column has received limited attention in the main body of the Red Sea and this study investigates the change in the community both spatially and temporally across 11 stations in the central Red Sea. Using molecular methods to target the v9 region of the 18S rRNA gene a total of approximately 11.5 million reads were sequenced resulting in 2528 operational taxonomic units (OTUs) at 97% similarity. The phylum Arthropoda dominated in terms of reads accounting for on average 86.2% and 65.3% for neuston nets and vertical multinets respectively. A reduction in the number of OTUs was noticed with depth for both total metazoa and Maxillopoda whilst there was also a significant change in the composition of the Maxillopoda community. The genus Corycaeus had a higher proportion of reads in the epipelagic zone with Pleuromamma becoming increasingly dominant with depth. No significant difference was observed in the community between night and day sampling however there was a significant difference in the zooplankton community between two sampling periods separated by 10 days.

Introduction

Zooplankton play a substantial role in marine food webs and biogeochemical pathways [1]. Alterations in the planktonic community composition have been shown to have impacts higher up the trophic chain [2]. However whole community studies of zooplankton are limited and this is especially true of the Red Sea, a semi confined sea that combines high temperature, salinity and oligotrophic conditions [3]. The majority of studies undertaken on zooplankton in the Red Sea have concentrated on its northern extent and especially the Gulf of Aqaba [4, 5, 6]. These have included reports on zooplankton distributions throughout the water column [6, 7] as well as those at the surface [8]. However, studies in the central and southern reaches of the Red Sea are lacking. Böttger Schnack [9, 10] studied micro and mesozooplankton focusing on non-calanoid copepods whereas Schneider et al [11] described the abundance and size structure of the community. Whilst Pearman et al [12] studied zooplankton communities associated with coral reefs. The complex taxonomy of tropical communities and limited number of studies has meant that the overall diversity of the zooplankton assemblage has been poorly described.

Molecular techniques have enabled biodiversity studies to be undertaken at a greater depth and scale. Environmental DNA sampling relies on the extraction of bulk DNA and amplicon sequencing of a common molecular maker, such as the small subunit ribosomal RNA gene (SSU rRNA) or mitochondrial cytochrome oxidase I (COI). Molecular techniques were first utilized to study microbial populations (e.g. [13, 14]), but have more recently been applied to marine metazoan communities (for example [12, 15, 16, 17]. It is a technique, which allows analysis of organisms across a broad taxonomic range based on sequence similarities to reference databases. Molecular techniques have the possibility to simplify the analysis of marine zooplankton communities by allowing the whole community to be assessed simultaneously without the need for specific experts in morphological taxonomy of each taxa. As well as this the identification of larval and juvenile stages, which are complicated to differentiate morphologically would be able to be identified using genetic methods.

In this study, using high throughput sequencing methods we investigate the diversity of zooplankton communities in the central Red Sea region. To target the whole of the zooplankton community we used general eukaryotic primers targeting the v9 hypervariable region of the 18S rRNA gene [18]. Clustering of Operational Taxonomic Units (OTUs) was undertaken at 97% similarity so that OTU richness was most closely aligned to species richness [16, 19]. To date there have been few studies that have investigated the whole zooplankton community in the main body of the Red Sea and therefore this study provides a basic understanding of the community structure and improves our knowledge of the marine ecosystem in the Red Sea.

Methods

Sampling

Samples were collected during March—April 2013 in the central Red Sea region (S1 Table; Fig 1) during the KAUST Red Sea Research Cruise 2013 aboard the RV Aegaeo. Permits for sampling in Saudi Arabian waters were obtained from the Saudi Arabian coastguard. No specific permissions were required, as the study did not involve endangered or protected species. A total of 11 stations were sampled in the region with the stations being repeated 10 days after the first sampling. Zooplankton samples were collected using two methods. Firstly depth profiles were obtained using a multinet system (Hydro-Bios, Kiel) fitted with a mesh size of 64 μm and a CTD profiler attached. Vertical hauls were carried out between 500 m and the surface with a towing speed of 1 m s-1. The depths at which nets were triggered are detailed in S1 Table. Secondly, surface zooplankton samples were obtained using a neuston net with a 64 μm mesh size towed at approximately 3 knots for 10 min. Subsequent to sampling the samples were stored at 4°C in absolute ethanol.

thumbnail
Fig 1. Sampling stations.

Positions of sampling stations showing a) Positions relative to the Red Sea and b) the position of individual stations. Figure plotted in R using the package maps [20].

https://doi.org/10.1371/journal.pone.0133487.g001

DNA extraction and sequencing

DNA was extracted based on the methods described in [21]. Ethanol was removed from the samples prior to the addition of ATL lysis buffer (Qiagen). The whole sample was transferred to a pestle and mortar and crushed until the biomass was liquidized. To complete the lysis of the sample the resulting cell suspension was incubated overnight at 55°C with 20 μl proteinase K (20 mg/mL). For samples with larger biomass the addition of proteinase K was scaled up. Phenol:chloroform:iso amyl alcohol (IAA) was added in equal volume to the sample, mixed and centrifuged for 5 min. The aqueous layer was removed and a second round of phenol:chloroform:IAA extraction was undertaken followed by a round of chloroform:IAA. The aqueous layer was removed and DNA precipitated using 2.5 volumes ethanol and 0.1 volumes sodium acetate at – 20°C. DNA was washed in 70% ethanol and resuspended in autoclaved MilliQ water.

Amplification of eukaryotic DNA was achieved using the general eukaryotic primers designed by [18]. The forward primer 1389F was combined with the reverse primer 1510R. All primers had an Iontorrent tag (instead of a Roche 454 tag as used by [18]) and a 10 bp barcode. PCRs were undertaken in duplicate in 50 μl reaction volumes containing 2.5U Taq polymerase (Invitrogen), 1XTaq reaction Buffer, 200μM dNTPs (Invitrogen), 1.5mM MgSO4, 0.05 mg Bovine Serum Albumin (BSA) and 0.2μM of each primer. PCR conditions were adapted from [18]. An initial denaturation step of 3 min at 94°C was followed by 30 cycles of 94°C for 30 sec, 57°C for 45 sec and 72°C for 45 sec. This is followed by a final extension at 72°C for 10 min. A no template negative control was also run. Duplicate PCR reactions were combined and the PCR band visualized on a 1% agarose gel before being excised and purified with Qiagen’s gel extraction kit. Purified samples were pooled, having been equalized to 50ng (measured using a dsDNA high sensitivity kit for a Qubit 2.0 Fluorometer (Life Technologies)), according to barcode. A maximum of 15 unique barcodes were pooled together and multiplexed samples were analyzed for purity using an Agilent 2100 bioanalyzer machine. Emulsion PCR for IonTorrent sequencing was undertaken, using an Ion PGM Template OT2 400 Kit on a Ion One Touch system 2.0 following manufacturers protocols, with 15 pM of template amplicons. After enrichment in a Ion one touch ES system the templates were loaded on a 318 V2 chip and sequenced at the King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia, core facility on an IonTorrent PGM machine using an Ion PGM Sequencing 400 Kit. Raw sequence reads were deposited in the NCBI Sequence Read Archive under the experiment accession number: SRP056532.

Sequence analysis and bioinformatics

Raw sequences were automatically demultiplexed, based on barcode, on the IonTorrent PGM machine and each sample was filtered individually using the QIIME (version 1.8) software package [22]. Parameters for quality filtering were minlength = 100, maxlength = 200, quality = 25, windowsize = 25, max homopolymers = 6 and max forward primer mismatch = 1, barcode = 0. Filtered sequences were concatenated and clustering into OTUs occurred through a two-step process. Firstly, using the trie function of the QIIME version of the CD-HIT [23] algorithm, sequences were preclusted and the longest representative sequence of each precluster was selected. Secondly de novo clustering of these representative sequences using the USEARCH (version 5.2.236) [24] algorithm at 97% identity was undertaken where clusters had to include at least 2 sequences. Representative sequences (longest) were selected for the USEARCH OTUs and checked for chimeras against a reference database (SILVA 119 [25] using UCHIME [26]). Taxonomic assignments were made in QIIME against the SILVA 119 database using the naïve Bayesian classifier rdp [27]. To assess the composition across different stations reads were rarefied multiple times (n = 100) at an even depth (20,000 reads) and the average taken. Diversity statistics were undertaken using QIIME and the R (v 3.0.2 (R development core team)) packages vegan [28] and phyloseq [29]. Weighted and unweighted UniFrac [30] distance matrices were constructed in phyloseq based on the QIIME OTU table and phylogenetic tree. Non-metric multidimensional scaling (NMDS) plots were constructed with phyloseq and statistical analysis of the distance matrices was calculated using PERMANOVA (adonis in vegan).

Results

A total of 11,693,628 reads were obtained across all the samples and clustering resulted in a total of 2528 metazoan OTUs at 97% similarity, subsequent to the removal of 110 OTUs which were designated as chimeras as well as those which were not taxonomically classified as metazoa. In order for samples to be compared the metazoan reads were rarefied multiple times (n = 100) at an even depth (20,000 reads). After rarefaction a total of 2079 OTUs were observed. The majority of these OTUs were grouped in the phylum Arthropoda (1285 OTUs; 62% of total OTUs), and particularly the class Maxillopoda (886 OTUs; 43% of total OTUs). Other phyla accounting for more than 100 OTUs include the Cnidaria (131 OTUs; 6.3% of total OTUs), Urochordata (118 OTUs; 5.7% of total OTUs) and Mollusca (104 OTUs; 5.0% of total OTUs). Comparing alpha diversity down the depth profile it was found that there was a significant decrease in the numbers of both metazoan (F = 9.485; p = 0.004) and Maxillopoda (F = 15.9; p < 0.001) OTUs with depth (Fig 2).

thumbnail
Fig 2. Association of OTU counts with depth.

Number of OTUs for both Metazoa and Maxillopoda against depth. The average depth at which each multinet was open was used.

https://doi.org/10.1371/journal.pone.0133487.g002

At all stations, samples taken using the neuston net were dominated by the phyla Arthropoda (86.2%) (Fig 3A and 3B). Smaller proportions of reads were attributed to Chaetognatha (average 4.9%), Cnidaria (average 2.2%), Urochordata (average 1.7%) and Mollusca (average 1.6%). Mollusca and Cnidaria were more prevalent in the second sampling period (Fig 3B). On a finer taxonomic scale, the neuston nets were dominated by the class Maxillopoda (phyla Arthropoda) (average 76.9%), with smaller contributions to the samples from Malacostraca (phyla Arthropoda) (average 8.2%). During the second sampling period the Maxillopoda genus Corycaeus was more prevalent than in the first, accounting for an average of 10.23% and 3.24% of the community respectively.

thumbnail
Fig 3. Zooplankton composition for a) neuston and b) neuston repeat.

Composition of zooplankton at the phyla level for a) the neuston surface samples and b) the neuston surface repeated samples. Arthro = Arthropoda; Chaeto = Chaetognatha. Other accounts for those phyla, which accounted for less than 3% of a sample or OTUs that were not assigned to the phyla level. For station designations refer to S1 Table. Figure plotted in R using the package maps [20].

https://doi.org/10.1371/journal.pone.0133487.g003

On average throughout the water column Arthropoda were the dominant component of the plankton averaging 65.3%, however Cnidaria also accounted for a substantial proportion (average 15.1%) whilst Urochordata also accounted for a sizeable proportion of reads in several stations (Fig 4). Table 1 shows the taxonomic designations (at genus level where possible) of the most dominant OTU for the metazoan phyla down the water column. Maxillopoda were the dominant class of Arthropoda throughout the depth profile, however the genera present differed through the water column. A higher proportion of reads were attributed to the genus Pleuromamma at deeper depths than at shallower ones whilst the genus Corycaeus was more prominent in the top 100m (Fig 5).

thumbnail
Fig 4. Metazoa composition in multinet samples.

Combined composition of metazoa at each station in the multinet samples. Arthro = Arthropoda; Urocho = Urochordata. Other accounts for those phyla, which accounted for less than 3% of a sample or OTUs, which were not classified at the phyla level. For station designations refer to S1 Table. Figure plotted in R using the package maps [20].

https://doi.org/10.1371/journal.pone.0133487.g004

thumbnail
Fig 5. Depth profile of Maxillopoda taxa.

Depth profile of the composition of the Maxillopoda genera in the multinet samples. Pies are plotted at the average depth for which the net was open except where stations are close in longitude were one station is plotted at the 1st quartile whilst the other is plotted at the 3rd quartile to avoid overlap. Background is temperature (°C). Others includes those genera with < 1% reads in at least one sample and those OTUs classified as Maxillopoda but which could not be assigned to genus level. Figure plotted with R package oce [31].

https://doi.org/10.1371/journal.pone.0133487.g005

thumbnail
Table 1. Genus assignments for the dominant OTUs.

Taxonomic assignation of the dominant genus (where identified) for the metazoa phyla at each sampling point. See S1 Table and Fig 1 for station designations. Blank designations are where the dominant OTU could not be assigned to the genus level.

https://doi.org/10.1371/journal.pone.0133487.t001

To assess statistical differences in the zooplankton community samples, those, which had a full depth profile (3 night stations and 3 day stations), were analyzed. Adonis analysis was undertaken to investigate the effects of depth and time of day. Both weighted and unweighted UniFrac results showed that there was a significant difference with depth (weighted: F = 3.17, p < 0.001; unweighted F = 2.588, p < 0.001) with Fig 6A + 6B showing a NMDS of all sampling points. There was no significant effect of time of day (weighted: F = 2.018, p = 0.057; unweighted F = 1.514, p = 0.058). A similar finding was found when just considering the dominant group Maxillopoda with depth being found to be significant (weighted: F = 12.353, p < 0.001; unweighted F = 3.084, p < 0.001) whilst time of day was not significant (weighted: F = 0.782, p = 0.428; unweighted F = 1.197, p = 0.262).

thumbnail
Fig 6. NMDS plot for a) weighted and b) unweighted UniFrac for multinet samples.

NMDS plot based on a) weighted and b) unweighted UniFrac distance matrices showing clustering by depth. Samples were categorized based on the average depth, which the net was opened.

https://doi.org/10.1371/journal.pone.0133487.g006

In the neuston nets there was also found to be no significant difference in the metazoan community depending on the time of day of sampling (after 9 samples were randomly chosen for both night and day) (weighted: F = 0.789, p = 0.579; unweighted: F = 0.886, p = 0.612). However a statistical difference was observed between the two sampling periods (weighted: F = 2.738, p = 0.014; unweighted: F = 2.289, p = 0.004) between the first 10 stations (Station 11 only had results from the repeat sampling period) (Fig 7A + 7B). No significant trends (p > 0.05) were observed between distance and similarity of the zooplankton community for either UniFrac distance measure or sampling period.

thumbnail
Fig 7. NMDS plot for a) weighted and b) unweighted UniFrac for neuston samples.

NMDS plot based on a) weighted and b) unweighted UniFrac distance matrices showing clustering of the sampling period. A is the sampling period between 25th March 2013 and 27th March 2013 whilst B is the sampling period between 3rd April 2013 and 5th April 2013.

https://doi.org/10.1371/journal.pone.0133487.g007

Discussion

The present study uses high throughput sequencing technology to provide basic information on the vertical profile of zooplankton in the central Red Sea. Morphological investigations often focus on the abundance of each species in the zooplankton community. However using next generation sequencing techniques the correlation between reads and abundance cannot be presumed. Previous genetic studies of the zooplankton community targeting the 18S rRNA gene showed that for bivalve and decapod larvae metagenetic reads correlated more with morphological biomass than abundance [16], whilst [15] found a correlation between copepod LSU reads and dry weight. Although this correlation could be affected by biases such as primer mismatches or gene copy numbers [32, 33] it can give an approximate representation of changes in the biomass of the zooplankton community. The dominant phyla in the current study were Arthropoda in both the neuston and vertical profiles with Cnidaria also accounting for a substantial proportion of reads (Figs 3 and 4). Urochordata also accounted for a large proportion of reads in the vertical nets with Chaetognatha and Mollusca observed in the neuston nets. Similar groups were found to be the most abundant in previous studies in the Red Sea (e.g. [4, 5, 34, 35]) and in other water bodies including the subtropical Pacific [36] and UK coastal waters [16]. Copepods were the dominant component of the Arthropoda in both the neuston and vertical net samples. This has been previously shown in a molecular study of zooplankton associated with coral reefs in the Red Sea [12], morphological studies undertaken in the Gulf of Aden and Red Sea (e.g. [6, 11, 35, 37]) and in other marine basins (e.g. Mediterranean: [38]; subtropical Pacific: [36]; Caribbean: [39]). The dominant genus identified in the epipelagic zone was the Poecilostomatoid copepod Corycaeus, which has also been found to be present in other Red Sea studies [37] although other genera such as Clausocalanus and Ctenocalanus [7, 40] are often more abundant. Lindeque [16] has previously reported that Clauso-/Ctenocalanus accounted for a higher proportion of morphologically identified samples compared with molecular techniques. The differences were attributed to high abundance relative to biomass of juveniles of these species and similar explanations could account for the low proportion of these genera in the current study with an increase in copepodit stages of Clausocalanus spp. having been shown in March/April in the Gulf of Aqaba [41]. The calanoid Pleuromamma was dominant in the mesopelagic zone, which is in agreement with previous studies [7, 40] (Fig 5).

In the current study, in agreement with previous studies [6] we found that there was no significant difference in both the metazoan and maxillopoda community between samples taken at night and those taken during the day. However the dominant Maxillopoda genus Pleuromamma was observed to have higher abundances of reads in the deeper samples, as has been observed previously in the Red Sea [40]. In agreement with the fact this genus undertakes diel vertical migrations [40, 42, 43] higher abundances were observed in the epipelagic zone of those stations where sampling occurred during the night. However, large proportions of reads were still observed in the mesopelagic zone suggesting either a proportion of the population was non-migratory as has been observed in various species of this genus [43] or different life cycle distribution patterns [44].

A significant decrease in the number of metazoan and Maxillopoda OTUs was noticed with depth (Fig 2), which is in agreement with previous investigations in subtropical regions [6, 40, 45, 46] but shows a different distribution from that observed for zooplankton in the Amundsen and Makarov Basins where the maximum was observed between 500–750m [47].

Whilst studies have shown changes in the abundance of zooplankton in the Red Sea on a seasonal scale [4, 5, 7, 37] we are able to show the dynamic nature of the planktonic community with statistically significant changes in the neuston community between sampling of the same stations 10 days apart. The difference in the community appeared to be due to the change in rare taxa with the unweighted UniFrac, which favors rare taxa being more significant than the weighted UniFrac. However, the genus Corycaeus accounted for 10.23% of the reads in the repeat sampling and only 3.24% in the first period of sampling. This suggests that studies incorporating short time periods between sampling efforts will be required to fully understand the changes in the plankton community in the Red Sea.

One advantage of molecular techniques over those based on morphology is that larval and juvenile specimens can be hard to identify morphologically [16] whilst there is no differentiation in life stage when using molecular techniques. El Sherbiny et al [5] found that adult copepods only constituted 22.3% of the total copepods meaning identification of a substantial proportion of the zooplankton community would be problematic. The focus on adult copepods is a fundamental issue in morphological studies and could mean that there are misrepresentations of zooplankton diversity with [16] suggesting that next generation sequencing technologies may reveal double the copepod diversity.

Whilst sequencing technologies have their advantages there are currently limitations for their use in ecological studies. For example, in the current study, on average 40% of the reads from the multinet samples classified as Maxillopoda could not be classified at a lower taxonomic level. Genera such as Oncaea, which have previously been shown to be a key component of the Red Sea zooplankton community [5], are absent from reference databases meaning that reads would fall into the unknown category. Lindeque et al [16] found that off the UK coast only half the possible local species of Oithona were present in the database and this limits the ability to classify species at low taxonomic levels and missing sequences of representative species from the Red Sea may partly account for the lower proportion of this genera observed in the current study despite high abundances previously reported [35]. This problem is likely to be exacerbated in regions, such as the Red Sea where the diversity of zooplankton has received considerably less attention.

In summary, this study has increased the amount of knowledge available on the diversity of zooplankton in the central Red Sea. Whilst there are still limitations in the use of molecular techniques to study the diversity of metazoans this study has shown that the techniques are able to discriminate changes in the structure of the zooplankton community not only in time but also down a depth profile at different taxonomic levels. In the future with the collaboration of molecular ecologists and morphological taxonomists the gaps in reference databases can be closed leading to an improvement in the speed and scale at which ecological studies can be undertaken.

Supporting Information

S1 Table. Ancillary data for neuston and multinet samples including temperature, conductivity, salinity, depth, latitude, longitude, day/night and date of sampling.

https://doi.org/10.1371/journal.pone.0133487.s001

(CSV)

Acknowledgments

The authors would like to thank all those who helped in the process of undertaking this study including the captain and crew of the RV Aegaeo. We would also like to thank Ioannis Georgakakis for helping with the nets and CTD rosette used during the sampling process. Also we would like to thank the two reviewers for their comments on the manuscript. The authors declare no conflicts of interest.

Author Contributions

Conceived and designed the experiments: JKP XI. Performed the experiments: JKP. Analyzed the data: JKP. Contributed reagents/materials/analysis tools: JKP XI. Wrote the paper: JKP XI.

References

  1. 1. Mitra A, Castellani C, Gentleman WC, Jonasdottir SH, Flynn KJ, Bode A, et al. Bridging the gap between marine biogeochemical and fisheries sciences; configuring the zooplankton link. Progress in Oceanography. 2014; 129: 176–99.
  2. 2. Beaugrand G, Brander KM, Lindley JA, Souissi S, Reid PC. Plankton effect on cod recruitment in the North Sea. Nature. 2003; 426(6967): 661–4 pmid:14668864
  3. 3. Raitsos DE, Pradhan Y, Brewin RJW, Stenchikov G, Hoteit I. Remote Sensing the Phytoplankton Seasonal Succession of the Red Sea. PLoS ONE. 2013; 8(6): e64909. pmid:23755161
  4. 4. Dorgham MM, Elsherbiny MM, Hanafi MH. Vertical distribution of zooplankton in the epipelagic zone off Sharm El-Sheikh, Red Sea, Egypt. Oceanologia. 2012; 54(3): 473–89.
  5. 5. El-Sherbiny MMH, Mahmoud H.; Aamer , Maher A. Monthly variations in abundance and species composition of the epipelagic zooplankton off Sharm El-Sheik, Northern Red Sea. Research Journal of Environmental Sciences. 2007; 1(5): 200–10.
  6. 6. Cornils A, Schnack-Schiel SB, Hagen W, Dowidar M, Stambler N, Plahn O, et al. Spatial and temporal distribution of mesozooplankton in the Gulf of Aqaba and the northern Red Sea in February/March 1999. Journal of Plankton Research. 2005; 27(6): 505–18.
  7. 7. Cornils A, Schnack-Schiel SB, Al-Najjar T, Badran MI, Rasheed M, Manasreh R, et al. The seasonal cycle of the epipelagic mesozooplankton in the northern Gulf of Aqaba (Red Sea). Journal of Marine Systems. 2007; 68(1–2): 278–92.
  8. 8. Al-Najjar TH, El-Sherbiny MM. Spatial and Seasonal Variations in Biomass and Size Structure of Zooplankton in Coastal Waters of the Gulf of Aqaba. Jordan Journal of Biological Sciences. 2008; 1(2): 55–9.
  9. 9. Böttger-Schnack R. Summer distribution of micro- and small mesozooplankton in the Red Sea and Gulf of Aden, with special reference to non-calanoid copepods. Marine Ecology Progress Series. 1995; 118(1–3): 81–102.
  10. 10. Böttger-Schnack R. The microcopepod fauna in the Eastern Mediterranean and Arabian Seas: A comparison with the Red Sea fauna. Hydrobiologia. 1994; 292–293: 271–82.
  11. 11. Schneider G, Lenz J, Rolke M. Zooplankton standing stock and community size structure within the epipelagic zone–A comparison between the central Red Sea and the Gulf of Aden. Marine Biology. 1994; 119(2): 191–8.
  12. 12. Pearman J, El-Sherbiny MM, Lanzén A, Al-Aidaroos AM, Irigoien X. Zooplankton diversity across three Red Sea reefs using pyrosequencing. Frontiers in Marine Science. 2014;1:27:
  13. 13. Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, et al. Microbial diversity in the deep sea and the underexplored "rare biosphere". Proceedings of the National Academy of Sciences of the United States of America. 2006; 103(32): 12115–20. pmid:16880384
  14. 14. Logares R, Audic S, Santini S, Pernice MC, de Vargas C, Massana R. Diversity patterns and activity of uncultured marine heterotrophic flagellates unveiled with pyrosequencing. ISME Journal. 2012; 6(10): 1823–33. pmid:22534609
  15. 15. Hirai J, Kuriyama M, Ichikawa T, Hidaka K, Tsuda A. A metagenetic approach for revealing community structure of marine planktonic copepods. Molecular ecology resources. 2015; 15(1): 68–80. pmid:24943089
  16. 16. Lindeque PK, Parry HE, Harmer RA, Somerfield PJ, Atkinson A. Next Generation Sequencing Reveals the Hidden Diversity of Zooplankton Assemblages. PLoS ONE. 2013;8(11). e81327: pmid:24244737
  17. 17. Thomsen PF, Kielgast J, Iversen LL, Moller PR, Rasmussen M, Willerslev E. Detection of a Diverse Marine Fish Fauna Using Environmental DNA from Seawater Samples. PloS ONE. 2012; 7(8): e41732. pmid:22952584
  18. 18. 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(7): e6372. pmid:19633714
  19. 19. Bucklin A, Frost BW, Bradford-Grieve J, Allen LD, Copley NJ. Molecular systematic and phylogenetic assessment of 34 calanoid copepod species of the Calanidae and Clausocalanidae. Marine Biology. 2003; 142: 333–343
  20. 20. Becker RA, and Wilks AR. Original S code. R version by Brownrigg R. Enhancements by Minka TP. tpminka@media.mit.edu maps: Draw Geographical Maps. R package version 2.3–9. http://CRAN.R-project.org/package=maps 2014
  21. 21. Bucklin A. Methods for population genetic analysis of zooplankton. 2000; In: The ICES Zooplankton Methodology Manual. International Council for the Exploration of the Sea. London: Academic Press; p. 533–570. Chapter 11
  22. 22. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al QIIME allows analysis of highthroughput community sequencing data. Nature methods. 2010; 7(5): 335–6. pmid:20383131
  23. 23. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006; 22(13): 1658–9. pmid:16731699
  24. 24. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010; 26(19): 2460–1 pmid:20709691
  25. 25. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Research. 2007; 35(21): 7188–96. pmid:17947321
  26. 26. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011; 27(16): 2194–200. pmid:21700674
  27. 27. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology. 2007; 73(16): 5261–7. pmid:17586664
  28. 28. Oksanen J, Blanchet F.G, Kindt R, Legendre P, Minchin PR, O’Hara RB et al vegan: Community Ecology Package. R package version 2.0–10. http://CRAN.R-project.org/package=vegan. 2013
  29. 29. McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013; 8(4). e61217. pmid:23630581
  30. 30. Lozupone C, Knight R, UniFrac: a new phylogenetic method for comparing microbial communities. Applied and Environmental Microbiology. 2005; 71; 8228–8235 pmid:16332807
  31. 31. Kelley D. oce: Analysis of Oceanographic data. R package version 0.9–14 http://CRAN.R-project.org/package=oce. 2014
  32. 32. Machida RJ, Knowlton N. PCR Primers for Metazoan Nuclear 18S and 28S Ribosomal DNA Sequences. PLoS ONE. 2012; 7(9): e46180. pmid:23049971
  33. 33. Prokopowich CD, Gregory TR, Crease TJ. The correlation between rDNA copy number and genome size in eukaryotes. Genome. 2003; 46(1): 48–50. pmid:12669795
  34. 34. Echelman T, Fishelson L. Surface zooplankton dynamics and community structure in the Gulf of Aqaba (Eliat), Red Sea. Marine Biology. 1990; 107(1): 179–90.
  35. 35. Böttger R. The Vertical Distribution of Micro- and Small Mesozooplankton in the Central Red Sea. Biological Oceanography. 1987; 4(4): 383–402.
  36. 36. Cornils A, Schulz J, Schmitt P, Lanuru M, Richter C, Schnack-Schiel SB. Mesozooplankton distribution in the Spermonde Archipelago (Indonesia, Sulawesi) with special reference to the Calanoida (Copepoda). Deep-Sea Research Part II-Topical Studies in Oceanography. 2010; 57(24–26): 2076–88.
  37. 37. Khalil MT, ElRahman NSA. Abundance and diversity of surface zooplankton in the Gulf of Aqaba, Red Sea, Egypt. Journal of Plankton Research. 1997; 19(7): 927–36.
  38. 38. Pancucci-Papadopoulou MA, Siokoufrangou I, Theocharis A, Georgopoulos D. Zooplankton vertical distribution in relation to the hydrology in the NW Levantine and the SE Aegean Seas (Spring 1986). Oceanologica Acta. 1992; 15(4): 365–81.
  39. 39. Moore E, Sander F. Quantitative and Qualitative Aspects of the Zooplankton and Breeding Patterns of Copepods at two Caribbean Coral Reef Stations. Estuarine and Coastal Marine Science. 1976; 4: 589–607.
  40. 40. Weikert H. The vertical distribution of zooplankton in relation to habitat zones in the area of the Atlantis II Deep, Central Red Sea. Marine Ecology Progress Series.1982. 1982; 8; 129–143
  41. 41. Cornils A, Niehoff B, Richter C, Al-Najjar T, Schnack-Schiel SB. Seasonal abundance and reproduction of clausocalanid copepods in the northern Gulf of Aqaba (Red Sea). Journal of Plankton Research. 2007; 29(1): 57–70
  42. 42. Farstey V, Lazar B, Genin A. Expansion and homogeneity of the vertical distribution of zooplankton in a very deep mixed layer. Marine Ecology Progress Series. 2002; 238: 91–100.
  43. 43. Haury LR. Vertical distribution of Pleuromamma (Copepoda, Metridinidae) across the eastern North Pacific Ocean. Hydrobiologia. 1988; 167: 335–42.
  44. 44. Longhurst A, Williams R. Materials for plankton modeling vertical distribution of Atlantic zooplankton in summer. Journal of Plankton Research. 1979; 1(1):1–28
  45. 45. Nishikawa J, Matsuura H, Castillo LV, Campos WL, Nishida S. Biomass, vertical distribution and community structure of mesozooplankton in the Sulu Sea and its adjacent watrs. Deep Sea Research II. 2007; 54: 114–130
  46. 46. Böttger-Schnack R, Schnack D, Hagen W. Microcopepod community structure in the Gulf of Aqaba and northern Red Sea, with special reference to the Oncaeidae. Journal of Plankton Research. 2008; 30(5): 529–550
  47. 47. Kosobokova K, Hirche H-J. Zooplankton distribution across the Lomonosov Ridge, Arctic Ocean: species inventory, biomass and vertical structure. Deep Sea Research I. 2000; 47; 2029–206