The brown sea cucumber Isostichopus fuscus (Ludwig, 1875) is an echinoderm of commercial importance with an elongated body and a soft, tough and thick border with irregular papules that inhabits rock bottoms, coral reefs and is occasionally found on sandy or soft bottoms, down to 61 m depth. This species is found from the Gulf of California to the Galapagos Islands and the coast of Ecuador (Deichmann, 1985; Maluf, 1988; Sonnenholzner et al., 2013; Toral-Granda & Martínez, 2007). Its spatial distribution depends on substrate complexity, depth, temperature, salinity, light, rugosity, and food availability (Maluf, 1988).

In Mexico, I. fuscus was intensively fished from 1988 to 1993, with catches that surpasses 1 000 tones shortly after, its abundance was so low that it was listed as an at-risk species by the Mexican government (NOM-059-SEMARNAT-1994) and its fishing was completely banned (SEMARNAT, 1994). In 2002, the ban was suspended but fishing is dependent upon stock assessments (SEMARNAT, 2002) and currently the species is listed under special protection (SEMARNAT, 2010). This situation has motivated studies on the reproduction and growth (Herrero-Pérezrul, Reyes-Bonilla, García-Domínguez, & Cintra-Buenrostro, 1999), population dynamics (Reyes-Bonilla & Herrero-Pérezrul, 2003) fishery (Herrero-Pérezrul & Chávez, 2005), weight-length relationships and condition index (Herrero-Pérezrul & Reyes-Bonilla, 2008) of the brown sea cucumber. However, all of these studies were conducted in the southern portion of the Baja California peninsula.

There are a number of studies of the reproduction and growth of the brown sea cucumber that suggest that sea surface temperature, light, dissolved oxygen concentration, pH, and the concentration of ammonia, nitrites and inorganic phosphorous are variables that influence the optimal development of this species (Asha & Muthiah, 2005). Five gonad stages have been described for this species: I) undifferentiated, II) gametogenesis, III) maturity, IV) spawning and V) post-spawning (Fajardo-León, Michel-Guerrero, Singh-Cabanillas, Vélez-Barajas, & Massó-Rojas, 1995). It has also been observed that spawning coincides with rising temperatures and increased light (Fajardo-León et al., 1995; Herrero-Pérezrul et al., 1999; Toral-Granda & Martínez, 2007). Reproduction follows a seasonal pattern in Baja California Sur, México (Herrero-Pérezrul et al., 1999). However, there are no reproductive studies of this species at the northernmost limits of its distribution. Therefore, the objective of this study was to describe the reproductive cycle of I. fuscus and to evaluate the potential effects of environmental conditions on the cycle. We hypothesized that rising temperature and increased chlorophyll a concentration, a proxy for primary production, would promote gonad development.

MATERIALS AND METHODS

Study Area: The gulf of California (GC) has been divided in three faunistic regions (Brusca & Brusca, 2002). The northern region is influenced by the Colorado river delta, with high salinity and sea surface temperatures (SST) ranging between 10°C in winter and up to 32°C in summer. The salinity and temperature distributions observed are driven by seasonal flows of heat and humidity, with strong tidal and convective mixing (Paden, Abbott, & Winant, 1993; Soto-Mardones, Marinone, & Parés-Sierra, 1999). The central region or Greater Islands region (Isla Tiburon and Angel de la Guarda) presents lower SST due to intense tidal mixing, with strong seasonal changes from 16°C in winter to 31°C in summer (Álvarez-Borrego, 2007). The southern region or mouth of the gulf, is characterized by high primary productivity due to high nutrient concentrations that come from wind-induced upwelling, tidal mixing and the exchange of water masses between the gulf and the Pacific Ocean (Álvarez-Borrego, 2007; Lavín & Marinone, 2003).

Field work: Sampling surveys were conducted monthly from October 2014 to December 2016 in 118 sites located in the northern and central regions of the GC. For logistic reasons the area was divided in three zones (Fig. 1). At each site SST, bottom temperature, salinity, dissolved oxygen and salinity were recorded with a sonde (YSI 85). Collection was performed by night because sea cucumbers hide during daylight (Reyes-Bonilla, Ramírez-Ortiz, Herrero-Pérezrul, & Calderón-Aguilera, 2016) at three depths: shallow (1-9 m), mid-water (10-17 m) and deep (18-27 m). All specimens found along 2 x 25 m transects (two transects per depth, i.e., six transects per site) were collected and taken on board. Once aboard all specimens were left resting and then weighed (OHAUSMR; precision 0.1 g) and sized with a flexible tape (1 mm) and returned on site. Only 30 specimens were kept for histological analysis in each monthly survey. The gonad was extracted and preserved in 10% formalin for three days and then transferred to 70% alcohol for histology.

Lab work: A portion of gonad was fixed in formaldehyde solution (10%) for approximately 24 h, dehydrated in a graded alcohol series, cleared with xylene and embedded in paraplast at 56oC. Histological sections (4 μm) were cut with a Leica 2040 Autocut microtome and stained with hematoxylin and eosin, following Humason (1962) and Bancroft, Cook, & Stirling (1988). Gametogenic development and sex ratios were determined by examination of histological samples at different magnifications (4X, 10X, 20X and 40X; Microscope Nikon Optiphot-II, coupled with a digital camera Sight DS-5M and the digital system Digital Sight DS-L1; Table 2 supplemental material).

Data analysis: Satellite data of monthly mean of chlorophyll a (Chl a) concentration and sea surface temperature (SST) with a 4 x 4 km pixel resolution from October 2014 to December 2016 were gathered from AQUAMODIS (Ocean Color Web http://oceancolor.gsfc.nasa.gov/cms/) and processed with SeaDAS (ver. 7.3.1). A time series from 2003-2016 of Chl a and SST were generated using NASA’s Ocean Color Radiometry Online Visualization and Analysis Tool utility (https://giovanni.sci.gsfc.nasa.gov/giovanni/). These series were standardized (µ = 0, SD = 1) to compare them with those observed during the study period. A cross-correlation analysis was performed to determine whether the timing of gametogenic development and spawning was correlated with Chl a concentration and SST using STATISTICA (Statsoft ver. 7.1).

The weight - length relationship was described with a power model of the form W = a∙Lb, where W is weight, L is length, a is the intercept and b is the allometry coefficient (Safran, 1992). The growth parameters were estimated through the Von Bertalanffy equation: Lt = L∞ (1 - [e –k (t - t0)]), where Lt is the size or mean length at age t and L∞ is the maximum theoretical length that an individual can reach, k is the growth rate and t0 is an initial fitting parameter. This was done using the ELEFAN routine (Electronic Frequency Analysis) from FiSAT (ver. 1.2.2) (© FAO 2006-2016), following Gayanilo & Pauly (1997). A χ2 test was performed to test if the sex ratio was 1:1. The oocyte area was estimated using ImageJ 1.46 (http://imagej.net/ImageJ). A total of 10 291 oocytes that were contained within a field and presented a well-defined germinal vesicle were measured. The oocytes theoretical diameter (OTD) was calculated through the equation OTD= √4A/π, where OTD is in μm, A is the oocyte area and π is a constant (Saout, Quéré, Donval, Paulet, & Samain, 1999). In order to identify if there was more than one cohort, a modal progression analysis of OTD (Hmida, Ayache, Haouas, & Romdhane, 2010) was performed using ELEFAN from FiSAT II.

RESULTS

A total of 2 808 specimens of Isostichopus fuscus from 25 months of sampling (October 2014-October 2015 and January-December 2016) were evaluated. The mean length observed (± SD) was 21.4 ± 6.0 cm (min: 5 cm, max: 41 cm; Table 1 Supplemental Material). From the length frequency distribution depicted in Figure 2, the minimum size (5 cm) was recorded in October 2014 and March 2015 and the maximum (41 cm) in December 2014. The mean weight was 375.6 ± 249.1 g (min: 4 g, recorded in May 2015, max: 1 200 g, recorded in June 2015 and in June, July and October 2015; Table 1 supplemental material) and the frequency distribution is shown in Figure 3. The best represented weight class was 301-400 g; on average, males were bigger than females, but the difference was not statistically significant for either length (t = -0.98, P > 0.05, df = 381) or weight (t = -1.79, P > 0.05, df = 379).

Five hundred and fifty-three gonads were obtained from 717 dissected organisms, from which 224 were females, 162 were males, 157 were undifferentiated and 10 were hermaphrodites (Fig. 4). The sex ratio was statistically different from 1:1 (χ2 = 36.63, P = 0.03, df = 23). Sexually active specimens accounted for 77% of the sample with a mean length of 27.5±4.06 cm (min: 16, max: 41 cm) and weight of 647 ± 224 g (min: 50 g, max: 1,200 g).

The weight-length relationship is described as W = 0.18L2.4 (r2 = 0.82). The Von Bertalanffy parameters are Lt = 39.9 cm and k = 0.18, supposing a t0 = 0. With those parameters, the size-at-first maturity is estimated to be between 16 and 18 cm.

Histological characteristics: Five stages of gonadal development were identified: undifferentiated, gametogenesis (oogenesis and spermatogenesis), maturity, spawning (spawning and expulsion), post-spawning (post-spawning and post-expulsion, Table 2 supplementary material, Fig. 6). For both sexes, the gonad consists of beams of numerous tubules of different lengths of white to beige coloration.

Undifferentiated

From all examined gonads, 28% were undifferentiated and they were present throughout the year but more frequently from December to March (Fig. 5). At this stage, it is impossible to distinguish sex, because there are no gametes, only abundant conjunctive tissue and a thick gonadal wall. In some cases it was possible to observe phagocytes; this stage may be present in juvenile specimens or specimens that have passed the post-spawning stage, as the gonad is reabsorbed.

Gametogenesis

Only 9% of the sample was found at this stage. Specimens in the gametogenesis stage were observed from February to October and more frequently from April to August (Fig. 5). At this stage there is a proliferation of female and male gametes.

Oogenesis: In the gonad wall, oogonia are embedded or attached to the germinal epithelium, proliferating towards the lumen. Primary oocytes have a defined nucleus and peripheral nucleoli. In some cases, small cells (two or three products of meiosis) are seen around the primary oocyte. A large amount of connective tissue is observed, as well as folds in the wall of the gonad, which are reduced as gametogenesis progresses (Fig. 6 a).

Spermatogenesis: In spermatogenesis, the gonadal wall has longitudinal folds that extend towards the lumen and on the periphery of the wall there are spermatogonia of spherical appearance. There is proliferation of primary spermatocytes of ovoid and lumpy aspect and spermatocytes are seen arranged in folds. Mature spermatocytes can be observed in the lumen. The connective tissue begins to decrease in thickness by the accumulation of spermatocytes (Fig. 6 b).

Maturity

Fifteen percent of the collected sample were classified as mature. Mature specimens were observed from June to November, with a peak in the months of July to October (Fig. 5). At this stage, the gametes fill most of the light in the follicles. In females, the mature oocyte is completely round, measuring approximately 70 μm. The nucleus of the oocyte is easily distinguishable and it is possible to observe the three peripheral nuclei, as well as developing oocytes. The gonadal wall is thin and no connective tissue is visible. Although uncommon, the presence of atresia in the lumen can be observed, surrounded by phagocytes. Phagocytes are observed inside and outside the lumen (Fig. 6 c). In males, the follicles are completely saturated with mature spermatozoa, arranged in dense layers. Mature spermatozoa have a spherical shape, flattened slightly in the dorsoventral direction. The gonadal wall is thin and completely deployed, on the periphery there are spermatocytes and spermatids (Fig. 6 d).

Spawning

Nineteen percent of the sample was in the spawning stage, often in the months of August to October (Fig. 5). In this stage, a decrease of the gametes inside the follicles is observed, as well as empty spaces in the lumen.

Spawning: In this stage, little folds are observed in the gonadal wall, as well as a thickening of the gonadal wall and increasing connective tissue. The residual oocytes are in maturity or oogenesis. A large quantity of phagocytes are observed inside and outside the lumen. In some cases it is possible to observe the gametes in atresia (Fig. 6 e).

Expulsion: A separation between the mature spermatocytes at the center of the lumen and the immature spermatocytes at the periphery of the follicles is observed. The gonadal wall remains thin, although a thin layer of connective tissue is observed. It is possible to observe phagocytes (Fig. 6 f).

Post-spawning

Twenty nine percent of the sample was in the post-spawning stage, mostly present in November and December (Fig. 5). At this stage, the lumen of the follicles is practically empty. The gonadal wall is thick and has a large amount of connective tissue. A large quantity of phagocytes are observed inside and outside of follicles. The gonad is practically reabsorbed by phagocytic activity.

Post-spawning: The oocytes present are amorphous and of smaller size, clearly distinguishable atresia in different sizes (Fig. 6 g).

Post-expulsion: Distinction between spermatozoa and phagocytes is difficult, since the remaining spermatozoa are being phagocytosed (Fig. 6 h).

Hermaphrodite

Two percent (n = 10) of the sexed specimens were hermaphroditic (Fig. 7). Nine of the ten hermaphroditic organisms were collected in San Luis Gonzaga in the months of September and October of 2015 (Fig. 5). Female and male follicles were found in different gonadal stages, and three specimens were in spawning stage. In four specimens, female follicles were spawning and male follicles were in gametogenesis and in two specimens, female follicles were mature and male follicles were in gametogenesis. The other hermaphroditic specimen was collected in San Francisquito in July 2016; in this specimen both female and male follicles were spawning.

Gametogenic stage and environmental variables: For the study period, Chl a SST anomalies were observed and compared to historical values. Even when SST followed a seasonal pattern (min 17°C in winter and 31°C max in summer), there was a large difference in October 2015 of 6.5°C higher than the historical average (Fig. 8). Chl a follows an irregular pattern, with peaks in February, March, April and November, reaching a maximum of 5.96 mg Chl a m-3 in February 2015, with means of 4.3 mg Chl a m-3 more than the historical average for that month. However, in 2016 there were negative anomalies in January, February, March, May, June, July, September and October compared to the historical average.

The reproductive peak occurs on June, coinciding with an increase in SST. The cross-correlation analysis shows a positive correlation between gametogenesis and SST two months before the observed reproductive peak (r = 0.42, P < 0.05). On the other hand, there is negative correlation between Chl a and gametogenesis (r = -0.65, P < 0.05) four months before gametogenic development. Gamete expulsion occurs from August to October, when SST is higher and this was corroborated by the cross-correlation analysis (r = 0.57, P < 0.05). Moreover, there was a negative correlation between Chl a and gamete expulsion (r = -0.71, P < 0.05) when lag = 0.

Oocyte theoretical diameter (OTD) and modal progression analysis (MPA): A total of 10 291 oocytes were measured (mean ± SD = 65.3 ± 19 μm). The mean OTD during gametogenesis was 66.0 ± 17 μm), 73.8 ± 14 μm during maturity, 60.4 ± 20 μm during spawning and 49.2 ± 20 μm during post-spawning. The MPA shows that the OTD increases 23%, i.e., 15 μm∙month-1 when passing from oogenesis to maturity, and decreases 9% (7 μm∙month-1) when passing from maturity to spawning; finally, OTD in post-spawning is, on average, 53 μm∙month-1 smaller than during spawning.

DISCUSSION

In a seminal paper on echinoderm reproduction, Mercier & Hamel (2009) affirm that every time there is more evidence that there are exogenous factors that affect the reproduction of these animals. This is the first work conducted at the northernmost limit of the distribution of Isostichopus fuscus, as well as the first to address the relationship between its reproductive cycle and environmental variables. This is of utmost importance considering that a population living at the limit of its distribution may be more sensitive to environmental variables and that may manifest as anomalies in gonadal development.

On average, sea cucumbers from the north of the gulf of California are smaller (21.4 cm) and lighter (375.6 g) than those from the central and southern regions of the gulf (27° a 29 °N; 25.1-32 cm and 458-562 g: Fajardo-León et al., 1995; Fajardo-León & Vélez-Barajas, 1996 and 23.3cm and 385.9g: Herrero-Pérezrul & Reyes-Bonilla, 2008) and Oaxaca, México (16° N; 396 g, Glockner-Fagetti & Benítez-Villalobos, 2017) but larger than in the Galápagos (0° N; 20 cm and 271 g: Sonnenholzner, 1997 and 20.8 cm: Toral-Granda & Martínez, 2007). Even when sizes may indicate a gradient from north to south, as proposed by Herrero-Pérezrul et al., (1999), it might also be the result of manipulation, as holothurids are very malleable. Conand, (1990) mentions that aspidoquirots commercial species such as Holothuria scabra, H. nobilis and H. fuscogilva have sizes larger than 28 cm and may weigh over 900 g. Size is very important for fisheries and commercial purposes as 90% of body weight is lost during dehydration (Herrero-Pérezrul & Reyes-Bonilla, 2008).

The total weight frequency distribution of I. fuscus was found to be between 4-1 200g, with specimens most frequent falling between 200-400g. Toral-Granda & Martínez (2007) consider juveniles to be those individuals below 150 g. Following that criteria, 19% of the sampled population in this study was juvenile. Moreover, Nuño-Hermosillo (2003) considered recruits to be specimens shorter than 13 cm. Therefore, recruitment at our study region takes place between October and February, when 7% of the population was between 4-12 cm of length.

The smallest reproductively active organism of I. fuscus was 16 cm long; however, 50% of the population was sexually active at ≅ 25 cm. In the south of the Baja California peninsula, size-at-fist-maturity was found at five years or at 21 cm long (Herrero-Pérezrul & Chávez, 2005; Herrero-Pérezrul & Reyes-Bonilla, 2008), which is very close to what was found in this study.

The weight-length parameters a and b have a biological meaning; in holothurids 35% of its weight is soft tissue (guts and gonads) and water; therefore, the weight-length relationship could be affected by manipulation, either by retention or expulsion of water. Our estimates of a (0.18) and b (2.4) are different from those of Fajardo-León et al., (1995; a = 0.60 and b = 1.6) and Herrero-Pérezrul & Reyes-Bonilla (2008; a = 1.6 and b = 1.8) for this species.

The 1:1 sex ratio is common in species with sexual reproduction (Hamel & Mercier, 1996; Pennington, 1985) and this is the case for most holothurids, but not always. In the case of I. fuscus there is contradictory evidence regarding the sex ratio. In Mexico, Herrero-Pérezrul et al., (1999) found a 1:1 sex ratio, whereas in Ecuador, Toral-Granda & Martínez (2007) found more females than males, as we found in this study. In another species of sea cucumber (Parastichopus parvimensis), Tapia Vázquez et al., (1996) observed that sex ratios varied monthly.

In our study area, I. fuscus exhibits an asynchronous reproductive cycle between sexes which is different from what other authors have observed in tropical regions (Conand, 1990; Fajardo-León et al., 1995; Tapia Vázquez et al., 1996). However, this same species may present continuous spawning (Toral-Granda & Martínez, 2007) and those differences are related to environmental conditions. One caveat of our study is that we performed monthly surveys and we may have lost more frequent reproductive pulses. This limitation in reproductive studies has already been pointed out by Mercier & Hamel (2009), but it is very difficult to overcome due to logistic and budgeting concerns. Nevertheless, 25 monthly surveys and the histological analysis of hundreds of gonads ensure that our study is representative of the gonad cycle of this species in the northern Gulf of California.

Gonad development starts in March becomes more active between April and July, when temperatures rise. Spawning takes place from June to November, with a peak during August and September, coinciding with higher temperatures. This relationship for holothurids has been observed by many authors (Drumm & Loneragan, 2005; Guzmán, Guevara, & Hernández, 2003; Hamel, Himmelman, & Dufresne, 1993; Hamel & Mercier, 1995; Muthiga, Kawaka, & Ndirangu, 2009; Omar, Abdel Razek, Abdel Rahman, & El Shimy, 2013; Ramofafia, Battaglene, Bell, & Byrne, 2000; Ramofafia, Byrne, & Battaglene, 2001).

The peak in Chl a concentration takes place in February, which might induce the reproductive process due to increased food availability for the system as a whole. Female gametogenesis starts one month earlier than in males, coinciding with the peak of primary productivity, perhaps because of the energetic cost of egg production. Hamel et al., (1993) observed that phytoplankton increase triggered spawning of Psolus fabricii; in our study, the increase in Chl a concentration triggers gametogenesis and spawning takes place during the warm season while undifferentiated individuals are present in the cold months.

The oocyte theoretical diameter (OTD) for other holothurids ranges between 8.93-119.83 µm (Conand, 1982). Costelloe, (1985) found that OTD for Aslia lefevrei was between 40-220 μm in oogenesis, 250-340 μm when active and spawning and between 20-40 μm during post-spawning. In our study, we found a mean OTD of 66 μm in oogenesis, 73.8 μm when active, 60.4 μm spawning and 49.2 μm post-spawning, which is similar to results by Costelloe, (1985).

In marine organisms, there are many factors such as genetics, contaminants and overfishing that can influence hermaphroditism (Ghiselin, 1969). In the case of holothurids, Herrero-Pérezrul, Reyes-Bonilla, & García-Dominguez, (1998) have reported very few cases: 2 out of 155 organisms for Holothuria atra, 1 out of 113 organisms for Peniagone azorica, 1 out of 30 organisms for P. diaphana and 1 out of 429 organisms for Cherbonniera utriculus. For I. fuscus this is the third case reported and the one with the largest proportion (2%) of hermaphrodites: 10 out of 553 organisms. Previously, Herrero-Pérezrul et al. (1999) found 2 out of 173 and organisms, Nuño-Hermosillo (2003) reported one out of 165 organisms. Interestingly, all hermaphrodites came from San Luis Gonzaga, a location with a mining influence and the presence of heavy metals and pesticides (Servicio Geológico Mexicano, 2003). Even when the region has low human activity, resuspension of sediments, mining and mixing due to upwelling, are known to be a source and a means of transport of heavy metals. Moreover, high concentrations of Mn, Al, Zn, Cd and Cu have been detected in mussels from the northern Gulf of California (Gutiérrez-Galindo, Villaescusa-Celaya, & Arreola-Chimal, 1999). Avise (2011) mentions that heavy metals and metaloids affect the endocrine system, reproduction and other biological process. In our study, the hermaphrodites were present spawning as females and present in gametogenesis as males. It has been observed in dioic species that hermaphroditism increases with overfishing (Borgia & Blick, 1981; Ghiselin, 1969). This last author claims that hermaphroditism may be an adaption to increase the chances of successful fertilization (Ghiselin, 1969). We suggest that mining, chemicals and fishing pressure might be driving hermaphroditism in the area.

On the other hand, asexual reproduction has been reported for holothurids (Dolmatov, 2014). Lawrence and Herrera (2000) mention that fission in holothurids seems to be related with habitat and it is induced environmentally and that is why some species draw on fission as an adaptation to sustain populations where recruitment is unlikely. In a recent work, Sonnenholzner, Searcy-Bernal, & Panchana-Orrala, (2017) achieved asexual reproduction of I. fuscus through transverse fission. This is highly relevant for mariculture and restocking under the current depauperate condition of this species in the region (Glockner-Fagetti, Calderon-Aguilera, & Herrero-Pérezrul, 2016).

One of the open questions of this work refers to the effect of fishing on reproduction, because current density is one order of magnitude less than ten years ago (Glockner-Fagetti et al., 2016). In dioic organisms it is not clear how they synchronize to release gametes to find mature conspecifics. In I. fuscus Shepherd, Martinez, Toral-Granda, & Edgar (2004) suggest that a density of 1.2 ind m-2 is needed in order to expect a 50% fertilization success. In our study area the mean density is 0.07 ind∙m-2, well below that density, so more studies on stock assessment and reproduction should be conducted.

acknowledgmentS

The first author received a scholarship during her Master’s studies (CONACYT CVU 599881). T. Rubilar assisted with the description of gonad stages. B. Rojas performed the histological processing of samples collected in 2016. FJ Ponce assisted with graphic preparations. This work was supported by the Fondo Sectorial de Investigación en Materias, Agrícola, Pecuaria, Acuacultura, Agrobiotecnología y Recursos Fitogenéticos (SAGARPA-CONACYT, grant 2013-2-0225235 to LECA). Comments of two anonymous reviewers improved previous versions of this work. Professional proofreading by A. L. McTavish.

RESUMEN

El pepino de mar Isostichopus fuscus es un recurso de gran valor comercial, por lo cual es altamente explotado. Sin embargo, no hay estudios sobre su ciclo reproductivo en la zona más norteña de su distribución. Con el fin de caracterizar su periodo reproductivo, se realizaron expediciones mensuales (Oct 2014-Dic 2016), en la costa oriental de Baja California. Se visitaron 118 sitios comprendidos en la zona costera e insular desde Bahía San Luis Gonzaga hasta el Paralelo 28 para colectar gónadas. Se colectaron 2 808 especímenes de I. fuscus, la talla media ± desviación estándar fue de 21 ± 6.7 cm y el peso medio fue 375.6 ± 249 g. Se disectaron 717 organismos, de los que se obtuvieron 553 gónadas. Mediante análisis histológico se identificó que 224 eran hembras, 162 machos, 157 indiferenciados y 10 hermafroditas. Se encontraron más hembras que machos (χ2 = 36.63, P = 0.03, gl = 23). No hubo diferencias significativas ni en la longitud ni en el peso entre hembras y machos (P > 0.05), sin embargo, las hembras fueron relativamente más grandes. La relación peso-longitud se describe con la ecuación P = 0.18L2.4, r2 = 0.82; la talla de primera madurez se calculó a los 16 cm. Se identificaron cinco estadios de desarrollo gonádico, 28% de las gónadas se catalogaron como indiferenciadas, 9% en gametogénesis, 15% en madurez, 19% en expulsión y 29% en postexpulsión. Se estimó el diámetro teórico del ovocito (DTO) a partir de la medición del área de 10 291 ovocitos. La media del DTO en ovogénesis fue de 65.3 ± 19.7 µm, (n=2 037), en madurez fue de 66.0 ± 17.8 µm (n=3 630), en desove fue de 73.8 ± 14.6 µm (n=3 756) y en postdesove fue de 49.18 ± 20.7 µm (n=868). El análisis de progresión modal indicó que los ovocitos al pasar de la ovogénesis a la madurez incrementan su DTO un 23%, mientras que de la madurez al desove disminuyen 9%, y en postdesove son en promedio 72% más chicos que en desove.

Palabras clave: reproducción; pepino de mar; Holothuroidea; diámetro de ovocitos.

REFERENCES

Álvarez-Borrego, S. (2007). Oceanografía de la región de las grandes islas. In G. D. Danemann & E. Ezcurra (Eds.), Bahía de los Ángeles: recursos naturales y comunidad: linea base 2007 (pp. 45-65). Instituto Nacional de Ecología.

Asha, P. S., & Muthiah, P. (2005). Effects of temperature, salinity and pH on larval growth, survival and development of the sea cucumber Holothuria spinifera Theel. Aquaculture, 250(3), 823-829.

Avise, J. C. (2011). Hermaphroditism: a primer on the biology, ecology, and evolution of dual sexuality. Columbia University Press.

Bancroft, J. D., Cook, H. C., & Stirling, R. W. (1988). Manual of histological techniques and their diagnostic application. Churchill Livingstone.

Borgia, G., & Blick, J. (1981). Sexual competition and the evolution of hermaphroditism. Journal of Theoretical Biology, 89(3), 523-532.

Brusca, R. C., & Brusca, G. J. (2002). Invertebrates (2 ed). Sinauer Associates Incorporated.

Conand, C. (1982). Reproductive cycle and biometric relations in a population of Actinopyga echinites (Echinodermata: Holothuroidea) from the lagoon of New Caledonia, western tropical Pacific. In International Echinoderm Conference, Belkema, Rotterdam (pp. 437-442).

Conand, C. (1990). The Fishery Resourse of Pacific Island Countries: Part 2. Food & Agriculture Org.

Costelloe, J. (1985). The annual reproductive cycle of the holothurian Aslia lefevrei (Dendrochirota: Echinodermata). Marine Biology, 88(2), 155-165.

Deichmann, E. (1985). The Holothurioidea collected by the Velero III during the years 1932 to 1938. University of Suouthern California Press.

Dolmatov, I. (2014). Asexual reproduction in holothurians. The Scientific World Journal, 2014, 1-13.

Drumm, D. J., & Loneragan, N. R. (2005). Reproductive biology of Holothuria leucospilota in the Cook Islands and the implications of traditional fishing of gonads on the population. New Zealand Journal of Marine and Freshwater Research, 39(1), 141-156.

Fajardo-León, M. del C., Michel-Guerrero, E., Singh-Cabanillas, J., Vélez-Barajas, J. A., & Massó-Rojas, J. A. (1995). Estructura poblacional y ciclo reproductor del pepino de mar Isostichopus fuscus en Santa Rosalía, BCS, México. Ciencia Pesquera, 11, 45-53.

Fajardo-León, M. del C., & Vélez-Barajas, A. (1996). Pesquería del pepino de mar. In M. Casas-Vladez & G. Ponce-Díaz (Eds.), Estudio del potencial pesquero y acuícola de Baja California Sur (Vol. 2, pp. 151-165). La Paz: SEMARNAT/CICIMAR.

© FAO 2006-2016. Fisheries and aquaculture software. FISAT II - FAO-ICLARM Stock Assessment Tool. In: FAO Fisheries and Aquaculture Department [online]. Rome.

Gayanilo, F. C. J., & Pauly, D. (1997). FAO-ICLARM Stock Assessment Tools (FiSAT). Reference Manual., 8, 261.

Ghiselin, M. T. (1969). The evolution of hermaphroditism among animals. The Quarterly Review of Biology, 44(2), 189-208.

Glockner-Fagetti, A., & Benítez-Villalobos, F. (2016). Spatio-temporal variation in density and size structure of the endangered sea cucumber Isostichopus fuscus in Huatulco National Park, Mexico. Marine Ecology, 38(1), e12386. http://doi.org/10.1111/maec.12386

Glockner-Fagetti, A., Calderon-Aguilera, L. E., & Herrero-Pérezrul, M. D. (2016). Density decrease in an exploited population of brown sea cucumber Isostichopus fuscus in a biosphere reserve from the Baja California peninsula, Mexico. Ocean & Coastal Management, 121, 49-59.

Gutiérrez-Galindo, E. A., Villaescusa-Celaya, J. A., & Arreola-Chimal, A. (1999). Bioacumulación de metales en Mejillones de cuatro sitios selectos de la región costera de Baja California. Ciencias Marinas, 25(4), 557-578.

Guzmán, H. M., Guevara, C., & Hernández, I. (2003). Reproductive cycle of two commercial species of sea cucumber (Echinodermata: Holothuroidea) from Caribbean Panama. Marine Biology, 142(2), 271-279.

Hamel, J.-F., Himmelman, J. H., & Dufresne, L. (1993). Gametogenesis and spawning of the sea cucumber Psolus fabricii (Duben and Koren). The Biological Bulletin, 184(2), 125-143.

Hamel, J.-F., & Mercier, A. (1995). Spawning of the sea cucumber Cucumaria frondosa in the St. Lawrence Estuary, eastern Canada. Lawrence Estuary, Eastern Canada. SPC Beche-de-Mer Information Bulletin, 7, 12-18.

Hamel, J.-F., & Mercier, A. (1996). Gonad morphology and gametogenesis of the sea cucumber Cucumaria frondosa. SPC Beche-de-Mer Information Bulletin, 8, 22-33.

Herrero-Pérezrul, M. D., & Chávez, E. A. (2005). Optimum fishing strategies for Isostichopus fuscus (Echinodermata: Holothuroidea) in the Gulf of California, México. Revista de Biología Tropical, 53(3), 357.

Herrero-Pérezrul, M. D., & Reyes-Bonilla, H. (2008). Weight-Length relationship and relative condition of the holothurian Isostichopus fuscus at Espíritu Santo Island, Gulf of California, México. Revista De Biologia Tropical, 56(3), 273-280.

Herrero-Pérezrul, M. D., Reyes-Bonilla, H., & García-Dominguez, F. (1998). Casual hermaphroditism in gonochoric Isostichopus fuscus (Ludwig, 1875)(Echinodermata: Holothuroidea) of the southern Gulf of California, Mexico. Bulletin of Marine Science, 63(3), 611-615.

Herrero-Pérezrul, M. D., Reyes-Bonilla, H., García-Domínguez, F., & Cintra-Buenrostro, C. E. (1999). Reproduction and growth of Isostichopus fuscus (Echinodermata: Holothuroidea) in the southern Gulf of California, México. Marine Biology, 135(3), 521-532.

Hmida, L., Ayache, N., Haouas, Z., & Romdhane, M. S. (2010). Oocyte cohort analysis: criteria for an evaluation of the reproductive cycle in Solen marginatus (Pennant, 1777),(Bivalvia: Solenacea) in southern Tunisia. Journal of Shellfish Research, 29(1), 129-134.

Humason, G. L. (1962). Animal Tissue Techniques. Animal tissue techniques. United States of America: W. H. Freeman and Company.

Lavín, M. F., & Marinone, S. G. (2003). An overview of the physical oceanography of the Gulf of California. In Nonlinear processes in geophysical fluid dynamics (pp. 173-204). Springer.

Lawrence, J. M., & Herrera, J. (2000). Stress and deviant reproduction in echinoderms. Zoological Studies, 39(3), 151-171.

Maluf, L. Y. (1988). Composition and distribution of the central eastern Pacific echinoderms (6th ed.). Los Angeles, California: Natural History Museum of Los Angeles County.

Mercier, A., & Hamel, J. (2009). Endogenous and exogenous control of gametogenesis and spawning in echinoderms. Advances in Marine Biology, 55, 1-320.

Muthiga, N. A., Kawaka, J. A., & Ndirangu, S. (2009). The timing and reproductive output of the commercial sea cucumber Holothuria scabra on the Kenyan coast. Estuarine, Coastal and Shelf Science, 84(3), 353-360.

Nuño-Hermosillo, A. (2003). Ecología poblacional, ciclo reproductivo e historia de la pesquería del pepino de mar Isostichopus fuscus (Ludwig, 1875) (Echinodermata: Holothuroidea) en Bahía Chamela, Jalisco, México. Centro Universitario de Ciencias Biológicas y Agropecuarias. Universidad de Guadalajara, México.

Omar, H. A., Abdel Razek, F. A., Abdel Rahman, S. H., & El Shimy, N. A. (2013). Reproductive periodicity of sea cucumber Bohadschia vitiensis (Echinodermata: Holothuroidea) in Hurghada area, Red Sea, Egypt. Egyptian Journal of Aquatic Research, 39, 115-123.

Paden, C. A., Abbott, M. R., & Winant, C. D. (1993). Tidal and atmospheric forcing of the upper ocean in the Gulf of California: 2. Surface heat flux. Journal of Geophysical Research, 98(11), 20091-20103.

Pennington, J. T. (1985). The ecology of fertilization of echinoid eggs: the consequences of sperm dilution, adult aggregation, and synchronous spawning. The Biological Bulletin, 169(2), 417–430.

Ramofafia, C., Battaglene, S. C., Bell, J. D., & Byrne, M. (2000). Reproductive biology of the commercial sea cucumber Holothuria fuscogilva in the Solomon Islands. Marine Biology, 136(6), 1045-1056.

Ramofafia, C., Byrne, M., & Battaglene, S. (2001). Reproductive biology of the intertidal sea cucumber Actinopyga mauritiana in the Solomon Islands. Journal of the Marine Biological Association of the UK, 81(3), 523-531.

Reyes-Bonilla, H., & Herrero-Pérezrul, M. D. (2003). Population parameters of an exploited population of Isostichopus fuscus (Holothuroidea) in the southern Gulf of California, México. Fisheries Research, 59(3), 423-430.

Reyes-Bonilla, H., Ramírez-Ortiz, G., Herrero-Pérezrul, M. D., & Calderón-Aguilera, L. E. (2016). Subestimación de la abundancia del pepino café Isostichopus fuscus (Holothuroidea: Echinodermata) en muestreos diurnos con respecto a nocturnos en el. Revista Mexicana de, 87(2), 519-522.

Safran, P. (1992). Theoretical analysis of the weight-length relationship in fish juveniles. Marine Biology, 112(4), 545-551.

Saout, C., Quéré, C., Donval, A., Paulet, Y.-M., & Samain, J.-F. (1999). An experimental study of the combined effects of temperature and photoperiod on reproductive physiology of Pecten maximus from the Bay of Brest (France). Aquaculture, 172(3), 301-314.

SEMARNAT. (1994). Secretaría de Medio Ambiente y Recursos Naturales. Norma Oficial Mexicana NOM-059-ECOL-1994. México: Diario Oficial de la Federación (DOF), miércoles 16 de marzo de 1994.

SEMARNAT. (2002). Secretaría de Medio Ambiente y Recursos Naturales. Norma Oficial Mexicana NOM-059-ECOL-2001. Diario Oficial de la Federación (DOF), miércoles 6 de marzo de 2002.

SEMARNAT. (2010). Secretaría de Medio Ambiente y Recursos Naturales. Norma Oficial Mexicana NOM-059-SEMARNAT-2010. Diario Oficial de la Federación (DOF), jueves 30 de diciembre de 2010.

Servicio-Geológico-Mexicano. (2003). Carta Geológico-Minera Ensenada, H11-2 Baja California, Esc. 1:250 000.

Shepherd, S. A., Martinez, P., Toral-Granda, M. V, & Edgar, G. J. (2004). The Galápagos sea cucumber fishery: management improves as stocks decline. Environmental Conservation, 31(2), 102-110.

Sonnenholzner, J. I. (1997). A brief survey of the commercial sea cucumber Isostichopus fuscus (Ludwig, 1875). SPC Beche-de-Mer Information Bulletin, (March), 1995-1998.

Sonnenholzner, J. I., Brandt, M., Francisco, V., Hearn, A., Luzuruaga, M., Guarderas, P., & Navarro, J. C. (2013). Echinoderms of Ecuador. In J. J. Alvarado & F. A. Solís-Marín (Eds.), Echinoderm Research and Diversity in Latin America (pp. 183-233). Springer-Verlag Berlin Heidelberg.

Sonnenholzner, J. I., Searcy-Bernal, R., & Panchana-Orrala, M. (2017). The potential for propagation of the commercial sea cucumber Isostichopus fuscus (Ludwig, 1875) by induced transverse fission. Regional Studies in Marine Science, 9, 35-42.

Soto-Mardones, L., Marinone, S. G., & Parés-Sierra, A. (1999). Variabilidad espariotemporal de la temperetura superficial del mar en el Golfo de California. Ciencias Marinas, 25(1), 1-30.

Tapia Vázquez, O. M., Castro González, J. J., & Valles Ríos, H. (1996). Madurez gonádica del pepino de mar Parastichopus parvimensis en la costa occidental de Baja California, México, en 1994. Ciencia Pesquera, 12, 5-12.

Toral-Granda, M. V., & Martínez, P. C. (2007). Reproductive biology and population structure of the sea cucumber Isostichopus fuscus (Ludwig, 1875) (Holothuroidea) in Caamaño, Galápagos Islands, Ecuador. Marine Biology, 151(6), 2091-2098.

Supplementary material

Table 1

Descriptive statistics of the size (TL = total length, cm and W = weight, g) of Isostichopus fuscus

on the east coast of Baja California

Month		N	Mean	Mode	SD	Min	Max		Month		N	Mean	Mode	SD	Min	Max
Oct 2014	TL	43	18.54	21	5.51	5	28		Jan 2016	TL	160	19.40	21	4.27	8	32
	W	45	280.23	259	137.24	6	616			W	158	231.08	250	132.26	20	760
November	TL	122	21.42	19	4.82	7	34		February	TL	93	22.67	25	5.41	12	36
	W	132	312.05	378	153.48	18	746			W	93	360.98	300	211.26	60	1 000
December	TL	32	22.03	21	6.86	9	41		March	TL	140	22.80	22	5.60	9	36
	W	32	360.16	297	194.57	28	739			W	140	325.79	400	197.90	20	920
Jan 2015	TL	154	15.71	11	5.9	6	35		April	TL	70	25.64	26	4.16	16	35
	W	154	189.97	91	181.51	13	1 188			W	70	627.71	500	211.61	150	1 080
February	TL	79	26.1	25	4.31	14	35		May	TL	136	19.43	18	4.39	9	31
	W	79	563.54	600	183.1	80	980			W	136	278.16	360	156.84	50	840
March	TL	122	18.2	12	6.44	5	35		June	TL	152	22.10	18	4.15	10	38
	W	122	266.31	400	214.1	20	800			W	152	373.03	400	180.64	50	1 200
April	TL	92	26.79	30	5.72	14	40		July	TL	163	21.03	17	5.80	4	36
	W	92	588.48	600	216.57	40	1 000			W	162	416.54	260	274.41	10	1 200
May	TL	74	17.7	18	4.72	7	29		August	TL	86	26.65	26	4.24	9	37
	W	74	184.38	100	183.14	4	800			W	86	745.23	700	206.17	30	1 080
June	TL	321	18.56	16	6.4	7	35		September	TL	148	19.65	21	4.53	6	31
	W	321	221.72	100	203.51	10	900			W	148	324.46	300	141.26	60	820
July	TL	200	23.36	21	5.69	6	38		October	TL	125	22.06	20	5.05	9	33
	W	200	374.22	300	215.51	8	1 200			W	125	589.52	800	268.43	60	1 200
August	TL	133	22.51	25	7.33	9	40		November	TL	121	23.08	27	5.27	7	34
	W	133	345.94	100	235.4	20	1 000			W	121	642.31	840	269.84	10	1 120
September	TL	40	28.58	26	4.65	16	36		December	TL	38	25.82	27	4.02	15	32
	W	40	556	500	133.22	50	800			W	38	427.89	440	153.84	120	800
October	TL	30	23.93	26	3.12	14	28		Pooled Average	TL		21.46	21	6.04	4	41
	W	30	399	400	46.19	280	500			W		375.64	400	249.18	4	1 200

Table 2

Description of the gonadal development stages of Isostichopus fuscus

Gonadal stage	Females	Males
II Gametogenesis	Oogenesis Connective tissue thick, lax and undulating Visible hemal sinus Activity in germinal epithelium Lumen empty or half empty Vitellogenic oocytes Primary oocytes	Spermatogenesis Connective tissue thick, lax and undulating Visible hemal sinus Activity in germinal epithelium Lumen empty or half empty Columns of primary spermatocytes orderly in wavy form
III Mature	Mature Connective tissue not visible Hemal sinus not visible Little or no activity of the germinal epithelium Lumen filled with mature oocytes	Mature Semi lax connective tissue, poorly visible Hemal sinus not visible Little or no activity of the germinal epithelium Lumen filled with mature spermatozoa
IV Spawning	Spawning Connective tissue visible, fine and lax. Sinus hemal barely visible. Null activity of the germinal epithelium. Half-empty lumen with oocyte concentration in the center. Residual oocytes. Phagocytes	Expulsion Connective tissue visible, fine and lax. Sinus hemal barely visible. Null activity of the germinal epithelium. Half-empty lumen, agglomeration of spermatozoa in the center and empty spaces in the periphery Phagocytes
V Post spawning	Post spawning Connective tissue thick, lax and undulating Visible hemal sinus Little or no activity of the germinal epithelium Lumen empty or half empty Oocytes residual and in lysis Phagocytes	Post expulsion Connective tissue thick, lax and undulating Visible hemal sinus Little or no activity of the germinal epithelium Lumen empty or half empty Spermatozoa residual and in lysis Phagocytes