The research, conducted by Armaan Raina and John Meitzen of North Carolina State University, focused on medium spiny neurons in the nucleus accumbens, a GABAergic output cell that serves as the key computational integrator for this hub of motivated and reward-driven behavior. Earlier work from the same laboratory had shown that the electrophysiological properties of these neurons shift across the four phases of the rodent estrous cycle, including action potential characteristics, passive membrane attributes, and the properties of miniature excitatory postsynaptic currents, the tiny synaptic events that reveal how individual glutamatergic synapses are functioning. The team hypothesized that if those hormone-driven shifts carry real information, then machine learning models trained on the data should be able to recover it, and the amount of information should differ depending on which type of electrical measurement is analyzed.

To test this idea, the researchers drew on a publicly available archive of whole-cell patch clamp recordings from adult female Sprague-Dawley rats, originally collected across diestrus, early proestrus, late proestrus, and estrus phases, with each phase assigned through careful vaginal cytology and confirmed by circulating sex steroid hormone levels. From these raw recordings they extracted three distinct feature sets. The action potential set captured metrics such as burst length, spikes per burst, and spike frequency from more than four thousand depolarization-evoked spike trains. The passive set included thirteen properties such as resting membrane potential, input resistance, membrane time constant, and inward rectification ratio, derived from over five thousand hyperpolarizing voltage responses. The miniature excitatory postsynaptic current set described more than sixty-one thousand synaptic events using measures like amplitude, rise time, decay, and baseline current.

Seven classical machine learning architectures were then trained on each feature set: random forest, gradient boosting, logistic regression, k-nearest neighbors, support vector classifier, multi-layer perceptron, and a single decision tree, producing twenty-one models in total. Crucially, the researchers deliberately avoided hyperparameter tuning, since their goal was not to squeeze out maximum accuracy but to compare how much phase-related information each data type inherently encodes. Each model was trained five times on fixed random seeds, with features Z-score normalized and the data split seventy-five to twenty-five between training and testing. Control experiments in which class labels were scrambled produced sharply lower accuracies, confirming that the models were learning genuine biological patterns rather than statistical artifacts.

The results were striking in their consistency. Across every model, the hierarchy of performance was the same: passive membrane features yielded the best classification, followed by action potential features, with synaptic mEPSC features trailing behind. The random forest classifier trained on passive properties reached 93.81 percent accuracy, and gradient boosting came in a close second at 93.55 percent, while logistic regression struggled at 68.53 percent. For action potential features, k-nearest neighbors led at 83.36 percent. Even the weakest feature set, the synaptic events, allowed the random forest to reach nearly 70 percent accuracy, well above chance. Statistical analysis confirmed that accuracy differed significantly by model type, feature set, and their interaction, and the substantial gap between linear models and ensemble methods on passive properties suggests that the hormone-driven changes carve out complex, non-linear boundaries in the data.

Feature importance analysis added biological depth to the machine learning results. For the passive set, input resistance stood out as by far the most discriminative property, nearly doubling the importance of the next feature, the membrane time constant. This makes intuitive sense given prior experiments showing that estradiol replacement in ovariectomized animals directly modulates nucleus accumbens neuron input resistance. In the action potential set, burst length was the most informative feature, followed by spike frequency and average spike amplitude. Most unexpectedly, for the synaptic set the strongest signal came not from the synaptic events themselves but from the baseline current measured just before each event, which likely reflects cycle-dependent differences in resting membrane potential rather than changes in synaptic transmission per se.

Per-class recall analysis then revealed which phases carry the most distinctive electrical signatures. Late proestrus, the brief window when estrogen and progesterone surges peak, proved the most easily identified phase of all, with recall values of 99.92 percent on passive data, 93.84 percent on action potential data, and 86.28 percent on synaptic data. This aligns neatly with the intense endocrine activity of that phase, when hormones act on the nucleus accumbens to increase exploratory behavior, reduce anxiety-like behavior, and boost motivation to mate. In contrast, the estrus phase was consistently the hardest to classify, with recall as low as 43.61 percent on action potential data, perhaps reflecting a transitional period when some hormone effects persist while others fade back toward the diestrus phenotype.

Latent space visualization using t-SNE independently confirmed the classifier results, showing clear separation of late proestrus from the other three phases in both the passive and action potential feature sets, while the synaptic features produced more diffuse, overlapping clusters. The authors also addressed a methodological concern: because more recordings were available for late proestrus, the dataset was imbalanced. Retraining the models on an equalized subset produced no significant change in accuracy, indicating the results were not an artifact of unequal class sizes. The main limitations acknowledged by the team include the a priori selection of electrophysiological features and the inherently low throughput of manual patch clamp recordings.

Beyond the immediate findings, the study signals a broader shift in how neuroscientists might handle the long-neglected variability introduced by hormonal cycles. Most published neuroscience studies document animal sex but not cycle phase, meaning vast archives of existing electrophysiological data may contain hidden hormonal structure that no one looked for. The authors suggest that models like theirs could one day be applied retrospectively to such datasets, potentially assigning cycle phases and dramatically increasing the value of previously collected data. They also envision extending the approach to other brain regions such as the caudate-putamen, hippocampus, and hypothalamus, to other rhythmic biological processes such as circadian and seasonal cycles, and to combined analyses incorporating morphology and transcriptomics. In demonstrating that machine learning can read a hormonal clock out of a single neuron’s electrical whisper, the work positions these algorithms not merely as classifiers but as interpretive instruments for probing how neuromodulation sculpts the fundamental computations of the brain.

Subject of Research: Machine learning classification of neural electrophysiological differences across estrous cycle phases in rat nucleus accumbens medium spiny neurons

Article Title: Application of Machine Learning Models to Identify Differences in Neural Electrophysiological Properties Across Estrous Cycle Phases

Article References: Raina, A., & Meitzen, J. (2026). Application of Machine Learning Models to Identify Differences in Neural Electrophysiological Properties Across Estrous Cycle Phases. Neuroinformatics, 24(3), Article 57. https://doi.org/10.1007/s12021-026-09814-0

Image Credits: AI Generated

DOI: 10.1007/s12021-026-09814-0

Keywords: machine learning, estrous cycle, nucleus accumbens, medium spiny neurons, electrophysiology, neuroinformatics, random forest, input resistance, action potentials, mEPSC, estradiol, neuromodulation