<figure><img src="image_1.png" /></figure>

HAL
open science

# Comparative Expression Study of the Endo-G Protein Coupled Receptor (GPCR) Repertoire in Human Glioblastoma Cancer Stem-like Cells, U87-MG Cells and Non Malignant Cells of Neural Origin Unveils New Potential Therapeutic Targets

Marie Fève, Jean-Michel Saliou, Maria Zeniou, Sarah Lennon, Christine Carapito, Jihu Dong, Alain van Dorsselaer, Marie-Pierre Junier, Hervé Chneiweiss, Sarah Cianférani, et al.

► To cite this version:

Marie Fève, Jean-Michel Saliou, Maria Zeniou, Sarah Lennon, Christine Carapito, et al. Comparative Expression Study of the Endo-G Protein Coupled Receptor (GPCR) Repertoire in Human Glioblastoma Cancer Stem-like Cells, U87-MG Cells and Non Malignant Cells of Neural Origin Unveils New Potential Therapeutic Targets. PLoS ONE, 2014, 9 (3), pp.e91519. 10.1371/journal.pone.0091519. hal-01358283

HAL Id: hal-01358283

https://hal.sorbonne-universite.fr/hal-01358283v1

Submitted on 31 Aug 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Comparative Expression Study of the Endo–G Protein Coupled Receptor (GPCR) Repertoire in Human Glioblastoma Cancer Stem-like Cells, U87-MG Cells and Non Malignant Cells of Neural Origin Unveils New Potential Therapeutic Targets

Marie Fève¹,², Jean-Michel Saliou²,³, Maria Zeniou¹,³, Sarah Lennon², Christine Carapito², Jihu Dong¹, Alain Van Dorsselaer², Marie-Pierre Junier³, Hervé Chneiweiss³, Sarah Cianférani², Jacques Haiech¹, Marie-Claude Kilhoffer¹*

¹ Laboratoire d'Innovation Thérapeutique, UMR7200, Laboratoire d'Excellence Medalis, CNRS, Université de Strasbourg, Faculté de Pharmacie, Illkirch, France,

²Laboratoire de Spectrométrie de Masse BioOrganique, UMR7178, CNRS, Institut Pluridisciplinaire Hubert Curien, Université de Strasbourg, Strasbourg, France,

³ Neuroscience Paris Seine, UMR8246, Inserm U1130, Institut de Biologie Paris Seine, CNRS, Université Pierre et Marie Curie, Paris, France

## Abstract

Glioblastomas (GBMs) are highly aggressive, invasive brain tumors with bad prognosis and unmet medical need. These tumors are heterogeneous being constituted by a variety of cells in different states of differentiation. Among these, cells endowed with stem properties, tumor initiating/propagating properties and particularly resistant to chemo- and radiotherapies are designed as the real culprits for tumor maintenance and relapse after treatment. These cells, termed cancer stem-like cells, have been designed as prominent targets for new and more efficient cancer therapies. G-protein coupled receptors (GPCRs), a family of membrane receptors, play a prominent role in cell signaling, cell communication and crosstalk with the microenvironment. Their role in cancer has been highlighted but remains largely unexplored. Here, we report a descriptive study of the differential expression of the endo-GPCR repertoire in human glioblastoma cancer stem-like cells (GSCs), U-87 MG cells, human astrocytes and fetal neural stem cells (f-NSCs). The endo-GPCR transcriptome has been studied using Taqman Low Density Arrays. Of the 356 GPCRs investigated, 138 were retained for comparative studies between the different cell types. At the transcriptomic level, eight GPCRs were specifically expressed/overexpressed in GSCs. Seventeen GPCRs appeared specifically expressed in cells with stem properties (GSCs and f-NSCs). Results of GPCR expression at the protein level using mass spectrometry and proteomic analysis are also presented. The comparative GPCR expression study presented here gives clues for new pathways specifically used by GSCs and unveils novel potential therapeutic targets.

**Citation:** Fève M, Saliou J-M, Zeniou M, Lennon S, Carapito C, et al. (2014) Comparative Expression Study of the Endo–G Protein Coupled Receptor (GPCR) Repertoire in Human Glioblastoma Cancer Stem-like Cells, U87-MG Cells and Non Malignant Cells of Neural Origin Unveils New Potential Therapeutic Targets. PLoS ONE 9(3): e91519. doi:10.1371/journal.pone.0091519

**Editor:** Roland Seifert, Medical School of Hannover, Germany

**Received** December 28, 2013; **Accepted** February 10, 2014; **Published** March 24, 2014

**Copyright:** © 2014 Fève et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported by La Ligue contre le Cancer, La Fondation pour la Recherche Médicale, the Agence Nationale pour la Recherche (ANR), the Proteomic French Infrastructure (ANR-10-INSB-08-03), the Région Alsace, Université de Strasbourg, CNRS, the French Ministère de l'enseignement supérieur et de la recherche, the Fondation ARC (SFI20121205939, HC/MPJ) and has been preformed within the LABEX ANR-10-LABX-0034_Medalis and received a financial support from French government managed by "Agence National de la Recherche" under "Programme d'investissement d'avenir". Marie Fève is indebted to French Ministère de l'enseignement supérieur et de la recherche for doctoral fellowship and Sarah Lennon to the Fondation ARC pour la Recherche sur le Cancer. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: mckilhoffer@unistra.fr

ⓘ These authors contributed equally to this work.

## Introduction

Glioblastomas (World Health Organization (WHO) grade IV astrocytomas) are highly aggressive, angiogenic and infiltrating brain tumors representing more than 50% of all gliomas [1]. Their outcome is poor, most treatments currently in use remaining inefficient on long term survival [2,3,4] and less than 5% of patients survive 5 years post diagnosis [1]. Numerous clinical trials, mainly in phase I/II, are ongoing worldwide to find treatments with increased efficacy (http://clinicaltrials.gov/, http://apps.

who.int/trialsearch/and for recent review see [5]). Improvement of the standard Stupp protocol [3,6], association to topoisomerase inhibitors, alkylating agents, tyrosine kinase inhibitors, intercalating agents or antibodies targeting VEGF or the EGF receptor, development of new combinations of chemotherapeutic agents and testing new delivery drug devices have been used in different protocols (for exhaustive review see [5]). However no therapy appears as the panacea. Seeking for new molecules acting on targets different from those explored so far and improvement of
multi-targeted protocols are necessary to decrease recurrence and morbidity of this brain cancer that afflicts humans of all ages.

A new paradigm in cancer therapies arose a decade ago with the identification, within tumors, of cells endowed with stem cell properties and able to propagate tumors in immunodeficient mice with high efficacy. Such cells were first identified in hematological cancers [7]. In 2001, Reya postulated that they may be an integral part of most if not all tumors [8]. Cells with stem properties, long term survival, able to self-renew, differentiate into various cells types and lead to tumor formation after serial transplantation of few cells to immuno-deficient mice were found in many solid tumors including glioblastomas (GBM) [9,10], breast [11], colon cancer [12,13], melanoma [14,15] and pancreatic cancer [16]. Recent experiments, using genetic techniques, allowing the tracing of cancer cells within tumors strongly argue in favor of the presence of tumor stem cells within GBM [17], skin [18] and adenoma [19] tumors in vivo. These cancer stem cells (CSCs), also termed tumor stem cells, tumor initiating cells, tumor propagating cells or stem-like cancer propagating cells are considered as the cells responsible for tumor maintenance and relapse after treatment. Classical radiotherapies and chemotherapies which often induce tumor regression are thought to act preferentially on the non-stem cancer cells sparing CSCs and thus causing cancer recurrence. In addition, CSCs are thought to be able to enter a quasi-quiescent state compared to the rest of the tumor cells and thus escape treatments often designed to act on highly proliferating cells. CSCs have thus become the cells to target in order to increase cancer therapy efficacy. Numerous reviews are highlighting this point (for the most recent ones see [20,21,22,23,24,25,26]).

In the present study, we focused our attention on G-protein coupled receptors' (GPCRs) expression in cancer cells with stem and tumor propagating properties isolated from human GBM biopsies [27] and called herein glioblastoma stem-like cells (GSCs). The approach relied on the pivotal role played by this large family of membrane receptors in most key physiological functions including cell communication, cell signalling, cell migration, cell proliferation and apoptosis. These seven trans-membrane domain receptors convey information from outside the cell to inside by binding a variety of physiological compounds of wide chemical diversity. Their paramount importance in cell activity places GPCRs among the most appealing pharmacological targets. This has been extensively exploited to develop major drugs of therapeutic interest for a variety of pathologies including cardiovascular diseases, psychological disorders, allergies, so that nowadays 30-50% of current drugs target GPCRs [28,29,30,31]. Recently, GPCRs also gained attention in anti-tumor drug discovery due to the prominent role they can play, directly or indirectly though transactivation of other receptors, in cancer cell pathophysiology including proliferation, migration/metastasis, angiogenesis and metabolism [31,32,33,34,35,36]. The SMO antagonist vismodegib has recently been FDA-approved in basal cell carcinoma [37] and several molecules targeting GPCRs are in clinical trials as adjuvants of chemotherapies (namely Bosentan, an Endothelin A/B subtype receptor antagonist for melanoma or BKT140, a CXCR4 antagonist for multiple myeloma). In addition, the possible use in cancer therapy of monoclonal antibodies targeting GPCRs or their ligands has been reviewed in the literature [33,38].

In order to get a deeper insight into the pathophysiology of GSCs and find new potential therapeutic targets, we undertook a descriptive study of the set of GPCRs expressed in CSCs isolated from patient primary GBMs in comparison with the GPCR expression repertoires of the GBM cell line U87-MG, non-cancerous human fetal neural stem cells (f-NSCs) and primary

human astrocytes (HA). This set of cells was chosen as a first step to compare cancerous and non-cancerous neural cells with stem properties or in a more differentiated state. The study was performed both at the transcriptomic and proteomic levels. The expression profile of this family of proteins appears finely regulated despite karyotype alterations and changes in DNA ploidy [27]. The exquisitely tuned GPCR expression changes observed in the different cell types tested (consisting of cancerous and non-cancerous cells, stem and non-stem cells) lead to a molecular signature evidenced by hierarchical clustering. Finally, the study discloses several GPCRs with high expression specificity in GSCs compared to non-cancerous stem cells or the U-87 GBM cell line. In addition, GPCRs related to stemness have been highlighted. Those GPCRs may give a hook to modify the fate of cancer stem-like cells and more differentiated cancer cells.

## Results

### Global analysis of endo-GPCR gene expression in human GSCs, fetal NSCs, astrocytes in primary culture and GBM U-87 MG cells

Expression of GPCRs was studied in five different human cell types comprising TG1 and OB1 GSCs isolated from GBM biopsies of two patients, the GBM cancer cell line U-87 MG, fetal neural stem cells (f-NSCs) and human astrocytes (HA) in primary cultures. Expression studies were performed by RT-QPCR using the GPCR low density microarrays (Applied Biosystems) in the 384-well plate format. These microarrays allow testing 356 different GPCRs, 2 members of the LANCL family (LANCL1 and LANCL2), which are not considered as GPCRs, 14 different housekeeping genes (with 18 s ribosomal RNA present in quadruplicate) and 11 other genes not related to GPCRs. The list of genes present on the microarray is presented in Table S1. To analyze GPCR expression in the different cell types, data were normalized taking RPLP0, one of the most stable housekeeping genes in our experiments, as reference (Figure S1). The use of RPLP0 alone gave results similar to those using geometric averaging of multiple internal control genes [39]. The results were also similar to those obtained after normalization against TBP previously reported as being an adequate housekeeping gene for expression studies in GBM [40]. In the present experiments, the cycle threshold (Ct) value of RPLP0 was 21.1±0.28. RPLP0 is expressed at a ~6000 fold lower level compared to the ribosomal RNA 18 s. The most highly expressed GPCRs in each cell type are given in Table 1. GPR56 is highly expressed in all tumor cells but also and even at a higher level, in f-NSCs, the non-tumor stem cells of fetal origin. Its expression in astrocytes is reduced (Ct = 31.2; ratio over RPLP0 = 0.0009). For comparison, low levels of GPR56 were also observed in HEK293 cells, another non-stem cell line (Ct = 30.6 (after normalization); ratio over RPLP0 = 0.001) (data not shown).

Considering the whole set of 356 GPCRs (list in Table S1), we focused our attention on those expressed at Ct values ≤31.5 (i.e. expressed down to $1.33×10^{-3}$ fold the level of RPLP0) in at least one of the five cell types tested. This stringent condition gives the best signal to noise ratio, although it may leave aside some interesting genes expressed at lower levels. Out of the 356 different GPCRs present on the array, 138 satisfy this criterion (see list and expression values in Table S1). Expression of these 138 GPCRs in the different cell types is represented in Figure 1, where GPCRs are listed, in alphabetical order, as a function of their expression level (expressed as $2^{-Ct} \times 10^{12}$ (arbitrary units)). Among these 138 GPCRs, 90 genes were expressed above the cutoff level in OB1 cells, 83 genes in TG1 cells, 77 genes in f-NSCs, 61 genes in
<table><caption>**Table 1.** GPCRs most highly expressed in the five cell types.</caption><thead><tr><th>Cell type</th><th>GPCR</th><th>Ct</th><th>Delta Ct (Ct<sub>GPCR</sub> - Ct<sub>RPLP0</sub>)</th><th>GPCR expression compared to RPLP0 (fold)</th></tr></thead><tbody><tr><td rowspan="2">TG1</td><td>LPHN2</td><td>24±0.3</td><td>2.88±0.58</td><td>7.36±1.49</td></tr><tr><td>GPR56</td><td>24.2±0.04</td><td>3.08±0.32</td><td>8.46±1.25</td></tr><tr><td rowspan="3">OB1</td><td>F2R</td><td>24.3±0.01</td><td>3.18±0.28</td><td>9.06±1.22</td></tr><tr><td>GPR56</td><td>24.4±0.01</td><td>3.28±0.29</td><td>9.71±1.22</td></tr><tr><td>FZD7</td><td>24.5±0.03</td><td>3.38±0.31</td><td>10.41±1.24</td></tr><tr><td rowspan="3">f-NSC</td><td>LPHN2</td><td>24.6±0.04</td><td>3.48±0.32</td><td>11.16±1.25</td></tr><tr><td>GPR56</td><td>22.6±0.08</td><td>1.48±0.36</td><td>2.79±1.28</td></tr><tr><td>FZD3</td><td>24.7±0.18</td><td>3.58±0.46</td><td>12.38±1.38</td></tr><tr><td rowspan="2">HA</td><td>GPRC5B</td><td>24.8±0.15</td><td>3.68±0.43</td><td>13.09±1.35</td></tr><tr><td>FZD1</td><td>25±0.08</td><td>3.88±0.36</td><td>14.72±1.28</td></tr><tr><td rowspan="3">U87</td><td>FZD7</td><td>25±0.1</td><td>3.88±0.38</td><td>14.72±1.3</td></tr><tr><td>BDKRB2</td><td>24.5±0.3</td><td>3.38±0.58</td><td>10.41±1.49</td></tr><tr><td>GPR56</td><td>24.9±0.13</td><td>3.78±0.41</td><td>13.74±1.33</td></tr></tbody></table>

GPCRs with Ct values ≤25 are indicated. TG1 and OB1 correspond to the GBM cancer stem-like cells. f-NSC, HA and U-87 stand for fetal neural stem cells, human astrocytes and the U-87 GBM cell line, respectively. Ct corresponds to GPCR cycle threshold normalized to the Ct of RPLP0 (Ct = 21.12). GPCR fold expression corresponds to the expression of a given GPCR in a given cell type compared to the expression of RPLP0 calculated as: fold = 2^(-Ct<sub>RPLP0</sub> - Ct<sub>GPCR</sub>). doi:10.1371/journal.pone.0091519.t001

human astrocytes and 57 genes in the U-87 MG GBM cell line (see expression data in Table S1).

## Hierarchical clustering of GSCs, f-NSCs, human astrocytes and GBM U-87 MG cells using endo-GPCR gene expression

The 138 GPCRs significantly expressed (with Ct ≤31.5 or expression level >329 units) in at least one of the five cell types studied were used to perform unsupervised hierarchical clustering. Three main clusters could be observed (Figure 2). A first cluster gathers the GSCs TG1 and OB1, the second assembles human astrocytes (HA) and GBM U-87 MG cells and the third cluster corresponds to f-NSCs, which appeared closer to the tumor stem cells than to the astrocytes and GBM U-87 MG cells. The same three clusters are obtained when using all GPCRs present on the microarray without filtering (data not shown). The clustering clearly points to sets of genes specifically expressed or overexpressed either in a given cell type or common to two or more cell types. GPCR expression pattern thus constitutes a signature specific to a given cell type and allows clustering cells as a function of their physiological or pathophysiological states. Of note, the GBM U-87 MG cells cluster apart from the GSCs TG1 and OB1.

Comparative analysis of the differences in GPCR expression among the five cell types, with emphasis on GPCRs specifically expressed in the GBM CSCs will be performed in the following paragraphs.

## Analysis of the core of GPCRs expressed in all of the five cell types

Among the 138 GPCRs retained in our study, 26 are present at levels above the cutoff value (expression level ≥329 units, i.e. Ct ≤31.5) in the five cell types (Figure 3A and 3B). Receptors of the family of adhesion GPCRs, including BAI2, CD97, CELSR1, GPR125, LPHN1, LPHN2 and GPR56, are the most highly represented (7/26), followed by members of the Frizzled receptor family (5/26) comprising FZD1, FZD2, FZD5, FZD7 and FZD8. The vomeronasal GPCR VN1R1, is also a member of this group of “common GPCRs”. Although expressed in all five cell types,

most of the genes present a differential expression level between the five cell types.

As aforementioned, GPR56 is one of the most highly expressed genes in the three cancer cell types tested (TG1, OB1 and U-87 MG). Figure 3 shows that this gene is also highly expressed in f-NSCs. Its expression level in this cell type is respectively 3, 3.6 and 4.8 fold higher compared to that in TG1, OB1 and U-87 MG cells and largely exceeds its expression level in human astrocytes, which is close to the limit fixed in our study. Besides GPR56, the thrombin receptor F2R/PAR1 and to a lesser extent BAI2, GABBR1 and the orphan receptor GPR153 show significantly higher expression in the cells with stem cell properties (f-NSCs, TG1 and OB1), compared to the two other cell types. Only one gene in this group, LPHN2, appears to be overexpressed in the GSCs TG1 and OB1, but not in f-NSCs, nor in the GBM U-87 MG cells or in human astrocytes. Its expression in these two cell types exceeds the mean expression in the three other cell types HA, f-NSC and U-87 MG by 13 fold. On the overall, astrocytes and U-87 MG cells express the “commonly expressed” GPCR genes at the lowest level (with the exception of FZD1 and FZD7 in astrocytes).

## GPCRs with expression specificity in TG1 and OB1 GSCs

Expression data analyses unveil nine GPCRs with expression specificity in TG1 and OB1 cells (Figure 4). LPHN2, one of these GPCRs, belongs to the group of “commonly expressed” genes (Figure 3 and previous section). However, LPHN2 appears significantly overexpressed in the two GSC cell types compared to the GBM U-87 MG cells or to the non-cancerous cells (Figure 4).

Seven other genes are expressed with specificity in TG1 and OB1 cells. These include CALCRL, GPR37, GPR73, HRH2, GPR103 and GPR128, the levels of the latter two being close to the limit set in our study (Figure 4). The calcitonin receptor-like receptor CALCRL (or CRLR or CLR) is expressed above, but close to the cutoff level in astrocytes and f-NSCs, whereas it is downregulated in U-87 MG cells. Six genes, GPR37, HRH2, GPR73, EDG8, GPR128 and GPR103 are expressed above the cutoff level only in TG1 and OB1 cells. GPR37 is expressed well
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>
**Figure 1. GPCR expression in TG1, OB1, f-NSC, HA and U-87-MG cells.** GPCR expression was determined by RT-QPCR. GPCRs are presented in alphabetical order. Expression levels are related to the cycle threshold value and expressed as $2^{-Ct} * 10^{-12}$, where Ct corresponds to the cycle threshold of a given GPCR. A Ct of 31.5 corresponds to an expression level of 329 (arbitrary units). Dotted lines indicate the cutoff level (Ct = 31.5). Expression levels are presented on a logarithmic scale. All data were normalized using RPLP0 as housekeeping gene.
doi:10.1371/journal.pone.0091519.g001

above the cutoff value in TG1 and OB1 cells and is strongly upregulated in these cells in comparison to f-NSCs, HA and U-87 MG cells. GPR73 expression is increased in both TG1 and OB1 cells, with highest amplitude in TG1 cells. The last three receptors of the group, EDG8, GPR128 and GPR103 are specifically expressed in TG1 and OB1 cells, but at levels close to the cutoff value fixed in the present study. GPR103, although expressed in the lower range level, shows total specificity for TG1 and OB1 cells.

Among the 138 selected GPCRs, only one, FZD6, appears downregulated in TG1 and OB1 cells. FZD6 shows a rather high expression in HA cells and f-NSCs. It is downregulated in U-87 MG cells and absent in the two GSC types. FZD6 belongs to the frizzled family of GPCRs that counts a total of ten members (FZD1-10). Nine FZD genes (FZD1-9) have been analyzed in the present study (Figure 1). FZD10 mRNA could not be detected at a significant level in any of the five cell types tested. Among the nine remaining FZD genes, two may be noticed. FZD3 which shows the highest expression in the cells with stem properties and FZD6, the expression of which is almost completely turned down in the GSCs tested.

## GPCRs showing expression specificity in cells with stem cell properties

Besides the nine genes with specific changes in TG1 and OB1 cells compared to f-NSCs, astrocytes and U-87 MG cells (Figure 4), other genes show more specific expression in cells with stem properties. Ten genes show expression higher or in the same order of magnitude in f-NSC compared to TG1 and OB1 cells. These include BAI2, CXCR4, F2R, GABBR1 and GPRC5B (Figure 5A) and EDNRB (Figure 5B), BAI3, CHRM4, GPR19, and GPR82 (Figure 5C). Among these genes, CXCR4 appears to be the most specifically expressed in cells with stem properties. Its expression is similar in f-NSCs and TG1 cells, increased in OB1 cells, lower in HA cells and largely downregulated in U-87 MG cells (Figure 5A). Similarly to CXCR4, F2R/PAR1 shows the highest expression in OB1 cells. F2R belongs to the Proteinase activated receptors (PAR) also called Thrombin receptors or Coagulation factor II receptors. Two out of the four PARs, namely FR2 (PAR-1) and F2RL1 (PAR-2) are significantly expressed in four (PAR-2) and all (PAR-1) of the five cell types tested, respectively (Figure 1C). Both F2R/PAR-1 and F2RL1/PAR-2 show overexpression in TG1 and OB1 cells, but F2R/PAR-1 has also increased expression in f-NSCs, whereas F2RL1 (PAR-2) has higher expression in HA cells.

Three genes, CHRM3, GPR171 and P2Y5, although expressed in f-NSCs, show higher expression in TG1 and OB1 cells (Figure 5C).

Three more genes LGR4, CNR1 and GPR51 (GABBR2), although expressed in TG1 cells, appear more specific to f-NSCs and OB1 cells (Figure 5B and C).

Finally, CELSR2 expression is also upregulated in these cells, with the highest expression in OB1 cells (Figure 5B).

## GPCRs with higher expression in the cancer cells TG1, OB1 and U-87 MG compared to HA and f-NSCs

We previously showed that GPR56 was overexpressed in TG1, OB1, U-87 MG and f-NSC cells (Figure 2). Further analysis of the GPCR transcriptome indicates that the expression of one GPCR,

CD97, appears more specifically related to the cancer phenotype. Among these cancer cells, CD97 expression appears upregulated in cancer cells with stem cell properties compared to U-87 MG cells.

## GPCRs characteristic of a given cell type

Analysis of Figure 1 also points to GPCRs that appear highly expressed in only one of the five cell types (Figure 6). Thus, HA cells are characterized by the specific expression of OXTR and DRD4, f-NSCs by BAI1, GPR37L1 and U-87 MG cells by BDKRB1 and to some extent BDKRB2 (although this GPCR is also expressed at low levels in HA cells).

## Proteomic analysis of GPCR expression in human GBM CSCs, f-NSCs, HA and GBM U-87 MG cells

One step further in our investigation was to study GPCR expression at the protein level. Antibodies are often used as tools to investigate protein expression and localization. However the use of antibodies is cumbersome if a great number of proteins have to be analyzed. Moreover, the use of antibodies and especially those targeting GPCRs is often hampered by their lack of specificity. This has been reported for several receptors [41] including α1 adrenergic [42], dopamine [43], beta adrenergic [44], muscarinic receptors [45] and histamine H4 receptors [46]. To overcome the tedious steps needed for generation and validation of reliable antibodies, other approaches are needed. In recent years, high throughput protein profiling, even in complex biological systems, has gained considerable interest thanks to the development of mass spectrometry techniques (MS) coupled to liquid chromatography and bioinformatics.

Thus, in parallel to transcriptomic analysis, we developed a mass spectrometry-based proteomic approach to investigate the expression of GPCRs at the protein level. This task appeared challenging as GPCRs are membrane proteins mostly expressed at low levels. Therefore, a dedicated membrane proteomic approach based on 1D gel protein separation, followed by reduction, alkylation and enzymatic digestion prior to nanoLC-MS/MS analysis of peptide mixtures was developed (see material and methods). With a 1% false discovery rate (FDR), between 3000 and 5000 unique proteins were identified for each cell type (more precisely 4090 for f-NSCs, 3997 for TG1, 3960 for OB1, 3758 for U-87 MG and 3216 for HA cells). This led to more than 7000 non-redundant protein identifications when the results from the five different cell types were combined. Among those proteins, about 20% were annotated as plasma membrane proteins according to their GO annotation: 1043 for f-NSCs, 976 for OB1, 966 for U-87 MG, 963 for TG1 and 837 for HA cells. A final list of seven GPCRs was obtained, among which five (GPR56, LPHN2, F2R/PAR1, GPRC5B and CD97) belong to the fifteen most expressed transcripts (expression level >20000 units) (Table 2). Interestingly, CELSR2 which is expressed at a lower mRNA level (7542 units maximum in OB1 cells) was identified with at least 2 peptides in f-NSCs, TG1 and OB1 cells. A detailed analysis of proteomic results was performed in order to explain why some highly abundant gene transcripts were not detected at the protein level, whereas less abundant ones were. Based on in depth GPCR amino acid sequence analysis, the number of MS detectable peptides was estimated (Table 2). To be
<figure><img src="image_12.png" /><figcaption>**Figure 2. Hierarchical clustering of TG1, OB1, f-NSC, HA and U-87 MG cells based on GPCR expression.** The set of 138 GPCRs with Ct values ≤31.5 (or expression levels expressed as 2⁻⁵CT * 10⁻¹²≥329 units) was used for the clustering. For a given cell type, mean expression values (m) from at least two experiments were used for the clustering. Clustering was performed using the DChip software [180]. doi:10.1371/journal.pone.0091519.g002</figcaption></figure>

best detected in MS, peptides should have: i) molecular weights comprised between 1000 and 3000 Da; ii) extracellular or cytoplasmic regions exposed to trypsin digestion and characterized by lower hydrophobicity than trans-membrane domains and iii) no posttranslational modifications. Considering those criteria, our analyses suggest that, grossly, GPCRs with more than 5 potentially detectable peptides and with gene expression levels over 20000 units were observed in the present proteomic study.

To obtain semi-quantitative information on proteins identified by mass spectrometry, we performed spectral and peptide counting analysis. For GPR56, one of the best identified GPCRs by proteomic analysis, spectral counting led to an expression pattern similar to that observed at the transcript level, i.e. expression of GPR56 being higher in f-NSCs than in TG1, OB1 and U-87 MG cells and not detected in HA (Table 2). LPHN2 which shows high transcript levels specifically in OB1 and TG1 cells, exhibits similar results at the protein level (LPHN2 is detected with a high number of peptides and spectra in both TG1 and OB1 cells and is not detected in f-NSCs, U-87 MG cells or in astrocytes) (Table 2). Interestingly, all the MS/MS selected peptides are located in the extracellular domain of LPHN2 even if 34 peptides could have been expected from the cytoplasmic domain. This observation could be consistent with cytoplasmic domain processing which eliminates part of the C-terminal domain of the receptor. CD97 is another GPCR for which proteomics follows the trend of transcriptomics. CD97 transcript was found only in the cancer cells with overexpression in GSCs. At the protein level CD97 appears highly overexpressed in TG1 (9 peptides), OB1 (10 peptides) and U-87 cells (6 peptides) while it was not identified in f-NSCs and HA cells (Table 2). Also, CELSR2 which showed a transcript level in OB1>f-NSC>TG1>HA was detected with the same expression hierarchy at the proteomic level with OB1>f-NSC>TG1>HA (Table 2). F2R (PAR-1), GPRC5B and BAI2 belong to the group of GPCRs overexpressed at the mRNA level in cells with stem properties. With only five potentially detectable peptides, F2R was identified by proteomic analysis only in OB1 cells that present the highest mRNA level for this GPCR (Table 2). Despite a lower level of transcript and few peptides detectable by MS/MS, GPCR5B was identified in f-NSCs, OB1 and TG1 cells in agreement with the RT-QPCR analysis (Table 2). Surprisingly, BAI2, a 173 kDa protein for which 40 unique peptides could be expected by MS analysis and with a mRNA level similar to that of CD97 (Table 2), was identified by a single peptide in f-NSCs (the cell type which showed the highest transcript level) (Table 2).

## Discussion

This paper presents the first exhaustive study of GPCR expression at both the transcriptomic and proteomic levels in GSCs. Analysis of the difference in GPCR expression between several cell types of neural origin highlights potential functioning specificities. Investigating the GPCR proteome was challenging as most of these membrane receptors are expressed at low levels. Besides the expression level, receptor accessibility to proteolytic enzymes, which can be related to the length of the extra- and intracellular domains and post-transcriptional modifications, must
<figure><img src="image_13.png" /><figcaption>**Figure 3. Ensemble of GPCRs expressed in all of the five cell types tested.** GPCRs expressed with a cycle threshold (Ct) ≤ to 31.5 in GSCs (TG1 and OB1) human fetal neural stem cells (f-NSCs), human astrocytes (HA) and U-87 MG GBM cells (U-87) are represented. Ct = 31.5 (or expression level expressed as 2⁻⁻Ct * 10⁻¹² = 329 units) was taken as the lowest limit for expression. All data were normalized using RPLP0 as housekeeping gene. Ordinate indicates the fold change in expression of a given GPCR relative to the Ct value of 31.5. Relative fold change is calculated as 2⁽³¹.₅-Ct_GPCR⁾, where Ct_GPCR corresponds to the cycle threshold of a given GPCR in a given cell type. Values correspond to the mean of at least 2 experiments (from 2 to 4). GPCRs are presented in alphabetical order. Note the ordinate scale change for GPR 56.
doi:10.1371/journal.pone.0091519.g003</figcaption></figure>

be considered. Seven GPCRs were detected by MS - based proteomics. For all GPCRs, but BAI2, detection by mass spectrometry correlates with their transcript level and number of potential proteolytic peptides. These results thus demonstrate that MS - based proteomic approaches nowadays allow both identification of proteins in complex mixtures of thousands of proteins and quantification with good confidence. It also reinforces transcriptomic data by showing (for the receptors that could be detected) a good correlation between levels evaluated by RT Q-PCR and mass spectrometry, i.e. the highest protein levels being observed for the highest mRNA levels. Due to the lack of selectivity of most commercial GPCR - directed antibodies, MS appears as an attractive methodology for GPCR expression profiling.

The different GPCRs presented in the result section, with expression specificity in either GBM cancer cells, GBM cancer cells and f-NSCs, GSCs and GSCs/f-NSC) are presented in Table 3 along with their known ligand(s), class, function and chromosome localization. In the discussion we will focus our attention on the most highly expressed genes detected by proteomics and on some GPCRs with specific expression in TG1 and OB1 cells. The relevant GPCRs will be discussed in order to highlight their main known functions and, if documented, their implications in glioma or, more generally, in cancer.

## Genes with specific overexpression in GBM cancer cells

One gene, CD97, emerges for its specificity in GBM cancer cells. This GPCR shows significant expression change in the brain cancer cells vs non-cancer cells and its expression is highest in cancer cells with stem cell properties (TG1 and OB1). This expression specificity is observed both at the mRNA and protein levels. CD97 is a group II adhesion GPCR, member of the

epidermal growth factor-seven transmembrane (EGF-TM7) family. Aberrant expression of this gene has been reported in various cancers, namely thyroid, gastric, pancreatic, and esophageal carcinomas [47] and expression seems to be related to tumor aggressiveness. Recently, CD97 was shown to be upregulated in three different GBM cell lines and to be one of the targets of the transcription factor WT-1 (Wilms tumor protein), the expression of which is also upregulated in gliomas [48]. CD97 was found to be involved in invasiveness and metastasis of prostate cancer, (where it forms heterodimers with the lysophosphatidic acid receptor LPAR1 [49,50]) and more recently of GBMs [48,51]. Its overexpression is correlated with decreased survival in patients [51]. The present study clearly shows the expression specificity of this receptor in GBM cancer cells compared to non-cancerous neural cells. Further studies need to be performed to establish which of the possible splice variants [52] are upregulated. So far there are no reports on the physiological role of CD97 in brain. The low expression of this receptor in astrocytes and f-NSCs compared to cancer cells places it as a good candidate for further investigation of its role in cancer in general and GSC cells in particular.

## Genes expressed in GBM cancer cells with and without stem properties and in f-NSCs

GPR56, another adhesion GPCR belonging to class VIII is expressed in GBM cells, but its expression is even higher in f-NSCs. GPR56 is known for its expression in cancer cells, including gliomas [53]. In several tumors (melanoma, prostate cancer) [54,55,56], the level of GPR56 expression appears to be inversely correlated to malignancy and metastatic potential. Physiologically, GPR56 plays a key role in brain development, being involved in cortical development and radial migration of neural progenitors
<figure><img src="image_14.png" /><figcaption>**Figure 4.** GPCRs preferentially expressed in TG1 and OB1 GSCs.
GPCR expression was determined by RT-QPCR. Expression levels are presented as 2^(ΔCt) * 10^(ΔΔCt), where Ct corresponds to the cycle threshold of a given GPCR. A) and B) present GPCRs overexpressed in TG1 and OB1 cells and GPCRs are presented as a function of decreasing expression level in TG1 cells; C) presents the GPCR specifically underexpressed in TG1 and OB1 GSCs. Data correspond to the mean of 4 experiments.
doi:10.1371/journal.pone.0091519.g004</figcaption></figure>

(NPC/NSC) from the lateral ventricles to the cortical plate. Collagen has been reported to be its major ligand for this function [57]. GPR56 has been shown to be highly expressed in nestin-positive neural stem/progenitor cells in the ventricular/subventricular zone of human and mouse fetal brains as well as in cultured neurospheres derived from both human and mouse fetal brains [58] [59]. Its expression is downregulated during differentiation [59].

The high expression of GPR56 in f-NSCs and its low expression in astrocytes are in agreement with the pre-cited studies performed on neural stem/progenitor cells. However the present study which compares the expression of GPR56 in GBM cells and f-NSCs

shows for the first time that GPR56 can be expressed at the same and even higher level in normal (non-cancerous) stem cells than in cancer cells.

## GPCRs showing expression specificity in cells with stem cells properties TG1, OB1 and f-NSC cells

At the transcriptomic level, 17 GPCRs show specific expression in cells with stem properties (Figure 5A–C; Table 3). Four of them, BAI2, F2R, GPCR5B and CELSR2 are detected at the protein level. Information relative to the role of these receptors in stem cells is limited and their role in CSCs remains to be explored.

BAI2 is a member of the brain angiogenesis inhibitor family (BAIs) which belongs to group VII of adhesion GPCRs and counts 3 members, BAI1 to 3, located on chromosomes 8, 1 and 6 respectively. BAIs are mainly expressed in brain. This group exhibits an N-terminal CUB (for complement C1r/C1s, Uegf, Bmp1) domain, 3 to 5 thrombospondin type 1 repeats (TSR) discovered to regulate the anti-angiogenic activity of thrombospondin-1, one hormone-binding domain, the evolutionary conserved GPCR proteolytic site (GPS) and multiple glycosylation sites. The three receptors are involved in neo-vascularization and have been reported to be involved in apoptosis and regulation of tumor progression (for review see [60]). The antiangiogenic activity of BAI1 in normal and pathological states is well documented [61,62,63]. It is triggered by vasculostatin-40, a 40 kDa fragment induced by matrix metalloprotease cleavage [63]. Tumor-suppressor activity has thus been attributed to BAI1 and its use in cancer gene therapy has been proposed [62]. The angiostatic effect of BAI2 (and BAI3) is less documented. Activation through N-terminal cleavage has been reported for BAI2 [64]. Concerning a role of this receptor family in high grade gliomas, data are scarce and concern global levels of expression measured on tissue samples composed of a mixture of cells [65]. The present study which gives a comparative overview of the transcription expression pattern in five different cell types shows that f-NSCs express BAI1-3 at the highest level (Figure 1). BAI2-3 are overexpressed in cells with stem properties compared to more differentiated cells. Such observations have not yet been reported in the literature. Only the lack of expression of BAI1 in various GBM cell lines, including the U-87 MG cells has been documented [61], which is in agreement with the results of the present study (Figure 1). The different members of the family and namely BAI2 and BAI3 may be involved in functions different from neovascularization, related to stemness maintenance or regulation of apoptosis. Further experiments are necessary to uncover the role of these receptors in GBM and in cancer, more generally. Of the three BAIs, the transcription level of BAI2 is highest. Despite its high molecular weight and high level of mRNA in f-NSCs, TG1 and OB1 cells, BAI2 is only scarcely detected (one peptide) at the proteomic level in f-NSCs (Table 2). This raises further questions concerning the expression of the gene at the protein level or receptor processing in the stem-like cancer cells.

F2R is a second GPCR detected by mass spectrometry with expression specificity in cells with stem properties. F2R belongs to the Proteinase activated receptors (PAR) also called Thrombin receptors or Coagulation factor II receptors which count four members. Two PAR receptors, namely FR2 (PAR-1) and F2RL1 (PAR-2) are significantly expressed in, respectively, four (PAR-2) and all (PAR-1) of the five cell types tested (Figure 1). Both F2R/ PAR-1 and F2RL1/PAR2 show overexpression in TG1 and OB1 cells. F2R, but not F2RL1 is also overexpressed in f-NSCs. PAR receptors have plural effects on neural cells. PAR receptors and namely F2R/PAR-1 and F2RL1/PAR-2 are expressed in most tumors (for review see [66,67]) and are known for their
<figure><img src="image_15.png" /><figcaption>**Figure 5.** GPCRs overexpressed in the GSCs TG1 and OB1 and the fetal neural stem cells (f-NSCs). The ordinate represents the expression level of the different GPCRs shown as (2⁻⁻⁰⁻⁰⁻¹²) where Ct corresponds to the cycle threshold of a given GPCR. Note the difference in scale between A, B and C. HA corresponds to human astrocytes and U-87 to the GBM cell line U-87 MG.
doi:10.1371/journal.pone.0091519.g005</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 6.** GPCRs specifically expressed in one of the five cell types. Expression of the given GPCR is compared between TG1 and OB1 GSCs, fetal neural stem cells (f-NSCs), human astrocytes (HA) and the GBM cell line U-87 MG. The ordinate represents the expression level of the different GPCRs shown as (2⁻⁻⁰⁻⁰⁻¹²) where Ct corresponds to the cycle threshold of a given GPCR. Note the difference in scale between A) and B). TG1 corresponds to GSCs, HA to human astrocytes and U-87 to the GBM cell line U-87 MG.
doi:10.1371/journal.pone.0091519.g006</figcaption></figure>

cooperation with EGFR signaling [68]. F2RL1/PAR-2 triggers proliferation, migration, invasion of GBM cells and activation of downstream MAPKs including ERKs after stimulation by the HIF-inducible tissue factor (TF) [69,70]. F2R/PAR-1, the principal thrombin-activated receptor involved in platelet aggregation and endothelial cell proliferation, has been shown to be overexpressed in various invasive and metastatic tumors and its expression levels directly correlate with the degree of invasiveness of the cancer. A link between angiogenesis and EGFR (via the ligand independent EGFRvIII mutant), tissue factor (TF) and PAR receptors has recently been established in GBM cell lines [71]. In melanoma, cell signaling through PAR releases several factors, such as cytokines, proteases, adhesion molecules and growth factors causing changes in the tumor microenvironment [72]. This has not been shown so far for GBMs. PAR receptors probably play a major role in GBM, and even more generally in cancer, linking the tumor microenvironment, hypoxia and angiogenesis to tumor cell proliferation, migration and tissue invasion. However, PARs also play a major role in normal cell functioning and are expressed, as shown in the present study, in astrocytes and fetal neural stem cells. The overexpression of these receptors in cancer stem-like cells has not been reported so far. F2R is detected by mass spectrometry in OB1 cells that have the highest mRNA level for this receptor (Table 2).
<table><caption>**Table 2.** GPCRs detected by proteomic analysis and comparison to those most highly expressed at the transcriptomic level.</caption><thead><tr><th rowspan="2">GPCR Name</th><th rowspan="2">Access number</th><th rowspan="2">Molecular weight (kDa)</th><th colspan="2">OB1</th><th colspan="2">TG1</th><th colspan="2">f-NSC</th><th colspan="2">HA</th><th colspan="2">U87</th></tr><tr><th>Nb of expected unique peptides</th><th>Transcript level (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th><th>Nb of unique Transcript peptide (nb level of spectrum) (2 c<sub>3</sub>)*10<sup>12</sup></th></tr></thead><tbody><tr><td>GPR56</td><td>Q91653</td><td>78</td><td>14</td><td>42820</td><td>4 (8)</td><td>52347</td><td>4 (10)</td><td>155121</td><td>8 (37)</td><td>625</td><td>nd</td><td>32369</td><td>3 (6)</td></tr><tr><td>LPHN2</td><td>O95490</td><td>163</td><td>50</td><td>37915</td><td>9 (15)</td><td>62684</td><td>12 (24)</td><td>2338</td><td>nd</td><td>5921</td><td>nd</td><td>3495</td><td>nd</td></tr><tr><td>F2R</td><td>P25116</td><td>44</td><td>5</td><td>48045</td><td>2 (6)</td><td>26751</td><td>nd</td><td>18912</td><td>nd</td><td>7701</td><td>nd</td><td>1259</td><td>nd</td></tr><tr><td>BDKRB2</td><td>P30411</td><td>44</td><td>0</td><td>78</td><td>nd</td><td>0</td><td>nd</td><td>1</td><td>nd</td><td>1655</td><td>nd</td><td>42664</td><td>nd</td></tr><tr><td>FZD7</td><td>O75084</td><td>61</td><td>3</td><td>39613</td><td>nd</td><td>10182</td><td>nd</td><td>14788</td><td>nd</td><td>30943</td><td>nd</td><td>4011</td><td>nd</td></tr><tr><td>GPRC5B</td><td>Q9NZH0</td><td>41</td><td>3</td><td>24853</td><td>1 (3)</td><td>24463</td><td>1 (4)</td><td>33650</td><td>3 (7)</td><td>2888</td><td>nd</td><td>248</td><td>nd</td></tr><tr><td>TM7SF1</td><td>O60478</td><td>45</td><td>3</td><td>15303</td><td>nd</td><td>3322</td><td>nd</td><td>29968</td><td>nd</td><td>2855</td><td>nd</td><td>1743</td><td>nd</td></tr><tr><td>CXCR4</td><td>P61073</td><td>40</td><td>2</td><td>29644</td><td>nd</td><td>14193</td><td>nd</td><td>11354</td><td>nd</td><td>120</td><td>nd</td><td>19</td><td>nd</td></tr><tr><td>BA12</td><td>O60241</td><td>173</td><td>40</td><td>19770</td><td>nd</td><td>15863</td><td>nd</td><td>29243</td><td>1 (1)</td><td>5050</td><td>nd</td><td>5478</td><td>nd</td></tr><tr><td>BDKRB1</td><td>P46663</td><td>40</td><td>2</td><td>5</td><td>nd</td><td>0</td><td>nd</td><td>9</td><td>nd</td><td>65</td><td>nd</td><td>27717</td><td>nd</td></tr><tr><td>GPR153</td><td>Q6NV75</td><td>65</td><td>1</td><td>8111</td><td>nd</td><td>26080</td><td>nd</td><td>7762</td><td>nd</td><td>2438</td><td>nd</td><td>342</td><td>nd</td></tr><tr><td>CD97</td><td>P48960</td><td>90</td><td>19</td><td>14593</td><td>10 (40)</td><td>24958</td><td>9 (31)</td><td>405</td><td>nd</td><td>1079</td><td>nd</td><td>5233</td><td>6 (21)</td></tr><tr><td>LGR4</td><td>Q8N537</td><td>104</td><td>8</td><td>22726</td><td>nd</td><td>4135</td><td>nd</td><td>19526</td><td>nd</td><td>1804</td><td>nd</td><td>51</td><td>nd</td></tr><tr><td>CNR1</td><td>P21554</td><td>53</td><td>3</td><td>20381</td><td>nd</td><td>3161</td><td>nd</td><td>14232</td><td>nd</td><td>190</td><td>nd</td><td>1</td><td>nd</td></tr><tr><td>CELSR2</td><td>O9HCU4</td><td>317</td><td>79</td><td>7543</td><td>2 (2)</td><td>1769</td><td>2 (2)</td><td>3501</td><td>8 (12)</td><td>601</td><td>nd</td><td>118</td><td>nd</td></tr></tbody></table>

List of the GPCRs with the highest transcription level in at least one of the five cell types (expressed as 2<sup>-ΔΔCt</sup>*10<sup>12</sup> values >20000) and their identification by proteomic analysis. GPCRs are listed from top to bottom as a function of their mRNA level (the most highly expressed in one or several cell types appears at the top). Number of expected unique peptides corresponds to the unique peptides potentially detectable by MS/MS analysis. Expected peptides were calculated after *in silico* trypsin digestion of a given GPCR taking into account different criteria: i) no miscleavage, ii) molecular weights comprised between 1000 and 3000 Da, iii) presence in extracellular or cytoplasmic regions and iv) no posttranslational modifications. GPCRs identified by proteomic analysis are indicated in italics. OB1 and TG1 correspond to the GSCs, f-NSC, HA and U-87 stand for fetal neural stem cells, human astrocytes and the U-87 MG GBM cell line, respectively. Ct=cycle threshold; nd=not detected).
doi:10.1371/journal.pone.0091519.t002
<table><caption>**Table 3.** Summary of the GPCRs expressed with the highest specificity in different subgroups of cells comprising at least the GSCs TG1 and OB1.</caption><thead><tr><th>Subgroup</th><th>Function</th><th>Ligand</th><th>Class</th><th>Chromosomal localisation</th></tr></thead><tbody><tr><td colspan="5"><strong>TOU</strong></td></tr><tr><td>CD97</td><td>tumor angogenesis, cell migration, EMT</td><td>CD55</td><td>B Adhesion/EGF-like family</td><td>19p13</td></tr><tr><td colspan="5"><strong>NTOU</strong></td></tr><tr><td>GPR56</td><td>brain development; tumor metastasis; cell plasticity</td><td>collagen III</td><td>B Adhesion</td><td>16q13</td></tr><tr><td colspan="5"><strong>TO</strong></td></tr><tr><td>CALCRL/CRLR</td><td>angiogenesis, cell proliferation, antiapoptotic, resistance to hypoxia</td><td>adrenomedullin, CGRP</td><td>B Calcitonin group</td><td>2q32.1</td></tr><tr><td>EDG8/S1P5</td><td>cell migration/proliferation</td><td>lysophosphatidic acid, sphingosine 1-phosphate</td><td>A</td><td>19p13.2</td></tr><tr><td>FZD6</td><td>tumor suppressor</td><td>Wnt</td><td>B Frizzled</td><td>8q22.3-q23.2</td></tr><tr><td>GPR103</td><td>pain. RF-amide Receptor</td><td>QRFP family peptides</td><td>A</td><td>4q27</td></tr><tr><td>GPR128</td><td>orphan no report</td><td>orphan</td><td>B Adhesion</td><td>3q12.2</td></tr><tr><td>GPR37</td><td>macroautophagy, glio and neuronal protection</td><td>prosaposin/Prosapeptide?</td><td>A</td><td>7q31</td></tr><tr><td>GPR73/PKR1</td><td>angiogenesis, cell proliferation, invloved in pain such as GPR103</td><td>prokinecitin-1 and 2</td><td>A</td><td>2p13.1</td></tr><tr><td>HRH2</td><td>angiogenesis, cell proliferation</td><td>histamine</td><td>A</td><td>5q35.2</td></tr><tr><td>LPHN2</td><td>EMT in heart development</td><td>lectomedin</td><td>B Adhesion</td><td>1p31.1</td></tr><tr><td colspan="5"><strong>NTO</strong></td></tr><tr><td>BAI2</td><td>unclear (neovascularization, apotosis)</td><td>orphan</td><td>B Adhesion</td><td>1p25.2</td></tr><tr><td>BAI3</td><td>synapse formation, denditric morphogenesis</td><td>C1q-like proteins (C1q1-4) ?</td><td>B Adhesion</td><td>6q12</td></tr><tr><td>CELSR2</td><td>planar cell plority</td><td>Cadherin domains - homophilic binding</td><td>B Adhesion</td><td>11p13.3</td></tr><tr><td>CHRM3</td><td>Cell proliferation,differentiation, survival. Link with EGFR activation. Implicated in cancer</td><td>acetylcholine</td><td>A</td><td>1q43</td></tr><tr><td>CHRM4</td><td>Cell proliferation, differentiation, survival. Implicated in cancer</td><td>acetylcholine</td><td>A</td><td>11p11.2</td></tr><tr><td>CNR1</td><td>energy metabolim; glucose oxidation, ketogenesis, apoptosis, autophagy; stem gene expression</td><td>anandamide, 2-arachidonyl-glycerol</td><td>A</td><td>6q14-q15</td></tr><tr><td>CXCR4</td><td>cell migration, angiogenesis, invasivness hypoxia controlled, member of the CRG (cooperation regulatory genes).</td><td>CXCL12</td><td>A</td><td>2q22.1</td></tr><tr><td>EDNRB</td><td>self renewal, migration, proapoptotic</td><td>endothelins 1 to 3</td><td>A</td><td>13q22.3</td></tr><tr><td>F2R/PAR-1</td><td>Various pathophysiological processes. Participates in EGFR signalling and thus involved in carcinoma.HIF inducible.</td><td>proteolytic cleavage of N-terminus</td><td>A</td><td>5q13.3</td></tr><tr><td>GABBR1</td><td>cell proliferation</td><td>GABA</td><td>C Metabotropic GABA(B)R</td><td>6p22.1</td></tr><tr><td>GABBR2 (GPR51)</td><td>cell proliferation</td><td>GABA</td><td>C Metabotropic GABA(B)R</td><td>9q22.1-22.3</td></tr><tr><td>GPR19</td><td>links with stemness</td><td>orphan</td><td>A</td><td>12p12.3</td></tr><tr><td>GPR82</td><td>orphan, food intake, not essential</td><td>orphan</td><td>A</td><td>Xp11.4</td></tr><tr><td>GPR171</td><td>Unknown; retinoic acid induced</td><td>orphan</td><td>A</td><td>3q25.1</td></tr><tr><td>GPRC5B/Retinoic acid-induced gene 2 protein2</td><td>Unknown; retinoic acid induced</td><td>orphan</td><td>C Glutamate metabotropic</td><td>16p12.3</td></tr><tr><td>LGR4</td><td>stem cell proliferation/differentiation; Invasivness metastasis</td><td>R-spondin</td><td>A</td><td>11p14-p13</td></tr><tr><td>P2RY5/LPAR2</td><td>Sparce information. Involved in cancer</td><td>oleoyl-L-alpha-lysophosphatidic acid</td><td>A</td><td>13q14.2</td></tr></tbody></table>

GPCRs with expression specificity in the different subgroups of cells are listed together with i) their reported function in cancer or their physiological implications (if known), ii) their ligand(s), iii) the GPCR class they belong to and iv) their chromosome localization. T, O, U, N, stands for TG1, OB1, U-87 -MG and f-NSC, respectively. GPCRs in TO thus stands for the receptors specifically expressed in TG1 and OB1 GSCs compared to their expression in the three other cell types; TOU for those expressed in TG1, OB1 and U-87 MG and so on.
doi:10.1371/journal.pone.0091519.t003
GPRC5B (or RAIG-2), a retinoic acid-inducible orphan GPCR, is a member of the group C of metabotropic glutamate receptor family proteins with short N-termini. It is expressed in human brain [35,73,74]. Studies performed on Xenopus suggest a signaling through the non-canonical Wnt/frizzled pathway and a possible cross-talk between the two pathways [75]. This receptor may play an important role in brain maturation and neural circuitry building during brain development. Expression variants, due to alternative splicing have been reported in mouse brain [76]. GPRC5B-v1 is ubiquitously expressed, whereas GPRC5B-v2 is more specific to brain. Nevertheless, besides its implication in behavioral abnormalities, information concerning this receptor family is limited. GPRC5B is almost absent from U-87 MG cells. Its mRNA level is highest in f-NSCs, TG1 and OB1 cells with f-NSC>OB1~TG1 (Figure 1). Interestingly, mass spectrometry confirms its presence at the protein level with the same relative abundance in the three cell types.

CELSR2 (Cadherin EFG LAG seven-pass G-type receptor 2) is a non-classical cadherin adhesion receptor (group IV of adhesion GPCRs), member of the flamingo subfamily (fmi) characterized by the presence of a seven transmembrane domain [77]. The three members of the family, CELSR1-3, have multiple functions in epithelial and nervous systems and are among the core planar cell polarity (PCP) receptors [78,79,80]. A role in maintaining quiescence of long term hematopoietic stem cells has been reported [81]. They play a pivotal role in contact mediated communication. PCP regulation appears to take place through cooperation with the Frizzled receptors and the non-canonical Wnt pathway [81,82]. Aberrant regulation of these signaling pathways has been suggested to be involved in cancer invasion and metastasis [83,84]. Although expressed in HA cells, CELSR2 mRNA is overexpressed in cells with stem properties with levels in OB1>f-NSC>TG1. The same order of expression is also found at the protein level. The overexpression of this receptor in CSCs deserves further investigation.

Several other highly expressed GPCRs, namely CXCR4, GABA(B)Rs, LGR4 and CNR1, although not detected by mass spectroscopy, will be discussed.

CXCR4 presents the same relative expression as F2R/PAR1 at the mRNA level (Figure 4A) and similarities in functional properties. CXCR4 is a member of the chemokine receptor family and its endogenous ligand CXCL12 has been implicated in tumor development, including gliomas, facilitating angiogenesis, infiltration and metastasis through modulation of VEGF signaling, PAI-1 (plasminogen activator inhibitor-1) production and membrane type-2 matrix metalloproteinase expression (reviewed in [85,86]). In agreement with the results of the present study, overexpression of CXCR4 in GBMs has been recognized as an attribute of cells with progenitor/stem properties and CXCL12 promotes specifically the proliferation of these cells compared to more differentiated cancer cells [87,88]. Of note, the micro-RNA cluster miR-302-367 has recently been shown to inhibit the GSCs' self-renewal and infiltrative properties through CXCR4 repression [89]. Recent studies also show a role of the CXCR4/CXCL12 axis in triggering interaction of cancer cells with microvascular endothelial cells and GBM cancer stem cell transdifferentiation into pericytes [90,91]. The prominent role of this axis in tumor growth and cancer stem-like cells' interaction with their microenvironment and the perivascular space, designed CXCR4 as a potential target for adjuvant chemotherapies. So far the CXCR4 antagonist BKT140 is in phase I/IIA of a clinical trial for multiple myeloma (http://www.clinicaltrials.gov NCT01010880; for review see [86]). Expression of CXCR4 in neural stem cells has also recently been reported [92]. The present

study shows that the expression level of CXCR4 in f-NSCs can reach that observed in GSCs. Using antibodies, CXCR4 was found to be expressed in U-87 GBM cells. Our study indicates low expression of the receptor at the mRNA level in U-87 MG cells (below the cutoff level used in the study (Table S1)). Despite its high transcription level in f-NSCs, TG1 and OB1 cells, the protein could not be detected by mass spectrometry.

GABBR1 and GABBR2 (GPR51) which belong to the metabotropic GABA(B) receptors (class C GPCRs) present a relatively high expression in f-NSCs (Figures 5A and 5C). These two receptors are largely distributed in the central nervous system (CNS) and are particularly present on neurons where they play a role in neurotransmission both at the presynaptic and postsynaptic level. The two receptors function by forming heteromeric complexes of GABA(B) R1/R2 subunits [93,94,95,96], but homodimers [97], higher order oligomers [98] and complexes with other GPCRs have also been reported (e.g. the closely related extracellular calcium sensing receptor CaR [99], the muscarinic receptor M2 [100], as well as with the non GPCR GABA A receptors {Balasubramanian, 2004 #242}. GABA(A)R and GABA(B)R have been shown to exert opposite effects on rodent NSCs. GABA(A)R activation restricts NSC proliferation [101,102], whereas GABA(B)R activation promotes in vitro neurospheres' development [103]. Besides neurons, GABA(B) R1 and R2 are also present in astrocytes [104] and cells outside the CNS [105,106,107]. GABA(B) receptors have been involved in various cancers and the GABA(B) receptor agonist baclofen was shown to reduce cancer development in various animal models [108,109,110,111]. However, GABA signaling remains extremely complex and the multiple possible interaction of GABA(B) subunits with other partners, as well as the regulatory role played by the GABA(B)R2 in receptor trafficking [112] renders assessing a clear cut role of these receptors in cancer difficult. The high expression of the two GABA(B) receptors in f-NSCs may suggest a role in stem cells.

Similarly to GABBR2 (GPR51), LGR4 and CNR1 present an increased mRNA level in f-NSCs and OB1 cells.

LGR4/GPR48, a class A receptor, belongs to the family of leucine-rich repeat (LRR)-containing G protein-coupled receptors (or LGRs) composed of six members. LGR1-3 are hormone receptors. LGR4-6 were deorphanized in 2011 when De Lau and Carmon [113,114], working on intestinal crypt cells, showed that these receptors bind R-spondin, a molecule secreted by the crypt cells and known as a positive regulator of Wnt signaling. In the study of Carmon, LGR4 and 5 were also shown to interact with the Wnt- binding LRP/Frizzled receptor complex. LGR4-6 have been found in various tissues [115,116,117,118], where they regulate stem cell proliferation and differentiation. Inactivation of the LGR4 coding gene decreases cell proliferation and strongly reduces terminal differentiation of Paneth cells in postnatal mouse intestinal crypts [119]. Physiologically, LGR4 appears to play a key role in organism and tissue development [100,120,121,122]. Expression of LGR4 was also associated with invasiveness and metastatic activity of colon and cervical cancer cells [123,124]. So far no information is available on LGR4 expression in gliomas. Of note, in our study, among the three R-spondin binding receptors, only LGR4 is expressed at a significant level in TG1 and OB1 cells.

CNR1, which encodes the CB1 cannabinoid receptor (CB1), exhibits an expression profile comparable to LGR4 (Figure 5B). CB1 has a widespread distribution in various brain regions [125] and has two known endogenous ligands (anandamide and 2-arachidonyl-glycerol). Several genetic studies have been performed to link CNR1 mutations to psychiatric but also metabolic disorders
[126,127,128]. In rat astrocytes CB1 plays pivotal roles in energy metabolism and activation of this receptor increases the rate of glucose oxidation and ketogenesis, two mechanisms involved in energy supply of the CNS (for review see [125]). CNR1 expression has been reported in various glioma cell lines and receptor agonists trigger apoptosis and autophagy of these cells [129,130]. However, cannabinoids' effect is complex and relies on both ligand concentrations and receptor levels [131]. In glioma stem like cells, CB1 agonists have been reported to modulate the expression of stem genes [132]. In human glioma tissues, increased expression of cannabinoid receptors (CB1 and CB2) has been reported to correlate with higher tumor grades [133]; but this remains controversial [134,135,136]. In other tumors, a role of CNR1 in invasion and metastasis has been reported [137,138] and the endocannabinoid system has been reported as a promising tool to improve the efficacy of steroids in colon cancer [139]. Globally, the cannabinoid system appears to present anti-proliferative, anti-angiogenic and pro-apoptotic properties. Nevertheless, its potential use for treating glioma and may be cancer more generally needs further analysis and is hindered by side effects such as anxiety and depression.

Several less abundant GPCRs of the subgroup of genes expressed in cells with stem properties (Table 3), show specific upregulation in GSCs. These include CHRM3, GPR171 and P2RY5. As for F2R/PAR-1 and CXCR4, CHRM3 activation may synergistically act on the EGFR signaling pathway. Much less is known for P2RY5 and GPR171 and their relation to cancer.

GPCRs specifically expressed/overexpressed in TG1 and
OB1 GSCs

At the transcriptomic level, eight GPCRs appeared specifically upregulated and one downregulated in TG1 and OB1 GSCs (Figure 4 and Table 3) compared to their expression in HA or f-NSCs.

LPHN2 shows the highest expression among the different GPCRs of this group. The specific upregulation of LPHN2 at the mRNA level is corroborated at the protein level and both transcriptomic and proteomic data are in agreement with a level in TG1>OB1. LPHN2 is an adhesion GPCR (group I) belonging to the family of latrophilins composed of three members (LPHN1-3) (for review see [140]). Expression of latrophilin 1 and 3 is mainly found in the central nervous system, whereas LPHN2 is more ubiquitous [140,141,142,143]. LPHN2 shares with the two other members a large N-terminus composed of lectin, olfactomedin and hormone binding domains. The receptors bind α-latrotoxin, the black widow spider venom and this binding triggers strong exocytosis at neuron synapses. LPHN2 however shows lower binding affinities than the LPHN1 and LPHN3 isoforms. Not much is known about the functional role of these receptors in vertebrates. Through their large N-terminus, latrophilins may establish interactions with the extracellular matrix. Such interactions have been reported recently for LPHN1 [144,145] and LPHN3 [146]. Information on the function of LPHN2 is scarce. It has been reported to play a role in the epithelial mesenchymal transition (EMT) necessary for heart valve development during embryogenesis [143]. In addition, one report indicates a higher expression of the gene in breast tumors [147]. More recently, in a genome wide association study approach aimed at identifying genes involved in the response of tumor cells to paclitaxel, a microtubule cytotoxic used in chemotherapies, a SNP located in the intronic part of the LPHN2 coding gene was identified as being associated with the cellular response to this compound [148]. However there is no further information as to how an intronic SNP in LPHN2 may interfere with the cellular effect of

paclitaxel. In the current state of knowledge, one may speculate that the overexpression of LPHN2 in TG1 and OB1 GSCs may be related to EMT, the developmental program often activated during cancer metastasis and regulating the equilibrium between CSCs and the non-stem cancer cells [149]. A recent study using in silico analysis and aiming at finding genes under p53 control, designed LPHN2 as a novel target gene of p53 [150]. LPHN1 and LPHN3 are also upregulated in TG1 and OB1 cells, but both receptors appear less specific than LPHN2 and we did not detect them by mass spectrometry.

Other receptors, namely the parkin-associated endothelin B-like receptor (GPR37 or PAEL-R), the calcitonin receptor-like receptor (CALCRL or CRL), the histamine receptor HRH2, the prokinecitin receptor (GPR73 or PKR1), the endothelial differentiation G-protein coupled receptor 8 (EDG8 or S1P5), GPR128 and GPR103, show specific increase in their mRNA levels in TG1 and OB1 GSCs (Figure 4 and Table 3). However their expression is too low to be detected by mass spectrometry.

Concerning the histamine receptors, only changes in HRH2 could be observed in the present study. HRH1, H3 and H4 were expressed below the threshold level used in the study. HRH2 has a large tissue distribution and is involved in a plethora of physiological functions as reviewed in [151,152]. As other members of the histamine receptor family [153] it has been reported to play a role in cell proliferation, differentiation energy metabolism and thus to be potentially involved in cancer. HRH2 are found in human mammary and gastric carcinoma cells and several human melanoma cell lines and the HRH2 antagonist cimetidine appears to be endowed with anticancer properties, although its action is not exclusively mediated through HRH2 [154,155]. Clinical trials using HRH2 antagonists as adjuvant treatment in resected colorectal cancer have also been reported ([156]; for review see [157]), but their real benefits need to be assessed. One additional report points to an antiproliferative activity of histamine mediated through HRH2 in pancreatic cells (PANC-1), with a G0/G1-phase arrest and a modulation of the Bcl-2 family proteins [158]. In our hands, two HRH2 agonists (anthamine and dimaprit) and five antagonists (cimetidine, zolandine, ranitidine, thiotidine and famotidine) tested at 5 and 50 µM for 24 hours did not show a significant effect on TG1 GSCs' proliferation and survival. The specific overexpression of this histamine receptor subtype in TG1 and OB1 GSCs suggests its potential involvement in the cells' pathophysiology that needs further and deeper investigation.

The function of the calcitonin receptor-like receptor CALCRL depends on the sub-type of RAMP (receptor-activity modulating protein) to which it is associated. It functions either as a receptor for calcitonin-gene-related peptide (CGRP) when associated to RAMP1 or for adrenomodullin (AM) when associated to RAMP2 or RAMP3 [159]. AM has vasorelaxation properties and its expression is controlled by the hypoxia inducible factor (HIF). Cells overexpressing AM exhibit higher resistance to hypoxia induced apoptosis [160]. AM induces overexpression of anti-apoptotic factors and decreases expression of pro-apoptotic factors [160,161]. The location of AM close to necrotic zones in tumors suggests a major role of this peptide in the resistance of tumor cells to hypoxia. Numerous data are available on AM and its expression in tumors. However data on its receptor are scarce. A recent work shows that CALCRL and its partners RAMP2 and RAMP3 are present at the mRNA and protein levels in GBM tumors isolated from several patients [162]. Moreover, systemic injections of AM receptor antibodies in animals were shown to inhibit both tumor growth and neo-vascularization of U87 xenografted cells in a dose dependent manner [163]. Nevertheless, information on the
distribution, regulation and function of endogenous AM receptors
remains poor. Overall, data from the literature suggest a role of
AM and its receptor in the survival of tumor cells within the
hypoxic regions of solid tumors. AM may induce neovasculariza-
tion and promote tumor growth. However, in the absence of
functional data, it remains difficult to assess the exact relationship
between AM, hypoxia, cell proliferation, cell survival and
CALCR.

Concerning EDG8/S1P5, the role of the ligand S1P (sphingosine 1-phosphate) in cancer has recently been reviewed [164]. For the other receptors of this subgroup, information relative to their implication in cancer is rare (one report for GPR37 [165]) or lacking.

Concluding Remarks

GPCRs, which mediate extracellular to intracellular signaling,
are key receptors in cell functioning in both physiological and
disease conditions. Over the years, a great number of GPCRs have
successfully been developed as pharmacological targets in various
diseases. However, the complexity of GPCR signaling offered by the
possible formation of homo- or hetero oligomers, interactions with a
wealth of other proteins, including growth factor receptors
[166,167,168,169] and downstream effectors, and its regulatory
role in several pathways including the Hippo pathway which
coordinates tissue homeostasis [170], unveils bewildering opportu-
nities to target new diseases and fuel drug discovery [31]. Oncology
is one of the sectors that may benefit from these new developments
as deregulation of GPCR-dependent cell signaling sustains and may
even drive tumorigenesis [32,35,36]. Exploration of the compara-
tive expression of 356 GPCRs (among which 138 were retained in
the present study) in GSCs and in non-cancer neural cells with stem
cell properties (human fetal neural stem cells) or without stem cell
properties (human astrocytes) and to the established GBM U-87
MG cell line, unveils differences in GPCR expression between the
five cell types of neural origin. Several genes, including GPR37,
GPR73 (PKR1), GPR103, GPR128, EDG8 (S1P5) and HRH2
appear as the most specifically enriched in TG1 and OB1 cells.
LPHN2, which is expressed in the five cell types, is predominantly
upregulated in TG1 and OB1 cells. LPHN2 is the only GPCR of
the group detected by mass spectrometry. The specific under-
expression of FZD6 in TG1 and OB1 cells is noteworthy as FZD6 is
the only GPCR showing such a profile. The GPCRs emerging from
our analysis are related to metabolic functions, pro-angiogenic
activities, cell proliferation, resistance to hypoxia and interactions
with the microenvironment (Table 3), an ensemble of characteristics
which endow these cancer cells with properties allowing them to be
maintained and develop in the specific microenvironment created
by the tumor and which may induce death of the normal neural
cells. GSCs may originate from differentiated cells that acquired
stem properties or from stem/immature progenitor cells that lost
their normal regulatory capacities to become cancerous. Whatever
the origin, the GSCs may have acquired or retained gene expression
found in normal stem/progenitor cells. In this GPCR transcriptome
study, several GPCRs show a differential expression pattern in f-
NSCs, TG1 and OB1 cells compared to HA or U-87 MG cells.
Some GPCRs are already known for their involvement in GBMs.
Others have been reported for their involvement in different
cancers, but not in brain cancers. To our knowledge, this is the first
exhaustive report presenting comparable expression profiles of the
whole set of GPCRs in cells from GBM tumors and non-cancerous
tissues at the transcriptomic and proteomic levels. Not many
GPCRs have so far been reported for their role in GSCs. For those
genes that are already documented, our findings mostly correlate

with the reported studies obtained using other techniques such as
use of agonists and antagonists, or knockout/knockdown experi-
ments. GPCRs specific to GSCs or U-87 MG cells constitute the
most appealing ones for cancer therapy. In addition to changes in
the level of expression, genomic modifications with mutations that
can alter the integrity of the receptor, namely its structure,
interaction with partner proteins, affinity for its ligand(s) and
downstream signaling have to be considered [171]. Protein altering
mutations affecting several GPCRs as well as G protein subunits,
have been reported in a study of 441 breast, lung, ovarian,
pancreatic, prostate cell types and subtypes [172]. Targeting the
finely tuned-GPCR dependent signaling with biased agonists or
antagonists or allosteric effectors [173] selected to act on the
identified deregulated pathways opens new avenues for more
efficient and personalized multi-target oriented chemotherapies.
Getting deeper insight into tumor pathophysiology, identifying
cancer subtypes and developing system biology approaches to
simulate cell intra- and intercellular communication networks will
pave the route for disruptive innovation in oncology. The present
descriptive presentation of GPCRs exhibiting specific expression
changes in GCSs is a first step that calls for future studies at the
functional level in order to decipher the precise role of these
receptors in sustaining the cancerous phenotype.

**Materials and Methods**

Ethics statement

The biomedical research was conducted according to the
declaration of Helsinki, to the French laws and was approved by
the institutional review board of Ste Anne Hospital, Paris. Patients
have given written informed consent. Isolation and characteriza-
tion of neural stem cells from human fetal brain at embryonic day
50-55 (Carnegie stage 19-22) were performed under ethical
approval from the University Paris-Descartes internal review
board using tissue donated with written informed consent after
elective termination of pregnancy.

Cell Cultures

Patients' derived cells were isolated from surgical resections of human primary GBMs provided by the Neurosurgery department of Sainte Anne Hospital in Paris. Patients' tumors were diagnosed as GBMs according to the World Health Organization (WHO) classification system [174] and as malignant glio-neuronal tumor cells according to the Sainte Anne's classification and grading system [175]. Malignant glio-neuronal tumor cells (TG1 and OB1 cells) with stem cell properties were isolated from those biopsies as described in [27]. Their characterization was previously reported [27,176]. *TP53* and *IDH1* were not mutated in these cells and SHH and NOTCH pathways were active. TG1 and OB1 cells were able to differentiate and to induce tumors when orthotopically grafted (<500 cells) in immune-compromised mice. Cells were expanded in Dulbecco's modified Eagle's: F-12 medium (1:1) supplemented with N2, G5 (containing bFGF and EGF) and B27 (all from Invitrogen, France) at a density of 3.10⁴ cells/cm². Human fetal neural stem cells (f-NSCs/NSC XXIV) were derived from human fetal brain and characterized as previously described [177]. Cells were cultured in "Neurocult"® NSC Proliferation Medium supplemented with basic fibroblast growth factor (10 ng/ml, bFGF, Peprotech, Rocky Hill, NJ) and with epidermal growth factor (20 ng/ml, EGF, Peprotech, Rocky Hill, NJ). Under the culture conditions used, TG1, OB1 and f-NSC cells grew as neurospheres. For cell expansion, once a week, neurospheres were mechanically dissociated into a single-cell suspension and 90% of their culture medium was renewed. Human Astrocyte primary
cells (HA cells, ScienCell Research Laboratories) were expanded in Astrocyte Medium (ScienCell Research Laboratories) according to the manufacturer's instructions. Human brain tumor cells U-87 MG (American Type Culture Collection, ATCC) were expanded in ATCC complete growth medium according to the manufacturer's instructions. For all cell lines, Master Cell Banks were established and Working Cell Banks were derived from these Master Cell Banks for further experiments.

## G-protein-coupled receptor transcriptome

Total RNA was isolated from $10^7$ cells of each cell type with TRI Reagent™ (Euromedex, France) according to the manufacturer's instructions and further purified with RNeasy mini kit columns (Qiagen). RNA purity was assessed by absorption spectra analysis using a NanoDrop ND-1000 apparatus (Labtech). Absorption ratios A260/A280 and A260/A230 were comprised between 1.8 and 2.1. RNA integrity and concentration were checked using an Agilent 2100 Bioanalyzer and the RNA 6000 LabChip kit. Only RNA with a RNA Integrity Number (RIN) higher than 8.5 was processed (2100 expert software, Agilent Technologies). 10 µg of total RNA were reverse transcribed to single-stranded cDNA using the High Capacity cDNA Reverse Transcription kit (Applied Biosystems). Real-time PCR analysis was performed using TaqMan® Low Density Arrays (Applied Biosystems). TaqMan® GPCR Gene Expression Arrays correspond to quantitative PCR (Q-PCR) assays spotted onto a 384-well card. cDNA (850 ng) was carefully mixed with 850 µL of TaqMan Universal PCR Master Mix (PE Applied Biosystems). Real-time RT-PCR was performed using an ABI Prism 7900HT apparatus (Applied Biosystems, Foster City, CA). Thermal cycling was carried out for 2 min at 50°C (to activate uracil-DNA glycosylase), 10 min at 94.5°C (activation) and 40 cycles of 30 s at 97°C and 60 s at 59.7°C. Each reaction well contained cDNA derived from about 1.5 ng of total RNA. Experiments were performed in triplicates for f-NSCs and duplicates for the four other cell types.

## Cell membrane preparation, protein extraction and separation

Cells grown as neurospheres were dissociated, washed in PBS and suspended in buffer A consisting of 10 mM Hepes buffer pH 7.4, 2 mM EDTA and protease inhibitors (complete proteases inhibitors cocktail from Roche). Cells were mechanically lysed by pipetting the solution up and down with a P1000 pipette and pressing the solution through a 22G needle mounted on a syringe several times (~15) until 90% of the cells were lysed. Cell debris were removed by a 10 min. centrifugation at 1560 g. The supernatant was collected and centrifuged for 90 min at 20 800 g. Cell ghosts were resuspended in buffer A. After denaturation at 100°C in 5% SDS, 5% βmercaptoethanol, 1 mM EDTA, 10% glycerol, 10 mM Tris pH 8 buffer for 3 min, protein samples were fractionated on 1D SDS-PAGE gradient gel (8–15%). Gels were stained with colloidal Coomassie Blue. 90 slices were cut and stored at -20°C prior to mass spectrometry analysis.

## Mass spectrometry proteomic analysis

*In gel* digestion of gel slices was performed as previously described [178]. The resulting peptide extracts were analysed by C18 reversed phase nanoHPLC on a nanoHPLC-Chip/MS system (Agilent Technologies, Palo Alto, USA) coupled to an ion trap amaZon (Bruker Daltonics, Bremen, Germany) mass spectrometer. Peptide mixtures were loaded on a 40 nl enrichment column and separated on a 75 µm × 150 mm column packed with Zorbax 300SB-C18 5 µm material (Agilent Technologies, Palo

Alto, USA). Trapping was performed during 3 min, at 3.75 µL/min, with solvent A (98% H₂O, 2% ACN and 0.1% FA). Elution was performed using two solvents, A (98% H₂O, 2% ACN and 0.1% FA) and B (98% ACN, 2% H₂O and 0.1% FA) at a flow rate of 300 nL/min, using a 2–24% B gradient over 45 min and a 24–40% B gradient over 16 min, followed by 70% B over 5 min.

For tandem MS experiments, the system was operated with automatic switching between MS and MS/MS modes. The 8 most abundant peptides, preferably doubly and triply charged ions, were selected on each MS spectrum for further isolation and CID fragmentation. Then, selected ions were excluded during 0.6 min. The MS/MS scanning was performed in the ultrascan resolution mode at a scan rate of 32 500 m/z per second. The complete system was fully controlled by Hystar 3.2 (Bruker Daltonics).

## Proteomic Data Analysis

Raw data collected during nanoLC-MS/MS analyses were processed and converted into \*.mgf peak list format with DataAnalysis 4.0 (Bruker Daltonics). MS/MS data were interpreted using two search engines: Mascot (version 2.3.02, Matrix Science, London, UK) installed on a local server and OMSSA (Open Mass Spectrometry Search Algorithm version 2.1.9; [179]) interfaced for grid computing on an in-house developed and freely available software suite (https://msda.unistra.fr). Searches were performed against a composite target-decoy UniProtKB-SwissProt human entries (TaxID = 9606, 10 July 2012, 40500 total entries). Searches were performed with a tolerance on mass measurement of 0.2 Da for precursor and 0.3 Da for fragment ions. Carbamidomethylation of cysteine residues, oxidation of methionine residues and propionamidation of cysteine residues were searched as variable modifications. Up to one trypsin missed cleavage was allowed. Scaffold software (version 3.6.5 Proteome Science, Portland, USA) was used for validation and false positive rate estimation of protein identifications. For each sample and each search engine, peptides were filtered out according to the cutoff set for protein hits with 1 or more peptides and a false positive identification rate less of than 1%. Further data manipulations were performed with Microsoft Excel. Semi-quantification was performed by peptide and spectral counting.

## Supporting Information

<caption>**Figure S1** Expression of 14 housekeeping genes present on the TaqMan GPCR Array Card (Applied Biosystems).</caption>

Expression distribution for each housekeeping gene, in the different experiments performed, is given as box plots. Ordinates represent the cycle threshold (Ct) of the quantitative PCR experiment for each housekeeping gene indicated. The middle line in the boxes corresponds to the median. The end of the wishers corresponds to the lowest and highest values within 1.5 IQR (interquartile range) of the lower and upper quartile. Outliers are represented by stars. Abbreviations: 18 s, 18 s ribosomal RNA; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; POLR2A, DNA dependent RNA polymerase 2A subunit; ACTB, actin β; PPIA, peptidylpropyl isomerase A; PGK1, phosphoglycerate kinase 1; B2M, β2-microglobulin; GUSB, glucuronidase β; HPRT1, hypoxanthine phosphoribosyltransferase 1; TBP, TATA-box binding protein; TFRC, transferrin receptor; HMBS, hydroxymethylbilane synthase; IPO8, importin 8; PHGDH, phosphoglycerate dehydrogenase. (DOCX)

<caption>**Table S1** Gene assays present on the GPCR microarray plates and their expression values in five cell types.</caption>

381 different assays are present on the GPCR microarray, allowing the analysis of the expression of 356 GPCRs, 14 different housekeeping genes (rRNA 18 s is in quadruplicate) and 11 other genes
(among which two genes of the LANCL family). Expression data are presented in four groups. Group I corresponds to GPCRs expressed with Ct values ≤31.5 (138 genes) in at least one of the five cell types, group II to GPCRs expressed at Ct>31.5 in all five cell types, group III to the housekeeping genes (14 genes) and group IV to eleven genes not related to housekeeping genes or GPCRs. In each group, genes are listed in alphabetical order. Expression values are indicated as cycle threshold values (Ct) and as 2^Ct * 10^12 (arbitrary units). The cutoff level was set at Ct ≤31.5 (i.e. expression level ≥329 units). For GPCRs which did not show amplification, Ct values were arbitrarily set at 40. (XLSX)

## References

1. Dolecek TA, Propp JM, Stroup NE, Kruchko C (2012) CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2005–2009. *Neuro Oncol* 14 Suppl 5: v1–49.

2. Schechter CB (1999) Re: Brain and other central nervous system cancers: recent trends in incidence and mortality. *J Natl Cancer Inst* 91: 2050–2051.

3. Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, et al. (2005) Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. *N Engl J Med* 352: 987–996.

4. Holdhoff M, Grossman SA (2011) Controversies in the adjuvant therapy of high-grade gliomas. *Oncologist* 16: 351–358.

5. Yin AA, Cheng JX, Zhang X, Liu BL (2013) The treatment of glioblastomas: a systematic update on clinical Phase III trials. *Crit Rev Oncol Hematol* 87: 265–282.

6. Stupp R, Hegi ME, Mason WP, van den Bent MJ, Taphoorn MJ, et al. (2009) Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. *Lancet Oncol* 10: 459–466.

7. Bonnet D, Dick JE (1997) Human acute myeloid leukemia is organized as a hierarchy that originates from a primitive hematopoietic cell. *Nat Med* 3: 730–737.

8. Reya T, Morrison SJ, Clarke MF, Weissman IL (2001) Stem cells, cancer, and cancer stem cells. *Nature* 414: 105–111.

9. Hemmati HD, Nakano I, Lazareff JA, Masterman-Smith M, Geschwind DH, et al. (2003) Cancerous stem cells can arise from pediatric brain tumors. *Proc Natl Acad Sci U S A* 100: 15178–15183.

10. Singh SK, Clarke ID, Terasaki M, Bonn VE, Hawkins C, et al. (2003) Identification of a cancer stem cell in human brain tumors. *Cancer Res* 63: 5821–5828.

11. Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF (2003) Prospective identification of tumorigenic breast cancer cells. *Proc Natl Acad Sci U S A* 100: 3983–3988.

12. Dalerba P, Cho RW, Clarke MF (2007) Cancer stem cells: models and concepts. *Annu Rev Med* 58: 267–284.

13. O'Brien CA, Pollett A, Gallinger S, Dick JE (2007) A human colon cancer cell capable of initiating tumour growth in immunodeficient mice. *Nature* 445: 106–110.

14. Schatton T, Frank MH (2008) Cancer stem cells and human malignant melanoma. *Pigment Cell Melanoma Res* 21: 39–55.

15. Schatton T, Murphy GF, Frank NY, Yamaura K, Waaga-Gasser AM, et al. (2008) Identification of cells initiating human melanomas. *Nature* 451: 345–349.

16. Hermann PC, Huber SL, Herrler T, Aicher A, Ellwart JW, et al. (2007) Distinct populations of cancer stem cells determine tumor growth and metastatic activity in human pancreatic cancer. *Cell Stem Cell* 1: 313–323.

17. Chen J, Li Y, Yu TS, McKay RM, Burns DK, et al. (2012) A restricted cell population propagates glioblastoma growth after chemotherapy. *Nature*.

18. Driessens G, Beck B, Caauwe A, Simons BD, Blanpain C (2012) Defining the mode of tumour growth by clinical analysis. *Nature*.

19. Schepers AG, Snippert HJ, Stange DE, van den Born M, van Es JH, et al. (2012) Lineage tracing reveals Lgr5+ stem cell activity in mouse intestinal adenomas. *Science* 337: 730–735.

20. Frank NY, Schatton T, Frank MH (2010) The therapeutic promise of the cancer stem cell concept. *J Clin Invest* 120: 41–50.

21. Alison MR, Lin WR, Lim SM, Nicholson LJ (2012) Cancer stem cells: in the line of fire. *Cancer Treat Rev* 38: 589–598.

22. Sharkis SJ, Jones RJ, Civin C, Jang YY (2012) Pluripotent stem cell-based cancer therapy: promise and challenges. *Sci Transl Med* 4: 127ps129.

23. Lee C, Dunn SE, Yip S (2012) Stem cells in brain tumour development and therapy- two-sides of the same coin. *Can J Neurol Sci* 39: 145–156.

24. Nduom EK, Hadjipanayis CG, Van Meir EG (2012) Glioblastoma cancer stem-like cells: implications for pathogenesis and treatment. *Cancer J* 18: 100–106.

25. Malik B, Nie D (2012) Cancer stem cells and resistance to chemo and radio therapy. *Front Biosci* (Elite Ed) 4: 2142–2149.

## Acknowledgments

We thank Dr René Frydman (Hopital Antoine Béclère, Clamart, France) and Dr Hervé Coffiny for providing fetal brain tissues and the Sainte Anne neurosurgery team for providing glioblastoma tissues.

## Author Contributions

Conceived and designed the experiments: JMS CC AVD HC SC JH MCK. Performed the experiments: MF JMS MZ SL JD. Analyzed the data: MF JMS MZ SL CC JD JH MCK. Contributed reagents/materials/ analysis tools: MPJ HC. Wrote the paper: JMS SL CC SC JH MCK. Revised the paper critically: HC AVD MPJ HC.

26. Persano L, Rampazzo E, Basso G, Viola G (2013) Glioblastoma cancer stem cells: role of the microenvironment and therapeutic targeting. *Biochem Pharmacol* 85: 612–622.

27. Patru C, Romao L, Varlet P, Coulombel L, Raponi E, et al. (2010) CD133, CD15/SSEA-1, CD34 or side populations do not resume tumor-initiating properties of long-term cultured cancer stem cells from human malignant glioma neuronal tumors. *BMC Cancer* 10: 66.

28. Drews J (2000) Drug discovery: a historical perspective. *Science* 287: 1960–1964.

29. Hopkins AL, Groom CR (2002) The druggable genome. *Nat Rev Drug Discov* 1: 727–730.

30. Jacoby E, Bouhelal R, Gerspacher M, Seuwen K (2006) The 7 TM G-protein-coupled receptor target family. *ChemMedChem* 1: 761–782.

31. Garland SL (2013) Are GPCRs Still a Source of New Targets? *J Biomol Screen*.

32. Lappano R, Maggiolini M (2011) G protein-coupled receptors: novel targets for drug discovery in cancer. *Nat Rev Drug Discov* 10: 47–60.

33. Herr DR (2012) Potential use of G protein-coupled receptor-blocking monoclonal antibodies as therapeutic agents for cancers. *Int Rev Cell Mol Biol* 297: 45–81.

34. Dorsam RT, Gutkind JS (2007) G-protein-coupled receptors and cancer. *Nat Rev Cancer* 7: 79–94.

35. Audigier Y, Picault FX, Chaves-Almagro G, Masri B (2013) G Protein-Coupled Receptors in cancer: biochemical interactions and drug design. *Prog Mol Biol Transl Sci* 115: 143–173.

36. Feigin ME (2013) Harnessing the genome for characterization of GPCRs in cancer pathogenesis. *FEBS J*.

37. Axelson M, Liu K, Jiang X, He K, Wang J, et al. (2013) U.S. Food and Drug Administration approval: vismodegib for recurrent, locally advanced, or metastatic basal cell carcinoma. *Clin Cancer Res* 19: 2289–2293.

38. Webb DR, Handel TM, Kretz-Rommel A, Stevens RC (2013) Opportunities for functional selectivity in GPCR antibodies. *Biochem Pharmacol* 85: 147–152.

39. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. *Genome Biol* 3: RESEARCH0034.

40. Valente V, Teixeira SA, Neder L, Okamoto OK, Oba-Shinjo SM, et al. (2009) Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR. *BMC Mol Biol* 10: 17.

41. Mackrill JJ (2004) Generation, use, and validation of receptor-selective antibodies. *Methods Mol Biol* 259: 47–65.

42. Jensen BC, Swigart PM, Simpson PC (2009) Ten commercial antibodies for alpha-1-adrenergic receptor subtypes are nonspecific. *Naunyn Schmiedebergs Arch Pharmacol* 379: 409–412.

43. Bodei S, Arrighi N, Spano P, Sigala S (2009) Should we be cautious on the use of commercially available antibodies to dopamine receptors? *Naunyn Schmiedebergs Arch Pharmacol* 379: 413–415.

44. Hamdani N, van der Velden J (2009) Lack of specificity of antibodies directed against human beta-adrenergic receptors. *Naunyn Schmiedebergs Arch Pharmacol* 379: 403–407.

45. Pradidarcheep W, Stallen J, Labruyere WT, Dabhoiwala NF, Michel MC, et al. (2009) Lack of specificity of commercially available antisera against muscarinic and adrenergic receptors. *Naunyn Schmiedebergs Arch Pharmacol* 379: 397–402.

46. Beermann S, Seifert R, Neumann D (2012) Commercially available antibodies against human and murine histamine H(4)-receptor lack specificity. *Naunyn Schmiedebergs Arch Pharmacol* 385: 125–135.

47. Aust G, Eichler W, Laue S, Lehmann I, Heldin NE, et al. (1997) CD97: a dedifferentiation marker in human thyroid carcinomas. *Cancer Res* 57: 1798–1806.

48. Chidambaram A, Fillmore HL, Van Meter TE, Dumur CI, Broaddus WC (2012) Novel report of expression and function of CD97 in malignant gliomas: correlation with Wilms tumor 1 expression and glioma cell invasiveness. *J Neurosurg* 116: 843–853.
49. Ward Y, Lake R, Yin JJ, Heger CD, Raffeld M, et al. (2011) LPA receptor heterodimerizes with CD97 to amplify LPA-initiated RHO-dependent signaling and invasion in prostate cancer cells. Cancer Res 71: 7301–7311.

50. Ward Y, Lake R, Martin PL, Killian K, Salerno P, et al. (2013) CD97 amplifies LPA receptor signaling and promotes thyroid cancer progression in a mouse model. Oncogene 32: 2726–2738.

51. Safae M, Clark AJ, Oh MC, Ivan ME, Bloch O, et al. (2013) Overexpression of CD97 confers an invasive phenotype in glioblastoma cells and is associated with decreased survival of glioblastoma patients. PLoS One 8: e62765.

52. Liu D, Trojanowicz B, Radestock Y, Fu T, Hammje K, et al. (2010) Role of CD97 isoforms in gastric carcinoma. Int J Oncol 36: 1401–1408.

53. Shashidhar S, Lorente G, Nagavarapu U, Nelson A, Kuo J, et al. (2005) GPR56 is a GPCR that is overexpressed in gliomas and functions in tumor cell adhesion. Oncogene 24: 1673–1682.

54. Xu L, Begum S, Hearn JD, Hynes RO (2006) GPR56, an atypical G protein-coupled receptor, binds tissue transglutaminase, TG2, and inhibits melanoma tumor growth and metastasis. Proc Natl Acad Sci U S A 103: 9023–9028.

55. Xu L, Begum S, Barry M, Crowley D, Yang L, et al. (2010) GPR56 plays varying roles in endogenous cancer progression. Clin Exp Metastasis 27: 241–249.

56. Yang L, Chen G, Mohanty S, Scott G, Fazal F, et al. (2011) GPR56 Regulates VEGF production and angiogenesis during melanoma progression. Cancer Res 71: 5558–5568.

57. Luo R, Jeong SJ, Jin Z, Strokes N, Li S, et al. (2011) G protein-coupled receptor 56 and collagen III, a receptor-ligand pair, regulates cortical development and lamination. Proc Natl Acad Sci U S A 108: 12925–12930.

58. Terskikh AV, Easterday MC, Li L, Hood L, Kornblum HI, et al. (2001) From hematopoiesis to neuropoiesis: evidence of overlapping genetic programs. Proc Natl Acad Sci U S A 98: 7934–7939.

59. Bai Y, Du L, Shen L, Zhang Y, Zhang L (2009) GPR56 is highly expressed in neural stem cells but downregulated during differentiation. Neuroreport 20: 918–922.

60. Cork SM, Van Meir EG (2011) Emerging roles for the BAI1 protein family in the regulation of phagocytosis, synaptogenesis, neurovasculature, and tumor development. J Mol Med (Berl) 89: 743–752.

61. Kaur B, Brat DJ, Calkins CC, Van Meir EG (2003) Brain angiogenesis inhibitor 1 is differentially expressed in normal brain and glioblastoma independently of p53 expression. Am J Pathol 162: 19–27.

62. Kang X, Xiao X, Harata M, Bai Y, Nakazaki Y, et al. (2006) Antiangiogenic activity of BAI1 in vivo: implications for gene therapy of human glioblastomas. Cancer Gene Ther 13: 385–392.

63. Cork SM, Kaur B, Devi NS, Cooper L, Saltz JH, et al. (2012) A proprotein convertase/MMP-14 proteolytic cascade releases a novel 40 kDa vasculostatin from tumor suppressor BAI1. Oncogene.

64. Okajima D, Kudo G, Yokota H (2010) Brain-specific angiogenesis inhibitor 2 (BAI2) may be activated by proteolytic processing. J Recept Signal Transduct Res 30: 143–153.

65. Kee HJ, Ahn KY, Choi KC, Won Song J, Heo T, et al. (2004) Expression of brain-specific angiogenesis inhibitor 3 (BAI3) in normal brain and implications for BAI3 in ischemia-induced brain angiogenesis and malignant glioma. FEBS Lett 569: 307–316.

66. Garcia-Lopez MT, Gutierrez-Rodriguez M, Herranz R (2010) Thrombin-activated receptors: promising targets for cancer therapy? Curr Med Chem 17: 109–128.

67. Elste AP, Petersen I (2010) Expression of proteinase-activated receptor 1–4 (PAR 1–4) in human cancer. J Mol Histol 41: 89–99.

68. Arora P, Cuevas BD, Russo A, Johnson GL, Trejo J (2008) Persistent transactivation of EGFR and ErbB2/HER2 by protease-activated receptor-1 promotes breast carcinoma cell invasion. Oncogene 27: 4434–4445.

69. Gessler F, Voss V, Dutzmann S, Seifert V, Gerlach R, et al. (2010) Inhibition of tissue factor/protease-activated receptor-2 signaling limits proliferation, migration and invasion of malignant glioma cells. Neuroscience 165: 1312–1322.

70. Dutzmann S, Gessler F, Harter PN, Gerlach R, Mittelbronn M, et al. (2010) The pro-migratory and pro-invasive role of the procoagulant tissue factor in malignant gliomas. Cell Adh Migr 4: 515–522.

71. Magnus N, Garnier D, Rak J (2010) Oncogenic epidermal growth factor receptor up-regulates multiple elements of the tissue factor signaling pathway in human glioma cells. Blood 116: 815–818.

72. Zigler M, Kamiya T, Brantley EC, Villares GJ, Bar-Eli M (2011) PAR-1 and thrombin: the ties that bind the microenvironment to melanoma metastasis. Cancer Res 71: 6561–6566.

73. Robbins MJ, Michalovich D, Hill J, Calver AR, Medhurst AD, et al. (2000) Molecular cloning and characterization of two novel retinoic acid-inducible orphan G-protein-coupled receptors (GPRC5B and GPRC5C). Genomics 67: 8–18.

74. Brauner-Osborne H, Krogsgaard-Larsen P (2000) Sequence and expression pattern of a novel human orphan G-protein-coupled receptor, GPRC5B, a family C receptor with a short amino-terminal domain. Genomics 65: 121–128.

75. Harada Y, Yokota C, Habas R, Slusarski DC, He X (2007) Retinoic acid-inducible G protein-coupled receptors bind to frizzled receptors and may activate non-canonical Wnt signaling. Biochem Biophys Res Commun 358: 968–975.

76. Cool BH, Chan GC, Lee L, Oshima J, Martin GM, et al. (2010) A flanking gene problem leads to the discovery of a Gprc5b splice variant predominantly expressed in C57Bl/6J mouse brain and in maturing neurons. PLoS One 5: e10351.

77. Usui T, Shima Y, Shimada Y, Hirano S, Burgess RW, et al. (1999) Flamingo, a seven-pass transmembrane cadherin, regulates planar cell polarity under the control of Frizzled. Cell 98: 585–595.

78. Boutin C, Goffinet AM, Tissir F (2012) Celsr1-3 cadherins in PCP and brain development. Curr Top Dev Biol 101: 161–183.

79. Tissir F, Goffinet AM (2006) Expression of planar cell polarity genes during development of the mouse CNS. Eur J Neurosci 23: 597–607.

80. Tissir F, Qu Y, Montcouqiol M, Zhou L, Komatsu K, et al. (2010) Lack of cadherins Celsr2 and Celsr3 impairs ependymal ciliogenesis, leading to fatal hydrocephalus. Nat Neurosci 13: 700–707.

81. Sugimura R, He XC, Venkatraman A, Arai F, Box A, et al. (2012) Noncanonical Wnt signaling maintains hematopoietic stem cells in the niche. Cell 150: 351–365.

82. Wada H, Tanaka H, Nakayama S, Iwasaki M, Okamoto H (2006) Frizzled3a and Celsr2 function in the neuroepithelium to regulate migration of facial motor neurons in the developing zebrafish hindbrain. Development 133: 4749–4759.

83. Katoh M (2005) WNT/PCP signaling pathway and human cancer (review). Oncol Rep 14: 1583–1588.

84. Huang H, Groth J, Sossey-Alaoui K, Hawthorn L, Beall S, et al. (2005) Aberrant expression of novel and previously described cell membrane markers in human breast cancer cell lines and tumors. Clin Cancer Res 11: 4357–4364.

85. Terasaki M, Sugita Y, Arakawa F, Okada Y, Ohshima K, et al. (2011) CXCL12/CXCR4 signaling in malignant brain tumors: a potential pharmacological therapeutic target. Brain Tumor Pathol 28: 89–97.

86. Domanska UM, Kruizinga RC, Nagengast WB, Timmer-Bosscha H, Huls G, et al. (2013) A review on CXCR4/CXCL12 axis in oncology: no place to hide. Eur J Cancer 49: 219–230.

87. Ehtesham M, Mapara KY, Stevenson CB, Thompson RC (2009) CXCR4 mediates the proliferation of glioblastoma progenitor cells. Cancer Lett 274: 305–312.

88. Schulte A, Gunther HS, Phillips HS, Kemming D, Martens T, et al. (2011) A distinct subset of glioma cell lines with stem cell-like properties reflects the transcriptional phenotype of glioblastomas and overexpresses CXCR4 as therapeutic target. Glia 59: 590–602.

89. Fareh M, Turchi L, Virolle V, Debruyne D, Almairac F, et al. (2012) The miR 302–367 cluster drastically affects self-renewal and infiltration properties of glioma-initiating cells through CXCR4 repression and consequent disruption of the SHH-GLI-NANOG network. Cell Death Differ 19: 232–244.

90. Cheng L, Huang Z, Zhou W, Wu Q, Donnola S, et al. (2013) Glioblastoma stem cells generate vascular pericytes to support vessel function and tumor growth. Cell 153: 139–152.

91. Rao S, Sengupta R, Choe EJ, Woerner BM, Jackson E, et al. (2012) CXCL12 mediates trophic interactions between endothelial and tumor cells in glioblastoma. PLoS One 7: e33005.

92. Li M, Chang CJ, Lathia JD, Wang L, Pacenta HL, et al. (2011) Chemokine receptor CXCR4 signaling modulates the growth factor-induced cell cycle of self-renewing and multipotent neural progenitor cells. Glia 59: 108–118.

93. Jones KA, Borowsky B, Tamm JA, Craig DA, Durkin MM, et al. (1998) GABA(B) receptors function as a heteromeric assembly of the subunits GABA(B)R1 and GABA(B)R2. Nature 396: 674–679.

94. White JH, Wise A, Main MJ, Green A, Fraser NJ, et al. (1998) Heterodimerization is required for the formation of a functional GABA(B) receptor. Nature 396: 679–682.

95. Kaupmann K, Malitschek B, Schuler V, Heid J, Froestl W, et al. (1998) GABA(B)-receptor subtypes assemble into functional heteromeric complexes. Nature 396: 683–687.

96. Billinton A, Ige AO, Wise A, White JH, Disney GH, et al. (2000) GABA(B) receptor heterodimer-component localisation in human brain. Brain Res Mol Brain Res 77: 111–124.

97. Villemure JF, Adam L, Bevan NJ, Gearing K, Chenier S, et al. (2005) Subcellular distribution of GABA(B) receptor homo- and hetero-dimers. Biochem J 388: 47–55.

98. Maurel D, Comps-Agrar L, Brock C, Rives ML, Bourrier E, et al. (2008) Cell-surface protein-protein interaction analysis with time-resolved FRET and snap-tag technologies: application to GPCR oligomerization. Nat Methods 5: 561–567.

99. Chang W, Tu C, Cheng Z, Rodriguez L, Chen TH, et al. (2007) Complex formation with the Type B gamma-aminobutyric acid receptor affects the expression and signal transduction of the extracellular calcium-sensing receptor. Studies with HEK-293 cells and neurons. J Biol Chem 282: 25030–25040.

100. Boyer SB, Clancy SM, Terunuma M, Revilla-Sanchez R, Thomas SM, et al. (2009) Direct interaction of GABAB receptors with M2 muscarinic receptors enhances muscarinic signaling. J Neurosci 29: 15796–15809.

101. Fernando RN, Eleuteri B, Abdelhady S, Nussenzweig A, Andang M, et al. (2011) Cell cycle restriction by histone H2AX limits proliferation of adult neural stem cells. Proc Natl Acad Sci U S A 108: 5837–5842.

102. Kriegstein AR (2005) GABA puts the brake on stem cells. Nat Neurosci 8: 1132–1133.
103. Fukui M, Nakamichi N, Yoneyama M, Ozawa S, Fujimori S, et al. (2008) Modulation of cellular proliferation and differentiation through GABA(B) receptors expressed by undifferentiated neural progenitor cells isolated from fetal mouse brain. *J Cell Physiol* 216: 507-519.

104. Oka M, Wada M, Wu Q, Yamamoto A, Fujita T (2006) Functional expression of metabotropic GABAB receptors in primary cultures of astrocytes from rat cerebral cortex. *Biochem Biophys Res Commun* 341: 874-881.

105. Mizuta K, Osawa Y, Mizuta F, Xu D, Emala CW (2008) Functional expression of GABAB receptors in airway epithelium. *Am J Respir Cell Mol Biol* 39: 296-304.

106. Tamura S, Watanabe M, Kanbara K, Yanagawa T, Watanabe K, et al. (2009) Expression and distribution of GABAergic system in rat knee joint synovial membrane. *Histol Histopathol* 24: 1009-1019.

107. Zhou Z, Sun H, Li X, Li Y, Zhao S, et al. (2010) A local GABAergic system is functionally expressed in human fallopian tube. *Biochem Biophys Res Commun* 398: 237-241.

108. Jiang X, Su L, Zhang Q, He C, Zhang Z, et al. (2012) GABAB Receptor Complex as a Potential Target for Tumor Therapy. *J Histochem Cytochem* 60: 269-279.

109. Roberts SS, Mendonca-Torres MC, Jensen K, Francis GL, Vasko V (2009) GABA receptor expression in benign and malignant thyroid tumors. *Pathol Oncol Res* 15: 645-650.

110. Wang T, Huang W, Chen F (2008) Baclofen, a GABAB receptor agonist, inhibits human hepatocellular carcinoma cell growth in vitro and in vivo. *Life Sci* 82: 536-541.

111. Tatsuta M, Ishi H, Baba M, Taniguchi H (1992) Attenuation by the GABA receptor agonist baclofen of experimental carcinogenesis in rat colon by azoxymethane. *Oncology* 49: 241-245.

112. Hannan S, Wilkins ME, Dehghani-Tafti E, Thomas P, Baddeley SM, et al. (2011) GABAB receptor internalisation is regulated by the R2 subunit. *J Biol Chem*.

113. Carmon KS, Gong X, Lin Q, Thomas A, Liu Q (2011) R-spondins function as ligands of the orphan receptors LGR4 and LGR5 to regulate Wnt/beta-catenin signaling. *Proc Natl Acad Sci U S A* 108: 11452-11457.

114. de Lau W, Barker N, Low TY, Koo BK, Li VS, et al. (2011) Lgr5 homologues associate with Wnt receptors and mediate R-spondin signalling. *Nature* 476: 293-297.

115. Birchmeier W (2011) Stem cells: Orphan receptors find a home. *Nature* 476: 287-288.

116. Sato T, van Es JH, Snippert HJ, Stange DE, Vries RG, et al. (2011) Paneth cells constitute the niche for Lgr5 stem cells in intestinal crypts. *Nature* 469: 415-418.

117. Snippert HJ, Clevers H (2011) Tracking adult stem cells. *EMBO Rep* 12: 113-122.

118. Barker N, van Es JH, Kuipers J, Kujala P, van den Born M, et al. (2007) Identification of stem cells in small intestine and colon by marker gene Lgr5. *Nature* 449: 1003-1007.

119. Mustata RC, Van Loy T, Lefort A, Libert F, Strollo S, et al. (2011) Lgr4 is required for Paneth cell differentiation and maintenance of intestinal stem cells ex vivo. *EMBO Rep* 12: 558-564.

120. Kato S, Matsubara M, Matsuo T, Mohri Y, Kazama I, et al. (2006) Leucine-rich repeat-containing G protein-coupled receptor-4 (LGR4, Gpr48) is essential for renal development in mice. *Nephron Exp Nephrol* 104: e63-75.

121. Oyama K, Mohri Y, Sone M, Nawa A, Nishimori K (2011) Conditional knockout of Lgr4 leads to impaired ductal elongation and branching morphogenesis in mouse mammary glands. *Sex Dev* 5: 205-212.

122. Wang Z, Jin C, Li H, Li C, Hou Q, et al. (2010) GPR48-Induced keratinocyte proliferation occurs through HB-EGF mediated EGFR transactivation. *FEBS Lett* 584: 4057-4062.

123. Gao Y, Kitagawa K, Hiramatsu Y, Kikuchi H, Isobe T, et al. (2006) Up-regulation of GPR48 induced by down-regulation of p27Kip1 enhances carcinoma cell invasiveness and metastasis. *Cancer Res* 66: 11623-11631.

124. Gao Y, Shan ZY, Wang H, Zhang HM, Teng WP (2009) Inhibitory effect of shRNA targeting GPR48 on invasion and metastasis of human cervical carcinoma cell line HeLa. *Ai Zheng* 28: 104-107.

125. Stella N (2010) Cannabinoid and cannabinoid-like receptors in microglia, astrocytes, and astrocytomas. *Glia* 58: 1017-1030.

126. Zhuang M, Yang Y, Cao F, Lu M, Wang X, et al. (2012) Associations of variants of CNR1 with obesity and obesity-related traits in Chinese women. *Gene* 495: 194-198.

127. Verde Z, Santiago C, Rodriguez Gonzalez-Moro JM, de Lucas Ramos P, Lopez Martin S, et al. (2011) 'Smoking genes': a genetic association study. *PLoS One* 6: e26668.

128. Smith SB, Maixner DW, Fillingim RB, Slade G, Gracely RH, et al. (2012) Large candidate gene association study reveals genetic risk factors and therapeutic targets for fibromyalgia. *Arthritis Rheum* 64: 584-593.

129. Carracedo A, Lorente M, Egia A, Blazquez C, Garcia S, et al. (2006) The stress-regulated protein p8 mediates cannabinoid-induced apoptosis of tumor cells. *Cancer Cell* 9: 301-312.

130. Salazar M, Carracedo A, Salanueva IJ, Hernandez-Tiedra S, Lorente M, et al. (2009) Cannabinoid action induces autophagy-mediated cell death through stimulation of ER stress in human glioma cells. *J Clin Invest* 119: 1359-1372.

131. Cudaback E, Marrs W, Moeller T, Stella N (2010) The expression level of CB1 and CB2 receptors determines their efficacy at inducing apoptosis in astrocytomas. *PLoS One* 5: e8702.

132. Aguado T, Carracedo A, Julien B, Velasco G, Milman G, et al. (2007) Cannabinoids induce glioma stem-like cell differentiation and inhibit gliomagenesis. *J Biol Chem* 282: 6854-6862.

133. Wu X, Han L, Zhang X, Li L, Jiang C, et al. (2012) Alteration of endocannabinoid system in human gliomas. *J Neurochem* 120: 842-849.

134. Held-Feindt J, Dorner L, Sahan G, Mehdorn HM, Mentlein R (2006) Cannabinoid receptors in human astroglial tumors. *J Neurochem* 98: 886-893.

135. Schley M, Stander S, Kerner J, Vajkoczy P, Schupfer G, et al. (2009) Predominant CB2 receptor expression in endothelial cells of glioblastoma in humans. *Brain Res Bull* 79: 333-337.

136. De Jesus ML, Hostalot C, Garibi JM, Salles J, Meana JJ, et al. (2010) Opposite changes in cannabinoid CB1 and CB2 receptor expression in human gliomas. *Neurochem Int* 56: 829-833.

137. Marshall AD, Lagutina I, Grosveld GC (2011) PAX3-FOXO1 induces cannabinoid receptor 1 to enhance cell invasion and metastasis. *Cancer Res* 71: 7471-7480.

138. Kenessey I, Banki B, Mark A, Varga N, Tovari J, et al. (2012) Revisiting CB1 Receptor as Drug Target in Human Melanoma. *Pathol Oncol Res*.

139. Proto MC, Gazzerro P, Di Croce L, Santoro A, Malfitano AM, et al. (2012) Interaction of endocannabinoid system and steroid hormones in the control of colon cancer cell growth. *J Cell Physiol* 227: 250-258.

140. Silva JP, Ushkaryov YA (2010) The latrophilins, “split-personality” receptors. *Adv Exp Med Biol* 706: 59-75.

141. Matsushita H, Lelianova VG, Ushkaryov YA (1999) The latrophilin family: multiply spliced G protein-coupled receptors with differential tissue distribution. *FEBS Lett* 443: 348-352.

142. Ichtchenko K, Bittner MA, Krasnoperov V, Little AR, Chepurny O, et al. (1999) A novel ubiquitously expressed alpha-latrotoxin receptor is a member of the CIRL family of G-protein-coupled receptors. *J Biol Chem* 274: 5491-5498.

143. Doyle SE, Scholz MJ, Greer KA, Hubbard AD, Darnell DK, et al. (2006) Latrophilin-2 is a novel component of the epithelial-mesenchymal transition within the atrioventricular canal of the embryonic chicken heart. *Dev Dyn* 235: 3213-3221.

144. Boucard AA, Ko J, Sudhof TC (2012) High Affinity Neurexin Binding to Cell Adhesion G-protein-coupled Receptor CIRL1/Latrophilin-1 Produces an Intercellular Adhesion Complex. *J Biol Chem* 287: 9399-9413.

145. Silva JP, Lelianova VG, Ermolyuk YS, Vysokov N, Hitchen PG, et al. (2011) Latrophilin 1 and its endogenous ligand Lasso/teneurin-2 form a high-affinity transsynaptic receptor pair with signaling capabilities. *Proc Natl Acad Sci U S A* 108: 12113-12118.

146. O'Sullivan ML, de Wit J, Savas JN, Comoletti D, Otto-Hitt S, et al. (2012) FLRT Proteins Are Endogenous Latrophilin Ligands and Regulate Excitatory Synapse Development. *Neuron* 73: 903-910.

147. White GR, Varley JM, Heighway J (1998) Isolation and characterization of a human homologue of the latrophilin gene from a region of 1p31.1 implicated in breast cancer. *Oncogene* 17: 3513-3519.

148. Eng L, Ibrahim-zada I, Jarjanazi H, Savas S, Meschian M, et al. (2011) Bioinformatic analyses identifies novel protein-coding pharmacogenomic markers associated with paclitaxel sensitivity in NCI60 cancer cell lines. *BMC Med Genomics* 4: 18.

149. Yang G, Quan Y, Wang W, Fu Q, Wu J, et al. (2012) Dynamic equilibrium between cancer stem cells and non-stem cancer cells in human SW620 and MCF-7 cancer cell populations. *Br J Cancer*.

150. Garritano S, Inga A, Gemignani F, Landi S (2013) More targets, more pathways and more clues for mutant p53. *Oncogenesis* 2: e54.

151. Del Valle J, Gantz I (1997) Novel insights into histamine H2 receptor biology. *Am J Physiol* 273: G987-996.

152. Hill SJ, Ganellin CR, Timmerman H, Schwartz JC, Shankley NP, et al. (1997) International Union of Pharmacology. XIII. Classification of histamine receptors. *Pharmacol Rev* 49: 253-278.

153. Medina VA, Rivera ES (2010) Histamine receptors and cancer pharmacology. *Br J Pharmacol* 161: 755-767.

154. Kubecova M, Kolostova K, Pinterova D, Kacprzak G, Bobek V (2011) Cimetidine: an anticancer drug? *Eur J Pharm Sci* 42: 439-444.

155. Zheng Y, Xu M, Li X, Jia J, Fan K, et al. (2013) Cimetidine suppresses lung tumor growth in mice through proapoptosis of myeloid-derived suppressor cells. *Mol Immunol* 54: 74-83.

156. Lisy K (2013) Histamine type 2 receptor antagonists as adjuvant treatment for resected colorectal cancer. *Clin J Oncol Nurs* 17: 211-212.

157. Deva S, Jameson M (2012) Histamine type 2 receptor antagonists as adjuvant treatment for resected colorectal cancer. *Cochrane Database Syst Rev* 8: CD007814.

158. Cricco G, Martin G, Medina V, Nunez M, Mohamad N, et al. (2006) Histamine inhibits cell proliferation and modulates the expression of Bcl-2 family proteins via the H2 receptor in human pancreatic cancer cells. *Anticancer Res* 26: 4443-4450.

159. McLatchie LM, Fraser NJ, Main MJ, Wise A, Brown J, et al. (1998) RAMPs regulate the transport and ligand specificity of the calcitonin-receptor-like receptor. *Nature* 393: 333-339.

160. Oehler MK, Norbury C, Hague S, Rees MC, Bicknell R (2001) Adrenomedullin inhibits hypoxic cell death by upregulation of Bcl-2 in endometrial
cancer cells: a possible promotion mechanism for tumour growth. Oncogene
20: 2937–2945.

161. Martinez A, Vos M, Guedez L, Kaur G, Chen Z, et al. (2002) The effects of adrenomedullin overexpression in breast tumor cells. J Natl Cancer Inst 94: 1226–1237.

162. Metellus P, Voutsinos-Porche B, Nanni-Metellus I, Colin C, Fina F, et al. (2011) Adrenomedullin expression and regulation in human glioblastoma, cultured human glioblastoma cell lines and pilocytic astrocytoma. Eur J Cancer 47: 1727–1735.

163. Kaafarani I, Fernandez-Sauze S, Berenguer C, Chinot O, Delfino C, et al. (2009) Targeting adrenomedullin receptors with systemic delivery of neutralizing antibodies inhibits tumor angiogenesis and suppresses growth of human tumor xenografts in mice. FASEB J 23: 3424–3435.

164. Pyne NJ, Tonelli F, Lim KG, Long JS, Edwards J, et al. (2012) Sphingosine 1-phosphate signalling in cancer. Biochem Soc Trans 40: 94–100.

165. Huang X, Wang Y, Nan X, He S, Xu X, et al. (2013) The role of the orphan G protein-coupled receptor 37 (GPR37) in multiple myeloma cells. Leuk Res.

166. Bhola NE, Grandis JR (2008) Crosstalk between G-protein-coupled receptors and epidermal growth factor receptor in cancer. Front Biosci 13: 1857–1865.

167. Daub H, Weiss FU, Wallasch C, Ullrich A (1996) Role of transactivation of the EGF receptor in signalling by G-protein-coupled receptors. Nature 379: 557–560.

168. Delcourt N, Bockaert J, Marin P (2007) GPCR-jacking: from a new route in RTK signalling to a new concept in GPCR activation. Trends Pharmacol Sci 28: 602–607.

169. Lappano R, De Marco P, De Francesco EM, Chimento A, Pezzi V, et al. (2013) Cross-talk between GPER and growth factor signaling. J Steroid Biochem Mol Biol.

170. Yu FX, Guan KL (2013) The Hippo pathway: regulators and regulations. Genes Dev 27: 355–371.

171. O'Hayre M, Vazquez-Prado J, Kufareva I, Stawiski EW, Handel TM, et al. (2013) The emerging mutational landscape of G proteins and G-protein-coupled receptors in cancer. Nat Rev Cancer 13: 412–424.

172. Kan Z, Jaiswal BS, Stinson J, Janakiraman V, Bhatt D, et al. (2010) Diverse somatic mutation patterns and pathway alterations in human cancers. Nature 466: 869–873.

173. Wootten D, Christopoulos A, Sexton PM (2013) Emerging paradigms in GPCR allostery: implications for drug discovery. Nat Rev Drug Discov 12: 630–644.

174. Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, et al. (2007) The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol 114: 97–109.

175. Varlet P, Soni D, Miquel C, Roux FX, Meder JF, et al. (2004) New variants of malignant glioneuronal tumors: a clinicopathological study of 40 cases. Neurosurgery 55: 1377-1391; discussion 1391-1372.

176. Galan-Moya EM, Le Guelte A, Lima Fernandes E, Thirant C, Dwyer J, et al. (2011) Secreted factors from brain endothelial cells maintain glioblastoma stem-like cell expansion through the mTOR pathway. EMBO Rep 12: 470–476.

177. Thirant C, Bessette B, Varlet P, Puget S, Cadusseau J, et al. (2011) Clinical relevance of tumor cells with stem-like properties in pediatric brain tumors. PLoS One 6: e16375.

178. Miguet L, Bechade G, Fornecker L, Zink E, Felden C, et al. (2009) Proteomic analysis of malignant B-cell derived microparticles reveals CD148 as a potentially useful antigenic biomarker for mantle cell lymphoma diagnosis. J Proteome Res 8: 3346–3354.

179. Geer LY, Markey SP, Kowalak JA, Wagner L, Xu M, et al. (2004) Open mass spectrometry search algorithm. J Proteome Res 3: 958–964.

180. Li C, Wong WH (2001) Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A 98: 31–36.