Until now, some aspects of the reproductive biology of H. leucospilota from other regions have been studied; however, there are differences from one population to another, as well as differences in the environmental conditions of their habitats. Therefore, it was necessary to investigate the Persian Gulf population in order to properly manage fisheries. According to Sewell (1992), reproductive biology may even differ among specimens of the same species collected at different latitudes.
H. leucospilota is a common species in the Persian Gulf that is gonochoric (having the sexes separate) and show no sexual dimorphism which is observed in many aspidochirotida holothurians such as H. fuscocinerea (Benitez-Villalobos et al. 2013), H. poli (Slimane-Tamacha et al. 2019) H. spinifera (Asha and Muthiah 2007) and H. (Platyperona) sanctori (Mezali et al. 2014). Drumm and Loneragan (2005) also reported no external sexual dimorphism for H. leucospilota feom Rarotonga. In the present study, the sex was determined using the mature gonad color (testes are beige and the ripe ovaries are pink/red (Purwati and Luongvan 2003; Drumm and Loneragan 2005)) and gonad histology.
Holothurians often exhibit an annual reproductive cycle (Chao et al. 1995); however, semiannual or even continuous reproductive cycles throughout the year are also common, especially in tropical regions (Conand 1993). Similar to H. poli from Lesvos Island in Greece (Bardanis and Batjakas 2018), the lack of hermaphroditic individuals and existence of an annual reproductive cycle were noted in the present study. In the present study, the annual reproductive cycle of H. leucospilota was determined using the histological and sex hormone analysis and the effect of environmental variables (including the salinity and temperature) on it was also evaluated. Based on the research, the main environmental stimuli that directly or indirectly affect the sexual cycles of organisms include changes in temperature, salinity, photoperiod, availability of food and habitat, and recruitment and survival of juveniles (Costelloe 1985). Among these, two factors (salinity and temperature) were investigated in the present study. The relationship between the reproductive cycle and natural environmental conditions has been investigated for a number of holothurians (Costelloe 1985; Sewell and Bergquist 1990; Hamel et al. 1993; Hopper et al. 1998; Asha and Muthiga 2008; Benitez-Villalobos et al. 2013; Kazanidis et al. 2014; Panola-Madrigal et al. 2017). Like H. Tubulosa and H. fuscocinerea, H. leucospilota reproduces in the warmer months when the days are longer (Benitez-Villalobos et al. 2013; Kazanidis et al. 2014). In the Persian Gulf H. leucospilota population, gamete release coincided with a significant increase in temperature (when temperatures reach 33.9°C in July-August. Similarly, H. tubulosa spawns during the warm months (from July to September) in the Adriatic Sea, when surface water temperature is at its warmest (Despalatovic et al. 2004). The main reason is that the larvae have more food available during this time (warm months) due to the phytoplankton proliferation (Kazanidis et al. 2014; Panola-Madrigal et al. 2017). According to Conand (1993), the seasonal reproductive cycle in holothurians is characteristic of most temperate species, whose spawning season occurs with increasing temperature, which is associated with increasing phytoplankton biomass and, in turn, more food source for sea cucumber. Unlike these species, the spawning period of some holothurians such as Aslia lefevrei (Costelloe 1985), H. forskali (Ballesteros et al. 2021)d whitmaei (Shiell and Uthicke 2006), occurs in colder months. Anyway, considering the stimulation of spawning by thermal changes, it seems that temperature is probably the most important factor in the release of gametes in many holothurian species, and for this reason, temperature changes are used abundantly to stimulate spawning in holothurian cultivation (Rakaj et al. 2018; Ballesteros et al. 2021). On the other hand, Vargas-Yanez et al. (2017) stated that temperature might not affect the gonadal cycle of Parastichopus regalis collected from depths of 60 to 150 m in the northwestern Mediterranean Sea, because the water temperature at these depths is quite constant and about 13–14 ◦C in winter, spring and summer when the reproductive cycle develops. In autumn, when the temperature decreases with depth (from 18°C at 50 m to 13°C at 150 m), no gonads were found in the studied P. regalis (Vargas-Yanez et al. 2017).
Salinity could be another factor affecting the reproductive cycle in holothurians. However, there was no correlation between gonadal development in H. leucospilota and salinity in the present study. Analysis of salinity fluctuations in the Persian Gulf based on 10-year reports shows that changes in salinity in the region have a constant trend. The average salinity is 39.12 PSU during the year; it decreases to a minimum of 38.68 PSU in summer and increases to a maximum of 39.34 PSU in winter (Omidi and Norinezhad 2009). Therefore, salinity does not seem to have a significant effect on the reproductive cycle of H. leucospilota.
In the present study, H. leucospilota did not show the characteristic pattern of a continuous reproductive cycle. This is to be expected given that it is a tropical species and considering that most cycles are timed so that larval or juvenile production coincides with periods favorable for feeding or survival (Giese and Pearse 1974). Slimane-Tamacha et al. (2019) reported the same result with H. poli that exhibits a single annual reproductive cycle. On the other hand, Gaudron et al. (2008) reported a biennial cycle consisting of a dominant warm-season spawning and a smaller and more variable secondary peak in H. leucospilota in the western Indian Ocean.
In the present study, the gonad morphology and gametogenesis stages were similar to those obtained with holothurians such as H. leucospilota in the western Indian Ocean (Gaudron et al. 2008). the gonad of H. leucospilota from the Persian Gulf was similar to H. fuscocinerea (Benitez-Villalobos et al. 2013)d poli (Slimane-Tamacha et al. 2019), consisting of two bundles of long tubules, one on each side of the dorsal mesentery. However, the results of the present study differed in some aspects from the results obtained by other researchers even on the same species in other regions. Such as H. fuscocinerea (Kazanidis et al. 2014), no sexual dimorphism was observed in the color of the gonads of H. leucospilota from the Persian Gulf. This was in contrast to Gaudron et al. (2008), who noted sexual dimorphism for gonad color in H. leucospilota in the western Indian Ocean, a feature usual in many Pacific holothurians. In addition, gonad tubules were absorbed after spawning in Stichopus californicus (Smiley 1988), while in H. leucospilota (such as H. leucospilota in the western Indian Ocean (Gaudron et al. 2008)d fuscocinerea (Kazanidis et al. 2014)) the tubules were not absorbed and progenitor cells were visible along the germinal epithelium of those at the end of next January and growth continued.
In all holothurian reproductive research, although researchers may have used different names to express the stages of gonad development (e.g. two stages in a stage 1), they all express the stages of rest, growth, maturation, spawning, and exhaustion. The annual reproductive cycle of H. leucospilota culminated in spawning in mid-summer. Such as Cucumaria frondosa (Hamel and Mercier 1996), gametogenesis in H. leucospilota started in late January, after a short period of recovery. Intense synthesis of gametes and their storage continued until spawning, when gametes were released from the mature gonad tubules, which had reached their maximum size in late July to mid-October. This growth period (on average 6–8 months) appears to be necessary for each tubule to reach maturity. Such results were quite different from those reported by Gaudron et al. (2008) for the same species in the western Indian Ocean. In the present study, male and female H. leucospilota were mainly in stages 4 and 5 from July to mid-October, while at the same time, the western Indian Ocean H. leucospilota were reported in stage 3 (Gaudron et al. 2008). A Higher percentage of H. leucospilota were in post-spawning stage by the end of October, while the western Indian Ocean H. leucospilota were reported post-spawning stage at the end of February (Gaudron et al. 2008). Some western Indian Ocean H. leucospilota remained at this stage until April followed by a second spawning event in May. Gaudron et al. (2008) suggested that the occurrence of the second spawning event could be because the females that failed to spawn in February spawned in suitable environmental conditions during the wet season of the Australian summer. In agreement with Gaudron et al (2008), in the present study, gametes were in the maturation stage and distribution of oocyte diameters was more homogeneous in mature male and female H. leucospilota. While, coexistence of different stages of developing gametes within one tubule was usual in male and female growing H. leucospilota from the Persian Gulf (The present study) and the western Indian Ocean (Gaudron et al. 2008). In the present study, about 75–90% of male and female H. leucospilota had gonads every month of the year. However, the gonads were very small and undistinguishable in lower than 25% samples during the year especially during December to mid-January. It is also common to find a high proportion of specimens with gonads for several months among other sea cucumber species. Christophersen et al. (2020) showed that 67–100% of Parastichopus tremulus samples had gonads every month of the year.
Based on the microscopic observations, the gonad development in H. leucospilota male and female conforms to a synchronous "tubule recruitment pattern" and all gonad tubules are at the same stages of development in this species. It seems to be due to changes in water temperature and photoperiod (Drumm and Loneragan 2005). This is not unusual and many other holothurian species do not follow the model of an asynchronous "tubule recruitment pattern", such as H. fuscocinerea (Benitez-Villalobos et al. 2013). The same result was reported by Gaudron et al. (2008) on H. leucospilota in the western Indian Ocean. On the other hand, the asynchronous model is usual in some other holothurians such as H. atra (Chao et al. 1994), H. fuscogilva and H. mauritiana (Ramofafia and Byrne 2001).
Body size is an important factor to explain changes in reproductive function and index, both interspecies and intraspecies. For most holothurian species, with the increase in body size and weight, the weight of the gonads and thus the reproductive index increases (Conand 1993). In the present study, female and male H. leucospilota with a body weight between 61 and 70 gr and GSI more than 12% had gonads at stage III (Growing stage) of gametogenic development indicating that they were approaching sexual maturity. In addition, in H. leucospilota with a body weight between 72 and 85 gr and GSI more than 19%, the gonads were fully mature in stage VI of sexual maturity. However, these are more than the mean body weight reported for growing and mature H. leucospilota in southern Vietnam (45 gr and 67 gr, respectively) (Nguyen and Britaev 1993). Drumm and Loneragan (2005) reported female and male H. leucospilota with body weights more than 37.5 g and 55 were in growing and mature gonad development stages, respectively. Furthermore, H. leucospilota had the lowest values of GSI (between 1.2 to 2.7%) in resting months (late-October to early-January), during which the ‘sexually undifferentiated’ individuals were observed. Benitez-Villalobos et al. (2013) also recorded the lowest values of GSI (< 2) for Holothuria fuscocinerea during October to January.
In many holothurian species, the mean gonad index (MGI) has been used to determine the gonad growth pattern (Conand et al. 2002; Foglietta et al. 2004; Christophersen et al. 2020). In the current study, the MGI was used to document the development pattern of gonads, and due to the relatively low dispersion of MGI between individuals in both males and female H. leucospilota, we recommend its use to obtain suitable information on the reproductive cycle. However, in some researches on holothurian reproduction, the use of this index has not been successful. For example, Ramon et al. (2022) did not suggest the use of this index due to the high dispersion of the results in P. regalis, especially in summer when the gonad reached its maximum size. During the reproductive cycle (except from August to the end of November), the mean gonadal index in female H. leucospilota was higher than that of males. Tuwo and Conand (1992) also stated that female H. forskali had higher gonadal indices than male. Inversely, the mean gonadal index was reported to be higher in male Cucumaria frondosa (Hamel and Mercier 1996) and Psolus fabricii (Hamel et al. 1993) than that of female, before the spawning event. On the other hand, there was no inequality between male and female gonadal indices in holothurians such as Parastichopus californicus (Cameron and Fankboner 1986) and Stichopus mollis (Sewell and Bergquist 1990). It seems that the female gonads had larger gonads and higher MGI due to the need for more energy to synthesize an equivalent amount of gametes, especially considering that no significant difference was observed in the body size of males and females during the reproductive cycle. Lawrence and Lane (1982) showed that ovaries are richer than testes even though they are of similar size. According to Hamel and Mercier (1996), during periods of low food availability, the energy content decreases in the female body wall, even though its mass does not differ significantly. This indicates that a metabolic activity that requires energy is taking place and probably at the same time, (the early stage of gametogenesis), the transfer of energy to the gonads occurs. Synthesis of female gamete probably initially needs more energy than males (Hamel and Mercier 1996).
The oocyte size is another important factor in the reproductive cycle. In the holothurians studied, the size of the mature oocytes varies between 65–150 µm (Conand 1993; Hamel and Mercier 1996). In the present study, the mean diameter of the mature oocytes (a parameter that is commonly used in the population reproduction study (Conand 1993)) in H. leucospilota reached its maximum (86.57 ± 18 µm) in stage VI of maturity (May to July). This result was comparable with the reports about mature oocyte diameter in H. atra (88 µm; Harriott 1982); H. parvula (93 µm; Emson and Mladenov 1987).
The present study showed the existence of estradiol and testosterone in the gonad tissues of H. leucospilota. This helps build evidence for steroids in invertebrates, including holothurians. Although the presence of steroids has been proven in a limited number of species of holothurians such as H. scabra (Thongbuakaew et al. 2021), the amount of these hormones in these species has not been determined. The results presented here are consistent with data reported on the presence of T and E2 in other echinoderms such as the sea urchins Paracentrotus lividus (Barbaglio et al. 2007)d lividus (Janer et al. 2005). Such as other species, a close relationship was observed between T or E2 levels and the degree of gonad growth and gametes maturation in H. leucospilota in the present study. A significant difference in the concentration of T and E2 was found among the months as well as between males and females in H. leucospilota. This result contrasted with the results reported by Barbaglio et al. (2007) who found no clear correlation between T and E2 levels and distribution of reproductive stages along the year in sea urchin P. lividus. On the other hand, as in the present study, Janer et al. (2005) reported that the T level was significantly lower in the testis of P. lividus in the developmental stage (during early and advanced spermatogenesis) than in the mature stage (the end of gametogenesis). This result shows the role of T in regulating late spermatogenesis and spawning. T levels also begin to increase immediately after the resting (recovery) phase (when the gonad is preparing to reorganize for the next reproductive cycle), indicating the involvement of T in the initiation of a new spermatogenesis, as has been reported for sea urchin Lytechinus variegatus (Wasson et al. 2000). Hines et al. (1992) stated that, T levels are also elevated during the developmental stage in ovaries, during vitellogenesis, as observed in sea star Asterias vulgaris, thus indicating involvement of T also in oogenesis. In the present study, mean E2 concentration was more than T levels in females than in males, possibly reflecting a more important role of this hormone in females. As also reported in L. variegatus (Wasson et al. 2000)d vulgaris (Hines et al. 1992), mean E2 levels increased in ovaries immediately after recovery stages and it was high in the growth phase, during villogenesis, which indicates the special role of this hormone in the regulation of oogenesis. According to the data, the significant role of T and E2 in the reproductive biology of echinoderms, especially sea cucumber, seems evident.