Abstract

Aims

Coarctation of the aorta (CoA) accounts for 4–8% of congenital heart defects (CHDs) and confers substantial morbidity despite treatment. It is increasingly recognized as a highly heritable condition. The aim of the study was to search for sequence variants that affect the risk of CoA.

Methods and results

We performed a genome-wide association study of CoA among Icelanders (120 cases and 355 166 controls) based on imputed variants identified through whole-genome sequencing. We found association with a rare (frequency = 0.34%) missense mutation p.Arg721Trp in MYH6 (odds ratio = 44.2, P = 5.0 × 10−22), encoding the alpha-heavy chain subunit of cardiac myosin, an essential sarcomere protein. Approximately 20% of individuals with CoA in Iceland carry this mutation. We show that p.Arg721Trp also associates with other CHDs, in particular bicuspid aortic valve. We have previously reported broad effects of p.Arg721Trp on cardiac electrical function and strong association with sick sinus syndrome and atrial fibrillation.

Conclusion

Through a population approach, we found that a rare missense mutation p.Arg721Trp in the sarcomere gene MYH6 has a strong effect on the risk of CoA and explains a substantial fraction of the Icelanders with CoA. This is the first mutation associated with non-familial or sporadic form of CoA at a population level. The p.Arg721Trp in MYH6 causes a cardiac syndrome with highly variable expressivity and emphasizes the importance of sarcomere integrity for cardiac development and function.

Introduction

Coarctation of the aorta (CoA) is the most common birth defect of the aorta with an incidence of about one per 2500 live births.1 It is defined by local narrowing of the proximal descending aorta and/or aortic arch, accompanied by bicuspid aortic valve (BAV) in more than 50% of cases and generally presenting as either neonatal heart failure or hypertension later in life.2 Surgical or interventional treatment considerably improves outcome but risk of premature cardiovascular morbidity and mortality remains despite appropriate therapy.3

Coarctation of the aorta is primarily a non-familial or sporadic disease.4 However, it has been shown to cosegregate in families with left-ventricular outflow tract obstruction (LVOTO) malformations, a mechanistically defined subgroup of congenital heart defects (CHDs) including CoA, BAV, congenital aortic stenosis, and hypoplastic left heart syndrome (HLHS).5 As a group, the LVOTO malformations are markedly heritable (0.71–0.90) and have a high relative risk for first-degree relatives (36.9).6 In addition, around 15% of individuals with CoA occur as part of a recognized genetic syndrome (e.g. 45, X, or Turner).7

Not much is known about genetic causes of non-syndromic CoA. Several studies have found mutations in families with LVOTO malformations and a few instances of sporadic CoA, both with and without concomitant CHDs. The most strongly implicated gene is NOTCH1,8–10 encoding a transmembrane receptor that regulates cell fate during development. Mutations in other genes, including MYH611,12,SMAD6,13,NKX2-5,14 and GATA5,15 have been found in one or few individuals with CoA. The MYH6 mutations were found in two families, one with predisposition to atrial septal defect (ASD)11 and the other to HLHS.12 Some individuals in these families presented with CoA. In addition, knockout in mice of several genes found within copy number variants in individuals with CoA, including MCTP2,16,MATR3,17 and FOXC1,18 have resulted in CoA-like phenotypes.

Methods

GWAS study design

Study samples

The CoA sample set included 120 Icelanders who received the discharge diagnosis of CoA at Landspitali, The National University Hospital (LUH) in Reykjavik, the only tertiary referral centre in Iceland, between 1984 and 2016. The individuals were diagnosed with CoA between the years 1950 and 2016, with most individuals (75%) diagnosed after 1990. The individuals diagnosed with CoA were identified either through diagnostic codes of CoA (ICD-9 code 747.1, ICD-10 code Q25.1) registered between 1990 and 2016 or procedure codes of CoA (WHO codes 1-273, 5-369, 5-382, and 5-387, NOMESCO codes FDJ 00, FDJ 10, FDJ 20, FDJ 30, FDJ 42, and FDJ 96) registered between 1984 and 2016. The diagnoses of CoA were confirmed and detailed phenotypic characteristics (Supplementary material online, Table S1) established through review of electronic and paper medical records at LUH. Coarctation of the aorta was defined as a non-syndromic congenital narrowing of the aorta, the diagnosis of which was confirmed by a cardiologist with echocardiography and/or cardiac catheterization. The individuals used as controls in the CoA GWAS analyses consisted of disease-free individuals randomly drawn from the Icelandic genealogical database and individuals from other genetic studies at deCODE.

In addition to CoA, we included in the study the following samples from the deCODE phenotype database: BAV, ASD, ventricular septal defect (VSD), patent ductus arteriosus (PDA), late onset aortic valve stenosis (AVS), sick sinus syndrome (SSS), atrial fibrillation (AF), heart failure (HF), ischaemic stroke (IS), hypertension (HTN), coronary artery disease (CAD), left atrial diameter (LAD), aortic root diameter (ARD), left ventricular end-diastolic diameter (LVEDD), electrocardiogram (ECG) data, thoracic aortic aneurysm, high-degree atrioventricular block, and hypertrophic cardiomyopathy (Supplementary material online, Supplementary methods).

All DNA samples used in the study are part of deCODE’s biobank established in 1996 and built up since then through various genetic studies at deCODE.

The study was approved by the Icelandic Data Protection Authority and the National Bioethics Committee of Iceland. Study approval numbers were VSN-15-053, VSN-15-016, VSN-15-056, VSN-15-058, VSN-15-114, VSN-15-057, and 10-009-S1. Written informed consent was obtained from all study participants. The study complies with the declaration of Helsinki.

Genotyping, whole-genome sequencing, and imputation

For chip genotyping, 151 677 samples were typed with the Illumina HumanHap300, HumanCNV370, HumanHap610, HumanHap1M, HumanHap660, Omni-1, Omni 2.5, or Omni Express bead chips at deCODE. Long range phasing of all chip-genotyped individuals was performed with methods previously described19,20 (Supplementary material online, Supplementary methods).

The whole genomes of 15 220 Icelanders were sequenced using Illumina technology to a median depth of 35X (Supplementary material online, Supplementary methods). The sequence variants identified in the 15 220 sequenced Icelanders were then imputed into 151 677 Icelanders who had been genotyped with various Illumina single nucleotide polymorphism chips and their genotypes phased using long-range phasing.19,20 The imputation of the sequence variants, identified thorough whole-genome sequencing (WGS), into the chip typed long-range phased individuals was performed with the same model as used by IMPUTE.21 The utilization of long-range phased haplotypes enables accurate imputation of variants with frequency down to approximately 0.02% in this data set. Using genealogic information on Icelanders from The Book of Icelanders,22 the sequence variants were imputed into first and second-degree relatives of chip genotyped individuals (genealogical imputation),23 to further increase the sample size for association analysis and to increase the power to detect associations. We identified 32.5 million high quality sequence variants (all with imputation information >0.8 that mapped to build hg38) that were tested for association with CoA under the multiplicative model.

Association analysis

In the association analysis were 120 individuals diagnosed with CoA and 355 116 individuals as controls, all with imputed genotypes. The sequence variants imputed were identified through WGS of 15 220 Icelanders (n = 33 individuals with CoA and n = 15 187 individuals as controls). Of the individuals diagnosed with CoA, 39 were chip typed and long-range phased and of the individuals who were controls 140 661 were chip typed and long-range phased. These were imputed with the same model as used by IMPUTE. The remaining individuals (n = 81 CoA cases and n = 214 412 controls), were not chip typed themselves but were first or second degree relatives of the chip typed individuals and imputed using genealogical imputation as described in Ref.23

To account for inflation in test statistics due to cryptic relatedness and stratification, we applied the method of linkage disequilibrium score regression24 (Supplementary material online, Supplementary methods). The estimated correction factor was 1.04 for the multiplicative model of the CoA association. To correct for multiple testing we used the weighted Holm–Bonferroni method25 to allocate family wise error rate of 0.05 equally between four annotation-based classes of sequence variants (Supplementary material online, Supplementary methods). When testing the association of p.Arg721Trp with several other cardiac phenotypes, the individuals that served as controls consisted of disease-free individuals randomly drawn from the Icelandic genealogical database and individuals from other genetic studies at deCODE.

Phenotypic differences between carriers and non-carriers of p.Arg721Trp

To analyse if CoA carriers of the p.Arg721Trp mutation differed clinically from non-carrier individuals with CoA, we evaluated the frequencies of various clinical characteristics in these two groups with CoA (see Supplementary material online, Table S2). Fisher’s exact test was used to test for significant difference in the mean frequency of the variants between non-carriers and carriers, and the odds ratio (OR) was calculated as [pa/(1−pa)]/[pc/(1−pc)], where pa and pc are the mean frequencies of the variants in non-carriers and carriers, respectively.

Results

A rare missense mutation in MYH6 associates with coarctation of the aorta

To search for sequence variants that associate with non-syndromic CoA, we performed a GWAS including 120 Icelanders with CoA and 355 116 Icelanders who served as population controls. We observed a genome-wide significant association with CoA at chromosome 14q11 (Figure 1), explained by a rare (allele frequency = 0.34%) missense mutation p.Arg721Trp (c.2161C>T) in MYH6, encoding the alpha myosin heavy chain subunit (αMHC) in cardiac muscle. Alpha myosin heavy chain subunit is a main component of the sarcomere, the basic contractile unit of cardiac muscle.26 P.Arg721Trp associates with CoA with an OR of 44.2 (95% confidence interval 20.5–95.5) and P = 5.01 × 10−22 (genome-wide significance threshold for missense variants was set at 6.5 × 10−8, see Methods)27 (Figure 2, Table 1). None of the genotyped individuals (N = 151 677) were homozygous for the mutation, consistent with its low frequency (1.8 homozygotes expected under Hardy–Weinberg equilibrium). Since we observed no homozygotes, we could not discriminate between the dominant and the multiplicative modes of inheritance.

Manhattan plot of coarctation of the aorta genome-wide association study in Iceland. The P values (−log10) are plotted against their respective positions on each chromosome.
Figure 1

Manhattan plot of coarctation of the aorta genome-wide association study in Iceland. The P values (−log10) are plotted against their respective positions on each chromosome.

Region plot for the association of variants on 14q11 with coarctation of the aorta. Shown is a 1 Mb region on chromosome 14. The strongest association is with the missense variant p.Arg721Trp in MYH6 located at position 23 396 970 on chromosome 14 (chr14: 23 396 970). The nine other variants shown are weakly correlated with p.Arg721Trp, r2 between 0.6–0.4 (green) and 0.4–0.2 (light blue).
Figure 2

Region plot for the association of variants on 14q11 with coarctation of the aorta. Shown is a 1 Mb region on chromosome 14. The strongest association is with the missense variant p.Arg721Trp in MYH6 located at position 23 396 970 on chromosome 14 (chr14: 23 396 970). The nine other variants shown are weakly correlated with p.Arg721Trp, r2 between 0.6–0.4 (green) and 0.4–0.2 (light blue).

The p.Arg721Trp mutation is located in exon 18 (out of 39 exons) of MYH6 and leads to an arginine to tryptophan alteration at amino acid 721 (full-length protein 1939 amino acid) (see Supplementary material online, Figure S1). It is located in the converter domain of αMHC (see Supplementary material online, Figure S1 and S2), a small domain crucial in conveying a conformational change from the active site to the lever arm upon adenosine triphosphate (ATP) hydrolysis.28 It is considered likely that the mutation alters protein function (SIFT = 0, PolyPhen = 0.99, MutationTaster = 0.93), probably by altering the folding of the converter domain.

There were 987 carriers of p.Arg721Trp among the 151 677 chip-typed Icelanders and eight of those (one per 123 carriers) were diagnosed with CoA. In line with low penetrance of the mutation for CoA, p.Arg721Trp carriers diagnosed with CoA did not cluster in families. However, 20% of the 39 chip-typed individuals with CoA carried p.Arg721Trp. Thus, while the penetrance of the mutation for CoA is low, it accounts for a large proportion of individuals with CoA in the Icelandic population.

The p.Arg721Trp mutation is not present in the Exome Variant Server, containing sequence data from 6503 individuals [Exome Variant Server, NHLBI Exome Sequencing Project (ESP), Seattle, WA, USA] (http://evs.gs.washington.edu/EVS/) (August 2016) and one copy was found in The Genome Aggregation Database (gnomAD), holding data from 126 216 exome sequences and 15 136 WGS unrelated individuals.29 The p.Arg721Trp mutation thus appears to be absent from or present at a very low frequency in other populations.

We show the phenotypic characteristics of individuals with CoA in Supplementary material online, Table S1. About half were diagnosed during the first month of life and three quarters during the first year of life. As expected,3 CoA was more common in males than in females (1.6:1). About three quarters of individuals with CoA had other CHDs, most commonly BAV and VSD. Similar to other studies,30 the aortic valve was bicuspid in about half of those with CoA.

To determine whether there are phenotypic differences between carriers (n = 24) and non-carriers (n = 96) of the p.Arg721Trp mutation within the CoA sample set, we evaluated the frequencies of various clinical characteristics in the two groups (see Supplementary material online, Table S2). Carriers were nominally more likely to present with mild rather than more critical and complex forms of CoA (OR = 4.2 and P = 0.023). We observed no other differences.

Association of p.Arg721Trp in MYH6 with other cardiac diseases

We have previously demonstrated that p.Arg721Trp associates strongly with SSS and AF, atrial arrhythmias that are common in the elderly and frequently coexist.31 With larger sample sizes, these associations have become stronger; however, previously reported association with thoracic aortic aneurysm is no longer significant (Table 1). In the context of assessing effects of AF risk variants on cardiac conduction, we have also recently shown that p.Arg721Trp associates with many ECG measures corresponding to a widespread effect on electrical function of the heart32 (see Supplementary material online, Figure S3).

Table 1

Association of p.Arg721Trp with congenital heart defects and various cardiac phenotypes

NaffNcontrOR/effect (95% CI)aP-value
Congenital heart defects
 Coarctation of the aorta120355 11644.2 (20.5 to 95.5)5.0 × 10−22
 Bicuspid aortic valve208293 34610.5 (2.6 to 38.0)7.3 × 10−8
 Ventricular septal defect715357 6414.4 (1.9 to 10.0)3.7 × 10−4
 Patent ductus arteriosus594357 7624.9 (2.1 to 11.6)2.3 × 10−4
 Atrial septal defect657353 0963.3 (1.5 to 7.1)0.0026
Cardiac conditions
 Sick sinus syndrome3310346 0828.7 (6.8 to 11.2)6.2 × 10−64
 Atrial fibrillation13 471374 9392.4 (1.9 to 3.0)1.1 × 10−14
 Aortic valve stenosis2457349 3422.7 (1.8 to 4.0)1.8 × 10−6
 Heart failure10 480353 5081.8 (1.4 to 2.3)2.3 × 10−6
 Ischaemic stroke8948369 6241.5 (1.1 to 2.0)0.0029
 High degree atrioventricular block1303361 9192.1 (1.2 to 3.5)0.0092
 Coronary artery disease37 782318 8451.2 (1.0 to 1.5)0.056
 Hypertrophic cardiomyopathy163239 2930.0 (0.0 to 4.5)0.15
 Thoracic aortic aneurysm353302 4581.8 (0.6 to 5.3)0.31
 Hypertension54 974324 8031.1 (0.9 to 1.3)0.44
Echocardiogram
 Left atrial diameter19 3800.3 (0.1 to 0.5)2.6 × 10−4
 Aortic root diameter19 506−0.1 (−0.2 to 0.1)0.41
 LVEDDb57010.0 (−0.3 to 0.3)0.93
NaffNcontrOR/effect (95% CI)aP-value
Congenital heart defects
 Coarctation of the aorta120355 11644.2 (20.5 to 95.5)5.0 × 10−22
 Bicuspid aortic valve208293 34610.5 (2.6 to 38.0)7.3 × 10−8
 Ventricular septal defect715357 6414.4 (1.9 to 10.0)3.7 × 10−4
 Patent ductus arteriosus594357 7624.9 (2.1 to 11.6)2.3 × 10−4
 Atrial septal defect657353 0963.3 (1.5 to 7.1)0.0026
Cardiac conditions
 Sick sinus syndrome3310346 0828.7 (6.8 to 11.2)6.2 × 10−64
 Atrial fibrillation13 471374 9392.4 (1.9 to 3.0)1.1 × 10−14
 Aortic valve stenosis2457349 3422.7 (1.8 to 4.0)1.8 × 10−6
 Heart failure10 480353 5081.8 (1.4 to 2.3)2.3 × 10−6
 Ischaemic stroke8948369 6241.5 (1.1 to 2.0)0.0029
 High degree atrioventricular block1303361 9192.1 (1.2 to 3.5)0.0092
 Coronary artery disease37 782318 8451.2 (1.0 to 1.5)0.056
 Hypertrophic cardiomyopathy163239 2930.0 (0.0 to 4.5)0.15
 Thoracic aortic aneurysm353302 4581.8 (0.6 to 5.3)0.31
 Hypertension54 974324 8031.1 (0.9 to 1.3)0.44
Echocardiogram
 Left atrial diameter19 3800.3 (0.1 to 0.5)2.6 × 10−4
 Aortic root diameter19 506−0.1 (−0.2 to 0.1)0.41
 LVEDDb57010.0 (−0.3 to 0.3)0.93

Shown are the number of affected individuals and control individuals used in the association analysis for each of the traits.

a

Estimated odds ratio (OR) or the effect in standard deviation and the 95% confidence interval (CI) for the association with p.Arg721Trp.

b

Left ventricular end-diastolic diameter.

Table 1

Association of p.Arg721Trp with congenital heart defects and various cardiac phenotypes

NaffNcontrOR/effect (95% CI)aP-value
Congenital heart defects
 Coarctation of the aorta120355 11644.2 (20.5 to 95.5)5.0 × 10−22
 Bicuspid aortic valve208293 34610.5 (2.6 to 38.0)7.3 × 10−8
 Ventricular septal defect715357 6414.4 (1.9 to 10.0)3.7 × 10−4
 Patent ductus arteriosus594357 7624.9 (2.1 to 11.6)2.3 × 10−4
 Atrial septal defect657353 0963.3 (1.5 to 7.1)0.0026
Cardiac conditions
 Sick sinus syndrome3310346 0828.7 (6.8 to 11.2)6.2 × 10−64
 Atrial fibrillation13 471374 9392.4 (1.9 to 3.0)1.1 × 10−14
 Aortic valve stenosis2457349 3422.7 (1.8 to 4.0)1.8 × 10−6
 Heart failure10 480353 5081.8 (1.4 to 2.3)2.3 × 10−6
 Ischaemic stroke8948369 6241.5 (1.1 to 2.0)0.0029
 High degree atrioventricular block1303361 9192.1 (1.2 to 3.5)0.0092
 Coronary artery disease37 782318 8451.2 (1.0 to 1.5)0.056
 Hypertrophic cardiomyopathy163239 2930.0 (0.0 to 4.5)0.15
 Thoracic aortic aneurysm353302 4581.8 (0.6 to 5.3)0.31
 Hypertension54 974324 8031.1 (0.9 to 1.3)0.44
Echocardiogram
 Left atrial diameter19 3800.3 (0.1 to 0.5)2.6 × 10−4
 Aortic root diameter19 506−0.1 (−0.2 to 0.1)0.41
 LVEDDb57010.0 (−0.3 to 0.3)0.93
NaffNcontrOR/effect (95% CI)aP-value
Congenital heart defects
 Coarctation of the aorta120355 11644.2 (20.5 to 95.5)5.0 × 10−22
 Bicuspid aortic valve208293 34610.5 (2.6 to 38.0)7.3 × 10−8
 Ventricular septal defect715357 6414.4 (1.9 to 10.0)3.7 × 10−4
 Patent ductus arteriosus594357 7624.9 (2.1 to 11.6)2.3 × 10−4
 Atrial septal defect657353 0963.3 (1.5 to 7.1)0.0026
Cardiac conditions
 Sick sinus syndrome3310346 0828.7 (6.8 to 11.2)6.2 × 10−64
 Atrial fibrillation13 471374 9392.4 (1.9 to 3.0)1.1 × 10−14
 Aortic valve stenosis2457349 3422.7 (1.8 to 4.0)1.8 × 10−6
 Heart failure10 480353 5081.8 (1.4 to 2.3)2.3 × 10−6
 Ischaemic stroke8948369 6241.5 (1.1 to 2.0)0.0029
 High degree atrioventricular block1303361 9192.1 (1.2 to 3.5)0.0092
 Coronary artery disease37 782318 8451.2 (1.0 to 1.5)0.056
 Hypertrophic cardiomyopathy163239 2930.0 (0.0 to 4.5)0.15
 Thoracic aortic aneurysm353302 4581.8 (0.6 to 5.3)0.31
 Hypertension54 974324 8031.1 (0.9 to 1.3)0.44
Echocardiogram
 Left atrial diameter19 3800.3 (0.1 to 0.5)2.6 × 10−4
 Aortic root diameter19 506−0.1 (−0.2 to 0.1)0.41
 LVEDDb57010.0 (−0.3 to 0.3)0.93

Shown are the number of affected individuals and control individuals used in the association analysis for each of the traits.

a

Estimated odds ratio (OR) or the effect in standard deviation and the 95% confidence interval (CI) for the association with p.Arg721Trp.

b

Left ventricular end-diastolic diameter.

To further explore the effect of the p.Arg721Trp MYH6 mutation, we tested it for association with additional cardiac phenotypes, including other CHDs, common heart diseases, and several echocardiogram variables (Table 1, Supplementary material online, Table S3 and Figure S3; significance threshold set at P < 0.003 (0.05/17 individual phenotypes tested). The p.Arg721Trp mutation associates with increased risk of several CHDs: BAV, VSD, ASD, and PDA (Table 1). As expected, the strongest association was with BAV (OR = 10.5 and P = 7.3 × 10−8). In addition, the mutation associates with late onset AVS. To assess if p.Arg721Trp associates with CHDs in the absence of diagnosed CoA, we re-tested for association after removing individuals with CoA from the analysis. Although the effect of p.Arg721Trp is consistently weaker, the associations remain (see Supplementary material online, Table S3). We cannot exclude the existence of undiagnosed CoA in these individuals. The mutation also associates with HF and IS and with LAD but not with other variables derived from the echocardiographic data such as ARD or LVEDD (Table 1). The p.Arg721Trp mutation did not associate with HTN or CAD.

Discussion

Through GWAS based on variants identified through WGS, we found a rare missense variant in the sarcomere gene MYH6 that has a strong effect on the risk of CoA in the Icelandic population and explains a substantial fraction of CoA in Icelanders. The same mutation also associates with other CHDs, in particular BAV. It has a widespread effect on cardiac electrical function and associates strongly with atrial arrhythmias, both SSS and AF. This is the first mutation shown to associate with non-familial or sporadic form of CoA at a population level. The p.Arg721Trp mutation appears to be absent from other populations or if present, at a very low frequency. The Icelandic population is a founder population in that a small number of ancestors account for a relatively large proportion of genetic diversity in the current population. Hence, sequence variants that are very rare in more outbred populations, like p.Arg721Trp, may thus be more frequent in Icelanders.33

Myosin is a major component of the sarcomere, the building block of the contractile system of cardiac muscle. Myosin is an ATPase cellular motor protein composed of two heavy chains and two pairs of light chains. The two heavy chains are αMHC and beta myosin heavy chain (βMHC) encoded by MYH6 and MYH7, respectively. Both αMHC and βMHC are expressed throughout the heart during embryonic cardiogenesis and βMHC continues to do so in the adult heart whereas αMHC expression becomes restricted to the atrium.34 Expression of MYH6 has not been detected in the aorta.35

The pathogenesis of CoA is not well understood. One of the main models of CoA pathogenesis, the haemodynamic theory,2,36 maintains that cardiac lesions resulting in decreased left ventricular outflow promote development of CoA by reducing blood flow through the Foetal aorta. The p.Arg721Trp mutation could predispose to CoA by reducing blood flow through the Foetal aorta because of diminished contraction of the developing heart. This hypothesis is supported by overexpression studies in rat cardiomyocytes showing that the p.Arg721Trp mutation impairs sarcomeric structure37 and by our ECG data demonstrating widespread effect of p.Arg721Trp on cardiac electrical function, including in the ventricles. Our hypothesis is compatible with the fact that MYH6 is expressed in the ventricles during the development of the heart but not in the aorta.

Very rare mutations in MYH6, other than p.Arg721Trp, have been linked to various CHDs,11,12,38 particularly familial ASD39,40 and both dilated and hypertrophic cardiomyopathy.41 In all instances, these mutations have been restricted to a few sporadic cases or too few families. In two of these families, one with predisposition to ASD11 and the other to HLHS,12 some of the affected family members had other cardiac defects, including CoA. p.Arg721Trp in MYH6 differs from these rare familial mutations in that it associates with CoA at the population level and explains about 20% of individuals with CoA in Iceland.

The main limitation of the study is the small size of the CoA sample set. A larger set might have facilitated detection of more variants associating with CoA with weaker effects than observed for p.Arg721Trp, and allowed a better estimate of the effect (OR), penetrance and the fraction of CoA cases explained by the mutation.

Conclusion

In conclusion, our findings give insights into the pathophysiology of CoA, supporting the haemodynamic theory of the pathogenesis. Moreover, the pleiotropic effect of p.Arg721Trp in MYH6 suggests it causes a cardiac syndrome with highly variable expressivity that is difficult to understand clinically without sequence information. Furthermore, these data emphasize the importance of sarcomere integrity for cardiac development and function.

Acknowledgements

We thank all the study subjects for their valuable participation as well as our colleagues who contributed to data collection, sample handling, and genotyping.

Funding

deCODE genetics/Amgen, Inc.

Conflict of interest: The following authors affiliated with deCODE Genetics/Amgen, Inc. are employed by the company: T.B., R.B.T., G.S., P.S., G.L.N., A.H., S.G., A.M., I.J., D.O.A., D.F.G., U.T., H.H., K.S.

Footnotes

See page 3250 for the editorial comment on this article (doi: 10.1093/eurheartj/ehy230)

References

1

Hoffman
JI
,
Kaplan
S.
The incidence of congenital heart disease
.
J Am Coll Cardiol
2002
;
39
:
1890
1900
.

2

Allen
HD
,
Driscoll
DJ
,
Shaddy
RE
,
Feltes
TF
.
Moss & Adams' Heart Disease in Infants, Children, and Adolescents: Including the Fetus and Young Adult
. 8th ed.
Philadelphia
:
Lippincott Wilkins & Wilkins
;
2013
.

3

Aboulhosn
J
,
Child
JS.
Left ventricular outflow obstruction: subaortic stenosis, bicuspid aortic valve, supravalvar aortic stenosis, and coarctation of the aorta
.
Circulation
2006
;
114
:
2412
2422
.

4

Keane
J
,
Fyler
D
,
Lock
J.
Nadaś Pediatric Cardiology
. 2nd ed.
Philadelphia
:
Saunders
;
2006
.

5

Loffredo
CA
,
Chokkalingam
A
,
Sill
AM
,
Boughman
JA
,
Clark
EB
,
Scheel
J
,
Brenner
JI.
Prevalence of congenital cardiovascular malformations among relatives of infants with hypoplastic left heart, coarctation of the aorta, and d-transposition of the great arteries
.
Am J Med Genet A
2004
;
124a
:
225
230
.

6

McBride
KL
,
Pignatelli
R
,
Lewin
M
,
Ho
T
,
Fernbach
S
,
Menesses
A
,
Lam
W
,
Leal
SM
,
Kaplan
N
,
Schliekelman
P
,
Towbin
JA
,
Belmont
JW.
Inheritance analysis of congenital left ventricular outflow tract obstruction malformations: segregation, multiplex relative risk, and heritability
.
Am J Med Genet A
2005
;
134a
:
180
186
.

7

van der Bom
T
,
Zomer
AC
,
Zwinderman
AH
,
Meijboom
FJ
,
Bouma
BJ
,
Mulder
BJ.
The changing epidemiology of congenital heart disease
.
Nat Rev Cardiol
2011
;
8
:
50
60
.

8

McBride
KL
,
Riley
MF
,
Zender
GA
,
Fitzgerald-Butt
SM
,
Towbin
JA
,
Belmont
JW
,
Cole
SE.
NOTCH1 mutations in individuals with left ventricular outflow tract malformations reduce ligand-induced signaling
.
Hum Mol Genet
2008
;
17
:
2886
2893
.

9

Kerstjens-Frederikse
WS
,
van de Laar
IM
,
Vos
YJ
,
Verhagen
JM
,
Berger
RM
,
Lichtenbelt
KD
,
Klein Wassink-Ruiter
JS
,
van der Zwaag
PA
,
Du Marchie Sarvaas
GJ
,
Bergman
KA
,
Bilardo
CM
,
Roos-Hesselink
JW
,
Janssen
JH
,
Frohn-Mulder
IM
,
van Spaendonck-Zwarts
KY
,
van Melle
JP
,
Hofstra
RM
,
Wessels
MW.
Cardiovascular malformations caused by NOTCH1 mutations do not keep left: data on 428 probands with left-sided CHD and their families
.
Genet Med
2016
;
18
:
914
923
.

10

Freylikhman
O
,
Tatarinova
T
,
Smolina
N
,
Zhuk
S
,
Klyushina
A
,
Kiselev
A
,
Moiseeva
O
,
Sjoberg
G
,
Malashicheva
A
,
Kostareva
A.
Variants in the NOTCH1 gene in patients with aortic coarctation
.
Congenit Heart Dis
2014
;
9
:
391
396
.

11

Arrington
CB
,
Bleyl
SB
,
Matsunami
N
,
Bonnell
GD
,
Otterud
BE
,
Nielsen
DC
,
Stevens
J
,
Levy
S
,
Leppert
MF
,
Bowles
NE.
Exome analysis of a family with pleiotropic congenital heart disease
.
Circ Cardiovasc Genet
2012
;
5
:
175
182
.

12

Tomita-Mitchell
A
,
Stamm
KD
,
Mahnke
DK
,
Kim
MS
,
Hidestrand
PM
,
Liang
HL
,
Goetsch
MA
,
Hidestrand
M
,
Simpson
P
,
Pelech
AN
,
Tweddell
JS
,
Benson
DW
,
Lough
JW
,
Mitchell
ME.
Impact of MYH6 variants in hypoplastic left heart syndrome
.
Physiol Genomics
2016
;
48
:
912
921
.

13

Tan
HL
,
Glen
E
,
Topf
A
,
Hall
D
,
O'Sullivan
JJ
,
Sneddon
L
,
Wren
C
,
Avery
P
,
Lewis
RJ
,
ten Dijke
P
,
Arthur
HM
,
Goodship
JA
,
Keavney
BD.
Nonsynonymous variants in the SMAD6 gene predispose to congenital cardiovascular malformation
.
Hum Mutat
2012
;
33
:
720
727
.

14

McElhinney
DB
,
Geiger
E
,
Blinder
J
,
Benson
DW
,
Goldmuntz
E.
NKX2.5 mutations in patients with congenital heart disease
.
J Am Coll Cardiol
2003
;
42
:
1650
1655
.

15

Bonachea
EM
,
Chang
SW
,
Zender
G
,
LaHaye
S
,
Fitzgerald-Butt
S
,
McBride
KL
,
Garg
V.
Rare GATA5 sequence variants identified in individuals with bicuspid aortic valve
.
Pediatr Res
2014
;
76
:
211
216
.

16

Lalani
SR
,
Ware
SM
,
Wang
X
,
Zapata
G
,
Tian
Q
,
Franco
LM
,
Jiang
Z
,
Bucasas
K
,
Scott
DA
,
Campeau
PM
,
Hanchard
N
,
Umana
L
,
Cast
A
,
Patel
A
,
Cheung
SW
,
McBride
KL
,
Bray
M
,
Craig Chinault
A
,
Boggs
BA
,
Huang
M
,
Baker
MR
,
Hamilton
S
,
Towbin
J
,
Jefferies
JL
,
Fernbach
SD
,
Potocki
L
,
Belmont
JW.
MCTP2 is a dosage-sensitive gene required for cardiac outflow tract development
.
Hum Mol Genet
2013
;
22
:
4339
4348
.

17

Quintero-Rivera
F
,
Xi
QJ
,
Keppler-Noreuil
KM
,
Lee
JH
,
Higgins
AW
,
Anchan
RM
,
Roberts
AE
,
Seong
IS
,
Fan
X
,
Lage
K
,
Lu
LY
,
Tao
J
,
Hu
X
,
Berezney
R
,
Gelb
BD
,
Kamp
A
,
Moskowitz
IP
,
Lacro
RV
,
Lu
W
,
Morton
CC
,
Gusella
JF
,
Maas
RL.
MATR3 disruption in human and mouse associated with bicuspid aortic valve, aortic coarctation and patent ductus arteriosus
.
Hum Mol Genet
2015
;
24
:
2375
2389
.

18

Sanchez-Castro
M
,
Eldjouzi
H
,
Charpentier
E
,
Busson
PF
,
Hauet
Q
,
Lindenbaum
P
,
Delasalle-Guyomarch
B
,
Baudry
A
,
Pichon
O
,
Pascal
C
,
Lefort
B
,
Bajolle
F
,
Pezard
P
,
Schott
JJ
,
Dina
C
,
Redon
R
,
Gournay
V
,
Bonnet
D
,
Le Caignec
C.
Search for rare copy-number variants in congenital heart defects identifies novel candidate genes and a potential role for FOXC1 in patients with coarctation of the aorta
.
Circ Cardiovasc Genet
2016
;
9
:
86
94
.

19

Kong
A
,
Masson
G
,
Frigge
ML
,
Gylfason
A
,
Zusmanovich
P
,
Thorleifsson
G
,
Olason
PI
,
Ingason
A
,
Steinberg
S
,
Rafnar
T
,
Sulem
P
,
Mouy
M
,
Jonsson
F
,
Thorsteinsdottir
U
,
Gudbjartsson
DF
,
Stefansson
H
,
Stefansson
K.
Detection of sharing by descent, long-range phasing and haplotype imputation
.
Nat Genet
2008
;
40
:
1068
1075
.

20

Kong
A
,
Steinthorsdottir
V
,
Masson
G
,
Thorleifsson
G
,
Sulem
P
,
Besenbacher
S
,
Jonasdottir
A
,
Sigurdsson
A
,
Kristinsson
KT
,
Jonasdottir
A
,
Frigge
ML
,
Gylfason
A
,
Olason
PI
,
Gudjonsson
SA
,
Sverrisson
S
,
Stacey
SN
,
Sigurgeirsson
B
,
Benediktsdottir
KR
,
Sigurdsson
H
,
Jonsson
T
,
Benediktsson
R
,
Olafsson
JH
,
Johannsson
OT
,
Hreidarsson
AB
,
Sigurdsson
G
,
Ferguson-Smith
AC
,
Gudbjartsson
DF
,
Thorsteinsdottir
U
,
Stefansson
K.
Parental origin of sequence variants associated with complex diseases
.
Nature
2009
;
462
:
868
874
.

21

Marchini
J
,
Howie
B
,
Myers
S
,
McVean
G
,
Donnelly
P.
A new multipoint method for genome-wide association studies by imputation of genotypes
.
Nat Genet
2007
;
39
:
906
913
.

22

Gudbjartsson
DF
,
Sulem
P
,
Helgason
H
,
Gylfason
A
,
Gudjonsson
SA
,
Zink
F
,
Oddson
A
,
Magnusson
G
,
Halldorsson
BV
,
Hjartarson
E
,
Sigurdsson
GT
,
Kong
A
,
Helgason
A
,
Masson
G
,
Magnusson
OT
,
Thorsteinsdottir
U
,
Stefansson
K.
Sequence variants from whole genome sequencing a large group of Icelanders
.
Sci Data
2015
;
2
:
150011.

23

Styrkarsdottir
U
,
Thorleifsson
G
,
Sulem
P
,
Gudbjartsson
DF
,
Sigurdsson
A
,
Jonasdottir
A
,
Jonasdottir
A
,
Oddsson
A
,
Helgason
A
,
Magnusson
OT
,
Walters
GB
,
Frigge
ML
,
Helgadottir
HT
,
Johannsdottir
H
,
Bergsteinsdottir
K
,
Ogmundsdottir
MH
,
Center
JR
,
Nguyen
TV
,
Eisman
JA
,
Christiansen
C
,
Steingrimsson
E
,
Jonasson
JG
,
Tryggvadottir
L
,
Eyjolfsson
GI
,
Theodors
A
,
Jonsson
T
,
Ingvarsson
T
,
Olafsson
I
,
Rafnar
T
,
Kong
A
,
Sigurdsson
G
,
Masson
G
,
Thorsteinsdottir
U
,
Stefansson
K.
Nonsense mutation in the LGR4 gene is associated with several human diseases and other traits
.
Nature
2013
;
497
:
517
520
.

24

Bulik-Sullivan
BK
,
Loh
PR
,
Finucane
HK
,
Ripke
S
,
Yang
J
,
Patterson
N
,
Daly
MJ
,
Price
AL
,
Neale
BM.
LD Score regression distinguishes confounding from polygenicity in genome-wide association studies
.
Nat Genet
2015
;
47
:
291
295
.

25

Holm
S.
A simple sequentially rejective multiple test procedure
.
Scand J Stat
1979
;
6
:
65
70
.

26

Fuster
V
,
Harrington
RA
,
Narula
J
,
Eapen
ZJ.
Hurst’s the Heart
. 14th ed.
New York City
:
McGraw-Hill Education
;
2017
.

27

Sveinbjornsson
G
,
Albrechtsen
A
,
Zink
F
,
Gudjonsson
SA
,
Oddson
A
,
Masson
G
,
Holm
H
,
Kong
A
,
Thorsteinsdottir
U
,
Sulem
P
,
Gudbjartsson
DF
,
Stefansson
K.
Weighting sequence variants based on their annotation increases power of whole-genome association studies
.
Nat Genet
2016
;
48
:
314
317
.

28

Sweeney
HL
,
Houdusse
A.
Myosin VI rewrites the rules for myosin motors
.
Cell
2010
;
141
:
573
582
.

29

Lek
M
,
Karczewski
KJ
,
Minikel
EV
,
Samocha
KE
,
Banks
E
,
Fennell
T
,
O'Donnell-Luria
AH
,
Ware
JS
,
Hill
AJ
,
Cummings
BB
,
Tukiainen
T
,
Birnbaum
DP
,
Kosmicki
JA
,
Duncan
LE
,
Estrada
K
,
Zhao
F
,
Zou
J
,
Pierce-Hoffman
E
,
Berghout
J
,
Cooper
DN
,
Deflaux
N
,
DePristo
M
,
Do
R
,
Flannick
J
,
Fromer
M
,
Gauthier
L
,
Goldstein
J
,
Gupta
N
,
Howrigan
D
,
Kiezun
A
,
Kurki
MI
,
Moonshine
AL
,
Natarajan
P
,
Orozco
L
,
Peloso
GM
,
Poplin
R
,
Rivas
MA
,
Ruano-Rubio
V
,
Rose
SA
,
Ruderfer
DM
,
Shakir
K
,
Stenson
PD
,
Stevens
C
,
Thomas
BP
,
Tiao
G
,
Tusie-Luna
MT
,
Weisburd
B
,
Won
HH
,
Yu
D
,
Altshuler
DM
,
Ardissino
D
,
Boehnke
M
,
Danesh
J
,
Donnelly
S
,
Elosua
R
,
Florez
JC
,
Gabriel
SB
,
Getz
G
,
Glatt
SJ
,
Hultman
CM
,
Kathiresan
S
,
Laakso
M
,
McCarroll
S
,
McCarthy
MI
,
McGovern
D
,
McPherson
R
,
Neale
BM
,
Palotie
A
,
Purcell
SM
,
Saleheen
D
,
Scharf
JM
,
Sklar
P
,
Sullivan
PF
,
Tuomilehto
J
,
Tsuang
MT
,
Watkins
HC
,
Wilson
JG
,
Daly
MJ
,
MacArthur
DG.
Analysis of protein-coding genetic variation in 60,706 humans
.
Nature
2016
;
536
:
285
291
.

30

Teo
LL
,
Cannell
T
,
Babu-Narayan
SV
,
Hughes
M
,
Mohiaddin
RH.
Prevalence of associated cardiovascular abnormalities in 500 patients with aortic coarctation referred for cardiovascular magnetic resonance imaging to a tertiary center
.
Pediatr Cardiol
2011
;
32
:
1120
1127
.

31

Holm
H
,
Gudbjartsson
DF
,
Sulem
P
,
Masson
G
,
Helgadottir
HT
,
Zanon
C
,
Magnusson
OT
,
Helgason
A
,
Saemundsdottir
J
,
Gylfason
A
,
Stefansdottir
H
,
Gretarsdottir
S
,
Matthiasson
SE
,
Thorgeirsson
GM
,
Jonasdottir
A
,
Sigurdsson
A
,
Stefansson
H
,
Werge
T
,
Rafnar
T
,
Kiemeney
LA
,
Parvez
B
,
Muhammad
R
,
Roden
DM
,
Darbar
D
,
Thorleifsson
G
,
Walters
GB
,
Kong
A
,
Thorsteinsdottir
U
,
Arnar
DO
,
Stefansson
K.
A rare variant in MYH6 is associated with high risk of sick sinus syndrome
.
Nat Genet
2011
;
43
:
316
320
.

32

Thorolfsdottir
RB
,
Sveinbjornsson
G
,
Sulem
P
,
Helgadottir
A
,
Gretarsdottir
S
,
Benonisdottir
S
,
Magnusdottir
A
,
Davidsson
OB
,
Rajamani
S
,
Roden
DM
,
Darbar
D
,
Pedersen
TR
,
Sabatine
MS
,
Jonsdottir
I
,
Arnar
DO
,
Thorsteinsdottir
U
,
Gudbjartsson
DF
,
Holm
H
,
Stefansson
K.
A missense variant in PLEC increases risk of atrial fibrillation
.
J Am Coll Cardiol
2017
;
70
:
2157
2168
.

33

Gudbjartsson
DF
,
Helgason
H
,
Gudjonsson
SA
,
Zink
F
,
Oddson
A
,
Gylfason
A
,
Besenbacher
S
,
Magnusson
G
,
Halldorsson
BV
,
Hjartarson
E
,
Sigurdsson
GT
,
Stacey
SN
,
Frigge
ML
,
Holm
H
,
Saemundsdottir
J
,
Helgadottir
HT
,
Johannsdottir
H
,
Sigfusson
G
,
Thorgeirsson
G
,
Sverrisson
JT
,
Gretarsdottir
S
,
Walters
GB
,
Rafnar
T
,
Thjodleifsson
B
,
Bjornsson
ES
,
Olafsson
S
,
Thorarinsdottir
H
,
Steingrimsdottir
T
,
Gudmundsdottir
TS
,
Theodors
A
,
Jonasson
JG
,
Sigurdsson
A
,
Bjornsdottir
G
,
Jonsson
JJ
,
Thorarensen
O
,
Ludvigsson
P
,
Gudbjartsson
H
,
Eyjolfsson
GI
,
Sigurdardottir
O
,
Olafsson
I
,
Arnar
DO
,
Magnusson
OT
,
Kong
A
,
Masson
G
,
Thorsteinsdottir
U
,
Helgason
A
,
Sulem
P
,
Stefansson
K.
Large-scale whole-genome sequencing of the Icelandic population
.
Nat Genet
2015
;
47
:
435
444
.

34

Franco
D
,
Lamers
WH
,
Moorman
AF.
Patterns of expression in the developing myocardium: towards a morphologically integrated transcriptional model
.
Cardiovasc Res
1998
;
38
:
25
53
.

35

GTEx Consortium
.
The genotype-tissue expression (GTEx) project
.
Nat Genet
2013
;
45
:
580
585
.

36

Rudolph
AM
,
Heymann
MA
,
Spitznas
U.
Hemodynamic considerations in the development of narrowing of the aorta
.
Am J Cardiol
1972
;
30
:
514
525
.

37

Ishikawa
T
,
Jou
CJ
,
Nogami
A
,
Kowase
S
,
Arrington
CB
,
Barnett
SM
,
Harrell
DT
,
Arimura
T
,
Tsuji
Y
,
Kimura
A
,
Makita
N.
Novel mutation in the alpha-myosin heavy chain gene is associated with sick sinus syndrome
.
Circ Arrhythm Electrophysiol
2015
;
8
:
400
408
.

38

Granados-Riveron
JT
,
Ghosh
TK
,
Pope
M
,
Bu'Lock
F
,
Thornborough
C
,
Eason
J
,
Kirk
EP
,
Fatkin
D
,
Feneley
MP
,
Harvey
RP
,
Armour
JA
,
David Brook
J.
Alpha-cardiac myosin heavy chain (MYH6) mutations affecting myofibril formation are associated with congenital heart defects
.
Hum Mol Genet
2010
;
19
:
4007
4016
.

39

Ching
YH
,
Ghosh
TK
,
Cross
SJ
,
Packham
EA
,
Honeyman
L
,
Loughna
S
,
Robinson
TE
,
Dearlove
AM
,
Ribas
G
,
Bonser
AJ
,
Thomas
NR
,
Scotter
AJ
,
Caves
LS
,
Tyrrell
GP
,
Newbury-Ecob
RA
,
Munnich
A
,
Bonnet
D
,
Brook
JD.
Mutation in myosin heavy chain 6 causes atrial septal defect
.
Nat Genet
2005
;
37
:
423
428
.

40

Posch
MG
,
Waldmuller
S
,
Müller
M
,
Scheffold
T
,
Fournier
D
,
Andrade-Navarro
MA
,
De Geeter
B
,
Guillaumont
S
,
Dauphin
C
,
Yousseff
D
,
Schmitt
KR
,
Perrot
A
,
Berger
F
,
Hetzer
R
,
Bouvagnet
P
,
Özcelik
C.
Cardiac alpha-myosin (MYH6) is the predominant sarcomeric disease gene for familial atrial septal defects
.
PLoS One
2011
;
6
:
e28872.

41

Carniel
E
,
Taylor
MR
,
Sinagra
G
,
Di Lenarda
A
,
Ku
L
,
Fain
PR
,
Boucek
MM
,
Cavanaugh
J
,
Miocic
S
,
Slavov
D
,
Graw
SL
,
Feiger
J
,
Zhu
XZ
,
Dao
D
,
Ferguson
DA
,
Bristow
MR
,
Mestroni
L.
Alpha-myosin heavy chain: a sarcomeric gene associated with dilated and hypertrophic phenotypes of cardiomyopathy
.
Circulation
2005
;
112
:
54
59
.

Author notes

Thorsteinn Bjornsson and Rosa B. Thorolfsdottir authors contributed equally to this work.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]

Supplementary data