Adult-neurogenesis allows for representational stability and flexibility in early olfactory system
Figures
Spiking network model of adult-neurogenesis.
(A1). Schematic of early olfactory circuit. M/T: mitral/tufted cells, GC: granule cells, abGCs: adult-born GCs, PCs: piriform cortical cell, FFIN: feedforward inhibitory neuron, FBIN: feedback inhibitory neuron. Example plus signs indicate excitatory synapses and minus signs indicate inhibitory synapses. (A2). Schematic of synaptic reshuffling of GCs due to adult-neurogenesis. On each day, 10% of GCs have their weights reshuffled, including the synaptic weights from M/T cells, other GCs or short-axon cells (SACs), and PCs to GC (feedback). (B1). Partial weight matrix on two example days between GCs and M/T cells (top), other GCs or SACs (middle), and feedback from PCs (bottom). (B2). Histograms of all synaptic weights across days. (B3). Weight matrix dissimilarity between Day-0 and each other day. (C1). An example model odor generated by stimulating different combinations of glomerular identity and timing of activation. (C2). Pair-wise correlations between model odors show the relative similarities and differences across all model odors.
Odor definition and synaptic weight changes.
(A) Glomerular activation latency for all 100 odors we defined. Each column corresponds to one odor. (B) Histogram of the number of activated glomeruli by all odors (n = 100 odors). (C) Histogram of correlation coefficients between each pair of odors, measuring the similarity between odors. (D) Sparsely sampled synaptic weights on Day-0 (x axis) and Day-10 (y axis). Each dot is a synapse. Left: M/T to GCs; middle: GC/SAC to GCs; right: PCs to GCs.
Adult-neurogenesis modulates individual M/T cell responses but preserves population representation.
(A1) Trail-averaged firing rate of two example M/T cells responding to the same odor (within-odor) on three different days (mean ± SD, n = 10 trials). (A2) Firing rate patterns of odor-activated M/T cells to the same odor on 3 days. M/T cells driven by different glomeruli are separated by a white dashed line. The two arrows correspond to the two example cells in (A1). (A3) Pairwise within-odor full-ensemble correlation of single-trial M/T cell responses between each day averaged across all odors (n = 100). Each small box separated by the dashed lines is a 10 × 10 matrix (n = 10 trials) corresponding to autocorrelation (same day on diagonal) and cross-correlation (different days off diagonal). (A4) Full-ensemble correlation of trial-averaged M/T cell responses for within-odor (black solid line with error bar: mean ± SD, n = 100 odors) and across-odor (black dashed line with error bar: mean ± SD, n = 10 pairs) between Day-0 and each other day. Gray dashed line with no error bar: same as the black solid line in (B4). (B1) Low-dimensional trajectories of M/T responses to two example odors (color) on two different days (darkness of color). Thin curves: single-trial trajectories, thick curves: trial-averaged trajectories, points: maximal-distance points on trajectory to the origin. (B2) Only maximal-distance points are shown for different odors on different days. (B3) Similar to (A3) but for reduced-space correlation computed using the low-dimensional M/T trajectories (single-trial). (B4) Similar to (A4) but for reduced-space correlation computed using the low-dimensional M/T trajectories (trial-averaged). Gray dashed line with no error bar: same as the black solid line in (A4).
Example M/T cell responses on Day-0 and Day-10.
Three example M/T cells (column) in response to the five example odors (row). Black: response on Day-0; red: response on Day-10. Some cells increase their responses while some cells decrease their responses from Day-0 to Day-10 (mean ± SD, n = 10 trials).
Adult-neurogenesis modulates both individual PCs’ responses and population representations.
(A1). Trail-averaged firing rate of two example PCs responding to the same odor on three different days (mean ± SD, n = 10 trials). (A2). Firing rate patterns of partial PCs responding to the same odor on 3 days. The two arrows correspond to the two example cells in (A1). (A3) Pairwise within-odor full-ensemble correlation of single-trial PCs’ responses between each day averaged across all odors (n = 100). Each small box separated by the dashed lines is a 10 × 10 matrix (n = 10 trials) corresponding to autocorrelation (same day on diagonal) and cross-correlation (different days off diagonal). (A4). Full-ensemble correlation of trial-averaged cell (PCs: purple; M/T: black) responses for within-odor (solid line with error bar: mean ± SD, n = 100 odors) and across-odor (dashed line with error bar: mean ± SD, n = 10 pairs) between Day-0 and each other day. The curves for M/T cells are the same lines as in Figure 2. (B1) Low-dimensional trajectories of PCs’ responses to two example odors (color) on two different days (darkness of color). Thin curves: single-trial trajectories, thick curves: trial-averaged trajectories, points: maximal-distance points on trajectory to the origin. (B2) Only maximal-distance points are shown for different odors on different days. (B3) Similar to (A3) but for reduced-space correlation computed using the single-trial principal component analysis (PCA) trajectories of PCs. (B4). Similar to (A4) but for reduced-space correlation computed using the trial-averaged PCA trajectories of PCs.
Example piriform cell responses on Day-0 and Day-10.
Three example piriform cells (column) in response to the five example odors (row). Black: response on Day-0; red: response on Day-10 (mean ± SD, n = 10 trials).
Example piriform cell responses (trial-averaged) on all 10 days.
(A) Another three example piriform cells in response to four odors from Day-0 and Day-10. Each curve is the response averaged across 10 trials (smoothed by a Gaussian kernel with 20-ms width). Note some piriform cells increase or decrease their responses, while some cells change their temporal profile (cell-3).
Cell response changes between Day-0 and Day-10.
(A) Histogram showing the differences between the trial-averaged response on Day-10 and the trial-averaged response on Day-0 for each odor-activated M/T cell (n = 100 odors). Positive values indicate that the responses on Day-10 are larger than Day-0. (B) Same as (A) but for piriform cells.
Odor manifold and geometric reshaping of odor representations in PCx.
Odor manifolds of M/T cells (top) and PCs (bottom) on two example days (red and blue) are plotted by clusters of transparent round dots in the principal component analysis (PCA) space. Each column contains the trial-averaged trajectories on two days (red and blue) evoked by three example odors. M/T trajectories stay close while PCs trajectories are reshaped in different ways: shifted (A), rotated (B), and warped (C).
Quantifying geometric reshaping and decoding analysis on M/T and PCs.
(A1) Cosine similarity to quantify the geometric reshaping using population firing rate as a function of intervals (number of days separated). Mean ± SD, n = 100 odors. (A2). Decoding accuracy by K-nearest neighbors algorithm trained on M/T population firing rate of each sample day and tested on M/T population firing rate of each reference day. (A3) Same as (A2) but using PCs population firing rate. (B). Same as (A) but using principal component analysis (PCA) trajectories. Mean ± SD, n = 10 odor pairs.
Experience-dependent plasticity enhances representational stability in PCx.
(A1) Different types of plastic synapses (blue) related to an abGC. Circle: inhibitory; triangle: excitatory. (A2) Dependence of synaptic modification on pre/post inter-spike interval used for the spike-timing-dependent plasticity (STDP). S1 and S2 are two example synapses with different modification – S1: facilitated; S2: suppressed. (B) Weight history of the two example synapses in A2 when an odor is repeatedly presented. Adult-neurogenesis happens on Day-4 which randomly resets the weights. The synaptic weights stay constant without plasticity (black) while change trial by trial with plasticity (blue). (C1) Left: reduced-space correlation of trial-averaged M/T trajectories between Day-0 and each day. Blue solid line with error bar: within-odor and with STDP (mean ± SD, n = 10 odors); black dashed line without error bar: within-odor and without STDP. Black dashed line with error bar: across-odors, mean ± SD, n = 10 pairs. Right: drift rate for M/T trajectories with and without STDP. (C2) Same as C1 but for PCs trajectories.
Synaptic weight history of spike-timing-dependent plasticity (STDP).
(A) Traces of five example synapses from M/T cells to GCs. STDP was applied when adult-neurogenesis happened at the start of the day. STDP drove the synaptic weight (y axis) to certain values through trial-by-trial exposure to certain odors. For some synapses, adult-neurogenesis only happened on 1 day while for others it happened on multiple days. (B) With STDP applied, three example piriform cells (column) in response to the three example odors (row). Black: response on Day-0; red: response on Day-10 (mean ± SD, n = 10 trials). (C) Slope of the synaptic weight traces over the last five trials of odor exposure with STDP applied. The slopes were all close to zero indicating the convergence of synaptic weights. Boxplot with the central mark indicating the median, and the bottom and top edges indicating the 25th and 75th percentiles, respectively. The whiskers extend to the most extreme data points not considered outliers. n = 100 randomly sampled synapses for each day.
Tables
Parameters of Izhikevich neuron model for different cell types.
| M/T cells | GCs | PCs | FFIs/FBIs | |
|---|---|---|---|---|
| 0.2 | 0.2 | 0.2 | 0.2 | |
| −65 | −65 | |||
| 2 | 2 |
Network parameters controlling the connectivity between cell types.
| Connection density | Average synaptic strength | References | |
|---|---|---|---|
| MCtoMC (same glomerulus) | 0.8 | 0.25 | Urban and Sakmann, 2002; Egger and Urban, 2006; Schoppa and Westbrook, 2002; Soucy et al., 2009 |
| MCtoMC (different glomeruli)* | 0 | 0 | |
| MC2GC | 0.3 | 0.25 | Schoppa and Urban, 2003; Kapoor and Urban, 2006; Burton and Urban, 2015; Willhite et al., 2006 |
| GC2MC | 0.02 | –0.4 | |
| GC2GC | 0.05 | –0.1 | Boyd et al., 2012; Burton and Urban, 2015 |
| MC2PC | 0.5 | 0.06 | Sosulski et al., 2011; Miyamichi et al., 2011; Stern et al., 2018; Luna and Schoppa, 2008 |
| MC2FFI | 0.2 | 0.2 | |
| PC2GC | 0.9 | 0.03 | Otazu et al., 2015; Nunez-Parra et al., 2013; Padmanabhan et al., 2018 |
| PC2PC | 0.01 | 0.1 | Stern et al., 2018; Franks et al., 2011; Suzuki and Bekkers, 2012; Large et al., 2016 |
| FFI2PC | 0.1 | –0.1 | |
| FBI2PC | 0.8 | –0.1 | |
| FFI2FFI | 0.01 | –1.0 | |
| PC2FBI | 0.02 | 0.3 | |
| FBI2FBI | 0.02 | –0.5 |
-
*
We assumed zero connection between M/T cells associated with different glomeruli.