What constant-depth threshold circuits say about what wetware is for
A team at INESC-ID and Instituto Superior Tecnico in Lisbon has designed the first neuromorphic circuit for exhaustive epistasis detection, one of the most combinatorially brutal problems in bioinformatics. The headline result is theoretical, but the method behind it is the part organoid intelligence should steal: define the hardware's real restrictions first, then derive the algorithm, not the other way around.
Source: Constant Depth Threshold Circuits For Exhaustive Epistasis Detection, Ribeiro, Ilic and Sousa, arXiv:2605.29719, 2026. Primary source. Read: the full arXiv HTML version, including all circuit constructions, the complexity analyses and the authors' own limitations.
What the work claims
Epistasis detection asks which combinations of genetic variants, usually SNPs, jointly influence a phenotype such as disease presence. An exhaustive search over k-way interactions must enumerate on the order of n^k combinations of n SNP positions, and in practice k runs from 2 to 4. GPUs handle this poorly: the best classical approaches burn power on memory transfers and synchronization while counting frequencies into contingency tables.
Ribeiro, Ilic and Sousa claim a neuromorphic circuit that computes all combination frequencies in a fully pipelined way, with total runtime O(n^k), that is, bounded by the number of combinations to calculate with no additional complexity overhead, using log-linear space in the relevant variables. The dataset is embedded in the hardware itself as addressable firing-pattern memories, and counting is done by constant-depth population-count circuits built from threshold gates. The authors are explicit that this is, to their knowledge, the first neuromorphic algorithm for the problem, and they derive it under stated hardware restrictions rather than idealized ones.1
How it works
Three ideas carry the construction. First, memory is co-located with computation: each sample's genotype is stored as a binarized firing pattern in a synaptic stack circuit, addressable and non-destructive, so no off-chip transfer ever occurs. Binarizing the three possible genotype values (0, 1, 2) into three bits means each memory neuron spikes about one third of the time.
Second, counting is done by threshold gates, the circuit-theory abstraction behind spiking neurons, rather than by arithmetic units. The paper builds parity, population count (POPC) and binary sum circuits in the class TC0, meaning constant depth. With unbounded fan-in these run in constant time; once you admit that a real neuron has a bounded number of input synapses, the honest depth becomes O(log_Fin n) where F_in is the fan-in limit. A repeater circuit generates the SNP combinations themselves, again from LIF-style neurons.
Third, everything is pipelined: one dataset entry is processed per timestep once the pipeline fills, and the pipeline latency is only O(log n). Counting synapses, the full circuit occupies O(mn + m log^3 m) synapses and O(mn) neurons for n SNPs and m samples. A parity sub-circuit on l bits uses linear neurons and quadratic connections; the binary-sum circuit's synapse count grows roughly as l^3 n^2, which the authors report openly.1
Where a skeptic should push
The single most load-bearing assumption is that this circuit can ever be built at the scale where the asymptotics matter. There is no chip run in the paper, no simulation, no energy measurement. Every energy statement is a counting argument: neurons spike a third of the time, high-frequency combinations add at most O(log m) spikes per combination, therefore the implementation is "favorable to low energy". That may be right, but it is a hypothesis, not a result, and the authors say as much in their conclusion, where they note that comparing against classical algorithms properly requires incorporating data-movement costs they did not model.
The space bill deserves real scrutiny. O(mn) neurons and synapses means the dataset lives in the substrate; for modern biobank scales (hundreds of thousands of samples, millions of SNPs) the neuron count would be astronomical. The scheme is honest about this, since co-located memory is precisely the point, but it means the construction is an existence proof for a regime where the data fits, not a drop-in replacement for a GPU cluster. There is also a subtlety the paper handles well and the reader should notice: the glamorous constant-depth results all assume unbounded fan-in, and the authors themselves warn that ignoring such restrictions produces unrealistic conclusions, citing the fact that rational-valued currents would let one neuron encode infinite information in a Cantor-style stack. The constant-depth headline quietly becomes logarithmic depth under hardware limits. That is correct science, and it is also a reminder that the boldest number in any neuromorphic claim is usually the one with the least realistic assumption under it.
The organoid lesson in native workloads
Organoid intelligence has no equivalent of this paper, and that absence is the story. Nobody has sat down and asked, with this level of rigor, which algorithms are native to a ball of living neural tissue: spikes as events, synapses as co-located analog memory, graded membrane state, severe fan-in limits, stochastic firing, and no off-chip anything. The field instead benchmarks organoids on tasks imported from machine learning, where the substrate loses by construction, or it makes substrate-free claims about "real-number computing in neurons" that are the exact analogue of the Cantor-encoding assumption this paper explicitly rejects.
The opportunity is methodological borrowing. A complexity program for organoid computing would do for wetware what this paper does for threshold silicon: define the restrictions honestly (fan-in, precision, memory capacity, metabolic power), then derive the workload class where those restrictions are advantages. Exhaustive enumeration with co-located memory is the exemplar: tissue stores information in exactly the co-located, analog, fan-in-limited form the theory assumes. It is conceivable that frequency counting, associative enumeration, or population-code aggregation are native wetware workloads in a way that matrix multiplication is not, and no benchmark culture will tell you that, because benchmarks measure the wrong task class.
The threat cuts both ways. Without a native-workload theory, the field is defenseless against two failure modes: dismissing organoids because they underperform GPUs at GPU tasks, and overselling them with fantasy assumptions about infinite analog precision or all-to-all connectivity. There is also an obsolescence risk. Neuromorphic silicon is building exactly this kind of native-workload library first; the Ion Channel Neuromorphic Platform and similar efforts are already circling the same co-located-memory territory. If wetware cannot name its native task class in the same rigorous language, it will be compared against silicon on silicon's terms and found redundant.
The bottom line
Established: a sound theoretical construction showing that exhaustive epistasis detection maps onto threshold-gate neuromorphic primitives with runtime bounded by the enumeration itself, under explicitly stated hardware restrictions. Hypothesis: that this yields real energy advantages over GPUs at practical scales. The paper is a model of intellectual honesty about which is which. For organoid intelligence, the transferable asset is not the circuit but the discipline: restrictions first, algorithms second, claims that survive both. What would confirm the approach is a hardware or large-scale simulation mapping SynOps to measured joules against a tuned GPU baseline; what would break it is evidence that dataset-scale space requirements swamp any asymptotic win.
Frequently asked questions
What is epistasis detection?
It is the search for combinations of genetic variants, typically SNPs, that jointly influence a trait such as a disease. Exhaustive detection checks all k-way combinations, which grows as n^k for n SNP positions and quickly becomes one of the most demanding combinatorial workloads in bioinformatics.
What does constant depth actually buy?
With idealized unbounded fan-in, threshold-gate circuits can compute parity, population count and binary sums in a constant number of time steps, which is what makes the fully pipelined O(n^k) enumeration possible. Under realistic fan-in limits the depth degrades gracefully to O(log n), which the pipeline latency already reflects.
Is this circuit running on real hardware?
No. The paper is a theory and circuit-design contribution with no fabrication, simulation, or measurement. The authors state that hardware comparison, including data-movement costs across chip cores, is future work.
Why does the neuron count matter for the claim?
Space scales as O(mn) neurons and O(mn + m log^3 m) synapses for m samples and n SNPs, because the dataset is stored in the hardware. The scheme only pays off when the dataset fits in the substrate, which confines it to regimes far smaller than modern biobanks.
What should organoid computing take from this?
The discipline of deriving algorithms from stated substrate restrictions instead of porting benchmarks. A rigorous account of which computations are native to spiking, analog, fan-in-limited tissue would give the field task classes it can actually win, and protect it from both hype and premature dismissal.
References
- A. Ribeiro, A. Ilic and L. Sousa. Constant Depth Threshold Circuits For Exhaustive Epistasis Detection. arXiv:2605.29719. 2026. https://arxiv.org/abs/2605.29719. Accessed 2026-09-18.