1. Introduction
Medullary thyroid carcinoma (MTC) is a rare tumor that arises from calcitonin-secreting parafollicular C-cells in the thyroid gland. Approximately 25-40% of MTC cases are hereditary, occurring as part of multiple endocrine neoplasia type 2A or 2B (MEN 2A or MEN 2B) syndromes, while the remaining 60-75% are sporadic cases [
1,
2].
Familial cases are attributed to germline-activating pathogenic variants in the
RET gene, located at 10q11.21 [
2]. Despite a well-established genotype-phenotype correlation between RET pathogenic variants and MEN 2 syndromes, there exists significant variability in the age of onset, disease aggressiveness, and outcomes both within and among families with the same RET pathogenic variant [
1,
2,
3,
4,
5,
6,
7]. Additionally, somatic pathogenic variants in the
RET gene, and less frequently in the RAS genes, are identified in sporadic MTCs [
8,
9,
10].
Both hereditary and sporadic forms of MTC often present with cervical lymph node involvement at diagnosis, occurring in approximately 50% of cases. The primary treatment is extensive and meticulous surgical resection. However, distant metastases are observed in 10-15% of patients, underscoring the need for thorough assessment of disease extent and severity [
1,
2].
In the absence of a standardized risk-stratification system, clinical-pathological prognostic factors -such as age at diagnosis, gender, tumor stage, tumor volume doubling time, regional lymph node involvement, and distant metastases -are used to identify high-risk patients who may benefit from more aggressive treatment [
2,
11].
For patients with progressive or symptomatic metastatic disease not amenable to surgery, radiotherapy and targeted systemic therapies are effective interventions. However, despite their efficacy and generally favorable tolerability, resistance mechanisms have been documented, highlighting the need of ongoing research and drug development to overcome these challenges [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25].
In our prior study, we discovered recurrent somatic copy number alterations (CNAs) affecting the DLK1 locus (14q32.2) in both sporadic and hereditary cases of MTC [
26]. A correlation was observed between DLK1 copy number alterations (CNA) and pathological features indicative of more aggressive disease and poorer outcomes (larger tumor size and advanced AJCC tumor stages).
Delta-like non-canonical Notch ligand 1 (DLK1) is a maternally imprinted and paternally expressed gene that belongs to the DLK1-DIO3 gene cluster on chromosome 14q34. It was initially identified as a non-canonical Notch ligand with established roles during development and supportive functions in several aggressive cancers [
27]. Elevated DLK1 expression has been documented across a spectrum of cancers affecting various organ systems, including the endocrine system, gastrointestinal tract, lungs, liver, kidneys, brain, breast, sarcomas, and both pediatric and adult blood cancers. Its increased expression often correlates with advanced disease stages and is associated with worse outcomes in terms of progression-free survival (PFS), recurrence-free survival (RFS), and overall survival (OS). Moreover, elevated DLK1 levels have been linked to increased self-renewal capacity and resistance to chemotherapy and other treatment modalities across diverse malignancies [
28,
29,
30].
In this study, we conducted a meticulous examination of the expression of a range of well-established stemness markers, alongside potential stem cell markers identified in the MTC literature [
20,
21,
31]. Given that multicellular spheroid generation is often used as an experimental model to measure the self-renewal, multidrug resistance, and multipotent nature of the cancer stem cell (CSC) subpopulation within a tumor or cancer cell line, we further investigated the ability of these cell lines to form multicellular spheroids and express multidrug resistance proteins. We next investigated DLK1 expression in medullary and non-medullary thyroid carcinoma cells, isolated subpopulations of MTC cells expressing DLK1 (DLK1+) and those lacking DLK1 expression (DLK1-) and subjected them to a detailed characterization of their expression of stem cell markers.
3. Discussion
Although MTC is less prevalent compared to differentiated thyroid carcinomas, its prognosis is usually worse. Approximately 50% of patients experience involvement of cervical lymph node metastasis, while about 30% develop distant metastasis, affecting organs such as the liver, lungs, bones, and occasionally the skin and brain [
32].
The primary treatment for MTC involves total thyroidectomy with neck dissection for lymph nodes. Due to their ineligibility for radioactive iodine therapy, patients with progressive disease may undergo radiotherapy or systemic therapy, which encompasses cytotoxic chemotherapy, RET-selective inhibitors, multi-tyrosine kinase inhibitors (TKIs), or immunotherapy [
33]. Although response rates vary, systemic chemotherapy typically yields partial or transient responses. Among cytotoxic drugs, doxorubicin, either as a monotherapy or combined with cisplatin, is the most employed in MTC patients [
34]. However, cytotoxic chemotherapy, while capable of controlling tumor burden, is not considered the first-line therapy for patients with persistent or recurrent MTC due to its low efficacy. Instead, targeted therapies such as tyrosine kinase inhibitors (e.g., vandetanib, cabozantinib) are typically preferred as first-line options for managing persistent or recurrent MTC in many countries [
32,
33,
34,
35,
36]. New Selective RET inhibitors (e.g., Selpercatinib or pralsetinib), have revolutionized the treatment of cancer have been approved as treatment option for cancer with
RET gene alterations. While targeted therapies have demonstrated superior efficacy and tolerability when compared to conventional chemotherapy in treating MTC, they may not achieve complete eradication of cancer cells. Regrettably, clinical responses for advanced MTC are often partial, and patients may eventually develop resistance to these drugs, leading to disease progression over time. Understanding the mechanisms of resistance is crucial for devising more effective and durable therapeutic strategies [
37,
38,
39,
40,
41,
42].
Emerging evidence suggests that drug resistance in cancers may arise from the survival of cancer cells displaying a stem cell phenotype. Indeed, CSCs have been shown to play a critical role in tumor initiation and progression, as evidenced by their ability to initiate new tumors when injected into animal models. CSCs represent a small fraction of cells within the tumor microenvironment, typically ranging from 0.1% to 30% depending on the tumor type and its stage [
43].
These CSCs express canonical markers linked to self-renewal and pluripotency, including core proteins, drug efflux transporters, and surface factors such as CD44 and CD133. Such molecular features significantly contribute to tumor progression and metastasis by enhancing the survival and proliferation of CSCs [
20,
21,
27,
44,
45,
46].
Moreover, CSCs exhibit the unique ability to form multicellular spheroids, potentially enabling their dissemination via bloodstream or lymphatic circulation, thereby facilitating development of distant metastases. This multifaceted role of CSCs underscores their pivotal significance in tumor biology and emphasizes the need for innovative therapeutic strategies targeting cancer progression and metastasis [
47,
48,
49] Consequently, several groups have suggested that CSC markers can be used to predict the prognosis of patients [
50,
51].
There is evidence that CSC markers can distinguish CSCs from other cancer cells, enabling the identification of circulating CSCs before metastasis formation through liquid biopsy. These markers also hold potential as targets for innovative therapeutic approaches [
52,
53].
In the field of thyroid research, substantial experimental evidence has provided support to the hypothesis CSCs are central in both the initiation of cancer and the formation of metastases [
54,
55,
56,
57,
58,
59,
60,
61,
62,
63]. The upregulation of cell surface markers such as CD44, CD133, and ALDH1A1 has been linked to a poorer prognosis and increased resistance to chemotherapy and radiotherapy in both differentiated (DTC) and undifferentiated thyroid carcinomas (UTC) [
21,
64,
65].
We have previously demonstrated recurrent somatic copy number gain of the region (14q32.2), comprising the gene
DLK1, in MTC, which correlates with larger tumor size and advanced tumor stage (III and IV) [
26]. Furthermore, our findings, corroborated by others, indicate that the
DLK1 gene is frequently overexpressed in MTC [
26,
66].
DLK1, located on chromosome band 14q32.2 in humans, is an imprinted gene encoding secreted and membrane-bound isoforms. While its expression is widespread in various human tissues during embryonic development, in adults, it is predominantly restricted to (neuro)endocrine tissues and other immature stem/progenitor cells. Notably,
DLK1 exhibits elevated expression levels in numerous common malignancies, including liver, breast, brain, pancreas, colon, and lung cancers, as well as endocrine-related cancers such as ovarian tumors and pheochromocytoma. Indeed,
DLK1’s elevated expression has been identified in pediatric cancers as well as aggressive and therapy-resistant adult cancers [
28,
67,
68,
69,
70].
Emerging evidence suggests that drug resistance in cancers may stem from the survival of cancer cells exhibiting a stem cell phenotype, and DLK1 may play a crucial role in regulating and maintaining this phenotype [
30,
71].
In this study, to assess their stem cell-like properties, we characterized two MTC cell lines for the expression of 15 known stem cell markers. These lines harbor pathogenic variants in the two most common MEN2A and MEN2B-associated codons (p.C634 and p.M918), enabling us to investigate the relationship between these RET variants and the stem cell phenotype. MZ-CRC-1 cells, which carry the pathogenic variant (p.M918T) associated with MEN2B syndrome, exhibited higher expression levels of most canonical markers compared to TT cells, which possess the p.C634W variant linked to MEN2A.
The protein-array analysis revealed that most markers that are fundamental for the initiation and maintenance of pluripotent states were expressed in both cell lines, with MZ-CRC-1 exhibiting higher levels of most markers. This suggests that MZ-CRC-1 possesses a robust stem cell phenotype, characterized not only by the presence of key pluripotency factors such as OCT3/4, SOX2, and NANOG but also by additional markers that may enhance stemness. Validation of the canonical stem cell markers (OCT3/4, SOX2, NANOG) by flow cytometry confirmed their expression in both cell lines and increased expression in MZ-CRC-1 cells.
Subsequent assessment of CD44, CD133, ALDH1, multidrug resistance proteins (MRP1 and MRP3), and DLK1 revealed varied expression patterns. CD133 and MRP proteins were higher in MZ-CRC-1, while CD44 and ALDH1A1 were elevated in TT cells. Notably, DLK1 was specifically identified in MTC cell lines and showed significantly higher expression in MZ-CRC-1 cells compared to TT cells, indicating its potential role as a marker for MTC.
These findings suggest that both cell lines harbor a subpopulation of cells with a CSC phenotype, potentially influencing tumor progression, aggressiveness, and drug resistance. Differences observed in specific markers like CD133 and MRP proteins in MZ-CRC-1 and CD44 in TT cells may be attributed to varying levels of DLK1 isoforms or interactions with signaling proteins such as Notch or RET, potentially influencing cell maintenance and self-renewal as well as treatment resistance.
Notably, MZ-CRC-1 cells exhibited not only increased DLK1 expression but also enhanced spheroid-forming ability compared to TT cells.
Further investigation involved sorting cells into DLK1-positive and negative populations and assessing stem cell marker levels in both groups. The DLK1-positive subpopulation in both cell lines expressed stemness markers, with notably higher expression in MZ-CRC-1 compared TT cells.
While this study uncovers novel insights, it’s acknowledged that further investigations, including gene knockdown experiments and in vivo analyses, are needed to elucidate DLK1’s role comprehensively. We also acknowledge that besides somatic copy number alteration, other mechanisms may exert an influence on DLK1 expression.
4. Materials and Methods
4.1. Cell Culture
This study utilized two RET-driven MTC cell lines: MZ-CRC-1 (harboring p.M918T mutation) and TT (harboring p.C634W mutation). Additionally, four non-MTC carcinoma cell lines were included: B-CPAP (Papillary thyroid carcinoma), XTC.UC1 (oncocytic cell carcinoma), as well as 8505c and KTC-2 (undifferentiated thyroid carcinoma).
The TT cell line was obtained from the ATCC (Cat# CRL-1803), while the MZ-CRC-1 cell line was generously provided by Prof. Barry Nelkin from Johns Hopkins University. The XTC.UC1 cell line was kindly donated by Prof. Ian Ganly from Memorial Sloan Kettering Cancer Center (New York, USA). The B-CPAP (ACC 273) and 8505c (ACC 219) cell lines were purchased from DSMZ (Germany). The TT cell line was cultured in F-12K medium supplemented with 10% FBS and 100 U/mL penicillin/streptomycin. The MZ-CRC-1 and XTC.UC1 cell lines were cultured in a 1:1 mixture of DMEM and Ham’s F12, supplemented with 10% FBS and 100 U/mL penicillin/streptomycin. B-CPAP and 8505c cell lines were grown in RPMI-1640 medium supplemented with 10% FBS and 100 U/mL penicillin/streptomycin. The KTC-2 cell line was grown in RPMI-1640 medium supplemented with 5% FBS and 100 U/mL penicillin/streptomycin. All cells were maintained at 37°C in a 5% CO2 humidified atmosphere. To ensure the reliability and authenticity of the thyroid cancer cell lines used in our study, we conducted a panel of short tandem repeat (STR) profiling. This step was essential to mitigate concerns about potential cross-contamination and to maintain the integrity of our experimental results.
4.2. Protein Array Analysis
To explore the expression of known CSC markers in the TT and MZ-CRC-1 cell lines, protein extracts isolated from the cells underwent analysis using the antibody-based human pluripotent stem cell proteomic array (Cat # ARY010; R&D Systems, Minneapolis, MN, USA). This array simultaneously detects 15 established stem cell markers (
Table 1 and
Figure S2). Approximately 10
7 TT and MZ-CRC-1 cells were used for protein lysate isolation, and quantification was carried out as previously described [
72].
Nitrocellulose membranes, each spotted with 15 different antibodies in duplicate, were incubated with 300 μg of cellular extract for 16 h at 4°C. After two washes, the membranes were incubated with streptavidin-HRP secondary antibody and revealed with a chemiluminescent reagent (R&D Systems). The intensity score of each duplicated spot was analyzed using the ImageQuant LAS4000 Analyzer (GE Healthcare, Chicago, Illinois, USA) and quantified using ImageQuant TL software (GE Healthcare), as described previously [
72]. The averaged intensity signal of the duplicate spots was calculated by subtracting the averaged background signal (negative control) from each spot, following the manufacturer’s instructions. Fold changes were then calculated based on the average density values of each protein expressed in the MZ-CRC-1 cells, divided by the average density values obtained in the TT cells.
4.3. Flow Cytometry Analysis
FC analysis was employed for validation, assessing the expression of recognized stem cell markers (OCT3/4, SOX2, NANOG), immunophenotyping cells for the expression of both recognized (ALDH1A1, CD44, and CD133), and a hypothetical (DLK1) stem cell markers in the thyroid, and measuring cell size at the single-cell level. FC was also utilized to measure the expression of multidrug resistance markers, MRP1 and MRP3.
To prepare for FC, cells were fixed with 4% paraformaldehyde, blocked with 5% BSA, and then incubated with the appropriate primary antibodies as detailed in
Table S1. Subsequently, cells were incubated with the secondary labeled antibody, goat anti-rabbit IgG-FITC (Thermo Fisher Scientific, Waltham, Massachusetts, USA), and washed in 1X PBS. Two independent experiments were performed, each performed in triplicate. For each flow cytometric analysis, 10.000 events were recorded and analyzed on a BD Accuri C6 flow cytometer with a filter in FL1 (BD Biosciences, Franklin Lakes, New Jersey, USA). The resulting files were exported and analyzed using FlowJo software (Tree Star Inc., Ashland, Oregon, USA).
Negative controls included cells without primary and secondary antibodies, as well as cells incubated solely with the secondary antibody in a blocking buffer. Statistical analyses were based on the median fluorescence intensity (MFI), with results presented as flow cytometry histograms and bar graphs, following the recommendations [
73].
4.4. Hoechst 33342 Efflux Assay
To evaluate the cells’ ability to efflux the fluorescent dye Hoechst 33342, approximately 8 x103 cells per well were cultured in a 96-well plate and exposed to 5 µg/mL of Hoechst 33342 (Thermo Fisher Scientific) for approximately 90 minutes at 37°C. Following the staining period, the cells were washed once with PBS, and the medium was replaced. The plate was then transferred to a SpectraMax Plus M3 microplate reader (Molecular Devices, San José, CA, USA) and the blue fluorescence of Hoechst 33342 was excited at 350 nm, and its fluorescence was measured at 450 nm (Time 0). Subsequently, fluorescence readings were taken at 30, 60, and 90-minute time points (cumulative time points). Two independent experiments were performed in triplicates. Statistical analysis was performed using GraphPad Prism software, and the results were presented graphically to illustrate the efflux kinetics of Hoechst 33342 from the cells over the specified time intervals.
4.6. Generation of Spheroids Cultures
To confirm the ability of MTC cell lines to form multicellular spheroids [
74], 2x10
4 cells were cultured in serum-free medium in 24-well ultra-low attachment plates (Corning, New York, USA). The spheroid formation efficiency (SFE) was monitored using a Nikon Eclipse TE2000 inverted microscope system (Nikon, Melville, New York, USA) equipped with optical system and a TCapture software. Images were captured at 0-, 24-, and 48-h post-seeding. The quantification of spheroids generated by MTC cells was conducted by counting the number and measuring the size of spheroids at specific time points. Two independent experiments were conducted, and each experiment performed in triplicate.
4.7. Western Blot Analysis
Cell homogenates were prepared from 2 medullary and four non-medullary thyroid cell lines and incubated on ice in radio-immunoprecipitation (RIPA) buffer supplemented with protease and phosphatase inhibitor cocktails (Merck, New Jersey, USA). After centrifugation at 14.000 rpm for 10 minutes, protein concentration was quantified using the BCA protein assay kit (Pierce Biotechnology, Rockford, Illinois, USA). Fifty micrograms of reduced protein were loaded onto a 10% polyacrylamide gel. The proteins were then transferred onto a nitrocellulose membrane (Bio-Rad, California, USA) and blocked for 1 hour in TBS with 0.1% Tween and 5% nonfat dry milk before probing with primary antibodies against DLK1 (Ab21682, 1:500) and β-actin (4967, Cell Signaling, 1:10,000). Primary antibodies were incubated overnight at 4°C in blocking buffer, followed by washing with TBST (150 mM NaCl, 2.7 mM KCl, 25 mM Tris, 0.1% Tween). The membrane was then incubated with anti-rabbit horseradish peroxidase-conjugated secondary antibodies (Dako, P0448, 1:1.000) for 1 hour at room temperature. Detection was performed using Immobilon Western Chemiluminescent HRP substrate (Millipore, Massachusetts, USA) and visualized with an ImageQuant LAS 4000 imaging system (GE Healthcare).
4.8. Fluorescence-Activated Cell Sorting (FACS)
To isolate a subpopulation of MZ-CRC-1 and TT cells expressing DLK1, approximately 2 x 107 cells/mL were blocked with 2% BSA, followed by incubation with the primary antibody targeting DLK1 for 1 hour. Subsequently, cells were fluorescently labeled with a secondary antibody before being subjected to FACS analysis using the BD FACS Aria II system (BD Bioscience). The DLK1-positive (DLK1+) and DLK1-negative (DLK1-) subpopulations were collected and subjected to protein isolation. Subsequently, the isolated proteins were analyzed using the Proteome Profiler Human Pluripotent Stem Cell array, as aforementioned. The obtained data was analyzed to determine differential protein expression patterns between the DLK1+ and DLK1- subpopulations.
4.9. Statistical Analysis
Protein array analysis involved calculating fold changes based on average density values of each protein expressed in MZ-CRC-1 cells divided by those in TT cells. Flow cytometry utilized the MFI from 10,000 events, analyzed with an unpaired t-test. The Hoechst 33342 assay was conducted in two independent experiments, each performed in triplicate, and analyzed using a t-test. Spheroid formation was quantified by counting and measuring spheroid number and size using TCapture software, with statistical analysis performed using a t-test in GraphPad Prism software. Relative expression levels between DLK1+ and DLK1- cells in both cell lines were compared after log-transforming the data (Log2) to enhance normality, maintain proportionality in fold changes (both positive and negative), and improve interpretability. A histogram was generated using RStudio software to visualize these comparisons.
Figure 1.
The figure provides a comprehensive analysis of stem cell markers in medullary thyroid carcinoma (MTC) cell lines. (A) shows results from antibody array membrane analysis, while (B) displays a heatmap illustrating the differential expression of 15 analysed proteins. Notably, the MZ-CRC-1 cell line, characterized by the RET p.M918T variant, shows elevated expression of most proteins associated with the stem cell phenotype compared to the TT cell line, characterized by the RET p.C634W variant. (C) (F was employed to assess the expression of recognized stem cell markers (OCT3/4, SOX2, NANOG) in MZ-CRC-1 and TT cell lines. Representative histograms illustrate the red peak indicating target protein expression, distinctly separated from the green peak (secondary antibody) and black peak (negative control). (D) For FC analysis, two independent experiments were conducted in triplicate. About 10,000 events were recorded and reported as median fluorescence intensity (MFI). (E) FC indicates that the percentage of MZ-CRC-1 and TT cells positive for OCT3/4 (40.3%, 26.5%), SOX2 (96.5%, 60.5%), and NANOG (63.8%, 1.7%), respectively. Statistical significance was determined using unpaired t-tests and denoted as follows: >0.05 (ns), *** p ≤ 0.001, and **** p ≤ 0.000.1.
Figure 1.
The figure provides a comprehensive analysis of stem cell markers in medullary thyroid carcinoma (MTC) cell lines. (A) shows results from antibody array membrane analysis, while (B) displays a heatmap illustrating the differential expression of 15 analysed proteins. Notably, the MZ-CRC-1 cell line, characterized by the RET p.M918T variant, shows elevated expression of most proteins associated with the stem cell phenotype compared to the TT cell line, characterized by the RET p.C634W variant. (C) (F was employed to assess the expression of recognized stem cell markers (OCT3/4, SOX2, NANOG) in MZ-CRC-1 and TT cell lines. Representative histograms illustrate the red peak indicating target protein expression, distinctly separated from the green peak (secondary antibody) and black peak (negative control). (D) For FC analysis, two independent experiments were conducted in triplicate. About 10,000 events were recorded and reported as median fluorescence intensity (MFI). (E) FC indicates that the percentage of MZ-CRC-1 and TT cells positive for OCT3/4 (40.3%, 26.5%), SOX2 (96.5%, 60.5%), and NANOG (63.8%, 1.7%), respectively. Statistical significance was determined using unpaired t-tests and denoted as follows: >0.05 (ns), *** p ≤ 0.001, and **** p ≤ 0.000.1.
Figure 2.
The figure describes FC analysis of ALDH1A1, CD44, CD133, MRP1, and MRP3. (A) Shows representative histograms of ALDH1A1, CD44, and CD133 expression, with the red peak indicating target protein expression. This distinct shift separates it from the secondary antibody signal (green peak) and the negative control (black peak). For FC analysis, two independent experiments were conducted in triplicate. About 10.000 events were recorded and reported as median fluorescence intensity (MFI). (B) Shows the percentage of MZ-CRC-1 and TT cells positive for ALDH1A1 (52.1%, 93%), CD44 (47.6%, 48.8%), and CD133 (18.2%, 1%), respectively. Statistical significance was determined using unpaired t-tests and denoted as follows: ** p ≤ 0.01, *** p ≤ 0.001, and **** p ≤ 0.0001.
Figure 2.
The figure describes FC analysis of ALDH1A1, CD44, CD133, MRP1, and MRP3. (A) Shows representative histograms of ALDH1A1, CD44, and CD133 expression, with the red peak indicating target protein expression. This distinct shift separates it from the secondary antibody signal (green peak) and the negative control (black peak). For FC analysis, two independent experiments were conducted in triplicate. About 10.000 events were recorded and reported as median fluorescence intensity (MFI). (B) Shows the percentage of MZ-CRC-1 and TT cells positive for ALDH1A1 (52.1%, 93%), CD44 (47.6%, 48.8%), and CD133 (18.2%, 1%), respectively. Statistical significance was determined using unpaired t-tests and denoted as follows: ** p ≤ 0.01, *** p ≤ 0.001, and **** p ≤ 0.0001.
Figure 3.
(A) Representative Western blot results showing DLK1 expression in TT (line 1) and MZ-CRC-1 cells (line 2) and BCPAP, XTC.UC1, KTC2, and 8505 cells (lines 1-6). β-actin was used as a loading control. (B) Representative histograms from FC analysis depict DLK1 expression in MZ-CRC-1 and TT cells. The red peak in the histograms indicates DLK1 protein expression, clearly distinguishable from the secondary antibody signal (green peak) and the negative control (black peak). (C) Illustrates the proportion of cells expressing DLK1 in MZ-CRC-1 (37%) and TT cells (34,5%). (D) FC analyses was performed in triplicate across two independent experiments and reported median fluorescence intensity (MFI). Statistical significance was assessed using unpaired t-tests and is denoted as >0.05 (ns).
Figure 3.
(A) Representative Western blot results showing DLK1 expression in TT (line 1) and MZ-CRC-1 cells (line 2) and BCPAP, XTC.UC1, KTC2, and 8505 cells (lines 1-6). β-actin was used as a loading control. (B) Representative histograms from FC analysis depict DLK1 expression in MZ-CRC-1 and TT cells. The red peak in the histograms indicates DLK1 protein expression, clearly distinguishable from the secondary antibody signal (green peak) and the negative control (black peak). (C) Illustrates the proportion of cells expressing DLK1 in MZ-CRC-1 (37%) and TT cells (34,5%). (D) FC analyses was performed in triplicate across two independent experiments and reported median fluorescence intensity (MFI). Statistical significance was assessed using unpaired t-tests and is denoted as >0.05 (ns).
Figure 4.
The figure discribed FC analysis of MRP1 and MRP3 and spheroid formation in MTC cells. (A) Representative histograms of FC analysis illustrate the red peak indicating target protein expression, distinctly separated from the green peak (secondary antibody) and black peak (negative control). FC analyses was performed in triplicate across two independent experiments and reported median fluorescence intensity (MFI).(B) Illustrates the concentration of Hoechst 33342 dye in cell supernatant at specified time points. Two independent experiments were conducted in triplicate. (C) Displays representative photographs of spheroid formation in MTC cells taken at distinct time points (24h and 48h after plating). (D) Two independent experiments were conducted in triplicate. The total numbers of multicellular spheroids at each time point are graphically represented. Statistical significance was assessed using an unpaired t-test (**** p < 0.000).
Figure 4.
The figure discribed FC analysis of MRP1 and MRP3 and spheroid formation in MTC cells. (A) Representative histograms of FC analysis illustrate the red peak indicating target protein expression, distinctly separated from the green peak (secondary antibody) and black peak (negative control). FC analyses was performed in triplicate across two independent experiments and reported median fluorescence intensity (MFI).(B) Illustrates the concentration of Hoechst 33342 dye in cell supernatant at specified time points. Two independent experiments were conducted in triplicate. (C) Displays representative photographs of spheroid formation in MTC cells taken at distinct time points (24h and 48h after plating). (D) Two independent experiments were conducted in triplicate. The total numbers of multicellular spheroids at each time point are graphically represented. Statistical significance was assessed using an unpaired t-test (**** p < 0.000).
Figure 5.
The figure illustrates the cell size (FSC) and granularity (SSC) of MZ-CRC-1 and TT cell lines sorted by DLK1 expression, either positive (DLK1+) or negative (DLK1-). The bar chart presents the mean cell size of each populations highlighting that the DLK1+ subpopulation (P2) exhibits significantly smaller size compared to DLK1- subpopulation (P3) (** p ≤ 0.01, ***p<0.001). .
Figure 5.
The figure illustrates the cell size (FSC) and granularity (SSC) of MZ-CRC-1 and TT cell lines sorted by DLK1 expression, either positive (DLK1+) or negative (DLK1-). The bar chart presents the mean cell size of each populations highlighting that the DLK1+ subpopulation (P2) exhibits significantly smaller size compared to DLK1- subpopulation (P3) (** p ≤ 0.01, ***p<0.001). .
Figure 6.
(A) presents results from an antibody array membrane analysis comparing DLK1-positive (DLK1+) and DLK1-negative (DLK1-) subpopulations from the MZ-CRC-1 and TT cell lines. (B) displays a heatmap illustrating the differential expression of 15 analysed proteins between the DLK1+ and DLK1- subpopulations (see
Table 2).
Figure 6.
(A) presents results from an antibody array membrane analysis comparing DLK1-positive (DLK1+) and DLK1-negative (DLK1-) subpopulations from the MZ-CRC-1 and TT cell lines. (B) displays a heatmap illustrating the differential expression of 15 analysed proteins between the DLK1+ and DLK1- subpopulations (see
Table 2).
Figure 7.
This figure presents the analysis of stem cell markers in DLK1-positive (DLK1+) compared to DLK1-negative (DLK1-) subpopulations in both MZ-CRC-1 (
A) and TT (
B) cell lines. The data presented in
Table 2 were log-transformed, and a histogram was generated using RStudio software to compare the relative expression levels between DLK1+ and DLK1- cells in both cell lines. The bar chart depicts the fold change (log2) on the y-axis, highlighting 15 stemness markers. Proteins significantly upregulated in the DLK1+ subpopulation are represented in light gray, while those showing downregulation are depicted in darker gray or black for each cell line.
Figure 7.
This figure presents the analysis of stem cell markers in DLK1-positive (DLK1+) compared to DLK1-negative (DLK1-) subpopulations in both MZ-CRC-1 (
A) and TT (
B) cell lines. The data presented in
Table 2 were log-transformed, and a histogram was generated using RStudio software to compare the relative expression levels between DLK1+ and DLK1- cells in both cell lines. The bar chart depicts the fold change (log2) on the y-axis, highlighting 15 stemness markers. Proteins significantly upregulated in the DLK1+ subpopulation are represented in light gray, while those showing downregulation are depicted in darker gray or black for each cell line.
Table 1.
Expression of Stem Cell Markers in Medullary Thyroid Carcinoma cell lines.
Table 1.
Expression of Stem Cell Markers in Medullary Thyroid Carcinoma cell lines.
Protein |
Average intensity in TT cells |
Average intensity in MZ-CRC-1 cells |
Fold Change* |
P-value |
OCT3/4 |
173 |
941 |
5.439306 |
0,009 |
SOX2 |
283 |
389 |
1.374558 |
0.0428 |
NANOG |
398 |
786 |
1.974874 |
0.031 |
OTX2 |
452 |
570 |
1.261062 |
0.0423 |
E-Cadherin |
7956 |
893 |
0.112242 |
0.0017 |
Snail |
846 |
589 |
0.696217 |
0.0059 |
GSC |
430 |
479 |
1.113953 |
NS |
VEGF |
339 |
662 |
1.952802 |
0.0139 |
SOX17 |
249 |
492 |
1.975904 |
0.0282 |
FOXA2 |
218 |
435 |
1.995413 |
NS |
TP63 |
272 |
341 |
1.253676 |
NS |
PDX-1 |
359 |
434 |
1.208914 |
0.0247 |
GATA-4 |
660 |
469 |
0.710606 |
NS |
AFP |
163 |
192 |
1.177914 |
NS |
HCG |
335 |
2048 |
6.113433 |
0.0035 |
Table 2.
Expression of Stem Cell Markers in Medullary Thyroid Carcinoma Cell Lines Stratified by DLK1 Status.
Table 2.
Expression of Stem Cell Markers in Medullary Thyroid Carcinoma Cell Lines Stratified by DLK1 Status.
Protein |
Average intensity in TT DLK1- cells |
Average intensity in TT DLK1+ cells |
P-value |
Average intensity in MZ-CRC-1 DLK1- cells |
Average intensity in MZ-CRC-1 DLK1+ cells |
P-value |
OCT-3/4 |
528 |
689 |
0.0255 |
544 |
771 |
0.0041 |
NANOG |
628 |
731 |
0.0119 |
569 |
859 |
0.0042 |
SOX2 |
563 |
639 |
0.0241 |
440 |
864 |
0.013 |
OTX2 |
810 |
818 |
NS |
995 |
936 |
NS |
E-cadherin |
1495 |
3194 |
0.038 |
985 |
1544 |
0.011 |
Snail |
781 |
703 |
NS |
708 |
844 |
0.0091 |
GSC |
669 |
695 |
NS |
884 |
926 |
NS |
VEGF |
699 |
1156 |
0.0033 |
783 |
1002 |
NS |
SOX17 |
703 |
700 |
NS |
716 |
907 |
NS |
FOXA2 |
803 |
1302 |
0.0245 |
856 |
941 |
NS |
TP63 |
650 |
691 |
NS |
887 |
792 |
NS |
PDX-1 |
618 |
660 |
NS |
679 |
811 |
NS |
GATA-4 |
1511 |
1495 |
NS |
2192 |
1587 |
0.0105 |
AFP |
603 |
533 |
NS |
619 |
608 |
NS |
HCG |
757 |
901 |
NS |
988 |
962 |
NS |