Abstract
The frontal cortex-basal ganglia circuit plays an important role in interval timing. We examined neuronal discharges in the dorsomedial and dorsolateral striatum (DMS and DLS) in rats performing a temporal categorization task and compared them with previously recorded neuronal activity in the medial prefrontal cortex (mPFC). All three structures conveyed significant temporal information, but striatal neurons seldom showed the prolonged, full-interval spanning ramping activity frequently observed in the mPFC. Instead, the majority fired briefly during sample intervals. Also, the precision of neural time decoding became progressively worse with increasing time duration in the mPFC, but not in the striatum. With the caveat that mPFC and striatal units were recorded from different animals, our results suggest that the striatum and mPFC convey temporal information via distinct neural processes.
Introduction
The frontal cortex-basal ganglia circuit has been strongly implicated in interval timing—the estimation of time intervals in the range of seconds to minutes (Buhusi and Meck, ; Meck et al., ). Brain imaging studies in humans have found enhanced blood-oxygen-level dependent signals in the frontal cortex and the striatum during various timing tasks (Rao et al., ; Lewis and Miall, ; Nenadic et al., ; Coull et al., ; Hinton and Meck, ; Jahanshahi et al., ; Grahn and Brett, ; Penney and Vaitilingam, ; Grahn and McAuley, ; Koch et al., ; Teki et al., ; Geiser et al., ). Local lesions or inactivation of the prefrontal cortex (PFC; Glickstein et al., ; Dietrich et al., ; Dietrich and Allen, ; Mangels et al., ; Onoe et al., ; Koch et al., ; Jones et al., ; Kim et al., ) or striatum (Artieda et al., ; Meck, , ; Gibbon et al., ; Rammsayer and Classen, ; Harrington et al., ; Harrington and Haaland, ; Chiba et al., ; Gouvêa et al., ) impair interval timing behavior in both humans and animals. Physiological studies have shown temporal changes in neuronal activity in the PFC, premotor cortex and striatum in monkeys and rats performing various interval timing tasks (Brody et al., ; Matell et al., ; Roux et al., ; Reutimann et al., ; Sakurai et al., ; Genovesio et al., , ; Oshio et al., ; Tanaka, ; Chiba et al., , ; Lebedev et al., ; Jin et al., ; Mita et al., ; Kim et al., ; Knudsen et al., ; Gouvêa et al., ; Mello et al., ; Murakami et al., ; Mendoza et al., ). In addition, optogenetic stimulation of D1 receptor-expressing neurons in the PFC alters interval timing behavior in mice (Narayanan et al., ).
Although a large body of evidence indicates the involvement of the frontal cortex and basal ganglia in interval timing, their respective roles in interval timing and underlying neural processes require further study. To this end, we sought to compare timing-related neural activity between connected regions of the frontal cortex and basal ganglia under the same behavioral setting. We showed previously that in rats performing a temporal categorization task, some medial PFC (mPFC) neurons convey temporal information in the form of monotonically changing (ramping) activity. In these rats, the activity of the recorded mPFC neuronal ensemble was tightly correlated with time interval discrimination behavior (Kim et al., ). Here, to better understand how the frontal cortex-basal ganglia circuit contributes to interval timing, we examined neuronal activity in the striatum in the same temporal categorization task we used in our previous study (Kim et al., ) so we could compare the results with our results in the mPFC.
The cortex-basal ganglia circuit is thought to consist of multiple parallel loops (Alexander et al., ; Alexander and Crutcher, ; Graybiel, ; Redgrave et al., ). In rats, the dorsomedial and dorsolateral striatum (DMS and DLS, respectively) appear to be parts of distinct cortico-basal ganglia loops; they have distinct anatomical connection patterns (McGeorge and Faull, ; Voorn et al., ; Balleine et al., ; Redgrave et al., ; Devan et al., ) and their inactivation or lesions yield dissociable changes in the animal’s choice behavior (Yin et al., , ; Balleine et al., ; White, ). The DMS receives direct projections from the mPFC and is part of the associative cortico-basal ganglia loop, while the DLS receives direct projections from the sensorimotor cortex and is part of the sensorimotor cortico-basal ganglia loop (McGeorge and Faull, ; Voorn et al., ; Balleine et al., ; Redgrave et al., ; Devan et al., ). We recorded neuronal discharges from both the DMS and DLS and compared them to neuronal discharges recorded from the mPFC. We found that both striatal and mPFC neurons carry information about elapsed time, but timing-related neural activity dynamics differ between the two regions.
Materials and Methods
Subjects
Three young male Sprague-Dawley rats (approximately 9–11 weeks old; 280–380 g) were used. After 1 week of extensive handling and water deprivation, their body weights fell to 80%–85% of their ad libitum weights. Once behavioral training began, they were restricted to 30 min of water access immediately after finishing their once-daily behavioral session. Experiments were performed in the dark phase of a 12 h light/dark cycle. All animal care and experimental procedures were performed in accordance with protocols approved by the directives of the Animal Care and Use Committee of the Korea Advanced Institute of Science and Technology (Daejeon, South Korea).
Behavioral Task
We used the same temporal categorization task described in our previous studies (Kim et al., , ). Briefly, the rats were required to discriminate six randomly presented sample intervals into short (3,018, 3,310 and 3,629 ms) or long (3,979, 4,363 and 4,784 ms) and navigate to the corresponding target sites (long-choice target, left in one rat and right in two rats) to obtain a 30 μl water reward (Figure 1A). The animals were required to come back from the target sites to the central stem via lateral alleys. The onset of a sample interval was triggered by the animal’s breaking of the central photobeam. The onset of the sample interval was accompanied by a brief tone (200 ms, 3.3 KHz, 90 dB) and its offset was marked by the lowering of the central bridge. Each session consisted of an average of 250.7 trials (SD 14.1). The animals were over-trained before unit recordings (30–49 and 14–24 days of training before and after electrode implantation, respectively). They also experienced 10 practice trials, five with the shortest interval and five with the longest interval, before each recording session.
Figure 1
Unit Recordings
Two sets of six tetrodes were implanted in the DMS (centered 0.5 mm anterior and 1.8 mm lateral to Bregma; 3.4–5.9 mm ventral to the brain surface) and DLS (centered 0.5 mm anterior and 3.8 mm lateral to Bregma; 3.4–5.9 mm ventral to the brain surface) of well-trained rats under deep anesthesia with sodium pentobarbital (50 mg/kg; Figure 2A). After >1 week of recovery from surgery, the tetrodes were gradually lowered to obtain isolated unit signals. Once unit recording began, the tetrodes were advanced by 100 μm after each daily recording session. Unit signals were amplified 10,000×, filtered between 600 and 6,000 Hz, digitized at 32 kHz, and stored on a personal computer using the Cheetah data acquisition system (Neuralynx, Bozeman, MT, USA). The animal’s head position was monitored by tracking at 60 Hz a set of light-emitting diodes mounted on the headstage. Spike clustering was performed offline using MClust (A.D. Redish). Only well-isolated unit clusters were included in the analysis (L ratio <0.2, isolation distance >15). After completing the recordings, small marking lesions were made by passing an electrolytic current (50 mA, 30 s, cathodal) through one channel of each tetrode. Then, the recording locations were verified histologically as previously described (Baeg et al., ).
Figure 2
Analysis
Choice Behavior
Animal choice data were subjected to the following logistic regression analysis:
where Plong is the probability to choose the long target (averaged across all sessions for each animal and each sample duration), T is the sample duration, and a and b are constants.
Unit Classification
Units with mean firing rates <5 Hz and spike widths >0.26 ms were classified as putative medium spiny neurons (MSNs) and the rest were classified as putative interneurons (Figure 2B).
Activity Duration
A spike density function was generated for each neuron for the longest sample interval (4784 ms) by applying a Gaussian kernel to each spike (σ = 100 ms). Activity half-duration was defined as the duration between the maximal and half-maximal activity of the spike density function. When a spike density function yielded two half-durations, the longer one was used to calculate the activity half-duration for that neuron.
Ramping Activity
The longest sample interval was divided into four equal-duration bins for each neuron. If neural firing rates of the four bins varied significantly (one-way ANOVA, p < 0.05) and neural firing rates changed monotonically across the four bins, the neuron was considered as showing full-interval spanning ramping activity.
Multiple Linear Regression Analysis
To determine whether linear or logarithmic functions better explain individual neuronal activity, we divided each sample interval into 10 equal bins. The relationship between the mean firing rate of a neuron within each time bin and the time since interval onset was described by the following equation:
where S is the trial-by-trial firing rate of a single neuron within a specific time bin, T is the raw or log-transformed time since the interval onset, X, Y and D indicate the animal’s mean lateral head position (X-position), mean vertical head position (Y-position), and overall displacement, respectively, in the corresponding analysis time window, PC denotes the animal’s goal choice in the previous trial, ε is the error term, and a0–a5 are regression coefficients.
Neural Decoding of Temporal Information
Neural decoding of temporal information was done as described previously (Kim et al., ). Briefly, we examined how well neuronal ensemble activity (simultaneously recorded neurons or all neurons pooled across sessions) during the last 500 ms of each sample interval classified sample intervals as short or long using a linear discriminant analysis. When analyzing error trials, we excluded the shortest and longest sample intervals from the analysis because of low error rates (9.7% and 7.9%, respectively). We removed a single trial and generated a linear discriminant function based on the neuronal ensemble activity in the remaining correct trials separated according to the correct target (short vs. long). We repeated this procedure for all correct and error trials and calculated the rate at which the model correctly classified sample interval lengths. The numbers of correct and error trials were equalized across sessions for each sample interval when decoding temporal information using all neurons pooled across sessions.
To assess the amount of temporal information conveyed by neuronal ensemble activity after controlling for the influence of movement-related variables, we used the partial residuals (ε1, Larsen and McCleary, ) of the following regression model for the classification of sample interval lengths:
where S, X, Y and D are neuronal firing rate, the animal’s mean X-position, mean Y-position and overall displacement, respectively, during the last 500 ms of each sample interval.
Statistical Analysis
Student’s t-tests were used to test the statistical significance of regression coefficients. The Kruskal-Wallis test was used to compare the activity durations of DMS, DLS and mPFC neurons. One-way ANOVA was used to compare the decoding performances of DMS, DLS and mPFC neurons. Two-tailed tests were employed for all statistical tests. A p-value < 0.05 was used as the criterion for a significant statistical difference. Data are expressed as means and SD.
Results
Behavior
Across all trials and all animals, we observed a 79.0% (SD 1.4) success rate in choosing the correct target. The probability of choosing the long interval target increased as a function of sample interval across all animals. A logistic regression model accounts well for the relationship between the probability of long target choice and interval duration (eq. 1; animal #1, R2 = 0.958; animal #2, R2 = 0.965; animal #3, R2 = 0.959; Figure 1B).
Neuronal Database
We obtained a total of 490 well-isolated units from the DMS and 577 from the DLS in 71 recording sessions (Figure 2A). In the DMS, 355 (72.4%) were classified as putative MSNs and 135 (27.6%) were classified as putative interneurons; in the DLS, 485 (84.1%) were classified as putative MSNs and 92 (15.9%) were classified as putative interneurons (Figure 2B). The mean discharge rates of the putative MSNs were 0.88 Hz (SD 0.90) in the DMS and 0.73 Hz (SD 0.80) in the DLS. The mean discharge rates of the putative interneurons were 10.82 Hz (SD 8.98) in the DMS and 15.14 Hz (SD 14.99) in the DLS. For comparison, we re-analyzed recordings of 993 well-isolated units from the mPFC of rats performing the same task (Kim et al., ). Of these, 791 putative pyramidal cells had a mean discharge rate of 3.30 Hz (SD 0.13) and 202 putative interneurons had a mean discharge rate of 12.12 Hz (SD 0.69). We included in the analysis only putative MSNs and putative pyramidal cells with mean discharge rates >0.1 Hz during the task (DMS, n = 311, 1.00 Hz (SD 0.90); DLS, n = 429, 0.82 Hz (SD 0.81); mPFC, n = 782, 3.34 Hz (SD 3.76)). We failed to find significant variations in timing-related striatal neural activity along the dorsoventral axis except activity duration (see below).
Activity Duration
In a large proportion of putative MSNs recorded from the DMS and DLS, we observed phasic activity patterns at different times during the sample intervals (examples shown in Figure 3). Figure 4A shows activity profiles during the longest sample interval (4,784 ms) for all the striatal MSNs and mPFC pyramidal neurons we analyzed. As shown, the DLS and DMS MSNs tended toward phasic discharges, with the DLS neurons showing more brief discharges. The putative pyramidal cells in the mPFC tended toward broader temporal activity profiles, with some showing prolonged ramping (i.e., gradually increasing or decreasing during the entire sample interval). To quantify the activity duration for each individual neuron during the longest sample interval, we measured the length of time between the maximum and half-maximum level of each neuron’s spike density function and plotted it as the “Activity half-duration” in Figure 4B. The mean activity half-duration was the longest for the mPFC neurons and shortest for the DLS neurons (DMS, 542.5 ms (SD 477.1); DLS, 417.6 ms (SD 348.9); mPFC, 744.9 ms (SD 578.5); Kruskal-Wallis test, F(2,1472) = 178.84, p < 0.001; Tukey-Kramer post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001; DLS vs. mPFC, p < 0.001; effect size, η2 = 0.108; Figure 4B). Also, the proportions of neurons showing full-interval spanning ramping activity (see “Materials and Methods” section) were 25.8, 16.5 and 12.6% in the mPFC, DMS and DLS, respectively, which deviated significantly from an equal distribution (χ2-test, χ2 = 31.5; p < 0.001; Mann-Whitney post hoc test with Bonferroni correction, DMS vs. DLS, p = 0.146; DMS vs. mPFC, p = 0.001; DLS vs. mPFC, p < 0.001).
Figure 3
Figure 4
To test whether the activity duration varies along the dorsoventral axis of the striatum, we divided the DMS and DLS units into tertiles according to their recording depths and compared their activity durations. More ventral units tended to show shorter activity half-durations in the DMS (dorsal, 605.0 ms (SD 495.7); middle, 580.0 ms (SD 536.5); ventral, 445.0 ms (SD 373.7); Kruskal-Wallis test, F(2,297) = 7.94, p = 0.019; Tukey-Kramer post hoc test, dorsal vs. middle, p = 0.594; dorsal vs. ventral, p = 0.015; middle vs. ventral, p = 0.168) as well as DLS (dorsal, 474.0 ms (SD 362.8); middle, 429.8 ms (SD 393.2); ventral, 348.9 ms (SD 269.0); p = 0.001; dorsal vs. middle, p = 0.096; dorsal vs. ventral, p < 0.001; middle vs. ventral, p = 0.239).
Temporal Information
To quantify the temporal information carried by a neuronal population, we asked how well neuronal ensemble activity during the last 500 ms of each sample interval classifies sample intervals into short and long ones using correct trials. Neural classification of sample-interval length based on simultaneously recorded units was well above chance level (50%) in all areas (ensemble size = 6 neurons, mPFC, 59.0% (SD 0.76) correct classification; t-test, p < 0.001; DMS, 61.4% (1.1), p < 0.001; DLS, 56.9% (SD 1.2), p < 0.001; Figure 5A). Performance increased when the length classification was based on all units pooled across sessions (ensemble size = 200 neurons, mPFC, 88.8% (SD 3.47), p < 0.001; DMS, 97.4% (SD 0.92), p < 0.001; DLS, 92.4% (SD 2.42), p < 0.001; Figure 5B). We then performed a neuron-dropping analysis to examine the relationship between ensemble size and temporal information. At equivalent ensemble sizes, DMS ensembles conveyed significantly more temporal information than the DLS and mPFC ensembles (one-way ANOVA, ensemble size = 100 neurons, F(2,297) = 1741.83, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001; ensemble size = 200 neurons, F(2,297) = 512.01, p < 0.001, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001; Figure 5C).
Figure 5
We have shown previously that timing-related mPFC neural activity in the current task cannot be accounted for by behavioral variations during sample presentation (Kim et al., ). We performed the same analysis for the striatal neurons. For this, we divided the behavioral sessions into quintiles according to temporal information carried by each or whole set of three behavioral variables (the animal’s X-position, Y-position and displacement) and repeated the neural decoding analysis for each group. These behavioral variables did convey temporal information in varying degrees (Figures 6A,B). However, the pattern of neural decoding was generally similar across the quintiles (Figures 6A,B) so that there was no significant variation in the percentage of correct neural classification across the quintiles (one-way ANOVA, DMS, X-position, F(4,25) = 1.94, p = 0.135; Y-position, F(4,25) = 1.38, p = 0.268; displacement, F(4,25) = 0.92, p = 0.469; all variables, F(4,25) = 0.47, p = 0.758; DLS, X-position, F(4,25) = 0.23, p = 0.919; Y-position, F(4,25) = 0.48, p = 0.751; displacement, F(4,25) = 0.70, p = 0.596; all variables, F(4,25) = 0.46, p = 0.767; Figures 6C,D). Thus, neural decoding did not vary significantly as a function of behavioral decoding.
Figure 6
To further address this matter, we used the partial residuals (Larsen and McCleary, ) of the regression model containing the movement-related variables (X-position, Y-position and displacement; eq. 3) instead of raw firing rates for the length classification. All three areas conveyed significant amounts of temporal information, although their performances were somewhat lower (asymptotic values, ~80% correct classification; Figures 6E,F) compared to those using raw firing rates (asymptotic values, ~90% correct classification; Figure 5C). At equivalent ensemble sizes, the DMS and DLS ensembles conveyed significantly more temporal information than the mPFC ensemble (one-way ANOVA, ensemble size = 100 neurons, F(2,297) = 286.05, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001; DLS vs. mPFC, p < 0.001). These results indicate relatively minor contributions of behavioral variations to the temporal information conveyed by striatal neural activity.
Error Trial Analysis
We next analyzed neuronal activity in error trials to determine whether striatal neural activity correlates with the animal’s judgment of time. We excluded the shortest and longest samples in this analysis because they had so few error trials (3,018 ms, n = 4.04 (SD 2.13); 4,784 ms, n = 3.31 (SD 1.92) per session). We then generated discriminant functions using the neuronal ensemble activity during the last 500 ms of each sample interval in correct trials and classified sample intervals into short and long ones using neuronal ensemble activity in error trials. The performance of all three ensembles fell significantly below the level expected by chance (50%; t-test, ensemble size = 100 neurons, DMS, t(198) = 24.14, p < 0.001; DLS, t(198) = 9.75, p < 0.001; mPFC, t(198) = 61.07, p < 0.001; ensemble size = 200 neurons, DMS, t(198) = 20.18, p < 0.001; DLS, t(198) = 9.46, p < 0.001; mPFC, t(198) = 62.13, p < 0.001). At equivalent ensemble sizes, the mPFC ensemble performance (% correct classification) was significantly lower than that of the other ensembles, and the performance of the DLS ensemble was significantly higher (one-way ANOVA, ensemble size = 100 neurons, F(2,297) = 166.05, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001, DLS vs. mPFC, p < 0.001; ensemble size = 200 neurons, F(2,297) = 133.95, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001, DLS vs. mPFC, p < 0.001; Figures 7A,B). Note that poor classification of sample interval length means good prediction of the animal’s target choice in error trials. Thus, mPFC neuronal activity was most well correlated with the animal’s judgment of sample duration, and DLS neuronal activity was only weakly correlated with duration judgment.
Figure 7
The animal may have committed an error because of imprecise timing-related neural activity even with properly maintained attention (“misclassification” trials). Alternatively, in some error trials, the animal may have totally lost track of time, such as due to poorly maintained attention, and chosen a target without referring to timing-related neural activity (“uncertainty” trials). Because the animal’s reaction time is expected to be relatively short and long in misclassification and uncertainty trials, respectively, we examined how reaction time (the time between sample interval offset and the animal’s arrival at a target) varied as a function of sample interval duration in error trials. The reaction time was relatively short in error trials with “difficult” samples intervals (i.e., those close to the classification boundary; 3,629 and 3,979 ms) and relatively long in error trials with “easy” sample intervals (the shortest and longest ones; 3,310 and 4,363 ms; one-way-ANOVA, mPFC, F(5,267) = 6.09; p < 0.001; striatum, F(5,414) = 4.68; p < 0.001; Figure 7C). This is likely due to a high proportion of misclassification for difficult sample intervals. To further explore the relationship between error-trial neural activity and the animal’s target choice, we divided all error trials into two groups according to the reaction time (longer or shorter than the mean reaction time in correct trials for each animal group; mPFC, 1,765.1 ms; striatum, 2,412.4 ms) and repeated the same decoding analysis. For the error trials with relative short reaction times, decoding performance was similarly low for the mPFC and DMS ensembles and close to the chance level for the DLS ensemble (one-way ANOVA, ensemble size = 50 neurons, F(2,2997) = 189.30, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p = 0.114, DLS vs. mPFC, p < 0.001). For the error trials with relatively long reaction times, the performance of the mPFC ensemble showed the lowest decoding performance, followed by the DMS and then the DLS ensembles (one-way ANOVA, ensemble size = 50 neurons, F(2,2997) = 462.18, p < 0.001, Tukey’s HSD post hoc test, DMS vs. DLS, p < 0.001; DMS vs. mPFC, p < 0.001, DLS vs. mPFC, p < 0.001; Figure 7D). Thus, the mPFC neuronal ensemble activity was well correlated with the animal’s duration judgment for all types of error trials, the DMS ensemble activity was well correlated with the animal’s duration judgment for only short reaction-time error trials, and the DLS neuronal ensemble activity was poorly correlated with the animal’s duration judgment for all types of error trials.
Temporal Precision as a Function of Time
We next examined the relationship between temporal discrimination precision and elapsed time. To avoid sensory cue-related neural activity artifacts, we excluded the first 500 ms of each sample interval from the analysis. We then divided each sample interval into 10 equal bins (the duration of each bin is 10% of a given sample duration) and calculated Mahalanobis distances between adjacent bins based on the neuronal ensemble activity in each bin. The Mahalanobis distance measures how far apart two groups of vectors are considering their centers and variances (McCune et al., ). To make it possible to compare the three brain regions, we calculated the Mahalanobis distances after matching ensemble size in each region to that of the smallest ensemble (i.e., DMS, n = 311 neurons). We accomplished this by randomly deleting 118 DLS MSNs and 471 mPFC neurons. We repeated this procedure 100 times and calculated the mean Mahalanobis distances for the DLS and mPFC neuronal ensembles. The Mahalanobis distances for the mPFC neuronal ensembles tended to decrease over time (Figure 8); the slopes of the linear regression curves were significantly different from zero for all sample intervals (t-test, t(7)-values > 7.54, p-values < 0.029) except the shortest one (t(7) = 3.89, p = 0.089). In contrast, for the DMS ensembles, the slope for the shortest sample interval (t(7) = 10.77, p = 0.013) was significantly different from zero, but the slopes for the other intervals were not (t(7)-values < 4.41, p-values > 0.074). For the DLS ensembles, none of the slopes were significantly different from zero (t(7)-values < 3.34, p-values > 0.110; Figure 8). These results are unlikely to reflect different time interval discrimination behavior between the two groups of animals. We found no significant difference in behavior between the two groups (striatum vs. mPFC animals, t-test, performance, 79.0 (SD 1.4) vs. 80.0% (SD 4.6) correct choices, t(4) = −1.11, p = 0.330; parameters of logistic regression analysis (eq. 1), slope, 2.83 (SD 0.23) vs. 3.01 (SD 0.10), t(4) = −1.63, p = 0.244; intercept, −10.86 (SD 0.78) vs. −11.47 (SD 0.36), t(4) = 1.60, p = 0.251; R2-values, 0.961 (SD 0.004) vs. 0.963 (SD 0.009), t(4) = −0.38, p = 0.739).
Figure 8
Linear vs. Logarithmic Activity Profile
We showed previously that mPFC neuronal activity profiles are better described by logarithmic than linear functions (Kim et al., ). This may explain why the precision of neural decoding of time progressively decreases over time. For comparison’s sake, we asked how well linear vs. logarithmic time scales explain the activity of individual striatal neurons. As in our previous study (Kim et al., ), we divided each sample interval into 10 equal bins and used a linear regression model (eq. 2) to examine the dependence of trial-by-trial neuronal activity within each bin on linear or log-transformed time. For the mPFC neurons, we found that the R2 values determined using the logarithmic model were significantly larger than those determined with the linear model for all sample interval durations (Table 1). In contrast, for the DLS and DMS neurons, we found no significant difference between the R2 values determined using the logarithmic or linear models for any sample interval duration (Table 1). Because most DMS and DLS MSNs show phasic activity, we repeated this analysis using only the neurons with relatively long activity profiles (half activity duration >400 ms, DMS, n = 116; DLS, n = 88; mPFC, n = 470). This analysis yielded similar results (Table 1). Thus, in contrast with mPFC neurons, there is no significant difference between the ability of logarithmic or linear models to account for the activity of DMS or DLS neurons.
Table 1
| DMS | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4,784 ms |
| R2 of linear function | 0.115 (0.100) | 0.123 (0.118) | 0.150 (0.134) | 0.145 (0.128) | 0.135 (0.123) | 0.124 (0.111) |
| R2 of log function | 0.116 (0.106) | 0.124 (0.126) | 0.152 (0.146) | 0.145 (0.141) | 0.136 (0.134) | 0.125 (0.116) |
| t-value | −0.371 | −0.382 | −0789 | −0256 | −0.443 | −0.301 |
| p-value | 0.711 | 0.703 | 0.431 | 0.798 | 0.658 | 0.763 |
| DLS | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4,979 ms |
| R2 of linear function | 0.131 (0.135) | 0.128 (0.133) | 0.141 (0.128) | 0.155 (0.151) | 0.130 (0.134) | 0.111 (0.107) |
| R2 of log function | 0.130 (0.134) | 0.127 (0.136) | 0.142 (0.137) | 0.156 (0.153) | 0.133 (0.143) | 0.110 (0.112) |
| t-value | 0.405 | 0.245 | −0.596 | −0275 | −1.307 | 0.737 |
| p-value | 0.686 | 0.807 | 0.552 | 0.783 | 0.192 | 0.461 |
| mPFC | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4,979 ms |
| R2 of linear function | 0.092 (0.102) | 0.103 (0.113) | 0.110 (0.112) | 0.111 (0.113) | 0.108 (0.114) | 0.108 (0.113) |
| R2 of log function | 0.096 (0.109) | 0.107 (0.119) | 0.113 (0.118) | 0.114 (0.122) | 0.110 (0.120) | 0.110 (0.119) |
| t-value | −4.890 | −4.248 | −3.461 | −2.997 | −2.590 | −1.988 |
| p-value | <0.001 | <0.001 | <0.001 | 0.003 | 0.010 | 0.047 |
| DMS | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4,784 ms |
| R2 of linear function | 0.120 (0.112) | 0.133 (0.125) | 0.154 (0.139) | 0.143 (0.126) | 0.131 (0.121) | 0.122 (0.109) |
| R2 of log function | 0.124 (0.126) | 0.136 (0.135) | 0.160 (0.159) | 0.147 (0.147) | 0.136 (0.138) | 0.125 (0.122) |
| t-value | −1.052 | −0.830 | −1.379 | −1.084 | −1.286 | −0.865 |
| p-value | 0.295 | 0.408 | 0.171 | 0.280 | 0.201 | 0.389 |
| DLS | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4,979 ms |
| R2 of linear function | 0.138 (0.144) | 0.144 (0.161) | 0.133 (0.131) | 0.162 (0.172) | 0.113 (0.110) | 0.102 (0.108) |
| R2 of log function | 0.143 (0.144) | 0.149 (0.165) | 0.134 (0.135) | 0.162 (0.170) | 0.116 (0.124) | 0.100 (0.111) |
| t-value | −0.972 | −1.525 | −0.451 | 0.091 | −0.909 | 0.536 |
| p-value | 0.334 | 0.131 | 0.653 | 0.928 | 0.366 | 0.593 |
| mPFC | 3,018 ms | 3,310 ms | 3,629 ms | 3,979 ms | 4,363 ms | 4979 ms |
| R2 of linear function | 0.087 (0.096) | 0.098 (0.110) | 0.104 (0.107) | 0.105 (0.108) | 0.102 (0.112) | 0.100 (0.107) |
| R2 of log function | 0.091 (0.102) | 0.102 (0.117) | 0.107 (0.112) | 0.108 (0.117) | 0.104 (0.116) | 0.102 (0.113) |
| t-value | −3.730 | −3.019 | −2.382 | −2.391 | −2.044 | −1.897 |
| p-value | <0.001 | 0.003 | 0.018 | 0.017 | 0.042 | 0.058 |
Comparison of R2 values for linear vs. logarithmic functions.
R2 values of all neurons (top) or only those with activity half-durations >400 ms (bottom, dorsomedial striatum, DMS, n = 116: dorsolateral striatum, DLS, n = 88; medial prefrontal cortex, mPFC = 470) were compared for the regression models containing linear vs. log-transformed time (eq. 2) for each sample interval duration using t-tests. Data are presented as means (SD).
Discussion
The frontal cortex-basal ganglia circuit is strongly implicated in interval timing (Buhusi and Meck, ; Meck et al., ), but the exact contributions of each circuit component to interval-timing behavior and their underlying neural processes remain unclear. We, therefore, compared neuronal activity patterns in the mPFC, DMS and DLS during identical behavioral tasks. We found that we were able to decode elapsed time from the neuronal ensemble activity recorded from all three regions, meaning all three regions carry substantial temporal information. The neuronal discharge characteristics we recorded from the DMS and DLS during sample interval presentation, however, were distinct from those recorded in the mPFC. DMS and DLS neurons seldom showed the prolonged ramping activity spanning the entire interval that was frequently observed in the mPFC. Instead, most DMS and DLS neurons showed periodic phasic discharges within the sample intervals. In addition, whereas mPFC neuronal activity tended to change logarithmically over time making the neural decoding of elapsed time less precise, this was not the case for DMS or DLS neurons.
Ramping and sequential discharges, which are both found ubiquitously in the brain (Durstewitz and Seamans, ; Buzsáki, ), are excellent candidate neural processes for conveying temporal information. Both ramping activity and sequential discharges are found in the mPFC as well as striatum (e.g., Kim et al., ; Emmons et al., ; Tiganj et al., ). However, in our task in which the confounding effect of movement-related neural activity is relatively small (Figure 6), activity durations were significantly shorter for striatal than mPFC neurons. Moreover, ramping activity spanning the entire duration of the longest sample interval (4874 ms) was less frequently found in the striatum. Previous studies also have found phasic discharges of striatal neurons under diverse behavioral conditions (e.g., Kimura et al., ; Apicella et al., ; Lau and Glimcher, ). These results raise the possibility that the mPFC and striatum may carry temporal information via distinct neural processes. Why is prolonged ramping observed rarely in the striatum? This may be due to a fundamental difference in the organization of the neural circuits in each region. Cortical circuits mainly consist of excitatory neurons interconnected by recurrent collaterals, whereas the principal striatal neurons (i.e., the MSNs) are inhibitory. Excitatory network activity may drive ramping activity in some PFC neurons (Mongillo et al., ). Alternatively, mPFC pyramidal neurons and striatal MSNs may have different channel compositions. In other words, mPFC neurons may exhibit ramping activity because of unique internal channel dynamics (Durstewitz, ) that are absent in striatal neurons. Additional studies will be necessary to clarify the mechanism underlying this difference.
There is disagreement about whether the brain represents time on a linear (Gibbon, ; Gibbon and Church, ; Roberts, ; Church and Gibbon, ; Gallistel, ; Wearden and Jones, ) or logarithmic (Church and Deluty, ; Staddon and Higa, ; Roberts, ; Yi, ) scale. Some studies have reported linear changes in neural activity over time (e.g., Komura et al., ; Machens et al., ), but we showed that activity profiles of mPFC neurons are better described by logarithmic than linear functions in the current task (Kim et al., ). The logarithmic encoding of time means that as time elapses, its representation becomes progressively less precise. This was indeed the case for mPFC neurons regardless of whether we analyzed the whole population (Kim et al., ), analyzed the phasic-firing “time cells” alone (Tiganj et al., ), or excluded the putative sensory-related responses (Figure 8). In contrast, when we excluded the putative sensory-related responses (i.e., the first 500 ms) from the analysis, the DMS and DLS showed similar levels of precision in temporal discrimination over time. Also, unlike with the mPFC neurons, logarithmic functions were no better at accounting for the temporal profiles of individual striatal neurons than linear functions. With the caveats that we recorded the striatal and mPFC neurons from different animals and that the results of comparing linear vs. logarithmic functions were equivocal for striatal neurons, our results raise the possibility that different brain areas may represent time on linear and logarithmic scales in parallel.
Our finding that DMS neurons convey more temporal information and predict animal choices better in error trials than DLS neurons is consistent with the anatomical connections between the mPFC and striatum (the mPFC directly projects to the DMS, but not the DLS; McGeorge and Faull, ; Voorn et al., ; Balleine et al., ; Redgrave et al., ; Devan et al., ) and recent physiological studies that suggest the transfer of temporal information from the PFC to the striatum (Emmons et al., ; Wang et al., ). Although it is unknown whether other cortical regions (e.g., the sensorimotor cortex) similarly convey temporal information to the DLS, our results suggest the mPFC may be an important source of temporal information for the DMS. Then a question arises as to why we found evidence for logarithmic encoding of time in the mPFC, but not in the striatum. It would be informative in the future to examine the effects of manipulating the mPFC on DMS timing-related neural activity in order to understand whether and how mPFC neural activity contributes to temporal information processing in the DMS (see Emmons et al., ; Wang et al., ). It is also unclear whether close interactions between the PFC and striatum are necessary for all types of interval timing behavior (Matell and Meck, ; Meck et al., ; Agostino and Cheng, ) or they may work independently to be in charge of controlling different types of interval timing behavior. A recent study has shown that striatal units predict the timing of licking responses better than orbitofrontal cortical units in head-fixed mice performing an appetitive trace conditioning task (Bakhurin et al., ). In our study, the mPFC neural activity was most well correlated with the animal’s target choice in all types of error trials. The mPFC and striatum may play more important roles in controlling where to navigate according to duration judgment (as in our task) and when to emit a specific motor response (such as licking), respectively, for example. In general, it is largely unknown how temporal information contained in temporally-changing neural activity in a given brain area contributes to behavioral control. Clearly, further studies will be necessary to clarify the functional roles of timing-related neural activity in different areas of the brain.
Statements
Author contributions
JK and MWJ designed the study. JK performed the experiments. JK, DK and MWJ analyzed the data. MWJ wrote the article with input from all authors.
Funding
This work was supported by the Institute for Basic Science (IBS-R002-G1; MWJ).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AgostinoP. V.ChengR.-K. (2016). Contributions of dopaminergic signaling to timing accuracy and precision. Curr. Opin. Behav. Sci.8, 153–160. 10.1016/j.cobeha.2016.02.013
2
AlexanderG. E.CrutcherM. D. (1990). Functional architecture of basal ganglia circuits: neural substrates of parallel processing. Trends Neurosci.13, 266–271. 10.1016/0166-2236(90)90107-l
3
AlexanderG. E.DeLongM. R.StrickP. L. (1986). Parallel organization of functionally segregated circuits linking basal ganglia and cortex. Annu. Rev. Neurosci.9, 357–381. 10.1146/annurev.ne.09.030186.002041
4
ApicellaP.LjungbergT.ScarnatiE.SchultzW. (1991). Responses to reward in monkey dorsal and ventral striatum. Exp. Brain Res.85, 491–500. 10.1007/bf00231732
5
ArtiedaJ.PastorM. A.LacruzF.ObesoJ. A. (1992). Temporal discrimination is abnormal in Parkinson’s disease. Brain115, 199–210. 10.1093/brain/115.1.199
6
BaegE. H.KimY. B.JangJ.KimH. T.Mook-JungI.JungM. W. (2001). Fast spiking and regular spiking neural correlates of fear conditioning in the medial prefrontal cortex of the rat. Cereb. Cortex11, 441–451. 10.1093/cercor/11.5.441
7
BakhurinK. I.GoudarV.ShobeJ. L.ClaarL. D.BuonomanoD. V.MasmanidisS. C. (2017). Differential encoding of time by prefrontal and striatal network dynamics. J. Neurosci.37, 854–870. 10.1523/JNEUROSCI.1789-16.2016
8
BalleineB. W.DelgadoM. R.HikosakaO. (2007). The role of the dorsal striatum in reward and decision-making. J. Neurosci.27, 8161–8165. 10.1523/jneurosci.1554-07.2007
9
BrodyC. D.HernandezA.ZainosA.RomoR. (2003). Timing and neural encoding of somatosensory parametric working memory in macaque prefrontal cortex. Cereb. Cortex13, 1196–1207. 10.1093/cercor/bhg100
10
BuhusiC. V.MeckW. H. (2005). What makes us tick? Functional and neural mechanisms of interval timing. Nat. Rev. Neurosci.6, 755–765. 10.1038/nrn1764
11
BuzsákiG. (2010). Neural syntax: cell assemblies, synapsembles, and readers. Neuron68, 362–385. 10.1016/j.neuron.2010.09.023
12
ChibaA.OshioK.InaseM. (2008). Striatal neurons encoded temporal information in duration discrimination task. Exp. Brain Res.186, 671–676. 10.1007/s00221-008-1347-3
13
ChibaA.OshioK.InaseM. (2015). Neuronal representation of duration discrimination in the monkey striatum. Physiol. Rep.3:e12283. 10.14814/phy2.12283
14
ChurchR. M.DelutyM. Z. (1977). Bisection of temporal intervals. J. Exp. Psychol. Anim. Behav. Process.3, 216–228. 10.1037/0097-7403.3.3.216
15
ChurchR. M.GibbonJ. (1982). Temporal generalization. J. Exp. Psychol. Anim. Behav. Process.8, 165–186. 10.1037/0097-7403.8.2.165
16
CoullJ. T.VidalF.NazarianB.MacarF. (2004). Functional anatomy of the attentional modulation of time estimation. Science303, 1506–1508. 10.1126/science.1091573
17
DevanB. D.HongN. S.McDonaldR. J. (2011). Parallel associative processing in the dorsal striatum: segregation of stimulus-response and cognitive control subregions. Neurobiol. Learn. Mem.96, 95–120. 10.1016/j.nlm.2011.06.002
18
DietrichA.AllenJ. D. (1998). Functional dissociation of the prefrontal cortex and the hippocampus in timing behavior. Behav. Neurosci.112, 1043–1047. 10.1037//0735-7044.112.5.1043
19
DietrichA.FrederickD. L.AllenJ. D. (1997). The effects of total and subtotal prefrontal cortex lesions on the timing ability of the rat. Psychobiol.25, 191–201.
20
DurstewitzD. (2003). Self-organizing neural integrator predicts interval times through climbing activity. J. Neurosci.23, 5342–5353. 10.1523/jneurosci.23-12-05342.2003
21
DurstewitzD.SeamansJ. K. (2006). Beyond bistability: biophysics and temporal dynamics of working memory. Neuroscience139, 119–133. 10.1016/j.neuroscience.2005.06.094
22
EmmonsE. B.De CorteB. J.KimY.ParkerK. L.MatellM. S.NarayananN. S. (2017). Rodent medial frontal control of temporal processing in the dorsomedial striatum. J. Neurosci.37, 8718–8733. 10.1523/JNEUROSCI.1376-17.2017
23
GallistelC. R. (1999). Can a decay process explain the timing of conditioned responses?J. Exp. Anal. Behav.71, 264–271. 10.1901/jeab.1999.71-264
24
GeiserE.NotterM.GabrieliJ. D. (2012). A corticostriatal neural system enhances auditory perception through temporal context processing. J. Neurosci.32, 6177–6182. 10.1523/JNEUROSCI.5153-11.2012
25
GenovesioA.TsujimotoS.WiseS. P. (2006). Neuronal activity related to elapsed time in prefrontal cortex. J. Neurophysiol.95, 3281–3285. 10.1152/jn.01011.2005
26
GenovesioA.TsujimotoS.WiseS. P. (2009). Feature- and order-based timing representations in the frontal cortex. Neuron63, 254–266. 10.1016/j.neuron.2009.06.018
27
GibbonJ. (1977). Scalar expectancy theory and Weber’s law in animal timing. Psychol. Rev.84, 279–279. 10.1037//0033-295x.84.3.279
28
GibbonJ.ChurchR. M. (1981). Time left: linear versus logarithmic subjective time. J. Exp. Psychol. Anim. Behav. Process.7, 87–108. 10.1037//0097-7403.7.2.87
29
GibbonJ.MalapaniC.DaleC. L.GallistelC. (1997). Toward a neurobiology of temporal cognition: advances and challenges. Curr. Opin. Neurobiol.7, 170–184. 10.1016/s0959-4388(97)80005-0
30
GlicksteinM.QuigleyW. A.StebbinsW. C. (1964). Effect of frontal and parietal lesions on timing behavior in monkeys. Psychon. Sci.1, 265–266. 10.3758/bf03342902
31
GouvêaT. S.MonteiroT.MotiwalaA.SoaresS.MachensC.PatonJ. J. (2015). Striatal dynamics explain duration judgments. Elife4:e11386. 10.1101/020883
32
GrahnJ. A.BrettM. (2007). Rhythm and beat perception in motor areas of the brain. J. Cogn. Neurosci.19, 893–906. 10.1162/jocn.2007.19.5.893
33
GrahnJ. A.McAuleyJ. D. (2009). Neural bases of individual differences in beat perception. Neuroimage47, 1894–1903. 10.1016/j.neuroimage.2009.04.039
34
GraybielA. M. (2008). Habits, rituals, and the evaluative brain. Annu. Rev. Neurosci.31, 359–387. 10.1146/annurev.neuro.29.051605.112851
35
HarringtonD. L.HaalandK. Y. (1999). Neural underpinnings of temporal processing: a review of focal lesion, pharmacological and functional imaging research. Rev. Neurosci.10, 91–116. 10.1515/revneuro.1999.10.2.91
36
HarringtonD. L.HaalandK. Y.KnightR. T. (1998). Cortical networks underlying mechanisms of time perception. J. Neurosci.18, 1085–1095. 10.1523/jneurosci.18-03-01085.1998
37
HintonS. C.MeckW. H. (2004). Frontal-striatal circuitry activated by human peak-interval timing in the supra-seconds range. Cogn. Brain Res.21, 171–182. 10.1016/j.cogbrainres.2004.08.005
38
JahanshahiM.JonesC. R.DirnbergerG.FrithC. D. (2006). The substantia nigra pars compacta and temporal processing. J. Neurosci.26, 12266–12273. 10.1523/JNEUROSCI.2540-06.2006
39
JinD. Z.FujiiN.GraybielA. M. (2009). Neural representation of time in cortico-basal ganglia circuits. Proc. Natl. Acad. Sci. U S A106, 19156–19161. 10.1073/pnas.0909881106
40
JonesC. R.RosenkranzK.RothwellJ. C.JahanshahiM. (2004). The right dorsolateral prefrontal cortex is essential in time reproduction: an investigation with repetitive transcranial magnetic stimulation. Exp. Brain Res.158, 366–372. 10.1007/s00221-004-1912-3
41
KimJ.GhimJ.-W.LeeJ. H.JungM. W. (2013). Neural correlates of interval timing in rodent prefrontal cortex. J. Neurosci.33, 13834–13847. 10.1523/jneurosci.1443-13.2013
42
KimJ.JungA. H.ByunJ.JoS.JungM. W. (2009). Inactivation of medial prefrontal cortex impairs time interval discrimination in rats. Front. Behav. Neurosci.3:38. 10.3389/neuro.08.038.2009
43
KimuraM.KatoM.ShimazakiH. (1990). Physiological properties of projection neurons in the monkey striatum to the globus pallidus. Exp. Brain Res.82, 672–676. 10.1007/bf00228811
44
KnudsenE. B.PowersM. E.MoxonK. A. (2014). Dissociating movement from movement timing in the rat primary motor cortex. J. Neurosci.34, 15576–15586. 10.1523/jneurosci.1816-14.2014
45
KochG.OliveriM.CaltagironeC. (2009). Neural networks engaged in milliseconds and seconds time processing: evidence from transcranial magnetic stimulation and patients with cortical or subcortical dysfunction. Philos. Trans. R. Soc. Lond. B Biol. Sci.364, 1907–1918. 10.1098/rstb.2009.0018
46
KochG.OliveriM.TorrieroS.CaltagironeC. (2003). Underestimation of time perception after repetitive transcranial magnetic stimulation. Neurology60, 1844–1846. 10.1212/wnl.60.11.1844
47
KomuraY.TamuraR.UwanoT.NishijoH. (2001). Retrospective and prospective coding for predicted reward in the sensory thalamus. Nature412, 546–549. 10.1038/35087595
48
LarsenW. A.McClearyS. J. (1972). The use of partial residual plots in regression analysis. Technometrics14, 781–790. 10.2307/1267305
49
LauB.GlimcherP. W. (2008). Value representations in the primate striatum during matching behavior. Neuron58, 451–463. 10.1016/j.neuron.2008.02.021
50
LebedevM. A.O’DohertyJ. E.NicolelisM. A. (2008). Decoding of temporal intervals from cortical ensemble activity. J. Neurophysiol.99, 166–186. 10.1152/jn.00734.2007
51
LewisP. A.MiallR. C. (2003). Distinct systems for automatic and cognitively controlled time measurement: evidence from neuroimaging. Curr. Opin. Neurobiol.13, 250–255. 10.1016/s0959-4388(03)00036-9
52
MachensC. K.RomoR.BrodyC. D. (2010). Functional, but not anatomical, separation of “what” and “when” in prefrontal cortex. J. Neurosci.30, 350–360. 10.1523/JNEUROSCI.3276-09.2010
53
MangelsJ. A.IvryR. B.ShimizuN. (1998). Dissociable contributions of the prefrontal and neocerebellar cortex to time perception. Cogn. Brain Res.7, 15–39. 10.1016/s0926-6410(98)00005-6
54
MatellM. S.MeckW. H. (2004). Cortico-striatal circuits and interval timing: coincidence detection of oscillatory processes. Cogn. Brain Res.21, 139–170. 10.1016/j.cogbrainres.2004.06.012
55
MatellM. S.MeckW. H.NicolelisM. A. (2003). Interval timing and the encoding of signal duration by ensembles of cortical and striatal neurons. Behav. Neurosci.117, 760–773. 10.1037/0735-7044.117.4.760
56
McCuneB.GraceJ.UrbanD. (2002). Analysis of Ecological Communities, Chapter 6.Gleneden Beach, OR: MJM Software Design.
57
McGeorgeA.FaullR. (1989). The organization of the projection from the cerebral cortex to the striatum in the rat. Neuroscience29, 503–537. 10.1016/0306-4522(89)90128-0
58
MeckW. H. (1996). Neuropharmacology of timing and time perception. Cogn. Brain Res.3, 227–242. 10.1016/0926-6410(96)00009-2
59
MeckW. H. (2006). Neuroanatomical localization of an internal clock: a functional link between mesolimbic, nigrostriatal, and mesocortical dopaminergic systems. Brain Res.1109, 93–107. 10.1016/j.brainres.2006.06.031
60
MeckW. H.PenneyT. B.PouthasV. (2008). Cortico-striatal representation of time in animals and humans. Curr. Opin. Neurobiol.18, 145–152. 10.1016/j.conb.2008.08.002
61
MelloG. B.SoaresS.PatonJ. J. (2015). A scalable population code for time in the striatum. Curr. Biol.25, 1113–1122. 10.1016/j.cub.2015.02.036
62
MendozaG.MendezJ. C.PerezO.PradoL.MerchantH. (2018). Neural basis for categorical boundaries in the primate pre-SMA during relative categorization of time intervals. Nat. Commun.9:1098. 10.1038/s41467-018-03482-8
63
MitaA.MushiakeH.ShimaK.MatsuzakaY.TanjiJ. (2009). Interval time coding by neurons in the presupplementary and supplementary motor areas. Nat. Neurosci.12, 502–507. 10.1038/nn.2272
64
MongilloG.AmitD. J.BrunelN. (2003). Retrospective and prospective persistent activity induced by Hebbian learning in a recurrent cortical network. Eur. J. Neurosci.18, 2011–2024. 10.1046/j.1460-9568.2003.02908.x
65
MurakamiM.ShteingartH.LoewensteinY.MainenZ. F. (2017). Distinct sources of deterministic and stochastic components of action timing decisions in rodent frontal cortex. Neuron94, 908.e7–919.e7. 10.1016/j.neuron.2017.04.040
66
NarayananN. S.LandB. B.SolderJ. E.DeisserothK.DiLeoneR. J. (2012). Prefrontal D1 dopamine signaling is required for temporal control. Proc. Natl. Acad. Sci. U S A109, 20726–20731. 10.1073/pnas.1211258109
67
NenadicI.GaserC.VolzH.-P.RammsayerT.HägerF.SauerH. (2003). Processing of temporal information and the basal ganglia: new evidence from fMRI. Exp. Brain Res.148, 238–246. 10.1007/s00221-002-1188-4
68
OnoeH.KomoriM.OnoeK.TakechiH.TsukadaH.WatanabeY. (2001). Cortical networks recruited for time perception: a monkey positron emission tomography (PET) study. Neuroimage13, 37–45. 10.1006/nimg.2000.0670
69
OshioK.ChibaA.InaseM. (2006). Delay period activity of monkey prefrontal neurones during duration-discrimination task. Eur. J. Neurosci.23, 2779–2790. 10.1111/j.1460-9568.2006.04781.x
70
PaxinosG.WatsonC. (1998). A Stereotaxic Atlas of the Rat Brain.New York, NY: Academic Press.
71
PenneyT. B.VaitilingamL. (2008). “Imaging time,” in Psychology of Time, ed. GrondinS. (Bingley: Emerald Group), 261–294.
72
RammsayerT.ClassenW. (1997). Impaired temporal discrimination in Parkinson’s disease: temporal processing of brief durations as an indicator of degeneration of dopaminergic neurons in the basal ganglia. Int. J. Neurosci.91, 45–55. 10.3109/00207459708986364
73
RaoS. M.HarringtonD. L.HaalandK. Y.BobholzJ. A.CoxR. W.BinderJ. R. (1997). Distributed neural systems underlying the timing of movements. J. Neurosci.17, 5528–5535. 10.1523/jneurosci.17-14-05528.1997
74
RedgraveP.RodriguezM.SmithY.Rodriguez-OrozM. C.LehericyS.BergmanH.et al. (2010). Goal-directed and habitual control in the basal ganglia: implications for Parkinson’s disease. Nat. Rev. Neurosci.11, 760–772. 10.1038/nrn2915
75
ReutimannJ.YakovlevV.FusiS.SennW. (2004). Climbing neuronal activity as an event-based cortical representation of time. J. Neurosci.24, 3295–3303. 10.1523/jneurosci.4098-03.2004
76
RobertsS. (1981). Isolation of an internal clock. J. Exp. Psychol. Anim. Behav. Process.7, 242–268. 10.1037//0097-7403.7.3.242
77
RobertsW. A. (2006). Evidence that pigeons represent both time and number on a logarithmic scale. Behav. Processes72, 207–214. 10.1016/j.beproc.2006.03.002
78
RouxS.CoulmanceM.RiehleA. (2003). Context-related representation of timing processes in monkey motor cortex. Eur. J. Neurosci.18, 1011–1016. 10.1046/j.1460-9568.2003.02792.x
79
SakuraiY.TakahashiS.InoueM. (2004). Stimulus duration in working memory is represented by neuronal activity in the monkey prefrontal cortex. Eur. J. Neurosci.20, 1069–1080. 10.1111/j.1460-9568.2004.03525.x
80
StaddonJ. E.HigaJ. J. (1999). Time and memory: towards a pacemaker-free theory of interval timing. J. Exp. Anal. Behav.71, 215–251. 10.1901/jeab.1999.71-215
81
TanakaM. (2007). Cognitive signals in the primate motor thalamus predict saccade timing. J. Neurosci.27, 12109–12118. 10.1523/jneurosci.1873-07.2007
82
TekiS.GrubeM.KumarS.GriffithsT. D. (2011). Distinct neural substrates of duration-based and beat-based auditory timing. J. Neurosci.31, 3805–3812. 10.1523/jneurosci.5561-10.2011
83
TiganjZ.JungM. W.KimJ.HowardM. W. (2017). Sequential firing codes for time in rodent medial prefrontal cortex. Cereb. Cortex27, 5663–5671. 10.1093/cercor/bhw336
84
VoornP.VanderschurenL. J.GroenewegenH. J.RobbinsT. W.PennartzC. M. (2004). Putting a spin on the dorsal-ventral divide of the striatum. Trends Neurosci.27, 468–474. 10.1016/j.tins.2004.06.006
85
WangJ.NarainD.HosseiniE. A.JazayeriM. (2018). Flexible timing by temporal scaling of cortical responses. Nat. Neurosci.21, 102–110. 10.1038/s41593-017-0028-6
86
WeardenJ.JonesL. A. (2007). Is the growth of subjective time in humans a linear or nonlinear function of real time?Q. J. Exp. Psychol.60, 1289–1302. 10.1080/17470210600971576
87
WhiteN. M. (2009). Some highlights of research on the effects of caudate nucleus lesions over the past 200 years. Behav. Brain Res.199, 3–23. 10.1016/j.bbr.2008.12.003
88
YiL. (2009). Do rats represent time logarithmically or linearly?Behav. Processes81, 274–279. 10.1016/j.beproc.2008.10.004
89
YinH. H.KnowltonB. J.BalleineB. W. (2004). Lesions of dorsolateral striatum preserve outcome expectancy but disrupt habit formation in instrumental learning. Eur. J. Neurosci.19, 181–189. 10.1111/j.1460-9568.2004.03095.x
90
YinH. H.OstlundS. B.KnowltonB. J.BalleineB. W. (2005). The role of the dorsomedial striatum in instrumental conditioning. Eur. J. Neurosci.22, 513–523. 10.1111/j.1460-9568.2005.04218.x
Summary
Keywords
interval timing, temporal categorization task, prefrontal cortex, striatum, rats
Citation
Kim J, Kim D and Jung MW (2018) Distinct Dynamics of Striatal and Prefrontal Neural Activity During Temporal Discrimination. Front. Integr. Neurosci. 12:34. doi: 10.3389/fnint.2018.00034
Received
11 May 2018
Accepted
24 July 2018
Published
13 August 2018
Volume
12 - 2018
Edited by
Mark Laubach, American University, United States
Reviewed by
Matthew S. Matell, Villanova University, United States; Marcelo S. Caetano, Universidade Federal do ABC, Brazil
Updates
Copyright
© 2018 Kim, Kim and Jung.
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: Min Whan Jung mwjung@kaist.ac.kr
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.