Copy Number Variation in CNP267 Region May Be Associated with Hip Bone Size

doi:10.1371/journal.pone.0022035

PLoS One. 2011; 6(7): e22035.

Published online 2011 July 15. doi: 10.1371/journal.pone.0022035

PMCID: PMC3137628

Copyright Liu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Copy Number Variation in CNP267 Region May Be Associated with Hip Bone Size

Shan-Lin Liu,^#¹ Shu-Feng Lei,^#¹^* Fang Yang,¹ Xi Li,¹ Rong Liu,¹ Shan Nie,¹ Xiao-Gang Liu,³ Tie-Lin Yang,³ Yan Guo,³ Fei-Yan Deng,² Qing Tian,² Jian Li,² Yao-Zhong Liu,² Yong-Jun Liu,² Hui Shen,^2,³ and Hong-Wen Deng^1,^2,^3,⁴^*

¹Laboratory of Molecular and Statistical Genetics and the Key Laboratory of Protein Chemistry and Developmental Biology of Ministry of Education, College of Life Sciences, Hunan Normal University, Changsha, Hunan, People's Republic of China

²Center of Bioinformatics and Genomics, School of Public Health and Tropical Medicine, Tulane University, New Orleans, Louisiana, United States of America

³School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, Shanxi, People's Republic of China

⁴Center of Systematic Biomedical Research, University of Shanghai for Science and Technology, Shanghai, People's Republic of China

Brian P. Chadwick, Editor

Florida State University, United States of America

^#Contributed equally.

* E-mail: hdeng2/at/tulane.edu (H-WD); Email: lei/at/hunnu.edu.cn (S-FL)

Conceived and designed the experiments: S-L. Liu S-F. Lei H-WD. Performed the experiments: X. Li X-G. Liu Yong-Jun Liu. Analyzed the data: F. Yang RL T-LY YG F-Y. Deng. Contributed reagents/materials/analysis tools: QT JL Yao-Zhong Liu HS. Wrote the paper: S-L. Liu S-F. Lei. Sample recruitment: SN.

Received March 8, 2011; Accepted June 13, 2011.

This article has been cited by other articles in PMC.

Abstract

Osteoporotic hip fracture (HF) is a serious global public health problem associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of key measurable risk factors for HF, independent of bone mineral density (BMD). Hip BS is highly genetically determined, but genetic factors underlying BS variation are still poorly defined. Here, we performed an initial genome-wide copy number variation (CNV) association analysis for hip BS in 1,627 Chinese Han subjects using Affymetrix GeneChip Human Mapping SNP 6.0 Array and a follow-up replicate study in 2,286 unrelated US Caucasians sample. We found that a copy number polymorphism (CNP267) located at chromosome 2q12.2 was significantly associated with hip BS in both initial Chinese and replicate Caucasian samples with p values of 4.73E-03 and 5.66E-03, respectively. An important candidate gene, four and a half LIM domains 2 (FHL2), was detected at the downstream of CNP267, which plays important roles in bone metabolism by binding to several bone formation regulator, such as insulin-like growth factor-binding protein 5 (IGFBP-5) and androgen receptor (AR). Our findings suggest that CNP267 region may be associated with hip BS which might influence the FHL2 gene downstream.

Introduction

Osteoporosis is a serious global public health problem especially in the elderly, which is characterized by low bone mineral density (BMD) and low trauma fracture. Osteoporotic hip fracture (HF) is the most serious low trauma fracture with high morbidity and mortality. It has been predicted that the number of osteoporotic fractures (OFs) worldwide will increase to 6.3 million in 2050, with most of these future fractures happening in Asian [1]. The economic burden is estimated at 17 billion dollars every year on OFs since 2005 in the United States alone [2]. Especially in recent years, an increasing aging Chinese population has resulted in even more OFs in China [3].

BMD, as an important risk factor for OFs, has been extensively investigated for identifying the genetic factors underlying BMD variation [4]. However, BMD is not the only risk factor for osteoporotic fracture, and recent studies suggested that bone size (BS) is another important risk factor for OFs independent of BMD [5], [6]. BS has more than 50% of heritability in both Chinese and Caucasians [7], [8]. Abnormal BS contributes significantly to the pathogenesis of OFs [5], [6]. Bigger bone size is a logical adaptation to enhance the mechanical competence of bone, because a larger cross-sectional area can bear larger compressive loads and cope more efficiently with bending loading. Spine BMD in females is comparable with that in males, however, females suffer a larger incidence of spine fractures than males partially attributable to the fact that female's spine bone sizes are 20–25% smaller than male's after adjusting for body size differences [9]. A genetically homogenous inbred mouse strain has higher bone mass but smaller bone size, and is less sensitive in adapting to mechanical loading to increase bone strength when compared with another inbred mouse strain [10]. Currently, only a few of genes such as VDR [11], [12], ER, COL1A2 [13], CYP17 [14], PTH [15], were tested in association with BS, and most of the genetic factors for BS largely remain unknown.

Copy number variation (CNV) is a type of DNA variation due to diversity in the number of copies of a DNA segment that may range from one kilobase to several megabases in size. Copy number polymorphisms (CNPs) refer to common CNVs that appear to involve the same affected genomic sequence and are therefore consistent with a model of a genetic polymorphism. CNVs may alter gene dosage, disrupt coding sequences, or exert long-range positional effects on gene expression pattern outside CNV region, consequently leading to phenotype variation [16], [17]. Individuals with lower copy number of CCL3L1 [18], FCGR3B [19] and DEFB4 [20] genes were predisposed to risk of AIDS, immunologically mediated glomerulonephritis and Crohn disease, respectively. Recently, our groups performed two genome-wide CNV association studies on BMD in Chinese and Caucasian and found CNV regions containing UGT2B17 [21] and VPS13B [22] genes were significantly associated with BMD, respectively. However, to the best of our knowledge, there is no reported genome-wide CNV association research that focused on hip BS variation.

This study performed an initial genome-wide common CNV association analysis in 1,627 Chinese Han subjects using Affymetrix GeneChip Human Mapping SNP 6.0 array and a follow-up replicate study in 2,286 unrelated US Caucasians sample, and identified CNP267 region may be associated with hip BS, which might influence the FHL2 gene downstream.

Materials and Methods

Ethics Statement

The study was approved by the necessary institutional review boards of all the participating institutions. Before entering the project, all the subjects signed the informed-consent documents.

Subjects and Measurement

We first performed an initial genome-wide CNV analysis for hip BS in 1,627 unrelated Chinese Han subjects. The most significant results identified in Chinese sample were further replicated in an unrelated Caucasian sample.

1) Initial study sample

The sample for the initial genome-wide CNV analysis consisted of 1,627 unrelated Chinese Han adults including 825 women and 802 men. All the subjects were recruited from the cities of Xi'an/Changsha and their vicinity with the purposes of searching for genetic factors for osteoporosis related phenotypes (such as BMD, BS, bone geometry). In order to minimize the confounding effects of environmental and therapeutic factors that may influence the skeletal system, individuals with chronic diseases and conditions that might potentially affect bone metabolism were excluded. These diseases/conditions included chronic disorders involving vital organs (heart, lung, liver, kidney, and brain), serious metabolic diseases (diabetes, hypo- and hyper-parathyroidism, etc.), other skeletal diseases (Paget disease, osteogenesis imperfecta, rheumatoid arthritis, etc.), chronic use of drugs affecting bone metabolism (hormone replacement therapy, corticosteroid therapy, anti-convulsant drugs), and malnutrition conditions (such as chronic diarrhea, chronic ulcerative colitis, etc.), etc. In addition, subjects taking anti-bone-resorptive or bone anabolic agents/drugs, such as bisphosphonates, were excluded from this study.

Hip BS (total area including femoral neck, trochanter and intertrochanter) was measured by a Hologic QDR 4500W dual-energy X-ray absorptiometry (DXA) scanner (Hologic Corp., Waltham, MA). The machine was calibrated daily, and the coefficient of variation (CV) for hip BS value, obtained from seven individuals repeatedly measured five times, of the DXA measurements was 1.18%. Height and weight were measured at the same visit for BS measurement.

2) Caucasian replication sample

To replicate the associations for CNPs at significant level of p<0.01, we performed another independent study in a sample including a total of 2,286 random Caucasian subjects (age: 51.4±13.8 yrs) that were recruited in Midwestern US in Kansas City, Missouri and Omaha, Nebraska for osteoporosis study. All identified subjects were US Caucasians of European origin. The BS values of total hip were measured using DXA machines (Hologic Inc., Bedford, MA, USA) that were calibrated daily. The CV of the DXA measurement for total hip BS was about 1.94%.

Genotyping

1) In initial study sample

Genomic DNA was extracted from peripheral blood leukocytes using standard protocols. Genotyping with Affymetrix Genome-Wide Human SNP Array 6.0, which features 1.8 million genetic markers, including more than 906,600 SNPs and more than 946,000 probes for detection of copy number variation, was performed using the standard protocol recommended by the manufacturer. Fluorescence intensities were quantified using an Affymetrix array scanner 30007G. Data management and analyses were performed using the Affymetrix® GeneChip® Command Console® Software (AGCC). Contrast quality control (QC) threshold was set at the default value of greater than 0.4 for sample quality control. The final average contrast QC across the entire sample reached the high level of 2.62. The Birdsuite package (http://www.broadinstitute.org/science/programs/medical-and-population-genetics/birdsuite/birdsuite-0) was used for genotype calling, genotyping quality control, and CNV identification.

For the QC for sample, we firstly measured the copy number estimates of each chromosome and genome-wide average (sum of all chromosomes), reported by the Birdseye Hidden Markov Model [23], and removed the subjects who showed excessively high or low estimate for copy number according to either genome-wide average or more than 2 chromosomes (>3 standard deviations). Then, we measured the variability of CNP and SNP probe intensities according to each chromosome and genome-wide average (sum of all chromosomes). We removed the subjects with excessive variability in probe intensity according to either genome-wide average or more than 2 chromosomes (>3 standard deviations). We kept the subjects who only had 1 or 2 chromosomes failing in copy number estimate QC and probe intensity QC, and treated the CNPs in the chromosomes of those subjects as missing data in further association analysis. As a result, 1,531 samples were used in CANARY software [23] for CNP call.

For the QC for CNPs, we removed the CNPs with more than 5% of uncertain or missing copy call (CC) or with less than 1% of allele frequency (AF) that refers to the total proportion of the subjects with copy number less or more than two in total sample. Of the initial full-set of 1,316 CNPs, there are 817 CNPs with AF <1%, 735 CNPs with CC greater than 0.05 and 430 CNPs with both. Finally 194 CNPs covering ~12 Mb with an average marker spacing of 64 kb remain in the final association analyses (Table S1).

2) In replication sample

Genomic DNA was extracted from human blood using a commercial isolation kit (Gentra systems, Minneapolis, MN, USA) following the protocols detailed in the kit. Genotyping and quality control procedures are the same as in Chinese.

Statistical analysis

We used stepwise regression model to screen significant covariates for each study cohort. Parameters including age, age², sex, age-sex, age²-sex, weight, height, BMI, birth year were tested for their association with hip BS. Significant (p≤0.05) parameters (weight, height in initial Chinese sample and age, sex, age-sex, age²-sex, height, birth year in replication Caucasian sample) were then included as covariates to adjust the raw BS values. EIGENSTRAT was employed to perform principal component analysis to correct for stratification in genome-wide association studies [24]. We used ~370,000 SNPs to calculate the principal components and the ten default main eigenvectors were used as covariates to adjust the raw BS values for correction of population stratification. The adjusted BS data, if not following normal distributions, were further subjected to BoxCox transformation into normal distribution. Finally, association analyses between CNPs and hip BS were performed using the PLINK software package (version 1.07) (http://pngu.mgh.harvard.edu/~purcell/plink/).

The CNPs at the significant level of p<0.01 identified in Chinese sample were selected in replication study in Caucasian. The processes of raw hip BS value adjustment and association analyses in replication study were the same as in initial discovery study in Chinese.

Results

The basic characteristics of the subjects used in association analyses and replication analyses, including sex, age, height, weight, hip BS were summarized in Table 1.

Table 1

Basic characteristics of the study subjects.

In the initial genome-wide CNV analysis in Chinese sample, four CNPs were associated with hip BS at significant level of p<0.01 (CNP11164 P

6.18E-04, CNP182 P

1.05E-03, CNP267 P

4.73E-03, CNP10799 P

7.23E-03) (Table 2). To find replication association evidence, we performed a subsequent replication analysis for these four CNPs in Caucasians that have evidently different genetic background from Chinese. Due to AF of CNP11164 (AF

0.0031) and CNP10799 (AF

0.003) is smaller than 1%, only CNP182 and CNP267 were analyzed for replicate association and reached p values at 2.39E-01 and 5.66E-03, respectively (Table 2). Meta-analysis (Fisher's combined p method [25]) was used to combine association tests in the two populations. The Fisher's combined p value (p

3.09E-04) for CNP267 is nearly approaching the significant level after multiple testing (0.05/194

2.58E-04).

Table 2

Characteristics of the interesting CNPs for association analysis.

CNP267 was located in 2q12.2 with physical position from 106,247,145 bp to 106,251,789 bp. The four and a half LIM domains 2 (FHL2) located at the downstream of CNP267, is an important candidate gene, which plays important roles in bone metabolism by binding to several bone formation regulator (details in Discussion). Of the subjects analyzed in the initial and replication samples, carriers of copy number (CN) of 0CN, 1CN, 2CN and 3CN, were 210, 669, 513, 5 in Chinese sample and 214, 910, 914, 7 in Caucasian sample, respectively. Cause the number of carriers of CN

3 is too small to be statistic significant, we excluded the data of CN

3 in the following analysis. For the effect of CNP267, we compared the adjusted hip BS for subjects with different CNs of CNP267 in both the initial and replication studies. As shown in Figure 1, CN was negatively associated with hip BS in Chinese sample (i.e., individuals with more CNs trend to have lower hip BS), but positively in Caucasian sample. The opposite effect of CN on BS variation may be caused by several reasons (details in Discussion). Regression analysis showed that CNP267 can explain ~0.97% of hip BS variation.

Figure 1

Hip BS values in groups with different copy number (CN) of CNP267 in the Chinese and Caucasians.

Discussion

CNV is a type of recently increasingly being recognized DNA variation associated with human complex diseases. CNV presumably alters gene dosage, disrupts coding sequences, or exerts long-range positional effects on the gene expression pattern outside the CNV region, and consequently leads to phenotype variation. We performed the first genome-wide CNV association study in a large Chinese sample for an important bone strength phenotype, hip BS, and a follow-up replicate study in 2,286 unrelated US Caucasians sample, and identified that CNP267 region may be associated with hip BS. We also tried to find if there are any BS candidate genes in the regions of significant CNPs in initial Chinese sample. Although there were no well known BS candidate genes, such as VDR, ER, COL1A2, CYP17, PTH et al., we suggested that these genes should deserve more attention in further studies and presented the results in Table S2.

CNP267 was located in 2q12.2 with physical position from 106,247,145 bp to 106,251,789 bp, where no important candidate gene has known function directly associated with osteoporosis or bone. However in the downstream of CNP267 there is an important candidate gene, the four and a half LIM domains 2 (FHL2), which plays important roles in bone metabolism by binding to several bone formation regulator. FHL2 is strongly expressed in all types of human osteoblasts [26]. FHL2 can bind to several regulation factors that are important during bone formation. For example, it binds to insulin-like growth factor-binding protein 5 (IGFBP-5) [26], the most abundant IGFBP stored in bone, which plays an important role in regulating bone formation [27], [28], [29], [30]. FHL2 also increases transcriptional activity of androgen receptor (AR) [31], which is important in regulation of bone formation and resorption [32], [33]. In addition, FHL2 can also influence bone formation as one of modulators of the Wnt signaling pathway [34], [35] and TNF receptor–associated factor 6 (TRAF6)-mediated receptor activator of NF-κB (RANK) signaling pathway [36]. Gene knockdown studies showed that FHL2-deficient mice have decreased BMD and BS due to decreased osteoblast activity [37], [38]. Conversely, enforced expression of FHL2 boosts mineralization in cell culture and increased bone mass in transgenic mice [37].

Genes located outside of CNV intervals can show disrupted expression patterns by long-range position effects. For example, transcription levels of some genes located in the diploid flanking sequence outside deletion region were altered in patients carrying the Williams-Beuren syndrome deletion [39]. CNV17p11.2 can change the expression of some flanking genes, and thus cause Smith-Magenis syndrome [40]. Whole genome association studies between CNVs and gene expression levels showed that more than half of identified genes significantly associated with CNVs were not located in the CNV intervals [41]. These observations may change the traditional view, i.e., CNVs regulate gene expression via directly altering copy number of dosage-sensitive genes located in CNV intervals. However CNVs may alter the expression of genes located outside of CNV intervals by disrupting local chromatin environment and/or gene regulatory elements such as enhancers, silencers, and insulators. Previous studies reported that the relative expression levels of genes were modified through changes in copy numbers as far as 2–7 Mb away from the breakpoints [41], [42].

In association study, the inconsistence of association directions between initial and replicate association studies is usually observed. There may be several reasons to address the opposite effect of CN on BS variation. First, it is well known that the same risk locus with diverse genetic background may generate different effects [43], [44]. Studies also showed that osteoporosis related phenotypes are clearly under ethnic specific genetic determination [45], [46], [47]. Chinese and Caucasians are two evidently different populations in genetic background. The comparison of frequencies of CNs in the two samples shows significantly different with P<0.01 in chi-square test. The ethnic difference may lead to the reverse effect of association in our study. Second, in most situations, significant association suggests strong linkage disequilibrium between marker and causal locus. If the causal loci are different in two samples, it is possible to cause discordant direction of association. For example, the risk allele of marker is negatively correlated with a risk allele at one causal locus in one population, but positively associated with a protective allele at another causal locus in another population. Czarnomska et al. also reported such observation, i.e., a set of genes control the impact of the ApcMin mutation in two organs but with opposite effects [48]. Finally, environmental covariates may influence the effect direction of gene variants. Eder et al. discovered opposite effects of CD14/−260 on serum IgE levels in children raised in different environments [49]. Reneland et al. found that rs1498608 in PDE4D gene showed an opposite relationship with BMD variation, which indicating that the variant's effect may be context-dependent [50]. The two samples come from different geographic regions, suggesting obviously different environment such as living condition and dietary habit. These differences may also influence the effect direction of the allele.

In this study, FHL2 may be recognized as an important candidate gene involved in the association between CNP267 and hip BS by the potential mechanism that CNP267 may influence the expression of FHL2 by causing changes in local chromatin structure or alterations in enhancer activity. This potential mechanism is waiting for being confirmed by functional studies.

In conclusion, this is the first genome-wide CNV association study for hip BS in Chinese Han sample. CNP267 was significantly associated with hip BS in both Chinese and Caucasian. Our findings suggest that CNP267 region may be associated with hip BS.

Supporting Information

Table S1

Information of other 190 CNPs with p>0.01 in initial Chinese association analysis.

(DOC)

Click here for additional data file.^{(241K, doc)}

Table S2

Candidate genes at the interesting CNP regions in initial Chinese sample.

(DOC)

Click here for additional data file.^{(38K, doc)}

Acknowledgments

We would like to thank Fan Yin and Quan-wei Liu for routine maintenance of all the instruments in our laboratory.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: The study was partially supported by Natural Science Foundation of China (NSFC) (31071097, 30600364, 30771222 and 30900810), NSFC-Canadian Institutes of Health Research (CIHR) Joint Health Research Initiative Proposal (30811120436), NSFC/Research Grants Council Joint Research Scheme (30731160618), Shanghai Leading Academic Discipline Project (S30501) and startup fund from Shanghai University of Science and Technology. H-WD was partially supported by grants from National Institutes of Health (P50AR055081, R01AG026564, R01AR050496, RC2DE020756, R01AR057049, and R03TW008221) and Franklin D. Dickson/Missouri Endowment. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Liu Z, Piao J, Pang L, Qing X, Nan S, et al. The diagnostic criteria for primary osteoporosis and the incidence of osteoporosis in China. J Bone Miner Metab. 2002;20:181–189. [PubMed]

2. Burge R, wson-Hughes B, Solomon DH, Wong JB, King A, et al. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005–2025. J Bone Miner Res. 2007;22:465–475. [PubMed]

3. Deng HW, Deng XT, Conway T, Xu FH, Heaney R, et al. Determination of bone size of hip, spine, and wrist in human pedigrees by genetic and lifestyle factors. J Clin Densitom. 2002;5:45–56. [PubMed]

4. Liu YJ, Shen H, Xiao P, Xiong DH, Li LH, et al. Molecular genetic studies of gene identification for osteoporosis: a 2004 update. J Bone Miner Res. 2006;21:1511–1535. [PMC free article] [PubMed]

5. Seeman E, Duan Y, Fong C, Edmonds J. Fracture site-specific deficits in bone size and volumetric density in men with spine or hip fractures. J Bone Miner Res. 2001;16:120–127. [PubMed]

6. Deng HW, Xu FH, Davies KM, Heaney R, Recker RR. Differences in bone mineral density, bone mineral content, and bone areal size in fracturing and non-fracturing women, and their interrelationships at the spine and hip. J Bone Miner Metab. 2002;20:358–366. [PubMed]

7. Jian WX, Long JR, Deng HW. High heritability of bone size at the hip and spine in Chinese. J Hum Genet. 2004;49:87–91. [PubMed]

8. Cooper C, Campion G, Melton LJ., III Hip fractures in the elderly: a world-wide projection. Osteoporos Int. 1992;2:285–289. [PubMed]

9. Gilsanz V, Boechat MI, Gilsanz R, Loro ML, Roe TF, et al. Gender differences in vertebral sizes in adults: biomechanical implications. Radiology. 1994;190:678–682. [PubMed]

10. Akhter MP, Cullen DM, Pedersen EA, Kimmel DB, Recker RR. Bone response to in vivo mechanical loading in two breeds of mice. Calcif Tissue Int. 1998;63:442–449. [PubMed]

11. Need AG, Horowitz M, Stiliano A, Scopacasa F, Morris HA, et al. Vitamin D receptor genotypes are related to bone size and bone density in men. Eur J Clin Invest. 1996;26:793–796. [PubMed]

12. Fang Y, van Meurs JB, Rivadeneira F, van Schoor NM, van Leeuwen JP, et al. Vitamin D receptor gene haplotype is associated with body height and bone size. J Clin Endocrinol Metab. 2007;92:1491–1501. [PubMed]

13. Lei SF, Deng FY, Xiao SM, Chen XD, Deng HW. Association and haplotype analyses of the COL1A2 and ER-alpha gene polymorphisms with bone size and height in Chinese. Bone. 2005;36:533–541. [PubMed]

14. Zmuda JM, Cauley JA, Kuller LH, Ferrell RE. A common promotor variant in the cytochrome P450c17alpha (CYP17) gene is associated with bioavailability testosterone levels and bone size in men. J Bone Miner Res. 2001;16:911–917. [PubMed]

15. Gong G, Johnson ML, Barger-Lux MJ, Heaney RP. Association of bone dimensions with a parathyroid hormone gene polymorphism in women. Osteoporos Int. 1999;9:307–311. [PubMed]

16. Sharp AJ. Emerging themes and new challenges in defining the role of structural variation in human disease. Hum Mutat. 2009;30:135–144. [PubMed]

17. Lupski JR, Stankiewicz P. Genomic disorders: molecular mechanisms for rearrangements and conveyed phenotypes. PLoS Genet. 2005;1:e49. [PMC free article] [PubMed]

18. Shao W, Tang J, Song W, Wang C, Li Y, et al. CCL3L1 and CCL4L1: variable gene copy number in adolescents with and without human immunodeficiency virus type 1 (HIV-1) infection. Genes Immun. 2007;8:224–231. [PubMed]

19. Fanciulli M, Norsworthy PJ, Petretto E, Dong R, Harper L, et al. FCGR3B copy number variation is associated with susceptibility to systemic, but not organ-specific, autoimmunity. Nat Genet. 2007;39:721–723. [PMC free article] [PubMed]

20. Boniotto M, Ventura M, Eskdale J, Crovella S, Gallagher G. Evidence for duplication of the human defensin gene DEFB4 in chromosomal region 8p22-23 and implications for the analysis of SNP allele distribution. Genet Test. 2004;8:325–327. [PubMed]

21. Yang TL, Chen XD, Guo Y, Lei SF, Wang JT, et al. Genome-wide copy-number-variation study identified a susceptibility gene, UGT2B17, for osteoporosis. Am J Hum Genet. 2008;83:663–674. [PubMed]

22. Deng FY, Zhao LJ, Pei YF, Sha BY, Liu XG, et al. Genome-wide copy number variation association study suggested VPS13B gene for osteoporosis in Caucasians. Osteoporos Int. 2009;21:579–87. [PMC free article] [PubMed]

23. Korn JM, Kuruvilla FG, McCarroll SA, Wysoker A, Nemesh J, et al. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet. 2008;40:1253–1260. [PMC free article] [PubMed]

24. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, et al. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–909. [PubMed]

25. Fisher KA, Miles R. Modeling the acoustic radiation force in microfluidic chambers. J Acoust Soc Am. 2008;123:1862–1865. [PubMed]

26. Amaar YG, Thompson GR, Linkhart TA, Chen ST, Baylink DJ, et al. Insulin-like growth factor-binding protein 5 (IGFBP-5) interacts with a four and a half LIM protein 2 (FHL2). J Biol Chem. 2002;277:12053–12060. [PubMed]

27. Mohan S, Nakao Y, Honda Y, Landale E, Leser U, et al. Studies on the mechanisms by which insulin-like growth factor (IGF) binding protein-4 (IGFBP-4) and IGFBP-5 modulate IGF actions in bone cells. J Biol Chem. 1995;270:20424–20431. [PubMed]

28. Richman C, Baylink DJ, Lang K, Dony C, Mohan S. Recombinant human insulin-like growth factor-binding protein-5 stimulates bone formation parameters in vitro and in vivo. Endocrinology. 1999;140:4699–4705. [PubMed]

29. Govoni KE, Baylink DJ, Mohan S. The multi-functional role of insulin-like growth factor binding proteins in bone. Pediatr Nephrol. 2005;20:261–268. [PMC free article] [PubMed]

30. Andress DL. IGF-binding protein-5 stimulates osteoblast activity and bone accretion in ovariectomized mice. Am J Physiol Endocrinol Metab. 2001;281:E283–E288. [PubMed]

31. Muller JM, Isele U, Metzger E, Rempel A, Moser M, et al. FHL2, a novel tissue-specific coactivator of the androgen receptor. EMBO J. 2000;19:359–369. [PubMed]

32. Yeh S, Tsai MY, Xu Q, Mu XM, Lardy H, et al. Generation and characterization of androgen receptor knockout (ARKO) mice: an in vivo model for the study of androgen functions in selective tissues. Proc Natl Acad Sci U S A. 2002;99:13498–13503. [PubMed]

33. Kawano H, Sato T, Yamada T, Matsumoto T, Sekine K, et al. Suppressive function of androgen receptor in bone resorption. Proc Natl Acad Sci U S A. 2003;100:9416–9421. [PubMed]

34. Labalette C, Renard CA, Neuveut C, Buendia MA, Wei Y. Interaction and functional cooperation between the LIM protein FHL2, CBP/p300, and beta-catenin. Mol Cell Biol. 2004;24:10689–10702. [PMC free article] [PubMed]

35. Martin B, Schneider R, Janetzky S, Waibler Z, Pandur P, et al. The LIM-only protein FHL2 interacts with beta-catenin and promotes differentiation of mouse myoblasts. J Cell Biol. 2002;159:113–122. [PMC free article] [PubMed]

36. Bai S, Kitaura H, Zhao H, Chen J, Muller JM, et al. FHL2 inhibits the activated osteoclast in a TRAF6-dependent manner. J Clin Invest. 2005;115:2742–2751. [PMC free article] [PubMed]

37. Gunther T, Poli C, Muller JM, Catala-Lehnen P, Schinke T, et al. Fhl2 deficiency results in osteopenia due to decreased activity of osteoblasts. EMBO J. 2005;24:3049–3056. [PubMed]

38. Govoni KE, Baylink DJ, Chen J, Mohan S. Disruption of four-and-a-half LIM 2 decreases bone mineral content and bone mineral density in femur and tibia bones of female mice. Calcif Tissue Int. 2006;79:112–117. [PMC free article] [PubMed]

39. Merla G, Howald C, Henrichsen CN, Lyle R, Wyss C, et al. Submicroscopic deletion in patients with Williams-Beuren syndrome influences expression levels of the nonhemizygous flanking genes. Am J Hum Genet. 2006;79:332–341. [PubMed]

40. Bi W, Yan J, Stankiewicz P, Park SS, Walz K, et al. Genes in a refined Smith-Magenis syndrome critical deletion interval on chromosome 17p11.2 and the syntenic region of the mouse. Genome Res. 2002;12:713–728. [PubMed]

41. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315:848–853. [PMC free article] [PubMed]

42. Reymond A, Henrichsen CN, Harewood L, Merla G. Side effects of genome structural changes. Curr Opin Genet Dev. 2007;17:381–386. [PubMed]

43. Miao X, Garcia-Barcelo MM, So MT, Tang WK, Dong X, et al. Lack of association between nNOS −84G>A polymorphism and risk of infantile hypertrophic pyloric stenosis in a Chinese population. J Pediatr Surg. 2010;45:709–713. [PubMed]

44. Chang XL, Mao XY, Li HH, Zhang JH, Li NN, et al. Functional parkin promoter polymorphism in Parkinson's disease: new data and meta-analysis. J Neurol Sci. 2011;302:68–71. [PubMed]

45. Dvornyk V, Liu PY, Long JR, Zhang YY, Lei SF, et al. Contribution of genotype and ethnicity to bone mineral density variation in Caucasians and Chinese: a test for five candidate genes for bone mass. Chin Med J (Engl) 2005;118:1235–1244. [PubMed]

46. Dvornyk V, Liu XH, Shen H, Lei SF, Zhao LJ, et al. Differentiation of Caucasians and Chinese at bone mass candidate genes: implication for ethnic difference of bone mass. Ann Hum Genet. 2003;67:216–227. [PubMed]

47. Lei SF, Chen Y, Xiong DH, Li LM, Deng HW. Ethnic difference in osteoporosis-related phenotypes and its potential underlying genetic determination. J Musculoskelet Neuronal Interact. 2006;6:36–46. [PubMed]

48. Czarnomska A, Krysiak E, Piskorowska J, Sitarz M, Pysniak K, et al. Opposite effects of modifiers of the ApcMin mutation in intestine and mammary gland. Cancer Res. 2003;63:4533–4537. [PubMed]

49. Eder W, Klimecki W, Yu L, von ME, Riedler J, et al. Opposite effects of CD 14/−260 on serum IgE levels in children raised in different environments. J Allergy Clin Immunol. 2005;116:601–607. [PubMed]

50. Reneland RH, Mah S, Kammerer S, Hoyal CR, Marnellos G, et al. Association between a variation in the phosphodiesterase 4D gene and bone mineral density. BMC Med Genet. 2005;6:9. [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science