Genome-wide association study identifies a single major locus contributing to survival into old age; the APOE locus revisited

By studying the loci that contribute to human longevity, we aim to identify mechanisms that contribute to healthy aging. To identify such loci, we performed a genome-wide association study (GWAS) comparing 403 unrelated nonagenarians from long-living families included in the Leiden Longevity Study (LLS) and 1670 younger population controls. The strongest candidate SNPs from this GWAS have been analyzed in a meta-analysis of nonagenarian cases from the Rotterdam Study, Leiden 85-plus study, and Danish 1905 cohort. Only one of the 62 prioritized SNPs from the GWAS analysis (P<1×10(-4) ) showed genome-wide significance with survival into old age in the meta-analysis of 4149 nonagenarian cases and 7582 younger controls [OR=0.71 (95% CI 0.65-0.77), P=3.39 × 10(-17) ]. This SNP, rs2075650, is located in TOMM40 at chromosome 19q13.32 close to the apolipoprotein E (APOE) gene. Although there was only moderate linkage disequilibrium between rs2075650 and the ApoE ε4 defining SNP rs429358, we could not find an APOE-independent effect of rs2075650 on longevity, either in cross-sectional or in longitudinal analyses. As expected, rs429358 associated with metabolic phenotypes in the offspring of the nonagenarian cases from the LLS and their partners. In addition, we observed a novel association between this locus and serum levels of IGF-1 in women (P=0.005). In conclusion, the major locus determining familial longevity up to high age as detected by GWAS was marked by rs2075650, which tags the deleterious effects of the ApoE ε4 allele. No other major longevity locus was found.     

Meta-analysis for genome-wide association study identifies multiple variants at the BIN1 locus associated with late-onset Alzheimer's disease

Recent GWAS studies focused on uncovering novel genetic loci related to AD have revealed associations with variants near CLU, CR1, PICALM and BIN1. In this study, we conducted a genome-wide association study in an independent set of 1034 cases and 1186 controls using the Illumina genotyping platforms. By coupling our data with available GWAS datasets from the ADNI and GenADA, we replicated the original associations in both PICALM (rs3851179) and CR1 (rs3818361). The PICALM variant seems to be non-significant after we adjusted for APOE e4 status. We further tested our top markers in 751 independent cases and 751 matched controls. Besides the markers close to the APOE locus, a marker (rs12989701) upstream of BIN1 locus was replicated and the combined analysis reached genome-wide significance level (p = 5E-08). We combined our data with the published Harold et al. study and meta-analysis with all available 6521 cases and 10360 controls at the BIN1 locus revealed two significant variants (rs12989701, p = 1.32E-10 and rs744373, p = 3.16E-10) in limited linkage disequilibrium (r² = 0.05) with each other. The independent contribution of both SNPs was supported by haplotype conditional analysis. We also conducted multivariate analysis in canonical pathways and identified a consistent signal in the downstream pathways targeted by Gleevec (P = 0.004 in Pfizer; P = 0.028 in ADNI and P = 0.04 in GenADA). We further tested variants in CLU, PICALM, BIN1 and CR1 for association with disease progression in 597 AD patients where longitudinal cognitive measures are sufficient. Both the PICALM and CLU variants showed nominal significant association with cognitive decline as measured by change in Clinical Dementia Rating-sum of boxes (CDR-SB) score from the baseline but did not pass multiple-test correction. Future experiments will help us better understand potential roles of these genetic loci in AD pathology.     

Comparison of the power between microsatellite and single-nucleotide polymorphism markers for linkage and linkage disequilibrium mapping of an electrophysiological phenotype

We performed linkage and linkage disequilibrium (LD) mapping analyses to compare the power between microsatellite and single nucleotide polymorphism (SNP) markers. Chromosome-wide analyses were performed for a quantitative electrophysiological phenotype, ttth1, on chromosome 7. Multipoint analysis of microsatellite markers using the variance component (VC) method showed the highest LOD score of 4.20 at 162 cM, near D7S509 (163.7 cM). Two-point analysis of SNPs using the VC method yielded the highest LOD score of 3.98 in the Illumina SNP data and 3.45 in the Affymetrix SNP data around 152-153 cM. In family-based single SNP and SNP haplotype LD analysis, we identified seven SNPs associated with ttth1. We searched for any potential candidate genes in the location of the seven SNPs. The SNPs rs1476640 and rs768055 are located in the FLJ40852 gene (a hypothetical protein), and SNP rs1859646 is located in the TAS2R5 gene (a taste receptor). The other four SNPs are not located in any known or annotated genes. We found the high density SNP scan to be superior to microsatellites because it is effective in downstream fine mapping due to a better defined linkage region. Our study proves the utility of high density SNP in genome-wide mapping studies.     

Genome-wide association analysis reveals putative Alzheimer's disease susceptibility loci in addition to APOE

Alzheimer's disease (AD) is a genetically complex and heterogeneous disorder. To date four genes have been established to either cause early-onset autosomal-dominant AD (APP, PSEN1, and PSEN2(1-4)) or to increase susceptibility for late-onset AD (APOE5). However, the heritability of late-onset AD is as high as 80%, (6) and much of the phenotypic variance remains unexplained to date. We performed a genome-wide association (GWA) analysis using 484,522 single-nucleotide polymorphisms (SNPs) on a large (1,376 samples from 410 families) sample of AD families of self-reported European descent. We identified five SNPs showing either significant or marginally significant genome-wide association with a multivariate phenotype combining affection status and onset age. One of these signals (p = 5.7 x 10(-14)) was elicited by SNP rs4420638 and probably reflects APOE-epsilon4, which maps 11 kb proximal (r2 = 0.78). The other four signals were tested in three additional independent AD family samples composed of nearly 2700 individuals from almost 900 families. Two of these SNPs showed significant association in the replication samples (combined p values 0.007 and 0.00002). The SNP (rs11159647, on chromosome 14q31) with the strongest association signal also showed evidence of association with the same allele in GWA data generated in an independent sample of approximately 1,400 AD cases and controls (p = 0.04). Although the precise identity of the underlying locus(i) remains elusive, our study provides compelling evidence for the existence of at least one previously undescribed AD gene that, like APOE-epsilon4, primarily acts as a modifier of onset age.     

Genome-wide association study of genetic determinants of LDL-c response to atorvastatin therapy: importance of Lp(a)

We carried out a genome-wide association study (GWAS) of LDL-c response to statin using data from participants in the Collaborative Atorvastatin Diabetes Study (CARDS; n = 1,156), the Anglo-Scandinavian Cardiac Outcomes Trial (ASCOT; n = 895), and the observational phase of ASCOT (n = 651), all of whom were prescribed atorvastatin 10 mg. Following genome-wide imputation, we combined data from the three studies in a meta-analysis. We found associations of LDL-c response to atorvastatin that reached genome-wide significance at rs10455872 (P = 6.13 × 10(-9)) within the LPA gene and at two single nucleotide polymorphisms (SNP) within the APOE region (rs445925; P = 2.22 × 10(-16) and rs4420638; P = 1.01 × 10(-11)) that are proxies for the ε2 and ε4 variants, respectively, in APOE. The novel association with the LPA SNP was replicated in the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) trial (P = 0.009). Using CARDS data, we further showed that atorvastatin therapy did not alter lipoprotein(a) [Lp(a)] and that Lp(a) levels accounted for all of the associations of SNPs in the LPA gene and the apparent LDL-c response levels. However, statin therapy had a similar effect in reducing cardiovascular disease (CVD) in patients in the top quartile for serum Lp(a) levels (HR = 0.60) compared with those in the lower three quartiles (HR = 0.66; P = 0.8 for interaction). The data emphasize that high Lp(a) levels affect the measurement of LDL-c and the clinical estimation of LDL-c response. Therefore, an apparently lower LDL-c response to statin therapy may indicate a need for measurement of Lp(a). However, statin therapy seems beneficial even in those with high Lp(a).     

A Multivariate Genome-Wide Association Analysis of 10 LDL Subfractions,and Their Response to Statin Treatment,in 1868 Caucasians

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log₁₀Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns - both in our original data and in a replication cohort - consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log₁₀Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.     

Influence of Genetic Variation on Plasma Protein Levels in Older Adults Using a Multi-Analyte Panel

Proteins, widely studied as potential biomarkers, play important roles in numerous physiological functions and diseases. Genetic variation may modulate corresponding protein levels and point to the role of these variants in disease pathophysiology. Effects of individual single nucleotide polymorphisms (SNPs) within a gene were analyzed for corresponding plasma protein levels using genome-wide association study (GWAS) genotype data and proteomic panel data with 132 quality-controlled analytes from 521 Caucasian participants in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort. Linear regression analysis detected 112 significant (Bonferroni threshold p = 2.44×10⁻⁵) associations between 27 analytes and 112 SNPs. 107 out of these 112 associations were tested in the Indiana Memory and Aging Study (IMAS) cohort for replication and 50 associations were replicated at uncorrected p<0.05 in the same direction of effect as those in the ADNI. We identified multiple novel associations including the association of rs7517126 with plasma complement factor H-related protein 1 (CFHR1) level at p<1.46×10⁻⁶⁰, accounting for 40 percent of total variation of the protein level. We serendipitously found the association of rs6677604 with the same protein at p<9.29×10⁻¹¹². Although these two SNPs were not in the strong linkage disequilibrium, 61 percent of total variation of CFHR1 was accounted for by rs6677604 without additional variation by rs7517126 when both SNPs were tested together. 78 other SNP-protein associations in the ADNI sample exceeded genome-wide significance (5×10⁻⁸). Our results confirmed previously identified gene-protein associations for interleukin-6 receptor, chemokine CC-4, angiotensin-converting enzyme, and angiotensinogen, although the direction of effect was reversed in some cases. This study is among the first analyses of gene-protein product relationships integrating multiplex-panel proteomics and targeted genes extracted from a GWAS array. With intensive searches taking place for proteomic biomarkers for many diseases, the role of genetic variation takes on new importance and should be considered in interpretation of proteomic results.     

SNPing away at complex diseases: analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease

There has been great interest in the prospects of using single-nucleotide polymorphisms (SNPs) in the search for complex disease genes, and several initiatives devoted to the identification and mapping of SNPs throughout the human genome are currently underway. However, actual data investigating the use of SNPs for identification of complex disease genes are scarce. To begin to look at issues surrounding the use of SNPs in complex disease studies, we have initiated a collaborative SNP mapping study around APOE, the well-established susceptibility gene for late-onset Alzheimer disease (AD). Sixty SNPs in a 1.5-Mb region surrounding APOE were genotyped in samples of unrelated cases of AD, in controls, and in families with AD. Standard tests were conducted to look for association of SNP alleles with AD, in cases and controls. We also used family-based association analyses, including recently developed methods to look for haplotype association. Evidence of association (P相似文献   

Association study of CRP gene polymorphisms with serum CRP level and cardiovascular risk in the NHLBI Family Heart Study

Recent epidemiological studies have indicated that baseline C-reactive protein (CRP) levels may have value in prediction of cardiovascular risk. Using six tag single-nucleotide polymorphisms (SNPs) selected from our complete list of SNPs on the CRP gene, we investigated the association of CRP genotypes with plasma CRP levels and cardiovascular risk in the National Heart, Lung, and Blood Institute (NHLBI) Family Heart Study cohort (1,296 Caucasians, 48.5% male, 54.7 +/- 12.8 yr old). There was a significant trend toward association of CRP haplotypes with CRP levels (P = 0.045). SNP analysis indicated a highly significant association of SNP -757 (rs3093059, P = 0.0004) and SNP -286 (rs3091244, P = 0.0065) and a borderline association of SNP -7180 (rs1341665, P = 0.06) with CRP levels. Neither CRP haplotypes nor individual SNP genotypes were associated with intima-media thickness of the common carotid or internal carotid artery or the bifurcation of the carotid arteries. These results indicated a strong impact of local SNPs of the CRP gene on plasma CRP levels, but there was no direct evidence that these genetically controlled CRP elevations by local CRP SNPs contributed to cardiovascular disease phenotypes.     

A genome-wide association study of BMI in American Indians

Numerous studies have been done to understand genetic contributors to BMI, but only a limited number of studies have been done in nonwhite groups such as American Indians. A genome-wide association study (GWAS) for BMI was therefore performed in Pima Indians. BMI measurements from a longitudinal study of 1,120 Pima Indians and 454,194 single-nucleotide polymorphisms (SNPs) from the 1 million Affymetrix SNP panel were used (35% of SNPs were excluded due to minor allele frequency <0.05). Data included BMI measured at multiple examinations collected from 1965 to 2004, as well as the maximum BMI at one of these visits. General and within-family tests were performed using a maximum-likelihood based mixed model procedure. No SNP reached a genome-wide significance level (estimated at P < 4.94 × 10(-7)). For repeated measures analyses, the strongest associations for general and within-family tests mapped to two different regions on chromosome 6 (rs9342220 (P = 1.39 × 10(-6)) and rs7758764 (P = 2.51 × 10(-6)), respectively). For maximum BMI, the strongest association for the general tests mapped to chromosome 4 (rs17612333; P = 1.98 × 10(-6)) and to chromosome 3 (rs11127958; P = 1.53 × 10(-6)) for the within-family tests. Further analysis is important because only a few of these regions have been previously implicated in a GWAS and genetic susceptibility may differ by ethnicity.     

Interleukin-15 and interleukin-15R alpha SNPs and associations with muscle, bone, and predictors of the metabolic syndrome

The aims of this study were to examine associations between two SNPs in the human IL-15 gene and three SNPs in the IL-15Ralpha gene with predictors of metabolic syndrome and phenotypes in muscle, strength, and bone at baseline and in response to resistance training (RT). Subjects were Caucasians who had not performed RT in the previous year and consisted of a strength cohort (n=748), volumetric cohort (n=722), and serum cohort (n=544). Subjects completed 12 weeks of unilateral RT of the non-dominant arm, using their dominant arm as an untrained control. ANCOVA analyses revealed gender-specific associations with: (1) IL-15 SNP (rs1589241) and cholesterol (p=0.04), LDL (p=0.02), the homeostasis model assessment (HOMA; p=0.03), and BMI (p=0.002); (2) IL-15 SNP (rs1057972) and the pre- to post-training absolute difference in 1RM strength (p=0.02), BMI (p=0.008), and fasting glucose (p=0.03); (3) IL-15Ralpha SNP (rs2296135) and baseline total bone volume (p=0.04) and the pre- to post-training absolute difference in isometric strength (p=0.01); and 4) IL-15Ralpha SNP (rs2228059) and serum triglycerides (p=0.04), baseline whole muscle volume (p=0.04), baseline cortical bone volume (p=0.04), and baseline muscle quality (p=0.04). All associations were consistent in showing a potential involvement of the IL-15 pathway with muscle and bone phenotypes and predictors of metabolic syndrome.     

ZBTB40基因与中国汉族人群和美国高加索人群骨密度的关联研究

最近的全基因组关联研究发现ZBTB40(zinc finger and BTB domain containing 40)基因是一个潜在的调节骨密度的新基因,为了可靠地验证该基因与骨密度的关联关系,采用精细定位关联研究来检测ZBTB40基因内及其附近的SNPs与骨密度的关系.首先在中国样本(1 627个不相关的汉族样本)和美国样本(2 286个不相关高加索样本)中对ZBTB40基因的50个SNPs进行基因分型,然后采用Plink软件检测ZBTB40基因与腰椎和髋部骨密度的关联关系.证实了以前报道的SNPs rs7524102与腰椎和髋部骨密度的关联关系.另外还发现5个SNPs(rs10917209、rs6426748、rs4433361、rs3856183和rs7550872)与腰椎和髋部骨密度相关联.其中最显著的SNP是位于ZBTB40基因上游区域的rs6426748,其与腰椎骨密度的关联P值为9.82×10-4,多重检验校正后仍然显著.连锁不平衡分析表明rs6426748与rs10917209,rs7524102和rs7550872呈现高度连锁关系,位于长度为27 kb的单体型块内,因此ZBTB40基因的功能致病位点可能是rs6426748或该单体型块内某一未分型的SNP.     

GWAS-Linked Loci and Neuroimaging Measures in Alzheimer’s Disease

Recently, 19 susceptibility loci for Alzheimer’s disease (AD) had been identified through AD genome-wide association studies (GWAS) meta-analysis. However, how they influence the pathogenesis of AD still remains largely unknown. We studied those loci with six MRI measures, abnormal glucose metabolism, and β-amyloid (Aβ) deposition on neuroimaging in a large cohort from Alzheimer’s Disease Neuroimaging Initiative (ADNI) database in order to provide clues of the mechanisms through which these genetic variants might be acting. As a result, single nucleotide polymorphisms (SNPs) at rs983392 within MS4A6A and rs11218343 within SOLR1 were both associated with the percentage of increase in the volume of left inferior temporal regions in the follow-up study. Meanwhile, rs11218343 at SORL1 and rs6733839 at BIN1 was associated with rate of volume change of left parahippocampal and right inferior parietal, respectively. Moreover, rs6656401 at CR1 and rs983392 at MS4A6A were both associated with smaller volume of right middle temporal at baseline. However, in addition to the APOE locus, we did not detect any influence on glucose metabolism and Aβ deposition. APOE ε4 allele was associated with almost all measures. Altogether, five loci (rs6656401 at CR1, rs983392within MS4A6A, rs11218343 at SORL1, rs6733839 at BIN1, and APOE ε4) have been detected to be associated with one or a few established AD-related neuroimaging measures.     

Identification of six novel susceptibility loci for dyslipidemia using longitudinal exome-wide association studies in a Japanese population

Recent genome-wide association studies have identified various dyslipidemia-related genetic variants. However, most studies were conducted in a cross-sectional manner. We thus performed longitudinal exome-wide association studies of dyslipidemia in a Japanese population. We used ~244,000 genetic variants and clinical data of 6022 Japanese individuals who had undergone annual health checkups for several years. After quality control, the association of dyslipidemia-related phenotypes with 24,691 single nucleotide polymorphisms (SNPs) was tested using the generalized estimating equation model. In total, 82 SNPs were significantly (P < 2.03 × 10^?6) associated with dyslipidemia phenotypes. Of these SNPs, four (rs74416240 of TCHP, rs925368 of GIT2, rs7969300 of ATXN2, and rs12231744 of NAA25) and two (rs34902660 of SLC17A3 and rs1042127 of CDSN) were identified as novel genetic determinants of hypo-HDL- and hyper-LDL-cholesterolemia, respectively. A replication study using the cross-sectional data of 8310 Japanese individuals showed the association of the six identified SNPs with dyslipidemia-related traits.     

GAB2 alleles modify Alzheimer's risk in APOE epsilon4 carriers

The apolipoprotein E (APOE) epsilon4 allele is the best established genetic risk factor for late-onset Alzheimer's disease (LOAD). We conducted genome-wide surveys of 502,627 single-nucleotide polymorphisms (SNPs) to characterize and confirm other LOAD susceptibility genes. In epsilon4 carriers from neuropathologically verified discovery, neuropathologically verified replication, and clinically characterized replication cohorts of 1411 cases and controls, LOAD was associated with six SNPs from the GRB-associated binding protein 2 (GAB2) gene and a common haplotype encompassing the entire GAB2 gene. SNP rs2373115 (p = 9 x 10(-11)) was associated with an odds ratio of 4.06 (confidence interval 2.81-14.69), which interacts with APOE epsilon4 to further modify risk. GAB2 was overexpressed in pathologically vulnerable neurons; the Gab2 protein was detected in neurons, tangle-bearing neurons, and dystrophic neuritis; and interference with GAB2 gene expression increased tau phosphorylation. Our findings suggest that GAB2 modifies LOAD risk in APOE epsilon4 carriers and influences Alzheimer's neuropathology.     

Relationship between five GLUT1 gene single nucleotide polymorphisms and diabetic nephropathy: a systematic review and meta-analysis

So far, case-control studies on the association between glucose transporter 1 (GLUT1) gene single nucleotide polymorphisms (SNPs) and diabetic nephropathy (DN) have generated considerable controversy. To clarify the linkage of GLUT1 SNPs on the risk of DN, a systematic review and meta-analysis was performed. A comprehensive literature search of electronic databases was conducted to obtain relative studies. Nine case-control studies were included. Significant differences were found between XbaI SNP (rs841853) and increased risk of DN in all genetic models. Subgroup analyses for Caucasians population and DN from both type 1 and type 2 diabetes also revealed positive results. For Enh2-1 SNP (rs841847), Enh2-2 SNP (rs841848) and HaeIII SNP (rs1385129), obvious linkages were demonstrated in recessive model. However, analysis for the association between HpyCH4V SNP (rs710218) and the susceptibility of DN showed no significance. Likewise, negative outcome was also found in the assessment for the influence of XbaI or Enh2-2 SNP on the pathogenesis progress of DN. The evidence currently available shows that XbaI, Enh2 and HaeIII SNPs, but not HpyCH4V SNP, in GLUT1 gene may be genetic susceptibility to DN. However, data does not support the association between either XbaI or Enh2-2 SNP and the severity of DN.     

Genetic Variants in the Fat and Obesity Associated (FTO) Gene and Risk of Alzheimer's Disease

Background

Recent studies showed that polymorphisms in the Fat and Obesity-Associated (FTO) gene have robust effects on obesity, obesity-related traits and endophenotypes associated with Alzheimer''s disease (AD).

Methods

We used 1,877 Caucasian cases and controls from the NIA-LOAD study and 1,093 Caribbean Hispanics to further explore the association of FTO with AD. Using logistic regression, we assessed 42 SNPs in introns 1 and 2, the region previously reported to be associated with AD endophenotypes, which had been derived by genome-wide screenings. In addition, we performed gene expression analyses of neuropathologically confirmed AD cases and controls of two independent datasets (19 AD cases, 10 controls; 176 AD cases, 188 controls) using within- and between-group factors ANOVA of log₁₀ transformed rank invariant normalized expression data.

Results

In the NIALOAD study, one SNP was significantly associated with AD and three additional markers were close to significance (rs6499640, rs10852521, rs16945088, rs8044769, FDR p-value: 0.05<p<0.09). Two of the SNPs are in strong LD (D′>0.9) with the previously reported SNPs. In the Caribbean Hispanic dataset, we identified three SNPs (rs17219084, rs11075996, rs11075997, FDR p-value: 0.009<p<0.01) that were associated with AD. These results were confirmed by haplotype analyses and in a metaanalysis in which we included the ADNI dataset. FTO had a significantly lower expresssion in AD cases compared to controls in two independent datasets derived from human cortex and amygdala tissue, respectively (p = 2.18×10−5 and p<0.0001).

Conclusions

Our data support the notion that genetic variation in Introns 1 and 2 of the FTO gene may contribute to AD risk.     

Identification of PLCL1 gene for hip bone size variation in females in a genome-wide association study

Osteoporosis, the most prevalent metabolic bone disease among older people, increases risk for low trauma hip fractures (HF) that are associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of the key measurable risk factors for HF. Although hip BS is highly genetically determined, genetic factors underlying the trait are still poorly defined. Here, we performed the first genome-wide association study (GWAS) of hip BS interrogating approximately 380,000 SNPs on the Affymetrix platform in 1,000 homogeneous unrelated Caucasian subjects, including 501 females and 499 males. We identified a gene, PLCL1 (phospholipase c-like 1), that had four SNPs associated with hip BS at, or approaching, a genome-wide significance level in our female subjects; the most significant SNP, rs7595412, achieved a p value of 3.72x10(-7). The gene's importance to hip BS was replicated using the Illumina genotyping platform in an independent UK cohort containing 1,216 Caucasian females. Two SNPs of the PLCL1 gene, rs892515 and rs9789480, surrounded by the four SNPs identified in our GWAS, achieved p values of 8.62x10(-3) and 2.44x10(-3), respectively, for association with hip BS. Imputation analyses on our GWAS and the UK samples further confirmed the replication signals; eight SNPs of the gene achieved combined imputed p values<10(-5) in the two samples. The PLCL1 gene's relevance to HF was also observed in a Chinese sample containing 403 females, including 266 with HF and 177 control subjects. A SNP of the PLCL1 gene, rs3771362 that is only approximately 0.6 kb apart from the most significant SNP detected in our GWAS (rs7595412), achieved a p value of 7.66x10(-3) (odds ratio = 0.26) for association with HF. Additional biological support for the role of PLCL1 in BS comes from previous demonstrations that the PLCL1 protein inhibits IP3 (inositol 1,4,5-trisphosphate)-mediated calcium signaling, an important pathway regulating mechanical sensing of bone cells. Our findings suggest that PLCL1 is a novel gene associated with variation in hip BS, and provide new insights into the pathogenesis of HF.     

Analysis of secondary phenotypes in multigroup association studies

Although case-control association studies have been widely used, they are insufficient for many complex diseases, such as Alzheimer's disease and breast cancer, since these diseases may have multiple subtypes with distinct morphologies and clinical implications. Many multigroup studies, such as the Alzheimer's Disease Neuroimaging Initiative (ADNI), have been undertaken by recruiting subjects based on their multiclass primary disease status, while extensive secondary outcomes have been collected. The aim of this paper is to develop a general regression framework for the analysis of secondary phenotypes collected in multigroup association studies. Our regression framework is built on a conditional model for the secondary outcome given the multigroup status and covariates and its relationship with the population regression of interest of the secondary outcome given the covariates. Then, we develop generalized estimation equations to estimate the parameters of interest. We use both simulations and a large-scale imaging genetic data analysis from the ADNI to evaluate the effect of the multigroup sampling scheme on standard genome-wide association analyses based on linear regression methods, while comparing it with our statistical methods that appropriately adjust for the multigroup sampling scheme. Data used in preparation of this article were obtained from the ADNI database.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.