PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosgenPLoS GeneticsSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Genet. 2010 October; 6(10): e1001167.
Published online 2010 October 21. doi:  10.1371/journal.pgen.1001167
PMCID: PMC2958814

Genetic Association Study Identifies HSPB7 as a Risk Gene for Idiopathic Dilated Cardiomyopathy

Greg Gibson, Editor

Abstract

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−3, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10−3, 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−13, 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.

Author Summary

Dilated cardiomyopathy is a severe disease of the heart muscle and often leads to chronic heart failure, eventually with the consequence of cardiac transplantation. Identification of genetic disease markers in at-risk persons could play an important role in preventive health care. Several mutations in familial forms of the disease are described. Here, we examine the role of common genetic variants on the sporadic form of dilated cardiomyopathy. By screening about 2,000 candidate genes previously related to cardiovascular disease in more than 1,900 cases and 3,600 controls, we show that a polymorphism in the HSPB7 gene (rs1739843) is strongly associated with susceptibility to dilated cardiomyopathy. We also show that the effect on disease risk is present in both German and French cohorts. Therefore, this study is an important step towards revealing insight in the genetic background of the sporadic form of dilated cardiomyopathy.

Introduction

Dilated cardiomyopathy (DCM) is a common form of heart muscle disease with a prevalence of 1[ratio]2,500 in the general population. It represents a major cause of cardiovascular morbidity and mortality and is characterized by systolic dysfunction as well as dilation and impaired contraction of the ventricles, often leading to chronic heart failure and eventually requiring cardiac transplantation [1]. In about 35% of cases DCM is a familial disease [2]. However, in the sporadic form of DCM, i. e. after exclusion of affected family members and all detectable causes (also called idiopathic DCM), a genetic component is discussed, but can thus far not be assigned to single gene defects. Knowledge of genetic risk factors for both, familial and non-familial forms of DCM is important to initiate treatment prior to symptomatic onset of the disease, to delay its occurrence or possibly halt its progression. To date, only a few common susceptibility alleles for sporadic DCM were identified from candidate-gene approaches, but could not be confirmed in replication samples [2], [3], this being a common problem of single gene based analyses [4]. In contrast, unbiased genome-wide association studies (GWAS) allow the identification of genetic risk factors even outside of known genes, but higher power is needed to compensate for multiple testing [5]. No comprehensive GWAS was performed to date on sporadic form of DCM.

The cardiovascular gene-centric 50K single nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array represents an established compromise between GWAS and hypothesis-driven candidate gene approach by analyzing polymorphisms in more than 2,000 genes known or predicted to be involved in cardiovascular phenotypes [6].

In this study, we conducted a screening based on the cardiovascular 50K SNP array with three independent replication studies to reveal insight in genetic contribution to idiopathic DCM. The four samples from Germany and France included 1,910 sporadic DCM cases and 3,630 healthy controls individuals. We identified a common intronic variant in HSPB7, encoding a cardiovascular small heat shock protein, to be associated with sporadic form of DCM.

Results

Screening stage

In our screening case-control sample, DCM cases were more likely men, were slightly younger and less frequently smokers, had a lower BMI and a higher prevalence of hypertension, hypercholesterolemia as well as type 2 diabetes (Table 1).

Table 1
Characteristics of DCM cases and controls used for initial screening.

After quality control, 30,920 SNPs were available for analysis with 23,307 independent markers (defined as SNPs with pairwise r2<0.8 based on linkage disequilibrium (LD) in the control group). Therefore, we set a significance threshold to 0.05/23,307 = 2.15×10−6 to account for the multiple testing. In association analyses of this stage 1 study applying logistic regression adjusted for gender, four SNPs, namely rs1739843 (HSPB7, intron 2), rs11701453 (RUNX1, intron 1), rs7597774 (ADD2, intron 1) and rs2229714 (RPS6KA1, 3′ untranslated region) showed a p-value below this threshold (3.16*10−8, 1.65*10−7, 2.05*10−7, and 1.51*10−6, respectively). Results were similar when additionally adjusting for age (e.g. for rs1739843 p = 2.40*10−8). None of the four polymorphisms showed deviation from Hardy-Weinberg equilibrium. The lowest p-value for association with DCM was observed for a SNP located in HSPB7 intron 2 (rs1739843) leading to a protective effect of the minor allele (OR = 0.67 [95% CI 0.58–0.77]). Analysis of the region around the SNP rs1739843 using HapMap data (release #22) revealed the presence of six genes and 27 polymorphisms in LD with the lead SNP (r2-value>0.5) (Figure 1A). Nine of these SNPs were present on the cardiovascular 50K array after quality control and were located in HSPB7 gene as well as two genes downstream, CLCNKA and CLCNKB (Figure 1B; Table S1).

Figure 1
Linkage disequilibrium (LD) structure of HSPB7 genomic region and association results.

In this sample, the genomic inflation factor λ was 1.285 for the highest 90% of the 30,920 observed p-values. When correcting rs1739843 for this λ factor, the p-value was 1.06*10−6 and OR = 0.67 [95% CI 0.57–0.79] (Table 2).

Table 2
Association of SNPs showing p-values < 2×10−6 in the initial screening sample and follow-up in three independent replication samples analyzed by logistic regression adjusted for gender.

Replication

The four SNPs with uncorrected p<2.15×10−6 in the initial scan (rs1739843, rs11701453, rs7597774 and rs2229714) were analyzed using logistic regression adjusted for gender in three independent replication samples. First, n = 564 additional German DCM patients and n = 981 controls were genotyped for the four SNPs. Marker rs1739843 showed strong association with DCM (p = 2.07*10−3, OR = 0.79 [95% CI 0.67–0.92]). Conversely, for rs2229714 (p = 0.075, OR = 1.20 [95% CI 0.98–1.47]), rs11701453 (p = 0.373, OR = 1.09 [95% CI 0.90–1.31]) and rs7597774 (p = 0.621, OR = 1.04 [95% CI 0.89–1.22]) the initial association results were not replicated. Second, a French replication sample (France 1) consisted of n = 433 cases and n = 395 controls. Only rs1739843 showed association with DCM after adjustment for gender (p = 3.73*10−3, OR = 0.74 [95% CI 0.60–0.91]). For the other SNPs, no significant association was seen in this sample. Third, in an independent second French replication sample (France 2), again only rs1739843 showed association with DCM after adjustment for gender (p = 2.26*10−4, OR = 0.63 [95% CI 0.50–0.81]). Replication results are summarized in Table 2. None of the four polymorphisms showed deviation from Hardy-Weinberg equilibrium in any replication samples.

In a combined analysis of the screening step, corrected for the λ factor of 1.285, and the three follow-up studies (n = 5,540), the SNP rs1739843 reached a p-value of 5.28*10−13 (OR = 0.72 [95% CI 0.65–0.78]) for association with idiopathic DCM (Table 2, Figure 2). There was no between-study heterogeneity for this effect (I2 = 6.9%, p = 0.36).

Figure 2
Forest plot for rs1739843 in initial screening sample (λ-corrected) and three replication samples (Germany, France 1, and France 2), together with results from the combined analysis.

Resequencing

To reveal potential causal variants, the coding region of HSPB7 was resequenced in a total of 48 DCM patients. We detected three known synonymous variants (rs945416, rs732286 and rs1739840). The synonymous variants rs945416 (position 19, serine) and rs732286 (position 33, alanine) are in high LD with rs1739843 (r2 = 0.96, HapMap CEU data release #24). SNP rs1739840 (position 117, threonine) is not available in HapMap. In the initial sample of 664 DCM patients, all three synonymous polymorphisms are in perfect LD to each other and to rs1739843 as shown by genotyping. Neither missense nor splice site de novo mutations were identified by sequencing. Synonymous SNP rs11807575, as well as non-synonymous variants rs77021870 and rs74626772 were listed in databases, but not found to be polymorphic in our sample.

Analysis of DCM candidate genes

Since the design of the 50K human gene-centric bead chip (IBC array) aims at a large-scale gene-based approach, we screened candidate genes which are known for or potentially involved in susceptibility to DCM in our initial screening sample utilizing the information on 30,920 SNPs. We established a list of previously reported genes for DCM by searching PubMed and OMIM databases (http://www.ncbi.nlm.nih.gov/) for “CARDIOMYOPATHY, DILATED” and “GENETIC”. A total of 315 SNPs including 234 independent SNPs (defined as SNPs with pairwise r2<0.8 based on LD in the control group) were located in or near (+/−10kb) the chosen candidate genes, representing 1.01% of array content. DCM association results for these SNPs were obtained from our screening study on 664 cases and 1,874 controls (Table 3, more details in Table S2). On a single candidate gene level, polymorphisms in or near ABCC9, DES, MYH6 and TPM1 showed nominal significance after Bonferroni correction for the number of SNPs tested in gene regions (p = 0.010, p = 0.022, p = 0.005 and p = 0.018, respectively). However, none of these SNPs remained significant after correction for the 234 independent markers tested in this DCM candidate gene approach. Our study was powered to detect moderate to large effects (e. g. for OR>1.3 and MAF = 30% or OR>1.5 and MAF = 20% or OR>1.7 and MAF = 10%, the power was 56%, 96% and 97% for two-sided p<0.05/234 = 2.14*10−4, respectively).

Table 3
Candidate gene approach on DCM causing or susceptibility genes.

Discussion

In the present case-control study, we evaluated the relationship of common SNPs with sporadic DCM using a large-scale screening approach. Our comprehensive strategy set out to analyze the human gene-centric 50K bead chip (IBC array), which focuses on loci with a potential functional link to cardiovascular disease (CVD) and covers more than 45,000 SNPs from about 2,000 genes [6].

Our study identified a polymorphism (rs1739843) in intron 2 of the HSPB7 gene being associated with susceptibility to DCM in a German case-control sample with three replication steps. Recently, Cappola et al. reported an association between rs1739843 and both, ischemic and non-ischemic heart failure, applying the same gene-centric 50K bead chip [7]. They found a protective effect of the minor allele, which is in conformity with our results on DCM. As DCM is a potential preliminary stage for non-ischemic heart failure, these independent findings point to a possible common pathophysiologic cascade. However, a second association signal for heart failure located in the FRMD4B region (rs6787362, minor allele frequency (MAF) 10.4%) identified by Cappola et al. [7] could not be detected in our DCM case-control sample (p = 0.64). Our study had a power of 99% to find a nominal association between DCM and rs6787362 with p<0.05 and an OR = 0.67.

The finding on HSPB7 is also in-line with a previously reported large-scale re-sequencing approach in four biologically relevant cardiac signaling genes, which detected HSPB7 sequence diversity in sporadic cardiomyopathy [8]. Our data together with the results from Cappola et al. [7] and Matkovich et al. [8], substantiate the importance of rs1739843 or related polymorphisms in the HSPB7 locus for DCM and heart failure and possibly underscore a common genetic basis for these related phenotypes.

Matkovich et al. further report that none of the detected HSPB7 gene variants altered amino acid sequence [8], which is also consistent with the fact that we found neither missense nor splice site mutations in the HSPB7 sequence. Therefore, the biological mechanism explaining the association between the polymorphism rs1739843 and DCM risk remains still unclear. The three detected synonymous variants (rs945416, rs732286 and rs1739840) are in high LD with each other as well as with our lead SNP rs1739843 and lie on one LD block. Therefore, it could be hypothesized that these SNPs represent causal risk factors for DCM, as described for the P-glycoprotein encoding gene MCP1 and affected drug and inhibitor interactions [9]. Synonymous SNPs lead to changes in codon usage and may cause functional implications by conformational changes in protein structure due to translation efficiency. Alternatively, a de novo splice site could be created by a SNP or other (unmapped) polymorphisms outside the HSPB7 coding region may alter its gene expression. Clearly, functional studies would be required to prove these hypotheses.

Besides the HSPB7 gene, where the lead SNP is located, also five genes (CLCNKA, CLCNKB, C1orf64, ZBTB17 and SPEN) lying on the same LD block may potentially be responsible for the association with DCM. CLCNKA and CLCNKB encode for two members of the family of voltage-gated chloride channels. These proteins are predominantly expressed in the kidney and participate in renal salt reabsorption [10]. The function of C1orf64 is currently unknown. ZBTB17, also known as MIZ-1, encodes a zinc finger protein involved in the regulation of c-myc [11]. SPEN (RBM15C or MINT) encodes a conserved transcriptional repressor that controls the expression of regulators in diverse signaling pathways [12], [13].

HSPB7, encoding the small heat shock protein cvHsp (also known as HspB7), is the functionally most plausible candidate gene in this genomic region. It is known to be expressed in cardiovascular and insulin-sensitive tissues [14]. In general, the expression and activation of heat shock proteins is influenced by elevated temperatures as well as ischemia, hypoxia and acute cellular stress [15], [16]. In the aging skeletal muscle increase of cvHsp protein content was observed [17]. cvHsp was shown to be constitutively localized under non-stressful conditions to nuclear splicing speckles and may influence mRNA processing [18]. Recent data suggest co-localization between cvHsp and α-B-crystallin in the z-band of cardiac tissue and interaction with other small heat shock proteins [19]. However, further investigations like genomic fine-mapping and subgroup analyses in the context of cardiomyopathies are needed.

Genetic analyses in familial forms of DCM led to the identification of risk loci showing X-linked, autosomal dominant or autosomal recessive patterns of inheritance [2], [20], [21]. Some of the DCM causing genes or plausible candidate genes were also covered by the 50K bead chip, wherefore we specifically tested those SNPs lying in risk gene regions (10 kb upstream and downstream, respectively). In these analyses, no significant association with any of the gene variants was found, indicating that in sporadic cases of DCM probably other pathways are involved than in familial DCM. However, less frequent variants may have been missed due to insufficient power of our screening sample. Furthermore, the distinction between familial and sporadic forms of DCM is, to a certain degree, somewhat arbitrary. Screening of family members is rarely done in clinical routine, but when carried out on a systematic basis, up to 7% of previously healthy first-degree relatives have reduced left ventricular function or dilation without presence of cardiac symptoms [22]. Therefore, it might be anticipated that genetic testing could help to identify individuals at risk in familial DCM but also in families of patients affected by so-called idiopathic forms of the disease.

Already known genetic factors account for only a fraction of DCM heritability [20]. Given a 1.5-fold increased risk of DCM among heterozygous subjects in our screening sample (48% in the general population-based KORA study) and a 2.25 times increased risk among homozygous subjects (34% in KORA), 49% of DCM cases would be attributable to the SNP rs1739843 (or correlated polymorphisms) with 19% attributable to heterozygous and 30% to homozygous carriers, respectively. Therefore, the genetic component seems to comprise a large proportion for this disease. However, with the prevalence of the idiopathic form of the disease being about 1[ratio]2,700 [23], a genetic screening of the general population would include four cases out of 10,000 screened persons and two of these would have the disease due to this SNP. Therefore, the great potential of this variant might rather be screening of high risk populations, or this pathway indicates potential drug targets. Further investigations should aim (1) to identify additional variants underlying DCM susceptibility with otherwise unknown etiology and (2) to analyze potential influence of these common alleles as modifiers for familial forms of DCM. Taken together for both, modifiers of familial forms and susceptibility alleles in idiopathic DCM, knowledge of genetic background will support preventive medical measures in the future.

Some limitations of our study should be mentioned. First, we conducted a large-scale SNP analysis focused on genes potentially involved in cardiovascular traits. Therefore, on the one hand we were able to detect associations between DCM and polymorphisms only in these pre-selected genes. On the other hand, the 50K human CVD bead chip allows comprehensive gene-based analysis with more than 2,000 well covered loci. Second, our sample size only allowed to detect moderate to large effects (e. g. for OR>1.3 and MAF = 30% or OR>1.5 and MAF = 20% or OR>1.7 and MAF = 10%, the power was 19%, 75% and 80% for p<2.15*10−6, respectively). Therefore, we may have overlooked real association signals in our screening step. Third, there could be some population stratification in our initial screen sample. However, the observed λ could also be caused - in part - by underlying association due to the analysis of pre-selected loci known or suggested to be involved in cardiovascular phenotypes. The fact that the association between rs1739843 in HSPB7 and idiopathic DCM was replicated in three independent samples strongly enhances the confidence in our results.

Materials and Methods

Ethics statement

The ethics committees of the participating study centers approved the study protocol and all participants gave their written informed consent. The study was in accordance with the principles of the current version of the Declaration of Helsinki.

Case-control samples and phenotyping

Cases for the initial German screening study were recruited from the German Heart Institute (Berlin), and controls were from a population-based German KORA study (follow-up survey F3, Augsburg) [24]. Phenotypic details are summarized in Table 1. Controls (n = 1,874) had no medical history for coronary artery disease (CAD), myocardial infarction or DCM; mean age was 62±11 years and slightly more women (n = 986) than men (n = 888) were present in the control group. Inclusion criteria for DCM cases were the following: reduced systolic function (left ventricular ejection fraction (LVEF) <45%) without angiographically assessed evidence of major CAD, significant valvular heart disease (>grade 2, i. e. such as mitral or aortic regurgitation), hypertensive heart disease, congenital heart disease, myocarditis (by endomyocardial biopsy, when available) or other secondary forms of heart failure. Patients with a positive family history were also excluded from this study. In DCM cases (n = 664), mean LVEF was 24±3% and mean age of disease diagnosis was 46±11 years.

For the first replication step, additional German DCM cases (mean age 53±13 years; n = 564, n = 440 men, n = 124 women) were recruited from different German study centers: Berlin, n = 64; Lübeck, n = 96 (Angio-Lueb); Marburg, n = 61 (EUROGENE); Münster, n = 101 (EUROGENE); Regensburg, n = 150 (EUROGENE); Regensburg, n = 92 (GoKard). Independent German KORA controls from surveys S1 and S2 (n = 981, n = 539 men, n = 442 women) had a mean age of 52±10 years [24]. Inclusion and exclusion criteria were identical to the initial case-control sample.

A second replication study (France 1) was recruited in France (CARDIGENE) [25], [26]. The French cases were of white European origin (all born in France, from parents born in France or neighboring countries) with a diagnosis of DCM, i. e. enlarged left ventricle end-diastolic volume/diameter >140 ml/m2 on ventriculography or >34 mm/m2 on echocardiography and LVEF ≤40% confirmed over a six-month period, in the absence of causal factors such as CAD or sustained hypertension, intrinsic valvular disease, documented myocarditis, congenital malformation, insulin-dependent diabetes. Only apparently sporadic DCM cases without additional (first degree) relative with DCM were included (but 8% were in fact with familial form after careful cardiac examination in relatives). Recruitment was performed in ten hospitals in six regions in France (Lille, Lyon, Nancy, Nantes, Paris-Ile de France, Strasbourg) from September 1994 to February 1996. A total of 433 patients (229 had undergone a cardiac transplantation) were included (n = 345 men, n = 88 women). Mean age of patients was 45±11 years, mean LVEF was 23±7% and mean end-diastolic volume was 195±67 ml/m2. Controls (n = 395) were age- and gender-matched (n = 310 men, n = 85 women).

The third replication sample was also of French origin (France 2). Inclusion criteria were identical to the France 1 sample. A total of 249 patients from EUROGENE and PHRC were included (n = 198 men, n = 51 women). Mean age of patients at diagnosis was 51±10 years. Controls (n = 380) were free of medical history for CAD, myocardial infarction or DCM and mean age was 46±11 years (n = 301 men, n = 79 women).

Genotyping

Initial genotyping was carried out using the 50K gene-centric human CVD bead chip version 1 (IBC v1 array) (Illumina, San Diego, CA, USA) [6] following the manufacturer's protocol. Data were analyzed (calling and sample clustering) and exported employing BeadStudio analysis software (Illumina). From the initial 45,707 SNPs, those markers with low call rates (<95%) or low frequency (MAF<1%) were excluded. Minimal call rate per individual was 90%. We used identity-by-descent methods to exclude unknown first-degree relation of participants.

Replication samples were taken from human CVD bead chip data or genotyped with 5′ exonuclease TaqMan technology (Applied Biosystems, Foster City, CA, USA) as previously described [27]. A by-design assay for rs1739843 was used with primer sequences 5′-CTCTGCCATCACCATCTCACA-3′ and 5′-GGCAGAGGGAGCCTGAG-3′ and probe sequences 5′-VIC-AGGGTGGGAGGTGACAG-NFQ-3′ and 5′-FAM-AGGGTGGGAGATGACAG-NFQ-3′ (site of rs1739843 is underlined; fluorescence dyes VIC and FAM on 5′ end and non-fluorescence quencher (NFQ) on 3′ end are indicated). All other assays were obtained pre-designed directly from Applied Biosystems. Detailed information on assays used in France 2 sample are available at http://genecanvas.ecgene.net/infusions/genecanvas/Polymorphisms/PolymorphismsList.php.

SNP rs1739843 was re-genotyped using the by-design TaqMan assay in initial case sample (n = 664) to check for discrepancies between human CVD bead chip and TaqMan genotypes. A >99.8% concordance of genotypes was found. For all genotyped samples a call rate >97% for each SNP assay was reached.

Resequencing

Polymerase chain reaction (PCR) primer were generated using Primer3Plus (http://www.bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi) [28] to cover the coding parts of the three HSPB7 exons (GenBank accession No. NM_14424.4). The primer sequences and PCR amplification products are listed in Table S3. Included intronic regions were 267 bp for 5′ end of intron 1, 156 bp for 3′ end of intron 1, 136 bp for 5′ end of intron 2, and 89 bp for 3′ end of intron 2, respectively. PCR cycling conditions consisted of an initial denaturation at 95°C for 9 min, followed by 40 cycles with denaturation at 95°C for 30 s, annealing at 60°C for 30 s, and elongation at 72°C for 30 s, with a final elongation step at 72°C for 7 min.

After PCR amplification, primers and dNTPs were removed using ExoSAP-IT (USB Europe, Staufen, Germany) following the manufacturer's instructions. The purified PCR products were directly sequenced using the ABI PRISM BigDye Terminator Cycle Sequencing Ready Reaction Kit Version 3.1 on the ABI 3730 (Applied Biosystems, Foster City, CA, USA).

Statistical analyses

For initial screening and replication analyses, logistic regression adjusted for gender was used. P-values, odds ratios (OR) and their 95% confidence intervals (CI) were reported. The inflation factor λ was computed in the 50K initial screening analysis for logistic regression analysis assuming a χ2 distribution with two degrees of freedom of the minus two-times logep measures (90% highest p-values). The p-values and CI from initial screening analysis were genomic-control corrected using this λ factor via standard errors (standard error[corrected] = sqrt(λ)*standard error) and beta estimates (95%CI beta[corrected] = beta±1.96*standard error[corrected]). Deviation from Hardy-Weinberg equilibrium was calculated with an exact test [29]. Statistical and association analyses were performed using JMP 7.0.2 (SAS Institute Inc, Cary, NC, USA) and PLINK v1.07 (http://pngu.mgh.harvard.edu/~purcell/plink/) [30], respectively. Power analysis was carried out using Quanto 1.2.4 (http://hydra.usc.edu/gxe/). We combined the initial scan results corrected for λ with the replication studies' results using a fixed effect model. Annotation of association results on a genome level was performed with WGAViewer software (http://people.genome.duke.edu/~dg48/WGAViewer/) [31]. LD patterns were calculated using HapMap releases #22 and #24 (http://www.hapmap.org/) [32].

Supporting Information

Table S1

Association results of SNPs in HSPB7 genomic region in initial screening sample (664 cases, 1,874 controls).

(0.05 MB DOC)

Table S2

Association of previously reported DCM causing or susceptibility genes in initial screening sample (664 cases, 1,874 controls).

(0.49 MB DOC)

Table S3

Primer sequences for PCR amplification of the HSPB7 coding region and product sizes.

(0.03 MB DOC)

Acknowledgments

We appreciate the invaluable contribution of participants of the studies. We gratefully acknowledge the excellent technical assistance of Anja Angelov, Françoise Gary, Claire Perret, Michaela Schmid, and Josef Simon. We thank Peter Oefner and Corinna Feuchtinger for sequencing support.

The KORA study group consists of H.-E. Wichmann (speaker), R. Holle, J. John, T. Illig, C. Meisinger, A. Peters, and their coworkers, who are responsible for the design and conduct of the KORA studies.

The members of the EUROGENE Heart Failure project consisted of the following persons: Coordinator: Richard Isnard (Paris, France). France: Michel Komajda, Philippe Charron, Marie-Claude Aumont (Paris), Olivier Dubourg (Boulogne). Germany: Christian Hengstenberg (Regensburg), Anette Richter, Bernhard Maisch (Marburg), Thomas Wichter (Münster). Italy: Alessandra Serio, Eloisa Arbustini, Luigi Tavazzi (Pavia). Portugal: Dulce Brito, Hugo Madeira (Lisbon). Spain: Ana Diaz, Enrique Galve (Barcelona). Sweden: Stellan Morner, Anders Waldenström (Umea). UK: John G.F. Cleland (Kingston-upon-Hull).

The PHRC French network (Programme Hospitalier de Recherche Clinique AOM04141) consisted of the following persons: Jean-Noël Trochu (Nantes), Laurent Fauchier (Tours), Pascal DeGroote (Lille), Olivier Dubourg (Boulogne), Marie-Claude Aumont (Paris-Bichat), Pierre Gibelin (Nice), Jean F. Aupetit (Lyon), Michel Komajda, Richard Isnard, Philippe Charron (all Paris-Salpêtrière).

Footnotes

The authors have declared that no competing interests exist.

Part of this work was financed by the German National Genome Research Network (NGFN-2 01GS0417 and NGFNplus Atherogenomics 01GS0832) and by the University Hospital of Regensburg, Germany. Part of this work was supported by grants from Assistance Publique - Hopitaux de Paris [PHRC: AOM04141], the Fondation LEDUCQ (Eurogene Heart Failure network and CAERUS network), the Societe Francaise de Cardiologie/Federation Francaise de Cardiologie (grant in 2006). The KORA research platform (KORA, Cooperative Research in the Region of Augsburg) was initiated and financed by the Helmholtz Center Munich - German Research Center for Environmental Health, which is funded by the German Federal Ministry of Education and Research (BMBF) and by the State of Bavaria. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Maron BJ, Towbin JA, Thiene G, Antzelevitch C, Corrado D, et al. Contemporary definitions and classification of the cardiomyopathies: an American Heart Association Scientific Statement from the Council on Clinical Cardiology, Heart Failure and Transplantation Committee; Quality of Care and Outcomes Research and Functional Genomics and Translational Biology Interdisciplinary Working Groups; and Council on Epidemiology and Prevention. Circulation. 2006;113:1807–1816. [PubMed]
2. Richard P, Villard E, Charron P, Isnard R. The Genetic Bases of Cardiomyopathies. Journal of the American College of Cardiology. 2006;48:A79–A89.
3. Rampersaud E, Kinnamon DD, Hamilton K, Khuri S, Hershberger RE, et al. Common Susceptibility Variants Examined for Association with Dilated Cardiomyopathy. Ann Hum Genet. 2010;74:110–116. [PMC free article] [PubMed]
4. Ioannidis JP. Why most published research findings are false. PLoS Med. 2005;2:e124. doi: 10.1371/journal.pmed.0020124. [PubMed]
5. Pearson TA, Manolio TA. How to interpret a genome-wide association study. JAMA. 2008;299:1335–1344. [PubMed]
6. Keating BJ, Tischfield S, Murray SS, Bhangale T, Price TS, et al. Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies. PLoS ONE. 2008;3:e3583. doi: 10.1371/journal.pone.0003583. [PMC free article] [PubMed]
7. Cappola TP, Li M, He J, Ky B, Gilmore J, et al. Common variants in HSPB7 and FRMD4B associated with advanced heart failure. Circ Cardiovasc Genet. 2010;3:147–154. [PMC free article] [PubMed]
8. Matkovich SJ, Van Booven DJ, Hindes A, Kang MY, Druley TE, et al. Cardiac signaling genes exhibit unexpected sequence diversity in sporadic cardiomyopathy, revealing HSPB7 polymorphisms associated with disease. J Clin Invest. 2010;120:280–289. [PMC free article] [PubMed]
9. Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, et al. A “silent” polymorphism in the MDR1 gene changes substrate specificity. Science. 2007;315:525–528. [PubMed]
10. Uchida S. In vivo role of CLC chloride channels in the kidney. Am J Physiol Renal Physiol. 2000;279:F802–F808. [PubMed]
11. Ikegaki N, Gotoh T, Kung B, Riceberg JS, Kim DY, et al. De novo identification of MIZ-1 (ZBTB17) encoding a MYC-interacting zinc-finger protein as a new favorable neuroblastoma gene. Clin Cancer Res. 2007;13:6001–6009. [PubMed]
12. Newberry EP, Latifi T, Towler DA. The RRM domain of MINT, a novel Msx2 binding protein, recognizes and regulates the rat osteocalcin promoter. Biochemistry. 1999;38:10678–10690. [PubMed]
13. Li J, Wang J, Yang X, Li J, Qin H, et al. The Spen homolog Msx2-interacting nuclear target protein interacts with the E2 ubiquitin-conjugating enzyme UbcH8. Mol Cell Biochem. 2006;288:151–157. [PubMed]
14. Krief S, Faivre JF, Robert P, Le DB, Brument-Larignon N, et al. Identification and characterization of cvHsp. A novel human small stress protein selectively expressed in cardiovascular and insulin-sensitive tissues. J Biol Chem. 1999;274:36592–36600. [PubMed]
15. Lindquist S, Craig EA. The heat-shock proteins. Annu Rev Genet. 1988;22:631–677. [PubMed]
16. Anckar J, Sistonen L. Heat shock factor 1 as a coordinator of stress and developmental pathways. Adv Exp Med Biol. 2007;594:78–88. [PubMed]
17. Doran P, Gannon J, O'Connell K, Ohlendieck K. Aging skeletal muscle shows a drastic increase in the small heat shock proteins alphaB-crystallin/HspB5 and cvHsp/HspB7. Eur J Cell Biol. 2007;86:629–640. [PubMed]
18. Vos MJ, Kanon B, Kampinga HH. HSPB7 is a SC35 speckle resident small heat shock protein. Biochim Biophys Acta. 2009;1793:1343–1353. [PubMed]
19. Brodehl A, Martin I, Gawlowski T, Stork I, Gummert J, et al. The small heat shock proteins cvHSP and alpha-B-Crystallin are colocalized and interact in human myocardial tissue. Clin Res Cardiol. 2010;99:P467.
20. Burkett EL, Hershberger RE. Clinical and genetic issues in familial dilated cardiomyopathy. J Am Coll Cardiol. 2005;45:969–981. [PubMed]
21. Osterziel KJ, Hassfeld S, Geier C, Perrot A. [Familial dilated cardiomyopathy]. Herz. 2005;30:529–534. [PubMed]
22. Michels VV, Olson TM, Miller FA, Ballman KV, Rosales AG, et al. Frequency of development of idiopathic dilated cardiomyopathy among relatives of patients with idiopathic dilated cardiomyopathy. Am J Cardiol. 2003;91:1389–1392. [PubMed]
23. Codd MB, Sugrue DD, Gersh BJ, Melton LJ., III Epidemiology of idiopathic dilated and hypertrophic cardiomyopathy. A population-based study in Olmsted County, Minnesota, 1975–1984. Circulation. 1989;80:564–572. [PubMed]
24. Wichmann HE, Gieger C, Illig T. KORA-gen–resource for population genetics, controls and a broad spectrum of disease phenotypes. Gesundheitswesen. 2005;67(Suppl 1):S26–S30. [PubMed]
25. Tesson F, Charron P, Peuchmaurd M, Nicaud V, Cambien F, et al. Characterization of a unique genetic variant in the beta1-adrenoceptor gene and evaluation of its role in idiopathic dilated cardiomyopathy. CARDIGENE Group. J Mol Cell Cardiol. 1999;31:1025–1032. [PubMed]
26. Charron P, Tesson F, Poirier O, Nicaud V, Peuchmaurd M, et al. Identification of a genetic risk factor for idiopathic dilated cardiomyopathy. Involvement of a polymorphism in the endothelin receptor type A gene. CARDIGENE group. Eur Heart J. 1999;20:1587–1591. [PubMed]
27. Stark K, Reinhard W, Neureuther K, Wiedmann S, Sedlacek K, et al. Association of common polymorphisms in GLUT9 gene with gout but not with coronary artery disease in a large case-control study. PLoS ONE. 2008;3:e1948. doi: 10.1371/journal.pone.0001948. [PMC free article] [PubMed]
28. Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, et al. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35:W71–W74. [PMC free article] [PubMed]
29. Wigginton JE, Cutler DJ, Abecasis GR. A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet. 2005;76:887–893. [PubMed]
30. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
31. Ge D, Zhang K, Need AC, Martin O, Fellay J, et al. WGAViewer: software for genomic annotation of whole genome association studies. Genome Res. 2008;18:640–643. [PubMed]
32. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–861. [PMC free article] [PubMed]
33. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, et al. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–2939. [PMC free article] [PubMed]
34. Bienengraeber M, Olson TM, Selivanov VA, Kathmann EC, O'Cochlain F, et al. ABCC9 mutations identified in human dilated cardiomyopathy disrupt catalytic KATP channel gating. Nat Genet. 2004;36:382–387. [PMC free article] [PubMed]
35. Olson TM, Michels VV, Thibodeau SN, Tai YS, Keating MT. Actin mutations in dilated cardiomyopathy, a heritable form of heart failure. Science. 1998;280:750–752. [PubMed]
36. Mohapatra B, Jimenez S, Lin JH, Bowles KR, Coveler KJ, et al. Mutations in the muscle LIM protein and alpha-actinin-2 genes in dilated cardiomyopathy and endocardial fibroelastosis. Mol Genet Metab. 2003;80:207–215. [PubMed]
37. Duboscq-Bidot L, Charron P, Ruppert V, Fauchier L, Richter A, et al. Mutations in the ANKRD1 gene encoding CARP are responsible for human dilated cardiomyopathy. Eur Heart J. 2009;30:2128–2136. [PubMed]
38. Moulik M, Vatta M, Witt SH, Arola AM, Murphy RT, et al. ANKRD1, the gene encoding cardiac ankyrin repeat protein, is a novel dilated cardiomyopathy gene. J Am Coll Cardiol. 2009;54:325–333. [PMC free article] [PubMed]
39. Knoll R, Hoshijima M, Hoffman HM, Person V, Lorenzen-Schmidt I, et al. The cardiac mechanical stretch sensor machinery involves a Z disc complex that is defective in a subset of human dilated cardiomyopathy. Cell. 2002;111:943–955. [PubMed]
40. Li D, Tapscoft T, Gonzalez O, Burch PE, Quinones MA, et al. Desmin mutation responsible for idiopathic dilated cardiomyopathy. Circulation. 1999;100:461–464. [PubMed]
41. Muntoni F, Cau M, Ganau A, Congiu R, Arvedi G, et al. Brief report: deletion of the dystrophin muscle-promoter region associated with X-linked dilated cardiomyopathy. N Engl J Med. 1993;329:921–925. [PubMed]
42. Towbin JA, Hejtmancik JF, Brink P, Gelb B, Zhu XM, et al. X-linked dilated cardiomyopathy. Molecular genetic evidence of linkage to the Duchenne muscular dystrophy (dystrophin) gene at the Xp21 locus. Circulation. 1993;87:1854–1865. [PubMed]
43. Posch MG, Posch MJ, Geier C, Erdmann B, Mueller W, et al. A missense variant in desmoglein-2 predisposes to dilated cardiomyopathy. Mol Genet Metab. 2008;95:74–80. [PubMed]
44. Schönberger J, Wang L, Shin JT, Kim SD, Depreux FF, et al. Mutation in the transcriptional coactivator EYA4 causes dilated cardiomyopathy and sensorineural hearing loss. Nat Genet. 2005;37:418–422. [PubMed]
45. Arimura T, Hayashi YK, Murakami T, Oya Y, Funabe S, et al. Mutational analysis of fukutin gene in dilated cardiomyopathy and hypertrophic cardiomyopathy. Circ J. 2009;73:158–161. [PubMed]
46. Friedrichs F, Zugck C, Rauch GJ, Ivandic B, Weichenhan D, et al. HBEGF, SRA1, and IK: Three cosegregating genes as determinants of cardiomyopathy. Genome Res. 2009;19:395–403. [PubMed]
47. Sylvius N, Tesson F, Gayet C, Charron P, Benaiche A, et al. A new locus for autosomal dominant dilated cardiomyopathy identified on chromosome 6q12-q16. Am J Hum Genet. 2001;68:241–246. [PubMed]
48. Fatkin D, MacRae C, Sasaki T, Wolff MR, Porcu M, et al. Missense mutations in the rod domain of the lamin A/C gene as causes of dilated cardiomyopathy and conduction-system disease. N Engl J Med. 1999;341:1715–1724. [PubMed]
49. Daehmlow S, Erdmann J, Knueppel T, Gille C, Froemmel C, et al. Novel mutations in sarcomeric protein genes in dilated cardiomyopathy. Biochem Biophys Res Commun. 2002;298:116–120. [PubMed]
50. Carniel E, Taylor MR, Sinagra G, Di LA, Ku L, et al. Alpha-myosin heavy chain: a sarcomeric gene associated with dilated and hypertrophic phenotypes of cardiomyopathy. Circulation. 2005;112:54–59. [PubMed]
51. Kamisago M, Sharma SD, DePalma SR, Solomon S, Sharma P, et al. Mutations in sarcomere protein genes as a cause of dilated cardiomyopathy. N Engl J Med. 2000;343:1688–1696. [PubMed]
52. Arimura T, Nakamura T, Hiroi S, Satoh M, Takahashi M, et al. Characterization of the human nebulette gene: a polymorphism in an actin-binding motif is associated with nonfamilial idiopathic dilated cardiomyopathy. Hum Genet. 2000;107:440–451. [PubMed]
53. Hassel D, Dahme T, Erdmann J, Meder B, Huge A, et al. Nexilin mutations destabilize cardiac Z-disks and lead to dilated cardiomyopathy. Nat Med. 2009;15:1281–1288. [PubMed]
54. Haghighi K, Kolokathis F, Pater L, Lynch RA, Asahi M, et al. Human phospholamban null results in lethal dilated cardiomyopathy revealing a critical difference between mouse and human. J Clin Invest. 2003;111:869–876. [PMC free article] [PubMed]
55. Schmitt JP, Kamisago M, Asahi M, Li GH, Ahmad F, et al. Dilated cardiomyopathy and heart failure caused by a mutation in phospholamban. Science. 2003;299:1410–1413. [PubMed]
56. Li D, Parks SB, Kushner JD, Nauman D, Burgess D, et al. Mutations of presenilin genes in dilated cardiomyopathy and heart failure. Am J Hum Genet. 2006;79:1030–1039. [PubMed]
57. Brauch KM, Karst ML, Herron KJ, de AM, Pellikka PA, et al. Mutations in ribonucleic acid binding protein gene cause familial dilated cardiomyopathy. J Am Coll Cardiol. 2009;54:930–941. [PMC free article] [PubMed]
58. Olson TM, Michels VV, Ballew JD, Reyna SP, Karst ML, et al. Sodium channel mutations and susceptibility to heart failure and atrial fibrillation. JAMA. 2005;293:447–454. [PMC free article] [PubMed]
59. Tsubata S, Bowles KR, Vatta M, Zintz C, Titus J, et al. Mutations in the human delta-sarcoglycan gene in familial and sporadic dilated cardiomyopathy. J Clin Invest. 2000;106:655–662. [PMC free article] [PubMed]
60. Bione S, D'Adamo P, Maestrini E, Gedeon AK, Bolhuis PA, et al. A novel X-linked gene, G4.5. is responsible for Barth syndrome. Nat Genet. 1996;12:385–389. [PubMed]
61. D'Adamo P, Fassone L, Gedeon A, Janssen EA, Bione S, et al. The X-linked gene G4.5 is responsible for different infantile dilated cardiomyopathies. Am J Hum Genet. 1997;61:862–867. [PubMed]
62. Hayashi T, Arimura T, Itoh-Satoh M, Ueda K, Hohda S, et al. Tcap gene mutations in hypertrophic cardiomyopathy and dilated cardiomyopathy. J Am Coll Cardiol. 2004;44:2192–2201. [PubMed]
63. Taylor MR, Slavov D, Gajewski A, Vlcek S, Ku L, et al. Thymopoietin (lamina-associated polypeptide 2) gene mutation associated with dilated cardiomyopathy. Hum Mutat. 2005;26:566–574. [PubMed]
64. Mogensen J, Murphy RT, Shaw T, Bahl A, Redwood C, et al. Severe disease expression of cardiac troponin C and T mutations in patients with idiopathic dilated cardiomyopathy. J Am Coll Cardiol. 2004;44:2033–2040. [PubMed]
65. Murphy RT, Mogensen J, Shaw A, Kubo T, Hughes S, et al. Novel mutation in cardiac troponin I in recessive idiopathic dilated cardiomyopathy. Lancet. 2004;363:371–372. [PubMed]
66. Olson TM, Kishimoto NY, Whitby FG, Michels VV. Mutations that alter the surface charge of alpha-tropomyosin are associated with dilated cardiomyopathy. J Mol Cell Cardiol. 2001;33:723–732. [PubMed]
67. Gerull B, Gramlich M, Atherton J, McNabb M, Trombitas K, et al. Mutations of TTN, encoding the giant muscle filament titin, cause familial dilated cardiomyopathy. Nat Genet. 2002;30:201–204. [PubMed]
68. Olson TM, Illenberger S, Kishimoto NY, Huttelmaier S, Keating MT, et al. Metavinculin mutations alter actin interaction in dilated cardiomyopathy. Circulation. 2002;105:431–437. [PubMed]

Articles from PLoS Genetics are provided here courtesy of Public Library of Science