Abstract
Background: Infiltrating lobular carcinoma (ILC) is the second most common histologicaI subtype of breast cancer, accounting for 10% of all cases. ILC has a characteristic genomic profile. ILC shows a high frequency of cadherin 1 (CDH1) mutations, along with loss of phosphatase and tensin homolog (PTEN), activation of alpha serine/threonine kinase (AKT), and mutations in T-box transcription factor (TBX3) and forkhead box protein A1 (FOXA1). We suspected that another gene, von Willebrand factor (VWF), might also be part of the profile, since coagulation tests reveal significant differences in patients with breast cancer. Materials and Methods: For newly-diagnosed breast cancer, the association between VWF and histology in the GDC Breast Cancer dataset in The Cancer Genome Atlas (TCGA) was evaluated. The following were used to access and analyze the data: Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/); Xena browser (https://xenabrowser.net); cBioportal (http://cbioportal.org); Oncomine (https://oncomine.org); and Prediction Analysis of Microarray 50 (PAM50). Results: Patients with ILC had higher VWF RNA expression than patients with infiltrating ductal carcinoma and other histology. The difference of expression was present to the same degree in both pre-menopausal and post-menopausal patients. Nine alterations in VWF and PTEN were significantly co-occurrent. Considering all histologies in 843 samples, Tukey's honest significant difference post hoc test showed that VWF RNA expression of the normal subtype was significantly greater than that of the other subtypes (p<0.001). Conclusion: Our finding of significantly higher VWF RNA expression in the PAM50 normal (non-basal-like) breast cancer subtype suggests that VWF protein measurement might complement or corroborate PAM50 results. VWF and PAM50 results both suggesting a low risk of recurrence might make the decision whether to give chemotherapy easier, especially if VWF protein were an independent predictor.
Infiltrating lobular carcinoma (ILC) is the second most common histological subtype of breast cancer, accounting for 10% of all cases. ILC differs from infiltrating ductal carcinoma (IDC) in its pathological features and responsiveness to systemic therapy. ILC derives a distinct benefit from systemic therapy compared to IDC (1).
ILC has a characteristic genomic profile, exhibiting a high frequency of cadherin 1 (CDH1) mutation, along with loss of phosphatase and tensin homolog (PTEN), activation of alpha serine/threonine kinase (AKT), and mutations in T-box transcription factor (TBX3) and forkhead box protein A1 (FOXA1) (2, 3). The genomic profile of lobular carcinoma in situ is similar (4).
We suspected that another gene, von Willebrand factor (VWF), might be part of the profile, since coagulation tests reveal significant differences in patients with breast cancer (5, 6). One report documents VWF up-regulation in ILC and suggests an effect on cell adhesion (7).
VWF, the largest human plasma protein, is a multimeric glycoprotein that mediates platelet adhesion to both the subendothelial matrix and endothelial surfaces and acts as a carrier for coagulation factor VIII in the circulation (6). Besides its essential role in hemostasis, VWF has an effect on tumors, mainly by inhibiting angiogenesis and apoptosis (6). Increased expression of VWF gene is associated with poorer survival in primary lower-grade gliomas (8).
Comparison of alterations in genes in the entire The Cancer Genome Atlas breast cohort. The von Willebrand factor (VWF) and phosphatase and tensin homolog (PTEN) gene pair had nine significantly co-occurrent alterations, while PTEN and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) had 25 mutually exclusive alterations.
In the current study, we used The Cancer Genome Atlas (TCGA) to investigate whether VWF might be part of the genetic profile of ILC.
Materials and Methods
For newly diagnosed breast cancer, this study evaluated the association between VWF gene expression and histology in the GDC Breast Cancer dataset in TCGA (9). The following were used to access and analyze the data: Genomic Data Commons Data Portal (https://portal.gdc.cancer.gov/); Xena browser (https://xenabrowser.net) (10); cBioportal (http://cbioportal.org) (11); Oncomine (https://oncomine.org) (12); and Prediction Analysis of Microarray 50 (PAM50), a Food and Drug Administration-approved tooI for intrinsic subtyping of breast cancer using quantitative reverse transcriptase-polymerase chain reaction that assays the expression of 50 genes. VWF Data from these assays were sourced from Oncomine. Comparisons were made only for VWF mRNA. In the case of the PAM50, we did not look at data from any of the 50 genes included in the test; we looked only at whether the PAM50 result was reported as normal or abnormal.
Gene expression is reported as fragments per kilobase of transcript per million mapped reads upper quartile (fpkm-uq), which is an RNA-Seq-based method for expression normalization (13).
Simple statistics were calculated to identify patterns of mutual exclusivity or co-occurrence. For each pair of query genes (e.g. VWF and PTEN), an odds ratio (OR) was calculated (Equation 1) that indicates the likelihood that the events in the two genes are mutually exclusive or co-occurrent across the selected cases
(Eq. 1)
where A was the number of cases in which both genes were altered; B was the number of cases in which VWF but not PTEN was altered; C was the number of cases in which PTEN but not VWF was altered; and D was the number of cases in which neither gene was altered. Each pair was then assigned to one of three categories indicative of a tendency toward mutual exclusivity, co-occurrence, or of no association, respectively. Only Fisher's exact test was performed to determine whether the identified relationship was significant for a gene pair. The test is used to support co-occurrence when the number of tumors with alterations in both genes is significantly higher than expected by chance. Likewise, it suggests mutual exclusivity when the number of tumors with alterations in both genes is significantly lower (11).
Distribution of age at diagnosis of women patients in the GDC The Cancer Genome Atlas breast cancer cohort.
Results
We analyzed the Breast Invasive Carcinoma (TCGA) cohort of 1,113 samples from women. The age distribution of patients is shown in Figure 1. Race information: 69.3% were White, 16.8% were Black, 5.5% Asian, 0.1% American Indian, 8.3% unclassified. Ethnicity information: 80% were non-Hispanic, 4% were Hispanic, 16% were unclassified (unclassified can refer to race or non-Hispanic, Hispanic classifications). Out of the tumors that were staged, 26% were stage T1, 58% stage T2, 12% stage T3, and 4% stage T4 (one tumor had no staging).
VWF mRNA expression was significantly related to histology. Patients with ILC had higher VWF mRNA expression than patients with IDC and other histology (Figure 2). The difference of VWF mRNA expression was present to the same degree in both pre-menopausal and post-menopausal breast cancer. Data from Zhao et al. (7) show a 2.331-fold change (ILC/IDC) in VWF mRNA expression that is significant (data shown in Figure 3; these data confirm the analysis in Figure 2).
Tukey's honest significant difference post hoc test confirmed that VWF mRNA expression in ILC was significantly increased over that of IDC (p<0.001), medullary carcinoma (p=0.038), and mucinous carcinoma (p=0.004). Excluding ILC, VWF mRNA expression in IDC was not significantly different from that of any other histology.
von Willebrand factor (VWF) mRNA expression according to histology. Patients with infiltrating lobular carcinoma (ILC) had higher VWF mRNA expression than patients with infiltrating ductal carcinoma (IDC) and other histologies. Tukey's honest significant difference post hoc test: ILC vs.: IDC, p<0.001; medullary carcinoma (Med.), p=0.038; and mucinous carcinoma (Muc.), p=0.004. Excluding ILC, VWF mRNA expression in IDC did not significantly differ from that of any other histology. Met: Metaplastic carcinoma; NOS: infiltrating carcinoma, not otherwise specified.
Because of the relationship between ABO blood group and serum VWF levels (14), we performed univariate analysis of variance, general linear model using VWF mRNA expression as the dependent variable, histological type as the fixed factor, and ABO mRNA expression as covariate. Variability of VWF mRNA expression was significant (p<0.001) and independent of ABO mRNA expression (p<0.001).
PTEN mutations are characteristic of ILC (2, 3). Nine alterations in VWF and PTEN were significantly co-occurrent (p=0.007; Table I and Figure 4). In contrast, 25 alterations in PTEN and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) were mutually exclusive.
VWF is part of a network of neighboring genes, including PTEN, PIK3CA, and erythroblastic leukemia viral oncogene homolog 2 (ERB2), implicated in the genesis of ILC (Figure 5).
Figure 6 shows the correlation of VWF and PTEN mRNA expression in 963 breast cancer cases. The positive correlation was significant (p=0.000627); although weak, the correlation is still important.
Data from Zhao et al. (24) show a 2.331-fold increase of von Willebrand factor (VWF) in infiltrating lobular carcinoma (ILC) compared to infiltrating ductal carcinoma (IDC). To show the fold change, microarray data was log2-transformed to generate a normal distribution of the data, a requirement for most statistical tests and a commonly used scale for microarray data. Oncomine displays log2, median centered values on the y-axis of bar charts, as above (https://oncomine.org).
Diagrammatic representation of nine significantly co-occurring alterations in von Willebrand factor (VWF) and phosphatase and tensin homolog (PTEN). VWF and PTEN were altered in 3% and 11%, respectively, of 996 sequenced cases (cBioportal.org).
von Willebrand factor (VWF) and its network of neighboring genes, including phosphatase and tensin homolog (PTEN), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) and erythroblastic leukemia viral oncogene homolog 2 (ERB2), implicated in the genesis of infiltrating lobular carcinoma. An arrow indicates a directed interaction and a line an undirected interaction. Green arrow: Control of phosphorylation; blue arrow: control of state of change; brown lines indicate targeting by a drug (cBioportal.org). FGF3,4,19: Fibroblast growth factor 3,4,19; YWHAZ: tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta; FGFR1: fibroblast growth factor receptor 1; PTK2: protein tyrosine kinase 2; TP53: tumor protein 53; COL14A1: collagen type XIV alpha 1 chain; ACN2: actinin alpha 2.
Mixed ILC/IDC breast cancers have distinct histopathological characteristics and are often diagnosed at a late stage, although survival does not vary significantly from ILC and IDC (15). In the TCGA cohort, 30 cases were recorded as having a mixed histology; however, TCGA does not specify what the mixture was.
Correlation of mRNA expression of von Willebrand factor (VWF) with that of phosphatase and tensin homolog (PTEN) in 963 breast cancer cases (p=0.000627).
VWF mRNA expression in ILC was not related to T-stage (p=0.2) nor the presence of metastases (p=0.97). Nor was VWF mRNA expression in IDC affected by T-stage or the presence of metastases (p=0.55). VWF RNA expression, copy number, or mutations had no effect on survival of patients with IDC or ILC.
VWF mRNA expression was related to PAM50 subtype (Figure 7). Considering all histologies in 843 samples, Tukey's honest significant difference post hoc test showed that VWF mRNA expression of the normal subtype was significantly greater than that of the other subtypes (p<0.001).
A receiver operating characteristic curve was constructed with the data used in Figure 7 to determine the optimum cut-off for VWF mRNA expression above which the PAM50 subtype was probably normal and below which it was probably abnormal (Figure 8). Using a cut-off of 12 fpkm-uq, sensitivity was 45%, specificity was 99% positive predictive value=78%, negative predictive value=97%, and accuracy=97% (p<0.001).
VWF mRNA expression according to the Prosigna Breast Cancer Prognostic Gene Signature Assay (PAM50) subtype in 840 samples, all histologies. Tukey's honest significant difference post hoc test showed that VWF mRNA expression of the normal subtype was significantly greater than that of the other subtypes (p<0.001). LumA: Luminal A; LumB: luminal B; HER2: human epidermal growth factor receptor 2.
Discussion
Armand Trousseau first reported the relationship of malignant tumors and coagulation in 1865. Trousseau diagnosed the syndrome in himself 2 years later, dying soon afterward of gastric cancer. Recent findings suggest that genetic pathways within tumor cells might trigger thrombotic phenomena, thereby worsening prognosis (16).
Von Willebrand disease (vWD) is the most common hereditary coagulation disorder in humans, and results from a deficiency in the quality or quantity of VWF. vWD presents with bleeding diathesis: easy bruising, frequent nosebleeds, and bleeding gums. Women may experience menorrhagia and hemorrhage during parturition (17).
Some studies have found that VWF might promote cancer dissemination (6). Tumor cells actively interact with coagulation cascade factors such as VWF, thereby enhancing their metastatic capacity (6). A survey of 92 patients with vWD from 54 Italian hemophilia treatment centers revealed that many cancers occur in people with vWD (106 carcinomas) (18). In addition, cancer-cell-derived VWF mediates gastric cancer metastasis, and VWF may be a new therapeutic target (19). Other studies indicate that VWF might protect against tumor dissemination [reviewed in (6)]. In patients with breast cancer, elevated serum VWF is associated with advanced stage and tumor progression (20).
Receiver operating characteristic curve of the data used in Figure 7 in order to determine the optimum cut-off for von Willebrand factor (VWF) mRNA expression above which the Prosigna Breast Cancer Prognostic Gene Signature Assay (PAM50) subtype was probably normal and below which it was probably abnormal determine. Area under curve=0.784 (95% confidence intervaI=0.682-0.887), standard error=0.052, p<0.001. Using a cut-off of 12 sensitivity was 45%, specificity was 99%. positive predictive value=78%, negative predictive value=97% accuracy 97%.
As noted above, patients with ILC had higher VWF mRNA expression than patients with IDC or other histology. Additional evidence of the involvement of VWF in ILC is the relationship of VWF to PTEN. PTEN is part of the ILC genomic profile (2, 3). PTEN and VWF were found to share nine co-occurrent alterations in breast cancer, and mRNA expression of the two genes significantly correlated.
The PTEN gene codes for an enzyme that acts as a tumor suppressor. The PTEN enzyme is a tyrosine phosphatase: it modifies other proteins and lipids by removing phosphate groups. PTEN tyrosine phosphatase signals cells to stop dividing and to undergo apoptosis. The PTEN tyrosine phosphatase helps control cell migration, the adhesion of cells to surrounding tissues, and angiogenesis. Additionally, PTEN tyrosine phosphatase maintains the stability of a cell's genetic information. All of these functions prevent uncontrolled cell growth and the formation of tumors (21). Loss of expression of the PTEN tyrosine phosphatase is associated with poor outcome in breast cancer (21).
Biological processes or pathways in cancer are often deregulated through different genes or by multiple different mechanisms. But cancer gene mutations usually do not occur at random. Mutations of certain cancer genes tend to co-occur, indicating that they may work in tandem to drive tumor formation and development (11). This appears to be the case with the nine co-occurring alterations of VWF and PTEN in ILC.
In contrast, mutations of other genes appear in a mutually exclusive fashion, suggesting that two genes may have highly similar downstream components. The concept of mutual exclusivity can be exploited to identify previously unknown mechanisms that contribute to oncogenesis and cancer progression. In mutual exclusivity, events in genes associated with a specific cancer tend to be mutually exclusive across a set of tumors – that is, each tumor is likely to have only one of the genetic events (11). The 25 mutually exclusive mutations in PTEN and PIK3CA are an example of mutual exclusivity presented here (22).
In a cohort of women with breast cancer with all histologies, recurrence and death of women with basal-like subtype was significantly higher than for those with normal (non-basal-like) subtype (23). GDC TCGA data analyzed here demonstrated the same relationship of basal-like subtype to early recurrence.
PAM50 intrinsic subtypes are not equally distributed across breast cancer histologies (24), and VWF is not one of the genes analyzed in PAM50, although VWF RNA data we analyzed here is part of the GDC TCGA breast cancer cohort. Our finding of significantly higher VWF mRNA expression in those with normal (non-basal-like) subtype suggests that serum VWF measurement might complement or corroborate PAM50 results. VWF and PAM50 results, both suggesting low risk of recurrence, might make the decision whether to give chemotherapy easier, especially if VWF was an independent predictor.
Footnotes
Authors' Contributions
All contributed equally in data collection, data analysis, manuscript writing, and revision.
This article is freely accessible online.
Conflicts of Interest
None declared.
- Received January 28, 2019.
- Revision received February 27, 2019.
- Accepted March 12, 2019.
- Copyright© 2019, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved