Study S10 · generated 2026-10-13T21:32:40+00:00 · 11.3 s

Susceptible vs resilient: neural correlates of the CSDS phenotype

Which per-animal network, activity and coding measures separate the lab-classified susceptible and resilient mice?

Question

Do susceptible and resilient mice differ in how their neuronal activity, functional networks and population codes respond to defeat, and in how those responses evolve across days?

Background & rationale

Only a subset of mice subjected to CSDS develop lasting social avoidance (susceptible); the rest remain resilient (Krishnan et al. 2007; Golden et al. 2011). Resilience is an active process with its own neural signatures. Pre-stress nucleus accumbens D1-MSN activity predicts later resilience (Muir et al. 2018), a brain-wide electrical network signature tracks vulnerability (Hultman et al. 2018), and susceptible and resilient mice differ in large-scale functional connectivity after CSDS (Mol Psychiatry 2025). Hippocampal social coding also differs between the phenotypes (Biol Psychiatry Glob Open Sci 2025). Here the lab classified 3 mice as resilient (4703, 6597, 6669) and 3 as susceptible (6445, 6672, 6675), so we can ask which single-neuron network measures, measured during the defeat period, differ between the groups.

Methodology

Per-animal features (days 1–9 without day 7). * Behavior (S01): defeat-response index (fraction of SD time submissive + freezing + escape) and the annotation-based SI proxy. * Activity and networks (S02, S04): for event rate, mean pairwise correlation, significant-edge density, clustering, global efficiency and modularity, the animal's baseline (BL) mean and the defeat response (SD − BL) and sensory-stress response (SS − BL), computed over days on which the animal has all of BL, SD and SS. * Reorganisation (S06): FC context specificity and the drift slope of SD networks (similarity per day of separation, from the pairwise session similarities). * Coding (S08): context-decoding accuracy and behavior-decoding accuracy above each animal's own null.

Group comparison. For each feature, the susceptible − resilient mean difference was tested with an exact permutation test over all 20 ways of labelling 3 of the 6 animals as susceptible (two-sided). With 3 vs 3 the smallest attainable p is 0.10, which happens only when the groups are completely separated. The study therefore reports Hedges' g and complete separation (every susceptible value above, or below, every resilient value) as the primary descriptors. Under the null, any given feature shows complete separation with probability 0.10. BH-FDR q-values are reported for completeness.

Trajectories. Session-level measures were averaged per phenotype for each recording day and condition to show whether the groups diverge as defeats accumulate.

Continuous coupling. Spearman ρ (exact, 720 permutations) between each neural feature and the defeat-response index.

Parameters

ParameterValue
groupslab classification (3 vs 3)
testexact permutation (20)
days1, 2, 3, 5, 9 (day 7 excluded)

Highlights

info

Phenotype groups

3 susceptible · 3 resilient

Lab classification: resilient 4703, 6597, 6669; susceptible 6445, 6672, 6675.

null

Features separating the groups

2/25

Complete separation of susceptible and resilient animals; ≈2.5 expected by chance (binomial p = 0.73, ignoring correlations among features).

info

event rate (Hz) · SS − BL

g = -3.22

Susceptible lower (-0.107 vs 0.0151), complete separation, exact p = 0.100 (minimum attainable).

info

event rate (Hz) · SD − BL

g = -1.90

Susceptible lower (-0.31 vs -0.0536), complete separation, exact p = 0.100 (minimum attainable).

caution

Power

n = 3 per group

No 3-vs-3 comparison can reach p < 0.05; these are candidate markers.

Results

Group separation. Of 25 per-animal features, 2 completely separated the 3 susceptible from the 3 resilient mice, against ≈2.5 expected by chance: event rate (Hz) · SS − BL (g = -3.22, susceptible < resilient); event rate (Hz) · SD − BL (g = -1.90, susceptible < resilient). The chance of seeing at least this many separating features is p = 0.73 (binomial, treating features as independent; they are correlated, so this is optimistic). The count is compatible with chance, so no single feature should be read as a phenotype marker. Each such feature has the minimum attainable exact p of 0.10. Table 1 lists all features ranked by separation and |Hedges' g|, and Figure 1 shows the strongest ones per animal.

Largest effects. event rate (Hz) · SS − BL: susceptible -0.107 vs resilient 0.0151 (g = -3.22); event rate (Hz) · SD − BL: susceptible -0.31 vs resilient -0.0536 (g = -1.90); SD network drift (Δsim/day): susceptible -0.0125 vs resilient -0.0372 (g = +1.35); SI proxy (annotations): susceptible 1.41 vs resilient 3.13 (g = -0.83); modularity Q · SD − BL: susceptible 0.0464 vs resilient 0.0149 (g = +0.80).

Trajectories. Figure 2 shows activity and network measures per recording day for each phenotype in BL, SD and SS sessions. Figure 3 shows the BL → SD → SS profile of each group.

Interpretation. With three animals per group, this analysis cannot establish a phenotype difference. Features that separate the groups and agree with prior literature are the natural pre-registered hypotheses for a larger cohort.

Figures

Features that best distinguish susceptible from resilient mice
Figure 1. Features that best distinguish susceptible from resilient mice. Per-animal values (labelled by ID) for the 8 features ranked highest by complete separation and |Hedges' g|. Bars = group means. p = exact permutation (minimum attainable 0.10).
Day-by-day trajectories by phenotype
Figure 2. Day-by-day trajectories by phenotype. Rows: BL, SD and SS sessions; columns: activity and network measures. Thin lines = animals, thick = phenotype mean ± SEM. Day 7 excluded.
Condition profiles by phenotype
Figure 3. Condition profiles by phenotype. Mean over days of each measure in BL, SD and SS for each animal (thin) and phenotype (thick).

Tables

Table 1. Susceptible vs resilient: per-animal features

Exact permutation over all 20 relabelings (two-sided; minimum attainable p = 0.10). Hedges' g = susceptible − resilient in pooled-SD units. Download CSV

featuresusceptible_meanresilient_meandifferencehedges_gcomplete_separationp_exactn_relabelingsrho_defeat_indexp_rhoq_fdr
event rate (Hz) · SS − BL-0.1070.0151-0.122-3.22True0.120-0.7710.1031
event rate (Hz) · SD − BL-0.31-0.0536-0.256-1.9True0.120-0.7710.1031
SD network drift (Δsim/day)-0.0125-0.03720.02461.35False0.220-0.1430.8031
SI proxy (annotations)1.413.13-1.73-0.825False0.4100.20.7831
modularity Q · SD − BL0.04640.01490.03150.801False0.3200.4860.3561
defeat-response index0.310.1990.1120.76False0.4201
modularity Q · SS − BL-0.02080.0343-0.0551-0.701False0.320-0.08570.9191
clustering · SS − BL0.0007520.03-0.0292-0.628False0.5200.3140.5641
FC context specificity0.001040.0091-0.00806-0.615False0.420-0.20.7141
event rate (Hz) · BL1.591.92-0.33-0.475False0.6200.5430.2971
decoding behavior SD (above null)0.06440.05550.008820.371False0.620-0.2570.6581
decoding behavior BL+SS (above null)0.1020.0980.003970.287False0.8200.3710.4971
global efficiency · SS − BL0.00219-0.006070.008260.287False0.7200.4860.3561
global efficiency · SD − BL0.01760.01070.006910.172False0.6200.4860.3561
clustering · SD − BL0.01610.01110.004950.154False0.9200.60.2421
mean pairwise r · SD − BL0.002810.00906-0.00625-0.149False120-0.4860.3561
significant-edge density · SS − BL-0.0507-0.06350.01280.112False120-0.08570.9191
global efficiency · BL0.2560.262-0.00502-0.111False0.820-0.4290.4191
clustering · BL0.480.4760.004360.11False0.920-0.08570.9191
mean pairwise r · BL0.4120.421-0.00868-0.0903False0.920-0.4290.4191
modularity Q · BL0.3560.361-0.00467-0.0546False0.920-0.6570.1751
decoding context BL/SD/SS (above null)0.1290.134-0.00509-0.0382False120-0.7140.1361
mean pairwise r · SS − BL-0.0298-0.028-0.00181-0.0235False0.920-0.08570.9191
significant-edge density · SD − BL0.004190.00482-0.000626-0.0107False1200.08570.9191
significant-edge density · BL0.7470.748-0.00131-0.00939False120-0.2570.6581

Table 2. Per-animal feature matrix

All features used in the comparison. Download CSV

subjectevent rate (Hz) · BLmean pairwise r · BLsignificant-edge density · BLclustering · BLglobal efficiency · BLmodularity Q · BLevent rate (Hz) · SD − BLmean pairwise r · SD − BLsignificant-edge density · SD − BLclustering · SD − BLglobal efficiency · SD − BLmodularity Q · SD − BLevent rate (Hz) · SS − BLmean pairwise r · SS − BLsignificant-edge density · SS − BLclustering · SS − BLglobal efficiency · SS − BLmodularity Q · SS − BLdefeat-response indexSI proxy (annotations)FC context specificitySD network drift (Δsim/day)decoding behavior BL+SS (above null)decoding behavior SD (above null)decoding context BL/SD/SS (above null)phenotype
47031.720.4790.7810.4490.2330.3470.06730.0499-0.003440.00872-0.01380.02440.03620.00687-0.01350.03740.007150.08160.0842.490.0212-0.02180.08770.04650.256resilient
64450.8390.4550.8190.5130.2990.453-0.2470.01780.0259-0.01390.005130.00495-0.0623-0.012-0.0384-0.0358-0.0152-0.06610.1780.454-0.00707-0.002680.09840.08570.175susceptible
65971.840.3960.7830.480.2990.323-0.201-0.03240.02840.03160.04730.02050.00304-0.00695-0.0247-0.0106-0.0339-0.05760.2151.960.0108-0.03420.1040.06420.156resilient
66692.210.3880.6810.4970.2530.413-0.02740.00974-0.0105-0.00702-0.00146-0.0002020.00596-0.0838-0.1520.06310.00850.07880.2974.96-0.0047-0.05560.1020.0559-0.00955resilient
66722.330.3020.5780.4890.2260.32-0.3820.01510.05330.04740.0540.0897-0.1330.03190.0475-4.4e-050.0269-0.00210.4330.00619-0.009720.09140.03630.0494susceptible
66751.620.480.8440.4390.2440.295-0.299-0.0245-0.06670.0148-0.00630.0444-0.125-0.109-0.1610.0381-0.005160.005910.322.360.00398-0.02520.1160.07110.162susceptible

Limitations

  • n = 3 per group. No feature can reach p < 0.05, and with ~30 features about 3 would be expected to separate the groups by chance. Results identify candidate markers for a confirmatory cohort and are not evidence of a phenotype difference.
  • Mouse 6597 has no sensory-stress block on day 9, so that animal's network features use days 1, 2, 3 and 5. Mouse 6669 lacks the SIT imaging block; day-1–9 features for that animal are complete. Day 7 is excluded for every mouse.
  • No non-defeated control group, so "resilient" differences cannot be separated from "less affected by the assay".

References

  1. Krishnan V. et al. (2007) Molecular adaptations underlying susceptibility and resistance to social defeat in brain reward regions. Cell 131:391–404.
  2. Golden S.A. et al. (2011) A standardized protocol for repeated social defeat stress in mice. Nat Protoc 6:1183–1191.
  3. Muir J. et al. (2018) In vivo fiber photometry reveals signature of future stress susceptibility in nucleus accumbens. Neuropsychopharmacology 43:255–263.
  4. Hultman R. et al. (2018) Brain-wide electrical spatiotemporal dynamics encode depression vulnerability. Cell 173:166–180.
  5. Resting-state fMRI reveals altered functional connectivity associated with resilience and susceptibility to CSDS (2025) Mol Psychiatry. doi:10.1038/s41380-025-02897-2
  6. Diminished social memory and hippocampal correlates of social interactions in chronic social defeat stress susceptibility (2025) Biol Psychiatry Glob Open Sci. https://www.bpsgos.org/article/S2667-1743(25)00009-6/fulltext
  7. Button K.S. et al. (2013) Power failure: why small sample size undermines the reliability of neuroscience. Nat Rev Neurosci 14:365–376.