Abstract
Dysregulation of tyrosine kinase receptor (RTK) signaling pathways play important roles in glioblastoma (GBM). However, therapies targeting these signaling pathways have not been successful, partially because of drug resistance. Increasing evidence suggests that tumor heterogeneity, more specifically, GBM-associated stem and endothelial cell heterogeneity, may contribute to drug resistance. In this perspective article, we introduce a high-throughput, quantitative approach to profile plasma membrane RTKs on single cells. First, we review the roles of RTKs in cancer. Then, we discuss the sources of cell heterogeneity in GBM, providing context to the key cells directing resistance to drugs. Finally, we present our provisionally patented qFlow cytometry approach, and report results of a “proof of concept” patient-derived xenograft GBM study.
Introduction
GBMs are the most frequent and lethal malignant primary adult brain tumor (Yadav et al., 2018), which presents a critical need to develop new therapeutics. Addressing the dysregulation of RTK signaling pathways offers promise in overcoming GBM lethality (Hawkins-Daarud et al., 2013; Cloughesy et al., 2014; Smith et al., 2016; Massey et al., 2018). RTK dysfunction has been observed in GBM, where these pathways are correlated with tumor cell proliferation (Johnson et al., 2012; Furnari et al., 2015), angiogenesis (Plate et al., 1994; Kuczynski et al., 2011), tumor invasiveness (Giannini et al., 2005; Sangar et al., 2014), and resistance to therapy (Murat et al., 2008; Lu and Bergers, 2013; Popescu et al., 2015). Moreover, these pathways are popular targets for small-molecule inhibitors (Rich and Bigner, 2004; Candolfi et al., 2011). Unfortunately, the clinical benefit of these targeted therapies is limited by drug resistance (De Witt Hamer, 2010; Szopa et al., 2017).
Increasing evidence suggests that drug resistance may be attributed to tumor heterogeneity (variations within an individual tumor) (Saunders et al., 2012; Furnari et al., 2015; Qazi et al., 2017). For example, a landmark study identified tumor subpopulations resistant to therapy prior to treatment by sequencing 4,645 single cells from 19 melanoma patients. This thorough analysis was enabled by single-cell technology, and may have been overlooked with ensemble sequencing (Tirosh et al., 2016). Additionally, a single-cell analysis of patient-derived xenografts (PDXs) of GBM39 also found higher heterogeneity in resistant tumors than in responsive tumors (Shi et al., 2012). In line with these single-cell measurements, we previously discovered, measured, and statistically described heterogeneity in breast cancer xenografts by quantifying vascular endothelial growth factor plasma membrane receptor (VEGFR) concentrations at the single-cell level (Imoukhuede and Popel, 2014). When we combined this quantitative analysis with computational modeling, we arrived at the prediction that tumors having “high” concentrations of plasma membrane VEGFR1 could be resistant to anti-VEGF drugs (angiogenesis inhibitors) (Weddell and Imoukhuede, 2014). Clinical work supports this prediction for colorectal cancer (Weickhardt et al., 2015), and application of this quantification and prediction should offer a new paradigm for biomarker discovery in cancer medicine.
To address the need for quantitative, single-cell analysis of GBM heterogeneity, we apply our optimized and provisionally patented VEGFR quantitative flow (qFlow) cytometry approach to GBM (Imoukhuede and Popel, 2011, 2014; Imoukhuede et al., 2013; Weddell and Imoukhuede, 2014; Chen et al., 2015, 2017; Imoukhuede and Chen, 2018). We describe expanded measurement to several RTKs critical to tumor development. To provide further context we, briefly, review the roles of RTKs in cancer and present connections between RTK heterogeneity and drug resistance. We then present our approach, qFlow cytometry, and report promising findings of a “proof of concept” PDX GBM study.
Roles of RTKs in cancer
RTKs are widely expressed transmembrane proteins (Cadena and Gill, 1992; Ferrara et al., 2003). Upon ligand binding, they are activated via canonical (Mac Gabhann and Popel, 2007; Sarabipour and Hristova, 2016) and non-canonical (Steinkamp et al., 2014; Chen et al., 2015; Pennock et al., 2016; Mamer et al., 2017) ligand-induced dimerization and tyrosine phosphorylation mechanisms. Importantly, unligated receptors can dimerize (Ruch et al., 2007; Chung et al., 2010; Low-Nam et al., 2011; Lin et al., 2012; Comps-Agrar et al., 2015; King et al., 2016; Sarabipour and Hristova, 2016; Sarabipour et al., 2016) and signal (Wu et al., 2010; Sarabipour et al., 2016; Kazlauskas, 2017), although ligand binding stabilizes the dimeric receptor structure. These receptor-initiated signaling events regulate cell survival, proliferation, differentiation, and motility (Hubbard and Miller, 2007; Volinsky and Kholodenko, 2013).
VEGFRs are upregulated in many cancers (Ferrara, 2002; Kut et al., 2007; Mac Gabhann and Popel, 2008). Signals through endothelial VEGFRs and the neuropilin (NRP) co-receptors (Imoukhuede and Popel, 2011, 2012, 2014; Imoukhuede et al., 2013; Gelfand et al., 2014) induce the sprouting angiogenic hallmarks of cell proliferation and cell migration (Simons et al., 2016). These sprouting angiogenesis hallmarks also sustain tumor growth and enable tumor metastasis (Hanahan and Weinberg, 2011; Shibuya, 2014). VEGF and other pro-angiogenic factors, may also regulate vascular growth and regression in tumors that co-opt pre-existing blood vessels (Holash et al., 1999; Jayson et al., 2016; Kuczynski et al., 2016).
In addition to these canonical pathways, cross-family signaling may also affect tumor vascularization. In this paradigm, ligands from one growth factor family bind to and signal through receptor(s) of another family. For instance, we have shown VEGF-mediated downregulation of PDGFRs (Chen et al., 2015), and discovered that both VEGF–PDGFR binding and PDGF–VEGFR binding is high affinity (Mamer et al., 2017). Other cross-family studies have identified VEGF–PDGFR binding and signaling (Ball et al., 2007; Pennock et al., 2016) and VEGFR–PDGFR dimerization in tumor associated pericytes (Greenberg et al., 2008). Altogether, these canonical and cross-family RTK mechanisms suggest several possible receptor activation landscapes that can contribute to tumor growth and drug resistance.
Gbm-associated cell heterogeneity: stem and endothelial
An accepted origin of tumor heterogeneity involves clonal evolution; an reiterative process of genetic mutation, clonal selection, and expansion, which drives the growth of single cancer cells into heterogeneous tumor masses (Greaves and Maley, 2012; Greaves, 2015; McGranahan and Swanton, 2017). In addition to cancer cells, other cell types within the tumor may also differentiate or transition as tumor develops. Some such cells include: tumor-associated fibroblasts, macrophages/monocytes, endothelial cells (ECs), and stem cells (Saunders et al., 2012). Here, we describe glioblastoma stem cells (GSCs) and ECs, which we focus on in our pilot study.
GSCs are an important tumor cell component, because despite their small number (~0.5–10%) (Pallini et al., 2011), GSCs are more resistant to radiotherapy and chemotherapy than other cancer cells (Schonberg et al., 2014; Seymour et al., 2015). Furthermore, their resistance can amplify tumor heterogeneity, because they have self-renewing and tumor-initiating capabilities (Lathia et al., 2015). GSCs are often identified by CD133 (Mak et al., 2011), which is associated with poor prognosis in a number of tumor types. There is controversy surrounding the usage of CD133 as a GSC marker (Golebiewska et al., 2013; Seymour et al., 2015; Bradshaw et al., 2016). Early studies showed a subpopulation of GBM cells expressing CD133 were able to form tumors (Singh et al., 2004) and further studies showed subpopulations of CD133− cells were also able to form tumors in vivo (Beier et al., 2007). While these studies do not negate the possible role of CD133 in identifying GSCs, they do highlight the importance of heterogeneity and the need for additional markers. Therefore, establishing a “barcode” of RTK plasma membrane concentrations on GSCs may help to identify novel markers, aiding in the isolation and understanding of these stem cells.
ECs, the primary structural unit of the vasculature, are an important contributor to GBM development. Unlike normal vessels, tumor vasculature is leaky, tortuous, and dilated (Jain, 2005; Aird, 2009). In addition to typical tumor vascular pathological features, brain tumor vasculature exhibits the loss of the important blood-brain-barrier feature of tight EC-EC junctions when tumor size grows beyond 1–2 mm in diameter (Jain et al., 2007). The close interaction between tumors and tumor vessels, and the observation of extensive EC heterogeneity supports the need for profiling tumor-associated ECs.
A paradigm shift in single-cell technologies: from gene-centric to proteomics
Studies characterizing GBM heterogeneity primarily focus on genetic and transcriptomic profiling (Verhaak et al., 2010; Snuderl et al., 2011; Dunn et al., 2012; Szerlip et al., 2012; Brennan et al., 2013; Patel et al., 2014; Ellis et al., 2015), which does not always correlate with functional changes (Simonson and Schnitzer, 2007; Feng et al., 2009; Taniguchi et al., 2010). Moreover, multiple studies show discordance between sequence data and protein expression in GBM, particularly with regards to epidermal growth factor receptor (EGFR) (Brennan et al., 2009) and PDGFR (Hermanson et al., 1992) gene vs. protein expression. Because proteins are the effectors of signaling toward functional response (Grecco et al., 2011; Imoukhuede et al., 2013; Chen et al., 2017), there is a need for increased protein-based, functional measurements.
qFlow cytometry offers a powerful tool for protein-based, single-cell measurements. It applies fluorescent calibration to traditional flow cytometry, converting signal to absolute protein concentrations (Lyer et al., 1997; Lee-Montiel and Imoukhuede, 2013; Chen et al., 2017). Absolute protein quantification allows detection of variations in proteins across published studies, tissues, replicates, and instrument settings (Wheeless et al., 1989; Rocha-Martins et al., 2012; Baumgartner et al., 2013; Nguyen et al., 2013; Vigelsø et al., 2015). Moreover, qFlow cytometry advances systems biology, providing the quantitative data needed for computational studies (Chen et al., 2014; Weddell and Imoukhuede, 2018). For example, using qFlow cytometry coupled with systems biology, we predicted that anti-VEGF efficacy depends on tumor endothelial VEGFR1 plasma membrane concentrations (Weddell and Imoukhuede, 2014). Furthermore, a receptor-internalization computational model recently predicted that small increases in plasma membrane RTK concentrations (< 1,000 receptors/cell) may double nuclear-based RTK signaling (Weddell and Imoukhuede, 2017), which further implicates RTK concentrations as a determinant of signal transduction. These predictions were only possible with the accurate experimental data offered by qFlow cytometry.
A new approach for examining GBM heterogeneity
We performed a “proof of concept” qFlow cytometry study on a PDX, GBM39 (Figure 1). GBM39 is known for EGFRvIII and low invasiveness, in vivo (Johnson et al., 2012; Wei et al., 2016). The xenograft was established with tumor tissue from patients undergoing surgical treatment at Mayo Clinic, Rochester, MN. Multiple studies characterize these PDX models and report maintenance of patient morphologic and molecular characteristics including EGFR amplification as well as tumor invasiveness (Giannini et al., 2005; Sarkaria et al., 2007).
Figure 1
Following dissociation, PDX cells were stained with Sytox Blue (a live/dead cell stain), CD45 (Patenaude et al., 2010), CD34 (Soares et al., 2007; Moghadam et al., 2015), and CD133 (Singh et al., 2004; Calabrese et al., 2007; Molina et al., 2014; Naujokat, 2014; Soeda et al., 2015) fluorophore-conjugated antibodies that target EC-like cells and GSCs, respectively (Figure 1). This labeling scheme excludes both dead cells and hematopoietic cells and enables identification of human tumor EC-like cells (hCD34+), mouse tumor EC-like cells (mCD34+), and GSCs (hCD133+) from the live CD45− pool (Figure 2A). To obtain reliable data, we obtained fluorescence signals from 2 to 3 samples/RTK with 10,000–35,000 live single cells collected per sample. As expected, the bulk GBM39 PDX sample was primarily non-EC, non-GSC cells (62.46%). In addition, we found 6-fold higher mouse tumor EC-like cells than human tumor EC-like cells (Figure 2B). This quantification aligned with prior studies of GBM xenograft showing ~7.1% EC population (CD45−CD31+CD34+). Consistent with our quantification of GSCs, a primary human study of 37 patients reported a range of 0.5–10% (Pallini et al., 2011), when identifying GSCs using the CD133 marker.
Figure 2
We labeled and screened 9 plasma membrane RTKs on these cells, which included two established GBM biomarkers, EGFR and insulin-like growth factor receptor (IGFR) (Sangar et al., 2014), and angiogenic signaling biomarkers: VEGFRs, PDGFRs, NRP1, and Tie2 (Carmeliet and Jain, 2000, 2011; Ferrara, 2002; Ferrara and Kerbel, 2005; Dudley, 2012). Using qFlow cytometry and statistical models, we quantitatively characterized GBM39 PDX via four patented metrics (Figure 1): cell composition, ensemble RTK concentration, cell-by-cell analysis with Gaussian mixture modeling, and heterogeneity analysis (Imoukhuede and Chen, 2018).
Percentage of gated cell populations were exported using FlowJo software (TreeStar). Ensemble RTK concentrations and cell-by-cell analysis were performed as previously described (Chen et al., 2015, 2017). We then applied Gaussian mixture modeling to identify log-normal sub-populations within each distribution, described by its mean, standard deviation, and density. We reduced the chance of overfitting the subpopulations by using Bayesian Information Criterion (BIC) (Raftery, 1995; Huedo-Medina et al., 2006). A detailed description of heterogeneity quantification is provided in section Quantification of Cell-RTK Heterogeneity.
Human tumor EC-like cells have high EGFR and IGFR on plasma membrane
EGFR and IGFR are expressed on tumor cells and contribute to tumor progression. Interestingly, the human tumor EC-like population had high plasma membrane EGFR and IGFR concentrations (~21,000/cell and ~20,000/cell, respectively) (Figure 2C), consistent with qualitative findings of higher EGFR on breast carcinoma-derived ECs compared to normal ECs (Amin et al., 2006). Our results of high EGFR on human tumor EC-like cells from GBM39 is also consistent with results of clinical GBM samples (Soda et al., 2011).
The mixture modeling revealed that 8% of human tumor EC-like subpopulations had a ~12-fold higher membrane localization of EGFRs than average. We found a similar pattern for IGFRs in human tumor EC-like subpopulations. Together, the ensemble-averaged data and the mixture modeling indicated significant plasma membrane localization of EGFR and IGFR on human tumor EC-like cells. High concentrations of EGFR and IGFR suggest an opportunity for targeted inhibition, which could be a mechanism for disrupting tumor vessels on GBMs with a similar profile.
Mouse tumor EC-like cells have similar plasma membrane VEGFR concentrations as healthy mouse ECs from skeletal muscle
VEGFRs are key regulators of tumor angiogenesis, so their quantification can offer insight into the tumor vasculature. Furthermore, as biomarkers of vasculature, these receptors have been proposed as diagnostic biomarkers of anti-angiogenic drug efficacy (Lambrechts et al., 2013; Wehland et al., 2013) with computational (Weddell and Imoukhuede, 2014) and clinical (Weickhardt et al., 2015) support to their use. We found that VEGFR1 and VEGFR2 had similar concentrations and ratios on mouse tumor EC-like cells (~3,100 VEGFR1/cell and ~1,000 VEGFR2/cell) as on healthy ECs obtained from mouse skeletal muscle (Imoukhuede and Popel, 2012) (Figure 2E). This finding of a low VEGFR2:VEGFR1 ratio aligns with a previous study on breast cancer xenografts (Imoukhuede and Popel, 2014); however, the receptor abundance we report here is much lower. These findings of EC-like cells from GBM39 having VEGFRs at levels similar to normal mouse skeletal muscle ECs suggests a need for further quantification of normal brain ECs VEGFR concentrations to establish tissue standards. Similarly, it suggests a need to examine other GBM specimens to identify whether this is a property of co-opted vessels or specific to this GBM strain.
We analyzed the human tumor EC-like population (5.20% of the population, Figure 2B), which should reflect the original tumor vessels from the patient (Figure 2C). We found similar plasma membrane VEGFR1 and VEGFR2 ratios (~3,600 VEGFR1/cell & ~5,800 VEGFR2/cell) as previous reports in vitro (Imoukhuede and Popel, 2011; Chen et al., 2015). However, these data show that not all tumors have the same concentrations or ratios of plasma membrane VEGFRs on their endothelium. Importantly, tumor EC-like cells display much greater heterogeneity than normal ECs with subpopulations that have high concentrations of VEGFRs. Indeed, cell-by-cell analysis and mixture modeling of human tumor EC-like cells reveals the existence of a high-VEGFR1 subpopulation (~10%) with ~41,000 VEGFR1/cell, while the highest VEGFR2 subpopulation is ~18,500 VEGFR2/cell, comprising ~35% of the total human tumor EC-like population (Figures 2D,G). The difference in VEGFR2:VEGFR1 ratio and receptor concentrations between human and mouse tumor EC-like population shows a significant level of endothelial heterogeneity. Such data may enable correlations between these tumor vessel regulators and anti-angiogenic drug efficacy.
Plasma membrane PDGFRs localize on tumor EC-like cells
PDGFRs serve important roles in supporting vasculature in tumor microenvironments (Andrae et al., 2008). We observed lower levels of PDGFRs on human tumor EC-like cell membranes than on mouse (Figures 2C,E). The cell-by-cell analysis and mixture modeling suggests that this ensemble average does not capture the subpopulations having high-PDGFR plasma membrane localization: 66 and 16% of mouse tumor EC-like cell membrane had ~23,400 PDGFRα and ~19,800 PDGFRβ, respectively (Figure 2F). This significant heterogeneity may be attributed to the use of the CD34 marker to designate EC-like cells, because it is also found on stem cells/precursors, mast cells, and neurons (Nielsen and McNagny, 2008; Imoukhuede and Popel, 2014; AbuSamra et al., 2017). PDGFRα is also considered an important mesenchymal stem cell marker (Farahani and Xaymardan, 2015). So, the co-labeling of PDGFRα and CD34 suggests these cells may be mesenchymal stem cells (Aguirre et al., 2010).
If these CD34+PDGFR+ cells are endothelial, then our data correlates with studies finding PDGFRs on tumor ECs (Hermansson et al., 1988; Werner et al., 1990; Plate et al., 1992). PDGFR localization on ECs is controversial, because it is characteristic of mural cells and not of ECs (Heldin et al., 1981; Bowen-Pope and Ross, 1982; Kazlauskas and DiCorleto, 1985; Raines et al., 1991; Battegay et al., 1994; Marx et al., 1994). However, they have been observed on monolayer microvascular ECs, in vitro (Bar et al., 1989; Marx et al., 1994) and on angiogenic ECs that formed sprout and tubes in vitro (Battegay et al., 1994). If we subscribe to the canonical PDGFR localization understanding, then these tumor vessels induce “non-conventional” PDGFR localization patterns.
GSCs have little-to-no surface EGFR or IGFR
Multiple studies suggest that a higher degree of GSC “stemness” is associated with EGFR amplification (Mazzoleni et al., 2010; Liffers et al., 2015); however, we observed ~13-fold lower EGFRs on GSC plasma membranes compared to the bulk PDX cells (Figure 2C). This trend was also seen with IGFR (Figure 2C). The low membrane EGFR concentrations on GSCs is concerning, given reports that EGFR signaling is necessary for GSC proliferation and tumor-sphere formation (Soeda et al., 2008; Griffero et al., 2009). Yet, this may explain the lower percentage of GSCs in the PDX sample (~0.9%) compared to the expected stem cell fraction (0.5–10%; Pallini et al., 2011). A possible explanation is that serially transplanted tumors can lose their EGFR overexpression, even in vivo (Liffers et al., 2015). Clearly, further investigation of both gene expression and protein quantification on other GBM PDX GSCs is necessary to understand their contribution to heterogeneity and drug-resistance.
Quantification of cell-RTK heterogeneity
To quantify heterogeneity of each cell subpopulation, we used two parameters: number of mixture components and Quadratic Entropy (QE). To quantitatively assess the number of subpopulations within each cell population, we fit each cell-by-cell RTK distribution with mixture models consisting of 1–9 log-normal Gaussian sub-distributions (mixture components); we then applied BIC as the criterion to select the mixture model with the lowest BIC. The number of mixture components is determined by how many log-normal Gaussian sub-distributions are in the mixture model. The number of mixture components, thus, is a measurement of cell heterogeneity. Generally, 1–2 mixture components are considered low heterogeneity (Chen et al., 2017), while more than two components is considered highly heterogeneous (Imoukhuede and Popel, 2014; Weddell and Imoukhuede, 2014).
Alternatively, QE requires equally spaced bins, here we chose 500 bins, from each cell-by-cell distribution (Figures 2D,F). QE then sums the weighted differences of the means between two bins (Rao, 1982; Pavoine and Dolédec, 2005; Zoltán, 2005). Thus, QE is a measurement of the increase in random variation in the cellular response. Because healthy ECs and human fibroblasts in vitro have shown QE within 0.2–0.7 (Chen et al., 2015), we describe QE < 0.7 as low heterogeneity and QE > 0.7 as high heterogeneity. QE provides a quantitative measure of the diversity of cellular phenotypes in cancer tissue sections for diagnostic applications (Potts et al., 2012) and drug discovery (Gough et al., 2014). Interestingly, human tumor EC-like cells showed lower QE and number of mixture components when compared to mouse tumor EC-like cells (Figure 2G). We suspect that the likely loss of human tumor-associated cells over time in a PDX model (Chao et al., 2017) may be the reason why human tumor EC-like cells present a more homogenous state than the mouse tumor EC-like cells.
Clinical implications of GBM heterogeneity
We envision that RTK quantification can identify ideal receptor targets across the bulk tumor specimen and on specific cell populations in the tumor. First, the ideal receptor target would be highly available (Rich and Bigner, 2004; Weis and Cheresh, 2011; Cloughesy et al., 2014): it would have high concentrations on a high percentage of bulk cells or specific cells. Next, the target RTK would exhibit low heterogeneity: it would have low QE in bulk cells or on the specific cell subpopulation (Jain et al., 2009; Heath et al., 2016). An ideal receptor target would also be highly specific to the tumor, which would manifest as higher receptor concentrations in the tumor vs. healthy tissue (Rich and Bigner, 2004).
Based on these guidelines, we offer possible targets on GBM39. If the goal is targeting tumor vessels, then VEGFR2 and PDGFRα are highly targetable: >70% target cells have > 6,000 VEGFR2 or PDGFRα/cell plasma membrane with QE = 0.20 or 0.32, respectively. Furthermore, they are likely targets, because they are more highly expressed in GBM specimens than health tissue (Chen et al., 2015): ~5-fold higher VEGFR2 and ~4-fold higher PDGFRα. Therefore, targeting VEGFR2 and PDGFRα should preferentially target the tumor.
Our work suggests that targeting EGFR and IGFR on tumors like GBM39 may not be effective by itself. Although, they have high concentrations on ~70–90% EC-like and non-EC-like GBM cells, their high GBM heterogeneity (QE = ~1.0) and high concentration on healthy tissue (2–2,000 × 103 EGFR/fibroblast or epithelial cell; 2.5 × 104 IGFR/NIH 3T3 mouse fibroblasts (Sorkin and Duex, 2010; Brennan et al., 2013; Weddell and Imoukhuede, 2017) may lower their targeting specificity, resulting in lower drug efficacy (Wheeler et al., 2010). Better drug delivery to the tumor site will likely improve targeting specificity without disrupting healthy tissue. An alternative strategy is to develop dual-inhibitors targeting both EGFR/IGFR and VEGFR2 to increase their specificity for tumor EC-like cells.
We believe our method can also identify cellular and molecular mechanisms underlying reduced response to drugs. For example, upregulation of alternative signaling pathways has been implicated in anti-VEGF drug resistance (Bergers and Hanahan, 2008; Shojaei and Ferrara, 2008). This mechanism of drug resistance is often accompanied by significant tumor heterogeneity (Snuderl et al., 2011; Szerlip et al., 2012; Lu and Bergers, 2013). Therefore, these alternative pathways may be overlooked in bulk studies if they are only present on small cell subpopulations. From this study, we suggest targeting RTKs that are localized on plasma membrane at high concentrations on small cell populations (< 10%) for combination therapy. For example, VEGFR1 and Tie2 on tumor ECs may become “alternative” RTKs for anti-VEGF treatment, because we found ~10% human tumor EC-like cell subpopulations had 41,000 VEGFR1 and ~8% had 65,700 Tie2 on the plasma membrane. Identifying alternative RTK pathways that contribute to resistance can provide tumor-specific drug targets for combination therapy.
Future direction in characterizing GBM heterogeneity
Our study of the GBM39 PDX model, arrived at 4 key findings and 2 recommendations: (1) tumor EC-like subpopulations have high concentrations of plasma membrane VEGFR1 and VEGFR2; (2) human vs. mouse tumor EC-like cells have inverted VEGFR2:VEGFR1 ratios; (3) tumor EC-like subpopulations have high plasma membrane EGFR, IGFR, and PDGFR concentrations; and (4) GSCs compose a low percentage of cells in the tumor and have little-to-no EGFRs and IGFRs on their plasma membranes.
Based on findings in this study and our RTK-targeting criteria, VEGFR2 or PDGFRα would be likely drug targets for GBM39. In addition, VEGFR1 and Tie2 are likely drug targets for combination therapy. The next step would be to test these targets in a GBM PDX model.
The results of this “proof of concept” study should be interpreted as such: it offers an approach for continued measurement of tumor samples, broadly, and GBM samples, specifically, with the GBM39 PDX sample as a first example. We present the novel method, qFlow cytometry, and show its application in characterizing GBM heterogeneity. Larger and well-powered samples are warranted to expand the current preliminary results, and to discover ideal drug targets and mechanisms underlying drug resistance.
Future opportunities for expanding this research lies in establishing protein concentration ranges on additional samples and continued development of biomimetic tumor models. Firstly, additional measurements of protein concentration on normal ECs and other cells would provide the baselines needed to compare to tumor. In establishing EC baselines, isolation of a pure EC population may be a challenge. Previous qFlow studies have identified ECs using both the CD34 and CD31 markers (Imoukhuede and Popel, 2012, 2014; Imoukhuede et al., 2013). However, it is important to note, that using multiple markers can bias cell collection: CD34 is a progenitor marker, so its use biases selection from more mature cells. Whereas, CD31 is a mature cell marker that is found on ECs, platelets, natural killer cells, monocytes, macrophages, and among other cells (Liu and Shi, 2012), so its use can lead to sample impurity. Here, we chose to bias toward progenitor-like ECs; however, expanded studies may determine if protein concentrations correlate with marker presentation (e.g., identifying whether progenitor-like cells having higher or lower protein concentrations).
Another opportunity for advancement lies in our quantitative single-cell RTK mapping, moving toward multiplexed measurement of RTKs. Toward multiplexed quantification, Lee-Montiel et al., developed a quantum dot method for receptor labeling and calibration (Lee-Montiel and Imoukhuede, 2013; Lee-Montiel et al., 2015) that can be translated to qFlow cytometry. Another approach could be to adapt receptor quantification to mass cytometry (CyTOF) (Spitzer and Nolan, 2016). Such advancements would provide multi-RTK, multi-cell insight into tumor heterogeneity.
In conclusion, cancer research is experiencing a paradigm shift from ensemble analysis to cell-to-cell variability (Niepel et al., 2009; Hoppe et al., 2014; Dar and Weiss, 2018) because of the increasing evidence correlating drug resistance with tumor heterogeneity. The perspective and work that we present here offers sensitive methods for heterogeneity characterization in tumors that should enable improved treatment. We believe that continued quantification of single-cell receptor heterogeneity is a new frontier that will offer significant clinical impact.
Statements
Author contributions
SC and PI designed the study. BH provided samples and contributed critical insight into GBM and discussion on the studies. SC performed the experiments and quantified RTK plasma membrane concentrations. SC, TL, and PI conducted and discussed data analysis. SC, TL, and PI prepared the manuscript.
Acknowledgments
This work was supported by grants from the National Science Foundation (1512598 and 1653925), and the American Heart Association (16SDG26940002). Research reported in this publication was also supported by the National Cancer Institute of the National Institutes of Health under Award Number R01 CA197488 (BACH). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. We would like to thank members of BH's lab, Dr. Sara Pedron and Jee-Wei Emily Chen at University of Illinois at Urbana-Champaign, and Dr. Jann Sarkaria at Mayo Clinic (Rochester, MN) for the generous supply of GBM PDX samples.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AbuSamraD. B.FajrA. A.AsmaS. A.HebaM. J. A.CheeJ. C.AymanF. A.et al. (2017). Not just a marker: cd34 on human hematopoietic stem/progenitor cells dominates vascular selectin binding along with CD44. Blood Adv.1, 2799–2816. 10.1182/bloodadvances.2017004317
2
AguirreA.PlanellJ. A.EngelE. (2010). Dynamics of bone marrow-derived endothelial progenitor cell/mesenchymal stem cell interaction in co-culture and its implications in angiogenesis. Biochem. Biophys. Res, Commun. 400, 284–291. 10.1016/j.bbrc.2010.08.073
3
AirdW. C. (2009). Molecular heterogeneity of tumor endothelium. Cell Tissue Res.335, 271–281. 10.1007/s00441-008-0672-y
4
AminD. N.HidaK.BielenbergD. R.KlagsbrunM. (2006). Tumor endothelial cells express epidermal growth factor receptor (EGFR) but not ErbB3 and are responsive to EGF and to EGFR |. Cancer Res.66, 2173–2180. 10.1158/0008-5472.CAN-05-3387
5
AndraeJ.GalliniR.BetsholtzC. (2008). Role of platelet-derived growth factors in physiology and medicine. Genes Dev.22, 1276–1312. 10.1101/gad.1653708
6
BallS. G.ShuttleworthC. A.KieltyC. M. (2007). Vascular endothelial growth factor can signal through platelet-derived growth factor receptors. J. Cell Biol.177, 489–500. 10.1083/jcb.200608093
7
BarR. S.BoesM.BoothB. A.DakeB. L.HenleyS.HartM. N. (1989). The effects of platelet-derived growth factor in cultured microvessel endothelial cells. Endocrinology124, 1841–1848. 10.1210/endo-124-4-1841
8
BattegayE. J.RuppJ.Iruela-ArispeL.SageE. H.PechM. (1994). PDGF-BB modulates endothelial proliferation and angiogenesis in vitro via PDGF beta-receptors. J. Cell Biol.125, 917–928. 10.1083/jcb.125.4.917
9
BaumgartnerR.UmlaufE.VeitingerM.GuterresS.RappoldE.BabelukR.et al. (2013). Identification and validation of platelet low biological variation proteins, superior to GAPDH, actin and tubulin, as tools in clinical proteomics. J. Proteomics94, 540–551. 10.1016/j.jprot.2013.10.015
10
BeierD.HauP.ProescholdtM.LohmeierA.WischhusenJ.OefnerP. J.et al. (2007). CD133(+) and CD133(-) glioblastoma-derived cancer stem cells show differential growth characteristics and molecular profiles. Cancer Res.67, 4010–4015. 10.1158/0008-5472.CAN-06-4180
11
BergersG.HanahanD. (2008). Modes of resistance to anti-angiogenic therapy. Nat. Rev. Cancer8, 592–603. 10.1038/nrc2442
12
Bowen-PopeD. F.RossR. (1982). Platelet-derived growth factor. II. Specific binding to cultured cells. J. Biol. Chem.257, 5161–5171.
13
BradshawA.WickremsekeraA.TanS. T.PengL.DavisP. F.ItinteangT. (2016). Cancer stem cell hierarchy in glioblastoma multiforme. Front. Surg.3:21. 10.3389/fsurg.2016.00021
14
BrennanC. W.VerhaakR. G. W.McKennaA.CamposB.NoushmehrH.SalamaS. R.ZhengS.et al. (2013). The somatic genomic landscape of glioblastoma. Cell155, 462–477. 10.1016/j.cell.2013.09.034
15
BrennanC.MomotaH.HambardzumyanD.OzawaT.TandonA.PedrazaA.et al. (2009). Glioblastoma subclasses can be defined by activity among signal transduction pathways and associated genomic alterations. PLoS ONE4:e7752. 10.1371/journal.pone.0007752
16
CadenaD. L.GillG. N. (1992). Receptor tyrosine kinases. FASEB J.6, 2332–2337. 10.1096/fasebj.6.6.1312047
17
CalabreseC.PoppletonH.KocakM.HoggT. L.FullerC.HamnerB.OhE. Y.et al. (2007). A perivascular niche for brain tumor stem cells. Cancer Cell11, 69–82. 10.1016/j.ccr.2006.11.020
18
CandolfiM.KroegerK. M.XiongW.LiuC.PuntelM.YagizK.MuhammadA. G.et al. (2011). Targeted toxins for glioblastoma multiforme: pre-clinical studies and clinical implementation. Anticancer. Agents Med. Chem.11, 729–738. 10.2174/187152011797378689
19
CarmelietP.JainR. K. (2000). Angiogenesis in cancer and other diseases. J. Nat.407, 249–257. 10.1038/35025220
20
CarmelietP.JainR. K. (2011). Molecular mechanisms and clinical applications of angiogenesis. Nature473, 298–307. 10.1038/nature10144
21
ChaoC.WidenS. G.WoodT. G.ZatarainJ. R.JohnsonP.GajjarA.GomezG.et al. (2017). Patient-derived xenografts from colorectal carcinoma: a temporal and hierarchical study of murine stromal cell replacement. Anticancer Res.37, 3405–3412. 10.21873/anticanres.11707
22
ChenS.AnsariA.SterrettW.HurleyK.KemballJ.WeddellJ. C.et al. (2014). Current state-of-the-art and future directions in systems biology. Progr. Commun. Sci.1, 12–26.
23
ChenS.GuoX.ImarenezorO.ImoukhuedeP. I. (2015). Quantification of VEGFRs, NRP1, and PDGFRs on endothelial cells and fibroblasts reveals serum, intra-family ligand, and cross-family ligand regulation. Cell. Mol. Bioeng.8, 383–403. 10.1007/s12195-015-0411-x
24
ChenS.WeddellJ. C.GuptaP.ConardG.ParkinJ.ImoukhuedeP. I. (2017). qFlow cytometry-based receptoromic screening: a high-throughput quantification approach informing biomarker selection and nanosensor development. Methods Mol. Biol.1570, 117–138. 10.1007/978-1-4939-6840-4_8
25
ChungI.AkitaR.VandlenR.ToomreD.SchlessingerJ.MellmanI. (2010). Spatial control of EGF receptor activation by reversible dimerization on living cells. Nature464, 783–787. 10.1038/nature08827
26
CloughesyT. F.CaveneeW. K.MischelP. S. (2014). Glioblastoma: from molecular pathology to targeted treatment. Ann. Rev. Pathol. Mech. Dis.9, 1–25. 10.1146/annurev-pathol-011110-130324
27
Comps-AgrarL.DunsheeD. R.EatonD. L.SonodaJ. (2015). Unliganded fibroblast growth factor receptor 1 forms density-independent dimers. J. Biol. Chem.290, 24166–24177. 10.1074/jbc.M115.681395
28
DarR. D.WeissR. (2018). Perspective: engineering noise in biological systems towards predictive stochastic design. APL Bioeng.2:020901. 10.1063/1.5025033
29
De Witt HamerP. C. (2010). Small molecule kinase inhibitors in glioblastoma: a systematic review of clinical studies. Neurooncology12, 304–316. 10.1093/neuonc/nop068
30
DudleyA. C. (2012). Tumor endothelial cells. Cold Spring Harb. Perspect. Med.2:a006536. 10.1101/cshperspect.a006536
31
DunnG. P.RinneM. L.WykoskyJ.GenoveseG.QuayleS. N.DunnI. F.AgarwallaP. K.et al. (2012). Emerging insights into the molecular and cellular basis of glioblastoma. Genes Dev.26, 756–784. 10.1101/gad.187922.112
32
EllisH. P.GreensladeM.PowellB.SpiteriI.SottorivaA.KurianK. M. (2015). Current challenges in glioblastoma: intratumour heterogeneity, residual disease, and models to predict disease recurrence. Front. Oncol.5:251. 10.3389/fonc.2015.00251
33
FarahaniR. M.XaymardanM. (2015). Platelet-derived growth factor receptor alpha as a marker of mesenchymal stem cells in development and stem cell biology. Stem Cells Int.2015:362753. 10.1155/2015/362753
34
FengY.MitchisonT. J.BenderA.YoungD. W.TallaricoJ. A. (2009). Multi-parameter phenotypic profiling: using cellular effects to characterize small-molecule compounds. Nat. Rev. Drug Discov.8, 567–578. 10.1038/nrd2876
35
FerraraN. (2002). VEGF and the quest for tumour angiogenesis factors. Nat. Rev. Cancer2, 795–803. 10.1038/nrc909
36
FerraraN.KerbelR. S. (2005). Angiogenesis as a therapeutic target. Nature438, 967–974. 10.1038/nature04483
37
FerraraN.GerberH.LeCouterJ. (2003). The biology of VEGF and its receptors. Nat. Med.9, 669–676. 10.1038/nm0603-669
38
FurnariF. B.CloughesyT. F.CaveneeW. K.MischelP. S. (2015). Heterogeneity of epidermal growth factor receptor signalling networks in glioblastoma. Nat. Rev. Cancer15, 302–310. 10.1038/nrc3918
39
GelfandM. V.HaganN.TataA.OhW.LacosteB.KangK.et al. (2014). Neuropilin-1 functions as a VEGFR2 co-receptor to guide developmental angiogenesis independent of ligand binding. ELife3:e03720. 10.7554/eLife.03720
40
GianniniC.SarkariaJ. N.SaitoA.UhmJ. H.GalanisE.CarlsonB. L.et al. (2005). Patient tumor EGFR and PDGFRA gene amplifications retained in an invasive intracranial xenograft model of glioblastoma multiforme. Neurooncology7, 164–176. 10.1215/S1152851704000821
41
GolebiewskaA.BougnaudS.StieberD.BronsN. H. C.VallarL.HertelF.et al. (2013). Side population in human glioblastoma is non-tumorigenic and characterizes brain endothelial cells. Brain136 (Pt 5), 1462–1475. 10.1093/brain/awt025
42
GoughA. H.ChenN.ShunT. Y.LezonT. R.BoltzR. C.ReeseC. E.WagnerJ.et al. (2014). Identifying and quantifying heterogeneity in high content analysis: application of heterogeneity indices to drug discovery. PLoS ONE9:e102678. 10.1371/journal.pone.0102678
43
GreavesM. (2015). Evolutionary determinants of cancer. Cancer Discov.5, 806–820. 10.1158/2159-8290.CD-15-0439
44
GreavesM.MaleyC. C. (2012). Clonal evolution in cancer. Nature481, 306–313. 10.1038/nature10762
45
GreccoH. E.SchmickM.BastiaensP. I. H. (2011). Signaling from the living plasma membrane. Cell144, 897–909. 10.1016/j.cell.2011.01.029
46
GreenbergJ. I.ShieldsD. J.BarillasS. G.AcevedoL. M.MurphyE.HuangJ.ScheppkeL.et al. (2008). A Role for VEGF as a negative regulator of pericyte function and vessel maturation. Nature456, 809–813. 10.1038/nature07424
47
GrifferoF.DagaA.MarubbiD.CapraM. C.MelottiA.PattarozziA.GattiM.et al. (2009). Different response of human glioma tumor-initiating cells to epidermal growth factor receptor kinase inhibitors. J. Biol. Chem.284, 7138–7148. 10.1074/jbc.M807111200
48
HanahanD.WeinbergR. A. (2011). Hallmarks of cancer: the next generation. Cell144, 646–674. 10.1016/j.cell.2011.02.013
49
Hawkins-DaarudA.RockneR. C.AndersonA. R. A.SwansonK. R. (2013). Modeling tumor-associated edema in gliomas during anti-angiogenic therapy and its impact on imageable tumor. Front. Oncol.3:66. 10.3389/fonc.2013.00066
50
HeathJ. R.RibasA.MischelP. S. (2016). Single-cell analysis tools for drug discovery and development. Nat. Rev. Drug Discov.15, 204–216. 10.1038/nrd.2015.16
51
HeldinC. H.WestermarkB.WastesonA. (1981). Specific receptors for platelet-derived growth factor on cells derived from connective tissue and glia. Proc. Natl. Acad. Sci. U.S.A.78, 3664–3668. 10.1073/pnas.78.6.3664
52
HermansonM.FunaK.HartmanM.Claesson-WelshL.HeldinC. H.WestermarkB.et al. (1992). Platelet-derived growth factor and its receptors in human glioma tissue: expression of messenger RNA and protein suggests the presence of autocrine and paracrine loops. Cancer Res.52, 3213–3219.
53
HermanssonM.NistérM.BetsholtzC.HeldinC. H.WestermarkB.FunaK. (1988). Endothelial cell hyperplasia in human glioblastoma: coexpression of mRNA for platelet-derived growth factor (PDGF) B chain and PDGF receptor suggests autocrine growth stimulation. Proc. Natl. Acad. Sci. U.S.A.85, 7748–7752. 10.1073/pnas.85.20.7748
54
HolashJ.MaisonpierreP. C.ComptonD.BolandP.AlexanderC. R.ZagzagD.et al. (1999). Vessel cooption, regression, and growth in tumors mediated by angiopoietins and VEGF. Science284, 1994–1998. 10.1126/science.284.5422.1994
55
HoppeP. S.CoutuD. L.SchroederT. (2014). Single-cell technologies sharpen up mammalian stem cell research. Nat. Cell Biol.16, 919–927. 10.1038/ncb3042
56
HubbardS. R.MillerW. T. (2007). Receptor tyrosine kinases: mechanisms of activation and signaling. Curr. Opin. Cell Biol.19, 117–123. 10.1016/j.ceb.2007.02.010
57
Huedo-MedinaT. B.Sánchez-MecaJ.Marín-MartínezF.BotellaJ. (2006). Assessing heterogeneity in meta-analysis: Q statistic or I2 index?Psychol. Methods11, 193–206. 10.1037/1082-989X.11.2.193
58
ImoukhuedeP. I.ChenS. (2018). Characterizing Heterogeneity via Single-Cell Receptor Quantification. UIUC2018-032-01(PRO).
59
ImoukhuedeP. I.PopelA. S. (2011). Quantification and cell-to-cell variation of vascular endothelial growth factor receptors. Exp. Cell Res.317, 955–965. 10.1016/j.yexcr.2010.12.014
60
ImoukhuedeP. I.PopelA. S. (2012). Expression of VEGF receptors on endothelial cells in mouse skeletal muscle. PLoS ONE7:e44791. 10.1371/journal.pone.0044791
61
ImoukhuedeP. I.PopelA. S. (2014). Quantitative fluorescent profiling of VEGFRs reveals tumor cell and endothelial cell heterogeneity in breast cancer xenografts. Cancer Med.3, 225–244. 10.1002/cam4.188
62
ImoukhuedeP. I.DokunA. O.AnnexB. H.PopelA. S. (2013). Endothelial cell-by-cell profiling reveals the temporal dynamics of VEGFR1 and VEGFR2 membrane localization after murine hindlimb ischemia. Am. J. Physiol. Heart Circ. Physiol.304, H1085–H1093. 10.1152/ajpheart.00514.2012
63
JainR. K. (2005). Normalization of tumor vasculature: an emerging concept in antiangiogenic therapy. Science307, 58–62. 10.1126/science.1104819
64
JainR. K.DudaD. G.WillettC. G.SahaniD. V.ZhuA. X.LoefflerJ. S.et al. (2009). Biomarkers of response and resistance to antiangiogenic therapy. Nat. Rev. Clin. Oncol.6, 327–338. 10.1038/nrclinonc.2009.63
65
JainR. K.TomasoE. D.DudaD. G.LoefflerJ. S.SorensenA. G.BatchelorT. T. (2007). Angiogenesis in brain tumours. Nat. Rev. Neurosci.8, 610–622. 10.1038/nrn2175
66
JaysonG. C.KerbelR.EllisL. M.HarrisA. L. (2016). Antiangiogenic therapy in oncology: current status and future directions. Lancet388, 518–529. 10.1016/S0140-6736(15)01088-0
67
JohnsonH.RosarioA. M. D.BrysonB. D.SchroederM. A.SarkariaJ. N.WhiteF. M. (2012). Molecular characterization of EGFR and EGFRvIII signaling networks in human glioblastoma tumor xenografts. Mol. Cell. Proteomics11, 1724–1740. 10.1074/mcp.M112.019984
68
KazlauskasA. (2017). PDGFs and their receptors. Gene614, 1–7. 10.1016/j.gene.2017.03.003
69
KazlauskasA.DiCorletoP. E. (1985). Cultured endothelial cells do not respond to a platelet-derived growth-factor-like protein in an autocrine manner. Biochim. Biophys. Acta846, 405–412. 10.1016/0167-4889(85)90013-8
70
KingC.StonemanM.RaicuV.HristovaK. (2016). Fully quantified spectral imaging reveals in vivo membrane protein interactions. Integr. Biol.8, 216–229. 10.1039/C5IB00202H
71
KuczynskiE. A.PattenS. G.CoomberB. L. (2011). VEGFR2 expression and TGF-β signaling in initial and recurrent high-grade human glioma. Oncology81, 126–134. 10.1159/000332849
72
KuczynskiE. A.YinM.Bar-ZionA.LeeC. R.ButzH.ManS.DaleyF.et al. (2016). Co-option of liver vessels and not sprouting angiogenesis drives acquired sorafenib resistance in hepatocellular carcinoma. J. Natl. Cancer Inst.108:djw030. 10.1093/jnci/djw030
73
KutC.GabhannF. M.PopelA. S. (2007). Where is VEGF in the body? A meta-analysis of VEGF distribution in cancer. Brit. J. Cancer97, 978–985. 10.1038/sj.bjc.6603923
74
LambrechtsD.LenzH.HaasS.CarmelietP.SchererS. J. (2013). Markers of response for the antiangiogenic agent bevacizumab. J. Clin. Oncol.31, 1219–1230. 10.1200/JCO.2012.46.2762
75
LathiaJ. D.MackS. C.Mulkearns-HubertE. E.ValentimC. L. L.RichJ. N. (2015). Cancer stem cells in glioblastoma. Genes Dev.29, 1203–1217. 10.1101/gad.261982.115
76
Lee-MontielF. T.ImoukhuedeP. I. (2013). Engineering quantum dot calibration standards for quantitative fluorescent profiling. J. Mater. Chem. B1, 6434. 10.1039/c3tb20904k
77
Lee-MontielF. T.LiP.ImoukhuedeP. I. (2015). Quantum dot multiplexing for the profiling of cellular receptors. Nanoscale7, 18504–18514. 10.1039/C5NR01455G
78
LiffersK.LamszusK.SchulteA. (2015). EGFR amplification and glioblastoma stem-like cells. Stem Cells Int.2015:427518. 10.1155/2015/427518
79
LinC. C.MeloF. A.GhoshR.SuenK. M.StaggL. J.KirkpatrickJ.et al. (2012). Inhibition of basal FGF receptor signaling by dimeric Grb2. Cell149, 1514–1524. 10.1016/j.cell.2012.04.033
80
LiuL.ShiG. P. (2012). CD31: beyond a marker for endothelial cells. Cardiovasc. Res.94, 3–5. 10.1093/cvr/cvs108
81
Low-NamS. T.LidkeK. A.CutlerP. J.RooversR. C.HenegouwenP. M. B. P.WilsonB. S.et al. (2011). ErbB1 dimerization is promoted by domain co-confinement and stabilized by ligand binding. Nat. Struct. Mol. Biol.18, 1244–1249. 10.1038/nsmb.2135
82
LuK. V.BergersG. (2013). Mechanisms of evasive resistance to anti-VEGF therapy in glioblastoma. CNS Oncol.2, 49–65. 10.2217/cns.12.36
83
LyerS.BishopJ.AbramsB.MainoV.WardA.ChristianT.et al. (1997). QuantiBRITE: A New Standard for PE Flourescence Quantitation. San Jose, CA: Becton Dickinson Immunocytometry Systems.
84
Mac GabhannF.PopelA. S. (2007). Dimerization of VEGF receptors and implications for signal transduction: a computational study. Biophys. Chem.128, 125–139. 10.1016/j.bpc.2007.03.010
85
Mac GabhannF. and Popel, A. S. (2008). Systems biology of vascular endothelial growth factors. Microcirculation. 15, 715–738. 10.1080/10739680802095964
86
MakA. B.BlakelyK. M.WilliamsR. A.PenttiläP.ShukalyukA. I.OsmanK. T.et al. (2011). CD133 protein N-glycosylation processing contributes to cell surface recognition of the primitive cell marker AC133 epitope. J. Biol. Chem.286, 41046–41056. 10.1074/jbc.M111.261545
87
MamerS. B.ChenS.WeddellJ. C.PalaszA.WittenkellerA.KumarM.et al. (2017). Discovery of high-affinity PDGF-VEGFR interactions: redefining RTK dynamics. Sci. Rep.7:16439. 10.1038/s41598-017-16610-z
88
MarxM.PerlmutterR. A.MadriJ. A. (1994). Modulation of platelet-derived growth factor receptor expression in microvascular endothelial cells during in vitro angiogenesis. J. Clin. Invest.93, 131–139. 10.1172/JCI116936
89
MasseyS. C.RockneR. C.Hawkins-DaarudA.GallaherJ.AndersonA. R. A.CanollP.et al. (2018). Simulating PDGF-driven glioma growth and invasion in an anatomically accurate brain domain. Bull. Math. Biol.80, 1292–1309. 10.1007/s11538-017-0312-3
90
MazzoleniS.PolitiL. S.PalaM.CominelliM.FranzinA.SergiL. S.FaliniA.et al. (2010). Epidermal growth factor receptor expression identifies functionally and molecularly distinct tumor-initiating cells in human glioblastoma multiforme and is required for gliomagenesis. Cancer Res.70, 7500–7513. 10.1158/0008-5472.CAN-10-2353
91
McGranahanN.SwantonC. (2017). Clonal heterogeneity and tumor evolution: past, present, and the future. Cell168, 613–628. 10.1016/j.cell.2017.01.018
92
MoghadamS. A.AbadiA. M.MokhtariS. (2015). Immunohistochemical analysis of CD34 expression in salivary gland tumors. J. Oral Maxillofac. Pathol.19, 30–33. 10.4103/0973-029X.157197
93
MolinaE. S.PillatM. M.Moura-NetoV.LahT. T.UlrichH. (2014). Glioblastoma stem-like cells: approaches for isolation and characterization. J. Cancer Stem Cell Res.1:1. 10.14343/JCSCR.2014.2e1007
94
MuratA.MigliavaccaE.GorliaT.LambivW. L.ShayT.HamouM.TriboletN.et al. (2008). Stem cell–related ‘self-renewal' signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma. J. Clin. Oncol.26, 3015–3024. 10.1200/JCO.2007.15.7164
95
NaujokatC. (2014). Monoclonal antibodies against human cancer stem cells. Immunotherapy6, 290–308. 10.2217/imt.14.4
96
NguyenR.PerfettoS.MahnkeY. D.ChattopadhyayP.RoedererM. (2013). Quantifying spillover spreading for comparing instrument performance and aiding in multicolor panel design. Cytometry83, 306–15. 10.1002/cyto.a.22251
97
NielsenJ. S.McNagnyK. M. (2008). Novel functions of the CD34 family. J. Cell Sci.121, 3683–3692. 10.1242/jcs.037507
98
NiepelM.SpencerS. L.SorgerP. K. (2009). Non-genetic cell-to-cell variability and the consequences for pharmacology. Curr. Opin. Chem. Biol.13, 556–561. 10.1016/j.cbpa.2009.09.015
99
PalliniR.Ricci-VitianiL.MontanoN.MollinariC.BiffoniM.CenciT.et al. (2011). Expression of the stem cell marker CD133 in recurrent glioblastoma and its value for prognosis. Cancer117, 162–174. 10.1002/cncr.25581
100
PatelA. P.TiroshI.TrombettaJ. J.ShalekA. K.GillespieS. M.WakimotoH.CahillD. P.et al. (2014). Single-cell RNA-Seq highlights intratumoral heterogeneity in primary glioblastoma. Science344, 1396–1401. 10.1126/science.1254257
101
PatenaudeA.ParkerJ.KarsanA. (2010). Involvement of endothelial progenitor cells in tumor vascularization. Microvasc. Res.79, 217–223. 10.1016/j.mvr.2010.01.007
102
PavoineS.DolédecS. (2005). The apportionment of quadratic entropy: a useful alternative for partitioning diversity in ecological data. Environ. Ecol. Stat.12, 125–138. 10.1007/s10651-005-1037-2
103
PennockS.KimL. A.KazlauskasA. (2016). Vascular endothelial cell growth factor A acts via platelet-derived growth factor receptor α to promote viability of cells enduring hypoxia. Mol. Cell. Biol.36, 2314–2327. 10.1128/MCB.01019-15
104
PlateK. H.BreierG.FarrellC. L.RisauW. (1992). Platelet-derived growth factor receptor-beta is induced during tumor development and upregulated during tumor progression in endothelial cells in human gliomas. Lab. Invest.67, 529–534.
105
PlateK. H.BreierG.WeichH. A.MennelH. D.RisauW. (1994). Vascular endothelial growth factor and glioma angiogenesis: coordinate induction of VEGF receptors, distribution of vegf protein and possible in vivo regulatory mechanisms. Int. J. Cancer.59, 520–529. 10.1002/ijc.2910590415
106
PopescuA. M.AlexandruO.BrindusaC.PurcaruS. O.TacheD. E.TataranuL. G.et al. (2015). Targeting the VEGF and PDGF signaling pathway in glioblastoma treatment. Int. J. Clin. Exp. Pathol.8, 7825–7837.
107
PottsS. J.KruegerJ. S.LandisN. D.EberhardD. A.YoungG. D.SchmechelS. C.et al. (2012). Evaluating tumor heterogeneity in immunohistochemistry-stained breast cancer tissue. Lab. Invest.92, 1342–1357. 10.1038/labinvest.2012.91
108
QaziM. A.VoraP.VenugopalC.SidhuS. S.MoffatJ.SwantonC.et al. (2017). Intratumoral heterogeneity: pathways to treatment resistance and relapse in human glioblastoma. Ann. Oncol.28, 1448–1456. 10.1093/annonc/mdx169
109
RafteryA. E. (1995). Bayesian model selection in social research. Sociol. Methodol.25, 111. 10.2307/271063
110
RainesE. W.Bowen-PopeD. F.RossR. (1991). Platelet-derived growth factor, in Peptide Growth Factors and Their Receptors I, eds SpornM. B.RobertsA. B. (New York, NY: Springer), 173–262.
111
RaoC. R. (1982). Diversity and dissimilarity coefficients: a unified approach. Theor. Popul. Biol.21, 24–43.
112
RichJ. N.BignerD. D. (2004). Development of novel targeted therapies in the treatment of malignant glioma. Nat. Rev. Drug Discov.3, 430–446. 10.1038/nrd1380
113
Rocha-MartinsM.NjaineB.SilveiraM. S. (2012). Avoiding pitfalls of internal controls: validation of reference genes for analysis by QRT-PCR and western blot throughout rat retinal development. PLoS ONE7:e43028. 10.1371/journal.pone.0043028
114
RuchC.SkiniotisG.SteinmetzM. O.WalzT.Ballmer-HoferK. (2007). Structure of a VEGF–VEGF receptor complex determined by electron microscopy. Nat. Struct. Mol. Biol.14, 249–250. 10.1038/nsmb1202
115
SangarV.FunkC. C.KusebauchU.CampbellD. S.MoritzR. L.PriceN. D. (2014). Quantitative proteomic analysis reveals effects of epidermal growth factor receptor (EGFR) on invasion-promoting proteins secreted by glioblastoma cells. Mol. Cell. Proteomics13, 2618–2631. 10.1074/mcp.M114.040428
116
SarabipourS.HristovaK. (2016). Mechanism of FGF receptor dimerization and activation. Nat. Commun.7:10262. 10.1038/ncomms10262
117
SarabipourS.Ballmer-HoferK.HristovaK. (2016). VEGFR-2 conformational switch in response to ligand binding. Elife5:e13876. 10.7554/eLife.13876
118
SarkariaJ. N.YangL.GroganP. T.KitangeG. J.CarlsonB. L.SchroederM. A.GalanisE.et al. (2007). Identification of molecular characteristics correlated with glioblastoma sensitivity to EGFR Kinase inhibition through use of an intracranial xenograft test panel. Mol. Cancer Ther.6, 1167–1174. 10.1158/1535-7163.MCT-06-0691
119
SaundersN. A.SimpsonF.ThompsonE. W.HillM. M.Endo-MunozL.LeggattG.et al. (2012). Role of intratumoural heterogeneity in cancer drug resistance: molecular and clinical perspectives. EMBO Mol. Med.4, 675–684. 10.1002/emmm.201101131
120
SchonbergD. L.LubelskiD.MillerT. E.RichJ. N. (2014). Brain tumor stem cells: molecular characteristics and their impact on therapy. Mol. Aspects Med.39, 82–101. 10.1016/j.mam.2013.06.004
121
SeymourT.NowakA.KakulasF. (2015). Targeting aggressive cancer stem cells in glioblastoma. Front. Oncol.5:159. 10.3389/fonc.2015.00159
122
ShiQ.QinL.WeiW.GengF.FanR.ShinY. S.et al. (2012). Single-cell proteomic chip for profiling intracellular signaling pathways in single tumor cells. Proc. Natl. Acad. Sci. U.S.A.109, 419–424. 10.1073/pnas.1110865109
123
ShibuyaM. (2014). VEGF-VEGFR signals in health and disease. Biomol. Ther.22, 1–9. 10.4062/biomolther.2013.113
124
ShojaeiF.FerraraN. (2008). Role of the microenvironment in tumor growth and in refractoriness/resistance to anti-angiogenic therapies. Drug Resist. Updates11, 219–230. 10.1016/j.drup.2008.09.001
125
SimonsM.GordonE.Claesson-WelshL. (2016). Mechanisms and regulation of endothelial VEGF receptor signalling. Nat. Rev. Mol. Cell Biol.17, 611–625. 10.1038/nrm.2016.87
126
SimonsonA. B.SchnitzerJ. E. (2007). Vascular proteomic mapping in vivo. J. Thromb. Haemostasis5, 183–187. 10.1111/j.1538-7836.2007.02551.x
127
SinghS. K.HawkinsC.ClarkeI. D.SquireJ. A.BayaniJ.HideT.et al. (2004). Identification of human brain tumour initiating cells. Nature432, 396–401. 10.1038/nature03128
128
SmithC. L.KilicO.SchiapparelliP.Guerrero-CazaresH.KimD.Sedora-RomanN. I.et al. (2016). Migration phenotype of brain-cancer cells predicts patient outcomes. Cell Rep.15, 2616–2624. 10.1016/j.celrep.2016.05.042
129
SnuderlM.FazlollahiL.LeL. P.NittaM.ZhelyazkovaB. H.DavidsonC. J.AkhavanfardS.et al. (2011). Mosaic amplification of multiple receptor tyrosine kinase genes in glioblastoma. Cancer Cell20, 810–817. 10.1016/j.ccr.2011.11.005
130
SoaresA. B.JulianoP. B.AraujoV. C.MetzeK.AltemaniA. (2007). Angiogenic switch during tumor progression of carcinoma ex-pleomorphic adenoma. Virchows Arch.451, 65–71. 10.1007/s00428-007-0438-z
131
SodaY.MarumotoT.Friedmann-MorvinskiD.SodaM.LiuF.MichiueH.PastorinoS.et al. (2011). Transdifferentiation of glioblastoma cells into vascular endothelial cells. Proc. Natl. Acad. Sci. U.S.A.108, 4274–4280. 10.1073/pnas.1016030108
132
SoedaA.HaraA.KunisadaT.YoshimuraS.IwamaT.ParkD. M. (2015). The evidence of glioblastoma heterogeneity. Sci. Rep.5:7979. 10.1038/srep07979
133
SoedaA.InagakiA.OkaN.IkegameY.AokiH.YoshimuraS.et al. (2008). Epidermal growth factor plays a crucial role in mitogenic regulation of human brain tumor stem cells. J. Biol. Chem.283, 10958–10966. 10.1074/jbc.M704205200
134
SorkinA.DuexJ. E. (2010). Quantitative analysis of endocytosis and turnover of epidermal growth factor (EGF) and EGF receptor. Curr. Protoc. Cell Biol. Chapter 15:Unit 15.14. 10.1002/0471143030.cb1514s46
135
SpitzerM. H.NolanG. P. (2016). Mass cytometry: single cells, many features. Cell165, 780–791. 10.1016/j.cell.2016.04.019
136
SteinkampM. P.Low-NamS. T.YangS.LidkeK. A.LidkeD. S.WilsonB. S. (2014). ErbB3 is an active tyrosine kinase capable of homo- and heterointeractions. Mol. Cell. Biol.34, 965–977. 10.1128/MCB.01605-13
137
SzerlipN. J.PedrazaA.ChakravartyD.AzimM.McGuireJ.FangY.OzawaT.et al. (2012). Intratumoral heterogeneity of receptor tyrosine kinases EGFR and PDGFRA amplification in glioblastoma defines subpopulations with distinct growth factor response. Proc. Natl. Acad. Sci. U.S.A.109, 3041–3046. 10.1073/pnas.1114033109
138
SzopaW.BurleyT. A.Kramer-MarekG.KasperaW. (2017). Diagnostic and therapeutic biomarkers in glioblastoma: current status and future perspectives. Biomed. Res. Int.2017:8013575. 10.1155/2017/8013575
139
TaniguchiY.ChoiP. J.LiG.ChenH.BabuM.HearnJ.et al. (2010). Quantifying E. coli proteome and transcriptome with single-molecule sensitivity in single cells. Science329, 533–538. 10.1126/science.1188308
140
TiroshI.IzarB.PrakadanS. M.WadsworthM. H.TreacyD.TrombettaJ. J.RotemA.et al. (2016). Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science352, 189–196. 10.1126/science.aad0501
141
VerhaakR. G.HoadleyK. A.PurdomE.WangV.QiY.WilkersonM. D.MillerC. R.et al. (2010). Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell17, 98–110. 10.1016/j.ccr.2009.12.020
142
VigelsøA.DybboeR.HansenC. N.DelaF.HelgeJ. W.GrauA. G. (2015). GAPDH and β-actin protein decreases with aging, making stain-free technology a superior loading control in western blotting of human skeletal muscle. J. Appl. Physiol.118, 386–394. 10.1152/japplphysiol.00840.2014
143
VolinskyN.KholodenkoB. N. (2013). Complexity of receptor tyrosine kinase signal processing. Cold Spring Harb. Perspect. Biol.5:a009043. 10.1101/cshperspect.a009043
144
WeddellJ. C.ImoukhuedeP. I. (2014). Quantitative characterization of cellular membrane-receptor heterogeneity through statistical and computational modeling. PLoS ONE9:e97271. 10.1371/journal.pone.0097271
145
WeddellJ. C.ImoukhuedeP. I. (2017). Integrative meta-modeling identifies endocytic vesicles, late endosome and the nucleus as the cellular compartments primarily directing RTK signaling. Integr. Biol.9, 464–484. 10.1039/C7IB00011A
146
WeddellJ. C.ImoukhuedeP. I. (2018). Computational systems biology for the VEGF family in angiogenesis. Encycloped. Cardiovasc. Res. Med. 659–676. 10.1016/B978-0-12-809657-4.99548-6
147
WehlandM.BauerJ.MagnussonN. E.InfangerM.GrimmD. (2013). Biomarkers for anti-angiogenic therapy in cancer. Int. J. Mol. Sci.14, 9338–9364. 10.3390/ijms14059338
148
WeiW.ShinY. S.XueM.MatsutaniT.MasuiK.YangH.IkegamiS.et al. (2016). Single-cell phosphoproteomics resolves adaptive signaling dynamics and informs targeted combination therapy in glioblastoma. Cancer Cell29, 563–573. 10.1016/j.ccell.2016.03.012
149
WeickhardtA. J.WilliamsD. S.LeeC. K.ChionhF.SimesJ.MuroneC.WilsonK.et al. (2015). Vascular endothelial growth factor D expression is a potential biomarker of bevacizumab benefit in colorectal cancer. Br. J. Cancer113, 37–45. 10.1038/bjc.2015.209
150
WeisS. M.ChereshD. A. (2011). Tumor angiogenesis: molecular pathways and therapeutic targets. Nat. Med.17, 1359–1370. 10.1038/nm.2537
151
WernerS.HofschneiderP. H.HeldinC. H.OstmanA.RothW. K. (1990). Cultured kaposi's sarcoma-derived cells express functional PDGF A-type and B-type receptors. Exp. Cell Res.187, 98–103.
152
WheelerD. L.DunnE. F.HarariP. M. (2010). Understanding resistance to EGFR inhibitors[mdash]impact on future treatment strategies. Nat. Rev. Clin. Oncol.7, 493–507. 10.1038/nrclinonc.2010.97
153
WheelessL. L.CoonJ. S.CoxC.DeitchA. J.WhiteR. W. D.KossL. G.MelamedM.et al. (1989). Measurement variability in DNA flow cytometry of replicate samples. Cytometry10, 731–738. 10.1002/cyto.990100610
154
WuF. T. H.StefaniniM. O.GabhannF. M.KontosC. D.AnnexB. H.PopelA. S. (2010). A systems biology perspective on SVEGFR1: its biological function, pathogenic role and therapeutic use. J. Cell. Mol. Med.14, 528–552. 10.1111/j.1582-4934.2009.00941.x
155
YadavV. N.AltshulerD.KadiyalaP.ZamlerD.CombaA.AppelmanH.et al. (2018). Molecular ablation of tumor blood vessels inhibits therapeutic effects of radiation and bevacizumab. Neuro-Oncology, April. 10.1093/neuonc/noy055
156
ZoltánB. (2005). Rao's quadratic entropy as a measure of functional diversity based on multiple traits. J. Veg. Sci.16, 533–540. 10.1111/j.1654-1103.2005.tb02393.x
Summary
Keywords
single-cell, glioblastoma, RTK, heterogeneity, VEGFR, EGFR, IGFR, stem cell
Citation
Chen S, Le T, Harley BAC and Imoukhuede PI (2018) Characterizing Glioblastoma Heterogeneity via Single-Cell Receptor Quantification. Front. Bioeng. Biotechnol. 6:92. doi: 10.3389/fbioe.2018.00092
Received
13 March 2018
Accepted
21 June 2018
Published
11 July 2018
Volume
6 - 2018
Edited by
Paul De Vos, University Medical Center Groningen, Netherlands
Reviewed by
Alessandro Poggi, Dipartimento delle Terapie Oncologiche Integrate, Ospedale Policlinico San Martino, Italy; Cristiana Tanase, Victor Babes National Institute of Pathology, Romania
Updates
Copyright
© 2018 Chen, Le, Harley and Imoukhuede.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: P. I. Imoukhuede pimoukhuede@wustl.edu
This article was submitted to Biomaterials, a section of the journal Frontiers in Bioengineering and Biotechnology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.