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.

Genome-wide association analysis in dilated cardiomyopathy reveals two new players in systolic heart failure on chromosomes 3p25.1 and 22q11.23

AIMS

Our objective was to better understand the genetic bases of dilated cardiomyopathy (DCM), a leading cause of systolic heart failure.

METHODS AND RESULTS

We conducted the largest genome-wide association study performed so far in DCM, with 2719 cases and 4440 controls in the discovery population. We identified and replicated two new DCM-associated loci on chromosome 3p25.1 [lead single-nucleotide polymorphism (SNP) rs62232870, P = 8.7 × 10-11 and 7.7 × 10-4 in the discovery and replication steps, respectively] and chromosome 22q11.23 (lead SNP rs7284877, P = 3.3 × 10-8 and 1.4 × 10-3 in the discovery and replication steps, respectively), while confirming two previously identified DCM loci on chromosomes 10 and 1, BAG3 and HSPB7. A genetic risk score constructed from the number of risk alleles at these four DCM loci revealed a 3-fold increased risk of DCM for individuals with 8 risk alleles compared to individuals with 5 risk alleles (median of the referral population). In silico annotation and functional 4C-sequencing analyses on iPSC-derived cardiomyocytes identify SLC6A6 as the most likely DCM gene at the 3p25.1 locus. This gene encodes a taurine transporter whose involvement in myocardial dysfunction and DCM is supported by numerous observations in humans and animals. At the 22q11.23 locus, in silico and data mining annotations, and to a lesser extent functional analysis, strongly suggest SMARCB1 as the candidate culprit gene.

CONCLUSION

This study provides a better understanding of the genetic architecture of DCM and sheds light on novel biological pathways underlying heart failure.

Genome wide association analysis in dilated cardiomyopathy revealed two new susceptibility loci for systolic heart failure

 

We conducted the largest Genome Wide Association Study performed so far in Dilated Cardiomyopathy (DCM), a leading cause of systolic heart failure and cardiovascular death. Using a discovery phase of 2,719 cases and 4,440 controls and a replication phase of 584 independent cases and 966 controls, we identified and replicated two new DCM-associated loci, one on chromosome 3p (meta-analysis OR = 1.36 [1.25 - 1.48], p=5.3 10–13) and the second on chromosome 22q (meta-analysis OR = 1.33 [1.22 - 1.46], p=5.0 10–10), while confirming the two previously identified DCM loci on chromosome 10 and 1, BAG3 and HSPB7. We estimated the global heritability to 31% ± 8%. The genetic risk score constructed from the number of lead risk-alleles at these 4 loci revealed a 27% risk increased in individuals with 8 risk-alleles compared to the 5 risk alleles reference group (OR = 1.27 [1.14–1.42]).

The two association signals were then fine-mapped by combining in silico and functional genomics investigations (as 4C-sequencing on iPSC-derived cardiomyocytes). While a few genes remain candidates at the second locus and deserve further investigations, our work clearly identified one gene as responsible for the association at the first locus whose role in the pathophysiology of DCM is supported by recent observations in human and mice. As the biological pathway in which this gene is involved is a potential target for pharmacological agents, our finding opens novel therapeutic perspectives for treating or preventing heart failure.

We are convinced that these results provide new findings that add both on the understanding of the genetic architecture of heart failure and on potential new players involved in the pathophysiology of this devastating disease.

Funding Acknowledgement

Type of funding source: Public grant(s) – National budget only. Main funding source(s): GENMED Laboratory of Excellence on Medical Genomics, DETECTIN-HF project (ERA-CVD framework), Assistance Publique-Hôpitaux de Paris, Délégation à la recherche clinique, the “Fondation LEDUCQ”, the PROMEX charitable foundation, the Société Française de Cardiologie/Fédération Française de Cardiologie, the Deutsche Forschungsgemeinschaft, The Federal Ministry of Education and Research and the Ministry of Cultural Affairs and the Social Ministry of the Federal State of Mecklenburg-West Pomerania, The German Center for Cardiovascular Research (DZHK), Hospitals NIHR Biomedical Research Centre, NWO VENI grant (no. 016.176.136)

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 (R²=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.

The Leu33/Pro Polymorphism (PIA1/PIA2 ) of the Glycoprotein IIIa (GPIIIa) Receptor Is Not Related to Myocardial Infarction in the ECTIM Study

The GPIIb/IIIa receptor complex may contribute to acute coronary syndromes by mediating platelet aggregation. The Leu33/Pro polymorphism (PlAl/PlA2) of the GPIIIa has recently been shown to be associated with CHD in a small case-control study. We have investigated this polymorphism in a large multicenter study of patients with myocardial infarction and controls and found no difference in the distribution of allele and genotype frequencies between cases and controls.

The 4G/5G Genetic Polymorphism in the Promoter of the Plasminogen Activator Inhibitor-1 (PAI-1) Gene Is Associated with Differences in Plasma PAI-1 Activity but not with Risk of Myocardial Infarction in the ECTIM Study

We have investigated the interrelationships of plasma PAI-1 activity, the PAI-1 4G/5G polymorphism and risk of myocardial infarction (MI) in the ECTIM study, a case-control study of MI based in Belfast, Lille, Strasbourg and Toulouse. Mean PAI-1 levels in cases were similar across all centres but in controls, levels in the French centres were significantly higher. Only in Belfast were PAI-1 levels higher in cases (11.7AU/ml) than controls (10.5AU/ml). The PAI-1 4G allele frequency was similar in cases and controls (0.55 and 0.54). In all groups, 4G homozygotes had the highest mean plasma PAI-1 level (4G4G vs 5G5G; cases overall: 14.2 vs 12.1 AU/ml; controls overall: 15.0 vs 12.6AU/ml), with the heterozygotes generally intermediate. The data from Belfast are consistent with the literature implicating PAI-1 level as an MI risk factor. In ECTIM, the PAI-1 4G/5G polymorphism is not a genetic risk factor for MI but is associated with PAI-1 activity. Thus homozygosity for the 4G allele may predispose to elevated PAI-1 and impaired fibrinolysis, perhaps requiring interaction with other genetic or environmental factors to influence MI risk.

Expression of a Paternal History of Premature Myocardial Infarction on Fibrinogen, Factor VIIC and PAI-1 in European Offspring - The EARS Study

To assess the role of genetic and environmental factors in the predisposition to atherosclerosis, 682 students whose father had suffered a myocardial infarction before the age of 55 (“cases”) and 1312 controls matched for age and sex, were recruited from 14 Universities in Europe. Fibrinogen, factor Vile and PAI-1 were compared between cases and controls across European regions. Fibrinogen and factor Vile were positively correlated with BMI, smoking and contraception. PAI-1 was positively and independently correlated with BMI and waist-to-hip ratio, and negatively with contraception. Factor Vile and PAI-1 were correlated with cholesterol and triglycerides, and fibrinogen was weakly correlated with LDL-cholesterol. After adjustment for covariates and lipids, fibrinogen level was significantly higher in male cases than in controls (2.38 vs 2.29, p <0.01). No such difference was found in females (2.59 vs 2.57 - NS). There was no significant case/control difference for factor Vile and PAI-1. These results support the hypothesis that fibrinogen is a transmissible risk factor of coronary artery disease in males. No such evidence was provided for factor Vile and PAI-1.

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.

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.

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.

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.

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

First case of a human RASGRP2 mutation affecting Rap1 activation in platelets and causing 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.

Genome-wide haplotype analysis of cis expression quantitative trait loci in monocytes

In order to assess whether gene expression variability could be influenced by several SNPs acting in cis, either through additive or more complex haplotype effects, a systematic genome-wide search for cis haplotype expression quantitative trait loci (eQTL) was conducted in a sample of 758 individuals, part of the Cardiogenics Transcriptomic Study, for which genome-wide monocyte expression and GWAS data were available. 19,805 RNA probes were assessed for cis haplotypic regulation through investigation of ~2,1 × 10(9) haplotypic combinations. 2,650 probes demonstrated haplotypic p-values >10(4)-fold smaller than the best single SNP p-value. Replication of significant haplotype effects were tested for 412 probes for which SNPs (or proxies) that defined the detected haplotypes were available in the Gutenberg Health Study composed of 1,374 individuals. At the Bonferroni correction level of 1.2 × 10(-4) (~0.05/412), 193 haplotypic signals replicated. 1000 G imputation was then conducted, and 105 haplotypic signals still remained more informative than imputed SNPs. In-depth analysis of these 105 cis eQTL revealed that at 76 loci genetic associations were compatible with additive effects of several SNPs, while for the 29 remaining regions data could be compatible with a more complex haplotypic pattern. As 24 of the 105 cis eQTL have previously been reported to be disease-associated loci, this work highlights the need for conducting haplotype-based and 1000 G imputed cis eQTL analysis before commencing functional studies at disease-associated loci.

Graphical modeling of gene expression in monocytes suggests molecular mechanisms explaining increased atherosclerosis in smokers

Smoking is a risk factor for atherosclerosis with reported widespread effects on gene expression in circulating blood cells. We hypothesized that a molecular signature mediating the relation between smoking and atherosclerosis may be found in the transcriptome of circulating monocytes. Genome-wide expression profiles and counts of atherosclerotic plaques in carotid arteries were collected in 248 smokers and 688 non-smokers from the general population. Patterns of co-expressed genes were identified by Independent Component Analysis (ICA) and network structure of the pattern-specific gene modules was inferred by the PC-algorithm. A likelihood-based causality test was implemented to select patterns that fit models containing a path “smoking→gene expression→plaques”. Robustness of the causal inference was assessed by bootstrapping. At a FDR ≤0.10, 3,368 genes were associated to smoking or plaques, of which 93% were associated to smoking only. SASH1 showed the strongest association to smoking and PPARG the strongest association to plaques. Twenty-nine gene patterns were identified by ICA. Modules containing SASH1 and PPARG did not show evidence for the “smoking→gene expression→plaques” causality model. Conversely, three modules had good support for causal effects and exhibited a network topology consistent with gene expression mediating the relation between smoking and plaques. The network with the strongest support for causal effects was connected to plaques through SLC39A8, a gene with known association to HDL-cholesterol and cellular uptake of cadmium from tobacco, while smoking was directly connected to GAS6, a gene reported to have anti-inflammatory effects in atherosclerosis and to be up-regulated in the placenta of women smoking during pregnancy. Our analysis of the transcriptome of monocytes recovered genes relevant for association to smoking and atherosclerosis, and connected genes that before, were only studied in separate contexts. Inspection of correlation structure revealed candidates that would be missed by expression-phenotype association analysis alone.

Powerful identification of cis-regulatory SNPs in human primary monocytes using allele-specific gene expression

A large number of genome-wide association studies have been performed during the past five years to identify associations between SNPs and human complex diseases and traits. The assignment of a functional role for the identified disease-associated SNP is not straight-forward. Genome-wide expression quantitative trait locus (eQTL) analysis is frequently used as the initial step to define a function while allele-specific gene expression (ASE) analysis has not yet gained a wide-spread use in disease mapping studies. We compared the power to identify cis-acting regulatory SNPs (cis-rSNPs) by genome-wide allele-specific gene expression (ASE) analysis with that of traditional expression quantitative trait locus (eQTL) mapping. Our study included 395 healthy blood donors for whom global gene expression profiles in circulating monocytes were determined by Illumina BeadArrays. ASE was assessed in a subset of these monocytes from 188 donors by quantitative genotyping of mRNA using a genome-wide panel of SNP markers. The performance of the two methods for detecting cis-rSNPs was evaluated by comparing associations between SNP genotypes and gene expression levels in sample sets of varying size. We found that up to 8-fold more samples are required for eQTL mapping to reach the same statistical power as that obtained by ASE analysis for the same rSNPs. The performance of ASE is insensitive to SNPs with low minor allele frequencies and detects a larger number of significantly associated rSNPs using the same sample size as eQTL mapping. An unequivocal conclusion from our comparison is that ASE analysis is more sensitive for detecting cis-rSNPs than standard eQTL mapping. Our study shows the potential of ASE mapping in tissue samples and primary cells which are difficult to obtain in large numbers.

Cytokine phenotype, genotype, and renal outcomes at cardiac surgery

BACKGROUND

Cardiac surgery modulates pro- and anti-inflammatory cytokine balance involving plasma tumour necrosis factor alpha (TNFα) and interleukin-10 (IL-10) together with urinary transforming growth factor beta-1 (TGFβ1), interleukin-1 receptor antagonist (IL1ra) and tumour necrosis factor soluble receptor-2 (TNFsr2). Effects on post-operative renal function are unclear. We investigated if following cardiac surgery there is a relationship between cytokine (a) phenotype and renal outcome; (b) genotype and phenotype and (c) genotype and renal outcome. Since angiotensin-2 (AG2), modulates TGFβ1 production, we determined whether angiotensin converting enzyme insertion/deletion (ACE I/D) genotype affects urinary TGFβ1 phenotype as well as renal outcome.

METHODS

In 408 elective cardiac surgery patients we measured pre- and 24 h post-operative urinary TGFβ-1, IL1ra and TNFsr2 and pre- and 2 h post-operative plasma TNFα and IL-10. Post-operative responses were compared for each cytokine in patients grouped according to presence or absence of renal dysfunction defined as a drop from baseline eGFR of greater than 25% (as calculated by the method of modification of diet in renal disease (MDRD)) occurring (1) within the first 24 and (2) 48 postoperative hours (early renal dysfunction), (3) on the fifth postoperative day (late renal dysfunction) or (4) at any time throughout the 5 day postoperative period (early and late combined). Patient genotype was determined for TNF/G-308A, TGFβ1-509 C/T, IL10/G-1082A and ACE I/D.

RESULTS

Post-operative plasma IL-10 and urinary TGFβ1 responses were significantly higher in patients who developed early renal dysfunction. IL1ra and TNFsr2 responses were significantly lower 24h post-operatively in patients who developed late renal dysfunction. Genotype did not alter cytokine phenotype or outcome. CONCLUSIONS/INFERENCES: Cytokine profiling may help predict early and late renal dysfunction. Genotypes studied did not alter phenotype or outcome.

Comprehensive exploration of the effects of miRNA SNPs on monocyte gene expression

We aimed to assess whether pri-miRNA SNPs (miSNPs) could influence monocyte gene expression, either through marginal association or by interacting with polymorphisms located in 3'UTR regions (3utrSNPs). We then conducted a genome-wide search for marginal miSNPs effects and pairwise miSNPs × 3utrSNPs interactions in a sample of 1,467 individuals for which genome-wide monocyte expression and genotype data were available. Statistical associations that survived multiple testing correction were tested for replication in an independent sample of 758 individuals with both monocyte gene expression and genotype data. In both studies, the hsa-mir-1279 rs1463335 was found to modulate in cis the expression of LYZ and in trans the expression of CNTN6, CTRC, COPZ2, KRT9, LRRFIP1, NOD1, PCDHA6, ST5 and TRAF3IP2 genes, supporting the role of hsa-mir-1279 as a regulator of several genes in monocytes. In addition, we identified two robust miSNPs × 3utrSNPs interactions, one involving HLA-DPB1 rs1042448 and hsa-mir-219-1 rs107822, the second the H1F0 rs1894644 and hsa-mir-659 rs5750504, modulating the expression of the associated genes.

As some of the aforementioned genes have previously been reported to reside at disease-associated loci, our findings provide novel arguments supporting the hypothesis that the genetic variability of miRNAs could also contribute to the susceptibility to human diseases.

Statistical colocalization of monocyte gene expression and genetic risk variants for type 1 diabetes

One mechanism by which disease-associated DNA variation can alter disease risk is altering gene expression. However, linkage disequilibrium (LD) between variants, mostly single-nucleotide polymorphisms (SNPs), means it is not sufficient to show that a particular variant associates with both disease and expression, as there could be two distinct causal variants in LD. Here, we describe a formal statistical test of colocalization and apply it to type 1 diabetes (T1D)-associated regions identified mostly through genome-wide association studies and expression quantitative trait loci (eQTLs) discovered in a recently determined large monocyte expression data set from the Gutenberg Health Study (1370 individuals), with confirmation sought in an additional data set from the Cardiogenics Transcriptome Study (558 individuals). We excluded 39 out of 60 overlapping eQTLs in 49 T1D regions from possible colocalization and identified 21 coincident eQTLs, representing 21 genes in 14 distinct T1D regions. Our results reflect the importance of monocyte (and their derivatives, macrophage and dendritic cell) gene expression in human T1D and support the candidacy of several genes as causal factors in autoimmune pancreatic beta-cell destruction, including AFF3, CD226, CLECL1, DEXI, FKRP, PRKD2, RNLS, SMARCE1 and SUOX, in addition to the recently described GPR183 (EBI2) gene.

Integrating genome-wide genetic variations and monocyte expression data reveals trans-regulated gene modules in humans

One major expectation from the transcriptome in humans is to characterize the biological basis of associations identified by genome-wide association studies. So far, few cis expression quantitative trait loci (eQTLs) have been reliably related to disease susceptibility. Trans-regulating mechanisms may play a more prominent role in disease susceptibility. We analyzed 12,808 genes detected in at least 5% of circulating monocyte samples from a population-based sample of 1,490 European unrelated subjects. We applied a method of extraction of expression patterns-independent component analysis-to identify sets of co-regulated genes. These patterns were then related to 675,350 SNPs to identify major trans-acting regulators. We detected three genomic regions significantly associated with co-regulated gene modules. Association of these loci with multiple expression traits was replicated in Cardiogenics, an independent study in which expression profiles of monocytes were available in 758 subjects. The locus 12q13 (lead SNP rs11171739), previously identified as a type 1 diabetes locus, was associated with a pattern including two cis eQTLs, RPS26 and SUOX, and 5 trans eQTLs, one of which (MADCAM1) is a potential candidate for mediating T1D susceptibility. The locus 12q24 (lead SNP rs653178), which has demonstrated extensive disease pleiotropy, including type 1 diabetes, hypertension, and celiac disease, was associated to a pattern strongly correlating to blood pressure level. The strongest trans eQTL in this pattern was CRIP1, a known marker of cellular proliferation in cancer. The locus 12q15 (lead SNP rs11177644) was associated with a pattern driven by two cis eQTLs, LYZ and YEATS4, and including 34 trans eQTLs, several of them tumor-related genes. This study shows that a method exploiting the structure of co-expressions among genes can help identify genomic regions involved in trans regulation of sets of genes and can provide clues for understanding the mechanisms linking genome-wide association loci to disease.

New gene functions in megakaryopoiesis and platelet formation

Platelets are the second most abundant cell type in blood and are essential for maintaining haemostasis. Their count and volume are tightly controlled within narrow physiological ranges, but there is only limited understanding of the molecular processes controlling both traits. Here we carried out a high-powered meta-analysis of genome-wide association studies (GWAS) in up to 66,867 individuals of European ancestry, followed by extensive biological and functional assessment. We identified 68 genomic loci reliably associated with platelet count and volume mapping to established and putative novel regulators of megakaryopoiesis and platelet formation. These genes show megakaryocyte-specific gene expression patterns and extensive network connectivity. Using gene silencing in Danio rerio and Drosophila melanogaster, we identified 11 of the genes as novel regulators of blood cell formation. Taken together, our findings advance understanding of novel gene functions controlling fate-determining events during megakaryopoiesis and platelet formation, providing a new example of successful translation of GWAS to function.

EDNRB gene variants and melanoma risk in two southern European populations

BACKGROUND

EDNRB gene variants were reported to be associated with melanoma risk in French patients, with the S305N variant showing the highest frequency.

AIM

To verify the S305N association with melanoma risk in an independent larger French population (378 patients, 389 controls); to investigate the role of EDNRB variants in melanoma risk in an Italian population (133 patients, 118 controls); and to explore the association of CDKN2A or CDK4 mutations with the S305N EDNRB variant in a subgroup of patients (59 French, 12 Italian) with a suspected hereditary predisposition to melanoma (familial melanoma, sporadic multiple primary melanoma or melanoma associated with pancreatic cancer).

METHODS

The S305N variant was genotyped in the French population, while the EDNRB gene in the Italian population was entirely sequenced.

RESULTS

Overall, there was no significant difference in the frequency of the S305N variant between patients with sporadic melanoma and controls in either the French or the Italian population. However, a significantly higher S305N allele frequency was detected in French patients with a suspected hereditary predisposition to melanoma compared with controls (P = 0.04). In addition, in this subgroup of patients, the S305N allele was also significantly associated with the presence of CDKN2A mutations (P = 0.04).

CONCLUSIONS

Our results showed no evidence of association of the S305N EDNRB polymorphism with sporadic melanoma risk in either the French or Italian populations, but there was an indication that EDNRB might be a melanoma-predisposing gene in French patients with a suspected hereditary predisposition to melanoma.

A genome-wide association study identifies two loci associated with heart failure due to dilated cardiomyopathy

AIMS

Dilated cardiomyopathy (DCM) is a major cause of heart failure with a high familial recurrence risk. So far, the genetics of DCM remains largely unresolved. We conducted the first genome-wide association study (GWAS) to identify loci contributing to sporadic DCM.

METHODS AND RESULTS

One thousand one hundred and seventy-nine DCM patients and 1108 controls contributed to the discovery phase. Pools of DNA stratified on disease status, population, age, and gender were constituted and used for testing association of DCM with 517 382 single nucleotide polymorphisms (SNPs). Three DCM-associated SNPs were confirmed by individual genotyping (P < 5.0 10(-7)), and two of them, rs10927875 and rs2234962, were replicated in independent samples (1165 DCM patients and 1302 controls), with P-values of 0.002 and 0.009, respectively. rs10927875 maps to a region on chromosome 1p36.13 which encompasses several genes among which HSPB7 has been formerly suggested to be implicated in DCM. The second identified locus involves rs2234962, a non-synonymous SNP (c.T757C, p. C151R) located within the sequence of BAG3 on chromosome 10q26. To assess whether coding mutations of BAG3 might cause monogenic forms of the disease, we sequenced BAG3 exons in 168 independent index cases diagnosed with familial DCM and identified four truncating and two missense mutations. Each mutation was heterozygous, present in all genotyped relatives affected by the disease and absent in a control group of 347 healthy individuals, strongly suggesting that these mutations are causing the disease.

CONCLUSION

This GWAS identified two loci involved in sporadic DCM, one of them probably implicates BAG3. Our results show that rare mutations in BAG3 contribute to monogenic forms of the disease, while common variant(s) in the same gene are implicated in sporadic DCM.

Multiple marker approach to risk stratification in patients with stable coronary artery disease

AIMS

multimarker approaches for risk prediction in coronary artery disease have remained inconsistent. We assessed multiple biomarkers representing distinct pathophysiological pathways in relation to cardiovascular events in stable angina.

METHODS AND RESULTS

we investigated 12 biomarkers reflecting inflammation [C-reactive protein, growth-differentiation factor (GDF)-15, neopterin], lipid metabolism (apolipoproteins AI, B100), renal function (cystatin C, serum creatinine), and cardiovascular function and remodelling [copeptin, C-terminal-pro-endothelin-1, mid-regional-pro-adrenomedullin (MR-proADM), mid-regional-pro-atrial natriuretic peptide (MR-proANP), N-terminal-pro-B-type natriuretic peptide (Nt-proBNP)] in 1781 stable angina patients in relation to non-fatal myocardial infarction and cardiovascular death (n = 137) over 3.6 years. Using Cox proportional hazards models and C-indices, the strongest association with outcome for log-transformed biomarkers in multivariable-adjusted analyses was observed for Nt-proBNP [hazard ratio (HR) for one standard deviation increase 1.65, 95% confidence interval (CI) 1.28-2.13, C-index 0.686], GDF-15 (HR 1.59, 95% CI 1.25-2.02, C-index 0.681), MR-proANP (HR 1.46, 95% CI 1.14-1.87, C-index 0.673), cystatin C (HR 1.39, 95% CI 1.10-1.75, C-index 0.671), and MR-proADM (HR 1.63, 95% CI 1.21-2.20, C-index 0.668). Each of these top single markers and their combination (C-index 0.690) added predictive information beyond the baseline model consisting of the classical risk factors assessed by C-index and led to substantial reclassification (P-integrated discrimination improvement <0.05).

CONCLUSION

comparative analysis of 12 biomarkers revealed Nt-proBNP, GDF-15, MR-proANP, cystatin C, and MR-proADM as the strongest predictors of cardiovascular outcome in stable angina. All five biomarkers taken separately offered incremental predictive ability over established risk factors. Combination of the single markers slightly improved model fit but did not enhance risk prediction from a clinical perspective.

Genetics and beyond--the transcriptome of human monocytes and disease susceptibility

Background

Variability of gene expression in human may link gene sequence variability and phenotypes; however, non-genetic variations, alone or in combination with genetics, may also influence expression traits and have a critical role in physiological and disease processes.

Methodology/Principal Findings

To get better insight into the overall variability of gene expression, we assessed the transcriptome of circulating monocytes, a key cell involved in immunity-related diseases and atherosclerosis, in 1,490 unrelated individuals and investigated its association with >675,000 SNPs and 10 common cardiovascular risk factors. Out of 12,808 expressed genes, 2,745 expression quantitative trait loci were detected (P<5.78×10⁻¹²), most of them (90%) being cis-modulated. Extensive analyses showed that associations identified by genome-wide association studies of lipids, body mass index or blood pressure were rarely compatible with a mediation by monocyte expression level at the locus. At a study-wide level (P<3.9×10⁻⁷), 1,662 expression traits (13.0%) were significantly associated with at least one risk factor. Genome-wide interaction analyses suggested that genetic variability and risk factors mostly acted additively on gene expression. Because of the structure of correlation among expression traits, the variability of risk factors could be characterized by a limited set of independent gene expressions which may have biological and clinical relevance. For example expression traits associated with cigarette smoking were more strongly associated with carotid atherosclerosis than smoking itself.

Conclusions/Significance

This study demonstrates that the monocyte transcriptome is a potent integrator of genetic and non-genetic influences of relevance for disease pathophysiology and risk assessment.

Phospholipolyzed LDL induces an inflammatory response in endothelial cells through endoplasmic reticulum stress signaling

Secreted phospholipases A2 (sPLA2s) are present in atherosclerotic plaques and are now considered novel attractive therapeutic targets and potential biomarkers as they contribute to the development of atherosclerosis through lipoprotein-dependent and independent mechanisms. We have previously shown that hGX-sPLA2-phospholipolyzed LDL (LDL-X) induces proinflammatory responses in human umbilical endothelial cells (HUVECs); here we explore the molecular mechanisms involved. Global transcriptional gene expression profiling of the response of endothelial cells exposed to either LDL or LDL-X revealed that LDL-X activates multiple distinct cellular pathways including the unfolded protein response (UPR). Mechanistic insight showed that LDL-X activates UPR through calcium depletion of intracellular stores, which in turn disturbs cytoskeleton organization. Treatment of HUVECs and aortic endothelial cells (HAECs) with LDL-X led to activation of all 3 proximal initiators of UPR: eIF-2alpha, IRE1alpha, and ATF6. In parallel, we observed a sustained phosphorylation of the p38 pathway resulting in the phosphorylation of AP-1 downstream targets. This was accompanied by significant production of the proinflammatory cytokines IL-6 and IL-8. Our study demonstrates that phospholipolyzed LDL uses a range of molecular pathways including UPR to initiate endothelial cell perturbation and thus provides an LDL oxidation-independent mechanism for the initiation of vascular inflammation in atherosclerosis.

A follow-up study of a genome-wide association scan identifies a susceptibility locus for venous thrombosis on chromosome 6p24.1

To identify genetic susceptibility factors conferring increased risk of venous thrombosis (VT), we conducted a multistage study, following results of a previously published GWAS that failed to detect loci for developing VT. Using a collection of 5862 cases with VT and 7112 healthy controls, we identified the HIVEP1 locus on chromosome 6p24.1 as a susceptibility locus for VT. Indeed, the HIVEP1 rs169713C allele was associated with an increased risk for VT, with an odds ratio of 1.20 (95% confidence interval 1.13-1.27, p = 2.86 x 10(-9)). HIVEP1 codes for a protein that participates in the transcriptional regulation of inflammatory target genes by binding specific DNA sequences in their promoter and enhancer regions. The current results provide the identification of a locus involved in VT susceptibility that lies outside the traditional coagulation/fibrinolysis pathway.

Contribution of SELP and PSGL-1 genotypes and haplotypes to the presence of coronary heart disease in Tunisians

P-selectin (SELP) and its counter-receptor, P-selectin glycoprotein ligand-1 (PSGL-1), play key role in the transient attachment of leukocytes to endothelial cells predisposing to coronary heart disease (CHD). In the current report, 293 angiographically proven CHD patients and 327 age, gender, and race-matched controls were included. Our aim was to evaluate the contribution to CHD of the following SNPs: C-2123G, G-1969A and T715P in SELP, Met62Ile and the VNTR variants in PSGL-1 gene in a North African population from Tunisia. While there were no significant differences in the distribution of SELP or PSGL-1 alleles or genotypes between patients and controls, a trend for a significant association of the C-2123G genotypes distribution with incident CHD was observed (P=0.06). Assuming an additive model of transmission, the risk was 74% higher among subjects carrying the GG genotypes in comparison to those carrying the CC genotype (OR=1.74 [1.01-2.98], P=0.04) and 80% higher in the recessive model (OR=1.80 [1.08-3.01], P=0.02). Haplotype analysis did not identify any specific SELP or PSGL-1 haplotypes to be associated with CHD. The present study demonstrated no evidence of association between individual SELP or PSGL-1 SNPs or haplotypes with incident CHD. However, this study replicates absence of association of the mostly studied SNP, T715P, previously reported in individuals with African origin.

Genome-wide association study of intracranial aneurysm identifies three new risk loci

Saccular intracranial aneurysms (IAs) are balloon-like dilations of the intracranial arterial wall; their hemorrhage commonly results in severe neurologic impairment and death. We report a second genome-wide association study with discovery and replication cohorts from Europe and Japan comprising 5,891 cases and 14,181 controls with ∼832,000 genotyped and imputed SNPs across discovery cohorts. We identified three new loci showing strong evidence for association with IA in the combined data set, including intervals near RBBP8 on 18q11.2 (OR=1.22, P=1.1×10^-12), STARD13/KL on 13q13.1 (OR=1.20, P=2.5×10^-9) and a gene-rich region on 10q24.32 (OR=1.29, P=1.2×10^-9). We also confirmed prior associations near SOX17 (8q11.23-q12.1; OR=1.28, P=1.3×10^-12) and CDKN2A/B (9p21.3; OR=1.31, P=1.5×10^-22). It is noteworthy that several putative risk genes play a role in cell-cycle progression, potentially affecting proliferation and senescence of progenitor cell populations that are responsible for vascular formation and repair.

Association of three polymorphisms selected from a genome-wide association study with coronary heart disease in the Tunisian population

Despite extensive exploration of many genes, strong evidence of a molecular genetic association with coronary heart disease (CHD) or myocardial infarction (MI) remains to be obtained. Recently, significant interest has emerged in mapping genetic susceptibility for complex traits through whole-genome studies association generating promoting data that will determine the genetic contribution to common human diseases such as coronary heart disease. The aim of the present case-control study including 324 healthy controls and 296 patients with coronary heart disease from Tunisia, was to assess relation between three polymorphisms previously reported to be strongly associated with coronary heart disease in the Welcome Trust Case Control Consortium (WTCCC) and the German myocardial infarction family studies: locus 9p21.3 (rs 1333049), locus 6q25.1 (rs6922269) and 2q36.3 (rs2943634). By single locus analysis, no differences in genotype distribution and allelic frequency were found between the two groups of study. The risk allele (C) for rs2943634 was less frequent among Tunisian population than in controls from the WTCCC and German studies (57% vs 65%). The three SNPs previously reported to be associated with CHD were not replicated in our small sample.