- Research article
- Open Access
- Open Peer Review
Synopsis and meta-analysis of genetic association studies in osteoporosis for the focal adhesion family genes: the CUMAGAS-OSTEOporosis information system
BMC Medicinevolume 9, Article number: 9 (2011)
Focal adhesion (FA) family genes have been studied as candidate genes for osteoporosis, but the results of genetic association studies (GASs) are controversial. To clarify these data, a systematic assessment of GASs for FA genes in osteoporosis was conducted.
We developed Cumulative Meta-Analysis of GAS-OSTEOporosis (CUMAGAS-OSTEOporosis), a web-based information system that allows the retrieval, analysis and meta-analysis (for allele contrast, recessive, dominant, additive and codominant models) of data from GASs on osteoporosis with the capability of update. GASs were identified by searching the PubMed and HuGE PubLit databases.
Data from 72 studies involving 13 variants of 6 genes were analyzed and catalogued in CUMAGAS-OSTEOporosis. Twenty-two studies produced significant associations with osteoporosis risk under any genetic model. All studies were underpowered (<50%). In four studies, the controls deviated from the Hardy-Weinberg equilibrium. Eight variants were chosen for meta-analysis, and significance was shown for the variants collagen, type I, α1 (COL1A1) G2046T (all genetic models), COL1A1 G-1997T (allele contrast and dominant model) and integrin β-chain β3 (ITGB3) T176C (recessive and additive models). In COL1A1 G2046T, subgroup analysis has shown significant associations for Caucasians, adults, females, males and postmenopausal women. A differential magnitude of effect in large versus small studies (that is, indication of publication bias) was detected for the variant COL1A1 G2046T.
There is evidence of an implication of FA family genes in osteoporosis. CUMAGAS-OSTEOporosis could be a useful tool for current genomic epidemiology research in the field of osteoporosis.
Osteoporosis is a common skeletal disease characterized by generalized reduction in bone mineral density and microarchitectural deterioration of bone tissue, leading to impaired skeletal strength and increased susceptibility to fracture . Genetic factors have long been recognized as playing an important role in osteoporosis . Major efforts are currently underway to identify the specific genes and allelic variants predisposing patients to this disease. The identification of genes is achieved by conducting hypothesis-free, genome-wide association studies (GWASs) and candidate gene association studies (GASs) . Candidate genes are typically chosen on the basis of having biological effects on bone metabolism or bone cell activity and whether they contribute to the risk of osteoporosis .
The focal adhesion (FA) gene family has emerged as a logical candidate for osteoporosis. Focal adhesions are specialized structures at the cellular-extracellular matrix contact points, where bundles of actin filaments are anchored to transmembrane receptors of the integrin family through a multimolecular complex of junctional plaque proteins. Some of the constituents of FA genes participate in the structural link between membrane receptors and the actin cytoskeleton, while others are signaling molecules [4, 5]. Although there are a number of genes constituting the pathway, only a small number of variants of these genes have been studied in GASs of osteoporosis. More recently, haplotype-based approaches and genome-wide genotyping platforms have enabled more comprehensive capture of genetic variation in these genes . The most studied gene is the collagen, type I, α1 (COL1A1) gene, especially its variant G2046T. Other genes studied in the FA family include the genes insulin-like growth factor I (IGF-I), integrin β-chain β3 (ITGB3), α-actinin-3 (ACTN3), COL1A2 and type 1 insulin-like growth factor receptor (IGF-IR). However, the results of the GASs involving genes of the FA pathway and osteoporosis are controversial and inconclusive, possibly because of methodological limitations, including inadequate sample size, patient selection, ethnicity of the populations studied and lack of adjustments for confounders .
To explore the involvement of FA family gene polymorphisms in osteoporosis susceptibility, we systematically searched for all available GASs of FA family genes and osteoporosis (as a binary phenotype) and created the Cumulative Meta-Analysis of Genetic Association Studies-OSTEOporosis (CUMAGAS-OSTEOporosis) information system. Then we catalogued all retrieved articles and estimated the risk effects of all individually investigated variants. Finally, the available data were synthesized using meta-analysis techniques to increase the power for detecting significant results and to decrease the uncertainty of the estimated genetic risks .
CUMAGAS-OSTEOporosis is a web-based database and an information system for cumulative meta-analysis of GASs [; see also [10, 11]]. CUMAGAS-OSTEOporosis performs meta-analysis for all genetic models (allele contrast, dominant, recessive, additive and codominant) and provides data on various covariates. Currently, CUMAGAS-OSTEOporosis operates for binary phenotypes (that is, osteoporosis: yes or no), but our study group is expanding the system to analyze continuous phenotype (bone mineral density).
CUMAGAS-OSTEOporosis is a dynamic system, since it has the capacity of continuous updating. Authors of published and unpublished studies may contribute their data by entering their studies' data into a prespecified data entry form (CUMAGAS-FORM) . Furthermore, authors may correct previously stored data or notify for missed studies by contacting the CUMAGAS investigators (email@example.com).
Selection of Studies
All studies published before June 2010 were identified by conducting extended computer-based searches of the PubMed and HuGE PubLit databases. The search criteria in the PubMed database included a combination of the following terms: Focal adhesion, ACP1, ACTN3, ADRB1, AKT1, COL1A1, COL1A2, COMP, CREBBP, CSNK1D, CTNNA3, CTNNB1, DRD2, FGFR1, GRB2, GRB2, IGF-I, IGF-IR, ITGA1, ITGA2, ITGA2B, ITGB3, KDR, PARD6A, PDGFRA, PDGFRB, PIK3CA, PRKACB, PTPN1, SMAD2, SMAD3, SMAD4, SPP1, TCF7, TGFBR1, TGFBR2, VEGFA, VWF, osteoporosis, gene, polymorphism, allele and variant. The bibliographies in the articles that these searches identified were used to find further references. The HuGE PubLit database  was searched for the disease term "osteoporosis" and for the gene terms listed above.
The eligible studies fulfilled the following inclusion criteria: (1) inclusion of cases with clinically diagnosed osteoporosis and controls free of osteoporosis, (2) information provided on genotype frequency or risk estimates (only studies that reported a particular variant were considered, and missing data were not imputed), (3) use of DNA-based analytical methods for genotyping and (4) studies of humans. Studies investigating disease progression, severity, phenotype modification, response to treatment or survival were excluded from our study. Case reports, editorials, review articles and non-English-language articles were also excluded. Finally, family-based studies were excluded because of different design settings. Abstracts of retrieved studies were independently read by two investigators (CD and EZ) to assess their appropriateness for this study. Full-text articles of the studies were evaluated (by CD and EZ) according to the inclusion criteria. The results were compared, and disagreements were resolved by consensus.
Published GWASs of osteoporosis cited in the HuGE PubLit database  and the National Human Genome Research Institute Catalog of Published Genome-Wide Association Studies  were screened for variants of the FA gene family. Open access databases for GWASs  were also searched. The variants tested in candidate gene studies were examined regardless of whether they had been included or tagged by proxy variants in the genotyping platforms used in the GWASs of osteoporosis .
From each article, the following information was extracted: first author, year of publication, ethnicity of the study population, study design, demographics and number of cases and controls for each genotype and effect size. The frequencies of the alleles and the genotypic distributions were extracted or calculated for both the cases and the controls. The reference single nucleotide polymorphism identification numbers (rs numbers), the chromosomal gene position and the nucleotide base changes for all genetic variants were identified by performing extended searches of bioinformatics databases [12, 14–16].
Data analysis and synthesis
Prior to meta-analysis, the risk effect of gene variants for the allele contrast and the dominant models were evaluated separately for each study. All associations were indicated as odds ratios (ORs) with the corresponding 95% confidence intervals (95% CIs). In the meta-analysis, the heterogeneity between studies was tested using the Q-statistic [17, 18], and it was quantified with the I 2 metric . Heterogeneity was considered significant when P Q < 0.10 (P Q is the P-value for Q-statistic). The pooled OR was estimated using the random effects (RE) model [; see also [8, 20]]. The RE model was chosen because it is more conservative than the alternative fixed effects model, which does not consider heterogeneity. The differential magnitude of effect in large versus small studies (that is, indication of publication bias) was tested using a modified linear regression test for funnel plot asymmetry proposed by Harbord et al. . This effect was considered significant when the P-value for Harbord's test was P H<0.05. The meta-analysis consisted of the main (that is, overall) analysis, which included all available data, as well as subgroup analyses by ethnicity, gender, age, menopausal status and sensitivity analysis which examined the effect of excluding specific studies [8, 20].
The distribution of each variant in the control group was tested for the Hardy-Weinberg equilibrium (HWE) . HWE indicates possible genotyping errors and/or population stratification . Sensitivity analysis was carried out for the studies that deviated from HWE and the studies for which HWE could not be tested (that is, the pooled OR was calculated after excluding these studies). The power of each study for the allele contrast was calculated assuming a 20% alteration in effect size (that is, modest effect), a significance level of 0.05 and a disease allele frequency equal to the one of the study population . Analyses were performed using the CUMAGAS-OSTEOporosis database  and Compaq Visual Fortran 90 software (Compaq Computer Corporation, Houston, Texas, US) with the International Mathematics and Statistics Library (Visual Nuemerics Inc, Houston, Texas, US).
The literature review identified 169 titles in the PubMed and HuGELit databases that met the search criteria. The search in HuGE PubLit and the databases for GWASs traced articles already identified by PubMed. After abstract selection, 132 articles remained. When an article provided data for different populations, then each population was considered as a different study [23–27]. Thirty-nine articles consisting of 72 studies that investigated the association between genetic variants of the FA family genes and osteoporosis fulfilled the inclusion criteria. Figure 1 presents a flowchart of retrieved studies and studies that were excluded, with specification of the reasons for inclusion or exclusion (a list of the excluded studies is provided in Additional file 1). Overall, 6 genes and 13 distinct variants of these genes investigated in the 72 gene-disease association studies were identified. The studies were published between 1996 and 2010.
Studies' characteristics and association results
The characteristics of each study and the association results of variants are shown in Additional file 2. In GWASs, none of the variants of the FA gene family were reported as significant, nor were the variants examined in the meta-analyses captured by the genotyping platforms used in the GWASs [28, 29].
Studies were conducted in various populations of different racial descent: Sixty-six studies involved solely Caucasians, two studies recruited East Asians, three studies involved Turks and one study involved a population of Mexican origin. Twenty-three studies involved only postmenopausal women. Twelve studies provided data for men and women. One study involved children.
The distribution of genotypes in the control group deviated from HWE in 4 studies, and in 20 studies the HWE deviation could not be tested. In all studies, the statistical power for detecting a significant risk effect was lower than 50%. In total, 22 studies produced significant associations with osteoporosis risk under any genetic model. The significant associations concerned the variants COL1A1 G2046T, COL1A1 RsaI intron 5, COL1A2 PvuII and ITGB3 T176C.
In total, eight variants were investigated in two or more studies, and their results were meta-analyzed: COL1A1 G2046T (Sp1 SS/ss, rs1800012), COL1A1 G-1997T (rs1107946), COL1A1 -1663T ins/del (rs2412298), COL1A1 MspI 26 kb upstream, COL1A1 RsaI intron 5, COL1A1 MnII exon 52, IGF-I 192 bp CA and ITGB3 T176C (rs5920). Table 1 shows the meta-analysis results for the association between the different variants and the risk of developing osteoporosis. Significant results are shown for the variants COL1A1 G2046T, COL1A1 G-1997T and ITGB3 T176C. Additional file 3 shows the associations of the individual studies and the meta-analysis results for the dominant model of the variant COL1A1 G2046T.
In particular, there was a significant overall association for the allele contrast of the variant COL1A1 G2046T (OR, 1.65; 95% CI, 1.39-1.94) with the heterogeneity between studies being significant (I 2 = 69%; P Q < 0.01). In subgroup analysis, a significant association was shown for Caucasians (OR, 1.60; 95% CI, 1.36-1.88), adults (OR, 1.48; 95% CI, 1.29-1.70), females (OR, 1.37; 95% CI, 1.20-1.57), males (OR, 1.85; 95% CI, 1.09-3.14) and postmenopausal women (OR, 1.34; 95% CI, 1.14-1.58). The recessive, dominant, additive and codominant models produced similar results (with the exception of subgroup analysis for males for the recessive and additive models). The sensitivity analysis for HWE did not alter the pattern of results.
Significant results were also shown for the variants COL1A1 G-1997T (allele contrast and dominant model) and ITGB3 T176C (recessive and additive models). However, these results were based on a small number of studies, and therefore safe conclusions could not be drawn. In the overall meta-analyses for the allele contrast, a differential magnitude of effect in large versus small studies was detected only for the variant COL1A1 G2046T (P H = 0.01).
In this project, the currently available data from GASs on human FA family genes in osteoporosis were catalogued. Then the data were synthesized, and the involvement of FA gene variants in disease susceptibility was assessed comprehensively. The eligible GASs were catalogued in a publicly available web-based database and information system called CUMAGAS-OSTEOporosis . In total, eight variants were meta-analyzed. Significant results were shown for the variants COL1A1 G2046T, COL1A1 G-1997T and ITGB3 T176C. Type I collagen is the major protein constituent of bone and is therefore a strong and plausible candidate gene for osteoporosis. The COL1A1 G2046T polymorphism is a single base pair substitution (G → T) within the regulatory region of the COL1A1 gene . The COL1A1 G-1997T polymorphism has been identified in the proximal promoter of COL1A1 at position -1997 and is in linkage disequilibrium with the COL1A1 G2046T polymorphism . There is evidence that the promoter polymorphisms are functional and have effects on DNA binding and gene transcription , but it is unclear to what extent this polymorphism is associated with the biomechanical properties of bone or susceptibility to fracture. Taking into consideration the abundance of COL1A1 in the bones and the fact that these variants seem to play a major role in the function of the gene, the significant association found in our results may encourage intensive research in this area. The ITGB3 T176C polymorphism changes the conformational structure of the β3-subunit of integrin . The integrin β3-subunit is known to play a key role in the resorptive function of osteoclasts, as shown in vitro and in transgenic animal studies . Genetic variation in integrin β3 may influence bone remodeling and subsequent bone loss and risk of osteoporotic fractures.
CUMAGAS-OSTEOporosis is an evidence-based information system for systematically searching, reviewing and synthesizing data for GASs of osteoporosis, with the capacity for continuous updating. CUMAGAS is being expanded to additional complex diseases such chronic lymphocytic leukemia, peripheral arterial disease, hypertension and osteoarthritis [10, 9, 35, 36]. CUMAGAS also has the capacity to incorporate data from GWASs subject to their public availability.
The GWASs of osteoporosis have not highlighted a significant role for the FA family of genes. However, the commercial genotyping platforms  may underrepresent the variants of the FA pathway and, of course, the variants included in the meta-analysis. In addition, the variants identified to date from the GWAS approach explain only a fraction of the disease heritability, and therefore the potential role for the variants of FA pathway may not be excluded . Furthermore, the analyses of GWASs have missed associations of multilocus variants involved in pathways with pathophysiological relevance to disease mechanisms [3, 36].
Most of the published GASs are underpowered for detecting the minor contribution of common alleles. For example, a sample size of more than 10,000 patients is needed to achieve >80% power to detect a significant OR between 1.1 and 1.5 (modest effect) for a polymorphic locus in association with a complex disease . Meta-analysis is a tool that allows for analysis with the potential for higher power by pooling the results of multiple studies . However, there is no formal, established methodology for calculating the power of a meta-analysis. In addition, power analysis may not be applicable in meta-analysis, since it is a retrospective, all-inclusive synthesis of published studies [8, 38]. Nevertheless, type II errors are expected to be less common in a meta-analysis than in single studies [8, 39]. Currently, no single institution alone is able to provide a sufficient number of patients, and therefore the creation of large databases from consortia where researchers share their data are required. However, this need for data sharing has been pointed out by previous initiatives [40, 41].
Two problems in human genome epidemiology research are that negative studies are frequently unpublished and some studies do not provide extractable data, which lead to the well-known phenomenon of publication bias . However, negative results should also have a venue for publication. Moreover, the inclusion of "negative" and unpublished data in meta-analyses of GASs as a means of reducing publication bias is commonly suggested and is believed to help in pointing out genetic effects . Thus, the establishment of an electronic information system to aid in performing cumulative meta-analyses of (published and unpublished) GASs of osteoporosis and identifying significant genetic variants could be a valuable tool for ongoing research in the field. Furthermore, CUMAGAS-OSTEOporosis will support rapid progress in human genome epidemiology of osteoporosis by identifying valid and replicable associations and making the overall effect of each variant from published and unpublished studies rapidly available to researchers. Concerning the retrieval of unpublished studies, the authors of unpublished studies will be able to submit their data to the CUMAGAS-OSTEOporosis database. Since these data will not have undergone peer review, a sensitivity analysis may be carried out (that is, a meta-analysis that examines the effect of excluding these studies). Finally, CUMAGAS-OSTEOporosis is an open access system, and it may support efforts to prevent publication bias .
Publication bias was tested using the method proposed by Harbord et al. , which is a modification of Egger's test , and it is appropriate for small-study effects. The visual inspection of funnel plots was avoided, since their validity is questionable [46, 47]. However, the statistical tests used to evaluate studies for publication bias actually compare the differential magnitude of effects in large versus small studies [8, 48].
The significance of risk effects in the GASs was assessed using the OR metric for various genetic models (dominant, recessive, additive and codominant) by merging genotypes. These models are not independent, and there is no a priori biological justification for their choice. Consequently, the interpretation of the results can be problematic, especially when all genetic contrasts are significant as in the case of COL1A1 G2046T. In these cases, the introduction of the recently proposed generalized OR (ORG) as an overall genetic risk effect may be a remedy . The ORG is a single statistic that summarizes the magnitude and significance of the association without considering the hash of possible contrasts, and thus the interpretation of the results is straightforward . The meta-analysis methodology in CUMAGAS-OSTEOporosis is expanding to incorporate the ORG metric and the continuous phenotype of osteoporosis (that is, bone mineral density).
In the meta-analysis, various genetic contrasts in different populations were explored, since there is no biological justification for choosing a specific contrast . Thus, adjustment for multiple testing is not strictly required for such an exploratory study [8, 50]. In addition, the adjustment for multiple testing might not be necessary, since the data were synthesized with the objective of reducing the uncertainty of effect size, without a prespecified hypothesis [50–52]. Furthermore, an appropriate multiple test adjustment might be difficult because the investigated contrasts are not independent and there is no clear structure in the multiple tests [49, 52]. Finally, an adjustment for multiple comparisons (for example, Bonferroni's correction) concerns a general null hypothesis that there is no association in all genetic contrasts simultaneously, which is not likely [53, 54].
The phenotypic heterogeneity of osteoporosis in the included studies in our synopsis makes the effort to combine the findings of GASs meaningfully in the complex field of osteoporosis a difficult task. Thus, in the presence of the large heterogeneity, the results should be interpreted with caution. There may be fundamental differences in the nature of genetic susceptibility to osteoporosis between postmenopausal and premenopausal women and between the two sexes or even among patients of different ethnicities. Furthermore, environmental factors, particularly nutrients, have to be accurately evaluated together with complex genotyping, to weigh their importance in revealing functional variants with respect to specific genetic background [55, 56]. Our analysis used the available study-level allele and genotype distributions, precluding adjusted analysis for potential gene-gene and gene-environment interactions, for which raw genotype data would be required. Failure to account for interactions may have reduced the power of our analysis but is unlikely to have inflated the number of positive results.
There is evidence implicating the activity of the FA family of genes in osteoporosis. Future studies designed to investigate epistatic and gene-environment interactions may help in deriving more conclusive claims about the role of these genes in osteoporosis. The CUMAGAS-OSTEOporosis information system can be a useful resource for reviewing and interpreting the findings of the accumulating genomic epidemiologic research in osteoporosis.
Stewart TL, Ralston SH: Role of genetic factors in the pathogenesis of osteoporosis. J Endocrinol. 2000, 166: 235-245. 10.1677/joe.0.1660235.
Zintzaras E, Rodopoulou P, Koukoulis GN: BsmI, TaqI, ApaI and FokI polymorphisms in the vitamin D receptor (VDR) gene and the risk of osteoporosis: a meta-analysis. Dis Markers. 2006, 22: 317-326.
Kitsios GD, Zintzaras E: Genome-wide association studies: hypothesis-"free" or "engaged"?. Transl Res. 2009, 154: 161-164. 10.1016/j.trsl.2009.07.001.
Zaidel-Bar R, Cohen M, Addadi L, Geiger B: Hierarchical assembly of cell-matrix adhesion complexes. Biochem Soc Trans. 2004, 32: 416-420. 10.1042/BST0320416.
Petit V, Thiery JP: Focal adhesions: structure and dynamics. Biol Cell. 2000, 92: 477-494. 10.1016/S0248-4900(00)01101-1.
Li WF, Hou SX, Yu B, Li MM, Férec C, Chen JM: Genetics of osteoporosis: accelerating pace in gene identification and validation. Hum Genet. 2010, 127: 249-285. 10.1007/s00439-009-0773-z.
Mann V, Ralston SH: Meta-analysis of COL1A1 Sp1 polymorphism in relation to bone mineral density and osteoporotic fracture. Bone. 2003, 32: 711-717. 10.1016/S8756-3282(03)00087-5.
Zintzaras E, Lau J: Synthesis of genetic association studies for pertinent gene-disease associations requires appropriate methodological and statistical approaches. J Clin Epidemiol. 2008, 61: 634-645. 10.1016/j.jclinepi.2007.12.011.
Cumulative Meta-Analysis of Genetic Association Studies-OSTEOporosis (CUMAGAS-OSTEOporosis) database. [http://biomath.med.uth.gr/]
Zintzaras E, Zdoukopoulos N: A field synopsis and meta-analysis of genetic association studies in peripheral arterial disease: the CUMAGAS-PAD database. Am J Epidemiol. 2009, 170: 1-11. 10.1093/aje/kwp094.
Zintzaras E, Kitsios G: Field synopsis and synthesis of candidate-gene association studies in chronic lymphocytic leukemia: the CUMAGAS-CLL information system. Am J Epidemiol. 2009, 170: 671-678. 10.1093/aje/kwp201.
Yu W, Gwinn M, Clyne M, Yesupriya A, Khoury MJ: A navigator for human genome epidemiology. Nat Genet. 2008, 40: 124-125. 10.1038/ng0208-124.
Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA: Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci USA. 2009, 106: 9362-9367. 10.1073/pnas.0903103106.
Johnson AD, O'Donnell CJ: An open access database of genome-wide association results. BMC Med Genet. 2009, 10: 6-10.1186/1471-2350-10-6.
Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, de Bakker PI: SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008, 24: 2938-2939. 10.1093/bioinformatics/btn564.
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29: 308-311. 10.1093/nar/29.1.308.
Ioannidis JP, Trikalinos TA, Zintzaras E: Extreme between-study homogeneity in meta-analyses could offer useful insights. J Clin Epidemiol. 2006, 59: 1023-1032. 10.1016/j.jclinepi.2006.02.013.
Zintzaras E, Ioannidis JP: Heterogeneity testing in meta-analysis of genome searches. Genet Epidemiol. 2005, 28: 123-137. 10.1002/gepi.20048.
DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7: 177-188. 10.1016/0197-2456(86)90046-2.
Trikalinos TA, Salanti G, Zintzaras E, Ioannidis JP: Meta-analysis methods. Adv Genet. 2008, 60: 311-334. full_text.
Harbord RM, Egger M, Sterne JA: A modified test for small-study effects in meta-analyses of controlled trials with binary endpoints. Stat Med. 2006, 25: 3443-3457. 10.1002/sim.2380.
Zintzaras E: Variance estimation of allele-based odds ratio in the absence of Hardy-Weinberg equilibrium. Eur J Epidemiol. 2008, 23: 323-326. 10.1007/s10654-008-9242-6.
Ralston SH, Uitterlinden AG, Brandi ML, Balcells S, Langdahl BL, Lips P, Lorenc R, Obermayer-Pietsch B, Scollen S, Bustamante M, Husted LB, Carey AH, Diez-Perez A, Dunning AM, Falchetti A, Karczmarewicz E, Kruk M, van Leeuwen JP, van Meurs JB, Mangion J, McGuigan FE, Mellibovsky L, del Monte F, Pols HA, Reeve J, Reid DM, Renner W, Rivadeneira F, van Schoor NM, Sherlock RE, GENOMOS Investigators, et al: Large-scale evidence for the effect of the COLIA1 Sp1 polymorphism on osteoporosis outcomes: the GENOMOS study. PLoS Med. 2006, 3: e90-10.1371/journal.pmed.0030090.
Langdahl BL, Ralston SH, Grant SF, Eriksen EF: An Sp1 binding site polymorphism in the COLIA1 gene predicts osteoporotic fractures in both men and women. J Bone Miner Res. 1998, 13: 1384-1389. 10.1359/jbmr.19220.127.116.114.
Välimäki S, Tähtelä R, Kainulainen K, Laitinen K, Löyttyniemi E, Sulkava R, Välimäki M, Kontula K: Relation of collagen type I alpha 1 (COLIA 1) and vitamin D receptor genotypes to bone mass, turnover, and fractures in early postmenopausal women and to hip fractures in elderly people. Eur J Intern Med. 2001, 12: 48-56.
McGuigan FE, Armbrecht G, Smith R, Felsenberg D, Reid DM, Ralston SH: Predication of osteoporotic fractures by bone densitometry and COLIA1 genotyping: a prospective, population-based study in men and women. Osteoporos Int. 2001, 12: 91-96. 10.1007/s001980170139.
Husted LB, Harsløf T, Gonzalez-Bofill N, Schmitz A, Carstens M, Stenkjaer L, Langdahl BL: Haplotypes of promoter and intron 1 polymorphisms in the COLIA1 gene are associated with increased risk of osteoporosis. Calcif Tissue Int. 2009, 84: 85-96. 10.1007/s00223-008-9199-1.
Liu YZ, Pei YF, Liu JF, Yang F, Guo Y, Zhang L, Liu XG, Yan H, Wang L, Zhang YP, Levy S, Recker RR, Deng HW: Powerful bivariate genome-wide association analyses suggest the SOX6 gene influencing both obesity and osteoporosis phenotypes in males. PLoS One. 2009, 4: e6827-10.1371/journal.pone.0006827.
Kiel DP, Demissie S, Dupuis J, Lunetta KL, Murabito JM, Karasik D: Genome-wide association with bone mass and geometry in the Framingham Heart Study. BMC Med Genet. 2007, 8 (Suppl 1): S14-10.1186/1471-2350-8-S1-S14.
Keen RW, Woodford-Richens KL, Grant SF, Ralston SH, Lanchbury JS, Spector TD: Association of polymorphism at the type I collagen (COL1A1) locus with reduced bone mineral density, increased fracture risk, and increased collagen turnover. Arthritis Rheum. 1999, 42: 285-290. 10.1002/1529-0131(199902)42:2<285::AID-ANR10>3.0.CO;2-3.
Jin H, van't Hof RJ, Albagha OM, Ralston SH: Promoter and intron 1 polymorphisms of COL1A1 interact to regulate transcription and susceptibility to osteoporosis. Hum Mol Genet. 2009, 18: 2729-2738. 10.1093/hmg/ddp205.
García-Giralt N, Enjuanes A, Bustamante M, Mellibovsky L, Nogués X, Carreras R, Díez-Pérez A, Grinberg D, Balcells S: In vitro functional assay of alleles and haplotypes of two COL1A1-promoter SNPs. Bone. 2005, 36: 902-908.
Teitelbaum SL: Osteoclasts, integrins, and osteoporosis. J Bone Miner Metab. 2000, 18: 344-349. 10.1007/s007740070007.
Honda S, Honda Y, Bauer B, Ruan C, Kunicki TJ: The impact of three-dimensional structure on the expression of PlA alloantigens on human integrin β3. Blood. 1995, 86: 234-242.
Zintzaras E, Kitsios GD, Ziogas DC, Rodopoulou P, Karachalios T: Field synopsis and synthesis of genetic association studies in osteoarthritis: the CUMAGAS-OSTEO information system. Am J Epidemiol. 2010, 171: 851-858. 10.1093/aje/kwq016.
Kitsios GD, Zintzaras E: Synopsis and data synthesis of genetic association studies in hypertension for the adrenergic receptor family genes: the CUMAGAS-HYPERT database. Am J Hypertens. 2010, 23: 305-313. 10.1038/ajh.2009.251.
Goldstein DB: Common genetic variation and human traits. N Engl J Med. 2009, 360: 1696-1698. 10.1056/NEJMp0806284.
Muncer S, Taylor S, Craigie M: Power dressing and meta-analysis: incorporating power analysis into meta-analysis. J Adv Nurs. 2002, 38: 274-280. 10.1046/j.1365-2648.2002.02177.x.
Zintzaras E, Ioannidis JP: Meta-analysis for ranked discovery datasets: theoretical framework and empirical demonstration for microarrays. Comput Biol Chem. 2008, 32: 38-46. 10.1016/j.compbiolchem.2007.09.003.
GAIN Collaborative Research Group, Manolio TA, Rodriguez LL, Brooks L, Abecasis G, Association Study of Psoriasis, Ballinger D, Daly M, Donnelly P, Faraone SV, International Multi-Center ADHD Genetics Project, Frazer K, Gabriel S, Gejman P, Molecular Genetics of Schizophrenia Collaboration, Guttmacher A, Harris EL, Insel T, Kelsoe JR, Bipolar Genome Study, Lander E, McCowin N, Mailman MD, Nabel E, Ostell J, Pugh E, Sherry S, Sullivan PF, Major Depression Stage 1 Genomewide Association in Population-Based Samples Study, Thompson JF, Warram J, Genetics of Kidneys in Diabetes (GoKinD) Study, Wholley D, Milos PM, Collins FS: New models of collaboration in genome-wide assosiation studies: the Genetic Associatiof Information Network. Nat Genet. 2007, 39: 1045-1051. 10.1038/ng2127.
Zhai G, van Meurs JB, Livshits G, Meulenbelt I, Valdes AM, Soranzo N, Hart D, Zhang F, Kato BS, Richards JB, Williams FM, Inouye M, Kloppenburg M, Deloukas P, Slagboom E, Uitterlinden A, Spector TD: A genome-wide association study suggests that a locus within the ataxin 2 binding protein 1 gene is associated with hand osteoarthritis: the Treat-OA consortium. J Med Genet. 2009, 46: 614-616. 10.1136/jmg.2009.067314.
Calnan M, Smith GD, Sterne JA: The publication process itself was the major cause of publication bias in genetic epidemiology. J Clin Epidemiol. 2006, 59: 1312-1318. 10.1016/j.jclinepi.2006.05.002.
Ioannidis JP: Journals should publish all "null" results and should sparingly publish "positive" results. Cancer Epidemiol Biomarkers Prev. 2006, 15: 186-201. 10.1158/1055-9965.EPI-05-0921.
Evans T, Gulmezoglu M, Pang T: Registering clinical trials: an essential role for WHO. Lancet. 2004, 363: 1413-1414. 10.1016/S0140-6736(04)16136-9.
Egger M, Davey Smith G, Schneider M, Minder C: Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997, 315: 629-634.
Munafò MR, Clark TG, Flint J: Assessing publication bias in genetic association studies: evidence from a recent meta-analysis. Psychiatry Res. 2004, 129: 39-44.
Lau J, Ioannidis JP, Terrin N, Schmid CH, Olkin I: The case of the misleading funnel plot. BMJ. 2006, 333: 597-600. 10.1136/bmj.333.7568.597.
Ioannidis JP, Trikalinos TA, Ntzani EE, Contopoulos-Ioannidis DG: Genetic associations in large versus small studies: an empirical assessment. Lancet. 2003, 361: 567-571. 10.1016/S0140-6736(03)12516-0.
Zintzaras E: The generalized odds ratio as a measure of genetic risk effect in the analysis and meta-analysis of association studies. Stat Appl Genet Mol Biol. 2010, 9: Article 21
Rothman KJ: No adjustments are needed for multiple comparisons. Epidemiology. 1990, 1: 43-46. 10.1097/00001648-199001000-00010.
Savitz DA, Olshan AF: Multiple comparisons and related issues in the interpretation of epidemiologic data. Am J Epidemiol. 1995, 142: 904-908.
Bender R, Lange S: Adjusting for multiple testing: when and how?. J Clin Epidemiol. 2001, 54: 343-349. 10.1016/S0895-4356(00)00314-0.
Feise RJ: Do multiple outcome measures require p-value adjustment?. BMC Med Res Methodol. 2002, 2: 8-10.1186/1471-2288-2-8.
Perneger TV: What's wrong with Bonferroni adjustments. BMJ. 1998, 316: 1236-1238.
Ioannidis JP, Ng MY, Sham PC, Zintzaras E, Lewis CM, Deng HW, Econs MJ, Karasik D, Devoto M, Kammerer CM, Spector T, Andrew T, Cupples LA, Duncan EL, Foroud T, Kiel DP, Koller D, Langdahl B, Mitchell BD, Peacock M, Recker R, Shen H, Sol-Church K, Spotila LD, Uitterlinden AG, Wilson SG, Kung AW, Ralston SH: Meta-analysis of genome-wide scans provides evidence for sex- and site-specific regulation of bone mass. J Bone Miner Res. 2007, 22: 173-183. 10.1359/jbmr.060806.
Hsu YH, Xu X, Terwedow HA, Niu T, Hong X, Wu D, Wang L, Brain JD, Bouxsein ML, Cummings SR, Rosen CJ, Xu X: Large-scale genome-wide linkage analysis for loci linked to BMD at different skeletal sites in extreme selected sibships. J Bone Miner Res. 2007, 22: 184-194. 10.1359/jbmr.061015.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/9/9/prepub
This study was partially supported by the Hellenic Society for the Study of Bone Metabolism.
The authors declare that they have no competing interests.
TK conceived the study. EZ and TK were the principal investigators. EZ, CD and TK retrieved and assessed the articles. EZ, AK and TM developed the algorithms and the information system. EZ and PR performed the analysis. EZ, CD and TK drafted the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.