International Journal of Environmental Research and Public Health

Article Analyzing the Response Behavior of variegatus (: ) to Different Concentrations of Copper Sulfate Based on Line Body Shape Detection and a Recurrent Self-Organizing Map

Chang Woo Ji 1, Young-Seuk Park 2,* , Yongde Cui 3, Hongzhu Wang 3, Ihn-Sil Kwak 1 and Tae-Soo Chon 4,*

1 Fisheries Science Institute, Chonnam National University, Yeosu 59626, Korea; [email protected] (C.W.J.); [email protected] (I.-S.K.) 2 Department of Biology and Department of Life and Nanopharmaceutical Sciences, Kyung Hee University, Seoul 02447, Korea 3 Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China; [email protected] (Y.C.); [email protected] (H.W.) 4 Ecology and Future Research Association (EnFRA), Dusil-ro 45 beon-gil 21, Geumjeong-gu, Busan 46228, Korea * Correspondence: [email protected] (Y.-S.P.); [email protected] (T.-S.C.); Tel.: +82-2-961-0946 (Y.-S.P.); +82-51-512-2262 (T.-S.C.)  Received: 30 March 2020; Accepted: 7 April 2020; Published: 11 April 2020 

Abstract: Point detection (e.g., the centroid of the body) of species has been conducted in numerous studies. However, line detection (i.e., the line body shape) of elongated species has rarely been investigated under stressful conditions. We analyzed the line movements of an Oligochaeta in response to treatments with a toxic chemical, copper sulfate, at low concentrations (0.01 mg/L and 0.1 mg/L). The automatic line-tracking system was devised to identify the movement of body segments (body length) and the movements of segments (i.e., the speed and angles between segments) were recorded before and after treatment. Total body length was shortened from 31.22 ( 5.18) mm to 20.91 ( 4.65) mm after the 0.1 mg/L treatment. The Shannon entropy index decreased ± ± from 0.44 ( 0.1) to 0.28 ( 0.08) after treatment. On the other hand, the body and movement segments ± ± did not significantly change after the 0.01 mg/L treatment. Sequential movements of test organisms were further analyzed with a recurrent self-organizing map (RSOM) to determine the pattern of time-series line movements. The RSOM made it feasible to classify sequential behaviors of indicator organisms and identify various continuous body movements under stressful conditions.

Keywords: Lumbriculus variegatus; line-tracking; response behavior; computational method; copper sulfate; risk assessment; Shannon entropy

1. Introduction Response behavior of indicator species is considered an efficient means of monitoring risk in aquatic ecosystems [1–4]. Real-time behavior measurements are able to fill the gap between macroscale (e.g., community surveys) and microscale (e.g., individual development) assessments [5,6]. Response behaviors of indicator species are sensitive to sublethal exposure of various toxic chemicals at low concentrations [1,7]. Response behaviors have been monitored in various taxa, including insects [8,9], crustaceans [2,10], snails [11], and fish [12–14].

Int. J. Environ. Res. Public Health 2020, 17, 2627; doi:10.3390/ijerph17082627 www.mdpi.com/journal/ijerph Int. J. Environ. Res. Public Health 2020, 17, 2627 2 of 15

Considering that indicator species are primarily either round (e.g., Daphnia), elongated, or fusiform (e.g., fish), point estimation methods (e.g., centroid of the body location) are conveniently used for behavioral monitoring [13–15]. However, it is also feasible to monitor species with long body shapes under stressful conditions, including oligochaetes [16] and nematodes (e.g., Caenorhabditis elegans)[17]. Compared with point movement, line movement provides an extra dimension to express the status of an indicator species. Line detection is primarily used with the nematode Caenorhabditis elegans [18–20]. C. elegans were also used to describe line movements of organisms by tracking the centroids of their x and y positions [21–23] and calculating the angle change rate of the body skeleton [19]. Line movements of C. elegans were used to differentiate between response behaviors of different mutants after exposure to benzene [24]. Computational methods, such as the Frenet–Serret formula, were also used to describe sequential line movements of C. elegans [20,25]. Recently, chemical engineering methods were developed to detect and remove copper ions from aqueous media [26–28]. These methods were relatively more convenient than behavioral monitoring methods. However, the potential environmental toxicity is not well characterized [29] and the risks they may pose to aquatic organisms and the environment are not certain. Oligochaeta, including Lumbriculus variegatus, are important fauna in the aquatic food chain or web and are widely used as bioindicators to pollution. In particular, L. variegatus was identified as an aquatic organism to assess test specimens of bioaccumulation by the US Environmental Protection Agency [30]. Although oligochaetes were used to indicate toxic effects [16,31,32], their line body shapes were not extensively studied or compared with those of C. elegans. Both nematodes and oligochaetes possess line body shapes, but their overall shapes are different. The body is more elongated in oligochaetes and the scope of movements is substantially more complex due to this body elongation. Although line movement shapes were not computationally analyzed for oligochaetes, various movement patterns were reported in response to stimuli. The body reverse and swimming movements of L. variegatus, for instance, were observed when the worms were touched by anterior segments [33]. The activities of Oligochaeta were recorded using a quadrupole impedance conversion technique under chemical stress [2,16,34], which showed an increase in locomotory activity of L. variegatus in an intermediate concentration of lead (1–5 mg/L) [16]. However, changes in body shapes of the oligochaetes per se were not extensively studied, especially in response to chemical stressors. Regarding line body shape methodology, various computational approaches were introduced to analyze complex movement behaviors [13,15,35]. A multi-layer perceptron, which is a supervised artificial neural network, was used to detect changes in movement patterns of medaka fish (Oryzias latipes) before and after treatment with diazinon [15]. However, the multi-layer perceptron requires templates a priori to recognize the movement patterns in a supervised manner. Meanwhile, a self-organizing map (SOM), which is an unsupervised artificial neural network, was proposed to characterize movement patterns in an unsupervised learning manner (i.e., without a priori knowledge) and was implemented to classify response behaviors of various indicator species, including cockroaches [9] and medaka fish [13,14]. Son et al. [36] used a recurrent self-organizing map (RSOM), an extended version of a SOM, to pattern sequential data, in this case, time-series movements of L. variegatus to detect the effects of treatments with copper. Because of the complexity of behavioral data, one chemical concentration was primarily studied when analyzing movement patterns [13,15,37]. The temporal network models were more feasible in learning time-series data than conventional methods based on linear (e.g., Autoregressive Model and Autoregressive Moving Average-Neural Network) and non-linear (e.g., Nonlinear Autoregressive Moving Average with the Exogenous Input and Multivariate Adaptive Regression Splines) statistical analyses [38]. In ecology, time-series data for community dynamics were predicted using the Elman network [39] and a recurrent network [40]. These networks are necessary for training templates (i.e., supervised learning). However, obtaining all the pre-determined patterns for behavioral data is difficult due to the numerous complex patterns. Subsequently, a huge amount of behavioral data is required to be prepared through continuous recording. Data mining provides the overall scope of behavioral data; in this regard, we utilized the Int. J. Environ. Res. Public Health 2020, 17, 2627 3 of 15 unsupervised network to accommodate the sequential line movement data by collecting data from continuous recording. The Temporal Kohonen Map (TKM) [38], which was derived from the Kohonen Self-Organizing Map [41], is considered as an efficient learning tool for temporal sequence processing. The involvement of the earlier input vectors in each unit occurred following the recursive difference in TKM. RSOM and an unsupervised temporal model provided more flexibility when dealing with sequential data. In RSOM, direct learning of the temporal context is possible, while TKM is unable to directly use the temporal information of input sequences in weight updating [42]. TKM allows the use of many data with only a little a priori knowledge, while RSOM provides constant results to deal with the classification of temporal data with little effort [42–44]. In this study, we focused on evaluating the effects of exposure to chemicals on the body segments of L. variegatus at different concentrations using computational methods, including recurrent SOM. The body segments (i.e., body length), which were only obtained from line body shape detection, were characterized according to different concentrations of copper sulfate. We demonstrated that the unsupervised temporal model feasibly patterned the sequential line-movement of oligochaetes after toxic substance treatment and efficiently characterized the stressful behavior of the specimens.

2. Materials and Methods

2.1. Test Organisms and Observation System

L. variegatus were collected at the long-term field survey station, Beopgi (35◦21044” N, 129◦06015” E, Yangsan, Republic of Korea), in August 2006, and used as a stock population in the Laboratory of Ecology and Behavior System, Department of Biological Sciences, Pusan National University, Republic of Korea. The sample site was located upstream in the Suyong River and was free from anthropogenic effects due to the protection of the area by the local government for drinking water. The stock population was reared in a glass tank (5 L) with an artificial dry diet (Tetramin®) at 20 1 C. Dry ± ◦ leaves (oak, Quercus sp.) gathered at the sample site were used as the substrate in the tank. Only mature organisms that were 4 cm ( 0.5 cm) in length were used for behavioral observations. ± The movement behaviors of test specimens were individually observed and recorded using the observation system for 6 h without treatment. Subsequently, copper sulfate (CuSO4) was used as a treatment for the same organisms, which were recorded for an additional 6 h after treatment. After being dissolved in distilled water, the chemical was applied to the observation aquarium at two concentrations (0.01 and 0.1 mg/L). The LC50 of copper sulfate was reported as 0.23 mg/L[45]. The response behaviors of 20 individuals were observed before and after chemical treatments at each concentration using the observation system. After consideration of noise and observed data quality due to the elongated body shape in recorded images, 12 of 20 observed individuals at each concentration were selected for further analyses.

2.2. Detection of Line Movement L. variegatus were individually placed in an observation aquarium (9 cm diameter) and their positions were scanned from the top view at 0.25 s intervals using a CCTV (Closed circuit television) camera (Kukjae Electronics Co. Ltd.; IVC-841). The time interval of 0.25 s was sufficiently short to observe the slow, overall responses of the line bodies of the oligochaetes. The backlight conditions (backlight; 0.2 Watt green diodes) were provided underneath the observation cage according to Son et al. [36] to enhance visual recognition of the organisms. Before placing the test organisms in the observation aquarium, a background image was acquired with an image size of 320 240 pixels. × Then, a new image with the test organism in the observation aquarium was compared with the background image. Considering the complexity in calculating line movements of with elongated body shapes [18,19], we adopted four processes to record the positions of test organisms before analyzing Int. J. Environ. Res. Public Health 2020, 17, 2627 4 of 15 movements. First, we obtained line skeletons of the elongated bodies of the oligochaetes. Second, the heads and tails of the line skeletons were determined by experienced observers after watching continuous movement sequences. Due to the difficulty in automatic recognition of the head and tail in the elongated bodies of the oligochaetes, direct observations were performed to head and tail positions. For instance, the head initiated the movement and the rest of the body followed the head movement in a complex manner, with the last movement usually ending in a tail movement. Additionally, a wide range of turning was observed in the head area in a short period (broad angular speed). However, direct observation was time-consuming and some segments were difficult to identify as either the head or tail, particularly when test specimens performed reversal swimming, where the position of the head and the tail on the elongated body rapidly switched several times. When the body was released again, it was difficult to identify either the head or the tail. The segments with unclear identification were excluded from analysis, accounting for approximately 40% of observed data. Approximately 60% of the segments with activity were chosen for analysis. Third, the body skeleton was divided evenly into 12 segments. Based on the position of each segment, length (12 data points) and angles between segments (11 data points) were obtained at every time interval (0.25 s) for an average per 1 min (23 data points for one movement segment). Fourth, changes in body segments were investigated by calculating movement parameters under stressful conditions. The length of each segment and the curvatures between neighboring segments were recorded at each time interval by image framing according to Son et al. [36].

2.3. Analysis of Movement Data Based on preliminary studies with test organisms [5,36], the following five parameters (average in 1 min) were selected to characterize the movement of the body segments of oligochaetes:

Speed (mm/s): Average of movement distance of the test organisms in 1 min; • Stop duration (s): Average time during which the test specimen did not move; • Stop number: Stop frequency during 1 min; • Turning rate (rad/s): The sum of angular changes in radians in absolute values divided by the • cumulated time duration of movement; Meander (rad/s): The sum of angular changes in radians in absolute values divided by the • cumulated movement distance.

In addition, the Shannon entropy index (H)[46] was calculated to represent diversity in movement activities in response to chemical stress, as follows:

Xn H = p(xi) log p(xi), − i=1 2 where p(xi) is the probability of observing the test organisms at each coordinate position (i). The entropy index represents the complexity of the movement activities of the test organisms. If the organisms were active with a higher degree of movement, the index was higher. An analysis of variance (ANOVA) and Tukey’s honestly significant difference test were conducted to compare differences in variables among clusters. Statistical differences in the variables before and after chemical treatment were evaluated based on the paired sample t-tests. These analyses were conducted with the software Statistica [47]. Recurrence of time series data for movement was applied to the RSOM [42]. For time-series data, the sequential line data consisting of 12 segments with 23 data points (including 12 body segment lengths and 11 angles between segments) were used as input data for the RSOM. In total, data for 5 min (25 s 12 = 300 s) were used to characterize the behavioral responses using the RSOM. Based on × a previous study [36], we used 48 (=8 6) of the RSOM output units to show the most biologically × relevant results. Int. J. Environ. Res. Public Health 2020, 17, 2627 5 of 15

The SOM is a vector quantization method to map patterns from an input space Vi onto a lower dimensional space VM of the map to determine the best matching unit b in time step t, according to the following equation:  x(t) w (t) = min x(t) w (t) , k − b k i k − i k where i V , x(t) is an input vector and w (t) is a weight vector of the units i on the map. Then, the ∈ M i weight vector of the best matching unit b is updated toward the given input vector x(t) in the SOM, according to the equation:

w (t + 1) = w (t) + γ(t)h (t)(x(t) w (t)), (1) b b b − b where γ(t), 0<γ(t) 1 is a learning rate and h (t) is the neighborhood function. ≤ b RSOM is similar to SOM, except it searches for the best matching unit as following [36,42]:

y (t) = (1 α)y (t 1) + α(x(t) w (t)), (2) i − i − − i where 0 < α 1 is a leaking coefficient, y (t) is a leaked difference vector, w (t) is the reference or weight ≤ i i vector in unit i, and x(t) is the input pattern for time step t. The best matching unit b at time step t is searched for by the equation:  y = min y (t) , (3) b i k i k where i V . The RSOM had the same weight updating process as the SOM, although the input ∈ M sequence was noted before the learning process [36]. The learning process of the RSOM was conducted with a program developed by the authors in the Matlab environment [48]. After the learning process, a hierarchical cluster analysis was conducted based on the Ward linkage method with Euclidean distance measures to define the cluster boundaries in RSOM units. A multi-response permutation procedure [49], which is a nonparametric procedure to test the hypothesis of no difference between two or more groups, was conducted to evaluate the significance of the clusters defined in the RSOM using PC-ORD [50].

3. Results

3.1. Effects on Body Length and Parameters The differences in the average body length of the organisms were evaluated before and after treatment with copper sulfate. Total length was not different before and after treatments at the 0.01 mg/L concentration (paired t-test, t = 1.62, P = 0.13), whereas it was significantly shortened from 31.22 5.18 (S.D.) mm to 20.91 4.65 (S.D.) mm after the 0.1 mg/L treatment (paired t-test, t = 4.352, ± ± P = 0.001) (Figure1a). The Shannon entropy index was not differentiated when the test organisms were treated with the 0.01 mg/L concentration; the values were 0.508 0.11 (SD) before treatment and 0.501 0.09 (SD) ± ± after treatment (paired t-test, t = 0.14, P = 0.89) (Figure1b). However, the Shannon entropy index significantly decreased after the 0.1 mg/L treatment (0.44 0.1 (SD) before treatment and 0.287 0.08 ± ± (SD) after treatment), although the difference was not statistically significant (paired t-test, t = 3.674, P = 0.036) (Figure1b). Considering the overall decrease in the entropy value, movement activities of test individuals became less diverse after treatment with 0.1 mg/L of copper. Subsequently, five parameters were evaluated according to the segments of line-shaped bodies before and after treatment at the two concentrations (Figure2). Overall, speed was lower in the middle part of the body (top panel, Figure2), whereas the meander was higher in the middle part of the body (bottom panel, Figure2) before and after treatments, regardless of concentration. Meander tended to be lower at both ends of the body, although this difference was not statistically significant (Figure2). Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 5 of 15

x(t) − wb (t) = mini{ x(t) − wi (t) },

where i ∈ VM, x(t) is an input vector and wi(t) is a weight vector of the units i on the map. Then, the weight vector of the best matching unit b is updated toward the given input vector x(t) in the SOM, according to the equation:

wb (t +1) = wb (t) +  (t)hb (t)(x(t) − wb (t)), (1)

where  (t), 0< (t) ≤ 1 is a learning rate and hb(t) is the neighborhood function. RSOM is similar to SOM, except it searches for the best matching unit as following [36,42]: y (t) = (1−)y (t −1) +(x(t) − w (t)) i i i , (2)

where 0 < α ≤ 1 is a leaking coefficient, yi(t) is a leaked difference vector, wi(t) is the reference or weight vector in unit i, and x(t) is the input pattern for time step t. The best matching unit b at time step t is searched for by the equation: y = min { y (t) } b i i , (3)

where i ∈ VM. The RSOM had the same weight updating process as the SOM, although the input sequence was noted before the learning process [36]. The learning process of the RSOM was conducted with a program developed by the authors in the Matlab environment [48]. After the learning process, a hierarchical cluster analysis was conducted based on the Ward linkage method with Euclidean distance measures to define the cluster boundaries in RSOM units. A multi-response permutation procedure [49], which is a nonparametric procedure to test the hypothesis of no difference between two or more groups, was conducted to evaluate the significance of the clusters defined in the RSOM using PC-ORD [50].

3. Results

3.1. Effects on Body Length and Parameters The differences in the average body length of the organisms were evaluated before and after treatment with copper sulfate. Total length was not different before and after treatments at the 0.01 mg/L concentration (paired t-test, t = 1.62, P = 0.13), whereas it was significantly shortened from 31.22 Int. J. Environ.± 5.18 Res. (S.D.) Public mm Health to 20.912020 ,±17 4.65, 2627 (S.D.) mm after the 0.1 mg/L treatment (paired t-test, t = 4.352, P = 6 of 15 0.001) (Figure 1a).

Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 6 of 15

The Shannon entropy index was not differentiated when the test organisms were treated with the 0.01 mg/L concentration; the values were 0.508 ± 0.11 (SD) before treatment and 0.501 ± 0.09 (SD) after treatment (paired t-test, t = 0.14, P = 0.89) (Figure 1b). However, the Shannon entropy index significantly decreased after the 0.1 mg/L treatment (0.44 ± 0.1 (SD) before treatment and 0.287 ± 0.08 (SD) after treatment), although the difference was not statistically significant (paired t-test, t = 3.674, P = 0.036) (Figure 1b). Considering the overall decrease in the entropy value, movement activities of test individuals became less diverse after treatment with 0.1 mg/L of copper. Subsequently, five parameters were evaluated according to the segments of line-shaped bodies Figure 1.FigureComparison 1. Comparison of total of total body body lengths lengths ofof LumbriculusLumbriculus variegatus variegatus (a) and(a) andShannon Shannon information information before and after treatment at the two concentrations (Figure 2). Overall, speed was lower in the entropy representing movement of test organisms (b) before and after treatments with copper sulfate entropymiddle representing part of the body movement (top panel, of test Figure organisms 2), whereas (b) beforethe meander and after was treatmentshigher in the with middle copper part sulfate of at concentrations of 0.01 mg/L and 0.1 mg/L. Paired sample t-tests were conducted to compare the atthe concentrations body (bottom of panel, 0.01 mgFigure/L and 2) before 0.1 mg and/L. after Paired treatments, sample tregardless-tests were of conductedconcentration. to compareMeander the differences in body length and Shannon entropy index before and after treatments. Error bars indicate differences in body length and Shannon entropy index before and after treatments. Error bars indicate tendedstandar to bed lower deviation. at both ends of the body, although this difference was not statistically significant standard(Figure deviation.2).

FigureFigu 2. Direff 2.erences Differences in parameters in parameters characterizing characterizing the the movement movement of body of body segments segments of Lumbriculus of Lumbriculus variegatus before and after copper sulfate treatments at two concentrations, 0.01 mg/L (a) and 0.1 mg/L variegatus before and after copper sulfate treatments at two concentrations, 0.01 mg/L(a) and 0.1 mg/L (b). Statistical differences (paired sample t-test, t0.05,(2),11 = 2.20) in the selected variables are represented (b). Statistical differences (paired sample t-test, t = 2.20) in the selected variables are represented by an asterisk (*: 0.01 < P < 0.05). Error bars indic0.05,(2),11ate standard deviation. by an asterisk (*: 0.01 < P < 0.05). Error bars indicate standard deviation.

Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 7 of 15 Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 7 of 15

TheThe chem chemicalical effects effects were were observed observed more more clearly clearly with with 0.1 0.1 mg/L mg/L compared compared with with 0.01 0.01 mg/L mg/L. .Stop Stop durationduration and and turning turning rate rate substantially substantially increased increased in in all all or or most most segments segments at at 0.1 0.1 mg/Lmg/L, ,respectively respectively (Figure(Figure 2). 2). Additionally, Additionally, a a significant significant increase increase in in stop stop number number was was observed observed in in the the posterior posterior segments.Int.segments. J. Environ. It ItRes. was was Public noteworthy noteworthy Health 2020 ,that 1that7, x only FORonly PEER the the lastREVIEWlast segment segment showed showed a a significant significant decrease decrease in in speed speed7 after ofafter 15 treatmenttreatment atat 0.10.1 mg/Lmg/L. . AtAt thethe lowlowerer concentrationconcentration of of 0.010.01 mg/Lmg/L, , the the differencedifference waswas not not as as clearlyclearly The chemical effects were observed more clearly with 0.1 mg/L compared with 0.01 mg/L. Stop observableobservable asas inin the the 0.10.1 mg/Lmg/L concentration.concentration. However,However, aa partialpartial difference difference waswas observedobserved inin speedspeed duration and turning rate substantially increased in all or most segments at 0.1 mg/L, respectively accordingaccording to to bod bodyy segments segments, ,with with a a decrease decrease in in both both the the anterior anterior (p1 (p1 and and p2) p2) and and posterior posterior segments segments (Figure 2). Additionally, a significant increase in stop number was observed in the posterior (p11(p11––p13)p13) (paired (paired t -ttest,-test, P P < < 0.05; 0.05; top top panel, panel, Figure Figure 2a), 2a), indicating indicating that that a a body body response response was was observed observed segments. It was noteworthy that only the last segment showed a significant decrease in speed after atat the the lowerlower concentration. concentration. AccordingAccording toto ANOVA ANOVA andand Tukey Tukey’s’s hhonestlyonestly s significantignificant d differenceifference test test treatment at 0.1 mg/L. At the lower concentration of 0.01 mg/L, the difference was not as clearly (Table(Table 1), 1), speed speed and and turning turning rate rate were were divided divided into into two two groups, groups, whereas whereas stop stop duration duration and and meander meander observable as in the 0.1 mg/L concentration. However, a partial difference was observed in speed werewereInt. J. not Environ.not different different Res. Public between between Health 2020the the, experimental 17experimental, 2627 con conditions.ditions. 7 of 15 according to body segments, with a decrease in both the anterior (p1 and p2) and posterior segments TheThe overall overall speed speed was was lower lower in in the the middle middle parts parts of of the the body. body. Speed Speed and and mea meandernder showed showed high high (p11–p13) (paired t-test, P < 0.05; top panel, Figure 2a), indicating that a body response was observed correlationscorrelations at at both both concentrations concentrations (r (r = = −−0.8150.815 and and − −0.915,0.915, respectively, respectively, for for 0.01 0.01 mg/L mg/L and and 0.1 0.1 mg/L mg/L at the Thelower chemical concentration. effects were According observed to moreANOVA clearly and with Tukey 0.1’s mg honestly/L compared significant with 0.01difference mg/L. test Stop concentrations;concentrations; P P < < 0.0001 0.0001 for for both both concentrations), concentrations), indicating indicating that that the the linear linear and and circular circular movements movements (Tableduration 1), speed and turning and turning rate substantiallyrate were divided increased into two in all groups, or most whereas segments stop at duration 0.1 mg/ L,and respectively meander werewere sensitive sensitive to to the the chemical chemical but but negatively negatively associated associated with with each each other. other. Based Based on on the the results results of of body body were(Figure not2 ).different Additionally, between a significant the experimental increase con in stopditions. number was observed in the posterior segments. sizesize (Figure (Figure 1) 1) and and parameters parameters obtained obtained from from body body segments segments (Figure (Figure 2, 2, Table Table 1), 1), the the overall overall effects effects of of It wasThe noteworthy overall speed that was only lower the lastin the segment middle showed parts of a the significant body. Speed decrease and inmea speednder after showed treatment high treatmentstreatments appeared appeared to to be be distinctive distinctive at at the the 0.1 0.1 mg/L mg/L concentration concentration and and partially partially det detectableectable in in speed speed correlationsat 0.1 mg/L. at At both the concentrations lower concentration (r = −0.815 of 0.01 and mg −0.915,/L, the respectively, difference was for not0.01 asmg/L clearly and observable 0.1 mg/L withinwithin the the body body shape shape at at the the lower lower concentration concentration of of 0.01 0.01 mg/L mg/L. . concentrations;as in the 0.1 mg P/ L< 0.0001 concentration. for both concentrations), However, a partial indicating difference that was the observedlinear and in circular speed movements according to bodyTable segments, 1. Tukey’s with multiple a decrease comparison in both test the for anteriorparameters (p1 between and p2) two and concentrations posterior segments (a: before (p11–p13)0.01 were sensitiveTable 1. Tukey’s to the chemicalmultiple comparison but negatively test for associated parameters with between each other.two concentrations Based on the (a: results before 0.01of body (pairedmg/Lt-test,; b: afterP < 0.010.05; mg/L top; c: panel, before Figure 0.1 mg/L2a),; d: indicating after 0.1 mg/L that). a Groups body response are divided was by observed upper lines at theand lower size (Figuremg/L; b: 1) after and 0.01 parameters mg/L; c: before obtained 0.1 mg/L from; d: body after segments 0.1 mg/L). (FigureGroups are2, Table divided 1), by the upper overall lines effects and of concentration.under lines based According on Tukey’s to ANOVA multiple and comparison Tukey’s test honestly (P < 0.05). significant difference test (Table1), speed treatmentsunder appearedlines based to on be Tukey’s distinctive multiple at thecomparison 0.1 mg/L test concentration (P < 0.05). and partially detectable in speed withinand turning the body rate shape were dividedat the lower into concentration two groups, whereas of 0.01 stopmg/L duration. and meander were not different between the experimental conditions. PositionPosition SpeedSpeed StopStop Duration Duration StopStop Number Number TurningTurning Rate Rate MeanderMeander Table 1. Tukey’s multiple comparison test for parameters between two concentrations (a: before 0.01 p1 n.s*. n.s. n.s. c b a d n.s. mg/LTablep1; b: 1. afterTukey’s 0.01n.s multiplemg/L*. ; c: before comparison 0.1n.s. mg/L test; d: for after parameters 0.1 mg/Ln.s. between). Groups two are concentrations dividedc b a d by upper (a: beforelinesn.s. and 0.01 mg/L; b: after 0.01 mg/L; c: before 0.1 mg/L; d: after 0.1 mg/L). Groups are divided by upper lines and underp2 lines baseda b onc d Tukey’s multiplen.s. comparison test (Pn.s. < 0.05). c b a d n.s. underp2 lines baseda b c on d Tukey’s multiplen.s. comparison test (Pn.s.< 0.05). c b a d n.s. Positionp3 aSpeed b c d Stop Durationn.s. Stop n.s.Number Turningc b a dRate Meandern.s. p3Position a b c Speedd Stopn.s. Duration Stopn.s. Number Turningc b a Rated Meandern.s.

p4p1p4 p1aa n.sb b c *c .d d n.s *.n.s.n.s. n.s. n.s. n.s. cc b b a a d d n.s.n.s.n.s. p2 a b c d n.s. n.s. n.s. p5p2p5 aa b b c c d d n.s.n.s. n.s.n.s. cc b b a a d d n.s.n.s. p3 a b c d n.s. n.s. n.s. p6 a b c d n.s. n.s. c b a d n.s. p3p6 p4 aa b b c c da d b c d n.s.n.s.n.s. n.s. n.s. cc b b a a d d n.s.n.s.n.s. p5 a b c d n.s. n.s. n.s. p7p4p7 aa b b c c d d n.s.n.s. n.s.n.s. cc b b a a d d n.s.n.s. p6 a b c d n.s. n.s. n.s. p8 a b c d n.s. n.s. c b a d n.s. p5p8 p7 aa b b c c da d b c d n.s.n.s.n.s. n.s. n.s. cc b b a a d d n.s.n.s.n.s. p8 a b c d n.s. n.s. n.s. p9p6p9 aa b b c c d d n.s.n.s. n.s.n.s. cc b b a a d d n.s.n.s. p9 a b c d n.s. n.s. n.s. p10 a b c d n.s. a b c d c b a d n.s. p7p10 p10 aa b b c c da d b c d n.s.n.s.n.s. an.s. b c d cc b b a a d d n.s.n.s.n.s. p11p8p11 p11 aa b b c c da d b c d n.s.n.s.n.s. aa n.s.b b c c d d cc b b a a d d n.s.n.s.n.s. p12 a c b d n.s. n.s. p12p9p12 aa cb c b cb d d n.s.n.s. aa n.s.b b c c d d cc b b a a d d n.s.n.s. p13 a c b d n.s. n.s. p13p10 a cb bc d n.s. a b c d c b a d n.s. p13 a c b d n.s.* n.s.: statistically nota significant b c d c b a d n.s. ** n.s.: n.s.:p11 statistically statistically a b notc not d significant significant n.s. a b c d c b a d n.s. The overall speed was lower in the middle parts of the body. Speed and meander showed high 3.2. Effectsp12 on Anglesa c of b Bodyd Segmentsn.s. a b c d c b a d n.s. 3.2.correlations Effects on atAngles both of concentrations Body Segments (r = 0.815 and 0.915, respectively, for 0.01 mg/L and 0.1 mg/L − − concentrations; P < 0.0001 for both concentrations), indicating that the linear and circular movements TheThep13 curvature curvature a c of ofb d the the body body waswasn.s. representedrepresented by bya ab a seriesc series d of of angles anglesc b a observed observedd between betweenn.s. every every consecutivewere sensitive pair to of the body chemical segments, but negativelycreating an associated extra dime withnsion each of movement other. Based patterns. on the results The patterns of body consecutive* n.s.: statistically pair of body not significant segments, creating an extra dimension of movement patterns. The patterns size (Figure1) and parameters obtained from body segments (Figure2, Table1), the overall e ffects of ofof angle angle changes changes were were different different in in the the different different concentrations concentrations (Figure (Figure 3). 3). The The curvature curvature of of the the head head (a1)treatments and tail appeared(a11) varied to beaccording distinctive to treatment at the 0.1 strength mg/L concentration of copper, as and shown partially for both detectable the 0.01 in mg/L speed 3.2.(a1) Effects and tail on (a11)Angles varied of Body according Segments to treatment strength of copper, as shown for both the 0.01 mg/L within the body shape at the lower concentration of 0.01 mg/L. (Figure(Figure 3a) 3a) and and 0.1 0.1 mg/L mg/L (Figure (Figure 3b) 3b) concentrations concentrations. .Curvature Curvature strongly strongly differed differed in in the the middle middle of of the the The curvature of the body was represented by a series of angles observed between every consecutive3.2. Effects onpair Angles of body of Body segments, Segments creating an extra dimension of movement patterns. The patterns of angleThe changes curvature were of different the body in the was different represented concentrations by a series (Figure of angles 3). The observed curvature between of the head every (a1)consecutive and tail (a11) pair ofvaried body according segments, to creating treatment an strength extra dimension of copper, of movementas shown for patterns. both the The 0.01 patterns mg/L (Figureof angle 3a) changes and 0.1were mg/L di (Figurefferent 3b) in the concentrations different concentrations. Curvature strongly (Figure3 differed). The curvature in the middle of the of head the (a1) and tail (a11) varied according to treatment strength of copper, as shown for both the 0.01 mg/L (Figure3a) and 0.1 mg /L (Figure3b) concentrations. Curvature strongly di ffered in the middle of the body at 0.1 mg/L. The curvature of the tail (a11) was sensitive to increases after the 0.1 mg/L treatment (0.389 0.033 and 0.494 0.088 for before and after treatments, respectively) (paired t-test, t = 2.247, ± ± P = 0.046). Curvatures were not different before or after the 0.01 mg/L treatment in any segment, Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 8 of 15 bodyInt. J. Environ.at 0.1 mg/L Res. Public. The Health curvature2020, 17 of, 2627 the tail (a11) was sensitive to increases after the 0.1 mg/L treatment 8 of 15 (0.389 ± 0.033 and 0.494 ± 0.088 for before and after treatments, respectively) (paired t-test, t = 2.247, P = 0.046). Curvatures were not different before or after the 0.01 mg/L treatment in any segment, includingincluding the the head head and and tail tail segments (paired t--test,test, t = 0.570.5733 and 0.503, P = 0.5770.577 and and 0.624 0.624 for for head head andand tail, tail, respectively). respectively). Curvature Curvature specifically specifically at at the the head head segment segment (a1) (a1) was was not not significantly significantly different different between the before (0.393 0.046) or after 0.1 mg/L treatment groups (0.397 0.055) (paired t-test, between the before (0.393 ±± 0.046) or after 0.1 mg/L treatment groups (0.397 ± 0.055)± (paired t-test, t = 0.091,t = 0.091, P = P0.929).= 0.929).

FigureFigure 3. 3. DifferencesDifferences in in curvatures curvatures of of the the head head (a1) (a1) and and tail tail (a11) (a11) at attwo two different different concentrations concentrations of copperof copper sulfate sulfate:: (a) ( 0.01a) 0.01 mg/L mg /andL and (b) ( b 0.1) 0.1 mg/L mg/.L. Head Head and and tail tail regions regions were were identified identified by by direct observation.observation. Statistical Statistical difference differencess (paired (paired sampled t-test,-test, tt0.05,(2),110.05,(2),11 = = 2.20)2.20) in in the the curvatures curvatures are representedrepresented by an asterisk (*: P < 0.05;0.05; **: **: PP << 0.01).0.01). Error Error bars bars indicate indicate standard standard deviation. deviation. 3.3. Patterning Sequential Line Movement 3.3. Patterning Sequential Line Movement Because the line movement was recorded in time series, continuous movements were patterned Because the line movement was recorded in time series, continuous movements were patterned before and after chemical treatment via RSOM. Considering the effects of the chemical were conspicuous before and after chemical treatment via RSOM. Considering the effects of the chemical were at high concentrations (Figures1 and2), we used data for the 0.1 mg /L copper sulfate treatment to conspicuous at high concentrations (Figures 1 and 2), we used data for the 0.1 mg/L copper sulfate pattern the sequential line movement. The movement patterns were accordingly grouped using RSOM treatment to pattern the sequential line movement. The movement patterns were accordingly (Figure4a). Di fferent sizes of circles indicated the number of patterns assigned in each SOM unit. grouped using RSOM (Figure 4a). Different sizes of circles indicated the number of patterns assigned A gradient was observed in the vertical position with a slight inclination to the diagonal. The movement in each SOM unit. A gradient was observed in the vertical position with a slight inclination to the patterns were grouped vertically into three clusters (I, II, and III) (Figure4b). The before-treatment diagonal. The movement patterns were grouped vertically into three clusters (I, II, and III) (Figure groups (white circles) were mainly placed in the upper right part of the map, whereas the after-treatment 4b). The before-treatment groups (white circles) were mainly placed in the upper right part of the groups (black circles) occupied a broader area in the bottom left part of the map. The movement map, whereas the after-treatment groups (black circles) occupied a broader area in the bottom left patterns before treatment were mainly grouped in cluster I, whereas the movement patterns after part of the map. The movement patterns before treatment were mainly grouped in cluster I, whereas treatment were grouped in cluster III. the movement patterns after treatment were grouped in cluster III. The body movements were correspondingly different in different clusters (Figure 5). Cluster I showed a sequence of large and straightened shapes, cluster II displayed folding shapes, and cluster III represented contracted and folding shapes. As stated above, body shapes in cluster II were mostly movements before the chemical treatment, whereas those in cluster III were mostly movements observed after the treatment. Therefore, body shapes after the chemical treatments were represented by short body lengths with a high degree of folding.

Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 9 of 15

Each movement parameter differed according to different clusters in the RSOM (Figure 6). To compare parameters, 1000 samples were randomly selected from each cluster. The total length in cluster I showed the greatest range among the clusters (Figure 6a), with length also being relatively greater in cluster II, another cluster belonging to the before-treatment group, compared with cluster III, which represented after treatment. The total lengths of each cluster significantly differed from Int. J. Environ. Res. Public Health 2020, 17, x FOR PEER REVIEW 9 of 15 each other (Tukey’s test, q0.05,3000,4 = 3.314; I ≠ II ≠ III) (Table 2). The angles between body segments were also significantly different among clusters. In cluster I, the angles were relatively greater at the Each movement parameter differed according to different clusters in the RSOM (Figure 6). To end of the body, however, curvatures were greater in the middle part of the body. Patterns in the compare parameters, 1000 samples were randomly selected from each cluster. The total length in RSOM were capable of characterizing the sequential line shapes of test organisms before and after clusterInt. J. Environ. I showed Res. Public the Healthgreatest2020 range, 17, 2627 among the clusters (Figure 6a), with length also being relatively 9 of 15 treatment. greater in cluster II, another cluster belonging to the before-treatment group, compared with cluster III, which represented after treatment. The total lengths of each cluster significantly differed from each other (Tukey’s test, q0.05,3000,4 = 3.314; I ≠ II ≠ III) (Table 2). The angles between body segments were also significantly different among clusters. In cluster I, the angles were relatively greater at the end of the body, however, curvatures were greater in the middle part of the body. Patterns in the RSOM were capable of characterizing the sequential line shapes of test organisms before and after treatment.

FigureFigure 4. 4. RecurrentRecurrent self self-organizing-organizing map map (RSOM (RSOM)) applied applied to to the the line line movements movements of of LumbriculusLumbriculus variegatus before ( ) and after ( ) copper sulfate treatments. (a) Classification of line movement tracks in variegatus before () and after (•· ) copper sulfate treatments. (a) Classification of line movement tracks b# inthe the RSOM RSOM and and ( )( dendrogramb) dendrogram for thefor classificationthe classification of RSOM of RSOM units units based based on the on Ward the linkageWard linkage method with Euclidean distance measures. Different sizes of circles indicate the number of patterns assigned to method with Euclidean distance measures. Different sizes of circles indicate the number of patterns each SOM unit. assigned to each SOM unit. The body movements were correspondingly different in different clusters (Figure5). Cluster I showedFigure a sequence 4. Recurrent of large self-organizing and straightened map (RSOM shapes,) applied cluster to II the displayed line movements folding shapes, of Lumbriculus and cluster III representedvariegatus before contracted () and and after folding (•) copper shapes. sulfate As treatments stated above,. (a) Classification body shapes of line in cluster movement II were tracks mostly movementsin the RSOM before and the (b)chemical dendrogram treatment, for the classification whereas those of RSOM in cluster units based III were on the mostly Ward movementslinkage observedmethod after with the Euclidean treatment. distance Therefore, measures. body Different shapes sizes after of the circles chemical indicate treatments the number were of patterns represented by shortassigned body to lengths each SOM with unit. ahigh degree of folding.

Figure 5. The sequential tracks of line movements of test organisms in different clusters defined in the RSOM. The number underneath each figure indicates the legend of a snapshot at 5 s intervals.

Figure 5. The sequentialsequential trackstracks ofof line line movements movements of of test test organisms organisms in in di ffdifferenterent clusters clusters defined defined in the in theRSOM. RSOM. The The number number underneath underneath each each figure figure indicates indicates the legendthe legend of a of snapshot a snapshot at 5 sat intervals. 5 s intervals.

Each movement parameter differed according to different clusters in the RSOM (Figure6). To compare parameters, 1000 samples were randomly selected from each cluster. The total length in cluster I showed the greatest range among the clusters (Figure6a), with length also being relatively Int. J. Environ. Res. Public Health 2020, 17, 2627 10 of 15 greater in cluster II, another cluster belonging to the before-treatment group, compared with cluster III, which represented after treatment. The total lengths of each cluster significantly differed from each other (Tukey’s test, q0.05,3000,4 = 3.314; I , II , III) (Table2). The angles between body segments were also significantly different among clusters. In cluster I, the angles were relatively greater at the end of the body, however, curvatures were greater in the middle part of the body. Patterns in the RSOM were

Int.capable J. Environ. of characterizingRes. Public Health 20 the20, sequential 17, x FOR PEER line REVIEW shapes of test organisms before and after treatment.10 of 15

FigureFigure 6. 6. VariationVariation in in body body segments segments of of LumbriculusLumbriculus variegatus variegatus accaccordingording to to subclusters subclusters matching matching the the RSOM.RSOM. ( (aa)) T Totalotal length and and ( (bb)) angle angle between between body body segments. segments. Error Error bars bars indicate indicate standard standard deviation. deviation. DifferentDifferent lettersletterson on error error bars bars indicate indicate significant significant differences differences between between clusters clusters based on based Tukey’s on honestlyTukey’s hsignificantonestly significant difference difference test (P < 0.05).test (P < 0.05).

TableTable 2. Tukey’s multiplemultiple comparison comparison test test for for homogeneous homogeneous groups groups defined defined in the in self-organizingthe self-organizing map map(SOM) (SOM for angles) for angles characterizing characterizing body segments body segments of Lumbriculus of Lumbriculus variegatus variegatusbefore and before after and 0.1 mgafter/L and0.1 mg/L0.01 mg and/L 0.01 copper mg/L treatments. copper treatments.

a1 a1a2 a2a3 a3 a4a4 a5 a5 F F100.89 100.89 35.16 35.16 25.16 25.16 10.6010.60 5.095.09 P <0.0001 <0.0001 <0.0001 <0.0001 0.0056 P <0.0001 <0.0001 <0.0001 <0.0001 0.0056 Group I , II , III I , II , III I , II , III I , II = III I , II = III Group I ≠ II ≠ III I ≠ II ≠ III I ≠ II ≠ III I ≠ II = III I ≠ II = III a6 a7 a8 a9 a10 a11 a6 a7 a8 a9 a10 a11 F 33.82 22.26 21.03 40.61 55.64 70.96 F P 33.82<0.0001 22.26<0.0001 <21.030.0001 <0.000140.61 <0.000155.64 <0.000170.96 P Group<0.0001 I , II = III<0.0001 I , II = III I<0.0001, II = III I ,<0.0001II = III I , II<0.0001, III I , II ,<0.0001III Group I ≠ II = III I ≠ II = III I ≠ II = III I ≠ II = III I ≠ II ≠ III I ≠ II ≠ III 4. Discussion 4. Discussion We demonstrated that the movement of segments of line-shaped animals (e.g., oligochaetes) could be analyzedWe demonstrated to detect chemical that the stress movement at low ofconcentrations. segments of line Detailed-shaped diff animalserences could (e.g., beoligochaetes) seen in the couldbody segmentsbe analyzed (Figures to detect1 and chemical3). Di ff stresserences at low in response concentrations. behaviors Detailed were also differences observed could according be seen to indi ff theerent body toxin segments concentrations (Figures (Figure 1 and2), 3). Shannon Differences entropy in indices response (Figure behaviors1b), and were angles also (Figure observed3), according to different toxin concentrations (Figure 2), Shannon entropy indices (Figure 1b), and angles (Figure 3), demonstrating that response behaviors were associated with the strength of the toxic impact. It was also noteworthy that differences in responses were represented by differences in angles according to different toxin concentrations. Although angle changes were higher at the center positions at higher concentrations with regard to the total body, the angle changes were differentiated at the end positions at low concentrations. Greater curvature was observed at the center position, indicating sinuousness of the body (Figure 3). However, parameters representing speed in body segments were not sufficiently differentiated to detect stress behaviors of test organisms.

Int. J. Environ. Res. Public Health 2020, 17, 2627 11 of 15 demonstrating that response behaviors were associated with the strength of the toxic impact. It was also noteworthy that differences in responses were represented by differences in angles according to different toxin concentrations. Although angle changes were higher at the center positions at higher concentrations with regard to the total body, the angle changes were differentiated at the end positions at low concentrations. Greater curvature was observed at the center position, indicating sinuousness of the body (Figure3). However, parameters representing speed in body segments were not su fficiently differentiated to detect stress behaviors of test organisms. Although detection of the response behaviors based on body segments may require extra computational resources (e.g., calculating body length and curvature) than movement parameters, body segments appeared to be efficient in determining detailed response behaviors, including partial movements of line-shaped bodies. The angle of the tail significantly changed after treatment with the toxic substance, indicating that the tail behavior of test organisms was more affected than the head region (Figure3b). Line-shaped body visualization would be useful to characterize line-movement patterns of test organisms due to the extra dimension included to express the complexity of response behaviors compared with point data. We used L. variegatus as the test organism in this study. Aquatic oligochaetes, including Lumbriculus, are important taxa in freshwater communities. Oligochaetes play an important role in aquatic communities, aiding in the decomposition of organic materials in sediment. In particular, L. variegatus is sensitive to exterior stimuli and has a diverse scope of behaviors [33,51], and was proposed as a standard organism for sediment bioaccumulation tests [52]. This organism has stereotypical behaviors that can be used for sublethal toxicology [16,53,54], and could be a biological indicator. The species was reported to be more sensitive to toxic chemicals, including copper, compared with other indicator species in Tubificidae, such as Tubifex tubifex [55]. However, line movements were not reported after exposure to stressors prior to this work. To the best of our knowledge, the pattern analysis of line data for oligochaetes, as shown in this study, was not previously conducted in response to toxic chemicals. Time-series analysis by RSOM could be a stepping-stone for the quantitative patterning of continuous movements of elongated animals, such as oligochaetes. We used copper sulfate as a toxicant in this study. Copper is commonly used for fertilizers and pesticides due to its antifungal properties; however, excessive levels can contaminate aquatic ecosystems, even though it is an essential trace element [56]. Monitoring and removal methods were developed due to its toxicity to organisms [57,58]. Regarding efficiency in chemical engineering tools, biomonitoring with observations of response behaviors were not able to measure the precise concentration of toxic materials. However, the risk of environmental toxicity of copper is varied in regard to freshwater organisms. We demonstrated that copper sulfate caused typical symptomatic behaviors in the temporal sequence of line movement of blackworms according to different concentrations. The toxic effects were mostly reported at one concentration in continuous behavioral monitoring [5,13,15]. However, physiological–behavioral studies are needed in the future to reveal the causality of relationships in the production of symptoms. Line tracking of test specimens is typically performed with C. elegans to characterize movement patterns according to different phenotypes [18,20,22]. Although C. elegans is a good indicator for the analysis of behavior according to phenotype, some drawbacks exist regarding monitoring in practical situations. This species is somewhat small and a microscope is usually required for monitoring. Consequently, it would be difficult for monitoring under field conditions, requiring more in situ facilities. In addition, the nematode’s life-span is too short for continuous monitoring under real conditions. However, aquatic Oligochaeta, including Lumbriculus, are reasonable species for monitoring under field conditions because of their large size and long life-span. In this study, the head and tail regions were not specifically differentiated through digital image processing. We manually identified the head and tail regions of oligochaetes in movement segments, and modified methods were used to confirm body parts. Feng et al. [22] distinguished the heads and tails of C. elegans manually on a monitor. In the preliminary study, we attempted to automatically Int. J. Environ. Res. Public Health 2020, 17, 2627 12 of 15 identify the head and tail regions of C. elegans, but it was difficult to automatically distinguish between the head and tail regions of oligochaetes. Because the body was much longer and movement behaviors of oligochaetes (e.g., body reverse swimming) were highly complex compared with those of C. elegans, further development of methods to recognize body parts in elongated animals, such as oligochaetes, is required in the future.

5. Conclusions Analysis of the body segment data of the line-shaped animals was useful for detecting responses. Line-shaped body movements with parameters expressing body segments broadened the scope of behavioral characterization of elongated animals, adding an extra dimension to reveal the complexity of line-shaped bodies compared with point detection in round-shaped species. RSOM is therefore helpful in addressing the sequential line movements of test organisms, especially line-shaped animals.

Author Contributions: Conceptualization, T.-S.C. and Y.-S.P.; methodology, C.W.J., Y.C., H.W., and I.-S.K.; formal analysis, C.W.J. and Y.-S.P.; investigation, C.W.J., T.-S.C., and Y.-S.P.; visualization, C.W.J. and Y.-S.P.; writing—original draft preparation, C.W.J.; writing—Review and editing, C.W.J., Y.C., H.W., I.-S.K., T.-S.C., and Y.-S.P. All authors have read and agreed to the published version of the manuscript. Funding: This study was supported by the National Research Foundation of Korea (NRF) funded by the Korean government (NRF-2016R1A2B4011801, NRF-2018R1A6A1A03024314, NRF-2019R1H1A2080107 and NRF-2019R1A2C1087099). Conflicts of Interest: The authors declare no conflict of interest.

References

1. Lemly, A.D.; Smith, R.J. A behavioral assay for assessing effects of pollutants of fish chemoreception. Ecotoxicol. Environ. Saf. 1986, 11, 210–218. [CrossRef] 2. Gerhardt, A.; Carlsson, A.; Ressemann, C.; Stich, K.P. New online biomonitoring system for Gammarus pulex (L.) (Crustacea): In situ test below a copper effluent in south Sweden. Environ. Sci. Technol. 1998, 32, 150–156. [CrossRef] 3. Teather, K.; Harris, M.; Boswell, J.; Gray, M. Effects of Acrobat MZ®and Tattoo C®on Japanese medaka (Oryzias latipes) development and adult male behavior. Aquat. Toxicol. 2001, 51, 419–430. [CrossRef] 4. Borcherding, J. Ten years of practical experience with the Dreissena-Monitor, a biological early warning system for continuous water quality monitoring. Hydrobiologia 2006, 556, 417–426. [CrossRef] 5. Ji, C.W.; Lee, S.H.; Choi, K.-H.; Kwak, I.-S.; Lee, S.G.; Cha, E.Y.; Lee, S.-K.; Chon, T.-S. Monitoring of movement behaviors of chironomid larvae after exposure to diazinon by using fractal dimension and Self-Organizing Map. Int. J. Ecodyn. 2007, 2, 27–38. [CrossRef] 6. Bae, M.-J.; Park, Y.-S. Biological early warning system based on the responses of aquatic organisms to disturbances: A review. Sci. Total Environ. 2014, 466–467, 635–649. [CrossRef][PubMed] 7. Dutta, H.; Marcelino, J.; Richmonds, C. Brain acetylcholinesterase activity and optomotor behavior in bluegills, Lepomis macrochirus exposed to different concentrations of diazinon. Arch. Int. Physiol. Biochim. Biophys. 1992, 100, 331–334. [PubMed] 8. Chon, T.-S.; Park, Y.-S.; Ross, M.H. Activity of German cockroach, Blattella germanica (L.) (Orthoptera : Blattellidae), at different microhabitats in semi-natural conditions when treated with sublethal doses of pesticides. J. Asia-Pac. Entomol. 1998, 1, 99–107. 9. Chon, T.-S.; Park, Y.-S.; Park, J.Y.; Choi, S.-Y.; Kim, K.T.; Cho, E.C. Implementation of computational methods to pattern recognition of movement behavior of Blattella germanica (Blattaria : Blattellidae) treated with Ca2+ signal inducing chemicals. Appl. Entomol. Zool. 2004, 39, 79–96. [CrossRef] 10. Roast, S.D.; Widdows, J.; Jones, M.B. Disruption of swimming in the hyperbenthic mysid Neomysis integer (Peracarida: Mysidacea) by the organophosphate pesticide chlorpyrifos. Aquat. Toxicol. 2000, 47, 227–241. [CrossRef] 11. Ibrahim, W.L.F.; Furu, P.; Ibrahim, A.M.; Christensen, N.Ø. Effect of the organophosphorous insecticide, chlorpyrifos (Dursban), on growth, fecundity and mortality of Biomphalaria alexandrina and on the production of Schistosoma mansoni cercariae in the snail. J. Helminthol. 1992, 66, 79–88. [CrossRef][PubMed] Int. J. Environ. Res. Public Health 2020, 17, 2627 13 of 15

12. Gray, M.A.; Teather, K.L.; Metcalfe, C.D. Reproductive success and behavior of Japanese Medaka (Oryzias latipes) exposed to 4-tert-octylphenol. Environ. Toxicol. Chem. 1999, 18, 2587–2594. 13. Park, Y.-S.; Chung, N.-I.; Choi, K.-H.; Cha, E.Y.; Lee, S.-K.; Chon, T.-S. Computational characterization of behavioral response of medaka (Oryzias latipes) treated with diazinon. Aquat. Toxicol. 2005, 71, 215–228. [CrossRef][PubMed] 14. Ji, C.W.; Lee, S.H.; Kwak, I.-S.; Cha, E.Y.; Lee, S.-K.; Chon, T.-S. Computational analysis of movement behaviors of medaka (Oryzias latipes) after the treatments of copper by using fractal dimension and artificial neural networks. In Environmental Toxicology; Kungolos, A.G., Brebbia, C.A., Samaras, C.P., Popov, V., Eds.; WIT Press: Southampton and Mykonos, UK, 2006; pp. 93–107. 15. Kwak, I.-S.; Chon, T.-S.; Kang, H.-M.; Chung, N.-I.; Kim, J.-S.; Koh, S.C.; Lee, S.-K.; Kim, Y.-S. Pattern recognition of the movement tracks of medaka (Oryzias latipes) in response to sub-lethal treatments of an insecticide by using artificial neural networks. Environ. Pollut. 2002, 120, 671–681. [CrossRef] 16. Gerhardt, A. Importance of exposure route for behavioural responses in Lumbriculus variegatus Müller (Oligochaeta: Lumbriculida) in short-term exposures to Pb. Environ. Sci. Pollut. Res. 2007, 14, 430–434. [CrossRef] 17. Clark, D.A.; Gabel, C.V.; Gabel, H.; Samuel, A.D.T. Temporal activity patterns in thermosensory neurons of freely moving Caenorhabditis elegans encode spatial thermal gradients. J. Neurosci. 2007, 27, 6083–6090. [CrossRef] 18. Baek, J.-H.; Cosman, P.; Feng, Z.; Silver, J.; Schafer, W.R. Using machine vision to analyze and classify Caenorhabditis elegans behavioral phenotypes quantitatively. J. Neurosci. Methods 2002, 118, 9–21. [CrossRef] 19. Geng, W.; Cosman, P.; Berry, C.C.; Feng, Z.; Schafer, W.R. Automatic tracking, feature extraction and classification of C. elegans phenotypes. IEEE Trans. Biomed. Eng. 2004, 51, 1811–1820. [CrossRef] 20. Stephens, G.J.; Johnson-Kerner, B.; Bialek, W.; Ryu, W.S. Dimensionality and dynamics in the behavior of C. elegans. PLoS Comput. Biol. 2008, 4, 10. [CrossRef] 21. Hardaker, L.A.; Singer, E.; Kerr, R.; Zhou, G.; Schafer, W.R. Serotonin modulates locomotory behavior and coordinates egg-laying and movement in Caenorhabditis elegans. J. Neurobiol. 2001, 49, 303–313. [CrossRef] 22. Feng, Z.; Cronin, C.J.; Wittig, J.H.J.; Sternberg, P.W.; Schafer, W.R. An imaging system for standardized quantitative analysis of C. elegans behavior. BMC Bioinform. 2004, 5, 1–6. [CrossRef][PubMed] 23. Ajie, B.C.; Estes, S.; Lynch, M.; Phillips, P.C. Behavioral degradation under mutation accumulation in Caenorhabditis elegans. Genetics 2005, 170, 655–660. [CrossRef][PubMed] 24. Eom, H.-J.; Kim, H.; Kim, B.-M.; Chon, T.-S.; Choi, J. Integrative assessment of benzene exposure to Caenorhabditis elegans using computational behavior and toxicogenomic analyses. Environ. Sci. Technol. 2014, 48, 8143–8151. [CrossRef][PubMed] 25. Bilbao, A.; Patel, A.K.; Rahman, M.; Vanapalli, S.A.; Blawzdziewicz, J. Roll maneuvers are essential for active reorientation of Caenorhabditis elegans in 3D media. Proc. Natl. Acad. Sci. USA 2018, 115, E3616. [CrossRef] [PubMed] 26. Awual, M.R. New type mesoporous conjugate material for selective optical copper(II) ions monitoring & removal from polluted waters. Chem. Eng. J. 2017, 307, 85–94. [CrossRef] 27. Awual, M.R.; Hasan, M.M.; Rahman, M.M.; Asiri, A.M. Novel composite material for selective copper(II) detection and removal from aqueous media. J. Mol. Liq. 2019, 283, 772–780. [CrossRef] 28. Awual, M.R.; Hasan, M.M.; Khaleque, M.A.; Sheikh, M.C. Treatment of copper(II) containing wastewater by a newly developed ligand based facial conjugate materials. Chem. Eng. J. 2016, 288, 368–376. [CrossRef] 29. Petersen, E.J.; Huang, Q.; Weber, W.J. Ecological uptake and depuration of carbon nanotubes by Lumbriculus variegatus. Environ. Health Perspect. 2008, 116, 496–500. [CrossRef] 30. US EPA. Methods for Measuring the Toxicity and Bioaccumulation of Sediment-associated Contaminants with Freshwater Invertebrates; U.S. Environmental Protection Agency: Washington, DC, USA, 2000. 31. Phipps, G.L.; Ankley, G.T.; Benoit, D.A.; Mattson, V.R. Use of the aquatic oligochaete Lumbriculus varieatus for assessing the toxicity and bioaccumulation of sediment-associated contaminants. Environ. Toxicol. Chem. 1993, 12, 269–279. [CrossRef] 32. O’Gara, B.A.; Bohannon, V.K.; Teague, M.W.; Smeaton, M.B. Copper-induced changes in locomotor behaviors and neuronal physiology of the freshwater Oligochaete, Lumbriculus variegatus. Aquat. Toxicol. 2004, 69, 51–66. [CrossRef] Int. J. Environ. Res. Public Health 2020, 17, 2627 14 of 15

33. Drewes, C.D. Helical swimming and body reversal behaviors in Lumbriculus variegatus (Annelida: : Lumbriculidae). Hydrobiologia 1999, 406, 263–269. [CrossRef] 34. Gerhardt, A.; de Bisthoven, L.J.; Penders, E. Quality control of drinking water from the River Rhine with the multispecies freshwater biomonitor. Aquat. Ecosyst. Health Manag. 2003, 6, 159–166. [CrossRef] 35. Kim, C.-K.; Kwak, I.-S.; Cha, E.-Y.; Chon, T.-S. Implementation of wavelets and artificial neural networks to detection of toxic response behavior of chironomids (Chironomidae: Diptera) for water quality monitoring. Ecol. Model. 2006, 195, 61–71. [CrossRef] 36. Son, K.H.; Ji, C.W.; Park, Y.-M.; Cui, Y.; Wang, H.Z.; Chon, T.-S.; Cha, E.Y. Recurrent Self-Organizing Map implemented to detection of temporal line-movement patterns of Lumbriculus variegates (Oligochaeta: Lumbriculidae) in response to the treatments of heavy metal. In Environmental Toxicology; Kungolos, A.G., Brebbia, C.A., Samaras, C.P., Popov, V., Eds.; WIT Press: Southampton/Mykonos, UK, 2006; pp. 93–107. 37. Ji, C.W.; Kim, H.; Park, Y.-S.; Chon, T.-S. Computational analysis of movement behaviors of medaka (Oryzias latipes) in response to chemical and thermal stressors. J. Korean Phys. Soc. 2012, 60, 570–575. [CrossRef] 38. Gershenfeld, N.; Weigend, A. The future of time series: Learning and understanding. In Time Series Prediction: Forecasting the Future and Understanding the Past; Weigend, A., Gershnfeld, N., Eds.; Addison-Wesley: Reading, MA, USA, 1993; pp. 1–70. 39. Chon, T.S.; Park, Y.S.; Cha, E.Y. Patterning of community changes in benthic macroinvertebrates collected from urbanized streams for the short time prediction by temporal artificial neuronal networks. In Artificial Neuronal Networks; Lek, S., Guégan, J.-F., Eds.; Springer: Berlin/Heidelberg, Germany, 2000; pp. 99–114. [CrossRef] 40. Chon, T.-S.; Kwak, I.-S.; Park, Y.-S.; Kim, T.-H.; Kim, Y. Patterning and short-term predictions of benthic macroinvertebrate community dynamics by using a recurrent artificial neural network. Ecol. Model. 2001, 146, 181–193. [CrossRef] 41. Kohonen, T. Self-Organizing Maps; Springer: Berlin, Germany, 2001. 42. Varsta, M.; Heikkonen, J.; Millan, J.d.R. Context learning with the self organizing map. In Proceedings of the Proceedings of WSOM 1997, Workshop on Self-Organizing Maps, Espoo, Finland, 4–6 June 1997; pp. 197–202. 43. Timo, K.; Varsta, M.; Heikkonen, J.; Kaski, K. Temporal sequence processing using recurrent SOM. In Proceedings of the 2nd International Conference on Knowledge-Based Intelligent Engineering Systmes, Adelaide, Australia, 21 April 1998; pp. 290–297. 44. Angeloviˇc,P. Time series prediction using RSOM and local models. In Proceedings of the IIT. SRC 2005 Student Research Conference, Bratislava, Slovakia, 27 May 2005; pp. 27–34. 45. Bailey, H.C.; Liu, D.H.W. Lumbriculus variegatus, a benthic Oligochaete, as a bioassay organism. In Aquatic Toxicology Hazard Assessment 3rd Symposium, ASTM STP 707; Eaton, J.C., Parrish, P.R., Hendricks, A.C., Eds.; ASTM: Philadelphia, PA, USA, 1980. 46. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423, 623–656. [CrossRef] 47. StatSoft Inc. STATISTICA (Data Analysis Software System), Version 7; StatSoft Inc.: Tulsa, Oklahoma, 2004. 48. The Mathworks. MATLAB, Version 7.0; The Mathworks, Inc.: Massachusetts, MA, USA, 2001. 49. Mielke, E.W.; Berry, K.J.; Johnson, E.S. Multiresponse permutation procedures for a priori classifications. Commun. Stat. Theory Methods 1976, 5, 1409–1424. [CrossRef] 50. McCune, B.; Mefford, M.J. PC-ORD. Multivariate Analysis of Ecological Data. Version 4.41; MjM Software: Gleneden Beach, OR, USA, 1999. 51. Lesiuk, N.M.; Drewes, C.D. Autotomy reflex in a freshwater oligochaete, Lumbriculus variegatus (Clitellata: Lumbriculidae). Hydrobiologia 1999, 406, 253–261. [CrossRef] 52. ASTM (American Society for Testing and Materials). Annual Book of ASTM Standards. Volume 11.05; ASTM: Philadelphia, PA, USA, 1995. 53. Rogge, R.W.; Drewes, C.D. Assessing sublethal neurotoxicity effects in the freshwater oligochaete, Lumbriculus variegatus. Aquat. Toxicol. 1993, 26, 73–90. [CrossRef] 54. Ding, J.; Drewes, C.D.; Hsu, W.H. Behavioral effects of ivermectin in a freshwater oligochaete, Lumbriculus variegatus. Environ. Toxicol. Chem. 2001, 20, 1584–1590. [CrossRef] 55. Chapman, K.K.; Benton, M.J.; Brinkhurst, R.O.; Scheuerman, P.R. Use of the aquatic oligochaetes Lumbriculus variegatus and Tubifex tubifex for assessing the toxicity of copper and cadmium in a spiked-artificial-sediment toxicity test. Environ. Toxicol. 1998, 14, 271–278. [CrossRef] Int. J. Environ. Res. Public Health 2020, 17, 2627 15 of 15

56. GERBE. Programme de Recherche: Toxicologie - Ecotoxicologie des Pesticides et des Metaux Lourds. Deuxieme Rapport d’Activites; Laboratory of Eco-Toxicology, Faculty of Sciences, University of Reims Champagne: Ardenne, France, 1998. 57. Awual, M.R. Novel ligand functionalized composite material for efficient copper(II) capturing from wastewater sample. Compos. Part B Eng. 2019, 172, 387–396. [CrossRef] 58. Awual, M.R. A novel facial composite adsorbent for enhanced copper(II) detection and removal from wastewater. Chem. Eng. J. 2015, 266, 368–375. [CrossRef]

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).