Free-living marine nematodes community structure in the conservation area (Chaojing Park) and its adjacent area of Keelung, Taiwan

Abstract

Studies conducted in the same seas or even study sites nearby each other, showed very different community structure, implying the patchiness of free-living marine nematodes which may be related to the sedimentary environment such as sediment type and food availability of the study area. This study was motivated by the concerns about the impacts of high level of anthropogenic activities on Chaojing Park (gazetted as Wanghaixiang’s Chao-Jing Bay Resource Conservation Area (WCJBRA) in 2016). The present study provides baseline knowledge of free-living marine nematode community structure in WCJBRA and identify potential marine nematodes as bioindicators to indicate possible impacts of the anthropogenic activities to the Chaojing Park. A total of 15 stations were selected in the subtidal zones of WCJBRA and its adjacent area. Marine nematode sample collection was carried out on the 13th and 14th of September 2019 using SCUBA diving technique. Results showed positive correlation between nematode density and medium sand (500μm-1.0mm). Presence of certain species such as Daptonema sp., Pomponema sp. and Innocuonema sp. indicates presence of disturbances in S12 and S13. Several species also showed potential to be introduced as indicator for healthy environment subjected to further studies on nematode-pollutants relationship, particularly on autecology as per se. Higher species diversity, H’ index of S1-S8 and S11 was categorised as Good Condition; followed by station with moderate species diversity index (S9, S10, S14 –Moderate Condition) zone; and lastly S12, S13 and S15 (Poor Condition).

Introduction

Free-living marine nematodes exist as epifauna (on the sediment) or infauna (in the sediment) [1]. They comprise a sizable percentage of marine benthic populations, accounting for four out of every five bottom-dwelling metazoans on the world [2]. One of its ecological functions is to regulate nitrogen remineralization [3] by grazing on decomposers and generating ammonium for use by plants and bacteria [1,4]. Additionally, they are capable of immobilising carbon inside their body structure and re-mineralizing it during energy production via respiration or decomposition [1]. Additionally, marine nematodes contribute to the marine food cycle by acting as prey for predators [1,3,4].

The majority of studies on free-living marine nematodes were ecological investigations [5–7] or environmental pollution impact studies [7,8]. These are due to their close association with sediment, short life cycle, high reproduction rate, abundance, and small size [9–12], allowing them to respond quickly and early in the history of the contamination site [7,13–15]. As a result, investigations conducted in the same oceans or even close study sites revealed highly diverse community structures, emphasising the patchiness of free-living marine nematodes that may be related to the sedimentary environment, such as sediment type and food availability of the study area. For example, Liu [16] found 232 species and 149 genera in the community structure of nematode in the southern Yellow Sea. In the southern Taiwan Straits, there are 75 genera of free-living marine nematodes [17] and 101 species (58 genera) in Hong Kong waters [18].

Keelung City Government designated Chaojing as Wanghaixiang’s Chao-Jing Bay Resource Conservation Area (WCJBRA) in 2016. The conserved area was to ensure the survival of coral reefs as well as to address overfishing issues [19]. Although the bottom substrate was mainly comprise of rock, siltation was reported in Chaojing area, suffocating the coral reefs [20]. Moreover, Lin [20] also stated that sea urchin and lobsters which are indicator organisms for overfishing in the Reef Check Survey were rarely spotted in the area. As a result, significant concerns had been raised about the anthropogenic activity in the bay, as well as the presence of a commercial fishing harbour right close to the conservation area. The discharge of domestic garbage from surrounding restaurants, houses, boats, and shipping activities at the commercial fishing port, as well as antifouling paint activities in the harbour, may have an impact on the marine protected area’s environmental health. Few studies on free-living marine nematodes in relation to the impact of anthropogenic disturbance from fishing harbours have been conducted [21–23], and they have revealed that the meiofauna community responds well to anthropogenic disturbances.

Although free-living marine nematodes have shown promise in quantifying environmental status [7,8], it is unfortunate that little research has been conducted on the community structure of free-living marine nematodes in Taiwan, especially in marine protected area. The most recent study conducted was in northwest Taiwan on the intertidal zone of Danshuei River estuary [24]. Thus, this is the first study of its kind undertaken in Taiwan on the community structure of free-living marine nematodes at WCJBRA and its adjacent area, providing a solid foundation for future research. The aims of this study are to (1) characterise the community structure of free-living marine nematodes in the WCJBRA and nearby area, and (2) determine the physico-chemical environmental variables that influence the community structure of marine nematodes in this area.

Methodology

On the 13th and 14th of September 2019, this project was carried out in the subtidal zones of the marine conservation area and its surrounding region in Chaojing Park, Keelung, Taiwan. There was no historical environmental data in this location, but anthropogenic activity such as a fishing harbour, restaurants, antifouling coating, and a housing community were seen. A total of 15 sampling stations were chosen, and the geographical coordinates for each are shown below in Table 1. In Fig 1, the WCJBRA was denoted by a red dashed line, whereas the fishing harbour was at S14.

[Figure omitted. See PDF.]

Map illustrated by digitizing the base map of OpenStreetMap (under Open Database Lisence). Details refer to www.opendatacommons.org/licenses/odbl).

[Figure omitted. See PDF.]

The collection of marine nematode samples was done using the subtidal zone quantitative approach provided by Somerfield and Warwick [25]. Subtidal sediments have been gathered via SCUBA diving directly. A modified corer (11.95cm²) was used to gather triplicate samples. It was inserted vertically into sediments around 5cm deep. Additional sediment samples for particle size analysis have been obtained (PSA). Physico-chemical water quality parameters of bottom waters including salinity, temperature, dissolve oxygen (DO) and pH were measured in situ using a YSI multiparameter probe (ProPlus).

Sample processing

Particle size samples were dried in an oven at 80°C for 24 hours to remove moist. Each sample weighing 100g was transferred to a beaker. Hydrogen peroxide were added to remove organic matter and left to stand overnight with sodium hexametaphosphate prior to dry sieving [26] to determine the particle size of the sediment from each station. Sediment retained on each sieve was weighed. Size of the sediments was graded according to the size scale suggested by [27].

Sediments collected for nematode analysis were sieved through a 500 μm mesh and decanted over a 25 μm mesh (Retch, ASTM E11) with tap water to remove macrofauna and microfauna prior to preservation for future extraction [23,28]. Sediments were washed into sample bottles and then preserved with 5% formalin [11]. Rose Bengal (0.5 g L^-1) was applied and allowed to dye the nematode specimens thoroughly overnight. Sediments were rinsed and transferred to a centrifugation tube, which was filled to the 80% mark with Ludox and centrifuged at a rate of 3000 rpm for 5 minutes at 4°C. [29]. Supernatants were preserved in 5% formalin [23] for further sorting, counting, and identifying. The nematodes were isolated using a wire loop and dehydrated in ethanol-glycerol solution according to [30] to create permanent mounts for identification under a compound microscope using the nematode pictorial keys [31–36]. Additionally, marine nematodes have been divided into four trophic functional feeding groups (FFGs) based on their buccal cavity morphology: selective deposit feeders are designated as 1A; non-selective deposit feeders are designated as 1B; epigrowth feeders are designated as 2A; and predators/omnivores are designated as 2B [37–40].

Data and statistical analyses

PRIMER v6 and SPSS (Statistical package for Social Science) v27 software were used for the statistical analyses in present study. DIVERSE analysis was used to determine the biological indices including average number of genera per sample (S), Total number of individuals (N), Shannon-Weiner Index (H’), Pielou’s Evenness Index (J’) and Margalef’s Index (d) from each station. The Maturity Index (MI) [41,42] was calculated as the weighted average of the individual colonizer-persistent (c-p) values to identify sites under stress. Threshold for nematode descriptors can be referred in Moreno’s study [13].Analysis of Similarity (ANOSIM) was used to determine the similarity of nematode communities between stations while Similarity Percentage (SIMPER) was used to identify the similarity percentage within and among each station. Multidimensional Scaling (MDS) with Bray-Curtis Similarity matric was used to illustrate the level of similarity of nematode species composition between stations. The data was square root transformed prior to the analysis. Nematode genera were classified according to Wieser [37] into four feeding groups: selective feeder (1A), nonselective feeder (1B), deposit feeders or epistrate feeders (2A), and predators or omnivores (2B) to investigate the trophic structure of the assemblages. CLUSTER analysis subjected to Euclidean distance was used to show the similarity of environmental parameters between stations. Pearson Correlation between environmental parameters and nematode density were conducted.

Results

Nematode composition and density

Across all sample sites in this research, a total of 111 nematode species from 25 families were identified and documented. DIVERSE analysis found that the northern study sites (S1–S10) had a higher biological indices value, while the southern research sites had a lower value (S11-S15). The maximum nematode density was found in S10 (N = 363 ± 91.65 ind./10cm²), while the lowest was found in S13 (N = 11 ± 8.89 ind./10cm²), followed by S12 (N = 15 ± 4.16 ind./10cm²). Stations situated inside WCJBRA, such as S2 (N = 115 ± 19.86 ind./10cm²) and S3 (N = 92 ± 17.95 ind./10cm²), had significantly higher nematode densities than other stations, but S1 and S4 had significantly lower densities (N = 50 ± 11.5 ind./10cm² and N = 43 ± 4.36 ind./10cm², respectively). On the other hand, S2 has the third largest total number of individuals (N), the most species (S = 31 ± 3.51), the most species richness (d = 6.19 ± 0.53), the highest species diversity (H’ = 3.11 ± 0.07%), and the second highest species evenness (J’ = 0.91 ± 0.01%). Additionally, S13 had the fewest species and individuals (S = 3 ± 2.08%; N = 11 ± 8.89 individuals/10cm²), the lowest species richness and diversity (d = 0.73 ± 0.25%; H’ = 0.79 ± 0.29%), but moderate high in species evenness (J’ = 0.87 ± 0.11%). In general, the nematode population in this research was dispersed uniformly (J’>0.7), except for S15 (J’ = 0.50 ± 0.33). Additionally, S15 has the second-lowest species richness and lowest diversity of species (d = 0.97 ± 0.30; H’ = 0.79 ±0.51%).

Nematode assemblages

S4 had the lowest MI value (1.59 ± 0.05) and S13 had the highest (4.34 ± 0.20), followed by S15 (3.17 ± 0.22). Despite its location inside the WCJBRA, S4 was classified as having a bad environmental state using Moreno’s (13) standards. On the other hand, S14, the location of the commercial fishing port, had excellent environmental conditions. C-p 5 (extreme persisters) was not detected in all stations, while c-p 1 (extreme colonisers) was detected only in S12 (6.06%), which was dominated by c-p 2 (45.65%) and c-p 3. (37.10%). C-p 2 (general colonisers) and C-p 3 were typically dominant in all studied areas except for S13, S14 and S15 (Table 2). The lowest percentages of c-p 2 (11.85%) and c-p 3 (9.59%) were identified in S13 and S15 respectively. S14 were dominated by c-p 4 (53.55%) with lowest percentage of c-p 2 at 12.57%. The highest c-p 4 value were 58.55% in S15 (provided by Eurystomina sp.2). S4 had the greatest concentration of c-p 2 at 73.94%, which was mostly attributed to Spilophorella sp. (27.62%) and Axonolaimus sp. (21.23%). Whilst, the highest concentration of c-p 3 were identified in S8 at 59.98%, contributed by Parapinnanema sp. (13.32%) and Innocuonema sp. (10.81%).

[Figure omitted. See PDF.]

The results of one-way ANOSIM revealed a 0.1% significant difference among groups with a R value of 0.748. SIMPER analysis showed that the marine nematode population in S10 had the highest percentage of similarity (58.47%), followed by S2 (55.44%) and then S14 (54.56%). Similar percentage in S10 was contributed by Spirinia sp. (26.6%), Microlaimus sp. (24.82%), and Pomponema sp. (10.24%) while S2 was contributed by Axonolaimus sp., Metadesmodora sp. and Meyersia sp. (9.95%, 6.91% and 6.63% respectively). On the other hand, the lowest similarity percentage within station was S11 for 20.45%. Most of the stations were highly dissimilar with S12, S13, and S15 in the community structure. In general, S12 has the highest dissimilarity percentage with S3 (96.82%), contributed by Spilophorella sp (12.23%)., Draconema sp. (11.23%), and Acticonema sp. (6.03%) whereas S13 has the highest dissimilarity with S1 (97.30%) on Spillophorella sp., Viscosia sp., and Chromadorita sp., (5.98%, 5.6% and 5.56% respectively). Besides that, Eurystomina sp.2 (16.06%), Innocuonema sp. (5.26%), and Axonolaimus sp. (5.18%) contributed the highest dissimilarity percentage between S7 and S15 with average dissimilarity of 97.95%. Eurystomina sp. 2 was recorded to be the dominant species in S15 was noted to have high dissimilarity percentage with other stations. Furthermore, S14 has the highest average dissimilarity with S15 (96.40%) due to the present of Eurystomina sp. (21.78%) and Axonolaimus sp. (7.8%) in S15, while Bathyeurystomina sp. with contribution of 12.14% in S14.

The result of Multi-dimensional Scaling (MDS) on the species composition after square root transformation and Bray-Curtis similarity matrices revealed that S1 and S2 were grouped together at 40% similarity, followed by S5, S6, S7, S8, and S11, and S9 and S10. The remaining six stations were grouped independently as an independent group (Fig 2). Nonetheless, when the similarity threshold was raised to 60%, all stations formed autonomous groups. In general, the findings indicated a low degree of resemblance between the marine nematode communities at the research sites. Each cluster has its own distinct mix of nematode species.

[Figure omitted. See PDF.]

Fig 3 shows the percentage of FFG in each station where most of the stations were dominated by 2A except S12, S14 and S15. In fact, S12 was dominated by 1B (55.00%) while both S14 and S15 were comprised mostly of 2B (S14: 54.17%; S15: 86.50%). On the other hand, the highest percentage of 2A was found in S13 at 76.92%.

[Figure omitted. See PDF.]

Sedimentary characteristics and Pearson correlation

During the current investigation, a total of four in situ physico-chemical water parameters were collected, as well as sediment for granulometry analysis (Table 1). The similarity of environmental factors across stations was shown by CLUSTER analysis using Euclidean distance (Fig 4). The stations were grouped into four groups when the similarity distance was set at 4. Station S9 and S10 were in one group, S11 and S12 were in another, S14 was in a separate group, and the rest of the stations were in a single group. The temperatures in S9 (25.17±0.26°C) and S10 (25.33±0.52°C) were slightly higher than in other stations, while the granulometry analysis showed S9 and S10 had coarser sediment environments. S9 was dominated by coarse sand (500μm-1.0mm) at 47.08% while S10 comprised mostly medium sand (250μm-500μm) at 37.25%. On the other hand, S11 and S12 were grouped together due to their extremely low salinity (S11: 16.55±0.45ppt; S12: 16.89±0.08ppt). In addition, S14 was in an independent group because of its lowest pH (6.37±0.10), highest DO (5.93±0.05mg/L) and highest silt percentage (<63μm) at 27.79% among other stations. Results of Pearson correlation showed that nematode density was positively correlated with medium sand (r = 0.593, p = 0.02).

[Figure omitted. See PDF.]

Discussion

Numerous studies have established correlations between environmental variables and the formation of nematode communities [43–46]. According to the majority of research, sediment grain size is the fundamental factor influencing the abundance and variety of marine nematodes [1,47,48]. The current investigation established a positive correlation between nematode density and medium sand. S9 and S10 were classified separately in CLUSTER analysis due to their coarser sediment conditions, with S9 being predominantly coarse sand and S10 being predominantly medium sand with the highest nematode density. The finer the grain size, the greater the abundance of marine nematodes, and vice versa, since the finer grain size gives a greater total surface area for organic matter adsorption [49] and so supports higher nematode community abundances [45,50].

Although S10 had the highest nematode density, S2 in the WCJBRCA had the highest number of species, species richness, and species diversity. The Shannon-Wiener diversity index (H’), which is extensively employed to quantify benthic community diversity, may also indicate sediment condition [51]. Indexes of species diversity with a value less than 2.0 indicate low species diversity and poor sediment condition, whereas those with a value greater than 3.0 indicate the reverse [52]. In general, stations located northeast of the commercial fishing harbour (S14) have a more than 2.0 species diversity value, whereas those located southeast of the commercial fishing harbour have a value less than 2.0. On the other hand, greater nematode species diversity may suggest a diversity of food supplies, whilst lower indices may indicate a scarcity or lack of food sources diversity [45,53].

S15 and S13 had the lowest sediment condition (diversity index < 2). According to residents, S13 was once a port and is now occupied with a few residential units. As a result, local inputs such as sewage or waste disposal were anticipated. Domestic sewage may contribute to localised eutrophication [54], altering the composition of nematode assemblages by organic matter enrichment and being related with pollution [55], resulting in a lack of nematode species diversity. According to Urban-Malinga [50], species diversity often decreases with pollution, while Sahraean [56] showed decreased nematode diversity and abundance in sewage-impacted areas of the city or residential complex. On the other hand, the low diversity index in S15 could be a result of the rainy weather, along with high wave conditions during sampling. Currents caused by the wind can suspend nematodes from sediment into the water column until they reach still water [57]. Additionally, Venekey [58] indicate that meiofauna may burrow down to ten centimetres into the sediment to avoid being suspended in the water column. Moreover, vertical movement of nematodes to deeper strata is conceivable in well-aerated habitats with sediment sizes greater than 125 m (S15 were primarily composed of very coarse sand, > 1 mm) and to mitigate the possibility of erosion caused by high flow velocity [59]. However, upward migration of some species may occur as a result of resource allocation and quality [60], possibly contribute to the result of certain species domination in S15.

Additionally, S10 and S14 exhibit poor sediment condition as measured by nematode diversity values, but S9, which is located between S10 and S14, exhibits moderate sediment condition. It is believed that discharge originates at S14 and collects at S10 as a result of the breakwater (Fig 1). The breakwater may have a substantial effect on the area’s currents and wave energy [61], perhaps depositing anthropogenic contaminants in S10 and explaining the moderate sediment condition in S11. Petroleum, sewage, and biomass combustion are only a few of the significant toxins found in fishing harbours that have an effect on the nematode community structure [62,63]. Additionally, elevated PAH levels in sediments suggest a high organic matter content [64], which influences the composition of marine nematodes [53,56,65].

Maturity Index was initially proposed for the study of terrestrial and freshwater habitats, and was then extended to marine and brackish ecosystems [41,66], yet it was being used conservatively [67]. Bongers [42] distinguished colonizers or c-p 1 which are more tolerant to environmental variations while presisters or c-p 5 which are more sensitive, thus, increase in number of persisters leads to an increase in MI value. Despites that it could be a good descriptor of stress as it is influenced lesser by natural variables (e.g. sediment granulometry) than other available indices such as H’ and J, there were cases with unsatisfactory results [13,68,69]. In current study, MI values showed significant contrasting results with H’ especially in S15 which was dominant by Eurystomina sp. 2 (c-p 4). Semprucci et al. [70] reported only a slight influence of granulometry on MI and c-p, thus, giving more support to the hypothesis of higher anthropogenic activity in the area. Furthermore, persisters nematodes in high abundance may either indicates the prolonged poor conditions, leading to nematode assemblage which is well adapted to the unfavourable conditions [71,72], or its ability to profit from the organisms killed by the extreme environmental conditions [13]. Therefore, this study further supports the application of H’ particularly in marine environment and H’ has been selected for the EcoQ class definition was selected for the EcoQ class definition due to its wide application [67].

In general, the MDS results indicated a low degree of similarity across the marine nematode community studied in this study. Coral reefs dominated the first category, which included S1 and S2 in the WCJBRA. Previous research has shown that the coral reef area contains a greater number of marine nematodes, either in terms of species or densities [73], in comparison to areas surrounding those subjected to anthropogenic activities [33,34], such as S14. Spilophorella sp., Viscosia sp. and Meyersia sp. were the top three nematode genera discovered in S1. Marine nematode, Spilophorella sp., was one of the dominant species observed in the study site of Semprucci [15] which has high H’ value (3.5–4.5), indicating good sediment quality. Although S1 have fairly high species diversity index, yet the index value was slightly below 3.0 which indicated the sediment quality still under a satisfactory level due to slight disturbed by the ecotourism activities in the WCJRCA. Several studies had noted the impact of tourism and diving activities on coral ecosystem [74–78]. Inexperienced divers with low buoyancy control skill have higher physical contact rate with coral reefs compare to experienced divers causing coral breakage [79]. This scenario not only will reduce the coral biomass but also cause the injured coral to suffer from slower growth rates and become more susceptible to disease and predation [80], altering the reefs state into a less structural complex state as well as impacting the benthic assemblage [81,82], thus, the lower species diversity index in S1.

In general, for most current sites, the community structure of marine nematodes was considerably different from S12, S13 and S15. The three stations were positioned most far from WCJRCA, indicating possible substantial anthropogenic disorders in household, port operations and exposure to the hydrodynamic processes. Previous studies have shown that several species of free-living marine nematodes may be used to identify pollutant levels and habitat or sediment quality [15,23,83], different environmental impacts [28,84] as well as depth of sediment collected [45]. Contaminated site evaluation normally incorporates water and sediment assessments to evaluate the overall level of contamination, but often does not indicate the toxicity of contaminants to biota [34,85]. The study by Fischer [86] shows that species of nematode showed ability and adaptive mechanisms to handle shifting environments to reduce their detrimental impacts. Some nematodes (for example, Enoplus brevis and Enoplus communis) have shown that their unique abilities to contains pollutants like Cd, Cu, Pb and Zn in their biotope [21].

The top three nematode abundances in S12 were Daptonema sp., Pomponema sp. and Cyatholaimus sp., while in S13 were Chromadorita sp., and Innocuonema sp. In previous studies, Daptonema sp., Sabatieria sp., Terschellingia sp., Theristus sp., and Paradontophora sp., were reported to be able to indicate pollution and environmental disturbances such as heavy metals [11,22,47,63], while Microlaimus sp. are known to inhabit anoxic environments in the study of Nanajkar [87]. Daptonema sp. was also proven by Mahmoudi [63] that this species showed a positive relationship with the concentration of diesel while this species was categorized as ‘opportunistic’ nematode species and elaborate their position in marine environmental biomonitoring. Thus, the presence of Daptonema sp. not only suggests potential low-level diesel contamination in that area but also shows that there are minimal levels of metal pollution [88]. Nematode species, Innocuonema sp., were found dominant in S13 and are known for their ability to tolerate different levels and types of disturbances in the study carried out by Semprucci [89]. Although the current study collected limited information on water quality parameters, previous studies had well documented the use of marine nematodes as bioindicators that could potentially be used as guidelines to estimate the potential impact of their presence in the study area.

Both Terschellingia sp. and Pomponema sp. were recorded as dominant species from S14. The former species was reported by Sahraeian [90] at study site with high polychlorinated biphenyls (PCBs) levels, a result of zinc production smelting and electrical power generation [91]. Terschellingia sp. also has a positive correlation with heavy metals such as Cadmium, Colbat, Chromium, Copper, Iron, Manganese, Nickel, Vadium, and Aluminium [92]. In addition, it has been linked to organic enrichment and as an indicator of poor ecological status due to its tolerance to pollution [13,93]. The latter species was documented in a previous study by Mahmoudi [63] as Pomponema sp. was recorded to be significantly affected by diesel contamination but not eliminated, thus, it is considered as the diesel-sensitive nematode. On the other hand, S15 was dominated by Eurystomina sp. with a minority of Axonolaimus sp. Previously, the presence of Axonolaimus sp. may potentially indicate the occurrence of anthropogenic disturbance in S15 as discussed by previous research [11,63,88,94] it is commonly found in areas with metal pollution and osmoregulatory stress on account of salinity values beyond the optimum range in that area.

Furthermore, different types of organic matter allow different nematode genera and species to coexist with distinct functional feeding groups (FFG) [39,40,45]. In general, epigrowth feeder (2A) dominated FFG in our investigation, implying that benthic microalgae were the most important food source coincided with the study by Leduc and Probert in 2011. On the other hand, S15 was discovered to be dominated by predators/omnivores (2B), which were made up of Eurystomina sp. 2. The presence of huge volumes of specific food kinds or nutrients, such as domestic discharge from neighbouring families, restaurants, or fishing vessels that potentially accumulate in the area that favours that specific group, could suggest the dominance of a given FFG [43,47].

The population structure of free-living marine nematodes in the WCJBRA and its surrounding areas was identified in this study. Daptonema sp., Innocuonema sp., Axonolaimus sp., and Pomponema sp. were discovered as nematode genera that could be used as bioindicators in the area, however more research is needed to undertake on water and sediment analyses. All the study stations were divided into three groups based on Shannon-Weiner diversity index, nematode species composition, and FFG distribution. The first group (good condition) included S1 to S8 and S11; the second group (moderate condition) included S9, S10, and S14; and the third group (poor condition) included S12, S13, and S15.

Conclusion

The present study only showed a positive correlation between nematode density and medium sand. This indicated the presence of confounding factors in characterising the nematode community in WCJBRA and its adjacent area. Findings on the species diversity index, together with the presence of potential bioindicator species, indicate that some stations are probably being affected by anthropogenic activities regardless of the level of disturbances. Certain genera such as Spilophorella sp. in the stations of GC in the present study have the potential to be introduced as indicators for a healthy environment. The presence of Daptonema sp., Innocuonema sp., Axonolaimus sp., and Pomponema sp., in S12, S13, and S15 indicated potential metal, diesel and organic contaminations corresponded to the anthropogenic activities in the area (domestic waste, diesel, antifouling coating and etc). The high abundance of Eurystomina sp. in S15 may be a potential bioindicator and symbolises the presence of a particular contaminant which was not measured in this study. Although data on contaminants are lacking, this research enables there to be a preliminary evaluation of the current ecological conditions in WCJBRA. It may also provide a baseline for the future monitoring in WCJBRA and for assessing long-term changes therein. Therefore, further studies on the detailed water quality profiling and nematode-pollution relationship are suggested, particularly on autecology as per se for future application in environment monitoring.

Acknowledgments

The current project is part of the research activity under the Memorandum of Understanding between Universiti Malaysia Sabah and National Museum of Marine Science and Technology, Taiwan facilitated this research collaboration. The authors would also like to thank Borneo Marine Research Institute for providing the facilities and for administrative and logistic support during fieldwork.

Citation: Ng W-L, Chen C-A, Mustafa S, Soo C-L, Liao Y-C, Shih T-W (2022) Free-living marine nematodes community structure in the conservation area (Chaojing Park) and its adjacent area of Keelung, Taiwan. PLoS ONE 17(5): e0268691. https://doi.org/10.1371/journal.pone.0268691

About the Authors:

Wei-Ling Ng

Roles: Formal analysis, Writing – original draft

Affiliation: Borneo Marine Research Institute, Universiti Malaysia Sabah, Kota Kinabalu, Sabah, Malaysia

Cheng-Ann Chen

Roles: Conceptualization, Investigation, Project administration, Resources, Software, Supervision, Writing – review & editing

E-mail: [email protected]

Affiliation: Borneo Marine Research Institute, Universiti Malaysia Sabah, Kota Kinabalu, Sabah, Malaysia

https://orcid.org/0000-0001-9168-1910

Saleem Mustafa

Roles: Conceptualization

Affiliation: Borneo Marine Research Institute, Universiti Malaysia Sabah, Kota Kinabalu, Sabah, Malaysia

Chen-Lin Soo

Roles: Investigation

Affiliation: Institute for Tropical Biology and Conservation, Universiti Malaysia Sabah, Kota Kinabalu, Sabah, Malaysia

Yun-Chih Liao

Roles: Conceptualization, Resources

Affiliation: Department of Earth and Life Science, University of Taipei, Taipei, Taiwan (R.O.C.)

Tung-Wei Shih

Roles: Conceptualization

Affiliation: National Museum of Marine Science & Technology, Keelung, Taiwan (R.O.C.)

References

1. Coull BC, Chander G. Meiobenthos. Nat Environ Res Counc. 2001;3: 1705–1711.

2. Bongers T, Ferris H. Nematode community structure as a bioindicator in environmental monitoring. Trends Ecol Evol. 1999;14: 224–228. pmid:10354624

3. Ferris H. Contribution of Nematodes to the Structure and Function of the Soil Food Web. J Nematol. 2010;42: 63–67. pmid:22736838

4. Neher DA. Role of Nematodes in Soil Health and Their Use as Indicators. J Nematol. 2001;33: 161–168. pmid:19265875

5. Forster S. Osmotic stress tolerance and osmoregulation of intertidal and subtidal nematodes. J Exp Mar Biol Ecol. 1997;54: 109–125.

6. Vermeeven H, Vanreusel A, Vanhove S. Species Distribution Within the Free-Living Marine Nematode Genus Dichromadora In the Weddell Sea and Adjacent Areas. Deep Res II. 2004;54: 1643–1664.

7. Chen C-A, Soo C-L, Balsamo M, Semprucci F. An approach based on nematode descriptors for the classification of ecological quality (EcoQ) of the Malaysian coasts. Mar Biodivers. 2018;48: 117–126.

8. Schratzberger M, Somerfield PJ. Effects of widespread human disturbances in the marine environment suggest a new agenda for meiofauna research is needed. Sci Total Environ. 2020;728: 138435. pmid:32570325

9. Singh R, Ingole BS, Nanajkar MR. The Life Cycle of The Free-Living Marine Nematode Innocuonema tentabunda De Man, 1890. Nematol Mediterr. 2009;37: 235–238.

10. Singh R, Ingole B. Life History of a Free-Living Marine Nematode Daptonema normandicum Reared in Laboratory. J Environ Biol. 2011;32: 147–152. pmid:21882647

11. Chen C, Long SM, Rosli NM. An Ecological Study of Free-living Marine Nematodes in Teluk Awar, Sarawak, Malaysia. 2012; 1–10.

12. Hong J-H, Semprucci F, Jeong R, Kim K, Lee S, Jeon D, et al. Meiobenthic nematodes in the assessment of the relative impact of human activities on coastal marine ecosystem. Environ Monit Assess. 2020;192: 81. pmid:31897611

13. Moreno M, Semprucci F, Vezzuli L, Balsamo M, Fabiano M, Albertelli G. The Use of Nematodes in Assessing Ecological Quality Status in The Mediterranean Coastal Ecosystems. Ecol Indic. 2011;11: 328–336.

14. Ürkmez D, Sezgin M, Bat L. Use of Nematode Maturity Index for The Determination of Ecological Quality Status: A Case Study from The Black Sea. Mediterr Environ. 2014;20: 96–107.

15. Semprucci F, Frontalini F, Sbrocca C. Meiobenthos and free-living nematodes as tools for biomonitoring environments affected by riverine impact. 2015.

16. Liu XS, Zhang ZN, Huang Y. Sublittoral Meiofauna With Particular Reference to Nematodes in The Southern Yellow Sea, China. Estuar Coast Shelf Sci. 2007;71: 616–628.

17. Cai L, Hong H, Zou C, Fang S, Lv X, Zhang Y. Studies on Species Composition of Marine Nematodes and Their Food Types in South Taiwan Strait. J Oceanogr Taiwan Strait. 2000;19: 212–217.

18. Liu XS, Cheung SG, Shin PKS. Meiofauna With Special Reference to Nematodes in Trawling Ground of Subtropical Hong Kong. Mar Pollut Bull. 2009;58: 607–615. pmid:19187945

19. Saving the Sea: Wang Ming-hsiang and His Squid Birthing Suites. Taiwan Panorama. 11 Nov 2019: 1. Available: https://nspp.mofa.gov.tw/nsppe/news.php?post=165423&unit=412&unitname=RELATED-READS&postname=Saving-the-Sea:-Wang-Ming-hsiang-and-His-Squid-Birthing-Suites.

20. Lin YZ. Taiwan Reef Check 2018 Reveals Reefs Covered by Silt. Green News Taiwan. 12 Sep 2018: 1. Available: https://e-info.org.tw/node/213943.

21. Fichet D, Bowher G, Radenac G, Miramand P. Concentration and Mobilisation of Cd, Cu, Pb And Zn By Meiofauna Populations Living in Harbour Sediment: Their Role in The Heavy Metal Flux from Sediment to Food Web. Sci Total Environ. 1999;243: 263–272. pmid:10635598

22. Moreno M, Vezzulli L, Marin V, Laconi P, Albertelli G, Fabiano M. The Use of Meiofauna Diversity as An Indicator of Pollution in Harbours. ICES J Mar Sci. 2008;68: 1428–1435.

23. Losi V, Ferrero TJ, Moreno M, Gaozza L, Rovere A, Firpo M, et al. The Use of Nematodes in Assessing Ecological Conditions in Shallow Waters Surrounding a Mediterranean Harbour Facility. Estuar Coast Shelf Sci. 2013;130: 209–221.

24. Cai L, Fu S, Zhou X, Tseng LC, Hwang JS. Benthic meiofauna with emphasis on nematode assemblage response to environmental variation in the intertidal zone of the Danshuei River estuary, northwest Taiwan. Ecol Res. 2020;35: 857–870.

25. Somerfield PJ, Warwick RM. Meiofauna Techniques. Methods Study Mar Benthos. 2013; 253–284.

26. Bale AJ, Kenny AJ. Sediment Analysis and Seabed Characterisation. Methods for the Study of Marine Benthos. Oxford, UK: Blackwell Science Ltd; 2005. pp. 43–86. https://doi.org/10.1002/9780470995129.ch2

27. Wentworth CK. A Scale of Grade and Class Terms for Clastic Sediments. J Geol. 1922;30: 377–392.

28. Semprucci F, Frontalini F, Covazzi-Harriague A, Coccioni R, Balsamo M. Meio- And Macrofauna In the Marine Area of The Monte St. Bartolo Natural Park (Central Adriatic Sea, Italy). Sci Mar. 2013;77: 189–199.

29. Burgess R. An improved protocol for separating meiofauna from sediments using colloidal silica sols. Mar Ecol Prog Ser. 2001;214: 161–165.

30. Seinhorst JW. A Rapid Method for The Transfer of Nematodes from Fixative to Anhydrous Glycerin. Nematologica. 1959;4: 67–69.

31. Platt HM, Warwick RM. Free-Living Marine Nematodes, Part I. British Enoplids. 1983.

32. Warwick RM, Platt HM, Sommerfield PJ. Free-Living Nematodes (Part III) Monhysterids. 1998.

33. Shabdin ML, Rosli NM, Chen CA. Free-Living Nematodes in Sarawak Coastal Waters. Penerbit UMT, Kuala Terengganu. 2013.

34. Chen CA, Mohd Long S. A new species of Setosabatieria Platt, 1985 (Nematoda: Comesomatidae) from Chek Jawa, Singapore, with a key to valid species of the genus. Mar Biol Res. 2015;11: 203–208.

35. Chen C-A, Mohd Long S. A New Marine Nematode Species, Metalinhomoeus ramsarensis (Linhomoeidae Filipjev, 1922) from Kuching Wetland National Park, Sarawak, Malaysia. Borneo J Resour Sci Technol. 2015;5: 11–15.

36. Chen CA, Nguyen DT, Smol N. Two new free-living marine nematode species from an intertidal sandyrocky shore on Pulau Ubin, Singapore with a key to the valid species of the genera Prooncholaimus and Acanthonchus. Raffles Bull Zool. 2015;2015: 68–74.

37. Wieser W. Die Beziehung Zwischen Mundh Ö Hlengestalt, Ern Ä Hrungsweise Und Vorkommen Bei Freilebenden Marinen Nematoden. Ark für Zool. 1953;4: 439–484.

38. Romeyn K, Bouwman L. Food Selection and Consumption by Estuarine Nematodes. Hydrobiol Bull. 1983;17: 103–109.

39. Jensen P. Feeding Ecology of Free-living Aquatic Nematodes. Mar Ecol Ser. 1987;35: 187–196.

40. Moens T, Vincx M. Observations on the Feeding Ecology of Estuarine Nematodes. J Mar Biol Assoc United Kingdom. 1997;77: 211–227.

41. Bongers T, Alkemade R, Yeates GW. Interpretation of disturbance-induced maturity decrease in marine nematode assemblages by means of the Maturity Index. Mar Ecol Prog Ser. 1991;76: 135–142.

42. Bongers T. The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia. 1990;83,: 14–19. pmid:28313236

43. Hua E, Zhang ZN, Zhang Y. Environmental Factors Affecting Nematode Community Structure in The Changjiang Estuary and Its Adjacent Waters. J Mar Biol Assoc United Kingdom. 2009;89: 109–117.

44. Leduc D, Pilditch C. Effect of A Physical Disturbance Event on Deep-Sea Nematode Community Structure and Ecosystem Function. J Exp Mar Bio Ecol. 2013;440: 35–41.

45. Liu XS, Xu M, Z J.H., Liu D, Li X. Community Structure and Biodiversity of Free-Living Marine Nematodes in The Northern South China Sea. Acta Oceanol Sin. 2015;34: 77–85.

46. Xu M, Liu QH, Zhang ZN, Liu XS. Response of Free-Living Marine Nematodes to The Southern Yellow Sea Cold Water Mass. Mar Pollut Bull. 2016;105: 58–64. pmid:26965091

47. Gyedu-Ababio TK, Furstenberg JP, Baird D, Vanreusel A. Nematodes as Indicators of Pollution: A Case Study from The Swartkops River System, South Africa. Hydrobiologia. 1999;397: 155–169.

48. Steyaert M, Vanaverbeke J, Vanreusel A, Barranguet C, Lucas C, Vincx M. The Importance of Fine-Scale, Vertical Profiles in Characterising Nematode Community Structure. Estuar Coast Shelf Sci. 2003;58: 353–366.

49. Ng W-L, Chen C-A, Kawi SD ak., Musta B, Chan T-Y. Effects of hydrogen peroxide treatment on the particle size distribution of hydrothermal vent sediments: A case study in Guishan Island, Taiwan. Borneo J Mar Sci Aquac. 2019;3: 52–56.

50. Urban-Malinga B. Meiobenthos In Marine Coastal Sediments. Sediment Coast Zo from High to Low Latitudes Similarities Differ. 2014;388: 59–78.

51. Clarke K, Warwick R. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation. Nat Environ Res Counc. 1994;144.

52. Lewis M. Sediment habitat assessment for targeted near-coastal areas. In Bortone S.A. (ed.),. Estuar Indic CRC Mar Sci Ser. 2005; 79–98.

53. Chen C, Long SM, Rosli NM. Spatial distribution of tropical estuarine nematode communities in Sarawak, Malaysia (Borneo). Raffles Bull Zool. 2012;60: 173–181.

54. Gawad E, Azab Ma, Lotfy M. Assessment of Organic Pollutants in Coastal Sediments, UAE. Env Geol. 2008;54: 1091–1102.

55. Pinto TK, Benvenuti C. Changes in Nematoda Community Structure Related to Sewage Discharge. J Coast Res. 2006;1: 1807–1810.

56. Sahraean N, Campinas T, Keivan B, Khanaghah E, Mosallanejad H, Ranst E Van, et al. Effects of pollution on nematode assemblage structure and diversity on beaches of the northern Persian Gulf. Hydrobiologia. 2017;799: 349–369.

57. Ullberg J, Ólafsson E. Free-living marine nematodes actively choose habitat when descending from the water column. Mar Ecol Prog Ser. 2003;260: 141–149.

58. Venekey V, Melo TPG, Rosa Filho JS. Effects of seasonal fluctuation of amazon river discharge on the spatial and temporal changes of meiofauna and nematodes in the amazonian coast. Estuar Coast Shelf Sci. 2019;227: 106330.

59. Traunspurger W, Threis I, Majdi N. Vertical and Temporal Distribution of Free-Living Nematodes Dwelling in Two Sandy-Bed Streams Fed by Helocrene Springs. Nematology. 2015;17: 923–940.

60. Neira C, King I, Mendoza G, Sellanes J, de Ley P, Levin LA. Nematode Community Structure Along A Central Chile Margin Transect Influenced by The Oxygen Minimum Zone. Deep Sea Res Part I Oceanogr Res Pap. 2013;78: 1–15.

61. Kudale MD. Impact of Port Development on The Coastline and The Need for Protection. Indian J Geo-Marine Sci. 2010;39: 597–604.

62. Zheng GJ, Richardsan B. Petroleum Hydrocarbons and Polycyclic Aromatic Hydrocarbons (Palls) in Hong Kong Marine Sediments. Chemosphere. 1999;38: 2625–2632.

63. Mahmoudi E, Essid N, Beyrem H, Hedfi A, Boufahja F, Vitiello P, et al. Effects of Hydrocarbon Contamination on A Free-Living Marine Nematode Community: Results from Microcosm Experiments. Mar Pollut Bull. 2005; 1197–1204. pmid:15896812

64. Yang GP. Polycyclic Aromatic Hydrocarbons in The Sediments of The South China Sea. Environ Pollut. 2000;108: 163–171. pmid:15092946

65. Gingold R, Moens T, Rocha-Olivares A. Assessing the Response of Nematode Communities to Climate Change-Driven Warming: A Microcosm Experiment. PLoS ONE,. 2013;8: e66653. pmid:23825552

66. Fu S, Rao Y, Chen X, Zhou X, Wu C, Li X, et al. Comparison of benthic nematode assemblages in native mangrove forest and exotic mangrove plantations (Sonneratia apetala Buch-Ham) along the South China Coast. Mar Pollut Bull. 2021;166: 112249. pmid:33756351

67. Semprucci F, Balsamo M. Free-living marine nematodes as bioindicators: past, present and future perspective. Environ Res J. 2014;6: 17–35.

68. Semprucci F, Balsamo M, Frontalini F. La comunidad de nematodos de una laguna costera (Laguna de Varano, Italia meridional): Patrones de la ecología y la biodiversidad. Sci Mar. 2014;78: 579–588.

69. Armenteros M, Ruiz-Abierno A, Fernández-Garcés R, Pérez-García JA, Díaz-Asencio L, Vincx M, et al. Biodiversity patterns of free-living marine nematodes in a tropical bay: Cienfuegos, Caribbean Sea. Estuar Coast Shelf Sci. 2009;85: 179–189.

70. Semprucci F, Boi P, Manti A, Harriague AC, Rocchi M, Colantoni P, et al. Benthic communities along a littoral of the Central Adriatic Sea (Italy). Helgol Mar Res. 2010;64: 101–115.

71. Moreno M, Ferrero TJ, Gallizia I, Vezzulli L, Albertelli G, Fabiano M. An assessment of the spatial heterogeneity of environmental disturbance within an enclosed harbour through the analysis of meiofauna and nematode assemblages. Estuar Coast Shelf Sci. 2008;77: 565–576.

72. Moreno M, Albertelli G, Fabiano M. Nematode response to metal, PAHs and organic enrichment in tourist marinas of the mediterranean sea. Mar Pollut Bull. 2009;58: 1192–1201. pmid:19386320

73. Pavlyuk ON, Trebukhova JA. Meiobenthos in Nha Trang Bay of the South China Sea (Vietnam). Ocean Sci J. 2006;41: 139–148.

74. Hawkins JP, Roberts CM. Effects of Recreational SCUBA Diving on Fore-Reef Slope Communities of Coral Reefs. Biol Conserv. 1992;62: 171–178.

75. Hawkins JP, Roberts CM, Van’T Hof T, De Meyer K, Tratalos J, Aldam C. Effects of Recreational Scuba Diving on Caribbean Coral and Fish Communities. Conserv Biol. 1999;13: 888–897.

76. Hawkins JP, Roberts CM, Kooistra D, Buchan K, White S. Sustainability of Scuba Diving Tourism on Coral Reefs of Saba. Ocean Coast Manag. 2005;33: 373–387.

77. Aeby G, Williams G, Franklin E, Haapkyla J, Drew Harvell C, Neale S, et al. Growth Anomalies on The Coral Genera Acropora And Porites Are Strongly Associated with Host Density and Human Population Size Across the Indo-Pacific. PLoS One. 2011;6: e16887. pmid:21365011

78. Altizer S, Ostfeld R, Johnson P, Kutz S, Harvell C. Climate Change and Infectious Diseases: From Evidence to A Predictive Framework. Science (80-). 2013;341: 514–519. pmid:23908230

79. Toyoshima J, Nadaoka K. Importance of Environmental Briefing and Buoyancy Control on Reducing Negative Impacts of SCUBA Diving on Coral Reefs. Ocean Coast Manag. 2015;116: 20–26.

80. Lamb JB, True JD, Piromvaragorn S, Willis BL. Scuba diving damage and intensity of tourist activities increases coral disease prevalence. Biol Conserv. 2014;178: 88–96.

81. Lyons PJ, Arboleda E, Benkwitt CE, Davis B, Gleason M, Howe C, et al. The Effect of Recreational SCUBA Divers on The Structural Complexity and Benthic Assemblage of a Caribbean Coral Reef. Biodivers Conserv. 2015;24: 3491–3504.

82. Giglio V, Luiz O, Ferreira C. Ecological Impacts and Management Strategies for Recreational Diving: A Review. J Environ Manage. 2020;256: 109949. pmid:31818747

83. Xuan QN, Vanreusel A, Thanh NV, Smol N. Biodiversity of meiofauna in the intertidal Khe Nhan Mudflat, can gio mangrove forest, Vietnam with special emphasis on free living nematodes. Ocean Sci J. 2007;42: 135–152.

84. Ansari K, Pattnaik A, Rastogi G, Bhadury P. Inventory of Free-Living Marine Nematodes from Asia’s Largest Coastal Lagoon, Chilika, India. Wetl Ecol Manag. 2015;23: 881–890.

85. Soo C-L, Chen C-A, Mohd-Long S. Assessment of Near-Bottom Water Quality of Southwestern Coast of Sarawak, Borneo, Malaysia: A Multivariate Statistical Approach. J Chem. 2017;2017: 1–12.

86. Fischer K, Kölzow N, Höltje H, Karl I. Assay Conditions in Laboratory Experiments: Is the Use of Constant Rather Than Fluctuating Temperatures Justified When Investigating Temperature-Induced Plasticity? Oecologia. 2011;166: 23–33. pmid:21286923

87. Nanajkar MR, Ingole B. Nematode Species Diversity as Indicator of Stressed Benthic Environment Along the Central West Coast of India. Conf Artic Divers Life Process from Ocean Land, Eds Desai, PV; Roy, R. 2007; 42–52.

88. Somerfield PJ, Gee JM, Warwick RM. No TitleSoft Sediment Meiofaunal Community Structure in Relation to A Long-Term Heavy Metal Gradient in The Fal Estuary System. Mar Ecol Prog Ser. 1994;105: 79–88.

89. Semprucci F, Colantoni P, Balsamo M. Is maturity index an efficient tool to assess the effects of the physical disturbance on the marine nematode assemblages?—A critical interpretation of disturbance-induced maturity successions in some study cases in Maldives. Acta Oceonologica Sin. 2016;35: 89–98.

90. Sahraeian N, Hosseinzadeh H, Mosallanejad H, Ingels J. Temporal and spatial variability of free-living nematodes in a beach system characterized by domestic and industrial impacts (Bandar Abbas, Persian. Ecol Indic. 2020;118: 106697.

91. Mohebbi Nozar S, Ismail WR, Zakaria MP. Distribution, Sources Identification, and Ecological Risk of PAHs and PCBs in Coastal Surface Sediments from the Northern Persian Gulf. Hum Ecol Risk Assess An Int J. 2014;20: 1507–1520.

92. Soko MI, Gyedu-Ababio T. Free-Living Nematodes as Pollution Indicator in Incomati River Estuary, Mozambique. Open J Ecol. 2019;9: 117–133.

93. Balsamo M, Semprucci F, Frontalini F, Coccioni R. Meiofauna As A Tool for Marine Ecosystem Biomonitoring. Cruzado A. (Ed) Marine Ecosystems. InTech Publ. 2012;4: 77–104.

94. Gyedu-ababio TK, Baird D. Response of meiofauna and nematode communities to increased levels of contaminants in a laboratory microcosm experiment. 2006;63: 443–450.