Genome-wide pharmacogenetics of anti-drug antibody response to bococizumab highlights key residues in HLA DRB1 and DQB1.

In this largest to-date genetic analysis of anti-drug antibody (ADA) response to a therapeutic monoclonal antibody (MAb), genome-wide association was performed for five measures of ADA status among 8844 individuals randomized to bococizumab, which targets PCSK9 for LDL-C lowering and cardiovascular protection. Index associations prioritized specific amino acid substitutions at the DRB1 and DQB1 MHC class II genes rather than canonical haplotypes. Two clusters of missense variants at DRB1 were associated with general ADA measures (residues 9, 11, 13; and 96, 112, 120, 180) and a third cluster of missense variants in DQB1 was associated with ADA measures including neutralizing antibody (NAb) titers (residues 66, 67, 71, 74, 75). The structural disposition of the missense substitutions implicates peptide antigen binding and CD4 effector function, mechanisms that are potentially generalizable to other therapeutic mAbs.Clinicaltrials.gov: NCT01968954, NCT01968967, NCT01968980, NCT01975376, NCT01975389, NCT02100514.

Advancing human genetics research and drug discovery through exome sequencing of the UK Biobank.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

Acute respiratory failure and the kinetics of neutrophil recovery in pediatric hematopoietic cell transplantation: a multicenter study.

In this multicenter study, we investigated the kinetics of neutrophil recovery in relation to acuity and survival among 125 children undergoing allogeneic hematopoietic cell transplantation (allo-HCT) who required invasive mechanical ventilation (IMV). Recovery of neutrophils, whether prior to or after initiation of IMV, was associated with a significantly decreased risk of death relative to never achieving neutrophil recovery. A transient increase in acuity (by oxygenation index and vasopressor requirements) occurred among a subset of the patients who achieved neutrophil recovery after initiation of IMV; 61.5% of these patients survived to discharge from the intensive care unit (ICU). Improved survival among patients who subsequently achieved neutrophil recovery on IMV was not limited to those with peri-engraftment respiratory distress syndrome. The presence of a respiratory pathogen did not affect the risk of death while on IMV but was associated with an increased length of IMV (p < 0.01). Among patients undergoing HCT who develop respiratory failure and require advanced therapeutic support, neutrophil recovery at time of IMV and/or presence of a respiratory pathogen should not be used as determining factors when counseling families about survival.

Correction to: Acute respiratory failure and the kinetics of neutrophil recovery in pediatric hematopoietic cell transplantation: a multicenter study.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Non-alcoholic fatty liver disease and risk of incident acute myocardial infarction and stroke: findings from matched cohort study of 18 million European adults.

OBJECTIVE: To estimate the risk of acute myocardial infarction (AMI) or stroke in adults with non-alcoholic fatty liver disease (NAFLD) or non-alcoholic steatohepatitis (NASH). DESIGN: Matched cohort study. SETTING: Population based, electronic primary healthcare databases before 31 December 2015 from four European countries: Italy (n=1 542 672), Netherlands (n=2 225 925), Spain (n=5 488 397), and UK (n=12 695 046). PARTICIPANTS: 120 795 adults with a recorded diagnosis of NAFLD or NASH and no other liver diseases, matched at time of NAFLD diagnosis (index date) by age, sex, practice site, and visit, recorded at six months before or after the date of diagnosis, with up to 100 patients without NAFLD or NASH in the same database. MAIN OUTCOME MEASURES: Primary outcome was incident fatal or non-fatal AMI and ischaemic or unspecified stroke. Hazard ratios were estimated using Cox models and pooled across databases by random effect meta-analyses. RESULTS: 120 795 patients with recorded NAFLD or NASH diagnoses were identified with mean follow-up 2.1-5.5 years. After adjustment for age and smoking the pooled hazard ratio for AMI was 1.17 (95% confidence interval 1.05 to 1.30; 1035 events in participants with NAFLD or NASH, 67 823 in matched controls). In a group with more complete data on risk factors (86 098 NAFLD and 4 664 988 matched controls), the hazard ratio for AMI after adjustment for systolic blood pressure, type 2 diabetes, total cholesterol level, statin use, and hypertension was 1.01 (0.91 to 1.12; 747 events in participants with NAFLD or NASH, 37 462 in matched controls). After adjustment for age and smoking status the pooled hazard ratio for stroke was 1.18 (1.11 to 1.24; 2187 events in participants with NAFLD or NASH, 134 001 in matched controls). In the group with more complete data on risk factors, the hazard ratio for stroke was 1.04 (0.99 to 1.09; 1666 events in participants with NAFLD, 83 882 in matched controls) after further adjustment for type 2 diabetes, systolic blood pressure, total cholesterol level, statin use, and hypertension. CONCLUSIONS: The diagnosis of NAFLD in current routine care of 17.7 million patient appears not to be associated with AMI or stroke risk after adjustment for established cardiovascular risk factors. Cardiovascular risk assessment in adults with a diagnosis of NAFLD is important but should be done in the same way as for the general population.

African-American, Puerto-Rican, and Other Hispanic Fathers' Differential Responses to a Parenting Intervention.

Despite the recognized importance of fathers to children's well-being, there is a lack of research exploring the impact of parenting interventions on young fathers. Further, little work has been done to identify whether fatherhood interventions differentially benefit specific subgroups of fathers, including Hispanic subgroups.This research examines a 15-week fatherhood intervention for African American, Puerto Rican, and non-Puerto Rican Hispanic young fathers. Data were collected from 312 fathers ages 15-24 at baseline, post-intervention (15 weeks), and at 8 months, 12 months, and 16 months post-baseline follow-up. Latent growth models were used to examine differential shapes of change for the three groups. Results suggest different intervention responses across the three groups on risky parenting attitudes, African American fathers in the study indicated more risky parenting attitudes than Hispanic and Puerto Rican fathers. The practice and research implications of disaggregating data for minority fathers, particularly for Latino subgroups, are discussed.

Risks and clinical predictors of cirrhosis and hepatocellular carcinoma diagnoses in adults with diagnosed NAFLD: real-world study of 18 million patients in four European cohorts.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is a common condition that progresses in some patients to steatohepatitis (NASH), cirrhosis and hepatocellular carcinoma (HCC). Here we used healthcare records of 18 million adults to estimate risk of acquiring advanced liver disease diagnoses in patients with NAFLD or NASH compared to individually matched controls. METHODS: Data were extracted from four European primary care databases representing the UK, Netherlands, Italy and Spain. Patients with a recorded diagnosis of NAFLD or NASH (NAFLD/NASH) were followed up for incident cirrhosis and HCC diagnoses. Each coded NAFLD/NASH patient was matched to up to 100 "non-NAFLD" patients by practice site, gender, age ± 5 years and visit recorded within ± 6 months. Hazard ratios (HR) were estimated using Cox models adjusted for age and smoking status and pooled across databases by random effects meta-analyses. RESULTS: Out of 18,782,281 adults, we identified 136,703 patients with coded NAFLD/NASH. Coded NAFLD/NASH patients were more likely to have diabetes, hypertension and obesity than matched controls. HR for cirrhosis in patients compared to controls was 4.73 (95% CI 2.43-9.19) and for HCC, 3.51 (95% CI 1.72-7.16). HR for either outcome was higher in patients with NASH and those with high-risk Fib-4 scores. The strongest independent predictor of a diagnosis of HCC or cirrhosis was baseline diagnosis of diabetes. CONCLUSIONS: Real-world population data show that recorded diagnosis of NAFLD/NASH increases risk of life-threatening liver outcomes. Diabetes is an independent predictor of advanced liver disease diagnosis, emphasising the need to identify specific groups of patients at highest risk.

An exploratory phenome wide association study linking asthma and liver disease genetic variants to electronic health records from the Estonian Biobank.

The Estonian Biobank, governed by the Institute of Genomics at the University of Tartu (Biobank), has stored genetic material/DNA and continuously collected data since 2002 on a total of 52,274 individuals representing ~5% of the Estonian adult population and is increasing. To explore the utility of data available in the Biobank, we conducted a phenome-wide association study (PheWAS) in two areas of interest to healthcare researchers; asthma and liver disease. We used 11 asthma and 13 liver disease-associated single nucleotide polymorphisms (SNPs), identified from published genome-wide association studies, to test our ability to detect established associations. We confirmed 2 asthma and 5 liver disease associated variants at nominal significance and directionally consistent with published results. We found 2 associations that were opposite to what was published before (rs4374383:AA increases risk of NASH/NAFLD, rs11597086 increases ALT level). Three SNP-diagnosis pairs passed the phenome-wide significance threshold: rs9273349 and E06 (thyroiditis, p = 5.50x10-8); rs9273349 and E10 (type-1 diabetes, p = 2.60x10-7); and rs2281135 and K76 (non-alcoholic liver diseases, including NAFLD, p = 4.10x10-7). We have validated our approach and confirmed the quality of the data for these conditions. Importantly, we demonstrate that the extensive amount of genetic and medical information from the Estonian Biobank can be successfully utilized for scientific research.

Real-world data reveal a diagnostic gap in non-alcoholic fatty liver disease.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is the most common cause of liver disease worldwide. It affects an estimated 20% of the general population, based on cohort studies of varying size and heterogeneous selection. However, the prevalence and incidence of recorded NAFLD diagnoses in unselected real-world health-care records is unknown. We harmonised health records from four major European territories and assessed age- and sex-specific point prevalence and incidence of NAFLD over the past decade. METHODS: Data were extracted from The Health Improvement Network (UK), Health Search Database (Italy), Information System for Research in Primary Care (Spain) and Integrated Primary Care Information (Netherlands). Each database uses a different coding system. Prevalence and incidence estimates were pooled across databases by random-effects meta-analysis after a log-transformation. RESULTS: Data were available for 17,669,973 adults, of which 176,114 had a recorded diagnosis of NAFLD. Pooled prevalence trebled from 0.60% in 2007 (95% confidence interval: 0.41-0.79) to 1.85% (0.91-2.79) in 2014. Incidence doubled from 1.32 (0.83-1.82) to 2.35 (1.29-3.40) per 1000 person-years. The FIB-4 non-invasive estimate of liver fibrosis could be calculated in 40.6% of patients, of whom 29.6-35.7% had indeterminate or high-risk scores. CONCLUSIONS: In the largest primary-care record study of its kind to date, rates of recorded NAFLD are much lower than expected suggesting under-diagnosis and under-recording. Despite this, we have identified rising incidence and prevalence of the diagnosis. Improved recognition of NAFLD may identify people who will benefit from risk factor modification or emerging therapies to prevent progression to cardiometabolic and hepatic complications.

Utilization of Positive and Negative Controls to Examine Comorbid Associations in Observational Database Studies.

BACKGROUND: Opportunities to leverage observational data for precision medicine research are hampered by underlying sources of bias and paucity of methods to handle resulting uncertainty. We outline an approach to account for bias in identifying comorbid associations between 2 rare genetic disorders and type 2 diabetes (T2D) by applying a positive and negative control disease paradigm. RESEARCH DESIGN: Association between 10 common and 2 rare genetic disorders [Hereditary Fructose Intolerance (HFI) and α-1 antitrypsin deficiency] and T2D was compared with the association between T2D and 7 negative control diseases with no established relationship with T2D in 4 observational databases. Negative controls were used to estimate how much bias and variance existed in datasets when no effect should be observed. RESULTS: Unadjusted association for common and rare genetic disorders and T2D was positive and variable in magnitude and distribution in all 4 databases. However, association between negative controls and T2D was 200% greater than expected indicating the magnitude and confidence intervals for comorbid associations are sensitive to systematic bias. A meta-analysis using this method demonstrated a significant association between HFI and T2D but not for α-1 antitrypsin deficiency. CONCLUSIONS: For observational studies, when covariate data are limited or ambiguous, positive and negative controls provide a method to account for the broadest level of systematic bias, heterogeneity, and uncertainty. This provides greater confidence in assessing associations between diseases and comorbidities. Using this approach we were able to demonstrate an association between HFI and T2D. Leveraging real-world databases is a promising approach to identify and corroborate potential targets for precision medicine therapies.

Identifying Cases of Type 2 Diabetes in Heterogeneous Data Sources: Strategy from the EMIF Project.

Due to the heterogeneity of existing European sources of observational healthcare data, data source-tailored choices are needed to execute multi-data source, multi-national epidemiological studies. This makes transparent documentation paramount. In this proof-of-concept study, a novel standard data derivation procedure was tested in a set of heterogeneous data sources. Identification of subjects with type 2 diabetes (T2DM) was the test case. We included three primary care data sources (PCDs), three record linkage of administrative and/or registry data sources (RLDs), one hospital and one biobank. Overall, data from 12 million subjects from six European countries were extracted. Based on a shared event definition, sixteeen standard algorithms (components) useful to identify T2DM cases were generated through a top-down/bottom-up iterative approach. Each component was based on one single data domain among diagnoses, drugs, diagnostic test utilization and laboratory results. Diagnoses-based components were subclassified considering the healthcare setting (primary, secondary, inpatient care). The Unified Medical Language System was used for semantic harmonization within data domains. Individual components were extracted and proportion of population identified was compared across data sources. Drug-based components performed similarly in RLDs and PCDs, unlike diagnoses-based components. Using components as building blocks, logical combinations with AND, OR, AND NOT were tested and local experts recommended their preferred data source-tailored combination. The population identified per data sources by resulting algorithms varied from 3.5% to 15.7%, however, age-specific results were fairly comparable. The impact of individual components was assessed: diagnoses-based components identified the majority of cases in PCDs (93-100%), while drug-based components were the main contributors in RLDs (81-100%). The proposed data derivation procedure allowed the generation of data source-tailored case-finding algorithms in a standardized fashion, facilitated transparent documentation of the process and benchmarking of data sources, and provided bases for interpretation of possible inter-data source inconsistency of findings in future studies.

The epigenetic signature of systemic insulin resistance in obese women.

AIMS/HYPOTHESIS: Insulin resistance (IR) links obesity to type 2 diabetes. The aim of this study was to explore whether white adipose tissue (WAT) epigenetic dysregulation is associated with systemic IR by genome-wide CG dinucleotide (CpG) methylation and gene expression profiling in WAT from insulin-resistant and insulin-sensitive women. A secondary aim was to determine whether the DNA methylation signature in peripheral blood mononuclear cells (PBMCs) reflects WAT methylation and, if so, can be used as a marker for systemic IR. METHODS: From 220 obese women, we selected a total of 80 individuals from either of the extreme ends of the distribution curve of HOMA-IR, an indirect measure of systemic insulin sensitivity. Genome-wide transcriptome and DNA CpG methylation profiling by array was performed on subcutaneous (SAT) and visceral (omental) adipose tissue (VAT). CpG methylation in PBMCs was assayed in the same cohort. RESULTS: There were 647 differentially expressed genes (false discovery rate [FDR] 10%) in SAT, all of which displayed directionally consistent associations in VAT. This suggests that IR is associated with dysregulated expression of a common set of genes in SAT and VAT. The average degree of DNA methylation did not differ between the insulin-resistant and insulin-sensitive group in any of the analysed tissues/cells. There were 223 IR-associated genes in SAT containing a total of 336 nominally significant differentially methylated sites (DMS). The 223 IR-associated genes were over-represented in pathways related to integrin cell surface interactions and insulin signalling and included COL5A1, GAB1, IRS2, PFKFB3 and PTPRJ. In VAT there were a total of 51 differentially expressed genes (FDR 10%); 18 IR-associated genes contained a total of 29 DMS. CONCLUSIONS/INTERPRETATION: In individuals discordant for insulin sensitivity, the average DNA CpG methylation in SAT and VAT is similar, although specific genes, particularly in SAT, display significantly altered expression and DMS in IR, possibly indicating that epigenetic regulation of these genes influences metabolism.

Body Mass Index and Risk of Nonalcoholic Fatty Liver Disease: Two Electronic Health Record Prospective Studies.

CONTEXT: The relationship between rising body mass index (BMI) and prospective risk of nonalcoholic fatty liver disease (NAFLD)/nonalcoholic steatohepatitis (NASH) is virtually absent. OBJECTIVE: Determine the extent of the association between BMI and risk of future NAFLD diagnosis, stratifying by sex and diabetes. DESIGN: Two prospective studies using Humedica and Health Improvement Network (THIN) with 1.54 and 4.96 years of follow-up, respectively. SETTING: Electronic health record databases. PARTICIPANTS: Patients with a recorded BMI measurement between 15 and 60 kg/m(2), and smoking status, and 1 year of active status before baseline BMI. Patients with a diagnosis or history of chronic diseases were excluded. INTERVENTIONS: None. MAIN OUTCOME MEASURE: Recorded diagnosis of NAFLD/NASH during follow-up (Humedica International Classification of Diseases, Ninth Revision code 571.8, and read codes for NAFLD and NASH in THIN). RESULTS: Hazard ratios (HRs) were calculated across BMI categories using BMI of 20-22.5 kg/m(2) as the reference category, adjusting for age, sex, and smoking status. Risk of recorded NAFLD/NASH increased linearly with BMI and was approximately 5-fold higher in Humedica (HR = 4.78; 95% confidence interval, 4.17-5.47) and 9-fold higher in THIN (HR = 8.93; 7.11-11.23) at a BMI of 30-32.5 kg/m(2) rising to around 10-fold higher in Humedica (HR = 9.80; 8.49-11.32) and 14-fold higher in THIN (HR = 14.32; 11.04-18.57) in the 37.5- to 40-kg/m(2) BMI category. Risk of NAFLD/NASH was approximately 50% higher in men and approximately double in those with diabetes. CONCLUSIONS: These data quantify the consistent and strong relationships between BMI and prospectively recorded diagnoses of NAFLD/NASH and emphasize the importance of weight reduction strategies for prevention and management of NAFLD.

Latest approaches for the treatment of obesity.

INTRODUCTION: Obesity is a body weight disorder characterized by excess adiposity that increases the risk for developing co-morbidities such as type 2 diabetes. A large medical need exists for new anti-obesity treatments capable of promoting 10% or greater weight loss, with minimal side effects. AREAS COVERED: The authors describe the application of monogenic forms of rare obesity and genome-wide association studies in selecting critical pathways for drug discovery. Furthermore, they review in detail several pathways and pharmacological targets in the central nervous system (e.g., the leptin-melanocortin axis, the opioid system, GLP-1/GLP-1 system, and FGF21/FGFR1c/ß-Klotho axis) that play an important role in the regulation of feeding behavior and energy homeostasis. Special focus is given to new strategies that engage well-known targets via novel mechanisms in order to circumvent issues seen with previous drug candidates that failed in the clinic. Finally, the authors discuss the recent developments around fixed-dose combinations, targeted polypharmacology, and non-traditional combinations of drugs and devices. EXPERT OPINION: The future for new weight-loss approaches to treat obesity looks promising. Current therapies have shown modest effects on weight loss in the general obese population but will have greater impact in smaller homogeneous sub-populations of obese subjects using personalized medicine. Drug combinations that target multiple, complementary pathways have the potential to promote double-digit weight loss in a broader, heterogeneous patient population. Furthermore, the development of advanced subcutaneous delivery technologies has opened up opportunities to develop breakthrough peptide and biologic agents for the treatment of obesity.

Development and application of spherically curved charge-coupled device imagers.

Operation of a CCD imager on a curved focal surface offers advantages to flat focal planes, especially for lightweight, relatively simple optical systems. The first advantage is that the modulation transfer function can approach diffraction-limited performance for a spherical focal surface employed in large field-of-view or large-format imagers. The second advantage is that a curved focal surface maintains more uniform illumination as a function of radius from the field center. Examples of applications of curved imagers, described here, include a small compact imager and the large curved array used in the Space Surveillance Telescope. The operational characteristics and mechanical limits of an imager deformed to a 15 mm radius are also described.

Genotype Correlations With Blood Pressure and Efficacy From a Randomized Phase III Trial of Second-Line Axitinib Versus Sorafenib in Metastatic Renal Cell Carcinoma.

BACKGROUND: In the phase III axitinib second-line (AXIS) trial, axitinib significantly prolonged progression-free survival (PFS) versus sorafenib in patients with previously treated metastatic renal cell carcinoma (mRCC). Analyses of associations between germline single-nucleotide polymorphisms (SNPs) and outcomes are reported. PATIENTS AND METHODS: DNA samples from blood were genotyped using TaqMan allelic discrimination. Logistic/Cox regression analyses were used to evaluate association of 15 SNPs in vascular endothelial growth factor (VEGF)-A, VEGF receptor (VEGFR)1, VEGFR2, or hypoxia-inducible factor (HIF)-1α with outcomes for blood pressure (BP; Grade ≥ 3 hypertension, diastolic BP > 90 mm Hg, and increase ≥ 15 mm Hg from baseline) and efficacy (independent review committee-assessed objective response rate and PFS, and overall survival [OS]). Multivariate analyses assessed SNPs and baseline characteristics as potential predictors of PFS and OS. RESULTS: Genotype data were available for 305 (42.7%) of 714 patients; 159 received axitinib and 146 sorafenib. After Bonferroni adjustment, no SNP was associated with BP outcomes. In axitinib-treated patients, VEGF-A rs699947 (A/A vs. C/C) and rs833061 (C/C vs. T/T) were associated with longer OS (27.0 vs. 13.4 months; hazard ratio [HR], 0.39; Padjusted = .015). In sorafenib-treated patients, VEGFR2 rs2071559 (G/G vs. A/A) was associated with longer OS (26.8 vs. 13.8 months; HR, 0.41; Padjusted = .030). In multivariate analyses, no SNP predicted axitinib efficacy; VEGFR2 rs2071559 predicted PFS (P = .0053) and OS (P = .0027) for sorafenib. Sensitivity/specificity of VEGFR2 rs2071559 for OS was < 80%. CONCLUSION: No SNP predicted axitinib outcomes. Although VEGFR2 rs2071559 predicted sorafenib efficacy in patients with mRCC, sensitivity/specificity limitations preclude its use for selecting individual patients for sorafenib treatment.

An integrated epigenomic analysis for type 2 diabetes susceptibility loci in monozygotic twins.

DNA methylation has a great potential for understanding the aetiology of common complex traits such as Type 2 diabetes (T2D). Here we perform genome-wide methylated DNA immunoprecipitation sequencing (MeDIP-seq) in whole-blood-derived DNA from 27 monozygotic twin pairs and follow up results with replication and integrated omics analyses. We identify predominately hypermethylated T2D-related differentially methylated regions (DMRs) and replicate the top signals in 42 unrelated T2D cases and 221 controls. The strongest signal is in the promoter of the MALT1 gene, involved in insulin and glycaemic pathways, and related to taurocholate levels in blood. Integrating the DNA methylome findings with T2D GWAS meta-analysis results reveals a strong enrichment for DMRs in T2D-susceptibility loci. We also detect signals specific to T2D-discordant twins in the GPR61 and PRKCB genes. These replicated T2D associations reflect both likely causal and consequential pathways of the disease. The analysis indicates how an integrated genomics and epigenomics approach, utilizing an MZ twin design, can provide pathogenic insights as well as potential drug targets and biomarkers for T2D and other complex traits.

Small molecule ghrelin receptor inverse agonists and antagonists.

Ghrelin is an endogenous peptide hormone secreted primarily by the stomach and is involved in a number of physiological processes including growth hormone secretion, food intake, as well as energy and glucose homeostasis. The physiological actions of ghrelin are mediated through the growth hormone secretagogue receptor 1a (ghrelin receptor), a peptidic G-protein-coupled receptor. This target has attracted much interest, as agents that block ghrelin's actions on its receptor are anticipated to be pharmaceutical interventions for a number of diseases. This review provides an overview of ghrelin biology with a focus on metabolic diseases and summarizes recent medicinal chemistry programs aimed at delivering small molecule ghrelin receptor antagonists and inverse agonists to the clinic.

Differential methylation of the TRPA1 promoter in pain sensitivity.

Chronic pain is a global public health problem, but the underlying molecular mechanisms are not fully understood. Here we examine genome-wide DNA methylation, first in 50 identical twins discordant for heat pain sensitivity and then in 50 further unrelated individuals. Whole-blood DNA methylation was characterized at 5.2 million loci by MeDIP sequencing and assessed longitudinally to identify differentially methylated regions associated with high or low pain sensitivity (pain DMRs). Nine meta-analysis pain DMRs show robust evidence for association (false discovery rate 5%) with the strongest signal in the pain gene TRPA1 (P=1.2 × 10(-13)). Several pain DMRs show longitudinal stability consistent with susceptibility effects, have similar methylation levels in the brain and altered expression in the skin. Our approach identifies epigenetic changes in both novel and established candidate genes that provide molecular insights into pain and may generalize to other complex traits.