Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 21 May 2021
Sec. Neurotrauma
This article is part of the Research Topic Recent Advances in the Understanding of Mild Traumatic Brain Injury and Post-Concussive Syndrome: Changing the Paradigm View all 9 articles

Behavioral and Histopathological Impairments Caused by Topical Exposure of the Rat Brain to Mild-Impulse Laser-Induced Shock Waves: Impulse Dependency

\nMotoyuki JitsuMotoyuki Jitsu1Katsuki NiwaKatsuki Niwa1Go SuzukiGo Suzuki1Takeyuki ObaraTakeyuki Obara1Yukiko IwamaYukiko Iwama1Kohsuke HagisawaKohsuke Hagisawa1Yukihiro TakahashiYukihiro Takahashi1Yoshitaro MatsushitaYoshitaro Matsushita1Satoru TakeuchiSatoru Takeuchi2Hiroshi NawashiroHiroshi Nawashiro2Shunichi SatoShunichi Sato3Satoko Kawauchi
Satoko Kawauchi3*
  • 1Military Medicine Research Unit, Japan Ground Self Defense Force, Tokyo, Japan
  • 2Department of Neurosurgery, National Defense Medical College, Tokorozawa, Japan
  • 3Division of Bioinformation and Therapeutic Systems, National Defense Medical College Research Institute, Tokorozawa, Japan

Although an enormous number of animal studies on blast-induced traumatic brain injury (bTBI) have been conducted, there still remain many uncertain issues in its neuropathology and mechanisms. This is partially due to the complex and hence difficult experimental environment settings, e.g., to minimize the effects of blast winds (tertiary mechanism) and to separate the effects of brain exposure and torso exposure. Since a laser-induced shock wave (LISW) is free from dynamic pressure and its energy is spatially well confined, the effects of pure shock wave exposure (primary mechanism) solely on the brain can be examined by using an LISW. In this study, we applied a set of four LISWs in the impulse range of 15–71 Pa·s to the rat brain through the intact scalp and skull; the interval between each exposure was ~5 s. For the rats, we conducted locomotor activity, elevated plus maze and forced swimming tests. Axonal injury in the brain was also examined by histological analysis using Bodian silver staining. Only the rats with exposure at higher impulses of 54 and 71 Pa·s showed significantly lower spontaneous movements at 1 and 2 days post-exposure by the locomotor activity test, but after 3 days post-exposure, they had recovered. At 7 days post-exposure, however, these rats (54 and 71 Pa·s) showed significantly higher levels of anxiety-related and depression-like behaviors by the elevated plus maze test and forced swimming test, respectively. To the best of the authors' knowledge, there have been few studies in which a rat model showed both anxiety-related and depression-like behaviors caused by blast or shock wave exposure. At that time point (7 days post-exposure), histological analysis showed significant decreases in axonal density in the cingulum bundle and corpus callosum in impulse-dependent manners; axons in the cingulum bundle were found to be more affected by a shock wave. Correlation analysis showed a statistically significant correlation between the depression like-behavior and axonal density reduction in the cingulum bundle. The results demonstrated the dependence of behavior deficits and axonal injury on the shock wave impulse loaded on the brain.

Introduction

The risk of suffering from blast-induced traumatic brain injury (bTBI) continues not only for military personnel but also for civilians due to frequent attacks using improvised explosive devices (IEDs) (14). Most of the patients lack any abnormality in conventional neuroimaging for diagnosis, usually computed tomography, and they have therefore been diagnosed as having mild bTBI (57). However, many of them have developed chronic, persistent neuropsychiatric symptoms (810) as well as many other symptoms such as headache (11) and sleep disorder (12), causing a serious decrease in their quality of life. Although extensive studies have been conducted, the neuropathology and mechanisms of bTBI have not been fully elucidated (13). This is partially due to the complex and hence difficult experimental environment settings. In animal studies, it is important to examine the primary mechanism (effects of a shock wave itself) since the above-mentioned bTBI-related symptoms have been observed in patients free from interactions with the secondary mechanism (effects of propelled debris and shrapnel) and the tertiary mechanism (effects of acceleration due to blast winds) (14). Thus, the effects of wind or a jet should be carefully minimized. Although bTBIs are injuries due to systemic exposure in most cases (14), the outcomes of systemic exposure are too complex to analyze. In animal studies, therefore, it is occasionally necessary to separate the effects of brain exposure from those of torso exposure to examine the effects of brain exposure alone.

We have been using laser-induced shock waves (LISWs) to investigate the mechanisms of shock wave-induced brain injury (1517). An LISW is a pure shock wave that is free from wind or a jet. In addition, the energy of an LISW is spatially well confined and its size is controllable by changing the laser spot size, enabling site-specific shock wave application. Thus, we can exclude both the effects of the tertiary mechanism and torso exposure in animal experiments by using LISWs.

The purpose of this study was to examine behavioral and neuropathological changes caused by exposure of the brain alone to LISWs as a function of shock wave impulse in rats. Impulse is defined as the time-integrated positive pressure component of a shock wave and is known to be one of the most important parameters to determine shock wave-induced tissue damage (18). Considering neuropsychiatric symptoms that have been reported in bTBI patients (19, 20), we assessed locomotor activity, anxiety-related behavior and depression-like behavior for the rats. We also examined axonal injury in the brain by histological analysis using Bodian silver staining; diffuse axonal injury has also been reported in bTBI patients (8, 21), although it is not clear how the axonal injury is unique to bTBI. Associations between the findings obtained from behavioral tests and those obtained from pathological assessment are discussed.

Materials and Methods

All requests for animals and procedures intended to be used in the present study were approved by the Ethics Committee of Animal Care and Experimentation, National Defense Medical College, Japan (Permission numbers: 10040, 13024, and 16036).

Animals

Mice and rats have been most widely used for studying the pathology and mechanisms of bTBI. While the genetic background has been characterized in more detail for mice, rats have a larger forebrain and show greater sophistication and sociability than mice do. Thus, it has been suggested that rats can provide more comparability with humans in neurobehavioral studies (22, 23). Considering this factor, we chose rats as model animals for the present study. Nien-week-old male Sprague-Dawley rats were obtained from Japan SLC, Inc. (Shizuoka, Japan) and were housed in standard laboratory cages on a 12:12-h light/dark cycle with free access to food and water after arrival. Rats at 10 weeks of age weighing 300–340 g (n = 100) were divided into a control group without LISW exposure (n = 16) and five LISW groups (n = 84); the shock wave conditions are described in detail later.

Generation and Characteristics of Laser-Induced Shock Waves

A method for generating an LISW was described previously (15) (Figure 1A). Briefly, a laser target, which was a light-absorbing material (0.5-mm-thick natural black rubber disk) covered with an optically transparent material (1.0-mm-thick polyethylene terephthalate [PET] sheet) was irradiated with a high-intensity laser pulse from a Q-switched Nd:YAG laser (pulse width, 6 ns FWHM; wavelength, 532 nm; Brilliant b; Quantel, Les Ulis Cedex, France). The laser pulse was absorbed by the black rubber to induce plasma and its expansion was accompanied by the generation of a shock wave (LISW). Peak pressure and impulse of an LISW can be precisely controlled by changing the laser fluence on the target. The size of an LISW can also be controlled by changing the spot size (diameter) of the laser beam.

FIGURE 1
www.frontiersin.org

Figure 1. Generation and characteristics of laser-induced shock waves (LISWs) and their placement in rat brains. (A) Setup for generating LISW. (B) Typical temporal waveform of LISW (~82 MPa) with a Friedlander waveform, P = Ps(e(−t/t*))(1–(t/t*)) (Ps = 90 MPa, t* = 0.8 μs). (C) Positions of four LISWs transcutaneously applied to the rat brain, which are shown with respect to the skull anatomy.

A problem in using LISWs for bTBI study is that the duration of positive pressure (hereafter simply called duration) of an LISW is as short as ~1 μs (Figure 1B), which is much shorter than typical IED-related shock waves (18, 24). In most of the shock tube-based bTBI studies using rodents, on the other hand, shock waves with durations of milliseconds have been used to mimic IED-related shock waves. It should be noted, however, that the interaction of real IED shock waves with a human brain and that with rodents' brains are totally different in terms of pressure gradient formed within the brain and boundary effects due to anatomical differences of the brains. The effects of shock wave duration was discussed in the Discussion section. In this study, we assumed that the impulse of the shock wave is a primary parameter to determine brain injury. Alley estimated impulses of IED-related shock waves to range from ~6.9 to ~100 Pa·s at propagation distances of 1–10 m (18). We produced 3-mm diameter LISWs with impulses ranging from 15 to 71 Pa·s by changing the laser fluence on the target in the range of 0.6~3.0 J/cm2. This impulse range was within the conditions that were shown to cause mild bTBI in rats by using a shock tube (<250 Pa·s) in a study by Mishra et al. (25). As described in section Neuropathology, the rats were transcardially perfused at 7 days post-exposure, and the brains were extracted. No brain hemorrhages were seen upon visual observation at any of the impulses used.

LISW Application

Rats were anesthetized with pentobarbital sodium (i.p. 50 mg/kg) and their heads were shaved. A rat in the LISW groups was subjected to a set of four LISWs at a given impulse on the following sites to cover the whole limbic system (Figure 1C): 5 mm bilateral from a point at the intersection of the auricular line with the mid-sagittal line and then 5 mm anterior from each point. The interval between each exposure was ~5 s. For each LISW application, a laser target was placed on the scalp with forceps; ultrasound gel was applied between the bottom of the target (rubber) and the scalp for acoustic impedance matching. Rats in the LISW groups were exposed to a set of LISWs generated at laser fluences of 0.6, 1.2, 1.8, 2.4, and 3.0 J/cm2, the corresponding impulses (peak pressures) of LISWs being 15 Pa·s (33 MPa), 29 Pa·s (52 MPa), 43 Pa·s (66 MPa), 54 Pa·s (76 MPa), and 71 Pa·s (82 MPa), respectively. Each rat underwent either behavioral tests or pathological assessments. The number of animals in each group is stated in the Results section. Rats in the control group (sham controls) were treated in the same way except for exposure to LISWs. They received the same anesthesia and preparation as those the rats of LISW groups did. LISWs were generated at the position ~10 cm away from the rats but were not applied to them.

Behavioral Assessments

To assess potential behavioral changes after LISW exposure, rats underwent (i) a locomotor activity test, (ii) an elevated plus-maze test, and (iii) a forced swimming test to assess spontaneous activity, anxiety-related behavior (26) and depression-like behavior (27), respectively. Experimenter blinding was not done, but all behavioral experiment data were collected by recording systems and analyzed with a dedicated software, ensuring the objectivity of our results.

Locomotor activity of rats was measured in a 12-h dark cycle using a Supermex apparatus (Muromachi Kikai Co., Ltd., Tokyo, Japan) for 5 consecutive days from 1 day post-exposure. With the apparatus, spontaneous motor activities of individual rats can be assessed by detecting heat radiation from the body with a far-infrared sensor attached to the cage (28).

At 7 days after LISW exposure, rats were subjected to an elevated plus maze test (Model EPM-04R, Muromachi Kikai Co., Ltd., Tokyo, Japan). The open arms and closed arms were both 50 cm in length and 10 cm in width; the wall height and stand height were 40 and 50 cm, respectively. After acclimatization to the experimental room for at least 1 h, each rat was placed on the center platform facing one of the two open arms and allowed to behave freely for 5 min. Behavior was recorded with a video tracking system (DVTrack Video Tracking System CompACT VAS/DV; Muromachi Kikai Co., Ltd., Tokyo, Japan), and the percentage of time spent in the open arms to the total time and total distance traveled were analyzed.

Thereafter, rats were allowed to rest for more than 1 h and were then subjected to a forced swimming test using a cylindrical tank 60 cm high and 30 cm in diameter. We based on a published protocol in which a single session was used and results for the first vigorous phase was excluded (29). Each rat was placed in 25°C water in a plastic cylinder (water depth, 35 cm; diameter, 30 cm) and allowed to behave freely for 15 min. For the last 10 min, the time spent immobile was measured and its percentage to the total time (10 min) was calculated, for which ImageJ PS1, a software package modified on the basis of the public domain ImageJ program (developed at the U.S. National Institutes of Health and available at: http://rsb.info.nih.gov/ij) (30) (O'Hara and Co., Ltd., Tokyo, Japan), was used.

Neuropathology

At 7 days post-exposure, rats were deeply anesthetized with pentobarbital sodium (i.p. 100 mg/kg) and were transcardially perfused with normal saline followed by 4% buffered paraformaldehyde. Brains were removed after perfusion and immersed in 4% buffered paraformaldehyde overnight and then embedded in paraffin. Five-μm-thick coronal sections at 3.3 mm posterior to the bregma were used for Bodian silver staining. Axonal densities were quantified by assessing light transmittance for digital images of the brain specimens after staining (31) using ImageJ software described above. Images were converted into 8-bit gray scale images and the mean light transmittance within the regions of interest was determined (arbitrary units, ranging from 0 to 255). High values (increased light transmittance) correspond to low staining intensity, indicating decreased axonal densities and hence axonal injury.

Statistical Analysis

Data are presented as means ± standard deviation (SD). The results of the spontaneous movement test were analyzed by two-way repeated measures ANOVA with the Bonferroni post-hoc test. The results of the elevated plus maze test and forced swimming test and the light transmittance of brain specimens were analyzed by one-way factorial ANOVA with Dunnett's post-hoc test. The correlation between the result of the light transmittance of brain specimens and the results of the behavioral tests, i.e., the elevated plus maze test and the forced swimming test, were analyzed by Pearson's correlation coefficients. These statistical analyses were performed using GraphPad Prism software (Prism 6 for Windows, GraphPad software, Inc., La Jolla, CA, USA). A difference was considered statistically significant when p ≤ 0.05. To assess the validity of the sample sizes used in behavioral tests, post-hoc power analysis was performed using a program G*Power 3 (32).

Results

Spontaneous Locomotor Activity

Figure 2 shows spontaneous motor activities during 12-h dark cycles at 1–5 days post-exposure for rats in the control group and LISW groups with exposure at five different impulses (n = 10 for control [0 Pa·s], 15, 29, and 43 Pa·s; n = 5 for 54 and 71 Pa·s). The values (mean and SD) and results of the statistical analysis are summarized in Tables 1, 2. The two-way repeated measures ANOVA showed that shock wave impulses and days post-exposure had significant effects on locomotor activities; significant interaction effects were also found (p = 0.0042). Bonferroni post-hoc test showed that spontaneous motor activities were significantly decreased at 1 and 2 days post-exposure at the higher impulses of 54 and 71 Pa·s compared to those of the control group. After 3 days post-exposure, however, there were no significant differences in spontaneous motor activities between any of the LISW groups and the control group, indicating recovery of locomotor activity for LISW groups with exposure at the higher impulses.

FIGURE 2
www.frontiersin.org

Figure 2. Spontaneous motor activities in the locomotor activity test as a function of post-exposure days (Day 1–5) for rats whose brains were exposed to LISWs with different impulses (n = 10 for 0 [control], 15, 29, and 43 Pa·s; n = 5 for 54 and 71 Pa·s). The data are presented as means ± SD. Asterisks (***p ≤ 0.001) indicate statistically significant differences (two-way repeated measures ANOVA with Bonferroni post-hoc test, vs. control).

TABLE 1
www.frontiersin.org

Table 1. Summary of the two-way repeated measures ANOVA for the results of locomotor activity test.

TABLE 2
www.frontiersin.org

Table 2. Summary of the results of locomotor activity test and statistical analysis with Bonferroni post-hoc test.

Anxiety-Related Behavior (Elevated Plus Maze Test)

Figure 3 shows the results of the elevated plus maze test for all groups of rats at 7 days post-exposure (n = 10 for control [0 Pa·s], 15, 29 and 43 Pa·s; n = 5 for 54 and 71 Pa·s), from which data for rats that fell from the arms were excluded (33). One-way factorial ANOVA showed a significant difference in the percentages of time spent in open arms (ANOVA F-value = 6.968; p = 0.0001) but not in the total distances traveled among the groups (ANOVA F-value = 0.6305; p = 0.6775). The post-hoc test results are summarized in Table 3. For rats in the control group, the percentage of time spent in the open arms was 32.7 ± 7.3% (n = 9). The percentages of time spent in the open arms for rats in the LISW groups with exposure at 54 and 71 Pa·s were significantly lower [6.1 ± 4.9% (n = 5) and 12.1 ± 3.9% (n = 5), respectively], indicating a significant increase in anxiety-related behavior for these groups of rats. Rats in the LISW groups with exposure at the lower impulses (15, 29, and 43 Pa·s) showed no significant decrease in the time spent in the open arms. On the other hand, no significance differences were found in the total distances traveled by all rat groups in elevated plus maze test, suggesting no significant changes in their locomotor activities.

FIGURE 3
www.frontiersin.org

Figure 3. (A) Percentage of time spent in open arms and (B) total distance traveled in the elevated plus maze test at 7 days post-exposure as a function of impulses of LISWs that were applied to the rat brains (mean ± SD). Rats that fell from the arms during the test were excluded; the numbers of those rats excluded were 1, 2, 2, 1, 0, and 0 in the sham control group, 15, 29, 43, 54, and 71 Pa·s group, respectively. Asterisks (***p ≤ 0.001) indicate statistically significant differences (one-way factorial ANOVA with Dunnett's post-hoc test, vs. control).

TABLE 3
www.frontiersin.org

Table 3. Summary of the results of elevated plus maze test and statistical analysis with Dunnett's post-hoc test.

Depression-Like Behavior (Forced Swimming Test)

Figure 4 shows the results of the forced swimming test for all groups of rats at 7 days post- exposure (n = 10 for control [0 Pa·s], 15, 29 and 43 Pa·s; n = 5 for 54 and 71 Pa·s), from which data for rats rescued from drowning during the test were excluded (34). No rats were excluded due to visually recognized health problems. One-way factorial ANOVA showed a significant difference in the percentages of immobility time among the groups (ANOVA F-value = 19.3; p < 0.0001). The results of the post-hoc tests are summarized in Table 4. For rats in the control group, the percentage of immobility time was 35.8 ± 6.2% (n = 9). The percentages of immobility time for rats in the LISW groups with exposure at 51 and 71 Pa·s were 56.6 ± 14.4% (n = 5) and 81.3 ± 15.7% (n = 5), respectively, which were significantly greater than that for the control group. These results indicate significantly higher levels of depression-like behavior in these groups of rats. Rats in the LISW groups with exposure at the lower impulses (15, 29, and 43 Pa·s) showed no significant decrease in immobility time.

FIGURE 4
www.frontiersin.org

Figure 4. Percentage of immobility time in the forced swimming test at 7 days post-exposure as a function of impulses of LISWs that were applied to the rat brains (mean ± SD). Rats that were rescued from drowning during the test were excluded; the numbers of those rats excluded were 1, 1, 1, 0, 0, and 0 in the sham control group, 15, 29, 43, 54, and 71 Pa·s group, respectively. Asterisks (*p < 0.001, ***p ≤ 0.001) indicate statistically significant differences (one-way factorial ANOVA with Dunnett's post-hoc test, vs. control).

TABLE 4
www.frontiersin.org

Table 4. Summary of the results of forced swimming test and statistical analysis with Dunnett's post-hoc test.

Histopathological Assessments

Figures 5a–f shows representative images of Bodian-stained brain sections at a boundary region of the cingulum bundle (CB) and the corpus callosum (CC) in the brains of all groups of rats. The whole image region of the control brain section is relatively dark (Figure 5a). For the sections of brains of rats in the LISW groups, on the other hand, the CB regions are lighter than the CC regions, indicating increased light transmittance and hence decreased axonal density in the CB regions (Figures 5b–f). For all of the images, a unidirectional fibrous structure is seen in the CC regions, while there are no such microstructures in the CB regions, indicating different fiber orientations in the CB and CC.

FIGURE 5
www.frontiersin.org

Figure 5. Representative images of Bodian silver-stained sections of the boundary region between the CB and the CC in brains exposed to LISWs with different impulses and their light transmittance 7 days post-exposure. (a) Control (0 Pa·s). (b) 15 Pa·s. (c) 29 Pa·s. (d) 43 Pa·s. (e) 54 Pa·s. (f) 71 Pa·s. Dashed lines indicate the boundary of the CB and CC. Scale bar indicates 50 μm. Light transmittance of Bodian silver-stained sections of the regions of (g) the CB and (h) the CC in rat brains exposed to LISWs as a function of impulse (n = 6 for control [0 Pa·s]; n = 8 for 15, 29, and 43 Pa·s; n = 10 for 54 and 71 Pa·s). Higher light transmittance indicates lower axonal density. The data are presented as means ± SD. Asterisks (*p ≤ 0.05, ***p ≤ 0.001) indicate statistically significant differences (one-way factorial ANOVA with Dunnett's post-hoc test, vs. control).

To evaluate axonal densities in the CB and CC regions in each image, light transmittance was quantified and the values are shown in Figures 5g,h, respectively (n = 6 for control [0 Pa·s]; n = 8 for 15, 29, and 43 Pa·s; n = 10 for 54 and 71 Pa·s). One-way factorial ANOVA showed a significant difference in light transmittances in CB (ANOVA F-value = 16.98; p < 0.0001) and CC (ANOVA F-value = 4.033; p = 0.0042). The results of the post-hoc test are summarized in Table 5. For the CB region, there is a tendency for increase in light transmittance, and hence axonal density decreases with increase in the impulse of LISWs. Light transmittances for all LISW groups were significantly higher than that for the control group (Figure 5g). At the highest impulse (71 Pa·s), the light transmittance was reduced by ~30% compared with the value for the control group. For the CC region, on the other hand, the extent of increased light transmittance (i.e., decreased axonal density) was smaller than that for the CB region; light transmittance was significantly higher only at the highest impulse (71 Pa·s) than that for the control (Figure 5h).

TABLE 5
www.frontiersin.org

Table 5. Summary of the results of light transmittance of Bodian-stained brain sections and statistical analysis with Dunnett's post-hoc test.

Correlations Between Changes in Axonal Density and Behavioral Outcomes

To investigate if the changes in axonal density in the CB and the CC were correlated with the behavioral outcomes, correlation analyses were performed for the data of mean light transmittance (LT) in CB and in CC, mean percentage of time spent in the open arms in the elevated plus maze test, and mean percentage of immobility time in the forced swimming test. Figure 6 shows relations between the LT in CB and the time spent in open arms (Figure 6A, r = −0.7581, p = 0.0807), between the LT in CB and the immobility time (Figure 6B, r = 0.8360, p = 0.0381), between the LT in CC and the time spent in open arms (Figure 6C, r = −0.4089, p = 0.4208) and between the LT in CC and the immobility time (Figure 6D, r = 0.6103, p = 0.1982). There was a statistically significant correlation only between the LT in CB and the immobility time in the forced swimming test, suggesting that a decreased axonal density was significantly correlated with an increased depression-like behavior.

FIGURE 6
www.frontiersin.org

Figure 6. Results of correlation analysis for the data of behavioral outcomes and axonal densities: Correlations (A) between the light transmittance (LT) in CB and the time spent in open arms in the elevated plus maze test (EPMT), (B) between the LT in CB and the immobility time in the forced swimming test (FST), (C) between the LT in CC and the time spent in open arms in EPMT, and (D) between the LT in CC and the immobility time in FST. Mean data shown in Tables 24 were plotted for the analysis. Correlation coefficient r and p-values are shown in the figure. An asterisk (*p ≤ 0.05) indicate a statistically significant difference.

Verification of Sample Sizes

Since the numbers of rats in LISW groups with 51 and 71 Pa·s exposure were smaller than that in the control group, post-hoc power analyses were performed to verify the sample sizes. The results are shown in Table 6. In all analyses, the power was >0.8, indicating that the sample sizes in this study were acceptable.

TABLE 6
www.frontiersin.org

Table 6. Results of post hoc power analysis for behavioral tests.

Discussion

In this study, behavioral and neuropathological changes were evaluated for rats whose brains were topically exposed to LISWs to examine the mechanisms of shock wave-related TBI. Due to the nature of LISWs, the results are free from the effects of acceleration (tertiary mechanism) and torso exposure, purely reflecting the effect of the primary mechanism on the brain. Our rat model showed significantly lower locomotor activity in the acute phase (1 and 2 days post-exposure) at the higher impulse conditions (54 and 71 Pa·s), indicating loss of basic physical function under the unconstrained conditions. However, their locomotor function was recovered after 3 days post-exposure, and the elevated maze and forced swim tests, which were conducted at 7 days post-exposure, are therefore considered not to be affected by change in the basic physical functions. Importantly, the rats exposed to LISWs at the higher impulse conditions showed significantly higher levels of anxiety-related and depression-like behaviors at 7 days post-exposure. These behavioral changes were shown to be correlated with decreased axonal densities in CB and CC.

Some studies reported the results of anxiety tests using an elevated plus maze or elevated zero maze test for rats exposed to blast overpressures (3539). Kovesdi et al. (35) reported that the time spent in open arms was not significantly decreased in rats at 9 days post-exposure, while it was significantly reduced (anxiety increased) in rats after 46 days in experiments using a shock tube (~20 psi [~138 kPa], single exposure). On the other hand, the experiments of Budde et al. that used a blast tube (100 and 450 kPa, single exposure) showed no significant decreases in open arm time either in the acute phase (1–4 days post-injury) or chronic phase (28–31 days post-injury) (36). Elder et al. (37) found that the rats exposed to a single 74.5-kPa blast per day for three consecutive days exhibited signs of anxiety at 40 days after exposures; rats with exposures moved less and spent significantly more time in the closed arms than in the open arms when compared with controls. Perez-Garcia et al. (38) showed for the same rat model that the blast-exposed rats spent more time in the closed arms even at 29–30 weeks after exposure. Furthermore, Blaze et al. (39) observed an anxiety-like phenotype for the same rat model at 1–1.5 months post-injury, but it was diminished by 12–13 months post-injury. The present results on elevated plus maze test can be featured by the fact that the increased rat anxiety appeared in the acute phase (7 days post-exposure). Unexpectedly, the results of the forced swimming test were limitedly reported for the rodents exposed to a blast. Kawa et al. (40) applied a blast (550 kPa, single exposure) to rats but did not observe an increase in rat immobility time either on day 1, 14, or 35. To the best of the authors' knowledge, there have been few reports on rat bTBI models that showed both anxiety-related and depression-like behaviors. A prospective, longitudinal clinical study on active-duty US military personnel who were diagnosed with mild bTBI showed that, in addition to PTSD-related symptoms, depression and anxiety were the most prominent psychiatric symptoms (20). Therefore, it can be said that anxiety-related and depression-like behaviors should be involved in important behavioral traits of bTBI animal models. Our rat model may be valuable for shock wave-related behavioral changes. The features of our model might be associated with the temporal pressure characteristic of LISWs, which was discussed below.

Goldstein et al. (41) reported that the mice exposed to a blast did not show deficits in learning and memory in the absence of rotational and/or acceleration-deceleration forces. As described above, our model was free from the effects of acceleration (tertiary mechanism) but showed behavioral alterations under high impulse conditions. This might be associated with the difference in shock wave durations: ~5 ms (shock tube) vs. ~1 μs (LISW). Although the animal species were different in the two studies, the pressure gradient for a ~1 μs-duration LISW was much steeper than that for a ~5 ms-duration shock wave, causing a greater local compression-tensile force and/or shear force in the brains.

The CB region, which was found to be affected by shock wave exposure, is the largest commissural white matter bundle in the brain and is involved in the neural circuit known as Papez circuit that controls emotional expression (42). The Papez circuit begins and ends with the hippocampus through the neural pathways via mammillary bodies, anterior thalamic nucleus and cingulum (42). Thus, axonal injury in the CB region can affect not only emotional expression but also memory and cognitive functions in individuals. Anatomically, the circuit exists in the limbic system, which should interact with the incident shock/pressure wave in the present exposure scheme. Accordingly, the observed depression-like behavior in the rats was shown to be significantly correlated with the axonal density reduction found in the CB region (Figure 6B). However, it should be noted that the current analysis is limited by the small sample size. Other structures, such as the frontal cortex and amygdala, are also involved in emotionality (43). Thus, further study is needed to fully understand the mechanisms of the observed behavioral changes.

It is interesting that the loss of axonal density was much more evident in the CB region than in the CC region (Figures 5g,h), although both the CB and CC form white matter bundles and they make contact with each other (the CB being located just above the CC). This may be associated with their axonal orientations perpendicular to each other; the axonal orientation of the CB is sagittal, while that of the CC is transverse (44). Thus, there would be a certain discontinuity of tissue mechanical properties at the CB-CC boundary, resulting in an anisotropic compliance and/or an acoustic impedance mismatching. Since LISWs were applied to the parietal region in the present study, an incident shock/pressure wave propagates from the top of the brain. This propagating pressure wave can be disturbed and reflected at the boundary, possibly generating a shear force or cavitation (2, 45). Such events might be associated with vulnerability of neuronal axons in the CB. An LISW (3 mm in diameter) propagating through the water was visualized by shadowgraphing. The diameter of the high-pressure wavefront was ~4 mm at a depth of 6 mm, indicating the relatively high directionality of LISW propagation (46). Thus, LISW application to other sites of the brain would be useful to examine the abovementioned hypothesis.

In our experiments, locomotor deficits were observed in the rats of LISW groups with exposure at 54 and 71 Pa·s only at 1 and 2 days post-exposure, while anxiety and depression-like behaviors were observed from the behavior tests conducted at 7 days post-exposure. We previously found that LISW application directly caused the occurrence of spreading depolarization, which was followed by a long-lasting oligemia/hypoxemia in the cortex (15). Since the motor cortex was involved in the region with such hemodynamic abnormality, it could affect the locomotion of rats in super-acute phase. An increase in the intracranial pressure was also observed due to the blood brain barrier disruption in this phase, which might also alter the locomotor function of rats.

There were limitations in this study. First, the observation time points were limited, being only in the acute phase. Traumatic neuropathology generally proceeds in a cascading manner, and the clinical study referred to above showed that the symptoms in the 5-year evaluation were worse than those in the 1-year evaluation (20). Thus, it is important to reveal time-dependent changes in neurobehaviors and neuropathology. Second, for neuropathology, we focused only on axonal injury, although we previously reported the increased oxidative stress in the cortex, hippocampus and cerebellum in the same model at 7 days post-exposure (17). Although the clinical neuropathology for bTBI is not fully understood, glial scar formation (47) and aberrant protein accumulation (41) have been reported on the basis of results of studies using postmortem brains of bTBI patients in addition to axonal injury. It is not clear how our neuropathological findings (axonal injury) are associated with those outcomes. Thus, comprehensive analysis of the neuropathology is needed for the rat brain in the chronic phase.

Another important issue to be addressed is the appropriate temporal characteristics of the shock wave to be applied to rodents' brains. In most of the bTBI studies using a blast tube or shock tube with rodents, shock waves with durations of milliseconds have been used to replicate actual IED-relevant shock waves (24). The positive pressure duration of a typical IED explosion-related shock wave ranges from 250 μs to several milliseconds (24). However, the interaction of the human brain with such IED-related shock waves, especially the brain to skull boundary effects, would not be reproduced in the brains of small animals when using similar shock wave duration. The size of the rat brain is roughly one-tenth of that of the human brain (in one-dimensional scale). On the basis of the scaling law, we think that the appropriate shock wave duration to be applied to the rat brain would be one-tenth of the ones described above, ranging from 25 μs to several hundred microseconds. The present LISW (~1-μs duration) was much shorter than this range, but it would still be valid to examine the brain injury associated with underwater blast-related shock waves, of which durations were reported to range from 10 to 200 μs (24). Recently, we have succeeded in controlling the duration of LISWs in a certain range. Thus, we plan to investigate the effects of different shock wave durations on the neuropathological and behavioral changes for the rats, enabling a further understanding of the results of the present study.

Conclusions

We applied a set of four LISWs to rat brains in the impulse range of 15–71 Pa·s and examined neurobehavioral and neuropathological changes in the rats. At 7 days post-exposure, we observed significant anxiety-related and depression-like behaviors in the rats in impulse-dependent manners. Decreased axonal density was observed in both the CB and CC regions, but axons in the CB region were found to be more affected by shock wave exposure. The depression-like behavior was found to be significantly correlated with the axonal density reduction in CB. The rats characterized in this study may be used as a new model to study blast-related TBI.

Data Availability Statement

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

Ethics Statement

The animal study was reviewed and approved by Ethics Committee of Animal Care and Experimentation, National Defense Medical College, Japan.

Author Contributions

MJ, KN, GS, TO, KH, YT, YM, and SS contributed to the conception, design, and administrative support of the study. MJ, KN, GS, TO, and YI performed the research and analyzed the data. MJ, ST, SS, and SK wrote the manuscript. All authors discussed and interpreted the results. All authors reviewed and approved the final manuscript.

Funding

This work was partially supported by the Research Promotion Program for Defense Medicine under the Ministry of Defense.

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.

Acknowledgments

We are grateful staff members of the Division of Behavioral Science and the Division of Environmental Medicine, the National Defense Medical College Research Institute for the use of equipment for behavioral tests.

References

1. Hoge CW, McGurk D, Thomas JL, Cox AL, Engel CC, Castro CA. Mild traumatic brain injury in U.S. Soldiers returning from Iraq. N Engl J Med. (2008) 358:453–63. doi: 10.1056/NEJMoa072972

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Rosenfeld JV, McFarlane AC, Bragge P, Armonda RA, Grimes JB, Ling GS. Blast-related traumatic brain injury. Lancet Neurol. (2013) 12:882–93. doi: 10.1016/S1474-4422(13)70161-3

CrossRef Full Text | Google Scholar

3. UA Army Asia-Pacific Counter-IED Fusion Center Monthly Report - Jane 2020 (2020). Available online at: https://www.usarpac.army.mil/apcied/Monthly/Current.pdf (accessed October 10, 2020).

4. Magnus D, Khan MA, Proud WG. Epidemiology of civilian blast injuries inflicted by terrorist bombings from 1970-2016. Def Technol. (2018) 14:469–76. doi: 10.1016/j.dt.2018.07.014

CrossRef Full Text | Google Scholar

5. Salat DH, Robinson ME, Miller DR, Clark DC, McGlinchey RE. Neuroimaging of deployment-associated traumatic brain injury (TBI) with a focus on mild TBI (mTBI) since 2009. Brain Inj. (2017) 31:1204–19. doi: 10.1080/02699052.2017.1327672

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Yamamoto S, DeWitt DS, Prough DS. Impact & blast traumatic brain injury: implications for therapy. Molecules. (2018) 23:245. doi: 10.3390/molecules23020245

CrossRef Full Text | Google Scholar

7. DePalma RG, Hoffman SW. Combat blast related traumatic brain injury (TBI): decade of recognition; promise of progress. Behav Brain Res. (2018) 340:102–5. doi: 10.1016/j.bbr.2016.08.036

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Mac Donald CL, Johnson AM, Cooper D, Nelson EC, Werner NJ, Shimony JS, et al. Detection of blast-related traumatic brain injury in U.S. military personnel. N Engl J Med. (2011) 364:2091–100. doi: 10.1056/NEJMoa1008069

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nelson NW, Davenport ND, Sponheim SR, Anderson CR. Chapter 32: Blast-related mild traumatic brain injury. In: Kobeissy FH, editor. Brain Neurotrauma: Molecular, Neuropsychological, and Rehabilitation Aspects. Boca Raton, FL: CRC Press/Taylor & Francis (2015) 6–7.

10. Cernak I. Understanding blast-induced neurotrauma: how far have we come?. Concussion. (2017) 2:CNC42. doi: 10.2217/cnc-2017-0006

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Couch JR, Stewart KE. Headache prevalence at 4-11 years after deployment-related traumatic brain injury in veterans of Iraq and Afghanistan wars and comparison to controls: a matched case-controlled study. Headache. (2016) 56:1004–21. doi: 10.1111/head.12837

CrossRef Full Text | Google Scholar

12. Good CH, Brager AJ, Capaldi VF, Mysliwiec V. Sleep in the United States Military. Neuropsychopharmacology. (2020) 45:176–91. doi: 10.1038/s41386-019-0431-7

CrossRef Full Text | Google Scholar

13. Cernak I, Stein DG, Elder GA, Ahlers S, Curley K, DePalma RG, et al. Preclinical modelling of militarily relevant traumatic brain injuries: challenges and recommendations for future directions. Brain Inj. (2017) 31:1168–76. doi: 10.1080/02699052.2016.1274779

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Cernak I, Noble-Haeusslein LJ. Traumatic brain injury: an overview of pathobiology with emphasis on military populations. J Cereb Blood Flow Metab. (2010) 30:255–66. doi: 10.1038/jcbfm.2009.203

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Sato S, Kawauchi S, Okuda W, Nishidate I, Nawashiro H, Tsumatori G. Real-time optical diagnosis of the rat brain exposed to a laser-induced shock wave: observation of spreading depolarization, vasoconstriction and hypoxemia-oligemia. PLoS One. (2014) 9:e82891. doi: 10.1371/journal.pone.0082891

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kawauchi S, Okuda W, Nawashiro H, Sato S, Nishidate I. Multispectral imaging of cortical vascular and hemodynamic responses to a shock wave: observation of spreading depolarization and oxygen supply-demand mismatch. J Biomed Opt. (2019) 24:1–17. doi: 10.1117/1.JBO.24.3.035005

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Maekawa T, Uchida T, Nakata-Horiuchi Y, Kobayashi H, Kawauchi S, Kinoshita M, et al. Oral ascorbic acid 2-glucoside prevents coordination disorder induced via laser-induced shock waves in rat brain. PLoS One. (2020) 15:e0230774. doi: 10.1371/journal.pone.0230774

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Alley M. Explosive Blast Loading Experiments for TBI Scenarios: Characterization and Mitigation. Thesis, West Lafayette, IN: Purdue University (2009).

Google Scholar

19. Spelman JF, Hunt SC, Seal KH, Burgo-Black AL. Post deployment care for returning combat veterans. J Gen Intern Med. (2012) 27:1200–9. doi: 10.1007/s11606-012-2061-1

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mac Donald CL, Barber J, Jordan M, Johnson AM, Dikmen S, Fann JR, et al. Early clinical predictors of 5-year outcome after concussive blast traumatic brain injury. JAMA Neurol. (2017) 74:821–9. doi: 10.1001/jamaneurol.2017.0143

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Magnuson J, Leonessa F, Ling GS. Neuropathology of explosive blast traumatic brain injury. Curr Neurol Neurosci Rep. (2012) 12:570–9. doi: 10.1007/s11910-012-0303-6

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Aitman T, Dhillon P, Geurts AM. A RATional choice for translational research? Dis Model Mech. (2016) 9:1069–72. doi: 10.1242/dmm.027706

CrossRef Full Text | Google Scholar

23. Ellenbroek B, Youn J. Rodent models in neuroscience research: is it a rat race? Dis Model Mech. (2016) 9:1079–87. doi: 10.1242/dmm.026120

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Courtney MW, Courtney AC. Working toward exposure thresholds for blast-induced traumatic brain injury: thoracic and acceleration mechanisms. Neuroimage. (2011) 54(Suppl 1):S55–S61. doi: 10.1016/j.neuroimage.2010.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Mishra V, Skotak M, Schuetz H, Heller A, Haorah J, Chandra N. Primary blast causes mild, moderate, severe and lethal TBI with increasing blast overpressures: experimental rat injury model. Sci Rep. (2016) 6:26992. doi: 10.1038/srep26992

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Pellow S, Chopin P, File SE, Briley M. Validation of open:closed arm entries in an elevated plus-maze as a measure of anxiety in the rat. J Neurosci Methods. (1985) 14:149–67. doi: 10.1016/0165-0270(85)90031-7

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Porsolt RD, Le Pichon M, Jalfre M. Depression: a new animal model sensitive to antidepressant treatments. Nature. (1977) 266:730–2. doi: 10.1038/266730a0

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Masuo Y, Matsumoto Y, Morita S, Noguchi J. A novel method for counting spontaneous motor activity in the rat. Brain Res Brain Res Protoc. (1997) 1:321–6. doi: 10.1016/S1385-299X(97)00005-6

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Swiergiel AH, Leskov IL, Dunn AJ. Effects of chronic and acute stressors and CRF on depression-like behavior in mice. Behav Brain Res. (2008) 186:32–40. doi: 10.1016/j.bbr.2007.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Miyata S, Kumagaya R, Kakizaki T, Fujihara K, Wakamatsu K, Yanagawa Y. Loss of glutamate decarboxylase 67 in somatostatin-expressing neurons leads to anxiety-like behavior and alteration in the Akt/GSK3β signaling pathway. Front Behav Neurosci. (2019) 13:131. doi: 10.3389/fnbeh.2019.00131

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Gouw AA, Seewann A, Vrenken H, van der Flier WM, Rozemuller JM, Barkhof F, et al. Heterogeneity of white matter hyperintensities in Alzheimer's disease: post-mortem quantitative MRI and neuropathology. Brain. (2008) 131(Pt 12):3286–98. doi: 10.1093/brain/awn265

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Faul F, Erdfelder E, Lang AG, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. (2007) 39:175–91. doi: 10.3758/BF03193146

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Schmitt U, Hiemke C. Strain differences in open-field and elevated plus-maze behavior of rats without and with pretest handling. Pharmacol Biochem Behav. (1998) 59:807–11. doi: 10.1016/S0091-3057(97)00502-9

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Yamada K, Kobayashi M, Mori A, Jenner P, Kanda T. Antidepressant-like activity of the adenosine A(2A) receptor antagonist, istradefylline (KW-6002), in the forced swim test and the tail suspension test in rodents. Pharmacol Biochem Behav. (2013) 114-115:23–30. doi: 10.1016/j.pbb.2013.10.022

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kovesdi E, Kamnaksh A, Wingo D, Ahmed F, Grunberg NE, Long JB, et al. Acute minocycline treatment mitigates the symptoms of mild blast-induced traumatic brain injury. Front Neurol. (2012) 3:111. doi: 10.3389/fneur.2012.00111

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Budde MD, Shah A, McCrea M, Cullinan WE, Pintar FA, Stemper BD. Primary blast traumatic brain injury in the rat: relating diffusion tensor imaging and behavior. Front Neurol. (2013) 4:154. doi: 10.3389/fneur.2013.00154

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Elder GA, Dorr NP, De Gasperi R, Gama Sosa MA, Shaughness MC, Maudlin-Jeronimo E, et al. Blast exposure induces post-traumatic stress disorder-related traits in a rat model of mild traumatic brain injury. J Neurotrauma. (2012) 29:2564–75. doi: 10.1089/neu.2012.2510

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Perez-Garcia G, Gama Sosa MA, De Gasperi R, Tschiffely AE, McCarron RM, Hof PR, et al. Blast-induced “PTSD”: evidence from an animal model. Neuropharmacology. (2019) 145(Pt B):220–9. doi: 10.1016/j.neuropharm.2018.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Blaze J, Choi I, Wang Z, Umali M, Mendelev N, Tschiffely AE, et al. Blast-related mild TBI alters anxiety-like behavior and transcriptional signatures in the rat amygdala. Front Behav Neurosci. (2020) 14:160. doi: 10.3389/fnbeh.2020.00160

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Kawa L, Arborelius UP, Yoshitake T, Kehr J, Hökfelt T, Risling M, et al. Neurotransmitter systems in a mild blast traumatic brain injury model: catecholamines and serotonin. J Neurotrauma. (2015) 32:1190–9. doi: 10.1089/neu.2014.3669

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Goldstein LE, Fisher AM, Tagge CA, Zhang XL, Velisek L, Sullivan JA, et al. Chronic traumatic encephalopathy in blast-exposed military veterans and a blast neurotrauma mouse model. Sci Transl Med. (2012) 4:134ra60. doi: 10.1016/j.jalz.2012.05.592

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Bubb EJ, Metzler-Baddeley C, Aggleton JP. The cingulum bundle: anatomy, function, and dysfunction. Neurosci Biobehav Rev. (2018) 92:104–27. doi: 10.1016/j.neubiorev.2018.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Banks SJ, Eddy KT, Angstadt M, Nathan PJ, Phan KL. Amygdala-frontal connectivity during emotion regulation. Soc Cogn Affect Neurosci. (2007) 2:303–12. doi: 10.1093/scan/nsm029

CrossRef Full Text | Google Scholar

44. Wedeen VJ, Rosene DL, Wang R, Dai G, Mortazavi F, Hagmann P, et al. The geometric structure of the brain fiber pathways. Science. (2012) 335:1628–34. doi: 10.1126/science.1215280

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Goeller J, Wardlaw A, Treichler D, O'Bruba J, Weiss G. Investigation of cavitation as a possible damage mechanism in blast-induced traumatic brain injury. J Neurotrauma. (2012) 29:1970–81. doi: 10.1089/neu.2011.2224

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Ando T, Sato S, Takano S, Ashida H, Obara M. Propagation characteristics of laser-induced stress wave in deep tissue for gene transfer. Appl Surf Sci. (2009) 255:9898–901. doi: 10.1016/j.apsusc.2009.04.115

CrossRef Full Text | Google Scholar

47. Shively SB, Horkayne-Szakaly I, Jones RV, Kelly JP, Armstrong RC, Perl DP. Characterisation of interface astroglial scarring in the human brain after blast exposure: a post-mortem case series. Lancet Neurol. (2016) 15:944–53. doi: 10.1016/S1474-4422(16)30057-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: diffuse axonal injury, elevated plus maze test, forced swimming test, spontaneous motor activity, blast-induced traumatic brain injury, laser-induced shock wave

Citation: Jitsu M, Niwa K, Suzuki G, Obara T, Iwama Y, Hagisawa K, Takahashi Y, Matsushita Y, Takeuchi S, Nawashiro H, Sato S and Kawauchi S (2021) Behavioral and Histopathological Impairments Caused by Topical Exposure of the Rat Brain to Mild-Impulse Laser-Induced Shock Waves: Impulse Dependency. Front. Neurol. 12:621546. doi: 10.3389/fneur.2021.621546

Received: 26 October 2020; Accepted: 23 April 2021;
Published: 21 May 2021.

Edited by:

Randall Reed Benson, Center for Neurological Studies (CNS), United States

Reviewed by:

Stuart J. McDonald, Monash University, Australia
Matthew J. Robson, University of Cincinnati, United States
Michael Shaughness, Uniformed Services University of the Health Sciences, United States

Copyright © 2021 Jitsu, Niwa, Suzuki, Obara, Iwama, Hagisawa, Takahashi, Matsushita, Takeuchi, Nawashiro, Sato and Kawauchi. 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: Satoko Kawauchi, satok-bits@ndmc.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.