Email updates

Keep up to date with the latest news and content from BMC Genetics and BioMed Central.

Open Access Highly Accessed Research article

A genome-wide scan for quantitative trait loci affecting limb bone lengths and areal bone mineral density of the distal femur in a White Duroc × Erhualian F2 population

Huirong Mao1, Yuanmei Guo12, Guangcheng Yang1, Bin Yang1, Jun Ren1, Sanfeng Liu1, Huashui Ai1, Junwu Ma1, Bertram Brenig2 and Lusheng Huang1*

Author Affiliations

1 Key Laboratory for Animal Biotechnology of Jiangxi Province and the Ministry of Agriculture of China, Jiangxi Agricultural University, Nanchang 330045, PR China

2 Institute of Veterinary Medicine, Georg-August-University of Göttingen, Burckhardtweg 2, 37077 Göttingen, Germany

For all author emails, please log on.

BMC Genetics 2008, 9:63  doi:10.1186/1471-2156-9-63

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/9/63


Received:4 February 2008
Accepted:8 October 2008
Published:8 October 2008

© 2008 Mao et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Limb bone lengths and bone mineral density (BMD) have been used to assess the bone growth and the risk of bone fractures in pigs, respectively. It has been suggested that limb bone lengths and BMD are under genetic control. However, the knowledge about the genetic basis of the limb bone lengths and mineralisatinon is limited in pigs. The aim of this study was to identify quantitative trait loci (QTL) affecting limb bone lengths and BMD of the distal femur in a White Duroc × Erhualian resource population.

Results

Limb bone lengths and femoral bone mineral density (fBMD) were measured in a total of 1021 and 116 F2 animals, respectively. There were strong positive correlations among the lengths of limb bones and medium positive correlations between the lengths of limb bones and fBMD. A whole-genome scan involving 183 microsatellite markers across the pig genome revealed 35 QTL for the limb bone lengths and 2 for femoral BMD. The most significant QTL for the lengths of five limb bones were mapped on two chromosomes affecting all 5 limb bones traits. One was detected around 57 cM on pig chromosome (SSC) 7 with the largest F-value of more than 26 and 95% confidence intervals of less than 5 cM, providing a crucial start point to identify the causal genes for these traits. The Erhualian alleles were associated with longer limb bones. The other was located on SSCX with a peak at 50–53 cM, whereas alleles from the White Duroc breed increased the bone length. Many QTL identified are homologous to the human genomic regions containing QTL for bone-related traits and a list of interesting candidate genes.

Conclusion

This study detected the QTL for the lengths of scapula, ulna, humerus and tibia and fBMD in the pig for the first time. Moreover, several new QTL for the pig femoral length were found. As correlated traits, QTL for the lengths of five limb bones were mainly located in the same genomic regions. The most promising QTL for the lengths of five limb bones on SSC7 merits further investigation.

Background

Bone length and bone mineral density (BMD) are generally regarded as two important parameters to assess the bone growth in pigs [1]. Individuals with longer limb bones usually have taller body heights or heights at shoulder. It has been shown that the body height or height at shoulder is negatively correlated with backfat thickness, and the height at shoulder is an important item in determining the yield of ham, loin, picnic shoulder and shoulder butt [2,3]. Visual selection for moderate length of the leg combined with appropriate body length can improve structural soundness and decrease the economic loss resulted from structural unsoundness for producers in pigs [4].

Longitudinal growth of the skeleton occurs through the action of chondrocytes in the proliferative and hypertrophic zones of the growth plate [5]. The cellular and biochemical processes influencing endochondral bone growth are complex and not yet fully elucidated. So far, there are few genetic studies on limb bone lengths in pigs. To our knowledge, only four quantitative trait loci (QTL) affecting femur dimensions has been detected on pig chromosomes (SSC) 2, 4, 16 and 17 [6]. Studies in mice [7-10] and chickens [11] have shown that the limb bone lengths are strongly controlled by genes. Elucidation of the genetic basis of bone growth will provide insight to understanding of the pathways and molecules involved in endochondral growth [12]. Bone mineral density (BMD) is a major determinant of risk for osteoporosis and bone fracture, which commonly affect the bone qualities of the distal forearm, thoracic and lumbar vertebrae, and proximal femur [13]. BMD is a complex trait with a high heritability, ranging from 0.5 to 0.9 in humans [14] and 0.6 to 0.7 in mice [15,16]. Furthermore, 60–85% and 80% of the phenotypic variance of BMD are genetically determined in humans [17] and mice [15], respectively. More recently, genome-wide scans for QTL affecting BMD have been widely performed in humans [18-20] and mice [15,16,21], and a number of QTL for BMD at different bone sites were mapped. In contrast, few QTL studies have been performed in farm animals [22], and none of QTL study was reported in pigs.

Because the physiologies are similar between pigs and humans, the pig is an ideal animal model for studying the genetic basis of bone growth and osteoporosis in humans, such as osteonecrosis of the femoral head, fractures of cartilage and bone, bone ingrowth [23,24]. The purpose of this study was to identify genomic regions affecting the lengths of limb bones (scapula, ulna, tibia, humerus and femur) and BMD of the distal femur in a large F2 resource population.

Results

Descriptive statistics of traits and correlations between traits

The descriptive statistics of the traits measured are listed in Table 1. Each limb bone was longer on average in castrated males than intact females. The correlation coefficients between the traits are shown in Table 2. All of the correlations were positive and highly significant (P < 0.001), especially for those among the lengths of limb bones. The BMD of the distal femur (fBMD) showed medium correlations with the lengths of limb bones.

Table 1. Descriptive statistics of limb bones lengths and BMD of the distal femur in the White Duroc × Erhualian intercross

Table 2. Phenotypic correlation coefficients among limb bone lengths and fBMDa

QTL for limb bone lengths

The results of QTL analysis are given in Table 3. The threshold values of 1% and 5% genome-wise significant and suggestive QTL for limb bone lengths were 7.53, 6.49, and 4.17, respectively. A total of 35 QTL were detected for the lengths of limb bones, and the number of QTL for each trait ranged from five (TL) to nine (SL). Due to the strong correlations among the limb bone lengths (Table 2), most of QTL regions showed associations with the lengths of two or more limb bones. Moreover, most of QTL had highly significant additive effects, but five QTL had significant dominance effects and a paternally expressed QTL affecting SL at 1% genome-wide significant level was found at 17 cM on SSC2 (Fig. 1A and Table 3).

Two QTL at 1% genome-wide significant level were found for the lengths of five limb bones (Table 3). One was located around 57 cM flanking by SW1856 and S0102 on SSC7, explaining 7 to 20% of the phenotypic variance (Fig. 1B). This QTL was the most significant QTL detected with small 95% confidence intervals (CI95) ranging from 2.5 (FL) to 5 cM (TL). The Erhualian alleles increased bone lengths at this locus. The other was in a region flanked by SW2456 and SW1943 on SSCX, explaining 2 to 6% of the phenotypic variance (Fig. 1C). At this locus, alleles from the White Duroc breed were associated with longer limb bones, and the CI95 varied from 12 (FL) to 59.5 cM (TL).

Three 1% genome-wide significant QTL each for SL, UL and TL were detected in an interval from 67 to 72 cM on SSC13 (Figure. 1D and Table 3). The favorable alleles were inherited from the White Duroc breed. A 1% genome-wide significant QTL for TL was found at the distal end of SSC1. In the same region, a 1% genome-wide significant QTL for UL and a 10% genome-wide significant QTL for HL were detected. Favorable QTL alleles at these loci were from the Erhualian breed. A suggestive QTL for SL and FL was evidenced at a different position on this chromosome (Fig. 1E and Table 3).

Table 3. The QTL mapping results for the lengths of limb bones and fBMD

thumbnailFigure 1. Evidence of significant QTL for lengths of the limb bone on pig chromosomes 2 (A), 7 (B), X (C), 13 (D), 1 (E) and 4 (F). The relative positions in cM on the linkage map are indicated in the x-axis, and the F values are given in the y-axis. Three lines are provided for 1% genome-wide (―), 5% genome-wide (----), and suggestive significance (– - – - – -) levels.

On SSC4 at 58 cM between SWR2179 and SW45, QTL was detected for HL and UL at 1% genome-wide significant level. At another region between SW45 and S0097 on this chromosome, two QTL were mapped for SL at 5% genome-wide significant level and for FL at suggestive level, respectively (Fig. 1F and Table 3). The White Duroc alleles were associated with longer HL, UL, SL and FL at these loci.

QTL for areal BMD

The critical values for 1%, 5% genome-wide and suggestive significance levels of fBMD QTL were 11.3, 8.16 and 4.56, respectively. Only two suggestive QTL were detected for this trait, one was located at 106 cM on SSC11 and the other was at 127 cM on SSC4 (Table 3). The two QTL jointly explained 25% of phenotypic variance, having significant additive and dominance effects whereas no imprinting effect. The Erhualian alleles were the favorable alleles for fBMD at the two loci.

Interaction between QTL and sex

Three of 37 QTL had significant interaction effects with sex (P < 0.05). One was a suggestive QTL for SL at 54 cM on SSC5 and the other two were 1% genome-wide significant QTL for UL at 58 cM and FL at 57 cM on SSC7, respectively (Table 4).

Table 4. The interactions between QTL and sex for limb bones

Discussion

In this study, significant segregation was observed in the bone-related traits measured in White Duroc × Erhualian intercross. Strong positive correlations among the lengths of five limb bones were found, indicating that causal genes have effects on general mechanisms of bone length growth. Accordingly, several QTL for the lengths of different limb bones were detected at the same positions, such as QTL on SSC7, SSCX and SSC13. A similar situation has been observed in the previous QTL analysis of limb bone lengths in mice [10].

A total of 37 QTL were evidenced including 35 QTL for the limb bone lengths and 2 QTL for fBMD. To our knowledge, it is the first time to report QTL for SL, HL, UL, TL, and fBMD in pigs. All QTL regions for FL, HL, TL and UL on autosomes are orthologous to the intervals harboring the loci for these long bones lengths in mouse models [7,9,10], whereas QTL for SL except that on SSC4 do not correspond to the SL QTL regions in mice [7]. QTL for femoral dimension have been detected on SSC2, 4, 16 and 17 in a Wild boar × Large White intercross [6], which was not repeated in this study. We instead found novel QTL for FL on SSC1, 3, 7, 10, 15 and X and a suggestive QTL for FL on SSC4 at a different position (35 cM vs 68 cM). The discordant QTL for FL are possibly caused by the different genetic background of founder animals in the two experimental populations or due to false discovery of QTL.

The QTL on SSC7 was the most significant QTL detected with a CI95 of less than 5 cM, providing a promising region to identify genes responsible for the length of limb bones. This genomic region is homologous to human chromosome (HSA) 6p21, which has significant effect on human stature by interaction with HSA2q21. A cluster of candidate genes for longitudinal or skeletal growth have been considered on HSA6p21, such as RUNX2/CBFA1 (runt-related transcription factor 2), COL11A2 (collagen, type XI, alpha 2) and RXRB (retinoid × receptor, beta) [25]. Another promising candidate gene is SCUBE3 (signal peptide-CUB-EGF-like domain-containing protein 3) that is proximal to the QTL peak and plays a critical role in bone cells with an exclusive expression in bones and osteoblasts [26]. Significant QTL for growth, carcass and fatness traits have been consistently mapped at positons close to the QTL region on SSC7 http://www.animalgenome.org/QTLdb/ webcite. It has been shown that body height has negative correlation with backfat thickness and has effect on the yield of ham, loin, picnic shoulder and shoulder butt [2,3]. The overlapping QTL region for limb bone length, growth and fatness traits indicates that there might be gene(s) that has pleiotropic effects on these traits in the region.

The QTL on SSCX is another region showing significant association with limb bone lengths. The region is homologous to HSAXq2.4, where a QTL for height and several QTL for syndromes of idiopathic short stature have been evidenced [27]. At the proximal end of the p-arm of SSC2, a 1% genome-wide significant QTL specific for SL was found with a significant maternal imprinting effect. Insulin-like growth factor 2 (IGF2) playing a major role in muscle growth and fat deposition is an interesting candidate gene for the QTL because its effect coincides with IGF2 imprinting effect. A single nucleotide substitution in IGF2 intron 3 has been identified to be the causative mutation for a major maternal-imprinting QTL in this region [28].

We only detected two suggestive QTLs for fBMD on SSC11 and SSC4, respectively. The QTL region on SSC4 is homologous to HSA1q21-23 encompassing a QTL for human spine BMD [19]. The QTL region on SSC11 is homologous to HSA13q32-34 harbouring a QTL for spine BMD [18] and a QTL for distal forearm areal BMD [20].

In this study, we performed QTL mapping by using combined-sex analyses. Taking the caution that QTL for these bone-related traits could be affected by sex, we analyzed the interaction effects between QTL and sex. Three QTL had significant interaction between QTL and sex. One was a suggestive QTL for SL on SSC5 that was expressed exclusively in intact females. The others were two genome-wide significant QTL each for UL and FL on SSC7, at which the dominant effect was expressed only in intact females. Nevertheless, these interaction effects only dropped the F-values without altering the evidence of QTL on this chromosome, confirming that the QTL mentioned above are reliable by the combined-sex analyses.

Conclusion

This study detected a total of 37 QTL for the lengths of scapula, ulna, humerus and tibia and BMD of the distal femur in pigs for the first time. Moreover, several new QTL for the femoral length in pigs were found. As correlated traits, many QTL for the lengths of five limb bones were located in the same genomic regions. The most promising QTL for each limb bone length were all located on SSC7 with the largest F-values and the smallest confidence intervals. It merits fine mapping of the QTL or the identification of positional candidate genes. Alleles from the White Duroc breed were not systematically favorable for longer length of limb bones.

Methods

Animals

A three-generation resource population was developed and managed as describled in our previous report [29]. Briefly, two White Duroc boars and 17 Erhualian sows were mated to produce 9 F1 boars and 59 F1 sows. These F1 animals were then intercrossed avoiding sister-brother mating to generate a total of 1912 F2 animals in 6 batches. In this study, 116 castrated males were recorded for fBMD and 1021 F2 animals including castrated males and intact female were measured for limb bone lengths. All F2 piglets were weaned at 45 d of age, and the males were castrated at the age of three months. All animals were housed in half-open pens and fed with mineral enforced corn-soybean diets.

Phenotype recording

Both forelimb and hind limb were removed from the right side carcass of F2 animals when slaughtered at 240 d of age. Five limb bones were dissected from the limbs and the lengths of these bones were measured with a large caliper: the scapula (the maximum straight line distance from the cavitas glenoidalis to the border of scapular cartilage), humerus (total length from the head to the trochlea), ulna (length from the olecranon process to the styloid process), femur (total length from the greater trochanter to the intercondyloid fossa) and the tibia (length from the intercondylar eminence to the medial malleolus). After that, the femur bones of the 116 F2 individuals in the first batch were stored at -20°C until utilized.

Areal bone mineral densities at the distal femur (fBMD) of the pigs were measured by a dual-energy X-ray absorptiometry (DXA) (Challenger, Montpellier, France) with a precision of ± 1.5%. At least 40% of each bone was scanned from the same end, and a common area was selected for fBMD measurement.

Genotyping and genetic map construction

Genome DNA was extracted from ear chip or spleen tissue, and a total of 183 informative microsatellites distributing 19 porcine chromosomes were genotyped for 1828 pigs including 19 founder animals, 68 F1 animals and 1741 F2 pigs in the White Duroc × Erhualian intercross. The number of markers on each chromosome varies from five (SSC18) to 24 (SSC13). Linkage analyses were carried out using the Crimap2.4 package [30]. Recombination units were transformed to map distances using the Haldane mapping function, and then the linkage maps were constructed.

Statistical analysis

Descriptive statistics of the traits measured in the White Duroc × Erhualian intercross and the correlation coefficients between traits were calculated by the MEANS and CORR procedures of SAS Version 9.0 (SAS Institute Inc., Cary, USA). The GLM procedure of SAS was employed to determine the fixed effects and the covariates in the final QTL mapping model. Factors having significant effect on traits were included in the QTL analysis model.

A least squares regression method was used to perform QTL analysis using QTLexpress at http://qtl.cap.ed.ac.uk webcite[31]. The underlying assumption of this method is that the alternative alleles at a given QTL are fixed in the founder breeds. The two QTL alleles were defined as D from White Duroc and E from Erhualian in this study. There are four possible genotypes (DD, DE, ED and EE) at each analysis point. The probabilities of the four QTL genotypes, viz. Prob(DD), Prob(DE), Prob(ED) and Prob(EE), for each F2 animal were inferred from the flanking markers at each centi-Morgan (cM) and the coefficients of additive, dominance and imprinting effects were equal to Prob(DD) minus Prob(EE), Prob(DE) plus Prob(ED) and Prob(DE) minus Prob(ED), respectively. Finally, the phenotypic values were regressed onto these coefficients to estimate the QTL effects at each cM across the genome. Family and sex were included in all QTL models as fixed effect. Batch and carcass length were considered respectively as fixed effect and covariate in the models for the lengths of limb bones. Additionally, some other special factors were also considered in the final model for QTL mapping as fixed effects or covariates including sex as fixed effect in analysis UL, live weight as a covariate in SL, FL and aBMD, and the age at slaughter as a covariate in SL. For analyses of the X chromosome, only the phenotypic data of female F2 progeny was considered in the QTL analysis like autosomal QTL scans.

To investigate whether the QTL had different effects between males and females, the interaction between QTL and sex was tested. Each detected QTL was fixed at the same position and reanalyzed with an interaction model. If the interaction model was significantly better than the non-interaction model, the interaction between QTL and sex was considered. The statistic of the test followed an F distribution, and the F value is the ratio of the deviation of the residual sum of squares (RSS) of the non-interaction model from the RSS of the interaction model to the residual mean square (MS) of the interaction model. The numerator and the denominator degree of freedom were equal to 3 (gender fixed in the non-interaction model) or 4 (gender not fixed in the non-interaction model) and the population size minus the number of effects was fixed in the interaction model, respectively.

A permutation method was used to get the empirical distribution of the test statistic by 1,000 times permutation [32], and the critical values were obtained based on the empirical distribution. The 95% confidence interval of QTL was constructed by a bootstrap method [33]. Based on 2,000 replacement resample, the empirical distribution of the QTL position was obtained and the 95% confidence interval was determined. Percentage of variance explained by each QTL was calculated using the following formula:

<a onClick="popup('http://www.biomedcentral.com/1471-2156/9/63/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2156/9/63/mathml/M1">View MathML</a>

Where MSfull, MSreduce1 and MSreduce were the mean squares of the models with all detected QTL, with the rest detected QTL except for the current focused one, and without all of the detected QTL, respectively.

Authors' contributions

HRM performed marker genotyping with BY and prepared the draft with YMG. YMG finished all of the data analysis, and recorded phenotypes with GCY and HSA. JR made a major revision of the manuscript. BB revised the draft. LSH conceived the research plan and made a final revision of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors give their appreciations to Xiao-mei Wang, form the Second Affiliated Hospital of Nanchang University, for her help in measurement of fBMD. This study was supported by the Natural Science Foundation of China (30425045) and the National 973 Program of China (2006CB708213).

References

  1. Mitchell AD, Scholz AM, Pursel VG: Total body and regional measurements of bone mineral content and bone mineral density in pigs by dual energy X-ray absorptiometry.

    J Anim Sci 2001, 79(10):2594-2604. PubMed Abstract | Publisher Full Text OpenURL

  2. Hetzer HO, Hankins OG, King JX, Zeller JH: Relationship between certain body measurements and carcass characteristics in swine.

    J Anim Sci 1950, 9:37-47. OpenURL

  3. Hetzer HO, Miller RH: Correlated responses of various body measurements in swine selected for high and low fatness.

    J Anim Sci 1972, 35(4):743-751. OpenURL

  4. Wood CM: Don't ignore feet and leg soundness in pigs.

    Livestock Update 2001. OpenURL

  5. Karsenty G, Wagner EF: Reaching a genetic and molecular understanding of skeletal development.

    Dev Cell 2002, 2(4):389-406. PubMed Abstract | Publisher Full Text OpenURL

  6. Andersson-Eklund L, Uhlhorn H, Lundeheim N, Dalin G, Andersson L: Mapping quantitative trait loci for principal components of bone measurements and osteochondrosis scores in a Wild Boar × Large White intercross.

    Genet Res 2000, 75(2):223-230. PubMed Abstract | Publisher Full Text OpenURL

  7. Christians JK, Bingham VK, Oliver FK, Heath TT, Keightley PD: Characterization of a QTL affecting skeletal size in mice.

    Mamm Genome 2003, 14(3):175-183. PubMed Abstract | Publisher Full Text OpenURL

  8. Drake TA, Hannani K, Kabo JM, Villa V, Krass K, Lusis AJ: Genetic loci influencing natural variations in femoral bone morphometry in mice.

    J Orthop Res 2001, 19(4):511-517. PubMed Abstract | Publisher Full Text OpenURL

  9. Kenney-Hunt JP, Vaughn TT, Pletscher LS, Peripato A, Routman E, Cothran K, Durand D, Norgard E, Perel C, Cheverud JM: Quantitative trait loci for body size components in mice.

    Mamm Genome 2006, 17(6):526-537. PubMed Abstract | Publisher Full Text OpenURL

  10. Leamy LJ, Pomp D, Eisen EJ, Cheverud JM: Pleiotropy of quantitative trait loci for organ weights and limb bone lengths in mice.

    Physiol Genomics 2002, 10(1):21-29. PubMed Abstract OpenURL

  11. Sharman PW, Morrice DR, Law AS, Burt DW, Hocking PM: Quantitative trait loci for bone traits segregating independently of those for growth in an F2 broiler × layer cross.

    Cytogenet Genome Res 2007, 117(1–4):296-304. PubMed Abstract | Publisher Full Text OpenURL

  12. Shiang R, Thompson LM, Zhu YZ, Church DM, Fielder TJ, Bocian M, Winokur ST, Wasmuth JJ: Mutations in the transmembrane domain of FGFR3 cause the most common genetic form of dwarfism, achondroplasia.

    Cell 1994, 78(2):335-342. PubMed Abstract | Publisher Full Text OpenURL

  13. Jin H, Ralston SH: Genetics of osteoporosis.

    Curr Rheumatol Rep 2005, 7(1):66-70. PubMed Abstract | Publisher Full Text OpenURL

  14. Huang QY, Kung AW: Genetics of osteoporosis.

    Mol Genet Metab 2006, 88(4):295-306. PubMed Abstract | Publisher Full Text OpenURL

  15. Beamer WG, Shultz KL, Churchill GA, Frankel WN, Baylink DJ, Rosen CJ, Donahue LR: Quantitative trait loci for bone density in C57BL/6J and CAST/EiJ inbred mice.

    Mamm Genome 1999, 10(11):1043-1049. PubMed Abstract | Publisher Full Text OpenURL

  16. Li X, Mohan S, Gu W, Wergedal J, Baylink DJ: Quantitative assessment of forearm muscle size, forelimb grip strength, forearm bone mineral density, and forearm bone size in determining humerus breaking strength in 10 inbred strains of mice.

    Calcif Tissue Int 2001, 68(6):365-369. PubMed Abstract | Publisher Full Text OpenURL

  17. Krall EA, Dawson-Hughes B: Heritable and life-style determinants of bone mineral density.

    J Bone Miner Res 1993, 8(1):1-9. PubMed Abstract OpenURL

  18. Deng HW, Xu FH, Huang QY, Shen H, Deng H, Conway T, Liu YJ, Liu YZ, Li JL, Zhang HT, Davies KM, Recker RR: A whole-genome linkage scan suggests several genomic regions potentially containing quantitative trait Loci for osteoporosis.

    J Clin Endocrinol Metab 2002, 87(11):5151-5159. PubMed Abstract | Publisher Full Text OpenURL

  19. Koller DL, Econs MJ, Morin PA, Christian JC, Hui SL, Parry P, Curran ME, Rodriguez LA, Conneally PM, Joslyn G, Peacock M, Johnston CC, Foroud T: Genome screen for QTLs contributing to normal variation in bone mineral density and osteoporosis.

    J Clin Endocrinol Metab 2000, 85(9):3116-3120. PubMed Abstract | Publisher Full Text OpenURL

  20. Niu T, Chen C, Cordell H, Yang J, Wang B, Wang Z, Fang Z, Schork NJ, Rosen CJ, Xu X: A genome-wide scan for loci linked to forearm bone mineral density.

    Hum Genet 1999, 104(3):226-233. PubMed Abstract | Publisher Full Text OpenURL

  21. Yu H, Mohan S, Edderkaoui B, Masinde GL, Davidson HM, Wergedal JE, Beamer WG, Baylink DJ: Detecting novel bone density and bone size quantitative trait loci using a cross of MRL/MpJ and CAST/EiJ inbred mice.

    Calcif Tissue Int 2007, 80(2):103-110. PubMed Abstract | Publisher Full Text OpenURL

  22. Rubin CJ, Brandstrom H, Wright D, Kerje S, Gunnarsson U, Schutz K, Fredriksson R, Jensen P, Andersson L, Ohlsson C, Mallmin H, Larsson S, Kindmark A: Quantitative trait loci for BMD and bone strength in an intercross between domestic and wildtype chickens.

    J Bone Miner Res 2007, 22(3):375-384. PubMed Abstract | Publisher Full Text OpenURL

  23. Pearce AI, Richards RG, Milz S, Schneider E, Pearce SG: Animal models for implant biomaterial research in bone: a review.

    Eur Cell Mater 2007, 13:1-10. PubMed Abstract | Publisher Full Text OpenURL

  24. Turner AS: Animal models of osteoporosis–necessity and limitations.

    Eur Cell Mater 2001, 1:66-81. PubMed Abstract | Publisher Full Text OpenURL

  25. Liu YZ, Guo YF, Xiao P, Xiong DH, Zhao LJ, Shen H, Liu YJ, Dvornyk V, Long JR, Deng HY, Li JL, Deng HW: Epistasis between loci on chromosomes 2 and 6 influences human height.

    J Clin Endocrinol Metab 2006, 91(10):3821-3825. PubMed Abstract | Publisher Full Text OpenURL

  26. Wu BT, Su YH, Tsai MT, Wasserman SM, Topper JN, Yang RB: A novel secreted, cell-surface glycoprotein containing multiple epidermal growth factor-like repeats and one CUB domain is highly expressed in primary osteoblasts and bones.

    J Biol Chem 2004, 279(36):37485-37490. PubMed Abstract | Publisher Full Text OpenURL

  27. Liu YZ, Xiao P, Guo YF, Xiong DH, Zhao LJ, Shen H, Liu YJ, Dvornyk V, Long JR, Deng HY, Li JL, Recker RR, Deng HW: Genetic linkage of human height is confirmed to 9q22 and Xq24.

    Hum Genet 2006, 119(3):295-304. PubMed Abstract | Publisher Full Text OpenURL

  28. Van Laere AS, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, Archibald AL, Haley CS, Buys N, Tally M, Andersson G, Georges M, Andersson L: A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig.

    Nature 2003, 425(6960):832-836. PubMed Abstract | Publisher Full Text OpenURL

  29. Ren J, Guo YM, Ma JW, Huang LS: Growth and meat quality QTL in pigs with special reference to a very large White Duroc × Erhualian resource population.

    Proc 8WCGALP Brazil poster ID 2006, 11-13. OpenURL

  30. Green P, Falls k, Crooks S: Documentation for CRIMAP, version 2.4. Washington University School of Medicine, St Louis; 1990.

  31. Haley CS, Knott SA, Elsen JM: Mapping quantitative trait loci in crosses between outbred lines using least squares.

    Genetics 1994, 136(3):1195-1207. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping.

    Genetics 1994, 138(3):963-971. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping.

    Genetics 1996, 143(2):1013-1020. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL