1. Introduction
Cushing’s disease (CD) is a serious endocrine disorder caused by an adrenocorticotropic hormone (ACTH)-secreting pituitary neuroendocrine tumor (PitNET) that stimulates the adrenal glands to overproduce cortisol [
1,
2,
3,
4]. The WHO renamed pituitary adenomas as PitNETs [
5]. While PitNETs have been defined as benign, implying that these tumors cause a disease that is not life threatening or harmful to health, in fact chronic exposure to excess cortisol has wide-ranging and detrimental effects on health. Hypercortisolism causes increased stroke rates, diabetes, obesity, depression, anxiety, and a three-fold increase in the risk of death from cardiovascular disease and cancer [
4,
6,
7,
8].
The first-line treatment for CD is pituitary surgery, which is followed by disease recurrence in 50% of patients during the 10-year follow-up period after surgery in the hands of an experienced surgeon [
9,
10,
11]. Studies have demonstrated that surgical failures and recurrences of CD are common, and despite multiple treatments, biochemical control is not achieved in approximately 30% of patients. This suggests that in routine clinical practice, initial and long-term disease remission is not achieved in a substantial number of CD patients [
7,
12]. Hence, medical therapy is often considered in the following situations: when surgery is contraindicated or fails to achieve remission, or when recurrence occurs after apparent surgical remission. While stereotactic radiosurgery treats incompletely resected or recurrent PitNETs, the main drawbacks include the longer time to remission (12–60 months) and the risk of hypopituitarism [
3,
13,
14]. There is an inverse relationship between disease duration and reversibility of complications associated with the disease, thus emphasizing the importance of identifying an effective medical strategy to rapidly normalize cortisol production by targeting the pituitary adenoma [
4,
7,
12]. Unfortunately, the lack of current standard of care treatments with low efficacy and tolerability makes CD a medical therapeutic challenge.
The overall goal of medical therapy for CD is to target the signaling mechanisms to lower cortisol levels in the body [
15,
16]. The drugs offered for treatment of CD vary in the mechanism of action, safety, tolerability, route of administration, and drug–drug interactions [
15,
16]. In the era of precision medicine [
17], where it is imperative to identify effective therapies early, there is an urgent need to accelerate the identification of therapies targeted to the ACTH-secreting pituitary tumor which are tailored for each individual patient. The absence of preclinical models that replicate the complexity of the PitNET microenvironment has prevented us from acquiring the knowledge to advance clinical care by implementing therapies specifically targeting the tumor, which would have a higher efficacy and tolerability for CD patients. In this instance, organoids can replicate much of the complexity of an tumor. An “organoid” is defined as a three-dimensional cell structure, grown from primary cells of dissociated pituitary tumors in Matrigel matrix, which proliferate, and differentiate in three dimensions, eventually replicating key biological properties of the tissue [
18]. While pituitary cell lines predominantly represent hormonal lineages, these cultures do not reproduce the primary pituitary tissue because of the tumor transformation and non-physiological 2D culture conditions [
19,
20,
21]. Pituitary tissue-derived organoids have been generated from mouse models [
22,
23]. While several human and rat pituitary spheroid/aggregate/tumoroid models have been reported, these cultures consist of poorly differentiated cells with high replicative potential which can affect drug response and produce data that poorly translate to the clinic [
24,
25]. In this study, we developed an organoid model derived from human PitNETs that replicated much of the cellular complexity and function of the patient’s tumor. Organoids derived from corticotroph PitNETs retained the genetic alterations of the patient’s primary tissue.
2. Materials and Methods
2.1. Generation and Culture of Human Pituitary Neuroendocrine Tumor (PitNET) Organoids
Patients with planned transsphenoidal surgery for pituitary tumors were identified in the outpatient neurosurgery clinics. Tissues were collected under the St. Joseph’s Hospital and Barrow Neurological Institute Biobank collection protocol PHXA-05TS038 and collection of outcomes data protocol PHXA-0004-72-29, with the approval of the Institutional Review Board (IRB) and patient consent. Samples were de-identified and shipped to the Zavros laboratory (University of Arizona) for processing.
Pituitary tumor tissue was collected in Serum-Free Defined Medium (SFDM) supplemented with ROCK inhibitor (Y27632, 10 µM), L-glutamine (2 mM), A83-01 (activin receptor-like kinase (Alk) 4/5/7 inhibitor, 0.5 mM), penicillin/streptavidin (1%), kanamycin (1%), amphotericin/gentamycin (0.2%), CHIR-98014 (4 mM), and thiazovivin (TZV, 2.5 mM). Tissues that contained red blood cells were incubated with Red Blood Cell (RBC) Lysis Buffer according to the manufacturer’s protocol (Thermo Fisher Scientific, San Fransisco, CA, USA). Tissues were dissected into small pieces, transferred to digestion buffer (DMEM/F12 supplemented with 0.4% collagenase 2, 0.1% hyaluronic acid, 0.03% trypsin-EDTA) and incubated for 5–10 min at 37 °C with gentle shaking. Tissue was further incubated with Accutase™ (Thermo Fisher Scientific) for 5 min at 37 °C. Enzymatically dissociated cells were pelleted and washed in DPBS supplemented with antibiotics at a 400 relative centrifugal force (RCF) for 5 min. Dissociated adenoma cells were resuspended in Matrigel™, and Matrigel™ domes containing the cells were then plated in culture dishes and overlaid with pituitary growth media (
Supplemental Table S1). The culture was maintained at 37 °C at a relative humidity of 95% and 5% CO
2. Organoid growth medium was replenished every 3–4 days and passaged after 15 days in culture.
2.2. Generation of Induced Pluripotent Stem Cells (iPSCs)
Induced pluripotent stem cell lines (iPSC lines) were generated from control individuals (no reported disease) or CD patients according to published protocols by the University of Arizona iPSC Core [
26]. All human iPSC lines were tested and found to be negative for mycoplasma contamination using the Mycoalert Mycoplasma testing kits (LT07-318, Lonza), and no karyotype abnormalities were found (KaryoStat+, Thermo).
2.3. Pituitary Organoids Generated from iPSCs
Six well culture plates were coated with 2 mL/well 0.67% Matrigel (diluted in E8 media, UA iPSC core, 151169-01) and incubated at 37 °C at a relative humidity of 95% and 5% CO
2 overnight. The iPSC lines were reprogrammed from the blood of either a healthy donor (JCAZ001) or a CD patient (iPSC
7 and iPSC
1063) at the University of Arizona iPSC Core. Passage 12 iPSCs were plated onto the coated plates and incubated at 37 °C at a relative humidity of 95% and 5% CO
2. At 70% confluency, cells were passaged to freshly coated 24 well plates at a ratio of 1:8 and grown to 85–90% confluency before beginning the directed differentiation schedule. From days 0 to 3, cells were cultured in E6 media supplemented with 1% penicillin/streptomycin, 10 μM SB431542, and 5 ng/mL BMP4. BMP4 was withdrawn from the culture at day 3. Starting on day 4, the cells were cultured in E6 media, supplemented with 10 μM SB431542, 30 ng/mL human recombinant SHH, 100 ng/mL FGF8b, 10 ng/mL FGF18, and 50 ng/mL FGF10. Fifteen days after culture, the cells were harvested in cold E6 media by pipetting and resuspended in Matrigel™ (20,000 cells/50 mL Matrigel™). Matrigel™ domes containing the cells were plated in culture dishes and overlaid with differentiation media containing E6 media which was supplemented with 10 μM Y-27632, 30 ng/mL human recombinant SHH, 100 ng/mL FGF8b, 10 ng/mL FGF18, and 50 ng/mL FGF10 (
Supplemental Table S2). Organoids were cultured for a further 15 days at 37 °C at a relative humidity of 95% and 5% CO
2.
2.4. Spectral Flow Cytometry (Cytek™ Aurora)
The multicolor flow cytometry panel was designed using the Cytek
® Full Spectrum Viewer online tool to calculate the similarity index (
Supplemental Figure S1). The organoids were harvested in cold SFDM media and centrifuged at 400× g for 5 min. Supernatant was discarded and organoids were dissociated to single cells using Accutase
® (Thermo Fisher Scientific 00-4555-56). The enzymatic reaction was stopped using prewarmed DPBS, and cells were then centrifuged at 400× g for 5 min and incubated with fluorochrome-conjugated/unconjugated primary surface or cytoplasmic antibodies (
Supplemental Figure S1) at 4 °C for 30 min. Cells were then washed with Cell Staining Buffer (BioLegend # 420-201) and incubated with secondary antibodies (
Supplemental Figure S1) at 4 °C for 30 min. Cells were fixed using Cytofix/Cytoperm™ Fixation/Permeabilization Solution (BD Biosciences # 554714) at 4 °C for 20 min, followed by washing with Fixation/Permeabilization wash buffer. Cells were labeled with fluorochrome-conjugated/unconjugated intracellular primary antibodies (
Supplemental Figure S1) at 4 °C for 30 min, then washed and incubated with secondary antibodies at 4 °C for 30 min. Cells were resuspended in cell staining buffer and fluorescence and measured using the Cytek Aurora 5 Laser Spectral Flow Cytometer. An unstained cell sample was fixed and used as a reference control. UltraComp eBeads™, Compensation Beads (Thermo Fisher Scientific # 01-2222-42) were stained with the individual antibodies and used as single stain controls for compensation and gating. Data were acquired using the Cytek™ Aurora and analyzed using Cytobank software (Beckman Coulter, Indianapolis, IN, USA).
2.5. Whole Mount Immunofluorescence
Organoids were immunostained using published protocols by our laboratory [
27,
28,
29]. Proliferation was measured by using 5-ethynyl-2′-deoxyuridine (EdU) incorporation according to the Manufacturer’s protocol (Click-IT EdU Alexa Fluor 555 Imaging Kit, Thermo Fisher Scientific C10338). Co-staining was performed by blocking fixed organoids with 2% donkey serum (Jackson Immuno Research, # 017-000-121) diluted in 0.01% PBST for 1hr at room temperature. Organoids were then incubated overnight at 4 °C with primary antibodies, followed by secondary antibodies and Hoechst (Thermo Fisher Scientific H1399, 1:1000 in 0.01% PBST) for 1 h at room temperature. Human specific primary antibodies used included: rabbit anti-ACTH (Thermo Fisher Scientific 701293, 1:250), rabbit anti-Synaptophysin (Thermo Fisher Scientific PA5-27286, 1:100), species PIT1 (Thermo Fisher Scientific PA5-98650, 1:50), rabbit anti-LH (Thermo Fisher Scientific PA5-102674, 1:100), mouse anti-FSH (Thermo Fisher Scientific MIF2709, 1:100), mouse anti-PRL (Thermo Fisher Scientific CF500720, 1:100), Alexa Flour conjugated GH (NB500-364AF647, 1:100), and mouse anti-CAM5.2 (SIGMA 452M-95, 1:250). The secondary antibodies used included Alexa Fluor 488 Donkey Anti Rabbit IgG (H+L) (Thermo Fisher Scientific A21206, 1:100) or Alexa Fluor 647 Donkey Anti Mouse IgG (H+L) (Thermo Fisher Scientific A31571, 1:100). Organoids were visualized and images were acquired by confocal microscopy using the Nikon CrestV2 Spinning Disk (Nikon, Melville, NY, USA). Fluorescence intensity and percentage of EdU positive cells of total cells, were calculated using Nikon Elements Software (Version 5.21.05, Nikon, Melville, NY, USA).
2.6. Nuclear Morphometric Analysis (NMA)
Nuclear Morphometric Analysis (NMA) using treated organoids was performed based on a published protocol that measures cell viability based on the changes in nuclear morphology of the cells, using nuclear stain Hoechst or DAPI [
30]. Images of organoid nuclei were analyzed using the ImageJ Nuclear Irregularity Index (NII) plugin for key parameters, which included cell area, radius ratio, area box, aspect, and roundness. Using the published spreadsheet template [
30], the NII of each cell was calculated with the following formula: NII = Aspect − Area Box + Radius Ratio + Roundness. The area vs. NII of vehicle-treated cells were plotted as a scatter plot using the template, and was considered as the normal cell nuclei. The same plots were generated for each condition, and the NII and area of treated cells were compared to the normal nuclei, and classified as one of the following NMA populations: Normal (N; similar area and NII), Mitotic (S; similar area, slightly higher NII), Irregular (I; similar area, high NII), Small Regular (SR; apoptotic, low area and NII), Senescent (LR; high area, low NII), Small Irregular (SI; low area, high NII), or Large Irregular (LI; high area, high NII). Cells classified as SR exhibited early stages of apoptosis, and cells classified as either I, SI, or LI exhibited significant nuclear damage. The percentage of cells in each NII classification category were calculated and plotted as a histogram using GraphPad Prism.
2.7. ELISA
Concentration of secreted ACTH in conditioned media that was collected from organoid cultures was measured using the Human ACTH ELISA Kit (Novus Biologicals, NBP2-66401), according to the manufacturer’s protocol. The enzyme–substrate reaction was measured spectrophotometrically (BioTek Gen5 Micro Plate Reader Version 3.11, Santa Clara, CA, USA) at a wavelength of 450 nm, and the ACTH concentration (pg/mL) was interpolated by a standard curve with a 4-parameter logistic regression analysis, using GraphPad Prism (Version 9.2.0, San Diego, CA, USA).
2.8. Drug Assay
Patient adenoma-derived pituitary organoids were grown in 96-well plates and treated with 147 small molecules taken from the NCI AOD9 compound library for 72 h. (
https://dtp.cancer.gov/organization/dscb/obtaining/available_plates.html (accessed on 22 August 2021)). Drugs were diluted from 10 mM DMSO stock plates into 100 M DMSO working stocks with a final concentration of 1μM. All vehicle controls were treated with 0.1% DMSO. Organoid proliferation was measured using a CellTiter 96
® AQueous One Solution Cell Proliferation Assay kit (MTS, Promega, G3582, Madison, WI, USA) according to the manufacturer’s instruction. Organoid death was calculated based on the absorbance readings at 490 nm, collected from the MTS assay relative to the vehicle controls. Drug screens were performed with biological replicates in the same screen. Drugs were selected based on their ability to target key signaling pathways as well as clinical relevance to the treatment. Drug sensitivity is represented by cell viability, and is significant at <0.5 suppressive effect of the drugs. The percent of cell viability relative to the vehicle control was calculated. Correlation coefficients across each organoid were calculated using the Pearson method to assess confidence in replication. The variance component was detected for each drug across all organoids. A random effect model was run with a single random factor for each drug, and estimated variance was calculated by rejecting the null hypothesis that variation was not present among samples. The drug responses were grouped by variance factor, into large (vc > 100), median (100 > vc > 50), and small (vc < 50). A heatmap was used to display the differential responses in cell viability for the drugs.
Drugs that clustered together and showed response within corticotrophs were investigated further based on their mode of action. Pathways (Kegg and Reactome) and gene ontology mapping were conducted for the genes that were being targeted by the drugs, in order to evaluate the key responses in cellular processes. A network was constructed in Cytoscape v 3.8.2 (San Diego, CA, USA) for the purpose of association between the drugs and genes.
2.9. Drug Dose Responses
Organoids were grown in Matrigel™ domes within 96-well round-bottom culture plates. Recombinant human SHH was removed from the pituitary organoid growth media, 24 h prior to drug treatment. Organoids were treated with either vehicle (DMSO), cabergoline (Selleckchem S5842), ketoconazole (Selleckchem S1353), roscovitine (Selleckchem S1153), GANT61 (Stemcell Technologies 73692), pasireotide (TargetMol TP2207), mifeprostone (Selleckchem S2606), etomidate (Selleckchem S1329), mitotane (Selleckchem S1732), metyropane (Selleckchem S5416), or osilodrostat (Selleckchem S7456) at concentrations of 0, 1, 10, 100, 1000, and 10,000 nM, for 72 h. The percentage of cell viability was measured using an MTS assay (Promega G3580). Absorbance was measured at 490 nm and normalized to the vehicle. Concentrations were plotted in a logarithmic scale, and a nonlinear dose response curve regression was calculated using GraphPad Prism. An IC50 value for each drug treatment was determined based on the dose response curve, using GraphPad Prism analysis software.
2.10. Calculation of Area under the Curve (AUC)
AUC (area under the curve) was determined by plotting the normalized % cell viability versus transformed concentration of the drugs, using a trapezoidal approximation for the area [31]. The formula was based on splitting the curve into trapezoids with bases equal to the % viability (V) and height equal to the interval length (difference in concentrations (C), and then summing the areas of each trapezoid:
2.11. Quantitative RT PCR (qRT-PCR)
RNA was collected from patient-derived organoid cultures using the RNeasy Mini Kit (Qiagen). cDNA was generated from the extracted RNA, and then pre-amplified using TaqMan PreAmp Master Mix (Thermo Fisher Scientific 391128). The primers used were human-specific GAPDH (Thermo Fisher Scientific, Applied Biosystems Hs02786624_g1), NR5A1 (SF1) (Thermo Fisher Scientific, Hs00610436_m1), PIT1 (Thermo Fisher Scientific, Hs00230821_m1), TPit (Thermo Fisher Scientific, Hs00193027), and POMC (Thermo Fisher Scientific, Hs01596743_m1). Each PCR reaction was performed using a final volume of 20 µL, composed of 20X TaqMan Expression Assay primers, 2X TaqMan Universal Master Mix (Applied Biosystems, TaqMan
® Gene Expression Systems), and a cDNA template. Amplification of each PCR reaction was conducted in a StepOne™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA), using the following PCR conditions: 2 min at 50 °C, 10 min at 95 °C, denaturing for 15 s at 95 °C, and annealing/extending for 1 min at 60 °C, for a total of 40 cycles. Relative fold change was calculated using the 2 − ∆∆Ct method [
32], where CT = threshold cycle. Results were analyzed as the average fold change in gene expression compared to the control, and GAPDH served as an internal control.
2.12. Whole Exome Sequencing
WES was performed by the University of Arizona Center for Applied Genetics and Genomic Medicine. Isolated DNA from patient adenoma tissue will be quantified using the Qubit quantitation system with standard curve, as per the supplier protocol (Thermo Fisher Scientific). All samples were further tested for quality using the Fragment Analyzer (Advanced Analytical), following the manufacturer-recommended protocols. Whole exome sequencing (WES) was performed by array capture and approximately 60 Mb of exome target sequence, using the SureSelectXT Human All Exon V6 enrichment (Agilent) or equivalent (which one was used). All exome library builds were quantified via qPCR and subsequently sequenced to a minimum 20X coverage, using paired-end chemistry on the Illumina NovaSeq platform. Whole exome sequencing (WES) was performed by hybridization capture of approx. 35 Mb of the exome target sequence, using the Swift Exome Hyb Panel (Swift Biosciences 83216). All exome library builds were quantified via qPCR and subsequently sequenced to a minimum 20X coverage, using paired-end chemistry on the Illumina NextSeq500 or NovaSeq platform (Illumina). DNA reads were trimmed, filtered by quality scores and aligned to the human genome (hg38) with Burrows–Wheeler Aligner with default parameters. Picard (
http://broadinstitute.github.io/picard (accessed on 22 December 2021)) was used to mark duplicates. Germline single nucleotide variants (SNV) were called using the Genome Analysis Tool Kit (GATK), using the given guidelines. Mutations were annotated using ANNOVAR for coding sequences. Variants that passed the quality filter were further investigated for similarity. Concordance between tissue and organoids was calculated using Jaccard similarity index (J
ij = M
ij/(M
i + M
j − M
ij) where M
i is the number of variants in tissues, M
j is the number of variants in organoids, and M
ij is the number of identical variants in both tissue and organoid.
2.13. Single Cell RNA Sequencing (scRNA-Seq)
Cultures were collected on day 15 of the pituitary directed differentiation schedule, and cells were dissociated into a single-cell suspension using Cell Dissociation Buffer (Thermo Fisher Scientific 13151014). Cells (15,000 cells/sample) were resuspended in the sample buffer (BD Biosciences 65000062), filtered using cell strainer (40 microns), and loaded into a BD Rhapsody cartridge (BD Biosciences 400000847) for single-cell transcriptome isolation. Based on the BD Rhapsody system whole-transcriptome analysis for single-cell whole-transcriptome analysis, microbead-captured single-cell transcriptomes were used to prepare a cDNA library. Briefly, double-stranded cDNA was first generated from the microbead-captured single-cell transcriptome in several steps, including reverse transcription, second-strand synthesis, end preparation, adapter ligation, and whole-transcriptome amplification (WTA). Then, the final cDNA library was generated from double-stranded full-length cDNA by random priming amplification using a BD Rhapsody cDNA Kit (BD Biosciences, 633773), as well as the BD Rhapsody Targeted mRNA and WTA Amplification Kit (BD Biosciences, 633801). The library was sequenced in PE150 mode (paired-end with 150-bp reads) on NovaSeq6000 System (Illumina). A total of 80,000 reads were demultiplexed, trimmed, mapped to the GRCh38 annotation, and quantified using the whole transcriptome analysis pipeline (BD Rhapsody™ WTA Analysis Pipeline v1.10 rev6, San Jose, CA, USA) on the Seven Bridges Genomics platform (
https://igor.sbgenomics.com (accessed on 4 April 2022)), prior to clustering analysis in Seurat. For QC and filtration, read counting and unique molecular identifier (UMI) counting were the principal gene expression quantification schemes used in this single-cell RNA-sequencing (scRNA-seq) analysis. The low-quality cells, empty droplets, cell doublets, or multiplets were excluded based on unique feature count (less than 200 or larger than 2500), as they may often exhibit either an aberrantly high gene count or very few genes. Additionally, the mitochondrial QC metrics were calculated, and the cells with >5% mitochondrial counts were filtered out, as the percentage of counts originating from a set of low-quality or dying cells often exhibit extensive mitochondrial contamination. After the removal of unwanted cells from the single cell dataset, the global-scaling normalization method LogNormalize was employed. This method normalizes the feature expression measurements for each cell by the total expression, multiplies this by a scale factor (10,000), and log-transforms the result. The molecules per gene per cell, based on RSEC error correction (RSEC_MolsPerCell file) matrix files from iPSC
ctrl and iPSC
CDH23 samples, were imported into Seurat v4, merged, and processed (as stated above) for UMAP reduction, cluster identification, and differential marker assessment using the FindAllMarkers function within Seurat.
2.14. Statistical Analyses
Sample size was based on assessment of power analysis using SigmaStat software. Data collected from each study from at least 4 in vitro technical replicates were analyzed by obtaining the mean ± standard error of the mean (SEM), unless otherwise stated. The significance of the results was then tested using commercially available software (GraphPad Prism, GraphPad software, San Diego, CA, USA).
4. Discussion
Our studies demonstrate the development of organoids generated from human PitNETs (hPITOs) can potentially be used to screen for the sensitivity and efficacy of responses to targeted therapies for CD patients that either fail to achieve remission or exhibit recurrence of disease after surgery. In addition, we have documented that induced pluripotent stem cells (iPSCs) generated from a CD patient expressing germline mutation CDH23 (iPSC
CDH23) reveals the disease pathogenesis under directed differentiation. Many early in vitro experiments have used pituitary cell lines, spheroids, aggregates, and/or tumoroids that do not replicate the primary PitNET microenvironment [
19,
20,
21], and lack a multicellular identity [
39,
40]. The development of PitNET tissue-generated organoids is limited to the use of transgenic mouse models as the source [
22,
23,
41]. The recent organoid cultures reported by Nys et al. [
42] have been generated from single stem cells isolated from PitNET tissue, and are claimed to be true organoids due to their clonality. However, multicellular complexity was not validated by the protein expression or hormone secretion from pituitary cell lineages in these cultures [
42]. According to the National Cancer Institute (NCI, NIH), an ‘organoid’ is defined as “a tiny, 3-dimensional mass of tissue that is made by growing stem cells (cells from which other types of cells develop) in the laboratory” [
43]. The hPITOs reported here begin from single and/or 3–4 cell clusters dissociated from the PitNET tissue that harbors the stem cells.
Supplemental Video S2 demonstrates a process of ‘budding,’ as well as lumen formation as organoids grow and differentiate. We document differentiation and function by comprehensive spectral flow cytometry, ELISA, and response to standard of care drugs. The growth of PitNET organoids reported in the current study is consistent with that of gastrointestinal tissue derived cultures that begin from cell clusters, crypts, or glands [
27,
44,
45].
Our studies report a PitNET tissue organoid culture with a multicellular identity consisting of differentiated cell lineages, stem/progenitor cells, and immune and stromal cell compartments, which replicates much of the patient’s own adenoma pathology, functionality, and complexity. We have also demonstrated that iPSCs, derived from the blood of a CD patient, can be directly differentiated into pituitary organoids that resemble similar characteristics to the tumor tissue. Many investigators have proposed the use of organoids in personalized medicine, but have focused these efforts on targeted treatment of cancers [
27,
46,
47,
48]. The findings reported in these studies are the first to implement this approach for the potential treatment of PitNETs. Collectively, we have developed a relevant human in vitro approach to potentially advance our knowledge as well as our approach to studies in the field of pituitary tumor research. Both the hPITOs and the iPSC
CDH23 may be implemented in studies that strive to (1) define the molecular and cellular events that are crucial for the development of PitNETs leading to CD, and (2) accelerate the identification of effective targeted therapies for patients with CD.
While published studies have advanced our understanding of the molecular mechanisms of the pathogenesis of corticotroph adenomas and elucidated candidate therapeutic targets for CD, these reports fall short of directly informing clinical decisions for patient treatment. Using organoids to screen potential drugs and compounds can potentially improve therapeutic accuracy.
Figure 3 demonstrated a variation in drug responsiveness amongst the organoid lines generated from individual patients. Importantly, there was further divergence in drug responsiveness amongst the individual organoid lines within each pathologically defined corticotroph subtype. For example, hPITOs generated from patients with sparsely granulated corticotroph adenomas (hPIT0s 10, 25, 34, 35) and Crooke’s cell adenomas (hPITOs 7, 33) showed variable responses regardless of similar pathologically defined subtypes. In addition, the response of the tumor cells within the organoids to the standard of care drugs that directly target the pituitary in the body, including mifepristone and cabergoline, was only 50% in hPITO34 and hPITO35, and almost 0% in the other lines, including hPITO7, 10, and 25. These data clearly demonstrate that the inherent patient difference to drug response that is often observed among CD patients is reflected in the organoid culture. This culture system may be an approach that will provide functional data revealing actionable treatment options for each patient. Patient-derived organoids from several tumors have served as a platform for testing the efficacy of anticancer drugs and predicting responses to targeted therapies in individual patients [
27,
46,
48,
49,
50]. An example of the use of organoids in identifying drug responsiveness within an endocrine gland is that of papillary thyroid cancer [
51]. Organoids developed from PTC patients were used as a preclinical model for studying responsiveness to anticancer drugs in a personalized approach [
51]. However, our study is the first report of the use of hPITOs for drug screening. Connecting genetic and drug sensitivity data will further categorize corticotroph subtypes associated with CD. WES analysis of each hPITO line was compared to the results for the corresponding primary adenoma tissues. We showed the concordance rate of exonic variants between the primary tumor tissues obtained from CD patients and the corresponding organoid line. On average, approximately 80% of the variants observed in the CD patients’ adenoma tissues were retained in the corresponding hPITOs.
Pituitary organoids were also developed from iPSCs generated from PBMCs of a CD patient expressing a germline genetic alteration in cadherin-related 23 CDH23 (iPSC
CDH23), a CD patient expressing an MEN1 mutation (iPSC
MEN1), and a healthy individual (iPSC
ctrl). Foundational studies performed by investigators at the genome level have revealed significant knowledge regarding the pathophysiology of CD [
36,
37,
52,
53]. In some instances, CD is a manifestation of genetic mutation syndromes that include multiple endocrine neoplasia type 1 (MEN1), familial isolated pituitary adenoma (FIPA), and Carney complex [
54,
55]. CDH23 syndrome is clinically associated with the development of Usher syndrome, deafness, and vestibular dysfunction [
56]. Several mutations in CDH23 are associated with inherited hearing loss and blindness [
57]. However, none of the variants found in this study were linked to any symptoms of deafness or blindness. A possible explanation is that deafness-related CDH23 mutations are caused by either homozygous or compound heterozygous mutations [
57]. In a study that linked mutations in CDH23 with familial and sporadic pituitary adenomas, it was suggested that these genetic alterations could play important roles in the pathogenesis of CD [
38]. Genomic screening in a total of 12 families with familial PitNETs, 125 individuals with sporadic pituitary tumors, and 260 control individuals showed that 33% of the families with familial pituitary tumors and 12% of individuals with sporadic pituitary tumors expressed functional or pathogenic CDH23 variants [
38]. Consistent with the expected pathology and function of a PitNET from a patient with CD, iPSC
CDH23 organoids exhibited hypersecretion of ACTH, and expression of transcription factors and cell markers were reported in the pathology report for corticotroph PitNETs. Collectively, these findings warrant further investigation to determine whether carriers of CDH23 mutations are at a high risk of developing CD and/or hearing loss. Specifically, clinical investigation is required to determine whether pituitary MRI scans should be adopted in the screening of CDH23-related diseases, including Usher syndrome and age-related hearing loss.
Pituitary organoids generated from iPSCs of a CD patient revealed the existence of cell populations that potentially contribute to the support of PitNET growth and disease progression, as well as an expansion of stem and progenitor cells that may be the targets for tumor recurrence. Organoids derived from both pituitary adenomas and iPSCs exhibited increased expression of stem cell and progenitor markers at both the protein and transcriptomic levels. Unique clusters that were proliferative in the iPSC
CDH23 organoids expressed a hybrid pituitary cell population which was in an epithelial/mesenchymal state (CK20+/VIM+/CXCR4+/Ki67+). In support of our findings, a similar report of a hybrid epithelial/mesenchymal pituitary cell has been made as part of the normal developmental stages of the human fetal pituitary [
58]. Previous studies have suggested that pituitary stem cells undergo an EMT-like process during cell migration and differentiation [
59,
60,
61]. Consistent with our findings are extensive studies using single cells isolated from human pituitary adenomas to show increased expression of stem cell markers SOX2 and CXCR4 [
22,
23,
41,
62,
63]. Within the clusters identified in the iPSC
CDH23 culture were cell populations expressing stem cell markers, including SOX2, NESTIN, CXCR4, KLF4, and CD34. The same iPSC
CDH23 cell clusters, 4, 8, 9, and 11, co-expressed upregulated genes of NOTCH, Hedgehog, WNT, and TGFβ signaling, which are pivotal not only in pituitary tumorigenesis and pituitary embryonic development, but also in ‘tumor stemness’ [
22,
23,
41,
62,
63,
64]. We also noted that clusters of cell populations 5 and 14 unique within the iPSC
CDH23 cultures expressed upregulated genes which were indicative of high proliferation. We observed upregulated expression of the E2F family of transcription factors (E2Fs) E2F1 and E2F7. These findings are of significance, given that there is evidence to show that upregulation of E2Fs is fundamental for tumorigenesis, metastasis, drug resistance, and recurrence [
65]. Within the pituitary adenoma microenvironment, whether these stem cells directly differentiate into pituitary tumors or support the growth of the adenoma is largely unknown. In addition, whether pituitary stem cell populations become activated in response to injury is also understudied. Although the role of stem cells has been identified using a mouse model through implantation of the cells within the right forebrain [
66], the identification of pituitary tumor-initiating stem cells using in vivo orthotopic transplantation models is impossible in mice. Pituitary tumors harboring the stem cells may require engraftment within the environment from which the cells are derived in order to enable growth and differentiation of the tumor. However, it is technically impossible to implant cells orthotopically in the murine pituitary. The pituitary tumor organoid cultures presented in these studies may offer an approach by which isolation, identification, and characterization of this stem cell population is possible. Therefore, we would gain knowledge on the mechanisms of pituitary tumor pathogenesis and reveal potential novel targets for therapeutic interventions by using the iPSC generated pituitary organoid culture.
PitNETs associated with the development of CD cause serious morbidity due to chronic cortisol exposure that dysregulates almost every organ system in the body. Overall, existing medical therapies remain suboptimal, with negative impact on health and quality of life, including considerable risk of therapy resistance and tumor recurrence. To date, little is known about the pathogenesis of PitNETs. Here, we present a human organoid-based approach that will allow us to acquire knowledge of the mechanisms underlying pituitary tumorigenesis. Such an approach is essential to identify targeted treatments and improve clinical management of patients with CD.
CushieBlog 