DNA-pools targeted-sequencing as a robust cost-effective method to detect rare variants: Application to dilated cardiomyopathy genetic diagnosis.

Dilated cardiomyopathy (DCM) is a heart disease characterized by left ventricular dilatation and systolic dysfunction. In 30% of cases, pathogenic variants, essentially private to each patient, are identified in at least one of almost 50 reported genes. Thus, while costly, exons capture-based Next Generation Sequencing (NGS) of a targeted gene panel appears as the best strategy to genetically diagnose DCM. Here, we report a NGS strategy applied to pools of 8 DNAs from DCM patients and validate its robustness for rare variants detection at 4-fold reduced cost. Our pipeline uses Freebayes to detect variants with the expected 1/16 allele frequency. From the whole set of detected rare variants in 96 pools we set the variants quality parameters optimizing true positives calling. When compared to simplex DNA sequencing in a shared subset of 50 DNAs, 96% of SNVs/InsDel were accurately identified in pools. Extended to the 384 DNAs included in the study, we detected 100 variants (ACMG class 4 and 5), mostly in well-known morbid gene causing DCM such as TTN, MYH7, FLNC, and TNNT2. To conclude, we report an original pool-sequencing NGS method accurately detecting rare variants. This innovative approach is cost-effective for genetic diagnostic in rare diseases.

Correction: Exome-wide association study reveals novel susceptibility genes to sporadic dilated cardiomyopathy.

[This corrects the article DOI: 10.1371/journal.pone.0172995.].

Genetic Reduction in Left Ventricular Protein Kinase C-α and Adverse Ventricular Remodeling in Human Subjects.

BACKGROUND: Inhibition of PKC-α (protein kinase C-α) enhances contractility and cardioprotection in animal models, but effects in humans are unknown. Genotypes at rs9912468 strongly associate with PRKCA expression in the left ventricle, enabling genetic approaches to measure effects of reduced PKC-α in human populations. METHODS AND RESULTS: We analyzed the cis expression quantitative trait locus for PRKCA marked by rs9912468 using 313 left ventricular specimens from European Ancestry patients. The forward strand minor allele (G) at rs9912468 is associated with reduced PKC-α transcript abundance (1.7-fold reduction in minor allele homozygotes, P=1×10-41). This association was cardiac specific in expression quantitative trait locus data sets that span 16 human tissues. Cardiac epigenomic data revealed a predicted enhancer in complete (R2=1.0) linkage disequilibrium with rs9912468 within intron 2 of PRKCA. We cloned this region and used reporter constructs to verify cardiac-specific enhancer activity in vitro in cardiac and noncardiac cells and in vivo in zebrafish. The PRKCA enhancer contains 2 common genetic variants and 4 haplotypes; the haplotype correlated with the rs9912468 PKC-α-lowering allele (G) showed lowest activity. In contrast to previous reports in animal models, the PKC-α-lowering allele is associated with adverse left ventricular remodeling (higher mass, larger diastolic dimension), reduced fractional shortening, and higher risk of dilated cardiomyopathy in human populations. CONCLUSIONS: These findings support PKC-α as a regulator of the human heart but suggest that PKC-α inhibition may adversely affect the left ventricle depending on timing and duration. Pharmacological studies in human subjects are required to discern potential benefits and harms of PKC-α inhibitors as an approach to treat heart disease.

Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy.

BACKGROUND: Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits, including cardiovascular diseases. RESULTS: Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 97 patients with dilated cardiomyopathy and 108 non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage, and allele-specific expression. Systematic annotation of genome-wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome-wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome-wide association signals in two independent cohorts. CONCLUSIONS: RNA transcription, splicing, and allele-specific expression are each important determinants of the dilated cardiomyopathy phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics.

Exome-wide association study reveals novel susceptibility genes to sporadic dilated cardiomyopathy.

AIMS: Dilated cardiomyopathy (DCM) is an important cause of heart failure with a strong familial component. We performed an exome-wide array-based association study (EWAS) to assess the contribution of missense variants to sporadic DCM. METHODS AND RESULTS: 116,855 single nucleotide variants (SNVs) were analyzed in 2796 DCM patients and 6877 control subjects from 6 populations of European ancestry. We confirmed two previously identified associations with SNVs in BAG3 and ZBTB17 and discovered six novel DCM-associated loci (Q-value<0.01). The lead-SNVs at novel loci are common and located in TTN, SLC39A8, MLIP, FLNC, ALPK3 and FHOD3. In silico fine mapping identified HSPB7 as the most likely candidate at the ZBTB17 locus. Rare variant analysis (MAF<0.01) demonstrated significant association for TTN variants only (P = 0.0085). All candidate genes but one (SLC39A8) exhibit preferential expression in striated muscle tissues and mutations in TTN, BAG3, FLNC and FHOD3 are known to cause familial cardiomyopathy. We also investigated a panel of 48 known cardiomyopathy genes. Collectively, rare (n = 228, P = 0.0033) or common (n = 36, P = 0.019) variants with elevated in silico severity scores were associated with DCM, indicating that the spectrum of genes contributing to sporadic DCM extends beyond those identified here. CONCLUSION: We identified eight loci independently associated with sporadic DCM. The functions of the best candidate genes at these loci suggest that proteostasis regulation might play a role in DCM pathophysiology.

Role of lipid phosphate phosphatase 3 in human aortic endothelial cell function.

AIMS: Lipid phosphate phosphatase 3; type 2 phosphatidic acid phosphatase ß (LPP3; PPAP2B) is a transmembrane protein dephosphorylating and thereby terminating signalling of lipid substrates including lysophosphatidic acid (LPA) and sphingosine-1-phosphate (S1P). Human LPP3 possesses a cell adhesion motif that allows interaction with integrins. A polymorphism (rs17114036) in PPAP2B is associated with coronary artery disease, which prompted us to investigate the possible role of LPP3 in human endothelial dysfunction, a condition promoting atherosclerosis. METHODS AND RESULTS: To study the role of LPP3 in endothelial cells we used human primary aortic endothelial cells (HAECs) in which LPP3 was silenced or overexpressed using either wild type or mutated cDNA constructs. LPP3 silencing in HAECs enhanced secretion of inflammatory cytokines, leucocyte adhesion, cell survival, and migration and impaired angiogenesis, whereas wild-type LPP3 overexpression reversed these effects and induced apoptosis. We also demonstrated that LPP3 expression was negatively correlated with vascular endothelial growth factor expression. Mutations in either the catalytic or the arginine-glycine-aspartate (RGD) domains impaired endothelial cell function and pharmacological inhibition of S1P or LPA restored it. LPA was not secreted in HAECs under silencing or overexpressing LPP3. However, the intra- and extra-cellular levels of S1P tended to be correlated with LPP3 expression, indicating that S1P is probably degraded by LPP3. CONCLUSIONS: We demonstrated that LPP3 is a negative regulator of inflammatory cytokines, leucocyte adhesion, cell survival, and migration in HAECs, suggesting a protective role of LPP3 against endothelial dysfunction in humans. Both the catalytic and the RGD functional domains were involved and S1P, but not LPA, might be the endogenous substrate of LPP3.

Preservation Analysis of Macrophage Gene Coexpression Between Human and Mouse Identifies PARK2 as a Genetically Controlled Master Regulator of Oxidative Phosphorylation in Humans.

Macrophages are key players involved in numerous pathophysiological pathways and an in-depth characterization of their gene regulatory networks can help in better understanding how their dysfunction may impact on human diseases. We here conducted a cross-species network analysis of macrophage gene expression data between human and mouse to identify conserved networks across both species, and assessed whether such networks could reveal new disease-associated regulatory mechanisms. From a sample of 684 individuals processed for genome-wide macrophage gene expression profiling, we identified 27 groups of coexpressed genes (modules). Six modules were found preserved (P < 10-4) in macrophages from 86 mice of the Hybrid Mouse Diversity Panel. One of these modules was significantly [false discovery rate (FDR) = 8.9 × 10-11] enriched for genes belonging to the oxidative phosphorylation (OXPHOS) pathway. This pathway was also found significantly (FDR < 10-4) enriched in susceptibility genes for Alzheimer, Parkinson, and Huntington diseases. We further conducted an expression quantitative trait loci analysis to identify SNP that could regulate macrophage OXPHOS gene expression in humans. This analysis identified the PARK2 rs192804963 as a trans-acting variant influencing (minimal P-value = 4.3 × 10-8) the expression of most OXPHOS genes in humans. Further experimental work demonstrated that PARK2 knockdown expression was associated with increased OXPHOS gene expression in THP1 human macrophages. This work provided strong new evidence that PARK2 participates to the regulatory networks associated with oxidative phosphorylation and suggested that PARK2 genetic variations could act as a trans regulator of OXPHOS gene macrophage expression in humans.

SASH1, a new potential link between smoking and atherosclerosis.

OBJECTIVE: We have previously reported that SASH1 expression is increased in circulating human monocytes from smokers and was positively correlated with the number of carotid atherosclerotic plaques. The aim of this study was to further validate the link between smoking, SASH1 and atherosclerosis within the vascular wall and to assess the impact of SASH1 expression on endothelial cell functions. METHOD: Human carotids with atherosclerotic plaques were obtained from 58 patients (45 of them with known smoking status: smoker, non-smoker, ex-smokers), and were processed for gene expression analyses and immunostaining. To investigate its function, SASH1 was silenced in human aortic endothelial cells (HAECs) using two different siRNA and subcellular localization of SASH1 was determined by immunostaining and subcellular fractionation. Subsequently the transcriptomic analyses and functional experiments (wound healing, WST-1 proliferation or Matrigel assays) were performed to characterize SASH1 function. RESULTS: SASH1 was expressed in human vascular cells (HAECs, smooth muscle cells) and in monocytes/macrophages. Its tissue expression was significantly higher in the atherosclerotic carotids of smokers compared to non-smokers (p < 0.01). In HAECs, SASH1 was expressed mostly in the cytoplasm and SASH1 knockdown resulted in an increased cell migration, proliferation and angiogenesis. Transcriptomic and pathway analyses showed that SASH1 silencing results in a decreased CYP1A1 expression possibly through the inhibition of TP53 activity. CONCLUSION: We showed that SASH1 expression is increased in atherosclerotic carotids in smokers and its silencing affects endothelial angiogenic functions; therefore we provide a potential link between smoking and atherosclerosis through SASH1 expression.

Runs of Homozygosity: Association with Coronary Artery Disease and Gene Expression in Monocytes and Macrophages.

Runs of homozygosity (ROHs) are recognized signature of recessive inheritance. Contributions of ROHs to the genetic architecture of coronary artery disease and regulation of gene expression in cells relevant to atherosclerosis are not known. Our combined analysis of 24,320 individuals from 11 populations of white European ethnicity showed an association between coronary artery disease and both the count and the size of ROHs. Individuals with coronary artery disease had approximately 0.63 (95% CI: 0.4-0.8) excess of ROHs when compared to coronary-artery-disease-free control subjects (p = 1.49 × 10(-9)). The average total length of ROHs was approximately 1,046.92 (95% CI: 634.4-1,459.5) kb greater in individuals with coronary artery disease than control subjects (p = 6.61 × 10(-7)). None of the identified individual ROHs was associated with coronary artery disease after correction for multiple testing. However, in aggregate burden analysis, ROHs favoring increased risk of coronary artery disease were much more common than those showing the opposite direction of association with coronary artery disease (p = 2.69 × 10(-33)). Individual ROHs showed significant associations with monocyte and macrophage expression of genes in their close proximity-subjects with several individual ROHs showed significant differences in the expression of 44 mRNAs in monocytes and 17 mRNAs in macrophages when compared to subjects without those ROHs. This study provides evidence for an excess of homozygosity in coronary artery disease in outbred populations and suggest the potential biological relevance of ROHs in cells of importance to the pathogenesis of atherosclerosis.

Contribution of Rare and Common Genetic Variants to Plasma Lipid Levels and Carotid Stiffness and Geometry: A Substudy of the Paris Prospective Study 3.

BACKGROUND: We assess the contribution of common and rare putatively functional genetic variants (most of them coding) present on the Illumina exome Beadchip to the variability of plasma lipids and stiffness of the common carotid artery. METHODS AND RESULTS: Measurements were obtained from 2283 men and 1398 women, and after filtering and exclusion of monomorphic variants, 32,827 common (minor allele frequency >0.01) and 68,770 rare variants were analyzed. A large fraction of the heritability of plasma lipids is attributable to variants present on the array, especially for triglycerides (fraction of variance attributable to measured genotypes: V(G)/V(p)=31.4%, P<3.1×10(-11)) and high-density lipoprotein cholesterol (V(G)/V(p)=26.4%, P<4.2×10(-12)). Plasma lipids were associated with common variants located in known candidate genes, but no implication of rare variants could be established. Gene sets for plasma lipids, blood pressure, and coronary artery disease were defined on the basis of recent meta-analyses of genome-wide association studies. We observed a strong association between the plasma lipids gene set and plasma lipid variables, but none of the 3 genome-wide association studies gene sets was associated with the carotid parameters. Significant V(G)/V(p) ratios were observed for external (14.5%, P<2.7×10(-5)) and internal diameter (13.4%, P<4.3×10(-4)), stiffness (12.5%, P<8.0×10(-4)), intima-media thickness (10.6%, P<7.9×10(-4)), and wall cross-sectional area (13.2%, P<2.4×10(-5)). A significant association was observed between the common rs2903692 polymorphism of the CLEC16A gene and the internal diameter (P<4.3×10(-7)). CONCLUSIONS: These results suggest an involvement of CLEC16A, a gene that has been reported to be associated with immune disorders, in the modulation of carotid vasodilatation.

Meta-analysis of 65,734 individuals identifies TSPAN15 and SLC44A2 as two susceptibility loci for venous thromboembolism.

Venous thromboembolism (VTE), the third leading cause of cardiovascular mortality, is a complex thrombotic disorder with environmental and genetic determinants. Although several genetic variants have been found associated with VTE, they explain a minor proportion of VTE risk in cases. We undertook a meta-analysis of genome-wide association studies (GWASs) to identify additional VTE susceptibility genes. Twelve GWASs totaling 7,507 VTE case subjects and 52,632 control subjects formed our discovery stage where 6,751,884 SNPs were tested for association with VTE. Nine loci reached the genome-wide significance level of 5 × 10(-8) including six already known to associate with VTE (ABO, F2, F5, F11, FGG, and PROCR) and three unsuspected loci. SNPs mapping to these latter were selected for replication in three independent case-control studies totaling 3,009 VTE-affected individuals and 2,586 control subjects. This strategy led to the identification and replication of two VTE-associated loci, TSPAN15 and SLC44A2, with lead risk alleles associated with odds ratio for disease of 1.31 (p = 1.67 × 10(-16)) and 1.21 (p = 2.75 × 10(-15)), respectively. The lead SNP at the TSPAN15 locus is the intronic rs78707713 and the lead SLC44A2 SNP is the non-synonymous rs2288904 previously shown to associate with transfusion-related acute lung injury. We further showed that these two variants did not associate with known hemostatic plasma markers. TSPAN15 and SLC44A2 do not belong to conventional pathways for thrombosis and have not been associated to other cardiovascular diseases nor related quantitative biomarkers. Our findings uncovered unexpected actors of VTE etiology and pave the way for novel mechanistic concepts of VTE pathophysiology.

Cigarette smoking reduces DNA methylation levels at multiple genomic loci but the effect is partially reversible upon cessation.

Smoking is a major risk factor in many diseases. Genome wide association studies have linked genes for nicotine dependence and smoking behavior to increased risk of cardiovascular, pulmonary, and malignant diseases. We conducted an epigenome wide association study in peripheral-blood DNA in 464 individuals (22 current smokers and 263 ex-smokers), using the Human Methylation 450 K array. Upon replication in an independent sample of 356 twins (41 current and 104 ex-smokers), we identified 30 probes in 15 distinct loci, all of which reached genome-wide significance in the combined analysis P < 5 × 10(-8). All but one probe (cg17024919) remained significant after adjusting for blood cell counts. We replicated all 9 known loci and found an independent signal at CPOX near GPR15. In addition, we found 6 new loci at PRSS23, AVPR1B, PSEN2, LINC00299, RPS6KA2, and KIAA0087. Most of the lead probes (13 out of 15) associated with cigarette smoking, overlapped regions of open chromatin (FAIRE and DNaseI hypersensitive sites) or/and H3K27Ac peaks (ENCODE data set), which mark regulatory elements. The effect of smoking on DNA methylation was partially reversible upon smoking cessation for longer than 3 months. We report the first statistically significant interaction between a SNP (rs2697768) and cigarette smoking on DNA methylation (cg03329539). We provide evidence that the metSNP for cg03329539 regulates expression of the CHRND gene located circa 95 Kb downstream of the methylation site. Our findings suggest the existence of dynamic, reversible site-specific methylation changes in response to cigarette smoking , which may contribute to the extended health risks associated with cigarette smoking.

Allelic expression mapping across cellular lineages to establish impact of non-coding SNPs.

Most complex disease-associated genetic variants are located in non-coding regions and are therefore thought to be regulatory in nature. Association mapping of differential allelic expression (AE) is a powerful method to identify SNPs with direct cis-regulatory impact (cis-rSNPs). We used AE mapping to identify cis-rSNPs regulating gene expression in 55 and 63 HapMap lymphoblastoid cell lines from a Caucasian and an African population, respectively, 70 fibroblast cell lines, and 188 purified monocyte samples and found 40-60% of these cis-rSNPs to be shared across cell types. We uncover a new class of cis-rSNPs, which disrupt footprint-derived de novo motifs that are predominantly bound by repressive factors and are implicated in disease susceptibility through overlaps with GWAS SNPs. Finally, we provide the proof-of-principle for a new approach for genome-wide functional validation of transcription factor-SNP interactions. By perturbing NFκB action in lymphoblasts, we identified 489 cis-regulated transcripts with altered AE after NFκB perturbation. Altogether, we perform a comprehensive analysis of cis-variation in four cell populations and provide new tools for the identification of functional variants associated to complex diseases.

Single nucleotide polymorphisms with cis-regulatory effects on long non-coding transcripts in human primary monocytes.

We applied genome-wide allele-specific expression analysis of monocytes from 188 samples. Monocytes were purified from white blood cells of healthy blood donors to detect cis-acting genetic variation that regulates the expression of long non-coding RNAs. We analysed 8929 regions harboring genes for potential long non-coding RNA that were retrieved from data from the ENCODE project. Of these regions, 60% were annotated as intergenic, which implies that they do not overlap with protein-coding genes. Focusing on the intergenic regions, and using stringent analysis of the allele-specific expression data, we detected robust cis-regulatory SNPs in 258 out of 489 informative intergenic regions included in the analysis. The cis-regulatory SNPs that were significantly associated with allele-specific expression of long non-coding RNAs were enriched to enhancer regions marked for active or bivalent, poised chromatin by histone modifications. Out of the lncRNA regions regulated by cis-acting regulatory SNPs, 20% (n = 52) were co-regulated with the closest protein coding gene. We compared the identified cis-regulatory SNPs with those in the catalog of SNPs identified by genome-wide association studies of human diseases and traits. This comparison identified 32 SNPs in loci from genome-wide association studies that displayed a strong association signal with allele-specific expression of non-coding RNAs in monocytes, with p-values ranging from 6.7×10(-7) to 9.5×10(-89). The identified cis-regulatory SNPs are associated with diseases of the immune system, like multiple sclerosis and rheumatoid arthritis.

Adrenomedullin and arterial stiffness: integrative approach combining monocyte ADM expression, plasma MR-Pro-ADM, and genome-wide association study.

BACKGROUND: Adrenomedullin (ADM) is a circulating vasoactive peptide involved in vascular homeostasis and endothelial function. Single nucleotide polymorphisms of the ADM gene are associated with blood pressure variability, and elevated levels of plasma midregional proadrenomedullin (MR-pro-ADM) are associated with cardiovascular diseases. METHODS AND RESULTS: We investigated the sources of variability of ADM gene expression and plasma MR-pro-ADM concentrations in the general population, and their relationship with markers of atherosclerosis. MR-pro-ADM levels were assessed in 4155 individuals who underwent evaluation of carotid intima-media thickness and arterial rigidity (reflection index and stiffness index). In a subsample of 1372 individuals, ADM gene expression was assessed as part of a transcriptomic study of circulating monocytes. Nongenetic factors explained 45.8% and 7.5% of MR-pro-ADM and ADM expression variability, respectively. ADM expression correlated with plasma C-reactive protein, interleukin-receptor 1A, and myeloperoxidase, whereas MR-pro-ADM levels correlated with C-terminal proendothelin-1, creatinine, and N-terminal pro-B-type natriuretic peptide. Genome-wide association study of ADM expression and MR-pro-ADM levels both identified a single locus encompassing the ADM gene. ADM expression was associated with 1 single nucleotide polymorphism rs11042717 (P=2.36×10(-12)), whereas MR-pro-ADM was associated with 2 single nucleotide polymorphisms with additive effects, rs2957692 (P=1.54×10(-13)) and rs2957717 (P=4.24×10(-8)). Reflection index was independently associated with rs11042717 (P<10(-4)) and ADM expression (P=0.0002) but not with MR-pro-ADM. Weaker associations were observed for stiffness index. Intima-media thickness was not related to ADM single nucleotide polymorphisms or expression. CONCLUSIONS: These results support an involvement of the ADM gene in the modulation of peripheral vascular tone.

Human CalDAG-GEFI gene (RASGRP2) mutation affects platelet function and causes severe bleeding.

The nature of an inherited platelet disorder was investigated in three siblings affected by severe bleeding. Using whole-exome sequencing, we identified the culprit mutation (cG742T) in the RAS guanyl-releasing protein-2 (RASGRP2) gene coding for calcium- and DAG-regulated guanine exchange factor-1 (CalDAG-GEFI). Platelets from individuals carrying the mutation present a reduced ability to activate Rap1 and to perform proper αIIbß3 integrin inside-out signaling. Expression of CalDAG-GEFI mutant in HEK293T cells abolished Rap1 activation upon stimulation. Nevertheless, the PKC- and ADP-dependent pathways allow residual platelet activation in the absence of functional CalDAG-GEFI. The mutation impairs the platelet's ability to form thrombi under flow and spread normally as a consequence of reduced Rac1 GTP-binding. Functional deficiencies were confined to platelets and megakaryocytes with no leukocyte alteration. This contrasts with the phenotype seen in type III leukocyte adhesion deficiency caused by the absence of kindlin-3. Heterozygous did not suffer from bleeding and have normal platelet aggregation; however, their platelets mimicked homozygous ones by failing to undergo normal adhesion under flow and spreading. Rescue experiments on cultured patient megakaryocytes corrected the functional deficiency after transfection with wild-type RASGRP2. Remarkably, the presence of a single normal allele is sufficient to prevent bleeding, making CalDAG-GEFI a novel and potentially safe therapeutic target to prevent thrombosis.

DNA methylation and body-mass index: a genome-wide analysis.

BACKGROUND: Obesity is a major health problem that is determined by interactions between lifestyle and environmental and genetic factors. Although associations between several genetic variants and body-mass index (BMI) have been identified, little is known about epigenetic changes related to BMI. We undertook a genome-wide analysis of methylation at CpG sites in relation to BMI. METHODS: 479 individuals of European origin recruited by the Cardiogenics Consortium formed our discovery cohort. We typed their whole-blood DNA with the Infinium HumanMethylation450 array. After quality control, methylation levels were tested for association with BMI. Methylation sites showing an association with BMI at a false discovery rate q value of 0·05 or less were taken forward for replication in a cohort of 339 unrelated white patients of northern European origin from the MARTHA cohort. Sites that remained significant in this primary replication cohort were tested in a second replication cohort of 1789 white patients of European origin from the KORA cohort. We examined whether methylation levels at identified sites also showed an association with BMI in DNA from adipose tissue (n=635) and skin (n=395) obtained from white female individuals participating in the MuTHER study. Finally, we examined the association of methylation at BMI-associated sites with genetic variants and with gene expression. FINDINGS: 20 individuals from the discovery cohort were excluded from analyses after quality-control checks, leaving 459 participants. After adjustment for covariates, we identified an association (q value ≤0·05) between methylation at five probes across three different genes and BMI. The associations with three of these probes--cg22891070, cg27146050, and cg16672562, all of which are in intron 1 of HIF3A--were confirmed in both the primary and second replication cohorts. For every 0·1 increase in methylation ß value at cg22891070, BMI was 3·6% (95% CI 2·4-4·9) higher in the discovery cohort, 2·7% (1·2-4·2) higher in the primary replication cohort, and 0·8% (0·2-1·4) higher in the second replication cohort. For the MuTHER cohort, methylation at cg22891070 was associated with BMI in adipose tissue (p=1·72 × 10(-5)) but not in skin (p=0·882). We observed a significant inverse correlation (p=0·005) between methylation at cg22891070 and expression of one HIF3A gene-expression probe in adipose tissue. Two single nucleotide polymorphisms--rs8102595 and rs3826795--had independent associations with methylation at cg22891070 in all cohorts. However, these single nucleotide polymorphisms were not significantly associated with BMI. INTERPRETATION: Increased BMI in adults of European origin is associated with increased methylation at the HIF3A locus in blood cells and in adipose tissue. Our findings suggest that perturbation of hypoxia inducible transcription factor pathways could have an important role in the response to increased weight in people. FUNDING: The European Commission, National Institute for Health Research, British Heart Foundation, and Wellcome Trust.

A meta-analysis of genome-wide association studies identifies ORM1 as a novel gene controlling thrombin generation potential.

Thrombin, the major enzyme of the hemostatic system, is involved in biological processes associated with several human diseases. The capacity of a given individual to generate thrombin, called the thrombin generation potential (TGP), can be robustly measured in plasma and was shown to associate with thrombotic disorders. To investigate the genetic architecture underlying the interindividual TGP variability, we conducted a genome-wide association study in 2 discovery samples (N = 1967) phenotyped for 3 TGP biomarkers, the endogenous thrombin potential, the peak height, and the lag time, and replicated the main findings in 2 independent studies (N = 1254). We identified the ORM1 gene, coding for orosomucoid, as a novel locus associated with lag time variability, reflecting the initiation process of thrombin generation with a combined P value of P = 7.1 × 10(-15) for the lead single nucleotide polymorphism (SNP) (rs150611042). This SNP was also observed to associate with ORM1 expression in monocytes (P = 8.7 × 10(-10)) and macrophages (P = 3.2 × 10(-3)). In vitro functional experiments further demonstrated that supplementing normal plasma with increasing orosomucoid concentrations was associated with impaired thrombin generation. These results pave the way for novel mechanistic pathways and therapeutic perspectives in the etiology of thrombin-related disorders.

Increased monocyte adhesion by endothelial expression of VCAM-1 missense variation in vitro.

OBJECTIVE: In whole genome and single gene analyses, genetic variation at the vascular cell adhesion molecule-1 (VCAM-1) locus has been associated with inflammatory disease and stroke in sickle cell anaemia. In the current work, we investigated the functional impact of VCAM-1 missense variants and their effect on cell-cell adhesion. METHODS AND RESULTS: To determine the functional in vitro relevance of five missense VCAM-1 variants (S318F; T384A; G413A; L555V; I716L), we generated wild type and single variant VCAM-1 forms [318F, 384A, 413A, 555V, 716L] in EA.hy926 endothelial cells. Real-time PCR, western blot and ELISA analyses revealed significant differences in mRNA and protein levels for VCAM-1 variants. Monocytic cell lines THP-1 and U937 showed significantly increased adhesion to endothelial cells overexpressing VCAM-1 forms 318F, 555V and 716L compared to those overexpressing wild type VCAM-1 (p < 0.05). CONCLUSIONS: VCAM-1-dependent cell adhesion to endothelial cells in vitro is significantly increased when expressing VCAM-1 missense mutations 318F, 555V and 716L. The underlying mechanism involves altered VCAM-1 protein levels and function. This observation may be of particular relevance for chronic inflammatory pathophysiologic conditions involving cell-cell adhesion such as atherosclerosis and other proinflammatory conditions.

GUESS-ing polygenic associations with multiple phenotypes using a GPU-based evolutionary stochastic search algorithm.

Genome-wide association studies (GWAS) yielded significant advances in defining the genetic architecture of complex traits and disease. Still, a major hurdle of GWAS is narrowing down multiple genetic associations to a few causal variants for functional studies. This becomes critical in multi-phenotype GWAS where detection and interpretability of complex SNP(s)-trait(s) associations are complicated by complex Linkage Disequilibrium patterns between SNPs and correlation between traits. Here we propose a computationally efficient algorithm (GUESS) to explore complex genetic-association models and maximize genetic variant detection. We integrated our algorithm with a new Bayesian strategy for multi-phenotype analysis to identify the specific contribution of each SNP to different trait combinations and study genetic regulation of lipid metabolism in the Gutenberg Health Study (GHS). Despite the relatively small size of GHS (n  =  3,175), when compared with the largest published meta-GWAS (n > 100,000), GUESS recovered most of the major associations and was better at refining multi-trait associations than alternative methods. Amongst the new findings provided by GUESS, we revealed a strong association of SORT1 with TG-APOB and LIPC with TG-HDL phenotypic groups, which were overlooked in the larger meta-GWAS and not revealed by competing approaches, associations that we replicated in two independent cohorts. Moreover, we demonstrated the increased power of GUESS over alternative multi-phenotype approaches, both Bayesian and non-Bayesian, in a simulation study that mimics real-case scenarios. We showed that our parallel implementation based on Graphics Processing Units outperforms alternative multi-phenotype methods. Beyond multivariate modelling of multi-phenotypes, our Bayesian model employs a flexible hierarchical prior structure for genetic effects that adapts to any correlation structure of the predictors and increases the power to identify associated variants. This provides a powerful tool for the analysis of diverse genomic features, for instance including gene expression and exome sequencing data, where complex dependencies are present in the predictor space.