Reproductive characteristics of the hermaphroditic four-finger threadfin, Eleutheronema tetradactylum (Shaw, 1804), in tropical coastal waters
BMC Zoology volume 8, Article number: 22 (2023)
This study investigated the reproductive traits of the hermaphroditic four-finger threadfin, Eleutheronema tetradactylum, along the coasts of Thailand during January to December 2021. Fish samples were collected from Pattani Bay, Thailand to assess the sex ratio, gonadosomatic index (GSI), maturity stage and fecundity. Additional fish samples were also collected from other areas to evaluate the length and weight at first sex change (Ls50 and Ws50) and length at first maturity (Lm50). The overall sex ratio for male and female was 1:0.69 with male being predominant throughout the year. Threadfin fish spawn the whole year round with peaks during moderate rainy and heavy rainy seasons. Histological examination confirmed its protandrous hermaphrodite posing multiple spawning habits. The average fecundity was 1.85 × 105 ± 1.05 × 105 eggs and positively related with standard length, body weight, gonad weight, and egg diameter (p < 0.05). The Ls50 and Ws50 were 27.58 cm and 419.39 g, and 29.71 cm and 457.28 g, for fish from Pattani Bay and Samut Prakan province, respectively. The Lm50 of male from Pattani Bay and Samut Prakan province were 25.78 cm and 25.56 cm, respectively, which were larger than those from Satun and Nakhon Sri Thammarat provinces. The Lm50 of females from Pattani Bay was smaller than that from Samut Prakan province. This study provided fundamental information on the reproductive characteristics of E. tetradactylum, which can be implemented to support management of natural fish stock and aquaculture development.
Information on fish reproductive biology is crucial for managing the fishery resources and aquaculture [1, 2]. Several parameters such as sex ratio, gonadosomatic index (GSI), maturity stage, fecundity, size and weight at first sex change, size at first maturity and relationship between various parameters with reproductive traits have been reported in many fish species, e.g., neotropical fish; Astyanax fasciatus, Oligosarcus robustus, Loricariichthys anus, Trachelyopterus lucenai  and other 1207 marine fish species from Brazil . The four-finger threadfin, Eleutheronema tetradactylum, is widely distributed along the coasts of India, Bangladesh, Sri Lanka, the Indo-West Pacific and the northern and western Australia [3,4,5]. Its habitats range from coastal waters to estuaries in the tropics [5,6,7,8]. It has long been a target fish species by fishermen in many parts of the world, including Thailand, with high value and market demand . The fishing gears used typically to capture this fish include gillnet, beach seine, longline, trap and trawl .
Many studies were conducted on the biology and fisheries of this species including its breeding , length at maturity size , fecundity , histology of gonads , stock structure and exploitation status [10, 14, 15], aquaculture , taxonomy [3,4,5] and compositions of parasitic organisms [17, 18]. However, most studies on reproductive aspects focused only on a particular reproductive trait viz., length and weight relationship [10, 11, 14], fecundity [8, 12], gonad histology , size at first maturity , and size at the first sex change . Additionally, some studies reported that this fish was found from late winter to mid-summer in Australia , March to September in Malaysia , January to April and July to September in India , and February to March and July to August in the Bay of Bengal . In Thailand, Sritakon et al.  found that this fish spawned throughout the whole year with three peak periods. Nevertheless, none of them simultaneously comprehensively investigated the reproductive trait from the wild population. Moreover, it is known that sex change is an evolutionary adaptation of organisms to enhance a successful reproduction [21, 22]. E. tetradactylum is a protandrous hermaphrodite fish, beginning as male and later changing to female. The sizes at the first sex change were reported from different areas. For instance, the Ls50 was 32.50 cm of total length (TL) from the north-western Australia , and 20.8 cm to 46.5 cm FL across the northern Australia of fork length (FL) from the northern Australia . However, no existing data on the size at first sex change from the Southeast-Asian region is currently available. This study was thus conducted to fulfill the gap and provide scientific knowledge on reproductive biology of E. tetradactylum. This is considered the first comprehensive investigation on the reproductive traits of this species.
Materials and methods
The peninsular of Thailand is situated in between the Andaman Sea on the west side and the Gulf of Thailand on the east side. Both coastal regions provide high diversity of fish and are regarded as important fishing grounds [23,24,25]. The Gulf of Thailand in this study can be divided into two parts: the inner Gulf of Thailand and the southern Gulf of Thailand. The former is mainly influenced by river runoff providing high diversity of fish while the latter is influenced by intrusion of the South China Sea. Four different study sites were selected from both the Gulf of Thailand and Andaman coasts, including the provinces of Pattani (PTN), Nakhon Sri Thammarat (NST) and Samut Prakan (SPR) and Satun (STN) (Fig. 1). Study area was illustrated by QGIS software, version 3.32.1 . Satun province is situated in the coast of Andaman Sea, while Samut Prakan is in the inner Gulf of Thailand, and Nakhon Sri Thammarat and Pattani are in the southern Gulf of Thailand. Among these, Pattani Bay situated in Pattani province is selected as the main study site since it is known as a major fishing ground by local fishermen especially for four-finger threadfin [9, 27, 28], mackerel fish [29,30,31], shrimps  and crabs . Three different seasons were classified in the bay: dry or pre-monsoon from January to April, moderate rainy season (SW monsoon) during May to September, and heavy rainy season or NE monsoon season during October to December .
Sampling was conducted monthly from January to December 2021 by a traditional monofilament gill net (mesh size of 4.5 cm stretched, 5 m deep, and 540 m long) in Pattani Bay. The nets were set and left floating for 60 min before being hauled onboard. The sampling period was conducted between 0500–0900 h. The captured specimens, which died immediately after capture, were kept in the ice box and transported to the laboratory at the Faculty of Science and Technology, Prince of Songkla University, Thailand. Approximately 30 specimens of E. tetradactylum were selected from a monthly sample for further analysis. Additionally, fish samples from the provinces of Satun (STN), Nakhon Sri Thammarat (NST) and Samut Prakan (SPR) were collected from the catches of local fishermen (Fig. 1).
Standard length (SL) and body weight (BW) were measured to the nearest 0.01 cm and weighted to 0.01 g, respectively. Consequently, fish samples were later categorized into four size classes: small (18.13 ± 1.71 cm), medium (25.74 ± 2.19 cm), large (32.96 ± 2.06 cm) and extra-large (42.07 ± 2.20 cm) based on the length frequency distribution measured by Past software .
A total of 408 fish samples, excluding 21 transitional sex fishes, from Pattani Bay were used to assess the sex ratio, gonadosomatic index (GSI), maturity stage, spawning period and fecundity. The length and weight at first sex change (Ls50 and Lw50) were calculated for fishes from Pattani Bay (n = 429) and Samut Prakan province (n = 39), whereas the length at first maturity (Lm50) of male and female fish (n = 721) were compared for all four sites (Table 1).
After removing the gonad from the body cavity, sex and maturity stage of fish samples were distinguished macroscopically with the adapted criteria proposed by Laevastu . It was assigned as male, female or transitional sex based on the presence of both testis and ovary with the relevant literatures [7, 13]. For transitional gonads, the testicular tissue appeared along the dorsal edge and in the inner lateral part of the gonads while the ovarian tissues were along the ventral side and outer part of the gonads. Characteristics of gonadal stages for male and female and transitional sex are given in Tables 2 and 3.
For histological examination, 95 gonads, accounting for 22% of total fish samples, were wiped with tissue paper and weighed to the nearest 0.01 g. They were immediately fixed with 10% buffered formalin and stored in bottle jars [37,38,39]. Immature individuals for which sex could not be identified with the naked eyes were excluded for this histological analysis. In the histological process, three sub-samples from gonad tissues were taken from three parts of the gonad: anterior, middle, and posterior, respectively. The protocol for further analysis was based on the methodologies described by Zamidi et al. , Shihab et al.  and Simon et al. . Gonad tissues were loaded in an automatic tissue processor (LEICA JUNG HISTOKINETTE 2000) with the process of 12 solutions consisting of propanol, chloroform and paraplast for 24 h. After cleaning and dehydrating processes, gonad tissues were embedded in the paraffin wax, shaped into blocks and trimmed. The tissues were then transversely sectioned into 4–5 µm slices by a semi-automatic microtome (LEICA Mod. 2035 Biocut). Those tissues were later defiled with hematoxylin and eosin, mounted on glass slides for light microscopy analysis. Finally, gonadal maturation was classified into eight developmental phases viz., stage I (immature), stage II (resting), stage III (developing), stage IV (maturing), stage V (mature), stage VI (spawning), stage VII (spent) and stage VIII (recovering) as indicated by Pember et al. , and Shihab et al. .
The terminology applied in this study was based on the criteria used by Pember et al. , Wallace et al.  and Brown‐Peterson et al. . In male, oocyte stages included connective tissue (ct), lumen (l), spermatocytes (sc), sperm duct (sd), spermatogonia (sg), sperm sinuses (ss), spermatids (st), spermatozoa (sz). In transitional gonads, they were connective tissue (ct), nucleolar oocytes (cn), perinucleolar oocytes (pn) and spermatozoa (sz). In female, they were atretic oocytes (ao), cortical alveolar oocytes (ca), nucleolar oocytes (cn), connective tissue (ct), lipid droplets (ld), ovigerous lamellae (ol), primary growth oocytes (pg), post ovary follicle (pof), perinucleolar oocytes (pn), and yolk granule oocytes (yg).
To define a monthly reproductive phase of fish, a gonadosomatic index (GSI) of male and female was calculated with the following equation :
where GSI is gonadosomatic index, GW is gonad weight (g) and TW is total weight (g). A non-parametric one-way ANOVA (Kruskal–Wallis test) was performed to examine the homogeneity of the monthly means of GSI and number of mature fish using the R program .
To determine spawning or breeding season, seasonal variation of percentage of male and female at different gonadal stages was classified as immature (stage I-IV) and mature (stage V-VIII) based on visual observation. The stages of V and above were grouped into mature fish as their readiness to spawn. A spawning season was determined by a high percentage of mature fish .
To determine fish fecundity, the posterior, middle and anterior parts of mature female gonads from preserved samples were gently dissected. Eggs were carefully put in a petri dish, weighed to the nearest 0.001 g, and poured with distilled water. The number of eggs was later counted by gravimetric method based on the weight of ovary and density of egg . A fecundity was estimated by a gravimetric method as:
where, n = number of eggs in three subsamples, OW = total weight of ovary (g) and WS = weight of three subsample (g). Egg diameter was measured to the nearest 0.1 µm by microscope with the assistance of Digimizer image analysis software version 4.3.4. Then, variances of fecundity among fish size and region were examined. Fish samples from Satun province (STN) were excluded due to an absence of developing/mature ovaries.
Relationship between fecundity (F) and standard length (SL), body weight (BW), gonad weight (GW) and egg diameter (ED) were tested by linear regression analysis . Prior to the analysis, raw data were transformed to Log (x + 1) to reduce non-normality. Then, relationship between fecundity (Y) and independent variables (X) were analyzed with PAST software  by the following equation:
where, a and b are the constant values from linear regression, F is fecundity and X is independent variables (SL, BW, GW and ED). Linear relationship of egg diameter (ED) to SL, BW, GW and fecundity were assessed.
To determine length at the first sex change (Ls50), a logistic regression was fitted to the proportions of both male and female in each of the 1 cm length class , then applied with PAST software  by the following function :
where, Pls is the proportion of male and female in each interval class S, S50 and S95 are the lengths at which 50% and 95% of the populations are females, respectively. In addition, a similar calculation was applied for weight at the first sex change (Lw50) in each of the 10 g of weight class. Due to insufficient samples, fish from STN and NST provinces were excluded for this analysis.
To estimate length at the 50% of fish maturity (Lm50), frequency distribution and proportion of both male and female were binned into the 1 cm length class . They were then analyzed by the following logistic function :
where, Plm is the proportion of mature fish in each 1 cm length class M, M50 and M95 are the length at which 50% and 95% of the populations are matured, respectively. Due to insufficient samples of female fish from STN and NST, the analysis of Lm50 of females for those sites were excluded.
Size and sex
Based on the examined samples, the smallest fish was 11.96 cm SL and 36.86 g BW, and the largest was 48.30 cm SL and 1900 g BW. The majority of fish standard length ranged between 20–30 cm SL, accounting for 49.51% of all samples (Table 4). Sex could be allocated into three categories viz., male, female and transitional sex. The average SL and BW of males, females and transitional sex were 22.36 ± 4.96 cm, 30.43 ± 5.29 cm and 28.49 ± 3.68 cm, and 229.69 ± 168.90 g, 520.92 ± 265.47 g and 436.91 ± 168.54 g, respectively. Spatially, the lowest SL (19.38 ± 2.86 cm) and BW (137.51 ± 67 g) were found in NST, while the highest SL (30.02 ± 4.86 cm) and BW (442.76 ± 205.4 g) were in SPR.
The overall sex ratio between male and female was 1:0.69 (χ2 = 14.16, p < 0.001) (Table 5). Monthly sex ratio varied sporadically, including in January, February, April, May and August (χ2 > 3.84, p < 0.05).
Gonadosomatic index (GSI) and maturity
The Kruskal–Wallis test indicated different GSI values among months (p < 0.05). Different patterns of mean GSI between male and female were observed (Fig. 2). The GSI of males was the highest in July and dropped gradually from September onwards. The GSI of females demonstrated five times of peak within a year, indicating multiple spawning behavior. The highest number of mature males was observed in July while mature females in May (Fig. 3). It coincided with the peak of male’s GSI which was observed in July as the observation of high preponderance of mature male. Kruskal–Wallis test indicated that there was no variation between the number of mature males and females among months (p > 0.05).
There were eight maturity stages observed in males but no resting stage (Stage II) found in females during the period of this study (Fig. 3). Overall, 91.62% of males and 73.66% of females were found to be from immature to maturing stages throughout the year. The highest percentage of mature males (46.15%) and females (80%) occurred in July and April, respectively (Fig. 4). Occurrence of mature stages of females (Stage V-VIII) during May to October and December indicated that this fish spawned throughout the whole year.
Stage of gonad maturity
Gonad maturation could be categorized into male (Fig. 5a), transitional sex (Fig. 5b) and female (Fig. 5c). The standard length of fish ranged from 11.96—44.32 cm SL in male, 22.07—37.57 cm in transitional sex and 22.37—48.30 cm in female. Histological analysis indicated that the presence of spermatogonia (sg) predominantly occurred in immature male gonads (Fig. 6a). It was distinguished by a slightly stained cytoplasm and homogenous nuclei. The connective tissue (ct), spermatocyte (sc), spermatids (st) and spermatozoa (sz) appeared in the resting, developing and maturing stages (Fig. 6b, c and d, respectively). These stages represented an active spermatogenesis, and a number of spermatozoa could be recognized. The presence of lumen (l) in male gonads appeared in the maturing stage (Fig. 6d). When E. tetradactylum became mature, other characteristics including sperm sinuses (ss) and sperm duct (sd) appeared in the gonad wall (Fig. 6e). Sperm duct and majority of the testes were full of spermatids (st) and spermatozoa (sz). The presence of crypts containing spermatozoa (sz) was usually captive in the outer perimeter of testis in the spawning stage (Fig. 6f). The largest numbers of spermatogonia (sg) and spermatozoa (sz) were found in the spent stage.
However, the spaces of empty sperm were observed within the connective tissue (Fig. 6g). In the recovering stage, testes are dominated by connective tissues (ct) consisting of extensive spaces with a few appearances of spermatozoa (sz) (Fig. 6h). This appearance of spermatozoa indicated an undergoing spermiation of the gonad. In transitional gonads, the testicular and ovarian tissues were separated by the connective tissue (ct). The presence of chromatin nucleolar oocyte (cn) and perinucleolar oocytes (pn) in early transitional gonads occurred with an interseptum of connective tissues (Fig. 7a). While in the late transitional period, size of testicular region was reduced and gonad assumed as a typical characteristic of ovaries with large numbers of perinucleolar oocytes (pn) (Fig. 7b). In female gonads, chromatin nucleolar oocyte (cn), perinucleolar oocytes (pn) and ovigerous lamellae (ol) are present in immature ovaries with small nucleoli (Fig. 8a and b). The cytoplasm was strongly basophilic in this stage. During the ovary developing stage, cortical alveolar oocytes (ca) presented clearly (Fig. 8c) and the ovaries began to develop but were not ready to spawn. Maturing and mature stages were notable by the presence of cortical alveolar (ca) and abundance of yolk granule oocytes (yg) with several lipid droplets (ld). The number of yolk granule oocytes (yg) were large (Fig. 8d and e, respectively). Primary growth oocytes (pg), atretic oocytes (ao) and some ripped eggs (rp) were found in the spawning stage (Fig. 8f). Remnant yolk granule oocytes (yg) were recorded in the spent stage and atretic oocyte (ao) and connective tissues (ct) were typically detected (Fig. 8g). Connective tissue (ct) and ovigerous lamellae (ol) appeared in the recovering stage (Fig. 8h) indicating the continued process of ovulation.
A total of 78 female fishes with developing and maturing ovaries, SL ranging from 22.37 to 48.30 cm, were estimated for fecundity. It was found that the fecundity varied between 0.49 × 105 to 4.96 × 105 eggs, with an average of 1.85 × 105 ± 1.05 × 105 eggs and increased along an increment of fish size (Table 6). The relationship between fish fecundity and egg diameter was observed in which the higher fecundity posed the greater egg diameter. Fecundity also differed between study sites (p < 0.05) where fish from SPR had the highest compared to fishes from the other two regions, PTN and NST (Table 6). Positive logarithmic relationships of fish fecundity (F) to SL, BW, OW and ED from PTN were observed (Fig. 9). Positive correlations of fecundity to gonad weight (R2 = 0.86) and ED to fecundity (R2 = 0.52) were greater than other body parameters (Figs. 9 and 10). Moreover, a positively higher correlation between ED and GW was found compared to other body parameters.
Length and weight at first sex change and size at first maturity
Results from logistic regression on the proportion of males and females indicated variations of Ls50 and Lw50 from two different regions. The Ls50 and Lw50 of fish from PTN and SPR were 27.58 cm SL and 419.39 g BW and 29.71 cm SL and 457.28 g BW, respectively (Fig. 11). The Lm50 of both males and females were estimated at 18.53 cm SL and 30.14 cm SL, respectively. By comparing among all study sites, results from logistic regressions indicated that Lm50 of males from NST matured at 18.52 cm SL, which was earlier than other sites. Female fish from PTN and SPR matured at 30.14 cm SL and 31.40 cm SL, respectively (Fig. 12).
This study provides an intensive aspect on the reproductive traits of E. tetradactylum along the coasts of Thailand. It was observed that the majority of fish samples collected by a traditional gill net, a selective fishing gear, were larger than 20.0 cm SL with the maximum size of 48.3 cm SL (Table 4). All fishes found in this study with the size of smaller than 22.0 cm SL are male. Theoretically, the sex ratio between male and female of the population is 1:1, in exception of the organism living under a particular circumstance . Some factors that may influence sex ratio of fish are fishing, mortality, migration, environmental conditions including reproductive behavior , fish size classes , growth rate , different period of maturity [12, 55] and different position in the water column . The variations of sex ratio in different fish species from different ecosystems have been reported widely. For instances, different sex ratio was observed in Lethrinids fish (Lethrinus miniatus) , puffer fish (Carinotetraodon travancoricus)  and sawtooth barracuda, Sphyraena putnamae . On the contrary, some studies reported a greater proportion of males than females  and some highlighted females predominated males [54, 59]. Different migratory patterns between male and female were also reported to influence sex ratio [53, 54, 58]. Moreover, size at first sex change affected sex ratio of the mature fish stock structure [14, 50]. Results from this study indicated that the overall sex ratio between male and female of E. tetradactylum from PTN was 1:0.69, which is highly predominated by males especially in January, February, April, May and August (p < 0.05). This circumstance, as a greater proportion of male found in dry and SW monsoon seasons (Table 5), is probably because this fish is a protandrous hermaphrodite, which initially grows as male and later changes to female at certain age, thus more male can be found in the coastal habitat where small sized-fish dominates this ecosystem. This is supported by a report from Zamidi et al.  who indicated that males fish dominated females from gillnets in Malaysia. Another study from the rearing pond by Cheng et al.  found that the ratio of females was 34.5% in the first year of culture and changed to 90.3% after three years in the rearing pond. An equal sex ratio in monsoon season might be related to seasonal reproduction and pattern of sex change. Future study on this particular aspect is essential.
Gonadosomatic index (GSI) indicates spawning or breeding season of fish [55, 56], as an increment of GSI value posing an increasing reproductive activity . The GSI normally increases along with the maturation of fish at peak during the spawning period and declines rapidly after spawning . This study found that the GSI value of females is higher than males reflecting a greater proportion of females’ bodies retained the larger gonad [52, 53, 56]. The significant monthly variation of mean GSI was found for both male and female (p < 0.05). The trend of GSI values found in this study indicates multiple spawning behavior of this fish species. It is related to the greater proportion of mature gonasd found along with high GSI values in an exception of the period without mature females in November. Moreover, the greatest number of mature males and females were also observed in moderate rainy season (SW monsoon) and heavy rainy season (NE monsoon), indicating higher chance to breed and spawn during these periods. Thus, it can be concluded that the spawning season of E. tetradactylum in PTN, based on the GSI value, proportion of mature gonad and number of mature fish, occurs the whole year round with the highest period in moderate rainy and heavy rainy seasons. This finding supports some earlier studies reported spawning periods of E. tetradactylum from different regions. This fish species was found to have a lengthy spawning period with peak in June in Malaysian coastal waters  and peak from October to December in the north-western part of Australia . Those findings also prove that habitat and season influence spawning patterns of fish. Similar to other fish species, environmental conditions during monsoon season are suitable for breeding activity due to strong water current, high turbidity, low temperature, low sunlight [52, 55, 56, 58] and less predation . As spawning strategy and maturity phase of male and female fishes can be determined by the GSI and histological examination of gonads [61, 62], results from gonad maturation indicate the synchronous and asynchronous spawning behavior of male and female, respectively. Male can be defined as a total spawner reproducing once a year, whilst female is an asynchronous fish or a multiple spawner reproducing several times a year [42, 63, 64]. This finding supports the study conducted by Cheng et al.  who reported an asynchronous reproductive pattern of this fish from an aquaculture pond in Taiwan.
Fecundity data is able to support an exploitation of fish stock, population dynamics and fisheries management practices [52, 53]. The four-finger threadfin is a benthic-pelagic species which tends to spawn with a small number of eggs but larger in size [42, 55, 64, 65]. The average fecundity from this study is 1.84 × 105 eggs, considerably lower than those reported from other areas. The earlier studies of fecundity from Malaysia, Bay of Bengal and Taiwan were reported at 3.41 × 105—11.14 × 105 eggs , 10.05 × 105—20.91 × 105 eggs  and 3.17 × 105—25.08 × 105 eggs , respectively. However, these earlier reports investigated the bigger size classes of fish compared to the present study, thus the variation of fish size may directly affect the fecundity [56, 58]. This study also confirms that fecundity increased proportionally with an increasing body and gonad weights [12, 40, 66]. Relationships between fecundity and body characteristics were used to assess corresponding among them. More accurate assessment of relationship when gonad weight was used as indicator to correlate with fecundity, instead of length and weight of fish . This study found that egg diameter is positively correlated with gonad weight and fecundity. This pattern of relationship, as the larger fish size corresponding to the heavier gonad and leading to the bigger egg size, was highlighted by Nesarul et al. . However, Thorsen et al.  found different results for the Atlantic cod in which fecundity decreased when egg diameter increased. Moreover, the present study also observed that fishes from different sites produce different numbers of fecundities, with the highest at SPR compared to PTN and NST.
Knowledge on sex change during ontogeny is useful for fisheries management practice [22, 68]. Several factors regulate the size at sex change such as body size [22, 49, 68], predation, fishing pressure [14, 68, 69], breeding season  and parental care . This study confirms that E. tetradactylum is a protandrous hermaphrodite with the length at first sex change (Ls50) of 27.58 cm TL and weight (Lw50) of 419.39 g. This size is considered smaller compared to the earlier reports from other areas. Thus, locality may affect the size at first sex change of this fish. This is coincident with a report by Shimizu et al.  who suggested that the estimated value of Ls50 differed across the locality where fishes resided. For instance, it was 32.5 cm TL in the north-western Australia, 40.0 cm TL in the western Australia , 20.8 cm to 46.5 cm FL across northern Australia . Pember et al.  found that E. tetradactylum in eastern Queensland grew up to 100 cm TL at the age of 8 years and sex change at 40 cm TL or equivalent to 2 years. They also suggested that the first sex change of this fish took at least a year to reveal the transitional sex after spawning season. Similarly, the sex change of Labrus bergylta occurred at the end of their breeding season . An estimation of the first maturation of fish (Lm50) is essential to apply for assessing a minimum permissible mesh size and avoid over exploitation of juveniles for fisheries management [54, 60, 71]. Environmental conditions can lead to changes of growth and onset of gonadal maturation [63, 64, 72]. Stock density and food availability can also influence sexual maturity of marine fish species . Generally, female fish grow larger than male but the male matures earlier than female [52, 57] due to different energy consumption . However, this is not the case for E. tetradactylum as it starts the earlier age as male and later changes to female during ontogeny. This study indicates that male and female fishes larger than 25.78 cm SL and 31.40 cm SL, respectively, can be defined as mature. It is slightly different from other areas such as 20.11 cm TL for male from Australia , 31.10 cm FL for pooled sex from India , and 39.60 cm FL for pooled sex from Indonesia . Hence, this study provides the first report on the Lm50 values for both male and female simultaneously.
The present study illustrates reproductive aspects of E. tetradactylum in the coastal waters of Thailand. Male fish predominated the annual stock with the overall sex ratio of 1:0.69. It spawned the whole year round with peaks in moderate rainy and heavy rainy seasons. Histological examination confirms the protandrous hermaphrodite posing multiple spawning habits. The average fecundity was 1.85 × 105 ± 1.05 × 105 and positively related with standard length, body weight, gonad weight and egg diameter. The Ls50 and Ws50 for fishes from Pattani Bay and Samut Prakan province were 27.58 cm, 419.39 g and 29.71 cm, 457.28 g, respectively. The Lm50 of male from Pattani Bay and Samut Prakan province were 25.78 cm and 25.56 cm. This study thus provides crucial scientific information on the reproductive biology of E. tetradactylum which is fundamental for the management of natural fish stock and aquaculture development of this species.
Availability of data and materials
The data that support in this study are accessible on request from the corresponding author.
Fontoura NF, Braun AS, Milani PCC. Estimating size at first maturity (L50) from Gonadosomatic Index (GSI) data. Neotrop Ichthyol. 2009;7(2):217–22.
Trindade-Santos I, Freire KMF. Analysis of reproductive patterns of fishes from three large marine ecosystems. Front Mar Sci. 2015;2(38):1–10.
Mohsin AKM, Ambak MA. Marine fishes and fisheries of Malaysia and neighboring countries. Serdang: Universiti Pertanian Malaysia Press; 1996.
Carpenter KE, Niem VH. FAO species identification guide for fishery purposes. In: The living marine resources of the Western Central Pacific, Bony fishes part 3 (Menidae to Pomacentridae). Rome: Food and Agriculture Organization of the Nations; 2001. p. 2791–3379.
Motomura H. Threadfins of the world. An annotated and illustrated catalogue of polynemid species known to date. FAO species catalogue for fishery purposes. Rome: Food and Agriculture Organization of the Nations; 2004.
Motomura H, Iwatsuki Y, Kimura S, Yoshino T. Revision of the Indo-West Pacific polynemid fish genus Eleutheronema (Teleostei: Perciformes). Ichthyol Res. 2002;49(1):47–61.
Pember MB, Newman SJ, Hesp SA, Young GC, Skepper CL, Hall NG, Potter IC. Biological parameters for managing the fisheries for blue and king threadfin salmons, estuary rockcod, Malabar grouper and mangrove jack in north-western Australia. Fisheries Report and Development Corporation Report. , Western Australia: Centre for Fish and Fisheries Research Murdoch University Murdoch; 2005.
Nesarul HM, Abu Hena MK, Saifullah SM, Idris MH. Breeding biology of Eleutheronema tetradactylum (Shaw, 1804) from the Bay of Bengal Indian Ocean. World Appl Sci J. 2014;30(2):240–4.
Soe KK, Pradit S, Jaafar Z, Hajisamae S. Effects of mesh size, fishing depth and season on the catch and discards of short mackerel Rastrelliger brachysoma gillnet fishery at the mouth of Pattani Bay Thailand. Fish Sci. 2022;88:15–27.
Tirtadanu T, Chodrijah U. Some population parameters and exploitation status of fourfinger threadfin (Eleutheronema tetradactylum Shaw, 1804) in Tarakan waters, north Kalimantan. Indones Fish Res J. 2018;24(2):91–8.
Karna SK, Sahoo DK, Panda S. Length-weight relationship (LWR), growth estimation and length at maturity of Eleutheronema tetradactylus in the Chilika Lagoon India. South Asian J Exp Biol. 2012;2(3):97–102.
Zamidi I, Samat A, Zaidi CC, Mazlan AG, Gazi MA, Abul Q, Simon KD. Fecundity and temporal reproductive cycle of four finger threadfin (Eleutheronema tetradactylum) in Malaysian coastal water. Asian J Ani Vet Adv. 2012;7:1100–9.
Shihab I, Gopalakrishnan A, Vineesh N, Muktha M, Akhilesh KV, Vijayagopal P. Histological profiling of gonads depicting protandrous hermaphroditism in Eleutheronema tetradactylum. J Fish Biol. 2017;90(6):2402–11.
Ballagh AC, Welch DJ, Newman SJ, Allsop Q, Stapley JM. Stock structure of the blue threadfin (Eleutheronema tetradactylum) across northern Australia derived from life-history characteristics. Fish Res. 2012;121–122:63–72.
Newman SJ, Pember MB, Rome BM, Mitsopoulos GEA, Skepper CL. Stock structure of blue threadfin Eleutheronema tetradactylum across northern Australia as inferred from stable isotopes in sagittal otolith carbonate. Fish Manag Ecol. 2011;18:246–57.
Hena MKA, Idris MH, Wong SK, Kibria MM. Growth and survival of Indian salmon Eleutheronema tetradactylum (Shaw, 1804) in brackish water pond. J Fish Aquat Sci. 2011;6:479–84.
Zischke MT, Cribb TH, Welch DJ, Sawynok W, Lester RJG. Stock structure of blue threadfin Eleutheronema tetradactylum on the Queensland east coast, as determined by parasites and conventional tagging. J Fish Biol. 2009;75(1):156–71.
Moore BR, Stapley J, Allsop Q, Newman SJ, Ballagh A, Welch DJ, Lester RG. Stock structure of blue threadfin Eleutheronema tetradactylum across northern Australia, as indicated by parasites. J Fish Biol. 2011;78(3):923–36.
Kagwade PV. The Polynemid fishes of India. Cochin: Bulletin Central Marine Fisheries Research Institute; 1970.
Sritakon T, Nooklum R, Vibunpant S, Songkaew N, Kongchai T. Reproductive biology of fourfinger threadfin (Eleutheronema tetradactylum (Shaw, 1804)) in Nakhon Si Thammarat and Songkhla provinces. Thailand: Department of Fisheries; 2019. Technical Paper No. 23/2019.
Chopelet J, Waples RS, Mariani S. Sex change and the genetic structure of marine fish populations. Fish. 2009;10:329–43.
Molloy PP, Goodwin NB, Côté IM, Gage MJG, Reynolds JD. Predicting the effects of exploitation on male-first sex-changing fish. Anim Conserv. 2007;10:30–8.
Satapoomin U. The fishes of southwestern Thailand, the Andaman sea – a review of research and a provisional checklist of species. Phuket Mar Biol Cent Res Bull. 2011;70:29–77.
Tint KK, Ngin K, Sapari A, Souliphone K, Suwannapoom S, Viron JG, Thi V, Thanh P. Enhancing the management of the Indo-Pacific mackerel resources in the Gulf of Thailand: a synthesis. In: Fish for the People, SEAFDEC, Bangkok, Thailand. 2020;18(1):14–19.
Wattayakorn, G. Environmental Issues in the Gulf of Thailand. In: Wolanski E, editor. The environment in asia pacific harbours. Dordrecht: Springer:2006.
QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation. 2023. URL http://qgis.org.
Hajisamae S, Ibrahim S. Seasonal and spatial variations of fish trophic guilds in a shallow, semi-enclosed tropical estuarine bay. Environ Biol Fishes. 2008;82:251–64.
Hajisamae S. Trophic ecology of bottom fishes assemblage along coastal areas of Thailand. Estuar Coast Shelf Sci. 2009;82(3):503–14.
Hajisamae S, Soe KK, Pradit S, Chaiyvareesajja J, Fazrul H. Feeding habits and microplastic ingestion of short mackerel, Rastrelliger brachysoma, in a tropical estuarine environment. Environ Biol Fishes. 2022;105:289–302.
Soe KK, Pradit S, Hajisamae S. Feeding habits and seasonal trophic guilds structuring fish community in the bay mouth region of a tropical estuarine habitat. J Fish Biol. 2021;99:1430–45.
Soe KK, Hajisamae S, Petchsupa N, Jaafar Z, Fazrul H, Pradit S. Reproductive biology of short mackerel, Rastrelliger brachysoma, off Pattani Bay, Lower Gulf of Thailand. Songklanakarin J Sci Technol. 2022;44(1):103–11.
Hajisamae S, Yeesin P. Do habitat, month and environmental parameters affect shrimp assemblage in a shallow semi-enclosed tropical bay, Thailand? Raffles Bull Zool. 2014;62:107–14.
Hisam F, Hajisamae S, Ikhwanuddin M, Pradit S. Distribution pattern and habitat shifts during ontogeny of the blue swimming crab, Portunus pelagicus (Linnaeus, 1758) (Brachyura, Portunidae). Crustaceana. 2020;93(1):17–32.
Chaiwanawut C, Hattha K, Duangmala P. Patterns of rainfall in Pattani province from 1982 to 2001. Songklanakarin J Sci Technol. 2005;27(1):116–76.
Hammer O, Harper DAT, Ryan, PD. PAST: Paleontological Statistics Software Package for Education and Data Analysis. Palaeontologia Electronica (version 4.03). 2001;4(1):1–9. http://palaeo-electronica.org/2001_1/past/issue1_01.htm.
Laevastu T. Manual of methods in fisheries biology. Rome: Food and Agriculture Organization of the Nations; 1965.
Bruch RM, Dick TA, Choudhury A. A field guide for the identification of stages of gonadal development in Lake sturgeon (Acipenser fulvescens: Rafinesque): with note on Lake sturgeon reproductive biology and management implications. Fond du Lac: Sturgeon for Tomorrow; 2001.
Blazer VS. Histopathological assessment of gonadal tissue in wild fishes. Fish Physiol Biochem. 2002;26:85–101.
Peterson TS, Spitsbergen JM, Feist SW, Kent ML. Luna stain, an improved selective stain for detection of microsporidian spores in histologic sections. Dis Aquat Org. 2011;95:175–80.
Simon KD, Bakar Y, Samat A, Zaidi CC, Aziz A, Mazlan AG. Population growth, trophic level, and reproductive biology of two congeneric archer fishes (Toxotes chatareus, Hamilton 1822 and Toxotes jaculatrix, Pallas 1767) inhabiting Malaysian coastal waters. J Zhejiang Univ Sci. 2009;10(12):902–11.
Wallace RA, Selman K. Cellular and dynamic aspects of oocyte growth in teleosts. Am Zool. 1981;21(2):325–343.
Brown-Peterson NJ, Wyanski DM, Saborido-Rey F, Macewicz BJ, Lowerre-Barbieri SK. A standardized terminology for describing reproductive development in fishes. Mar Coast Fish. 2011;3:52–70.
Oliveira MR, Silva NB, Yamamoto ME, Chellappa S. Gonad development and reproduction of the ballyhoo half beak, Hemiramphus brasiliensis from the coastal waters of Rio Grande do Norte. Brazil Braz J Biol. 2015;75(2):324–30.
Sokal RR, Rohlf FJ. Introduction to biostatistics. 2nd ed. New York: Dover Press; 1981.
Nikolsky GV. The Ecology of fishes. New York: Academic Press; 1963.
R Core Team. R: A Language and Environment for Statistical Computing; Version 3.6.2.; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/. Accessed 1 May 2023.
King M. Fisheries biology, assessment and management. Oxford: Oxford University Press; 1995.
Wolfert DR. Maturity and fecundity of walleyes from the eastern and western basins of Lake Erie. J Fish Res Board Can. 1969;26:1877–88.
Gulland JA. Fish stock assessment. A manual of basic method. Rome: Food and Agriculture Organization of the Nations; 1983.
Williams AJ, Davies CR, Mapstone BD. Regional patterns in reproductive biology of Lethrinus miniatus on the Great Barrier Reef. Mar Freshw Res. 2006;57(4):403–14.
Norman GR, Streiner DL. Biostatistics: The bare essentials. Shelton: People’s Medical Publishing House; 2008.
Radhakrishnan RC, Kuttanelloor R, Balakrishna MK. Reproductive biology of the endemic cyprinid fish Hypselobarbus thomassi (Day, 1874) from Kallada River in the Western Ghats India. J Appl Ichthyol. 2020;36:604–12.
Fujimoto S, Takeda S, Yagi M, Yamahira K. Seasonal change in male reproductive investment of a fish. Environ Biol Fishes. 2021;104:107–18.
Wang SB, Hsu TL, Joung SJ. The reproductive biology of the hermaphroditic goldlined seabream Rhabdosargus sarba in the central waters of the Taiwan Strait. Fish Sci. 2020;86:793–805.
Simon KD, Bakar Y, Mazlan AG, Zaidi CC, Samat A, Arshad A, Brown-Peterson NJ. Aspects of the reproductive biology of two archer fishes Toxotes chatareus, (Hamilton 1822) and Toxotes jaculatrix (Pallas 1767). Environ Biol Fishes. 2012;93:491–503.
Roshni K, Renjithkumar CR, Kurup BM. Reproductive biology of the endemic fish Etroplus suratensis (Cichlidae) from a tropical estuary in southern India. J Ichthyol. 2021;61:460–6.
Anupama KM, Hari Sankar HS, Rithin Raj M, Harikrishnan M. Reproductive biology of malabar pufferfish Carinotetraodon travancoricus (Tetraodontidae). J Appl Ichthyol. 2019;59(4):545–54.
Ghosh S, Satishkumar M, Manas HM, Rohit P, Abdussamad EM, Gopalakrishnan A. Reproductive biology of the sawtooth barracuda, Sphyraena putnamae, in the Bay of Bengal. India Environ Biol Fishes. 2021;104:421–1432.
Rodríguez SA, Grau A, Castro-Fernández J, Castejón I, Terrados J, Morales-Nin B, Arechavala-Lopez P. Reproductive biology of pipefish Syngnathus typhle and S. abaster (Syngnathidae) from Western Mediterranean Sea. J Ichthyol. 2021;61:608–615.
Cheng AC, Liu PC, Chin TS, Wang SB, Liu CH. Growth and reproductive biology of the four finger threadfin, Eleutheronema tetradactylum (Shaw), in an aquaculture pond, Pingtung Taiwan. J Fish Society Taiwan. 2017;44(1):23–35.
Fernandes JFF, Freitas J, de Araújo SA, de Santana TC, Lobato RS, Figueiredo MB. Reproductive biology of the lane snapper, Lutjanus synagris (Linnaeus 1758) (Perciformes, Lutjanidae), in the Maranhão continental shelf Northeast of Brazil. Environ Biol Fishes. 2022;105:1033–50.
Salas-Singh C, Morales-Bojorquez E, Aguirre-Villaseñor H. Reproductive biology of the bullseye puffer Sphoeroides annulatus: Gonadosomatic index and its suitability for estimating length at maturity. J Fish Biol. 2022;101(5):1119–33.
Muchlisin ZA. A general overview on some aspects of fish reproduction. Aceh Int J Sci Technol. 2014;3(1):43–52.
Murua H, Kraus G, Saborido-Rey F, Witthames PRPR, Thorsen A, Junquera S. Procedures to estimate fecundity of marine fish species from field samples in relation to reproductive strategy. J Northwest Atl Fish Sci. 2003;33:33–54.
Pensinger LG, Brown-Peterson NJ, Green CC, Midway SR. Reproductive biology of hardhead catfish Ariopsis felis: evidence for overwintering oocytes. J Fish Biol. 2021;99(2):308–20.
Hossain MD, Rahman MM, Abdallah EM. Relationships between body size, weight, condition and fecundity of the threatened fish Puntius ticto (Hamilton, 1822) in the Ganges River Northwestern Bangladesh. Sains Malays. 2012;41(7):803–14.
Thorsen A, Witthames PR, MarteinsdÓttir G, Nash RDM, Kjesbu OS. Fecundity and growth of Atlantic cod (Gadus morhua L.) along a latitudinal Gradient. Fish Res. 2010;104:45–55.
Mariani S, Sala-Bozano M, Chopelet J, Benvenuto C. Spatial and temporal patterns of size-at-sex-change in two exploited coastal fish. Environ Biol fishes. 2013;96:535–41.
Muncaster S, Norberg B, Andersson E. Natural sex change in the temperate protogynous Ballan wrasse Labrus bergylta. J Fish Biol. 2013;82:1858–70.
Shimizu S, Endo S, Kihara S, Sunobe T. Size, age, and social control of protogynous sex change in the labrid fish Pteragogus aurigarius. Ichthyol Res. 2022;69:75–81.
Havimana L, Ohtomi J, Masuda Y, Vazquez AM. The reproductive biology of female crimson sea bream Evynnis tumifrons off the southwestern coast of Kyushu. Japan Fish Sci. 2020;86:65–75.
Silva JP, Costa MRD, Gomes ID, Araújo FG. Gonadal development and fecundity of the smooth weakfish Cynoscion leiarchus (Teleostei: Perciformes: Sciaenidae) in a tropical Brazilian bay. Zoologia (Curitiba). 2016;33(6): e20160032.
We thank the Faculty of Science and Technology, Prince of Songkla University for allowing using all laboratories and equipment and providing scholarship to Mr. Teuku Haris Iqbal to undertake his Ph. D program under this project. A special acknowledgement is given to TUYF Charitable Trust, Hong Kong to support Mr. Teuku Haris Iqbal and all research costs. This research was also supported by Prince of Songkla University under the postdoctoral fellowship program. We thank the crew of fisherman from the provinces of Pattani, Satun, Nakhon Sri Thammarat and Samut Prakan who helped in fish sampling, as well as Wasina Rungruang, Arun Lohhem and Sofiyudin Maae for field sampling and laboratory work.
Consent to participate
The grant for research costs was supported by TUYF Charitable Trust, Hong Kong with the collaboration of Faculty of Science and Technology, Prince of Songkla University for Sukree Hajisamae as principal investigator and supervisor. Postdoctoral scholarship was granted by Prince of Songkla University.
Ethics approval and consent to participate
Ethical approval was waived by the institutional animal care and use committee, Prince of Songkla University (MHESI 68014/605; Ref. Ex001/2022). All sampling and laboratory processing involving animals have been followed according to national and/or international guidelines.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Soe, K.K., Iqbal, T.H., Lim, A. et al. Reproductive characteristics of the hermaphroditic four-finger threadfin, Eleutheronema tetradactylum (Shaw, 1804), in tropical coastal waters. BMC Zool 8, 22 (2023). https://doi.org/10.1186/s40850-023-00181-w
- Gonadosomatic index
- Size at first maturity
- Size and weight at first sex change
- Spawning season