Skip to main content

Contribution of sarcomere gene mutations to left atrial function in patients with hypertrophic cardiomyopathy

Abstract

Background

Left atrial (LA) enlargement and dysfunction are related to clinical course in patients with hypertrophic cardiomyopathy (HCM). We aimed to investigate genetic contribution to LA structural and functional remodeling.

Methods

Two hundred twelve patients were consecutively enrolled, and echocardiography and extensive genetic analysis were performed. Cardiac magnetic resonance (CMR) was performed in 135 patients. Echocardiography was also performed in controls (n = 30).

Results

Patients with HCM had lower late-diastolic mitral annular velocity (a’) and higher LA volume index (LAVI) than controls. Patients with pathogenic or likely pathogenic sarcomere gene mutations (PSM, n = 67, 32%) had higher LAVI and lower CMR-derived LA total emptying fraction (37.0 ± 18.5 vs. 44.2 ± 12.4%, p = 0.025). In patients without AF (n = 187), the PSM had lower a’ (6.9 ± 2.0 vs. 7.8 ± 1.9 cm/s, p = 0.004) than others. The PSM had higher prevalence and amount of late gadolinium enhancement (LGE) in the left ventricle (LV). In multivariate analysis, PSM was significantly related to lower a’ independent of E/e’, LV mass index, and LAVI. However, the relation significantly attenuated after adjustment for the extent of LGE in the LV, suggesting common myopathy in the LV and LA. In addition, PSM was significantly related to lower LA total emptying fraction independent of age, E/e’, s’, LV ejection fraction, LV myocardial global longitudinal strain and %LGE mass.

Conclusions

PSM was related to LA dysfunction independent of LV filling pressure and LAVI, suggesting its contribution to atrial myopathy in HCM.

Peer Review reports

Background

Left atrial (LA) enlargement and dysfunction are related to the development of atrial fibrillation (AF) and exercise intolerance in patients with hypertrophic cardiomyopathy (HCM) [1,2,3]. LA dysfunction is usually related to chronic structural remodeling owing to long-standing pressure and volume overload. In addition, LA structural remodeling such as enlargement or fibrosis is a marker of development of AF in HCM [4]. Recent studies also showed that LA function measured by strain imaging is an independent predictor of poor outcome in patients with HCM [1, 2]. However, it has been reported that in HCM, LA dysfunction can occur as a myopathic process itself, irrespective of loading conditions [4, 5]. In this regard, genetic contribution to LA structural and functional remodeling has not been well investigated. Tissue Doppler imaging-based late mitral annular velocity (a’) and LA phasic function have been validated as good LA functional parameters reflecting active contraction in several studies and are widely used [6]. Herein, we sought to evaluate the relationship between related gene mutations and LA function in patients with HCM. In addition, we aimed to study the direct effect on LA function as a process of LA myopathy, independent of LA loading conditions.

Materials and methods

Study population

A total of 432 patients treated at a single center were enrolled in an HCM Registry from 2006 to 2014. Among them, 220 patients were excluded owing to insufficient data, follow-up loss, or declining study enrollment. Finally, 212 patients underwent genetic testing. The patients had maximal left ventricular (LV) hypertrophy greater than 13 mm and a ratio of maximal thickness to posterior wall thickness greater than 1.3 without an underlying cause of hypertrophy, such as uncontrolled hypertension or aortic stenosis. Patterns of LV hypertrophy were classified as ApHCM and non-ApHCM (asymmetrical hypertrophy, diffuse hypertrophy, and focal segmental hypertrophy). All patients underwent screening for Fabry disease and were confirmed negative for the galactosidase alpha variant. For comparison, conventional echocardiography was performed in controls. The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki, and it was approved by institutional review board of Gangnam Severance Hospital (3–2015-0019). Written informed consent was obtained from all participants.

Genetic testing and analysis

HCM gene panel (nDNA) and mtDNA design

HCM genes consisted of 8 validated sarcomere genes and 25 putative HCM genes [7]. A comprehensive HCM-specific panel, consisting of 82 nuclear DNAs (nDNAs: 33 sarcomere-associated genes, 5 phenocopy genes, and 44 nuclear genes linked to mitochondrial cardiomyopathy) and 37 mitochondrial DNAs (mtDNAs), was analyzed (Supplementary Table S1).

DNA preparation, library construction and sequencing of HCM gene panel and mtDNA

The details are described in Supplementary Method 1 to 3.

Identification of potential pathogenic mtDNA variants

Mitochondrial genome databases, including MITOMAP [8] and Human Mitochondrial Genome Database (mtDB) [9] and Phylotree [10] were referred to validate the detected variants. Novel and rare non-haplogroup-associated variants were further evaluated for their potential pathogenicity based on the variant’s location, changes in the amino acid sequence, and evolutionary conservation [11]. We have assessed potential pathogenicity using multiple software programs including Polyphen2, Fathmmw, Mutation Assessor, and PROVEAN. When the majority of computational evidence supported a deleterious effect, we have assigned novel and rare non-haplogroup-associated variants as damaging mtDNA variants [12]. Data analysis of mitochondrial genome are described in Supplementary Method 4.

Data analysis of the HCM gene panel

The Burrows–Wheeler aligner algorithm with default parameters was used to align reads to the human reference genome sequence GRCh37 [13]. SAMTools was used to convert the sequence alignment map file to the BAM format [14]. Duplicates were sorted and removed using the Picard tool (http://broadinstitute.github.io/picard/). The Genome Analysis Toolkit was used for indel realigning and base-quality score re-calibration [15]. Variants were further filtered with altered allele frequency > 30%, 50× coverage. For 33 HCM genes, annotated variants using ANNOVAR [16] and Variant Effect Predictor (http://asia.ensembl.org/info/docs/tools/vep/index.html) were classified as pathogenic and likely pathogenic based on refined American College of Medical Genetics and Genomics (ACMG) standards and guidelines for inherited cardiac conditions [12, 17,18,19]. For 44 mitochondria-related nuclear DNA genes (recessive conditions), annotated variants were classified as pathogenic and likely pathogenic based on ACMG guidelines [17]. And we adapted gnomAD AF cutoff 0.01% as the moderate level of evidence supporting pathogenicity (ACMG/AMP criterion PM2) based on maximum credible population AF (http://cardiodb.org/allelefrequencyapp/) [17, 18].

Echocardiographic analysis

Comprehensive echo-Doppler evaluation was performed according to the current American Society of Echocardiography guidelines [20]. LV ejection fraction was measured by biplane Simpson’s method. LA volume was measured at the end-systole by the ellipsoidal method, and LA volume index was calculated as LA volume/body surface area (BSA). Peak early (E) and late (A) diastolic mitral inflow velocities were measured in apical four-chamber view. The filter was set to exclude high frequency signal, and the Nyquist limit was adjusted to a range of 15 to 20 cm/s. Gain and sample volume were minimized to allow for a clear tissue signal with minimal background noise. Systolic (s’), early (e’) and late (a’) diastolic velocities of the mitral annulus were measured from the apical 4-chamber view with a 2- to 5-mm sample volume placed at the septal corner of the mitral annulus. (Fig. 1) The ratio of E/e′ was calculated. LV wall thickness was measured in all cross-sectional planes. Continuous wave Doppler was used to measure the peak velocity across the LV outflow tract (LVOT), and the pressure gradient was calculated using the Bernoulli equation, as follows: 4 × (peak velocity across the LVOT)2 [21]. It was measured at rest and during Valsalva maneuver. LVOT obstruction was defined as a systolic pressure gradient ≥30 mmHg.

Fig. 1
figure 1

Tissue Doppler index (right panel) and phasic left atrial function measurement with cardiac magnetic resonance imaging (right panel) in a patient with pathogenic sarcomere gene mutation (TNNI3 mutation, a, b) and without sarcomere mutations (c, d)

Cardiac magnetic resonance imaging (CMR) and LV chamber performance assessment

CMR was performed using a 1.5-T scanner (Magnetom Avanto®; Siemens Medical Solutions, Erlangen, Germany) with a phased array body coil. (Supplementary Method 5).

Extent of LGE assessment

From the LGE images, the LV was divided into 16 segments [22]. Presence of LGE involvement in each segment and the total number of LGE-involving segments were measured. In addition, the percentage of LGE in LV mass were measured using dedicated quantitative analysis software (Qmass®MR 8.1, Medis, Leiden, Netherland) using LGE images with PSIR sequence. To improve the reproducibility, experienced radiologist and cardiologist with more than 10 years of experience analyzed the LGE sizes. In each short-axis slice image, the boundaries of contrast-enhanced areas were automatically traced. On LGE-MR images, myocardium with abnormal enhancement was defined as an area of hyper-enhancement more than five standard deviations from the remote myocardium. Remote myocardium was defined as non-enhanced myocardium, opposite of the hyper-enhanced myocardium [23]. The maximum signal was determined by computer-assisted window thresholding of the enhanced area. Obvious artifacts, such as those caused by motion, were excluded using a tool from the software package. Total LGE amount was calculated by summation of all slice volumes of enhancement.

Myocardial strain analysis using feature tracking CMR

Myocardial strain analysis using feature tracking CMR was performed in 135 patients with semi-automated software (Qstrain®MR 2.0, Medis, Leiden, Netherland). LV endocardial borders were manually drawn at a reference frame. LV endocardial and epicardial borders were manually traced in 2-, 3-, 4-chamber long-axis views at end-systolic and end-diastolic phase. LV global longitudinal strain (GLS) was obtained from averaging longitudinal strains of apical 2-, 3-, and 4-chamber view. LA endocardial border was manually traced in 4-chamber long-axis view using LV end-diastole as reference phase. LA global longitudinal strain is defined as the average peak strain value. LA maximal, pre-contraction (pre-A in cases without AF) and minimal volume were measured [24]. LA total emptying fraction was calculated as (LA maximal volume – LA minimal volume) / LA maximal volume. LA reservoir fraction, as (LA maximal volume – LA minimal volume / LA minimal volume), LA conduit fraction as (LA maximal volume – LA pre-A volume) / LA maximal volume, and LA active emptying fraction as (LA pre-A volume – LA minimal volume) / LA pre-A volume [6]. (Fig. 1).

Statistical analysis

Continuous variables with normal distributions are reported as the mean ± standard deviation or 95% confidence interval. Student’s t-tests were used to compare the means of continuous variables that were approximately normally distributed between the two groups. Normality was determined using the Shapiro–Wilk test. Categorical variables are reported as counts (or percentages) and were compared using chi-square tests. For comparison of more than two groups, analysis of variance was performed with post-hoc analysis (LSD) for subgroup comparison. For the correlates of LA function, Pearson’s correlation coefficients were determined and Pearson’s correlation analyses were performed. For the multivariable analysis, linear regression analysis or logistic regression analysis was performed with variables with p < 0.05 in univariate analysis to check the independence of the variables. All statistical analyses were performed using SPSS version 25.0 (IBM Corp., Armonk, NY, USA). A two-sided p-value < 0.05 was considered statistically significant.

Results

Baseline characteristics

The mean age of the patients was 59 ± 14 years, and 63 (30%) of them were female. Of the total patients, 49 (23%) had obstructive HCM; 100 patients (47%) had ApHCM, and 64 (64%, 64/100) of these had pure-type ApHCM. Patients with HCM had higher LA volume index and lower a’ than those of age- and sex- matched controls. CMR was performed in 135 patients (LGE in 133 patients). The mean LV mass index was 85.7 ± 22.9 g/m2; 92 (69%) patients showed evidence of LGE. The average number of LGE segments and %LGE mass were 3.5 ± 3.4 and 8.0 ± 9.7%, respectively (Table 1).

Table 1 Clinical, echocardiographic, electrocardiographic, and cardiac magnetic resonance imaging findings in controls and patients with hypertrophic cardiomyopathy according to sarcomere gene mutations

Genetic characteristics

Based on the American College of Medical Genetics and Genomics guidelines, [19] 67 of 212 (32%) patients had 71 pathogenic or likely pathogenic variants in 33 sarcomere-associated genes (33 MYBPC3, 19 MYH7, 14 TNNI3, 2 MYH6, 1 JPH2, 1 TNNC1, and 1 MYL3). Four patients harbored more than one variants in HCM genes. Homozygous or compound heterozygous variants in MYBPC3 were identified in one patient and co-variants were identified in three patients (two had MYBPC3 and MYH7, and one had MYBPC3 and JPH2). In total, 26 (12%) patients had probably damaging mtDNA variants (one patient was missed due to non-analysis of mtDNA), 15 (7%) had mitochondria-related nDNA variants, and 1 had a pathogenic variant in GAA. Seven patients (3%) had both pathogenic sarcomere gene variants and mitochondria-related variants. Detailed genetic variants of detected sarcomere genes, mitochondria-related nDNAs and damaging mtDNAs are shown in Supplementary Table S2, S3 and S4 [12]. The patients with pathogenic or likely pathogenic sarcomere gene mutations (PSM) had higher prevalence of LGE (90% vs. 63% vs. 50%, p = 0.001) and number of LGE-involved segments (4.9 ± 2.8 vs. 2.9 ± 3.3 vs. 2.6 ± 3.8, p = 0.005) than patients without any mutations and those with mitochondria-related mutations.

Effects of sarcomere mutations on LA function

The PSM (n = 67, 32%) had higher LA volume index (41.0 ± 22.0 vs. 34.9 ± 15.0 mL/m2, p = 0.045) and tendency of higher prevalence of AF (16.4% vs. 9.7%, p = 0.156). Even in patients without AF (n = 187), PSM had significantly lower A (61.9 ± 19.6 vs. 72.4 ± 20.8 cm/s, p = 0.001) and a’ (6.9 ± 2.0 vs. 7.8 ± 1.9 cm/s, p = 0.004). PSM had lower LA total emptying fraction by CMR than others (44.2 ± 12.4 vs. 37.0 ± 18.5%, p = 0.025). However, patients with only mitochondria-related mutations had higher a’, compared with PSM (8.6 ± 3.3 vs. 6.9 ± 2.0 cm/s, p < 0.001) (Fig. 2). Within PSM, those with TNNI3 mutation and MYH6 had lower a’, A and higher LA volume index than others (Fig. 3). Patients with thick filament mutations (n = 46) had a lower A (64.0 ± 19.0 vs. 71.0 ± 21.2 cm/s, p = 0.048) and a tendency to lower a’ (7.13 ± 1.96 vs. 7.68 ± 1.96 cm/s, p = 0.099) compared with those without thick filament mutations (n = 141).

Fig. 2
figure 2

Comparison of prevalence of atrial fibrillation (a), left atrial volume index (b), and a’ velocity (c) and left atrial total emptying fraction (d) between sarcomere gene mutation-positive group and mutation-negative group. Comparison of left atrial volume index (e), and a’ (f) velocity between sarcomere- and mitochondria-related gene mutation groups*. HCM, hypertrophic cardiomyopathy; AF, atrial fibrillation, LA, left atrial, TDI, tissue Doppler imaging; a’, peak late-diastolic septal mitral annular; CMR, cardiac magnetic resonance imaging; bar indicates standard error; *one patient was missed due to non-analysis of mtDNA

Fig. 3
figure 3

Comparisons of a’ (a), A (b) velocity and left atrial volume index (c) within patients with individual sarcomere gene mutations

Correlates and independent relationship of sarcomere mutations to LA function

Patients with ApHCM had higher a’ than that of patients without ApHCM. In univariate analysis, a’ was significantly correlated to LA volume index, E/e’, e’, and s’. In addition, LV mass index, LA minimal volume, LA total emptying fraction, LA active emptying fraction, 5-year risk of sudden cardiac death, and the presence and number of LGE segments were significantly correlated to a’ but not with LV global longitudinal strain. In multivariate analysis, PSM was significantly related to a’ independent of the presence of ApHCM, E/e’, s’, and LA volume index. When using the definition of a’ < 6 cm/s (less than 3rd quartile value) as LA dysfunction, PSM was also significantly correlated to LA dysfunction independent of ApHCM, E/e’, s’ and LA volume index by echocardiography (Table 2). However, the relation significantly attenuated (p > 0.05) after adjustment for the presence of LGE or number of LGE segment in LV, which suggests common myopathy in LV and LA. In subgroup who underwent CMR, LA total emptying fraction was significantly correlated to age, E/e’, s’, a’, %LGE mass, LV ejection fraction by CMR, LV myocardial GLS and presence of sarcomere mutations. In multivariate analysis, presence of sarcomere mutations was significantly correlated to LA total emptying fraction independent of age, E/e’, s’, LV ejection fraction by CMR, LV myocardial GLS and %LGE mass (Table 3).

Table 2 Univariate and multivariate analysis for a’ velocity
Table 3 Univariate and multivariate analysis for left atrial total emptying fraction

Discussion

This study shows the following major findings. First, LA function was significantly reduced in patients with HCM, compared with controls despite same age and sex distribution. Second, among patients with HCM, the PSM had significantly reduced LA function (measured by a’ in echo and LA total emptying fraction by CMR), compared with patients without sarcomere gene mutations. Third, worse LA function in the PSM was independent of the LA volume index, LV mass index, and E/e’. Fourth, LA function, measured by a’, was significantly correlated to the extent of LGE in LV, and the PSM had significantly higher extent of LGE. Although, the relationship between sarcomere gene mutation and a’ was significantly attenuated after adjustment for the extent of LGE in LV, the LA total emptying fraction was significantly correlated to sarcomere gene mutation independent of %LGE mass, which suggests that LA dysfunction is a unique finding of HCM as an LA myopathy related to sarcomere gene mutations independent of LA loading conditions. Although we did not measure LA fibrosis by LGE quantification, this finding suggests that sarcomere gene mutation may contribute to both LV and LA myopathy, and may then induce LA dysfunction. This sarcomere gene mutation- related LA dysfunction is not just through elevated LV filling pressure, because even after adjustment for E/e’, LA volume index, and LV mass index, sarcomere gene mutation-related LA dysfunction remained significant. In addition, the findings that PSM were significantly younger than those in the non-sarcomere-related group and the LA function of patients with HCM was significantly lower than controls despite same age and sex distribution, suggest that LA dysfunction in HCM is a unique myopathic process and that sarcomere gene mutation significantly contribute to LA dysfunction. Previous basic studies showed that sarcomere gene mutation directly induces myocyte hypertrophy, myocyte disarray, and fibrosis in the LV in HCM mouse models [25]. Some studies showed that TGF β1-medicated LA fibrosis also developed in sarcomere gene mutation-derived mouse HCM models only with dysfunctional LV [25], which supports our speculation of a common myopathic process in the LV and LA. This finding also suggests that current anti-hypertrophy and anti-fibrotic agents [26] would affect LA structural and functional remodeling, thereby reducing AF. In our study, mitochondria-related mutations were not significantly related to LA dysfunction and LV fibrosis; therefore, phenotypic contribution and myocardial fibrosis is mainly contributed by pathogenic sarcomere gene mutations. Despite the small number of subjects with validated pathogenic sarcomere gene mutations in our study, the degree of genetic contribution to LA function seemed to vary between mutations. According to our study, patients with TNNI3 mutation had worse LA function than patients with MYBPC3 mutation did, whereas mitochondria-related gene mutations showed more benign phenotypes in terms of LA function. Therefore, gene-targeted therapy needs to be individualized.

Prevention of sudden cardiac death due to ventricular tachycardia is a major goal of management of HCM [27]. However, several recent studies showed that the general prognosis of HCM is better than that in previous reports, and active primary prevention has reduced the rate of sudden deaths related to HCM [28]. With the advancing age of patients with HCM, heart failure due to advanced diastolic dysfunction or development of AF become a major concern [29]. In fact, several studies showed that concomitant or newly developed AF is a major prognostic factor for heart failure admission, new-onset stroke, and mortality in HCM [30, 31]. Therefore, prevention of AF should be one of the goals of management of HCM. In this regard, several new attempts for upstream genetic modulation using CRISPR/Cas9, [32] a sarcomere-targeted anti-hypertrophic or anti-fibrotic agent, [26] could have potential preventive effects on the development of AF and heart failure.

This study has several limitations. First, we did not evaluate the myopathic process of the LA itself, such as direct tissue characterization or CMR-derived LA fibrosis measurement [33]. However, measurement of LA-LGE involves extensive variability, and 3D-LGE is not routinely measured in most institutions. Hence, use of the most popular parameters with less measurement variability would be justified. Further studies are warranted to develop a more accurate and reproducible method for assessing the LA myopathic process in HCM. Second, owing to lack of long-term follow-up data, the relationship between LA dysfunction and future development of new-onset AF or heart failure was not proved. However, several previous studies showed that LA structural or functional remodeling was significantly related to poor outcomes, especially to the development of new-onset AF and heart failure [1, 2]. The findings of this study could help elucidate this relationship. In addition, our results attempt to explain the reason behind the worse prognosis of PSM among those with HCM. Third, due to limited number of CMR, LA active emptying fraction did not reach the statistical differences between sarcomere groups. However, LA total emptying fraction was significantly lower in sarcomere mutation groups even in the limited patients, which supports our observation measured by a’ velocity. Larger number of echocardiographic or CMR studies with phasic LA functional analyses warrant further investigation.

Conclusion

Presence of pathogenic sarcomere mutation was significantly related to LA dysfunction independent of LV filling pressure and LAVI, suggesting its contribution to atrial myopathy in HCM. Therefore, our study results may provide an explanation why HCM patients with sarcomere gene mutations had poorer prognosis, regarding higher prevalence of AF and hospitalization due to heart failure [34].

Availability of data and materials

The datasets of the current study are available from the corresponding author on reasonable request.

Abbreviations

LA:

Left atrial

AF:

Atrial fibrillation

HCM:

Hypertrophic cardiomyopathy

a’:

Tissue Doppler imaging-based late mitral annular velocity

LV:

Left ventricular

ApHCM:

Apical hypertrophic cardiomyopathy

nDNA:

Nuclear DNA

mtDNA:

Mitochondrial DNA

CMR:

Cardiac magnetic resonance imaging

LGE:

Late gadolinium enhancement

PSM:

Patients with pathogenic or likely pathogenic sarcomere gene mutations

E:

Early diastolic transmitral inflow velocity

s’:

Tissue Doppler imaging-based systolic mitral annular tissue velocity

GLS:

Global longitudinal strain

References

  1. Hinojar R, Zamorano JL, Fernandez-Mendez M, Esteban A, Plaza-Martin M, Gonzalez-Gomez A, Carbonell A, Rincon LM, Nacher JJJ, Fernandez-Golfin C. Prognostic value of left atrial function by cardiovascular magnetic resonance feature tracking in hypertrophic cardiomyopathy. Int J Cardiovasc Imaging. 2019.

  2. Vasquez N, Ostrander BT, Lu DY, Ventoulis I, Haileselassie B, Goyal S, Greenland GV, Vakrou S, Olgin JE, Abraham TP, et al. Low left atrial strain is associated with adverse outcomes in hypertrophic cardiomyopathy patients. J Am Soc Echocardiography. 2019.

  3. Debonnaire P, Joyce E, Hiemstra Y, Mertens Bart J, Atsma Douwe E, Schalij Martin J, Bax Jeroen J, Delgado V, Marsan Nina A. Left atrial size and function in hypertrophic cardiomyopathy patients and risk of new-onset atrial fibrillation. Circ Arrhythm Electrophysiol. 2017;10(2):e004052.

    Article  PubMed  Google Scholar 

  4. Sivalokanathan S, Zghaib T, Greenland GV, Vasquez N, Kudchadkar SM, Kontari E, Lu DY, Dolores-Cerna K, van der Geest RJ, Kamel IR, et al. Hypertrophic cardiomyopathy patients with paroxysmal atrial fibrillation have a high burden of left atrial fibrosis by cardiac magnetic resonance imaging. JACC Clinical electrophysiology. 2019;5(3):364–75.

    Article  PubMed  Google Scholar 

  5. Ko T. Left atrium as an active component of the pathophysiology in HCM. Int Heart J. 2018;59(5):906–8.

    Article  PubMed  Google Scholar 

  6. Blume GG, Mcleod CJ, Barnes ME, Seward JB, Pellikka PA, Bastiansen PM, Tsang TSM. Left atrial function: physiology, assessment, and clinical implications. Eur Heart J - Cardiovascular Imaging. 2011;12(6):421–30.

    Google Scholar 

  7. Walsh R, Buchan R, Wilk A, John S, Felkin LE, Thomson KL, Chiaw TH, Loong CC, Pua CJ, Raphael C, et al. Defining the genetic architecture of hypertrophic cardiomyopathy: re-evaluating the role of non-sarcomeric genes. Eur Heart J. 2017.

  8. Brandon MC, Lott MT, Nguyen KC, Spolim S, Navathe SB, Baldi P, Wallace DC. MITOMAP: a human mitochondrial genome database--2004 update. Nucleic Acids Res. 2005;33(Database issue):D611–3.

    Article  CAS  PubMed  Google Scholar 

  9. Ingman M. Gyllensten U: mtDB: human mitochondrial genome database, a resource for population genetics and medical sciences. Nucleic Acids Res. 2006;34(Database issue):D749–51.

    Article  CAS  PubMed  Google Scholar 

  10. Chaitanya L, Ralf A, van Oven M, Kupiec T, Chang J, Lagace R, Kayser M. Simultaneous whole mitochondrial genome sequencing with short overlapping amplicons suitable for degraded DNA using the ion torrent personal genome machine. Hum Mutat. 2015;36(12):1236–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Zaragoza MV, Brandon MC, Diegoli M, Arbustini E, Wallace DC. Mitochondrial cardiomyopathies: how to identify candidate pathogenic mutations by mitochondrial DNA sequencing, MITOMASTER and phylogeny. European journal of human genetics : EJHG. 2011;19(2):200–7.

    Article  CAS  PubMed  Google Scholar 

  12. Chung H, Kim Y, Cho SM, Lee HJ, Park CH, Kim JY, Lee SH, Min PK, Yoon YW, Lee BK, et al. Differential contributions of sarcomere and mitochondria-related multigene variants to the endophenotype of hypertrophic cardiomyopathy. Mitochondrion. 2020;53:48–56.

    Article  CAS  PubMed  Google Scholar 

  13. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics (Oxford, England). 2009;25(14):1754–60.

    Article  CAS  Google Scholar 

  14. Yang J, Ding X, Sun X, Tsang SY, Xue H. SAMSVM: a tool for misalignment filtration of SAM-format sequences with support vector machine. J Bioinforma Comput Biol. 2015;13(6):1550025.

    Article  CAS  Google Scholar 

  15. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. 2015;17(5):405–24.

  18. Whiffin N, Minikel E, Walsh R, O'Donnell-Luria AH, Karczewski K, Ing AY, Barton PJR, Funke B, Cook SA, MacArthur D, et al. Using high-resolution variant frequencies to empower clinical genome interpretation. Genetics in medicine : official journal of the American College of Medical Genetics. 2017;19(10):1151–8.

    Article  Google Scholar 

  19. Whiffin N, Walsh R, Govind R, Edwards M, Ahmad M, Zhang X, Tayal U, Buchan R, Midwinter W, Wilk AE, et al. CardioClassifier: disease- and gene-specific computational decision support for clinical genome interpretation. Genetics in medicine : official journal of the American College of Medical Genetics. 2018;20(10):1246–54.

    Article  CAS  Google Scholar 

  20. Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, Flachskampf FA, Foster E, Goldstein SA, Kuznetsova T, et al. Recommendations for Cardiac Chamber Quantification by Echocardiography in Adults: An Update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. J Am Soc Echocardiography. 2015;28(1):1–39 e14.

    Article  Google Scholar 

  21. Sasson Z, Yock PG, Hatle LK, Alderman EL, Popp RL. Doppler echocardiographic determination of the pressure gradient in hypertrophic cardiomyopathy. J Am Coll Cardiol. 1988;11(4):752–6.

    Article  CAS  PubMed  Google Scholar 

  22. Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart. A Statement for Healthcare Professionals From the Cardiac Imaging Committee of the Council on Clinical Cardiology of the American Heart Association. 2002;105(4):539–42.

    Google Scholar 

  23. Park CH, Chung H, Kim Y, Kim JY, Min PK, Lee KA, Yoon YW, Kim TH, Lee BK, Hong BK, et al. Electrocardiography based prediction of hypertrophy pattern and fibrosis amount in hypertrophic cardiomyopathy: comparative study with cardiac magnetic resonance imaging. Int J Cardiovasc Imag. 2018;34(10):1619–28.

    Article  Google Scholar 

  24. Li L, Chen X, Yin G, Yan W, Cui C, Cheng H, Lu M, Zhao S. Early detection of left atrial dysfunction assessed by CMR feature tracking in hypertensive patients. Eur Radiol. 2020;30(2):702–11.

    Article  PubMed  Google Scholar 

  25. Nakajima H, Nakajima Hisako O, Salcher O, Dittiè Andrea S, Dembowsky K, Jing S, Field Loren J. Atrial but not ventricular fibrosis in mice expressing a mutant transforming growth factor-β1 transgene in the heart. Circ Res. 2000;86(5):571–9.

    Article  CAS  PubMed  Google Scholar 

  26. Green EM, Wakimoto H, Anderson RL, Evanchik MJ, Gorham JM, Harrison BC, Henze M, Kawas R, Oslob JD, Rodriguez HM, et al. A small-molecule inhibitor of sarcomere contractility suppresses hypertrophic cardiomyopathy in mice. Science (New York, NY). 2016;351(6273):617–21.

    Article  CAS  Google Scholar 

  27. Marian AJ, Braunwald E. Hypertrophic cardiomyopathy: genetics, pathogenesis, clinical manifestations, diagnosis, and therapy. Circ Res. 2017;121(7):749–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Weissler-Snir A, Adler A, Williams L, Gruner C, Rakowski H. Prevention of sudden death in hypertrophic cardiomyopathy: bridging the gaps in knowledge. Eur Heart J. 2016;38(22):1728–37.

    Google Scholar 

  29. members ATF, Elliott PM, Anastasakis A, Borger MA, Borggrefe M, Cecchi F, Charron P, Hagege AA, Lafont A, Limongelli G, et al. 2014 ESC guidelines on diagnosis and management of hypertrophic cardiomyopathy: the task force for the diagnosis and Management of Hypertrophic Cardiomyopathy of the European Society of Cardiology (ESC). Eur Heart J. 2014;35(39):2733–79.

    Article  Google Scholar 

  30. Kubo T, Kitaoka H, Okawa M, Hirota T, Hayato K, Yamasaki N, Matsumura Y, Yabe T, Takata J, Doi YL. Clinical impact of atrial fibrillation in patients with hypertrophic cardiomyopathy. Results from Kochi RYOMA study. Circulation journal : official journal of the Japanese Circulation Society. 2009;73(9):1599–605.

    Article  Google Scholar 

  31. Guerrero JCC, Jimenez-Baena E, Martinez-Martinez A, Valle-Racero J, Francisco L-P, Lopez-Haldon J, Urbano-Moral J. REVISITING PROGNOSTIC IMPACT OF ATRIAL FIBRILLATION IN HYPERTROPHIC CARDIOMYOPATHY. 2016;67(13 Supplement):1515.

  32. Mosqueira D, Mannhardt I, Bhagwan JR, Lis-Slimak K, Katili P, Scott E, Hassan M, Prondzynski M, Harmer SC, Tinker A, et al. CRISPR/Cas9 editing in human pluripotent stem cell-cardiomyocytes highlights arrhythmias, hypocontractility, and energy depletion as potential therapeutic targets for hypertrophic cardiomyopathy. Eur Heart J. 2018;39(43):3879–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Oakes Robert S, Badger Troy J, Kholmovski Eugene G, Akoum N, Burgon Nathan S, Fish Eric N, Blauer Joshua JE, Rao Swati N, DiBella Edward VR, Segerson Nathan M, et al. Detection and quantification of left atrial structural remodeling with delayed-enhancement magnetic resonance imaging in patients with atrial fibrillation. Circulation. 2009;119(13):1758–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ho CY, Day SM, Ashley EA, Michels M, Pereira AC, Jacoby D, Cirino AL, Fox JC, Lakdawala NK, Ware JS, et al. Genotype and lifetime burden of disease in hypertrophic cardiomyopathy: insights from the Sarcomeric human cardiomyopathy registry (SHaRe). Circulation. 2018;138(14):1387–98.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

None.

Funding

This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2014R1A1A2055872).

Author information

Authors and Affiliations

Authors

Contributions

Hyemoon Chung:Writing - Original Draft, Formal analysis, Investigation; Yoonjung Kim: Investigation, Formal analysis Writing - Original Draft; Chul-Hwan Park: Methodology, Validation; In-Soo Kim, Jong-Youn Kim, Pil-Ki Min, Young Won Yoon: Writing - Review & Editing, Supervision; Tae Hoon Kim: Methodology; Byoung Kwon Lee, Bum-Kee Hong, Se-Joong Rim, Hyuck Moon Kwon:Supervision, Investigation; Kyung-A Lee: Methodology, Investigation, Formal analysis, Writing - Review & Editing; Eui-Young Choi: Conceptualization, Formal analysis, Writing - Review & Editing, Funding acquisition. The author(s) read and approved the final manuscript.

Corresponding authors

Correspondence to Kyung-A Lee or Eui-Young Choi.

Ethics declarations

Ethics approval and consent to participate

The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki, and it was approved by institutional review board of Gangnam Severance Hospital (3–2015-0019). Written informed consent was obtained from all participants.

Consent for publication

All the authors read final version of manuscript and approved publication.

Competing interests

The authors declare that there is no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Method S1.

DNA preparation. Method S2. Library construction and sequencing of the HCM gene panel. Method S3. Library construction and mtDNA sequencing. Method S4. Data analysis of the mitochondrial genome. Method S5. CMR. Table S1. Summary of 82 genes associated with hypertrophic cardiomyopathy. Supplementary Table S2. Likely pathogenic or Pathogenic variants in the 33 sarcomere associated genes classified according to the refined American College of Medical Genetics and Genomics (ACMG) standards and guidelines for inherited cardiac conditions. Supplementary Table S3. Likely pathogenic or pathogenic variants in the 6 non-sarcomere genes and the 44 mitochondria-related nuclear genes. Supplementary Table S4. Non-haplogroup-associated variants with a GenBank frequency < 0.1%.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chung, H., Kim, Y., Park, C. et al. Contribution of sarcomere gene mutations to left atrial function in patients with hypertrophic cardiomyopathy. Cardiovasc Ultrasound 19, 4 (2021). https://doi.org/10.1186/s12947-020-00233-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12947-020-00233-y

Keywords