1. Neuroscience
Download icon

Extensive and spatially variable within-cell-type heterogeneity across the basolateral amygdala

  1. Timothy P O'Leary
  2. Kaitlin E Sullivan
  3. Lihua Wang
  4. Jody Clements
  5. Andrew L Lemire
  6. Mark S Cembrowski  Is a corresponding author
  1. Department of Cellular and Physiological Sciences, Life Sciences Institute, University of British Columbia, Canada
  2. Janelia Research Campus, Howard Hughes Medical Institute, United States
  3. Djavad Mowafaghian Centre for Brain Health, University of British Columbia, Canada
  4. School of Biomedical Engineering, University of British Columbia, Canada
Tools and Resources
Cite this article as: eLife 2020;9:e59003 doi: 10.7554/eLife.59003
9 figures, 1 video, 2 tables, 2 data sets and 1 additional file

Figures

Workflow for assessing heterogeneity within the basolateral amygdala complex.

(A) Atlas schematic of the basolateral amygdala complex (BLA), schematizing microdissection for scRNA-seq. (B) Overview of scRNA-seq data, as visualized through t-SNE dimensionality reduction. Data, along with analysis and visualization tools, available at http://scrnaseq.janelia.org/amygdala. (C) Analysis of coarse discrete heterogeneity within the BLA (see Figures 2 and 3). (D) Analysis of fine continuous heterogeneity within the BLA (see Figures 4 and 5). (E) Spatial registration of discrete and continuous heterogeneity (see Figure 6). (F) Higher order discrete and continuous heterogeneity within the BLA (see Figures 7 and 8).

Figure 2 with 2 supplements
scRNA-seq analysis of coarse, discrete transcriptomic differences.

(A) Overview of scRNA-seq data, as visualized through t-SNE dimensionality reduction and colored according to cluster identity. The number of enriched genes for each cluster is provided. (B) As in (A), but with projections onto the first two principal components. (C) Random forest classification of subsampled data. (D) Expression of known control genes. Expression is colored from low (white) to high (red). Inset numerical values denote maximum CPM value across all cells. (E) As in (D), but with novel cluster-enriched marker genes. (F) Functionally relevant neuronal genes that are differentially expressed between clusters. Numerical values denote CPM values of right tick mark.

Figure 2—figure supplement 1
UMAP and t-SNE dimensionality reduction recapitulate the same general organization.

(A) UMAP visualization of scRNA-seq, with two coarse clusters. (B) t-SNE visualization of scRNA-seq data, with two coarse clusters, included for comparison to UMAP reduction. (C,D) As in (A,B), but for six finer subpopulations.

Figure 2—figure supplement 2
Comparison with previously published BLA scRNA-seq data.

(A) t-SNE visualization of scRNA-seq from Zeisel et al., 2018 excitatory basolateral amygdala neurons (‘TEGLU22’). Example enriched genes listed, which tend to correlate with cell health (e.g. Malat1, mt-Atp6) or with neuronal populations outside of the BLA (e.g. Zbtb20, Strip2; see panels C-Q). (B) As in (A), but for projection onto first two principal components. (C) Expression of Kcng1, a marker restriction to the BLA, in the Allen Mouse Brain Atlas (AMBA) (Lein et al., 2007). (D) Expression of Kcng1 in the TEGLU22 scRNA-seq dataset. Most cells (76%) do not express this BLA marker. (E) Expression of Kcng1 in the dataset of our manuscript. Most cells (86%) express this BLA marker. (F-H) As in (C–E), but for the LA marker Negr1. (I-K) As in (C–E), but for the BA marker Cplx1. L-N. As in (C–E), but for the TEGLU22 marker Strip2, which is expressed outside the BLA. O-Q. As in (L–N), but for Zbtb20.

Figure 3 with 2 supplements
Differentially expressed genes show spatially discrete variation within the BLA.

(A) Three-dimensional coronal volumetric rendering of the BLA in the mouse brain. Image modified from Allen Brain Explorer (Lau et al., 2008). (B,C) As in (A), but for sagittal (B) and horizontal (C) views of the mouse brain. Coronal sections defined in (D) are shown. (D) Coronal sections highlighting the anterior, intermediate, and posterior geometry of the BLA. Atlas definitions of LA vs. BA are shown (modified from Paxinos and Franklin, 2004). (E) Spatial expression of Slc17a7 via chromogenic ISH across the anterior, intermediate, and posterior BLA. Scale bar: 500 μm. Inset shows scRNA-seq expression. (F) Segmentation of Slc17a7-expressing cells from the BLA. (G,H) As in (E,F), but for Negr1 and Ddit4l, novel markers of cluster 1. For reference, locations of segmented cells are superimposed on the location of excitatory (Slc17a7-expressing) neurons within the BLA. (I,J) As in (G,H), but for Cplx1 and Lynx1, novel markers of cluster 2.

Figure 3—figure supplement 1
Spatial organization of inhibitory marker gene expression, and additional cluster-specific marker gene expression.

(A) Representative ISH images for inhibitory neuron marker genes (Gad1, Slc32a1) along the A/P axis of the BLA (left to right). Inserts show t-SNE visualization of scRNA-seq data for each gene with points colored based on expression levels within cells (white to red). (B) Locations of segmented cells for marker genes along the A/P axis of the BLA (anterior, intermediate, posterior). For reference, locations of cells are superimposed on the location of excitatory neurons (Slc17a7-positive) within the BLA. (C,D) As in (A,B), except for an additional pair of marker genes (Cpne8, Gpr88) for the LA cluster (cluster 1). (E,F) As in (A and B), except for except for an additional pair (Rspo2, Cdh9) of marker genes for the BA cluster (cluster 2).

Figure 3—figure supplement 2
Cplx1 expression in the lateral amygdala is associated with inhibitory neurons.

(A) Representative mFISH images of Cplx1 and Gad1 expression in the LA. Scale bar: 200 µm. (B) As in (A), but for an expansion of the LA. Scale bar: 20 µm. (C) Summary of Cplx1-expressing cells that co-expressed Gad1 in the LA. Numbers within each bar denote total number of cells in each category (n = 8 sections from n = 5 animals, cells counted manually).

scRNA-seq analysis of fine, graded transcriptomic differences.

(A) Overview of scRNA-seq data, as visualized through t-SNE dimensionality reduction, at a relatively fine clustering resolution. (B) As in (A), but with projections onto the first two principal components. (C) Random forest classification of subsampled data. For comparison, results from the coarse 2-cluster scheme are also provided in light grey. Inset illustrates rates of misclassification on a per-cell basis, with misclassified cells typically found at the interface of fine clusters. Inset results depict averages across 1000 simulations for 800 training cells. (D) t-SNE illustration of fine-scale clustering within the LA, with number of enriched genes shown for each subpopulation. (E) Expression of example marker genes for the two subpopulations of LA neurons. Note high gene-to-gene variation in expression for marker genes associated with each subpopulation. Inset numerical values denote maximum CPM value across all cells. (F) Expression of functionally relevant genes that are enriched in each LA subpopulation. (G,H,I) As in (D,E), but for subpopulations of the BA.

Graded spatial variation, within the LA and BA, of subpopulation-enriched marker genes.

(A) Representative ISH images for selected LA subpopulation-enriched genes. Scale bar: 500 μm. Insets show t-SNE visualization of scRNA-seq data for each gene. (B) Locations of segmented cells for pairs of genes for each LA subpopulation across anterior, intermediate, and posterior sections. For reference, locations of cells are superimposed on the location of excitatory (Slc17a7-expressing) neurons within the BLA. (C,D) As in (A,B), except for pairs of enriched genes for BA subpopulations.

Figure 6 with 3 supplements
Multiplexed FISH simultaneously maps discrete and continuous heterogeneity.

(A) Overview of mFISH in the anterior amygdala. Scale bar: 400 µm. (B) As in (A), but for intermediate amygdala. Inset provides representative expansion, scale bar: 20 µm. (C) As in (B), but for posterior amygdala. (D) Overview of discrete heterogeneity in the BLA, obtained by assessing relative expression of Negr1 and Cplx1 across cells. (E) Histogram of phenotype identities for (D). (F) Illustration of gene-expression correlations for all cells, relative to the Negr1 phenotype seed cell highlighted in magenta. (G) Correlations computed and averaged across space, relative to the magenta seed cell shown in (F). Magenta and green points represent other cells within the posterior BLA with a binarized discrete cell-type identity. Blue and orange lines indicate binned average correlation across all cells at 100 µm spacing for the same and opposite phenotypes. (H,I) As in (F,G), but for a Cplx1 phenotype seed cell. (J) As in (G,I), but pooling and analyzing averages across all cells in the BLA. Blue and orange lines illustrate mean correlations to other cells of the same and opposite phenotypes, respectively, with spread about mean indicating one standard deviation.

Figure 6—figure supplement 1
Individual mFISH marker-gene expression in the BLA.

(A) Overview of smFISH signals for the 12 BLA marker genes used in this study. Scale bar: 300 μm. (B) Expanded views of marker gene expression, for the respective boxed areas shown in (A). Scale bar: 40 μm.

Figure 6—figure supplement 2
Clustering of mFISH excitatory neurons.

(A) Points denoting individual excitatory neurons in mFISH data, colored according to respective coronal section. (B) t-SNE visualization of mFISH data for the cells shown in (A). (C) Spatial location of cells comprising two clusters. (D) t-SNE visualization for two excitatory neuron clusters, colored according to cluster identity from (C). (E,F) As in (C,D), but for six clusters.

Figure 6—figure supplement 3
Interneuron subtypes resolved by scRNA-seq and chromogenic ISH.

(A) t-SNE visualization of 51 interneurons, colored according to cluster identity (Int1 and Int2). Inset numeric values denote number of enriched genes in each cluster, relative to the other interneuron cluster. Inset t-SNE visualization shows expression of interneuron marker Gad1. (B) Examples of novel interneuron cluster-enriched genes obtained from scRNA-seq. (C) Summary of all genes differentially expressed between the two interneuron clusters. Functional categories include calcium handling (Hpca, Wfs1), synapses (Gabra1, Cplx1, Syt5), and regulation of transcription (Tshz1). (D) Chromogenic ISH images for differentially expressed genes. Genes enriched in Int1 show diffuse spatial patterning across the BLA and minimal labeling in the intercalated region. Genes enriched in Int2, conversely, show dense expression in the capsular region. Scale bar: 300 μm.

Figure 7 with 1 supplement
BLA projection neurons show spatially discrete and continuous cellular variation.

(A) Atlas schematic showing retrograde viral injection site for the ACX (modified from Paxinos and Franklin, 2004). (B) Representative maximum intensity projections of labeled cells within the BLA. Scale bar: 500 μm. The boundary of the BLA is shown (solid line) along with the division between the lateral and basal sub-regions (dashed line) (C). Locations of all virally labeled cells along the A/P axis of the BLA (left to right). For reference, virally labeled cells (black) are superimposed over the position of excitatory neurons (Slc17a7-expressing cells), with cells colored to illustrate LA (magenta) and BA (green) per atlas definitions. (D) Percentage of total labeled cells located within the LA and BA, depicted as mean ± SEM. Grey lines indicate range of values within 95% confidence intervals for the expected number of cells in the LA and BA, via Monte Carlo simulations of random selection of excitatory neurons. E-X. As in (A–D), except for the NAC (E–H), PFC (I–L), VHC (M–P), MEC (Q–T) and RSC (U–X) injection sites.

Figure 7—figure supplement 1
Spatial cellular gradients of projection neurons are present in the BLA and differ based on projection site.

(A) Illustration of the BLA showing the LA (magenta) and sub-regional parcellation of the BA (green; B-M: basal medial; B-L: basal lateral). (B) Percentage of virally labeled NAC projection neurons located within the lateral BLA, and sub-regions of the BA. Individual points and connecting lines represent values from single animals. Bars depict means ± SEM. Solid black line indicates percentage of excitatory cells (Slc17a7 positive) expected in each region. (C-F) As in (B), except for PFC, VHC, MEC and RSC projection neurons. (G) As in (A), except showing parcellation of the LA (L-M: lateral medial; L-L: lateral lateral). (H) Percentage of ACX projection cells located within the BA and sub-regions of the LA.

LA and BA neurons vary by morphology.

(A) Representative images of labeled cells projecting to the ACX, NAC, PFC, VHC, MEC, and RSC. Scale bars: 25 μm. (B) The cell-body area of virally labeled cells for each projection site, depicted as mean ± SEM. (C) Top row: overview of mFISH signals for Slc17a7, Negr1, and Cplx1 in the intermediate BLA. Scale bars: 300 μm. Middle, bottom rows: expansions for the LA and BA, respectively. Scale bars: 40 μm. (D) Left: Cell-body area of Slc17a7-expressing cells segmented from (C). Middle, right: Cell-body area of Slc17a7-expressing cells that also express Negr1 or Cplx1, respectively. (E) Cell-body area for Slc17a7-expressing cells that also express Negr1 (magenta) or Cplx1 (green), depicted as mean ± SEM.

Summary of heterogeneity within the basolateral amygdala.

Examples of spatially variable gene-expression, circuit, and cellular properties of excitatory neurons in the basolateral amygdala, as resolved in this study.

Videos

Video 1
Overview of mFISH marker gene expression in the posterior BLA.

Tables

Table 1
List of image series used from the Allen Mouse Brain Atlas.
GeneFigureImage series
Kcng12-S2C77340480
Negr12-S2F, 3D692
Cplx12-S2I, 3F, 6-S1D67752308
Strip22-S2L72283809
Zbtb202-S2O79568020
Slc17a73B70436317
Ddit4l3D71836878
Lynx13F655
Gad13-S1B79556706
Slc32a13-S1B72081554
Cpne83-S1E73520974
Gpr883-S1E79567811
Rspo23-S1H71016632
Cdh93-S1H72472764
Rorb5B79556597
Myl45B72129251
Cdh135B79490066
Otof5B73788043
Prss235E70634118
S100b5E79591593
Nnat5E77887874
Slc24A25E71924238
Adamts25E71924385
Etv15E72119595
Bdnf5E79587720
Nr4a25E732
Slc6a16-S1D79591685
Tshz16-S1D72129289
Meis26-S1D1231
Key resources table
Reagent type
(species) or
resource
DesignationSource or
reference
IdentifiersAdditional
information
Sequence-based reagentCplx1 ISH probeAdvanced Cell Diagnostics482531-C3Two-color FISH
Sequence-based reagentGad1 ISH probeAdvanced Cell Diagnostics400951-C2Two-color FISH
Sequence-based reagentNr4a2 ISH probeAdvanced Cell Diagnostics423351-T1Multiplexed FISH
Sequence-based reagentOtof ISH probeAdvanced Cell Diagnostics485671-T2Multiplexed FISH
Sequence-based reagentCdh13 ISH probeAdvanced Cell Diagnostics443251-T3Multiplexed FISH
Sequence-based reagentRorb ISH probeAdvanced Cell Diagnostics444271-T4Multiplexed FISH
Sequence-based reagentAdamts2 ISH probeAdvanced Cell Diagnostics806371-T5Multiplexed FISH
Sequence-based reagentPrss23 ISH probeAdvanced Cell
Diagnostics
447921-T6Multiplexed FISH
Sequence-based reagentBdnf ISH probeAdvanced Cell Diagnostics424821-T7Multiplexed FISH
Sequence-based reagentSlc5a5 ISH probeAdvanced Cell Diagnostics487721-T8Multiplexed FISH
Sequence-based reagentSlc17a7 ISH probeAdvanced Cell Diagnostics416631-T9Multiplexed FISH
Sequence-based reagentNnat ISH probeAdvanced Cell Diagnostics432631-T10Multiplexed FISH
Sequence-based reagentNegr1 ISH probeAdvanced Cell Diagnostics806361-T11Multiplexed FISH
Sequence-based reagentCplx1 ISH probeAdvanced Cell Diagnostics482531-T12Multiplexed FISH
Software,algorithmRhttps://www.r-project.orgSCR_001905-
Software,algorithmSeurathttps://satijalab.org/seurat/RRID:SCR_007322-
Software,algorithmFijihttps://imagej.net/FijiRRID:SCR_002285-
Software, algorithmPrismhttps://www.graphpad.com/scientific-software/prism/RRID:SCR_002798-
OtherAAV-SL1-CAG-tdTJanelia Viral Core--
OtherAAV-SL1-CAG-GFPJanelia Viral Core--
Strain, strain background Mus musculusPenk-creJanelia Research Campus--

Data availability

Raw and processed scRNA-seq datasets have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus under GEO: GSE148866. Data underlying the results of this manuscript have been provided on FigShare (https://doi.org/10.6084/m9.figshare.c.5108165).

The following data sets were generated
  1. 1
    NCBI Gene Expression Omnibus
    1. M Cembrowski
    (2020)
    ID GSE148866. Discrete and continuous heterogeneity in the basolateral amygdala.
  2. 2

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)