Replication of 6 Obesity Genes in a Meta-Analysis of Genome-Wide Association Studies from Diverse Ancestries

Obesity is a major public health problem with a significant genetic component. Multiple DNA polymorphisms/genes have been shown to be strongly associated with obesity, typically in populations of European descent. The aim of this study was to verify the extent to which 6 confirmed obesity genes (FTO, CTNNBL1, ADRB2, LEPR, PPARG and UCP2 genes) could be replicated in 8 different samples (n = 11,161) and to explore whether the same genes contribute to obesity-susceptibility in populations of different ancestries (five Caucasian, one Chinese, one African-American and one Hispanic population). GWAS-based data sets with 1000 G imputed variants were tested for association with obesity phenotypes individually in each population, and subsequently combined in a meta-analysis. Multiple variants at the FTO locus showed significant associations with BMI, fat mass (FM) and percentage of body fat (PBF) in meta-analysis. The strongest association was detected at rs7185735 (P-value = 1.01×10−7 for BMI, 1.80×10−6 for FM, and 5.29×10−4 for PBF). Variants at the CTNNBL1, LEPR and PPARG loci demonstrated nominal association with obesity phenotypes (meta-analysis P-values ranging from 1.15×10−3 to 4.94×10−2). There was no evidence of association with variants at ADRB2 and UCP2 genes. When stratified by sex and ethnicity, FTO variants showed sex-specific and ethnic-specific effects on obesity traits. Thus, it is likely that FTO has an important role in the sex- and ethnic-specific risk of obesity. Our data confirmed the role of FTO, CTNNBL1, LEPR and PPARG in obesity predisposition. These findings enhanced our knowledge of genetic associations between these genes and obesity-related phenotypes, and provided further justification for pursuing functional studies of these genes in the pathophysiology of obesity. Sex and ethnic differences in genetic susceptibility across populations of diverse ancestries may contribute to a more targeted prevention and customized treatment of obesity.


Introduction
Obesity, a state in which excess lipids accumulate in various body fat depots due to a chronic imbalance between energy intake and energy expenditure, is associated with many diseases such as type 2 diabetes mellitus, hypertension, coronary heart disease, and some cancers [1]. Although the true causes of the accelerating obesity epidemic have not been fully clarified, the prevalence of obesity continues to increase around the world. If current trends continue, it is estimated that by the year 2030 nearly 50% of adults in the United States will be clinically obese [2], and the world population will include 1.12 billion obese individuals [3]. Although the impact of environmental factors is likely to be significant, many studies have shown that body weight and obesity are strongly influenced by genetic factors, with heritability estimates often in excess of 50% particularly when derived from comparisons of identical and fraternal twins [4,5].
Most genetic studies of obesity have focused on body mass index (BMI) to evaluate whether a person is obese or not. However, BMI cannot distinguish fat mass from fat free mass [6,7]. Alternative measurements such as percentage of body fat (PBF) and fat mass (FM) are more homogeneous and reflect body fat content more accurately than BMI [8]. Few genetic studies of obesity, however, have utilized these latter measures of body fat content.
In the past few years, extensive efforts had focused on the detection of obesity genes. In this regard, relatively few genes identified through preliminary linkage scans or candidate gene approaches have been confirmed to be truly associated with obesity by replication studies and other methodologies [9]. Among potential reasons for the failure to replicate most candidate genes or linkage peaks, one can cite small sample sizes and limited numbers of DNA variants upon which these studies were based. More recently, under the common variant-common disease hypothesis, several genome-wide association studies (GWAS) and large-scale meta-analyses of multiple GWAS on obesity (mostly using BMI as phenotype) have been conducted. Thus far, a total of 58 genetic loci, all with small effect sizes, were found to be robustly associated with obesity-related traits in multiple populations ( [10]; see also: http://www.ncbi.nlm.nih.gov/gap/phegeni/). Among these genes, fat mass-and obesity-associated gene (FTO) stands out as the gene with the strongest significant association with obesity and it has been found to be associated with obesity in virtually all populations in which replication was attempted. Nevertheless, all genomic markers identified along with their putative genes have only been shown to have very small effects on BMI or the risk of obesity. Cumulatively, these genetic loci identified through GWAS account for less than 5% of the total heritability of BMI [11], leaving the vast majority of heritability yet unidentified. Presumably, additional variants/loci will eventually be detected with larger sample sizes combined with incorporation of rare variants, copy number variation markers, and other genomic and epigenomic features.
The leptin receptor (LEPR), a single-transmembrane-domain receptor of the cytokine receptor family [12], acts with an adipocyte-specific hormone leptin that regulates adipose-tissue mass through hypothalamic effects on satiety and energy metabolism. Peroxisome proliferator-activated receptor gamma(P-PARG) is a transcription factor expressed abundantly in adipose tissue which involved in adipogenesis by activating adipocyte differentiation and mediating the expression of fat cell-specific genes [13]. Adrenoceptor beta 2, surface (ADRB2) is a major lipolytic receptor in human fat cells which plays a key role in regulating energy balance through both thermogenesis and lipid mobilization from adipose tissues [14]. Uncoupling protein 2 (UCP2) is an inner mitochondrial membrane transporter which dissipates the proton gradient of inner mitochondrial membranes and releases stored energy as heat, thus has an important role in energy expenditure [15]. FTO gene encoding a nucleic acid demethylase plays a role in controlling feeding behavior and energy expenditure [16]. Catenin (cadherin-associated protein), blike 1 (CTNNBL1) encodes a protein homolog to b-catenin responsible for cell-to-cell adhesion and Wnt-signalling [17]. Replication of previous findings in well-designed and statistically powered studies and fine-mapping the causal variants are essential to elucidate the importance of these six genes on obesity. The present study represents a meta-analysis of samples with diverse ancestries in which we attempt to replicate associations between sequence variants in six important candidate genes of obesity and body fatness phenotypes and to explore whether the same genes contribute to obesity-susceptibility in populations of different ancestries in large study samples.

Phenotype Measurements
Several obesity-related phenotypes were measured. These include BMI, body composition (FM and PBF) measured by dual-energy X-ray absorptiometry (DXA) scanners (either Lunar Corp., Madison, WI, USA, or Hologic Inc., Bedford, MA, USA) following the manufacturer's protocols or by underwater weighing with corrections for pulmonary residual volume. Covariates, including sex, age, age 2 , weight, height, and scanner ID (in WHI-AA and WHI-HIS), were screened with a stepwise linear regression model. Raw measurements were adjusted for significant covariates. To correct for potential population stratification, principal components (PCs) were computed and the first five PCs (i.e., PC1-PC5, explained .95% variation in each population.) derived from genome-wide genotype data were also included as covariates. Residual scores for each phenotype were normalized by inverse quantile of the standard normal distribution to impose a normal distribution on phenotypes which were then subjected to further analysis.

Genotype Imputation
To combine data across different genotyping platforms and to achieve a higher genome coverage, extensive genotype imputation was performed. Briefly, haplotypes of individual GWAS were first phased by a Markov Chain Haplotyping algorithm (MACH) (http://genome.sph.umich.edu/wiki/MaCH) [22], and untyped genotypes were then imputed by Minimac (http://genome.sph. umich.edu/wiki/Minimac), based on phased haplotypes, using the freely available haplotype data of the 1000 Genomes Project (as of August, 2010) as reference panels. Reference samples included 283 individuals of European ancestry, 193 individuals of Asian ancestry, and 174 individuals of African ancestry. Imputation was performed by comparing the respective panel with the closest ancestry. For each GWAS, genotypes for untyped SNPs were imputed based on relevant population's haplotype reference panel. SNPs with imputation quality score (as assessed by r2.hat by Minimac) greater than 0.3 were retained in at least two studies, and with minor allele frequency (MAF) .0.05 in at least one study, were included for subsequent analyses. Prior to genotype imputation, strand orientations were checked and inconsistencies were resolved. Imputation results are summarized as an 'allele dosage' defined as the expected number of copies of the coded allele at that SNP (i.e., a fractional value between 0 and 2) for each genotype. In total, 2,954 genotyped or imputed autosomal SNPs spanning the FTO(1,275 SNPs), ADRB2(119 SNPs), CTNNBL1(483 SNPs), LEPR(550 SNPs), PPARG(435 SNPs) and UCP2 (92 SNPs) genes were analyzed.

Association Tests
Association was tested in each study between directly-typed or imputed SNPs and obesity phenotypes under an additive genetic model. For each study of unrelated subjects (i.e., OOS, KCOS, COS, WHI-AA, and WHI-HIS), association was examined by fitting a linear regression model using MACH2QTL (http://www. sph.umich.edu/csg/abecasis/MACH/download/), in which allele dosage was used as a predictor of phenotype. For family-based samples (i.e., QFS, FHS and IFS), a mixed linear model was used in which the effect of genetic relatedness within each pedigree was taken into consideration [23].

Meta-analysis
Conventional meta-analysis of individual studies was performed with weights proportional to the square root of the sample size using METAL software (http://www.sph.umich.edu/csg/ abecasis/metal/) [24], and Cochran's Q statistic and I 2 were calculated as measures of between-study heterogeneity. Randomeffect meta-analyses were performed particularly for SNPs with Q statistic P-value,0.05 or I 2 value.50%. Meta-analysis of effect size (regression slope) of the identified candidate SNP and the forest plot were performed using Review Manager (RevMan) (http://ims.cochrane.org/revman).
Regional association plots of the most significant SNPs were generated using LocusZoom [25]. Wright's F-statistics (F ST ) was calculated with the R package to assess genetic differences among populations using the differences in allele frequencies. We

Trans-ethnic Meta-analysis
Given the multiple ancestry groups included in the study, we further performed a trans-ethnic meta-analysis using a recent MANTRA (Meta-Analysis of Trans-ethnic Association studies) software [26] as an alternative and robust analytic approach that accounts for the trans-ethnic nature of the various cohorts. In contrast to the traditional meta-analysis, the trans-ethnic approach in MANTRA accommodates between-population heterogeneity of associated variants and their effect sizes by allowing for allelic effects to be most similar between the most closely related populations. MANTRA adopts a Bayesian framework and assumes that studies from closely related populations are more likely to share a common true effect size, and the true effect size is allowed to vary across different population clades. Default settings were used in MANTRA. Evidence in favor of association of the trait with the variant was assessed by means of a Bayes' factor (BF).

Results
Baseline characteristics of the subjects in each cohort are presented in Table 1. A total of 2,954 SNPs (77.2%-91.1% imputed SNPs) covering the six targeted genes were metaanalyzed. A total of 45 SNPs located in the FTO gene demonstrated significant associations with obesity phenotypes using a Bonfferoni corrected significance level of P,1.68610 25 . Linkage disequilibrium (LD) analysis revealed that these 45 SNPs were in almost complete LD (r 2 .0.80) and were located within the same LD block of approximately 50 kb ( Figure 1). Another 125 SNPs out of the 1,275 total SNPs spanning the FTO gene in the present study showed marginal associations with obesity phenotypes (meta-analysis P-values ranging from 1.71610 25 to 4.99610 22 ) (Results not shown). Table 2 shows the top 15 FTO SNPs associated with obesity traits. For most variants, little heterogeneity among cohorts was observed, except for rs9922708 for which moderate heterogeneity was observed (I 2 = 43% and Q P-value = 0.091). The direction of effects was consistent across studies except for three SNPs (rs9922708, rs17817449 and rs7206790). The trans-ethnic metaanalysis using MANTRA closely mirrored the results from the standard meta-analysis using METAL. The most significant association was observed for FTO rs7185735 (P-value = 1.01610 27 for BMI, 1.80610 26 for FM, and 5.29610 24 for PBF) ( Table 2), and allele-specific OR (95% CI) reached 1.20 (1.12-1.28) for BMI, 1.18(1.10-1.28) for FM and 1.13 (1.05-1.21) for PBF ( Figure 2). The frequency of minor allele G at this imputed SNP ranged from 0.11 (Chinese) to 0.47 (Caucasian) in different ethnic groups. The forest plot for this SNP indicated that its strongest association was observed in QFS (OR = 1.51 for BMI, 1.49 for FM and 1.37 for PBF). Carriers of the G allele had higher BMI, FM and PBF values ( Figure 2).
To identify potential sex-specific loci, we performed a series of meta-analyses stratified by sex. In female specific samples, 25 FTO SNPs were significantly associated with BMI and FM, all of which overlapped with the findings obtained with the combined sample. There was no evidence of significant associations in the male specific samples (Table 3). However, a nominally significant association was revealed at rs16952725 of FTO in males (Pvalue = 6.87610 24 for BMI and P-value = 8.55610 23 for FM), but not in the female specific sample (P-value = 0.55 for BMI and P-value = 0.43 for FM) ( Table 3). We estimated the statistical power at various locus heritabilities for the total sample, males and females ( Figure 3). Our total sample and female sample have . 90% power to detect a variant explaining .0.40% of heritability. However, the male sample only have 30% power to detect a variant explaining .0.40% of heritability.
Finally, we examined potential differences of SNP association with obesity phenotypes among the four ethnic groups (Chinese, Caucasians, African Americans and Hispanic-Americans). In Caucasians, we found that 35 FTO SNPs were significantly associated with BMI and FM (P,1.68610 25 ); 28 of these 35 SNPs overlapped with those derived from the analysis performed on the combined samples. The most significant association was observed for rs17817288 in FTO (P-value = 1.96610 27 for BMI, 2.76610 25 for FM, and 2.39610 23 for PBF). However, these significant FTO SNPs in Caucasians showed no or weak evidence of association in non-Caucasian populations (Table 4, Figure 1). Therefore, the association results of FTO obtained with the combined samples were mostly contributed by the Caucasian group. For example, Pvalue was 1.01610 27 for the most significant SNP rs7185735 in the combined samples for BMI, 4.90610 27 in Caucasians, 7.47610 22 in Chinese, 6.05610 21 in African-Americans and 4.45610 21 in Hispanic-Americans. In non-Caucasian populations, we only found marginally significant associations. Two FTO SNPs (rs7201444 and rs13335146) specific in Hispanic-Americans, 14 FTO SNPs and one LEPR SNP(rs9436299) specific in African Americans and one CTNNBL1 SNP (rs45500793) specific in Chinese showed marginally significant associations (P-value, 1610 23 ) (Table 4, Figure 1). Table 4 listed 2 SNPs specific in Hispanic-Americans, 15 SNPs specific in African Americans, 1 SNP specific in Chinese and top 15 SNPs specific in Caucasians. To assess genetic differences for these ethnic-specific SNPs (listed in Table 4) among the eight populations, F ST were calculated. Nine FTO variants (rs1861869, rs7186521, rs17817288, rs8044769, rs11075987, rs9935401, rs8051591, rs7193144, rs8043757) and CTNNBL1 SNP (rs45500793) showed little population differentiation (F ST = 0-0.05). F ST of other SNPs ranged from 0.05 to 0.31, suggesting moderate/great genetic differentiation. FTO rs7201444 showed very great differentiation with the frequency of minor allele A ranging from 0.001 (Chinese) to 0.365 (African-Americans).
In addition, we examined the association of previously reported SNPs of the six genes among our dataset ( rs7202116, rs9939609 and rs9930506) showed significant/nominal associations with BMI, FM and/or PBF. Two PPARG SNPs (rs1801282 and rs3856806) and 3 CTNNBL1 variants (rs6013029, rs16986921 and rs6020712) showed nominal associations with obesity traits.

Discussion
In this study we investigated the associations between 6 genes and obesity phenotypes, including direct measures of body fatness, using a dense set of variants genotyped or imputed in 11,161 subjects from four different ethnic backgrounds (Chinese, Caucasian, African-Americans and Hispanic-Americans). Our findings in this multiethnic population confirmed the importance of the FTO gene for obesity risk in humans. A total of 45 SNPs located in the FTO gene showed significant associations with the obesity phenotypes. The FTO protein affects demethylation of nuclear RNA in vitro [27], but whether the efficiency of this process depends on the FTO genotype or how this may be related to the observed effects on BMI or body fatness is not clear.
At present, the strongest associations between FTO SNPs and BMI belong to intronic SNPs, which might have a role in the regulation of FTO and/or nearby genes. It is critical to recognize, however, that associated SNPs are not necessarily causal SNPs underlying the association, and that the functional variants are still unknown. A recent study reported that rs7202116 G allele creates a CpG site along with other variants in perfect linkage disequilibrium with it [28], and these risk alleles may potentially have increased DNA methylation. Bioinformatic analyses also revealed that six FTO SNPs (rs11642015, rs17817497, rs3751812, rs17817964, rs62033408, and rs1421085) were located within candidate intronic regulatory elements and that two SNPs (rs11642015 and rs1421085) were predicted to have allele-specific binding affinities for different transcription factors [29]. Specifically, the T allele at rs11642015 binds Paired box protein 5 (PAX5) while the C allele at rs1421085 is predicted to have a substantially reduced binding affinity for Cut-like homeobox 1 (CUX1). Further investigation is warranted to identify potentially functional SNPs and the mechanisms by which various alleles at FTO influence the level of adiposity.
We detected moderate associations for multiple variants from the CTNNBL1, PPARG and LEPR genes, but did not detect association with any variants in ADRB2 and UCP2. Our data confirmed the associations with adiposity for three previously reported CTNNBL1 variants (rs6013029, rs16986921 and rs6020712) [17]. This finding is consistent with the results from Andreasen and co-workers, who in a study comprising 18,014 Danish participants found that the CTNNBL1 rs6013029 T-allele    and the rs6020846 G-allele confer an increased risk of developing obesity, especially morbid obesity [30]. Interestingly, a German population-based study (KORA) [31] failed to demonstrate association of CTNNBL1 variant rs6013029 with obesity or BMI. The discrepancy might be due to population-specific differences and/or insufficient statistical power. For the most frequently studied coding variants of UCP2 [32], ADRB2 [33], PPARG [34] and LEPR [35] genes, only rs1801282 and rs3856806 in PPARG gene showed nominal associations with BMI and/or PBF. There are substantial ethnic differences in the prevalence of excess body weight and obesity. Although differences in lifestyle are likely to account for some of the observed differences, genetic variability could also play a major role. Among the four ethnic groups (Chinese, Caucasian, African-Americans and Hispanic populations), we found that 35 SNPs in intron 1 of FTO gene were significantly associated with indicators of obesity in Caucasians. All of the SNPs overlapped with those reported in a previous metaanalysis study [36]. There was very limited or no evidence for associations between these SNPs and adiposity in the other ethnic groups. On the other hand, two FTO SNPs (rs7201444 and rs13335146) in intron 8 specific in Hispanic Americans, 11 FTO SNPs (rs8043738, rs8063472, rs8058460, rs1077129, rs8048396, rs1861869, rs17217144, rs2892469, rs1861868, rs7186521, rs7184874) in intron 1 specific in African Americans showed moderate evidence for associations that were not reported by previous studies. While several previous studies reported association between FTO SNPs and obesity-related phenotypes in Hispanic Americans [37][38][39], African Americans [29,40] and Asian populations [41][42][43][44], these studies were relatively small in sample sizes and showed mixed results. Wing et al. [39] genotyped 26 SNPs in intron 1 of FTO in 373 Hispanic Americans and observed associations between BMI and several SNPs that were previously reported to be associated with obesity (rs9939609, rs8050136, rs1121980, rs1421085, rs17817449 and rs3751812), and four other SNPs (rs8047395, rs10852521, rs8057044 and rs8044769). Song et al. [37] and Scuteri et al. [38] replicated associations of BMI with the well-known rs9939609 [37] and rs9930506 [38] in Hispanic Americans. However, their findings were not confirmed in our Hispanic Americans. The two Hispanic-American specific SNPs found in our study were in intron 8 of FTO. Interestingly, a detailed mapping (262 tag SNPs across the entire FTO gene) in individuals of African American descent demonstrated significant association for seven SNPs (rs708262, rs11076017, rs16952725, rs9932411, rs7191513, rs2689269, rs16952987) in intron 8 [40]. Our study together with Adeyemo et al. suggests the intron 8 of FTO is a second site, in addition to intron 1, playing a role in the association between FTO and indicators of obesity.
Previous studies in African Americans showed either very limited or no evidence for associations with the SNPs initially reported in European populations, such as rs9939609, rs1121980, rs17817449 and rs8050136 [29]. More comprehensive evaluations of the FTO locus identified several other variants, rs56137030 [29], rs3751812 [45,46], rs1108102 [39] and rs8057044 [47] that showed associations with BMI in African Americans. Among these variants, only rs1108102 (P-value = 6.15610 23 ) replicated in our African American sample. Among the 14 FTO SNPs specific to African Americans, rs1861868 was previously found to be the strongest SNP associated with BMI in Old Order Amish individuals with low physical activity (P,0.001) [48]. rs1861869  Table 4. P values of ethnic-specific SNPs identified in the meta-analysis.  and rs7186521 were found to be associated with weight and waist circumference in 843 unrelated individuals from an island population in the eastern Adriatic coast of Croatia [49]. The allele frequencies for most of the SNPs indicated moderate population differentiation, and thus could potentially lead to variable genetic impact on obesity phenotypes across populations. In summary, our study replicated the associations of FTO intron 1 variants with BMI, FM and PBF in Caucasians, and confirmed the importance of variants in intron 8 as well. Therefore, our results suggest that there are ethnic differences with regard to the effects of FTO on obesity and body fatness.
Marked sex differences in the prevalence of obesity between women and men suggest sex-specific genetic impacts on obesity risk, even though lifestyle and dietary factors are also likely to contribute to these differences. Furthermore, genetic factors may interact with levels of physical activity to modify obesity risk [50].
To determine if the effect of genetic variants differed by sex, we stratified by sex and found that 25 FTO SNPs were specifically associated with obesity in females only. This is consistent with a few studies that have shown sex-differences in the heritability of BMI and fat percentage [51,52]. In contrast, a nominally significant association was revealed at rs16952725 of the FTO gene in males only, but not in females. Our findings are further supported by a recent study in children which found that the FTO variant rs9939609 showed association with obesity and BMI among girls but not among boys [53]. A recent genome-wide association study modeling the effect of genotype-by-sex interaction on obesity phenotypes demonstrated sex-influenced associations between genetic variation at the LYPLAL1 locus and obesityrelated traits [54]. This study suggests that FTO may be a gene playing a role in the commonly observed sex-dimorphism in adiposity.
The majority of obesity loci have been discovered through GWAS in individuals of European descent, and, more recently, in Asians as well. Okada et al. performed a GWAS with BMI in 62,245 East Asian subjects and observed a significant association with FTO rs12149832 (P-value = 4.8610 222 ) [55]. However, our study failed to replicate the association of FTO rs12149832. Several other studies replicated FTO rs9939609 [41][42][43]56,57] and rs8050136 [44] with obesity and BMI in Chinese or Asians. Our study showed weak evidence of association with FM and PBF for rs9939609 and rs8050136 (P-values range from 9.53610 23 to 3.23610 22 ), but not with BMI. Lu and Loos reviewed the transferability of 36 GWAS identified BMI-associated SNPs between European and East Asian ancestry populations using both SNP-to-SNP and locus-wide comparisons [32]. SNP-to-SNP comparisons suggest that of the 32 SNPs found in European populations, 8 SNPs were non-polymorphic and another 6 showed no convincing evidence of association in East Asians and two of the four loci identified in East Asians showed some evidence of transferability to European populations [32]. However, locus-wide analyses suggested the more extensive transferability. For example, CDKAL1 rs9356744 showed genome-wide significant association with BMI in East Asians but not in European populations (Pvalue = 0.19). When examining the locus surrounding rs9356744, however, other SNPs that are not in LD with rs9356744 were actually associated with BMI in European populations (P-values, 10 23 ). Replication studies in other populations that scan the entire gene (rather than just one or a few SNPs) may lead to the discovery of other important genetic variants. Studying populations of different ancestries (especially those with smaller LD) could help fine-map disease or trait loci, eventually pinpointing the causal variant (s) and gene(s). Our study used a comprehensive approach by thoroughly examining six genes with dense SNP coverage to search for population-specific and/or shared obesity loci in multiple ancestry groups.
An important limitation of the current study is the relatively small numbers of male subjects (29% of total subjects) and subjects from the three non-Caucasian ethnic groups (Chinese, African and Hispanic populations, 25% of the total subjects). These low numbers may have contributed to the nominal associations detected in these sub-groups. Nevertheless, the present study provided additional evidence supporting the presence of ethnicand sex-differences for some prominent obesity variants and perhaps genes.
In conclusion, we have attempted to replicate previously reported associations between multiple DNA common variants relating to 6 obesity genes in populations representing four ethnicities. We were able to find confirmatory evidence for contributions of FTO, CTNNBL1, LEPR and PPARG related genomic variants to human variation in adiposity. In particular, FTO variants showed sex-specific and ethnic-specific associations with adiposity traits. between Illumina and Boston University. Funding support for the Framingham Whole Body DXA dataset was provided by NIH grants R01 AR/AG 41398. The datasets used for the analyses described in this manuscript were obtained from dbGaP at http://www.ncbi.nlm.nih.gov/ sites/entrez?db = gap through dbGaP accession phs000007.v14.p6.
Funding support for the Genetic Determinants of Bone Fragility was provided through the NIA Division of Geriatrics and Clinical Gerontology. Genetic Determinants of Bone Fragility is a genome-wide association studies funded as part of the NIA Division of Geriatrics and Clinical Gerontology. Assistance with phenotype harmonization and genotype cleaning, as well as with general study coordination, was provided by the NIA Division of Geriatrics and Clinical Gerontology and the NIA Division of Aging Biology. Support for the collection of datasets and samples were provided by the parent grant, Genetic Determinants of Bone Fragility (P01-AG018397). Funding support for the genotyping which was performed at the Johns Hopkins University Center for Inherited Diseases Research was provided by the NIH NIA. The datasets used for the analyses described in this manuscript were obtained from dbGaP at http://www. ncbi.nlm.nih.gov/sites/entrez?db = gap through dbGaP accession phs000138.v2.p1.
The WHI program is funded by the National Heart, Lung, and Blood Institute, National Institutes of Health, U.S. Department of Health and Human Services through contracts N01WH22110, 24152, 32100-2, 32105-6, 32108-9, 32111-13, 32115, 32118-32119, 32122, 42107-26, 42129-32, and 44221. This manuscript was not prepared in collaboration with investigators of the WHI, has not been reviewed and/or approved by the Women's Health Initiative (WHI), and does not necessarily reflect the opinions of the WHI investigators or the NHLBI. The datasets used for the analyses described in this manuscript were obtained from dbGaP at http://