Diversity and Distribution of Polychaetes in Mangroves of East Coast of India Diversity and Distribution of Polychaetes in Mangroves of East Coast of India

This research article reports an exhaustive account on the mangrove-associated poly- chaetes. Polychaetes are an important component in marine benthic communities and they play a major ecological role in mangrove ecosystem. This article gives an overview of polychaete diversity associated to five major mangrove forests of east coast of India (Muthupettai, Pichavaram, Coringa, Bhitarkanika and Sundarban). The results of this sur - vey indicated that the physicochemical parameters did not vary much except a few parameters that showed only marginal variations. With regard to the macrobenthic organisms, the polychaetes topped the list. Crustaceans were found to be the next dominant group in the order of abundance and followed by gastropods and bivalves of the total benthic organ- isms collected. The results of the statistical analysis revealed that the parameters such as salinity, pH, silt, clay, total organic carbon (TOC), total nitrogen (TN) and total phosphate (TP) were manifested as best match in determining benthic fauna distributions followed by TOC, slit, clay and TP. The maximum number of polychaete species was recorded from Sundarban mangroves (68 species) and minimum in Muthupettai mangroves (39 species).


Introduction
Mangroves are unique coastal ecosystem contributing as a rich store house of biodiversity. Mangrove forests are extremely important coastal resources [1] which play a pivotal role in socio-economic development. It also plays a major role as nursery ground for juveniles of a plethora of fin and shell fishes. A total of 54 mangrove species belonging to 20 genera and detection of pollution and are considered as the taxonomic group with the highest level of sensitivity to perturbation of the soft substrata [17].
No comprehensive study has been undertaken so far on benthic biodiversity in general and polychaete taxonomy in particular in the mangroves of east coast of India. Taking cognizance of the facts stated above, a case study on the diversity and distribution pattern of polychaetes in five major mangroves of east coast of India is posted in this article.

Study area
For the present investigation, survey was conducted in five different mangrove ecosystems of east coast of India. The description of the study area is detailed in the following section (Table 1 & Figure 1).
The water, sediment and macrofaunal samples were collected seasonally from five major mangrove ecosystem of east coast of India during 2013-2014. In each mangrove, three stations representing i) Land ward zone, ii) core mangrove, and iii) Seaward zone, were fixed and thus altogether 15 stations were sampled:

Collection of water and sediment samples
The environmental parameters such as pH, salinity, temperature and dissolved oxygen (DO was measured following the modified Winkler's method [18] in the site itself. The sediment nutrient parameters such as total nitrogen (TN) was estimated by following the method of Strickland and Parsons [18], total phosphorous (TP) by following the method of Menzel and Corwin [19]; and total organic carbon (TOC) by following the standard method of El Wakeel and Riley [20].

Biological sample (field and lab routines)
In each station, three replicate samples were collected using Peterson grab. This type of grab is considered to be the most efficient gear in obtaining the good penetrative samples in shallow water environments. The grab employed was found to take a sample covering an area of 0.1m 2 .
The procedure adopted for sampling was following the method of Mackie [21]. After collecting the samples, they were emptied into a plastic tray. The larger organisms were handpicked immediately from the sediments and then sieved through 0.5 mm mesh screen. The organisms retained by the sieve were placed in a labeled container and fixed in 5-7% formalin. Subsequently, the organisms were stained with Rose Bengal solution (0.1 g in 100 ml of distilled water) for greater visibility during sorting. All the species were sorted, enumerated and identified to the advanced possible level with the consultation of available literature. The works of Fauvel [22] and Day [23] and http://www.marinespecies.org/polychaeta/ were referred for identification.

Statistical analyses
The data were approached to various statistical methods namely univariate, graphical/distributional and multivariate methods available in PRIMER (Ver. 7.) statistical software [24]. The data were analyzed for diversity index (H′) using the method of Shannon -Wiener's formula [25]; for species richness (d) using the formula of Margalef [26] and species evenness (J') using Pielou [27].
Cluster analysis was done to find out the similarities between the samples/stations/regions. The most commonly used clustering technique is the hierarchical agglomerative method. MDS (non -metric Multi-Dimensional Scaling) [28,29], was used to find out the similarities (or dissimilarities) between each pair of entities to produce a 'map', which would ideally show the interrelationships of all.
The principal component analysis-Bi-plot (PCA-Bi-plot), a multivariate procedure capable of providing a data reduction and easy visualization through the Pearson correlation between the physicochemical parameters and sampling stations were performed using XLSTAT-Pro version 5.1.4. Canonical Correspondence Analysis (CCA) was also done to relate the abundance of benthic species with linear combination of environmental variables [30,31].
Canonical Correspondence Analysis (CCA) allows to obtaining a simultaneous representation of the sites, the objects, and the variables in two or three dimensions that is optimal for a variance criterion [30]. To confirm the results obtained through CCA, BIO-ENV procedure [32] was also employed. A weighted Spearman rank correlation coefficient (ρω) was used to determine the harmonic rank correlation between the biological variable and all possible combinations of the environmental variables.

Environmental variables
The mean values of physicochemical parameters recorded at each sampling station are summarized in

Principal component analysis
The PCA was performed using physicochemical parameters to set a well defined distinction between the stations and the parameters. The PCA drawn for five mangroves showed 85.67% variance of the total axis wherein the first axis (F1) explained up to 62.47% of the total variance and F2 axis explained only 23.20% of the total variance. When the results were viewed, the parameters such as salinity, pH, Silt, Clay, TN, TP and TOC got positively correlated with MUT-1, PIC-1, BIT-2, PIC-2, SUN-2 and SUN-3 and MUT-1 while water temperature, DO and sand were negatively correlated with stations MUT-3, PIC-3, BIT-1, BIT-3, SUN-1, COR-1, COR-2 and COR-3 (Figure 2).

Species composition of macrofauna
In the present study, organisms of the following five groups were recorded in the benthic samples collected: 1. polychaetes, 2. crustaceans, 3. bivalves, 4. gastropods and 5. 'others.' As many as 97 species of macrofauna were recorded from 5 mangrove ecosystems of the present  study. Of these species, polychaetes were found to be the largest component in the collection with 68 species. Crustaceans emerged as next dominant group in the order of abundance with 11 species. The bivalves and gastropods came next in the order with 8 and 6 species respectively and the group 'others' came last in the order with 4 species.
In Muthupettai mangroves, a total of 69 species were recorded. Among these, 39 species belonged to polychaetes, 10 species to crustaceans, 8 species each to bivalves and gastropods and 4 species to group 'others.' With respect to Pichavaram mangroves, a total of 88 species of macrofauna were recorded. Among these, there were 59 species of polychaetes, 10 species were crustaceans, 8 and 7 species were bivalves and gastropods respectively and 4 species of 'others.' Regarding Coringa, 77 species of macrofauna were found. Among these, 50 species of polychaetes, 9 species of crustaceans and 8 and 7 species of bivalves and gastropods and 3 species of 'others' were recorded. Coming to Bhitarkanika mangroves, 81 species of macrofauna were found. Among these, 54 species of polychaetes, 10 species of crustaceans and 7 species each of bivalves and gastropods and 3 species of 'others' were recorded.
Coming to Sundarban mangroves, 97 species of macrofauna were found. Of these, 68 species of polychaetes, 11 species of crustaceans and 8 and 6 species of bivalves and gastropods respectively, and 4 species of 'Others' were recorded.     in their occurrence in the entire mangrove ecosystem. With respect to bivalves, Anadara rhombea, Crassostrea madrasensis, Katelysia opima, Meretrix meretrix, Meretrix casta, Perna indica, and similarly among gastropods, Cerithidea cingulata, Nassarius stollatus, Turritella acutangula and Murex trapa were recorded frequently. Group "others" constitute fish larvae, sea urchins, crab and foraminiferans. The common macro benthic species recorded in various stations of five mangrove ecosystems is shown in Table 3 & Figure 3.

Population density of macrofauna
The results of population density recorded in five mangroves are given in the following section: In Muthupettai mangroves, the population density of benthic macrofauna varied from 417 to 3545nos/m 2 with the maximum was noticed during summer and minimum during monsoon.
Coming to Pichavaram mangroves, the density of benthic organisms varied between 451 and 5645 nos/m 2 with during summer and minimum during monsoon. Regarding Coringa mangroves, the density of benthic organisms ranged from 386 to 4262 nos/m 2 with maximum during summer and minimum during monsoon. Coming to Bhitarkanika mangroves, the density of benthic organisms varied between 433 and 4862 nos/m 2 with maximum during summer and minimum during monsoon. With respect to Sundarban mangroves, the density of organisms varied from 511 to 6845 nos/m 2 . The minimum density was recorded monsoon and maximum during summer. Among the mangroves, the maximum density of macrofauna was recorded in Sundarbans (6845 nos/m 2 ) during summer and minimum in Muthupettai (3545nos/m 2 ) during monsoon (Figure 4).

Percentage composition of benthos
The percentage composition of macrofauna recorded in five different mangroves ecosystems are given below: In Muthupettai, when the results of percentage composition of benthic fauna were viewed, polychaetes constituted the maximum with 54% of the total benthic organisms followed by crustaceans with 14%, bivalves with 12%, gastropods with 13% each and group 'others' with 7% to the samples collected in Muthupettai mangroves. With respect to Pichavaram mangroves, polychaetes continued to emerge as the dominant group in terms of abundance with a percentage occurrence of 56%. Crustaceans ranked second with a percentage contribution of 15%. Gastropods, bivalves contributed 11%, 12% respectively and 'others' with 6% to the total benthic organisms recorded.
Regarding Coringa, as in other mangroves, polychaetes continued to be the dominant group with 61%, followed by crustaceans, bivalves, gastropods and 'others' with 13%, 12%, 9% and 5% respectively. Coming to Bhitarkanika mangroves, polychaetes remained as the dominant group with a percentage contribution of 53%. Crustaceans were found to be the next dominant group with a percentage contribution of 13%. Gastropods, bivalves and 'others' contributed 8%, 11% and 5% respectively to the total benthic organisms collected. In Sundarban mangroves, polychaetes topped the list in terms of abundance with a percentage of 62%. Crustaceans formed second dominant group with a percentage contribution of 15%. Gastropods, bivalves contributed with 7% and 10% respectively and 'others' with 6% of the total benthic organisms ( Figure 5).

Diversity indices
The Diversity indices (mean value) recorded at each sampling station is summarized in Table 4. The species diversity varied from 3.018 to 4.476 with maximum in Sundarbans and minimum in Muthupettai mangroves; species richness fluctuated from 3.216 to 4.194 with maximum in Sundarbans and minimum in Coringa mangroves; with respect to Pielou's evenness, it varied from 0.852 to 0.991 with maximum in Bhitarkanika and minimum in Coringa mangroves.

Cluster analysis
The seaward stations (MUT-1, PIC-1, COR-1, BIT-1 and SUN-1) in all the mangroves got grouped at the highest level of similarity followed by stations of core mangrove zone (MUT-2, PIC-2, COR-2, BIT-2 & SUN-2) and stations of landward zone (MUT-3, PIC-3, COR-3, BIT-3 & SUN-3) got grouped to form cluster based on the species composition with the exception of a few outliers (stations), which might be due to the species commonality between zones. This fact was further confirmed through MDS, and the results also revealed the same pattern of groupings as recognized in cluster analysis (Figure 6).

Canonical correspondence analysis (CCA)
Canonical correspondence analysis (CCA) was done to ascertain the relationship between the physicochemical parameters and benthic faunal density. The CCA drawn for five mangrove ecosystem showed 91.43% variance of the total axis wherein the F1 axis showed 74.56% and F2 axis 16.87% of the total variance. The environmental parameters such as salinity, Silt, Clay, TOC, TP and TN were showing strong correlation with the benthic faunal diversity, while other parameters like water temperature, depth, sand and DO had weak correlation with the benthic faunal distribution (Figure 7).

BIO-ENV (biota-environment matching)
In the BIO-ENV procedure, which was employed to measure the agreement between the rank correlations of the biological (Bray-Curtis similarity) and environmental (Euclidean distance) matrices, ten environmental variables were allowed to match the biota. The results of best  combinations are given in Table 5. In this case, as evidenced in CCA plot, salinity, silt, clay, TOC, total nitrogen and total phosphorous were featured as the major variables explaining the best match (0.90) with faunal distributions followed by pH, TOC and total nitrogen were also got manifested in the second best variable combinations in determining the faunal distribution in the mangrove ecosystems.

Discussion
Composition of benthic communities and their role varies from one habitat to another depending upon the water and sediment characteristics of the mangroves. The distribution of mangrove fauna in relation to water quality has been described quantitatively [33]. Among the five mangroves, the maximum temperature was recorded at Muthupettai during summer and minimum in Sundarbans, which could be ascribed to the effect of atmospheric cooling. Similar conclusion was also drawn earlier by Bolam et al. [34] in UK continental shelf waters and in shelf waters of southeast coast of India [35]. The temperature levels recorded presently are comparable with the study made by Kathiresan [36] who reported the temperature range of 28-31°C.
The high salinity values observed during summer compared to other seasons is might be due to low rain fall and the rise in atmospheric temperature resulting in high evaporation rate of the surface water. Similar seasonal variations were observed by Manokaran [35] in the inshore waters of Parangipettai and Cuddalore; by Murugesan et al. [37] in Tuticorin coastal waters and Rahaman et al. [38] in Sundarbans mangroves; Sivaraj et al. [39] in Vellar-Coleroon estuarine system.
In the present study, the maximum pH of 8.23 was recorded during summer and minimum of 7.1 was recorded during wet season. Hydrogen-ion concentration was found to vary among the five mangroves and was alkaline throughout the study period. Higher pH observed in summer season could be attributed to the removal of CO 2 by the photosynthetic organisms and the lower pH during monsoon season could be due to the dilution of saline water with fresh-water inflow from nearby sources as has been reported by Murugesan et al. [37].
Coming to dissolved oxygen, (DO) it varied from 3.80 to 7.27 mg/l with the maximum (7.27) during wet season and minimum 3.80 was recorded during dry season. All the stations of various mangroves showed the similar seasonal pattern in the distribution of dissolved oxygen with minimum value during dry months and maximum during wetter months. The relatively low DO values observed in the summer are attributed to the entry of high saline waters in to the mangroves, as well as fluctuations in temperature and salinity, which in turn affect the dissolution of oxygen [40]. This fact is in close agreement with earlier studies done elsewhere [38,41].
Mangrove ecosystems are able to store large amounts of organic carbon [42]. In the present study, the maximum TOC of 16.52mgC/g was recorded at SUN-12 during dry season and minimum of 6.45mgC/g was recorded at COR-13 during wet season. As noticed in temperature and salinity, all the stations showed similar seasonal pattern in the distribution of organic carbon content with maximum value during dry months and minimum during wet months. Similarly, Hasrizal et al. [43] studied the seasonal changes of organic carbon content in the surface sediments of the Terengganu near shore coastal area of Malaysia with maximum value during postmonsoon and summer seasons and they also opined that the sediment characteristics and the organic carbon concentration are largely influenced by southwest and northeast monsoons.
In the present study, total nitrogen content showed striking seasonal variation with maximum TN (5.98 μg/g) was recorded during monsoon and minimum (3.48 μg/g) during dry season. Likewise, the maximum TP (1.73 μg/g) was recorded during wet season and minimum (0.88 μg/g) was recorded during dry season. The maximum values in wet season might be attributed to the higher amount of rainfall and river runoff as has been reported earlier by Sreedevi [44]. Similarly Kamykowski and Zentoura [45] also opined that the accumulation of nitrite in the near bottom samples depends on diffusion from sediments as well as mechanisms such as nitrification near the sediment and water interface. Similar observation was made by Gouda and Panigrahy [46] in Rushikulya estuary, Orissa, east coast of India. Manikoth and Salih [47] recorded high nitrogen concentration during monsoon season in the Vembanad estuarine complex, southwest coast of India. Joshi and Ghose [48] studied nutrient characteristics of Sundarban mangroves. Martin et al. [49] studied on the benthic fauna in a tropical estuary of Cochin backwaters and Sekar et al. [50] in Pichavaram and Muthupettai mangroves in relation to nutrient characteristics.
Studies on the sediment composition are of paramount importance in benthic ecology. The comprehensive knowledge on the sediment composition is a pre-requisite and inevitable one to understand the benthic ecology [51]. The nature of the substratum has a profound effect on the bottom fauna and conversely, the benthos can influence the sediment characteristics. Gray and Snelgrove and Butman [52,53] posted the information regarding the relationship between sediments and benthic organisms. They also pointed out that the grain size distribution of the sediments is of great importance in determining the distribution of benthos. Snelgrove and Butman [53] also concluded that the relationship was a complex interaction of the seabed flow and sediment characteristics and that could explain the distribution of organisms across all sedimentary habitats.
The correlation between the physicochemical parameters and benthic faunal density for the surveyed five mangrove ecosystem showed that the environmental parameters such as salinity, Silt, Clay, TOC, TP and TN were showing strong correlation with the benthic faunal diversity, while other parameters like water temperature, depth, sand and DO had weak correlation with the benthic faunal distribution. Similar variables combination were reported earlier by Sundaray et al. in Mahanadi River [54]; Satheeshkumar et al. [55] in Pondicherry coast; Sivaraj et al. [56] in Nandgoan coastal waters; Sivaraj et al. [41] in Vellar-Coleroon estuarine system.
Percentage contribution of benthic species composition of the present study showed in the order of polychaetes, crustaceans, bivalves, gastropods and groups 'others'. The dominance of polychaetes in terms of density and species composition in diverse ecological niche is due to their high degree of adaptability to a wide range of environmental factors. Similar preponderance of polychaetes has been observed earlier by Kumar [32] in Cochin backwaters; Prabha Devi [57] in Coleroon estuary, and Ansari et al. [58], in Mandovi estuary. Athalye and Gokhale [59] reported the dominance of polychaetes followed by gastropods, bivalves, and hermit crabs in Thane creek, Mumbai. The dominance of polychaetes might be due to the fact that firm substrate provided by roots and dense canopy of the mangroves which also provide protection against desiccation [60]. Similar dominance of polychaetes was also reported in other tropical waters [61,62].
In a study conducted by Harkantra and Parulekar [63], polychaetes outnumbered the other faunal groups where the substratum was mainly composed of mud. Bhat and Neelakandan [64] also observed maximum number of polychaetes in the clayey-silty substratum, the fine particles of mud and clay substratum, which retains more water than coarse particles (sand and gravel). Such fine deposits or particles are commonly composed of decomposable organic constituents. As the organic content represents an important direct or indirect food source for benthic organisms, elevated organic matter may result in an enhancement of benthic faunal diversity [52,65]. Therefore, it is clear that polychaetes abound in finer sediments as noticed by the above referred researchers. This fact also corroborates the results of present study. The population density of macrofauna is governed by various environmental variables such as temperature, salinity, sediment type, organic carbon level in the sediments besides tidal action [66]. Monsoon months registered low density followed by gradual increase in postmonsoon and peaked during summer season, which are in agreement with the results of Sekar et al. [38,50].
The population density recorded presently is comparable with the following studies made in the back waters along the east and west coasts of India: Harkantra et al. [ In the present study, a marked seasonal variation in the Shannon diversity was found with minimum diversity value (3.018) in Muthupettai mangroves during monsoon and maximum (4.476) in Sundarbans mangroves during dry season. Similar range of diversity values was recorded earlier in Vellar estuary [71]. Shillabeer and Tapp [72] stated that the estuarine and mangrove environment is far more dynamic than the fully marine and therefore, there may be a wide range of variations in the benthic diversity of an estuary.
As in the species diversity, species richness values were also low during wet season and high during dry season, which might be due to adaptability to high salinities at high temperatures than at low temperatures [73], as a result more marine forms are able to flourish in tropical waters [74].
The trend with respect to richness values of the present study is evident in the studies made by Raveenthiranath Nehru [14] in Coleroon estuary and Sebastin Raja [14] in Sunnambar estuary; Palanisamy and Anisa [51] in Pondicherry coastal waters. With respect to evenness (J'), it largely followed the trend of species diversity.
With respect to classification and ordination techniques, the stations of marine zone (seaward) grouped at the highest level of similarity followed by stations of core mangrove zone and stations of fresh water zone (landward zone) grouped to form clusters based on the species composition. The physicochemical parameters such as salinity, Silt, Clay, TOC, TP and TN in landward zone and core mangrove were found relatively similar and it highly influenced the benthic faunal diversity, while in seaward zone the trends of the same parameters varied significantly and it didn't affect the distribution and diversity of the benthic fauna. The MDS results also largely followed the trend of dendrogram. Investigation similar to this was carried out by Sivaraj et al. [41] who made a comparative study of Vellar-Coleroon estuarine system using macrobenthic communities through cluster analysis. The stress value observed in MDS plot is comparable with the studies [75][76][77].
Canonical Correspondence Analysis (CCA) was done to ascertain the relationship between the physicochemical parameters and benthic faunal density. Similar combinations of environmental variables influencing benthic faunal distribution was reported in Nandgaon coastal waters, Maharashtra, India [56]; Sivaraj et al. [41] in Vellar-Coleroon estuarine system. This fact was further confirmed through BIO-ENV, which yielded the combinations of six environmental entities (salinity-silt-clay-TOC-TN-TP) as best match 'defining' the faunal distributions. The associated coefficient of environmental to biotic similarity was 0.90. True to this, studies [39,71] reported the similar combinations of environmental variables influencing the benthic faunal distribution. Clarke and Ainsworth [62] also reported the organic carbon-sediment particle size, to constitute the best match explaining the distribution of meiobenthic organisms. Similarly, Mackie et al. [78,79] reported the combination as silt-clay-organic carbon forming the best match in explaining the faunal distribution. The combinations recognized in the above referred studies corroborate the results of the present study.
Comparing our own data with the studies made elsewhere in mangroves of other Asian countries, a few inferences could be drawn. In our study, as many as 68 species of polychaetes were recorded from 5 mangrove ecosystems of the present study. The density and number of species recorded presently is comparable with the works carried out in mangroves of other Asian countries barring a few variations in their density and diversity which might be due to the dynamic nature of the mangrove environment. Shillabeer and Tapp [72] stated that the mangrove environment is far more dynamic than the fully marine and therefore, there is every possibility in the variations in the occurrence of species. Similarly, there was no pronounced variation with respect to commonality in the species occurrence between our data and data of others. With regard to representation of polychaete families, by and large the representatives from Errant polychaetes were found to outnumber compared to sedentary counterparts. The similar dominance of errant polychaetes could be seen invariably in the works done in the mangroves of other Asian countries.

Conclusion
Based on the foregoing account, it is concluded that the present study yielded quite a good amount of information on the benthic biodiversity in general and polychaete taxonomy in particular in the mangroves of east coast of India. As there was no comprehensive report on the polychaetes of mangroves of east coast of India, comparison was done only based on the available sporadic reports and thus a clear -cut inference could not be drawn.
On the other hand, studies related to taxonomy of benthic fauna is limited as the researchers worldwide did not evince much interest in this line besides the enrolment of a new generation of benthic taxonomists has also been poor in the recent past. There are several reasons for this: (i) indifferent attitudes, both in society and educational systems, and (ii) organisms that are "invisible" from the perspective of immediate economic and medical interest to man and more importantly poor funding from the Government. To achieve this, an intensive collaboration of benthic researchers among the Asian countries is need of the hour, as it will throw an important beam of light on the Polychaete taxonomy in the mangroves with a view to formulate management strategies and also to arrive at meaningful conclusions for the policy makers.