Protein-altering and regulatory genetic variants near GATA4 implicated in bicuspid aortic valve

Abstract

Bicuspid aortic valve (BAV) is a heritable congenital heart defect and an important risk factor for valvulopathy and aortopathy. Here we report a genome-wide association scan of 466 BAV cases and 4,660 age, sex and ethnicity-matched controls with replication in up to 1,326 cases and 8,103 controls. We identify association with a noncoding variant 151 kb from the gene encoding the cardiac-specific transcription factor, GATA4, and near-significance for p.Ser377Gly in GATA4. GATA4 was interrupted by CRISPR-Cas9 in induced pluripotent stem cells from healthy donors. The disruption of GATA4 significantly impaired the transition from endothelial cells into mesenchymal cells, a critical step in heart valve development.

ARTICLE

Received 5 Oct 2016 | Accepted 31 Mar 2017 | Published 25 May 2017

DOI: 10.1038/ncomms15481 OPEN

Protein-altering and regulatory genetic variants near GATA4 implicated in bicuspid aortic valve

Bo Yang1,2,*, Wei Zhou3,*, Jiao Jiao1,*, Jonas B. Nielsen4, Michael R. Mathis5, Mahyar Heydarpour6,

Guillaume Lettre7,8, Lasse Folkersen9,10, Siddharth Prakash11, Claudia Schurmann12, Lars Fritsche13,14, Gregory A. Farnum3, Maoxuan Lin4, Mohammad Othman15, Whitney Hornsby2, Anisa Driscoll2, Alexandra Levasseur2, Marc Thomas2, Linda Farhat2, Marie-Pierre Dub7,8, Eric M. Isselbacher6, Anders Franco-Cereceda16, Dong-chuan Guo11, Erwin P. Bottinger12, G. Michael Deeb1,2, Anna Booher2,4, Sachin Kheterpal5, Y. Eugene Chen1,2,4, Hyun Min Kang13, Jacob Kitzman3,17, Heather J. Cordell18,

Bernard D. Keavney19,20, Judith A. Goodship18, Santhi K. Ganesh2,4,17, Gonalo Abecasis13, Kim A. Eagle2,4, Alan P. Boyle3,17, Ruth J.F. Loos12,21, Per Eriksson9,**, Jean-Claude Tardif7,8,**, Chad M. Brummett5,**, Dianna M. Milewicz11,**, Simon C. Body6,** & Cristen J. Willer2,3,4,17,**

Bicuspid aortic valve (BAV) is a heritable congenital heart defect and an important risk factor for valvulopathy and aortopathy. Here we report a genome-wide association scan of 466 BAV cases and 4,660 age, sex and ethnicity-matched controls with replication in up to 1,326 cases and 8,103 controls. We identify association with a noncoding variant 151 kb from the gene encoding the cardiac-specic transcription factor, GATA4, and near-signicance forp.Ser377Gly in GATA4. GATA4 was interrupted by CRISPR-Cas9 in induced pluripotent stem cells from healthy donors. The disruption of GATA4 signicantly impaired the transition from endothelial cells into mesenchymal cells, a critical step in heart valve development.

1 Department of Cardiac Surgery, University of Michigan, Ann Arbor, Michigan 48109, USA. 2 Frankel Cardiovascular Center, University of Michigan,Ann Arbor, Michigan 48109, USA. 3 Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 48109, USA.

4 Division of Cardiovascular Medicine, Department of Internal Medicine, University of Michigan, Ann Arbor, Michigan 48109, USA. 5 Department of Anesthesiology, University of Michigan, Ann Arbor, Michigan 48109, USA. 6 Department of Anesthesiology, Perioperative, and Pain Medicine, Brigham and Womens Hospital, Harvard Medical School, Boston, Massachusetts 02115, USA. 7 Montreal Heart Institute, Montreal, Quebec, Canada HIT 1C8.

8 Department of Medicine, Universit de Montral, Montreal, Quebec, Canada QC H3T 1J4. 9 Cardiovascular Medicine Unit, Center for Molecular Medicine, Department of Medicine, Karolinska University Hospital Solna, Karolinska Institutet, Stockholm SE-171 76, Sweden. 10 Center for Biological Sequence Analysis, Technical University of Denmark, Copenhagen DK-2800, Denmark. 11 Department of Internal Medicine, Division of Medical Genetics, University of Texas Health Science Center at Houston McGovern Medical School, Houston, Texas 77030, USA. 12 The Charles Bronfman Institute for Personalized Medicine, The Icahn School of Medicine at Mount Sinai, New York, New York 10029, USA. 13 Department of Biostatistics, University of Michigan, Ann Arbor, Michigan 48109, USA. 14 Norwegian University of Science and Technology, Trondheim 7491, Norway. 15 Department of Ophthalmology and Visual Sciences, University of Michigan, Ann Arbor, Michigan 48105, USA. 16 Cardiothoracic Surgery Unit, Department of Molecular Medicine and Surgery, Karolinska University Hospital Solna, Karolinska Institutet, Stockholm SE-171 76, Sweden. 17 Department of Human Genetics, University of Michigan, Ann Arbor, Michigan 48109, USA. 18 Institute of Genetic Medicine, Newcastle University, Newcastle Upon Tyne NE1 3BZ, UK. 19 Division of Cardiovascular Sciences, Faculty of Biology, Medicine and Health, The University of Manchester, Manchester M13 9PL, UK. 20 Manchester Heart Centre, Central Manchester University Hospitals NHS Foundation Trust, Manchester M13 9WL, UK. 21 The Mindich Child Health Development Institute, The Icahn School of Medicine at Mount Sinai, New York, New York 10029, USA. * These authors contributed equally to this work. ** These authors jointly supervised this work. Correspondence and requests for materials should be addressed to B.Y. (email: mailto:[email protected]

Web End [email protected] ) or to C.J.W. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481

Bicuspid aortic valve (BAV) is a congenital aortic valve defect characterized by fusion of two of the normal three leaets. With a prevalence of B1% in the population and a feature

of some rare connective-tissue syndromes, BAV is the most common cardiovascular malformation in humans1,2. BAV is associated with serious consequences: 3070% of those with BAV will develop dilated thoracic aorta3; 1571% of BAV patients develop aortic valve stenosis depending on age group and individuals with BAV have a 50-fold higher risk of severe aortic valve stenosis4; and up to 47% of BAV patients develop aortic valve incompetence5. The presence of a BAV confers an eightfold increased risk of aortic dissection, which carries very high mortality6. Overall, 27% of BAV patients will require surgical intervention to either replace their aortic valve or aorta for aortic aneurysm and dissection7. BAV accounts for B40% of the 450,000 aortic valve replacements performed in the United States each year8.

BAV is moderately heritable, with estimates ranging from 20 to 89% (refs 911). Despite the prevalence, importance and heritability of BAV, its genetic origins remain elusive. Previous genetic studies of BAV have focused primarily on linkage analysis in families9,12 or sequencing candidate genes in cases13 under a hypothesis of Mendelian inheritance. Only one previous genome-wide association study (GWAS) for BAV has been published in a limited number of cases (n 68; ref. 14), which did not identify

any genome-wide signicant results. The only gene in which variants have been identied to cause BAV in multiple families is NOTCH1, but o6% of BAV cases are accounted for by NOTCH1 variation13. It is clear that BAV is not a simple Mendelian trait9,15, but is indeed heritable, and, therefore, we applied genetic association methods typically used for complex traits.

With a goal of identifying genetic variants associated with BAV, leading to biological insight of the underlying causes, here we perform an unbiased genome scan in a large study of BAV cases (n 466) and controls (n 4,660), with replication in

additional samples of up to 1,326 cases and 8,103 controls. We identify two genetic variants that reached or were near genome-wide signicance levels (Po5 10 8). One is a low-

frequency intergenic variant rs6601627 (odds ratio (OR) 2.38,

Pafter-replication 3 10 15) with a substantially higher frequency

in BAV cases (8.3%) than in controls (4.2%), and the other one is an independent association signal at a common protein-altering variant p.Ser377Gly (rs3729856) in GATA4, which encodes a cardiac-specic transcription factor that is 151 kilobases(kb) away from the rst variant (Pafter-replication 8.8 10 8). Induced

pluripotent stem cells (iPSCs) with GATA4 disrupted by CRISPR-Cas9 demonstrate impaired transition of endothelial into mesenchymal cells (EndoMT), a critical step in valve formation16.

ResultsGWAS of BAV. To discover the underlying genetic basis of BAV, we successfully genotyped 498,075 genetic variants with enrichment of protein-altering variants (43.8% of variants examined) for 466 BAV cases and 4,660 controls. Imputation from the Haplotype Reference Consortium (HRC) panel17 enabled examination of a total of 12,320,487 variants. Clinical characteristics for BAV cases are summarized in Supplementary Table 1. Following a genome-wide association scan (Supplementary Figs 1 and 2), we examined three variants in replication cohorts with a combined total of up to 1,326 additional cases and 8,103 controls.

Variants near GATA4. The strongest result from the genome-wide discovery for BAV was observed for a genotyped low-frequency variant, rs6601627, in an intergenic region of

chromosome 8 (rs6601627, minor allele frequency (MAF) 4.1%,

OR 1.9, P

3.0 10 15; Table 1, Fig. 1 and

Supplementary Fig. 2). Ninety-seven imputed variants in this region also reached genome-wide signicance (Po5 10 8).

The two nearest genes are not obvious functional candidates (CTSB and DEFB135); however, this variant is 151 kb from the 30 end of the GATA4 gene (Fig. 1).

We also observed a common missense variant p.Ser377Gly in GATA4 (rs3729856) that was also associated with BAV (Pdiscovery 3.2 10 4, MAF 14.5%; Table 1 and Fig. 1). We

selected GATA4 p.Ser377Gly for in silico replication because it is a protein-altering variant and because it was located in the local genomic region of the most signicant variant (Supplementary Fig. 3). The p.Ser377Gly variant reached near genome-wide signicance after including in silico replication data (OR 1.31,

P 8.8 10 8; Supplementary Fig. 4), and exceeded the typical

signicance level used for exome-wide studies of coding variation (typically Po2 10 7; ref. 18). This suggests that GATA4 may

be the functional gene at this GWAS locus; however, further experiments will be needed to demonstrate which gene(s) causes BAV. The two variants at 8p23.1 (rs6601627 and rs3729856) appear to be independent of each other, since they were not in linkage disequilibrium (LD r2 0.013) and reciprocal conditional

association analysis maintained nominal signicance for both (Pcond rs6601627 8.92 10 9, Pcond rs3729856 0.012). After

including in silico replication data, the reciprocal conditional association analysis still maintained nominal signicance (Pmeta

rs6601627 1.52 10 9, Pmeta rs3729856 8.17 10 3;

Supplementary Fig. 5). The non-additive association tests showed that both variants appear to have dominant effects on risk of BAV (Supplementary Table 2).

The ExAC database characterizes protein-altering variants in 60,706 multiethnic individuals with whole-exome sequences19. ExAC lists 96 missense variants in GATA4 (95 of them have MAFo1%), a decit compared to the 140 variants predicted based on gene size. In addition, 9.4 loss-of-function (LoF) variants are predicted and only 1 was observed (p.Lys365Ter) out of 60,706 individuals with deep exome sequences. The probability that the gene is intolerant to LoF, a measure of the relative importance of gene function, is high (loss intolerance probability (pLI) score 0.8 where 40.9 is considered extremely intolerant).

Moreover, the missense GATA4 variant rs3729856 is predicted as benign or tolerated by PolyPhen-2 (ref. 20) and SIFT21, and it has a combined annotation dependent depletion (CADD) score 9.418 (in top 11% of deleterious variants in the human genome)22. These results suggest the importance of the GATA4 genes function, although the missense variant rs3729856 itself may not be signicantly deleterious.

We hypothesized that the functional BAV gene at this CTSB/ GATA4 locus would demonstrate high expression in heart or vascular tissue. Using the GTEx portal23, we examined mRNA expression levels of all genes within the 200 kb surrounding the noncoding associated variant rs6601627 and found that GATA4 showed strong expression in the heart (atrial appendage and left ventricle) and coronary artery, and also in the ovary, testis, pancreas and liver (Supplementary Fig. 6)23. The other genes in the region (NEIL2, FDFT1 and CTSB) showed ubiquitous expression levels across all tissues (Supplementary Fig. 6). Examination of all GTEx association results did not identify any signicant expression quantitative trait locus with the noncoding rs6,601,627 (Po10 5; ref. 23). We propose that this noncoding variant, or a variant tagged by it, inuences GATA4 expression in a manner not detectable by GTEeither exerting an inuence on gene expression levels only in the developing fetal heart or with a relatively modest effect that was not detectable in the current GTEx sample size.

combined

2 NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481 ARTICLE

Table 1 | Genetic variants associated with BAV.

Variants Discovery Replication Combined

Chr:posrsid Protein

change

Freq

Case

Ctrl (%)

N Case/

Ctrl

OR P value N Case/

Ctrl

Freq

Case/

Ctrl (%)

OR P value N Case/

Ctrl

OR P value

15

8:11778803rs6601627 Intergenic 8.2/3.7 466/

4,660

2.38 (1.813.13) 1.5 10

10 1,021/

5,357

7.2/4.2 1.73 (1.422.12) 1.1 10

7 1,487/

10,017

1.93 (1.642.27) 3.0 10

8

8:11614575rs3729856 p.S377G

GATA4

18.2/14.1 466/

4,660

1.39 (1.171.66) 3.2 10

4 1,326/

8,103

15.3/12.7 1.28 (1.141.45) 5.3 10

5 1,792/

12,763

1.31 (1.191.45) 8.8 10

8

16:72146374rs137867582 p.T1221M

DHX38

0.9/0.1 466/

4,660

13.14 (5.3932.04) 1.5 10

8 720/

5,831

0.37/0.15 2.87 (18.22) 5.0 10

2 1,186/

10,491

7.13 (3.6314) 1.2 10

BAV, bicuspid aortic valve; Ctrl, control; Freq, frequency; OR, odds ratio.

After detecting the association with both coding and noncoding variants at the GATA4 locus, we sought to examine the role of GATA4 in the development of the aortic valve. In the primitive heart tube, heart valves develop from endocardial cushions, which are formed by mesenchymal cells derived from endothelial cells (ECs) through a process called EndoMT16. Despite the critical role in heart valve development, the mechanism of EndoMT is not well understood2426. GATA4 was previously shown to be essential for heart formation and for endocardial cushion development in mice2729. Here we evaluated the impact of disruption of GATA4 on human iPSC differentiation into mesenchymal cells through EndoMT to examine the role of GATA4 in the development of aortic valves in humans.

The GATA4 knockout mouse is embryonically lethal between embryonic day (E) 7.0 and E9.5, and lacks a primitive heart tube28. Deletions of GATA4 in humans have been associated with congenital heart defects (CHDs)30,31, and a missense variantp.Gly296Ser was identied in a family with atrial and ventricular septal defects32. A mouse model of the p.Gly296Ser missense change is also embryonically lethal by E11.5, but a subset of these mice demonstrate semilunar valve stenosis and small defects of the atrial septum, thought to be resultant from defects in cardiomyocyte proliferation during embryogenesis33. Previous studies observed missense variants in GATA4 in patients with septal defects34, CHDs35 and Tetralogy of Fallot (ToF)35, but have not been tested in casecontrol models. The frequency of GATA4 variants in healthy controls is not clear from these studies and their pathogenicity is unknown. The co-appearance of congenital heart disease and testicular anomalies was found in a family with a GATA4 p.Gly221Arg mutation, thought to disrupt interaction with FOG2 and/or NR5A1, important factors for gonadal development36. GATA5 has 46% homology with GATA4. GATA5 sequence variants have been identied in humans with BAV37,38, and GATA5 knockout mice and zebrash demonstrate high rates of cardiac abnormalities39.

Chromatin conformation at the GATA4 locus. We attempted to evaluate the hypothesis that noncoding variants in LD with rs6601627 have an impact on expression of GATA4 during a critical stage of development. This is supported by prior evidence that GATA4 dosage has an impact on cardiac formation40. We rst identied potentially functional variants using RegulomeDB and HaploReg in the region near rs6601627 or variants in high LD (r240.6; ranging from rs112197605 to rs117851931;

hg19:chr8:1177495211838697)41,42. After examining local chromatin states, DNase-hypersensitive regions and transcription factor-binding sites, we identied rs118065347 as a variant likely to be functional because it co-localizes with binding regions for multiple transcription factors (including KAP1, CCNT2, CJUN, C-MYC, GATA2, HDAC2, HMGN3, JUND, MAX, SP1, TAL1, YY1 and ZBTB7A) and is in a known

enhancer active in fetal heart, left and right ventricle, right atrium, as well as other tissues42,43. This variant disrupts the binding motif for a variety of transcription factors including PAX6 (ref. 44). There are other candidate functional variants in high LD with the index variant, and molecular experiments will be required to denitively identify the functional variant(s) and the mechanism of action on aortic valve development.

We next asked which genes in the locus may interact with this candidate enhancer region. We identied chromatin interaction loops in K562 and GM12878 cells using chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) and high-throughput sequencing (Hi-C) data (Fig. 2)45,46. rs118065347 falls near the edge of a topologically associated domain spanning from hg19:chr8:11250000 to 11825000 dened by Hi-C in both cell lines. The variant falls inside a ChIA-PET loop connecting to a region also annotated as an enhancer 30 of GATA4 and C8orf49 and 50 of NEIL2. These data indicate that this distal region is brought in close proximity to GATA4 and disruption of this region may have direct impact on GATA4 expression. Further molecular experiments will be needed to clarify the gene(s) that has an impact on BAV.

Phenotypic characteristics of BAV cases in discovery sample. Among our 466 non-syndromic BAV cases, 93 (20%) reported one or more family members also having BAV (Supplementary Table 1). This suggests a high recurrence risk and supports the hypothesis of large-effect variants, but not necessarily Mendelian inheritance47. Majority of BAV cases were recruited from a cardiac surgery clinic at the University of Michigan Frankel Cardiovascular Center (FCVC), where patients are referred to cardiac surgery for aneurysm repair or valve replacement; thus, we found a high proportion of patients with thoracic aortic aneurysm (TAA; 83%). However, at these two loci, we saw no evidence for heterogeneity between BAV cases with or without TAA (Supplementary Table 3) and between BAV cases with or without a positive family history of BAV and/or TAA (Supplementary Table 4), suggesting that BAV probably has an impact on the risk of TAA because of altered haemodynamic blood ow and aortopathy from different mechanisms instead of sharing molecular mechanisms with TAA that has an impact on both aorta and valve tissue47. In addition, we did not nd evidence for heterogeneity in the association results at GATA4 and BAV subtypes (Supplementary Table 5) and among male and female subjects (Supplementary Table 6).

Implication of rs6601627 and rs3729856 in other CHDs. To investigate whether the two variants at GATA4 that we report are involved in development of other and more severe CHDs, we tested for association with 806 cases of ToF along with 5,029 matched controls and performed association tests for the two variants as described previously48. In an additive genetic model, we did not nd

NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481

100

Recombination rate (cM/Mb)Recombination rate (cM/Mb)

10

80

8

log 10(Pvalue)

log 10(Pvalue)

60

6

4

40

2

20

0

0

100

10

80

8

60

6

4

40

2

20

0

0

11.4 11.6 11.8 12 12.2 12.4 Position on chr8 (Mb)

Figure 1 | Regional association plot of the chr8 association region near GATA4in the discovery cohort. Genome-wide single variant association tests were performed on 466 BAV cases and 4,660 controls. The upper panel shows all variants that were directly genotyped in the chip array in this region. A missense variant (rs3729856, p.S377G) within GATA4 was observed to be associated with BAV with P 3.2 10 4, that reached P 8.8 10 8

following replication in 1,326 BAV cases and 8,103 controls. The bottom panel shows results after genotypes imputed from the HRC reference17. Coding

variants are represented by triangles and noncoding variants are represented by squares.

evidence for association between the noncoding rs6601627 and ToF (MAF 0.03, OR 0.89, 95% condence interval (CI) 0.671.20,

P 0.46); however, for GATA4 p.Ser377Gly, the association was

nominally signicant (MAF 0.11, OR 1.24, 95% CI 1.061.45.

P 0.007). This suggests that the regulatory variant associated

with BAV may act in a highly tissue or developmentally controlled manner to cause only BAV and not other CHDs, whereas GATA4 missense changes may have a more broad impact on other CHDs.

Missense variant in DHX38. The GWAS highlighted a rare missense variant (0.14% frequency in controls) in DHX38(p.Thr1221Met) with a strong association with BAV (OR 13.14,

95% CI 5.3932.04, P 1.5 10 8) in the discovery sample

(Supplementary Fig. 2). After genotyping of this variant in 720 cases and 5,831 controls, only 22 copies of the rare allele were identied (5 in cases and 17 in controls), providing a replication P value of 0.05. Additional large studies will be needed to conrm this rare variant association with BAV.

GATA4 deciency impairs EndoMT in iPSC-derived cells. We investigated the biological impact of GATA4 in the EndoMT process required for human valve formation. Human iPSCs were generated from peripheral blood mononuclear cells of a donor with normal trileaet aortic valve, using non-integrated DNA vectors containing OCT4, SOX2, C-MYC and KLF4 (ref. 49). The pluripotency of iPSCs was conrmed by expression of OCT4, SOX2, NANOG and SSEA4, TRA-1-60 and TRA-1-81 (Supplementary Fig. 7a,b). In addition, iPSCs generated teratoma-containing tissues from three germ layers, demonstrating their pluripotency in vivo (Supplementary Fig. 7c). In a previous study, wild-type GATA4 localized completely in the nucleus, whereas GATA4 mutant p.Ser377Gly (a C-terminal mutant) was shown to be partially distributed to the cytoplasm, indicating a LoF mutation50. To evaluate whether disruption of GATA4 may result in a LoF phenotype, iPSCs were electrotransfected with plasmid containing Cas9, GATA4 single guide RNA (sgRNA) and green uorescent protein (GFP) as an indicator for transfection51. As control, iPSCs were transfected with plasmid containing Cas9 and GFP. Transfected cells were enriched by ow cytometry

4 NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481 ARTICLE

11100000

11750000

11300000

11700000

11500000

Figure 2 | Chromatin interactions between associated variants and GATA4. The topological domain region containing associated variants (orange vertical lines), genes (green bars), chromatin interactions by Hi-C (blue loops) and ChIA-PET (purple loops), and chromatin state (outer ring and standard colours from ref. 68 but of signicance here: yellow as enhancers, red as promoters, green as transcribed, blue as CTCF and grey as inactive). All data are from K562 cells. rs3729856 is indicated as falling within a coding exon of GATA4. rs6601627 was identied as the associated variant to BAV and rs118065347 is the putative functional variant in linkage. rs11865347 overlaps an annotated enhancer as well as a ChIA-PET loop connecting to a region 30 of GATA4.

sorting based on GFP positivity (Supplementary Fig. 8ac). iPSCs were differentiated into ECs with efciency above 90% (Supplementary Fig. 8d).

The GATA4 level was signicantly lower in ECs from the GATA4 sgRNA-transfected group than control (Fig. 3a), indicating successful targeting to GATA4. When EndoMT was induced by TGFb2 and BMP2 in ECs, smooth muscle actin (SMA), a mesenchymal marker gene was upregulated in control cells (Fig. 3b). Noticeably, the GATA4 sgRNA group showed signicantly lower SMA levels (Fig. 3b). ECs were also explanted to collagen gel to induce EndoMT27. The GATA4 sgRNA group showed signicantly fewer mesenchymal cells migrating out after

3 days than control cells (Fig. 3c). Cells undergoing EndoMT express SMA and CD31 simultaneously at a certain point27. Immunouoresence staining of SMA and CD31, markers of EndoMT27, also showed signicantly less SMA and CD31 double-positive cells in the GATA4 sgRNA group (Fig. 3d). These results indicate that EndoMT was impaired by disruption of GATA4 with GATA4 sgRNA.

DiscussionIn this study we nd variants associated with BAV that reach genome-wide signicance. We identied association with a low-

NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481

frequency noncoding variant 151 kb from GATA4, as well as a common missense variant in GATA4. Although we cannot yet conrm the mechanism of action of the noncoding variant(s) on chromosome 8 on aortic valve development, chromatin conformation experiments suggest that the region near the associated variants appears to loop and physically interact with regions intronic to GATA4. This hypothesis could be tested in future functional experiments to investigate whether the noncoding BAV-associated variants identied here affect expression of GATA4 at a critical time in heart development. This could possibly disrupt EndoMT, a process important for normal trileaet aortic valve formation.

GATA4, a zinc-nger transcription factor, is one of three major transcription factors, together with Nkx2.5 and TBX5, that are critical for heart differentiation52. Although not previously associated with BAV, the GATA4 gene is a plausible biological candidate. The missense GATA4 mutation G296S disrupts the transcriptional cooperativity between GATA4 and TBX5, resulting in abnormal cellular functions related to morphogenetic defects53. Many mutations in GATA4 have been previously reported to be found in different kinds of CHDs atrial septal defect3234,5458, ventricular septal defect32,34,57,59,60 and ToF34,57,61, although mostly tested in family studies. Furthermore, the GATA4 mutations have been identied in CHD patients with various ancestries: European3234,58, Asian32,54,55,57,5961 and Native and Hispanic American34. In addition, GATA4 knockout mice are

embryonically lethal with heart defects. Mice that are missing GATA5 also develop BAV39, and rare GATA5 mutations have been identied in humans with BAV37.

We have provided evidence for the complexity of the BAV phenotype, with multiple genetic variants of incomplete penetrance contributing to susceptibility. To assess whether the two variants that we report are specic for BAV or whether they are also implicated in other CHDs, we studied cases of ToF, characterized by several cardiac malformations including an over-riding aorta, pulmonic stenosis, ventricular septal defect and right ventricular hypertrophy. We found that the common coding variant in GATA4 (p.Ser377Gly) was associated with an increased risk of ToF, whereas the low-frequency noncoding variant (rs6601627) was not associated. We speculate that the low-frequency noncoding variant disrupts a regulatory element that plays a critical role in regulating GATA4 expression in a precise time of cardiac embryogenesis that may have an impact on the valve more specically, whereas the common GATA4 missense variant might disrupt GATA4 function more generally and increase the risk of several cardiac malformations, including ToF. The frequencies of the associated variants at the GATA4 locus (variants with r240.6 in EUR samples of 1000 G) vary among different populations62. For example, among the noncoding variant rs6601627 and its 115 correlated variants, 108 variants have MAFo0.01 and 73 are monomorphic in East Asians62.

a

b

ECs

+ +

EndoMT ECs

Control

sgRNA

GATA4

GATA4 sgRNA

GATA4

50

SMA

50

GAPDH

37

GAPDH

37

MW (kDa) marker

**

Marker MW (kDa)

*

GATA4/GAPDH

1.5

1.0

0.5

0.0

SMA/GAPDH

1.5

1.0

0.5

0.0

*

+ +

EndoMT ECs

Control

Control

GATA4 sgRNA

GATA4 sgRNA

GATA4

sgRNA ECs

c

d

SMA/CD31/DAPI

Number of mesenchymal cells

relative to control

1.5

1.0

0.5

**

0.5

Collagen gel assay

Figure 3 | EndoMT is a key process in aortic valve development and is impaired by GATA4 deciency. (a) Western blot of GATA4 and GAPDH from control and GATA4 sgRNA ECs. GATA4 sgRNA ECs were differentiated from iPSCs transfected with px458 with GATA4 sgRNA and enriched by GFP. Control ECs were derived from iPSCs with px458 and enriched by GFP. An uncropped version is presented in Supplementary Fig. 9. Lower panel: quantication of western blot data. The data were normalized to control ECs. Experiments were repeated three times; averages and standard derivations were plotted.(b) Western blot of SMA and GAPDH from control ECs, control ECs undergoing EndoMT, GATA4 sgRNA ECs and GATA4 sgRNA ECs undergoing EndoMT. An uncropped version is presented in Supplementary Fig. 10. Lower panel: quantication of western blot data. The data were normalized to control ECs undergoing EndoMT. Experiments were repeated three times; averages and standard derivations were plotted. (c) Numbers of mesenchymal cells from control and GATA4 sgRNA in collagen gel assay. The data were normalized to control. Experiments were repeated three times; averages and standard derivations were plotted. (d) Immunouorescence staining of SMA and CD31 of the control and GATA4 sgRNA undergoing EndoMT. Scale bars, 50 mm. EC, endothelial cell;

EndoMT, endothelial-to-mesenchymal transition; iPSC, induced pluripotent stem cell; kDa, kilodalton; MW, molecular weight. *Po0.05; **Po0.01.

Control

GATA4 sgRNA

6 NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481 ARTICLE

Association studies in other populations will be critical for determining whether the association exists in other populations and may be helpful at narrowing the associated interval.

To investigate the possible role of GATA4 in aortic valve development, we used sgRNA-guided Cas9 to disrupt GATA4 in iPSCs from a healthy human donor with normal tricuspid aortic valves (TAVs). The iPSCs were differentiated into ECs and then induced to mesenchymal cells through EndoMT. We demonstrated that deciency of GATA4 impaired the EndoMTs, a critical step in valve formation16 (Fig. 3). This indicates that GATA4 is required for aortic valve formation and that disruption of the GATA4 gene, either by noncoding or protein-altering variants, may affect aortic valve formation.

Methods

GWAS genotyping and genotype imputation. We performed genotyping of a combined set of 498,075 GWAS variants, including 217,957 protein-altering variants, using a GWAS exome chip array (Illumina Human CoreExome). To avoid

any potential batch effects, cases and controls were genotyped using the same array in the same genotyping centre (Sequencing and Genotyping Core at the University of Michigan). Genotype calling was performed using GenTrain version 2.0 in GenomeStudio V2011.1 (Illumina) using identical cluster les for cases and controls. Samples with o98% genotype calls, evidence of gender discrepancy, duplicates as well as individuals with non-European ancestry identied by plotting the rst 10 genotype-driven principal components were excluded from further analysis. We performed variant-level quality control (QC) by excluding 22,983 variants that met any of the following criteria: variants with a cluster separation scoreo0.3, o98% genotype call rate or deviation from HardyWeinberg equilibrium(Po1 10 5). We phased the autosomal genotype data using SHAPEIT2 (ref. 63)

and imputed variants from the HRC v1 reference panel17 using minimac3 (ref. 64). We excluded poorly imputed variants with imputation r2o0.3, and then merged the genotyped variants and the successfully imputed variants to a combined data set, which contains 12,320,487 variants in total.

Description of cases in discovery cohort. We collected DNA from consented individuals with BAV from the FCVC at the University of Michigan as part of the University of Michigan BAV registry or the Cardiovascular Health Improvement Project (CHIP). All repository projects utilized for this study are approved by the University of Michigan, Medical School, Institutional Review Board (IRBMED), and informed consent was obtained from study participants. Patients were typically seen in clinic for aortic valve replacement or aortic aneurysm. Diagnoses of BAV were made by cardiac surgeons upon visual inspection of the aortic valve during open surgery for aneurysm repair or valve replacement. BAV cases with major syndromic connective-tissue disorders (for example, Marfan syndrome) were excluded. DNA was isolated from peripheral blood lymphocytes.

Description and selection of controls in discovery cohort. We identied potential controls from a surgical-based biobank, the Michigan Genomics Initiative (MGI), that were genotyped with the same GWAS array (Illumina Human CoreExome). After excluding those with possible aortic disease (n 1,586,

Supplementary Table 7), we were left with 15,642 potential controls with GWAS data. We performed age-matching by requiring controls to have a birth year within 5

and 10 years of the case. From the available controls in the appropriate age and sex

category, we selected the best ethnic match for each case and repeated the greedy algorithm until a control was selected for each case. We repeated the entire process so that 10 controls were selected for each case. We opted for this approach to provide the best ancestry-matching between cases and controls, to reduce the potential for false-positives due to ethnicity mismatch and to also provide the most power for rare variants that increase the risk of BAV by including the highest number of matching controls. All MGI research subjects provided informed consent.

Statistical analyses. In the discovery cohort, we performed association testing for the BAV status using logistic regression with single genetic variants (295,759 with MAF41%), with age, sex and the rst four principal components as covariates using PLINK for hard call genotypes65 and the EPACTS software (http://csg.sph.umich.edu//kang/epacts/

Web End =http://csg.sph.umich.edu//kang/epacts/) for imputed dosages. We identied two genetic variants that were directly genotyped and reached genome-wide signicance levels (Po5 10 8). We observed no evidence for ination because

of population stratication (l 1.033, Supplementary Fig. 1). We observed a

genotyped missense variant within 200 kb (rs3729856) of one of the signicant noncoding variants (rs6601627) and selected this third variant for follow-up in additional samples.

Association with ToF. A total of 835 unrelated ToF cases and 5,159 controls were genotyped and imputed from 1,000 Genomes Phase 3 for the region 1112 MB on

chromosome 8 using IMPUTE2 (refs 48,66). The association tests were performed using logistic regression of the best-guess genotypes for all imputed single-nucleotide polymorphisms (SNPs) with IMPUTE2 info score Z0.5 and with MAF Z0.01 in controls using SNPTEST67. This study has been approved by Newcastle and North Tyneside NHS Research Ethics Committee.

Incorporating gene expression and epigenetics data. We assessed expression levels of relevant genes using the GTEx server (http://www.gtexportal.org/home/

Web End =http://www.gtexportal.org/home/ )23. We obtained the Hi-C interaction calls from Rao et al.45,46 and ChIA-PET interactions from Phanstiel et al.45,46 available through the ENCODE DCC accessions ENCSR000FDB and ENCSR752QCX (https://www.encodeproject.org/

Web End =https://www.encodeproject.org/). ChromHMM data are displayed from the K562 Genome Segmentation by ChromHMM from ENCODE/Analysis available at https://genome.ucsc.edu/

Web End =https://genome.ucsc.edu/ and is created through chromatin segmentation using eight histone modications, CTCF, Pol2 and open chromatin annotations68.

IPSC generation and culture. The procedure of iPSC derivation was performed according to methods we described49: Peripheral blood mononuclear cells (PBMC) were separated from human peripheral blood with lymphocyte seperation medium (LSM) (MP Biomedicals LLC.), cultured in medium containing Iscoves modied Dulbeccos medium (IMDM) (Life Technologies Corp.), 10% fetal bovine serum (Life technologies Corp.), thrombopoietin (TPO), stem-cell factor (SCF) and FMS-like tyrosine kinase-3 (FLT-3) at a nal concentration of 100 ng ml 1, granulocytemacrophage colony-stimulating factor and interleukin-3 at a nal concentration of 10ng ml 1 (Peprotech Inc.) and penicillinstreptomycin(Life Technologies Corp.), and electrotransfected with episomal DNA plasmids containing OCT4, SOX2, KLF4 and C-MYC using Nucleofector 2 Device(Lonza Corp.). At around day 30 post infection, the colonies became compact. The colonies were mechanically picked up from the culture dishes and rst cultured with mouse embryonic broblasts for three passages69 and transited to TesRE8 medium (Stemcells Inc.) on matrigel-coated (BD Corp.) dishes. iPSCs were passaged every 46 days with Versene (Life Technologies Corp.). In addition, iPSCs from passages 2535 were used in experiments.

Teratoma formation in immune-decient mice. Conduction of Animal experiments was in compliance with the regulations of the Unit for Laboratory Animal Medicine at the University of Michigan. Two million iPSCs were injected subcutaneously into each ank of the recipient male, 68-week-old nonobese diabetic-severe combined immunodecient mice (Jackson Laboratory, Bar Harbor, Maine). Three to ve weeks after injection, teratomas were collected from the mouse anks and xed with formalin (Thermo Corp.) for 2 days. The tumours were then embedded in parafn, and sections were prepared with a microtome (Leica Corp.) and stained with haematoxylin and eosin staining solutions from Thermo Corp. The slides were examined and micrographs were taken under brighteld with microscope (Nikon Corp.).

GATA4 sgRNA design and electrotransfection of iPSCs. sgRNAs were designed to target GATA4 exon2 (the rst coding exon) with sgRNA design tool (http://www.genome-engineering.org

Web End =http:// http://www.genome-engineering.org

Web End =www.genome-engineering.org ) developed by Zhang and co-workers51. Sequence of GATA4 sgRNA was 50-CGCGCCGTGCATGAAGGCGCCGG-30. Target site was chr8:-11565888. Quality score was 93. Minimal number of mismatch nucleotides in offsite targets was 3. SgRNAs were cloned into PX458, which contains SpCas9-2AEGFP, using AgeI and EcoRI at 50 and 30 cloning sites51. One million iPSCs were electrotransfected with constructed 5 mg PX458 containing GATA4 sgRNA, using the Lonza Human Stem Cell Nucleofector Kit 2 with programme U-023 on Nuclefector 2 device (Lonza Ltd.). Another one million iPSCs were electrotransfected with the PX458 vector as control under the same conditions.

EC differentiation from iPSCs. To differentiation iPSCs into ECs, iPSCs were dissociated with Versene (Life Technologies Corp.) into single cells and seeded at2 104 cells per cm2 with the TesRE8 (Stemcell Technology Inc.) medium supple

mented with Rocki (Y27632, Stemgent Inc.). When the cells reached a conuence of 2030%, the medium was changed into a differentiation medium, which contained DMEM-F12 (Life Technologies Corp.), B27 supplement without vitamin A (Life Technologies Corp.), L-glutamine (Life Technologies Corp.), penicillinstreptomycin (Life Technologies Corp.), 400 mM 1-thioglycerol (Sigma Corp.), 50 mg ml 1 ascorbic acid (Sigma Corp.), 25 ng ml 1 BMP4 (R&D Systems Corp.) and 6 mM GSK3 inhibitor CHIR99021 (Sigma Corp.). Differentiation medium was refreshed daily for 3 days. Then, cells were dissociated with Accutase (Life Technologies Corp.) and seeded at 1 104 cells per cm2 on matrigel (BD Corp.)-coated dishes with an EC medium

containing Stempro34(Life Technologies Corp.), Stempro34 supplement (Life Technologies Corp.), L-glutamine (Life Technologies Corp.), penicillinstreptomycin (Life Technologies Corp.) and 50 ng ml 1 vascular endothelial growth factor (Peprotech

Inc.). Medium was refreshed every 2 days for 13 days.

Immunouorescence staining and ow cytometry. Immunouorescence staining and ow cytometry were performed as follows: rst, cells were xed in 4%

NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481

formaldehyde (Thermo Corp.) for 1 h at room temperature, and then the cells were washed with DPBS (Thermo Corp.) once and incubated with primary antibodies for 2 h at room temperature70. The following primary antibodies were used: anti-OCT4 (mouse IgG, dilute 500 times upon usage, sc-5,279, Santa Cruz Biotechnology Inc.), anti-SOX2 (mouse IgG, dilute 500 times upon usage, sc-365964, Santa Cruz Biotechnology Inc.), anti-NANOG (rabbit polyclonal, dilute 500 times upon usage, REC-RCAB004PF, Cosmo Inc.), anti-SSEA4 (mouse IgG, dilute 100 times upon usage, 60062, Stemcell Technology Inc.), anti-TRA-1-60 (mouse IgM, dilute 100 times upon usage, 60,064, Stemcell Technology Inc.), anti-TRA-1-81 (mouse IgM, dilute 100 times upon usage, 60,065, Stemcell Technology Inc.), anti-CD31 (rabbit polyclonal, dilute 500 times upon usage, ab28364, Abcam Inc.) and anti-SMA (mouse IgG, dilute 1,000 times upon usage, A5228, Sigma Corp.). Cells were washed three times with DPBS (Thermo Corp.), and then incubated with secondary antibodies for 1 h at room temperature. The following uorochrome-conjugated secondary antibodies were used: Alexa Fluor 488 goat anti-rabbit IgG (goat, dilute 1,000 times upon usage, A11034, Thermo Corp.), Alexa Fluor 488 goat anti-mouse IgG (goat, dilute 1,000 times upon usage, A32723, Thermo Corp.) and Alexa Fluor 594 goat anti-mouse IgG (goat, dilute 1,000 times upon usage, A11032, Thermo Corp.). Slides were mounted with anti-fade mounting media containing 4,6-diamidino-2-phenylindole (Prolong gold, Life Technologies Corp.), and were observed on a Nikon A1 confocal microscope (Nikon Corp.). In a ow cytometry study, electrotransfected iPSCs were dissociated into single cells with Accutase (Stemcell Technology Inc.), and applied to the MoFlo Astrios (Beckman Coulter Inc.) ow cytometry machine.

Western blot analysis. Whole-cell extracts were prepared using RIPA buffer (1% NP-40, 1% sodium deoxycholate, 0.1% SDS, 0.15 M NaCl, 0.01 M sodium phosphate, 2 mM EDTA, 50 mM sodium uoride, 0.2 mM Na3VO4.2H2O, 100 U ml 1 protease inhibitor), resolved on SDSPAGE gels and transferred to acetate cellulose membranes. Primary antibodies used were anti-GATA4 (rabbit IgG, diluted500 times upon usage, 36,968, Cell Signaling Technology Inc.), anti-SMA (mouse IgG, diluted 300 times upon usage, A5228 Sigma) and anti-GAPDH (rabbit IgG, diluted 2,000 times upon usage, sc25778, Santa Cruz Inc.). Secondary antibodies used were IRDye800CW Donkey anti-Mouse (92532212), IRDye680LT Donkey anti-Rabbit (92568023), IRDye800CW Donkey anti-Rabbit (92532213;all secondary antibodies were diluted 500 times upon usage and purchased from Licor Inc.). The Licor western blot detection system was used for the dual-colour imaging. Uncropped versions of western blots are presented in Supplementary Figs 9 and 10. ImageJ was used for quantication of bands. Each band was normalized by GAPDH. Experiments were repeated three times. Average value and s.d. were plotted.

Endothelial-to-mesenchymal transition and collagen gel assay. EndoMT was induced by changing medium to an EndoMT-inducing medium that contained stempro34 medium with stempro34 supplement (Life Science Technology Corp.), L-glutamine (Life Technologies Corp.), penicillinstreptomycin (Life Technologies Corp.), 200 ng ml 1 BMP2 (Peprotech Inc.) and 50 ng ml 1 TGFb2 (Peprotech

Inc.). Non-EndoMT control groups were kept in EC medium containing Stem-pro34 (Life Technologies Corp.), Stempro34 supplement (Life Technologies Corp.), L-glutamine (Life Technologies Corp.), penicillinstreptomycin (Life Technologies Corp.) and 50 ng ml 1 vascular endothelial growth factor (Peprotech Inc.). Cells were collected 3 days after induction.

Type I collagen (Sigma Corp.) at 1 mg ml 1 (nal concentration) was mixed with stempro34 medium, stempro34 supplement (Life Science Technology Corp.)

and 50 mM NaOH (Sigma Corp.). The mixture was poured into 24-well tissue culture plates (0.5 ml per well) and allowed to gel in 5% CO2 incubator at 37 C for 30 min. And then 0.5 ml EndoMT-inducing medium was added. After 3 days, pictures of cells were taken with 100 times magnicance under the Eclipse Ti-U inverted research microscope (Nikon Corp.). Mesenchymal cells that migrated out in three pictures from different elds were counted. Experiments were repeated three times. Average value and standard derivation were plotted.

Genetic association replication cohorts. CHIP. In the CHIP replication cohort, an additional 140 BAV cases from the University of Michigan FCVC biobank were collected. These samples were genotyped using the same GWAS array as the discovery cohort, but were only examined for the three variants described here. The association was tested using PLINK65 with 1,400 age-, sex- and ancestry-matched controls from the MGI study, which were independent samples from previously used controls. Informed consent was obtained from all participants and approval was obtained from the Institutional Review Board of the University of Michigan Medical School.

Montreal Heart Institute. In the Montreal Heart Institute (MHI) biobank, 305 BAV cases and 2,746 controls were collected and genotyped on the Illumina Core Exome array at the MHI Pharmacogenomic Centre. Controls were selected by excluding those with myocardial infarction (MI), percutaneous coronary intervention (PCI), Angina, congestive heart failure (CHF), valve defects, heart surgeries, heart arrest, atrial brillation and sudden cardiac death. Genotyping was performed with the Illumina HumanExome array. Association analysis was performed in PLINK65 using a logistic regression model correcting for sex, age and

principal components of ancestry 110. The project has been approved by the Ethics Committee of the MHI and informed consent was obtained from study participants.

Partners HealthCare. In the Partners HealthCare cohort, 452 Caucasian BAV cases were identied from the electronic medical records (EMRs) of Partners HealthCare (Boston, MA). Individual echocardiographic images were reviewed to conrm BAV diagnosis. Whole-blood DNA was genotyped using the Illumina Omni2.5 Beadchip. The Framingham Heart Study dbGaP cohort, genotyped using the Illumina Omni5.0 Beadchip, was used as controls. QC and population stratication of the genotype data were performed in PLINK65. SNPs with MAF less than 1%, without physical map reference, not in HardyWeinberg equilibrium (Po10 4), differential missingness (Po10 5), were removed. Related individuals (PI_HAT40.25) were excluded. Genome-wide IBD and IBS were used to detect outliers and clusters. After merging case and control genotypes, additional genotypes have been imputed against the 1,000 Genome reference (phase3) and HRC (Michigan University) panels using SHAPEIT2 (ref. 63) and IMPUTE2 (ref. 66). After QC, 452 cases and 1,634 controls (1,094 males 992 females) with

7.5 million markers were analysed using an additive logistic regression model accounting for gender, age and principal components. This study has been approved by Partners HealthCare Human Research Committee, and informed consent was obtained from study participants.

University of Texas Health Science Centre. In the University of Texas Health Science Centre cohort, 765 patients with sporadic TAAs or aortic dissections were collected and genotyped. In all, 874 genotypes from dbGAP (NINDS Neurologically Normal control collection) were used as controls. QC and population stratication of the genotype data were performed in PLINK65. SNPs with MAF less than 1% or missing more than 1% of genotypes were excluded. Multidimensional scaling was used to detect and exclude population outliers.

We imputed additional genotypes against 1,000 Genomes Phase3 using SHAPEIT2 (ref. 63) and IMPUTE2 (ref. 66). After QC, a total of 152 BAV cases and 633 TAV cases or 874 controls were analysed using an additive logistic regression model accounting for gender and principal components. This study has been approved by the Committee for the Protection of Human Subjects atUT Health Science Center at Houston, and informed consent was obtained from study participants.

ASAPARTISTPOLCAOlivia cohort. Three cohorts ASAP, ARTIST and POLCA were included in this replication group with a total of 275 BAV cases and 1,686 controls used for analysis. The POLCA/Olivia cohort is a merged cohort with a total of 1,295 individuals. The ASAP cohort consists of 429 patients genotyped on Illumina 610wQuad beadchips. Approximately 588,400 SNPs were provided after QC. The Artist cohort consists of 406 samples genotyped with Omni2.5 Quad beadchips on 2,443,180 SNPs. In POLCA, 625 control samples were genotyped on Illumina 610kwQuad, and in Olivia 670 control samples were genotyped on Illumina 1M-genotyping arrays. The vast majority of included samples are of Scandinavian ancestry. For the ASAP database, where ancestry is specically registered, this corresponds to 495% of the individuals, supported by PCA plots of genotype clustering. Imputation was performed using Impute2 from 1,000G phase1 v3 (ref. 66). Analysis was performed using SNPTEST67, with age, sex and rst 10 principal components as covariates. This study was approved by the Regional Ethical Committee of Stockholm, and informed consent was obtained from study participants.

BioMe. The Mount Sinai BioMe Biobank (BioMe) is an ongoing, prospective, hospital- and outpatient-based population research programme operated by The Charles Bronfman Institute for Personalized Medicine at Mount Sinai and has enroled over 33,000 participants since September 2007. BioMe is an EMR-linked biobank that integrates research data and clinical care information for consented patients at The Mount Sinai Medical Center, which serves diverse local communities of upper Manhattan with broad health disparities. BioMe populations include 25% of African ancestry (AA), 36% of Hispanic Latino ancestry (HL), 30% of white European ancestry (EA) and 9% of other ancestry. The BioMe disease burden is reective of health disparities in the local communities. BioMe operations are fully integrated in clinical care processes, including direct recruitment from clinical sites waiting areas and phlebotomy stations by dedicated recruiters independent of clinical care providers, prior to or following a clinician standard of care visit. Recruitment currently occurs at a broad spectrum of over 30 clinical care sites. Information on BAV status, age and sex was derived from participants EMRs. BAV cases were dened as BioMe participants with the ICD-9 code 746.4 (Congenital insufciency of aortic valve). In total, there were 41 BAV cases with available genotyping data (8 AA and 13 HL BAV cases genotyped on the Innium Multi-Ethnic Global BeadChip from Illumina as well as 13 additional HL and 7 EA BAV cases genotyped on the Illumina HumanOmniExpressExome-8 v1.0 BeadChip. For each case, three controls were selected by genetically matching using the rst two genetic principal components and stratication by age and sex. Logistic regression was performed in PLINK for the three SNPs in the four groups65. We performed analyses both including and excluding the BioME non-European samples, and results were highly similar. We present results in this study excluding the non-European samples since there were few cases. This study has been approved by Icahn School of Medicine IRB and informed consent was obtained from study participants.

Data availability. The data that support the ndings of this study are available from the corresponding author upon reasonable request.

8 NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481 ARTICLE

References

1. Hoffman, J. I. & Kaplan, S. The incidence of congenital heart disease. J. Am. Coll. Cardiol. 39, 18901900 (2002).

2. Tutar, E., Ekici, F., Atalay, S. & Nacar, N. The prevalence of bicuspid aortic valve in newborns by echocardiographic screening. Am. Heart J. 150, 513515 (2005).

3. Losenno, K. L., Goodman, R. L. & Chu, M. W. Bicuspid aortic valve disease and ascending aortic aneurysms: gaps in knowledge. Cardiol. Res. Pract. 2012, 145202 (2012).

4. Ward, C. Clinical signicance of the bicuspid aortic valve. Heart 83, 8185 (2000).

5. Siu, S. C. & Silversides, C. K. Bicuspid aortic valve disease. J. Am. Coll. Cardiol. 55, 27892800 (2010).

6. Michelena, H. I. et al. Incidence of aortic complications in patients with bicuspid aortic valves. JAMA 306, 11041112 (2011).

7. Michelena, H. I. et al. Natural history of asymptomatic patients with normally functioning or minimally dysfunctional bicuspid aortic valve in the community. Circulation 117, 27762784 (2008).

8. Roberts, W. C. & Ko, J. M. Frequency by decades of unicuspid, bicuspid, and tricuspid aortic valves in adults having isolated aortic valve replacement for aortic stenosis, with or without associated aortic regurgitation. Circulation 111, 920925 (2005).

9. Ellison, J. W. et al. Evidence of genetic locus heterogeneity for familial bicuspid aortic valve. J. Surg. Res. 142, 2831 (2007).

10. Garg, V. Molecular genetics of aortic valve disease. Curr. Opin. Cardiol. 21, 180184 (2006).

11. Cripe, L., Andelnger, G., Martin, L. J., Shooner, K. & Benson, D. W. Bicuspid aortic valve is heritable. J. Am. Coll. Cardiol. 44, 138143 (2004).

12. Martin, L. J. et al. Evidence in favor of linkage to human chromosomal regions 18q, 5q and 13q for bicuspid aortic valve and associated cardiovascular malformations. Hum. Genet. 121, 275284 (2007).

13. Foffa, I. et al. Sequencing of NOTCH1, GATA5, TGFBR1 and TGFBR2 genes in familial cases of bicuspid aortic valve. BMC Med. Genet. 14, 44 (2013).14. Wooten, E. C. et al. Application of gene network analysis techniques identies AXIN1/PDIA2 and endoglin haplotypes associated with bicuspid aortic valve. PLoS ONE 5, e8830 (2010).

15. McBride, K. L. et al. NOTCH1 mutations in individuals with left ventricular outow tract malformations reduce ligand-induced signaling. Hum. Mol. Genet. 17, 28862893 (2008).

16. Lin, C. J., Lin, C. Y., Chen, C. H., Zhou, B. & Chang, C. P. Partitioning the heart: mechanisms of cardiac septation and valve development. Development 139, 32773299 (2012).

17. McCarthy, S. et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat. Genet. 48, 12791283 (2016).

18. Sveinbjornsson, G. et al. Rare mutations associating with serum creatinine and chronic kidney disease. Hum. Mol. Genet. 23, 69356943 (2014).

19. Lek, M. et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285291 (2016).

20. Adzhubei, I., Jordan, D. M. & Sunyaev, S. R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. Chapter 7, Unit 7.20 (2013).

21. Kumar, P., Henikoff, S. & Ng, P. C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 4, 10731081 (2009).

22. Kircher, M. et al. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46, 310315 (2014).

23. GTEx Consortium. The Genotype-Tissue Expression (GTEx) project. Nat. Genet. 45, 580585 (2013).

24. de Lange, F. J. et al. Lineage and morphogenetic analysis of the cardiac valves. Circ. Res. 95, 645654 (2004).

25. Lincoln, J., Aleri, C. M. & Yutzey, K. E. Development of heart valve leaets and supporting apparatus in chicken and mouse embryos. Dev. Dyn 230, 239250 (2004).

26. Wirrig, E. E. & Yutzey, K. E. Conserved transcriptional regulatory mechanisms in aortic valve development and disease. Arterioscler. Thromb. Vasc. Biol. 34, 737741 (2014).

27. Rivera-Feliciano, J. et al. Development of heart valves requires Gata4 expression in endothelial-derived cells. Development 133, 36073618 (2006).

28. Molkentin, J. D., Lin, Q., Duncan, S. A. & Olson, E. N. Requirement of the transcription factor GATA4 for heart tube formation and ventral morphogenesis. Genes Dev. 11, 10611072 (1997).

29. Kuo, C. T. et al. GATA4 transcription factor is required for ventral morphogenesis and heart tube formation. Genes Dev. 11, 10481060 (1997).

30. Pehlivan, T. et al. GATA4 haploinsufciency in patients with interstitial deletion of chromosome region 8p23.1 and congenital heart disease. Am. J. Med. Genet. 83, 201206 (1999).

31. Kennedy, S. J., Teebi, A. S., Adatia, I. & Teshima, I. Inherited duplication, dup(8) (p23.1p23.1) pat, in a father and daughter with congenital heart defects. Am.J. Med. Genet. 104, 7980 (2001).

32. Garg, V. et al. GATA4 mutations cause human congenital heart defects and reveal an interaction with TBX5. Nature 424, 443447 (2003).

33. Sarkozy, A. et al. Spectrum of atrial septal defects associated with mutations of NKX2.5 and GATA4 transcription factors. J. Med. Genet. 42, e16 (2005).

34. Tomita-Mitchell, A., Maslen, C. L., Morris, C. D., Garg, V. & Goldmuntz, E. GATA4 sequence variants in patients with congenital heart disease. J. Med. Genet. 44, 779783 (2007).

35. Zhang, W. et al. GATA4 mutations in 486 Chinese patients with congenital heart disease. Eur. J. Med. Genet. 51, 527535 (2008).

36. Lourenco, D. et al. Loss-of-function mutation in GATA4 causes anomalies of human testicular development. Proc. Natl Acad. Sci. USA 108, 15971602 (2011).

37. Bonachea, E. M. et al. Rare GATA5 sequence variants identied in individuals with bicuspid aortic valve. Pediatr. Res. 76, 211216 (2014).

38. Padang, R., Bagnall, R. D., Richmond, D. R., Bannon, P. G. & Semsarian, C. Rare non-synonymous variations in the transcriptional activation domains of GATA5 in bicuspid aortic valve disease. J. Mol. Cell Cardiol. 53, 277281 (2012).

39. Laforest, B., Andelnger, G. & Nemer, M. Loss of Gata5 in mice leads to bicuspid aortic valve. J. Clin. Invest. 121, 28762887 (2011).

40. Pu, W. T., Ishiwata, T., Juraszek, A. L., Ma, Q. & Izumo, S. GATA4 is a dosage-sensitive regulator of cardiac morphogenesis. Dev. Biol. 275, 235244 (2004).

41. Ward, L. D. & Kellis, M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 40, D930D934 (2012).

42. Boyle, A. P. et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 22, 17901797 (2012).

43. Bernstein, B. E. et al. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 5774 (2012).

44. Piper, J. et al. Wellington: a novel method for the accurate identication of digital genomic footprints from DNase-seq data. Nucleic Acids Res. 41, e201 (2013).

45. Phanstiel, D. H., Boyle, A. P., Heidari, N. & Snyder, M. P. Mango: a bias-correcting ChIA-PET analysis pipeline. Bioinformatics 31, 30923098 (2015).

46. Rao, S. S. et al. A 3D map of the human genome at kilobase resolution reveals principles of chromatin looping. Cell 159, 16651680 (2014).

47. Pasta, S. et al. Difference in hemodynamic and wall stress of ascending thoracic aortic aneurysms with bicuspid and tricuspid aortic valve. J. Biomech. 46, 17291738 (2013).

48. Cordell, H. J. et al. Genome-wide association study identies loci on 12q24 and 13q32 associated with tetralogy of Fallot. Hum. Mol. Genet. 22, 14731481 (2013).

49. Su, R. J. et al. Efcient generation of integration-free ips cells from human adult peripheral blood using BCL-XL together with Yamanaka factors. PLoS ONE 8, e64496 (2013).

50. Wang, E. et al. Identication of functional mutations in GATA4 in patients with congenital heart disease. PLoS ONE 8, e62138 (2013).

51. Ran, F. A. et al. Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. 8, 22812308 (2013).

52. Huang, W. Y., Cukerman, E. & Liew, C. C. Identication of a GATA motif in the cardiac alpha-myosin heavy-chain-encoding gene and isolation of a human GATA-4 cDNA. Gene 155, 219223 (1995).

53. Ang, Y. S. et al. Disease model of GATA4 mutation reveals transcription factor cooperativity in human cardiogenesis. Cell 167, 17341749.e1722 (2016).

54. Hirayama-Yamada, K. et al. Phenotypes with GATA4 or NKX2.5 mutations in familial atrial septal defect. Am. J. Med. Genet. A 135, 4752 (2005).

55. Yang, Y. Q. et al. Mutation spectrum of GATA4 associated with congenital atrial septal defects. Arch. Med. Sci. 9, 976983 (2013).

56. LaHaye, S. et al. Utilization of whole exome sequencing to identify causative mutations in familial congenital heart disease. Circ. Cardiovasc. Genet. 9, 320329 (2016).

57. Mattapally, S., Nizamuddin, S., Murthy, K. S., Thangaraj, K. & Banerjee, S. K.c.620C4T mutation in GATA4 is associated with congenital heart disease in South India. BMC Med. Genet. 16, 7 (2015).58. Posch, M. G. et al. Mutations in GATA4, NKX2.5, CRELD1, and BMP4 are infrequently found in patients with congenital cardiac septal defects. Am. J. Med. Genet. A 146a, 251253 (2008).

59. Yang, Y. Q. et al. Novel GATA4 mutations in patients with congenital ventricular septal defects. Med. Sci. Monit. 18, CR344CR350 (2012).

60. Wang, J. et al. A novel GATA4 mutation responsible for congenital ventricular septal defects. Int. J. Mol. Med. 28, 557564 (2011).

61. Yang, Y. Q. et al. GATA4 loss-of-function mutations underlie familial tetralogy of fallot. Hum. Mutat. 34, 16621671 (2013).

NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15481

62. Auton, A. et al. A global reference for human genetic variation. Nature 526, 6874 (2015).

63. Delaneau, O., Marchini, J. & Zagury, J. F. A linear complexity phasing method for thousands of genomes. Nat. Methods 9, 179181 (2012).

64. Fuchsberger, C., Abecasis, G. R. & Hinds, D. A. minimac2: faster genotype imputation. Bioinformatics 31, 782784 (2015).

65. Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet.Hum. Genet. 81, 559575 (2007).

66. Howie, B. N., Donnelly, P. & Marchini, J. A exible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 5, e1000529 (2009).

67. 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. 39, 906913 (2007).

68. Ernst, J. & Kellis, M. ChromHMM: automating chromatin-state discovery and characterization. Nat. Methods 9, 215216 (2012).

69. Jiao, J. et al. Modeling Dravet syndrome using induced pluripotent stem cells (iPSCs) and directly converted neurons. Hum. Mol. Genet. 22, 42414252 (2013).

70. Jiao, J. et al. Promoting reprogramming by FGF2 reveals that the extracellular matrix is a barrier for reprogramming broblasts to pluripotency. Stem Cells 31, 729740 (2013).

Acknowledgements

We thank all participants of the CHIP and BAV registry at the University of Michigan for their contribution to research. We appreciate the valuable efforts from the CHIP and BAV registry collection team (W.H., A.D., A.L., M.T., L.F., Jennifer McNamara, Joyvina Evans, Daniel Ferman, Neelima Goyal, Joanna Hider, Tom Jones, John Kruszewski, Brian Kulick, Jamie Love-Nichols, Joe Luciano, Alison Lunau, Phillip Nulty, Kevin Packard, Ian Patterson, Bailey Pearce, Mike Ranella, Alexander Shevchik, Carolyn Sommer, Christina White). The University of Michigan Health System CHIP was supported by the FCVC. We appreciate the Aikens Fund for Aortic Research (C.J.W. and B.Y.) and McKay research award for supporting this project (B.Y.). B.Y. is supported by American Association of Thoracic Surgery (AATS) Graham Foundation and Thoracic Surgery Foundation of Research and Education (TSFRE). C.J.W. is supported by HL109946, HL130705 and HL127564. S.K.G. is supported by Doris Duke Charitable Foundation, grant #2013104 and HL122684. J.B.N. is supported by grants from The Danish Heart Foundation, the Lundbeck Foundation and the A.P. Mller Foundation for the Advancement of Medical Science. We thank all participants and staff of the Andr and France Desmarais MHI Biobank, and in particular John D. Rioux and Simon de Denus who participated in the design of the study. The MHI replication study was funded by the MHI Foundation and the Canada Research Chair Program.

We thank all participants and physician investigators of the International Bicuspid Aortic Valve Consortium (BAVCon) for their scientic collaboration. Investigators of the International Bicuspid Aortic Valve Consortium (BAVCon): Eduardo Bossone, Azienda Ospedaliera Universitaria Salerno; Rodolfo Citro, Azienda Ospedaliera Universitaria Salerno; Stefano Nistri, CMSR Veneto Medica, Vicenza; Dan Gilon, Hadassah-Hebrew University Medical Center; Ronen Durst, Hadassah-Hebrew University Medical Center; Simon Body, Harvard University; Thoraf M. Sundt, Harvard University; J. Daniel Muehlschlegel, Harvard University; Carlo de Vincentis, IRCCS Policlinico San Donato; Francesca R. Pluchinotta, IRCCS Policlinico San Donato; Hector I. Michelena, Mayo Clinic; Maurice Sarano, Mayo Clinic; Bo Yang, University of Michigan; Kim Eagle, University of Michigan; Cristen J. Willer, University of Michigan; Giuseppe Limongelli, Monaldi Hospital; Malenka M. Bissell, Oxford University; Alessandro DellaCorte, Second University of Naples; Amalia Forte, Second University of Naples; Gordon Huggins, Tufts University; Victor Dayan, Universidad de la Repblica, Uruguay; Yohan Boss, Universit Laval; Evaldas Girdsaukas, University of Hamburg; Ashutosh Hardikar, University of Tasmania; Thomas Marwick, University of Tasmania; Joseph Bavaria, University of Pennsylvania; Rita C. Milewski, University of Pennsylvania;

Dianna M. Milewicz, University of Texas Health Science Center at Houston;

Siddarth K. Prakash, University of Texas Health Science Center at Houston; Arturo Evangelista, Vall dHebron University Hospital; Joshua C Denny, Vanderbilt University; and Edward Hulten, Walter Reed National Military Medical Center. The BAVCon efforts at Partners HealthCare (Boston, MA) are supported by HL114823 to S.C.B. The ASAP study was supported by the Leducq Foundation (MIBAVA). The Mount Sinai BioMe Biobank is supported by The Andrea and Charles Bronfman Philanthropies. BioMe is partially funded by a grant from NIH (NHGRI U01HG007417). Analyses of BioMe data were supported in part through the computational resources and staff expertise provided by the Department of Scientic Computing at the Icahn School of Medicine at Mount Sinai. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Ofce of the Director of the National Institutes of Health. Additional funds were provided by the NCI, NHGRI, NHLBI, NIDA, NIMH and NINDS. Additional funding was provided to the Broad Institute (HHSN268201000029C), the van Andel Institute (10ST1035, HHSN261200800001E), University of Miami (DA006227, DA033684 and N01MH000028), the University of Geneva (MH090941 and MH101814), the University of Chicago (MH090951, MH090937, MH101820 and MH101825), the University of North CarolinaChapel Hill (MH090936 and MH101819), Harvard University (MH090948), Stanford University (MH101782), Washington University St Louis (MH101810) and the University of Pennsylvania (MH101822). The data used for the analyses described in this manuscript were obtained from the GTEx Portal on 11 16 November 2015.

Author contributions

C.J.W. and B.Y. designed the study. C.J.W., B.Y., W.Z. and J.J. drafted the manuscript. B.Y., J.B.N., M.M., M.H., G.L., L.F., S.P., C.S., L.F., M.L., W.H., A.D., A.L., M.T., L.F.,M.-P.D., E.M.I., A.F.-C., D.-C.G., E.P.B., G.M.D., A.B., S.K., H.C., B.K., J.G., S.G. and K.E. characterized phenotypes of samples. M.O., J.K., Y.E.C. and B.Y. contributed to the discussion of biological mechanisms. W.Z., J.B.N., M.H., G.L., L.F., S.P., C.S., L.F., M.L., H.C. and H.M.K. performed the statistical analysis of association data. G.A.F. and A.P.B. performed the chromatin conformation analysis. J.J. and B.Y. performed experiments of CRISPR-Cas9 on iPSCs and experiments of EndoMT. R.J.F.L., P.E., J.-C.T., C.M.B., D.M., S.C.B., C.J.W. and B.Y. are Principal Investigators of cohorts.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Yang, B. et al. Protein-altering and regulatory genetic variants near GATA4 implicated in bicuspid aortic valve. Nat. Commun. 8, 15481 doi: 10.1038/ncomms15481 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

10 NATURE COMMUNICATIONS | 8:15481 | DOI: 10.1038/ncomms15481 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications