# PHACTR1 Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance

Soto-Romuald Kiando, Nathan R. R Tucker, Luis-Jaime Castro-Vega, Alexander Katz, Valentina d'Escamard, Cyrielle Treard, Daniel Fraher, Juliette Albuisson, Daniella Kadian-Dodov, Zi Ye, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Soto-Romuald Kiando, Nathan R. R Tucker, Luis-Jaime Castro-Vega, Alexander Katz, Valentina d'Escamard, et al. PHACTR1 Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance. PLoS Genetics, 2016, 12 (10), pp.e1006367. 10.1371/journal.pgen.1006367. hal-02379438

HAL Id: hal-02379438

https://hal.science/hal-02379438v1

Submitted on 25 Nov 2019

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# PHACTR1 Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance

Soto Romuald Kiando¹,², Nathan R. Tucker³, Luis-Jaime Castro-Vega¹,², Alexander Katz⁴, Valentina D'Escamard⁵, Cyrielle Tréard¹,², Daniel Fraher³, Juliette Albuisson¹,²,⁶,⁷, Daniella Kadian-Dodov⁵, Zi Ye⁸, Erin Austin⁸, Min-Lee Yang⁴, Kristina Hunker⁴, Cristina Barlassina⁹, Daniele Cusi¹⁰, Pilar Galan¹¹, Jean-Philippe Empana¹,², Xavier Jouven¹,²,¹², Anne-Paule Gimenez-Roqueplo¹,²,⁷, Patrick Bruneval¹,², Esther Soo Hyun Kim¹³, Jeffrey W. Olin⁵, Heather L. Gornik¹³, Michel Azizi²,¹⁴,¹⁵, Pierre-François Plouin¹,²,¹⁴, Patrick T. Ellinor³, Iftikhar J. Kullo⁸, David J. Milan³, Santhi K. Ganesh⁴, Pierre Boutouyrie¹,²,¹⁶, Jason C. Kovacic⁵, Xavier Jeunemaitre¹,²,⁶,⁷, Nabila Bouatia-Naji¹,²*

<figure><img src="image_3.png" /></figure>

CrossMark
click for updates

1. INSERM, UMR970 Paris Cardiovascular Research Center (PARCC), Paris F-75015, FRANCE.
2. Paris-Descartes University, Sorbonne Paris Cité, Paris 75006, FRANCE.
3. Cardiovascular research Center, Massachusetts General Hospital, Charlestown, MA 02114, USA.
4. Department of Internal Medicine and Department of Human Genetics, University of Michigan, Ann Arbor, MI 48109, USA.
5. The Zena and Michael A. Wiener Cardiovascular Institute, Icahn School of Medicine, Marie-Josée and Henry R. Kravis Cardiovascular Health Center at Mount Sinai, One Gustave L. Levy Place, Box 1030 New York, NY 10029, New York, NY, USA.
6. Assistance Publique-Hôpitaux De Paris, Referral Center for Rare Vascular Diseases, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE.
7. Assistance Publique-Hôpitaux De Paris, Department of Genetics, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE.
8. Department of Cardiovascular Diseases, Mayo Clinic, Rochester, MN 55905, USA.
9. Dept. of Health Sciences, Genomic and Bioinformatics Unit, Viale Ortles 22/4, Milano, Chair and Graduate School of Nephrology, University of Milano, Division of Nephrology, San Paolo Hospital, Milano, 20142, ITALY.
10. Institute of Biomedical Technologies, Italian National Centre of Research, Via F.lli Cervi 93, 20090 Segrate - Milano.
11. Nutritional Epidemiology Research Group, Sorbonne-Paris-Cité, UMR University of Paris 13/Inserm U-557/INRA U-1125/CNAM, Bobigny, France F-93017, Bobigny, FRANCE.
12. Assistance Publique-Hôpitaux De Paris, Department of Cardiology, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE.
13. Department of Cardiovascular Medicine, Cleveland Clinic Heart and Vascular Institute, Cleveland, OH 44195, USA.
14. Assistance Publique-Hôpitaux De Paris, Department of Hypertension, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE.
15. INSERM, Clinical Investigation Center CIC1418, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE.
16. Assistance Publique-Hôpitaux De Paris, Department of Pharmacology, Hôpital Européen Georges Pompidou, Paris, F-75015, FRANCE

* nabila.bouatia-naji@inserm.fr

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Kiando SR, Tucker NR, Castro-Vega L-J, Katz A, D'Escamard V, Tréard C, et al. (2016) *PHACTR1* Is a Genetic Susceptibility Locus for Fibromuscular Dysplasia Supporting Its Complex Genetic Pattern of Inheritance. PLoS Genet 12(10): e1006367. doi:10.1371/journal.pgen.1006367

**Editor:** Ruth McPherson, University of Ottawa
Heart Institute, CANADA

**Received:** May 4, 2016

**Accepted:** September 16, 2016

**Published:** October 28, 2016

**Copyright:** © 2016 Kiando et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was supported by a Young Investigator funding GENET-FMD (ANR-13-JSV1-0002-01 to NBN), ANR GDPM2 (2010-BLANC-1129-02 to XJ), National French clinical grant (PHRC ARCADIA-PROFILE to PFP) and Fondation pour la recherche Médicale support (FRM MEDYA to MA and Equipe FRM 2015-2018 DEQ20150331716 to XJ and NBN), in part by a European funding FP7 HYPERGENES (FP7-

## Abstract

Fibromuscular dysplasia (FMD) is a nonatherosclerotic vascular disease leading to stenosis, dissection and aneurysm affecting mainly the renal and cerebrovascular arteries. FMD is often an underdiagnosed cause of hypertension and stroke, has higher prevalence in females (~80%) but its pathophysiology is unclear. We analyzed ~26K common variants (MAF>0.05) generated by exome-chip arrays in 249 FMD patients and 689 controls. We replicated 13 loci (P<10⁻⁴) in 402 cases and 2,537 controls and confirmed an association between FMD and a variant in the phosphatase and actin regulator 1 gene (*PHACTR1*).
HEALTH-F4-2007-201550 to XJ and DC), the National Institutes of Health (R01-HL127692 to DJM, and T32HL007824 to JK and VDE, and 1R01HL092577, R01HL128914, K24HL105780 to PTE), the Doris Duke Charitable Foundation (Grant #:2013104 to SKG), an Established Investigator Award from the American Heart Association (13EIA14220013 to PTE) and by the Fondation Leducq (14CVD01 to PTE). Dr. Gornik and the Cleveland Clinic FMD Biorepository has been supported in part by the National Institutes of Health, National Center for Research Resources, CTSA 1UL1RR024989, Cleveland, Ohio. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

Three additional case control cohorts including 512 cases and 669 replicated this result and overall reached the genomic level of significance (OR = 1.39, P = 7.4×10⁻¹⁰, 1,154 cases and 3,895 controls). The top variant, rs9349379, is intronic to *PHACTR1*, a risk locus for coronary artery disease, migraine, and cervical artery dissection. The analyses of geometrical parameters of carotids from ~2,500 healthy volunteers indicate higher intima media thickness (P = 1.97×10⁻⁴) and wall to lumen ratio (P = 0.002) in rs9349379-A carriers, suggesting indices of carotid hypertrophy previously described in carotids of FMD patients. Immunohistochemistry detected PHACTR1 in endothelium and smooth muscle cells of FMD and normal human carotids. The expression of *PHACTR1* by genotypes in primary human fibroblasts showed higher expression in rs9349379-A carriers (N = 86, P = 0.003). *Phactr1* knockdown in zebrafish resulted in dilated vessels indicating subtle impaired vascular development.

We report the first susceptibility locus for FMD and provide evidence for a complex genetic pattern of inheritance and indices of shared pathophysiology between FMD and other cardiovascular and neurovascular diseases.

## Author Summary

Fibromuscular Dysplasia (FMD) is a vascular disease characterized by a succession of occlusions and dilatation of medium-sized arteries (e.g renal, carotid or coronary arteries) with important health consequences, mainly resistant hypertension and stroke. FMD is an atypical vascular disease because it is not associated with overweight or dyslipidemia and 80% of patients are early middle aged women. Our genetic study conducted in >1100 patients and >3800 controls demonstrate that a common variant rs9349379 located on chromosome 6 in the phosphatase and actin regulator 1 gene (*PHACTR1*) increases by ~40% the risk of FMD. This is the first time a genetic risk factor is reported for FMD because it has been longtime considered rare and potentially under a Mendelian mode of inheritance. We also show that rs9349379 correlates with the expression of *PHACTR1* in fibroblasts from FMD patients and controls. Interestingly, the same allele that increases the risk of FMD is at risk for cervical artery dissection and migraine, often reported in FMD patients but protective from myocardial infarction and coronary disease, where atherosclerosis is more common. The clear role of PHACTR1 in maintaining vascular well integrity is not fully elucidated. Using a specific antibody we detected PHACTR1 both on endothelial and smooth muscle cells of human FMD and control carotids, which suggests that PHACTR1 may have multiple functions depending on the cell type and the degree of atherosclerosis of the arteries.

## Introduction

Fibromuscular dysplasia (FMD) is a non-atherosclerotic and non-inflammatory vascular disease leading to stenosis, aneurysm, dissection, and/or occlusion of medium-sized arteries, in particular the renal and extracranial cerebrovascular arteries [1-3]. FMD predisposes to hypertension, transient ischemic attack and stroke [2,3]. Intriguingly, 75% to 90% of FMD patients are women [3,4] and FMD is increasingly considered to be a silent and under-diagnosed condition [5].
The angiography-based classification of renal FMD distinguishes between patients with multifocal stenoses, including the "string-of-beads" FMD pattern, and unifocal (or focal) FMD with corresponding differences in sex ratio, median age of diagnosis and smoking status [6,7].

The pathogenesis of FMD is unknown and there are strong arguments in favor of a genetic origin. We have previously reported familiality of FMD in ~10% of patients [8]. The US FMD registry has also described a family history in first or second-degree relatives for FMD (7%) and aneurysm (23%) [3]. Nonetheless, the unknown status of most of family members compromised heritability estimates and made the assessment of the genetic mode of inheritance difficult. Recently, we investigated the coding genomes from 16 affected siblings and excluded the existence of a major mutated gene in familial FMD [9]. The complexity of the diagnosis based on computed tomography angiography (CTA) and/or magnetic resonance angiography (MRA) and a general lack of awareness among affected patients and clinicians result in an under-diagnosis bias.

To advance our understanding of the etiology of FMD, we performed a genetic association study that identifies a first genetic susceptibility locus for FMD. In addition to the high prevalence of asymptomatic FMD (~3–6%) [3,10] and the existence of environmental modifiers (e.g female hormones, lifetime mechanical stress) our study provides genetic and functional evidence supporting for the first time a complex genetic basis for FMD.

## Results

### Association of rs9349379 in PHACTR1 with FMD

We carried-out a multi-stage genetic association study, including a discovery and four validation case control cohorts, to identify genetic determinants of FMD. All FMD patients and controls were of European ancestry and have similar overall clinical features (Table 1). First, we analyzed 25,606 common (minor allele frequency ≥ 0.05) genetic variants in 249 FMD cases and 689 controls. Despite the small sample size, we also performed stratified analyses including only females (193 patients and 416 controls). Associations across chromosomes are summarized in Manhattan Plots (S2 Fig). In the global analysis, no SNP achieved the adjusted significance threshold for multiple testing (P = 1.95×10⁻⁶, S2 Fig). Nevertheless, the strongest association signal located on chromosome 6 (rs9349379, effect allele frequency (EAF) in

<table><caption>Table 1. Clinical features of the study populations.</caption><thead><tr><th>Cohorts</th><th>N</th><th>Females (%)</th><th>Age at exam (yrs)</th><th>Multifocals (%)</th><th>HTN (%)</th><th>HTN diagnosis (yrs)</th><th>Study Design</th></tr></thead><tbody><tr><td>RVDRC Cases</td><td>249</td><td>193 (78%)</td><td>44.12 ± 14.31</td><td>164 (66%)</td><td>228 (92%)</td><td>36.6 ± 13.65</td><td>Clinical recruitment</td></tr><tr><td>SU.VI.MAX controls</td><td>689</td><td>416 (60%)</td><td>49.80 ± 6.20</td><td>NR</td><td>NA</td><td>NA</td><td>Population based</td></tr><tr><td>ARCADIA cases</td><td>402</td><td>319 (79%)</td><td>42.46 ± 15.68</td><td>273 (68%)</td><td>341 (85%)</td><td>34.9 ± 15.15</td><td>Clinical recruitment</td></tr><tr><td>PPS3 controls</td><td>2,537</td><td>1012 (40%)</td><td>58.73 ± 5.94</td><td>NR</td><td>0</td><td>NR</td><td>Population based</td></tr><tr><td>Mayo Cases</td><td>143</td><td>119 (83%)</td><td>61.48 ± 13.43</td><td>NA</td><td>108 (76%)</td><td>58.75 ± 10.31</td><td>Clinical recruitment</td></tr><tr><td>Mayo Controls</td><td>333</td><td>286 (86%)</td><td>65.78 ± 10.50</td><td>NR</td><td>169 (51%)</td><td>62.52 ± 10.25</td><td>Clinical recruitment</td></tr><tr><td>UM/Cleveland Cases</td><td>304</td><td>295 (97%)</td><td>53.93 ± 10.45</td><td>254 (84%)</td><td>193 (63%)</td><td>NA</td><td>Clinical recruitment</td></tr><tr><td>UM/Cleveland Controls</td><td>289</td><td>280 (97%)</td><td>55.13 ± 9.97</td><td>NR</td><td>NA</td><td>NA</td><td>Clinical recruitment</td></tr><tr><td>DEFINE-FMD Cases</td><td>65</td><td>65 (100%)</td><td>57.03 ± 10.13</td><td>65 (100%)</td><td>41 (63%)</td><td>42.68 ± 13.05</td><td>Clinical recruitment</td></tr><tr><td>DEFINE-FMD Controls</td><td>47</td><td>47 (100%)</td><td>50.70 ± 9.38</td><td>NR</td><td>4 (9%)</td><td>55.25 ± 13.05</td><td>Clinical recruitment</td></tr></tbody></table>

HTN, Hypertension

NR, Not Relevant

NA, Not Available.

doi:10.1371/journal.pgen.1006367.t001
cases = 0.70, odds ratio (OR) = 1.65, *P* = 1.47×10⁻⁵) surpasses the adjusted threshold in the female only analysis (OR = 1.99, *P* = 8.16 × 10⁻⁷, S2 Fig, S1 Table).

In addition to rs9349379, we selected for follow-up loci that showed suggestive association with FMD (*P*<10⁻³) and were located in or near either biological candidate genes (e.g extracel- lular matrix degradation) or previous cardiovascular genome-wide association (GWAS) sig- nals. We also prioritized four SNPs located within one megabase interval around rs9349379, the top associated variant. The first follow-up study included 402 FMD patients from the ARCADIA registry, who had similar clinical characteristics to the patients of the discovery stage and were compared to 2,537 controls from PPS3 (Table 1). Of the 16 SNPs selected, 13 passed genotyping QC criteria in cases and controls. Three SNPs, all located in the *PHACTR1* locus, showed replicated association with FMD (rs9349379-A, *P* = 7.21 × 10⁻⁴, rs9369640-C, *P* = 8.45 × 10⁻⁴ and rs1332844-C, *P* = 1.72 × 10⁻³, S1 Table). Conditioned regression analyses for rs9349379 in the case control cohort including the discovery and first follow-up study sam- ples indicated association signal redundancy for rs9369640 (*P* = 0.21) and rs1332844 (*P* = 0.24) with rs9349379 being the most statistically significant. Next, we investigated the association of rs9349379 with FMD in three additional and independent case-control cohorts from the USA. Overall, of the four follow-up studies, three showed a significant effect of rs9349379 on the risk of FMD; all studies showed consistent direction of effect of the FMD risk allele rs9349379-A being more prevalent in FMD cases (Table 2). The association meta-analysis included 1,154 FMD patients and 3,895 controls and indicated an overall OR of 1.39 for rs9349379 (EAF = 0.69 in global cases sample), with a high level of significance (*P* = 7.36×10⁻¹⁰) that is below the genomic threshold and no indices for heterogeneity between studies (*P* = 0.574, Table 2).

As for rare variants (MAF < 5%), we observed 62,767 polymorphic variants in the discovery case control analyses that we analyzed using the SKAT-O gene-based association test. None of the 9,967 genes covered with at least 2 polymorphic variants showed significant association

<table><caption>Table 2. Association of rs9349379 in *PHACTR1* with FMD in five independent case-control cohorts.</caption><thead><tr><th rowspan="2">Study</th><th rowspan="2">Cohorts</th><th rowspan="2">n</th><th colspan="4">Additive model</th><th colspan="2">Recessive model</th></tr><tr><th>AA</th><th>AG</th><th>GG</th><th>EAF (A)</th><th>OR (95% CI)<sup>a</sup></th><th>P Value<sup>a</sup></th></tr></thead><tbody><tr><td rowspan="2">Discovery</td><td>RVDRC Cases</td><td>249</td><td>124</td><td>103</td><td>22</td><td>0.70</td><td>1.65 (1.32–2.07)</td><td>1.47 × 10<sup>-5</sup></td></tr><tr><td>SU.VI.MAX controls</td><td>689</td><td>237</td><td>341</td><td>111</td><td>0.59</td><td></td><td>1.87 (1.39–2.52)</td></tr><tr><td rowspan="8">Follow-up</td><td>ARCADIA</td><td>393</td><td>183</td><td>170</td><td>40</td><td>0.68</td><td>1.32 (1.12–1.54)</td><td>7.21 × 10<sup>-4</sup></td></tr><tr><td>PPS3 controls</td><td>2537</td><td>982</td><td>1174</td><td>381</td><td>0.62</td><td></td><td>1.38 (1.11–1.71)</td></tr><tr><td>Mayo Cases</td><td>143</td><td>62</td><td>65</td><td>16</td><td>0.66</td><td>1.37 (1.02–1.83)</td><td>0.034</td></tr><tr><td>Mayo Controls</td><td>333</td><td>116</td><td>159</td><td>58</td><td>0.59</td><td></td><td>1.94 (0.91–4.16)</td></tr><tr><td>UM/Cleveland Cases</td><td>304</td><td>145</td><td>130</td><td>29</td><td>0.69</td><td>1.31 (1.03–1.66)</td><td>0.026</td></tr><tr><td>UM/Cleveland Controls</td><td>289</td><td>116</td><td>135</td><td>38</td><td>0.63</td><td></td><td>1.43 (0.96–2.14)</td></tr><tr><td>DEFINE-FMD Cases</td><td>65</td><td>41</td><td>18</td><td>6</td><td>0.77</td><td>1.39 (0.77–2.49)</td><td>0.276</td></tr><tr><td>DEFINE-FMD Controls</td><td>47</td><td>22</td><td>22</td><td>3</td><td>0.70</td><td></td><td>1.38 (0.99–1.91)</td></tr><tr><td rowspan="2">Meta-analysis<sup>b</sup></td><td>All Cases</td><td>1154</td><td>555</td><td>486</td><td>113</td><td>0.69</td><td>1.39 (1.25–1.54)</td><td>7.36 × 10<sup>-10</sup></td></tr><tr><td>All Controls</td><td>3895</td><td>1473</td><td>1831</td><td>591</td><td>0.61</td><td></td><td>1.50 (1.3–1.73)</td></tr></tbody></table>

OR, Odds Ratio

CI, Confidence Interval

EAF, Effect Allele Frequency.

<footnote><sup>a</sup>Odds Ratio (OR) and P values were computed by logistic regression under the additive and recessive genetic model.</footnote>

<footnote><sup>b</sup>Meta-analysis was performed using inverse variance-weighted method. Heterogeneity between cohorts was tested using Cochran's Q statistics and was not significant (*P*<sub>add</sub> = 0.574; *P*<sub>rec</sub> = 0.483).</footnote>
with FMD after Bonferroni correction (P = 5.02 × 10⁻⁶, strongest gene P = 1.42 × 10⁻⁴), or bio-
logical candidacy among best scoring genes, which discourages follow-up genotyping or
sequencing efforts in larger samples.

Association of rs9349379 with carotid parameters in PPS3

Here we aimed to characterize the association of rs9349379 with artery thickness and stiffness using a non-invasive method in 2,458 healthy volunteers. Genetic association analyses indicate significant association between the FMD increasing risk allele rs9349379-A and greater intima-media thickness (IMT) (β<sub>add</sub> = 11.65 µm, P<sub>add</sub> = 1.65×10<sup>-4</sup>, Fig 1) and wall to lumen ratio (WLR) (β<sub>add</sub> = 0.004, P<sub>add</sub> = 0.002, Fig 1), and decreased circumferential wall stress at DBP (β<sub>add</sub> = -0.72, P<sub>add</sub> = 0.004, S2 Table). Interestingly, despite their smaller number (N = 975) compared to males (N = 1,483), females present more accentuated effect of the rs9349379-A allele on IMT (β<sub>add</sub> = 14.83 µm, P<sub>add</sub> = 0.001) and WLR (β<sub>add</sub> = 0.006, P<sub>add</sub> = 0.003). However, we detected no significant interaction with sex. rs9349379-A also associated with increased wall cross-sectional area (β<sub>add</sub> = 0.24 mm², P<sub>add</sub> = 8.67×10<sup>-4</sup>, S2 Table), peripheral SBP (β<sub>add</sub> = 0.70 mmHg, P<sub>add</sub> = 0.009, S2 Table) and central pulse pressure (β<sub>add</sub> = 0.62 mmHg, P<sub>add</sub> = 0.002, Fig 1).

PHACTR1 expression studies

Variant rs9349379 resides in the fourth intron of the phosphatase actin regulator 1 gene (*PHACTR1*) in a putative regulatory sequence. We assessed *PHACTR1* expression using quantitative real time PCR in mRNAs from primary cultured human fibroblasts. *PHACTR1* expression did not differ in fibroblasts from FMD patients compared to age and sex matched controls (Fig 2). However, stratifying by rs9349379 genotype indicated increased expression in individuals carrying the FMD risk allele rs9349379-A (*P*<sub>add</sub> = 0.003, Fig 2). At the protein level, using a specific antibody, we found that PHACTR1 is expressed in endothelial and medial smooth muscle cells of carotid arteries from both normal and FMD patients (Fig 2).

Functional role of *Phactr1* in vascular development

Given the unclear role for *PHACTR1* in vascular development and maintenance, we assessed the effect of *PHACTR1* perturbation in the zebrafish. Morpholino injected embryos had undetectable *phactr1* transcript levels at 72 hours post fertilization while overall morphology was unchanged (S3 Fig). Compared to control injected zebrafish, *phactr1* knockdown resulted in a marked disorganization of the developing hepatic portal vein and the segmental vessels in the developing trunk (Fig 3 and S3 Fig). Further analysis of the diameter of three major peripheral vessels, the dorsal aorta, caudal artery and posterior cardinal vein, demonstrated a nearly 8% dilatation (*P*<0.05) of the posterior cardinal vein in *phactr1* suppressed embryos.

Discussion

Our study describes for the first time the genetic association of rs9349379, a common variant in *PHACTR1*, with arterial fibromuscular dysplasia (FMD). We demonstrate that this common variant increases by ~40% the risk of FMD in five independent case-control studies. These findings are based on genetic data from 1,154 FMD patients and 3,895 controls, the largest investigation conducted so far to elucidate the genetic basis of this intriguing vascular disease with unknown etiological origin and challenging clinical features.

To date, the genetic investigation on FMD was limited to the screening of candidate genes involved in rare vascular and arterial syndromes, and/or the study of underpowered series of
Intima-media thickness (mm x 10⁻³)

<figure><img src="image_5.png" /></figure>

Wall to Lumen Ratio

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

**Fig 1. Association between rs9349379 and artery thickness in healthy controls.** Intima-media thickness and wall to lumen ratio are presented by genotype in all (A,D) women (B,E) and men (C,F) participants. We indicated p-values for the linear regression analyses under the additive model adjusted for age, sex and body surface area.

doi:10.1371/journal.pgen.1006367.g001

patients.[9,11,12] For long considered a rare disease, FMD has also been hypothesized to be under the genetic control of highly penetrant genetic defects. Our study provides the first evidence for a complex genetic pattern of inheritance for FMD, involving a common genetic allele (frequency = 0.60 in the general population). The identification of a genetic susceptibility locus for FMD supports the concept that this disease is controlled by a large number of genetic
<figure><img src="image_11.png" /><figcaption>**Fig 2.** *PHACTR1 mRNA expression and immunostaining.* (A) *PHACTR1* mRNA expression in human fibroblasts. *PHACTR1* mRNA levels were determined by RT-qPCR in cultured fibroblasts from controls (N = 39) and FMD cases (N = 51). *GAPDH* expression was used as control for normalization and data are expressed as mean fold change of *PHACTR1* relative to *GAPDH*. There is a non-significant trend toward overexpression of *PHACTR1* in FMD cases compared to controls (P = 0.61, Mann-Whitney test) whereas significant differences are uncovered after stratification by genotypes (P = 0.003, Kruskal-Wallis test). (B) Immunostaining of normal and FMD internal carotid using anti-PHACTR1 antibody. PHACTR1 was detected in endothelium and smooth muscle cells in the media. Staining is mostly cytoplasmic with regular alignment in normal carotids and typical disorganized cellular structure in media of FMD carotid.</figcaption></figure>

doi:10.1371/journal.pgen.1006367.g002

determinants in strong interaction with various environmental factors, including female sex and mechanical stresses.

We provide genetic evidence for increased IMT, narrowed artery lumen and no change in arterial stiffness, even after adjusting for blood pressure, in healthy volunteers carrying the rs9349379-A, the FMD risk allele. These genetic associations support the observation of carotid concentric hypertrophy in FMD risk carriers, especially females. This finding is consistent with previous arterial features reported in FMD patients when compared to age, sex and SBP matched controls.[13] Circumferential wall stress was decreased in rs9349379-A carriers, showing that hypertrophy may overcome the moderate increase in blood pressure. Of note, the effect size of the FMD risk allele on IMT corresponds to approximately five years of aging for the IMT, when reported to sex and age reference values for arterial geometry,[14] suggesting potentially an accelerated arterial aging in carriers of the rs9349379-A FMD risk allele.

*PHACTR1* was previously identified by several genome-wide association studies (GWAS) as a risk locus for cardiovascular and neurovascular diseases. *PHACTR1* is a confirmed susceptibility locus for coronary artery disease (CAD) and myocardial infarction (MI),[15,16] migraine [17,18] and more recently cervical artery dissection (CeAD), a rare condition defined as a mural hematoma in a carotid or vertebral artery and a cause of stroke.[19] Of note, the association of rs9349379 is in the opposite direction for FMD, CeAD and migraine, with rs9349379-A at risk, when compared to CAD and MI with rs9349379[G] allele at risk. Migraine and CeAD
<figure><img src="image_12.png" /><figcaption>**Fig 3.** *phactr1* modulation in zebrafish affects vascular dimensions and patterning. Two-dimensional projections obtained from z-series confocal images in the trunk of control and phactr1 knockdown (KD) zebrafish embryos (two representative images per condition). Green represents the vascular endothelium as marked by EGFP. Greyscale represents the corresponding DIC bright field image of the fish trunk. DA: Dorsal aorta, SV: Segmental vessel, PCV: Posterior cardinal vein. Quantification of inner vascular diameter for the dorsal aorta (DA), posterior cardinal vein (PCV) and caudal artery (CA). (*) represents *P*<0.05.</figcaption></figure>

doi:10.1371/journal.pgen.1006367.g003

share several clinical features with FMD. Migraine is also more prevalent in females than in males and it is reported by a third of FMD patients.[3] CeAD is an important risk factor for subarachnoid hemorrhage stroke, as is cerebrovascular FMD.[3] The investigation of a large series of CeAD patients indicated that cervical FMD is reported in 5.6% of the patients[20] and carotid dissection is the presenting manifestation in 12.1% of patients of the US FMD registry.
[3] In contrast to migraine and CeAD, the implication of the same common variant in *PHACTR1* in CAD/MI and FMD is rather unexpected, though it involved a different allele at the same genetic variant. The clinical link between FMD and CAD/MI is less obvious, except for the high proportion of FMD in rare forms of MI involving spontaneous coronary artery dissection, that also present in young females without atherosclerotic risk factors.[21–23] By
definition FMD does not involve atherosclerotic stenosis that is concomitant with CAD and MI through pathogenesis processes (e.g dyslipidemia and inflammation). Future comprehensive genetic investigation by full GWAS for FMD will allow the assessment of the putative consistency and opposition of effects between risk loci for these cardiovascular and neurovascular diseases.

The rs9349379 is intronic to *PHACTR1* in a putative noncoding regulatory sequence. *PHACTR1* encodes a phosphatase and actin regulator protein and its function is not fully elucidated. Here we describe significant correlation between rs9349379 genotypes and *PHACTR1* expression in human fibroblasts from FMD patients and controls. This association is consistent with the expression quantitative trait loci (eQTL) data from several artery beds of the Genotype-Tissue Expression (GTEx) project and a recent study where increased expression is reported in human coronary arteries from donors carrying the rs9349379-A.[16] Using a genome-editing technique with CRISPR-Cas9 applied to human umbilical vein endothelial cells, the deletion of the sequence containing rs9349379 caused a 35% decrease in *PHACTR1* expression and impaired the fixation of the myocyte enhancer factor-2 (MEF2), supporting a functional regulatory effect of this genetic variant *in vitro*. [16] However, other transcription factors may also bind to this site, as the knockdown of the expression of two members of the MEF2 family did not change *PHACTR1* expression.[24] We note that ENCODE ChIP-Seq data derived from genomes extracted from two immortalized cell lines after stimulation by estradiol indicates several consistent peaks for ESR1 suggesting a putative regulation of PHACTR1 by this female hormone. This putative regulation that needs to be established experimentally is consistent with the high proportion (~75–90%) of women among FMD patients.

Molecular studies have also linked PHACTR1 to cell adhesion and migration in angiogenesis via vascular endothelial growth factor (VEGF) stimulation.[25,26] In contrast to other studies,[26,27] Beaudoin et al reported a lack of evidence for the induction of PHACTR1 expression in endothelial cells with pro-angiogenic (VEGF), pro-inflammatory stimulations or shear stress.[16] Our immunohistochemistry staining of PHACTR1 indicates the presence of PHACTR1 both in endothelial and smooth muscle cells of normal and FMD carotid arteries. The genetic implication of *PHACTR1* in FMD, the subtle impairment in the development of vasculature in zebrafish, in addition to the evidence from recent GWAS describing an increasing number of loci with genes involved in vessel wall biology in CAD and MI[28] all support PHACTR1 plays a key etiological role in vascular structure. Recent work on *PHACTR1* regulation in atherosclerosis showed strong expression in human atherosclerotic plaque macrophages lipid-laden foam cells, adventitial lymphocytes and endothelial cells.[29] This group also describes the absence of PHACTR1 in SMCs of healthy and atherosclerotic aortas, which is not supported by our immunostaining on medium-sized arteries using a different anti-body, though they have detected the expression of an intermediate transcript in SMCs. Further investigation is required to understand this discrepancy that might reflect *PHACTR1* is much more abundantly expressed in macrophages from atherosclerotic plaques, where this population of cells is highly represented and biologically active, which is different from the FMD and healthy vessels without atherosclerosis. Thus, according to rs9349379 genotype, the manifestation of one specific vascular disease depends on multiple and particular environmental triggers (e.g. female specific hormonal context, mechanical movements of medium size arteries in the case of FMD) that are still to be specifically determined for each of these cardiovascular and neurovascular diseases.

In summary, our study provides genetic and functional evidence to support *PHACTR1* as a first susceptibility locus for FMD. Further functional exploration of this locus and more comprehensive genetic investigation through genome-wide association will provide additional predisposing loci to FMD and help understanding the etiological mechanisms of non-atherosclerotic arterial stenosis.
# Methods

## Ethics Statement

We obtained individual written informed consent from all participants included in France and the USA case control studies.

French FMD patients (RVDRC and ARCADIA) from the European Hospital Georges Pompidou (HEGP) are part of the ARCADIA/PROFILE protocol that was approved by the Ile-De-France research ethics committee (Comité de Protection des Personnes: CPP d'île de France) on 03/04/2009 (ID: 2009-A00288-49). The ethics committee of the Paris-Cochin hospital approved the protocols of the SU.VI.MAX study (CCPPRB No. 706) and PPS3 (CPP No. 2007-A01386-47). PPS3 is declared at ClinicalTrials.gov (Identifier: NCT00741728). Mayo Clinic case control study is part of the Vascular Diseases Biorepository study and was approved by the Mayo I Mayo Clinic Institutional Review Boards (IRB # 08-008355). UM/Cleveland case control study was approved by each institution IRB protocols: the University of Michigan IRB number HUM00044507 and the Cleveland Clinic IRB number 10-318. The DEFINE-FMD study was approved by the Human Research Ethics Committee of the Icahn School of Medicine at Mount Sinai (Study ID: HS#13-00575/GCO#13-1118 and is registered with ClinicalTrials.gov Identifier: NCT01967511.

## Patients and control populations

We used a three-stage association design. First, we performed an exome-chip based genetic association in 249 French FMD patients (RVDRC cohort) and 689 controls from SU.VI.MAX [30]. Second, we followed-up 13 loci in an independent set of 402 French patients (ARCADIA registry) and 2,537 controls from PPS3 [31]. Three additional studies from the USA totaling 512 patients and 669 controls were used for further replication: Mayo Clinic cohort [32], University of Michigan (UM)/Cleveland Clinic cohort and the DEFINE-FMD study. All participants are of European ancestry and presented similar clinical characteristics (Table 1) and homogeneous diagnosis, exclusion and inclusion criteria.

**RVDRC and ARCADIA cases.** We analyzed unrelated FMD followed-up at the Rare Vascular Diseases Reference Center (RVDRC) of the European Hospital Georges Pompidou (HEGP), Paris, France. We ascertained patients from the ARCADIA (Assessment of Renal and Cervical Artery Dysplasia) register, an ongoing national FMD registry at the HEGP, Paris. The diagnosis of FMD in RVDRC and ARCADIA patients was established using clinical information from the medical history, the interpretation of angiography and/or computed tomography scan of arterial beds after the exclusion of other causes of arterial stenosis such as atherosclerosis, Takayasu disease and Elhers Danlos syndrome and neurofibromatosis type 1. Given the complexity of the interpretation of imaging of vascular diseases, a local panel of experts including clinicians from the departments of hypertension, radiology, vascular medicine and medical genetics validated the diagnosis of FMD.

**SU.VI.MAX and PPS3 controls.** Controls in discovery stage are French Europeans from the SU.VI.MAX study as previously described [30]. PPS3 controls were 2,537 participants ascertained from the Paris Prospective Study 3 (PPS3), an ongoing observational French prospective study evaluating the possible implication of numerous vascular health parameters in cardiovascular disease in healthy subjects aged from 50 to 75 years undergoing a standard health examination in a large preventive medical center that is subsidized by the French national insurance system for salaried workers [31].

**Mayo Case-Control study.** Potential FMD cases were identified from the Mayo vascular disease bio-repository. Electronic health records of 501 patients were manually reviewed by
clinicians and the diagnosis was confirmed for 191 patients according to criteria in [1], of whom 159 had DNA available and 143 (83% females) had high-density genotyping data. Patients were matched for sex and age (+/- 2 years) to 333 controls without known FMD or other atherosclerotic vascular diseases.

**UM/Cleveland Clinic case control study.** Patients were recruited at the University of Michigan (UM) and/or the Cleveland Clinic. UM FMD cases were recruited into an IRB-approved study through a referral clinic at the University of Michigan and through self-referral to the study. The Cleveland Clinic cases were enrolled among consecutive patients seen at a dedicated FMD referral clinic. Clinical diagnosis of FMD was ascertained by a vascular medicine specialist after review of diagnostic imaging and prior to blood sample collection. DNAs from healthy controls without vascular disease were obtained from the Cleveland Clinic GeneBank, which was approved by the Cleveland Clinic IRB. Genomic DNA was isolated from a peripheral blood sample and analyzed as described[33].

**The DEFINE-FMD Study.** Eligible cases were females with an imaging-confirmed diagnosis of multifocal FMD and who fulfilled other accepted diagnostic criteria [1]. Healthy controls were matched to FMD cases according to age and sex, required to be receiving ≤ 2 blood pressure medications, have a body mass index < 35kg/m2 and to be non-smokers. Healthy controls underwent physical exam and those with bruits; unexplained hypertension or other cardiovascular findings were excluded. Exclusion criteria for all subjects included male gender, multifocal FMD, use of immunosuppressive agents, major comorbidities, diseases that may confound genetic/genomic analyses (i.e. Crohn's disease, multiple sclerosis etc.) or any other form of heritable vascular disease (i.e. Ehlers-Danlos, Marfan, Loeys-Dietz).

## Genotyping

Genotypes were generated using the Illumina-HumanExome-12v1 array in RVDRC, SU.VI. MAX and PPS3,[34] Illumina-Human-Omni-Express-Exome in the DEFINE-FMD Study participants, Illumina-Infinium-Human-CoreExome in Mayo Clinic cohorts and by individual genotyping in ARCADIA (KASP technology) and UM/Cleveland Clinic cohorts (Taqman). We applied quality control (QC) filters to the discovery cases and controls as recommended [35] (individual's call rate < 97%, extreme heterozygosity, sex-discordant, duplicate or relatedness unsing PLINK (version 1.07).[36] Individuals with non-European ancestry were detected and excluded (54 cases and 9 controls) using the EIGENSTRAT program[37] and visualized by principal components analysis (PCA) including HapMap phase 3 samples. Cases and controls displayed a comparable distribution after ancestry QC (S1 Fig). From 240,748 successfully genotyped SNPs, we excluded monomorphic variants (n = 129,890), call rate < 99%, deviation from Hardy-Weinberg Equilibrium in cases and/or controls (P < 10⁻⁵) and minor allele frequency (MAF)<0.05 in controls. The final analysis included 25,606 variants (25,138 autosomal and 468 X-linked) in a sample of 249 FMD cases and 689 controls. Comparable individuals and markers QC was applied to follow-up cohorts using arrays (PPS3, Mayo Clinic case control cohorts, the DEFINE-FMD study) or individual genotyping (ARCADIA and UM/Cleveland cases control cohort).

## Carotid geometry and stiffness measurements in PPS3 cohort

Carotid parameters were measured and calculated as previously reported [34]. Briefly, a 10 MHz 128 transducer linear array probe was positioned on the carotid area. Measurements were performed on a 4 cm segment of the right common carotid artery, 1 cm proximal to the bifurcation/sinus throughout the cardiac cycle for 6 seconds. A longitudinal section showing clear interfaces for blood/intima and media/adventitia was obtained. The system allows real-
time radiofrequency signal analysis with operator-independent determination of external diameter (Dext), internal diameter (Dint), and intima-media thickness (IMT) on 128 lines throughout the cardiac cycle. Distension was measured on 14 lines at high-pulsed radiofre- quency (600 Hz). The axial resolution was 34 µm for diameter, 17 µm for IMT, and 1.7 µm for distension.[38] Aortic blood pressure was estimated from the distension waveform according to van Bortel et al.[39] The distensibility coefficient, representing the elastic properties of the artery as a hollow structure, was calculated as dLCSA/(LCSA×central PP), where LCSA is the lumen cross-sectional area and central PP is the aortic pulse pressure. Carotid stiffness (Cstif) was calculated as DC⁻⁰.⁵ and circumferential wall stress as diastolic blood pressure×Dint/ 2×IMT. WCSA is wall cross-sectional area. Young's elastic modulus is calculated as DC⁻¹x3(1 +LCSA/WCSA). Here we examined the association of one genetic variant (rs9349379) with 14 interdependent traits in a sub-sample of 2458 participants (975 females and 1483 males). Hypertensive participants (BP over 140 and/or 90 mmHg, and/or use of antihypertensive treat- ments) were excluded to avoid confounding with genetic effects.

Statistical analyses

We tested the association with FMD using logistic regression under the additive genetic model as implemented in PLINK[36] (version 1.07) in discovery and follow-up studies. In the discov- ery analysis, we included the first five principal components axes as co-variates to control for hidden population stratification. The Bonferroni adjusted threshold for significance was set to P = 1.95×10⁻⁶ to account for multiple testing of 25,606 common variants. We used the inverse variance-weighted method for meta-analysis implemented in Metal [40]. Gene-based analyses of rare variants was performed using SKAT-O method as described in.[41] Heterogeneity was assessed with Cochran's Q statistics. We tested the association of rs9349379 with carotid parame- ters in 2,458 normotensive subjects (975 females and 1483 males) from PPS3 using a linear regression on an additive genetic model including age, sex, body surface area (BSA), and mean blood pressure (MBP) as covariates when relevant. Before performing the analyses, all parameters (all quantitative traits) were quantile-transformed to a standard normal distribution. The 14 traits fall into three main carotid parameters categories: geometry (intima-media-thickness, external and internal carotid diameters, wall to lumen ratio and circumferential wall stress), arterial stiff- ness (stiffness, cross-sectional distensibility, Young's elastic modulus, cross-sectional compliance, and wall cross-sectional area) and central blood pressure (systolic, diastolic and aortic blood pres- sure). Considering that carotid parameters are interdependent, we applied a Bonferroni correc- tion (P<adj = 0.017) for the three main categories as recommended [34].

PHACTR1 expression in fibroblasts

We analyzed the expression of PHACTR1 from fibroblasts of skin biopsy samples of 104 indi-
viduals (51 FMD, 39 controls, 12 undetermined). Genotypes and detectable expression levels
were available in 86 individuals.

Fibroblasts of the DEFINE-FMD Study were derived from skin biopsy samples using stan-
dard explant techniques from 103 individuals (46 FMD and 57 controls). Briefly, skin biopsy
samples were dissected into small pieces and cultured under cover slips in DMEM/F12 media
containing 20% fetal bovine serum, 1% antibiotic-antimycotic solution, 1% 200 mM L-gluta-
mine, 1% 100mM Sodium pyruvate, and 1% MEM Non-Essential Amino Acids (all from Life
Technologies, Grand Island, NY, USA) at 37°C in 5% CO2. Cell culture medium was replaced
every 48 hours and confluent cells from passages 2–3 were used for RNA extractions. RNA was
then converted to cDNA using iscript cDNA synthesis kit (Bio Rad, Hercules, CA, USA) and
quantitative real time PCR (qRT-PCR) was performed using SYBR green fastmix low rox
(Quanta Biosciences, Gaithersburg, MD, USA). Fibroblasts from the University of Michigan were obtained through the Coriell Biorespository (Camden, NJ, USA) for 20 subjects with FMD and apparently 20 healthy control subjects and matched for passage number. DNAs and genotypes were not available. Confluent cells from passages 6–8 were used to extract RNA and generate cDNAs to assess *PHACTR1* expression according to culture and extraction protocol described for the DEFINE-FMD study. Fibroblasts from the Department of genetic, HGEP study were derived from skin biopsy samples from patients followed up for unknown arterial diseases where FMD and Elhers Danlos syndrome were discarded. FMD status is not available. SNP rs9349379 was genotyped by direct forward and reverse sequencing (BigDye Terminator kit v3.1 cycle sequencing kit) and run on an ABI Prism 3730XL DNA Analyzer Sequencer (Perkin Elmer Applied Biosystems, Foster City, CA). Cell culture and RNA extraction conditions and reagent were identical to the DEFINE-FMD study. Informed consent was obtained from all patients of all centers.

In total, we analyzed the expression of *PHACTR1* from fibroblasts of skin biopsy samples of 104 individuals (51 FMD, 39 controls, 12 undetermined), and 86 individuals had genotypes and relevant expression levels. The expression of PHACTR1 transcript variant 1 (NCBI Reference Sequence: NM_030948.2) and *GAPDH*, which was used as normalization control, was determined by RT-qPCR following the MIQE guidelines.[42] All expression studies from the three laboratories used the same following primers: PHACTR1 fwd ATGACCGCAGGGCAG ATAAG, rev TTCGGATGGCAGCTTTGTCT; GAPDH fwd GGGTGTGAACCATGAGAAG TATGA, rev GGTGCAGGAGGCATTGCT. The efficacy and efficiency of amplifications for PHACTR1 and GAPDH were equivalent as determined by the linearity tests: (*PHACTR1* efficacy 1.97 and efficiency 98.7%); (GAPDH efficacy 2.04 and efficiency 102%). The specificity of amplification products was verified by the presence of a single peak in the melting temperature curve analysis. Outliers displaying *GAPDH* Cq values different from the mean or a ΔCq value >14 (Cq value PHACTR1 –Cq value *GAPDH*) were filtered out. Good correlations for RT and technical qPCR replicates were obtained. Given that data obtained from three independent laboratories was equivalent: mean Cq for *GAPDH* 17,41 +/- 0,65 and mean Cq for *PHACTR1* was 30.46 +/- 1.19, we pooled them as indicated to gain power, especially for the correlation by genotypes where few samples per study were GG carriers. Data is presented as the mean fold change in the expression of PHACTR1 relative to GAPDH, calculated using the formula $2^{-\Delta Cq}$, as recommended.[43] Statistical analyses were performed using a non-parametric Mann-Whitney test for single comparison between FMD cases and controls and Kruskal Wallis test for comparisons according to rs9349379 genotypes.

## PHACTR1 Immunostaining

Protein staining for PHACTR1 was performed on paraffin embedded artery samples using primary anti-PHACTR1 antibody (Sigma-Aldrich, St Louis, MI, USA) and revealed using an ABC peroxidase kit with diaminobenzidine (Vector laboratories, Burlingame, CA, USA). The paraffin blocks of arterial tissues were obtained from surgical pathology archives of the HEGP as remnants of the regular diagnostic procedure (2 renal arteries, 1 internal carotid, 1 popliteal and 1 femoral artery). FMD arteries belonged to patients who had surgical reparation for aneurysm resection (1 carotid, 1 radial, 1 coeliac, 1 ulnar and 1 pancreaticoduodenal artery). The arterial tissues were fixed in formalin and embedded in paraffin. Antigen retrieval were performed by incubating tissue sections in alkaline solution (Dako, Trappes, France) for 40 minutes at 94°C in a hot water bath. The sections were then incubated for 60 minutes with the primary anti-PHACTR1 rabbit polyclonal antibody diluted at 1/50. For reparation we used ABC peroxidase kit with diaminobenzidine (Vector laboratories, Burlingame, CA, USA).
## Zebrafish experiments

Anti-*phactr1* morpholino was designed to target the exon7/intron7 boundary of zebrafish *phactr1* transcript. A non-targeting morpholino of equivalent length but differing nucleotide composition was injected at equivalent concentrations as a control. For injections, 0.8nL of control of anti-*phactr1* morpholino was injected into single cell embryos of the Tubingen/AB strain. For confirmation of morpholino efficacy, primers spanning the targeted boundary were used to amplify cDNA constructed from the isolated RNA of injected 72 hour post fertilization embryos. All the morphological analyses as well as the knockout status of the fish were performed blinded to the experimenter. For imaging studies, fish were reared under standard conditions with addition of phenylthiourea to media in order to inhibit pigmentation. At 72hpf, high-speed videography was performed as described[44] used to determine the diameter of major vessels. Measurement of vessels was manually performed using ImageJ. For confocal reconstructions of the embryonic vasculature, embryos obtained from lines expressing EGFP driven by the *flk* promoter were used. Confocal images were obtained on a Nikon A1SiR confocal. Two-dimensional projections of maximal intensity (head and torso) or averaged intensity (trunk) were used to visualize the overall vascular architecture using ImageJ.

## Supporting Information

**S1 Fig. Principal component analysis using genotypes from FMD cases and controls of the discovery stage.** Data used is post quality control filtration for non-European ancestry origin and relatedness
(TIFF)

**S2 Fig. Association between Single-Nucleotide Polymorphisms (SNPs) and Fibromuscular Dysplasia (FMD).** (A) Association between SNPs and FMD in global analysis. (B) Association between SNPs and FMD in females only. The red line indicates the Bonferroni adjusted significance threshold for 25,606 common variants tested set to $P = 1.95 \times 10^{-6}$.
(TIFF)

**S3 Fig. Phenotypic evaluation of *phactr1* knockdown zebrafish.** (A) RT-PCR evaluation of splice alteration observed following microinjection of phactr1 morpholino (PHACTR1 KD). (B) Brightfield micrographs of overt morphology at 60 hours post fertilization. (C) Two-dimensional projections obtained from z-series confocal images in the head and trunk of control and phactr1 knockdown zebrafish embryos. Green represents the vascular endothelium as marked by EGFP. Greyscale represents the corresponding DIC brightfield image of the fish head and trunk region. HPV indicates the relative position of the developing hepatic portal vein
(TIFF)

**S1 Table. Association of SNPs selected for follow-up in ARCADIA and PPS3 case control analyses.** SNP, Single Nucleotide polymorphism; Chr, Chromosome; EA, Effect allele; OR, Odds-Ratio; CI, Confidence Interval. <sup>a</sup>For genic SNPs, the relevant gene is listed; for intergenic SNPs nearest upstream and downstream genes are listed. <sup>b</sup>SNPs added despite a P value in the discovery stage were above the threshold given their location in or near *PHACTR1*, the most associated locus.
(DOCX)

**S2 Table. Association of rs9349379 in *PHACTR1* with carotid artery parameters in PPS3 controls.** <sup>a</sup>All (N = 2,458), Female (N = 975) and Males only (N = 1,483) participants of the PPS3. $\beta_{q.t.}$ is the effect size estimated using quantile-transformed values and $\beta$ is the effect size estimated using untransformed values. P values indicated are for the analyses of quantile
transformed values. Age, sex, body surface area (BSA), smoking status and cholesterol were included as covariables in all regression models. MBP was an additional covariable for the following carotid geometry (IMT, Dext, Dint and WLR) and arterial stiffness (Stiffness, Dist, Comp and WCSA) parameters.

(DOCX)

## Acknowledgments

We thank all patients who participate in this study in all the participating centers. We also thank clinicians who participated in the recruitment of patients. From the HEGP: Laurence Amar, Arshid Azarine, Guillaume Bobrie, Antoine Chédid, Béatrice Fiquet, Michael Frank, Agnes Labatide-Alanore, Olivier Steinchen. We also acknowledge the contribution to clinical data acquisition by Alessandra Giavarini, University of Milano and Sebastien Savard, Hôtel-Dieu de Québec Hospital. From the ARCADIA register: Jean-Philippe Baguet, Louis Boyer, Béatrice Duly-Bouhannick, Philippe Gosse, Claire Mounier-Vehier, Alexandre Persu, François Sihol, Frédéric Thony and Emmanuel Touzé. We acknowledge Susan Blackburn and James C. Stanley, University of Michigan and Stanley Hazen, Cleveland Clinic and Carin Smith, Mayo Clinic, in the recruitment of USA patients and controls. We thank Louise Marais, Inserm UMR970, Paris, for the acquisition of echo-tracking data, Patrick Lacolley, Inserm UMR_S 1116, Nancy and François Cambien, Inserm UMR_S 1166, Paris who contributed to the genetic investigation in PPS3 and Mark Lathrop, Genome Quebec who contributed to the generation of genotypes in SU.VI.MAX. We thank Xavier Loyer, Siham Bennis, Lama Cheikh-Hussein and Giulia Hardouin for assistance performing *in silico* analyses and expression experiments.

## Author Contributions

**Conceived and designed the experiments:** MA PFP PTE DJM SKG PB JCK XJe NBN.

**Performed the experiments:** SRK NRT LJCV AK VD CT DF JA ZY MLY KH CB PB NBN.

**Analyzed the data:** SRK NRT LJCV EA CB PB NBN.

**Contributed reagents/materials/analysis tools:** DC PG JPE XJo PTE DJM SKG JCK XJe NBN.

**Wrote the paper:** SRK NRT LJCV DJM SKG JCK NBN.

**Examined and followed-up patients and established the diagnosis:** JA DKD ESYK PG JPE XJo APGR JWO HLG MA PFP IJK SKG JCK XJe.

## References

1. Olin JW, Sealove BA (2011) Diagnosis, management, and future developments of fibromuscular dysplasia. J Vasc Surg 53: 826-836 e821. doi: 10.1016/j.jvs.2010.10.066 PMID: 21236620

2. Plouin PF, Perdu J, La Batide-Alanore A, Boutouyrie P, Gimenez-Roqueplo AP, et al. (2007) Fibromuscular dysplasia. Orphanet J Rare Dis 2: 28. doi: 10.1186/1750-1172-2-28 PMID: 17555581

3. Olin JW, Froehlich J, Gu X, Bacharach JM, Eagle K, et al. (2012) The United States Registry for Fibromuscular Dysplasia: results in the first 447 patients. Circulation 125: 3182-3190. doi: 10.1161/CIRCULATIONAHA.112.091223 PMID: 22615343

4. Cragg AH, Smith TP, Thompson BH, Maroney TP, Stanson AW, et al. (1989) Incidental fibromuscular dysplasia in potential renal donors: long-term clinical follow-up. Radiology 172: 145-147. doi: 10.1148/radiology.172.1.2662248 PMID: 2662248

5. Hendricks NJ, Matsumoto AH, Angle JF, Baheti A, Sabri SS, et al. (2014) Is fibromuscular dysplasia underdiagnosed? A comparison of the prevalence of FMD seen in CORAL trial participants versus a
single institution population of renal donor candidates. Vasc Med 19: 363-367. doi: 10.1177/
1358863X14544715 PMID: 25082538

6. Persu A, Giavarini A, Touze E, Januszewicz A, Sapoval M, et al. (2014) European consensus on the diagnosis and management of fibromuscular dysplasia. J Hypertens 32: 1367-1378. doi: 10.1097/
HJH.0000000000000213 PMID: 24842696

7. Savard S, Steichen O, Azarine A, Azizi M, Jeunemaitre X, et al. (2012) Association between 2 angio-
graphic subtypes of renal artery fibromuscular dysplasia and clinical characteristics. Circulation 126:
3062–3069. doi: 10.1161/CIRCULATIONAHA.112.117499 PMID: 23155180

8. Pannier-Moreau I, Grimbert P, Fiquet-Kempf B, Vuagnat A, Jeunemaitre X, et al. (1997) Possible familial origin of multifocal renal artery fibromuscular dysplasia. J Hypertens 15: 1797-1801. PMID: 9488242

9. Kiando SR, Barlassina C, Cusi D, Galan P, Lathrop M, et al. (2015) Exome sequencing in seven fami-
lies and gene-based association studies indicate genetic heterogeneity and suggest possible candi-
dates for fibromuscular dysplasia. J Hypertens 33: 1802-1810. doi: 10.1097/HJH.0000000000000625
PMID: 26147384

10. McKenzie GA, Oderich GS, Kawashima A, Misra S (2013) Renal artery fibromuscular dysplasia in 2,640 renal donor subjects: a CT angiography analysis. J Vasc Interv Radiol 24: 1477-1480. doi: 10.1016/j.jvir.2013.06.006 PMID: 23911200

11. Perdu J, Gimenez-Roqueplo AP, Boutouyrie P, Beaujour S, Laloux B, et al. (2006) Alpha1-antitrypsin gene polymorphisms are not associated with renal arterial fibromuscular dysplasia. J Hypertens 24: 705-710. doi: 10.1097/01.hjh.0000217853.97369.42 PMID: 16531799

12. Poloskey SL, Kim E, Sanghani R, Al-Quthami AH, Arscott P, et al. (2012) Low yield of genetic testing for known vascular connective tissue disorders in patients with fibromuscular dysplasia. Vasc Med 17: 371-378. doi: 10.1177/1358863X12459650 PMID: 23064905

13. Boutouyrie P, Gimenez-Roqueplo AP, Fine E, Laloux B, Fiquet-Kempf B, et al. (2003) Evidence for carotid and radial artery wall subclinical lesions in renal fibromuscular dysplasia. J Hypertens 21: 2287-2295. doi: 10.1097/01.hjh.0000098143.70956.c1 PMID: 14654749

14. Engelen L, Ferreira I, Stehouwer CD, Boutouyrie P, Laurent S, et al. (2013) Reference intervals for common carotid intima-media thickness measured with echotracking: relation with risk factors. Eur Heart J 34: 2368-2380. doi: 10.1093/eurheartj/ehs380 PMID: 23186808

15. Coronary Artery Disease Genetics C (2011) A genome-wide association study in Europeans and South Asians identifies five new loci for coronary artery disease. Nat Genet 43: 339-344. doi: 10.1038/ng.782 PMID: 21378988

16. Beaudoin M, Gupta RM, Won HH, Lo KS, Do R, et al. (2015) Myocardial Infarction-Associated SNP at 6p24 Interferes With MEF2 Binding and Associates With PHACTR1 Expression Levels in Human Coronary Arteries. Arterioscler Thromb Vasc Biol 35: 1472-1479. doi: 10.1161/ATVBAHA.115.305534 PMID: 25838425

17. Freilinger T, Anttila V, de Vries B, Malik R, Kallela M, et al. (2012) Genome-wide association analysis identifies susceptibility loci for migraine without aura. Nat Genet 44: 777-782. doi: 10.1038/ng.2307 PMID: 22683712

18. Anttila V, Winsvold BS, Gormley P, Kurth T, Bettella F, et al. (2013) Genome-wide meta-analysis iden-
tifies new susceptibility loci for migraine. Nat Genet 45: 912-917. doi: 10.1038/ng.2676 PMID:
23793025

19. Debette S, Kamatani Y, Metso TM, Kloss M, Chauhan G, et al. (2015) Common variation in PHACTR1 is associated with susceptibility to cervical artery dissection. Nat Genet 47: 78-83. doi: 10.1038/ng.3154 PMID: 25420145

20. Bejot Y, Aboa-Eboule C, Debette S, Pezzini A, Tatlisumak T, et al. (2014) Characteristics and out-
comes of patients with multiple cervical artery dissection. Stroke 45: 37-41. doi: 10.1161/
STROKEAHA.113.001654 PMID: 24326451

21. Saw J, Aymong E, Sedlak T, Buller CE, Starovoytov A, et al. (2014) Spontaneous coronary artery dis-
section: association with predisposing arteriopathies and precipitating stressors and cardiovascular
outcomes. Circ Cardiovasc Interv 7: 645-655. doi: 10.1161/CIRCINTERVENTIONS.114.001760
PMID: 25294399

22. Prasad M, Tweet MS, Hayes SN, Leng S, Liang JJ, et al. (2015) Prevalence of extracoronary vascular abnormalities and fibromuscular dysplasia in patients with spontaneous coronary artery dissection. Am J Cardiol 115: 1672-1677. doi: 10.1016/j.amjcard.2015.03.011 PMID: 25929580

23. Michelis KC, Olin JW, Kadian-Dodov D, d'Escamard V, Kovacic JC (2014) Coronary artery manifesta-
tions of fibromuscular dysplasia. J Am Coll Cardiol 64: 1033-1046. doi: 10.1016/j.jacc.2014.07.014
PMID: 25190240
24. Turner AW, McPherson R (2015) PHACTR1: Functional Clues Linking a Genome-Wide Association Study Locus to Coronary Artery Disease. Arterioscler Thromb Vasc Biol 35: 1293–1295. doi: 10.1161/ATVBAHA.115.305680 PMID: 25995042

25. Allain B, Jarray R, Borriello L, Leforban B, Dufour S, et al. (2012) Neuropilin-1 regulates a new VEGF-induced gene, Phactr-1, which controls tubulogenesis and modulates lamellipodial dynamics in human endothelial cells. Cell Signal 24: 214–223. doi: 10.1016/j.cellsig.2011.09.003 PMID: 21939755

26. Jarray R, Allain B, Borriello L, Biard D, Loukaci A, et al. (2011) Depletion of the novel protein PHACTR-1 from human endothelial cells abolishes tube formation and induces cell death receptor apoptosis. Biochimie 93: 1668–1675. doi: 10.1016/j.biochi.2011.07.010 PMID: 21798305

27. Jarray R, Pavoni S, Borriello L, Allain B, Lopez N, et al. (2015) Disruption of phactr-1 pathway triggers pro-inflammatory and pro-atherogenic factors: New insights in atherosclerosis development. Biochimie.

28. Consortium CAD (2015) A comprehensive 1000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet 47: 1121–1130. doi: 10.1038/ng.3396 PMID: 26343387

29. Reschen ME, Lin D, Chalisey A, Soilleux EJ, O’Callaghan CA (2016) Genetic and environmental risk factors for atherosclerosis regulate transcription of phosphatase and actin regulating gene PHACTR1. Atherosclerosis 250: 95–105. doi: 10.1016/j.atherosclerosis.2016.04.025 PMID: 27187934

30. Hercberg S, Galan P, Preziosi P, Bertrais S, Mennen L, et al. (2004) The SU.VI.MAX Study: a randomized, placebo-controlled trial of the health effects of antioxidant vitamins and minerals. Arch Intern Med 164: 2335–2342. doi: 10.1001/archinte.164.21.2335 PMID: 15557412

31. Empana JP, Bean K, Guibout C, Thomas F, Bingham A, et al. (2011) Paris Prospective Study III: a study of novel heart rate parameters, baroreflex sensitivity and risk of sudden death. Eur J Epidemiol 26: 887–892. doi: 10.1007/s10654-011-9618-x PMID: 21964902

32. Ye Z, Kalloo FS, Dalenberg AK, Kullo IJ (2013) An electronic medical record-linked biorepository to identify novel biomarkers for atherosclerotic cardiovascular disease. Glob Cardiol Sci Pract 2013: 82–90. doi: 10.5339/gcsp.2013.10 PMID: 24689004

33. Tang WH, Hazen SL (2014) The contributory role of gut microbiota in cardiovascular disease. J Clin Invest 124: 4204–4211. doi: 10.1172/JCI72331 PMID: 25271725

34. Proust C, Empana JP, Boutouyrie P, Alivon M, Challande P, et al. (2015) Contribution of Rare and Common Genetic Variants to Plasma Lipid Levels and Carotid Stiffness and Geometry: A Substudy of the Paris Prospective Study 3. Circ Cardiovasc Genet 8: 628–636. doi: 10.1161/CIRCGENETICS. 114.000979 PMID: 26160806

35. Guo Y, He J, Zhao S, Wu H, Zhong X, et al. (2014) Illumina human exome genotyping array clustering and quality control. Nat Protoc 9: 2643–2662. doi: 10.1038/nprot.2014.174 PMID: 25321409

36. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575. doi: 10.1086/519795 PMID: 17701901

37. Patterson N, Price AL, Reich D (2006) Population structure and eigenanalysis. PLoS Genet 2: e190. doi: 10.1371/journal.pgen.0020190 PMID: 17194218

38. Meinders JM, Brands PJ, Willigers JM, Kornet L, Hoeks AP (2001) Assessment of the spatial homogeneity of artery dimension parameters with high frame rate 2-D B-mode. Ultrasound Med Biol 27: 785–794. PMID: 11516538

39. Van Bortel LM, Balkestein EJ, van der Heijden-Spek JJ, Vanmolkot FH, Staessen JA, et al. (2001) Non-invasive assessment of local arterial pulse pressure: comparison of applanation tonometry and echo-tracking. J Hypertens 19: 1037–1044. PMID: 11403351

40. Willer CJ, Li Y, Abecasis GR (2010) METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26: 2190–2191. doi: 10.1093/bioinformatics/btq340 PMID: 20616382

41. Lee S, Emond MJ, Bamshad MJ, Barnes KC, Rieder MJ, et al. (2012) Optimal unified approach for rare-variant association testing with application to small-sample case-control whole-exome sequencing studies. Am J Hum Genet 91: 224–237. doi: 10.1016/j.ajhg.2012.06.007 PMID: 22863193

42. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, et al. (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55: 611–622. doi: 10.1373/clinchem.2008.112797 PMID: 19246619

43. Schmittgen TD, Livak KJ (2008) Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc 3: 1101–1108. PMID: 18546601

44. Dina C, Bouatia-Naji N, Tucker N, Delling FN, Toomer K, et al. (2015) Genetic association analyses highlight biological pathways underlying mitral valve prolapse. Nat Genet 47: 1206–1211. doi: 10.1038/ng.3383 PMID: 26301497