Thalamus and claustrum control parallel layer 1 circuits in retrosplenial cortex

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.


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.
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 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 'lowrheobase' 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. Figure 1. Anterior thalamic input controls LR, but not RS, neurons in superficial granular retrosplenial cortex. (A) Representative example of the firing properties of LR neurons. Top trace: sustained high-frequency firing of action potentials with little spike frequency adaptation at high current inputs in response to 200 pA current injection. Bottom trace: Delay to first spike at threshold current inputs (50 pA) and little to no sag when hyperpolarized (À100 pA). Right inset is a zoomed in view of the first spike in response to a near-threshold current input. Cell is held at À65 mV resting potential. (B) Figure 1 continued on next page 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 mm 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 Figure 1 continued Same as A, but for a representative RS neuron. Note the spike frequency adaptation with high current inputs, delay to first spike at near-threshold current inputs, and sag potential seen with hyperpolarization. Cell is held at À65 mV resting potential. (C) Population data comparisons between LR and RS neurons of the following intrinsic properties: input resistance (p=1.31e-26; Wilcoxon rank sum), rheobase (p=1.24e-17; Wilcoxon rank sum), and adaptation ratio (p=2.94e-24; Wilcoxon rank sum). Error bars are standard error of the mean (SEM). (D) Left: Schematic showing the injection of optogenetic viral construct into the anterodorsal/anteroventral (ADAV) nuclei of the anterior thalamus. Middle: DIC image of RSG. Right: EYFP fluorescence image of RSG showing expressing anterior thalamic axons and terminal arbors. (E) Schematic of optogenetic stimulation of TC axons and terminal arbors over the patched cell body in RSG L3. (F) LR (purple) and RS (blue) cell responses to 10 Hz light pulses (see Materials and methods) at maximum LED power. LR cell shows spiking responses, while the RS cell has very small EPSPs. Both cells were held at À65 mV before stimulation. Blue triangles indicate light pulses. (G) Optogenetic activation of thalamic axons in RSG increases the firing rate of the LR (purple) but not RS (blue) neuron in a simultaneously patched pair. Neurons were held in ACSF +picrotoxin at a constant firing rate of 10-30 Hz via a 2 s current injection, and a 1 ms LED pulse was delivered 500 ms into the spike train (see Materials and methods). Top: Raster plots for all sweeps of a representative LR (purple) and RS (blue) example. Middle: Firing trace of one sweep for that same LR and RS pair. Bottom: PSTH plots evaluating all sweeps for the same LR and RS example cells. (H) Population analysis of the spike ratio for all tested LR (n = 7) and RS (n = 6) cells showing a significant increase in firing rate post-LED pulse in LR cells compared to RS cells (p=0.0012; Wilcoxon rank sum test). Error bars are SEM. See Figure 1-source code 1 and Figure 1-source data 1 for MATLAB code and source data used to plot bar graphs in panel C and associated supplements. The online version of this article includes the following source data, source code and figure supplement(s) for figure 1: Source data 1. Cell-type specific intrinsic physiology. Source code 1. Code to generate intrinsic physiology comparisons.  Right: confocal image of layers 1-5 of RSG showing lamination of anterior thalamic axons and terminal arbors (green) and cell membrane marker, NeuroTrace (cyan). (C) Left: LR reconstruction superimposed on patch location in an RSG slice with thalamic projections (green) and NeuroTrace (cyan). Right: Projection density plot showing density of LR dendrites (n = 10; mean ± SEM shaded) in purple and density of ADAV axon expression (n = 5; mean ± SEM shaded) in green plotted as distance from the pia (mm). (D) Same as C for RS cells with dendrites (n = 5) plotted in blue. (E) Top: Schematic of L1a-targeted optogenetic stimulation. Bottom: Significantly larger EPSP amplitude at high LED power for LR (n = 15) cells compared to RS (n = 8) cells (p=0.0185; Wilcoxon rank sum test). Error bars are SEM. (F) Left: Confocal image of RSG slice with anterior thalamic axons and terminal arbors (green). Layers are demarcated, and schematic LR and RS neurons are placed in their representative locations within the superficial layers. Blue circles indicate targeting of minimum Figure 2 continued on next page 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 L5 sup (p=0.420), while L5 dp 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, L5 sup , and L5 dp 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 mM) or sodium channel blocker TTX (1 mM)+4 AP (100 mM) 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 Representative trace examples of an LR (purple) and RS (blue) neurons' responses to the LED stimulation at each location. Green shading is representative of the projection density of thalamic axons and terminal arbors across the shown layers. (G) Population analysis for LR and RS cells' normalized EPSP responses to minimum LED power optogenetic stimulations at L1a, L1/2 boundary, cell body in L3, superficial L5 (L5 sup ), and deep L5 (L5 dp ) locations. Top: LR cells have significantly larger responses at L1a (n = 58, p=0.00001, paired t-test) and cell body (n = 53, p=0.0006) compared to L1/2 (n = 59). Responses at L5 dp (n = 10, p=0.0462) are significantly smaller than L1/2, and responses at L5 sup (n = 12, p=0.420) do not significantly differ from L1/2. Bottom: EPSP amplitude in RS cells does not significantly differ at any stimulation location compared to L1/2 stimulation (L1a: n = 27, p=0.289; cell body: n = 25, p=0.702; L5 sup : n = 12, p=0.197; L5 dp : n = 11, p=0.261; L1/2: n = 27). Error bars are SEM. The paired t-test was used for all statistical comparisons. See Figure 2-source data 1 for source data. The online version of this article includes the following source data and figure supplement(s) for figure 2: Source data 1. Laminar-specific responses to ADAV input.  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 5 sup 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 L5 sup (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 mM) or sodium channel blocker TTX (1 mM)+4 AP (100 mM) 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 mM) or sodium channel blocker TTX (1 mM)+4 AP (100 mM) 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 neurons' responses to the LED stimulation at each location. Red shading is representative of the projection density of CLA axons and terminal arbors across the layers. (G) Population analysis for LR and RS normalized EPSP responses to minimum LED power stimulations at L1a, L1/2, cell body, L5 sup , and L5 dp locations. Top: EPSP amplitude in LR cells does not significantly differ at any stimulation location compared to L1a stimulation (L1/2: n = 24, p=0.5151; L3: n = 23, p=0.3276; L5 sup : n = 10, p=0.3434; L5 dp : n = 7, p=0.4348; L1a: n = 24). Bottom: RS cells have significantly larger responses at L1/2 (n = 18, p=0.00725) and L5 sup (n = 9, p=0.0301) compared to L1a (n = 18), while normalized response amplitudes at cell body (n = 17, p=0.1377) and L5 dp (n = 8, p=0.1355) did not differ from responses at L1a. The paired t-test was used for all statistical analyses. Error bars are SEM. See Figure 3source data 1 for source data. The online version of this article includes the following source data for figure 3: Source data 1. Laminar-specific responses to CLA input. 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.
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 anticorrelated (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- are demarcated, and schematic LR and RS neurons are placed in their representative locations within the superficial layers. Blue circles indicate targeting of the minimum LED power stimulations, beginning at midline and extending into layer 5, stimulating every 20-30 mm (see Materials and methods). Right: Representative trace examples of an LR (purple) and RS (blue) neurons' responses to the LED stimulation at each location. Magenta shading is representative of the projection density of ACC axons and terminal arbors across the layers. (G) Population analysis for LR and RS normalized EPSP responses to minimum LED power stimulations at L1a, L1/2, cell body, L5 sup , and L5 dp locations. Top: EPSP amplitude in LR cells does not significantly differ at any stimulation location compared to L1a stimulation (L1/2: n = 19, p=0.1126; cell body: n = 19, p=0.35; L5 sup : n = 13, p=0.3423; L5 dp : n = 13, p=0.7166; L1a: n = 19). Bottom: RS cells have significantly larger responses at L1/2 (n = 14, p=0.000037) and cell body (n = 14, p=0.000525) compared to L1a (n = 14) and did not significantly differ at L5 sup (n = 11, p=0.0607) or L5 dp (n = 10, p=0.1340) compared to L1a. The paired t-test was used for all statistical analyses. Error bars are SEM. See Figure 4-source data 1 for source data. The online version of this article includes the following source data for figure 4: Source data 1. Laminar-specific responses to ACC input. Blue circles indicate targeting of the minimum LED power stimulations, beginning at midline and extending into layer 5, stimulating every 20-30 mm Figure 5 continued on next page 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 6figure 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 (

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  (blue) neurons' responses to the LED stimulation at each location. Yellow shading is representative of the projection density of DS axons and terminal arbors across the layers. (G) Population analysis for LR and RS normalized EPSP responses to minimum LED power stimulations at L1a, L1/2, cell body, L5 sup , and L5 dp locations. Top: LR cells have significantly larger responses at L1/2 (n = 8, p=0.0132) and cell body (n = 8, p=0.0135) compared to responses at L1a (n = 8), while responses at L5 sup (n = 8, p=0.9146) and L5 dp (n = 7, p=0.4902) do not differ from those at L1a. Bottom: EPSP amplitude in RS cells does not significantly differ at any stimulation location compared to L1a stimulation (L1/2: n = 7, p=0.8811; cell body: n = 7, p=0.3304; L5 sup : n = 7, p=0.2260; L5 dp : n = 6, p=0.0816; L1a: n = 7). The paired t-test was used for all statistical analyses. Error bars are SEM. See Figure 5-source data 1 for source data. The online version of this article includes the following source data for figure 5: Source data 1. Laminar-specific responses to DS input. 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 8figure 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 9figure 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 Figure 6 continued terminal arbors in RSG. Note distinct sublaminar distribution of axons from ADAV, ACC, and CLA in layer 1. Right: Projection density of LR (purple; n = 10) and RS (blue; n = 5; mean ± SEM shaded for all) neurons. Note the distinct difference between LR and RS dendrite lamination across the layers. (C) Correlation matrix of means for all axon-axon, axon-dendrite, and dendrite-dendrite comparisons. Note that LR dendrites are significantly positively correlated with ADAV and DS axons but negatively correlated with CLA and ACC axons, while RS dendrites are significantly positively correlated with CLA and ACC axons but negatively correlated with ADAV and DS axons. (D) Summary schematic showing selective control of LR neurons by inputs from ADAV and DS, including head direction signals. In contrast, RS neurons are preferentially controlled by inputs from the CLA and ACC. This precise organization of principal cell dendrites and afferent axons forms two parallel circuits in superficial RSG. The online version of this article includes the following figure supplement(s) for figure 6:  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   We modeled a heterogeneous population of HD cells providing input to a single LR neuron via either depressing (green) or non-depressing (orange) synapses. Presynaptic HD cells varied in their tuning width, maximum firing rate, background firing rate, and anticipatory time interval (ATI). Tuning curves of a randomly selected subset of 100 (out of 7500) HD cells in the simulated ensemble are depicted here. Boxplot on the right depicts the empirical distribution of ATIs of all HD cells in the presynaptic population, with the mean ATI chosen as 50 ms in accordance with Taube, 2010. (B) Sample traces spanning ten seconds of simulation time. Green: firing rate of postsynaptic LR cell receiving HD input via depressing synapses; Orange: firing rate of postsynaptic LR cell receiving HD input via non-depressing synapses; Gray: head turning speed. Note that for the non-depressing inputs, the firing rate of the LR cell remains approximately constant throughout. In contrast, depressing synapses produce a firing rate whose fluctuations visibly reflect fluctuations in the angular head speed. (C) Firing rate of postsynaptic LR cell plotted against head speed 24 ms in the past shows a clear monotonic relationship. A similar relationship exists between firing rate and current angular head speed (see panels E, F). The inset shows firing rate as a function of angular head Figure 8 continued on next page 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 velocity. The error bars indicate the standard deviation of all firing rates in each AHV bin; note that these variabilities are low, emphasizing the consistency and reliability of this code. (D) The gray histogram depicts the correlations between input HD cell firing rate and angular head speed during active (>10 Hz firing) windows. Since presynaptic firing rate was only explicitly modulated by head direction and not head speed, these values are low, with a mean value of 0.08. The green bar depicts correlation between postsynaptic firing rate and head speed at the optimal time lag (24 ms, see panels E, F). Hence, LR cells can utilize depressing inputs from HD cells to compute de novo head speed. (E) Cross-correlation between LR cell firing rate and head speed. We computed the Pearson correlation between the firing rate of the postsynaptic LR cell and the head speed L ms in the future for varying values of L, ranging from À70 to 70 ms. Postsynaptic firing rate was maximally correlated with head speed 24 ms in the past. (F) Same as E but for cross-mutual information. Cross-mutual information identifies the same optimal lag as cross-correlation. The online version of this article includes the following figure supplement(s) for figure 8:   . Encoding of future head direction in thalamus helps to better encode present head speed in retrosplenial LR cells. (A) The schematic depicts tuning curves for an HD cell with a preferred direction of 180 degrees. Top: the tuning curves of this cell if it displayed no anticipatory firing (ATI = 0 ms). Note that clockwise and counterclockwise turns produce identical tuning curves in this case. Bottom: the tuning curves of this cell if it had an ATI = 50 ms. Note that now, during head turns in either direction, the cell will fire 50 ms prior to when the animal faces 180 degrees. Our convention takes positive angular head velocity to denote counterclockwise turning. (B) Anticipatory firing of presynaptic HD cells improves speed coding in the postsynaptic LR cell. To quantify strength of speed coding independent of latency between head speed and postsynaptic firing rate, we used the maximum correlation between head speed and postsynaptic firing rate across all time lags. Inset shows a similar relationship, now for the maximum of cross-mutual information. (C) Anticipatory firing of presynaptic HD cells improves the latency between current head speed and postsynaptic firing rate, enabling more temporally precise speed coding. The online version of this article includes the following figure supplement(s) for figure 9:  (ii) Response properties of the three LR cells whose preferred direction distributions are depicted in Ai,Bi,Ci. As the preferred direction distribution becomes less uniform (from A to C), directional encoding improves (left graph) while speed coding becomes less precise (right graph). (D) For each simulated LR cell, we quantified the heterogeneity of its input HD distribution using the standard metric for this measure (Wasserstein distance from a uniform distribution; WD). Plot shows the cumulative distribution of WD values over the entire simulated LR population. Markers A,B,C correspond to the three example LR cells shown above. (E) Scatterplot showing each LR cell's correlation with angular head speed versus its correlation with head direction. Both individually and as a population, LR cells can show conjunctive encoding of head direction and speed, with a clear tradeoff between the two. The color-coded WD is also shown for each cell. (F) Same as E but for mutual information.
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). 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 lowrheobase (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 projectionspecific 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.

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), Camk2a Cre (Jackson Laboratories, 005359), Grp-KH288 Cre (RRID:MMRRC 037585-UCD), Pvalb Cre (PV-IRES-Cre; Jackson Laboratories, 008069), and Pvalb Cre 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
Surgical 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 Grp Cre 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 mL total virus volume at each depth were given via a picospritzer at 0.05-0.07 mL/ 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
Slices 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 mm 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
During 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 mm and resistances between 2 and 5 MW. 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
Channelrhodopsin-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 (L5 sup : defined as 100 mm from the L3/5 boundary), and deep layer 5 (L5 dp ; defined as 200 mm from the L3/5 boundary). For LED-based stimulation over L1a, the objective was centered approximately 30 mm 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 L5 sup , the objective was centered 100 mm into L5, and for L5 dp , the objective was centered 200 mm 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 shortterm 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
To 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
Z-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
Representative 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 mm at 0.5 mm 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
Multiple 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 (R in ), input capacitance (C in ), and membrane time constant (t 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 ISI last /ISI first . 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 (R in ), input capacitance (C in) , and membrane time constant (t m ) were all measured from a protocol that delivered a series of small negative current steps resulting in~4 mV deflection in membrane potential. R in was calculated using Ohm's law, mean voltage change divided by mean current. t 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, C in was calculated using the formula t m = R in x C in . 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
To 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
All 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 N thalamic HD cells, via synapses exhibiting short-term depression.

Modeling the presynaptic HD population
Throughout, let T t ð Þ 2 0; 2p ½ denote the head direction at time t ms. Each HD cell has a preferred angle (PA) at which its firing rate is maximized. The firing rate at time t of an HD cell with PA takes the form: where f is a fixed function describing the shape of the cell's tuning curve. In all simulations, we utilize a Gaussian tuning curve: Here, f bg is the background firing rate, f max 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 r ; t ð Þ at time t. 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 a; b; loc; scale ð Þ is given by where G 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 r ; t ð Þ at time t. We additionally impose a post-spike refractory period of 4ms.

Anticipatory firing
Thalamic 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 À !A. Equivalently, we can think of HD cell PAs as static and instead say anticipatory firing effectively transforms the heading trajectory itself: if T t ð Þ describes the true angular trajectory, the effective trajectory experienced by our chosen HD cell is:

In vivo head direction recordings
To simulate the angular trajectory T t ð Þ, 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 In 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 x; u; g ð Þ, specified as: . x t ð Þ: fraction of vesicles remaining in the pool at time t.
. u t ð Þ: fraction of pool released into cleft upon presynaptic spike arrival at time t (i.e. release probability).
. g t ð Þ: synaptic conductance at time t.
Here, t is measured in milliseconds, g in umho, and the variables x and u are dimensionless. Depression results from depletion of x due to a spike train, followed by slow recovery of the vesicle pool; facilitation results from increase in u due to a spike train, followed by slow decay of this release probability. Formally, these dynamical variables satisfy the system: Here, t sp denotes an arbitrary presynaptic spike time, and the d term specifies a discontinuous modification to the variable upon arrival of a presynaptic spike. We use the notation f AE t ð Þ ¼ lim h!0 AE f t þ h ð Þ, and thus u þ and u À satisfy the relation u þ ¼ u À þ U 1 À u À ð Þ. The parameters t d , t f , and t g set the timescale of depression recovery, facilitation decay, and synaptic conductance decay, respectively. Note that u þ ! U always, so U represents the minimum fraction of x contributing to the conductance on each presynaptic spike arrival. Finally, the parameter A is just a tunable synaptic weight (with units of umho).
In the regime of t d ) t f , depression dominates the synaptic dynamics; conversely, in the regime of t f ) t d , facilitation dominates. Accordingly, the limits t f ! 0 and t d ! 0 describe depressiononly and facilitation-only dynamics, respectively. On the other hand, when t f and t d are on the same order of magnitude, the synapse may exhibit combined depressing-facilitating dynamics.

Synapse implementation and parameter choices
We 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: t d ¼ 270ms; t f ¼ 40ms; U ¼ 0:28 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
We model the postsynaptic LR cell using the morphologically detailed model originally presented in Brennan et al., 2020. This model, implemented in NEURON (Carnevale andHines, 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
Each postsynaptic cell receives input from N 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 0; 2p ½ . 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 2p periodicity of the density). We then used inverse transform sampling to draw from the generated distributions.

Data analysis
To quantify the extent to which postsynaptic activity reflects head turning speed, we utilized two similarity measures M: Pearson correlation and mutual information. In fact, we used the time-lagged counterparts of these measures: for a pair of time series x t , y t , their l-lagged similarity measure is M x t ; y tÀl ð Þ. Varying l, we obtain the cross-correlation and cross-mutual information between the two variables. Throughout the paper, when referring to an instance of measure M between time series x t , y t , we either report the function l7 !M x t ; y tÀl ð Þ or just the maximum value of the function: l max M x t ; y tÀl ð Þ 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.
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.  Tables. contains Table 1, Table 2, and Table 3. Table 1 presents mean ± SEM values of intrinsic properties of LR, L3 RS, and L5 RS cells. Table 2 presents the exact p-values evaluating differences between the three cell types for each intrinsic property examined. Source data files for these intrinsic properties are found in Figure 1-source data 1. Table 3 presents exact parameter distributions used for the model (see 'Modeling the Presynaptic HD Population').

Decision letter and Author response
. Transparent reporting form