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.].

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.

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.

Dysfunctional nitric oxide signalling increases risk of myocardial infarction.

Myocardial infarction, a leading cause of death in the Western world, usually occurs when the fibrous cap overlying an atherosclerotic plaque in a coronary artery ruptures. The resulting exposure of blood to the atherosclerotic material then triggers thrombus formation, which occludes the artery. The importance of genetic predisposition to coronary artery disease and myocardial infarction is best documented by the predictive value of a positive family history. Next-generation sequencing in families with several affected individuals has revolutionized mutation identification. Here we report the segregation of two private, heterozygous mutations in two functionally related genes, GUCY1A3 (p.Leu163Phefs*24) and CCT7 (p.Ser525Leu), in an extended myocardial infarction family. GUCY1A3 encodes the α1 subunit of soluble guanylyl cyclase (α1-sGC), and CCT7 encodes CCTη, a member of the tailless complex polypeptide 1 ring complex, which, among other functions, stabilizes soluble guanylyl cyclase. After stimulation with nitric oxide, soluble guanylyl cyclase generates cGMP, which induces vasodilation and inhibits platelet activation. We demonstrate in vitro that mutations in both GUCY1A3 and CCT7 severely reduce α1-sGC as well as ß1-sGC protein content, and impair soluble guanylyl cyclase activity. Moreover, platelets from digenic mutation carriers contained less soluble guanylyl cyclase protein and consequently displayed reduced nitric-oxide-induced cGMP formation. Mice deficient in α1-sGC protein displayed accelerated thrombus formation in the microcirculation after local trauma. Starting with a severely affected family, we have identified a link between impaired soluble-guanylyl-cyclase-dependent nitric oxide signalling and myocardial infarction risk, possibly through accelerated thrombus formation. Reversing this defect may provide a new therapeutic target for reducing the risk of myocardial infarction.

Expression pattern in human macrophages dependent on 9p21.3 coronary artery disease risk locus.

OBJECTIVE: Genome-wide association studies identified a risk haplotype on chromosome 9p21.3 to be associated with coronary artery disease (CAD) and myocardial infarction (MI). Since this region does not contain a clear candidate gene with known pathophysiology, we performed a haplotype-specific expression study in human macrophages during pro-inflammatory stimulation to investigate the locus-dependent expression patterns in a model of atherosclerosis, the underlying cause for CAD and MI. METHODS: Blood samples were taken from 40 male stable MI patients either homozygous for 9p21.3 risk (n = 20) or non-risk haplotype (n = 20) as well as from 28 healthy male individuals (n = 14 for each haplotype). Monocytes were isolated by density gradient centrifugation followed by differentiation into macrophages via M-CSF. Macrophages were either incubated with a pro-inflammatory IFNγ-LPS cocktail or kept untreated as controls. After 24 h, RNA was isolated and applied to Affymetrix Human Exon 1.0 ST Arrays. RESULTS: Macrophages from MI patients and controls stratified for 9p21.3 haplotypes, exhibited marked differences in gene expression. Most pronounced differences were found in inflammatory mediators, like the chemokines CCL8 and CCL2 and the lectines CLEC4E and CLEC5A. Differences in expression changes could be seen most obviously during inflammatory stimulation for both, the interleukins IL12B and IL1B, and members of metallothionein gene family. CONCLUSION: These findings show that gene expression is different in 9p21.3 haplotype-stratified macrophages. While these effects are relatively small in our in vitro model of atherosclerosis, these biological effects may contribute to a long term effect in risk haplotype carriers increasing susceptibility to CAD and MI.

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06 × 10⁻6, OR  = 0.67 [95% CI 0.57-0.79] for the minor allele T). Three more SNPs showed p < 2.21 × 10⁻5. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n =564, n = 981 controls, p = 2.07 × 10⁻³, OR = 0.79 [95% CI 0.67-0.92]), France 1 (n = 433 cases, n = 395 controls, p =3.73 × 10⁻³, OR  = 0.74 [95% CI 0.60-0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26 × 10⁻4, OR  = 0.63 [95% CI 0.50-0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28 × 10⁻¹³, OR= 0.72 [95% CI 0.65-0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.