The Horizontal Distribution of Siliceous Planktonic Radiolarian Community in the Eastern Indian Ocean

Abstract

The plankton radiolarian community was investigated in the spring season during the two-month cruise ‘Shiyan1’ (10 April–13 May 2014) in the Eastern Indian Ocean. This is the first comprehensive plankton tow study to be carried out from 44 sampling stations across the entire area (80.00°–96.10° E, 10.08° N–6.00° S) of the Eastern Indian Ocean. The plankton tow samples were collected from a vertical haul from a depth 200 m to the surface. During the cruise, conductivity–temperature–depth (CTD) measurements were taken of temperature, salinity and chlorophyll a from the surface to 200 m depth. Shannon–Wiener’s diversity index (H’) and the dominance index (Y) were used to analyze community structure. There was a total of 168 plankton species, composed of Acantharia, Phaeodaria, Polycystina, Collodaria and Taxopodida (monospecific—Sticholonche zanclea, Hertwig is the only recognized species). Hence, it included both celestine-based and siliceous organisms, which are also described here for the first time from this region. Total radiolarians ranged from 5 to 5500 ind/m⁻³, dominated by co-occurrences of Sphaerozoum punctatum and Stichonche zanclea species at the south-equator zone (SEQ)-transect 80° E and equator zone (EQ)-transect Lati-0. The possible environmental variables were tested through RDA analysis; although no result was obtained for the full species dataset, the samples from the equatorial transect related strongly to mixed-layer chlorophyll a concentration and those of a north–south transect to surface silicate concentrations or mixed-layer nitrate were significantly correlated (p < 0.01) to the radiolarian community. Our results indicate that the silicate and chlorophyll-a concentrations are the two major factors affecting the radiolarian distribution along two of the investigated transects (southern equator and equator) in the study area.

Highlights:

- First comprehensive plankton-net sampling of Eastern Indian Ocean.

- Apparent dependence of some observed taxa of particular oceanic environmental variables.

- Explanatory environmental variables for plankton abundances in two cruise transects.

1. Introduction

Radiolaria are the siliceous-shelled, single-cell, planktonic organisms frequently studied for palaeoceanographic and environmental reconstruction. They are globally distributed and widely collected from plankton tow and sediment samples around the world’s oceans [1]. They are one of the most diverse group found in all water depths of the ocean [2] and a good indicator of water masses [3,4]. Radiolarian protists have opaline silica (SiO₂-nH₂O), strontium sulfate or calcium carbonate skeletons and are characterized by central capsule shells and axopodia spines [5]. They also contribute significantly to the deposition of silicate from the ocean’s upper layer to the bottom sediments [6,7]. They are the second most important producers of biogenic opal and siliceous oozes in the oceans, exceeding the effect of diatoms in equatorial regions [8]. Usually, they contribute 14 to 31% of the total bSiO2 export in the Ocean [9,10].

The biogeographical distribution of radiolaria has been reported from the Atlantic Ocean [11,12,13,14], the Pacific Ocean [15,16,17,18,19,20,21,22], the Arctic Sea [23] and the Indian Ocean [24,25,26,27]. The Eastern Indian Ocean (EIO) is part of one of the three largest oceans globally and is geographically located at the equator, covering 1000 km²–5000 km² area [28,29]. It is divided into two semi-enclosed basins from the Bay of Bengal (BOB) and the Arabian Sea [28]. The Eastern Indian Ocean has complex oceanography due to semiannual prevailing currents [28,30]. It is characterized by four mainstream ocean currents (Figure 1), which are known as BBR (Bay of Bengal runoffs), EICC (east India coastal current), WICC (west India coastal currents), EK (Ekman transport), SEC (south equatorial counter-current), SJS (south Java current), SMC (northeast and southwest monsoon currents) and WJ (Wyrtki jets) currents [28,29,31,32]. The EICC flows equatorward along the western boundary of the BOB before crossing to the east of Sri Lanka. Along the Sri-Lankan coast, high nutrient concentrations are entrained in the EICC [31]. The WJs are eastward surface currents that only occur in the equatorial EIO during the spring (April–May) and result in a deep upper mixed layer [33,34]. On the eastern side of the basin, the WJs cause a decline in the thermocline and nitracline, leading to decreases in phytoplankton primary production in May [33,34,35]. The conductivity–temperature–depth (CTD) profiles from 0 to 200 m and nutrients from the World Ocean Atlas (WOA13) during the boreal spring season (April–May) in the Eastern the Indian Ocean, described from the three transect areas, south-equator zone (SEQ), north-equator zone (NEQ) and equator zone (EQ) (Figure 2 and Figure 3a–l). The temperatures ranged from 15 to 30 °C and salinity varied from 35.0 to 36.8 psu, maximum SST and SSS are reached in the south-equator during this season, compared to the northern BOB and Sumatran coastal area [29], which becomes stratified. The nitrate and silicate concentrations are 6.09 µM to 40.38 µM, which has no obvious different distribution trend (Figure 3). The concentrations of nitrite and silicate are higher in the northern part of the NEQ transect area (Figure 3a,c), while the concentration near the equator is low (Figure 3e,g). Phosphate values are 0.38 to 2.5 µM and the maximum is reached in the equator and Sumatran station (Figure 3b,f). The dissolved oxygen (DO) values are 0.1–4.50 mg/L and the maximum also occurs in the equatorial region (Figure 3h).

Only a few studies have been conducted in the Eastern Indian Ocean and most of them are based on surface sediment samples and the investigation of polycystine radiolarians [24,25,27]. This is the first comprehensive plankton tow sampling study in the Eastern Indian Ocean. The research objectives were to estimate the diversity and horizontal distribution of the radiolarian community from 200 m to the surface of the Eastern Indian Ocean and find out the possible explanatory environmental factors between three transects; NEQ in the north; SEQ in the south-equator and EQ in the central region of the Eastern Indian Ocean. Radiolarian are key players in biogeochemical processes such as transformation of silica and carbon. Their composition and distribution characteristics in the EIO are of great significance for future ecological and molecular diversity research.

2. Materials and Methods 2.1. Sampling Sites and Collections

Plankton-net tow sampling was conducted in the Eastern Indian Ocean (80.00°–96.10° E, 10.08° N–6.00° S). The study area covered the ca. 1000–1500 km distance from the Bay of Bengal (BOB) and parallel to the Sumatran coast of Indonesia to equatorial section and southern equatorial zone (Figure 2). The first transect (hereafter “Long-90”) consisting of the 17 stations along the NNW line from 5.5° S, 96.1° E (station I511) to 10.08° N 88.1° E (station IQ2), the second (hereafter, “Long-80”) of 15 stations along two adjacent north–south lines from 2.42° N, 80.0° E (station I302) to 6.0° S 80.0° E (station I318) and the third, which partially overlapped the other two 11 stations (hereafter “Lat-0”), closely followed the equator (Table 1).

Using a modified Indian Ocean Standard Net (20 µm mesh size, 0.57-m diameter, 470 cm, with a mechanical flow meter), samples were collected from a vertical haul from a depth 200 m to the surface. About 1000 mL of sample was collected in each PM bottle and fixed with 2% formalin preservative. Preserved samples in the PM bottles were returned to Dr. Sun’s laboratory (Tianjin University of Marine Science and Technology, China) for analyses. 2.2. CTD Data Collection

Temperature and salinity data were collected using a Seabird conductivity–temperature–depth (CTD) system (SeaBird SBE 911 plus V2). Using the Strickland and Parson [36] method, a standard protocol to measure the chlorophyll a, was employed using 500 mL of seawater filtered through a 25-mm diameter Whatman filter paper and kept in 5 mL of acetone (93%) at 20 °C for 24 h, then analyzed using a fluorometer (Trilogy, Turner Designs, Sunnyvale, CA) for estimation of chlorophyll a.

2.3. Satellite Data

Satellite remote sensing on sea-surface chlorophyll a data (MODIS-aqua level 5 km resolution) was obtained during the spring season (April–May 2014) to evaluated evidence for Eastern Indian Ocean from NASA’s World view center (NASA’s Worldview website); (https://worldview.earthdata.nasa.gov?v=57.9508133990794216.300628602547988,126.0113303398581,17.116840099466998;t=2014-04-23-T15%3A13%3A42Z;l=MODIS_Aqua_Chlorophyll_A,MODIS_Terra_Chlorophyll_A).

2.4. Nutrients Concentration

The nutrients, silicate, nitrate, phosphate and dissolved oxygen values were derived from the interpolation of the objectively-analyzed values during boreal Spring data (April to June) collected by World Ocean Atlas 2013 (https://www.nodc.noaa.gov/OC5/woa13/woa13data.html) at the 25 m, 50 m, 75 m, 100 m, 150 m and 200 m below sea level (Figure 3a–l). The concentration of nitrate ranged from 6.9 to 16.7 µM), silicate 0.1–11.7 µM, phosphate 0.1–1.69 µM, dissolved oxygen 0.1–4.50 mg/L (Figure 3). The maximum concentrations of nitrate and silicate was observed at the BOB stations (HFA01; I601) at 90° E (Figure 3a,c) which decreases to the equator at stations (I312, I409, I101) (Figure 3e,g) and the maximum phosphate concentration was also found from the equator—EQ station (I410) (Figure 3f), Sumatran station (I511) and BOB stations (I607, HFA-01) (Figure 3b). The dissolved oxygen (DO) was maximum in two transects at 80° E and 90° E (Figure 3d,l) but the oxygen saturation was high 34 mg/L recorded in the equator-EQ (Figure 3h).

2.5. Community Analysis and Microscopic Observation

To enumerate the community structure and for primary identification, fixed samples were examined at 20× and 40× magnification under an inverted light microscope (Motic, AE 2000, Xiamen, China) which was equipped with a digital Camera (Moticam 2506, 5.0 m, Pixel, Xiamen, China). About 2 mL of each sample were examined to count cells in a settling chamber using a modified Utermöhl method [37]. Species identification was confirmed by the use of scanning electron microscopy (SEM). Samples were pre-screened through a net again to remove large zooplanktons and subsamples were filtered on 0.6 µm, pore-size Millipore filters and air-dried for 24 h, then mounted on a stub, sputter-coated and examined under the SEM microscope (JEOL JSM-IT300LV, Hitachi, Japan) at Biological building (13), Tianjin University of science and Technology, TEDA, China. Species were identified using published taxonomic references [13,38,39,40].

2.6. Data Analysis and Statistical Methods

Shannon–Wiener diversity index (^H′ ), Pielou evenness index (J) and Dominance index (Y) has been used by the modified equation [37] to analyze the composition and richness.

^H′=−∑i=1S_{Pi log2 Pi}

Y=_niN_fi

J=^H′/_log2S

where S is the phytoplankton species of each sample [41]. N is the total number of cells of all species counted; ni is the cell number of species i; Pi = ni/N represents the relative abundance of a given species; fi is the frequency of occurrence of species i in each sample. Species with Y > 0.02 was defined as dominant species [42].

2.7. Correspondence Statistical Analysis

The Hellinger-transformed Legendre and Gallagher [43] of species data were investigated by redundancy analysis (“RDA”) using the R software [44] RDA was chosen to allow the use of ordiR2step (Blanchet et al. [45] and an equivalent of Canoco (ter Braak and Šmilauer, [46] for the coefficient freedom value plus Pearson correlation values)).

3. Results 3.1. Environmental Dataset of Spring 2014 in Eastern Indian Ocean

An environmental dataset was constructed comprising the values for temperature, salinity and chlorophyll a measured during the cruise at the depths (0 m, 50 m, 75 m, 100 m and 200 m) were derived from CTD data (Figure 4a–i).

3.2. Vertical and Horizontal Profile of Temperature and Salinity

The vertical distribution of temperature and salinity in the three transects are summarized in Figure 4). The temperature and salinity in the depth 0–50 m ranged from 12.82 and 32.55 psu and in the depth 100–200 m ranged from 31.82 °C and 32.5 to 35.5 psu, respectively (Figure 4). Horizontally, maximum temperature value (>31.82 °C) was observed in the BOB samples (HFA01, HFA05, I601, I603, IQ2) and (<30 °C) was observed in the equator samples (I401-I413) (Figure 5a). The maximum salinity was observed also in the equator stations (I308, I401, I403, I406) at the transect 80° E and minimum value of salinity (>32 psu) was observed in BOB sampling areas (IQ2, HFA01-HFA05) (Figure 5b). The low salinity values in the BOB could be influenced by riverine runoff inputs and aggravate stratification. However, the salinity remains steady in the area nearby the Sumatran coast of Indonesia (Figure 4b and Figure 5b), which probably contributed to an increase in salinity gradients due to riverine influences from these coastal locales.

3.3. Vertical and Horizontal Profile of Chlorophyll a

The vertical distribution of Chlorophyll a in the Eastern Indian Ocean are shown in (Figure 4g–i). Chl-a ranged from 0.005 to 0.89 µg/L, which was higher at the 50 m and 75 m depth and minimum Chl- a was observed at 200 m (Figure 4c,f,i). Horizontally, the high chlorophyll a was observed in the equator region at the stations (I403, I308, I410-I412) and the Sumatran coast at the stations (I501, I503) (Figure 5c) respectively.

3.4. Community Structure and Diversity Pattern

In total, 168 taxa were identified and classified into 28 families and four orders (Table S1). Table S1 includes the number of individual (ind/m⁻³) of each species/sample. The identified taxa were composed of Polycystinea groups, e.g., Collodaria (15 spp.; Figure 6 and Figure 7), Spumellaria (44 spp.; Figure 8), Nassellaria (75 spp.; Figure 9), followed by Acantharia (8 spp.) and Taxopodida (1 spp.; Figure 10) and Phaeodaria; (7 spp.; Figure 11).

Out of 168 taxa, the 5 to 40 species were horizontally distributed in the study area (Figure 12a). The species richness was highest at stations (I308, I815, I817) and (I509, I607, I609) between the two transects area at the 80° E and 90° E (Figure 12a). A high diversity index (2.5 H’) was observed in stations (I304, I310, I412, I413, I809, I813, I817) at the transect 80° E and (HF-05) at the BOB zone. The low diversity (1–1.5 H’) was observed in the stations (I410, I807, I611, HF-01) (Figure 12b).

3.5. Plankton-Net Tows Abundance of Radiolarian

Radiolarian abundance varied from 5 to 5838 ind/m⁻³ across the sampling area. The highest abundance was observed at the station (I815), followed by (2200 ind/m⁻³) at stations (I310, I308) and low abundance (247 ind/m⁻³) at station (I410), (174 ind/m⁻³) at station (I312), (168 ind/ m⁻³) at station (HFA03), (157 ind/m⁻³) at station (I506) and (109–102 ind/m⁻³) at stations (I609, I306, I815), respectively (Figure 12c). Collodaria and Taxopodida were the most abundant group (90% and 70% of the total radiolarian abundance), respectively (Figure 12d,e), followed by Acantharia (35%) (Figure 12f), Spumellaria and Nassellaria (15% and 6%) (Figure 12g,h) and Phaeodaria (8%) (Figure 12i), respectively. The top three species, S. punctatum, S. fuscum, S. zanclea contributed 56.6% and other fifteen species contributed 10–25% to total radiolarian (Table 2).

3.6. Horizontal Distribution of Radiolarian in EIO The abundance, diversity and taxonomic patterns within the three transect zone show the following characteristics.

• NEQ transect—90° E: Radiolarian abundance ranged from 6 to 1080 ind/m⁻³ which recorded from 17 stations at transect 90° E (Figure 12c). These samples contained a total of 89 species. Each species has heterogeneous distribution, 2–23 species were found from each station and maximum assemblage occurred at the stations (I609) (Figure 12b). The maximum abundance was recorded at the BOB stations (AH-05), followed by 661 ind/m⁻³ at stations (I109, I609) which were dominated by the occurrences of Stylodictya multispina, C. inerme and Siphonosphaera polysiphonia (Figure 13; Table 2) and low abundance <250 ind/m⁻³ was recorded by the occurrences of S. punctatum and Callimitra sociliatorchia, <190 ind/m⁻³ by the occurrences of S. zanclea and Didymocyrtis tetrathalamus (Figure 13), respectively.

• Equatorial transect—0°: The abundance of radiolarian in the equatorial samples area ranged from 14 to 996 ind/m⁻³ and collected from 11 stations (Figure 12c). These samples contain total of (39 species) and approximately, 2–15 species were found from each station (Figure 12b). The diversity increases maximum of 15 species at stations (I403, I412, I413) but subsequently drops to two species at other stations. The assemblage in the equatorial region is dominated by S. punctatum, S. cf ovodimare, L. pentagona pentagone and S. zanclea species (Figure 13).

• SEQ transect—80° E. The abundance of radiolarian in the transect 80° E (SEQ) ranged from 1000 to 5257 ind/m⁻³ and recorded from 15 stations (Figure 12c). The high species richness numbers (41 species) occurred at station (I308) and is immediately followed by low species numbers (27 species) at station (I509) and (23–21 species) occurred at stations (I310, I312) (Figure 12b), respectively. The maximum abundance >5257 ind/m⁻³ was observed at the station (I817) and assemblages are dominated by S. punctatum, followed by >1900–1780 ind/m⁻³ observed at stations (I310, I807) dominated by S. fuscum and S. zanclea (Figure 13), respectively.

3.7. Selected Dominant Taxa in the Eastern Indian Ocean

We selected 20 dominant species which identified by dominant index (Y) and frequency distribution (fi) (Table 2). Overall, relative contribution of S. punctatum (22%), S. fuscum (15%), D. tetrathalamus (3.5%), S. zanclea (7%), S. polysiphonia (2.3%), L. pentagona pentagone (2.3%), Eucytidium hexastichum (32%), C. sociliatorchia (2%), Pterocanium pretextum (1.2%), Hexacontium armatum-hostile, E. gegenbauri, Anthocyrtidium cineraria and Arachnocorys circumtexta (0.8%), Lophophaena dodecantha (0.7%), Amphirrhopalum ypsilon, Acanthsphaera actinote, Euchitonia elegans and Spongodiscus resurgens (0.6%) (Table 2), respectively.

4. Canonical Ordination Analysis 4.1. Identification of Potentially Explanatory Environmental Variables

The RDA-generated inertia and the first ten eigenvalues appear as Table 3.

Two tests were employed to determine how many, if any, variables the eigenvalues from the species dataset might explain. The first test used the Guttman–Kaiser Criterion (“GKC”) (Yeomans and Golder [47]): Fourteen eigenvalues exceeded the mean eigenvalue, although two only marginally, suggesting many explanatory variables. The second, more stringent, test employed the “Broken Stick” Model (“BSM”) Jackson [48] and this indicated that there might be about two explanatory variable. The test results are plotted as Figure 14a.

An environmental dataset was constructed comprising the values for temperature, salinity and chlorophyll measured during the cruise at the ocean’s surface and depths of 25, 50, 75, 100, 150 and 200 m below sea level (bsl) plus the World Ocean Atlas 2013 objectively-analysed values for nitrate, phosphate, silicate and dissolved oxygen for the season April–June at the same depths [49,50]. RDA was applied to the species data and each variable in the environment dataset in turn and Analysis of Variance (ANOVA) run on the output: Those environmental variables where the probability of the F statistic was greater than 0.005 and the ratio of the first RDA axis to the first unconstrained axis was greater than one was to be published. No significant results were found. Because the species data were drawn from a relatively large ocean region with a range of environmental conditions, three transects were selected to examine whether plankton abundances were affected by different environmental variables in different parts of the study area (Figure 2 and Table 1).

Application of the GKC to the transect species data suggested that Long-90 and Long-80 might each have about four explanatory variables and Lat-0 three. The BSM indicated that no explanatory variable would be revealed. RDA and ANOVA were applied to each of the transects with the same publication criteria as for the full species dataset. Transect Long-80 gave five and Lat-0 gave two results significant at the 0.1% level; Long-90 yields no significant results—See Table 4.

An R program was written to use ordiR2step to find the most likely additional explanatory environmental variables given the five already found for transect Long-80 and the two for Lat-0. The program’s results suggested some secondary explanatory variables but, when strongly correlated pairs (see Table 5) were eliminated, none remained. The most likely explanatory environmental variables (based on the Pr(>F) statistic) are, therefore, for the Long-80 transect is silicate and for Lat-0 chlorophyll a, both at the sea surface.

4.2. Selected Dominant Taxa and Environmental Factor

The dominant species were explained by RDA analysis with six environmental variables (Figure 15; Table 6). The species showed the positively correlation to temperature and salinity, e.g., S. zanclea, S. polysiphonia, S. fuscum, H. armatum-hostile, D. tetrathalamus, E. elegans, S. resurgens and C. sociliatorchia (Table 6). The species, S. zanclea, D. tetrathalamus, A. hexagonalis, E. elegans, L. pentagona pentagona and L. dodecantha were positively correlated to chlorophyll a (Figure 15; Table 6), but only L. pentagona pentagona was significantly correlated (p < 0.05). The species, S. zanclea, S. polysiphonia, S. punctatum, A. actinote, H. armatum-hostile, A. ypsilon and E. elegans was positively correlated to phosphate (Figure 15; Table 6) and only S. polysiphonia was significantly correlated (p < 0.05) (Table 6). The species, S. zanclea, S. polysiphonia, S. punctatum, E. elegans, A. circumtexta, C. sociliatorchia and L. pentagona pentagone were also positively correlated to silicate (Figure 15; Table 6) and S. polysiphonia showed the significance level (p < 0.01) (Table 6).

5. Discussion

There are only a few sediment sampling studies of radiolaria have been conducted in the Eastern Indian Ocean so far [24,25,26]. We describe a comprehensive study of the planktonic radiolarian community and its horizontal distribution on the three transects in the Eastern Indian Ocean. We made no distinction between living and dead radiolarian specimens, the samples were immediately fixed with formalin after collection and no evidence was gathered to determine how many were living at the time of collection, therefore our samples represent both the living and dead individuals that were in the water column from sea surface to 200 m depth at the time of plankton net sampling.

5.1. Radiolarian Communities in the Eastern Indian Ocean

Preliminary studies reported only the 35 taxa from the Western Indian Ocean [24], 74–261 taxa from the Eastern Indian Ocean [25,27], 40 taxa from the BOB [26]. In this study about 168 taxa of the major radiolaria group, such as Polycystinea, Acantharia, Collodaria, Taxopodida and Phaeodaria were collected from plankton tow samples from the surface to euphotic zone of EIO (Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11).

Most of these groups were tropical or sub-tropical Pacific and Atlantic Ocean dwelling species [3,16]. On the base of previous geographical records, the distribution of radiolarian species in the EIO does not reflect their Pacific and Atlantic counterparts, was described here for the first time in the EIO (Table 7). However, the species of Polycystinea group were also retrieved, which were occasionally important (Supplemental Table S1; Figure 8 and Figure 9). The Collodaria and Taxopodida group were mainly represented by the high frequency distribution 70–90% which were dominant in our samples (Figure 12d,e), have a widespread distribution and contribute high fluxes within phytoplankton communities in open waters [51]. Nonetheless, this is the first report on the co-occurrence of Collodaria and Taxopodida in the Eastern Indian Ocean. However, variations of Collodarian contribution to the radiolarians were lower in the Antarctic and Indian Ocean [51]. The species of Spumellaria (15%) and Nassellaria (6%) were found in relatively low abundance in all samples (Figure 12g,h). The Spumellaria and Nassellaria species were highly abundant in the BOB [26] and in the Indian Ocean [1]. Spumellaria and Nassellaria have a global distribution in the surface and sediments waters of the Atlantic and Pacific Ocean [1]. More recently, a high percentage of Spumellaria were correlated with a low percentage of Collodaria and Nassellaria in the South China Sea [52]. We additionally found the Acantharia and Phaeodarian in low numbers in our samples, contributing only about 8–13% to the total radiolarian population (Figure 12f,i) and notably occurred in the twilight zones of the Ocean [51]. Acantharia are usually found in the tropical region of the Atlantic Ocean and East China Sea (Table 7) and Phaeodarian are predominantly found in the Adriatic Sea and Arctic Sea [53]. However, Gupta et al. [26] reported the few Phaeodarian species from sediments of BOB. Despite their low abundances, they are a significant component because they are celestine and opaline silica protists. We can conclude that more work is needed to examine the sampling depth per station and to distinguish between living and dead radiolaria in the Eastern Indian Ocean.

5.2. Characteristics of Radiolarian Abundance, Horizontal Distribution and Environmental Factor

Diversity indices and abundance was higher around the equator and SEQ transect regions (Figure 12). One of the reasons for high diversity was the occurrence of two species, S. punctatum and S. zanclea, in these two transect area. The diversity of Collodaria species such as S. punctatum and S. polysiphonia can be influenced by the trophic status of oceanic provinces and increased towards more oligotrophic regions [51]. As a result, high abundance in the equatorial region at stations (I815, I310, I312) can be predicted by the Wyrtki jets in spring season when WJs carries the high salinity and high nutrient waters from west to eastward [34]. However, low diversity was observed at the northern BOB and near Sumatran coastal area due to less saline waters due to runoff from the Ganges river and southwest upwelling region from these coastal locales. Radiolaria are rare in coastal areas, bays, estuarine and fjord areas due to less saline water [21,25]. This result is contrary to recent research in the Pacific Ocean, the radiolarian abundance and distribution in the continental shelf area of the Pacific Ocean was reported under low salinity and high temperature [66].

Radiolarian abundance in the Eastern Indian Ocean during this study was ranged from 5 to 5850 ind/m⁻³ (Figure 12). The high abundance at the one station (I815) from the south equator transect, where the colonies of S. punctatum and S. fuscum were very abundant (Figure 6d,f). Over all, abundance was below the 1500 ind/m⁻³. Our results generally agree with previous information that radiolarian densities are higher in the shallow water [1,16]. Plankton net-studies in the Pacific and Atlantic Ocean also show that the highest abundance of radiolarian occurs in the upper ocean layer [21]. The abundance of radiolarians was similar to the other studies from the Atlantic Ocean [13] and Indian Ocean [1]. It ranged from 1 to 100 m⁻³ and 500 m⁻³ and polycytinea radiolarian numbers are ~1200 times greater at 0–150 m than at 500–1000 m [13]. At the depths (100 m), the maximum densities such as (1400 specimen/m⁻³) and low densities (250–500 ind/m⁻³) can be observed and reported from the Pacific Ocean [4,16]. Boltovoskoy and Correa, [16] also recorded the maximum abundance such as (~200,000 specimens/m⁻²) at 200 m and (~300,000 specimens/m⁻²) at 0–400 m from 30° N to 30° S. Our results recorded the Taxopodida species, S. zanclea at maximum abundance (1780 ind/m⁻³) and Collodarian species, S. punctatum (5500 ind/m⁻³) and S. polysiphonia (661 ind/m⁻³) are the most dominant species from integrated depth 0 to 200 m (Figure 13a–c). S. punctatum and S. polysiphonia are the cosmopolitan species form large colonies in surface waters to 100 m depth, estimated at 30 to 20,000 colonies m⁻³ [67]. However, S. polysiphonia species are only encountered along the transect longitude-90° E (Figure 13c) and are restricted to shallow photic layers of the ocean [5].

The radiolarian assemblage in our data consisted of a high abundance of photosynthetic symbionts species and subtropical species such as S. punctatum, S. zanclea, E. hexastichum, D. tetrathalamus, A. circumtexta, L. pentagona pentagone, A. actinote, A. cineraria, H. armatum-hostile, E. gegenbauri. Most of the polycystinea radiolarian species are the common dominant species in the EIO [24,25,26], the Pacific Ocean [1,16,66]. Among of them, the S. punctatum and S. zanclea has new geographical distribution records in the EIO. High percentages of S. punctatum, S. zanclea, D. tetrathalamus, are common in the modern tropical-subtropical the Pacific Ocean and Atlantic Ocean [2,18,21,68]. The horizontal distribution of S. punctatum, S. fuscum and S. zanclea, with high abundance and L. pentagona pentagone, D. tetrathalamus, E. hexastichum, C. sociliatorchia low abundance were present in the transect 80° E region which can be indicator of the SCM (SEC: South Equatorial Countercurrent) (Figure 13), while S. punctatum L. pentagona pentagone and S. zanclea dominated in the equator region (Figure 13a–c,e–i) and are the indicator species of WJs (Wyrtki jets currents), While S. polysiphonia, S. punctatum, C. sociliatorchia, H. armatum-hostile, A. hexagonalis, S. zanclea and D. tetrathalamus were commonly found in the transect 90° E (Figure 13a,b), are the indicator species of BBR (Bay of Bengal Runoffs) and SJC: South Java Current. Our RDA results show that the abundances of S. punctatum and A. cineraria are positively correlated to temperature (Pearson rank correlation coefficient r > 0.31, p < 0.05), S. polysiphonia negatively correlated with salinity (r > −0.431, p < 0.001) and with phosphate (r > −0.431, p < 0.001) and L. dodocantha with Chl a concentration (r > 0.004, p < 0.005). D. tetrathalamus and H. armatum-hostile, abundance was negatively correlated with silicate concentration (r = −0.318, p < 0.005) and (r = −0.045, p < 0.005) (Table 6).

5.3. Factors Influencing on the Radiolarian Distribution

The vertical water structure of 0–50 m has the characteristic of high temperature and low salinity and high nutrients. At the depth below 100–200 m, it is characterized by high salinity and low temperature as shown in Figure 4. BOB is located in the north which has the characteristics of intrusive river runoff with low salinity and high temperature in the EIO [31], Sumatra is open coastal zone with conjunction of the BOB, characterized by low salinity and oxygenated water and the equator area is characterized by warm water and mixed layer of Chlorophyll a [30].

In our study area, the temperature and salinity of the equator are high and the boreal-spring nutrient (WOA13) data shows that the equatorial-nutrient concentration is low (Figure 3e–g). The RDA analysis gave most of the negative but significant correlation values between the environmental variables (Table 5). RDA explains the positive correlation at level (p < 0.05) between the Chlorophyll a, nitrate and temperature and between the salinity, silicate and phosphate (Table 5). Dissolved Oxygen was negatively correlated, indicating the minimum OMZ condition. This suggests that more than one variable, e.g., Chlorophyll a, temperature and silicate can shape the radiolarian community in the EIO [69,70].

5.4. Factor Influence on Radiolarian in Three Transect Area, EIO

Salinity and temperature are two significant factors that are closely related to the radiolarian community at the depth of 200 m [25,26,27,52,70,71]. High temperature can be associated with high radiolarian diversity [52]. Rogers and De Deccker, [27] clearly demonstrate that, temperature and salinity between the 100 m and 300 m and nutrients at the 100 m in the Eastern and Southern Indian Ocean, have strongly impacted on the distribution of the radiolarian community. In the present study, the relationships between the all species dataset, within each transect, with seven environmental variables often did not provide any significant correlations (Figure 14a), but significant relationships were also found with the dominant taxa at the transect 80° E and equator (Figure 14). The ratio of first RDA axis to the first unconstrained axis was greater than 1 and probability of F statistic was significant 0.005, but they did not provide a significant correlation at the transect long-90° E. The difference between the initial and subsequent unconstrained eigenvalues was not great, an attempt was made to strengthen the effect by excluding those species with an RDA score of less than two standard deviations from the mean score in any of the first six components—That is using only the dominant species. The ratio of the first to the second unconstrained eigenvalues did increase (Table 3—labelled “Dominant”) but application of GKC, BSM and RDA failed to produce any better results than using all the species data. It was noted that the species data matrix was very sparse (approximately 90% of the cells contained zeros). A dense species dataset was generated by deleting the sparsest species records and the battery of tests described above applied. No significant results were obtained. The variance in the full species data cannot be explained by a small number of variables probably because the study area is large enough to cover a number of ocean environments. Analysis of two sections of the total species data revealed possible explanatory variables: (a) The variance in species from a north–south transect at approximately 80° E (“Long-80”) could be due to the surface silicate concentrations but, alternatively, might result from nitrate concentrations or dissolved oxygen in the mixed layer. Silicate, nitrate and oxygen are necessary for radiolarian growth so the connection with the variance is likely to be causal. (b) The variance in an equatorial transect (“Lat-0”) is strongly related to the chlorophyll a level in the mixed layer but it is not certain that chlorophyll a is necessary for radiolarian growth. Thus, the connection may be indirect such as radiolarians may consume chlorophyll-producing organisms or may simply require a similar environment. However, neither temperature nor salinity showed significant correlations with the radiolarian distribution (Table 4). As salinity and temperature enable radiolarians to thrive in different water masses [16]; its effect on the vertical distribution is greater than that on the horizontal distribution [23].

5.5. Factors Influencing the Radiolarian Communities

Temperature and salinity are not sufficient explanatory variables: Chlorophyll a and nutrients and dissolved organic matter (DOM) can stimulate heterotrophic communities such as the radiolarian [4,27,71,72], tintinnid [29] and diazotroph communities [73,74]. These environmental variables play an important role in circulatory changes of EIO in the spring season. The nutrients (silicate, phosphate and nitrate) may be depleted throughout the EIO, except in coastal areas, when coastal upwelling cause nutrients to elevate [31]. In this study, the Chlorophyll a concentration in the surface layer ranged from 0.005 to 0.89 µg/L, which is similar to the remote sensing Chlorophyll-a data and high chlorophyll a of the same period can be seen in this region (Figure 4; Supplemantry Figure S1). The Chlorophyll a value is significantly higher in spring season in same region [75]. This strengthens the possibility of using Chl-a remote sensing data to estimate the geographic distributions of radiolarian in the Ocean. These concentration values often lead to low nutrients in the DCM layer in spring [29,31]. Our findings also confirmed earlier reports that low nutrients at the equator correspond to the DCM (deep Chl a concentration maximum) layer from 50 to 100 m in a same study area [29]. This DCM chlorophyll a maximum has caused the dominance of nano and picoplankton biomass [29,31] including some heterotrophs bacteria [73] in the region. The silicate, nitrate or oxygen concentrations including chlorophyll a in surface water are mostly highly correlated at the transects of 80° E and equator (Table 4). Among environmental variables, the only significant positive correlation was observed between the S. punctatum and A. cineraria with temperature and L. dodecantha with Chl a concentration (p < 0.005) (Table 6). Differently, a negative correlation was observed for S. polysiphonia with phosphate at a high significant value (p < 0.001) and salinity, D. tetrathalamus and H. armatum-hostile, with silicate concentration (Table 6). The significant correlation between radiolarian species with SST and SSS also has been reported from the Pacific Ocean [66].

These were the major environmental factors affecting on the distribution of silicified dominant species including Collodarian in the EIO (Figure 15). This finding supports previous research in Australia [76], which suggested that the symbionts of such protists can take up nutrients from their hosts, so nutrients affecting their distribution in that area and so their presence at low latitudes and absence from high ones are to be expected. The symbiotic species of Collodarian were dominant in the study area, which is associated with high radiolarian densities in numerous studies of the tropical ocean [2,5]. The symbiotic radiolarians are key player for the carbon export from nutrient-poor and subtropical marine areas [2,10]. Recently, radiolarian is contributing 29% to the meso-zooplankton in the ocean, cause the high biomass 0.089 pg C in the upper column 200 m depth and about 5.2% of total carbon biomass [2]. Photosymbiotic Collodarian species often found with the high abundant at the 100 m depths and with Acantharia contributed only 0.18% primary production in the oligotrophic waters [2]. However, Polycystinea radiolarian at the 200–400 m constituted of high biomass to the total Rhizarian [4]. The energy obtained by their primary production can enhance carbon uptake and increase carbon fluxes due to the high rates of primary production (1400–41,000 ng carbohydrate colony⁻¹ h⁻¹) and result in high carbon fixation (91.16 µg) [77,78]. Earlier studies indicated that colonial radiolarians in the subsurface have a strong correlation with chlorophyll a maximum at the subsurface layers (25–75 m) and cause high fluxes in the Gulf of Mexico [76] and also in the Pacific Ocean [4].

Our studies suggest that gelatinous packed species belonging to Collodarian were abundant in the euphotic zone and dominate the radiolarian community. Many species occurring in the surface areas are particularly herbivorous and omnivorous planktons which can be indirectly related to chlorophyll a, similar to other studies [72,79]. Recently, bacteria [28], tintinnids [29], coccolithophores [30], phytoplankton studies [31] have been conducted at the same sampling area and these can be nutrient source of radiolarian [23]. The radiolarian depicting in the surface for feeding the chlorophyll a producing organism such as bacterial feeders or crustaceans or it is possible that the species may be required similar environment.

6. Conclusions The results presented here describe the first comprehensive plankton-net sampling of the Eastern Indian Ocean. Planktonic radiolarian within five classes, e.g., Acantharia, Taxopodida, Collodaria, Phaeodaria and Polycystinea, including 60 newly reported taxa, were described in this study area. Community composition, abundance and distribution of radiolaria are different in each transect but, in each case, are dominated by the combined occurrences of Collodaria and Taxopodida. The highest abundance was recorded in the south-equator region, and the horizontal distribution appears to be controlled by three explanatory variables; chlorophyll a, silicate, and nitrate dominate the radiolarian species:

(a) Based on the RDA analysis, the radiolarian abundance in the north–south transect at approximately 80° E (“Long-80”) is mainly controlled by the surface silicate concentrations or might be possible for high nitrate concentrations or dissolved oxygen in the mixed layer. So, their connection with the growth of radiolarian is casual and dominates only the S. punctatum, S. fuscum and S. zanclea species.

(b) Radiolarian abundance in an equatorial transect (“Lat-0”) is strongly related to the chlorophyll a level in the mixed layer which dominates the S. punctatum, L. pentagona pentagone and S. zanclea but it is not certain that chlorophyll a is necessary for radiolarian growth. Thus, the connection may be indirect, such as radiolarians may consume chlorophyll-producing organisms or may simply require a similar environment.

(c) The radiolarian abundance in a north-equator transect (“Long-90”) is controlled by the high salinities and silicate concentrations intrude via the BBR runoff flow from BOB and the open coastal zone of Sumatra, which dominates the S. multispina, C. inerme and S. polysiphonia, S. punctatum, C. sociliatorchia, S. zanclea and D. tetrathalamus.

NEQ-Long-90			SEQ-Long-80			EQ-Lat-0
Station	Longitude	Latitude	Station	Longitude	Latitude	Station	Longitude	Latitude
HFA01	6.502	88.002	I302	2.415	80.003	I101	0.485	92.685
HFA03	7.505	88.003	I304	1.502	80	I401	0.002	80.004
HFA03	8.5	88	I308	−0.997	80.006	I403	0	82.007
I101	0.485	92.685	I310	−1.996	80.004	I404	0	83.016
I103	1.5	92.233	I312	−2.999	80	I406	0	85.002
I105	2.499	91.632	I314	−4	80.003	I407	0	86.002
I107	3.5	91.03	I316	−5	80	I410	0	89.002
I109	4.5	90.434	I318	−6	80.002	I412	0	91.003
I501	−0.5	93.103	I401	0.002	80.004	I413	0	92.014
I503	−1.5	93.81	I807	−0.5	81.003	I501	−0.5	93.103
I505	−2.5	94.337	I809	−1.499	81.002	I807	−0.5	81.003
I507	−3.5	94.935	I811	−2.499	80.998
I509	−4.5	95.532	I813	−3.5	80.999
I511	−5.499	96.098	I815	−4.499	81
I601	6	89.5	I817	−5.5	81
I603	5.999	88.502
IQ2	10.078	88.133

Species	Frequency (fi)/%	Relative Abundance (P)/%	Dominance Index Y
Sphaerozoum punctatum	52	22.2	5.1
Sphaerozoum fuscum	50	15.8	3.49
Didymocyrtis tetrathalamus	52	3.5	0.8
Sticholonche zanclea	39	7.7	1.32
Siphonsphaera polysiphonia	32	3.3	0.46
Lophospyris pentagona pentagone	36	2.3	0.37
Eucytidium hexastichum	36	1.9	0.31
Anthocyrtidium cineraria	20	0.8	0.07
Pterocanium pretextum	25	0.8	0.09
Eucecryphalus gegenbauri	20	1.2	0.11
Lophophaena dodecantha	20	0.7	0.07
Arachnocorys circumtexta	23	1.2	0.12
Callimitra sociliatorchia	27	2	0.24
Amphirrhopalum ypsilon	20	0.7	0.06
Hexacontium armatum-hostile	23	1.4	0.14
Acanthsphaera actinote	23	0.7	0.07
Euchitonia elegans	18	0.6	0.05
Spongodiscus resurgens	18	0.6	0.05

Source	Dataset	Inertia	PC1	PC2	PC3	PC4	PC5	PC6	PC7	PC8	PC9	PC10
All	Hellinger	0.786	0.098	0.074	0.054	0.048	0.042	0.038	0.033	0.029	0.027	0.026
All	Dominant	0.786	0.124	0.099	0.081	0.069	0.050	0.048	0.045	0.033	0.030	0.026
Long-90	Hellinger	0.755	0.127	0.099	0.073	0.064	0.056	0.052	0.044	0.038	0.036	0.030
Long-90	Dominant	0.755	0.147	0.112	0.085	0.076	0.055	0.044	0.040	0.035	0.021	0.017
Long-80	Hellinger	0.774	0.164	0.108	0.095	0.078	0.065	0.055	0.052	0.040	0.033	0.028
Long-80	Dominant	0.774	0.185	0.124	0.102	0.077	0.064	0.047	0.023	0.020	0.016	0.008
Lat-0	Hellinger	0.820	0.154	0.123	0.111	0.093	0.090	0.078	0.067	0.046	0.033	0.023
Lat-0	Dominant	0.820	0.176	0.124	0.106	0.091	0.089	0.063	0.048	0.029	0.016	0.011

Transect	Explanatory Variable	Depth (m bsl)	Total Variance	Variance Explained	% Variance Explained	F Statistic	Pr(>F)	Λ1/Λ2
Lat-0	Chlorophyll a	0	0.8085	0.1399	17.3	1.884	0.0001	1.13
Lat-0	Chlorophyll a	25	0.8085	0.1312	16.2	1.743	0.0015	1.05
Long-80	Nitrate	75	0.7939	0.1020	12.9	1.917	0.0013	1.03
Long-80	Diss. Oxygen	25	0.7939	0.1077	13.6	2.041	0.0006	1.06
Long-80	Diss. Oxygen	50	0.7939	0.1077	13.6	2.042	0.0008	1.09
Long-80	Silicate	0	0.7939	0.1065	13.4	2.014	0.0003	1.06

Correlation Coefficients (43 Degrees of Freedom)	Nitrate µM	Phosphate µM	Salinity psu	Silicate µM	Temperature °C	Dissolved Oxygen mg/L
Chlorophyll a µg/L	0.131	−0.326	−0.164	−0.126	0.149	−0.124
Nitrate µM		−0.540	−0.339	−0.556	0.280	0.129
Phosphate µM			0.430	0.662	−0.465	−0.159
Salinity psu				0.672	−0.551	−0.406
Silicate µM					−0.676	−0.622
Temperature °C						0.590
Dissolved Oxygen mg/L

Figure in bold italics are for those combinations significant at the 1% level.

Dominant Species	Temperature °C	Salinity psu	Chlorophyll a µg/L	Phosphate µM	Silicate µM	Nitrate µM
Sticholonche zanclea	0.118	0.130	0.208	0.031	0.121	−0.020
Siphonosphaera polysiphonia	0.177	−0.471**	−0.213	0.004*	0.130	0.068
Sphaerozoum punctatum	−0.312*	−0.043	−0.187	0.085	0.092	−0.048
Sphaerozoum fuscum	0.007	0.117	0.027	−0.014	−0.237	−0.097
Acanthosphaera actinota	−0.062	0.089	−0.094	0.059	−0.100	0.055
Astrosphaera hexagonalis	−0.028	0.065	0.136	−0.058	−0.291	−0.191
Hexacontium armatum-hostile	0.051	0.164	−0.015	0.066	−0.318*	−0.272
Didymocyrtis tetrathalamus	0.102	0.290	0.226	−0.290	−0.158	−0.070
Amphirhopalum ypsilon	−0.058	0.036	−0.087	0.054	0.019	−0.066
Euchitonia elegans	0.104	0.071	0.117	0.032	−0.035	−0.111
Spongodiscus resurgens	0.261	0.207	−0.094	−0.100	−0.022	−0.119
Arachnocorys circumtexta	−0.132	−0.200	−0.067	−0.028	0.117	0.195
Callimitra sociliatorchia	0.030	−0.123	−0.131	−0.018	0.178	0.061
Lophospyris pentagona pentagone	−0.026	0.061	0.100*	−0.045	0.045*	0.062
Lophophaena dodecantha	−0.066	0.099	0.360	−0.063	−0.309	−0.163
Anthocyrtidium cineraria	−0.308*	0.154	−0.098	−0.196	−0.050	−0.049
Pterocanium pretextum	−0.029	0.054	−0.109	−0.060	0.081	−0.103
Eucytidium hexastichum	−0.013	0.285	0.171	−0.123	−0.056	−0.275
Eucecryphalus gegenbauri	−0.080	0.043	−0.158	−0.077	−0.043	0.207

Figure in bold italics are for those combinations significant at the * p < 0.05 and ** p <0.01 level.

Class	Taxa	Worldwide Distribution Records	References
Acantharia
	Acanthochiasma fusiforme	Mediterranean Sea, Atlantic and Pacific Ocean	[54]
	Acanthostaurus conacanthus	South Atlantic Ocean	[54]
	Trizona brandti	Mediterranean Sea, the Red Sea, the English Channel and the West Pacific Ocean	[38]
	Dictyacantha tabulate	Pacific Ocean	[38]
	Amphibelone anomala	East China Sea, Australia	[38,55]
	Amphilonche concreta	Mediterranean Sea, Atlantic and Pacific Ocean	[38]
	Diploconus cylindrus	North Atlantic Ocean, Pacific Ocean	[14,38]
	Diploconus faces	North Pacific Ocean, southward Pacific, Atlantic and Indian Ocean	[1]
Taxopodida
	Sticholonche zanclea	East China Sea; Mediterranean Sea	[56,57]
Collodaria	Solenosphaera zanguebarica	Atlantic, Pacific Ocean, eastern Indian Ocean, China Sea	[25,58]
	Collosphaera huxleyi	Atlantic, Pacific Ocean	[25,58]
	Collosphaera macropora	China Sea	[25,58]
	Collosphaera tuberosa	China Sea	[25,58]
	Siphonosphaera polysiphonia	Pacific, Atlantic Ocean	[25,58]
	Siphonosphaera socialis	Gulf of California, Atlantic Ocean	[39]
	Siphonosphera magnisphera	Atlantic Ocean	[39]
	Sphaerozoum punctatum	Gulf of California; Southern Indian Ocean	[27]
	Sphaerozoum cf. ovodimare	Mediterranean Sea, Naples, Messina, Atlantic, Canary Islands, Cape Verde Islands, West Coast of Africa	[59]
	Sphaerozoum fuscum	East China Sea	[60]
	Collozoum inerme	Northern Norwegian Sea, the South Atlantic Ocean	[59]
	Thalassoxanthium cervicorne	Atlantic Ocean	[14]
	Thalassoxanthium octoceras	Madagascar, Rabbe Island, Indian Ocean	[38]
Phaeodaria	Aularia ternaria	North Pacific Ocean, south Atlantic and Pacific Ocean	[39]
	Aulacantha scolymantha	Pacific, Atlantic, Indian and Mediterranean Sea	[14,61]
	Auloceros arborescens subelegens	Atlantic and Arctic Ocean	[14]
	Coelodendrum ramosissimum	China Sea, Eastern Mediterranean Sea	[58]
	Castanella longispinium	South Atlantic Ocean, Gulf of Oman	[62]
	Conchopsis compressa	North Pacific (Japan and San Francisco), South Pacific Ocean	[14]
	Conchellium capsula	Atlantic Ocean and Pacific Ocean	[14]
	Conchidium caudatum	Atlantic Ocean	[14]
	Challengeron radians	Pacific, Atlantic Ocean and Japan Sea	[39]
	Pharyngella gastrula	Pacific and Atlantic Ocean	[39]
Spumellaria	Acanthosphaera dodecastyla	Atlantic	[39]
	Acanthosphaera pinchuda	Equatorial Pacific Ocean	[13]
	Actinosphaera capillacea	Pacific and Atlantic Ocean; East China Sea	[21,39]
	Astrosphaera hexagonalis	Pacific and Atlantic Ocean	[39]
	Centrocubus cladostylus	Gulf of California, Pacific Ocean	[1]
	Cladococcus scoparius	Central Pacific Ocean	[38]
	Cyrtidosphaera reticulata	Japan	[63]
	Spongurus pylomaticus	Gulf of California	[64]
	Streblacantha circumtexta	Arctic Sea, Nordic Sea; North Atlantic Ocean	[16,65]
Nassellaria	Archibursa tripodiscus	Central Pacific Ocean	[38]
	Callimitra solocicribrata	Central Pacific Ocean	[38]
	Cephalospyris cancellata	Central Pacific Ocean	[38]
	Cladoscenium limbatum	Nordic Seas	[39]
	Conarachnium parabolicum	Pacific and Atlantic Ocean	[39]
	Dictyocodon palladius	Central Pacific Ocean	[59]
	Eucecryphalus clinatus	Gulf of California	[39]
	Lampromitra schultzei	Pacific and Atlantic Ocean	[14,39]
	Litharachnium tentorium	Mediterranean Sea	[38]
	Pterocanium korotnevi	Central Equatorial Pacific Ocean	[1]
	Pteroscenium pinnatum	Pacific and Atlantic Ocean	[13,38]
	Pseudodictyophimus gracilipes	Arctic Ocean; Gulf of California; Nordic Sea; North Atlantic; North Pacific Ocean	[63,64]
	Phormacantha hystrix	Arctic Ocean, Gulf of California, Nordic Sea	[64,65]
	Stichopilium bicorne	Gulf of California,	[64]
	Tetraplecta pinigera	Central Pacific Ocean	[38]
	Theocorys veneris	Pacific and Atlantic Ocean	[39]

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/2073-4441/12/12/3502/s1. Table S1: The raw data as (ind/m-3) of all radiolarian taxa that were recorded in Eastern Indian Ocean cruise 'Shiyan I' in 2014. Figure S1. The Satellite Chlorophyll a data profile derived from NASA during April 2014 (a) and May 2014 (b) in the Eastern Indian Ocean.

Author Contributions

J.S. designed the cruise study for community analysis in EIO. C.D. collected the samples on board during cruise "Shiyan I" in 2014. S.M. did all the laboratory work including microscopy examination for identification, counting and wrote the manuscript. J.R. helped to modify all statistical analyses and drafted the statistical portion of the manuscript. X.Z. did chlorophyll a extraction. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Nature Science Foundation of China, grants # (41876134, 41676112 and 41276124), the Tianjin 131 Innovation Team Program (20180314), and the Changjiang Scholar Program of Chinese Ministry of Education (T2014253) to Jun Sun. The study was also partly funded by the National Science Foundation project (41406155) to Changling Ding.

Acknowledgments

We would also like to thank the Open Cruise Project in Eastern Indian Ocean of National Nature Science Foundation of China (NORC2014-10) for sharing their ship time, which carried out by the RV "Shiyan I". We are particularly grateful for the corrections and editing by O. Roger Anderson (Sr. Research scientist, Lamont-Doherty Earth Observatory, Columbia University).

Conflicts of Interest

The authors declare no conflict of interest.