Adult skull bone marrow is an expanding and resilient haematopoietic reservoir

feature-image

Play all audios:

Loading...

ABSTRACT The bone marrow microenvironment is a critical regulator of haematopoietic stem cell self-renewal and fate1. Although it is appreciated that ageing, chronic inflammation and other


insults compromise bone marrow function and thereby negatively affect haematopoiesis2, it is not known whether different bone compartments exhibit distinct microenvironmental properties and


functional resilience. Here we use imaging, pharmacological approaches and mouse genetics to uncover specialized properties of bone marrow in adult and ageing skull. Specifically, we show


that the skull bone marrow undergoes lifelong expansion involving vascular growth, which results in an increasing contribution to total haematopoietic output. Furthermore, skull is largely


protected against major hallmarks of ageing, including upregulation of pro-inflammatory cytokines, adipogenesis and loss of vascular integrity. Conspicuous rapid and dynamic changes to the


skull vasculature and bone marrow are induced by physiological alterations, namely pregnancy, but also pathological challenges, such as stroke and experimental chronic myeloid leukaemia.


These responses are highly distinct from femur, the most extensively studied bone marrow compartment. We propose that skull harbours a protected and dynamically expanding bone marrow


microenvironment, which is relevant for experimental studies and, potentially, for clinical treatments in humans. SIMILAR CONTENT BEING VIEWED BY OTHERS RESILIENT ANATOMY AND LOCAL


PLASTICITY OF NAIVE AND STRESS HAEMATOPOIESIS Article Open access 20 March 2024 SKULL BONE MARROW CHANNELS AS IMMUNE GATEWAYS TO THE CENTRAL NERVOUS SYSTEM Article 23 November 2023 SKULL


BONE MARROW AND SKULL MENINGES CHANNELS: REDEFINING THE LANDSCAPE OF CENTRAL NERVOUS SYSTEM IMMUNE SURVEILLANCE Article Open access 29 January 2025 MAIN Niche microenvironments composed of


multiple cell types and molecular signals regulate the self-renewal and the ultimate fate of haematopoietic stem cells (HSCs) in bone marrow1,3 (BM). The skull, a flat bone encasing the


brain, develops via intramembranous ossification and fundamentally differs in both composition and remodelling from long bones, which undergo endochondral ossification during development and


sustain load-bearing mechanical stress4,5,6. The thin and flat mouse calvarium, the main component of the skull roof, has been successfully used as a platform for the imaging of HSC


dynamics7,8. Potential functional differences between BM compartments, such as those in long and flat bone, however, remain little understood. It has been shown that HSCs are homogenously


distributed throughout calvarium and long bone (namely diaphysis and epiphysis) under homeostatic conditions, whereas chimerism in mice transplanted with congenic haematopoietic stem and


progenitor cells (HSPCs) is initially higher in calvarium and epiphysis relative to diaphysis9. The discovery of designated channels between the calvarial BM and the dura mater of the


meninges indicates that skull BM acts as a haematopoietic reservoir for meningeal and central nervous system immunity10,11. Ageing causes an overall decline in haematopoietic function, with


compromised lymphopoiesis and a bias toward dysfunctional myelopoiesis2,12. Again, it remains unknown whether all BM compartments are equally affected by ageing processes and therefore


collectively contribute to haematopoietic decline. Here we show that calvarial BM is fundamentally different to the BM in long bone. Calvarial BM continuously expands during adulthood, is


resistant to most major hallmarks of BM niche ageing, and gradually increases its systemic haematopoietic contribution throughout adult life and ageing. Our work also establishes that BM


compartments exhibit unexpected functional differences, with implications for a range of physiological conditions such as pregnancy and ageing, as well as pharmacological treatments and


disease settings. SKULL BM EXPANDS DURING ADULT LIFE To characterize changes in calvarial BM throughout adult life, we collected skulls from young adult (10–14 weeks; hereafter referred to


as ‘young’), middle-aged (31–37 weeks), old (52–75 weeks) and geriatric (95+ weeks) mice. Initial stereoscopic observation showed that blood-filled structures, representing blood vessels and


surrounding BM, occupy only limited areas of the frontal and parietal bone in young adults, whereas interparietal bone is mostly filled by blood cells (Fig. 1a). Although angiogenic growth


of blood vessels is normally confined to development, regeneration or certain pathological conditions13,14, substantial expansion of the vasculature can be seen in calvaria from middle-aged


mice (Fig. 1a–d). Further increases in vascular area and vessel diameter can be observed in old and geriatric skulls (Fig. 1a–d). Flow cytometric analyses show that the expansion of


calvarial vasculature is accompanied by a profound increase in total haematopoietic cells, HSCs, HSPCs, committed haematopoietic progenitors and stromal cells (Supplementary Data 1; gating


strategies shown in Supplementary Data 2). Skull cross-sections from young and geriatric mice show that the enlargement of calvarial BM occurs more rapidly in female mice than in males (Fig.


1a and Extended Data Fig. 1a,b). Expansion of calvarial BM results in an increase in total skull thickness (of approximately 70% between young and geriatric mice of both sexes) without


decreasing the thickness of the cortical bone tables (Extended Data Fig. 1a,b). Immunofluorescence analysis of skull sections shows a profound increase in V-type proton ATPase-positive


activated osteoclasts lining the surface of the inner and outer cortical bone tables in skulls from old relative to young mice (Extended Data Fig. 1c,d), indicating that bone resorption by


osteoclasts contributes to BM expansion. To address whether skull BM expansion also occurs in ageing humans, we analysed computed tomography (CT) head scan images obtained from 36 human


patients evaluated for small cerebral aneurysms. This revealed significantly increased area enclosed by the inner and outer cortical bone tables in older (61–69 years of age) compared with


younger (21–40 years of age) individuals (Extended Data Fig. 2a–c). A larger difference was seen in female skulls (83.2% increase versus 24.6% in males), showing that this sexual dichotomy


occurs in both mice and humans. These results demonstrate a substantial expansion of the BM volume in both mouse and human skull during adulthood and ageing. VESSEL GROWTH DURING SKULL BM


EXPANSION Previous studies have shown that arterial and sinusoidal endothelial cells are important components of the HSC microenvironment15,16. To gain more insight into the organization of


the calvarial vasculature, we first utilized an in vivo immunofluorescence staining method with CD31 antibody in _Flk1_-GFP knock-in mice17 (_Flk1_ is also known as _Kdr_). This approach


revealed profound differences in vascular architecture in young adult frontal, parietal and interparietal skull bones (Extended Data Fig. 3a,c). Specifically, larger FLK1+CD31+ sinusoidal


vessels with distal sprouts dominate in frontal bone, thinner FLK1+CD31+ sinusoidal vessel loops connect to FLK1low/−CD31+ arterioles in parietal bone, and interparietal bone is completely


filled by a dense network of sinusoidal vessels and arterioles (Extended Data Fig. 3a,c). Vascular leakage assays with Evans Blue corroborated local differences in vascular properties in


different skull parts and revealed significantly higher leakage in frontal bone compared with the other two skull bones (Extended Data Fig. 3b,d). To gain more insight into the organization


of the calvarial vasculature and the changes during adulthood and ageing, we performed in vivo labelling with fluorescent CD31 and endomucin antibodies. In young adult mice, only


interparietal bone was largely filled by a dense vascular network, whereas substantial and sustained vascular expansion was observed in frontal and parietal bones of middle-aged and old mice


(Fig. 1b–d). Vascular area increased by 6.4- and 3.6-fold in frontal and parietal bones, respectively, between young and geriatric mice. Endomucin, which decorates calvarial sinusoidal


vessels but not arteries and smaller arterioles, has been implicated in bone maintenance and HSC function18,19. Vascular endomucin expression significantly increased until the geriatric


stages (3.1-fold increase in frontal bone and 7.0-fold increase in parietal bone). We next focused on parietal bone, which showed the largest BM expansion and a transition from thin looping


sinusoidal vessels in young adult mice to a dense network at later stages (Fig. 1b,c). Sinus diameter increased 4.8-fold and dorsal–ventral vertical extension of the vascular network can be


observed at the old stage (Fig. 1c,d and Extended Data Fig. 3e). Immunofluorescence staining of caveolin-119,20 showed a 2.0-fold increase of interconnecting arterioles during ageing


(Extended Data Fig. 3f,g). In vivo immunofluorescence staining of haematopoietic cells using fluorescent CD45 antibody and flow cytometric analyses show that vascular growth is associated


with a substantial increase in haematopoietic cells (13.1-fold increase between young and geriatric mice) within the expanding BM (Extended Data Fig. 4a,b and Supplementary Data 1). CD45+


cells were predominantly associated with large-calibre sinusoidal vessels, which were highly abundant in old and geriatric skull (Extended Data Fig. 4a). The link between vessel calibre and


BM was further confirmed by a genetic approach; labelling of haematopoietic cells in the _Vav-cre;Rosa26-mTmG_ background21,22 shows that haematopoietic cells expressing GFP surround


large-calibre sinusoidal vessels (Extended Data Fig. 4c). Similarly, CD3e+ T lymphocytes, B220+ B lymphocytes and CD11b+ cells (including monocytes/macrophages, granulocytes and natural


killer cells) are associated with large-calibre sinusoids in skull BM (Extended Data Fig. 4d). With the exception of interparietal bone, BM-associated sinusoidal vessels with strong


endomucin immunostaining are comparably sparse in young adult mice and increase substantially throughout adult life and ageing (Fig. 1b–d). Interparietal bone, which is largely filled by


vessels and BM at the young adult stage, shows minimal age-related changes in vascular area (24.3% increase between young and geriatric samples) and sinus diameter (9.3% decrease between


young and geriatric samples) but a significant increase in endomucin expression (2.0-fold increase between young and geriatric samples) (Fig. 1b,d). By contrast, vascular density in the


metaphyseal region of the femoral BM decreased by 45.2% between the middle-aged and old stages, but rebounded at geriatric stages with the emergence of very thin vessels displaying high


endomucin expression (Extended Data Fig. 5a,b), which have been recently characterized as dysfunctional sinusoidal vessels23. Electron micrographs of the luminal surface of endothelial cells


in young and old skull BM show comparable patterns of endothelial fenestrations, which are a hallmark of vessels in primary lymphoid organs including BM. By contrast, endothelial cells in


old femur BM show a highly irregular pattern of fenestrae and the emergence of larger pores in the endothelial surface (Extended Data Fig. 5c). In contrast to the increasing vasculature of


skull marrow, vessel density in the dura mater, which is adjacent to the skull and connected directly to the skull BM via specialized channels10, decreases significantly during ageing


(Extended Data Fig. 5d,e). Collectively, these results show continuous expansion of the BM and vasculature in the adult and ageing calvarium but not in the adjacent dura mater. Moreover,


ageing femoral BM displays features of compromised vascular integrity, which are not seen in ageing skull BM. SKULL BM RESPONSES IN PATHOPHYSIOLOGY High demand for certain blood cells leads


to strongly increased haematopoiesis in a range of physiological but also pathological conditions1,24. Pregnancy, for example, induces extramedullary haematopoiesis and the rapid expansion


of maternal blood volume and red blood cells25,26. Skull BM and associated vasculature increased substantially during pregnancy, namely 1.8-fold at 17 days postcoitum (dpc) and 2.3-fold at 2


days postpartum (dpp), which involves substantial vessel enlargement in frontal and parietal bone together with a strong increase in total haematopoietic cells, HSCs and HSPCs, committed


haematopoietic progenitors and stromal cells (Fig. 2a,b and Supplementary Data 1). As reported previously27, CD31hiendomucinhi type H vessels in femur increased strongly at late-stage


pregnancy and decreased postpartum, but overall vascular density in the metaphyseal region showed more limited changes (26.2% increase at 17 dpc and 8.5% decrease at 2 dpp) relative to skull


(Extended Data Fig. 6a,b). Although the total number of haematopoietic cells in femur was not significantly changed by pregnancy, flow cytometry showed increases in HSCs, HSPCs, certain


committed haematopoietic progenitors and stromal cells (Supplementary Data 1), consistent with enhanced blood cell production. Ischaemic stroke in a transient mid-cerebral artery occlusion


model (tMCAO) was shown to activate BM haematopoietic stem cells28. Consistent with this finding, tMCAO induced a striking expansion of skull BM and associated vasculature (2.1-fold


increase), filling a substantial portion of parietal bone within 7 days (Fig. 2c,d). These changes are accompanied by strong increases in calvarial total haematopoietic cells, HSCs, HSPCs,


committed haematopoietic progenitors and stromal cells (Supplementary Data 1). By contrast, vascular density in femoral BM decreased by 26.0% and haematopoietic cells and subpopulations were


either unchanged or reduced in the stroke condition (Extended Data Fig. 6c,d and Supplementary Data 1). Haematological malignancies are known to lead to BM alterations, which, in turn,


facilitate disease progression29. A chronic myeloid leukaemia (CML) model showed a modest but significant expansion of the skull BM and associated blood vessels (39.1% increase), which


involved an increase in haematopoietic and stromal cells (Fig. 2e,f and Supplementary Data 1). By contrast, substantial loss of vascularity (50.4% decrease), severe vascular rarefaction and


limited changes in total haematopoietic cell number were observed in CML femoral BM (Extended Data Fig. 6e,f and Supplementary Data 1). Biologically active fragments of parathyroid hormone


(PTH) are used to treat osteoporosis30, and also lead to HSPC expansion31. We found that 4 weeks of daily PTH treatment resulted in substantial skull BM expansion (4.6-fold increase),


increased large-calibre vessels and expansion of vessel-associated BM area (Fig. 2g,h). Flow cytometry confirmed the increase in total haematopoietic cell number, HSCs, HSPCs and certain


committed progenitors, as well as in endothelial and stromal cells (Supplementary Data 1). By contrast, only very limited changes were observed in femoral BM immunostaining and flow


cytometry (Extended Data Fig. 6g,h and Supplementary Data 1). These results demonstrate an unexpected dynamic plasticity of skull BM in a range of pathophysiological contexts and uncover


differences between calvarial and femoral BM. SKULL HSPCS DRIVE BM GROWTH DURING AGEING HSPCs residing in BM change drastically during ageing, both in quantity and differentiation


potential32. To address whether such cells might promote skull BM expansion, we isolated lineage-negative (Lin−) cells from young or old skull and femur and transplanted them into lethally


irradiated young recipients (Fig. 3a). In vivo immunofluorescence staining showed substantial expansion (2.2-fold) of skull BM vessels in recipient mice that received cells from old donor


skull BM compared with those that received cells from young donors (Fig. 3b,c). Transplantation of old donor femur BM also resulted in the expansion of skull BM vessels in recipient mice but


to a lesser degree than old donor skull BM. Notably, vascular density in femur was not significantly different between recipients of transplants from young and old donors (Extended Data


Fig. 7a,b), further suggesting divergent responsiveness of different BM compartments. Flow cytometric analyses of HSPCs showed that both long-term HSCs (LT-HSCs) and Lin−Sca1+KIT+ (LSK)


HSPCs increase in number in the ageing skull BM (Fig. 3d,e), similar to reports for ageing long bone33,34. To further test whether HSPCs contribute to BM expansion and vascular growth in


skull, we treated mice with prostaglandin E235,36 (PGE2) or AMD3100 (ref. 37) to modulate HSPC number. PGE2 treatment increased the number of HSPCs in the skull BM by 92.3% (Extended Data


Fig. 7c,d), resulting in significant expansion of calvarial vasculature and vessel-associated BM relative to control mice (Fig. 3f,g and Supplementary Data 1). AMD3100 treatment mobilized


HSPCs into the peripheral blood by 11.2-fold (Extended Data Fig. 7e,f) and induced a 60.0% reduction of skull BM and vasculature relative to control mice (Fig. 3h,i). Although both PGE2 and


AMD3100 treatments led to changes in vascular density in femoral BM, these alterations were confined to existing marrow without changing the BM area (Extended Data Fig. 7g–j). These results


indicate that treatments acting on HSPCs can dynamically modulate the skull BM and its vasculature. VEGFA DRIVES SKULL BM AND VESSEL GROWTH VEGFA is a master regulator of vascular


growth38,39, and enzyme-linked immunosorbent assays (ELISAs) showed that levels of the growth factor were significantly increased in the skull BM in response to pregnancy, stroke and PTH


treatment, whereas no changes could be detected in the femoral BM under the same conditions (Extended Data Fig. 8a). VEGFA was significantly reduced in the femoral BM of mice with CML but


was increased in the skulls of the same mice (Extended Data Fig. 8a). Previous studies showed that HSPCs isolated from long bone express high levels of _Vegfa_ transcripts compared with


other BM cell types40. We performed quantitative PCR with reverse transcription (RT–qPCR) analyses of various haematopoietic and stromal cell populations after fluorescence-activated cell


sorting (FACS) from enzymatically digested BM. Lin−Sca1+KIT+ HSPCs from skull indeed showed higher _Vegfa_ expression compared with Lin−Sca1−KIT+ cells (4.72-fold), Lin−Sca1−KIT− stromal


cells (5.97-fold) or Lin+ mature haematopoietic cells (5.92-fold) in young mice (Fig. 4a). HSPCs isolated from aged skull show higher _Vegfa_ expression compared with Lin−Sca1−KIT+ cells


(8.39-fold), Lin−Sca1−KIT− stromal cells (10.59-fold) and Lin+ mature haematopoietic cells (9.86-fold) (Fig. 4a). By contrast, HSPCs isolated from aged femur show significantly lower _Vegfa_


expression (5.5-fold) compared with the same population from young mice (Extended Data Fig. 9a). These results are mirrored by the expression of transcripts encoding VEGFA120 and VEGFA164


isoforms, which are highest in HSPCs from old skull, but also substantially higher in HSPCs from young mice compared with other cell populations (Extended Data Fig. 9b). At the protein


level, ELISA of whole-BM lysates shows that VEGFA levels are initially highest in young femur, probably reflecting expression by non-haematopoietic sources such as chondrocytes, and decrease


by 40.2% during ageing (Fig. 4b). By contrast, the amount of VEGFA in calvarial BM increases 2.2-fold from young to geriatric stages (Fig. 4b). Immunofluorescence of skull sections confirms


high expression of VEGFA by resident KIT+ HSPCs (Fig. 4c). As hypoxia is a critical regulator of VEGFA expression in angiogenesis41,42, we administered Hypoxyprobe (pimonidazole


hydrochloride) to analyse age-related changes in the oxygenation of skull and femur (Extended Data Fig. 9c,d). The Hypoxyprobe signal was significantly higher in skull BM relative to femoral


BM at old (3.3-fold) and geriatric stages (4.8-fold) (Fig. 4d,e and Extended Data Fig. 9c,d), consistent with the observed increase of VEGFA in ageing skull. Next, we manipulated


VEGFA-dependent VEGFR2 signalling in vivo to directly address the role of this pathway in the expansion of skull vasculature and BM. For gain-of-function experiments, a cDNA encoding a


bone-homing version of VEGFA (Methods) was cloned into the vector pLIVE, which enables constitutive protein expression in liver after hydrodynamic tail vein injection43. _Vegfa_


overexpression induced a 2.5-fold increase of skull BM, including HSCs and other haematopoietic and stromal cells, and associated vasculature (Fig. 4f,g and Supplementary Data 1).


Conversely, inhibition of VEGFR2 signalling with the blocking antibody DC101 inhibited BM expansion and vascular growth during ageing, as shown by immunostaining and flow cytometry (Fig.


4h,i and Supplementary Data 1). By contrast, although we observed vascular density changes within the confines of the femoral BM, these alterations did not involve any substantial volumetric


changes of the BM as in the skull (Extended Data Fig. 9e–j). These results demonstrate that VEGFA–VEGFR2 signalling regulates skull BM expansion during ageing and that the increase of HSPCs


in ageing skull is likely to be a relevant source of VEGFA. RESILIENCE OF AGEING SKULL BM The ageing BM microenvironment disrupts normal self-renewal and other essential functions of


HSCs44. To assess whether hallmarks of ageing affect calvarial and femoral BM equally, we characterized adipogenesis and inflammation in both compartments in young and geriatric mice.


Adipocytes have previously been reported to accumulate in the ageing BM and negatively regulate HSC function45. Whole-mount or cross-sectional staining of skull or femur, respectively, with


BODIPY shows minimal changes in the number of lipid-filled mature adipocytes in the skull, whereas a substantial 5.9-fold increase in adipocytes can be observed in the metaphyseal region of


the ageing femoral BM (Fig. 5a,b). Next, we used a bead-based multiplex array to quantify 12 pro-inflammatory cytokines in lysates from skull and femoral BM isolated from young and old mice.


Of note, only IFNγ was significantly upregulated in old skull, whereas 7 out of 12 cytokines, including TNF and IL-6, showed increased expression in old femur relative to young femur (Fig.


5c). It should also be noted that the expression of most pro-inflammatory cytokines in old femur was substantially higher than in skull, further supporting the notion that detrimental


processes associated with ageing are more prominent in femur. Ageing HSCs exhibit a myeloid bias at the expense of lymphopoiesis46. To assess whether such myeloid bias occurs and determine


which stages of the myeloid differentiation cascade are affected in the ageing skull, we analysed myeloid-committed progenitors and mature myeloid cells in the geriatric skull or femoral BM


using flow cytometry. We detected significant proportional increases of Lin−Sca1−KIT+ multipotent, common myeloid and granulocyte-monocyte progenitors, but not megakaryocyte-erythroid


progenitors, in the femur compared with the skull (Fig. 5d). Furthermore, RT–qPCR analyses of the key myeloid determination factors _Cebpa_ and _Spi1_ (which encodes the transcription factor


PU.1) show substantial upregulation in the geriatric femur BM versus skull BM (Fig. 5e). Finally, to compare the myeloid output from skull versus hindlimb BM into peripheral blood, we


selectively shielded these areas during lethal irradiation of young and old mice (Fig. 5f). Whereas there was no difference in myelopoiesis between head or leg shielding in young mice,


peripheral blood output from hindlimbs of old mice showed a strong myeloid bias compared with skull (Fig. 5f,g). Together, these results show that calvarial BM is remarkably resilient


against key features of ageing and thereby maintains a healthy microenvironment for haematopoiesis even in old mice. OLD SKULL PRESERVES FUNCTIONAL BM NICHES It has been established that


skull BM contains functional HSCs that can fully reconstitute the entire haematopoietic system in a transplantation setting9 and is the source of meningeal immune cells11,47. The effect of


ageing on haematopoietic output by different BM compartments, however, remains unknown. Following lethal irradiation in a non-transplantation setting, all mice with whole-body exposure died


within 12 days of BM failure (Fig. 5h). However, all young adult mice that were protected by shielding of the head or hindlimbs survived past 16 weeks without transplantation of exogenous


BM, demonstrating that both bone compartments contain sufficient HSCs for survival (Fig. 5h). Secondary haematopoietic reconstitution following an additional round of lethal irradiation with


head or hindlimb shielding also demonstrated the long-term reconstitution potential of HSCs derived from both BM compartments (Extended Data Fig. 10a,b). Whereas shielding of the head of


old mice was also fully compatible with long-term survival after irradiation, the same was not the case for hindlimb shielding, and none of the mice in this group survived beyond 200 days


(Fig. 5h). The calvarial vasculature is also highly resilient against irradiation-induced changes. Consistent with our previous findings40, vascular density (10.2-fold) and diameter


(3.3-fold) were substantially increased in femoral BM at 7 days after irradiation. By contrast, vessels in skull showed minimal changes both with and without BM transplantation (Extended


Data Fig. 10c–e). To explore the potential mechanisms that contribute to BM heterogeneity in ageing skull and femur, we performed single-cell RNA sequencing (scRNA-seq) of KIT+ HSPCs, and


stromal and endothelial cells isolated from either BM compartment in young or geriatric mice (Extended Data Fig. 11a–d). We found significantly higher expression of factors associated with


stress and inflammation such as _Jund_, _Fosb_ and _Dusp1_, as well as myeloid differentiation factors such as _Ngp_, _Ltf_ and _S100a8_, in geriatric femur HSCs compared with skull HSCs


(Extended Data Fig. 12a–f). These differentially expressed genes correlate strongly with gene expression changes reported for IL-1β-treated HSCs or the _Tet2-_knockout model of clonal


haematopoiesis23,48. To further determine whether the observed differences in BM function between ageing skull and femur are owing to HSC- or HSPC-intrinsic or extrinsic, microenvironmental


factors, we performed colony-forming unit (CFU) assays of HSPCs isolated from skull or femur of young or geriatric mice (Extended Data Fig. 13a,b). Although HSPCs derived from both geriatric


skull and femur initially showed higher colony-forming potential than their younger counterparts in primary CFU plating, there were no significant differences in secondary plating or in


HSPC self-renewal or myeloid differentiation potential in both rounds of plating, indicating a role for niche-derived factors. scRNA-seq analysis of endothelial cells from geriatric BM


revealed differences in the expression of inflammatory and myeloid differentiation factors such as _Hotairm1_, _Lgals3_, _Il6_ and _S100a6_, which were significantly upregulated in geriatric


femur compared with skull (Extended Data Fig. 14a–g). Conversely, endothelial cells from skull exhibited higher expression of HSC maintenance factors such as _Ptn_, _Gpr182_ and _Spp1_


(Extended Data Fig. 14h). Finally, to assess the haematopoietic function of skull BM in a physiological context devoid of irradiation injury, we expressed the green-to-red photoconvertible


fluorescent protein, Kikume green–red49 (KikGR) in haematopoietic cells under control of _Vav_-_cre_21 (Fig. 5i). Following in vivo photoconversion of haematopoietic cells in calvarial BM,


flow cytometric analysis of peripheral blood revealed a marked age-dependent increase in systemic haematopoietic contribution from photoconverted skull BM (Fig. 5j). These results show that


systemic haematopoietic output from skull increases relative to other BM compartments during ageing. In summary, our findings reveal an unexpected level of heterogeneity between two major BM


compartments that are well-established model systems in experimental studies. In contrast to nearly all other organ systems and long bone, which become fully vascularized during embryonic


and postnatal development50, blood vessel growth in skull persists throughout adult life. Lifelong angiogenesis in skull involves the continued formation of large-calibre sinusoidal vessels,


which are strongly associated with a range of haematopoietic cell populations and thus represent a central landmark in the skull BM environment. Vascular growth and BM expansion in skull


rely on VEGFA and its receptor VEGFR2, which are key regulators of developmental and regenerative angiogenesis in many organs39,51,52. Consistent with previous studies40, skull BM-resident


HSPCs express more VEGFA than Lin+ differentiated haematopoietic cells, but other relevant sources of VEGFA might exist, as has been shown for chondrocytes and osteoprogenitors in long


bone53. Our experiments also establish that circulation-derived VEGFA can promote calvarial vessel growth and BM expansion, which is consistent with the finding that systemic VEGF can


increase lifespan and delay ageing processes in mice54. Conversely, HSPC behaviour is controlled by the BM endothelium and associated reticular cells, which provide critical niche signals,


such as stem cell factor55 (also known as SCF or KIT ligand) or CXCL1256 (also known as SDF-1). Several studies have shown that haematopoietic cell production is enhanced in response to


various physiological and pathological stimuli, such as pregnancy and ischaemic stroke57. Our present findings indicate that these conditions lead to the expansion of calvarial BM but induce


only limited alterations in load-bearing long bone, where BM growth may not be compatible with mechanical loading and structural integrity. Furthermore, our results establish that calvarial


BM is more resistant to ageing-related degenerative processes, namely adipogenesis, the upregulation of inflammatory cytokines and compromised vascular integrity. This might, in part,


reflect that calvarial BM is biologically ‘younger’ owing to its formation during adult life. Alternatively, unknown cellular or molecular features might differentially influence


haematopoiesis in different BM compartments. Consistent with the latter, recent work has shown that skull BM has a distinct molecular profile in health and neurological disorders58. Notably,


our findings also suggest that BM expansion occurs in human skull, which is consistent with previous observations59,60. Thus, our findings regarding the plasticity and resilience of the BM


in skull may have wider relevance for ageing and disease processes in humans as well as pharmacological treatments that target angiogenic signalling pathways. METHODS ANIMAL MODELS C57BL/6J


mice were used for all experiments involving wild-type mice and pharmacological treatments. Mice at the age of 10–14 weeks, 31–37 weeks, 52–75 weeks and >95 weeks were chosen for young


adult, middle-aged, old and geriatric groups, respectively. Both female and male mice of all age groups were used for initial BM expansion analyses, while only female mice were used for


remaining experiments. _Flk1-_GFP reporter mice17 were used for initial blood vessel characterization. For genetic labelling of haematopoietic cells, _Vav1-cre_ mice21 were interbred with


_ROSA26-mTmG_ reporter mice22 to generate _Vav-mTmG_ mice. For photoconversion of haematopoietic cells, _Vav1-cre_ mice were interbred with _ROSA26_-CAG-_loxP_-stop-_loxP_-KikGR knock-in


mice49 to generate _Vav-KikGR_ mice. For pregnancy experiments, 10-week-old C57BL/6J female mice were paired with 10- to 12-week-old C57BL/6J male mice and the onset of pregnancy was


determined by the presence of a vaginal plug in the morning. Ten-week-old mice received daily intraperitoneal injections of PTH (1–34) (Bachem, 0.1 mg kg−1 for 28 days), PGE2 (Cayman


Chemical, 2 mg kg−1 for 7 days), AMD3100 (Abcam, 5 mg kg−1 for 14 days) before they were euthanized. For DC101 treatment, 10-week-old mice received intraperitoneal injections of DC101


(BioXCell, 40 mg kg−1) every 2 days for 12 weeks. Mice were kept in individually ventilated cages, with constant access to food and water under a 12 h light and 12 h dark cycle regime. Air


flow, temperature (21–22 °C) and humidity (55–60%) were controlled by an air management system. Mice were checked daily and maintained in specific pathogen-free conditions. Sufficient


nesting material and environmental enrichment was provided. All animal experiments were performed according to the institutional guidelines and laws, approved by local animal ethical


committee and were conducted at the Max Planck Institute for Molecular Biomedicine (84-02.04.2016.A160, 81-02.04.2018.A171, 81-02.04.2020.A212, 81-02.04.2020.A416 and 81-02.04.2022.A198),


Universitätsmedizin Berlin (G0220/17), Georg-Speyer-Haus (F123/2017) and the University Medical Center Mainz Institute of Transfusion Medicine (G23-1-067 A1TE) under the indicated


permissions granted by the Landesamt für Natur, Umwelt und Verbraucherschutz (LANUV) of North Rhine-Westphalia, the State Office for Health and Social Affairs Berlin, Regierungspräsidium


Darmstadt and the Landesuntersuchungsamt Rheinland-Pfalz, Germany. HUMAN SUBJECTS, CT ACQUISITION AND DATA ANALYSIS The study was approved by the local ethics committee and the institutional


review board (IRB) of Asan Medical Center, and the requirement for informed consent was waived due to the retrospective nature of the study (IRB number: 2023-0658). The study population


consisted of 36 patients, divided into 4 groups according to age (between 20 and 40 years, over 60 years) and sex (male, female), with 9 individuals in each group. Patients who underwent CT


for evaluation of small cerebral aneurysm from April to May 2023 were eligible. Participants were excluded if they had a previous history of surgery or radiation therapy to the head and


neck, vascular or bone-related medical implants, or a suspicious disease other than small cerebral aneurysm. All human patients underwent CT examinations on the same 128-channel


multidetector CT system (Somatom Definition Edge; Siemens). Imaging variables were as follows: 100 kV; 100 effective mAs; axial scan mode; section thickness, 0.5 mm; display FOV, 20.5 cm;


pitch, 1; gantry rotation time, 0.5 s; pixel matrix, 512 × 512. Images were obtained from the vertex to first cervical spine, without an intravenous injection of contrast media. The CT data


were digitally transferred to a personal computer and processed with ImageJ software (http://rsb.info.nih.gov/ij/). A representative image was selected on a coronal CT image perpendicular to


the outermost convex area on an axial CT image. After whole-bone segmentation of the parietal bone, cortical bone and BM were defined by attenuation densities on CT scan: cortical bone as


over 850 Hounsfield units and BM as less than 850 Hounsfield units and their areas were calculated. SAMPLE PROCESSING AND IMMUNOSTAINING Mice were euthanized by transcardial perfusion of PBS


and 4% paraformaldehyde (PFA), skulls and femur were collected and fixed immediately in ice-cold 4% PFA for 6–8 h under gentle agitation. Bones were decalcified in 0.5 M EDTA for 3 days


(for skulls) or 7 days (for femurs) at 4 °C under gentle shaking agitation, washed 5 times in PBS in 5 min intervals, followed by overnight incubation in cryoprotectant solution (20%


sucrose, 2% polyvinylpyrrolidone) and embedding in bone embedding medium (8% gelatin, 20% sucrose, 2% polyvinylpyrrolidone). Samples were stored overnight at −80 °C. 80-μm-thick cryosections


were prepared for immunofluorescence staining. Bone sections were washed in PBS and permeabilized with 0.3% Triton X-100 in PBS for 10 min at room temperature. Samples were incubated in


blocking solution (5% heat-inactivated donkey serum in 0.3% Triton X-100) for 1 h at room temperature. Primary antibodies (rat monoclonal anti-endomucin (V.7C7) (Santa Cruz, sc-65495, 1:200


dilution), rabbit monoclonal anti-vATPaseB1/B2 (Abcam, 200839, 1:200 dilution), goat polyclonal anti-osteopontin (R&D Systems, AF808, 1:200 dilution), goat polyclonal anti-CD31 (R&D,


AF3628, 1:200 dilution), rabbit polyclonal anti-caveolin-1 (Cell Signaling, 3238, 1:100), goat polyclonal anti-VEGF164 (R&D Systems, AF-493-NA, 1:200 dilution), and rat monoclonal


APC-conjugated anti-CD117 (KIT) (BD Biosciences, 553356, 1:100 dilution) were diluted in PBS with 5% donkey serum and incubated overnight at 4 °C. Next, slides were washed 3–5 times in PBS


in 5 min intervals. Species-specific Alexa Fluor-conjugated secondary antibodies Alexa Fluor 488 (Thermo Fisher Scientific, A21208), Alexa Fluor 594 (Thermo Fisher Scientific, A21209), Alexa


Fluor 647 (Thermo Fisher Scientific, A31573 or A21447) diluted 1:500 in PBS with 5% donkey serum were added and incubated overnight at 4 °C. Slides were washed 3–5 times in PBS in 5 min


intervals. Nuclei were counterstained with DAPI (Sigma-Aldrich, D9542, 1:1,000 dilution). Coverslips were mounted with FluoroMount-G (Southern Biotech, 0100-01). IN VIVO IMMUNOSTAINING AND


EVANS BLUE LEAKAGE ASSAY Rat monoclonal anti-CD31 (BD Biosciences, 553708) was conjugated to Alexa Fluor 647 using the Alexa Fluor 647 Antibody Labeling Kit (Thermo Fisher Scientific,


A20186) according to the manufacturer’s instructions. For blood vessel immunostaining, the conjugated anti-CD31 antibody and rat monoclonal PE-conjugated anti-endomucin (V.7C7) (Santa Cruz,


65495 PE) were diluted 1:10 in 200 μl PBS and injected intravenously into the tail vein. For haematopoietic cell immunostaining, rat monoclonal FITC-conjugated anti-CD45 (eBioscience,


11-0451-82), hamster monoclonal FITC-conjugated anti-CD3e (eBioscience, 16-0031-82), rat monoclonal PE-conjugated anti-CD45R/B220 (BD Biosciences, 553090), rat monoclonal FITC-conjugated


anti-CD11b (BD Biosciences, 553310) were diluted 1:10 in PBS and injected intravenously into the tail vein. Mice were euthanized 1 h after injection with transcardial perfusion with PBS and


4% PFA and bones were collected and fixed immediately in ice-cold 4% PFA for 6–8 h under gentle agitation. The dura mater was carefully removed from the skull with forceps. Bones were


decalcified in 0.5 M EDTA for 1 day (for skulls) or 7 days (for femurs) at 4 °C under gentle shaking agitation, and washed 5 times in PBS in 5 min intervals. Skulls were counterstained with


DAPI (1:500 dilution) for 1 h and trimmed down to the calvarium before mounting with iSpacers (Sunjin Lab, IS011) in PBS. Femurs were cryosectioned, counterstained and mounted as described


above. For the Evans Blue leakage assay, mice were anaesthetized immediately prior to tail vein injection of 200 μl Evans Blue solution (Sigma-Aldrich, E2129, 1% v/w). Mice were euthanized


via transcardial perfusion 5 min after injection as described above. In order to distinguish vascular leakage in the dura mater from the calvarial BM, dura mater tissues were separated from


the calvarial bone before overnight decalcification. Immunostained samples were imaged with a Zeiss LSM980 (Carl Zeiss). Images were analysed, quantified and processed using ZEN Black (Carl


Zeiss, v2.3), ImageJ (NIH, v2.0.0) and IMARIS (Bitplane, v10.0.1). Tilescan overview images of skull BM were superimposed on top of a black background, filling empty corners without image


data. Vessel diameter was measured by selecting the _z_-plane image with the widest vessel diameter from the _z_-stack of the individual vessel. SCANNING ELECTRON MICROSCOPY Skull and femur


from 12-week-old and 73-week-old mice were isolated and submerged in 4% PFA, 0.5% glutaraldehyde, 2 mM MgCl2, 2 mM CaCl2 in 0.1 M cacodylate buffer, pH 7.4, under agitation for 2 h at room


temperature. Samples were fixed further overnight in 2% glutaraldehyde, 2 mM MgCl2, 2 mM CaCl2 in 0.1 M cacodylate buffer, pH 7.4 at 4 °C. Bones were then decalcified over 12 days, changing


solution every other day in 5% EDTA in 0.1 M cacodylate buffer, pH 7.4 under rotation at 4 °C. Subsequently, 150 μm sections were generated with a vibratome (VT 1200, Leica). Sections were


post-fixed in 1% OsmO4, containing 2.5% PFA–glutaraldehyde mixture buffered with 0.1 M phosphate (pH 7.2) for 5 h and then were placed in graded ethanol for critical-point drying using E3000


(Polaron) critical-point dryer. Critical-point-dried bones were placed on a piece of carbon tape and sputter coated with gold in a SC502 Sputter Coater (Polaron). Specimens were imaged on a


Quanta 250 Field Emission Scanning Electron microscope (FEI Quanta 250 FEG, FEI, Hillsboro, OR) installed at the Korea Research Institute of Bioscience and Biotechnology. DURA MATER


WHOLE-MOUNT IMMUNOSTAINING Mice were euthanized by transcardial perfusion of PBS and 4% PFA, skulls were collected and fixed immediately in ice-cold 4% PFA for 6–8 h under gentle agitation.


Skulls were decalcified in 0.5 M EDTA for 24 h at 4 °C under gentle shaking agitation, washed 5 times in PBS in 5 min intervals, trimmed down to the calvarium, and incubated in blocking


solution (5% heat-inactivated donkey serum in 0.3% Triton X-100) for 1 h at room temperature. Goat polyclonal anti-CD31 (R&D, AF3628, 1:100 dilution) was diluted in PBS with 5% donkey


serum and incubated overnight at 4 °C with gentle agitation. Samples were washed 3–5 times in PBS in 10 min intervals. Alexa Fluor 647 (Thermo Fisher Scientific, A21447) diluted 1:500 in PBS


with 5% donkey serum was added and incubated overnight at 4 °C with gentle agitation. Samples were washed 3–5 times in PBS in 10 min intervals and mounted with iSpacers (Sunjin Lab, IS011)


in PBS. TMCAO Sixteen-week-old female C57BL6/J mice were used throughout the experiments. Mice were anaesthetized by intraperitoneal injection of a mixture of 10 mg kg−1 xylazine (cp-pharma)


and 90 mg kg−1 ketamine hydrochloride (cp-pharma). Throughout the whole procedure and during recovery, body temperature was maintained at 37 °C via a heating pad. After ligation of the left


proximal common carotid artery and external carotid artery, a 7.0-nylon monofilament (Doccol) with a 0.23-mm coated tip was introduced into the distal internal carotid artery via an


incision in the ligated common carotid artery. The monofilament was advanced distal to the carotid bifurcation to occlude the middle cerebral artery. Arter topical application of the local


anaesthetic lidocaine hydrochloride (Xylocain Spray 2%, Aspen) the neck wound was closed temporarily for a 45 min ischaemic period. At reperfusion, the monofilament was withdrawn from the


carotid artery and the wound was stitched with 4-0 non-resorbable sutures (Ethibond Excel, Ethicon) and the single s.c. injection of Penicillin G 20 000 U (Benzylpenicillin-Natrium,


InfectoPharm) was given. The mouse was returned to its cage to recover under observation. CHRONIC MYELOID LEUKAEMIA Six-week-old female C57BL/6 mice were purchased from Charles River


Laboratories and were used as donors and recipients in all transplants. The transplantation experiments were performed as previously described61. In brief, to induce CML-like


myeloproliferative neoplasia, donor BM cells from donor mice pre-treated with 5-fluorouracil (200 mg kg−1 intravenously; 4 days prior to collection) were pre-stimulated overnight in medium


containing SCF (50 ng ml−1), IL-6 (10 ng ml−1) and IL-3 (6 ng ml−1) and transduced on two consecutive days with murine stem cell virus (MSCV)-IRES-GFP-BCR-ABL1 to induce CML or MSCV-IRES-GFP


control virus. Subsequently, transduced cells were intravenously transplanted (2.5 × 105 cells per mouse) into sublethally irradiated (900 cGy) recipient mice. Mice were euthanized 14 days


after transplantation. LINEAGE DEPLETION AND TRANSPLANTATION In order to transplant lineage-negative BM cells, young (10- to 14-week-old) or old (52- to 75-week-old) donor mice were


euthanized and skulls were collected. The calvarium was first chopped with scissors in FACS buffer (PBS with 2% fetal calf serum), then crushed with a mortar and pestle. Cell suspension was


filtered through a 40-μm mesh filter (Falcon, 352340), resuspended in RBC lysing buffer (Sigma-Aldrich, R7757) for red blood cell lysis and washed with FACS buffer. Cells were resuspended in


FACS buffer and incubated with a biotinylated anti-haematopoietic lineage antibody cocktail (Miltenyi-Biotec, 130-092-613, 1:10 dilution), followed by washing with FACS buffer and


incubation with R-PE-conjugated streptavidin secondary antibody (Invitrogen, S866, 1:50 dilution). DAPI (1:1,000 dilution) was added to resuspended cells to distinguish live and dead cells


and were FACS-sorted for live lineage-negative cells on a FACSAria Fusion (BD Biosciences). Sorted cells were intravenously transplanted (5×105 cells/mouse) into lethally irradiated (12 Gy,


Best Theratronics, Gammacell 40 Exactor) recipient mice (12-week-old). Mice were euthanized 14 days after transplantation. FACS ANALYSIS OF BM AND PERIPHERAL BLOOD Mice from each age group


were euthanized and skull and femur were collected. Skulls were chopped with scissors in FACS buffer before crushed with mortar and pestle; femurs were crushed without chopping. BM stromal


samples were dissociated with Collagenase I (Gibco, 17100-017, 2 mg ml−1) and Collagenase IV (Gibco, 17104-019, 2 mg ml−1) in PBS for 20 min at 37 °C with intermittent shaking. Cell


suspensions were strained through a 40-μm mesh filter, resuspended in RBC lysing buffer (when applicable) and washed with FACS buffer. Cells were resuspended and incubated with the following


primary antibodies in FACS buffer: biotinylated rat monoclonal anti-haematopoietic lineage antibody cocktail (Miltenyi-Biotec, 130-092-613, 1:50 dilution), APC-conjugated rat monoclonal


anti-CD117 (BD Biosciences, 553356, 1:100 dilution), FITC-conjugated rat monoclonal anti-CD117 (Biolegend, 105806, 1:100 dilution), FITC-conjugated rat monoclonal anti-Ly6A/E (Sca1)


(eBioscience, 11-5981-85, 1:100 dilution), PerCP-Cy5.5-conjugated rat monoclonal anti-Ly-6A/E (Invitrogen, 45-5981, 1:100), APC-Cy7-conjugated hamster monoclonal anti-CD48 (BD Biosciences,


561242, 1:100 dilution), PE-conjugated rat monoclonal anti-CD150 (SLAM) (Biolegend, 115904, 1:100 dilution), Alexa Fluor 647-conjugated rat monoclonal anti-CD150 (Biolegend, 115918, 1:100


dilution), PE-Cy7-conjugated rat monoclonal anti-CD45 (eBioscience, 25-0451-82, 1:100 dilution), BV421-conjugated rat monoclonal anti-TER-119 (Biolegend, 116234, 1:100 dilution),


FITC-conjugated rat monoclonal anti-CD71 (Biolegend, 113806, 1:100 dilution), Alexa Fluor 647-conjugated rat monoclonal anti-CD31 (BD Biosciences, 553708, conjugation described above, 1:200


dilution), PE-conjugated rat monoclonal anti-endomucin (Santa Cruz, 65495 PE, 1:100 dilution), PE-Cy7-conjugated rat monoclonal anti-CD16/32 (eBioscience, 25-0161, 1:100 dilution), eFluor


450-conjugated rat monoclonal anti-CD34 (eBioscience, 48-0341, 1:100 dilution), PE-conjugated rat monoclonal anti-CD127 (eBioscience, 12-1271, 1:100 dilution), BV711-conjugated rat


monoclonal anti-CD41 (BD Biosciences, 740712, 1:100 dilution), PE-conjugated rat monoclonal anti-CD105 (eBioscience, 12-1051-82, 1:100 dilution), APC-conjugated hamster monoclonal anti-CD3e


(eBioscience, 17-0031, 1:100 dilution), PE-conjugated rat monoclonal anti-CD45R/B220 (BD Biosciences, 553090, 1:100 dilution), FITC-conjugated rat monoclonal anti-CD11b (BD Biosciences,


553310, 1:100 dilution). Cells were washed, resuspended in FACS buffer with Alexa Fluor 405-conjugated (Invitrogen, S32351, 1:100 dilution) or APC-Cy7-conjugated (BD Biosciences, 554063,


1:100 dilution) streptavidin secondary antibody, washed again before analysis with a FACSymphony A5 Cell Analyzer (BD Biosciences). Peripheral blood was collected from the submandibular vein


with lancets (Medipoint) into EDTA-coated tubes. Blood was resuspended in RBC lysing buffer and washed with FACS buffer. Cells were resuspended and incubated with the following primary


antibodies in FACS buffer: biotinylated rat monoclonal anti-haematopoietic lineage antibody cocktail (Miltenyi-Biotgec, 130-092-613, 1:50 dilution), APC-conjugated rat monoclonal anti-CD117


(BD Biosciences, 553356, 1:100 dilution), FITC-conjugated rat monoclonal anti-Ly6A/E (Sca1) (eBioscience, 11-5981-85, 1:100 dilution), Pacific Blue-conjugated mouse monoclonal anti-CD45.2


(Biolegend, 109820, 1:100 dilution), APC-conjugated hamster monoclonal anti-CD3e (eBioscience, 17-0031, 1:100 dilution), PE-conjugated rat monoclonal anti-CD45R/B220 (BD Biosciences, 553090,


1:100 dilution), FITC-conjugated rat monoclonal anti-CD11b (BD Biosciences, 553310, 1:100 dilution). Cells were washed and resuspended in FACS buffer before analysis with a FACSymphony A5


Cell Analyzer (BD Biosciences). RNA EXTRACTION AND QUANTITATIVE PCR FACS-sorted cells from 10-week-old mouse skulls were lysed and RNA was extracted using a Monarch Total RNA Miniprep Kit


(New England BioLabs, T2010S). Extracted RNA concentration was measured with a NanoDrop 8000 Spectrophotometer (Thermo Fisher Scientific) and cDNA was generated with a LuncaScript RT


SuperMix Kit (New England BioLabs, E3010L). Quantitative PCR with reverse transcription was performed with a BioRad CFX96 real-time PCR system using FAM-conjugated Taqman probes for _Vegfa_


(Mm00437306_m1) or using PowerUp SYBR Green Master Mix (Applied Biosystems, A25742) with primers designed using Pimer-BLAST or adopted from previously published studies: _Vegfa__120_


(5′-AACGATGAAGCCCTGGAGTG-3′; 5′-TGAGAGGTCTGGTTCCCGA-3′); _Vegfa__164_ (5′-AACGATGAAGCCCTGGAGTG; 5′-GACAAACAAATGCTTTCTCCG-3′); _Vegfa__188_ (5′-AACGATGAAGCCCTGGAGTG-3′;


5′-AACAAGGCTCACAGTGAACG-3′). Gene expression levels were normalized to the endogenous VIC-conjugated _Gapdh_ probe (44326317E) as control. ELISA Mice from each age group were euthanized,


bones were collected. Skulls were chopped before being crushed with a mortar and pestle in ice-cold RIPA lysis buffer; femurs were crushed without chopping. Supernatants of centrifuged


lysates were further concentrated using an Ultra-0.5 Centrifugal Filter Unit with a 3 kDa cutoff (Millipore, UFC500396), resulting concentrations were measured using a Pierce BCA Protein


Assay Kit (Thermo Fisher Scientific, 23225), and the concentrations of VEGFA in tissue extracts were measured using a Mouse VEGFA Quantikine ELISA Kit (R&D Systems, MMV00-1). HYPOXIA


ANALYSIS Hypoxic cells were detected with the hypoxia probe pimonidazole (Pimo, Hypoxyprobe) according to the manufacturer’s instructions. Mice were intraperitoneally injected with 60 mg 


kg−1 1 h before analysis. VEGFA PLASMID CONSTRUCTION AND OVEREXPRESSION To generate the pLIVE-VEGFA165-HA-MP-Asp8x bone-homing protein containing VEGF165 fused to a HA tag, metalloprotease


and 8x Asp peptide sequences, a cDNA fragment encoding amino acids 1–191 of human _VEGFA_ was amplified via PCR using the following oligonucleotide primers: VEGFA-AscI-Fwd:


5′-ATGAACTTTCTGCTGTCT-3′ and VEGFA-XhoI-Rev: 5′-CCGCCTCGGCTTGTCACATCTGCA-3′ and annealed with the NEBuilder Assembly Cloning Kit. Ten-week-old mice were used for hydrodynamic tail vein


injection. Mice were injected with 0.5 μg g−1 (plasmid/body weight) pLIVE-Vegfa plasmid suspended in TransIT-EE hydrodynamic delivery solution (Mirus, MIR5340). The appropriate amount of


plasmid was suspended in an injection volume of 10% of the body weight and injected into each individual mouse via the tail vein in 5–7 s as previously reported62. ADIPOCYTE ANALYSIS To


stain for neutral lipids, the entire calvarium or femur cryosections were incubated in BODIPY 493/503 (Invitrogen, D3922; 1:1,000 dilution) for 1 h at room temperature with gentle agitation


(only calvarium). Samples were washed with PBS 3–5 times at 5 min intervals before mounting. ANALYSIS OF INFLAMMATORY CYTOKINES Mice from each age group were euthanized and bones were


collected. Skulls were chopped before being crushed with a mortar and pestle in ice-cold RIPA lysis buffer; femurs were crushed without chopping. Supernatants of centrifuged lysates were


further concentrated using an Ultra-0.5 Centrifugal Filter Unit with a 3 kDa cutoff (Millipore, UFC500396), resulting concentrations were measured using a Pierce BCA Protein Assay Kit


(Thermo Fisher Scientific, 23225), and concentrations of inflammatory cyotokines were measured with LEGENDplex Mouse Inflammation Panel (13-plex) with V-bottom plates (Biolegend, 740446).


Analysis on a FACSymphony (BD Biosciences) and quantification were performed according to the manufacturer’s protocol. Data analysis was performed using software provided by Biolegend.


Manual gating was used to define beads A and B, and automatic gating strategy was used to gate individual cytokines in the APC–PE plot. IRRADIATION WITH PARTIAL SHIELDING Mice were


anaesthetized with ketamine (100 mg kg−1) and xylazine (10 mg kg−1) prior to irradiation. For partial shielding, the entire head or both legs of a mouse were inserted into the opening of the


cylindrical 1-inch-thick lead shield (JRT Associates, PTI-50-P) and exposed to lethal irradiation (12 Gy). The mouse was returned to its cage to recover under observation. SAMPLE


PREPARATION FOR SCRNA-SEQ Mice from each age group were euthanized and skull and femur were collected. Skulls were chopped with scissors in FACS buffer before crushed with mortar and pestle;


femurs were crushed without chopping. BM stromal samples were dissociated with Collagenase I (Gibco, 17100-017, 2 mg ml−1) and Collagenase IV (Gibco, 17104-019, 2 mg ml−1) in PBS for 20 min


at 37 °C with intermittent shaking. Cell suspensions were strained through a 40-μm mesh filter and washed with FACS buffer. Cells were resuspended and incubated with biotinylated rat


monoclonal anti-haematopoietic lineage antibody cocktail (Miltenyi-Biotec, 130-092-613, 1:50 dilution). Cells were washed, resuspended in FACS buffer with mouse monoclonal anti-Biotin


MicroBeads (Miltenyi-Biotec, 130-105-637, 1:50 dilution) and incubated before being loaded into a magnetic-associated cell sorting (MACS) column (Miltenyi-Biotec, 130-042-201) for lineage


depletion. Lin− cells were further incubated with rat monoclonal anti-CD45 MicroBeads (Miltenyi-Biotec, 130-052-301, 1:50 dilution), rat monoclonal anti-CD117 MicroBeads (Miltenyi-Biotec,


130-091-224, 1:50 dilution), biotinylated rat monoclonal anti-CD71 (Biolegend, 113803, 1:100). Cells were washed, resuspended in FACS buffer with mouse monoclonal anti-Biotin MicroBeads


(Miltenyi-Biotec, 130-105-637, 1:50 dilution) and incubated before being loaded into a magnetic-associated cell sorting (MACS) column (Miltenyi-Biotec, 130-042-201) for further


haematopoietic depletion. Single-cell suspensions were processed with BD Rhapsody and scRNA-seq libraries were evaluated and quantified by Agilent Bioanalyzer using High Sensitivity DNA Kit


(Agilent Technologies, 5067-4626) and Qubit (Thermo Fisher Scientific, Q32851). Individual libraries were diluted to 4 nM and pooled for sequencing. Pooled libraries were sequenced by using


High Output Kit (Illumina, TG-160-2002) with a NextSeq500 sequencer (Illumina). SCRNA-SEQ Preprocessing: STAR version 2.7.10a (PMID: 23104886) was used to generate a reference genome index


for GRCm39, with Gencode annotations vM29, subset to lncRNA and protein-coding genes. FASTQ reads were mapped against the reference genome index using STAR with the settings “--soloType


CB_UMI_Complex --soloCellFilter None --outSAMtype BAM SortedByCoordinate --soloFeatures GeneFull_Ex50pAS --soloCBmatchWLtype 1MM --soloUMIlen 8 --soloCBwhitelist BD_CLS1.txt BD_CLS2.txt


BD_CLS3.txt --runRNGseed 1 --soloMultiMappers EM --readFilesCommand zcat --outSAMattributes NH HI AS nM NM MD jM jI MC ch CB UB GX GN sS CR CY UR UY”. Libraries using standard BD Rhapsody


beads were mapped using the adapter parameters “--soloAdapterSequence NNNNNNNNNACTGGCCTGCGANNNNNNNNNGGTAGCGGTGACA --soloCBposition 2_0_2_8 2_21_2_29 3_1_3_9 --soloUMIposition 3_10_3_17”,


libraries with BD Rhapsody enhanced beads with --soloAdapterSequence NNNNNNNNNGTGANNNNNNNNNGACA --soloCBposition 2_0_2_8 2_13_2_21 3_1_3_9 --soloUMIposition 3_10_3_17. Raw counts were


imported as AnnData63 objects. We removed low complexity barcodes with the knee plot method, and further filtered out cells with a mitochondrial mRNA content, as well as unusually high total


and gene counts using manually determined cutoffs for each sample. Doublets were scored with scrublet64. Finally, each sample’s gene expression matrix was normalized using scran65 (1.22.1)


with Leiden clustering66 input at resolution 0.5. G2M and S phase scores were assigned to each cell using gene lists from ref. 67 and the scanpy68 (1.9.6) sc.tl.score_genes_cell_cycle


function. Embedding, clustering and annotation: different combinations of samples and cell populations (all, ECs, HSCs), were used as input for 2D embedding and clustering: the corresponding


expression matrix was subset to the 2,000 most highly variable genes (sc.pp.highly_variable_genes, flavour “seurat”). The top 50 principal components were calculated, and batch-corrected


using Harmony69 (0.0.9). The principal components served as basis for _k_-nearest neighbour calculation (sc.pp.neighbors, n_neighbors=30), which were used as input for UMAP70 layout


(sc.tl.umap, min_dist=0.3). Cell populations were clustered using scanpy.tl.leiden, and a suitable resolution was chosen for a first-pass annotation. Here, contaminating cell populations,


including multiplet clusters, were removed, and clustering was repeated. Cluster marker genes were calculated using a pseudobulk approach, comparing aggregate counts with 2 pseudoreplicates


for each cluster to all remaining cells (pyDeSEQ2 0.4.8). Finally, expression of select marker genes was plotted using Matplotlib71 (3.8.4) imshow, and clusters were annotated accordingly.


Differential expression analysis: Differentially expressed genes were calculated using a pseudobulk approach, comparing aggregate counts with two pseudoreplicates for each condition


(pyDeSEQ2 0.4.8). SKULL BM PHOTOCONVERSION _Vav1-KikGR_ mice were anaesthetized with ketamine (100 mg kg−1) and xylazine (10 mg kg−1). A skin flap was generated to expose the calvarium, as


previously described72. Each exposed area of the calvarium was then exposed to UV light from a Zeiss Axio Imager (Zeiss Microscopy) for 60 s, confirmed for photoconversion from green to red


fluorescence, before exposing another area. The skin flap was sutured back together and peripheral blood was analysed by flow cytometry, as described above, to check for the presence of


photoconverted cells, which were non-existent in the peripheral blood immediately after photoconversion. One week after photoconversion, peripheral blood was drawn, stained for Alexa


Fluor-conjugated rat monoclonal anti-CD45, and was analysed by flow cytometry for CD45+ photoconverted haematopoietic cells derived from the skull BM. STATISTICAL ANALYSIS No statistical


methods were used to predetermine sample size. The experiments were randomized and investigators were blinded to allocation during experiments and outcome analyses. All values are presented


as mean ± s.d. Statistical significance was determined by the two-tailed unpaired Student’s _t_-test between two groups or the Tukey multiple comparison test (one-way ANOVA) for


multiple-group comparison. Statistical analyses were performed using GraphPad Prism 9.0 (GraphPad Software). Statistical significance was set at _P_ < 0.05. REPORTING SUMMARY Further


information on research design is available in the Nature Portfolio Reporting Summary linked to this article. DATA AVAILABILITY The scRNA-seq data generated in this study have been deposited


in the Gene Expression Omnibus under accession number GSE275179. The mouse reference genome GRCm39 with GENCODE M26 annotation (https://www.gencodegenes.org/mouse/release_M26.html) was used


for mapping the reads in this study. All individual mouse lines used in this study are commercially available at The Jackson Laboratory. Plasmid constructs are available, upon request,


through the corresponding author. Source data are provided with this paper. CHANGE HISTORY * _ 27 NOVEMBER 2024 In the Methods, two instances of “cGy” have been updated to instead read as


“Gy”. _ REFERENCES * Hoggatt, J., Kfoury, Y. & Scadden, D. T. Hematopoietic stem cell niche in health and disease. _Annu. Rev. Pathol._ 11, 555–581 (2016). CAS  PubMed  Google Scholar  *


Verovskaya, E. V., Dellorusso, P. V. & Passegue, E. Losing sense of self and surroundings: hematopoietic stem cell aging and leukemic transformation. _Trends Mol. Med._ 25, 494–515


(2019). PubMed  PubMed Central  Google Scholar  * Comazzetto, S., Shen, B. & Morrison, S. J. Niches that regulate stem cells and hematopoiesis in adult bone marrow. _Dev. Cell_ 56,


1848–1860 (2021). CAS  PubMed  PubMed Central  Google Scholar  * van den Bos, T., Speijer, D., Bank, R. A., Bromme, D. & Everts, V. Differences in matrix composition between calvaria and


long bone in mice suggest differences in biomechanical properties and resorption: Special emphasis on collagen. _Bone_ 43, 459–468 (2008). PubMed  Google Scholar  * Sivaraj, K. K. &


Adams, R. H. Blood vessel formation and function in bone. _Development_ 143, 2706–2715 (2016). CAS  PubMed  Google Scholar  * Mills, W. A., Coburn, M. A. & Eyo, U. B. The emergence of


the calvarial hematopoietic niche in health and disease. _Immunol. Rev._ 311, 26–38 (2022). CAS  PubMed  PubMed Central  Google Scholar  * Lo Celso, C. et al. Live-animal tracking of


individual haematopoietic stem/progenitor cells in their niche. _Nature_ 457, 92–96 (2009). ADS  CAS  PubMed  Google Scholar  * Christodoulou, C. et al. Live-animal imaging of native


haematopoietic stem and progenitor cells. _Nature_ 578, 278–283 (2020). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Lassailly, F., Foster, K., Lopez-Onieva, L., Currie, E. &


Bonnet, D. Multimodal imaging reveals structural and functional heterogeneity in different bone marrow compartments: functional implications on hematopoietic stem cells. _Blood_ 122,


1730–1740 (2013). CAS  PubMed  Google Scholar  * Herisson, F. et al. Direct vascular channels connect skull bone marrow and the brain surface enabling myeloid cell migration. _Nat.


Neurosci._ 21, 1209–1217 (2018). CAS  PubMed  PubMed Central  Google Scholar  * Cugurra, A. et al. Skull and vertebral bone marrow are myeloid cell reservoirs for the meninges and CNS


parenchyma. _Science_ 373, eabf7844 (2021). CAS  PubMed  PubMed Central  Google Scholar  * Nakamura-Ishizu, A., Ito, K. & Suda, T. Hematopoietic stem cell metabolism during development


and aging. _Dev. Cell_ 54, 239–255 (2020). CAS  PubMed  PubMed Central  Google Scholar  * Potente, M., Gerhardt, H. & Carmeliet, P. Basic and therapeutic aspects of angiogenesis. _Cell_


146, 873–887 (2011). CAS  PubMed  Google Scholar  * Chung, A. S. & Ferrara, N. Developmental and pathological angiogenesis. _Annu. Rev. Cell Dev. Biol._ 27, 563–584 (2011). CAS  PubMed 


Google Scholar  * Ramasamy, S. K. et al. Regulation of hematopoiesis and osteogenesis by blood vessel–derived signals. _Annu. Rev. Cell Dev. Biol._ 32, 649–675 (2016). CAS  PubMed  Google


Scholar  * Ramalingam, P., Butler, J. M. & Poulos, M. G. Vascular regulation of hematopoietic stem cell homeostasis, regeneration, and aging. _Current Stem Cell Rep._ 7, 194–203 (2021).


Google Scholar  * Xu, Y. et al. Neuropilin-2 mediates VEGF-C–induced lymphatic sprouting together with VEGFR3. _J. Cell Biol._ 188, 115–130 (2010). CAS  PubMed  PubMed Central  Google


Scholar  * Kusumbe, A. P. et al. Age-dependent modulation of vascular niches for haematopoietic stem cells. _Nature_ 532, 380–384 (2016). ADS  CAS  PubMed  PubMed Central  Google Scholar  *


Langen, U. H. et al. Cell–matrix signals specify bone endothelial cells during developmental osteogenesis. _Nat. Cell Biol._ 19, 189–201 (2017). CAS  PubMed  PubMed Central  Google Scholar 


* Yu, J. Direct evidence for the role of caveolin-1 and caveolae in mechanotransduction and remodeling of blood vessels. _J. Clin. Invest._ 116, 1284–1291 (2006). CAS  PubMed  PubMed Central


  Google Scholar  * de Boer, J. et al. Transgenic mice with hematopoietic and lymphoid specific expression of Cre. _Eur. J. Immunol._ 33, 314–325 (2003). PubMed  Google Scholar  * Muzumdar,


M. D., Tasic, B., Miyamichi, K., Li, L. & Luo, L. A global double-fluorescent Cre reporter mouse. _Genesis_ 45, 593–605 (2007). CAS  PubMed  Google Scholar  * Mitchell, C. A. et al.


Stromal niche inflammation mediated by IL-1 signalling is a targetable driver of haematopoietic ageing. _Nat. Cell Biol._ 25, 30–41 (2023). CAS  PubMed  PubMed Central  Google Scholar  *


Olson, O. C., Kang, Y. A. & Passegue, E. Normal hematopoiesis is a balancing act of self-renewal and regeneration. _Cold Spring Harb. Perspect. Med._ 10, a035519 (2020). CAS  PubMed 


PubMed Central  Google Scholar  * Nakada, D. et al. Oestrogen increases haematopoietic stem-cell self-renewal in females and during pregnancy. _Nature_ 505, 555–558 (2014). ADS  CAS  PubMed


  PubMed Central  Google Scholar  * Oguro, H. et al. 27-Hydroxycholesterol induces hematopoietic stem cell mobilization and extramedullary hematopoiesis during pregnancy. _J. Clin. Invest._


127, 3392–3401 (2017). PubMed  PubMed Central  Google Scholar  * Rodrigues, J. et al. Estrogen enforces the integrity of blood vessels in the bone during pregnancy and menopause. _Nat.


Cardiovasc. Res._ 1, 918–932 (2022). PubMed  PubMed Central  Google Scholar  * Courties, G. et al. Ischemic stroke activates hematopoietic bone marrow stem cells. _Circ. Res._ 116, 407–417


(2015). CAS  PubMed  Google Scholar  * Mendez-Ferrer, S. et al. Bone marrow niches in haematological malignancies. _Nat. Rev. Cancer_ 20, 285–298 (2020). CAS  PubMed  PubMed Central  Google


Scholar  * Reid, I. R. & Billington, E. O. Drug therapy for osteoporosis in older adults. _Lancet_ 399, 1080–1092 (2022). CAS  PubMed  Google Scholar  * Calvi, L. M. et al. Osteoblastic


cells regulate the haematopoietic stem cell niche. _Nature_ 425, 841–846 (2003). ADS  CAS  PubMed  Google Scholar  * Akunuru, S. & Geiger, H. Aging, clonality, and rejuvenation of


hematopoietic stem cells. _Trends Mol. Med._ 22, 701–712 (2016). CAS  PubMed  PubMed Central  Google Scholar  * Morrison, S. J., Wandycz, A. M., Akashi, K., Globerson, A. & Weissman, I.


L. The aging of hematopoietic stem cells. _Nat. Med._ 2, 1011–1016 (1996). CAS  PubMed  Google Scholar  * Colom Diaz, P. A., Mistry, J. J. & Trowbridge, J. J. Hematopoietic stem cell


aging and leukemia transformation. _Blood_ 142, 533–542 (2023). CAS  PubMed  PubMed Central  Google Scholar  * North, T. E. et al. Prostaglandin E2 regulates vertebrate haematopoietic stem


cell homeostasis. _Nature_ 447, 1007–1011 (2007). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Porter, R. L. et al. Prostaglandin E2 increases hematopoietic stem cell survival and


accelerates hematopoietic recovery after radiation injury. _Stem Cells_ 31, 372–383 (2013). CAS  PubMed  Google Scholar  * Hoggatt, J. et al. Rapid mobilization reveals a highly engraftable


hematopoietic stem cell. _Cell_ 172, 191–204.e110 (2018). CAS  PubMed  Google Scholar  * Simons, M., Gordon, E. & Claesson-Welsh, L. Mechanisms and regulation of endothelial VEGF


receptor signalling. _Nat. Rev. Mol. Cell Biol._ 17, 611–625 (2016). CAS  PubMed  Google Scholar  * Apte, R. S., Chen, D. S. & Ferrara, N. VEGF in signaling and disease: beyond discovery


and development. _Cell_ 176, 1248–1264 (2019). CAS  PubMed  PubMed Central  Google Scholar  * Chen, Q. et al. Apelin+ endothelial niche cells control hematopoiesis and mediate vascular


regeneration after myeloablative injury. _Cell Stem Cell_ 25, 768–783.e766 (2019). CAS  PubMed  PubMed Central  Google Scholar  * Schipani, E., Wu, C., Rankin, E. B. & Giaccia, A. J.


Regulation of bone marrow angiogenesis by osteoblasts during bone development and homeostasis. _Front. Endocrinol._ 4, 85 (2013). Google Scholar  * Wong, B. W., Marsch, E., Treps, L., Baes,


M. & Carmeliet, P. Endothelial cell metabolism in health and disease: impact of hypoxia. _EMBO J._ 36, 2187–2203 (2017). CAS  PubMed  PubMed Central  Google Scholar  * Suda, T. &


Liu, D. Hydrodynamic gene delivery: its principles and applications. _Mol. Ther._ 15, 2063–2069 (2007). CAS  PubMed  Google Scholar  * Ho, Y.-H. & Méndez-Ferrer, S. Microenvironmental


contributions to hematopoietic stem cell aging. _Haematologica_ 105, 38–46 (2020). CAS  PubMed  PubMed Central  Google Scholar  * Ambrosi, T. H. et al. Adipocyte accumulation in the bone


marrow during obesity and aging impairs stem cell-based hematopoietic and bone regeneration. _Cell Stem Cell_ 20, 771–784.e776 (2017). CAS  PubMed  PubMed Central  Google Scholar  * Ho, Y.


H. et al. Remodeling of bone marrow hematopoietic stem cell niches promotes myeloid cell expansion during premature or physiological aging. _Cell Stem Cell_ 25, 407–418.e406 (2019). CAS 


PubMed  PubMed Central  Google Scholar  * Brioschi, S. et al. Heterogeneity of meningeal B cells reveals a lymphopoietic niche at the CNS borders. _Science_ 373, eabf9277 (2021). CAS  PubMed


  PubMed Central  Google Scholar  * McClatchy, J. et al. Clonal hematopoiesis related TET2 loss-of-function impedes IL1β-mediated epigenetic reprogramming in hematopoietic stem and


progenitor cells. _Nat. Commun._ 14, 8102 (2023). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Tomura, M. et al. Tracking and quantification of dendritic cell migration and antigen


trafficking between the skin and lymph nodes. _Sci. Rep._ 4, 6030 (2014). CAS  PubMed  PubMed Central  Google Scholar  * Maes, C. Role and regulation of vascularization processes in


endochondral bones. _Calcif. Tissue Int._ 92, 307–323 (2013). CAS  PubMed  Google Scholar  * Carmeliet, P. & Jain, R. K. Molecular mechanisms and clinical applications of angiogenesis.


_Nature_ 473, 298–307 (2011). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Hooper, A. T. et al. Engraftment and reconstitution of hematopoiesis is dependent on VEGFR2-mediated


regeneration of sinusoidal endothelial cells. _Cell Stem Cell_ 4, 263–274 (2009). CAS  PubMed  PubMed Central  Google Scholar  * Hu, K. & Olsen, B. R. Osteoblast-derived VEGF regulates


osteoblast differentiation and bone formation during bone repair. _J. Clin. Invest._ 126, 509–526 (2016). PubMed  PubMed Central  Google Scholar  * Grunewald, M. et al. Counteracting


age-related VEGF signaling insufficiency promotes healthy aging and extends life span. _Science_ 373, eabc8479 (2021). CAS  PubMed  Google Scholar  * Ding, L., Saunders, T. L., Enikolopov,


G. & Morrison, S. J. Endothelial and perivascular cells maintain haematopoietic stem cells. _Nature_ 481, 457–462 (2012). ADS  CAS  PubMed  PubMed Central  Google Scholar  * Greenbaum,


A. et al. CXCL12 in early mesenchymal progenitors is required for haematopoietic stem-cell maintenance. _Nature_ 495, 227–230 (2013). ADS  CAS  PubMed  PubMed Central  Google Scholar  *


Heidt, T. et al. Chronic variable stress activates hematopoietic stem cells. _Nat. Med._ 20, 754–758 (2014). CAS  PubMed  PubMed Central  Google Scholar  * Kolabas, Z. I. et al. Distinct


molecular profiles of skull bone marrow in health and neurological disorders. _Cell_ 186, 3706–3725.e3729 (2023). CAS  PubMed  PubMed Central  Google Scholar  * Sabanciogullari, V. et al.


Diploe thickness and cranial dimensions in males and females in mid-Anatolian population: an MRI study. _Forensic Sci. Int._ 219, 289.e281–287 (2012). Google Scholar  * Lillie, E. M., Urban,


J. E., Lynch, S. K., Weaver, A. A. & Stitzel, J. D. Evaluation of skull cortical thickness changes with age and sex from computed tomography scans. _J. Bone Miner. Res._ 31, 299–307


(2016). PubMed  Google Scholar  * Krause, D. S. et al. Differential regulation of myeloid leukemias by the bone marrow microenvironment. _Nat. Med._ 19, 1513–1517 (2013). CAS  PubMed  PubMed


Central  Google Scholar  * Liu, F., Song, Y. & Liu, D. Hydrodynamics-based transfection in animals by systemic administration of plasmid DNA. _Gene Ther._ 6, 1258–1266 (1999). CAS 


PubMed  Google Scholar  * Virshup, I., Rybakov, S., Theis, F. J., Angerer, P. & Wolf, F. A. anndata: annotated data. Preprint at _bioRxiv_ https://doi.org/10.1101/2021.12.16.473007v1


(2021). * Wolock, S. L., Lopez, R. & Klein, A. M. Scrublet: computational identification of cell doublets in single-cell transcriptomic data. _Cell Syst._ 8, 281–291.e9 (2019). CAS 


PubMed  PubMed Central  Google Scholar  * Lun, A. T. L., Bach, K. & Marioni, J. C. Pooling across cells to normalize single-cell RNA sequencing data with many zero counts. _Genome Biol._


17, 75 (2016). PubMed  Google Scholar  * Traag, V. A., Waltman, L. & van Eck, N. J. From Louvain to Leiden: guaranteeing well-connected communities. _Sci. Rep._ 9, 5233 (2019). ADS  CAS


  PubMed  PubMed Central  Google Scholar  * Tirosh, I. et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. _Science_ 352, 189–196 (2016). ADS  CAS 


PubMed  PubMed Central  Google Scholar  * Wolf, F. A., Angerer, P. & Theis, F. J. SCANPY: large-scale single-cell gene expression data analysis. _Genome Biol._ 19, 15 (2018). PubMed 


PubMed Central  Google Scholar  * Korsunsky, I. et al. Fast, sensitive and accurate integration of single-cell data with Harmony. _Nat. Methods_ 16, 1289–1296 (2019). CAS  PubMed  PubMed


Central  Google Scholar  * McInnes, L., Healy, J. & Melville, J. UMAP: uniform manifold approximation and projection for dimension reduction. Preprint at


https://doi.org/10.48550/arXiv.1802.03426 (2018). * Hunter, J. D. Matplotlib: a 2D graphics environment. _Comput. Sci. Eng._ 9, 90–95 (2007). Google Scholar  * Koh, B. I. et al. VEGFR2


signaling drives meningeal vascular regeneration upon head injury. _Nat. Commun._ 11, 3866 (2020). ADS  PubMed  PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS The


authors thank M. Stehling (MPI Münster), J. Kirsch (FCCF, DRFZ) and C. Conche (HI-TRON Mainz, DKFZ) for FACS technical assistance; F. Berkenfeld for supporting animal colony maintenance; and


K. Müller for sequencing technical assistance. V.M. is a part of the Cells in Motion-International Max Planck Research School (CiM-IMPRS) and University of Münster, Germany. R.H.A. is


supported by the Max Planck Society, the European Research Council (AdG 101139772, PROTECT), the DFG (CRC 1366, project no. 394046768), and the Leducq Foundation. FUNDING Open access funding


provided by Max Planck Society. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Department of Tissue Morphogenesis, Max Planck Institute for Molecular Biomedicine, Münster, Germany Bong Ihn


Koh, Vishal Mohanakrishnan, Hongryeol Park, Susanne Adams, M. Gabriele Bixel & Ralf H. Adams * Sequencing Core Facility, Max Planck Institute for Molecular Biomedicine, Münster, Germany


Hyun-Woo Jeong * Bioinformatics Service Unit, Max Planck Institute for Molecular Biomedicine, Münster, Germany Kai Kruse * Department of Radiology and Research Institute of Radiology, Asan


Medical Center, University of Ulsan College of Medicine, Seoul, Republic of Korea Young Jun Choi * Department of Neurosurgery, Charité-Universitätsmedizin Berlin, corporate member of Freie


Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany Melina Nieminen-Kelhä & Peter Vajkoczy * Institute of Transfusion Medicine, Transfusion Center, University Medicine


Mainz, Mainz, Germany Rahul Kumar & Daniela S. Krause * Georg-Speyer-Haus Institute for Tumor Biology and Experimental Medicine and Goethe University Frankfurt, Frankfurt, Germany


Raquel S. Pereira * Center for Vascular Research, Institute for Basic Science, Daejeon, Republic of Korea Hyuek Jong Lee Authors * Bong Ihn Koh View author publications You can also search


for this author inPubMed Google Scholar * Vishal Mohanakrishnan View author publications You can also search for this author inPubMed Google Scholar * Hyun-Woo Jeong View author publications


You can also search for this author inPubMed Google Scholar * Hongryeol Park View author publications You can also search for this author inPubMed Google Scholar * Kai Kruse View author


publications You can also search for this author inPubMed Google Scholar * Young Jun Choi View author publications You can also search for this author inPubMed Google Scholar * Melina


Nieminen-Kelhä View author publications You can also search for this author inPubMed Google Scholar * Rahul Kumar View author publications You can also search for this author inPubMed Google


Scholar * Raquel S. Pereira View author publications You can also search for this author inPubMed Google Scholar * Susanne Adams View author publications You can also search for this author


inPubMed Google Scholar * Hyuek Jong Lee View author publications You can also search for this author inPubMed Google Scholar * M. Gabriele Bixel View author publications You can also


search for this author inPubMed Google Scholar * Peter Vajkoczy View author publications You can also search for this author inPubMed Google Scholar * Daniela S. Krause View author


publications You can also search for this author inPubMed Google Scholar * Ralf H. Adams View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS


B.I.K. and R.H.A. designed the study and interpreted results. B.I.K. designed, conducted and analysed the majority of the experiments, including all pharmacological treatments, skull sample


processing, immunostainings, transplantations, FACS sorting and analyses, RNA and protein analyses, and photoconversion surgeries. V.M. processed femur samples and provided technical


assistance. H.-W.J. and H.P. provided technical assistance. K.K. performed scRNA-seq analyses. Y.J.C. conducted and analysed human head CT scans. M.N.-K. and P.V. designed and performed


transient middle cerebral artery occlusion surgeries. R.K., R.S.P. and D.S.K. generated the CML mouse model and contributed to project design. S.A. generated pLIVE plasmids. H.J.L. conducted


scanning electron microscopy sample processing and imaging. M.G.B. contributed to project design and data interpretation. R.H.A. supervised the project. CORRESPONDING AUTHORS Correspondence


to Bong Ihn Koh or Ralf H. Adams. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. PEER REVIEW PEER REVIEW INFORMATION _Nature_ thanks the anonymous


reviewer(s) for their contribution to the peer review of this work. Peer review reports are available. ADDITIONAL INFORMATION PUBLISHER’S NOTE Springer Nature remains neutral with regard to


jurisdictional claims in published maps and institutional affiliations. EXTENDED DATA FIGURES AND TABLES EXTENDED DATA FIG. 1 AGE-DEPENDENT EXPANSION OF BM IN SKULL. A, B. DAPI staining and


quantification of mouse skull coronal cryosections showing expansion of BM cellular content during adulthood and aging (n = 4 mice/group from three independent experiments). Blue arrowheads


indicate location of magnified inset. C, D. IF staining and quantification of skull coronal cryosections showing increased vATPase+ activated osteoclasts attached to Osteopontin+ bone


surfaces in old versus young specimen (n = 4 mice/group from two independent experiments). Scale bars, 1 mm. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple


comparison test (one-way ANOVA) and two-tailed unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 2 AGE-RELATED EXPANSION OF DIPLOIC SPACE IN HUMAN SKULL. A. Method for


quantification of whole bone, cortex and bone marrow areas in human patient CT scans. B, C. Representative coronal CT images (B) and quantification (C) of young (21–40 years old) and old


(61–69 years old) female and male human skulls showing enlargement of diploic space with aging (n = 9 patients/group). Scale bars, 2 cm. Vertical bars indicate mean ± SD. _P_ values were


calculated using two-tailed unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 3 REGIONAL DIFFERENCES AND EXPANSION OF SKULL BM. A, C. In vivo IF staining (A) and quantification of


BM vessels (C) in _Flk1-GFP_ reporter mice showing distinct blood vessel architecture in frontal (F), parietal (P) and interparietal (i-P) skull. Arrowheads indicate Flk1low/− CD31+


arterioles. n = 6 mice/group for vascular area, n = 20 randomly selected vessels from all samples/group for vessel diameter from two independent experiments. Boxed areas in overview image


(left) indicate location of magnified images. Scale bar, 1 mm. B, D. Imaging (B) and quantification (D) of Evans Blue vascular leakage in different skull bone parts (n = 4 mice/group from


three independent experiments). Boxed areas in overview image (left) indicate location of magnified images. Scale bar, 1 mm. E. IF staining of skull BM blood vessels showing increase in


vascular network complexity with aging. Dotted lines demarcate dorsal and ventral skull boundaries. Representative images from four independent experiments. Scale bars, 100 μm. F, G. IF


staining and quantification of caveolin-1-positive (Cav1+) blood vessels showing an increase with aging (n = 8 mice/group from two independent experiments). Dotted lines demarcate dorsal and


ventral boundaries of the skull. Scale bars, 100 μm. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple comparison test (one-way ANOVA) and two-tailed


unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 4 ASSOCIATION OF VESSELS AND HEMATOPOIETIC CELLS IN SKULL BM. A, B. In vivo IF staining (A) and quantification (B) of CD45+


hematopoietic cells in skull BM during aging (n = 4 mice/group from three independent experiments). Green arrowheads indicate dense clusters of CD45+ cells. Scale bars, 500 μm. C. Genetic


labeling of hematopoietic cells in _Vav-Cre Rosa26-mTmG_ mice (Vav-mTmG) shows association (arrowheads) of GFP-expressing cells with large CD31+ vessels. Representative images from three


independent experiments. Scale bars, 150 μm. D. Confocal images showing association of CD3ɛ+ T lymphocytes (green), B220+ B lymphocytes (red) and CD11b+ cells (green) with large-caliber


CD31+ BM vessel in aged skull. Representative images from three independent experiments. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple comparison test


(one-way ANOVA). Scale bars, 500 μm. Source Data EXTENDED DATA FIG. 5 AGE-RELATED CHANGES IN THE BM VASCULATURE. A, B. In vivo IF staining (A) and quantification (B) of femoral BM vessels


showing vascular deterioration during aging by comparing young (Y), middle-aged (M), old (O), and geriatric (G) specimen (n = 5 mice/group for vascular density and Endomucin intensity, n = 


10 randomly selected vessels from all samples/group from three independent experiments). Scale bars, 500 μm. C. Electron micrographs showing the surface of BM endothelial cells in young vs.


old skull or femur, as indicated. Note regular pattern of fenestrations (black arrowheads) in young and old skull but disorganized pattern with larger gaps (yellow arrowheads) in old femur.


Representative images from three independent experiments. Scale bars, 1 μm. D, E. IF staining (D) and quantification (E) of dural blood vessels revealing substantial decrease in vascular


density with aging. n = 3 (young), n = 6 (old) mice/group from three independent experiments. Scale bars, 250 μm. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey


multiple comparison test (one-way ANOVA) and two-tailed unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 6 PATHOPHYSIOLOGICAL CHANGES IN THE FEMORAL BM VASCULATURE. A, B. In vivo


IF staining (A) and quantification (B) of femoral BM blood vessels in pregnant (17 dpc) and post-partum (2 dpp) female mice. n = 3 mice/group from three independent experiments. C, D. In


vivo IF staining (C) and quantification (D) of femoral BM blood vessels in mice 7 days after transient mid-cerebral artery occlusion (tMCAO). n = 6 (sham), n = 8 (tMCAO) from three


independent experiments. E, F. In vivo IF staining (E) and quantification (F) of femoral BM blood vessels in mice with chronic myeloid leukemia. n = 4 mice/group from three independent


experiments. G, H. In vivo IF staining (G) and quantification (H) of femoral BM blood vessels in mice with 28-day sustained parathyroid hormone (PTH) treatment. n = 4 mice/group from three


independent experiments. Scale bars, 1 mm. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple comparison test (one-way ANOVA) and two-tailed unpaired Student’s


_t_-test. Source Data EXTENDED DATA FIG. 7 EFFECT OF TREATMENTS AFFECTING HSPCS. A, B. In vivo IF staining (A) and quantification (B) of femoral BM blood vessels in recipient mice receiving


hematopoietic lineage-negative cells isolated from young or old donors (young skull (YS), old skull (OS), young femur (YF), old femur (OF)). n = 4 mice/group from three independent


experiments. Scale bars, 1 mm. C, D. Diagram showing experimental scheme for PGE2 treatment for HSPC expansion (C) and quantification (D) of HSPCs in skull BM and spleen weight with PGE2


treatment (n = 3 mice/group from three independent experiments). E, F. Diagram showing experimental scheme for AMD3100 treatment for HSPC mobilization (E) and quantification (F) of HSPCs in


peripheral blood and spleen weight with AMD3100 treatment (n = 5 mice/group from three independent experiments). G-J. In vivo IF staining (G, I) and quantification (H, J) of femoral BM blood


vessels in mice treated with PGE2 (G, H) or AMD3100 (I, J) for HSPC expansion or mobilization, respectively. n = 6 (PGE2 Control), n = 5 (PGE2), n = 4 (all other groups) mice/group from


three independent experiments. Scale bars, 1 mm. Vertical bars indicate mean ± SD. _P_ values were calculated using two-tailed unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 8


VEGF-A RESPONSE TO PATHOPHYSIOLOGICAL CONDITIONS. A. VEGF-A protein concentrations in total BM lysates from skull or femur isolated from mice in various pathophysiological conditions,


including pregnancy, stroke, chronic myeloid leukemia and sustained PTH treatment. n = 4 (Pregnancy, Stroke and PTH) and n = 5 (CML) mice/group from two independent experiments. Vertical


bars indicate mean ± SD. _P_ values were calculated using two-tailed unpaired Student’s _t_-test. Source Data EXTENDED DATA FIG. 9 VEGF EXPRESSION AND HYPOXIA IN BM COMPARTMENTS. A. qRT-PCR


analyses of _Vegfa_ mRNA in FACS-sorted LSK (LIN- Sca1+ cKit+), cKit+ (LIN- Sca1- cKit+), Lin neg (LIN- Sca1- cKit-) and Lin pos (LIN+) populations isolated from femur BM of young or old


mice (n = 5 mice pooled/sample from three independent experiments). B. qRT-PCR analyses of _Vegfa_ mRNA isoforms in populations shown in (A) isolated from skull or femur BM of young or old


mice (n = 5 mice pooled/sample from three independent experiments). C, D. IF staining (C) and quantification (D) of IV-injected Hypoxyprobe in young, old and geriatric skull and femoral BM,


as indicated. Note increase of Hypoxyprobe signal in skull but not in femur (arrowheads, n = 5 mice/group from two independent experiments). Scale bars, 1 mm. Vertical bars indicate mean ± 


SD. **_P_ < 0.01, ***_P_ < 0.001 versus Young LSK (a) or Young (d), #_P_ < 0.05 versus Old by Tukey multiple comparison test (one-way ANOVA). E, H. Diagram showing experimental


scheme for _Vegfa_ overexpression or anti-VEGFR2 blocking antibody DC101 treatment. F-J. In vivo IF staining (F, I) and quantification (G, J) of femoral BM blood vessels in mice after


hydrodynamic injection of bone-homing _Vegfa_ construct (n = 4 mice/group from three independent experiments) (F, G) or treatment with anti-VEGFR2 blocking antibody (DC101) for 12 weeks (n =


 5 mice/group from two independent experiments) (I, J). Scale bars, 1 mm. Vertical bars indicate mean ± SD. _P_ values were calculated using two-tailed unpaired Student’s _t_-test. Source


Data EXTENDED DATA FIG. 10 LONG-TERM RECONSTITUTION POTENTIAL OF SKULL HSCS. A. Diagram showing experimental scheme for serial hematopoietic reconstitution following serial partial


irradiation with shielding. B. Monthly FACS analysis of T, B, and myeloid cells in peripheral blood isolated from animals with head versus leg shielding during primary and secondary


hematopoietic reconstitution (n = 5 (Head shield), n = 4 (Leg shield) mice/group from two independent experiments). Vertical bars indicate mean ± SD. C-E. Fluorescence imaging (C) and


quantification (D, E) of skull or femoral BM blood vessels in Cdh5-mTnG (mTomato, red) reporter mice at 7 days after lethal irradiation with (+) or without (-) bone marrow transplantation


(BMT), as indicated. Note minimal vascular alterations in the skull BM compared to substantial vascular changes in femoral BM (n = 6 mice/group for vascular area/density, n = 10 randomly


selected vessels from all samples/group for vessel diameter from three independent experiments). Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple comparison


test (one-way ANOVA). Source Data EXTENDED DATA FIG. 11 SINGLE-CELL RNA-SEQUENCING ANALYSIS OF HSPCS AND STROMAL CELLS IN SKULL AND FEMUR BONE MARROW DURING AGING. A. UMAP plot showing


color-coded merged cell clusters of FACS-sorted cKit+ HSPCs and hematopoietic lineage-depleted bone marrow cells. Hematopoietic stem cell (HSC), multi-potent progenitor (MPP), myeloid


progenitor (MyPro), lymphoid progenitor (LyPro), megakaryocyte-erythrocyte progenitor (MEP), pre-B-cell (Pre-B), sinusoidal endothelial cell (sEC), arterial endothelial cell (aEC), bone


marrow mesenchymal stromal cell (bmMSC), metaphyseal mesenchymal stromal cell (mpMSC), smooth muscle cell (SMC), osteoblast (Ob). B. Heat map showing top cell marker genes of each cell


population shown in (A). C. Group-selected color-coded cell clusters from young or old, skull or femur as indicated. D. Frequency plots of color-coded HSPCs (top) and endothelial cells


(bottom) from young skull (YS), old skull (OS), young femur (YF), old femur (OF). EXTENDED DATA FIG. 12 SINGLE-CELL RNA-SEQUENCING ANALYSIS OF HSCS IN SKULL AND FEMUR BONE MARROW DURING


AGING. A. UMAP plot showing color-coded merged cell clusters of HSCs identified in Extended Data Fig. 11. B. Heat map showing top cell marker genes of each cell population shown in (A). C.


Frequency plots of color-coded HSCs from young skull (YS), old skull (OS), young femur (YF), old femur (OF). D. Differentially expressed genes (DEG) analysis comparing HSCs from geriatric


skull versus femur. E. DEG analysis comparing each HSC subcluster from geriatric skull versus femur. F. Violin plots showing the expression of selected genes with log normalized values for


inflammatory and myeloid determination factors in each HSC subcluster identified in (A). Skull HSC1 (n = 6216 cells), HSC2 (n = 3713 cells), HSC3 (n = 768 cells) and Femur HSC1 (n = 6726


cells), HSC2 (n = 3006 cells), HSC3 (n = 719 cells) isolated from n = 5 mice from two independent experiments. The box of each boxplot starts in the first quartile and ends in the third,


with the line inside representing the median. Blue dots, _p_-adjusted value < 0.01 and Log2 fold change <0.5; orange dots, _p_-adjusted value < 0.01 and Log2 fold change > 0.5 by


two-tailed unpaired Student’s _t_-test. EXTENDED DATA FIG. 13 HSPCS ARE DEPENDENT ON MICROENVIRONMENT DURING AGING. A, B. Brightfield image (A) and quantification (B) of primary and


secondary colony forming unit (CFU) assays of 500 FACS-sorted LSK cells isolated from young or geriatric skull or femur bone marrow (n = 4 mice/group from two independent experiments). Scale


bar, 5 mm. Vertical bars indicate mean ± SD. _P_ values were calculated using Tukey multiple comparison test (one-way ANOVA). Source Data EXTENDED DATA FIG. 14 SINGLE-CELL RNA-SEQUENCING


ANALYSIS OF ECS IN SKULL AND FEMUR BONE MARROW DURING AGING. A. UMAP plot showing color-coded merged cell clusters of ECs identified in Extended Data Fig. 11. B. Heat map showing top cell


marker genes of each cell population shown in (A). C. Frequency plots of color-coded ECs from young skull (YS), old skull (OS), young femur (YF), old femur (OF). D, E. Violin plots showing


the expression of selected genes with log normalized values for vascular signaling factors in ECs (D) and each EC subcluster (E) from geriatric skull versus femur. Skull ECs (n = 2414


cells), Arterial (n = 1070 cells), Sinusoidal 1 (n = 329 cells), Sinusoidal 2 (n = 833 cells), Venous (n = 182 cells) and Femur ECs (n = 594 cells), Arterial (n = 430 cells), Sinusoidal 1 (n


 = 19 cells), Sinusoidal 2 (n = 58 cells), Venous (n = 87 cells) isolated from n = 5 mice from two independent experiments. The box of each boxplot starts in the first quartile and ends in


the third, with the line inside representing the median. F. Differentially expressed genes (DEG) analysis comparing ECs from geriatric skull versus femur. G, H. Heat map showing the


expression of selected genes with log normalized values for inflammatory factors (G) and HSC-regulating factors (H) identified in (F). Blue dots, _p_-adjusted value < 0.01 and Log2 fold


change <0.5; orange dots, _p_-adjusted value < 0.01 and Log2 fold change > 0.5 by two-tailed unpaired Student’s _t_-test. SUPPLEMENTARY INFORMATION REPORTING SUMMARY SUPPLEMENTARY


DATA 1 Flow cytometric analyses of skull or femur BM in ageing, pregnancy, stroke, CML, PTH treatment, young versus old transplantation, PGE2 treatment, AMD3100 treatment, pLIVE-Vegfa


overexpression and DC101 treatment. SUPPLEMENTARY DATA 2 Gating strategies for all FACS isolation and flow cytometric analyses performed in the study. PEER REVIEW FILE SOURCE DATA SOURCE


DATA FIGS. 1–5 AND SOURCE DATA EXTENDED DATA FIGS. 1–10 AND 13. RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License,


which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link


to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless


indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or


exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. Reprints


and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Koh, B.I., Mohanakrishnan, V., Jeong, HW. _et al._ Adult skull bone marrow is an expanding and resilient haematopoietic reservoir.


_Nature_ 636, 172–181 (2024). https://doi.org/10.1038/s41586-024-08163-9 Download citation * Received: 10 October 2023 * Accepted: 07 October 2024 * Published: 13 November 2024 * Issue Date:


05 December 2024 * DOI: https://doi.org/10.1038/s41586-024-08163-9 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a


shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative