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.

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.

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.

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.

Effect of polymorphisms in selected genes involved in pituitary-testicular function on reproductive hormones and phenotype in aging men.

CONTEXT: Polymorphisms in genes involved in regulation, biosynthesis, metabolism, and actions of testicular sex hormones may influence hormone balance and phenotype of aging men. OBJECTIVE: We investigated the relationships between polymorphisms in genes related to pituitary-testicular endocrine function and health status. DESIGN AND SETTING: Using cross-sectional baseline data, we conducted a multinational prospective cohort observational study consisting of a population survey of community-dwelling men. PARTICIPANTS: A total of 2748 men, aged 40-79 (mean +/- sd, 60.2 + 11.2) yr, were randomly recruited from eight European centers. Forty-three polymorphisms were genotyped in the following genes: androgen receptor (AR), estrogen receptor-alpha and -beta (ESR1 and ESR2), steroid 5alpha-reductase type II (SRD5A2), 17alpha-hydroxylase/17,20-lyase (CYP17A1), aromatase (CYP19A1), sex hormone-binding globulin (SHBG), LH beta-subunit (LHB), and LH receptor (LHCGR). MAIN OUTCOME MEASURES: We measured the associations between gene polymorphisms and endocrine, metabolic, and phenotypic parameters related to aging and sex hormone action. RESULTS: Several polymorphisms in SHBG, ESR2, AR, CYP19A1, and LHB were significantly associated with circulating levels of SHBG, LH, total, free, and bioavailable testosterone and estradiol, the LH x testosterone product, and indices of insulin sensitivity. Apart from several previously reported associations between genes affecting estrogen levels and heel ultrasound parameters, no associations existed between polymorphisms and nonhormonal variables (anthropometry, blood lipids, blood pressure, hemoglobin, prostate symptoms, prostate-specific antigen, sexual dysfunction, cognition). CONCLUSION: In aging men, polymorphisms in genes related to the pituitary-testicular endocrine function significantly influence circulating LH, testosterone, and estradiol levels, but the downstream effects may be too small to influence secondary phenotypic parameters.

Genetic variation in sex hormone genes influences heel ultrasound parameters in middle-aged and elderly men: results from the European Male Aging Study (EMAS).

Genes involved in sex hormone pathways are candidates for influencing bone strength. Polymorphisms in these genes were tested for association with heel quantitative ultrasound (QUS) parameters in middle-aged and elderly European men. Men 40-79 yr of age were recruited from population registers in eight European centers for the European Male Aging Study (EMAS). Polymorphisms were genotyped in AR, ESR1, ESR2, CYP19A1, CYP17A1, SHBG, SRD5A2, LHB, and LHCGR. QUS parameters broadband ultrasound attenuation (BUA) and speed of sound (SOS) were measured in the heel and used to derive BMD. The relationships between QUS parameters and polymorphisms were assessed using linear regression adjusting for age and center. A total of 2693 men, with a mean age of 60.1 +/- 11.1 (SD) yr were included in the analysis. Their mean BUA was 80.0 +/- 18.9 dB/Mhz, SOS was 1550.2 +/- 34.1 m/s, and BMD was 0.542 +/- 0.141 g/cm(2). Significant associations were observed between multiple SNPs in a linkage disequilibrium (LD) block within CYP19A1, peaking at the TCT indel with the deletion allele associating with reduced ultrasound BMD in heterozygotes (beta =-0.016, p = -0.005) and homozygotes (beta = -0.029, p = 0.001). The results for BUA and SOS were similar. Significant associations with QUS parameters were also observed for the CAG repeat in AR and SNPs in CYP17A1, LHCGR, and ESR1. Our data confirm evidence of association between bone QUS parameters and polymorphisms in CYP19A1, as well as modest associations with polymorphisms in CYP17A1, ESR1, LHCGR, and AR in a population sample of European men; this supports a role for genetically determined sex hormone actions in influencing male bone health.

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.

Endotoxin exposure, CD14, and allergic disease: an interaction between genes and the environment.

RATIONALE: High endotoxin exposure may reduce the risk of allergic sensitization. OBJECTIVE: To determine the relationship between a promoter polymorphism in the CD14 gene (CD14/-159 C to T) and endotoxin exposure in relation to the development of allergic sensitization, eczema, and wheeze within the setting of a birth cohort. METHODS: We genotyped 442 children (CD14/-159 C to T; rs2569190). We assessed children for allergic sensitization (IgE > 0.2 kU/L to at least one of seven allergens), eczema (physical examination), and parentally reported wheeze. Endotoxin was measured in house dust. MAIN RESULTS: Genotype frequencies were consistent with other populations (TT, 25%; CT, 47%; CC, 28%). Sensitization (present in 33% of children) was not associated with genotype. For children with TT and CT genotypes, there was no association between endotoxin and sensitization (odds ratio [OR], 0.95; 95% confidence interval [CI], 0.71-1.23; p = 0.7; and OR, 0.90; 95% CI, 0.77-1.04; p = 0.16, respectively) or endotoxin and eczema (OR, 0.99; 95% CI, 0.81-1.20; p = 0.89; and OR, 1.38; 95% CI, 0.83-2.30; p = 0.22, respectively). In children with the genotype CC, increasing endotoxin load was associated with a marked and significant reduction in the risk of sensitization (OR, 0.70; 95% CI, 0.55-0.89; p = 0.004) and eczema (OR, 0.73; 95% CI, 0.56-0.95; p = 0.02). However, we observed an increased risk of nonatopic wheeze with increasing endotoxin exposure in children with the CC genotype (OR, 1.42; 95% CI, 1.01-1.99; p = 0.04) but not other genotypes. No effect was seen for atopic wheeze. CONCLUSIONS: Increasing endotoxin exposure is associated with reduced risk of allergic sensitization and eczema but with increased risk of nonatopic wheeze in children with the CC genotype at -159 of the CD14 gene. The impact of environmental endotoxin may be enhanced in individuals with this genotype.

Polymorphisms in a disintegrin and metalloprotease 33 (ADAM33) predict impaired early-life lung function.

RATIONALE: Asthma commonly originates in early life in association with impaired lung function, which tracks to adulthood. OBJECTIVES: Within the context of a prospective birth cohort study, we investigated the association between single nucleotide polymorphisms (SNPs) in a disintegrin and metalloprotease 33 (ADAM33) gene and early-life lung function. METHODS: Children were genotyped for 17 SNPs in ADAM33. Lung function at age 3 (n = 285) and 5 years (n = 470) was assessed using plethysmographic measurement of specific airway resistance (sRaw). At age 5, we also measured FEV(1). SNPs were analyzed individually using logistic regression, followed by linkage disequilibrium mapping to identify the causal locus. MAIN RESULTS: Carriers of the rare allele of F+1 SNP had reduced lung function at age 3 years (p = 0.003). When the recessive model was considered, four SNPs (F+1, S1, ST+5, V4) showed association with sRaw at age 5 years (p < 0.04). Using linkage disequilibrium mapping, we found evidence of a significant causal location between BC+1 and F1 SNPs, at the 5' end of the gene. Four SNPs were associated with lower FEV(1) (F+1, M+1, T1, and T2; p < or = 0.04). The risk of transient early wheezing more than doubled among children homozygous for the A allele of F+1 (odds ratio, 2.39; 95% confidence intervals, 1.18-4.86; p = 0.02), but there was no association between any SNP and allergic sensitization or physician-diagnosed asthma. CONCLUSIONS: Polymorphisms in ADAM33 predict impaired early-life lung function. The functionally relevant polymorphism is likely to be at the 5' end of the gene.