Sally John

HLA-B*5701 genotype is a major determinant of drug-induced liver injury due to flucloxacillin

Drug-induced liver injury (DILI) is an important cause of serious liver disease. The antimicrobial agent flucloxacillin is a common cause of DILI, but the genetic basis for susceptibility remains unclear. We conducted a genome-wide association (GWA) study using 866,399 markers in 51 cases of flucloxacillin DILI and 282 controls matched for sex and ancestry. The GWA showed an association peak in the major histocompatibility complex (MHC) region with the strongest association (P = 8.7 × 10−33) seen for rs2395029[G], a marker in complete linkage disequilibrium (LD) with HLA-B*5701. Further MHC genotyping, which included 64 flucloxacillin-tolerant controls, confirmed the association with HLA-B*5701 (OR = 80.6, P = 9.0 × 10−19). The association was replicated in a second cohort of 23 cases. In HLA-B*5701 carrier cases, rs10937275 in ST6GAL1 on chromosome 3 also showed genome-wide significance (OR = 4.1, P = 1.4 × 10−8). These findings provide new insights into the mechanism of flucloxacillin DILI and have the potential to substantially improve diagnosis of this serious disease.

An atlas of genetic influences on human blood metabolites.

Genome-wide association scans with high-throughput metabolic profiling provide unprecedented insights into how genetic variation influences metabolism and complex disease. Here we report the most comprehensive exploration of genetic loci influencing human metabolism thus far, comprising 7,824 adult individuals from 2 European population studies. We report genome-wide significant associations at 145 metabolic loci and their biochemical connectivity with more than 400 metabolites in human blood. We extensively characterize the resulting in vivo blueprint of metabolism in human blood by integrating it with information on gene expression, heritability and overlap with known loci for complex disorders, inborn errors of metabolism and pharmacological targets. We further developed a database and web-based resources for data mining and results visualization. Our findings provide new insights into the role of inherited variation in blood metabolic diversity and identify potential new opportunities for drug development and for understanding disease.

Whole-Genome Scan, in a Complex Disease, Using 11,245 Single-Nucleotide Polymorphisms: Comparison with Microsatellites

Despite the theoretical evidence of the utility of single-nucleotide polymorphisms (SNPs) for linkage analysis, no whole-genome scans of a complex disease have yet been published to directly compare SNPs with microsatellites. Here, we describe a whole-genome screen of 157 families with multiple cases of rheumatoid arthritis (RA), performed using 11,245 genomewide SNPs. The results were compared with those from a 10-cM microsatellite scan in the same cohort. The SNP analysis detected HLA*DRB1, the major RA susceptibility locus (P=.00004), with a linkage interval of 31 cM, compared with a 50-cM linkage interval detected by the microsatellite scan. In addition, four loci were detected at a nominal significance level (P<.05) in the SNP linkage analysis; these were not observed in the microsatellite scan. We demonstrate that variation in information content was the main factor contributing to observed differences in the two scans, with the SNPs providing significantly higher information content than the microsatellites. Reducing the number of SNPs in the marker set to 3,300 (1-cM spacing) caused several loci to drop below nominal significance levels, suggesting that decreases in information content can have significant effects on linkage results. In contrast, differences in maps employed in the analysis, the low detectable rate of genotyping error, and the presence of moderate linkage disequilibrium between markers did not significantly affect the results. We have demonstrated the utility of a dense SNP map for performing linkage analysis in a late-age-at-onset disease, where DNA from parents is not always available. The high SNP density allows loci to be defined more precisely and provides a partial scaffold for association studies, substantially reducing the resource requirement for gene-mapping studies.

Genetically determined P2X7 receptor pore formation regulates variability in chronic pain sensitivity

Chronic pain is highly variable between individuals, as is the response to analgesics. Although much of the variability in chronic pain and analgesic response is heritable, an understanding of the genetic determinants underlying this variability is rudimentary. Here we show that variation within the coding sequence of the gene encoding the P2X7 receptor (P2X7R) affects chronic pain sensitivity in both mice and humans. P2X7Rs, which are members of the family of ionotropic ATP-gated receptors, have two distinct modes of function: they can function through their intrinsic cationic channel or by forming nonselective pores that are permeable to molecules with a mass of up to 900 Da. Using genome-wide linkage analyses, we discovered an association between nerve-injury–induced pain behavior (mechanical allodynia) and the P451L mutation of the mouse P2rx7 gene, such that mice in which P2X7Rs have impaired pore formation as a result of this mutation showed less allodynia than mice with the pore-forming P2rx7 allele. Administration of a peptide corresponding to the P2X7R C-terminal domain, which blocked pore formation but not cation channel activity, selectively reduced nerve injury and inflammatory allodynia only in mice with the pore-forming P2rx7 allele. Moreover, in two independent human chronic pain cohorts, a cohort with pain after mastectomy and a cohort with osteoarthritis, we observed a genetic association between lower pain intensity and the hypofunctional His270 (rs7958311) allele of P2RX7. Our findings suggest that selectively targeting P2X7R pore formation may be a new strategy for individualizing the treatment of chronic pain.

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 (r2 = 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.

Association between rheumatoid arthritis and polymorphism of tumor necrosis factor receptor II, but not tumor necrosis factor receptor I, in Caucasians

OBJECTIVE Tumor necrosis factor (TNF) is a powerful mediator of inflammation in rheumatoid arthritis (RA). In vivo, its acute effects are limited by binding to soluble receptors (TNFR), suggesting that TNFR genes could be important candidate risk factors. The present study was undertaken to investigate association of polymorphisms of TNFRI and TNFRII with RA in subjects in the UK. METHODS Unrelated Caucasian RA patients (n = 291) and healthy Caucasian controls (n = 143) were genotyped for A/G polymorphism in exon 1 of TNFRI. From this sample, 240 of the patients and 137 controls were also typed for a single-nucleotide polymorphism (SNP) in exon 6 of the TNFRII gene. In followup studies, DNA samples from UK Caucasian RA patients with a positive family history (n = 149) and UK Caucasian patients with sporadic RA (n = 208) were also typed for the exon 6 TNFRII polymorphism. RESULTS TNFRI polymorphism was not associated with RA (odds ratio [OR] for GG genotype 0.93, 95% confidence interval [95% CI] 0.54-1.60). For TNFRII, in the initial study group, patients with RA were significantly more likely to be positive for both the G allele and GG genotype than were controls (OR for GG genotype 2.55, 95% CI 1.11-5.86). The association appeared to be confined to those with a family history of RA. This finding was replicated in an independent cohort of patients with familial RA. CONCLUSION The results of this study provide evidence of association between an SNP in the TNFRII gene and RA, the strongest association being observed in patients with a family history. No evidence of association between RA and TNFRI was demonstrated.

Increased Estrogen Rather Than Decreased Androgen Action Is Associated with Longer Androgen Receptor CAG Repeats.

CONTEXT The individual variability in the waning androgenic-anabolic functions of aging men may be influenced by the CAG repeat polymorphism in exon 1 of the androgen receptor (AR), affecting androgen sensitivity. However, findings on its phenotypic effects are inconclusive. OBJECTIVE The aim was to investigate the relationships between health status, various reproductive hormones, and the AR CAG repeat length. DESIGN We conducted a multinational prospective cohort observational study with cross-sectional baseline data. SETTING This was a population survey of community-dwelling men. PARTICIPANTS Men (40-79 yr old; n = 3369) were randomly recruited from centers in eight European countries; CAG repeat analysis was performed in 2878 men. MAIN OUTCOME MEASURES We measured the correlations of the CAG repeat length with selected endocrine, metabolic, and phenotypic parameters related to aging and sex hormone action. RESULTS Only minor differences were found in CAG repeat lengths between the eight European countries. They showed significant positive association with total, free, and bioavailable levels of testosterone (T) and estradiol. FSH but not LH correlated inversely with CAG repeat length. Significant associations were found with bone ultrasound parameters at the calcaneus. Negative correlation was found with triglycerides, but not with other blood lipids or with anthropometry, blood pressure, hemoglobin, insulin sensitivity, or sexual and prostatic functions. CONCLUSIONS The AR CAG repeat length correlates significantly with serum T and estradiol of aging men. Weaker transcriptional activity of the AR with longer CAG-encoded polyglutamine repeats appears to be totally or nearly totally compensated for by higher T levels. The residual phenotypic correlations may reflect differences in estrogen levels/actions after aromatization of the higher T levels.

Metabolomic markers reveal novel pathways of ageing and early development in human populations

Background Human ageing is a complex, multifactorial process and early developmental factors affect health outcomes in old age. Methods Metabolomic profiling on fasting blood was carried out in 6055 individuals from the UK. Stepwise regression was performed to identify a panel of independent metabolites which could be used as a surrogate for age. We also investigated the association with birthweight overall and within identical discordant twins and with genome-wide methylation levels. Results We identified a panel of 22 metabolites which combined are strongly correlated with age (R2 = 59%) and with age-related clinical traits independently of age. One particular metabolite, C-glycosyl tryptophan (C-glyTrp), correlated strongly with age (beta = 0.03, SE = 0.001, P = 7.0 × 10−157) and lung function (FEV1 beta = −0.04, SE = 0.008, P = 1.8 × 10−8 adjusted for age and confounders) and was replicated in an independent population (n = 887). C-glyTrp was also associated with bone mineral density (beta = −0.01, SE = 0.002, P = 1.9 × 10−6) and birthweight (beta = −0.06, SE = 0.01, P = 2.5 × 10−9). The difference in C-glyTrp levels explained 9.4% of the variance in the difference in birthweight between monozygotic twins. An epigenome-wide association study in 172 individuals identified three CpG-sites, associated with levels of C-glyTrp (P < 2 × 10−6). We replicated one CpG site in the promoter of the WDR85 gene in an independent sample of 350 individuals (beta = −0.20, SE = 0.04, P = 2.9 × 10−8). WDR85 is a regulator of translation elongation factor 2, essential for protein synthesis in eukaryotes. Conclusions Our data illustrate how metabolomic profiling linked with epigenetic studies can identify some key molecular mechanisms potentially determined in early development that produce long-term physiological changes influencing human health and ageing.

A Combined Functional Annotation Score for Non-Synonymous Variants

Aims: Next-generation sequencing has opened the possibility of large-scale sequence-based disease association studies. A major challenge in interpreting whole-exome data is predicting which of the discovered variants are deleterious or neutral. To address this question in silico, we have developed a score called Combined Annotation scoRing toOL (CAROL), which combines information from 2 bioinformatics tools: PolyPhen-2 and SIFT, in order to improve the prediction of the effect of non-synonymous coding variants. Methods: We used a weighted Z method that combines the probabilistic scores of PolyPhen-2 and SIFT. We defined 2 dataset pairs to train and test CAROL using information from the dbSNP: ‘HGMD-PUBLIC’ and 1000 Genomes Project databases. The training pair comprises a total of 980 positive control (disease-causing) and 4,845 negative control (non-disease-causing) variants. The test pair consists of 1,959 positive and 9,691 negative controls. Results: CAROL has higher predictive power and accuracy for the effect of non-synonymous variants than each individual annotation tool (PolyPhen-2 and SIFT) and benefits from higher coverage. Conclusion: The combination of annotation tools can help improve automated prediction of whole-genome/exome non-synonymous variant functional consequences.

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.