MarrowQuant Across Aging and Aplasia: A Digital Pathology Workflow for Quantification of Bone Marrow Compartments in Histological Sections

The bone marrow (BM) exists heterogeneously as hematopoietic/red or adipocytic/yellow marrow depending on skeletal location, age, and physiological condition. Mouse models and patients undergoing radio/chemotherapy or suffering acute BM failure endure rapid adipocytic conversion of the marrow microenvironment, the so-called "red-to-yellow" transition. Following hematopoietic recovery, such as upon BM transplantation, a "yellow-to-red" transition occurs and functional hematopoiesis is restored. Gold Standards to estimate BM cellular composition are pathologists' assessment of hematopoietic cellularity in hematoxylin and eosin (H&E) stained histological sections as well as volumetric measurements of marrow adiposity with contrast-enhanced micro-computerized tomography (CE-μCT) upon osmium-tetroxide lipid staining. Due to user-dependent variables, reproducibility in longitudinal studies is a challenge for both methods. Here we report the development of a semi-automated image analysis plug-in, MarrowQuant, which employs the open-source software QuPath, to systematically quantify multiple bone components in H&E sections in an unbiased manner. MarrowQuant discerns and quantifies the areas occupied by bone, adipocyte ghosts, hematopoietic cells, and the interstitial/microvascular compartment. A separate feature, AdipoQuant, fragments adipocyte ghosts in H&E-stained sections of extramedullary adipose tissue to render adipocyte area and size distribution. Quantification of BM hematopoietic cellularity with MarrowQuant lies within the range of scoring by four independent pathologists, while quantification of the total adipocyte area in whole bone sections compares with volumetric measurements. Employing our tool, we were able to develop a standardized map of BM hematopoietic cellularity and adiposity in mid-sections of murine C57BL/6 bones in homeostatic conditions, including quantification of the highly predictable red-to-yellow transitions in the proximal section of the caudal tail and in the proximal-to-distal tibia. Additionally, we present a comparative skeletal map induced by lethal irradiation, with longitudinal quantification of the "red-to-yellow-to-red" transition over 2 months in C57BL/6 femurs and tibiae. We find that, following BM transplantation, BM adiposity inversely correlates with kinetics of hematopoietic recovery and that a proximal to distal gradient is conserved. Analysis of in vivo recovery through magnetic resonance imaging (MRI) reveals comparable kinetics. On human trephine biopsies MarrowQuant successfully recognizes the BM compartments, opening avenues for its application in experimental, or clinical contexts that require standardized human BM evaluation.