Research | Open | Published:
Bone marker gene expression in calvarial bones: different bone microenvironments
Journal of Biological Research-Thessalonikivolume 24, Article number: 9 (2017)
In calvarial mice, mesenchymal stem cells (MSCs) differentiate into osteoprogenitor cells and then differentiate into osteoblasts that differentiate into osteocytes, which become embedded within the bone matrix. In this case, the cells participating in bone formation include MSCs, osteoprogenitor cells, osteoblasts and osteocytes. The calvariae of C57BL/KaLwRijHsD mice consist of the following five bones: two frontal bones, two parietal bones and one interparietal bone. This study aimed to analyse some bone marker genes and bone related genes to determine whether these calvarial bones have different bone microenvironments.
C57BL/KaLwRijHsD calvariae were carefully excised from five male mice that were 4–6 weeks of age. Frontal, parietal, and interparietal bones were dissected to determine the bone microenvironment in calvariae. Haematoxylin and eosin staining was used to determine the morphology of different calvarial bones under microscopy. TaqMan was used to analyse the relative expression of Runx2, OC, OSX, RANK, RANKL, OPG, N-cadherin, E-cadherin, FGF2 and FGFR1 genes in different parts of the calvariae.
Histological analysis demonstrated different bone marrow (BM) areas between the different parts of the calvariae. The data show that parietal bones have the smallest BM area compared to frontal and interparietal bones. TaqMan data show a significant increase in the expression level of Runx2, OC, OSX, RANKL, OPG, FGF2 and FGFR1 genes in the parietal bones compared with the frontal and interparietal bones of calvariae.
This study provides evidence that different calvarial bones, frontal, parietal and interparietal, contain different bone microenvironments.
C57BL/KalwRijHsD mice develop a high frequency of monoclonal proliferative B cell disorders. Most C57BL/KalwRijHsD mice have a monoclonal gammopathy of undetermined significance (MGUS) that is similar to humans. The 5TMM series of myeloma models originate from spontaneously developed multiple myeloma (MM) in ageing C57BL/KalwRijHsD mice and have many of the features of the human disease [1,2,3]. Several murine models of MM exist, but 5T2, 5T33 and 5TGM1 are the best characterized and have been used in most recent studies [4,5,6,7,8,9,10]. Al-Amer (2015) recently used a confocal laser scanning microscope LSM 510, multiphoton microscopy, to trace the homing of DiD+ myeloma cells and determine the distribution of these cells in calvarial bone BM. The data showed myeloma cells arrive to all parts of calvariae after 3 days of injection . However, in another study, the same author showed myeloma tumours do not grow in all parts of the calvariae . These findings raised the question of whether parts of calvariae have different bone microenvironments.
Osteoblasts are mononucleate cells derived from putative MSCs that are responsible for bone formation and regulation of osteoclast differentiation. Bone morphogenic proteins (BMPs), Wnt/β-catenin and Notch signalling play a major role in osteoblastogenesis [13,14,15]. In addition, Runt-related transcription factor 2 (RUNX2) has been shown to be an essential transcription factor for osteoblastogenesis. It was found that Runx2−/− mice failed to develop calvarial bones, and the mice died at birth. Runx2 also regulates another transcription factor that is necessary for osteoblastogenesis, osterix (OSX) . Mature osteoblasts are identified by their cuboidal structure and location on the endosteal surface. Type I collagen is the major product of the bone-forming osteoblast in addition to osteocalcin (OC) .
Bone remodelling is a complex process in which mature bone tissue is removed by a process called bone resorption, and new bone tissue is formed by a process called bone formation. In this process, cellular activity and molecular mechanisms are closely coordinated to ensure the bone resorption-formation sequence is performed at a mutual location, sustaining the bone mass. Osteoclastic bone resorption is controlled by receptor activator of nuclear factor kappa-B ligand (RANKL). The binding of RANKL to RANK starts the osteoclastogenesis process and stimulates osteoclast activity and bone resorption . Osteoprotegerin (OPG), a soluble tumour necrosis factor (TNF) receptor superfamily, is produced by many cell types, such as bone-marrow stromal cells and osteoblasts. OPG binds to RANKL and protects the bone from excessive resorption by preventing the RANK/RANKL interaction [19,20,21]. It was determined that overexpression of OPG reduces osteoclastogenesis in mice; however, a lack of OPG has been shown to accelerate osteoclastogenesis, and OPG Ko mice (OPG−/− mice) develop severe osteoporosis .
Neural (N)-cadherin is an adhesion molecule that has been implicated in the localization of haematopoietic stem cells (HSCs) to ‘niches’ containing osteoblasts on endosteal bone surfaces. HSCs home into the endosteal surfaces and attach to N-cadherin+ osteoblasts. This attachment keeps the HSCs in a quiescent state in the HSC niche [23, 24]. Epithelial (E-) cadherin plays an essential role in Ca2+-dependent cell-to-cell adhesion. A striking feature of the type I (classic) cadherins is a highly conserved His-Ala-Val (HAV) motif in the first extracellular domain that mediates homophilic interactions between cadherin molecules expressed on adjacent cells . On the other side of the cell membrane, the cadherin intracellular domain forms a protein complex with catenin molecules, which anchor the cadherins to the cytoskeleton and provide a unique signal transduction pathway for cadherin interactions .
Fibroblast growth factor receptors (FGFRs) play major roles in skeletogenesis, and activating mutations of the human FGFR1, FGFR2 and FGFR3 genes cause premature fusion of the skull bones (craniosynostosis) . This genetic evidence establishes a role for FGFR1 in skeletal development and suggests that the FGFR1, FGFR2 and FGFR3 signalling pathways may have similar or redundant functions. In addition to premature fusion of cranial sutures, several of the classic craniosynostosis syndromes have associated phenotypes that affect long bone development in the appendicular skeleton. This genetic evidence establishes a role for FGFR1 in skeletal development and suggests that the FGFR1, FGFR2 and FGFR3 signalling pathways may have similar or redundant functions. In addition to premature fusion of cranial sutures, several of the classic craniosynostosis syndromes have associated phenotypes that affect long bone development in the appendicular skeleton . FGFRs 1–3 are expressed in the developing and mature skeleton in patterns suggestive of both unique and redundant function .
In the present study, I evaluated the relative expression levels of osteoblast markers (Runx2, OC and OSX) and the regulation of bone resorption molecules (RANK, RANKL and OPG), adhesion molecules (N-cadherin and E-cadherin) and skeletogenesis markers (FGF2 and FGFR1) in different parts of the calvariae. The osteoblast markers (Runx2, OC, and OSX) and the regulation of bone resorption markers (RANKL and OPG) and skeletogenesis markers (FGF2 and FGFR1) are significantly highly expressed in parietal bone compared with frontal and interparietal bones. However, there were no differences in the relative expression levels of the adhesion molecules (N-cadherin and E-cadherin) between the frontal, parietal and interparietal bones. This study provides evidence that different calvarial bones contain different bone microenvironments.
Male C57BL/KalwRiJHsD mice, aged 4–6 weeks, were purchased from Harlan, Netherlands and from the University of Leeds, UK. Mice were housed in the University of Sheffield biological services laboratory. All animals were provided with food and water ad libitum as well as light, and all procedures were performed under a personal license (40/10118).
C57BL/KaLwRijHsd calvariae were carefully excised from five male mice that were 4–6 weeks old. Mice were handled carefully. The frontal, parietal, and interparietal bones were carefully dissected, and sections of each calvarial part were cut for three levels. Each level cut was 200 μm after the previous level. The total sample size was 45, including 15 samples from frontal bones, 15 from parietal bones and 15 from interparietal bones. Each sample was used in triplicate during TaqMan analysis.
Haematoxylin and eosin staining
Haematoxylin and eosin stain (H&E) is a standard staining method in histological studies. Haemalum is formed from aluminium ions and oxidized haematoxylin stains nuclei with a blue colour, while eosin stains the cytoplasm and extracellular matrix with varying degrees of pink [30, 31]. Bone sections were dipped twice in xylene for 5 min to remove wax and rehydrated twice in absolute alcohol, 95% alcohol and 70% alcohol for 5 min before they were briefly washed in water. Gill’s II haematoxylin (VWR, UK) was used for 120 s to stain nuclei. Sections were washed in water. Then, 1% aqueous eosin (VWR, UK) with 1% calcium carbonate (Sigma, UK) was used for 5 min to stain the cytoplasm. Sections were dehydrated in 70% alcohol for 10 s, 95% alcohol for 10 s and twice in absolute alcohol for 30 s each. Xylene was used for 1 min to remove alcohol and for 3 min for clearance. Slides were mounted and covered using coverslips (VWR, UK). In the endo-cortical and trabecular surfaces of the calvariae, the nuclei were stained blue and the cytoplasm and extracellular matrices were stained pink.
High-resolution micro-computed tomography (μCT)
A high resolution μCT scanner (model 1172; Skyscan, Belgium) is a non-destructive technique that provides three-dimensional bone microstructure. A μCT scanner was used to scan calvarial bones. Image were captured every 0.35° through a 360° rotation with a pixel size of 4.1. Scanned images were reconstructed using Skyscan NRecon software (version 220.127.116.11, Skyscan, Belgium).
Leitz DMRB microscopy
Osteomeasure software version 4.10 (Osteomerics Incorporated, Atlanta, USA) was used with a Leitz DMRB microscope and a drawing tablet (CalComp Drawing Broad III) to quantify the BM area of calvariae. This microscopy was used to visualize non-fluorescent cells and sections.
Total RNA was isolated from primary osteoblasts at different stages of development using TRIzol. RNAs were transcribed to cDNA using SuperScript™ II reverse transcriptase (Invitrogen, UK). TaqMan probes for Runx2, OC, OSX, RANK, RANKL, OPG, N-cadherin, E-cadherin, FGF2 and FGFR1 were used to quantify different genes with TaqMan Universal Master Mix II (Applied Biosystems, UK) and normalized with the level of housekeeping gene, β-actin. An ABI PRISM >900 sequence detector was used to quantify the expression of Runx2, OC, OSX, RANK, RANKL, OPG, N-cadherin, E-cadherin, FGF2 and FGFR1 in different calvarial bones.
Statistical significance was determined using Tukey’s test using GraphPad Prism 6 (GraphPad Software Inc., California, USA). Data were considered significant for p ≤ 0.05 (* p < 0.05, ** p < 0.01, and *** p < 0.001). All data were presented as the mean ± standard deviation (±SD).
Huge BM area in interparietal bone compared with frontal and parietal bones
The BM area in C57BL/KaLwRijHsD calvariae was determined. Cortical bone was identified as an outside solid bone with a regular structure, while trabecular bone was identified as a bone inside the BM with an irregular structure (Fig. 1). The endo-cortical bone surface is a surface of the cortical bone inside the BM. Histological analysis demonstrated statistically significant increases in the BM area in interparietal bones compared to frontal and parietal bones (Fig. 2).
Relative level expression of Runx2, OC and OSX in calvarial bones
The expression levels of osteoblast markers, Runx2, OC and OSX, were analysed using TaqMan. Expression of these genes was normalized to the expression of a housekeeping gene, β-actin. Figure 3a shows a significant increase in the relative expression of Runx2 gene in the parietal bones compared with the frontal and interparietal bones. Figure 3b shows a significant increase in the relative expression of the OC gene in the parietal bones compared with the frontal and interparietal bones. Figure 3c shows a significant increase in the relative expression of the OSX gene in the parietal bones compared with the frontal and interparietal bones.
Relative expression levels of RANK, RANKL and OPG in calvarial bones
The expression levels of bone resorption regulation molecules, RANK, RANKL and OPG, were analysed using TaqMan. The expression of these genes was normalized to the expression of a housekeeping gene, β-actin. Figure 4a shows that RANK was only expressed in interparietal bone, but it was not expressed in frontal or parietal bone. Figure 4b shows that RANKL was expressed in parietal and interparietal bones, but it was not expressed in frontal bone. Data show a significant increase in the relative expression of the RANKL gene in the parietal bones compared with the frontal and interparietal bones. Figure 4c shows a significant increase in the relative expression of the OPG gene in the parietal bones compared with the frontal and interparietal bones.
Relative expression levels of N-cadherin and E-cadherin in calvarial bones
The expression levels of adhesion molecules, N-cadherin and E-cadherin, were analysed using TaqMan. The expression of these genes was normalized to the expression of a housekeeping gene, β-actin. Both genes are expressed in all parts of calvariae. Figure 5a shows a non significant increase in the relative expression of the N-cadherin gene in the parietal bones compared with the frontal and interparietal bones. Figure 5b shows an increase in the relative expression of the E-cadherin gene in the frontal bones compared with the parietal and interparietal bones, but it was not significant, although p value was close to significant (p = 0.06).
Relative expression levels of FGF2 and FGFR1 in calvarial bones
The expression levels of skeletogenesis markers, FGF2 and FGFR1, were analysed using TaqMan. The expression of these genes was normalized to the expression of a housekeeping gene, β-actin. Figure 6a shows a significant increase in the relative expression of the FGF2 gene in the frontal and parietal bones compared with the interparietal bones, but there was no significant different in the relative expression between the frontal and parietal bones. Figure 6b shows a significant increase in the relative expression of the FGFR1 gene in the parietal bones compared with the frontal and interparietal bones.
In this study, different calvariae—frontal, parietal and interparietal—were analysed. The aim of this study was to test the hypothesis that calvarial bones contain different bone microenvironments. In previous studies, C57BL/KaLwRijHsd mice were inoculated with myeloma cells via their tail vein. Multi-photon microscopy demonstrated a significant increase in the number of myeloma cells that homed to the interparietal bones 3 and 21 days post inoculation compared to frontal and parietal bones . In a separate study, C57BL/KaLwRijHsd mice were also inoculated with myeloma cells via their tail veins. Using the Lightools image analysis system with ImageJ software demonstrated that a large tumour was formed in interparietal bones and not in frontal or parietal bones after 3 weeks .
To date, no publication has demonstrated whether different calvarial bones—frontal, parietal or interparietal—contain different bone microenvironments in C57BL/KaLwRijHsD mice. Parfitt et al. (1987) summarized a unified system of terminology for bone histomorphometry. The bone histomorphometry community decided to measure the BM volume as the marrow area (Ma.Ar) . According to this terminology, the BM areas of the frontal, parietal and interparietal regions of calvariae were analysed in this study. It has been reported that calvariae of C57BL/KaLwRijHsD mice consist of the following five bones: two frontal bones, two parietal bones and one interparietal bone . In this study, the BM area of different parts of C57BL/KaLwRijHsd calvariae was analysed under a microscope. Interestingly, data showed that there was a significant increase in the BM area of interparietal bones compared to frontal and parietal bones. Together, these studies suggest the possibility that there are different bone microenvironments in the calvariae.
Healthy bone is regulated and coordinated by a bone remodelling cycle consisting of bone formation and bone resorption throughout the skeleton. One of the most important consequences of cancer metastasis in the bone is the induced imbalance in bone remodelling between osteoblasts and osteoclasts, which affects bone resorption and bone formation. This imbalance between osteoclasts and osteoblasts is responsible for an increase in bone resorption and a reduction in bone formation, resulting in osteolytic bone lesions [33,34,35,36]. C57BL/KalwRij mice are a pre-clinical model used to study MM bone disease. Most C57BL/KalwRij mice have a monoclonal gammopathy of undetermined significance (MGUS) similar to humans [1,2,3]. To support our understanding of the bone microenvironment in the pre-clinical model to study MM bone disease, the relative expression levels of osteoblast markers (Runx2, OC and OSX) and the regulation of bone resorption molecules (RANK, RANKL and OPG), adhesion molecules (N-cadherin and E-cadherin) and skeletogenesis markers (FGF2 and FGFR1) in different calvariae were analysed.
It has been reported that Runx2, OC and OSX play an important role in osteoblastogenesis, and a deficiency in one of them plays a critical role in murine bone development [16, 17]. On the other hand, it has been reported that the RANKL/RANK/OPG system plays an important role in bone resorption. RANKL/RANK signalling promotes osteoclast differentiation and activation, leading to bone resorption. OPG binds to RANKL and protects bone from excessive resorption [18,19,20,21,22]. For osteoblast markers (Runx2, OC and OSX), this study demonstrated that there was a significant increase in the relative expression levels of Runx2, OC and OSX in parietal bones compared with frontal and interparietal bones. For bone resorption regulation molecules (RANK, RANKL and OPG), this study demonstrated that there was a significant increase in the relative expression levels of RANKL and OPG in parietal bones compared to frontal and interparietal bones. However, the expression of RANK was significantly higher in interparietal bones compared with frontal and parietal bones. These findings show that the parietal bone has a high expression level of osteoblast markers (Runx2, OC and OSX) and RANKL/OPG, blocking bone resorption. This could explain the finding that there is a small area of bone marrow in parietal bone that shows more bone formation and blocks bone resorption. On the other hand, interparietal bone has a low expression level of osteoblast markers (Runx2, OC and OSX) in addition to a high expression level of RANK. This could explain the finding that there is a large bone marrow area in interparietal bone where there is more bone resorption and less bone formation. This study also showed that there is a significant increase in skeletogenesis markers (FGF2 and FGFR1) in parietal bones. On the other hand, this study showed that there were no differences in the relative expression of the adhesion molecules (N-cadherin and E-cadherin) between the frontal, parietal and interparietal bones.
This study provides evidence that different calvarial bones—frontal, parietal and interparietal—contain different bone microenvironments. This provides platforms for pre-clinical investigations to study myeloma colonization and growth in the skull (calvarial bone) as well as anti-myeloma therapies, particularly those targeting myeloma bone disease.
Radl J, Hollander CF, van den Berg p, de Glopper E. Idiopathic paraproteinaemia. I. Studies in an animal model—the ageing C57BL/KaLwRij mouse. Clin Exp Immunol. 1978;33:395–402.
Radl J, De Glopper ED, Schuit HR, Zurcher C. Idiopathic paraproteinemia. II. Transplantation of the paraprotein-producing clone from old to young C57BL/KaLwRij mice. J Immunol. 1979;122:609–13.
Vanderkerken K, De Raeve H, Goes E, Van Meirvenne S, Radl J, Van Riet I, et al. Organ involvement and phenotypic adhesion profile of 5T2 and 5T33 myeloma cells in the C57BL/KaLwRij mouse. Br J Cancer. 1997;76:451–60.
Vanderkerken K, De Greef C, Asosingh K, Arteta B, De Veerman M, Vande Broek I, et al. Selective initial in vivo homing pattern of 5T2 multiple myeloma cells in the C57BL/KaLwRij mouse. Br J Cancer. 2000;82:953–9.
Croucher PI, De Hendrik R, Perry MJ, Hijzen A, Shipman CM, Lippitt J, et al. Zoledronic acid treatment of 5T2MM-bearing mice inhibits the development of myeloma bone disease: evidence for decreased osteolysis, tumor burden and angiogenesis, and increased survival. J Bone Miner Res. 2003;18:482–92.
Alici E, Konstantinidis KV, Aints A, Dilber MS, Abedi-Valugerdi M. Visualization of 5T33 myeloma cells in the C57BL/KaLwRij mouse: establishment of a new syngeneic murine model of multiple myeloma. Exp Hematol. 2004;32:1064–72.
Olson DL, Burkly LC, Leone DR, Dolinski BM, Lobb RR. Anti-alpha4 integrin monoclonal antibody inhibits multiple myeloma growth in a murine model. Mol Cancer Ther. 2005;4:91–9.
Edwards CM, Lwin ST, Fowler JA, Oyajobi BO, Zhuang J, Bates AL, et al. Myeloma cells exhibit an increase in proteasome activity and an enhanced response to proteasome inhibition in the bone marrow microenvironment in vivo. Am J Hematol. 2009;84:268–72.
Fowler JA, Mundy GR, Lwin ST, Edwards CM. Bone marrow stromal cells create a permissive microenvironment for myeloma development: a new stromal role for Wnt inhibitor Dkk1. Cancer Res. 2012;72:2183–9.
Van Valckenborgh E, Matsui W, Agarwal P, Lub S, Dehui X, De Bruyne E, et al. Tumor-initiating capacity of CD138− and CD138+ tumor cells in the 5T33 multiple myeloma model. Leukemia. 2012;26:1436–9.
Al-Amer O. Tracking of individual myeloma cell homing to the calvarial bone marrow using myeloma murine models. Int Res J Med Med Sci. 2015;3:44–50.
Al-Amer O. Homing of myeloma cells into calvarial bone using myeloma murine models. Int J Health Sci Res. 2015;5:357–63.
Chen D, Zhao M, Mundy GR. Bone morphogenetic proteins. Growth Factors. 2004;22:233–41.
Milat F, Ng KW. Is Wnt signalling the final common pathway leading to bone formation? Mol Cell Endocrinol. 2009;310:52–62.
Engin F, Lee B. NOTCHing the bone: insights into multi-functionality. Bone. 2010;46:274–80.
Lian JB, Stein GS. Runx2/Cbfa1: a multifunctional regulator of bone formation. Curr Pharm Des. 2003;9:2677–85.
Di Lullo GA, Sweeney SM, Korkko J, Ala-Kokko L, San Antonio JD. Mapping the ligand-binding sites and disease-associated mutations on the most abundant protein in the human, type I collagen. J Biol Chem. 2002;277:4223–31.
Boyle WJ, Simonet WS, Lacey DL. Osteoclast differentiation and activation. Nature. 2003;423:337–42.
Lacey DL, Timms E, Tan HL, Kelley MJ, Dunstan CR, Burgess T, et al. Osteoprotegerin ligand is a cytokine that regulates osteoclast differentiation and activation. Cell. 1998;93:165–76.
Yasuda H, Shima N, Nakagawa N, Mochizuki SI, Yano K, Fujise N, et al. Identity of osteoclastogenesis inhibitory factor (OCIF) and osteoprotegerin (OPG): a mechanism by which OPG/OCIF inhibits osteoclastogenesis in vitro. Endocrinology. 1998;139:1329–37.
Grimaud E, Redini F, Heymann D. Osteoprotegerin: a new therapeutic agent for the treatment of bone disease. Drug Discov Today. 2001;6:1241–2.
Bucay N, Sarosi I, Dunstan CR, Morony S, Tarpley J, Capparelli C, et al. Osteoprotegerin-deficient mice develop early onset osteoporosis and arterial calcification. Genes Dev. 1998;12:1260–8.
Zhang J, Niu C, Ye L, Huang H, He X, Tong W-G, et al. Identification of the haematopoietic stem cell niche and control of the niche size. Nature. 2003;425:836–41.
Xie Y, Yin T, Wiegraebe W, He XC, Miller D, Stark D, et al. Detection of functional haematopoietic stem cell niche using real-time imaging. Nature. 2009;457:97–101.
Takeichi M. Morphogenetic roles of classic cadherins. Curr Opin Cell Biol. 1995;7:619–27.
Gumbiner BM, Yamada KM. Cell-to-cell contact and extracellular matrix. Curr Opin Cell Biol. 1995;7:615–8.
Iseki S, Wilkie AO, Morriss-Kay GM. Fgfr1 and Fgfr2 have distinct differentiation- and proliferation-related roles in the developing mouse skull vault. Development. 1999;126:5611–20.
Wilkie AO, Patey SJ, Kan SH, van den Ouweland AM, Hamel BC. Fgfs, their receptors, and human limb malformations: clinical and molecular correlations. Am J Med Genet. 2002;112:266–78.
Ornitz DM, Marie PJ. FGF signaling pathways in endochondral and intramembranous bone development and human genetic disease. Genes Dev. 2002;16:1446–65.
Fischer AH, Jacobson KA, Rose J, Zeller R. Hematoxylin and eosin staining of tissue and cell sections. CSH Protoc. 2008;2008:pdb prot4986.
Llewellyn BD. Nuclear staining with alum hematoxylin. Biotech Histochem. 2009;84:159–77.
Parfitt AM, Drezner MK, Glorieux FH, Kanis JA, Malluche H, Meunier PJ, et al. Bone histomorphometry: standardization of nomenclature, symbols, and units. Report of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res. 1987;2:595–610.
Heider U, Fleissner C, Zavrski I, Kaiser M, Hecht M, Jakob C, et al. Bone markers in multiple myeloma. Eur J Cancer. 2006;42:1544–53.
Esteve FR, Roodman GD. Pathophysiology of myeloma bone disease. Best Pract Res Clin Haematol. 2007;20:613–24.
Datta HK, Ng WF, Walker JA, Tuck SP, Varanasi SS. The cell biology of bone metabolism. J Clin Pathol. 2008;61:577–87.
Edwards CM, Zhuang J, Mundy GR. The pathogenesis of the bone disease of multiple myeloma. Bone. 2008;42:1007–13.
Thanks to The University of Sheffield for accepting me as a visiting researcher to conduct my own personal research and to meet with other academics to exchange ideas and research findings.
The author declares no competing interests.
Availability of data and materials
The raw data is not shared since they are saved in specific format of Power graph software.
Ethics approval and consent to participate
All procedures involving animals were approved by the Home Office (PIL 40/10118) and the University of Sheffield’s Animal Ethics Committee.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.