Contribution of molecular analysis to the typification of the non-functioning pituitary adenomas

Abstract

Aim

The WHO Classification of Tumours of Endocrine Organs considers the inmunohistochemical characterization of pituitary adenomas (PA) as mandatory for patient diagnosis. Recent advances in the knowledge of the molecular patterns of these tumours could complement this classification with gene expression profiling.

Methods

Within the context of the Spanish Molecular Registry of Pituitary Adenomas (REMAH), a multicentre clinical-basic research project, we analysed the molecular phenotype of 142 PAs with complete IHC and clinical information. Gene expression levels of all pituitary hormones, type 1 corticotrophin-releasing hormone receptor, dopamine receptors and arginine vasopressin receptor 1b were measured by quantitative real-time polymerase chain reaction. In addition, we used three housekeeping genes for normalization and a pool of nine healthy pituitary glands from autopsies as calibration reference standard.

Results

Based on the clinically functioning PA (FPA: somatotroph, corticotroph, thyrotroph and lactotroph adenomas), we established the interquartile range of relative expression for all genes studied in each PA subtype. That allowed molecularly the different PA subtypes, including the clinically non-functioning PA (NFPA). Afterwards, we estimated the concordance of the molecular and immunohistochemical classification with clinical diagnosis in FPA and between them in NFPA. The kappa values were higher in molecular than in immunohistochemical classification in FPA and showed a bad concordance in all NFPA subtypes.

Conclusions

According to these results, the molecular characterization of the PA complements the IHC analysis, allowing a better typification of the NFPA.

About the Authors:

Laura Sanchez-Tejada

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Writing – original draft, Writing – review & editing

Affiliation: Research Unit, Institute for Health and Biomedical Research (Isabial-Fisabio Foundation), Hospital General Universitario de Alicante, Alicante, Spain

ORCID http://orcid.org/0000-0001-5563-8923

Ruth Sanchez-Ortiga

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Endocrinology Department, Hospital General Universitario de Alicante, Alicante, Spain

Cristina Lamas

Contributed equally to this work with: Cristina Lamas, Rosa Camara, Pedro Riesgo, Carmen Fajardo

Roles Supervision, Validation, Visualization

Affiliation: Endocrinology Department, Hospital General Universitario de Albacete, Albacete, Spain

Rosa Camara

Contributed equally to this work with: Cristina Lamas, Rosa Camara, Pedro Riesgo, Carmen Fajardo

Roles Supervision, Validation, Visualization

Affiliation: Endocrinology Department, Hospital Universitario y Politecnico La Fe, Valencia, Spain

Pedro Riesgo

Contributed equally to this work with: Cristina Lamas, Rosa Camara, Pedro Riesgo, Carmen Fajardo

Roles Supervision, Visualization

Affiliation: Neurosurgery Department, Hospital de La Ribera, Alzira, Valencia, Spain

Carmen Fajardo

Contributed equally to this work with: Cristina Lamas, Rosa Camara, Pedro Riesgo, Carmen Fajardo

Roles Supervision, Validation, Visualization

Affiliation: Endocrinology Department, Hospital de La Ribera, Alzira, Spain

Francisco Ignacio Aranda

Roles Methodology, Validation, Visualization

Affiliation: Pathology Department, Hospital General Universitario de Alicante, Alicante, Spain

Antonio Pico

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Endocrinology Department, Hospital General Universitario de Alicante, Alicante, Spain

ORCID http://orcid.org/0000-0002-7549-7563

Introduction

Pituitary adenomas (PA) constitute 10%–15% of intracranial neoplasms. They are currently classified according to their clinical, biochemical, radiological and inmunohistochemical (IHC) characteristics as well as their tumour behaviour [1, 2, 3].

Clinically, PAs are divided into two main types: functioning pituitary adenomas (FPAs) and non-functioning pituitary adenomas (NFPAs). FPAs are accompanied by a clinical syndrome related to hormone hypersecretion: Cushing’s disease, acromegaly, galactorrhea-hypogonadism or hyperthyroidism. NFPAs include a cluster of pituitary tumours without endocrine manifestations of hormone overproduction, which are diagnosed incidentally or due to neuro-ophthalmological symptoms. They comprise approximately 30%–35% of pituitary tumours [4, 5].

Radiologically, two classification systems are currently used. The Hardy classification divides PAs into four grades based on their size and the invasiveness in the sella turcica [6]. The Knops classification takes into consideration the tumour invasion of the cavernous sinus according to coronal sections of magnetic resonance imaging (MRI) scans, with the readily detectable internal carotid artery serving as the radiological landmark [7].

Pathologically, different proposals have aimed to classify pituitary tumours over the years, ranging from the cellular staining properties to the latest World Health Organization (WHO) classification based on IHC and ultrastructural criteria [8]. This classification recognises five PA subtypes, depending on the aberrant hormone immunostaining: prolactin tumour (PRL), growth hormone tumour (GH), adrenocorticotropic hormone tumour (ACTH), thyroid-stimulating hormone tumour (TSH, α-subunit (SU)) and gonadotroph hormone tumour (follicle-stimulating hormone (FSH), luteinising hormone (LH) or SU) [8]. Among these, those non-gonadotroph and non-functioning tumours are known as silent (S) adenomas. Some of them are multihormonal. Finally, pituitary adenomas without IHC-detected aberrant hormone production are called null cell (NC) adenomas [3, 9].

Currently, translational research is based mainly on genetic and epigenetic characteristics. As experience accumulates in this field, molecular studies will gradually become integrated into the standardised clinical management of PA, similarly to other tumours [10,11].

Silent adenomas have been attributed to the release of pituitary hormone isoforms without biological activity or to the reduced or absent secretion of the normal pituitary hormones [3]. Null cell adenomas show ultrastructural features of a particular pituitary cell type or express cell specific transcription factors but do not secrete proteins. In these scenarios, the IHC technique is unable to define the cellular origin of the tumour. Conversely, the mRNA analysis of pituitary hormones or their specific receptors could identify the specific cell clone from which the tumour arises, allowing a molecular typification of these tumours. The clinical implications of this knowledge have yet to be clarified.

The aim of this study was to characterize molecularly the silent and immunohistochemically null cell adenomas in a large series of pituitary adenomas and to study the concordance between the IHC and molecular typification of PA.

Matherial and methods

The study was performed within the context of the Spanish Molecular Registry of Pituitary Adenomas (REMAH), a Spanish Multicentre Clinical-Basic Project [12]. The study complies with the Declaration of Helsinki and other applicable laws and received approval from the Local Ethics Committee (CEIC Hospital General Universitario de Alicante). None of the donors was from a vulnerable population and all donors or next of kin provided written informed consent that was freely given. One of the patients was 14 years old boy and his parents signed his informed consent.

Participants

We collected 142 PA specimens with complete clinical and IHC data and sufficient quantity and quality of RNA for molecular analysis. Additionally, we obtained nine normal pituitary glands from autopsies and used them as calibration reference samples.

The tumours came from four university referral hospitals for pituitary surgery. Clinical, pathological, and radiological data were collected anonymously for each sample from the REMAH database. The information retrieved included sex; age; endocrine syndrome (Cushing's syndrome, acromegaly, hyperthyroidism or amenorrhea-hypogonadism); hormonal data; previous medical, surgical or radiotherapy treatment; tumour maximum diameter (TMD); radiological MRI extension of the tumour, and IHC results. Tumours were classified as invasive (Hardy’s grade IV) or non-invasive (Hardy’s grades I–III). We present the clinical baseline characteristics in Table 1.

[Figure omitted. See PDF.]

Table 1. Clinical Baseline characteristics of patients.

https://doi.org/10.1371/journal.pone.0180039.t001

Clinical identification of pituitary adenomas

All cases presenting clinical and biochemical characteristics of a recognised endocrine syndrome were classified as FPA. We considered participants to have NFPA when they did not present any clinical or biochemical characteristics of a recognised endocrine syndrome and were diagnosed according to neuro-ophthalmological or radiological signs.

IHC identification of NFPA

According to the IHC data, we defined the following subtypes of NFPA: 1) IHC gonadotroph adenomas (IHC-GT): tumours with positive immunostaining for LH, FSH or SU; 2) IHC silent corticotroph adenomas (IHC-sCT): tumours with positive immunostaining for ACTH without Cushing’s syndrome; 3) IHC silent somatotroph adenomas (IHC-sST): tumours with positive immunostaining for GH without acromegaly; 4) IHC silent thyrotroph adenomas (IHC-sTT): tumours with positive immunostaining for TSH without hyperthyroidism; 5) IHC silent lactotroph adenomas (IHC-sLT): tumours with positive immunostaining for PRL without clinical amenorrhea-hypogonadism; 6) IHC silent multihormonal adenomas (IHC-sMHA): tumours with positive immunostaining for more than one hormone without predominance of one of them, not associated to a clinical syndrome; 7) IHC null cell adenomas (IHC-NC): tumours without positive immunostaining for any pituitary hormone.

Molecular identification of pituitary adenomas

First, we defined the interquartile ranges (IQR) of gene expression for each clinically diagnosed FPA subtype (Tables 2 and 3). Afterward, we assessed all cases with clinical diagnosis of NFPA (74/142) individually and reclassified them according to their dominant gene expression, taking as a reference the 25th percentile of their corresponding FPA, as follows: 1) molecular gonadotroph adenomas (M-GT): tumours with a dominant expression of LH, FSH or S; 2) molecular silent corticotroph adenomas (M-sCT): tumours with a dominant gene expression of POMC, AVPR1b and CRH-R1 without Cushing syndrome [13,14]; 3) molecular silent somatotroph adenomas (M-sST): tumours with a dominant gene expression of GH without acromegaly syndrome; 4) molecular silent thyrotroph adenomas (M-sTT): tumours with dominant expression of TSH without hyperthyroidism; 5) molecular silent lactotroph adenomas (M-sLT): dominant expression of PRL without amenorrhea-hypogonadism syndrome; 6) molecular silent multihormonal adenomas (M-sMHA): tumours with a dominant expression for more than one hormone without predominance of any of them, not associated to a clinical syndrome; 7) molecular null cell adenomas (M-NC): tumours that do not express any pituitary hormone gene or their receptors.

[Figure omitted. See PDF.]

Table 2. Molecular classification of Pituitary Adenomas (PA).

https://doi.org/10.1371/journal.pone.0180039.t002

[Figure omitted. See PDF.]

Table 3. Ranges of relative expression (p25-p75) of the pituitary genes.

https://doi.org/10.1371/journal.pone.0180039.t003

Although we did not find any molecular sST, sLT or sMHA in our series, we considered these categories because they have been described by others authors [15]. The expression of each hormone and receptor involved in the molecular classification of our series is shown in Fig 1.

[Figure omitted. See PDF.]

Fig 1. Relative gene expression of hormone and receptor implicated in the molecular classification of pituitary adenomas (n = 142).

Figure shows the bar chart of the gene expression levels analyzed for each subtype. Silent lactotroph and silent somatotroph are not represented because there was none in our series.

https://doi.org/10.1371/journal.pone.0180039.g001

Immunohistochemical studies

The IHC studies were performed in the respective pathology departments of the four participating hospitals. The application of immunohistochemical techniques was performed by equal in the four hospital on complete sections of formalin-fixed tissue and embedded in paraffin in automated form with Autostainer Link48 (Dako-Agilent) equipment, with the Envision high sensitivity display system (Dako). The tissue sections were subjected to dewaxing by incubation in xylol and rehydration using decreasing solutions of ethanol and water. Antigenic recovery was performed with the PT-Link apparatus, treating the sections with target retrieval solution Envision Flex, with high pH (9.0). Primary antibodies were incubated in the Autostainer for 20 to 30 minutes, depending on the setting of the antibodies. Antibodies used and their dilutions in each centre are detailed in “S1 Table”. In all cases immunostaining against the growth hormone, prolactin, adrenocorticotrophin, follicle-stimulating hormone, luteinizing hormone, thyroid-stimulating hormone, and alpha-subunit of pituitary hormones were performed. Antibodies were detected using Envision kit in an automated immunostainer. Depending on the hospital, the IHC results were expressed as either a percentage of positive cells or in semi-quantitative categories (−, +, ++ or +++). To normalize these results, two researchers blinded to clinical information revised the pathology reports and performed a pathological diagnosis, which was used in the study of concordance between molecular and IHC techniques. The criteria used to define the pathological diagnosis was: 1) in case of the results were expressed as a percentage of positive cells, only those with more than 5% were considered positive staining, and 2) in case of semi-quantitative categories only was considered the positive staining regardless of the intensity.

Molecular analysis

All the molecular studies of the tumours were performed in the Research Unit of the General University Hospital of Alicante.

All specimens were preserved immediately after surgery in RNAlater® solution at 4°C for 24–48 h and then frozen at −80°C until molecular analysis. Molecular methods were adjusted and set up to obtain the highest possible quantity of material from the samples, preserving the quality of the nucleic acids.

Before analysing the study samples, we performed efficiency analyses with two non-valuable samples and two commercial controls of total RNA from brain and total human RNA. The chosen method for extraction was the AllPrep DNA-RNA-Protein (Qiagen). We evaluated the quality and quantity of total RNA in the Nanodrop (Thermo Scientific) and the 2100 Bioanalyzer (Agilent Technologies). The protocol of retrotranscription, selection of primers and quantitative real-time-PCR (qRT-PCR) conditions were adjusted to the material and methods published by Luque et al [12]. Specifically, gene expression levels of all pituitary hormones, type 1 corticotrophin-releasing hormone receptor, arginine vasopressin receptor 1b, and five types of dopamine receptors (DR1, DR4, DR5, DR2 and its isoform DR2 long) were measured by qRT-PCR. The results of the dopamine receptors are detailed in “S2 Table”. For normalisation, three housekeeping genes (GAPDH, ACTB and HPRT-1) were analysed simultaneously and with the same protocol. The reactions were performed in 7500 Fast Real Time PCR System (Applied Biosystems), and results were analysed using SDS software (Applied Biosystems®), taking into account the three housekeeping genes. All samples and normal pituitary autopsy samples were analysed in duplicate and following the same protocol. Once normal pituitary autopsy samples were analysed separately and their concordance tested, we mixed them in equal parts into a pool that was analysed as another sample that served as calibrator. A relative quantification was established based on the ΔΔCt method.

Statistical analysis

Statistical analysis was performed with SPSS 15.0 software (IBM Software). Participant ages and tumour diameters were expressed as mean ± standard deviation (SD). Qualitative variables, including PA subtypes, were expressed as absolute and relative frequencies. Molecular variables showed a non-normal distribution (Kolmogorov-Smirnov test), therefore we chose percentiles as the measure of distributions and reported them as median plus interquartile range (IQR: p25-p75). Cohen's kappa coefficient was calculated to measure inter-classifications agreement (κ = 1 representing complete concordance and κ = 0 null concordance). P values < 0.05 were considered statistically significant.

Results

Prevalence of the different subtypes of pituitary adenomas according to identification technique

Clinical and biochemically, 70 patients (49.3%) in our series had a hormonal profile compatible with an endocrine syndrome and were considered as having FPAs, while 72 patients (50.7%) showed normal pituitary function or hypopituitarism and were considered as having NFPAs. The prevalence of the different subtypes according to the IHC or molecular identification is shown in Table 4.

[Figure omitted. See PDF.]

Table 4. Prevalence of the different pituitary subtypes according to their inmunohistochemical or molecular identification in the series studied (N = 142).

https://doi.org/10.1371/journal.pone.0180039.t004

Reliability of the molecular classification of PA

Once we molecularly identified the different subtypes of PA, we assessed the contribution of the molecular study to the identification of the cellular origin of PA, comparing it with participants’ biochemical and IHC profiles. We contrasted these techniques against the established clinical diagnosis on an individual basis (one per one in all participants).

In the total series (N = 142), the concordance between the three techniques was generally high in all subtypes, with the biochemical study showing the highest kappa values, followed by the molecular analysis (Table 5). The concordance between inmunohistochemical and molecular results was different depending on the different centres as shown in “S3 Table”.

[Figure omitted. See PDF.]

Table 5. Concordance between clinical diagnosis and independent diagnosis with each technique.

https://doi.org/10.1371/journal.pone.0180039.t005

When we considered only NFPAs (n = 74), the concordance between IHC and molecular profiles was generally low, with the highest level of concordance (κ = 0.519) shown for the sCT subtype and the lowest one (κ = 0.259) for the NC subtype (Table 6). It was not possible to calculate the concordance between IHC and molecular profiles on sST, sLT and sMHA subtypes because we did not identify these molecular subtypes in our series.

[Figure omitted. See PDF.]

Table 6. Concordance between immunohistochemical and molecular classifications for the 74 clinically NFPA.

https://doi.org/10.1371/journal.pone.0180039.t006

Discussion

In this study, we show the molecular typification of a large series of PA. The main strength of the article is the complete clinical, biochemical and immunohistochemical studies of all patients. Moreover, all molecular studies were performed in the same molecular laboratory. Conversely, the main limitation of the study is that the immunohistochemical studies were performed locally in the four University hospitals participants. As immunohistochemistry is a semi-quantitative technique and depends largely of the antibodies chosen, this fact biases significantly the inmunohistochemical results. Indeed, we observed an important variability in the IHC results between the four participating centresmeanwhile in the Lyon’s pathological series, using sensitive IHC techniques, the percentage of NC adenomas were reduced from 10% in 1992 to 1% in 2012 [16]. The variability in the IHC results between the four participating centres is detailed in “S3 Table”.

According to our molecular results, the most prevalent pituitary tumours in our series were GT followed by CT, ST, NC, TT, LT and MHAs. These results were very different from those obtained in the IHC study, where the most prevalent adenomas were the ST, followed by MHA, NC, CT, GT, LT and TT (Table 4). The prevalence of LT, both in the IHC and molecular analysis, is underestimated due to the origin of our study sample, as most LTs are not removed surgically. As previously mentioned, the estimated prevalence of pituitary adenoma subtypes has changed frequently, depending on the different series published but also on the technique used in their classification [2].

The classification of PA has also evolved with the advances in pathology. The most important development with clinical implications has been the application of IHC techniques to the characterisation of the tumours, linking the PAs with their protein secretion. IHC techniques have given support to the currently used WHO 2004 classification of PA [8]. However, while there is good concordance between clinical and IHC diagnosis in FPA, the information provided by IHC is more limited in NFPA. Indeed, our results demonstrate that many tumours classified as MHA and NC adenomas immunohistochemically may be reclassified by the molecular study.

NFPAs are frequently macroadenomas and sometimes invasive. The first line of treatment is always surgery, followed when necessary by radiotherapy. In contrast to FPAs, medical treatment of NFPA is mostly ineffective [3]. The prognostic and follow-up of NFPA depends on the results of the mentioned treatments. The overwhelming majority of PAs is benign and slow growing, so if there are no tumour remains after surgery or they have been irradiated, the frequency of clinical follow-up is usually low. However, some NFPA subtypes behave more aggressively than others and should be monitored more closely. This fact has prompted the development of classifications that tailor the prognosis of the PA to the subtype, size, degree of invasiveness and the result of the immunostaining of Ki-67 and p53 [17].

The IHC classification of PA has some limitations. While it represents a spectacular advance in the classification of these tumours, the technique is semi-quantitative and therefore observer-dependent [4,18]. Moreover, immunohistochemistry is an immunometric method that depends on the antibody used. As some PAs secrete different hormone isoforms, they cannot be recognised by specific antibodies. Indeed, Petrus et al. demonstrated that analytical methods using monoclonal antibodies could underestimate FSH levels when the balance of FSH isoforms varied [19]. That could explain why 33 tumours molecularly identified as GT (71.7% of the M-GT), were diagnosed immunohistochemically as NC or other subtypes in our series.

The term ‘null cell adenomas’ has been an important topic of discussion in the literature. Some authors prefer the term ‘non-immunoreactive tumours’ [20], speculating that IHC negativity could correspond to defects in hormone secretion or the production of biologically inactive or insufficient quantities of hormone. In our series, the prevalence of NC adenomas decreased from 12.7% in the IHC classification to 7% in the molecular one. This could be attributed to the higher sensitivity of the molecular technique and the detection of all pituitary hormone isoforms at mRNA level, especially in those tumours without enough protein translation to be detected with the IHC technique.

Similarly, the identification of silent pituitary adenomas was very different between the molecular and the IHC analyses. In our series, a considerable group of NFPA (17/74, 23%) were subclassified as M-sCT, but only 7 (9.5%) of these tumours were confirmed by IHC. The remaining 10 were classified as MHA (n = 4), NC (n = 1), GT (n = 1) and sST adenomas (n = 4). Previous reports in surgical series showed prevalences of sCT between 2.9–5.7% [21,22] similar to our IHC results but clearly inferior to that of the molecular study. The clinical significance of this could be important because it has been suggested that sCT are more aggressive, with higher rates of recurrence, than other PA [23].The silencing of ACTH has been attributed to the production of dysfunctional or high molecular weight ACTH isoforms [24]. Furthermore, Tateno et al. reported a repression of the activity of prohormone convertase (PC) in silent corticotrophs adenomas compared with the functioning CT [25]. PC is the enzyme responsible for the proteolytic processing of POMC in ACTH and other hormones. Therefore, its lack of activity would lead to a defect in the production of ACTH and would explain the negative results of immunostaining and hence the underestimation of the real prevalence of silent corticotropinomas [25,26,27]. Silent corticotropinomas also overexpress the AVPR1b and the CRH-R1 genes, similar to functioning CT[13,14]. The overexpression of these two genes together with the POMC gene in the sCT of our series (Fig 1 and Table 3) corroborated the correct molecular identification of this pituitary tumour subtype.

In the same way, gene expression identified five NFPAs as M-sTT, only one of which was confirmed by IHC. Three cases were identified as MHA (two were positive for TSH, GH and PRL, and one for TSH, GH and FSH). The last one is a discordant case, which we will discuss later.

Eleven sST were identified by the IHC method but not by the molecular analysis. The primer chosen for GH had specificity for all four transcripts of GH, so it is unlikely that the absence of M-sST represented a false negative of the molecular study. Rather, the IHC has a low specificity, attributed to the dilution and the time of incubation with the antibody and to the immunoperoxidase reaction. Authors have highlighted the high number of false positive obtained with this technique [28]. Indeed, one of the IHC-sST was clinically and biochemically diagnosed as Cushing’s syndrome and molecularly identified as M-CT. Nine IHC-sST were identified as M-GT and one as M-NC.

Similarly, four sLT were identified by the IHC method but not by the molecular analysis. Of these, one tumour was molecularly identified as M-CT, and indeed, the patient had Cushing’s syndrome. This case represents, as above, a false positive from the IHC analysis. The other three IHC-sLT were molecularly identified as M-GT.

The molecular identification of PAs proposed in the present article has shown to be highly robust to the cellular origin of these tumours. The molecular characterisation is quantitative and easily standardisable technique, so it has lower variability and higher specificity and sensitivity than the IHC technique. Indeed, the molecular classification of PA was always more consistent with the clinical manifestations than IHC (Table 4). Moreover, the molecular analysis allowed the re-classification of an important number of silent and IHC-NC tumours.

The clinical significance of molecular PA identification has not been yet established. As previously said, it is possible that M-sCT could behave with the same aggressiveness as IHC-sCT. Moreover, molecular silent tumours expressing PRL, TSH or GH genes could reduce their size if be treated with cabergoline or somatostatin analogues. Besides, the number of anti-angiogenic drugs is increasing daily, and we have recently reported that the angiogenic pathway largely depends on the cellular origin of the PA [29]. In this scenario, the complete identification of the cellular origin of PA could be essential in achieving optimal results.

The molecular classification of PA tumours also has some limitations. The mechanisms behind the silencing of some types of tumours and its clinical relevance are still unclear. Moreover, the pre-operative treatment of the tumour can affect gene expression. While in our series no tumour was irradiated before the molecular study, some patients had received prior medical treatment (Table 1). We observed that all LT treated previously with cabergoline showed lower than expected PRL gene expression levels (0.01 FC to 2.00 FC). In fact, in the present series there were some discrepant cases between the molecular, clinical and IHC diagnosis:

1. One patient with NFPA with positive immunostaining for ACTH and PRL was molecularly identified as M-TT. The patient suffered from primary hypothyroidism in substitution treatment with levothyroxine. Although previous studies have reported that patients with non-controlled primary hypothyroidism could develop TT [30], hormonal replacement usually reverses this condition. Therefore, we consider that this tumour is a sTT and that the IHC results were false positives.

2. Three patients with gene expression profile of M-GT and clinically diagnosed as NFPA had positive immunostaining for PRL. Biochemically, only one of these patients had mild hyperprolactinaemia (42.5 ng/mL), which was attributed to a lesion of the dopaminergic pathways. Therefore, we consider again that all three of these cases are probably false positives of the IHC technique.

3. One patient with clinical, biochemical and IHC diagnosis of LT did not express the PRL gene. However, the expression of dopamine receptors DR2 and DR4 in this tumour resembled that of the LT adenomas. Moreover, the tumour overexpressed the receptors DR1 and DR5 (50 and 21 FC respectively). This patient had been previously treated with high doses of cabergoline (3.5 mg/week) for four years. The most plausible explanation was that the long-term treatment with cabergoline downregulated the PRL gene and increased the DR1 and DR5 gene expression, as reported elsewhere [31,32].

4. One patient with clinical and biochemical diagnosis of LT expressed neither PRL protein nor the PRL gene and was erroneously identified via IHC and molecularly as having NC adenoma. Once again, the patient had been treated with cabergoline for a long time (more than four years). It is possible that different degrees of downregulation have sequential effects on the gene and the protein expression. In fact, other authors had shown that treatment with dopamine agonists can decrease PRL and its receptor mRNA levels [33]. Because of the complexity of the regulatory role of dopamine agonists on PRL cell function, many more studies are needed to elucidate the mechanisms by which they affect to PRL mRNA or protein levels.

5. One patient molecularly classified as M-ST had a clinical LT and showed positive immunostaining for all pituitary hormones. He was a 14-year-old patient with macroprolactinoma, in treatment with cabergoline for four years. Despite the normalisation of the PRL levels, the boy had a growth and puberty delay. The tumour expressed the GH gene at high levels (more than 3 FC) without gene expression of PRL. The most plausible explanation is that the sample analysed was normal pituitary, with a high growth hormone expression corresponding to adolescence and a lack of PRL gene expression due to the cabergoline treatment. This case highlights the need for a pathologist to select the tissue for analysis before performing any molecular study.

Conclusions

In this large series, molecular profiling showed higher accuracy than the immunohistochemical study in the identification of pituitary adenomas subtypes. This finding entails changes in the known prevalence of these tumours with possible clinical repercussions. Indeed, the prevalence of silent CT, a tumour PA subtype that often behaves aggressively, increases significantly with the molecular compared to IHC techniques. Moreover, gene expression allows reclassification of a large part of the IHC null cell adenomas. Therefore, the molecular profiling complements the pathological information of the IHC in the study of PA, although previous medical treatments may act as confounders by modifying the gene expression. Even though the not centralized performance of IHC analysis in this study could bias the IHC results, this fact increases the value of the molecular typification of PA as complement of the IHQ information in NFPA. However, it is worth highlighting the importance that an expert pathologist carefully identifies the used tissue for molecular analysis.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Antibodies used in immunohistochemical studies performed in pathology departments of the four participating hospitals.

https://doi.org/10.1371/journal.pone.0180039.s001

(DOC)

S2 Table. Relative expression [median (p25-p75)] of dopamine receptors (DR) DR1, DR4, DR5, DR2 and its isoform DR2 long, according to molecular subtypes of pituitary adenomas.

This is the S2 Table legend: PA: pituitary adenomas, DR: dopamine receptors, NC: null cell adenomas, GT: gonadotroph adenomas, CT: corticotroph adenomas, sCT: silent corticotroph adenomas, ST: somatotroph adenomas, LT: lactotroph adenomas, TT: thyrotroph adenomas, sTT: silent thyrotroph adenomas.

https://doi.org/10.1371/journal.pone.0180039.s002

(DOC)

S3 Table. Concordance between clinical and immunohistochemical / molecular diagnosis in each hospital participating in the study.

This is the S3 Table legend: Values show Cohen's kappa coefficient (κ = 1 represents complete concordance and κ = 0 the null concordance). All p-values were 0,000. IHC: Immunohistochemistry; H: Hospital; NFPA: Nonfunctioning Pituitary Adenomas; NC: non-calculable.

https://doi.org/10.1371/journal.pone.0180039.s003

(DOCX)

Acknowledgments

This work is the result of analysis of the REMAH project in its Alicante node, sponsored by the Spanish Society of Endocrinology and Nutrition (SEEN) with a grant from Novartis Oncology.The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Citation: Sanchez-Tejada L, Sanchez-Ortiga R, Lamas C, Camara R, Riesgo P, Fajardo C, et al. (2017) Contribution of molecular analysis to the typification of the non-functioning pituitary adenomas. PLoS ONE 12(7): e0180039. https://doi.org/10.1371/journal.pone.0180039

References

1. Greenman Y, Stern N (2009) Non-functioning pituitary adenomas. Best Pract Res ClinEndocrinol Metab 23: 625–638.

2. Ezzat S, Asa SL, Couldwell WT, Barr CE, Dodge WE, Lee Vance ML et al. (2004) The prevalence of pituitary adenomas: a systematic review. Cancer 101: 613–619. pmid:15274075

3. Chanson P, Raverot G, Castinetti F, Cortet-Rudelli C, Galland F, Salenave S. (2015) Management of clinically non-functioning pituitary adenoma. Ann Endocrinol (Paris) 76: 239–247.

4. Asa SL (2008) Practical pituitary pathology: what does the pathologist need to know? Arch Pathol Lab Med 132: 1231–1240. pmid:18684022

5. Al-Brahim NY, Asa SL (2006) My approach to pathology of the pituitary gland. J Clin Pathol 59: 1245–1253. pmid:17142570

6. Vezina JL, Hardy J, Yamashita M (1975) [Microadenomas and hypersecreting pituitary adenomas]. Arq Neuropsiquiatr 33: 119–127. pmid:167695

7. Knosp E, Steiner E, Kitz K, Matula C (1993) Pituitary adenomas with invasion of the cavernous sinus space: a magnetic resonance imaging classification compared with surgical findings. Neurosurgery 33: 610–617; discussion 617–618. pmid:8232800

8. DeLellis RAL, R V; Heltz Ph U and Eng CH. (2004) Tumors of Endocrine Organs. In: Press I, editor. Pathology and genetics of tumours of endocrine organs. Lyon: IARC. pp. 10–47.

9. Trouillas J, Roy P, Sturm N, Dantony E, Cortet-Rudelli C, Viennet G. (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 126: 123–135. pmid:23400299

10. Sepulveda AR, Hamilton SR, Allegra CJ, Grody W, Cushman-Vokoun AM,Funkhouser WK, et al. (2017) Molecular Biomarkers for the Evaluation of Colorectal Cancer: Guideline From the American Society for Clinical Pathology, College of American Pathologists, Association for Molecular Pathology, and the American Society of Clinical Oncology. Am J Clin Pathol 147: 221–26011.

11. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees ChA, et al. (2000) Molecular portraits of human breast tumours. Nature 406: 747–752. pmid:10963602

12. Luque RM, Ibanez-Costa A, Sanchez-Tejada L, Rivero-Cortes E, Robledo M, Madrazo-Atutxa A, et al. (2016) The Molecular Registry of Pituitary Adenomas (REMAH): A bet of Spanish Endocrinology for the future of individualized medicine and translational research. Endocrinol Nutr 63: 274–284. pmid:27091627

13. Tateno T, Izumiyama H, Doi M, Yoshimoto T, Shichiri M, Inoshita N, et al. (2007) Differential gene expression in ACTH -secreting and non-functioning pituitary tumors. Eur J Endocrinol 157: 717–724. pmid:18057378

14. de Keyzer Y, Rene P, Beldjord C, Lenne F, Bertagna X (1998) Overexpression of vasopressin (V3) and corticotrophin-releasing hormone receptor genes in corticotroph tumours. Clin Endocrinol (Oxf) 49: 475–482.

15. Cooper O, Melmed S (2012) Subclinical hyperfunctioning pituitary adenomas: the silent tumors. Best Pract Res Clin Endocrinol Metab 26: 447–460. pmid:22863387

16. Trouillas J (2014) In search of a prognostic classification of endocrine pituitary tumors. Endocr Pathol 25: 124–132. pmid:24839229

17. Raverot G, Jouanneau E, Trouillas J (2014) Management of endocrine disease: clinicopathological classification and molecular markers of pituitary tumours for personalized therapeutic strategies. Eur J Endocrinol 170: R121–132. pmid:24431196

18. Gomez-Hernandez K, Ezzat S, Asa SL, Mete O (2015) Clinical Implications of Accurate Subtyping of Pituitary Adenomas: Perspectives from the Treating Physician. Turk Patoloji Derg 31 Suppl 1: 4–17.

19. Petrus M, Rittie JL, Causse JE, Moulie N, Rhabbour M, Netter JC, et al. (1996) [Turner's syndrome and isoforms of LH and FSH. Value of polyclonal enzymatic techniques]. Arch Pediatr 3: 245–247. pmid:8785562

20. Mahta A, Haghpanah V, Lashkari A, Heshmat R, Larijani B, Tavangar SM et al. (2007) Non-functioning pituitary adenoma: immunohistochemical analysis of 85 cases. Folia Neuropathol 45: 72–77. pmid:17594597

21. Horvath E, Kovacs K, Killinger DW, Smyth HS, Platts ME, Singer W (1980) Silent corticotropic adenomas of the human pituitary gland: a histologic, immunocytologic, and ultrastructural study. Am J Pathol 98: 617–638. pmid:6244736

22. Baldeweg SE, Pollock JR, Powell M, Ahlquist J (2005) A spectrum of behaviour in silent corticotroph pituitary adenomas. Br J Neurosurg 19: 38–42. pmid:16147581

23. Cho HY, Cho SW, Kim SW, Shin CS, Park KS, Kim SY (2010) Silent corticotroph adenomas have unique recurrence characteristics compared with other nonfunctioning pituitary adenomas. Clin Endocrinol (Oxf) 72: 648–653.

24. Yamakita N, Murai T, Kawamura S, Teramachi H, Matsuhisa T, Hirata T, et al. (1999) High molecular weight corticotropin measured with immunoradiometric assay in a patient with asymptomatic pituitary corticotropinoma. Endocr J 46: 563–571. pmid:10580749

25. Tateno T, Izumiyama H, Doi M, Akashi T, Ohno K, Hirata Y. (2007) Defective expression of prohormone convertase 1/3 in silent corticotroph adenoma. Endocr J 54: 777–782. pmid:17917309

26. Ohta S, Nishizawa S, Oki Y, Yokoyama T, Namba H (2002) Significance of absent prohormone convertase 1/3 in inducing clinically silent corticotroph pituitary adenoma of subtype I—immunohistochemical study. Pituitary 5: 221–223. pmid:14558669

27. Raverot G, Wierinckx A, Jouanneau E, Auger C, Borson-Chazot F, Lachuer J, et al. (2010) Clinical, hormonal and molecular characterization of pituitary ACTH adenomas without (silent corticotroph adenomas) and with Cushing's disease. Eur J Endocrinol 163: 35–43. pmid:20385723

28. True LD (2008) Quality control in molecular immunohistochemistry. Histochem Cell Biol 130: 473–480. pmid:18648842

29. Sanchez-Ortiga R, Sanchez-Tejada L, Moreno-Perez O, Riesgo P, Niveiro M, Pico A. (2013) Over-expression of vascular endothelial growth factor in pituitary adenomas is associated with extrasellar growth and recurrence. Pituitary 16: 370–377. pmid:22990332

30. Sugimoto K, Mori K, Uchida K, Kobayashi D, Itoi K (2007) Quantitative analysis of thyroid-stimulating hormone messenger RNA and heterogeneous nuclear RNA in hypothyroid rats. Brain Res Bull 74: 142–146. pmid:17683800

31. Gruszka A, Culler MD, Melmed S (2012) Somatostatin analogs and chimeric somatostatin-dopamine molecules differentially regulate human growth hormone and prolactin gene expression and secretion in vitro. Mol Cell Endocrinol 362: 104–109. pmid:22705877

32. Gruszka A, Ren SG, Dong J, Culler MD, Melmed S (2007) Regulation of growth hormone and prolactin gene expression and secretion by chimeric somatostatin-dopamine molecules. Endocrinology 148: 6107–6114. pmid:17656461

33. Jin L, Qian X, Kulig E, Scheithauer BW, Calle-Rodrigue R, Abboud Ch, et al. (1997) Prolactin receptor messenger ribonucleic acid in normal and neoplastic human pituitary tissues. J Clin Endocrinol Metab 82: 963–968. pmid:9062514