Skip to main content

Genetic relevance and determinants of mitral leaflet size in hypertrophic cardiomyopathy

Abstract

Background

Whether mitral leaflet elongation is a primary phenotype of hypertrophic cardiomyopathy (HCM) is controversial. We investigated the genetic relevance and determinants of mitral leaflet size by performing extensive gene analyses in patients with HCM.

Methods

Anterior mitral leaflet (AML) lengths were measured in HCM patients (n = 211) and age- and sex-matched controls (n = 30) using echocardiography with hemodynamic and chamber geometric assessments. We analyzed 82 nuclear DNA (8 sarcomeric genes, 74 other HCM-associated genes) and mitochondrial DNA. Cardiac magnetic resonance imaging (CMR) was performed in the 132 HCM patients.

Results

Average indexed AML was significantly longer for HCM than for controls (17.2 ± 2.3 vs. 13.3 ± 1.6 mm/m2, P <  0.001). Average AML length correlated with body surface area (BSA), left ventricular (LV) end-systolic volume (P <  0.001) and LV mass by CMR (P < 0.001). Average indexed AML by BSA of pure-apical HCM was significantly shorter than other typed HCM (16.6 ± 2.0 vs. 17.4 ± 2.4 mm/m2, P = 0.025). Indexed AML was independently correlated with left atrial wall stress. The thin filament mutation group showed larger average AML (31.9 ± 3.8 vs. 29.6 ± 3.8 mm, P = 0.045), but this was not significant with the indexed value. No difference in AML size among subgroups was observed based on the presence of sarcomere protein or mitochondria-related gene variants (P > 0.05).

Conclusion

AML elongation was a unique finding of HCM. However, the leaflet size was more related to chamber geometry and hypertrophy pattern rather than genetic factors within overt HCM.

Peer Review reports

Introduction

Mitral leaflet elongation is related to obstructive hypertrophic cardiomyopathy (HCM) [1]. Although, there are still some controversies whether classical pathogenic sarcomere gene variants independently affect mitral leaflet enlargement [2,3,4]. Previous basic studies have shown that some sarcomere genes affect valvular growth [5, 6]. In addition, a recent study further confirmed that leaflet elongation is a primary phenotype in subclinical HCM patients with pathogenic sarcomere gene mutations [7]. However, these investigations mainly focused on variants in classical pathogenic sarcomere genes such as MYBPC3 and MYH7, regardless of the morphologic phenotype of the hypertrophy pattern [2,3,4, 7]. The contributions of various gene variant groups other than sarcomere genes have not been extensively investigated. In addition, about 40–60% of patients with overt HCM do not have any pathogenic sarcomere gene mutations [4]. Furthermore, regarding mitral leaflet size, there was no comparison study between pathogenic variant group and non-variant group within overt hypertrophy patients. Therefore, it is still unclear whether mitral leaflet size in overt HCM is mainly determined by genotype or by additional geometric or hemodynamic factors. In this study, to reveal the genetic relevance and effects of LV geometric change to mitral leaflet size, an extensive HCM gene panel comprising 82 nuclear DNA (nDNA) genes (8 sarcomere genes, 74 other HCM-associated genes) and mitochondrial DNA (mtDNA) genes was analyzed along with diverse morphologic phenotype definitions and hemodynamic factors.

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 due to insufficient data, follow-up loss, or declining study enrollment. Of these patients, one had poor image quality that precluded measuring mitral leaflet size in any view. Finally, 211 patients underwent genetic testing. The patients enrolled in the study 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. Cases of abnormal papillary muscle insertion to the apex with thickening were also included. Patterns of LV hypertrophy were classified as apical HCM (ApHCM) and non-ApHCM (asymmetrical hypertrophy, diffuse hypertrophy, and focal segmental hypertrophy). ApHCM was classified as pure apical type (hypertrophy confined below the papillary muscle level) and mixed type (mixed pure apical and asymmetrical septal hypertrophy but maximal thickness in the apex) according to echocardiographic findings. All patients underwent screening analysis for Fabry disease and were confirmed negative for the galactosidase alpha (GLA) variant. For comparison, anterior mitral leaflet (AML) size was measured in 30 age- and sex-matched controls. The study protocol was approved by our institutional review board (3–2015-0019), and written informed consent was obtained for each subject.

Genetic testing and analysis

HCM gene panel (nDNA) design

A literature search of the PubMed database was performed for targeted gene selection for the comprehensive HCM-specific panel. It included 82 nDNA genes (33 sarcomere protein genes, 5 phenocopy genes, and 44 nuclear genes linked to mitochondrial cardiomyopathy, Fig. 1 and Additional file 2 Table S1). HCM genes consisted of 8 validated sarcomere genes and 25 putative HCM genes [8].

Fig. 1
figure 1

The comprehensive hypertrophic cardiomyopathy (HCM)-specific panel consisted of 82 nuclear genes including (a) 33 sarcomere-associated genes, (b) 5 phenocopy genes, and (c) 44 mitochondria-related nuclear DNA (nDNA) genes

DNA preparation

The details are described at Additional file 1 Method S1.

Library construction and sequencing of the HCM gene panel and mtDNA

The details are described at Additional file 1 Method S2.

Identification of potential pathogenic mtDNA variants

Non-haplogroup-associated novel and rare variants were evaluated for their potential pathogenicity based on variant location, amino acid change, and evolutionary conservation [9].

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 [10]. SAMTools was used to convert the sequence alignment map file to the BAM format [11]. Sorting and removal of duplications were performed using the Picard tool (http://broadinstitute.github.io/picard/). GATK was used to perform indel realigning and base-quality score re-calibration [12]. Variants were annotated with ANNOVAR [13] and Variant Effect Predictor (http://asia.ensembl.org/info/docs/tools/vep/index.html). Only bases meeting the minimum base quality (≥ 20) from reads are considered. Variants were further filtered with altered allele frequency > 30%, 50× coverage, and population frequency < 0.01 in the 1000 Genome Project, ESP6500, ExAC databases, and the Korean Reference Genome Database (http://152.99.75.168/KRGDB), which was constructed with whole-genome sequencing data from 1100 Korean individuals. Prediction of the potential pathogenicity was performed using Alamut® Visual software (Interactive Biosoftware, Rouen, France), Human Gene Mutation Database professional version, release 2016.1 (http://www.hgmd.org/), ClinVar, and the Atlas of Cardiac Genetic Variation. The impact of missense change was predicted with Align GVGD, SIFT, PolyPhen and MutationTaster. Visual inspection of candidate variants were performed using the Integrative Genomics Viewer (IGV) [14] and we performed Sanger confirmation of pathogenic variants which showed unacceptable quality metrics. Variants were classified based on the American College of Medical Genetics and Genomics standards and guidelines [15, 16].

Data analysis of the mitochondrial genome

The details are described in Additional file 1 Method S3.

Echocardiographic analysis

Comprehensive echo-Doppler evaluation was performed according to current American Society of Echocardiography guidelines [17]. A routine standard echocardiography study was performed to measure systolic and diastolic parameters as follows: LV end-diastolic volume and LV end-systolic volume (LVESV) were measured with the biplane Simpson’s method, and LV ejection fraction was calculated. Left atrial (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. Tissue Doppler interrogation was performed in the septal mitral annulus in apical four-chamber view, followed by measurement of the peak systolic mitral annulus velocity (s′) and early diastolic mitral annulus peak velocity (e′). The ratio of E/e′ was calculated. LV wall thickness was measured in all cross-sectional planes. Continuous wave Doppler was used to measure 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) [18]. It was measured at resting and during Valsalva maneuver. LVOT obstruction was defined as a systolic pressure gradient ≥30 mmHg. Contrast echocardiography was performed in patients with a poorly defined LV border. Mitral regurgitation (MR) degree was defined as trivial, I, II, III, or IV according to regurgitation area.

Cardiac magnetic resonance imaging (CMR)

CMR was performed using a 1.5-T scanner (Magnetom Avanto®; Siemens Medical Solutions, Erlangen, Germany) with a phased array body coil. The LV 2-, 3-, 4-chamber, and short axis views were obtained using cine images with steady-state free precession sequence. The endocardial and epicardial borders were contoured using a semi-automated method (Argus®; Siemens, Germany or Qmass® MR 8.1; Medis, Leiden, the Netherlands); subsequently, the LVEDV and LVESV were measured. To determine the end-diastolic LV mass, the difference between the epicardial and endocardial areas for all slices was multiplied by the slice thickness and section gap, and then multiplied by the specific gravity of the myocardium (1.05 g/mL). Papillary muscle mass was included in the LV cavity and excluded from the LV mass measurements.

Mitral leaflet length measurement

AML lengths were measured in the parasternal long axis (PLX) and apical three-chamber (3CH) views, the same measurements from age- and sex-matched controls were taken for comparison. In both views, AML lengths were measured in diastole, with the leaflets maximally extended and parallel to the anterior septum, and defined as the distance from the most distal extent of the anterior leaflet to its insertion into the posterior aortic wall (Fig. 2). The measurements were repeated in 10 randomly selected patients to assess reproducibility.

Fig. 2
figure 2

Schematic illustration of anterior mitral leaflet distance measurement in the parasternal long axis view (a) and apical 3-chamber view (b)

LV geometry and mitral valve wall stress assessment

For end-systolic mitral valve wall stress assessment, the index from LV side was roughly calculated as (LV end-systolic pressure) × LVESV, where LV end-systolic pressure was calculated as aorta end-systolic pressure (2 × systolic blood pressure + diastolic blood pressure/3) + peak systolic LVOT pressure gradient. E/e′ × LA volume was used for measurement of diastolic mitral valve wall stress index from the LA side (LAWS).

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 reproducibility test, paired sample t-tests and Pearson’s correlation coefficients were determined. All clinical statistical analyses were performed using SPSS version 23.0 statistical package (IBM Corp., Armonk, NY, USA). A two-sided p-value < 0.05 was considered statistically significant.

Results

Baseline characteristics

The mean age of patients was 59 ± 14 years, and 63 (30%) were female. Of them, 49 (23%) had obstructive HCM; 100 patients (47%) had ApHCM, and 64 (64%, 64/100) of these patients had pure-type ApHCM. The AML lengths of patients with HCM, measured in both PLX and 3CH, were significantly longer than those of controls (32.1 ± 4.5 vs. 26.4 ± 2.9 mm in PLX; 28.3 ± 3.6 vs. 23.6 ± 3.1 mm in 3CH; both P < 0.001). Even after being indexed by BSA, patients with HCM had longer AMLs than control (18.3 ± 2.8 vs. 14.0 ± 2.0 mm/m2 in PLX; 16.1 ± 2.4 vs. 12.5 ± 1.8 mm/m2 in 3CH, P < 0.001) (Fig. 3).

Fig. 3
figure 3

Comparison of indexed AML between sarcomere gene mutation positive and mutation negative patients along with comparison to control (a). Comparison of AML-average between thick sarcomere gene variants and thin or regulatory sarcomere gene variants (b). Comparison of indexed AML (c) and prevalence of sarcomere gene variants (d) between pure apical HCM and mixed or non-apical HCM. AML, anterior mitral leaflet length; iAML, indexed AML, HCM, hypertrophic cardiomyopathy; bar represents mean ± standard error

Genetic characteristics

Based on the American College of Medical Genetics and Genomics guideline, [16] 67 of 211 (31.8%) cases 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 variant in HCM genes. We identified homozygous or compound heterozygous variants in MYBPC3 in one patient and co-variants in three patients (two had MYBPC3 and MYH7, and one had MYBPC3 and JPH2). Twenty-seven (13%) patients had probably damaging mtDNA variants, 15 (7%) had mitochondria-related nDNA variants, and 1 had a pathogenic variant in GAA (Additional file 2 Table S2 and Table S3). Patients with ApHCM had a lower prevalence of sarcomere protein gene variants (41.4% vs. 20.8%, P = 0.001) than patients with non-ApHCM.

Correlation with mitral leaflet size

We observed significant correlations between average AML length and BSA (r = 0.329, P < 0.001), LV end-diastolic volume, LVESV (r = 0.341, P < 0.001), LV mass by CMR (r = 0.338, P < 0.001), LV end-systolic wall stress index, LA volume, and maximal thickness. Of these parameters, BSA showed significant and independent correlation with average AML length (β = 0.292, P < 0.001); hence, AML lengths indexed by BSA (iAML) were used for further analyses. The average iAML of pure ApHCM (n = 64) was significantly shorter than that of mixed or non-ApHCM (n = 147; 16.6 ± 2.0 vs. 17.4 ± 2.4 mm/m2, P = 0.025) along with a lower prevalence of sarcomere protein gene mutation (19% vs. 36%, P = 0.006) (Table 1 and Fig. 3). Although the iAML of obstructive HCM was not significantly higher than non-obstructive HCM (17.8 ± 2.9 vs. 17.0 ± 2.1 mm/m2, P = 0.103), iAML was significantly correlated with the LVOT peak pressure gradient (resting and during Valsalva maneuver, r = 0.254, P = 0.001). Average iAML was significantly correlated with age, LA volume, LAWS (r = 0.275, P < 0.001), MR grade (r = 0.373, P < 0.001), and E/e′. In multivariable analysis, LAWS was independently associated with the average iAML (β = 0.212, P = 0.012) (Table 1). The relationships between AML size and echocardiographic parameters remained significant even after exclusion of pure-type ApHCM.

Table 1 Comparison between sarcomere protein gene variant group and non-variant group

Genetic relevance to mitral leaflet length

No significant difference was observed in iAML lengths between the patients with and without pathogenic sarcomere protein gene variants. This non-significance remained after exclusion of patients with pure-type ApHCM. Among patients with sarcomere protein gene variants, there was no significant difference between the MYBPC3 (n = 32) and MYH7 (n = 17) groups (average iAML, 16.8 ± 2.6 vs. 17.8 ± 2.2 mm/m2, P = 0.184). The MYBPC3 mutation group has shorter AML-PLX values (32.4 ± 4.4 vs. 30.6 ± 5.0 mm, P = 0.042), but iAML was not significantly different. The sarcomere protein gene variant group was further divided into thick filament gene variant (n = 52) and thin filament or regulatory gene variant (n = 15) groups. The thin filament group showed larger AML size (31.9 ± 3.8 vs. 29.6 ± 3.8 mm/m2, P = 0.045 with average AML), but this was not significant with the indexed value (P > 0.05) (Table 2). No significant difference in mitral leaflet size was observed upon classification of patients into non-variant, only mitochondrial-related variant, only sarcomere protein gene variant, and both sarcomere and non-sarcomere gene variant groups (P > 0.05) (Table 3).

Table 2 Univariate and multivariable analysis of correlation for anterior mitral leaflet lengths
Table 3 Clinical, anatomic, and hemodynamic characteristics according to the presence of genetic variant

Measurement reproducibility of mitral leaflet length

With regard to intraobserver measurement variability, there was no difference between the first and second AML length measurements in PLX (29.8 ± 1.8 vs. 30.4 ± 2.1 mm, P = 0.161) or 3CH (27.8 ± 2.5 vs. 27.0 ± 2.6 mm, P = 0.138). Correlation between measurements was good for AML length in PLX (r = 0.835, P = 0.003) and 3CH (r = 0.827, P = 0.003).

Discussion

In this study, we confirmed that mitral leaflet size was significantly larger in patients with HCM than in controls. Among these patients, significant and independent correlations were observed between mitral leaflet size and BSA, LV volume, LV mass by CMR, degree of MR and LAWS. An extensive HCM panel showed no independent relationships of genetic factors with mitral leaflet size.

Chamber remodeling and hemodynamic load to mitral leaflet length

Previous studies have shown that hemodynamic load induces mitral leaflet development. In addition, adaptive mitral leaflet growth may be observed in adults [19]. Our results also support the hypothesis that LV end-systolic valve stress from the LV side and LA end-diastolic valve wall stress from the LA side correlate with AML size. We confirmed that leaflet lengths contribute to the development of a trans-LVOT pressure gradient. Our data showing a significant correlation between MR grade and mitral leaflet length suggest that hemodynamic load contributes to leaflet elongation. However, whether increase in the trans-LVOT pressure gradient and MR is a cause or result of AML elongation is unclear. There are principally two types of cells found in mitral leaflet tissue, namely, endothelial cells that cover the surface of the cusps and interstitial cells (ICs) that form a network within the extracellular matrix (ECM) in the body of the cusp [6]. Both cell types exhibit unique functions that are different from those of other endothelial cells and ICs found throughout the body. Valve ICs express a complex pattern of cell surface, cytoskeletal, and muscle proteins. These are able to bind to- and communicate with each other and the ECM. Although, endothelial cells on the outflow and inflow surfaces of the valve differ from one another, [6, 20] valve stresses on both sides would promote leaflet elongation even in overt HCM status as supported by our study results.

Genetic contribution to mitral leaflet length

In our results, patients with HCM had larger iAMLs than age- and sex-matched controls, indicating that genetic factors contribute to mitral leaflet elongation. However, we could not find a direct relationship between leaflet size and sarcomere genes, phenocopy genes, mitochondria-related nDNA or mtDNA in the stage of overt HCM. Previous some studies showed that sarcomere gene mutations promote mitral leaflet elongation in subclinical HCM [2, 7], however after development of myocardial hypertrophy and LV remodeling, hemodynamic and geometrical factors predominantly affect mitral leaflet elongation. Therefore, mitral leaflet elongation could be a secondary phenotype following myocardial hypertrophy. These findings are supported by the previous study results of larger leaflets size in overt HCM patients than in subclinical HCM, in spite of same sarcomere gene mutations [2, 7].

Embryologically, the mitral leaflet originates from the atrio-ventricular annulus, which is adjacent to the basal portion of the left ventricle [20]. According to our study, leaflet size was smaller for pure ApHCM than for others, including asymmetrical basal septal hypertrophy, suggesting that hypertrophy signals affect adjacent AML growth through hemodynamic load or genetic predisposition to basal myocardial hypertrophy. Why the AML is larger in classical HCM than in ApHCM might be explained by genetic contribution, because of more prevalent in mitochondria-related gene or mtDNA mutations and less prevalent in sarcomere gene mutations in ApHCM. Troponin I and T genes, included in thin sarcomere protein gene, were shown to be related to valve IC activation [5]. Thin filament or regulatory genes were more related to mitral leaflet elongation in the present and a previous study [4], which suggests that mitral leaflet size is determined by multiple genes and modified by chronic hemodynamic load and chamber remodeling. However, this relationship was significantly attenuated after controlling for body size or LV size.

A recent study revealed that mutations in DCHS1 gene cause mitral valve prolapse with elongation [21], which suggests that MV specific gene might directly affect mitral leaflet enlargement even in HCM. Therefore, studies for revealing genetic factors not only sarcomere genes but also genes coding MV ICs need to be performed in HCM. This effort would open new therapeutic option for preventing obstructive HCM by genetic manipulation.

Limitations

This study has several limitations. First, only AML size was measured with echocardiography. Although cardiovascular magnetic resonance (CMR) has been used good imaging modality for both AML and posterior mitral leaflet, a review of MV leaflet size showed that CMR cannot perfectly measure the distance. AML usually represents the main pathology in HCM, and leaflet length measurements were performed in two different views, with the average value used to reduce measurement errors. Second, our extensive HCM gene panel and mtDNA analysis did not assess rare variant contribution to mitral leaflet size. Third, the enrolled study population number was not large enough to reach statistical power. The cost of next-generation sequencing is continuously decreasing; hence, extensive gene panel analyses to evaluate multigenic contributions to leaflet size in a large HCM registry is warranted.

Conclusion

Mitral leaflet elongation was a unique finding of HCM. However, the leaflet size was more related to LV geometry and hypertrophy pattern rather than genetic factors within overt HCM, suggesting that mitral leaflet elongation is due to gene–hemodynamic interaction and chamber remodeling.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

3CH:

three-chamber view

AML:

Anterior mitral leaflet

ApHCM:

Apical hypertrophic cardiomyopathy

CMR:

Cardiac magnetic resonance imaging

ECM:

Extracellular matrix

HCM:

Hypertrophic cardiomyopathy

IC:

Interstitial cell

LA:

Left atrial

LV:

Left ventricular

LVESV:

Left ventricular end-systolic volume

LVOT:

Left ventricular outflow tract

MR:

Mitral regurgitation

PLX:

Parasternal long axis view

References

  1. Klues HG, Proschan MA, Dollar AL, Spirito P, Roberts WC, Maron BJ. Echocardiographic assessment of mitral valve size in obstructive hypertrophic cardiomyopathy. Anatomic validation from mitral valve specimen. Circulation. 1993;88(2):548–55.

    Article  CAS  Google Scholar 

  2. Maron MS, Olivotto I, Harrigan C, Appelbaum E, Gibson CM, Lesser JR, et al. Mitral valve abnormalities identified by cardiovascular magnetic resonance represent a primary phenotypic expression of hypertrophic cardiomyopathy. Circulation. 2011;124(1):40–7.

    Article  CAS  Google Scholar 

  3. Tarkiainen M, Sipola P, Jalanko M, Helio T, Laine M, Jarvinen V, et al. Cardiovascular magnetic resonance of mitral valve length in hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2016;18(1):33.

    Article  Google Scholar 

  4. Weissler-Snir A, Hindieh W, Gruner C, Fourey D, Appelbaum E, Rowin E, et al. Lack of Phenotypic Differences by Cardiovascular Magnetic Resonance Imaging in MYH7 (beta-Myosin Heavy Chain)- Versus MYBPC3 (Myosin-Binding Protein C)-Related Hypertrophic Cardiomyopathy. Circ Cardiovasc Imaging. 2017;10(2).

  5. Roy A, Brand NJ, Yacoub MH. Molecular characterization of interstitial cells isolated from human heart valves. J Heart Valve Dis. 2000;9(3):459–64 discussion 64-5.

    CAS  PubMed  Google Scholar 

  6. Chester AH, Taylor PM. Molecular and functional characteristics of heart-valve interstitial cells. Philos Trans R Soc Lond Ser B Biol Sci. 2007;362(1484):1437–43.

    Article  CAS  Google Scholar 

  7. Cirino AL, Groarke JD, Lakdawala NK, Galazka PZ, Ho CY, Thune JJ, et al. Intrinsic mitral valve alterations in hypertrophic cardiomyopathy sarcomere mutation carriers. Eur Heart J Cardiovasc Imaging. 2018;19(10):1109–16.

    Article  Google Scholar 

  8. Walsh R, Buchan R, Wilk A, John S, Felkin LE, Thomson KL, et al. Defining the genetic architecture of hypertrophic cardiomyopathy: re-evaluating the role of non-sarcomeric genes. Eur Heart J. 2017. https://doi.org/10.1093/eurheartj/ehw603.

  9. 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. Eur J Hum Genet. 2011;19(2):200–7.

    Article  CAS  Google Scholar 

  10. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  Google Scholar 

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

  12. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, 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  Google Scholar 

  13. 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  Google Scholar 

  14. Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.

    Article  CAS  Google Scholar 

  15. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, 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. Genet Med. 2015;17(5):405–24.

    Article  Google Scholar 

  16. Whiffin N, Walsh R, Govind R, Edwards M, Ahmad M, Zhang X, et al. CardioClassifier: disease- and gene-specific computational decision support for clinical genome interpretation. Genet Med. 2018;20(10):1246–54.

    Article  CAS  Google Scholar 

  17. Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, 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 Echocardiogr. 2015;28(1):1–39.e14.

    Article  Google Scholar 

  18. 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  Google Scholar 

  19. Rausch MK, Tibayan FA, Craig Miller D, Kuhl E. Evidence of adaptive mitral leaflet growth. J Mech Behav Biomed Mater. 2012;15:208–17.

    Article  Google Scholar 

  20. Levine RA, Hagege AA, Judge DP, Padala M, Dal-Bianco JP, Aikawa E, et al. Mitral valve disease--morphology and mechanisms. Nat Rev Cardiol. 2015;12(12):689–710.

    Article  CAS  Google Scholar 

  21. Durst R, Sauls K, Peal DS. deVlaming a, Toomer K, Leyne M, et al. mutations in DCHS1 cause mitral valve prolapse. Nature. 2015;525:109.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Informed consent

Informed consent was obtained from all individual participants included in the study.

Funding

The authors had full access to the database and take full responsibility for its integrity. All authors have read and agreed to the manuscript as written. 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

HC, YK, KL and EYC made the study design and wrote the manuscript. YK and KL analyzed genetic test. CHP, EYC and THK analyzed CMR images. HC, EYC and SJR analyzed the echocardiography. HC and EYC collected the echocardiographic and clinical data. YWY, JYK, BKL, PKM, BKH, and HMK helped collecting clinical data and provide critical comments to the manuscript. All authors 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

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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.

Additional file 2. Table S1.

Summary of 82 genes associated with hypertrophic cardiomyopathy; Table S2. Likely pathogenic or Pathogenic variants in the 33 sarcomere associated genes classified according to the 2015 American College of Medical Genetics guidelines; Table S3. Likely pathogenic or Pathogenic variants in the 6 non-sarcomere genes and the 44 mitochondria-related nuclear genes.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chung, H., Kim, Y., Park, CH. et al. Genetic relevance and determinants of mitral leaflet size in hypertrophic cardiomyopathy. Cardiovasc Ultrasound 17, 21 (2019). https://doi.org/10.1186/s12947-019-0171-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12947-019-0171-1

Keywords