Human hepatic vascular cavernomas, the most common benign tumor of the liver, were described in the mid-1800s, yet the mechanisms for their formation and effective treatments remain unknown. Here, we demonstrate gain-of-function mutations in KRAS or BRAF genes within liver endothelial cells as a causal mechanism for hepatic vascular cavernomas. We identified gain-of-function mutations in KRAS or BRAF genes in pathological liver tissue samples from patients with hepatic vascular cavernomas. Mice expressing these human KRASG12D or BRAFV600E mutations in hepatic endothelial cells recapitulated the human hepatic vascular cavernoma phenotype of dilated sinusoidal capillaries with defective branching patterns. KRASG12D or BRAFV600E induced “zipper-like” contiguous expression of junctional proteins at sinusoidal endothelial cell–cell contacts, switching capillaries from branching to cavernous expansion. Pharmacological or genetic inhibition of the endothelial RAS–MAPK1 signaling pathway rescued hepatic vascular cavernoma formation in endothelial KRASG12D- or BRAFV600E-expressing mice. These results uncover a major cause of hepatic vascular cavernomas and provide a road map for their personalized treatment.
Sporadic cavernous hemangiomas or vascular cavernomas of the liver, classified as slow-flow venous malformations, are the most common benign tumor of the liver (Tait et al., 1992; Graivier et al., 1967). These cavernomas are composed of dilated endothelial-lined vascular spaces that are irregular in arrangement and sizes (Hamilton and Holmes, 1950). Large hepatic vascular cavernomas, >5 cm in diameter, often cause lethal complications such as hepatic rupture, consumption coagulopathy, and cardiac failure (Hendrick, 1948; Graivier et al., 1967; Ribeiro et al., 2010). Although Virchow (1863) and Frerichs (1860) described hepatic vascular cavernomas as a distinct clinical entity in the mid-1800s, their genetic etiology, molecular mechanism, and effective treatment remain undefined.
Human patients with hepatic vascular cavernomas exhibit somatic gain-of-function mutations in KRAS or BRAF genes in pathological liver tissue samples
Using next-generation deep sequencing and Sanger sequencing, along with bridged nucleic acid (BNA) clamp–mediated PCR assay, we sequenced hepatic vascular cavernoma tissue samples from 39 patients diagnosed with hepatic cavernous hemangioma or hemangioma (Table 1 and Fig. S1 A). Our analyses revealed KRAS and/or BRAF gain-of-function pathogenic mutations in ∼36% of cases (14 of 39; Table 1 and Fig. 1, A and B). Activating KRAS mutations (n = 10) were found at hotspot codon 12 and/or 13: p.Gly-Asp and/or p.Gly-Ser with an allele frequency ranging from ∼0.1 to 0.75% (Table 1 and Fig. 1 A). Among these, four cases exhibited mosaic KRAS mutations. BRAF gain-of-function mutations (n = 3) were identified at hotspot codon 600: p.Val600Met with an allele frequency ∼0.2% (Table 1 and Fig. 1 B). Interestingly, one case revealed activating mutations in both KRAS (p.Gly12Ser) and BRAF (p.Val600Met) with a variant allele load of ∼0.2% (Table 1). These hepatic lesions exhibited cavernous spaces filled with erythrocytes and lined by lymphatic vessel endothelial hyaluronan receptor (LYVE1)+/Cadherin 5 (CDH5)+/ platelet and endothelial cell adhesion molecule 1 (PECAM1)+ endothelial cells, malformed sinusoidal capillary vessels, and normal endothelial cell proliferation (Fig. 1 C and Fig. S1, B and C).
Endothelial KrasG12D or BrafV600E gain-of-function mutations cause hepatic vascular cavernomas in mice
Dysregulation of endothelial cells underlies the formation of vascular anomalies. To investigate the endothelial cell–specific role of KRAS or BRAF gain-of-function mutations in hepatic vascular cavernomas, we generated KrasG12D F/+; Cdh5CreERT2 and BrafV600E F/+; Cdh5CreERT2 mice, in which tamoxifen administration activates the KrasG12D or BrafV600E allele expression specifically in Cdh5+ endothelial cells. Tamoxifen-treated KrasG12D F/+; Cdh5CreERT2 and BrafV600E F/+; Cdh5CreERT2 mice rapidly developed sinusoidal capillary dilation and hepatic vascular cavernomas with complete penetrance at perinatal and adult stages (Fig. 2, A and B; and Fig. S2, A–D). Histological and immunological analyses of tamoxifen-treated KrasG12D F/+; Cdh5CreERT2 and BrafV600E F/+; Cdh5CreERT2 mice revealed cavernous, malformed, erythrocyte-filled liver sinusoidal spaces lined by Lyve1+/Cdh5+/Pecam1+ endothelial cells, and normal endothelial cell proliferation consistent with human phenotype (Fig. 2, A and B; Fig. S2, A–D; and Fig. S3, A–D). In addition, Podoplanin (Pdpn)+ lymphatic vessels appear normal and do not contribute to hepatic vascular cavernomas in tamoxifen-treated KrasG12D F/+; Cdh5CreERT2 and BrafV600E F/+; Cdh5CreERT2 mice (Fig. 2, A and B). Further, we expressed KRASG12D or BRAFV600E allele in liver sinusoidal endothelial cells of developing murine embryos using Lyve1-Cre (KrasG12D F/+; Lyve1Cre and BrafV600E F/+; Lyve1Cre) mice (Fig. 2, C–E). These embryos exhibited complete embryonic lethality, reduced liver size, and fully penetrant hepatic vascular cavernomas phenotype without major changes in endothelial cell proliferation (Fig. 2, C–E; Fig. S2, E–H; and Fig. S3, E–H). KrasG12D F/+; Lyve1Cre and BrafV600E F/+; Lyve1Cre embryonic livers revealed hypoxia and apoptosis in clusters of nonendothelial hepatic cells near vascular cavernomas, thus limiting growth and expansion of fetal liver (Fig. 2, F and G). Tissue-resident macrophages express Lyve1; hence, we expressed KRASG12D or BRAFV600E allele in macrophages of developing murine embryos using LysM-Cre (KrasG12D F/+; LysMCre and BRAF-V600ELysM-Cre). KrasG12D F/+; LysMCre and BRAF-V600ELysM-Cre embryos appeared normal and displayed normal sinusoidal and hepatic development (Fig. S4, A and B). Similarly, tamoxifen-treated KrasG12D F/+; Cdh5CreERT2 and BrafV600E F/+; Cdh5CreERT2 mice did not reveal sinusoidal capillary dilation and vascular cavernomas in spleen, an organ with fenestrated endothelium (Fig. S4, C and D). Together, these results suggest that KRAS or BRAF gain-of-function mutations within sinusoidal endothelial cells cause hepatic vascular cavernomas.
Endothelial KrasG12D or BrafV600E gain-of-function mutations aberrantly enhance “zipper-like” contiguous expression of adherens junctional proteins at sinusoidal endothelial cell–cell contacts
To define transcriptional changes, we performed an Affymetrix Clariom D assay on embryonic BrafV600E F/+; Lyve1Cre livers (Fig. 3, A–C; and Fig. S5, A and B). Among 5,317 differentially regulated transcripts with National Center for Biotechnology Information identifiers, 62% were reduced, while 38% were increased in BrafV600E F/+; Lyve1Cre compared with control (Fig. 3, A–C; and Fig. S5 A). Gene set enrichment analysis showed enrichment in cell adhesion molecules, focal adhesion, tight junction, and MAPK signaling pathways (Fig. 3, A–C; and Fig. S5, A and B). These enriched categories contained almost exclusively up-regulated genes, suggesting that the BRAF gain-of-function mutation promotes endothelial cell–cell adhesion (Fig. 3, B and C). We observed a corresponding increase in protein expression of critical endothelial adherens junctional proteins Cdh5, Fms-related receptor tyrosine kinase 4 (Flt4), and Pecam1 and aberrant zipper-like contiguous expression pattern of Cdh5 in BrafV600E F/+; Lyve1Cre, KrasG12D F/+; Lyve1Cre, BrafV600E F/+; Cdh5CreERT2, KrasG12D F/+; Cdh5CreERT2, and human hemangioma livers (Fig. 1 C; Fig. 3, D and E; Fig. S3, A, B, E, and F; and Fig. S5, C and D), suggesting abnormal contiguous adhesion between the hepatic sinusoidal endothelial cells.
Pharmacological inhibition of BrafV600E-Map2k or genetic ablation of Mapk1 rescue hepatic vascular cavernomas in BrafV600E F/+; Lyve1Cre and KrasG12D F/+; Lyve1Cre mice
Gain-of-function mutations in KRAS or BRAF hyperactivate the downstream canonical MAPK kinase (MAP2K)–MAPK signaling pathway in endothelial cells (Al-Olabi et al., 2018; Nikolaev et al., 2018). To investigate the requirement of Map2k hyperactivation in the pathogenesis of hepatic vascular cavernomas, we treated BrafV600E F/+; Lyve1Cre murine embryos with pharmacological inhibitors of Map2k (trametinib) and BrafV600E (dabrafenib; Fig. 4 A). Inhibition of Map2k activity rescued hepatic vascular cavernomas, liver size, embryonic lethality, and aberrant zipper-like contiguous expression of sinusoidal endothelial cell adhesion proteins in developing BrafV600E F/+; Lyve1Cre murine embryos (Fig. 4 A). Similarly, genetic ablation of endothelial Mapk1 prevented embryonic lethality, hepatic vascular cavernomas, and zipper-like contiguous expression of adherens junctional proteins in BrafV600E F/+; Lyve1Cre and KrasG12D F/+; Lyve1Cre murine embryos (Fig. 4, B and C; and Fig. S5, E and F). Together, these results suggest that the hyperactivation of KRAS–BRAF–MAP2K–MAPK1 signaling pathway in sinusoidal endothelial cells enhances aberrant zipper-like cell–cell adhesion, promoting cavernous expansion of hepatic sinusoidal capillaries (Fig. 4 D).
This study reveals the first and unexpected causal link between KRAS or BRAF mutations and hepatic vascular cavernomas. Mutations in RAS, the first human oncogene, and BRAF genes drive ∼30% of human tumors (Alcantara et al., 2019; Hunter et al., 2015; Lu et al., 2016; Smith et al., 2013; Der et al., 1982; Zhong et al., 1999; Davies et al., 2002; Thomas et al., 2007), yet no reports describe RAS or RAF mutations in hepatic cavernous hemangiomas, the most common benign tumor of the liver. We identified somatic activating KRAS and/or BRAF mutations in human hepatic tissue samples with sporadic hemangioma lesion (Fig. 1, A and B; and Table 1). Low allele frequencies of KRAS and BRAF mosaic-activating mutations in these sporadic lesions are consistent with the cellular heterogeneity of vascular malformations (Nikolaev et al., 2018; Al-Olabi et al., 2018). This study is the first to model human hepatic cavernous vascular malformation in mice. Our murine models, for the first time, established that sinusoidal endothelial cells expressing gain-of-function KRASG12D or BRAFV600E mutations give rise to hepatic vascular cavernomas (Fig. 2 and Fig. S2). This completely penetrant phenotype distinguishes KRASG12D and BRAFV600E as a somatic driver-mutations in slow-flow hepatic vascular cavernomas. In contrast, known genetic causes of human slow-flow malformations, such as germline mutations in KRIT1, CCM2, PDCD10, GLMN, RASA1, or TEK, may require a somatic “second hit” to drive vascular malformations (Macmurdo et al., 2016; Plummer et al., 2004; Vikkula et al., 1996).
Our study highlights a nononcogenic role of KRASG12D and BRAFV600E mutations in hepatic vascular cavernomas. These oncogenic forms of KRAS and BRAF drive tumorigenic growth of human cancer cells (Dankort et al., 2007; Jackson et al., 2001). Here we demonstrate that KRASG12D or BRAFV600E mutations cause hepatic vascular cavernomas without modulating endothelial cell proliferation, suggesting a context-dependent role for KRAS and BRAF mutations in the sinusoidal endothelium (Fig. 1 C; Fig. S1, B and C; and Fig. S3). Consistent with this interpretation, somatic activating KRAS or BRAF mutations do not cause cancer in some tissues, such as brain, uterus, and spinal cord (Nikolaev et al., 2018; Al-Olabi et al., 2018). In addition, our study highlights overlapping phenotypic consequences of gain-of-function KRAS and BRAF mutations in hepatic vascular cavernomas. The Catalogue of Somatic Mutations in Cancer reveals ∼99% mutation frequency in KRAS at three amino acids, G12 (∼89%), G13 (∼9%), and Q61 (∼1%). Among these mutations, G12D (∼36%) is the most prevalent, followed by G12C (∼14%), G12V (∼23%), and G13D (∼7%; Waters and Der, 2018; Prior et al., 2012). Similarly, V600 mutation represents ∼50% of all BRAF mutations in human tumors, with striking differences in replaced amino acids, i.e., V600E (∼70–90%), V600R (∼1%), V600M (∼0.3%), and V600D (∼0.1%; Pollock et al., 2003; Davies et al., 2002; Wan et al., 2004). Consistent with this, human hepatic tissue samples with sporadic hemangioma lesions revealed G12D, G12S, G13D, and/or G13S mutations in KRAS and/or V600M mutation in BRAF (Fig. 1, A and B; and Table 1). Despite these variations at hotspot amino acids G12/13 or V600, pathological human liver samples exhibited remarkable phenotypic resemblance, such as malformed sinusoidal capillary vessels, cavernous spaces filled with erythrocytes and lined by LYVE1+/CDH5+/PECAM1+ endothelial cells, and normal endothelial cell proliferation (Fig. 1 C and Fig. S1, B and C). Consistent with the human phenotype, mice expressing KRASG12D or BRAFV600E in endothelial cells rapidly developed sinusoidal capillary dilation and hepatic vascular cavernomas without altering endothelial cell proliferation (Fig. 2, Fig. S2, and Fig. S3). These data and previous literature reports suggest structural, biochemical, and functional similarities among KRAS G12/13 and BRAF V600 mutations (Yao et al., 2015; Smith et al., 2013; Lu et al., 2016; Hunter et al., 2015; Hobbs et al., 2016; Hobbs and Der, 2019). For example, KRAS G12D, G12C, and G12V mutations exhibit similar crystal structure and hyperactivate RAS, leading to sustained activation of the RAF–MAP2K–MAPK1 signaling pathway (Li et al., 2018; Muñoz-Maldonado et al., 2019; Johnson et al., 2019; Ihle et al., 2012; Hobbs et al., 2016; Hunter et al., 2015; Smith et al., 2013; Lu et al., 2016). Similarly, both BRAF V600E and V600M mutations form homodimers in a RAS-dependent manner to activate the MAP2K–MAPK1 signaling pathway (Yao et al., 2015; Pratilas et al., 2009; Rajakulendran et al., 2009). These structural and biochemical studies have supported the development of pharmacologic inhibitors directly targeting KRAS G12/13 mutations, BRAF V600 mutations, MAP2K, and MAPK1 in human tumors (Ostrem and Shokat, 2016; Shi et al., 2012; Waters and Der, 2018; Mullard, 2019; Cox et al., 2014; Samatar and Poulikakos, 2014; Ross et al., 2017).
Here we identify the RAS–RAF–MAP2K–MAPK1 signaling pathway as a critical regulator of the size, shape, and branching pattern of hepatic sinusoidal capillary vessels. The formation and maintenance of a 3D network of sinusoidal capillaries with fine adaptations in branching and vessel diameter is essential for liver growth and function (Géraud et al., 2017; Si-Tayeb et al., 2010). The endothelial lining of sinusoidal capillaries exhibits a unique morphology, with fenestrations grouped into sieve plates, specialized junctional complexes, and an incomplete basement membrane, facilitating transfer of solutes and large molecules between the bloodstream and hepatocytes. Sinusoidal endothelial cells exhibit a dynamic cellular rearrangement during sinusoid patterning and maintenance (Zapotoczny et al., 2019); however, the cellular and molecular mechanisms underlying hepatic vascular cavernomas remain unknown. In developing murine embryos, “button-like” discontiguous expression of adherens junctional complexes, comprising Pecam1, Cdh5, Vegfr2, and Flt4, between the endothelial cells drives the patterning and polarized cellular migration during vascular morphogenesis (Szymborska and Gerhardt, 2018; Bentley et al., 2014). Our data demonstrate that constitutive activation of the KRAS–BRAF–MAP2K–MAPK1 signaling pathway in sinusoidal endothelial cells promotes zipper-like contiguous expression of adherens junctional proteins, such as Cdh5, Pecam1, and Flt4, switching hepatic sinusoidal capillaries from branching to cavernous expansion (Fig. 4 D). Our findings support a model in which discontiguous adhesion drives intercalation, and thus sinusoidal capillary branching, whereas contiguous adhesion disrupts intercalation, and thereby expansion of sinusoidal capillary diameter (Fig. 4 D). These data support prior observations that Pecam1, Cdh5, Flt4, and Vegfr2 complex, classically associated with endothelial cell proliferation, regulates integrity of endothelial cell–cell junctions (Tzima et al., 2005; Privratsky and Newman, 2014; Lampugnani et al., 2018; Dejana and Vestweber, 2013; Coon et al., 2015). For example, Pecam1 interacts with Cdh5 at endothelial cell–cell contacts to regulate junctional integrity and movement of endothelial cells (Privratsky and Newman, 2014; Tzima et al., 2005). Mice lacking Cdh5 exhibit embryonic lethality owing to defective Vegfr signaling and endothelial cell–cell interactions, suggesting a functional relationship between these receptors (Carmeliet et al., 1999). Consistent with this interpretation, mice with loss of Flt4 function show defective integrity of endothelial cell–cell junctions due to reduced junctional Cdh5 localization and increased Vegfr2 expression (Heinolainen et al., 2017; Tammela et al., 2008). Taken together, these findings demonstrate that activity of Flt4 and Vegfr2 is essential for Pecam1-Cdh5 adherens junction integrity at endothelial cell–cell contacts.
This is the first report, to our knowledge, of an effective molecular therapy for hepatic vascular cavernomas in mice. KRASG12D or BRAFV600E mutations drive constitutive activation of the classic MAP2K–MAPK1 signaling pathway (Dankort et al., 2007; Van Meter et al., 2007). Food and Drug Administration–approved drugs dabrafenib and trametinib inhibit BRAFV600E and MAP2K, respectively (Ribas and Flaherty, 2011), two kinases within the RAS–RAF–MAP2K–MAPK1 signaling pathway. In our murine models, concurrent administration of dabrafenib and trametinib inhibited the growth of BRAFV600E-dependent hepatic vascular cavernomas (Fig. 4 A). This is consistent with recent randomized, multicenter clinical trials (NCT01584648 and NCT01597908) demonstrating that dabrafenib plus trametinib improve progression-free survival and overall survival of patients with BRAFV600E-positive tumors (Robert et al., 2019). The present study takes this therapeutic approach one step further by revealing that Mapk1 functions as a critical downstream factor to promote KrasG12D- or BrafV600E-dependent hepatic vascular cavernomas. Genetic ablation of Mapk1 in BrafV600E F/+; Lyve1Cre and KrasG12D F/+; Lyve1Cre mice rescued the hepatic vascular cavernomas phenotype, suggesting Mapk1 as a novel therapeutic target (Fig. 4, B and C). Together, our murine models provide a platform to test novel pharmacological compounds to prevent formation and enlargement of RAS–MAPK1 signaling pathway–dependent hepatic vascular cavernomas.
Materials and methods
A total of 43 patients were initially included in the study after searching the database covering a period of 10 yr for diagnostic keywords “liver hemangioma.” For comparison, initially included as controls were a total of 10 patients whose liver tissue specimens were available in the tissue archive for non–vascular-related pathology and did not have known KRAS- or BRAF-driven cancer pathology. Tissue blocks containing liver hemangioma lesions and blocks containing liver from controls were retrieved, and a trained pathologist reviewed H&E-stained sections from the blocks. Patient liver specimens (n = 4) that did not contain hemangioma lesions or control liver samples affected with marked cirrhosis/fibrosis were excluded from further analysis. Review of clinical and pathological information and study of archived tissue samples were approved with waiver for informed consent (Health Insurance Portability and Accountability Act waiver) by the University of Massachusetts Medical School Institutional Review Board. Because of the large number of subjects and the time frame and retrospective nature of the study, the research could not practicably be performed without the waiver of informed consent. Patient data and specimens with or without preexisting pathology were collected for clinical reasons as per the Institutional Review Board approval. A data collection sheet was used to record data using a deidentified, separate numbering system. The data, including slides, were deidentified and given a separate number. The deidentified data are stored in locked investigator offices and password protected, with access to only authorized study personnel.
DNA extraction from formalin-fixed paraffin-embedded (FFPE) tissue
Sections of FFPE tissue blocks from patients and controls were used for DNA extraction; where necessary, patient sections were microdissected to improve lesion coverage. DNA extraction was performed using the QIAamp DNA FFPE Tissue Kit according to manufacturer’s instructions. Briefly, sections were deparaffinized by sequential washes in xylene and 100% ethanol. Samples were then centrifuged, air-dried, and incubated overnight in tissue lysis buffer and proteinase K. Finally, samples were loaded in QIAcube as per the instructions in the instrument manual, and DNA was eluted into the buffer provided with the QIAamp DNA FFPE Tissue Kit. Eluted DNA was quantified using nanodrop, aliquoted, and stored at −80°C until downstream analyses.
Next-generation sequencing was performed as previously described (Kamionek et al., 2016). 37 primer pairs were used to cover the known hotspot mutations in KRAS, NRAS, HRAS (codons 12, 13, 61, and 146), BRAF (codons 460–470 and 600), and the entire coding sequence of the MAP2K1 gene. For the purposes of the current study, only mutations in the hotspot codons of KRAS (codons 12 and 13) and BRAF (codon 600) were further analyzed in depth. Amplicon libraries were created from genomic DNA (10 ng), according to the manufacturer’s protocol (Ion AmpliSeq Library Kit 2.0; Thermo Fisher Scientific). The Ion CHEF Template System (Thermo Fisher Scientific) was used for emulsion PCR to amplify library DNA onto IonSphere Particles and loaded on 318v2 chips. Sequencing was performed on an Ion Torrent PGM (IC200 Sequencing Kit; Thermo Fisher Scientific) with coverage of ×5,000 to ×10,000. Raw sequencing data files were analyzed using two pipelines, Variant Caller (Life Technologies) and NextGene (SoftGenetics), and output was compared using a laboratory-developed visual basic Excel program. Variants with an allele frequency >0.1% detected at the hotspot codons of KRAS and BRAF from either pipeline were called positive. Further, all positive calls were manually verified using NextGene (SoftGenetics) alignment viewer.
Mutational analysis by PCR clamp assay and Sanger sequencing
KRAS codon 12/13 or BRAF codon 600 was analyzed using BNA clamp real-time quantitative PCR, and specific primer sets (IDT) were used (BRAF-F, 5′-GTAAAACGACGGCCAGTAAACTCTTCATAATGCTTGCTCTG-3′, BRAF-R, 5′-CAGGAAACAGCTATGACGGACCCACTCCATCGAGA-3′; KRAS-F, 5′-GTAAAACGACGGCCAGTGTACTGGTGGAGTATTTGATAGTG-3′; and KRAS-R, 5′-CAGGAAACAGCTATGACATCGTCAAGGCACTCTTGCCTAC-3′). M13F or M13R sequencing primers are underlined. SsoFast Evagreen Supermix (Bio-Rad) and DNA in a total volume of 20 µl were used with or without a BNA clamp (BRAF BNA, 5′-GATTTCACTGTAGC-spacerC3-3′, and KRAS BNA, 5′-mUA+CGCCACCAGCT-spacerC3-3′; BioSynthesis). A Bio-Rad CFX96 real-time thermal cycler was used. The PCR conditions were as follows: 98°C denaturation for 2 min, then 60 cycles of 60°C for 5 s, and 98°C for 5 s. Briefly, real-time quantitative PCR was performed with and without a BNA clamp designed to block amplification of wild-type sequences. The delta threshold count (ΔCT) was calculated from the difference between the Ct values obtained from PCR performed without clamp and with clamp. The RKO human colon cell line, with c. 1799T>A (p.V600E) BRAF mutation, was used as positive control (ATCC). Specimens with a ΔCT ≥ the negative control were deemed K-RAS and/or BRAF wild type. The remaining specimens were regarded as positive, and mutation status was confirmed by BigDye Sanger methodology (Life Technologies) using capillary gel electrophoresis (ABI 3500xL). Results were interpreted using Mutation Surveyor V4.0.8 (SoftGenetics). This assay detects all known mutations in KRAS codons 12/13 and between BRAF codons 598 and 602 with a detection limit of <0.1%.
KrasG12D, BRAFV600E, Mapk1Flox, Lyve1Cre, LysMCre, and Gt(ROSA)26Sor-LacZ reporter mice and Cdh5CreERT2 have been described previously (Clausen et al., 1999; Dankort et al., 2007; Jackson et al., 2001; Pham et al., 2010; Samuels et al., 2008; Sörensen et al., 2009; Soriano, 1999). All mice were maintained on a mixed genetic background. Mice were backcrossed for at least five generations to reference strains, and littermates were used as controls for all experiments. Neonatal mice were administered tamoxifen once daily orally on postnatal days 0, 1, and 2 or postnatal days 3, 4, and 5 at a dose of 50 µg/d, and all mice were analyzed at postnatal day 12. Adult mice were treated with tamoxifen once daily for 3 d consecutively at 6 wk of age with a dose of 1 mg/10 g body weight, and mice were euthanized at 16 wk (KrasG12D F/+; Cdh5CreERT2/+) or 8–9 wk (BrafV600E F/+; Cdh5CreERT2/+) of age for tissue analyses. For pharmacologic rescue experiments, pregnant mice were treated once daily by oral gavage with either vehicle (0.2% Tween 80 and 0.5% hydroxymethyl cellulose dissolved in water) or trametinib (0.5 µg/g body weight/d) and dabrafenib (25 µg/g body weight/d) dissolved in vehicle solution for 5 consecutive days from embryonic days 9.5 to 13.5. University of Massachusetts Medical School is accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International and follows the Public Health Service Policy for the Care and Use of Laboratory Animals. Animal care was provided in accordance with the procedures outlined in the Guide for the Care and Use of Laboratory Animals (National Research Council Committee for the Update of the Guide for the Care and Use of Laboratory Animals, 2011). University of Massachusetts Medical School Institutional Animal Care and Use Committee (Worcester, MA) approved all animal use protocols.
Murine embryos and tissues were fixed in 2% paraformaldehyde at 4°C overnight, ethanol dehydrated, embedded in paraffin, and sectioned at 8-µm thickness using a Leica fully motorized rotary microtome. FFPE tissue sections were deparaffinized in xylene and rehydrated through an ethanol gradient, followed by 2-min Harris-modified hematoxylin and 30-s eosin-Y staining. Slides were dehydrated with ethanol, cleared with xylene, and mounted with Vectashield glass mounting medium.
Immunohistochemistry was performed as previously described (Janardhan et al., 2017). Briefly, slides with sections were deparaffinized and immersed in sodium citrate buffer (10 mM sodium citrate and 0.05% Tween-20, pH 6) and placed in a 2100 Antigen Retriever (Aptum Biologics) for heat-induced antigen retrieval. After antigen retrieval, sections were blocked in 10% donkey serum, 0.1% BSA, and 0.3% Triton X-100 in PBS for 1 h at room temperature. Sections were then washed in PBS and incubated with primary antibodies (Table S1) in 10% donkey serum and PBS overnight at 4°C. Finally, slides were washed in PBS, incubated in Alexa Fluor 488 or 546/568 conjugated secondary antibodies (1:500) with Hoechst (1:1,000) for 1 h at room temperature, rinsed in PBS, and mounted with Vectashield glass mounting medium.
Images of dissected embryos, mice, and tissue samples were captured using a Leica MZ10 F fluorescence stereomicroscope equipped with a 0.7× C-mount, Achromat 1.0 × 90-mm objective, a SOLA light engine, a DS-Fi1 color camera (Nikon), and NIS-Elements Basic Research software (Nikon). H&E-stained section images were captured using a Nikon Eclipse 80i microscope equipped with CFI Plan Fluor 4×/10×/20×/40× objective lenses, a DS-Fi1 color camera, and NIS-Elements Basic Research software. Immunofluorescent-stained slides were imaged using Plan-Apochromat objectives with DIC (63×/1.4 oil, and 20×/0.8) on Zeiss LSM800 Airyscan inverted digital spectral confocal microscope equipped with Definite Focus 2.0, three confocal GaAsP detectors including Airyscan detector with 6,000 × 6,000-pixel resolution, and solid-state laser module with 405-, 488-, 561-, and 640-nm beam splitter. Image stacks of vertical projections were assembled using Zeiss Zen 2.5 imaging software.
Livers from embryonic day 11.5 embryos were dissected, cleared of blood, snap frozen in liquid nitrogen, and stored at −80°C until further processing. After genotyping, eight liver samples each from control and mutant embryos were respectively pooled before RNA isolation. Total RNA was isolated using Qiagen RNAeasy microkit as per the manufacturer’s instructions. Eluted RNA was then aliquoted into three tubes each for control and mutants. All samples were checked for quality using Agilent 2100 Bioanalyzer and had an RNA integrity number >9.5. Microarray analysis was performed by the University of Massachusetts Genomics Core Facility using Mouse Clariom D array (Thermo Fisher Scientific). Raw microarray .cel files were imported and processed using Transcriptome Analysis Console version 4.0. Full microarray datasets can be accessed at the Gene Expression Omnibus under accession no. GSE147848.
Statistical significance was determined based on a two-tailed Student’s t test or χ2 test (GraphPad Prism 8.0). A P value < 0.05 was considered significant.
Online supplemental material
Fig. S1 shows hepatic vascular cavernomas with normal endothelial cell proliferation in human patients with somatic gain-of-function mutations in KRAS or BRAF genes. Fig. S2 shows hepatic vascular cavernomas and embryonic lethality in mice with endothelial KrasG12D or BrafV600E gain-of-function mutations. Fig. S3 shows normal proliferation of hepatic endothelial cells in mice with endothelial KrasG12D or BrafV600E gain-of-function mutations. Fig. S4 shows normal sinusoidal and hepatic development in mice with macrophage KrasG12D or BrafV600E gain-of-function mutations. In addition, mice with endothelial KrasG12D or BrafV600E gain-of-function mutations show normal spleen morphology and sinusoids. Fig. S5 shows increased expression of adherens junctional proteins in mice with endothelial KrasG12D or BrafV600E gain-of-function mutations. In addition, genetic ablation of endothelial Mapk1 rescues lethality at embryonic day 14.5 in embryos with endothelial KrasG12D or BrafV600E gain-of-function mutations. Table S1 lists the antibodies used.
We gratefully acknowledge Benjamin Chen, MD, PhD (University of Massachusetts Medical School, Worcester, MA) for his help with examination of hemangioma patient samples. We thank the University of Massachusetts Medical School Genomics Core Facility for help with the microarray experiment.
This work was supported by National Institutes of Health grant HL141377 (to C.M. Trivedi).
Author contributions: H.P. Janardhan and C.M. Trivedi conceived the study, designed experiments, and wrote the manuscript. H.P. Janardhan, X. Meng, K. Dresser, L. Hutchinson, and C.M. Trivedi performed experiments and acquired and analyzed data. All authors reviewed results and approved the final version of the manuscript.
Disclosures: The authors declare no competing interests exist.