Tie-Lin Yang

Genome-wide Copy-Number-Variation Study Identified a Susceptibility Gene, UGT2B17, for Osteoporosis

Osteoporosis, a highly heritable disease, is characterized mainly by low bone-mineral density (BMD), poor bone geometry, and/or osteoporotic fractures (OF). Copy-number variation (CNV) has been shown to be associated with complex human diseases. The contribution of CNV to osteoporosis has not been determined yet. We conducted case-control genome-wide CNV analyses, using the Affymetrix 500K Array Set, in 700 elderly Chinese individuals comprising 350 cases with homogeneous hip OF and 350 matched controls. We constructed a genomic map containing 727 CNV regions in Chinese individuals. We found that CNV 4q13.2 was strongly associated with OF (p = 2.0 x 10(-4), Bonferroni-corrected p = 0.02, odds ratio = 1.73). Validation experiments using PCR and electrophoresis, as well as real-time PCR, further identified a deletion variant of UGT2B17 in CNV 4q13.2. Importantly, the association between CNV of UGT2B17 and OF was successfully replicated in an independent Chinese sample containing 399 cases with hip OF and 400 controls. We further examined this CNVs relevance to major risk factors for OF (i.e., hip BMD and femoral-neck bone geometry) in both Chinese (689 subjects) and white (1000 subjects) samples and found consistently significant results (p = 5.0 x 10(-4) -0.021). Because UGT2B17 encodes an enzyme catabolizing steroid hormones, we measured the concentrations of serum testosterone and estradiol for 236 young Chinese males and assessed their UGT2B17 copy number. Subjects without UGT2B17 had significantly higher concentrations of testosterone and estradiol. Our findings suggest the important contribution of CNV of UGT2B17 to the pathogenesis of osteoporosis.

Genome-wide Association and Follow-Up Replication Studies Identified ADAMTS18 and TGFBR3 as Bone Mass Candidate Genes in Different Ethnic Groups

To identify and validate genes associated with bone mineral density (BMD), which is a prominent osteoporosis risk factor, we tested 379,319 SNPs in 1000 unrelated white U.S. subjects for associations with BMD. For replication, we genotyped the most significant SNPs in 593 white U.S. families (1972 subjects), a Chinese hip fracture (HF) sample (350 cases, 350 controls), a Chinese BMD sample (2955 subjects), and a Tobago cohort of African ancestry (908 males). Publicly available Framingham genome-wide association study (GWAS) data (2953 whites) were also used for in silico replication. The GWAS detected two BMD candidate genes, ADAMTS18 (ADAM metallopeptidase with thrombospondin type 1 motif, 18) and TGFBR3 (transforming growth factor, beta receptor III). Replication studies verified the significant findings by GWAS. We also detected significant associations with hip fracture for ADAMTS18 SNPs in the Chinese HF sample. Meta-analyses supported the significant associations of ADAMTS18 and TGFBR3 with BMD (p values: 2.56 x 10(-5) to 2.13 x 10(-8); total sample size: n = 5925 to 9828). Electrophoretic mobility shift assay suggested that the minor allele of one significant ADAMTS18 SNP might promote binding of the TEL2 factor, which may repress ADAMTS18 expression. The data from NCBI GEO expression profiles also showed that ADAMTS18 and TGFBR3 genes were differentially expressed in subjects with normal skeletal fracture versus subjects with nonunion skeletal fracture. Overall, the evidence supports that ADAMTS18 and TGFBR3 might underlie BMD determination in the major human ethnic groups.

Robust and Comprehensive Analysis of 20 Osteoporosis Candidate Genes by Very High-Density Single-Nucleotide Polymorphism Screen Among 405 White Nuclear Families Identified Significant Association and Gene–Gene Interaction

Many “novel” osteoporosis candidate genes have been proposed in recent years. To advance our knowledge of their roles in osteoporosis, we screened 20 such genes using a set of high‐density SNPs in a large family‐based study. Our efforts led to the prioritization of those osteoporosis genes and the detection of gene–gene interactions.

Polymorphisms of the Low-density lipoprotein receptor-related protein 5 (LRP5) gene are associated with obesity phenotypes in a large family-based association study

Background: The low-density lipoprotein receptor-related protein 5 (LRP5) gene, essential for glucose and cholesterol metabolism, may have a role in the aetiology of obesity, an important risk factor for diabetes. Participants and methods: To investigate the association between LRP5 polymorphisms and obesity, 27 single-nucleotide polymorphisms (SNPs), spacing about 5 kb apart on average and covering the full transcript length of the LRP5 gene, were genotyped in 1873 Caucasian people from 405 nuclear families. Obesity (defined as body mass index (BMI) >30 kg/m2) and three obesity-related phenotypes (BMI, fat mass and percentage of fat mass (PFM)) were investigated. Results: Single markers (12 tagging SNPs and 4 untaggable SNPs) and haplotypes (5 blocks) were tested for associations, using family-based designs. SNP4 (rs4988300) and SNP6 (rs634008) located in block 2 (intron 1) showed significant associations with obesity and BMI after Bonferroni correction (SNP4: p<0.001 and p = 0.001, respectively; SNP6: p = 0.002 and 0.003, respectively). The common allele A for SNP4 and minor allele G for SNP6 were associated with an increased risk of obesity. Significant associations were also observed between common haplotype A–G–G–G of block 2 with obesity, BMI, fat mass and PFM with global empirical values p<0.001, p<0.001, p = 0.003 and p = 0.074, respectively. Subsequent sex-stratified analyses showed that the association in the total sample between block 2 and obesity may be mainly driven by female subjects. Conclusion: Intronic variants of the LRP5 gene are markedly associated with obesity. We hypothesise that such an association may be due to the role of LRP5 in the WNT signalling pathway or lipid metabolism. Further functional studies are needed to elucidate the exact molecular mechanism underlying our finding.

Genome-wide association study identifies ALDH7A1 as a novel susceptibility gene for osteoporosis.

Osteoporosis is a major public health problem. It is mainly characterized by low bone mineral density (BMD) and/or low-trauma osteoporotic fractures (OF), both of which have strong genetic determination. The specific genes influencing these phenotypic traits, however, are largely unknown. Using the Affymetrix 500K array set, we performed a case-control genome-wide association study (GWAS) in 700 elderly Chinese Han subjects (350 with hip OF and 350 healthy matched controls). A follow-up replication study was conducted to validate our major GWAS findings in an independent Chinese sample containing 390 cases with hip OF and 516 controls. We found that a SNP, rs13182402 within the ALDH7A1 gene on chromosome 5q31, was strongly associated with OF with evidence combined GWAS and replication studies (P = 2.08×10−9, odds ratio = 2.25). In order to explore the target risk factors and potential mechanism underlying hip OF risk, we further examined this candidate SNPs relevance to hip BMD both in Chinese and Caucasian populations involving 9,962 additional subjects. This SNP was confirmed as consistently associated with hip BMD even across ethnic boundaries, in both Chinese and Caucasians (combined P = 6.39×10−6), further attesting to its potential effect on osteoporosis. ALDH7A1 degrades and detoxifies acetaldehyde, which inhibits osteoblast proliferation and results in decreased bone formation. Our findings may provide new insights into the pathogenesis of osteoporosis.

Genome-wide association study suggested copy number variation may be associated with body mass index in the Chinese population

Obesity is a major public health problem characterized with high body mass index (BMI). Copy number variations (CNVs) have been identified to be associated with complex human diseases. The effect of CNVs on obesity is unknown. In this study, we explored the association of CNVs with BMI in 597 Chinese Han subjects using Affymetrix GeneChip Human Mapping 500K Array Set. We found that one CNV at 10q11.22 (from 46.36 Mb to 46.56 Mb) was associated with BMI (the raw P=0.011). The CNV contributed 1.6% of BMI variation, and it covered one important obesity gene—pancreatic polypeptide receptor 1(PPYR1). It was reported that PPYR1 was a key regulator of energy homeostasis. Our findings suggested that CNV might be potentially important for the BMI variation. In addition, our study suggested that CNV might be used as a genetic marker to locate genes associated with BMI in Chinese population.

Genome-wide Association and Replication Studies Identified TRHR as an Important Gene for Lean Body Mass

Low lean body mass (LBM) is related to a series of health problems, such as osteoporotic fracture and sarcopenia. Here we report a genome-wide association (GWA) study on LBM variation, by using Affymetrix 500K single-nucleotide polymorphism (SNP) arrays. In the GWA scan, we tested 379,319 eligible SNPs in 1,000 unrelated US whites and found that two SNPs, rs16892496 (p = 7.55 x 10(-8)) and rs7832552 (p = 7.58 x 10(-8)), within the thyrotropin-releasing hormone receptor (TRHR) gene were significantly associated with LBM. Subjects carrying unfavorable genotypes at rs16892496 and rs7832552 had, on average, 2.70 and 2.55 kg lower LBM, respectively, compared to those with alternative genotypes. We replicated the significant associations in three independent samples: (1) 1488 unrelated US whites, (2) 2955 Chinese unrelated subjects, and (3) 593 nuclear families comprising 1972 US whites. Meta-analyses of the GWA scan and the replication studies yielded p values of 5.53 x 10(-9) for rs16892496 and 3.88 x 10(-10) for rs7832552. In addition, we found significant interactions between rs16892496 and polymorphisms of several other genes involved in the hypothalamic-pituitary-thyroid and the growth hormone-insulin-like growth factor-I axes. Results of this study, together with the functional relevance of TRHR in muscle metabolism, support the TRHR gene as an important gene for LBM variation.

Pathway-based genome-wide association analysis identified the importance of regulation-of-autophagy pathway for ultradistal radius BMD.

Wrist fracture is not only one of the most common osteoporotic fractures but also a predictor of future fractures at other sites. Wrist bone mineral density (BMD) is an important determinant of wrist fracture risk, with high heritability. Specific genes underlying wrist BMD variation are largely unknown. Most published genome‐wide association studies (GWASs) have focused only on a few top‐ranking single‐nucleotide polymorphisms (SNPs)/genes and considered each of the identified SNPs/genes independently. To identify biologic pathways important to wrist BMD variation, we used a novel pathway‐based analysis approach in our GWAS of wrist ultradistal radius (UD) BMD, examining approximately 500,000 SNPs genome‐wide from 984 unrelated whites. A total of 963 biologic pathways/gene sets were analyzed. We identified the regulation‐of‐autophagy (ROA) pathway that achieved the most significant result (p = .005, qfdr = 0.043, pfwer = 0.016) for association with UD BMD. The ROA pathway also showed significant association with arm BMD in the Framingham Heart Study sample containing 2187 subjects, which further confirmed our findings in the discovery cohort. Earlier studies indicated that during endochondral ossification, autophagy occurs prior to apoptosis of hypertrophic chondrocytes, and it also has been shown that some genes in the ROA pathway (e.g., INFG) may play important roles in osteoblastogenesis or osteoclastogenesis. Our study supports the potential role of the ROA pathway in human wrist BMD variation and osteoporosis. Further functional evaluation of this pathway to determine the mechanism by which it regulates wrist BMD should be pursued to provide new insights into the pathogenesis of wrist osteoporosis.

Molecular genetic studies of gene identification for osteoporosis: the 2009 update.

Osteoporosis is a complex human disease that results in increased susceptibility to fragility fractures. It can be phenotypically characterized using several traits, including bone mineral density, bone size, bone strength, and bone turnover markers. The identification of gene variants that contribute to osteoporosis phenotypes, or responses to therapy, can eventually help individualize the prognosis, treatment, and prevention of fractures and their adverse outcomes. Our previously published reviews have comprehensively summarized the progress of molecular genetic studies of gene identification for osteoporosis and have covered the data available to the end of September 2007. This review represents our continuing efforts to summarize the important and representative findings published between October 2007 and November 2009. The topics covered include genetic association and linkage studies in humans, transgenic and knockout mouse models, as well as gene-expression microarray and proteomics studies. Major results are tabulated for comparison and ease of reference. Comments are made on the notable findings and representative studies for their potential influence and implications on our present understanding of the genetics of osteoporosis.

Whole Genome Distribution and Ethnic Differentiation of Copy Number Variation in Caucasian and Asian Populations

Although copy number variation (CNV) has recently received much attention as a form of structure variation within the human genome, knowledge is still inadequate on fundamental CNV characteristics such as occurrence rate, genomic distribution and ethnic differentiation. In the present study, we used the Affymetrix GeneChip® Mapping 500K Array to discover and characterize CNVs in the human genome and to study ethnic differences of CNVs between Caucasians and Asians. Three thousand and nineteen CNVs, including 2381 CNVs in autosomes and 638 CNVs in X chromosome, from 985 Caucasian and 692 Asian individuals were identified, with a mean length of 296 kb. Among these CNVs, 190 had frequencies greater than 1% in at least one ethnic group, and 109 showed significant ethnic differences in frequencies (p<0.01). After merging overlapping CNVs, 1135 copy number variation regions (CNVRs), covering approximately 439 Mb (14.3%) of the human genome, were obtained. Our findings of ethnic differentiation of CNVs, along with the newly constructed CNV genomic map, extend our knowledge on the structural variation in the human genome and may furnish a basis for understanding the genomic differentiation of complex traits across ethnic groups.