Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Feb 17;185(4):654-671.e22.
doi: 10.1016/j.cell.2021.12.031. Epub 2022 Jan 21.

A functional cellular framework for sex and estrous cycle-dependent gene expression and behavior

Affiliations

A functional cellular framework for sex and estrous cycle-dependent gene expression and behavior

Joseph R Knoedler et al. Cell. .

Erratum in

Abstract

Sex hormones exert a profound influence on gendered behaviors. How individual sex hormone-responsive neuronal populations regulate diverse sex-typical behaviors is unclear. We performed orthogonal, genetically targeted sequencing of four estrogen receptor 1-expressing (Esr1+) populations and identified 1,415 genes expressed differentially between sexes or estrous states. Unique subsets of these genes were distributed across all 137 transcriptomically defined Esr1+ cell types, including estrous stage-specific ones, that comprise the four populations. We used differentially expressed genes labeling single Esr1+ cell types as entry points to functionally characterize two such cell types, BNSTprTac1/Esr1 and VMHvlCckar/Esr1. We observed that these two cell types, but not the other Esr1+ cell types in these populations, are essential for sex recognition in males and mating in females, respectively. Furthermore, VMHvlCckar/Esr1 cell type projections are distinct from those of other VMHvlEsr1 cell types. Together, projection and functional specialization of dimorphic cell types enables sex hormone-responsive populations to regulate diverse social behaviors.

Keywords: aggression; deep sequencing; estrous cycle; maternal behavior; mating; menstrual cycle; sex differences; social behaviors; synaptic plasticity; transcriptomically defined cell types.

PubMed Disclaimer

Conflict of interest statement

Declaration of interests K.D. is a member of Cell’s Advisory Board. D.W.B., S.I., J.R.K., and N.M.S. are named Inventors on a patent application to be filed on the subject matter.

Figures

Figure 1:
Figure 1:. TRAPseq identification of sex differences in gene expression.
A. Schematic of TRAPseq workflow. B. Scatter plots of sDEGs in different Esr1+ populations. Dots represent DEGs with >1.5-fold change and FDR-adjusted p value < 0.05 (colored) or genes that did not meet both criteria (gray). Colored numbers show number of sDEGs upregulated in that condition and comparison for each region and (far right) for all regions combined. In total, TRAPseq identified 890 and 651 sDEGs in M v FR and M v FNR comparisons. v, versus; TPM, transcripts/million. n = 3/condition. C. ISH for sDEGs in coronal sections through regions enclosed within red dotted areas on Nissl-stained sections on left. ISHs confirm TRAPseq data showing higher expression of Pappa, Auts2, Criml, and Nripl in M and other genes upregulated in FR or FNR mice. Heatmap at far right shows scaled expression of DEGs that were visualized by ISH. Scale = mean z-scored expression of DEGs centered at zero for individual comparisons between conditions (sex or estrous state) and regions; analogous heatmaps provided for ISH studies in all figures. Coronal plane as well as dorsoventral and mediolateral orientations preserved for histological panels in all figures. n = 2/condition/probe. Red, mRNA; blue, DAPI. Scale bars = 100 μm. See also Fig. S1.
Figure 2:
Figure 2:. TRAPseq identification of estrous stage-dependent differences in gene expression.
A. Scatter plots of eDEGs in different Esr1+ populations. Dots represent DEGs with >1.5-fold change and FDR-adjusted p value < 0.05 (colored) or genes that did not meet both criteria (gray). Colored numbers show number of eDEGs upregulated in that condition and comparison for each region and (far right) for all regions combined. In total, TRAPseq identified 770 eDEGs. n = 3. B. ISH for eDEGs. ISHs confirm TRAPseq data showing higher expression of Pgr15l in FNR and other genes upregulated in FR mice. n = 2/condition/probe. Scale bars = 100 μm. C. Heatmap of p values of individual DEGs for the relevant pairwise comparison (with darker green colors indicating more significant p values) illustrating that most DEGs are restricted to one Esr1+ population. D. Heatmap of p values of individual DEGs for the relevant pairwise comparison (with darker green colors indicating more significant p values) illustrating that most DEGs within an Esr1+ population are specific to one comparison between sexes or estrous states. See also Fig. S2.
Figure 3:
Figure 3:. ASD-association and distribution of DEGs in vivo.
A. Heatmap of expression of ASD risk-conferring DEGs. Scale = z-scored expression of DEGs centered at zero for individual comparisons between conditions (sex or estrous state) and regions. See Fig. 1C for heatmap and ISH for Auts2, an ASD-risk conferring DEG. B. ISH for ASD-risk conferring DEGs confirm TRAPseq data showing upregulation of Tcf4, Ptchdl, Sox5, Rfx3, and Ube3a in males. n = 2/condition/probe. Scale bar = 100 μm. C. ISHs for DEGs show that they are distributed in a few (sparse) or most (dense) Esr1+ neurons in each region. Quantification of these ISHs revealed that ‘sparse’ and ‘dense’ DEGs were expressed in ≤10% and ≥70% of Esr1+ cells. n = 2/condition/probe. Scale bar = 100 μm. See also Fig. S3.
Figure 4:
Figure 4:. Categorizing Esr1+ populations into cell types with snRNAseq.
A-D. Violin plots classifying Esr1+ cell types (rows) by virtue of expression of sex hormone receptors, major neuronal neurotransmitter type (excitatory and inhibitory), and enriched genes in the BNSTpr (A), MeA (B), POA (C), and VMHvl (D). CPM, counts per million. n = 10,392 (BNSTpr), 15,929 (MeA), 17,784 (POA), and 5671 (VMHvl) Esr1+ neurons. See also Fig. S4 and Tables S5, S6.
Figure 5:
Figure 5:. Distributive representation of DEGs across Esr1+ cell types.
A. UMAP representation of cell types in Esr1+ populations, with numbers corresponding to those in Fig. 4. B. UMAP representation of Esr1+ cell types by sex or estrous state shows that most of them are proportionally equivalently distributed. C. UMAP representation of excitatory and inhibitory cell types. D. UMAP representation shows that each cell type expresses ≥20% of DEGs in the respective Esr1+ population. E. UMAP representation with scaled distribution of DEGs. Individual cell types express both sDEGs and eDEGs, with different cell types expressing proportionally more DEGs upregulated in a given condition (sex or estrous stage). Gray cell types (scaled value = “0”) express equivalent proportion of DEGs upregulated in both conditions being compared. See also Fig. S5 and Table S2.
Figure 6:
Figure 6:. The BNSTprTac1/Esr1 cell type, among all BNSTprEsr1 cell types, is essential for sex recognition, mating, and aggression in males.
A. Violin plots of a subset of DEGs in the BNSTprTac1/Esr1 cell type. Table S7 lists all DEGs in this cell type. B. Salt and pepper distribution of Tac1 mRNA in Esr1+ BNSTpr neurons. Scale bar = 100 μm. C. Schematic of intersectional chemogenetic strategy to inhibit BNSTprTac1/Esr1 (top) or BNSTprTac1−/Esr1 cell types. D. Inhibition of the BNSTprTac1/Esr1 cell type, but not BNSTprTac1−/Esr1 cell types, abolishes male preference for female urine. E. Inhibition of the BNSTprTac1/Esr1 cell type reduces the probability of resident males initiating mating toward receptive females as well as the number of mount or intromission events/test. F. Inhibition of the BNSTprTac1/Esr1 cell type reduces the probability of resident males attacking intruder males as well as the number of attacks/test. G. Inhibition of BNSTprTac1−/Esr1 cell types does not alter mating of resident males with receptive females. H. Inhibition of BNSTprTac1−/Esr1 cell types does not alter aggression of resident males toward intruder males. Mean ± SEM. n = 2 (B), 8 (D-H). *p < 0.05, **p < 0.01, ***p < 0.001. See also Fig. S6.
Figure 7:
Figure 7:. The VMHvlCckar/Esr1 and VMHvlCckar−/Esr1 cell types are required for female mating and maternal aggression, respectively.
A. Violin plot showing FR-specific expression of Cckar in the VMHvlCckar/Esr1 cell type (#20). B. Representative eDEGs significantly enriched in VMHvlCckar/Esr1 cell type compared to VMHvlCckar−/Esr1 cell types. C. Cckar expression is restricted to the ventrolateral component of the Esr1+ VMHvl population, in agreement with previous work (Hashikawa et al., 2017; Xu et al., 2012). Scale bar = 100 μm. D. Schematic of intersectional chemogenetic strategy to inhibit VMHvlCckar/Esr1 (top) or VMHvlCckar−/Esr1 (bottom) cell types. E. Inhibition of VMHvlCckar/Esr1 cell type in Fr mice significantly reduces lordosis, receptivity index and increases rejection behavior. Lordosis quotient, (lordosis events/# of mounts); Receptivity index (# of intromissions/# of mounts). F. Inhibition of the VMHvlCckar−/Esr1, but not VMHvlCckar/Esr1, cell types abrogates maternal aggression. G. Inhibition of VMHvlCckar−/Esr1 cell types in FR mice does not disable mating behavior. H. Schematic of strategy to label presynaptic termini of VMHvlCckar/Esr1 and Cckar neurons. I. Higher density of mCherry+ termini of VMHvlCckar/Esr1 neurons in AVPV of FR compared to FNR mice. For all findings related to mCherry+ termini in this and subsequent Figures, density of termini is normalized to the number of mCherry+ soma in the VMHvl as described previously (Inoue et al., 2019). J. No change in density of mCherry+ termini of VMHvlCckar− neurons in AVPV between FR and FNR mice. Mean ± SEM. n = 2 (B), 6 (F), ≥7 (G), 8 (H), ≥7 (J), and 6 (K). *p < 0.05. See also Fig. S7 and Table S7.

References

    1. Abrahams BS, Arking DE, Campbell DB, Mefford HC, Morrow EM, Weiss LA, Menashe I, Wadkins T, Banerjee-Basu S, and Packer A (2013). SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol Autism 4, 36. - PMC - PubMed
    1. Alexa A, and Rahnenfuhrer J (2021). topGO: Enrichment Analysis for Gene Ontology.
    1. Arnold AP (2009). The organizational-activational hypothesis as the foundation for a unified theory of sexual differentiation of all mammalian tissues. Horm Behav 55, 570–578. - PMC - PubMed
    1. Asahina K, Watanabe K, Duistermars BJ, Hoopfer E, González CR, Eyjólfsdóttir EA, Perona P, and Anderson DJ (2014). Tachykinin-Expressing Neurons Control Male-Specific Aggressive Arousal in Drosophila. Cell 156, 221–235. - PMC - PubMed
    1. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. (2000). Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 25, 25–29. - PMC - PubMed

Publication types