Next Article in Journal
Characterization of an Endolysin Targeting Clostridioides difficile That Affects Spore Outgrowth
Previous Article in Journal
Dysregulation of Muscle-Specific MicroRNAs as Common Pathogenic Feature Associated with Muscle Atrophy in ALS, SMA and SBMA: Evidence from Animal Models and Human Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interspecies Behavioral Variability of Medaka Fish Assessed by Comparative Phenomics

1
Department of Chemistry, Chung Yuan Christian University, Chung-Li 320314, Taiwan
2
Department of Bioscience Technology, Chung Yuan Christian University, Chung-Li 320314, Taiwan
3
Department of Applied Chemistry, National Pingtung University, Pingtung 900391, Taiwan
4
Faculty of Pharmacy and The Graduate School, University of Santo Tomas, Manila 1008, Philippines
5
Department of Aquatic Biosciences, National Chiayi University, 300 University Rd., Chiayi 600, Taiwan
6
Center for Nanotechnology, Chung Yuan Christian University, Chung-Li 320314, Taiwan
7
Research Center for Aquatic Toxicology and Pharmacology, Chung Yuan Christian University, Chung-Li 320314, Taiwan
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(11), 5686; https://doi.org/10.3390/ijms22115686
Submission received: 3 May 2021 / Revised: 23 May 2021 / Accepted: 23 May 2021 / Published: 26 May 2021
(This article belongs to the Section Molecular Biology)

Abstract

:
Recently, medaka has been used as a model organism in various research fields. However, even though it possesses several advantages over zebrafish, fewer studies were done in medaka compared to zebrafish, especially with regard to its behavior. Thus, to provide more information regarding its behavior and to demonstrate the behavioral differences between several species of medaka, we compared the behavioral performance and biomarker expression in the brain between four medaka fishes, Oryzias latipes, Oryzias dancena, Oryzias woworae, and Oryzias sinensis. We found that each medaka species explicitly exhibited different behaviors to each other, which might be related to the different basal levels of several biomarkers. Furthermore, by phenomics and genomic-based clustering, the differences between these medaka fishes were further investigated. Here, the phenomic-based clustering was based on the behavior results, while the genomic-based clustering was based on the sequence of the nd2 gene. As we expected, both clusterings showed some resemblances to each other in terms of the interspecies relationship between medaka and zebrafish. However, this similarity was not displayed by both clusterings in the medaka interspecies comparisons. Therefore, these results suggest a re-interpretation of several prior studies in comparative biology. We hope that these results contribute to the growing database of medaka fish phenotypes and provide one of the foundations for future phenomics studies of medaka fish.

1. Introduction

Medaka is a small, oviparous freshwater teleost fish distributed in East Asia that is often found in rice fields and, therefore, has been called ‘ricefish’ [1]. Medaka, especially Japanese medaka (Oryzias latipes), has been used as a model organism in basic fish biology and behavioral studies. Furthermore, in recent years, the use of medaka as a model organism has highly contributed to the knowledge in various research fields, such as genetics, toxicology, and behavior science [2,3,4,5]. Moreover, O. latipes has been proposed by the OECD as the standard fish for toxicology tests [6,7]. As an experimental organism, medaka possesses characteristics that are similar to those of the zebrafish (Danio rerio), including its small size (adult 2–4 cm in length), high fecundity, short generation time (2–3 months), simple dietary and habitat requirements, transparency of embryos, and availability of genomic information [3,8]. Thus, it serves as a complementary model to the well-established zebrafish in many fields, including cancer research [1,9]. Also, medaka holds several advantages over zebrafish, which come from its species-specific features. These features include hardiness, availability of highly polymorphic inbred strains, smaller genome size (800 Mb), and the adaptation of photoperiod, higher salinity, and a wide temperature range (6–40 °C) [9]. These traits contribute to the vast usage of medaka as an aquatic toxicological model [3].
Nowadays, behavioral profiles of fish have become a widely used approach in toxicology and pharmacology due to several reasons, including their high sensitivity to various chemicals during their early developmental stages and providing less cost, time, and space-consumption for in vivo drug screening compared to traditional animal models. Furthermore, juvenile and adult fish exhibit more complex behaviors, such as social interactions, learning, and memory that may directly influence the survival of individuals and future population structures, which can provide more options for the study of aquatic neurotoxicity [10,11]. It has recently become evident that medaka, especially O. latipes, possesses complex social and visually evoked behaviors, such as aggressive behavior, predation, social learning, startle response, mating preference, and shoaling, which represents the complex interaction of fishes moving together. These features emphasize medaka as an emerging model for neurobehavioral research [3]. Moreover, due to their high genetic resemblance to humans, their genetic examination progresses rapidly, opening up new approaches for studying the genetic control of behavior and establishing their role as a highly valuable human disorders study model [12,13,14,15,16,17]. Additionally, medaka also affords many beneficial traits for behavioral studies, including its central nervous system (CNS), which has a basic structure that resembles those of amniotes [18]. Finally, yet importantly, medaka also can provide invaluable data for comparative research on zebrafish [19]. However, the current behavioral profiles based on medaka mainly focus on their early stages, or only involve simple tests in aquatic toxicology research or ecotoxicological risk assessment [10]. Therefore, systemically comparing their behavioral features from simple to complex tasks is required to obtain full toxicity information of several toxicants.
Many animals exhibit inter- and intra-specific behavioral diversities, generated not only by genetic factors, but also by learning and development [20,21]. These behavioral diversities in some animal groups are thought to be a significant factor contributing to the emergence of social organization, influencing fitness, and are suggested to be under natural selection [22]. Thus, in evolutionary ecology, this type of biological variation in behavioral characteristics has attracted the interest of many researchers [18,23,24]. These variations are widely recognized in both clinical populations and animal models of human disorders. It is prevalent for the CNS phenotypes, such as neuropharmacological, behavioral, and toxic responses [25,26]. However, the reasons causing the individual differences that are maintained in natural populations, and the connection between genetic polymorphisms and behavioral diversities, remain unsolved. While interspecific behavioral diversities have been relatively well studied in rodents and primates, aquatic models are far less characterized in this case [27]. Initially, fish behavior was considered ‘simple’ and instinctive; however, it is recognized as homologous to mammals, complex, adaptive, context-specific, and highly variable nowadays [28].
Unfortunately, even though medaka is a small freshwater fish that is commonly used as an animal model for aquatic toxicology research, it has received relatively little attention in behavioral studies, specifically in interspecific behavioral diversities, compared to zebrafish [10,29]. Thus far, few reports describe the differences in the behavioral traits among the medaka species and inbred strains, including a prior study by Hyodo-Taguchi that found some inbred strains of medaka tend to be attached to humans. In contrast, other inbred strains tend to avoid interactions with humans, and these responses may be mediated by a visual stimulus [30]. However, this difference has not been quantitatively analyzed based on paradigmatic experiments [18]. In addition, a prior study showed the great phenotypic diversity of Oryzias woworae, and studied the causative gene underlying divergence in sexually selected traits [31]. Later, based on the red pectoral fins that are unique to the male of this fish, csf1 was revealed as a causative gene for red pectoral fins that can contribute to male reproductive success by the integration of genomic analysis and genome-editing technology. This result demonstrates that integrating genomic and phenomic approaches enables the identification of the causative genes underlying selected traits.
In the present study, we examine the interspecific variation in behavioral responses in adult medaka fishes with a series of behavioral assays to determine whether these fishes can be used in behavioral studies in a way similar to how zebrafish are used, or even in a more in-depth behavioral study considering the toughness of medaka over zebrafish in several extreme environmental conditions. Furthermore, biochemical studies were also carried out to help in elucidating the behavior results. We hypothesized that each species has its species-specific features, and the results of this study may provide information for other researchers to decide which medaka fish is more suitable for the study of behavioral responses, with a particular focus on neuroscience, pharmacology, and toxicology. Additionally, the present result can contribute to the growing database of phenotypical differences between several medaka fish species. The overview of the experimental design in the present study can be found in Figure A1.

2. Results

2.1. Novel Tank Assay Performance Comparison of Four Species of Medaka Fish

The novel tank assay is a test to inspect fish locomotor activity and its exploration ability to respond to the new environment [32]. Typically, fish spend most of the time at the bottom of the tank when they are introduced into a new environment, and they expand their swimming area to higher portions of the test tank after acclimating over time [33,34]. Four behavioral endpoints were used to observe their locomotor activity, which were average speed, freezing movement time ratio, swimming movement time ratio, and rapid movement time ratio. Meanwhile, regarding the exploratory behavior, the time in top duration, number of entries to the top, latency to enter the top, total distance traveled in the top, and average distance to the center of the tank (thigmotaxis) were quantified. During the test, each medaka group displayed a significantly different locomotor activity level to every other group. The highest level of locomotor activity was exhibited by O. woworae among the medaka groups. This phenomenon was shown by a high average speed, swimming, rapid movement time ratios, and a low level of freezing movement time ratio (Figure 1A–D). However, zebrafish still displayed higher locomotor activity than this medaka fish. Following O. woworae, a relatively high level of locomotor activity was also observed in Oryzias dancena. Even though not as high as O. woworae, O. dancena also displayed a higher average speed and rapid movement ratio than the other two medaka fishes (Figure 1A,D). Furthermore, a similar level of locomotor activity was detected in O. latipes and O. sinensis. While their average speed and rapid movement ratio were not statistically different from each other, different movement types were observed in these fishes, supported by the low level of freezing time movement ratio and high level of swimming time movement ratio of O. latipes over Oryzias sinensis (Figure 1B,C). Next, regarding the exploratory behavior, a different response to the novel environment from the zebrafish was displayed by all of the medaka groups. These differences were supported by the significant difference in all of the exploratory behavior-related endpoints between the zebrafish and all of the medaka fishes observed during the test (Figure 1E–I). Interestingly, O. sinensis showed a longer time in top duration than the other medaka fishes and even higher than the zebrafish (Figure 1E). However, a relatively short distance was traveled in the top indicating that this medaka swam slower than the zebrafish in the top portion of the test tank, which was also supported by a low level of average speed displayed by this fish (Figure 1A,I). Taken together, O. woworae possessed the highest locomotor activity among all of the medaka fishes, while the zebrafish locomotor activity was still higher than that of the medaka fish. Also, each medaka fish showed a specific and unique exploratory behavior, and this behavior was significantly different from the zebrafish. The detailed statistical analysis results of this test can be found in Table A1.

2.2. Aggressiveness Comparison of Four Species of Medaka Fish

Next, to evaluate the aggressiveness level of the fishes, a mirror biting assay was carried out. This measurement was performed by counting the relative interaction time of the fish with its mirror-reflecting image. Generally, to drive away from the potential intruder, fish immediately display mirror biting behavior when introduced into a tank with a mirror [35]. Similar to the novel tank test result, each medaka fish displayed a different level of aggressiveness. The most pronounced aggressive behavior was shown by O. sinensis, followed by O. woworae. This finding was indicated by a significantly higher mirror biting time percentage and the longest duration on the mirror side percentage of these two medaka fishes than the other two medaka fishes, O. dancena and O. latipes (Figure 2A,B). Interestingly, their level of aggressiveness was also found to be significantly higher than Danio rerio. On the other hand, O. dancena and O. latipes displayed a comparable level regarding their aggressiveness, which was also similar to D. rerio (Figure 2A,B). The detailed statistical analysis results of this test can be found in Table A2.

2.3. Comparison of Predator Avoidance Test Performance for Four Species of Medaka Fish

Afterward, we evaluated the fish’s fear level when facing their predators in the predator avoidance test [36]. The convict cichlid (Amatitlania nigrofasciata) was used in this study as a stimulus fish to promote the fear level of the tested fish based on the previous protocol [37]. The least predator avoidance behavior was displayed by O. dancena, followed by O. latipes. A high level of approaching predator time percentage and a low level of the average distance to the predator’s separator exhibited by these medaka fishes during the test indicated that they were not as fearful as the other fishes, including the zebrafish (Figure 2C,D). However, O. woworae and O. sinensis still displayed quite clear predator avoidance behavior on a similar level with the zebrafish. The detailed statistical analysis results of this test can be found in Table A2.

2.4. Conspecific Social Interaction Comparison of Four Species of Medaka Fish

Later, a conspecific social interaction test based on a similar rodent paradigm was conducted to evaluate the fishes social behaviors. This test is conducted by observing their interactions with the conspecifics, and it is one of the useful assays to study fish social phenotypes [35]. From the results, O. latipes displayed the slightest interest to interact with their conspecific. This phenomenon was supported by a low level of conspecific interaction time percentage and the longest conspecific interaction percentage, and a high level of average distance to the conspecific separator (Figure 3A–C). Following O. latipes, O. dancena and O. sinensis showed a slightly more pronounced conspecific social interaction during the test. The most profound conspecific interaction among the tested medaka fishes was observed in the O. woworae group. In addition, this medaka fish also exhibited a comparable level of this social behavior with the D. rerio. The detailed statistical analysis results of this test can be found in Table A2.

2.5. Comparison of Shoaling Behavior for Four Species of Medaka Fish

Shoaling, an innate behavior for several fish to swim together, was observed in each medaka fish. Generally, this behavior was intended to reduce anxiety and the risk of being captured by the predators [38,39]. O. latipes and O. dancena formed a quite tight shoal during the test, shown by the low levels of all of the behavioral endpoints, which are the average inter-fish distance, average shoal area, average nearest neighbor distance, and average farthest neighbor distance (Figure 4A–D). Meanwhile, O. woworae and O. sinensis displayed a slightly looser shoal, which is interesting since a similar shoal size was also observed in D. rerio. The detailed statistical analysis results of this test can be found in Table A2.

2.6. Comparison of Circadian Locomotor Activity Rhythms for Four Medaka Species

Next, we assessed their circadian locomotor activity rhythm since this daily rhythm of gross locomotor activity is frequently used as an assay of the circadian rhythmicity of animals, especially mammals [40]. From the results, all the medaka fishes displayed a significantly different circadian locomotor activity pattern compared to the zebrafish (Figure 5A). The high levels of locomotor activity showed these differences during both the day and night cycles. Meanwhile, regarding the medaka fish results, differences between each species were also observed during the experiment. In the day cycle, the highest locomotor activity was exhibited by O. sinensis, followed by O. dancena, O. latipes, and O. sinensis, which was consistent with the novel tank test results, shown by the differences in average speed, average angular velocity, freezing, swimming, and rapid movement time ratios (Figure 5B–G). Also, in terms of movement orientation, more pronounced zig-zag-like movement was shown by every medaka fish compared to the zebrafish. This phenomenon was shown by the high level of meandering measured during the day cycle, which was also observed during the night cycle (Figure 5J). Furthermore, similar results regarding locomotor activity were also found in the night cycle, where O. woworae possessed the highest locomotor activity among the rest of the medaka fishes, indicated by the differences in the average speed, average angular velocity, freezing, swimming, and rapid movement time ratios (Figure 5H,I,K–M). The detailed statistical analysis results of this test can be found in Table A2.

2.7. Biochemical Assay of Biomarker Expression in the Brain

Since behavior variations between each medaka fish were observed in previous behavioral tests, it was intriguing to investigate the contents of several neurotransmitters, the antioxidant activity, and oxidative stress that might play roles in these differences. Interestingly, based on the neurotransmitter levels, we found that O. woworae and O. dancena showed significantly higher levels of all of the tested neurotransmitters (serotonin (5-HT), dopamine (DA), and norepinephrine (NE)), except for acetylcholinesterase (AChE), than O. latipes and O. sinensis. Furthermore, a similar result was also observed in the stress hormone cortisol (Table 1). On the other hand, while O. woworae still showed the highest level of AChE, O. dancena possessed a comparable level of this neurotransmitter to O. sinensis, followed by O. latipes that also displayed a similar level of AChE to O. sinensis. In addition, this pattern was also found in the oxidative stress-related markers, which are catalase (CAT) and reactive oxygen species (ROS) (Table 1). The detailed statistical analysis results of this test can be found in Table A3.

2.8. PCA Analysis and Hierarchical Clustering Analysis of Several Medaka Fish Behavioral Endpoints

To explore the behavioral phenomics between several different medaka fish species, principal component analysis (PCA), hierarchical clustering, and heatmap comparison were performed after all of the behavioral tests. This process is also vital to reduce the data dimension and complexity. As the outgroup, the AB strain zebrafish (D. rerio) behavioral data were included to conduct a more profound study about the behavioral pattern differences between the tested medaka fishes. Two major clusters were generated from the hierarchical clustering result, separating O. latipes, O. dancena, and O. sinensis in one cluster, with O. woworae and D. rerio in another cluster (Figure 6A,B). The definition of all of the behavioral endpoints was described in Appendix A Table A4.

2.9. Phylogenetic Relationships of Four Medaka Species

To verify the relationship between each medaka fish species tested, we constructed a phylogenetic tree of medaka and D. rerio as the outgroup based on the NADH dehydrogenase subunit two gene (nd2). The medaka fish were classified into three groups from the phylogenetic tree (Figure 7), including celebensis, javanicus, and latipes species groups, as described by Kinoshita et al. [9].

3. Discussion

This is the first study to demonstrate the comparison of the innate behaviors between several medaka fish species to the best of our knowledge. As we predicted, each species has its species-specific features in each of the behaviors tested. O. woworae displayed the highest locomotor activity, while each medaka exhibited a unique exploratory behavior to each other. Moreover, this phenomenon was also obviously shown in the circadian locomotor activity rhythm test. Furthermore, this fish, together with O. sinensis, exhibited a significantly high aggressiveness level than D. rerio. Meanwhile, regarding predator avoidance behavior, these fishes showed similar fear-like behavior to D. rerio, while this behavior response to the predator was not clearly observed in O. dancena and O. latipes. O. woworae also showed a more distinguished social behavior in the same manner with D. rerio than other medaka fishes. Lastly, while all of the medaka fishes showed interest in forming a shoal with their conspecifics, a quite tight shoal was displayed by O. dancena and O. latipes.
The novel tank test has been widely used to study animal habituation. Generally, habituation is defined as a change in locomotor activity and exploratory behavior over time to sustain an animal’s survivability [41]. Based on the prior study, this behavior test on the zebrafish provided results comparable to other similar studies in other animal species, including medaka [42]. The typical behavioral pattern changes of all of the medaka, in a stereotypical manner, were demonstrated by the current results. This phenomenon indicated that the animals became familiar with the novel location as the exposure time increased. Consistent with this result, a previous study demonstrated a similar pattern change during the habituation of O. latipes in an open-field test [19]. In terms of their differences from the zebrafish, all of the medaka fishes exhibited a relatively lower locomotor activity than the zebrafish most of the time. This result is somewhat similar to a prior study in another species of medaka (Oryzias javanicus), which found that the medaka fish were passive compared to the D. rerio who were more active and aggressive [43]. Moreover, a high thigmotaxis level, a preference of animals towards the periphery of a novel arena and avoiding the center area, during the test was displayed by all of the medaka fish [44]. This result is plausible since zebrafish tend to spend time in the center area of a tank after it is habituated, while this behavior is not shown in medaka [45]. All of the medaka fishes displayed similar phases of habituation in the novel environment. This finding confirmed the usefulness of medaka and novel tank tests to investigate the habituation phenotype even though each species possessed different degrees of locomotor activity and exploratory behavior.
Similar to zebrafish, mirror approaching behavior in medaka is a well-established fish paradigm that reflects their aggressive behavior [46]. When a single medaka is placed in a tank, it swims freely in all directions; however, the same individual will tend to swim close to the mirror when a mirror is placed on one side of the tank [47]. Unfortunately, while this paradigm is well studied in medaka, most of the studies were done in O. latipes [3,47]. Thus, this is the first study that revealed the mirror approaching behavior in other medaka species. Surprisingly, in this study, while O. latipes displayed a similar mirror approaching behavior to D. rerio, other medaka fishes exhibited more robust mirror approaching behaviors. Since mirror approaching behavior in medaka is considered a simple and robust model of socially induced anxiety, the high level of mirror biting time might indicate a higher stress level in these fishes than O. latipes and D. rerio. Interestingly, cortisol, one of the primary fish hormones, was found to be lower in O. latipes compared to other medaka fishes, especially O. dancena and O. latipes. These observations indicate that cortisol might affect the differences in their aggressive behavior since it is involved in the medaka fish stress response [48]. Also, this speculation was based on the prior study in diazepam and fluoxetine, which are anxiolytic agents that decreased the mirror biting time and the stress response measured by the cortisol level [3,47,49]. Furthermore, based on the ELISA results on ROS, the molecules responsible for the signaling stress response, we hypothesize that different basal levels of ROS might also influence the differences in their behavior. This possibility was based on the previous research in adult zebrafish that showed anxiety-like or stress-like behavior, which was likely contributed to by higher ROS levels [50,51]. Besides ROS, this speculation is also supported by the variation in the catalase basal levels of each medaka fish. Catalase is primarily a peroxisomal enzyme that catalyzes the enzymatic decomposition of H2O2, and in rats it is associated with depression-like behavior associated with Alzheimer’s disease improvement [52,53]. Lastly, aggression also has been linked to serotonergic function in a variety of invertebrate and vertebrate species. Generally, a high level of serotonin activity clearly shown in the results is associated with low levels of aggressive behavior [54].
Avoiding predation is an essential behavioral reaction critical for survival, and may have a significant fitness component. Many animals, including fish, sense a danger of predation through a multitude of methods, such as alarm signals, calls, and chemical cues. To avoid predation, the animals have to respond to stimuli that represent the presence of danger properly [36,55]. It not only helps the researchers to elucidate the evolution and ecology of the studied species, but analysis of such fear responses may have clinical relevance [56]. Unfortunately, regardless of the promising advantages, only a few studies of this predation response were done in medaka. Moreover, most of these studies only used O. latipes as an animal model [55,57]. However, in the current study, we found that each medaka species responded to the presence of A. nigrofasciata as the fear stimulus in different magnitudes. Interestingly, while the O. latipes fear responses were clearly demonstrated in a prior study, this behavior was not shown in this experiment. This difference might be due to the different stimuli used to elicit the alarm response. While we stimulated the fear response by exploiting the visual cue from another fish, the prior study used a conspecific skin extract as the chemical cue of alarm substance [55]. Thus, each medaka fish might respond differently to other fish species or a different stimulus upon their first exposure to these fish, which needs to be confirmed in future studies. In addition, a more robust reaction in response to the predator stimulus was observed in O. woworae and O. sinensis. These results might indicate that these medaka fishes displayed a fear response in the presence of the stimulus since their response was similar to D. rerio, which is already well studied and proven to elicit a fear response in a similar situation [37,58,59]. One possibility that caused the differences is related to the specific basal serotonin content of each medaka fish, as shown in the current study. As mentioned above, serotonin is a good physiological indicator of various types of stress in fish. Mosienko et al. demonstrated that serotonin deficiency in certain regions of the brain, exhibited by tph2-knockout medaka, was associated with elevated stress and fear-related behaviors [60].
Generally, the development of social behavior is observed among various vertebrates, from fish to mammals, and coordinated with social factors such as the presence of conspecific, physiologic, and environmental factors [61,62]. In medaka fish, especially O. latipes, this social behavior had been well addressed previously. A prior study discovered that medaka fish robustly and reproducibly maintained proximity to a single target conspecific fish since it is attracted by conspecifics’ biological motion [63,64]. However, whether other medaka species also possess this social behavior by the presence of conspecific remains unknown. Even though O. latipes was found to maintain proximity to its conspecific, other medaka species also displayed a similar behavior, even in a higher magnitude and more robust manner. These results are intriguing since some of these medaka species, especially O. woworae, displayed social behavior as clear as D. rerio, a well-established animal model for this behavior [37,65]. Furthermore, shoaling, a simple form of affective behavior displayed in social fish species, is usually observed in small fish, including medaka (O. latipes) and zebrafish. This behavior may directly influence the survival of individuals and future population structures [11]. In the previous study, medaka was demonstrated to form shoals only with conspecifics and not with other species [66]. In this study, all medaka fishes formed a shoal that was even tighter than D. rerio for some medaka species, which were O. dancena and O. latipes. These results are plausible since medaka are known to have high visual acuity and exhibit a strong tendency to form shoals [67]. The differences in the social behaviors of each medaka fish might be related to the differences in their dopamine and norepinephrine levels, as shown in the ELISA results. These neurotransmitters are involved in the organization of the stress response in vertebrates, and together with serotonin, they are involved in the mediation of behavioral and neuroendocrine responses to social stress [68,69]. A prior study in teleost shows that the brain dopamine level increased after being stressed. Thus, the current results indicate that some medaka fishes might have different levels of stress compared to the others, which is consistent with the mirror biting test results. The stress is plausible since recent research discovered that grouping and pairing are more stressful to male medaka than isolation [48]. Regarding the norepinephrine content, the previous finding found that juvenile lake sturgeon (Acipenser fulvescens) held in isolation had a significantly longer norepinephrine response than fish held with conspecifics [70]. In addition, another study in zebrafish found that the whole-brain dopamine level was also found to be associated with the development of shoaling [71]. In a prior study, the dopamine level was found to have rapidly risen in the brain of adult zebrafish in response to social stimuli [72].
The time of oviposition in medaka fish, especially in O. latipes, was well-accounted several years ago. Interestingly, this rhythm is controlled by a circadian pacemaker entrained by light–dark (LD) cycles since it was gradually disappeared under continuous light [73]. The locomotor activity rhythm seems to be dependent on the light intensity during dark periods [74]. Here, typical diurnal circadian locomotor activity rhythms were displayed by all of the medaka fishes. Even though they were not as pronounced as in D. rerio, they are consistent with their natural behaviors [75]. This phenomenon might be explained by the circadian rhythm trait of medaka fish, especially O. latipes, which can be entrained by feeding time [76]. Thus, the feeding schedule that occurred during the day cycle in the present study might play a role in affecting this behavior since a prior study demonstrated that this entrainable rhythm is shown to be a persistent behavior. In their study, it was discovered that after a three-day fast, a meal-feeding readily entrained circadian rhythm of agonistic behavior of O. latipes remained fixed to the feeding time [77]. However, the reader should not conclude that all medaka species are capable of being entrained to feeding. Thus, future studies are needed to be conducted to verify this speculation. Furthermore, O. woworae was found to be the most active fish regarding their locomotor activity levels compared to other medaka fishes, which is consistent with the novel tank test results. This phenomenon could be related to the high level of serotonin and AChE measured in the fish. This speculation was taken since serotonin has been reported to play a positive role in regulating locomotor activity in many animals, while AChE is an enzyme responsible for the breakdown of acetylcholine, which plays an important role in modulating neuromuscular activity in neural synapses [78]. This trend is also observed in other medaka fishes.
In summary, the behavioral differences between each medaka fish might be associated with their different basal level of neurotransmitters. This phenomenon has been briefly noted in a previous study. However, the comparison only covered three kingdoms of living organisms [79]. Nevertheless, even though this difference was not fully investigated, it has been mentioned in several prior studies, including in fish. In normoxia conditions, three fish species, goldfish, tilapia, and carp, possessed different energy status levels, lactic acid accumulation, and amino acid patterns in brain tissue [80]. In addition, by using a combination of immunohistochemistry and confocal imaging, a previous experiment confirmed the differences in the serotonin and acetylcholine contents in the gill’s filament and lamellae of various fish species, including goldfish (Carasius auratus), sockeye salmon (Oncorhynchus nerka), trout, and O. latipes [81]. Furthermore, several prior findings also found some variations in the levels of adrenaline and noradrenaline stored within the chromaffin tissue of various species of fish, such as cyclostomes, dipnoans, elasmobranchs, ganoids, and teleosts [82]. Lastly, Aprison et al. also found some differences regarding glycine concentration, an amino acid that has potent inhibitory actions on the vertebrate nervous system, at various levels of the neuraxis in five different vertebrates [83]. Additionally, these phenomena might be related to the different genes underlying the divergences in selected traits between species, particularly in vertebrates, which needs further studies to confirm this hypothesis [84]. Overall, this is the first study comparing the basal levels of neurotransmitters between each species of medaka fishes to the best of our knowledge. Thus, further studies are required to deeply investigate this phenomenon.
Next, the phenomic-based hierarchical cluster generated two major clusters, which were O. latipes, O. dancena, and O. sinensis in one cluster, and O. woworae and D. rerio in another cluster (Figure 6A,B). In the first cluster, O. latipes and O. dancena were closer to each other than O. sinensis. This categorization was plausible since these medaka fishes exhibited several behaviors that were not observed in other fish, such as less pronounced predator avoidance behavior (endpoint 3-1) and tightened shoal formation (endpoint 5-1 to 5-4). Meanwhile, even though O. sinensis was found to be in this cluster, it showed a different behavioral pattern to the two medaka fishes in several vital endpoints. The most evident difference was in one of the novel tank test endpoints, the time in top duration (endpoint 1-1-5 and 1-2-5). As already mentioned in Figure 1, a significantly high level of this endpoint was displayed by this medaka fish during the whole 30 min of the novel tank test. Next, it was found that O. woworae and D. rerio belonged in the same cluster. After further investigation, we discovered that these two fishes exhibited similar behavior patterns, especially in the locomotor activity-related endpoints, such as the average speed and rapid movement ratio. A strong conspecific interaction shown by these fishes also became one of the factors that defined this grouping. Afterward, to precisely determine the phylogenetic position of each medaka fish and D. rerio, we estimated the phylogenetic relationship using published whole-genome datasets as references. The tree indicates that O. sinensis, O. latipes, and O. dancena formed a monophyletic group, which is consistent with previous trees based on a concatenated mitochondrial sequence matrix [85]. Interestingly, this genetic-based grouping follows with the salinity tolerances of these fishes, which are euryhaline for O. dancena, O. latipes, and O. sinensis [86,87,88], and stenohaline for O. woworae and D. rerio [89,90]. This result is plausible since O. sinensis is described as the subspecies O. latipes sinensis Chen. These species have approximately the same length of anal–fin rays and chromosome arms numbering 58 or more, distinguishing them from putative close relatives [88,91]. Furthermore, O. dancena and O. woworae were phylogenetically clearly separated from O. sinensis and O. latipes, which is also in agreement with other phylogenetic trees from previous studies that are based on the nuclear tyrosinase, mitochondrial 12S and 16s rRNA genes, and several characters. Based on these studies, the Oryzias species itself is ramified into three monophyletic groups. While O. dancena belongs to the javanicus clade, both O. sinensis and O. latipes belong to the latipes clade, which can elucidate the phylogenetic tree result of the present study. The javanicus clade is a monoarmed chromosome group and possesses subtelo- and acro-centric chromosomes that distinguish it from the other clades [9,85,88,92]. Next, as reported in a phylogeny based on the 3440-bp concatenated mitochondrial and nuclear sequences by Mokodongan and Yamahira, O. woworae was genetically distinct from other medaka fishes, such as O. latipes and O. javanicus. This grouping might also be caused by the different wider salinity tolerance between the medaka fishes. While medaka fishes, including O. latipes, O. sinensis, and O. dancena, have a wide salinity tolerance, this trait was not possessed by O. woworae, which might also explain the close distance between this fish and D. rerio in the phylogenetic tree. However, based on the preliminary phylogenomic analyses, the distance between zebrafish and medaka is due to the differences in the 20 nuclear protein-coding genes [93,94,95,96]. Interestingly, the phylogenetic tree and phenomic-based hierarchical cluster generated in this study displayed several resemblances, including the interspecies relationship between zebrafish and the medaka fishes. However, a slightly different result to the phylogenetic tree was observed in the medaka interspecies comparisons. The phenomic-based hierarchical clustering result showed that O. latipes and O. dancena appeared to be more closely related to each other than either of them was to O. sinensis. Thus, this is evidence that genetic-based phylogenetic analyses might have a different relationship among the species to the phenomic-based cluster analysis. However, a prior study has already addressed those differences regarding the relationship between zebrafish, medaka, pufferfish, and cichlids. In their study, the molecular data supported a close relationship between atherinomorphs (including the medaka and platy) and putatively more derived perform fish such as the cichlids. At the same time, previously the Atherinomorpha have been historically placed in an intermediate position among the other branches of the acanthomorph tree since they share several putative ‘primitive’ morphological features with more basal teleosts. Thus, their phylogenomic analysis revealed an unexpected relationship among the other three species, contrary to traditionally held systematic views based on morphology [96].
Next, to observe their reproducibility, we also calculated the coefficient of variation of each medaka fish species and D. rerio from every behavioral test. Based on the result in Table A5, O. woworae displayed the lowest average coefficient of variation between all of the medaka fishes, even though it was still higher than D. rerio. This result indicates the stronger reproducibility of O. woworae than other medaka fish in these behavior tests. A high coefficient of variation may impact the animal usage of an experiment since this condition is a problem to obtain significant results. Thus, a large sample size has become mandatory to overcome this problem. Moreover, this problem also means that the researchers might face ethical-related issues [29,97].

4. Materials and Methods

4.1. Animal Husbandry

Indian ricefish (Oryzias dancena) and Daisy’s ricefish (Oryzias woworae) were purchased from a local pet store, while Chinese ricefish (Oryzias sinensis) and Japanese ricefish (Oryzias latipes) were obtained from the Freshwater Bioresource Center at National Chiayi University. The AB strain zebrafish were obtained from the Taiwan Zebrafish Core Facility at Academia Sinica (http://icob.sinica.edu.tw/tzcas/, accessed on 26 November 2020). All fish in each experiment were mixed gender of 6–10 months old adults in healthy condition. One month prior to the experiment, all tested fish species were reared in the centralized fish facility with a recirculating aquatic system at 28 ± 1 °C and a 10/14-h dark/light cycle. This process was important to eliminate some external factors that might affect their behavior performance, such as the stressful condition during the transfer [98]. The conductivity of the circulating system’s water was kept between 300–1500 µS with pH 7.0–7.5. Ultraviolet (UV) light was utilized to filter the water constantly. All fish were held and raised in a trapezoid plastic tank with 34 cm at the top, 23 cm along the bottom, 19 cm along the diagonal side, 18 cm high, and 27 cm wide filled with 8 L of filtered water. Feed was given twice a day (09:00 and 17:00) with either lab-grown brine shrimp or commercial dry food. The general maintenance procedures and housing conditions were as previously described by Avdesh et al. [99].

4.2. Animal Ethics and Behavioral Tests

All fish experiments were performed following the guidelines issued by the Institutional Animal Care and Use Committees (IACUCs) of Chung Yuan Christian University (application number: CYCU106025, issue date 6 May 2018). All behavioral tests were conducted within the morning until afternoon (10:00 to 16:00), except for circadian locomotor activity rhythm, in a temperature-controlled room (26 ± 1 °C). After the acclimation process, a battery of behavioral tests, which were the novel tank, mirror biting, predator avoidance, conspecific social interaction, and shoaling tests based on the previous method, was conducted in all the fish groups [37]. A total number of 142 fish were used in behavior tests (n = 30 for D. rerio, O. dancena, O. latipes, and O. woworae; n = 22 for O. sinensis). The slightly lower number of O. sinensis used compared to other fishes was caused by the low availability of this fish in our region during the time period. In addition, the differences in sample size number of fish groups in several behavior tests were due to the unknown death of the fish during the rest period between the behavior tests. The entire tests were done using a trapezoid tank with 22 cm along the bottom, 28 cm at the top, 15.2 cm high, and 15.9 cm along the diagonal side, and later, this tank was filled with ~1.25 L of circulating filtered water. In avoiding stress conditions in the tested fish, all the behavior tests, excluding the circadian rhythm locomotor activity test, were divided into three sessions and carried out on a single week. The first session consisted of the novel tank and shoaling tests. A novel tank test was conducted as the first behavior test since the tested fish had not been introduced to the test tank previously, thus, their behavioral responses toward the novel environment would be more genuine and reliable. The mirror biting and social interaction tests belonged to the second session, followed by the predator avoidance test in the last session since in this test, the convict cichlid (Amatitlania nigrofasciata) was supposed to induce the fear response on the tested fish and since fear can be possessed by the fish for some period, it could affect their behavior in the following behavior tests. Before all of the tests, except the novel tank test, 3–5 min of acclimation in the test tank was applied to all tested fishes. Afterward, a 5-min video recording was conducted. Meanwhile, in the novel tank test, fish behavior was recorded for 1 min as soon as the fish were exposed to the test tank. The behaviors were recorded at 5 min intervals for 31 min. Next, the behavior tests were continued with the circadian locomotor activity rhythm test on the following week, which monitored its locomotor activity for 24 h [100]. This test was conducted in a 30 × 30 × 7.5 cm acrylic tank filled with ~3 L of filtered water. A lightbox, which consisted of two types of a light source (chip-on-board (COB) light-emitting diode (LED) and 940 nm infrared LED), was placed below the tank. Later, the fish locomotion during 24-h test was recorded for 1 min every hour. Next, idTracker (http://www.idtracker.es/, accessed on 25 May 2020), open-source software, was used to collect and convert the fish movement data to trajectories from all of the recorded videos [101]. All fish behavior tests were done in triplicate.

4.3. Brain Tissue Preparation, Total Protein Determination, and Quantification of Neurotransmitters, Stress Hormones, and Oxidative Stress Markers

In obtaining the brain tissues, immediate anesthesia and euthanasia were performed in medaka fishes by immersing them in tricaine solution (A5040, Sigma, St. Louis, MO, USA). Later, biochemical analyses were conducted on their whole-brain tissue extract. Ice-cold phosphate-buffered saline (PBS) in volumes of 10 (v/w) was used to standardize a single homogenate of two to three whole medaka fish brains at pH 7.2. Afterward, to homogenize the tissue a bullet blender (Next Advance, Inc., Troy, NY, USA) was utilized. After 15 min of 13,000 rpm centrifugation, the supernatant was transferred to a sterilized microtube and stored at −20 °C. Subsequently, Pierce BCA protein assay kit (23225, Thermo Fisher Scientific, Massachusetts, MA, USA) was applied to measure the total brain tissue’s protein level. After the color was formed, it was analyzed by using a microplate reader (Multiskan GO, Thermo Fisher Scientific, Waltham, MA, USA) at 562 nm. Subsequently, all medaka fish brain tissues were analyzed to compare the differences in neurotransmitters, stress hormone, and oxidative stress marker levels between each species of the medaka fishes. Several neurotransmitters, including serotonin (5-HT, ZGB-E1572), acetylcholine esterase (AChE, ZGB-E1637), dopamine (DA, ZGB-E1573), and norepinephrine (ZGB-E1571), were measured by target-specific ELISA kits. Meanwhile, tissue oxidative and anti-oxidative stress markers, reactive oxygen species (ROS) and catalase (CAT), respectively, and cortisol, one of the stress hormones, were also quantified by ELISA kits (ZGB-E1561, ZGB-E1598, and ZGB-E1575 Zgenebio Inc., Taipei, Taiwan). Later, a microplate reader (Multiskan GO, Thermo Fisher Scientific, Waltham, MA, USA) was utilized to measure the absorbance at 450 nm and compared it to the standard curve to quantify the relative concentration of the target protein. Ten biological and three technical replicates were applied in the analysis (n = 30, except O. latipes (n = 29)).

4.4. Statistical Analyses

In determining the statistical difference between each group, statistical analyses were carried out using GraphPad Prism (GraphPad Software version 8 Inc., La Jolla, CA, USA). For all behavioral and biochemical data analyses except the novel tank test, Kruskal–Wallis with uncorrected Dunn’s test was conducted to find the statistical differences in each group with every other group since the data are not normally distributed [102]. Meanwhile, two-way ANOVA with Geisser–Greenhouse correction continued, and uncorrected Fisher’s LSD test was used for the novel tank test. Data for each fish group are expressed as either median with interquartile range, the mean with a 95% confidence interval (CI), or standard deviation (SD). The statistic details for each behavioral and biochemical test are summarized in Table A2 and Table A3, respectively.

4.5. PCA, Heatmap, and Clustering Analysis

All the behavioral endpoint values from all the tested fish in every behavior test were input into a comma-separated values type file (.csv) using Microsoft Excel. All of the essential behavioral endpoints based on the prior study are listed and explained in Table A4 [37]. Next, the .csv file was uploaded to ClustVis (https://biit.cs.ut.ee/clustvis, accessed on 19 November 2020), a web tool designed to visualize and cluster multivariate data. Later, unit variance scaling for each row was carried out in order to treat each variable equally. Furthermore, singular value decomposition (SVD) with the imputation method was used to calculate principal components since there were no missing values in the dataset [103].

4.6. Phylogenetic Tree Construction

The complete mitochondrial NADH dehydrogenase subunit 2 gene (nd2) sequences for each species of medaka fish were downloaded from the NCBI database with the accession numbers NC_012976 (O. dancena), NC_004387 (O. latipes), LC051726 (O. woworae), and NC_013434 (O. sinensis). The gene sequences of nd2 of D. rerio (NC_002333), representing outgroup species, were also downloaded from the NCBI database. All nucleotide sequences nd2 of all species described above were then aligned using Geneious software (Biomatters, Auckland, New Zealand). Subsequently, phylogenetic analysis was also conducted with Geneious software with default parameter settings.

5. Conclusions

To sum up, the current study demonstrated that the phenotypes of four medaka fish, which were O. latipes, O. dancena, O. woworae, and O. sinensis, were different in their behavior in various behavioral tests. These differences might be related to the different basal levels of several neurotransmitters, the stress hormone, and oxidative stress markers measured. Furthermore, while our phenomic-based hierarchical cluster results showed a similar interspecies relationship between zebrafish and medaka fishes to the phylogenetic tree, the interspecies comparison between each medaka fish of these two approaches demonstrated a slight difference in the relationship result between three medaka fishes. Our findings also suggest that many challenges remain to be addressed before robustly using the behavioral endpoints for ecological hazard evaluation, since the instability of some behavior performances can be observed in the medaka fishes. We believe that our results contribute to the growing database of phenotypical differences between several medaka fish species and provide one of the foundations for future phenomics studies of medaka fish. Moreover, comparative studies involving zebrafish and medaka are remarkably informative for identifying the highly conserved genetic control mechanisms [104,105]. Furthermore, considering the abundance of medaka fish strains listed in the National BioResource Project (NBRP) Medaka, it is intriguing to study the differences between each medaka strain in every aspect, including the relationship between their genetics, behavior, and biochemistry [106].

Author Contributions

Conceptualization, G.A. and C.-D.H.; methodology, G.A.; software, P.S.; validation, G.A., K.H.-C.C., and J.-C.H.; formal analysis, P.S.; investigation, G.A.; resources, H.-T.L.; data curation, G.A.; writing—original draft preparation, G.A. and C.-D.H.; writing—review and editing, M.J.M.R., K.H.-C.C., and J.-C.H.; visualization, G.A.; supervision, C.-D.H., and H.-T.L.; project administration, C.-D.H.; funding acquisition, C.-D.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology, grant number MOST 108-2313-B-033-001-MY3 and MOST 107-2622-B-033-001-CC2 to C.-D.H and grant number MOST 107-2633-B-153-001- to J.-C.H.

Institutional Review Board Statement

The study was conducted according to the guidelines approved by the Institutional Animal Care and Use Committees (IACUCs) of Chung Yuan Christian University (application number: CYCU106025, issue date 6 May 2018).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Appendix A and further inquiries can be directed to the corresponding authors.

Acknowledgments

The authors appreciate Freshwater Bioresource Centre at National Chiayi University for providing the medaka fishes and Taiwan Zebrafish Core Facility at Academia Sinica (TZCAS) for providing the zebrafish AB strain.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Overview of the experimental design in the present study.
Figure A1. Overview of the experimental design in the present study.
Ijms 22 05686 g0a1
Table A1. The individual P-value of each fish behavioral endpoint comparison in the novel tank test. The P-values were obtained from a two-way ANOVA test with Geisser–Greenhouse correction continued with uncorrected Fisher’s LSD test. p-values less than 0.05 are written in bold.
Table A1. The individual P-value of each fish behavioral endpoint comparison in the novel tank test. The P-values were obtained from a two-way ANOVA test with Geisser–Greenhouse correction continued with uncorrected Fisher’s LSD test. p-values less than 0.05 are written in bold.
Behavioral Endpoints O. dancenaO. latipesO. woworaeO. sinensisD. rerio
Average SpeedO. dancena 0.00970.0005<0.0001<0.0001
O. latipes0.0097 <0.00010.1134<0.0001
O. woworae0.0005<0.0001 <0.0001<0.0001
O. sinensis<0.00010.1134<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Freezing Time Movement RatioO. dancena 0.2654<0.00010.0155<0.0001
O. latipes0.2654 0.00020.0002<0.0001
O. woworae<0.00010.0002 <0.0001<0.0001
O. sinensis0.01550.0002<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Swimming Time Movement RatioO. dancena 0.10900.00020.0286<0.0001
O. latipes0.1090 0.0171<0.0001<0.0001
O. woworae0.00020.0171 <0.0001<0.0001
O. sinensis0.0286<0.0001<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Rapid Time Movement RatioO. dancena <0.0001<0.00010.0007<0.0001
O. latipes<0.0001 <0.00010.3643<0.0001
O. woworae<0.0001<0.0001 <0.00010.4803
O. sinensis0.00070.3643<0.0001 <0.0001
D. rerio<0.0001<0.00010.4803<0.0001
Time in TopO. dancena 0.06360.0595<0.00010.0164
O. latipes0.0636 0.9472<0.0001<0.0001
O. woworae0.05950.9472 <0.0001<0.0001
O. sinensis<0.0001<0.0001<0.0001 <0.0001
D. rerio0.0164<0.0001<0.0001<0.0001
Number of Entries to The TopO. dancena 0.0050<0.00010.41430.0619
O. latipes0.0050 <0.00010.0461<0.0001
O. woworae<0.0001<0.0001 <0.0001<0.0001
O. sinensis0.41430.0461<0.0001 0.0057
D. rerio0.0619<0.0001<0.00010.0057
Latency to Enter the TopO. dancena 0.26740.4705<0.0001<0.0001
O. latipes0.2674 0.0648<0.0001<0.0001
O. woworae0.47050.0648 <0.00010.0002
O. sinensis<0.0001<0.0001<0.0001 0.5059
D. rerio<0.0001<0.00010.00020.5059
Total Distance Traveled in The TopO. dancena 0.12330.22730.0753<0.0001
O. latipes0.1233 0.00610.0009<0.0001
O. woworae0.22730.0061 0.5744<0.0001
O. sinensis0.07530.00090.5744 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
ThigmotaxisO. dancena 0.0162<0.0001<0.0001<0.0001
O. latipes0.0162 <0.0001<0.00010.0005
O. woworae<0.0001<0.0001 0.2250<0.0001
O. sinensis<0.0001<0.00010.2250 <0.0001
D. rerio<0.00010.0005<0.0001<0.0001
Table A2. The individual P-value of each fish behavioral endpoint comparison in behavior tests. The p-values were obtained from the Kruskal–Wallis test continued with uncorrected Dunn’s LSD test. p-values less than 0.05 are written in bold.
Table A2. The individual P-value of each fish behavioral endpoint comparison in behavior tests. The p-values were obtained from the Kruskal–Wallis test continued with uncorrected Dunn’s LSD test. p-values less than 0.05 are written in bold.
Behavior TestsBehavioral Endpoints O. dancenaO. latipesO. woworaeO. sinensisD. rerio
Mirror Biting TestMirror Biting Time PercentageO. dancena 0.09880.30110.00400.1159
O. latipes0.0988 0.0073<0.00010.9375
O. woworae0.30110.0073 0.05420.0091
O. sinensis0.0040<0.00010.0542 <0.0001
D. rerio0.11590.93750.0091<0.0001
Longest Duration in the Mirror Side PercentageO. dancena 0.76560.04660.00230.4318
O. latipes0.7656 0.09070.00540.2782
O. woworae0.04660.0907 0.22100.0055
O. sinensis0.00230.00540.2210 0.0002
D. rerio0.43180.27820.00550.0002
Predator Avoidance TestApproaching Predator Time PercentageO. dancena 0.71110.16690.00160.0104
O. latipes0.7111 0.31150.00440.0285
O. woworae0.16690.3115 0.04480.2384
O. sinensis0.00160.00440.0448 0.3064
D. rerio0.01040.02850.23840.3064
Average Distance to the Predator’s SeparatorO. dancena 0.13600.01530.00170.1091
O. latipes0.1360 <0.0001<0.00010.0020
O. woworae0.0153<0.0001 0.26280.4109
O. sinensis0.0017<0.00010.2628 0.0710
D. rerio0.10910.00200.41090.0710
Social Interaction TestConspecific Interaction Time PercentageO. dancena 0.02040.06750.50320.0155
O. latipes0.0204 <0.00010.1432<0.0001
O. woworae0.0675<0.0001 0.01870.5531
O. sinensis0.50320.14320.0187 0.0038
D. rerio0.0155<0.00010.55310.0038
Average Distance to the Conspecific’s SeparatorO. dancena 0.20370.03490.49020.0294
O. latipes0.2037 0.00070.63180.0006
O. woworae0.03490.0007 0.00850.9450
O. sinensis0.49020.63180.0085 0.0071
D. rerio0.02940.00060.94500.0071
Longest Conspecific Interaction PercentageO. dancena 0.09690.03600.70530.2989
O. latipes0.0969 0.00020.25060.0070
O. woworae0.03600.0002 0.02110.2902
O. sinensis0.70530.25060.0211 0.1823
D. rerio0.29890.00700.29020.1823
Shoaling TestAverage Inter-fish DistanceO. dancena 0.3414<0.00010.0014<0.0001
O. latipes0.3414 0.00210.01930.0023
O. woworae<0.00010.0021 0.64980.9773
O. sinensis0.00140.01930.6498 0.6685
D. rerio<0.00010.00230.97730.6685
Average Shoal AreaO. dancena 0.18110.00430.00360.0010
O. latipes0.1811 0.12940.08880.0503
O. woworae0.00430.1294 0.74470.6597
O. sinensis0.00360.08880.7447 0.9410
D. rerio0.00100.05030.65970.9410
Average Nearest Neighbor DistanceO. dancena 0.0654<0.00010.0001<0.0001
O. latipes0.0654 0.00560.02970.0066
O. woworae<0.00010.0056 0.73500.9572
O. sinensis0.00010.02970.7350 0.7721
D. rerio<0.00010.00660.95720.7721
Average Farthest Neighbor DistanceO. dancena 0.5503<0.00010.0001<0.0001
O. latipes0.5503 0.00060.00090.0008
O. woworae<0.00010.0006 0.84180.9270
O. sinensis0.00010.00090.8418 0.7774
D. rerio<0.00010.00080.92700.7774
Circadian Locomotor Activity Rhythm—Day CycleAverage SpeedO. dancena 0.04320.0403<0.0001<0.0001
O. latipes0.0432 <0.0001<0.0001<0.0001
O. woworae0.0403<0.0001 <0.0001<0.0001
O. sinensis<0.0001<0.0001<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Average Angular VelocityO. dancena 0.40150.9343<0.0001<0.0001
O. latipes0.4015 0.3569<0.0001<0.0001
O. woworae0.93430.3569 <0.0001<0.0001
O. sinensis<0.0001<0.0001<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
MeanderingO. dancena 0.06030.00970.5545<0.0001
O. latipes0.0603 <0.00010.2761<0.0001
O. woworae0.0097<0.0001 0.0037<0.0001
O. sinensis0.55450.27610.0037 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Freezing Movement Time RatioO. dancena 0.33570.3036<0.0001<0.0001
O. latipes0.3357 0.0464<0.0001<0.0001
O. woworae0.30360.0464 <0.0001<0.0001
O. sinensis<0.0001<0.0001<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Swimming Movement Time RatioO. dancena 0.13370.2348<0.0001<0.0001
O. latipes0.1337 0.0072<0.00010.0086
O. woworae0.23480.0072 <0.0001<0.0001
O. sinensis<0.0001<0.0001<0.0001 <0.0001
D. rerio<0.00010.0086<0.0001<0.0001
Rapid Movement Time RatioO. dancena 0.6290<0.00010.0087<0.0001
O. latipes0.6290 <0.00010.0022<0.0001
O. woworae<0.0001<0.0001 0.0066<0.0001
O. sinensis0.00870.00220.0066 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Circadian Locomotor Activity Rhythm—Night CycleAverage SpeedO. dancena 0.06020.0013<0.0001<0.0001
O. latipes0.0602 <0.00010.0043<0.0001
O. woworae0.0013<0.0001 <0.0001<0.0001
O. sinensis<0.00010.0043<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Average Angular VelocityO. dancena 0.17340.03310.2380<0.0001
O. latipes0.1734 0.00050.9701<0.0001
O. woworae0.03310.0005 0.0020<0.0001
O. sinensis0.23800.97010.0020 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
MeanderingO. dancena 0.49480.0310<0.0001<0.0001
O. latipes0.4948 0.00450.0002<0.0001
O. woworae0.03100.0045 <0.0001<0.0001
O. sinensis<0.00010.0002<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Freezing Movement Time RatioO. dancena 0.01190.0208<0.0001<0.0001
O. latipes0.0119 <0.00010.0012<0.0001
O. woworae0.0208<0.0001 <0.0001<0.0001
O. sinensis<0.00010.0012<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Swimming Movement Time RatioO. dancena 0.01060.0226<0.0001<0.0001
O. latipes0.0106 <0.00010.0009<0.0001
O. woworae0.0226<0.0001 <0.0001<0.0001
O. sinensis<0.00010.0009<0.0001 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Rapid Movement Time RatioO. dancena 0.0834<0.00010.0279<0.0001
O. latipes0.0834 <0.00010.5159<0.0001
O. woworae<0.0001<0.0001 0.0005<0.0001
O. sinensis0.02790.51590.0005 <0.0001
D. rerio<0.0001<0.0001<0.0001<0.0001
Table A3. The individual P-value of each medaka fish’s neurotransmitters, antioxidant activity, and oxidative stress content comparisons in the biochemical tests. The p-values were obtained from a Kruskal–Wallis test continued with uncorrected Dunn’s test. p-values less than 0.05 are written in bold.
Table A3. The individual P-value of each medaka fish’s neurotransmitters, antioxidant activity, and oxidative stress content comparisons in the biochemical tests. The p-values were obtained from a Kruskal–Wallis test continued with uncorrected Dunn’s test. p-values less than 0.05 are written in bold.
BiomarkersMedaka FishesO. dancenaO. latipesO. woworaeO. sinensis
5-HTO. dancena 0.04070.15690.0488
O. latipes0.0407 0.00050.939
O. woworae0.15690.0005 0.0007
O. sinensis0.04880.9390.0007
CATO. dancena 0.0030.86330.0093
O. latipes0.0030 0.00520.7163
O. woworae0.86330.0052 0.0151
O. sinensis0.00930.71630.0151
NEO. dancena 0.00360.50320.0583
O. latipes0.0036 0.00030.3107
O. woworae0.50320.0003 0.0104
O. sinensis0.05830.31070.0104
DAO. dancena 0.00590.31990.0322
O. latipes0.0059 0.00020.5405
O. woworae0.31990.0002 0.0017
O. sinensis0.03220.54050.0017
CortisolO. dancena 0.01160.21380.0168
O. latipes0.0116 0.00020.8935
O. woworae0.21380.0002 0.0003
O. sinensis0.01680.89350.0003
AChEO. dancena 0.03890.13080.3293
O. latipes0.0389 0.00030.2756
O. woworae0.13080.0003 0.0129
O. sinensis0.32930.27560.0129
ROSO. dancena 0.02170.06920.0692
O. latipes0.0217 <0.00010.6325
O. woworae0.0692<0.0001 0.0003
O. sinensis0.06920.63250.0003
Table A4. Summary of fish behavioral endpoints measured in each behavioral test during this experiment.
Table A4. Summary of fish behavioral endpoints measured in each behavioral test during this experiment.
IndexBehavior Endpoints (Units)DefinitionApplied to
1-1-1 & 1-2-1Average speed (cm s−1)Total distance traveled by fish divided by total time durationNovel Tank Test
(1-1 for 0–11 min and 1-2 for 15–31 min)
1-1-2 & 1-2-2Freezing time movement ratio (%)Total percentage of time when fish’s speed less than 1 cm s−1
1-1-3 & 1-2-3Swimming time movement ratio (%)Total percentage of time when fish’s speed is between 1 and 10 cm s−1
1-1-4 & 1-2-4Rapid movement ratio (%)Total percentage of time when fish’s speed more than 10 cm s−1
1-1-5 & 1-2-5Time in top duration (%)Total time spent in the top portion of the novel tank in percentage
1-1-6 & 1-2-6Number of entries to the topTotal times fish enter the upper half of the tank
1-1-7 & 1-2-7Latency to enter the top (s)The amount of time it takes the fish to cross into the upper half of the tank
1-1-8 & 1-2-8Total distance traveled in the top (cm)Total distance traveled in the top portion of the novel tank
1-1-9 & 1-2-9Thigmotaxis (cm)The average distance of the fish from the center of the tank
2-1Mirror biting time percentage (%)Total percentage of time when fish stayed in the mirror biting zoneMirror Biting Test
2-2Longest duration in the mirror side percentage (%)Total percentage of fish longest duration stayed in the mirror biting zone
3-1Approaching predator time percentage (%)Total percentage of time when fish stayed in the approaching predator zonePredator Avoidance Test
3-2Average distance to the predator’s separator (cm)Average distance of fish to the predator’s separator
4-1Conspecific interaction time percentage (%)Total percentage of time when fish stayed in the conspecific interaction zoneSocial Interaction Test
4-2Average distance to the conspecific’s separator (cm)Average distance of fish to the conspecific’s separator
4-3Longest conspecific interaction percentage (%)Total percentage of fish longest duration stayed in the conspecific interaction
5-1Average inter-fish distance (cm)Average distance between the body center of every member of the shoalShoaling Test
5-2Average shoal area (cm2)Average size of the shoal
5-3Average nearest neighbor distance (cm)Distance for the body center of each fish to the closest neighboring fish
5-4Average farthest neighbor distance (cm)Distance for the body center of each fish to the farthest neighboring fish
6-1-1 & 6-2-1Average speed (cm s−1)Total distance traveled by fish divided by total time durationCircadian Locomotor Activity Rhythm
(6-1 for day cycle and 6-2 for night cycle)
6-1-2 & 6-2-2Average Angular Velocity (°/s)Total turning angle of fish divided by total time duration
6-1-3 & 6-2-3Meandering (°/m)Total turning angle of fish divided by total distance traveled by fish
6-1-4 & 6-2-4Freezing time movement ratio (%)Total percentage of time when fish’s speed less than 1 cm s−1
6-1-5 & 6-2-5Swimming time movement ratio (%)Total percentage of time when fish’s speed is between 1 and 10 cm s−1
6-1-6 & 6-2-6Rapid movement ratio (%)Total percentage of time when fish’s speed more than 10 cm s−1
Table A5. Summary table of each strain coefficient of variation in each behavior endpoint.
Table A5. Summary table of each strain coefficient of variation in each behavior endpoint.
IndexCoefficient of Variation (%)
Danio rerioOryzias dancenaOryzias latipesOryzias woworaeOryzias sinensis
1-1-130.1876.4381.7870.8364.87
1-1-285.2867.3550.1666.5450.09
1-1-311.7250.7539.6240.1251.07
1-1-4119.26183.56300.39127.78151.28
1-1-594.04125.63163.79130.8369.09
1-1-682.26122.4156.09125.4290.9
1-1-7150.18102.7874.590.45180.32
1-1-899.42127.26173.59139.0681.21
1-1-933.233.0233.8131.9237.59
1-2-136.264.8257.8143.0561.88
1-2-271.5980.0870.0772.6757.6
1-2-38.8347.6632.2619.3646.96
1-2-4153.56149.89223.11107.58173.84
1-2-561.4592.2183.5169.7162.6
1-2-678.0591.7189.1367.0388.95
1-2-7164.92136.08142.32156.62207.24
1-2-870.0396.2106.174.8870.96
1-2-935.5836.7430.8731.2528.83
2-147.4260.3164.6936.6242.34
2-283.05104.588.4471.9179.35
3-113796.47128.2160.5226
3-213.8432.4528.2325.1721.64
4-133.8645.181.9543.768.01
4-255.2445.4631.0372.4555.06
4-396.13109.6131.694.3192.97
5-120.7347.4831.4627.8739.97
5-246.0266.5356.8554.2376.82
5-320.7351.5228.6827.4944.87
5-419.3334.872725.2528.02
6-1-143.7847.2564.574561.46
6-1-244.3345.9651.8650.7171.5
6-1-371.446.451.8348.3558.84
6-1-4205.747.1252.7952.0336.53
6-1-542.5137.8448.9138.673.53
6-1-649.05265293.3177.9171.9
6-2-170.457.5869.3576.3842.71
6-2-262.5247.0349.1244.5350.04
6-2-364.6351.8651.7464.8354.78
6-2-468.432.228.9647.3111.73
6-2-538.8877.8191.6273.2572.81
6-2-6170.9217.7225.2217.1220.2
Average70.5268381.7709889.9095174.1607379.91122

References

  1. Hilgers, L.; Schwarzer, J. The untapped potential of medaka and its wild relatives. eLife 2019, 8, 8. [Google Scholar] [CrossRef] [PubMed]
  2. Nakao, R.; Kitagawa, T. Differences in the behavior and ecology of wild type medaka (Oryzias latipes complex) and an orange commercial variety (himedaka). J. Exp. Zool. Part A Ecol. Genet. Physiol. 2015, 323, 349–358. [Google Scholar] [CrossRef] [PubMed]
  3. Ansai, S.; Hosokawa, H.; Maegawa, S.; Kinoshita, M. Chronic fluoxetine treatment induces anxiolytic responses and altered social behaviors in medaka, Oryzias latipes. Behav. Brain Res. 2016, 303, 126–136. [Google Scholar] [CrossRef] [PubMed]
  4. Shima, A.; Mitani, H. Medaka as a research organism: Past, present and future. Mech. Dev. 2004, 121, 599–604. [Google Scholar] [CrossRef] [PubMed]
  5. Lucon-Xiccato, T.; Conti, F.; Loosli, F.; Foulkes, N.S.; Bertolucci, C. Development of Open-Field Behaviour in the Medaka, Oryzias latipes. Biology 2020, 9, 389. [Google Scholar] [CrossRef] [PubMed]
  6. Khalil, F.; Kang, I.J.; Undap, S.; Tasmin, R.; Qiu, X.; Shimasaki, Y.; Oshima, Y. Alterations in social behavior of Japanese medaka (Oryzias latipes) in response to sublethal chlorpyrifos exposure. Chemosphere 2013, 92, 125–130. [Google Scholar] [CrossRef]
  7. OECD. Test No. 203: Fish, Acute Toxicity Test; OECD Guidelines for the Testing of Chemicals, Section 2; OECD: Paris, France, 1992. [Google Scholar]
  8. Lee, W.; Yang, K.-L. Using medaka embryos as a model system to study biological effects of the electromagnetic fields on development and behavior. Ecotoxicol. Environ. Saf. 2014, 108, 187–194. [Google Scholar] [CrossRef]
  9. Kinoshita, M.; Murata, K.; Naruse, K.; Tanaka, M. Medaka: Biology, Management, and Experimental Protocols; John Wiley & Sons: Honoken, NJ, USA, 2009. [Google Scholar]
  10. Hong, X.; Zha, J. Fish behavior: A promising model for aquatic toxicology research. Sci. Total. Environ. 2019, 686, 311–321. [Google Scholar] [CrossRef]
  11. Nakayama, K.; Oshima, Y.; Hiramatsu, K.; Shimasaki, Y.; Honjo, T. Effects of Polychlorinated Biphenyls on the Schooling Behavior of Japanese Medaka (Oryzias Latipes). Environ. Toxicol. Chem. 2005, 24, 2588–2593. [Google Scholar] [CrossRef]
  12. Ju, Z.; Wells, M.C.; Heater, S.J.; Walter, R.B. Multiple tissue gene expression analyses in Japanese medaka (Oryzias latipes) exposed to hypoxia. Comp. Biochem. Physiol. Part C Toxicol. Pharmacol. 2007, 145, 134–144. [Google Scholar] [CrossRef]
  13. Naruse, K.; Tanaka, M.; Takeda, H. Medaka: A Model for Organogenesis, Human Disease, and Evolution; Springer: Japan, Tokyo, 2011. [Google Scholar]
  14. Naruse, K.; Tanaka, M.; Mita, K.; Shima, A.; Postlethwait, J.; Mitani, H. A Medaka Gene Map: The Trace of Ancestral Vertebrate Proto-Chromosomes Revealed by Comparative Gene Mapping. Genome Res. 2004, 14, 820–828. [Google Scholar] [CrossRef] [Green Version]
  15. Ishikawa, Y.; Yoshimoto, M.; Ito, H. A brain atlas of a wild-type inbred strain of the medaka, Oryzias latipes. Fish Biol. J. MEDAKA 1999, 10, 1–26. [Google Scholar]
  16. Naruse, K.; Fukamachi, S.; Mitani, H.; Kondo, M.; Matsuoka, T.; Kondo, S.; Hanamura, N.; Morita, Y.; Hasegawa, K.; Nishigaki, R.; et al. A Detailed Linkage Map of Medaka, Oryzias latipes: Comparative Genomics and Genome Evolution. Genetics 2000, 154, 1773–1784. [Google Scholar] [CrossRef]
  17. Khorasani, M.Z.; Hennig, S.; Imre, G.; Asakawa, S.; Palczewski, S.; Berger, A.; Hori, H.; Naruse, K.; Mitani, H.; Shima, A.; et al. A first generation physical map of the medaka genome in BACs essential for positional cloning and clone-by-clone based genomic sequencing. Mech. Dev. 2004, 121, 903–913. [Google Scholar] [CrossRef]
  18. Tsuboko-Ishii, S.; Kimura, T.; Shinya, M.; Suehiro, Y.; Okuyama, T.; Shimada, A.; Takeda, H.; Naruse, K.; Kubo, T.; Takeuchi, H. Genetic Control of Startle Behavior in Medaka Fish. PLoS ONE 2014, 9, e112527. [Google Scholar] [CrossRef] [Green Version]
  19. Matsunaga, W.; Watanabe, E. Habituation of medaka (Oryzias latipes) demonstrated by open-field testing. Behav. Process. 2010, 85, 142–150. [Google Scholar] [CrossRef]
  20. Réale, D.; Reader, S.M.; Sol, D.; McDougall, P.T.; Dingemanse, N.J. Integrating animal temperament within ecology and evolution. Biol. Rev. 2007, 82, 291–318. [Google Scholar] [CrossRef] [Green Version]
  21. Bouchard, T.J.; Loehlin, J.C. Genes, evolution, and personality. Behav. Genet. 2001, 31, 243–273. [Google Scholar] [CrossRef]
  22. Smith, B.R.; Blumstein, D. Fitness consequences of personality: A meta-analysis. Behav. Ecol. 2008, 19, 448–455. [Google Scholar] [CrossRef] [Green Version]
  23. Carter, A.J.; Feeney, W.E.; Marshall, H.H.; Cowlishaw, G.; Heinsohn, R. Animal personality: What are behavioural ecologists measuring? Biol. Rev. 2013, 88, 465–475. [Google Scholar] [CrossRef] [Green Version]
  24. Sih, A.; Bell, A.; Johnson, J. Behavioral syndromes: An ecological and evolutionary overview. Trends Ecol. Evol. 2004, 19, 372–378. [Google Scholar] [CrossRef] [Green Version]
  25. Koolhaas, J.; de Boer, S.; Coppens, C.; Buwalda, B. Neuroendocrinology of coping styles: Towards understanding the biology of individual variation. Front. Neuroendocrinol. 2010, 31, 307–321. [Google Scholar] [CrossRef]
  26. Saunders, B.T.; Robinson, E.T. Individual Variation in the Motivational Properties of Cocaine. Neuropsychopharmacology 2011, 36, 1668–1676. [Google Scholar] [CrossRef]
  27. Conrad, J.L.; Weinersmith, K.L.; Brodin, T.; Saltz, J.B.; Sih, A. Behavioural syndromes in fishes: A review with implications for ecology and fisheries management. J. Fish Biol. 2011, 78, 395–435. [Google Scholar] [CrossRef]
  28. Demin, K.A.; Lakstygal, A.M.; Alekseeva, P.A.; Sysoev, M.; de Abreu, M.S.; Alpyshov, E.T.; Serikuly, N.; Wang, D.; Wang, M.; Tang, Z.; et al. The role of intraspecies variation in fish neurobehavioral and neuropharmacological phenotypes in aquatic models. Aquat. Toxicol. 2019, 210, 44–55. [Google Scholar] [CrossRef]
  29. Audira, G.; Siregar, P.; Strungaru, S.-A.; Huang, J.-C.; Hsiao, C.-D. Which Zebrafish Strains Are More Suitable to Perform Behavioral Studies? A Comprehensive Comparison by Phenomic Approach. Biology 2020, 9, 200. [Google Scholar] [CrossRef]
  30. Hyodo-Taguchi, Y. Inbred Strains of the Medaka and their Characteristics; Egami, N., Yamagami, K., Shima, A., Eds.; Tokyo University Press: Tokyo, Japan, 1990; pp. 129–142. [Google Scholar]
  31. Ansai, S.; Mochida, K.; Fujimoto, S.; Mokodongan, D.F.; Sumarto, B.K.A.; Masengi, K.W.A.; Hadiaty, R.K.; Nagano, A.J.; Toyoda, A.; Naruse, K.; et al. Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes. Nat. Commun. 2021, 12, 1–13. [Google Scholar] [CrossRef] [PubMed]
  32. Blaser, R.E.; Rosemberg, D.B. Measures of Anxiety in Zebrafish (Danio rerio): Dissociation of Black/White Preference and Novel Tank Test. PLoS ONE 2012, 7, e36931. [Google Scholar] [CrossRef] [Green Version]
  33. Levin, E.D.; Bencan, Z.; Cerutti, D.T. Anxiolytic effects of nicotine in zebrafish. Physiol. Behav. 2007, 90, 54–58. [Google Scholar] [CrossRef]
  34. Cachat, J.; Stewart, A.; Grossman, L.; Gaikwad, S.; Kadri, F.; Chung, K.M.; Wu, N.; Wong, K.; Roy, S.; Suciu, C.; et al. Measuring behavioral and endocrine responses to novelty stress in adult zebrafish. Nat. Protoc. 2010, 5, 1786–1799. [Google Scholar] [CrossRef]
  35. Pham, M.; Raymond, J.; Hester, J.; Kyzar, E.; Gaikwad, S.; Bruce, I.; Fryar, C.; Chanin, S.; Enriquez, J.; Bagawandoss, S. Assessing social behavior phenotypes in adult zebrafish: Shoaling, social preference, and mirror biting tests. Zebrafish Protoc. Neurobehav. Res. 2012, 66, 231–246. [Google Scholar] [CrossRef]
  36. Ahmed, O.; Séguin, D.; Gerlai, R. An automated predator avoidance task in zebrafish. Behav. Brain Res. 2011, 216, 166–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Audira, G.; Sampurna, B.P.; Juniardi, S.; Liang, S.-T.; Lai, Y.-H.; Hsiao, C.-D. A Versatile Setup for Measuring Multiple Behavior Endpoints in Zebrafish. Inventions 2018, 3, 75. [Google Scholar] [CrossRef] [Green Version]
  38. Engeszer, R.E.; Da Barbiano, L.A.; Ryan, M.J.; Parichy, D.M. Timing and plasticity of shoaling behaviour in the zebrafish, Danio rerio. Anim. Behav. 2007, 74, 1269–1275. [Google Scholar] [CrossRef] [Green Version]
  39. Miller, N.; Gerlai, R. Quantification of shoaling behaviour in zebrafish (Danio rerio). Behav. Brain Res. 2007, 184, 157–166. [Google Scholar] [CrossRef]
  40. Gwinner, F. Testosterone Induces “Splitting” of Circadian Locomotor Activity Rhythms in Birds. Science 1974, 185, 72–74. [Google Scholar] [CrossRef]
  41. Baker, M.R.; Goodman, A.C.; Santo, J.B.; Wong, R.Y. Repeatability and reliability of exploratory behavior in proactive and reactive zebrafish, Danio rerio. Sci. Rep. 2018, 8, 1–9. [Google Scholar] [CrossRef]
  42. Wong, K.; Elegante, M.; Bartels, B.; Elkhayat, S.; Tien, D.; Roy, S.; Goodspeed, J.; Suciu, C.; Tan, J.; Grimes, C.; et al. Analyzing habituation responses to novelty in zebrafish (Danio rerio). Behav. Brain Res. 2010, 208, 450–457. [Google Scholar] [CrossRef]
  43. Sataa, N.S.A.M.; Bakar, N.A.; Hodin, N.A.S.; Ramlan, N.F.; Ibrahim, W.N.W. Behavioral Responses of Javanese Medaka (Oryzias Javanicus) Versus Zebrafish (Danio Rerio) in Open Field Test. Res. Sq. 2020, 1–24. [Google Scholar] [CrossRef]
  44. Schnörr, S.; Steenbergen, P.; Richardson, M.; Champagne, D. Measuring thigmotaxis in larval zebrafish. Behav. Brain Res. 2012, 228, 367–374. [Google Scholar] [CrossRef]
  45. Peitsaro, N.; Kaslin, J.; Anichtchik, O.V.; Panula, P. Modulation of the histaminergic system and behaviour by α-fluoromethylhistidine in zebrafish. J. Neurochem. 2004, 86, 432–441. [Google Scholar] [CrossRef]
  46. Nakayasu, T.; Yasugi, M.; Shiraishi, S.; Uchida, S.; Watanabe, E. Three-dimensional computer graphic animations for studying social approach behaviour in medaka fish: Effects of systematic manipulation of morphological and motion cues. PLoS ONE 2017, 12, e0175059. [Google Scholar] [CrossRef]
  47. Tsubokawa, T.; Saito, K.; Kawano, H.; Kawamura, K.; Shinozuka, K.; Watanabe, S. Pharmacological effects on mirror approaching behavior and neurochemical aspects of the telencephalon in the fish, medaka (Oryzias latipes). Soc. Neurosci. 2009, 4, 276–286. [Google Scholar] [CrossRef]
  48. Otsuka, A.; Inahata, M.; Shimomura, Y.; Kagawa, N. Physiological changes in response to social isolation in male medaka fish. Fish. Sci. 2020, 86, 1–7. [Google Scholar] [CrossRef]
  49. De Abreu, M.S.; Koakoski, G.; Ferreira, D.; Oliveira, T.A.; Da Rosa, J.G.S.; Gusso, D.; Giacomini, A.C.V.; Piato, A.L.; Gil Barcellos, L.J. Diazepam and Fluoxetine Decrease the Stress Response in Zebrafish. PLoS ONE 2014, 9, e103232. [Google Scholar] [CrossRef]
  50. Dickinson, B.C.; Chang, C.J. Chemistry and biology of reactive oxygen species in signaling or stress responses. Nat. Chem. Biol. 2011, 7, 504–511. [Google Scholar] [CrossRef] [Green Version]
  51. Mocelin, R.N.; Herrmann, A.P.; Marcon, M.; Rambo, C.L.; Rohden, A.; Bevilaqua, F.; de Abreu, M.S.; Zanatta, L.; Elisabetsky, E.; Barcellos, L.J.; et al. N-acetylcysteine prevents stress-induced anxiety behavior in zebrafish. Pharmacol. Biochem. Behav. 2015, 139, 121–126. [Google Scholar] [CrossRef]
  52. Polavarapu, R.; Spitz, D.R.; Sim, J.E.; Follansbee, M.H.; Oberley, L.W.; Rahemtulla, A.; Nanji, A.A. Increased lipid peroxidation and impaired antioxidant enzyme function is associated with pathological liver injury in experimental alcoholic liver disease in rats fed diets high in corn oil and fish oil. Hepatology 1998, 27, 1317–1323. [Google Scholar] [CrossRef]
  53. Abdel-Salam, O.M.E.; Youness, E.R.; Morsy, F.A.; Mahfouz, M.M.; Kenawy, S.A. Study of the effect of antidepressant drugs and donepezil on aluminum-induced memory impairment and biochemical alterations in rats. Comp. Haematol. Int. 2015, 24, 847–860. [Google Scholar] [CrossRef]
  54. Lillesaar, C. The serotonergic system in fish. J. Chem. Neuroanat. 2011, 41, 294–308. [Google Scholar] [CrossRef]
  55. Mathuru, A.S. Conspecific injury raises an alarm in medaka. Sci. Rep. 2016, 6, 36615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Reznick, D.N.; Ghalambor, C.K.; Crooks, K. Experimental studies of evolution in guppies: A model for understanding the evolutionary consequences of predator removal in natural communities. Mol. Ecol. 2008, 17, 97–107. [Google Scholar] [CrossRef] [PubMed]
  57. Eisenberg, M.; Dudai, Y. Reconsolidation of fresh, remote, and extinguished fear memory in Medaka: Old fears don’t die. Eur. J. Neurosci. 2004, 20, 3397–3403. [Google Scholar] [CrossRef] [PubMed]
  58. Ladu, F.; Bartolini, T.; Panitz, S.G.; Chiarotti, F.; Butail, S.; Macrí, S.; Porfiri, M. Live Predators, Robots, and Computer-Animated Images Elicit Differential Avoidance Responses in Zebrafish. Zebrafish 2015, 12, 205–214. [Google Scholar] [CrossRef]
  59. Bass, S.L.; Gerlai, R. Zebrafish (Danio rerio) responds differentially to stimulus fish: The effects of sympatric and allopatric predators and harmless fish. Behav. Brain Res. 2008, 186, 107–117. [Google Scholar] [CrossRef]
  60. Mosienko, V.; Bert, B.; Beis, D.; Matthes, S.; Fink, H.; Bader, M.; Alenina, N. Exaggerated aggression and decreased anxiety in mice deficient in brain serotonin. Transl. Psychiatry 2012, 2, e122. [Google Scholar] [CrossRef] [Green Version]
  61. Ball, G.F.; Balthazart, J.; McCarthy, M.M. Is it useful to view the brain as a secondary sexual characteristic? Neurosci. Biobehav. Rev. 2014, 46, 628–638. [Google Scholar] [CrossRef]
  62. Mason, W. Ontogeny of social behavior. In Social Behavior and Communication; Springer: Berlin/Heidelberg, Germany, 1979; pp. 1–28. [Google Scholar]
  63. Isoe, Y.; Konagaya, Y.; Yokoi, S.; Kubo, T.; Takeuchi, H. Ontogeny and Sexual Differences in Swimming Proximity to Conspecifics in Response to Visual Cues in Medaka Fish. Zool. Sci. 2016, 33, 246–254. [Google Scholar] [CrossRef]
  64. Shibai, A.; Arimoto, T.; Yoshinaga, T.; Tsuchizawa, Y.; Khureltulga, D.; Brown, Z.P.; Kakizuka, T.; Hosoda, K. Attraction of posture and motion-trajectory elements of conspecific biological motion in medaka fish. Sci. Rep. 2018, 8, 8589. [Google Scholar] [CrossRef] [Green Version]
  65. Saverino, C.; Gerlai, R. The social zebrafish: Behavioral responses to conspecific, heterospecific, and computer animated fish. Behav. Brain Res. 2008, 191, 77–87. [Google Scholar] [CrossRef] [Green Version]
  66. Imada, H.; Hoki, M.; Suehiro, Y.; Okuyama, T.; Kurabayashi, D.; Shimada, A.; Naruse, K.; Takeda, H.; Kubo, T.; Takeuchi, H. Coordinated and Cohesive Movement of Two Small Conspecific Fish Induced by Eliciting a Simultaneous Optomotor Response. PLoS ONE 2010, 5, e11248. [Google Scholar] [CrossRef]
  67. Nakayasu, T.; Watanabe, E. Biological motion stimuli are attractive to medaka fish. Anim. Cogn. 2013, 17, 559–575. [Google Scholar] [CrossRef] [Green Version]
  68. Winberg, S.; Nilsson, G.E. Roles of brain monoamine neurotransmitters in agonistic behaviour and stress reactions, with particular reference to fish. Comp. Biochem. Physiol. Part C Pharmacol. Toxicol. Endocrinol. 1993, 106, 597–614. [Google Scholar] [CrossRef]
  69. Øverli, Ø.; Harris, C.A.; Winberg, S. Short-Term Effects of Fights for Social Dominance and the Establishment of Dominant-Subordinate Relationships on Brain Monoamines and Cortisol in Rainbow Trout. Brain Behav. Evol. 1999, 54, 263–275. [Google Scholar] [CrossRef]
  70. Hare, A.J.; Waheed, A.; Hare, J.F.; Anderson, W.G. Cortisol and catecholamine responses to social context and a chemical alarm signal in juvenile lake sturgeon, Acipenser fulvescens. Can. J. Zool. 2015, 93, 605–613. [Google Scholar] [CrossRef]
  71. Shams, S.; Amlani, S.; Buske, C.; Chatterjee, D.; Gerlai, R. Developmental social isolation affects adult behavior, social interaction, and dopamine metabolite levels in zebrafish. Dev. Psychobiol. 2018, 60, 43–56. [Google Scholar] [CrossRef]
  72. Saif, M.; Chatterjee, D.; Buske, C.; Gerlai, R. Sight of conspecific images induces changes in neurochemistry in zebrafish. Behav. Brain Res. 2013, 243, 294–299. [Google Scholar] [CrossRef]
  73. Egami, N. Effect of artificial photoperiodicity on time of oviposition in the fish, Oryzias latipes. Anim. Behav. 1954, 27, 57–62. [Google Scholar]
  74. Ueda, M.; Oishi, T. Circadian oviposition rhythm and locomotor activity in the medaka, Oryzias latipes. J. Interdiscip. Cycle Res. 1982, 13, 97–104. [Google Scholar] [CrossRef]
  75. Karigo, T.; Kanda, S.; Takahashi, A.; Abe, H.; Okubo, K.; Oka, Y. Time-of-Day-Dependent Changes in GnRH1 Neuronal Activities and Gonadotropin mRNA Expression in a Daily Spawning Fish, Medaka. Endocrinology 2012, 153, 3394–3404. [Google Scholar] [CrossRef] [Green Version]
  76. Ali, M.A. Rhythms in Fishes; Springer: Berlin/Heidelberg, Germany, 1992. [Google Scholar]
  77. Weber, D.N.; Spieler, R.E. Effects of the light-dark cycle and scheduled feeding on behavioral and reproductive rhythms of the cyprinodont fish, Medaka, Oryzias latipes. Cell. Mol. Life Sci. 1987, 43, 621–624. [Google Scholar] [CrossRef]
  78. Thi, N.H.B.; Audira, G.; Siregar, P.; Liang, S.-T.; Huang, J.-C.; Hsiao, C.-D.; Thi, N.A.N.; Thi, N.B.; Thi, N.N. Chronic Exposure to Low Concentration Lead Chloride-Induced Anxiety and Loss of Aggression and Memory in Zebrafish. Int. J. Mol. Sci. 2020, 21, 1844. [Google Scholar] [CrossRef] [Green Version]
  79. Roshchina, V.V. Evolutionary Considerations of Neurotransmitters in Microbial, Plant, and Animal Cells. Microb. Endocrinol. 2010, 17–52. [Google Scholar] [CrossRef]
  80. Van Ginneken, V.; Nieveen, M.; Van Eersel, R.; Thillart, G.V.D.; Addink, A. Neurotransmitter levels and energy status in brain of fish species with and without the survival strategy of metabolic depression. Comp. Biochem. Physiol. Part A Physiol. 1996, 114, 189–196. [Google Scholar] [CrossRef]
  81. Porteus, C.S.; Brink, D.L.; Milsom, W.K. Neurotransmitter profiles in fish gills: Putative gill oxygen chemoreceptors. Respir. Physiol. Neurobiol. 2012, 184, 316–325. [Google Scholar] [CrossRef]
  82. Reid, S.G.; Bernier, N.J.; Perry, S.F. The adrenergic stress response in fish: Control of catecholamine storage and release. Comp. Biochem. Physiol. Part C Pharmacol. Toxicol. Endocrinol. 1998, 120, 1–27. [Google Scholar] [CrossRef] [Green Version]
  83. Aprison, M.; Shank, R.; Davidoff, R.; Werman, R. The distribution of glycine, a neurotransmitter suspect in the central nervous system of several vertebrate species. Life Sci. 1968, 7, 583–590. [Google Scholar] [CrossRef]
  84. Wilkinson, G.S.; Breden, F.; Mank, J.; Ritchie, M.G.; Higginson, A.D.; Radwan, J.; Jaquiery, J.; Salzburger, W.; Arriero, E.; Barribeau, S.M.; et al. The locus of sexual selection: Moving sexual selection studies into the post-genomics era. J. Evol. Biol. 2015, 28, 739–755. [Google Scholar] [CrossRef] [Green Version]
  85. Yoon, M.; Kim, K.-Y.; Bang, I.-C.; Nam, Y.K.; Kim, D.S. Complete mitogenome sequence of the Chinese medaka Oryzias sinensis (Teleostei: Beloniformes) and its phylogenetic analysis. Genes Genom. 2011, 33, 307–312. [Google Scholar] [CrossRef]
  86. Sakamoto, T.; Kozaka, T.; Takahashi, A.; Kawauchi, H.; Ando, M. Medaka (Oryzias latipes) as a model for hypoosmoregulation of euryhaline fishes. Aquaculture 2001, 193, 347–354. [Google Scholar] [CrossRef]
  87. Kim, Y.K.; Lee, S.Y.; Kim, B.S.; Kim, D.S.; Nam, Y.K. Isolation and mRNA expression analysis of aquaporin isoforms in marine medaka Oryzias dancena, a euryhaline teleost. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2014, 171, 1–8. [Google Scholar] [CrossRef] [PubMed]
  88. Parenti, L.R. A phylogenetic analysis and taxonomic revision of ricefishes, Oryzias and relatives (Beloniformes, Adrianichthyidae). Zool. J. Linn. Soc. 2008, 154, 494–610. [Google Scholar] [CrossRef] [Green Version]
  89. Juo, J.-J.; Kang, C.-K.; Yang, W.-K.; Yang, S.-Y.; Lee, T.-H. A Stenohaline Medaka, Oryzias woworae, Increases Expression of Gill Na, K-ATPase and Na, K, 2Cl-Cotransporter 1 to Tolerate Osmotic Stress. Zool. Sci. 2016, 33, 414–425. [Google Scholar] [CrossRef] [PubMed]
  90. Hsu, H.-H.; Lin, L.-Y.; Tseng, Y.-C.; Horng, J.-L.; Hwang, P.-P. A new model for fish ion regulation: Identification of ionocytes in freshwater- and seawater-acclimated medaka (Oryzias latipes). Cell Tissue Res. 2014, 357, 225–243. [Google Scholar] [CrossRef] [PubMed]
  91. Chen, Y. Taxonomy and distribution of the genus Oryzias in Yunnan, China. Acta Taxonom. Sin. 1989, 14, 239–246. [Google Scholar]
  92. Takehana, Y.; Naruse, K.; Sakaizumi, M. Molecular phylogeny of the medaka fishes genus Oryzias (Beloniformes: Adrianichthyidae) based on nuclear and mitochondrial DNA sequences. Mol. Phylogenet. Evol. 2005, 36, 417–428. [Google Scholar] [CrossRef]
  93. Kang, C.-K.; Tsai, S.-C.; Lee, T.-H.; Hwang, P.-P. Differential expression of branchial Na+/K+-ATPase of two medaka species, Oryzias latipes and Oryzias dancena, with different salinity tolerances acclimated to fresh water, brackish water and seawater. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2008, 151, 566–575. [Google Scholar] [CrossRef]
  94. Cui, L.; Fan, M.; Belanger, S.; Li, J.; Wang, X.; Fan, B.; Li, W.; Gao, X.; Chen, J.; Liu, Z. Oryzias sinensis, a new model organism in the application of eco-toxicity and water quality criteria (WQC). Chemosphere 2020, 261, 127813. [Google Scholar] [CrossRef]
  95. Myosho, T.; Takahashi, H.; Yoshida, K.; Sato, T.; Hamaguchi, S.; Sakamoto, T.; Sakaizumi, M. Hyperosmotic tolerance of adult fish and early embryos are determined by discrete, single loci in the genus Oryzias. Sci. Rep. 2018, 8, 6897. [Google Scholar] [CrossRef]
  96. Chen, W.-J.; Orti, G.; Meyer, A. Novel evolutionary relationship among four fish model systems. Trends Genet. 2004, 20, 424–431. [Google Scholar] [CrossRef] [Green Version]
  97. Colón-Cruz, L.; Kristofco, L.; Crooke-Rosado, J.; Acevedo, A.; Torrado, A.; Brooks, B.W.; Sosa, M.A.; Behra, M. Alterations of larval photo-dependent swimming responses (PDR): New endpoints for rapid and diagnostic screening of aquatic contamination. Ecotoxicol. Environ. Saf. 2018, 147, 670–680. [Google Scholar] [CrossRef]
  98. Spence, R.; Gerlach, G.; Lawrence, C.; Smith, C. The behaviour and ecology of the zebrafish, Danio rerio. Biol. Rev. 2007, 83, 13–34. [Google Scholar] [CrossRef]
  99. Avdesh, A.; Chen, M.; Martin-Iverson, M.T.; Mondal, A.; Ong, D.; Rainey-Smith, S.; Taddei, K.; Lardelli, M.; Groth, D.M.; Verdile, G.; et al. Regular Care and Maintenance of a Zebrafish (Danio rerio) Laboratory: An Introduction. J. Vis. Exp. 2012, e4196. [Google Scholar] [CrossRef] [Green Version]
  100. Audira, G.; Sampurna, B.P.; Juniardi, S.; Liang, S.-T.; Lai, Y.-H.; Han, L.; Hsiao, C.-D. Establishing simple image-based methods and cost-effective instrument for toxicity assessment on circadian rhythm dysregulation in fish. Biol. Open 2019, 8, 041871. [Google Scholar] [CrossRef] [Green Version]
  101. Pérez-Escudero, A.; Vicente-Page, J.; Hinz, R.C.; Arganda, S.; De Polavieja, G.G. idTracker: Tracking individuals in a group by automatic identification of unmarked animals. Nat. Methods 2014, 11, 743–748. [Google Scholar] [CrossRef]
  102. Adams, D.C.; Anthony, C.D. Using randomization techniques to analyse behavioural data. Anim. Behav. 1996, 51, 733–738. [Google Scholar] [CrossRef] [Green Version]
  103. Metsalu, T.; Vilo, J. ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015, 43, W566–W570. [Google Scholar] [CrossRef]
  104. Furutani-Seiki, M.; Wittbrodt, J. Medaka and zebrafish, an evolutionary twin study. Mech. Dev. 2004, 121, 629–637. [Google Scholar] [CrossRef]
  105. Wittbrodt, J.; Shima, A.; Schartl, M. Medaka—A model organism from the far East. Nat. Rev. Genet. 2002, 3, 53–64. [Google Scholar] [CrossRef]
  106. Sasado, T.; Tanaka, M.; Kobayashi, K.; Sato, T.; Sakaizumi, M.; Naruse, K. The National BioResource Project Medaka (NBRP Medaka): An Integrated Bioresource for Biological and Biomedical Sciences. Exp. Anim. 2010, 59, 13–23. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Novel tank behavior endpoints comparison between each tested fish species (medaka and zebrafish). (A) Average speed, (B) freezing time movement ratio, (C) swimming time movement ratio, (D) rapid movement time ratio, (E) time in top duration, (F) number of entries to the top, (G) latency to enter the top, (H) total distance traveled in the top, and (I) average distance to the center of the tank were analyzed. The data were analyzed by the two-way ANOVA test with Geisser–Greenhouse correction continued with uncorrected Fisher’s LSD test. Different letters (a, b, c, d) on the error bars represent a significant statistical difference (p < 0.05). The data are expressed as the median with interquartile range (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 22 for O. sinensis).
Figure 1. Novel tank behavior endpoints comparison between each tested fish species (medaka and zebrafish). (A) Average speed, (B) freezing time movement ratio, (C) swimming time movement ratio, (D) rapid movement time ratio, (E) time in top duration, (F) number of entries to the top, (G) latency to enter the top, (H) total distance traveled in the top, and (I) average distance to the center of the tank were analyzed. The data were analyzed by the two-way ANOVA test with Geisser–Greenhouse correction continued with uncorrected Fisher’s LSD test. Different letters (a, b, c, d) on the error bars represent a significant statistical difference (p < 0.05). The data are expressed as the median with interquartile range (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 22 for O. sinensis).
Ijms 22 05686 g001
Figure 2. Mirror biting and predator avoidance behavior comparisons between each tested fish species (medaka and zebrafish). (A) Mirror biting time percentage and (B) longest duration on the mirror side percentage were analyzed in the mirror biting assay (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 22 for O. sinensis). (C) Approaching predator time percentage and (D) average distance to the predator’s separator were analyzed in the predator avoidance assay (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 16 for O. sinensis). The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c) on the error bars represent a significant statistical difference (p < 0.05). The violin plot’s median and interquartile were labeled with the bold line colored with cyan and yellow.
Figure 2. Mirror biting and predator avoidance behavior comparisons between each tested fish species (medaka and zebrafish). (A) Mirror biting time percentage and (B) longest duration on the mirror side percentage were analyzed in the mirror biting assay (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 22 for O. sinensis). (C) Approaching predator time percentage and (D) average distance to the predator’s separator were analyzed in the predator avoidance assay (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 16 for O. sinensis). The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c) on the error bars represent a significant statistical difference (p < 0.05). The violin plot’s median and interquartile were labeled with the bold line colored with cyan and yellow.
Ijms 22 05686 g002
Figure 3. Conspecific social interaction behavior comparisons between each tested fish species (medaka and zebrafish). (A) Conspecific interaction time percentage, (B) average distance to the conspecifics separator, and (C) longest conspecific interaction percentage were analyzed. The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c, d) on the error bars represent a significant statistical difference (p < 0.05) (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 21 for O. sinensis). The median and interquartile for the violin plot were labeled with the bold line colored with cyan and yellow.
Figure 3. Conspecific social interaction behavior comparisons between each tested fish species (medaka and zebrafish). (A) Conspecific interaction time percentage, (B) average distance to the conspecifics separator, and (C) longest conspecific interaction percentage were analyzed. The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c, d) on the error bars represent a significant statistical difference (p < 0.05) (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 21 for O. sinensis). The median and interquartile for the violin plot were labeled with the bold line colored with cyan and yellow.
Ijms 22 05686 g003
Figure 4. Shoaling behavior comparisons between each tested fish species (medaka and zebrafish). (A) Average inter-fish distance, (B) average shoal area, (C) average nearest neighbor distance, and (D) average farthest neighbor distance were analyzed. Groups of three fish were tested for shoaling behavior. The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b) on the error bars represent a significant statistical difference (p < 0.05) (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 21 for O. sinensis). The median and interquartile for the violin plot were labeled with the bold line colored with cyan and yellow.
Figure 4. Shoaling behavior comparisons between each tested fish species (medaka and zebrafish). (A) Average inter-fish distance, (B) average shoal area, (C) average nearest neighbor distance, and (D) average farthest neighbor distance were analyzed. Groups of three fish were tested for shoaling behavior. The data were analyzed by the Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b) on the error bars represent a significant statistical difference (p < 0.05) (n = 30 for zebrafish, O. dancena, O. latipes, and O. woworae; n = 21 for O. sinensis). The median and interquartile for the violin plot were labeled with the bold line colored with cyan and yellow.
Ijms 22 05686 g004
Figure 5. The circadian locomotor activity rhythm of medaka fish and AB strain zebrafish as the outgroup. (A) Comparison of the average speed between each tested fish species (medaka and zebrafish) during the day and night cycles. Comparisons of (B,H) average speed, (C,I) average angular velocity, (D,J) meandering, (E,K) freezing movement time ratio, (F,L) swimming movement time ratio, and (G,M) rapid movement time ratio between each tested fish species (medaka and zebrafish) in the day and night cycles, respectively. Data are presented as median with interquartile range, except for (G,M), which are presented as mean with a 95% confidence interval (CI) since some median of the data in those figures are 0. Data were analyzed by Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c, d, e) on the error bars represent a significant statistical difference (p < 0.05) (n = 18 for O. dancena, O. latipes, O. woworae, and D. rerio; n = 12 for O. sinensis).
Figure 5. The circadian locomotor activity rhythm of medaka fish and AB strain zebrafish as the outgroup. (A) Comparison of the average speed between each tested fish species (medaka and zebrafish) during the day and night cycles. Comparisons of (B,H) average speed, (C,I) average angular velocity, (D,J) meandering, (E,K) freezing movement time ratio, (F,L) swimming movement time ratio, and (G,M) rapid movement time ratio between each tested fish species (medaka and zebrafish) in the day and night cycles, respectively. Data are presented as median with interquartile range, except for (G,M), which are presented as mean with a 95% confidence interval (CI) since some median of the data in those figures are 0. Data were analyzed by Kruskal–Wallis test continued with uncorrected Dunn’s test. Different letters (a, b, c, d, e) on the error bars represent a significant statistical difference (p < 0.05) (n = 18 for O. dancena, O. latipes, O. woworae, and D. rerio; n = 12 for O. sinensis).
Ijms 22 05686 g005
Figure 6. (A) Principal component analysis and (B) hierarchical clustering analysis of multiple behavior endpoints in several different medaka species and zebrafish. In (A), two major clusters from hierarchical clustering analysis results were marked with the red color (1st cluster) and blue (2nd cluster) circle. The behavioral data from zebrafish were also included as the outgroup to conduct a more in-depth study about their behavior differences pattern.
Figure 6. (A) Principal component analysis and (B) hierarchical clustering analysis of multiple behavior endpoints in several different medaka species and zebrafish. In (A), two major clusters from hierarchical clustering analysis results were marked with the red color (1st cluster) and blue (2nd cluster) circle. The behavioral data from zebrafish were also included as the outgroup to conduct a more in-depth study about their behavior differences pattern.
Ijms 22 05686 g006
Figure 7. Phylogenetic relationships among four Oryzias species and Danio rerio inferred from the datasets of NADH dehydrogenase subunit 2 gene (nd2). The phylogenetic relationships were analyzed by the maximum likelihood method. The numbers beside the branches indicate the bootstrap values. Based on the prior publication by Murata et al. (2019), O. latipes and O. sinensis genetically belong to the latipes species group while O. dancena and O. woworae belong to javanicus and celebensis species groups, respectively.
Figure 7. Phylogenetic relationships among four Oryzias species and Danio rerio inferred from the datasets of NADH dehydrogenase subunit 2 gene (nd2). The phylogenetic relationships were analyzed by the maximum likelihood method. The numbers beside the branches indicate the bootstrap values. Based on the prior publication by Murata et al. (2019), O. latipes and O. sinensis genetically belong to the latipes species group while O. dancena and O. woworae belong to javanicus and celebensis species groups, respectively.
Ijms 22 05686 g007
Table 1. Comparison of neurotransmitters, antioxidant activity, and oxidative stress contents in each medaka species brain tissue measured using enzyme-linked immunosorbent assay (ELISA). The data are expressed as the mean with SD. Kruskal–Wallis test continued with uncorrected Dunn’s test was used to analyze the data. Different SSD letters (a, b, c) represent significant statistical differences (p < 0.05; n = 10 for all groups; * SSD = statistically significant difference).
Table 1. Comparison of neurotransmitters, antioxidant activity, and oxidative stress contents in each medaka species brain tissue measured using enzyme-linked immunosorbent assay (ELISA). The data are expressed as the mean with SD. Kruskal–Wallis test continued with uncorrected Dunn’s test was used to analyze the data. Different SSD letters (a, b, c) represent significant statistical differences (p < 0.05; n = 10 for all groups; * SSD = statistically significant difference).
BiomarkersO. dancenaO. latipesO. woworaeO. sinensisUnit
ConcentrationSSD * LetterConcentrationSSD LetterConcentrationSSD LetterConcentrationSSD Letter
5-HT (Serotonin)36.21 ± 11.25a25.95 ± 5.634b40.59 ± 7.444a27.16 ± 7.525bng/total protein (mg)
CAT (Catalase)15.83 ± 5.214a, c10.16 ± 2.194b14.35 ± 2.980a10.77 ± 2.829b, cng/total protein (mg)
NE (Norepinephrine)1.512 ± 0.6271a0.8307 ± 0.2756b1.653 ± 0.5085a1.011 ± 0.4499bng/total protein (mg)
DA (Dopamine)22.86 ± 9.282a12.84 ± 3.738b24.93 ± 5.888a14.45 ± 4.942bpg/total protein (mg)
Cortisol228.1 ± 69.08a135.5 ± 44.60b283.3 ± 68.17a146.9 ± 68.29bpg/total protein (mg)
AChE (Acetylcholinesterase)29.82 ± 6.940a, c17.09 ± 9.654b35.92 ± 7.873a23.44 ± 15.40b, cU/total protein (mg)
ROS (Reactive Oxygen Species)133.5 ± 44.48a, c80.81 ± 31.65b199 ± 69.21a89.78 ± 44.00b, cIU/total protein (mg)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Audira, G.; Siregar, P.; Chen, K.H.-C.; Roldan, M.J.M.; Huang, J.-C.; Lai, H.-T.; Hsiao, C.-D. Interspecies Behavioral Variability of Medaka Fish Assessed by Comparative Phenomics. Int. J. Mol. Sci. 2021, 22, 5686. https://doi.org/10.3390/ijms22115686

AMA Style

Audira G, Siregar P, Chen KH-C, Roldan MJM, Huang J-C, Lai H-T, Hsiao C-D. Interspecies Behavioral Variability of Medaka Fish Assessed by Comparative Phenomics. International Journal of Molecular Sciences. 2021; 22(11):5686. https://doi.org/10.3390/ijms22115686

Chicago/Turabian Style

Audira, Gilbert, Petrus Siregar, Kelvin H.-C. Chen, Marri Jmelou M. Roldan, Jong-Chin Huang, Hong-Thih Lai, and Chung-Der Hsiao. 2021. "Interspecies Behavioral Variability of Medaka Fish Assessed by Comparative Phenomics" International Journal of Molecular Sciences 22, no. 11: 5686. https://doi.org/10.3390/ijms22115686

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop