| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Genetics (L.M.H., M.C.M., L.A.C., J.R.) and Southwest National Primate Research Center (M.C.M., J.R.), Southwest Foundation for Biomedical Research, San Antonio, Texas 78245-0549; Southwest Fisheries Science Center (P.A.M.), La Jolla, California 92037; and Ernest Gallo Clinic and Research Center (G.J.), Emeryville, California 94608
Address all correspondence and requests for reprints to: Dr. Lorena M. Havill, Department of Genetics, Southwest Foundation for Biomedical Research, P.O. Box 760549, San Antonio, Texas 78245-0549. E-mail: lhavill{at}darwin.sfbr.org.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Although environmental factors (e.g. age, sex, diet, and exercise) are important contributors to individual variation in BMD, statistical genetic analyses in humans and nonhuman animal models clearly indicate that genes also play a substantial role (7, 8), accounting for 4090% of BMD variance (9). Studies provide increasing evidence that gene effects vary by skeletal site (9, 10) and that multiple genes are involved (7, 8).
Several whole-genome linkage screens in humans have resulted in significant evidence for quantitative trait loci (QTLs) affecting BMD (8), including one on chromosome 11q12-13 that has received considerable attention. This region harbors the gene for low-density lipoprotein receptor-related protein 5 (LRP5), mutations of which are responsible for an autosomal dominant high bone mass phenotype (11, 12) and an autosomal recessive low bone mass phenotype in the case of osteoporosis-pseudoglioma syndrome (13). Two other Mendelian disorders, autosomal recessive osteopetrosis (14) and autosomal dominant osteopetrosis (15), also map to this region.
Equivocal, but mounting, evidence indicates that the 11q12-13 region may also harbor a gene or genes that affect normal variation in BMD in the general population (16, 17, 18, 19). This effect may be due to LRP5 polymorphisms (20) or to other genes in the region or both. Further research is warranted to elucidate the relationship between 11q12-13 and normal variation in BMD.
Animal models can be used to facilitate and expedite the search for genes influencing normal variation in osteoporosis risk factors. Indeed, much of the current knowledge of the genetics of BMD and bone turnover stems from animal models such as mice, which, relative to humans, allow for enhanced experimental manipulation and environmental control. However, extrapolation of results concerning genetic influences from mice to humans can be problematic because of the phylogenetic distance between the species and the consequent differences in skeletal maintenance and repair (21).
In contrast, the baboon, a nonhuman primate of the taxonomic family Cercopithecidae (Old World monkeys), shares physiological and developmental characteristics that make it particularly well-suited to studies of skeletal maintenance and turnover. This animal shares with humans a relatively long lifespan, bone loss with advancing age (22, 23, 24), and physiological and endocrine correlates of reproductive biology and senescence, including menopause (25, 26). More extensive characterization of the baboon model, including normal variation and age changes in aBMD in the baboon, can be found in Havill et al. (24). The similarities between baboons and humans inspire confidence that study results will be substantively relevant to humans. In addition, when QTLs are identified in the baboon, the existence of a genetic linkage map of the baboon genome consisting primarily of human microsatellite markers allows us to easily identify the corresponding chromosomal region in humans (27).
We conducted the study presented here to characterize the impact of genetic variation on normal phenotypic variation in forearm aBMD using the baboon. The specific aims of this study were to detect and quantify the effects of genes on variation in aBMD in the baboon forearm and to localize chromosomal regions that harbor loci that contribute to genetic effects. Our results reveal a QTL in a region corresponding to human 11q12-13 that affects normal variation in this baboon population. The existence of QTLs in the same location across these two primate species indicates that there is a gene (or genes) in this region that is important to normal variation in BMD in both human and nonhuman primates, thereby attesting to the importance of this region to primate skeletal biology and, ultimately, to osteoporosis risk in humans.
| Materials and Methods |
|---|
|
|
|---|
After sedation by means of im ketamine followed by iv rompun, atropine, acepromazine, and ketamine to effect relaxation, aBMD (grams per centimeter squared) was assessed at six different skeletal sites via DXA (Lunar DPX +6529, Lunar Corp., Madison, WI). All DXA measurements were obtained using the manufacturers software for adults, version 3.65. Animals were scanned in a standard supine position. The scan width varied for each animal. We used Auto Width (Intelligent Scan Mode) to automatically narrow the path of the scan after locating bone masses. Speed was set manually depending on the thickness of the animal. Medium Mode was used for most animals, but Slow Mode was used for larger ones (>26-cm chest thickness). Both modes involve a scan interval of 1.2 x 1.2 mm. Three forearm sites were measured: ultradistal radius, diaphyseal radius at a point 33% proximal to the styloid process (radius 1/3), and diaphyseal ulna at a point 33% proximal to the styloid process (ulna 1/3). Five percent of the monkeys were scanned twice. The resulting coefficients of variation for ultradistal radius, radius 1/3, and ulna 1/3 were 2.3, 2.5, and 2.5%, respectively.
Data on additional variables possibly related to aBMD were also collected for each animal at the time of DXA. These variables include body weight, height (crown-rump length), and spinal osteophytosis score [related to development of spinal arthritis and diffuse idiopathic spinal hyperostosis in these animals assessed subjectively from visual inspection of radiographs and assigned a score on a four-point ordinal scale (0 = no osteophytosis, 3 = most severe osteophytosis)].
Baboon genotyping
Microsatellite markers were amplified from baboon genomic DNA by PCR. The forward primer for each microsatellite carried a fluorescent label. Amplification reactions included 2 U TaKaRa (Shiga, Japan) Taq Polymerase, 1x TaKaRa buffer, 0.33 mM dNTP mix, 1 µM forward and reverse primers, and 50 ng baboon genomic DNA. MgCl2 concentrations were optimized for each primer pair. PCR was performed with an initial denaturing step of 5 min at 94 C followed by 35 cycles of 40-sec denaturation at 94 C, 30-sec annealing, and 30-sec extension at 72 C followed by a final 5-min extension at 72 C. Annealing temperatures were optimized for each primer pair. Genotypes of baboons were determined through gel electrophoresis of the fluorescently labeled PCR products in ABI 377 automated sequencers with Gene Scan software (PerkinElmer Life and Analytical Sciences, Boston, MA) and analyzed using Genotyper software (PerkinElmer Life and Analytical Sciences). Detailed information about marker properties has been published elsewhere (27) and can be found at http://www.snprc.org/baboon/genome/index.html.
Baboon whole-genome linkage map
Genotype data from many of the same animals for whom aBMD data were obtained for this study were used to develop a first generation genetic linkage map of the baboon genome (27). Genotype data for 325 human microsatellite loci [short tandem repeats (STRs)], and six novel baboon microsatellites were used in marker-to-marker linkage analyses, facilitated by the expert system program Multimap (30, 31), which implements routines of the computer program CRIMAP (30), to produce a meiotic recombination map covering all 20 baboon autosomes. The average maker density is 7.2 cM. The most recent version of this map, used for multipoint interval mapping of aBMD QTLs in this study, contains 275 of these STRs that have been placed in unique positions at 1000:1 odds and eight more that have been placed in unique positions on the map at 100:1 odds.
Direct comparison among homologous (orthologous) loci reveals large regions of synteny in which the human marker order is conserved [seven autosomes with no major rearrangement, 15 with one or more rearrangements; please see Rogers et al. (27) for further details]. Given the degree of similarity, we have chosen to refer to baboon chromosomes by the number usually assigned their human homolog throughout the manuscript (with the baboon chromosome number in parentheses) to facilitate comprehension of our results and their comparison with the human situation. For example, chromosome 6 (Pha 4) designates the syntenic grouping of microsatellite marker loci that map to human chromosome 6, a syntenic grouping that comprises the loci on chromosome 4 in the baboon (Papio hamadryas; Pha). Using similar logic, chromosome 7_21 (Pha 3) designates a baboon chromosome that represents (relative to the human condition) a fusion of what are two different syntenic groups in humans, i.e. chromosomes 7 and 21.
Baboon pedigree
The baboons are assigned to 11 extended pedigrees that include the following relative pairs: 412 parent-offspring, 381 siblings, 32 grandparent-grandchild, 72 avuncular, 5951 half-siblings, 1212 half-avuncular, 4 first cousins, 179 half-first cousins, 40 half-first cousins once removed, 2 half-siblings and first cousins, 54 half-siblings and half-first cousins, 13 half-siblings and half-avuncular, and six double half-avuncular. Figure 1
shows one of the smaller baboon pedigrees to illustrate the general complexity.
|
We used a variance decomposition approach implemented in SOLAR (Sequential Oligogenic Linkage Analysis Routines) (32) to test for evidence of linkage between QTLs for aBMD and the 283 autosomal STR loci. This method, described in detail elsewhere (32), entails specification of the genetic covariance between arbitrary relatives as a function of the identity-by-descent (IBD) relationships at a given marker locus. The covariance matrix for a pedigree is modeled as the sum of the additive genetic covariance attributable to the QTL, the additive genetic covariance due to the effects of loci other than the QTL, and the variance due to unmeasured environmental factors.
We estimated marker locus-specific IBD probabilities for the pedigrees using a pair-wise maximum likelihood-based procedure (32). To permit multipoint analysis for QTL mapping, we employed an extension to the method of Fulker et al. (33) to estimate IBD probabilities at 1-centimorgan (cM) intervals along each chromosome. An LOD score evaluation was performed every centimorgan along each chromosome.
We tested the hypothesis of linkage by comparing the likelihood of a restricted model in which variance due to the QTL was constrained to zero (no linkage) with that of a model in which it was estimated. The LOD score of classical linkage analysis was obtained as the quotient of the difference between the 2 ln likelihoods divided by ln10 (34).
Significance of the mean effect of each covariate [age, age2, sex, sex x age2, height (crown-rump length in baboons), and weight] was assessed using a likelihood ratio test. This test compares differences in the likelihoods of a restricted model, in which the mean effect of the covariate to be tested is constrained to zero, and a general model, in which all parameters are estimated, to the
2 distribution with 1 df. If the model including the covariate effect was marginally significantly better than the model without it (P < 0.10), then the covariate was included in the final model. Reported h2 is residual heritability (that part of the variance that is attributable to the additive effects of genes after covariate effects are removed).
To control for the overall false positive rate given the finite marker locus density in the baboon genome linkage map, we estimated the LOD score associated with a genome-wide P value of 0.05 by means of a method suggested by Feingold et al. (35) At LOD more than 2.69, the genome-wide P value was <0.05. Therefore, we employ LOD equalling 2.7 as the threshold for genome-wide significance (at
= 0.05).
| Results |
|---|
|
|
|---|
|
|
Figures 24![]()
![]()
present the genome-wide linkage results for ultradistal radius, radius 1/3, and ulna 1/3, respectively. The highest LOD score for each phenotype (3.11 for ultradistal radius, 3.00 for radius 1/3, and 1.95 for ulna 1/3) occurred on chromosome 11q (Pha 14) at points approximately 6061 cM from the pter-most marker (see Fig. 5
), with a peak in the area of markers D11S1329, D11S1338, and D11S916. For our baboon map, the 95% confidence interval includes an area 0.8340 below the peak LOD score. Focusing on the two peaks that reach statistical significance (those for the two radius sites), the 95% confidence interval for the ultradistal radius is a region extending from approximately 51100 cM from the pter-most marker. This 49-cM region includes the following markers [physical location in humans (megabase pairs from pter) in parentheses]: D11S4200 (34), D11S904 (26), D11S1349 (11), D11S1329 (10), D11S1338 (5), D11S916 (73), D11S1902 (74), D11S1352 (78), D11S2002 (79), and D11S1366 (96). For the radius 1/3, the 95% confidence interval is a 16-cM interval extending from approximately 4766 cM from the pter-most marker. This region includes all of the markers listed for the ultradistal radius 95% confidence interval, with the exception of D11S1366. The estimated effect sizes for the QTLs (h2Q) are 0.29 for the ultradistal radius, 0.30 for the radius 1/3, and 0.16 for the ulna 1/3. These values are 83, 71, and 89% of the entire additive genetic effect for the ultradistal radius, the radius 1/3, and the ulna 1/3, respectively. A post hoc pedigree-specific LOD score evaluation reveals that no single pedigree or sire line accounts for more than 38% of the total significant LOD scores for these traits; at least seven of the 11 pedigrees contribute positively to their LODs, with the highest LOD coming from a different pedigree in each case.
|
|
|
|
Although no other chromosomal regions yielded a LOD greater than 2.69 (our genome-wide significance level at P = 0.05), several regions showed LOD scores greater than 1.0 including chromosomes 2q (Pha 12) (ultradistal radius, LOD 1.11; radius 1/3, LOD 1.91), 4p (Pha 5p) (ultradistal radius, LOD 1.20; radius 1/3, LOD 1.65), 7_21p (Pha 3p) (ultradistal radius, LOD 1.40; radius 1/3, LOD 1.06), and 7_21q (Pha 3q) (radius 1/3, LOD 1.47; ulna 1/3, LOD 1.16).
| Discussion |
|---|
|
|
|---|
This first genome-wide linkage scan in a nonhuman primate for QTLs affecting aBMD provides cross-species replication, in the form of significant evidence for linkage, of a QTL on chromosome 11q (Pha 14), a region that has received considerable attention in humans. This region, specifically 11q12-13, harbors the LRP5 gene and the T-cell immune regulator 1 (TCIRG1) gene, mutations of which are associated with the Mendelian disorders of the skeleton mentioned in the introduction. Naturally, this has prompted study into whether or not these same genes, or others in the region, might also contribute significantly to normal population level variation in BMD. Indeed, there is mounting evidence that this may be the case. Although not all published investigations of linkage between 11q12-13 and normal variation in bone phenotypes have found such (36), several studies report suggestive to significant evidence for linkage or association with genetic markers in this region. These include a QTL affecting femoral neck and lumbar spine aBMD in a Caucasian and African-American sample (16, 19), suggestive evidence of linkage with low aBMD in an Irish sample (18), evidence for association with radiographically assessed BMD of hand phalanges in a Chuvasha population (17), and suggestive evidence for association with metacarpal cortical index in a Chuvasha sample (37). Investigations of specific genes in this area have only just begun. A recent report indicates that particular LRP5 polymorphisms are associated with lumbar spine bone mineral content, bone area, and vertebral bone growth, particularly in males (20). Koay et al. (38) identified three single-nucleotide polymorphisms in the LRP5 gene that show significant associations with aBMD in their Caucasian sample. Another study found evidence of linkage between femoral neck aBMD and TCIRG1 but no evidence of association with polymorphisms in this gene (39).
Recent investigations of the relationship between LRP5 and mechanosensation are interesting in light of our finding of an 11q12-13 QTL in the forearm, a load-bearing site in the baboon. We did not find a QTL in this location for BMD of the spine (our unpublished data), a site that is subject to less mechanical loading in this locomotor quadruped. LRP5 is a coreceptor for Wnt, and recent research indicates that the role of Wnt/LRP5 signaling in the regulation of bone mass is linked to activation of this pathway through mechanical loading (40). Intensive and ongoing studies of the molecular mechanism by which a G171V mutation (41) in the LRP5 gene (using transgenic mice) regulates bone mass indicate that the high bone mass in individuals with this mutation results from decreased osteoblast and osteocyte apoptosis, increased bone formation in response to mechanical loading due to a lower threshold of sensitivity to mechanosensation, and potentially to reduced osteoclastogenesis (40). If LRP5 is indeed the candidate gene responsible for our linkage signal, the mediation of bone mass by LRP5 through increased bone formation in areas of greater mechanical loading could explain our detection of this effect in the baboon forearm.
Our results demonstrate that there is a gene (or genes) in this area that is very important to normal variation in forearm aBMD in this population of baboons. The effect size of the 11q12-13 QTL we report accounts for a very large portion of the total genetic effect (83, 71, and 89% for the ultradistal radius, radius 1/3, and ulna 1/3, respectively). Moreover, because our post hoc pedigree-specific LOD analysis shows that the LOD score we obtained stems from multiple pedigrees in the sample, rather than coming primarily from one or two pedigrees, we conclude it is unlikely that the detected QTL represents a single mutation with a major effect segregating through the baboon population. This QTL is in the same location that has been found to influence normal variation in BMD in humans at other skeletal sites (femoral neck and spine). This apparent conservation of function to the same genomic region attests to the overall importance of this region to human skeletal biology and, indeed, to primate skeletal biology in general. The fact that a QTL affecting BMD in inbred mice also maps to 11q12-13 may mean that this region is important beyond the primate order (42).
We cannot, of course, at this point, attribute the QTL effect to any particular gene. Although we consider LRP5 and TCIRG1 to be likely candidate genes responsible for the linkage to 11q12-13, there are other potential candidates in this region. These include FOS-like antigen 1 (FOSL1), latent transforming growth factor ß-binding protein (LTBP3), and fibronectin like 2 (FNL2) protein. Mice without FOS develop osteopetrosis due to a differentiation block in the osteoclast lineage (43). FOSL1 showed the highest rescue activity in a study of retroviral gene transfer of FOS genes that resulted in rescue of differentiation blocks in vitro, and a transgene expressing FOSL1 rescued the osteopetrosis of the mice lacking FOS (44). LTBP3 is a part of the biologically inactive large latent complex in which TGF-ß, a known stimulator of bone formation (45), is secreted, and affects TGF-ß bioavailability (46). FNL2 has not itself been directly shown to affect bone tissue, but fibronectin is an osteoid scaffolding protein and a downstream direct transcriptional target of the Wnt signaling pathway (47) and is found at elevated levels in serum in a high bone mass kindred (41).
Ultimately, the value of this study lies in its contribution to a growing body of data indicating that intensive investigation of genes in the 11q12-13 region is warranted because a gene or genes in this area affect not only extreme disorders of BMD but also normal variation in this fracture risk-related phenotype. Only after we successfully identify the gene and polymorphisms in this region that affect BMD can we screen individuals for the functional polymorphisms. Such screening would provide a mechanism for earlier identification of those at greatest risk of osteoporosis, allowing for earlier and more efficacious treatment and intervention. Given its current state of genetic characterization, the baboon provides an animal model for generating and testing hypotheses involving human skeletal genetics that is more directly biologically relevant than rodent models. Also, relative to humans, the baboon is a more accessible and less expensive model for testing the effects of potential treatment and intervention strategies on bone maintenance and turnover.
| Acknowledgments |
|---|
| Footnotes |
|---|
First Published Online March 8, 2005
Abbreviations: aBMD, Areal bone mineral density; DXA, dual-energy x-ray absorptiometry; IBD, identity-by-descent; LRP5, lipoprotein receptor-related protein 5; QTL, quantitative trait locus; STR, short tandem repeat.
Received August 19, 2004.
Accepted March 2, 2005.
| References |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |