Case framing at a high level
A list of epilepsy-associated genes can be strongly brain expressed without revealing where their effects converge. Neuronal genes are often detected across multiple cortical populations, and disease annotations frequently share broad terms. The initial analysis therefore produced a biologically plausible but experimentally vague conclusion: the module was neuronal.
We reframed the decision around developmental timing. The question became whether a reproducible subset of genes aligned with a transition that could be modeled—progenitor proliferation, neuronal specification, migration, excitatory-neuron maturation, inhibitory-neuron development, or synaptic network assembly.
Phenotype decomposition
Seizure type and onset were separated from global developmental delay, intellectual disability, regression, tone abnormalities, movement findings, sleep disturbance, behavioral features, and structural brain abnormalities. Exact and ancestor-expanded representations were analyzed in parallel. Leave-one-term-out analysis identified which clinical observations carried the module definition.
This reduced circularity. Genes were not considered convergent merely because all were already labeled with a broad epilepsy term. The analysis required the phenotype-defined module to show additional developmental structure in independent expression data.
Trajectory-aware cellular analysis
Developmental cortical atlases were mapped to conservative lineage stages while preserving original annotations. Module activity was summarized per donor and state. The analysis compared broad lineages for transfer and fine states for mechanism specificity. Pseudotime was treated as an inferred ordering, not literal chronological age.
Matched null sets controlled for gene length, baseline expression, annotation degree, and generic neuronal membership. Cell-cycle, immediate-early, ribosomal, and general synaptic programs were scored alongside the target module. A developmental window was retained only if the target carried information beyond these controls.
Cross-atlas validation
The primary result was repeated after excluding each atlas, donor group, and phenotype branch. We tracked direction, peak trajectory position, leading-edge genes, and neighborhood overlap. Exact cluster names differed across references, so replication was defined by compatible developmental ordering and cell lineage rather than lexical identity.
A subset of candidate genes shifted out of the stable leading edge when one broad developmental term was removed or when expression matching was applied. Those genes remained relevant to epilepsy biology but did not support the narrower shared-window conclusion.
Mechanistic interpretation
The stable module was most coherent after excitatory-neuron identity emerged and during maturation of synaptic, ion-channel, and activity-regulated programs. This does not imply one molecular pathway. It indicates that distinct upstream perturbations may become phenotypically legible during a shared cellular transition.
The analysis did not support a universal neuronal vulnerability claim. Inhibitory-neuron, progenitor, and glial programs remained relevant for subsets of genes, but they did not carry the same integrated module after resampling and phenotype refinement. The correct conclusion is modular convergence with a defined boundary.
Experimental decision
The chosen handoff is a staged perturbation design in which representative genes are edited before and after the predicted transition. Readouts include progression through state markers, neurite and synapse morphology, intrinsic excitability, network activity, and recovery under gene correction or mechanism-specific rescue.
The shared-window model predicts stronger and more coherent effects near the nominated maturation stage than in an undifferentiated progenitor state. Failure of that timing prediction would require revision even if a general cellular phenotype remained.
What the analysis establishes
Brain expression was not the finding
Broad neuronal detection confirmed tissue relevance but did not identify a mechanistic window.
A maturation-stage module survived controls
A bounded leading edge remained coherent during excitatory-neuron differentiation and maturation across references.
The convergence was partial
Progenitor, inhibitory-neuron, and glial contexts remained plausible for subsets and were not forced into one model.
Timing became a falsifiable variable
The analysis generated a prediction about when perturbation should have its strongest state-level effect.
What we conclude
We conclude that the most useful shared context for this phenotype-linked module is a defined excitatory-neuron maturation window, not the brain or neurons in general. The conclusion is supported by trajectory position, leading-edge stability, and cross-atlas direction after matched controls.
A valid follow-up must test timing. Perturbations introduced at multiple developmental stages, with rescue and orthogonal functional readouts, can determine whether the shared window represents mechanism or merely correlated expression. Gene-specific deviations should be preserved rather than interpreted as experimental failure.
Limitations
- Cortical atlases sample different regions, gestational stages, and protocols.
- Inferred trajectories are not direct lineage tracing and do not represent exact developmental time.
- Epilepsy phenotypes and gene annotations are heterogeneous and unevenly curated.
- Organoid and induced-neuron models may incompletely reproduce in-vivo maturation and circuit context.