In this study, we acquired a head image of individual inbred and non-inbred medaka and extracted a distribution pattern of dark spots on a mesencephalic region over time. Comparing the patterns of the six inbred medaka, we found different distribution patterns over 34 weeks (Fig. 1). Although part of the pattern changed gradually in some individuals, we could attribute each pattern to each individual (Fig. 1). Then, we expanded the analysis to 30 individuals and found that all distribution patterns could be distinguished (Fig. 2) and that pattern characteristics were maintained after four weeks (Fig. 3). Moreover, in the blind identification test, all three examinees showed 100% correct answer rates (Table 1). From these findings, we concluded that individual inbred medaka could be identified over a long time period based on the characteristics of its pattern of dark spots. In addition, in contrast to work in wild animals, in laboratory animals we can determine the time interval for image acquisition at will.
The appearance of the dark spots on the mesencephalic region changed over time in the same individual (Figs. 1–3). Dark spot recognition, however, was affected by image quality (Supplementary Figs. 1, 2), so not-discernible spots in an image does not mean nonexistent dark spot. Although the image quality varied during the long observation period, it would be reasonable to assume that a dark spot discerned at two consecutive time points was real. The case of vanishing and reemerging in Fig. 1, since it only vanished at one time point, this could be an artifact due to the image quality. In Fig. 1, four different time intervals (4, 6, 8 and 12 weeks) were used. At those time intervals, the appearance change was largest at 12 weeks and almost none at 4 weeks. Although the time interval was not short enough to determine whether the dark spot itself appeared suddenly or gradually, the distribution pattern as a whole changed gradually depending on interval length; therefore, the probability of changing pattern was minimized in a time window of 4 weeks. In fact, the observed change in distribution pattern in Figs. 2–3 was small enough to connect the pattern at the second time point with that at the first time point. Therefore, we could conclude that a 4-week interval was short enough for individual identification of adult inbred medaka.
We increased the number of medaka from 6 to 30 for the variation and stability verification test (Figs. 2, 3) since according to the central limit theorem in probability theory, the mean of random samples tends to be normally distributed when the number of samples is ≥ 30 even if the random sample itself does not follow a normal distribution [25]. Therefore, we utilized 30 medaka to investigate the distribution pattern, although we did not compute the arithmetic mean and statistical significance of the derived properties. Moreover, in most tank systems as well as ours, considering the tank size and fish density, medaka would be maintained as < 30 individuals in a tank, so verifying the patterns in 30 medaka would show enough robustness in the individual identification based on the distribution pattern against the individual variation in ordinary rearing conditions [26].
In non-inbred medaka, we also found six distinguishable distribution patterns with overall constant characteristics after four weeks (Fig. 4). Differences in the size and number of dark spots between inbred and non-inbred strains could be due to body coloration. On the other hand, in the non-inbred strain, the distribution pattern seemed to change more rapidly at four weeks. At the first time point in the non-inbred strain case, medaka fish were 10 weeks of age, corresponding to a rapid body-growth period. By 15 weeks, medaka’s body size has increased substantially and since sexual maturity has been reached, the growth rate slows down [27]. During the rapid body-growth period, the distance between dark spots would increase faster and changes in distance and direction might not be homogeneous. Therefore, the larger changes observed in the non-inbred strain were probably due to the larger growth rate between time points.
In the blind identification test of non-inbred medaka, two examinees showed 100% correct answer rate but one examinee only 67% (Table 1); however, this result confirmed individual identification even in non-inbred medaka. In the test, Examinee 3 failed in the identification of Individuals #39 and #42, in which the size of dark spots changed between time points. Since the dark spot patterns maintained the same characteristics, for individual identification it would be better to pay attention to individual dark spots and the overall pattern. Nevertheless, a large number of dark spots and fast growth rate would make difficult tracing individuals by their distribution patterns. Thus, a shorter interval of image acquisition might be better. In addition, the patterns in other body parts such as the dorsal side of the body just above the spinal cord could be utilized as supplementary information (Fig. 5). Despite the difficulty of focusing on a wider region, the broader distribution pattern of the dark spots might help tracking individuals. Moreover, a guide in extracted images of the mesencephalic region might help recognizing a characteristic pattern. In the enlarged image of inbred medaka, a vertical and horizontal line guide were placed at the automatically-chosen center (Figs. 2, 3, Supplementary File 1). In the blind test of non-inbred medaka, however, the guide was removed (Fig. 4, Supplementary File 2), suggesting that this guide helped perceive the relative position of the dark spots, even if not at the same anatomical position in both two time points.
There are other inbred medaka strains aside from the HdrR strain utilized in this study. Given differences in craniofacial morphology between inbred medaka strains [28], the distribution patterns of dark spots in other inbred medaka strains can differ from the ones in this study. Although a large and small number of dark spots were investigated in this study, verification in other inbred strains is needed. On the other hand, the genetic variation within the same inbred strain can be estimated as well. In this study, the F97 of the inbred strain was utilized. At F97, the probability of heterozygote occurrence can be calculated at 1.06 × 10− 9. Considering that the genomic size of medaka is 700.4 Mb [1] and assuming that all variations come from the heterozygote occurrence, individual variation in the genome of the same inbred strain can be estimated to exist in 0.74 bases (700.4 × 106 × 1.06 × 10− 9). The estimated ~ 0.74 base number may represent a single-nucleotide polymorphism in the medaka genome. Although the existence of single-nucleotide polymorphism changes might increase susceptibility to a disease [29], future research is required to unveil the causes for the observed individual variations of the distribution pattern of dark spots in the inbred strain.
In this study, a digital compact camera with digital microscopic mode was used to capture an enlarged view of medaka’s mesencephalic region. To increase image quality, a stereomicroscope with a digital camera could be utilized, although it would involve higher cost and possible need for anesthesia. In addition, qualitative measures rely on recognition of distribution pattern characteristics. Using a pattern recognition method in machine learning [30] might allow quantifying characteristics and automate identification. It might also lead to more precise identification and studies in the original environment, as in ecological studies. Furthermore, in blind identification tests, a one-to-many identification could be performed only in the beginning. Since all images of individual medaka at both time points were available, the method of elimination could be employed. In laboratory animals this is not a problem, since similar conditions to those in this study would be used. In addition, individuals with difficult distribution patterns to distinguish could be separately maintained after the initial observation.
Despite the potential for future improvements to the current set-up, our proposed method of individual identification based on the distribution pattern of dark spots would serve as animal biometrics relevant to medical and neurological research in medaka. In a time-series analysis, the proposed method would make it possible to return medaka to its original rearing environment after one experiment without imposing any burden. Therefore, it can be used in the long term observation needed for research on neurological and neurodegenerative diseases and on drugs with a long administration period.
In conclusion, we showed that the distribution pattern of dark spots on the head of inbred and non-inbred medaka differed in each individual and maintained its characteristics over time. In blind identification tests, three examinees scored high correct answer rates; thus, in the typical rearing environment, this pattern, recorded at 4-week intervals for adult and shorter time intervals for young medaka, could be utilized as biometric identifier, and play an important role in long term observations in medical and neurological studies.