OBJECTIVE: Methotrexate (MTX) is the drug of first choice for the treatment of rheumatoid arthritis (RA), but is effective only in around 60% of the patients. Identification of genetic markers to predict response is essential for effective treatment within a critical window period of 6 months after diagnosis, but have been hitherto elusive. In this study, we used genome-wide genotype data to identify the potential risk variants associated with MTX (poor)response in a north Indian RA cohort. MATERIALS AND METHODS: Genome-wide genotyping data for a total of 457 RA patients [297 good (DAS28-3≤3.2) and 160 poor (DAS28-3≥5.1) responders] on MTX monotherapy were tested for association using an additive model. Support vector machine and genome-wide pathway analysis were used to identify additional risk variants and pathways. All risk loci were imputed to fine-map the association signals and identify causal variant(s) of therapeutic/diagnostic relevance. RESULTS: Seven novel suggestive loci from genome-wide (P≤5×10(-5)) and three from support vector machine analysis were associated with MTX (poor)response. The associations of published candidate genes namely DHFR (P=0.014), FPGS (P=0.035), and TYMS (P=0.005) and purine and nucleotide metabolism pathways were reconfirmed. Imputation, followed by bioinformatic analysis indicated possible interaction between two reversely oriented overlapping genes namely ENOSF1 and TYMS at the post-transcriptional level. CONCLUSION: In this first ever genome-wide analysis on MTX treatment response in RA patients, 10 new risk loci were identified. These preliminary findings warrant replication in independent studies. Further, TYMS expression at the post-transcriptional level seems to be probably regulated through an antisense-RNA involving the 6-bp ins/del marker in the overlapping segment at 3'UTR of TYMS-ENOSF1, a finding with impending pharmacogenetic applications.

Fasting plasma glucose (FPG) has been recognized as an important indicator for the overall glycemic state preceding the onset of metabolic diseases. So far, most indentified genome-wide association loci for FPG were derived from populations with European ancestry, with a few exceptions. To extend a thorough catalog for FPG loci, we conducted meta-analyses of 13 genome-wide association studies in up to 24,740 nondiabetic subjects with East Asian ancestry. Follow-up replication analyses in up to an additional 21,345 participants identified three new FPG loci reaching genome-wide significance in or near PDK1-RAPGEF4, KANK1, and IGF1R. Our results could provide additional insight into the genetic variation implicated in fasting glucose regulation.

Circulating homocysteine levels (tHcy), a product of the folate one carbon metabolism pathway (FOCM) through the demethylation of methionine, are heritable and are associated with an increased risk of common diseases such as stroke, cardiovascular disease (CVD), cancer and dementia. The FOCM is the sole source of de novo methyl group synthesis, impacting many biological and epigenetic pathways. However, the genetic determinants of elevated tHcy (hyperhomocysteinemia), dysregulation of methionine metabolism and the underlying biological processes remain unclear. We conducted independent genome-wide association studies and a meta-analysis of methionine metabolism, characterized by post-methionine load test tHcy, in 2,710 participants from the Framingham Heart Study (FHS) and 2,100 participants from the Vitamin Intervention for Stroke Prevention (VISP) clinical trial, and then examined the association of the identified loci with incident stroke in FHS. Five genes in the FOCM pathway (GNMT [p = 1.60 × 10(-63)], CBS [p = 3.15 × 10(-26)], CPS1 [p = 9.10 × 10(-13)], ALDH1L1 [p = 7.3 × 10(-13)] and PSPH [p = 1.17 × 10(-16)]) were strongly associated with the difference between pre- and post-methionine load test tHcy levels (ΔPOST). Of these, one variant in the ALDH1L1 locus, rs2364368, was associated with incident ischemic stroke. Promoter analyses reveal genetic and epigenetic differences that may explain a direct effect on GNMT transcription and a downstream affect on methionine metabolism. Additionally, a genetic-score consisting of the five significant loci explains 13% of the variance of ΔPOST in FHS and 6% of the variance in VISP. Association between variants in FOCM genes with ΔPOST suggest novel mechanisms that lead to differences in methionine metabolism, and possibly the epigenome, impacting disease risk. These data emphasize the importance of a concerted effort to understand regulators of one carbon metabolism as potential therapeutic targets.

Apolipoprotein A-IV (apoA-IV) is a major component of HDL and chylomicron particles and is involved in reverse cholesterol transport. It is an early marker of impaired renal function. We aimed to identify genetic loci associated with apoA-IV concentrations and to investigate relationships with known susceptibility loci for kidney function and lipids. A genome-wide association meta-analysis on apoA-IV concentrations was conducted in five population-based cohorts (n = 13,813) followed by two additional replication studies (n = 2,267) including approximately 10 M SNPs. Three independent SNPs from two genomic regions were significantly associated with apoA-IV concentrations: rs1729407 near APOA4 (P = 6.77 × 10 -  44), rs5104 in APOA4 (P = 1.79 × 10-24) and rs4241819 in KLKB1 (P = 5.6 × 10-14). Additionally, a look-up of the replicated SNPs in downloadable GWAS meta-analysis results was performed on kidney function (defined by eGFR), HDL-cholesterol and triglycerides. From these three SNPs mentioned above, only rs1729407 showed an association with HDL-cholesterol (P = 7.1 × 10 -  07). Moreover, weighted SNP-scores were built involving known susceptibility loci for the aforementioned traits (53, 70 and 38 SNPs, respectively) and were associated with apoA-IV concentrations. This analysis revealed a significant and an inverse association for kidney function with apoA-IV concentrations (P = 5.5 × 10-05). Furthermore, an increase of triglyceride-increasing alleles was found to decrease apoA-IV concentrations (P = 0.0078). In summary, we identified two independent SNPs located in or next the APOA4 gene and one SNP in KLKB1 The association of KLKB1 with apoA-IV suggests an involvement of apoA-IV in renal metabolism and/or an interaction within HDL particles. Analyses of SNP-scores indicate potential causal effects of kidney function and by lesser extent triglycerides on apoA-IV concentrations.

Progranulin is a secreted protein with important functions in processes including immune and inflammatory response, metabolism and embryonic development. The present study aimed at identification of genetic factors determining progranulin concentrations. We conducted a genome-wide association meta-analysis for serum progranulin in three independent cohorts from Europe: Sorbs (N = 848) and KORA (N = 1628) from Germany and PPP-Botnia (N = 335) from Finland (total N = 2811). Single nucleotide polymorphisms (SNPs) associated with progranulin levels were replicated in two additional German cohorts: LIFE-Heart Study (Leipzig; N = 967) and Metabolic Syndrome Berlin Potsdam (Berlin cohort; N = 833). We measured mRNA expression of genes in peripheral blood mononuclear cells (PBMC) by micro-arrays and performed mRNA expression quantitative trait and expression-progranulin association studies to functionally substantiate identified loci. Finally, we conducted siRNA silencing experiments in vitro to validate potential candidate genes within the associated loci. Heritability of circulating progranulin levels was estimated at 31.8% and 26.1% in the Sorbs and LIFE-Heart cohort, respectively. SNPs at three loci reached study-wide significance (rs660240 in CELSR2-PSRC1-MYBPHL-SORT1, rs4747197 in CDH23-PSAP and rs5848 in GRN) explaining 19.4%/15.0% of the variance and 61%/57% of total heritability in the Sorbs/LIFE-Heart Study. The strongest evidence for association was at rs660240 (P = 5.75 × 10-50), which was also associated with mRNA expression of PSRC1 in PBMC (P = 1.51 × 10-21). Psrc1 knockdown in murine preadipocytes led to a consecutive 30% reduction in progranulin secretion. In conclusion, the present meta-GWAS combined with mRNA expression identified three loci associated with progranulin and supports the role of PSRC1 in the regulation of progranulin secretion.

OBJECTIVE: The authors conducted a genome-wide association study of anorexia nervosa and calculated genetic correlations with a series of psychiatric, educational, and metabolic phenotypes. METHOD: Following uniform quality control and imputation procedures using the 1000 Genomes Project (phase 3) in 12 case-control cohorts comprising 3,495 anorexia nervosa cases and 10,982 controls, the authors performed standard association analysis followed by a meta-analysis across cohorts. Linkage disequilibrium score regression was used to calculate genome-wide common variant heritability (single-nucleotide polymorphism [SNP]-based heritability [h2SNP]), partitioned heritability, and genetic correlations (rg) between anorexia nervosa and 159 other phenotypes. RESULTS: Results were obtained for 10,641,224 SNPs and insertion-deletion variants with minor allele frequencies >1% and imputation quality scores >0.6. The h2SNP of anorexia nervosa was 0.20 (SE=0.02), suggesting that a substantial fraction of the twin-based heritability arises from common genetic variation. The authors identified one genome-wide significant locus on chromosome 12 (rs4622308) in a region harboring a previously reported type 1 diabetes and autoimmune disorder locus. Significant positive genetic correlations were observed between anorexia nervosa and schizophrenia, neuroticism, educational attainment, and high-density lipoprotein cholesterol, and significant negative genetic correlations were observed between anorexia nervosa and body mass index, insulin, glucose, and lipid phenotypes. CONCLUSIONS: Anorexia nervosa is a complex heritable phenotype for which this study has uncovered the first genome-wide significant locus. Anorexia nervosa also has large and significant genetic correlations with both psychiatric phenotypes and metabolic traits. The study results encourage a reconceptualization of this frequently lethal disorder as one with both psychiatric and metabolic etiology.

BACKGROUND & AIMS: Hereditary hemochromatosis (HH) is the most common form of genetic iron loading disease. It is mainly related to the homozygous C282Y/C282Y mutation in the HFE gene that is, however, a necessary but not a sufficient condition to develop clinical and even biochemical HH. This suggests that modifier genes are likely involved in the expressivity of the disease. Our aim was to identify such modifier genes. METHODS: We performed a genome-wide association study (GWAS) using DNA collected from 474 unrelated C282Y homozygotes. Associations were examined for both quantitative iron burden indices and clinical outcomes with 534,213 single nucleotide polymorphisms (SNP) genotypes, with replication analyses in an independent sample of 748 C282Y homozygotes from four different European centres. RESULTS: One SNP met genome-wide statistical significance for association with transferrin concentration (rs3811647, GWAS p value of 7×10(-9) and replication p value of 5×10(-13)). This SNP, located within intron 11 of the TF gene, had a pleiotropic effect on serum iron (GWAS p value of 4.9×10(-6) and replication p value of 3.2×10(-6)). Both serum transferrin and iron levels were associated with serum ferritin levels, amount of iron removed and global clinical stage (p<0.01). Serum iron levels were also associated with fibrosis stage (p<0.0001). CONCLUSIONS: This GWAS, the largest one performed so far in unselected HFE-associated HH (HFE-HH) patients, identified the rs3811647 polymorphism in the TF gene as the only SNP significantly associated with iron metabolism through serum transferrin and iron levels. Because these two outcomes were clearly associated with the biochemical and clinical expression of the disease, an indirect link between the rs3811647 polymorphism and the phenotypic presentation of HFE-HH is likely.

BACKGROUND: Thiazide and thiazide-like diuretics are first-line medications for treating uncomplicated hypertension. However, their use has been associated with adverse metabolic events, including hyperglycemia and incident diabetes mellitus, with incompletely understood mechanisms. Our goal was to identify genomic variants associated with thiazide-like diuretic/chlorthalidone-induced glucose change. METHODS AND RESULTS: Genome-wide analysis of glucose change after treatment with chlorthalidone was performed by race among the white (n=175) and black (n=135) participants from the PEAR-2 (Pharmacogenomic Evaluation of Antihypertensive Responses-2). Single-nucleotide polymorphisms with P<5×10-8 were further prioritized using in silico analysis based on their expression quantitative trait loci function. Among blacks, an intronic single-nucleotide polymorphism (rs9943291) in the HMGCS2 was associated with increase in glucose levels following chlorthalidone treatment (ß=12.5; P=4.17×10-8). G-allele carriers of HMGCS2 had higher glucose levels (glucose change=+16.29 mg/dL) post chlorthalidone treatment compared with noncarriers of G allele (glucose change=+2.80 mg/dL). This association was successfully replicated in an independent replication cohort of hydrochlorothiazide-treated participants from the PEAR study (ß=5.54; P=0.023). A meta-analysis of the 2 studies was performed by race in Meta-Analysis Helper, where this single-nucleotide polymorphism, rs9943291, was genome-wide significant with a meta-analysis P value of 3.71×10-8. HMGCS2, a part of the HMG-CoA synthase family, is important for ketogenesis and cholesterol synthesis pathways that are essential in glucose homeostasis. CONCLUSIONS: These results suggest that HMGCS2 is a promising candidate gene involved in chlorthalidone and Hydrochlorothiazide (HCTZ)-induced glucose change. This may provide insights into the mechanisms involved in thiazide-induced hyperglycemia that may ultimately facilitate personalized approaches to antihypertensive selection for hypertension treatment. CLINICAL TRIAL REGISTRATION: URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00246519 and NCT01203852.

Genome-wide association studies (GWAS) in coronary artery disease (CAD) had identified 66 loci at 'genome-wide significance' (P < 5 × 10-8) at the time of this analysis, but a much larger number of putative loci at a false discovery rate (FDR) of 5% (refs. 1,2,3,4). Here we leverage an interim release of UK Biobank (UKBB) data to evaluate the validity of the FDR approach. We tested a CAD phenotype inclusive of angina (SOFT; ncases = 10,801) as well as a stricter definition without angina (HARD; ncases = 6,482) and selected cases with the former phenotype to conduct a meta-analysis using the two most recent CAD GWAS. This approach identified 13 new loci at genome-wide significance, 12 of which were on our previous list of loci meeting the 5% FDR threshold, thus providing strong support that the remaining loci identified by FDR represent genuine signals. The 304 independent variants associated at 5% FDR in this study explain 21.2% of CAD heritability and identify 243 loci that implicate pathways in blood vessel morphogenesis as well as lipid metabolism, nitric oxide signaling and inflammation.

PURPOSE: Acute pancreatitis is one of the common causes of asparaginase intolerance. The mechanism is unknown, and genetic predisposition to asparaginase-induced pancreatitis has not been previously identified. METHODS: To determine clinical risk factors for asparaginase-induced pancreatitis, we studied a cohort of 5,185 children and young adults with acute lymphoblastic leukemia, including 117 (2.3%) who were diagnosed with at least one episode of acute pancreatitis during therapy. A genome-wide association study was performed in the cohort and in an independent case-control group of 213 patients to identify genetic risk factors. RESULTS: Risk factors associated with pancreatitis included genetically defined Native American ancestry (P < .001), older age (P < .001), and higher cumulative dose of asparaginase (P < .001). No common variants reached genome-wide significance in the genome-wide association study, but a rare nonsense variant rs199695765 in CPA2, encoding carboxypeptidase A2, was highly associated with pancreatitis (hazard ratio, 587; 95% CI, 66.8 to 5166; P = 9.0 × 10(-9)). A gene-level analysis showed an excess of additional CPA2 variants in patients who did versus those who did not develop pancreatitis (P = .001). Sixteen CPA2 single-nucleotide polymorphisms were associated (P < .05) with pancreatitis, and 13 of 24 patients who carried at least one of these variants developed pancreatitis. Biologic functions that were overrepresented by common variants modestly associated with pancreatitis included purine metabolism and cytoskeleton regulation. CONCLUSION: Older age, higher exposure to asparaginase, and higher Native American ancestry were independent risk factors for pancreatitis in patients with acute lymphoblastic leukemia. Those who inherit a nonsense rare variant in the CPA2 gene had a markedly increased risk of asparaginase-induced pancreatitis.