PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
J Invest Dermatol. Author manuscript; available in PMC 2014 April 1.
Published in final edited form as:
PMCID: PMC3971866
NIHMSID: NIHMS404715

Copy Number Variation analysis in 98 individuals with PHACE syndrome

Abstract

PHACE syndrome is the association of large segmental facial hemangiomas and congenital anomalies, such as posterior fossa malformations, cerebral arterial anomalies, coarctation of the aorta, eye anomalies and sternal defects. To date, the reported cases of PHACE syndrome have been sporadic suggesting that PHACE may have a complex pathogenesis. We report here genomic copy number variation (CNV) analysis in 98 individuals with PHACE syndrome as a first step in deciphering a potential genetic basis of PHACE syndrome. A total of 3,772 CNVs (2,507 duplications and 1,265 deletions) were detected in 98 individuals with PHACE syndrome. CNVs were then eliminated if they failed to meet established criteria for quality, spanned centromeres or did not contain genes. CNVs were defined as “rare” if not documented in the Database of Genomic Variants. Ten rare CNVs were discovered (size range: 134–406 kb), located at 1q32.1, 1q43, 3q26.32-3q26.33, 3p11.1, 7q33, 10q24.32, 12q24.13, 17q11.2, 18p11.31, and Xq28. There were no rare CNV events that occurred in greater than one subject. Therefore, further study is needed to determine the significance of these CNVs in the pathogenesis of PHACE syndrome.

Keywords: PHACE syndrome, PHACES syndrome, Pascual-Castroviejo II syndrome, aortic arch anomaly, hemangioma, copy number variation

Introduction

Children with large infantile hemangiomas (IH) on the face and scalp are at risk for associated birth defects. In 1978, Pascual-Castroviejo described a female patient with a facial hemangioma and Dandy-Walker syndrome with abnormal intracranial cerebral arteries(Pascual-Castroviejo, 1978b). He subsequently described seven additional female patients with hemangioma and neurovascular anomalies (Pascual-Castroviejo, 1978a). In 1985, the term sternal malformation/vascular dysplasia syndrome was used to describe a child with sternal cleft and cutaneous, craniofacial hemangiomas as well as airway and liver hemangiomas (Hersh et al., 1985). In 1996, Frieden et al proposed the acronym PHACE (posterior fossa malformations, infantile hemangiomas, arterial anomalies of the great cerebral vessels, coarctation of the aorta, eye anomalies and sternal defects) to describe this condition (OMIM: 606519) (Frieden et al., 1996). Children with PHACE are at risk for cerebral vasculopathy that can lead to arterial ischemic stroke in a subset of individuals (Drolet et al., 2006; Siegel et al., 2012). Hemangiomas on the lower body can be associated with congenital anomalies in LUMBAR syndrome (lower body hemangioma, urogenital anomalies, ulceration, myelopathy, bone deformities, anorectal or arterial anomalies and renal anomalies)(Iacobas et al., 2010). Recently, identification of the genes for the overgrowth syndromes associated with vascular malformations have been rapidly unfolding (Kurek et al., 2012; Riviere et al., 2012), however the genes for the hemangioma syndromes have yet to be elucidated.

All reported cases of PHACE to date have been sporadic, yet the consistent phenotype and the characteristic brain and aortic anomalies provide evidence for a genetic etiology. PHACE has a significant female predominance suggesting that there may be a contribution from the X chromosome or hormonal influences (Haggstrom et al., 2007). X-linked inheritance has been suggested, although males do not have increased severity (Levin and Kaler, 2007; Metry et al., 2008). Preliminary X-inactivation studies in this PHACE cohort did not identify significant skewing in most affected females or their mothers. Although linkage studies have been performed on solitary IH identifying a locus at 5q, and subsequent studies have pointed to the involvement of VEGFR2 (Jinnin et al., 2008; Walter et al., 1999), there have been no linkage studies performed on PHACE specifically. There have been no significant prenatal environmental factors or toxins identified in association with PHACE (Metry et al., 2006).

There is growing evidence that small structural chromosomal abnormalities contribute to a broad array of complex congenital malformations (Southard et al., 2012). High-density arrays now provide a powerful method for evaluating genome-wide copy number variation (CNV) and have been used to advance knowledge of the genetic basis of numerous disorders including thoracic aortic aneurysms (Prakash et al., 2010), autism spectrum disorders (Glessner et al., 2009) and developmental brain disorders (Aldinger et al., 2009; Pinto et al., 2010). The objective of this study is to comprehensively evaluate the contribution of copy number variation in PHACE syndrome.

Results

Ninety-eight individuals diagnosed with PHACE were genotyped and analyzed for CNVs. Of these, 94 individuals with PHACE were genotyped using the Affymetrix Human Mapping GeneChip 6.0 arrays and analyzed with the Genotyping Console™ Software (Figure 1: Flowchart of CNV Analysis and Identification). A collaborating group also provided an additional cohort of 4 individuals with PHACE which were run and analyzed with Roche Nimblegen programs and these results were included in this analysis (see Supplemental methods). Eighty-two individuals with PHACE were female and 16 male. Additional demographic data are summarized in Table 1. A control group included 91 individuals without PHACE from an unrelated research cohort (see methods for a description of this cohort).

Figure 1
Flowchart of CNV Analysis and Identification
Table 1
Clinical Characteristics of 98 individuals with PHACE Syndrome

CNV in 98 PHACE syndrome cases

First, we looked at CNV burden and found a total of 3,696 CNVs (2,473 duplications and 1,223 deletions) in 94 individuals with PHACE using Genotyping Console, averaging 39 CNVs per individual (18–83 CNVs) with a mean length of 328 kb (range: 25 kb-23.8 MB). In the control cohort of 91 individuals, a total of 3,484 CNVs (2,506 duplications and 978 deletions) were identified, averaging 38 CNVs per individual (18–91 CNVs) with a mean length of 267 kb (25 kb-23.4 MB). The average number of CNVs per individual was not found to be statistically different between subjects and controls (p=0.91). In data provided by a collaborating group from an additional cohort of 4 individuals with PHACE a total of 76 CNVs (34 duplications and 42 deletions) with a mean length of 373 kb (range: 27 kb-3.9 Mb) were detected using the Roche NimbleGen array. Thus, a total of 3,772 CNVs (2,507 duplications and 1,265 deletions) were detected in 98 individuals with PHACE syndrome. CNVs were then eliminated from further analysis if they failed to meet established criteria for quality, spanned centromeres or did not contain genes. Of the CNVs meeting our criteria as described in Figure 1, only two CNVs were excluded based on identification in the control cohort (1q34 and 2p13.2). In these 98 individuals with PHACE, we identified rare CNVs greater than 100kb at 1q32.1, 1q43, 3q26-3q26.33, 3p11.1, 7q33, 10q24.32, 12q24.13, 17q11.2, 18p11.31, and Xq28. The genes contained in the reported CNVs are shown in Supplemental Table 1 and clinical phenotypes are detailed in Supplemental Table 2. There were no rare CNV events that occurred in greater than one individual with PHACE.

Validation

Taqman qPCR was used to validate 9 discovered rare CNV regions. The 7q33 region was validated with the Nimblegen chip. Taqman probes for selected genes in each genomic region were designed as follows: 1) KLHL12, CYB5R1, and PPFIA4 (1q32.1), 2) FH and KMO (1q43), 3) PIK3CA and MFN1 (3q26.32-3q26.33), 4) EPHA3 (3q11.1), 5) ARL3 (10q24.32), 6) TPCN1 and SLC24A6 (12q24.13), 7) TNFAIP1 and VTN (17q11.2), 8) EMILIN2 and MYOM1 (18p11.31), and 9) AFF2 (Xq28). Additional candidate regions (tested but not listed above) that failed to validate with Taqman qPCR were excluded from further analysis.

De novo versus Inherited CNVs

The de novo versus inherited status was determined by testing the mother and father for the identified CNV using the Taqman qPCR analysis. Parental DNA for either one or both parents was not available for 5 cases. Inherited CNVs included 1q32.1, 3p11.1, 10q24.32, 17q11.2, and 18p11.31 (Supplemental Table 1). De novo status for all other validated CNVs was indeterminate as one or both parental DNA samples were unavailable for analysis

Functional network pathway analysis

We screened the CNV regions for known genes important in vasculogenesis such as TGFBR1, TGFBR2, VEGF, NOTCH1, and NOTCH2; however, CNVs containing these specific genes were not identified in this cohort. As a next step in the analysis, validated CNVs were analyzed with the Ingenuity program to identify gene networks and enrich for common pathways (see methods). The analysis identified direct and indirect interactions between the genes in the CNV and the AKT/PIK3/VEGF pathways (Figure 2). A duplication containing the gene PIKC3A is of interest because PIKC3A has indirect interactions with mTOR, NOS2, EGF and HIF-1α and a direct interaction with ERK and AKT1. MYOG, CYP17A1, and CHI3L1 also showed interactions with these pathways.

Figure 2
Network of interacting PHACE syndrome CNV genes

Discussion

As a first step in determining the genetic influence in the development of PHACE, we utilized genome wide CNV analysis in a discovery cohort of 98 individuals with PHACE to evaluate CNV regions. This is the largest series to date analyzing the contribution of copy number variation in PHACE. We report here several rare CNVs (see methods for selection criteria) which were not found in the Database of Genomic Variants. We have further reviewed the CNVs in the Database of Genomic Structural Variation (dbVar) containing a dataset with copy number variants from 15,767 children with developmental delay. There was no overlap between the CNVs in individuals with PHACE and the CNVs reported as pathogenic in the publication from Cooper et al(Cooper et al., 2011). Overlap was noted with the CNVs recorded in dbVar from the Cooper cohort, however the majority of the reported overlapping CNVs were much larger, encompassing several additional genes. In the dbVar database, there was no reported clinical interpretation, validation or testing of inheritance for these CNVs. As these variants have not been previously linked with a specific disease phenotype other than developmental delay, and in rare cases craniofacial dysmorphism or unspecified abnormalities of the cardiovascular system, they should be considered 'variants of unknown significance’ until larger datasets for PHACE are investigated or functional studies are performed. In this study we did not find any large, rare CNVs occurring in more than one individual with PHACE.

The CNVs reported in this paper contain genes that interact with the HIF-1α, VEGF, AKT and mTOR pathways (Figure 2) which have been implicated in angiogenesis and hemangioma pathogenesis (Boscolo et al., 2011; Chen et al., 2005; Karar and Maity, 2011). The CNV gain at 3q26.32 contains the PIK3CA gene. PIK3CA (phosphoinositide 3-kinase catalytic subunit α) interacts with AKT1, VEGF and HIF-1α (Figure 2). PIK3CA, specifically, has been demonstrated to promote angiogenesis via VEGF with a positive correlation between PIK3CA and VEGF expression in ovarian cancer cells (Zhang et al., 2003). The mutations in the PIK3CA gene have recently been implicated in a condition characterized by megalencephaly and growth dysregulation with variable asymmetry, skin capillary malformations, distal limb malformations, and variable cortical malformation (MCAP) (Mirzaa et al., 2012; Riviere et al., 2012).

In addition, PIK3CA gene mutations have been identified in CLOVES syndrome (congenital lipomatous overgrowth with vascular, epidermal, and skeletal anomalies)(Kurek et al., 2012). The 1q32.1 region also contains the liprin-α4 gene promoter (PPFIA4) which is directly activated by binding of the hypoxia-inducible factor 1α (HIF-1α) (Mattauch et al., 2010). HIF-1α functions as a transcriptional regulator of the adaptive response to hypoxia and plays a role in embryonic vascularization. A portion of the MYOM1 gene was present in both a 26 kb loss and a 247 kb gain, in two unrelated individuals with PHACE. The significance of this finding is uncertain. The myomesin genes play an important structural and functional role in striated muscle and myomesin mutations have been identified in a family with hypertrophic cardiomyopathy (Siegert et al., 2011).

Genes involved with embryonic development

The candidate gene VTN in region 17q11.2 interacts with a biological network involved in embryonic cranial mesenchyme and vascular tissue development, as well as endothelial-mesenchymal interactions in development (Arciniegas et al., 2006). Such interactions with the cranial mesenchyme have been shown to be important in normal cerebellar development and disruption can lead to Dandy-Walker malformation (Aldinger et al., 2009). This CNV in this female proband was inherited from the mother; therefore we consider this to be a variant of unknown significance. A gain of the ephrin receptor tyrosine kinase (EPHA3) gene at chromosome 3p11.1 was identified in an individual with a sternal cleft and supraumbilical raphe. EPHA3 is involved in the guidance and coordination of the migration of motor and sensory axons during development and in cardiac valve development (Bevins et al., 2011; Jayasena et al., 2005; Oates et al., 1999; Stephen et al., 2007). A gain of the ADORA1 gene was found in a patient with arterial anomalies and coarctation of the aorta. The adenosine A1 receptor (ADORA1) gene has been shown to participate in axon guidance for the retino-tectal patterning through the engrailed homeoprotein transcription factors and ephrin A5 (Stettler et al., 2012). It is notable that a gain in the EMILIN2 gene at 18p11.31 was identified in a patient with a posterior fossa malformation and both carotid and cerebral arterial anomalies. EMILIN2 protein expression has been demonstrated in the mouse heart and aorta (Doliana et al., 2001; Sa and Hoover-Plow, 2011). Fibrosis and decreased elastic fibers seen in aorta tissue after coarctation repair has been reported, suggesting that abnormalities in elastic fibers could potentially play a role in the pathogenesis of PHACE (Prada et al., 2010).

Summary

PHACE is a complex condition which is likely to have multifactorial risk factors. Although several rare CNVs were identified in this cohort, none were found in more than one patient. Five discovered CNVs, 1q32.1, 3p11.1, 10q24.32, 17q11.2, and 18p11.31 were inherited, although there was no family history of either PHACE or IH in any of these individuals. It is interesting to note the chromosomal regions, 1q32.1, 3p11.1, 10q24.32, 17q11.2, and 18p11.31 were enriched for genes involved in the HIF-1α, VEGF and mTOR signaling pathways. Larger datasets will be necessary to determine the significance of the CNVs reported in this manuscript. These findings are consistent with our current understanding of PHACE as a complex condition and more sophisticated approaches, including next-generation sequencing of both germline and tissue DNA, will likely be necessary to determine if there are genetic factors in the pathogenesis of this disorder.

Limitations

The greatest challenges with CNV analysis are discriminating the effects of common and rare CNV polymorphisms and the disease relevance of structural chromosomal abnormalities. We used the existing databases of published CNVs including the University of California Santa Cruz (UCSC) Genome Browser, the Database of Genomic Variants (DGV) and others to filter variants that appear to be common in the general population. We recognize a limitation of this study is that we may have excluded some of the CNVs with a potential influence in PHACE based on a previously reported association with a different disorder. When available, we used a family-based approach to determine the de novo status. Array CNV analysis is limited by its resolution to detect events. It is also possible that the coronary artery disease control group may have shared some common CNV with the PHACE cohort, wherein the causative CNV could have been masked. We recognize this as another limitation of this study.

Materials and Methods

For additional methods for the samples run on Nimblegen array, see the supplemental methods.

Study Populations

The study has ongoing approval from the Institutional Review Board (IRB) at Children’s Hospital of Wisconsin (CHW) and was approved previously by the IRBs at Oregon Health and Science University (OHSU), the University of Utah, and the Committee for Human Research (CHR) at the University of California at San Francisco (UCSF). Written informed consent was obtained from all study subjects upon enrollment into the PHACE International Clinical Registry and Genetic Repository and includes information about DNA sharing and Database of Genotypes and Phenotypes and adherence to the Helinski Guidelines. Ninety-four individuals with PHACE were genotyped on the Affymetrix Human Mapping GeneChip 6.0 arrays and analyzed with the Genotyping Console™ Software-Affymetrix program. The diagnosis of definite or possible syndrome was made at participating institutions using previously described criteria (Metry et al., 2009) and confirmed by central review of magnetic resonance imaging (MRI), magnetic resonance angiography (MRA), or echocardiogram reports. Internal controls for the study included 91 subjects with coronary artery disease who did not carry a diagnosis of PHACE. Individuals in the control cohort were consented under an unrelated research study which allows DNA information sharing for research related purposes. Parent and control cohorts were genotyped on the same platform and subjected to equivalent QC measures and CNV analysis used for affected cases.

DNA and Tissue Repository

The DNA and tissue repository for the PHACE syndrome International Clinical Registry and Genetic Repository has been maintained in the Children’s Research Institute at MCW. DNA was extracted from saliva obtained from buccal swabs and blood collected from each subject (4 ml in an EDTA tube). The quantity and quality of the DNA was measured with the Nanodrop (Thermo Scientific).

Genotyping and CNV Data Analysis and Novelty Assessment

Two hundred eighteen samples (118 PHACE subjects, and 100 control samples) were hybridized to the Affymetrix 6.0 GeneChips. Genotyping was performed with the Affymetrix Human Mapping GeneChip 6.0 arrays to assay a minimum of roughly 1.9M probes for analysis per standard protocol (see Supplemental methods). CEL files were exported to Genotyping Console (Affymetrix) for analysis per standard protocols. Stringent quality parameters were applied. Samples with greater than 92 CNV segments, out-of-bound, and duplicate patient samples were excluded from further analysis, resulting in a total of 94 individuals with PHACE and 91 controls for the study. Regions of CNV in our participants with PHACE were compared to CNVs previously reported in public CN databases (Toronto Database of Genomic Variants (DGV)) to determine regions that were novel to individuals with PHACE. An additional search for CNV regions in Decipher was conducted to look for any associated phenotypes and in OMIM for known diseases. Candidate genes within discovered genomic regions were determined for its potential role in disease pathogenesis based on information gathered from sources in PubMed, UCSC Genome Browser, and Ingenuity. We used the following criteria for CNV detection: 1) Size greater than 100 kb, 2) less than 5% overlap with reported CNVs in the DGV or our control cohort, and 3) we only included regions containing genes. CNVs with a less than 100 kb copy loss or gain were noted if found within a larger CNV region discovered in another individual. An average distance between markers of ≤ 2 kb was required to determine adequate probe density. Once rare CNVs were discovered in our patient cohort, validation was performed with the Applied Biosystems TaqMan quantitative polymerase chain reaction (qPCR) and de novo status was ascertained when parental genotypes were available. Finally, the validated CNV were cross-referenced against the published neurodevelopmental CNVs at http://www.ncbi.nlm.nih.gov/dbvar/studies/nstd54/ (Cooper et al., 2011).

Quantitative Polymerase Chain Reaction Validation

Findings from our CNV analysis with genotyping console were validated with the Applied Biosystems TaqMan qPCR. Primers and probes used to detect CNVs with qPCR were ordered from Applied Biosystems (pre-mixed assay) (see Supplemental methods).

Functional Analysis: Ingenuity Network Pathways

Candidate genes within discovered genomic regions from CNV analysis with Genotyping Console were imported into the Ingenuity Pathways Analysis (IPA) Web server for functional analysis (maximum 35 genes). The IPA software utilizes the IngenuityR Knowledge Base (http://ingenuity.com, accessed Jan 2012) which contains information on millions of biological interactions and functional annotations derived from published data. This referenced information is used to model relationships between genes, proteins, complexes, cells, cellular components, and diseases amongst other molecules. The Path Explorer function of the IPA program searches for interactions between genes and was utilized to create the shortest interconnected pathway for genes contained in detected CNV regions. Shared pathways were subsequently noted.

Supplementary Material

01

Acknowledgements

We would like to thank Shawna Joachim for assisting with study coordination and manuscript preparation and Gina Im for assisting with data management. We are grateful to the families for their participation in the study and the physicians who referred patients for this study.

Funding:

This publication was made possible with support from the Oregon Clinical and Translational Research Institute (OCTRI), grant# UL1 RR024140 from the National Center for Research Resources (NCRR) and the National Center for Advancing Translational Sciences (NCATS) of the National Institutes of Health (NIH). Additional funding sources include Society for Pediatric Dermatology (DHS), Dermatology Foundation (DHS), Greater Milwaukee Foundation (DHS, BD), Families, Society of Pediatric Otolaryngology, Public Health Services research grants UL1-RR025764 and C06-RR11234 from the NCRR, the Children’s Health Research Center and Clinical Genetics Research Program at the University of Utah, American Society of Pediatric Otolaryngology (FG, PT, SM, DHS). JS was funded by NIH/NHLBI K08HL092970. The project was supported in part by funds from Children’s Research Institute, Medical College of Wisconsin (UB).

List of abbreviations

CNV
copy number variation
PHACE
an acronym for Posterior fossa brain malformation, large facial Hemangiomas, Arterial anomalies, Cardiac anomalies and aortic coarctation, and Eye abnormalities
DGV
Database of Genomic Variants
dbVar
Database of Genomic Structural Variation

Footnotes

References

  • Aldinger KA, Lehmann OJ, Hudgins L, et al. FOXC1 is required for normal cerebellar development and is a major contributor to chromosome 6p25.3 Dandy-Walker malformation. Nat Genet. 2009;41:1037–1042. [PMC free article] [PubMed]
  • Arciniegas E, Neves YC, Carrillo LM. Potential role for insulin-like growth factor II and vitronectin in the endothelial-mesenchymal transition process. Differentiation. 2006;74:277–292. [PubMed]
  • Bevins N, Lemke G, Reber M. Genetic dissection of EphA receptor signaling dynamics during retinotopic mapping. J Neurosci. 2011;31:10302–10310. [PMC free article] [PubMed]
  • Boscolo E, Mulliken JB, Bischoff J. VEGFR-1 mediates endothelial differentiation and formation of blood vessels in a murine model of infantile hemangioma. Am J Pathol. 2011;179:2266–2277. [PubMed]
  • Chen D, Lin XX, Li W. The relationship between the expression of HIF-1alpha and the angiogenesis in infancy hemangioma. Zhonghua Zheng Xing Wai Ke Za Zhi. 2005;21:115–118. [PubMed]
  • Cooper GM, Coe BP, Girirajan S, et al. A copy number variation morbidity map of developmental delay. Nat Genet. 2011;43:838–846. [PMC free article] [PubMed]
  • Doliana R, Bot S, Mungiguerra G, et al. Isolation and characterization of EMILIN-2, a new component of the growing EMILINs family and a member of the EMI domain-containing superfamily. J Biol Chem. 2001;276:12003–12011. [PubMed]
  • Drolet BA, Dohil M, Golomb MR, et al. Early stroke and cerebral vasculopathy in children with facial hemangiomas and PHACE association. Pediatrics. 2006;117:959–964. [PubMed]
  • Frieden IJ, Reese V, Cohen D. PHACE syndrome. The association of posterior fossa brain malformations, hemangiomas, arterial anomalies, coarctation of the aorta and cardiac defects, and eye abnormalities. Arch Dermatol. 1996;132:307–311. [PubMed]
  • Glessner JT, Wang K, Cai G, et al. Autism genome-wide copy number variation reveals ubiquitin and neuronal genes. Nature. 2009;459:569–573. [PMC free article] [PubMed]
  • Haggstrom AN, Drolet BA, Baselga E, et al. Prospective study of infantile hemangiomas: demographic, prenatal, and perinatal characteristics. J Pediatr. 2007;150:291–294. [PubMed]
  • Hersh JH, Waterfill D, Rutledge J, et al. Sternal malformation/vascular dysplasia association. Am J Med Genet. 1985;21:177–186. 201–202. [PubMed]
  • Iacobas I, Burrows PE, Frieden IJ, et al. LUMBAR: association between cutaneous infantile hemangiomas of the lower body and regional congenital anomalies. J Pediatr. 2010;157:795–801. e1–e7. [PubMed]
  • Jayasena CS, Flood WD, Koblar SA. High EphA3 expressing ophthalmic trigeminal sensory axons are sensitive to ephrin-A5-Fc: implications for lobe specific axon guidance. Neuroscience. 2005;135:97–109. [PubMed]
  • Jinnin M, Medici D, Park L, et al. Suppressed NFAT-dependent VEGFR1 expression and constitutive VEGFR2 signaling in infantile hemangioma. Nat Med. 2008;14:1236–1246. [PMC free article] [PubMed]
  • Karar J, Maity A. PI3K/AKT/mTOR Pathway in Angiogenesis. Front Mol Neurosci. 2011;4:51. [PMC free article] [PubMed]
  • Kurek KC, Luks VL, Ayturk UM, et al. Somatic Mosaic Activating Mutations in PIK3CA Cause CLOVES Syndrome. Am J Hum Genet. 2012;90:1108–1115. [PubMed]
  • Levin JH, Kaler SG. Non-random maternal X-chromosome inactivation associated with PHACES. Clin Genet. 2007;72:345–350. [PubMed]
  • Mattauch S, Sachs M, Behrens J. Liprin-alpha4 is a new hypoxia-inducible target gene required for maintenance of cell-cell contacts. Exp Cell Res. 2010;316:2883–2892. [PubMed]
  • Metry D, Heyer G, Hess C, et al. Consensus Statement on Diagnostic Criteria for PHACE Syndrome. Pediatrics. 2009;124:1447–1456. [PubMed]
  • Metry DW, Haggstrom AN, Drolet BA, et al. A prospective study of PHACE syndrome in infantile hemangiomas: demographic features, clinical findings, and complications. Am J Med Genet A. 2006;140:975–986. [PubMed]
  • Metry DW, Siegel DH, Cordisco MR, et al. A comparison of disease severity among affected male versus female patients with PHACE syndrome. J Am Acad Dermatol. 2008;58:81–87. [PubMed]
  • Mirzaa GM, Conway RL, Gripp KW, et al. Megalencephaly-capillary malformation (MCAP) and megalencephaly-polydactyly-polymicrogyria-hydrocephalus (MPPH) syndromes: two closely related disorders of brain overgrowth and abnormal brain and body morphogenesis. Am J Med Genet A. 2012;158A:269–291. [PubMed]
  • Oates AC, Lackmann M, Power MA, et al. An early developmental role for eph-ephrin interaction during vertebrate gastrulation. Mech Dev. 1999;83:77–94. [PubMed]
  • Pascual-Castroviejo I. Vascular and nonvascular intracranial malformation associated with external capillary hemangiomas. Neuroradiology. 1978a;16:82–84. [PubMed]
  • Pascual-Castroviejo I. Vascular changes in cerebellar developmental defects. Neuroradiology. 1978b;16:58–60. [PubMed]
  • Pinto D, Pagnamenta AT, Klei L, et al. Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010;466:368–372. [PMC free article] [PubMed]
  • Prada F, Mortera C, Bartrons J, et al. Complex aortic coarctation and PHACE syndrome. Rev Esp Cardiol. 2010;63:1367–1370. [PubMed]
  • Prakash SK, LeMaire SA, Guo DC, et al. Rare copy number variants disrupt genes regulating vascular smooth muscle cell adhesion and contractility in sporadic thoracic aortic aneurysms and dissections. Am J Hum Genet. 2010;87:743–756. [PubMed]
  • Riviere JB, Mirzaa GM, O'Roak BJ, et al. De novo germline and postzygotic mutations in AKT3, PIK3R2 and PIK3CA cause a spectrum of related megalencephaly syndromes. Nat Genet. 2012;44:934–940. [PMC free article] [PubMed]
  • Sa Q, Hoover-Plow JL. EMILIN2 (Elastin microfibril interface located protein), potential modifier of thrombosis. Thromb J. 2011;9:9. [PMC free article] [PubMed]
  • Siegel DH, Tefft KA, Kelly T, et al. Stroke in children with posterior fossa brain malformations, hemangiomas, arterial anomalies, coarctation of the aorta and cardiac defects, and eye abnormalities (PHACE) syndrome: a systematic review of the literature. Stroke. 2012;43:1672–1674. [PubMed]
  • Siegert R, Perrot A, Keller S, et al. A myomesin mutation associated with hypertrophic cardiomyopathy deteriorates dimerisation properties. Biochem Biophys Res Commun. 2011;405:473–479. [PubMed]
  • Southard AE, Edelmann LJ, Gelb BD. Role of copy number variants in structural birth defects. Pediatrics. 2012;129:755–763. [PubMed]
  • Stephen LJ, Fawkes AL, Verhoeve A, et al. A critical role for the EphA3 receptor tyrosine kinase in heart development. Dev Biol. 2007;302:66–79. [PubMed]
  • Stettler O, Joshi RL, Wizenmann A, et al. Engrailed homeoprotein recruits the adenosine A1 receptor to potentiate ephrin A5 function in retinal growth cones. Development. 2012;139:215–224. [PubMed]
  • Walter JW, Blei F, Anderson JL, et al. Genetic mapping of a novel familial form of infantile hemangioma. Am J Med Genet. 1999;82:77–83. [PubMed]
  • Zhang L, Yang N, Katsaros D, et al. The oncogene phosphatidylinositol 3'-kinase catalytic subunit alpha promotes angiogenesis via vascular endothelial growth factor in ovarian carcinoma. Cancer Res. 2003;63:4225–4231. [PubMed]