We address the problem of inferring the number of independently blinking fluorescent light emitters, when only their combined intensity contributions can be observed. This problem occurs regularly in light microscopy of objects smaller than the diffraction limit, where one wishes to count the number of fluorescently labeled subunits. Our proposed solution directly models the photophysics of the system, as well as the blinking kinetics of the fluorescent emitters as a fully differentiable hidden Markov model, estimating a posterior distribution of the total number of emitters. We show that our model is more accurate and increases the range of countable subunits by a factor of 2 compared to current state-of-the-art methods. Furthermore, we demonstrate that our model can be used to investigate the effect of blinking kinetics on counting ability and therefore can inform optimal experimental conditions.

Synaptic plasticity alters neuronal connections in response to experience, which is thought to underlie learning and memory. However, the loci of learning-related synaptic plasticity, and the degree to which plasticity is localized or distributed, remain largely unknown. Here we describe a new method, DELTA, for mapping brain-wide changes in synaptic protein turnover with single-synapse resolution, based on Janelia Fluor dyes and HaloTag knock-in mice. During associative learning, the turnover of the ionotropic glutamate receptor subunit GluA2, an indicator of synaptic plasticity, was enhanced in several brain regions, most markedly hippocampal area CA1. More broadly distributed increases in the turnover of synaptic proteins were observed in response to environmental enrichment. In CA1, GluA2 stability was regulated in an input-specific manner, with more turnover in layers containing input from CA3 compared to entorhinal cortex. DELTA will facilitate exploration of the molecular and circuit basis of learning and memory and other forms of plasticity at scales ranging from single synapses to the entire brain.

The spatial distribution of diverse biomolecules in multicellular organisms is essential for their physiological functions. High-throughput in situ mapping of biomolecules is crucial for both basic and medical research, and requires high scanning speed, spatial resolution, and chemical sensitivity. Here, we developed a Tissue Expansion method compatible with matrix-assisted laser desorption/ionization Mass spectrometry Imaging (TEMI). TEMI reaches single-cell spatial resolution without sacrificing voxel throughput and enables the profiling of hundreds of biomolecules, including lipids, metabolites, peptides (proteins), and N-glycans. Using TEMI, we mapped the spatial distribution of biomolecules across various mammalian tissues and uncovered metabolic heterogeneity in tumors. TEMI can be easily adapted and broadly applied in biological and medical research, to advance spatial multi-omics profiling.

Limited color channels in fluorescence microscopy have long constrained spatial analysis in biological specimens. Here, we introduce cycle Hybridization Chain Reaction (HCR), a method that integrates multicycle DNA barcoding with HCR to overcome this limitation. cycleHCR enables highly multiplexed imaging of RNA and proteins using a unified barcode system. Whole-embryo transcriptomics imaging achieved precise three-dimensional gene expression and cell fate mapping across a specimen depth of ~310 μm. When combined with expansion microscopy, cycleHCR revealed an intricate network of 10 subcellular structures in mouse embryonic fibroblasts. In mouse hippocampal slices, multiplex RNA and protein imaging uncovered complex gene expression gradients and cell-type-specific nuclear structural variations. cycleHCR provides a quantitative framework for elucidating spatial regulation in deep tissue contexts for research and potentially diagnostic applications.

bioRxiv preprint: 10.1101/2024.05.17.594641

Effective classification of neuronal cell types requires both molecular and morphological descriptors to be collected in situ at single cell resolution. However, current spatial transcriptomics techniques are not compatible with imaging workflows that successfully reconstruct the morphology of complete axonal projections. Here, we introduce a new methodology that combines tissue clearing, submicron whole-brain two photon imaging, and Expansion-Assisted Iterative Fluorescence In Situ Hybridization (EASI-FISH) to assign molecular identities to fully reconstructed neurons in the mouse brain, which we call morphoFISH. We used morphoFISH to molecularly identify a previously unknown population of cingulate neurons projecting ipsilaterally to the dorsal striatum and contralaterally to higher-order thalamus. By pairing whole-brain morphometry, improved techniques for nucleic acid preservation and spatial gene expression, morphoFISH offers a quantitative solution for discovery of multimodal cell types and complements existing techniques for characterization of increasingly fine-grained cellular heterogeneity in brain circuits.Competing Interest StatementThe authors have declared no competing interest.

We address the problem of inferring the number of independently blinking fluorescent light emitters, when only their combined intensity contributions can be observed at each timepoint. This problem occurs regularly in light microscopy of objects that are smaller than the diffraction limit, where one wishes to count the number of fluorescently labelled subunits. Our proposed solution directly models the photo-physics of the system, as well as the blinking kinetics of the fluorescent emitters as a fully differentiable hidden Markov model. Given a trace of intensity over time, our model jointly estimates the parameters of the intensity distribution per emitter, their blinking rates, as well as a posterior distribution of the total number of fluorescent emitters. We show that our model is consistently more accurate and increases the range of countable subunits by a factor of two compared to current state-of-the-art methods, which count based on autocorrelation and blinking frequency, Further-more, we demonstrate that our model can be used to investigate the effect of blinking kinetics on counting ability, and therefore can inform experimental conditions that will maximize counting accuracy.

Expansion microscopy (ExM) is a powerful technique to overcome the diffraction limit of light microscopy that can be applied in both tissues and cells. In ExM, samples are embedded in a swellable polymer gel to physically expand the sample and isotropically increase resolution in x, y, and z. By systematic exploration of the ExM recipe space, we developed a novel ExM method termed Ten-fold Robust Expansion Microscopy (TREx) that, as the original ExM method, requires no specialized equipment or procedures. TREx enables ten-fold expansion of both thick mouse brain tissue sections and cultured human cells, can be handled easily, and enables high-resolution subcellular imaging with a single expansion step. Furthermore, TREx can provide ultrastructural context to subcellular protein localization by combining antibody-stained samples with off-the-shelf small molecule stains for both total protein and membranes.

The fluorescent glutamate indicator iGluSnFR enables imaging of neurotransmission with genetic and molecular specificity. However, existing iGluSnFR variants exhibit low in vivo signal-to-noise ratios, saturating activation kinetics and exclusion from postsynaptic densities. Using a multiassay screen in bacteria, soluble protein and cultured neurons, we generated variants with improved signal-to-noise ratios and kinetics. We developed surface display constructs that improve iGluSnFR's nanoscopic localization to postsynapses. The resulting indicator iGluSnFR3 exhibits rapid nonsaturating activation kinetics and reports synaptic glutamate release with decreased saturation and increased specificity versus extrasynaptic signals in cultured neurons. Simultaneous imaging and electrophysiology at individual boutons in mouse visual cortex showed that iGluSnFR3 transients report single action potentials with high specificity. In vibrissal sensory cortex layer 4, we used iGluSnFR3 to characterize distinct patterns of touch-evoked feedforward input from thalamocortical boutons and both feedforward and recurrent input onto L4 cortical neuron dendritic spines.

The central nucleus of the amygdala (CEA) is a brain region that integrates external and internal sensory information and executes innate and adaptive behaviors through distinct output pathways. Despite its complex functions, the diversity of molecularly defined neuronal types in the CEA and their contributions to major axonal projection targets have not been examined systematically. Here, we performed single-cell RNA-sequencing (scRNA-Seq) to classify molecularly defined cell types in the CEA and identified marker-genes to map the location of these neuronal types using expansion assisted iterative fluorescence in situ hybridization (EASI-FISH). We developed new methods to integrate EASI-FISH with 5-plex retrograde axonal labeling to determine the spatial, morphological, and connectivity properties of ∼30,000 molecularly defined CEA neurons. Our study revealed spatio-molecular organization of the CEA, with medial and lateral CEA associated with distinct cell families. We also found a long-range axon projection network from the CEA, where target regions receive inputs from multiple molecularly defined cell types. Axon collateralization was found primarily among projections to hindbrain targets, which are distinct from forebrain projections. This resource reports marker-gene combinations for molecularly defined cell types and axon-projection types, which will be useful for selective interrogation of these neuronal populations to study their contributions to the diverse functions of the CEA.

Cells regulate function by synthesizing and degrading proteins. This turnover ranges from minutes to weeks, as it varies across proteins, cellular compartments, cell types, and tissues. Current methods for tracking protein turnover lack the spatial and temporal resolution needed to investigate these processes, especially in the intact brain, which presents unique challenges. We describe a pulse-chase method (DELTA) for measuring protein turnover with high spatial and temporal resolution throughout the body, including the brain. DELTA relies on rapid covalent capture by HaloTag of fluorophores that were optimized for bioavailability in vivo. The nuclear protein MeCP2 showed brain region- and cell type-specific turnover. The synaptic protein PSD95 was destabilized in specific brain regions by behavioral enrichment. A novel variant of expansion microscopy further facilitated turnover measurements at individual synapses. DELTA enables studies of adaptive and maladaptive plasticity in brain-wide neural circuits.