Beta-2 microglobulin (B2M) is a component of the major histocompatibility complex (MHC) class I molecule and has been studied as a biomarker of kidney function, cardiovascular diseases and mortality. Little is known about the genes influencing its levels directly or through glomerular filtration rate (GFR). We conducted a genome-wide association study of plasma B2M levels in 6738 European Americans from the Atherosclerosis Risk in Communities study to identify novel loci for B2M and assessed its association with known estimated GFR (eGFR) loci. We identified 2 genome-wide significant loci. One was in the human leukocyte antigen (HLA) region on chromosome 6 (lowest p value = 1.8 × 10(-23) for rs9264638). At this locus, 6 index SNPs accounted for 3.2 % of log(B2M) variance, and their association with B2M could largely be explained by imputed classical alleles of the MHC class I genes: HLA-A, HLA-B, or HLA-C. The index SNPs at this locus were not associated with eGFR based on serum creatinine (eGFRcr). The other locus of B2M was on chromosome 12 (rs3184504 at SH2B3, beta = 0.02, p value = 3.1 × 10(-8)), which was previously implicated as an eGFR locus. In conclusion, although B2M is known to be a component of MHC class I molecule, the association between HLA class I alleles and plasma B2M levels in a community-based population is novel. The identification of the two novel loci for B2M extends our understanding of its metabolism and informs its use as a kidney filtration biomarker.

Refractive error is the most common eye disorder worldwide and is a prominent cause of blindness. Myopia affects over 30% of Western populations and up to 80% of Asians. The CREAM consortium conducted genome-wide meta-analyses, including 37,382 individuals from 27 studies of European ancestry and 8,376 from 5 Asian cohorts. We identified 16 new loci for refractive error in individuals of European ancestry, of which 8 were shared with Asians. Combined analysis identified 8 additional associated loci. The new loci include candidate genes with functions in neurotransmission (GRIA4), ion transport (KCNQ5), retinoic acid metabolism (RDH5), extracellular matrix remodeling (LAMA2 and BMP2) and eye development (SIX6 and PRSS56). We also confirmed previously reported associations with GJD2 and RASGRF1. Risk score analysis using associated SNPs showed a tenfold increased risk of myopia for individuals carrying the highest genetic load. Our results, based on a large meta-analysis across independent multiancestry studies, considerably advance understanding of the mechanisms involved in refractive error and myopia.

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged≥25 y and women≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF)≤12 µg/L (cases) and iron replete controls (SF>100 µg/L in men, SF>50 µg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10(-7) for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P=7.0×10(-9), corrected P=0.012) was replicated within the VA samples (observed P=0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.

OBJECTIVE: We aimed to identify novel, functional, replicable and genome-wide significant risk regions specific for alcohol dependence using genome-wide association studies (GWASs). METHODS: A discovery sample (1409 European-American cases with alcohol dependence and 1518 European-American controls) and a replication sample (6438 European-Australian family subjects with 1645 alcohol dependent probands) underwent association analysis. Nineteen other cohorts with 11 different neuropsychiatric disorders served as contrast groups. Additional eight samples underwent expression quantitative locus (eQTL) analysis. RESULTS: A genome-wide significant risk gene region (NKAIN1-SERINC2) was identified in a meta-analysis of the discovery and replication samples. This region was enriched with 74 risk SNPs (unimputed); half of them had significant cis-acting regulatory effects. The distributions of -log(p) values for the SNP-disease associations or SNP-expression associations in this region were consistent throughout eight independent samples. Furthermore, imputing across the NKAIN1-SERINC2 region, we found that among all 795 SNPs in the discovery sample, 471 SNPs were nominally associated with alcohol dependence (1.7×10(-7)≤p≤0.047); 53 survived region- and cohort-wide correction for multiple testing; 92 SNPs were replicated in the replication sample (0.002≤p≤0.050). This region was neither significantly associated with alcohol dependence in African-Americans, nor with other non-alcoholism diseases. Finally, transcript expression of genes in NKAIN1-SERINC2 was significantly (p<3.4×10(-7)) associated with expression of numerous genes in the neurotransmitter systems or metabolic pathways previously associated with alcohol dependence. CONCLUSION: NKAIN1-SERINC2 may harbor a causal variant(s) for alcohol dependence. It may contribute to the disease risk by way of neurotransmitter systems or metabolic pathways.

Dietary intake of macronutrients (carbohydrate, protein, and fat) has been associated with risk of chronic conditions such as obesity and diabetes. Family studies have reported a moderate contribution of genetics to variation in macronutrient intake. In a genome-wide meta-analysis of a population-based discovery cohort (n = 33 533), rs838133 in FGF21 (19q13.33), rs197273 near TRAF family member-associated NF-kappa-B activator (TANK) (2p24.2), and rs10163409 in FTO (16q12.2) were among the top associations (P < 10(-5)) for percentage of total caloric intake from protein and carbohydrate. rs838133 was replicated in silico in an independent sample from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE) Nutrition Working Group (n = 38 360) and attained genome-wide significance in combined analysis (Pjoint = 7.9 × 10(-9)). A cytokine involved in cellular metabolism, FGF21 is a potential susceptibility gene for obesity and type 2 diabetes. Our results highlight the potential of genetic variation for determining dietary macronutrient intake.

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

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.

Hydrochlorothiazide (HCTZ) is one of the most widely prescribed antihypertensive medications. Although it is well known that HCTZ is associated with hyperglycemia and hypertriglyceridemia, the mechanisms underlying these adverse effects are not well understood. We performed a genome-wide association study and meta-analysis of the change in fasting plasma glucose and triglycerides in response to HCTZ from two different clinical trials: the Pharmacogenomic Evaluation of Antihypertensive Responses and the Genetic Epidemiology of Responses to Antihypertensive studies. Two single-nucleotide polymorphisms (rs12279250 and rs4319515 (r(2)=0.73)), located at 11p15.1 in the NELL1 gene, achieved genome-wide significance for association with change in fasting plasma triglycerides in African Americans, whereby each variant allele was associated with a 28 mg dl(-1) increase in the change in triglycerides. NELL1 encodes a cytoplasmic protein that contains epidermal growth factor-like repeats and has been shown to represses adipogenic differentiation. These findings may represent a novel mechanism underlying HCTZ-induced adverse metabolic effects.

BACKGROUND AND AIMS: Adiponectin is an adipose-secreted protein that has been linked to changes in insulin sensitivity, high-density lipoprotein cholesterol levels, and inflammatory patterns. Although fenofibrate therapy can raise adiponectin levels, treatment response is heterogeneous and heritable, suggesting a role for genetic mediators. This is the first genome-wide association study of fenofibrate effects on circulating adiponectin. METHODS AND RESULTS: Plasma adiponectin was measured in participants of the Genetics of Lipid Lowering Drugs and Diet Network (n = 793) before and after a 3-week daily treatment with 160 mg of fenofibrate. Associations between variants on the Affymetrix Genome-Wide Human SNP Array 6.0 and adiponectin were assessed using mixed linear models, adjusted for age, sex, site, and family. We observed a statistically significant (P = 5 × 10⁻⁸) association between rs2384207 in 12q24, a region previously linked to several metabolic traits, and the fenofibrate-induced change in circulating adiponectin. Additionally, our genome-wide analysis of baseline adiponectin levels replicated the previously reported association with CDH13 and suggested novel associations with markers near the PCK1, ZBP1, TMEM18, and SCUBE1 genes. The findings from the single marker tests were corroborated in gene-based analyses. Biological pathway analyses suggested a borderline significant association between the EGF receptor signaling pathway and baseline adiponectin levels. CONCLUSIONS: We present preliminary evidence linking several biologically relevant genetic variants to adiponectin levels at baseline and in response to fenofibrate therapy. Our findings provide support for fine-mapping of the 12q24 region to investigate the shared biological mechanisms underlying levels of circulating adiponectin and susceptibility to metabolic disease.

BACKGROUND: Beta-trace protein (BTP), measured in serum or plasma, has potential as a novel biomarker for kidney function. Little is known about the genes influencing BTP levels. METHODS: We conducted a genome-wide association study of log-transformed plasma BTP levels in 6720 European Americans (EAs) and replicated the significant associations in 1734 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study. RESULTS: We identified a genome-wide significant locus in EA upstream of Prostaglandin D2 synthase (PTGDS), the gene encoding BTP. Each copy of the A allele at rs57024841 was associated with 5% higher BTP levels (P = 1.2 × 10(-23)). The association at PTGDS was confirmed in AAs (6% higher BTP for each A allele at rs57024841, P = 1.9 × 10(-7)). The index single nucleotide polymorphisms (SNPs) in EAs and AAs explained ∼1.1% of the log(BTP) variance within each population and explained over 30% of the difference in log(BTP) levels between EAs and AAs. The index SNPs at the PTGDS locus in the two populations were not associated with the estimated glomerular filtration rate (eGFR) or the urine albumin creatinine ratio (P > 0.05). We further tested for the associations of BTP with 16 known loci of the eGFR in EA, and BTP was associated with 3 of 16 tested. CONCLUSIONS: The identification of a novel BTP-specific (non-renal related) locus and the confirmation of several genetic loci of the eGFR with BTP extend our understanding of the metabolism of BTP and inform its use as a kidney filtration biomarker.

Studying genetic determinants of intermediate phenotypes is a powerful tool to increase our understanding of genotype-phenotype correlations. Metabolic traits pertinent to the central nervous system (CNS) constitute a potentially informative target for genetic studies of intermediate phenotypes as their genetic underpinnings may elucidate etiological mechanisms. We therefore conducted a genome-wide association study (GWAS) of monoamine metabolite (MM) levels in cerebrospinal fluid (CSF) of 414 human subjects from the general population. In a linear model correcting for covariates, we identified one locus associated with MMs at a genome-wide significant level (standardized β=0.32, P=4.92 × 10(-8)), located 20 kb from SSTR1, a gene involved with brain signal transduction and glutamate receptor signaling. By subsequent whole-genome expression quantitative trait locus (eQTL) analysis, we provide evidence that this variant controls expression of PDE9A (β=0.21; P unadjusted=5.6 × 10(-7); P corrected=0.014), a gene previously implicated in monoaminergic transmission, major depressive disorder and antidepressant response. A post hoc analysis of loci significantly associated with psychiatric disorders suggested that genetic variation at CSMD1, a schizophrenia susceptibility locus, plays a role in the ratio between dopamine and serotonin metabolites in CSF. The presented DNA and mRNA analyses yielded genome-wide and suggestive associations in biologically plausible genes, two of which encode proteins involved with glutamate receptor functionality. These findings will hopefully contribute to an exploration of the functional impact of the highlighted genes on monoaminergic transmission and neuropsychiatric phenotypes.

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

Smoking is a major public health problem, but the genetic factors associated with smoking behaviors are not fully elucidated. Here, we have conducted an integrated genome-wide association study to identify common copy number polymorphisms (CNPs) and single nucleotide polymorphisms (SNPs) associated with the number of cigarettes smoked per day (CPD) in Japanese smokers (N = 17,158). Our analysis identified a common CNP with a strong effect on CPD (rs8102683; P=3.8 x 10(-42)) in the 19q13 region, encompassing the CYP2A6 locus. After adjustment for the associated CNP, we found an additional associated SNP (rs11878604; P=9.7 x 10(-30)) located 30 kb downstream of the CYP2A6 gene. Imputation of the CYP2A6 locus revealed that haplotypes underlying the CNP and the SNP corresponded to classical, functional alleles of CYP2A6 gene that regulate nicotine metabolism and explained 2% of the phenotypic variance of CPD (ANOVA F-test P=9.5 x 10(-52)). These haplotypes were also associated with smoking-related diseases, including lung cancer, chronic obstructive pulmonary disease and arteriosclerosis obliterans.

The study of the genetic regulation of metabolism in human serum samples can contribute to a better understanding of the intermediate biological steps that lead from polymorphism to disease. Here, we conducted a genome-wide association study (GWAS) to discover metabolic quantitative trait loci (mQTLs) utilizing samples from a study of prostate cancer in Swedish men, consisting of 402 individuals (214 cases and 188 controls) in a discovery set and 489 case-only samples in a replication set. A global nontargeted metabolite profiling approach was utilized resulting in the detection of 6,138 molecular features followed by targeted identification of associated metabolites. Seven replicating loci were identified (PYROXD2, FADS1, PON1, CYP4F2, UGT1A8, ACADL, and LIPC) with associated sequence variants contributing significantly to trait variance for one or more metabolites (P = 10(-13) -10(-91)). Regional mQTL enrichment analyses implicated two loci that included FADS1 and a novel locus near PDGFC. Biological pathway analysis implicated ACADM, ACADS, ACAD8, ACAD10, ACAD11, and ACOXL, reflecting significant enrichment of genes with acyl-CoA dehydrogenase activity. mQTL SNPs and mQTL-harboring genes were over-represented across GWASs conducted to date, suggesting that these data may have utility in tracing the molecular basis of some complex disease associations.

Bone and muscle, two major tissue types of musculoskeletal system, have strong genetic determination. Abnormality in bone and/or muscle may cause musculoskeletal diseases such as osteoporosis and sarcopenia. Bone size phenotypes (BSPs), such as hip bone size (HBS), appendicular bone size (ABS), are genetically correlated with body lean mass (mainly muscle mass). However, the specific genes shared by these phenotypes are largely unknown. In this study, we aimed to identify the specific genes with pleiotropic effects on BSPs and appendicular lean mass (ALM). We performed a bivariate genome-wide association study (GWAS) by analyzing ~690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) followed by a replication study in 2,286 unrelated US Caucasians (558 males and 1,728 females). We identified 14 interesting single nucleotide polymorphisms (SNPs) that may contribute to variation of both BSPs and ALM, with p values <10(-6) in discovery stage. Among them, the association of three SNPs (rs2507838, rs7116722, and rs11826261) in/near GLYAT (glycine-N-acyltransferase) gene was replicated in US Caucasians, with p values ranging from 1.89 × 10(-3) to 3.71 × 10(-4) for ALM-ABS, from 5.14 × 10(-3) to 1.11 × 10(-2) for ALM-HBS, respectively. Meta-analyses yielded stronger association signals for rs2507838, rs7116722, and rs11826261, with pooled p values of 1.68 × 10(-8), 7.94 × 10(-8), 6.80 × 10(-8) for ALB-ABS and 1.22 × 10(-4), 9.85 × 10(-5), 3.96 × 10(-4) for ALM-HBS, respectively. Haplotype allele ATA based on these three SNPs was also associated with ALM-HBS and ALM-ABS in both discovery and replication samples. Interestingly, GLYAT was previously found to be essential to glucose metabolism and energy metabolism, suggesting the gene's dual role in both bone development and muscle growth. Our findings, together with the prior biological evidence, suggest the importance of GLYAT gene in co-regulation of bone phenotypes and body lean mass.

Compressive strength index (CSI) is a newly established index for predicting hip fracture, the most serious consequence of osteoporosis. Appendicular lean mass (ALM), which influences skeletal strength of the lower limbs, is another trait associated with the risk of hip fracture. In this study, we performed a bivariate genome-wide association study (GWAS) to identify new candidate genes responsible for both CSI and ALM. In our discovery sample of 1627 unrelated Chinese subjects (802 males and 825 females), we scanned 909,509 SNPs using the Affymetrix Human Genome SNP 6.0 genotyping array. We successfully replicated our results in a sample of 2286 Caucasian subjects (558 males and 1728 females). The results indicated that five SNPs (rs174583, rs174577, rs174549, rs174548, rs7672337) in the FADS1, FADS2, and DCHS2 genes had significant bivariate associations with CSI and ALM in male subjects for both the GWAS discovery (with P<8.42×10(-6)) and the Caucasian sample (with P<0.07). We performed further replication analysis in a 2nd Caucasian sample with 501 Caucasian male subjects, using Affymetrix 500k arrays, and found that two of the above SNPs (rs174548 and rs174549, P=0.07) had bivariate associations with both CSI and ALM in males; the other 3 SNPs were not typed with the 500k array. The above findings suggest that the 3 genes, FADS1, FADS2, and DCHS2, containing these SNPs might play dual roles influencing both CSI and ALM in males. Our findings provide new insights into our understanding of the genetic basis of bone metabolism and the pathogenesis of osteoporosis.

BACKGROUND: Capecitabine, an oral 5-fluorouracil (5-FU) prodrug, is widely used in the treatment of breast, colorectal, and gastric cancers. To guide the selection of patients with potentially the greatest benefit of experiencing antitumor efficacy, or, alternatively, of developing toxicities, identifying genomic predictors of capecitabine sensitivity could permit its more informed use. METHODS: The objective of this study was to perform capecitabine sensitivity genome-wide association studies (GWAS) using 503 well genotyped human cell lines from individuals representing multiple different world populations. A meta-analysis that included all ethnic populations then enabled the identification of novel germline determinants (single nucleotide polymorphisms [SNPs]) of capecitabine susceptibility. RESULTS: First, an intrapopulation GWAS of Caucasian individuals identified reference SNP 4702484 (rs4702484) (within adenylate cyclase 2 [ADCY2]) at a level reaching genome-wide significance (P = 5.2 × 10(-8) ). This SNP is located upstream of the 5 methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene, and it is known that the enzyme for MTRR is involved in the methionine-folate biosynthesis and metabolism pathway, which is the primary target of 5-FU-related compounds, although the authors were unable to identify a direct relation between rs4702484 and MTRR expression in a tested subset of cells. In the meta-analysis, 4 SNPs comprised the top hits, which, again, included rs4702484 and 3 additional SNPs (rs8101143, rs576523, and rs361433) that approached genome-wide significance (P values from 1.9 × 10(-7) to 8.8 × 10(-7) ). The meta-analysis also identified 1 missense variant (rs11722476; serine to asparagine) within switch/sucrose nonfermentable-related, matrix-associated, actin-dependent regulator of chromatin (SMARCAD1), a novel gene for association with capecitabine/5-FU susceptibility. CONCLUSIONS: Toward the goal of individualizing cancer chemotherapy, the current study identified novel SNPs and genes associated with capecitabine sensitivity that are potentially informative and testable in any patient regardless of ethnicity.