MDPI

Article A Comprehensive Investigation of Steroidogenic Signaling in Classical and New Experimental Cell Models of Adrenocortical Carcinoma

Sandra Sigala 1,+D, Christina Bothou 2,+D, David Penton 3, Andrea Abate 1D, Mirko Peitzsch 4D ,

Deborah Cosentini 5, Guido A. M. Tiberio 6D, Stefan R. Bornstein 2,4,7,8,9,10, Alfredo Berruti 5,#[D and Constanze Hantel 2,4,*,

check for updates

Citation: Sigala, S .; Bothou, C .; Penton, D .; Abate, A .; Peitzsch, M .; Cosentini, D .; Tiberio, G.A.M .;

Bornstein, S.R .; Berruti, A .; Hantel, C. A Comprehensive Investigation of Steroidogenic Signaling in Classical and New Experimental Cell Models of Adrenocortical Carcinoma. Cells 2022, 11, 1439. https://doi.org/ 10.3390/cells11091439

Academic Editor: Marcin Rucinski

Received: 14 March 2022

Accepted: 22 April 2022

Published: 24 April 2022

Corrected: 14 September 2023

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affil- iations.

CC

İ

BY

Copyright: @ 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

1 Section of Pharmacology, Department of Molecular and Translational Medicine, University of Brescia, 25124 Brescia, Italy; sandra.sigala@unibs.it (S.S.); a.abate005@unibs.it (A.A.)

2 Department of Endocrinology, Diabetology and Clinical Nutrition, University Hospital Zurich (USZ) and University of Zurich (UZH), 8091 Zürich, Switzerland; christina.bothou@usz.ch (C.B.); stefan.bornstein@uniklinikum-dresden.de (S.R.B.)

3 Electrophysiology Facility (e-phac), Department of Molecular Life Sciences, University of Zurich (UZH), 8057 Zürich, Switzerland; david.pentonribas@uzh.ch

4 Medizinische Klinik und Poliklinik III, University Hospital Carl Gustav Carus Dresden, 01307 Dresden, Germany; mirko.peitzsch@uniklinikum-dresden.de

5 Medical Oncology Unit, Department of Medical and Surgical Specialties, Radiological Sciences, and Public Health, University of Brescia at ASST Spedali Civili di Brescia, 25124 Brescia, Italy; deborah.cosentini@gmail.com (D.C.); alfredo.berruti@unibs.it (A.B.)

6 Surgical Clinic, Department of Clinical and Experimental Sciences, University of Brescia at ASST Spedali Civili di Brescia, 25124 Brescia, Italy; guido.tiberio@unibs.it

7 Diabetes and Nutritional Sciences, King’s College London, London WC2R 2LS, UK

8 Center for Regenerative Therapies, Technische Universität Dresden, 01307 Dresden, Germany

9 Paul-Langerhans-Institute Dresden, Helmholtz Center Munich, University Hospital Carl Gustav Carus, Faculty of Medicine, Technische Universität Dresden, 01307 Dresden, Germany

10 Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore 636921, Singapore

* Correspondence: constanze.hantel@usz.ch; Tel .: +41-43-253-3008

t These authors contributed equally to this work.

# These authors contributed equally to this work.

Abstract: Adrenocortical carcinoma is a heterogeneous and aggressive cancer that originates from steroidogenic cells within the adrenal cortex. In this study, we have assessed for the preclinical gold standard NCI-H295 in direct comparison with the more recently established MUC-1 and a here newly reported ACC cell line (TVBF-7) the mutational status of important driver genes (TP53, MEN1, PRKAR1A, CTNNB1, APC, ZNRF-3, IGF-2, EGFR, RB1, BRCA1, BRCA2, RET, GNAS and PTEN), Wnt-signaling specificities (CTNNB1 mutation vs. APC mutation vs. wildtype), steroidogenic- (CYP11A1, CYP17A1, HSD3B2, HSD17B4, CYP21A2, CYP11B1,CYP11B2, MC2R, AT1R) and nuclear- receptor-signaling (AR, ER, GCR), varying electrophysiological potentials as well as highly individual hormone secretion profiles (Cortisol, Aldosterone, DHEA, DHEAS, Testosterone, 17-OH Progesterone, among others) which were investigated under basal and stimulated conditions (ACTH, AngII, FSK). Our findings reveal important genetic and pathophysiological characteristics for these three cell lines and reveal the importance of such cell-line panels reflecting differential endocrine functionalities to thereby better reflect clinically well-known ACC patient heterogeneities in preclinical studies.

Keywords: adrenocortical carcinoma cell lines; steroidogenesis; electrophysiology; genotype; NCI-H295; MUC-1; TVBF-7

1. Introduction

Adrenocortical carcinoma (ACC) is a rare and aggressive cancer that originates from steroidogenic cells within the adrenal cortex with an estimated incidence of 0.7-1 cases

per million of the population per year [1]. ACCs demonstrate heterogeneous steroid secretion patterns producing hormones such as mineralocorticoids, glucocorticoids and androgens or secrete precursor metabolites, deriving from intermediate steps along the three major adrenocortical biosynthetic pathways ([2], Figure 1). Only in rare cases ACCs are biochemically inactive and the elucidation of the steroid profile of ACCs is important in diagnostic, but also in therapeutic level.

Figure 1. Schematic representation of steroidogenesis pathway from stimulation to hormone pro- duction. Figures have been created modifying an image set from Servier Medical Art (SMART) http://smart.servier.com/ (accessed on 19 July 2021).

Stimulation

K+ CI-

Angiotensin II

ACTH

Forskolin

CaV.

AT1R

MC2R

GNRHR

Ca2+

Secondary Messengers

Hormones

ER HSP

AR HSP

GCR HSP

Transcription Factors

Steroidogenesis

SF-1

Mitochondria

smooth ER

Steroidogenic Enzymes

Nucleus

Glucocorticoid precursors

HSD11B1

11-Dehydro corticosterone

Mineralocorticoids

Cholesterol

StAR

Pregnenolone

HSD3B2

Progesterone

CYP21A2

11-Deoxy corticosterone

CYP11B2

Corticosterone

CYP11A1

CYP1182

18-OH Corticosterone

CYP11B2

Aldosterone

CYP17A1

CYP17A1

11-Deoxy cortisol

CYP11B1

HSD1181/2

Glucocorticoids

Cortisone

17-OH Pregnenolone

HSD3B2

17-OH Progesterone

CYP21A2

Cortisol

CYP11B1

21-Deoxy cortisol

CYP21A2

CYP11B2

18-OH Cortisol

CYP11B2

18-Oxo- cortisol

CYP17A1

Androgen precursors

Androgens

DHEA

HSD3B2

Androstenedione

HSD17B

Testosterone

5a-reductase

Dihydro testosterone

SULT2A1

CYP19

CYP19

Estrogens

HSD17B

Estrone

Estradiol

HSD17B

DHEAS

Estriol

In recent years, steroid profiling of patient plasma or 24-h urine samples have been de- veloped and provide steroid metabolomic information via liquid chromatography tandem mass spectrometry (LC-MS/MS) or gas chromatography mass spectrometry (GC-MS) [3,4]. Those tools combined with machine-learning based approaches have resulted in the identifi- cation of distinct malignant steroid “fingerprints” for ACC, which could differentiate benign from malignant adrenal tumor [2,3], but also be used as biomarkers and as a screening tool for early identification of tumor recurrence [5]. On the other hand, steroidogenesis is the direct or indirect target of action of drugs used as chemotherapeutics in patients diagnosed with ACC, with main representative the gold standard medication: mitotane [6,7].

To date preliminary data obtained for mineralocorticoids and glucocorticoids show that the steroidogenic pattern could be related with the disease outcome [8-10]. However, detailed mechanistic insights are lacking and as indicated above, also androgens and androgen receptor signalling is commonly present in ACC. The androgen receptor (AR) is a ligand-activated transcription factor that plays an important role in the context of various severe diseases. The AR is already a therapeutic target for prostate cancer and is also emerging as a new marker and potential therapeutic target for breast cancer. The availability of selective AR inhibitors (e.g., bicalutamide, enzalutamide, apalutamide) approved for the treatment of prostate cancer might have, thus, potential to be translated to other endocrine cancers. However, for breast and prostate cancer AR-signalling appears to have different functions according to the specific subtype [11], and disease stage [12]. Steroid hormones other than the classical ligands testosterone and dihydrotestosterone may be in this context also of interest considering their potential role in the development of resistance mechanisms [12].

However, important pre-requisite for detailed mechanistic and therapeutic studies are preclinical models. Considering that primary cell cultures from adrenal tumors are, due to rarity of this tumor type, overall limited, many groups have attempted to establish cell lines from ACCs and currently different ACC cell lines are meanwhile available [13-17]. The most implemented model for the study of steroidogenic gene expression and chemotherapeutic responsiveness in ACC until today is NCI-H295, which was established in the early 1990s from Gazdar et al. from a female patient diagnosed with a primary ACC [14]. Extensive in vitro studies have shown that the original line NCI-H295 expresses in baseline all the enzymes participating in normal human adrenal steroidogenic gene expression, including the enzymes that catalyse the rate- limiting steroidogenic step (STAR and CYP11A1) and all the major biosynthetic steroidogenic enzymes (CYP17, HSD3B2, CYP21, CYP11B2, CYP11B1 and HSD17B4) and retain, thus, the general potential to produce all major adrenal steroids. Considering that adrenocortical cells are, normally, excitable cells, the precise control of the membrane voltage is very important for the initiation of steroid synthesis. However, in context of steroidogenesis it should be also mentioned that clinical ACC cases of combined mineralocorticoid and glucocorticoid secretion, as reflected by NCI-H295, are comparably rare. Thus, from a clinical point of view, these cells do maybe not represent the most common phenotype in terms of ACC functionality [18]. Likewise, combined mutations in the two ACC driver genes TP53 and CTNNB1, as described for NCI-H295, is for ACC clinically rather rarely observed [19]. However, due to a lack of ACC models for many years, comprehensive data from additional human models reflecting diverse subtypes are nowadays for many aspects still lacking.

A recently developed, highly emerging model in this field is MUC-1, which has been established as patient-derived tissue xenograft and cell line from an ACC neck metastasis of a male patient [15]. MUC-1 was previously presented to be SF-1 and 3ßHSD positive and mice bearing MUC-1 xenografts had increased plasma cortisol [15]. Meanwhile, MUC-1 has been implemented in a variety of preclinical ACC studies and demonstrated repeatedly a different response pattern and clinically frequently observed drug resistance phenotype compared to NCI-H295R [20-29].

Here, we describe the development and characterization of a new ACC cell line, named TVBF-7, which is carrying a nonsense APC and represents functional signs of autonomous cortisol secretion. In a comprehensive study, we investigate TVBF-7in direct comparison

to NCI-H295 and MUC-1 regarding mutational status of main driver genes, steroidogenic signaling, electrophysiological properties as well as secretion profiles and reveal patterns of genetical, steroidogenic distinct and highly relevant pathophysiological ACC sub-types.

2. Materials and Methods

2.1. Cancer Cell Lines

NCI-H295R cells were originally obtained from ATCC while MUC-1 cells were previ- ously established by our group [15]. Both cell lines were authenticated again in December 2020 by Microsynth (Balgach, Switzerland) and maintained as previously described [15]. Additionally, the novel cell line described as TVBF-7 was authenticated for multiple pas- sages at the universities Zurich and Brescia (Balgach, Switzerland, and BMR Genomics, Padova, Italy). The cells were cultured in advanced D-MEM-F12 (#12634010, Gibco, Waltham, MS, USA), 10% FBS (#10082147, Gibco), 1% penicillin-streptomycin (#15070063, Gibco), 1% amphotericin B (#15290026, Gibco) and 2 mM glutamine (#A2916801, Gibco).

2.2. Whole Genome Sequencing and Genomic Data Presentation

Whole genome sequencing (WGS) for NCI-H295R, MUC-1 and TVBF-7 cells has been performed and standard bioinformatic analysis has been carried out (BGI, Shenzhen, Guang- dong, China). The filtered reads were aligned to the human reference genome (UCSC build HG19) using Burrows-Wheeler Aligner (BWA) software. Single Nucleotide Polymorphisms (SNP), insertions and deletions (InDel) and copy number variations (CNVs) have been annotated (SnpEff tool: http://snpeff.sourceforge.net/SnpEff_manual.html (accessed on 19 July 2021) and VEP tool: https://asia.ensembl.org/info/docs/tools/vep/index.html (accessed on 19 July 2021)) and filtered as presented in detail in the Supplementary File. The data were delivered post-analysis in vcf files, which were further assessed by the use of EmEditor (Washington, DC, USA) for the generation of the information provided in Figure 2B and in the Supplementary Table S3.

A

TVBF-7 STR Profiling
Marker519631 - D1
AMELx
D3S135817
D1S165616
D2S44112
D10S124814
D13S31712
Penta E5
D1655398
D1855119
D2S133819
CSF1PO12
Penta D9 12
TH016
vWA17
D2151129
D7582011 12
D5581812
TPOX8
DYS391-
D85117915
D12539118 20
D19543313
FGA21
D22S104515

B

Gene NameCell lineType of alterationMutation
TP53NCI-H295Rlarge deletionr.783_993del
TP53MUC-1deletion/ frameshiftNM_000546: c.1024deIC (p.R342fs)
CTNNB1NCI-H295RmissenseNM_001098209: c.T133C (p.S45P)
APCTVBF-7nonsenseNM_000038.5: c.739C>T (p.Q247*)
BRCA2NCI-H295RmissenseNM_000059: c.C9349G (p.H3117D)

Figure 2. STR profiling of the newly presented TVBF-7 cells (A). Most important pathogenic mutations of cancer driver genes as detected by the whole genome sequencing of all three cell lines (B).

2.3. Stimulation Experiments

On the treatment starting day, the cell medium was removed, and the cells were washed with PBS and then treated with the following agents: potassium chloride (KCI) (5, 10, 20 mM concentrations added on top to the baseline K+ concentration contained in the medium calculated to around 5 mM), (#P5405, Sigma-Aldrich, Buchs, Switzerland), AngII

,

(20, 100, 250 nM) (#A9525, Sigma-Aldrich), ACTH (5, 15, 25 nM) (#A0423, Sigma-Aldrich) or FSK (5, 10, 25 LM) (#F3918, Sigma-Aldrich). Each concentration of the samples and each control was included in triplicate. The plates were incubated at 37 ℃ and 5% CO2 for 24 h, and then either processed for quantitative Real-Time PCR in 24-well plates (TPP #92024; NCI-H295R (200,000/well), MUC-1 cells (85,000/well), TVBF-7 (180,000/well)) or for steroid measurements in the cell supernatant with simultaneous protein quantifi- cation within the cells in 6-well plates (TPP#92006; NCI-H295R (1,000,000/well), MUC-1 (425,000/well) and TVBF-7 (750,000/well)).

2.4. Quantitative Real-Time PCR

After 24h of stimulation, cells were prepared for total RNA isolation, genomic DNA removal and cDNA generation as previously described [20]. For real-time PCR analy- sis, EvaGreen® reaction mix (#1725200, Bio-Rad, Hercules, CA, USA) in QuantStudio5 (applied biosystems, Waltham, MS, USA) was used. The primers used are described in Supplementary Table S2. Differences in the threshold cycle (Ct) values between the GAPDH housekeeping gene and the gene of interest (Ct) were then calculated as an indicator of difference in the amount of mRNA expressed, corrected for the efficiency of the reaction previously acquired via standard curve.

2.5. Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS) Steroid Measurements

After 24h of stimulation, cell supernatant was collected, directly frozen on dry ice and stored at. - 80 ℃. A whole panel of steroid metabolomics was determined by LC-MS/MS as described previously (Peitzsch et al., 2015). The results were provided in concentration (ng/ml). As blank the respective medium supernatant for similarly treated but cell free wells was used for untreated and stimulated conditions. For this experiment, the following concentrations have been used: KCl 10 mM, AngII 100 nM, FSK 10 µM and a low (l) and a high ACTH (h) concentration corresponding to 10 and 50 nM, respectively.

2.6. Protein Quantification

Protein samples from the wells used for supernatant collection were acquired, following the same order. More specifically, cell proteins were extracted in RIPA buffer (50 mM Tris pH 8.0, 150 mM NaCl, 0.01 v/v NP-40 #74385, Sigma-Aldrich, St. Louis, MI, USA), 0.005 v/v sodium deoxycholate (#D6750, Sigma-Aldrich), and 0.001 w/v SDS (#2326.3, Roth, Karlsruhe, Germany) supplemented with a complete protease inhibitor cocktail (#11836170001, Roche, NY, USA) and phosphatase inhibitor cocktail (#P5726, Sigma-Aldrich). The homogenized lysate was centrifuged at 16,000x g for 15 min and protein concentration was quantified by Pierce BCA Protein Assay (#23225, Thermo Scientific™M, Reinach, Switzerland) following the manufacturer’s recommendations.

2.7. Electrophysiological Studies

For whole-cell automated patch clamp, NCI-H295R, MUC-1 and TVBF-7 cells were grown until ~80% confluency. Cells were used immediately after detaching with Accutase (Sigma-Aldrich, Buchs, Switzerland) and resuspended in an extracellular solution contain- ing (in mM) 135 NaCl, 1.8 MgCl2, 1.8 CaCl2, 10 HEPES, 5 KCl, pH 7.4 adjusted NaOH/HCl. 3.0 x 106 cells/mL and then directly added to the centrifuge tube of the QPatchII auto- mated patch clamp platform (Sophion Bioscience, Ballerup, Denmark). 48X single hole QChips with a resistance of ~2 MQ were used for experiments with an intracellular solu- tion containing (in mM) 95 K-gluconate, 30 KCl, 4.8 Na2HPO4, 1.2 NaH2PO4, 5 glucose, 2.38 MgCl2, 0.73 CaCl2, 1 EGTA, 3 Mg-ATP, pH 7.2 adjusted KOH/HCl. No seal enhancer and no liquid junction potential correction (11 mV measured for the solution pair) was used. Only cells maintaining a seal resistance of >0.5 G2, and with stable access resistance throughout the experiment, were used. After obtaining the whole-cell patch clamp, the extracellular solution was exchanged twice with solution containing 5, 15, 50 mM K+ and a protocol including 20 mv voltage steps (-100 to +80 mV, holding potential -60 mV, 200 ms)

was applied every 50 s. The osmolality of the extracellular solution was kept constant by proportionally decreasing the concentration of NaCl for increasing concentrations of KCl (e.g., in mM 135 NaCl + 5 KCl or 90 NaCl + 50 KCl). Current density was calculated using steady state current (average of last 100 ms of each step) divided by the cell capacitance measured immediately before each protocol. For the effect of AngII (20 nM), FSK (5 LM) and ACTH (5 nM), cells were washed twice with extracellular solution containing 5 mM K+ before the application of the effector. Voltage protocols including 20 mv voltage steps (-100 to +80 mV, holding potential -60 mV, 200 ms) was applied. For each treatment, currents and reversal potential were measured.

2.8. Statistical Analysis and Graphical Designs

Statistical analysis and graphical representation of the data was carried out using GraphPad Prism software (version 8, GraphPad Software, La Jolla, CA, USA). If not stated otherwise, comparison between control group and two or more treatment groups or between cell lines (mean of each) were performed by one-way ANOVA followed by Dunnett’s multiple comparisons test. The data are presented in column graphs depicting the mean ± SEM. The statistical significance is denoted as stars in the graphs (* p < 0.05; ** p < 0.01; *** p < 0.001).

For the representation of the graphical abstract of Figure 1, the figures created by modifying an image set from Servier Medical Art (SMART) http://smart.servier.com/ (accessed on 19 July 2021), which is cited appropriately.

3. Results

3.1. Establishment of a Novel Cancer Cell Line, TVBF-7

Recently, ACC primary cells (primary culture ACC115m [27]) have been isolated from a lymph node metastasis from a male patient without clinically obvious signs of steroid excess, while no functional testing has been performed. When the cells were found to be continuously passageable, the cells were further cultured for multiple passages at the university of Brescia (to date P 34) and remained stable as confirmed by cell authentication performed via STR profiling in passages 7, 14, 16, 23, and 29 (Supplementary Table S1). After transfer to the University of Zurich, the cells were cultured during the subsequently described experiments and again authenticated by STR-Analysis (Figure 2A). The newly established and characterized ACC cell line has been named TVBF-7.

3.2. Mutational Status of Important Driver Genes

The mutational status of TVBF-7 in main driver genes (TP53, MEN1, PRKAR1A, CTNNB1, APC, ZNRF-3, IGF-2, EGFR, RB1, BRCA1, BRCA2, RET, GNAS and PTEN) has been assessed vs. NCI-H295R and MUC-1. In Figure 2B, the most critical mutations are presented, a full panel of findings is provided in the Supplementary Materials. The main differences were detected regarding Wnt-signaling pathway regulators, as our analyses revealed for TVBF-7 a nonsense mutation in APC (Figure 2B). For NCI-H295R, the known CTNNB1 mutation was confirmed while MUC-1 represented the wild type for both genes. Moreover, the analysis revealed TP53 WT for TVBF-7 compared to the TP53 mutated genotypes of NCI-H295R and MUC-1. For BRCA2 a missense mutation was found for NCI-H295R exclusively.

3.3. Baseline Gene Expression Levels and Electrophysiological Properties

Interestingly, comparative characterizations of baseline gene expressions revealed for TVBF-7 much higher levels of Melanocortin 2 Receptor (MC2R) compared to NCI-H295R (NCI-H295R: 100.0 ± 3.9% vs. MUC-1: 0.4 ± 0.1% vs. TVBF-7: 301.4 ± 18.1, p < 0.001 for all comparisons), but lower for other hormonal receptors (Angiotensin II receptor type 1-AT1R, Estrogen Receptor 1-ER1, Androgen Receptor-AR and Gonadotropin-Releasing-Hormone Receptor-GNRHR). MUC-1 demonstrate overall rather low or intermediate (Glucocorticoid receptor) baseline expression (Figure 3A). In accordance, also CYP11B1-expression was

extraordinary high in TVBF-7 (NCI-H295R: 100.0 ± 3.2%, MUC-1: 0.0 ± 0.0% mV and TVBF-7: 20555.1 ± 991.4%, p < 0.001 for all comparisons) and CYP11A1, while for MUC-1 the basal gene expression was found to be overall low or undetectable (Figure 3B). For CYP17A1, HSD3B2, HSD17B4, CYP21A2 and CYP11B2 baseline gene-expression was demonstrated to be highest in NCI-H295R.

Figure 3. Baseline gene expression levels of different hormone receptors (MC2R, AT1R, AR, GNRH, ER1 and GCR) of unstimulated cells (A). Baseline gene expression levels of steroidogenic enzymes (CYP11A1, CYP17A1, HSD3B2, HSD17B4, CYP21A2, CYP11B2 and CYP11B1) (B). Comparative depiction of the reversal potential in baseline (5mM KCl) and upon stimulation with total 15 and 50 mM KCl for all three

A

B

Hormone Receptor Baseline Expression

reversal potential

NCI-H295R

400

mM KCI

Hormone Receptor/GAPDH

MUC-1

5

15

50

5

15

50

5

15

50

300

TVBF-7

0

gene expression

200

150

reversal potential

-10

(mV)

-20

100

NCI-H295R

-30

MUC-1

-40-

TVBF-7

50

-50

0

MC2R

AT1R

AR

GNRHR

ER1

GCR

C

D

Steroidogenic Enzyme Baseline Expression

25000

NCI-H295R

Currents at 80mV

80

ns

ns

ns

Ster. Enzyme/GAPDH gene expression

15000

MUC-1

250

TVBF-7

current density (pA/pF)

NCI-H295R

60

MUC-1

200

TVBF-7

150

40

100

20

50

0

5

15

50

5

15 50

5

15

50

0

CYP11A1

CYP17A1

HSD3B2

HSD17B4

CYP21A2

CYP11B2

CYP11B1

mM KCI

E

CYP11A1

F

HSD3B2

1500

NCI-H295

1300

NCI-H295

CYP11A1/GAPDH gene expression

1000

MUC-1

1100

TVBF-7

HSD3B2/GAPDH

gene expression

900

MUC-1

500

700

TVBF-7

150

500

100

200

50

100

0

0

0

5

10

25

0

5

10

25

0

5

10

25

0

5

10

25

0

5

10

25

0

U

10

25

μM FSK

µM FSK

G

H

CYP17A1

17-OH Progesterone

12000

200

NCI-H295

NCI-H295

CYP17A1/GAPDH

gene expression

10000

MUC-1

17 OH-Progesterone (ng/mg protein)

150

MUC-1

8000

TVBF-7

100

75

TVBF-7

6000

300

200

35

100

0

0

5

10

25

0

5

10

25

0

5

10

25

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (I)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

UM FSK

cell lines (C). Mean current densities in baseline (5mM KCl) and upon stimulation with total 15 and 50mM KCl for NCI-H295R (n = 51), MUC-1 (n = 39) and TVBF-7 (n = 26) cells at 80 mV (D). Stimulation of the first steroidogenesis steps among the different cell lines and upon different stimulations. CYP11A1 (E), HSD3B2 (F) and CYP17A1 (G) gene expression levels upon FSK stimulation per cell line. Hormonal production of 17-OH progesterone (H) in ng per mg of total protein for unstimulated cells in comparison with KCl, AngII, FSK and two different ACTH stimulations. Stars represent significance vs. untreated cells ( ** p < 0.01; *** p < 0.001).

Next, the baseline cell sensitivity to potassium (K+) was tested electrophysiologically. All cell lines responded to increased concentrations of KCl with depolarization, as shown by the changes in the reversal potential, used as a marker of the resting membrane voltage (Figure 3C, NCI-H295R (n=28-43), MUC-1 (n=25-28) and TVBF-7 (n=24-28). In the overall presentation for all cell lines, these results are indicated as not-significant for MUC-1, but when the data were analyzed with paired t-tests and not as pooled means, the significances were highlighted also in case of MUC-1 (data not shown). Of note, TVBF-7 tend to be the most hyperpolarized (NCI-H295R: - 32.11 ± 4.91 mV, MUC-1: - 34.74 ± 3.42 mV and TVBF-7: 41.06 ± 4.52, p > 0.05 for all comparisons). Considering that, in the presenting data, the liquid junction potential has not been subtracted (calculated around 11 mV for the current setting), the resting membrane voltage of the TVBF-7 cells demonstrates overall the closest value to the known hyperpolarized standard for adrenal cells (around -80 mV). For the same experimental setup, the membrane currents for different KCl concentrations tested at 80 mV remained insignificant for NCI-H295R (n =51), MUC-1 (n = 39) and TVBF-7 (n =25) cells (Figure 3D).

3.4. Stimulation of the First Steps of Steroidogenesis upon Known Stimuli: TVBF-7 Are Unresponsive

Next, we investigated if the cells respond to known stimulators of steroidogenesis such as KCl, AngII, ACTH and FSK, using FSK as a known inducer of overall steroidogenesis by increasing intracellular cAMP levels. Indeed, FSK stimulation resulted in induction of gene expression of the rate-limiting and first main steps of steroidogenesis (CYP11A1, HSD3B2 and CYP17A1) in NCI-H295R and also MUC-1 cells, while TVBF-7 remained unresponsive (Figure 3E-G). Interestingly, MUC-1 cells responded upon FSK stimulation here even more intensely than NCI-H295R cells. In agreement with these findings, the investigation of 17-OH- Progesterone detected by LC-MS/MS revealed significant increases for NCI-H295R and MUC-1, but not for TVBF-7. Progesterone levels increased for MUC-1 only (NCI-H295R: untreated 0.0 ± 0 vs. FSK 0.0 ± 0 ng/mg protein; MUC-1: untreated 1.6± 0.03 vs. FSK 6.8 ± 0.12 ng/mg protein, p <0.001; TVBF-7: untreated 0.12 ± 0.025 vs. FSK 0.13 ± 0.009 ng/mg protein, p > 0.999; Supplemental Figure S3).

3.5. Mineralocorticoid Pathway Stimulation Reveals Dysregulated and Distinct Patterns among the Cell Lines

Next, we have assessed the response of candidates and products of the mineralocorti- coid pathway. Hence, stimulation with KCl, AngII and ACTH induced CYP11B2 expression (Figure 4A-C) in NCI-H295R cells but not in TVBF-7 cells while the enzyme expression was undetected in MUC-1. ATR1 gene expression was induced only in NCI-H295R cells upon AngII stimulation (Figure 4D) but not in the other cell lines and not upon KCl or ACTH stimulation (data not shown). The induction towards mineralocorticoid phenotype, though, was not depicted in the hormonal secretion profile within 24 h, since an aldosterone increase was detected for none of the cell lines and 18-OH corticosterone production in NCI-H295R cells upon FSK stimulation only (data are provided in the Supplementary Materials).

AngII stimulation of cells for 5 min (300 s) did furthermore not result in statistically significant changes in reversal potential (Supplemental Figure S4), but regarding membrane currents, various cells responding at different time points and/or consecutive responses of the same cell were detected (change in current density >2 fold of baseline; Figure 4E). Having isolated the responders, we could furthermore observe NCI-H295R (n = 3/20) increasing their currents at 80 mV after the first 50 s of stimulation and again, with less

intensity in later time points, while MUC-1 (3/16) and TVBF-7 (2/17) responded at 100 s or later (Figure 4F).

3.6. Responsive, Non-Producing and Autonomous Secretion Detected Regarding the Glucocorticoid Production

Next, the glucocorticoid production, respective steroidogenic gene and receptor stimu- lation upon ACTH and FSK stimulation was investigated. While the cells did not show any clear response in terms of reversal potential in the conditions measured (data are in the Supplementary Materials), ACTH stimulation revealed some electrophysiological respon- ders with excessive change in the currents at 80 mV, surprisingly mostly deriving from the MUC-1 (n = 3/8) and TVBF-7 (n = 4/18) and not from NCI-H295R cells (n =0/15) cells. However, a cell-to-cell analysis did not reveal a specific response pattern (Figure 4G,H). In contrast, upon FSK stimulation, NCI-H295R (n = 2/16) and MUC-1 (n =4/15) responded with increases in the current density reaching to a maximum 200 s upon the beginning of FSK stimulation (Figure 4I,J).

Figure 4. Activation of the mineralocorticoid pathway upon different stimulations. CYP11B2 gene

A

B

CYP11B2

CYP11B2

500


NCI-H295

250

NCI-H295

CYP11B2/GAPDH

gene expression

400

MUC-1

CYP11B2/GAPDH

gene expression

200

MUC-1

300

TVBF-7


x

150

TVBF-7

150

100

T

100

50-

50-

0

0

20

100 250

0

20 100 250

0

20 100 250

0

0

5

10

25

0

10 25

0

5

10

25

nM Angli

nM ACTH

C

D

CYP11B2

AT1R

800

200

NCI-H295


NCI-H295

MUC-1

CYP11B2/GAPDH

gene expression

600

MUC-1

AT1R/GAPDH

gene expression

150

TVBF-7

400

TVBF-7

200

100

150

100

50-

50

0

0

0

0

6

10

20

0

5

10

20

0

5

10

20

0

20 100 250

0

20 100 250

0

20

100 250

mM KCI

nM Angll

E

Angll at 80 mV (all cells)

F

Ang, 80mV

12-

50

current fold change (normalized for time point 0)

current fold change (normalized for time point 0)

30

- NCI-H295R

10-

MUC-1

TVBF-7

8-

10-

6-

4.

5-

2.

0

0

0

50

100

150

200

250

300

0

50

100

150

200

250

300

time (sec)

time (sec)

G

ACTH at 80 mV (all cells)

H

AC TH, 80mV

50-

21

current fold change (normalized for time point 0)

30

current fold change (normalized for time point 0)

MUC-1

20

TVBF-7

X

8-

10-

6-

4-

5-

2

0

0

0

50

100

150

200

250

300

0

50

100

150

200

250

300

time (sec)

time (sec)

FSK at 80 mV (all cells)

FSK, 80mV

5-

5-

NCI-H295R

current fold change (normalized for time point 0)

current fold change (normalized for time point 0)

MUC-1

*

1

W

3

2

2

1.

1

0

0

50

100

150

200

250

300

0

0

50

100

150

200

250

300

time (sec)

time (sec)

expression upon AngII (A), ACTH (B) and KCl (C) per cell line as well as AT1R expression upon AngII stimulation (D). Scatter dot representation of the fold change of the current densities versus baseline at 80 mV for 5 min (300 s) of AngII stimulation cumulatively for all cells per time point aiming to divide the responders (>2 fold change, marked with red) from the non-responders (E). Representation of the fold change of the currents of the responder cells in conjunction with time for the AngII stimulation (F). Scatter dot representation of the fold change of the current densities versus baseline at 80 mV for 5 min (300 s) of ACTH stimulation cumulatively for all cells per time point aiming to divide the responders (>2 fold change, marked with red) from the non-responders (G). Representation of the fold change of the currents of the responder cells in conjunction with time for the ACTH stimulation (H). Scatter dot representation of the fold change of the current densities versus baseline at 80 mV for 5 min (300 s) of FSK stimulation cumulatively for all cells per time point aiming to divide the responders (>2 fold change, marked with red) from the non-responders (I). Representation of the fold change of the currents of the responder cells in conjunction with time for the FSK stimulation (J). Stars represent significance vs. untreated cells ( ** p < 0.01; *** p < 0.001).

NCI-H295R responded upon ACTH and NCI-H295R and MUC-1 upon FSK with increase of the MC2R and CYP21A2 gene expression (Figure 5B,D,F) but only in NCI-H295R this was also accompanied by an increase of CYP11B1 gene expression as well as elevation in cortisol and 21-deoxycortisol secretion (Figure 5A,C,E,H). As already indicated before, TVBF-7 demonstrated basally a profile represented by high basal glucocorticoid expression with upregulated MC2R and CYP11B1 expression levels which was again confirmed by extraordinary high basal cortisol (Figure 5E) and 21-Deoxycortisol (Figure 5H) levels. However, in contrast to NCI-H295R they remained unresponsive upon ACTH and FSK stimulation, which corresponds to an autonomous glucocorticoid secretion profile.

3.7. MUC-1 Induce Androgen Production and Upregulation of the AR

While NCI-H295R demonstrated basally very high androstendione (Figure 6C) and TVBF-7 DHEAS levels (Figure 6B), only MUC-1 represented a remarkable stimulability and downstream response in the context of androgens (Figure 6A,C,D) which was then by far exceeding DHEA and also Testosterone levels of NCI-H295R and TVBF-7 (Figure 6A,D). Apart from the above-mentioned distinct and around 200 s also specifically patterned electrophysiological response of MUC-1 cells upon FSK stimulation, this was accompanied by highly significant increases in expression of HSD17B4 (Figure 6E), androgen receptor gene (AR, Figure 6F), gonadotropin releasing hormone receptor gene (GNRH, Figure 6J), SF-1 (Figure 6I) and tendencies towards ER1-upregulation (Figure 6G) for MUC-1. Inter- estingly, in the context of nuclear receptors also the glucocorticoid receptor was under these conditions strongly upregulated for MUC-1 only. In contrast, HSD17B4 remained unchanged, SF-1 was upregulated and AR and ER1 were conversely down-regulated for NCI-H295R. For TVBF-7, again, all levels remained unchanged.

Figure 5. Activation of the glucocorticoid pathway upon different stimulations. CYP11B1 and MC2R gene expression upon ACTH (A,B) and FSK (C,D), respectively, per cell line. Hormonal production of cortisol (E) in ng per mg of total protein for unstimulated cells in comparison with KCI, AngII, FSK and two different ACTH stimulations. Activation of intermediate steroidogenesis steps upon different stimulations. CYP21A2 gene expression upon FSK stimulation (F) per cell line. Hormonal production of 11- (G) and 21-deoxycortisol (H) in ng per mg of total protein for unstimulated cells in comparison with KCI, AngII, FSK and two different ACTH stimulations. Stars represent significance vs. untreated cells (* p < 0.05; ** p < 0.01; *** p < 0.001).

A

CYP11B1

B

MC2R

800

NCI-H295

300

NCI-H295

CYP11B1/GAPDH

gene expression

600

MUC-1

MUC-1

400

TVBF-7

MC2R/GAPDH gene expression

TVBF-7

150

200

100

100

50

0

0

5

10

25

0

5

10

25

0

5

10

25

0

nM ACTH

0

10

25

0

5

10

25

0

5

10

25

nM ACTH

C

CYP11B1

D

MC2R

30000

NCI-H295

1500

NCI-H295

CYP11B1/GAPDH

gene expression

20000

MUC-1

MC2R/GAPDH

gene expression

1000

MUC-1

TVBF-7

TVBF-7

150

500

150

100

100

50

50

0

0

0

5

10

25

0

5

10

25

0

5

10

25

0

U

10

25

0

5

10

25

0

5

10

25

μM FSK

μM FSK

E

NCI-H295

F

Cortisol

MUC-1

CYP21A2

TVBF-7

1500

4000

3000

NCI-H295

1000

Cortisol (ng/mg protein)

CYP21A2/GAPDH gene expression

2000

MUC-1

500

1000

TVBF-7

150

75

100

35

50

Control

KCI

Angli

FSK ACTH (I)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

0

ACTH (1)

ACTH (h)

0

U1

10

25

0

5

10

25

0

5

10

25

μM FSK

G

H

11-Deoxycortisol

21-Deoxycortisol

NCI-H295

400

NCI-H295

TVBF-7

MUC-1

150

MUC-1

11-Deoxycortis ol (ng/mg protein)

TVBF-7

200

21-Deoxycortisol (ng/mg protein)

75

5

50

18

15

2

0.2

12

1.5

0.1

1.0

0.5

0.0

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

0.0

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (I) ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Figure 6. Androgenic pathway activation. Hormonal production of DHEA (A), DHEAS (B), An- drostenedione (C) and Testosterone (D) in ng per mg of total protein for unstimulated cells in comparison with KCI, AngII, FSK and two different ACTH stimulations. HSD17B4 gene expression upon FSK (E) per cell line. Differential response of different hormonal receptor expression per cell 36 line upon FSK stimulation: expression levels of AR (F), ER1 (G), GCR (H), SF-1 (I) and GNRHR (J) Stars represent significance vs. untreated cells (* p<0.05; ** p < 0.01; *** p < 0.001).

A

B

DHEA

NCI-H295

DHEAS

NCI-H295

MUC-1

26

MUC-1

TVBF-7

6000

TVBF-7

DHEA (ng/mg protein)

8~

DHEAS (ng/mg protein)

4000

2000

200

150

1

100

0

H

50

ACTH (h)

Control

FSK

ACTH (h)

Control

KCI

Angli

FSK

0

Control

KCI

Angli

FSK

ACTH (I)

KCI

Angli

ACTH (I)

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

CI

Angli

FSK

ACTH (1)

ACTH (h)

C

Androstenedione

NCI-H295

D

Testosterone

MUC-1

NCI-H295

250

TVBF-7

0.20

MUC-1

150

TVBF-7

Androstenedione (ng/mg protein)

50

Testosterone (ng/mg protein)

0.15

5

4

0.10-

3

1.0

0.05

0.5

0.0

0.00

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (I)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (1)

ACTH (h)

Control

KCI

Angli

FSK ACTH (1)

ACTH (h)

Control

CI

Angli

FSK

ACTH (I)

ACTH (h)

Control

KCI

Angli

FSK

ACTH (I)

ACTH (h)

E

F

HSD17B4

NCI-H295

AR

NCI-H295

150

MUC-1

200

MUC-1

HSD17B4/GAPDH gene expression

TVBF-7

TVBF-7

100

AR/GAPDH

gene expression

150

100

50-

50.

0

0

5

10

25

0

5

0

10

25

0

5

10

25

0

U

10

25

0

5

10

25

0

5

10

25

μM FSK

HM FSK

G

H

ER1

NCI-H295

GCR

NCI-H295

600-

MUC-1

200-

MUC-1

TVBF-7

TVBF-7

ER1/GAPDH gene expression

400

GCR/GAPDH

gene expression

150

200

150

100

100

50-

50-

0

0

0

U

10

25

0

5

10

25

0

5

10

25

0

5

10

25

0

5

10

25

0

5

10

25

HM FSK

UM FSK

SF-1

GNRHR

NCI-H295

NCI-H295

8000-


MUC-1

400

MUC-1

TVBF-7


TVBF-7

SF-1/GAPDH gene expression

GNRHR/GAPDH

gene expression

6000-

300

4000

200

250

200

150

100

100

50

0

0

0

5

10

25

0

5

10

25

0

U

10

25

0

5

10

25

0

5

10

25

0

5

10

25

UM FSK

HM FSK

4. Discussion

In our study, we have established a novel cell line, namely TVBF-7, which we character- ized in direct comparison to the gold standard NCI-H295R and the highly emerging MUC-1 cells. In a first step, we investigated a large panel of known driver genes. Our data demon- strate heterogenous genotypes. For TP53, e.g., we confirmed the known mutations for NCI-H295R and MUC-1 [13] and newly characterized TVBF-7 to be a wildtype. Our studies confirmed also the previously described activating CTNNB1 mutation for NCI-H295R [30]. For TVBF-7, we identified a non-sense APC mutation, while MUC-1 represented wildtype in both genetic loci. Overall, our data suggest a cell panel that includes important tools for preclinical mechanistic and therapeutic studies in the future. The Wnt/ -catenin pathway is a target of many novel treatments suggested for ACC on a preclinical level [31]. Even though a focus of great research, detailed mechanistic and therapeutic studies in a panel of human cell lines naturally reflecting different genotypes are still lacking. Similar considera- tions can be taken into account for TP53 or other potential candidates. Overall, the initial genomic data revealed that the three cell lines cover important aspects of the mutational heterogeneity met in patients with ACC [32].

In two therapeutic studies, NCI-H295R, MUC-1 and TVBF-7 (as primary culture ACC115m) have been already implemented. Interestingly, independent from genetic status both metastatic models often shared rather less therapeutic responsiveness compared to NCI-H295R [27,29], a therapeutic phenotype which is already well known from further studies implementing MUC-1 [20-23,25,26,28]. It should be mentioned in this context that, in contrast to the primary derived chemo-naive NCI-H295R, both metastatic models were obtained from EDP-M treated patients [15,27]. However, patients with advanced, metastatic and pre-treated disease are the ones urgently requiring ongoing multi-chemotherapeutic treatments, while early stages with a low risk of recurrence undergo in first-instance surgery [1]. Thus, in the best case, preclinical platforms provide such patient tumor hetero- geneities as various stages, therapeutic responsiveness, genotypes, hormonal phenotypes and gender, among others.

Consequently, we went on with the investigation of other important aspects such as steroidogenic signaling, secretion and, electrophysiological responsiveness using an innovative and unbiased automated high-through-put patch-clamp system. Alongside maintenance in characteristics of steroid producing cells such as the most hyperpolarized reversal potential and sensitivity to increased potassium concentrations [33], our studies revealed for TVBF-7 extremely high expression of MC2R, CYP11B1, cortisol and DHEAS secretion at baseline. The cells were furthermore unresponsive to a selection of known physiological stimuli. A profile comparable to autonomous cortisol secretion could be found [34]. Hypercortisolism in such patients is often mild, and most patients lack typical clinical features of overt Cushing’s syndrome. Likewise, TVBF-7 was derived from a patient without obvious clinical signs of hormone excess [27]. However, a correlation of APC inactivating mutations and elevated cortisol levels has been already previously reported for patients with adrenal tumors [30,35].

In contrast, for MUC-1 which was also derived from a male patient tumor originally diagnosed as hormonally diffuse with no clinical signs of steroid excess [15,36] at baseline no parameter of steroid excess was detectable. Fittingly, MUC-1 demonstrated a profile with comparably low steroidogenic activity and secretion. However, upon FSK stimulation MUC-1 demonstrated the capacity to induce impressively strong main steroidogenic genes such as CYP11A1, HSD3B2, CYP17A1 and also the upregulation of genes involved in glucocorticoid metabolism such as MC2R and CYP21A2. Electrophysiologically, MUC-1 responded furthermore to FSK showing a clear pattern, which is in alignment with the respective outcome for gene expression and hormonal levels. Interestingly, in previous studies with MUC-1 xenografts and first passages after explantation from a female host, MUC-1 demonstrated increased cortisol levels [15]. Hormonal patterns in ACC might be, thus, conditional and inducible by a variety of local factors. Tumors, clinically considered as not functional, might be thereby still able to lead to profound hormonal changes within

the tumor-microenviroment, even if not leading to whole body excess. This finding could be of high clinical relevance in the context of therapeutic treatments and responsiveness. Moreover, also gender aspects could be involved. Data which are scarce in the literature concerning the effects of gender on adrenal tumors [37]. However, it is known that cortisol- secreting adrenal tumors are more often diagnosed in female patients [38]. Of note, SF-1 is a key regulator of human sex determination [39] and its activation might lead to different downstream effects in tissues of male and female origin.

Androgens, estrogens, and progestins are known as sex steroids. For prostate and breast cancer mechanisms of sex-steroid hormone-regulated DNA damage repair is already well known [40]. Recently, our group reported marked differences in DNA damage repair leading to a rather drug-resistant phenotype for MUC-1 vs. NCI-H295R [20]. Together with the data obtained from the current study, it is prudent to speculate about a potential involvement of AR and its ligands in the drug resistant phenotype of MUC-1. Of note, the AR was markedly differently regulated in NCI-H295R (down), MUC-1 (up) and TVBF-7 (unchanged). The main AR ligand testosterone was within the observed time-frame de- tectable for MUC-1 exclusively, as well as an extraordinary increase in DHEA. DHEA is known to augment AR levels, and, furthermore, is able to compete with dihydrotestos- terone indicating intrinsic androgenic activity that is potentially independent of metabolic conversion to other androgens and can affect gene function through the AR [41]. In contrast, further studies demonstrate that androgenicity of DHEAS, which is instead rather elevated in NCI-H295R and even more in TVBF-7, is negligible [42]. However, also other nuclear receptors (ER and GCR) underly differential regulation in the various models. In this context, it is relevant to mention that previous studies demonstrated at baseline higher HSP90-abundance in MUC-1 compared to NCI-H295R, which was again correlated with less therapeutic responsiveness towards therapeutic HSP90 inhibition for MUC-1 [22]. Overall, further studies will be required, to elucidate the underlying mechanisms of action under varying conditions in more detail.

Interestingly, also HSD17B4 a gene with a dual role in steroidogenesis and fatty acid oxidation is specifically upregulated for MUC-1 under FSK-conditions. HSD17B4 expression is known to increase in castration resistant prostate cancer, leading to metabolic re-programming resulting in AR-stimulation and poor prognosis [43]. Furthermore, a recent study reported a main role in the context of mitochondrial and peroxisomal fatty acid oxidation and that the latter one can serve as a compensatory mechanism in case of mitochondrial defects or overload [44]. This is in the specific context of great interest, as we have recently reported significant differences in cholesterol storage between mi- totane sensitive (NCI-H295R) and resistant (MUC-1) adrenocortical cells [21]. Furthermore, upcoming data indicate strong differences in storage of cholesterylester-containing lipid droplets (NCI-H295R) vs. triacyglycerol-LDs (MUC-1) and the appropriate differential lipid-metabolism in this context [45].

In sum, the aim of the current study was to characterize NCI-H295R, MUC-1 and TVBF-7 regarding their underlying heterogenic geno- and phenotypes, to thereby reveal varying steroidogenic and down-stream signalling/ secretion capacities and highlight their significant, but not interchangeable values for new therapeutic and mechanistic studies by representing various clusters known for ACC patients.

Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/cells11091439/s1, Figure S1: AT1R gene expression upon increasing KCl (A) and ACTH stimulation (B) and MC2R gene expression upon increasing KCl (C) and AngII stimulation (D). Stars represent significance vs. untreated cells (*, p < 0.05; ** , p < 0.01; *** , p < 0.001); Figure S2: Hormonal production of Aldosterone (A), 18-oxo-cortisol (B), 18-OH-cortisol (C) and 11-deoxycorticosterone (D) in ng per mg of total protein for unstimulated cells in comparison with KCl, AngII, FSK and two different ACTH stimulations. Stars represent significance vs. untreated cells (*, p< 0.05; ** , p < 0.01; *** , p < 0.001); Figure S3: Hormonal production of 18-OH-corticosterone (A), Pregnenolone (B), Progesterone (C) and cortison (D) in ng per mg of total protein for unstimulated cells in comparison with KCl, AngII, FSK and two different ACTH stimulations. Stars represent

significance vs. untreated cells (*, p < 0.05; ** , p < 0.01; *** , p < 0.001); Figure S4: Reversal potential changes over time upon stimulation with AngII (A), ACTH (B) and Forskolin (FSK) (C) for all three cell lines. Table S1: STR Profile at the University of Brescia; Table S2: Primers used in the Quantitative Real-Time PCR; Table S3: Single-nucleotide polymorphism (SNP) findings for certain ACC and general cancer genes filtered for mutations located in coding (exonic) regions, with exclusion of the synonymous (silence) mutations.

Author Contributions: Conceptualization, S.S., A.B. and C.H .; Data curation, C.B., D.P. and A.A .; Formal analysis, S.S., C.B., D.P., A.A. and C.H .; Funding acquisition, S.R.B., A.B. and C.H .; Inves- tigation, S.S., C.B., A.A., M.P., D.C., G.A.M.T., A.B. and C.H .; Methodology, S.S., C.B., D.P., A.A., M.P., D.C., G.A.M.T., A.B. and C.H .; Resources, G.A.M.T., S.R.B. and A.B .; Supervision, S.S., D.P., A.B. and C.H .; Validation, S.S., C.B. and C.H .; Writing-original draft, S.S., C.B., A.A., A.B. and C.H .; Writing-review and editing, S.S., D.P., A.A., M.P., D.C., G.A.M.T., S.R.B., A.B. and C.H. All authors have read and agreed to the published version of the manuscript.

Funding: This work was supported by the Uniscientia Foundation (keyword tumor model) to C.H .; by the Heuberg Foundation (working title: “Novel multidimensional models for the adrenal gland and adrenal tumors”) to C.H. and S.B .; by the Italian Association for Cancer Research to A.B. (funding code: IG23009).

Institutional Review Board Statement: The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of Brescia (protocol NP1924; 27 January 2015).

Informed Consent Statement: Informed consent was obtained from all patients involved in the study.

Data Availability Statement: All data that were needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. The cell lines MUC-1 and TVBF-7 can be provided by C. Hantel and S. Sigala/ A. Berruti, respectively, pending scientific review and completed MTAs.

Acknowledgments: We thank Sereina Gut for her excellent technical support. The authors thank Sophion Bioscience for their support with automated patch clamp.

Conflicts of Interest: A.B. received research funds for ACC studies from Janssen-Cilag and Sanofi. All other authors declare that there is no potential conflict of interest.

References

1. Fassnacht, M .; Dekkers, O.M .; Else, T .; Baudin, E .; Berruti, A .; de Krijger, R .; Haak, H.R .; Mihai, R .; Assie, G .; Terzolo, M. European Society of Endocrinology Clinical Practice Guidelines on the management of adrenocortical carcinoma in adults, in collaboration with the European Network for the Study of Adrenal Tumors. Eur. J. Endocrinol. 2018, 179, G1-G46. [CrossRef]

2. Schweitzer, S .; Kunz, M .; Kurlbaum, M .; Vey, J .; Kendl, S .; Deutschbein, T .; Hahner, S .; Fassnacht, M .; Dandekar, T .; Kroiss, M. Plasma steroid metabolome profiling for the diagnosis of adrenocortical carcinoma. Eur. J. Endocrinol. 2019, 180, 117-125. [CrossRef] [PubMed]

3. Arlt, W .; Biehl, M .; Taylor, A.E .; Hahner, S .; Libe, R .; Hughes, B.A .; Schneider, P .; Smith, D.J .; Stiekema, H .; Krone, N .; et al. Urine steroid metabolomics as a biomarker tool for detecting malignancy in adrenal tumors. J. Clin. Endocrinol. Metab. 2011, 96, 3775-3784. [CrossRef] [PubMed]

4. Bancos, I .; Taylor, A.E .; Chortis, V .; Sitch, A.J .; Jenkinson, C .; Davidge-Pitts, C.J .; Lang, K .; Tsagarakis, S .; Macech, M .; Riester, A .; et al. Urine steroid metabolomics for the differential diagnosis of adrenal incidentalomas in the EURINE-ACT study: A prospective test validation study. Lancet Diabetes Endocrinol. 2020, 8, 773-781. [CrossRef]

5. Chortis, V .; Bancos, I .; Nijman, T .; Gilligan, L.C .; Taylor, A.E .; Ronchi, C.L .; O’reilly, M.W .; Schreiner, J .; Asia, M .; Riester, A .; et al. Urine Steroid Metabolomics as a Novel Tool for Detection of Recurrent Adrenocortical Carcinoma. J. Clin. Endocrinol. Metab. 2020, 105, e307-e318. [CrossRef] [PubMed]

6. Lin, C.W .; Chang, Y.H .; Pu, H.F. Mitotane exhibits dual effects on steroidogenic enzymes gene transcription under basal and cAMP-stimulating microenvironments in NCI-H295 cells. Toxicology 2012, 298, 14-23. [CrossRef]

7. Lehmann, T.P .; Wrzesinski, T .; Jagodzinski, P.P. The effect of mitotane on viability, steroidogenesis and gene expression in NCIH295R adrenocortical cells. Mol. Med. Rep. 2013, 7, 893-900. [CrossRef] [PubMed]

8. Suzuki, S .; Minamidate, T .; Shiga, A .; Ruike, Y .; Ishiwata, K .; Naito, K .; Ishida, A .; Deguchi, H .; Fujimoto, M .; Koide, H .; et al. Steroid metabolites for diagnosing and predicting clinicopathological features in cortisol-producing adrenocortical carcinoma. BMC Endocr. Disord. 2020, 20, 173. [CrossRef] [PubMed]

9. Abiven, G .; Coste, J .; Groussin, L .; Anract, P .; Tissier, F .; Legmann, P .; Dousset, B .; Bertagna, X .; Bertherat, J. Clinical and biological features in the prognosis of adrenocortical cancer: Poor outcome of cortisol-secreting tumors in a series of 202 consecutive patients. J. Clin. Endocrinol. Metab. 2006, 91, 2650-2655. [CrossRef] [PubMed]

10. Kendrick, M.L .; Curlee, K .; Lloyd, R .; Farley, D.R .; Grant, C.S .; Thompson, G.B .; Rowland, C .; Young, W.F .; van Heerden, J. Aldosterone-secreting adrenocortical carcinomas are associated with unique operative risks and outcomes. Surgery 2002, 132, 1008-1011; discussion 1012. [CrossRef] [PubMed]

11. Salvi, S .; Bonafe, M .; Bravaccini, S. Androgen receptor in breast cancer: A wolf in sheep’s clothing? A lesson from prostate cancer. Semin. Cancer Biol. 2020, 60, 132-137. [CrossRef]

12. Snaterse, G .; Visser, J.A .; Arlt, W .; Hofland, J. Circulating steroid hormone variations throughout different stages of prostate cancer. Endocr. Relat. Cancer 2017, 24, R403-R420. [CrossRef]

13. Pinto, E.M .; Kiseljak-Vassiliades, K .; Hantel, C. Contemporary preclinical human models of adrenocortical carcinoma. Curr. Opin. Endocr. Metab. Res. 2019, 8, 139-144. [CrossRef] [PubMed]

14. Gazdar, A.F .; Oie, H.K .; Shackleton, C.H .; Chen, T.R .; Triche, T.J .; Myers, C.E .; Chrousos, G.P .; Brennan, M.F .; Stein, C.A .; La Rocca, R.V. Establishment and characterization of a human adrenocortical carcinoma cell line that expresses multiple pathways of steroid biosynthesis. Cancer Res. 1990, 50, 5488-5496.

15. Hantel, C .; Shapiro, I .; Poli, G .; Chiapponi, C .; Bidlingmaier, M .; Reincke, M .; Luconi, M .; Jung, S .; Beuschlein, F. Targeting heterogeneity of adrenocortical carcinoma: Evaluation and extension of preclinical tumor models to improve clinical translation. Oncotarget 2016, 7, 79292-79304. [CrossRef] [PubMed]

16. Kiseljak-Vassiliades, K .; Zhang, Y .; Bagby, S.M .; Kar, A .; Pozdeyev, N .; Xu, M .; Gowan, K .; Sharma, V .; Raeburn, C.D .; Albuja-Cruz, M .; et al. Development of new preclinical models to advance adrenocortical carcinoma research. Endocr. Relat. Cancer 2018, 25, 437-451. [CrossRef] [PubMed]

17. Landwehr, L.S .; Schreiner, J .; Appenzeller, S .; Kircher, S .; Herterich, S .; Sbiera, S .; Fassnacht, M .; Kroiss, M .; Weigand, I. A novel patient-derived cell line of adrenocortical carcinoma shows a pathogenic role of germline MUTYH mutation and high tumour mutational burden. Eur. J. Endocrinol. 2021, 184, 823-835. [CrossRef]

18. Vicennati, V .; Repaci, A .; di Dalmazi, G .; Rinaldi, E .; Golfieri, R .; Giampalma, E .; Minni, F .; Marrano, N .; Santini, D .; Pasquali, R. Combined aldosterone and cortisol secretion by adrenal incidentaloma. Int. J. Surg. Pathol. 2012, 20, 316-319. [CrossRef] [PubMed]

19. Zheng, S .; Cherniack, A.D .; Dewal, N .; Moffitt, R.A .; Danilova, L .; Murray, B.A .; Lerario, A.M .; Else, T .; Knijnenburg, T.A .; Ciriello, G .; et al. Comprehensive Pan-Genomic Characterization of Adrenocortical Carcinoma. Cancer Cell 2016, 29, 723-736. [CrossRef]

20. Bothou, C .; Sharma, A .; Oo, A .; Kim, B .; Perge, P .; Igaz, P .; Ronchi, C.L .; Shapiro, I .; Hantel, C. Novel Insights into the Molecular Regulation of Ribonucleotide Reductase in Adrenocortical Carcinoma Treatment. Cancers 2021, 13, 4200. [CrossRef] [PubMed]

21. Warde, K.M .; Schoenmakers, E .; Ribes Martinez, E .; Lim, Y.J .; Leonard, M .; Lawless, S.J .; O’Shea, P .; Chatterjee, V.K .; Gurnell, M .; Hantel, C .; et al. Liver X receptor inhibition potentiates mitotane-induced adrenotoxicity in ACC. Endocr. Relat. Cancer 2020, 27, 361-373. [CrossRef]

22. Siebert, C .; Ciato, D .; Murakami, M .; Frei-Stuber, L .; Perez-Rivas, L.G .; Monteserin-Garcia, J.L .; Nölting, S .; Maurer, J .; Feuchtinger, A .; Walch, A.K .; et al. Heat Shock Protein 90 as a Prognostic Marker and Therapeutic Target for Adrenocorti- cal Carcinoma. Front. Endocrinol. 2019, 10, 487. [CrossRef] [PubMed]

23. Hasanovic, A .; Ruggiero, C .; Jung, S .; Rapa, I .; Signetti, L .; Ben Hadj, M .; Terzolo, M .; Beuschlein, F .; Volante, M .; Hantel, C .; et al. Targeting the multidrug transporter Patched potentiates chemotherapy efficiency on adrenocortical carcinoma in vitro and in vivo. Int. J. Cancer 2018, 143, 199-211. [CrossRef] [PubMed]

24. Fragni, M .; Palma Lopez, L.P .; Rossini, E .; Abate, A .; Cosentini, D .; Salvi, V .; Vezzoli, S .; Poliani, P.L .; Bosisio, D .; Hantel, C .; et al. In Vitro cytotoxicity of cabazitaxel in adrenocortical carcinoma cell lines and human adrenocortical carcinoma primary cell cultures. Mol. Cell Endocrinol. 2019, 498, 110585. [CrossRef]

25. Liang, R .; Weigand, I .; Lippert, J .; Kircher, S .; Altieri, B .; Steinhauer, S .; Hantel, C .; Rost, S .; Rosenwald, A .; Kroiss, M .; et al. Targeted Gene Expression Profile Reveals CDK4 as Therapeutic Target for Selected Patients With Adrenocortical Carcinoma. Front. Endocrinol. 2020, 11, 219. [CrossRef]

26. Abate, A .; Rossini, E .; Bonini, S.A .; Fragni, M .; Cosentini, D .; Tiberio, G.A.M .; Benetti, D .; Hantel, C .; Laganà, M .; Grisanti, S .; et al. Cytotoxic Effect of Trabectedin in Human Adrenocortical Carcinoma Cell Lines and Primary Cells. Cancers 2020, 12, 928. [CrossRef]

27. Rossini, E .; Tamburello, M .; Abate, A .; Beretta, S .; Fragni, M .; Cominelli, M .; Cosentini, D .; Hantel, C .; Bono, F .; Grisanti, S .; et al. Cytotoxic Effect of Progesterone, Tamoxifen and Their Combination in Experimental Cell Models of Human Adrenocortical Cancer. Front. Endocrinol. 2021, 12, 669426. [CrossRef] [PubMed]

28. Cantini, G .; Fei, L .; Canu, L .; Lazzeri, E .; Sottili, M .; Francalanci, M .; Angelotti, M.L .; De Filpo, G .; Ercolino, T .; Gelmini, S .; et al. Stimulated Expression of CXCL12 in Adrenocortical Carcinoma by the PPARgamma Ligand Rosiglitazone Impairs Cancer Progression. J. Pers. Med. 2021, 11, 1097. [CrossRef] [PubMed]

29. Abate, A .; Rossini, E .; Tamburello, M .; Lagana, M .; Cosentini, D .; Grisanti, S .; Fiorentini, C .; Tiberio, G.A .; Scatolini, M .; Grosso, E .; et al. Ribociclib Cytotoxicity Alone or Combined with Progesterone and/or Mitotane in In Vitro Adrenocortical Carcinoma Cells. Endocrinology 2022, 163, bqab248. [CrossRef]

30. Hosogi, H .; Nagayama, S .; Kanamoto, N .; Yoshizawa, A .; Suzuki, T .; Nakao, K .; Sakai, Y. Biallelic APC inactivation was responsible for functional adrenocortical adenoma in familial adenomatous polyposis with novel germline mutation of the APC gene: Report of a case. Jpn. J. Clin. Oncol. 2009, 39, 837-846. [CrossRef] [PubMed]

31. Altieri, B .; Ronchi, C.L .; Kroiss, M .; Fassnacht, M. Next-generation therapies for adrenocortical carcinoma. Best Pract. Res. Clin. Endocrinol. Metab. 2020, 34, 101434. [CrossRef] [PubMed]

32. Gara, S.K .; Lack, J .; Zhang, L .; Harris, E .; Cam, M .; Kebebew, E. Metastatic adrenocortical carcinoma displays higher mutation rate and tumor heterogeneity than primary tumors. Nat. Commun. 2018, 9, 4172. [CrossRef]

33. Dierks, A .; Lichtenauer, U.D .; Sackmann, S .; Spyroglou, A .; Shapiro, I .; Geyer, M .; Manonopoulou, J .; Reincke, M .; Hantel, C .; Beuschlein, F. Identification of adrenal genes regulated in a potassium-dependent manner. J. Mol. Endocrinol. 2010, 45, 193-206. [CrossRef]

34. Ueland, G.A .; Grinde, T .; Methlie, P .; Kelp, O .; Lovas, K .; Husebye, E.S. Diagnostic testing of autonomous cortisol secretion in adrenal incidentalomas. Endocr. Connect. 2020, 9, 963-970. [CrossRef] [PubMed]

35. Gaujoux, S .; Pinson, S .; Gimenez-Roqueplo, A.P .; Amar, L .; Ragazzon, B .; Launay, P .; Meatchi, T .; Libé, R .; Bertagna, X .; Audebourg, A .; et al. Inactivation of the APC gene is constant in adrenocortical tumors from patients with familial adenomatous polyposis but not frequent in sporadic adrenocortical cancers. Clin. Cancer Res. 2010, 16, 5133-5141. [CrossRef]

36. Assie, G .; Letouze, E .; Fassnacht, M .; Jouinot, A .; Luscap, W .; Barreau, O .; Omeiri, H .; Rodriguez, S .; Perlemoine, K .; René-Corail, F .; et al. Integrated genomic characterization of adrenocortical carcinoma. Nat. Genet. 2014, 46, 607-612. [CrossRef] [PubMed]

37. Audenet, F .; Mejean, A .; Chartier-Kastler, E .; Roupret, M. Adrenal tumours are more predominant in females regardless of their histological subtype: A review. World J. Urol. 2013, 31, 1037-1043. [CrossRef]

38. Alesina, P.F .; Walz, M.K. Adrenal Tumors: Are Gender Aspects Relevant? Visc. Med. 2020, 36, 15-19. [CrossRef]

39. Sekido, R .; Lovell-Badge, R. Sex determination involves synergistic action of SRY and SF1 on a specific Sox9 enhancer. Nature 2008, 453, 930-934. [CrossRef]

40. Ganguly, S .; Naik, D .; Muskara, A .; Mian, O.Y. The Nexus of Endocrine Signaling and Cancer: How Steroid Hormones Influence Genomic Stability. Endocrinology 2021, 162, bqaa177. [CrossRef] [PubMed]

41. Lu, S.F .; Mo, Q .; Hu, S .; Garippa, C .; Simon, N.G. Dehydroepiandrosterone upregulates neural androgen receptor level and transcriptional activity. J. Neurobiol. 2003, 57, 163-171. [CrossRef]

42. Mo, Q .; Lu, S.F .; Hu, S .; Simon, N.G. DHEA and DHEA sulfate differentially regulate neural androgen receptor and its transcriptional activity. Brain Res. Mol. Brain Res. 2004, 126, 165-172. [CrossRef] [PubMed]

43. Ko, H.K .; Berk, M .; Chung, Y.M .; Willard, B .; Bareja, R .; Rubin, M .; Sboner, A .; Sharifi, N. Loss of an Androgen-Inactivating and Isoform-Specific HSD17B4 Splice form Enables Emergence of Castration-Resistant Prostate Cancer. Cell Rep. 2018, 22, 809-819. [CrossRef] [PubMed]

44. Violante, S .; Achetib, N .; van Roermund, C.W.T .; Hagen, J .; Dodatko, T .; Vaz, F.M .; Waterham, H.R .; Chen, H .; Baes, M .; Yu, C .; et al. Peroxisomes can oxidize medium- and long-chain fatty acids through a pathway involving ABCD3 and HSD17B4. FASEB J. 2019, 33, 4355-4364. [CrossRef]

45. Warde, K .; Lim, Y .- J .; Beuschlein, F .; Hantel, C .; Dennedy, M.C. Investigating the Role of Cholesterol and Lipid Trafficking in Mitotane Resistance in Adrenocortical Carcinoma. J. Endocr. Soc. 2021, 5 (Suppl. 1), A70. [CrossRef]