Influences of Salinity on Embryonic and Larval Development of Striped Catfish Pangasianodon hypophthalmus

Salinity intrusion in coastal areas due to climate change is alarming. In this study, the effects of salinity on embryonic and larval development of striped catfish (Pangasianodon hypophthalmus) were studied experimentally. Embryos and larvae were exposed to seven salinity treatments (0, 2, 4, 6, 8, 10, and 12 ppt), each with three replications. Considerable survivability of embryos was recorded up to 6 ppt salinity. Mortality of embryos significantly increased at 8 and 10 ppt salinity, and 100% mortality was displayed within 12 h of exposure at 12 ppt salinity. The rate of hatching was significantly reduced at 8 and 10 ppt salinity. The 24 h lethal concentration (LC50) value of salinity for embryo was 11.24 ppt. Different types of deformities, such as undeveloped yolk sac, elongated gastrula yolk sac, and yolk sac bud, were highest at 10 ppt salinity. Similar to the embryo, considerable survivability of larvae was recorded up to 6 ppt salinity, and 100% mortalities were found within 24 h of exposure at 12 ppt salinity. The 24 and 48 h LC50 values of salinity for larvae were 10.63 and 8.48 ppt, respectively. Several types of deformities, including yolk sac ulceration, spine scoliosis, tail bent, yolk sac edema, and compromised swim bladder inflation, were highest at 10 ppt salinity after 48 h of exposure. Within 24 h of exposure, about 80% yolk sac of the larvae was absorbed at 8 and 10 ppt salinity, while 30%–50% yolk sac was absorbed at 0–6 ppt salinity. Growth rates in terms of length and weight were higher at 0, 2, and 4 ppt salinity and moderate at 6 and 8 ppt salinity. Overall, the current findings define the limits to optimize hatchery procedures for the culture of this species in low saline brackish water.


INTRODUCTION
Climate changes induced salinity intrusions in the coastal areas are cause for serious concern in tropical and subtropical countries (Rabbani et al., 2013). For instance, over the last four decades, the water salinity level increased about 26% in the coastal areas of Bangladesh (Alam et al., 2017). This salinity intrusion significantly affected the fertility of coastal areas that act as a barrier in sustainable agriculture that contributes to food security (Islam et al., 2005). An increase in water salinity also changed the physicochemical parameters of the aquatic environments, which ultimately affected the aquatic flora and fauna (Ahmed, 2013). In this case, efficient management of brackish water environments is needed to overcome the climate-caused obstacle. It may be necessary to revise the culture practices of salineintolerant species as well as to find innovative techniques for other existing species.
Salinity is among the most important environmental factors that influence the physiological processes of aquatic organisms. Saravanan et al. (2018) reported that fluctuations in salinity greatly influenced the physiology of common carp, Cyprinus carpio. Salinity exposure affected the growth performance of catfish, Clarias gariepinus (Sarma et al., 2013). Fecundity and fertilization and embryonic and larval growth of fish are significantly influenced by water salinity ranges (Boeuf and Payan, 2001;Akitas et al., 2004). Although acclimation of freshwater fish to seawater is species-specific, it is also subject to variation in different developmental stages (Watanabe et al., 1985;Villegas, 1990;Suresh and Lin, 1992). For these reasons, acute salinity tolerance tests in the early stages are of value to forecast the sensitivity of freshwater species (Matern, 2001;Paavola et al., 2005). Embryos and larvae are very sensitive to any changes during all stages of fry production (Akitas et al., 2004;Shahjahan et al., 2017;Rahman et al., 2020;Ashaf-Ud-Doulah et al., 2021).
Striped catfish (P. hypophthalmus) is a popular freshwater aquaculture species. Its growth rate is high and grows well in higher density. Although this fish can survive in a wide range of environmental conditions, its tolerance for osmotic change is limited. For this reason, it should be considered essentially stenohaline in its environmental preferences (Kumar et al., 2017). This air-breathing species might be an important candidate for farming in brackish water on both small and commercial scales. In our previous study, we found better growth performance up to 8 ppt salinity of this fish species in brackish water (Jahan et al., 2019). However, it is important to examine the influence of salinity on embryonic and larval specimens to understand the likelihood of survival, growth, and physiological stability of this species throughout their life cycle in coastal aquaculture systems. A suitable salinity level for the early developmental stages and growth of fish is essential for adaptation in a new saline environment. Although adaptation to saline water is regulated by osmoregulatory mechanisms in different adult fish (Boeuf and Payan, 2001), those adaptive mechanisms may not be available during the early life stages for some fish species (Varsamos et al., 2005;Yang and Chen, 2006;Bodinier et al., 2010). Therefore, this study investigates the effect of variation in salinity on the embryonic and early larval development of striped catfish in a controlled laboratory condition. The present investigations may provide insights into the optimum salinity level for successful egg incubation and early development, which might be helpful to produce fingerlings in brackish water for the sustainable coastal aquaculture of this species.

Experimental Animal
Healthy and active striped catfish (P. hypophthalmus) broods were procured for experimentation from Bangladesh Fisheries Research Institute (BFRI), Mymensingh, Bangladesh. The average body weights of male and female brood fish were 8.16 ± 0.68 and 12.34 ± 0.44 kg, respectively. The experiment was conducted in compliance with relevant guidelines followed by the Animal and Ethical Committee, Bangladesh Agricultural University (Assent No.: BAU-FoF/2020/004). The legitimate authority of the university approved the study including all the methodologies and fishes for their use for scientific purposes.

Induced Breeding and Fertilization of Eggs
For induced breeding, two male and two female brood fish were randomly selected. Male and female fish were differentiated by observation of slender genital papillae and swollen abdomens with a reddish color, respectively. They were placed in separate tanks for 6 h with a mild flow of tap water to induce spawning. Both male and female fish were given intramuscular injection of carp pituitary gland (PG) hormone at dosages of 2 and 6 mg/kg of body weight, respectively. After injection, both male and female fish were placed in a spawning tank for 8 h. The fish were placed one by one in a plastic bucket for stripping. Male and female were stripped by applying gentle pressure on the abdomen to collect milt and eggs in a plastic bowl. Milt and eggs were mixed thoroughly by using clean and soft poultry feather to accomplish fertilization in well-aerated freshwater.

Experimental Scheme
Fertilized eggs and larvae were exposed to seven salinity treatments (0, 2, 4, 6, 8, 10, and 12 ppt), each with three replications. Natural brine solution containing 260 ppt of salinity was obtained from Marine Fisheries and Technology Station, Cox's Bazar, Bangladesh. To prepare desired concentration of saline water, both brine and dechlorinated tap water were filtered through 65 µm of filter paper to remove debris and other suspended solid particles. The different concentrations of test solution were prepared at 0, 2, 4, 6, 8, 10, and 12 ppt salinity by adding the required amount of dechlorinated tap water to dilute the stock solution.

Effects of Salinity on Embryonic Stages of Striped Catfish
Immediately after fertilization, the fertilized eggs were released into previously prepared salinity (0, 2, 4, 6, 8, 10, and 12 ppt)controlled 20 L aquaria. Eggs were quantified in order to stock the same number into each aquarium, at a density of 100 eggs/20 L water. They were sampled randomly from each aquarium repeatedly from the time of stocking until hatching. Sampled eggs were preserved in 1% formalin solution in glass vials for immediate observation under a microscope. A minimum of five eggs from each treatment was observed to monitor the effects of salinity on the progress of embryonic development in each set of samples. The effects of salinity were observed using a digital microscope (MCX 100, Micros Austria). Images were obtained using a video-tracking image camera (Model-AmScope1000) connected between microscope and computer. The developmental stages of striped catfish were determined according to the developmental stages described by Puvaneswari et al. (2009).

Incubation Period and Hatchability of Striped Catfish
Distributed fertilized eggs in different salinities (0, 2, 4, 6, 8, 10, and 12 ppt) were monitored carefully every hour. Dead eggs and anomalies in eggs were recorded and removed immediately to maintain the water quality free from deterioration. After hatching, the early larvae were counted, and hatching rates were calculated in percentages using the following formula: Hatching rate = (No. of hatched fry/No. of fertilized egg) × 100.

Effects of Salinity on the Larval Stages of Striped Catfish
After hatching in normal (fresh) water, larvae were stocked into seven salinity treatments (0, 2, 4, 6, 8, 10, and 12 ppt) at a density of 50 larvae/20 L of water, each with three replications. Aeration was maintained in the aquaria throughout the experimental period. Five larval samples were collected at 12 h intervals until the completion of yolk sac absorption for microscopic examination as described above.

Determination of the Median Lethal Concentration (LC50) for Embryos and Larvae of Striped Catfish
Median lethal concentrations for both embryos and larvae were determined at 24 h intervals. The dead embryos and larvae were collected and recorded from each treatment, and the probit regression analysis was performed using IBM SPSS software version 20.0 (Chicago, IL, United States).

Effect of Salinity on Yolk Sac Absorption
During the experiment, yolk sac absorption performance was monitored in response to each of the salinity treatments. For this, the samples were taken at early hatching of striped catfish larvae and when the yolk sac was absorbed completely. Then, the yolk sac volume was measured using a digital microscope (MCX 100, Micros Austria). The yolk sac volume and yolk sac absorption were calculated according to Heming and Buddington (1988): where g is the yolk sac absorption rate, V 0 is the early volume of the yolk sac (mm 3 ), V t is the final volume of the yolk sac (mm 3 ), and t is time (h).

Effect of Salinity on O 2 Consumption Rate
Starved 1-day-old larvae (mean length = 4.5 ± 0.01 mm; mean weight = 0.029 ± 0.02 mg) were incubated and hatched at 0, 2, 4, 6, 8, 10, and 12 ppt salinity in order to estimate the oxygen consumption rates. For this, 20 larvae were randomly collected from each treatment and stocked into each salinity treatment in 200 ml glass bottles. Aeration was provided to maintain the dissolved oxygen saturation, and larvae were allowed to acclimatize the condition for 6 h. After the acclimation period, whole water was exchanged without any disturbance to the larvae, and the bottles were closed. All bottles were maintained under constant illumination and isolation to eliminate external disturbances. Bottles without larvae of each treatment were used as controls. After 48 h of incubation, larvae were sampled for the measurement of standard and total length. The larvae dry weight was taken 48 h after holding specimens at 60 • C. The oxygen concentration was determined using Winkler titration method (Strickland and Parsons, 1972): , W is the dry weight of larvae, and T is the incubation time.

Growth Performance Indices
Different growth performance indices were calculated using the following formula: 1. Yolk consumed (YC) was calculated as follows: YC (dry weight) (mg) = yolk at hatch (dry weight) -yolk at absorption (dry weight).
2. Yolk sac utilization efficiency: where, α is the long growth rate and g is the yolk sac absorption rate.
3. Length growth rate was determined using the following formula: LGR where, length f is the length final and length i is the length initial.  4. Weight growth rate was determined using the following formula: where, wt f is the weight final and wt i is the weight initial.

Data Analysis
For statistical analysis, all data were collected and noted carefully in the Microsoft Excel Worksheet 2010 during the whole experimental period. All analyzed data were exhibited as mean ± standard deviation (SD). One-way analysis of variance (ANOVA) was used followed by Tukey's test to analyze data for checking statistical significant difference among different salinity groups. The significant level of p < 0.05 value (SPSS version 20.0 for windows SPSS Inc., Chicago, IL, United States) was used to perform all the statistical analyses.

Effect of Salinity on Embryonic Development of Striped Catfish
Considerable survivability of embryos was recorded at up to 6 ppt salinity. Mortality of embryos significantly (p < 0.05) increased at 8 and 10 ppt salinity, and 100% mortality was displayed at 12 ppt salinity within 12 h of exposure (Table 1). Similarly, the rate of hatching success was significantly reduced (p < 0.05) at 8 and 10 ppt salinity ( Table 1). The incubation period depends on the salinity and incubation duration increased with increasing salinity ( Table 1). The 24 h LC50 value of salinity for the embryo was 11.24 ppt. Different types of deformities, such as damaged cell, undeveloped yolk sac, coagulated embryo in glyphosate, elongated gastrula yolk sac, streaks within the chorion, and yolk sac bud, were recorded in embryos at different time intervals in response to different salinities (Figure 1). The highest rate of deformities of embryos was recorded at 10 ppt salinity ( Table 2).

Effect of Salinity on the Larval Development of Striped Catfish
Similar to the embryo, substantial survivability of larvae was observed up to 6 ppt salinity. Larval mortality increased significantly (p < 0.05) at and above 6 ppt salinity, and 100% mortalities were found within 24 h of exposure at 12 ppt salinity ( Table 3). The 24 and 48 h LC50 values of salinity for larvae were 10.63 and 8.48 ppt, respectively. Several types of deformities, including yolk sac ulceration, the bent axis in either direction, impairment of swim bladder inflation, pericardial edema, spine scoliosis, tail bent, and yolk sac edema, were observed in larvae at different time intervals in different salinities (Figure 2). The rate of larval abnormalities was significantly (p < 0.05) increased in response to increasing salinity after 24 and 48 h of exposure (Table 4). In response to the control condition, larvae exhibited no larval deformities. The highest percentage of deformities was recorded at 10 ppt salinity after 48 h of exposure.

Effect of Salinity on Yolk Sac Absorption, Yolk Sac Utilization Efficiency, and O 2 Consumption of Larvae of Striped Catfish
Yolk sac absorption, yolk sac utilization efficiency, and O 2 consumption of larvae of striped catfish exposed to different salinities are shown in Table 5. Rapidly increasing yolk sac absorption rates (%) were observed at increasing salinity. Within 24 h of exposure, about 80% of the larval yolk sac was absorbed at 8 and 10 ppt salinity, while 30-50% yolk sac was absorbed at 0-6 ppt salinity. As a result, weight at yolk sac absorption (mg) of the larvae was significantly (p < 0.05) highest at 10 ppt salinity. Accordingly, significantly (p < 0.05) higher yolk sac absorption efficiency was recorded at 0, 2, 4, and 6 ppt salinity. After 48 h of exposure, the highest larval O 2 consumption (p < 0.05) was observed at 8 and 10 ppt salinity, in comparison with controls ( Table 5).
Effects of Salinity on the Length Growth Rate, Weight Growth Rate, and Developmental Index of Larvae of Striped Catfish Length growth rate (LGR), weight growth rate (WGR), and developmental index (KD) of larvae of striped catfish exposed to different salinities are presented in Table 6. Significantly (p < 0.05) higher LGR and WGR were recorded at 0, 2, and 4 ppt salinity. Moderate LGR and WGR were observed at 6 and 8 ppt salinity, whereas significantly (p < 0.05) reduced LGR and WGR were found at 10 ppt salinity. The KD was significantly (p < 0.05) highest at 10 ppt salinity, whereas significantly (p < 0.05) lowest KD was observed in control specimens.

DISCUSSION
The salinity tolerance test among newly hatched larvae reveals an important determinant of the likelihood of successful introduction of freshwater species in estuarine environments (Boeuf and Payan, 2001;Matern, 2001;Paavola et al., 2005). This study demonstrated considerable survivability and hatching success of embryos at up to 6 ppt salinity with a 24 h LC50 value of 11.24 ppt salinity. The embryo of some strains of Nile tilapia can tolerate up to 20 ppt salinity (despite reductions in hatching success; Fridman et al., 2012), which exceeds the tolerance rate seen in this study. The 24 h LC50 value of salinity for the embryo of a freshwater aquaculture species pabda, Ompok pabda, was 15.07 ppt with lower hatching success (Alam et al., 2020), which is higher than the tolerance rate seen in this study. In another study, reduced hatching success was recorded in different salinities (0, 5, 10, 15, 20, 25, and 32 ppt) after the transfer of freshwater spawned eggs of Nile tilapia, Oreochromis niloticus (Watanabe et al., 1985). However, mortality of embryos significantly increased at 8 and 10 ppt salinity, and 100% mortality was displayed after 12 h    Frontiers in Marine Science | www.frontiersin.org of exposure to an ambient 12 ppt salinity. Higher mortalities in higher salinities might be due to changes in the osmotic and ionic influences of water on eggs and embryos (Alderdice, 1988). Limited tolerance for higher salinity is a characteristic of closely related adult pangas fish, such as Pangasius hypothalamus. This species displays higher growth rates and lower mortality rates in freshwater than in brackish water (Mandal et al., 2017). We observed various ontogenetic deformities during early development, such as damaged cells, undeveloped yolk sac, coagulated embryo in glyphosate, elongated gastrula yolk sac, streaks within the chorion, and yolk sac bud, were maximal in response to 10 ppt salinity in this study. Even slight variation from optimal salinity can induce extensive deformities, reflecting osmoregulatory challenges encountered by embryos at higher salinity (Cook et al., 2005;Brown et al., 2010;Ekrem et al., 2012;Hayes et al., 2014). In the aquatic environment, salinity has huge influences on the cellular stability of an organism. Thus, osmotic pressure was changed in case of need for a huge amount of energy.
As a result, cellular anomalies and cellular death occurred within the body of the organism (Cañedo-Argüelles et al., 2019). Higher rates of survival of larvae of striped catfish were recorded up to 6 ppt salinity, although 100% mortalities were found at 12 ppt salinity within 24 h of exposure. The 24 and 48 h LC50 values of salinity for larvae exhibited at 10.63 and 8.48 ppt salinity, respectively. Higher mortalities in higher salinity might be due to an inability to adapt to osmoregulatory challenges just after hatching and an insufficient capacity to maintain ionic and osmotic balance during the yolk sac period. Salinity plays a significant role in freshwater aquatic organisms to regulate the osmotic stability among ion concentrations and internal body fluids (Cañedo-Argüelles et al., 2019). As salinity enhances osmotic pressure, consequently, it has adverse effects on the body fitness and survivability of the organisms. The living organism in the freshwater aquatic environment requires an osmotic balance between ion concentrations and fluids within the body cells. Thus, salt ions producing toxicity have direct negative impacts  not only on physiology of the individual species but also on species population (Nielsen et al., 2003). Salinity adaptation in fully euryhaline species is reliant on the gene expression of compounds that are involved directly and indirectly in ion transport (Lee et al., 2020;Bonzi et al., 2021). A comparison of closely related medaka species revealed a relationship of cotransporter gene expression with the ability to tolerate osmotic challenges (Juo et al., 2016). Some selection favoring morphological change has been associated with the tolerance of variable osmotic conditions in the genus Fundulus (Styga et al., 2019).
However, some fish larval strains and species can maintain ionic and osmotic balance between their internal and external environments during their yolk sac stages (Alderdice, 1988;Kaneko et al., 1995), which is consistent with the survivability of the larvae exposed to salinity up to 6 ppt in this study. For instance, the larvae of the marine rabbit fish, Siganus guttatus can tolerate 14-37 ppt salinity range after transfer from 10 to 45 ppt salinity (Young and Dueñas, 1993). This species is known to be at least semi-euryhaline, inhabiting marine and estuarine waters and occasionally rivers as adults (Food Agriculture Organization [FAO], 2021). Likewise, a longer LC50 was recorded for the red snapper, Lutjanus argentimaculatus, larvae in lower salinity when transferred from 32 ppt salinity (Estudillo et al., 2000). In this study, several types of deformities including yolk sac ulceration, axis-curve or bent axis in either direction, swim bladder inflate, spine scoliosis, tail bent, and yolk sac edema were increased in response to an increase in salinity after 24 and 48 h of exposure. 6 | Effects of salinity on length growth rate (LGR), weight growth rate (WGR), and developmental index (KD) of larvae of striped catfish (n = 150 larvae).

Salinity
LGR ( In addition, ulceration of the yolk sac occurred tremendously as this is the early stage of hatching. This stage is exceedingly sensitive to environmental factors as well as toxic elements. In our previous studies, we observed similar findings that yolk sac ulceration occurred in rohu, Labeo rohita, due to the effect of temperature (Ashaf-Ud-Doulah et al., 2021) and ulceration of yolk sac in zebrafish, Danio rerio, caused by sumithion (Rahman et al., 2020). Likewise, several researchers exhibited yolk sac edema in African catfish, C. gariepinus (Marimuthu et al., 2013); yolk sac malformation in Nile tilapia, Oreochromis niloticus (El-Greisy et al., 2016); yolk sac deformities in zebrafish, Danio rerio (Zoupa and Machera, 2017) in their previous experiment. However, all of these deformities are probably due to complications faced by embryos adapting to osmoregulatory challenges in different salinities (Ekrem et al., 2012;Hayes et al., 2014). Similar responses to salinity including higher mortality, slower development, and a variety of abnormalities have been observed in the larvae of Nile tilapia (Fridman et al., 2012), catfish, C. gariepinus (Jomori et al., 2012) and pabda, O. pabda (Alam et al., 2020). Some recent investigation suggests that alterations in transepithelial potential (TEP) may serve in a signaling capacity during osmotic acclimation or recovery from injury (Kennard and Theriot, 2020). The increased survivability of larvae of striped catfish at up to 6 ppt salinity probably reflects an increase in the numbers of chloride cells and development of the branchial system, both of which help to address the osmoregulatory challenge. It has been reported that the first appearance of branchial chloride cells was observed at 8 days post-hatch in Japanese flounder, Paralichthys  (Hiroi et al., 1998). Li et al. (1995) identified chloride cells on branchial epithelia during yolk sac stages in Mozambique tilapia, Oreochromis mossambicus. These chloride cells are also observed in the membrane of the yolk sac as well as other body surfaces of embryos and larvae. However, it is assumed that the generalization of extrabranchial chloride cells poses a central crucial role in osmoregulation until the gill chloride cells to be functional (Kaneko et al., 2002). During ontogeny, a distinct shift in ionocyte from the body surface to branchial areas has been reported in the killifish, Fundulus hereroclitus (Katoh et al., 2000), the European sea bass, Dicentrarchus labrax (Varsamos et al., 2002), the Mozambique tilapia (Yanagie et al., 2009), and the Nile tilapia (Fridman et al., 2011). Salinity adaptations not only increase in the abundance of chloride cells but also require intracellular redistribution of ion transport molecules (i.e., Na+/K+-ATPase; McCormick et al., 2003). These ion transport molecules assist to regulate the osmotic steadiness and membrane potential in the cells (Geering, 2008;Kopec et al., 2014). The sodium and potassium move against the concentration gradients. The Na+ K+-ATPase pump maintains the gradient of a higher concentration of extracellular sodium and a higher level of intracellular potassium (Pivovarov et al., 2018). Normally kidney used this sodium gradient to regulate the pH and electrolyte levels in the blood, reabsorptions of amino acids and glucose, and overall filtering the waste product in the blood (El Mernissi et al., 1991). The rate of oxygen consumption has served as an indicator of metabolism in fish (Cech, 1990). In this study, O 2 consumption of larvae was found significantly highest at 8 and 10 ppt salinity. This is might be due to higher metabolic rates and additional physiological stress at non-isosmotic salinity, increasing osmoregulatory costs and oxygen demand. Oxygen consumption rates were increased with increasing salinity during the yolk sac period in O. niloticus, the Nile tilapia (Fridman et al., 2012). Similar results have been reported for juvenile and adult sea bream, Sparus sarba (Woo and Kelly, 1995) and O. mossambicus x Oreochromis hornorum hybrids (Febry and Lutz, 1987). Moreover, metabolic rates are strongly influenced by the developmental stage (Swanson, 1996). Linear relationships between age during embryogenesis and oxygen consumption have been recorded in common carp, C. carpio (Kaushik et al., 1982), striped mullet, Mugil cephalus (Walsh and Lund, 1989), and milkfish, Chanos Febri (Swanson, 1996).
In this study, within 24 h of exposure, about 80% of the larval yolk sac was absorbed at 8 and 10 ppt salinity, while 30-50% yolk sac was absorbed at 0-6 ppt salinity. Accordingly, significantly higher yolk sac absorption efficiency was recorded at 0, 2, 4, and 6 ppt salinity. This might be due to energy consumption in response to the osmoregulatory challenge just after hatching and to the need of maintaining ionic and osmotic balance during the yolk sac period. Yolk serves as the initial energy source for larvae, providing protein and free amino acids required for the early development (Dabrowski, 1983;Conceição et al., 1997). Overall, significantly higher LGR and WGR of larvae of striped catfish at 0, 2, and 4 ppt salinity, and moderate but tolerable LGR and WGR at 6 and 8 ppt salinity in this study, warrant the striped catfish as recommended candidates for coastal aquaculture.

CONCLUSION
In conclusion, for introducing a species in coastal areas, a suitable environmental salinity level is essential for the early development and growth of fish larvae. This study has been conducted and provides detailed information about the embryonic and larval development of striped catfish in different salinity concentrations. The present findings provide useful information on salinity effects on the early life development of striped catfish. The results of this study make species candidates for larval rearing in controlled conditions for restocking purposes. The results of this study have important practical implications for handling the hatchery production method properly and for the future potential for aquaculture of this species in brackish water.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by the Animal and Ethical Committee, Bangladesh Agricultural University.

AUTHOR CONTRIBUTIONS
FH performed the experiment, data collection, and drafting of the manuscript. SI and MA-U-D assisted in the data collection and analysis. MA facilitated the brood fish to get embryos and larvae. MI and CB edited the manuscript. MS conceived, designed, supervised, and analyzed the data, and edited the manuscript. All authors reviewed and approved the final manuscript.