Thalamus and claustrum control parallel layer 1 circuits in retrosplenial cortex
Abstract
The granular retrosplenial cortex (RSG) is critical for both spatial and non-spatial behaviors, but the underlying neural codes remain poorly understood. Here, we use optogenetic circuit mapping in mice to reveal a double dissociation that allows parallel circuits in superficial RSG to process disparate inputs. The anterior thalamus and dorsal subiculum, sources of spatial information, strongly and selectively recruit small low-rheobase (LR) pyramidal cells in RSG. In contrast, neighboring regular-spiking (RS) cells are preferentially controlled by claustral and anterior cingulate inputs, sources of mostly non-spatial information. Precise sublaminar axonal and dendritic arborization within RSG layer 1, in particular, permits this parallel processing. Observed thalamocortical synaptic dynamics enable computational models of LR neurons to compute the speed of head rotation, despite receiving head direction inputs that do not explicitly encode speed. Thus, parallel input streams identify a distinct principal neuronal subtype ideally positioned to support spatial orientation computations in the RSG.
eLife digest
Sitting in your car, about to drive home after a long day at work, you realize you have no idea which way to go: you recognize where you are right now, and you remember the name of the street your house is on, but you cannot figure out how to get there. This spatial disorientation happens to people with damage to a brain region called the retrosplenial cortex, whose role and inner workings remain poorly understood. Recent evidence has shown that this area contains ‘low-rheobase’ neurons which are not seen anywhere else in the brain, but what do these neurons do?
Brennan, Jedrasiak-Cape, Kailasa et al. decided to explore the role of these neurons, focusing on the brain regions they are connected to. Experiments were conducted in mice using optogenetics, a technique that activates neurons using pulses of light. This revealed that brain areas involved in processing information about direction and position preferentially communicate with low-rheobase neurons rather than with nearby, more standard neurons in the retrosplenial cortex. The way these spatial signals are sent to the low-rheobase neurons allows these cells to ‘calculate’ how fast a mouse is turning its head using only information about which direction the mouse is facing. Essentially, this neuron can turn directional compass-like signals into a gyroscope signal that can track both direction and speed of head movement. These unique neurons may therefore be ideally suited to combine information about direction and space, suggesting that they may have evolved specifically to support spatial navigation.
Individuals with Alzheimer’s disease show exactly the same type of spatial disorientation as individuals with direct damage to the retrosplenial cortex. This region is also one of the first to show altered activity in Alzheimer’s disease. Exploring whether these unique retrosplenial neurons and their communication patterns are altered in Alzheimer’s disease models could help to understand and potentially treat this debilitating condition.
Introduction
Activity in the granular retrosplenial cortex (RSG) is correlated with a wide variety of behaviors, including spatial navigation, learning, memory, fear conditioning, imagination, and planning for the future (Alexander et al., 2020; Alexander and Nitz, 2017; Chang et al., 2020; Chrastil, 2018; Hinman et al., 2018; Mao et al., 2018; Miller et al., 2021; Miller et al., 2019). The RSG is also among the most densely connected regions of the brain, integrating inputs from a bevy of cortical and subcortical sources and serving as part of the default mode network (Greicius et al., 2009; Kaboodvand et al., 2018; Liu et al., 2019; Whitesell et al., 2021). While these anatomical connections are well-documented (Van Groen and Wyss, 2003; van Groen and Wyss, 1990; Whitesell et al., 2021; Wyss and Van Groen, 1992), few studies have examined their functional nature (Nitzan et al., 2020; Yamawaki et al., 2019b; Yamawaki et al., 2019a; Yamawaki et al., 2016a), and cell-type specificity has yet to be sufficiently explored. Such knowledge is critical to develop a mechanistic understanding of how the RSG integrates information from multiple sources to carry out its spatial and non-spatial functions.
Cortical regions, including the RSG, receive inputs from one or more thalamic nuclei (Herkenham, 1986; Herkenham, 1980; Jones, 2001; Jones, 1998; Peters, 1979; Van Groen and Wyss, 1995). Thalamic relay cells (Clascá et al., 2012) are grouped into three distinct classes: core, intralaminar, and paralaminar/ventral midline (matrix) nuclei (Hanbery and Jasper, 1953; Herkenham, 1986; Jones, 2001; Jones, 1998; Morison and Dempsey, 1941; Rubio-Garrido et al., 2009). Thalamic matrix nuclei project predominantly to layer 1 of many cortical regions (Herkenham, 1986; Herkenham, 1980), resulting in subtype-specific activation of cortical inhibitory (Anastasiades et al., 2021; Cruikshank et al., 2012; Delevich et al., 2015) and excitatory neurons (Collins et al., 2018; Cruikshank et al., 2012; Guo et al., 2018; Rodriguez-Moreno et al., 2020; Rubio-Garrido et al., 2009; Van der Werf et al., 2002; Yamawaki et al., 2019b). Matrix thalamocortical (TC) inputs from the anterior thalamus, which contains the highest proportion of head direction cells in the brain (Taube and Bassett, 2003), influence layer 5 pyramidal cells in RSG (Yamawaki et al., 2019b). However, potential subtype-specific responses and their functional impact on the encoding of directional information in RSG have yet to be examined.
The cortex, including the RSG, also receives widespread inputs from the densely connected claustrum (Brown et al., 2017; Crick and Koch, 2005; Goll et al., 2015; Jackson et al., 2018; Kim et al., 2016; Narikiyo et al., 2020; Wang et al., 2017; White and Mathur, 2018). As with TC signaling, the claustrum innervates both excitatory (da Costa et al., 2010; Narikiyo et al., 2020) and inhibitory (Jackson et al., 2018; Narikiyo et al., 2020; Salerno et al., 1984) neurons via precise laminar targeting (Wang et al., 2017) and projects to many of the same cortical regions as the thalamus (Burman et al., 2011). Claustrocortical (ClaC) projections are thought to regulate complex functions (Crick and Koch, 2005; Goll et al., 2015; Jackson et al., 2018; Kitanishi and Matsuo, 2017; Narikiyo et al., 2020; Renouard et al., 2015; Smith et al., 2012; White and Mathur, 2018), including the coordination and amplification of correlated signals across cortical regions (Kim et al., 2016; Smith et al., 2012; Smythies et al., 2012). Surprisingly, the RSG is one of few regions that does not return reciprocal connections to the claustrum (Zingg et al., 2018). Despite the unique unidirectional nature of this connection, no studies have investigated the functional nature or subtype-specificity of ClaC inputs to RSG neurons.
Each cortical region also receives inputs from select other cortical areas. The RSG receives direct innervation from the secondary motor cortex (Yamawaki et al., 2016a), contralateral RSG (Sempere-Ferràndez et al., 2018; Wyss et al., 1990), and anterior cingulate cortex (ACC; Van Groen and Wyss, 2003; van Groen and Wyss, 1990), among others. The RSG and ACC, in particular, have been implicated in several shared limbic and cognitive functions, including fear conditioning (Frankland et al., 2004; Han et al., 2003; Yamawaki et al., 2019a; Yamawaki et al., 2019b), but the cellular targets of ACC inputs to RSG have not been identified.
Thalamic, claustral, and cortical inputs all converge in the RSG, but whether they do so via parallel or integrative control of target excitatory neurons remains unknown. Within the RSG, the superficial layers contain two subtypes of principal excitatory pyramidal neurons that have strikingly distinct physiology, morphology, and computational capabilities: the low-rheobase (LR) neuron and the regular-spiking (RS) neuron (Brennan et al., 2020; Kurotani et al., 2013; Yousuf et al., 2020). Here, we use channelrhodopsin (ChR2)-assisted circuit mapping (CRACM) to study the subtype-specific circuits formed by anterior thalamic, claustral, anterior cingulate, and dorsal subicular inputs to the RSG. Our results reveal a double dissociation in the precise organization of afferent axons and principal cell dendrites that can support parallel processing in superficial RSG. Specifically, we find that anterior thalamus and dorsal subiculum preferentially control LR neurons, whereas claustrum and anterior cingulate control RS cells. We show that the synaptic dynamics of thalamic inputs allow for the robust encoding of angular speed by LR cell models, even though these inputs only explicitly contain information about head direction, not speed. Together, our results suggest that LR neurons are ideally positioned to contribute to the initial spatial orientation computations performed in superficial RSG.
Results
Cell-type-specific thalamic control of layer three pyramidal cells in granular retrosplenial cortex
Layer 3 of the granular retrosplenial cortex (RSG) contains two types of principal pyramidal neurons: the low-rheobase (LR) neuron and the regular-spiking (RS) neuron (Brennan et al., 2020; Kurotani et al., 2013; Yousuf et al., 2020). Using whole-cell patch clamp recordings, we again distinguished LR from RS neurons by their unique intrinsic properties, predominantly the low rheobase, high input resistance, lack of adaptation, and narrow spike width (Figure 1A–C, Figure 1—figure supplement 1, Supplementary file 1 - Table 1, and Supplementary file 1 - Table 2; Brennan et al., 2020). We used channelrhodopsin-assisted circuit mapping (CRACM) to interrogate inputs to these two cell types. Consistent with previous studies (Odagiri et al., 2011; Van Groen and Wyss, 2003; van Groen and Wyss, 1990; Yamawaki et al., 2019b), virus injections into the anterodorsal and anteroventral thalamic nuclei (ADAV) resulted in ChR2-eYFP expression in thalamocortical (TC) axons and terminals predominantly in layer 1a (L1a; the most superficial part of layer 1) and in layer 3 (L3) of superficial RSG (Figure 1D). To test whether these thalamocortical inputs target the LR and RS neurons in layer 3, we stimulated the thalamic axons and terminals with 1 ms LED pulses at a 10 Hz frequency for 1 s at high LED power over the cell body of the patched neurons (Figure 1E; see Materials and methods). From resting membrane potentials of approximately −65 mV, LR neurons were strongly driven past spike threshold by these thalamic inputs while RS neurons exhibited far smaller EPSPs (Figure 1F). To quantify the impact of TC inputs on LR and RS spike rates, we next depolarized the cells until they were spiking at 10–25 Hz and delivered a single LED pulse (see Materials and methods). To rule out effects of feed-forward inhibition in these experiments, the GABAA antagonist, picrotoxin (50 µM), was added to the bath. Compared to their baseline activity, this stimulation evoked a significant increase in spiking in LR neurons (p=0.002, Wilcoxon rank sum test) but not RS neurons (p=0.19, Wilcoxon rank sum test; Figure 1G). The LR neurons thus had a significantly larger increase in spike rate evoked by the optical pulse than RS neurons (p=0.0012, Wilcoxon rank sum test; Figure 1H). This suggests that LR neurons receive stronger activation from TC inputs than their neighboring L3 RS neurons.
Precise overlap of anterior thalamic afferents with LR, but not RS, apical dendrites
LR neurons are morphologically distinct from their neighboring RS neurons (Brennan et al., 2020; Kurotani et al., 2013; Yousuf et al., 2020). Given the precise laminar pattern of anterior thalamic inputs to RSG (Figure 1D and Figure 2B), we next investigated the overlap between anterior thalamic afferents and the dendritic morphologies of reconstructed LR (n = 10) and RS (n = 5) cells (see Materials and methods). Apical dendrites of LR neurons were most densely localized to L1a, with basal dendrites remaining confined to L3 (Figure 2C). In contrast, RS apical dendrites rarely entered L1a, instead existing in the deeper subdivisions of layer 1 (L1b/c), with basal dendrites primarily in upper layer 5 (L5; Figure 2D). Projection density analysis showed localization of ADAV axons and terminal arbors in L1a and weakly in L3 (Figure 2C&D). Thus, dendrites of LR, but not RS, cells selectively co-localize within layers 1a and 3 with ADAV TC afferents (Figure 2C&D). This anatomical colocalization of LR, but not RS, dendrites with TC axons coupled with the increase in LR, but not RS, spiking in response to TC optical stimulation indicates that the anterior thalamic-retrosplenial circuit at least partially obeys Peters’ rule, which states that neuronal populations with anatomically overlapping axonal and dendritic arbors are more likely to show functional connectivity Peters, 1979; Peters and Feldman, 1976.
To further verify this, we examined the responses of L3 LR and RS neurons to ADAV input at L1a, where the highest projection density of ADAV inputs exist (see Methods). When delivering 1 ms high LED power light pulses, we indeed found that LR neurons had a significantly larger excitatory postsynaptic potential (EPSP) amplitude in response to the L1a LED stimulation compared to RS cell responses (p=0.0185, Wilcoxon rank sum test; Figure 2E). Notably, in response to layer one stimulation, LR cells also showed significantly larger EPSP amplitudes compared to L5 RS neurons (Figure 2—figure supplement 1), which have been previously shown to respond to TC input (Yamawaki et al., 2019b). This suggests that LR neurons, rather than layer 3 or 5 RS neurons, are the predominant RSG cell type receiving and processing head direction inputs from the anterior thalamus.
To further examine the effects of lamination of anterior thalamic inputs to LR and RS neurons on their resulting responses, we conducted a high-resolution protocol that used minimum thresholded LED power to optically stimulate every 20–30 µm along the neuron’s longitudinal axis (Figure 2F; see Materials and methods). We found that LR neurons have the largest EPSP amplitude when the LED pulse is targeting layer 1a, corresponding with the peak projection density of ADAV axons and terminal arbors at this lamina. LR response amplitude then decreases as the LED moves away from the pia before reaching no response at the layer 1/2 border. As LED stimulation entered layer 3, the LR EPSP amplitude increased, consistent with the increase in projection density of the ADAV axons and terminal arbors, before again decreasing until their response amplitude reached 0 mV in layer 5 (Figure 2F). In contrast, the RS neuron showed much smaller EPSP amplitude with no significant variation in response across the layers (Figure 2F&G). To verify this functional sublaminar correlation between thalamocortical axons and pyramidal cell dendrites, we ran a population analysis comparing the EPSP amplitudes within LR and RS cell groups from all layers relative to L1/2, which represents our minimum projection density of ADAV inputs in the superficial layers. We found that within the LR cell population, LR cells have significantly larger responses at L1a (p=0.00001) and cell body (p=0.0006) compared to L1/2 (paired t-test; Figure 2G). This lamination of functional responses in the LR population precisely aligns with the lamination of ADAV axons and terminal arbors within these superficial layers. When stimulating in layer 5, there was no significant difference between EPSP amplitude at L1/2 and L5sup (p=0.420), while L5dp responses were significantly smaller than those at L1/2 (p=0.0462; Figure 2G), as would be expected given the complete lack of LR dendrites in L5. In contrast, RS EPSP amplitudes at layers 1a, cell body, L5sup, and L5dp did not differ from responses at L1/2 (p>0.05 for all; paired t-test), indicating that there is no significant lamination of functional responses to thalamic inputs in the RS population, likely due to their relative lack of response to these inputs in general (Figure 2G). These relationships persisted when inhibition blocker picrotoxin (50 µM) or sodium channel blocker TTX (1 µM)+4 AP (100 µM) were added to the bath (data not shown; see Materials and methods). Thus, the connectivity between anterior thalamic TC arbors and L3 principal neurons in RSG highlights a circuit that depends on the precise sublaminar colocalization of both local principal cells’ dendrites and incoming thalamocortical axons. This demonstrates that LR neurons are uniquely anatomically adapted to respond to incoming inputs from the anterior thalamus, where ~ 60% of neurons are head direction cells (Taube, 1995; Taube and Bassett, 2003), due to a precise overlap of their apical dendrites with the dense thalamic axons and terminals in the uppermost sublamina of L1 (L1a).
Claustrum and anterior cingulate selectively control retrosplenial RS cells while avoiding LR cells
Two other main sources of input to RSG are the claustrum (CLA; Van Groen and Wyss, 2003; van Groen and Wyss, 1990) and anterior cingulate cortex (ACC; Shibata and Naito, 2008; Van Groen and Wyss, 2003; van Groen and Wyss, 1990). Both claustrocortical (ClaC) and ACC corticocortical (CC) projections have been anatomically shown to exhibit distinct laminar organization within its target cortical layers (Van Groen and Wyss, 2003; Vogt and Miller, 1983; Wang et al., 2017), but the precise functional targets of claustral and cingulate inputs to superficial RSG have, to our knowledge, never been examined. To address this, we next used CRACM, this time to examine CLA inputs to RSG LR and RS neurons and then compared these projections to RSG thalamic inputs. ChR2 injections into CLA (Figure 3A) resulted in expression of CLA axons and terminals in layers 1 c, 2, and 5 (Figure 3B). This expression pattern is distinct from the thalamic expression localized to layers 1a and 3 (Figure 2B–D). Dendritic lamination plots of LR and RS dendrites and projection density plots of CLA axons and terminal arbors show that RS, but not LR, apical dendrites anatomically overlap with CLA arbors in superficial layers 1 c and two and upper layer 5 (Figure 3D). Using the same targeted optogenetic stimulation approach as previously described (see Materials and methods), we found that RS neurons have significantly larger responses to CLA input compared to LR neurons when stimulated at layer 1/2 (p=0.0000035, Wilcoxon rank sum test; Figure 3E), the area of the strongest superficial projection density of CLA arbors. Similarly, layer 5 RS neurons exhibited significantly larger EPSP amplitudes in response to CLA inputs compared to LR cells (Figure 2—figure supplement 1). This again indicates that while LR neurons may be optimally positioned to receive thalamic input, RS neurons in both RSG L3 and L5 are instead the primary recipients of CLA inputs.
To further examine the lamination of ClaC inputs to LR and RS neurons and their resulting responses, we again conducted the high-resolution CRACM protocol using minimum LED power (Figure 3F; see Materials and methods). We found that RS neurons have the largest EPSP amplitude when the LED pulse is stimulating both layer 1c and 5sup and weakest EPSP amplitude at layer 1a, corresponding with, respectively, the strongest and the weakest projection density of CLA afferents at these laminar locations. In contrast, the complementary LR cell had no response to LED stimulation at any of the laminar locations (Figure 3F). Population analysis showed that within the RS cell population, RS cells have significantly larger responses at L1/2 (p=0.00725) and L5sup (p=0.0301) compared to those at L1a (paired t-test; Figure 3G). This lamination of functional responses in the RS cell precisely aligns with the lamination of ClaC axons and terminal arbors within these layers. In contrast, LR EPSP amplitudes across all layers were much lower in magnitude and did not differ from responses at L1a (p>0.05 for all; paired t-test; Figure 3G), indicating no significant lamination of functional responses to ClaC inputs in the LR population. These relationships persisted when inhibition blocker picrotoxin (50 µM) or sodium channel blocker TTX (1 µM)+4 AP (100 µM) were added to the bath (data not shown; see Materials and methods).
Using the same CRACM approach to examine ACC inputs to RSG (Figure 4A), we found that corticocortical inputs from the anterior cingulate target L1b/c and, to a lesser extent, L5 (Figure 4B), partially resembling the laminar pattern seen with CLA arbors (Figure 3B) and overlapping precisely with RS, but not LR, dendrites (Figure 3C&D). Indeed, when stimulating ACC arbors at L1/2, RS cells had significantly larger EPSP amplitudes compared to LR cells (p=0.0019). When examining the lamination of these inputs with higher spatial resolution, we found that RS cells have the largest EPSP response to ACC inputs at L1/2, while the same minimum LED stimulation power elicits no response in layer 1a (Figure 4F). Population analysis of RS cells revealed significantly larger responses to optogenetic stimulation of ACC inputs at the L1/2 (p=0.000037) and cell body (p=0.000525) compared to L1a stimulation (Figure 4G), reflecting the precise overlap of ACC arbors and RS dendrites. In contrast, the much smaller LR EPSP amplitudes across all layers did not differ from responses at L1a (p>0.05 for all; paired t-test; Figure 4G). Again, the same results were seen when inhibition blocker picrotoxin (50 µM) or sodium channel blocker TTX (1 µM)+4 AP (100 µM) were added to the bath (data not shown; see Materials and methods). Thus, both claustrocortical and corticocortical inputs to RSG target RS apical dendrites at the lower divisions of L1 (L1c and L1b/c, respectively), whereas thalamocortical inputs preferentially target the apical dendrites of LR cells in upper L1 (L1a). This again indicates that parallel circuits in RSG process TC versus ClaC and CC information, and this parallel processing is enabled by precise sublaminar organization of afferent axons and layer three principal cell apical dendrites.
Dorsal subiculum selectively controls LR, but not RS, cells
Previous work has shown that dorsal subiculum (DS), which serves to transmit allocentric spatial information such as axis cell and boundary vector signals (Lever et al., 2009; Derdikman, 2009; Olson et al., 2017; Simonnet and Brecht, 2019; Bicanski and Burgess, 2020), also precisely targets layer 2/3 principal neurons in RSG (Nitzan et al., 2020; Yamawaki et al., 2019a). These inputs overlap with LR cell bodies and basal dendrites and have been shown to evoke larger excitatory postsynaptic currents in superficial compared to layer 5 pyramidal neurons (Nitzan et al., 2020; Yamawaki et al., 2019a). These results suggest that DS projections to RSC may be targeting LR cells, potentially resulting in converging spatial information from both ADAV and DS onto LR cells. However, it remains unknown whether LR or layer 3 RS cells are the predominant target of these DS projections. To investigate this, we repeated our CRACM protocols and examined the functional responses of LR and layer 3 RS cells to DS inputs (Figure 5A&B). We found that DS afferents to layer 3 of RSG precisely overlap with LR basal and, to a lesser degree, proximal apical dendrites (Figure 5C) but not RS basal dendrites (Figure 5D). Correspondingly, LR neurons responded with significantly larger EPSPs compared to RS neurons when DS inputs were stimulated at the cell body in layer 3 (p=0.004; Figure 5E). Our high-resolution protocol and population analysis also confirmed that LR cells had significantly larger responses at layers 1/2 and 3 compared to L1a (L1/2 p=0.0132, L3 p=0.0135; Figure 5F&G), suggesting the possibility of integration of synchronous anterior thalamic and dorsal subicular inputs by LR neurons. In contrast, RS cells showed no significant laminal differences in response to DS input (p>0.05 for all; Figure 5G). Thus, inputs from both the dorsal subiculum and anterior thalamus target LR cells, reflecting their precise overlap with LR dendrites, while the claustrum and anterior cingulate cortex instead target RS cells.
Precise anatomical overlap of LR versus RS dendrites with distinct incoming axons facilitates parallel circuits in superficial RSG
We next extended our physiological examination of superficial RSG circuits to include correlational analyses and dual injection experiments (Figure 6A). We found that LR dendrites are strongly and significantly positively correlated with both ADAV (r = 0.88, p<0.001) and DS (r = 0.69, p=1.41e-254) axons. In contrast, RS dendrites are strongly correlated with CLA (r = 0.29, p=4.82e-36) and ACC (r = 0.43, p=5.73e-82) axons. Notably, LR dendrites are significantly negatively correlated with both CLA (r = −0.49, p=1.16e-110) and ACC (r = −0.15, p=1.00e-10) axons, and RS dendrites are negatively correlated with ADAV (r = −0.55, p=2.47e-145) and DS (r = −0.44, p=7.74e-85) axons (Figure 6B&C). These anatomical results mirror our functional findings that LR neurons are preferentially driven by both ADAV and DS inputs, while RS neurons are targeted by CLA and ACC inputs.
We also compared laminar expression of afferents from the four input regions and found that ADAV and DS inputs to RSG, which both at least in part target layer 3, were significantly positively correlated with one another (r = 0.55, p=2.93e-143). Similarly, CLA and ACC inputs, which both target the lower divisions of L1 as well as L5, were also significantly positively correlated (r = 0.23, p=2.3e-23). However, ADAV and DS inputs were significantly negatively correlated with both CLA and ACC inputs, confirming that these two streams target separate sublayers within RSG (Figure 6B,C). Dendrite-dendrite correlations of LR and RS cell populations were also strongly anti-correlated (r = −0.65, p=3.00e-215; Figure 6C), highlighting the distinct anatomical organization of these two neighboring principal cell types. Taken together, these axon-axon, dendrite-dendrite, and axon-dendrite correlations indicate that precise, fine-grained laminar organization of LR and RS dendrites and afferent axons creates a parallel circuit in which LR neurons are selectively optimized to integrate incoming spatial information (Figure 6D).
Supplemental examination of contralateral RSG (cRSG) inputs to RSG showed lamination and resulting cell-type-specific targeting resembling that of projections from CLA and ACC (Figure 6—figure supplement 1). Specifically, cRSG inputs targeted L1c, L2, and L5, resulting in strong activation of RS cells, with almost no responses in LR cells (Figure 6—figure supplement 1). This again suggests that inputs from CLA, ACC and cRSG follow a structured pattern by which these inputs preferentially control RS cells, while ADAV and DS projections target LR cells. Importantly, laminar dichotomies similar to those seen here in RSG also exist between TC and ClaC/corticocortical projections to other cortical regions such as the medial prefrontal cortex (Cruikshank et al., 2012), as also revealed by our examination of anatomical datasets from the Allen Brain Institute (Figure 6—figure supplement 2). Thus, our findings may also highlight a more universal framework by which other cortical regions integrate and process thalamic versus claustral and cortical inputs.
Anterior thalamic input to LR neurons is uniquely depressing
To investigate the functional implications of these parallel circuits, we next examined the short-term dynamics of each input onto the RSG principal cells. LR short-term dynamics were examined in response to ADAV and DS inputs, while RS short-term dynamics were examined in response to CLA and ACC inputs (see Materials and methods). In contrast to previously documented thalamocortical matrix inputs to both superficial pyramidal cells and interneurons (Anastasiades et al., 2021; Cruikshank et al., 2012), 10 Hz anterior thalamic inputs from ADAV to RSG LR neurons were uniquely and significantly depressing (p=1.02e-5, Wilcoxon rank sum test; Figure 7A&B). In contrast, 10 Hz DS inputs to LR neurons were not depressing and exhibited weak facilitation (Figure 7A&B). Both CLA and ACC inputs to RS cells were also weakly facilitating (Figure 7C&D). Thus, the synaptic depression of anterior thalamic inputs to LR neurons is unique among the inputs examined here and also distinct from anterior thalamic inputs to superficial principal neurons in PFC and ACC (Cruikshank et al., 2012). As expected, we found that ADAV inputs to LR cells do not show synaptic depression at the much slower stimulation frequency of 0.1 Hz (Figure 7—figure supplement 1 panel C), but 40 Hz inputs resulted in stronger depression than 10 Hz (p=0.0022, Wilcoxon rank sum test; Figure 7—figure supplement 1 panel B).
Short-term depression of anterior thalamic inputs enables encoding of angular head speed by LR cells
We next used computational modeling to elucidate the functional role of synaptic depression in this circuit, especially with respect to the processing of head direction input (Taube, 1995; Taube and Bassett, 2003). Up to 60% of cells in the anterior thalamus are classical head direction (HD) cells (Taube and Bassett, 2003). Each HD cell has a unique preferred direction at which its firing rate is highest when the head is facing that direction. The preferred directions of all cells in the HD ensemble span the full range of compass directions (Taube, 1998). Our modeling setup consisted of a postsynaptic RSG LR cell receiving input from an ensemble of 7500 presynaptic HD cells via depressing synapses (for full details, see Materials and methods). In the initial simulations shown in Figure 8, the morphologically realistic model LR cell (Brennan et al., 2020) received uniform inputs (all HDs being equally likely) via depressing synapses with Tsodyks-Markram short-term dynamics (Tsodyks et al., 1998). The synaptic parameters were fit to match the experimentally observed response of LR cells to optogenetic stimulation of thalamic afferents reported above (Figure 7).
We found that the firing rate of the postsynaptic LR cell receiving HD input via depressing synapses was strongly correlated with angular head speed, giving rise to symmetrical angular head velocity tuning (Figure 8B&C). Identical HD inputs transmitted via non-depressing synapses resulted in LR firing rates that were uncorrelated with angular head speed (Figure 8B,E–F). It is important to note that in our model, firing rates of input HD cells were not explicitly modulated by angular head velocity; therefore, the observed speed tuning resulted strictly from the depressing synaptic dynamics. For the parameter set whose results are depicted, LR cell firing rate was optimally correlated with head speed 24 ms in the past, confirmed by both cross-correlation and mutual information analyses (Figure 8E&F). We obtained the same results when modeling an ensemble of only 2500 presynaptic HD cells, indicating that our results are robust to the size of the presynaptic population (Figure 8—figure supplement 1). Taken together, these results suggest that synaptic depression of HD ensemble inputs introduces an angular head speed signal into the LR population, producing neurons that are more responsive during faster head turns, and potentially supporting the spatial orientation encoding functions attributed to the RSG (Epstein, 2008; Ino et al., 2007; Milczarek et al., 2018; Miller et al., 2019).
Anticipatory firing of thalamic head direction cells improves postsynaptic speed encoding by retrosplenial LR neurons
Anterior thalamic HD cells display anticipatory firing (Blair et al., 1997), a phenomenon where an HD cell becomes most active at a fixed time interval before the animal is facing that cell’s preferred direction (Figure 9A). This temporal offset associated with each cell is called its anticipatory time interval (ATI). In anterior thalamus, the mean ATI has been reported as 25–50 ms (Taube, 2010). In our initial simulations, shown in Figure 8, we drew the ATI of each HD cell randomly from distributions matching these known in vivo ranges. To more systematically understand if and how the ATI of HD cells influences angular head speed coding by LR neurons, we next performed a series of simulations using various fixed ATI values for the entire HD population (Figure 9A) and analyzed the resulting postsynaptic LR response.
Increasing ATI from 0 ms (no anticipation by the HD cell) to higher values (HD cell firing prior to facing the preferred HD) improved the lag between LR firing rate and angular head speed, with larger ATIs resulting in shorter latency between the head movement and LR coding of that angular head speed (Figure 9C). Remarkably, larger ATIs also improved the angular head speed tuning of LR neurons independent of lag, as quantified by the maximum value (over the full range of lags) of cross-correlation or cross-mutual information (Figure 9B). Thus, anticipatory firing of HD cells may constitute a powerful coding principle in the thalamo-retrosplenial circuit, helping LR cells to not only encode the current head speed with minimal lag, but also better encode the head speed independent of the lag.
In order to better understand this surprising effect, as well as how speed tuning generally arises, we analytically studied a simplified mean-field model of this thalamo-retrosplenial circuit (Figure 9—figure supplement 1). Our analysis allowed us to mathematically derive that LR activity should encode head speed with depressing synapses from simplified HD cells (see Appendix). Moreover, the analysis showed that anticipatory firing compensates for the lag introduced by integration time, leading to a theoretical parameter regime for which postsynaptic speed coding can be essentially perfect (where LR firing rate reflects head speed exactly with minimal latency). The observed improvement in quality of speed coding with increases in ATI corresponds to moving closer to this regime. Finally, our analysis showed that postsynaptic activity, at least at lower rotational speeds, is proportional to the square of head speed, thereby explaining the concave-up, parabolic shape of the speed-firing rate curve (Figure 8C and Figure 8—figure supplement 2). Thus, both the simulations and the analytical theory show that synaptic depression of HD cell input onto LR cells should result in angular speed coding in the LR cell firing rate, and that anticipatory firing improves both the precision and quality of this coding.
Non-uniform HD inputs can allow LR cells to encode both head direction and speed, with a tradeoff
In the simulations above, we utilized a uniform distribution of HD inputs to each LR cell, such that each preferred direction was equally represented. However, in practice, it is likely that there will be heterogeneity in the preferred HDs of the cells providing inputs to any given retrosplenial LR cell. To study the extent to which speed coding persists under these conditions, we simulated a population of 200 LR cells, each of which received inputs from a randomly selected distribution of preferred directions for the presynaptic HD cells (Figure 10Ai,,Bi,Ci). We then analyzed the information content of each LR cell’s firing rate, assessing the correlation and mutual information with both head direction and speed. We found that single cell activity in the LR population could encode both head direction and head speed (Figure 10Aii,Bii,Cii), but there was a clear tradeoff between the two types of information.
This tradeoff between directional versus speed encoding was explained by the level of heterogeneity in the presynaptic preferred direction distribution (Figure 10D–F), making intuitive sense. If an LR cell receives inputs that represent each HD equally, then it is not biased towards any one HD and encodes direction poorly. However, it encodes head speed very strongly in this case because it is able to sample the speed across all directions—there are no directional blind spots where the LR cell cannot detect the speed of head rotation from its presynaptic inputs. On the other hand, consider an LR cell that, by chance, gets inputs from a higher proportion of HD cells encoding a direction around 180 degrees. Now, the output of this LR cell is better correlated with direction, acting as a broadly tuned HD cell itself and preferring 180 degrees, but also firing at most other directions. However, this LR cell’s encoding of head speed gets worse because it under-samples many of the head directions faced by an animal. Thus, non-uniform HD inputs via depressing thalamocortical synapses represent a simple mechanism by which LR cells can conjunctively encode both head direction and head speed. Overall, our results predict in vivo LR neuronal spike trains that are likely to correlate with both the head direction and rotational speed of an animal, with a tradeoff between direction versus speed coding seen across individual LR cells (Figure 10E).
Discussion
We have identified a circuit that supports parallel processing of information as it first reaches the superficial granular retrosplenial cortex. Layer 3 of RSG contains two distinct pyramidal neuronal subtypes: small, excitable low rheobase (LR) neurons and regular spiking (RS) cells (Brennan et al., 2020). Here, we show that LR cells are driven by inputs from the anterior thalamus and dorsal subiculum but are essentially unresponsive to inputs from the claustrum, anterior cingulate, and contralateral retrosplenial cortex. (Figures 1–6, Figure 6—figure supplement 1). RS cells show precisely the opposite relationship: they are strongly driven by claustral, anterior cingulate, and contralateral retrosplenial inputs, but very weakly activated by anterior thalamic or dorsal subicular inputs (Figures 1–6, Figure 2—figure supplement 1, Figure 6—figure supplement 1). This dichotomy can be explained in large part by the precise overlap of LR versus RS apical dendrites with distinct afferent axons, including clear sublamination of both axons and dendrites within layer 1 (Figure 6). Anterior thalamic and dorsal subicular cells show robust directional and spatial modulation (Barry et al., 2006; Olson et al., 2017; Stewart et al., 2014; Taube, 1998; Taube and Bassett, 2003). Thus, the parallel circuit described here allows LR cells to selectively process this spatially relevant information and makes them ideally suited to support the spatial orientation computations carried out by the RSG.
LR cells are morphologically and computationally unique
The small, excitable LR cells are unique among pyramidal cells in RSG, and there is no evidence, to our knowledge, for such cells in any other cortical regions (Brennan et al., 2020; Holmgren et al., 2003; Jiang et al., 2015; Kurotani et al., 2013; Wyss et al., 1990). Indeed, the ‘granular’ in granular retrosplenial cortex refers to the appearance of these dense small cell bodies in superficial RSG (Sripanidkulchai and Wyss, 1987; Vogt and Peters, 1981). Thus, LR cells are the defining morphological feature of RSG.
The distinct electrophysiological properties of LR neurons, specifically, play an integral role in the RSG’s unique capacity to process anterior thalamic head direction inputs. We have previously reported that the low rheobase and lack of spike frequency adaptation of these neurons are their defining computational features (Brennan et al., 2020). For this reason, we refer to them as low-rheobase (LR) neurons. They have also been called late-spiking (LS) neurons (Kurotani et al., 2013; Yousuf et al., 2020), but we argue this name is not optimal for several reasons. First, many other neurons in superficial RSG are also late spiking (Brennan et al., 2020), with L2/3 fast-spiking (FS) interneurons and also some RS neurons exhibiting a substantial delay to first spike, likely due to increased Kv1.1 or 1.2 channel expression in this region compared to other cortical areas (Kurotani et al., 2013). Thus, late-spiking is neither a distinctive feature of LR neurons nor excitatory neurons in this region. Second, the name ‘late-spiking (LS) neuron’ is already the accepted and widely used nomenclature referring to a distinct inhibitory neuronal subtype, the neurogliaform cells that are located throughout the cortical layers in many regions of the brain, including within layer 1 of RSG (Cruikshank et al., 2012; Hestrin and Armstrong, 1996; Kawaguchi and Kubota, 1997; Overstreet-Wadiche and McBain, 2015). Instead, the defining computational feature of the small pyramidal cells in L2/3 of the RSG is their low rheobase, and this nomenclature serves to provide an unambiguous name for these hyperexcitable cells (Brennan et al., 2020).
LR cells respond most strongly to inputs from spatially relevant regions
Previous work has reported TC-evoked responses in layer 5 RS neurons at their apical dendrites (Yamawaki et al., 2019b), with some variation in the magnitude of responses across cells. Our results reproduce these observations (Figure 2—figure supplement 1). Layer 5 is not a homogenous layer, often divided into sublayers 5A and 5B (Sempere-Ferràndez et al., 2019; Sempere-Ferràndez et al., 2018; Sigwald et al., 2020; Sripanidkulchai and Wyss, 1987; Yamawaki et al., 2016a) and contains a variety of pyramidal cells, including the thin-tufted and thick-tufted neurons who have been shown to exhibit different responses to external inputs (Sempere-Ferràndez et al., 2018; Wyss et al., 1990). Further work is needed to determine whether any particular subtype of layer five pyramidal neurons may exhibit substantial TC-evoked responses, but our results, reporting significantly larger responses to TC input by LR neurons, suggest LR cells are the predominant spatial information-encoding subtype within RSG.
The dorsal subiculum, which serves to transmit allocentric spatial information such as axis and boundary vector signals (Lever et al., 2009; Derdikman, 2009; Olson et al., 2017; Simonnet and Brecht, 2019; Bicanski and Burgess, 2020), also precisely targets LR neurons via projections predominantly to layer 3 of RSG (Figure 5; Nitzan et al., 2020; Yamawaki et al., 2019a). These inputs overlap with LR cell bodies, basal dendrites, and apical dendrites in L1b/c and evoke larger excitatory postsynaptic currents in superficial cells compared to layer 5 pyramidal neurons (Figure 5; Nitzan et al., 2020; Yamawaki et al., 2019a). Thus, the distinct morphology and intrinsic properties of LR neurons make them ideally suited to integrate head direction input from the thalamus and spatial inputs from the hippocampal formation via the dorsal subiculum. Indeed, this integration of various types of directional and distance information is often stated as the key computational function of the retrosplenial cortex as a whole (van Wijngaarden et al., 2020; Burgess et al., 2001; Byrne et al., 2007; Epstein, 2008; Ichinohe et al., 2003; Maguire, 2001), further highlighting the defining role that LR cells are likely to serve in the RSG.
Thalamic inputs to LR cells show short-term depression: implications for angular head velocity and head direction coding in the RSG
We have shown that thalamocortical synapses onto LR cells are depressing (Figure 7). Our modeling results demonstrate that such depressing synapses allow LR neurons to compute angular head speed from head direction input, leading to robust encoding of angular speed in the postsynaptic LR firing rate (Figure 8). This result is in line with previous studies (Abbott et al., 1997; de la Rocha and Parga, 2008; Puccini et al., 2007) that use synaptic depression to implement neural circuits that perform rate of change computations. We also find that LR encoding of angular head speed is improved by the fact that input HD cells display anticipatory firing (Figure 9; Blair et al., 1997). This improvement not only allows LR cells to encode current head speed with a shorter time lag (Figure 9C), but also strengthens the correlation with head speed regardless of lag (Figure 9B). This proposed mechanism by which depressing synapses compute head speed is fully generalizable, and thus may also be of relevance to other brain circuits employing HD-like population codes (Maunsell and Van Essen, 1983). Furthermore, we have shown that this thalamocortical synapse can support conjunctive encoding of both head direction and angular head speed by individual LR cells (Figure 10).
Our results make specific predictions about the coding properties of LR cells that can be tested with in vivo recordings. First, based on their intrinsic properties (Figure 1—figure supplement 1; Brennan et al., 2020), LR cells in superficial RSG are likely to show high firing rates and spike widths intermediate to the standard ‘broad/narrow’ criteria currently used to distinguish putative RS and FS cells in most cortical recordings. Thus, future recordings will need a large sampling of cells in L2/3 of RSG to provide enough data to identify the predicted three distinct clusters of spike shapes (RS/LR/FS). Second, we expect LR cell activity to be correlated with angular head speed with a relatively small lag (Figures 8 and 9). Third, we predict that the experimentally observed head velocity tuning curves will be concave up and exhibit approximately quadratic scaling at low rotational speeds, as our analytical work suggests (Figure 8; Appendix). Fourth, we anticipate a significant proportion of LR cells will show conjunctive encoding, where both angular head speed and direction are encoded. A key aspect of this conjunctive encoding, as predicted by our simulations (Figure 10), is that individual cells will show a tradeoff between directional and speed encoding. As a population, this is likely to allow LR cells to comprehensively encode conjunctive speed and directional information and provide this information to downstream targets, helping with spatial orientation computations within RSG. This predicted conjunctive encoding will also necessitate that the test of predictions 2 and 3 be carried out in a head-direction delimited manner.
Our modeling work so far has focused on understanding the computational implications of thalamic inputs to LR cells. Future work, beyond the scope of the present study, will be needed to understand how these thalamic inputs interact precisely with dorsal subiculum inputs. Given the diversity of spatial and directional encoding seen in the dorsal subiculum (Lever et al., 2009; Derdikman, 2009; Olson et al., 2017; Simonnet and Brecht, 2019; Bicanski and Burgess, 2020), computational models will allow for the rigorous undertanding of how each possible type of subicular input (e.g. axis versus boundary vector cells) to LR cells interacts with thalamic head direction inputs. Very recent work on understanding the tuning properties of subicular cells in a projection-specific manner will help to narrow the functional range of possible dorsal subicular inputs to LR cells (Kitanishi et al., 2021). We anticipate that one key impact of dorsal subicular inputs to LR cells will be to impose a spatial filter on the angular speed and directional predictions stated above. In particular, the precise organization of layer 1 dendrites and axons may make LR cells most responsive to the near-synchronous arrival of inputs from both thalamic head direction cells (at distal apical dendrites in L1a) and dorsal subicular boundary vector or axis cells (at more proximal apical dendrite locations in lower L1). Similar cooperative dendritic activation patterns are also a hallmark of integration of entorhinal and CA3 inputs by CA1 pyramidal cells, where they help to overcome the effects of strong dendritic inhibition and shape the firing of place cells (Ahmed and Mehta, 2009; Golding and Spruston, 1998; Kamondi et al., 1998; Takahashi et al., 2012).
Downstream targets of LR cells and potential role of RS neurons in retrosplenial spatial computations
LR neurons do not synapse onto neighboring excitatory neurons within layers 2/3 (Brennan et al., 2020), and no (0/9) pairs tested in this study between LR and layer 5 RS cells exhibited connectivity (data not shown). Thus, LR cells likely do not provide any substantial input to other LR cells or any RS cells located in layers 2–5 of ipsilateral RSG. Instead, LR neurons send their axons into the corpus callosum and may likely target homotopic regions of contralateral RSG (Brennan et al., 2020; Kurotani et al., 2013). Indeed, retrograde injections into one RSG hemisphere have been shown to predominantly label small contralateral L2/3 cells (Sripanidkulchai and Wyss, 1987; Van Groen and Wyss, 2003), most likely corresponding to the LR cells that constitute the majority of cells in L2/3 (Brennan et al., 2020). Our supplementary results (Figure 6—figure supplement 1) show that inputs from contralateral RSG preferentially and strongly target RS neurons, with almost no recruitment of LR neurons. This is consistent with the observed distribution of afferent axons from contralateral RSG: these axons and terminals avoid L1a and L3 where all LR dendrites are located (Figure 6—figure supplement 1). Thus, we propose the following hypotheses for LR downstream targets and possible function: LR neurons in one hemisphere receive unilateral inputs from directional cells in the anterior thalamus (Figures 1, 2 and 6). These inputs are utilized to encode information about both the speed and direction of head rotations and are subsequently transmitted to RS cells in contralateral, but not ipsilateral, RSG (Figure 6—figure supplement 1). Thus, LR cells may serve to perform a hemispheric switch of rotational speed and direction information. Given the known bidirectional connectivity between secondary motor cortex and RSG (Yamawaki et al., 2016a), RS cells may be able to utilize inputs from contralateral LR cells to perform comparisons of the actual head speed and direction signal with motor efference signals, helping to sharpen the code for spatial orientation in the RSG. However, extensive future investigations of this contralateral circuit are needed to explore these hypotheses.
Potential role of claustral and anterior cingulate inputs
In contrast to thalamic and DS inputs to RSG, projections from both the claustrum and anterior cingulate cortex avoid laminar overlap with LR apical dendrites, instead projecting to the lower divisions of L1/2 and L5 to target RS apical and basal dendrites (Figures 3 and 4, and 6). Both the CLA and ACC are heavily implicated in a range of non-spatial behaviors often described as being ‘cognitive’ or ‘mnemonic’ in function (Botvinick, 2007; Brown and Braver, 2005; Carter et al., 1998; Devinsky et al., 1995; Kim et al., 2016; Mohanty et al., 2007; Smith et al., 2012; Smythies et al., 2012; Stevens et al., 2011), but there is also evidence that these regions exhibit some degree of spatial information processing (Guterstam et al., 2015; Jankowski and O'Mara, 2015; Sutherland et al., 1988; Whishaw et al., 2001). Indeed, 7–15% of cells in the claustrum were found to encode place or boundary-related information (Jankowski and O'Mara, 2015), though this percentage is relatively low compared to the 60% of anterior thalamic (Taube, 1998) and 61% of dorsal subicular (Kitanishi et al., 2021) neurons that are known to clearly encode elemental direction, place, and other spatial features. Regardless, retrosplenial RS neurons may be receiving direct spatial inputs from the CLA and ACC in addition to the cognitive and mnemonic information relayed by these two regions. We also speculate that retrosplenial RS neurons may receive directional and spatial information that has already been heavily integrated by LR neurons in the contralateral hemisphere (Figure 6—figure supplement 1), as discussed above. This convergence of disparate signals onto RS cells could then serve to support functions that require both spatial and non-spatial integration, such as contextual fear conditioning, a key function attributed to the RSC (de Sousa et al., 2019; Keene and Bucci, 2008a; Keene and Bucci, 2008b; Kwapis et al., 2015; Yamawaki et al., 2019a).
Indeed, a recent study reported that selective ablation of putative retrosplenial RS cells located near the boundary and within upper layer 5 resulted in contextual fear amnesia (Sigwald et al., 2020), indicating that these cells are involved in the fear-related functions of the retrosplenial cortex (Lukoyanov and Lukoyanova, 2006; Yamawaki et al., 2019a; Yamawaki et al., 2019b). The laminar location of these cells corresponds well to the position of the RS cells recorded in our study that we show to preferentially receive CLA and ACC inputs. Both the CLA (Cho et al., 2017; Ipser et al., 2013; Vakolyuk et al., 1980; Vetere et al., 2017; Zingg et al., 2018) and ACC Frankland et al., 2004; Han et al., 2003; Steenland et al., 2012 have been heavily implicated in fear-associated behaviors, suggesting again that CLA and ACC inputs to L3 RS cells help to support the role of these neurons in contextual fear conditioning. Future work that selectively inactivates CLA versus ACC axon terminals in RSC will be necessary to causally confirm the role of each of these synapses in fear conditioning.
The role of feedforward inhibition
Another important consideration within the circuits examined here is the role of feedforward inhibition. Matrix TC inputs evoke feedforward inhibition to regulate cortical signaling through both disinhibition (Anastasiades et al., 2021; Delevich et al., 2015) and the disynaptic inhibition of pyramidal neurons (Cruikshank et al., 2012; Yamawaki et al., 2019a; Yamawaki et al., 2019b). Similarly, ClaC inputs invoke feedforward inhibition in several cortical regions, particularly through recruitment of neuropeptide Y (NPY) and, to a lesser degree, FS cells (Jackson et al., 2018). Long-range inhibitory signals from CA1 of the hippocampus also target apical dendrites in retrosplenial layer 1a, converging with excitatory anterior thalamic inputs, to precisely regulate these TC inputs and establish a hippocampo-thalamo-retrosplenial network (Yamawaki et al., 2019b). This fine-grained laminar overlap strongly suggests that these inhibitory signals may regulate responses to matrix TC input more strongly for LR than layer 5 RS cells and should be investigated further, particularly in the context of integration of spatial information processing. Future work will examine the effect of feedforward inhibition on the subpopulations of retrosplenial principal neurons in response to these converging inputs in order to establish a thorough understanding of the circuitry facilitating retrosplenial information processing, particularly with regard to head direction signals.
Conclusions
In summary, our results highlight a superficial retrosplenial circuit enabled by the precise sublaminar organization of distinct principal neuronal subtypes and axonal afferents. LR neurons receive directional and spatial inputs from the anterior thalamus and dorsal subiculum. The synaptic dynamics of thalamic inputs to LR cells can give rise to rate coding of angular head speed in LR cells. In contrast, neighboring RS neurons in RSG respond very weakly to these directional inputs, as their apical and basal dendrites distinctly avoid laminar overlap with anterior thalamic and subicular afferents. Instead, RS neurons respond to claustrocortical and anterior cingulate inputs. Determining how these two parallel streams of information are integrated in downstream neurons within both the ipsilateral and contralateral granular and dysgranular retrosplenial cortices is the next critical step toward a mechanistic understanding of retrosplenial computations.
Materials and methods
Animals
All procedures and use of animals were approved by the University of Michigan Institutional Animal Care and Use Committee. The following mouse lines were used for the CRACM experiments: Ai14 (Jackson Laboratories, 007914), Ai32 (Jackson Laboratories, 024109), C57BL6 wildtypes (Charles River, stock #027), Camk2αCre (Jackson Laboratories, 005359), Grp-KH288Cre (RRID:MMRRC 037585-UCD), PvalbCre (PV-IRES-Cre; Jackson Laboratories, 008069), and PvalbCre x Ai14 (crossed in house). A combined total of 70 mice of both sexes between the ages of P41–P184 were used in this study.
Experimental procedures
Surgical procedures
Request a detailed protocolSurgical anesthesia was induced via isoflurane inhalation at 4% and then maintained on a surgical anesthetic plane at 2–3% isoflurane. Upon induction, atropine was injected subcutaneously at 0.05 mg/kg. A Physitemp (Clifton, New Jersey, USA) controller monitored and maintained body temperature at 37°C. Ophthalmic ointment was placed on the eyes. The incision site was prepared using 1:40 Nolvasan followed with isopropyl alcohol before being subcutaneously injected with 1% lidocaine. The skull was then leveled, and bregma was identified. Using a digital stereotaxic coordinate system, the following injection target sites were identified: right anterodorsal thalamic nucleus (AP = −0.6 mm, ML = +0.78 mm, DV = −2.5, –3.25 mm), right claustrum (AP = +1.25 mm, ML = +2.6 mm, DV = −3.3 mm), right dorsal subiculum (AP = −3.0 mm, ML = +1.5 mm, DV = −1.5, –1.8 mm), right anterior cingulate cortex (AP = +0.3 mm, ML = +0.3 mm, DV = −0.7, –1.0 mm), and contralateral RSG (see Figure 6—figure supplement 1; AP = −2.25 mm, ML = −0.2 mm, DV = −0.7 mm). Burr holes were then drilled through the skull at the identified sites, and dura removed. Micropipettes were lowered under stereotaxic guidance into the target injection site containing the ChR2 viral construct (AAV2-EF1a-DIO-hChR2(h134R)-eYFP for GrpCre mice and AAV2-hsyn-ChR2(H134R)-eYFP or AAV2-hsyn-ChR2(H134R)-mCherry for all other lines, UNC Gene Therapy Vector Core). For dual injection surgeries, the ChR2 viral construct was injected into the functional target, and a fluorescent tag (AAV2-hsyn-mCherry or AAV2-hsyn-eYFP) was injected into the other region for comparison of various afferent axonal and terminal arborizations within the same animal.
Injections of 0.5 µL total virus volume at each depth were given via a picospritzer at 0.05–0.07 µL/min with a 5-min pause between injecting the more dorsal and ventral DV coordinates, when applicable. After a region was injected, there was a 10-min period before removing the micropipette from the brain. Enrofloxacin was administered at 8.0 mg/kg after injections. Burr holes were sealed with bone wax, and the incision was closed with VetBond with antibiotic ointment placed under skin edges. Isoflurane was tapered down prior to removal. After removal from isoflurane, carprofen was administered at 5 mg/kg. The mice were kept warm through an artificial heat source during the recovery period. Mice then recovered for 2–14 weeks post-injection before being used for slice experiments.
Slice preparation
Request a detailed protocolSlices were prepared as described previously (Brennan et al., 2020). Briefly, mice were deeply anesthetized using isoflurane before decapitation. Brains were removed and placed in a carbogen-saturated ice-cold high-sucrose slicing solution within 30 s of decapitation. Using a Leica 1200VT or Leica 1000S vibratome, 300 µm coronal slices were cut and placed in a high-magnesium artificial cerebrospinal fluid (ACSF) solution at 32°C. After resting in this solution for 20 min, the entire bath was moved to room temperature where the slices rested for the remainder of the experiment.
Whole-cell recordings
Request a detailed protocolDuring whole-cell recordings, slices were submerged in a recording chamber with a 2 mL/minute flow of body-temperature ACSF (126 mM NaCl, 1.25 mM NaH2PO4, 26 mM NaHCO3, 3 mM KCl, 10 mM dextrose, 1.20 mM CaCl2, and 1 mM MgSO4). They were visualized using an Olympus BX51WI microscope, Olympus 60x water immersion lens, and Andor Neo sCMOS camera (Oxford Instruments, Abingdon, Oxfordshire, UK). Patch pipettes had a diameter of 2–4 µm and resistances between 2 and 5 MΩ. The potassium gluconate internal solution used in these experiments contained 130 mM K-gluconate, 2 mM NaCl, 4 mM KCl, 10 mM HEPES, 0.2 mM EGTA, 0.3 mM GTP-Tris, 14 mM phosphocreatine-Tris, and 4 mM ATP-Mg and had a pH of 7.25 and osmolarity of 290 mOsm.
Current-clamp recordings were conducted using the Multiclamp 700B and Digidata 1550B (Molecular Devices). Patched neurons were adjusted for series resistances and held at resting potentials around −65 mV. Recordings were not corrected post hoc for liquid junction potential. Resting membrane potential, defined as recorded potential within 30 s of break-in, was recorded, and cells with resting potentials more depolarized than −50 mV were not included in this study. Intrinsic neuronal properties were calculated using a set of protocols detailed below and measured with either Clampfit or simple custom MATLAB routines (Figure 1–source code 1). Pharmacological agents were prepared prior to experiments as outlined by the manufacturer and added to the ACSF. Agents were applied for at least 10 min before conducting experiments, and verification of effect was visualized (no spiking for strong, suprathreshold current injections for TTX-4AP; no feedforward inhibition for picrotoxin).
Channelrhodopsin-assisted circuit mapping
Request a detailed protocolChannelrhodopsin-assisted circuit mapping (CRACM) experiments (Yamawaki et al., 2016b) were conducted under the same rig set-up as described above while using a 5500K white light-emitting diode (LED; Mightex; maximum power of 14.47 mW measured at the slice focal plane). Synaptic responses to optical stimulation of the ChR2-expressing axons were measured from postsynaptic retrosplenial neurons recorded under whole-cell current-clamp conditions.
To examine the effect of afferent axonal input to retrosplenial neurons, two optogenetic stimulation protocols were used. First, a targeted optogenetic stimulation approach was used to examine the distribution of afferent axons and terminal arbors across the cortical layers in which the postsynaptic retrosplenial neurons resided. The main areas of interest included layer 1a (L1a; defined as the top 1/3 of L1 adjacent to the midline), layer 1/2 boundary (L1/2; defined as the border between the sparse L1 and densely-packed L2), the location of the patched cell body (all patched cells were located in layers 3 or 5), upper layer 5 (L5sup: defined as 100 µm from the L3/5 boundary), and deep layer 5 (L5dp; defined as 200 µm from the L3/5 boundary). For LED-based stimulation over L1a, the objective was centered approximately 30 µm from the pia in L1. For LED-based stimulation over L1/2, the objective was centered over the border between L1 and L2. For LED-based stimulation over the cell body, the objective was centered on the patched cell body at its widest diameter. For L5sup, the objective was centered 100 µm into L5, and for L5dp, the objective was centered 200 µm into L5. All LED-based stimulations were aligned in the cortical column with the patched cell body. LED intensity was kept constant at each stimulation location and was set at the minimum intensity necessary to evoke the smallest visible response when stimulating at the minimum location (see below). For most cells, this minimal response was between 0.25 and 5 mV, while strongly activated cells had larger responses with the smallest possible LED intensity. The minimum location for ADAV experiments was L1/2, while the minimum location for CLA, ACC, and DS experiments was L1a, corresponding with the lowest projection density of inputs within the superficial layers of RSG. All layer test optogenetic stimulation protocols consisted of a 1 s 10 Hz train of 1 ms LED pulses. For short-term dynamics analyses, an additional optogenetic stimulation protocol was used that consisted of a 1 s 40 Hz train of 1 ms LED pulses.
Second, a current step protocol was used to examine the effect of thalamic input on the postsynaptic retrosplenial neurons’ spike trains. In this protocol, a 2 s current step was delivered to the postsynaptic cell that elicited a 10–30 Hz spike train with a 1 ms LED pulse delivered at 500 ms into the current step. The objective for this protocol was centered directly over the cell body.
Morphological investigations and reconstructions
Cell filling and visualization
Request a detailed protocolTo analyze the morphology of the cells, biocytin (5 mg/ml) was added to the recording solution immediately before recording. Biocytin was allowed to diffuse into the cell for no less than 20 min. Shortly before removing the patch pipette from the neuron, ten current pulses (1-3 nA at 1 Hz) were applied to aid the diffusion process (Jiang et al., 2015). After the fill process was complete, the patch pipette was retracted from the cell slowly to allow the membrane resealing. One to four cells were filled per slice, and then the slice was transferred from the recording chamber to 4% PFA for overnight fixation. The next day, slices were washed in PBS and incubated for 24–48 hr in streptavidin conjugated Alexa Fluor (488, 594, or 647) with 0.2% TritonX added to permeabilize the cells. For a subset of cells, fluorescent nissl stain, NeuroTrace, was added to the Alexa Fluor incubation step (at 1:200 dilution). After this incubation, slices were washed in PBS, mounted on slides, and cover slipped using Fluoromount-G.
Cell reconstruction
Request a detailed protocolZ-stacks of each filled cell were acquired using the Leica SP5 confocal microscope with a dry 40x lens or the Zeiss Axio Image M2 confocal microscope with 20x lens. Reconstructions from z-stacks were performed using NeuroLucida software in user-guided mode or NeuTube software (Feng et al., 2015). Additional reconstructions previously generated were also added to the analyses described below (Brennan et al., 2020).
Imaging of slice expression
Request a detailed protocolRepresentative images of axonal projection to RSG were obtained from mice injected with anterograde AAV2 into ADAV, CLA, ACC, DS, and cRSG. Slices were prepared as for electrophysiology recordings. After use, slices with ADAV, DS, and cRSG expression from AP range −1.7 to −2.06 mm were fixed overnight in 4% PFA and then stained with the fluorescent Nissl stain, NeuroTrace 435–455, to aid laminar demarcation. ACC and CLA slices underwent an additional signal amplification protocol. Those slices were first incubated overnight in PBS containing appropriate normal blocking serums with 0.2% TritonX. On the following day, anti-GFP (for eYFP) or anti-mCherry primary antibodies (Abcam ab13970 and ab167453 respectively, diluted 1:2000) were added and slices were incubated at 4°C on a shaker for up to 24 hr. Next, slices were washed in PBS and incubated on a shaker for 3 hr at room temperature in PBS containing 2% bovine serum albumin (BSA), 0.2% TritonX and biotin-conjugated secondary antibodies (Jackson ImmunoResearch AB_2313596 and AB_2340593, dilution 1:300) then again washed in PBS. Finally, slices were incubated on a shaker in PBS with 0.2% TritonX and streptavidin-conjugated Alexa Fluor (488 for eYFP and 594 for mCherry amplification, dilution 1:300) and NeuroTrace 435–455 (dilution 1:200) for another 3 hr at room temperature, and then washed in PBS. All slices were then mounted on slides using Fluoromount-G and allowed to dry overnight at room temperature. Z-stack images (7 µm at 0.5 µm z-steps) of the slices were obtained the following day using confocal microscope Zeiss Axio Image M2 with 20x dry objective.
Quantification and statistical analysis
Neuronal analysis and statistics
Request a detailed protocolMultiple intrinsic neuronal properties were calculated as previously reported (Brennan et al., 2020). Briefly, spike threshold, spike amplitude, spike width, spike frequency adaptation ratio, latency to first spike, rheobase, input resistance (Rin), input capacitance (Cin), and membrane time constant (τm) were measured. Spike amplitude, threshold, and width were calculated from the average of all spikes in a 600 ms current step that elicited a 5 Hz spike train. Amplitude was calculated as the voltage difference from threshold to the peak of the spike, threshold from the peak of the third derivative of membrane potential (Cruikshank et al., 2012), and width as the full width at half-maximum spike amplitude. Spike frequency adaptation ratio was calculated from the same 600 ms current step protocol using the first sweep that elicited a 10 Hz spike train using the equation ISIlast/ISIfirst. Latency to first spike and rheobase were measured using a 1 s current step protocol that began below threshold and increased by 1–5 pA steps until at least three sweeps post-threshold. Latency to first spike was calculated as the time from the onset of the current step to the peak of the first spike. Rheobase was measured as the minimum current needed to cause at least one spike in the 1 s current step. Lastly, input resistance (Rin), input capacitance (Cin), and membrane time constant (τm) were all measured from a protocol that delivered a series of small negative current steps resulting in ~4 mV deflection in membrane potential. Rin was calculated using Ohm’s law, mean voltage change divided by mean current. τm was measured by fitting a single exponential to the average of the initial 60 ms of the negative voltage deflection, ignoring the first 20 ms. Lastly, Cin was calculated using the formula τm = Rin x Cin. Statistical significance of the differences in intrinsic properties between the retrosplenial neuronal subtypes was calculated using the Wilcoxon rank sum test.
ChR2-assisted circuit mapping analysis and statistics
Request a detailed protocolTo quantify the effect of TC and ClaC input across the layers, the amplitude of the resulting postsynaptic responses from each laminar location of LED stimulation (L1a, L1/2, and cell body) were measured with Clampfit and Matlab. The average of the response to the first pulse in the 10 Hz train was calculated for each cell at each LED laminar location. Significant differences between cell types and/or laminar locations were calculated using the Wilcoxon rank sum test. From the current step protocol, raster plots and PSTHs were plotted for each cell in MATLAB. The percent increase in spike rate was calculated by subtracting the number of spikes in a 100 ms window after the LED pulse from the number of spikes in a 100 ms window before the LED pulse, then dividing this value by the number of spikes in the 100 ms window before the LED pulse. Significant differences between spike count pre- and post-stimulation for both LR and RS cells and LR and RS percent increase in spike rates were calculated using the Wilcoxon rank sum test.
Short-term dynamics were calculated from 10 Hz or 40 Hz LED pulses protocols (see Channelrhodopsin-assisted circuit mapping subsection of the Experimental Procedures section of Methods). To analyze 0.1 Hz short-term dynamics, the first response to the first pulse of each sweep from the 10 Hz protocols was used. Average response amplitudes of each pulse were measured for each cell. Amplitudes were then normalized using the equation stimulus n/stimulus 1, and normalized postsynaptic ratios were plotted. Cells with first pulse responses greater than 0.5 mV but not spiking were included. Significance of resulting short term synaptic depression of LR neuron responses to 10 Hz and 40 Hz thalamic stimulation were calculated using the Wilcoxon rank sum test.
Morphological analysis and statistics
Request a detailed protocolAll z-stack images of axonal projections and their corresponding NeuroTrace channel were first collapsed in ImageJ (using Z project option with max intensity). A section of image spanning layers 1–6 from each image was chosen to serve as a representative sample of expression pattern in that image. In each sample, boundaries of retrosplenial layers 1, 2, 3, and 5 were measured as distance from pia based on NeuroTrace expression, and the information was saved. L1 was distinguished by a low cell density, L2 by a thin, densely packed band of small pyramidal cell bodies, L3 by a lower density (compared to L2) of small pyramidal neurons, L5 as having larger cell bodies relative to L3, and L6 as having smaller cell bodies. Images of fills of the reconstructed cells underwent a similar process (with the exception of LR neurons, which only had layers 1, 2, and 3 boundaries measured, and the average layer width for L5 was applied to them by default). After all images (both cell reconstructions and sample projections into RSG) were measured and processed, they were turned into binary black and white thresholded versions using ImageJ. For the images of reconstructed cells, cell bodies were removed from the black and white images to allow the analysis of dendrites only.
Dendrites and axons were analyzed as a function of cortical depth (Yamawaki et al., 2014). White pixels in every image (axonal expression or reconstructed cell) were counted, and their total counts per row were divided by the total numbers of white pixels in the image, resulting in a projection density distribution vector along the layers of RSG. Average layer widths were calculated from the previously acquired measurements and used to set the number of bins per layer to uniformly represent layers. Density plots were then visualized from the adjusted distribution vectors using the following Python modules: pandas, matplotlib, scipy, and Seaborn. For correlation analysis, area from pia to layer 5A, corresponding to the end of RS basal dendrites, was chosen. Spearman’s correlation matrices were created using the pandas and scipy Python modules and plotted as a heatmap. Images of axonal projection patterns to prefrontal cortex (PFC) (Figure 6—figure supplement 2) were obtained from the Allen Mouse Brain Connectivity Atlas (https://connectivity.brain-map.org; Oh et al., 2014).
Modeling methodology
Our model circuit consists of an RSG LR cell receiving input from an ensemble of thalamic HD cells, via synapses exhibiting short-term depression.
Modeling the presynaptic HD population
Request a detailed protocolThroughout, let denote the head direction at time ms. Each HD cell has a preferred angle (PA) at which its firing rate is maximized. The firing rate at time of an HD cell with PA takes the form:
where is a fixed function describing the shape of the cell’s tuning curve. In all simulations, we utilize a Gaussian tuning curve:
Here, is the background firing rate, is the maximum firing rate, and sets the tuning curve width. We draw these parameters from Beta distributions with statistics in accordance with the ADN HD cell statistics reported in Taube, 2010; see Supplementary file 1 - Table 3 for exact quantities. The spike train of a single HD cell with PA follows an inhomogeneous Poisson process, with intensity given by at time . We additionally impose a post-spike refractory period of 4ms.
We draw these parameters from Beta distributions with statistics in accordance with the ADN HD cell statistics reported in Taube, 2010 (see Supplementary file 1 - Table 3 for exact quantities). The probability density function of the Beta distribution with parameters is given by with
where denotes the standard Gamma function.
Finally, the spike train of a single HD cell with PA follows an inhomogeneous Poisson process, with intensity given by at time . We additionally impose a post-spike refractory period of 4ms.
Anticipatory firing
Request a detailed protocolThalamic HD cells exhibit an anticipatory time interval, or ATI: for each cell, there exists some constant duration Ams such that the cell tends to spike Ams prior to when the animal’s heading equals the preferred direction of the cell. We define this notion more precisely following a previous study (Blair et al., 1997). Restrict attention to a specific HD cell, with ATI Ams, and let denote its PA when the angular head velocity is zero. The ATI is defined via the property that when angular head velocity equals , the PA of this cell shifts from to . Equivalently, we can think of HD cell PAs as static and instead say anticipatory firing effectively transforms the heading trajectory itself: if describes the true angular trajectory, the effective trajectory experienced by our chosen HD cell is:
In vivo head direction recordings
Request a detailed protocolTo simulate the angular trajectory , we employ in vivo mouse head tracking recordings from the data set CRCNS-th1 (Peyrache et al., 2015). These in vivo heading recordings were originally sampled at a relatively low rate of approximately 40 Hz; in order to facilitate use of this data at the finer time scale of HD cell spiking, we linearly interpolated the raw data. Angular head velocity was computed using the standard centered approximation to the derivative:
Modeling short-term synaptic plasticity
Request a detailed protocolIn the Tsodyks-Markram (TM) model (Tsodyks et al., 1998) of short-term synaptic plasticity, each synapse has a finite pool from which it releases vesicles into the cleft upon arrival of a presynaptic spike. The state of the synapse at any moment in time is characterized by the triple , specified as:
: fraction of vesicles remaining in the pool at time .
: fraction of pool released into cleft upon presynaptic spike arrival at time (i.e. release probability).
: synaptic conductance at time .
Here, is measured in milliseconds, in umho, and the variables and are dimensionless. Depression results from depletion of due to a spike train, followed by slow recovery of the vesicle pool; facilitation results from increase in due to a spike train, followed by slow decay of this release probability. Formally, these dynamical variables satisfy the system:
Here, denotes an arbitrary presynaptic spike time, and the term specifies a discontinuous modification to the variable upon arrival of a presynaptic spike. We use the notation , and thus and satisfy the relation . The parameters , , and set the timescale of depression recovery, facilitation decay, and synaptic conductance decay, respectively. Note that always, so represents the minimum fraction of contributing to the conductance on each presynaptic spike arrival. Finally, the parameter is just a tunable synaptic weight (with units of umho).
In the regime of , depression dominates the synaptic dynamics; conversely, in the regime of , facilitation dominates. Accordingly, the limits and describe depression-only and facilitation-only dynamics, respectively. On the other hand, when and are on the same order of magnitude, the synapse may exhibit combined depressing-facilitating dynamics.
Synapse implementation and parameter choices
Request a detailed protocolWe utilized an implementation of the Tsodyks-Markram model as a NEURON mechanism, freely available for download from Model DB (https://senselab.med.yale.edu/ModelDB/showmodel.cshtml?model=3815). For all simulations, we use the parameters:
in the TM model. We obtained this parameter set by fitting the TM model response to recorded EPSP amplitudes evoked via thalamocortical stimulation of RSG LR cells with a 10 Hz pulse train (Figures 2 and 7). This was accomplished by searching parameter space exhaustively for the parameter set minimizing a least-squares loss function. The resultant TM dynamics were moderately depressing, in accordance with the experimental results. Non-depressing synapses were modeled as standard exponential synapses with instantaneous rise time.
Modeling the postsynaptic LR population
Request a detailed protocolWe model the postsynaptic LR cell using the morphologically detailed model originally presented in Brennan et al., 2020. This model, implemented in NEURON (Carnevale and Hines, 2006), consists of Hodgkin-Huxley conductances representing a fast sodium current, a delayed rectifier potassium current, and a Kv1.1-mediated slow potassium current. We place each HD input → LR synapse at a random location on a distal apical dendrite of the LR cell; varying this synaptic placement does not qualitatively change results (data not shown).
PA distribution of a postsynaptic cell
Request a detailed protocolEach postsynaptic cell receives input from presynaptic HD cells, each of which has a unique PA. When studying postsynaptic speed coding, we initially drew the PA of each HD cell uniformly at random from . Later on, we examined the effect of introducing nonuniform PA distributions. In this case, we assumed that the PA distribution of each LR cell was continuous as a function of angle, and thus modeled the PA distribution probability densities randomly as resulting from Brownian bridge processes (to ensure periodicity of the density). We then used inverse transform sampling to draw from the generated distributions.
Data analysis
Request a detailed protocolTo quantify the extent to which postsynaptic activity reflects head turning speed, we utilized two similarity measures : Pearson correlation and mutual information. In fact, we used the time-lagged counterparts of these measures: for a pair of time series , , their -lagged similarity measure is . Varying , we obtain the cross-correlation and cross-mutual information between the two variables. Throughout the paper, when referring to an instance of measure between time series , , we either report the function or just the maximum value of the function:
In order to compute mutual information, we discretized the time series for head turning speed and postsynaptic firing rate into 40 bins each and calculated the mutual information between these discrete proxies.
Appendix 1
Mean-field model analysis
In order to better understand the observed speed tuning effects, we study an analytically tractable mean-field model of the simulation circuit. We show for this model that the aggregate postsynaptic activity does indeed reflect head turning speed – in fact, at low speeds it is proportional to the squared head speed. Moreover, we show that anticipatory firing compensates for lag introduced by the depression recovery timescale, enabling sharper and lower latency speed encoding. Finally, we briefly comment on model predictions when the HD synaptic input exhibits facilitating or combined facilitating-depressing dynamics.
Model description
Let denote the true head direction at time . For our mean-field analysis, we model all HD cells sharing preferred angle (PA) as a lumped subpopulation, whose firing rate at time is given by for the tuning curve function
Here denotes the indicator function of the interval . Note that this is not the same tuning function as utilized in simulations; here, the subpopulation of HD cells with PA is active precisely at times when the current head direction is within of . This simplification, in which HD cell tuning curves are discontinuous step functions rather than continuous Gaussians, gives qualitatively analogous results to the Gaussian case (data not shown). In the mean-field model, we further assume the background firing rate , maximum firing rate , and tuning curve width are constants independent of .
The synapses follow simplified Tsodyks-Markram (TM) dynamics (Tsodyks et al., 1998), in which the dynamical quantities represent population-averaged values of the standard TM synaptic variables (see Materials and methods for details of the standard TM model) over all synapses corresponding to HD cells having a particular . In what follows, we explicitly model synaptic depression only, although similar equations describe the full TM model with combined depression and facilitation. Thus, we are interested in the quantity , representing the average value of the depression variable for all synapses corresponding to HD cells with PA , at time . The dynamics of are given by the equation:
Here denotes the time constant of recovery from depression, and represents the synaptic release probability. When , the HD cells with PA are very depressed; when the firing rate lets up, recovers to its steady-state value of 1. Similarly, denotes the average conductance of all synapses corresponding to HD cells with PA , at time . We consider the limit in which the timescale of synaptic decay is much shorter than the depression recovery timescale , and thus we can make the approximation .
Now, let denote the density function of the presynaptic PA distribution. We are interested in dynamics of the -averaged quantities
In particular, is proportional to the sum of over all , i.e. the total conductance of the postsynaptic cell, and thus reflects the activity level of the postsynaptic cell.
Equation for dynamics of
In this section, we derive an explicit equation for the dynamics of in the limit of low background firing (). Crucially, the choice of step function tuning enables us to close the relevant system of ODEs to obtain low-dimensional dynamics. Although we will not study this case in any further depth, we note that when is nonneglible, the quantities and form a coupled planar system.
Proposition
In the limit of low background firing (), the mean synaptic conductance satisfies the equation
Here we have
and , with defined as the average firing rate of a presynaptic HD cell at time , i.e.
Proof
From the defining equation (Equation 2) for , we see that satisfies
We have for fixed that
For our specific choice of step function , we have . Setting and , we thus have
Now averaging over gives
Setting and , we can rewrite this as
where denotes the Dirac delta function and . It follows that
Finally setting gives the proposition.
Analyzing Contributions To
In the limit, the mean conductance is driven by fluctuations in , between which it relaxes to the value with time constant . Accordingly, to understand the behavior of we should study the behavior of and . For all of what follows, we assume .
Baseline activity:
When , i.e. the true head direction is not changing, we have following a transient. Thus, we can interpret as the baseline activity level of the postsynaptic cell. In the case , i.e. that of a uniform presynaptic PA distribution, the baseline activity is constant in time, since:
On the other hand, if is nonuniform, then the integral defining is not constant in time, taking larger values when is at angles around which the PA distribution has more mass. Since , exhibits the same effect, thereby explaining the observed modulation of baseline firing rate by head direction.
Contribution of depression:
We now assume in all of what follows that the PA distribution is uniform, i.e. . In this case, the function reduces to
We can interpret as the contribution of synaptic depression to , since it is the only term in the equation for that incorporates any effect of the depression variable . To understand the behavior of , we differentiate to obtain
where and
Note that for our step function tuning, is not actually defined. However, we can treat as if it is continuous at . One way to justify this is to imagine has been replaced with a continuous function that is close to in the sense that
is small; it is always possible to make this quantity arbitrarily small with a suitable choice of continuous .
The equation above further reduces the problem of understanding to that of understanding , which distorts the velocity information appearing in (Equation 7). To proceed, we make the approximation , where is the time-average of . This approximation will be valid if is typically positive and not too variable about its mean. Indeed, at each moment in time, we have a depression profile describing, for PA , the average depression of synapses from the HD subpopulation having PA . The graph of this profile has a well-like shape, centered close to . At low speeds, the slope will be positive, and the slope will be negative; it follows that the difference is itself typically positive, becoming smaller only during high speed head turns. For now, we assume also has low variability, and that therefore the approximation holds; later on, we will discuss the validity of this assumption. Continuing with the approximation, we have the simplified equation
It follows that tracks the variable , i.e, we have the approximate relation . Since is driven by the product , it follows that approximately . Thus, when the approximation for is valid, we conclude the postsynaptic conductance is proportional to the square of angular head speed.
However, these proportionality relations are only approximate, and are degraded by the integration times in (Equation 9) and in (Equation 4). In the next section, we will see that anticipatory firing compensates for these lags and sharpens speed coding.
Anticipatory firing improves speed coding
As mentioned above, the lags and degrade the approximate relation . Indeed, as shown in Figure 10, at ATI = 0 the postsynaptic activity is not especially informative about head speed, relative to the ATI = 50 ms and ATI = 100 ms cases. Intuitively, a positive ATI compensates for these lags; to make this statement precise, we note the following lemma, for which we omit the straight forward proof.
Lemma
Let be a differentiable function with . If , then the solution to the initial value problem
has the property that following a transient (more precisely, exponentially fast as )
—
Recall that introducing an ATI of A ms is equivalent to replacing the true angular trajectory with the effective trajectory . Applying the lemma to (Equation 9) with and , we see that when the ATI equals ms, the relation becomes exact (with a proportionality constant ). Moreover, in this case, the equation for simplifies to
Since we have assumed is uniform, and hence is constant in time, we can apply the lemma again taking . If we assume , then there is the exact equality
Thus, the condition for exact speed coding is
Finally, we note that in (Equation 8) the most parsimonious assumption to make is that ; this amounts to assuming we have replaced the step function tuning with a continuous function whose steepest portions are centered at and . It follows that indeed , and therefore, it is possible to be in the regime of exact speed coding. This analysis shows that the ATI optimizing speed coding in our model is . For the parameter set used in simulations (τ = 270ms, U = 0.28ms, fmax = 0.07 spikes/ms), we find that
Accordingly, the optimal ATI is . Indeed, simulations utilizing the above parameter set and step-function tuning curves exhibit optimal speed coding at an ATI of precisely 85 ms, thereby validating the analysis (Figure 8—figure supplement 1).
Scaling of postsynaptic activity with speed
As mentioned earlier, validity of the approximation depends on the variability of about its mean. When this approximation is valid, we have the scaling . Before describing the empirical validity of this approximation, we first note that even when is fairly variable, the function remains strongly correlated with , and consequently the analysis of ATI in the previous section remains approximately valid. Indeed, in the case of step-function tuning, we find that is actually fairly variable (CV = 0.6), but postsynaptic activity is still very strongly correlated with speed at the optimal ATI computed by the formula above (Figure 9—figure supplement 1). However, the distortion by of the velocity information in results in a speed-postsynaptic activity relation whose best-fit power law exponent is approximately 1.5-1.7, instead of the predicted value of 2. Thus, for step function tuning, the quadratic scaling does not hold.
On the other hand, for Gaussian tuning curves, the approximation is valid. Although the above derivations have assumed step-function tuning, we may heuristically imagine that the equations are approximately valid even with Gaussian tuning, for some suitable choice of tuning curve width . Proceeding as such, we find that with Gaussian tuning, is relatively invariable about its mean (CV = 0.2), especially for low head speeds. We note that the tuning curve shape directly affects the shape of the depression profile , and consequently affects the variability of . This suggests that for low speeds and Gaussian tuning, the scaling will be accurate. Mean-field simulation results do yield that for low speeds, the best-fit power law exponent is approximately 1.9-2.1.
This raises the question of whether this scaling is observable in data from simulations of our spiking model, since spiking simulations utilize the Gaussian tuning, which is more accurate to the experimentally observed shape of HD cell tuning curves. We fit results from the spiking model (postsynaptic firing rates binned by head speed) to a power law when including data points in the speed-firing rate curve only up to deg/s, for varying values of (Figure 8—figure supplement 2). We found that best-fit residuals were lowest for low speeds and that the best-fit exponents plateaued in a neighborhood of 2 before falling off as the considered range of speeds widened, suggesting that the approximate quadratic scaling is potentially visible in experimental data in this circuit and others utilizing the same synaptic-depression-based mechanism for computing rate of change.
Data availability
The LR neuron model utilized here is available on ModelDB (https://modeldb.yale.edu/260192), and the implementation of the Tsodyks-Markram model utilized here is also available from ModelDB (https://senselab.med.yale.edu/ModelDB/showmodel.cshtml?model=3815). Reconstructions will be available at http://neuromorpho.org/KeywordSearch.jsp and can be found by searching for the "Ahmed" archive. Source experimental data for figures 1-5, 7, Figure 1-Figure supplement 1, and Supplementary File 1-Table 1 are provided.
References
-
The hippocampal rate code: anatomy, physiology and theoryTrends in Neurosciences 32:329–338.https://doi.org/10.1016/j.tins.2009.01.009
-
Egocentric boundary vector tuning of the retrosplenial cortexScience Advances 6:eaaz2322.https://doi.org/10.1126/sciadv.aaz2322
-
The boundary vector cell model of place cell firing and spatial memoryReviews in the Neurosciences 17:71.https://doi.org/10.1515/REVNEURO.2006.17.1-2.71
-
Neuronal vector coding in spatial cognitionNature Reviews Neuroscience 21:453–470.https://doi.org/10.1038/s41583-020-0336-9
-
Conflict monitoring and decision making: reconciling two perspectives on anterior cingulate functionCognitive, Affective, & Behavioral Neuroscience 7:356–366.https://doi.org/10.3758/CABN.7.4.356
-
New breakthroughs in understanding the role of functional interactions between the neocortex and the claustrumThe Journal of Neuroscience 37:10877–10881.https://doi.org/10.1523/JNEUROSCI.1837-17.2017
-
Memory for events and their spatial context: models and experimentsPhilosophical Transactions of the Royal Society of London. Series B: Biological Sciences 356:1493–1503.https://doi.org/10.1098/rstb.2001.0948
-
Subcortical projections to the frontal pole in the marmoset monkeyEuropean Journal of Neuroscience 34:303–319.https://doi.org/10.1111/j.1460-9568.2011.07744.x
-
BookThe NEURON BookCambridge: Cambridge University Press.https://doi.org/10.1017/CBO9780511541612
-
Coordinated activities of retrosplenial ensembles during resting-state encode spatial landmarksPhilosophical Transactions of the Royal Society B: Biological Sciences 375:20190228.https://doi.org/10.1098/rstb.2019.0228
-
Brain-wide maps of fos expression during fear learning and recallLearning & Memory 24:169–181.https://doi.org/10.1101/lm.044446.116
-
Heterogeneity in human retrosplenial cortex: a review of function and connectivityBehavioral Neuroscience 132:317–338.https://doi.org/10.1037/bne0000261
-
Unveiling the diversity of thalamocortical neuron subtypesEuropean Journal of Neuroscience 35:1524–1532.https://doi.org/10.1111/j.1460-9568.2012.08033.x
-
What is the function of the claustrum?Philosophical Transactions of the Royal Society B: Biological Sciences 360:1271–1279.https://doi.org/10.1098/rstb.2005.1661
-
Thalamic control of layer 1 circuits in prefrontal cortexJournal of Neuroscience 32:17813–17823.https://doi.org/10.1523/JNEUROSCI.3231-12.2012
-
The synaptic organization of the claustral projection to the cat's visual cortexJournal of Neuroscience 30:13166–13170.https://doi.org/10.1523/JNEUROSCI.3122-10.2010
-
Thalamocortical transformations of periodic stimuli: the effect of stimulus velocity and synaptic short-term depression in the vibrissa-barrel systemJournal of Computational Neuroscience 25:122–140.https://doi.org/10.1007/s10827-007-0068-0
-
Are the boundary-related cells in the subiculum boundary-vector cells?Journal of Neuroscience 29:13429–13431.https://doi.org/10.1523/JNEUROSCI.4176-09.2009
-
Parahippocampal and retrosplenial contributions to human spatial navigationTrends in Cognitive Sciences 12:388–396.https://doi.org/10.1016/j.tics.2008.07.004
-
Attention: the claustrumTrends in Neurosciences 38:486–495.https://doi.org/10.1016/j.tins.2015.05.006
-
Anterolateral motor cortex connects with a medial subdivision of ventromedial thalamus through cell Type-Specific circuits, forming an excitatory Thalamo-Cortico-Thalamic loop via layer 1 apical tuft dendrites of layer 5B pyramidal tract type neuronsThe Journal of Neuroscience 38:8787–8797.https://doi.org/10.1523/JNEUROSCI.1333-18.2018
-
Independence of diffuse thalamo-cortical projection system shown by specific nuclear destructionsJournal of Neurophysiology 16:252–271.https://doi.org/10.1152/jn.1953.16.3.252
-
Morphology and physiology of cortical neurons in layer IThe Journal of Neuroscience 16:5290–5300.https://doi.org/10.1523/JNEUROSCI.16-17-05290.1996
-
Neural mechanisms of navigation involving interactions of cortical and subcortical structuresJournal of Neurophysiology 119:2007–2029.https://doi.org/10.1152/jn.00498.2017
-
Pyramidal cell communication within local networks in layer 2/3 of rat neocortexThe Journal of Physiology 551:139–153.https://doi.org/10.1113/jphysiol.2003.044784
-
Meta-analysis of functional brain imaging in specific phobiaPsychiatry and Clinical Neurosciences 67:311–322.https://doi.org/10.1111/pcn.12055
-
Dynamics of place, boundary and object encoding in rat anterior claustrumFrontiers in Behavioral Neuroscience 9:250.https://doi.org/10.3389/fnbeh.2015.00250
-
The thalamic matrix and thalamocortical synchronyTrends in Neurosciences 24:595–601.https://doi.org/10.1016/S0166-2236(00)01922-6
-
Dendritic spikes are enhanced by cooperative network activity in the intact HippocampusThe Journal of Neuroscience 18:3919–3928.https://doi.org/10.1523/JNEUROSCI.18-10-03919.1998
-
Neurotoxic lesions of retrosplenial cortex disrupt signaled and unsignaled contextual fear conditioningBehavioral Neuroscience 122:1070–1077.https://doi.org/10.1037/a0012895
-
Synaptic organization of the neuronal circuits of the claustrumThe Journal of Neuroscience 36:773–784.https://doi.org/10.1523/JNEUROSCI.3643-15.2016
-
Robust information routing by dorsal subiculum neuronsScience Advances 7:eabf1913.https://doi.org/10.1126/sciadv.abf1913
-
Organization of the Claustrum-to-Entorhinal cortical connection in miceThe Journal of Neuroscience 37:269–280.https://doi.org/10.1523/JNEUROSCI.1360-16.2016
-
Pyramidal neurons in the superficial layers of rat retrosplenial cortex exhibit a late-spiking firing propertyBrain Structure and Function 218:239–254.https://doi.org/10.1007/s00429-012-0398-1
-
The retrosplenial cortex is involved in the formation of memory for context and trace fear conditioningNeurobiology of Learning and Memory 123:110–116.https://doi.org/10.1016/j.nlm.2015.06.007
-
Boundary vector cells in the subiculum of the hippocampal formationThe Journal of Neuroscience 29:9771–9777.https://doi.org/10.1523/JNEUROSCI.1319-09.2009
-
Retrosplenial cortex lesions impair acquisition of active avoidance while sparing fear-based emotional memoryBehavioural Brain Research 173:229–236.https://doi.org/10.1016/j.bbr.2006.06.026
-
The retrosplenial contribution to human navigation: a review of lesion and neuroimaging findingsScandinavian Journal of Psychology 42:225–238.https://doi.org/10.1111/1467-9450.00233
-
Spatial memory Engram in the mouse retrosplenial cortexCurrent Biology 28:1975–1980.https://doi.org/10.1016/j.cub.2018.05.002
-
A study of thalamo-cortical relationsAmerican Journal of Physiology-Legacy Content 135:281–292.https://doi.org/10.1152/ajplegacy.1941.135.2.281
-
The claustrum coordinates cortical slow-wave activityNature Neuroscience 23:741–753.https://doi.org/10.1038/s41593-020-0625-7
-
Subiculum neurons map the current Axis of travelNature Neuroscience 20:170–172.https://doi.org/10.1038/nn.4464
-
Neurogliaform cells in cortical circuitsNature Reviews Neuroscience 16:458–468.https://doi.org/10.1038/nrn3969
-
Thalamic input to the cerebral cortexTrends in Neurosciences 2:183–185.https://doi.org/10.1016/0166-2236(79)90074-2
-
SoftwareCRCNS.org. Extracell Rec From Multi-Site Silicon Probes Anterior Thalamus Subicular Form Free Mov MiceCRCNS.org. Extracell Rec From Multi-Site Silicon Probes Anterior Thalamus Subicular Form Free Mov Mice.
-
Callosal responses in a retrosplenial columnBrain Structure & Function 223:1051–1069.https://doi.org/10.1007/s00429-017-1529-5
-
Organization of anterior cingulate and frontal cortical projections to the retrosplenial cortex in the ratThe Journal of Comparative Neurology 506:30–45.https://doi.org/10.1002/cne.21523
-
Burst firing and spatial coding in subicular principal cellsThe Journal of Neuroscience 39:3651–3662.https://doi.org/10.1523/JNEUROSCI.1656-18.2019
-
Rat claustrum coordinates but does not integrate somatosensory and motor cortical informationJournal of Neuroscience 32:8583–8588.https://doi.org/10.1523/JNEUROSCI.1524-12.2012
-
Hypotheses relating to the function of the claustrumFrontiers in Integrative Neuroscience 6:53.https://doi.org/10.3389/fnint.2012.00053
-
Anterior cingulate cortex: unique role in cognition and emotionThe Journal of Neuropsychiatry and Clinical Neurosciences 23:121–125.https://doi.org/10.1176/jnp.23.2.jnp121
-
Boundary coding in the rat subiculumPhilosophical Transactions of the Royal Society B: Biological Sciences 369:20120514.https://doi.org/10.1098/rstb.2012.0514
-
Contributions of cingulate cortex to two forms of spatial learning and memoryThe Journal of Neuroscience 8:1863–1872.https://doi.org/10.1523/JNEUROSCI.08-06-01863.1988
-
Head direction cells recorded in the anterior thalamic nuclei of freely moving ratsThe Journal of Neuroscience 15:70–86.https://doi.org/10.1523/JNEUROSCI.15-01-00070.1995
-
Head direction cells and the neurophysiological basis for a sense of directionProgress in Neurobiology 55:225–256.https://doi.org/10.1016/S0301-0082(98)00004-5
-
Interspike interval analyses reveal irregular firing patterns at short, but not long, intervals in rat head direction cellsJournal of Neurophysiology 104:1635–1648.https://doi.org/10.1152/jn.00649.2009
-
Persistent neural activity in head direction cellsCerebral Cortex 13:1162–1172.https://doi.org/10.1093/cercor/bhg102
-
Neural networks with dynamic synapsesNeural Computation 10:821–835.https://doi.org/10.1162/089976698300017502
-
Spontaneous electrical activity of the claustrumNeurophysiology 12:113–120.https://doi.org/10.1007/BF01065305
-
Connections of the retrosplenial granular a cortex in the ratThe Journal of Comparative Neurology 300:593–606.https://doi.org/10.1002/cne.903000412
-
Projections from the anterodorsal and anteroventral nucleus of the thalamus to the limbic cortex in the ratThe Journal of Comparative Neurology 358:584–604.https://doi.org/10.1002/cne.903580411
-
Connections of the retrosplenial granular b cortex in the ratThe Journal of Comparative Neurology 463:249–263.https://doi.org/10.1002/cne.10757
-
Cortical connections between rat cingulate cortex and visual, motor, and postsubicular corticesThe Journal of Comparative Neurology 216:192–210.https://doi.org/10.1002/cne.902160207
-
Form and distribution of neurons in rat cingulate cortex: areas 32, 24, and 29The Journal of Comparative Neurology 195:603–625.https://doi.org/10.1002/cne.901950406
-
Organization of the connections between claustrum and cortex in the mouseJournal of Comparative Neurology 525:1317–1346.https://doi.org/10.1002/cne.24047
-
Frontal cortical control of posterior sensory and association cortices through the claustrumBrain Structure and Function 223:2999–3006.https://doi.org/10.1007/s00429-018-1661-x
-
Dendritic bundling in layer I of granular retrosplenial cortex: intracellular labeling and selectivity of innervationThe Journal of comparative neurology 295:33–42.https://doi.org/10.1002/cne.902950104
-
A Corticocortical Circuit Directly Links Retrosplenial Cortex to M2 in the MouseThe Journal of Neuroscience 36:9365–9374.https://doi.org/10.1523/JNEUROSCI.1099-16.2016
-
Combining Optogenetics and Electrophysiology to Analyze Projection Neuron CircuitsCold Spring Harbor Protocols 2016:pdb.prot090084–.090847.https://doi.org/10.1101/pdb.prot090084
-
Heterogeneity of neuronal firing type and morphology in retrosplenial cortex of male F344 ratsJournal of Neurophysiology 123:1849–1863.https://doi.org/10.1152/jn.00577.2019
-
Input-output organization of the mouse claustrumJournal of Comparative Neurology 526:2428–2443.https://doi.org/10.1002/cne.24502
Article and author information
Author details
Funding
University of Michigan (Startup Funds)
- Omar J Ahmed
National Institute of Neurological Disorders and Stroke (NS121745)
- Omar J Ahmed
Whitehall Foundation (Research Grant)
- Omar J Ahmed
National Institute on Aging (5P30AG053760)
- Omar J Ahmed
National Science Foundation (Graduate Research Fellowship)
- Ellen KW Brennan
National Institute of Neurological Disorders and Stroke (T-32-NS076401)
- Ellen KW Brennan
- Sharena P Rice
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Scott Cruikshank, Roger Albin, and Clare Weiland for their comments on the manuscript, Vaughn Hetrick for his technical assistance, Michael Roberts for his help with the NeuroLucida software, Victoria Booth for her assistance with the modeling work, and Danny Siu for his assistance with MATLAB-based analyses. This work was supported by lab startup funds from the University of Michigan (OJA), a grant from the Whitehall Foundation (OJA), NIH grant NS121745 (OJA), pilot funds from the NIH/NIA funded Michigan Alzheimer’s Disease Research Center 5P30AG05376 (OJA), NIH T-32-NS076401 (EKWB and SPR), and by an NSF graduate research fellowship (EKWB).
Ethics
Animal experimentation: All procedures and use of animals were approved by the University of Michigan Institutional Animal Care and Use Committee (Protocol #9572). This study was performed in strict adherence with guidelines by the Guide for the Care and Use of Laboratory Animals and the Institutional Animal Care and Use Committee. Surgical procedures were performed under isoflurane anesthesia, and analgesics and other care was used to minimize animal discomfort. Slicing protocols also used deep isoflurane anesthesia to eliminate any pain or discomfort felt by the animal.
Copyright
© 2021, Brennan et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,071
- views
-
- 458
- downloads
-
- 34
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
Diffusional kurtosis imaging (DKI) is a methodology for measuring the extent of non-Gaussian diffusion in biological tissue, which has shown great promise in clinical diagnosis, treatment planning, and monitoring of many neurological diseases and disorders. However, robust, fast, and accurate estimation of kurtosis from clinically feasible data acquisitions remains a challenge. In this study, we first outline a new accurate approach of estimating mean kurtosis via the sub-diffusion mathematical framework. Crucially, this extension of the conventional DKI overcomes the limitation on the maximum b-value of the latter. Kurtosis and diffusivity can now be simply computed as functions of the sub-diffusion model parameters. Second, we propose a new fast and robust fitting procedure to estimate the sub-diffusion model parameters using two diffusion times without increasing acquisition time as for the conventional DKI. Third, our sub-diffusion-based kurtosis mapping method is evaluated using both simulations and the Connectome 1.0 human brain data. Exquisite tissue contrast is achieved even when the diffusion encoded data is collected in only minutes. In summary, our findings suggest robust, fast, and accurate estimation of mean kurtosis can be realised within a clinically feasible diffusion-weighted magnetic resonance imaging data acquisition time.
-
- Neuroscience
The retina transforms patterns of light into visual feature representations supporting behaviour. These representations are distributed across various types of retinal ganglion cells (RGCs), whose spatial and temporal tuning properties have been studied extensively in many model organisms, including the mouse. However, it has been difficult to link the potentially nonlinear retinal transformations of natural visual inputs to specific ethological purposes. Here, we discover a nonlinear selectivity to chromatic contrast in an RGC type that allows the detection of changes in visual context. We trained a convolutional neural network (CNN) model on large-scale functional recordings of RGC responses to natural mouse movies, and then used this model to search in silico for stimuli that maximally excite distinct types of RGCs. This procedure predicted centre colour opponency in transient suppressed-by-contrast (tSbC) RGCs, a cell type whose function is being debated. We confirmed experimentally that these cells indeed responded very selectively to Green-OFF, UV-ON contrasts. This type of chromatic contrast was characteristic of transitions from ground to sky in the visual scene, as might be elicited by head or eye movements across the horizon. Because tSbC cells performed best among all RGC types at reliably detecting these transitions, we suggest a role for this RGC type in providing contextual information (i.e. sky or ground) necessary for the selection of appropriate behavioural responses to other stimuli, such as looming objects. Our work showcases how a combination of experiments with natural stimuli and computational modelling allows discovering novel types of stimulus selectivity and identifying their potential ethological relevance.