- 1Department of Regulation Biochemistry, Graduate School of Medical Sciences, Kitasato University, Sagamihara, Japan
- 2Department of Behavioral Medicine, Tohoku University School of Medicine, Sendai, Japan
Leucine-rich repeat kinase 2 (LRRK2) is a molecule associated with familial and sporadic Parkinson’s disease. It regulates many central neuronal functions, such as cell proliferation, apoptosis, autophagy, and axonal extension. Recently, it has been revealed that LRRK2 is related to anxiety/depression-like behavior, implying an association between LRRK2 and stress. In the present study, we investigated for the first time the stress pathway and its relationship to gastrointestinal motility in LRRK2-knockout (KO) mice. The mice were subjected to acute restraint stress, and analyzed for activation of the paraventricular nucleus of the hypothalamus (PVN) using an immunohistochemical approach. Phosphorylation of extracellular signal-regulated kinase 1/2 (ERK1/2) was assessed by Western blotting. The KO mice showed a lower number of c-Fos-positive cells and disruption of the ERK signaling pathway in the PVN in the presence of restraint stress. Stress responses in terms of both upper and lower gastrointestinal motility were alleviated in the mice, accompanied by lower c-Fos immunoreactivity in enteric excitatory neurons. Our present findings suggest that LRRK2 is a newly recognized molecule regulating the stress pathway in the PVN, playing a role in stress-related gastrointestinal dysmotility.
Introduction
Leucine-rich repeat kinase 2 (LRRK2) was originally defined by linkage analysis of familial Parkinson’s disease (PD), and its association with inflammatory bowel disease (IBD) and leprosy has been demonstrated recently (Funayama et al., 2002; Paisán-Ruíz et al., 2004; Barrett et al., 2008; Zhang et al., 2009). LRRK2 is a complex kinase comprising Leucine-rich repeat (LRR), Ras of complex (ROC), C-terminal of Roc (COR), kinase, and WD40 domains (Meylan and Tschopp, 2005). It has already been reported that LRRK2 is expressed in neurons, astroglia, and microglia in the central nervous system (CNS), and several types of immune cells (Higashi et al., 2007; Kubo et al., 2010; Liu et al., 2011; Moehle et al., 2012). Accumulated evidence suggests that LRRK2 plays a key role in endocytosis, axonal extension, autophagy, proliferation, and survival of neurons (Esteves et al., 2014; Maekawa et al., 2016a, b). While many studies have been focusing on the role of LRRK2 in the pathogenesis of motor dysfunction in PD, some groups have reported an association between LRRK2 and the subtle non-motor phenotypes of PD. It has been shown that LRRK2-knockout (KO) mice display anxiety/depression-like behavior (Hinkle et al., 2012). In addition, transgenic mice harboring the G2019S (glycine to serine substitution at amino acid 2019) mutation have been shown to display anxiety/depression-like behavior, despite the paradoxical phenotype of G2019S mutation gain of function and loss of LRRK2 in KO mice (Lim et al., 2018). These findings suggest an association between LRRK2 and stress.
Physiological or psychological stress activates several brain regions, including the amygdala, paraventricular nucleus (PVN), and bed nucleus of the stria terminals, followed by activation of the pituitary gland (Rutherford et al., 2011). This activation induces peripheral stress responses, such as increases in gluconeogenesis, vasoconstriction, and heart rate together with decreases in growth, immunity, and digestion via the sympathetic system and hormonal pathway (Panagiotakopoulos and Neigh, 2014). Restraint stress in rodents is a well-known model of typical stress that activates PVN followed by development of peripheral responses (Kwon et al., 2006; Zheng et al., 2009; Campos-Rodríguez et al., 2013).
Gastrointestinal dysmotility is a feature of PD as well as stress and psychiatric disease (North et al., 2007; Konturek et al., 2011; Pellegrini et al., 2015). While many clinical studies have been reported, few investigations have focused on the role of enteric neurons and brain-gut interaction in PD-related gastrointestinal dysmotility (Derkinderen et al., 2011). Furthermore, in contrast to accumulated data on the function of LRRK2 in the CNS, it is unclear whether LRRK2 is associated with regulation of the stress signaling pathway and gastrointestinal motility.
In the present study using LRRK2-KO mice, we have revealed for the first time that the extracellular signal-regulated kinase (ERK) signaling pathway is disrupted in the PVN and that stress suppresses gastrointestinal dysmotility, thus highlighting the role of LRRK2 in the regulation of stress signaling and gastrointestinal motility.
Materials and Methods
Animals
Male LRRK2-KO mice on C57BL/6J, developed by Hinkle et al. (2012), and littermate wild-type (WT) mice aged 8–14 weeks were used. The mice were housed in a light- and temperature-controlled room with water and food available ad libitum. For sacrifice, the mice were euthanatized by cervical dislocation. All procedures had been approved by the Animal Experimentation and Ethics Committee of Kitasato University.
Restraint Stress
Single exposure to restraint stress for 60 min was performed by placing the animal in a 2.9 × 11.5-cm bottle, adjusting it with tape on the outside so that the animal was unable to move. 1 h after the end of stress exposure, the mice were sacrificed for each type of analysis.
Immunostaining of Brain Sections
Mice were subjected to flush-perfusion with phosphate-buffered saline (PBS) followed by perfusion-fixation with 4% paraformaldehyde (PFA). The removed brain was immersed in 4% PFA overnight, and subsequently in 30% sucrose for 48 h at 4°C. After the brain had been cut into coronal 30-μm sections, H2O2-inactivation of endogenous peroxidase activity was performed for immunohistochemistry, followed by treatment with 2% bovine serum albumin (BSA) in PBS-0.2% Triton X-100 for 60 min at room temperature to block non-specific protein binding. The sections were incubated with a rabbit monoclonal antibody against c-Fos (Cell Signaling Technology, Beverly, MA, United States) and a mouse monoclonal antibody against FOX3 overnight at 4°C. After incubation with a HRP-conjugated secondary antibody, the sections were then treated with 3,3′-diaminobenzidine (DAB) for immunohistochemistry. An Alexa Fluor 594-labeled donkey polyclonal antibody against rabbit IgG (H + L), F (ab’)2 fragment and an Alexa Fluor 488-labeled donkey polyclonal antibody against goat IgG F (ab’)2 fragment were used for immunofluorescence. Fluorescence was observed using a Nikon C2 Si confocal microscope system.
Western Blotting
Frozen brain sections (30 μm) were stained with toluidine blue, and then the PVN was collected by laser microdissection (LMD) in accordance with the mouse brain map. The collected PVNs were lysed in 2 x SDS sample buffer (EzApply, ATTO, Tokyo, Japan) and boiled at 100°C for 5 min. The boiled lysate was then subjected to sodium dodecyl sulfonyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) using a 5–20% e-PAGEL gradient gel (ATTO, Tokyo, Japan), and blotted onto polyvinylidene fluoride (PDNF) membranes. The membranes were blocked with 2% skim milk in PBS−0.1% Tween 20 for 60 min at room temperature, and then incubated with a rabbit monoclonal antibody against LRRK2 (MJFF2, Abcam, Cambridge, MA, United States), a rabbit monoclonal antibody against p44/42 MAPK (Erk1/2, Cell Signaling Technology, Beverly, MA, United States), a rabbit monoclonal antibody against phospho-p44/42 MAPK (Thr202/Tyr204, Cell Signaling Technology, Beverly, MA, United States), and a rabbit monoclonal antibody against glyceraldehyde-3-phosphate dehydrogenase (GAPDH, Cell Signaling Technology, Beverly, MA, United States) overnight at 4°C. After incubation with the secondary antibody (HRP-conjugated donkey anti-rabbit IgG antibody) for 30 min at room temperature, protein bands were detected using a Luminata Crescendo (Merck Millipore, Guyancourt, France). The signals were captured by an Odyssey Fc Imaging System (LI-COR, Lincoln, NE, United States) followed by quantification of the signals with Image Studio 5.0.
Quantitative Reverse Transcription (RT)-PCR
The PVN collected by LMD was homogenized in TRIzol Reagent (Thermo Fisher Scientific, Waltham, MA, United States). RNA was isolated in accordance with the manufacturer’s instructions. cDNA synthesis was performed using a ThermoScript RT-PCR System (Thermo Fisher Scientific, Waltham, MA, United States). Real-time PCR was performed using SYBR Green PCR Master Mix and a 7500 Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, United States). PCR primers used were as follows: CRF forward 5′−TGATCCGCATGGGTGAAGAAT ACTTCCTC−3′ and reverse 5′−CCCGATAATCTCCATCAGT TTCCTGTTGCTG−3′, and GAPDH forward 5′−GAGGCCG GTGCTGAGTATGTCGTG−3′ and reverse 5′−TCGGCAGAA GGGGCGGAGAT−3′. The threshold cycle (Ct) value was normalized by reference to GAPDH.
Measurement of Adrenocorticotropic Hormone
Blood samples were collected into tubes with EDTA 1 h after the stress exposure, and plasma samples were obtained by centrifugation at 1,500 rpm for 10 min at 4°C. The amounts of plasma ACTH were measured using an ACTH Rat/Mouse ELISA Kit (Phoenix Pharmaceuticals, Burlingame, CA, United States) in accordance with the manufacturer’s protocol.
Measurement of Metabolism
Body weight and metabolism were measured daily using a metabolic cage for 3 days in a row following habituation to the cage for 7 days. The amounts of food intake, feces, water intake, and urine volume were measured as weight in grams.
Intestinal Transit Assay
Intestinal motility in mice was evaluated from the distribution of orally administered fluorescein isothiocyanate (FITC)-conjugated dextran. Mice were exposed to the restraint stress just before oral administration of the dextran. After transportation of the dextran along the gastrointestinal tract for 3 h, the whole gut was collected after euthanasia. The small intestine and the colon were equally cut into ten segments and four segments, respectively. The fluorescence intensity of each piece was measured using FLUOstar OPTIMA (BMG LABTECH, Saitama, Japan). The intestinal geometric center was used for evaluation of the distribution of dextran, which was calculated using the following equation:
Bead Expulsion Assay
Colonic motility was evaluated by the bead expulsion assay. Following anesthesia with isoflurane (Wako, Osaka, Japan), a zirconia bead (3-mm diameter) was inserted into the rectum at a distance of 2-cm from the anal verge. The time required for expulsion of the bead beyond the anal verge was then measured. Each mouse was subjected to the same experiment again after the restraint stress. For comparison between WT and KO mice, we calculated the reduction in time until excretion of the bead after the restraint stress.
Whole-Mount Staining of the Longitudinal Muscle-Myenteric Plexus
The mouse colon was opened and pinned on a silicone-sheeted glass dish. Each preparation was fixed overnight with 4% PFA, and then the mucosa, submucosa, and circular muscle layer were removed under a dissecting microscope. Preparations were incubated in PBS-1% Triton X-100 for 30 min, followed by incubation with rabbit monoclonal antibody against c-Fos (Cell Signaling Technology, Beverly, MA, United States), goat polyclonal antibody against choline acetyltransferase (ChAT, Merck Millipore, Guyancourt, France), and goat polyclonal antibody against neuronal nitric oxide synthase (nNOS, Abcam, Cambridge, MA, United States) for 24 h at 4 C. And subsequently with Alexa Fluor 594-labeled donkey polyclonal antibody against rabbit IgG (H + L), F (ab’)2 fragment (Cell Signaling Technology, Beverly, MA, United States) or Alexa Fluor 488-labeled donkey polyclonal antibody against goat IgG F (ab’)2 fragment (Abcam, Cambridge, MA, United States) for 3 h at 4°C. Nuclear staining was performed with 2′-(4-ethoxyphenyl)-5-(4-methyl-1-piperazinyl)-2,5′-bi-1H-benzimidazole (Hoechst 33342, Kumamoto, Japan). Fluorescence was observed using a Nikon C2 Si confocal microscope system.
Statistical Analysis
All data are expressed as means + SEM. Significance of differences was assessed by Student’s t test for two comparisons and two-way ANOVA for multiple comparisons. Differences at p < 0.05 were considered to be significant. Sample size n represents the number of animals.
Results
Attenuated Reactivity of the PVN in Response to Restraint Stress in LRRK2-KO Mice
The PVN of the hypothalamus is one of the major regions involved in reactivity to stress, and restraint is known to be a major stressor that activates neurons (Kwon et al., 2006; Zheng et al., 2009; Rutherford et al., 2011; Campos-Rodríguez et al., 2013). Therefore, to explore the association between LRRK2 and stress, expression of c-Fos (a marker of neuronal activation) under restraint stress was analyzed in LRRK2-KO mice. Activated PVN neurons were clearly visualized by immunostaining with an antibody against c-Fos only after restraint stress in both WT and KO mice (Figures 1A,B). Interestingly, the number of c-Fos positive neurons was strikingly decreased in KO mice relative to WT mice, indicating attenuation of neuronal reactivity to the stress in KO mice (Figure 1C). We confirmed that reactivity for c-Fos after the stress was detected only in neurons under our restraint stress conditions (Figure 2).
Figure 1. Immunohistochemistry using an antibody against c-Fos in the PVN. (A) Section representing a level –1.94 mm posterior to the bregma. Yellow area indicates the PVN of the hypothalamus. (B) Representative immunohistochemistry using an antibody against c-Fos in the PVN of WT and LRRK2-KO mice in the absence or presence of restraint stress (RS). Methyl green was used for counter staining. Scale bar = 100 μm. (C) The total number of c-Fos-positive cells counted in the PVN. n = 5 per group. In all graphical representations, data are expressed as mean + SEM and were assessed by Student’st test (WT vs. KO), ∗p < 0.05.
Figure 2. Immunofluorescence of c-Fos positive cells in the PVN. Confocal images of double immunostaining using antibodies against NeuN and c-Fos in the PVN of WT and LRRK2-KO mice in the absence or presence of restraint stress (RS). Scale bar = 100 μm.
Disruption of the ERK Signaling Pathway in the PVN of LRRK2-KO Mice
It has been suggested that ERK signaling is the pathway affected by LRRK2 through its kinase activity, and that it is also crucial for stress signaling under stress (Monje et al., 2005; Liou et al., 2008). To investigate whether the ratio of phosphorylated ERK1/2 is changed in LRRK2-KO mice, the protein levels of LRRK2, pERK1/2, ERK1/2, and GAPDH in the PVN were quantified by Western blotting. A PVN-protein lysate was extracted from pieces of tissue collected by LMD (Figure 3A). Western blotting analysis detected specific bands at each of the predicted molecular weights in both WT and KO mice, and equivalent protein loading was confirmed by comparison with the signal intensity of GAPDH (Figure 3B). Unchanged levels of LRRK2 and ERK1/2 total protein after the stress were confirmed by GAPDH normalization (Figure 3C). It is noteworthy that the ratio of ERK1/2 phosphorylation was significantly increased in LRRK2-KO mice even in the absence of stress. Moreover, phosphorylation of ERK1/2 was increased by the stress in WT mice, but not in KO mice. These results indicated that LRRK2 deficiency dramatically affects the phosphorylation of ERK1/2 in the PVN, suggesting that LRRK2 acts as a molecule in stress pathways. Corticotropin releasing hormone (CRF) is the most important transmitter involved in the regulation of stress responses within the PVN (Lightman, 2008). Therefore, we investigated whether CRF expression was increased following restraint stress. Samples dissected from the PVN were subjected to quantitative RT-PCR using a CRF primer and a GAPDH primer as a target and internal control, respectively. Restraint stress induced up-regulation of the CRF-mRNA level in WT mice following the stress (Figure 4). On the other hand, the mRNA level in KO mice was not up-regulated following the stress.
Figure 3. Analysis of the PVN signaling pathway by Western blotting. (A) Representative toluidine blue staining at the level of the posterior PVN before (upper panel) and after (lower panel) LMD. Scale bar = 100 μm. (B) Representative images of Western blotting using antibodies against LRRK2, pERK1/2 (Thr202/Tyr204), ERK1/2, and GAPDH in LRRK2-KO mice. (C) Quantitative analysis of the density of protein bands. Each signal intensity was normalized against that of GAPDH or pan-ERK1/2. n = 5 per group. In all graphical representations, data are expressed as mean + SEM and were assessed by two-way ANOVA, ∗p < 0.05.
Figure 4. Down-regulated CRF-mRNA level in the PVN in the presence of restraint stress in LRRK2-KO mice. Total mRNA prepared from microdissected PVN in the absence or presence of restraint stress was subjected to quantitative PCR using a specific primer for CRF. n = 6 per group. In all graphical representations, data are expressed as mean + SEM and were assessed by two-way ANOVA, ∗p < 0.05.
Down-Regulation of ACTH Reactivity in LRRK2-KO Mice
An increase in the plasma level of ACTH is a major response following PVN activation caused by stress. Quantification of plasma ACTH revealed small amounts of ACTH in the absence of stress, without any significant difference between WT and KO mice (Figure 5). This increase in the level of ACTH was detected in the presence of stress in both WT and KO mice, but was significantly lower in the latter than in the former, in line with decreased activation of the PVN.
Figure 5. Lower plasma level of ACTH in the presence of restraint stress in LRRK2-KO mice. Quantitative analysis of plasma ACTH by ELISA in the absence or presence of restraint stress. n = 6 per group. In all graphical representations, data are expressed as mean + SEM and were assessed by two-way ANOVA, ∗∗p < 0.01, ∗p < 0.05.
Aberrant Gastrointestinal Response to Restraint Stress in LRRK2-KO Mice
It is known that gastrointestinal motility is affected by stress, sometimes resulting in functional gastrointestinal disorder (FGID). To investigate gastrointestinal motility under restraint stress, we performed the intestinal transit assay using FITC-conjugated dextran as an indicator of upper gastrointestinal motility and the bead expulsion assay as an indicator of colonic motility. Before these assays, the amounts of food intake, feces production, water intake, and urine volume per day were measured using metabolic cages. Accompanied by normal weight fluctuation, the metabolism of LRRK2-KO mice was comparable to that of WT mice (Supplementary Figure S1). The intestinal transit assay demonstrated a gradient distribution throughout the gastrointestinal tract with a peak point (Figure 6A). The peak of dextran distribution shifted to the oral side in the presence of stress relative to that in its absence, indicating delayed upper gastrointestinal motility in WT mice. The geometric center of dextran distribution also showed that upper gastrointestinal motility was delayed under restraint stress in WT mice (Figure 6B). In contrast, no alteration in the distribution of dextran was observed in KO mice, as indicated by an unchanged geometric center in the presence of stress. The bead expulsion assay for analysis of colonic motility revealed that the time required to excrete a bead that had been inserted into the rectum was shorter in the presence of stress than its absence in WT mice, indicating that the stress had induced colonic dysmotility (Figure 6C). This stress-induced colonic dysmotility was not observed in KO mice, suggesting disruption of stress reactivity in terms of colonic as well as upper gastrointestinal motility in KO mice.
Figure 6. Aberrant gastrointestinal response to restraint stress in LRRK2-KO mice. (A) Intestinal transit assay using FITC-conjected dextran in the absence or presence of restraint stress. Histogram represents the proximal to distal distribution of orally administered FITC-conjugated dextran along the gastrointestinal tract (1–14). n = 6 per group. The data represent the percentage distribution relative to the whole amount of dextran. (B) Intestinal geometric center of the distribution of dextran. (C) Bead expulsion assay for analysis of colonic motility. The time required to expel the bead was measured. Data represent the reduction in the time period until excretion due to stress. n = 6 per group. In all graphical representations, data are expressed as mean + SEM and were assessed by two-way ANOVA for multiple comparison and Student’s t test for two comparison (WT vs. KO), ∗p < 0.05.
Attenuation of c-Fos Immunoreactivity of Enteric Excitatory Neurons in Response to Restraint Stress in LRRK2-KO Mice
Since gastrointestinal motility is orchestrated by enteric neuronal circuits concurrently with signals from the CNS, enteric neurons activated in response to restraint stress were characterized by immunostaining for LMMP using antibodies against ChAT or nNOS, which visualize enteric excitatory or inhibitory neurons, respectively, and c-Fos. Activated neurons were simultaneously visualized using the antibody against c-Fos. This revealed that excitatory neurons were predominant in enteric neuronal ganglia, whereas inhibitory neurons were sparse (Figures 7A,B). Immunoreactivity for c-Fos was evident in the nuclei of neuronal cells only after the stress, and c-Fos positivity was more abundant in ChAT-positive than in nNOS-positive cells (Figure 7C). In KO mice, double-immunopositive cells (ChAT+/c-Fos+) were significantly decreased in comparison with WT mice; however, there was no difference in the number of nNOS+/c-Fos+ cells between WT and KO mice. These results suggested that disruption of gastrointestinal dysmotility in LRRK2-KO mice is due to lower activation of enteric excitatory neurons.
Figure 7. Attenuated reactivity of enteric excitatory neurons in response to restraint stress in LRRK2-KO mice. (A) Confocal images of double immunostaining using antibodies against ChAT and c-Fos in the LMMP of the colon after restraint stress. Nuclei were visualized using Hoechst 33342. Scale bar = 100 μm. (B) Confocal images of double immunostaining using antibodies against nNOS and c-Fos in the LMMP of the colon after restraint stress. (C) Graphical data represent the number of double-immunopositive cells (ChAT+/c-Fos+ or nNOS+/c-Fos+) per area in the colon. n = 6 per group in the presence of restraint stress. In all graphical representations, data are expressed as mean + SEM and were assessed by Student’s t test (WT vs. KO), ∗∗p < 0.01.
Discussion
In this study, we investigated the role of LRRK2 in the PVN stress pathway and the resulting gastrointestinal dysmotility for the first time. LRRK2-KO mice exhibited alleviation of gastrointestinal dysmotility caused by restraint stress, along with disruption of the ERK signaling pathway. Furthermore, peripheral stress responses – in terms of an increased plasma ACTH level and gastrointestinal dysmotility – were suppressed in KO mice. Our study has highlighted LRRK2 as a new candidate molecule regulating PVN reactivity in stress and gastrointestinal dysmotility.
It has been shown that the PVN plays a crucial role in the brain as a stress sensor (Herman, 2012). Given that LRRK2-KO mice exhibited a lower number of c-Fos-positive neurons in the presence of restraint stress than was the case in WT mice, there is a possibility that LRRK2 up-regulates neuronal reactivity in the PVN in response to stress (Figure 1C). Transgenic mice harboring LRRK2 mutation (R1441G, R1441C, G2019S, I2020T) might be helpful for understanding the association between stress and PD. We further revealed that abundant CRF expression following the stress was not observed in KO mice (Figure 3). Since CRF is the important transmitter involved in stress responses within the PVN, we considered that the loss of LRRK2 signaling ensures disturbance of stress responses (Wamsteeker et al., 2013; Jiang et al., 2018). Interestingly, linkage disequilibrium analysis has revealed that CRF receptor 1 is a molecule associated with PD (Oliveira et al., 2004). Accordingly, there is a possibility that the CRF-related pathway is associated with not only stress signaling but also PD development in tandem with LRRK2. From another viewpoint, the lower neuronal reactivity in the PVN of LRRK2-KO mice might be due to defective release and/or binding of CRF, as LRRK2 has a role in vesicle trafficking and receptor trafficking (Dodson et al., 2012; Matta et al., 2012; MacLeod et al., 2013; Beilina et al., 2014; Cho et al., 2014; Yun et al., 2015; Steger et al., 2016; Rassu et al., 2017). In addition to neurons expressing CRF, the PVN contains some types of neuron expressing thyrotropin-releasing hormone, arginine vasopressin, oxytocin, or somatostatin (Biag et al., 2012). It is problematic to clarify the expression of LRRK2 in the PVN because the reactivity of anti-LRRK2 antibody is limited, particularly under the conditions used for immunostaining (Davies et al., 2013). Only mRNA data are available with regard to LRRK2 expression in the PVN in rat and mouse brain (Simón-Sánchez et al., 2006; Taymans et al., 2006). Genetic approaches, such as the use of LRRK2-reporter mice, could be one way to resolve this issue.
In terms of c-Fos immunoreactivity, astroglia also have positive immunoreactivity under stress conditions (Fuente-Martín et al., 2016). Although the c-Fos-positive cells were neurons in our stress model, astroglia are another type of regulatory cell that can control stress reactivity through PVN neurons.
A number of intercellular signaling events have been attributed to LRRK2 kinase activity, and in addition LRRK2 acts as a scaffolding protein by interacting with ASK1, MKK3/6, and p38MAPK (Lobbestael et al., 2012; Yoon et al., 2017). Since the level of phosphorylation at Ser1292 has been proposed as a direct marker of LRRK2 kinase activity, we tried to detect phosphorylated LRRK2 (Ser1292) using PVN samples but were unable to find a clear specific band. Our present results were unable to indicate whether kinase activity or scaffolding function is substantially affected. Further in vitro experimental data will be necessary to clarify the predominant role of LRRK2 in the stress-signaling pathway.
LRRK2-KO mice also showed lower reactivity to stress and higher expression of phosphorylated ERK1/2 in the PVN, which appeared to suppress any stress-induced effects. As protein kinase A (PKA) activation is required for phosphorylation of ERK1/2 and LRRK2 negatively regulates PKA, augmented phosphorylation of ERK1/2 could be due to higher PKA activity caused by LRRK2 depletion (Impey et al., 1998; Zanassi et al., 2001; Chen et al., 2010; Parisiadou et al., 2014; Greggio et al., 2017; Russo et al., 2018). Kubo et al. (2016) have reported a higher degree of ERK1/2 phosphorylation in B2 B-cells of LRRK2-KO mice, thus highlighting the importance of further investigation into the interaction between LRRK2 and ERK1/2 (Kubo et al., 2016).
The neuronal activity of the PVN is regulated via the negative feedback action of corticosterone. Osterlund et al. have reported that ERK1/2 activation was inhibited by the tonic but not phasic feedback action of corticosterone (Osterlund et al., 2011). They suggested that ERK1/2 activation of the PVN may be associated with altered activity of stress-dependent neuronal inputs to the PVN, and that alteration of stress-dependent intracellular signaling mechanisms within those neurons might be associated with ERK1/2 activation. Therefore, this feedback mechanism might induce persistent activation of ERK1/2 and defective ERK1/2 activation in KO mice.
Although restraint stress suppresses upper gastrointestinal motility and promotes colonic motility, the response of LRRK2-KO mice to this form of stress was deficient and atypical (Taché et al., 1999). Although the lower ACTH level in KO mice indicates disruption of the hypothalamic-pituitary-adrenal (HPA) axis, other pathways can transmit stress signals from the brain to peripheral organs. The autonomic nervous pathway, prolactin release, gonadal steroid release, and IL-6 release may also have possible roles in regulation of the stress reaction (Herman et al., 2016). As the HPA axis plays an important role in stress responses, it may be necessary to consider several aspects of the stress reactive pathway, in order to integrally understand its association with LRRK2.
Activation of cholinergic enteric neurons in the presence of stress was suppressed in KO mice relative to WT mice, suggesting that disruption of the stress response in KO mice, in terms of altered gastrointestinal motility, was due to reduced reactivity of cholinergic neurons. Recently, a PET study has shown that the abnormal cholinergic function of individuals with LRRK2 mutations, both with or without PD, was due to increased rates of acetylcholinesterase hydrolysis (Liu et al., 2018). In addition, Higashi et al. (2007) have revealed that LRRK2 is expressed in ChAT-positive, but not nNOS-positive neurons (Higashi et al., 2007). These reports indicate that LRRK2 shows higher expression in cholinergic than in nitrergic neurons, and that the cholinergic system might be more affected by LRRK2 aberration. Our recent study showed that, like brain neurons, enteric neurons also express LRRK2 (Maekawa et al., 2016a, b). However, the characteristics of these neurons have not been revealed. Identification of the cells predominantly expressing LRRK2 in the enteric nervous system (ENS) could be the key to elucidating the mechanism of dysmotility in LRRK2-KO mice. More specific characterization of LRRK2-expressing cells in the ENS would help to shed light on cholinergic dysfunction in PD.
Enteric glial cells (EGCs) have emerged as another component of gastrointestinal motility control (Neunlist et al., 2014). Fujikawa et al. (2015) have investigated the association between EGCs and gastrointestinal dysmotility in maternally separated rats (Fujikawa et al., 2015). Also, a lower level of glial fibrillary acidic protein (GFAP) phosphorylation has been observed in colonic biopsy samples from PD patients (Clairembault et al., 2014). Given that neurons and EGCs comprise a huge complex circuit, and that astroglia express LRRK2 in the CNS, it is important to further analyze the ENS structure of LRRK2-KO mice, including their degree of EGC complexity (Coelho-Aguiar et al., 2015).
Gastrointestinal dysfunction is common in both idiopathic and familial PD (Marras et al., 2011). However, Gaiq et al. (2014) have recently reported that LRRK2 G2019S patients show no symptomatic gastrointestinal dysmotility (Gaiq et al., 2014). Since our study using LRRK2-KO mice also demonstrated normal gastrointestinal motility in the absence of stress, some additional risk factor, for example gene mutation, or environmental and/or psychological stress, might be necessary for induction of gastrointestinal dysfunction in LRRK2-related PD.
Taken together, our present findings suggest that LRRK2 is a newly recognized molecule regulating the stress pathway in the PVN, thus playing a role in stress-related gastrointestinal dysmotility. These findings will be of help for further studies of LRRK2 in the context of both the gut and the brain, especially in relation to PD.
Data Availability
All datasets generated for this study are included in the manuscript and/or the Supplementary Files.
Ethics Statement
All procedures had been approved by the Animal Experimentation and Ethics Committee of Kitasato University.
Author Contributions
TM, HT, and MI performed the experiments with support from FK, RK, MK, and TI. TM designed the study and wrote the manuscript with support from TI. All authors have read and approved the manuscript.
Funding
This study was supported by grant in Kitasato University School of Allied Health Sciences (Grant-in-Aid for Research Project 2015–2017).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Ms. Yamane and Ms. Sato for technical assistance in the analysis of PVN.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2019.00905/full#supplementary-material
References
Barrett, J. C., Hansoul, S., Nicolae, D. L., Cho, J. H., Duerr, R. H., Rioux, J. D., et al. (2008). Genome-wide association defines more than 30 distinct susceptibility loci for Crohn’s disease. Nat. Genet. 40, 955–962. doi: 10.1038/ng.175
Beilina, A., Rudenko, I. N., Kaganovich, A., Civiero, L., Chau, H., Kalia, S. K., et al. (2014). Unbiased screen for interactors of leucine-rich repeat kinase 2 supports a common pathway for sporadic and familial Parkinson disease. Proc. Natl. Acad. Sci. U.S.A. 111, 2626–2631. doi: 10.1073/pnas.1318306111
Biag, J., Huang, Y., Gou, L., Hintiryan, H., Askarinam, A., Hahn, J. D., et al. (2012). Cyto- and chemoarchitecture of the hypothalamic paraventricular nucleus in the C57BL/6J male mouse: a study of immunostaining and multiple fluorescent tract tracing. J. Comp. Neurol. 520, 6–33. doi: 10.1002/cne.22698
Campos-Rodríguez, R., Godínez-Victoria, M., Abarca-Rojano, E., Pacheco-Yépez, J., Reyna-Garfias, H., Barbosa-Cabrera, R. E., et al. (2013). Stress modulates intestinal secretory immunoglobulin A. Front. Integr. Neurosci. 7:86. doi: 10.3389/fnint.2013.00086
Chen, W. K., Liu, I. Y., Chang, Y. T., Chen, Y. C., Chen, C. C., Yen, C. T., et al. (2010). Ca(v)3.2 T-type Ca2+ channel-dependent activation of ERK in paraventricular thalamus modulates acid-induced chronic muscle pain. J. Neurosci. 30, 10360–10368. doi: 10.1523/JNEUROSCI.1041-10.2010
Cho, H. J., Yu, J., Xie, C., Rudrabhatla, P., Chen, X., Wu, J., et al. (2014). Leucine-rich repeat kinase 2 regulates Sec16A at ER exit sites to allow ER-golgi export. EMBO. J. 33, 2314–2331. doi: 10.15252/embj.201487807
Clairembault, T., Kamphuis, W., Leclair-Visonneau, L., Rolli-Derkinderen, M., Coron, E., Neunlist, M., et al. (2014). Enteric GFAP expression and phosphorylation in Parkinson’s disease. J. Neurochem. 130, 805–815. doi: 10.1111/jnc.12742
Coelho-Aguiar, J. D. M., Bon-Frauches, A. C., Gomes, A. L., Veríssimo, C. P., Aguiar, D. P., Matias, D., et al. (2015). The enteric glia: identity and functions. Glia 63, 921–935. doi: 10.1002/glia.22795
Davies, P., Hinkle, K. M., Sukar, N. N., Sepulveda, B., Mesias, R., Serrano, G., et al. (2013). Comprehensive characterization and optimization of anti-LRRK2 (leucine-rich repeat kinase 2) monoclonal antibodies. Biochem. J. 453, 101–113. doi: 10.1042/BJ20121742
Derkinderen, P., Rouaud, T., Lebouvier, T., Bruley des Varannes, S., Neunlist, M., and De Giorgio, R. (2011). Parkinson disease: the enteric nervous system spills its guts. Neurology 77, 1761–1767. doi: 10.1212/WNL.0b013e318236ef60
Dodson, M. W., Zhang, T., Jiang, C., Chen, S., and Guo, M. (2012). Roles of the drosophila LRRK2 homolog in Rab7-dependent lysosomal positioning. Hum. Mol. Genet. 21, 1350–1363. doi: 10.1093/hmg/ddr573
Esteves, A. R., Swerdlow, R. H., and Cardoso, S. M. (2014). LRRK2, a puzzling protein: insights into Parkinson’s disease pathogenesis. Exp. Neurol. 261, 206–216. doi: 10.1016/j.expneurol.2014.05.025
Fuente-Martín, E., García-Cáceres, C., Argente-Arizón, P., Díaz, F., Granado, M., Freire-Regatillo, A., et al. (2016). Ghrelin Regulates Glucose and Glutamate Transporters in Hypothalamic Astrocytes. Sci. Rep. 6:23673. doi: 10.1038/srep23673
Fujikawa, Y., Tominaga, K., Tanaka, F., Tanigawa, T., Watanabe, T., Fujiwara, Y., et al. (2015). Enteric glial cells are associated with stress-induced colonic hyper-contraction in maternally separated rats. Neurogastroenterol. Motil. 27, 1010–1023. doi: 10.1111/nmo.12577
Funayama, M., Hasegawa, K., Kowa, H., Saito, M., Tsuji, S., and Obata, F. (2002). A new locus for Parkinson’s disease (PARK8) maps to chromosome 12p11.2-q13.1. Ann. Neurol. 51, 296–301. doi: 10.1002/ana.10113
Gaiq, C., Vilas, D., Infante, J., Sierra, M., García-Gorostiaga, I., Buongiorno, M., et al. (2014). Nonmotor symptoms in LRRK2 G2019S associated Parkinson’s disease. PLoS One 9:e108982. doi: 10.1371/journal.pone.0108982
Greggio, E., Bubacco, L., and Russo, I. (2017). Cross-talk between LRRK2 and PKA: implication for Parkinson’s disease? Biochem. Soc. Trans. 45, 261–267. doi: 10.1042/BST20160396
Herman, J. P. (2012). Neural pathways of stress integration: relevance to alcohol abuse. Alcohol Res. 34, 441–447.
Herman, J. P., McKlveen, J. M., Ghosal, S., Kopp, B., Wulsin, A., Makinson, R., et al. (2016). Regulation of the hypothalamic-pituitary-adrenocortical stress response. Compr. Physiol. 6, 603–621. doi: 10.1002/cphy.c150015
Higashi, S., Moore, D. J., Colebrooke, R. E., Biskup, S., Dawson, V. L., Arai, H., et al. (2007). Expression and localization of Parkinson’s disease-associated leucine-rich repeat kinase 2 in the mouse brain. J. Neurochem. 100, 368–381. doi: 10.1111/j.1471-4159.2006.04246.x
Hinkle, K. M., Yue, M., Behrouz, B., Dächsel, J. C., Lincoln, S. J., Bowles, E. E., et al. (2012). LRRK2 knockout mice have an intact dopaminergic system but display alterations in exploratory and motor co-ordination behaviors. Mol. Neurodegener. 7:25. doi: 10.1186/1750-1326-7-25
Impey, S., Obrietan, K., Wong, S. T., Poser, S., Yano, S., Wayman, G., et al. (1998). Cross talk between ERK and PKA is required for Ca2+ stimulation of CREB-dependent transcription and ERK nuclear translocation. Neuron 21, 869–883. doi: 10.1016/s0896-6273(00)80602-9
Jiang, Z., Rajamanickam, S., and Justice, N. J. (2018). Local Corticotropin-releasing factor signaling in the hypothalamic paraventricular nucleus. J. Neurosci. 38, 1874–1890. doi: 10.1523/JNEUROSCI.1492-17.2017
Konturek, P. C., Brzozowski, T., and Konturek, S. J. (2011). Stress and the gut: pathophysiology, clinical consequences, diagnostic approach and treatment options. J. Physiol. Pharmacol. 62, 591–599.
Kubo, M., Kamiya, Y., Nagashima, R., Maekawa, T., Eshima, K., Azuma, S., et al. (2010). LRRK2 is expressed in B-2 but not in B-1 B cells, and downregulated by cellular activation. J. Neuroimmunol. 229, 123–128. doi: 10.1016/j.jneuroim.2010.07.021
Kubo, M., Nagashima, R., Ohta, E., Maekawa, T., Isobe, Y., Kurihara, M., et al. (2016). Leucine-rich repeat kinase 2 is a regulator of B cell function, affecting homeostasis, BCR sig, IgA production, and TI antigen responses. J. Neuroimmunol. 292, 1–8. doi: 10.1016/j.jneuroim.2016.01.005
Kwon, M. S., Seo, Y. J., Shim, E. J., Choi, S. S., Lee, J. Y., and Suh, H. W. (2006). The effect of single or repeated restraint stress on several signal molecules in paraventricular nucleus, arcuate nucleus and locus coeruleus. Neuroscience 142, 1281–1292. doi: 10.1016/j.neuroscience.2006.07.027
Lightman, S. L. (2008). The neuroendocrinology of stress: a never ending story. J. Neuroendocrinol. 20, 880–884. doi: 10.1111/j.1365-2826.2008.01711.x
Lim, J., Bang, Y., Choi, J. H., Han, A., Kwon, M. S., Liu, K. H., et al. (2018). LRRK2 G2019S Induces anxiety/depression-like behavior before the onset of motor dysfunction with 5-HT1A Receptor Upregulation in Mice. J. Neurosci. 38, 1611–1621. doi: 10.1523/JNEUROSCI.4051-15.2017
Liou, A. K., Leak, R. K., Li, Y., and Zigmond, M. J. (2008). Wild-type LRRK2 but not its mutant attenuates stress-induced cell death via ERK pathway. Neurobiol. Dis. 32, 116–124. doi: 10.1016/j.nbd.2008.06.016
Liu, S. Y., Wile, D. J., Fu, J. F., Valerio, J., Shahinfard, E., McCormick, S., et al. (2018). The effect of LRRK2 mutations on the cholinergic system in manifest and premanifest stages of Parkinson’s disease: a cross-sectional PET study. Lancet Neurol. 17, 309–316. doi: 10.1016/S1474-4422(18)30032-2
Liu, Z., Lee, J., Krummey, S., Lu, W., Cai, H., and Lenardo, M. J. (2011). The kinase LRRK2 is a regulator of the transcription factor NFAT that modulates the severity of inflammatory bowel disease. Nat. Immunol. 12, 1063–1070. doi: 10.1038/ni.2113
Lobbestael, E., Baekelandt, V., and Taymans, J. M. (2012). Phosphorylation of LRRK2: from kinase to substrate. Biochem. Soc. Trans. 40, 1102–1110. doi: 10.1042/BST20120128
MacLeod, D. A., Rhinn, H., Kuwahara, T., Zolin, A., Di Paolo, G., McCabe, B. D., et al. (2013). RAB7L1 interacts with LRRK2 to modify intraneuronal protein sorting and Parkinson’s disease risk. Neuron 77, 425–439. doi: 10.1016/j.neuron.2012.11.033
Maekawa, T., Sasaoka, T., Azuma, S., Ichikawa, T., Melrose, H. L., Farrer, M. J., et al. (2016a). Leucine-rich repeat kinase 2 (LRRK2) regulates α-synuclein clearance in microglia. BMC Neurosci. 17:77. doi: 10.1186/s12868-016-0315-2
Maekawa, T., Shimayama, H., Tsushima, H., Kawakami, F., Kawashima, R., Kubo, M., et al. (2016b). LRRK2: an emerging new molecule in the enteric neuronal system that quantitatively regulates neuronal peptides and iga in the gut. Dig. Dis. Sci. 62, 903–912. doi: 10.1007/s10620-017-4476-4473
Marras, C., Schüle, B., Munhoz, R. P., Rogaeva, E., Langston, J. W., Kasten, M., et al. (2011). Phenotype in parkinsonian and nonparkinsonian LRRK2 G2019S mutation carriers. Neurology 77, 325–333. doi: 10.1212/WNL.0b013e318227042d
Matta, S., Van Kolen, K., da Cunha, R., van den Bogaart, G., Mandemakers, W., Miskiewicz, K., et al. (2012). LRRK2 controls an EndoA phosphorylation cycle in synaptic endocytosis. Neuron 75, 1008–1021. doi: 10.1016/j.neuron.2012.08.022
Meylan, E., and Tschopp, J. (2005). The RIP kinases: crucial integrators of cellular stress. Trends. Biochem. Sci. 30, 151–159. doi: 10.1016/j.tibs.2005.01.003
Moehle, M. S., Webber, P. J., Tse, T., Sukar, N., Standaert, D. G., DeSilva, T. M., et al. (2012). LRRK2 inhibition attenuates microglial inflammatory responses. J. Neurosci. 32, 1602–1611. doi: 10.1523/JNEUROSCI.5601-11.2012
Monje, P., Hernández-Losa, J., Lyons, R. J., Castellone, M. D., and Gutkind, J. S. (2005). Regulation of the transcriptional activity of c-Fos by ERK. a novel role for the prolyl isomerase PIN1. J. Biol. Chem. 280, 35081–35084. doi: 10.1074/jbc.C500353200
Neunlist, M., Rolli-Derkinderen, M., Latorre, R., Van Landeghem, L., Coron, E., Derkinderen, P., et al. (2014). Enteric glial cells: recent developments and future directions. Gastroenterology 147, 1230–1237. doi: 10.1053/j.gastro.2014.09.040
North, C. S., Hong, B. A., and Alpers, D. H. (2007). Relationship of functional gastrointestinal disorders and psychiatric disorders: implications for treatment. World J. Gastroenterol. 13, 2020–2027.
Oliveira, S. A., Scott, W. K., Zhang, F., Stajich, J. M., Fujiwara, K., Hauser, M., et al. (2004). Linkage disequilibrium and haplotype tagging polymorphisms in the Tau H1 haplotype. Neurogenetics 5, 147–155. doi: 10.1007/s10048-004-0180-5
Osterlund, C. D., Jarvis, E., Chadayammuri, A., Unnithan, R., Weiser, M. J., and Spencer, R. L. (2011). Tonic, but not phasic corticosterone, constrains stress activatedextracellular-regulated-kinase 1/2 immunoreactivity within the hypothalamic paraventricular nucleus. J. Neuroendocrinol. 23, 1241–1251. doi: 10.1111/j.1365-2826.2011.02220.x
Paisán-Ruíz, C., Jain, S., Evans, E. W., Gilks, W. P., Simón, J., van der Brug, M., et al. (2004). Cloning of the gene containing mutations that cause PARK8-linked Parkinson’s disease. Neuron 44, 595–600. doi: 10.1016/j.neuron.2004.10.023
Panagiotakopoulos, L., and Neigh, G. N. (2014). Development of the HPA axis: where and when do sex differences manifest? Front. Neuroendocrinol. 35:285–302. doi: 10.1016/j.yfrne.2014.03.002
Parisiadou, L., Yu, J., Sgobio, C., Xie, C., Liu, G., Sun, L., et al. (2014). LRRK2 regulates synaptogenesis and dopamine receptor activation through modulation of PKA activity. Nat. Neurosci. 17, 367–376. doi: 10.1038/nn.3636
Pellegrini, C., Antonioli, L., Colucci, R., Ballabeni, V., Barocelli, E., Bernardini, N., et al. (2015). Gastric motor dysfunctions in Parkinson’s disease: current pre-clinical evidence. Parkinsonism Relat. Disord. 21, 1407–1414. doi: 10.1016/j.parkreldis.2015.10.011
Rassu, M., Del Giudice, M. G., Sanna, S., Taymans, J. M., Morari, M., Brugnoli, A., et al. (2017). Role of LRRK2 in the regulation of dopamine receptor trafficking. PLoS One 12:e0179082. doi: 10.1371/journal.pone.0179082
Russo, I., Di Benedetto, G., Kaganovich, A., Ding, J., Mercatelli, D., Morari, M., et al. (2018). Leucine-rich repeat kinase 2 controls protein kinase a activation state through phosphodiesterase 4. J. Neuroinflamm. 15:297. doi: 10.1186/s12974-018-1337-8
Rutherford, H. J., Williams, S. K., Moy, S., Mayes, L. C., and Johns, J. M. (2011). Disruption of maternal parenting circuitry by addictive process: rewiring of reward and stress systems. Front. Psychiatry 2:37. doi: 10.3389/fpsyt.2011.00037
Simón-Sánchez, J., Herranz-Pérez, V., Olucha-Bordonau, F., and Pérez-Tur, J. (2006). LRRK2 is expressed in areas affected by Parkinson’s disease in the adult mouse brain. Eur. J. Neurosci. 23, 659–666. doi: 10.1111/j.1460-9568.2006.04616.x
Steger, M., Tonelli, F., Ito, G., Davies, P., Trost, M., Vetter, M., et al. (2016). Phosphoproteomics reveals that Parkinson’s disease kinase LRRK2 regulates a subset of Rab GTPases. Elife 5:e12813. doi: 10.7554/eLife.12813
Taché, Y., Martinez, V., Million, M., and Rivier, J. (1999). Corticotropin-releasing factor and the brain-gut motor response to stress. Can. J. Gastroenterol. 13(Suppl. A), 18A–25A. doi: 10.1155/1999/375916
Taymans, J. M., Van den Haute, C., and Baekelandt, V. (2006). Distribution of PINK1 and LRRK2 in rat and mouse brain. J. Neurochem. 98, 951–961. doi: 10.1111/j.1471-4159.2006.03919.x
Wamsteeker, J. I. C., Füzesi, T., Watts, A. G., and Bains, J. S. (2013). Characterization of corticotropin-releasing hormone neurons in the paraventricular nucleus of the hypothalamus of Crh-IRES-Cre mutant mice. PLoS One. 8:e64943. doi: 10.1371/journal.pone.0064943
Yoon, J. H., Mo, J. S., Kim, M. Y., Ann, E. J., Ahn, J. S., Jo, E. H., et al. (2017). LRRK2 functions as a scaffolding kinase of ASK1-mediated neuronal cell death. Biochim. Biophys. Acta. 1864, 2356–2368. doi: 10.1016/j.bbamcr.2017.09.001
Yun, H. J., Kim, H., Ga, I., Oh, H., Ho, D. H., Kim, J., et al. (2015). An early endosome regulator, Rab5b, is an LRRK2 kinase substrate. J. Biochem. 157, 485–495. doi: 10.1093/jb/mvv005
Zanassi, P., Paolillo, M., Feliciello, A., Avvedimento, E. V., Gallo, V., and Schinelli, S. (2001). cAMP-dependent protein kinase induces cAMP-response element-binding protein phosphorylation via an intracellular calcium release/ERK-dependent pathway in striatal neurons. J. Biol. Chem. 276, 11487–11495. doi: 10.1074/jbc.M007631200
Zhang, F. R., Huang, W., Chen, S. M., Sun, L. D., Liu, H., Li, Y., et al. (2009). Genomewide association study of leprosy. N. Engl. J. Med. 361, 2609–2618. doi: 10.1056/NEJMoa0903753
Keywords: LRRK2, gut dysmotility, stress, paraventricular, signaling/signaling pathway
Citation: Maekawa T, Tsushima H, Kawakami F, Kawashima R, Kodo M, Imai M and Ichikawa T (2019) Leucine-Rich Repeat Kinase 2 Is Associated With Activation of the Paraventricular Nucleus of the Hypothalamus and Stress-Related Gastrointestinal Dysmotility. Front. Neurosci. 13:905. doi: 10.3389/fnins.2019.00905
Received: 15 February 2019; Accepted: 13 August 2019;
Published: 29 August 2019.
Edited by:
Pascal Derkinderen, Université de Nantes, FranceReviewed by:
Jean-Marc Taymans, Institut National de la Santé et de la Recherche Médicale (INSERM), FranceIsabella Russo, University of Brescia, Italy
Copyright © 2019 Maekawa, Tsushima, Kawakami, Kawashima, Kodo, Imai and Ichikawa. 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: Tatsunori Maekawa, maekawa@kitasato-u.ac.jp