Developmental single-cell transcriptomics of hypothalamic POMC neurons reveal the genetic trajectories of multiple neuropeptidergic phenotypes

  1. Hui Yu  Is a corresponding author
  2. Marcelo Rubinstein  Is a corresponding author
  3. Malcolm J Low  Is a corresponding author
  1. Department of Molecular and Integrative Physiology, University of Michigan Medical School, United States
  2. Instituto de Investigaciones en Ingeniería Genética y Biología Molecular, Consejo Nacional de Investigaciones Científicas y Técnicas, Argentina
  3. Departamento de Fisiología, Biología Molecular y Celular, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Argentina

Abstract

Proopiomelanocortin (POMC) neurons of the hypothalamic arcuate nucleus are essential to regulate food intake and energy balance. However, the ontogenetic transcriptional programs that specify the identity and functioning of these neurons are poorly understood. Here, we use single-cell RNA-sequencing (scRNA-seq) to define the transcriptomes characterizing Pomc-expressing cells in the developing hypothalamus and translating ribosome affinity purification with RNA-sequencing (TRAP-seq) to analyze the subsequent translatomes of mature POMC neurons. Our data showed that Pomc-expressing neurons give rise to multiple developmental pathways expressing different levels of Pomc and unique combinations of transcription factors. The predominant cluster, featured by high levels of Pomc and Prdm12 transcripts, represents the canonical arcuate POMC neurons. Additional cell clusters expressing medium or low levels of Pomc mature into different neuronal phenotypes featured by distinct sets of transcription factors, neuropeptides, processing enzymes, cell surface, and nuclear receptors. We conclude that the genetic programs specifying the identity and differentiation of arcuate POMC neurons are diverse and generate a heterogeneous repertoire of neuronal phenotypes early in development that continue to mature postnatally.

Editor's evaluation

This study investigates the developmental origins of functionally distinct neuronal populations in the arcuate nucleus of the hypothalamus. This study reports the transcriptional program of develop of multiple subtypes of POMC neurons, which are important for multiple physiological and behavioral functions, including appetite. The paper uses cutting edge methods will be important for research studying the arcuate nucleus including regulating food intake and metabolism.

https://doi.org/10.7554/eLife.72883.sa0

Introduction

Proopiomelanocortin (POMC)-expressing neurons located in the arcuate nucleus of the hypothalamus play an essential role in the regulation of food intake by maintaining a melanocortin-dependent anorexigenic tone. The critical importance of hypothalamic melanocortins is apparent in mice lacking Pomc that develop hyperphagia and early-onset extreme obesity (Bumaschny et al., 2012). This phenotype is closely mirrored in humans carrying biallelic null mutations in POMC (Krude et al., 1998). In addition, arcuate POMC neurons release the Pomc-encoded opioid peptide β-endorphin, which is critical in mediating normal food intake (Appleyard et al., 2003) and stress-induced analgesia (Rubinstein et al., 1996).

Postmitotic POMC neurons acquire their phenotypic identity in the developing hypothalamus days earlier than any other peptidergic neurons born in the presumptive arcuate nucleus (e.g. NPY/AGRP, GHRH, SST, KISSPEPTIN). Pomc mRNA is initially observed in the tuberal portion of the prospective mouse hypothalamus of E10.5 embryos, following the concurrent advent of the essential transcription factors (TFs) ISL1 (Nasif et al., 2015), NKX2-1 (Orquera et al., 2019), and PRDM12 (Hael et al., 2020) that act together in a combinatorial manner. In fact, the early ablation of either Isl1, Nkx2-1, or Prdm12 disrupts the onset of hypothalamic Pomc expression in conditional mutant mouse embryos (Nasif et al., 2015; Orquera et al., 2019; Hael et al., 2020). Given that the number of arcuate neurons co-expressing Nkx2-1, Isl1, and Prdm12 greatly exceeds that of POMC neurons in this area, it is expected that additional and still unknown TFs integrate a core regulatory complex completing the genetic program that specifies the identity of POMC neurons.

Other molecular markers of arcuate POMC neurons such as transporters, receptors, channels, and co-transmitters exhibit great variability indicating that a heterogeneous pool of diverse subpopulations of these neurons are involved in multiple circuits and functions. In fact, previous studies have suggested that subsets of arcuate Pomc-expressing neurons may act as precursors of terminally differentiated neurons of alternative neuropeptidergic phenotypes such as NPY/AGRP and KISSPEPTIN/NEUROKININ-B/DYNORPHIN (KNDY) neurons in the adult hypothalamus (Padilla et al., 2010; Sanz et al., 2015). Initial attempts to characterize the transcriptome of adult POMC neurons at the single-cell level confirmed the heterogeneous nature of arcuate POMC neurons and revealed the existence of distinct clusters including those co-expressing the aforementioned neuropeptide genes (Campbell et al., 2017; Chen et al., 2017; Lam et al., 2017; Huisman et al., 2020). However, the molecular signatures of the various hypothalamic Pomc-expressing lineages are still lacking.

Here, we track the origin and maturation of arcuate POMC neurons in the developing and early postnatal hypothalamus by performing single-cell RNA-sequencing (scRNA-seq) transcriptomics of fluorescently labeled cells taken from Pomc-TdDiscomaRed-Sv40PolyA (Pomc-TdDsRed) transgenic mice at different embryonic (E11.5, E13.5, E15.5, and E17.5) and early postnatal (P5 and P12) ages. Our data indicate that POMC neurons give rise to distinct transcriptional trajectories and show that a heterogeneous population of Pomc-expressing neurons is established early in the developing arcuate nucleus. Combining translating ribosome affinity purification with RNA-sequencing (TRAP-seq) of Pomc-eGFPL10a transgenic mice at P12 and P60, we identified transcriptional programs that emerge from distinct populations of early embryonic POMC neurons that later mature into fully functional POMC neurons. We believe that these datasets provide a valuable resource for comprehensively understanding the genetic bases underlying the ontogeny and development of POMC neurons as well as for the rational design of further functional studies of brain circuits involving POMC neurons.

Results

Cluster analysis of gene expression in POMC positive hypothalamic cells integrated across six development ages

We captured the single-cell transcriptomes of hypothalamic Pomc-expressing cells at four embryonic days (E11.5, E13.5, E15.5, and E17.5) and two early postnatal days (P5 and P12) when the maturing neurons already develop axons and innervate key target nuclei distal to the Arc (Bouret, 2004a, Bouret et al., 2004b, Bouret et al., 2012). Medial basal hypothalami from mice that express the reporter transgene Pomc-TdDsRed selectively in POMC cells were pooled from multiple embryos or postnatal pups at the indicated ages. Dissociated single-cell suspensions were sorted for DsRed fluorescence and used for library preparation with the 10× Genomics platform (Figure 1A). After a rigorous screening procedure to remove substandard scRNA-seq data (see Materials and methods), we used a combination of computational platforms to analyze the transcriptomes of 13,962 cells that had Pomc counts of ≥1 unique molecular identifiers (UMIs).

Figure 1 with 3 supplements see all
Single-cell sequencing of PomcDsRed cells from developing mouse hypothalami at E11.5, E13.5, E15.5, E17.5, P5, and P12 reveals eight distinct neuronal clusters and three non-neuronal clusters.

(A) Schematic diagram showing the overall experimental procedure. (B) UMAP (Uniform Manifold Approximation and Projection) plot showing the distribution of major cell clusters and their corresponding numbers of cells at each developmental stage. (C) Violin plots showing the expression of three signature genes and Pomc for each cluster. (D) UMAP plots showing the enrichment of feature genes in each cluster. Genes representing each cluster were colored and highlighted to show the cluster-specific enrichment. The color intensity corresponds to the normalized gene expression from Seurat analysis, where gene expression for each cell was normalized by the total expression, multiplied by a scale factor of 10,000, and then log-transformed. VLMC, vascular leptomeningeal cell.

Figure 1—source data 1

Average gene expression obtained from the Seurat analysis of each of the 11 proopiomelanocortin (Pomc) clusters integrated across all six developmental stages (corresponding to Figure 1).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig1-data1-v2.xlsx
Figure 1—source data 2

Feature genes that define each of the 11 proopiomelanocortin (Pomc) clusters (corresponding to Figure 1) integrated across all six developmental stages.

https://cdn.elifesciences.org/articles/72883/elife-72883-fig1-data2-v2.xlsx

A Seurat analysis of all cells integrated across the six time points revealed 11 distinct cell clusters projected onto a UMAP plot (Figure 1B). The same data are represented in a heatmap corresponding to the top 30 expressed genes per cluster (Figure 1—figure supplement 1). Clusters 1–8 (10,828 cells) were all highly enriched for a panel of three neuronal marker gene transcripts. The embedded table in Figure 1B lists the number of neurons of each cluster at each developmental age. The other three clusters (3134 cells) were enriched in transcripts characteristic of either microglia, vascular leptomeningeal cells (VLMC), or astrocytes (Figure 1C). Gene ontology analyses for biological processes further confirmed the neuronal or non-neuronal identities of the 11 individual clusters (data not shown).

By definition, all cells in the neuronal and non-neuronal clusters contained Pomc transcripts (UMI ≥1), however there were significant differences in the average expression levels per cluster, ranging from 1.3 to 33.8 [exp(log Normalized counts) –1]. Nomenclature of the eight neuronal clusters is based on a combination of relatively high (>20), medium (5–20), or low (<5) Pomc levels and high expression of one differential key feature gene as follows: no. 1, Pomc(high)/Prdm12; no. 2 Pomc(med)/Ebf1; no. 3, Pomc(med)/Nr4a2; no. 4, Pomc(low)/Otp; no. 5, Pomc(low)/Tac2; no. 6, Pomc(low)/Dlx5; no. 7, Pomc(low)/Prdm13; and no. 8, Pomc(low)/Pitx2. The cluster distribution and relative expression of those feature genes are represented by individual UMAP plots (Figure 1D). Average expression levels of all genes, including Pomc, in the 11 clusters and complete lists of all feature genes that define the individual clusters are presented in Figure 1—source data 1 and Figure 1—source data 2, respectively.

In order to validate our experimental approach of isolating Pomc-expressing cells from the developing hypothalamus by fluorescence-activated cell sorting (FACS) for the transgenic surrogate marker Pomc-TdDsRed, we performed an independent unsupervised cluster analysis of the transcriptomes from the 5724 cells that had DsRed counts of ≥1 UMI (Figure 1—figure supplement 2A-C and ). A projection analysis (Figure 1—figure supplement 2D-F) of the differentially expressed feature genes defined nine DsRed clusterswhich closely corresponded to 9 of the 11 clusters based on Pomc counts of ≥1 UMI. The two exceptions were the smallest clusters 7. Pomc(low)/Prdm13 and astrocytes (Figure 1—figure supplement 3D-F). Furthermore, the number of cells at each developmental age from the corresponding neuronal clusters was remarkably similar in the Pomc UMI ≥ 1 dataset compared to the DsRed UMI ≥ 1 dataset (compare Figure 1B to Figure 1—figure supplement 3A). The different population of DsRed fluorescent sorted cells showing UMI ≥ 1 for Pomc (13,962) or DsRed (5724) may be due to differential half-lives of the fluorescent protein DSRED used during cell sorting and its transcript used for barcoding the RNA-seq libraries.

Temporal gene expression analysis of the eight neuronal clusters at individual developmental stages

Next, we analyzed temporal gene expression patterns for each of the eight neuronal clusters at the six distinct developmental ages (Figure 2—source data 1). The progressive changes in feature genes across time are responsible for the observed within cluster heterogeneity on the heatmap in Figure 1—figure supplement 1.

Cluster 1

Pomc(high)/Prdm12 cells were the most abundant among the eight neuronal clusters and consistently had the highest levels of Pomc expression at all six developmental ages. We consider them to be the cell lineage that matures into the canonical POMC neurons in the adult hypothalamus (Figure 2A–C). Other feature genes included a set of four TFs known to be critical for the initiation and/or maintenance of Arc Pomc expression, Isl1, Nkx2-1, Prdm12, and Tbx3 (Nasif et al., 2015; Orquera et al., 2019; Quarta et al., 2019, Hael et al., 2020; Figure 2A). They were expressed in similar developmental patterns to each other, with the exception of Tbx3 that was significantly expressed only from day E15.5. Notably, there was a lack of the pituitary-specific Tbx19 TF, confirming that no pituitary tissue was included in the hypothalamic dissections. Two homologues of the Drosophila sine oculus homeobox gene, Six3 and Six6, which have not been described previously in the context of Arc POMC neuron development, were also expressed robustly at all six time points (Figure 2A). Co-expression of their mRNAs with Pomc by RNAScope at ages E12.5 and E16.5 confirmed the scRNA-seq analysis, although the two TFs were expressed in additional hypothalamic domains and the pituitary gland (Figure 2B). Vim1, Sox1, -2, -3, and -14, and Myc, whose expression is typically associated with neuronal precursors or immature neurons, tended to have higher expression at the earlier embryonic ages (Figure 2A). The nuclear receptor gene Nr5a1, also known as steroidogenic factor 1 (Sf1), was a key feature gene expressed almost exclusively at age E11.5 (Figure 2A), but only minimally expressed at later embryonic time points when NR5A1 is critically involved in the development of the hypothalamic ventromedial nucleus (Ikeda et al., 1995). This differential spatio-temporal presence of NR5A1 within or outside POMC neurons was confirmed by immunofluorescence on hypothalamic sections taken from E11.5 and E17.5 mouse embryos (Figure 2B). Interestingly, Nr5a1, Nr0b1 (Dax1), and Shh transcripts were also present in cluster 1 at the two earliest ages, and are factors known to interact with the Wnt/β-catenin signaling pathway in early cell fate determination in the adrenal and pituitary glands (Luo et al., 1994; Meeks et al., 2003; Mizusaki et al., 2003). A combination of pseudotime and gene ontology analyses for the top expressed genes (Log2FoldChange > 0.25, p < 0.05) in cluster 1 revealed that the genes expressed at the beginning of pseudotime function in the organization of cell projections, axonogenesis, and the regulation of cellular component size, whereas genes expressed at the middle and end of the pseudotime function are related to neuron maturation and neuron-specified functions such as the formation of myelin sheath and regulation of synaptic vesicle cycling (Figure 2C).

Figure 2 with 2 supplements see all
Temporal gene expression patterns of the Pomc(high)/Prdm12 cluster and Pomc(med)/Ebf1 clusters.

(A and D) Violin plots showing the expression of signature genes at each developmental stage in the Pomc(high)/Prdm12 cluster and Pomc(med)/Ebf1 cluster, respectively. (B) Fluorescence in situ hybridization showing the co-localization of Pomc (green) and Six6 (red), the co-localization of Pomc (green) and Six3 (red) at indicated developmental stages. Immunofluorescence showing the co-localization TDTOMATO (red) NR5A1 (green) at the indicated developmental stages. (C and F) Heatmaps showing the gene ontology analysis based on the top marker genes in the order of pseudotime in Pomc(high)/Prdm12 cluster and Pomc(med)/Ebf1 cluster, respectively. (E) Fluorescence in situ hybridization showing the co-localization of Pomc (green) and Ebf1 (red) at E11.5. (G) Reclustering cells from the Pomc(med)/Ebf1 cluster reveals two subclusters denoted as 2–1 Pomc(high) and 2–2 Pomc(low). Violin plots showing the signature gene expression in 2–1 Pomc(high) and 2–2 Pomc(low) clusters. Insets are magnified views of the indicated boxes. Arrows indicate co-expressing neurons in the merged panels. Scale bar: 50 µm. Image orientation: left, posterior(P); right, anterior(A). RP: Rathke’s pouch; Pit: pituitary gland. Arc: arcuate nucleus; VMH: ventromedial hypothalamus; V: ventricular zone; Mes: mesenchyme.

Figure 2—source data 1

Average gene expression obtained from the Seurat analysis of six developmental stages of each cluster (corresponding to Figure 2, Figure 3, and their linked Figure 2—figure supplement 2).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig2-data1-v2.xlsx
Figure 2—source data 2

Average gene expression obtained from the Seurat analysis of the two subclusters 3–1 and 3–2 derived from reclustering of cluster 3.

Pomc(med)/Nr4a2 and the two subclusters 2–1 and 2–2 derived from reclustering of cluster 2. Pomc(med)/Ebf1 (corresponding to Figure 2G, cluster 2 and Figure 2—figure supplement 1, cluster 3).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig2-data2-v2.xlsx
Figure 2—source data 3

Feature genes that define each of the two subclusters 3–1 and 3–2 derived from reclustering of cluster 3.

Pomc(med)/Nr4a2 and the two subclusters 2–1 and 2–2 derived from reclustering of cluster 2. Pomc(med)/Ebf1 Ebf1 (corresponding to Figure 2G, cluster 2 and Figure 2—figure supplement 1, cluster 3).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig2-data3-v2.xlsx

Cluster 2

Pomc(med)/Ebf1 neurons made up the second highest abundance cluster of POMC neurons, particularly at E11.5 and E13.5. After these time points, the fraction of all neurons constituting cluster 2 dropped by approximately two-thirds from ages E15.5 through P12. In contrast, the number of neurons in clusters 3–8 increased over time with peaks at different embryonic ages (Figure 1B). Cluster 2 neurons, on average, had medium levels of Pomc transcripts that were divided distinctly from their lowest at ages E11.5–E13.5 to highest at later embryonic and postnatal ages (Figure 2D; Figure 2—source data 1). Major cluster 2 feature markers included members of the non-basic helix-loop-helix early B-cell TF family EBF (Figure 1—source data 2). Ebf1 is typically co-expressed with either Ebf2 or Ebf3 and pairs of these factors were shared not only in cluster 2, but also in clusters 3 and 8 (Figure 2—figure supplement 2). Dual RNAScope in situ hybridization showed overlapping co-expression of Ebf1 and Pomc mRNAs in the ventricular zone of the developing hypothalamus at an early embryonic stage (E11.5) (Figure 2E). However, transcript levels of Ebf1 are much lower than what was found in other areas of the same sections, such as the mesenchyme (MES) (Figure 2E). Ebf1 mRNA was also detected in Pomc-expressing cells at E13.5, and then was barely detectable at E17.5 (Figure 2—figure supplement 1A and B). A similar expression pattern was found with another cluster 2 featured gene, Pou4f1 (POU Class 4 Homeobox 1). Dual RNAScope in situ hybridization confirmed the overlap of Pou4f1 and Pomc mRNAs at E13.5, but not at E17.5 (Figure 2—figure supplement 1C and D). These results, together, suggest that Pomc mRNA levels peak only after expression of both Ebf1 and Pou4f1 decrease. Two other prominent feature genes in cluster 2 were Crabp1 and Crabp2, which encode cellular retinoic acid binding proteins. A majority of the marker genes for cluster 2 are predicted to be expressed at the beginning of pseudotime and the main functions are associated with the maintenance of basic cell biological processes such as nucleosome organization, axonogenesis, microtubule transport, and intracellular transport (Figure 2F). The genes ordered in the middle and end of pseudotime function for the regulation of synapse structure and activities.

Visual examination of the UMAP plots in Figure 1D suggested that cluster 2 was composed of two compartments, one with high and one with low Pomc abundance. Unsupervised reclustering of the differentially expressed genes from integrated cluster 2 confirmed the existence of two subclusters (Figure 2G). Subcluster 2–1 contained two-thirds of the neurons, which were POMC(high) whereas subcluster 2–2 contained the other one-third of cluster 2 neurons, which in contrast, were Pomc(low). As shown in Figure 2D, subcluster 2–2 contains cells mainly from the early developmental stages E11.5 to E15.5, whereas, subcluster 2–1 contains cells from E17.5 to P12. Consistently, known activating TFs for Pomc expression such as Nkx2-1, Isl1, and Tbx3 were limited to subcluster 2–1 while relatively high levels of Ebf and Crabp transcripts were characteristic of subcluster 2–2 (Figure 2—source data 2 and Figure 2—source data 3). Another distinctive feature of cells in cluster 2 (both subclusters) is the absence of Prdm12 transcripts, one main component of cluster 1 (Figure 2A and D).

Cluster 3

Pomc(med)/Nr4a2 neurons displayed high expression of its major feature gene Nr4a2 (Nuclear Receptor Subfamily 4 Group A Member 2) (Figure 2—figure supplement 2). Nr4a2 mRNA was also expressed at relatively high levels in clusters 2. Pomc(med)/Ebf1 and 8. Pomc(low)/Pitx2. Nr4a2 (Nurr1) has been implicated in the regulation of Pomc expression in the pituitary corticotropic cell line AtT20 (Kovalovsky et al., 2002). Furthermore, Nr4a2 is a pioneer transcriptional regulator important for the differentiation and maintenance of meso-diencephalic dopaminergic (mdDA) neurons during development (Perlmann and Wallén-Mackenzie, 2004). It is crucial for expression of a set of genes including Slc6a3 (plasma membrane dopamine transporter, DAT), Slc18a2 (synaptic vesicular monoamine transporter, VMAT), Th (tyrosine hydroxylase), Foxa1 and Foxa2 (forkhead family TFs), and Drd2 (dopamine receptor 2) that together define the dopaminergic neuronal phenotype (Hong et al., 2014; Jankovic et al., 2005; Saucedo-Cardenas et al., 1998; Smits et al., 2003). An additional key factor involved in the generation and maintenance of mdDA neurons, Lmx1a (Lim- and homeodomain factor 1a) (Andersson et al., 2006; Hong et al., 2014), was also a key feature gene of clusters 3 and 8 (Figure 2—figure supplement 2A and D) suggesting that POMC and DA neurons share at least part of their genetic differentiation programs.

Similar to cluster 2, reclustering of the differentially expressed genes from cluster 3 revealed two subclusters (Figure 2—figure supplement 2C). Subcluster 3–1 constituted 72% of the main cluster 3 and was characterized by Pomc(low) neurons with strong Nr4a2 expression at all developmental ages. Feature genes for subcluster 3–1 also included Lmx1a, Slc6a3, Slc18a2, Foxa2, and Th, all feature genes for developmental subclusters E13.5–10 and E15.5–11 derived from cluster 3 indicating that a fraction of cluster 3 neurons were dopaminergic and existed transiently during a 3-day developmental window (Figure 1—source data 2, Figure 2—source data 1). In contrast, subcluster 3–2 neurons were Pomc(high) together with characteristic expression of the known Pomc transcriptional activator genes and were present only at ages P5 and P12. The dissociation of high Pomc expression with low Nr4a2 expression in subcluster 3–2 suggests that the latter is not necessary for hypothalamic Pomc expression, unlike its role in pituitary corticotrophs (Murphy and Conneely, 1997).

Cluster 8

Pomc(low)/Pitx2 is more closely related to clusters 2. Pomc(med)/Ebf1 and 3. Pomc(med)/Nr4a2 than to any of the others (Figure 2—figure supplement 2D). In addition to Pitx2, major feature genes were Barhl1, Foxa1, Foxp2, Lmx1a/b, and Tac1 (Figure 1—source data 2). Intriguingly, virtually all of the top cluster 8 defining genes (Figure 1—source data 2) are identical to those that characterize both the mesencephalic dopamine neurons and glutamatergic, non-dopaminergic, subthalamic nucleus neurons (Kee et al., 2017). Dual RNAscope in situ hybridization for Pitx2 and Pomc showed limited overlap at E12.5 (Figure 2—figure supplement 2E). Cluster 8 is the only cluster with the neuropeptide Cck as a feature gene at developmental ages P5 and P12 (Figure 1—source data 1 and Figure 2—source data 1).

Cluster 4

Pomc(low)/Otp cells uniquely expressed Otp at all developmental ages (Figure 3A). Pomc expression was moderate at E11.5 and then dropped to consistently lower levels through age P12 despite continued expression of Isl1, Nkx2-1, and Tbx3. Interestingly, expression of the other transcriptional modulator important for Pomc expression, Prdm12, was detected only at E11.5. Other characteristic features of cluster 4 were increasing gradients of Agrp, Npy, and Calcr expression over the course of development (Figure 3A). Minimal co-expression of Pomc and Npy was confirmed by fluorescence in situ hybridization (FISH) at age E15.5 (Figure 3B). Sst was highly expressed in this cluster at ages E13.5, P5, and P12 with a temporary drop off at E15.5 and E17.5. Altogether, these patterns of gene expression indicate that a subpopulation of POMC neurons expressing Otp further differentiates into AGRP/NPY neurons as suggested previously by studies based on transgenic lineage tracing (Padilla et al., 2010). Most of the top enriched genes of this cluster are expressed in the middle and end portions of pseudotime prediction highlighting functions such as regulation of inhibitory synapses, ribosomal biogenesis, and synaptic vesicle cytoskeletal transport (Figure 3C). A map of gene expression in individual cells within the cluster 4 UMAP plot demonstrated that Npy was expressed in virtually all Agrp neurons, but few Sst neurons. Sst and Agrp neurons were essentially separate cell populations (Figure 3D). The GABAergic marker Slc32a1 encoding the vesicular inhibitory amino acid transporter was expressed at the second highest level in cluster 4 relative to all other clusters (Figure 1—source data 1), consistent with the GABAergic phenotype of mature AGRP/NPY neurons (Figure 3—figure supplement 1).

Figure 3 with 2 supplements see all
Temporal gene expression patterns of the Pomc(low)/Otp cluster and Pomc(low)/Tac2 cluster.

(A and E) Violin plots showing the expression of signature genes at each developmental stage in Pomc(low)/Otp cluster and Pomc(low)/Tac2 cluster, respectively. (B) Fluorescence in situ hybridization showing the co-localization of Pomc (green) and Npy (red) at E15.5. (C and G) Heatmaps showing the gene ontology analysis based on the top marker genes in the order of pseudotime in Pomc(low)/Otp cluster and Pomc(low)/Tac2 cluster, respectively. (D) Uniform Manifold Approximation and Projection (UMAP) plots showing cells in the Pomc(low)/Otp cluster expressing Sst, Agrp, Npy, Sst/Agrp, Sst/ Npy, or Agrp/Npy. Cells expressing Sst and cells expressing both Agrp/Npy are from two distinct cell populations. (F) Fluorescence in situ hybridization showing the co-localization of Pomc (green) and Tac2 (red), the co-localization of Pomc (green) and Tacr3 (red), and the co-localization of Pomc (green) and Sox14 (red) at the indicated developmental stages. Immunofluorescence showing the co-localization of TDTOMATO (red) and ESR1 (green) at the indicated developmental stages. Insets are magnified views of the indicated boxes. Arrows indicate co-expressing neurons in the merged panels. Scale bar: 50 µm. Image orientation: left, posterior; right, anterior; Pit: pituitary gland, RP: Rathke’s pouch.

Cluster 5

Pomc(low)/Tac2 represents a distinct neuronal population from cluster 1 in its patterns of TFs associated with the development of mature POMC neurons and its low level of Pomc transcripts (Figure 3E). In particular, Isl1 was highly expressed only at E11.5 and Prdm12 at E11.5 and E13.5 in cluster 5. Starting from E15.5 through P12 there was a gradual onset of expression for a set of genes including Tac2, Ar, Esr1, Tacr3, Prlr, Pdyn, and Kiss1/Gm28040 isoforms (Figure 3E). Confirmations of co-expression for Pomc with either Tac2, Tacr3 and Sox14 or TDTOMATO with ESR1 at several embryonic ages are shown in Figure 3F. Interestingly, Sox14 has been reported to be necessary for expression of Kiss1 (Huisman et al., 2019). Most of the top enriched genes for this cluster are inferred to be expressed at the middle to end portions of the pseudotime prediction, and their functions are related to reproduction, synapse assembly/organization, pre/post-synaptic specialization, and neuron migration (Figure 3G). Taken together, this ensemble of genes is characteristic of adult Arc KNDY neurons that play an essential role in the control of reproductive physiology. A previous publication based on transgene lineage tracing also concluded that a significant subpopulation of KNDY neurons is derived from POMC progenitor cells (Sanz et al., 2015).

A prominent and unique characteristic of cluster 6. Pomc(low)/Dlx5 neurons was the expression of members of the Dlx family of homeobox TFs (Figure 3—figure supplement 2A and B). The six family members are usually expressed in three linked pairs of isoforms: 1 and 2, 3 and 4, or 5 and 6. Cluster 6 was enriched in Dlx1/2, and Dlx5/6 together with Dlxos1. The Dlx1/2 pair is known to activate expression of Ghrh in the mouse Arc (Lee et al., 2018). However, there were no transcripts for Ghrh expressed in cluster 6 or any of the other Pomc neuronal clusters suggesting that Dlx1/2 are not sufficient to induce Ghrh expression. Moreover, it has been recently shown that Dlx1/2 binding to Otp regulatory elements inhibits OTP production, subsequently reducing downstream Agrp expression (Lee et al., 2018). Cluster 6 had the highest expression of the opioid prepronociceptin (Pnoc) compared to the other integrated clusters, but low levels of preprodynorphin (Pdyn), similar to cluster 5 KNDY neurons.

Cluster 6

It had the highest Gad1, Gad2, and Slc32a1 (Vgat, vesicular inhibitory amino acid transporter) expression density of all clusters with virtually no Slc17a6 (Vglut2, vesicular glutamate transporter 2) expression, suggesting that these differentiating neurons were exclusively GABAergic. A comparison of GABAergic and glutamatergic markers among the eight neuronal clusters is shown in Figure 3—figure supplement 1. Cluster 4 neurons were also uniformly GABAergic. Cluster 1 neurons contained a mixture of glutamatergic and GABAergic markers with the exception of undetectable Slc32a1. However, they did express low levels of Slc18a2 (Vmat2, vesicular monoamine transporter 2), which possibly functions as an alternative vesicular GABA transporter in some dopamine neurons (German et al., 2015). This combination of mixed glutamatergic and GABAergic features in cluster 1 is characteristic of postnatal POMC neurons (Wittmann et al., 2013; Jones et al., 2019) and was also present to a limited extent in cluster 2. The remaining neuronal clusters 3, 5, 7, and 8 were all uniformly glutamatergic.

Cluster 7

Pomc(low)/Prdm13 contained the smallest number of cells among the eight neuronal clusters, and the cells were primarily limited to developmental ages E13.5–E17.5 (Figure 1B). Feature genes included Prdm13, Nr5a1, Adcypap1, Cnr1, Fam19a1, Rbp1, and Pdyn (Figure 3—figure supplement 2C). Prdm13 was reported recently (Chen et al., 2020) to be highly expressed in E15.5 POMC cells.

We further analyzed all of the neurons and non-neuronal cells by reclustering them based on their individual developmental ages (Figure 4) rather than for their gene expression profiles integrated across all ages. Average gene expression levels for each of these developmental subclusters are listed in Figure 4—source data 1 and the feature genes defining the developmental subclusters are listed in Figure 4—source data 2. Based on these data, it was then possible to trace the temporal continuity between developmental subclusters from age E11.5 to P12 relative to the original integrated cell cluster identities (Figure 4—figure supplement 1). These data are important because they connected the dynamic nature of gene expression profiles at each embryonic stage leading to the eventual transcriptomes present in the postnatal time points.

Figure 4 with 3 supplements see all
Uniform Manifold Approximation and Projection (UMAP) plots showing the distribution of major cell clusters at each developmental stage.

(A) E11.5, (B) E13.5, (C) E15.5, (D) E17.5, (E) P5, and (F) P12. The developmental subclusters derived from each primary cluster are colored to corresponding cells in the UMAP plot in Figure 1B.

Figure 4—source data 1

Average gene expression obtained from the Seurat analysis of each subcluster unique to the six developmental stages (corresponding to Figure 4 and Figure 4—figure supplement 1).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig4-data1-v2.xlsx
Figure 4—source data 2

Feature genes that define each subcluster at each of the six developmental stages (corresponding to Figure 4 and Figure 4—figure supplement 1).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig4-data2-v2.xlsx

Comparison of transcriptional profiles for neuropeptides, neuropeptide processing enzymes, neuroendocrine secretory proteins, and GPCRs in the eight neuronal POMC clusters at each developmental age

We identified 19 unique neuropeptide prohormone genes that were clearly differentially expressed features in at least one of the neuronal clusters at one or more developmental ages (Figure 4—figure supplement 2, Figure 2—source data 1). In addition to the neuropeptide genes already mentioned as characteristic features of certain clusters, additional neuropeptides were featured in other clusters primarily at postnatal ages P5 and P12. Notable examples are Gal, Tac1, Adcyap1, Cartpt, and Pthlh. Unlike the mosaic of neuropeptide gene expression profiles across clusters and developmental ages, the neuroendocrine secretory proteins, characteristic of dense core granules, were expressed ubiquitously in all eight neuronal clusters with similar gradients across all developmental ages (Figure 4—figure supplement 2). These included members of the chromogranin (Chga and Chgb) and secretogranin (Scg2, Scg3, and Scg5) families. Among the prohormone processing enzymes, Pcsk1 was prominently expressed only at P5 and P12 in most neuronal clusters while Pcsk2 was expressed at both embryonic and postnatal ages in the same clusters as Pcsk1. Pam was expressed in the same cluster/developmental age pattern as Pcsk1. Only Cpe was expressed in virtually all neurons at all developmental ages.

Split plots of G-protein-coupled receptor (GPCR) gene expression revealed a wide range of cluster- and age-specific patterns (Figure 4—figure supplement 3). Cnr1 (cannabinoid receptor 1) was an important feature primarily of cluster 7 and the two interrelated clusters 3 and 8. The two metabotropic GABA receptors Gabbr1 and Gabbr2 were strongly expressed in all eight neuronal clusters but differed in their developmental gradients. Gabbr1 was expressed at all developmental ages while Gabbr2 transcripts were largely present only in more mature neurons at postnatal ages P5 and P12. The opioid receptor Oprl1 that is selectively activated by nociception/orphanin FQ (Toll et al., 2016) stood out for its consistently strong expression in all eight neuronal clusters, particularly at postnatal days P5 and P12.

The prokineticin receptor 1 (Prokr1) was selectively expressed together with the GPCR accessory protein gene Mrap2 in cluster 1 Pomc(high)/Prdm12 neurons. Although MRAP2 was initially identified as an activating modulator of melanocortin receptor signaling, consequently reducing food intake, it was subsequently shown to be promiscuous in its interaction with additional GPCRs (Srisai et al., 2017). Unlike its action on MC4R signaling, MRAP2 inhibits PROKR1 signaling to promote food intake in mice independently of its interaction with MC4R (Chaly et al., 2016). Therefore, the co-expression of Prokr1 and Mrap2 in anorexigenic POMC neurons suggests an additional mechanism for PROKR1’s regulation of energy homeostasis by modulating the release of melanocortins from POMC neurons. Finally, the high expression of Npy2r (neuropeptide Y receptor 2) in both clusters 1 and 4 matches well with its expression in adult anorexigenic POMC neurons and orexigenic AGRP/NPY neurons, respectively.

TRAP-seq analysis of POMC-expressing cells at P12 and P60

To define the important genetic programs that direct the transition from early postnatal to adulthood and to compare the differences of the transcriptional programs that guide embryonic vs. postnatal/adult development, we performed Translating Ribosome Affinity Purification TRAP-seq using compound Pomc-CreERT2; Rosa26-eGFPL10a transgenic mice where POMC neurons are labeled with a eGFPL10a tag after tamoxifen administration at specific developmental stages. Colocalization of POMC and GFP immunoreactivity confirms the expression of eGFPL10a in Arc POMC neurons upon tamoxifen injection (Figure 5A). Anti-GFP-conjugated beads were used to pull down actively translating RNAs bound to the eGFPL10a ribosomes from hypothalamic extracts. Both the pull-down mRNA samples and the resultant supernatant mRNA samples were subjected to RNA sequencing. Principal component (PC) analysis shows distinct separation of pull-down and supernatant samples at both ages P12 and P60 from three independent experiments, suggesting good quality of these samples and the intrinsic differences between pull-down and supernatant samples (Figure 5B). Both Pomc and GFP were highly expressed in TRAP pull-down samples relative to the supernatants (Figure 5C), further validating the specificity of the TRAP-seq method. Gene enrichment analysis identified 1143 and 1047 highly enriched genes (p < 0.05) from Pomc- eGFPL10a P12 and P60 TRAP-seq pull-downs, respectively (Figure 5D, Figure 5—source data 1 and Figure 5—source data 2), from which 653 genes expressed in both datasets (Figure 5E). We hypothesize that the commonly expressed genes may be associated with the maintenance of POMC functional identity, whereas the uniquely expressed genes at P12 or P60 may be related to distinct age-dependent biological phenomena. To test this hypothesis, we performed gene ontology analysis on the 653 co-expressed genes, 394 P60 specifically expressed genes and 490 P12 specifically expressed genes. The results demonstrate that the commonly expressed genes are associated with the establishment and maintenance of neuron structure and basic neuronal functions. The top three functional annotations include the organization of synapses, regulation of neuron differentiation, and cytoskeleton-dependent intracellular transport. The ontology annotations on the uniquely expressed genes from P12 and P60 are dramatically different from each other. The functions of P12 uniquely enriched genes emphasize basic cellular metabolism, cellular modeling, and biochemical changes, whereas, the P60 enriched genes are related to the maturation of POMC neurons such as neurotransmitter development and POMC-specific function establishment including the responses to nutrient, and the construction of feeding behavioral circuits (Figure 5F). We next compared the TRAP-seq commonly expressed genes to our developmental scRNA-seq data to examine their expression levels in each cluster and each age. Approximately one-third (241) of the genes shown in the heatmaps on the left are abundantly expressing in the scRNA-seq dataset across all the clusters. Notably, compared to the four embryonic stages, a constantly higher transcript abundance was observed at both P5 and P12 ages (gray outlined boxes in Figure 5G), indicating great genetic similarities between P5 and P12. Intriguingly, the patterns of expression levels are very similar in both the P12 and P60 TRAP-seq datasets (Figure 5G, right heatmaps). To gain a better visualization of the distribution of TRAP-seq common genes expression at a single-cell level, we grouped the 241 genes as one module, calculated module scores based on the normalized counts and projected the scores to a UMAP plot of the eight neuronal clusters (Figure 5H). Most cells from cluster 1 (Pomc(high)/Prdm12), cluster 5 (Pomc(low)/Tac2), cluster 3 (Pomc(med)/Nr4a2), and a subcluster of cluster 2 (Pomc(med)/Ebf1) present a higher module score (high gene expression levels). Assessment together with Figure 6B confirms that most of the high scoring cells are derived from the two postnatal stages (P5 and P12).

Translating ribosome affinity purification with RNA-sequencing (TRAP-seq) showing gene enrichment and gene expression profiles at P12 and P60.

(A) Schematic diagram showing the generation of Pomc-CreERT2; ROSA26eGFP-L10a mice for TRAP-seq experiment. Immunofluorescence validated the co-localization of proopiomelanocortin (POMC) and eGFP. (B) Principal component analysis showing the separation of the RNA sequencing data from beads pull-down vs. supernatant at P12 (purple) and P60 (green), respectively. (C) Heatmap showing both Pomc and eGFP were highly expressed in beads pull-down samples. Rows are the biological replicates of each sample. Sup: Supernatant samples; PD: pull-down samples. Data were presented as scaled counts per million (CPM). (D) Gene enrichment plots showing 1143 genes and 1047 genes were significantly expressed in beads pull-down samples at P12 and P60, respectively (p < 0.05). (E) Venn diagram showing the number of genes highly enriched in both P12 (purple) and P60 (green) beads pull-down samples. (F) Gene ontology analysis showing the top 15 biological processes that were represented in P12 and P60 co-expressed genes, P12 uniquely expressed genes and P60 uniquely expressed genes. (G) Expression profile of the top enriched genes from both P12 and P60 TRAP-seq datasets across eight neuronal clusters; gray boxes indicate the higher expression of these genes in the two postnatal stages. The two heatmaps (right) indicating the top genes expression in P12 and P60 TRAP-seq datasets. (H) Uniform Manifold Approximation and Projection (UMAP) plot showing the distribution of the top enriched genes from both P12 and P60 TRAP-seq datasets in eight neuronal cell clusters. CPM: counts per million; org.: organization; reg.: regulation; pos.:positive; neg.: negative.

Figure 5—source data 1

Differentially expressed genes by RNA translating ribosome affinity purification with RNA-sequencing (TRAP-seq) between pull-down vs. supernatant at age P12 (corresponding to Figure 5).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig5-data1-v2.xlsx
Figure 5—source data 2

Differentially expressed genes by RNA translating ribosome affinity purification with RNA-sequencing (TRAP-seq) between pull-down vs. supernatant at age P60 (corresponding to Figure 5).

https://cdn.elifesciences.org/articles/72883/elife-72883-fig5-data2-v2.xlsx
Figure 6 with 1 supplement see all
RNA velocity and diffusion maps analyses illustrating the developmental trajectories of hypothalamic neuronal clusters.

(A) RNA velocity analysis showing the multiple origins of embryonic proopiomelanocortin (POMC) cells in eight neuronal clusters according to Seurat analysis cell embedding. (B) RNA velocity analysis showing the multiple origins of POMC cells correspond to the early embryonic stages. (C–D) Diffusion map showing the cell lineages of POMC neurons during embryogenesis. Cells are colored on the basis of their developmental stage (C) or based on their previously defined clusters (D). (E) Diffusion maps showing the expression of selected genes according to previous characterization of the feature genes representing each cluster. The color intensity corresponds to the normalized gene expression from Seurat analysis.

Developmental relationships of gene expression patterns in the eight neuronal POMC clusters at four embryonic ages

To investigate the gene expression trajectories of developing POMC neurons, we performed an RNA velocity analysis of cells from the eight neuronal clusters focusing on the four embryonic stages. This RNA velocity analysis revealed a major origin of POMC neurons captured by the Pomc(high)/Prdm12 cluster and a minor origin depicted by the Pomc(med)/Ebf1 cluster (Figure 6A and B), which is consistent with the early appearance of these two clusters at E11.5 (Figure 6—figure supplement 1). Further investigation of cells from E13.5, E15.5, and E17.5 also suggest a dual origin from Pomc(high)/Prdm12 or Pomc(med)/Ebf1 clusters (Figure 6—figure supplement 1).

To further validate the hypothesis of a dual origins of POMC neurons and to better visualize cell trajectories in three-dimensional Euclidean space, we constructed diffusion maps of all neurons from the four embryonic stages (Figure 6C). The image depicts a triangular beveled star polygon with rays 1 and 2 derived primarily from distinct groups of neurons at ages E11.5 and E13.5, a center vertex containing individual neurons from the intermediate ages E13.5 and E15.5 and a third ray with neurons from ages E15.5 and E17.5. The same data are alternatively color-coded for each neuron’s cluster identity (Figure 6D) showing that ray 1 is primarily derived from cluster 2, ray 2 from cluster 1, the center vertex from a convergence of all clusters 1–8 and ray 3 from clusters 1, 2, 4, 5, and 6. Compared to ray 2, cells in ray 1 are not continuously ordered, but instead show an abrupt transition at the junction of cells from E11.5 and E13.5.

We then overlaid the geometric location of each neuron with the expression levels of selected feature genes from each of the eight clusters (Figure 6E). Consistent with the diffusion maps analyses, the more granular data for individual genes and neurons illustrates the essentially identical patterns of high Pomc expression with the two originally discovered cognate transactivating factors Isl1 and Nkx2-1 across embryonic development (rays 2 and 3). However, Isl1 was also expressed in numerous Pomc(low) neurons at E11.5 and E13.5 (ray 1). Prdm12 and Tbx3 exhibited complementary temporal patterns of high expression at early (ray 2) vs. late (ray 3) embryonic ages, respectively. The combined developmental and geometric patterns of Six3, Six6, and Sox14 expression were all very similar to both Pomc and Nkx2-1, suggesting that all four of these TF genes might be necessary to define the cell fate of a subpopulation of hypothalamic neuronal progenitors into mature differentiated POMC neurons.

In contrast to those genes with expression patterns that defined the continuum of neuronal identity from ray 2 to ray 3, another set of genes including Ebf1, Pou4f1, and Nr4a2 was similarly highly expressed only at the two earliest embryonic ages in neurons constituting ray 1. The expression of Pitx2, which characterizes cluster 8 and a portion of cluster 3, was high in ray 1 as it merges into the center vertex and was mostly derived from age E13.5. Finally, the temporal and geometric patterns of high expression levels for Otp, Prdm13, Tac2, and Dlx5 were distinct from each other and the aforementioned genes, consistent with their ontogeny into the distinct differentiated cell fates of cluster 4. AGRP/NPY/GABA, cluster 5. KNDY, and cluster 6. GABA/PNOC neurons, respectively.

Discussion

In this study, we tracked the genetic programs that give rise to Pomc-expressing neurons in the developing hypothalamus and followed their progression into terminally differentiated POMC neurons and alternative phenotypic destinies. To this end, we performed an scRNA-seq study specifically of Pomc-expressing neurons present in developing mouse hypothalami at four embryonic time points (E11.5, E13.5, E15.5, E17.5) and two early postnatal days (P5 and P12). We also performed a complementary TRAP-seq study on affinity purified translating RNAs derived from hypothalamic POMC cells obtained at ages P12 and P60. These developmental stages were selected based on the ontogeny of Pomc expression in the presumptive hypothalamus (Japón et al., 1994), peak of neurogenesis of POMC neurons (Nasif et al., 2015; Orquera et al., 2019), subsequent developmental maturation (Padilla et al., 2010; Sanz et al., 2015), and terminal differentiation (Quarta et al., 2019).

Our study revealed the existence of a major group of POMC neurons (cluster 1. Pomc(high)/Prdm12) that showed by far the highest levels of Pomc UMIs relative to all other neuronal clusters at all ages. We believe that this cluster constitutes the canonical POMC neurons that further mature and integrate into functional circuits in the postnatal hypothalamus. Cluster 1 is featured by a unique combinatorial set of TFs that includes Isl1, Nkx2-1, and Prdm12. In the last few years, our laboratories demonstrated that these three TFs are present in neurons that start expressing Pomc in the developing arcuate nucleus at E10.5, and their expression continues at later developmental ages and adulthood in hypothalamic POMC neurons. Indeed, we showed that ISL1 (Nasif et al., 2015), NKX2-1 (Orquera et al., 2019), and PRDM12 (Hael et al., 2020) are absolutely necessary to specify the identity of POMC neurons. Moreover, we also demonstrated previously that adult mice lacking any of these three TFs exclusively from POMC cells display greatly reduced levels of hypothalamic Pomc mRNA together with increased adiposity and body weight. However, our results also indicated that the sole presence of ISL1, NKX2-1, and PRDM12 at E10.5–E11.5 is insufficient to determine the early identity of hypothalamic POMC neurons, suggesting that the combinatorial set of TFs necessary for neuronal-specific Pomc expression has additional early components, yet to be discovered. In this regard, our scRNA-seq study revealed several novel candidates in cluster 1 cells including Six3, Six6, Sox1/2/3, Sox14, and other TFs (Figure 2). Future molecular and functional genetic studies will be needed to assess their potential contribution to Arc Pomc expression and the functioning of POMC neurons.

Six3 and Six6 are two evolutionarily related TF genes known to play critical roles in the development of the eyes, forebrain, and pituitary (Liu and Cvekl, 2017; Liu et al., 2010; Suh et al., 2010). The expression patterns of Six3 and Six6 highly overlap during early embryogenesis but segregate into different territories during late embryonic development, indicating possible distinct postnatal functions (Geng et al., 2008). Indeed, Six3 knockout mice die at birth whereas mice lacking Six6 exhibit decreased numbers of GNRH neurons and impaired fertility (Larder et al., 2011). Although Six3 and Six6 are abundantly expressed in Arc POMC neurons across all examined developmental stages (Figure 2), their participation in the control of Pomc expression and the patterning and function of POMC neurons requires further investigation. Regarding Sox genes, a family of TFs known to regulate cell fate by transactivating cell-specific genes, in silico analysis of conserved SOX protein binding motifs suggests the potential binding of SOX proteins including SOX1, SOX2, SOX3, and SOX14 to Pomc neuronal enhancers nPE1 and nPE2 (Chen et al., 2020). All four of these Sox genes were highly enriched in Pomc(high)/Prdm12 cluster 1 (Figure 2 and Figure 1—source data 2) as well as in Pomc(low)/Tac2 cluster 5 and Pomc(low)/Prdm13 cluster 7. Furthermore, mice lacking Sox14 exhibit a great reduction in KNDY neurons (Huisman et al., 2019).

A particular feature gene present in cluster 1 was Tbx3, a T-box TF that has been recently reported to play a key role in Arc Pomc expression (Quarta et al., 2019). Unlike Isl1, Nkx2.1, and Prdm12, which are all strictly required at the onset of Pomc expression (E10.5), Tbx3 expression commences several days later (E15.5) in the Arc and, therefore, apparently does not participate in the early specification of POMC neurons. However, Tbx3 plays an important role in the functional maturation of these neurons as a transcriptional booster of arcuate Pomc expression. It was recently demonstrated in mice specifically lacking Tbx3 from POMC neurons that Pomc expression was diminished leading to hyperphagia and obesity (Quarta et al., 2019).

The TRAP-seq results add another layer to our understanding of the genetic machinery that leads to mature functional POMC neurons. Genes that are uniquely enriched at P12 are associated with the basic cellular biology process and cell metabolism whereas genes specifically expressed at P60 are involved in the establishment of POMC neuronal functions. The co-expressed genes from TRAP-seq are highly expressed at P5 and P12 in the single-cell sequencing datasets, suggesting that the genetic programs that dominate embryonic POMC neuronal development are substantially different from the programs that guide postnatal POMC neuronal differentiation. Although only 1 week apart, the cells at P5 share more genetic similarities to the cells at P12 rather than the cells at E17.5.

An unanticipated finding of our study is the discovery of an early population of Arc Pomc-expressing neurons characterized by the expression of the TF Ebf1 and Pou4f1, in the absence of Prdm12 transcripts. The expression profile of cluster 2. Pomc(med)/Ebf1 differs from cluster 1. Pomc(high)/Prdm12 most particularly at the earliest developmental ages E11.5–E15.5 and may both represent alternative genetic routes leading to mature arcuate neurons. The early expression of Ebf1 suggests that this TF may regulate the migration of differentiating neurons out of the ventricular zone to the mantle, as has been reported in the developing striatum (Garel et al., 1999). Further studies are needed to understand the functional significance of the non-canonical population of Pomc-expressing neurons captured by cluster 2 cells.

Another main finding of our study is the identification of multiple subpopulations of Pomc-expressing neurons that originate in the developing arcuate nucleus. In addition to the previously characterized fertility-regulating KNDY neurons (Pomc(low)/Tac2) (Sanz et al., 2015) and orexigenic AGRP/NPY neurons (Pomc(low)/Otp) (Padilla et al., 2010), our data indicated that, as early as E11.5, Pomc-expressing neurons may also give rise to four other less characterized neuronal populations noted as Pomc(med)/Nr4a2, Pomc(low)/Pitx2, Pomc(low)/Dlx5, and Pomc(low)/Prdm13. Moreover, subclustering cells from E11.5 revealed a distinct cell population featured by Otp expression (E11.5–9 cluster in Figure 4), suggesting that the presence of Otp transcripts at E11.5 may already be dictating the differentiation pathway of this neuronal cluster toward an AGRP phenotype, a hypothesis that requires further investigation. In addition, RNA-velocity developmental trajectory analysis (Figure 6A and Figure 6—figure supplement 1B) suggested that the early cells in Pomc(high)/Prdm12 cluster project to Pomc(low)/Otp cluster.

Pomc(med)/Nr4a2 and Pomc(low)/Pitx2 are related clusters sharing 20% of their top 150 feature genes being one of the most salient differences the presence of neuropeptide Cck transcripts only in the Pomc(low)/Pitx2 cluster. Comparing our dataset with a previous one obtained using FACS from adult hypothalami of POMC-EGFP mice (Lam et al., 2017) in which Nr4a2 showed to be expressed in 18.5% of cells in cluster 4, it is likely that embryonic Pomc(med)/Nr4a2 cells captured by our study keep their identity through adulthood. Similarly, signature genes for the Pomc(low)/Dlx5 cluster including Dlx1, Dlx2, Dlx5, and Dlx6 were also found to be highly enriched and specific for cluster 1 in Lam et al.’s study, suggesting that the Pomc(low)/Dlx5 cluster detected in our study becomes a subset of adult POMC neurons. Finally, our identification of the Pomc(low)/Prdm13 cluster matches with a recent study showing that Prdm13 is highly enriched in FACS cells from E15.5 POMC-EGFP embryos analyzed by whole transcriptome RNA sequencing (Chen et al., 2020).

It is worth mentioning that despite some extent of similarities between embryonic and adult POMC cells, the majority of embryonic and early postnatal cells profiled in our study are distantly related to mature adult POMC neurons, indicating that terminal differentiation of POMC neurons is a long-lasting maturation process that proceeds along multiple life stages. For example, most adult POMC cells express proprotein convertase PC1/3 (Pcsk1) and leptin receptor (Lepr), whereas only 22% and 2.7% of cells from our datasets express the two genes, respectively, suggesting critical functional changes in the transcriptome of POMC neurons during postnatal maturation after early postnatal ages.

Due to the limitations of the scRNA-seq techniques, some interesting questions cannot be answered from our study at the functional level. For example, are all TFs present in cluster 1. Pomc(high)/Prdm12 at E11.5 critical for dictating the identity of the cannonical POMC neurons? what TFs are critical for maintaining POMC neuron identity throughout lifetime? Is Otp the only critical gene that determines the transition of Pomc-expressing cells into AGRP/NPY neurons? Which are the TFs that give rise to the lineage of KNDY neurons during development? We believe that our study provides a valuable resource to tackle these and many other questions aimed at a better understanding of the regulation and functioning of POMC and other related neurons in the mammalian arcuate nucleus of the hypothalamus.

Concluding remarks

In summary, this study is the first to comprehensively characterize the transcriptomes of Arc hypothalamic POMC cells during embryonic and early postnatal development. This dataset extends our understanding of the diversification of early POMC cells and provides a valuable resource for further elucidating the regulatory mechanism for Pomc expression and neuronal maturation. The identification of new marker genes in each subpopulation from scRNA-seq and uniquely expressed genes at P12 or P60 from the TRAP-seq study will potentially lead to new directions for future functional studies of POMC neurons.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (Mus musculus)Pomc -TdDsRedPMID:19864580Dr Malcolm Low (University of Michigan)
Genetic reagent (Mus musculus)Pomc-CreERT2PMID:24177424RRID:MGI:5569339
Genetic reagent (Mus musculus)Rosa26eGFP-L10aJackson LaboratoryStock # 024750;RRID: IMSR_JAX:024750
AntibodyAnti-Nr5a1(Rabbit Polyclonal)Dr Gary Hammer, University of MichiganIF (1:750–1:1000)
AntibodyAnti-td-Tomato [16D7] (Rat Monoclonal)KerafastCat # EST203;RRID:AB_2732803IF (1:750)
AntibodyAnti-Pomc (Rabbit Polyclonal)Phoenix PharmaceuticalsCat # H-029–30;RRID:AB_2307442IF(1:1000)
AntibodyAnti-Rabbit Alexa 488 (Goat polyclonal)InvitrogenCat # A-11034;RRID:AB_2576217IF(1:500–1:1000)
AntibodyAnti-Rat Alexa Fluor 555 (Goat polyclonal)InvitrogenCat # A-21434;RRID:AB_2535855IF(1:500–1:1000)
AntibodyAnti-Esr1 (Rabbit Polyclonal)Dr Sue Moenter, University of MichiganIF(1:10000)
AntibodyAnti-GFP (Chicken polyclonal)AbcamCat # ab13970;RRID:AB_300798IF(1:1000)
Chemical compound, drugPapain Dissociation SystemWorthington Biochemical CorporationLK003178
Chemical compound, drugDNaseWorthington Biochemical CorporationLK003172
Software, algorithmRR Project for Statistical Computinghttps://www.r-project.org/
Software, algorithmSeurat 3.1.5Satija Labhttps://satijalab.org/seurat/articles/install.html
Software, algorithmMonocle 2Cole-trapnell labhttp://cole-trapnell-lab.github.io/monocle-release/docs/
Software, algorithmDestiny 2.14.0Carsten Marr & Florian Buettner labhttps://theislab.github.io/destiny/
Sequence-based reagentRNAscope Multiplex Fluorescent V2 AssayACDbioCat # 323110
Sequence-based reagentMm-Six3Advanced Cell DiagnosticsCat # 412941-C3
Sequence-based reagentMm-Six6Advanced Cell DiagnosticsCat # 574291
Sequence-based reagentMm-Tac2Advanced Cell DiagnosticsCat # 446391-C3
Sequence-based reagentMm-NpyAdvanced Cell DiagnosticsCat # 313321
Sequence-based reagentMm-Ebf1Advanced Cell DiagnosticsCat # 433411
Sequence-based reagentMm-Pou4f1Advanced Cell DiagnosticsCat # 414671-C3
Sequence-based reagentMm-PomcAdvanced Cell DiagnosticsCat # 314081-C2
Sequence-based reagentMm-Sox14Advanced Cell DiagnosticsCat # 516411
Sequence-based reagentMm-Prdm13Advanced Cell DiagnosticsCat # 543551-C2
Sequence-based reagentMm-Tacr3Advanced Cell DiagnosticsCat # 481671
Sequence-based reagentMm-Dlx5Advanced Cell DiagnosticsCat # 478151
Sequence-based reagentMm-Pitx2Advanced Cell DiagnosticsCat # 412841
Sequence-based reagentMm-PomcAdvanced Cell DiagnosticsCat # 314081
Sequence-based reagentMm-Positive Control ProbeAdvanced Cell DiagnosticsCat # 320881
Sequence-based reagentMm-Negative Control ProbeAdvanced Cell DiagnosticsCat # 320871
Sequence-based reagentRNAscope Probe DiluentAdvanced Cell DiagnosticsCat # 300041
Chemical compound, drugOpal 520 Reagent PackAKOYA BiosciencesCat # NC16018771:750–1:1500
Chemical compound, drugOpal 620 Reagent PackAKOYA BiosciencesCat # NC16120591:750–1:1500
Chemical compound, drugOpal 690 Reagent PackAKOYA BiosciencesCat # NC16050641:750–1:1500
OtherDAPIBD BiosciencesCat # 564907FACS (500–1000 ng/mL)
OtherProLong Gold Antifade MountantThermoFisherP36930
OtherEarle’s Balanced SaltsSigma-AldrichE2888

Mice

All procedures were performed in accordance with the Institutional Animal Care and Use Committee (IACUC) at the University of Michigan and followed the Public Health Service guidelines for the humane care and use of experimental animals. Mice were housed in ventilated cages under controlled temperature and photoperiod (12 hr light/12 hr dark cycle, lights on from 6:00 AM to 6:00 PM), with free access to tap water and laboratory chow (5L0D, LabDiet). Breeding mice were fed with the breeder chow diet (5008, LabDiet). Transgenic mice expressing the fluorescent protein Tdimer-Discosoma red (TdDsRed) in POMC neurons were generated previously (Hael et al., 2020). The vast majority of Pomc-TdDsRed cells were validated as authentic POMC neurons based on co-localization of ACTH immunostaining with the TdDsRed fluorophore (Hentges et al., 2009) Pomc-CreERT2 and Rosa26eGFP-L10a (Stock no. 024750, The Jackson Laboratory) mice were generated as previously described (Liu et al., 2014). To induce POMC neurons-specific expression of eGFP-L10a at P12 or P60, tamoxifen (50 mg/kg) was injected intraperitoneally for 5 constitutive days from P6 to P10 or P50 to P54.

Generation of single-cell suspension for scRNA-seq

Request a detailed protocol

Adult male Pomc-TdDsRed/+ mice were bred with female Pomc-TdDsRed/+ mice and the day of copulation plug detection was counted as embryonic day 0.5 (E0.5). Embryonic hypothalami were dissected at E11.5, E13.5, E15.5, and E17.5. Postnatal hypothalami were isolated from pups at postnatal day 5 (P5) and P12. In each dissection, hypothalami from at least six pups were pooled together to acquire enough cells for FACS. Pooled tissue samples were digested in papain solution supplemented with DNAse for 20 min (E11.5 and E13.5), 30 min (E15.5), 40 min (E17.5), or 1 hr (P5 and P12) with gentle agitation based on a published protocol (Garel et al., 1999). The isolated cells were stained with DAPI as an indicator for cell viability and sorted by FACS for the transgenic red fluorophore. The collected cells were counted and subjected to scRNA-seq library preparation and sequencing (10× Genomics).

scRNA-seq data processing

Request a detailed protocol

A total of 37,053 single cells (E11.5: 5329, E13.5: 7051, E15.5: 5148, E17.5: 2488, P5: 7551, and P12: 9486) from six different developmental stages were processed using the 10× Genomics Chromium system. The libraries were sequenced on Illumina HiSeq 4000 and NovaSeq platforms. We obtained a total of over 2.6 billion reads with an average of 70,7896 reads per cell. Over 85% of reads mapped confidently to the mouse genome across all six developmental stages. Raw reads were processed with Cell Ranger (v.2.2 and v.3.0). Seurat package (v.3.1.5) (Butler et al., 2018) was used for downstream analysis. Since the number of genes and UMI counts varied across developmental stage, we applied different criteria at each stage to filter out possible doublets and low-quality cells. Specifically, at E11.5, we removed outlier cells that had UMI counts <5000 or >60,000 and gene counts <1000 or >8000 (determined by the visualization of UMI counts and gene count distributions). At E13.5, we removed outlier cells that had UMI counts <2500 or >30,000 and gene counts >1500. At E15.5, we removed outlier cells that had UMI counts >40,000 and gene counts <1000 or >6000. At E17.5, we removed outlier cells that had UMI counts >25,000 and gene counts <1000 or > 5000. At P5, we removed outlier cells that had UMI counts <1500 or >60,000 and gene counts <1500 or >8000. At P12, we removed outlier cells that had UMI counts <2000 or >60,000 and gene counts <1500 or >8000. Moreover, cells with high proportions of mitochondrial genes (>10%) or hemoglobin genes (>10%) were filtered out. Finally, we removed all the cells without any Pomc transcript UMI counts. A total of 13,953 cells passed the criteria (E11.5: 1498, E13.5: 1796, E15.5: 2078, E17.5: 1139, P5: 3909 and P12: 3533) for downstream analysis.

The gene expression level for each cell was normalized by the total expression, multiplied by a scaling factor of 10,000 followed by log-transformation (Butler et al., 2018). The top 2000 most variable genes were selected from each age and 60 integration anchors were used to combine all the datasets together. After scaling and centering the integrated dataset, we performed PCs analysis on the data matrix and the top 50 PCs were selected based on JackStraw and Elbow plots from Seurat package (Hentges et al., 2009) for data visualization using Uniform Manifold Approximation and Projection (UMAP) technique. We next constructed a shared nearest neighbor (SNN) graph by setting an expected number of neighbors to 50. In order to cluster the cells, an SNN modularity optimization technique within a function FindClusters was used to group cells together with a resolution parameter 0.1. The identity for each cluster was assigned based on the level of Pomc transcripts and prior knowledge of marker genes. The similar analysis was applied to define TdDsRed cell clusters (Figure 1—figure supplement 2 and Figure 1—figure supplement 3). In addition to the above criteria applied to remove unwanted cells, we chose cells with at least one TdDsRed transcript UMI count and used the top 2000 highly variable genes for PCs analysis. The top 30 PCs were selected for unsupervised clustering with a resolution of 0.3.

Unsupervised cell clustering from each developmental age (Figure 4) was conducting using the following PC numbers (E11.5: 25; E13.5: 24, E15.5: 26, E17.5: 19, P5: 40, P12: 34) with their corresponding resolutions (E11.5: 0.2; E13.5: 0.6, E15.5: 0.5, E17.5: 0.5, P5: 0.5, P12: 0.5). The UMAP plots, violin plots, feature plots, and dot plots were all generated in the Seurat package. The average expression levels for each gene within each developmental subcluster (Figure 4—source data 1) were obtained from Seurat and calculated by the average of [exp(logNormalized counts)–1]. The differentially expressed genes in each identified cluster was identified by the comparison of gene expression levels in a specific cluster to all the other clusters.

Cluster identity comparisons

Request a detailed protocol

To compare cluster identities between POMC cells and DSRED cells (Figure 1—figure supplement 2), we normalized the DSRED cells data and acquired the top 2000 most variable genes. We then use the standard datasets integration method from Seurat packages with the top 50 anchors for data integration and top 50 anchors to transfer cell identities from POMC cells to DSRED cells. We also applied this approach to conduct age-pairwise comparisons (Figure 4—figure supplement 1) with default setting from Seurat V3.1.5 package.

Cell lineage construction

Request a detailed protocol

For RNA velocity analysis, spliced and unspliced matrices of reads were summarized using velocyto (v.0.17.16) with default parameters (La Manno et al., 2018). Low complexity and repeat regions were downloaded from the UCSC browser (rmsk table from mm10). scVelo was performed for RNA velocity analysis (v.0.2.2) (Bergen et al., 2020). The cell embedding information was acquired from Seurat analysis. The top 4000 genes, the top 20 PCs, and the top 20 neighbors under stochastic mode were used to generate velocity graph.

Cell lineage was constructed using the Destiny (v.2.14.0) package, which implements the formulation of diffusion maps (Angerer et al., 2016). Diffusion maps are a spectral method for non-linear dimension reduction, which is especially suitable for analyzing single-cell gene expression data from different time-courses. We removed VLMC, glial cells, and astrocytes (clusters 9, 10, and 11 in Figure 1) and kept only neurons at embryonic stages (E11.5, E13.5, E15.5, and E17.5), resulting in a total of 6005 cells for this analysis. The log-transformed normalized data from Seurat data slot and the corresponding annotation information were imported to construct an expression matrix. The diffusion maps were constructed under the default setting for the Gaussian kernel width sigma (σ) and 300 nearest neighbors. The first two diffusion components (DCs) were used to visualize the results. The 3D plots were produced using rgl package (v.0.100.50) with the top three DCs.

To perform pseudotime gene expression analysis, cells from cluster 1, cluster 2, cluster 4, and cluster 5 were extracted from the Seurat object, respectively. Raw counts were acquired to construct Monocle (v.2.4.0) (Qiu et al., 2017a, Trapnell et al., 2014; Qiu et al., 2017b) datasets. Cells were ordered along the pseudotime by setting E11.5 as the root state. Differentially expressed genes (Log2FoldChange > 0.25 and p < 0.05) in each cluster from the Seurat analysis were plotted along the pseudotime using ‘plot_pseudotime_heatmap’ function. Gene ontology analysis (Biological Process) was performed using ClusterProfiler (v.3.18.1) (Yu et al., 2012).

TRAP-seq analysis

Request a detailed protocol

Mice (Pomc-CreERT2; ROSA26eGFP-L10a) were euthanized and decapitated. The brain was removed from the skull and the arcuate nucleus was isolated from 2 mm thick coronal slices using a brain matrix. Tissues were then homogenized and the subsequent tissue lysate was subjected to immunopurification steps according to the previous protocol (Heiman et al., 2014). The ribo-depletion kit (RiboGone, Takara, CA) and SMARTer Stranded total RNA sample prep kit were used to remove excess ribosomal RNA and synthesize cDNA library. Samples were sequenced on a 50-cycle single end run on a HiSeq 4000 (Illumina) according to manufacturer’s protocols.

Raw sequencing data was processed at the University of Michigan Bioinformatics core. Briefly, the quality of the raw reads data was checked using FastQC (v.0.11.3) and the filtered reads were aligned to reference genome (UCSC mm10) using TopHat (v.2.0.13) and Bowtie2 (v.2.2.1) with default parameters. The HTSeq/DEseq2 method was used for differential expression analysis with paired-samples and treatment (pull-down vs. supernatant) as the main effects.

Tissue collection and immunofluorescence staining

Request a detailed protocol

Depending on developmental stage, embryos at E10.5–E17.5 were fixed in 4% paraformaldehyde in phosphate buffer solution (PBS) at 4°C for a various period of time. Specifically, embryos at E10.5–E12.5 were fixed for 1 hr; heads from embryos at E15.5–E16.6 and E17.5 were fixed for 2 and 4 hr, respectively. Tissues were stabilized in 10% sucrose/10% gelatin in PBS at 37°C for 30 min prior to embedding in OCT compound as described previously (Butler et al., 2018). Embryonic brains were sectioned sagittally on a cryostat (Leica CM1950) at 20 μm thickness. For dual immunofluorescence, after washing excess OCT compound with PBS, a heat-induced antigen retrieval process was performed using citrate buffer (10 mM anhydrous citric acid and 0.05% Tween-20, pH 6.0) at 80°C for 30 min followed by two washes in PBS. Sections were then blocked with 3% normal goat serum and 0.1% Triton X-100 for 1 hr and incubated with primary antibodies at room temperature overnight. After PBS washes, sections were incubated with goat secondary antibodies for 2 hr at room temperature. Nuclei were stained with DAPI (1 mg/L) for 10 min and the slides washed five times with PBS before mounting with ProLong Gold Antifade Mountant (ThermoFisher). Sections were imaged on a Nikon 90i fluorescence microscope with NIS-Elements software. Information on the sources and dilution of antibodies are listed in Key resources table.

In situ hybridization

Request a detailed protocol

FISH was performed using RNAscope Multiplex Fluorescent V2 Assay (Advanced Cell Diagnostics) according to the manufacturer’s instructions with slight modifications. Specifically, after washing with PBS, sections were post-fixed with 4% paraformaldehyde for 20 min, and dehydrated in a series of ethanol solutions (50%, 70%, and 100%). After drying, dehydrated sections were incubated in RNAscope Hydrogen Peroxide solution for 10 min followed by the protease treatment using RNAscope Protease IV for 20 min. Sections were then hybridized to target RNA probes for 2 hr in the HybEZ II hybridization oven. Hybridized probes were amplified using a cascade of signal amplification solutions (AMP 1–3) followed by standard signal developing protocol as described in the manufacturer’s brochure. Detailed information on RNA probes and dilutions are listed in Key resources table. Confocal images were obtained using a Nikon Instruments A1 Confocal Laser Microscope with NIS-Elements software.

Data availability

All raw data have been deposited in the Gene Expression Omnibus under accession numbers GSE154153 and GSE181539.

The following data sets were generated
    1. Hui Y
    2. Marcelo R
    3. Malcolm L
    (2021) NCBI Gene Expression Omnibus
    ID GSE154153. Single cell RNA-seq of the developing mouse hypothalamus.
    1. Hui Y
    2. Marcelo R
    3. Malcolm L
    (2021) NCBI Gene Expression Omnibus
    ID GSE181539. Developmental single-cell transcriptomics of hypothalamic POMC progenitors reveal the genetic trajectories of multiple neuropeptidergic phenotypes.

References

Decision letter

  1. Joel K Elmquist
    Reviewing Editor; University of Texas Southwestern Medical Center, United States
  2. Marianne E Bronner
    Senior Editor; California Institute of Technology, United States

Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "Developmental single-cell transcriptomics of hypothalamic POMC progenitors reveal the genetic trajectories of multiple neuropeptidergic phenotypes" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Marianne Bronner as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) Showing in situ hybridization data for E11.5 markers in a coronal plane, where the ventricular and mantle layer can de distinguished, would help to resolve the question of whether Ebf-hi expression represents newly differentiated neurons that are in the process of migrating out of the ventricular zone. (This appears to be case from the E13.5 section in Figure 1F, but it is hard to tell.)

2) Addition of scRNA-seq data from earlier embryonic stages (i.e. E10.5 and/or E11.0) would help to identify the precursor to the population of immature POMC neurons that is captured by cluster 1.POMC-hi/Prdm12.

3) One question raised by these data are the presence of an inhibitory population of GABA releasing (ie Vgat-expressing) POMC neurons. This is controversial with some studies reporting this and others indicating that there are no GABAergic POMC neurons. Given the evidence for a developmental cluster that is evident even at P12, it is important to resolve this issue in the context of the authors report of a bona fide transcriptional cluster with these co-expression patterns. The authors should examine the colocalization of these Pomc/Vgat neurons using RNA-FISH and report their distribution across development.

4) Finally, several conclusions need to be tempered. Nonetheless, the paper will be an outstanding resource for the field. The discussion should be written to reflect the important resource of this data set.

Reviewer #1 (Recommendations for the authors):

There seemed to be a substantial subset of microglia expressing Pomc (Figure 1), did the authors remove these from the dataset? How do they determine a bona fied Pomc neuron vs a neuron with Pomc mRNA contamination.

A supplemental model figure/cartoon integrating the various developmental trajectories with their defining transcriptional factors would be useful to visualize the key findings of the paper.

Reviewer #2 (Recommendations for the authors):

• Showing in situ hybridization data for E11.5 markers in a coronal plane, where the ventricular and mantle layer can de distinguished, would help to resolve the question of whether Ebf-hi expression represents newly differentiated neurons that are in the process of migrating out of the ventricular zone. (This appears to be case from the E13.5 section in Figure 1F, but it is hard to tell.)

• Addition of scRNA-seq data from earlier embryonic stages (i.e. E10.5 and/or E11.0) would help to identify the precursor to the population of immature POMC neurons that is captured by cluster 1.POMC-hi/Prdm12.

• The number of total neurons changes dramatically across gestation, in large part because cells that adopt non-POMC fates turn off Pomc expression and are no longer captured by FACS. In Figure 1, reporting the absolute number of neurons in each cluster, instead expressing it as a percentage, would make it easier to interpret.

• Analyses shown in Figure 3 Supplemental Figures 3-4 are critical to resolve issues related to POMC progenitor heterogeneity. This information should be moved to the main figures. Maintaining a consistent color scheme between Figure 1 and Figure 3 Supplement 3 would make it much easier to interpret this data. The chart in Supplemental Figure 4 could then remain as a supplemental figure to reinforce this point.

• Trajectory analyses shown in Figure 5 should include separate analyses of 2.Pomc-med/Ebf1 and 3.Pomc-med/Nra2, as they appear to map onto different populations in the mature ARH (Supplemental Figure 15)

• There are discrepancies between the text (lines 489-497) and Supplementary file 14, which do not link 3.Pomc-med/Nr4a2 neurons to mature POMC populations and Supplementary file 15, which does, making it difficult to interpret these data.

Reviewer #3 (Recommendations for the authors):

This manuscript by Yu et al. beautifully captures the transcriptional heterogeneity of mouse POMC neurons across hypothalamic development. This study unifies multiple other observations about the role for other neuron al cell types that express POMC transiently during development. The paper is an important contribution to understanding of the diversity of POMC neuron classes and their relationship to other cell types in the arcuate nucleus.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "Developmental single-cell transcriptomics of hypothalamic POMC neurons reveal the genetic trajectories of multiple neuropeptidergic phenotypes" for further consideration by eLife. Your revised article has been reviewed by 2 peer reviewers and the evaluation has been overseen by Marianne Bronner as the Senior Editor, and a Reviewing Editor.

The revisions by the authors improved the manuscript. These data will be a valuable resource for the field. However, there are a few important points that should be corrected or clarified in a final version.

1. Characterization of Ebf1 expression

The claim that Ebf1 expression at E11.5 is localized to the mantle layer (line 201) and the labels in the accompanying Figure 1E are not accurate. Expression of Pomc and Ebf1 are located within the ventricular zone (VZ), albeit toward the lateral edge, as indicated by the dense DAPI stain. The mantle layer "ML" is label is not accurate; there is almost no "mantle layer" in the mediobasal hypothalamus at E11.5. Also, the differentiating zone indicated at the lateral edge of the VZ in Figure 1E is not striatal.

Note that Figure 2-Figure Suppl 1A shows co-expression of Pomc and Ebf1 at the border between the VZ and ML at E13.5. There appears to be more co-expression on the VZ side than the ML side. The expression pattern shown here is consistent with reports in other regions of the developing CNS, where EBFs have been implicated in regulating migration out of the VZ.

2. Lines 474-5 "The diffusion maps analysis appears to confirm the multiple origins of POMC progenitor cells demonstrated by the RNA velocity study."

This is a confusing sentence because "POMC progenitors" is not the accurate terminology to use. It is easier to just eliminate this sentence.

3. Discussion of Ebf1+/cluster 2 neurons

Lines 614-6 "The expression profile of cluster 2.Pomc(med)614 /Ebf1 differs from cluster 1.Pomc(high)615 /Prdm12 most particularly at the earliest developmental ages E11.5 to E15.5 and may both represent alternative genetic routes leading to mature POMC neurons."

There is no evidence that the Ebf1+/cluster 2 neurons develop into mature neurons of any sort. The analysis in Figure 6 depicts the Ebf1+/cluster 2 neurons as a dead end in Ray 1 (i.e. retains immature properties). Since the transcriptomics performed here relied on the expression of a Pomc-GFP transgene, it is possible that cluster 2 neurons differentiate into mature non-POMC cell types. The sentence should be corrected by swapping "arcuate" for "POMC".

4. Discussion of developmental origins of NPY and KNDy neurons Lines 619-626 are confusing.

Timing question: With regard to NPY and KNDy neurons, which are known to derive from Pomc-expressing neurons, the authors argue that they arise "earlier than expected "(presumably at E11.5). And yet, according to the developmental clustering (Figure 4) and developmental trajectory analyses (Figure 6 supplement 1), meaningful numbers of Otp+/cluster 4 (future NPY) and Tac2+/cluster 5 (future KNDy) neuron precursors emerge at E13.5, consistent with reports in the literature.

Developmental lineage: These analyses also support the idea that Otp+/cluster 4 (future NPY) neurons share the same developmental trajectory as Prdm12+/cluster 1 (future POMC) neurons (Figure 6) as well as a common precursor (E13.5-4, Figure 4). This should be explicitly discussed. These studies do not provide clues into the origin of the Tac2+/cluster 5 (future KNDy) neurons.

https://doi.org/10.7554/eLife.72883.sa1

Author response

Essential revisions:

1) Showing in situ hybridization data for E11.5 markers in a coronal plane, where the ventricular and mantle layer can de distinguished, would help to resolve the question of whether Ebf-hi expression represents newly differentiated neurons that are in the process of migrating out of the ventricular zone. (This appears to be case from the E13.5 section in Figure 1F, but it is hard to tell.)

As suggested by reviewer #2, we performed in situ hybridization studies using antisense RNAscope probes to simultaneously label Pomc and Ebf1 expression in coronal sections of E11.5 embryos at the level of the developing basal hypothalamus. Using this additional embryonic time point and plane sectioning (in the original manuscript we had shown sagittal sections of E13.5 mouse embryos), we confirmed that Ebf1 and Pomc co-express in many cells located in the mantle zone of the developing hypothalamus, as can be observed in the newly submitted Figure 2E and Figure 2—figure supplement 1. The previously shown Figure 2F is now included in Figure 2—figure supplement 1. As a consequence of including this new supplement figure, the previous Figure 2—figure supplement 1 is now named Figure 2—figure supplement 2 and its reference changed in the Results section (line 184-195). It is worth mentioning that the level of Ebf1 expression in the developing hypothalamus is much lower than that found in other areas of the same sections, such as the striatal differentiating zone (SDZ) and mesenchyme (MES) as clearly observed in Figure 2E. Our scRNA-seq study shows that the expression level of Ebf1 in POMC cells is highest at E11.5, diminishes at E13.5, is even further reduced at E15.5 and completely disappears from E17.5 and beyond (Figure 2D). These results are confirmed by in situ hybridization (Figure 2E and Figure 2—figure supplement 1). Because Ebf1 has been involved in the control of cell differentiation and migration of early neurons from the subventricular zone to the mantle in the developing striatum (Garel, Marin, Grosschedl, and Charnay, 1999), it is tempting to speculate that Ebf1 plays a similar role in early Pomc-expressing hypothalamic cells. However, our data precludes us from making such a conclusion at this point and further studies will be necessary to evaluate the functional role of Ebf1 during the ontogeny of hypothalamic POMC neurons.

2) Addition of scRNA-seq data from earlier embryonic stages (i.e. E10.5 and/or E11.0) would help to identify the precursor to the population of immature POMC neurons that is captured by cluster 1.POMC-hi/Prdm12.

The reviewer points to a very interesting question which in our hands turned up to be technically more challenging than expected. At the initial design of this study, we sought to FACS cells from E10.5 developing hypothalami to capture the transcriptome at the onset of Pomc expression. However, the number of fluorescent cells in Pomc-dsRed transgenic E10.5 embryos was very small probably due to the existing delay between the onset of Pomc-driven transcription of the transgene and the translation of the early dsRed mRNAs into the fluorescent protein Tdimer2. In our hands, we need at least 17,000 cells to finally obtain 10,000 FACS cells for sequencing. In our current study, primary cells from POMC-TdDsRed/+ mice were first FACS sorted for the transgenic red fluorophore. The sorting percentage for POMC-TdDsRed cells is only 1-2%. Thus, it is very challenging to collect enough cells for FACS followed by scRNA-seq library preparation. We have previously reported using either an anti- ACTH antibody or anti-EGFP antibody in Pomc-EGFP E10.5 embryos that the number of cells expressing proteins from the Pomc gene or Pomc transgenes are very limited at this stage (Nasif et al., 2015, PNAS, Figure 1A and Hael et al., 2020 Mol Metab, Figure 1) (Hael, Rojo, Orquera, Low, and Rubinstein, 2020; Nasif et al., 2015).

3) One question raised by these data are the presence of an inhibitory population of GABA releasing (ie Vgat-expressing) POMC neurons. This is controversial with some studies reporting this and others indicating that there are no GABAergic POMC neurons. Given the evidence for a developmental cluster that is evident even at P12, it is important to resolve this issue in the context of the authors report of a bona fide transcriptional cluster with these co-expression patterns. The authors should examine the colocalization of these Pomc/Vgat neurons using RNA-FISH and report their distribution across development.

Typically, neurons are classified as either glutamatergic based on their expression of one or more vesicular glutamate transporters VGLUT1, VGLUT2 and VGLUT3 mRNAs (Slc17a7, Slc17a6 and Slc17a8, respectively) or GABAergic based on their expression of the vesicular GABA and glycine transporter VGAT (Slc32a1). Alternatively, GABAergic neurons are identified by expression of the synthetic enzymes glutamate decarboxylase 1 or 2 (Gad67 and Gad65, respectively). None of these three last mentioned markers can definitively identify a GABAergic neuron and we do not report the coexpression of Pomc and Slc32a1 in our manuscript. Moreover, triple in situ hybridization studies have demonstrated that subpopulations of POMC neurons express both Slc17a6 and Gad65 and the proportions of these neurons change during early postnatal life in the mouse (Jones et al., 2019; Wittmann, Hrabovszky, and Lechan, 2013). Synaptic storage and release of GABA from neuronal terminals after appropriate stimulation is in our view the definitive means to identify a neuron as GABAergic. Several publications have shown by E-physiology that some subpopulations of POMC neurons induce ipsps that can be blocked by GABA receptor antagonists (Hentges, Otero-Corchon, Pennock, King, and Low, 2009). Although VGAT expression has indeed not been found in POMC neurons, they do express VMAT as shown in Figure 3—figure supplement 1 and VMAT can act as an alternative synaptic GABA vesicular transporter in dopamine neurons (German, Baladi, McFadden, Hanson, and Fleckenstein, 2015).

4) Finally, several conclusions need to be tempered. Nonetheless, the paper will be an outstanding resource for the field. The discussion should be written to reflect the important resource of this data set.

As suggested by the reviewer we shortened some paragraphs of the Discussion to focus on the most important and less speculative matters. We have also tempered some conclusions and highlighted the importance and usefulness of this resource for future studies by many investigators worldwide. All changes in the revised Discussion are indicated in the revised manuscript.

Reviewer #1 (Recommendations for the authors):

There seemed to be a substantial subset of microglia expressing Pomc (Figure 1), did the authors remove these from the dataset? How do they determine a bona fied Pomc neuron vs a neuron with Pomc mRNA contamination.

A general caveat of the scRNA sequencing technique is the potential contamination with RNA molecules from other cells present in the sample or even from the environment. Although this possibility can never be completely ruled out, in our study we used fluorescence assisted cell sorting from medial basal hypothalami collected from POMC-DsRed transgenic mice to enrich the population of cells to be analyzed at the single cell level. In addition, we took several levels of precaution to limit, monitor and discard potential cases of cross contamination events: (1) We applied stringent filtering criteria to exclude doublets and low quality cells as described in the Materials and methods section; (2) We performed an independent unsupervised analysis on cells expressing the transgenic red fluorophore (DsRed) while ignoring Pomc UMI criteria per cell (Figure1—figure supplement 2 and Figure 1—figure supplement 3) and we found similar clustering at each developmental age compared to clustering based on Pomc UMI criteria of > = 1. (3) To validate each cluster, we performed dual RNAscope in situ hybridization to confirm the colocalization of Pomc and feature genes from each cluster (Figure 2, Figure 3 and their linked supplemental figures). The glial cell clusters in this study were defined by showing simultaneous expression of high levels of several marker genes selectively representing microglia, VLMC, or Astrocytes (Figure 1C), and the absence of neuronal markers whereas a reciprocal criterion was used to define neuronal clusters. Because the current study focuses on Pomc expressing neurons, we have only studied and reported neuronal clusters (1-8) in greater detail.

A supplemental model figure/cartoon integrating the various developmental trajectories with their defining transcriptional factors would be useful to visualize the key findings of the paper.

The reviewer’s suggestion is a genuine interest we had for a long time. Unfortunately, we have been unable to capture in a single schematic the design and most meaningful results of this complex work without duplicating graphs already shown in Figures 1, 2, 3 and 6 of this manuscript.

Reviewer #2 (Recommendations for the authors):

• Showing in situ hybridization data for E11.5 markers in a coronal plane, where the ventricular and mantle layer can de distinguished, would help to resolve the question of whether Ebf-hi expression represents newly differentiated neurons that are in the process of migrating out of the ventricular zone. (This appears to be case from the E13.5 section in Figure 1F, but it is hard to tell.)

This question was answered in the Essential Revision 1

• Addition of scRNA-seq data from earlier embryonic stages (i.e. E10.5 and/or E11.0) would help to identify the precursor to the population of immature POMC neurons that is captured by cluster 1.POMC-hi/Prdm12.

This question was answered in the Essential Revision 2.

• The number of total neurons changes dramatically across gestation, in large part because cells that adopt non-POMC fates turn off Pomc expression and are no longer captured by FACS. In Figure 1, reporting the absolute number of neurons in each cluster, instead expressing it as a percentage, would make it easier to interpret.

We appreciate the reviewer’s suggestion and have edited Figure 1 accordingly.

• Analyses shown in Figure 3 Supplemental Figures 3-4 are critical to resolve issues related to POMC progenitor heterogeneity. This information should be moved to the main figures. Maintaining a consistent color scheme between Figure 1 and Figure 3 Supplement 3 would make it much easier to interpret this data. The chart in Supplemental Figure 4 could then remain as a supplemental figure to reinforce this point.

We thank the reviewer’s suggestion and have edited Figure 3 Supplemental Figure3-4 accordingly in this revised version of the manuscript.

• Trajectory analyses shown in Figure 5 should include separate analyses of 2.Pomc-med/Ebf1 and 3.Pomc-med/Nra2, as they appear to map onto different populations in the mature ARH (Supplemental Figure 15).

We agree with the reviewer that the interpretation of original Supplementary Figure15 is not straightforward and may be misleading. Therefore, we decided to remove this comparative analysis from the Results section and its corresponding figures. Based on the reviewer’s comment, the original Figure 5 became Figure 6.

• There are discrepancies between the text (lines 489-497) and Supplementary file 14, which do not link 3.Pomc-med/Nr4a2 neurons to mature POMC populations and Supplementary file 15, which does, making it difficult to interpret these data.

This part has been removed from the Results section.

References

Garel, S., Marin, F., Grosschedl, R., & Charnay, P. (1999). Ebf1 controls early cell differentiation in the embryonic striatum. Development, 126(23), 5285-5294.

German, C. L., Baladi, M. G., McFadden, L. M., Hanson, G. R., & Fleckenstein, A. E. (2015). Regulation of the Dopamine and Vesicular Monoamine Transporters: Pharmacological Targets and Implications for Disease. Pharmacol Rev, 67(4), 1005-1024. doi:10.1124/pr.114.010397

Hael, C. E., Rojo, D., Orquera, D. P., Low, M. J., & Rubinstein, M. (2020). The transcriptional regulator PRDM12 is critical for Pomc expression in the mouse hypothalamus and controlling food intake, adiposity, and body weight. Mol Metab, 34, 43-53. doi:10.1016/j.molmet.2020.01.007

Hentges, S. T., Otero-Corchon, V., Pennock, R. L., King, C. M., & Low, M. J. (2009). Proopiomelanocortin expression in both GABA and glutamate neurons. J Neurosci, 29(43), 13684-13690. doi:10.1523/JNEUROSCI.3770-09.2009

Jones, G. L., Wittmann, G., Yokosawa, E. B., Yu, H., Mercer, A. J., Lechan, R. M., & Low, M. J. (2019). Selective Restoration of Pomc Expression in Glutamatergic POMC Neurons: Evidence for a Dynamic Hypothalamic Neurotransmitter Network. eNeuro, 6(2). doi:10.1523/ENEURO.0400-18.2019

Nasif, S., de Souza, F. S., Gonzalez, L. E., Yamashita, M., Orquera, D. P., Low, M. J., & Rubinstein, M. (2015). Islet 1 specifies the identity of hypothalamic melanocortin neurons and is critical for normal food intake and adiposity in adulthood. Proc Natl Acad Sci U S A, 112(15), E1861-1870. doi:10.1073/pnas.1500672112

Wittmann, G., Hrabovszky, E., & Lechan, R. M. (2013). Distinct glutamatergic and GABAergic subsets of hypothalamic pro-opiomelanocortin neurons revealed by in situ hybridization in male rats and mice. J Comp Neurol, 521(14), 3287-3302. doi:10.1002/cne.23350

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

The revisions by the authors improved the manuscript. These data will be a valuable resource for the field. However, there are a few important points that should be corrected or clarified in a final version.

1. Characterization of Ebf1 expression

The claim that Ebf1 expression at E11.5 is localized to the mantle layer (line 201) and the labels in the accompanying Figure 1E are not accurate. Expression of Pomc and Ebf1 are located within the ventricular zone (VZ), albeit toward the lateral edge, as indicated by the dense DAPI stain. The mantle layer "ML" is label is not accurate; there is almost no "mantle layer" in the mediobasal hypothalamus at E11.5. Also, the differentiating zone indicated at the lateral edge of the VZ in Figure 1E is not striatal.

Note that Figure 2-Figure Suppl 1A shows co-expression of Pomc and Ebf1 at the border between the VZ and ML at E13.5. There appears to be more co-expression on the VZ side than the ML side. The expression pattern shown here is consistent with reports in other regions of the developing CNS, where EBFs have been implicated in regulating migration out of the VZ.

In agreement with the reviewer’s comment, we have modified the sentence (line 186-189) indicating now that the co-expression of Ebf1 and Pomc occurs at the ventricular zone. Accordingly, we have modified the labels in Figure 2E. Concerning the reviewer’s suggestion that Ebf1 may regulate cell migration out of the VZ during neuronal development, we have now include this hypothesis in the Discussion section (lines 514 – 517).

2. Lines 474-5 "The diffusion maps analysis appears to confirm the multiple origins of POMC progenitor cells demonstrated by the RNA velocity study."

This is a confusing sentence because "POMC progenitors" is not the accurate terminology to use. It is easier to just eliminate this sentence.

Following the reviewer’s suggestion we have deleted this sentence in the newly revised manuscript.

3. Discussion of Ebf1+/cluster 2 neurons

Lines 614-6 "The expression profile of cluster 2.Pomc(med)614 /Ebf1 differs from cluster 1.Pomc(high)615 /Prdm12 most particularly at the earliest developmental ages E11.5 to E15.5 and may both represent alternative genetic routes leading to mature POMC neurons."

There is no evidence that the Ebf1+/cluster 2 neurons develop into mature neurons of any sort. The analysis in Figure 6 depicts the Ebf1+/cluster 2 neurons as a dead end in Ray 1 (i.e. retains immature properties). Since the transcriptomics performed here relied on the expression of a Pomc-GFP transgene, it is possible that cluster 2 neurons differentiate into mature non-POMC cell types. The sentence should be corrected by swapping "arcuate" for "POMC".

Following the reviewer’s suggestion, we have replaced “POMC” with “arcuate” in the newly revised manuscript.

4. Discussion of developmental origins of NPY and KNDy neurons Lines 619-626 are confusing.

Timing question: With regard to NPY and KNDy neurons, which are known to derive from Pomc-expressing neurons, the authors argue that they arise "earlier than expected "(presumably at E11.5). And yet, according to the developmental clustering (Figure 4) and developmental trajectory analyses (Figure 6 supplement 1), meaningful numbers of Otp+/cluster 4 (future NPY) and Tac2+/cluster 5 (future KNDy) neuron precursors emerge at E13.5, consistent with reports in the literature.

We thank the reviewer’s constructive suggestion and have revised this sentence accordingly (lines 520 – 523).

We would like to remark that POMC cells from cluster 4 display relatively high levels of Otp already at E11.5. Although Npy and Agrp transcripts only peak up days later, the early presence of Otp at E11.5 is likely an active marker for the differentiation pathway leading to the AGRP phenotype. We would also like to point out that developmental subclustering of POMC cells only at E11.5 revealed a small but distinct group of related cells (E11.5-9) featuring Otp (Figure 4-Source Data2). This subcluster shares high similarity to Pomc(low)/Otp cluster in regards to the gene expression profiles (comparing Figure 4-Source Data 2 to Figure 1-Source Data 2). In situ hybridization studies shown by the Allen Developmental Mouse Brain Atlas also confirmed the expression of Otp in the developing hypothalamus at E11.5 and it is likely these early Pomc(low)/Otp cells are the precursors of AGRP/NPY neurons at later embryonic stages, although this hypothesis requires further investigation.

Developmental lineage: These analyses also support the idea that Otp+/cluster 4 (future NPY) neurons share the same developmental trajectory as Prdm12+/cluster 1 (future POMC) neurons (Figure 6) as well as a common precursor (E13.5-4, Figure 4). This should be explicitly discussed. These studies do not provide clues into the origin of the Tac2+/cluster 5 (future KNDy) neurons.

We thoroughly examined our raw data. We agree that the Pomc(low)/Otp neurons share the common precursors as Pomc (high) / Prdm12 cluster, as shown by the RNA-velocity trajectory analysis in Figure 6A and Figure 6—figure supplement 1B and by the transcriptomic analysis at individual developmental stages (Figure 4—figure supplement 1). We have now addressed this issue in the revised manuscript (Line 528-530).

The reviewer pointed out the existence of common precursor (E13.5-4, Figure 4), indicating the Otp+/cluster 4 (future NPY) neurons share the same developmental trajectory as Prdm12+/cluster 1 (future POMC) neurons. We have carefully re-examined the Figure 4 and its associated supplement figures and source data and confirmed that the gene profiles in E11.5-9 and E13.5-4 subclusters are highly similar to that of Pomc(low)/Otp cluster with the feature genes – Otp, Calcr, and Sst. We also found these two clusters share genes that are highly expressed in Pomc (high) / Prdm12, such as Dlk1, Six3, Six6, Cited1, Isl1 and Peg3. That’s why we consider E11.5-9 and E13.5-4 subclusters as the common origins for both Pomc(low)/Otp and Pomc (high) / Prdm12 clusters (Figure 4—figure supplement 1). It worth noting that AGRP neurons share several similarities with POMC neurons in terms of transcription profiles although their physiological functions are opposite. For example, Tbx3 is expressed in both AGRP and POMC neurons, although its role appears to be different in each cell type based on the phenotype observed in the individual conditional mouse mutants.

https://doi.org/10.7554/eLife.72883.sa2

Article and author information

Author details

  1. Hui Yu

    Department of Molecular and Integrative Physiology, University of Michigan Medical School, Ann Arbor, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing - original draft, Writing - review and editing
    For correspondence
    huiyuz@umich.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5249-0193
  2. Marcelo Rubinstein

    1. Department of Molecular and Integrative Physiology, University of Michigan Medical School, Ann Arbor, United States
    2. Instituto de Investigaciones en Ingeniería Genética y Biología Molecular, Consejo Nacional de Investigaciones Científicas y Técnicas, Buenos Aires, Argentina
    3. Departamento de Fisiología, Biología Molecular y Celular, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina
    Contribution
    Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing - original draft, Writing - review and editing
    For correspondence
    mrubins@dna.uba.ar
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7500-7771
  3. Malcolm J Low

    Department of Molecular and Integrative Physiology, University of Michigan Medical School, Ann Arbor, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Resources, Supervision, Writing - original draft, Writing - review and editing
    For correspondence
    mjlow@umich.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9900-3708

Funding

National Institute of Diabetes and Digestive and Kidney Diseases (DK068400)

  • Malcolm J Low
  • Marcelo Rubinstein

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr Graham L Jones, Dr Zoe Thompson, and Charles Keane for their assistance in animal dissection, tissue sectioning, and data analysis, Dr Joshua Welch and Tongyu Liu for advice on data analysis, the Flow Cytometry Core, the Advanced Genomics Core, the Bioinformatics Core and the Microscopy, Imaging and Cellular Physiology Core (supported by NIH Grant: P30DK020572) and the University of Michigan Center for Gastrointestinal Research (supported by NIH grant: P30DK034933) at University of Michigan for data collection. Funding: This study was supported by NIH grant DK068400 (MJL and MR)

Ethics

All procedures were performed in accordance with the Institutional Animal Care and Use Committee (IACUC) protocol (PRO00008570) at the University of Michigan and followed the Public Health Service guidelines for the humane care and use of experimental animals.

Senior Editor

  1. Marianne E Bronner, California Institute of Technology, United States

Reviewing Editor

  1. Joel K Elmquist, University of Texas Southwestern Medical Center, United States

Publication history

  1. Preprint posted: May 13, 2021 (view preprint)
  2. Received: August 7, 2021
  3. Accepted: January 18, 2022
  4. Accepted Manuscript published: January 19, 2022 (version 1)
  5. Version of Record published: February 1, 2022 (version 2)

Copyright

© 2022, Yu et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,511
    Page views
  • 324
    Downloads
  • 3
    Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.

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)

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

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

  1. Hui Yu
  2. Marcelo Rubinstein
  3. Malcolm J Low
(2022)
Developmental single-cell transcriptomics of hypothalamic POMC neurons reveal the genetic trajectories of multiple neuropeptidergic phenotypes
eLife 11:e72883.
https://doi.org/10.7554/eLife.72883
  1. Further reading

Further reading

    1. Developmental Biology
    Seok Hee Lee, Xiaowei Liu ... Paolo F Rinaudo
    Research Article Updated

    In vitro fertilization (IVF) has resulted in the birth of over 8 million children. Although most IVF-conceived children are healthy, several studies suggest an increased risk of altered growth rate, cardiovascular dysfunction, and glucose intolerance in this population compared to naturally conceived children. However, a clear understanding of how embryonic metabolism is affected by culture condition and how embryos reprogram their metabolism is unknown. Here, we studied oxidative stress and metabolic alteration in blastocysts conceived by natural mating or by IVF and cultured in physiologic (5%) or atmospheric (20%) oxygen. We found that IVF-generated blastocysts manifest increased reactive oxygen species, oxidative damage to DNA/lipid/proteins, and reduction in glutathione. Metabolic analysis revealed IVF-generated blastocysts display decreased mitochondria respiration and increased glycolytic activity suggestive of enhanced Warburg metabolism. These findings were corroborated by altered intracellular and extracellular pH and increased intracellular lactate levels in IVF-generated embryos. Comprehensive proteomic analysis and targeted immunofluorescence showed reduction of lactate dehydrogenase-B and monocarboxylate transporter 1, enzymes involved in lactate metabolism. Importantly, these enzymes remained downregulated in the tissues of adult IVF-conceived mice, suggesting that metabolic alterations in IVF-generated embryos may result in alteration in lactate metabolism. These findings suggest that alterations in lactate metabolism are a likely mechanism involved in genomic reprogramming and could be involved in the developmental origin of health and disease.

    1. Developmental Biology
    2. Neuroscience
    Matthias Blanc, Giovanni Dalmasso ... Cristina Pujades
    Tools and Resources

    Reconstruction of prototypic three-dimensional (3D) atlases at the scale of whole tissues or organs requires specific methods to be developed. We have established a digital 3D-atlas maker (DAMAKER) and built a digital 3D-atlas to monitor the changes in the growth of the neuronal differentiation domain in the zebrafish hindbrain upon time. DAMAKER integrates spatial and temporal data of cell populations, neuronal differentiation and brain morphogenesis, through in vivo imaging techniques paired with image analyses and segmentation tools. First, we generated a 3D-reference from several imaged hindbrains and segmented them using a trainable tool; these were aligned using rigid registration, revealing distribution of neuronal differentiation growth patterns along the axes. Second, we quantified the dynamic growth of the neuronal differentiation domain by in vivo neuronal birthdating experiments. We generated digital neuronal birthdating 3D-maps and revealed that the temporal order of neuronal differentiation prefigured the spatial distribution of neurons in the tissue, with an inner-outer differentiation gradient. Last, we applied it to specific differentiated neuronal populations such as glutamatergic and GABAergic neurons, as proof-of-concept that the digital birthdating 3D-maps could be used as a proxy to infer neuronal birthdate. As this protocol uses open-access tools and algorithms, it can be shared for standardized, accessible, tissue-wide cell population atlas construction.