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.
Published in the journal:
. PLoS Genet 8(8): e32767. doi:10.1371/journal.pgen.1002858
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1002858
Summary
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.
Introduction
Studies of extreme phenotypes in humans have been instrumental in identifying molecular mechanisms underlying rare single gene disorders as well as common and chronic diseases including diabetes and obesity. Such studies have resulted in novel treatments that revolutionize the lives of affected individuals [1]–[4]. Collection of suitable cohorts, however, is expensive and takes many years to achieve, and progress has been limited to conditions in which simple and quantitative phenotypes can be defined [4]–[7]. By analogy, we hypothesized that an organ-specific extreme phenotype screen in knockout mice would more rapidly identify new genetic determinants of disease and also provide in vivo models to elucidate their molecular basis. The International Knockout Mouse Consortium (IKMC) has now established an ideal resource of mutant ES cells to test this hypothesis [8], [9].
The skeleton represents a paradigm organ system and osteoporosis is an important global disease ideally suited to such an approach. Osteoporosis is the commonest skeletal disorder affecting hundreds of millions of people worldwide and costing tens of billions of pounds each year [10]. Between 50 and 85% of the variance in bone mineral density (BMD) is genetically determined [11], but only 3% is accounted for by known genetic variation [12] and the vast majority of genes involved remain to be identified [13]. Current treatments reduce fracture risk by only 25–50% [14], [15] and thus there is urgent need to define new pathways that regulate bone turnover and strength in order to identify novel therapeutic targets. Accordingly, application of an extreme phenotype approach to study skeletal disorders in humans has already led to discovery of SOST (ENSG00000167941) and LRP5 (ENSG00000162337) as critical regulators of Wnt signaling in bone [5], [7], [16] and resulted in development of new drugs to stimulate bone formation [17].
The Wellcome Trust Sanger Institute Mouse Genetics Project (MGP) is undertaking high-throughput production of knockout mice using targeted ES cells generated by the IKMC. Knockout mice are generated using a knockout-first conditional gene targeting strategy (Figure S1), in which expression from the targeted allele can be investigated by X-gal staining for LacZ gene expression [9] (Figure S2). Each mouse undergoes a broad-based primary screen to identify developmental, anatomical, physiological and behavioral phenotypes [18]–[20]. A critical challenge now is to enhance this initial screening by developing organ- or disease-specific approaches [21] that are essential to identify biologically significant and functionally relevant phenotypes rapidly and cost-effectively for the benefit of the scientific community [18], [19], [21].
We, therefore, developed high-throughput skeletal phenotyping methods and prospectively studied 100 consecutive unselected mutant strains from the MGP. Using this approach, we discovered nine new genetic determinants of bone mass and strength and identified a novel functional classification of bone structure. These conditional knockout mice [9] can now be used to investigate cell-specific gene function and identify new regulatory pathways in the skeleton. The strategy can be applied to other physiological systems and complex diseases, thus realizing the full potential of the International Mouse Phenotyping Consortium.
Results
Broad primary phenotype screening
Mice generated by the MGP pipeline undergo a broad primary phenotype screen followed by terminal collection of blood and tissue at 16 weeks of age [20]. The screen is conducted on viable homozygote mutants, or heterozygotes in cases of embryonic lethality, and reports 233 variables relating to 28 physiological systems that include embyrogenesis; reproduction; growth; neurological; behaviour; sensory; skeleton; muscle; gastrointestinal and hepatobiliary; cardiovascular; endocrine; adipose; metabolism; haematopoietic; immune; skin and pigmentation; respiratory; and renal. Parameters relevant to the skeleton include body length, x-ray skeletal survey, dual energy x-ray absorptiometry (DEXA) analysis of BMD and biochemical measures of mineral metabolism. Tissues in which the targeted gene is expressed are determined by staining for lacZ reporter gene expression in heterozygous mice (Figure S2). To extend this broad initial screen, we incorporated novel imaging, statistical and biomechanical approaches for the specific and sensitive detection of functionally important skeletal abnormalities (Figure S3).
Generation of normal reference data
In order to establish strain-specific reference ranges for these new approaches, limbs from 16 week-old female C57BL/6 (B6Brd;B6Dnk;B6N-Tyrc-Brd) wild-type mice (n = 77) were obtained from 18 control cohorts. Normal ranges for six independent parameters of bone structure were obtained using Faxitron x-ray point projection microradiography and micro-computed tomography (micro-CT) (Fig. 1 and Figure S4). Bone mineral content (BMC), bone length and cortical bone thickness were determined by x-ray microradiography and measures of trabecular bone volume per tissue volume (BV/TV), trabecular number (Tb.N) and trabecular thickness (Tb.Th) by micro-CT (Fig. 1). Reference data were also obtained for six biomechanical parameters (Fig. 2 and Figure S5). The yield, maximum and fracture loads, stiffness and the proportions of energy dissipated prior to maximum load and fracture were determined from load displacement curves obtained in destructive 3-point bend tests (Fig. 2).
Validation of phenotyping methods
Limbs from 16 week-old female knockout mice in an identical C57BL/6 genetic background were obtained prospectively from the MGP pipeline (n = 100 unselected knockout strains, 2–6 mice per strain) and analyzed by x-ray microradiography, micro-CT and 3-point bend testing. Serendipitously, one of the 100 unselected strains was a homozygous knockout of Sparc (ENSMUSG00000018593), which encodes the extracellular matrix glycoprotein osteonectin. Deletion of Sparc is known to cause low bone turnover osteopenia resulting in weak and brittle bone with a higher mineral-to-matrix ratio due to reduced bone matrix content [22], [23]. Thus, Sparc knockout mice represented a well-characterized positive control for validation of our approach. Consistent with the reported phenotype [22], [23], we identified that Sparc knockout mice had reduced BMD and BMC with loss of trabecular bone but preservation of cortical bone thickness (Fig. 1 and Table S1), resulting in weak and brittle bone of reduced stiffness (Fig. 2). We also identified short stature in Sparc knockout mice (Fig. 1A), a parameter not investigated in previous studies. These findings validate the use of complementary and multi-parameter imaging together with biomechanical methods as a rapid and specific phenotyping approach to identify biologically significant and functionally relevant skeletal abnormalities using a minimal number of animals (n = 2).
Phenotyping of 100 knockout strains
To identify new genetic determinants of bone mass and strength, limbs from 100 knockout strains were analyzed for each structural and biomechanical variable. X-ray microradiography and micro-CT imaging identified 19 knockout strains in which at least one structural parameter was >2.0 standard deviations (SD) from the reference mean (Fig. 3, Figure S4 and Table S1). To ensure that significant abnormal phenotypes resulting from simultaneous but smaller variances in any of the six parameters were not overlooked, Mahalanobis distances were calculated as detailed in the methods and principal component analysis performed to identify multivariate outliers [24]–[27]. These studies identified 40 strains with outlier Mahalanobis distances (P<0.025), 24 of which had not been identified by analysis of individual x-ray microradiography or micro-CT values alone (Fig. 3 and Table S1). The MGP broad primary phenotype screen independently annotated 17 of these knockout strains with skeletal abnormalities (Table S1). Nine of the strains were also identified as outliers by x-ray microradiography, micro-CT or statistical methods whereas 8 did not display any abnormalities (Fig. 3).
The biomechanical significance of the 43 outlier phenotypes identified by imaging (n = 19, Faxitron and micro-CT) and statistical (n = 24, Mahalanobis analysis but excluding Faxitron and micro-CT) approaches, together with the 8 additional strains identified only in the MGP primary screen, was investigated (Fig. 3). Twelve of the 51 strains had at least one biomechanical parameter >2.0 SD from the reference mean (Fig. 3, Figure S5 and Table S1). However, 2 of these 12 strains had only minor abnormalities of bone morphology in the primary screen and were normal when investigated by x-ray microradiography, micro-CT and principal component analysis. Destructive 3-point bend testing of bones from the remaining 49 strains identified a further 7 with a single outlier biomechanical parameter but no other abnormality (Fig. 3 and Table S1).
In summary, the broad primary phenotype screen together with x-ray microradiography, micro-CT and statistical analysis identified knockout strains with at least one abnormal bone-related parameter. The addition of functional biomechanical testing demonstrated that 10 of these strains had major phenotypes affecting both the structure and strength of bone. Three of these carried heterozygous mutations (Asxl1 (ENSMUSG00000042548), Setdb1 (ENSMUSG00000015697) and Trim45 (ENSMUSG00000033233)) while the rest were homozygotes (Bbx (ENSMUSG00000022641), Cadm1 (ENSMUSG00000032076), Fam73b (ENSMUSG00000026858), Prpsap2 (ENSMUSG00000020528), Slc38a10 (ENSMUSG00000061306), Sparc and Spns2 (ENSMUSG00000040447)).
The primary phenotype database (http://www.sanger.ac.uk/mouseportal/) was interrogated for each of the strains identified with major skeletal phenotypes (Table 1).
Accuracy of rapid throughput skeletal phenotyping
The sensitivities, specificities and predictive values of each phenotyping method at a statistical threshold of >2.0 SD and 95% confidence limit were calculated to determine their ability to identify the 10 strains with major phenotypes. X-ray microradiography was the most accurate, identifying 8 of the 10 strains. Additional use of micro-CT and Mahalanobis analysis was required to identify the remaining 2 strains (Table S2). The MGP primary phenotype screen identified 5 of the 10 abnormal strains. Thus, addition of organ-specific imaging, statistical and biomechanical analyses to the primary phenotype screen resulted in increased sensitivity and specificity (Table S2), demonstrating the advantage of a complementary multi-parameter approach.
Correlations between imaging and biomechanical parameters were also determined to investigate relationships between bone structure and strength (Figure S6). Bone strength correlated strongly with BMC, cortical thickness and bone length but not with BV/TV, Tb.Th or Tb.N. Furthermore, there was no significant relationship between bone structural parameters determined by x-ray microradiography and micro-CT, thus demonstrating independence of the two techniques. Consistent with sensitivity, specificity and predictive value data (Table S2) these findings demonstrate that mineralization and cortical bone parameters determined by x-ray microradiography are excellent predictors of bone strength determined by 3-point bend testing, whereas trabecular bone structure determined by micro-CT is not. Importantly, trabecular bone parameters determined by micro-CT may represent good predictors of bone strength at sites of predominantly cancellous bone such as the vertebra or in response to age-related bone loss, although these possibilities were not investigated.
Functional classification of bone structure
Biomechanical analysis of wild-type mice and the 10 knockout strains with major phenotypes enabled four distinct biomechanical categories to be identified (Fig. 4). Normal bone had a stiffness of 30.2±4.1 N/mm (mean ± SD) with the capability to resist loading up to a yield load of 8.9±0.9 N (mean ± SD) and maximum load of 10.7±0.9 N (mean ± SD), prior to fracture at a load of 6.4±1.7 N (mean ± SD). These material properties of normal bone represent an optimised compromise between strength and flexibility that allows dissipation of 86.2% of energy prior to failure and limits the structural damage at fracture (Table S1).
Bones from strains with major phenotypes clustered into three abnormal biomechanical categories. Bones from Bbx, Cadm1 and Fam73b knockout mice were weak but flexible with reduced maximum load but the capability to bend and dissipate energy prior to fracture (Fig. 4B). Sparc, Prpsap2 and Slc38a10 bones were weak and brittle with reduced maximum load and lacking the capability to bend and dissipate energy prior to fracture (Fig. 4C). Asxl1, Trim45, Spns2 and Setdb1 bones were strong but brittle with an increased maximum load but were unable to bend and dissipate energy prior to fracture (Fig. 4D). A plot of the proportion of energy dissipated prior to fracture versus maximum load clearly separates these categories of bone strength (Fig. 4E). Further analysis demonstrated that Bbx, Cadm1, Fam73b Sparc, Prpsap2 and Slc38a10 bones have low BMC, whereas Asxl1, Trim45, Spns2 and Setdb1 bones have high BMC (Fig. 5). To determine whether these three categories of abnormal bone strength were related to a further morphological parameter known to have an important role for the biomechanical properties of long bones, we investigated their relationship with mid-diaphyseal diameter (Fig. 6). As expected [28], in WT mice mid-diaphyseal diameter correlated with fracture load and the proportion of energy dissipated prior to fracture. However, bones from the mutant strains identified with major skeletal phenotypes clustered into the same three abnormal categories, further demonstrating the validity of this functional classification.
Investigation of the biological activities of the 10 genes and their possible roles in bone indicates a broad diversity of function that was not clearly related to phenotype (Table 1), thus reinforcing the importance of an unbiased screening approach for identification of gene function in both homozygous and heterozygous mutants. In summary, a functional classification of four categories of bone structure was defined. Normal bone is strong and flexible with a normal mineral content, whereas abnormal bone is either (i) weak but flexible with low BMC, (ii) weak and brittle with low BMC or (iii) strong but brittle with high BMC.
Discussion
We have identified 10 genes with diverse and unrelated functions, the deletion of which resulted in major skeletal abnormalities. By adopting a multi-parameter phenotyping approach, we identified a new functional classification of bone structure based on its mineral content, strength and ductility that clarifies understanding of skeletal physiology and pathology, and which maps directly to human disease. As a result of evolutionary pressure, bone structure represents an optimal compromise between strength and flexibility that requires contributions from many diverse genes. Continuous bone remodeling enables the skeleton to adjust this compromise in response to changing physiological and environmental pressures [29], [30]. The current studies demonstrate that loss of function of individual genes can disrupt this optimal compromise resulting in skeletal phenotypes that cluster into three functionally distinct categories.
Postmenopausal osteoporosis is characterized by weak but flexible bone with low mineral content [31] and three of the identified knockout strains (Bbx, Cadm1, Fam73B) had phenotypes in this category. Bbx encodes a conserved transcription factor that contains a SOX-TCF HMG-box [32]–[36]. Family members include SRY-related Sox genes that are implicated in skeletal dysplasias [37], [38] and the TCF/LEF transcription factors that mediate Wnt/β-catenin signaling [39], a key pathway implicated in osteoporosis and osteoarthritis [40]–[43]. Cadm1 encodes a trans-membrane glycoprotein adhesion molecule of the immunoglobulin superfamily [44] for which a number of disparate functions have been reported including; tumor suppression [45], synapse development [46], behavioral regulation [47], T cell adhesion [48], mast cell interactions [49], and spermatogenesis [50]. However, no function in the skeleton has been reported. Fam73B encodes a conserved membrane protein of unknown function. These findings indicate that deletion of genes encoding proteins with diverse and unrelated functions can result in similar defects of bone strength and mineralization.
Disorders of bone matrix as typified by osteogenesis imperfecta [51] are characterized by bone that is weak and brittle with low BMC, and three of the strains (Prpsap2, Slc38a10, Sparc) displayed this phenotype. Prpsap2 encodes the non-catalytic inhibitory subunit of phosphoribosylpyrophosphate synthetase [52], and is required for synthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, and the coenzyme nicotinamide adenine dinucleotide [53]. Its function in the skeleton is unknown, although a recent study proposed PRPSAP2 as a candidate oncogene in osteosarcoma tumorigenesis [54]. Slc38a10 encodes a proposed sodium-coupled neutral amino acid membrane transporter [55] that may act as a cell volume regulator in mesenchyme [56]. The severe growth defect in Slc38a10 knockout mice suggests a critical function in chondrocytes, which mediate linear growth by cell volume expansion during hypertrophic differentiation [57]. Furthermore, related transporters have already been implicated in human skeletal disease. SLC35D1 is critical for chondroitin sulphate synthesis and mutations cause Schneckenbecken skeletal dysplasia [58]. Mutations in SLC26A2, cause four distinct chondrodysplasia syndromes [59] and emphasize the key role of these transporters in endochondral ossification. Sparc encodes the well-described extracellular matrix glycoprotein osteonectin and its deletion resulted in the characteristic and expected phenotype [22], [23] of weak and brittle bone with low BMC. These findings highlight the importance of enzymes, transporters and structural proteins to the functional integrity of bone matrix.
Diseases of high bone mass are rare and include sclerosteosis due to deletion of SOST [60] and autosomal dominant high bone mass due to gain-of-function mutations in LRP5 [61]. They are characterized by bone that is strong but brittle with high BMC, and four of the knockout strains (Asxl1, Setdb1, Spns2, Trim45) displayed such a phenotype. Asxl1 encodes a polycomb protein that interacts with heterochromatin protein-1 [62] and is required for regulation of Hox genes during axial patterning [63], suggesting a role in skeletal development [64], [65]. Indeed, ASXL1 heterozygous nonsense mutations were recently described to cause Bohring-Opitz syndrome [66], a developmental disorder characterized by mental retardation, impaired intrauterine growth, trigonocephaly and wrist and metacarpophalangeal joint abnormalities. Although the disease mechanism is unknown, craniofacial defects identified in homozygous Asxl1 knockout mice suggest that mutations in Bohring-Opitz syndrome result in a mutant protein with dominant-negative activity. Setdb1 encodes a histone H3 methyltransferase that regulates gene silencing [67], [68]. Although found to be expressed in cartilage but not bone in the primary phenotype screen, other studies demonstrated Setdb1 expression in osteoblasts and suggested a role in lineage commitment and differentiation [69], [70]. Spns2 encodes a sphingosine 1-phosphate (S1P) transporter [71] that is essential for S1P secretion. S1P binds to the G-protein coupled receptors, S1PR1 and S1PR2, and regulates osteoclast [72], [73] and osteoblast [74] precursor cell recruitment and migration. Thus, control of S1P secretion by Spns2 represents a novel mechanism that couples bone resorption and formation [75]. Trim45 is a member of the tripartite protein family, many of which act as ubiquitin or SUMO E3 ligases [76]–[78]. Although restricted to brain and testis in the primary phenotype screen, human studies demonstrate that Trim45 is more widely expressed [79]. Little is known about its function, although one study indicates Trim45 interacts with AP-1 and inhibits activity of the MAP kinase pathway [79]. The physiological significance of these findings and the role of Trim45 in the skeleton are unknown, although AP-1 proteins are key regulators of osteoblast and osteoclast differentiation and function [80], [81]. These findings emphasize the importance of lineage commitment, control of cell differentiation and coupling of both osteoblasts and osteoclasts in high bone mass disorders.
In the context of osteoporosis, our identification of many new genes that determine bone strength, and which otherwise could not be predicted, is consistent with studies indicating that diverse genetic polymorphisms result in small effects on phenotype [11], [12], [82], [83]. Accordingly, and in line with current understanding that only 3% of the heritability of BMD is accounted for by known genetic variation [12], none of the genes identified in this study have been recognized in osteoporosis genome-wide association studies [84]. We hypothesize, therefore, that unbiased multi-parameter and functional phenotyping of knockout mice has the power to identify many of the major genes that determine bone strength. Ultimately, this approach is likely to identify several genes from a single signaling pathway with an important role in the control of bone mass and strength. This has the advantage of independently confirming critical pathways and the potential to identify several alternative therapeutic targets. Importantly, however, the approach has limitations. The study of knockout mice can only identify phenotypes that result from gene deletion but cannot identify genes that only cause abnormalities when they harbor gain-of-function or dominant-negative mutations. Furthermore, the strategy does not include challenges such as ageing that may reveal additional phenotypes. However, if such provocative challenges were to be incorporated into screening approaches they would inevitably increase costs and limit throughput.
Our findings resulted from development and refinement of a rapid-throughput phenotyping algorithm to identify knockout mice with major abnormalities of bone structure and strength (Figure S3). The methods require bones from only two knockout mice, which first undergo digital point projection x-ray microradiography and micro-CT determination of six parameters of bone structure. Mahalanobis distance calculations and principal component analysis is performed and strains with at least one structural parameter >2.0 SD from the reference mean plus those with outlier Mahalanobis distances (95% confidence limit) are selected for biomechanical studies. Bones from selected strains undergo destruction 3-point bend testing to determine six measures of bone strength. Application of this unbiased approach to 100 consecutive knockout strains from the MGP pipeline identified 10% with major phenotypes affecting bone strength. Subsequent consideration of the results of primary phenotype screening and biological plausibility (Figure S3) allowed selection of mice to be refined.
Inherent in this approach is the capability to alter the statistical stringency threshold of analyses such that the number of strains for subsequent functional studies can be adjusted according to phenotype severity. For example, if the threshold for structural parameters is increased from 2.0 to 3.0 SD, then 9 outlier strains (rather than 19) are identified. Furthermore, if the confidence limit for Mahalanobis distance is increased from 95 to 99.7% then 21 multivariate outliers (rather than 40) would be identified. Of note, Trim45, which was recognized as an outlier only by Mahalanobis analysis, would still be identified if the confidence limit were to be increased to 99.7%, thus emphasizing the importance of a robust statistical method to ensure that all functional outliers are captured. Biomechanical analysis following application of these more stringent thresholds would detect 8 outlier strains including Trim45 and Sparc, resulting in the identification of 7 novel determinants of bone mass and strength rather than 9. The intrinsic flexibility of such a bespoke approach facilitates its application to other biological systems or polygenic diseases.
Materials and Methods
Ethics statement
All mouse studies were undertaken by Wellcome Trust Sanger Institute Mouse Genetics Project as part of the International Knockout Mouse Consortium and licensed by the UK Home Office in accordance with the Animals (Scientific Procedures) Act 1986 and the recommendations of the Weatherall report.
Primary phenotype screen
All mice generated by the MGP undergo a broad primary phenotype screen (http://www.sanger.ac.uk/mouseportal/) that includes measurement of body length, x-ray skeletal survey, DEXA analysis of bone mineral density and biochemical measures of mineral metabolism performed between 14–16 weeks of age, and determination of the normal tissue expression pattern of the targeted gene in 6–12 week old mice. Following primary phenotyping, lower limbs were fixed in 70% ethanol.
LacZ reporter gene expression
The pattern of LacZ reporter gene expression was determined in whole mount tissue preparations from heterozygous knockout mice between 6 and 12 weeks of age. Under terminal anaesthesia, mice were perfused with fresh cold 4% paraformaldehyde (PFA). Tissues were fixed for a further 30 min in 4% PFA, rinsed in phosphate buffered saline and stained with 0.1% X-gal for 48 hours at 4°C. Samples were subsequently fixed overnight in 4% PFA at 4°C, cleared with 50% glycerol and transferred to 70% glycerol. Specific and non-specific staining was determined in 41 tissues (Figure S2). A panel of 27 standardized images were recorded if expression was widespread (www.sanger.ac.uk/mouseportal).
Faxitron point projection digital x-ray microradiography
Bones from 16 week-old mice were fixed in 70% ethanol. Soft tissue was removed from the fixed bones and digital X-ray images were recorded at 10 µm pixel resolution using a Faxitron MX20 variable kV point projection x-ray source and digital image system (Qados, Cross Technologies plc, Sandhurst, Berkshire, UK) operating at 26 kV and 5× magnification [85]. Magnifications were calibrated by imaging a digital micrometer. Bone mineral content, bone length and cortical bone thickness were determined with coefficients of variation (CV) of 1.7%, 2.0% and 5.1%, respectively. The relative mineral content of calcified tissues was determined by comparison with standards included in each image frame, which comprised: a 1 mm steel plate; a 1 mm diameter spectrographically pure aluminum wire; and a 1 mm diameter polyester fiber. 2368×2340 16 bit DICOM images were converted to 8 bit Tiff images using ImageJ and the histogram stretched from the polyester (grey level 0) to steel (grey level 255) standards. Increasing gradations of mineralization density were represented in 16 equal intervals by a pseudocolor scheme. Cortical bone thickness was determined in at least 10 locations at the mid-femoral diaphysis. Bone length was determined using ImageJ 1.41 software (http://rsb.info.nih.gov/ij/).
Micro-CT
Tibias were analyzed by micro-CT (Skyscan 1172a, Skyscan, Belgium) at 50 kV and 200 µA using a 0.5 mm aluminum filter and a detection pixel size of 4.3 µm2 [85]. Images were captured every 0.7°, with 2× averaging, through 180° rotation of each bone and reconstructed using Skyscan NRecon software. A volume of 1 mm3 of trabecular bone was selected as the region of interest, 0.2 mm from the growth plate. Trabecular bone volume as proportion of tissue volume (BV/TV, %, CV 18.4%), trabecular thickness (Tb.Th, mm, CV 11.1%) and trabecular number (Tb. N, mm−1, CV 17.4%) were analyzed [86] using Skyscan CT analysis software.
Destructive 3-point bend testing
Bones were stored and tested in 70% ethanol. Destructive 3-point bend tests were performed on an Instron 5543 materials testing load frame (Instron Limited, High Wycombe, Buckinghamshire, UK) using custom built mounts with rounded supports that minimize cutting and shear loads [85]. Bones were positioned horizontally and centered on custom supports with the anterior surface upward. Load was applied vertically to the mid-shaft with a constant rate of displacement of 0.03 mm/second until fracture. A span of 12 mm was used. Load-displacement curves were plotted and yield load, maximum load and fracture load determined. Stiffness, the slope of the linear (elastic) part of the load-displacement curve, was calculated by the “least squares” method. Work energy was calculated from the area under curve at both maximum load and fracture. Elastic stored energy at maximum load was determined by calculating the area of a right angled triangle with the vertex at the point of maximum load and hypotenuse with a slope equal to that of the linear phase of the load-displacement curve. Elastic stored energy at fracture was similarly calculated but with the vertex of the triangle at the point of fracture (Fig. 2). Energy dissipated at maximum load or fracture was calculated by subtracting the elastic stored energy from the work energy at maximum load or fracture. CVs for each parameter were as follows: yield load (9.8%), maximum load (8.5%), fracture load (26.6%), stiffness (13.6%), the ratio of energy dissipated at maximum load to elastic stored energy at maximum load (25.1%), and the ratio of energy dissipated prior to fracture to elastic stored energy at fracture (11.0%).
Calculation of Mahalanobis distances
C57BL/6 reference ranges were generated for all Faxitron and micro-CT measures. Outliers in multivariate data were identified using robust Mahalanobis distances [27], which measure how far each observation is from the center of a data cluster, taking into account the shape of the cluster [87]. Robust Mahalanobis distances
were calculated for the vector of multivariate observations as described [27]. Here is a robust (i.e. relatively unaffected by outliers) estimate of the mean vector and is a robust estimate of the covariance matrix of the data set . Under the assumption of multivariate normality, if mouse i is from the same population as the rest of the data then has a chi-squared distribution with p degrees of freedom (where p is the number of variables). Robust estimates of the mean and covariance matrix are used so that potential outliers are not masked. The masking effect, by which outliers do not necessarily have a large Mahalanobis distance, can be caused by a small cluster of outliers that attract the mean and inflate the covariance in its direction. By replacing the sample mean and covariance with a robust estimate, the influence of these outliers is removed and the Mahalanobis distance is able to expose all outliers. Robust estimates of the mean and covariance matrix were calculated using the minimum volume ellipsoid method [87]. Given observations and variables, the minimum volume ellipsoid method seeks an ellipsoid containing
points of minimum volume. All analysis was conducted in the statistical computing package R (http://www.R-project.org).
Principal Component Analysis
Principal component analysis was used as a method to visualize multivariate data and reveal outliers by describing variation in a set of correlated variables in terms of a new set of uncorrelated variables. These new variables or principal components are linear combinations of the original variables derived in decreasing order of importance so that the first component accounts for the most variation of all possible linear combinations. The second component is then selected so that it accounts for as much of the remaining variance as possible (subject to it being uncorrelated with the first component), and so on [24]. Since the first few principal components often contain most of the variation in the data set they can be used as a lower-dimensional summary of the original variables.
Statistics
Normally distributed data were analyzed by Student's t test, or ANOVA followed by Tukey's multiple comparison post-hoc test. Relationships between bone structure and biomechanical measures, and between individual bone structure parameters, were determined by Pearson correlation. P values<0.05 were considered significant. Frequency distributions of bone mineral densities obtained by Faxitron were compared using the Kolmogorov-Smirnov test, in which P values for the D statistic in 1024 pixel data sets were D = >6.01 P<0.05, D = >7.20 P<0.01, and D = >8.62 P<0.001.
Supporting Information
Zdroje
1. FarooqiIS, JebbSA, LangmackG, LawrenceE, CheethamCH, et al. (1999) Effects of recombinant leptin therapy in a child with congenital leptin deficiency. N Engl J Med 341: 879–884.
2. MontagueCT, FarooqiIS, WhiteheadJP, SoosMA, RauH, et al. (1997) Congenital leptin deficiency is associated with severe early-onset obesity in humans. Nature 387: 903–908.
3. PearsonER, FlechtnerI, NjolstadPR, MaleckiMT, FlanaganSE, et al. (2006) Switching from insulin to oral sulfonylureas in patients with diabetes due to Kir6.2 mutations. N Engl J Med 355: 467–477.
4. YamagataK, OdaN, KaisakiPJ, MenzelS, FurutaH, et al. (1996) Mutations in the hepatocyte nuclear factor-1alpha gene in maturity-onset diabetes of the young (MODY3). Nature 384: 455–458.
5. BrunkowME, GardnerJC, Van NessJ, PaeperBW, KovacevichBR, et al. (2001) Bone dysplasia sclerosteosis results from loss of the SOST gene product, a novel cystine knot-containing protein. Am J Hum Genet 68: 577–589.
6. ChandrasekharappaSC, GuruSC, ManickamP, OlufemiSE, CollinsFS, et al. (1997) Positional cloning of the gene for multiple endocrine neoplasia-type 1. Science 276: 404–407.
7. GongY, SleeRB, FukaiN, RawadiG, Roman-RomanS, et al. (2001) LDL receptor-related protein 5 (LRP5) affects bone accrual and eye development. Cell 107: 513–523.
8. CollinsFS, RossantJ, WurstW (2007) A mouse for all reasons. Cell 128: 9–13.
9. SkarnesWC, RosenB, WestAP, KoutsourakisM, BushellW, et al. (2011) A conditional knockout resource for the genome-wide study of mouse gene function. Nature 474: 337–342.
10. JohnellO, KanisJA (2006) An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int 17: 1726–1733.
11. RalstonSH, UitterlindenAG (2010) Genetics of osteoporosis. Endocr Rev 31: 629–662.
12. RivadeneiraF, StyrkarsdottirU, EstradaK, HalldorssonBV, HsuYH, et al. (2009) Twenty bone-mineral-density loci identified by large-scale meta-analysis of genome-wide association studies. Nat Genet 41: 1199–1206.
13. CirulliET, GoldsteinDB (2010) Uncovering the roles of rare variants in common disease through whole-genome sequencing. Nat Rev Genet 11: 415–425.
14. BlackDM, ThompsonDE, BauerDC, EnsrudK, MuslinerT, et al. (2000) Fracture risk reduction with alendronate in women with osteoporosis: the Fracture Intervention Trial. FIT Research Group. J Clin Endocrinol Metab 85: 4118–4124.
15. MackeyDC, BlackDM, BauerDC, McCloskeyEV, EastellR, et al. (2011) Effects of antiresorptive treatment on non-vertebral fracture outcomes. J Bone Miner Res 26: 2411–2418.
16. BoydenLM, MaoJ, BelskyJ, MitznerL, FarhiA, et al. (2002) High bone density due to a mutation in LDL-receptor-related protein 5. N Engl J Med 346: 1513–1521.
17. RachnerTD, KhoslaS, HofbauerLC (2011) Osteoporosis: now and the future. Lancet 377: 1276–1287.
18. BrownSD, ChambonP, de AngelisMH (2005) EMPReSS: standardized phenotype screens for functional annotation of the mouse genome. Nat Genet 37: 1155.
19. JusticeMJ (2008) Removing the cloak of invisibility: phenotyping the mouse. Dis Model Mech 1: 109–112.
20. KarpNA, BakerLA, GerdinAK, AdamsNC, Ramirez-SolisR, et al. (2010) Optimising experimental design for high-throughput phenotyping in mice: a case study. Mamm Genome 21: 467–476.
21. Hardisty-HughesRE, ParkerA, BrownSD (2010) A hearing and vestibular phenotyping pipeline to identify mouse mutants with hearing impairment. Nat Protoc 5: 177–190.
22. BoskeyAL, MooreDJ, AmlingM, CanalisE, DelanyAM (2003) Infrared analysis of the mineral and matrix in bones of osteonectin-null mice and their wildtype controls. J Bone Miner Res 18: 1005–1011.
23. DelanyAM, AmlingM, PriemelM, HoweC, BaronR, et al. (2000) Osteopenia and decreased bone formation in osteonectin-deficient mice. J Clin Invest 105: 915–923.
24. Everitt B, editor. An R and S-PLUS Companion to Multivariate Analysis. London: Springer.
25. MitchellAFS, KrzanowskiWJ (1985) The Mahalanobis distance and elliptic distributions. Biometrika 72: 464–467.
26. RousseeuwPJ (1985) Multivariate estimation with high breakdown point;. W. Grossmann GP, I. Vinceze, and W. Wertz, Eds., editor The Netherlands: Reidel. 283–297.
27. RousseeuwPJ, Van ZomerenBC (1990) Unmasking Multivariate Outliers and Leverage Points. Journal of the American Statistical Association 85: 633–639.
28. AmmannP, RizzoliR (2003) Bone strength and its determinants. Osteoporos Int 14 Suppl 3: S13–18.
29. BoyleWJ, SimonetWS, LaceyDL (2003) Osteoclast differentiation and activation. Nature 423: 337–342.
30. HaradaS, RodanGA (2003) Control of osteoblast function and regulation of bone mass. Nature 423: 349–355.
31. DickensonRP, HuttonWC, StottJR (1981) The mechanical properties of bone in osteoporosis. J Bone Joint Surg Br 63-B: 233–238.
32. Marchler-BauerA, AndersonJB, ChitsazF, DerbyshireMK, DeWeese-ScottC, et al. (2009) CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res 37: D205–210.
33. DaileyL, BasilicoC (2001) Coevolution of HMG domains and homeodomains and the generation of transcriptional regulation by Sox/POU complexes. J Cell Physiol 186: 315–328.
34. LaudetV, StehelinD, CleversH (1993) Ancestry and diversity of the HMG box superfamily. Nucleic Acids Res 21: 2493–2501.
35. LoveJJ, LiX, CaseDA, GieseK, GrosschedlR, et al. (1995) Structural basis for DNA bending by the architectural transcription factor LEF-1. Nature 376: 791–795.
36. PohlerJR, NormanDG, BramhamJ, BianchiME, LilleyDM (1998) HMG box proteins bind to four-way DNA junctions in their open conformation. EMBO J 17: 817–826.
37. BernardP, TangP, LiuS, DewingP, HarleyVR, et al. (2003) Dimerization of SOX9 is required for chondrogenesis, but not for sex determination. Hum Mol Genet 12: 1755–1765.
38. MundlosS, OlsenBR (1997) Heritable diseases of the skeleton. Part I: Molecular insights into skeletal development-transcription factors and signaling pathways. FASEB J 11: 125–132.
39. MacDonaldBT, TamaiK, HeX (2009) Wnt/beta-catenin signaling: components, mechanisms, and diseases. Dev Cell 17: 9–26.
40. BalemansW, Van HulW (2007) The genetics of low-density lipoprotein receptor-related protein 5 in bone: a story of extremes. Endocrinology 148: 2622–2629.
41. BaronR, RawadiG (2007) Wnt signaling and the regulation of bone mass. Curr Osteoporos Rep 5: 73–80.
42. ChanA, van BezooijenRL, LowikCW (2007) A new paradigm in the treatment of osteoporosis: Wnt pathway proteins and their antagonists. Curr Opin Investig Drugs 8: 293–298.
43. LuytenFP, TylzanowskiP, LoriesRJ (2009) Wnt signaling and osteoarthritis. Bone 44: 522–527.
44. GiangrecoA, JensenKB, TakaiY, MiyoshiJ, WattFM (2009) Necl2 regulates epidermal adhesion and wound repair. Development 136: 3505–3514.
45. KuramochiM, FukuharaH, NobukuniT, KanbeT, MaruyamaT, et al. (2001) TSLC1 is a tumor-suppressor gene in human non-small-cell lung cancer. Nat Genet 27: 427–430.
46. BiedererT, SaraY, MozhayevaM, AtasoyD, LiuX, et al. (2002) SynCAM, a synaptic adhesion molecule that drives synapse assembly. Science 297: 1525–1531.
47. TakayanagiY, FujitaE, YuZ, YamagataT, MomoiMY, et al. (2010) Impairment of social and emotional behaviors in Cadm1-knockout mice. Biochem Biophys Res Commun 396: 703–708.
48. GiangrecoA, HosteE, TakaiY, RosewellI, WattFM (2012) Epidermal Cadm1 expression promotes autoimmune alopecia via enhanced T cell adhesion and cytotoxicity. J Immunol 188: 1514–1522.
49. ItoA, HagiyamaM, OonumaJ (2008) Nerve-mast cell and smooth muscle-mast cell interaction mediated by cell adhesion molecule-1, CADM1. J Smooth Muscle Res 44: 83–93.
50. FujitaE, KourokuY, OzekiS, TanabeY, ToyamaY, et al. (2006) Oligo-astheno-teratozoospermia in mice lacking RA175/TSLC1/SynCAM/IGSF4A, a cell adhesion molecule in the immunoglobulin superfamily. Mol Cell Biol 26: 718–726.
51. RauchF, GlorieuxFH (2004) Osteogenesis imperfecta. Lancet 363: 1377–1385.
52. TatibanaM, KitaK, TairaM, IshijimaS, SonodaT, et al. (1995) Mammalian phosphoribosyl-pyrophosphate synthetase. Adv Enzyme Regul 35: 229–249.
53. KatashimaR, IwahanaH, FujimuraM, YamaokaT, ItakuraM (1998) Assignment of the human phosphoribosylpyrophosphate synthetase-associated protein 41 gene (PRPSAP2) to 17p11.2–p12. Genomics 54: 180–181.
54. BothJ, WuT, BrasJ, SchaapGR, BaasF, et al. (2012) Identification of novel candidate oncogenes in chromosome region 17p11.2–p12 in human osteosarcoma. PLoS One 7: e30907.
55. SundbergBE, WaagE, JacobssonJA, StephanssonO, RumaksJ, et al. (2008) The evolutionary history and tissue mapping of amino acid transporters belonging to solute carrier families SLC32, SLC36, and SLC38. J Mol Neurosci 35: 179–193.
56. Franchi-GazzolaR, Dall'AstaV, SalaR, VisigalliR, BevilacquaE, et al. (2006) The role of the neutral amino acid transporter SNAT2 in cell volume regulation. Acta Physiol (Oxf) 187: 273–283.
57. KronenbergHM (2003) Developmental regulation of the growth plate. Nature 423: 332–336.
58. HiraokaS, FuruichiT, NishimuraG, ShibataS, YanagishitaM, et al. (2007) Nucleotide-sugar transporter SLC35D1 is critical to chondroitin sulfate synthesis in cartilage and skeletal development in mouse and human. Nat Med 13: 1363–1367.
59. Superti-FurgaA, HastbackaJ, RossiA, van der HartenJJ, WilcoxWR, et al. (1996) A family of chondrodysplasias caused by mutations in the diastrophic dysplasia sulfate transporter gene and associated with impaired sulfation of proteoglycans. Ann N Y Acad Sci 785: 195–201.
60. LiX, OminskyMS, NiuQT, SunN, DaughertyB, et al. (2008) Targeted deletion of the sclerostin gene in mice results in increased bone formation and bone strength. J Bone Miner Res 23: 860–869.
61. AkhterMP, WellsDJ, ShortSJ, CullenDM, JohnsonML, et al. (2004) Bone biomechanical properties in LRP5 mutant mice. Bone 35: 162–169.
62. LeeSW, ChoYS, NaJM, ParkUH, KangM, et al. (2010) ASXL1 represses retinoic acid receptor-mediated transcription through associating with HP1 and LSD1. J Biol Chem 285: 18–29.
63. FisherCL, LeeI, BloyerS, BozzaS, ChevalierJ, et al. (2010) Additional sex combs-like 1 belongs to the enhancer of trithorax and polycomb group and genetically interacts with Cbx2 in mice. Dev Biol 337: 9–15.
64. WellikDM (2007) Hox patterning of the vertebrate axial skeleton. Dev Dyn 236: 2454–2463.
65. WilliamsJA, KondoN, OkabeT, TakeshitaN, PilchakDM, et al. (2009) Retinoic acid receptors are required for skeletal growth, matrix homeostasis and growth plate function in postnatal mouse. Dev Biol 328: 315–327.
66. HoischenA, van BonBW, Rodriguez-SantiagoB, GilissenC, VissersLE, et al. (2011) De novo nonsense mutations in ASXL1 cause Bohring-Opitz syndrome. Nat Genet 43: 729–731.
67. SchultzDC, AyyanathanK, NegorevD, MaulGG, RauscherFJ3rd (2002) SETDB1: a novel KAP-1-associated histone H3, lysine 9-specific methyltransferase that contributes to HP1-mediated silencing of euchromatic genes by KRAB zinc-finger proteins. Genes Dev 16: 919–932.
68. YangL, XiaL, WuDY, WangH, ChanskyHA, et al. (2002) Molecular cloning of ESET, a novel histone H3-specific methyltransferase that interacts with ERG transcription factor. Oncogene 21: 148–152.
69. NifujiA, IdenoH, OhyamaY, TakanabeR, ArakiR, et al. (2010) Nemo-like kinase (NLK) expression in osteoblastic cells and suppression of osteoblastic differentiation. Exp Cell Res 316: 1127–1136.
70. TakadaI, MiharaM, SuzawaM, OhtakeF, KobayashiS, et al. (2007) A histone lysine methyltransferase activated by non-canonical Wnt signalling suppresses PPAR-gamma transactivation. Nat Cell Biol 9: 1273–1285.
71. HisanoY, KobayashiN, KawaharaA, YamaguchiA, NishiT (2011) The sphingosine 1-phosphate transporter, SPNS2, functions as a transporter of the phosphorylated form of the immunomodulating agent FTY720. J Biol Chem 286: 1758–1766.
72. IshiiM, EgenJG, KlauschenF, Meier-SchellersheimM, SaekiY, et al. (2009) Sphingosine-1-phosphate mobilizes osteoclast precursors and regulates bone homeostasis. Nature 458: 524–528.
73. IshiiT, ShimazuY, NishiyamaI, KikutaJ, IshiiM (2011) The role of sphingosine 1-phosphate in migration of osteoclast precursors; an application of intravital two-photon microscopy. Mol Cells 31: 399–403.
74. PedersonL, RuanM, WestendorfJJ, KhoslaS, OurslerMJ (2008) Regulation of bone formation by osteoclasts involves Wnt/BMP signaling and the chemokine sphingosine-1-phosphate. Proc Natl Acad Sci U S A 105: 20764–20769.
75. RyuJ, KimHJ, ChangEJ, HuangH, BannoY, et al. (2006) Sphingosine 1-phosphate as a regulator of osteoclast differentiation and osteoclast-osteoblast coupling. EMBO J 25: 5840–5851.
76. ChuY, YangX (2011) SUMO E3 ligase activity of TRIM proteins. Oncogene 30: 1108–1116.
77. NapolitanoLM, JaffrayEG, HayRT, MeroniG (2011) Functional interactions between ubiquitin E2 enzymes and TRIM proteins. Biochem J 434: 309–319.
78. OzatoK, ShinDM, ChangTH, MorseHC3rd (2008) TRIM family proteins and their emerging roles in innate immunity. Nat Rev Immunol 8: 849–860.
79. WangY, LiY, QiX, YuanW, AiJ, et al. (2004) TRIM45, a novel human RBCC/TRIM protein, inhibits transcriptional activities of ElK-1 and AP-1. Biochem Biophys Res Commun 323: 9–16.
80. KomoriT (2006) Regulation of osteoblast differentiation by transcription factors. J Cell Biochem 99: 1233–1239.
81. WagnerEF (2010) Bone development and inflammatory disease is regulated by AP-1 (Fos/Jun). Ann Rheum Dis 69 Suppl 1: i86–88.
82. RichardsJB, RivadeneiraF, InouyeM, PastinenTM, SoranzoN, et al. (2008) Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study. Lancet 371: 1505–1512.
83. StyrkarsdottirU, HalldorssonBV, GretarsdottirS, GudbjartssonDF, WaltersGB, et al. (2008) Multiple genetic loci for bone mineral density and fractures. N Engl J Med 358: 2355–2365.
84. EstradaK, StyrkarsdottirU, EvangelouE, HsuYH, DuncanEL, et al. (2012) Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat Genet 44: 491–501.
85. BassettJH, BoydeA, HowellPG, BassettRH, GallifordTM, et al. (2010) Optimal bone strength and mineralization requires the type 2 iodothyronine deiodinase in osteoblasts. Proc Natl Acad Sci U S A 107: 7604–7609.
86. ParfittAM, DreznerMK, GlorieuxFH, KanisJA, MallucheH, et al. (1987) Bone histomorphometry: standardization of nomenclature, symbols, and units. Report of the ASBMR Histomorphometry Nomenclature Committee. J Bone Miner Res 2: 595–610.
87. Venables WN, Ripley BD (2002) Modern Applied Statistics with S New York, NY, USA: Springer.
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2012 Číslo 8
Nejčtenější v tomto čísle
- Dissecting the Gene Network of Dietary Restriction to Identify Evolutionarily Conserved Pathways and New Functional Genes
- It's All in the Timing: Too Much E2F Is a Bad Thing
- Variation of Contributes to Dog Breed Skull Diversity
- The PARN Deadenylase Targets a Discrete Set of mRNAs for Decay and Regulates Cell Motility in Mouse Myoblasts