Although new drug discoveries are revolutionizing cancer treatments, repurposing existing drugs would accelerate the timeline and lower the cost for bringing treatments to cancer patients. Our goal was to repurpose CPI211, a potent and selective antagonist of the thromboxane A2-prostanoid receptor (TPr), a G-protein–coupled receptor that regulates coagulation, blood pressure, and cardiovascular homeostasis. To identify potential new clinical indications for CPI211, we performed a phenome-wide association study (PheWAS) of the gene encoding TPr, TBXA2R, using robust deidentified health records and matched genomic data from more than 29,000 patients. Specifically, PheWAS was used to identify clinical manifestations correlating with a TBXA2R single-nucleotide polymorphism (rs200445019), which generates a T399A substitution within TPr that enhances TPr signaling. Previous studies have correlated 200445019 with chronic venous hypertension, which was recapitulated by this PheWAS analysis. Unexpectedly, PheWAS uncovered an rs200445019 correlation with cancer metastasis across several cancer types. When tested in several mouse models of metastasis, TPr inhibition using CPI211 potently blocked spontaneous metastasis from primary tumors, without affecting tumor cell proliferation, motility, or tumor growth. Further, metastasis following intravenous tumor cell delivery was blocked in mice treated with CPI211. Interestingly, TPr signaling in vascular endothelial cells induced VE-cadherin internalization, diminished endothelial barrier function, and enhanced transendothelial migration by tumor cells, phenotypes that were decreased by CPI211. These studies provide evidence that TPr signaling promotes cancer metastasis, supporting the study of TPr inhibitors as antimetastatic agents and highlighting the use of PheWAS as an approach to accelerate drug repurposing.

This article is featured in Highlights of This Issue, p. 2407

The use of deidentified human genetic data tied to robust electronic health records (EHR) enables discovery of human genotypes correlating with specific clinical conditions. This approach, known as phenome-wide association study (PheWAS; refs. 1, 2), can discover clinical “phenotypic” indications associating with a genotypic variation, enabling discoveries of novel pathophysiologic gene functions. From a translational perspective, PheWAS could predict potential opportunities for repurposing an existing catalog of molecularly targeted therapies (3), diminishing the risk, cost, and time needed for bringing treatments to patients.

To test this principle, we used PheWAS to investigate repurposing opportunities for the drug CPI211 (4–6), a potent and highly selective small-molecule inhibitor of thromboxane A2 and prostanoid receptor (TPr). CPI211 (PubChem ID 3037233) is safe and well tolerated, as shown in a randomized, dose-escalating (≤1,000 mg), placebo-controlled clinical trial (7). Decades of basic and translational science determined that TPr, expressed in platelets, endothelial cells (EC), and smooth muscle cells, regulates different aspects of cardiovascular homeostasis. TPr dysregulation contributes to venous hypertension, pulmonary hypertension, thrombosis, and asthma, steering the development of CPI211 and other TPr antagonists toward these clinical indications (8, 9). However, aberrant activation of TPr may have yet undiscovered pathologic consequences for which TPr inhibitors might be effective. Because pathologically increased TPr expression can be caused by naturally occurring single-nucleotide polymorphisms (SNP) within the TPr gene, TBXA2R, we propose that PheWAS analysis using these previously identified and characterized TBXA2R SNPs as a genomic proxy for increased TPr expression could identify clinical manifestations correlating with increased TPr expression (10).

Several previously described SNPs within TBXA2R are known to alter TPr expression, some that diminish TPr and thus correlate with decreased platelet function (11), and others that increase TPr expression. As might be expected, the SNPs that increase TPr expression are known to correlate with advanced cardiovascular disease (12–14) or asthma (15). We investigated one TBXA2R SNP that increases TPr expression, rs200445019. This SNP was selected based on its occurrence within the TBXA2R coding region, generating a threonine-to-alanine alteration at residue 399 (T399A) within the C-terminus of the TPr splice variant, TPr-β (https://www.ncbi.nlm.nih.gov/snp/rs200445019/), the domain of TPr-β driving TPr desensitization, internalization, trafficking, and proteosomal degradation (16–19). Following ligand activation, TPr signaling activates protein kinase C (PKC), a serine threonine kinase that then phosphorylates TPr-β T399, inducing TPr internalization and degradation. TPr T399 mutation causes loss of PKC-mediated TPr phosphorylation, and loss of ligand-induced receptor downregulation. Thus, the TPr-β T399A substitution generated by rs200445019 results in greater TPr expression and signaling, supporting use of rs200445019 as a genomic proxy for increased TPr in PheWAS analyses. Here, we report that PheWAS analysis of rs200445019 identified a novel and unexpected correlation between rs200445019 and metastatic cancer diagnoses, and that the TPr inhibitor CPI211 blocked metastasis in multiple cancer models.

The BioVU biorepository at Vanderbilt University Medical Center (VUMC) contains >250,000 deidentified DNA samples extracted from excess patient blood samples collected during routine clinical testing, linked to corresponding, longitudinal clinical and demographic data derived from Synthetic Derivative, a deidentified EHR built for research purposes (20, 21). BioVU data were utilized in accordance with VUMC IRB# 151121. The allelic variant of focus was the TBXA2R missense SNP rs200445019, selected based on the presence of meaningful validation signals in the PheWAS results to support inference of variant effects in vivo, including hypertension and cardiovascular indications, phenotypes expected to be observed if TPr function was enhanced.

PheWAS using previously reported methods (22) focused on 29,722 patients of European ancestry genotyped using the Illumina Infinium Exomechip (Illumina, Inc.), containing ∼250,000 SNPs across genomic protein coding regions discovered through exome- and whole-genome sequencing in more than 12,000 individuals. All distinct phenotypic data (ICD9 billing codes representing approximately 1,660 phenotypes) were captured from EHRs and translated into corresponding phenotype groupings. Cases for each phenotype (condition) were defined as patients having ≥ 2 ICD9 codes documented on separate days and mapping to a specific PheCode. The PheCode ‘malignant secondary neoplasm’ (PheCode 198) included ICD9 codes for secondary malignant neoplasm of respiratory and digestive systems, secondary malignant neoplasms of other specified sites, malignant pleural effusion, and several others. Controls were defined as all patients who did not have any of the ICD9 codes defining a case. Only phenotypes occurring in ≥25 cases (0.42% of genotyped patients) were included. Physicians at VUMC with relevant qualifications and content-specific knowledge performed evidence reviews in cases of secondary malignancy. The PheWAS algorithm (22) was applied to calculate case and control genotype distribution, χ2 distribution, associated allelic P value and allelic odds ratio (OR). For χ2 distribution cell counts below five, Fisher exact test was used to calculate P value using R software (http://www.r-project.org/).

CPI211

CPI211 (Supplementary Fig. S1A, PubChem ID 3037233) was generously provided by Cumberland Pharmaceuticals, Inc. The chemical formula of CPI211 is C25H32N2O5, with a molecular weight of 440.5 g/mol. CPI211 is also known by the following aliases: Ifetroban, 143443-90-7, UNII-E833KT807K, and BMS-180291.

Cell culture

Cells purchased from American Tissue Type Collection [4T1: (CRL-2539) HEK-293T/17 (CRL-11268), MDA-MB-231 (HTB-26), MiaPaCa-2 (CRL-1420), A549 (CRM-CCL-185)] were maintained in growth medium defined by supplier. Cells constitutively expressing luciferase were generated by lentiviral particle (LVP) transduction with CAG-Luciferase LVP (GenTarget Inc.) followed by puromycin selection. LVP encoding human TPr-β and TPr-β(T399A) were custom synthesized in Lentifect (GeneCopoeia). 293T cells were transduced with empty LVP or those expressing TPrβ(WT) and TPr-β(T399A), followed by blasticidin selection. Mouse pulmonary microvascular ECs (MPMEC) were harvested and cultured as described (23). For MPMEC growth assays, cells (0.5 × 104) were seeded on gelatin-coated 35-mm dishes, and cultured for 7 days in growth media, U46619 (0.5 μmol/L), CPI211 (0.1 μmol/L), or 4T1 conditioned media (serum-free media conditioned for 48-hour culture of 4T1 cells, passed through a 0.2-μm syringe filter). Cells were stained with 0.01% crystal violet, rinsed with water and dried, and then imaged on a flatbed scanner. Human umbilical vein ECs (HUVEC) were purchased from Lonza and maintained in EGM-2 medium (Lonza). All cells used were cultured for <15 passages.

Platelet preparation

Human platelets were collected from blood donated by healthy anonymous volunteers, and prepared as previously described (24). The platelet pellet was suspended in Tyrode's buffer containing 0.1% bovine serum albumin and counted on a Beckman Z1 Coulter counter.

Western blotting

Cells were homogenized in lysis buffer [50 mmol/L Tris pH 7.4, 100 mmol/L NaF, 400 mmol/L NaCl, 0.5% NP-40, 100 μmol/L Na3VO4, 1X protease inhibitor cocktail (Roche)] and cleared by centrifugation. Platelets were homogenized with 0.6 N perchloric acid and cleared by centrifugation. For immunoprecipitation, cell lysates (1000 μg) cleared with protein A/G+ agarose were incubated with anti-VE-cadherin antibody (5 μg) and protein A/G+ agarose. Samples were resolved on 4%–12% polyacrylamide gels (Novex), transferred to nitrocellulose membranes (iBlot), and Western analysis was performed as described (25) using the following primary antibodies: [α-actin (Sigma-Aldrich); TPr (Abcam); from Cell Signaling Technologies: protein kinase C phospho-substrate, phospho-myosin light chain (MLC), MLC, Vascular Endothelial (VE)-Cadherin, phospho-tyrosine horseradish peroxidase (PTyr-100)].

Animals

All animals were housed under pathogen-free conditions. Experiments were performed in accordance with AAALAC guidelines and Vanderbilt University Institutional Animal Care and Use Committee approval. Where indicated, CPI211 was delivered daily at 50 mg/kg by oral gavage in 25 μL vehicle (4% sucrose in sterile water). 1 × 106 tumor cells (4T1, MDA-MB-231, A549, and MiaPaCa2) were injected via tail vein of athymic (nu/nu) Balb/C female mice. 4T1 tumor cells (1 × 106) were injected into the left inguinal mammary fat pad of athymic (nu/nu) Balb/C or WT Balb/C female mice. Tumor dimensions were measured using calipers. Tumor volume was calculated as: volume = length × width2 × 0.52). For bioluminescence imaging, mice were injected intraperitoneally with D-luciferin monosodium substrate (Thermo Fisher, 150 mg/kg) and imaged on IVIS Lumina III (Xenogen Corporation) 30 minutes after injection.

Histologic analysis

Tissue processing, staining, and immunohistochemistry (IHC) were performed by the Vanderbilt Translational Pathology Shared Resource. IHC on 5-μm paraffin-embedded sections was performed as described previously (26) using anti-CD31 (Abcam) and anti-Ki67 (Santa Cruz Biotechnologies). Photomicrographs were acquired on an Olympus CK40 inverted microscope through an Optronics DEI-750C camera using CellSens software for capture and morphometric analyses (27).

Transendothelial cell migration assay

MPMECs or HUVECs (1 × 105) were seeded on Matrigel-coated transwell lower surfaces and cultured 1 week before labeling with CellTracker Red (Molecular Probes; 0.5 mg/mL), transferred to 24-well dishes in 2% serum ± U46619 (0.5 μmol/L) and/or CPI211 (0.1 μmol/L) for 6 hours. 4T1-GFP cells (2 × 105) were seeded into the upper chamber for 16 hours, swabbed from upper transwell surfaces and then visualized on the lower filter surface using fluorescence microscopy.

Statistical analysis

Treatment groups were compared using either two-tailed Student t test or one-way ANOVA test coupled with Tukey means comparison test, where a P value < 0.05 was deemed representative of a significant difference between treatment groups.

Unbiased PheWAS correlates a naturally occurring TBXA2R SNP with metastasis

The TPr-T399A alteration generated by rs200445019 reportedly impairs TPr desensitization and degradation following ligand-induced stimulation. To confirm this, we expressed wild-type TPr-β (TPrβ-WT) and TPr-β T399A (TPrβ-T399A) in 293T cells, which express low endogenous TPr levels (Supplementary Fig. S1B). Under serum-starved conditions, similar protein expression of TPr was seen in cells expressing TPrβ-WT and TPrβ-T399A. As expected, addition of the TPr agonist U46619 (0.5 μmol/L) increased PKC activity and MLC phosphorylation (Fig. 1A), which occurred at similar levels in cells expressing TPrβ-WT and TPrβ-T399A. However, TPrβ-WT protein levels diminished within 5 minutes of U46619 treatment, whereas TPrβ-T399A levels remained unchanged, consistent with previous reports that TPrβ T399 drives ligand-induced TPr downregulation. This was confirmed by treating cells with U46619 for 5 minutes, then chasing with U46619-free media for up to 60 minutes, revealing a sustained downregulation of TPrβ-WT, but persistently elevated TPrβ-T399A (Fig. 1B). These findings confirm use of rs200445019 as a proxy for increased TPr expression.

Figure 1.

PheWAS and expression analyses correlate TBXA2R with metastasis and decreased progression-free survival. A and B, Western analysis of serum-starved cells treated 5 minutes with U46619 (0.5 μmol/L) using the antibodies shown to the left of each panel. U46619-treated cells were washed after 5 minutes, and then cultured up to 60 minutes with serum-free media (B). C, Distribution of cancer types producing secondary malignancies in PheWAS of >29,000 patients. D, Kaplan–Meier analysis of METABRIC-curated breast cancers assessing the relationship between TBXA2R and patient OS using cBio software (left). Cutoff between low/high: 2 SD above the average TBXA2R. Kaplan–Meier meta-analysis of NSCLC (center) and gastric cancer (right) gene-expression array data assessing the relationship between PFS and TBXA2R using Kmplot software. Cutoff between low/high: median TBXA2R expression. E, Western analysis of cells treated with CPI211 (100 nmol/L) and U46619 (0.5 μmol/L). F, Western analysis of washed human platelets treated ex vivo with CPI211 (0.1 nmol/L) and U46619 (0.5 μmol/L). G, Aggregometry of washed human platelets treated ex vivo with U46619 (0.5 μmol/L) ± CPI211 (5–500 nmol/L).

Figure 1.

PheWAS and expression analyses correlate TBXA2R with metastasis and decreased progression-free survival. A and B, Western analysis of serum-starved cells treated 5 minutes with U46619 (0.5 μmol/L) using the antibodies shown to the left of each panel. U46619-treated cells were washed after 5 minutes, and then cultured up to 60 minutes with serum-free media (B). C, Distribution of cancer types producing secondary malignancies in PheWAS of >29,000 patients. D, Kaplan–Meier analysis of METABRIC-curated breast cancers assessing the relationship between TBXA2R and patient OS using cBio software (left). Cutoff between low/high: 2 SD above the average TBXA2R. Kaplan–Meier meta-analysis of NSCLC (center) and gastric cancer (right) gene-expression array data assessing the relationship between PFS and TBXA2R using Kmplot software. Cutoff between low/high: median TBXA2R expression. E, Western analysis of cells treated with CPI211 (100 nmol/L) and U46619 (0.5 μmol/L). F, Western analysis of washed human platelets treated ex vivo with CPI211 (0.1 nmol/L) and U46619 (0.5 μmol/L). G, Aggregometry of washed human platelets treated ex vivo with U46619 (0.5 μmol/L) ± CPI211 (5–500 nmol/L).

Close modal

To detect phenotypes associated with rs200445019, we leveraged the BioVU database (21) to query 29,722 deidentified patient EHRs and their associated genotyping, scanning for novel associations between PheWAS codes and rs200445019. As expected, our unbiased PheWAS approach identified associations of rs200445019 with chronic venous hypertension, pulmonary heart disease, and primary pulmonary hypertension (Table 1; ref. 14). Unexpectedly, PheWAS also identified a novel association between rs200445019 and metastatic disease (coded as Secondary Malignancy) at multiple tissue sites, including lymph nodes, respiratory organs, digestive systems, brain/spine, and others. Chart review by medical experts confirmed the diagnoses of secondary malignancies, and revealed that these metastatic lesions were derived from multiple solid tumor types, including breast, colon, lung, head and neck, renal, gastric, ovarian cancers, melanomas, and others (Fig. 1C). However, PheWAS did not detect a specific correlation of rs200445019 with any primary cancer diagnosis, per se, suggesting that increased TPr expression and/or activity might have a selective role in tumor metastasis.

Table 1.

PheWAS of 29,722 patients of European ancestry in BioVU genotyped using the Illumina Infinium Exomechip and phenotyped for approximately 1,660 ICD9 codes show expected associations of rs200445019 with hypertension, pulmonary hypertension, and pulmonary heart disease.

ConditionPheWAS codeP valuesORCase carriersTotal casesControls
Chronic venous hypertension (validation phenotype) 453 9.1E−06 31.89 15 22,276 
Pulmonary heart disease 415 0.003 2.11 18 1,512 24,713 
Primary pulmonary hypertension 415.21 0.01 3.28 261 24,713 
Secondary malignant neoplasm 198 3.8E−03 1.95 32 3,568 20,663 
Secondary malignancy of respiratory organs 198.2 6.5E−03 2.42 13 1,190 20,663 
Secondary malignancy of lymph nodes 198.1 7.0E−03 2.20 17 1,698 20,663 
Secondary malignant neoplasm of digestive system 198.3 1.2E−02 2.85 520 20,663 
Secondary malignancy of brain/spine 198.5 2.1E−02 2.68 587 20,663 
ConditionPheWAS codeP valuesORCase carriersTotal casesControls
Chronic venous hypertension (validation phenotype) 453 9.1E−06 31.89 15 22,276 
Pulmonary heart disease 415 0.003 2.11 18 1,512 24,713 
Primary pulmonary hypertension 415.21 0.01 3.28 261 24,713 
Secondary malignant neoplasm 198 3.8E−03 1.95 32 3,568 20,663 
Secondary malignancy of respiratory organs 198.2 6.5E−03 2.42 13 1,190 20,663 
Secondary malignancy of lymph nodes 198.1 7.0E−03 2.20 17 1,698 20,663 
Secondary malignant neoplasm of digestive system 198.3 1.2E−02 2.85 520 20,663 
Secondary malignancy of brain/spine 198.5 2.1E−02 2.68 587 20,663 

Note: Novel associates of rs200445019 with secondary malignancy (i.e., metastasis) were identified. A case was defined as an EHR with two or more ICD9 codes mapping to one of the phenotypes used for this analysis. A case carrier is a defined as a case testing positive for rs200445019. Patients with records lacking the identified ICD9 codes were categorized as controls. The PheWAS algorithm was applied to calculate case and control genotype distribution, the χ2 distribution, associated allelic P value, and allelic OR. For those χ2 distributions in which observed cell counts fell below five, Fisher exact test was used to calculate the P value using the R statistical package (http://www.r-project.org/). Only phenotypes occurring in a minimum of 25 cases (0.42% of genotyped patients) were included in the analysis.

Next, we assessed the relationship between TBXA2R expression and cancer patient outcome. Analysis of METABRIC-curated invasive breast cancers assessed by gene-expression array (N = 2,051; ref. 28) revealed decreased overall survival (OS) of patients whose tumors harbored high TBXA2R (Fig. 1D). High TBXA2R expression also correlated with decreased OS in lung adenocarcinomas and gastric cancers (29). We extended this analysis across 14 additional cancer types, using publicly available tumor RNA-seq data tied to disease-free progression (DFP) outcomes (ref. 29; Supplementary Table S1), finding that 12 of 15 solid cancer types displayed an inverse correlation between TBXA2R expression and DFP, with strong statistical significance in 11 of these, including breast cancer, ovarian cancer, renal papillary cell carcinoma, sarcoma, and others. Interestingly, 2 of 14 cancer types assessed showed increased DFP in patients expressing increased TBXA2R [bladder cancer and hepatocellular carcinoma (HCC)]. The reasons underlying why bladder cancer and HCC might show improved survival upon increased TBXA2R expression are unclear, but may relate to unique pathophysiologies of these cancers, routes of metastatic spread, or tissue-selective roles of TPr. For the 12 cancers demonstrating a correlation between TBXA2R expression and worse outcome, these findings support the notion that TBXA2R/TPr may have an underappreciated role in tumor metastasis, and that targeted TPr inhibition might suppress metastasis.

To confirm TPr inhibition by the TPr antagonist CPI211, we assessed 293T cells expressing exogenous TPr-β. Western blot analysis confirmed that CPI211 blocked MLC phosphorylation in U46619-treated cells expressing TPrβ-WT and TPrβ-T399A (Fig. 1E). We also confirmed TPr-α inhibition by CPI211 using washed human platelets, which almost exclusively express TPr-α. U46619-induced P-MLC was blocked robustly by CPI211 (Fig. 1F). Dilutions of CPI211 down to 5 nmol/L blocked U46619-induced platelet aggregation (Fig. 1G), confirming potent TPrα inhibition by CPI211.

TPr inhibition reduces cancer metastasis in several preclinical cancer models

We assessed the impact of TPr inhibition on hematogenous metastasis using tumor cells [4T1 (mouse mammary cancer), MDA-MB-231 (human breast cancer), MiaPaCa2 (human pancreatic cancer), and A549 (human lung cancer)] delivered directly into circulation by intravenous (i.v.) injection. Western blot analysis revealed only very low TPr expression in each tumor cell line, and modest or no P-MLC in U46619-treated cells (Fig. 2A). Athymic (nu/nu) mice were pretreated with CPI211 (50 mg/kg) before tumor cell delivery, continuing treatment for the following 28 days (Fig. 2B). 4T1 lung metastases were seen in 100% of vehicle-treated mice, but only 70% of CPI211-treated mice (N = 10; Fig. 2C). Similarly, CPI211 treatment decreased the percentage of mice harboring MDA-MB-231 lung metastases from 90% to 20% (N = 10), and mice with A549 lung metastases from 60% to 10% (N = 10). Although 90% of mice inoculated with MiaPaca2 cells developed lung lesions regardless of vehicle or CPI211 treatment, the average number of MiaPaca2 metastases per mouse was diminished in CPI211-treated mice by more than half. These results suggest that TPr inhibition decreases hematogenous metastasis of multiple cancer types.

Figure 2.

The TPr inhibitor CPI211 decreases metastasis in mouse metastasis models. A, Western blot analysis of lysates collected from cells treated with U46619 (0.5 μmol/L). B, Schematic timeline used for treating mice with CPI211 (50 mg/kg) for 2 days prior to, through 28 days after, tail-vein delivery of tumor cells. C, Lungs were assessed histologically for metastases. Values shown are the number of lung lesions per mouse. Midlines are average, and error bars are SD. P values, Student t test. D, Treatment schematic for 4T1-Luc tumor–bearing mice with CPI211 for 28 days. Primary tumors were resected on treatment day 12. E, Left: average lung bioluminescence was determined weekly by IVIS. Each bar shows average bioluminescence per group. Error bars, SD. Center: raw bioluminescence (photons/second) for each mouse on treatment days 14 and 28. Right: values shown are number of lung metastases in each mouse on day 28. Midline represents the average; error bars represent SD. n.s., not significant.

Figure 2.

The TPr inhibitor CPI211 decreases metastasis in mouse metastasis models. A, Western blot analysis of lysates collected from cells treated with U46619 (0.5 μmol/L). B, Schematic timeline used for treating mice with CPI211 (50 mg/kg) for 2 days prior to, through 28 days after, tail-vein delivery of tumor cells. C, Lungs were assessed histologically for metastases. Values shown are the number of lung lesions per mouse. Midlines are average, and error bars are SD. P values, Student t test. D, Treatment schematic for 4T1-Luc tumor–bearing mice with CPI211 for 28 days. Primary tumors were resected on treatment day 12. E, Left: average lung bioluminescence was determined weekly by IVIS. Each bar shows average bioluminescence per group. Error bars, SD. Center: raw bioluminescence (photons/second) for each mouse on treatment days 14 and 28. Right: values shown are number of lung metastases in each mouse on day 28. Midline represents the average; error bars represent SD. n.s., not significant.

Close modal

Because venous tumor cell delivery allows tumor cells to bypass steps of intravasation, we used a rigorous model of spontaneous metastasis from primary orthotopic 4T1 mouse mammary tumors. Mice harboring luciferase-expressing 4T1 tumors were randomized into treatment groups receiving vehicle or CPI211 when tumors reached 50 mm3 (Fig. 2D). Mice were treated daily for 28 days, although tumors were surgically resected on treatment day 12. By treatment day 21, luminescence was detected in the thoracic region of vehicle-treated mice, but not in CPI211-treated mice (Fig. 2E). By day 28, average luminescence in vehicle-treated mice was approximately 18-fold higher than in CPI211-treated mice. Notably, 3 of 10 mice treated with CPI211 had no measurable luminescence above baseline (Fig. 2E, right). Histologic examination of lungs confirmed reduced lung metastases in CPI211-treated mice versus vehicle-treated mice, including an absence of metastases in 3 of 10 CPI211-treated mice.

TPr inhibition did not affect tumor cell proliferation

We measured volume of tumors grown in athymic mice and treated for 12 days with CPI211 beginning when tumors were 50 mm3, revealing no difference in tumor size between vehicle and CPI211-treated groups (Fig. 3A). IHC for the cellular proliferation marker Ki67 revealed similar Ki67 staining in tumors harvested from vehicle-treated and CPI211-treated mice (Fig. 3B and C). These findings were supported by measurements of tumor cell growth in culture, showing that neither TPr activation nor TPr inhibition affected growth of 4T1 or MDA-MB-231 cells (Fig. 3D).

Figure 3.

The TPr inhibitor CPI211 does not affect primary tumor growth but decreases tumor vessels. 4T1 mouse mammary tumors were grown in athymic mice until reaching 50 mm3. Mice were treated daily with CPI211 (50 mg/kg) for 12 days. A, Tumor volume was measured on treatment day 12. Each data point is the tumor volume in a single mouse. Midlines, average tumor volume; error bars, SD. B, Tumor sections were stained with hematoxylin and eosin (H&E) or by IHC to detect Ki67. C, The percentage of Ki67+ tumor cells was quantitated in 5 random fields per sample. Midline, average of biological replicates. D, 4T1 and MDA-MB-231 cells were treated with CPI211 and/or U46619 in serum-free media for 3 days. 10% serum was included as a positive control for growth. Value points represent the average of three technical replicates ± SD; N = 3. E, IHC to detect CD31+ vessels in 4T1 tumors. Inset represents the area shown in higher power. F and G, Quantitation of total CD31+ vessel area (F) and number of CD31+ vessels with a lumen area > 50 μm2 (G). Data points, average value measured in three random fields per section; N, 4. Midlines, average ± SD. n.s., not significant.

Figure 3.

The TPr inhibitor CPI211 does not affect primary tumor growth but decreases tumor vessels. 4T1 mouse mammary tumors were grown in athymic mice until reaching 50 mm3. Mice were treated daily with CPI211 (50 mg/kg) for 12 days. A, Tumor volume was measured on treatment day 12. Each data point is the tumor volume in a single mouse. Midlines, average tumor volume; error bars, SD. B, Tumor sections were stained with hematoxylin and eosin (H&E) or by IHC to detect Ki67. C, The percentage of Ki67+ tumor cells was quantitated in 5 random fields per sample. Midline, average of biological replicates. D, 4T1 and MDA-MB-231 cells were treated with CPI211 and/or U46619 in serum-free media for 3 days. 10% serum was included as a positive control for growth. Value points represent the average of three technical replicates ± SD; N = 3. E, IHC to detect CD31+ vessels in 4T1 tumors. Inset represents the area shown in higher power. F and G, Quantitation of total CD31+ vessel area (F) and number of CD31+ vessels with a lumen area > 50 μm2 (G). Data points, average value measured in three random fields per section; N, 4. Midlines, average ± SD. n.s., not significant.

Close modal

CPI211 decreased tumor vasculature

An abundance of large vascular structures seen in vehicle-treated tumors was substantially diminished in CPI211-treated tumors, as shown by staining with the vessel endothelial marker CD31 (Fig. 3E). Software-based morphometric analyses of CD31-stained tumor sections revealed a decreased total area of CD31+ vessels in CPI211-treated tumors (Fig. 3F) and fewer CD31+ structures exceeding 50 μm2 (Fig. 3G). Thus, it is possible that TPr signaling affects the tumor vasculature in a manner that promotes metastasis.

TPr signaling supports increased transendothelial migration of tumor cells

We investigated the impact of TPr signaling on ECs using primary cultures of mouse pulmonary microvascular ECs (MPMEC), which abundantly expressed TPr (Fig. 4A). U46619 induced PKC activity in MPMECs. As expected, TPr expression was not diminished in U46619-treated MPMECs, because mouse cells express only TPr-α, which, unlike TPr-β, is not internalized in response to ligand activation (19). Importantly, PKC substrate phosphorylation was potently inhibited by CPI211. Similar PKC substrate phosphorylation was seen in Western blot analyses of primary HUVECs treated with U46619, which was blocked upon treatment with CPI211. Notably, TPr expression diminished in U46619-treated HUVECs, which abundantly express TPr-β (30). Agonist-induced TPr diminution was blocked by CPI211 in HUVECs.

Figure 4.

The TPr inhibitor CPI211 increases barrier function of vascular ECs. A, Western analysis of MPMECs and HUVECs treated with CPI211 (100 nmol/L) ± U46619 (500 nmol/L). B, Crystal violet staining and quantitation measuring growth of MPMECs treated with U46619 (500 nmol/L) and CPI211 (100 nmol/L). C, Transendothelial tumor cell migration assay. ECs (CellTracker Red+) monolayers on the underside of Matrigel-coated transwells were treated for 6 hours with U46619 (500 nmol/L) and CPI211 (100 nmol/L) before adding GFP+ tumor cells to upper chambers. GFP+ cell migration to the lower transwell surface was visualized by fluorescence microscopy. D, Representative images of transendothelial migration of GFP+ 4T1 and MDA-MB-231 across mouse MPMECs and human HUVECs, respectively. E, Quantitation of GFP+ 4T1 (N = 8) and MDA-MB-231 (N = 4) cells migrating across an endothelial barrier treated with U46619 ± CPI211. Data points are average ± SD. F, Albumin leakage across an MPMEC monolayer grown on Matrigel-coated transwell; N = 3, each assessed in triplicate. Average ± SD is shown. Student t test. G, IF staining for VE-Cadherin in MPMECs grown on a Matrigel-coated transwell filter and treated for 6 hours with U46619. H, Western analysis of whole-cell lysates (left) or VE-Cadherin IP (right) from MPMECs treated for 6 hours with U46619 ± CPI211.

Figure 4.

The TPr inhibitor CPI211 increases barrier function of vascular ECs. A, Western analysis of MPMECs and HUVECs treated with CPI211 (100 nmol/L) ± U46619 (500 nmol/L). B, Crystal violet staining and quantitation measuring growth of MPMECs treated with U46619 (500 nmol/L) and CPI211 (100 nmol/L). C, Transendothelial tumor cell migration assay. ECs (CellTracker Red+) monolayers on the underside of Matrigel-coated transwells were treated for 6 hours with U46619 (500 nmol/L) and CPI211 (100 nmol/L) before adding GFP+ tumor cells to upper chambers. GFP+ cell migration to the lower transwell surface was visualized by fluorescence microscopy. D, Representative images of transendothelial migration of GFP+ 4T1 and MDA-MB-231 across mouse MPMECs and human HUVECs, respectively. E, Quantitation of GFP+ 4T1 (N = 8) and MDA-MB-231 (N = 4) cells migrating across an endothelial barrier treated with U46619 ± CPI211. Data points are average ± SD. F, Albumin leakage across an MPMEC monolayer grown on Matrigel-coated transwell; N = 3, each assessed in triplicate. Average ± SD is shown. Student t test. G, IF staining for VE-Cadherin in MPMECs grown on a Matrigel-coated transwell filter and treated for 6 hours with U46619. H, Western analysis of whole-cell lysates (left) or VE-Cadherin IP (right) from MPMECs treated for 6 hours with U46619 ± CPI211.

Close modal

Growth of MPMECs seeded at low density over 7 days was unaffected by U46619 and/or CPI211 (Fig. 4B). Interestingly, serum-free 4T1 tumor cell–conditioned media increased MPMEC growth, although this, too, was unaffected by CPI211.

Accumulating evidence shows that EC barrier function becomes compromised in cancers (31–33). We measured tumor cell migration across an endothelial barrier using a modified transwell approach (Fig. 4C). CellTracker Red-labeled MPMECs seeded on Matrigel on the lower side of a transwell filter were grown into a barrier monolayer, and then treated with U46619 prior to adding 4T1-GFP cells into the upper chamber, assessing transendothelial migrations by fluorescence microscopy. These studies revealed increased 4T1 transendothelial migration across U46619-treated MPMECs as compared with vehicle-treated MPMECs (Fig. 4D and E). Similarly, an increased number of MDA-MB-231 cells migrated across U46619-treated HUVECs as compared with vehicle-treated HUVECs. In both cases, CPI211 blocked tumor cell migration across the EC barrier. However, in the absence of ECs, 4T1 cell migration across Matrigel-coated transwells was unaffected by U46619 or by CPI211 (Supplementary Fig. S2A), as was MDA-MB-231 migration (Supplementary Fig. S2B). These findings suggest that TPr signaling may not affect tumor cell migration, per se, but may enhance transendothelial migration of tumor cells.

Next, we measured the barrier function of MPMEC monolayers cultured on transwells and treated with U46619 and/or CPI211. Transwell filters were transferred to fresh wells containing only PBS after 6 hours, at which point albumin (50 μg/mL) was added to the upper transwell chamber. After 2 minutes, the lower chamber was sampled for albumin distribution (Fig. 4F). U46619-treated MPMECs showed increased albumin passage to the lower chamber, which was blocked by CPI211. These results suggest that TPr signaling may increase permeability of an endothelial barrier, and are consistent with previous reports suggesting that thromboxane signaling may decrease vascular barrier function (34, 35).

VE-cadherin localization at EC junctions is critical for endothelial barrier function. Notably, VE-cadherin tyrosine phosphorylation induces VE-cadherin internalization, compromising endothelial barrier function (36), enabling metastatic dissemination of tumor cells (37, 38). Using immunofluorescence (IF) to assess VE-cadherin localization in MPMECs cultured on transwell filters, we identified punctate patterns of VE-cadherin staining along cell membranes at or near cell–cell contacts in vehicle-treated MPMECs (Fig. 4G), whereas U46619-treated MPMECs displayed diffuse intracellular VE-cadherin staining. Although U46619 did not affect total VE-cadherin expression (Fig. 4H), U46619-induced VE-cadherin tyrosine phosphorylation, which was blocked by CPI211 (Fig. 4H), suggesting that TPr signaling may induce VE-cadherin tyrosine phosphorylation and internalization, leading to decreased barrier function.

Neoadjuvant and adjuvant models of CPI211 treatment diminish metastatic spread from primary tumors in immune-competent mice

Next, we established models of spontaneous metastasis from orthotopic 4T1 tumors grown in immune-competent isogenic Balb/C mice, with two separate CPI211 treatment schemes assessed. The first approach modeled neoadjuvant CPI211 treatment (i.e., without tumor resection), in which tumor-bearing mice were randomized into treatment groups when tumors reached 50 mm3. Mice were treated daily for 28 days with vehicle or CPI211 (Fig. 5A). Although lung metastases were seen in all vehicle-treated mice (N = 10), only 60% of CPI211-treated mice harbored lung metastases (Fig. 5B). Of the mice that developed lung metastases, fewer lung metastases were seen CPI211-treated samples versus vehicle, demonstrating that TPr inhibition diminishes spontaneous metastases in the neoadjuvant setting.

Figure 5.

CPI211 diminishes spontaneous metastases. A and B, 4T1 tumor–bearing wild-type (Balb/C) mice were treated 28 days with neoadjuvant CPI211 (50 mg/kg daily) or vehicle when tumors reached 50 mm3 (A). Representative low-power images (20×) of H&E-stained sections of lungs collected at day 28 (B, left). Lung metastases were enumerated (B, right). Data points are the number of metastases in each mouse. Midlines are average ± SD; N = 10. C and D, 4T1 tumor–bearing wild-type (Balb/C) mice were treated with adjuvant CPI211 (50 mg/kg daily) or vehicle beginning after surgical resection of 300 mm3 tumors (C). Representative low-power images (20×) of H&E-stained sections of lungs collected at day 28 are shown (D, left). The inset shows the area of increased magnification. Lung metastases were enumerated (D, right). Midlines are the average ± SD; N = 8.

Figure 5.

CPI211 diminishes spontaneous metastases. A and B, 4T1 tumor–bearing wild-type (Balb/C) mice were treated 28 days with neoadjuvant CPI211 (50 mg/kg daily) or vehicle when tumors reached 50 mm3 (A). Representative low-power images (20×) of H&E-stained sections of lungs collected at day 28 (B, left). Lung metastases were enumerated (B, right). Data points are the number of metastases in each mouse. Midlines are average ± SD; N = 10. C and D, 4T1 tumor–bearing wild-type (Balb/C) mice were treated with adjuvant CPI211 (50 mg/kg daily) or vehicle beginning after surgical resection of 300 mm3 tumors (C). Representative low-power images (20×) of H&E-stained sections of lungs collected at day 28 are shown (D, left). The inset shows the area of increased magnification. Lung metastases were enumerated (D, right). Midlines are the average ± SD; N = 8.

Close modal

We next tested a model of adjuvant TPr inhibition following surgical resection of 300 mm3 tumors. Mice were randomized for treatment 24 hours after tumor resection (Fig. 5C), receiving daily treatment with vehicle or CPI211 for 21 days. Although 100% of vehicle-treated mice and 90% of CPI211-treated mice developed lung metastases (Fig. 5D), the number of lung metastases per mouse was reduced in CPI211-treated mice to half of what was seen in controls. These studies reveal a potential benefit for reducing metastatic burden through TPr inhibition, even after tumor cells may have entered circulation.

We identified a novel and unexpected correlation between TPr and tumor metastasis using a combination of PheWAS analyses in large population data sets, gene-expression studies in clinical cancer data sets, and pharmacologic blockade of TPr activity in models of metastasis (Supplementary Fig. S3). Importantly, the TPr inhibitor CPI211 decreased both experimental (hematogenous) and spontaneous tumor metastasis.

The metastatic process requires tumor cells to overcome a series of physiologic barriers (2). Tumor cells must invade through primary tumor extracellular matrix to reach the vasculature, migrate across the vascular endothelium into circulation, survive in circulation, evade immune recognition, and adhere to and migrate across the endothelium into distant tissue, where tumor cells must survive and grow (39). Findings shown herein highlight an underappreciated role for TPr signaling in disrupting the endothelial barrier in the metastatic process, facilitating tumor metastasis. These studies build upon published studies demonstrating that platelets, a rich source of the TPr ligand TxA2, are potent mediators of metastasis through a range of pathways involving interactions between tumor cells, ECs, and platelets (40–42). Platelet-CTC aggregation may enhance metastasis by providing important survival signals to the tumor cell, by shielding CTCs from immune surveillance, and by promoting adhesion to the vascular endothelium. In fact, previous studies report that TBXA2R-null mice display reduced metastases in experimental (tail-vein) models of metastasis (42). Thus, CPI211 may reduce platelet-assisted metastatic processes. Although we have not tested this hypothesis specifically, we confirmed that CPI211 inhibits platelet aggregation (Fig. 1).

Although platelet TPr likely has a prominent role in facilitating tumor metastasis, we show that endothelial TPr signaling increases transendothelial tumor cell migration, even in the absence of platelets (Fig. 4). Endothelial TPr signaling diminished VE-cadherin at endothelial cell–cell junctions. Interestingly, published studies show that hypoxia contributes to VE-cadherin endocytosis in ECs, resulting in endothelial barrier disruption and vascular dysfunction, while increasing TPr expression and sensitivity (38, 43), together suggesting that tumor hypoxia might increase endothelial TPr signaling, thus facilitating metastasis, although this hypothesis remains to be tested.

Recent studies highlighted a role for COX-1 in generation of a premetastatic niche within the vasculature of colorectal cancers, in large part through its ability to increase TxA2 production, thus enhancing platelet–tumor cell aggregation upon vessel walls. Inhibition of COX-1 using aspirin decreased metastasis through downregulation of platelet TxA2 (41), implicating TPr signaling in metastasis. These studies hold great clinical significance, as the Add-Aspirin study (NCT02804815) is currently testing if daily aspirin use (100 mg, 300 mg, or placebo, for 5 years) after SOC treatment prevents recurrence and improves survival in patients lacking clinically evident metastases at the conclusion of SOC treatment.

It is important to note that COX-1 produces both TxA2 and gut-protecting prostaglandin E2 (PGE2). As such, aspirin may have dangerous consequences for some patients. It is exciting to speculate TPr inhibition could provide the antimetastatic benefits of COX1 pathway inhibition, but without the potentially harmful side effects in aspirin-sensitive patients. However, given the potent antagonism of platelet coagulation by CPI211, studies going forward should proceed with awareness of bleeding or other potential adverse events that may arise in CPI211-treated patients.

Because TPr inhibition diminished metastasis without impacting tumor growth, it is anticipated that TPr inhibition could be effective in combination with SOC cytotoxic therapies in the adjuvant or neoadjuvant setting. Data presented herein suggest that TPr inhibitors could provide greatest benefit as a preventive treatment, a clinical course of action not frequently used. However, this approach is being tested in at least one clinical trial (NCT03694249), in which patients diagnosed with aggressive cancers with a high risk of recurrence (e.g., gastric cancer; triple-negative breast cancer and others) but lacking clinically identifiable metastases will be randomized post-surgically into treatment groups receiving CPI211 or placebo for 12 months, and followed for an additional 12 months. This landmark clinical trial will test the feasibility of a pharmaco-preventive strategy that may diminish metastatic burden in certain cancer patients.

The notion that TPr inhibitors will not eliminate existing metastases must be considered. However, TPr inhibition may reduce spread of tumor cells from secondary to tertiary tumor sites. Timing TPr inhibitor administration requires additional study, as does testing of other TPr inhibitors in clinical development (terutroban, ridogrel, terbogrel, and others; ref. 8). Further, it will be essential to identify which cancer types are responsive to TPr inhibition, necessitating a greater understanding of how different cancers use different metastatic routes.

In summary, we introduce TPr as an antimetastatic target across many types of cancers. Findings herein support the continued investigation of agents interfering with steps of the metastatic cascade, and support the use of PheWAS analysis to accelerate the pace of drug repurposing.

CPI211 used in this study was provided by Cumberland Pharmaceuticals, Inc. (CPI). CPI had no role in study design, experimental conduct, data collection and analysis, or preparation of the manuscript. Additionally, CPI provided no direct financial support for this work.

T.A. Werfel: Investigation, writing–original draft. D.J. Hicks: Investigation. B. Rahman: Investigation. W.E. Bindeman: Investigation. M.T. Duvernay: Investigation. J.G. Maeng: Data curation, investigation. H. Hamm: Resources, supervision. R.R. Lavieri: Conceptualization, data curation, investigation, writing–review and editing. M.M. Joly: Conceptualization, investigation, writing–review and editing. J.M. Pulley: Conceptualization, resources, funding acquisition. D.L. Elion: Investigation, writing–review and editing. D.M. Brantley-Sieders: Validation, investigation, writing–review and editing. R.S. Cook: Conceptualization, resources, formal analysis, supervision, funding acquisition, investigation, writing–original draft, project administration, writing–review and editing.

We kindly thank Cumberland Pharmaceuticals, Inc., and Dr. Ines Macias-Perez for providing CPI211. We acknowledge Drs. Ingrid Meier and Justin Balko for thoughtful discussion. We acknowledge Vanderbilt Shared Resources: the VICC Breast SPORE Pathology (Dr. Melinda Sanders), Translational Pathology, Digital Histology, and VANTAGE shared resources. This work was supported by grants NIH P50 CA098131, NIH P30 CA68485, National Center for Advancing Translational Sciences CTSA UL1TR000445, and Congressionally Directed Medical Research Program W81XWH-161-0063.

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Denny
JC
,
Ritchie
MD
,
Basford
MA
,
Pulley
JM
,
Bastarache
L
,
Brown-Gentry
K
, et al
PheWAS: demonstrating the feasibility of a phenome-wide scan to discover gene-disease associations
.
Bioinformatics
2010
;
26
:
1205
10
.
2.
Pulley
JM
,
Jerome
RN
,
Ogletree
ML
,
Bernard
GR
,
Lavieri
RR
,
Zaleski
NM
, et al
Motivation for launching a cancer metastasis inhibition (CMI) program
.
Target Oncol
2018
;
13
:
61
8
.
3.
Pulley
JM
,
Shirey-Rice
JK
,
Lavieri
RR
,
Jerome
RN
,
Zaleski
NM
,
Aronoff
DM
, et al
Accelerating precision drug development and drug repurposing by leveraging human genetics
.
Assay Drug Dev Technol
2017
;
15
:
113
9
.
4.
Schumacher
WA
,
Steinbacher
TE
,
Youssef
S
,
Ogletree
ML
. 
Antiplatelet activity of the long-acting thromboxane receptor antagonist BMS 180,291 in monkeys
.
Prostaglandins
1992
;
44
:
389
97
.
5.
Jones
RL
,
Giembycz
MA
,
Woodward
DF
. 
Prostanoid receptor antagonists: development strategies and therapeutic applications
.
Br J Pharmacol
2009
;
158
:
104
45
.
6.
Ogletree
ML
,
Harris
DN
,
Schumacher
WA
,
Webb
ML
,
Misra
RN
. 
Pharmacological profile of BMS 180,291: a potent, long-acting, orally active thromboxane A2/prostaglandin endoperoxide receptor antagonist
.
J Pharmacol Exp Ther
1993
;
264
:
570
8
.
7.
95th annual meeting of the American Society for Clinical Pharmacology and Therapeutics. New Orleans, Louisiana, March 30-April 1, 1994. Abstracts
.
Clin Pharmacol Ther
1994
;
55
:
123
218
.
8.
Fontana
P
,
Zufferey
A
,
Daali
Y
,
Reny
JL
. 
Antiplatelet therapy: targeting the TxA2 pathway
.
J Cardiovasc Transl Res
2014
;
7
:
29
38
.
9.
Nakahata
N
. 
Thromboxane A2: physiology/pathophysiology, cellular signal transduction and pharmacology
.
Pharmacol Ther
2008
;
118
:
18
35
.
10.
Jerome
RN
,
Pulley
JM
,
Roden
DM
,
Shirey-Rice
JK
,
Bastarache
LA
,
G
RB
, et al
Using human ‘experiments of nature’ to predict drug safety issues: an example with PCSK9 inhibitors
.
Drug Saf
2018
;
41
:
303
11
.
11.
Mundell
SJ
,
Mumford
A
. 
TBXA2R gene variants associated with bleeding
.
Platelets
2018
;
29
:
739
42
.
12.
Yi
X
,
Lin
J
,
Zhou
Q
,
Huang
R
,
Chai
Z
. 
The TXA2R rs1131882, P2Y1 rs1371097 and GPIIIa rs2317676 three-loci interactions may increase the risk of carotid stenosis in patients with ischemic stroke
.
BMC Neurol
2019
;
19
:
44
.
13.
Milanowski
L
,
Pordzik
J
,
Janicki
PK
,
Kaplon-Cieslicka
A
,
Rosiak
M
,
Peller
M
, et al
New single-nucleotide polymorphisms associated with differences in platelet reactivity and their influence on survival in patients with type 2 diabetes treated with acetylsalicylic acid: an observational study
.
Acta Diabetol
2017
;
54
:
343
51
.
14.
Bauer
J
,
Ripperger
A
,
Frantz
S
,
Ergun
S
,
Schwedhelm
E
,
Benndorf
RA
. 
Pathophysiology of isoprostanes in the cardiovascular system: implications of isoprostane-mediated thromboxane A2 receptor activation
.
Br J Pharmacol
2014
;
171
:
3115
31
.
15.
Takeuchi
K
,
Mashimo
Y
,
Shimojo
N
,
Arima
T
,
Inoue
Y
,
Morita
Y
, et al
Functional variants in the thromboxane A2 receptor gene are associated with lung function in childhood-onset asthma
.
Clin Exp Allergy
2013
;
43
:
413
24
.
16.
Sasaki
M
,
Sukegawa
J
,
Miyosawa
K
,
Yanagisawa
T
,
Ohkubo
S
,
Nakahata
N
. 
Low expression of cell-surface thromboxane A2 receptor beta-isoform through the negative regulation of its membrane traffic by proteasomes
.
Prostaglandins Other Lipid Mediat
2007
;
83
:
237
49
.
17.
Hamelin
E
,
Theriault
C
,
Laroche
G
,
Parent
JL
. 
The intracellular trafficking of the G protein-coupled receptor TPbeta depends on a direct interaction with Rab11
.
J Biol Chem
2005
;
280
:
36195
205
.
18.
Laroche
G
,
Rochdi
MD
,
Laporte
SA
,
Parent
JL
. 
Involvement of actin in agonist-induced endocytosis of the G protein-coupled receptor for thromboxane A2: overcoming of actin disruption by arrestin-3 but not arrestin-2
.
J Biol Chem
2005
;
280
:
23215
24
.
19.
Parent
JL
,
Labrecque
P
,
Orsini
MJ
,
Benovic
JL
. 
Internalization of the TXA2 receptor alpha and beta isoforms. Role of the differentially spliced cooh terminus in agonist-promoted receptor internalization
.
J Biol Chem
1999
;
274
:
8941
8
.
20.
Roden
DM
,
Pulley
JM
,
Basford
MA
,
Bernard
GR
,
Clayton
EW
,
Balser
JR
, et al
Development of a large-scale de-identified DNA biobank to enable personalized medicine
.
Clin Pharmacol Ther
2008
;
84
:
362
9
.
21.
Pulley
J
,
Clayton
E
,
Bernard
GR
,
Roden
DM
,
Masys
DR
. 
Principles of human subjects protections applied in an opt-out, de-identified biobank
.
Clin Transl Sci
2010
;
3
:
42
8
.
22.
Carroll
RJ
,
Bastarache
L
,
Denny
JC
. 
R PheWAS: data analysis and plotting tools for phenome-wide association studies in the R environment
.
Bioinformatics
2014
;
30
:
2375
6
.
23.
Brantley-Sieders
DM
,
Fang
WB
,
Hwang
Y
,
Hicks
D
,
Chen
J
. 
Ephrin-A1 facilitates mammary tumor metastasis through an angiogenesis-dependent mechanism mediated by EphA receptor and vascular endothelial growth factor in mice
.
Cancer Res
2006
;
66
:
10315
24
.
24.
Duvernay
M
,
Young
S
,
Gailani
D
,
Schoenecker
J
,
Hamm
HE
. 
Protease-activated receptor (PAR) 1 and PAR4 differentially regulate factor V expression from human platelets
.
Mol Pharmacol
2013
;
83
:
781
92
.
25.
Elion
DL
,
Jacobson
ME
,
Hicks
DJ
,
Rahman
B
,
Sanchez
V
,
Gonzales-Ericsson
PI
, et al
Therapeutically active RIG-I agonist induces immunogenic tumor cell killing in breast cancers
.
Cancer Res
2018
;
78
:
6183
95
.
26.
Qu
S
,
Rinehart
C
,
Wu
HH
,
Wang
SE
,
Carter
B
,
Xin
H
, et al
Gene targeting of ErbB3 using a Cre-mediated unidirectional DNA inversion strategy
.
Genesis
2006
;
44
:
477
86
.
27.
Youngblood
V
,
Wang
S
,
Song
W
,
Walter
D
,
Hwang
Y
,
Chen
J
, et al
Elevated Slit2 activity impairs VEGF-induced angiogenesis and tumor neovascularization in EphA2-deficient endothelium
.
Mol Cancer Res
2015
;
13
:
524
37
.
28.
Curtis
C
,
Shah
SP
,
Chin
SF
,
Turashvili
G
,
Rueda
OM
,
Dunning
MJ
, et al
The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups
.
Nature
2012
;
486
:
346
52
.
29.
Gyorffy
B
,
Surowiak
P
,
Budczies
J
,
Lanczky
A
. 
Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in non-small-cell lung cancer
.
PLoS One
2013
;
8
:
e82241
.
30.
Ashton
AW
,
Ware
JA
. 
Thromboxane A2 receptor signaling inhibits vascular endothelial growth factor-induced endothelial cell differentiation and migration
.
Circ Res
2004
;
95
:
372
9
.
31.
De Bock
K
,
Cauwenberghs
S
,
Carmeliet
P
. 
Vessel abnormalization: another hallmark of cancer? Molecular mechanisms and therapeutic implications
.
Curr Opin Genet Dev
2011
;
21
:
73
9
.
32.
Wieland
E
,
Rodriguez-Vita
J
,
Liebler
SS
,
Mogler
C
,
Moll
I
,
Herberich
SE
, et al
Endothelial Notch1 activity facilitates metastasis
.
Cancer Cell
2017
;
31
:
355
67
.
33.
Azzi
S
,
Hebda
JK
,
Gavard
J
. 
Vascular permeability and drug delivery in cancers
.
Front Oncol
2013
;
3
:
211
.
34.
Zhao
Z
,
Hu
J
,
Gao
X
,
Liang
H
,
Yu
H
,
Liu
S
, et al
Hyperglycemia via activation of thromboxane A2 receptor impairs the integrity and function of blood-brain barrier in microvascular endothelial cells
.
Oncotarget
2017
;
8
:
30030
8
.
35.
Xie
X
,
Sun
W
,
Wang
J
,
Li
X
,
Liu
X
,
Liu
N
. 
Activation of thromboxane A2 receptors mediates endothelial dysfunction in diabetic mice
.
Clin Exp Hypertens
2017
;
39
:
312
8
.
36.
Dejana
E
,
Orsenigo
F
,
Lampugnani
MG
. 
The role of adherens junctions and VE-cadherin in the control of vascular permeability
.
J Cell Sci
2008
;
121
:
2115
22
.
37.
Cai
J
,
Jiang
WG
,
Mansel
RE
. 
Phosphorylation and disorganization of vascular-endothelial cadherin in interaction between breast cancer and vascular endothelial cells
.
Int J Mol Med
1999
;
4
:
191
5
.
38.
Mazzone
M
,
Dettori
D
,
de Oliveira
RL
,
Loges
S
,
Schmidt
T
,
Jonckx
B
, et al
Heterozygous deficiency of PHD2 restores tumor oxygenation and inhibits metastasis via endothelial normalization
.
Cell
2009
;
136
:
839
51
.
39.
Talmadge
JE
,
Fidler
IJ
. 
AACR centennial series: the biology of cancer metastasis: historical perspective
.
Cancer Res
2010
;
70
:
5649
69
.
40.
Quail
DF
,
Joyce
JA
. 
Microenvironmental regulation of tumor progression and metastasis
.
Nat Med
2013
;
19
:
1423
37
.
41.
Lucotti
S
,
Cerutti
C
,
Soyer
M
,
Gil-Bernabe
AM
,
Gomes
AL
,
Allen
PD
, et al
Aspirin blocks formation of metastatic intravascular niches by inhibiting platelet-derived COX-1/thromboxane A2
.
J Clin Invest
2019
;
129
:
1845
62
.
42.
Matsui
Y
,
Amano
H
,
Ito
Y
,
Eshima
K
,
Suzuki
T
,
Ogawa
F
, et al
Thromboxane A(2) receptor signaling facilitates tumor colonization through P-selectin-mediated interaction of tumor cells with platelets and endothelial cells
.
Cancer Sci
2012
;
103
:
700
7
.
43.
Sawada
J
,
Urakami
T
,
Li
F
,
Urakami
A
,
Zhu
W
,
Fukuda
M
, et al
Small GTPase R-Ras regulates integrity and functionality of tumor blood vessels
.
Cancer Cell
2012
;
22
:
235
49
.

Supplementary data