First Assessment of the Benthic Meiofauna Sensitivity to Low Human-Impacted Mangroves in French Guiana

: Bioindicators assess the mangroves ecological state according to the types of pressures but they differ with the ecosystem’s speciﬁcities. We investigated benthic meiofauna diversity and structure within the low human-impacted mangroves in French Guiana (South America) in response to sediment variables with various distances to the main city. Contaminant’s concentrations differed among the stations, but they remained below toxicity guidelines. Meiofauna structure (Foraminifera, Kinorhyncha, Nematoda) however varied accordingly. Nematode’s identiﬁcation brought details on the sediment’s quality. The opportunistic genus Paraethmolaimus (Jensen, 1994) strongly correlated to the higher concentrations of Hg, Pb. Anoxic sediments were marked by organic enrichment in pesticides, PCB, and mangrove litter products and dominance of two tolerant genus, Terschellingia (de Man, 1888) and Spirinia (Gerlach, 1963). In each of these two stations, we found many Desmodora individuals (de Man, 1889) with the presence of epibionts highlighting the nematodes decreased ﬁtness and defenses. Oxic sediments without contaminants were distinguished by the sensitive genera Pseudocella (Filipjev, 1927) and a higher diversity of trophic groups. Our results suggested a nematodes sensitivity to low contaminants concentrations. Further investigations at different spatio-temporal scales and levels of deterioration, would be necessary to use of this group as bioindicator of the mangroves’ ecological status.


Introduction
Mangroves are among the most productive tropical environments. Over the past two decades, 35% of the world's mangrove was degraded and lost [1], mainly due to deforestation, erosion, urbanization/pollution, shrimp aquaculture, and tropical cyclones [2]. Mangroves of tropical developing countries with a high demographic pressure are the most degraded [3][4][5]. Key ecological functions are thus weakened (i.e., feeding areas, nurseries, blue carbon, sediment and contaminants retention [6,7]). Many countries need to develop policies for improving mangroves conservation and sustainability through management programs. Many mangroves restoration strategies failed, however, caused by multiple factors, e.g., lack of knowledge on the mangroves functioning and its spatial and temporal variability, conducive to a low selection of sites and to monospecific plantations that reduced animal and vegetal species interactions and diversity of ecological functions, and lack of incentives policies for a long-term management of restored areas that impeded strong mangroves monitoring [8][9][10]. In addition to anthropogenic pressures, it is now widely admitted that mangroves will suffer to the global change (i.e., increases in surface temperature and sea level, changes in precipitation regimes and in coastal sedimentation) that will impact the latitudinal mangrove ranges, vegetal, and animal species composition and the mangroves productivity [11,12]. Objectives of mangrove restoration may thus vary depending on sources of damage to mangroves that includes knowing the diagnosis of the causes of the deterioration and setting a baseline [13]. Bioindicators are among the proposed tools for assessing mangroves ecological state according to the types of pressures and the levels of deterioration but they differ with the ecosystems specificities [14,15]. The bioindicators approach has highlighted recently the needs to develop integrated studies of the mangroves ecosystem in preparation of its successful management [16][17][18]. Such an approach could be applied to some impact studies and management programs, in order for example, to reach the 2030 Agenda for Sustainable Development and the Paris Agreement, or the European Water Framework Directive which ones does not include yet specific Biological Quality Elements as the mangroves.
As an example, in temperate areas, benthic bioindicators have been widely used to monitor human disturbance as well as to establish management targets: Indicator species, diversity indices (Shannon-Wiener index, Pielou index, or Simpson index), ecological groups (AZTI Marine Biotic Index-AMBI, Ecological Quality Ratio-RQE) or trophic groups (Trophic Index of Endofauna; [19]). These indices are based on benthic macrofauna (>1 mm) composition and diversity, are easier to identify [20], but not on the benthic meiofauna ones (<1 mm). Several studies demonstrate that meiofauna (32 µm-1 mm), besides being fundamental to understand the structure and functioning of marine communities, can be used as an additional proxy for responses of benthic communities to environmental changes and can be a useful tool to investigate anthropogenic impacts, reflect spatial and temporal changes, and describe and document the good ecological status of soils and sediments ( [21] and references therein). Meiofauna is the most abundant and diversified benthic component into mangroves sediments providing important ecosystem benefits such as sediment bioturbation activities stimulating the recycling of organic matter [22][23][24][25] and trophic food webs [26,27]. Meiofauna abundance and community structure are shaped by sediments physical and biogeochemical characteristics and are rapidly affected by environmental changes, making this benthic component potentially suitable for biomonitoring [21]. Once mangroves are widely degraded, recovery of mangroves meiofauna communities occurs after 5 to 10 years reforestation [28]. It has been shown that meiofauna identification, even on a low taxonomic level (family, class or phylum) could provide sufficient information to characterize the community and the ecosystems health [20,29], especially in mangrove sediments [30]. A recent review about the use of meiofauna to assess environmental impact proposed four taxonomic groups sensitive to different kind of disturbances: Foraminifera, Nematoda, Copepoda, and Ostracoda [21]. Although some taxa are negatively impacted by disturbances such as Copepoda, others are favored by disturbed conditions such as Nematoda [31][32][33]. Knowing the importance of meiofauna in mangrove benthic ecosystems and their evidence as bioindicator in other marine systems [34], mangrove meiofauna, dominated by nematodes, could be potentially a good indicator of environmental changes in this peculiar ecosystem.
Among all these groups of the meiofauna (Nematodes, Foraminifera, Copepoda, Ostracods), nematodes are doubtless the most abundant metazoan group reaching impressive density in mangrove sediments. Furthermore Foraminifers and copepods are very sensitive to changing environmental conditions, while some nematode species are particularly tolerant to stress or they show spectacular physiological adaptations [21,34]. Thus, mangrove nematodes are known from a large variety of other different ecosystems worldwide, characterized by large fluctuations in environmental conditions [34]. If any nematode endemic genera are reported for mangroves, some genera are considered typically dominant in mangroves worldwide [34,35]. Previous studies have shown that nematode assemblages differ according to the type of disturbance due to certain nematode-specific tolerances [36,37]. In particular, nematode diversity is usually suppressed under environmental disturbances while some species become dominant [38]. Concerning metal pollution such as lead and zinc, studies on an environment with high concentrations of contaminants, exceeding tolerance levels, have been carried out [38][39][40][41]. Specific mangrove nematodes are thus potentially present reflecting different environmental conditions, including most impacted ones.
This study aimed to initiate works on the meiofauna bioindicators, firstly by assessing ecological status of the low human impacted mangroves in French Guiana (FG). French tropical coasts are composed of more than 100,000 hectares of mangroves, and almost 70,000 hectares are located in French Guiana [42]. FG coastal mangroves, dominated by Avicennia germinans, are very dynamic because they are exposed to important sedimentary inputs from the Amazon River generating an alternation of accreted and eroded mud banks, which determine the mangrove development or destruction. On the contrary, FG estuarine mangroves show different dynamics and structure with a much higher stability because of their distance to the mobile mud banks [43,44]. These estuarine mangroves are represented by simultaneous presence of A. germinans and Rhizophora spp. at the adult stand-age, spreading upstream at several kilometers inland along the riverbanks to the limit of tidal influence being the polyhaline area [45,46]. FG mangroves showed, at first sight, only a small anthropogenic disturbance [5] although no major study ever detailed their level of deterioration or pollution state [47]. These mangroves appear in a good initial state for impact studies before being threatened. Indeed, estuarine mangroves near major cities like Cayenne will be impacted by urban waste waters and rapid increase of industrialization and urbanization in the future. As an example, in 2017, the Cayenne's population was 61,268 habitants (2596,1 hab.m −2 ) with a birthrate of 28‰ since 2012 highlighting a rapid demographic increasing over short period [48]. In the current study, a complete set of environmental variables (physical and chemical sediment variables, organic contaminants, metals and metalloids, sediment organic biomarkers, prokaryotic biomasses) were identified and quantified at three stations progressively distanced from Cayenne city, along the Cayenne estuary. For each station, density, biomass, composition, and diversity of the mangrove meiofauna communities and specifically Nematoda were simultaneously identified and compared between the three stations in responses to mangroves sediment characteristics.

Study Area
The study area was located in the Cayenne estuary, French Guiana (South America, Figure 1). The climate is subtropical with two major rainfall periods extending from late March to early July and from mid-November to early March. Tides are semidiurnal and mesotidal with a mean tidal range of 1.68 m. FG is part of the French oversea departments and regions and is characterized by a 320 km-long mangrove coast, largely dominated by A. germinans, Rhizophora spp., and Laguncularia racemosa (C.F. Gaertn., 1807). A. germinans constitutes the often monospecific and homogeneous in adult ages stands, which characterize most of Guiana's coastal mangrove [45]. The three mangrove species spread upstream several kilometers inland, following the riverbanks, to the limit of tidal influence [45], and characterize the estuarine mangroves. In FG, the low slope of the coast enhances large salty areas at the estuary mouths where mangroves grow several kilometers inland. Despite the large amount of Amazonian waters, the FG coast exhibits high salinity [46]. In fact, the freshwater itself remains far from the shore during the dry season, following the north Brazilian current, while only the mud plume arrives at the coast. The three stations investigated in this study were located on the edges of the estuary (<5 m from the river channel) colonized by the adult mangroves, at an increasing distance from the Cayenne city ( Figure 1). The three stations exhibited the same mean sediment elevations (1.02 ± 0.08 m) as revealed by the Digital Elevation Model ( Figure S1A). They were all characterized by a similar adult mangrove stand-age with the same mean canopy height above the terrestrial reference ( Figure S1B), and all of them were dominated by A. germinans and Rhizophora spp., being proxies of constant salinity range over estuary at a similar tidal time [49]. station 1 was located near a wastewater treatment plant, which drained the waters from industrial, commercial, and urban activity at "Crique Fouillée" coming out of Cayenne city (4 • 54 53.208 N; 52 • 20 15.9324 W, Figure 1) [47,50]. Station 2 was located at the intersection between two rivers, Cayenne and Montsinnery, near to agricultural parcels (4 • 53 49.2288 N; 52 • 22 27.714 W, Figure 1). Station 3 was established 10km from the estuary mouth and from the agricultural and urban environments, located upstream the main "Cayenne" River, along the "Petit Cayenne" River (4 • 51 31.9716 N; 52 • 23 59.5248 W, Figure 1), making it the reference station. These stations were all located within the estuarine polyhaline zone, with the same salinity variations ranges as indicated by previous studies [46]. At the time of our sampling, at the rising tides of the spring tides during the dry season, the stations, all located in the same polyhaline area, presented a similar range of salinity for surface water comprised (p > 0.05) between 30 ± 2 (station 1) and 28 ± 2 (station 2) to 25 ± 3 (station 3).

Field Sampling
The stations were sampled between 19 and 21 of November 2017, during the spring tides, two hours before the high tide. High tide made the sampling stations easily accessible. However, the sampling sites were not inundated at high tide. Only station 2 presented thin water layer above sediment, which did not prevent the sampling. Three sediment cores (plexiglass tubes of 10.4 cm internal diameter, 20 cm height each) were manually collected from each station, each core being distanced of about 2m. Sediment cores were sliced horizontally from 0 to 2 cm (L1), 2 to 10 cm (L2), and 10 to 16 cm depth (L3), making three main sediments layers (L1, L2, L3) for each sediment core. Redox potential (Eh) and pH were measured directly in the core for each sediment layer using a multi-parameter system coupled to soil specific probes WTW Multi 3500i. Instrument precision was ±0.01 Eh and ±0.001 pH, the latter calibrated with NBS Buffers. Within each sediment layer, six sediment subsamples were extracted for the analysis of grain size, total carbon (TC) and nitrogen (TN), inorganic and organic contaminant concentrations, pigments content, prokaryotic biomass, and meiofauna, respectively. Inox spatulas and pyrolyzed glass bottles (450 • C, 3 h) were used to sample sediments for analysis of organic contaminants. Plastic spoons and 15 mL falcon tubes were used to sample sediments for metal analyses. Subsamples for grain size analysis were kept in the fridge (4 • C), whereas sediments for biogeochemistry analyses were directly frozen (−80 • C). Subsamples with a sub-surface area of 1.77 cm 2 were taken in each sediment layer for meiofauna analysis. The meiofauna samples were fixed in 4% diluted formaldehyde buffered with sodium bicarbonate and stored cold until analysis. main "Cayenne" River, along the "Petit Cayenne" River (4°51′31.9716′′ N; 52°23′59.5248′′ W, Figure 1), making it the reference station. These stations were all located within the estuarine polyhaline zone, with the same salinity variations ranges as indicated by previous studies [46]. At the time of our sampling, at the rising tides of the spring tides during the dry season, the stations, all located in the same polyhaline area, presented a similar range of salinity for surface water comprised (p > 0.05) between 30 ± 2 (station 1) and 28 ± 2 (station 2) to 25 ± 3 (station 3).   Sediments were freeze-dried (24 h) and gently crushed to powder and homogenized for bulk sediment analyses. Sedimentary grain size distribution was determined using a laser beam diffraction analyzer (Partica LA-950V2; Horiba Instruments, Inc.). Total Carbon (TC) and Total Nitrogen (TN) were analyzed by combustion at 930 • C using a CHN carbon analyzer (Flash-2000; Thermo Fisher Scientific Inc., Milan, Italy). Due to the lack of carbonates in the FG region [51], the organic fractions largely dominated or equally corresponded to the total carbon. The total carbon and nitrogen were thus used as proxies of the organic matter, and the molar C:N ratio was calculated to infer a proxy of the refractory versus labile nature of the organic matter.

Pigment Analysis
Lipophilic pigments were analyzed by high-performance liquid chromatography (HPLC) [52]. Mangrove sediments were incubated with 95% methanol (buffered with 2% ammonium acetate) during 15 min, at −20 • C in the dark. Extracts were then filtered with 0.2µm PTFE syringe filters and analyzed within 16h using an Agilent 1260 Infinity HPLC composed of a quaternary pump (VL 400 bar), a UV-VIS photodiode array detector (DAD 1260 VL, 190-950nm), a fluorescence detector (FLD 1260 excitation: 425 nm, emission: 655nm), and a 100µL automatic sample injector refrigerated at 4 • C in the dark. Chromatographic separation was carried out using a C18 column for reverse phase chromatography (Supelcosil, 25 cm long, 4.6 mm inner diameter). The solvents used were: 0.5M ammonium acetate in methanol and water (85:15, v:v), acetonitrile and water (90:10, v:v), and 100% ethyl acetate. The solvent gradient was set [52], with a 0.5 mL min −1 flow rate. Identification and calibration of the HPLC peaks were performed with antheraxanthin, ββ-carotene, canthaxanthin, chlorophyll a, chlorophyll b, chlorophyll c2, diatoxanthin, diadinoxanthin, fucoxanthin, and pheophytin a standards. We identified all detected peaks by their absorption spectra and relative retention times using the Agilent OpenLab software. Quantification was performed using standard calibration curves built with repeated injections of standards over a range of dilutions. Xanthophylls, carotenes, and chlorophyll b and c were quantified at 470 nm, and chlorophyll a and their derivatives as well as pheopigments were quantified at 665 nm. The relative abundance of each pigment (%) was calculated from its respective concentration in the sample (µg mg −1 ).

Metals and Metalloids Analysis
Sediments were digested using an aqua regia and microwave assistant. Then, the obtained solutions were filtered through 0.2 µm (cellulose nitrate, Sartorius) and diluted for analysis. Concentrations of minor/major/trace elements were measured by High-Resolution Inductively Coupled Plasma-Mass Spectrometer (HR-ICP/MS, Element 2, Finnigan) and validated by using SLRS-5 (River water certified reference material, National Research Council of Canada-NRC), LGC (River certified reference), and PACS2 (NRC) for sediments. The analytical recovery was within 10% compared to certified concentration. Particulate Hg was analyzed on aliquots of 50 mg dry sediment using a cold vapor atomic absorption spectrometry (CV-AAS; Leco Ama 254) equipped with a low-level optical cell. The method detection limit, estimated as three times standard deviation of the blank samples, was 2 pg. The certified reference material MESS-4 (National Research Council Canada) was run several times per analytical batch and constantly before starting the measurements, to check the accuracy of the measurements. The measured values were always within ±5% of the recommended values.

Bacterial and Archaeal Abundance
Total DNA was extracted from composite sediment samples representative of each layer (0.25-0.30 gdry sediment) using the PowerSoil DNA isolation kit (MoBio) according to the manufacturer's recommendations. The bacterial and archaeal abundance was measured by real-time quantitative PCR (qPCR), by the determination of bacterial and archaeal 16S rRNA gene copy numbers, respectively [56]. The resulting conditions lead to a quantitative PCR efficiency higher than 99% (R 2 = 0.99). 16S rRNA genes abundances were standardized by the mass of DNA recovered per g dry sediment.

Meiofauna Analysis
Meiofauna organisms were extracted from sediment samples [57]. Particles were conserved in Falcon tube containing a solution of formaldehyde with drops of Rose Bengal. The meiofauna was identified at the highest taxonomic level possible under a binocular lens (Stemi 508 binocular lens) and were measured with a micrometer screw (length and width). 150 nematodes per sample were collected; they were mounted on perma-Forests 2021, 12, 338 7 of 23 nent slides according to the formalin-ethanol-glycerol treatment and identified under an optical microscope [57,58]. Nematodes were classified in groups according to their trophic groups. Four groups were identified based on morphology characters: Absence or fine buccal cavity-bacterivorous; large buccal cavity without teeth-non-selective detritivores; buccal cavity with scraping teeth-grazers; and buccal cavity with large jaws-omnivorous-predators [59,60]. The maturity index (MI) assigned values from 1 for extremely opportunistic nematodes to 5 for extremely sensitive nematodes and according to life-history traits and their adaptability to a new environment [61]. The MI (Bongers, 1990) for a nematode assemblage was calculated as the weighted mean of the individual taxon scores: where v i was the c-p value of taxon I and f i was the frequency of that taxon.
Meiofauna density and biomass were calculated at all sampling stations and along a vertical profile into the sediment. Density was expressed as number of individuals per 10cm 2 . Biomass calculation was done in µg of carbon per 10cm 2 [57].

Statistical Analyses
Multiple factor analyses (MFA), which are extension of principal component analyses (PCA), are devoted to the study of tables in which groups of variables are introduced to the analysis (instead of using all the variables at once; [62]). This method can deal with tables that have many explicative variables and share the same statistical individuals. The variables are weighted so that in case of unbalanced groups of variables, the analysis will not be influenced by the largest group. The weights are identical for the variables of the same group (and vary from one group to another). They are such that the maximum axial inertia of a group is equal to 1. We thus used MFA analyses with the R package "ade4" [63] in order to find the links between several groups (matrix) of variables measured on the same set of statistical individuals (i.e., biological, physical, biogeochemical, and meiofauna variables on the same sediment core, and/or at the same sediment depth). The overall analysis consisted here in introducing several active groups (for instance: Physico-chemical parameters, metals, contaminants, prokaryotic biomass, organic matter, pigments) in a factorial analysis in order to analyze similarities and discrepancies between groups [64]. Two MFAs were carried out, one with environmental factors and the other with the density and biomass of the meiofauna. On the other hand, the studied entities came from the same samples (in the same sample were extracted several data), which allowed us to compare the two MFA between them thereafter. Separated PCAs were then performed on each group of variables: One with environmental factors, one with the density of nematodes, and one with the biomass of nematodes and then a weight equal to the inverse of the first Eigen value of the PCA was assigned to each group before concatenating them in a single dataset. Additionally, in order to explore the weight of a given factor (e.g., sediment depth or the interaction between sampling station and sediment depth), we performed a supervised Between Class Analysis (BCA) by decomposing the total inertia of the MFA dataset according to a given instrumental variable (i.e., sampling stations, and/or sediment depth). Percentage of total inertia explained by the instrumental variable was systematically calculated ( [63]; hereafter called bca-ratio). The above mentioned statistical procedure was previously detailed [65] and was performed (i) on the sediment parameters measured at each sampling stations (i.e., environmental variables: 3 variables (granulometry, pH, redox potential); pigments: 31 variables; prokaryotic biomass: 2 variables (archaeal and bacterial); organic matter quality: 2 variables (C:N ratio and organic carbon content); contaminants: 18 variables (organic and inorganic contaminants); (ii) on the nematodes density and iii) on the nematodes biomass, grouped by their feeding behavior (i.e., bacterivorous, detritivores, grazers, omnivorous predators). A threshold based on the square cosine of the coordinates of the variables was applied for highlighting the most discriminating variables. A Monte Carlo test was systematically performed to test the significance of the BCA ordination.
We have verified the significance of the most discriminating sediment variables, between the stations, by testing the station (station 1, 2, and 3) and depth effects (L1, L2, L3) on them with variances analysis. When the normality (Shapiro-Wilk test) and homogeneity of the variances were met (Bartlett test), Anova (analysis of variance) followed by Tukey (for pairwise comparisons) tests were carried out. Otherwise, non-parametric tests of Sheirer-Ray-Hare (SRH) test with Wilcoxon (for pairwise comparisons) were carried out [66].
Since the meiofauna occurred mainly over the first 5 cm depths, we tested the influence of station on the overall meiofauna diversity variables (density, biomass, trophic, and maturity indexes) by using a non-parametric one-way analysis of variance (Kruskal-Wallis test) following by Wilcoxon (for pairwise comparisons). The same analysis was done to test the station effect on each taxonomic unit of the meiofauna community (Foraminifera, Acarina, Nematoda, Copepoda, Kinorhyncha, Plathyelminths, Ostracoda).
A permutational multivariate analysis of variance, using a Bray-Curtis dissimilarity index (PERMANOVA; [67]), based on density and biomass per matrix of nematode genus group, was computed to test the multivariate response of infauna assemblages to the station location after verifying the multivariate homogeneity of group dispersions. Correlation analyses were then used to quantify the relatedness between nematodes genus and environmental variables, using the non-parametric Spearman coefficient ρ as data normality was not verified. These analyses were done using the software R (version 3.5.1) and Rcommander.

Environmental Characteristics
The first axis of the Multiple factor analyses (MFA), based on the environmental variables at each depth and station, explained 36.6% and the second axis explained 19.7% of the variance (Figure 2A). MFA analysis showed that the three layers of each station (L1, L2, L3) were grouped together ( Figure 2A). Then, the first axis of the supervised Between Class Analysis (BCA) on the same environmental variables dataset explained 75.6% and the second axis explained 24.4% of the variance ( Figure 2B). station 1 was separated from the other stations by the axis 1 and stations 2 and 3 were separated by the axis 2 ( Figure 2B). The supervised analysis (BCA; Figure 2C) showed that 39.8% of total inertia (i.e., bcaratio) was explained by the sampling stations. The first axis was defined by variations in organic contaminants and metals types, prokaryotic biomasses (archaeal versus bacterial), granulometry (silt and clay versus sand), and organic matter (C:N) while physico-chemical parameters (pH and redox potential) and pigments (Pheophorbide a-Pda versus all other pigments) defined the second axis.
Among these physical and biogeochemical sediment variables defining the two axis, we have defined the most discriminant ones with a square cosine of the coordinates superior to 0.4 ( Figure 2C; cos 2 > 0.4) and the most significant ones ( Figure S2; p < 0.005). Thus, we have determined that the significant higher proportion of sand and the higher contents of metals (i.e., mercury (Hg), molybdenum (Mo), bismuth (Bi), cadmium (Cd), arsenic (As), and lead (Pb)), number of copies of the bacterial 16S, total organic carbon (TOC), and C:N ratio characterized station 1 (Table 1, Figure 2C, Figure S2). Station 2 was distinguished by higher proportion of clay, negative redox potential (Eh), and pigments increase, i.e., Violaxanthin (Vio), Beta-carotene (Bb), Antheraxanthin (An), Lutein (L) (Table 1, Figure 2C, Figure S2). The significant higher proportions in PCB and pesticides also qualified stations 1 and 2 (Table 1, Figure 2C, Figure S2). Finally, significant higher concentrations of PAHs and specific metals, i.e., Zn, Ba (Table 1, Figure 2C, Figure S2) mostly differentiated station 3. We found the highest proportions of silt and higher quantities of Vanadium (V), Chromium (Cr), and Manganese (Mn) in stations 2 and 3.

Total Meiofauna Community
Total densities of the overall benthic meiofauna community varied from 2651 ± 121 individuals per 10 cm 2 at station 1, to 7431 ± 2141 individuals per 10 cm 2 at station 2 and 6712 ± 1619 individuals per 10cm 2 at station 3 ( Figure 3), but they were not significantly different. However, densities of certain taxa were different among stations. Foraminifera were found in higher abundances at station 2 (Kruskal-Wallis test; p = 0.03) whereas Acarina showed higher density values at station 1 (Kruskal-Wallis tests; p = 0.03). Nematoda was the most abundant taxon at all sampling stations: 76% at station 2, 90% at station 1, and 93.5% at station 3 ( Figure 3).
Total meiofaunal biomass ranged from 210 ± 43 µ gC per 10 cm 2 at station 1, 404 ± 138 µ gC per 10 cm 2 at station 2, and 375 ± 105 µ gC per 10 cm 2 at station 3 ( Figure 3), but they were not significantly different.  Table 1. Mean values of the three cores over the entire sedimentary column from each station for the physical, chemical, and biological parameters (mean ± SD; n = 3).

Total Meiofauna Community
Total densities of the overall benthic meiofauna community varied from 2651 ± 121 individuals per 10 cm 2 at station 1, to 7431 ± 2141 individuals per 10 cm 2 at station 2 and 6712 ± 1619 individuals per 10cm 2 at station 3 ( Figure 3), but they were not significantly different. However, densities of certain taxa were different among stations. Foraminifera were found in higher abundances at station 2 (Kruskal-Wallis test; p = 0.03) whereas Acarina showed higher density values at station 1 (Kruskal-Wallis tests; p = 0.03). Nematoda was the most abundant taxon at all sampling stations: 76% at station 2, 90% at station 1, and 93.5% at station 3 ( Figure 3).
The two first axes of the unsupervised analysis (MFA) explained 24% of the total variance ( Figure 4A). The MFA analysis showed a differentiation of the three stations based on L1 whereas L2 and L3 are grouped. The supervised BCA analysis showed that Total meiofaunal biomass ranged from 210 ± 43 µgC per 10 cm 2 at station 1, 404 ± 138 µgC per 10 cm 2 at station 2, and 375 ± 105 µgC per 10 cm 2 at station 3 ( Figure 3), but they were not significantly different.
The two first axes of the unsupervised analysis (MFA) explained 24% of the total variance ( Figure 4A). The MFA analysis showed a differentiation of the three stations based on L1 whereas L2 and L3 are grouped. The supervised BCA analysis showed that 14% of total inertia was explained by the sampling stations and showed a differentiation of the three stations based on the nematodes genera relative abundance: station 1 along the axis 1 and stations 2 and 3 along the axis 2 ( Figure 4B,C). The cos 2 threshold of 0.2 highlighted the most structuring genus of nematodes ( Figure 4C). The PERMANOVA analysis (not shown here) allows us to affirm that nematodes composition is different between stations (p < 0.001).   Seventy-three nematode genera were observed (all stations combined) and 18 were present at all three stations, including the Desmodora, Terschellingia, Paracomesoma (Schuurmans Stekhoven, 1950), and Microlaimus (de Man, 1880) ( Figure S3).
Thirty-seven nematode genera were identified at station 1, and the Desmodora and Spirinia genera accounted for more than 35% of the total nematode abundance at this station. However, the genus Paraethmolaimus significantly discriminated the station 1 (Table 2, Figure 4C).
At station 2, we identified 45 nematode genera and Terschellingia genus alone accounted for 27% of the nematode abundance at this station. The presence of Terschellingia significantly distinguished station 2 from the two other stations ( Figure 4C, Table 2). Table 2. Mean total nematode densities (individuals per 10cm 2 ; ±SD) and mean relative nematode densities (%; ±SD) at each station (here reported nematodes with at least 2% density on one station). The nematoda exhaustive list including those with a relative density < 2% is presented in Figure S3 with their appartenance trophic groups. In addition, we observed the presence of epibionts on nematodes of the genus Desmodora at stations 1 (i.e., 2.2%) and 2 (i.e., 23.5%).
Spirinia genus was significantly less present at station 3, as well as the Halalaimus genus at station 2, whereas Ptycholaimellus (Cobb, 1920) genus was less observed at station 1 ( Table 2). The same analysis based on the relative biomass of nematodes was carried out and show roughly the same results ( Figure S4).
Regarding the trophic groups of nematodes based on the nematodes densities, grazers were the most represented in all the stations with a higher grazing density at station 1 and predators the least ones (p < 0.05; Figure 5; Figure S3). The trophic diversity index, based on the nematode trophic groups, was higher at station 1 (0.48 ± 0.05) than at stations 2 and 3 (0.36 ± 0.03) (p < 0.03). However, we observed no difference in the proportion of the different trophic groups of nematodes based on the nematodes biomass.
The nematodes maturity was based on their adaptability to a new environment, which makes it possible to classify nematodes in very opportunistic (and therefore very adaptable) genera up to very sensitive (and not very adaptable) genera. The nematode maturity index did not differ between stations (2.5 ± 0.05 to 2.9 ± 0.02). However, we noticed that opportunistic and intermediate nematodes were the most abundant at each station    Relative to the total nematode densities within each station, we observed that stations 1 and 2 were made up of a higher proportion of opportunistic and intermediate nematodes (p < 0.05; Figure 5B).
The differentiation of station 3 took place at the level of the biomass of extremely opportunistic nematodes, which was the lowest (p < 0.02; Figure 5B).
Regarding the differentiation of stations at the level of nematode biomasses per maturity groups, only the biomass values of extremely opportunistic nematodes varied from one station to another. Station 2 showed more extremely opportunistic individuals (p < 0.02; Figure 5B).

Environmental Characteristics of the Three Stations
The overall studied stations, although located within similar mangrove stand age and along the polyhaline estuarine area, have shown distinct and significant different sediment variables. Although the dominant mode of grain size remained the silt for each station over all the sedimentary column, we think, however, station 1 was likely more impacted by the proximity to the seafront since the sand mode was more important at this station. Both stations 1 and 3 were dominated by aerobic redox conditions (Eh~100-200 mV) while the sediments of station 2 showed negative dominant anaerobic patterns with negative redox potential (Eh~−100mV), an indicator of the less oxidized sediments. We have found a dominance of the mangroves litter pigments (i.e., Antheraxanthin, Lutein, Betacarotene, and Violaxanthin) at station 2, thus highlighting an organic enrichment. On the contrary, although C:N ratio of 16-20 indicated more refractory material at station 1, the 16S DNA increase was the sign of numerous bacterial cells and suggested a higher bacterial OM degradation and mineralization at this station. However, archaeal abundances were about 10 (stations 2 and 3) to 1000-fold (station 1) lower than those of bacteria (Table 1), but they were not considered as a discriminant variable (low cos 2 < 0.4) in the multiple factor analysis (MFA) and between class analysis (BCA), and our interpretation on the prokaryotic community remained based on the bacterial abundances.
We have noticed a strong presence of pesticides and PCB at stations 1 and 2 that were relatively absent from station 3. Such organic contaminants, likely anthropogenic, might come from the wastewater treatment plant ("Crique Fouillée") located near station 1, and from agricultural parcels and by river runoff having a greater impact on station 2 [47,50]. Among the organic chemical components, we have also measured 4 to 5-fold higher concentrations of PAHs at stations 2 and 3 than station 1. By drawing a comparison with polluted mangrove areas where PAHs concentrations levels are up to several thousand ng/g [68], we have concluded that PAHs concentrations measured in the Cayenne estuary remained low. We think most of the PAHs found in coastal sediments of French Guiana are probably from pyrolytic origin (wood burning [69], maintenance of coastal savannas [21]) due to a low industrialization and a still-low human population density. In our case, we think that PAHs can more easily accumulate in finer sediments (station 2 and 3) but also when sediment organic content is higher (station 1) [70]. Here, differences were most probably due to differences in the number of atmospheric inputs and, perhaps, to a somehow easier accumulation in the sediments.
By referring to the official guideline concentrations for assessing toxicity hazards (i.e., effects range low (ERL) and effects range median (ERM; NOAA, 2012)), we have seen that organic contaminants and metals concentrations remained low. The most toxic PCB in our analyses were dioxin-like polychlorinated biphenyl 169 with a total concentration of 19.1 µg.kg −1 being just below the official used threshold [71] which does however depend on various intrinsic sediments parameters as the OM content. The most abundant pesticide was hexachlorobenzen found at station 2 at a concentration below the level of 17 ng/g that is calculated to be a predicted no effect HCB concentration for chronic exposure (PNEC chronic ) for organisms living in sediments [72]. We think that the different natural physical and biogeochemical context between the stations explained the high variations of metal concentrations. As example, we know that metals origin might be lithogenic with a succession of co-precipitation and co-complexification geochemical processes as already observed in another Guianese mangroves area [73]. However, the metals Pb, Hg, Mo, and Bi clearly characterized the station 1 likely due to the proximity to a wastewater treatment plant ("Crique Fouillée"). In this station, Pb (39.95 mg/kg) and Hg (0.64 mg/kg) concentrations were just above the thresholds from which metals can potentially affect benthic invertebrates but below lethal thresholds [74][75][76]. Station 3, located downstream of the Cayenne estuary still in the polyhaline area, received less contaminants and its more oxidized sediments revealed reactive biogeochemical functioning. At a similar tidal level to the station 1, metals as well as organic contaminants were less present at station 3 due to its distance from the Cayenne city, urban wastewater, and agricultural impacts. We thus think, as expected, this reference station is more preserved than the other two stations from anthropogenic effects.
Interestingly, we found the presence of a degradation pigment (Pheophorbide a-Pda) at station 3 and pheophytins a at station 2. The first step of chlorophyll a degradation in plants is hydrolysis of the phytyl ester bond by chlorophyllase to form chlorophylide a and phytol. In the second step, chlorophylide a is demetallized to form Pheophorbide a [77]. Pheophorbide a is thus a degradation product of chlorophylide a while pheophytin a is a degradation product of chlorophyll a. The presence of these two pigments could indicate different degradation processes between the two stations (stations 2 and 3) either in terms of metabolism or in terms of source of organic matter for invertebrates. Chlorophyll degradation could be advantageous in plants by preventing the accumulation of excited free chlorophyll or if protein amino acids must be redeployed for alternative needs in case of nutritional deficiency [78].

Overall Meiofauna Community
Meiofauna is generally more resistant than macrofauna and increases in density in the presence of organic and industrial contaminants [79]. However, this is true until a threshold level, after that a decrease in density is usually observed [80]. A previous study [81] had highlighted that mangrove meiofauna abundance is higher in undisturbed sediments compared to sediments disturbed by anthropogenic activities (e.g., from 2.6 ± 1.1 to 1.6 ± 4.41 individuals per 10 cm 2 ).
In our study, meiofauna diversity did not vary significantly between stations. Nevertheless, its community structure was significantly different among the stations. Higher Foraminifera density at station 2 confirmed the presence of less oxygenated sediment [82][83][84]. In low oxygen concentrations, Foraminifera set up specific mechanisms, in particular decreasing their size [85,86]. In the present study, Kinorhyncha reached almost 5% (station 2) of the overall meiofauna, which is in the range of values reported from stressful environments (e.g., wide fluctuations of salinities [87][88][89], while they are usually less than <5% in shallow waters and deep-sea ecosystems [30].

Nematodes Community
The deeper identification we made on the nematodes genera has brought more detailed indications of sediment conditions, as also confirmed by other experimental works [33]. Previous studies shown that salinity differences among the four haline zones of an entire estuary (i.e., euhaline, polyhaline, mesohaline, oligohaline; [90]) can control the distribution of the nematoda assemblages [91]. In our case, we excluded the salinity effect since the three stations were all located within the same polyhaline zone of the estuary. We sampled during the dry season, but during the wet season the estuary salinity decrease [92] might impact the nematoda assemblages over a year for the same stations, as observed in some Vietnamese estuaries [93]. In a study based on the metals impacts (e.g., Hg, Pb, As, Cd, Cu, and Cr) on meiofauna assemblages [40], authors supported the evidence of some nematode tolerance against metals loads whose concentrations exceeded the National Environment Quality. However, in our study, lower trace metal concentrations, supposed to be mainly of lithogenic origin, might affect nematodes community. At station 1, the genus Paraethmolaimus was dominant and significantly correlated to the presence of Hg, Pb, Cd, and Mo. In fact, recent work has underscored nematodes strong sensibility to low metals concentrations [94].
In our study, the genus Terschellingia accounted for 27% of the nematodes at station 2 ( Table 2) where it was correlated with the presence of PCB contaminants, but also a site with the highest densities and biomass of opportunistic nematoda. The dominance of a nematode genus is typical during a disturbance, with tolerant genera taking over sensitive organisms. The genus Terschellingia is known to be highly tolerant to organic contaminants [95], and particularly to PCB [96], and accounts for more than 10% of the nematode community in disturbed stations [97,98]. In addition to being dominant in organic-rich sediments [99], this genus is known to be adapted to the low oxidized muddy environments [100], which was the case for station 2 as indicated by the highly negative redox potential. Although high abundance of Terschellingia can be used as an indicator of polluted environments [101], other studies showed that this genus was unable to adapt to high metal concentrations [39] which might be the case in station 2 where metals concentrations remained low. Even if Terschellingia spp. tolerates organic enrichment, their activity levels and functions in benthic habitats can decline under oxygen limitation [102]. This case illustrates how the species in this Genus, characterized by different life cycles, respond to a wide range of organic contamination/enrichment regimes.
At station 3, the five abundant genera (Southerniella Allgén, 1932; Pseudocella, Deontolaimus, Halalaimus, and Anoplostoma) accounted for more than 50% of the total nematode abundance. The genera Pseudocella and Halalaimus are not tolerant to environmental disturbances and variations [61]. In our study, the genus Pseudocella accounted for 10% of nematode abundance. Thus, its presence indicated an absence or weak environmental disturbance at station 3. In addition, we noticed a low presence of the genus Spirinia, which demonstrated its stress tolerance, physiological, and behavioral adaptability to environmental disturbances in previous studies [36,99,103]. As the sediments of stations 3 were well oxidized without organic-rich material neither contaminants, we concluded that station 3 presented stable environmental conditions for those sensitive and highly sensitive nematodes.
The Desmodora genus was present at the three stations in similar proportions. This genus is known to adapt to different environmental conditions [95,99]. Its presence in the three stations can be explained by its adaptability. However, we observed Desmodora with epibionts only at stations 1 and 2 ( Figure 6). These epibionts, classified as suctorian ciliates [104], were found almost exclusively on the genus Desmodora: station 1 has lower epibionts values (2.2%) than station 2 (23.5%). A co-occurrence between epibionts and organic contaminants may support the fact that this epibiosis could be the result of a decrease in the immunity of organisms induced by stress, such as disrupted environmental conditions. Ciliate epibiosis is a common phenomenon in marine nematodes. As recently reported [104], the presence of epibionts on nematodes cuticle (which is involved in vital functions such as exchange of gases and nutrients, info-chemicals, and defense metabolites) can have a high cost for the host. Two to three epibionts already cause stress for the basibiont whereas more than five epibionts increase the energetic demand for locomotion [105]. Furthermore, it has been shown that in a variety of marine ecosystems, the presence of epibionts on nematodes was reported only in the most stressful condition sites (both due to anthropogenic actions or naturally extreme conditions) [104]. The occurrence of epibionts in a mangrove ecosystem could be due to habitat heterogeneity and the nature of anthropogenic disturbance [106]. Environmental stressful conditions can promote epibiosis because nematodes are more vulnerable and easily infested ( [104]; and references therein) promoting the use of epibiosis as ecological indicator in the assessment of ecological mangrove status. values (2.2%) than station 2 (23.5%). A co-occurrence between epibionts and organic contaminants may support the fact that this epibiosis could be the result of a decrease in the immunity of organisms induced by stress, such as disrupted environmental conditions. Ciliate epibiosis is a common phenomenon in marine nematodes. As recently reported [104], the presence of epibionts on nematodes cuticle (which is involved in vital functions such as exchange of gases and nutrients, info-chemicals, and defense metabolites) can have a high cost for the host. Two to three epibionts already cause stress for the basibiont whereas more than five epibionts increase the energetic demand for locomotion [105]. Furthermore, it has been shown that in a variety of marine ecosystems, the presence of epibionts on nematodes was reported only in the most stressful condition sites (both due to anthropogenic actions or naturally extreme conditions) [104]. The occurrence of epibionts in a mangrove ecosystem could be due to habitat heterogeneity and the nature of anthropogenic disturbance [106]. Environmental stressful conditions can promote epibiosis because nematodes are more vulnerable and easily infested ( [104] ; and references therein) promoting the use of epibiosis as ecological indicator in the assessment of ecological mangrove status. We classified the majority of nematodes as epistrate feeders, as previously reported [107]. The high grazers abundance in our study has indicated the presence of algae and vegetal debris accessible to nematodes. This is typically the case of mangroves where organisms have direct access to mangroves leaves and products of their degradation. We also observed that most of the representative nematodes were grazers [107], at the stations with a higher organic-rich sediment, highlighted by the dominance of mangroves litter pigments and of the degradation pigments. Station 2 was dominated by selective deposit feeders (e.g., Terschellingia) and epistrate feeders (e.g., Desmodora, Spirinia, Microlaimus). The presence of selective and non-selective deposit feeder nematodes indicated the presence of organic detritus that could result from the release of wastewater into the environment [108]. Opportunistic and Intermediate nematodes were the most abundant at all sampling stations. These life strategies are common between nematodes genera that are able to colonize impacted and unstable environments but also undisturbed ones. Our data can be partially explained by the high natural environmental variability typical of mangroves. A previous study [109] also showed similar maturity index values in both disturbed and preserved environments. However, the fact that station 3 hosted higher abundances and biomass of very sensitive nematodes confirmed that impacts in this station were lower allowing the surviving of sensitive organisms. We classified the majority of nematodes as epistrate feeders, as previously reported [107]. The high grazers abundance in our study has indicated the presence of algae and vegetal debris accessible to nematodes. This is typically the case of mangroves where organisms have direct access to mangroves leaves and products of their degradation. We also observed that most of the representative nematodes were grazers [107], at the stations with a higher organic-rich sediment, highlighted by the dominance of mangroves litter pigments and of the degradation pigments. Station 2 was dominated by selective deposit feeders (e.g., Terschellingia) and epistrate feeders (e.g., Desmodora, Spirinia, Microlaimus). The presence of selective and non-selective deposit feeder nematodes indicated the presence of organic detritus that could result from the release of wastewater into the environment [108]. Opportunistic and Intermediate nematodes were the most abundant at all sampling stations. These life strategies are common between nematodes genera that are able to colonize impacted and unstable environments but also undisturbed ones. Our data can be partially explained by the high natural environmental variability typical of mangroves. A previous study [109] also showed similar maturity index values in both disturbed and preserved environments. However, the fact that station 3 hosted higher abundances and biomass of very sensitive nematodes confirmed that impacts in this station were lower allowing the surviving of sensitive organisms.

Conclusions
We showed that the use of nematodes can provide accurate information to characterize the quality of the mangroves environment rather the investigation of the overall meiofauna community. An accurate identification of Foraminifera and Kinorhyncha in FG also should help in the future. This study confirmed that the structure of nematodes community varied according to the sediments characteristics. The presence of organic chemical components (PAHs, pesticides, PCB) and metals, even at low concentrations, in addition to redox state of the sediments, might influence meiofauna composition. The different concentrations of metals found at the three stations, which were below ERM and for most of the metals close to ERL, still resulted in a response from meiofauna and nematodes. Our results show that Paraethmolaimus, Terschellingia, Spirinia, and Desmodora and associated epibionts might be the first indicators of a pollution onset, while the sensible genus Pseudocella could be an indicator of absence or very low contamination rate. This information needs to be acquired in areas where human impacts are absent or low and before environmental threats are present. These data sets and information can inform the preparation of pollution management programs in case where the input of sediments with high organic and metal concentrations are occurring or expected. However, further studies of mangrove meiofauna at different spatial and temporal scales (e.g., seasons) from a large range of geographical locations with different levels of deterioration are essential for strengthening our conclusions on the use of Nematoda as bioindicator taxon of mangroves ecological status. Because the response of Nematoda can differ according to the type of disturbance, future work dedicated to mangrove management programs or restoration strategies should also include speciation analysis of metals (e.g., methylmercury, As (III/V)) and organic chemical components (e.g., pesticides, PCB) within mangroves sediments and surrounding waters.
Supplementary Materials: The following are available online at https://www.mdpi.com/1999-4 907/12/3/338/s1. Figure S1: (A) Digital elevation model of the studied area. Data were obtained using a Digital Elevation Model (DEM) distributed by "Guyane SIG -Plateforme Territoriale" which can be retrieved at https://catalogue.guyane-sig.fr/mapfishapp/map/9262ecd32a82ebc7308129df3 0603b2d. The DEM was derived from a LIDAR survey made in 2015 with a ground pixel size of 1 m 2 . (B) For each of the three study sites, 100 elevation data points were extracted inside square plots and basic statistics were conducted (Min, Max, Mean and STD to the mean canopy elevation). Figure S2: Differences in the most discriminant physical, chemical and biological sediment variables with a cos 2 > 0.4, at the three stations of Cayenne estuary (S1, S2, S3) for each sediment depth (L1, L2, L3) after parametric two-way ANOVAs (F and p-values) or non-parametric Scheirer-Ray-Hare (H and p-values) tests depending on the homogeneity of variances. Significant differences (p < 0.05) are indicated by a different letter (a, b, c). Last column indicates the pairwise comparisons between the three stations. Figure S3: Comprehensive list of the Nematode families and genus found at the three stations of the Cayenne estuary. This table presents the trophic groups for each genus (1A: bacterivorous; 1B: non-selective detritivores; 2A: grazers; 2B: omnivorous-predators) and their mean density (individuals per 10cm 2 ; ± SD) and mean relative densities (%; ±SD) at each station. Figure S4: Multiple Factor Analysis (MFA; A) and Between Class Analysis with a cos 2 threshold of 0.2 (BCA; B: stations distribution and C: variables distribution) based on relative percentages of nematodes biomasses: distribution of nematodes and their trophic group at each station (black: bacterivorous; red: non-selective detritivores; green: grazers; blue: omnivorous-predators). Figure S5: Table of correlation coefficients (Spearman) between significant sediment and nematodes variables discriminating the stations.