Next Article in Journal
“I Just Want You to Hear That Term”: Characterizing Language Used in Fetal Cardiology Consultations
Previous Article in Journal
Shared Decision Making and Cardioneuroablation Allow Discontinuation of Permanent Pacing in Patients with Vagally Mediated Bradycardia
Previous Article in Special Issue
Influence of Chromosome 9p21.3 rs1333049 Variant on Telomere Length and Their Interactive Impact on the Prognosis of Coronary Artery Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JCDD
Volume 10
Issue 9
10.3390/jcdd10090393
jcdd-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways

by
Tamiel N. Turley
1,2,
Jeanne L. Theis
2,
Jared M. Evans
3,
Zachary C. Fogarty
3,
Rajiv Gulati
4,
Sharonne N. Hayes
4,
Marysia S. Tweet
4 and
Timothy M. Olson
2,4,5,*
1
Molecular Pharmacology and Experimental Therapeutics Track, Mayo Clinic Graduate School of Biomedical Sciences, Mayo Clinic, Rochester, MN 55905, USA
2
Cardiovascular Genetics Research Laboratory, Mayo Clinic, Rochester, MN 55905, USA
3
Department of Quantitative Health Sciences, Division of Computational Biology, Mayo Clinic, Rochester, MN 55905, USA
4
Department of Cardiovascular Medicine, Mayo Clinic, Rochester, MN 55905, USA
5
Department of Pediatric and Adolescent Medicine, Division of Pediatric Cardiology, Mayo Clinic, Rochester, MN 55905, USA
*
Author to whom correspondence should be addressed.
J. Cardiovasc. Dev. Dis. 2023, 10(9), 393; https://doi.org/10.3390/jcdd10090393
Submission received: 11 August 2023 / Accepted: 8 September 2023 / Published: 12 September 2023
(This article belongs to the Special Issue 10th Anniversary of JCDD—'Genetics' Section)
Browse Figures
Versions Notes

Abstract

:
Rare familial spontaneous coronary artery dissection (SCAD) kindreds implicate genetic disease predisposition and provide a unique opportunity for candidate gene discovery. Whole-genome sequencing was performed in fifteen probands with non-syndromic SCAD who had a relative with SCAD, eight of whom had a second relative with extra-coronary arteriopathy. Co-segregating variants and associated genes were prioritized by quantitative variant, gene, and disease-level metrics. Curated public databases were queried for functional relationships among encoded proteins. Fifty-four heterozygous coding variants in thirteen families co-segregated with disease and fulfilled primary filters of rarity, gene variation constraint, and predicted-deleterious protein effect. Secondary filters yielded 11 prioritized candidate genes in 12 families, with high arterial tissue expression (n = 7), high-confidence protein-level interactions with genes associated with SCAD previously (n = 10), and/or previous associations with connective tissue disorders and aortopathies (n = 3) or other vascular phenotypes in mice or humans (n = 11). High-confidence associations were identified among 10 familial SCAD candidate-gene-encoded proteins. A collagen-encoding gene was identified in five families, two with distinct variants in COL4A2. Familial SCAD is genetically heterogeneous, yet perturbations of extracellular matrix, cytoskeletal, and cell–cell adhesion proteins implicate common disease-susceptibility pathways. Incomplete penetrance and variable expression suggest genetic or environmental modifiers.

1. Introduction

Spontaneous coronary artery dissection (SCAD) has gained recognition over the past decade as an important cause of acute coronary syndrome, myocardial infarction, and sudden death in predominantly young to middle-aged women with few atherosclerotic risk factors [1]. Complete or partial occlusion of the coronary vessel lumen by an expanding medial hematoma and/or an intimal tear are the diagnostic features of SCAD on coronary angiography. Coexistence of fibromuscular dysplasia (FMD) in >50% of individuals with SCAD [1], identification of a shared genetic risk allele for FMD and SCAD [2], and association with non-FMD extra-coronary aneurysms or dissections suggest that SCAD is a clinical manifestation of a systemic arteriopathy. Although a majority of cases are sporadic, the lack of clinical risk factors for vascular disease together with reports of rare familial associations implicate a genetic predisposition to SCAD [3]. In a Mayo Clinic cohort, whole-genome sequencing (WGS) and genome-wide associations studies (GWAS) previously revealed a rare co-segregating variant in TLN1 with major effect in a familial case, and common risk-conferring variants in sporadic cases, respectively [4,5,6]. In a large cohort from the United Kingdom, WGS and adjudication of rare variants by American College of Medical Genetics (ACMG) criteria revealed seven genes that harbored pathogenic or likely-pathogenic variants, yet they represented only 3.6% of cases [7]. Indeed, SCAD is becoming increasing recognized as a genetically heterogeneous and complex disorder, necessitating comprehensive and multidisciplinary approaches for discovery and validation of susceptibility genes [8]. In the current study, WGS in 15 families provided a unique opportunity to identify rare co-segregating variants and candidate genes to advance understanding of SCAD pathogenesis.

2. Materials and Methods

2.1. Study Subjects

Probands from each family were consecutively enrolled in the Mayo Clinic SCAD registry after diagnostic confirmation of SCAD by review of coronary angiograms and abstraction of detailed demographic and clinical data from questionnaires and medical records from 30 August 2011 to 30 June 2020 [9]. The diagnostic criteria for SCAD were as previously described [5,10]. First-degree and more distant relatives with a clinical diagnosis of SCAD or other arteriopathy, which included known FMD, arterial aneurysm or dissection, and cerebral infarctions, were invited to participate in the study.

2.2. Next-Generation Sequencing and Bioinformatics Analysis

WGS and whole-exome sequencing (WES) were performed by the Mayo Clinic Medical Genome Facility, and variant call format (VCF) files were created by the Bioinformatics Core. For WGS of 15 probands and 22 family members, 150 base pair (bp) paired-end sequencing was performed on either the Illumina HiSeq4000 or NovaSeq 6000 (Illumina, San Diego, CA, USA). For 10 family members, whole-exome 100 bp or 150 bp paired-end sequencing was performed on either the Illumina HiSeq2000 or 4000 platforms, respectively, using the Agilent SureSelect Human All Exon capture kits (V4 + UTR or V5 + UTR; Agilent technologies, Santa Clara, CA, USA).
WGS and WES reads were aligned to the hg38 or hg19 reference genome with either BWA-MEM [11] or Novoalign [12], respectively. Alignments were followed by the sorting and marking of duplicate reads using Picard [13]. Local realignment of insertions/deletions (INDELs) and base quality score recalibration were performed using the Genome Analysis Toolkit version 3.4 [14]. Single-nucleotide variants (SNVs) and INDELs were called across all samples simultaneously using the Genome Analysis Toolkit’s UnifiedGenotyper [15] (WES reads) or GATK HaplotypeCaller [16] and GenotypeGVCFs [14] (WGS reads), followed by variant quality score recalibration. WGS yielded on average 1.1 billion 150 base paired-end reads, and WES yielded 101 million 100- and 150-base paired-end reads.

2.3. Genetic Variant Filtering and Prioritization

VCF files with SNV and INDEL calls from WGS and WES were uploaded into Ingenuity® Variant Analysis™ (IVA; QIAGEN, Redwood City, CA, USA). To retain only high-confidence variants, a quality control filter was applied that required a base call quality of at least 20, genotype quality of ≥30, and a read depth of ≥10. Variants within a simple repeat region, defined as 1–50 bp in unit length with ≥5 repetitions, were excluded using the Tandem Repeats Finder [17] algorithm.
Annotated variants that passed quality control were prioritized by filtering for rarity, evidence of co-segregation with disease, impact on protein, tolerance of the corresponding gene to variation, mRNA expression pattern in arterial tissues, and gene–disease association. Variants with a minor allele frequency (MAF) of <0.1% from the Genome Aggregation Database [18] (gnomAD) v2.1 and v3 across all races were included. Variants with an MAF of <1% were used for recessive inheritance modeling. Co-segregation with disease required that a variant was present in all affected family members and obligate carriers. Except for parents of SCAD probands, unaffected family members were not included in the study. Co-segregation modeling was based on autosomal dominant transmission of a heterozygous variant with incomplete penetrance and variable expression and/or inheritance of homozygous or compound heterozygous recessive alleles for the affected sibling cases.
To predict the consequence of SNVs and INDELs, coding variants were annotated utilizing the Ensembl Variant Effect Predictor [19] tool. Variants classified as missense, frameshift, start-loss, stop-gain, stop-loss, and canonical splice site loss variants were included. The predicted deleterious effect of missense and canonical splice site variants was evaluated using a Combined Annotation Dependent Depletion [20] (CADD) score, which was then indexed as a percentile of genome-wide variants according to variant type. CADD percentiles were computed using a custom script that generated a histogram of CADD PHRED scores for all variants with the same variant consequence. Missense and canonical splice site variants each with a CADD score of ≤24.5, corresponding to the lower quartile, were excluded. Conservation of substituted amino acids was determined from the NCBI HomoloGene database and multialigned with MUSCLE version 3.6 [21]. The IVA software platform was used to identify and prioritize variants that were classified as deleterious by the Human Gene Mutation Database (HGMD) [22].
Non-coding variants were predicted to impact gene regulation based on location within a microRNA (miRNA) binding/coding site or a transcription factor binding site (TFBS) with a Position Weight Matrix score >0.75 from the Factorbook database [23]. Variants were prioritized by filtering for CADD score ≥7 corresponding to the upper quartile specific to regulatory variants, and rank score of ≤3 in the RegulomeDB V2.0 database [24]. To eliminate false positives, in silico visual assessment of all variants was performed using the Integrative Genomics Viewer [25].

2.4. Candidate Gene Prioritization

Publicly available databases and in silico tools were used to predict the biological impact of each identified variant. For a coding variant, the gene’s intolerance to variation was determined by metrics in gnomAD [18]. A probability of loss-of-function (pLI) and/or an observed over expected (o/e) score applied to stop gain/loss and frameshift variants, and a Z score applied to missense and canonical splice site variants. All candidate genes were required to have a pLI score ≥0.9 (range 0.0–1.0) and/or an o/e of <0.35, or a Z score ≥ 2. Regulatory variants were required to be predicted as pathogenic with a score ≥0.5 according to the Functional Analysis through Hidden Markov Models with eXtended Features (FATHMM-XF) V2.3 database [26].
Secondary filters to prioritize genes that harbored rare, predicted-deleterious, co-segregating variants were based on high arterial tissue expression ranking, a candidate gene shared by two or more familial SCAD (F-SCAD) cases, known or predicted functional associations with SCAD candidate genes previously identified in our cohort (TLN1, EDN1, PHACTR1, ADAMTSL4, C1orf54, ECM1, FBN1, LINC00310, LRP1, MRPS21, AFAP1) [2,4,5,6], an established gene for monogenic connective tissue disorders or aortopathies, or a gene associated with vascular disease risk. High arterial tissue messenger ribonucleic acid (mRNA) expression was defined by 1–5 ranking for one or more arterial tissues (aorta, coronary, or tibial), among 54 distinct tissue types in the Genotype-Tissue Expression (GTEx) database [27]. Protein–protein interaction networks were generated with the STRING v11 database [28] and visualized with Cytoscape v3.8.2 software [29]. The commercially available Heritable Disorders of Connective Tissue Panel (57 genes (2019); GeneDx, Gaithersburg, MD, USA) and Aortopathy Comprehensive Panel (27 genes (2019); Invitae, San Francisco, CA, USA) were used to identify known disease genes. A link to a vascular disease phenotype was identified based upon (1) a genetic mouse model reported in Mouse Genome Informatics (MGI) [30], (2) an established genetic link to vascular disease in humans reported in HGMD [22], or (3) a vascular disease-associated trait annotated in the GWAS Catalog [31]. The secondary filtering process was non-iterative, and top candidate gene(s) selection for each family was based on fulfilling the highest number of secondary filtering parameters.

2.5. In Silico Analysis of Shared Biological Networks and Processes

A network propagation analysis was performed on 11 prioritized candidate genes using STRING [28] v11 with a high-confidence threshold of 0.7 and visualized with Cytoscape [29] v3.8.2 to identify shared predicted and known functional associations among F-SCAD candidate genes. The network was built by inputting the previously reported F-SCAD gene TLN1 [4], known TLN1-interacting genes (18 alpha-integrins, 9 beta-integrins, and 7 actin isoforms), and the prioritized candidate genes. Interactions included both direct (physical) and indirect (functional) associations with high-confidence thresholds of >0.7. Enriched functional pathways and gene ontologies were investigated using biological molecular databases that included the Kyto Encyclopedia of Genes and Genomes [32] (KEGG), Reactome [33], National Center for Biotechnology Information BioSystems [34], and GeneALaCart meta-analysis tool from GeneCards [35].

3. Results

3.1. Clinical Characteristics

Clinical data and DNA samples were obtained from fifteen probands with SCAD who had a relative with SCAD, eight of whom had a second relative with extra-coronary arteriopathy (Figure 1, Table 1). Families were comprised of two affected mother–daughter pairs (SCAD-12, 14), two affected sister pair–parent trios (SCAD-13, 15), one affected sibling trio (SCAD-06), six affected sister pairs (SCAD-01, 02, 03, 04, 05, 07), one affected sibling pair–daughter trio (SCAD-10), and three families with distantly-related affected relatives (SCAD-08, 09, 11).
Among the 15 familial SCAD (F-SCAD) kindreds, SCAD was confirmed by review of angiograms or medical records (n = 28) or reported by history (n = 2). Probands and relatives with SCAD had a mean age at first event of 43.4 ± 7.2 years (range 32 to 59 years), with the left anterior descending coronary artery branch most affected. Associated phenotypes included FMD of the renal, iliac, splenic, or carotid arteries in twelve individuals who were screened; coronary artery tortuosity in five; peripartum status in three; and migraine headaches in five. Probands in three families had recurrent SCAD in a location separate from the initial event, occurring 2–11 years after the previous dissection. Recent episodes of emotional stress were reported in seven individuals and physical exertion in three individuals as potential triggers for SCAD. Similar to findings in previous studies of SCAD, typical risk factors for atherosclerotic coronary artery disease were less frequent in F-SCAD. Only five individuals had hypertension, nine had hyperlipidemia, eight had active or prior tobacco use, and one had diabetes mellitus. Overall, demographic and clinical characteristics of our F-SCAD cases were indistinguishable from our previously reported SCAD cohort (484 individuals) [5], with the exception of the proportion with migraine headaches (36% vs 15%; p = 0.019) and hypertension (32% vs 15%; p = 0.038). A majority of patients with SCAD or another arteriopathy were women (34 of 38; 89%), all of white European descent, with a mean age at diagnosis of 44.9 ± 10.1 years (range 25 to 71 years). The phenotypes of relatives with other arteriopathies (n = 8) included abdominal, iliac, and/or cerebral aneurysms; carotid or superior mesenteric artery dissections; coronary tortuosity; and stroke.

3.2. Whole-Genome Sequencing Reveals Disease-Associated Candidate Genes in Familial SCAD

WGS was performed on genomic DNA samples from affected individuals (n = 31) in each of the 15 families (Figure 1). Unaffected parents underwent WGS or WES (n = 16) to determine transmission of identified candidate variants. All genome and exome sequencing passed standard quality control metrics for coverage. Variant call format files were uploaded and analyzed for each individual with Ingenuity Variant Analysis software.
For each of the 15 families, over six million variants in 55,000 genes were identified. A systematic filtering approach based on quantitative metrics was applied to identify the most plausible candidate gene(s) for F-SCAD in 15 kindreds (Figure 1 and Figure 2). The primary variant-level filter yielded 54 candidate variants among families, which were further scrutinized by a secondary gene and disease-level filter. The 12 identified variants within each of the 11 prioritized candidate genes were heterozygous nucleotide substitutions or deletions that resulted in amino acid substitutions or frameshifts, all within highly conserved regions (Figure 3). All variants fulfilled quality control metrics, were rare in the general population with a minor allele frequency of <0.1% (gnomAD), co-segregated with disease; and were predicted deleterious based on CADD score, loss of function variant, and/or HGMD link (Table 2).
At the gene level, candidates were required to be intolerant to variation with a probability of loss-of-function score >0.9 and/or an observed over expected score of >0.35 for truncating variants, or a Z-score of >2 for missense variants. For each family, 0–9 candidate genes fulfilled all primary variant and gene-constraint metrics. Among the six families with two or more affected siblings, no homozygous or compound heterozygous variants were identified that passed the primary filters, excluding autosomal recessive inheritance. All families were thus considered to have autosomal dominant inheritance of a heterozygous variant with incomplete penetrance and, in families with extra-coronary arteriopathy, variable expression. A secondary filter of seven additional gene- and disease-level metrics was applied to identify the top candidate gene(s) in each family. Seven genes ranked highly in arterial tissue expression (rank 1–5) among 54 distinct tissues, none of which showed differential expression in females (Table 2, Figure 4).
Our previous genomic studies of familial and sporadic SCAD revealed TLN1, FBN1, LRP1, PHACTR1, ECM1, ADAMTSL4, LINC00310, C1orf54, MRPS21, EDN1, and AFAP1 as candidate genes for disease susceptibility [4,5,6]. Accordingly, potential protein–protein interactions of these 11 genes with the 11 F-SCAD candidate genes identified in the current study were investigated by STRING, revealing 10 high-confidence first- or second-degree functional associations. Four high-confidence first-degree interactions were identified: IQGAP1–TLN1, COL3A1–FBN1, COL5A2–FBN1, and LRP1–LRP2. Two of the previously identified SCAD candidate genes—TLN1 and FBN1—each had eight high-confidence first- and/or second-degree interactions (Figure 5, Table 2).
To further determine a gene’s candidacy, a biological context filter was applied to identify genes that were previously associated with aortopathies and connective tissue disorders (CTD). This comprehensive analysis revealed three genes that overlapped both CTD and aortopathy gene panels (COL3A1, COL5A2, NOTCH1; Table 2). Candidate genes were scrutinized for associations with reported murine or human vascular abnormalities in the Mouse Genome Informatics and HGMD databases or a GWAS vascular trait in the GWAS Catalog. From this analysis, 11 genes with known vascular phenotype-gene associations were identified. Overall, 12 heterozygous nucleotide substitutions or deletions that resulted in amino acid substitutions or frameshifts within highly conserved regions of 11 genes passed the systematic filtering metrics.

3.3. Potential Associations among Familial SCAD Candidate Genes

For five of the twelve families, the top candidate was a collagen-encoding gene, each of which fulfilled five or more secondary filtering parameters. Moreover, two families had distinct heterozygous missense variants in COL4A2, which encodes the collagen type IV alpha 2 chain (SCAD-07 and SCAD-14; Figure 1, Table 2). To identify additional associations among F-SCAD candidate genes, interactions at the protein level were investigated by creating a protein–protein network of TLN1 and the 11 other candidate gene proteins in STRING. Among the 11 top-ranked F-SCAD candidates, 10 genes had high-confidence, first- or second-degree associations with another F-SCAD-candidate-gene-encoded protein (Figure 6).
Genetic perturbation of the actin cytoskeleton has been proposed as one pathway in SCAD pathogenesis [4,5,6]. Accordingly, shared pathways and biological processes of the actin cytoskeleton and blood vessel formation were investigated using curated pathway and Gene Ontology databases. Among the 11 candidate genes, 10 genes shared pathways and processes involved in the cytoskeleton structure, remodeling and signaling, and/or angiogenic processes (Figure 7).

3.4. Ancillary Investigation of Non-Coding Variants

Our comprehensive and stringent variant filtering and gene prioritization scheme identified a top-ranked candidate gene in 12 of the 15 families in analyses confined to variants in coding regions that altered protein structure. To explore the possibility that rare non-coding variants predicted to alter gene expression could co-segregate with SCAD, an ancillary investigation of defined regulatory domains was conducted in the three families in which a candidate gene was not identified. Following the same primary analysis used for coding variants, implementation of secondary filters identified rare heterozygous variants within a microRNA binding site in PDGFRA (SCAD-08) and transcription factor binding site in ATF4 (SCAD-10) (Figures S1 and S2, Table S1). Both candidates demonstrated second-degree protein interactions with both previously identified SCAD-associated genes and other F-SCAD candidate genes with coding variants identified in the current study (Figures S3 and S4).

4. Discussion

SCAD is becoming increasingly recognized as a complex, genetically heterogeneous disorder [8]. To discover susceptibility genes for SCAD and potential shared disease pathways, we performed WGS in 15 rare F-SCAD cases and applied a systematic filtering approach that utilized quantitative metrics at variant, gene, and disease levels. Consistent with the low yield of pathogenic variants identified in a large sporadic SCAD cohort [7], only two top-ranked variants in NOTCH1 and KCNK3 in the current study were classified as pathogenic by ACMG criteria. However, co-segregation with SCAD and extra-coronary arteriopathy proved to be a powerful prioritization filter, reducing the number of rare, predicted-deleterious variants by 36–96% among families. Stringent gene constraint thresholds were used to identify variants with the highest likelihood of perturbing protein function, rendering short lists of 0–9 candidate genes in each family. Further prioritization of these genes was based on arterial tissue expression, interaction with known SCAD-associated genes, and genetic association with murine and human vascular disorders. The 12 identified heterozygous nucleotide substitutions and deletions exhibited autosomal dominant inheritance with incomplete penetrance and variable expression, including in the six families with 2–3 affected siblings for which no variants fitting an autosomal recessive mode of inheritance were identified. Collagen-encoding genes were the most common class among the 11 top-ranked candidate genes in our F-SCAD cohort, a finding consistent with an emerging body of evidence from WES and WGS in sporadic SCAD cohorts [36,37,38].
Three of the top-ranked F-SCAD candidate genes (COL3A1, COL5A2, and NOTCH1) are disease genes for inherited aortopathy and CTDs. While FMD is the most common trait associated with SCAD, 5.1–8.2% of sporadic SCAD cases have been found to carry deleterious variants within genes associated with heritable CTDs and aortopathy, including vascular Ehlers–Danlos syndrome (COL3A1), Marfan syndrome (FBN1), and Loeys–Dietz syndrome (SMAD3) [7,39,40]. Accordingly, familial variants within these well-established arteriopathy and connective tissue disorder genes may underlie non-syndromic SCAD susceptibility risk by increasing arterial fragility.
In addition to rare co-segregating variants, GWAS have implicated common risk-conferring variants and FBN1, LRP1, PHACTR1, ECM1, ADAMTSL4, LINC00310, C1orf54, MRPS21, and AFAP1 as candidate genes for SCAD [4,5,6], highlighting the heterogeneity of the disorder. Notwithstanding, we previously identified high-confidence second-degree protein–protein interactions between TLN1 and PHACTR1, FBN1 and LRP1 [5]. Moreover, among the 11 F-SCAD candidate-gene-encoded proteins identified in the current study, 10 had direct or indirect functional associations with TLN1 and/or common-risk gene-encoded proteins FBN1, LRP1, EDN1, and AFAP1 responsible for regulating cell adhesion, migration, and signaling pathways that orchestrate organization of the actin cytoskeleton [41,42,43,44,45]. These interactions suggest that both rare and common risk-conferring SCAD genes can have a shared role in regulation of the actin cytoskeleton, and variants within these genes may underlie SCAD pathogenesis.
Aberrant protein–protein interactions and shared pathways have been implicated in intracranial and thoracic aortic aneurysms and dissections [46,47]. In one family (SCAD-01), Val3698Glu substitution was identified in LRP2, a gene involved in coronary artery morphogenesis. The LRP2-encoded protein has a secondary interaction with two familial gene-encoded proteins (NOTCH1 and CTNNB1). LRP2, a major receptor of lipocalin-2, promotes angiogenesis through iron- and reactive-oxygen-species-related pathways in rodent endothelial cells [48].
Ten F-SCAD candidate-gene-encoded proteins had interactions with TLN1 and/or another familial gene-encoded protein and shared processes involved in cell adhesion. The formation and disassembly of adhesions through the binding of a cell to the extracellular matrix or other cells is essential for maintaining and repairing tissue architecture. Aberrant changes in cell adhesion dynamics have played a crucial role in the pathobiology of thoracic aortic dissections and aneurysms [49,50]. Interactions between adhesion complexes and actin, microtubules and intermediate filaments, elements that make up cytoskeleton networks, profoundly influence cell shape, motility, and remodeling [51].

4.1. Cell–Cell Adhesion

A Ser701Pro substitution in CDH4 was identified in SCAD-06. CDH4 is important for retinal vascular development [52]. CDH4 belongs to the type I classical cadherins responsible for mediating calcium-dependent adheren junctions [53]. Cell–cell junctions allow the cell to respond to changes within their microenvironment through connections with the actin cytoskeleton [54]. Moreover, within endothelial cells, these junctional complexes are responsible for regulating blood vessel maintenance and angiogenesis [55]. In SCAD-12, a Asn287Ser substitution was identified in CTNNB1, which is primarily responsible for regulating and coordinating cell–cell adhesion by anchoring cadherins to the actin cytoskeleton and involved in Wnt signaling [56]. Inactivation of Ctnnb1 in murine endothelial cells leads to a decrease in cell–cell adhesions and alterations in vascular morphogenesis [57]. Moreover, in rat carotid arteries, CDH4 expression and CTNNB1 signaling were associated with increased cyclin D1 expression and vascular smooth muscle proliferation, implicating both in vascular disease [56]. IQGAP1 is involved in regulating the cytoskeletal architecture through interactions with actin and microtubules [58]. Within this gene, we identified an Asn1066Asp substitution in SCAD-04. Acting as a scaffolding protein, IQGAP1 interacts with both catenins and cadherins and regulates cadherin-mediated cell–cell adhesions [59]. Further, IQGAP1 directly interacts with active vascular endothelial growth factor receptor-2, essential in endothelial migration and proliferation contributing to angiogenesis [60].

4.2. Cell–Extracellular Matrix Adhesion

A particularly notable finding in our study was the identification of collagen genes as top candidates for F-SCAD in five families (SCAD-09, -03, -07, -14, -15), each harboring a rare missense variant leading to substitution of a highly conserved amino acid: COL3A1 (Arg446Cys), COL4A1 (Pro530Ser), COL4A2 (Gly151Ser and Gly413Arg), and COL5A2 (Gln902Pro). While only three of the family members affected with SCAD had undergone formal evaluation by a medical geneticist (SCAD-03 proband; SCAD-07 proband and her sister), none had been suspected of Marfan syndrome, Ehlers–Danlos syndrome, or other multi-system CTD. Moreover, brain imaging performed in all but one affected family member (SCAD-09 proband) revealed subtle vascular disease in only one individual with a tiny cerebral aneurysm (SCAD-07 proband), whose maternal aunt also had a cerebral aneurysm by history.
Collagens are a large family of extracellular matrix proteins responsible for a variety of functions, including tissue scaffolding, cell adhesion, and cell migration [61]. Collagens contain three polypeptide (α) chains that make up a triple helix, with a repeating Gly-X-Y triplet, in which every third residue is a glycine [61]. X and Y positions are frequently proline and 4-hydroxyproline residues that provide stability to the triple helix [61]. Each of our F-SCAD collagen variants were located within or near a triple helix region, and three were glycine or proline residue substitutions (Gly151Ser, Gly413Arg, and Pro530Ser). At the center of a triple helix where the α chains come together, a glycine residue is essential, as space precludes any other larger residue from fitting within the region [62]. Collagens have been associated with a wide spectrum of diseases. A large number of deleterious variants have been found in 22 collagen genes that range from early lethal phenotypes to common risk-conferring variants [62]. Included in heritable collagen disorders are Ehlers–Danlos syndrome (COL3A1, COL5A2) [63,64], brain small vessel disease-2 (COL4A2) [65], and arterial aneurysms (COL3A1, COL4A1) [66,67]. Of note, COL3A1 was one of seven genes with pathogenic variants identified in the first large-scale WGS study of sporadic SCAD [7]. Moreover, enrichment of rare variants in collagen genes has been reported subsequently in three independent WES and WGS studies of sporadic SCAD [36,37,38].

4.3. Adhesion Signaling

Cell adhesive interactions and networks that exist between both cell–cell and cell–matrix adhesions are controlled by signaling networks. Signaling networks are responsible for regulating downstream functions that include cell adhesion, migration, and mechanotransduction, thereby integrating signals from the outside and within the cell [68]. In EPHB4, we identified a Gly462Arg substitution in SCAD-05. EPHB4-mediated signaling, by activation from ephrinB2 ligand, is involved in regulating cell adhesion and migration [69]. Inactivation of Ephb4 in murine endothelial cells leads to rupturing of cardiac capillaries and cardiomyocyte hypertrophy [70]. An Arg1279Cys substitution in NOTCH1 was identified in SCAD-13, one of two instances in our study of a likely pathogenic variant based on ACMG criteria. Moreover, NOTCH1 fulfilled five secondary parameters in our filtering workflow. NOTCH1 receptor signaling is involved in a wide variety of processes, including regulating adheren junction assembly and vascular barrier function [71]. Deleterious variants within NOTCH1 have been associated with bicuspid aortic valve disease, thoracic aortic aneurysms, and pulmonary hypertension [72,73]. In a recent case report, a Arg1438Cys substitution was identified by WES in a patient with recurrent and pregnancy-associated SCAD [74].
An Arg261fs* deletion in KCNK3 was identified in SCAD-02, the second instance of a likely pathogenic variant. KCNK3 is responsible for regulating K+ conductance and resting membrane potential in pulmonary artery smooth muscle cells [75]. Deleterious variants within KCNK3 have been linked to pulmonary arterial hypertension [76].

4.4. Limitations

The co-segregation filter proved to be a powerful tool for generating relatively short lists of potential SCAD-susceptibility genes in each family. However, the families in this study were not powered for genome-wide linkage analysis; co-segregation with SCAD was inferred in two families (SCAD-09, SCAD-15), for which clinical and genetic data were lacking in the relative with SCAD, and our candidate-gene-ranking workflow imposed strict gene-constraint requirements and relied on existing bioinformatics databases, predictive metrics, and known protein–protein and gene–disease associations. Consequently, additional co-segregating variants could have a role in SCAD pathogenesis. Predicting the potential biological impact of non-coding region variants poses a particular challenge. Functional experiments and identification of additional rare coding or non-coding variants in top candidate genes within other SCAD cohorts would lend further support to our findings.

5. Conclusions

F-SCAD is a genetically heterogeneous disorder, yet collagen genes were top candidates in five families in our cohort. Moreover, several genes had high-confidence primary or secondary protein–protein interactions with both TLN1 and FBN1, and our study further implicates perturbation of structural proteins that function in actin cytoskeleton and cell–cell adhesion in disease pathogenesis. The well-established predilection of SCAD for women was also observed in most of our F-SCAD kindreds. However, our top-ranked candidate genes did not exhibit sex-specific differences in coronary artery mRNA expression, suggesting that risk allele effects may be accentuated by a generalized increase in vascular fragility conferred by female hormone milieu [77]. Indeed, incomplete penetrance and variable expression of rare, predicted-deleterious, co-segregating variants implicate genetic or environmental modifiers of this episodic disorder.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/jcdd10090393/s1. The following supporting information can be found in the Supplementary File, Table S1: Top candidate genes harboring non-coding variants in familial SCAD. Figure S1: Non-coding variant filtering and gene prioritization workflow identifies additional candidate genes for familial SCAD. Figure S2: Non-coding variants in familial SCAD alter transcription factor and microRNA binding sites [78,79]. Figure S3: Functional association between individual F-SCAD non-coding candidates and previously reported SCAD genes. Figure S4: Network interactions of functional associations among F-SCAD non-coding and coding candidates.

Author Contributions

Conceptualization, T.N.T. and T.M.O.; methodology, T.N.T., J.M.E. and Z.C.F.; formal analysis, T.N.T.; investigation, T.N.T. and T.M.O.; resources, M.S.T., S.N.H. and R.G.; data curation, T.N.T., J.M.E. and Z.C.F.; writing—original draft preparation, T.N.T. and T.M.O.; writing—review and editing, T.N.T., T.M.O., J.L.T., M.S.T. and S.N.H.; visualization, T.N.T.; supervision, T.M.O.; project administration, T.N.T. and T.M.O.; funding acquisition, T.M.O., M.S.T. and S.N.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funding from SCAD Research Inc., the Heller Foundation Inc., and National Institutes of Health (NIH) T32 GM72474.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Mayo Clinic (protocol 11–000160 approved 25 May 2011).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data from this study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank Rhianna Sundsbak for technical assistance, Brenda Speltz and Jill Boyum for patient recruitment, and the patients who volunteered for this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACMG, American College of Medical Genetics; CTD, connective tissue disorder; F-SCAD, familial SCAD; FMD, fibromuscular dysplasia; GWAS, genome-wide association study; HGMD, Human Gene Mutation Database; SCAD, spontaneous coronary artery dissection; WES, whole-exome sequencing; WGS, whole-genome sequencing.

References

  1. Hayes, S.N.; Tweet, M.S.; Adlam, D.; Kim, E.S.H.; Gulati, R.; Price, J.E.; Rose, C.H. Spontaneous coronary artery dissection: Jacc state-of-the-art review. J. Am. Coll. Cardiol. 2020, 76, 961–984. [Google Scholar] [CrossRef] [PubMed]
  2. Adlam, D.; Olson, T.M.; Combaret, N.; Kovacic, J.C.; Iismaa, S.E.; Al-Hussaini, A.; O’Byrne, M.M.; Bouajila, S.; Georges, A.; Mishra, K.; et al. Association of the PHACTR1/EDN1 genetic locus with spontaneous coronary artery dissection. J. Am. Coll. Cardiol. 2019, 73, 58–66. [Google Scholar] [CrossRef] [PubMed]
  3. Goel, K.; Tweet, M.; Olson, T.M.; Maleszewski, J.J.; Gulati, R.; Hayes, S.N. Familial spontaneous coronary artery dissection: Evidence for genetic susceptibility. JAMA Intern. Med. 2015, 175, 821–826. [Google Scholar] [CrossRef]
  4. Turley, T.N.; Theis, J.L.; Sundsbak, R.S.; Evans, J.M.; O’Byrne, M.M.; Gulati, R.; Tweet, M.S.; Hayes, S.N.; Olson, T.M. Rare missense variants in tln1 are associated with familial and sporadic spontaneous coronary artery dissection. Circ. Genom. Precis. Med. 2019, 12, e002437. [Google Scholar] [CrossRef] [PubMed]
  5. Turley, T.N.; O’Byrne, M.M.; Kosel, M.L.; de Andrade, M.; Gulati, R.; Hayes, S.N.; Tweet, M.S.; Olson, T.M. Identification of susceptibility loci for spontaneous coronary artery dissection. JAMA Cardiol. 2020, 5, 929–938. [Google Scholar] [CrossRef] [PubMed]
  6. Turley, T.N.; Kosel, M.L.; Bamlet, W.R.; Gulati, R.; Hayes, S.N.; Tweet, M.S.; Olson, T.M. Susceptibility locus for pregnancy-associated spontaneous coronary artery dissection. Circ. Genom. Precis. Med. 2021, 14, e003398. [Google Scholar] [CrossRef]
  7. Carss, K.J.; Baranowska, A.A.; Armisen, J.; Webb, T.R.; Hamby, S.E.; Premawardhana, D.; Al-Hussaini, A.; Wood, A.; Wang, Q.; Deevi, S.V.V.; et al. Spontaneous coronary artery dissection: Insights on rare genetic variation from genome sequencing. Circ. Genom. Precis. Med. 2020, 13, e003030. [Google Scholar] [CrossRef]
  8. Amrani-Midoun, A.; Adlam, D.; Bouatia-Naji, N. Recent advances on the genetics of spontaneous coronary artery dissection. Circ. Genomic Precis. Med. 2021, 14, e003393. [Google Scholar] [CrossRef]
  9. Hayes, S.N.; Kim, E.S.H.; Saw, J.; Adlam, D.; Arslanian-Engoren, C.; Economy, K.E.; Ganesh, S.K.; Gulati, R.; Lindsay, M.E.; Mieres, J.H.; et al. Spontaneous coronary artery dissection: Current state of the science: A scientific statement from the american heart association. Circulation 2018, 137, e523–e557. [Google Scholar] [CrossRef]
  10. Tweet, M.S.; Hayes, S.N.; Pitta, S.R.; Simari, R.D.; Lerman, A.; Lennon, R.J.; Gersh, B.J.; Khambatta, S.; Best, P.J.M.; Rihal, C.S.; et al. Clinical features, management, and prognosis of spontaneous coronary artery dissection. Circulation 2012, 126, 579–588. [Google Scholar] [CrossRef]
  11. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997. [Google Scholar] [CrossRef]
  12. Novocraft Technologies. NovoAlign. Available online: http://www.novocraft.com/ (accessed on 18 January 2021).
  13. Picard Tools—By Broad Institute. Available online: http://broadinstitute.github.io/picard/ (accessed on 6 October 2021).
  14. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The genome analysis toolkit: A mapreduce framework for analyzing next-generation dna sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed]
  15. DePristo, M.A.; Banks, E.; Poplin, R.; Garimella, K.V.; Maguire, J.R.; Hartl, C.; Philippakis, A.A.; del Angel, G.; Rivas, M.A.; Hanna, M.; et al. A framework for variation discovery and genotyping using next-generation dna sequencing data. Nat. Genet. 2011, 43, 491–498. [Google Scholar] [CrossRef] [PubMed]
  16. Poplin, R.; Ruano-Rubio, V.; DePristo, M.; Fennell, T.; Carneiro, M.; Van der Auwera, G.; Kling, D.; Gauthier, L.; Levy-Moonshine, A.; Roazen, D.; et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv 2018. [Google Scholar] [CrossRef]
  17. Benson, G. Tandem repeats finder: A program to analyze dna sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [PubMed]
  18. Karczewski, K.J.; Francioli, L.C.; Tiao, G.; Cummings, B.B.; Alföldi, J.; Wang, Q.; Collins, R.L.; Laricchia, K.M.; Ganna, A.; Birnbaum, D.P.; et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 2020, 581, 434–443. [Google Scholar] [CrossRef]
  19. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.S.; Thormann, A.; Flicek, P.; Cunningham, F. The ensembl variant effect predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef]
  20. Rentzsch, P.; Witten, D.; Cooper, G.M.; Shendure, J.; Kircher, M. CADD: Predicting the deleteriousness of variants throughout the human genome. Nucleic Acids Res. 2019, 47, D886–D894. [Google Scholar] [CrossRef]
  21. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef]
  22. Stenson, P.D.; Mort, M.; Ball, E.V.; Shaw, K.; Phillips, A.; Cooper, D.N. The human gene mutation database: Building a comprehensive mutation repository for clinical and molecular genetics, diagnostic testing and personalized genomic medicine. Hum. Genet. 2014, 133, 1–9. [Google Scholar] [CrossRef]
  23. Wang, J.; Zhuang, J.; Iyer, S.; Lin, X.-Y.; Greven, M.C.; Kim, B.-H.; Moore, J.; Pierce, B.G.; Dong, X.; Virgil, D.; et al. Factorbook.Org: A wiki-based database for transcription factor-binding data generated by the encode consortium. Nucleic Acids Res. 2013, 41, D171–D176. [Google Scholar] [CrossRef] [PubMed]
  24. Boyle, A.P.; Hong, E.L.; Hariharan, M.; Cheng, Y.; Schaub, M.A.; Kasowski, M.; Karczewski, K.J.; Park, J.; Hitz, B.C.; Weng, S.; et al. Annotation of functional variation in personal genomes using regulomeDB. Genome Res. 2012, 22, 1790–1797. [Google Scholar] [CrossRef] [PubMed]
  25. Thorvaldsdóttir, H.; Robinson, J.T.; Mesirov, J.P. Integrative genomics viewer (IGV): High-performance genomics data visualization and exploration. Brief. Bioinform. 2013, 14, 178–192. [Google Scholar] [CrossRef]
  26. Shihab, H.A.; Gough, J.; Cooper, D.N.; Stenson, P.D.; Barker, G.L.A.; Edwards, K.J.; Day, I.N.M.; Gaunt, T.R. Predicting the functional, molecular, and phenotypic consequences of amino acid substitutions using hidden markov models. Hum. Mutat. 2013, 34, 57–65. [Google Scholar] [CrossRef] [PubMed]
  27. GTEx Consortium human genomics. the genotype-tissue expression (gtex) pilot analysis: Multitissue gene regulation in humans. Science 2015, 348, 648–660. [Google Scholar] [CrossRef]
  28. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef]
  29. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  30. Smith, C.L.; Blake, J.A.; Kadin, J.A.; Richardson, J.E.; Bult, C.J. Mouse genome database group mouse genome database (mgd)-2018: Knowledgebase for the laboratory mouse. Nucleic Acids Res. 2018, 46, D836–D842. [Google Scholar] [CrossRef]
  31. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The NHGRI-EBI GWAS catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef]
  32. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  33. Jassal, B.; Matthews, L.; Viteri, G.; Gong, C.; Lorente, P.; Fabregat, A.; Sidiropoulos, K.; Cook, J.; Gillespie, M.; Haw, R.; et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2020, 48, D498–D503. [Google Scholar] [CrossRef] [PubMed]
  34. Geer, L.Y.; Marchler-Bauer, A.; Geer, R.C.; Han, L.; He, J.; He, S.; Liu, C.; Shi, W.; Bryant, S.H. The NCBI BioSystems database. Nucleic Acids Res. 2010, 38, D492–D496. [Google Scholar] [CrossRef] [PubMed]
  35. Harel, A.; Inger, A.; Stelzer, G.; Strichman-Almashanu, L.; Dalah, I.; Safran, M.; Lancet, D. GIFtS: Annotation landscape analysis with genecards. BMC Bioinform. 2009, 10, 348. [Google Scholar] [CrossRef] [PubMed]
  36. Zekavat, S.M.; Chou, E.L.; Zekavat, M.; Pampana, A.; Paruchuri, K.; Lino Cardenas, C.L.; Koyama, S.; Ghazzawi, Y.; Kii, E.; Uddin, M.M.; et al. Fibrillar collagen variants in spontaneous coronary artery dissection. JAMA Cardiol. 2022, 7, 396–406. [Google Scholar] [CrossRef] [PubMed]
  37. Tarr, I.; Hesselson, S.; Iismaa, S.E.; Rath, E.; Monger, S.; Troup, M.; Mishra, K.; Wong, C.M.Y.; Hsu, P.-C.; Junday, K.; et al. Exploring the genetic architecture of spontaneous coronary artery dissection using whole-genome sequencing. Circ. Genom. Precis. Med. 2022, 15, e003527. [Google Scholar] [CrossRef] [PubMed]
  38. Wang, Y.; Starovoytov, A.; Murad, A.M.; Hunker, K.L.; Brunham, L.R.; Li, J.Z.; Saw, J.; Ganesh, S.K. Burden of rare genetic variants in spontaneous coronary artery dissection with high-risk features. JAMA Cardiol. 2022, 7, 1045–1055. [Google Scholar] [CrossRef]
  39. Henkin, S.; Negrotto, S.M.; Tweet, M.S.; Kirmani, S.; Deyle, D.R.; Gulati, R.; Olson, T.M.; Hayes, S.N. Spontaneous coronary artery dissection and its association with heritable connective tissue disorders. Heart 2016, 102, 876–881. [Google Scholar] [CrossRef]
  40. Kaadan, M.I.; MacDonald, C.; Ponzini, F.; Duran, J.; Newell, K.; Pitler, L.; Lin, A.; Weinberg, I.; Wood, M.J.; Lindsay, M.E. Prospective cardiovascular genetics evaluation in spontaneous coronary artery dissection. Circ. Genom. Precis. Med. 2018, 11, e001933. [Google Scholar] [CrossRef]
  41. Gough, R.E.; Goult, B.T. The tale of two talins—Two isoforms to fine-tune integrin signalling. FEBS Lett. 2018, 592, 2108–2125. [Google Scholar] [CrossRef]
  42. Bax, D.V.; Bernard, S.E.; Lomas, A.; Morgan, A.; Humphries, J.; Shuttleworth, C.A.; Humphries, M.J.; Kielty, C.M. Cell adhesion to fibrillin-1 molecules and microfibrils is mediated by alpha 5 beta 1 and alpha v beta 3 integrins. J. Biol. Chem. 2003, 278, 34605–34616. [Google Scholar] [CrossRef]
  43. Wujak, L.; Böttcher, R.T.; Pak, O.; Frey, H.; El Agha, E.; Chen, Y.; Schmitt, S.; Bellusci, S.; Schaefer, L.; Weissmann, N.; et al. Low density lipoprotein receptor-related protein 1 couples β1 integrin activation to degradation. Cell. Mol. Life Sci. CMLS 2018, 75, 1671–1685. [Google Scholar] [CrossRef] [PubMed]
  44. Koyama, Y.; Baba, A. Endothelin-induced cytoskeletal actin re-organization in cultured astrocytes: Inhibition by c3 adp-ribosyltransferase. Glia 1996, 16, 342–350. [Google Scholar] [CrossRef]
  45. Cunnick, J.M.; Kim, S.; Hadsell, J.; Collins, S.; Cerra, C.; Reiser, P.; Flynn, D.C.; Cho, Y. Actin filament-associated protein 1 is required for csrc activity and secretory activation in the lactating mammary gland. Oncogene 2015, 34, 2640–2649. [Google Scholar] [CrossRef] [PubMed]
  46. Kim, J.H.; Na, C.-Y.; Choi, S.Y.; Kim, H.W.; Du Kim, Y.; Kwon, J.B.; Chung, M.Y.; Hong, J.M.; Park, C.B. Integration of gene-expression profiles and pathway analysis in ascending thoracic aortic aneurysms. Ann. Vasc. Surg. 2010, 24, 538–549. [Google Scholar] [CrossRef]
  47. Bi, S.; Liu, R.; He, L.; Li, J.; Gu, J. Bioinformatics analysis of common key genes and pathways of intracranial, abdominal, and thoracic aneurysms. BMC Cardiovasc. Disord. 2021, 21, 14. [Google Scholar] [CrossRef] [PubMed]
  48. Wu, L.; Du, Y.; Lok, J.; Lo, E.H.; Xing, C. Lipocalin-2 Enhances angiogenesis in rat brain endothelial cells via reactive oxygen species and iron-dependent mechanisms. J. Neurochem. 2015, 132, 622–628. [Google Scholar] [CrossRef]
  49. Wei, X.; Sun, Y.; Wu, Y.; Zhu, J.; Gao, B.; Yan, H.; Zhao, Z.; Zhou, J.; Jing, Z. Downregulation of talin-1 expression associates with increased proliferation and migration of vascular smooth muscle cells in aortic dissection. BMC Cardiovasc. Disord. 2017, 17, 162. [Google Scholar] [CrossRef]
  50. Guo, D.-C.; Pannu, H.; Tran-Fadulu, V.; Papke, C.L.; Yu, R.K.; Avidan, N.; Bourgeois, S.; Estrera, A.L.; Safi, H.J.; Sparks, E.; et al. Mutations in smooth muscle alpha-actin (ACTA2) lead to thoracic aortic aneurysms and dissections. Nat. Genet. 2007, 39, 1488–1493. [Google Scholar] [CrossRef]
  51. Weber, G.F.; Bjerke, M.A.; DeSimone, D.W. Integrins and cadherins join forces to form adhesive networks. J. Cell Sci. 2011, 124, 1183–1193. [Google Scholar] [CrossRef]
  52. Dorrell, M.I.; Aguilar, E.; Friedlander, M. Retinal Vascular development is mediated by endothelial filopodia, a preexisting astrocytic template and specific r-cadherin adhesion. Investig. Ophthalmol. Vis. Sci. 2002, 43, 3500–3510. [Google Scholar]
  53. Oda, H.; Takeichi, M. Evolution: Structural and functional diversity of cadherin at the adherens junction. J. Cell Biol. 2011, 193, 1137–1146. [Google Scholar] [CrossRef] [PubMed]
  54. Kemler, R. From cadherins to catenins: Cytoplasmic protein interactions and regulation of cell adhesion. Trends Genet. TIG 1993, 9, 317–321. [Google Scholar] [CrossRef] [PubMed]
  55. Lampugnani, M.G.; Dejana, E. Adherens junctions in endothelial cells regulate vessel maintenance and angiogenesis. Thromb. Res. 2007, 120 (Suppl. 2), S1–S6. [Google Scholar] [CrossRef] [PubMed]
  56. Slater, S.C.; Koutsouki, E.; Jackson, C.L.; Bush, R.C.; Angelini, G.D.; Newby, A.C.; George, S.J. R-Cadherin:Beta-Catenin complex and its association with vascular smooth muscle cell proliferation. Arterioscler. Thromb. Vasc. Biol. 2004, 24, 1204–1210. [Google Scholar] [CrossRef] [PubMed]
  57. Cattelino, A.; Liebner, S.; Gallini, R.; Zanetti, A.; Balconi, G.; Corsi, A.; Bianco, P.; Wolburg, H.; Moore, R.; Oreda, B.; et al. The conditional inactivation of the beta-catenin gene in endothelial cells causes a defective vascular pattern and increased vascular fragility. J. Cell Biol. 2003, 162, 1111–1122. [Google Scholar] [CrossRef] [PubMed]
  58. Fukata, M.; Watanabe, T.; Noritake, J.; Nakagawa, M.; Yamaga, M.; Kuroda, S.; Matsuura, Y.; Iwamatsu, A.; Perez, F.; Kaibuchi, K. Rac1 and Cdc42 capture microtubules through IQGAP1 and CLIP-170. Cell 2002, 109, 873–885. [Google Scholar] [CrossRef]
  59. Fukata, M.; Kuroda, S.; Nakagawa, M.; Kawajiri, A.; Itoh, N.; Shoji, I.; Matsuura, Y.; Yonehara, S.; Fujisawa, H.; Kikuchi, A.; et al. Cdc42 and Rac1 regulate the interaction of IQGAP1 with Beta-Catenin. J. Biol. Chem. 1999, 274, 26044–26050. [Google Scholar] [CrossRef]
  60. Yamaoka-Tojo, M.; Ushio-Fukai, M.; Hilenski, L.; Dikalov, S.I.; Chen, Y.E.; Tojo, T.; Fukai, T.; Fujimoto, M.; Patrushev, N.A.; Wang, N.; et al. IQGAP1, a novel vascular endothelial growth factor receptor binding protein, is involved in reactive oxygen species--dependent endothelial migration and proliferation. Circ. Res. 2004, 95, 276–283. [Google Scholar] [CrossRef]
  61. Kadler, K.E.; Baldock, C.; Bella, J.; Boot-Handford, R.P. Collagens at a glance. J. Cell Sci. 2007, 120, 1955–1958. [Google Scholar] [CrossRef]
  62. Myllyharju, J.; Kivirikko, K.I. Collagens and collagen-related diseases. Ann. Med. 2001, 33, 7–21. [Google Scholar] [CrossRef]
  63. Frank, M.; Albuisson, J.; Ranque, B.; Golmard, L.; Mazzella, J.-M.; Bal-Theoleyre, L.; Fauret, A.-L.; Mirault, T.; Denarié, N.; Mousseaux, E.; et al. The Type of Variants at the COL3A1 Gene Associates with the Phenotype and Severity of Vascular Ehlers-Danlos Syndrome. Eur. J. Hum. Genet. EJHG 2015, 23, 1657–1664. [Google Scholar] [CrossRef] [PubMed]
  64. Watanabe, M.; Nakagawa, R.; Naruto, T.; Kohmoto, T.; Suga, K.; Goji, A.; Kagami, S.; Masuda, K.; Imoto, I. A Novel Missense Mutation of COL5A2 in a patient with ehlers–danlos syndrome. Hum. Genome Var. 2016, 3, 16030. [Google Scholar] [CrossRef]
  65. Yoneda, Y.; Haginoya, K.; Arai, H.; Yamaoka, S.; Tsurusaki, Y.; Doi, H.; Miyake, N.; Yokochi, K.; Osaka, H.; Kato, M.; et al. De Novo and inherited mutations in COL4A2, encoding the Type IV collagen A2 chain cause porencephaly. Am. J. Hum. Genet. 2012, 90, 86–90. [Google Scholar] [CrossRef] [PubMed]
  66. van de Luijtgaarden, K.M.; Heijsman, D.; Maugeri, A.; Weiss, M.M.; Verhagen, H.J.M.; IJpma, A.; Brüggenwirth, H.T.; Majoor-Krakauer, D. First genetic analysis of aneurysm genes in familial and sporadic abdominal aortic aneurysm. Hum. Genet. 2015, 134, 881–893. [Google Scholar] [CrossRef] [PubMed]
  67. Plaisier, E.; Gribouval, O.; Alamowitch, S.; Mougenot, B.; Prost, C.; Verpont, M.C.; Marro, B.; Desmettre, T.; Cohen, S.Y.; Roullet, E.; et al. COL4A1 mutations and hereditary angiopathy, nephropathy, aneurysms, and muscle cramps. N. Engl. J. Med. 2007, 357, 2687–2695. [Google Scholar] [CrossRef]
  68. Huveneers, S.; Danen, E.H.J. Adhesion signaling—Crosstalk between integrins, Src and Rho. J. Cell Sci. 2009, 122, 1059–1069. [Google Scholar] [CrossRef]
  69. Sakamoto, H.; Zhang, X.-Q.; Suenobu, S.; Ohbo, K.; Ogawa, M.; Suda, T. Cell Adhesion to Ephrinb2 Is induced by EphB4 independently of its kinase activity. Biochem. Biophys. Res. Commun. 2004, 321, 681–687. [Google Scholar] [CrossRef]
  70. Luxán, G.; Stewen, J.; Díaz, N.; Kato, K.; Maney, S.K.; Aravamudhan, A.; Berkenfeld, F.; Nagelmann, N.; Drexler, H.C.; Zeuschner, D.; et al. Endothelial EphB4 maintains vascular integrity and transport function in adult heart. eLife 2019, 8, e45863. [Google Scholar] [CrossRef]
  71. Polacheck, W.J.; Kutys, M.L.; Yang, J.; Eyckmans, J.; Wu, Y.; Vasavada, H.; Hirschi, K.K.; Chen, C.S. A Non-canonical notch complex regulates adherens junctions and vascular barrier function. Nature 2017, 552, 258–262. [Google Scholar] [CrossRef]
  72. McKellar, S.H.; Tester, D.J.; Yagubyan, M.; Majumdar, R.; Ackerman, M.J.; Sundt, T.M. Novel NOTCH1 mutations in patients with bicuspid aortic valve disease and thoracic aortic aneurysms. J. Thorac. Cardiovasc. Surg. 2007, 134, 290–296. [Google Scholar] [CrossRef]
  73. Dabral, S.; Tian, X.; Kojonazarov, B.; Savai, R.; Ghofrani, H.A.; Weissmann, N.; Florio, M.; Sun, J.; Jonigk, D.; Maegel, L.; et al. Notch1 signalling regulates endothelial proliferation and apoptosis in pulmonary arterial hypertension. Eur. Respir. J. 2016, 48, 1137–1149. [Google Scholar] [CrossRef] [PubMed]
  74. Bai, B.; Zhang, M.; Zhuang, Y.; Zhu, J.; Li, W.; Ma, W.; Chen, H. The pregnancy-associated spontaneous coronary artery dissection in a young woman with a novel missense mutation in NOTCH1: A case report. BMC Med. Genet. 2020, 21, 119. [Google Scholar] [CrossRef] [PubMed]
  75. Olschewski, A.; Li, Y.; Tang, B.; Hanze, J.; Eul, B.; Bohle, R.M.; Wilhelm, J.; Morty, R.E.; Brau, M.E.; Weir, E.K.; et al. Impact of TASK-1 in human pulmonary artery smooth muscle cells. Circ. Res. 2006, 98, 1072–1080. [Google Scholar] [CrossRef] [PubMed]
  76. Navas Tejedor, P.; Tenorio Castaño, J.; Palomino Doza, J.; Arias Lajara, P.; Gordo Trujillo, G.; López Meseguer, M.; Román Broto, A.; Lapunzina Abadía, P.; Escribano Subía, P. An homozygous mutation in kcnk3 is associated with an aggressive form of hereditary pulmonary arterial hypertension. Clin. Genet. 2017, 91, 453–457. [Google Scholar] [CrossRef] [PubMed]
  77. Tweet, M.S.; Miller, V.M.; Hayes, S.N. The evidence on estrogen, progesterone, and spontaneous coronary artery dissection. JAMA Cardiol. 2019, 4, 403–404. [Google Scholar] [CrossRef] [PubMed]
  78. Zuo, C.; Shin, S.; Keles, S. atSNP: Transcription factor binding affinity testing for regulatory SNP detection. Bioinformatics 2015, 31, 3353–3355. [Google Scholar] [CrossRef]
  79. McGeary, S.; Lin, K.; Shi, C.; Pham, T.; Bisaria, N.; Kelley, G.; Bartel, D. The biochemical basis of microRNA targeting efficacy. Science 2019, 366, eaav1741. [Google Scholar] [CrossRef]
Jcdd 10 00393 g001
Figure 1. Familial SCAD cohort and prioritized candidate genes. Candidate genes: CDH4, cadherin 4; COL3A1, collagen type III alpha I chain; COL4A1, collagen type IV alpha I chain; COL4A2, collagen type IV alpha 2 chain; COL5A2, collagen type V alpha 2 chain; CTNNB1, catenin beta 1; EPHB4, ephrin type-B receptor 4; IQGAP1, IQ-motif-containing GTPase-activating protein 1; KCNK3, potassium channel subfamily K member 3; LRP2, LDL-receptor-related protein 2; NOTCH1, neurogenic locus notch homolog protein 1.
Figure 1. Familial SCAD cohort and prioritized candidate genes. Candidate genes: CDH4, cadherin 4; COL3A1, collagen type III alpha I chain; COL4A1, collagen type IV alpha I chain; COL4A2, collagen type IV alpha 2 chain; COL5A2, collagen type V alpha 2 chain; CTNNB1, catenin beta 1; EPHB4, ephrin type-B receptor 4; IQGAP1, IQ-motif-containing GTPase-activating protein 1; KCNK3, potassium channel subfamily K member 3; LRP2, LDL-receptor-related protein 2; NOTCH1, neurogenic locus notch homolog protein 1.
Jcdd 10 00393 g001
Jcdd 10 00393 g002
Figure 2. Genetic variant filtering and gene prioritization workflow identified top candidate genes for familial SCAD. A variant filtering scheme was applied to all single-nucleotide variants and insertions/deletions, identifying 11 genes in 12 families. Top candidate gene(s) selection for each family was based on fulfilling the highest number of secondary filtering parameters. Abbreviations: CADD, Combined Annotation Dependent Depletion; CTD, connective tissue disorder; GWAS, genome-wide association study; HGMD, Human Gene Mutation Database; MAF, minor allele frequency; pLI, probability of loss-of-function; SCAD, spontaneous coronary artery dissection.
Figure 2. Genetic variant filtering and gene prioritization workflow identified top candidate genes for familial SCAD. A variant filtering scheme was applied to all single-nucleotide variants and insertions/deletions, identifying 11 genes in 12 families. Top candidate gene(s) selection for each family was based on fulfilling the highest number of secondary filtering parameters. Abbreviations: CADD, Combined Annotation Dependent Depletion; CTD, connective tissue disorder; GWAS, genome-wide association study; HGMD, Human Gene Mutation Database; MAF, minor allele frequency; pLI, probability of loss-of-function; SCAD, spontaneous coronary artery dissection.
Jcdd 10 00393 g002
Jcdd 10 00393 g003
Figure 3. Predicted-damaging variants in familial SCAD altered highly conserved amino acids. Alignment of variants and surrounding residues demonstrated conservation across species. Identical residues are indicated by a dot, and a dash indicates a gap.
Figure 3. Predicted-damaging variants in familial SCAD altered highly conserved amino acids. Alignment of variants and surrounding residues demonstrated conservation across species. Identical residues are indicated by a dot, and a dash indicates a gap.
Jcdd 10 00393 g003
Jcdd 10 00393 g004
Figure 4. Sex-specific arterial tissue expression of prioritized candidate genes. Violin plots were generated by GTEx, demonstrating no differences in sex-specific candidate gene expression in coronary, aorta, and tibial arteries.
Figure 4. Sex-specific arterial tissue expression of prioritized candidate genes. Violin plots were generated by GTEx, demonstrating no differences in sex-specific candidate gene expression in coronary, aorta, and tibial arteries.
Jcdd 10 00393 g004
Jcdd 10 00393 g005
Figure 5. Functional associations between individual F-SCAD candidates and previously reported SCAD genes. Talin 1 (red) and/or GWAS-derived candidate gene proteins EDN1, LRP1, FBN1, and AFAP1 (blue) interacted with 10 of the prioritized candidates for familial SCAD (yellow) directly or through a STRING second shell interactor protein (grey). STRING analysis was set at a high confidence score (0.7).
Figure 5. Functional associations between individual F-SCAD candidates and previously reported SCAD genes. Talin 1 (red) and/or GWAS-derived candidate gene proteins EDN1, LRP1, FBN1, and AFAP1 (blue) interacted with 10 of the prioritized candidates for familial SCAD (yellow) directly or through a STRING second shell interactor protein (grey). STRING analysis was set at a high confidence score (0.7).
Jcdd 10 00393 g005
Jcdd 10 00393 g006
Figure 6. Network interactome of functional associations among F-SCAD candidates. (A) Protein–protein interactions in 10 of the familial prioritized candidates (yellow) and talin 1 (red) directly or a STRING 2nd shell interactor (grey). (B) Shared protein–protein interactions had a high confidence score of <0.9.
Figure 6. Network interactome of functional associations among F-SCAD candidates. (A) Protein–protein interactions in 10 of the familial prioritized candidates (yellow) and talin 1 (red) directly or a STRING 2nd shell interactor (grey). (B) Shared protein–protein interactions had a high confidence score of <0.9.
Jcdd 10 00393 g006
Jcdd 10 00393 g007
Figure 7. Familial candidate genes share pathways and biological processes involved in the cytoskeleton and angiogenesis. Abbreviations: FAK1, focal adhesion kinase 1; GO, Gene Ontology. Asterisks denotes the previously identified TLN1.
Figure 7. Familial candidate genes share pathways and biological processes involved in the cytoskeleton and angiogenesis. Abbreviations: FAK1, focal adhesion kinase 1; GO, Gene Ontology. Asterisks denotes the previously identified TLN1.
Jcdd 10 00393 g007
Table 1. Familial SCAD demographics and clinical characteristics.
Table 1. Familial SCAD demographics and clinical characteristics.
FamilyIDSexRaceAge at Event/
Diagnosis, Yrs		Clinical FeaturesAffected Coronary Artery
SCAD-01II.1FW42SCAD, MHMV: LM, LCx, OM2, D1
II.2F47SCAD, FMD, P, ESOM
SCAD-02II.1FW35SCAD, CT, ESLAD
II.2F45SCADLAD
SCAD-03II.1FW43SCAD, FMD, CTLAD
II.2F45SCAD, MHLAD
SCAD-04II.1FW43SCAD, ESMV: LCx, OM
II.2F59SCADLAD
SCAD-05II.1FW34SCAD, FMD, P, PEOM2
II.2F36SCAD, FMDLAD
SCAD-06II.1FW50SCAD, FMD, ICAAMV: LAD, LCx
II.2M47SCADDA
II.3M40SMAD—
SCAD-07III.1FW41SCAD, FMD, ES, CALAD
III.2F37SCAD, FMD, CTLAD
II.3F65CA, ICAA—
SCAD-08III.1FW50, 59 SCAD, R, ES1. LAD, 2. OM1
III.2F50SCAD, CTOM1
I.1M71CA—
SCAD-09III.2FW42SCAD, ESMV: RCA, PDA
III.3FUnkSCAD, FMD
III.1M25CVA—
SCAD-10II.3FW42, 45, 56SCAD, R, ES1. LCx, 2. MV: LAD, LM, 3. RCA
III.1F39SCAD, MHLAD
II.2FUnkCT—
SCAD-11III.1FW36SCAD, FMDOM1
III.2F44SCAD, CT, MH, PERCA
SCAD-12II.3FW47SCAD, FMDLCx
III.1F38SCAD, POM
II.1FUnkCA—
SCAD-13II.1FW44SCADLAD
I.2F 54SCADLAD
II.2F 45CD—
SCAD-14I.2FW59, 61SCAD, FMD, R1. MV: LCx, OM2, 2. LAD
II.1F32SCAD, PELAD
SCAD-15II.1FW33SCAD, FMD, MHLAD
II.2FUnkSCAD
I.1F65AAA, IAA, PAA—
AAA, abdominal aortic aneurysm; CA, cerebral aneurysm; CD, carotid artery dissection; CT, coronary tortuosity; CVA, cerebrovascular accident; DA, diagonal artery; ES, emotional stress; F, female; FMD, fibromuscular dysplasia; ICAA, internal carotid artery aneurysm; IAA, iliac artery aneurysm; LAD, left anterior descending coronary artery; LCx, left circumflex coronary artery; LM, left main coronary artery; M, male; MH, migraine headache; MV, multivessel; OM, obtuse marginal artery; P, pregnancy-associated; PAA, popliteal artery aneurysm; PE, physical exertion; R, recurrent; RCA, right coronary artery; SMAD, superior mesenteric artery dissection; and W, white.
Table 2. Top candidate genes in familial SCAD. Minor allele frequencies are based on all populations. AA, aortic aneurysm; Aor, aortopathy; AAA, abdominal aortic aneurysm; AD, aortic dissection; CAC, coronary artery calcification; CAD, coronary artery disease; CADD, combined annotation-dependent depletion; ceAD, cervical artery dissection; CoA, coarctation of the aorta; CTD, connective tissue disorder; CTHD, conotruncal heart defects; CSVD, cerebral small vessel disease; CVD, cardiovascular disease; gnomAD, Genome Aggregation Database; GWAS, genome-wide association study; HGMD, Human Genome Mutation Database; IA, intracranial aneurysm; ICH, intracerebral hemorrhage; LVOTO, left ventricular outflow tract obstruction; MGI, Mouse Genome Informatics; PAH, pulmonary arterial hypertension; TAD, thoracic aortic disorder; TGA, transposition of the great arteries; TOF, tetralogy of fallot; VAD, intracranial vertebral artery dissection.
Table 2. Top candidate genes in familial SCAD. Minor allele frequencies are based on all populations. AA, aortic aneurysm; Aor, aortopathy; AAA, abdominal aortic aneurysm; AD, aortic dissection; CAC, coronary artery calcification; CAD, coronary artery disease; CADD, combined annotation-dependent depletion; ceAD, cervical artery dissection; CoA, coarctation of the aorta; CTD, connective tissue disorder; CTHD, conotruncal heart defects; CSVD, cerebral small vessel disease; CVD, cardiovascular disease; gnomAD, Genome Aggregation Database; GWAS, genome-wide association study; HGMD, Human Genome Mutation Database; IA, intracranial aneurysm; ICH, intracerebral hemorrhage; LVOTO, left ventricular outflow tract obstruction; MGI, Mouse Genome Informatics; PAH, pulmonary arterial hypertension; TAD, thoracic aortic disorder; TGA, transposition of the great arteries; TOF, tetralogy of fallot; VAD, intracranial vertebral artery dissection.
Gene Symbol
(Name)		rsIDVariantTranslation
Impact		gnomAD MAF
(%)		CADD
Score
(Percentile)		Deleterious Classification
(HGMD)		Gene
Constraint
(Z-Score/pLI †)		Arterial Tissue Expression—RankPPI
Talin 1 = 1,
GWAS = 2		CTD/
Aor
Gene		Abnormal Vascular Phenotype
(MGI)		Vascular Phenotype
(HGMD)		GWAS
Trait
(GWAS Cat)
CoronaryAortaTibial
CDH4
(cadherin 4)		rs1305825960c.2101T > C
p.S701P		Missense0.000425.1 (75)—2.024746432———CAC, ICH
COL3A1
(collagen type III alpha 1 chain)		rs1238066761c.1336C > T
p.R446C		Missense0.000432 (95)—4.096571, 2YesAD, aorta smooth muscle morphologyTAD, AA, CSVD, AAA, SCAD—
COL4A1
(collagen type IV alpha 1 chain)		rs145172612c.1588C > T
p.P530S		Missense0.0329.4 (90)CM18181943.021321, 2—ICH, retinal vascular morphologyceAD, CSVD, ICH, CAD, TGAArterial stiffness, CVD, CAD
COL4A2
(collagen type IV alpha 2 chain)		rs746743018c.451G > A
p.G151S		Missense0.00224.8 (75)—2.192131, 2—Aorta stenosis,
cranial blood
vascular morphology		CSVD, ICHArterial stiffness, carotid artery thickness, CAC, CAD
rs1464563247c.1237G > A
p.G413R		Missense0.00126.6 (85)—
COL5A2
(collagen type V alpha 2 chain)		—c.2705A > C
p.Q902P		MissenseNR27.5 (90)—2.446481, 2YesVascular congestion, cardiovascular system physiologyceAD, CTHD, AD, IA, SCAD—
CTNNB1
(catenin beta 1)		rs35288908c.860A > G
p.N287S		Missense0.0921.2 (30)CM0437573.85131531, 2—Pulmonary artery, vascular endothelial cell and vitelline morphology——
EPHB4
(ephrin type-B receptor 4)		rs146674844c.1384G > A
p.G462R		Missense0.0324.8 (75)—2.30303738——Angiogenesis, vascular branching, pulmonary artery and vitelline morphologyPAH, LVOTO—
IQGAP1
(IQ-motif-containing GTPase-activating protein1)		—c.3196A > G
p.N1066D		MissenseNR27.7 (90)—2.445431, 2——TOF and other cardiac abnormalities—
KCNK3
(potassium channel subfamily K
member 3)		—c.781_794delCGCGCGCTGCTCAC
p.R261fsTer		Frameshift *NR——0.90 †5137———PAHCVD, ICH
LRP2
(LDL-receptor-related protein 2)		—c.11093T > A
p.V3698E		MissenseNR24.6 (75)—2.073531342—Aortic arch morphologyPAH—
NOTCH1
(neurogenic locus notch homolog protein 1)		rs182330532c.3835C > T
p.R1279C		Missense *0.0634 (99)—3.452628221, 2YesAngiogenesis, vasculogenesis,
vitelline morphology		TAD, PAH, VAD, CoA, TOF, ceADCAC
* Predicted likely pathogenic by American College of Medical Genetics and Genomics criteria. † pLI, probability of loss of function intolerance score.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Turley, T.N.; Theis, J.L.; Evans, J.M.; Fogarty, Z.C.; Gulati, R.; Hayes, S.N.; Tweet, M.S.; Olson, T.M. Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways. J. Cardiovasc. Dev. Dis. 2023, 10, 393. https://doi.org/10.3390/jcdd10090393

AMA Style

Turley TN, Theis JL, Evans JM, Fogarty ZC, Gulati R, Hayes SN, Tweet MS, Olson TM. Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways. Journal of Cardiovascular Development and Disease. 2023; 10(9):393. https://doi.org/10.3390/jcdd10090393

Chicago/Turabian Style

Turley, Tamiel N., Jeanne L. Theis, Jared M. Evans, Zachary C. Fogarty, Rajiv Gulati, Sharonne N. Hayes, Marysia S. Tweet, and Timothy M. Olson. 2023. "Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways" Journal of Cardiovascular Development and Disease 10, no. 9: 393. https://doi.org/10.3390/jcdd10090393

APA Style

Turley, T. N., Theis, J. L., Evans, J. M., Fogarty, Z. C., Gulati, R., Hayes, S. N., Tweet, M. S., & Olson, T. M. (2023). Identification of Rare Genetic Variants in Familial Spontaneous Coronary Artery Dissection and Evidence for Shared Biological Pathways. Journal of Cardiovascular Development and Disease, 10(9), 393. https://doi.org/10.3390/jcdd10090393

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop