An arrhythmogenic metabolite in atrial fibrillation.

BACKGROUND: Long-chain acyl-carnitines (ACs) are potential arrhythmogenic metabolites. Their role in atrial fibrillation (AF) remains incompletely understood. Using a systems medicine approach, we assessed the contribution of C18:1AC to AF by analysing its in vitro effects on cardiac electrophysiology and metabolism, and translated our findings into the human setting. METHODS AND RESULTS: Human iPSC-derived engineered heart tissue was exposed to C18:1AC. A biphasic effect on contractile force was observed: short exposure enhanced contractile force, but elicited spontaneous contractions and impaired Ca2+ handling. Continuous exposure provoked an impairment of contractile force. In human atrial mitochondria from AF individuals, C18:1AC inhibited respiration. In a population-based cohort as well as a cohort of patients, high C18:1AC serum concentrations were associated with the incidence and prevalence of AF. CONCLUSION: Our data provide evidence for an arrhythmogenic potential of the metabolite C18:1AC. The metabolite interferes with mitochondrial metabolism, thereby contributing to contractile dysfunction and shows predictive potential as novel circulating biomarker for risk of AF.

Genome-wide association study of a semicontinuous trait: illustration of the impact of the modeling strategy through the study of Neutrophil Extracellular Traps levels.

Over the last years, there has been a considerable expansion of genome-wide association studies (GWAS) for discovering biological pathways underlying pathological conditions or disease biomarkers. These GWAS are often limited to binary or quantitative traits analyzed through linear or logistic models, respectively. In some situations, the distribution of the outcome may require more complex modeling, such as when the outcome exhibits a semicontinuous distribution characterized by an excess of zero values followed by a non-negative and right-skewed distribution. We here investigate three different modeling for semicontinuous data: Tobit, Negative Binomial and Compound Poisson-Gamma. Using both simulated data and a real GWAS on Neutrophil Extracellular Traps (NETs), an emerging biomarker in immuno-thrombosis, we demonstrate that Compound Poisson-Gamma was the most robust model with respect to low allele frequencies and outliers. This model further identified the MIR155HG locus as significantly (P = 1.4 × 10-8) associated with NETs plasma levels in a sample of 657 participants, a locus recently highlighted to be involved in NETs formation in mice. This work highlights the importance of the modeling strategy for GWAS of a semicontinuous outcome and suggests Compound Poisson-Gamma as an elegant but neglected alternative to Negative Binomial for modeling semicontinuous outcome in the context of genomic investigations.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood.

BACKGROUND: Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa. OBJECTIVES: To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions. METHODS: We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation. RESULTS: We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant. CONCLUSION: We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors.

BACKGROUND: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources. METHODS: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations. RESULTS: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis. CONCLUSIONS: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Common and Rare 5'UTR Variants Altering Upstream Open Reading Frames in Cardiovascular Genomics.

High-throughput sequencing (HTS) technologies are revolutionizing the research and molecular diagnosis landscape by allowing the exploration of millions of nucleotide sequences at an unprecedented scale. These technologies are of particular interest in the identification of genetic variations contributing to the risk of rare (Mendelian) and common (multifactorial) human diseases. So far, they have led to numerous successes in identifying rare disease-causing mutations in coding regions, but few in non-coding regions that include introns, untranslated (UTR), and intergenic regions. One class of neglected non-coding variations is that of 5'UTR variants that alter upstream open reading frames (upORFs) of the coding sequence (CDS) of a natural protein coding transcript. Following a brief summary of the molecular bases of the origin and functions of upORFs, we will first review known 5'UTR variations altering upORFs and causing rare cardiovascular disorders (CVDs). We will then investigate whether upORF-affecting single nucleotide polymorphisms could be good candidates for explaining association signals detected in the context of genome-wide association studies for common complex CVDs.

CRIP1 expression in monocytes related to hypertension.

Hypertension is a complex and multifactorial disorder caused by lifestyle and environmental factors, inflammation and disease-related genetic factors and is a risk factor for stroke, ischemic heart disease and renal failure. Although circulating monocytes and tissue macrophages contribute to the pathogenesis of hypertension, the underlying mechanisms are poorly understood. Cysteine rich protein 1 (CRIP1) is highly expressed in immune cells, and CRIP1 mRNA expression in monocytes associates with blood pressure (BP) and is up-regulated by proinflammatory modulation suggesting a link between CRIP1 and BP regulation through the immune system. To address this functional link, we studied CRIP1 expression in immune cells in relation to BP using a human cohort study and hypertensive mouse models. CRIP1 expression in splenic monocytes/macrophages and in circulating monocytes was significantly affected by angiotensin II (Ang II) in a BP-elevating dose (2 mg/kg/day). In the human cohort study, monocytic CRIP1 expression levels were associated with elevated BP, whereas upon differentiation of monocytes to macrophages this association along with the CRIP1 expression level was diminished. In conclusion, CRIP1-positive circulating and splenic monocytes seem to play an important role in hypertension related inflammatory processes through endogenous hormones such as Ang II. These findings suggest that CRIP1 may affect the interaction between the immune system, in particular monocytes, and the pathogenesis of hypertension.

Epigenome-wide association study identifies DNA methylation markers for asthma remission in whole blood and nasal epithelium.

BACKGROUND: Asthma is a chronic respiratory disease which is not curable, yet some patients experience spontaneous remission. We hypothesized that epigenetic mechanisms may be involved in asthma remission. METHODS: Clinical remission (ClinR) was defined as the absence of asthma symptoms and medication for at least 12 months, and complete remission (ComR) was defined as ClinR with normal lung function and absence of airway hyperresponsiveness. We analyzed differential DNA methylation of ClinR and ComR comparing to persistent asthma (PersA) in whole blood samples (n = 72) and nasal brushing samples (n = 97) in a longitudinal cohort of well characterized asthma patients. Significant findings of whole blood DNA methylation were tested for replication in two independent cohorts, Lifelines and Epidemiological study on the Genetics and Environment of Asthma (EGEA). RESULTS: We identified differentially methylated CpG sites associated with ClinR (7 CpG sites) and ComR (129 CpG sites) in whole blood. One CpG (cg13378519, Chr1) associated with ClinR was replicated and annotated to PEX11 (Peroxisomal Biogenesis Factor 11 Beta). The whole blood DNA methylation levels of this CpG were also different between ClinR and healthy subjects. One ComR-associated CpG (cg24788483, Chr10) that annotated to TCF7L2 (Transcription Factor 7 Like 2) was replicated and associated with expression of TCF7L2 gene. One out of seven ClinR-associated CpG sites and 8 out of 129 ComR-associated CpG sites identified from whole blood samples showed nominal significance (P < 0.05) and the same direction of effect in nasal brushes. CONCLUSION: We identified DNA methylation markers possibly associated with clinical and complete asthma remission in nasal brushes and whole blood, and two CpG sites identified from whole blood can be replicated in independent cohorts and may play a role in peroxisome proliferation and Wnt signaling pathway.

Correction to: Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

Following publication of the original article [1], the authors reported an error in Additional file 1.

bioSyntax: syntax highlighting for computational biology.

BACKGROUND: Computational biology requires the reading and comprehension of biological data files. Plain-text formats such as SAM, VCF, GTF, PDB and FASTA, often contain critical information which is obfuscated by the data structure complexity. RESULTS: bioSyntax ( https://biosyntax.org/ ) is a freely available suite of biological syntax highlighting packages for vim, gedit, Sublime, VSCode, and less. bioSyntax improves the legibility of low-level biological data in the bioinformatics workspace. CONCLUSION: bioSyntax supports computational scientists in parsing and comprehending their data efficiently and thus can accelerate research output.

Genetically defined elevated homocysteine levels do not result in widespread changes of DNA methylation in leukocytes.

BACKGROUND: DNA methylation is affected by the activities of the key enzymes and intermediate metabolites of the one-carbon pathway, one of which involves homocysteine. We investigated the effect of the well-known genetic variant associated with mildly elevated homocysteine: MTHFR 677C>T independently and in combination with other homocysteine-associated variants, on genome-wide leukocyte DNA-methylation. METHODS: Methylation levels were assessed using Illumina 450k arrays on 9,894 individuals of European ancestry from 12 cohort studies. Linear-mixed-models were used to study the association of additive MTHFR 677C>T and genetic-risk score (GRS) based on 18 homocysteine-associated SNPs, with genome-wide methylation. RESULTS: Meta-analysis revealed that the MTHFR 677C>T variant was associated with 35 CpG sites in cis, and the GRS showed association with 113 CpG sites near the homocysteine-associated variants. Genome-wide analysis revealed that the MTHFR 677C>T variant was associated with 1 trans-CpG (nearest gene ZNF184), while the GRS model showed association with 5 significant trans-CpGs annotated to nearest genes PTF1A, MRPL55, CTDSP2, CRYM and FKBP5. CONCLUSIONS: Our results do not show widespread changes in DNA-methylation across the genome, and therefore do not support the hypothesis that mildly elevated homocysteine is associated with widespread methylation changes in leukocytes.

Homocysteine levels associate with subtle changes in leukocyte DNA methylation: an epigenome-wide analysis.

AIM: Homocysteine (Hcy) is a sensitive marker of one-carbon metabolism. Higher Hcy levels have been associated with global DNA hypomethylation. We investigated the association between plasma Hcy and epigenome-wide DNA methylation in leukocytes. METHODS: Methylation was measured using Illumina 450 k arrays in 2035 individuals from six cohorts. Hcy-associated differentially methylated positions and regions were identified using meta-analysis. RESULTS: Three differentially methylated positions cg21607669 (SLC27A1), cg26382848 (AJUBA) and cg10701000 (KCNMA1) at chromosome 19, 14 and 10, respectively, were significantly associated with Hcy. In addition, we identified 68 Hcy-associated differentially methylated regions, the most significant of which was a 1.8-kb spanning domain (TNXB/ATF6B) at chromosome 6. CONCLUSION: We identified novel epigenetic loci associated with Hcy levels, of which specific role needs to be further validated.

Blood triglyceride levels are associated with DNA methylation at the serine metabolism gene PHGDH.

Efficient interventions to reduce blood triglycerides are few; newer and more tolerable intervention targets are needed. Understanding the molecular mechanisms underlying blood triglyceride levels variation is key to identifying new therapies. To explore the role of epigenetic mechanisms on triglyceride levels, a blood methylome scan was conducted in 199 individuals from 5 French-Canadian families ascertained on venous thromboembolism, and findings were replicated in 324 French unrelated patients with venous thromboembolism. Genetic context and functional relevance were investigated. Two DNA methylation sites associated with triglyceride levels were identified. The first one, located in the ABCG1 gene, was recently reported, whereas the second one, located in the promoter of the PHGDH gene, is novel. The PHGDH methylation site, cg14476101, was found to be associated with variation in triglyceride levels in a threshold manner: cg14476101 was inversely associated with triglyceride levels only when triglyceride levels were above 1.12 mmol/L (discovery P-value = 8.4 × 10-6; replication P-value = 0.0091). Public databases findings supported a functional role of cg14476101 on PHGDH expression. PHGDH catalyses the first step in the serine biosynthesis pathway. These findings highlight the role of epigenetic regulation of the PHGDH gene in triglyceride metabolism, providing novel insights on putative intervention targets.

Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

BACKGROUND: The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples. RESULTS: We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases. CONCLUSIONS: Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome.

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.

Long-range epigenetic regulation is conferred by genetic variation located at thousands of independent loci.

The interplay between genetic and epigenetic variation is only partially understood. One form of epigenetic variation is methylation at CpG sites, which can be measured as methylation quantitative trait loci (meQTL). Here we report that in a panel of lymphocytes from 1,748 individuals, methylation levels at 1,919 CpG sites are correlated with at least one distal (trans) single-nucleotide polymorphism (SNP) (P<3.2 × 10(-13); FDR<5%). These trans-meQTLs include 1,657 SNP-CpG pairs from different chromosomes and 262 pairs from the same chromosome that are >1 Mb apart. Over 90% of these pairs are replicated (FDR<5%) in at least one of two independent data sets. Genomic loci harbouring trans-meQTLs are significantly enriched (P<0.001) for long non-coding transcripts (2.2-fold), known epigenetic regulators (2.3-fold), piwi-interacting RNA clusters (3.6-fold) and curated transcription factors (4.1-fold), including zinc-finger proteins (8.75-fold). Long-range epigenetic networks uncovered by this approach may be relevant to normal and disease states.

Genome-wide investigation of DNA methylation marks associated with FV Leiden mutation.

In order to investigate whether DNA methylation marks could contribute to the incomplete penetrance of the FV Leiden mutation, a major genetic risk factor for venous thrombosis (VT), we measured genome-wide DNA methylation levels in peripheral blood samples of 98 VT patients carrying the mutation and 251 VT patients without the mutation using the dedicated Illumina HumanMethylation450 array. The genome-wide analysis of 388,120 CpG probes identified three sites mapping to the SLC19A2 locus whose DNA methylation levels differed significantly (p<3 10-8) between carriers and non-carriers. The three sites replicated (p<2 10-7) in an independent sample of 214 individuals from five large families ascertained on VT and FV Leiden mutation among which 53 were carriers and 161 were non-carriers of the mutation. In both studies, these three CpG sites were also associated (2.33 10-110.05). In conclusion, our work clearly illustrates some promises and pitfalls of DNA methylation investigations on peripheral blood DNA in large epidemiological cohorts. DNA methylation levels at SLC19A2 are influenced by SNPs in LD with FV Leiden, but these DNA methylation marks do not explain the incomplete penetrance of the FV Leiden mutation.

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.