Really beautiful rigorous work from the Liston and Levitz Labs combining circuit, molecular, and behavioral approaches to dissect the role of μ-opioid receptors (MOR/OPRM1) in ketamine’s antidepressant effects
One of the central findings is that MORs are enriched in somatostatin (SST) interneurons in cortex, and that these cells play a key role in mediating ketamine’s behavioral effects
There is a lot to like here: the study is technically sophisticated, the circuit logic is compelling, and the SST-dependent mechanism is supported by multiple independent experiments
But stepping back — it’s also critical that we get the cell-type distribution of OPRM1 correct, because this has direct implications for how we think about:
• opioid analgesia
• opioid use disorder
• cortical circuit modulation
• and the design of next-generation MOR-targeting therapeutics
And on that specific point, I think there is a non-trivial dataset interpretation issue worth discussing
The key claim relies heavily on older SMART-seq datasets (ACA/ALM; ~5K cells) from the Allen Institute
These datasets were incredibly important at the time — but they were not designed as unbiased quantitative cellular censuses
A few technical considerations that matter a lot for a gene like Oprm1:
1) Sampling design (not a census)
These datasets rely on:
• targeted dissections (ACA/ALM only)
• FACS-based sorting
• heavy use of transgenic driver lines
This means the data reflect what was selected, not necessarily the true population distribution
2) Cell numbers and statistical power
~5,000 ACA cells total. Once subdivided across excitatory subclasses, power drops quickly
For low-abundance GPCR transcripts like Oprm1, this creates:
• dropout sensitivity
• threshold artifacts
• unstable “presence/absence” calls
• inflated apparent enrichment in small populations
3) SMART-seq recovery biases
Requires intact dissociated cells:
• large projection neurons underrepresented
• pyramidal neurons more fragile
• interneurons often more recoverable
4) Quantification annotation (2018-era pipelines)
Older gene models isoform collapsing exon/intron handling can all affect detection of isoform-complex GPCRs like Oprm1
Now contrast that with newer datasets:
The Allen Brain Cell (ABC) Atlas includes ~4 million cells across the entire mouse brain, with:
• orders-of-magnitude larger sampling
• stable estimates across cortical subclasses
• improved taxonomy
• spatial integration
And importantly — across this dataset, OPRM1 signal is not restricted to SST interneurons
We see the same pattern across:
• our own single-nucleus RNA-seq (100k cells)
• Allen ABC Atlas (4 million cells)
• multiplexed FISH
• immunohistochemistry
• n=4 Oprm1-Cre mouse lines
• MORp viral promoter strategies
Across all of these orthogonal approaches, the result is highly consistent:
→ OPRM1/MOR is enriched in glutamatergic cortical populations, not exclusively confined to SST interneurons
So how do we reconcile this?
The most parsimonious explanation is not that either dataset is “wrong,” but that they are answering different questions under different constraints:
SMART-seq (2018):
→ high depth, small N, targeted sampling
→ vulnerable for sparse genes
ABC-scale (modern):
→ massive N, robust population estimates
→ better suited for cell-type distribution
For genes like Oprm1 — low abundance, heterogeneous, biologically critical —
scale sampling design cross-modal validation are decisive
The broader point:
As a field, we need to be careful about making strong claims about cell-type specificity of neuromodulatory receptors based on early-generation datasets — especially when:
• atlas-scale data now exist
• multiple orthogonal methods converge
• and the implications extend to therapeutics
None of this detracts from the importance of SST interneurons in MOR-dependent circuit function — that biology may be very real and important
But the global cortical distribution of OPRM1 appears broader, and heavily includes glutamatergic neurons
And that distinction matters for new appraoches for depression, pain and OUD