Skip to main content

ORIGINAL RESEARCH article

Front. Neurosci., 20 September 2022
Sec. Perception Science
This article is part of the Research Topic Dynamic Cognition: Empirical and Theoretical Insights From Nonhuman Primate Studies in Naturalistic Conditions View all 5 articles

The effect of ketamine on eye movement characteristics during free-viewing of natural images in common marmosets

\r\nZlata PolyakovaZlata Polyakova1Masao IwaseMasao Iwase2Ryota HashimotoRyota Hashimoto3Masatoshi Yoshida*Masatoshi Yoshida1*
  • 1Center for Human Nature, Artificial Intelligence, and Neuroscience, Hokkaido University, Sapporo, Japan
  • 2Department of Psychiatry, Osaka University Graduate School of Medicine, Suita, Japan
  • 3Department of Pathology of Mental Diseases, National Institute of Mental Health, National Center of Neurology and Psychiatry, Kodaira, Japan

Various eye movement abnormalities and impairments in visual information processing have been reported in patients with schizophrenia. Therefore, dysfunction of saccadic eye movements is a potential biological marker for schizophrenia. In the present study, we used a pharmacological model of schizophrenia symptoms in marmosets and compared the eye movement characteristics of marmosets during free-viewing, using an image set identical to those used for human studies. It contains natural and complex images that were randomly presented for 8 s. As a pharmacological model of schizophrenia symptoms, a subanesthetic dose of ketamine was injected intramuscularly for transient and reversible manipulation. Eye movements were recorded and compared under a ketamine condition and a saline condition as a control. The results showed that ketamine affected eye movement characteristics during free-viewing. Saccades amplitude and scanpath length were significantly reduced in the ketamine condition. In addition, the duration of saccades was longer under the ketamine condition than under the saline condition. A similar tendency was observed for the duration of fixations. The number of saccades and fixations tended to decrease in the ketamine condition. The peak saccades velocity also decreased after ketamine injection whereas there was no difference in the main sequence relationship between saccades amplitude and peak velocity. These results suggest that ketamine affected visual exploration but did not affect the oculomotor aspect of saccades in marmosets, consistent with studies in patients with schizophrenia. Therefore, we conclude that the subanesthetic dose of ketamine is a promising pharmacological model of schizophrenia symptoms in common marmosets and can be used in combination with free-viewing paradigms to establish “translatable markers” for schizophrenia in primates.

Introduction

Schizophrenia is known to affect one percent of the population and is classified as a mental disorder with psychosis (Andreasen, 1995). The brain mechanism of this disorder remains largely unknown and there is a need to identify behavioral markers that can be used in pathophysiology at the neuronal level using animal models. Human studies have reported dysfunction of eye movements in subjects with schizophrenia in various tasks (Holzman et al., 1973; Kojima et al., 1992; Levy et al., 1993; Arolt et al., 1998; Campana et al., 1999). Impaired visual information processing and abnormal eye movements during free-viewing of naturalistic images have been shown to be possible neurophysiological biomarkers of schizophrenia (Bestelmeyer et al., 2006; Benson et al., 2007; Miura et al., 2014; Yoshida et al., 2015). Common features of schizophrenic subjects associated with atypical visual exploration of images have been reported: prolonged fixation duration, decreased saccade amplitude, decreased saccade frequency, and decreased scanpath length (Morita et al., 2017). Increased salience-guided eye movements were reported in subjects with schizophrenia (Yoshida et al., 2015). Therefore, we hypothesize that the visuo-cognitive impairments during free-viewing reported in patients with schizophrenia may be used as an indicator in animal model of schizophrenia symptoms.

A promising non-human primate model for investigating cognitive dysfunction in neuroscientific studies is a common marmoset (Callithrix jacchus) model (Okano et al., 2012). The possibility of genetic manipulations and brain structure (e.g., lissencephalic cortex) of this species has multiple advantages for electrophysiological and optical imaging neuroscientific studies (Hung et al., 2015; Sadakane et al., 2015; Yamada et al., 2016; Ebina et al., 2018; Kondo et al., 2018). Marmosets, like macaques and humans, have the potential to be widely adopted for oculomotor research due to their similar saccade kinematics and their ability to perform psychophysical tasks (Mitchell et al., 2014; Chen et al., 2021). Thus, establishing an experimental paradigm for studying schizophrenia using common marmosets has multiple scientific benefits.

The motivation for the choice of free-viewing as a behavioral task is based on the fact that a significant amount of human studies on subjects with schizophrenia reported changes in eye movement characteristics during the observation of images. In addition, this type of task does not require pre-training or high cognitive capacity. Visual behavior of marmosets is comparable to that of other primates (Mitchell and Leopold, 2018) and eye movement monitoring experiments under head restraint are already well developed (Mitchell et al., 2014, 2015; Chen et al., 2021; Kaneko et al., 2022). Therefore, in the current study, we used a free-viewing task with natural and complex images presented to head-restrained marmosets. To date, there are no studies that directly compare eye movement characteristics of marmosets as a pharmacological model of schizophrenia symptoms with those of patients with schizophrenia, because of differences in visual stimuli. In the present study, we addressed this issue by choosing the same natural and complex images used in the human studies (Miura et al., 2014; Morita et al., 2017, 2019).

There are multiple approaches to establish animal models of schizophrenia (Carpenter and Koenig, 2008), which can be classified into four categories: developmental, pharmacological, lesion, or genetic manipulation (Jones et al., 2011). In the current study, we prioritized the presence of a reversible effect, which is possible in drug-induced models. Human data have shown that N-Methyl-D-aspartate (NMDA) antagonists such as ketamine can produce clinical symptoms of schizophrenia (Krystal et al., 1994). Pharmacological models with ketamine have been reported to produce both positive and negative symptoms of schizophrenia (Beck et al., 2020). Early human studies showed that high doses of ketamine strongly affected characteristics of eye movements (Radant et al., 1998), and it was subsequently shown that even low doses reproduce some of the abnormalities seen in subjects with schizophrenia, such as reduced eye acceleration during initiation, reduced closed loop gain (Weiler et al., 2000). Injections of subanesthetic dose of ketamine 0.3–0.5 mg/kg in non-human primates reproduce various sensory and cognitive dysfunctions of schizophrenia, such as dysfunction of the visual working memory (Sawagashira and Tanaka, 2021) and impaired visual contextual integration (Schielke and Krekelberg, 2021) in macaques, deficits in executive function in marmosets (Kotani et al., 2016). Recently, it has been reported in marmosets that ketamine also induced disruption of scan paths, with limited effects on saccade motor control during free-viewing of face images (Selvanayagam et al., 2021). However, there is no study to date examining the effect of ketamine during free-viewing of natural and/or complex images with various categories. To fill this gap, the present study examined whether subanesthetic doses of ketamine reproduce schizophrenia-like oculomotor symptoms in common marmosets, and directly compared changes in eye movement characteristics with data from human studies conducted using the same task and image dataset (Miura et al., 2014; Yoshida et al., 2015).

Materials and methods

Animals

The experimental protocols were approved by the Institutional Animal Care and Use Committee of the National Institutes of Natural Sciences and Hokkaido University Animal Care and Use Committee, and all experiments were performed in accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals. Five marmosets (Callithrix jacchus) were subjects in this study. Marmoset U (8 years, 360 g), Marmoset E (3 years, 390 g), Marmoset G (5 years, 370 g), and Marmoset J (4 years, 410 g) were males. Marmoset F was a 9 years old female and weighed 460 g. When animals were not in active use, they were housed in individual cages and had ad libitum access to food and water. During the experiments, the animals were under the close supervision of university veterinarians and/or experimenters.

Marmosets’ preparation

First, each marmoset was trained to sit in a marmoset chair (MMC-3505, O’Hara & Co., Tokyo, Japan) quietly for 15 min. During this first stage of chair training experimenter periodically rewarded the animal with a small piece of Baumkuchen cake. After animals were conditioned to be handled by the experimenter and to the movement restrictions, each marmoset underwent an aseptic surgical operation under isoflurane anesthesia. The head-holder implant (Reinforced plastic) was mounted using Estecem II (Tokuyama Dental) and used to stabilize the head during eye-tracking data collection. Then, several sessions of chair training under head-fixed conditions were repeated to give time for animals to adapt to a new restrained condition and stay calm for 30 min while licking continuously administered sweetened water. After these procedures, a free-viewing task was conducted.

Behavioral task

Free-viewing stimuli

In the present study, we used the same test pictures data set as for human studies (Yoshida et al., 2015; Morita et al., 2017), which allows us to perform a direct comparison of eye movement characteristics reported for human subjects with schizophrenia and data collected on marmoset with the pharmacological model of the disorder. The dataset of 56 natural and complex images in seven categories (nature, animal, face, building, daily, food, geometric, noise) was used in the study. Face images were selected from Matsumoto and Ekman database (Matsumoto and Ekman, 1988) and the rest from the International Affective Pictures System (IAPS) (Lang et al., 2008), which are not allowed to be published in scientific journals. Test images were randomly presented for 8 s, interleaved with calibration pictures with marmoset faces in different positions presented for 4 s (Figure 1A). The gaze calibration was performed using images with four small pictures of marmoset faces located in the corners or smaller versions of the same pictures grouped in the corner/center (Figure 1B). Task conditions were similar to the one used in the human studies (Miura et al., 2014; Morita et al., 2017) with slight minor modifications (e.g., calibration procedure) for studying non-human primates.

FIGURE 1
www.frontiersin.org

Figure 1. Methods. (A) Experimental design. Eye tracking in the free-viewing task was implemented in head-restrained marmosets (n = 5). Injection of saline or ketamine to marmosets was followed by a short time gap (approximately 4 min) to avoid eye nystagmus in the case of drug application. The data collection was performed in four time periods. Each period (approximately 5 min) includes a set of 22 calibration and 22 test images presented in the queue order. (B) Examples of gaze (blue traces) in Marmoset G during free-viewing of images. (1) Two representative examples of calibration images (see section “Materials and methods”) used for the eye position offset and gain adjustment. (2) Example of eye tracking data during free-viewing of the same test image from the category of faces after saline (left) and ketamine (right) injections. The test images are replaced with a gray rectangle. See section “Materials and methods” for detail.

During data collection, marmosets sat in the marmoset chair with its head fixed using the head-holder. The LCD monitor (Dell P1917S, 375 mm × 300 mm, 60 Hz) was placed at a distance of 40cm in front of the marmosets. The resolution of the screen was 1280 × 1024 pixels. Images were presented in the full-screen mode. The marmosets continuously received drops of sweetened water in order to maintain an arousal state and attention.

Eye tracking

Eye movements during free-viewing were binocularly recorded at 500 Hz with the EyeLink 1000 Plus eye-tracker (SR Research, Mississauga, ON, Canada). The eye position was detected based on pupil tracking and corneal reflection information available in the Eyelink system. The blink segments detected automatically by Eyelink were recorded and used for further exclusion of the data with artifacts.

Injections

We injected a subanesthetic dose of ketamine 0.5 mg/kg (Daiichi-Sankyo Co., Tokyo, Japan) intramuscularly for transient and reversible manipulation of schizophrenia symptoms. The dosage of ketamine was determined based on previous pharmacological studies in monkeys (Pouget et al., 2010; Sawagashira and Tanaka, 2021) and marmosets (Selvanayagam et al., 2021), so that marmosets could perform the free-viewing task. To avoid the development of the drug tolerance and cumulative effect, injections of ketamine were performed once a week. The saline injections were performed on different days in between ketamine injections and were used as a control condition for further comparison. The recording of eye movements was started after the strong effect of nystagmus was diminished and covered approximately 25 min time interval after the injection (Figure 1A).

Data analysis

The aim was to evaluate the oculomotor behavior of marmosets during free-viewing in a pharmacological model of schizophrenia symptoms. We performed a detailed analysis of pupil changes, gaze behavior, and saccade kinematics in control and modeled schizophrenic states.

Pupil changes

Shortly after the ketamine injection, we observed an increase in pupil size that lasts approximately 10–15 min. Since the pupil effect did not affect the eye detection quality, data reordered during pupil dilation were included in the analysis and used as an indicator of a successful drug administration.

Eye movement detection

EyeLink recorded gaze position and pupil diameter, which were used for the detailed analysis of eye movement characteristics. First, preprocessing of recorded data was done using MATLAB code written in-house. It includes a selection of time intervals without blinks and strong nystagmus caused by ketamine injection. Then, data were examined using GUI as described elsewhere (Chen and Hafed, 2013). Briefly, the eye movement velocities and accelerations were calculated using a fourth-order 21-Hz differentiating Savitzky-Golay smoothing FIR (Yoshida et al., 2012). Then we identified time samples of the movement in which radial eye velocity exceeded a threshold of 15°/s to exclude the noise. The start and end of saccades were defined based on the acceleration threshold of 450°/s2 (Chen et al., 2021; Selvanayagam et al., 2021). In addition, the minimal duration of saccades was set 15 ms and the minimal intersaccadic interval 10 ms. Finally, we manually verified all detected saccades. In addition, intervals of noise artifacts were marked by the experimenter and removed from the data. After all procedures mentioned above, the fixations were identified using self-written code as intervals of data without artifacts, blinks, and saccades with a threshold of minimal duration 50 ms.

Calculation of eye movement characteristics

In the present study, we evaluated parameters of saccades kinematics in the control saline condition and modeled schizophrenic state. The number of fixations/saccades was evaluated as a number of events within test image presentation time (8 s) excluding the time of blinks and noise. Duration of fixations/saccades similarly was normalized by eight seconds time interval of active viewing. Saccades amplitude in visual degrees and peak velocity (degrees per second) were calculated for each saccade and used further for the investigation of the main sequence relationship. Scanpath length was evaluated for each test image as the summation of amplitudes of saccades detected within the image presentation. The described parameters were averaged for each test image presentation case and then the median value with median absolute deviation was calculated among all images presented for each animal under the specific condition (saline, ketamine).

Main sequence

To characterize the main sequence, we fitted the amplitude versus the peak velocity of the saccades collected during free-viewing with the following equation (Leigh and Zee, 2015): V = Vmax × {1 – exp. (–A/C)}, where (A) is the amplitude and (V) the peak eye velocity of the saccades. Parameters of the main sequence Vmax and C were optimized.

Statistics

All statistical analyses were performed using MATLAB (Mathworks, Natick, MA, USA). Comparisons between median values of eye movement characteristics with saline and ketamine injections were performed using the two-sample Wilcoxon rank-sum test. The effect size of the estimated differences was quantified by Cliff’s delta (Δ). The general guidelines for interpretation was the following: negligible for |Δ| < 0.147, small for 0.147 ≤ |Δ| < 0.33, medium for 0.33 < = |Δ| ≤ 0.474, and large for |Δ| ≥ = 0.474 (Romano et al., 2006). A false discovery rate (FDR) with α = 0.05 was used for the post hoc correction (Benjamini and Hochberg, 1995).

For a statistical test of the main sequence relationship, a bootstrap resampling method was used. To test the null hypothesis in which the Vmax for ketamine condition and saline condition are the same, the data set was resampled from merged population and fitting was done to obtain Vmax for the resampled ketamine population and the resampled saline population. By repeating resampling procedures 999 times, the distribution of the difference between Vmax for the resampled ketamine population and the resampled saline population was obtained. Then p value was calculated as the probability of resampled differences higher (or lower) than that of actual difference. To make the p values for a two-sided test, the p values were multiplied by two.

Differences among parameters for all five marmosets were evaluated using two repeated measures analyses of variance (ANOVA) with factors: marmoset and injection type (saline, ketamine). To evaluate the overall effect of injection type across all animals, a model without interaction term was used. Changes related to the time period were evaluated using three-way ANOVA with marmoset, injection type, and time period factors. In the case of image category dependency, we performed a three-way ANOVA with factors: marmoset, injection type, and image category.

Results

Overview of recorded data

In the present study, we recorded free-viewing data during the presentation of natural and/or complex images from five common marmosets (Figure 1A). The number of recording sessions per injection condition is indicated in Table 1. Each recording session included four time periods (one period ∼ 5 min; Figure 1A), however in some of them the eyes detection was not reliable due to marmosets’ behavior or the effect of injection, thus was excluded from the analysis (Table 1). Moreover, the ketamine injection induced rapid, involuntary eye movements (nystagmus) within 2–4 min which was excluded from the data analysis. The total number of images with active viewing for five animals was 1006 in the control saline case and 1046 in the ketamine pharmacological model condition. In the analysis of eye movement characteristics described below, the values were averaged for each image. Saccades were reliably detected in all animals before and after ketamine injection. The number of saccades is summarized in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Database overview.

Example from one animal

The effect of ketamine on the distribution of saccades amplitudes is shown for one animal – Marmoset G (Figure 2, left). The example data demonstrates the decrease in saccades amplitudes after ketamine injection, in comparison to saline injections. On the other hand, the main sequence relationship, which is represented by plotting the peak velocity and the saccade amplitude of each saccade, was overall unaffected by ketamine injection (Figure 2, right).

FIGURE 2
www.frontiersin.org

Figure 2. Example of the effect of ketamine (magenta) vs. saline (blue) on eye movements for Marmoset G. Saccade kinematics during free-viewing of natural images. (Left) Saccades amplitudes for two conditions plotted on a log scale. The median of the data distribution for two conditions is indicated by a shaded vertical line. (Right) Main sequence relationship. The saccade amplitude vs. the peak velocity is plotted separately for saline and ketamine cases. Dots indicate all saccades collected during the presentation of images. The curve represents the fitting of the data to the equation described in section “Materials and methods.”

Statistical analysis of eye movement characteristics

The analysis of major eye movement characteristics was performed for the saline control case and ketamine pharmacological model condition in five marmosets (Table 2 and Figure 3).

TABLE 2
www.frontiersin.org

Table 2. Eye movement characteristics in the control condition and the ketamine condition in five marmosets.

FIGURE 3
www.frontiersin.org

Figure 3. Eye movement characteristics in control saline case and ketamine pharmacological model condition plotted for five marmosets. Median values for two conditions plotted (black solid line) for the following characteristics: (A) Number of saccades; (B) Amplitudes of saccades; (C) Scanpath length. The values averaged per picture are indicated as a dot on a violin plot. Error bars denote the first and the third quartile. Significant cases after FDR post hoc correction are indicated by an asterisk.

Number of fixations/saccades

We compared the number of fixations/saccades in the saline and ketamine injection conditions. The number of fixations significantly decreased in three marmosets (Marmosets F, E, J; Table 2). Similar tendency was observed for the number of saccades (Marmosets F, U, E; Table 2 and Figure 3A) in the ketamine condition. A two-way ANOVA was conducted that examined the effect of ketamine and saline injections on the number of fixations/saccades for the data from all animals. The significant difference was found in both cases of number of fixations [F(1,2126) = 4.05, p = 4.42 × 10–2] and number of saccades [F(1,2152) = 18.07, p = 2.22 × 10–5].

Duration of fixations/saccades

Duration of fixations had a tendency to increase after ketamine injection in four animals and was significant in three (Marmosets F, U, E; Table 2). Duration of saccades was increased significantly in four marmosets out of five (Marmosets U, E, G, J; Table 2). We conducted repeated measures ANOVA and observed significant effect of ketamine over saline for both duration of fixations [F(1,2002) = 13.07, p = 3.07 × 10–4] and saccades [F(1,2046) = 15.06, p = 1.07 × 10–4].

Saccades amplitude and scanpath length

The analysis of eye movements revealed that the mean saccade amplitudes were reduced in four marmosets after injection of ketamine (Marmosets U, E, G, J; Table 2 and Figures 2, 3B). We confirmed the presence of a significant decrease in saccade amplitudes after ketamine injection using two-way ANOVA for data from five marmosets [F(1,2034) = 66.35, p = 6.51 × 10–16].

Scanpath length was reduced significantly in four animals (Marmosets F, U, E, G; Table 2 and Figure 3C). A two-way ANOVA for all marmosets confirmed the presence of significant differences in scanpath length between ketamine and saline injection conditions [F(1,2139) = 36.61, p = 1.7 × 10–9].

Peak velocity and main sequence relationship

Peak velocity significantly decreased in four animals (Marmosets U, E, G, J; Table 2) in the ketamine condition. ANOVA results for the full dataset supported this finding [F(1,2042) = 58.84, p = 2.64 × 10–14].

The slope of the main sequence relationship between the peak velocities and the amplitudes of saccades was not altered by ketamine injection (Figure 2, right). Estimated parameters of the main sequence relationship Vmax and C were similar across marmosets in the case of saline injection and not significantly different from parameters optimized for results after ketamine injection (Table 3).

TABLE 3
www.frontiersin.org

Table 3. Optimized parameters of the main sequence relationship (see section “Materials and methods”) for each marmoset.

All eye movement characteristics analyzed in the present study significantly changed after administration of ketamine to marmosets, except for main sequence relationship parameters. The effect size (Cliff’s delta) for the majority of significant changes was within a range of small and medium guideline intervals (see section “Materials and methods”). For other cases, there were no significant differences.

Duration of the ketamine effect

Our experimental data demonstrate that the major differences in eye movement characteristics were in place within three time periods after ketamine injection and the effect was reduced in the fourth period (Figure 4). Three-way ANOVA demonstrated presence of statistical significance in the interaction between injection type and time period, for all characteristics related to saccades [duration of saccades F(3,2024) = 3.14, p = 2.46 × 10–2; saccades amplitudes F(3,2008) = 12.04, p = 8.18 × 10–8; peak velocity F(3,2020) = 11.03, p = 3.5 × 10–7; scanpath length F(3,2117) = 12.3, p = 5.59 × 10–8], except for the number of saccades. This finding indicates that the effect of ketamine is altered across time periods. The example data for saccades amplitudes and scanpath length in five animals illustrates that the effect of ketamine was reduced in the fourth period (Figure 4).

FIGURE 4
www.frontiersin.org

Figure 4. Changes in the eye movement characteristics within time periods after the injection (saline vs. ketamine). Example of alterations in amplitudes of saccades (log values) and scanpath length. (Up) the difference between median values of ketamine and saline cases indicated (black circle) for each time period and animal (gray line). (Down) Cliff’s delta was used to quantify the effect size of the estimated difference between median values in two injection conditions. For interpretation, we followed the general guidelines described in section “Materials and methods.”

Effect of image category

We checked whether there is a relation between an image category (nature, animal, face, building, daily, food, geometric, noise) and eye movement characteristics in the free-viewing task. The example of saccades amplitudes in saline and ketamine conditions demonstrates higher median values for the animal image category compared to noise, as well as a decrease in saccades amplitudes after drug application (Figure 5A). Cliff’s delta (the effect size) for each animal and image category was calculated for each characteristic. There was no simple relationship found for all eye movement characteristics. The example for amplitudes of saccades and scanpath length illustrates this finding (Figure 5B). Three-way ANOVA with the interaction between injection type and image category revealed no statistical significance in eye movement characteristics. This finding suggests that the effect of ketamine is general and does not depend on specific image categories. Thus, the effect of ketamine in our study is general enough and does not depend on the choice of images. In the model without interaction, image category is statistically significant in all cases related to saccades [number of saccades F(7,2131) = 2.27, p = 2.67 × 10–2; saccades amplitudes F(7,2023) = 5.14, p = 8.14 × 10–6; peak velocity F(7,2031) = 10.11, p = 1.74 × 10–12; scanpath length F(7,2132) = 7.12, p = 1.98 × 10–8], except for saccades duration. The major eye movement characteristics are dependent on image category.

FIGURE 5
www.frontiersin.org

Figure 5. The effect of image categories on the eye movement characteristics. (A) Example of saccades amplitudes after saline or ketamine injection for ‘animal’ and ‘noise’ image categories. Median values for a log of saccades amplitudes are indicated (black circle) and connected in two conditions for each animal (gray line). Note: Monkey J has not been presented with ‘noise’ images. (B) Effect sizes (Cliff’s delta, see section “Materials and methods”) for the test of difference in amplitudes of saccades (log values) and scanpath length median values in each image category are indicated by a symbol for each marmoset: F, asterisk; U, triangle; E, circle; G, square; J, cross.

Discussion

Our main purpose was to directly compare eye movement characteristics in patients with schizophrenia and the marmoset pharmacological model of schizophrenia symptoms. The analysis of eye movement characteristics after a subanesthetic dose of ketamine (intramuscular injection) applied in marmosets revealed similar tendencies of decreased number of fixations and saccades, reduced saccade amplitude, peak velocity, and scanpath length, and increased fixation duration. In addition, the parameters of the main sequence relationship (saccade amplitude vs. saccade peak velocity) were not affected by ketamine application. These results suggest that ketamine affected visual exploration including some attentional processes, but not oculomotor activities, as argued below. These results are consistent with data reported for subjects with schizophrenia in eye tracking studies using the same image dataset (Miura et al., 2014; Morita et al., 2017, 2019) and provide evidence for the validity of the ketamine pharmacological model. The difference was observed only for saccade duration, which was increased in marmosets but decreased in humans. Overall, our results support our hypothesis that the combination of the ketamine pharmacological model in marmosets and the free-viewing task is useful to study schizophrenia.

Visual exploration vs. oculomotor function

After the ketamine administration, there was no significant difference in the main sequence relationship between saccades amplitude and peak velocity whereas there was a disturbance in scanpath length. Whereas eye movement characteristics may reflect a variety of functions such as covert attention and target selection (Schiller and Tehovnik, 2005), the main sequence relationship has been used to assess the oculomotor dysfunctions (Gibaldi and Sabatini, 2021). These previous findings lead us to the conclusion that visual exploration possibly predominates over oculomotor functions in the ketamine model of schizophrenia symptoms. Further studies examining visual salience may help to clarify how visual attention is affected in schizophrenia, as described below.

Free-viewing of naturalistic stimuli

In the present study, natural and complex images were used to test the pharmacological model of schizophrenia symptoms in marmosets. The dataset was identical to that used for the human studies (Miura et al., 2014; Morita et al., 2017, 2019). Naturalistic stimuli are beneficial for human and animal studies because they mimic “real life” visual inputs in the experiment. This approach retains higher ecological validity, in contrast to analogs using tightly controlled stimuli (e.g., moving dots, oriented lines, etc.). In addition, free-viewing of natural images is suitable for extending research toward the analysis of visual salience and processes behind this phenomenon because various visual features can be examined with naturalistic images. Thus, the combination of salience analysis and neurophysiology can be beneficial in studying visual attention and its pathophysiological changes (Yoshida et al., 2012).

Ketamine as a model of schizophrenia symptoms

In the present study, subanesthetic doses of ketamine were applied as a model of schizophrenia symptoms in marmosets. Ketamine was chosen because its effect is reversible, whereas it is irreversible in other drugs such as phencyclidine (Jentsch and Roth, 1999) or methamphetamine (Featherstone et al., 2007). This was confirmed in the present study, with the effect of ketamine decreasing in the fourth period, corresponding to 20–25 min after injection (Figure 4). The three-way ANOVA also revealed a significant interaction between drug effect and period, suggesting that the effect of ketamine varied with periods. Because the effects of ketamine are reversible, it is possible to track changes in neuronal activity before and after ketamine injection, which has a potential advantage for future neurophysiological applications.

Ketamine is also known to have antidepressant effects which last up to 7 days (Zarate et al., 2006). Since the effect of ketamine on eye movement characteristics in the present study was transient and diminished in 25 min, it is unlikely that the present findings reflect the antidepressant effects of ketamine.

The effect of ketamine on image category

As shown in Figure 5 and in the accompanying analysis, the main effect of image category was statistically significant in all cases, i.e., the eye movement characteristics depend on categories. However, the interaction between drug and image category was not statistically significant, indicating that the effect of ketamine was not affected by image category. This suggests that the effect of ketamine is robust and independent of specific image types. This finding provides flexibility in the selection of image databases for future experiments and could be important based on the literature on gaze metrics in the schizophrenia research (Hashimoto, 2021). We also note that the effect of ketamine can depend on visual salience of images and further experiments are necessary to clarify this possibility.

Future directions

An extension of this study is to examine visual salience at gazes using the Itti-Koch saliency computational model (Itti et al., 1998; Yoshida et al., 2012). This approach may provide a better understanding of the abnormalities in attentional processes found in schizophrenic subjects and marmoset ketamine models may replicate the effect on visual salience reported in human studies (Yoshida et al., 2015). In addition to static images, movie stimulus can be adopted (Chen et al., 2021). Because movie stimuli have been reported to act simultaneously and reliably on multiple cortical processing systems (Bartels and Zeki, 2004), this approach is beneficial for neurophysiological studies of visual salience. Electrophysiological and Ca-imaging techniques that are already established for marmosets (Ebina et al., 2018; Kondo et al., 2018; Ghahremani et al., 2019; Selvanayagam et al., 2019), in combination with our experimental protocol, can be used to study neuronal activities at the cell level in pharmacological models of schizophrenia. The application of multiple approaches and the creation of a database for researchers and medical doctors can facilitate the development of both neuroscience and psychiatry (Onitsuka et al., 2022b).

Clinical significance

Direct comparison of eye movement characteristics in marmoset and human studies conducted on the identical image data set strengthened the validity of the pharmacological model of schizophrenia. The experimental protocol developed here can be used in various neuroscientific studies on animals to provide fundamental insight into pathophysiology. Based on the findings from animal studies, new, effective eye tracking tests could be created for the early diagnosis of psychiatric disorders, including schizophrenia (Onitsuka et al., 2022a). In addition, research with non-human primates can suggest effective pharmacological manipulation strategies and targets, which can facilitate drug development and consequently contribute to psychiatry.

Conclusion

In conclusion, our findings demonstrate that the subanesthetic administration of ketamine is a promising pharmacological model of schizophrenia in common marmosets and can be used in combination with free-viewing paradigms to establish “translatable markers” (Okano et al., 2016) of schizophrenia in non-human primates. Our experimental protocol can be further used for electrophysiological or optical imaging studies in simulated schizophrenic states and can provide valuable insights into the involvement of certain brain regions in cognitive functions and visual information processing. In addition, details of changes in neuronal activity in pathophysiological conditions could lead to a better understanding of the underlying processes behind neuropsychiatric disorders.

Data availability statement

The original contributions presented in this study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The animal study was reviewed and approved by Hokkaido University Animal Care and Use Committee.

Author contributions

MY designed the research. MY and ZP performed the research and analyzed the data and wrote the manuscript. MI and RH provided the images. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by AMED (JP22dm0207069) and JSPS KAKENHI (JP22H02936 and JP20H00001) to MY.

Conflict of interest

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Andreasen, N. (1995). Symptoms, signs, and diagnosis of schizophrenia. Lancet 346, 477–481. doi: 10.1016/S0140-6736(95)91325-4

CrossRef Full Text | Google Scholar

Arolt, V., Teichert, H.-M., Steege, D., Lencer, R., and Heide, W. (1998). Distinguishing schizophrenic patients from healthy controls by quantitative measurement of eye movement parameters. Biol. Psychiatry 44, 448–458. doi: 10.1016/S0006-3223(97)00479-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartels, A., and Zeki, S. (2004). Functional Brain Mapping during Free Viewing of Natural Scenes. Hum. Brain Mapp. 21, 75–85. doi: 10.1002/hbm.10153

PubMed Abstract | CrossRef Full Text | Google Scholar

Beck, K., Hindley, G., Borgan, F., Ginestet, C., McCutcheon, R., Brugger, S., et al. (2020). Association of Ketamine with Psychiatric Symptoms and Implications for Its Therapeutic Use and for Understanding Schizophrenia: A Systematic Review and Meta-analysis. JAMA Netw. Open 3, 1–20. doi: 10.1001/jamanetworkopen.2020.4693

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Benson, P. J., Leonards, U., Lothian, R. M., St Clair, D. M., and Merlo, M. C. G. (2007). Visual scan paths in first-episode schizophrenia and cannabis-induced psychosis. J. Psychiatry Neurosci. 32, 267–274.

Google Scholar

Bestelmeyer, P. E. G., Tatler, B. W., Phillips, L. H., Fraser, G., Benson, P. J., and St.Clair, D. (2006). Global visual scanning abnormalities in schizophrenia and bipolar disorder. Schizophr. Res. 87, 212–222. doi: 10.1016/j.schres.2006.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Campana, A., Duci, A., Gambini, O., and Scarone, S. (1999). An artificial neural network that uses eye-tracking performance to identify patients with schizophrenia. Schizophr. Bull. 25, 789–799. doi: 10.1093/oxfordjournals.schbul.a033419

PubMed Abstract | CrossRef Full Text | Google Scholar

Carpenter, W. T., and Koenig, J. I. (2008). The Evolution of Drug Development in Schizophrenia. Neuropsychopharmacology 33, 2061–2079. doi: 10.1038/sj.npp.1301639

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. Y., and Hafed, Z. M. (2013). Postmicrosaccadic enhancement of slow eye movements. Ann. Intern. Med. 158, 5375–5386. doi: 10.1523/JNEUROSCI.3703-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. Y., Matrov, D., Veale, R., Onoe, H., Yoshida, M., Miura, K., et al. (2021). Properties of visually guided saccadic behavior and bottom-up attention in marmoset, macaque, and human. J. Neurophysiol. 125, 437–457. doi: 10.1152/jn.00312.2020

PubMed Abstract | CrossRef Full Text | Google Scholar

Ebina, T., Masamizu, Y., Tanaka, Y. R., Watakabe, A., Hirakawa, R., Hirayama, Y., et al. (2018). Two-photon imaging of neuronal activity in motor cortex of marmosets during upper-limb movement tasks. Nat. Commun. 9:1879. doi: 10.1038/s41467-018-04286-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Featherstone, R. E., Kapur, S., and Fletcher, P. J. (2007). The amphetamine-induced sensitized state as a model of schizophrenia. Prog. Neuro Psychopharmacol. Biol. Psychiatry 31, 1556–1571. doi: 10.1016/j.pnpbp.2007.08.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghahremani, M., Johnston, K. D., Ma, L., Hayrynen, L. K., and Everling, S. (2019). Electrical microstimulation evokes saccades in posterior parietal cortex of common marmosets. J. Neurophysiol 122, 1765–1776. doi: 10.1152/jn.00417.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibaldi, A., and Sabatini, S. P. (2021). The saccade main sequence revised: A fast and repeatable tool for oculomotor analysis. Behav. Res. Methods 53, 167–187. doi: 10.3758/s13428-020-01388-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashimoto, R. (2021). Do eye movement abnormalities in schizophrenia cause Praecox Gefühl? Psychiatry Clin. Neurosci. 75, 79–80. doi: 10.1111/pcn.13197

PubMed Abstract | CrossRef Full Text | Google Scholar

Holzman, P. S., Proctor, L. R., and Hughes, D. W. (1973). Eye-Tracking Patterns in Schizophrenia. Science 181, 179–181. doi: 10.1126/science.181.4095.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Hung, C. C., Yen, C. C., Ciuchta, J. L., Papoti, D., Bock, N. A., Leopold, D. A., et al. (2015). Functional mapping of face-selective regions in the extrastriate visual cortex of the marmoset. J. Neurosci. 35, 1160–1172. doi: 10.1523/JNEUROSCI.2659-14.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Itti, L., Koch, C., and Niebur, E. (1998). A model of saliency-based visual attention for rapid scene analysis. IEEE Trans. Pattern. Anal. Mach. Intell. 20, 1254–1259. doi: 10.1109/34.730558

CrossRef Full Text | Google Scholar

Jentsch, J. D., and Roth, R. H. (1999). The neuropsychopharmacology of phencyclidine: From NMDA receptor hypofunction to the dopamine hypothesis of schizophrenia. Neuropsychopharmacology 20, 201–225. doi: 10.1016/S0893-133X(98)00060-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, C., Watson, D., and Fone, K. (2011). Animal models of schizophrenia. Br. J. Pharmacol. 164, 1162–1194. doi: 10.1111/j.1476-5381.2011.01386.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaneko, T., Komatsu, M., Yamamori, T., Ichinohe, N., and Okano, H. (2022). Cortical neural dynamics unveil the rhythm of natural visual behavior in marmosets. Commun. Biol. 5:108. doi: 10.1038/s42003-022-03052-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kojima, T., Matsushima, E., Ando, K., Ando, H., Sakurada, M., Ohta, K., et al. (1992). Exploratory eye movements and neuropsychological tests in schizophrenic patients. Schizophr. Bull. 18, 85–94. doi: 10.1093/schbul/18.1.85

PubMed Abstract | CrossRef Full Text | Google Scholar

Kondo, T., Saito, R., Otaka, M., Yoshino-Saito, K., Yamanaka, A., Yamamori, T., et al. (2018). Calcium Transient Dynamics of Neural Ensembles in the Primary Motor Cortex of Naturally Behaving Monkeys. Cell Rep. 24, 2191–2195. doi: 10.1016/j.celrep.2018.07.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotani, M., Enomoto, T., Murai, T., Nakako, T., Iwamura, Y., Kiyoshi, A., et al. (2016). The atypical antipsychotic blonanserin reverses (+)-PD-128907- and ketamine-induced deficit in executive function in common marmosets. Behav. Brain Res. 305, 212–217. doi: 10.1016/j.bbr.2016.02.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Krystal, J. H., Karper, L. P., Seibyl, J. P., Freeman, G. K., Delaney, R., Bremner, J. D., et al. (1994). Subanesthetic Effects of the Noncompetitive NMDA Antagonist. Ketamine, in Humans. Arch Gen Psychiatry 51, 199–214. doi: 10.1001/archpsyc.1994.03950030035004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, P. J., Bradley, M. M., and Cuthbert, B. N. (2008). International Affective Picture System (IAPS): Digitized Photographs, Instruction Manual and Affective Ratings, Technical Report A-8. Gainesville, FL: University of Florida.

Google Scholar

Leigh, R. J., and Zee, D. S. (2015). The Neurology of Eye Movements. Oxford: Oxford University Press. doi: 10.1093/med/9780199969289.001.0001

CrossRef Full Text | Google Scholar

Levy, D. L., Holzman, P. S., Matthysse, S., and Mendell, N. R. (1993). Eye tracking dysfunction and schizophrenia: A critical perspective. Schizophr. Bull. 19, 461–536. doi: 10.1093/schbul/19.3.461

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumoto, D., and Ekman, P. (1988). Japanese and Caucasian facial expressions of emotion (JACFEE) [Slides]. San Francisco, CA: San Francisco State University.

Google Scholar

Mitchell, J. F., and Leopold, D. A. (2018). The marmoset as a model for visual neuroscience. Neurosci. Res. 93, 20–46. doi: 10.1016/j.neures.2015.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, J. F., Priebe, N. J., and Miller, C. T. (2015). Motion dependence of smooth pursuit eye movements in the marmoset. J. Neurophysiol. 113, 3954–3960. doi: 10.1152/jn.00197.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, J. F., Reynolds, J. H., and Miller, C. T. (2014). Active vision in marmosets: A model system for visual neuroscience. J. Neurosci. 34, 1183–1194. doi: 10.1523/JNEUROSCI.3899-13.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Miura, K., Hashimoto, R., Fujimoto, M., Yamamori, H., Yasuda, Y., Ohi, K., et al. (2014). An integrated eye movement score as a neurophysiological marker of schizophrenia. Schizophr. Res. 160, 228–229. doi: 10.1016/j.schres.2014.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Morita, K., Miura, K., Fujimoto, M., Yamamori, H., Yasuda, Y., Iwase, M., et al. (2017). Eye movement as a biomarker of schizophrenia: Using an integrated eye movement score. Psychiatry Clin. Neurosci. 71, 104–114. doi: 10.1111/pcn.12460

PubMed Abstract | CrossRef Full Text | Google Scholar

Morita, K., Miura, K., Fujimoto, M., Yamamori, H., Yasuda, Y., Kudo, N., et al. (2019). Eye movement abnormalities and their association with cognitive impairments in schizophrenia. Schizophr. Res. 209, 255–262. doi: 10.1016/j.schres.2018.12.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Okano, H., Hikishima, K., Iriki, A., and Sasaki, E. (2012). The common marmoset as a novel animal model system for biomedical and neuroscience research applications. Semin. Fetal. Neonatal. Med. 17, 336–340. doi: 10.1016/j.siny.2012.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Okano, H., Sasaki, E., Yamamori, T., Iriki, A., Shimogori, T., Yamaguchi, Y., et al. (2016). Brain/MINDS: A Japanese National Brain Project for Marmoset Neuroscience. Neuron 92, 582–590. doi: 10.1016/j.neuron.2016.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Onitsuka, T., Hirano, Y., Nemoto, K., Hashimoto, N., Kushima, I., Koshiyama, D., et al. (2022b). Trends in big data analyses by multicenter collaborative translational research in psychiatry. Psychiatry Clin. Neurosci. 76, 1–14. doi: 10.1111/pcn.13311

PubMed Abstract | CrossRef Full Text | Google Scholar

Onitsuka, T., Hirano, Y., Nakazawa, T., Ichihashi, K., Miura, K., Inada, K., et al. (2022a). Toward recovery in schizophrenia: Current concepts, findings, and future research directions. Psychiatry Clin. Neurosci. 76, 282–291. doi: 10.1111/pcn.13342

PubMed Abstract | CrossRef Full Text | Google Scholar

Pouget, P., Wattiez, N., Rivaud-Péchoux, S., and Gaymard, B. (2010). Rapid development of tolerance to sub-anaesthetic dose of ketamine: An oculomotor study in macaque monkeys. Psychopharmacology 209, 313–318. doi: 10.1007/s00213-010-1797-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Radant, A. D., Bowdle, T. A., Cowley, D. S., Kharasch, E. D., and Roy-Byrne, P. P. (1998). Does ketamine-mediated N-methyl-D-aspartate receptor antagonism cause schizophrenia-like oculomotor abnormalities? Neuropsychopharmacology 19, 434–444. doi: 10.1016/S0893-133X(98)00030-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Romano, J., Kromrey, J. D., Coraggio, J., Skowronek, J., and Devine, L. (2006). “Exploring methods for evaluating group differences on the NSSE and other surveys: Are the t-test and Cohen’s d indices the most appropriate choices?,” in Paper Presented at the Annual Meeting of the Southern Association for Institutional Research, (Arlington, VA), 14–17.

Google Scholar

Sadakane, O., Masamizu, Y., Watakabe, A., Terada, S. I., Ohtsuka, M., Takaji, M., et al. (2015). Long-Term Two-Photon Calcium Imaging of Neuronal Populations with Subcellular Resolution in Adult Non-human Primates. Cell Rep. 13, 1989–1999. doi: 10.1016/j.celrep.2015.10.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawagashira, R., and Tanaka, M. (2021). Ketamine-induced alteration of working memory utility during oculomotor foraging task in monkeys. eNeuro 8:ENEURO.0403-20.2021. doi: 10.1523/ENEURO.0403-20.2021

PubMed Abstract | CrossRef Full Text | Google Scholar

Schielke, A., and Krekelberg, B. (2021). N-methyl d-aspartate receptor hypofunction reduces visual contextual integration. J. Vis. 21:9. doi: 10.1167/jov.21.6.9

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiller, P. H., and Tehovnik, E. J. (2005). Neural mechanisms underlying target selection with saccadic eye movements. Prog. Brain Res. 149, 157–171. doi: 10.1016/S0079-6123(05)49012-3

CrossRef Full Text | Google Scholar

Selvanayagam, J., Johnston, K. D., Schaeffer, D. J., Hayrynen, L. K., and Everling, S. (2019). Functional localization of the frontal eye fields in the common marmoset using microstimulation. J. Neurosci. 39, 9197–9206. doi: 10.1523/JNEUROSCI.1786-19.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Selvanayagam, J., Johnston, K. D., Wong, R. K., Schaeffer, D., and Everling, S. (2021). Ketamine disrupts gaze patterns during face viewing in the common marmoset. J. Neurophysiol. 126, 330–339. doi: 10.1152/jn.00078.2021

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiler, M. A., Thaker, G. K., Lahti, A. C., and Tamminga, C. A. (2000). Ketamine effects on eye movements. Neuropsychopharmacology 23, 645–653. doi: 10.1016/S0893-133X(00)00156-1

CrossRef Full Text | Google Scholar

Yamada, Y., Matsumoto, Y., Okahara, N., and Mikoshiba, K. (2016). Chronic multiscale imaging of neuronal activity in the awake common marmoset. Sci. Rep. 6:35722. doi: 10.1038/srep35722

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Itti, L., Berg, D. J., Ikeda, T., Kato, R., Takaura, K., et al. (2012). Residual attention guidance in blindsight monkeys watching complex natural scenes. Curr. Biol. 22, 1429–1434. doi: 10.1016/j.cub.2012.05.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, M., Miura, K., Hashimoto, R., Fujimoto, M., Yamamori, H., Yasuda, Y., et al. (2015). Saliency-guided eye movement during free-viewing in schizophrenic patients. J. Vis. 15:61. doi: 10.1167/15.12.61

CrossRef Full Text | Google Scholar

Zarate, C. A., Singh, J. B., Carlson, P. J., Brutsche, N. E., Ameli, R., Luckenbaugh, D. A., et al. (2006). A Randomized Trial of an N-methyl-D-aspartate Antagonist in Treatment-Resistant Major Depression. Arch. Gen. Psychiatry 63:856. doi: 10.1001/archpsyc.63.8.856

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: pharmacological model, ketamine, schizophrenia, scanpath, saccades, non-human primates

Citation: Polyakova Z, Iwase M, Hashimoto R and Yoshida M (2022) The effect of ketamine on eye movement characteristics during free-viewing of natural images in common marmosets. Front. Neurosci. 16:1012300. doi: 10.3389/fnins.2022.1012300

Received: 05 August 2022; Accepted: 29 August 2022;
Published: 20 September 2022.

Edited by:

Cheng Xue, The University of Chicago, United States

Reviewed by:

Gongchen Yu, National Institutes of Health (NIH), United States
Bart Krekelberg, Rutgers University–Newark, United States

Copyright © 2022 Polyakova, Iwase, Hashimoto and Yoshida. 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) and the copyright owner(s) 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: Masatoshi Yoshida, myoshi@chain.hokudai.ac.jp

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.