Skip to main content

Left ventricular diastolic function associated with common genetic variation in ATP12Ain a general population

Abstract

Background

Left ventricular (LV) function depends on the activity of transmembrane electrolyte transporters. Failing human myocardium has lower Na+/K+ ATPase expression and higher intracellular sodium concentrations. The ATP12A gene encodes a catalytic subunit of an ATPase that can function as a Na+/K+ pump. We, therefore, investigated the association between LV function and common genetic variants in ATP12A.

Methods

A random sample of 1166 participants (53.7% women; mean age 49.5 years, 44.8% hypertensive) was recruited in Belgium, Poland, Italy and Russia. We measured transmitral early and late diastolic velocities (E and A) by pulsed wave Doppler, and mitral annular velocities (e' and a') by tissue Doppler. Using principal component analysis, we summarized 7 Doppler indexes - namely, E, A, e' and a' velocities, and their ratios (E/A, e'/a', and E/e') - into a single diastolic score. We genotyped 5 tag SNPs (rs963984, rs9553395, rs10507337, rs12872010, rs2071490) in ATP12A. In our analysis we focused on rs10507337 because it is located within a transcription factor binding site.

Results

In the population-based analyses while adjusting for covariables and accounting for family clusters and country, rs10507337 C allele carriers had significantly higher E/A (P = 0.003), e' (P = 5.8×10-5), e'/a' (P = 0.003) and diastolic score (P = 0.0001) compared to TT homozygotes. Our findings were confirmed in the haplotype analysis and in the family-based analyses in 74 informative offspring.

Conclusions

LV diastolic function as assessed by conventional and tissue Doppler indexes including a composite diastolic score was associated with genetic variation in ATP12A. Further experimental studies are necessary to clarify the role of ATP12A in myocardial relaxation.

Peer Review reports

Background

Left ventricular (LV) diastolic dysfunction refers to a condition in which abnormalities in LV function are present during diastole and is characterized by impaired relaxation and/or filling of the heart. LV diastolic dysfunction is associated with common risk factors such as hypertension and it can progress to symptomatic heart failure [1]. The prevalence of subclinical diastolic dysfunction in the general population increases with age and is as high as 25.1% [2]-[5]. Echocardiographic measurement of diastolic Doppler velocities of the transmitral blood flow (E and A peaks) and the mitral annular movement (e' and a' peaks) during early and late diastole allows the non-invasive assessment of LV diastolic function [6]. Moreover, low e' velocity measured by Tissue Doppler Imaging (TDI) significantly and independently from other cardiovascular risk factors predicted higher risk of fatal and nonfatal cardiovascular events in patients with hypertension [7] or heart failure [8] and in general population [9]. Recently we demonstrated significant heritability of the diastolic Doppler indexes [10]. Understanding to what extent genetic factors along with anthropometric, hemodynamic factors and lifestyle influence diastolic Doppler indexes is an important issue in view of the relation of LV diastolic dysfunction with outcome.

Cardiomyocyte contraction and relaxation depend on the balance of electrolytes (calcium, sodium, potassium, etc.) across the cellular membranes [11],[12]. The electrolyte gradients are maintained by transmembrane channels and adenosine triphosphate (ATP) dependent pumps. Experimental [13] and clinical studies [14] demonstrated that intracellular sodium concentration (Na+) is increased in failing cardiomyoctes as compared to normal myocardium. The possible cause for this observation might be related to changes in Na+/K+ ATPase expression and/or function [11],[12]. Among the genes that encode the catalytic alpha subunit of Na+/K+-ATPase [15],[16], ATP12A belongs to the family of P-type cation transport ATPases that can function as Na+/K+ ATPase in human cells [17]. Taking together, genetic variability in ATP12A might influence myocardial Na+ handling and consequently myocardial function. We recently genotyped top SNPs that might be associated with cardiovascular phenotypes and that have been identified in a recent genome-wide association study [18]. Among these SNPs common genetic variants in ATP12A were genotyped. Therefore, in the Flemish Study on Environment, Genes and Health Outcomes (FLEMENGHO) and the European Project On Genes in Hypertension (EPOGH) we investigated whether echocardiographic variables reflecting LV diastolic function are associated with common genetic variants in ATP12A.

Methods

Study participants

From August 1985 until December 2005, we randomly recruited a family-based population sample (the FLEMENGHO cohort) from a geographically defined area in northern Belgium as described in previous publications [4]. EPOGH recruited nuclear families from 1999 until 2001. The EPOGH investigators were trained at the Studies Coordinating Centre in Leuven, Belgium, and applied the same protocol, questionnaires and procedures, as used in FLEMENGHO [3]. All study participants provided a signed informed consent and the study was approved by the Ethical Committees of the University of Leuven, the University of Padova, Jagiellonian University Medical College, and the Novosibirsk Institute of Internal and Preventive Medicine. All clinical investigations were conducted according to the principles expressed in the Declaration of Helsinki. The initial response rate at enrolment was 75.0%.

In the FLEMENGHO study, from May 2005 until January 2010, we invited 1055 former participants from a previously identified random population for a follow-up examination at the field centre, including echocardiography. Of those, 828 renewed their written consent. In the EPOGH study, from January 2007 until September 2009, we invited 631 former participants for an examination including echocardiography. Of those, 561 gave their consent in writing. We excluded 223 subjects from analysis, because DNA was missing or of bad quality (n = 217), atrial fibrillation (n = 1) or because LV diastolic Doppler indexes could not be reliably measured (n = 5). Thus, 1166 subjects were analyzed including 777 FLEMENGHO participants (Noorderkempen, Belgium) and 389 EPOGH participants from Krakow, Poland (n = 143), Mirano, Italy (n = 109) and Novosibirsk, Russia (n = 137).

Echocardiography

The participants refrained from smoking, heavy exercise, and drinking alcohol or caffeine-containing beverages for at least 3 hours before echocardiography.

Data acquisition

In each center one experienced physician did the ultrasound examination, using a Vivid7 Pro (GE Vingmed, Horten, Norway) interfaced with a 2.5- 3.5 MHz phased-array probe, according to a standardized protocol as published elsewhere [3],[4]. With the subjects in partial left decubitus and breathing normally, the observer obtained images, together with a simultaneous ECG signal, from the parasternal long and short axes and from the apical 4- and 2-chamber long-axis views. M-mode echocardiograms of the LV were recorded from the parasternal long-axis view under control of the 2D image. The ultrasound beam was positioned just below the mitral valve at the level of the posterior chordae tendineae. To record mitral flow velocities from the apical window, we positioned the Doppler sample volume at the mitral valve tips.

Using TDI, the observer recorded low-velocity, high-intensity myocardial velocity at a high frame rate (>190 FPS), while adjusting the imaging angle to ensure parallel alignment of the ultrasound beam with the myocardial segment of interest. From the apical window, a 5 mm Doppler sample was placed at the septal, lateral, inferior and posterior sites of the mitral annulus to obtain the pulsed wave TDI velocities.

Off-line analysis

All echocardiographic recordings included at least 5 cardiac cycles and were digitally stored for off-line analysis. One experienced observer (TK) analyzed the digitally stored images from all centers, using the EchoPac software, version 4.0.4 (GE Vingmed, Horten, Norway), selecting and averaging 3 cardiac cycles obtained at the end-expiration. The observer was blinded to the genetic results. LV internal diameter and interventricular septal and posterior wall thickness were measured at end-diastole from the 2 dimensionally guided M-mode tracing according to the guidelines. When optimal orientation of M-mode ultrasound beam could not be obtained (in approximately 5% of echocardiographic recordings), the reader performed linear measurements on correctly oriented two-dimensional images. End-diastolic left ventricular dimensions were used to calculate LV mass by an anatomically validated formula. LV mass was indexed to body surface area (BSA). LV end-systolic and end-diastolic volumes were measured off-line using the standard biplane Simpson's method.

We assessed LV diastolic function using recordings of conventional blood flow and tissue Doppler velocities. Pulsed-wave Doppler signals of transmitral blood flow were used to measure peak early (E) and late (A) diastolic velocities. From the pulsed wave TDI recordings, we measured the early (e') and late (a') peak diastolic velocities of the mitral annulus displacement, and the e'/a' ratio at the 4 acquisition sites. We calculated the E/e' ratio by dividing transmitral E peak by e' averaged from the 4 acquisition sites. As reported previously [3], the inter-observer intra-session reproducibility across the four sampling sites ranged from 4.48% to 5.34% for e' velocities and from 3.96% to 4.52% for a' velocities.

Other measurements

We administered a standardized questionnaire to collect detailed information on subjects' medical history, smoking and drinking habits, and intake of medications. The conventional blood pressure was the average of five consecutive auscultatory readings obtained with the subject in the seated position. Hypertension was defined as a blood pressure of at least 140 mm Hg systolic or 90 mm Hg diastolic or as use of antihypertensive drugs. Body mass index was weight in kilograms divided by the square of height in meters.

Genotypes

We extracted DNA from white blood cells. For genotyping of ATP12A SNPs we used 15 K Illumina Infinium custom chip (Illumina Inc, San Diego, CA), designed for the HYPERGENES project [18] to cover genes deemed to be relevant for hypertension and related target organ damage. Detailed information on the procedure of tag SNP selection for the custom chip is provided in Additional file 1. The human ATP12A gene lies within chromosome 13q12.1-q12.3 (http://www.ncbi.nlm.nih.gov/gene/479) and spans approximately 32 Kb. It contains 23 exons and 22 introns [19]. The selected 8 tag SNPs were annotated according to the Genome Reference Consortium Human Build 37 hg19 (Additional file 1: Table S1). Three SNPs were excluded from the statistical analysis because of minor allele frequency (MAF) = 0.02 (rs2289909), genotyping call-rate <99% (rs1867767) or bad cluster according to Illumina Genome Studio (rs1001806, Additional file 1: Table S1). Of remaining 5 SNPs, rs963984, rs12872010 and rs2071490 lie within coding regions, but do not alter the amino acid sequence of the coded protein (synonymous). Two other SNPs (rs9553395 and rs10507337) are located in the 5' flanking region, with rs10507337 within a transcription factor binding site.

Statistical methods

For database management and statistical analysis, we used SAS software, version 9.3 (SAS Institute, Cary, NC). We compared means and proportions by the z test and by the χ2 test, respectively. We tested minimum allele frequency, proportion of missing genotypes and Hardy-Weinberg equilibrium using JMP Genomics, version 6.1 (SAS Institute, Cary, NC). We reconstructed haplotypes using the PROC HAPLOTYPE procedure implemented in the genetics module of the SAS software.

Principal component analysis

The diastolic Doppler velocities are measured using similar technique (Doppler) and are highly intra-correlated. Therefore, some redundancy in the cumulative variance of the Doppler velocities exists. This redundancy can be decreased by a mathematical procedure, the principal component analysis. This thechnique enables summarization of intra-correlated variables into an artificial score that carries information from all the contributing variables. Individual values of the artificial score are calculated by summarization of optimally weighted values of the contributing measured indexes. We, therefore, summarized transmitral E and A blood flow velocities, averaged TDI e' and a' mitral annulus velocities and their ratios (E/A, e'/a', E/e') into a single composite diastolic score. Normality of diastolic score distribution was evaluated by Shapiro-Wilk's statistic and skewness by computation of the coefficient of skewness. The association between the composite diastolic score and velocities was described using factors loadings. We performed stepwise linear regression to identify correlates of the composite diastolic score with sex, age, height, weight, body mass index, heart rate, systolic and diastolic blood pressures, left ventricular mass index, and ejection fraction. We set the P value for variables to enter and stay in the regression models at 0.10.

Population-based analysis

First, we performed the association analyses of dependent variables (LV phenotypes) with the genotypes or haplotypes of interest by use of a mixed model in each examination center. Covariables with known relevance for LV structure and function [3] were included in the models as fixed effects, while family clusters was modeled as a random effect. Then we performed the association analyses in all centers combined while examination centers were modeled as random effects. We tested whether observed effect sizes are homogeneous across examination centers. The heterogeneity test statistic (Cochran's Q) was computed by summing the squared deviations of each center estimate from the overall meta-analytic estimate, weighting each center contribution in the same manner as in the meta-analysis. We also reported the P-values corrected for multiple testing using the Bonferroni method.

Family-based analysis

We performed the transmission disequilibrium test for quantitative traits (QTDT). We evaluated the within- and between-family components of phenotypic variance using the orthogonal model as implemented by Abecasis et al.[20] in the QTDT software (version 2.6.1; http://www.sph.umich.edu/csg/abecasis/QTDT). In this model the between-family component is sensitive to population structure whereas the within-family component is significant in the presence of transmission disequilibrium.

Results

Characteristics of participants by center, ATP12A alleles and genotypes and haplotypes frequencies

Table'1 summarizes baseline characteristics of participants by examination center. Overall, the 1166 white European participants included 626 (53.7%) women and 522 (44.8%) hypertensive patients of whom 315 (27.0%) were on antihypertensive drug treatment. Mean age (±SD) was 49.5 ± 15.2 years. Polish participants were younger, had higher systolic blood pressure and less frequently reported a daily alcohol consumption of ≥5 gram than participants from other centers (Table'1). Composite diastolic Doppler score were highest in the Polish center compared to the remaining centers (Table'1). Additional file 1: Table S2 presents baseline characteristics of participants by sex.

Table 1 Characteristics of participants by country

Additional file 1: Table S1 lists the chromosome position, location type, allele frequencies, call frequencies and tagged SNPs of the 8 selected ATP12A SNPs. In our analyses we focused on rs10507337 because it is located within a transcription factor binding site and, therefore, it might influence the ATP12A gene expression. The r2 reflecting linkage disequilibrium between rs10507337 and rs12872010 was 0.83. On the other hand, we observed no correlation between rs10507337 and three other ATP12A SNPs, rs963984, rs9553395 and rs2071490, with the r2 values ranged from 0.008 to 0.245, respectively. The genotype frequencies of rs10507337 complied with Hardy-Weinberg equilibrium (P = 0.97). Center specific and cumulative rs10507337 genotype distributions are listed in Table'2. As there were only 5 (0.4%) rs10507337 C allele homozygotes (all from the Belgian center), we contrasted rs10507337 C allele carriers with TT homozygotes and, therefore, use a recessive association model. Additional file 1: Table S3 lists the clinical and echocardiographic characteristics of the participants by rs10507337 genotypes. In unadjusted analysis, rs10507337 C allele carriers had lower heart rate, higher transmitral E, E/A ratio, TDI e' and e'/a' ratio compared to TT homozygotes. There were no other differences in clinical or echocardiographic parameters between these two groups.

Table 2 rs10507337 genotype distribution by country

We used the selected 5 tag SNPs of ATP12A (rs9553395, rs10507337, rs2071490, rs12872010 and rs963984; Additional file 1: Table S1) to reconstruct haplotypes. Haplotype frequencies for TTTCC, TTCCA, TTCCC, CTTCC, CTCCA, CCTTC were 67.4%, 7.1%, 8.5%, 6.9%, 2.3% and 5.3%, respectively.

Composite diastolic score

Because of the high intra-correlation of Doppler diastolic velocities, we summarized these traits into a single diastolic score using principal component analysis. The first principal component accounted for 61.5% of the overall variance in the 7 contributing Doppler indexes. Therefore, we used the first principal component as a normally distributed summary score (mean 0; standard deviation 1; Additional file 1: Figure S2). Figure'1 shows factor loadings for each of the diastolic Doppler measures which we used to calculate the composite diastolic score. Participants with better diastolic function profile (higher transmitral and tissue Doppler velocities during early diastole and lower E/e') had higher values of the composite score. In stepwise multiple regression, the composite diastolic score significantly and independently decreased with age, body mass index, heart rate, diastolic blood pressure and left ventricular mass index (Table'3). The explained variance totaled 79.2%.

Figure 1
figure1

Loading of composite diastolic Doppler score. Absolute values of bivariate correlations (factor loadings) between the composite diastolic score and each of the diastolic parameters. The correlations of transmitral A, TDI a' and E/e' ratio to the composite score were negative and are represented by non-filled points.

Table 3 Correlates of the composite diastolic Doppler score selected in stepwise regression

Population-based association study

We adjusted the model for important covariables such as sex, age, body mass index, heart rate and diastolic blood pressure as well as for non-independence of observations within families and centers. Table'4 lists the adjusted Doppler diastolic indexes and composite score by rs10507337 genotype by country and in all centers combined. When all centers were pooled (Table'4), subjects carrying at least one rs10507337 C allele had significantly higher values of TDI e' (P = 5.8×10-5), transmitral E/A (P = 0.003), TDI e'/a' (P = 0.003) and composite diastolic score (P = 0.0001). The heterogeneity tests (P = 0.04) suggested that there was borderline but significant heterogeneity among the centers for peaks E and e' , e'/a' ratio and composite diastolic score (Table'4). After excluding from the analyses the Polish participants, who were younger than the participants from other centers, we did not observe a heterogeneity among the remaining centers (P > 0.40). Haplotype analysis further confirmed that in all subjects combined, composite diastolic score (P = 0.001; Table'5) and related Doppler parameters such as TDI e', E/A, e'/a' (P ± 0.009; Table'6) were significantly higher in CCTTC carriers than in noncarriers. This is the only haplotype that included rs10507337 C allele.

Table 4 Adjusted Doppler diastolic indexes and composite score by rs10507337 by country and all centers combined
Table 5 Composite diastolic Doppler score by ATP12A haplotypes in all centers combined
Table 6 Adjusted Doppler diastolic indexes and composite score by CCTTC haplotype in all centers combined

Our findings remained consistent after exclusion of subjects with history of coronary heart disease or valve abnormalities (n = 67; Additional file 1: Table S4) or subjects on antihypertensive drugs (n = 315; Additional file 1: Table S5). For other LV phenotypes (Additional file 1: Table S6), none of the phenotype-genotype associations reached statistical significance. With adjustment for multiple testing applied, we observed that TDI e', e'/a' ratio and composite diastolic score were significantly (P ± 0.008) higher in carriers of rs12872010 T allele compared to non-carriers (Additional file 1: Table S7). These associations were somewhat expected given the high correlation between ATP12A rs10507337 and rs12872010. We did not observe any significant associations with other tested SNPs (Additional file 1: Tables S8, S9 and S10) or haplotypes (Table'5).

Family-based association study

Our study population (n = 1166) included 749 founders and 417 offspring from 357 families. In our family-based analyses we used 74 informative offspring belonging to 35 families (mean age 34.6 ± 8.2 years; 56.9% women). The number of offspring per informative family amounted to one in 10 families, two in 17 families, three in 4 families, and more than three in 4 families. We adjusted the QTDT analyses as described in the previous section. For LV diastolic function indexes in relation to rs10507337 genotypes, the orthogonal model did not demonstrate population stratification (P ± 0.08). Transmitral E/A, TDI e' , TDI e'/a' and composite diastolic score significantly increased with transmission of at least one C allele to offspring (Table'7). The effect sizes of the within-family components averaged 0.17 (χ2 = 6.61, P = 0.010) for E/A, 1.19 cm/s (χ2 = 8.61, P = 0.003) for e', 0.26 (χ2 = 9.13, P = 0.003) for e'/a' and 0.62 (χ2 = 9.04, P = 0.003) for composite score (Table'7).

Table 7 Results of the QTDT analyses for the rs10507337 C allele in relation to Doppler diastolic indexes and composite score

Discussion

The main finding of the present study was that LV diastolic function as assessed by conventional and tissue Doppler indexes was associated with genetic variation in the ATP12A promoter. Carriers of rs10507337 C allele had better pattern of myocardial relaxation as compared to non-carriers. Because rs10507337 C allele occurred only in one reconstructed haplotype (CCTTC), we found also the significant associations between this haplotype and LV diastolic function indexes. The family-based analyses included only 74 informative offspring but nevertheless confirmed that transmission of at least one C allele to offspring was associated with higher transmitral and tissue Doppler velocities during early diastole and composite diastolic score. On the contrary, we did not observe any association of the genetic variants in ATP12A and other LV phenotypes including LV mass index and ejection fraction.

The gold standard for assessing diastolic function remains the pressure-volume relationship, but this requires an invasive approach. Doppler measurements of mitral inflow and the TDI technique open up the possibility of evaluating non-invasively diastolic function [6]. In our study, we assessed LV diastolic function non-invasively using the transmitral flow and the TDI mitral annular velocities. Previous studies validated these indexes versus invasive measures of diastolic function [6]. LV diastolic dysfunction is defined as functional abnormalities that exist during LV relaxation and filling. Impaired myocardial relaxation is characterized by decreased early (E peak), but enhanced atrial LV filling (A peak) as well as less vigorous mitral annulus motion during early diastole (TDI e' peak). Moreover, e' peak velocity along the LV longitudinal axis is less susceptible to the effects of an increased preload and therefore provides a more direct measure of myocardial relaxation than, for instance, the transmitral E peak velocity. In addition, combining transmitral flow velocity with annular velocity (E/e' ratio) might be a tool for assessing the LV filling pressure, which combines the influence of the transmitral driving pressure and myocardial relaxation [6].

To our knowledge, no previous study reported the association of cardiovascular phenotypes with the ATP12A gene. So far, only a whole DNA array survey in spontaneously hypertensive rats identified the ATP12A gene among probable candidate genes for hypertension [21]. Kinoshita et al.[21] demonstrated that lower ATP12A expression might contribute to development of hypertension in spontaneously hypertensive rats. There is a high level of structural similarity (up to 86%) in the ATP12A gene between humans and rodents [16]. The ATP12A gene product is a catalytic alpha (ATP12a') subunit of an ATP dependent transmembrane electrolyte pump family (X+-K+-ATPases) [16],[19] which is expressed in human myocardium (http://www.genecards.org/cgi-bin/carddisp.pl?gene=ATP12A). Functional X+-K+-ATPases are heterodimers, composed of alpha and beta subunits. Alpha subunits of ATPases have catalytic and ion-binding properties [22]. Previous studies in Xenopus oocytes and human HEK293 cells demonstrated that ATP12a' could be assembled to ATPase that functions as a Na+-K+-pump [17],[23]. Moreover, ATP12a' is sensitive to cardiac glycoside ouabain, specific inhibitor of Na+-K+-ATPases [24].

Small perturbations in Na+ concentration due to changes of Na+-K+-ATPase expression and activity, might lead to important changes of Ca2+ extrusion, a process central for myocardial contractility and relaxation [11],[12]. Schwinger et al.[25],[26] demonstrated that failing human myocardium had low level of Na+-K+-ATPase concentration and, a consequence, had elevated intracellular Na+ levels [12],[13]. Moreover, studies in rodent hearts showed that decrease in ATP production effects Na+-K+-ATPase activity and leads to elevated intracellular Na+ concentration [27].

Involvement of Na+-K+-ATPase in electrolyte homeostasis makes ATP12A a possible candidate gene that might modify cardiac function. In our study, we investigated the genetic variation which is located in a transcription factor binding site of the ATP12A promoter area. Several transcription factors can bind this promoter region including CTCF (CCCTC-binding factor) and NF-kB that both act in NF-kB cascade [19]. Of notice, levels of NF-kB were increased in failing human myocardium [28]. Therefore, in conditions related to disturbances in ATP production such as ischemia or heart failure, these transcription factors might modify expression of Na+-K+-ATPase.

In our study we observed that subjects with minor C allele in the ATP12A promoter area had significantly higher transmitral and tissue Doppler velocities during early diastole and, therefore, demonstrated an enhanced early myocardial relaxation compared to overall population mean. Overall, these findings might be indicative of the functional importance of the described genetic variation in ATP12A. However, we do not know whether rs10507337 C allele is associated with increased or decreased expression of the ATP12A gene. The precise mechanisms underlying the association of LV diastolic function with ATP12A remain to be elucidated.

Our echocardiographic estimation of diastolic function included measurement of diastolic blood flow (A, E) and mitral annular movement (a', e'), and calculation of several ratios (E/A, e'/a', E/e'). Since each of these Doppler indexes represents a somewhat distinct feature of diastolic function, none of them is sufficient to stand alone [1]. On the other hand, because diastolic Doppler indexes are measuring the same construct, e.g. velocities of transmitral blood flow and LV mitral annulus movement during diastole and its ratios, it is possible to reduce the measured phenotypes into a principal component (composite score) that will account for most of the variance in the variables. Indeed, all used diastolic indexes are intra-correlated and, therefore, some redundancy is expected. This continuously and normally distributed diastolic Doppler score accounted for 61.5% of the variance in the contributing indexes. In our study, subjects with better diastolic function had higher values of the composite score. Similarly, the composite diastolic score was significantly higher in rs10507337 C allele carriers and in CCTTC haplotype carriers. Because of the linkage disequilibrium between rs10507337 and rs12872010, this association was also observed in rs12872010 minor allele carriers. Hence, our findings remained consistent and we may conclude that the composite score of diastolic Doppler velocities might be used in genomic studies as a phenotype describing LV diastolic function.

The present study must be interpreted within the context of its limitations and strengths. All participants were white Europeans. Thus, the association cannot be generalized to other ethnic or racial groups. Transmitral and TDI diastolic velocities are quantitative echocardiographic traits that arise from complex interaction between multiple genes, hemodynamic, and environmental factors and are prone to measurement error. For the present study, one experienced observer in each center performed all echocardiographic examinations. Digitally stored images were post-processed centrally by one observer with high reproducibility [4]. There was also a high degree of internal consistency between the results of the population-based and family-based analyses. The between-family components for rs10507337 were not statistically significant, which makes it unlikely that our results are driven by population stratification.

Conclusion

Left ventricular diastolic function as assessed by conventional and tissue Doppler indexes was associated with genetic variation in the ATP12A promoter. We observed that carriers of rs10507337 C allele or CCTTC haplotype had better pattern of myocardial relaxation as compared to non-carriers. Further studies are necessary to clarify the functional significance of this genetic variation.

Additional file

Abbreviations

A:

the transmitral peak late diastolic velocity

a':

the peak late diastolic mitral annular velocity

ATP:

Adenosine triphosphate

ATP12A :

Gene that encode the catalytic alpha subunit of Na+/K+-ATPase

BSA:

Body surface area

E:

the transmitral peak early diastolic velocity

E/A ratio:

E velocity/A velocity ratio

e':

the peak early diastolic mitral annular velocity

e'/a' ratio:

e' velocity/a' velocity ratio

E/e' ratio:

E velocity/e' velocity ratio

EPOGH:

European Project on Genes in Hypertension

EF:

Ejection fraction

FLEMENGHO:

Flemish Study on Environment, Genes and Health Outcomes

LV:

Left ventricle

MAF:

Minor allele frequency

QTDT:

Quantitative transmission disequilibrium test

SD:

Standard deviation

TDI:

Tissue Doppler Imaging

References

  1. 1.

    Kuznetsova T, Herbots L, Yin Y, Stolarz-Skrzypek K, Staessen J: Systolic and diastolic left ventricular dysfunction: from risk factors to overt heart failure. Expert Rev Cardiovasc Ther. 2010, 8: 251-258. 10.1586/erc.10.3.

    Article  PubMed  Google Scholar 

  2. 2.

    Abhayaratna W, Marwick T, Smith W, Becker N: Characteristics of left ventricular diastolic dysfunction in the community: an echocardiographic survey. Heart. 2006, 92: 1259-1264. 10.1136/hrt.2005.080150.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Kloch-Badelek M, Kuznetsova T, Sakiewicz W, Tikhonoff V, Ryabikov A, González A, López B, Thijs L, Jin Y, Malyutina S, Stolarz-Skrzypek K, Casiglia E, Díez J, Narkiewicz K, Kawecka-Jaszcz K, Staessen JA: Prevalence of left ventricular diastolic dysfunction in European populations based on cross-validated diagnostic thresholds. Cardiovasc Ultrasound. 2012, 10: 10-10.1186/1476-7120-10-10.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Kuznetsova T, Herbots L, López B, Jin Y, Richart T, Thijs L, González A, Herregods MC, Fagard RH, Díez J, Staessen JA: Prevalence of left ventricular diastolic dysfunction in a general population. Circ Heart Fail. 2009, 2: 105-112. 10.1161/CIRCHEARTFAILURE.108.822627.

    Article  PubMed  Google Scholar 

  5. 5.

    Redfield M, Jacobsen S, Burnett JJ, Mahoney D, Bailey K, Rodeheffer R: Burden of systolic and diastolic ventricular dysfunction in the community: appreciating the scope of the heart failure epidemic. JAMA. 2003, 289: 194-202. 10.1001/jama.289.2.194.

    Article  PubMed  Google Scholar 

  6. 6.

    Nagueh SF, Appleton CP, Gillebert TC, Marino PN, Oh JK, Smiseth OA, Waggoner AD, Flachskampf FA, Pellikka PA, Evangelista A: Recommendations for the evaluation of left ventricular diastolic function by echocardiography. J Am Soc Echocardiogr. 2009, 22: 107-133. 10.1016/j.echo.2008.11.023.

    Article  PubMed  Google Scholar 

  7. 7.

    Wang M, Yip GW, Wang AY, Zhang Y, Ho PY, Tse MK, Yu CM, Sanderson JE: Tissue Doppler imaging provides incremental prognostic value in patients with systemic hypertension and left ventricular hypertrophy. J Hypertens. 2005, 23: 183-191. 10.1097/00004872-200501000-00029.

    Article  PubMed  Google Scholar 

  8. 8.

    Dokainish H, Zoghbi WA, Lakkis NM, Ambriz E, Patel R, Quinones MA, Nagueh SF: Incremental predictive power of B-type natriuretic peptide and tissue Doppler echocardiography in the prognosis of patients with congestive heart failure. J Am Coll Cardiol. 2005, 45: 1223-1226. 10.1016/j.jacc.2005.01.025.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Kuznetsova T, Thijs L, Knez J, Herbots L, Zhang Z, Staessen J: Prognostic value of left ventricular diastolic dysfunction in a general population. J Am Heart Assoc. 2014, 3 (3): e000789-10.1161/JAHA.114.000789.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Kloch-Badelek M, Knez J, Tikhonoff V, Thijs L, Sakiewicz W, Ryabikov A, Stolarz-Skrzypek K, Jin Y, Malyutina S, Casiglia E, Narkiewicz K, Czarnecka D, Kawecka-Jaszcz K, Staessen JA, Kuznetsova T: Heritability and other determinants of left ventricular diastolic function in the family-based population study. J Hypertens. 2014, 32: 1854-1861. 10.1097/HJH.0000000000000256.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Despa S, Bers D: Na+ transport in the normal and failing heart - remember the balance. J Mol Cell Cardiol. 2013, 61: 2-10. 10.1016/j.yjmcc.2013.04.011.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Pieske B, Houser S: [Na+]i handling in the failing human heart. Cardiovasc Res. 2003, 57: 874-886. 10.1016/S0008-6363(02)00841-6.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Pogwizd S, Sipido K, Verdonck F, Bers D: Intracellular Na in animal models of hypertrophy and heart failure: contractile function and arrhythmogenesis. Cardiovasc Res. 2003, 57: 887-896. 10.1016/S0008-6363(02)00735-6.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Pieske B, Maier L, Piacentino V, Weisser J, Hasenfuss G, Houser S: Rate dependence of [Na+] i and contractility in nonfailing and failing human myocardium. Circulation. 2002, 106: 447-453. 10.1161/01.CIR.0000023042.50192.F4.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Blanco G, Mercer R: Isozymes of the Na-K-ATPase: heterogeneity in structure, diversity in function. Am J Physiol. 1998, 275: 633-650.

    Google Scholar 

  16. 16.

    Jaisser F, Beggah A: The nongastric H + -K + -ATPases: molecular and functional properties. Am J Physiol. 1999, 276: 812-824.

    Google Scholar 

  17. 17.

    Grishin A, Caplan M: ATP1AL1, a member of the non-gastric H, K-ATPase family, functions as a sodium pump. J Biol Chem. 1998, 273: 27772-27778. 10.1074/jbc.273.43.27772.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Salvi E, Kutalik Z, Glorioso N, Benaglio P, Frau F, Kuznetsova T, Arima H, Hoggart C, Tichet J, Nikitin YP, Conti C, Seidlerova J, Tikhonoff V, Stolarz-Skrzypek K, Johnson T, Devos N, Zagato L, Guarrera S, Zaninello R, Calabria A, Stancanelli B, Troffa C, Thijs L, Rizzi F, Simonova G, Lupoli S, Argiolas G, Braga D, D'Alessio MC, Ortu MF, et al: Genomewide association study using a high-density single nucleotide polymorphism array and case'control design identifies a novel essential hypertension susceptibility locus in the promoter region of endothelial NO synthase. Hypertension. 2012, 59: 248-255. 10.1161/HYPERTENSIONAHA.111.181990.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Sverdlov V, Kostina M, Modyanov N: Genomic organization of the human ATP1AL1 gene encoding a ouabain-sensitive H, K-ATPase. Genomics. 1996, 32: 317-327. 10.1006/geno.1996.0125.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Abecasis G, Cardon L, Cookson W: A general test of association for quantitative traits in nuclear families. Am J Hum Genet. 2000, 66: 279-292. 10.1086/302698.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Kinoshita K, Ashenagar M, Tabuchi M, Higashino H: Whole rat DNA array survey for candidate genes related to hypertension in kidneys from three spontaneously hypertensive rat substrains at two stages of age and with hypotensive induction caused by hydralazine hydrochloride. Exp Ther Med. 2011, 2: 201-212.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Geering K: The functional role of beta subunits in oligomeric P-type ATPases. J Bioenerg Biomembr. 2001, 33: 425-438. 10.1023/A:1010623724749.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Crambert G, Horisberger J, Modyanov N, Geering K: Human nongastric H+−K + −ATPase: transport properties of ATP1al1 assembled with different beta-subunits. Am J Physiol Cell Physiol. 2002, 283: 305-314. 10.1152/ajpcell.00590.2001.

    Article  Google Scholar 

  24. 24.

    Modyanov N, Mathews P, Grishin A, Beguin P, Beggah AT, Rossier BC, Horisberger JD, Geering K: Human ATP1AL1 gene encodes a ouabain-sensitive H-K-ATPase. Am J Physiol. 1995, 269: 992-997.

    Google Scholar 

  25. 25.

    Schwinger R, Wang J, Frank K, Müller-Ehmsen J, Brixius K, McDonough AA, Erdmann E: Reduced sodium pump alpha1, alpha3, and beta1-isoform protein levels and Na+, K + −ATPase activity but unchanged Na'+ −Ca2++ exchanger protein levels in human heart failure. Circulation. 1999, 99: 2105-2112. 10.1161/01.CIR.99.16.2105.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Schwinger R, Bundgaard H, Muller-Ehmsen J, Kjeldsen K: The Na, K-ATPase in the failing human heart. Cardiovasc Res. 2003, 57: 913-920. 10.1016/S0008-6363(02)00767-8.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Jansen M, Shen H, Zhang L, Wolkowicz P, Balschi J: Energy requirements for the Na'+ gradient in the oxygenated isolated heart: effect of changing the free energy of ATP hydrolysis. Am J Physiol Heart Circ Physiol. 2003, 285: 2437-2445.

    Article  Google Scholar 

  28. 28.

    Wong S, Fukuchi M, Melnyk P, Rodger I, Giaid A: Induction of cyclooxygenase-2 and activation of nuclear factor-kB in myocardium of patients with congestive heart *failure. Circulation. 1998, 98: 100-103. 10.1161/01.CIR.98.2.100.

    CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgments

The European Union (grants IC15-CT98-0329-EPOGH, LSHM-CT-2006-037093-InGenious HyperCare, HEALTH-F4-2007-201550-HyperGenes, HEALTH-2011-278249-EU-MASCARA, HEALTH-F7-305507-HoMAGE, and ERC Advanced Grant-2011-294713-EPLORE) supported the Research Unit Hypertension and Cardiovascular Epidemiology (Leuven, Belgium). The Research Unit Hypertension and Cardiovascular Epidemiology also received grants from the Fonds voor Wetenschappelijk Onderzoek Vlaanderen, Ministry of the Flemish Community, Brussels, Belgium (grants G.0734.09, G.0880.13 and G. 0881.13). The Department of Health, University of Milano and Genomics and Bioinformatics Platform, Fondazione Filarete are supported by InterOmics (PB05 MIUR-CNR Italian Flagship Project).

Author information

Affiliations

Authors

Corresponding author

Correspondence to Tatiana Kuznetsova.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TK had full access to all of the data and takes responsibility for the integrity of the data and the accuracy of the statistical analysis. Study concept and design: JK, ES, JAS, TK. Acquisition of the phenotypic data: JK, VT, KSS, AR, MKK, SM, EC, DC, KKJ, TK. Molecular analysis: ES, DB, DC, TN. Statistical analysis and manuscript preparation: JK, ES, LT, TK. Critical revision of the manuscript for important intellectual content: VT, KSS, AR, DB, MKK, SM, EC, DC, KKJ, DC, TN, JAS. Obtained funding: SM, EC, DC, KKJ, DC, JAS, TK. All authors read and approved the final manuscript.

Electronic supplementary material

a: Supplemental methods.

Additional file 1: b: Figure S1. Local plot of ATP12A region mapped by selected SNPs. c: Figure S2. Distribution and probability plot for the first principal component Describing LV diastolic function. d: Table S1. Selected SNPs for ATP12A. e: Table S2. Characteristics of participants by sex. f: Table S3. Characteristics of participants by rs10507337 genotypes. g: Table S4. Adjusted Doppler diastolic indexes and composite score by ATP12A rs10507337 in subjects without a previous history of coronary heart disease or valve abnormalities. h: Table S5. Adjusted Doppler diastolic indexes and composite score by ATP12A rs10507337 in untreated subjects. i: Table S6. Adjusted LV echocardiographic phenotypes by rs10507337. j: Table S7. Adjusted Doppler diastolic indexes and composite score by rs12872010 by country and in all centres combined. k: Table S8. Adjusted Doppler diastolic indexes and composite score by rs9553395 by country and in all centres combined. l: Table S9. Adjusted Doppler diastolic indexes and composite score by rs2071490 by country and in all centres combined. m: Table S10. Adjusted Doppler diastolic indexes and composite score by rs963984 by country and in all centres combined. (DOCX 132 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Knez, J., Salvi, E., Tikhonoff, V. et al. Left ventricular diastolic function associated with common genetic variation in ATP12Ain a general population. BMC Med Genet 15, 121 (2014). https://doi.org/10.1186/s12881-014-0121-6

Download citation

Keywords

  • Epidemiology
  • Echocardiography
  • Diastolic function
  • ATP12A