- 1Kagawa School of Pharmaceutical Sciences, Tokushima Bunri University, Sanuki, Japan
- 2Department of Biology, Waseda University, Tokyo, Japan
- 3Research Organization for Nano and Life Innovation, Waseda University, Tokyo, Japan
- 4Hotchkiss Brain Institute, University of Calgary, Calgary, AB, Canada
- 5WASEDA Bioscience Research Institute in Singapore, Singapore, Singapore
- 6Lipid Science and Aging Research Center and Center for Stem Cell Research, Kaohsiung Medical University, Kaohsiung, Taiwan
The pond snail Lymnaea stagnalis learns and forms long-term memory (LTM) following both operant conditioning of aerial respiratory behavior and classical conditioning of taste aversive behavior. In the present study, we examined whether there are interstrain differences in the ability to form LTM following these two types of conditioning. A strain of Lymnaea (TC1) collected in Alberta, Canada exhibits superior memory-forming ability following aerial respiratory operant conditioning compared to a laboratory-reared strain of Lymnaea from Netherlands known as the Dutch strain. We asked whether the offspring of the Canadian TC1 and Dutch snails (i.e., filial 1 (F1) cross snails) would have the superior memory ability and found, rather, that their memory ability was average like the Dutch snails. That is, the Canadian TC1 snails have superior ability for LTM formation following aerial respiratory operant conditioning, but the Dutch and the generated F1 cross have average ability for memory forming. We next examined the Canadian TC1, Dutch and F1 cross snails for their ability to learn and form memory following conditioned taste aversion (CTA). All three populations showed similar associative CTA responses. However, both LTM formation and the ratio of good-to-poor performers in the memory retention test were much better in the Dutch snails than the Canadian TC1 and F1 cross snails. The memory abilities of the Canadian TC1 and F1 cross snails were average. Our present findings, therefore, suggest that snails of different strains have different memory abilities, and the F1 cross snails do not inherit the memory ability from the smart strain. To our knowledge, there have been a limited number of studies examining differences in memory ability among invertebrate strains, with the exception of studies using mutant flies.
Introduction
The pond snail Lymnaea stagnalis is a useful animal model for investigating the causal, neuronal mechanisms of learning and memory from the behavioral to molecular levels (Lukowiak and Syed, 1999; Benjamin et al., 2000; Hatakeyama et al., 2006; Kojima et al., 2015; Sunada et al., 2017a,b). Lymnaea can be both operantly (aerial respiratory behavior) and classically (feeding behavior) conditioned (Kemenes and Benjamin, 1989; Kojima et al., 1996; Lukowiak et al., 1996). When snails placed in an aquatic, hypoxic environment received a tactile stimulus to the pneumostome area every time they attempted to breathe, the snails learned not to breathe, and the number of stimuli received in a memory test session was significantly lower than that in the training session (Lukowiak et al., 1996, 2000; Lukowiak and Syed, 1999). The early studies performed by the Lukowiak laboratory utilized an inbred laboratory-reared strain of Lymnaea from Netherlands (hereinafter, the Dutch strain; see Van Der Steen et al., 1969) for the aerial respiratory operant conditioning. In the Dutch strain, two 0.5 h training sessions with a 1 h interval between them were required for the snails to form long-term memory (LTM) persisting for at least 24 h (Lukowiak et al., 2000). However, some naturally occurring strains of Lymnaea form memory faster and better than the Dutch strain (Orr et al., 2008, 2009; Dalesman and Lukowiak, 2012). In these so-called smart snails (e.g., Canadian TC1 strain; see below for details), LTM is formed after only a single 0.5 h training session (Shymansky et al., 2017).
Having shown that there are differences in learning and memory-forming abilities following aerial respiratory operant conditioning, in a recent study we sought to begin to determine the cost of being smart. Hughes et al. (2017) showed that one apparent cost of possessing the smart phenotype was that the smart snails were less able to cope with stress than average snails. The stressors used in the Hughes et al. (2017) study were thermal, resource restriction combined with food odor, predator detection and tissue injury (shell damage). Their results suggested that a stressor or a combination of stressors act to enhance memory in average snails but obstruct memory formation in smart snails. These results are consistent with the Yerkes-Dodson/Hebb law, resulting that smart snails are more easily stressed than average snails (Yerkes and Dodson, 1908; Hebb, 1955; Ito et al., 2015).
In the present study, we ask whether there are other costs or benefits of possessing the smart phenotype as determined by the studies of operant conditioning of aerial respiration (e.g., Hughes et al., 2016; Shymansky et al., 2017). Previously, attempts have been made to study this question using the fruit fly Drosophila. However, almost all these studies utilized different mutants (Tully, 1996). It may be better to examine this question using different strains of a single animal species rather than mutants. This is because we are able to perform behavioral experiments on different known strains and then to find the causes, for example genetic drift. Genetic drift has been suggested to account for significant differences in European Lymnaea stagnalis populations (e.g., how they respond to different levels of copper in the water) that occur (Puurtinen et al., 2004a,b, 2007; Côte et al., 2015). Changes in learning and memory forming ability due to genetic drift possibly have a more important biological meaning than those due to mutation. Thus, we used different strains of Lymnaea to better answer these questions.
Two established strains of Lymnaea were used in the present study: the Dutch strain and the Canadian TC1 strain. As a third population, the filial 1 (F1) generation resulting from a cross of the Dutch and Canadian TC1 snails was also used. These three populations were bred in a laboratory (i.e., not freshly collected). The Dutch strain was isolated in the 1950s from freshly collected snails (see Van Der Steen et al., 1969). The Canadian TC1 strain was collected in 2010 (Braun et al., 2012) and has been bred in Alberta, Canada and in Kagawa and Tokyo for many generations. The F1 cross population was bred in Tokyo, Japan.
Studies on the classical conditioning of taste aversive behavior (i.e., conditioned taste aversion, CTA) in Dutch snails were previously performed at the Ito laboratory. In CTA experiments, an appetitive stimulus (e.g., sucrose) is used as the conditioned stimulus (CS), and an application of the CS to the lips increases the feeding response (i.e., the number of bites). An aversive stimulus (e.g., electric chock) is used as the unconditioned stimulus (US), and the US causes the snails to immediately cease feeding. In the CTA-training procedure, the CS is paired with the US. After repeated CS-US temporal pairings, the CS no longer elicits feeding, and this CTA persists for at least a month (Kojima et al., 1996; Ito et al., 1999, 2013).
Recent studies at the Ito laboratory have used an automated learning apparatus for CTA (see “Materials and Methods” Section for its details; Takigami et al., 2016). We have obtained good results in Dutch snails using a 5-s presentation of a 100 mM sucrose (the concentration in the reservoir) as the CS and a 0.2-s presentation of a 9 V electric shock as the US. The pond water always flowed at the speed of 200 ml/min, whereas the sucrose solution flowed at 1 ml/s. Thus, the concentration of sucrose in the reservoir is higher than the concentration sensed by the snails because the sucrose solution was diluted with the continuous pond water flow. On the other hand, no attempts have been made to examine the learning ability of CTA in the Canadian TC1 snails.
Lymnaea stagnalis is a species distributed across Europe and North America (Mozley, 1939). A mitochondrial ribosomal RNA study suggested that there are genetically distinct populations of Lymnaea stagnalis, including in Germany and Italy (Remigio and Blair, 1997). The sequence divergence and genetic distance between these two populations (i.e., German and Italian populations) were greater than those between some separate species of snails. Support for such a conclusion came from a series of studies by a Finnish group studying Lymnaea stagnalis in eight different ponds/shallow bays of larger lakes within 20 km of each other in Finland (Puurtinen et al., 2004a,b, 2007). They concluded that there is detectable genetic variability between these populations. This conclusion was also supported in a more recent study concerning the genetic variation in copper tolerance of Lymnaea stagnalis (Côte et al., 2015). Random genetic drift was hypothesized to explain the observed genetic divergence. A similar conclusion was drawn in a study of 14 different populations of Lymnaea stagnalis collected in Belgium, Netherlands and Germany (Bouétard et al., 2014). Thus, rather than studying differences in learning and memory in mutants (Tully, 1996), we propose to study interstrain differences that are thought to have arisen through normal genetic drift.
In the present study, we both operantly condition and classically condition the above-mentioned two strains and their common offspring (F1 cross snails) of Lymnaea stagnalis: the Dutch strain, the Canadian TC1 strain and the F1 cross snails resulting from their mating. We measure their abilities to learn and form memory, and we examine whether the strain that exhibits better performance in the aerial respiratory operant conditioning task also exhibits better ability to form LTM following the CTA classical conditioning. That is, we examine whether the same strain is “smarter” in terms of both operant and classical conditioning. Further, we established their common offspring (F1 cross snails) to examine whether the F1 cross snails inherit the memory ability from the smart strain or not.
Materials and Methods
Snails
We used two known strains and an F1 cross obtained by mating snails from the two strains (F1 cross snails) of Lymnaea stagnalis (Linnaeus, 1758) in the present study. The first was the Dutch strain, the standard laboratory strain used in studies performed worldwide. This strain was established from snails collected in a polder in Utrecht, Netherlands in the 1950s and continues to be maintained at the Vrije Universiteit Amsterdam in Netherlands (see Van Der Steen et al., 1969). The populations of Dutch snails used in the present study were established in Canada and Japan using snails derived from the Vrije Universiteit Amsterdam colony in 1980s and 1990s, respectively. The second strain used was the Canadian TC1 strain, which was collected from a pond adjacent to the Trans-Canada Highway, Alberta, Canada (at 51.07 N, 114.39 W; see Braun et al., 2012). This Canadian TC1 population was transferred to Japan and maintained separately in dechlorinated tap water as a substitute for pond water under a 12:12 light-dark cycle at around 20°C. Finally, the third was derived from eggs laid by mating the Dutch and Canadian TC1 snails (see below). We termed this population the F1 cross. All the snails were fed ad libitum on a kind of turnip leaf (Brassica rapa var. peruviridis: Komatsuna [in Japanese]) and a spiral shell food (Nisso, Saitama, Japan) every other day (Otsuka et al., 2013). The snails with a <5 mm shell were fed with goldfish food originated from vegetables (Kyorin: Hiraki, Himeji, Japan). The snails with a ca. 20 mm shell length were used for the behavioral experiments (Yamagishi et al., 2015). The generation time for these two strain snails and their offspring snails is almost 2 years.
F1 Cross Snails
The Dutch and the Canadian TC1 snails were reared in separate aquaria as described above. One immature Lymnaea individual with a ca. 10 mm shell length was selected from each tank, and these two snails were reared together in a separate aquarium. Because snails do not become sexually mature until obtaining a shell length of 15 mm, these two snails sexually immature (Yamanaka et al., 1999; Sadamoto et al., 2000). Copulation between the snails was observed, and selfing was not observed. Thus, the eggs laid were the result of the mating between the two strains. The hatchlings in these aquaria were designated the F1 cross snails. Further, we obtained the F1 cross snails from the combination of different Dutch snails and different Canadian TC1 snails. That is, we purposely obtained F1 cross snails with great genetic variable. In preliminary experiments, we had already confirmed that F1 cross snails can bear eggs and that the fry hatched from these eggs can be normally grown.
Aerial Respiratory Operant Conditioning
The snails were kept in aerated artificial pond water (a 1% concentration of Instant Ocean; Spectrum Brands, Middleton, WI, USA; with 80 mg/l CaCl2) for the aerial respiratory operant conditioning. The training procedure was identical to that described previously (Dalesman and Lukowiak, 2012). To increase the occurrence of aerial respiration, the artificial pond water in a beaker was made hypoxic (≤5% O2) by vigorously bubbling with N2 for 20 min before training commenced; bubbling was then continued at a reduced rate throughout training to maintain hypoxic conditions (Figure 1A). Snails were replaced into the breaker and allowed to acclimate for 10 min before the training session. This acclimation period was then followed by a 0.5 h training period. During the training period, each time a snail attempted to open its pneumostome (i.e., breathing tube) at the water’s surface, the pneumostome was gently poked using a sharpened wooden stick. This caused the snail to withdraw its pneumostome, but not its whole body to withdraw. To test for LTM, the protocol used for 0.5 h training session was repeated at 24 h after the training session.
Figure 1. Training procedures for aerial respiratory operant conditining and conditioned taste aversion (CTA) in Lymnaea. (A) Aerial respiratory operant conditining. The pond water was made hypoxic (≤5% O2) by bubbling with N2 before training, and bubbling was then continued throughout training. Snails were acclimatized for 10 min before the training session. During the 0.5 h training session, each time a snail attempted to open its pneumostome at the water’s surface, the pneumostome was gently poked. For memory test, the protocol used for 0.5 h training session was repeated at 24 h after the training session. (B,C) CTA. (B) For forward conditioning (i.e., CTA training), snails were subjected to 10 pairings of sucrose administration (100 mM, conditioned stimulus, CS) for 5 s and electric shock (9 V, 0.4 μA, unconditioned stimulus, US) for 0.2 s. The inter-trial interval was 65.2 s. For backward conditioning, the stimuli were reversed. That is, 10 pairings of application of electric shock (US) for 0.2 s and application of sucrose (CS) for 5 s were used. The inter-trial interval was also 65.2 s in the backward confitioning. For naive control experiments, the CS and the US were replaced with the application of DW. The timng was the same as that of the forward training. (C) The time schedule is shown for training. Snails were deprived of food for 1 day before training (except Figure 5). A pretest and two post-tests were performed to count the number of feeding responses to sucrose (bites/min) for 1 min after the application of CS for 5 s. One post-test was performed 10 min after training, and the other post-test was done 24 h after training.
CTA Classical Conditioning
We used a de novo, automatic training apparatus that was developed at the Sakakibara laboratory of Tokai University (Takigami et al., 2016), with slight modifications. In this conditioning apparatus, the pond water, which was propelled by a water pump, constantly flowed on snails at a speed of 200 ml per min. We could thus easily supplemented the flowing water with 100 mM sucrose (the CS) and delivered it to the snails for 5 s. The flow of sucrose (1 ml/s) into the test chamber began immediately when the pump was turned on. The amount of sucrose administered to the mouth of an individual snail in this manner was estimated to be 20–30 mM. Fluids were drained from the chamber via an overflow pathway. The US consisted of a 9 V, 0.4 μA electric shock delivered by a dry cell battery and applied for 0.2 s via a 2-channel electronic stimulator (DPS-1200D; Dia Medical System, Tokyo, Japan) or a microcomputer (Arduino Holdings, New York, NY, USA). The US was strong enough to immediately terminate the feeding behavior before the snail withdrew into its shell, but it was not harmful to the snail (Ito et al., 2015). Both the CS and US were delivered to the snail from signals generated from an electronic stimulator. This stimulator can generate any arbitral timing signals with any duration, delay and repetition number. The snail’s feeding behavior was observed by a human experimenter. Each test chamber consisted of a 50 ml culture flask bottle equipped with an in and out flow path.
Snails were first challenged with a pretest by the application of sucrose (CS; Murakami et al., 2013). The number of bites was recorded in the 1-min interval in pond water after the 5-s application of the CS. Twenty-four hours later, the snails were again subjected to the conditioning and control procedures in the same apparatus as used for the pretest. In the CTA training procedure, we paired the CS with the US, and the snails received 10 paired presentations of the CS-US (forward; Figure 1B). Controls included a backward-conditioned (US-CS) group and a naive group to validate associative learning (Figure 1B). For the naive control group, only pond water was applied to the lips instead of the CS and US. In every instance, the inter-trial interval was set as 65.2 s, because snails required this time to recover from the withdrawal response caused by the electrical stimulus (i.e., the US). The conditioning procedure was completed within 1 day. In the post-test sessions, snails were again challenged with the CS, and the number of bites was counted in the 1-min interval in pond water after the 5-s application of the CS. The time schedule of these protocols is shown in Figure 1C.
Based on our previous findings (Sugai et al., 2007), we set a performance cutoff in the post-test sessions to distinguish between good and poor performers. A snail possessing good learning and LTM (i.e., a good performer) will not open its mouth (i.e., a bite) following presentation of the CS. However, some snails bite by chance (i.e., spontaneously) in the absence of any delivered stimulus (Kojima et al., 1996). Such spontaneous openings occur at a rate of about one per minute. Thus, we defined a good performer as a snail that made 0–1 bites/min during the post-test session in response to presentation of the CS. Poor performers were thus defined as snails that made ≥2 bites/min in response to the CS during the post-test session. In other words, if a snail completely learns and forms memory of CTA, it does not bite ≥2 times/min in response to the sucrose CS. We thus set the threshold at 1 bites/min to distinguish a good performer from a poor performer. The behavioral experiments were performed in the morning, because it has been shown that snails exhibit better learning ability in the morning than at other times (Wagatsuma et al., 2004).
For both types of conditioning (i.e., aerial respiratory operant conditioning and CTA classical conditioning), the training procedures used would allow us to easily distinguish between snails exhibiting enhanced cognitive ability (Shymansky et al., 2017). Those that did not exhibit such enhancement were termed average snails (Hughes et al., 2017).
Data Analysis
The data are expressed as the mean ± SEM. Significant differences were examined at P < 0.05. A paired t-test was performed between the number of attempted openings in the training session and that in the memory test in the operant conditioning of aerial respiration. For there to be memory present, the number of attempted pneumostome openings in the memory test session had to be significantly less (i.e., P < 0.05) than the number of attempted openings in the single 0.5 h training session. Following the CTA training procedure, a one-way ANOVA followed by a post hoc Scheffé’s test or a two-way repeated measures ANOVA having one between subjects factor (strain: Dutch snails, Canadian TC1 snails and F1 cross snails) and one within-subjects factor (time course: pretest, 10 min post-test and 24 h post-test) followed by a post hoc Holm’s test was used to determine whether memory was present (i.e., whether there was a significant difference in the number of bites compared to the pretest response; P < 0.05). When the variances of within-subjects factor were not equal, the degrees of freedom were adjusted by Greenhouse-Geisser’s epsilon. The ratio of good to poor performers was estimated by a χ2 test to detect significant differences at P < 0.05. The computer software used was R (version 3.3.1).
Results
Aerial Respiratory Operant Conditioning
First, we performed experiments to show the difference between the smart Canadian TC1 snails and the average Dutch snails in forming memory following operant conditioning of aerial respiratory behavior. In each cohort of snails (N = 60 each), the snails were subjected to a 0.5 h training session, then an interval of 24 h and a 0.5 h memory test. The Dutch strain was unable to form LTM following the single training session (i.e., the number of attempted openings in the memory test was not significantly different than that in the training session), whereas the Canadian TC1 strain exhibited LTM (i.e., there was a significant difference in the number of attempted openings in the memory test compared to the training session; P < 0.01). To easily compare the data, we plotted the ratios of a change as a percentage of the initial responses in the training session. These data are presented in Table 1 and Figure 2 and are consistent with the previous published data (e.g., Braun et al., 2012).
Table 1. Comparison of learning and memory-forming abilities in two strain snails and their offsprings in Lymnaea stagnails.
Figure 2. Comparison of the pneumostome opening attempts between the training session (training) and the memory test (test) in the Dutch snails (black, N = 60), the Canadian TC1 snails (white, N = 60) and the F1 cross snails (hatched, N = 28) in the operant conditioning of aerial respiratory behavior. The number of pneumostome opening attempts was counted over 0.5 h during the training session and the memory test 24 h later. The data are shown in percentages, indicating that the results of F1 cross snails resembled those of the Dutch snails. We calculated the statistical significant using the raw data. **P < 0.01; N.S., not-significant. The statistical analysis was performed usign a paired t-test.
We next examined the behavioral phenotype of the F1 cross snails following aerial respiratory operant conditioning. We found that the F1 cross snails did not show LTM formation following the single 0.5 h training session (P > 0.05; Figure 2). That is, the F1 cross snails showed an average ability for aerial respiratory operant conditioning. However, as can be seen, the number of attempted openings in the training session of the F1 cross cohort of snails (N = 28) was much smaller than the number of attempted openings in either the Dutch or the Canadian TC1 cohorts. These data show that whereas the Canadian TC1 snails exhibited the smart phenotype, neither the Dutch nor the F1 cross snails met the criteria to be considered smart.
CTA Classical Conditioning
We next examined how the three cohorts of snails (i.e., the Dutch, the Canadian TC1 and the F1 cross snails) performed in a CTA trial. Previous reports demonstrated that the best learning and memory scores were obtained when the snails were food-deprived for 1 day before training (Sugai et al., 2007; Mita et al., 2014). Thus, in the first experiment using the CTA training procedure, we food-deprived all the snails for 1 day. In the pretest, the number of bites was recorded after application of the sucrose CS. There were no significant differences (P > 0.05) in the feeding responses elicited by the CS among the three cohorts of snails (i.e., the Dutch, the Canadian TC1 and the F1 cross snails) in the pretest session (Figure 3A). That is, each of the cohorts responded to the sucrose stimulus in a similar fashion. We then tested the three cohorts 10 min after CTA training (i.e., the 10 min post-test). The CS-elicited feeding responses in all three cohorts were significantly decreased in comparison with those of the pretest (P < 0.01). In addition, there were no significant differences in the number of feeding responses elicited by the CS among the three cohorts (P > 0.05, Figure 3A). That is, all three cohorts showed similar associative learning following CTA conditioning.
Figure 3. Comparison of the feeding responses to the sucrose CS among the Dutch snails (black), the Canadian TC1 snails (white) and the F1 cross snails (hatched) in the classical conditioning trial of taste aversion behavior. The feeding response was examined at pretest, 10 min post-test and 24 h post-test in response to the CS. (A) Forward conditioning (i.e., CTA training) for the three cohorts (i.e., the Dutch, the Canadian TC1 and the F1 cross snails) was performed. The results showed that (1) there was a significant difference in the feeding responses shown in (A) (two-way repeated measures ANOVA: for strains F(2,91) = 3.84, P < 0.05; for time course F(2,182) = 118.18, P < 0.01; for interaction F(4,182) = 2.66, P < 0.05). The simple effect test showed a significangt difference for strains at 24 h post-test (F(2,91) = 5.04, P < 0.01). Further, a Holm’s multiple comparison showed *P < 0.05 between the Dutch and F1 cross snails; *P < 0.05 between the Dutch and Canadian TC1; and P > 0.05 between the Canadian TC1 and F1 cross snails. (2) There were no significant differences in the feeding response to the sucrose CS at the pretest (F(2,91) = 6.07, P > 0.05). (3) The feeding responses to the CS were significantly suppressed in all three cohorts (for example, P < 0.01 for the Canadian TC1 snails), but there were no significant differences among the three cohorts at the 10 min post-test (F(2,91) = 1.27, P > 0.05). (B) Backward conditioning (pairings of the US first and then the CS) for the Dutch, the Canadian TC1 and the F1 cross snails showed that (1) there were no significant differences among the three cohorts at the pretest, the 10 min post-test and the 24 h post-test (two-way repeated measures ANOVA: for strains F(2,52) = 0.42, P > 0.05); but (2) there were significant differeces in the feeding responses shown in B (two-way repeated measures ANOVA: for time course F(2,79.32) = 9.22, P < 0.01. Further, a Holm’s multiple comparison showed *P < 0.05 between the pretest and 10 min post-test; *P < 0.05 between the pretest and 10 min post-test; *P < 0.05 between the 10 min post-test and 24 h post-test). There were no significant interaction (two-way repeated measures ANOVA: for interaction F(3.05,79.32) = 0.81, P > 0.05). (C) The naive snails also showed that (1) there were no significant differences among the three cohorts at the pretest, the 10 min post-test and the 24 h post-test (two-way repeated measures ANOVA: for strain F(2,52) = 0.81, P > 0.05). (2) There were no significant differences in the feeding responses (two-way repeated measures ANOVA: for time course F(1.57,81.69) = 2.95, P > 0.05). There were no significant interaction (two-way repeated measures ANOVA: for interaction F(3.05,79.32) = 0.81, P > 0.05).
When we tested for LTM 24 h after CTA training (i.e., the 24 h post-test), we found that memory was present in all three cohorts. That is, in all three cohorts the number of bites elicited in the 24 h post-test was significantly smaller than that in the pretest (P < 0.01, Figure 3A). Thus, all the cohorts were capable of forming CTA-LTM following the conditioning procedure. However, we found that the feeding response in the Dutch snails in the 24 h post-test was significantly smaller than those of the Canadian TC1 and the F1 cross snails (P < 0.05, Figure 3A). That is, the Dutch snails formed better CTA-LTM than the other two cohorts of snails. The results indicated that the Dutch snails exhibited the smart phenotype, whereas the Canadian TC1 and the F1 cross snails showed average intelligence in the CTA classical conditioning trial.
In the control behavioral experiments (i.e., the backward conditioning and the naive control procedures), there were no significant differences in the feeding responses to the sucrose CS among the three cohorts at either the pretest, the 10 min post-test or the 24 h post-test (Figures 3B,C). These results support the validity of the CTA training shown in Figure 3A and highlight the different learning abilities of the Dutch and the Canadian TC1 snails at the 24 h post-test after the CTA training procedure.
Good and Poor Performers in the CTA Trial
To further confirm that the Dutch snails had better LTM than the other two cohorts, we examined the ratio between the good and poor performers in response to the CS after both the 10 min and 24 h post-test sessions in the Dutch, Canadian TC1 and F1 cross snails (Figure 4). As described in the “Materials and Methods” Section, we defined a good performer as a snail that made 0–1 bites/min during the post-test session in response to the CS. A poor performer was thus defined as a snail that made ≥2 bites/min in response to the CS during the post-test session. At the 10 min post-test after the CTA training, there were no significant differences in the learning performance among the three cohorts (Figure 4A, P > 0.05). Thus, each of the three cohorts performed equally well on the 10 min post-test. This result was in agreement with the results of the feeding responses (Figure 3A).
Figure 4. Comparison of the ratio between the good and poor performers at the memory retention tests among the Dutch, the Canadian TC1 and the F1 cross snails. Pie charts express the ratio (percentage) of the good and poor performers. (A) At 10 min post-test, no significant differences were found (P > 0.05) in the ratios of good performers (black) or poor performers (white) among the three cohorts. (B) At 24 h post-test, there was a significant difference among the three cohorts (P < 0.01). That is, the Canadian TC1 and the F1 cross snails exhibited poor formation of CTA-long-term memory (LTM).
We also compared the ratio between the good and poor performers in response to the sucrose CS at the 24 h post-test in the three cohorts (Figure 4B). The ratio of good and poor performers in the Canadian TC1 and the F1 cross snails was significantly lower than that of the Dutch snails (P < 0.01, Figure 4B). Again, these data showed that the Dutch snails were more successful at forming memory than either of the other two cohorts.
Food Deprivation as a Stressor Affecting CTA
Finally, we examined the possibility that 1-day food deprivation was a greater stressor for the Canadian TC1 snails (Hughes et al., 2017). According to the Yerkes-Dodson/Hebb inverted-U law, a strong stressor impairs learning and memory, and thus there is an optimal level of arousal for learning and memory (Yerkes and Dodson, 1908; Hebb, 1955; Ito et al., 2015). However, our results using satiated Canadian TC1 snails showed that the feeding behavior was not suppressed at the post-tests after the CTA training (Figure 5). As demonstrated previously in Dutch snails (Aonuma et al., 2016, 2017), the scores obtained in a satiated state are not as good as those obtained in a modestly hungry state (i.e., after 1 day of food deprivation). Thus, 1 day of food deprivation was not a strong stressor for either the Canadian TC1 snails or the Dutch snails.
Figure 5. CTA in the Canadian TC1 snails without 24-h food deprivation. If the food deprivation before CTA training had been stressful to the Canadian TC1 snails, we expected that the memory scores would have declined. We therefore trained the food-satiated snails by the CTA training paradigm. However, there was no significant difference among the pretest (black), the 10 min post-test (white) and the 24 h post-test (hatched). That is, the learning and memory were worse in snails in the satiated state than in snails with modest food deprivation (i.e., food deprivation of 1 day), as previously reported.
Discussion
We here showed that there were significant differences among three distinct populations of Lymnaea stagnalis (i.e., the Dutch, the Canadian TC1, and the F1 cross snails) with respect to learning and memory-formation abilities following either an operant or a classical conditioning procedure. The Dutch and Canadian TC1 strains exhibited both strengths and weaknesses in their respective ability to form memory with the two different training procedures. The Canadian TC1 snails possessed better memory-forming ability in the aerial respiratory operant conditioning trial than the Dutch and the F1 cross snails, whose learning phenotype was average. On the other hand, all three cohorts exhibited associative learning following CTA conditioning. However, the LTM formation in the Dutch snails was much better than that exhibited by the Canadian snails and the F1 cross snails. Thus, whereas the Canadian TC1 snails had a superior memory-forming capability following the operant conditioning procedure, this superiority did not translate into forming the best memory on the CTA task. Similarly, the Dutch snails performed better on the CTA task but did not exhibit superior LTM following the operant conditioning procedure. Finally, snails of the F1 cross population inherited neither the smart phenotype for operant conditioning from the Canadian TC1 snails, nor the smart phenotype for classical conditioning from the Dutch snails.
To better understand the present findings, we should consider two points: genetic differences between the Dutch and the Canadian TC1 strains; and genetic drift caused by these differences. The genetic differences were the basis for the different memory phenotypes shown above. With respect to the genotype of this species, Mozley (1939) described three sub-species of Lymnaea stagnalis in North America based on shell morphology (Lymnaea stagnalis sanctaemariae, Lymnaea stagnalis lillianae and Lymnaea stagnalis wasatchensis). He believed that these three sub-species appeared to be “specially adapted to the conditions under which they lived”. In more recent studies using newer techniques, Remigio concluded that there are genetically distinct populations within Lymnaea stagnalis (Remigio and Blair, 1997; Remigio, 2002). This conclusion was also supported by data showing a genetic variation in the tolerance of Lymnaea to copper (Côte et al., 2015). In that study, five naturally occurring freshly collected strains from specific ponds as well as three laboratory inbred strains were used. One of those inbreed strains was the Dutch strain that we used here. Further support for the genetically distinct comes from the work performed in Finland on Lymnaea stagnalis collected from a number of small lakes, ponds and the bays of larger lakes (Puurtinen et al., 2004a,b, 2007). The authors found detectable genetic variability between those populations of Lymnaea stagnalis. Genetic drift also occurred in Lymnaea stagnalis populations in outdoor closed mesocosms (Coutellec and Caquet, 2011).
We are uncertain why the F1 cross snails performed more poorly than the Canadian TC1 and the Dutch strains in both the aerial respiratory operant conditioning and the CTA classical conditioning trials. There might be genetic variable in F1 cross snails, because we obtained the F1 cross snails from the combination of different Dutch snails and different Canadian TC1 snails. That is, we purposely obtained the F1 cross snails with great genetic variable. Nevertheless, the F1 cross snails showed the average ability for learning. In the case of the CTA, one possibility is that there are different levels of dopamine or octopamine in the central nervous system (CNS) of the three populations. We posit this idea because we have previously shown that the memory performance in the Dutch strain following CTA training was negatively correlated with dopamine/octopamine levels in the CNS (Aonuma et al., 2016, 2017). It may therefore be that the F1 cross snails have higher levels of CNS dopamine/octopamine compared to the Dutch strains. Further experiments will be necessary to confirm this hypothesis and possibly to elucidate other reasons for the poor performance of this population.
Understanding why different strains of the same species, or why different individuals of the same strain, exhibit different cognitive properties is one of the most important questions to be answered in the field of neuroscience. Such differences have long been noted with inbred mouse strains (Wahlsten et al., 2003; Tipps et al., 2014). In mice, the usual procedure for establishing inbred strains is mating of brother-sister pairs for a minimum of 20 generations, resulting in lines that are roughly 99% genetically identical (Green, 1981). Many attempts are going on to find out some differences in genetics among inbred mouse strains (Graybeal et al., 2014). However, because of the complexity of the mammalian brain compared to the snail CNS and the more complicated behaviors being studied in rodents compared to the simpler behavioral studies in snails, it may be that the answers to these questions will initially be made in the snail. The interstrain difference is rare in learning abilities in invertebrates as shown in the present study, with the exception of the studies using a parasitoid wasp (Froissart et al., 2017) and mutants in Drosophila (Tully, 1996).
As noted above, genetic drift has been shown to occur in laboratory-inbreed strains of Lymnaea stagnalis (Côte et al., 2015). It is possible that over time the laboratory-bred Canadian TC1 snails will become behaviorally more like the Dutch snails in the aerial respiratory operant conditioning as well as in CTA. We will, however, be able to monitor this both behaviorally and at the single neuron level, as Braun et al. (2012) have previously demonstrated the differences in the activity of a single neuron, RPeD1, between smart and average freshly collected snails. Thus, over time a smart phenotype snail may disappear in laboratory-bred populations for the aerial respiratory operant conditioning.
We consider the results shown for the backward conditioning snails and the naive snails for CTA training (Figures 3B,C). Both the backward conditioning and the naive snails showed that there were significant differences in the feeding responses for time course (P < 0.01). That is, siginificant differences were observed between the pretest, the 10 min post-test and the 24 h post-test. This decrease in the feeding response was considered to be due not to the learning effect but rather to the fatigue, because both the backward conditioning snails and the naive snails showed a tendency toward a slight decrease in the feeding response.
The studies using parasitoid wasps by (Hoedjes et al., 2011) showed variation in learning rate and memory dynamics, and described “learning has several ecological costs”. In our present study, however, we have used the same rearing systems and performed the same experiments under the same environmental conditions for the different strains of Lymnaea. We should also note that the elapsed time for rearing in our laboratories was so long (at least 7 years). That is, we should consider genetic shifts rather than ecological factors to well understand our present results. Further, the innate preference in different strains of Lymnaea for both “air” for the aerial respiratory operant conditioning and “sucrose” for the CTA classical conditioning was not observed. Thus, we do not have the ability to determine if there are ecological costs to the different strains.
In conclusion, the different strains of Lymnaea stagnalis have different learning and memory-forming capabilities. Strains that exhibit intelligence on one test will not necessarily do so on another test. The F1 cross snails may inherit a worse ability for both aerial respiratory operant conditioning and CTA classical conditioning. The difference in learning and memory-forming ability among the strains of this species may henceforth attract the attention of researchers in learning and memory.
Author Contributions
HS conducted the CTA experiment and prepared the figures. YT conducted the aerial respiratory operant conditioning and the CTA experiments, analyzed the data and prepared the figures. RN conducted the CTA experiment. MS planned, designed the experiments and completed the manuscript. KL planned, designed the experiments, conducted the aerial respiratory operant conditioning and completed the manuscript. EI planned, designed the experiments, conducted the aerial respiratory operant conditioning, analyzed the data and held the main responsibility for the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by KAKENHI grants from the Japanese Society for the Promotion of Science (24657055 and 25291074 to EI), Waseda University grants for Specific Research Projects in Japan (2016B-068, 2016B-069 and 2016S-037 to EI), and a grant from the Natural Sciences and Engineering Research Council of Canada (227993-2013 to KL).
References
Aonuma, H., Kaneda, M., Hatakeyama, D., Watanabe, T., Lukowiak, K., and Ito, E. (2016). Relationship between the grades of a learned aversive-feeding response and the dopamine contents in Lymnaea. Biol. Open 5, 1869–1873. doi: 10.1242/bio.021634
Aonuma, H., Kaneda, M., Hatakeyama, D., Watanabe, T., Lukowiak, K., and Ito, E. (2017). Weak involvement of octopamine in aversive taste learning in a snail. Neurobiol. Learn. Mem. 141, 189–198. doi: 10.1016/j.nlm.2017.04.010
Benjamin, P. R., Staras, K., and Kemenes, G. (2000). A systems approach to the cellular analysis of associative learning in the pond snail Lymnaea. Learn. Mem. 7, 124–131. doi: 10.1101/lm.7.3.124
Bouétard, A., Côte, J., Besnard, A.-L., Collinet, M., and Coutellec, M.-A. (2014). Environmental versus anthropogenic effects on population adaptive divergence in the freshwater snail Lymnaea stagnalis. PLoS One 9:e106670. doi: 10.1371/journal.pone.0106670
Braun, M. H., Lukowiak, K., Karnik, V., and Lukowiak, K. (2012). Differences in neuronal activity explain differences in memory forming abilities of different populations of Lymnaea stagnalis. Neurobiol. Learn. Mem. 97, 173–182. doi: 10.1016/j.nlm.2011.11.005
Côte, J., Bouétard, A., Pronost, Y., Besnard, A. L., Coke, M., Piquet, F., et al. (2015). Genetic variation of Lymnaea stagnalis tolerance to copper: a test of selection hypotheses and its relevance for ecological risk assessment. Environ. Pollut. 205, 209–217. doi: 10.1016/j.envpol.2015.05.040
Coutellec, M. A., and Caquet, T. (2011). Heterosis and inbreeding depression in bottlenecked populations: a test in the hermaphroditic freshwater snail Lymnaea stagnalis. J. Evol. Biol. 24, 2248–2257. doi: 10.1111/j.1420-9101.2011.02355.x
Dalesman, S., and Lukowiak, K. (2012). How stress alters memory in ‘smart’ snails. PLoS One 7:e32334. doi: 10.1371/journal.pone.0032334
Froissart, L., Giurfa, M., Sauzet, S., and Desouhant, E. (2017). Cognitive adaptation in asexual and sexual wasps living in contrasted environments. PLoS One 12:e0177581. doi: 10.1371/journal.pone.0177581
Graybeal, C., Bachu, M., Mozhui, K., Saksida, L. M., Bussey, T. J., Sagalyn, E., et al. (2014). Strains and stressors: an analysis of touchscreen learning in genetically diverse mouse strains. PLoS One 9:e87745. doi: 10.1371/journal.pone.0087745
Green, E. L. (1981). Genetics and Probability in Animal Breeding Experiments. New York, NY: Oxford University Press.
Hatakeyama, D., Sadamoto, H., Watanabe, T., Wagatsuma, A., Kobayashi, S., Fujito, Y., et al. (2006). Requirement of new protein synthesis of a transcription factor for memory consolidation: paradoxical changes in mRNA and protein levels of C/EBP. J. Molec. Biol. 356, 569–577. doi: 10.1016/j.jmb.2005.12.009
Hebb, D. O. (1955). Drives and the C. N. S. (conceptual nervous system). Psychol. Rev. 62, 243–254. doi: 10.1037/h0041823
Hoedjes, K. M., Kruidhof, H. M., Huigens, M. E., Dicke, M., Vet, L. E. M., and Smid, H. M. (2011). Natural variation in learning rate and memory dynamics in parasitoid wasps: opportunities for converging ecology and neuroscience. Proc. Biol. Sci. 278, 889–897. doi: 10.1098/rspb.2010.2199
Hughes, E., Shymansky, T., Sunada, H., and Lukowiak, K. (2016). Qualitatively different memory states in Lymnaea as shown by differential responses to propranolol. Neurobiol. Learn. Mem. 136, 63–73. doi: 10.1016/j.nlm.2016.09.013
Hughes, E., Shymansky, T., Swinton, E., Lukowiak, K. S., Swinton, C., Sunada, H., et al. (2017). Strain-specific differences of the effects of stress on memory in Lymnaea. J. Exp. Biol. 220, 891–899. doi: 10.1242/jeb.149161
Ito, E., Kobayashi, S., Kojima, S., Sadamoto, H., and Hatakeyama, D. (1999). Associative learning in the pond snail, Lymnaea stagnalis. Zool. Sci. 6, 711–723. doi: 10.2108/zsj.16.711
Ito, E., Kojima, S., Lukowiak, K., and Sakakibara, M. (2013). From likes to dislikes: conditioned taste aversion in the great pond snail (Lymnaea stagnalis). Can. J. Zool. 91, 405–412. doi: 10.1139/cjz-2012-0292
Ito, E., Yamagishi, M., Takigami, S., Sakakibara, M., Fujito, Y., and Lukowiak, K. (2015). The Yerkes-Dodson law and appropriate stimuli for conditioned taste aversion in Lymnaea. J. Exp. Biol. 218, 336–339. doi: 10.1242/jeb.113266
Kemenes, G., and Benjamin, P. R. (1989). Appetitive learning in snails shows characteristics of conditioning in vertebrates. Brain Res. 489, 163–166. doi: 10.1016/0006-8993(89)90019-x
Kojima, S., Sunada, H., Mita, K., Sakakibara, M., Lukowiak, K., and Ito, E. (2015). Function of insulin in snail brain in associative learning. J. Comp. Physiol. A Neuroethol. Sens. Neural Behav. Physiol. 201, 969–981. doi: 10.1007/s00359-015-1032-5
Kojima, S., Yamanaka, M., Fujito, Y., and Ito, E. (1996). Differential neuroethological effects of aversive and appetitive reinforcing stimuli on associative learning in Lymnaea stagnalis. Zool. Sci. 13, 803–812. doi: 10.2108/zsj.13.803
Lukowiak, K., Adatia, N., Krygier, D., and Syed, N. (2000). Operant conditioning in Lymnaea: evidence for intermediate- and long-term memory. Learn. Mem. 7, 140–150. doi: 10.1101/lm.7.3.140
Lukowiak, K., Ringseis, E., Spencer, G., Wildering, W., and Syed, N. (1996). Operant conditioning of aerial respiratory behaviour in Lymnaea stagnalis. J. Exp. Biol. 199, 683–691.
Lukowiak, K., and Syed, N. (1999). Learning, memory and a respiratory central pattern generator. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 124, 265–274. doi: 10.1016/s1095-6433(99)00114-2
Mita, K., Okuta, A., Okada, R., Hatakeyama, D., Otsuka, E., Yamagishi, M., et al. (2014). What are the elements of motivation for acquisition of conditioned taste aversion? Neurobiol. Learn. Mem. 107, 1–12. doi: 10.1016/j.nlm.2013.10.013
Mozley, A. (1939). The variation of Lymnaea (Linné). J. Mollusc. Stud. 23, 267–269. doi: 10.1093/oxfordjournals.mollus.a064365
Murakami, J., Okada, R., Sadamoto, H., Kobayashi, S., Mita, K., Sakamoto, Y., et al. (2013). Involvement of insulin-like peptide in long-term synaptic plasticity and long-term memory of the pond snail Lymnaea stagnalis. J. Neurosci. 33, 371–383. doi: 10.1523/JNEUROSCI.0679-12.2013
Orr, M. V., Hittel, K., and Lukowiak, K. (2008). Comparing memory-forming capabilities between laboratory-reared and wild Lymnaea: learning in the wild, a heritable component of snail memory. J. Exp. Biol. 211, 2807–2816. doi: 10.1242/jeb.020172
Orr, M., Hittel, K., Lukowiak, K. S., Han, J., and Lukowiak, K. (2009). Differences in LTM-forming capability between geographically different strains of Alberta Lymnaea stagnalis are maintained whether they are trained in the lab or in the wild. J. Exp. Biol. 212, 3911–3918. doi: 10.1242/jeb.024281
Otsuka, E., Matsunaga, M., Okada, R., Yamagishi, M., Okuta, A., Lukowiak, K., et al. (2013). Increase in cyclic AMP concentration in a cerebral giant interneuron mimics part of a memory trace for conditioned taste aversion of the pond snail. Biophysics 9, 161–166. doi: 10.2142/biophysics.9.161
Puurtinen, M., Emily Knott, K., Suonpää, S., Nissinen, K., and Kaitala, V. (2007). Predominance of outcrossing in Lymnaea stagnalis despite low apparent fitness costs of self-fertilization. J. Evol. Biol. 20, 901–912. doi: 10.1111/j.1420-9101.2007.01312.x
Puurtinen, M., Hytönen, M., Knott, K. E., Taskinen, J., Nissinen, K., and Kaitalia, V. (2004a). The effects of mating system and genetic variability on susceptibility to trematode parasites in a freshwater snail, Lymnaea stagnalis. Evolution 58, 2747–2753. doi: 10.1554/04-465
Puurtinen, M., Emily Knott, K., Suonpää, S., van Ooik, T., and Kaitalia, V. (2004b). Genetic variability and drift load in populations of an aquatic snail. Evolution 58, 749–756. doi: 10.1554/03-434
Remigio, E. A. (2002). Molecular phylogenetic relationships in the aquatic snail genus Lymnaea, the intermediate host of the causative agent of fascioliasis: insights from broader taxon sampling. Parasitol. Res. 88, 687–696. doi: 10.1007/s00436-002-0658-8
Remigio, E. A., and Blair, D. (1997). Molecular systematics of the freshwater snail family Lymnaeidae (Pulmonata: Basommatophora) utilising mitochondrial ribosomal DNA sequences. J. Mollusc. Stud. 63, 173–185. doi: 10.1093/mollus/63.2.173
Sadamoto, H., Yamanaka, M., Hatakeyama, D., Nakamura, H., Kojima, S., Yamashita, M., et al. (2000). Developmental study of anatomical substrate for conditioned taste aversion in Lymnaea stagnalis. Zool. Sci. 17, 141–148. doi: 10.2108/zsj.17.141
Shymansky, T., Protheroe, A., Hughes, E., Swinton, C., Swinton, E., Lukowiak, K. S., et al. (2017). Juveniles of Lymnaea ‘smart’ snails do not perseverate and have the capacity to form LTM. J. Exp. Biol. 220, 408–413. doi: 10.1242/jeb.147884
Sugai, R., Azami, S., Shiga, H., Watanabe, T., Sadamoto, H., Kobayashi, S., et al. (2007). One-trial conditioned taste aversion in Lymnaea: good and poor performers in long-term memory acquisition. J. Exp. Biol. 210, 1225–1237. doi: 10.1242/jeb.02735
Sunada, H., Lukowiak, K., and Ito, E. (2017a). Cerebral giant cells are necessary for the formation and recall of memory of conditioned taste aversion in Lymnaea. Zool. Sci. 34, 72–80. doi: 10.2108/zs160152
Sunada, H., Watanabe, T., Hatakeyama, D., Lee, S., Forest, J., Sakakibara, M., et al. (2017b). Pharmacological effects of cannabinoids on learning and memory in Lymnaea. J. Exp. Biol. 220, 3026–3038. doi: 10.1242/jeb.159038
Takigami, S., Sunada, H., Lukowiak, K., Ito, E., and Sakakibara, M. (2016). An automated learning apparatus for classical conditioning of Lymnaea stagnalis. J. Neurosci. Methods 259, 115–121. doi: 10.1016/j.jneumeth.2015.10.008
Tipps, M. E., Raybuck, J. D., Buck, K. J., and Lattal, K. M. (2014). Delay and trace fear conditioning in C57BL/6 and DBA/2 mice: issues of measurement and performance. Learn. Mem. 21, 380–393. doi: 10.1101/lm.035261.114
Tully, T. (1996). Discovery of genes involved with learning and memory: an experimental synthesis of Hirschian and Benzerian perspectives. Proc. Natl. Acad. Sci. U S A 93, 13460–13467. doi: 10.1073/pnas.93.24.13460
Van Der Steen, W. J., van der Hoven, N. P., and Jager, J. C. (1969). A method for breeding and studying freshwater snails under continuous water change, with some remarks on growth and reproduction in Lymnaea stagnalis (L.). Netherlands J. Zool. 19, 131–139. doi: 10.1163/002829669x00170
Wagatsuma, A., Sugai, R., Chono, K., Azami, S., Hatakeyama, D., Sadamoto, H., et al. (2004). The early snail acquires the learning. Comparison of scores for conditioned taste aversion between morning and afternoon. Acta Biol. Hung. 55, 149–155. doi: 10.1556/abiol.55.2004.1-4.18
Wahlsten, D., Metten, P., Phillips, T. J., Boehm, S. L. II, Burkhart-Kasch, S., Dorow, J., et al. (2003). Different data from different labs: lessons from studies of gene-environment interaction. J. Neurobiol. 54, 283–311. doi: 10.1002/neu.10173
Yamagishi, M., Watanabe, T., Hatakeyama, D., and Ito, E. (2015). Effects of serotonin on the heartbeat of pond snails in a hunger state. Biophysics 11, 1–5. doi: 10.2142/biophysics.11.1
Yamanaka, M., Sadamoto, H., Hatakeyama, D., Nakamura, H., Kojima, S., Kimura, T., et al. (1999). Developmental changes in conditioned taste aversion in Lymnaea stagnalis. Zool. Sci. 16, 9–16. doi: 10.2108/zsj.16.9
Keywords: aerial respiratory operant conditioning, conditioned taste aversion, F1 cross, Lymnaea, interstrain differences
Citation: Sunada H, Totani Y, Nakamura R, Sakakibara M, Lukowiak K and Ito E (2017) Two Strains of Lymnaea stagnalis and the Progeny from Their Mating Display Differential Memory-Forming Ability on Associative Learning Tasks. Front. Behav. Neurosci. 11:161. doi: 10.3389/fnbeh.2017.00161
Received: 10 May 2017; Accepted: 16 August 2017;
Published: 11 September 2017.
Edited by:
Martin Giurfa, UMR5169 Centre de Recherches sur la Cognition Animale (CRCA), FranceReviewed by:
Zoltán Serfozo, Balaton Limnological Institute (MTA), HungaryEmmanuel Desouhant, University of Lyon, France
Copyright © 2017 Sunada, Totani, Nakamura, Sakakibara, Lukowiak and Ito. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Etsuro Ito, ZWl0b0B3YXNlZGEuanA=
† These authors have contributed equally to this work.