A. K. Shakur Ahammad1*†
Neaz A. Hasan2†
Mohammad Mahfujul Haque2
Abul Bashar2†
Md. Borhan Uddin Ahmed1Mohammad Ashraful Alam3
Md. Asaduzzaman4Md. Abul Bashar3Yahia Mahmud3- 1Department of Fisheries Biology and Genetics, Bangladesh Agricultural University, Mymensingh, Bangladesh
- 2Department of Aquaculture, Bangladesh Agricultural University, Mymensingh, Bangladesh
- 3Bangladesh Fisheries Research Institute, Mymensingh, Bangladesh
- 4Department of Marine Bioresource Science, Chattogram Veterinary and Animal Sciences University, Chattogram, Bangladesh
In recent years, attaining gonadal maturation in smaller Hilsa (Tenualosa ilisha) has become a burning issue for Hilsa fishery of Bangladesh. Causes of early maturation are not yet clearly understood. Along with environmental parameters, genetic differentiation within the population was hypothesized as the main driver, and therefore, assessing the correlation between gonadosomatic index (GSI) and environmental factors and analyzing genetic diversity were set as objectives of the present study. To address these complex issues, six diverse habitats across Bangladesh were chosen for Hilsa sample collection. For GSI, gonad was dissected from fresh fish and preserved in Bouin’s fluid for histological observation. Water quality parameters such as temperature, dissolved oxygen, pH, and salinity were also assessed. 35 fish from each habitat were used to extract and amplify DNA through the PCR technique, and genetic diversity was examined. Further, to draw a firm conclusion, the phylogenetic tree of the Hilsa population was developed by the unweighted pair-group method of arithmetic mean method based on the Cyt b gene of mitochondrial DNA. Results of GSI studies revealed that peak spawning months of T. ilisha were in October and February, where October showed the highest values in all six habitats. Histological examination showed different stages of gonadal development in different sizes and ages of Hilsa. Among all sampling sites, no statistical difference was observed for GSI value; however, smaller sized and aged Hilsa being ripped were evident in Gaglajur Haor and Kali River. Among the observed water quality parameters, temperature correlated with GSI strongly. Increased GSI was observed with temperature augmentation from downstream to upper stream, irrespective of body size and age. A perplex correlation between dissolved oxygen of observed habitats and GSI was executed. Other physico-chemical parameters viz. pH and salinity exhibited weak and moderate positive association with the GSI, respectively. Haplotype diversity of mitochondrial DNA divided the Hilsa population into three possible sub-populations, where the most distant group (Gaglajur Haor and Kali River) was subjected to early gonadal maturity. Results of this study make clear conclusions regarding the role of environmental and genetic factors on early gonadal maturations, pointing fingers at the curse of climate change and anthropogenic stressors for the migration of the Hilsa fishery of Bangladesh.
Introduction
Bangladesh has one of the most extensive and active deltas, fed by three mighty rivers: the Padma, the Meghna, and the Jamuna. This contributes to a high potential for fresh and brackish-water fish, in addition to the vast marine resources (Hossain, 2001; Haldar and Islam, 2003). There are almost 293 fish species where 18 fish from Clupeiformes were reported from Bangladesh territory by Hossain et al. (2012), and Hilsa, Tenualosa ilisha (Hamilton, 1822), is the dominating species of the group. Hilsa is mainly an anadromous fish that migrates from saline water to freshwater for the purpose of breeding. During the breeding period, it comes to the river channels of the adjacent sea countries and returns to the sea after spawning. After hatching, the juvenile Hilsa remain in freshwater for a few months for their development and then return back to the sea. It owns the largest-species fishery in Bangladesh, contributing 12% (0.53 million MT) of the total fish production (Department of Fisheries, 2019).
With a wide geographical distribution of Hilsa ranging from the Persian Gulf to the Arabian Sea, this transboundary fish takes shelter in the shallow water of the Bay of Bengal (Blaber, 2000; Arai and Amalina, 2014). For breeding purposes (with a further intention of nursing the juvenile), Hilsa follow their main migration route [Meghna river (lower part)-Meghna river (upper stream)-Padma (lower stream of the Ganges)-Jamuna (New Brahmaputra River)] to reach in upper stream rivers, which facilitate them with an uncompromising favorable environment in purpose of their breeding success. Apart from this main migration route, there is also evidence that Hilsa takes entry in freshwater through some other unpopular and secondary coastal migration route (Pashur River-Madhumati River-Gorai River-Padma River; reviewed in Sarker et al., 2021). With an unknown number of mysterious issues about the migration intention of Hilsa, the fish have also been found to use unfocused water bodies, including marshy wetland ecosystems known as Haor and hill stream rivers (e.g., Someshwari River) bordering the northeastern part of Bangladesh as their travel route (Sarker et al., 2021).
Though it can be found in almost every river, estuary, and marine environment in the country (Ahmed et al., 2008), the number of small sized Hilsa with matured gonad in natural water bodies are increasing alarmingly (Almukhtar et al., 2016). The reason behind this menacing problem observed in the Hilsa of Bangladesh is still unknown. Whether these small sized Hilsa are adults, or environmental factors (e.g., climate, hydrological, and ecological changes) or genetical factors have forced the fish to develop premature gonad in the juvenile stage is a big question (Almukhtar et al., 2016). Environmental factors are very much related to climate change, and Bangladesh, a low-lying floodplain country bordering the GBM basins (Ganges, Brahmaputra, and Meghna), is considered as one of the most vulnerable countries to the consequences of climate change (Agrawala et al., 2003; World Bank, 2010). Intensive rainfall during monsoons and other natural disasters (for example, cyclones, floods, storm surges, etc.) becomes a regular incidence in Bangladesh due to geographical south attachment with Bay of Bengal and north attachment with Himalayas (Ferdous and Baten, 2011). With a sharp increase of 0.103°C temperature per decade of the country in the last four decades (Shahid, 2010), it has been projected that the temperature will continue to increase by 1°C by 2030, 1.4°C by 2050, and 2.4°C by 2100 (IPCC, 2007). Rising temperature has also been altered by the seasonal variability in rains with their spatial distribution (Shahid and Khairulmaini, 2009) by augmenting air moisture-holding capacity.
Usually, the old literature indicates that Hilsa stock of the Indian subcontinent became mature between 16 and 17 cm in case of males and for females the range was 19–20 cm, when their age was 2+ or 3+ (Pillay, 1958). However, in an alteration of seasonal pattern with a combined effect of vigorous nourishment during their migratory journey, the age range has diverged for the sexual maturity of Hilsa, which now can be appear at different sizes and ages of fish (Zhang et al., 2009; Mohamed and Qasim, 2014). Furthermore, genetic factors, feed competition, or their persistence of stay in the river may stimulate early gonadal development. It may also be possible for some small Hilsa to reside in the upper stream instead of returning to the downstream and thus, attain gonadal maturity in freshwater (Ahmed et al., 2018). However, these hypotheses need intensive research for confirmation. Moreover, determining the gonadosomatic index (GSI), observation on gonadal maturation, genetic diversity of the selected stocks as well as assessment of the water quality parameters of the habitat of Hilsa are reasonably necessary for the validation of these hypotheses. Therefore, it was the center of the attention of the present study.
Reproduction is a continuous developmental process throughout ontogeny, requiring energetic, ecological, anatomical, biochemical, and endocrinological adaptations (Caputo et al., 2000). The GSI is a good indicator of reproductive activity being used to determine gonadal maturation stages (Hismayasari et al., 2015). It represents details regarding the reproductive status and determining reproduction time of fish (Mohan and Jhajhria, 2001; Shankar and Kulkarni, 2005). GSI also provides data on seasonal variations in multiple spawning and variation of food consumption (Pradhan, 1965) and facilitates the effective conservation of any species (Ghaedi et al., 2013; Rahman, 2014). Zhang et al. (2009) confirmed that the spawning activities of fish can be determined based on the monthly variations of GSI, and it is proportionate to each maturity stage. Bladon et al. (2019) in their study, confirmed that a significant drop of GSI from high to low is a signal of egg depositing activity in fish and that is why the GSI is a convenient method to estimate the spawning period in fish. Variations in the GSI are related to the changes in the water quality parameters and food abundances, reflecting the size variation of Hilsa (Hasan et al., 2016). In addition, water quality parameters control the maturation of gonads in fish (Hazard and Eddy, 1951). Therefore, fish nutritional status and gonad maturation are a result of their relationship with environmental factors (Ribeiro-Silva et al., 2018). Once again, histological study on the gonadal changes of fish is an essential tool for identifying reproduction level, which includes different developmental phases. Unver and Saraydın (2012) found a strong correlation between histological change in gonad and the sexual maturation of Hilsa. Therefore, a combined study of GSI and gonad histology can precisely determine the maturation timing of any fish species. Study of vitellogenic stage (VG) of oocyte and increased gonad weight (GW) and histology of gonad are the important aspects of reproductive biology that indicate the maturation level and can estimate approximate spawning time (Brown-Peterson et al., 2011). Knowledge on the fish size and age at first sexual maturity is essential to estimate the spawning stock’s size (Neja, 1992). Hence, size and age at first sexual maturity are closely linked to lifetime fecundity of fish (Stearns, 1992; Bernardo, 1993). In this regard, Ahmed et al. (2020) initiated a good attempt to determine the size and age at first sexual maturity of T. ilisha by otolith examination and confirmed that Hilsa of 1.2 years with a total length (TL) of 26 cm is sexually matured, which is a clear indication of early gonadal maturation.
Water quality parameters such as water temperature, pH, and oxygen are the crucial drivers of the gonad development for several fish species in the tropical region (Bailly et al., 2008). For instance, the water temperature can trigger the reproductive processes of many fish species and, therefore, their gonadal development may be connected to rising sea water temperature (Mondal, 2012). Mylonas and Zohar (2007) reported that a change in water temperature had a high correlation with the GW of fish and GSI and sexual maturation. Seasonal variations of physico-chemical factors of water also directly influence the gonad maturation of fish since they affect the availability of food and trigger a number of neuro-endocrinological processes in fish (Araujo et al., 2019).
Roomiani et al. (2014) again reconfirmed that the differences in the spawning maturation of fish in different areas may be due to genetical and environmental factors like temperature, pH, DO, etc. After maturation for spawning purpose Hilsa shad follow two types of migration pattern- (1) fluvial potamodromous (fish remain in freshwater habitats and breed therein) and (2) marine ecotype (fish migrate to salt water for breeding purposes) in Bangladesh (Raja, 1985). Based on these, different levels of genetic variations and diversities have been recognized by previous studies (Majumder and Alam, 2009; Asaduzzaman et al., 2019, 2020; Sarker et al., 2021). However, interpretation on whether these variations are caused by different migratory patterns or is the cause of reciprocal migration is not yet understood. The genetic architecture of vertebrates, including fish, is a ubiquitous and potent caliber that shapes ovarian maturation, leading to defining the next generation (Barson et al., 2015; Mobley et al., 2021). It is well established that the gonadal maturation of migratory Atlantic Salmon is strongly influenced by genetic factors affecting the span of different life stages (Healy et al., 2019). However, the relationship between genetic diversity and early gonadal maturity of Hilsa shad remains unboxed. Due to retaining maternal inheritance throughout the generation, preventing recombination, and showing a better alteration rate, mitochondrial DNA is considered more appropriate to measure genetic deviation than nuclear DNA (Majumder and Alam, 2009). To represent genetic differentiation among the six different habitats and find an underlying relationship with gonadal maturation, the Cyt b gene of mitochondrial DNA were analyzed to calculate haplotypes and genetic diversity.
In compliance with the above research, a better understanding of early gonadal development in T. ilisha during the natural reproductive cycle is crucial to determine the actual cause of early gonadal maturation of Hilsa in Bangladesh. Therefore, the present study was aimed to (1) determine the seasonal variation of GSI, (2) assess the gonad maturity of T. ilisha across six diverse habitats of Bangladesh through histological observation, (3) determine the variation in the water quality parameters in those habitats during the reproductive period, and (4) analyze the genetic diversity to identify whether the selected Hilsa stock maintained single or multiple sub-populations.
Materials and Methods
Sample Collection, Age, and Growth Determination of T. ilisha
Fresh Hilsa were collected randomly from six different locations of Bangladesh, namely the Meghna River Estuary, Chandpur; the Bay of Bengal, Cox’s Bazar; the Kali River, Kishoreganj; the Tetulia River, Barisal; the Padma River, Mawa, Munshiganj; and the Gaglajur Haor, Mohanganj, Netrokona (Figure 1). Fish were collected with the help of local fishers by using fishing boats from January to December, 2018. All samples were carried to the Fisheries Biology and Genetics Laboratory of Bangladesh Agricultural University (BAU), Mymensingh. Thirty samples from each habitat were collected for analysis. TL were taken to the nearest mm using a measuring scale attached to a board. The individual whole fish was measured as body weight (BW) using an electronic digital balance (MONOBIOC, B2002-S). The age and growth pattern of different sized Hilsa were determined based on the lunar ring in the otolith and length-weight relationships, respectively (Ahmed et al., 2020). Thereafter, the samples which were under normal maturation were used for the present study. To visualize the rings in the otolith, they were ground and polished using a combination of the acid contact method (Sinha and Jones, 1967) and the wet-dry sandpaper method (Ahmed et al., 2020). The use of too much acid was avoided to hinder contact with acid to the edge of the otolith. During this grinding procedure, the otolith was sporadically examined under a microscope to check for the appearance of the rings. The growth pattern of Hilsa was determined by the length weight-relationship for each of the four habitats. The mean values of TL, weight, age, and their growth pattern were recorded for tabulation.
Figure 1. Graphical picture of six sampling locations of Hilsa for the experiment. (A) the Gaglajur Haor, Mohanganj, and Netrokona; (B) the Kali River, Kishoreganj; (C) the Padma River, Mawa, and Munshiganj; (D) the Meghna River Estuary, Chandpur; (E) the Bay of Bengal, Cox’s Bazar; and (F) the Tetulia River, Barisal.
Determination of Gonadosomatic Index
The fish’s abdomens were opened by cutting the abdomen with a sharp scissor, and the muscle, fat tissue, digestive organs, blood veins, etc. were appropriately removed. Both left and right gonads were measured together as GW. The collected gonads were kept in Bouins fixative about 8–10 h. The ratio of gonad and fixative was 1:10. After 8–10 h, gonads were rinsed with 10% formalin. At the end, these gonads were preserved with 10% formalin in small vials. The GSI of each fish specimen was calculated by the following formula (Howaida et al., 1998):
Where, GSI = Gonadosomatic Index, GW = Gonad Weight, and BW = Body Weight.
Relationship Between Body Weight, Total Length, Age, and Gonadosomatic Index
For establishing the relationship of GSI with different body parameters viz: TL, BW, and age, all data were recorded with sophisticated instruments. TL and BW were measured as described in section “Sample Collection, Age and Growth Determination of T. ilisha” in this manuscript. Age of each fish were counted from their otolith. GSI were calculated by the method as described in section “Determination of Gonadosomatic Index.” All data were processed in Microsoft excel to generate a comparative graph between GSI and TL, BW, age of Hilsa.
Histological Observation of the Ovary
The histological observation was performed with several steps such as dehydration, cleaning, infiltration, embedding, trimming, sectioning, staining, and mounting. The dehydration, clearing, and infiltration of the sample were done using an automatic tissue processor in a series of increasing alcohol solutions, according to Johnston et al. (2009). After that, two changes of xylene and three series of changes in paraffin wax were completed in the processor. Finally, samples were left overnight in the laboratory embedded in and covered by paraffin mass, which eventually solidified and supported the segment during sectioning. Paraffin-embedded blocks were cut by a microtome knife at 4–5 μm size and left the sections into a water bath at a temperature of 37°C as described by Ahammad et al. (2021). the side of the paraffin blocks were trimmed by a sharp knife to obtain a suitable cuboidal shape in which the undesirable wax (paraffin) layer of the embedded blocks was removed. Sectioning was done by using a manual microtome machine (Biobase, Bk-MT260M). After obtaining sections in ribbons, it was correctly placed in a hot water bath having hot water of 50°C. At the end, the slide was stained with xylene, hematoxylin followed by periodic acid Schiff reaction. All images from histological observation were processed following the method of Moallem et al. (2015), using a high magnification microscope (CX 41, Olympus) attached with a high-resolution camera.
Observation of the Water Quality Parameters
Water samples were collected in glass bottles throughout the year from six selected habitats of Hilsa. Samples were triplicated from each location for minimizing standard error. After sampling, the bottles were marked with the respective identifying numbers. Water quality parameters such as temperature, dissolved oxygen, and pH were recorded immediately, and salinity was measured in the laboratory. Water temperature (°C) from each habitat was recorded by using a Celsius thermometer (SMART SENSOR∗ AR867, China). Dissolved oxygen (mg/l) and pH were measured using a digital DO meter (HANNA, Model- HI 9146, Romania) and a direct reading digital pH meter (Milwaukee PH meter, Model-PH55/PH56, United States), respectively. Salinity was measured by a Hand Refractometer (ERMA Hand Refractometer, Model-MASTER-S10a, Mumbai, India).
Determination of Genetic Diversity of Selected Hilsa Populations
Genomic DNA was extracted from the fin samples of Hilsa (N = 210; 35 samples/populations) using GeneJET Genomic DNA Purification Kit (Thermo Fisher Scientific, United Kingdom). Next, DNA was measured by GENEVA NANO micro-volume spectrophotometer, and amplification was performed by PCR using the universal primer set for the Cyt b gene of the mitochondrial DNA. Then amplified PCR products were electrophoresed and visualized on 1.0% agarose gels containing 2 μl ethidium bromide. The PCR products were purified according to FavorPrep GEL/PCR Purification Kit (Favorgen Biotech Corp) protocol and sequenced bidirectionally using an ABI 3730 Genetic Analyzer following the manufacture’s protocol. The raw nucleotide sequences of all the specimens were trimmed and aligned among themselves using the ClustalW software in MEGA 7.0 (Kumar et al., 2018). Haplotypes were calculated using the Kimura-2-Parameter model (Kimura, 1980). The complete aligned dataset was calculated for polymorphic sites (Ps), the number of haplotypes (H), haplotype diversity, and nucleotide diversity for each population using DnaSP 5.10 (Librado and Rozas, 2009). We estimated Tajima’s D (Tajima, 1989) and Fu’s Fs (Fu, 1997) to detect signatures of population bottlenecks or expansions and deviations from the pattern of polymorphism expected from the neutral model of evolution. The deviation from neutrality assessed by both D and Fs indices was determined in 1,000 simulations. The observed and expected distributions of the number of pair-wise mutational differences between haplotypes were calculated to their relative frequencies (mismatch distributions) under the sudden expansion model (Rogers and Harpending, 1992).
Statistical Analysis
Growth and GSI data of Hilsa were analyzed using MS Excel 2016. Comparison between TL, weight, age and GSI were established using MS Excel 2016. Significant difference of GSI between the habitats was done by one-way analysis of variance followed by duncan’s multiple range test at 95% level of confidence. Scatter plot representing the water quality parameters against GSI of T. ilisha was attained from the linear model. The Spearman rank correlation test was executed to measure the linear or monotonic association (Al-jabery et al., 2020) between water quality parameters and GSI of T. ilisha. The statistical tests were two-sided, and statistical significance was considered at P < 0.05. Statistical analyses were performed in R (RStudio Team, 2018). The genetic distance values were used to construct the unweighted pair-group method of arithmetic mean (UPGMA) dendrogram using the software POPGENE version 1.32 and MEGA version 10.0.5 (2018).
Results
Age, Size, and Growth Pattern
Mean TL was found highest in the case of Hilsa from the Tetulia River. However, the highest standard deviation (in the case of mean weight) indicates that the size variation was very high in this habitat, and the fish were not in uniform size. On the other hand, the TL and weight was lowest in the case of fish from the Gaglajur Haor. Mean age was highest in the fish from the Tetulia River, which was on average 41.8 months. The length-weight relationship indicates a positive allometric growth for the Meghna River, the Bay of Bengal, the Tetulia River, and the Padma River. On the other hand, the Kali River and the Gaglajur Haor showed isometric and negative allometric growth, respectively (Table 1).
Table 1. Age and growth parameter of Hilsa sample collected from six different habitats.
Gonadosomatic Index During Spawning Season
In all the habitats, GSI values were found apparently similar (P > 0.05; Supplementary Table 1). Hilsa from the Meghna River, the Padma River, the Tetulia River, the Bay of Bengal, the Gaglajur Haor, and the Kali River showed GSI values ranged from 4.9–11.6, 4–12.4, 4.9–13.3, 4.7–12.4, 4.5–11.2, and, 4.7–11.6, respectively, recorded from January to December. The highest GSI values were recorded in October, and the lowest in April. However, in the case of the Bay of Bengal the lowest GSI was recorded in June. Mean GSI of the Tetulia River showed higher values in August (8.4) and October (13.3) than fish from other habitats in the same months. Clearly, the inconsistencies for the Gaglajur Haor and the Kali River were visible because they both showed lower curves than the other four habitats. It indicates that fish from these two habitats consist of more negligible differences between their body and GW (P > 0.05). The mean GSI value is higher in February and October in all the habitats indicating these 2 months as peak GSI and thereby the spawning time for Hilsa (Figure 2). Nevertheless, between these 2 months, GSI in October was higher than February, indicating October as the major spawning season for Hilsa in all six habitats.
Figure 2. Monthly representation of GSI of Hilsa collected from six different habitats.
Length-Weight, Age, and Gonadosomatic Index Relationship
Body weight, TL, and age with respect to GSI showed a regular and uniform pattern in all Hilsa samples from six habitats (P > 0.05). However, in the case of Hilsa from the Kali River and Gaglajur Haor, during the time of peak spawning season, the BW was the lowest (Figure 3). On the other hand, the BW of Hilsa from other habitats showed the highest value in different seasons of the year. In relation to GSI, all graphs showed a proportional relationship with TL, BW, and age as these three parameters increased with GSI in all six samples. In all the collected samples, fish were older and held the highest GSI in October compared to other months of the year.
Figure 3. Monthly variation of GSI in respect to the body weight of Hilsa in six different habitats.
Ovarian Maturity
In all the six habitats, at the first stage, oocytes were spherical, oval, and multi-faceted, and the large nucleus occupies most of the cell. Prominent nucleoli were found sporadically in the nucleoplasm. The cytoplasm was strongly basophilic and stained with hematoxylin (Figure 4A). Nucleolus and Yolk nucleus were clearly visible in the pronucleus stage (Figure 4B). Small alveoli were formed around the nucleus, which increased in number and size during the cortical alveoli stage (Figure 4C). In primary vitellogenesis stage, vitellin granules around the nucleus clearly appeared (Figure 4D). Fat vacuoles were seen in large numbers around the pore. The number of granules increased and moved to the center of the cell in secondary vitellogenesis (Figure 4E). Cell layers were thickened, and zona radiata was stained with eosin heavily (Figure 4F). In the tertiary vitellogenesis stage, oocytes were seen with an unaided eye (Figure 4G), and they were more spherical in shape (Figure 4H).
Figure 4. Different stages of Oocyte development in Tenualosa ilisha (A–H). (A) Section of immature ovary of T. ilisha. Stage of chromatin nucleus (green arrow). Stage of pronucleues (yellow arrow), Nucleus (N), Nucleolus (n), Yolk nucleolus (ny), and Oogonia (og; H & Ex 10); (B) Oocyte in pronucleues stage, Nucleus (N), Nucleolus (n), and Yolk nucleolus (ny; H & Ex 40); (C) Oocyte in cortical alveolus stage. Cortical alveolus stage (Yellow arrow), Cortical alveolus (ca), Pronucleus (Blue arrow), and Oogonia (og). (H & Ex 4); (D) Oocyte in primary vitellogenesis stage. An Oocyte in primary vitellogenesis (Yellow arrow with germinal vesicle in around), pronucleus (Blue arrow), cortical alveolus (ca), Nucleus (N), and Yolk globule (g; H & Ex 40); (E) Oocyte in secondary vitellogenesis, (H & Ex 4); (F) The magnified view of boxed areas of panel (E), Theca (Th), granulosa (Gr), and Zona radiate (Zr; H & Ex 40); (G) Oocyte in tertiary vitellogenesis stage, (H & Ex 4); and (H) The magnified view of boxed areas of panel (G) where gray arrow indicate the tertiary vitellogenesis stage, (H & Ex 40).
In the comparison of ovarian development with age from the previous study by Ahmed et al. (2020), Hilsa from the Meghna River, the Bay of Bengal, the Padma River, and the Tetulia River showed regular ovarian development in respect to their age. No fish were found with a fully matured ovary at the age below 2 years. However, in the case of the Kali River and the Gaglajur Haor, irregular ovarian development was recorded, and fully developed gonads were confirmed by identifying the tertiary vitellogenesis stage in some small size Hilsa with an age of 14 months.
Continuation of the Correlation Between the Gonadosomatic Index and Water Quality Parameters
Temperature
The uphill pattern of the regression line with R > 0.476 in habitats except Gaglajur Haor indicates a moderate positive relationship between temperature and GSI of T. ilisha. The relation between temperature and GSI of T. ilisha in Gaglajur Haor was found to have a strong positive relationship (R > 0.93; Figure 5). In all habitats, p-values were lower than 0.05, indicating that there is a very strong correlation between the temperature and gonadal maturity of T. ilisha.
Figure 5. Scatter plot and linear regression for the year-round means of Gonadosomatic index (GSI) of Tenualosa ilisha against mean temperatures for six diverse habitats across Bangladesh: (A) Kali river, (B) Gaglajur Haor, (C) Padma river, (D) Meghna river, (E) Tetulia river, and (F) Bay of Bengal. Pearson correlation coefficient is reported as “R” while probability value is reported as “P.”
Dissolved Oxygen
The downhill pattern of the regression line with R < −0.968 in four habitats from downstream except for the Meghna River indicates a negative relationship between dissolved oxygen and GSI of T. ilisha. The negative relationship found in the Padma River was significant (P < 0.05) and in the other two habitats, the Tetulia River, and the Bay of Bengal, the negative association was non-significant (P > 0.05; Figure 6). For the Meghna River and other two habitats from upper stream (the Kali River and the Gaglajur Haor), the uphill pattern of the regression line with a R < 0.476 indicates a moderate positive relationship between dissolved oxygen and GSI of T. ilisha. In the case of a positive relationship, the Meghna River was found to have a significant impression (P < 0.05), where in both habitats from the upper stream, a positive but non-significant relationship was found (P > 0.05; Figure 6).
Figure 6. Scatter plot and linear regression for the year-round means of Gonadosomatic index (GSI) of Tenualosa ilisha against mean dissolved oxygen (DO) for six diverse habitats across Bangladesh: (A) Kali river, (B) Gaglajur Haor, (C) Padma river, (D) Meghna river, (E) Tetulia river, and (F) Bay of Bengal. Pearson correlation coefficient is reported as “R” while probability value is reported as “P.”
pH
Among the four downstream habitats, the uphill pattern of the regression line with R < 0.476 in the Padma River, the Meghna River, and the Bay of Bengal apprise the moderate positive relation within pH and GSI of T. ilisha. The relationship between pH and GSI for the Tetulia River and other observed habitats downstream exhibited very weak positive relationships for having R∼0.07. And these positive relationships for all four habitats downstream were found to have non-significant impressions (P > 0.05). The relationship between GSI of T. ilisha and pH for Kali River was significantly moderate negative (R = −0.67, P < 0.05). In another habitat, Gaglajur Haor, from the upper stream, the relationship was very weakly negative (R∼−0.13) and non-significant (P ≥ 0.05; Figure 7).
Figure 7. Scatter plot and linear regression for the year-round means of Gonadosomatic index (GSI) of Tenualosa ilisha against mean pH for six diverse habitats across Bangladesh: (A) Kali river, (B) Gaglajur Haor, (C) Padma river, (D) Meghna river, (E) Tetulia river, and (F) Bay of Bengal. Pearson correlation coefficient is reported as “R” while probability value is reported as “P.”
Salinity
The relationship between GSI of T. ilisha and salinity in two saline habitats (Tetulia River and Bay of Bengal) from downstream was found to be moderately positive (R < 0.476) and non-significant (P > 0.05; Figure 8).
Figure 8. Scatter plot and linear regression for the year-round means of Gonadosomatic index (GSI) of Tenualosa ilisha against mean salinity for two diverse habitats across Bangladesh: (A) Tetulia river, (B) Bay of Bengal. For other four habitats (Kali river, Gaglajur Haor, Padma river, and Meghna river) scatter plot and linear regression was avoided for their observed data recorded of zero salinity. Pearson correlation coefficient is reported as “R” while probability value is reported as “P.”
Spearman Rank Test
Data regarding the GSI of T. ilisha and water quality parameters of all observed habitats were included in the Spearman rank correlation test to obtain their specific correlation coefficient (SCC) values. In the case of the Spearman rank test, the habitats were categorized into upper stream (Kali River and Gaglajur Haor) and downstream (Padma River, Meghna River, Tetulia River, and Bay of Bengal) and their average values were considered for both GSI and all water quality parameters. Among the average water quality parameters of the upper stream, only the average temperature showed a strong correlation (SCC: 0.615), with a very high statistical significance (P < 0.05) to GSI of upper stream T. ilisha. For the same pair of variables (average GSI and average temperature), the same strong (SCC: 0.581, P < 0.05) correlations were observed in T. ilisha downstream (Table 2). In the same dataset, the average temperature was found to have strong but negative (SCC: −0.628, P < 0.05) correlation with moderate dissolved oxygen content for upper stream. Moreover, the average temperature downstream had a strong correlation (SCC: 0.693, P < 0.05) with average pH (Table 2).
Table 2. Correlation between the upper stream and downstream water quality parameters with GSI of T. ilisha.
Genetic Diversity of Selected Populations of Hilsa Shad Using mtDNA Cyt b Gene Sequence
The genetic diversity of six different populations of Hilsa shad, T. ilisha was presented in Table 3. However, the highest number of haplotypes were observed in the Bay of Bengal population (24) followed by the populations of Tetulia River (18), Meghna River (16), Padma River (14), Gaglajur Haor (11), and Kali River (9). Among the populations, all samples were showed the highest haplotype diversity (Hd = 1.000) except Gaglajur Haor (0.92105) and Kali River (0.86104). The highest number of Ps was recorded in the Hilsa shad samples of Bay of Bengal (55) followed by Tetulia River (52), Meghna River (46), Padma River (44), Gaglajur Haor (31), and Kali River (22), respectively. The nucleotide diversity (Pi) of the Bay of Bengal (0.40006) was found to be highest, followed by the Tetulia River (0.29500), Meghna River (0.04846), Padma River (0.02189), Gaglajur Haor (0.01191), and Kali River (0.01010). Still, Tajima’s D and Fu’s F tests were not significant (P > 0.05) in all the Hilsa shad populations (Table 3).
Table 3. Genetic diversity among the six different populations of hilsa shad, T. ilisha.
From the pairwise mismatch distributions of mtDNA haplotypes, we observed smooth and unimodal (bell-shaped) distribution of all Hilsa populations (Bay of Bengal, Tetulia River, Meghna River, Padma River, Gaglajur Haor, and Kali River; Figure 9). No mismatch distribution was found in the wild population.
Figure 9. Observed pairwise differences (red line) and the expected mismatch distributions under the sudden expansion model (green line) for Tenualosa ilisha of (A) Kali river, (B) Gaglajur Haor, (C) Padma river, (D) Meghna river, (E) Tetulia river, and (F) Bay of Bengal. The X-axis show the number of pairwise differences and the Y-axis shows the frequency of the pairwise comparisons.
The UPGMA dendrogram showed that the six Hilsa shad populations and wild T. ilisha species had been grouped into two major clusters, as shown in the dendrogram (Figure 10). The wild T. ilisha was in one cluster and separated by the highest genetic distance of D = 0.80. In contrast, the remaining cluster for four rivers (i.e., Tetulia River, Meghna River, Padma River, and Kali River) including the Bay of Bengal and Gaglajur Haor populations was subdivided into three clusters and separated by the genetic distance of D = 0.15.
Figure 10. Unrooted phylogenetic tree of T. ilisha produced by the UPGMA method based on mtDNACyt b region.
Discussion
Gonadosomatic index is an indicator of the state of gonadal development of fish and can be estimated during spawning season (Duarte et al., 2007; Oso et al., 2011; Amtyaz et al., 2013). Narejo et al. (2008) studied GSI for female T. ilisha from River Indus, Pakistan and recorded GSI data separately for summer (from April to September) and winter (from October to March). They observed that in the females, GSI index was higher in July and February, which is inconsistent with the results obtained in the present study. This might be due to the different geographical habitats where the studies were conducted. Flura et al. (2015) conducted a study on GSI of T. ilisha from the Meghna River, Chandpur, Bangladesh where they found the peak GSI value for females in June, October, and February. The highest mean GSI value was found in October; hence it determined October as the main spawning season of T. ilisha. In the present study, from October to February is reported as peak spawning period for Hilsa, which is supportive of the findings of Flura et al. (2015). For all six habitats, GSI values started to increase gradually from August and reached their highest peak in October and then decreased, which indicates the spent stage of the fish. Thereafter, the GSI indices increased again in the month of February, which is determined as the winter spawning period. This report of one major spawning period in October and one minor spawning period in February corroborated with the previous studies conducted by Flura et al. (2015) and Bladon et al. (2019) in the case of Hilsa. However, it is evident from the present study that no early or late peak spawning season has been occurring in these six habitats.
In all sampling sites, the maturation of the ovary was found linear with the age and size of the Hilsa. But in the case of Kali River and Gaglajur Haor, smaller sized and aged Hilsa were found with matured ovaries. In general, Hilsa gets its first sexual maturity in 31–33 cm of TL (Hossain et al., 2018). But this length of fish does not necessarily represent the age of fish in some stock (Ahmed et al., 2020). In a previous work with the same fish sample, Ahmed et al. (2020) conducted a study from the same six locations using the same fish sample as in the present study. In that study, the author reported age of Hilsa in respect to their length and weight. In the present study, histological observation was carried out using the same ovary sample from the fish used by Ahmed et al. (2020) to determine age through otolith examination where some Hilsa were found with matured ovary before reaching their usual age of sexual maturity. Except for the Kali River and Gaglajur Haor, Hilsa showed no incongruity in the case of ovarian development in respect to their size and age. Mean GSI of Kali River and Gaglajur Haor were higher than that of the rest four sampling site, indicating that ovaries of kali and Gaglajur Haor sized larger when compared to their respective body size. Data showing that smaller sized and aged fish bear larger GSI, elucidates the evidence of early gonadal maturity in these two habitats. However, ripended fish ranging smaller than 25 cm (declared as Jhatka- juvenile Hilsa by Department of Fisheries under Ministry of Fisheries and Livestock, Bangladesh) provokes alarm from considering sustainable fisheries management.
A few fish were found to have developed gonad at small size but were mature in age. This irregularity in ovarian development may occur due to their separate environments because these two-sampling sites are situated in the upper stream of the country, which may cause difference in availability of nutrients. Some environmental factors also may drive this ovarian development in Hilsa. According to Nair and Nancollas (1958) and Mondal (2012), sea water temperature rise is an important factor that may push early gonadal maturity in Hilsa. As the temperature of the Bay of Bengal is rising day by day (Kay et al., 2018), early maturity of Hilsa in this habitat is possible. This matured fish may migrate to the upper streams, the Kali River and the Gaglajur Haor. According to Hasan et al. (2016), anomaly of current velocity in the upper stream areas also influence sexual maturity in Hilsa. These variations could be attributable to genetic makeup, such as inbreeding, genetic depression, or changes to the expression levels of genes related to gonad development of Hilsa (Ahmed et al., 2020). Variation in gonad development of fish may also be possible due to climate change, food availability, chemicals, agricultural fertilizers, industrials and urban waste water, various endocrine disrubters or a combination of these factors (Yoneda and Wright, 2005; Caballero-Gallardo et al., 2016; Servili et al., 2020). It can be assumed that the Hilsa in the Kali River and Gaglajur Haor could not migrate to the sea after spawning and were somehow trapped in these water bodies. However, this phenomenon needs further investigation for confirmation by transcriptomic analysis of gonad, kidney, brain, liver, and muscle samples of different sized Hilsa.
Environmental temperature variation of aquatic habitats over the year is fundamentally important to fish for their ectothermic physiological process to which they are ambient (Neuheimer and Taggart, 2007; Buisson et al., 2008; Durance and Ormerod, 2009). Successful reproduction of T. ilisha is influenced by different extrinsic factors such as temperature, pH, salinity, and other properties of water (Jahan et al., 2017) and limits to attain successful artificial reproduction. Non-uniform water temperature across different habitats throughout the year leaves a deep impression on the early reproductive process of T. ilisha as observed in this study (Supplementary Tables 2–7). This alteration in temperature may influence the gamete maturation, ovulation, and spawning of T. ilisha like other fish (Hilder and Pankhurst, 2003; Pankhurst and King, 2010). In the case of female T. ilisha, the increased or substandard temperature of their ambient environment possibly results in conformational change of cytoskeleton dynamics like Danio rerio (Nonnis et al., 2021) and this cytoskeleton participates in oocyte maturation following ooplasm segregation (Ivanenkov et al., 1987). The maturation accelerated by an aggregation of microtubules in high temperature and halted by breakdown of microtubules in decreased temperature (Lessman, 1987; Nogales et al., 1995). These interfere in nucleus migration for active participation of actin microfilaments in cytoplasmic polarization (Becker and Hart, 1996), formation of the first polar body (Maro and Verlhac, 2002), fertilization and embryo morphogenesis (Klymkowsky and Karnovsky, 1994) entangled with final oocyte maturation. In contrast to the existing relationship between temperature and gonadal maturity found in this study, Jonsson et al. (1991) proved no influence of differing temperature on the maturity of Salmo trutta. Meanwhile, geomorphological and ecological changes in the habitats altering the reproduction pattern of Hilsa. Moreover, temperature augmentation of sea water have initiated the maturation of Hilsa (Mondal, 2012). In previous years it took about minimum 2 years to reach gonadal maturation for Hilsa, however, in recent years they are maturing within 1 year (Mondal, 2012). This might be due to the factors that are adversely influencing the upstream migration are siltation in habitats, newly constructed manmade barriers on the rivers and change in temperature variation over the years between downstream and upper stream (Mondal, 2012). This improves our understanding of the largest GSI in the smallest fish of the youngest group.
The present study reveals weak and moderate positive relationships of GSI with pH and salinity, respectively. Hossain et al. (2016) reported the 6–7.5, the optimum pH for different stages of T. ilisha; though in this study, the pH ranged from 6.4 to 9.65 in different habitats. In the case of Padma River, Meghna River, and Bay of Bengal, the pH was very much close to the mentioned optimum range by Hossain et al. (2016;Supplementary Tables 2–7). For other habitats this relationship ranged between weak positive to moderate negative possibly for the reason that high pH driven waterbodies acidification differentially affect physiology, reproduction, biochemistry, and survival of organisms (McClellan-Green et al., 2007; Ellis et al., 2014; Lane et al., 2015) of those habitats. Especially the female T. ilisha possibly becomes vulnerable to stress from alteration of their resource energy location and high production cost of eggs compared to sperm like of other fish (Marčeta et al., 2020). Moderate positive relationships between GSI and salinity of two saline habitats (Tetulia River and Bay of Bengal) of this study suggest that gonadal development of T. ilisha was promoted by salinity. In Eriocheir sinensis (Long et al., 2017) and Crassostrea gasar (Paixão et al., 2013), gonadal development was stimulated by their subservient habitats salinity. There are a number of potential explanations that could have contributed to this relationship; (1) the resemblance of osmolality of brackish water and T. ilisha hemolymph spare less energy consumption on osmoregulation (Dutta et al., 2019) make the way to use nutrient/energy for gonadal development; (2) gonadal development of T. ilisha is substantially fueled by calcium ion enriched brackish water. This fact has been substantiated from poor gonadal development of Crangon crangon (Gelin et al., 2001), and Grapsid crabs (Bas and Spivak, 2000) in low salinity.
We observed that freshwater populations of Hilsa especially Kali River, Gaglajur Haor, Meghna River, Padma River, and Tetulia River showed the lower number of haplotype, Ps, haplotype, and nucleotide diversity compared to the Bay of Bengal population (see Table 3 and Figures 9, 10), indicating low genetic diversity existing in the Hilsa stock of freshwater habitats. Majumder and Alam (2009) also observed similar results, reporting the lowest genetic variability occurred in the Meghna River compared to the Bay of Bengal population of Hilsa shad in Bangladesh, while using D-Loop region of mitochondrial DNA. Their results suggested that Hilsa shad have at least two, one marine and the other riverine-estuarine, or possibly three (riverine, estuarine, and marine), stocks of Hilsa in Bangladesh. Dahle et al. (1997) revealed three sub-populations of Hilsa shad in Bangladesh waters, such as Freshwater, Brackish water, salt water (Bay of Bengal), when using RAPD markers. Preliminary genetic differentiation is more evident in distant groups (Gaglajur Haor, Kali River, Padma River, Meghna River, and Tetulia River) when we consider the Bay of Bengal as natal habitat and can opine that migration patterns of Hilsa might influence the genetic variation. Hilsa from the Gaglajur Haor and Kali River might not follow the complete migration to the marine environment due to environmental and/or anthropogenic stressors and return to their habitat from middle streams. On the other hand, Hilsa near the coastal region (lower part of Meghna, Padma, and Tetulia rivers) might have easier access to marine habitats during the juvenile stage and return during the breeding season. As a result, while Hilsa from Padma, Meghna, and Tetulia rivers showed close genetic diversity to the Bay of Bengal, fish from Gaglajur Haor, and Kali River showed abnormal diversity, getting off from wild population. With this anomaly, the sub-population of Hilsa as distinct and the distant groups (Kali River, Gaglajur Haor, upper part of Meghna and Padma Rivers) were found to be subjected to early gonadal maturation. Unrooted phylogenetic tree based on Cyt b gene of mitochondrial DNA clearly indicates that early gonadal maturation in Hilsa shad of Bangladesh might be controlled by genetic factors. However, this incidence needs more clear investigation of upregulatory or downregulatory genes involved in early gonad maturity of Hilsa shad through transcriptomic analysis.
Conclusion
For stock assessment of any fish, the evidence-based calculation is necessary. From the aforementioned results, the obvious and caliber dominance of environmental factors, particularly temperature and genetic diversity on the early gonadal maturity of Hilsa shad can be concluded. While GSI is strongly correlated with body size and age of Hilsa shad, early gonadal maturation evident in distant population groups (Gaglajur Haor and Kali River) is linked to the haplotype diversity. Premature gonad ripening of Hilsa, a menacing problem in Bangladesh, is clearly unpacked in the present research. It will support the researchers in developing further genetic research and assist policymakers in developing more applicable strategies to strengthen production and conserve future stocks.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The animal study was reviewed and approved by the Animal Welfare and Bangladesh Agricultural University Ethics Committee (Ref. no. BAURES/ESRC/Fish/10).
Author Contributions
AA: conceptualization, methodology, investigation, writing—original draft, visualization, and fund acquisition. NH: formal analysis, data curation, writing—original draft, and visualization. MH: conceptualization, validation, writing—review, and editing. AB: methodology, data arrangement, and writing—original draft. MAh: sample collection, methodology, validation, and data analysis. MAl: sample collection, writing—review, and editing and data analysis. MAs: writing—review, and editing and data analysis. MB: funding acquisition and writing—review, and editing. YM: funding acquisition and writing—review, and editing. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
Authors would like to acknowledge Bangladesh Fisheries Research Institute (BFRI) for the financial support to carry out the research.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.758868/full#supplementary-material
References
Agrawala, S. T. O., Ahmed, A. U., Smith, J., and Aalst, M. V. (2003). Development and Climate Change in Bangladesh: Focus on Coastal Flooding and the Sunderbans. Paris: Organisation for Economic Co-operation and Development (OECD).
Ahammad, A. S., Asaduzzaman, M., Ahmed, M. B. U., Akter, S., Islam, M. S., Haque, M. M., et al. (2021). Muscle cellularity, growth performance and growth-related gene expression of juvenile climbing perch Anabas testudineus in response to different eggs incubation temperature. J. Therm. Biol. 96:102830. doi: 10.1016/j.jtherbio.2020.102830
Ahmed, M. B. U., Ahammad, A. K. S., Shahjahan, Datta, B. K., Rabbi, F., Alam, M. A., et al. (2018). Age analysis and growth patterns of Tenualosa ilisha using otolith examination and length-weight relationships from different regions of Bangladesh. bioRxiv [Preprint]. doi: 10.1101/471664 bioRxiv:471664,
Ahmed, M. B. U., Ahammad, A. K. S., Shahjahan, M., Rabbi, M. F., Alam, M. A., Sakib, M. N., et al. (2020). Age, growth and maturity of the Indian Shad, Tenualosa ilisha through otolith examination from different habitats in Bangladesh. Egypt. J. Aquat. Biol. Fish. 24, 343–359. doi: 10.21608/ejabf.2020.119368
Ahmed, M. S., Sharif, A. S. M., and Latifa, G. A. A. (2008). Age, Growth and Mortality of hilsa shad, Tenualosa ilisha in the River Meghna, Bangladesh. Asian J. Biol. Sci. 1, 69–76. doi: 10.3923/ajbs.2008.69.76
Al-jabery, K. K., Obafemi-Ajayi, T., Olbricht, G. R., and Wunsch, D. C. II (2020). “Data preprocessing,” in Computational Learning Approaches to Data Analytics in Biomedical Applications (Cambridge, MA: Elsevier Academic Press), 7–27. doi: 10.1016/B978-0-12-814482-4.00002-4
Almukhtar, M. A., Jasim, W., and Mutlak, F. (2016). Reproductive biology of hilsa shad Tenualosa ilisha (Teleostei: Clupeidae) during spawning migration in the Shatt Al Arab River and Southern Al Hammar Marsh, Basra, Iraq. J. Fish. Aquat. Sci. 11, 43–55. doi: 10.3923/jfas.2016.43.55
Amtyaz, K. M., Khan, M. Z., and Hashmi, M. U. (2013). Studies on gonadosomatic index & stages of gonadal development of striped piggy fish, Pomadasys stridens (Forsskal, 1775) (Family; Pomadasyidae) of Karachi Coast, Pakistan. J. Entomol. Zool. Study 1, 28–31.
Arai, T., and Amalina, R. (2014). New record of a tropical shad Tenualosa ilisha (Teleostei: Clupeidae) in Malaysian waters. Mar. Biodivers. Rec. 7, 1–4. doi: 10.1017/S1755267214000736
Araujo, G. S., Abessa, D. M. S., Soares, A. M. V. M., and Loureiro, S. (2019). Multi-generational exposure to Pb in two monophyletic Daphnia species: individual, function and population related endpoint. Ecotox. Environ. Saf. 173, 77–85. doi: 10.1016/j.ecoenv.2019.02.001
Asaduzzaman, M., Wahab, M., Rahman, M., Nahiduzzaman, M., Dickson, M., Igarashi, Y., et al. (2019). Fine-scale population structure and ecotypes of anadromous hilsa shad (Tenualosa ilisha) across complex aquatic ecosystems revealed by NextRAD genotyping. Sci. Rep. 9:16050. doi: 10.1038/s41598-019-52465-2
Asaduzzaman, M., Wahab, M., Rahman, M. M., Mariom, Nahiduzzaman, M., Rahman, M., et al. (2020). Morpho-genetic divergence and adaptation of anadromous hilsa shad (Tenualosa ilisha) along their heterogenic migratory habitats. Front. Mar. Sci. 7:554. doi: 10.3389/fmars.2020.00554
Bailly, D., Agostinho, A. A., and Suzuki, H. I. (2008). Influence of the flood regime on the reproduction of fish species with different reproductive strategies in the Cuiabá River, Upper Pantanal, Brazil. River Res. Appl. 24, 1218–1229. doi: 10.1002/rra.1147
Barson, N., Aykanat, T., Hindar, K., Baranski, M., Bolstad, G. H., Fiske, P., et al. (2015). Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon. Nature 528, 405–408. doi: 10.1038/nature16062
Bas, C. C., and Spivak, E. D. (2000). Effect of salinity on embryos of two Southwestern Atlantic estuarine grapsid crab species cultured in vitro. J. Crustac. Biol. 20, 647–656. doi: 10.1163/20021975-99990088
Becker, K. A., and Hart, N. H. (1996). The cortical actin cytoskeleton of unactivatedzebrafish eggs: spatial organization and distribution of filamentous actin, nonfilamentous actin, and myosin-II. Mol. Reprod. Dev. 43, 536–547. doi: 10.1002/(SICI)1098-2795(199604)43:4<536::AID-MRD17<3.0.CO;2-X
Bernardo, J. (1993). Determinants animals of maturation in trends. J. Ecol. Evol. 8, 166–173. doi: 10.1016/0169-5347(93)90142-c
Blaber, S. J. M. (2000). Tropical Estuarine Fishes: Ecology, Exploitation and Conservation. Oxford: Blackwell Science.
Bladon, A., Myint, K. T., Ei, T., Khine, M., Aye, P. T., Thwe, T. L., et al. (2019). Spawning Seasonality of Hilsa (Tenualosa ilisha) in Myanmar’s Ayeyarwady Delta. IIED Working Paper. London: IIED, 1–36.
Brown-Peterson, N. J., Wyanski, D. M., Saborido-Rey, F., Macewicz, B. J., and Lowerre-Barbieri, S. K. (2011). A standardized terminology for describing reproductive development in fishes. Mar. Coast. Fish. 3, 52–70. doi: 10.1080/19425120.2011.555724
Buisson, L., Blanc, L., and Grenouillet, G. (2008). Modelling stream fish species distribution in a river network: the relative effects of temperature versus physical factors. Ecol. Freshw. Fish 17, 244–257. doi: 10.1111/j.1600-0633.2007.00276.x
Caballero-Gallardo, K., Olivero-Verbel, J., and Freeman, J. L. (2016). Toxicogenomics to evaluate endocrine disrupting effects of environmental chemicals using the zebrafish model. Curr. Genom. 17, 515–527. doi: 10.2174/1389202917666160513105959
Caputo, V., Candi, G., and La Mesa, M. (2000). Pattern of gonad maturation and the question of semelparity in the paedomorphic goby Aphia minuta. J. Fish Biol. 58, 656–669. doi: 10.1111/j.1095-8649.2001.tb00520.x
Dahle, G., Rahman, M., and Erikson, A. G. (1997). RAPD fingerprinting used for discriminating among three populations of Hilsa shad (Tenualosa ilisha). Fish. Res. 32, 263–269. doi: 10.1016/s0165-7836(97)00058-1
Department of Fisheries (2019). National Fish Week, Compendium (In Bengali). Dhaka: Department of Fisheries, Ministry of Fisheries and Livestock, 1–124.
Duarte, S., Araujo, F. G., Sales, A., and Bazzoli, N. (2007). Morphology of gonads, maturity and spawning season of Loricarfichthys spixii (Siluriformes, Loricariidae) in a suptropical reservoir. Braz. Arch. Biol. Tech. 50, 1019–1032. doi: 10.1590/s1516-89132007000700013
Durance, I., and Ormerod, S. J. (2009). Trends in water quality and discharge confound long-term warming effects on river macroinvertebrates. Freshw. Biol. 54, 388–405. doi: 10.1111/j.1365-2427.2008.02112.x
Dutta, S., Ray, S. K., Pailan, G. H., Suresh, V. R., and Dasgupta, S. (2019). Alteration in branchial NKA and NKCC ion-transporter expression and ionocyte distribution in adult Hilsa during up-river migration. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 189, 69–80. doi: 10.1007/s00360-018-1193-y
Ellis, R. P., Spicer, J. I., Byrne, J. J., Sommer, U., Viant, M. R., White, D. A., et al. (2014). 1H NMR metabolomics reveals contrasting response by male and female mussels exposed to reduced seawater pH, increased temperature, and a pathogen. Environ. Sci. Technol. 48, 7044–7052. doi: 10.1021/es501601w
Ferdous, M. G., and Baten, M. A. (2011). Climatic variables of 50 years and their trends over Rajshahi and Rangpur division. J. Environ. Sci. Nat. Resour. 4, 147–150. doi: 10.3329/jesnr.v4i2.10165
Flura, M. Z., Rahman, B. S., Rahman, M. A., Alam, M. A., and Pramanik, M. H. (2015). Length-weight relationship and GSI of Hilsa, Tenualosa ilisha (hamilton, 1822) fishes in Meghna river, Bangladesh. Int. J. Nat. Soc. Sci. 2, 82–88.
Fu, Y. X. (1997). Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147, 915–925. doi: 10.1093/genetics/147.2.915
Gelin, A., Crivelli, A. J., Rosecchi, E., and Kerambrun, P. (2001). Can salinity changes affect reproductive success in the brown shrimp Crangon crangon? J. Crustac. Biol. 21, 905–911. doi: 10.1651/0278-0372(2001)021[0905:cscars]2.0.co;2
Ghaedi, A., Kabir, M. A., and Hashim, R. (2013). Oocyte development and fecundity of Snakehead Murrel, Channa striatus (Bloch, 1793) in Captivity. Asian Fish. Sci. 26, 39–51.
Haldar, G. C., and Islam, M. R. (2003). Hilsa Fisheries Conservation, Development and Management Technique (in Bengali). Dhaka: Department of Fisheries (DoF), 1–89.
Hasan, K. M. M., Ahmed, Z. F., Wahab, M., and Mohammed, E. Y. (2016). Food and Feeding Ecology of Hilsa (Tenualosa ilisha) in Bangladesh’s Meghna River Basin. IIED Working Paper. London: IIED, doi: 10.13140/RG.2.1.5031.2727
Hazard, T. P., and Eddy, R. E. (1951). Modification of the sexual cycle in brook trout (Salvelinus fontinalis) by control of light. Trans. Am. Fish. Soc. 80, 158–162. doi: 10.1577/1548-8659(1950)80[158:motsci]2.0.co;2
Healy, K., Ezard, T. H. G., Jones, O. R., Salguero-Gómez, R., and Buckley, Y. M. (2019). Animal life history is shaped by the pace of life and the distribution of age-specific mortality and reproduction. Nat. Ecol. Evol. 3, 1217–1224. doi: 10.1038/s41559-019-0938-7
Hilder, M. L., and Pankhurst, N. W. (2003). Evidence that temperature change cues reproductive development in the spiny damselfish, Acanthochromis polyacanthus. Environ. Biol. Fishes 66, 187–196. doi: 10.1023/A:1023601729203
Hismayasari, I. B., Marhendra, A. P. W., Saidin, S. R., and Supriyadi, S. D. (2015). Gonadosomatic index (GSI), hepatosomatic index (HSI) and proportion of oocytes stadia as an indicator of rainbowfish Melanotaenia boesemani spawning season. Int. J. Fish. Aquat. Stud. 2, 359–362.
Hossain, M. A. R., Das, I., Genevier, L., Hazra, S., Rahman, M., and Barange, M. (2018). Biology and fisheries of hilsa shad in Bay of Bengal. Sci. Total Environ. 651, 1720–1734. doi: 10.1016/j.scitotenv.2018.10.034
Hossain, M. A. R., Wahab, M. A., and Belton, B. (2012). The Checklist of the Riverine Fishes of Bangladesh. Dhaka: The world Fish Center, Bangladesh and South Asia Office.
Hossain, M. S. (2001). Biological aspects of the coastal and marine environment of Bangladesh. Ocean Coast. Manag. 44, 261–282. doi: 10.1016/s0964-5691(01)00049-7
Hossain, M. S., Sharifuzzaman, S. M., Chowdhury, S. R., and Sarker, S. (2016). Habitats across the life cycle of Hilsa shad (Tenualosa ilisha) in aquatic ecosystem of Bangladesh. Fish. Manag. Ecol. 23, 450–462. doi: 10.1111/fme.12185
Howaida, R. G., Hanlon, T. H., Mahmoud, H. H., and Eletreby, S. G. (1998). Maturation, fecundity, seasonality and reproduction of two commercially valuable cuttlefish, Sephia pharaonis and S. dollfusi in Suez Canal. J. Fish. Res. 36, 99–115. doi: 10.1016/s0165-7836(98)00107-6
IPCC (2007). “Climate change: Impacts, Adaptation and Vulnerability”, Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press.
Ivanenkov, V. V., Minin, A. A., Meshcheryakov, V. N., and Martynova, L. E. (1987). The effect of local cortical microfilament disorganization on ooplasmic segregation in the loach (Misgurnus fossilis) egg. Cell Differ. 22, 19–28. doi: 10.1016/0045-6039(87)90410-6
Jahan, I., Ahsan, D., and Farque, M. H. (2017). Fishers’ local knowledge on impact of climate change and anthropogenic interferences on hilsa fishery in South Asia: evidence from Bangladesh. Environ. Dev. Sustain. 19, 461–478. doi: 10.1007/s10668-015-9740-0
Johnston, I. A., Lee, H.-T., Macqueen, D. J., Paranthaman, K., Kawashima, C., Anwar, A., et al. (2009). Incubation temperature affects muscle fiber recruitment in adult zebrafish: genomewide changes in gene and microRNA expression associated with the transition from hyperplastic to hypertrophic growth phenotypes. J. Exp. Biol. 212, 1781–1793. doi: 10.1242/jeb.029918
Jonsson, N., Jonsson, B., and Hansen, L. P. (1991). Energetic cost of spawning in male and female Atlantic salmon (Salmo salar L.). J. Fish Biol. 39, 739–744. doi: 10.1111/j.1095-8649.1991.tb04403.x
Kay, S., Caesar, J., and Janes, T. (2018). “Marine dynamics and productivity in the Bay of Bengal,” in Ecosystem Services for Well-Being in Deltas, eds R. Nicholls, C. Hutton, W. Adger, S. Hanson, M. Rahman, and M. Salehin (Cham: Palgrave Macmillan), doi: 10.1007/978-3-319-71093-8_14
Kimura, M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16, 111–120. doi: 10.1007/BF01731581
Klymkowsky, M. W., and Karnovsky, A. (1994). Morphogenesis and the cytoskeleton: studies of the Xenopus embryo. Dev. Biol. 165, 372–384. doi: 10.1006/dbio.1994.1260
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Lane, A., Campanati, C., Dupont, S., and Thiyagarajan, V. (2015). Trans-generational responses to low pH depend on parental gender in a calcifying tubeworm. Sci. Rep. 5, 1–7. doi: 10.1038/srep10847
Lessman, C. A. (1987). Germinal vesicle migration and dissolution in Rana pipiens oocytes: effect of steroids and microtubule poisons. Cell Differ. 20, 239–251. doi: 10.1016/0045-6039(87)90469-6
Librado, P., and Rozas, J. (2009). DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452. doi: 10.1093/bioinformatics/btp187
Long, X., Wu, X., Zhao, L., Ye, H., Cheng, Y., and Zeng, C. (2017). Effects of salinity on gonadal development, osmoregulation and metabolism of adult male Chinese mitten crab, Eriocheir sinensis. PLoS One 12:e0179036. doi: 10.1371/journal.pone.0179036
Majumder, S., and Alam, M. S. (2009). Genetic Structure and stock discrimination of hilsa (Tenualosa ilisha) by RFLP analysis of mtDNA. Evol. Genet. Genet. Mol. Biol. 32:2009. doi: 10.1590/S1415-47572009005000023
Marčeta, T., Matozzo, V., Alban, S., Badocco, D., Pastore, P., and Marin, M. G. (2020). Do males and females respond differently to ocean acidification? An experimental study with the sea urchin Paracentrotus lividus. Environ. Sci. Pollut. Res. 27, 39516–39530. doi: 10.1007/s11356-020-10040-7
Maro, B., and Verlhac, M. H. (2002). Polar body formation: new rules for asymmetric divisions. Nat. Cell Biol. 4, E281–E283. doi: 10.1038/ncb1202-e281
McClellan-Green, P., Romano, J., and Oberdörster, E. (2007). Does gender really matter in contaminant exposure? A case study using invertebrate models. Environ. Res. 104, 183–191. doi: 10.1016/j.envres.2006.09.008
Moallem, Z., Abdi, R., Movahedinia, A., Shirali, S., and Salati, A. P. (2015). Gonad histology and gonadosomatic index variations during gonadal development of wild female Tenualosa ilisha. J. Persian Gulf 6, 53–58.
Mobley, K. B., Aykanat, T., Czorlich, Y., House, A., Kurko, J., Miettinen, A., et al. (2021). Maturation in Atlantic salmon (Salmo salar, Salmonidae): a synthesis of ecological, genetic, and molecular processes. Rev. Fish Biol. Fish. 31, 523–571. doi: 10.1007/s11160-021-09656-w
Mohamed, A. R. M., and Qasim, A. M. (2014). Stock assessment and management of Hilsa shad (Tenualosa ilisha) in Iraqi marine waters, northwest Arabian Gulf. Int. J. Fish. Aquat. Stud. 1, 1–7. doi: 10.21276/ijaq.2015.12.1.1
Mohan, D., and Jhajhria, A. (2001). Comparative growth and production studies of carps in Skekhawat pond Jodhpur. J. Nate 13, 123–130.
Mondal, B. K. (2012). Delicious Hilsa Turns ‘Bland’ due to Early Maturity. National Fisheries Development Board of India. Cited by Niyogi, S. 2012. Available online at: http://timesofindia.indiatimes.com/topic (accessed April 22, 2021).
Mylonas, C. C., and Zohar, Y. (2007). “Promoting oocyte maturation, ovulation and spawning in farmed fish,” in The Fish Oocyte, eds P. J. Babin, J. Cerdà, and E. Lubzens (Dordrecht: Springer), 437–474. doi: 10.1007/978-1-4020-6235-3_15
Nair, V. S. K., and Nancollas, G. H. (1958). 742. Thermodynamics of ion association. Part IV. Magnesium and zinc sulphates. J. Chem. Soc. (Resumed) 3706–3710. doi: 10.1039/jr9580003706
Narejo, N. T., Lashari, P. K., and Jafri, S. I. H. (2008). Morphometric and meristic differences between two types of Palla, Tenualosa ilisha (Hamilton) from River Indus, Pakistan. Pak. J. Zool. 40, 31–35.
Neja, Z. (1992). Maturation and fecundity of mackerel (Scomber scombrus L.) in Northwest Atlantic. Acta Ichthyol. Piscat. 22, 125–140. doi: 10.3750/aip1992.22.1.11
Neuheimer, A. B., and Taggart, C. T. (2007). The growing degree-day and fish size-at-age: the overlooked metric. Can. J. Fish. Aquat. Sci. 64, 375–385. doi: 10.1139/f07-003
Nogales, E., Medrano, F. J., Diakun, G. P., Mant, G. R., Towns-Andrews, E., and Bordas, J. (1995). The effect of temperature on the structure of vinblastine-induced polymers of purified tubulin: detection of a reversible conformational change. J. Mol. Biol. 254, 416–430. doi: 10.1006/jmbi.1995.0628
Nonnis, S., Angiulli, E., Maffioli, E., Frabetti, F., Negri, A., Cioni, C., et al. (2021). Acute environmental temperature variation affects brain protein expression, anxiety and explorative behaviour in adult zebrafish. Sci. Rep. 11:2521. doi: 10.1038/s41598-021-81804-5
Oso, J. A., Idowu, E. O., Fagbuaro, O., Olaniran, T. S., and Ayorinde, B. E. (2011). Fecundity, condition factor and gonado-somatic index of hepsetusodoe (African Pike) in a tropical reservoir, Southwest Nigeria. World J. Fish Mar. Sci. 3, 112–116.
Paixão, L., Ferreira, M. A., Nunes, Z., Fonseca-Sizo, F., and Rocha, R. (2013). Effects of salinity and rainfall on the reproductive biology of the mangrove oyster (Crassostrea gasar): implications for the collection of broodstock oysters. Aquaculture 38, 6–12. doi: 10.1016/j.aquaculture.2012.11.019
Pankhurst, N. W., and King, H. R. (2010). Temperature and salmonid reproduction: implications for aquaculture. J. Fish Biol. 76, 69–85. doi: 10.1111/j.1095-8649.2009.02484.x
Pillay, T. V. R. (1958). Biology of the Hilsa, Hilsa ilisha (Hamilton) of the river hooghly. Indian J. Fish. 5, 201–257.
Pradhan, M. J. (1965). Observation on the maturity and spawning of Psettodes erumei. Indian J. Fish. 9, 580–589.
Rahman, M. A. (2014). Habitat Ecology and Reproductive Biology of Gagora Catfish (Arius gagora) in Monuriver, Sylhet. MS Thesis. Sylhet: Sylhet Agricultural University, 71.
Raja, B. T. A. (1985). “A review of the biology and fisheries of the Marine Fisheries of Hilsa ilisha in the upper Bay of Bengal,” in Proceedings of the Marine Fisheries resourse Management. Bay of Bengal (Colombo: FAO/UNDP BOBP/WP/37), 66.
Ribeiro-Silva, R. D. C., Fiaccone, R. L., Conceição-Machado, M. E. P. D., Ruiz, A. S., Barreto, M. L., and Santana, M. L. P. (2018). Body image dissatisfaction and dietary patterns according to nutritional status in adolescents. J. Pediatr.(Rio J) 94, 155–161. doi: 10.1016/j.jpedp.2017.08.011
Rogers, A. R., and Harpending, H. (1992). Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569.
Roomiani, L., Sotudeh, A. M., and Hakimi, R. (2014). Reproductive biology of Hilsa shad (Tenualosa ilisha) in coastal Waters of the Northwest of Persian Gulf. Iran. J. Fish. Sci. 13, 201–215.
Sarker, A., Jiang, J., Naher, H., Huang, J., Sarker, K. K., Yin, G., et al. (2021). Cross-species gene enrichment revealed a single population of hilsa shad (Tenualos ailisha) with low genetic variation in Bangladesh waters. Sci. Rep. 11:11560. doi: 10.1038/s41598-021-90864-6
Servili, A., Canario, A. V. M., Mouchel, O., and Muñoz-Cueto, J. A. (2020). Climate change impacts on fish reproduction are mediated at multiple levels of the brain-pituitary-gonad axis. Gen. Comp. Endocrinol. 291:113439. doi: 10.1016/j.ygcen.2020.113439
Shahid, S. (2010). Impact of climate change on irrigation water demand of dry season boro rice in northwest Bangladesh. Clim. Chang. 105, 433–453. doi: 10.1007/s10584-010-9895-5
Shahid, S., and Khairulmaini, O. S. (2009). Spatial and temporal variability of rainfall in Bangladesh. Asia Pac. J. Atmos. Sci. 45, 375–389.
Shankar, D. S., and Kulkarni, R. S. (2005). Somatic condition of the fish, Notopterus notopterus (pallas) during phases of the reproductive cycle. J. Environ. Biol. 26, 49–53.
Sinha, V. R. P., and Jones, J. W. (1967). On the age and growth of the freshwater eel (Anguilla anguilla). J. Zool. 153, 99–117. doi: 10.1111/j.1469-7998.1967.tb05033.x
Tajima, F. (1989). Statistical methods to test for nucleotide mutation hypothesis by DNA polymorphism. Gene 123, 585–595. doi: 10.1093/genetics/123.3.585
Unver, B., and Saraydın, S. U. (2012). Histological examination of ovarium development of shemaya Chalcalburnus chalcoides living in Lake Todurge (Sivas/Turkey). Folia Zool. 53, 99–106.
Yoneda, M., and Wright, P. J. (2005). Effect of temperature and food availability on reproductive investment of first-time spawning male Atlantic cod, Gadus morhua. ICES J. Mar. Sci. 62, 1387–1393. doi: 10.1016/j.icesjms.2005.04.018
Keywords: early gonadal maturity, genetic diversity of Hilsa, small Hilsa, mitochondrial DNA, phylogenetic tree of Hilsa, GSI, Cyt b gene, Hilsa migration
Citation: Ahammad AKS, Hasan NA, Haque MM, Bashar A, Ahmed MBU, Alam MA, Asaduzzaman M, Bashar MA and Mahmud Y (2021) Environmental Factors and Genetic Diversity as Drivers of Early Gonadal Maturation: A Gonadosomatic Index Based Investigation on Indian Shad, Tenualosa ilisha Population of Bangladesh. Front. Mar. Sci. 8:758868. doi: 10.3389/fmars.2021.758868
Received: 15 August 2021; Accepted: 24 September 2021;
Published: 21 October 2021.
Edited by:
Amin Golpour Dehsari, Centre for Biology, Academy of Sciences of the Czech Republic (ASCR), CzechiaReviewed by:
Shahdat Hossain, National Institute of Biotechnology (NIB), BangladeshDavood Zargham, Iranian Fisheries Science Research Institute, Agricultural Research, Education and Extention Organization, Iran
Copyright © 2021 Ahammad, Hasan, Haque, Bashar, Ahmed, Alam, Asaduzzaman, Bashar and Mahmud. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: A. K. Shakur Ahammad, c2FoYW1tYWQwOUB5YWhvby5jb20=; c2FoYW1tYWQwOUBiYXUuZWR1LmJk
†These authors have contributed equally to this work and share first authorship