Research Article

Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region

See allHide authors and affiliations

Science  16 Nov 2018:
Vol. 362, Issue 6416, eaau5324
DOI: 10.1126/science.aau5324

Mapping the brain, one neuron at a time

Spatial transcriptomics can link molecularly described cell types to their anatomical positions and functional roles. Moffitt et al. used a combination of single-cell RNA-sequencing and MERFISH (multiplexed error-robust fluorescence in situ hybridization) to map the identity and location of specific cell types within the mouse preoptic hypothalamus and surrounding areas of the brain (see the Perspective by Tasic and Nicovich). They related these cell types to specific behaviors via gene activity. The approach provides an unbiased description of cell types of the preoptic area, which are important for sleep, thermoregulation, thirst, and social behavior.

Science, this issue p. eaau5324; see also p. 749

Structured Abstract

INTRODUCTION

A mechanistic understanding of brain function requires the identification of distinct cell types in the brain at a molecular, spatial, and functional level. The preoptic region of the hypothalamus comprises multiple nuclei and controls many social behaviors and homeostatic functions. Discrete neuronal types within the preoptic region have been associated with specific hypothalamic behaviors and homeostatic controls, yet the organizational principles of the underlying circuits remain elusive. Further progress requires methods that can identify molecularly distinct cell types and map their spatial and functional organization in the tissue.

RATIONALE

Single-cell RNA sequencing (scRNA-seq) has revolutionized the understanding of many tissues by allowing a systematic, genome-wide molecular identification of cell types. However, scRNA-seq requires cell dissociation, leading to a loss of spatial context that is essential to understand the cellular architecture of brain circuits. Image-based approaches to single-cell transcriptomics enables gene expression profiling of individual cells within their native tissue and offers opportunities for simultaneous in situ cell-type identification and spatial mapping, as well as functional characterization when combined with activity marker imaging. The combination of these complementary techniques would allow us to generate a molecular inventory of neuronal types while mapping their spatial and functional organization.

RESULTS

We combined scRNA-seq and multiplexed error robust fluorescence in situ hybridization (MERFISH), a single-cell transcriptome imaging method, to investigate the molecular, spatial, and functional organization of the mouse hypothalamic preoptic region. We profiled ~31,000 cells using scRNA-seq and imaged ~1.1 million cells within intact tissues using MERFISH. Our data revealed a remarkable diversity of neurons in this region, comprising ~70 different neuronal populations, many of which were previously unknown. These neuronal types exhibited distinct neuromodulatory signatures and revealed a striking heterogeneity within cell populations that were previously thought to be functionally unitary. MERFISH measurements further allowed us to map the spatial organization of these neuronal types, determine the cellular composition of distinct nuclei, and provide insights into the functional organization of neuron populations, including topographical relationships that underlie sex hormone signaling.

Last, we combined MERFISH with immediate-early-gene expression imaging to identify specific neuronal populations activated by social behaviors, including parenting, mating, and aggression. Several neuronal populations were selectively activated in each of these behaviors, supporting the notion that transcriptionally distinct neuronal types control specific hypothalamic functions. We identified a core neuronal population activated in all animals that exhibit parenting, as well as cell populations differentially activated in mothers and fathers, providing insights into how physiological state may affect parental behavior. Moreover, we identified cells associated with sexual behavior in males and females as well as male aggression toward infants and conspecific males.

CONCLUSION

By combining MERFISH with scRNA-seq, we have revealed the molecular, spatial, and functional organization of neurons within the hypothalamic preoptic region. These results provide a framework for mechanistic investigation of behavior circuits with high molecular and spatial resolution and opens avenues for identifying and mapping cell types in a diverse range of tissues and organisms.

In situ single-cell profiling reveals the molecular and cellular organization of the hypothalamic preoptic region.

The combination of MERFISH with scRNA-seq to profile the gene expression of 1 million cells in situ revealed ~70 neuronal populations in the preoptic region, each with distinct molecular signatures and spatial organizations, providing insights into neuromodulatory signaling pathways. Further combination with activity marker imaging led to the identification of discrete neuronal types activated by key social behaviors, including parenting, aggression, and mating.

Abstract

The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei—including the molecular identity, spatial organization, and function of distinct cell types—is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.

A mechanistic understanding of brain function requires a systematic assessment of cell types and their spatial organization, connectivity, and functional properties. A case in point is the preoptic region of the hypothalamus, which comprises multiple nuclei and controls essential social behaviors such as parenting, mating, and aggression as well as homeostatic functions such as thermoregulation, thirst, and sleep (1, 2). Because these are evolutionarily conserved functions, it has been proposed that the associated neural circuits are genetically defined and thus composed of transcriptionally distinct neuronal types (13). Indeed, several neuronal populations within the preoptic region, each defined by discrete molecular markers, have been linked to distinct behavioral and homeostatic functions (411). However, the number of cell types present in the preoptic region as well as their molecular signatures, spatial organizations, and functional roles remain unclear, hampering our ability to investigate the underlying neural circuits.

Single-cell RNA-sequencing (scRNA-seq) provides a powerful means for the identification of cell types and cell states through genome-wide expression profiling of individual cells, offering rich insights into the cellular diversity of many tissues, including the brain (1215). However, scRNA-seq requires cell dissociation and thus results in the loss of the spatial context of cells that is critical for understanding tissue function (15, 16). Recently, image-based single-cell transcriptomic approaches have been developed that quantify gene expression by directly imaging individual RNA molecules within intact cells and tissues with multiplexed fluorescence in situ hybridization (FISH) or in situ sequencing (15, 1722). These approaches offer new opportunities to identify cell populations within complex tissues while simultaneously mapping their spatial organization and uncovering their functions by combining gene expression profiling with imaging of activity markers, such as the induction of immediate early genes (IEGs) (22, 23). Among these, multiplexed error-robust FISH (MERFISH) detects individual RNA molecules with single-molecule FISH (smFISH) (24, 25) and uses error-robust barcoding, combinatorial labeling, and sequential imaging to multiplex smFISH measurements, enabling transcriptome-scale RNA imaging of individual cells in situ (20, 26).

We developed a MERFISH-based imaging and analysis platform for in situ cell-type identification and mapping and used this approach, in combination with scRNA-seq, to create a cell atlas of the preoptic region of the mouse hypothalamus. We used scRNA-seq to catalog cell populations and identify their marker genes. We then performed MERFISH imaging of these marker genes together with genes of known functional importance to identify cell populations and map their spatial organization in situ. Last, we combined MERFISH with measurements of IEG expression in order to identify discrete cell populations activated by specific social behaviors—including parenting, aggression, and mating—in both sexes and different physiological states.

Results

scRNA-seq of the preoptic region

We dissected a rostral part of the mouse hypothalamus that contains the preoptic region (Fig. 1A)—the medial preoptic area (MPOA) and surrounding nuclei (~2.5 by 2.5 by 1.1 mm, Bregma +0.5 to –0.6)—from adult female and male brains and dissociated the tissue using a custom protocol that improved cell survival and capture (fig. S1). We collected scRNA-seq profiles from 31,299 cells across three replicates of each sex using droplet-based scRNA-seq (2729).

Fig. 1 scRNA-seq of the preoptic region in the mouse hypothalamus.

(A) Schematic of the preoptic region of the hypothalamus. Magenta boxes indicate the area dissected for scRNA-seq (Bregma +0.5 to –0.6). (B) t-distributed stochastic neighbor embedding (tSNE) for all cells and inhibitory and excitatory neurons, with cells colored by cluster. Numbers superimposed on the tSNE indicate the cluster ID. Total cell numbers for each tSNE plot are indicated. NFO, newly formed oligodendrocytes; OPC, oligodendrocyte progenitor cells; MO, mature oligodendrocytes. (C) Heat map of z-scores of expression for select genes within inhibitory neuronal clusters. Clusters are organized on the basis of the hierarchical tree constructed with expression in principal component space, with some of the genes differentially expressed between branches indicated (blue). The nomenclature of clusters uses a numeric indicator of excitatory or inhibitory cluster followed by one or two marker genes, with the first marker typically a neuromodulator (29). Inhibitory and excitatory clusters that lack a notable neuromodulator marker gene were designated as Gaba and Glut, respectively, with an additional marker gene to help differentiate among these clusters when possible. Cluster names are colored according to the first gene. Predicted anatomical locations for the clusters are listed on the tree, and the unlabeled lines indicate that such prediction was not possible. Thick black lines underscore clusters grouped by common neuropeptide expression. (D) As in (C) but for excitatory neurons. The hybrid neuronal clusters h1/h2 and h3 are listed in (C) and (D), respectively, because they were initially classified as inhibitory and excitatory, respectively. (E) –log10(P value) for the enrichment of gene categories in differentially expressed genes that mark neuronal clusters calculated based on a gene-set enrichment analysis as shown in fig. S6. *P < 0.05.

We used unsupervised, graph-based, community-detection methods (28, 30, 31) modified by us (fig. S2) to cluster cells (29). This led to the delineation of major cell classes, including inhibitory and excitatory neurons, microglia, astrocytes, immature oligodendrocytes (newly formed oligodendrocytes and oligodendrocyte progenitor cells), mature oligodendrocytes, ependymal cells, endothelial cells, fibroblasts, macrophages, and mural cells, as well as subdivisions within these cell classes (Fig. 1B and table S1).

Further clustering of inhibitory neurons (15,042 cells) and excitatory neurons (3511 cells) separately revealed 43 and 23 subpopulations, respectively (Fig. 1B; fig. S3, A and B; and tables S1 and S2). Hereafter, we denote excitatory and inhibitory neuronal clusters as e1, e2, …, and i1, i2, …, respectively. We also provide specific names for these clusters based on marker genes (Fig. 1, C and D, and figs. S4 and S5, the latter emphasizing neuropeptide expression) (29).

Although the majority of the identified clusters expressed either excitatory or inhibitory neuronal markers, we observed expression of the γ-aminobutyric acid (GABA) synthetic genes Gad1 and Gad2 in many excitatory neuronal clusters classified on the basis of expression of Vglut2 (Slc17a6), with Gad2 expression being particularly widespread (fig. S3C). By contrast, very few Slc17a6-positive clusters expressed the GABA transporter gene Vgat (Slc32a1). These data suggest that Slc17a6 and Slc32a1 are better discriminators for excitatory versus inhibitory neurons, corroborating evidence from other brain areas (32). Cells in two neuronal clusters originally designated as inhibitory and one originally designated as excitatory coexpressed Slc17a6 (or Slc17a8, vGlut3) and Slc32a1. These cells were unlikely to be a clustering artifact because individual cells coexpressed both markers, nor did they correspond to doublets (29); hence, they potentially represent hybrid neurons capable of GABA/glutamate corelease, as characterized in the hypothalamus and a few other brain regions (3234). We denote these clusters as h1, h2, and h3 (Figs. 1, C and D, and fig. S3C).

To determine the gene categories that best discriminate neuronal clusters, we examined the top five most differentially expressed genes in each cluster and observed enrichment for neuropeptides and molecules involved in neuromodulator production and transport, as well as for transcription factors, but not for neuromodulator (neuropeptide and hormone) receptors. Quantitative analyses of enrichment profiles of these three gene classes among differentially expressed genes further support this notion (Fig. 1E, fig. S6, and table S3). Neuromodulator receptors did discriminate some clusters (for example, Npr1, Rxfp1, Brs3, and Drd1) (Fig. 1, C and D, and figs. S4 and S5). However, on average, neuromodulator receptors were expressed more widely and at lower levels than neuromodulators and transcription factors, limiting their use as potential markers for functional studies. Most clusters were discriminated by combinations of genes rather than by single markers.

Hierarchical tree analyses (29) showed that inhibitory neuronal clusters that express a common neuromodulator were often grouped together on the tree—for example, clusters expressing Avp, Gal, Crh, Tac1, and Sst (Fig. 1C)—suggesting potential functional or developmental commonality among them. By contrast, neuromodulators largely failed to group excitatory neuronal clusters (Fig. 1D). Instead, predicted locations of individual clusters on the basis of spatial expression patterns of their marker genes observed in the Allen Brain Atlas (35) and our own in situ hybridization data (fig. S7) suggest that excitatory clusters tended to be grouped on the tree by anatomical structures or nuclei (Fig. 1D). For example, markers of clusters e4, e2, e21, h3, and e17 defined a node in the tree located in the PVN and adjacent nuclei (MPN, PaAP, BAC, and BNST), markers of node-sharing clusters e13 and e7 placed these populations in the MnPO/AvPe/VMPO region, whereas markers of e12, e6, e5, and e1 placed these cells in the MPN/MPA region (Fig. 1D) (full names of the nuclei described in this work are provided in table S4). We thus hypothesize that excitatory neuron types tend to be spatially segregated in distinct anatomical structures of the preoptic region, in a manner similar to the spatial segregation of different types of excitatory neurons in various layers of the cortex (36). Similar analysis with the inhibitory neuronal tree suggests that although some groups of clusters were defined by spatially restricted transcription factor expression—for example, Six6 marking the SCN (Fig. 1C)—such spatial grouping of transcriptionally similar clusters appeared to be less pronounced than with excitatory clusters. Additionally, transcription factors tended to mark groups of neuronal clusters further subdivided by neuromodulator expression (Fig. 1, C and D), which is consistent with earlier reports of hypothalamic parcellation by transcription factors during early development (37).

Specific neuronal clusters identified with scRNA-seq

Previous studies of the preoptic region have defined cell populations associated with the regulation of specific homeostatic and behavioral functions on the basis of the expression of one or more marker genes (table S5). Clusters that express these marker combinations were identified in our scRNA-seq data (figs. S4 and S5), together with many previously unknown cell populations. Moreover, we uncovered a high level of molecular heterogeneity among a number of previously reported singular cell types, thus partitioning them into multiple distinct populations, as illustrated below on specific examples.

The neuropeptide galanin (Gal) has been associated with behaviorally relevant cell populations of the preoptic region (4, 5, 38) in the MPOA (parenting and feeding) (5, 38) and VLPO (sleep) (4). Our scRNA-seq data revealed seven neuronal clusters that were statistically enriched in Gal expression, each characterized by distinct marker genes (Fig. 2A) validated with two-color in situ hybridization (fig. S7A). These clusters were each associated with different hormonal modulations, ranging from cluster i20:Gal/Moxd1, predicted to lie in the sexually dimorphic nucleus of the POA (Fig. 1C) and expressing a wide range of sex steroid and neuropeptide receptors, to cluster e24:Gal/Rxfp1, expressing no sex steroid receptor (Fig. 2A).

Fig. 2 scRNA-seq identifies subdivisions of cells that express markers previously associated with single neuronal populations.

(A to C) Expression distributions of selected marker genes and genes of interest in all neuronal clusters that are statistically enriched [Model-based Analysis of Single-cell Transcriptomics (MAST) (75), false discovery rate <0.01] in (A) galanin (Gal), (B) tyrosine hydroxylase (Th), or (C) Bdnf and Adcyap1. Gene names in black indicate differentially expressed genes for each selected neuronal cluster. Gene names in blue indicate inhibitory (Gad1, Gad2, Slc32a1) and excitatory (Slc17a6) neuronal markers, as well as dopaminergic markers (Ddc, Slc6a3, and Slc18a2). Gene names in green indicate sex hormone receptors. The y axis on each violin plot depicts the log transformed counts with the range set to the 95% expression quantile of the cluster with the highest expression (29). The sizes of red, cyan, and yellow circles correspond to the cell abundance of the inhibitory, excitatory, and hybrid clusters, respectively.

Second, cells that express tyrosine hydroxylase (Th), a key enzyme involved in catecholamine synthesis, have been viewed as a single population involved in several social behaviors (6, 39). We identified six Th-enriched neuronal clusters (Fig. 2B and fig. S7B), among which only i16:Gal/Th and i38:Kiss1/Th expressed both Dopa decarboxylase (Ddc) and the vesicular monoamine transporter Vmat2 (Slc18a2), genes required for dopaminergic function (Fig. 2B).

Last, the neuropeptide adenylate cyclase activating polypeptide 1 (Adcyap1) and brain-derived neurotrophic factor (Bdnf) have recently been identified as combined markers for preoptic neurons sensing warm temperature (8). Our data revealed nine Adcyap1- and Bdnf-enriched clusters (Fig. 2C). Although the warm-sensitive neurons have been previously considered as inhibitory neurons on the basis of their functional properties and expression of Gad2, all nine Adcyap1- and Bdnf-enriched clusters identified here coexpressed Gad2 and Slc17a6, and only one of them also expressed Slc32a1 (Fig. 2C), identifying these clusters as excitatory or hybrid neurons. We further identified one of these clusters as representing warm-sensitive neurons with the help of MERFISH. A recent study has revealed that a neuronal population that controls thirst-motivated behavior also expresses Adcyap1 and Bdnf (10), further supporting the notion that Adcyap1 and Bdnf are imperfect markers for warm-sensitive cells.

MERFISH measurements of the preoptic region

Next, we performed MERFISH measurements of the preoptic region (1.8 by 1.8 by 0.6 mm, Bregma +0.26 to –0.34), within the area characterized with scRNA-seq, targeting a set of 155 genes (Fig. 3A and table S6) (29). These genes were composed of two groups: (i) 85 preselected genes that were either known markers for major cell classes or relevant to neuronal functions of the hypothalamus, such as neuropeptides and neuromodulator receptors, and (ii) 70 additional genes that were identified with scRNA-seq as neuronal cluster markers but not already included in the 85 preselected genes. Among these 155 genes, 135 genes were imaged by using combinatorial smFISH with an error-robust barcoding scheme, as demonstrated previously for MERFISH (20, 26, 40). The remaining 20 genes were relatively short and/or expressed at high levels, which is challenging for combinatorial smFISH detection, and hence were measured in sequential rounds of multicolor FISH after the combinatorial run. The sexually dimorphic expression previously reported for 11 genes (41, 42) was confirmed here (fig. S8).

Fig. 3 Major cell classes and their spatial organizations in the preoptic region as revealed with MERFISH.

(A) (Left) Schematic of the MERFISH measurements. Combinatorial smFISH imaging was used to identify 135 genes, followed by sequential rounds of two-color FISH to identify 20 additional genes. Total polyadenylated mRNA and nuclei costains then allowed cell boundary segmentation. (Top right) Pseudo-colored dots marking localizations of individual molecules of eight example RNA species, each marking a distinct major cell class, in a 10-μm-thick, 1.8- by 1.8-mm slice. (Bottom right) Magnification of the white boxed region (left) and the total mRNA image and the segmented cell boundaries of the same region (right). The raw and decoded MERFISH images of the same field of view (FOV) for all 135 genes measured by using combinatorial smFISH are shown in fig. S9; the total mRNA and nuclei costain images and segmented cell boundaries for the same FOV are shown in fig. S10. The segmented cell boundaries represent the boundaries of the cell soma (29). A subset of identified RNA molecules fell outside these boundaries and are thus candidates for RNAs in neuronal or glial processes. (B) Expression of all genes measured with MERFISH for ~500,000 cells imaged in multiple naïve animals. Expression for each gene is normalized to the 95% expression quantile for that gene across all cells. Cells are grouped by major classes, and markers of each major cell class are listed on the right. OD, oligodendrocytes. (C) tSNE plot of these cells. (D) Pairwise Pearson correlation coefficients between the average expression profiles (in z-scores) of individual cell classes identified with MERFISH and scRNA-seq. (E) (Top) Spatial distribution of all major cell classes across sections at different anterior-posterior positions from a single female mouse. Cells are marked with cell segmentation boundaries and colored by cell classes as indicated. Six of the twelve 1.8- by 1.8-mm imaged slices are shown. The 0, 100, 200, 300, 400, and 500 μm labels indicate the distance from the anterior position (Bregma +0.26). (Bottom) Enlarged image of the slice at 400 μm from the anterior position (left) and a further magnified image of the region shown in the gray dashed box (right). Scale bars, 500 μm (left), 250 μm (right). (F) Spatial distributions of individual cell classes are shown as colored dots on the background of all cells shown as gray dots. Dashed ovals indicate several specific hypothalamic nuclei and are colored identically to the nuclei abbreviations listed to the right. BNST, bed nucleus of the stria terminalis; MPN, medial preoptic nucleus; MnPO, median preoptic nucleus; Pe, periventricular hypothalamic nucleus; AvPe, anteroventral periventricular nucleus; VMPO, ventromedial preoptic nucleus; VLPO, ventrolateral preoptic nucleus; PVA, paraventricular thalamic nucleus; PaAP, paraventricular hypothalamic nucleus, anterior parvicellular.

We sectioned the preoptic region into 60 evenly spaced slices along the anterior-posterior axis and performed three-dimensional MERFISH imaging on every fifth slice (29). Individual RNA molecules were clearly detected and identified (fig. S9), and individual cells were segmented based on 4′,6-diamidino-2-phenylindole (DAPI) and total mRNA staining (fig. S10) (29). In total, we profiled >400,000 cells from three to four replicates in naïve male and female animals, as well as >500,000 additional cells from three to five replicates of animals subjected to behavioral stimuli (29). MERFISH expression data showed high reproducibility between replicates (fig. S11A), good correlation with bulk RNA-seq data of the preoptic region (43) (fig. S11B), and a low false-detection rate (fig. S11C). For the targeted genes, MERFISH detected on average six- to eightfold more transcript copies per cell than did scRNA-seq (fig. S12, A to D), underscoring the high sensitivity of MERFISH.

We used an unsupervised, community-detection–based clustering approach similar to that applied to scRNA-seq data to identify transcriptionally distinct cell populations in MERFISH data (Fig. 3, B and C, and table S7) (29). MERFISH identified all major cell classes (Fig. 3, B and C), except for macrophages and fibroblasts, potentially because the corresponding marker genes were not included in the MERFISH gene library. The expression profiles of cell classes measured with MERFISH were strongly correlated with those determined by using scRNA-seq (Fig. 3D). However, the relative abundance of cells in various cell classes differed in the two datasets (fig. S12E). In particular, astrocytes, endothelial cells, and ependymal cells were depleted in our scRNA-seq data, presumably because of cell loss during tissue dissociation.

MERFISH also provided a direct measurement of the spatial distribution of major cell classes. As expected, mature oligodendrocytes were enriched in the anterior commissure and the fornix–major myelinated fiber tracts of the rostral hypothalamus, whereas immature oligodendrocytes, astrocytes, microglia, and endothelial cells were dispersed throughout (Fig. 3, E and F). Ependymal cells formed a single layer lining the more caudal aspects of the third ventricle, and mural cells were organized in vermiform structures that resemble blood vessels (Fig. 3, E and F). Notably, inhibitory and excitatory neurons exhibited distinct distributions (Fig. 3, E and F). Inhibitory neurons, the more abundant neuronal type in the preoptic region, were widely dispersed across this region but enriched in specific posterior nuclei, including the BNST and MPN. By contrast, excitatory neurons were specifically enriched in a few nuclei anteriorly but became more dispersed posteriorly and, in agreement with previous reports (44), were depleted in the posterior BNST.

MERFISH analyses of specific neuronal types

Clustering analyses of inhibitory neurons and excitatory neurons separately identified ~40 inhibitory and ~30 excitatory neuronal populations (Fig. 4, A and B, and tables S7 and S8). We investigated the impact of the number of genes used to cluster cells in MERFISH data and found that ~90% of the identified neuronal clusters were recovered by using the ~75 genes that were most informative among the 155 (fig. S13). Beyond this point, cluster recovery increased more slowly with the number of genes added (fig. S13). Hereafter, we denote excitatory and inhibitory neuronal clusters identified with MERFISH as E-1, E-2, … and I-1, I-2, …, respectively, and the one identified hybrid cluster as H-1.

Fig. 4 Neuronal clusters in the preoptic region as revealed with MERFISH.

(A and B) z-scores of expression profiles for (A) inhibitory and (B) excitatory neuronal clusters identified with MERFISH. Depicted are 100 random cells from each cluster. The neuronal clusters are organized on the basis of similarity in their expression profiles, as depicted by the dendrogram. The sizes of red, cyan, and yellow circles indicate the abundance of neuronal clusters, and only clusters with more than 100 cells are depicted. H-1 is grouped with the inhibitory clusters because it was initially classified as inhibitory neurons. (C) The pairwise Pearson correlation coefficients between the expression profile (in z-score) of the MERFISH and scRNA-seq clusters. The order of the clusters in (C) is not the same as in (A) and (B). (D) As in (C) but with only scRNA-seq cluster(s) most similar to each MERFISH cluster shown, identified as the cluster(s) with the highest Pearson correlation coefficient(s) (fig. S14 and table S9) (29). When multiple scRNA-seq clusters show statistically indistinguishable, highest correlation coefficients to a MERFISH cluster (29), all of them are indicated. scRNA-seq clusters outside the region imaged with MERFISH, as assessed by the expression patterns of the marker genes in the Allen Brain Atlas (35) and our own in situ data (fig. S7) (29), are excluded from this analysis (29). (E) Same as (D) but for clusters enriched in galanin (Gal).

The expression profiles of most neuronal clusters determined with MERFISH correlated well with those of scRNA-seq clusters (Fig. 4C and fig. S14, A and B). This observation allowed us to infer, for each MERFISH cluster, the putative corresponding or most similar scRNA-seq cluster(s), defined as the cluster(s) with the highest correlation coefficient(s) (Fig. 4D, fig. S14C, and table S9) (29), which could help expand our knowledge of the expression profiles of the MERFISH clusters. Similar correspondence was observed by using a neural network classifier (fig. S14, D and E). Correlations between MERFISH and scRNA-seq clusters were only moderately weaker than those between scRNA-seq clusters derived from bootstrapped replicates (fig. S14, F and G). Many MERFISH clusters had a distinct, most similar scRNA-seq cluster (Fig. 4D, fig. S14, and table S9). However, in some instances, multiple MERFISH clusters exhibited the highest correlation to the same scRNA-seq cluster; in addition, a small fraction of MERFISH clusters lacked a statistically significant correlation to any scRNA-seq cluster (Figs. 4D, fig. S14, and table S9). Both of these scenarios suggest that some clusters identified with MERFISH were not discriminated by scRNA-seq. Conversely, a small fraction of scRNA-seq clusters lacked a statistically significant correlation to any MERFISH cluster (Fig. 4D, fig. S14, and table S9), suggesting that these clusters were not identifiable by the MERFISH gene panel or were located outside the MERFISH-imaged area (29).

As a specific illustration, MERFISH identified 10 clusters enriched in Gal expression, some showing one-to-one correspondences to Gal-enriched scRNA-seq clusters (Fig. 4E; fig. S15, A and B; and table S9). We also observed instances in which two Gal-enriched MERFISH clusters putatively corresponded to the same Gal-enriched scRNA-seq cluster (for example, I-14 and I-16 to i16) (Fig. 4E and table S9), suggesting that MERFISH resolved subpopulations within the scRNA-seq cluster. Indeed, we identified two subsets of cells within the scRNA-seq cluster i16, each respectively expressing markers of the MERFISH clusters I-14 and I-16 (fig. S15, C and D). The calcitonin receptor (Calcr)– and bombesin receptor (Brs3)–positive I-14 and the Th-positive I-16 were found to be differentially activated in specific social behaviors, as described later (table S9), supporting the resolution of these cells into two distinct populations. We also observed a similar resolution of i8 into I-7 and I-31 (Fig. 4E; fig. S15, E and F; and table S9). A Gal-enriched scRNA-seq cluster could also be split into Gal-enriched and non–Gal-enriched MERFISH clusters [for example, i20 into Gal-enriched cluster I-34 and non–Gal-enriched clusters I-2 and I-32 (fig. S15, G and H, and table S9)].

Examination of the MERFISH or scRNA-seq clusters that were not discriminated by the other method showed several trends. Some of the MERFISH clusters not detected with scRNA-seq had relatively low abundance and thus might not be sufficiently represented in our scRNA-seq data, which profiled 4% as many neurons as we did in MERFISH. Some of the MERFISH clusters not discriminated by scRNA-seq had lowly expressed marker genes, which may not be reliably detected with scRNA-seq. Conversely, some scRNA-seq clusters not identified with MERFISH had marker genes that were not included in the MERFISH gene library. Some of the extremely low-abundance MERFISH or scRNA-seq clusters that lack correspondence may not represent well-identified clusters. These results thus demonstrate the complementary nature of MERFISH and scRNA-seq and an increased ability to characterize cell populations when both approaches are combined. Nevertheless, some clusters still exhibited heterogeneity in gene expression associated with distinct spatial locations (fig. S16), suggesting either spatial gradients in gene expression within the same cluster or the presence of unresolved cell subpopulations.

Spatial organization of specific neuronal cell types

Next, we examined the spatial distributions of individual neuronal clusters (Fig. 5A and figs. S17 and S18) within the framework of major anatomically defined nuclei of the preoptic region as depicted in Fig. 5B (45). About 30% of the MERFISH clusters were enriched primarily in a single nucleus (Fig. 5C, pink shading)—for example, cluster I-5 primarily in the VLPO and cluster E-9 in the PVA (Fig. 5, A to C)—whereas approximately half of the clusters were distributed over a few (two to four), often physically contiguous nuclei (Fig. 5C, unshaded clusters)—for example, cluster I-3 in the BNST-mv, PaAP, PS, and SHy and cluster I-12 in the StHy, MPA, and MPN (Fig. 5, A to C). This anatomical dispersion may reflect a similar function of the same cell type across distinct nuclei or a developmental relationship of spatially distinct cells. By contrast, a small fraction of the neuronal clusters were dispersed and not enriched in any given nucleus, such as I-21 and E-22 (Fig. 5A). Whereas most nuclei were populated by both excitatory and inhibitory neurons, the PVA and BAC only contained excitatory clusters (Fig. 5C, green shaded row), and the BNST-p and BNST-mv contained only inhibitory clusters (Fig. 5C, blue shaded rows), which is consistent with previous observations of high expression of Slc17a6 and Slc32a1 in these regions, respectively (35, 44).

Fig. 5 The spatial organization of neuronal clusters in the preoptic region.

(A) Spatial distribution of example neuronal clusters that are localized (top and middle) or dispersed (bottom). Depicted are six of the 12 slices imaged from a female mouse. Colored markers indicate cells of the specified neuronal clusters, and gray markers indicate all other neurons. Nuclei boundaries depicted in light gray are drawn according to (45) and aligned to the tissue slices according to the locations of landmarks, such as the anterior commissure, fornix, and ventricle. The 0, 100, 200, 300, 400, and 500 μm labels indicate the distance from the anterior position (Bregma +0.26). (B) Illustration of major hypothalamic nuclei spanning the imaged region and colored according to legend on the right (45). Nuclei abbreviations are as defined in Fig. 3F, and additionally, BAC, bed nucleus of the anterior commissure; LPO, lateral preoptic area; MPA, medial preoptic area; PS, parastrial nucleus; StHy, striohypothalamic nucleus; SHy, septohypothalamic nucleus; ACA, anterior commissure; Fx, fornix; 3V, third ventricle. Bregma locations are listed on top and the map at Bregma –0.22 is duplicated. (C) Summary of nuclei in which inhibitory (blue) or excitatory (green) neuronal clusters are enriched. Translucent horizontal bars indicate nuclei that contain only inhibitory (blue) or excitatory (green) clusters. Vertical pink bars highlight clusters primarily enriched in single nuclei. BNST-mal, BNST, medial division, anterolateral part; BNST-mv, BNST, medial division, ventral part; BNST-p, BNST, posterior part. (D and E) Analysis of spatial mixing of distinct neuronal clusters. We define the complexity of the neighborhood surrounding any given neuron as the number of distinct neuronal clusters present within that neighborhood, and the purity of that neighborhood as the fraction of all cells within the given neighborhood that are part of the most abundant cluster. Probability distributions of the complexity (D) and purity (E) of the 100-μm-radius neighborhood surrounding any given neuron are depicted.

Neuronal clusters of the preoptic region appeared highly intermixed, with multiple clusters occupying any given nucleus. To quantify the degree of intermixing, we calculated the neighborhood composition for each neuron. This analysis showed that each neighborhood contained multiple clusters and was typically not dominated by a single cell population (Fig. 5, D and E).

These direct spatial measurements allowed us to provide an anatomy-based taxonomy for the identified neuronal clusters, except for the dispersed clusters, which we named on the basis of marker genes (table S9). The putative correspondence between these MERFISH clusters and scRNA-seq clusters allowed us to further assess the spatial locations of scRNA-seq clusters and compare them with our earlier predictions shown in Fig. 1, C and D. In nearly all cases, the predicted locations matched or partially overlapped those of the corresponding MERFISH clusters (table S9), lending additional support to our earlier observations on the spatial relationship between transcriptionally similar clusters (Fig. 1, C and D). However, we caution that the predicted scRNA-seq cluster locations represent rough approximations because of the relatively low resolution of available in situ hybridization data and suffer from occasional ambiguity in the spatial patterns of some marker genes. Moreover, because of the remaining heterogeneity in some of the identified scRNA-seq and MERFISH clusters, the putative correspondence may only represent similarity between subsets of cells within these clusters.

Spatial and molecular organization of socially relevant cell populations

Neuromodulators, hormones, and associated signaling pathways play critical roles in hypothalamic functions, but analyses with cellular resolution have been limited owing to the low expression level of many of the corresponding receptors. MERFISH enabled us to examine the distribution of these genes throughout the preoptic region, providing functional insights into the associated cell populations.

Sex steroid hormones are essential to the development and modulation of social behaviors and reproduction. We examined the distribution of enzymes and receptors essential for steroid hormone signaling in the preoptic region. The enzyme aromatase (Cyp19a1) converts testosterone to estrogen and thus modulates steroid function (39). MERFISH revealed aromatase-enriched clusters with distinct repertoires of sex hormone receptors (Fig. 6, A and B). Several of these clusters (such as I-2, I-13, I-24, I-32, and E-12) expressed both androgen receptor (Ar) and estrogen receptor α (Esr1), suggesting that in these cell populations, circulating testosterone can be converted into estrogen and thus affect gene expression in a cell-autonomous manner through Esr1 activation. In addition, the aromatase-enriched clusters E-12 and I-2 substantially overlapped with locations of several Esr1-enriched clusters (Fig. 6, C to F), suggesting that estrogen synthesized by these aromatase-expressing cells may also act in a paracrine manner on cells of the nearby Esr1-enriched clusters, in addition to the autocrine signaling mode described above (Fig. 6G).

Fig. 6 Spatial and molecular organization of neuronal clusters enriched in genes relevant to social behaviors.

(A) Expression distributions of selected marker genes and genes of interest for neuronal clusters enriched in aromatase (Cyp19a1). Expression distributions are calculated as in Fig. 2. (B) Spatial distributions of neuronal clusters depicted in (A). Two of the 12 slices from a female mouse sample are depicted. Nuclei boundaries depicted in light gray are as defined in Fig. 5A. (C and D) As in (A) and (B) but for clusters enriched in estrogen receptor α (Esr1). (E) Schematic showing the nuclei spanned by individual clusters, as indicated by the color subdivisions of the rectangles, colored identically to the nuclei abbreviations listed below. The nuclei abbreviations are as defined in Figs. 3F and 5B. I-17 is not colored because it was found at the edge of our imaged region and falls outside of the boundaries of the nearest imaged nuclei, the VLPO (table S9). (F) Average overlap fraction between aromatase-enriched clusters and Esr1-enriched clusters for all measured animals. Cluster I-24 is enriched in both aromatase and Esr1 and listed only once. (G) Models of autocrine and paracrine signaling. Circulating testosterone (gray T) can activate cells expressing androgen receptor (AR) or, in cells expressing aromatase, can be converted to estrogen (orange E). Autocrine: In cells co-expressing aromatase and Esr1, estrogen produced in these cells can activate estrogen receptor (ERa) in the same cells. Paracrine: Estrogen produced by aromatase-enriched cells can activate ERa in nearby cells enriched with Esr1. (H) Comparison of the fraction of cells that belong to the specified neuronal clusters (I-15 or I-2) for all male (blue) and female (red) replicates as a function of the anterior-posterior position of the slices. Above each panel are the spatial distribution of the cluster in four slices from a single female (red) and male (blue) replicate. (I and J) As in (A) and (B) but for clusters enriched in oxytocin receptor (Oxtr). (K and L) As in (A) and (B) but for a cluster enriched in gonadotropin releasing hormone 1 (Gnrh1). MERFISH revealed 8, 15, and 19 aromatase-, Esr1- and Oxtr-enriched clusters, respectively, with only the seven most enriched clusters depicted for each.

Some Esr1-enriched and aromatase-enriched clusters exhibited differences in cell abundance of varying extent between males and females. For example, the Esr1-enriched cluster I-15 showed an appreciable enrichment in female animals, whereas the aromatase-enriched cluster I-2 showed a more modest male enrichment (Fig. 6H). I-2 overlapped with the sexually dimorphic nucleus of the preoptic area (SDN-POA), and its marker gene Cplx3 was co-expressed with MoxD1 [a canonical SDN-POA marker (46) not in our MERFISH gene library] in cells of the SDN-POA and BNST (fig. S19). However, I-2:BNST/StHy/MPN spatially extended beyond the boundaries of the SDN-POA, suggesting that it may contain unresolved subpopulations, which may partially mask a sexual dimorphism in SDN-POA cell populations.

We next considered clusters enriched in the expression of oxytocin receptor (Oxtr), an important modulator of social behaviors that exerts its effects broadly throughout the brain (47). Although the low expression level of Oxtr has previously made it challenging to identify oxytocin targets, the high sensitivity of MERFISH allowed us to detect enrichment of Oxtr in multiple clusters (Fig. 6, I and J). For example, the Oxtr-enriched BNST cluster I-24 coexpressed multiple sex hormone receptors as well as the neuropeptide Tac2 implicated in social isolation stress (48), aggression, and fear (48, 49), suggesting the involvement of oxytocin signaling in these functions. This cluster was specifically activated after pup-directed aggression by virgin males as described below, corroborating studies that implicate the BNST in this function (50). This highlights a seemingly paradoxical role for Oxtr in agonistic pup encounters versus its known roles in affiliative behaviors (47). Oxtr was also found in the VLPO cluster I-5, implying that oxytocin might have a role in the modulation of VLPO functions such as sleep or temperature sensing (4, 5153).

The high throughput of MERFISH measurements allowed us to identify some extremely rare cell types. GnRH-expressing cells (E-30) represent a rare cell population dispersed within the preoptic area and basal forebrain that integrates and orchestrates peripheral and central aspects of reproduction (5456). Only a few cells were found to express appreciable GnRH in our scRNA-seq data and did not form a distinct cluster. By contrast, MERFISH identified a GnRH-enriched cluster (Fig. 6, K and L), which expressed remarkably low levels of Esr1 and Ar (Fig. 6K), suggesting that GnRH neurons may receive indirect feedback from circulating hormones within the hypothalamic–pituitary–gonadal (HPG) axis, potentially through synaptic input from hormone-responsive cell-types such as Kisspeptin (Kiss)–expressing cells (57).

Partition of previously defined cell types into multiple cell populations

Our MERFISH data also partitioned a number of previously reported single cell types—for example, Gal-expressing and Adcyap1-expressing neurons—into multiple distinct cell populations (Fig. 7A). We observed 10 Gal-enriched MERFISH clusters, several of which were scattered across multiple nuclei, such as cluster I-14:MPA/MPN/StHy (Fig. 7, A to C). I-14 was strongly activated during parenting, as shown below, revealing that molecularly and functionally defined cell types can spread across multiple nuclei.

Fig. 7 Subdivisions of neuronal populations expressing Gal or Adcyap1 revealed with MERFISH.

(A) MERFISH subdivides galanin- and Adycap1-expressing cells into multiple transcriptionally and spatially distinct clusters. Color subdivision of the rectangles shows the nuclei spanned by individual clusters, colored identically to the nuclei abbreviations listed on the right. The nuclei abbreviations are as defined in Figs. 3F and 5B. (B) Expression distributions of selected marker genes and genes of interest for all neuronal clusters enriched in galanin (Gal). Expression distributions are calculated as in Fig. 2. (C) Spatial distributions of all inhibitory and excitatory Gal-enriched clusters. (D and E) As in (B) and (C) but for Adcyap1- and Bdnf-enriched clusters. The seven most enriched of the 14 Adcyap1- and Bdnf-enriched clusters are shown. (F) in situ hybridization images of cFos (red), Sncg (green), and overlay of an anterior slice of the preoptic region taken from a heat-stressed animal. The blue boxed region is magnified and shown on the right. Sncg is a marker for the scRNA-seq cluster e13 that corresponds to the MERFISH cluster E-3 (table S9).

Similarly, MERFISH identified 14 clusters enriched in Adcyap1 and Bdnf (Fig. 7, A, D, and E), which were previously designated as markers for warm-sensitive neurons (8). Yet only one of these clusters, E-3:AvPe/Pe/VMPO/VLPO, displayed the established spatial location of warm-sensitive cells (Fig. 7E) (8). Indeed, upon heat stress (29), a high level of the IEG cFos was expressed in cells in the region covered by E-3, and Sncg, a marker gene of E-3’s corresponding scRNA-seq cluster e13 (table S9), was highly enriched in these cFos-positive cells (Fig. 7F). E-3 expressed the leptin and prolactin receptors (Lepr and Prlr) (Fig. 7D), suggesting a mechanism by which metabolic and reproductive states may modulate thermoregulation. Preoptic cells receiving projections from Arcuate Nucleus Kiss1 cells were recently implicated in the regulation of hormonally induced hot flashes via the activation of the receptor Tacr3 (58). E-3 and e13 expressed both Kiss1 receptor and Tacr3, implicating the warm-sensitive cluster in the generation of hot flashes.

Neuronal cell types activated by key social behaviors

To investigate the role of specific neuronal populations in discrete social behaviors, we included cFos in MERFISH measurements and characterized animals after parenting, aggression, or mating. We performed clustering analysis of these behavioral samples together with naïve samples not subjected to behavioral stimuli and did not observe any cluster that was present only in behavioral samples. For each behavior, only a few neuronal clusters, each characterized by key markers, exhibited a statistically significant enrichment in cFos-positive cells (Fig. 8, A and B, and fig. S20). In addition, many other clusters showed a small fraction of cFos-positive cells, which together accounted for a substantial subset of all cFos-positive cells (Fig. 8A and fig. S20B). Activated neurons in all tested behaviors were predominately inhibitory.

Fig. 8 Neuronal clusters activated during specific social behaviors revealed with MERFISH.

(A) Enrichment in cFos-positive cells within each neuronal cluster observed in males or females after displaying a given social behavior. Red bars marked with asterisks are clusters with statistically significant enrichment in cFos-positive cells, as compared with the fraction of cFos-positive cells in all cells (binomial test; false-discovery rate < 5%). Error bars represent standard error of the mean (n = 3 to 5 replicates). We measured fewer slices in behaviorally stimulated animals than in naive animals (4 versus 12 slices per animal) (29), and only clusters in which at least 10 cells are present in two or more behavior replicates are depicted. (B) Expression distributions of selected marker genes and genes of interest for neuronal clusters enriched in cFos-positive cells in the tested social behaviors. Expression distributions are calculated as in Fig. 2. (C) Representative in situ hybridization images of 16-μm-thick sections from the preoptic region showing cFos expression in cells expressing markers of neuronal clusters activated during parenting, in virgin females, mothers, and fathers. Regions in blue dashed boxes are magnified and shown on the right. Red, green, and blue mark the listed genes, and white (or yellow for I-2) indicates coexpression in the merged images. Clusters that cannot be distinguished by a combination of two marker genes plus their spatial location (I-27 and I-10) were not tested. (D) Venn diagrams summarizing the clusters that were activated during specific behaviors in different sexes or physiological states.

We first examined clusters preferentially activated after pup exposure, which elicits parenting behaviors in virgin females, mothers, and fathers but triggers aggression in virgin males (5, 59). We identified, with MERFISH, preferentially activated clusters—clusters enriched in cFos-positive cells (referred to as cFos enrichment hereafter) (Fig. 8, A and B, and fig. S20)—and validated, by use of two- or three-color in situ hybridization, all clusters that could be specifically defined by two marker genes and spatial location (Fig. 8C). In all animal groups that display parenting, I-14:MPA/MPN/StHy had substantially higher cFos enrichment than that of any other cluster (Fig. 8A and fig. S20). The expression of Gal and vasopressin receptor (Avpr1a) within this cluster (Fig. 8B) is consistent with the established role of Gal neurons in parenting (5) and their recently documented vasopressinergic input (60). Moreover, this cluster expressed a large set of hormone and peptide receptors (Fig. 8B), substantiating the complex neuromodulation of parenting (61). In addition to I-14, we observed a preferential activation of clusters by pup exposure in a state-dependent manner. In mothers and fathers, cFos enrichment was observed in Oxtr-expressing cluster I-10:MPA/PaAP/SHy, and mothers additionally showed modest cFos enrichment in I-27:BNST and E-1:AvPe/Pe/MnPO/VMPO (Fig. 8 and fig. S20). Fathers additionally showed cFos enrichment in clusters I-2:BNST/StHy/MPN and in I-16:AvPe/Pe/SHy (Fig. 8 and fig. S20).

By contrast, I-14 was not preferentially activated in virgin males exposed to pups, which is consistent with their aggressive responses toward pups (5). Instead, I-16:AvPe/Pe/SHy and I-24:BNST exhibited cFos enrichment after pup-directed aggression (Fig. 8 and figs. S20 and S21). I-16 was also preferentially activated in virgin males that display inter-male aggression, as was I-2 (Fig. 8 and figs. S20 and S21), suggesting that I-16 is broadly involved in aggressive responses, whereas I-24 and I-2 may mediate differential responses to pups and adults. I-16 was also activated in fathers during parenting, which might indicate that the switch from pup-directed aggression to parenting (59, 61) occurs in circuit nodes downstream of I-16, or that the role of I-16 in aggression is inhibited by pro-parenting circuits. The Gal- and Th-enriched cluster I-16 expressed Vmat2 (Slc18a2) and showed correspondence to the scRNA-seq cluster i16:Gal/Th (table S9), which additionally expressed Ddc, suggesting that this cell population is dopaminergic. The activation of a dopaminergic neuronal population during aggression may provide a cellular basis to understand the observations that dopamine is released during aggression and modulates aggressive behavior (6). I-16 also expressed the opioid receptors Oprd1 and Oprk1 (Fig. 8B), supporting the effects of opioid receptor ligands on aggressive encounters in mice and other rodents (62).

Next, we examined clusters activated by successful mating (29) to capture neural activity associated with appetitive and consummatory aspects of sexual behavior. The cluster I-15:AvPe/Pe/VMPO, which displayed a cell abundance enrichment in female mice compared with males (Fig. 6H), was preferentially activated in females and to a lesser extent in males after mating (Fig. 8 and figs. S20 and S21). Both I-14 and I-15—activated by parenting and mating, respectively—expressed Esr1 (Fig. 8B), which is consistent with recent findings on the involvement of Esr1-expressing cells in parenting and mating behaviors (63, 64). However, our data showed that cells activated by parenting and mating belong to two distinct cell populations localized to distinct preoptic nuclei. More generally, we found Esr1 expression in nearly all behaviorally activated clusters, suggesting that this gene alone cannot define specific behaviorally relevant cell types. In addition to I-15, which was activated in both sexes after mating, we also observed a few clusters that exhibited sexually dimorphic cFos enrichment, such as I-16 in female mating and I-2, I-11, I-14, I-33, E-8, and E-15 in male mating (Fig. 8 and figs. S20 and S21). The weak activation of the parenting cluster I-14 after male sexual behavior is consistent with our previous finding that a small subset of Gal neurons are activated in both mating and parenting (5) and may suggest a mechanism underlying the mating-dependent switch to parental behavior in virgin male mice. Intriguingly, the Th-enriched cluster I-16 was activated by different behaviors in animals of different sexes, mating in females and aggression in males, similarly to the functional sexual dimorphism observed in a recent study of AvPe Th cells (6).

Discussion

Here, we combined the power of scRNA-seq and MERFISH to create a spatially resolved and functionally aware cell atlas of the preoptic region of the mouse hypothalamus. These methods identified major cell classes and neuronal subpopulations with correlated gene expression profiles, providing cross-validations for both methods. Moreover, the two methods are complementary: scRNA-seq measured more genes than MERFISH and helped define marker genes for MERFISH, whereas MERFISH provided spatial context of cells at high resolution as well as more accurate detection and quantification of weakly expressed genes, including functionally important genes such as neuropeptide and hormone receptors. As a result, the combined data provided a more complete picture of the transcriptional diversity and spatial organization of individual cells in the preoptic region.

We observed a remarkable diversity of neurons in this region, comprising ~70 different neuronal clusters. Transcription factors and cell-surface markers have been observed as markers of neuronal identity in the mouse spinal cord (65) and cortex (66) and the Drosophila olfactory system (67). In the mouse preoptic region, genes discriminating neuronal clusters were enriched for neuropeptides and molecules involved in neuromodulator synthesis and transport and for transcription factors. By contrast, neuromodulator receptors were weaker discriminators of these neuronal populations or were expressed at low levels, providing a useful note of caution in using these genes for targeted functional studies in the preoptic region. Many of these neuronal populations were defined by a combination of multiple genes, indicating that genetic intersectional approaches will be most useful in the functional interrogation of specific cell types (67).

The list of cell populations identified in this work substantially expands and further defines previously reported cell types in this region. We observed specific clusters enriched in genes previously identified as markers of functionally important preoptic cells (for example, cells involved in parenting, aggression, thermoregulation, sleep, or thirst), such as Galanin, Th, Adcyap1, Nts, Crh, Tac1, Cck, Agtr1a, Nos1, and aromatase (tables S5 and S9). Our data also resolved many of the previously described singular cell types into multiple cell populations. We observed good correlation between our data and recent scRNA-seq analyses of the whole mouse hypothalamus (68) and of the whole mouse brain (69). In the latter study, which we compared with ours in more detail because it had a larger whole hypothalamic dataset, 15 neuronal clusters were identified from ~2000 profiled hypothalamic neurons (fig. S22) (69). However, because we used ~10 times more neuronal scRNA-seq profiles to characterize about one-fifth of the whole hypothalamus (the preoptic region), we were able to analyze the preoptic region with a greater depth and thus gain finer delineation of cell populations (fig. S22B). The neuronal populations that we uncovered in the preoptic region also largely differed from cell clusters described in scRNA-seq studies of other hypothalamic areas (33, 70), perhaps suggesting a molecular and cellular distinctiveness of this brain area.

MERFISH further allowed us to map the spatial organization of cell populations. Structural features of the hypothalamus are not as visibly apparent as in laminated parts of the brain, and hypothalamic nuclei have largely been defined by subtle differences in neuronal density together with specific connectivity and functional roles (1). However, the differences versus similarities in the cell-type composition and function of distinct nuclei in the hypothalamus remains unclear (71). MERFISH allowed us to examine the organization of distinct cell populations within individual hypothalamic nuclei, providing a framework with which to explore the molecular basis of their anatomical segregation. The spatial organizations of neuronal clusters were diverse: Many of the neuronal clusters were each primarily enriched in one or a few nuclei, whereas several clusters were substantially more dispersed. Moreover, individual nuclei were composed of multiple neuronal clusters. We also observed specific topographical organizations that can support defined modes of function; for example, the spatial proximity of aromatase- and Esr1-expressing cells may support paracrine estrogen signaling. Although aromatase- and Esr1-enriched cell populations regulate sex hormone production and signaling, rarely any of them appeared to be exclusively expressed in either sex, which is consistent with behavioral evidence that males and females are capable of exhibiting behaviors typical of the opposite sex (5, 72).

Last, the ability of MERFISH to interrogate intact tissue allowed us to include activity-dependent IEGs in our measurements, allowing the identification of neuronal populations activated by specific behaviors. Using this approach to study several social behaviors—including parenting, aggression, and mating—we observed that only a small number of neuronal clusters displayed statistically significant enrichment in cFos-positive cells after each behavior. This observation supports a model in which genetically encoded circuits composed of transcriptionally distinct neuronal types control specific hypothalamic functions. However, in all three behaviors, we also observed widespread activation of many neuronal clusters at a substantially weaker level, suggesting a secondary role for many different neuronal types in these behaviors and possibly reflecting necessary cross-talk between different behavior circuits. We caution that the large range of cFos expression levels seen in our samples suggests that some activated clusters with low levels of activity-dependent cFos induction may not have been identified if their cFos levels were below the background noise in our measurements.

This study also extended our previous work on circuits that underlie parenting behavior (5, 60) by resolving preoptic Gal neurons into several distinct subpopulations, with only one of them involved in both male and female parenting. In addition, we identified distinct cell populations that were differentially activated in mothers and fathers during parenting, providing insights into how physiological state may affect parental behavior. Moreover, we identified cell populations associated with sexual behavior in males and females as well as those involved in male aggression toward infants and conspecific males. Together, our data defined functionally relevant cell populations that underlie social behavior with a high molecular and spatial resolution.

Overall, our study demonstrates the power of combining scRNA-seq and MERFISH to map cell types and their organization in the brain, reveal their functional roles in diverse behaviors, and generate hypotheses about structure-function relationships in neural circuits. The identification of marker gene combinations and spatial locations defining the neuronal populations in the preoptic region provides necessary tools for the precise targeting and perturbation of these neurons, thus enabling future functional studies. As an imaging-based approach, we envision that MERFISH can be combined with diverse imaging methods for anatomical tracing and functional interrogation to provide insights into how distinct cell types communicate to form functional circuits in the healthy and diseased brain, as well as in other tissues.

Methods summary

scRNA-seq of the preoptic region was performed by using protocols modified from (70) to increase neuronal survival. Tissue fixation and sectioning as well as MERFISH probe construction, staining, and imaging were performed by using established protocols (26). We imaged 155 genes in MERFISH measurements, with 135 genes imaged by using combinatorial smFISH measurements and 20 additional genes imaged by using sequential rounds of noncombinatorial FISH. The sequences of all probes used for MERFISH are provided in tables S10 and S11. Individual cells were segmented with a seeded watershed algorithm by using DAPI and total mRNA costains (29). Cell clusters were identified by using Louvain community detection on a nearest-neighbor graph built on the statistically significant principle components of gene expression (28, 31, 73) modified to allow an optimized choice of the number of nearest neighbors in the graph. Behavioral stimuli were applied to animals by using established protocols (5, 72, 74), and only animals scored as displaying the desired behavior were used for MERFISH measurements.

Supplementary Materials

www.sciencemag.org/content/362/6416/eaau5324/suppl/DC1

Materials and Methods

Figs. S1 to S22

Tables S1 to S11

References (7798)

References and Notes

  1. Materials and methods are available as supplementary materials.
Acknowledgments: We thank Y. Sigal and H. Babcock for advice on sample preparation and assistance with instrumentation and J. Fan for assistance in the gene set enrichment analysis. Funding: This work is in part supported by the National Institutes of Health (R01MH113094 to X.Z and C.D.; R01MH111502 to C.D., X.Z., and A.R.; and R01HD082131 to C.D.). D.B.-M. and K.S. are supported in part by the National Institutes of Health Pathway to Independence Award (K99HD092542 to D.B-M. and K99EY028625 to K.S.). A.R. was supported by the Klarman Cell Observatory. A.R., C.D., and X.Z. are Howard Hughes Medical Institute Investigators. Author contributions: J.R.M., D.B-M., C.D., and X.Z. conceived the study; J.R.M., D.B-M., and S.W.E. performed experiments with help from J.H. and J.D.P.; J.R.M., D.B.-M., S.W.E., and E.V. analyzed data, with help from K.S. and N.D.R.; K.S. and A.R. provided some of the analysis tools; J.R.M., D.B-M., S.W.E., E.V., C.D., and X.Z. wrote the manuscript, with input from A.R., K.S., J.D.P., N.D.R., and J.H. Competing interests: J.R.M. and X.Z. are inventors on patents applied for by Harvard University related to MERFISH. A.R. is a founder of Celsius Therapeutics and a SAB member of Thermo Fisher Scientific, Driver Group, and Syros Pharmaceuticals. Data and materials availability: scRNA-seq data are available at Gene Expression Omnibus (GEO) (GSE113576). MERFISH data are available at Dryad (76), and the analysis software is available at https://github.com/ZhuangLab/MERFISH_analysis.
View Abstract

Navigate This Article