Next Article in Journal
Identification, Elucidation and Deployment of a Cytoplasmic Male Sterility System for Hybrid Potato
Previous Article in Journal
Adenylyl Cyclase in Ocular Health and Disease: A Comprehensive Review
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 13
Issue 6
10.3390/biology13060446
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Travelling in Microphis (Teleostei: Syngnathidae) Otoliths with Two-Dimensional X-ray Fluorescence Maps: Twists and Turns on the Road to Strontium Incorporation

by
Clara Lord
1,*,
Vincent Haÿ
1,
Kadda Medjoubi
2,
Sophie Berland
1 and
Philippe Keith
1
1
UMR 8067, Biologie Des Organismes Et Écosystèmes Aquatiques (BOREA), Sorbonne Université, Muséum National d’Histoire Naturelle, Université de Caen Normandie, Université Des Antilles, CNRS, IRD, CP26, 43 Rue Cuvier, 75005 Paris, France
2
Synchrotron SOLEIL, 91192 Saint-Aubin, France
*
Author to whom correspondence should be addressed.
Biology 2024, 13(6), 446; https://doi.org/10.3390/biology13060446
Submission received: 16 May 2024 / Revised: 4 June 2024 / Accepted: 7 June 2024 / Published: 18 June 2024
(This article belongs to the Section Zoology)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Otoliths, calcified biomineralised structures in bony fish inner ear, are remarkable objects as keepers of individual life, both in terms of time and environment. Otolith microchemical analysis was undertaken with a synchrotron-based scanning X-ray fluorescence 2D high-resolution mapping method and we developed analytical imaging processes to retrieve biological information. Strontium incorporation in the otoliths of four pipefish species from Pacific island streams was analysed, as its presence in otoliths is a sign of life in a marine environment, while its absence corresponds to life in freshwater. Their life cycle is expected to be diadromous, allowing for species sustainability in labile and fragmented habitats. From a biological point of view, we uncovered the existence of plasticity in the diadromous life cycle, enhancing our knowledge of these poorly studied species, which are highly threatened by anthropogenic actions. In addition, our 2D hyperfine-scale images outlined the heterogeneity of strontium distribution in the otolith to a greater extent than what is generally assumed. Our results led us to uncover a trade-off in expression between the biological and chemical driving forces that govern element incorporation in the otoliths. Our methods and results provide a strong basis for future works and prospects in fish biology, ecology, and conservation.

Abstract

Indo-Pacific tropical island streams are home to freshwater pipefish (Microphis spp., Syngnathidae). Otoliths were used to uncover life history traits in four species, including a New Caledonian endemic. All four species present the same methodological challenge: their otoliths are small, fragile and mute for growth marks using basic observation tools. Strontium (Sr) is calcium substituent in the mineral lattice, driven by salinity conditions, and thus useful to study diadromous migrations. Synchrotron-based scanning X-ray fluorescence 2D high-resolution mapping allowed us to tackle the global and hyperfine strontium (Sr) distribution. We developed analytical imaging processes to retrieve biological information from otoliths from the data generated via synchrotron analysis. We uncovered plasticity in the life cycle: all species were amphidromous, apart from some freshwater residents from New Caledonia. Understanding life cycle modalities is crucial to categorize species distribution limits and to implement adapted conservation measures, especially when endemic species are at stake. 2D fine-scale images outlined the heterogeneity of Sr distribution: in addition to the trivial Sr incorporation driven by environmental ionic conditions, there is an unusual mosaic arrangement of Sr distribution and we hypothesize that biological control, especially growth during the early life stages, may sometimes overrule stoichiometry. This shows that it is worth studying otolith formation and element integration at imbricated scales, and our methods and results provide a strong basis for future works and prospects in otolith science.

1. Introduction

Teleost fish otoliths are paired calcium carbonate structures located in the inner ear. Approximately 98% of the otolith is a CaCO3 mineral, usually in the mineral form of aragonite, while the other ~2% is composed of an organic matrix and minor and trace elements [1].
Otoliths have long been used as recorders of individual life as they grow throughout the fish’s life and are not subjected to remodelling [2]. Their incremental growth allows for time estimations; minor and trace elements are preserved in the otolith matrix and their concentrations can reflect environments with different chemical signatures, such as marine or freshwater environments in which the fish has been. These trace elements may also provide information on the physiological state of the fish [3,4].
The study of otolith chemical composition can be used to reconstruct individual migratory and environmental histories (e.g., [5,6,7]). Indeed, Sr concentrations are considered to be stable and much higher in the marine environment. Otolith Sr content is correlated with an ambient Sr concentration in the surrounding water [8]. The strontium/calcium ratio (Sr:Ca) in the otolith, in relation to the transition between water bodies of different salinities, is commonly reported in the range of 2 to 4 mmol·mol−1 in freshwater and >6 mmol·mol−1 in estuarine and sea water [9,10,11]. This ratio has been commonly used to reconstruct migration patterns in several diadromous species: anadromous salmonid species [9,11], the catadromous genus Kuhlia [12], and the eel Anguilla marmorata, Quoy and Gaimard, 1824 [13], as well as amphidromous Sicydiinae gobies [5,14].
Unlike other diadromous species, the migration between fresh and marine waters of amphidromous species does not involve reproduction. Adults live, grow and reproduce in freshwater; after hatching, individuals drift downstream to sea where they spend a variable amount of time before returning to rivers [15]. In tropical island rivers, fish are mainly amphidromous (especially Gobiidae and Eleotridae) [16] as the highest likelihood for colonisation of such environments is by organisms coming from the sea and able to tackle the transition between the two biomes [17].
Syngnathidae, including seahorses, seadragons and pipefish, are a charismatic teleost family including about 50 genera and 300 species [18]. They are mainly found in tropical marine areas; however, about 30 pipefish species inhabit freshwater [19] and are essentially distributed in tropical island rivers of the Indo-Pacific. Freshwater pipefish are highly endangered species [20] and their protection is impeded by an important lack of knowledge.
The choice to study freshwater pipefish is not only driven by conservation purposes; it is also motivated by the challenge of studying reluctant otoliths. Pipefish have tricky otoliths: very small in size (<400 µm), fragile and lacking any sort of growth marks using classic observation methods [7].
Nowadays, the most commonly used method is 1D transect or point-by-point methods such as laser ablation coupled to an inductively coupled plasma mass spectrometer (LA-ICPMS). This method has many uses such as in the study of the amphidromous life cycle of Sicydiinae gobies [5], stock assessment of the migratory porcupine fish by Sr marking [21] or the study of the migratory behaviour of commercial fish [22]. With this method, researchers can focus on the early growth region [23] or can obtain life history profiles along a transect over a chosen growth axis of the otolith; these types of transects are often used to detect environmental changes during the individual’s life such as migrations between water masses of varying salinity [5,12,24]. However, interpretation of the results with this powerful but destructive technique (since the laser beam ablates matter leaving a few μm to several tens of μm wide trench on the sample) is carried out under the assumption that the analysed zone chosen is representative of the individual life history. However, caution must be taken since physiological jamming and biomineralization processing may lead to spatial heterogeneity in elemental incorporations [1,10,25]. Adding a second dimension, i.e., 2D mapping of elements, allows refining the study of element incorporation patterns. Elemental 2D maps can be obtained using a UV high-repetition-rate femtosecond laser ablation system coupled to a high-resolution inductively coupled plasma mass spectrometer. This method consists in ablating the entire surface of the otolith repeatedly; each transverse section ablated produces an elemental image [26]. The width of the ablation beam is usually 20–30 μm in diameter. Electron microprobe analysis (EPMA) was also used for otolith analysis. It is fundamentally the same as scanning electronic microscopy, with the added capability of chemical analysis. EPMA allows acquiring precise, quantitative elemental analyses at very small “spot” sizes (c. 1–2 microns), primarily by wavelength-dispersive spectroscopy. But EPMA otolith analysis is usually coupled to LA-ICPMS techniques: the EPMA detection limit of elements is quite high (100 μg/g) and only a few elements in the otolith can be analysed; the LA-ICPMS detection limit is around 0.1–1 μg/g. Researchers use these two complementary techniques, with an overlap of elements that can be measured with both, ensuring the correctness of the results [27,28]. X-ray-fluorescence-based methods have also been developed such as proton-induced energy emission spectroscopy (micro-PIXE) or scanning X-ray fluorescence microscopy (SXFM). These two methods create X-ray fluorescence by exciting atoms either with accelerated protons of high energy (PIXE) or high energy X-rays produced by a synchrotron accelerator (SXFM). At de-excitation, atoms emit X-ray photons characteristic of the element present and the number of X-ray photons emitted gives the quantity present. The size of the spots can vary depending on the technique used and has an impact on element quantification [29]. To analyse elements at the microscale, synchrotron X-ray fluorescence (XRF) is a promising method for producing fine-scale 2D elemental maps and can simultaneously produce multi-element maps of biomineralized structures [7,29,30,31,32]. XRF maps allow studying element incorporation dynamics throughout the otolith at a fine scale (presently down to 0.5 μm) [7,9], even for elements present at low concentrations (lower limit of detection: 1–10 ppm; see [33]).
In this paper, we used XRF 2D mapping because of the high resolution offered by this method to understand better Sr incorporation in pipefish otoliths. The underlying biological question is the diadromous life cycle of tropical freshwater pipefish. Sr XRF 2D mapping was recently used on the otoliths of two tropical freshwater species [7], Microphis brachyurus (Bleeker, 1854) and Microphis nicoleae Haÿ et al., 2023 [19]; it was shown that these two species navigate between water masses of different salinity during their life, a strong suggestion that they are amphidromous. In this paper, we further analysed Microphis brachyurus (Sri Lanka to French Polynesia) and M. nicoleae (Papua New Guinea and the Solomon Islands) and added two species, Microphis torrentius Jordan and Seale, 1906 (Samoa, Vanuatu and New Caledonia), and Microphis cruentus Dawson and Fourmanoir, 1981 (New Caledonian endemic). We used hyperfine 2D methods and developed analytical imaging processes to reveal otolith Sr incorporation patterns and to retrieve information on life history traits.

2. Materials and Methods

2.1. Sample Preparation

Otoliths were obtained from the Muséum national d’Histoire naturelle (MNHN, Paris, France) collections: fish were collected during field campaigns and authorization was given by the curators to retrieve the otoliths. Freshwater Syngnathidae were sampled from different locations in the Indo-Pacific realm and preserved in 95% ethanol, a preservation medium that has been shown not to hamper microchemical analysis [34]. We used a total of 32 specimens: 17 Microphis brachyurus, 8 Microphis nicoleae, 3 Microphis torrentius and 4 Microphis cruentus. The standard length of each individual was consistent with syngnathids in adulthood (Figure 1, Table 1).
The otoliths (sagittae) were removed from the saccules of the inner ear of each sample, thoroughly rinsed with MilliQ water and kept dry until they were individually embedded in epoxy resin (Araldite 2020, Escil, Chassieu, France). The embedded otolith was cut on either side of its core along the longest transverse axis using a low-speed diamond-bladed saw (Buehler, Leinfelden-Echterdingen, Germany) and ground-on carbide silicon abrasive discs of decreasing grain size (Escil, Chassieu, France). A slice less than 200 µm thick was obtained exposing a frontal section down to the core. The slice was checked by light microscopy.
Utmost care was taken to prevent impurity deposition on the sample, especially metal contamination, at every step of the process. The resin was also used for mechanical fixation of the sample to avoid movement during the beam time.

2.2. Synchrotron-Induced XRF Scanning and Data Acquisition

X-ray fluorescence spectrometry (XRF) was performed by scanning the full-size area of the otolith cross sections at the Nanoscopium Nanoprobe CX3-scanning spectromicroscopy beamline of Synchrotron Soleil (Saint-Aubin, France).
The beamline (approx. 155 m long) allows focusing high energy and coherent X-ray beam from a U18 X-ray source undulator down to 100 nanometres. The incident monochromatic photon beam was set to 16.2 keV energy and focused by a Kirckpatrick–Baez mirror. The multilength scale capability of the beamline allows for adjusting the resolution of the map from 1 × 1 down to 0.3 × 0.3 (H × V µm2) for our experiment [35]. The sample was mounted on an X and Z translation stage allowing displacement with a positioning accuracy down to 100 nm and with a travel range up to 1 cm. Two silicon drift detectors (SDD VITUS H50, KETEK GmbH, Munich, Germany) were mounted at 120° to the incident X-ray beam to measure the X-ray fluorescence emitted by the sample.
Spot-by-spot data was collected under continuous scanning by using the available FLYSCAN mode at the beamline. In XRF, being a full-spectral technique, the signal was simultaneously collected from all elements that fluoresce under the experimental conditions. Quantitative reconstruction of the data was conducted by an in-house process of raw data optimization by reducing the noise, removing artefact signals and measuring the discrete energy intensities of the target elements [36].
Energy calibration, a crucial step for element identification, is performed in a standard way from each data detector apart by selecting at least two known X-ray lines [37] and specifying the energy confidence limits allowed to contain the characteristic X-ray peak for each element. Then, each element derivation of the XRF spectrum is obtained by summing the energy signals from all the channels over which the characteristic peak of this element extends.
The output dataset is a collection of 2D elemental distribution encoded in Tag Image File Format (.TIF). In this study, the samples were subjected to the same conditions to obtain an X-ray energy intensity distribution of calcium (Ca) and strontium (Sr) with a pixel size of 0.5 µm and a dwell time of 40 ms. The acquisition time ranged from 1 to 3 h depending on the size of the surface area.
Sr:Ca analysis between samples was deduced from Sr and Ca X-ray signal intensity ratio. In order to retrieve the actual Sr mass fraction in the CaCO3 matrix in each sample, a calibration process was performed using PyMca stand-alone open-source software, institutionally supported and updated by the European Synchrotron Radiation Facility (ESRF) for spectral analysis (PyMca 5.9.2 http://pymca.sourceforge.net/, accessed on 6 June 2024) [38]. The corrected mass fraction was then translated to a molar fraction of Sr:Ca.

2.3. Best Use of Synchrotron Facility Allocation Time

The access to synchrotron beam time is not commonplace, hence most of the time allotted at the synchrotron facility site was dedicated to beam time use for data acquisition. Alongside the incompressible scan time, the large amount of raw data generated is imported reduced, corrected and pre-processed within the time frame of the experiment through an in-house built workflow pipeline for on-site data viewing. Then, a data reconstruction process makes the data outputs portable through regular storage devices thanks to a user-friendly interface and in the ready-to-use format (.txt, .TIF). All of the further data treatment steps are then feasible by the researchers on standard working stations.

2.4. Sr:Ca Ratio 2D Raster

Ca and Sr XRF data in TIF format was loaded as images using Fiji (Fiji Is Just ImageJ) software, the collaborative open source with a self-update project for scientific image analysis (http://fiji.sc/cgi-bin/gitweb.cgi, accessed on 6 June 2024) [39]. Each 2D raster image is displayed scaled in grey values, by assigning the black colour to the pixels with the highest fluorescence intensity signal, white for the pixel with the lowest intensity, and scaling the in-between intensities accordingly. Nevertheless, the encoded pixel value still calls the actual energy intensity figure in any of the calculations made by image processing, provided the high dynamic range (32 bit) is maintained. Sr:Ca 2D rasters were obtained through image processing using the ‘divide function’ of the ‘image calculator’ package. The tool performs the arithmetic operations pixel by pixel between the two ‘images’.
The outcome is a 2D raster of Sr:Ca X-ray energy ratio map, which is then either tuned to pixel-by-pixel colour allocation or used as is for their value dataset. In each raster, the Sr:Ca values along with their spatialization by addressing their cartesian coordinates were exported and saved as .csv (comma-separated values) in order to analyse zonation and intensity frequency distribution.

2.5. Colour Map Tuning

‘Red-to-Slate Grey’ Sr:Ca Ratio Permutation

Sr:Ca colour maps are displayed by assigning a colour to the pixel’s value following a home-designed ‘red-yellow-brown-cyan-slate grey’ palette for permutation (scale range provided by the inset in the figures). An in-house built lookup table (LUT) was set using the LUT tool of the Fiji package library. The rational colour permutation criteria were normalised using the range of Sr:Ca ratios for all samples. To do so, a normalization step was performed by setting manually a few square pixels in each image background to the maximum Sr:Ca ratio value retrieved in the whole sampling. Then, the in-house lookup table was applied to the image scaled in grey values. The colour maps obtained reflected the differences in intensity matching the Sr:Ca ratio range in the sample and provided a visualization of the phenotypic variations between the samples.

2.6. 3-3-2 RGB Image Conversion

The 3-3-2 RGB conversion process was applied to the Sr:Ca 2D rasters using the dedicated LUT tool of the Fiji software. Tabulation transforms the pixel value by addressing three times the value to red colour, three times to green and two times to blue. This mathematical transformation allows visualization of neighbouring information. The closer the adjoining Sr:Ca values are, the smoother the offset in the converted image is. This method provides an interesting addition to the basic Sr:Ca colour scale mapping as it can unmask the presence of subtle local gradients in the range of the raw values or, conversely, highlight sustained granular zones.

2.7. Post-Treatment Computation for Sr:Ca Ratio Frequency Distribution

As described in a previous study [7], zonal changes in Sr XRF signal intensity depict three concentric zones in the otoliths. This configuration is used to allocate each raster coordinate as belonging to one of three zones, which are named here—from the core to the periphery—Zone 1 (Z1), Zone 2 (Z2) and Zone 3 (Z3).
This was carried out by using Fiji software to adjust boundaries on the signal output intensity transitions and select the discrete areas of which the corresponding pixel cartesian coordinates (x,z) were allocated by means of specific transformation in the result (.csv) file and export process.
Then, the Sr:Ca values retrieved from the 2D raster were subjected to zone allocation by matrix superimposition of the values to the cartesian coordinate sets using basic tools of the collaborative open source R project available on the Comprehensive R Archive Network (https://cran.r-project.org/, accessed on 6 June 2024). All analyses were carried out with R version 4.2.3. Data tagged with allocation was then gathered and processed using ‘dplyr’ package and ‘summarise’ function in order to aggregate and summarise the Sr:Ca value frequencies. The figures were set up by running the ggplot2 package geom_area function for normalized frequencies (reported to surface area) above the 1% threshold value.

3. Results

For each otolith analysed, 2D maps show that Sr content is distributed in a concentric manner. The visualization of the Sr:Ca ratio reveals an integration of Sr within the otolith matrix in every growth direction.

3.1. Output Process

Figure 2 sketches the key steps of the data analysis. Prior to the analysis, we used the following starting postulates: (i) Sr is a major mineral component of the otolith and an element that marks water salinity, (ii) otoliths record time and chemical information and (iii) synchrotron-based 2D XRF is a semi-quantitative method. Figure 2a shows a zonal pattern on the Sr distribution image when Ca’s was smooth and uniform [7] indicating that the spatial distribution of the Sr:Ca ratio was strikingly driven by the Sr variation. The transect drawn from edge to edge through the core (Figure 2a) shows steep changes in the Sr:Ca ratio during the fish’s life (Figure 2b). Three main zones are highlighted by the use of Sr:Ca thresholding (Figure 2c): there is an outermost zone of lowest Sr:Ca intensity (called Z3); Z3 shows a clear break with Z2, of higher Sr:Ca intensity. Finally, the last zone identified, Z1, is a circular zone of c. 20 to 30 μm in diameter with the core at its centre (Figure 2d). Variations in Sr:Ca in the so-called zones were compared between samples by imaging (Figure 2e) and computation, allowing the interpretation of the results both in terms of life cycle hypothesis and Sr incorporation dynamics in the otolith.

3.2. Otolith Sr:Ca Distribution of Discrete Zones

Figure 3 represents the 2D Sr:Ca (mmol·mol−1) maps with the visualisation of the three concentric zones for each otolith. As observed with the 2D maps, over the zonal Sr:Ca pattern, the distribution of Sr:Ca ratios seems to be heterogeneous in the otolith. Figure 4 represents the frequencies of the Sr:Ca ratio in each zone defined in the otolith in order to characterize the variation in the elemental incorporation. The latter representation allows seeing the spread of values in each zone, and apart from a few individuals (see further in the results), Z3 is always clearly separated from the other two zones with homogeneously low Sr:Ca concentrations, consistent with a freshwater environment (around 2 mmol·mol−1) and a low spread of value centred on a peak. Z1 and Z2 always overlap, but Z1 shows a low spread of values centred on a peak unlike Z2, which shows a heterogeneous distribution, a very diffuse spread of values and no peak.

3.3. Environment Transition Signature in M. brachyurus and M. nicoleae

All specimens show a concentric Sr:Ca pattern consistent with transitions between environments of varying salinity (Figure 3a and Figure 4a; Figure 3b and Figure 4b). There is always a clear break between the outermost Z3 and Z2, with a very low Sr:Ca in Z3, consistent with life in freshwater (around 1 to 2 mmol·mol−1) and a Z2 area with Sr:Ca main maximum concentration around 5 to 9 mmol·mol−1, consistent with the marine environment. For most specimens, [Sr:Ca]Z1 < [Sr:Ca]Z2 >> [Sr:Ca]Z3, with a discrete innermost Z1 region (e.g., JAP 8094). This pattern will be called ‘classic pattern’ hereafter. Based on this classic pattern there are a number of noticeable variations: there are specimens for which [Sr:Ca]Z1 = [Sr:Ca]Z2 >> [Sr:Ca]Z3 (e.g., NC 8099 or POL 3912); several specimens (e.g., 8091, 3917 and 3922) show many areas in Z2 with a low salinity signature, but with the occurrence of patches of high concentrations (>6 mmol·mol−1) giving a marine signature, suggesting that these individuals with a ‘patchy’ Z2 migrated between environments of different salinities. The results suggest that these two species from all locations (Japan, PNG, Solomon Islands, New Caledonia and French Polynesia) are migrating species. The four M. brachyurus specimens from French Polynesia show a very heterogeneous Z2 with patches of very high concentration (up to about 16 mmol·mol−1; see specimen 3915), which do not concur with a concentric growth pattern. These specimens show the highest concentrations of Sr:Ca in Z2, compared to fish caught in other regions (maximum 7 mmol·mol−1 for specimens caught in Japan, 9 mmol·mol−1 in PNG, 12 mmol·mol−1 in the Solomon Islands and around 8 mmol·mol−1 in New Caledonia).

3.4. M. torrentius

For specimen 8073, Z1 signs for an environment under the influence of high salinity with Sr:Ca up to 6 mmol·mol−1. Z2 mode is around 3 mmol·mol−1 (freshwater) and Z3 is less than 1 mmol·mol−1. This specimen does not seem to migrate to the sea during the juvenile phase represented by Z2. The other M. torrentius specimens (8084, 8077) show a ‘patchy’ Z2 with only some areas of high Sr:Ca concentration (between 6 and 12 mmol·mol−1), testifying to a passage in the marine environment (Figure 3c and Figure 4c).

3.5. M. cruentus

For this endemic species to New Caledonia (Figure 3d and Figure 4d), the three zones Z1, Z2 and Z3 are clearly visible, but the variability of Sr:Ca between the different areas in the otolith is less marked, with overlapping values for the three zones delimited in the otolith indicating that this species probably stays in water bodies of low salinity throughout its life. For all specimens except 8068, Sr:Ca concentrations remain low and representative of freshwater residency (<4 mmol·mol−1). Although heterogeneously distributed in Z2, 8068 shows an overall high Sr:Ca (mean 6 mmol·mol−1 for Z2), meaning that this fish may have passed through a marine environment during the juvenile growth phase.

3.6. Added Value of 3-3-2 RGB Image Conversion Method

The visualisation of Sr:Ca mathematically transformed on interleaved blocks of 3_square_pixels allows disentangling focal organisation in patches and local gradients and provides valuable additional information to the Sr:Ca maps (Figure 5). For all samples, Z2 shows a granular pattern but there are ring patterns at the Z2/Z3 boundary, which marks the return to freshwater after the marine phase. Interestingly, the ring zonation often extends in Z2. We can invoke that the ring patterns revealed by this method indicate stabling in estuaries or wandering at the river mouth. A likely explanation for the granular pattern could be the result of physiological and growth rate parameters overriding environmental chemistry at early times in the life cycle.

4. Discussion

In the present work, we analysed the Sr:Ca ratio of the otoliths of four species of freshwater pipefish. Diadromy was our entry point to study Sr incorporation in the otolith using 2D XRF mapping. Strontium (Sr) mapping revealed three alternating concentric zones consistent with the fish’s life stages (Figure 2, Figure 3, Figure 4 and Figure 5). This Sr integration is compatible with a circumperipheric continuous growth pattern as Sr is integrated in all growth directions. This Sr integration has already been shown for the European hake and the European sea bass: Sr:Ca does not show proximo-distal differences [26].

4.1. M. brachyurus and M. nicoleae: Amphidromous Species

The ‘classic pattern’ was the most common among M. brachyurus and M. nicoleae. It strongly suggests that these individuals are amphidromous: Z1 corresponds to hatching postulated to take place in a freshwater environment; Z2 with higher Sr:Ca corresponds to a marine phase (between 6 and 10 mmol·mol−1 and even higher) and Z3 with very low Sr concentration corresponds to life back in freshwater (<2 mmol·mol−1). These three life phases correspond to a typical amphidromous life cycle [5,7]. Adults live in freshwater; during field missions, adult fish are caught using an electric portative fishing device, which only functions in freshwater. Pipefish live upriver, as far up as the first impassable barrier, so often several hundred meters or several kilometres from the river mouth, far from the influence of salt water [40]. Mature adults are regularly observed in the rivers with the eggs borne by the adult male (pers. obs). The eggs hatch in freshwater and the juvenile phase, the phase where growth is very strong, takes place at sea. Juveniles then return to rivers where they colonise the adult habitat. Haÿ et al. [7] actually validated the amphidromous life cycle for the two species cited here and gave each zone a specific name: Z1 corresponds to hatching in freshwater (hFW), Z2 corresponds to the marine phase in sea water (SW) and Z3 corresponds to the adult phase in freshwater (AdFW). However, the Z1 region often exhibited Sr concentrations typical of environments of higher salinity than expected in freshwater (>5 mmol·mol−1). This could be explained by the fact that adults sometimes migrate to the estuary for the eggs to hatch closer to the sea, making it easier for juveniles to reach the marine environment.
M. brachyurus is the most widespread species studied in this paper, as is known from Sri Lanka to French Polynesia. Figure 3a and Figure 4a show that all the specimens studied are amphidromous, which is consistent with this species’ wide distribution. Marine Z2 signatures are quite different between specimens caught in the West Pacific Ocean (Japan, PNG, Solomon Islands and New Caledonia) and French Polynesia. These differences may be due to exposure to different environmental conditions (sea water signature may be different around French Polynesia than around Japan or New Caledonia) and/or to variations in the individual integration pathways. These different signatures may be used in the future to study migration routes for these amphidromous species.

4.2. New Caledonian Microphis: A Versatile Lifestyle?

M. torrentius is known from New Caledonia and Vanuatu, two close archipelagos, and Samoa further to the East. The three specimens studied in this paper were caught in New Caledonia. The maps and Sr:Ca profiles showed results comparable to what was analysed for M. brachyurus and M. nicoleae with two specimens (Figure 3c and Figure 4c: 8077 and 8084) with quite low Sr:Ca Z1 signature, a Z2 with patches of higher Sr:Ca signing for a marine environment and a Z3 of very low Sr:Ca, signing for freshwater. These two specimens are thought to be amphidromous. Specimen 8073, however, does not seem to migrate to sea during the juvenile growth phase. Indeed, Z1 shows a higher Sr:Ca (up to 6 mmol·mol−1) showing that the individual was under the influence of higher salinity at birth, probably as the adult migrated to the estuary during reproduction. Z2 then has a low Sr:Ca, signing for freshwater residency (around 2–3 mmol·mol−1) during the juvenile growth phase, and Z3 is typical of a freshwater environment (<1 mmol·mol−1). This specimen, 8073, stayed all its life in the river and did not undertake its juvenile development at sea. It did not have an amphidromous life cycle.
For M. cruentus (Figure 3d) the three concentric zones, Z1, Z2 and Z3, driven by Sr incorporation are clearly visible with Z1 and Z2 being higher in Sr concentration than Z3. M. cruentus 8068 has a ‘patchy’ Z2 with only some areas of high Sr:Ca concentration (up to 8 mmol·mol−1) (Figure 3d and Figure 4d), testifying to a passage in a marine environment during juvenile growth. For the other three M. cruentus specimens, the Sr:Ca molar ratio in all three zones reflects that of a freshwater environment (<4 mmol·mol−1) or at least a less pronounced transition to an environment of higher salinity. Maps shown in Figure 3d (NC 8066, 8071 and 8069) and graphs in Figure 4d suggest that these individuals are freshwater residents and that juvenile development took place in the river.
These results suggest that some specimens of M. torrentius and M. cruentus do not seem to carry out their marine migration while others go to sea. Both species have a versatile lifestyle with a facultative amphidromous life cycle. Facultative amphidromy has been observed in other species, usually endemic or with a restricted geographic distribution. New Caledonia is known to harbour many endemic species, both plants and animals and both terrestrial and aquatic [41]. This is partially explained by the fact that New Caledonia is a particular geological formation with a metamorphic unit in the northeastern part of the island and peridotite nappes with an ultramafic nickel-rich soil in the southern part of the island [16]. Although the amphidromous life cycle is clear for the freshwater goby endemic to New Caledonia Sicyopterus sarasini Weber & Beaufort, 1915 [5], it is not obvious for the New Caledonian endemic goby Protogobius attiti Watson and Pöllabauer, 1988 [42]. The preference for an ultramafic substrate by New Caledonian species may explain the loss of an amphidromous life cycle. This observation is also relevant for the endemic M. cruentus for which facultative amphidromy is observed here. The number of specimens should be increased in order to state more precisely the amphidromous versatile modalities. Facultative amphidromy is known for fish and has been already observed in Hawaiian or Asian endemic species [43,44,45]. According to Closs et al. [46], the absence of complete migration between freshwater and marine environments in a few populations of amphidromous species (e.g., Cottus, Rhinogobius and Galaxias) could be explained by their specialization in the upstream habitats. Freshwater residency shows advantages: reduced mortality due to predation or the loss of larvae in ocean currents as well as reduced energy costs caused by the environmental transition [47]. The fragility of newly hatched pelagic individuals is a major cost of this life cycle, as they are unable to withstand starvation and other adverse physical and chemical conditions likely to be encountered during their downstream migration to reach their pelagic nursery habitat [48]. While some individuals may lose their ability to disperse, amphidromy is not lost on an evolutionary time scale, as it remains advantageous at the population or species levels for the colonisation of new habitats and for genetic mixing [47]. Therefore, variation in the amphidromous life cycle could be an adaption to particular environmental conditions. Since M. torrentius has a broader distribution than the endemic M. cruentus (New Caledonia, Vanuatu and Samoa), it was surprising to uncover the possibility of freshwater residency. Maybe it is caused by the particular conditions in New Caledonia. In any case, additional otoliths of this species from its entire distribution area should be analysed in order to quantify better the importance of facultative amphidromy for M. torrentius.

4.3. Advantages of 2D XRF Mapping When Addressing Biological Questions

The amphidromous life cycle and its versatility could have been uncovered using any 2D technic described in the introduction. 1D transect would even have allowed validating the amphidromous life cycle for otoliths with a ‘classic pattern’. But XRF methods allowed visualizing Sr integration within the calcified lattice at a hyperfine scale, finer than any of the other 2D methods, which led to surprising results. As our method mixes element quantification and 2D imaging, we were able to see that Sr concentration, especially in the Z2 area, is not at all homogeneously distributed. The spread of Sr:Ca values more or less important observed for Z2 (Figure 4) is consistent with the non-homogeneous integration of elements during the growth of the otolith during this phase representing an early life phase. Z2 has, for many specimens a patchy distribution, which does not fit with the otolith’s concentric growth even though hollow ring patterning can be unravelled underneath apparent gross features. This means that within a continuous growth band, heterogeneities occur. The heterogeneity of elemental distribution within the otolith matrix is greater than what is generally assumed [1,10,25,29]. Sr integration in the otolith may depend on the occurrence of vaterite inclusions, one of the three natural polymorphs of CaCO3 (calcite, aragonite and vaterite). Vaterite was shown to trap less Sr than aragonite [25,49,50]. In the present case, the four species studied were previously verified for the occurrence of CaCO3 polymorphs in their otoliths by using Raman spectroscopy. All the otoliths analysed only exhibited aragonite spectra (Supplementary Material Figure S1). As a result, we can hypothesize that this unusual pattern of Sr distribution is due to biomineralisation processes and the underlying control mechanisms. Most studies support that Sr is incorporated in equilibrium with environmental concentration [1,51,52,53], but that patchy distribution shows that Sr incorporation is also controlled by other additional factors, e.g., temperature, ontogeny, growth, metabolism, and this may particularly be true during the juvenile stage. Most studies find no effect of the temperature on Sr incorporation [53,54], but some experimental studies show that temperatures close to the species’ thermal tolerance limits may affect element uptake, transport in blood plasma or biomineralisation processes [55,56]. Temperatures in the rivers usually are around 23 to 25 °C (pers. obs). Sea surface temperatures around PNG, the Solomon Islands and French Polynesia are around 30 °C (https://www.ospo.noaa.gov, accessed on 22 February 2024). For amphidromous species, migrating from freshwater at hatching to the sea to spend their juvenile phase, the increase in water temperature of 5 to 8 °C may have a strong impact on individual metabolism and otolith growth via increased accretion [57]. In terms of endogenous and physiologically driven processes, ontogeny has a great impact on otolith growth. Although the environmental component in Sr uptake and incorporation is strong, it has been shown that ontogenetic changes alter fish growth and ion processing; therefore, important physiological and morphological changes in the fish’s life can affect Sr otolith incorporation, in addition to the strong environmental component [58,59]. There are very few studies on Sr incorporation independent of habitat use but some note that Sr can exhibit an ontogenetic decrease during larval or juvenile stages [60,61] while other studies show an increase in Sr with fish age [62,63]. For Microphis species, regardless of habitat use, Sr incorporation may be higher during the early life and juvenile stages due to physiological and genetic controls [64]. That may explain why even non-amphidromous individuals have higher Sr concentrations in Z1 and Z2 than in Z3.
Amphidromous species have to adapt twice to environmental changes affecting osmoregulation. The passage to marine life as juveniles is physiologically and metabolically challenging [65] and it impacts element incorporation within the otolith matrix, which may be reflected by the patchy pattern.
The rules of physico-chemistry may also be challenged by biological control [66]. For instance, the granular pattern could reflect the biological or genetic factors acting at a certain period in the life cycle on the growth directions of the otoliths and controlling their shape. The 3-3-2 RGB transformation allowed the enhancement of Sr incorporation information and visualisation of the internal growth dynamic of the mineral (Figure 4): the core zone Z1 and the first growth phase of Z2 have a granular aspect, which does not reflect the otolith growth by concentric accretion, and Sr incorporation is more discreet (mathematically speaking). Concentricity appears mid-Z2 or at the end of Z2, with a rhythmicity of Sr incorporation, particularly in the transition between Z2 and Z3. This transition is likely to reflect the transition from the sea to freshwater, and to a greater extent, the transition from a fast-growing juvenile phase to a slower growing adult phase (since this phenomenon is also seen on non amphidromous individuals). This shows that Sr incorporation during the early juvenile phase is under a biological driver and is influenced by ontogeny and growth.
The patchy effect can be very minimalistic in some specimens as patches of higher Sr concentration can only appear either on only one side of the otolith (towards the sulcus or towards the distal area) or on very small areas of Z2 and was uncovered only because hyperfine-scale XRF mapping was used (see Figure 3a M. brachyurus 8083; Figure 3b—M. nicoleae PNG 3922, 3917, SOL 3919, 3920; Figure 3c—M. torrentius NC 8077, NC 8084; Figure 3d—M. cruentus NC 8068). For these individuals, adding a second dimension to the analysis has proven to be a highly effective method. Typically, 1D analysis could have led to misinterpretation of the biology of these species as a transect, or points, chosen in a patchless area would have resulted in concluding to freshwater residency, while patches of high Sr are proof of environmental transition through a marine phase and show that the individual has undergone amphidromous migrations. We then conclude that all specimens analysed are amphidromous apart from three M. cruentus specimens, which remain their entire life in freshwater and M. torrentius specimens for which passage in marine environment is questionable.

5. Conclusions

Even though otolith microchemical composition is mainly driven by environmental elemental composition, biotic factors clearly act on element incorporation in the otolith lattice [1]. Synchrotron-based scanning X-ray fluorescence 2D high-resolution mapping allowed us to tackle global and hyperfine strontium heterogeneous distribution. It led to the discovery of plasticity in the freshwater pipefish life cycle as well as the trade-off expression between biological and chemical driving forces that govern element incorporation. This study illustrates the value of otolith chemistry to track fish movements, but caution must be taken against the routine interpretation of fish movement from otolith Sr:Ca variation under the assumption of a strict correlation with salinity. The analytical imaging process we developed and our results provide a strong basis for future works and prospects in fish biology and ecology. Fundamental knowledge of plasticity in the life cycle of amphidromous species has important implications in terms of biodiversity surveys, as varying life history traits are the motor of speciation and evolution. Furthermore, understanding life cycle modalities is crucial to categorize species distribution limits and to implement adapted conservation measures, especially when endemic species are at stake.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology13060446/s1, Figure S1: Otolith Raman analysis.

Author Contributions

Conceptualization, C.L., S.B. and P.K.; Methodology, S.B. and K.M.; Investigation, V.H. and S.B.; writing—original draft preparation, C.L. and S.B.; writing—review and editing, C.L. and S.B. and P.K.; supervision, P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available at https://data.indores.fr/dataset.xhtml?persistentId=doi:10.48579/PRO/2REAXP, accessed on 6 June 2024.

Acknowledgments

We acknowledge SOLEIL for the provision of synchrotron radiation facilities and we would like to thank K. Medjoubi and A. Somogyi for assistance in using beamline “NANOSCOPIUM” and monitoring Proposal 20210141. The authors would like to acknowledge assistance from the colleagues and support groups of Soleil for their contributions, continuous help, and outstanding technical support. Data were obtained on the Nanoscopium beamline. For these special thanks go to the beamline crew K. Medjoubi and A. Somogyi for their advice and constructive suggestions. For French Polynesia, we would like to thank V. Mazel (Ichtyo-Pacific). For PNG, thanks to G. Kaipu (PNG NRI) for the research permit (N°018196), N. Gowep (CEPA) for the export permit and our friends, P. Amick, B. Ruli (Live & Learn), J. Anamiato (National Museum and art Gallery of PNG) and D. Vaghelo (WNB Provincial Govt, Environment section). We want to thank all the responsible Chiefs of the areas investigated for their kind permission, and the village communities who have always heartily received us and helped us in our prospecting of rivers. For the Solomon Islands, we would like to thank D. Boseto (ESSI) for his invaluable help on the field, the landowners and tribes of the Solomon Islands for allowing the expedition team to enter their customary lands, and the Government for the support and facilitation of the legal process. For New Caledonia, we would like to thank the New Caledonian Government and the New Caledonian North and South Provinces, the DAFE (C. Fort), J.J. Cassan, and N. Charpin from the NGO “Vies d’Ô Douce” for his help on the field. We also would like to thank K. Maeda and H. Kobayashi for specimens from Japan. The authors gratefully acknowledge Oulfa Belhadj, head of vibrational spectroscopy equipment set (Raman, FTIR) for her scientific advice and her technical assistance, as well as CRC MNHN laboratory for Raman beamline time.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hüssy, K.; Limburg, K.E.; De Pontual, H.; Thomas, O.R.; Cook, P.K.; Heimbrand, Y.; Blass, M.; Sturrock, A.M. Trace Element Patterns in Otoliths: The Role of Biomineralization. Rev. Fish. Sci. Aqua. 2021, 29, 445–477. [Google Scholar] [CrossRef]
  2. Campana, S.E.; Thorrold, S.R. Otoliths, Increments, and Elements: Keys to a Comprehensive Understanding of Fish Populations? Can. J. Fish. Aquat. Sci. 2001, 58, 30–38. [Google Scholar] [CrossRef]
  3. Teichert, N.; Valade, P.; Grondin, H.; Trichet, E.; Sardenne, F.; Gaudin, P. Pelagic Larval Traits of the Amphidromous Goby Sicyopterus lagocephalus Display Seasonal Variations Related to Temperature in La Réunion Island. Ecol. Freshw. Fish. 2016, 25, 234–247. [Google Scholar] [CrossRef]
  4. Taddese, F.; Reid, M.; Heim-Ballew, H.; Jarvis, M.G.; Closs, G.P. Otolith Chemistry of Triplefin Forsterygion nigripenne Indicates Estuarine Residency. Fish. Sci. 2021, 87, 271–281. [Google Scholar] [CrossRef]
  5. Lord, C.; Tabouret, H.; Claverie, F.; Pécheyran, C.; Keith, P. Femtosecond Laser Ablation ICP-MS Measurement of Otolith Sr:Ca and Ba:Ca Composition Reveal Differential use of Freshwater Habitats for Three Amphidromous Sicyopterus (Teleostei: Gobioidei: Sicydiinae) species. J. Fish Biol. 2011, 79, 1304–1321. [Google Scholar] [CrossRef] [PubMed]
  6. Reis-Santos, P.; Gillanders, B.M.; Sturrock, A.M.; Izzo, C.; Oxman, D.S.; Lueders-Dumont, J.A.; Walther, B.D. Reading the Biomineralized Book of Life: Expanding Otolith Biogeochemical Research and Applications for Fisheries and Ecosystem-Based Management. Rev. Fish. Biol. Fisheries. 2023, 33, 411–449. [Google Scholar] [CrossRef]
  7. Haÿ, V.; Berland, S.; Medjoubi, K.; Somogyi, A.; Mennesson, M.I.; Keith, P.; Lord, C. Unmasking Pipefish Otolith Using Synchrotron-Based Scanning X-ray Fluorescence. Sci. Rep. 2023, 13, 4794. [Google Scholar] [CrossRef] [PubMed]
  8. Secor, D.; Rooker, J. Is Otolith Strontium a Useful Scalar of Life Cycles in Estuarine Fishes? Fish. Res. 2000, 46, 359–371. [Google Scholar] [CrossRef]
  9. Limburg, K.E.; Huang, R.; Bilderback, D.H. Fish Otolith Trace Element Maps: New Approaches with Synchrotron Microbeam X-ray Fluorescence. X-ray Spectrom. J. 2007, 36, 336–342. [Google Scholar] [CrossRef]
  10. Zimmerman, C.E. Relationship of otolith strontium-to-calcium ratios and salinity: Experimental validation for juvenile salmonids. Can. J. Fish. Aquat. Sci. 2005, 62, 88–97. [Google Scholar] [CrossRef]
  11. Zimmerman, C.E.; Swanson, H.K.; Volk, E.C.; Kent, A.J. Species and Life History Affect the Utility of Otolith Chemical Composition for Determining Natal Stream of Origin for Pacific Salmon. Trans. Am. Fish. Soc. 2013, 142, 1370–1380. [Google Scholar] [CrossRef]
  12. Feutry, P.; Castelin, M.; Ovenden, J.R.; Dettaï, A.; Robinet, T.; Cruaud, C.; Keith, P. Evolution of Diadromy in Fish: Insights from a Tropical Genus (Kuhlia Species). Am. Nat. 2013, 181, 52–63. [Google Scholar] [CrossRef]
  13. Arai, T. Early Life History and Recruitment Processes of a Tropical Anguillid Eel Anguilla marmorata to the Pacific Coast, as Revealed by Otolith Sr:Ca Ratios and Microstructure. Biology 2022, 11, 803. [Google Scholar] [CrossRef] [PubMed]
  14. Taillebois, L.; Tabouret, H.; Pecheyran, C.; Keith, P. Inputs from microchemistry to the understanding of three Sicydiinae species’ life cycle. Life. Environ. 2015, 65, 73–84. [Google Scholar]
  15. McDowall, R.M. On Amphidromy, a Distinct Form of Diadromy in Aquatic Organisms. Fish. Fish. 2007, 8, 1–13. [Google Scholar] [CrossRef]
  16. Keith, P.; Lord, C. Tropical Freshwater Gobies: Amphidromy as a Life Cycle. In The Biology of Gobies; Patzner, R.A., Van Tassel, J.L., Kapoor, B.G., Eds.; Science Publishers: Jersey, UK; Enfield, NH, USA, 2011; pp. 243–277. [Google Scholar]
  17. Keith, P. Biology and Ecology of Amphidromous Gobiidae of the Indo-Pacific and the Caribbean Regions. J. Fish. Biol. 2003, 63, 831–847. [Google Scholar] [CrossRef]
  18. Van der Laan, R.; Eschmeyer, W.N.; Fricke, R. Family-Group Names of Recent Fishes. Zootaxa 2014, 3882, 1–230. [Google Scholar] [CrossRef] [PubMed]
  19. Haÿ, V.; Mennesson, M.I.; Lord, C.; Keith, P. Why Several When One Can Unite Them All? Integrative Taxonomic Revision of Indo-Pacific Freshwater Pipefish (Nerophinae). J. Linn. Soc. London Zool. 2023, 198, 923–956. [Google Scholar] [CrossRef]
  20. Vincent, A.C.J.; Foster, S.J.; Koldewey, H.J. Conservation and Management of Seahorses and Other Syngnathidae. J. Fish. Biol. 2011, 78, 1681–1724. [Google Scholar] [CrossRef]
  21. Gu, L.; Zhang, H.; Feng, G.; Liu, Y.; Han, Z.; Zhao, F.; Ye, Q.; Hu, W.; Song, C. Evaluation of the Effectiveness of Marking Juvenile Takifugu obscurus Otoliths With Strontium. Fishes 2022, 7, 371. [Google Scholar] [CrossRef]
  22. Ferreira, I.; Daros, F.A.; Moreira, C.; Feijó, D.; Rocha, A.; Mendez-Vicente, A.; Pisonero, J.; Correia, A.T. Is Chelidonichthys lucerne (Linnaeus, 1758) a Marine Estuarine-Dependent Fish? Insights from Saccular Otolith Microchemistry. Fishes 2023, 8, 383. [Google Scholar] [CrossRef]
  23. Kerr, L.A.; Campana, S.E. Chemical Composition of Fish Hard Parts as a Natural Marker of Fish Stocks. In Stock Identification Methods, 2nd ed.; Cadrin, S.X., Kerr, L.A., Mariani, S., Eds.; Elsevier: Boston, MA, USA, 2014; pp. 205–234. [Google Scholar]
  24. Walther, B.D.; Limburg, K.E. The Use of Otolith Chemistry to Characterize Diadromous Migrations. J. Fish. Biol. 2012, 81, 796–825. [Google Scholar] [CrossRef] [PubMed]
  25. Tzeng, W.N.; Chang, C.W.; Wang, C.H.; Shiao, J.C.; Iizuka, Y.; Yang, Y.J.; You, C.F.; Lozys, L. Misidentification of the Migratory History of Anguillid Eels by Sr/Ca Ratios of Vaterite Otoliths. Mar. Ecol. Prog. Ser. 2007, 348, 285–295. [Google Scholar] [CrossRef]
  26. de Pontual, H.; MacKenzie, K.M.; Tabouret, H.; Daverat, F.; Mahé, K.; Pecheyran, C.; Hüssy, K. Heterogeneity of Otolith Chemical Composition from Two-Dimensional Mapping: Relationship With Biomineralization Mechanisms and Implications for Microchemistry Analyses. J. Fish. Biol. 2024, 104, 20–33. [Google Scholar] [CrossRef] [PubMed]
  27. Thorrold, S.R.; Shuttleworth, S. In Situ Analysis of Trace Elements and Isotope Ratios in Fish Otoliths Using Laser Ablation Sector Field Inductively Coupled Plasma Mass Spectrometry. Can. J. Fish. Aquat. Sci. 2000, 57, 232–242. [Google Scholar] [CrossRef]
  28. FitzGerald, J.L.; Thorrold, S.R.; Bailey, K.M.; Brown, A.L.; Severin, K.P. Elemental Signatures in Otoliths of Larval Walleye Pollock (Theragra chalcogramma) From the Northeast Pacific Ocean. Fish. Bull. 2004, 102, 604–616. [Google Scholar]
  29. Limburg, K.E.; Elfman, M. Insights From Two-Dimensional Mapping of Otolith Chemistry. J. Fish. Biol. 2017, 90, 480–491. [Google Scholar] [CrossRef] [PubMed]
  30. Dimasi, E.; Sarikaya, M. Synchrotron X-ray Microbeam Diffraction From Abalone Shell. J. Mater. Res. 2004, 19, 1471–1476. [Google Scholar] [CrossRef]
  31. Bottini, C.; Dapiaggi, M.; Erba, E.; Faucher, G.; Rotiroti, N. High Resolution Spatial Analyses of Trace Elements in Coccoliths Reveal New Insights into Element Incorporation in Coccolithophore Calcite. Sci. Rep. 2020, 10, 9825. [Google Scholar] [CrossRef]
  32. Cuif, J.P.; Dauphin, Y.; Lo, C.; Medjoubi, K.; Saulnier, D.; Somogyi, A. Synchrotron-Based HR-Fluorescence and Mineralogical Mapping of the Initial Growth Stages of Polynesian Cultivated Pearls Disprove the ‘Reversed Shell’ Concept. Minerals 2022, 12, 172. [Google Scholar] [CrossRef]
  33. Adams, F.; Janssens, K.; Snigirev, A. Microscopic X-ray Fluorescence Analysis and Related Methods with Laboratory and Synchrotron Sources. J. Anal. At. Spectrom. 1998, 13, 319–331. [Google Scholar] [CrossRef]
  34. Hedges, K.J.; Ludsin, S.A.; Frayer, B.J. Effects of Ethanol Preservation on Otolith Microchemistry. J. Fish Biol. 2004, 64, 923–937. [Google Scholar] [CrossRef]
  35. Somogyi, A.; Medjoubi, K.; Baranton, G.; Le Roux, V.; Ribbens, M.; Polack, F.; Philippot, P.; Samama, J.P. Optical Design and Multi-Scale Scanning Spectro-Microscopy Possibilities at the Nanoscopium Beamline of Synchrotron Soleil. J. Synchrotron Radiat. 2015, 22, 1118–1129. [Google Scholar] [CrossRef] [PubMed]
  36. Bergamaschi, A.; Medjoubi, K.; Messaoudi, C.; Marco, S.; Somogyi, A. MMX-I: Data-Processing Software for Multimodal X-Ray Imaging and Tomography. J. Synchrotron Rad. 2016, 23, 783–794. [Google Scholar] [CrossRef] [PubMed]
  37. Van Grieken, R.; Markowicz, A. (Eds.) Handbook of X-ray Spectrometry; CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
  38. Solé, V.A.; Papillon, E.; Cotte, M.; Walter, P.; Susini, J.A.A. Multiplatform Code for the Analysis of Energy-Dispersive X-Ray Fluorescence Spectra. Spectroc. Acta Part B At. Spectrosc. 2007, 62, 63–68. [Google Scholar] [CrossRef]
  39. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An Open-Source Platform for Biological-Image Analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [PubMed]
  40. Keith, P.; Boseto, D.; Lord, C. Freshwater Fish of the Solomon Islands; Société Française d’Ichtyologie: Paris, France, 2021; pp. 1–173. [Google Scholar]
  41. Caesar, M.; Grandcolas, P.; Pellens, R. Outstanding micro-endemism in New Caledonia: More Than One Out of Ten Animal Species Have a Very Restricted Distribution Range. PLoS ONE 2017, 12, e0181437. [Google Scholar] [CrossRef] [PubMed]
  42. Tabouret, H.; Tomadin, M.; Taillebois, L.; Iida, M.; Lord, C.; Pécheyran, C.; Keith, P. Amphidromy and Marine Larval Phase of Ancestral Gobioids Rhyacichthys guilberti and Protogobius attiti (Teleostei: Rhyacichthyidae). Mar. Freshw. Res. 2014, 65, 776–783. [Google Scholar] [CrossRef]
  43. Heim-Ballew, H.; Moody, K.N.; Blum, M.J.; McIntyre, P.B.; Hogan, J.D. Migratory Flexibility in Native Hawai’ian Amphidromous fishes. J. Fish Biol. 2020, 96, 456–468. [Google Scholar] [CrossRef]
  44. Liao, T.Y.; Huang, W.C.; Iizuka, Y.; Chou, M.T.; Shiao, J.C. Facultative Amphidromy and Pelagic Larval Duration Plasticity of Rhinogobius formosanus (Teleostei, Gobioidei). ZooKeys 2020, 951, 91. [Google Scholar] [CrossRef]
  45. Lisi, P.J.; Hogan, J.D.; Holt, G.; Moody, K.N.; Wren, J.L.; Kobayashi, D.R.; Blum, M.J.; McIntyre, P.B. Stream and Ocean Hydrodynamics Mediate Partial Migration Strategies in an amphidromous Hawaiian goby. Ecology 2022, 103, e3800. [Google Scholar] [CrossRef] [PubMed]
  46. Closs, G.P.; Hicks, A.S.; Jellyman, P.G. Life Histories of Closely Related Amphidromous and Non-Migratory Fish Species: A Trade-off Between Egg Size and Fecundity. Freshw. Biol. 2013, 58, 1162–1177. [Google Scholar] [CrossRef]
  47. Hogan, J.D.; Blum, M.J.; Gilliam, J.F.; Bickford, N.; McIntyre, P.B. Consequences of Alternative Dispersal Strategies in a Putatively Amphidromous Fish. Ecology 2014, 95, 2397–2408. [Google Scholar] [CrossRef]
  48. Iida, M.; Watanabe, S.; Tsukamoto, K. Oceanic Larval Duration and Recruitment Mechanism of the Amphidromous Fish Sicyopterus japonicus (Gobioidei: Sicydiinae). Reg. Stud. Mar. Sci. 2015, 1, 25–33. [Google Scholar] [CrossRef]
  49. Brown, R.J.; Severin, K.P. Elemental Distribution within Polymorphic Inconnu (Stenodus leucichthys) Otoliths is Affected by Crystal Structure. Can. J. Fish. Aquat. Sci. 1999, 56, 1898–1903. [Google Scholar] [CrossRef]
  50. Melancon, S.; Fryer, B.J.; Ludsin, S.A.; Gagnon, J.E.; Yang, Z. Effects of Crystal Structure on the Uptake of Metals by Lake Trout (Salvelinus namaycush) Otoliths. Can. J. Fish. Aquat. Sci. 2005, 62, 2609–2619. [Google Scholar] [CrossRef]
  51. Elsdon, T.S.; Gillanders, B.M. Relationship Between Water and Otolith Elemental Concentrations in Juvenile Black Bream Acanthopagrus butcheri. Mar. Ecol. Prog. Ser. 2003, 260, 263–272. [Google Scholar] [CrossRef]
  52. Elsdon, T.S.; Gillanders, B.M. Consistency of Patterns Between Laboratory Experiments and Field Collected Fish in Otolith Chemistry: An Example and Applications for Salinity Reconstructions. Mar. Freshw. Res. 2005, 56, 609–617. [Google Scholar] [CrossRef]
  53. Reis-Santos, P.; Tanner, S.E.; Elsdon, T.S.; Cabral, H.N.; Gillanders, B.M. Effects of Temperature, Salinity and Water Composition on Otolith Elemental Incorporation of Dicentrarchus labrax. J. Exp. Mar. Biol. Ecol. 2013, 446, 245–252. [Google Scholar] [CrossRef]
  54. Walther, B.D.; Kingsford, M.J.; O’Callaghan, M.D.; McCulloch, M.T. Interactive Effects of Ontogeny, Food Ration and Temperature on Elemental Incorporation in Otoliths of a Coral Reef Fish. Environ. Biol. Fish 2010, 89, 441–451. [Google Scholar] [CrossRef]
  55. Thomas, O.R.B.; Ganio, K.; Roberts, B.R.; Swearer, S.E. Trace Element–Protein Interactions in Endolymph From the Inner Ear of Fish: Implications for Environmental Reconstructions Using Fish Otolith Chemistry. Metallomics 2017, 9, 239–249. [Google Scholar] [CrossRef] [PubMed]
  56. Izzo, C.; Reis-Santos, P.; Gillanders, B.M. Otolith Chemistry Does not Just Reflect Environmental Conditions: A Meta-Analytic Evaluation. Fish Fish. 2018, 19, 441–454. [Google Scholar] [CrossRef]
  57. Fablet, R.; Pecquerie, L.; de Pontual, H.; Høie, H.; Millner, R.; Mosegaard, H.; Kooijman, S.A. Shedding Light on Fish Otolith Biomineralization Using a Bioenergetic Approach. PLoS ONE 2011, 6, e27055. [Google Scholar] [CrossRef] [PubMed]
  58. Jessop, B.; Cairns, D.; Thibault, I.; Tzeng, W. Life History of American Eel Anguilla rostrata: New Insights From Otolith Microchemistry. Aquat. Biol. 2008, 1, 205–216. [Google Scholar] [CrossRef]
  59. Brown, R.J.; Severin, K.P. Otolith Chemistry Analyses Indicate That Water Sr:Ca is the Primary Factor Influencing Otolith Sr:Ca for Freshwater And Diadromous Fish But Not for Marine Fish. Can. J. Fish. Aquat. Sci. 2009, 66, 1790–1808. [Google Scholar] [CrossRef]
  60. de Pontual, H.; Lagardere, F.; Amara, R.; Bohn, M.; Ogor, A. Influence of Ontogenetic and Environmental Changes in the Otolith Microchemistry of Juvenile Sole (Solea solea). J. Sea Res. 2003, 50, 199–211. [Google Scholar] [CrossRef]
  61. Morales-Nin, B.; Swan, S.C.; Gordon, J.D.M.; Palmer, M.; Geffen, A.J.; Shimmield, T.; Sawyer, T. Age-Related Trends in Otolith Chemistry of Merluccius merluccius from the North-Eastern Atlantic Ocean and the Western Mediterranean Sea. Mar. Freshw. Res. 2005, 56, 599–607. [Google Scholar] [CrossRef]
  62. Hughes, J.M.; Stewart, J.; Gillanders, B.M.; Collins, D.; Suthers, I.M. Relationship Between Otolith Chemistry and Age in a Widespread Pelagic Teleost Arripis trutta: Influence of Adult Movements on Stock Structure and Implications for Management. Mar. Freshw. Res. 2016, 67, 224–237. [Google Scholar] [CrossRef]
  63. Grammer, G.L.; Morrongiello, J.R.; Izzo, C.; Hawthorne, P.J.; Middleton, J.F.; Gillanders, B.M. Coupling Biogeochemical Tracers with Fish Growth Reveals Physiological and Environmental Controls on Otolith Chemistry. Ecol. Monogr. 2017, 87, 487–507. [Google Scholar] [CrossRef]
  64. Söllner, C.; Burghammer, M.; Busch-Nentwich, E.; Berger, J.; Schwarz, H.; Riekel, C.; Nicolson, T. Control of Crystal Size and Lattice Formation by Starmaker In Otolith Biomineralization. Science 2003, 302, 282–286. [Google Scholar] [CrossRef]
  65. Taillebois, L.; Keith, P.; Valade, P.; Torres, P.; Baloche, S.; Dufour, S.; Rousseau, K. Involvement of Thyroid Hormones in the Control of Larval Metamorphosis in Sicyopterus lagocephalus (Teleostei: Gobioidei) at the Time of River Recruitment. Gen. Comp. Endocrinol. 2011, 173, 281–288. [Google Scholar] [CrossRef] [PubMed]
  66. Otter, L.M.; Agbaje, O.; Kilburn, M.R.; Lenz, C.; Henry, H.; Trimby, P.; Hoppe, P.; Jacob, D.E. Insights into Architecture, Growth Dynamics, and Biomineralization From Pulsed Sr-Labelled Katelysia rhytiphora Shells (Mollusca, Bivalvia). Biogeosciences 2019, 16, 3439–3455. [Google Scholar] [CrossRef]
Biology 13 00446 g001
Figure 1. Collection sites in the Western and Central Pacific Ocean region (red squares). .br—Microphis brachyurus, .ni—Microphis nicoleae, .to—Microphis torrentius, .cr Microphis cruentus.
Figure 1. Collection sites in the Western and Central Pacific Ocean region (red squares). .br—Microphis brachyurus, .ni—Microphis nicoleae, .to—Microphis torrentius, .cr Microphis cruentus.
Biology 13 00446 g001
Biology 13 00446 g002
Figure 2. Workflow stages for the analysis of 2D Sr:Ca distribution. (a) Zonal pattern is highlighted by the XRF signal in the otolith section down to the core (red arrow); (b) Transition between zones are shown by steep slopes on the graph (red dotted lines) of Sr:Ca XRF signal ratio along a transect (orange transect) through the core (red arrow). (c) Three main zones can be defined by thresholding the Sr:Ca signal ratio values: Z1, quite circular zone of c. 20 to 30 μm centred around the core, then Z2 and Z3 from the core to the periphery, respectively. (d) Each pixel on the otolith section is assigned a tag addressing the zone it belongs to; it will then be used to extract the spatialized data. (e) Phenotypic variations between samples were shown using a colour gradient for the Sr:Ca ratio.
Figure 2. Workflow stages for the analysis of 2D Sr:Ca distribution. (a) Zonal pattern is highlighted by the XRF signal in the otolith section down to the core (red arrow); (b) Transition between zones are shown by steep slopes on the graph (red dotted lines) of Sr:Ca XRF signal ratio along a transect (orange transect) through the core (red arrow). (c) Three main zones can be defined by thresholding the Sr:Ca signal ratio values: Z1, quite circular zone of c. 20 to 30 μm centred around the core, then Z2 and Z3 from the core to the periphery, respectively. (d) Each pixel on the otolith section is assigned a tag addressing the zone it belongs to; it will then be used to extract the spatialized data. (e) Phenotypic variations between samples were shown using a colour gradient for the Sr:Ca ratio.
Biology 13 00446 g002
Biology 13 00446 g003aBiology 13 00446 g003b
Figure 3. (a–d) Fine and global Sr:Ca ratio intensities highlighted in colour scale on the 2D rasters. (a). M. brachyurus—(b). M. nicoleae—(c). M. torrentius—(d). M. cruentus. (PNG: Papua New Guinea). Colour scale: Sr:Ca mmol·mol−1.
Figure 3. (a–d) Fine and global Sr:Ca ratio intensities highlighted in colour scale on the 2D rasters. (a). M. brachyurus—(b). M. nicoleae—(c). M. torrentius—(d). M. cruentus. (PNG: Papua New Guinea). Colour scale: Sr:Ca mmol·mol−1.
Biology 13 00446 g003aBiology 13 00446 g003b
Biology 13 00446 g004aBiology 13 00446 g004b
Figure 4. (a–d) Frequency distribution of the Sr:Ca molar ratio values computed from each individual 2D raster after the data has been assigned to the spatial zones of the otolith, Z1 (light blue), Z2 (orange) and Z3 (slate grey). The horizontal axis plots the Sr:Ca in molar concentration converted from XRF energy ratio values. Frequency (y axis) was normalized to the depicted zone surface (pixel unit). JAP: Japan, POL: French Polynesia, PNG: Papua New Guinea, SOL: Solomon Islands, NC: New Caledonia. (a). .br: M. brachyurus—(b). .ni. M. nicoleae—(c). .to M. torrentius—(d). .cr: M. cruentus.
Figure 4. (a–d) Frequency distribution of the Sr:Ca molar ratio values computed from each individual 2D raster after the data has been assigned to the spatial zones of the otolith, Z1 (light blue), Z2 (orange) and Z3 (slate grey). The horizontal axis plots the Sr:Ca in molar concentration converted from XRF energy ratio values. Frequency (y axis) was normalized to the depicted zone surface (pixel unit). JAP: Japan, POL: French Polynesia, PNG: Papua New Guinea, SOL: Solomon Islands, NC: New Caledonia. (a). .br: M. brachyurus—(b). .ni. M. nicoleae—(c). .to M. torrentius—(d). .cr: M. cruentus.
Biology 13 00446 g004aBiology 13 00446 g004b
Biology 13 00446 g005
Figure 5. Colour-scaled Sr:Ca signal ratio raster pictures (left side) and their corresponding 3-3-2-RGB conversion (right side). The 3-3-2-RGB conversion method allows for enhanced offset information of each pixel with its immediate neighbours. br: M. brachyurus, ni: M. nicoleae, to: M. torrentius, cr: M. cruentus—JAP: Japan, POL: French Polynesia, PNG: Papua New Guinea, SOL: Solomon Islands, NC: New Caledonia—Scale bar: 50 µm—3-3-2-RGB conversion images were processed with filter for easier colour blindness reading.
Figure 5. Colour-scaled Sr:Ca signal ratio raster pictures (left side) and their corresponding 3-3-2-RGB conversion (right side). The 3-3-2-RGB conversion method allows for enhanced offset information of each pixel with its immediate neighbours. br: M. brachyurus, ni: M. nicoleae, to: M. torrentius, cr: M. cruentus—JAP: Japan, POL: French Polynesia, PNG: Papua New Guinea, SOL: Solomon Islands, NC: New Caledonia—Scale bar: 50 µm—3-3-2-RGB conversion images were processed with filter for easier colour blindness reading.
Biology 13 00446 g005
Table 1. Sampling: species name, ID: field number; SL: Standard length in mm; collection location details (country, island, river); XRF otolith scan number; Museum national d’Histoire naturelle (MNHN, Paris, France) collection number.
Table 1. Sampling: species name, ID: field number; SL: Standard length in mm; collection location details (country, island, river); XRF otolith scan number; Museum national d’Histoire naturelle (MNHN, Paris, France) collection number.
SpeciesIDSLSexCountryIslandRiverScanCollection Number
M. brachyurusJP48175.75FJapanOkinawaGenka8067MNHN-IC-2023-0276
M. brachyurusJP44187.81FJapanOkinawaGenka8079MNHN-IC-2021-0320
M. brachyurusJP35173.66MJapanOkinawaTeima8080MNHN-IC-2023-0274
M. brachyurusJP32165.12FJapanOkinawaTeima8091MNHN-IC-2023-0275
M. brachyurusJP33180.69MJapanOkinawaTeima8094MNHN-IC-2023-0275
M. brachyurus17765129.38MPapua New GuineaNew BritainRangihi3928MNHN-IC-2021-0318
M. brachyurus17766122.49FPapua New GuineaNew BritainRangihi3943MNHN-IC-2021-0318
M. brachyurus19190101.98FSolomon IslandsKolombangaraVaqe3916MNHN-IC-2021-0331
M. brachyurusRTNC033124.32FNew CaledoniaGrande TerreCreek Baie Nord8072MNHN-IC-2023-0273
M. brachyurusRTNC117168.60MNew CaledoniaGrande TerrePwé hiit8075MNHN-IC-2023-0272
M. brachyurusRTNC145120.41FNew CaledoniaGrande TerreGarana8081MNHN-IC-2023-0277
M. brachyurusRTNC139105.57FNew CaledoniaGrande TerreGarana8083MNHN-IC-2021-0322
M. brachyurusRTNC109131.88MNew CaledoniaGrande TerrePwé hiit8099MNHN-IC-2023-0272
M. brachyurusPFV52125.32FFrench PolynesiaTahitiPapenoo3912MNHN-IC-2023-0050
M. brachyurusPFV06105.02FFrench PolynesiaTahitiPapenoo3913MNHN-IC-2023-0051
M. brachyurusPFV39109.36MFrench PolynesiaTahitiPapenoo3915MNHN-IC-2023-0052
M. brachyurusPFV42104.15FFrench PolynesiaTahitiPapenoo3933MNHN-IC-2023-0053
M. nicoleae1905582.26MPapua New GuineaNew BritainWalindi3917MNHN-IC-2023-0045
M. nicoleae1918369.41FPapua New GuineaNew BritainGavuvu3939MNHN-IC-2023-0446
M. nicoleae1917690.86MPapua New GuineaNew BritainGavuvu3940MNHN-IC-2023-0446
M. nicoleae1769393.66MPapua New GuineaNew BritainHoskins3922MNHN-IC-2021-0338
M. nicoleae1496264.37FSolomon IslandsRanonggaOvana3918MNHN-IC-2021-0336
M. nicoleae1496183.33FSolomon IslandsRanonggaOvana3919MNHN-IC-2021-0036
M. nicoleae1825397.65MSolomon IslandsSanta IsabelRakata3920MNHN-IC-2023-0047
M. nicoleae18268103.90MSolomon IslandsSanta IsabelRakata3924MNHN-IC-2023-0048
M. torrentiusRTNC11984.82FNew CaledoniaGrande TerreGarana8073MNHN-IC-2021-0335
M. torrentiusRTNC11889.69MNew CaledoniaGrande TerreGarana8077MNHN-IC-2021-0335
M. torrentiusRTNC114C86.15MNew CaledoniaGrande TerrePwé hiit8084MNHN-IC-2021-0334
M. cruentusRTNC080112FNew CaledoniaGrande TerreCarénage8068MNHN-IC-2021-0348
M. cruentusRTNC107149.66MNew CaledoniaGrande TerreOuenghi8066MNHN-IC-2021-0349
M. cruentusRTNC106110.23FNew CaledoniaGrande TerreOuenghi8069MNHN-IC-2021-0349
M. cruentusRTNC108168.60FNew CaledoniaGrande TerreOuenghi8071MNHN-IC-2021-0349
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lord, C.; Haÿ, V.; Medjoubi, K.; Berland, S.; Keith, P. Travelling in Microphis (Teleostei: Syngnathidae) Otoliths with Two-Dimensional X-ray Fluorescence Maps: Twists and Turns on the Road to Strontium Incorporation. Biology 2024, 13, 446. https://doi.org/10.3390/biology13060446

AMA Style

Lord C, Haÿ V, Medjoubi K, Berland S, Keith P. Travelling in Microphis (Teleostei: Syngnathidae) Otoliths with Two-Dimensional X-ray Fluorescence Maps: Twists and Turns on the Road to Strontium Incorporation. Biology. 2024; 13(6):446. https://doi.org/10.3390/biology13060446

Chicago/Turabian Style

Lord, Clara, Vincent Haÿ, Kadda Medjoubi, Sophie Berland, and Philippe Keith. 2024. "Travelling in Microphis (Teleostei: Syngnathidae) Otoliths with Two-Dimensional X-ray Fluorescence Maps: Twists and Turns on the Road to Strontium Incorporation" Biology 13, no. 6: 446. https://doi.org/10.3390/biology13060446

APA Style

Lord, C., Haÿ, V., Medjoubi, K., Berland, S., & Keith, P. (2024). Travelling in Microphis (Teleostei: Syngnathidae) Otoliths with Two-Dimensional X-ray Fluorescence Maps: Twists and Turns on the Road to Strontium Incorporation. Biology, 13(6), 446. https://doi.org/10.3390/biology13060446

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop