- Open Access
Intratumoural spatial distribution of S100B + folliculostellate cells is associated with proliferation and expression of FSH and ERα in gonadotroph tumours
Acta Neuropathologica Communications volume 10, Article number: 18 (2022)
Folliculostellate cells are S100B-expressing cells with numerous functions in the normal anterior pituitary. These cells have also been identified in pituitary neuroendocrine tumours (PitNETs), where their precise role remains elusive. Here, we aimed to build a refined cartography of S100B-expressing cells to characterise their interpatient and intratumoural spatial distribution, and to start identifying their potential functions in PitNETs. High-throughput histological analysis of S100B-stained tumour sections of 54 PitNETs revealed a significant decrease in S100B + cells in PitNETs compared to the normal anterior pituitary. A Ki67 index ≥ 3, a mitosis count > 2/10 per high power fields, and a proliferative status, were all associated with fewer S100B + cells in gonadotroph tumours. Gonadotroph tumours also showed interpatient and intratumoural heterogeneity in the spatial distribution of S100B + cells. The existence of an intratumoural heterogeneity was further confirmed by the incorporation to our spatial analysis of additional markers: Ki67, FSH, LH, ERα and SSTR2. The tumour areas with fewer S100B + cells displayed a higher percentage of Ki67 + cells, whereas strong positive correlations were observed between S100B + , FSH + , and ERα + cells. Such spatial associations suggest that S100B + folliculostellate cells could play a role in gonadotroph tumorigenesis, and may contribute to the maintenance of tumour cells in a low proliferating, FSH + /ERα + differentiated state. Albeit, further in-depth functional studies are required to decipher the mechanisms underlying these spatial associations and to potentially identify a therapeutic use.
Pituitary neuroendocrine tumours (PitNETs) are the second most common adult intracranial tumours. They account for around 15% of intracranial neoplasms and represent a major health issue with a prevalence of 1 in 1000 people [1,2,3]. Their diagnosis and the prediction of their behaviour are challenging due to differences in their cell-of-origin within the anterior pituitary gland (APG) (gonadotroph, corticotroph, and cells of the Pit1 lineage, namely somatotroph, lactotroph, and thyrotroph cells) and their secretory capacity (functioning or non-functioning lesions). Though patients suffering from Pit1 and corticotroph tumours have benefited from the development of effective therapeutic options, gonadotroph tumours, which originate from follicle-stimulating hormone (FSH)- and luteinizing hormone (LH)-producing cells, and account for up to 35% of all PitNETs [4, 5], currently lack effective medical treatments . In addition, gonadotroph tumours may present multiple recurrences with an unpredicted time-course following their surgical resection (relapse at 5 years being seen in 10–40% of patients) [7, 8]. Moreover, a small percentage of PitNETs displays an aggressive behaviour that remains difficult to predict [9, 10], characterised by resistance and/or multiple recurrences despite multimodal treatments [9, 11]. Therefore, the identification of factors associated with and/or leading to an aggressive phenotype, as well as the identification of new therapeutic targets and novel treatment options are of particular importance.
The tumour microenvironment (TME) that includes blood and lymph vessels, immune cells, cancer-associated fibroblasts, as well as extracellular matrix components, and soluble molecules , represents an appealing prognostic tool and a therapeutic target, and has been extensively described in numerous cancers . However, the knowledge on the TME of PitNETs is still limited [14, 15]. Folliculostellate (FS) cells are non-endocrine APG-resident cells, identified based on their morphology and immunoreactivity to the S100B protein [16, 17]. FS cells account for ~ 5% of the normal APG [18, 19], where they have numerous roles, including a supportive role , the regulation of hormone secretion , phagocytic properties , and the production of numerous cytokines and growth factors such as interleukin 6 (IL-6), follistatin, basic fibroblast growth factor, transforming growth factor β, and vascular endothelial growth factor . FS cells were also described in PitNETs [23,24,25], suggesting they might play a role in these tumours. However, their function on tumorigenesis-related processes and their usefulness as a biomarker or as a therapeutic target remain underexplored.
Here, we assessed the presence and distribution of S100B + cells within the TME of PitNETs in order to provide a refined cartography of S100B + cells within tumours and between patients. We developed a high-throughput pipeline of single-cell histological analysis based on whole tumour serial sections obtained from a series of 54 surgically-resected PitNETs. We uncover an interpatient and an intratumoural heterogenous distribution of S100B + cells in gonadotroph tumours. In addition, our findings highlight the spatial relationships between the presence of FS cells and the proliferative capacity, as well as FSH and oestrogen receptor alpha (ERα) immunoreactivity of gonadotroph tumour cells, suggesting potential roles of S100B + cells in gonadotroph tumorigenesis.
Material and methods
Tumour and normal APG samples
All pituitary tumours (n = 54) were collected by the same surgeon (E.J.), and subsequently diagnosed in the Pathology Department of Groupement Hospitalier Est, Hospices Civils de Lyon, Lyon, France by a single pathologist (A.V.). Normal APG samples (n = 4) used as controls were obtained from the Pathology Department of Targu Mures Emergency Hospital, Targu Mures, Romania. These were collected during autopsies performed at 6–20 h after death on adult patients who have died suddenly, without prior known endocrine dysfunction. Those patients did not have traumatic head injury, documented prolonged multiple organ failure or cerebral ischemia prior to death, and they did not have evident signs of cerebral autolysis or incidental sellar/parasellar lesions.
S100B immunostaining was prospectively incorporated in the histopathological and immunohistochemistry (IHC) diagnostic routines. Clinicopathological data were collected from the medical files of patients included. The IHC diagnosis of different PitNET subtypes was based on hormone immunoreactivity: FSH, LH, prolactin, growth hormone, thyroid-stimulating hormone, and adrenocorticotropic hormone. When no hormone expression was identified, IHC for transcription factors (steroidogenic factor 1 (SF1), pituitary specific transcription factor (Pit1), and T-box transcription factor TBX19 (T-pit)), completed the diagnosis [26, 27]. For statistical analysis, tumours of the Pit1 lineage were considered together. Tumours were classified as functioning or non-functioning based on FSH, LH, and sex hormone levels; prolactin levels; growth hormone levels during the glucose tolerance test and insulin-like growth factor 1; thyroid-stimulating hormone and free thyroxine levels; free urinary cortisol, the response to dexamethasone suppression tests, and adrenocorticotropic hormone levels. Tumours were classified as proliferative (i.e., “proliferative status”) when at least two of three criteria of proliferation, namely Ki67 index ≥ 3%, mitoses n > 2/10 high power fields (HPFs) (400 × magnification), and positive p53 (> 10 strongly positive nuclei/10 HPFs), were present . Invasion was defined as cavernous and/or sphenoid sinus invasion detected on magnetic resonance imaging and/or during surgery; for the cavernous sinus invasion detected on magnetic resonance imaging, the revised Knosp classification was used . The grading of tumours was based on the published five-tiered clinicopathological classification (1a, non-invasive and non-proliferative; 1b, non-invasive, but proliferative; 2a, invasive, but not proliferative; and 2b, invasive and proliferative; 3, metastatic tumour) .
For IHC, staining was conducted on 4 µm serial sections obtained from formalin-fixed paraffin-embedded (FFPE) tissues through a routine protocol performed on a BenchMark® ULTRA immunostainer (Ventana Medical Systems Inc, Roche Diagnostics, Tucson, Arizona, USA). Detection of positive staining was done using diaminobenzidine (ultraView Universal DAB Detection Kit, Ventana Medical Systems Inc, Tucson, Arizona, USA) followed by haematoxylin counterstaining (Bluing Reagent, Ventana Medical Systems Inc, Roche Diagnostics, Tucson, Arizona, USA). Primary antibodies for IHC were as follow: S100B (Agilent Dako, #Z0311), FSH (Beckman Coulter Immunotech, #0373), LH (Agilent Dako, #M3502), Ki67 (Agilent Dako, #M7240), ERα (clone SP1, Ventana Medical Systems Inc.), and somatostatin receptor type 2 (SSTR2) (Abcam, #ab134152). For immunofluorescence (IF), staining was performed using standard protocols. Briefly, 4 µm-thick FFPE sections were first deparaffinised at 55 °C for 30 min prior to 5 min of incubation in xylene. Tissue sections were subsequently rehydrated in decreasing concentrations of ethanol (100%, 95% and 70%). Heat-induced epitope retrieval was performed using a citrate-based unmasking solution at pH 6.0 (Vector Laboratories, California, USA, #H-3300) prior to a 1 h blocking step and an overnight incubation with S100B (Abcam, #ab52642) and FSH (Beckman Coulter Immunotech, #0373) primary antibodies. Secondary antibodies coupled to fluorochromes (anti-mouse-AlexaFluor488, cat#712-546-151, and anti-rabbit-Cy3, cat#711–166-152, Jackson ImmunoResearch, UK) were then incubated for 30 min before mounting the sections with a DAPI counterstaining solution (Vector Laboratories, California, USA, #H-1200).
Slide scanning and image analysis
Stained whole sections were scanned at 20 × using a Pannoramic Scan (3D HISTECH Ltd., Budapest, Hungary) or a Zeiss Axio Scan.Z1 (Carl Zeiss Microscopy GmbH, Jena, Germany). Scanned images were visualised with CaseViewer software (3D HISTECH Ltd., Budapest, Hungary) and Zen2 software (blue edition). Image analysis, subsequent segmentation, and staining quantification, were performed with HALO® software v.2.3.2089.27 (Indica Labs, New Mexico, USA). Our pathologist (A.V.) manually annotated the areas corresponding to tumour tissue and to adjacent APG tissue (when present). The regions with staining artefacts were excluded prior to performing a single-cell segmentation with HALO® CytoNuclear module v1.6. Segmentation was then double-checked by eye for each sample to ensure that cells were correctly identified before quantifying the percentage of positive cells out of the total number of cells present within the area analysed. Finally, spatial plots that illustrate each stained cell as an individual point, were generated.
Statistical analysis was performed in R commander/R package version 3.5.1 (R-project, Vienna, Austria) and GraphPad Prism 5 (GraphPad, California, USA). Mann–Whitney test, Kruskal–Wallis test followed by Dunn’s post-test, Wilcoxon matched paired test, Student's t-test with Welch’s correction, and Spearman’s correlation were applied as indicated. p < 0.05 (adjusted p < 0.05 in case of multiple comparisons) was considered significant. Graphs were generated with GraphPad Prism 5 (GraphPad, California, USA).
High-throughput histological single-cell analysis based on whole-slide quantification reveals a loss of S100B + cells in PitNETs in comparison to normal APG
Quantification of immunoreactive cells on tissue sections is generally conducted either semi-quantitatively or by quantifying a limited number of selected fields. For the purpose of our study, we deemed such approaches to be imprecise, as the heterogeneous distribution of cells among tissues analysed can be overlooked. Thus, we developed a pipeline of IHC analysis to quantify single-cell staining on whole tumour sections (Fig. 1A). Analysis of FSH and SF1 expression of a gonadotroph tumour confirmed the feasibility of our technical approach, giving us the possibility to quantify and address the spatial distribution of immunoreactive cells in a large number (hundreds of thousands) of cells per section (Fig. 1B, C).
Following this validation, we analysed the number and the spatial distribution of FS cells through their S100B immunoreactivity, in a cohort of 54 patients with PitNETs (26 gonadotroph tumours, 20 tumours of the Pit1 lineage, and 8 corticotroph tumours) and 4 normal APG samples. The clinicopathological characteristics of the tumours analysed are detailed in Table 1. Quantification of S100B + cells on scanned whole-sections (Additional file 1: Fig. S1) confirmed that S100B + cells could be detected in normal APG, as well as in gonadotroph, corticotroph, and Pit1 tumours (Fig. 2A). Quantification of S100B + cells at the whole tumour level highlighted that tumour tissue showed a significant decrease in S100B + cells (median 0.18%) compared to normal APG (median 3.95%) (Fig. 2B). To refine this observation, we used histological sections in which both tumour and adjacent APG areas were present, and comparatively quantified the percentage of S100B + cells in individual patients. As shown in Fig. 2C, this quantification confirmed that the percentage of S100B + cells was lower in the tumour tissue compared to adjacent APG tissue, independently of the PitNET subtype. Interestingly, whereas most of the tumours analysed displayed a very low percentage of S100B + cells (43 tumours had < 1% of S100B + cells and 19 tumours < 0.1% S100B + cells), a small fraction of them had a higher percentage of S100B + cells reaching a maximum of 9.55% (Fig. 2B). Taken together, these observations suggest an interpatient heterogeneity regarding the percentage of intratumoral S100B + cells, which, nonetheless, is decreased in the majority of PitNETs in comparison to the normal APG.
The presence of intratumoural S100B + cells is associated with lower proliferative properties and positively correlated with ERα expression in gonadotroph tumours
We next investigated whether the percentage of intratumoural S100B + cells at the whole tumour level was associated with the clinicopathological characteristics of PitNETs analysed. We found that tumours with a Ki67 index ≥ 3% were associated with significantly fewer S100B + cells (p = 0.01) (Fig. 3A), whereas there was no association between the percentage of S100B + cells and the tumour invasive trait (Fig. 3B). We then wondered whether these observations were related to specific PitNET subtypes. Although the percentage of intratumoural S100B + cells was not statistically associated with the tumour subtype itself (Fig. 3C), a Ki67 index ≥ 3%, a mitosis count > 2/10 HPFs, and a proliferative status, were all significantly associated with a lower percentage of S100B + cells exclusively in gonadotroph tumours (Fig. 3D–F), but not in Pit1, nor in corticotroph tumours (Additional file 2: Fig. S2A-C and Fig. S2E-G). Moreover, S100B + cells were more widely spread in gonadotroph tumours (range 0.01–9.55%) compared to Pit1 (range 0–2.17%) and corticotroph tumours (range 0.007–3.42%). No association between invasiveness and the percentage of intratumoural S100B + cells was obtained when considering PitNET subtypes separately (Fig. 3G, Additional file 2: Fig. S2D and S2H). In addition, no statistical significance between the percentage of S100B + cells and other clinicopathological characteristics (age, gender, tumour dimension, grade, and the percentage of SSTR2 + or SSTR5 + cells) was reached in gonadotroph, corticotroph, or Pit1 tumours (Fig. 3H, I and data not shown); of note, SSTR5 was analysed only for corticotroph and Pit1 tumours. Interestingly, the percentage of S100B + cells showed a moderate positive correlation with the percentage of ERα + cells in gonadotroph tumours (Fig. 3J). These observations suggest that intratumoural S100B + cells may play a role in gonadotroph tumorigenesis-related processes.
Spatial analysis of S100B + cell distribution highlights an intratumoural heterogeneity of S100B + cells
Given that our results indicated a possible role for S100B + cells in gonadotroph tumorigenesis, we aimed at refining our cartography of this rare population of TME cells, by incorporating a spatial resolution to our analysis. In order to investigate the intratumoural distribution of S100B + cells in our cohort of 26 gonadotroph tumours, we defined for each sample multiple tumour areas on the same tumour slide. The percentage of S100B + cells was then quantified in each area, and the results subsequently confirmed by visually inspecting the stained areas (Fig. 4A). Analysis of the 26 gonadotroph tumours confirmed interpatient heterogeneity, with a subset of tumours displaying almost no S100B expression within the areas quantified and a group of lesions presenting areas with variable percentages of S100B + cells (Fig. 4A, B). Among the latter, the percentage of S100B + cells was heterogeneous between tumours, but also within the different areas of the same tumour, highlighting an intratumoural spatial heterogeneity of S100B + cells in gonadotroph PitNETs (Fig. 4).
The intratumoural spatial distribution of S100B + cells is associated with FSH, ERα and Ki67 immunoreactivity
The identification of a heterogeneous spatial distribution of S100B + cells in individual gonadotroph tumours led us to question whether the spatial distribution of S100B + cells may be of functional relevance in these tumours. To address this question, we built a cartography based on the IHC analysis of multiple markers on serial tumour sections, including five new markers, namely Ki67 (a proliferation marker), FSH and LH (the two hormones produced by gonadotroph cells), SSTR2 (a potential marker of response to treatment), and ERα (a transcription factor and a potential marker of response to treatment) [26, 30] (Fig. 5A). Spatial quantification of each marker was performed on 20 independent areas of 5 gonadotroph tumours detailed in Table 2 (4 areas per tumour). We chose the 5 gonadotroph tumours based on the percentages of S100B + cells: 4 tumours having a heterogeneous distribution, including the tumour having the widest spread of the percentages of S100B + cells—gonadotroph tumour #3, and one tumour having homogeneous percentages of S100B + cells—gonadotroph tumour #18.
The spatial quantification of each marker shown in Fig. 5B revealed their heterogeneous distribution between the different tumours. In these 5 tumours, the percentage of ERα-expressing and LH-expressing cells seemed to show less intratumoural heterogeneity than S100B-expressing, FSH-expressing, Ki67-expressing, and SSTR2-expressing cells (Fig. 5B). For each area, we further questioned the spatial association of the tested markers with the percentage of S100B + cells. Areas expressing low and high levels of S100B were defined based on the median percentage of S100B + cells within the 20 areas (n = 10 for each group). S100BHigh areas had significantly more FSH + and ERα + cells (p = 0.0001, and p = 0.008, respectively) (Fig. 5C, D, Fig. 6A, Additional file 3: Fig. S3A) and fewer Ki67 + cells (p = 0.02) (Fig. 5E, Additional file 3: Fig. S3B). The percentage of LH + and SSTR2 + cells was not significantly different between S100BHigh and S100BLow tumour areas (Additional file 4: Fig. S4). Correlation analysis of the 20 tumour areas further confirmed that the percentage of S100B + cells was positively and strongly correlated with the percentage of FSH + cells (Spearman’s rho = 0.78, adjusted p = 0.0007) and of ERα + cells (Spearman’s rho = 0.76, adjusted p = 0.001) (Table 3). A strong positive correlation was also found between the percentage of FSH + cells and ERα + cells (Spearman’s rho = 0.67, adjusted p = 0.01) (Table 3), suggesting the existence of tumour areas simultaneously rich and poor in FSH + , ERα + , and S100B + cells (Additional file 3: Fig. S3C).
To further validate these associations, we confirmed the existence of gonadotroph tumour areas presenting either an enrichment in S100B, ERα, and FSH immunoreactivity, or a lack of immunoreactivity for the three markers (Fig. 6B). Double IF S100B/FSH staining also revealed that S100B + cells were mainly located in areas enriched in FSH immunoreactivity (Fig. 6C). Hence, these observations confirmed the heterogeneous nature of gonadotroph tumours, and suggest that the heterogeneous distribution of S100B + cells might be of functional relevance within the TME of gonadotroph PitNETs.
The concepts of intratumoural, intertumoural and interpatient heterogeneity are intensively studied and acknowledged to be of tremendous importance in cancer biology [31, 32]. In PitNETs, tumour heterogeneity has only started to be unravelled by recent studies [30, 33, 34]. Here, our work contributes to extending this knowledge by providing an insight into the heterogeneity of PitNETs through the cartography of a rare TME component, namely S100B-expressing FS cells. Our work presents a novel high-throughput histological approach to analyse the spatial distribution and the association of multiple markers on scanned whole sections. To our knowledge, this report is the first to investigate the intratumoural spatial distribution of FS cells, and even gonadotroph intratumoural heterogeneity in general. As illustrated by our results, this kind of approach may refine our knowledge of tumour heterogeneity and favour the identification of functional associations between TME components and tumour cells.
Using whole-scan quantification of S100B IHC staining, we first confirmed an overall loss of S100B + FS cells in the majority of our PitNETs compared to normal APG and to adjacent APG, with 80% of these tumours having < 1% S100B + cells, and 35% having < 0.1% S100B + cells. This result is consistent with previous work that reported a loss of S100B + cells in PitNETs in comparison to normal APG collected from autopsies or from APG regions adjacent to tumour tissue [19, 35]. Interestingly, our analysis further revealed an interpatient heterogeneity with percentages of S100B + cells ranging from 0 to 9.55% in the 54 PitNETs analysed. Among PitNET subtypes, gonadotroph tumours had the most variable S100B distribution, though this may be partly due to our cohort composition.
Retrospective investigation of clinicopathological characteristics revealed that PitNETs with a Ki67 index ≥ 3% were associated with fewer S100B + cells. Further analysis comparing PitNET subtypes, pointed this proliferative trait to be driven by the gonadotroph subtype. Indeed, in gonadotroph tumours, not only was a Ki67 index ≥ 3% significantly associated with fewer S100B + cells, but also a mitosis count > 2/10 per HPFs, and the proliferative status. While we also found the overall percentage of S100B + cells to positively correlate with the percentage of ERα + cells in gonadotroph tumours, the overall percentage of S100B + cells was not associated with age, gender, invasion or tumour dimension in any of the PitNET subtypes, a result consistent with previous work by Voit et al. in acromegaly patients . Moreover, the percentage of intratumoural S100B + cells was not statistically associated with the grade, and did not correlate with the percentage of SSTR2 + cells in any of the PitNET subtypes. In Pit1 and in corticotroph tumours, this percentage was not associated with proliferation markers, and did not correlate with the percentage of SSTR5 + cells, either. Of note, in gonadotroph tumours, it would have been interesting to also compare S100B + cells to SSTR3 + cells. Although the absence of statistically significant associations might be due to small sample sizes, these observations support nevertheless the hypothesis that S100B + cells may have distinct functions in different PitNET subtypes, and that the loss of these cells could favour the proliferative status of gonadotroph tumours.
Spatial analysis of cell distribution in different tumour areas of 26 gonadotroph PitNETs further revealed an intratumoural heterogeneity of S100B + cells. While we identified a first group of gonadotroph tumours with almost no S100B + cells, we found a second group that presented a wide distribution of S100B + cells ranging from 0.005 to 18.85% within different areas analysed. The presence of an intratumoural heterogeneity and the existence of low- and high-expressing S100B areas was intriguing and made us question the functional consequences of such a spatial cell distribution in gonadotroph tumours. To address this question, we added to our spatial analysis 5 new markers: Ki67, FSH, LH, ERα and SSTR2. Analysing 20 areas from 5 independent gonadotroph tumours, we found that higher percentages of S100B + cells were associated with lower percentages of Ki67. This observation was consistent with our previous observation based on the overall percentage of S100B + cells per tumour, and it further confirmed the functional relationship that may exist between gonadotroph tumour proliferation and the loss of S100B + cells.
More interestingly, through our spatial analysis, we observed strong positive correlations between the percentage of S100B + cells and FSH + cells, and between S100B + cells and ERα + cells, respectively. A similar correlation was also found between the percentage of FSH + cells and ERα + cells, suggesting the occurrence of tumour areas simultaneously rich and poor in FSH + , ERα + , and S100B + cells. Consistently, a positive correlation between FSH and ERα at the whole tumour level has already been reported in gonadotroph tumours [30, 36], suggesting that our technical approach is appropriate to investigate spatial associations of marker co-expression. Given the fact that ERα is also expressed in the vast majority of normal human gonadotroph cells which produce, secrete and stain for FSH [37, 38], we hypothesize that gonadotroph tumour cells expressing ERα and FSH might be closer to a normal phenotypic state than gonadotroph tumour cells lacking ERα and FSH expression. Identification of an enrichment in S100B + folliculostellate cells in areas also enriched in ERα-expressing and FSH-expressing cells may suggest a role for these cells in maintaining a more differentiated, functioning or closer to the normal phenotype state of gonadotroph tumour cells. This hypothesis appears to be coherent with our other findings, as the loss of intratumoural S100B + cells was associated with increased proliferation. Alternatively, FS cells may secrete cofactors that modulate the signalling of oestrogens through ERα  leading to FSH production, or FS cells may secrete molecules or physically interact with gonadotroph tumour cells leading to FSH production independently of ERα. Such molecules might be IL-6 or transforming growth factor β, both of which were shown to be produced by FS cells [22, 40], and to stimulate the secretion of FSH by the normal APG [40, 41]. Indeed, given that FS cells have been shown to modulate the hormone secretion of the normal APG  and that relationships between S100B + cells and the preoperative levels of growth hormone and of prolactin were revealed in acromegaly patients , it would not be surprising if FS cells also have an effect on hormone production/secretion by gonadotroph tumour cells.
Knowing that FS cells produce a large numbers of secreted cytokines, it would be tempting to speculate that FS cells could modulate tumour cell behaviour through paracrine signalling, by acting, for example, as a senescence inducer and/or proliferation controller. Interestingly, the secretion of the senescence-associated IL-6 cytokine by FS cells was shown to inhibit the proliferation of endocrine cells in the normal APG . In pituitary tumours, it has been shown that IL-6 was mainly produced by tumour cells [40, 42], and that its silencing in primary cultures of human pituitary tumours decreased tumour senescence . To the same extent, FS cells also produce follistatin , which may be another interesting candidate since dysregulation of follistatin function was shown to be linked to multiple cancers, including pituitary adenomas . Our results highlight therefore the need for refining our knowledge on the secretome of intratumoural S100B + cells since it could have potent actions in PitNETs.
Aside from their possible function in gonadotroph tumours and other PitNETs, one may question whether FS cells may act as “gatekeeper” cells involved in maintaining normal structure and function of the APG, and whether their loss might favour the onset of pituitary tumours. Studies on the normal human APG have shown an increased presence of S100B + cells with age, especially after the age of 70 [45, 46], and proposed that the increased presence of these cells represents the successful modality of APG aging . Given the multiple roles FS cells have in the normal APG through paracrine [17, 22] and physical actions [47, 48], it appears indeed possible that these cells are necessary to maintain a normal tissue structure and function. Their loss might therefore shift the balance from normal tissue preservation to adenoma formation/progression. And, reciprocally, a normal APG tissue might be required to maintain the presence of FS cells.
We conducted an interpatient and an intratumoural cartography of S100B + cells to gain insight into the potential functions these cells might have in gonadotroph tumours. Concomitantly, we demonstrated the feasibility of analysing whole tumour sections at the single-cell resolution in human PitNETs, and of spatially analysing multiple markers on serial sections. The overall percentage of S100B + cells was associated with lower proliferative properties and positively correlated with ERα expression in gonadotroph tumours. However, the overall percentage of S100B + cells was not associated/correlated with age, gender, invasion, grade, tumour dimension and percentage of SSTR2 + cells in any of the PitNET subtypes. In Pit1 and in corticotroph tumours, this percentage was not associated with proliferation markers, and did not correlate with the percentage of SSTR5 + cells, either. Our work highlights a potential role for intratumoural S100B + cells in modulating the proliferation, the differentiation, and the hormonal production of gonadotroph tumour cells. It also suggests that S100B + cells may contribute to maintaining a phenotype closer to the normal tissue in certain tumour areas. However, the importance of S100B + cell loss in pituitary tumorigenesis remains to be demonstrated. Our work paves the way towards a refined view of gonadotroph tumour heterogeneity and towards a better understanding of the role of S100B + folliculostellate cells in gonadotroph tumours. It also emphasizes the need to perform high-resolution analysis of the TME at the single-cell level to unravel novel mechanisms that drive tumorigenesis-associated processes, which will hopefully lead to the discovery of new diagnostic and prognostic markers and of new and/or personalized treatments.
Availability of data and materials
Anterior pituitary gland
Formalin-fixed paraffin embedded
High power fields
Oestrogen receptor alpha
Pituitary neuroendocrine tumours
Somatostatin receptor type 2
Steroidogenic factor 1
Asa SL, Casar-Borota O, Chanson P, Delgrange E, Earls P, Ezzat S et al (2017) From pituitary adenoma to pituitary neuroendocrine tumor (PitNET): an International Pituitary Pathology Club proposal. Endocr Relat Cancer 24:C5-8
Melmed S (2020) Pituitary-tumor endocrinopathies. N Engl J Med 382:937–950
Raverot G, Ilie MD, Lasolle H, Amodru V, Trouillas J, Castinetti F et al (2021) Aggressive pituitary tumours and pituitary carcinomas. Nat Rev Endocrinol 17:671–684
Nishioka H, Inoshita N, Mete O, Asa SL, Hayashi K, Takeshita A et al (2015) The complementary role of transcription factors in the accurate diagnosis of clinically nonfunctioning pituitary adenomas. Endocr Pathol 26:349–355
Drummond J, Roncaroli F, Grossman AB, Korbonits M (2019) Clinical and pathological aspects of silent pituitary adenomas. J Clin Endocrinol Metab 104:2473–2489
Ilie MD, Raverot G (2020) Treatment options for gonadotroph tumors: current state and perspectives. J Clin Endocrinol Metab 105:e3507–e3518
Raverot G, Dantony E, Beauvy J, Vasiljevic A, Mikolasek S, Borson-Chazot F et al (2017) Risk of recurrence in pituitary neuroendocrine tumors: a prospective study using a five-tiered classification. J Clin Endocrinol Metab 102:3368–3374
Raverot G, Vasiljevic A, Jouanneau E (2018) Prognostic factors of regrowth in nonfunctioning pituitary tumors. Pituitary 21:176–182
Raverot G, Burman P, McCormack A, Heaney A, Petersenn S, Popovic V et al (2018) European Society of Endocrinology Clinical Practice Guidelines for the management of aggressive pituitary tumours and carcinomas. Eur J Endocrinol 178:G1-24
Asa SL, Ezzat S (2016) Aggressive pituitary tumors or localized pituitary carcinomas: defining pituitary tumors. Exp Rev Endocrinol Metab 11:149–162
McCormack A, Dekkers OM, Petersenn S, Popovic V, Trouillas J, Raverot G et al (2018) Treatment of aggressive pituitary tumours and carcinomas: results of a European Society of Endocrinology (ESE) survey 2016. Eur J Endocrinol 178:265–276
Balkwill FR, Capasso M, Hagemann T (2012) The tumor microenvironment at a glance. J Cell Sci 125:5591–5596
Hui L, Chen Y (2015) Tumor microenvironment: sanctuary of the devil. Cancer Lett 368:7–13
Ilie MD, Vasiljevic A, Raverot G, Bertolino P (2019) The microenvironment of pituitary tumors—biological and therapeutic implications. Cancers 11:1605
Marques P, Grossman AB, Korbonits M. The tumour microenvironment of pituitary neuroendocrine tumours. Front Neuroendocrinol. 2020;58:100852.
Nakajima T, Yamaguchi H, Takahashi K (1980) S100 protein in folliculostellate cells of the rat pituitary anterior lobe. Brain Res 191:523–531
Morris J, Christian H (2011) Folliculo-stellate cells: paracrine communicators in the anterior pituitary. Open Neuroendocrinol J 4:77–89
Allaerts W, Vankelecom H (2005) History and perspectives of pituitary folliculo-stellate cell research. Eur J Endocrinol 153:1–12
Hofler H, Walter GF, Denk H (1984) Immunohistochemistry of folliculo-stellate cells in normal human adenohypophyses and in pituitary adenomas. Acta Neuropathol (Berl) 65:35–40
Delfin L, Mete O, Asa SL (2021) Follicular cells in pituitary neuroendocrine tumors. Hum Pathol 114:1–8
Devnath S, Inoue K (2008) An insight to pituitary folliculo-stellate cells. J Neuroendocrinol 20:687–691
Denef C (2007) Paracrinicity: the story of thirty years of cellular pituitary crosstalk. J Neuroendocrinol 20:1–70
Vajtai I, Kappeler A, Sahli R (2007) Folliculo-stellate cells of “true dendritic” type are involved in the inflammatory microenvironment of tumor immunosurveillance of pituitary adenomas. Diagn Pathol 2(20):1–4
Farnoud M, Kujas M, Derome P, Racadot J, Peillon F, Li J (1994) Interactions between normal and tumoral tissues at the boundary of human anterior pituitary adenomas. An immunohistochemical study. Virchows Arch 424:75–82
Voit D, Saeger W, Lüdecke DK (1999) Folliculo-stellate cells in pituitary adenomas of patients with acromegaly. Pathol Res Pract 195:143–147
Osamura RY, Grossman A, Korbonits M, Kovacs K, Lopes MBS, Matsuno A, et al. Pituitary adenoma. WHO classif tumours endocr organs chapter 1 tumors Pituit Gland. Lyon: IARC Press; 2017. p. 14–8.
Villa C, Vasiljevic A, Jaffrain-Rea ML, Ansorge O, Asioli S, Barresi V et al (2019) A standardised diagnostic approach to pituitary neuroendocrine tumours (PitNETs): a European Pituitary Pathology Group (EPPG) proposal. Virchows Arch 475:687–692
Trouillas J, Roy P, Sturm N, Dantony E, Viennet G, Bonneville J-F et al (2013) A new prognostic clinicopathological classification of pituitary adenomas: a multicentric case–control study of 410 patients with 8 years post-operative follow-up. Acta Neuropathol (Berl) 126:123–135
Micko ASG, Wöhrer A, Wolfsberger S, Knosp E (2015) Invasion of the cavernous sinus space in pituitary adenomas: endoscopic verification and its correlation with an MRI-based classification. J Neurosurg 122:803–811
Ilie MD, Vasiljevic A, Louvet C, Jouanneau E, Raverot G (2020) Gonadotroph tumors show subtype differences that might have implications for therapy. Cancers 12:1012
Grzywa TM, Paskal W, Włodarski PK (2017) Intratumor and intertumor heterogeneity in melanoma. Transl Oncol 10:956–975
Liu J, Dang H, Wang XW (2018) The significance of intertumor and intratumor heterogeneity in liver cancer. Exp Mol Med 50:e416–e416
Hickman RA, Bruce JN, Otten M, Khandji AG, Flowers XE, Siegelin M et al (2021) Gonadotroph tumours with a low SF-1 labelling index are more likely to recur and are associated with enrichment of the PI3K-AKT pathway. Neuropathol Appl Neurobiol 47:415–427
Neou M, Villa C, Armignacco R, Jouinot A, Raffin-Sanson M-L, Septier A et al (2020) Pangenomic classification of pituitary neuroendocrine tumors. Cancer Cell 37:123-134.e5
Iwaki T, Kondo A, Takeshita I, Nakagaki H, Kitamura K, Tateishi J (1986) Proliferating potential of folliculo-stellate cells in human pituitary adenomas: immunohistochemical and electron microscopic analysis. Acta Neuropathol (Berl) 71:233–242
Kolnes AJ, Øystese KAB, Olarescu NC, Ringstad G, Berg-Johnsen J, Casar-Borota O et al (2020) FSH levels are related to E-cadherin expression and subcellular location in nonfunctioning pituitary tumors. J Clin Endocrinol Metab 105:2587–2594
Friend KE, Chiou YK, Lopes MB, Laws ER, Hughes KM, Shupnik MA (1994) Estrogen receptor expression in human pituitary: correlation with immunohistochemistry in normal tissue, and immunohistochemistry and morphology in macroadenomas. J Clin Endocrinol Metab 78:1497–1504
Shupnik MA, Pitt LK, Soh AY, Anderson A, Lopes MB, Laws ER (1998) Selective expression of estrogen receptor α and β isoforms in human pituitary tumors 1. J Clin Endocrinol Metab 83:3965–3972
Shang Y, Hu X, DiRenzo J, Lazar MA, Brown M (2000) Cofactor dynamics and sufficiency in estrogen receptor-regulated transcription. Cell 103:843–852
Renner U, Gloddek J, Pereda MP, Arzt E, Stalla GK (1998) Regulation and role of intrapituitary IL-6 production by folliculostellate cells. Domest Anim Endocrinol 15:353–362
Ying S-Y, Becker A, Baird A, Ling N, Ueno N, Esch F et al (1986) Type beta transforming growth factor (TGF-β) is a potent stimulator of the basal secretion of follicle stimulating hormone (FSH) in a pituitary monolayer system. Biochem Biophys Res Commun 135:950–956
Ueta Y, Levy A, Chowdrey HS, Lightman SL (1995) S?100 antigen-positive folliculostellate cells are not the source of IL-6 gene expression in human pituitary adenomas. J Neuroendocrinol 7:467–474
Sapochnik M, Haedo MR, Fuertes M, Ajler P, Carrizo G, Cervio A et al (2017) Autocrine IL-6 mediates pituitary tumor senescence. Oncotarget 8:4690–4702
Bilezikjian LM, Justice NJ, Blackler AN, Wiater E, Vale WW (2012) Cell-type specific modulation of pituitary cells by activin, inhibin and follistatin. Mol Cell Endocrinol 359:43–52
Pavlović M, Jovanović I, Ugrenović S, Vasović L, Krstić M, Bakić M et al (2013) Morphometric analysis of the human anterior pituitary’s folliculostellate cells during the aging process. Ann Anat Anat Anz 195:231–237
Čukuranović Kokoris J, Jovanović I, Pantović V, Krstić M, Stanojković M, Milošević V et al (2017) Morphometric analysis of the folliculostellate cells and luteinizing hormone gonadotropic cells of the anterior pituitary of the men during the aging process. Tissue Cell 49:78–85
Fauquier T, Guerineau NC, McKinney RA, Bauer K, Mollard P (2001) Folliculostellate cell network: a route for long-distance communication in the anterior pituitary. Proc Natl Acad Sci 98:8891–8896
Hodson DJ, Legros C, Desarménien MG, Guérineau NC (2015) Roles of connexins and pannexins in (neuro)endocrine physiology. Cell Mol Life Sci 72:2911–2928
The authors thank CRB HCL—Tissue Tumorothèque Est for its assistance with biological samples, and Dr Brigitte MANSHIP for English proofreading.
This work has been supported by the Exchange in Endocrinology Expertise (3E) program award of the European Union of Medical Specialists (UEMS), Section and Board of Endocrinology, Fondation ARC pour la recherche sur le cancer (PJA 20181207935 and ARCMD-DOC12020020001361), La Ligue contre le cancer, Région Rhone-Alpes-Auvergne, and La Société Française d’Endocrinologie (FRIEMM). These funding bodies had no role in the design of the study, nor in the collection, analysis, and interpretation of data, or in writing the manuscript.
Ethics approval and consent to participate
The study was performed in line with the principles of the Declaration of Helsinki and its later amendments, and was approved by the bioethics committees of Hospices Civils de Lyon (No. 21_439) and of Targu Mures Hospital (No. 5894). Patient information was recorded in a local database (PITUICARE-Lyon, registered with the French data protection agency CNIL, 16–021, and clinicaltrials.org, NCT 02854228). Informed consent was obtained from all individual patients included in the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1
: Fig. S1. Representative S100B immunohistochemical (IHC) staining and segmentation. A, Representative immunostaining patterns of S100B in tumour tissue and adjacent anterior pituitary. Left panel shows a scanned S100B immunostained section of a gonadotroph tumour. Right panels are magnified views of the observed staining in the adjacent anterior pituitary and tumour tissue present on the same slide (adjacent anterior pituitary is annotated in blue). Scale bar = 50 µm. B, Representative example of a single-cell segmentation performed with the HALO® software (Indica Labs, New Mexico, USA): S100B IHC staining (top), and the resulting segmentation (bottom) are shown. Nuclei are segmented in blue, while the positive staining is segmented in red.
Additional file 2
: Fig. S2. Association between the percentage of S100B+ cells and clinicopathological traits in Pit1 and corticotroph tumours. A, Percentage of S100B+ cells in Pit1 tumours with a Ki67 index < 3% (n = 16, median = 0.14%) versus a Ki67 index ≥ 3% (n = 4, median = 0.19%). B, Percentage of S100B+ cells in Pit1 tumours with a number of mitosis ≤ 2/10 HPFs (n = 18, median = 0.11%) versus a number of mitosis > 2/10 HPFs (n = 2, median = 0.39%). C, Percentage of S100B+ cells in non-proliferative (n = 16, median = 0.14%) versus proliferative Pit1 tumours (n = 4, median = 0.19%). D, Percentage of S100B+ cells in non-invasive (n = 9, median = 0.09%) versus invasive Pit1 tumours (n = 9, median = 0.35%). E, Percentage of S100B+ cells in corticotroph tumours with a Ki67 index < 3% (n = 5, median = 0.83%) versus a Ki67 index ≥ 3% (n = 3, median = 0.13%). F, Percentages of S100B+ cells in corticotroph tumours with a number of mitosis ≤ 2/10 HPFs (n = 7, median = 0.64%) versus a number of mitosis > 2/10 HPFs (n = 1, median = 1.11%). G, Percentage of S100B+ cells in non-proliferative (n = 7, median = 0.64%) versus proliferative corticotroph tumours (n = 1, median = 1.11%). H, Percentage of S100B+ cells in non-invasive (n = 2, median = 0.07%) versus invasive corticotroph tumours (n = 4, median = 0.74%). Graphs show median with interquartile range. Statistical test: Mann-Whitney. Abbreviations: non-significant (ns), high power fields (HPFs).
Additional file 3
: Fig. S3. Spatial association between S100B+, ERα+, Ki67+, and FSH+ cells. A, Graph showing the percentage of S100B+ cells and ERα+ cells in 20 matching tumour areas (Z1 to Z4) of 5 gonadotroph tumours (Gon#3, #5, #11, #12 and #18). The areas represented on the right side (red background) show the highest percentage of ERα+ cells. Black and red lines delineate the higher versus lower percentages of S100B+ cells, and ER + cells, respectively. The higher versus lower percentages of S100B+ cells are separated by the median. B, Graph showing the percentage of S100B+ cells and Ki67+ cells in 20 matching tumour areas (Z1 to Z4) of 5 gonadotroph tumours (Gon#3, #5, #11, #12 and #18). The areas represented on the right side (green background) show the highest percentages of Ki67+ cells. Black and green lines delineate the higher versus lower percentages of S100B+ cells, and Ki67+ cells, respectively. The higher versus lower percentages of S100B+ cells are separated by the median. C, Graph showing the percentage of S100B+ cells, FSH+ cells, and ERα+ cells, in 20 matching tumour areas (Z1 to Z4) of 5 gonadotroph tumours (Gon#3, #5, #11, #12 and #18). The 10 areas represented on the right side (grey background) show the highest percentage of S100B+ cells (the higher versus lower percentages of S100B+ cells are split by the median). Blue and red lines delineate the higher versus lower percentages of FSH+, and ERα+ cells, respectively. Abbreviations: follicle-stimulating hormone (FSH), oestrogen receptor alpha (ERα)
Additional file 4
: Fig. S4. The percentage of LH+ and SSTR2+ cells was not associated with the percentage of S100B+ cells in gonadotroph tumours. A, Percentage of LH+ cells in S100BHigh (n = 9, median percentage of LH+ cells = 12.32%) versus S100BLow areas (n = 10, median percentage of LH+ cells = 2.01%). B, Percentage of SSTR2+ cells in S100BHigh (n = 10, median percentage of SSTR2+ cells = 3.13%) versus S100BLow areas (n = 10, median percentage of SSTR2+ cells = 0.08%). S100BHigh and S100BLow expressing areas were separated by the median percentage of S100B+ cells for the 20 areas. Graphs show median with range. Statistical test: Mann-Whitney, ns=non-significant. Abbreviations: luteinizing hormone (LH), somatostatin receptor type 2 (SSTR2).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ilie, M.D., Vasiljevic, A., Chanal, M. et al. Intratumoural spatial distribution of S100B + folliculostellate cells is associated with proliferation and expression of FSH and ERα in gonadotroph tumours. acta neuropathol commun 10, 18 (2022). https://doi.org/10.1186/s40478-022-01321-y
- Tumour heterogeneity
- Folliculostellate cells
- S100B + cells
- Tumour microenvironment
- Gonadotroph adenoma