Lateral–Medial Dissociation in Orbitofrontal Cortex–Hypothalamus Connectivity
Satoshi Hirose1,2† 
Takahiro Osada1,2† 
Akitoshi Ogawa1 
Masaki Tanaka1 
Hiroyuki Wada3 Yasunori Yoshizawa3 Yoshio Imai3 Toru Machida3,4 Masaaki Akahane3 
Ichiro Shirouzu3 
Seiki Konishi1,2*- 1Department of Neurophysiology, Juntendo University School of Medicine, Tokyo, Japan
- 2Department of Physiology, The University of Tokyo School of Medicine, Tokyo, Japan
- 3Department of Radiology, NTT Medical Center Tokyo, Tokyo, Japan
- 4International University of Health and Welfare, Tokyo, Japan
The orbitofrontal cortex (OFC) is involved in cognitive functions, and is also closely related to autonomic functions. The OFC is densely connected with the hypothalamus, a heterogeneous structure controlling autonomic functions that can be divided into two major parts: the lateral and the medial. Resting-state functional connectivity has allowed us to parcellate the cerebral cortex into putative functional areas based on the changes in the spatial pattern of connectivity in the cerebral cortex when a seed point is moved from one voxel to another. In the present high spatial-resolution fMRI study, we investigate the connectivity-based organization of the OFC with reference to the hypothalamus. The OFC was parcellated using resting-state functional connectivity in an individual subject approach, and then the functional connectivity was examined between the parcellated areas in the OFC and the lateral/medial hypothalamus. We found a functional double dissociation in the OFC: the lateral OFC (the lateral orbital gyrus) was more likely connected with the lateral hypothalamus, whereas the medial OFC (the medial orbital and rectal gyri) was more likely connected with the medial hypothalamus. These results demonstrate the fundamental heterogeneity of the OFC, and suggest a potential neural basis of the OFC–hypothalamic functional interaction.
Introduction
Although the orbitofrontal cortex (OFC) has been well-documented in relation to cognitive functions (Mishkin, 1964; Robbins, 1996; Kringelbach and Rolls, 2004; O’Doherty, 2004; Schoenbaum et al., 2009; Rushworth et al., 2011), its involvement in autonomic functions has been less highlighted. Electric stimulation to the OFC leads to autonomic responses such as changes in blood pressure, heart rate, and respiration rate (Kaada et al., 1949). The OFC is connected, both directly and indirectly, with the hypothalamus (Ongur and Price, 2000; Barbas, 2007). The hypothalamus is recognized as the highest-level center of autonomic functions, and synthesizes the sympathetic/parasympathetic nervous system and the endocrine system (Kandel et al., 2013). The hypothalamus is a small structure (approximately 1 cm3 per hemisphere in humans) that contains several nuclei, and can conventionally be divided into two parts: the lateral and the medial. The lateral part roughly corresponds to one relatively large area, the lateral hypothalamic area, and is known as the hunger center (Anand and Brobeck, 1951; Delgado and Anand, 1953). The medial part of the hypothalamus contains several nuclei, including the ventromedial nucleus known as the satiety center (Hetherington and Ranson, 1940; Miller, 1960), and is related to various autonomic functions such as cardiovascular regulation, hormonal release, and circadian rhythms (Kandel et al., 2013). Thus the OFC may be involved in autonomic functions via its connections with the hypothalamus. Although the OFC–hypothalamus connections have been well described in animal studies (Ongur et al., 1998; Rempel-Clower and Barbas, 1998), the precise anatomical architecture in humans is not well understood.
Functional MRI has allowed us to reveal connectivity-based areal dissociations among small brain structures such as the locus coeruleus in comparison with the ventral tegmental area/substantia nigra pars compacta (Zhang et al., 2016) and the basal nucleus of Meynert in comparison with the ventral striatum (Li et al., 2014). Moreover, areal parcellation using clustering methods has been developed using diffusion-weighted imaging (Behrens et al., 2003; Johansen-Berg et al., 2004; Anwander et al., 2007; Klein et al., 2007; Mars et al., 2011) and resting-state functional connectivity (Yeo et al., 2011; Kahnt et al., 2012; Zhang et al., 2012b; Hutchison et al., 2013, 2015; Onoda and Yamaguchi, 2013; Shen et al., 2013; Long et al., 2014; Eickhoff et al., 2015; Finn et al., 2015; Poldrack et al., 2015; Wang et al., 2015; Reid et al., 2016). In particular, Kahnt et al. (2012) revealed six functionally distinct regions in the OFC using clustering method in resting-state functional connectivity. More recently, resting-state functional connectivity has allowed us to map boundaries of putative functional areas (Margulies et al., 2007; Cohen et al., 2008; Biswal et al., 2010; Nelson et al., 2010; Hirose et al., 2012, 2013; Milham, 2012; Wig et al., 2014a,b; Gordon et al., 2016). For the hypothalamus, despite pioneering attempts to parcellate the human hypothalamus using diffusion-weighted imaging (Schindler et al., 2013; Schönknecht et al., 2013), it is difficult to parcellate the human hypothalamus into individual sub-nuclei with the spatial resolution provided by the current MRI technique.
As one of early steps toward understanding the human OFC–hypothalamic interaction, in the present study, we investigated the connectivity-based architecture of the OFC by applying resting-state functional MRI with high spatial resolution of 2 mm cubic voxels (Hirose et al., 2012, 2013) in an individual-subject approach (Laumann et al., 2015; Wang et al., 2015). The OFC was parcellated into modular areas using the boundary mapping method, and the resting-state functional connectivity, which is known to primarily reflect anatomic connections (Buckner et al., 2008), was then examined between the parcellated OFC areas and the lateral/medial hypothalamus, a conventionally common separation of the hypothalamus (Figure 1).
FIGURE 1. An overview of image data analysis in the present study. Functional MRI images were processed for (1) areal parcellation of the OFC, (2) division of the hypothalamus into the lateral and medial, and (3) calculation of functional connectivity between each of the parcellated areas in the OFC and the lateral/medial hypothalamus.
Materials and Methods
Subjects and MRI Procedures
Twelve healthy right-handed subjects (aged 20–39 years; seven males and five females) participated in the experiments after they gave written informed consent in accordance with the Declaration of Helsinki. The procedures of fMRI experiments were approved by the Institutional Review Board of The University of Tokyo School of Medicine. The experiments were conducted using a 3T scanner (Philips Achieva X 3T Rel. 2.6, Best, The Netherlands). T1-weighted structural images were collected (resolution = 0.81 mm × 0.81 mm × 1.20 mm). Functional imaging was conducted by using gradient-echo echo-planar sequences (TR = 9.0 s, TE = 35 ms, flip angle = 90 deg, FOV = 192 mm × 192 mm, matrix size = 96 × 96, resolution = 2.0 mm × 2.0 mm × 2.0 mm, 75 contiguous transverse slices, ascending interleaved order). The data were sampled using the cubic voxels of 2 mm to minimize signal contamination from the other bank of a sulcus. The higher spatial resolution required a TR of 9 s, but the long TR will not influence the outcome because signals in a lower frequency range of the temporal filter (0.009–0.08 Hz; Fox et al., 2005; Fair et al., 2007), which are known to be predominant in functional connectivity (Salvador et al., 2005), were spared for the resting-state fMRI analysis. During the functional imaging, the subjects were instructed to passively view a fixation point on the screen. One run took about 5 min (35 volumes after discarding the first three volumes), and each subject underwent 90 runs. Typically, one session consisted of 15 runs, and two sessions were administered in 1 day, and subjects underwent these sessions in three separate days.
Figure 2A shows group-averaged functional images of the OFC used in the present study. In the MNI space, the OFC ranged approximately from Z = -30 to -10 and Y = 20–60 (Chiavaras et al., 2001). Despite moderate signal loss due to the sinuses, the orbital surface of the functional images taken in this MRI system appeared relatively well-preserved. As a reference, Figure 2B shows group-averaged functional images of twenty subjects randomly sampled from the data platform of Human Connectome Project1. Although the data platform provides a standard large sample of high-resolution functional images, the large signal dropout centered around the medial OFC in either one hemisphere is less suitable for the boundary mapping analysis employed in the present study.
FIGURE 2. (A) Horizontal sections of group-averaged functional and structural images of the OFC used in the present study. (B) Horizontal sections of group-averaged functional and structural images of the OFC taken from the data platform in the Human Connectome Project. (C) The left and right hypothalami delineated using anatomical landmarks. The hypothalamus was divided into the lateral and medial parts such that their volumes of were matched with one another. (D) Voxel-wise statistical maps of correlation with the whole left and right hypothalami overlaid on top of the anatomic images. The color scale indicates the correlation value. (E) Voxel-wise statistical maps of correlation with the lateral or the medial hypothalamus. The format is similar to that in (D).
Analysis Overview
An overview of the image data analysis of the present study is shown in Figure 1. Acquired functional images were used for (1) areal parcellation in the OFC, (2) delineation of the hypothalamus and its division into the lateral and medial parts, and (3) calculation of functional connectivity between the OFC and the hypothalamus. More specifically, the OFC was parcellated into functional areas, and their centers were defined as regions of interest (ROIs) in the OFC. The hypothalamus was defined based on anatomical landmarks, and was divided into the lateral and medial parts such that their volumes were matched with one another. The functional connectivity was then calculated between each of the centers of the parcellated areas in the OFC and the lateral/medial hypothalamus.
Preprocessing for Resting-State Functional Connectivity
Functional images were preprocessed for the resting-state functional connectivity analysis, and were used later for areal parcellation in the OFC and for calculation of functional connectivity between the OFC and the hypothalamus. Images were realigned and were slice-timing corrected using SPM8 (Worsley and Friston, 1995). Neither spatial normalization nor spatial smoothing was used in order to maintain the higher spatial resolution of the functional images. The following preprocessing procedures are essentially the same as those used in previous literatures of resting-state functional MRI (Fox et al., 2005; Fair et al., 2007). Temporal filters (0.009 Hz < f < 0.08 Hz) were applied to the functional images using FSL (Smith et al., 2004). A general linear model (Worsley and Friston, 1995; Miezin et al., 2000) was used to regress out nuisance signals that correlated with head motion, whole-brain global signal, averaged ventricular signal, and averaged white matter signal (Fox et al., 2005; Fair et al., 2007). Both the global and white matter signals were removed together as covariates of no-interest, and the correlation between the two signals will not affect the results. To examine the effect of global signal regression on the results, the left and right OFC of the representative subject was analyzed without the global signal regression (Supplementary Figures S1 and S2). The results confirmed that the global signal regression had little effect on the analyses of the present study. To avoid phase shifting, we employed frame-wise displacement (FD; Power et al., 2012, 2014) and evaluated the amount of head motion. The FD is instantaneous head motion that can be calculated as locational difference between two successive images in a run. We excluded those runs from analysis where large FD (>0.6 mm) occurred by more than 10% of images in the runs. There were 79 ± 16 runs (mean ± SD) per subject included in the analysis, and the FD was 0.18 ± 0.05 mm (mean ± SD). These runs were concatenated in subsequent correlation analyses.
Areal Parcellation
The probabilistic boundary maps were generated based on the boundary mapping method (Cohen et al., 2008; Nelson et al., 2010; Hirose et al., 2012, 2013; Gordon et al., 2016). The OFC was flattened into the 2D space based on a cortical surface-based analysis (Fischl et al., 1999) using Caret (Van Essen et al., 2001). Each pixel in the 2D flattened cortical space was used as the seed for calculation of correlation with the target voxels. The target voxels were restricted to only those located in the corresponding region in the contralateral hemisphere (Hirose et al., 2012, 2013), without calculation of correlations with all the voxels in the whole brain in the original parcellation method. This modification of the original method has been validated by the observation that a region has the strongest functional connectivity with the corresponding region in the contralateral hemisphere (Stark et al., 2008). It has also been shown that the modified method is efficient in a signal to noise ratio, with a minimal calculation time (Hirose et al., 2012). More specifically, for each seed voxel in the OFC, a spherical region (radius: 2 mm) was generated in the contralateral hemisphere, and the spherical regions were ‘or’-combined across the seed voxels, to form a collection of target voxels. A voxel-wise correlation map in the target region was generated for each seed voxel, and the correlation coefficient was then converted to the Fisher’s z (Fox et al., 2005; Fair et al., 2007).
The analysis procedures after correlation map generation are basically the same as those described in the previous studies (Cohen et al., 2008; Nelson et al., 2010). It is assumed that the pixels at which the spatial pattern of the correlation maps changes drastically represent the boundaries between functional areas. The changes of the spatial pattern of the correlation maps were quantified using the similarity (index of eta2) of the correlation maps between the seeds. The Canny edge detection method (Canny, 1986) was used for the eta2 maps to create a gradient map and also to detect edges. Averaging across the entire sets of binary edge maps generated a probabilistic boundary map where intensity of the pixel represents the probability of the pixel being an edge. The local minima in the gradient map were also detected, and the binary local minimum maps were generated. The binary maps were then averaged to generate a probabilistic ‘center’ map.
We next estimated the overall distance between the centers of the adjacent parcellated areas. Since the probabilistic center maps exhibited a spatially periodic pattern of probability, an autocorrelation analysis was applied to the probabilistic center maps (Hirose et al., 2013). The 2D autocorrelation was calculated in the probabilistic center map using the Signal Processing Toolbox of MATLAB, after spatial smoothing (FWHM: 4 mm). A sliding window of a 50 mm × 50 mm space was moved across the probabilistic center map of OFC (70 mm × 70 mm). The autocorrelation was computed repeatedly for each window position and was averaged across the positions. The autocorrelation value was further clumped along the same radial coordinates to neglect directions.
Hypothalamus Division
The masks that delineated the hypothalamus of individual subjects were created in MNI space after spatial normalization of functional images, based on anatomical landmarks in structural and functional images described in previous MRI studies of the hypothalamus (Schindler et al., 2013; Schönknecht et al., 2013; see Table 1 in Schönknecht et al., 2013 for the detail of definition). The hypothalamus masks were then divided into the lateral and medial parts by a parasagittal plane in the MNI space (Figure 2C). Ideally, the hypothalamus should be divided based on anatomical landmarks, but there are very few anatomical landmarks that can be seen inside the hypothalamus in structural and functional images. So the hypothalamus was split into two parts by a parasagittal plane such that the volumes of the lateral and medial parts of the hypothalamus were matched. Matching volumes helps to minimize the effect of their volumes on the significance level of functional connectivity with the OFC.
The hypothalamus mask in the MNI space created for each subject was converted back to the original individual brain to calculate the functional connectivity in the original subject space. The volume of the hypothalamus in the original subject space was 754 ± 121 mm3 and 779 ± 124 mm3 (mean ± SD) in the left and right hypothalamus, respectively, and was consistent with those reported in the previous studies (Schindler et al., 2013; Schönknecht et al., 2013; Gabery et al., 2015).
Orbitofrontal–Hypothalamic Interaction
The functional connectivity was calculated between each of the parcellated areas in the OFC and the lateral and medial parts of the hypothalamus. The spherical ROIs (radius: 3 mm) were placed on the centers of the parcellated areas in the OFC. The significance level for functional connectivity was set at P < 0.05 for each parcellated area, Bonferroni corrected by the number of centers in the OFC for each subject/hemisphere (35.8 ± 7.0 and 32.9 ± 6.8, mean ± SD, in the left and right OFC).
To investigate the connectivity pattern of many parcellated areas in the OFC with the hypothalamus, the orbital surface of the prefrontal cortex in each subject was sectioned into five parts (the lateral orbital gyrus, the anterior orbital gyrus, the posterior orbital gyrus, the medial orbital gyrus, and the rectal gyrus), based on the sulcus landmarks (the lateral orbital sulcus, the medial orbital sulcus, the transverse orbital sulcus and the olfactory sulcus; Chiavaras and Petrides, 2000; Rolls et al., 2015). The positions of the sulci were provided by the sulcal depth maps in Caret software.
To investigate the connectivity pattern in the OFC in another way, the orbital surface was divided into the lateral and medial parts in a simple geometrical manner. The boundaries surrounding the orbital surface were first delineated, and the three vertices in the orbital surface were determined as the anterior vertex, the lateral posterior vertex, and the medial posterior vertex. The orbital surface was divided into the lateral and the medial parts by the line connecting the anterior vertex and the middle point of the line connecting the lateral posterior and the medial posterior vertices.
As a standard analysis of functional connectivity, a voxel-wise second-level analysis of Fisher’s z maps was separately conducted (Figures 2D,E). After preprocessing including spatial normalization to the standard template and spatial smoothing (FWHM: 6 mm), functional connectivity was calculated between the lateral/medial part of the hypothalamus and each of the voxels in the whole brain, and was converted to Fisher’s z. The resultant correlation z images for individual subjects were entered into a second-level group analysis using a random effect model.
Results
As a standard analysis, the functional connectivity was calculated between the bilateral entire hypothalamus and each voxel in the whole brain. The group analysis showed correlation in the OFC, and major correlation clusters were located in the medial part of the OFC (Figure 2D). The hypothalamus was divided into the lateral and medial parts, and these parts of the hypothalamus were used as separate seed regions. The OFC showed correlation with both the lateral and medial parts of the hypothalamus primarily in the medial part of the OFC, but the difference in the correlation between the lateral and the medial hypothalamus was detected only in (6, 42, -28; medial > lateral, t = 5.4; Figure 2E).
Areal boundaries and centers were then calculated, in an individual subject approach. Figure 3A shows a modular pattern of parcellated areas in the OFC in one representative subject. To confirm the reliability of the areal parcellation, we divided the whole image data in one representative subject into two halves, and each half was applied to the areal parcellation method. The observed pattern was successfully replicated across the two halves of data (Figure 3A). We next measured the overall average distance between adjacent parcellated areas based on the autocorrelation analysis. The probabilistic center map (Figure 3B) was used to calculate the autocorrelation (Figure 3C), and the autocorrelogram was averaged along circles of different radii (Figure 3D). The local maxima was observed at the radius of 11 and 13 mm in the left and right hemispheres respectively (average: 12 mm), indicating that the overall distance between the centers of the parcellated areas was similar to that observed in the posterior inferior frontal cortex (12 mm) reported previously (Hirose et al., 2013).
FIGURE 3. (A) Probabilistic boundary and center maps of the OFC generated by the areal parcellation method. Probabilistic boundaries are shown in blue, with a brightness scale that reflects probability of positive binary edges. White dots indicate the centers of the parcellated areas. Purple curve indicates anatomical landmarks surrounding the orbital surface of the prefrontal cortex. The boundary pattern was preserved even when the whole images were divided into two halves. (B) A probabilistic center map of the OFC. The color scale reflects the probability of local minima in the gradient maps. (C) Autocorrelogram of the probabilistic center maps of the OFC. To measure the overall distance between adjacent centers of the parcellated areas, periodicity of the center map was quantified using an autocorrelation analysis. The color scale indicates the correlation values at particular spatial lags. (D) The correlation values averaged along circles with different radii in the autocorrelogram. Red curves indicate standard error of means of the 12 subjects.
To analyze the functional connectivity pattern of parcellated areas in the OFC with the hypothalamus, the orbital surface was sectioned into five parts in the 2D space based on sulcus landmarks (Chiavaras and Petrides, 2000; Rolls et al., 2015). Figure 4A demonstrates the functional connectivity in the OFC centers in one representative subject that were correlated significantly more with the lateral than the medial hypothalamus, or were correlated significantly more with the medial than the lateral hypothalamus. Table 1 shows the summary of the number of the areas in the five OFC parts that were correlated differentially more with the lateral or the medial hypothalamus in all 12 subjects. The group data of the left and right OFC (Figure 4B) showed that the areas in the lateral orbital gyrus were more likely connected with the lateral hypothalamus, whereas the areas in the medial orbital gyrus and the rectal gyrus were more likely connected with the medial hypothalamus [the lateral orbital gyrus: lateral > medial, t(11) = 3.7, P < 0.01; the medial orbital gyrus: medial > lateral, t(11) = 3.2, P < 0.01; the rectal gyrus: medial > lateral, t(11) = 2.7, P < 0.05]. The anterior and posterior orbital gyri, located in-between the lateral and medial orbital gyri, did not show the differential pattern.
FIGURE 4. (A) Five parts of the OFC in one representative subject that were sectioned using sulcus landmarks. The sections were delineated in the 2D space using the lateral orbital sulcus, the medial orbital sulcus, the transverse orbital sulcus, and the olfactory sulcus, based on Chiavaras and Petrides (2000) and Rolls et al. (2015). The sulcal depth map of Caret in the representative subject was shown in a gray scale. The areal centers that were correlated significantly more with the lateral than the medial hypothalamus are shown in blue, and those was correlated significantly more with the medial than the lateral hypothalamus are shown in purple. LOG, lateral orbital gyrus; AOG, anterior orbital gyrus; POG, posterior orbital gyrus; MOG, medial orbital gyrus; RG, rectal gyrus; los, lateral orbital sulcus; mos, medial orbital sulcus; tos, transverse orbital sulcus; olf, olfactory sulcus. (B) Percentages of areal centers in the left and right OFC that were connected differentially with the lateral and medial hypothalamus. The percentages were group averaged across the 12 subjects in each of the five sections. The bars in the graphs indicate the standard errors of means. ∗P < 0.05, ∗∗P < 0.01.
TABLE 1. The number of areal centers in the sulcus-based OFC sections that were connected differentially with the lateral or the medial hypothalamus.
The functional connectivity pattern in the OFC was analyzed further by sectioning the orbital surface geometrically into the lateral and medial parts in the 2D space. Figure 5A demonstrates the functional connectivity in the OFC centers in the same representative subject, similarly to Figure 4A. Table 2 shows the summary of the number of the areal centers in the lateral and medial parts of the OFC. The group data of the left and right OFC (Figure 5B) showed that the lateral OFC contained the greater number of areas with greater connectivity with the lateral hypothalamus, whereas the medial OFC contained the greater number of areas with greater connectivity with the medial hypothalamus [the lateral OFC: lateral > medial, t(11) = 2.7, P < 0.05; the medial OFC: medial > lateral, t(11) = 3.1, P < 0.01]. The interaction in a two-way ANOVA with the lateral/medial OFC and the lateral/medial hypothalamus as main effects was also significant [F(1,11) = 34.0, P < 0.001].
FIGURE 5. (A) The lateral and medial parts of the orbital surface that were sectioned geometrically in the same representative subject. The sections were delineated in the 2D space using three vertices shown in orange (anterior, lateral posterior, and medial posterior). The format is similar to that in Figure 4A. (B) Percentages of areal centers in the left and right OFC that were connected differentially with the lateral and medial hypothalamus. The percentages were group averaged across the 12 subjects in the lateral and medial OFC. The format is similar to that in Figure 4B. ∗P < 0.05, ∗∗P < 0.01, ∗∗∗P < 0.001.
TABLE 2. The number of areal centers in the geometry-based OFC sections that were connected differentially with the lateral or the medial hypothalamus.
Discussion
The areal parcellation method applied to high-resolution resting-state functional MRI data of individual human subjects demonstrated a basic anatomical architecture of the OFC–hypothalamus interaction: the lateral OFC (the lateral orbital gyrus) was more likely connected with the lateral hypothalamus, whereas the medial OFC (the medial orbital gyrus and rectal gyrus) was more likely connected with the medial hypothalamus. The results provide an important step toward understanding the functional interaction between the OFC and the hypothalamus.
Since the hypothalamus is a small structure of approximately 1 cm3 (Schindler et al., 2013; Schönknecht et al., 2013; Gabery et al., 2015), it is difficult to parcellate the hypothalamus into its individual sub-nuclei. The present resting-state fMRI study only reports the lateral–medial dissociation in the OFC–hypothalamus functional connectivity based on the lateral–medial division of the hypothalamus. Since the resting-state functional connectivity primarily reflects anatomical connections (Buckner et al., 2008), the present study also suggests underlying anatomic connectivity between the OFC and hypothalamus. The indirect connections via the third region may drive the functional connectivity, but it has been demonstrated that the majority of functional connectivity is derived from the direct connections (Honey et al., 2009; Adachi et al., 2012). Moreover, the lateral–medial dissociation is consistent with previous anatomic tracer studies of macaque monkeys (Ongur et al., 1998; Rempel-Clower and Barbas, 1998). These studies revealed differential connection patterns between the lateral hypothalamus and the OFC, and between the medial hypothalamus and the medial prefrontal cortex. Although the present study revealed a dissociation within the OFC, the orbital to medial prefrontal trend that goes along the lateral to medial hypothalamus of the monkeys is consistent with the present dissociation pattern of the lateral to medial OFC trend.
A standard group analysis revealed functional connectivity with the lateral or the medial hypothalamus primarily in the medial OFC (Figure 2E), consistent with the previous study of OFC parcellation (Kahnt et al., 2012). The standard group analysis did not reveal a double dissociation of the lateral–medial differential connectivity (Figure 2E), presumably due to larger individual variations of anatomical architecture in the association cortex (Fischl et al., 2008). One recent study has successfully reported that a right medial OFC region is connected differentially more with the medial than lateral hypothalamus (Kullmann et al., 2014), consistent with the results of the present study (Figure 2E). The effectiveness of the boundary mapping method in defining ROIs with homogeneous connectivity patterns (Gordon et al., 2016) may have helped reveal the lateral–medial double dissociation of the OFC–hypothalamus connectivity in the present study. It is important to discriminate signals of one bank of a sulcus from those of the other bank, and cubic voxels of 2 mm seem sufficient for this purpose. The size of the parcellated areas in the OFC shown in the present study was approximately 12 mm, and the 2 mm voxels also seem sufficient for the parcellation pattern. To examine smaller structures, however, smaller voxels will be useful, but with the need to collect a greater amount of data to compensate for lower signal to noise ratio. Future studies would be required with higher spatial resolution to parcellate the hypothalamus further, as has been done in previous studies of larger subcortical structures such as the thalamus (Behrens et al., 2003), the amygdala (Bzdok et al., 2012), and the striatum (Choi et al., 2012; Janssen et al., 2015).
The lateral–medial dissociation of the OFC–hypothalamus connectivity may explain some of autonomic aspects of the OFC in terms of hypothalamic functions. For example, the somatic marker (physiological arousal) that supports decision making, as measured with skin conductance response (Bechara et al., 2000; Zhang et al., 2012a, 2014, 2015), might be provided from the medial hypothalamus, such as the dorsomedial and the anterior nuclei. These nuclei regulate the sympathetic nervous system, and may explain the loss of skin conductance response following lesions to the OFC (Bechara et al., 2000). Impaired extinction following selective lesions to the medial OFC (Butter, 1969) might be caused by excessive desire for food reward, after loss of interaction with the ventromedial nucleus of the hypothalamus, the satiety center (Hetherington and Ranson, 1940; Miller, 1960). Further understanding of precise anatomical architecture of the OFC–hypothalamus connectivity may reveal some aspects in cognitive tasks that have largely been ignored in neuroscience literatures.
Author Contributions
SH, TO, and SK designed the study, HW, YY, YI, TM, MA, IS, and SK collected data, SH, TO, AO, MT, and SK analyzed the data, SH, TO, and SK wrote the manuscript.
Funding
This work was supported by a Grant-in-Aid for Scientific Research B (22300134) and a grant from Takeda Science Foundation to SK.
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.
Acknowledgment
We thank Ms. Suzuki for technical assistance of the MRI acquisition.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnhum.2016.00244
FIGURE S1 | (A) Probabilistic boundary and center maps with global signal regression (GSR). The right OFC of case 1 is shown. The format is similar to that of Figure 3A. (B) Similar to (A) but without GSR. (C) The correlation values averaged along circles with different radii in the autocorrelogram with GSR. The format is similar to that of Figure 3D. (D) Similar to (C) but without GSR. (E) Five parts of the OFC that were sectioned using sulcus landmarks. Correlation with the lateral/medial hypothalamus was calculated with GSR. The format is similar to that of Figure 4A. (F) Similar to (E) but without GSR. (G) The lateral and medial parts of the orbital surface that were sectioned geometrically. The format is similar to that of Figure 5A. (H) Similar to (G) but without GSR.
FIGURE S2 | The same analyses in Supplementary Figure S1 were performed on the left OFC of case 1. The format is similar to that of Supplementary Figure S1.
Footnotes
References
Adachi, Y., Osada, T., Sporns, O., Watanabe, T., Matsui, T., Miyamoto, K., et al. (2012). Functional connectivity between anatomically unconnected areas is shaped by collective network-level effects in the macaque cortex. Cereb. Cortex 22, 1586–1592. doi: 10.1093/cercor/bhr234
Anand, B. K., and Brobeck, J. R. (1951). Hypothalamic control of food intake in rats and cats. Yale J. Biol. Med. 24, 123–140.
Anwander, A., Tittgemeyer, M., von Cramon, D. Y., Friederici, A. D., and Knösche, T. R. (2007). Connectivity-based parcellation of Broca’s area. Cereb. Cortex 17, 816–825. doi: 10.1093/cercor/bhk034
Barbas, H. (2007). Flow of information for emotions through temporal and orbitofrontal pathways. J. Anatomy. 211, 237–249. doi: 10.1111/j.1469-7580.2007.00777.x
Bechara, A., Damasio, H., and Damasio, A. R. (2000). Emotion, decision making and the orbitofrontal cortex. Cereb. Cortex 10, 295–307. doi: 10.1093/cercor/10.3.295
Behrens, T. E., Johansen-Berg, H., Woolrich, M. W., Smith, S. M., Wheeler-Kingshott, C. A., Boulby, P. A., et al. (2003). Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat. Neurosci. 6, 750–757. doi: 10.1038/nn1075
Biswal, B. B., Mennes, M., Zuo, X. N., Gohel, S., Kelly, C., Smith, S. M., et al. (2010). Toward discovery science of human brain function. Proc. Natl. Acad. Sci. U.S.A. 107, 4734–4739. doi: 10.1073/pnas.0911855107
Buckner, R., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function and relevance to disease. Ann. N. Y. Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011
Butter, C. M. (1969). Perseveration in extinction and in discrimination reversal tasks following selective frontal ablations in Macaca mulatta. Physiol. Behav. 4, 163–171. doi: 10.1016/0031-9384(69)90075-4
Bzdok, D., Laird, A. R., Zilles, K., Fox, P. T., and Eickhoff, S. B. (2012). An investigation of the structural, connectional, and functional subspecialization in the human amygdala. Hum. Brain Mapp. 34, 3247–3266. doi: 10.1002/hbm.22138
Canny, J. (1986). A computational approach to edge detection. IEEE Trans. Pattern Anal. Mach. Intell. PAMI. 8, 679–698. doi: 10.1109/TPAMI.1986.4767851
Chiavaras, M. M., LeGoualher, G., Evans, A., and Petrides, M. (2001). Three-dimensional probabilistic atlas of the human orbitofrontal sulci in Standardized Stereotaxic Space. Neuroimage 13, 479–496. doi: 10.1006/nimg.2000.0641
Chiavaras, M. M., and Petrides, M. (2000). Orbitofrontal sulci of the human and macaque monkey brain. J. Comp. Neurol. 422, 35–54. doi: 10.1002/(SICI)1096-9861(20000619)422:1<35::AID-CNE3>3.0.CO;2-E
Choi, E. Y., Yeo, B. T., and Buckner, R. L. (2012). The organization of the human striatum estimated by intrinsic functional connectivity. J. Neurophysiol. 108, 2242–2263. doi: 10.1152/jn.00270.2012
Cohen, A. L., Fair, D. A., Dosenbach, N. U., Miezin, F. M., Dierker, D., Van Essen, D. C., et al. (2008). Defining functional areas in individual human brains using resting functional connectivity MRI. Neuroimage 41, 45–57. doi: 10.1016/j.neuroimage.2008.01.066
Delgado, J. M. R., and Anand, B. K. (1953). Increase of food intake induced by electrical stimulation of the lateral hypothalamus. Am. J. Physiol. 172, 162–168.
Eickhoff, S. B., Thirion, B., Varoquaux, G., and Bzdok, D. (2015). Connectivity-based parcellation: critique and implications. Hum. Brain Mapp. 36, 4771–4792. doi: 10.1002/hbm.22933
Fair, D. A., Dosenbach, N. U. F., Church, J. A., Cohen, A. L., Brahmbhatt, S., Miezin, F. M., et al. (2007). Development of distinct control networks through segregation and integration. Proc. Natl. Acad. Sci. U.S.A. 104, 13507–13512. doi: 10.1073/pnas.0705843104
Finn, E. S., Shen, X., Scheinost, D., Rosenberg, M. D., Huang, J., Chun, M. M., et al. (2015). Functional connectome fingerprinting: identifying individuals using patterns of brain connectivity. Nat. Neurosci. 18, 1664–1674. doi: 10.1038/nn.4135
Fischl, B., Rajendran, N., Busa, E., Augustinack, J., Hinds, O., Yeo, B. T., et al. (2008). Cortical folding patterns and predicting cytoarchitecture. Cereb. Cortex 18, 1973–1980. doi: 10.1093/cercor/bhm225
Fischl, B., Sereno, M. I., and Dale, A. M. (1999). Cortical surface-based analysis: II: inflation, flattening, and a surface-based coordinate system. Neuroimage 9, 195–207. doi: 10.1006/nimg.1998.0396
Fox, M. D., Snyder, A. Z., Vincent, J. L., Corbetta, M., Van Essen, D. C., and Raichle, M. E. (2005). The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. U.S.A. 102, 9673–9678. doi: 10.1073/pnas.0504136102
Gabery, S., Georgiou-Karistianis, N., Lundh, S. H., Cheong, R. Y., Churchyard, A., Chua, P., et al. (2015). Volumetric analysis of the hypothalamus in huntington disease using 3T MRI: the IMAGE-HD study. PLoS ONE 10:e0117593. doi: 10.1371/journal.pone.0117593
Gordon, E. M., Laumann, T. O., Adeyemo, B., Huckins, J. F., Kelley, W. M., and Petersen, S. E. (2016). Generation and evaluation of a cortical area parcellation from resting-state correlations. Cereb. Cortex 26, 288–303. doi: 10.1093/cercor/bhu239
Hetherington, A. W., and Ranson, S. W. (1940). Hypothalamic lesions and adiposity in the rat. Anat. Rec. 78, 149–172. doi: 10.1111/j.1753-4887.1983.tb07169.x
Hirose, S., Watanabe, T., Jimura, K., Katsura, M., Kunimatsu, A., Abe, O., et al. (2012). Local signal time-series during rest used for areal boundary mapping in individual human brains. PLoS ONE 7:e36496. doi: 10.1371/journal.pone.0036496
Hirose, S., Watanabe, T., Wada, H., Imai, Y., Machida, T., Shirouzu, I., et al. (2013). Functional relevance of micromodules in the human association cortex delineated with high-resolution fMRI. Cereb. Cortex 23, 2863–2871. doi: 10.1093/cercor/bhs268
Honey, C. J., Sporns, O., Cammoun, L., Gigandet, X., Thiran, J. P., Meuli, R., et al. (2009). Predicting human resting-state functional connectivity from structural connectivity. Proc. Natl. Acad. Sci. U.S.A. 106, 2035–2040. doi: 10.1073/pnas.0811168106
Hutchison, R. M., Culhamd, J. C., Flanagan, J. R., Everling, S., and Gallivan, J. P. (2015). Functional subdivisions of medial parieto-occipital cortex in humans and nonhuman primates using resting-state fMRI. Neuroimage 116, 10–29. doi: 10.1016/j.neuroimage.2015.04.068
Hutchison, R. M., Womelsdorf, T., Gati, J. S., Everling, S., and Menon, R. S. (2013). Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques. Hum. Brain Mapp. 34, 2154–2177. doi: 10.1002/hbm.22058
Janssen, R. J., Jylänki, P., Kessels, R. P. C., and van Gerven, M. A. J. (2015). Probabilistic model-based functional parcellation reveals a robust, fine-grained subdivision of the striatum. Neuroimage 119, 398–405. doi: 10.1016/j.neuroimage.2015.06.084
Johansen-Berg, H., Behrens, T. E., Robson, M. D., Drobnjak, I., Rushworth, M. F., Brady, J. M., et al. (2004). Changes in connectivity profiles define functionally distinct regions in human medial frontal cortex. Proc. Natl. Acad. Sci. U.S.A. 101, 13335–13340. doi: 10.1073/pnas.0403743101
Kaada, B. R., Pribram, K. H., and Epstein, J. A. (1949). Respiratory and vascular responses in monkeys from temporal pole, insula, orbital surface and cingulate gyrus. J. Neurophysiol. 12, 347–356.
Kahnt, T., Chang, L. J., Park, S. Q., Heinzle, J., and Haynes, J. D. (2012). Connectivity-based parcellation of the human orbitofrontal cortex. J. Neurosci. 32, 6240–6250. doi: 10.1523/JNEUROSCI.0257-12.2012
Kandel, E. R., Schwartz, J. H., Jessell, T. M., Siegelbaum, S. A., and Hudspeth, A. J. (2013). Principles of Neural Science, 5th Edn. New York, NY: McGraw Hill.
Klein, J. C., Behrens, T. E., Robson, M. D., Mackay, C. E., Higham, D. J., and Johansen-Berg, H. (2007). Connectivity-based parcellation of human cortex using diffusion MRI: establishing reproducibility, validity and observer independence in BA 44/45 and SMA/pre-SMA. Neuroimage 34, 204–211. doi: 10.1016/j.neuroimage.2006.08.022
Kringelbach, M. L., and Rolls, E. T. (2004). The functional neuroanatomy of the human orbitofrontal cortex: evidence from neuroimaging and neuropsychology. Prog. Neurobiol. 72, 341–372. doi: 10.1016/j.pneurobio.2004.03.006
Kullmann, S., Heni, M., Linder, K., Zipfel, S., Haring, H.-U., Veit, R., et al. (2014). Resting-state functional connectivity of the human hypothalamus. Hum. Brain Mapp. 35, 6088–6096. doi: 10.1002/hbm.22607
Laumann, T. O., Gordon, E. M., Adeyemo, B., Snyder, A. Z., Joo, S. J., Chen, M. Y., et al. (2015). Functional system and areal organization of a highly sampled individual human brain. Neuron 87, 657–670. doi: 10.1016/j.neuron.2015.06.037
Li, C. S., Ide, J. S., Zhang, S., Hu, S., Chao, H. H., and Zaborszky, L. (2014). Resting state functional connectivity of the basal nucleus of Meynert in humans: in comparison to the ventral striatum and the effects of age. Neuroimage 97, 321–332. doi: 10.1016/j.neuroimage.2014.04.019
Long, X., Goltz, D., Margulies, D. S., Nierhaus, T., and Villringer, A. (2014). Functional connectivity-based parcellation of the human sensorimotor cortex. Eur. J. Neurosci. 39, 1332–1342. doi: 10.1111/ejn.12473
Margulies, D. S., Kelly, A. M., Uddin, L. Q., Biswal, B. B., Castellanos, F. X., and Milham, M. P. (2007). Mapping the functional connectivity of anterior cingulate cortex. Neuroimage 37, 579–588. doi: 10.1016/j.neuroimage.2007.05.019
Mars, R. B., Jbabdi, S., Sallet, J., O’Reilly, J. X., Croxson, P. L., Olivier, E., et al. (2011). Diffusion-weighted imaging tractography-based parcellation of the human parietal cortex and comparison with human and macaque restingstate functional connectivity. J. Neurosci. 31, 4087–4100. doi: 10.1523/JNEUROSCI.5102-10.2011
Miezin, F. M., Maccotta, L., Ollinger, J. M., Petersen, S. E., and Buckner, R. L. (2000). Characterizing the hemodynamic response: effects of presentation rate, sampling procedure, and the possibility of ordering brain activity based on relative timing. Neuroimage 11, 735–759. doi: 10.1006/nimg.2000.0568
Milham, M. P. (2012). Open neuroscience solutions for the connectome-wide association era. Neuron 73, 214–218. doi: 10.1016/j.neuron.2011.11.004
Mishkin, M. (1964). “Perseveration of central sets after frontal lesions in monkeys,” in The Frontal Granular Cortex and Behavior, eds J. M. Warren and K. Akert (New York, NY: McGraw-Hill), 219–241.
Nelson, S. M., Cohen, A. L., Power, J. D., Wig, G. S., Miezin, F. M., Wheeler, M. E., et al. (2010). A parcellation scheme for human left lateral parietal cortex. Neuron 67, 156–170. doi: 10.1016/j.neuron.2010.05.025
O’Doherty, J. P. (2004). Reward representations and reward-related learning in the human brain: insights from neuroimaging. Curr. Opin. Neurobiol. 14, 769–776. doi: 10.1016/j.conb.2004.10.016
Ongur, D., An, X., and Price, J. L. (1998). Prefrontal cortical projections to the hypothalamus in macaque monkeys. J. Comp. Neurol. 401, 480–505. doi: 10.1002/(SICI)1096-9861(19981130)401:4<480::AID-CNE4>3.3.CO;2-6
Ongur, D., and Price, J. L. (2000). The organization of networks within the orbital and medial prefrontal cortex of rats, monkeys and humans. Cereb. Cortex 10, 206–219. doi: 10.1093/cercor/10.3.206
Onoda, K., and Yamaguchi, S. (2013). Small-worldness and modularity of the resting-state functional brain network decrease with aging. Neurosci. Lett. 556, 104–108. doi: 10.1016/j.neulet.2013.10.023
Poldrack, R. A., Laumann, T. O., Koyejo, O., Gregory, B., Hover, A., Chen, M. Y., et al. (2015). Long-term neural and physiological phenotyping of a single human. Nat. Commun. 6:8885. doi: 10.1038/ncomms9885
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Power, J. D., Mitra, A., Laumann, T. O., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2014). Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage 84, 320–341. doi: 10.1016/j.neuroimage.2013.08.048
Reid, A. T., Lewis, J., Bezgin, G., Khundrakpamb, B., Eickhoff, S. B., McIntosh, A. R., et al. (2016). A cross-modal, cross-species comparison of connectivity measures in the primate brain. Neuroimage 125, 311–331. doi: 10.1016/j.neuroimage.2015.10.057
Rempel-Clower, N. L., and Barbas, H. (1998). Topographic organization of connections between the hypothalamus and prefrontal cortex in the rhesus monkey. J. Comp. Neurol. 398, 393–419. doi: 10.1002/(SICI)1096-9861(19980831)398:3<393::AID-CNE7>3.0.CO;2-V
Robbins, T. W. (1996). Dissociating executive functions of the prefrontal cortex. Philos. Trans. R. Soc. Lond. B Biol. Sci. 351, 1463–1470. doi: 10.1098/rstb.1996.0131
Rolls, E. T., Joliot, M., and Tzourio-Mazoyer, N. (2015). Implementation of a new parcellation of the orbitofrontal cortex in the automated anatomical labeling atlas. Neuroimage 122, 1–5. doi: 10.1016/j.neuroimage.2015.07.075
Rushworth, M. F., Noonan, M. P., Boorman, E. D., Walton, M. E., and Behrens, T. E. (2011). Frontal cortex and reward-guided learning and decision-making. Neuron 70, 1054–1069. doi: 10.1016/j.neuron.2011.05.014
Salvador, R., Suckling, J., Schwarzbauer, C., and Bullmore, E. (2005). Undirected graphs of frequency-dependent functional connectivity in whole brain networks. Philos. Trans. R. Soc. B 360, 937–946. doi: 10.1098/rstb.2005.1645
Schindler, S., Schönknecht, P., Schmidt, L., Anwander, A., Strauß, M., Trampel, R., et al. (2013). Development and evaluation of an algorithm for the computer-assisted segmentation of the human hypothalamus on 7-tesla magnetic resonance images. PLoS ONE 8:e66394. doi: 10.1371/journal.pone.0066394
Schoenbaum, G., Roesch, M. R., Stalnaker, T. A., and Takahashi, Y. K. (2009). A new perspective on the role of the orbitofrontal cortex in adaptive behavior. Nat. Rev. Neurosci. 10, 885–892. doi: 10.1038/nrn2753
Schönknecht, P., Anwander, A., Petzold, F., Schindler, S., Knösche, T. R., Möller, H. E., et al. (2013). Diffusion imaging-based subdivision of the human hypothalamus: a magnetic resonance study with clinical implications. Euro. Arch. Psychiatry Clin. Neurosci. 263, 497–508. doi: 10.1007/s00406-012-0389-5
Shen, X., Tokoglu, F., Papademetris, X., and Constable, R. T. (2013). Groupwise whole-brain parcellation from resting-state fMRI data for network node identification. Neuroimage 82, 403–415. doi: 10.1016/j.neuroimage.2013.05.081
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, 208–219. doi: 10.1016/j.neuroimage.2004.07.051
Stark, D. E., Margulies, D. S., Shehzad, Z. E., Reiss, P., Kelly, A. M., Uddin, L. Q., et al. (2008). Regional variation in interhemispheric coordination of intrinsic hemodynamic fluctuations. J. Neurosci. 51, 13754–13764. doi: 10.1523/JNEUROSCI.4544-08.2008
Van Essen, D. C., Dickson, J., Harwell, J., Hanlon, D., Anderson, C. H., and Drury, H. A. (2001). An integrated software suite for surface-based analyses of cerebral cortex. J. Am. Med. Inform. Assoc. 8, 443–459. doi: 10.1136/jamia.2001.0080443
Wang, D., Buckner, R. L., Fox, M. D., Holt, D. J., Holmes, A., Stoecklein, S., et al. (2015). Parcellating cortical functional networks in individuals. Nat. Neurosci. 18, 1853–1860. doi: 10.1038/nn.4164
Wig, G. S., Laumann, T. O., Cohen, A., Power, J. D., Nelson, S. M., Glasser, M. F., et al. (2014a). Parcellating an individual subject’s cortical and subcortical structures using snowball sampling of resting-state correlations. Cereb. Cortex 24, 2036–2054. doi: 10.1093/cercor/bht056
Wig, G. S., Laumann, T. O., and Petersen, S. E. (2014b). An approach for parcellating human cortical areas using resting-state correlations. Neuroimage 93, 276–291. doi: 10.1016/j.neuroimage.2013.07.035
Worsley, K. J., and Friston, K. J. (1995). Analysis of fMRI time-series revisited again. Neuroimage 2, 173–181. doi: 10.1006/nimg.1995.1023
Yeo, B. T., Krienen, F. M., Sepulcre, J., Sabuncu, M. R., Lashkari, D., Hollinshead, M., et al. (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 106, 1125–1165. doi: 10.1152/jn.00338.2011
Zhang, S., Hu, S., Chao, H. H., Ide, J. S., Luo, X., Farr, O. M., et al. (2014). Ventromedial prefrontal cortex and the regulation of physiological arousal. Soc. Cogn. Affect. Neurosci. 9, 900–908. doi: 10.1093/scan/nst064
Zhang, S., Hu, S., Chao, H. H., and Li, C. S. (2016). Resting-state functional connectivity of the locus coeruleus in humans: in comparison with the ventral tegmental area/substantia nigra pars compacta and the effects of age. Cereb. Cortex doi: 10.1093/cercor/bhv172 [Epub ahead of print].
Zhang, S., Hu, S., Chao, H. H., Luo, X., Farr, O. M., and Li, C. S. (2012a). Cerebral correlates of skin conductance responses in a cognitive task. Neuroimage 62, 1489–1498. doi: 10.1016/j.neuroimage.2012.05.036
Zhang, S., Ide, J. S., and Li, C. S. (2012b). Resting-state functional connectivity of the medial superior frontal cortex. Cereb. Cortex 22, 99–111. doi: 10.1093/cercor/bhr088
Keywords: functional MRI, neuroanatomy, brain mapping, frontal lobe, orbitofrontal cortex
Citation: Hirose S, Osada T, Ogawa A, Tanaka M, Wada H, Yoshizawa Y, Imai Y, Machida T, Akahane M, Shirouzu I and Konishi S (2016) Lateral–Medial Dissociation in Orbitofrontal Cortex–Hypothalamus Connectivity. Front. Hum. Neurosci. 10:244. doi: 10.3389/fnhum.2016.00244
Received: 01 April 2016; Accepted: 09 May 2016;
Published: 26 May 2016.
Edited by:
Tetsuo Kida, National Institute for Physiological Sciences, JapanCopyright © 2016 Hirose, Osada, Ogawa, Tanaka, Wada, Yoshizawa, Imai, Machida, Akahane, Shirouzu and Konishi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Seiki Konishi, skonishi@juntendo.ac.jp
†These authors have contributed equally to this work.