Abstract
Hypothalamic neural circuits are known to regulate energy homeostasis and feeding behavior, but how these circuits are established during development is not well understood. Here we report that embryonic neural progenitors that express the transcription factor OLIG1 contribute neurons to the ventral hypothalamus including the arcuate nucleus (ARH), a center that regulates feeding behavior. Ablation of bone morphogenetic protein receptor 1a (BMPR1A) in the OLIG1 lineage resulted in hypophagia, hypoglycemia, and weight loss after the second postnatal week with death by week 4. Differentiation and specification of inhibitory hypothalamic neurons contributing to melanocortin and dopaminergic systems were abnormal in the BMPR1A-deficient ARH. Although the hypophagia promoted expression of the orexigenic neuropeptide agouti related protein (AgRP) in the BMPR1A-deficient ARH, there was a profound decrease of AgRP+ axonal terminals in the mutant ARH targets including dorsomedial and paraventricular hypothalamic nuclei. Projection of AgRP+ neurons to these nuclei is known to be regulated by leptin. Leptin injection in neonatal mice increased bone morphogenic protein (BMP) signaling in the ventral hypothalamus, and blocking BMP signaling prevented leptin-induced neurite outgrowth in ARH explant cultures. These findings suggest that BMPR1A signaling is critical for postnatal establishment of leptin-responsive orexigenic fibers from ARH to multiple hypothalamic nuclei. More generally these observations indicate that BMPR1A signaling regulates postnatal establishment of OLIG1 lineage-derived ARH neuronal circuits that are critical for leptin-mediated feeding behavior.
Introduction
The hypothalamus regulates many aspects of feeding behavior and energy homeostasis in response to signals from the periphery. For example, leptin (Zhang et al., 1994), an adipocyte-derived hormone, inhibits food intake and promotes catabolic responses through neural circuits involving multiple nuclei of the ventral hypothalamus. Central to the circuit are two distinct neuronal populations located in the arcuate nucleus of hypothalamus (ARH) that express either proopiomelanocortin (POMC) or agouti-related protein (AgRP)/neuropeptide Y (NPY). AgRP+/NPY+ neurons are orexigenic and inhibited by increased circulating levels of leptin whereas POMC+ neurons are anorexigenic and stimulated by leptin. These ARH neurons project directly to nuclei including the dorsomedial hypothalamus, paraventricular nucleus of the hypothalamus (PVH), and lateral hypothalamus area, which all play important roles in leptin-regulated feeding behavior (Cone, 2005; Gautron and Elmquist, 2011).
Despite advances in the understanding of hypothalamic regulation of feeding behavior, molecular mechanisms involved in development of ARH neural circuits remain limited. During embryonic development, POMC+ and NPY+ ARH neurons are generated from Nkx2.1+ progenitors (Yee et al., 2009), and several other transcription factors have been implicated in hypothalamic neuroendocrine development and function (Caqueret et al., 2005; McNay et al., 2006; Shimogori et al., 2010). Hypothalamic neuronal differentiation is regulated by sequential actions of floor plate-secreted signaling ligands sonic hedgehog (SHH) and bone morphogenetic protein 7 (BMP7) (Ohyama et al., 2005; 2008). BMP signaling is mediated by a heterodimeric complex of type I (BMPR1A, BMPR1B, and ALK2) and type II (BMPR2) receptors, which canonically lead to the phosphorylation and nuclear translocation of the SMAD family effector proteins (Shi and Massagué, 2003). Disruption of signaling through BMPR1A in cells that express the transcription factor OLIG1 causes mutant mice to fail to gain weight normally and die by 4 weeks of age (Samanta et al., 2007). OLIG1 is expressed throughout the ventral neural tube including developing telencephalon and diencephalon, where OLIG1-expressing progenitors differentiate into interneurons and oligodendrocytes (Takebayashi et al., 2000; Lu et al., 2002; Samanta et al., 2007; Mukhopadhyay et al., 2009). Although OLIG1 lineage neurons were not known to populate hypothalamic nuclei, the feeding abnormalities in mice with OLIG1 lineage deletion of BMPR1A suggested to us disruption of hypothalamic function.
Here we show that the OLIG1 lineage gives rise to three neuronal subtypes in the ARH. Ablation of Bmpr1a in the OLIG1 lineage resulted in hypophagia-induced AgRP upregulation and POMC downregulation in ARH at postnatal day 21 (P21). Despite increased expression of AgRP in cell bodies, AgRP+ axon terminals were greatly reduced at the target nuclei. Since projection of ARH neurons is regulated by leptin, this suggested that BMP signaling may be necessary for leptin-mediated axon outgrowth. In support of this hypothesis, we found that leptin injection in neonatal mice increased BMP signaling in the ventral hypothalamus, and blocking BMP signaling with the antagonist noggin prevented leptin-induced neurite outgrowth in ARH explant cultures. These findings suggest that BMPR1A signaling is critical for postnatal establishment of leptin-responsive orexigenic fibers from ARH to multiple hypothalamic nuclei.
Materials and Methods
Generation and maintenance of mouse lines.
OLIG1-cre (OLIG1cre/cre) mice were maintained in a C57BL/6 background while Bmpr1a floxed (Bmpr1afx/fx), Bmpr1a+/−, and Bmpr1b+/− were maintained in a 129Sv:C57BL/6 mixed background. The generation of the OLIG1-cre, Bmpr1afx/fx, Bmpr1b +/−, and Bmpr1b+/− has already been reported (Mukhopadhyay et al., 2009). The OLIG1cre/cre mice were bred to Bmpr1a+/− to get OLIG1cre/+; Bmpr1a+/−, which were then bred to Bmpr1b+/− to generate the OLIG1cre/cre; Bmpr1a+/−; Bmpr1a+/−. The Bmpr1afx/fx animals were bred to Bmpr1b+/− to generate the Bmpr1afx/fx; Bmpr1b+/−, which were then bred to the ROSA lacZ/lacZ (Jackson Laboratories) to generate Bmpr1afx/fx; Bmpr1b+/−; ROSA lacZ/lacZ. In the final mating OLIG1cre/cre; Bmpr1a+/−; Bmpr1b+/− animals were bred to Bmpr1afx/fx; Bmpr1b+/−; ROSA lacZ/lacZ to get OLIG1cre/+; Bmpr1afx/+; Bmpr1b+/−; ROSA lacZ/+ (control), OLIG1cre/+; Bmpr1afx/−; Bmpr1b+/−; ROSA lacZ/+ (R1a KO), OLIG1cre/+; Bmpr1afx/+; Bmpr1b−/−; ROSA lacZ/+ (R1b KO), and OLIG1cre/+; Bmpr1afx/−; Bmpr1b−/−; ROSA lacZ/+ (R1a/b KO) animals. Male mice of each genotype were selected for all postnatal body weight and cellular analyses. The excision of the Bmpr1a floxed allele by the OLIG1-cre expression has already been reported (Samanta et al., 2007). All animal experiments were conducted in conformation to the regulatory standards laid down by the Northwestern University Animal Care and Use Committee.
Immunohistochemistry, in situ hybridization, and EdU labeling.
Mice were transcardially perfused with cold HBSS medium (Invitrogen) and the brains were fixed in 4% paraformaldehyde (PFA) for 2 h. The fixed brains were cryoprotected in 30% sucrose solution overnight and frozen in Tissue-Tek OCT before being cryosectioned for immunohistochemistry. The tissue sections were incubated with primary antibody overnight at 4°C in PBS containing 1% bovine serum albumin (BSA), 0.25% Triton X-100, and 2% goat serum. After a second set of washes the cell and the tissue sections were incubated with fluorophore-conjugated Alexa Fluor goat secondary antibodies (1:500) (Invitrogen) in a solution containing 1% BSA, 0.25% Triton X-100, and DAPI nuclear stain (Invitrogen) for 1 h. After a second set of washes, the slides were mounted in ProLong Gold antifade reagent (Invitrogen). The primary antibodies used were mouse anti NKX2.1 (1:200; Thermo Fisher), rabbit anti-PAX6 (1:2000; Millipore), rabbit anti-OLIG2 (1: 500; Millipore), goat anti-OLIG2 (1:500; R&D), rabbit anti-ASCL1 (1:1000; BD Pharmingen), guinea pig anti-DLX2 (1:3000), rabbit anti-cleaved caspase 3 (1:1000; Cell Signaling Technology), rabbit anti-tyrosine hydroxylase (TH; 1:1000; Pel Freez), mouse anti-TH (1:500; Sigma), rabbit anti-POMC (1:1000; Phoenix Pharmaceuticals), rabbit anti-NPY (1:1000; Phoenix Pharmaceuticals), rabbit anti-AgRP (1:2000; Phoenix Pharmaceuticals), guinea pig anti-AgRP (1:2000; Abnova), rabbit anti-MC4R (1:250; Alomone Labs), sheep anti-α-MSH (1:10,000, Millipore), rabbit anti-β-galactosidase (1:2000, Covance), chick anti-β-galactosidase (1:500, Fitzgerald Scientific), rabbit anti-BMPR1A (Abgent AP2004a, 1:200), rabbit anti-phospho-SMAD1/5/8 (1:300; Cell Signaling Technology), rabbit anti-phospho-STAT3 (1:1000, Cell Signaling Technology). For mRNA expression analysis of POMC and AgRP, plasmids containing full-length POMC and AgRP cDNA were purchased (Open Biosystems) and used to generated digoxigenin-labeled antisense RNA probes to use for in situ hybridization as previously described (Mukhopadhyay et al., 2009). For EdU labeling pregnant dams were give one injection (10 ml/kg) of EdU (Invitrogen) solution (2.5 mg/ml) intraperitoneally. To assess proliferation and cell cycle exit, the embryos were collected 1 or 24 h post-EdU injection, respectively. For phospho-SMAD1/5/8 and phospho-STAT3 staining, recombinant leptin (5 mg/kg; R&D) or vehicle (PBS) is injected intraperitoneally into wild-type mice at P7 and killed at 45 min or 3 h postinjection. All tissue was processed as described above for immunohistochemistry. Immunohistochemistry was followed by staining for EdU using Clik-IT EdU Imaging Kit (Invitrogen) following manufacturer's instructions.
X-gal staining.
Tissue sections were washed in rinse buffer (PBS with 0.01% deoxycholate, 5 mm EGTA, 0.02% NP40, and 2 mm MgCl2) for 10 min twice. After the washes the slides were incubated in staining buffer (PBS with 5 mm K3Fe(CN)6, 5 mm K4Fe(CN)6, 0.01% deoxycholate, 5 mm EGTA, 0.02% NP40, 2 mm MgCl2, and 1 mg/ml X-gal) at 37°C until sufficient staining was detected. The slides were washed with PBS and incubated in primary antibody overnight at 4°C in PBS containing 1% BSA, 0.25% Triton X-100, and 2% goat serum. The following day the slides were washed in PBS and incubated in biotinylated secondary antibodies, which were then peroxidase conjugated as per manufacturer's instructions in VECTASTAIN Elite ABC Reagent Kit (Vector Laboratories). A red stain was developed using the AEC substrate kit for peroxidase (Vector Laboratories).
Blood glucose and serum hormone measurements.
Blood glucose was measured from a drop of blood from the tail vein using Ascensia Contour Blood Glucose Monitoring System (Bayer) as per manufacturer's instructions. For serum hormone measurements, blood samples were left at room temperature for 45 min to allow coagulation and spun at 14,000 rpm for 10 min. The supernatant was collected, aliquoted, and stored at −80°C until further use. Ghrelin, leptin, and insulin concentrations were determined by multiplex ELISA using the Milliplex MAP Mouse Gut Hormone Panel (Millipore) following manufacturer's instructions.
ARH explant culture.
Slice and explant culture was performed based on a modified protocol from the Simerly and Hevner labs (Bouret et al., 2004; Daza et al., 2007). Briefly, brains of P7 C57BL6 pups were dissected, embedded in 4% low melting point agarose, and sectioned in ice-cold KREB buffer to 200-μm-thick slices using a vibratome (Leica VT1000S). The ARH was isolated under a dissecting microscope and placed inside a 30 μl droplet of 4 mg/ml rat tail collagen solution (Millipore) on glass coverslips. Collagen matrix was allowed to set for 30 min before serum ad libitum growth media with either vehicle, leptin (100 ng/ml; R&D 498-OB), or leptin plus Noggin (250 ng/ml; R&D 719-NG) was added to the wells (n = 6–8 per condition). Slices were cultured at 37°C with 5% CO2 for 3 d with fresh media exchanged daily, and then fixed with 4% PFA before immunostaining with β-III tubulin (1:1000; Sigma) and analyzed on a confocal microscope (Leica SP5). For neurite outgrowth quantification, a z-series stack through the entire slice was merged, binarized, and skeletonized using National Institutes of Health (NIH) ImageJ. Per pixel-integrated density of the image field, excluding the explant itself was calculated and compared between experimental conditions.
Cell counts and statistical analysis.
Coronal adult mouse brain sections (10 μm) close to bregma −0.1 mm and bregma +1.7 mm, based on landmarks defined by the Paxinos Atlas of the Mouse Brain, were used for analysis of the rostral hypothalamus (RH) and the tuberal hypothalamus (TuH), respectively. Images of the different hypothalamic nuclei were obtained from four brain sections separated by 40 μm using a Zeiss Axiovert epifluorescence microscope and counted manually using NIH ImageJ software. The ImageJ software was also used to measure the area of the hypothalamic nuclei from which the cells were counted. To study colocalization, images were acquired using a Zeiss UV LSM510 confocal scanning microscope. All counts and image acquisitions were performed blinded. In experiments with only two conditions Student's t test assuming unequal variance was used while for multicondition experiments one-way ANOVA followed by Tukey's post hoc text was used to determine statistical significance. All data are presented as mean ± SEM.
Results
Mice lacking BMPR1A signaling in OLIG1 lineage cells are hypophagic and lean
We previously reported that mice with OLIG1 lineage-specific deletion of Bmpr1a weigh less than littermate controls at P20 and do not survive past P25 (Samanta et al., 2007). To further investigate the role of BMP receptors (BMPRs) in development of OLIG1 lineage cells, we generated mice carrying a lacZ reporter and mutant for Bmpr1a (R1a KO) (OLIG1cre/+;Bmpr1afx/;Bmpr1b+/;ROSAlacZ/+), Bmpr1b (R1b KO), or both (R1a/b KO) (OLIG1cre/+;Bmpr1afx/−;Bmpr1b−/−,ROSAlacZ/+). The body weights of newborn R1a KO mice were comparable to control littermates (OLIG1cre/+;Bmpr1afx/+;Bmpr1b+/−,ROSAlacZ/+) or R1b KO mice until P10 when their growth rates slowed and ultimately turned negative after P16 (Fig. 1C). R1a/b KO mice had only slightly slower growth rates until P10 when they slowed further and paralleled R1a KO mice thereafter (Fig. 1C). At P21, R1a KO and R1a/b KO mice were distinguishable from control littermates by their smaller sizes (Fig. 1A), hypoactivity, and poorly groomed coats. At P21 the body lengths (naso-anal) of R1a KO and R1a/b KO mice were reduced ∼15 and 18%, respectively, compared with control (Fig. 1B; Control: 7.7 ± 0.1 cm; R1a KO: 6.6 ± 0.2 cm; R1b KO: 7.6 ± 0.1 cm; R1a/b KO: 6.4 ± 0.1 cm; p < 0.01), and R1a KO and R1a/b KO mice weighed ∼45 and 55% less than controls (Control: 12.3 ± 0.3 g; R1a KO: 6.8 ± 0.6 g; R1b KO: 12.0 ± 0.4 g; R1a/b KO: 5.6 ± 0.6 g; p < 0.01). Since growth of R1b KO mice was not affected, the weight loss was attributable to loss of Bmpr1a in OLIG1 lineage cells.
The severe anorexia in R1a KO and R1a/b KO mice was accompanied by ∼70 and 56 reductions, respectively, in blood glucose levels at P21 (Fig. 1D; Control: 122.0 ± 4.5 mg/dl; R1a KO: 36.4 ± 7.1 mg/dl; R1b KO: 118.6 ± 7.8 mg/dl; R1a/b KO: 53.9 ± 7.6 mg/dl; p < 0.01). Mutations in TGF-β signaling, which is closely related to BMP signaling, cause hypoglycemia due to hyperinsulinemia (Bouret and Simerly, 2007). However, serum insulin and orexigenic hormone ghrelin concentrations of the mutant mice were comparable to controls (Fig. 1E,F). In contrast, plasma levels of the anorexigenic hormone leptin were reduced ∼70% both in R1a KO and R1a/b KO mice (Control: 5940.2 ± 1194.2 pg/ml; R1a KO: 1960.9 ± 302.4 pg/ml; R1b KO: 6003.1 ± 1954.1 pg/ml; R1a/b KO: 1910.8 ± 81.3 pg/ml; p < 0.05) at P21 (Fig. 1G). To exclude the possibility that R1A deletion in the OLIG1 lineage leads to global reduction of leptin production, we examined plasma leptin levels before the onset of hypophagia at P7, and found no significant differences between R1a KO and control (Fig. 1H). These observations suggest that leptin is produced normally in R1a KO mice, but is subsequently downregulated at P21 as the result of starvation.
Bmpr1a mutant mice exhibit hypophagia-induced neuropeptide alterations in the ARH
Targets of circulating leptin include ARH neurons that form the central melanocortin pathway (Tartaglia et al., 1995; Cheung et al., 1997; Cowley et al., 2001). Subpopulations of ARH neurons upregulate POMC or downregulate NPY/AgRP expression in response to leptin signaling to modulate energy balance (Stephens et al., 1995; Erickson et al., 1996; Thornton et al., 1997). To investigate if hypoglycemia in R1a KO mice reflected inability of ARH neurons to respond to suppressed anorexigenic signals, we first examined POMC and NPY protein expression in ARH. Consistent with lower levels of circulating leptin, we found that anorexigenic peptide POMC immunoreactivity is reduced in R1a KO and R1a/b KO mice at P21 (Fig. 2A–E; 255.2 ± 21.6 cells/mm2 for control; 47.1 ± 16.3 cells/mm2 for R1a KO; 193.5 ± 23.4 cells/mm2 for R1b KO; 86.8 ± 23.0 cells/ mm2 for R1a/b KO; p < 0.01). We also observed a reciprocal increase of soma-localized orexigenic peptide NPY immunoreactivity in R1a KO and R1a/b KO mice at P21 (Fig. 2F–J; 26.4 ± 3.4 cells/mm2 for control; 201.7 ± 37.5 cells/mm2 for R1a KO; 19.7 ± 1.3 cells/mm2 for R1b KO; 168.4 ± 29.4 cells/mm2 for R1a/b KO; p < 0.01). A similar increase is also observed for AgRP peptide expression (data not shown; see Fig. 4H).
Since AgRP and NPY peptides are normally not enriched in the soma, quantification of immunoreactive cell bodies may be biased by detectability. To address this issue, we performed AgRP and POMC mRNA expression analysis by in situ hybridization in R1a KO and littermate control mice at P21. Similar to protein expression analysis, POMC mRNA-expressing cells are significantly reduced in R1a KO mice (Fig. 2P,Q; 205.5 ± 12.8 cells/mm2 for control; 75.5 ± 5.15 cells/mm2 for R1a KO; p < 0.01). More importantly, examination of AgRP mRNA revealed the absolute abundance of AgRP+ neurons and further confirmed an increase in both number and transcript level of AgRP in R1a KO mice (Fig. 2P,Q; 195.5 ± 11.13 cells/mm2 for control; 366 ± 18.43 cells/mm2 for R1a KO; p < 0.01). These changes in melanocortin pathway peptide expression are known responses to hypoglycemia, suggesting that ARH circuit appropriately suppressed anorexigenic peptide expression and upregulated orexigenic peptide expression in response to low leptin levels in R1a KO mice.
Hypothalamic but not midbrain dopaminergic neurons are lost in R1a KO mice
Hyperphagia observed in leptin mutant mice is dependent upon dopamine release (Szczypka et al., 2000). The absence of hyperphagia in response to low leptin levels in R1a KO mice suggested possible dysregulation of dopamine signaling. We therefore analyzed the ARH cell population expressing the catecholamine biosynthetic enzyme, TH, situated dorsal to the region containing POMC+ and NPY+ cells (Fig. 2K), and found ∼58 and 64 reductions in R1a KO and R1a/b KO mice (Fig. 2K–N; 241.6 ± 22.1 cells/mm2 for control; 101.2 ± 29.1 cells/mm2 for R1a KO; 215.2 ± 20.5 cells/mm2 for R1b KO; 87.0 ± 12.8 cells/mm2 for R1a/b KO; p < 0.01). Examination of TH expression in other hypothalamic nuclei also revealed reductions in the anterior preoptic nucleus (APN) and periventricular preoptic nucleus (PPN) in RH (Fig. 3A–E), but no changes in zona incerta (ZI) and dorsomedial hypothalamic nucleus (DMH) in the TuH (Fig. 3K–J). In addition to the hypothalamus, loss of TH activity in dopaminergic neurons in the substantia nigra pars compacta (SNc) causes an aphagic phenotype in mice (Szczypka et al., 2001; Hnasko et al., 2006). However, immunostaining for TH at P21 shows that total dopaminergic neurons in both SNc and ventral tegmental area (VTA) in R1a and R1a/b KO mice did not differ from control (Fig. 3K–R,W). The presence of a ROSAlacZ reporter in the BMPR mutant mice also allowed us to visualize OLIG1-expressing cells and their progeny by β-gal immunohistochemistry. Examination of the OLIG1 lineage contribution to SNc and VTA showed <10% of TH+ neurons in SNc were from OLIG1+ progenitors (Fig. 3S; 9.25 ± 1.27%), and even fewer colabeled cells were observed in VTA (Fig. 3U; 1.35 ± 0.5%). In addition, there were also no significant differences in the number of OLIG1 lineage-derived dopaminergic neurons between control (Fig. 3S,U,M) and mutant mice (Fig. 3U,V,M). Thus it is unlikely that the anorexic phenotype observed here is mediated by midbrain dopaminergic neurons.
OLIG1 lineage contributes to BMPR1A-expressing ARH neuronal subtypes
To determine whether altered BMPR function in ARH neurons contributes to the changes in neuropeptide expression observed in R1a KO mice, we first asked whether OLIG1+ progenitors give rise to neurons in the ARH by double immunostaining for β-gal with markers of the three ARH neuronal subtypes at E14 and P21 in the control mice. We found that a significant number of ARH neurons express β-gal at E14, including neurons that express POMC (Fig. 4A,A′), AgRP (Fig. 4B,B′), and TH (Fig. 4C,C′). Since the majority of ARH neurons are born between E13 and E14 (data not shown), this finding confirms the OLIG1+ origin of ARH neurons. At P21, β-gal expression colocalized with all three ARH neuronal subtype markers (Fig. 4D—F′) but with labeling of smaller percentages of ARH neurons (Fig. 4H; 35.38 ± 3.61% of AgRP+, 36.76 ± 3.93% of POMC+, 21.25 ± 2.8% of TH+ neurons express β-gal). This validates that OLIG1+ progenitor-derived ARH neurons participate in the postnatal ARH feeding circuit, but suggests that there may also be contributions from spatially or temporally distinct progenitors other than the OLIG1+ population. Next we asked whether BMPR1A is expressed by OLIG1+ lineage-derived neurons in the ARH. Immunostaining with a BMPR1A-specific antibody in the control mice revealed a broad distribution of BMPR1A+ cells in the ARH at P21 (Fig. 4G,G′), including the majority of β-gal-expressing OLIG1+ lineage-derived neurons (Fig. 4G,G″, arrows). Since BMP signaling is known to regulate neural progenitor maintenance and differentiation (Gross et al., 1996; Li et al., 1998; Bonaguidi et al., 2008), we investigated whether deletion of BMPR1A alters OLIG1+ lineage contribution to the ARH neurons. Despite an increase of AgRP expression in the R1a KO (Fig. 4I), examination of β-gal+ cells in the mutant and control ARH shows similar numbers and distribution at both E14 and P21 (Fig. 4I, data not shown), suggesting that changes in neuropeptide expression do not result from altered contribution of OLIG1+ lineage to the ARH in the absence of BMPR1A.
Progenitor patterning and proliferation are unaffected in R1a KO mice
To investigate how BMP signaling affects OLIG1-derived ARH neurons during development, we first examined patterning and proliferation of neural progenitors that generate ARH neurons. In TuH, PAX6, OLIG2, and NKX2.1 transcription factors label progressively ventral progenitor domains at E13 (Fig. 5A,B). NKX2.1 is expressed in the most ventral progenitor domain that gives rise to ARH neurons (Yee et al., 2009). Although OLIG1 and OLIG2 expression typically overlap, ARH OLIG1+ lineage cells extend ventrally beyond the OLIG2+ domain and overlap with the NKX2.1+ domain (Fig. 5A,C, inset). These findings suggest that ARH neurons likely differentiate from NKX2.1+/OLIG1+ expressing progenitors. We examined proliferation of progenitor subtypes by colabeling with the thymidine analog, EdU (Fig. 5A,B′), and found that the number of cells double labeled with progenitor markers and EdU in R1a KO mice was comparable to the control (Fig. 5D). To investigate a possible global change in cell cycle exit, we double labeled with EdU and cell proliferation marker Ki67 to identify cells that exited cell cycle within a 24 h period. We found no significant changes in the total or the percentages of EdU+Ki67− cells in the ARH (Fig. 5K). These observations suggest that loss of BMPR1A signaling does not affect progenitor patterning and proliferation and cell cycle exit in the mutant hypothalamus.
Differentiation of specified neuronal precursors is reduced in R1a KO mice
Next we asked whether differentiation of ARH neurons is altered in R1a KO mice. ASCL1, a transcription factor involved in BMP-mediated neuronal specification and differentiation (Lo et al., 1997; 1998; Shou et al., 1999), has been shown to regulate ARH neuron cell fate (McNay et al., 2006). We found at E13 that most ASCL1+ cells within the ARH Nkx2.1+ domain are derived from the OLIG1 lineage, and R1a deletion did not affect the overall distribution of ASCL1+ cells (Fig. 5E,E′). However, there was a significant reduction in ASCL1+ cell numbers in mutants (Fig. 5G,G′,L; 1776.6.16 ± 184.05.34 cells/mm2 vs control 2278.8 ± 121.74 cells/mm2; p = 0.031) but no change in numbers of ASCL1+EdU+ cells (Fig. 5F,F′,L; 536.33 ± 46.72 cells/mm2 vs control 522.41 ± 14.89 cells/mm2; p = 0.8).
There are several possible explanations for the reduction in ASCL1+ cell numbers despite normal numbers of ASCL1+EdU+ cells in R1a KO embryos. One possibility is increased apoptosis of ASCL1+ precursors. However, cleaved-caspase 3 and TUNEL staining were unchanged around the ARH (Fig. 5J,J′,K; data not shown). Alternatively, mutant ASCL1+ precursors may prematurely downregulate ASCL1 expression and adopt a different cell fate. The transcription factor DLX2 is expressed in the developing ARH and maintained in some differentiated neurons (Yee et al., 2009). We found a 23% reduction of total DLX2+ cells in mutant ARH at E13 (Fig. 5H,H′,K; 2597.53 ± 118.8 cells/mm2; littermate controls: 3213.54 ± 137.22 cells/mm2; p = 0.0057). Similar to the ASCL1+ ARH population, numbers of dividing (EdU+ labeled) DLX2+ cells (Fig. 5F,F′,K) were unchanged in the mutants (499.01 ± 14.05 cells/mm2 vs controls 502.81 ± 14.11 cells/mm2, p = 0.21). Interestingly, despite decreases in cells labeled both for DLX2 and ASCL1, there was a small but significant increase in double-labeled cells (Fig. 5M) in the mutant ARH (282.08 ± 45.4 cells/mm2 vs controls 162.54 ± 14.52 cells/mm2, p = 3.92 E06). The majority of DLX2+ASCL1+ cells do not colabel with EdU (Fig. 5M) but reside within the subventricular zone, suggesting that they may be transient precursors that exited cell cycle but failed to further differentiate in the absence of BMPR1A.
ARH neurons in R1a KO mice show altered embryonic cell fate specification
ASCL1 and DLX2 both promote GABAergic cell fate commitment in ventral telencephalon (Anderson et al., 1997; 1999; Letinic et al., 2002; Stühmer et al., 2002; Long et al., 2009), a process that is also known to be regulated by BMP signaling (Yung et al., 2002; Gulacsi and Lillien, 2003; Samanta et al., 2007; Mukhopadhyay et al., 2009). Since DLX2+, AgRP+, and TH+ neurons in the ARH all express GAD67 at birth (Yee et al., 2009), we hypothesized that loss of BMPR1A signaling with the resultant decrease in numbers of ASCL1 and DLX2 neuronal precursors would affect specification of these three ARH neuronal subtypes. We therefore examined the development of DLX2+, AgRP+, and TH+ neurons in the ARH by immunohistochemistry and quantified cell numbers at E14, E17, P0, P7, and P14. We found that in the embryonic ARH, DLX2 and AgRP are expressed by distinct neuronal populations that show opposing changes in cell numbers in R1a KO mice (Fig. 6A,B,E′,F′,H′). Specifically, DLX2+ neurons at E14 are significantly reduced in R1a KO mice (Fig. 7A, E′',H′; 1257 ± 32 cells/mm2; controls: 1545.42 ± 33.11 cells/mm2, p = 4.57 E05), whereas AgRP+ neurons are slightly increased (Fig. 6B,F′,H′; 196.06 ± 11.66 cells/mm2; controls:154.07 ± 13.97cells/mm2; p = 0.046). In contrast, TH+ neurons in E14 ARH are found to always coexpress DLX2, but cell numbers in R1a KO mice are comparable to the control (Fig. 6C,G′,H′; 203.02 ± 27.54 cells/mm2, control: 197.33 ± 6.88 cells/mm2, p = 0.85), suggesting that the observed reduction in DLX2 expression is restricted to TH− neurons. Similar cell number changes are observed for all three neuronal subtypes at E17 (Fig. 6A,B,C; data not shown). These findings confirm that BMPR1A is necessary for the initial differentiation and specification of ARH neuronal subtypes during embryonic development.
ARH neurons in R1a KO mice have different postnatal neuropeptide properties
Postnatally, changes in DLX2+ and AgRP+ neuron numbers observed in R1a KO embryos are no longer present at P0 and P7 (Fig. 6A,B; data not shown). However, a significant increase of DLX2+ neurons is found at P14 in the ARH of R1a KO mice (Fig. 6I; 1252.27 ± 55.24 cells/mm2, control: 694.44 ± 29.27 cells/mm2; p = 2.29 E8), along with a small increase of AgRP+ neurons (Fig. 6J; 197.32 ± 15.18 cells/mm2; control: 137.15 ± 7.59 cells/mm2; p = 0.023). In contrast to their opposing embryonic expression, most AgRP+ neurons at P14 also expressed DLX2 (Fig. 6D,L; 139.45 ± 7.02), and numbers of AgRP+/DLX2+ cells were increased in R1a KO mice (Fig. 6D′,L, 197.32 ± 23.03; p = 0.026). In addition to changes in cell number, we also observed a change in AgRP protein distribution where fibrous staining found in the control is significantly reduced in mutant, accompanied by elevated soma staining intensity (Fig. 6J,J′). Surprisingly, postnatal TH+ neuron numbers were similar in R1a KO and control ARH at all time points analyzed except at P21 when they were reduced (Fig. 6C,K,K′). Although total TH+ neuron numbers were unaffected at P14, the number of TH+/DLX2+ double-labeled cells significantly increased in R1a KO mice (Fig. 6D,L,L′; 225.78 ± 9.16; control: 143.25 ± 11.51; p = 7.72 E05). Since DLX2 function is correlated with GABAergic cell fate specification, increased DLX2 expression in AgRP and TH in R1a KO mice supporting the idea that ARH neurons may acquire more inhibitory characteristics in the absence of BMPR1A signaling.
BMPR1A signaling is necessary for postnatal ARH neurite outgrowth
To investigate how BMP signaling regulates the postnatal establishment of ARH circuits, we focused on developmental events transpiring when weight loss in R1a KO mice is first observed. ARH neurons project to target nuclei during P6 and P12, during which time circulating leptin is elevated and necessary to promote axonal growth (Ahima et al., 1998; Bouret et al., 2004). Since this correlates temporally with the onset of weight loss in R1a KO mice, we asked whether ARH neurons project to appropriate targets in the mutants. Since AgRP+ cell bodies are only found in the ARH, fiber-associated AgRP expression in ARH target nuclei was used as an indicator of ARH-originated neuronal projections. AgRP+ fibers in R1a KO mice (Fig. 7A″) were greatly reduced at P21 in DMH (Fig. 7A–A″), the earliest target of ARH projections, despite increased expression in cell bodies (Fig. 7A, arrows). Similarly, AgRP expression in the PVH (Fig. 7B—B″), an essential component of the hypothalamic feeding circuit, was drastically reduced (Fig. 7B″). TH+ fibers were also reduced in DMH and PVH, consistent with the loss of TH+ cell bodies in R1a KO ARH (Fig. 7A,B). In addition, the loss of POMC+ cell bodies in R1a KO ARH resulted in fewer fibers expressing αMSH in PVH, where the receptor MC4R resides and remains intact in mutants (Fig. 7C—C″). Together, these findings suggest that the leptin-mediated orexigenic signaling which increased AgRP expression in the ARH was unable by itself to enable ARH neurons to project to targets necessary to promote feeding in R1a KO mice.
We therefore asked whether BMP signaling is necessary for ARH neurons to extend neurites in response to leptin. We microdissected the ARH from P7 hypothalamus and cultured it with vehicle, leptin, or leptin plus the BMP antagonist noggin for 72 h. Leptin treatment significantly increased the number of neurites extending from explants (Fig. 7D,E; leptin treatment sample fiber density = 59.28 ± 6.53; vehicle treatment sample fiber density: 12.51 ± 2.05; p = 0.014). However, leptin-induced neurite outgrowth was profoundly inhibited by noggin (fiber density = 15.27 ± 2.81; p = 0.016). Interestingly, when long neurites were occasionally observed in noggin-treated explants, they had markedly reduced branching (Fig. 7D). To examine specifically whether BMPR1A participates the leptin-induced neurite outgrowth, we cultured ARH explants from P7 R1a KO mice and littermate control mice with leptin, and observed a similar reduction of neurite outgrowth in R1a KO ARH explants (Fig. 7F,G; control fiber density = 80.04 ± 11.45; R1a KO fiber density = 26.55 ± 6.29; p = 0.004). These findings suggest that BMP signaling is necessary for the leptin-induced neurite outgrowth from ARH neurons during the onset of establishment of feeding circuits.
Phosphorylation of signal transducer and activator of transcription-3 (STAT3) mediates leptin-regulated adult feeding behavior in the adult (Vaisse et al., 1996; Bates et al., 2003; Xu et al., 2007) as well as postnatal development of ARH circuits (Bouret et al., 2012). Additional signaling cascades, including phosphorylation of mitogen-activated protein kinase and phosphatidylinositol 3 kinase, are responsive to leptin receptor activation (Niswender et al., 2001; Bouret et al., 2012). To determine the effect of leptin on BMP signaling during the period of hypothalamic feeding circuit establishment, we injected leptin into P7 wild-type mice and examined BMP signaling activation via phospho-SMAD 1/5/8 (pSMAD) immunostaining. The efficacy of the leptin injection was confirmed by the increased number of phospho-STAT3 (pSTAT3)-expressing ARH neurons at 45 min postinjection, with many AgRP+ neurons expressing pSTAT3 (data not shown). In the vehicle-injected controls pSMAD-expressing (pSMAD+) cells were found sparsely distributed across the ventral hypothalamus, including the ARH. However, the number of pSMAD+ cells in the ventral hypothalamus, including the ARH, increased significantly at 45 min postleptin injection (Fig. 8A,B; vehicle = 103.72 ± 7.95 cells/mm2; leptin+ 45 min = 158.75 ± 11.57 cells/mm2; p = 4.43 E05). Total pSMAD+ cell numbers returned close to the control levels at 3 h postleptin injection, suggesting that the effect of leptin on overall pSMAD activation is transient. Interestingly, detailed examination of pSMAD expression in AgRP+ neurons in the ARH revealed a small but significant increase at 3 h after leptin injection (Fig. 8A, inset, C; vehicle = 11.38 ± 1.98 cells/mm2; leptin + 3 h = 24.87 ± 3.61 cells/mm2; p = 0.002). These findings suggest that in addition to STAT3 pathway activation, leptin also promotes both short- and long-term activation of BMP signaling during the development of feeding circuits in the ARH.
Discussion
We find that conditional ablation of BMPR1A in the OLIG1 lineage causes severe preweaning anorexia in mice. Genetic mutations previously known to cause anorexia in mice close to weaning include mutations in contactin, anorexic (anx) gene product and Th (Maltais et al., 1984; Zhou and Palmiter, 1995; Berglund et al., 1999). Contactin mutant mice appear anorexic and have increased NPY and AgRP expression in ARH similar to R1a KO mice (Berglund et al., 1999; Fetissov et al., 2005). However, contactin mutant mice exhibit severe ataxia ∼P10 that may lead to anorexia (Berglund et al., 1999). Unlike contactin mutants, R1a KO mice do not exhibit motor deficits, suggesting dissimilar mechanisms underlying the anorexic phenotype. Anx/anx and Th mutant mice exhibit phenotypes more similar to R1a KO. Anx/anx mice begin losing weight around the same time, and show similar reductions in serum leptin levels and changes in hypothalamic NPY and POMC expression (Johansen et al., 2003). Ablation of TH specifically in dopaminergic neurons causes an anorexic phenotype with a time line similar to R1a KO mice (Zhou and Palmiter, 1995), which prompted us to investigate dopaminergic neurons in R1a KO mice. There were no changes in dopaminergic neuron numbers in SNc or VTA, most likely due to the very limited contribution from OLIG1+ lineage to these nuclei.
We also examined the possibility that dopamine dysregulation in the hypothalamus might lead to the anorexic phenotype, as it has been shown that dopamine signaling increases feeding in the medial hypothalamus and inhibits feeding in the lateral hypothalamus (Robert et al., 1990; Fetissov et al., 2000). Based on TH expression, we found that ARH dopaminergic neurons are born and can survive past 2 weeks after birth in R1a KO mice. However, an increased percentage of TH+ neurons at P14 express DLX2, a transcription factor that promotes GABAergic traits. This suggests that despite normal TH expression, alterations in other neurotransmitters in R1a KO mutants could affect how dopaminergic neurons function in neural circuits. We also observed a reduction of ARH TH+ neurons at P21 in R1a KO mice. However, loss of ARH dopaminergic neurons is likely not the cause of hypophagia, since R1a KO mice stop gaining weight at P14 when TH+ neuron numbers were unchanged. Furthermore, loss of ARH dopaminergic neurons has been observed in Dlx1 and Six6 mutant mice, which are both viable into adulthood (Yee et al., 2009; Larder et al., 2011). Thus changes in dopaminergic neurons in ARH could contribute to, but are not primarily responsible for, the anorexic phenotype in R1a KO mice.
Oligodendrocytes and somatic motor neurons arise from the OLIG1+ lineage (Zhou et al., 2000; Lu et al., 2002). Mice mutant for Bmpr1a in the OLIG1+ lineage have increased numbers of oligodendrocytes but no changes in spinal motor neurons (Samanta et al., 2007). The increased number of oligodendrocytes is unlikely to cause a hypophagic phenotype as most oligodendrocytes in postnatal brain are derived from an EMX1+ lineage (Kessaris et al., 2006) and mice carrying Bmpr1a mutations in the EMX1+ lineage are not hypophagic and are viable (Yuhki et al., 2004). In hindbrain, the OLIG1+ lineage contributes to hypoglossal neurons that innervate the tongue, but we found no changes in neuron numbers in the hypoglossal nucleus (data not shown). Further, mice with deletion of Bmpr1a in the HB9+ lineage that generates all somatic motor neurons are viable into adulthood (E. Laufer, personal communication). Thus loss of BMPR1A function in somatic motor neurons, including hypoglossal neurons, does not underlie the hypophagia.
Cholinergic neurons of the basal forebrain (BFCN) are also generated from OLIG1/2+ lineage precursors and cholinergic neuron differentiation is affected in OLIG2 mutant brains (Furusho et al., 2006). Moreover, cholinergic differentiation is induced by BMP9 (López-Coviella et al., 2000). However, there were no significant changes in BFCN differentiation in R1a KO animals (data not shown). This likely reflects the presence of intact ALK1, the major type I receptor for BMP9 (Upton et al., 2009). Thus the hypophagia and weight loss in R1a KO mice does not involve disruption of cholinergic neurons. OLIG1 progenitors also generate cortical interneurons (Samanta et al., 2007; Mukhopadhyay et al., 2009). Although unlikely, changes in forebrain interneurons could contribute to the hypophagia observed in R1a KO mice.
BMP7 is required for specification of NKX2.1+ progenitor-derived arcuate neurons (Ohyama et al., 2005). Consistent with reports that BMP7 acts on postmitotic cells that are ventralized by SHH to induce dopaminergic differentiation, we observed no changes in progenitor proliferation in R1a KO animals. Our observation of reduced ASCL1 and DLX2 expression in ARH suggests a role for BMPR1A signaling in neuronal differentiation. However, the reduced number of DLX2+ cells in embryos does not lead to reduced numbers of postnatal DLX2+ neurons. It is possible that reduced numbers of DLX2+ neurons in R1a KO mice are still functionally sufficient after birth. Alternatively, different circuit properties resulting from DLX2+ neuron reduction may elicit compensatory changes, including recruitment of other ARH neurons to become DLX2+/GABAergic. This hypothesis is supported by our observation that DLX2+ neuron numbers increase at P14 and P21. However, it is unclear whether this increase is a result of hypophagia or contributes to the behavioral phenotype. GABAergic signaling has recently been shown to be critical for the orexigenic function of AgRP neurons(Wu et al., 2009; Vong et al., 2011), and the observed DLX2 expression increase correlates well with the need to increase GABAergic/orexigenic signaling as result of hypophagia.
Nevertheless, the profound reduction of AgRP+ terminals in ARH target nuclei of R1a KO mice likely prevented the transduction of orexigenic signals. Our in vitro findings strongly suggest that this reduction is the result of defective neurite outgrowth of ARH neurons, as inhibition of BMP signaling by noggin or R1A deletion prevented leptin-induced neurite outgrowth of cultured ARH explants. It has been suggested that leptin preferentially activates pSTAT3 and increases the density of orexigenic AgRP+ fibers during postnatal development (Blackshaw et al., 2010; Bouret et al., 2012). These findings raise the possibility that reduced neurite outgrowth in ARH explants observed with BMP signaling inhibition might be specific to AgRP+ neurons and contributes to hypophagia. However, postnatal ablation of AgRP neurons has been shown to not affect feeding behavior in the adult (Luquet et al., 2005), suggesting that compensatory mechanisms are likely in place to relay orexigenic signals originally mediated by the lost AgRP neurons. The Olig1 lineage also contributes to other hypothalamic nuclei including the ventral medial hypothalamus (VMH), which is known to participate in leptin-regulated feeding behavior(Dhillon et al., 2006; Bingham et al., 2008; Zhang et al., 2008). Our observation of leptin-induced BMP signaling activation in the VMH suggests BMPR1A may also regulate the development of VMH-mediated circuits, which in turn contribute to hypophagia observed in R1a KO mice.
Increased SMAD1/5/8 phosphorylation in ventral hypothalamus 45 min after leptin injection suggests that BMPR activation rapidly responds to elevated leptin concentrations. The leptin receptor has been shown to physically interact with LRP1 (Liu et al., 2011), a low-density lipoprotein receptor-related protein that functions as a TGF-β coreceptor in the vascular system (Boucher et al., 2007). However, direct interaction between leptin receptor and BMPR1A has not been reported. In the ARH, AgRP+ neurons show a small but significant increase of pSMAD1/5/8 3 h after leptin injection, following the increase of pSTAT3 observed at 45 min after leptin injection. Sequential activation of pSTAT3 and pSMAD has been reported during astrogliogenesis, during which leukemia inhibitory factor activated pSTAT3 promotes BMP2 expression that leads to commitment of neural stem cells to the astrocyte lineage (Fukuda et al., 2007). Further examination of pSTAT3 transcriptional targets in the ARH neurons would help address the possibility that pSTAT3 and pSMAD signaling are linked.
There are precedents for BMP promoting neurite outgrowth (Li et al., 1998; Zhang et al., 1998; Iwasaki et al., 1999). Developmentally, BMP7 regulates the guidance and extension of commissural axons in embryonic spinal cord (Butler and Dodd, 2003) Unlike our observation in the ARH, these effects of BMP7 are mediated by BMPR1B instead of BMPR1A (Yamauchi et al., 2008). Expression of BMP7 is present in the embryonic ventral hypothalamus, and is found postnatally in meningeal cells adjacent to the ARH during leptin-dependent neurite outgrowth (Söderström and Ebendal, 1999; Ohyama et al., 2008; Choe et al., 2012). These findings suggest that BMP7 may act as the ligand for BMPR1A-expressing ARH neurons. Our study again highlighted the functional diversity of type 1 BMPR subtypes, which has been shown to perform distinct physiological roles in other systems (Sahni et al., 2010). BMPR II, through SMAD independent regulation of LIM kinase and Cofilin activation, has also been implicated in axon turning and extension (Wen et al., 2007; Phan et al., 2010). While it is currently unclear what downstream target of BMPR1A is responsible for the observed phenotype, our findings clearly demonstrate that BMPR1A signaling is essential for the perinatal development and establishment of ARH neuronal circuits.
Footnotes
This project was supported by National Institutes of Health Grants NS 20013 and NS 20778. We thank Yuji Mishina for providing us the Bmpr1a, David Rowitch for the OLIG1-cre, and Karen Lyons and Henry Kronenberg for the Bmpr1b mouse lines. We thank Lixin Kan and Jayshree Samanta for technical assistance and manuscript discussion.
The authors declare no competing financial interests.
- Correspondence should be addressed to Chian-Yu Peng, Department of Neurology, Northwestern University's Feinberg School of Medicine, Chicago, IL 60611. c-peng{at}northwestern.edu