Apoptotic cell fragments locally activate tingible body macrophages in the germinal center.

Germinal centers (GCs) that form within lymphoid follicles during antibody responses are sites of massive cell death. Tingible body macrophages (TBMs) are tasked with apoptotic cell clearance to prevent secondary necrosis and autoimmune activation by intracellular self antigens. We show by multiple redundant and complementary methods that TBMs derive from a lymph node-resident, CD169-lineage, CSF1R-blockade-resistant precursor that is prepositioned in the follicle. Non-migratory TBMs use cytoplasmic processes to chase and capture migrating dead cell fragments using a "lazy" search strategy. Follicular macrophages activated by the presence of nearby apoptotic cells can mature into TBMs in the absence of GCs. Single-cell transcriptomics identified a TBM cell cluster in immunized lymph nodes which upregulated genes involved in apoptotic cell clearance. Thus, apoptotic B cells in early GCs trigger activation and maturation of follicular macrophages into classical TBMs to clear apoptotic debris and prevent antibody-mediated autoimmune diseases.

Osteoclasts recycle via osteomorphs during RANKL-stimulated bone resorption.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases.

ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate.

MOTIVATION: Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory. RESULTS: We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation. AVAILABILITY AND IMPLEMENTATION: ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

Tracking the clonal dynamics of SARS-CoV-2-specific T cells in children and adults with mild/asymptomatic COVID-19.

Children infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) develop less severe coronavirus disease 2019 (COVID-19) than adults. The mechanisms for the age-specific differences and the implications for infection-induced immunity are beginning to be uncovered. We show by longitudinal multimodal analysis that SARS-CoV-2 leaves a small footprint in the circulating T cell compartment in children with mild/asymptomatic COVID-19 compared to adult household contacts with the same disease severity who had more evidence of systemic T cell interferon activation, cytotoxicity and exhaustion. Children harbored diverse polyclonal SARS-CoV-2-specific naïve T cells whereas adults harbored clonally expanded SARS-CoV-2-specific memory T cells. A novel population of naïve interferon-activated T cells is expanded in acute COVID-19 and is recruited into the memory compartment during convalescence in adults but not children. This was associated with the development of robust CD4+ memory T cell responses in adults but not children. These data suggest that rapid clearance of SARS-CoV-2 in children may compromise their cellular immunity and ability to resist reinfection.

Spliceosome mutations are associated with clinical response in a phase 1b/2 study of the PLK1 inhibitor onvansertib in combination with decitabine in relapsed or refractory acute myeloid leukemia.

PLK1 is overexpressed in acute myeloid leukemia (AML). A phase 1b trial of the PLK1 inhibitor onvansertib (ONV) combined with decitabine (DAC) demonstrated initial safety and efficacy in patients with relapsed/refractory (R/R) AML. The current study aimed to identify molecular predictors of response to ONV + DAC in R/R AML patients. A total of 44 R/R AML patients were treated with ONV + DAC and considered evaluable for efficacy. Bone marrow (BM) samples were collected at baseline for genomic and transcriptomic analysis (n = 32). A 10-gene expression signature, predictive of response to ONV + DAC, was derived from the leading-edge genes of gene set enrichment analyses (GSEA). The gene signature was evaluated in independent datasets and used to identify associated mutated genes. Twenty percent of the patients achieved complete remission, with or without hematologic count recovery (CR/CRi), and 32% exhibited a ≥50% reduction in bone marrow blasts. Patients who responded to treatment had elevated mitochondrial function and OXPHOS. The gene signature was not associated with response to DAC alone in an independent dataset. By applying the signature to the BeatAML cohort (n = 399), we identified a positive association between predicted ONV + DAC response and mutations in splicing factors (SF). In the phase 1b/2 trial, patients with SF mutations (SRSF2, SF3B1) had a higher CR/CRi rate (50%) compared to those without SF mutations (9%). PLK1 inhibition with ONV in combination with DAC could be a potential therapy in R/R AML patients, particularly those with high OXPHOS gene expression and SF mutations.

Regenerating zebrafish scales express a subset of evolutionary conserved genes involved in human skeletal disease.

BACKGROUND: Scales are mineralised exoskeletal structures that are part of the dermal skeleton. Scales have been mostly lost during evolution of terrestrial vertebrates whilst bony fish have retained a mineralised dermal skeleton in the form of fin rays and scales. Each scale is a mineralised collagen plate that is decorated with both matrix-building and resorbing cells. When removed, an ontogenetic scale is quickly replaced following differentiation of the scale pocket-lining cells that regenerate a scale. Processes promoting de novo matrix formation and mineralisation initiated during scale regeneration are poorly understood. Therefore, we performed transcriptomic analysis to determine gene networks and their pathways involved in dermal scale regeneration. RESULTS: We defined the transcriptomic profiles of ontogenetic and regenerating scales of zebrafish and identified 604 differentially expressed genes (DEGs). These were enriched for extracellular matrix, ossification, and cell adhesion pathways, but not in enamel or dentin formation processes indicating that scales are reminiscent to bone. Hypergeometric tests involving monogenetic skeletal disorders showed that DEGs were strongly enriched for human orthologues that are mutated in low bone mass and abnormal bone mineralisation diseases (P< 2× 10-3). The DEGs were also enriched for human orthologues associated with polygenetic skeletal traits, including height (P< 6× 10-4), and estimated bone mineral density (eBMD, P< 2× 10-5). Zebrafish mutants of two human orthologues that were robustly associated with height (COL11A2, P=6× 10-24) or eBMD (SPP1, P=6× 10-20) showed both exo- and endo- skeletal abnormalities as predicted by our genetic association analyses; col11a2Y228X/Y228X mutants showed exoskeletal and endoskeletal features consistent with abnormal growth, whereas spp1P160X/P160X mutants predominantly showed mineralisation defects. CONCLUSION: We show that scales have a strong osteogenic expression profile comparable to other elements of the dermal skeleton, enriched in genes that favour collagen matrix growth. Despite the many differences between scale and endoskeletal developmental processes, we also show that zebrafish scales express an evolutionarily conserved sub-population of genes that are relevant to human skeletal disease.

Gene expression predicts dormant metastatic breast cancer cell phenotype.

BACKGROUND: Breast cancer can recur months to decades after an initial diagnosis and treatment. The mechanisms that control tumor cell dormancy remain poorly understood, making it difficult to predict which patients will recur and thus benefit from more rigorous screening and treatments. Unfortunately, the extreme rarity of dormant DTCs has been a major obstacle to their study. METHODS: To overcome this challenge, we developed an efficient system to isolate and study rare dormant breast cancer cells from metastatic organs including bones, which represent a major site of metastasis. After isolation of cells from the long bones, we used single cell RNA-sequencing (scRNA-seq) to profile proliferative and dormant PyMT-Bo1 breast cancer cells. We also compared this signature to dormant versus proliferative tumor cells isolated from the lungs. Finally, we compared our dormant signature to human datasets. RESULTS: We identified a group of genes including Cfh, Gas6, Mme and Ogn that were highly expressed in dormant breast cancer cells present in the bone and lung. Expression of these genes had no impact on dormancy in murine models, but their expression correlated with disease-free survival in primary human breast cancer tumors, suggesting that these genes have predictive value in determining which patients are likely to recur. CONCLUSIONS: Dormant breast cancer cells exhibit a distinct gene expression signature regardless of metastatic site. Genes enriched in dormant breast cancer cells correlate with recurrence-free survival in breast cancer patients.

Prostate cancer cell-intrinsic interferon signaling regulates dormancy and metastatic outgrowth in bone.

The latency associated with bone metastasis emergence in castrate-resistant prostate cancer is attributed to dormancy, a state in which cancer cells persist prior to overt lesion formation. Using single-cell transcriptomics and ex vivo profiling, we have uncovered the critical role of tumor-intrinsic immune signaling in the retention of cancer cell dormancy. We demonstrate that loss of tumor-intrinsic type I IFN occurs in proliferating prostate cancer cells in bone. This loss suppresses tumor immunogenicity and therapeutic response and promotes bone cell activation to drive cancer progression. Restoration of tumor-intrinsic IFN signaling by HDAC inhibition increased tumor cell visibility, promoted long-term antitumor immunity, and blocked cancer growth in bone. Key findings were validated in patients, including loss of tumor-intrinsic IFN signaling and immunogenicity in bone metastases compared to primary tumors. Data herein provide a rationale as to why current immunotherapeutics fail in bone-metastatic prostate cancer, and provide a new therapeutic strategy to overcome the inefficacy of immune-based therapies in solid cancers.

A niche-dependent myeloid transcriptome signature defines dormant myeloma cells.

The era of targeted therapies has seen significant improvements in depth of response, progression-free survival, and overall survival for patients with multiple myeloma. Despite these improvements in clinical outcome, patients inevitably relapse and require further treatment. Drug-resistant dormant myeloma cells that reside in specific niches within the skeleton are considered a basis of disease relapse but remain elusive and difficult to study. Here, we developed a method to sequence the transcriptome of individual dormant myeloma cells from the bones of tumor-bearing mice. Our analyses show that dormant myeloma cells express a distinct transcriptome signature enriched for immune genes and, unexpectedly, genes associated with myeloid cell differentiation. These genes were switched on by coculture with osteoblastic cells. Targeting AXL, a gene highly expressed by dormant cells, using small-molecule inhibitors released cells from dormancy and promoted their proliferation. Analysis of the expression of AXL and coregulated genes in human cohorts showed that healthy human controls and patients with monoclonal gammopathy of uncertain significance expressed higher levels of the dormancy signature genes than patients with multiple myeloma. Furthermore, in patients with multiple myeloma, the expression of this myeloid transcriptome signature translated into a twofold increase in overall survival, indicating that this dormancy signature may be a marker of disease progression. Thus, engagement of myeloma cells with the osteoblastic niche induces expression of a suite of myeloid genes that predicts disease progression and that comprises potential drug targets to eradicate dormant myeloma cells.

Expression of the chemokine receptor CCR1 promotes the dissemination of multiple myeloma plasma cells in vivo

Multiple myeloma (MM) disease progression is dependent on the ability of MM plasma cells (PCs) to egress from the bone marrow (BM), enter the circulation and disseminate to distal BM sites. Expression of the chemokine CXCL12 by BM stromal cells is crucial for MM PC retention within the BM. However, the mechanisms which overcome CXCL12-mediated retention to enable dissemination are poorly understood. We have previously identified that treatment with the CCR1 ligand CCL3 inhibits the response to CXCL12 in MM cell lines, suggesting that CCL3/CCR1 signalling may enable egress of MM PC from the BM. Here, we demonstrated that CCR1 expression was an independent prognostic indicator in newly diagnosed MM patients. Furthermore, we showed that CCR1 is a crucial driver of dissemination in vivo, with CCR1 expression in the murine MM cell line 5TGM1 being associated with an increased incidence of bone and splenic disseminated tumours in C57BL/KaLwRij mice. Furthermore, we demonstrated that CCR1 knockout in the human myeloma cell line OPM2 resulted in a >95% reduction in circulating MM PC numbers and BM and splenic tumour dissemination following intratibial injection in NSG mice. Therapeutic targeting of CCR1 with the inhibitor CCX9588 significantly reduced OPM2 or RPMI-8226 dissemination in intratibial xenograft models. Collectively, our findings suggest a novel role for CCR1 as a critical driver of BM egress of MM PCs during tumour dissemination. Furthermore, these data suggest that CCR1 may represent a potential therapeutic target for the prevention of MM tumour dissemination.

Inhibiting the osteocyte-specific protein sclerostin increases bone mass and fracture resistance in multiple myeloma.

Multiple myeloma (MM) is a plasma cell cancer that develops in the skeleton causing profound bone destruction and fractures. The bone disease is mediated by increased osteoclastic bone resorption and suppressed bone formation. Bisphosphonates used for treatment inhibit bone resorption and prevent bone loss but fail to influence bone formation and do not replace lost bone, so patients continue to fracture. Stimulating bone formation to increase bone mass and fracture resistance is a priority; however, targeting tumor-derived modulators of bone formation has had limited success. Sclerostin is an osteocyte-specific Wnt antagonist that inhibits bone formation. We hypothesized that inhibiting sclerostin would prevent development of bone disease and increase resistance to fracture in MM. Sclerostin was expressed in osteocytes from bones from naive and myeloma-bearing mice. In contrast, sclerostin was not expressed by plasma cells from 630 patients with myeloma or 54 myeloma cell lines. Mice injected with 5TGM1-eGFP, 5T2MM, or MM1.S myeloma cells demonstrated significant bone loss, which was associated with a decrease in fracture resistance in the vertebrae. Treatment with anti-sclerostin antibody increased osteoblast numbers and bone formation rate but did not inhibit bone resorption or reduce tumor burden. Treatment with anti-sclerostin antibody prevented myeloma-induced bone loss, reduced osteolytic bone lesions, and increased fracture resistance. Treatment with anti-sclerostin antibody and zoledronic acid combined increased bone mass and fracture resistance when compared with treatment with zoledronic acid alone. This study defines a therapeutic strategy superior to the current standard of care that will reduce fractures for patients with MM.

A Microsatellite Analysis Used to Identify Global Pathways of Movement of Phytophthora cinnamomi and the Likely Sources of Wildland Infestations in California and Mexico.

The genetic structure of a sample of isolates of the oomycete plant pathogen Phytophthora cinnamomi from natural and agricultural outbreaks and the long-distance movement of individual genotypes were studied using four microsatellite markers to genotype 159 isolates of Californian, Mexican, and worldwide origins. Allelic profiles identified 75 multilocus genotypes. A STRUCTURE analysis placed them in three groups characterized by different geographic and host ranges, different genic and genotypic diversity, and different reproductive modes. When relationships among genotypes were visualized on a minimum spanning network (MSN), genotypes belonging to the same STRUCTURE group were contiguous, with rare exceptions. A putatively ancestral group 1 had high genic diversity, included all A1 mating type isolates and all Papuan isolates in the sample, was rarely isolated from natural settings in California and Mexico, and was positioned at the center of the MSN. Putatively younger groups 2 and 3 had lower genic diversity, were both neighbors to group 1 but formed two distinct peripherical sectors of the MSN, and were equally present in agricultural commodities and natural settings in Mexico and California. A few genotypes, especially in groups 2 and 3, were isolated multiple times in different locations and settings. The presence of identical genotypes from the same hosts in different continents indicated that long-distance human-mediated movement of P. cinnamomi had occurred. The presence of identical genotypes at high frequencies in neighboring wildlands and agricultural settings suggest that specific commodities may have been the source of recent wild infestations caused by novel invasive genotypes.

The skeletal cell-derived molecule sclerostin drives bone marrow adipogenesis.

The bone marrow niche is a dynamic and complex microenvironment that can both regulate, and be regulated by the bone matrix. Within the bone marrow (BM), mesenchymal stromal cell (MSC) precursors reside in a multi-potent state and retain the capacity to differentiate down osteoblastic, adipogenic, or chondrogenic lineages in response to numerous biochemical cues. These signals can be altered in various pathological states including, but not limited to, osteoporotic-induced fracture, systemic adiposity, and the presence of bone-homing cancers. Herein we provide evidence that signals from the bone matrix (osteocytes) determine marrow adiposity by regulating adipogenesis in the bone marrow. Specifically, we found that physiologically relevant levels of Sclerostin (SOST), which is a Wnt-inhibitory molecule secreted from bone matrix-embedded osteocytes, can induce adipogenesis in 3T3-L1 cells, mouse ear- and BM-derived MSCs, and human BM-derived MSCs. We demonstrate that the mechanism of SOST induction of adipogenesis is through inhibition of Wnt signaling in pre-adipocytes. We also demonstrate that a decrease of sclerostin in vivo, via both genetic and pharmaceutical methods, significantly decreases bone marrow adipose tissue (BMAT) formation. Overall, this work demonstrates a direct role for SOST in regulating fate determination of BM-adipocyte progenitors. This provides a novel mechanism for which BMAT is governed by the local bone microenvironment, which may prove relevant in the pathogenesis of certain diseases involving marrow adipose. Importantly, with anti-sclerostin therapy at the forefront of osteoporosis treatment and a greater recognition of the role of BMAT in disease, these data are likely to have important clinical implications.

Homozygous Dkk1 Knockout Mice Exhibit High Bone Mass Phenotype Due to Increased Bone Formation.

Wnt antagonist Dkk1 is a negative regulator of bone formation and Dkk1 +/- heterozygous mice display a high bone mass phenotype. Complete loss of Dkk1 function disrupts embryonic head development. Homozygous Dkk1 -/- mice that were heterozygous for Wnt3 loss of function mutation (termed Dkk1 KO) are viable and allowed studying the effects of homozygous inactivation of Dkk1 on bone formation. Dkk1 KO mice showed a high bone mass phenotype exceeding that of heterozygous mice as well as a high incidence of polydactyly and kinky tails. Whole body bone density was increased in the Dkk1 KO mice as shown by longitudinal dual-energy X-ray absorptiometry. MicroCT analysis of the distal femur revealed up to 3-fold increases in trabecular bone volume and up to 2-fold increases in the vertebrae, compared to wild type controls. Cortical bone was increased in both the tibiae and vertebrae, which correlated with increased strength in tibial 4-point bending and vertebral compression tests. Dynamic histomorphometry identified increased bone formation as the mechanism underlying the high bone mass phenotype in Dkk1 KO mice, with no changes in bone resorption. Mice featuring only Wnt3 heterozygosity showed no evident bone phenotype. Our findings highlight a critical role for Dkk1 in the regulation of bone formation and a gene dose-dependent response to loss of DKK1 function. Targeting Dkk1 to enhance bone formation offers therapeutic potential for osteoporosis.

Parathyroid Hormone (PTH) Increases Skeletal Tumour Growth and Alters Tumour Distribution in an In Vivo Model of Breast Cancer.

Breast cancer cells colonize the skeleton by homing to specific niches, but the involvement of osteoblasts in tumour cell seeding, colonization, and progression is unknown. We used an in vivo model to determine how increasing the number of cells of the osteoblast lineage with parathyroid hormone (PTH) modified subsequent skeletal colonization by breast cancer cells. BALB/c nude mice were injected for five consecutive days with PBS (control) or PTH and then injected with DiD-labelled breast cancer cells via the intra-cardiac route. Effects of PTH on the bone microenvironment and tumour cell colonization and growth was analyzed using bioluminescence imaging, two-photon microscopy, and histological analysis. PTH treatment caused a significant, transient increase in osteoblast numbers compared to control, whereas bone volume/structure in the tibia was unaffected. There were no differences in the number of tumour cells seeding to the tibias, or in the number of tumours in the hind legs, between the control and PTH group. However, animals pre-treated with PTH had a significantly higher number of tumour colonies distributed throughout skeletal sites outside the hind limbs. This is the first demonstration that PTH-induced stimulation of osteoblastic cells may result in alternative skeletal sites becoming available for breast cancer cell colonization.

Myeloma impairs mature osteoblast function but causes early expansion of osteo-progenitors: temporal changes in bone physiology and gene expression in the KMS12BM model.

Myeloma bone disease results from an uncoupling of osteoclastic resorption and osteoblastic bone formation, but early changes in osteogenic function remain poorly defined. We used the KMS12BM xenograft model to investigate cellular and molecular events at early and late stages of disease. Lytic lesions and changes in osteoblast and osteoclast numbers occur late (8 weeks), however, micro-computed tomography of femora revealed significant reduction in bone volume at earlier disease stages (3 weeks) when tumour burden is low. Calcein labelling demonstrated reduced mineralization and bone formation at 3 weeks, suggesting functional impairment despite preserved osteoblast numbers. Osteo-progenitors from compact bone increased early (1 week), but fell at 3 weeks and were profoundly suppressed by 8 weeks. Exposure of osteoblast progenitors to multiple myeloma (MM) cells in vitro induced cell cycling, suggesting a mechanistic basis for early expansion of osteo-progenitors. We observed temporal changes in chemokine, osteogenic and osteoclastogenic genes in the stromal compartment. Notably, an early rise in CCL3 may underlie functional changes in mature osteoblasts at 3 weeks. Our data indicate that MM has distinct effects on mature osteoblasts and immature osteo-progenitors. Our findings argue for early clinical intervention to prevent bone changes that ultimately lead to the development of osteolytic disease.

Mitotic quiescence, but not unique "stemness," marks the phenotype of bone metastasis-initiating cells in prostate cancer.

This study aimed to identify subpopulations of prostate cancer cells that are responsible for the initiation of bone metastases. Using rapidly dividing human prostate cancer cell lines, we identified mitotically quiescent subpopulations (<1%), which we compared with the rapidly dividing populations for patterns of gene expression and for their ability to migrate to the skeletons of athymic mice. The study used 2-photon microscopy to map the presence/distribution of fluorescently labeled, quiescent cells and luciferase expression to determine the presence of growing bone metastases. We showed that the mitotically quiescent cells were very significantly more tumorigenic in forming bone metastases than fast-growing cells (55 vs. 15%) and had a unique gene expression profile. The quiescent cells were not uniquely stem cell like, with no expression of CD133 but had the same level expression of other putative prostate stem cell markers (CD44 and integrins α2/ß1), when compared to the rapidly proliferating population. In addition, mitotic quiescence was associated with very high levels of C-X-C chemokine receptor type 4 (CXCR4) production. Inhibition of CXCR4 activity altered the homing of quiescent tumor cells to bone. Our studies suggest that mitotic dormancy is a unique phenotype that facilitates tumor cell colonization of the skeleton in prostate cancer.

OPG-Fc inhibits ovariectomy-induced growth of disseminated breast cancer cells in bone.

Dormant disseminated tumour cells can be detected in the bone marrow of breast cancer patients several years after resection of the primary tumour. The majority of these patients will remain asymptomatic, however, ∼ 15% will go on to develop overt bone metastases and this condition is currently incurable. The reason why these dormant cells are stimulated to proliferate and form bone tumours in some patients and not others remains to be elucidated. We have recently shown that in an in vivo model, increasing bone turnover by ovariectomy stimulated proliferation of disseminated tumour cells, resulting in formation of bone metastasis. We now show for the first time that osteoclast mediated mechanisms induce growth of tumours from dormant MDA-MB-231 cells disseminated in the bone. We also show that disruption of RANK-RANKL interactions following administration of OPG-Fc inhibits growth of these dormant tumour cells in vivo. Our data support early intervention with anti-resorptive therapy in a low-oestrogen environment to prevent development of bone metastases.

Rapid-throughput skeletal phenotyping of 100 knockout mice identifies 9 new genes that determine bone strength.

Osteoporosis is a common polygenic disease and global healthcare priority but its genetic basis remains largely unknown. We report a high-throughput multi-parameter phenotype screen to identify functionally significant skeletal phenotypes in mice generated by the Wellcome Trust Sanger Institute Mouse Genetics Project and discover novel genes that may be involved in the pathogenesis of osteoporosis. The integrated use of primary phenotype data with quantitative x-ray microradiography, micro-computed tomography, statistical approaches and biomechanical testing in 100 unselected knockout mouse strains identified nine new genetic determinants of bone mass and strength. These nine new genes include five whose deletion results in low bone mass and four whose deletion results in high bone mass. None of the nine genes have been implicated previously in skeletal disorders and detailed analysis of the biomechanical consequences of their deletion revealed a novel functional classification of bone structure and strength. The organ-specific and disease-focused strategy described in this study can be applied to any biological system or tractable polygenic disease, thus providing a general basis to define gene function in a system-specific manner. Application of the approach to diseases affecting other physiological systems will help to realize the full potential of the International Mouse Phenotyping Consortium.