Difference in miRNA Expression in Functioning and Silent Corticotroph Pituitary Adenomas Indicates the Role of miRNA in the Regulation of Corticosteroid Receptors

Abstract

Corticotroph pituitary adenomas commonly cause Cushing’s disease (CD), but some of them are clinically silent. The reason why they do not cause endocrinological symptoms remains unclear. We used data from small RNA sequencing in adenomas causing CD (n = 28) and silent ones (n = 20) to explore the role of miRNA in hormone secretion and clinical status of the tumors. By comparing miRNA profiles, we identified 19 miRNAs differentially expressed in clinically functioning and silent corticotroph adenomas. The analysis of their putative target genes indicates a role of miRNAs in regulation of the corticosteroid receptors expression. Adenomas causing CD have higher expression of hsa-miR-124-3p and hsa-miR-135-5p and lower expression of their target genes NR3C1 and NR3C2. The role of hsa-miR-124-3p in the regulation of NR3C1 was further validated in vitro using AtT-20/D16v-F2 cells. The cells transfected with miR-124-3p mimics showed lower levels of glucocorticoid receptor expression than control cells while the interaction between miR-124-3p and NR3C1 3′ UTR was confirmed using luciferase reporter assay. The results indicate a relatively small difference in miRNA expression between clinically functioning and silent corticotroph pituitary adenomas. High expression of hsa-miR-124-3p in adenomas causing CD plays a role in the regulation of glucocorticoid receptor level and probably in reducing the effect of negative feedback mediated by corticosteroids.

1. Introduction

Pituitary adenomas (also referred to as pituitary neuroendocrine tumors, PitNETs) represent about 10–20% of intracranial neoplasms in adults. They may originate from different kinds of secretory pituitary cells including corticotroph ACTH-secreting cells. Corticotroph adenomas commonly cause ACTH-dependent Cushing’s disease, but a significant proportion of these tumors are endocrinologically non-functioning and classified as subclinical/silent corticotroph adenomas (SCAs) [1].
CD-causing ACTH tumors are commonly small microadenomas with approximately 50% being smaller than 5 mm, which is challenging for MRI diagnostics [2]. In contrary, SCAs are commonly diagnosed due to neurological symptoms related to tumor mass at the stage of large macroadenomas. Frequently they show invasive growth and increased proliferation index [1]. According to current recommendations, SCAs are now referred to as “high-risk” pituitary adenomas which refers to their fast and invasive growth, high risk of recurrence and resistance to medical therapy [3,4]. They are recognized to be more aggressive than other clinically nonfunctioning pituitary tumors such as those of gonadotroph origin or null-cell adenomas [5].
The mechanism underlying the difference in secretory activity of CD-causing and subclinical tumors is unclear and only a few studies focused on this issue were published. The results indicated a role of the expression levels of particular genes/proteins involved in the regulation of POMC expression and pro-hormone conversion into ACTH as well as genes involved in pituitary differentiation [6,7,8,9,10,11,12,13]. However, it also appears that both active and silent corticotroph adenomas share a similar overall gene expression profile [14,15].
The aim of this study was to compare the profiles of microRNA (miRNA) expression in clinically functioning and silent corticotroph adenomas and to identify miRNAs that play a role in different ACTH secretory activity.

2. Results

2.1. Patients Characteristics

The study included 28 patients with CD and 20 patients suffering from SCA. All patients with CD had clear clinical signs and symptoms of hypercortisolism verified according to biochemical criteria including elevated midnight cortisol levels and 24 h urinary free cortisol (UFC). Patients with SCA had no clinical or biochemical signs of hypercortisolism and showed normal levels of midnight cortisol and 24 h UFC. Patients with CD had significantly higher morning serum cortisol levels than patients with SCAs (p = 0.0002) while no significant difference was observed in the morning serum ACTH levels. No difference in cortisol/ACTH ratio was observed between CD and SCA patients.
All the adenoma samples were ACTH-positive upon immunohistochemical staining against pituitary hormones (ACTH, GH, TSH, FSH, LH, α-subunit) and had characteristic ultrastructural features of corticotroph adenoma. Forty-one adenomas were positive only for ACTH, while seven ACTH-positive adenomas showed additional moderate/weak immunoreactivity for α-subunit. Increased proliferation assessed by Ki67 index ≥ 3% was observed in a similar proportion of CD and SCA patients, seven tumors causing CD and five SCAs. A higher proportion of sparsely vs. densely granulated adenomas was observed in SCAs than in CD-related adenomas, but the difference did not cross a significance threshold (p = 0.0787). No difference in the proportion of invasive/noninvasive adenomas was observed in clinically functioning and silent corticotroph adenomas.
All SCAs were macroadenomas, while tumors causing CD included 17 macroadenomas and 11 microadenomas. No significant differences in preoperative clinical parameters, including 24 h UFC, morning serum ACTH level, morning and midnight serum cortisol level, cortisol/ACTH ratio, were observed between CD patients with micro- and macroadenomas. Irrespectively, a correlation between tumors size and ACTH level (Spearman R= 0.4678; p = 0.0121) and a negative correlation between cortisol/ACTH ratio (Spearman R= −0.4015; p = 0.0342) was observed in CD patients.
No correlation was found between the remaining biochemical parameters and tumor size. Overall, the patients’ characteristics are presented in Table 1, while details including both the clinical and histopathological data are shown in Supplementary Table S1.
Table 1. Summary of clinical features of patients with Cushing’s disease and silent corticotroph adenomas.
Table

2.2. Identification of miRNAs Differentially Expressed in Corticotroph Adenomas Causing CD and Subclinical Cortiotroph Adenomas

NGS data on miRNA expression of 48 corticotroph adenomas from previous investigation were used to compare miRNA expression levels between adenomas causing CD (n = 24) and subclinical corticotroph adenomas (n = 20). Sequencing of small RNA libraries produced approximately 2,497,367 reads per sample, which were mapped to the human genome (hg19) and used for quantification of expression levels of known miRNAs, according to miRBase 22 release. Sequencing reads were annotated to 1917 miRNAs. Measurements of 1902 mature miRNAs expression were included in the analysis, after filtering out the miRNAs with low expression.
When miRNA profiles of adenomas causing CD and SCAs were compared, a total of 19 differentially expressed miRNAs were found that met the criteria of adjusted p-value < 0.05. This set included 16 miRNAs with higher expression in tumors causing CD: hsa-miR-129-2-3p, hsa-miR-129-5p, hsa-miR-124-3p, hsa-miR-132-5p, hsa-miR-129-1-3p, hsa-miR-135b-5p, hsa-miR-27a-3p, hsa-miR-10b-5p, hsa-miR-9-3p, hsa-miR-6506-3p, hsa-miR-6864-5p, hsa-let-7b-5p, hsa-miR-670-3p, hsa-miR-22-5p, hsa-miR-346 and hsa-miR-9-5p, Three miRNAs with lower expression in CD patients were found: hsa-miR-1909-3p, hsa-miR-4319 and hsa-miR-181b-3p. Details are presented in Table 2 and Figure 1A,B.
Ijms 23 02867 g001 550
Figure 1. MiRNA expression profiling in corticotroph adenomas. (A). Difference in miRNA expression between functioning and silent corticotroph adenomas. Volcano plot showing differentially expressed miRNAs. Significance and fold change thresholds are marked with dashed lines. (B). Heat map representing the expression of differentially expressed miRNAs and clustering the samples of adenomas causing Cushing’s disease (CD) and silent corticotroph adenomas (SCA). (C). The correlation between the expression levels of differentially expressed miRNAs and POMC expression or hormonal laboratory measurements in patients: morning plasma ACTH level, morning and midnight plasma cortisol levels and 24 h urinary free cortisol; * indicate p-value < 0.05; ** indicate p-value < 0.01; *** indicate p-value < 0.001
Table 2. The list of miRNAs differentially expressed in corticotroph pituitary adenomas causing CD and silent corticotroph adenomas.
Table

2.3. The Correlation of miRNA Expression and Patients’ Clinical Data

Since the clustering of the tumors based on the expression of differentially expressed miRNAs did not clearly separate functioning and silent adenomas, we determined whether the expression of the identified differentially expressed miRNAs is directly related to the results of patients’ laboratory tests as well as POMC expression, measured in tumor samples with qRT-PCR. For this purpose, Spearman’s correlation was applied to calculate a correlation matrix. We observed a significant positive correlation between 13 miRNAs out of 19 differentially expressed miRNAs and at least one of clinical laboratory parameters: serum ACTH, morning cortisol level, midnight cortisol level or 24 h UFC. For 11 miRNAs, with higher expression in patients with CD a positive correlation was observed, while a negative correlation was observed for 3 miRNAs that have lower expression in patients with CD. Four of the differentially expressed miRNAs, hsa-miR-9-3p, hsa-miR-9-5p, hsa-miR-27a-3p and hsa-miR-6506-3p, are correlated with POMC expression level in tumor tissue. The absolute value of correlation coefficient ranged between 0.31 and 0.55 which indicates a weak/moderate relationship. Details are presented in Figure 1C.

2.4. Funtional Enrichment Analysis of Differentially Expressed miRNAs

To investigate the possible functional role of the identified miRNAs with different expression levels in CD tumors and SCAs, we used the information on experimentally validated miRNA targets gathered in the miRtarbase release 8.0 database. High confidence known miRNA targets that were validated with luciferase reporter assay, reported in miRtarbase, were included in the analysis. The enrichment of the genes reported as miRNA targets of our 19 miRNAs of interest was determined with gene set over-representation analysis (GSOA) based on Gene Ontology (GO) Molecular Function and GO Biological Processes. The list of all the genes reported in miRTarbase as validated with reporter gene assay was used as reference. As a result, we found 30 GO Molecular Function terms and 293 GO Biological Processes terms as significantly enriched with genes that are targets of the 19 differentially expressed miRNAs. Top 10 enriched terms were related mainly to steroid hormone activity, regulation of transcription and regulation of stem cell differentiation, as shown in Figure 2. Details are presented in Supplementary Table S2. We paid special attention to the terms that refer to steroid hormone action, i.e., steroid hormone receptor activity (GO:0003707), nuclear receptor activity (GO:0004879), ligand-activated transcription factor activity (GO:0098531), as well as steroid hormone-mediated signaling pathway (GO:0043401) and hormone-mediated signaling pathway (GO:0009755). Importantly, the miRNA target genes that were overrepresented in these terms included NR3C1 and NR3C2 that encode for adrenal hormones glucocorticoid receptor (GR) and mineralocorticoid receptor (MR), respectively. According to the miRtarbase 9.0 database, hsa-miR-124-3p is a negative regulator of NR3C1 gene [16] while both hsa-miR-124-3p and hsa-miR-135b-5p downregulate MR [17].
Ijms 23 02867 g002 550
Figure 2. Gene set over-representation analysis of putative target genes of miRNAs differentially expressed in clinically functioning and silent corticotroph adenomas.
Using the PubMed search, we found additional evidence strongly supporting the role of hsa-miR-124-3p in the regulation of NR3C1 [18,19,20,21] as well as the role of hsa-miR-135b-5p in downregulating NR3C2 [22,23].

2.5. Comparison of the Expression of NR3C1 and NR3C2 in Corticotroph Adenomas Causing CD and Silent Adenomas

We determined the expression levels of NR3C1 and NR3C2 in corticotroph adenomas with qRT-PCR. We observed a significantly lower expression of both genes in samples from CD patients (n = 24) as compared to SCAs (n = 24); fold change (FC) 0.49 p = 0.0166 and FC 0.37 p = 0.0132, for NR3C1 and NR3C2, respectively. However, the observed difference is rather slight and a notable dispersion of the results was observed (Figure 3). The differences in NR3C1 and NR3C2 expression correspond to the differences in hsa-miR-124-3p and hsa-miR-135b-5p levels. Patients with CD have higher levels of both miRNAs and lower levels of NR3C1 and NR3C2 mRNA (Figure 3). Unfortunately, we did not find a direct correlation between the expression levels of hsa-miR-124-3p and NR3C1 or hsa-miR-135b-5p and NR3C2.
Ijms 23 02867 g003 550
Figure 3. The expression levels of NR3C1 and NR3C2 measured with qRT-PCR as well as hsa-miR-124-3p and hsa-miR-135b-5p measured with small RNA sequencing in tumor samples from CD patients and silent corticotroph adenomas; * indicate p-value < 0.05

2.6. Investigtion of miRNA-Related Regulation of NR3C1 In Vitro

Transfecting the cultured cells with miRNA mimics is the commonly used approach of in vitro validation of specific miRNA–mRNA interaction. We used mice corticotroph tumor AtT-20/D16v-F2 cells for in vitro experiment and initially verified whether these cells do express Nr3c1 and Nr3c2 genes using deposited RNAseq data from a previous experiment on AtT-20 cells (GSE132324; Gene Expression Omnibus) and qRT-PCR. This showed that the AtT-20/D16v-F2 have relatively high expression of Nr3c1 but do not express Nr3c2. Thus, we focused on the regulatory role of miR-124-3p on Nr3c1 expression. We used miRBase [24] and Targetscan [25] to determine whether miR-124-3p is evolutionarily conserved in humans and mice and whether it targets NR3C1 in both species. It confirmed that miR-124-3p is broadly conserved and it shares the same sequence of mature miRNA in humans and mice. Importantly, GR is among highly rated miR-124-3p predicted targets in both humans and mice and two highly conserved miR-124-3p binding motifs in 3′UTR of this gene were identified in these two species (Figure 4A).
Ijms 23 02867 g004 550
Figure 4. Role of mir-124-3p in regulation of glucocorticoid receptor gene. (A). Putative hsa-mir-124-3p target sites in 3′UTR of NR3C1. (B). Reduced expression of Nr3c1 gene expression and glucocorticoid receptor (GR) protein level in AtT-20/D16v-F2 cells treated with hsa-miR-124-3p mimics. (C). Results of luciferase reporter gene assay, showing the interaction between Nr3c1 3′UTR site 2 and mir-124-3p; * indicate p-value < 0.05; ns—not significant.
When we transfected AtT-20/D16v-F2 cells with miR-124-3p miRNA mimic and unspecific negative control miRNA mimic, we observed a significant decrease in Nr3c1 expression in cells treated with miR-124-3p miRNA mimic (Figure 4B). It was significantly lower than in cells treated with unspecific miRNA mimic. This difference was also clearly visible at the protein level. The GR level was reduced in cells treated with miR-124-3p miRNA mimic as compared to control (Figure 4B).
Two fragments of Nr3c1 3′UTR including each of putative miR-124-3p binding motifs were cloned in plasmid vector into 3′ region of the firefly luciferase gene. AtT-20/D16v-F2 cells were transfected with empty vector, vector with miR-124-3p binding site 1 and vector miR-124-3p binding site 2. Each of the three variants of the cells were cotransfected with miR-124-3p miRNA mimic or unspecific miRNA mimic that served as a negative control. Luminescence was developed 48 h after transfection and detected with microplate reader. As a result, we observed a significant decrease in luminescence in the cells with introduced plasmid with miR-124-3p binding site 2 treated with miR-124-3p mimic as compared to the cells transfected with the same plasmid construct but with control miRNA mimic. This observation confirms the interaction between miR-124-3p and 3′ UTR of Nr3c1 at putative binding site 2 (Figure 4C). The experiment did not confirm an interaction between miR-124-3p and 3′ UTR of Nr3c1 at binding site 1 since no significant difference of luminescence was found in cells transfected with plasmid vector harboring this binding motif treated with miR-124-3p mimic and the same cells treated with negative miRNA mimic (Figure 4C).

3. Discussion

Based on the clinical manifestation and biochemical tests results, pituitary corticotroph adenomas can be divided into functioning adenomas causing Cushing’s disease and SCAs. These two subtypes of tumors also differ in terms of some characteristics in MRI [2,26] and pathological features [27]. In contrast to CD-causing adenomas which are commonly small microadenomas, SCAs are diagnosed as macroadenomas due to neurological symptoms related to tumor mass. They are characterized by invasive growth, high risk of recurrence and resistance to medical therapy and are therefore referred to as “high-risk” pituitary adenomas according to current classification [3,4]. In our study, the SCAs were larger than functioning counterparts, as expected. A clear prevalence of women is observed among CD patients according to literature data [28], while it is not observed in patients suffering from SCAs. Our SCA group contained near equal representation of women and men as in previous reports [29,30]; however, some studies indicated female prevalence in SCAs [31].
Comparing functioning and silent corticotroph adenomas, we did not observe difference in patients’ age as well as differences in invasive growth status, ratio of adenomas with increased proliferation index and proportions of sparsely and densely granulated adenomas that may suggest the lack of difference in the tumors’ “aggressiveness”. Importantly, limitations for generalization of our results should be noted. The number of patients included in the analysis is relatively low and the group is not representative of the general population, especially in the case of patients suffering from Cushing’s disease. Since the main goal of our study was a molecular profiling of tumor tissue, we intentionally preselected large adenomas, which allowed us to have enough tissue for DNA/RNA isolation and successful molecular procedures.
In our investigation, we observed a negative correlation between cortisol/ACTH ratio and tumor volume in functioning corticotroph adenomas as described previously [32]. However, we did not observe any difference between micro- and macroadenomas causing CD as compared to SCAs (data not shown) as was found in previous studies [12].
The reason why some of corticotroph adenomas exhibit excessive hormone secretion and the others remain clinically silent is unclear and only few attempts have been made to determine the possible molecular mechanism underlying this difference in secretory activity. They were mainly focused on investigating the expression of the selected genes or proteins by comparing subclinical and functioning corticotroph adenomas. These studies indicated different expression levels of prohormone convertase 1/3 POMC, genes encoding somatostatin receptors, corticotropin releasing hormone receptor 1, vasopressin receptor (V1BR), corticosteroid 11-beta-dehydrogenase as well as NEUROD1 and TPIT [6,7,8,9,10,11,12,13]. However, whole transcriptome studies indicated that adenomas causing CD and subclinical corticotroph adenomas share a very common gene expression profile and a very low number of differentially expressed genes can be found by comparing transcriptome of silent and CD-causing ACTH tumors [14,15].
In our study, we determined the miRNA expression profile of 28 clinically functioning adenomas and 20 SCAs with next-generation sequencing of small RNA fraction. This allowed for the quantification of over 1900 miRNA annotated to current version of miRbase database and comparing their expression in two groups of tumor samples. We found a significant difference only in the expression levels of 19 miRNAs, that represent less than 1% of the miRNAs included in the analysis. This result resembles the observation from previous comparison of whole transcriptome profiles in functioning adenomas and SCAs where only 34 differentially expressed genes were found. Generally, both observations indicate a very common molecular profile of corticotroph adenomas, regardless of the functional status.
In our study, the expression levels of 13 out of 19 identified differentially expressed miRNAs were also correlated with peripheral ACTH/cortisol levels, further supporting the role of these miRNAs in secretory activity of corticotroph adenomas.
The possible role of miRNA in subclinical nature of SCAs was addressed in only one previous study by García-Martínez A et al. [33]. The authors compared the expression of 5 miRNAs in 24 functioning and 23 silent adenomas and observed a difference in hsa-miR-200a and hsa-miR-103 levels [33]. Their results were not confirmed by our investigation since these two miRNAs were not found among differentially expressed miRNAs. In our data, very a similar expression level of hsa-miR-200a was observed in clinically functioning and silent adenomas. In turn, a slightly higher expression of hsa-miR-103a-3p was observed in SCAs as previously reported, but the difference did not cross the significance threshold level. We should note that different methods were used for these two studies and technical and analytical differences could result in this discrepancy.
Since miRNAs play a role in gene regulation, their effect should be investigated in the context of the function of targeted genes. The interaction between miRNA and its target mRNA 3′UTR can be predicted with in silico tools. Unfortunately, prediction results can be very difficult to interpret since a huge number of predicted interactions can be found for some miRNAs. For example, when using the Targetescan (http://www.targetscan.org; accessed on 28 February 2022) prediction tool [25], over 4000 target genes were predicted for each hsa-miR-9-3p, hsa-miR-1909-3p, hsa-miR-22-5p and hsa-miR-181b-3p that we found as differentially expressed in CD and SCA. Therefore, to investigate a possible functional relevance of differentially expressed miRNAs we used a database of experimentally validated miRNA targets [34]. Gene set over-representation analysis of miRNA target genes indicated their enrichment in the pathways of steroid hormone nuclear receptors functioning. This result indicates that miRNAs that have different expression levels in CD and SCAs play a role in the regulation of expression of genes involved in steroid hormone signaling at hormone receptor level. It is especially interesting since this group of compounds includes adrenal hormones that play a role in the regulation of the hypothalamic–pituitary–adrenal (HPA) axis.
The particular enriched miRNA target genes included NR3C1 and NR3C2 that encode for corticosteroid hormone receptors (GR and MR, respectively). Both receptors are located in the cytoplasm where they bind glucocorticoids. Upon ligand binding, they are translocated to nucleus where they form dimers on DNA at glucocorticoid response elements (GREs). Glucocorticoid and mineralocorticoid receptors directly regulate the expression of target genes and/or influence the expression indirectly through the interaction with other transcription factors [35].
Glucocorticoids play a role in the basic mechanism of negative feedback of HPA axis. They act on hypothalamus, where high cortisol levels reduce secretion of corticotropin-releasing hormone (CRH), thus they directly reduce stimulation of ACTH secretion by anterior pituitary lobe. Glucocorticoids also inhibit the activity of pituitary cells indirectly. Corticotroph cells express GRs and their activation results in the reduction of POMC expression and secretion of ACTH [36,37]. In pituitary corticotroph adenomas, NR3C1 point mutations and loss of heterozygosity in NR3C1 locus were identified [38]. These mutations seem to affect the secretory activity and result in tumor resistance to corticosteroids [39]. Reduced expression of corticosteroid receptors in corticotroph adenomas has been reported in patients with resistance to high doses of dexamethasone [40]. These data indicate a role of GR in secretory activity of clinically functioning corticotroph adenomas. The expression of corticosteroid genes was previously investigated in CD-causing tumors and SCAs and no significant differences were found. However, it is worth noting that a low number of SCA patients was included in these studies: n = 9 [13], n = 8 [11] and n = 2 [41].
According to previously published results, hsa-miR-124-3p is a negative regulator of NR3C1 [16,18,19,20,21]. This was observed in acute lymphoblastic leukemia [19], adipocytes [20] and human embryonic kidney cells [21], where the reduced expression of NR3C1 upon an increase in hsa-miR-124-3p as well as a direct interaction between this miRNA and 3′UTR of GR gene were observed. Some additional clinical observations also suggest the role of hsa-miR-124-3p in the regulation of the response to cortiosteroids in patients with acute-on-chronic liver failure [18] and lymphoblastic leukemia [19]. Hsa-miRNA-124 also mediates corticosteroid resistance in T-cells of sepsis patients through the downregulation of GR [42].
Our analysis of the expression level of NR3C1 in corticotroph adenomas showed that tumors causing CD have lower gene expression and accordingly they exhibit higher levels of hsa-miR-124-3p. Subsequently, the role of hsa-miR-124-3p in NR3C1 downregulation was confirmed in mice AtT-20/D16v-F2 corticotroph cells using miRNA mimics and reporter gene assay. Transfection of AtT-20/D16v-F2 cells with hsa-miR-124-3p mimics resulted in reduced NR3C1 mRNA expression and GR protein level. We also confirmed the interaction between hsa-miR-124-3p and one of two predicted binding motifs in 3′UTR of NR3C1 with luciferase reporter gene assay. Since sequences of hsa-miR-124-3p and target sequence in 3′UTR of NR3C1 mRNA are the same in mice and in humans, we believe that results showing the regulation of the GR-encoding gene in mice AtT-20/D16v-F2 cells are also relevant to humans. Together, the available data indicate that in pituitary corticotrophs, hsa-miR-124-3p downregulates the expression of the GR gene. Since this receptor mediates the response of pituitary cells to cortisol, the expression of hsa-miR-124-3p appears to be an important element in the regulation of secretory activity of corticotroph cells. Based on these results, we can hypothesize that in CD, a high level of hsa-miR-124-3p contributes to lowering of GR expression and in consequence it plays a role in lowering the effect of glucocorticoid feedback on the activity of corticotroph adenoma. Hsa-miR-124-3p and hsa-miR-135b-5p can downregulate the expression level of MR, as proven in model HeLa cells [17]. Expression of both miRNAs is higher in corticotroph adenomas causing CD which corresponds to the lower expression of the NR3C2 gene in these tumors as compared to SCAs. Since the role of the MR receptor expression in pituitary cells is poorly understood, the functional implication of this observation is much less clear than in the case of GR downregulation. MR and GR have similar amino acid sequences, especially in DNA-binding domain, but they differ in affinity to corticosteroids. MR is specific for both mineralocorticoids and glucocorticoids while GR is specific predominantly for glucocorticoids. MRs have much higher affinity for glucocorticoids than GRs and are activated at basal glucocorticoid conditions, while GR occupancy is increased when glucocorticoid levels rise during the circadian peak or stress. Due to these differences, these two receptors play slightly different roles, despite the fact that they share a number of target genes [43]. MR expression is considered more tissue-specific than GR and was reported to be the most prevalent in kidney and adipose tissue but also in the hippocampus and hypothalamus [44]. However, the available databases of human expression pattern such as the Genotype-Tissue Expression project (https://gtexportal.org; accessed on 10 December 2021) or Protein atlas (https://www.proteinatlas.org; accessed on 10 December 2021) indicate that MR is widely expressed in multiple human tissues and organs including the pituitary gland. Unfortunately, a role of MR receptor in pathogenesis of pituitary tumors remains unknown.
AtT-20 cells, which are the only available cell line model of corticotroph adenoma, do not express MR receptor, thus the procedure of experimental validation of the role of miRNA in NR3C2 silencing is not applicable. With a lack of experimental data on the exact role of MR, we can only hypothesize that miRNA-mediated silencing of NR3C2 may have the similar effect on HPA axis feedback as silencing of NR3C1. It may enhance ACTH secretion by reducing the direct inhibitory effect of glucocorticoids on neoplastic pituitary corticotrophs.
The difference in expression of hsa-miR-124-3p and hsa-miR-135b-5p between subclinical and CD-causing adenomas is not big, thus we suppose that high expression of these miRNAs is not the only cause of difference in ACTH secretion. Presumably this is one of the mechanisms in the regulation of corticotrophs’ secretory activity. The model of miRNA-based corticosteroid receptor regulation does not undermine the role of previously described differences in the expression of convertase 1/3, POMC, somatostatin receptors or corticotropin releasing hormone receptor 1 or genes involved in differentiation of pituitary cells [6,7,8,9,10,11,12,13]. When considering the complex nature of the regulation of ACTH secretion, it can be assumed that multiple mechanisms may be involved in the silent character of subclinical adenomas. The low number of identified differentially expressed miRNAs or genes in silent and clinically functioning adenomas probably results from the intertumoral molecular heterogeneity of SCAs. This is also in line with clinical evidence indicating that some silent corticotroph adenomas can transform into clinically functioning ones while the others remain silent [1].
The misregulation of GR expression or NR3C1 mutation may have important therapeutical implications in CD patients. Non-selective GR antagonist Mifepristone was officially approved for treatment in patients with Cushing’s syndrome [45] while another new GR inhibitor, Relacorilant (CORT125134), is under clinical investigation for its use in this group of patients [46]. The further studies will be required to assess the role of GR abnormalities in response to GR-targeting treatment in CD.
In our study, we focused mainly on the role of hsa-miR-124-3p and hsa-miR-135b-5p in the regulation of corticosteroid receptors, but the role of other differentially expressed miRNAs can also be elucidated, based on the function of putative target genes. In the pathways enrichment analysis of the putative targets, molecular functions related to transcriptional regulation were found among the top processes. Interestingly, five miRNAs, i.e., hsa-miR-132-5p, hsa-miR-135b-5p, hsa-miR-27a-3p, hsa-miR-9-3p and hsa-miR-9-5p, were previously reported to downregulate the expression of FOXO1 transcription factor [47,48,49,50,51]. FOXO1 plays an important role in the differentiation of pituitary cells [52] and secretion of gonadotropic hormones [53,54] and prolactin [55]. The role of FOXO1 in pituitary corticotroph cells was not investigated but it was shown to regulate POMC expression in POMC hypothalamic neurons [56]. In POMC, neurons of arcuate nucleus FOXO1 directly suppresses POMC expression. A similar mechanism was also observed in prolactin pituitary adenomas where FOXO1 suppresses the promoter of PRL gene [55]. It is possible that high expression of hsa-miR-132-5p, hsa-miR-135b-5p, hsa-miR-27a-3p, hsa-miR-9-3p and hsa-miR-9-5p in pituitary corticotroph adenomas reduces the level of FOXO1 and eventually contributes to the upregulation of POMC expression. In our data from corticotroph adenomas, we observed the correlation between levels of hsa-miR-9-3p/hsa-miR-9-5 and POMC expression, which also supports this concept, but the exact role of miRNAs in possible FOXO1-related regulation of secretory activity of corticotroph cells requires further functional investigation.

4. Materials and Methods

4.1. Patients and Tissue Samples

Pituitary tumor samples from 48 patients were collected during transsphenoidal surgery. Formalin-fixed and paraffin-embedded (FFPE) tissue samples, including 28 samples from patients with Cushing’s disease and 20 samples of SCA were used for the study. Diagnosis of hypercortisolism was based on standard hormonal criteria: increased UFC in three 24 h urine collections, disturbances of cortisol circadian rhythm, increased serum cortisol levels accompanied by increased or not suppressed plasma ACTH levels at 8.00 and a lack of suppression of serum cortisol levels to <1.8 µg/dL during an overnight dexamethasone suppression test (1 mg at midnight). The pituitary etiology of Cushing’s disease was confirmed based on the serum cortisol levels or UFC suppression < 50% with a high-dose dexamethasone suppression test (2 mg q.i.d. for 48 h) or a positive result of a corticotrophin-releasing hormone stimulation test (100 mg i.v.) and positive pituitary magnetic resonance imaging.
ACTH levels were assessed using IRMA (ELSA-ACTH, CIS Bio International, Gif-sur-Yvette Cedex, France). The analytical sensitivity was 2 pg/mL (reference range: 10–60 pg/mL). Serum cortisol concentrations were determined by the Elecsys 2010 electrochemiluminescence immunoassay (Roche Diagnostics, Mannheim, Germany). Sensitivity of the assay was 0.02 μg/dL (reference range: 6.2–19.4 μg/dL). UFC was determined after extraction (liquid/liquid with dichloromethane) by electrochemiluminescence immunoassay (Elecsys 2010, Roche Diagnostics)—reference range: 4.3–176 μg/24 h.
All the tumors underwent detailed histopathological diagnosis including immunohistochemical staining with antibodies against particular pituitary hormones (ACTH, GH, TSH, FSH, LH, α-subunit) and Ki67 as well as ultrastructural analysis with electron microscopy.
The SCAs were characterized by the following clinicopathological criteria: positive immunohistochemical staining for ACTH, lack of signs and symptoms of hypercortisolism (Cushing’s syndrome), negative hormonal evaluation and non-compliance with diagnostic criteria of the CD.
Macroadenoma was defined as an adenoma with at least one diameter exceeding 10 mm, and the tumor volume was assessed with the diChiro Nelson formula (height × length × width × π/6). Invasive growth of the tumors was evaluated using Knosp grading [57]. Adenomas with Knosp grades 0, 1 and 2 were considered non-invasive, while those with Knosp 3 and 4 were considered invasive.
Forty-three patients had a clear history of not using any drugs that control the overproduction of the cortisol or ACTH (ketoconazole, mitotane, metyrapone, osilodrostat, mifepristone, pasireotide) before surgical treatment. The information on preoperative pharmacological treatment was not available for 5 patients.
Tumor tissue content of each FFPE sample ranged between 80 and 100% (median 99%), as assessed with histopathological examination. Patients’ characteristics are presented in Table 1 and details on each patient’s data are available in Supplementary Table S1.
The study was approved by the local Ethics Committee of Maria Sklodowska-Curie National Research Institute of Oncology in Warsaw, Poland. Each patient provided informed consent for the use of tissue samples for scientific purposes.
Total RNA from FFPE samples was purified with RecoverAll™ Total Nucleic Acid Isolation Kit for FFPE tissue (Thermo Fisher Scientific, Waltham, MA, USA) and measured using NanoDrop 2000 (Thermo Fisher Scientific). RNA was stored at −70 °C.

4.2. Micro RNA Expression Profiling

For comparing the miRNA expression profiles in CD-causing and clinically silent adenomas, NGS data from our previous investigation of miRNA expression in corticotroph adenomas were used. The dataset is available at Gene Expression Omnibus, accession no GSE166279. Sequencing of small RNA fraction was performed in 48 tumor samples (28 CD patients and 20 SCA patients) with ion semiconductor sequencing technology, as described previously [58]. Briefly, Ion Total RNA-Seq Kit v2 (Thermo Fisher Scientific) was used for sequencing library construction, Ion Xpress™ RNA-Seq Barcode Kit was used for hybridization and ligation of RNA adapters. RNA reverse transcription and subsequent cDNA purification and library size selection were performed using Nucleic Acid Binding Beads and verified using Bioanalyzer 2100 with High Sensitivity DNA Kit (Agilent, Santa Clara, CA, USA). Ion Chef instrument, with Ion PI™ Hi-Q™ Chef Kit (Thermo Fisher Scientific) and Ion Proton sequencer (Thermo Fisher Scientific) were used for library preparation and sequencing, respectively.
BamToFastq package was applied for converting unmapped bam files into fastq files. miRDeep2 was applied for read mapping to known human miRNAs (according to miRBase release 22) and reads quantification. Data normalization and differential expression analysis were performed using DESeq2. Filtration for low-expression miRNAs was applied as described previously. FC of expression calculated as the ratio of the normalized read-count value in CD-causing and silent adenomas was used as a measure of expression difference. Adjusted p-value < 0.05 was used as significance threshold. MiRtarbase release 9.0 database [34] was used to identify known miRNA target genes. PANTHER (http://pantherdb.org; accessed on 10 December 2021) [59] was used for gene set over-representation analysis.

4.3. qRT-PCR gene Expression Analysis

One microgram of RNA was subjected to reverse transcription with Transcriptor First Strand cDNA Synthesis Kit (Roche Diagnostics). qRT-PCR reaction was carried out in 384-well format using 7900HT Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) and Power SYBR Green PCR Master Mix (Thermo Fisher Scientific) in a volume of 5 μL, containing 2.25 pmol of each primer. The samples were amplified in triplicates. GAPDH was used as reference gene. Delta Ct method was used to calculate the relative expression level. PCR primers’ sequences are presented in Supplementary Table S3.

4.4. Cell Line Culture and miRNA Mimic Transfection

AtT-20/D16v-F2 cells were purchased from ATCC collection and cultured in DMEM medium supplemented with 10% FBS, as recommended. MiRCURY LNA miRNA Mimics including hsa-miR-124-3p mimic (YM00471256, Qiagen, Hilden, Germany) and negative control mimic (YM00479902-ADB, Qiagen) were used. AtT-20/D16v-F2 cells were seeded at 5 × 104 per well of a 24-well plate in culture medium and transfected with 50 nM miRNA with 1% (v/v) HiPerFect Transfection Reagent (Qiagen), according to the manufacturer’s instructions. The next day, the culture medium was changed. In total, 48 h after transfection the cells were harvested and subjected to isolation of total RNA with RNeasy Mini Kit (Qiagen). The expression of the putative hsa-miR-124-3p target gene was determined with qRT-PCR.

4.5. Luciferase Reporter Gene Assay

Hsa-miR-124-3p target sites in 3′UTR of NR3C1 were determined with Targetscan [25]. Each of two predicted hsa-miR-124-3p target sites were cloned into pmirGLO Dual-Luciferase miRNA Target Expression Vector (Promega, Madison, WI, USA). AtT-20/D16v-F2 cells (2 × 104/well) were seeded onto a 96-well plate in 100 µL culture medium. The next day, the cells were transfected with 100 ng of each plasmid vector, independently using 0.25% (v/v) lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA) in 10 µL of DMEM. The cells were subsequently transfected with either hsa-miR-124-3p mimic (YM00471256, Qiagen) or negative control mimic (YM00479902-ADB, Qiagen) in a final concentration of 50 nM using HiPerfectReagent (Qiagen). Culture medium was changed on the next day. Luciferase activity was measured with One-Glo Luciferase Assay System (Promega) 48 h after transfection.

4.6. Western Blotting

Cells were lysed in ice cold RIPA buffer, incubated for 30 min in 4 °C and centrifuged at 12,500× g rpm for 20 min at 4 °C. Samples were resolved using SDS-PAGE and electrotransferred to polyvinylidene fluoride membranes (PVDF) (Thermo Fisher). GR protein was detected with monoclonal anti-Glucocorticoid Receptor antibody (ab183127, Abcam, Cambridge, UK), and secondary anti-rabbit antibody conjugated to HRP (#7074, Cell Signaling, Beverly, MA, USA). Glyceraldehyde-3-Phosphate Dehydrogenase (#MAB374, Millipore, Bedford, MA, USA) detected with mouse HRP-conjugated antibody (#7076 Cell Signaling) served as control. Visualization was performed with SuperSignal West Pico Chemiluminescent Substrate (Thermo Fisher Scientific) and CCD digital imaging system Alliance Mini HD4 (UVItec Limited, Cambridge, UK).

4.7. Statistical Analysis

A two-sided Mann–Whitney U-test was used for analysis of continuous variables. The Spearman correlation method was used for correlation analysis. Significance threshold of α = 0.05 was adopted. Data were analyzed using GraphPad Prism 6.07 (GraphPad Software, La Jolla, CA, USA). Hierarchical clustering analysis was carried out with Cluster 3.0, and the results were visualized using TreeView 1.6 software (Stanford University School of Medicine, Stanford, CA, USA).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23052867/s1.

Author Contributions

Conceptualization, M.M. and M.B.; Methodology, M.B. and B.J.M.; Software, J.B.; Formal analysis, P.K., B.J.M. and M.B.; Investigation, B.J.M., P.K., N.R., M.B. and M.P.; Resources, J.K., G.Z., A.S. and T.M.; Data curation, J.B., B.J.M. and M.B.; Writing—original draft preparation, M.B., P.K. and B.J.M.; Writing—review and editing, all the authors; Visualization, M.B. and B.J.M.; Supervision, M.M.; Project administration M.B.; Funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Science Centre, Poland, grant number 2021/05/X/NZ5/01874.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the local Ethics Committee of Maria Sklodowska-Curie Institute—Oncology Center in Warsaw, Poland; approval no. number 44/2018, date of approval 26 July 2018.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data from next-generation sequencing of small RNA fraction of 48 corticotroph adenoma samples are available at Gene Expression Omnibus, accession no GSE166279.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ben-Shlomo, A.; Cooper, O. Silent Corticotroph Adenomas. Pituitary 2018, 21, 183–193. [Google Scholar] [CrossRef] [PubMed]
  2. Vitale, G.; Tortora, F.; Baldelli, R.; Cocchiara, F.; Paragliola, R.M.; Sbardella, E.; Simeoli, C.; Caranci, F.; Pivonello, R.; Colao, A. Pituitary Magnetic Resonance Imaging in Cushing’s Disease. Endocrine 2017, 55, 691–696. [Google Scholar] [CrossRef] [PubMed]
  3. Kontogeorgos, G.; Thodou, E.; Osamura, R.Y.; Lloyd, R.V. High-Risk Pituitary Adenomas and Strategies for Predicting Response to Treatment. Hormones 2022, 1, 3. [Google Scholar] [CrossRef] [PubMed]
  4. Osamura, R.Y.; Grossman, A.; Korbonits, M.; Kovacs, K.; Lopes, M.B.S.; Matsuno, A.; Trouillas, J. WHO Classification of Tumours of Endocrine Organs, 4th ed.; Lloyd, R.V., Osamura, R.Y., Rosari, J., Eds.; IARC Press: Lyon, France, 2017. [Google Scholar]
  5. Jiang, S.; Chen, X.; Wu, Y.; Wang, R.; Bao, X. An Update on Silent Corticotroph Adenomas: Diagnosis, Mechanisms, Clinical Features, and Management. Cancers 2021, 13, 6134. [Google Scholar] [CrossRef]
  6. Tateno, T.; Kato, M.; Tani, Y.; Oyama, K.; Yamada, S.; Hirata, Y. Differential Expression of Somatostatin and Dopamine Receptor Subtype Genes in Adrenocorticotropin (ACTH)- Secreting Pituitary Tumors and Silent Corticotroph Adenomas. Endocr. J. 2009, 56, 579–584. [Google Scholar] [CrossRef]
  7. Gabalec, F.; Beranek, M.; Netuka, D.; Masopust, V.; Nahlovsky, J.; Cesak, T.; Marek, J.; Cap, J. Dopamine 2 Receptor Expression in Various Pathological Types of Clinically Non-Functioning Pituitary Adenomas. Pituitary 2012, 15, 222–226. [Google Scholar] [CrossRef]
  8. Righi, A.; Faustini-Fustini, M.; Morandi, L.; Monti, V.; Asioli, S.; Mazzatenta, D.; Bacci, A.; Foschini, M.P. The Changing Faces of Corticotroph Cell Adenomas: The Role of Prohormone Convertase 1/3. Endocrine 2017, 56, 286–297. [Google Scholar] [CrossRef]
  9. Ohta, S.; Nishizawa, S.; Oki, Y.; Yokoyama, T.; Namba, H. Significance of Absent Prohormone Convertase 1/3 in Inducing Clinically Silent Corticotroph Pituitary Adenoma of Subtype I—Immunohistochemical Study. Pituitary 2002, 5, 221–223. [Google Scholar] [CrossRef]
  10. Tateno, T.; Izumiyama, H.; Doi, M.; Akashi, T.; Ohno, K.; Hirata, Y. Defective Expression of Prohormone Convertase 1/3 in Silent Corticotroph Adenoma. Endocr. J. 2007, 54, 777–782. [Google Scholar] [CrossRef]
  11. Tateno, T.; Izumiyama, H.; Doi, M.; Yoshimoto, T.; Shichiri, M.; Inoshita, N.; Oyama, K.; Yamada, S.; Hirata, Y. Differential Gene Expression in ACTH -Secreting and Non-Functioning Pituitary Tumors. Eur. J. Endocrinol. 2007, 157, 717–724. [Google Scholar] [CrossRef]
  12. Nagaya, T.; Seo, H.; Kuwayama, A.; Sakurai, T.; Tsukamoto, N.; Nakane, T.; Sugita, K.; Matsui, N. Pro-Opiomelanocortin Gene Expression in Silent Corticotroph-Cell Adenoma and Cushing’s Disease. J. Neurosurg. 1990, 72, 262–267. [Google Scholar] [CrossRef] [PubMed]
  13. Raverot, G.; Wierinckx, A.; Jouanneau, E.; Auger, C.; Borson-Chazot, F.; Lachuer, J.; Pugeat, M.; Trouillas, J. Clinical, Hormonal and Molecular Characterization of Pituitary ACTH Adenomas without (Silent Corticotroph Adenomas) and with Cushing’s Disease. Eur. J. Endocrinol. 2010, 163, 35–43. [Google Scholar] [CrossRef] [PubMed]
  14. Neou, M.; Villa, C.; Armignacco, R.; Jouinot, A.; Raffin-Sanson, M.L.; Septier, A.; Letourneur, F.; Diry, S.; Diedisheim, M.; Izac, B.; et al. Pangenomic Classification of Pituitary Neuroendocrine Tumors. Cancer Cell 2020, 37, 123–134.e5. [Google Scholar] [CrossRef] [PubMed]
  15. Bujko, M.; Kober, P.; Boresowicz, J.; Rusetska, N.; Paziewska, A.; Dabrowska, M.; Piaścik, A.; Pȩkul, M.; Zieliński, G.; Kunicki, J.; et al. USP8 Mutations in Corticotroph Adenomas Determine a Distinct Gene Expression Profile Irrespective of Functional Tumour Status. Eur. J. Endocrinol. 2019, 181, 615–627. [Google Scholar] [CrossRef]
  16. Vreugdenhil, E.; Verissimo, C.S.L.; Mariman, R.; Kamphorst, J.T.; Barbosa, J.S.; Zweers, T.; Champagne, D.L.; Schouten, T.; Meijer, O.C.; De Ron Kloet, E.; et al. MicroRNA 18 and 124a Down-Regulate the Glucocorticoid Receptor: Implications for Glucocorticoid Responsiveness in the Brain. Endocrinology 2009, 150, 2220–2228. [Google Scholar] [CrossRef]
  17. Sõber, S.; Laan, M.; Annilo, T. MicroRNAs MiR-124 and MiR-135a Are Potential Regulators of the Mineralocorticoid Receptor Gene (NR3C2) Expression. Biochem. Biophys. Res. Commun. 2010, 391, 727–732. [Google Scholar] [CrossRef]
  18. Wang, X.; Xu, H.; Wang, Y.; Shen, C.; Ma, L.; Zhao, C. MicroRNA-124a Contributes to Glucocorticoid Resistance in Acute-on-Chronic Liver Failure by Negatively Regulating Glucocorticoid Receptor Alpha. Ann. Hepatol. 2020, 19, 214–221. [Google Scholar] [CrossRef]
  19. Liang, Y.N.; Tang, Y.L.; Ke, Z.Y.; Chen, Y.Q.; Luo, X.Q.; Zhang, H.; Huang, L. Bin MiR-124 Contributes to Glucocorticoid Resistance in Acute Lymphoblastic Leukemia by Promoting Proliferation, Inhibiting Apoptosis and Targeting the Glucocorticoid Receptor. J. Steroid Biochem. Mol. Biol. 2017, 172, 62–68. [Google Scholar] [CrossRef]
  20. Liu, K.; Zhang, X.; Wei, W.; Liu, X.; Tian, Y.; Han, H.; Zhang, L.; Wu, W.; Chen, J. Myostatin/Smad4 Signaling-Mediated Regulation of Mir-124-3p Represses Glucocorticoid Receptor Expression and Inhibits Adipocyte Differentiation. Am. J. Physiol. Endocrinol. Metab. 2019, 316, E635–E645. [Google Scholar] [CrossRef]
  21. Roy, B.; Dunbar, M.; Shelton, R.C.; Dwivedi, Y. Identification of MicroRNA-124-3p as a Putative Epigenetic Signature of Major Depressive Disorder. Neuropsychopharmacology 2017, 42, 864–875. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Che, X.; Yang, N.; Bai, Z.; Wu, Y.; Zhao, L.; Pei, H. MiR-135b-5p Promotes Migration, Invasion and EMT of Pancreatic Cancer Cells by Targeting NR3C2. Biomed. Pharmacother. 2017, 96, 1341–1348. [Google Scholar] [CrossRef] [PubMed]
  23. Huang, Y.; Wang, Y.; Ouyang, Y. Elevated MicroRNA-135b-5p Relieves Neuronal Injury and Inflammation in Post-Stroke Cognitive Impairment by Targeting NR3C2. Int. J. Neurosci. 2021, 132, 58–66. [Google Scholar] [CrossRef] [PubMed]
  24. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. MiRBase: From MicroRNA Sequences to Function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef] [PubMed]
  25. McGeary, S.E.; Lin, K.S.; Shi, C.Y.; Pham, T.M.; Bisaria, N.; Kelley, G.M.; Bartel, D.P. The Biochemical Basis of MicroRNA Targeting Efficacy. Science 2019, 366, eaav1741. [Google Scholar] [CrossRef] [PubMed]
  26. Kasuki, L.; Antunes, X.; Coelho, M.C.A.; Lamback, E.B.; Galvão, S.; Silva Camacho, A.H.; Chimelli, L.; Ventura, N.; Gadelha, M.R. Accuracy of Microcystic Aspect on T2-Weighted MRI for the Diagnosis of Silent Corticotroph Adenomas. Clin. Endocrinol. 2020, 92, 145–149. [Google Scholar] [CrossRef] [PubMed]
  27. Thodou, E.; Argyrakos, T.; Kontogeorgos, G. Galectin-3 as a Marker Distinguishing Functioning from Silent Corticotroph Adenomas. Hormones 2007, 6, 227–232. [Google Scholar]
  28. Valassi, E.; Santos, A.; Yaneva, M.; Tóth, M.; Strasburger, C.J.; Chanson, P.; Wass, J.A.H.; Chabre, O.; Pfeifer, M.; Feelders, R.A.; et al. The European Registry on Cushing’s Syndrome: 2-Year Experience. Baseline Demographic and Clinical Characteristics. Eur. J. Endocrinol. 2011, 165, 383–392. [Google Scholar] [CrossRef]
  29. Langlois, F.; Shao, D.; Lim, T.; Yedinak, C.G.; Cetas, I.; Mccartney, S.; Cetas, J.; Dogan, A.; Fleseriu, M. Predictors of Silent Corticotroph Adenoma Recurrence; a Large Retrospective Single Center Study and Systematic Literature Review. Pituitary 2018, 21, 32–40. [Google Scholar] [CrossRef]
  30. Jahangiri, A.; Wagner, J.R.; Pekmezci, M.; Hiniker, A.; Chang, E.F.; Kunwar, S.; Blevins, L.; Aghi, M.K. A Comprehensive Long-Term Retrospective Analysis of Silent Corticotrophic Adenomas vs Hormone-Negative Adenomas. Neurosurgery 2013, 73, 8–17. [Google Scholar] [CrossRef]
  31. Strickland, B.A.; Shahrestani, S.; Briggs, R.G.; Jackanich, A.; Tavakol, S.; Hurth, K.; Shiroishi, M.S.; Liu, C.-S.J.; Carmichael, J.D.; Weiss, M.; et al. Silent Corticotroph Pituitary Adenomas: Clinical Characteristics, Long-Term Outcomes, and Management of Disease Recurrence. J. Neurosurg. 2021, 135, 1–8. [Google Scholar] [CrossRef]
  32. Machado, M.C.; Alcantara, A.E.E.; Pereira, A.C.L.; Cescato, V.A.S.; Castro Musolino, N.R.; de Mendonça, B.B.; Bronstein, M.D.; Fragoso, M.C.B.V. Negative Correlation between Tumour Size and Cortisol/ ACTH Ratios in Patients with Cushing’s Disease Harbouring Microadenomas or Macroadenomas. J. Endocrinol. Invest. 2016, 39, 1401–1409. [Google Scholar] [CrossRef] [PubMed]
  33. García-Martínez, A.; Fuentes-Fayos, A.C.; Fajardo, C.; Lamas, C.; Cámara, R.; López-Muñoz, B.; Aranda, I.; Luque, R.M.; Picó, A. Differential Expression of MicroRNAs in Silent and Functioning Corticotroph Tumors. J. Clin. Med. 2020, 9, 1838. [Google Scholar] [CrossRef] [PubMed]
  34. Huang, H.Y.; Lin, Y.C.D.; Li, J.; Huang, K.Y.; Shrestha, S.; Hong, H.C.; Tang, Y.; Chen, Y.G.; Jin, C.N.; Yu, Y.; et al. MiRTarBase 2020: Updates to the Experimentally Validated MicroRNA-Target Interaction Database. Nucleic Acids Res. 2020, 48, D148–D154. [Google Scholar] [CrossRef] [PubMed]
  35. Koning, A.S.C.A.M.; Buurstede, J.C.; van Weert, L.T.C.M.; Meijer, O.C. Glucocorticoid and Mineralocorticoid Receptors in the Brain: A Transcriptional Perspective. J. Endocr. Soc. 2019, 3, 1917–1930. [Google Scholar] [CrossRef]
  36. Nakai, Y.; Usui, T.; Tsukada, T.; Takahashi, H.; Fukata, J.; Fukushima, M.; Senoo, K.; Imura, H. Molecular Mechanisms of Glucocorticoid Inhibition of Human Proopiomelanocortin Gene Transcription. J. Steroid Biochem. Mol. Biol. 1991, 40, 301–306. [Google Scholar] [CrossRef]
  37. Drouin, J.; Charron, J.; Gagner, J.-P.; Jeannotte, L.; Nemer, M.; Plante, R.K.; Wrange, Ö. Pro-opiomelanocortin Gene: A Model for Negative Regulation of Transcription by Glucocorticoids. J. Cell. Biochem. 1987, 35, 293–304. [Google Scholar] [CrossRef]
  38. Huizenga, N.A.T.M.; De Lange, P.; Koper, J.W.; Clayton, R.N.; Farrell, W.E.; Van Der Lely, A.J.; Brinkmann, A.O.; De Jong, F.H.; Lamberts, S.W.J. Human Adrenocorticotropin-Secreting Pituitary Adenomas Show Frequent Loss of Heterozygosity at the Glucocorticoid Receptor Gene Locus. J. Clin. Endocrinol. Metab. 1998, 83, 917–921. [Google Scholar] [CrossRef]
  39. Miao, H.; Liu, Y.; Lu, L.; Gong, F.; Wang, L.; Duan, L.; Yao, Y.; Wang, R.; Chen, S.; Mao, X.; et al. Effect of 3 NR3C1 Mutations in the Pathogenesis of Pituitary ACTH Adenoma. Endocrinology 2021, 162, bqab167. [Google Scholar] [CrossRef]
  40. Mu, Y.M.; Takayanagi, R.; Imasaki, K.; Ohe, K.; Ikuyama, S.; Yanase, T.; Nawata, H. Low Level of Glucocorticoid Receptor Messenger Ribonucleic Acid in Pituitary Adenomas Manifesting Cushing’s Disease with Resistance to a High Dose-Dexamethasone Suppression Test. Clin. Endocrinol. 1998, 49, 301–306. [Google Scholar] [CrossRef]
  41. Ebisawa, T.; Tojo, K.; Tajima, N.; Kamio, M.; Oki, Y.; Ono, K.; Sasano, H. Immunohistochemical Analysis of 11-β-Hydroxysteroid Dehydrogenase Type 2 and Glucocorticoid Receptor in Subclinical Cushing’s Disease Due to Pituitary Macroadenoma. Endocr. Pathol. 2008, 19, 252–260. [Google Scholar] [CrossRef]
  42. Ledderose, C.; Möhnle, P.; Limbeck, E.; Schütz, S.; Weis, F.; Rink, J.; Briegel, J.; Kreth, S. Corticosteroid Resistance in Sepsis Is Influenced by MicroRNA-124-Induced Downregulation of Glucocorticoid Receptor-α. Crit. Care Med. 2012, 40, 2745–2753. [Google Scholar] [CrossRef] [PubMed]
  43. Spencer, R.L.; Deak, T. A Users Guide to HPA Axis Research. Physiol. Behav. 2017, 178, 43–65. [Google Scholar] [CrossRef] [PubMed]
  44. Han, F.; Ozawa, H.; Matsuda, K.I.; Nishi, M.; Kawata, M. Colocalization of Mineralocorticoid Receptor and Glucocorticoid Receptor in the Hippocampus and Hypothalamus. Neurosci. Res. 2005, 51, 371–381. [Google Scholar] [CrossRef]
  45. Castinetti, F.; Fassnacht, M.; Johanssen, S.; Terzolo, M.; Bouchard, P.; Chanson, P.; Cao, D.; Morange, I.; Picó, A.; Ouzounian, S.; et al. Merits and Pitfalls of Mifepristone in Cushing’s Syndrome. Eur. J. Endocrinol. 2009, 160, 1003–1010. [Google Scholar] [CrossRef] [PubMed]
  46. Fleseriu, M.; Laws, E.R.; Witek, P.; Pivonello, R.; Ferrigno, R.; de Martino, M.C.; Simeoli, C.; di Paola, N.; Pivonello, C.; Barba, L.; et al. Medical Treatment of Cushing’s Disease: An Overview of the Current and Recent Clinical Trials. Front. Endocrinol. 2020, 11, 648. [Google Scholar] [CrossRef]
  47. Senyuk, V.; Zhang, Y.; Liu, Y.; Ming, M.; Premanand, K.; Zhou, L.; Chen, P.; Chen, J.; Rowley, J.D.; Nucifora, G.; et al. Critical Role of MiR-9 in Myelopoiesis and EVI1-Induced Leukemogenesis. Proc. Natl. Acad. Sci. USA 2013, 110, 5594–5599. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, D.Z.; Chang, B.; Li, X.D.; Zhang, Q.H.; Zou, Y.H. MicroRNA-9 Promotes the Proliferation, Migration, and Invasion of Breast Cancer Cells via down-Regulating FOXO1. Clin. Transl. Oncol. 2017, 19, 1133–1140. [Google Scholar] [CrossRef] [PubMed]
  49. Guttilla, I.K.; White, B.A. Coordinate Regulation of FOXO1 by MiR-27a, MiR-96, and MiR-182 in Breast Cancer Cells. J. Biol. Chem. 2009, 284. [Google Scholar] [CrossRef]
  50. Xu, Y.; Zhao, S.; Cui, M.; Wang, Q. Down-Regulation of MicroRNA-135b Inhibited Growth of Cervical Cancer Cells by Targeting FOXO1. Int. J. Clin. Exp. Pathol. 2015, 8, 10294–10304. [Google Scholar]
  51. Lau, P.; Bossers, K.; Janky, R.; Salta, E.; Frigerio, C.S.; Barbash, S.; Rothman, R.; Sierksma, A.S.R.; Thathiah, A.; Greenberg, D.; et al. Alteration of the MicroRNA Network during the Progression of Alzheimer’s Disease. EMBO Mol. Med. 2013, 5, 1613–1634. [Google Scholar] [CrossRef]
  52. Kapali, J.; Kabat, B.E.; Schmidt, K.L.; Stallings, C.E.; Tippy, M.; Jung, D.O.; Edwards, B.S.; Nantie, L.B.; Raeztman, L.T.; Navratil, A.M.; et al. Foxo1 Is Required for Normal Somatotrope Differentiation. Endocrinology 2016, 157, 4351–4363. [Google Scholar] [CrossRef] [PubMed]
  53. Garrel, G.; Denoyelle, C.; L’Hôte, D.; Picard, J.Y.; Teixeira, J.; Kaiser, U.B.; Laverrière, J.N.; Cohen-Tannoudji, J. GnRH Transactivates Human AMH Receptor Gene via Egr1 and FOXO1 in Gonadotrope Cells. Neuroendocrinology 2019, 108, 65–83. [Google Scholar] [CrossRef] [PubMed]
  54. Skarra, D.V.; Arriola, D.J.; Benson, C.A.; Thackray, V.G. Forkhead Box O1 Is a Repressor of Basal and GnRH-Induced Fshb Transcription in Gonadotropes. Mol. Endocrinol. 2013, 27, 1825–1839. [Google Scholar] [CrossRef] [PubMed]
  55. Xiao, Z.; Wang, Z.; Hu, B.; Mao, Z.; Zhu, D.; Feng, Y.; Zhu, Y. MiR-1299 Promotes the Synthesis and Secretion of Prolactin by Inhibiting FOXO1 Expression in Drug-Resistant Prolactinomas. Biochem. Biophys. Res. Commun. 2019, 520, 79–85. [Google Scholar] [CrossRef]
  56. Benite-Ribeiro, S.A.; de Lima Rodrigues, V.A.; Machado, M.R.F. Food Intake in Early Life and Epigenetic Modifications of Pro-Opiomelanocortin Expression in Arcuate Nucleus. Mol. Biol. Rep. 2021, 48, 3773–3784. [Google Scholar] [CrossRef]
  57. Knosp, E.; Steiner, E.; Kitz, K.; Matula, C.; Parent, A.D.; Laws, E.R.; Ciric, I. Pituitary Adenomas with Invasion of the Cavernous Sinus Space: A Magnetic Resonance Imaging Classification Compared with Surgical Findings. Neurosurgery 1993, 33, 610–618. [Google Scholar] [CrossRef]
  58. Bujko, M.; Kober, P.; Boresowicz, J.; Rusetska, N.; Zeber-Lubecka, N.; Paziewska, A.; Pekul, M.; Zielinski, G.; Styk, A.; Kunicki, J.; et al. Differential MicroRNA Expression in USP8-Mutated and Wild-Type Corticotroph Pituitary Tumors Reflect the Difference in Protein Ubiquitination Processes. J. Clin. Med. 2021, 10, 375. [Google Scholar] [CrossRef]
  59. Mi, H.; Ebert, D.; Muruganujan, A.; Mills, C.; Albou, L.P.; Mushayamaha, T.; Thomas, P.D. PANTHER Version 16: A Revised Family Classification, Tree-Based Classification Tool, Enhancer Regions and Extensive API. Nucleic Acids Res. 2021, 49, D394–D403. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Covid-19 and Cushing’s Disease in a Patient with ACTH-secreting Pituitary Carcinoma

Abstract

Summary

The pandemic caused by severe acute respiratory syndrome coronavirus 2 is of an unprecedented magnitude and has made it challenging to properly treat patients with urgent or rare endocrine disorders. Little is known about the risk of coronavirus disease 2019 (COVID-19) in patients with rare endocrine malignancies, such as pituitary carcinoma. We describe the case of a 43-year-old patient with adrenocorticotrophic hormone-secreting pituitary carcinoma who developed a severe COVID-19 infection. He had stabilized Cushing’s disease after multiple lines of treatment and was currently receiving maintenance immunotherapy with nivolumab (240 mg every 2 weeks) and steroidogenesis inhibition with ketoconazole (800 mg daily). On admission, he was urgently intubated for respiratory exhaustion. Supplementation of corticosteroid requirements consisted of high-dose dexamethasone, in analogy with the RECOVERY trial, followed by the reintroduction of ketoconazole under the coverage of a hydrocortisone stress regimen, which was continued at a dose depending on the current level of stress. He had a prolonged and complicated stay at the intensive care unit but was eventually discharged and able to continue his rehabilitation. The case points out that multiple risk factors for severe COVID-19 are present in patients with Cushing’s syndrome. ‘Block-replacement’ therapy with suppression of endogenous steroidogenesis and supplementation of corticosteroid requirements might be preferred in this patient population.

Learning points

  • Comorbidities for severe coronavirus disease 2019 (COVID-19) are frequently present in patients with Cushing’s syndrome.
  • ‘Block-replacement’ with suppression of endogenous steroidogenesis and supplementation of corticosteroid requirements might be preferred to reduce the need for biochemical monitoring and avoid adrenal insufficiency.
  • The optimal corticosteroid dose/choice for COVID-19 is unclear, especially in patients with endogenous glucocorticoid excess.
  • First-line surgery vs initial disease control with steroidogenesis inhibitors for Cushing’s disease should be discussed depending on the current healthcare situation.

Background

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had a significant impact on the health care systems to date. The clinical presentation of coronavirus disease 2019 (COVID-19) is diverse, ranging from asymptomatic illness to respiratory failure requiring admission to the intensive care unit (ICU). Risk factors for severe course include old age, male gender, comorbidities such as arterial hypertension, diabetes mellitus, chronic lung-, heart-, liver- and kidney disease, malignancy, immunodeficiency and pregnancy (1). Little is known about the risk of COVID-19 in patients with rare endocrine malignancies, such as pituitary carcinoma.

Case presentation

This case concerns a 43-year-old man with adrenocorticotrophic hormone (ACTH)-secreting pituitary carcinoma (with cerebellar and cervical drop metastases) with a severe COVID-19 infection. He had previously received multiple treatment modalities including surgery, radiotherapy, ketoconazole, pasireotide, cabergoline, bilateral (subtotal) adrenalectomy and temozolomide chemotherapy as described elsewhere (2). His most recent therapy was a combination of immune checkpoint inhibitors consisting of ipilimumab (3 mg/kg) and nivolumab (1 mg/kg) (anti-CTLA-4 and anti-PD-1, respectively) every 3 weeks for four cycles, after which maintenance therapy with nivolumab (240 mg) every 2 weeks was continued. Residual endogenous cortisol production was inhibited with ketoconazole 800 mg daily. He had stabilized disease with a decrease in plasma ACTH, urinary free cortisol and stable radiological findings (2). Surgical resection of the left adrenal remnant was planned but was not carried out due to the development of a COVID-19 infection.

In March 2021, he consulted our emergency department for severe respiratory complaints. He had been suffering from upper respiratory tract symptoms for one week, with progressive dyspnoea in the last three days. He tested positive for SARS-CoV-2 the day before admission. On examination, his O2 saturation was 72%, with tachypnoea (40/min) and bilateral pulmonary crepitations. His temperature was 37.2°C, blood pressure 124/86 mmHg and pulse rate 112 bpm. High-flow oxygen therapy was initiated but yielded insufficient improvement (O2 saturation of 89% and tachypnoea 35/min). He was urgently intubated for respiratory exhaustion.

Investigation

Initial investigations showed type 1 respiratory insufficiency with PaO2 of 52.5 mmHg (normal 75–90), PaCO2 of 33.0 mmHg (normal 36–44), pH of 7.47 (normal 7.35–7.45) and a P/F ratio of 65.7 (normal >300). His inflammatory parameters were elevated with C-reactive protein level of 275.7 mg/L (normal <5·0) and white blood cell count of 7.1 × 10⁹ per L with 72.3% neutrophils. His most recent morning plasma ACTH-cortisol level (measured using the Elecsys electrochemiluminescence immunoassays on a Cobas 8000 immunoanalyzer [Roche Diagnostics]) before his admission was 213 ng/L (normal 7.2–63) and 195 µg/L (normal 62–180) respectively, while a repeat measurement 3 weeks after his admission demonstrated increased cortisol levels of 547 µg/L (possibly iatrogenic due to treatment with high-dose hydrocortisone) and a decreased ACTH of 130 ng/L.

Treatment

On admission, he was started on high-dose dexamethasone therapy for 10 days together with broad-spectrum antibiotics for positive sputum cultures containing Serratia, methicillin-susceptible Staphylococcus aureus and Haemophilus influenzae. Thromboprophylaxis with an intermediate dose of low molecular weight heparin (tinzaparin 14 000 units daily for a body weight of 119 kg) was initiated. A ‘block-replacement’ regimen was adopted with the continuation of ketoconazole (restarted on day 11) in view of his endocrine treatment and the supplementation of hydrocortisone at a dose depending on the current level of stress. The consecutive daily dose of hydrocortisone and ketoconazole is shown in Fig. 1.

Figure 1View Full Size
Figure 1
‘Block-replacement’ therapy with ketoconazole and hydrocortisone/dexamethasone. Dexamethasone 10 mg daily was initially started as COVID-19 treatment, followed by hydrocortisone at a dose consistent with current levels of stress. Ketoconazole was restarted on day 11 and titrated to a dose of 800 mg daily to suppress endogenous glucocorticoid production.

Citation: Endocrinology, Diabetes & Metabolism Case Reports 2022, 1; 10.1530/EDM-21-0182

Outcome and follow-up

He developed multiple organ involvement, including metabolic acidosis, acute renal failure requiring continuous venovenous hemofiltration, acute coronary syndrome type 2, septic thrombophlebitis of the right jugular vein, and critical illness polyneuropathy. He was readmitted twice to the ICU, for ventilator-associated pneumonia and central line-associated bloodstream infection respectively. He eventually recovered and was discharged from the hospital to continue his rehabilitation.

Discussion

We describe the case of a patient with severe COVID-19 infection with active Cushing’s disease due to pituitary carcinoma, who was treated with high-dose dexamethasone followed by ‘block-replacement’ therapy with hydrocortisone in combination with off-label use of ketoconazole as a steroidogenesis inhibitor. His hospitalization was prolonged by multiple readmissions to the ICU for infectious causes. Our case illustrates the presence of multiple comorbidities for a severe and complicated course of COVID-19 in a patient with active Cushing’s disease.

Dexamethasone was initially chosen as the preferred corticosteroid therapy, in analogy with the RECOVERY trial, in which dexamethasone at a dose of 6mg once daily (oral or i.v.) resulted in lower 28-day mortality in hospitalized patients with COVID-19 requiring oxygen therapy or invasive mechanical ventilation (3). However, the optimal dose/choice of corticosteroid therapy is unclear, especially in a patient population with pre-existing hypercortisolaemia. A similar survival benefit for hydrocortisone compared to dexamethasone has yet to be convincingly demonstrated. This may be explained by differences in anti-inflammatory activity but could also be due to the fact that recent studies with hydrocortisone were stopped early and were underpowered (45).

Multiple risk factors for a complicated course of COVID-19 are present in patients with Cushing’s syndrome and might increase morbidity and mortality (67). These include a history of obesity, arterial hypertension and impaired glucose metabolism. Prevention and treatment of these pre-existing comorbidities are essential.

Patients with Cushing’s syndrome also have an increased thromboembolic risk, which is further accentuated by the development of severe COVID-19 infection (67). Thromboprophylaxis with low molecular weight heparin is associated with lower mortality in COVID-19 patients with high sepsis‐induced coagulopathy score or high D-dimer levels (8) and is presently widely used in the treatment of severe COVID-19 disease (9). Subsequently, this treatment is indicated in hospitalized COVID-19 patients with Cushing’s syndrome. It is unclear whether therapeutic anticoagulation dosing could provide additional benefits (67). An algorithm based on the International Society on Thrombosis and Hemostasis-Disseminated Intravascular Coagulation score was proposed to evaluate the ideal anticoagulation therapy in severe/critical COVID-19 patients, with an indication for therapeutic low molecular weight heparin dose at a score ≥5 (9).

Furthermore, the chronic cortisol excess induces suppression of the innate and adaptive immune response. Patients with Cushing’s syndrome, especially when severe and active, should be considered immunocompromised and have increased susceptibility for viral and other (hospital-acquired) infections. Prophylaxis for Pneumocystis jirovecii with trimethoprim/sulfamethoxazole should therefore be considered (67).

Additionally, there is a particular link between the pathophysiology of COVID-19 and Cushing’s syndrome. The SARS-CoV-2 virus (as well as other coronaviruses) enter human cells by binding the ACE2 receptor. The transmembrane serine protease 2 (TMPRSS2), expressed by endothelial cells, is additionally required for the priming of the spike-protein of SARS-CoV-2, leading to viral entry. TMPRSS2 was studied in prostate cancer and found to be regulated by androgen signalling. Consequently, the androgen excess frequently associated with Cushing’s syndrome might be an additional risk factor for contracting COVID-19 via higher TMPRSS2 expression (10), especially in women, in whom the effect of excess androgen would be more noticeable compared to male patients with Cushing’s syndrome.

Treating Cushing’s syndrome with a ‘block-replacement’ approach, with suppression of endogenous steroidogenesis and supplementation of corticosteroid requirements, is an approach that should be considered, especially in severe or cyclic disease. The use of this method might decrease the need for monitoring and reduce the occurrence of adrenal insufficiency (7). Our patient was on treatment with ketoconazole, which was interrupted at initial presentation and then restarted under the coverage of a hydrocortisone stress regimen. Ketoconazole was chosen because of its availability. Advantages of ketoconazole over metyrapone include its antifungal activity with the potential for prevention of invasive pulmonary fungal infections, as well as its antiandrogen action (especially in female patients) and subsequent inhibition of TMPRSS2 expression (10). Regular monitoring of the liver function (every month for the first 3 months, at therapy initiation or dose increase) is necessary. Caution is needed due to its inhibition of multiple cytochrome P450 enzymes (including CYP3A4) and subsequently greater risk of drug-drug interactions vs metyrapone (710). Another disadvantage of ketoconazole is the need for oral administration. In our patient, ketoconazole was delivered through a nasogastric tube. i.v. etomidate is an alternative in case of an unavailable enteral route.

Finally, as a general point, the first-line treatment of a patient with a novel diagnosis of Cushing’s disease is transsphenoidal surgery. Recent endocrine recommendations pointed out the possibility of initial disease control with steroidogenesis inhibitors in patients without an indication for urgent intervention during a high prevalence of COVID-19 (7). This would allow the optimalization of metabolic parameters; emphasizing that the short-to mid-term prognosis is related to the cortisol excess and not its cause. Surgery could then be postponed until the health situation allows for safe elective surgery (7). This decision depends of course on the evolution of COVID-19 and the healthcare system in each country and should be closely monitored by policymakers and physicians.

Declaration of interest

The authors declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.

Funding

This work did not receive any specific grant from any funding agency in the public, commercial, or not-for-profit sector.

Patient consent

Written informed consent for publication of their clinical details and/or clinical images was obtained from the patient.

Author contribution statement

J M K de Filette is an endocrinologist-in-training and was the main author. All authors were involved in the clinical care of the patient. All authors contributed to the reviewing and editing process and approved the final version of the manuscript.

References

Vitamin D Deficiency in Cushing’s Disease: Before and After Its Supplementation

1
Department of Health Promotion, Maternal-Infantile Care, Excellence Internal and Specialist Medicine “G. D’Alessandro” [PROMISE], Section of Endocrine Disease and Nutrition, University of Palermo, 90127 Palermo, Italy
2
Biochemistry Head CQRC Division (Quality Control and Biochemical Risk), Department of Health Promotion, Maternal-Infantile Care, Excellence Internal and Specialist Medicine “G. D’Alessandro” [PROMISE], University of Palermo, 90127 Palermo, Italy
Author to whom correspondence should be addressed.
Academic Editor: Edgard Delvin
Nutrients 202214(5), 973; https://doi.org/10.3390/nu14050973

Abstract

Background: The primary objective of the study was to assess serum 25-hydroxyvitamin D [25(OH)D] values in patients with Cushing’s disease (CD), compared to controls. The secondary objective was to assess the response to a load of 150,000 U of cholecalciferol. Methods: In 50 patients with active CD and 48 controls, we evaluated the anthropometric and biochemical parameters, including insulin sensitivity estimation by the homeostatic model of insulin resistance, Matsuda Index and oral disposition index at baseline and in patients with CD also after 6 weeks of cholecalciferol supplementation. Results: At baseline, patients with CD showed a higher frequency of hypovitaminosis deficiency (p = 0.001) and lower serum 25(OH)D (p < 0.001) than the controls. Six weeks after cholecalciferol treatment, patients with CD had increased serum calcium (p = 0.017), 25(OH)D (p < 0.001), ISI-Matsuda (p = 0.035), oral disposition index (p = 0.045) and decreased serum PTH (p = 0.004) and total cholesterol (p = 0.017) values than at baseline. Multivariate analysis showed that mean urinary free cortisol (mUFC) was independently negatively correlated with serum 25(OH)D in CD. Conclusions: Serum 25(OH)D levels are lower in patients with CD compared to the controls. Vitamin D deficiency is correlated with mUFC and values of mUFC > 240 nmol/24 h are associated with hypovitaminosis D. Cholecalciferol supplementation had a positive impact on insulin sensitivity and lipids.

1. Introduction

Vitamin D is the precursor of a hormone with pleiotropic effects. Its deficiency has been largely investigated and shown to be associated with many complications including diabetes mellitus, adrenal insufficiency, cardiovascular disease, neurological disorders and other endocrinopathies [1,2,3].
Vitamin D, also known as cholecalciferol, is first formed in the skin by the photolysis of 7-dehydrocholesterol and after hydroxylated in the liver to 25-hydroxyvitamin D [25(OH)D]. It is further transformed in the kidney into 1,25-dihydroxyvitamin D3 (1,25(OH)2D3) (calcitriol) that is the active form [4].
Cushing’s disease (CD) is characterized by a cortisol excess due to autonomous pituitary ACTH secretion. Patients with CD show many comorbidities such as cardiovascular disease, metabolic disease, diabetes mellitus, metabolic syndrome, dyslipidemia, obesity, osteoporosis/osteopenia and infections that contribute to increasing the mortality risk for these patients [5,6,7,8,9,10,11]. Indeed, GCs are key regulators of intermediary metabolism promoting hepatic gluconeogenesis and glycogenosis and on lipid metabolism favouring the deposition of fat to the upper trunk and the face [12]. They stimulate water diuresis, glomerular filtration rate and renal plasma flow and these effects result in arterial hypertension and atherosclerosis. GCs reduce bone remodelling, augment urinary calcium excretion and decrease the intestinal calcium absorption. In addition, they act on immune and hematological systems inhibiting the secretion of interleukins and increasing the red blood cell count, respectively [12].
An interesting relationship exists between glucocorticoids (GCs) and vitamin D values [13,14,15,16]. Indeed, exogenous steroid therapy has been reported to be associated with vitamin deficiency [13]. The mechanism by which GCs reduce 25(OH)D levels is not direct, but indirect, regulating vitamin D receptor expression in many tissues and cells [17,18]. Some authors have shown that treatment with dexamethasone in mice was associated with a decrease in 1α-hydroxylase which is involved in the conversion from 25(OH)D3 to the active metabolite 1,25(OH)2D3 and an increase in 24-hydroxylase, able to break down the active form of calcitriol, in inactive, reducing circulating 25(OH)D levels [19]. In a clinical setting, controversial data have been reported on GCs effects on serum 1,25(OH)2D concentrations [20,21,22,23]. A likely reason for these discrepancies might be the marked heterogeneity of the studied groups. Some of these studies were performed in humans [23,24,25,26], and others in animal models [27,28], but only a few studies were conducted in subjects with endogenous hypercortisolism.
Low serum 25(OH)D levels have significant skeletal and extra-skeletal consequences such as myopathy, high risk of fractures and also affect the immune system and metabolism. All of these systems are impaired in patients with hypercortisolism and a vitamin D deficiency may provide a further aggravation of CD comorbidities. Indeed, it may cause a reduced intestinal calcium absorption resulting in secondary hypocalcemia and hyperparathyroidism leading to a bone demineralization. Its deficiency can contribute to obesity and metabolic syndrome due to the lack of antiadipogenic effect of vitamin D and to cardiovascular disease by a deregulation of the renin–angiotensin–aldosterone system, cardiac contractility and increase in cytokine release [29]. In the end, vitamin D deficiency causes impaired insulin sensitivity and immune system [30].
The discrepancies that emerge in the above-mentioned studies suggest a need to investigate the role of 25(OH)D in patients with CD. Therefore, the primary objective of the study was to evaluate serum 25(OH)D levels in patients with CD, compared to a control group matched for age, BMI and gender, and search for a possible correlation with the degree of hypercortisolism. The secondary objective was to evaluate the response to a course of 150,000 U of cholecalciferol on metabolic and hormonal parameters 6 weeks after the administration in patients with CD.

2. Materials and Methods

2.1. Subjects and Study Design

Fifty patients with active CD, 43 of them women (86%) and 7 of them men (20%) (mean age 50.9 ± 17.4 years; mean duration of disease 32.5 ± 22.4 years), followed from January 2016 to December 2020, by the Endocrinology of the University of Palermo, were included in the current study. Clinical practice guidelines and a recent consensus statement were used to diagnose CD [31,32].
We recruited a control group matched for age, BMI and gender in the same temporal period. It was composed of 48 patients, 33 women (82.5%) and 7 men (17.5%) (mean age 48.5 ± 13.4 years) were evaluated by our team for a suspicion not biochemically confirmed of Cushing’s syndrome (CS).
In all patients, we evaluated phenotypic characteristics including moon face, facial rubor, dorsal fat pad or buffalo hump, defined as a fatty tissue deposit between the shoulders, purple striae, defined as wide, reddish-purple streaks, and myopathy defined as muscle weakness at the proximal level.
We also assessed cardiovascular, metabolic and bone comorbidities. The diagnosis of metabolic syndrome was based on National Cholesterol Education Program Adult Treatment Panel (NCEP ATP III) criteria, while the diagnosis of diabetes mellitus and prediabetes were based on the American Diabetes Association (ADA, Arlington, VA, USA) criteria [33,34].
Among patients with diabetes mellitus (18 out of 50), 16 were treated with metformin alone, while 2 were treated with a combination of metformin and GLP-1 agonist receptors. Metformin and GLP-1 agonist receptors were discontinued 24 h and 2 weeks before metabolic evaluations, respectively, to avoid any interference with metabolic parameters. Diabetic patients were on good metabolic control (HbA1c ≤ 7%). Both CD patients and the controls were naïve to cholecalciferol.
In CD and the controls, BMI and waist circumference (WC), fasting serum lipids (total cholesterol (TC), HDL cholesterol, LDL cholesterol and triglycerides (TG), HbA1c, glycaemia, insulinaemia, albumin corrected calcium, phosphorus and parathyroid hormone (PTH) were assessed. To avoid seasonal influences, serum 25(OH)D levels were only assayed between winter and spring seasons (November–April). We evaluated urinary free cortisol (UFC) as the mean of three 24 h urine collections (mUFC), cortisol after a low dose of dexamethasone suppression test and plasma ACTH. We defined patients with mild hypercortisolism when mUFC levels not exceeding twice the upper limit of normal (ULN), moderate hypercortisolism by a level of mUFC more than 2 to 5 times the ULN and severe hypercortisolism by a mUFC level more than 5 times the ULN, as previously reported [35].
As defined by the Endocrine Society guidelines, we considered 25(OH)D deficiency for values < 20 ng/mL (50 nmol/L), insufficiency as levels of 20–30 ng/mL (50–75 nmol/L) and sufficiency for values ≥ 30 ng/mL (≥75 nmol/L) [36]. In addition, severe 25(OH)D deficiency was defined by levels < 10 ng/mL (<25 nmol/L) [37].
As markers of insulin sensitivity, we calculated the homeostatic model of insulin resistance (HOMA2-IR) [38], and in 32 patients with CD and in 40 controls who had no previous diagnosis of diabetes, we also evaluated the Matsuda index of insulin sensitivity (ISI-Matsuda) [39], the oral disposition index (DIo) [40] and the area under the curve for insulin (AUC2h insulinemia) and glucose (AUC2h glycaemia).
At the baseline visit, we assessed patients’ lifestyle habits: physical activity level, balanced diet (consumption of dairy products, meat, coffee, soft drinks), exposure to ultraviolet (UV) radiation, smoking status and alcohol use.
We excluded patients with adrenal-dependent hypercortisolism, pregnancy, taking oral contraceptives, liver or renal disease, cholecalciferol supplementation within 3 months before the study, malabsorption syndrome and exposure to ultraviolet (UV) radiation (solarium and sunscreen usage).
Patients with CD received an oral load dose of cholecalciferol of 150,000 UI [41,42] and biochemical parameters (metabolic and hormonal) were assayed 6 weeks after administration.
The study protocol was approved by the Ethics Committee of the Policlinico Paolo Giaccone hospital. All patients signed a written informed consent.

2.2. Assays

Biochemical parameters were measured by standard methods (Modular P800, Roche, Milan, Italy), as previously reported [9].
Hormonal parameters were measured by electrochemiluminescence immunoassay (ECLIA, Elecsys, Roche, Milan, Italy) following the manufacturer’s instructions, as previously reported [9].
Mean UFC was measured by mass spectrometry, as previously reported [35].
Normal values for hormonal markers were defined as follows: ACTH 2.2–14 pmol/L and UFC 59–378 nmol/24 h.

2.3. Statistical Analysis

We used statistical Packages for Social Science SPSS version 19 (SPSS, Inc., Chicago, IL, USA) for data analysis. The normality of quantitative variables was tested with the Shapiro–Wilk test. We calculated mean ± SD for continuous variables and rates and proportions for categorical variables. The differences between paired continuous variables (CD vs. controls) were analysed using one-way ANOVA. We used univariate Pearson correlation to evaluate the relations with the outcome parameters. For those variables which were significant at univariate correlation, we performed multiple linear regression analysis to identify independent predictors of the dependent variable 25(OH)D. A p-value of 0.05 was considered statistically significant. A receiver operating characteristic (ROC) analysis was performed to investigate the diagnostic ability of significantly associated risk factors to predict 25(OH)D deficiency. The ROC curve is plotted as sensitivity versus 1-specificity. The area under the ROC curve (AUC) was estimated to measure the overall performance of the predictive factors for serum 25(OH)D deficiency.

3. Results

At baseline, patients with CD had a higher frequency of arterial hypertension (p = 0.009), osteoporosis/osteopenia (p = 0.002), hypercholesterolemia (p = 0.002), diabetes mellitus (p = 0.026), myopathy (p < 0.001), facial rubor (p = 0.005), buffalo hump (p = 0.002) and hypovitaminosis deficiency (p = 0.001) than the controls (Table 1).
Table 1. Comorbidities of patients with CD and controls at baseline.
Table
By contrast, the controls had a higher frequency of vitamin D sufficiency (p = 0.004). Patients with CD also had higher WC (p = 0.031), PTH (p = 0.003), glycaemia (p = 0.010), HbA1c (p = 0.004), total cholesterol (p < 0.001), LDL cholesterol (p = 0.002), ACTH (p < 0.001), mUFC (p = 0.001), cortisol after a low dose of dexamethasone suppression test (p = 0.001) and lower 25(OH)D (p < 0.001), ISI-Matsuda (p = 0.007) and DIo (p = 0.003) than the controls (Table 2).
Table 2. Anthropometric and biochemical parameters of patients with CD and controls at baseline.
Table
Six weeks after cholecalciferol treatment, patients with CD showed increased serum calcium (p = 0.017), 25(OH)D (p < 0.001), ISI-Matsuda (p = 0.035), DIo (p = 0.045) and a decrease in PTH (p = 0.004) and total cholesterol (p = 0.017) levels than at baseline (Table 3).
Table 3. Anthropometric and biochemical parameters at baseline and 6 weeks after cholecalciferol supplementation in patients with CD.
Table
Considering the degree of hypercortisolism, in patients with severe hypercortisolism we observed 25(OH)D deficiency in 73.1% of cases (53.8% of them had a severe deficiency), insufficiency in 12.5% of cases and sufficiency in 6.3% of cases. In patients with moderate hypercortisolism, we observed 25(OH)D deficiency in 64.7% of cases (29% of them had a severe deficiency), insufficiency in 23.5% of cases and sufficiency in 11.8% of cases. In patients with mild hypercortisolism, we observed deficiency in 52.9% of cases (20% of them had a severe deficiency), insufficiency in 41.1% of cases and sufficiency in 6% of cases.
At univariate correlation, in patients with CD at baseline, serum 25(OH)D was inversely correlated with glycaemia (r = −0.385, p = 0.019), HbA1c (r = −0.391, p = 0.017), WC (r = −0.373, p = 0.023), mUFC (r = −0.466, p = 0.033) and cortisol after a low dose of dexamethasone suppression test (r = −0.299, p = 0.049) (Table 4). In the controls, at baseline, 25(OH)D was inversely correlated with WC (r = −0.130, p = 0.042) (Table 4).
Table 4. Correlation of serum 25-hydroxyvitamin D [25(OH)D] levels at baseline in patients with Cushing’s disease and controls.
Table
Multivariate analysis showed that mUFC was independently inversely associated with 25(OH)D (p = 0.010) in patients with CD (Figure 1). In the controls, no significant associations were found.
Nutrients 14 00973 g001 550
Figure 1. Independent variables associated with serum 25(OH)D in patients with active CD at multivariate analysis. mUFC: mean urinary free cortisol.
The ROC analysis showed that a cut-off of mUFC > 240 nmol/24 h was associated with 25(OH)D deficiency with a specificity of 100% and a sensitivity of 56.9%, AUC 0.803 (Figure 2).
Nutrients 14 00973 g002 550
Figure 2. 25(OH)D status and mUFC. ROC curve showed that a cut-off of mUFC > 240 nmol/24 h could be associated with 25(OH)D deficiency. Statistical analysis was performed using the chi-square test and receiver operator characteristic (ROC) curve analysis.

4. Discussion

The present study shows that patients with active CD have lower serum 25(OH)D values than the controls and that serum 25(OH)D levels are inversely correlated with mUFC in CD. In addition, a cholecalciferol load is associated after 6 weeks from the administration with an improvement of serum 25(OH)D and glycometabolic and lipid parameters in patients with CD. Furthermore, we found that higher values of mUFC than 240 nmol/24 h are predictive of 25(OH)D deficiency. The degree of hypercortisolism evaluated by UFC levels is a useful parameter to quantify the “amount” of cortisol secretion, even though it is not sufficiently exhaustive to assess the aggressiveness of the disease [35]. Indeed, a combination of several factors, including the degree of hypercortisolism, but also the duration of the disease, age and other individual predisposing factors, contribute to the aggressiveness of the disease.
Long-standing studies were conducted on vitamin D levels in patients with CD. Patients with CD, with and without osteopenia, were compared before and after oral calcium load showing that serum 1,25 (OH)2D3 plasma levels were higher in subjects with osteopenia than in those without it, likely due to a secondary increase in PTH levels as an effect of hypercortisolism [19]. Another study investigated the effect of hypercortisolism and eucortisolism, showing a reduction in serum 25(OH)D levels, but not in 1,25 (OH)2D3 in patients with hypercortisolism. By contrast, two other studies found normal serum 25(OH)D values in patients with CD [23,24]. However, all the above-mentioned studies were conducted on a small sample of patients. Recently, a meta-analysis conducted on the studies that evaluated serum 25(OH)D levels in patients treated with GCs reported lower serum 25(OH)D levels in these patients compared to healthy subjects [16]. A hypothetical reason was that patients with CD had low 24-hydroxylase levels than the controls, causing an alteration of vitamin D catabolism.
An interesting in vitro study in NCI-H295R cells found that treatment with 1,25(OH)2D3 decreased corticosterone secretion without affecting cortisol levels [43].
As expected, in the current study, we showed that treatment with cholecalciferol is associated with an improvement in insulin sensitivity and total cholesterol values in patients with CD. Indeed, cholecalciferol supplementation has been reported to be associated with improved peripheral insulin sensitivity and secretion in patients at high risk of diabetes or with type 2 diabetes [44]. A recent meta-analysis on 41 randomized controlled studies showed a significant improvement in total cholesterol levels after cholecalciferol supplementation. In addition, this improvement was more pronounced in patients with vitamin D deficiency [45,46].
A recent study compared the metabolism of vitamin D in patients with CD and controls after cholecalciferol treatment, showing that patients with CD had a higher 25(OH)D/24,25(OH)2D ratio than healthy controls, likely due to a decrease in 24-hydroxylase activity. The authors concluded that this alteration of vitamin D catabolism might have an influence on the effectiveness of cholecalciferol therapy in CD [47].
There are some limitations in the current study. First, the study is not randomized. Second, the dose of cholecalciferol administered is the same independently of the baseline serum 25(OH)D values. Third, we did not register the intake of milk and dairy products of the patients included in the study.
In conclusion, serum 25(OH)D levels are lower in subjects with active CD compared to controls matched for age, BMI and gender. Vitamin D deficiency is correlated with mUFC and values of mUFC > 240 nmol/24 h are predictive of 25(OH)D deficiency. In addition, cholecalciferol supplementation has a positive impact on insulin sensitivity and lipids and therefore should be considered part of the treatment of patients with CD at diagnosis, in order to improve the comorbidities. However, further studies are needed to evaluate a possible effect of cholecalciferol supplementation on the aggressiveness of CD.

Author Contributions

Conceptualization, V.G. and F.D.G.; methodology, V.G.; software, V.G.; validation, V.G., F.D.G. and C.G.; formal analysis, V.G.; investigation, V.G.; resources, F.D.G.; data curation, V.G.; writing—original draft preparation, V.G.; writing—review and editing, V.G.; visualization, V.G.; supervision, C.G.; project administration, C.G.; funding acquisition, C.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and was approved by the Institutional Review Board (or Ethics Committee) of Policlinico Paolo Giaccone (number 1, approved on the 17 January 2022).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

Data are available on demand at corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Muscogiuri, G.; Altieri, B.; Annweiler, C.; Balercia, G.; Pal, H.B.; Boucher, B.J.; Cannell, J.J.; Foresta, C.; Grübler, M.R.; Kotsa, K.; et al. Vitamin D and chronic diseases: The current state of the art. Arch. Toxicol. 2017, 91, 97–107. [Google Scholar] [CrossRef] [PubMed]
  2. Marino, R.; Misra, M. Extra-skeletal effects of Vitamin D. Nutrients 2019, 11, 1460. [Google Scholar] [CrossRef] [PubMed]
  3. Zendehdel, A.; Arefi, M. Molecular evidence of role of vitamin D deficiency in various extraskeletal diseases. J. Cell. Biochem. 2019, 120, 8829–8840. [Google Scholar] [CrossRef] [PubMed]
  4. Bikle, D.; Christakos, S. New aspects of vitamin D metabolism and action-addressing the skin as source and target. Nat. Rev. Endocrinol. 2020, 16, 234–252. [Google Scholar] [CrossRef]
  5. Pivonello, R.; Isidori, A.; De Martino, M.C.; Newell-Price, J.; Biller, B.M.K.; Colao, A. Complications of Cushing’s syndrome: State of the art. Lancet Diabetes Endocrinol. 2016, 4, 611–629. [Google Scholar] [CrossRef]
  6. Guarnotta, V.; Ferrigno, R.; Martino, M.; Barbot, M.; Isidori, A.M.; Scaroni, C.; Ferrante, A.; Arnaldi, G.; Pivonello, R.; Giordano, C. Glucocorticoid excess and COVID-19 disease. Rev. Endocr. Metab. Disord. 2020, 22, 703–714. [Google Scholar] [CrossRef]
  7. Giordano, C.; Guarnotta, V.; Pivonello, R.; Amato, M.C.; Simeoli, C.; Ciresi, A.; Cozzolino, A.; Colao, A. Is diabetes in Cushing’s syndrome only a consequence of hypercortisolism? Eur. J. Endocrinol. 2014, 170, 311–319. [Google Scholar] [CrossRef]
  8. Drey, M.; Berr, C.M.; Reincke, M.; Fazel, J.; Seissler, J.; Schopohl, J.; Bidlingmaier, M.; Zopp, S.; Reisch, N.; Beuschlein, F.; et al. Cushing′s syndrome: A model for sarcopenic obesity. Endocrine 2017, 57, 481–485. [Google Scholar] [CrossRef]
  9. Guarnotta, V.; Prinzi, A.; Pitrone, M.; Pizzolanti, G.; Giordano, C. Circulating irisin levels as a marker of osteosarcopenic-obesity in Cushing’s disease. Diabetes Metab. Syndr. Obes. 2020, 13, 1565–1574. [Google Scholar] [CrossRef]
  10. Hakami, O.A.; Ahmed, S.; Karavitaki, N. Epidemiology and mortality of Cushing′s syndrome. Best Pr. Res. Clin. Endocrinol. Metab. 2021, 35, 101521. [Google Scholar] [CrossRef]
  11. Javanmard, P.; Duan, D.; Geer, E.B. Mortality in patients with endogenous Cushing′s Syndrome. Endocrinol. Metab. Clin. North Am. 2018, 47, 313–333. [Google Scholar] [CrossRef] [PubMed]
  12. McKay, L.I.; Cidlowski, J.A. Physiologic and pharmacologic effects of corticosteroids. In Holland-Frei Cancer Medicine, 6th ed.; Kufe, D.W., Pollock., R.E., Weichselbaum, R.R., Eds.; BC Decker: Hamilton, ON, Canada, 2003. [Google Scholar]
  13. Tirabassi, G.; Salvio, G.; Altieri, B.; Ronchi, C.L.; Della Casa, S.; Pontecorvi, A.; Balercia, G. Adrenal disorders: Is there any role for Vitamin D? Rev. Endocr. Metab. Disord. 2016, 18, 355–362. [Google Scholar] [CrossRef] [PubMed]
  14. Skversky, A.L.; Kumar, J.; Abramowitz, M.K.; Kaskel, F.J.; Melamed, M.L. Association of glucocorticoid use and low 25-Hydroxyvitamin D levels: Results from the National Health and Nutrition Examination Survey (NHANES): 2001–2006. J. Clin. Endocrinol. Metab. 2011, 96, 3838–3845. [Google Scholar] [CrossRef] [PubMed]
  15. Muscogiuri, G.; Altieri, B.; Penna-Martinez, M.; Badenhoop, K. Focus on Vitamin D and the Adrenal Gland. Horm. Metab. Res. 2015, 47, 239–246. [Google Scholar] [CrossRef]
  16. Davidson, Z.E.; Walker, K.Z.; Truby, H. Clinical review: Do Glucocorticosteroids alter Vitamin D status? A systematic review with meta-analyses of observational studies. J. Clin. Endocrinol. Metab. 2012, 97, 738–744. [Google Scholar] [CrossRef]
  17. Hidalgo, A.A.; Trump, D.L.; Johnson, C.S. Glucocorticoid regulation of the vitamin D receptor. J. Steroid Biochem. Mol. Biol. 2010, 121, 372–375. [Google Scholar] [CrossRef]
  18. Hidalgo, A.A.; Deeb, K.K.; Pike, J.W.; Johnson, C.S.; Trump, D.L. Dexamethasone enhances 1α,25-Dihydroxyvitamin D3 effects by increasing Vitamin D receptor transcription. J. Biol. Chem. 2011, 286, 36228–36237. [Google Scholar] [CrossRef]
  19. Favus, M.J.; Kimberg, D.V.; Millar, G.N.; Gershon, E. Effects of cortisone administration on the metabolism and localization of 25-Hydroxycholecalciferol in the rat. J. Clin. Investig. 1973, 52, 1328–1335. [Google Scholar] [CrossRef]
  20. Kugai, N.; Koide, Y.; Yamashita, K.; Shimauchi, T.; Nagata, N.; Takatani, O. Impaired mineral metabolism in Cushing’s syndrome: Parathyroid function, vitamin D metabolites and osteopenia. Endocrinol. Jpn. 1986, 33, 345–352. [Google Scholar] [CrossRef]
  21. Aloia, J.F.; Roginsky, M.; Ellis, K.; Shukla, K.; Cohn, S. Skeletal metabolism and body composition in Cushing’s Syndrome. J. Clin. Endocrinol. Metab. 1974, 39, 981–985. [Google Scholar] [CrossRef]
  22. Findling, J.W.; Adams, N.D.; Lemann, J., Jr.; Gray, R.W.; Thomas, C.J.; Tyrrell, J.B. Vitamin D metabolites and parathyroid hormone in Cushing’s Syndrome: Relationship to calcium and phosphorus homeostasis. J. Clin. Endocrinol. Metab. 1982, 54, 1039–1044. [Google Scholar] [CrossRef] [PubMed]
  23. Seeman, E.; Kumar, R.; Hunder, G.G.; Scott, M.; Heath, H., 3rd; Riggs, B.L. Production, degradation, and circulating levels of 1,25-dihydroxyvitamin D in health and in chronic glucocorticoid excess. J. Clin. Investig. 1980, 66, 664–669. [Google Scholar] [CrossRef] [PubMed]
  24. Klein, R.G.; Arnaud, S.B.; Gallagher, J.C.; DeLuca, H.F.; Riggs, B.L. Intestinal calcium absorption in exogenous Hypercortisonism. J. Clin. Investig. 1977, 60, 253–259. [Google Scholar] [CrossRef] [PubMed]
  25. Chaiamnuay, S.; Chailurkit, L.-O.; Narongroeknawin, P.; Asavatanabodee, P.; Laohajaroensombat, S.; Chaiamnuay, P. Current daily glucocorticoid use and serum creatinine levels are associated with lower 25(OH) Vitamin D levels in Thai patients with systemic lupus erythematosus. JCR J. Clin. Rheumatol. 2013, 19, 121–125. [Google Scholar] [CrossRef]
  26. Slovik, D.M.; Neer, R.M.; Ohman, J.L.; Lowell, F.C.; Clark, M.B.; Segre, G.V.; Potts, J.T., Jr. Parathyroid hormone and 25-hydroxyvitamin D levels in glucocorticoid-treated patients. Clin. Endocrinol. 1980, 12, 243–248. [Google Scholar] [CrossRef]
  27. Lindgren, J.U.; Merchant, C.R.; DeLuca, H.F. Effect of 1,25-dihydroxyvitamin D3 on osteopenia induced by prednisolone in adult rats. Calcif. Tissue Res. 1982, 34, 253–257. [Google Scholar] [CrossRef]
  28. Corbee, R.; Tryfonidou, M.; Meij, B.; Kooistra, H.; Hazewinkel, H. Vitamin D status before and after hypophysectomy in dogs with pituitary-dependent hypercortisolism. Domest. Anim. Endocrinol. 2012, 42, 43–49. [Google Scholar] [CrossRef]
  29. Park, J.E.; Pichiah, P.T.; Cha, Y.-S. Vitamin D and metabolic diseases: Growing roles of Vitamin D. J. Obes. Metab. Syndr. 2018, 27, 223–232. [Google Scholar] [CrossRef]
  30. Medrano, M.; Carrillo-Cruz, E.; Montero, I.; Perez-Simon, J.A. Vitamin D: Effect on Haematopoiesis and immune system and clinical applications. Int. J. Mol. Sci. 2018, 19, 2663. [Google Scholar] [CrossRef]
  31. Fleseriu, M.; Auchus, R.; Bancos, I.; Ben-Shlomo, A.; Bertherat, J.; Biermasz, N.R.; Boguszewski, C.L.; Bronstein, M.D.; Buchfelder, M.; Carmichael, J.D.; et al. Consensus on diagnosis and management of Cushing’s disease: A guideline update. Lancet Diabetes Endocrinol. 2021, 9, 847–875. [Google Scholar] [CrossRef]
  32. Nieman, L.K.; Biller, B.M.K.; Findling, J.W.; Newell-Price, J.; Savage, M.O.; Stewart, P.M.; Montori, V. The diagnosis of Cushing’s Syndrome: An endocrine society clinical practice guideline. J. Clin. Endocrinol. Metab. 2008, 93, 1526–1540. [Google Scholar] [CrossRef] [PubMed]
  33. Expert Panel on Detection, Evaluation, Treatment of High Blood Cholesterol in Adults. Executive summary of the third report of the National Cholesterol Education Program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (Adult Treatment Panel III). JAMA 2001, 285, 2486–2497. [Google Scholar] [CrossRef] [PubMed]
  34. American Diabetes Association. Classification and diagnosis of diabetes: Standards of medical care in diabetes—2021. Diabetes Care 2021, 44 (Suppl. S1), S15–S33. [Google Scholar] [CrossRef] [PubMed]
  35. Guarnotta, V.; Amato, M.C.; Pivonello, R.; Arnaldi, G.; Ciresi, A.; Trementino, L.; Citarrella, R.; Iacuaniello, D.; Michetti, G.; Simeoli, C.; et al. The degree of urinary hypercortisolism is not correlated with the severity of cushing’s syndrome. Endocrine 2016, 55, 564–572. [Google Scholar] [CrossRef] [PubMed]
  36. Holick, M.F.; Binkley, N.C.; Bischoff-Ferrari, H.A.; Gordon, C.M.; Hanley, D.A.; Heaney, R.P.; Murad, M.H.; Weaver, C.M. Endocrine Society. Evaluation, treatment, and prevention of Vitamin D deficiency: An endocrine society clinical practice guideline. J. Clin. Endocrinol. Metab. 2011, 96, 1911–1930. [Google Scholar] [CrossRef]
  37. Fiebrich, H.-B.; Berg, G.V.D.; Kema, I.P.; Links, T.P.; Kleibeuker, J.H.; Van Beek, A.P.; Walenkamp, A.M.E.; Sluiter, W.J.; De Vries, E.G.E. Deficiencies in fat-soluble vitamins in long-term users of somatostatin analogue. Aliment. Pharmacol. Ther. 2010, 32, 1398–1404. [Google Scholar] [CrossRef]
  38. Matthews, D.R.; Hosker, J.P.; Rudenski, A.S.; Naylor, B.A.; Treacher, D.F.; Turner, R.C. Homeostasis model assessment: Insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 1985, 28, 412–419. [Google Scholar] [CrossRef]
  39. Matsuda, M.; DeFronzo, R.A. Insulin sensitivity indices obtained from oral glucose tolerance testing: Comparison with the euglycemic insulin clamp. Diabetes Care 1999, 22, 1462–1470. [Google Scholar] [CrossRef]
  40. Utzschneider, K.M.; Prigeon, R.L.; Faulenbach, M.V.; Tong, J.; Carr, D.B.; Boyko, E.J.; Leonetti, D.L.; McNeely, M.J.; Fujimoto, W.Y.; Kahn, S.E. Oral disposition index predicts the development of future diabetes above and beyond fasting and 2-h glucose levels. Diabetes Care 2009, 32, 335–341. [Google Scholar] [CrossRef]
  41. Glendenning, P.; Zhu, K.; Inderjeeth, C.; Howat, P.; Lewis, J.R.; Prince, R.L. Effects of three-monthly oral 150,000 IU cholecalciferol supplementation on falls, mobility, and muscle strength in older postmenopausal women: A randomized controlled trial. J. Bone Miner. Res. 2011, 27, 170–176. [Google Scholar] [CrossRef]
  42. Kearns, M.D.; Alvarez, J.A.; Tangpricha, V. Large, single-dose, oral Vitamin D supplementation in adult populations: A systematic review. Endocr. Pract. 2014, 20, 341–351. [Google Scholar] [CrossRef] [PubMed]
  43. Lundqvist, J.; Norlin, M.; Wikvall, K. 1α,25-Dihydroxyvitamin D3 affects hormone production and expression of steroidogenic enzymes in human adrenocortical NCI-H295R cells. Biochim. Biophys. Acta 2010, 1801, 1056–1062. [Google Scholar] [CrossRef] [PubMed]
  44. Lemieux, P.; Weisnagel, S.J.; Caron, A.Z.; Julien, A.-S.; Morisset, A.-S.; Carreau, A.-M.; Poirier, J.; Tchernof, A.; Robitaille, J.; Bergeron, J.; et al. Effects of 6-month vitamin D supplementation on insulin sensitivity and secretion: A randomised, placebo-controlled trial. Eur. J. Endocrinol. 2019, 181, 287–299. [Google Scholar] [CrossRef] [PubMed]
  45. Li, Y.; Tong, C.H.; Rowland, C.M.; Radcliff, J.; Bare, L.A.; McPhaul, M.J.; Devlin, J.J. Association of changes in lipid levels with changes in vitamin D levels in a real-world setting. Sci. Rep. 2021, 11, 21536. [Google Scholar] [CrossRef] [PubMed]
  46. Dibaba, D.T. Effect of vitamin D supplementation on serum lipid profiles: A systematic review and meta-analysis. Nutr. Rev. 2019, 77, 890–902. [Google Scholar] [CrossRef] [PubMed]
  47. Povaliaeva, A.; Bogdanov, V.; Pigarova, E.; Zhukov, A.; Dzeranova, L.; Belaya, Z.; Rozhinskaya, L.; Mel’Nichenko, G.; Mokrysheva, N. Assessment of Vitamin D metabolism in patients with Cushing’s disease in response to 150,000 IU cholecalciferol treatment. Nutrients 2021, 13, 4329. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Assessment of Vitamin D Metabolism in Patients with Cushing’s Disease

Endocrinology Research Centre, 117292 Moscow, Russia
*
Author to whom correspondence should be addressed.
Academic Editor: Spyridon N. Karras
Nutrients 202113(12), 4329; https://doi.org/10.3390/nu13124329
Received: 12 November 2021 / Revised: 26 November 2021 / Accepted: 27 November 2021 / Published: 30 November 2021

Abstract

In this study we aimed to assess vitamin D metabolism in patients with Cushing’s disease (CD) compared to healthy individuals in the setting of bolus cholecalciferol treatment. The study group included 30 adults with active CD and the control group included 30 apparently healthy adults with similar age, sex and BMI. All participants received a single dose (150,000 IU) of cholecalciferol aqueous solution orally. Laboratory assessments including serum vitamin D metabolites (25(OH)D3, 25(OH)D2, 1,25(OH)2D3, 3-epi-25(OH)D3 and 24,25(OH)2D3), free 25(OH)D, vitamin D-binding protein (DBP) and parathyroid hormone (PTH) as well as serum and urine biochemical parameters were performed before the intake and on Days 1, 3 and 7 after the administration. All data were analyzed with non-parametric statistics. Patients with CD had similar to healthy controls 25(OH)D3 levels (p > 0.05) and higher 25(OH)D3/24,25(OH)2D3 ratios (p < 0.05) throughout the study. They also had lower baseline free 25(OH)D levels (p < 0.05) despite similar DBP levels (p > 0.05) and lower albumin levels (p < 0.05); 24-h urinary free cortisol showed significant correlation with baseline 25(OH)D3/24,25(OH)2D3 ratio (r = 0.36, p < 0.05). The increase in 25(OH)D3 after cholecalciferol intake was similar in obese and non-obese states and lacked correlation with BMI (p > 0.05) among patients with CD, as opposed to the control group. Overall, patients with CD have a consistently higher 25(OH)D3/24,25(OH)2D3 ratio, which is indicative of a decrease in 24-hydroxylase activity. This altered activity of the principal vitamin D catabolism might influence the effectiveness of cholecalciferol treatment. The observed difference in baseline free 25(OH)D levels is not entirely clear and requires further study.

1. Introduction

Cushing’s disease (CD) is one of the disorders associated with endogenous hypercortisolism and is caused by adrenocorticotropic hormone (ACTH) hyperproduction originating from pituitary adenoma [1]. Skeletal fragility is a frequent complication of endogenous hypercortisolism, and fragility fractures may be the presenting clinical feature of disease. The prevalence of osteoporosis in endogenous hypercortisolism as assessed by dual-energy X-ray absorptiometry (DXA) or incidence of fragility fractures has been reported to be up to 50%. Osteoporosis in CD patients has a complex multifactorial pathogenesis, characterized by a low bone turnover and severe suppression of bone formation [2]. Exogenous glucocorticoids are used in the treatment of a wide range of diseases and it is estimated that 1–2% of the population is receiving long-term glucocorticoid therapy. As a consequence, glucocorticoid-induced osteoporosis is the most common secondary cause of osteoporosis [3].
Native vitamin D (in particular D3, or cholecalciferol) and its active metabolites (such as alfacalcidol) are universally considered as the essential components of the osteoporosis management [4,5]. The search for the optimal treatment of bone complications during chronic exposure to glucocorticoid excess provoked the investigation of vitamin D metabolism in this state. Early studies on this topic were focused predominantly on the general vitamin D status (assessed as 25(OH)D level) and on the levels of the active vitamin D metabolite (1,25(OH)2D). These studies showed inconsistent results, reporting that the chronic excess of glucocorticoids decreased [6,7,8,9], increased [10,11,12] or did not change [13,14,15] the levels of 25(OH)D or 1,25(OH)2D. A likely reason for such inconsistency might have been the high heterogeneity of the studied groups. Some of these studies were performed in humans [6,7,9,10,11,12,13,15] and some in animal models [8,14], and only several of them included subjects with specifically endogenous hypercortisolism [10,12,14,15]. Only two studies assessed both the levels of the active (1,25(OH)2D) and the inactive (24,25(OH)2D) vitamin D metabolites in endogenous hypercortisolism. One of them lacked control group and reported low-normal 24,25(OH)2D levels in patients with Cushing’s syndrome [10]. The second study by Corbee et al. reported similar circulating concentrations of 25(OH)D, 1,25(OH)2D and 24,25(OH)2D in studied groups of dogs regardless of either the presence of CD or hypophysectomy status [14].
Several experimental studies were performed to evaluate the impact of glucocorticoid excess on the enzymes involved in vitamin D metabolism. In mouse kidney glucocorticoid treatment increased 24-hydroxylase expression [16] and 24-hydroxylase activity [17]. An increased expression of 24-hydroxylase was also shown in rat osteoblastic and pig renal cell cultures treated with 1,25(OH)2D [18]. Dhawan and Christakos showed that 1,25(OH)2D-induced transcription of 24-hydroxylase was glucocorticoid receptor-dependent [19]. However, some works showed conflicting results. In particular, the steroid and xenobiotic receptor (SXR) which is activated by glucocorticoids [20], repressed 24-hydroxylase expression in human liver and intestine in work by Zhou et al. [21]. Lower 24-hydroxylase expression was observed in the brain and myocardium of glucocorticoid-treated rats [22] as well as in human osteosarcoma cells and human osteoblasts [23].
Nevertheless, based on experimental data, it has been suggested that the acceleration of 25(OH)D catabolism in the presence of glucocorticoid excess may predispose to vitamin D deficiency. Yet, relatively recent meta-analysis of the studies assessing 25(OH)D levels in chronic glucocorticoid users showed that serum 25(OH)D levels in these patients were suboptimal and lower than in healthy controls, but similar to steroid-naive disease controls [24].
Glucocorticoids also affect calcium and phosphorus homeostasis. In particular, they were shown to reduce gastrointestinal absorption by antagonizing vitamin D action (reducing the expression of genes for proteins involved in calcium transport—epithelial Ca channel TRPV6 and calcium-binding protein calbindin-D9K) [25]. Glucocorticoids increased fractional calcium excretion due to mineralocorticoid receptor-mediated action on epithelial sodium channels [26]. Hypercalciuria is highly prevalent in people with CD [27]. These effects might result in a negative calcium balance, although plasma ionized calcium was normal in people and dogs with hypercortisolism compared to control subjects [12,28]. Glucocorticoids also reduced tubular phosphate reabsorption by inhibiting tubular expression of the sodium gradient-dependent phosphate transporter, and induced phosphaturia [29], which was accompanied by phosphate lowering in humans [12].
Overall, current data on vitamin D status in hypercortisolism are conflicting and need clarification. In particular, clinical data on the state of vitamin D metabolism in the state of glucocorticoids excess are quite scarce. Studies were very heterogeneous in design, some lacked a control group, and the absolute majority of the studies were performed before the introduction of vitamin D measurement standardization [30]. Nevertheless, determining the optimal vitamin D treatment regimen in these high-risk patients is fairly relevant.
The aim of this study was to assess vitamin D metabolism in patients with CD compared to healthy individuals particularly in the setting of cholecalciferol treatment.

2. Materials and Methods

2.1. Study Population and Design

The study group included 30 adult patients with CD admitted for inpatient treatment at a tertiary pituitary center. Diagnosis of CD was established in accordance with the federal guidelines [31]. All patients were confirmed to be positive for endogenous hypercortisolism in at least two of the following tests: 24-h urine free cortisol (UFC) greater than the normal range for the assay and/or serum cortisol > 50 nmol/L after the 1-mg overnight dexamethasone suppression test and/or late-night salivary cortisol greater than 9.4 nmol/L). All patients also had morning ACTH ≥ 10 pg/mL and pituitary adenoma ≥ 6 mm identified by magnetic resonance imaging (MRI) or a positive for CD bilateral inferior petrosal sinus sampling (BIPSS). MRI was performed using a GE Optima MR450w 1.5T with Gadolinium (Boston, MA, USA). BIPSS was performed according to the standard procedure described elsewhere [32,33].
The control group included 30 apparently healthy adult individuals recruited from the staff and the faculty of the facility.
Inclusion criteria were age from 18 to 60 for both groups and the presence of the disease activity for the study group (defined as the presence of endogenous hypercortisolism at the time of participation in the study). Exclusion criteria for both groups were: vitamin D supplementation for 3 months prior to the study; severe obesity (body mass index (BMI) ≥ 35 kg/m2); pregnancy; the presence of granulomatous disease, malabsorption syndrome, liver failure; decreased GFR (less than 60 mL/min per 1.73 m2); severe hypercalcemia (total serum calcium > 3.0 mmol/L); allergic reactions to vitamin D medications; 25(OH)D level more than 60 ng/mL (determined by immunochemiluminescence analysis). All patients were recruited in the period from October 2019 to April 2021. The study protocol (ClinicalTrials.gov Identifier: NCT04844164) was approved by the Ethics Committee of Endocrinology Research Centre, Moscow, Russia on 10 April 2019 (abstract of record No. 6), all patients signed informed consent to participate in the study.
All participants received standard therapeutic dose (150,000 IU) of an aqueous solution of cholecalciferol (Aquadetrim®, Medana Pharma S.A., Sieradz, Poland) orally as a single dose [34]. Blood and urine samples were obtained before the intake as well as on days 1, 3 and 7 after administration; time points of sample collection were determined based on the authors’ previous work evaluating changes in 25(OH)D levels after a therapeutic dose of cholecalciferol [35]. The assessment included serum biochemical parameters (total calcium, albumin, phosphorus, creatinine, magnesium), parathyroid hormone (PTH), vitamin D-binding protein (DBP), vitamin D metabolites (25(OH)D3, 25(OH)D2, 1,25(OH)2D3, 3-epi-25(OH)D3 and 24,25(OH)2D3), free 25(OH)D and urine biochemical parameters (calcium- and phosphorus-creatinine ratios in spot urine).

2.2. Socio–Demographic and Anthropometric Data Collection

At the baseline visit, patients underwent a questionnaire aimed to assess their lifestyle: the presence of unhealthy habits, physical activity level, balanced diet (consumption of dairy products, meat, coffee, soft drinks), exposure to ultraviolet (UV) radiation (solarium and sunscreen usage, traveling south and the number of daytime walks in the sunny weather in the 3 months preceding study participation). Smoking status was classified as current smoker, former smoker and non-smoker; current and former smokers were collectively referred to as total smokers. A unit of alcohol was defined as a glass of wine, a bottle of beer or a shot of spirits, approximating 10–12 g ethanol. Serving of dairy products was defined as 100 g of cottage cheese, 200 mL of milk, 125 g of yogurt or 30 g of cheese. Patients’ weight was measured in light indoor clothing with a medical scale to the nearest 100 g, and their height with a wall-mounted stadiometer to the nearest centimeter. BMI was calculated as weight in kilograms divided by height in meters squared.

2.3. Laboratory Measurements

Morning ACTH (reference range 7–66 pg/mL), serum cortisol after a low-dose dexamethasone suppression test (cutoff value for suppression, 50 nmol/L [36]), late-night salivary cortisol (reference range 0.5–9.4 nmol/L [37]) were assayed by electrochemiluminescence assay using a Cobas 6000 Module e601 (Roche, Rotkreuz, Switzerland). The 24-h UFC (reference range 60–413 nmol/24 h) was measured by an immunochemiluminescence assay (extraction with diethyl ether) on a Vitros ECiQ (Ortho Clinical Diagnostics, Raritan, NJ, USA).
Total 25(OH)D levels (25(OH)D2 + 25(OH)D3; reference range 30–100 ng/mL) at the baseline visit were determined by the immunochemiluminescence analysis (Liaison, DiaSorin, Saluggia, Italy). PTH levels were evaluated by the electrochemiluminescence immunoassay (ELECSYS, Roche, Basel, Switzerland; reference range for this and subsequent laboratory parameters are given in the Results section for easier reading). Biochemical parameters of blood serum and urine were assessed by the ARCHITECT c8000 analyzer (Abbott, Chicago, IL, USA) using reagents from the same manufacturer according to the standard methods. Serum DBP and free 25(OH)D levels were measured by enzyme-linked immunosorbent assay (ELISA) using commercial kits. The assay used for free 25(OH)D levels assessment (DIAsource, ImmunoAssays S.A., Ottignies-Louvain-la-Neuve, Belgium) has <6.2% intra- and inter-assay coefficient of variation (CV) at levels 5.8–9.6 pg/mL. The assay used for DBP levels assessment (Assaypro, St Charles, MO, USA) has 6.2% average intra-assay CV and 9.9% average inter-assay CV.
The levels of vitamin D metabolites (25(OH)D3, 25(OH)D2, 1,25(OH)2D3, 3-epi-25(OH)D3 and 24,25(OH)2D3) in serum were determined by ultra-high performance liquid chromatography in combination with tandem mass spectrometry (UPLC-MS/MS) using an in-house developed method, described earlier [38]. With this technique, the laboratory participates in DEQAS quality assurance program (lab code 2388) and the results fall within the target range for the analysis of 25(OH)D and 1,25(OH)2D metabolites in human serum (Supporting Information, Figures S1 and S2). All UPLC-MS/MS measurements were made after the first successful completion (5/5 samples within the target range) of the DEQAS distributions for both analytes simultaneously. Each batch contained control samples (analytes in blank serum) with both high and low analyte concentrations. The samples were barcoded and randomized prior to the measurements to eliminate analyst-related errors.
Serum samples (3 aliquots) collected at each visit were either transferred directly to the laboratory for biochemical analyzes, total 25(OH)D and PTH measurement (1 aliquot) or were stored at −80 °C avoiding repeated freeze-thaw cycles for measurement of DBP, free 25(OH)D and vitamin D metabolites at a later date (2 aliquots).
Albumin-adjusted serum calcium levels were calculated using the formula [39]: total plasma calcium (mmol/L) = measured total plasma calcium (mmol/L) + 0.02 × (40 − measured plasma albumin (g/L)).
Baseline free 25(OH)D levels were also calculated using the formula introduced by Bikle et al. [40,41]. The affinity constant for 25(OH)D and albumin binding (Kalb) used for the calculation was equal 6 × 105 M−1, and affinity constant for 25(OH)D and DBP binding (KDBP) was equal 7 × 108 M−1.

Free 25(OH)D=total 25(OH)D1+Kalbalbumin+KDBPDBP

2.4. Statistical Analysis

Statistical analysis was performed using Statistica version 13.0 (StatSoft, Tulsa, OK, USA). All data were analyzed with non-parametric statistics and expressed as median [interquartile range] unless otherwise specified. Mann-Whitney U-test and Fisher’s exact two-tailed test were used for comparisons between two groups. Friedman ANOVA was performed to evaluate changes in indices throughout the study and pairwise comparisons using Wilcoxon test with adjustment for multiple comparisons (Bonferroni) were also made if the Friedman ANOVA was significant. Spearman rank correlation method was used to obtain correlation coefficients among indices. A p-value of less than 0.05 was considered statistically significant. When adjusting for multiple comparisons, a p-value greater than the significance threshold, but less than 0.05 was considered as a trend towards statistical significance.

3. Results

The groups were similar in terms of age, sex and BMI (p > 0.05). Both groups consisted predominantly of young and middle-aged women and the majority of patients were overweight or moderately obese (Table 1). Patients from the study group presented with lower screening levels of total 25(OH)D (p < 0.05).
Table 1. General characteristics of the patients at the baseline visits. For detailed description of the data format please refer to the Section 2.
The features of the underlying disease course in the study group are listed in Table 2. 15 patients (50%) had diabetes mellitus with an almost compensated state at the time of participation in the study, and 7 patients (23%) reported a history of low-energy fractures.
Table 2. Characteristics of the patients with Cushing’s disease (CD) in terms of the underlying disease.
The groups did not differ significantly in the reported smoking status, the level of daily physical activity, dietary habits and UV exposure (p > 0.05) and although there was a slight difference in alcohol consumption (p < 0.05), the absolute values were minor in both groups (Table 3).
Table 3. Questionnaire results.

3.1. Baseline Laboratory Evaluation

Detailed results of laboratory studies are presented in Table 4 and Table 5.
Table 4. Changes in the levels of the biochemical parameters and parathyroid hormone (PTH) during the study.
Table 5. Changes in the levels of free 25(OH)D, vitamin D-binding protein (DBP) and vitamin D metabolites during the study.
Patients with CD had several alterations in biochemical parameters, in particular, lower baseline serum creatinine and albumin levels, while magnesium levels were higher than in the control group (p < 0.05). They also had higher levels of urine phosphorus-creatinine ratio (p < 0.05). The rest of the studied biochemical parameters did not show significant difference between the groups (p > 0.05). 3 patients (10%) from the study group and 5 patients (17%) from the control group had secondary hyperparathyroidism, one patient with CD (3%) was diagnosed with mild primary hyperparathyroidism.
As for the assessment of vitamin D metabolism, unexpectedly the levels of 25(OH)D3 occurred to be equal in the groups (p > 0.05), with only two patients (7%) from the study group and one patient (3%) from the control group having sufficient vitamin D levels, according to the Endocrine Society and the Russian Association of Endocrinologists guidelines (≥30 ng/mL [34,42]). The levels of the active vitamin D metabolite—1,25(OH)2D3—were equal between the groups as well (p > 0.05), whereas the levels of 3-epi-25(OH)D3 and 24,25(OH)2D3 were lower in CD patients. Further calculation of 25(OH)D3/24,25(OH)2D3 and 25(OH)D3/1,25(OH)2D3 ratios corresponded to the observed levels of metabolites: 25(OH)D3/24,25(OH)2D3 ratio was higher in the study group (p < 0.05) assuming lower 24-hydroxylase activity and 25(OH)D3/1,25(OH)2D3 ratio was equal between the groups (p > 0.05).
Levels of free 25(OH)D were lower in CD patients (p < 0.05) and the levels of DBP did not differ between the groups (p > 0.05). Although calculated free 25(OH)D showed prominent positive correlation with the measured free 25(OH)D in both groups (r = 0.63 in the study group, r = 0.87 in the control group, p < 0.05), the association appeared to be weaker in the study group. In the control group, DBP levels correlated with both measured and calculated 25(OH)D levels (r = −0.48, p < 0.05 and r = −0.69, p < 0.05 respectively), while in patients with CD there was no association with measured free 25(OH)D levels (r = 0.04, p > 0.05 and r = −0.50, p < 0.05 respectively).
Correlation with 24-h UFC in CD patients was observed for serum albumin level (r = −0.37, p < 0.05) and urine calcium-creatinine ratio (r = 0.51, p < 0.05) among assessed biochemical parameters, and only with 25(OH)D3/24,25(OH)2D3 ratio among the parameters of vitamin D metabolism (r = 0.36, p < 0.05).

3.2. Laboratory Evaluation after the Intake of Cholecalciferol

All patients from the study group and 28 patients (93%) from the control group completed the study.
The observed baseline differences in biochemical parameters mostly preserved during the follow-up. In the study group there was an increase in serum phosphorus levels by Day 1 (p = 0.006) and a tendency to an increase in the urine phosphorus-creatinine ratio by Day 7 (p = 0.02). Patients from the control group showed a clinically insignificant increase in serum creatinine levels by Day 1 (p = 0.002) and a non-significant trend towards an increase in serum total and albumin-adjusted calcium (p = 0.01 for both measurements). No change in PTH levels was observed in patients with CD during the follow-up (p > 0.05), while in the control group there was a tendency for PTH to decrease by Day 3 (p = 0.02). There were no new cases of hypercalcemia in both groups during the follow-up. One patient from the study group and one patient from the control group had persistently increased urine calcium-creatinine ratio throughout the study. Four patients from the study group (13%) and none from the control group developed hypercalciuria during the follow-up, however these patients had no clinical manifestations during the observation period.
By Day 7, 25 patients (83%) from the study group and 22 patients (79%) reached sufficient 25(OH)D3 levels (≥30 ng/mL). Levels of 25(OH)D3 continued to increase by Day 3 in both groups (p < 0.001), after which tended to decrease in the study group (p = 0.01) and remained stable in the control group (p = 0.65). The increase in 25(OH)D3 after cholecalciferol intake was equal between the groups (18.5 [15.9; 22.5] ng/mL in the study group vs. 16.6 [13.1; 19.8] ng/mL in the control group, p > 0.05). In the presence of obesity, Δ25(OH)D3 was higher in the CD patients than in the control group (18.3 [14.2; 23.0] vs. 12.1 [10.0; 13.1] ng/mL, p < 0.05), while in non-obese patients no difference was observed (p > 0.05).
Obese and non-obese patients with CD had equal Δ25(OH)D3 (18.3 [14.2; 23.0] vs. 19.6 [16.0; 21.5] ng/mL, p > 0.05), while in the control group it was significantly lower in obese patients (12.1 [10.0; 13.1] vs. 18.3 [15.3; 21.4] ng/mL, p < 0.05). BMI showed significant correlation with Δ25(OH)D3 only in the control group (r = −0.47, p < 0.05), while in CD patients there was no such association (r = −0.06, p > 0.05) (Figure 1).
Figure 1. Relationship between Δ25(OH)D3 and BMI in groups.
1,25(OH)2D3 levels increased in CD patients by Day 1 and were stable during the follow-up in the control group. The rest of the studied parameters of vitamin D metabolism changed in a similar way between groups: 3-epi-25(OH)D3 levels increased until the Day 3, after which they decreased by the Day 7; 24,25(OH)2D3 levels showed more graduate elevation throughout the follow-up. In both groups 25(OH)D3/24,25(OH)2D3 ratios increased by Day 1, after which they decreased by Day 7, and 25(OH)D3/1,25(OH)2D3 ratios increased by Day 1, after which they remained stable. DBP levels didn’t change and free 25(OH)D levels showed an increase in both groups during the follow-up. The levels of 25(OH)D2 did not exceed 0.5 ng/mL in all examined individuals throughout the study. Among assessed parameters of vitamin D metabolism, higher 25(OH)D3/24,25(OH)2D3 ratios in the study group was the only difference between the groups which remained significant throughout the observation period (p < 0.05) (Figure 2).
Figure 2. Dynamic evaluation of 25(OH)D3/24,25(OH)2D3 ratios in groups.

4. Discussion

The main goal of our study was to evaluate the 25(OH)D3 levels and its response to the therapeutic dose of cholecalciferol in patients with CD as compared to healthy individuals. We observed no difference in baseline 25(OH)D3 assessed by UPLC-MS/MS between groups. Similar to our data were obtained in most studies conducted specifically in the state of endogenous hypercortisolism in humans [12,15] and dogs [14]. The study by Kugai et al. lacked control group and reported plasma levels of 25(OH)D corresponding to the vitamin D deficiency in most of the examined patients [10], while in our study only 2/3 of the patients with CD had 25(OH)D levels below 20 ng/mL. As for exogenous hypercortisolism, the meta-analysis aimed to explore serum 25(OH)D levels in glucocorticoid users showed lower levels than in healthy controls, but similar to steroid-naive disease controls, thus causing concern regarding the influence of the disease status on 25(OH)D levels [24]. Somewhat surprisingly, we obtained significantly discordant results in the study group when screening total 25(OH)D by ELISA and when measuring baseline 25(OH)D3 by UPLC-MS/MS, since the initial difference between the groups revealed by ELISA data with lower total 25(OH)D levels in the study group was not replicated by UPLC-MS/MS. It should be noted that our ELISA method did not participate in an external quality control program at the time of the study unlike UPLC-MS/MS; furthermore, a lower analytical performance was previously described for this technique with tendency for low specificity and lower measurement results [45].
When assessing other parameters of vitamin D metabolism, the most significant finding was the higher 25(OH)D3/24,25(OH)2D3 ratio in CD patients, both initially and during the observation after the intake of the cholecalciferol loading dose, indicating consistently reduced activity of 24-hydroxylase, the main enzyme of vitamin D catabolism. Earlier clinical and experimental studies also suggested altered activity of enzymes of vitamin D metabolism in hypercortisolism. However, these studies were heterogeneous and aimed predominantly at studying the activity of 1α-hydroxylase [7,8,10,11,12,14], which was not altered in patients with CD as compared to healthy individuals in our study. In the setting of the short-term glucocorticoid administration, Lindgren et al. showed transient increase in 24,25(OH)2D3 levels in rats [8], while in the study of Hahn et al. there was no change in 24,25(OH)2D3 levels [11]. Dogs with CD had similar 24,25(OH)2D3 levels before and after hypophysectomy as well as compared to control dogs [14]. The only study of considerably similar design by Kugai et al. reported low-normal 24,25(OH)2D3 in patients with Cushing’s syndrome [10], which is consistent with our result, as well as some experimental works indicative of suppression on CYP24A1 expression by glucocorticoids in human osteoblasts [23], liver and intestine [21] and in rat brain and myocardium [22]. However, in the present work, the activity of 24-hydroxylase in patients with hypercortisolism was for the first time evaluated by calculating the 25(OH)D3/24,25(OH)2D3 ratio, which has recently emerged as a new tool for vitamin D status assessment [46,47]. Given the correlation of this parameter with laboratory marker of the underlying disease activity (24-h UFC), a direct effect of cortisol overproduction on 24-hydroxylase activity might be assumed. Interestingly, it seems that the decreased activity of 24-hydroxylase observed in CD influenced the effectiveness of cholecalciferol treatment, decreasing the negative effect of obesity, as patients with CD had similar increase in 25(OH)D3 in obese and non-obese state and lacked correlation between Δ25(OH)D3 and BMI, as opposed to the control group. Moreover, the increase in 25(OH)D3 in obese patients from the control group was lower not only than in non-obese controls, but also than in obese patients with CD.
Another intriguing finding was lower levels of free 25(OH)D observed in patients with CD despite similar DBP levels and lower albumin levels, which, on the contrary, allows one to expect higher values of free 25(OH)D. Considering the weaker correlation between the measured and calculated free 25(OH)D in patients with CD, as well as the lack of correlation of the measured 25(OH)D with the main transport protein, an altered affinity of DBP might be suspected. One possible explanation is protein glycosylation as a consequence of diabetes mellitus, which was present in half of the patients [38,48,49]. After cholecalciferol intake, which was accompanied by an increase in free 25(OH)D, the differences between the groups were leveled; therefore, another suggested explanation might be competitive binding to the ligand. Since actin binds DBP with high affinity [50] and considering catabolic action of glucocorticoids on muscle tissue [51], actin is a presumable competing ligand candidate. Although this is mostly speculative, as far as the authors are aware, the present work was the first to assess free vitamin D in the glucocorticoid excess, so the described findings require verification of reproducibility and further evaluation.
The obtained discrepancies in the biochemical parameters characterizing calcium and phosphorus metabolism were generally consistent with the data of early studies discussed in the introduction [12,25,26,27,28,29], except for similar to controls serum phosphorus levels and lower prevalence of hypercalciuria. An interesting observation was the complete absence of the PTH decrease in patients with CD after receiving a loading dose of cholecalciferol. The mechanism of this phenomenon is not entirely clear, we tend to agree with the earlier hypothesis that this may be an adaptation to chronic urinary calcium loss [52].
Our research is distinguished by a number of important strengths: a prospective design, substantial sample of patients with CD, accounting for social and behavioral factors affecting vitamin levels D, comprehensive spectrum of vitamin D metabolism parameters investigated and participation in an external quality control program for vitamin D metabolites measurement.
Nevertheless, the study also had several limitations: the amount of dietary vitamin D and phosphorus, as well as possible differences in DBP affinity to vitamin D metabolites due to genetic isoforms of DBP [53] or other possible involved parameters (e.g., fibroblast growth factor-23) were not taken into account. A few patients from both groups received therapy with possible impact on vitamin D and calcium metabolism within 3 months preceding the participation in the study (spironolactone, diuretics, proton pump inhibitors, oral contraceptives, antifungal treatment, antidepressants, barbiturates, antiepileptic drugs). The groups had a trend for differences in sex and BMI (p = 0.07 for both parameters). Also, the study lacked a study group of patients with remission of CD to test the hypotheses put forward, however, this is a promising direction for further research.

5. Conclusions

We report that patients with endogenous ACTH-dependent hypercortisolism of pituitary origin have a consistently higher 25(OH)D3/24,25(OH)2D3 ratio than healthy controls, which is indicative of a decrease in 24-hydroxylase activity. This altered activity of the principal vitamin D catabolism might influence the effectiveness of cholecalciferol treatment. There is also a lack of clarity regarding the lower levels of free 25(OH)D observed in patients with CD, which require further study. To test the proposed hypotheses and to develop specialized clinical guidelines for these patients, longer-term randomized clinical trials are needed.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13124329/s1, Method validation against DEQAS, Figure S1: Comparison between DEQAS data for 25(OH)D scheme and our lab results, Figure S2: Comparison between DEQAS data for 1,25(OH)2D scheme and our lab results.

Author Contributions

Conceptualization, L.R., E.P., A.P. and A.Z.; methodology, V.B., Z.B., L.R. and G.M.; formal analysis, A.P.; investigation, A.P., V.B., E.P., L.D. and A.Z.; data curation, A.P. and V.B.; writing—original draft preparation, A.P.; writing—review and editing, V.B., E.P., A.Z., Z.B., L.R.; visualization, V.B.; supervision, L.D., L.R., G.M. and N.M.; project administration, L.R. and N.M.; funding acquisition, L.R. and N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant number 19-15-00243.

Institutional Review Board Statement

This study was performed in line with the principles of the Declaration of Helsinki. Approval was granted by the Ethics Committee of Endocrinology Research Centre, Moscow, Russia on 10 April 2019 (abstract of record No. 6).

Informed Consent Statement

Written informed consent was obtained from all individual participants included in the study.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We express our deep gratitude to our colleagues: Natalya M. Malysheva, Vitaliy A. Ioutsi, Larisa V. Nikankina for the help with the laboratory research.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Nishioka, H.; Yamada, S. Cushing’s disease. J. Clin. Med. 20198, 1951. [Google Scholar] [CrossRef]
  2. Mazziotti, G.; Frara, S.; Giustina, A. Pituitary Diseases and Bone. Endocr. Rev. 201839, 440–488. [Google Scholar] [CrossRef] [PubMed]
  3. Compston, J. Glucocorticoid-induced osteoporosis: An update. Endocrine 201861, 7–16. [Google Scholar] [CrossRef] [PubMed]
  4. Buckley, L.; Guyatt, G.; Fink, H.A.; Cannon, M.; Grossman, J.; Hansen, K.E.; Humphrey, M.B.; Lane, N.E.; Magrey, M.; Miller, M. 2017 American College of Rheumatology Guideline for the prevention and treatment of glucocorticoid-induced osteoporosis. Arthritis Care Res. 201769, 1095–1110. [Google Scholar] [CrossRef]
  5. Belaya, Z.; Rozhinskaya, L.y.; Grebennikova, T.; Kanis, J.; Pigarova, E.; Rodionova, S.; Toroptsova, N.; Nikitinskaya, O.; Skripnikova, I.; Drapkina, O.; et al. Summary of the Draft Federal Clinical Guidelines for Osteoporosis. Osteoporos. Bone Dis. 202023, 4–21. [Google Scholar] [CrossRef]
  6. Klein, R.G.; Arnaud, S.B.; Gallagher, J.C.; Deluca, H.F.; Riggs, B.L. Intestinal calcium absorption in exogenous hypercortisolism. J. Clin. Investig. 197760, 253–259. [Google Scholar] [CrossRef] [PubMed]
  7. Seeman, E.; Kumar, R.; Hunder, G.G.; Scott, M.; Iii, H.H.; Riggs, B.L. Production, degradation, and circulating levels of 1,25-dihydroxyvitamin D in health and in chronic glucocorticoid excess. J. Clin. Investig. 198066, 664–669. [Google Scholar] [CrossRef]
  8. Lindgren, J.U.; Merchant, C.R.; DeLuca, H.F. Effect of 1,25-dihydroxyvitamin D3 on osteopenia induced by prednisolone in adult rats. Calcif. Tissue Int. 198234, 253–257. [Google Scholar] [CrossRef]
  9. Chaiamnuay, S.; Chailurkit, L.O.; Narongroeknawin, P.; Asavatanabodee, P.; Laohajaroensombat, S.; Chaiamnuay, P. Current daily glucocorticoid use and serum creatinine levels are associated with lower 25(OH) vitamin D levels in thai patients with systemic lupus erythematosus. J. Clin. Rheumatol. 201319, 121–125. [Google Scholar] [CrossRef]
  10. Kugai, N.; Koide, Y.; Yamashita, K.; Shimauchi, T.; Nagata, N.; Takatani, O. Impaired mineral metabolism in Cushing’s syndrome: Parathyroid function, vitamin D metabolites and osteopenia. Endocrinol. Jpn. 198633, 345–352. [Google Scholar] [CrossRef]
  11. Hahn, T.J.; Halstead, L.R.; Baran, D.T. Effects of short term glucocorticoid administration on intestinal calcium absorption and circulating vitamin D metabolite concentrations in man. J. Clin. Endocrinol. Metab. 198152, 111–115. [Google Scholar] [CrossRef]
  12. Findling, J.W.; Adams, N.D.; Lemann, J.; Gray, R.W.; Thomas, C.J.; Tyrrell, J.B. Vitamin D metabolites and parathyroid hormone in Cushing’s syndrome: Relationship to calcium and phosphorus homeostasis. J. Clin. Endocrinol. Metab. 198254, 1039–1044. [Google Scholar] [CrossRef] [PubMed]
  13. Slovik, D.M.; Neer, R.M.; Ohman, J.L.; Lowell, F.C.; Clark, M.B.; Segre, G.V.; Potts, J.T., Jr. Parathyroid hormone and 25-hydroxyvitamin D levels in glucocorticoid-treated patients. Clin. Endocrinol. 198012, 243–248. [Google Scholar] [CrossRef]
  14. Corbee, R.J.; Tryfonidou, M.A.; Meij, B.P.; Kooistra, H.S.; Hazewinkel, H.A.W. Vitamin D status before and after hypophysectomy in dogs with pituitary-dependent hypercortisolism. Domest. Anim. Endocrinol. 201242, 43–49. [Google Scholar] [CrossRef] [PubMed]
  15. Aloia, J.F.; Roginsky, M.; Ellis, K.; Shukla, K.; Cohn, S. Skeletal metabolism and body composition in Cushing’s syndrome. J. Clin. Endocrinol. Metab. 197439, 981–985. [Google Scholar] [CrossRef]
  16. Van Cromphaut, S.J.; Stockmans, I.; Torrekens, S.; Van Herck, E.; Carmeliet, G.; Bouillon, R. Duodenal calcium absorption in dexamethasone-treated mice: Functional and molecular aspects. Arch. Biochem. Biophys. 2007460, 300–305. [Google Scholar] [CrossRef] [PubMed]
  17. Akeno, N.; Matsunuma, A.; Maeda, T.; Kawane, T.; Horiuchi, N. Regulation of vitamin D-1-hydroxylase and -24-hydroxylase expression by dexamethasone in mouse kidney. J. Endocrinol. 2000164, 339–348. [Google Scholar] [CrossRef] [PubMed]
  18. Kurahashi, I.; Matsunuma, A.; Kawane, T.; Abe, M.; Horiuchi, N. Dexamethasone enhances vitamin D-24-hydroxylase expression in osteoblastic (UMR-106) and renal (LLC-PK 1) cells treated with 1a, 25-dihydroxyvitamin D3. Endocrine 200217, 109–118. [Google Scholar] [CrossRef]
  19. Dhawan, P.; Christakos, S. Novel regulation of 25-hydroxyvitamin D3 24-hydroxylase (24(OH)ase) transcription by glucocorticoids: Cooperative effects of the glucocorticoid receptor, C/EBPb, and the vitamin D receptor in 24(OH)ase transcription. J. Cell. Biochem. 2010110, 1314–1323. [Google Scholar] [CrossRef] [PubMed]
  20. Luo, G.; Cunningham, M.; Kim, S.; Burn, T.; Lin, J.; Sinz, M.; Hamilton, G.; Rizzo, C.; Jolley, S.; Gilbert, D.; et al. CYP3A4 induction by drugs: Correlation between a pregnane X receptor reporter gene assay and CYP3A4 expression in human hepatocytes. Drug Metab. Dispos. 200230, 795–804. [Google Scholar] [CrossRef]
  21. Zhou, C.; Assem, M.; Tay, J.C.; Watkins, P.B.; Blumberg, B.; Schuetz, E.G.; Thummel, K.E. Steroid and xenobiotic receptor and vitamin D receptor crosstalk mediates CYP24 expression and drug-induced osteomalacia. J. Clin. Investig. 2006116, 1703–1712. [Google Scholar] [CrossRef] [PubMed]
  22. Jiang, P.; Xue, Y.; Li, H.; Liu, Y. Dysregulation of vitamin D metabolism in the brain and myocardium of rats following prolonged exposure to dexamethasone. Psychopharmacology 2014231, 3445–3451. [Google Scholar] [CrossRef] [PubMed]
  23. Zayny, A.; Almokhtar, M.; Wikvall, K.; Ljunggren, Ö.; Ubhayasekera, K.; Bergquist, J.; Kibar, P.; Norlin, M. Effects of glucocorticoids on vitamin D3-metabolizing 24-hydroxylase (CYP24A1) in Saos-2 cells and primary human osteoblasts. Mol. Cell. Endocrinol. 2019496, 110525. [Google Scholar] [CrossRef] [PubMed]
  24. Davidson, Z.E.; Walker, K.Z.; Truby, H. Do glucocorticosteroids alter vitamin D status? A systematic review with meta-analyses of observational studies. J. Clin. Endocrinol. Metab. 201497, 738–744. [Google Scholar] [CrossRef] [PubMed]
  25. Huybers, S.; Naber, T.H.J.; Bindels, R.J.M.; Hoenderop, J.G.J. Prednisolone-induced Ca2+ malabsorption is caused by diminished expression of the epithelial Ca2+ channel TRPV6. Am. J. Physiol.-Gastrointest. Liver Physiol. 2007292, 92–97. [Google Scholar] [CrossRef] [PubMed]
  26. Ferrari, P.; Bianchetti, M.G.; Sansonnens, A.; Frey, F.J. Modulation of renal calcium handling by 11β-hydroxysteroid dehydrogenase type 2. J. Am. Soc. Nephrol. 200213, 2540–2546. [Google Scholar] [CrossRef]
  27. Faggiano, A.; Pivonello, R.; Melis, D.; Filippella, M.; Di Somma, C.; Petretta, M.; Lombardi, G.; Colao, A. Nephrolithiasis in Cushing’s disease: Prevalence, etiopathogenesis, and modification after disease cure. J. Clin. Endocrinol. Metab. 200388, 2076–2080. [Google Scholar] [CrossRef] [PubMed]
  28. Ramsey, I.K.; Tebb, A.; Harris, E.; Evans, H.; Herrtage, M.E. Hyperparathyroidism in dogs with hyperadrenocorticism. J. Small Anim. Pract. 200546, 531–536. [Google Scholar] [CrossRef]
  29. Freiberg, J.M.; Kinsella, J.; Sacktor, B. Glucocorticoids increase the Na+-H+ exchange and decrease the Na+ gradient-dependent phosphate-uptake systems in renal brush border membrane vesicles. Proc. Natl. Acad. Sci. USA 198279, 4932–4936. [Google Scholar] [CrossRef] [PubMed]
  30. Sempos, C.T.; Heijboer, A.C.; Bikle, D.D.; Bollerslev, J.; Bouillon, R.; Brannon, P.M.; DeLuca, H.F.; Jones, G.; Munns, C.F.; Bilezikian, J.P.; et al. Vitamin D assays and the definition of hypovitaminosis D: Results from the First International Conference on Controversies in Vitamin D. Br. J. Clin. Pharmacol. 201884, 2194–2207. [Google Scholar] [CrossRef]
  31. Melnichenko, G.A.; Dedov, I.I.; Belaya, Z.E.; Rozhinskaya, L.Y.; Vagapova, G.R.; Volkova, N.I.; Grigor’ev, A.Y.; Grineva, E.N.; Marova, E.I.; Mkrtumayn, A.M.; et al. Cushing’s disease: The clinical features, diagnostics, differential diagnostics, and methods of treatment. Probl. Endocrinol. 201561, 55–77. [Google Scholar] [CrossRef]
  32. Machado, M.C.; De Sa, S.V.; Domenice, S.; Fragoso, M.C.B.V.; Puglia, P.; Pereira, M.A.A.; De Mendonça, B.B.; Salgado, L.R. The role of desmopressin in bilateral and simultaneous inferior petrosal sinus sampling for differential diagnosis of ACTH-dependent Cushing’s syndrome. Clin. Endocrinol. 200766, 136–142. [Google Scholar] [CrossRef]
  33. Findling, J.W.; Kehoe, M.E.; Raff, H. Identification of patients with Cushing’s disease with negative pituitary adrenocorticotropin gradients during inferior petrosal sinus sampling: Prolactin as an index of pituitary venous effluent. J. Clin. Endocrinol. Metab. 200489, 6005–6009. [Google Scholar] [CrossRef] [PubMed]
  34. Pigarova, E.A.; Rozhinskaya, L.Y.; Belaya, J.E.; Dzeranova, L.K.; Karonova, T.L.; Ilyin, A.V.; Melnichenko, G.A.; Dedov, I.I. Russian Association of Endocrinologists recommendations for diagnosis, treatment and prevention of vitamin D deficiency in adults. Probl. Endocrinol. 201662, 60–84. [Google Scholar] [CrossRef]
  35. Petrushkina, A.A.; Pigarova, E.A.; Tarasova, T.S.; Rozhinskaya, L.Y. Efficacy and safety of high-dose oral vitamin D supplementation: A pilot study. Osteoporos. Int. 201627, 512–513. [Google Scholar] [CrossRef]
  36. Pivonello, R.; De Martino, M.C.; De Leo, M.; Lombardi, G.; Colao, A. Cushing’s Syndrome. Endocrinol. Metab. Clin. N. Am. 200837, 135–149. [Google Scholar] [CrossRef]
  37. Belaya, Z.E.; Iljin, A.V.; Melnichenko, G.A.; Rozhinskaya, L.Y.; Dragunova, N.V.; Dzeranova, L.K.; Butrova, S.A.; Troshina, E.A.; Dedov, I.I. Diagnostic performance of late-night salivary cortisol measured by automated electrochemiluminescence immunoassay in obese and overweight patients referred to exclude Cushing’s syndrome. Endocrine 201241, 494–500. [Google Scholar] [CrossRef]
  38. Povaliaeva, A.; Pigarova, E.; Zhukov, A.; Bogdanov, V.; Dzeranova, L.; Mel’nikova, O.; Pekareva, E.; Malysheva, N.; Ioutsi, V.; Nikankina, L.; et al. Evaluation of vitamin D metabolism in patients with type 1 diabetes mellitus in the setting of cholecalciferol treatment. Nutrients 202012, 3873. [Google Scholar] [CrossRef] [PubMed]
  39. Thode, J.; Juul-Jørgensen, B.; Bhatia, H.M.; Kjaerulf-Nielsen, M.; Bartels, P.D.; Fogh-Andersen, N.; Siggaard-Andersen, O. Comparison of serum total calcium, albumin-corrected total calcium, and ionized calcium in 1213 patients with suspected calcium disorders. Scand. J. Clin. Lab. Investig. 198949, 217–223. [Google Scholar] [CrossRef]
  40. Bikle, D.D.; Siiteri, P.K.; Ryzen, E.; Haddad, J.G.; Gee, E. Serum protein binding of 1, 25-Dihydroxyvitamin D: A reevaluation by direct measurement of free metabolite levels. J. Clin. Endocrinol. Metab. 198561, 969–975. [Google Scholar] [CrossRef]
  41. Bikle, D.D.; Gee, E.; Halloran, B.; Kowalski, M.A.N.N.; Ryzen, E.; Haddad, J.G. Assessment of the Free Fraction of 25-Hydroxyvitamin D in Serum and Its Regulation by Albumin and the Vitamin D-Binding Protein. J. Clin. Endocrinol. Metab. 198663, 954–959. [Google Scholar] [CrossRef]
  42. Holick, M.F.; Binkley, N.C.; Bischoff-Ferrari, H.A.; Gordon, C.M.; Hanley, D.A.; Heaney, R.P.; Murad, M.H.; Weaver, C.M. Evaluation, treatment, and prevention of vitamin D deficiency: An Endocrine Society clinical practice guideline. J. Clin. Endocrinol. Metab. 201196, 1911–1930. [Google Scholar] [CrossRef] [PubMed]
  43. Dirks, N.F.; Martens, F.; Vanderschueren, D.; Billen, J.; Pauwels, S.; Ackermans, M.T.; Endert, E.; den Heijer, M.; Blankenstein, M.A.; Heijboer, A.C. Determination of human reference values for serum total 1,25-dihydroxyvitamin D using an extensively validated 2D ID-UPLC–MS/MS method. J. Steroid Biochem. Mol. Biol. 2016164, 127–133. [Google Scholar] [CrossRef] [PubMed]
  44. Tang, J.C.Y.; Nicholls, H.; Piec, I.; Washbourne, C.J.; Dutton, J.J.; Jackson, S.; Greeves, J.; Fraser, W.D. Reference intervals for serum 24,25-dihydroxyvitamin D and the ratio with 25-hydroxyvitamin D established using a newly developed LC–MS/MS method. J. Nutr. Biochem. 201746, 21–29. [Google Scholar] [CrossRef]
  45. Máčová, L.; Bičíková, M. Vitamin D: Current challenges between the laboratory and clinical practice. Nutrients 202113, 1758. [Google Scholar] [CrossRef] [PubMed]
  46. Ginsberg, C.; Hoofnagle, A.N.; Katz, R.; Hughes-Austin, J.; Miller, L.M.; Becker, J.O.; Kritchevsky, S.B.; Shlipak, M.G.; Sarnak, M.J.; Ix, J.H. The Vitamin D Metabolite Ratio Is Associated With Changes in Bone Density and Fracture Risk in Older Adults. J. Bone Miner. Res. 2021, 1–8. [Google Scholar] [CrossRef]
  47. Cavalier, E.; Huyghebaert, L.; Rousselle, O.; Bekaert, A.C.; Kovacs, S.; Vranken, L.; Peeters, S.; Le Goff, C.; Ladang, A. Simultaneous measurement of 25(OH)-vitamin D and 24,25(OH)2-vitamin D to define cut-offs for CYP24A1 mutation and vitamin D deficiency in a population of 1200 young subjects. Clin. Chem. Lab. Med. 202058, 197–201. [Google Scholar] [CrossRef]
  48. Rondeau, P.; Bourdon, E. The glycation of albumin: Structural and functional impacts. Biochimie 201193, 645–658. [Google Scholar] [CrossRef]
  49. Soudahome, A.G.; Catan, A.; Giraud, P.; Kouao, S.A.; Guerin-Dubourg, A.; Debussche, X.; Le Moullec, N.; Bourdon, E.; Bravo, S.B.; Paradela-Dobarro, B.; et al. Glycation of human serum albumin impairs binding to the glucagon-like peptide-1 analogue liraglutide. J. Biol. Chem. 2018293, 4778–4791. [Google Scholar] [CrossRef]
  50. McLeod, J.F.; Kowalski, M.A.; Haddad, J.G. Interactions among serum vitamin D binding protein, monomeric actin, profilin, and profilactin. J. Biol. Chem. 1989264, 1260–1267. [Google Scholar] [CrossRef]
  51. Gupta, Y.; Gupta, A. Glucocorticoid-induced myopathy: Pathophysiology, diagnosis, and treatment. Indian J. Endocrinol. Metab. 201317, 913. [Google Scholar] [CrossRef] [PubMed]
  52. Smets, P.; Meyer, E.; Maddens, B.; Daminet, S. Cushing’s syndrome, glucocorticoids and the kidney. Gen. Comp. Endocrinol. 2010169, 1–10. [Google Scholar] [CrossRef] [PubMed]
  53. Arnaud, J.; Constans, J. Affinity differences for vitamin D metabolites associated with the genetic isoforms of the human serum carrier protein (DBP). Hum. Genet. 199392, 183–188. [Google Scholar] [CrossRef] [PubMed]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Updated Cushing’s disease guideline highlights new diagnosis, treatment ‘roadmap’

An updated guideline for the treatment of Cushing’s disease focuses on new therapeutic options and an algorithm for screening and diagnosis, along with best practices for managing disease recurrence.

Despite the recent approval of novel therapies, management of Cushing’s disease remains challenging. The disorder is associated with significant comorbidities and has high mortality if left uncontrolled.

Adrenal transparent _Adobe
Source: Adobe Stock

“As the disease is inexorable and chronic, patients often experience recurrence after surgery or are not responsive to medications,” Shlomo Melmed, MB, ChB, MACP, dean, executive vice president and professor of medicine at Cedars-Sinai Medical Center in Los Angeles, and an Endocrine Today Editorial Board Member, told Healio. “These guidelines enable navigation of optimal therapeutic options now available for physicians and patients. Especially helpful are the evidence-based patient flow charts [that] guide the physician along a complex management path, which usually entails years or decades of follow-up.”

Shlomo Melmed

The Pituitary Society convened a consensus workshop with more than 50 academic researchers and clinical experts across five continents to discuss the application of recent evidence to clinical practice. In advance of the virtual meeting, participants reviewed data from January 2015 to April 2021 on screening and diagnosis; surgery, medical and radiation therapy; and disease-related and treatment-related complications of Cushing’s disease, all summarized in recorded lectures. The guideline includes recommendations regarding use of laboratory tests, imaging and treatment options, along with algorithms for diagnosis of Cushing’s syndrome and management of Cushing’s disease.

Updates in laboratory, testing guidance

If Cushing’s syndrome is suspected, any of the available diagnostic tests could be useful, according to the guideline. The authors recommend starting with urinary free cortisol, late-night salivary cortisol, overnight 1 mg dexamethasone suppression, or a combination, depending on local availability.

If an adrenal tumor is suspected, the guideline recommends overnight dexamethasone suppression and using late-night salivary cortisol only if cortisone concentrations can also be reported.

The guideline includes several new recommendations in the diagnosis arena, particularly on the role of salivary cortisol assays, according to Maria Fleseriu, MD, FACE, a Healio | Endocrine Today Co-editor, professor of medicine and neurological surgery and director of the Pituitary Center at Oregon Health & Science University in Portland.

Maria Fleseriu

“Salivary cortisol assays are not available in all countries, thus other screening tests can also be used,” Fleseriu told Healio. “We also highlighted the sequence of testing for recurrence, as many patients’ urinary free cortisol becomes abnormal later in the course, sometimes up to 1 year later.”

The guideline states combined biochemical and imaging for select patients could potentially replace petrosal sinus sampling, a very specialized procedure that cannot be performed in all hospitals, but more data are needed.

“With the corticotropin-releasing hormone stimulation test becoming unavailable in the U.S. and other countries, the focus is now on desmopressin to replace corticotropin-releasing hormone in some of the dynamic testing, both for diagnosis of pseudo-Cushing’s as well as localization of adrenocorticotropic hormone excess,” Fleseriu said.

The guideline also has a new recommendation for anticoagulation for high-risk patients; however, the exact duration and which patients are at higher risk remains unknown.

“We always have to balance risk for clotting with risk for bleeding postop,” Fleseriu said. “Similarly, recommended workups for bone disease and growth hormone deficiency have been further structured based on pitfalls specifically related to hypercortisolemia influencing these complications, as well as improvement after Cushing’s remission in some patients, but not all.”

New treatment options

The guideline authors recommended individualizing medical therapy for all patients with Cushing’s disease based on the clinical scenario, including severity of hypercortisolism. “Regulatory approvals, treatment availability and drug costs vary between countries and often influence treatment selection,” the authors wrote. “However, where possible, it is important to consider balancing cost of treatment with the cost and the adverse consequences of ineffective or insufficient treatment. In patients with severe disease, the primary goal is to treat aggressively to normalize cortisol concentrations.”

Fleseriu said the authors reviewed outcomes data as well as pros and cons of surgery, repeat surgery, medical treatments, radiation and bilateral adrenalectomy, highlighting the importance of individualized treatment in Cushing’s disease.

“As shown over the last few years, recurrence rates are much higher than previously thought and patients need to be followed lifelong,” Fleseriu said. “The role of adjuvant therapy after either failed pituitary surgery or recurrence is becoming more important, but preoperative or even primary medical treatment has been also used more, too, especially in the COVID-19 era.”

The guideline summarized data on all medical treatments available, either approved by regulatory agencies or used off-label, as well as drugs studied in phase 3 clinical trials.

“Based on great discussions at the meeting and subsequent emails to reach consensus, we highlighted and graded recommendations on several practical points,” Fleseriu said. “These include which factors are helpful in selection of a medical therapy, which factors are used in selecting an adrenal steroidogenesis inhibitor, how is tumor growth monitored when using an adrenal steroidogenesis inhibitor or glucocorticoid receptor blocker, and how treatment response is monitored for each therapy. We also outline which factors are considered in deciding whether to use combination therapy or to switch to another therapy and which agents are used for optimal combination therapy.”

Future research needed

The guideline authors noted more research is needed regarding screening and diagnosis of Cushing’s syndrome; researchers must optimize pituitary MRI and PET imaging using improved data acquisition and processing to improve microadenoma detection. New diagnostic algorithms are also needed for the differential diagnosis using invasive vs. noninvasive strategies. Additionally, the researchers said the use of anticoagulant prophylaxis and therapy in different populations and settings must be further studied, as well as determining the clinical benefit of restoring the circadian rhythm, potentially with a higher nighttime medication dose, as well as identifying better markers of disease activity and control.

“Hopefully, our patients will now experience a higher quality of life and fewer comorbidities if their endocrinologist and care teams are equipped with this informative roadmap for integrated management, employing a consolidation of surgery, radiation and medical treatments,” Melmed told Healio.