Emergent Topological Brain Organoids

---
name: emergent-topological-brain-organoids
description: "Apply persistent homology to MEA recordings of spontaneous brain-organoid activity. Detects H1 loops and H2 voids in correlation networks, compares them to rate-preserving null models, and identifies non-redundant loop-carrying cores."
metadata:
  arxiv_id: "2607.16517"
  published: "2026-07-21"
  authors: "Eve Bodnia, Margaux Basart, Sofie Hai, Lenzie Ford, Nina Miolane, Kenneth S. Kosik, Dirk Bouwmeester, Lincoln D. Carr"
  tags: [neuroscience, brain-organoid, persistent-homology, topological-data-analysis, MEA, neural-dynamics, correlation-network]
license: Complete terms in LICENSE.txt
---

# Emergent Topological Structure in Brain-Organoid Activity

Apply persistent homology (Vietoris–Rips filtration) to microelectrode-array (MEA) recordings of spontaneous activity from human and mouse cortical organoids, to detect loops (H1) and voids (H2) that exceed a rate- and population-preserving null model.

## When to Use

- Analyzing spontaneous multi-unit spike recordings from organoids, slices, or dense in-vivo arrays.
- You want to test whether correlation networks carry structured topology beyond firing rate and population bursting.
- You need to identify the non-redundant core of units that carry topological features.
- You want to know how many nodes are required for persistent homology to resolve H1 and H2.

## Core Methodology

### 1. Build the Correlation Network

- Spike-sort MEA data; use sorted single units as network nodes.
- Gaussian-smooth spike trains with a 50 ms kernel.
- For each pair of units (a, b), compute the maximum normalized cross-correlation over a short lag window (0, 10, or 20 ms):

  C(a, b) = max_j  ( Σ_i a_i b_{i+j} ) / ( sqrt(Σ_i a_i²) * sqrt(Σ_i b²) )

- Convert to a dissimilarity: d_ij = 1 - C_ij.

### 2. Vietoris–Rips Filtration and Persistent Homology

- Build a Vietoris–Rips filtration on d_ij using Ripser (or similar).
- Report density-indexed Betti curves β0(ρ), β1(ρ), β2(ρ), where ρ is edge density.
- Integrate Betti curves (e.g., ∫β1(ρ)dρ) for a threshold-free summary.

### 3. Rate- and Population-Preserving Null Model

- Generate raster-marginals surrogates: preserve every unit’s total spike count and every time bin’s total population activity, while destroying higher-order co-firing.
- Use 10^5 binary swaps per surrogate; create ~100 surrogates.
- Recompute correlations and topology on each surrogate.
- Compare integrated β1 / β2 from data to the surrogate distribution; report empirical rank p-values.

### 4. Identify Non-Redundant Core

- Compute H1 cocycle representatives for the persistent loops.
- Randomly remove 10% of units and recompute topology (average over 100 trials).
- Targeted removal: remove the 10% of units most frequently appearing in H1 cocycle representatives.
- Compare retained integrated β1 and bottleneck distance between original and post-removal persistence diagrams.

### 5. Size-Scaling Analysis

- Repeat for networks with varying N (e.g., subsample units).
- H1 becomes resolvable above roughly N ≈ 100 in the organoid datasets.
- H2 emerges significantly only in larger networks (N ≥ 119 in the reported data).

## Key Findings

- H1 loops exceeded the null in 14 of 18 organoid datasets (p ≤ 0.05).
- Loops peak at low edge density (ρ ≲ 0.15) and are carried by strongly co-active units, not weak correlations.
- Random 10% removal retained 92.5% of integrated β1; targeted removal retained 72.6% and produced larger bottleneck distances (median ratio 1.48 vs random).
- H2 voids appeared in 6 datasets, all with N ≥ 119, suggesting larger recordings are needed to resolve higher-order structure robustly.
- Loop structure reflects co-firing, not electrode layout; geometry is topologically trivial in the datasets tested.

## Practical Workflow

1. Load spike-sorted MEA data; create spike trains and quality-control units.
2. Gaussian-smooth with 50 ms kernel; compute pairwise maximum cross-correlation at 0–20 ms lags.
3. Build correlation distance matrix d = 1 - C.
4. Run Ripser on d; obtain persistence diagrams and Betti curves indexed by density.
5. Generate 100 raster-marginals surrogates; compute null distribution of integrated β1 and β2.
6. Report p-values and z-scores for data vs null.
7. (Optional) Compute cocycle representatives and test robustness to random vs targeted node removal.
8. (Optional) Subsample units to estimate topology-resolvable node count.

## Parameters and Defaults

| Parameter | Default | Notes |
|-----------|---------|-------|
| Smoothing kernel | 50 ms Gaussian | Matches tens-of-ms synaptic timescale |
| Lag window | 0 ms (also 10, 20 ms) | Zero lag is symmetric for undirected VR |
| Surrogate swaps | 10^5 per surrogate | Sufficient for destroying higher-order co-firing |
| Surrogate count | 100 | For empirical null distribution |
| Removal fraction | 10% | For robustness / core analysis |
| Density range | ρ ∈ [0, 1] | Threshold-free comparison |

## Pitfalls

- **Small N**: networks below ~50 units may not resolve H1 above surrogate noise; report this as an insufficiency rather than absence of structure.
- **Smoothing too wide or narrow**: a kernel far from the synaptic timescale can suppress or distort co-firing structure.
- **Using Erdős–Rényi as null**: this only tests structure against random connectivity, not against rate/population effects. Use raster-marginals surrogates.
- **Directed vs undirected**: undirected VR is symmetric; for firing-order information, use directed flag complexes separately.
- **Loop count alone is misleading**: a hub deletion can create new loops. Use bottleneck distance to measure topological disruption.

## Activation Keywords

- brain organoid topology
- persistent homology neural recordings
- MEA correlation network loops
- H1 H2 neural topology
- rate-preserving null model
- loop-carrying core
- topological data analysis neuroscience

## References

- Giusti, Pastalkova, Curto, Itskov (2015). Clique topology reveals intrinsic geometric structure in neural correlations. *PNAS*.
- Schreiber et al. (2003). A new correlation-based measure of spike timing reliability. *Neurocomputing*.
- Okun et al. (2012). Diverse coupling of neurons to populations in sensory cortex. *Nature*.
- Ripser: https://ripser.scikit-tda.org
- Lancaster & Pasca organoid protocols (see paper references).