Abstract
The frontal cortex-striatum circuit plays a pivotal role in adaptive goal-directed behaviors. However, it remains unclear how decision-related signals are mediated through cross-regional transmission between the medial frontal cortex and the striatum by neuronal ensembles in making decision based on outcomes of past action. Here, we analyzed neuronal ensemble activity obtained through simultaneous multiunit recordings in the secondary motor cortex (M2) and dorsal striatum (DS) in rats performing an outcome-based left-or-right choice task. By adopting tensor component analysis (TCA), a single-trial–based unsupervised dimensionality reduction approach, for concatenated ensembles of M2 and DS neurons, we identified distinct three spatiotemporal neural dynamics (TCA components) at the single-trial level specific to task-relevant variables. Choice-position–selective neural dynamics reflected the positions chosen and was correlated with the trial-to-trial fluctuation of behavioral variables. Intriguingly, choice-pattern–selective neural dynamics distinguished whether the incoming choice was a repetition or a switch from the previous choice before a response choice. Other neural dynamics was selective to outcome and increased within-trial activity following response. Our results demonstrate how the concatenated ensembles of M2 and DS process distinct features of decision-related signals at various points in time. Thereby, the M2 and DS collaboratively monitor action outcomes and determine the subsequent choice, whether to repeat or switch, for action selection.
- frontal corticostriatal ensemble
- outcome-based decision-making
- single-trial analysis
- tensor component analysis
Significance Statement
We analyzed neuronal ensemble activity simultaneously recorded in the secondary motor cortex (M2) and dorsal striatum (DS) to show how M2-DS circuit mediates decision-relevant signal through cross-regional transmission in decision-making. Decomposed cross-regional neural dynamics exhibited distinct characteristics related to choice position, switch/repetitive choice, and outcome of action at various points in time within trial. These results indicate M2-DS ensemble collaboratively process multiplicate decision-related signals.
Introduction
Animals can select an appropriate action based on sensory cues and outcomes of past actions to adapt flexibly to changing circumstances (Dolan and Dayan, 2013). The underlying neural circuit is widely believed to involve the frontal cortex-basal ganglia circuit (Ragozzino, 2007; Balleine and O’Doherty, 2010; Hikosaka and Isoda, 2010). However, it remains unclear how population neuronal dynamics interact between the frontal cortex and basal ganglia for the adaptive choice behavior. Synchronous neuronal activity across the frontal cortex and downstream subcortical striatum is correlated with skilled motor learning (Koralek et al., 2013; Lemke et al., 2019) and flexible behaviors based on different rules in a T-maze (Oberto et al., 2022), as well as outcomes of action (Handa et al., 2021). Neural trajectories of the secondary motor cortex (M2), a part of the medial frontal cortex, and dorsal striatum (DS) concomitantly represent choice and outcome information during an outcome-based two-alternative choice task. Precise spike synchrony between M2 and DS neurons becomes more prominent during periods of improved task performance (Handa et al., 2021), suggesting cross-regional coactivation of neuronal population correlated with behavioral variables. However, little is known about how cross-regional neuronal dynamics (i.e., M2-DS ensemble) contribute to decision-making during adaptive outcome-based action selection.
Large-scale neuronal activity can be analyzed by adopting dimensionality reduction methods to profile the neuronal ensemble activity (Cunningham and Yu, 2014). Furthermore, conducting trial-by-trial analyses of ensemble activity provides a means to assess the relationship between alterations in the physiological aspects of neuronal ensemble activity and emerging behavioral variables within a single session (Quian Quiroga and Panzeri, 2009; Musall et al., 2019; Veuthey et al., 2020). Motivational states of an animal regulate behavior and such internal states in the brain can be altered over trials within a single session (Burton et al., 1976; Berridge, 2004; Allen et al., 2019). For example, in the initial trials, a thirsty animal is notably motivated to engage in a behavioral task to attain a drop of water as a reward, whereas its motivated performance may diminish in the later trials. Neural representation in the brain can be altered in accordance with such behavioral changes (Allen et al., 2019). We wonder whether the M2-DS ensemble adaptively processes decision-related information in correlation with behavioral variables during a single behavioral session.
To address these questions, we carried out simultaneous electrophysiological recordings in the M2 and DS while rats performed an outcome-based two-alternative choice task. We analyzed cross-region ensemble spike activity at the single-trial level by applying tensor component analysis (TCA) to the concatenated ensemble activity of M2 and DS neurons. TCA provides a crucial advantage over dimensional reduction approaches, allowing us to quantify trial-to-trial variations in neural activity. Additionally, because TCA is an unsupervised method, unexpected features of ensemble activity can be unveiled (Williams et al., 2018). We observed that choice-position–selective neural dynamics were altered over trials, and such trial-by-trial alterations were correlated with trial-basis behavioral fluctuations. Certain neural dynamics revealed differential states between choice patterns, repetitive choices, and switch choices, even when the incoming motor responses and outcomes were the same. Other neural dynamics could discriminate between the rewarded and unrewarded outcomes following action selection. Choice-pattern–selective within-trial activity differed temporally from outcome-selective within-trial activity. These results suggest that trial-basis fluctuations in M2-DS ensembles could be attributed to behavioral variables and that the M2-DS ensemble was implemented in outcome monitoring and continued decision-making for action selection.
Materials and Methods
Animal preparation
All animal procedures were performed in accordance with the Animal Experiment Plan of the Animal Experiment Committee of RIKEN (approval number: H25-2-234[1]). The multiunit recording results during task performance presented in this study were obtained from the reanalysis of previously collected behavioral and electrophysiological data (Handa et al., 2021). Male Long–Evans rats (N = 14, 6 weeks, 200–220 g, Japan SLC) were employed. Home cages were situated in a temperature- and humidity-controlled environment with lights maintained on a 12 h light/dark cycle.
Stereotaxic surgery
All surgical procedures were performed under sterile conditions. Rats were anesthetized with 2% isoflurane, and their body temperature was monitored with a rectal probe and maintained at 37°C on a heating pad during the surgery. A sliding head attachment (Narishige) was implanted in the skull using a dual-curing resin cement (Panavia, Kuraray Noritake Dental) and dental resin (Unifast II, GC). Reference and grounding electrodes (Teflon-coated silver wire, A-M Systems) were positioned over the dura mater above the cerebellum. Following recovery from surgery, rats were deprived of water in their home cages, with water used as a reward for behavioral task execution; however, food was available ad libitum. Rats obtained ∼10 ml of water at the task chamber when they engaged in the task performance, whereas they were supplied 10 ml of water at the cage when the behavioral experiment was not performed. To confirm recording sites within regions showing corticostriatal projections from the M2 to DS, a retrograde tracer, Fluoro-Gold (FG; Fluorochrome), was injected into the DS 3 d before the electrophysiological recording experiment. A glass micropipette filled with 2% FG dissolved in 0.1 M cacodylic acid was installed on a micromanipulator angled medially by 27°. The pipette was inserted through a small burr hole drilled in the skull over the left hemisphere (AP: +1.5 mm to the bregma, ML: 1.0 mm to midline, 4.3 mm traveling distance). The pipette tip reached the dorsocentral part of the striatum (AP: +1.5 mm anterior to the bregma, approximately ML: 3.0 mm to midline, ∼3.8 mm ventral to pia mater) based on previous anatomical evidence (Reep et al., 2003). FG was iontophoretically infused using an iontophoresis pump (BAB-501; Kation Scientific). After the completion of training sessions, two cranial windows were created above the DS and M2 of the left hemisphere (AP: +1.0 and +3.0 mm to the bregma, ML: 3.0 and 1.0 mm to midline for DS and M2, respectively), and their dura maters were removed for electrophysiological recordings.
Behavioral task
Rats underwent training to perform an outcome-based two-choice task using a customized multiple-rat training system (O’Hara & Co.), facilitating parallel learning of the task paradigm for multiple rats simultaneously (Handa et al., 2021). The behavioral task was controlled using a custom-written software in LabVIEW (National Instruments). Individual rats were secured in a body-supporting cylinder, and their heads were rigidly and painlessly fixed using a sliding head holder on a stereotaxic frame (Fig. 1A). Spouts were linked to a syringe on a single-syringe pump (AL-1000; World Precision Instruments) using silicon tubing. Water delivery from each spout was regulated by a pinch valve and syringe pump triggered by the TTL signal. The trial initiated with a pure tone presentation (3 kHz, 1 s, 60 dB SPL; Fig. 1A, “Start”). Rats were instructed to refrain from licking any spouts from the start cue until the appearance of another auditory cue (10 kHz, 0.2 s, 60 dB SPL; “Go”). If rats licked any spouts during the delay period (“1st Delay”), the trial was promptly aborted. The pseudorandom delay period ranged from 0.7 to 2.3 s. Following the onset of the Go cue, rats could lick either the left or right spouts within a response window (5 s). The first lick was considered as a choice response (“Choice”). If the chosen spout location aligned with the ongoing reward location, 0.1% saccharin water was delivered as a reward after a pseudorandom delay period ranging between 0.3 and 0.7 s (“2nd Delay”). The next trial commenced after an outcome period (4 s, “Outcome”). However, when rats chose a no-reward spout, they received no sensory feedback but had an additional time of 5 s after the outcome period. The subsequent trial began after a timeout. Once the accumulated total number of rewarded trials reached 10 within each block, the reward-associated spout position reversed without any feedback, such as sensory or physical differences in the task. Block reversal occurred after ∼11–12 trials if the rats frequently repeated the rewarded choice. On the other hand, if the rats frequently selected the incorrect spout, the block reversal occurred though >11–12 trials because the block reversal occurred when the rat achieved the criterion (10 rewarded choices within the block). Therefore, rats could not anticipate block reversal without experiencing forthcoming trials. We trained all rats over 21 training sessions. If the performance of rat reached an achievement level (the reward acquisition probability of 75%) within the series of training sessions, we defined the remaining sessions as overtraining sessions. Otherwise, there was no overtraining day. After the training sessions, the rats were subject to recording experiments regardless of the learning achievement.
Electrophysiological recordings during task performance
Following the training sessions, each animal underwent two daily recording experiments. Multineuron activity was simultaneously recorded from the M2 and DS of the left hemisphere using two 32-channel silicon probes. These probes consisted of four shanks (0.4 mm shank separation), each featuring tetrode-like electrode sites spaced vertically by 0.5 mm (A4×2-tet-7/5mm-500-400-312, NeuroNexus Technologies). Each probe was connected to a custom-made headstage on one of two fine micromanipulators (1760-61, David Kopf Instruments) mounted on a stereotaxic frame (SR-8N, Narishige). The silicon probe was vertically inserted (depth from the pia mater: 1.2 mm) into M2 (at the center of probe: +3.0–3.6 mm to the bregma, 1.0–1.4 mm to midline), with the shanks aligned along the midline (Fig. 1D). Another silicon probe, angled posteriorly by 6°, was inserted into DS through a cranial window (at the center of probe: +0.6–1.0 mm to the bregma, 2.7–3.1 mm to midline, 4.0 mm traveling distance), with the shanks aligned along the coronal suture (Fig. 1D). Multiunit signals were amplified by the headstages before being fed into main amplifiers (Nihon Kohden) with a bandpass filter (0.5 Hz–10 kHz). All neural data were sampled at 20 kHz using two hard disk recorders (LX-120, TEAC), capturing the time of the task and licking events for each spout (left and right).
Histology
After the recording sessions, rats were deeply anesthetized with urethane (2–3 g/kg, i.p.) and subsequently perfused intracardially with chilled saline followed by 4% paraformaldehyde (PFA) dissolved in 0.1 M phosphate buffer (PB). The fixed brains were stored in 4% PFA overnight and then placed in a 30% sucrose solution in 0.1 M PB for 2 weeks. Postfixed brains were frozen and coronally sliced into 50-μm-thick serial sections using a microtome cryostat (HM500OM, Microm). The brain sections were stored in 0.1 M PB at 4°C overnight. For fluorescent visualization of FG-labeled neurons, brain sections were incubated with an anti-FG antibody from rabbit (AB153, 1:3,000, Millipore) at 4°C overnight. Subsequently, they underwent incubation with goat anti-rabbit IgG conjugated with Alexa Fluor 594 (A11012, 1:500, Invitrogen) for 2 h. Fluorescence images were acquired using a fluorescence microscope (Olympus, AX70) to confirm the presence of FG-labeled neurons around the silicon probe recording locations in the M2 and near the FG injection site in the DS. To verify the silicon probe track, the slices were counterstained with neutral red Nissl. Recording locations in the M2 and DS as well as AP coordinates were determined in accordance with the rat brain atlas (Paxinos and Watson, 2009).
Data analysis
All behavioral and neuronal data were analyzed by custom-written MATLAB scripts (The MathWorks).
Spike sorting, clustering, and refining
For each tetrode, spikes were isolated from multiunit activity by a custom-made semiautomatic spike-sorting program EToS (12 feature dimensions for four channels; high-pass filter at 300 Hz; time resolution at 20 kHz; spike-detection interval >0.5 ms; Takekawa et al., 2010, 2012). The sorted spike clusters were combined, divided, and discarded manually to refine single-neuron clusters by Klusters (Hazan et al., 2006). To avoid overlapping of detection of same units recorded between distinct tetrodes, we checked cross-correlations of spike times among isolated units across all of tetrodes. If there was a high correlation peak only at zero time between a pair of units, one of the units was excluded from further analyses because those spikes which originated from the same neuron were presumably recorded through different tetrodes.
Dataset
We reanalyzed behavioral and multineuron spike datasets previously recorded, sorting them in 14 recording sessions with 12 rats using a new analytical approach (Handa et al., 2021). Dataset for this study is summarized in Table 1, where the experimental ID, number of overtraining sessions, reward acquisition probability (the fraction of correct choices) in the recording session, and number of M2 and DS units are listed.
Tensor component analysis
To generate an original neural data tensor X for each recording session, we computed trial-based perievent time histograms aligned at the choice response, using a 200 ms sliding window with a 50 ms step for individual M2 and DS neurons. To obtain substantial firing activity in population of neurons, we chose the 200 ms bin for the sliding window because the signal-to-noise ratio in characteristics of TCA components using 200 ms bin was better than that analyzed using 50 ms bin for the sample size of our data. Subsequently, we obtained a third-order tensor (neuron, time, and trial) for the M2-DS ensemble through TCA (Fig. 2A). To investigate the single-trial dynamics of the M2-DS ensemble activity, we applied TCA to tensor X using a MATLAB-based toolbox (Tensor Toolbox for MATLAB, version 3.2.1, https://www.tensortoolbox.org/; Bader and Kolda, 2008). In this analysis, the M2-DS ensemble activity was decomposed into a third-order tensor
Choice-position selectivity while trials progress
We assessed significant differences in neural dynamics between left and right choice trials in a given session by applying a two-sample t test (p < 0.05) to trial factors for each TCA component. If the differences were statistically significant, we referred to the TCA component as “choice-position” selective TCA component. We computed the discrimination degree d’ as a choice-position selectivity index in a window of 30 trials by sliding the window with one trial step. Instantaneous d’ was computed as follows:
Correlation of choice-position selectivity with behavioral variables
To investigate whether the change in choice-position selectivity of the TCA component was correlated with the change in behavioral variables, we calculated two behavioral variables: reaction times (RTs) and the number of licks after the choice response. In each trial, RT was calculated as the duration between Go cue onset and the time at which the first lick occurred after Go cue onset and the number of licks that emerged within 2 s after the choice response. We computed the average RTs and number of licks by sliding the analysis window of 30 trials in one trial step. We subsequently computed Pearson's correlation between the trial series of choice-position selectivity d’ and each behavioral variable. If the p value was <0.01, we defined the TCA component as significantly correlated with the behavioral variable.
Trial series of intertrial intervals
We measured a duration between the lick choice at trial t and the onset of Start cue onset at the trial t + 1, in which the rat engaged in task performance (lick choice) after Go cue onset. We defined this duration as an intertrial interval (ITI). If the rat did not make the lick choice at the trial t + 1, we repeated the same trial until the animal made response and thus the ITI got longer. The ITIs were averaged by sliding the analysis window of 30 trials in one trial step as treated for trial series of RTs and the number of licks above.
Quantification of choice-pattern selectivity in M2-DS combined ensemble activity
To quantify the extent to which the neural dynamics in switch choice trials differed from those in repetitive-choice trials, we computed the standard score (Z-score) of trial factors based on their mean and SD within individual blocks. In the case of the left block, during which the left spout was associated with reward delivery, we calculated the Z-score of trial factors in individual left blocks, ranging from the trial where the rat switched its choice from the preceding unrewarded choice (switching from right unrewarded choice to left rewarded choice) to the last trial before the reversal of reward position. The Z-scores in the switch choice trial (Switch) and in the three continued repetitive-choice trials (Rep-1, Rep-2, and Rep-3) were computed as follows:
Quantification of outcome selectivity in M2-DS combined ensemble activity
To quantify the level of trial factor differences between rewarded and unrewarded choice trials in M2-DS combined ensemble activity, Z-scores of trial factors were computed, similar to the aforementioned analysis of choice-pattern selectivity. For the left block, trial factors were collected in a sequence of trials beginning from a trial where the rat correctly switched its choice from the preceding unrewarded choice (right to left) up to one trial before the initial switch choice trial (left to right) in the next block. This chunk of trials encompassed both repetitive-choice but unrewarded trials, occurring due to the reversal of the reward position. Z-scores were computed for the initial three repetitive choice and reward trials (Rep-1 & rwd, Rep-2 & rwd, and Rep-3 & rwd), along with those of the repetitive choices and unrewarded trials (Rep & unrwd), as described earlier. If any unrewarded choice trials were interspersed before achieving the four consecutive rewarded choices, the block was excluded from the analysis. Assuming that Rep-1 and rwd, Rep-2 and rwd, and Rep-3 & rwd and Rep and unrwd represented the same choice pattern (repetitive choice) but different outcome conditions (Rep-1 & rwd, Rep-2 & rwd, and Rep-3 & rwd: rewarded; Rep & unrwd: unrewarded). Statistical significance was determined through one-way ANOVA among the four choices, followed by post hoc multiple-comparisons Dunnett test (p < 0.05). If all three pairs (Rep & unrwd vs Rep-1 & rwd, Rep & unrwd vs Rep-2 & rwd, and Rep & unrwd vs Rep-3 & rwd) exhibited significance, the TCA component was categorized as “outcome”-selective TCA component. A parallel computation was conducted for right choice trials separately.
Quantification of choice-pattern and outcome selectivity in M2- and DS-alone ensemble activities
TCA was applied separately to the M2- and DS-alone ensemble activities using the same dataset as the M2-DS ensemble activity. Z-scores were computed for the analyses of choice-pattern and outcome selectivity, following the previously described methods. To compare the TCA components based on the M2- and DS-alone ensemble with those based on M2-DS ensemble, we randomly resampled the equivalent number of neurons in the M2-DS ensemble relative to the number of neurons in M2- and DS-alone ensemble, respectively. TCA was applied to the M2-DS ensemble to get TCA component. The number of TCA components and their Z-scores were compared with the TCA data based on M2-/DS-alone ensemble. When we compared Z-scores between M2-DS ensemble and M2-/DS-alone ensemble, we used two-sample t test to obtain t-statistics. We repeated this procedure 10 times. If t statistic was positively larger than 0, the Z-score of the M2-DS ensemble was larger than that of the M2-/DS-alone ensemble and vice versa.
Quantification of regional contribution in M2-DS combined ensemble activity
To quantify the contribution of DS and M2 to the TCA component based on M2-DS combined activity, we computed the regional contribution index using neuron factor. The values of each neuron factor were separated into values for M2 and DS cells, and their root mean square (RMS) was computed for each region. The regional contribution index C was computed as follows.
Results
Rats exhibit retrospective outcome-based choices
Head-restrained rats engaged in an outcome-based two-choice task, selecting between two spouts to make their choice (Fig. 1A). A reward was received when the chosen spout matched the ongoing reward-spout position. The action–outcome association was systematically reversed without sensory feedback after accumulating 10 rewarded trials in each block. We analyzed the behavioral data from 14 recording sessions with 12 rats. Across all sessions, the average number of trials per session was 589 ± 124 (mean ± SD). Rats discerned the reversal of choice-reward contingency within subsequent trials, responding to the experience of no-reward events, time-out, and reward acquisition (Fig. 1B). The averaged trial number per block was 12.5 ± 0.650 (mean ± SD). In both the left and right choice trials, rats repeatedly selected the same spout as that in the preceding rewarded trial but switched their choice following one to several unrewarded trials (Fig. 1B,C; one-way ANOVA: p < 10−21, post hoc Tukey–Kramer test: last rewarded trials vs last unrewarded trials, left choice, p < 0.001, right choice, p < 0.001). No statistically significant differences were observed in the repetitive choice probabilities between the left and right choice trials regardless of the outcome conditions in the preceding trial (post hoc Tukey–Kramer test: last rewarded trials, left choice vs right choice, p = 0.973; last unrewarded trials, left choice vs right choice, p = 0.201). The rats demonstrated an inability to predict the reversal of the reward block, choosing nearly no reward-associated spouts in some trials postreversal. In response to the block reversal, rats then retrospectively switched to another choice (Fig. 1C, right). In essence, the choice pattern of the rats closely mirrored the win-stay and lose-shift strategy, an optimal approach for maximizing reward acquisition in the current task, considering that the rats did not anticipate the block reversal.
Simultaneous recording of ensemble neural activity from M2 and DS
To investigate whether frontal corticostriatal ensembles are encoded at the single-trial level, we conducted simultaneous recordings of multineuron activity in both the M2 and DS of the left hemisphere using two multielectrode probes (Fig. 1D). Across 14 recording sessions, a considerable number of units were recorded in both M2 (mean ± SD, 37.3 ± 16.5 units) and DS (38.8 ± 15.5 units). The rodent M2 (or the rostral agranular medial cortex) is one of the primary cortical areas projecting to the dorsocentral region of the striatum (Cheatwood et al., 2003; Reep et al., 2003; Hintiryan et al., 2016). To confirm the presence of such corticostriatal projections from the recording site in M2 to that in the DS, we iontophoretically infused a retrograde tracer FG into the central part of the DS before recording sessions (Materials and Methods). FG-labeled corticostriatal neurons were primarily observed in layers 3 and 5 of the M2 (Fig. 1E). The probe track for M2 and DS recordings aligned with the FG-labeled region in the M2 or near the injection site in the DS (Fig. 1D,E), confirming the recording of the multineuron activity from the directly connected subregions of the M2 and DS.
Trial-by-trial changes in choice-position–selective activity of M2-DS ensembles are correlated with behavioral variables
To investigate the dynamic changes in the characteristics of M2-DS ensembles across trials, we utilized TCA (Williams et al., 2018) for both dimensionality reduction of high-dimensional neuronal activity and quantification of trial-to-trial fluctuations in ensemble activity (Materials and Methods). Here, we refer to the collective firing activity of multiple M2 and DS neurons as the M2-DS ensemble. TCA decomposes the M2-DS ensemble activity into a third order tensor (Fig. 2A). Each component comprises three vectors, including a “neuron factor” vector (Fig. 2B), a “temporal factor” vector (Fig. 2C), and a “trial factor” vector (Fig. 2D).
We hypothesized that the TCA components of the M2-DS ensemble would exhibit choice-position selectivity that changes across trials, consistent with our previous findings of choice-position selectivity in single-unit and population activity in both the M2 and DS during this task (Handa et al., 2021). Indeed, the trial factor of this TCA component revealed significant differences between the left and right choices (two-sample t test; p < 10−126). Additionally, we observed that the difference in trial factor between choice positions was small in the initial trials but gradually increased in the middle and later parts of this session (Fig. 2D). To quantify the gradual changes in trial factors across trials, we calculated the degree of discrimination across trials as choice-position selectivity (Materials and Methods). The discrimination degree increased across trials and significantly differed from the discrimination degree derived by shuffling trial order (permutation test; p < 0.01; Materials and Methods; Fig. 2E).
We examined whether the gradual changes in the discrimination degree were associated with alterations in behavioral variables across trials. Specifically, we investigated the relationship between the discrimination degree and two behavioral variables: reaction time (RTs) and the number of licks after the choice response. The trial series of the number of licks exhibited a higher correlation with the degree of discrimination of the TCA component (Pearson's correlation: r = −0.60; p < 0.0001) than with the trial series of RTs (r = −0.17; p < 0.0001; Fig. 2F). Although RTs (Fig. 2G1) and the number of licks (Fig. 2G2) varied across the 14 sessions, an average reduction in the number of licks was observed across trials. The slope of linear regression model, which was fitted to the trial series of the number of licks for each session, was negative (slope, mean ± SD = −0.00192 ± 0.00101) and significantly different from 0 (t test; p = 8.20 × 10−6). Among the 166 choice-position–selective TCA components, 112 (67.4%) and 125 (75.3%) were significantly correlated with RTs and the number of licks, respectively (Pearson's correlation; p < 0.01; Fig. 2G1,G2, right). Among these significantly correlated TCA components, 87 revealed a significant correlation with both RTs and the number of licks. The magnitudes of the significant correlation coefficients were larger in the correlation with the number of licks than in the correlation with RTs (two-sample t test; p < 0.001; Fig. 2H). These behavior-correlated TCA components were observed in all recording sessions. Their correlation coefficients were not significantly correlated with the number of overtraining sessions on this task (Pearson's correlation test: correlation with RTs: r = −0.160, p = 0.583; correlation with the number of licks: r = 0.024, p = 0.933), suggesting that the properties of trial factors in the choice-position–selective TCA were not affected by the overtraining for this task.
To quantify how much each brain region contributed to the choice-position–selective TCA components which correlated to RTs or the number of licks, we computed the regional contribution index using neuron factor (Materials and Methods). The regional contribution indices broadly ranged (Fig. 2I1,I2, respectively). On average, the contribution index was significantly shifted to positive values, suggesting the DS was more of a contributor than the M2 (t test: TCA components correlated to RTs: p = 2.81 × 10−6, TCA components correlated to the number of licks: p = 1.03 × 10−5).
In contrast to the reduction tendency of number of licks at the late trials in the sessions, we found that ITIs, which were intervals between consecutive trials where animals were engaged in lick response post Go cue, got longer at the late trials in the sessions (Fig. 2J; slope, mean ± SD = 0.00659 ± 0.00806, t test: p = 0.00912). If the rat did not respond to Go cue within the response window, the same trial was repeated until the rat made choice, resulting in longer ITIs. This increasing tendency of ITIs at the late trials indicates the reduction of motivation to engage in task performance. Thus, the ITIs could be a kind of behavioral variable to quantify how constantly animals were motivated to engage in choice action to obtain a reward. The number of licks were negatively correlated with the ITIs (Pearson's correlation coefficient: mean ± SD = −0.521 ± 0.20751, t test: p = 3.62 × 10−7; Fig. 2J,K), but RTs were not (mean ± SD = 0.0307 ± 0.262, t test: p = 0.676; Fig. 2K). Therefore, the behavioral results suggest that this reduction in the number of licks over the trials may reflect changes in the motivational state of the rats to engage in task performance.
These results indicate that TCA unveils the dynamics of choice-position selectivity in M2-DS ensembles across trials, which may be linked to changes in behavioral variables such as motor preparation and/or motivational state.
TCA unveils activity patterns of M2-DS ensembles distinguishing between repetitive and switch choices
In this study, TCA not only confirmed the anticipated choice-position–selective activity type of M2-DS ensembles, as revealed in our previous analysis (Handa et al., 2021), but also unveiled an unexpected trial factor pattern that differentially altered depending on the choice pattern—switch choice versus repetitive choice (Fig. 3A–C). A representative TCA component exhibited an increase and decrease in within-trial activity (temporal factor) before and after the response, respectively (Fig. 3B). The trial factor displayed distinct patterns between repetitive and switch choices in right choice trials but not in left choice trials (Fig. 3C). In right choice trials, the trial factors in switch trials significantly deviated from those in the repetitive trials (Fig. 3C,D). In contrast, in left choice trials, the trial factor fluctuated similarly in both switch and repetitive choices. To quantify the deviation within a single block of trials, we calculated the Z-scores of trial factors within each block to compare switch choice, first, second, and third repetitive choices, where the movement direction and outcome condition (rewarded) were the same (Fig. 3E, top; Materials and Methods). The Z-scores in switch choice trials significantly differed from those in all other repetitive choice trials in the right (contralateral) choice condition (one-way ANOVA: p < 10−13, post hoc Dunnett test, Switch vs Rep-1: p < 0.0001, Switch vs Rep-2: p < 0.0001, Switch vs Rep-3: p < 0.0001), whereas the Z-score was not significantly different in the left (ipsilateral) choice condition (one-way ANOVA: p = 0.102; Fig. 3E). If the difference in Z-scores between switch and repetitive choices reflected difference in previous choice positions (or previous outcomes), the difference in Z-score should be observed in both left and right choice conditions. Therefore, this lateralized difference between repetitive and switch choices could not be attributed to the differences in the previous choices or outcomes between the choice patterns, suggesting that this difference may reflect rather a lateralized cognitive function for switch or repetitive choices. Choice-pattern–selective TCA components were more frequently observed in right (contralateral) choice trials (N = 21) than in left (ipsilateral) choice trials (N = 13). Four TCA components displayed choice-pattern–selective activity in both left and right choice trials. The right choice-preferred TCA component (11 sessions) was detected in more recording sessions than the left choice-preferred TCA component (six sessions; Fig. 3F). The TCA component revealing choice-pattern selectivity in both left and right choices was observed in three sessions. The population of Z-scores for the choice-pattern–selective TCA trial factor was significantly higher in switch choice trials than in continued repetitive choice trials (Fig. 3G). The regional contribution indices broadly ranged (Fig. 3H). On average, the contribution index was not different from 0 (t test: p = 0.135 for left choice block, p = 0.0911 for right choice block), suggesting the M2 and DS equally contributed to the choice-pattern–selective TCA components.
This result suggests that the M2-DS ensembles differentially encode incoming choice information between switch choice (left to right) and repetitive choice (right to right) at the single-trial level.
TCA component exhibits a differential magnitude depending on rewarded and unrewarded choices
In our previous study, both M2 and DS ensembles displayed outcome-related activity (Handa et al., 2021). As expected, the TCA of the M2-DS ensemble revealed a change in trial factors depending on the outcome (rewarded and unrewarded events) at the single-trial level (Fig. 4A–C). A representative TCA component showcased an increase in within-trial activity (temporal factor) following the response (Fig. 4B). The trial factor highly deviated in unrewarded choice trials without a bias of laterality (Fig. 4C,D; left choice: one-way ANOVA: p < 10−18, post hoc Dunnett test, Rep-1 & rwd vs Rep & unrwd: p < 0.0001, Rep-2 & rwd vs Rep & unrwd: p < 0.0001, Rep-3 & rwd vs Rep & unrwd: p < 0.0001; right choice: one-way ANOVA: p < 10−27, post hoc Dunnett test, Rep-1 & rwd vs Rep & unrwd: p < 0.0001, Rep-2 & rwd vs Rep & unrwd: p < 0.0001, Rep-3 & rwd vs Rep & unrwd: p < 0.0001). Outcome-selective TCA components were frequently observed in both left and right choice trials (Fig. 4E), in contrast to the choice-pattern–selective TCA component discussed earlier. The population of Z-scores for the outcome-selective TCA trial factor was significantly higher in the repetitive but unrewarded choice trials than in the repetitive and rewarded choice trials (Fig. 4F). The regional contribution indices broadly ranged in both left and right choice blocks (Fig. 4G). The mean contribution index was not different from 0 (t test: p = 0.138 for left choice block, p = 0.0700 for right choice block), suggesting the M2 and DS equally contributed to the outcome-selective TCA components.
This finding suggests that the M2-DS ensembles differentially encoded outcome information in a given trial.
Differential within-trial activity in choice-pattern–selective and outcome-selective TCA components
The functionally distinct TCA components likely reflect characteristics revealing different roles of the M2-DS ensemble in adaptive outcome-based decision-making. To address this question, we examined the differences in within-trial activity (temporal factor) between the choice-pattern–selective and outcome-selective TCA components. Most choice-pattern–selective TCA temporal factors were activated before the response and exhibited decreased activity after the response (Fig. 5A). Conversely, >50% of the outcome-selective TCA temporal factors showed an increase in ensemble activity following the response (Fig. 5B). On average, these functionally distinct TCA components displayed opposite trends in temporal factors before and after the response (Fig. 5C). We quantified these temporal differences between choice-pattern–selective and outcome-selective TCA components by comparing the peak time of the temporal factor. The peak time of the temporal factor of choice-pattern–selective TCA components (mean ± SD: left choice trials, −0.473 ± 1.78 s, right choice trials, −0.608 ± 1.51 s) was significantly earlier than that of outcome-selective TCA components (left choice trials, 0.696 ± 1.61 s, right choice trials, 0.525 ± 1.51 s; Mann–Whitney U test, left choice trials: p = 0.0189, right choice trials: p = 0.00401; Fig. 5D).
This finding suggests that the M2-DS ensemble plays a temporally distinct role in adaptive choice behavior by detecting outcomes after the choice response and flexibly making action decisions (such as repeating or switching choice) based on the outcome information.
Comparison of TCA components based on M2- and DS-alone ensembles with TCA components based on M2-DS ensemble
If the features of ensemble activity of M2 and DS neurons were heterogeneous, or spike activity was not well coordinated between M2 and DS, the application of TCA to the concatenated ensemble activity of M2 and DS neurons (the M2-DS ensemble) could potentially attenuate task-related signals compared with TCA on ensemble activity in each region separately (Runyan et al., 2017). To explore this possibility, we investigated whether the choice-pattern–selective and outcome-selective TCA components derived from the M2-DS ensemble provided worse task-related signals than those obtained from ensembles in M2 or DS alone. TCA was separately applied to the M2-alone and DS-alone ensembles over the same data used for the M2-DS ensemble analysis. Choice-pattern–selective and outcome-selective TCA components were observed for both M2-alone and DS-alone ensembles. To control the number of neurons between M2-DS combined ensemble and M2-alone/DS-alone ensemble, we randomly resampled the equivalent number of neurons in the M2-DS ensemble relative to the number of neurons in M2-alone and DS-alone ensemble, respectively (Materials and Methods).
Regarding choice-pattern–selective TCA components, in left (ipsilateral) choice trials, the total number of detected TCA components was significantly fewer in the M2-alone ensemble (n = 4; t test: p = 4.32 × 10−5; Fig. 6A1), whereas that of the DS-alone ensemble (n = 15) was more than that of the M2-DS ensemble (t test: p = 0.0334; Fig. 6B1). In right (contralateral) choice trials, the number of choice-pattern–selective TCA components were comparable with that of the M2 (n = 15) and DS (n = 14) ensembles alone (t test: p = 0.244 for M2, p = 0.421 for DS; Fig. 6C1,D1). To compare the selectivity strength in the choice-pattern–selective TCA components between M2-DS ensemble and M2-/DS-alone ensembles, we compared the absolute Z-scores of the trial factors in switch choice trials between the M2-DS and M2-alone ensembles and between the M2-DS and DS-alone ensembles by two-sample t test to obtain t statistic. These values exhibited no significant difference by comparison with M2-alone ensemble (Fig. 6A2,C2). However, we found that the Z-scores of the trial factor were larger in M2-DS ensemble than those in DS-alone ensemble in right choice blocks although the statistical significance was detected in 4 out of 10 cases (Fig. 6D2, black symbols), but not in the left choice blocks (Fig. 6B2).
Regarding outcome-selective TCA components, the total number of TCA components based on the M2-DS ensemble were comparable with that based on the M2-alone ensemble (n = 43) in the left choice block (t test: p = 0.938; Fig. 6E1), whereas they were fewer than that based on the M2-alone ensemble (n = 55) in right choice block (t test: p = 0.0198; Fig. 6G1). Regarding the comparison with DS-alone ensembles, the number of TCA components in M2-DS ensemble were more than that in DS-alone ensemble (left/right choice: n = 40/39) in both blocks (t test: p = 0.00279 for left choice block; Fig. 6F1; p = 0.000103 for right choice block; Fig. 6H1). The selectivity strengths were not different between M2-DS ensemble and M2-/DS-alone ensembles in left and right choice blocks (Fig. 6E2,G2,F2,H2).
This finding indicates that TCA using the combined M2 and DS ensembles did not attenuate task-related signals in most cases; instead, it may not only reflect the contribution of either the M2 or DS ensemble but also the cooperative contribution of M2 and DS ensemble activity.
Discussion
In this study, we demonstrated distinct representations of decomposed M2-DS ensemble activity for various task-relevant behaviors at the trial level in rats performing an outcome-based choice task using the TCA approach. TCA is an unsupervised method and, in contrast to linear discriminant analysis and linear regression, it automatically identifies choice-position, choice-pattern, and outcome-specific patterns. Choice-position–selective TCA components (the choice-position–specific spatiotemporal neural dynamics) revealed dynamic changes in selectivity over trials, which were correlated with changes in behavioral variables, such as RTs, and/or the number of licks after the response. The TCA revealed neural dynamics showing selectivity for choice patterns (repeat and switch choices), even when the choice position and outcome were identical. Choice-pattern–selective and outcome-selective neural dynamics revealed functionally distinct within-trial activity. Choice-pattern–selective within-trial activity increased activity earlier (before the choice response), whereas outcome-selective within-trial activity increased activity later (after the choice response). TCA application on M2-DS ensemble activity tended to yield more task-related neural dynamics than TCA application on M2 or DS alone.
The rodent M2 integrates sensory and outcome information to make decisions regarding motor planning (Barthas and Kwan, 2017). The activity of M2 neurons represents laterality in the body (Erlich et al., 2011; Sul et al., 2011), forelimb (Soma et al., 2017), and tongue movements (Siniscalchi et al., 2016; Handa et al., 2017, 2021; Kurikawa et al., 2018). Additionally, the M2 is implicated in outcome evaluation for action (Sul et al., 2011; Gremel and Costa, 2013; Handa et al., 2021) and in the neuronal encoding of outcomes, such as reward and nonreward events (Sul et al., 2011; Handa et al., 2021). The rodent DS serves a significant gateway of the basal ganglia, receiving major synaptic inputs from the cerebral cortex (Kincaid et al., 1998; Cheatwood et al., 2003; Reep et al., 2003; Wall et al., 2013; Hintiryan et al., 2016). Neuronal activity in DS similarly shows the laterality of the choice response (Kim et al., 2009; Handa et al., 2021) and the value of the outcome (Kim et al., 2009; Nonomura et al., 2018). Therefore, similar behavior-related signals are observed in the two regions. To gain a deeper understanding of the mechanisms underlying cross-region transmission between M2 and DS in outcome-based decision-making, a previous study used Fisher's linear discriminant (FLD) for dimensionality reduction of ensemble activity. The dynamic characteristics of the ensemble activity in each region were detected, and these characteristics were compared between M2 and DS. Neural trajectories in the M2 and DS revealed temporally similar task-related signals, such as choice position and outcome. Precise spike synchrony between M2 and DS neurons emerges more frequently when task performance is superior (Handa et al., 2021). Consistent with this result, another study demonstrated that cortical representation is topographically reflected in the striatal subregion (Peters et al., 2021). Our findings demonstrate that the M2-DS ensemble activity can be deconstructed into distinct functions—specifically, choice position, choice pattern, and outcome—without compromising the representation of decision-related signals in each region. This evidence suggests that the ensemble activity of interconnected regions (the M2-DS ensemble) effectively processes similar decision-related signals in a cooperative manner for further processing through downstream structures in the basal ganglia for motor selection.
Although the previous study successfully identified temporally parallel processing of information related to choice position and outcome between the M2 and DS using FLD, this approach encountered challenges in identifying choice-pattern–selective activity at the population activity level (Handa et al., 2021). Similar to conventional principal component analysis, FLD utilizes trial-averaged data to compute a vector of hyperplanes, maximizing the degree of discrimination between two conditions (left and right choices; Bishop, 2006). The assumption of trial averaging is that trial-by-trial variability is a task-irrelevant noise. Conversely, TCA considers trial-by-trial variability to extract features of population activity at a single-trial level in an unsupervised manner (Williams et al., 2018). Our current results not only demonstrated choice-position and outcome-selective M2-DS ensemble activity but also revealed choice-pattern–selective activity based on trial-by-trial analysis. These choice-pattern–selective neural dynamics were observed even in neural dynamics based on M2 or DS ensemble activity alone, suggesting that the choice pattern is commonly processed in both M2 and DS. In this case, dimensionality reduction in an unsupervised manner enables the uncovering of the latent features of the neural ensemble. Choice-pattern–selective neural dynamics of M2-DS ensemble activity were found more frequently in contralateral choice trials than in ipsilateral choice trials. In the contralateral choice trials, the selectivity strength of choice-pattern–selective TCA neural dynamics was larger in M2-DS ensembles than that in DS-alone ensemble. This result may reflect the cooperative contribution of M2-DS ensembles to the lateralized function. This intriguing component of the M2-DS ensemble activity may reflect cognitive features for action selection based on outcome rather than motor commands as the incoming choice position is identical, but the difference lies in selecting the spout position by deciding whether to repeat the previous choice or switch from the preceding choice in the outcome-based choice task. In support of this interpretation, previous studies have indicated that M2 is involved in cognitive switching between sensory-guided and automated choice rules (Siniscalchi et al., 2016), as well as flexible visual categorization (Wang et al., 2020). The DS is additionally implicated in the processing of action selection based on outcome probability or value (Nonomura et al., 2018; Cox and Witten, 2019). Our findings are consistent with these neuronal functions in both the M2 and DS, extending to M2-DS ensembles.
Choice-position–selective TCA neural dynamics revealed trial-to-trial fluctuations in choice-position selectivity, which correlated with variable behavioral parameters. The correlation with the number of licks may indicate the significance of the M2-DS ensemble activity in relation to changes in motivation to participate in the task, given the decrease in the number of licks observed in the later period of the session. Our finding suggests that the DS could be more of a contributor in the M2-DS circuit activity for the choice-position–selective TCA neural dynamics, whereas the M2 and DS contribute equally to the choice-pattern–selective and outcome-selective TCA dynamics. This newly observed result is also attributed to the trial-based TCA approach although we could not clarify the internal correlation between M2 and DS within M2-DS ensemble activity. The analysis of simultaneously recorded population activity is useful to interpret neural function related to behavioral or cognitive variables at the single-trial level (Kiani et al., 2014). Recent study demonstrates the relationship between trial-wise variability in choices and variability in value signals, which are related to decision-making, decoded from neuronal population activity in the nonhuman primate orbitofrontal cortex (McGinty and Lupkin, 2023). For the rodents, simultaneous recording of spiking activity of population of neurons from multiple brain regions could be applicable thanks to the recent technical advance in high-density electrode probes (Steinmetz et al., 2021). The read-out of neuronal functions across different brain regions by using such large-scale spike activity is important to interpret the functions of cross-regional neural circuits by means of the approach of dimensionality reduction together with cross-regional correlation analysis at the single trial level (Veuthey et al., 2020; Gokcen et al., 2022; Kondapavulur et al., 2022) in future studies.
A recent study demonstrated that synchronous spike activations across regions, including the medial prefrontal cortex and the downstream dorsomedial and ventral striatum, emerge with behavioral correlations contingent on task demands during the T-maze task. Different firing assemblies were observed at various times within the task (Oberto et al., 2022). In our study, the within-trial activity revealed distinct characteristics, such as choice-pattern–selective neural dynamics and outcome-selective neural dynamics, which varied temporally. This suggests that the M2-DS ensemble serves different functions at different time points. Choice-pattern–selective neural dynamics may be modulated before the choice response, influencing the decision to repeat or switch, whereas the outcome-selective component could be activated after the response to monitor the outcome. In essence, the M2-DS ensemble contributes to temporally distinct roles in adaptive choice behavior by detecting outcomes after choice response and flexibly making decisions of action (repeating or switching choices) based on the outcome information. Growing evidence indicates that frontal cortex-striatal ensemble coactivation emerges at specific times and plays specific roles in ongoing goal-directed behavior.
Footnotes
The authors declare no competing financial interests.
This work was partially supported by Grants-in-Aid for Scientific Research (KAKENHI) from MEXT (nos. 22K06485 to T.H., 23H05476 to T.F., and 20K07716 to T.K.). We express our gratitude to R. Harukuni for her invaluable support in behavioral training and histological analysis.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.