3D melanoma spheroid model for the development of positronium biomarkers

Scientific Reports volume 13, Article number: 7648 (2023) Cite this article

299 Accesses

1 Altmetric

Metrics details

It was recently demonstrated that newly invented positronium imaging may be used for improving cancer diagnostics by providing additional information about tissue pathology with respect to the standardized uptake value currently available in positron emission tomography (PET). Positronium imaging utilizes the properties of positronium atoms, which are built from the electrons and positrons produced in the body during PET examinations. We hypothesized that positronium imaging would be sensitive to the in vitro discrimination of tumor-like three-dimensional structures (spheroids) built of melanoma cell lines with different cancer activities and biological properties. The lifetime of ortho-positronium (o-Ps) was evaluated in melanoma spheroids from two cell lines (WM266-4 and WM115) differing in the stage of malignancy. Additionally, we considered parameters such as the cell number, spheroid size and melanoma malignancy to evaluate their relationship with the o-Ps lifetime. We demonstrate pilot results for o-Ps lifetime measurement in extracellular matrix-free spheroids. With the statistical significance of two standard deviations, we demonstrated that the higher the degree of malignancy and the rate of proliferation of neoplastic cells, the shorter the lifetime of ortho-positronium. In particular, we observed the following indications encouraging further research: (i) WM266-4 spheroids characterized by a higher proliferation rate and malignancy showed a shorter o-Ps lifetime than WM115 spheroids characterized by a lower growth rate. (ii) Both cell lines showed a decrease in the lifetime of o-Ps after spheroid generation on day 8 compared to day 4 in culture, and the mean o-Ps lifetime was longer for spheroids formed from WM115 cells than for those formed from WM266-4 cells, regardless of spheroid age. The results of this study revealed that positronium is a promising biomarker that may be applied in PET diagnostics for the assessment of the degree of cancer malignancy.

Over the past few decades, three-dimensional (3D) cell cultures have been widely used as in vitro models, which can bridge the gap between in vitro and in vivo cell conditions1. The comparison of the 3D cell culture to a cell monolayer revealed some specific physiological and morphological characteristics, such as cell-to-cell and cell–matrix interactions, cell signaling, proliferation and necrosis. Unlike a monolayer cell culture, a 3D spheroid is an appropriate model to mimic the real tumor cell environment and nutrient diffusion rate between cells. Multicellular tumor spheroids allow us to study the biochemical mechanism of cell growth, enzymatic reactions, and various treatment modalities2,3,4,5.

Melanoma is a prevalent type of skin cancer that has been categorized as one of the most lethal cancers. Melanoma, as the most lethal cancer, is a multifactorial disease in which both genetic susceptibility and environmental exposure, predominantly to ultraviolet light, play important roles6. The appointed environmental risk factors for melanoma cancer are exposure to solar ultraviolet radiation (UVB) and sunburns (UVA). UVB irradiation causes direct DNA damage, which leads to DNA strand breaks. UVB also promotes melanoma cell survival, angiogenesis, and invasion due to increases in the penetration of macrophages and neutrophils into skin cells. UVA causes DNA damage through the production of free radicals, which induce oxidative stress in melanocytes. The risk for melanoma can also be associated with inherited mutations and somatic mutations, but the genetic tendency is responsible for only a small number of cases. As the effectiveness of melanoma treatment at advanced stages is low, there is a constant need for the development of new targeted therapies, immunotherapy, and combined therapies7 (Fig. 1).

Schematic view of melanocytes and the formation of a melanoma lesion in the epidermal layer. (A) Normal melanocytes have a regular dendritic shape and form many branching processes extending between numerous keratinocytes (left). One melanocyte reaches between 36 and 40 keratinocytes, forming an epidermal melanin unit (EMU), with a balanced proportion of one melanocyte to each of 8–10 keratinocytes in the basal layer of the epidermis. In a melanoma lesion (right), melanocytes lose their dendricity and become malignant and amoeboid in shape, changing their cell–cell contacts by expressing different adhesion molecules (N-cadherin instead of E-cadherin expressed by melanocytes)7. Melanosomes store the melanin granules produced by melanocytes, which can be distributed among surrounding keratinocytes to protect them from UV radiation damage. (B) Pictorial illustrations of positron annihilations in the melanin molecule. Carbon, oxygen, hydrogen, nitrogen, p-Ps and o-Ps atoms are indicated in colours explained in the legend. Dashed arrows indicate photons from direct electron–positron annihilation. Red arrows indicate photons from o-Ps atom self-annihilation. Blue arrows show photons from p-Ps self-annihilation. Green and brown arrows indicate o-Ps decay via the pick-off and conversion processes, respectively.

Although numerous studies have been dedicated to diagnosing melanoma in the early stages, more investigation and follow-up are needed. The rate of metastasis in patients due to melanoma is high, with 350,000 new cases per year, and the number of new cases is significantly increasing each year8,9,10.

In this study, a 3D spheroid model of two melanoma cell lines, WM266-4 and WM115, with different stages of malignancy was evaluated. The melanoma cell lines WM266-4 and WM115 have various biological properties that make them different in proliferation, migration rate, and cell size. The activation of advanced glycation end products (AGEs), which are produced during the combination of lipids and proteins with sugar in human metabolism, and proteins such as the receptor for AGEs (RAGE) and c-Jun N-terminal kinases (JNK) in the WM266-4 line makes this cell line more invasive than the WM115 line11. Moreover, the expression of the SLC7A5 gene (variant of Large Amino-acid Transporter 1) in WM115 is higher than in WM266-4, in contrast to the expression of the SLC7A8 gene, which is lower in WM115 than in WM266-412. These genes are involved in amino acid import and are essential for cell pigmentation, and the WM115 cell line is darker than the WM266-4 cell line.

In this study, we hypothesized that the difference between the grade of malignancy of WM115 and WM266-4 melanoma cancers present at the level of cell physiology can be probed by positronium biomarkers. A positronium is an exotic atom built of positrons and electrons13. The positronium is also formed in the intramolecular spaces during PET diagnosis13. In the tissue, the positronium may be formed and trapped in the free voids of the intramolecular spaces, as shown pictorially in Fig. 114. The positron emitted from the radionuclide (e.g., 18F in the PET diagnosis or 22Na in the typical Positron Annihilation Lifetime Spectroscopy (PALS) experiments) penetrates the object, and after losing the energy, annihilates with an electron from the molecules constituting the cells. The electron–positron annihilation into photons may occur directly (e + e− → photons (black dashed arrows in Fig. 1B)) or viaa positronium atom [e + e− → positronium → photons (solid arrows in Fig. 1B)]. In a quarter of cases, the positronium is formed as a short-lived (125 ps) state called the para-positronium (p-Ps), and in three-quarters of cases, it is formed as a long-lived (142 ns) ortho-positronium (o-Ps). Para-positronium (indicated in blue in Fig. 1B) decays predominantly into two photons (blue arrows), and ortho-positronium decays in vacuum predominantly into three photons (red arrows). However, in the intramolecular voids, ortho-positronium undergoes processes such as pick-off (the positron from o-Ps annihilates into two photons (green arrows) with an electron from the surrounding molecule) and conversion into para-positronium via interaction with molecules such as oxygen molecules. The resulting para-positronium decays into two photons (brown arrows). The range of the o-Ps mean lifetime variation is significant, and it changes from 142 ns in vacuum to 1.8 ns in water15,16. Therefore, the ortho-positronium lifetime depends strongly on the molecular environment: the nanomolecular structure and concentration of the bioactive molecules13,14,15,16,17,18,19. The properties of positronium in biological samples have been scarcely studied thus far. Only recently were the first studies of positronium in 3D cell structures performed by culturing cells on collagen scaffolds20. Other studies on skin cancer cells (basal and squamous cell carcinomas) were performed using a low-energy positron beam that limited skin penetration21,22,23.

The first in vitro research on the positronium lifetime in tissues from patients indicated promising results in view of the application of positronium as a biomarker for the diagnosis of uterine cancer and myxoma cancer and as a biomarker of hypoxia24,25,26,27. Recently, positronium was proposed as a novel biomarker for the in vivo assessment of tissue pathology13,25 that can be imaged using a newly developed positronium imaging method15,19,26,28 when applying prompt photon radionuclides28,29,30,31 and high-sensitivity PET scanners28. Moreover, the advent of total-body PET systems characterized by high sensitivity32,33,34,35 will enable the simultaneous application of positronium imaging to standard metabolic imaging during positron emission tomography29,36,37.

The aim of this study was to determine whether positronium imaging would be sensitive to the in vitro discrimination of tumour-like three-dimensional structures (spheroids) built from melanoma cell lines with different cancer activities and biological properties.

Two melanoma cell lines, WM266-4, a human malignant melanoma cancer cell line, and WM115, a cell line from a primary melanoma, were purchased from ESTDAB Melanoma Cell Bank (Tübingen, Germany) and cultured as we previously described38. A Luna-II™ automated cell counter (Logos Biosystems, Inc.) was used to determine the cell counts and viability before cell seeding.

Both cell lines, WM266-4 and WM115, were seeded in 5D spherical plates (5D sp5dplate, Kugelmeiers, Switzerland) to form spheroids39. The 5D microplate has 24 wells, 12 wells with nanocoated surface for spheroid formation, and 12 wells as a control. In our study, control wells were not used. Each well contains 750 microcavities (with 509 µm diameter and 320 µm depth) that are separated from each other by sharp borders, and these borders prevent cell migration from one microcavity to another; thus, 9000 spheroids with uniform shapes and diameters can be cultured in a single plate.

For cell seeding, 0.5 ml of complete medium was added to each well. Then, an extra 0.5 ml medium including 1,125,000 cells/well (1500 cells/microwell) was added to each well of the plate. The cells in the Falcon tube were resuspended to distribute them in the whole medium and then added to the wells.

After cell seeding, plates were maintained in an incubator at 37 °C and 5% CO2 with humidity. The cell culture medium was renewed every day after the spheroids were formed. Two suitable time points of spheroid growth were chosen to measure the mean o-Ps lifetime. Spheroid morphology and proliferation rate were also determined under an optical microscope (Olympus, IX-LWPO, T2, Japan) on days 4 and 8 after cell seeding. Image analysis was conducted by ImageJ software.

First, spheroids were dissociated to a single-cell suspension with trypsin/EDTA (Gibco, cat. no. 25200072, Waltham, MA USA). Then, 100 µl of trypsin/EDTA Spheroids was incubated for 10–20 min at 37 °C and pipetted several times to separate the cells. Then, the cells were centrifuged at 300g for 3 min, the supernatant was removed, and 100 µl of medium (Gibco, cat. no. 10010056, Waltham, MA, USA) was added to the cells and pipetted several times to separate the cells completely. In the final step, 10 µl of cells was added to 10 µl trypan blue and counted (the Luna-II™ cell counter, Logos Biosystems, Aligned Genetics, Inc). The workflow is shown in Fig. 2. Each viability test was performed on 4 different samples in each measurement under the same conditions.

Workflow for investigation of spheroid viability after harvesting 3D spheroids from 5D microplates.

After 15 min in the incubator, the cells were checked to ensure that there were no aggregates. After the viability test using the cell counter, the cell cluster map was evaluated to determine the percentage of clusters. In our experiments, cluster maps showed a high percentage of dissociation, above 90%. Afterward, the size of cells for both cell lines in 3D shape was checked by cell counter and microscopic investigations. For estimation of the cell size,

A Luna II ™ cell counter was used. Regarding the histogram of size distribution, the mean diameter of cells was estimated. For comparison, the size of cells in 2D cultures was also determined. For this purpose, spheroids were removed from each well using a 3 ml pipette with a large bore to avoid destroying the spheroid structure. Then, spheroids were poured onto a 15 ml Falcon tube and centrifuged at 500g for 7 min. In the next step, the supernatant was removed, and 1 ml fresh medium was added to the spheroids.

To evaluate glucose uptake, a fluorescence-emitting d-glucose probe 2-NBDG (2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxyglucose) (Thermo Fisher Scientific, cat. no N13195) was used. To assess the oxygen uptake rate distribution in the spheroids, WM266-4 and WM115 cells were seeded at densities of 1000 and 2000 cells/drop. On days 4 and 8, the spheroids were transferred to glass-bottom dishes, and 20 μl of 2-NBDG (200 μM) was added to each spheroid. The spheroids were incubated at 37 °C for 45 min. Finally, spheroids were washed with PBS, and spheroid images were taken with a Nikon Eclipse Ti-E microscope coupled to an A1 scanning confocal system (Nikon, Japan) at a 494 nm/551 nm fluorescent wavelength.

Oxygen uptake rate distribution was determined using a hypoxia kit, Image-IT™ Green Hypoxia Reagent (Thermo Fisher Scientific, cat. no I14834). To assess the oxygen uptake rate distribution in the spheroids, WM266-4 and WM115 cells were seeded at densities of 1000 and 2000 cells/drop. On days 4 and 8, the spheroids were transferred to glass-bottom dishes, and 20 μl of Image-IT™ Hypoxia Reagent (10 μM) was added to each spheroid. The spheroids were incubated at 37 °C for 1 h. The hypoxia dye was exchanged with fresh growth medium, and the spheroids were placed again in a cell culture incubator for the next 4 h. Spheroid images were taken with a Nikon Eclipse Ti-E microscope coupled to an A1 scanning confocal system (Nikon, Japan) at a 488 nm/520 nm fluorescent wavelength.

To obtain reliable and precise results of the positronium lifetime in spheroids, we used spheroids without any medium, supernatant, or chemical compounds. In this method, harvested spheroids from the microplates were prepared for positronium lifetime measurement after centrifugation and removal of the medium completely from the spheroids.

The positronium lifetime was measured using a spectrometer consisting of two vertically arranged BaF2 plastic scintillators (SCIONIX, Holland) and two photomultipliers with serial numbers SBO696 and SBO697 (Hamamatsu, Japan) powered by a high voltage power supply (CAEN SY4527). The signals from photomultipliers were read out and analyzed by the 6000A data analyzer oscilloscope (Le Croy). Spheroids were irradiated by positrons emitted from the 22Na radionuclide (with an activity of approximately 1 MBq) placed in Kapton foil. A dedicated aluminum chamber was used as a container for spheroids, and a holder was connected to a heater to keep the cells at 37 °C. The measurement setup is shown pictorially in Fig. 3.

Spheroids surrounding the 22Na source are located in the dedicated chamber between two BaF2 detectors. Spheroids are adjacent to the 22Na source, and there is no space or bubble between them. For each measurement, 106 events with the coincident registration of 511 keV photons and 1274 keV photons were collected.

The 22Na radionuclide after emission of the positron transforms to the excited state of the 22Ne isotope that deexcites (on average after 2.6 ps) via emission of the 1274 keV gamma quantum29,30. The positron loses energy while passing through the cells and eventually annihilates with the electron into two back-to-back 511 keV gamma quanta. Positron–electron annihilation may proceed directly or via creation of the positronium. The time between the emission of the positron and its annihilation is measured by the registration of the 1274 keV deexcitation photon and one of the 511 keV annihilation photons. The sample with the source is positioned in a way (see Fig. 3) that enables coincident registration of 1274 keV gamma and one 511 keV photon, and it prevents the coincident registration of both 511 keV photons flying back-to-back. An exemplary lifetime spectrum determined as a result of the measurement is shown in Fig. 4.

Example positronium lifetime spectrum. Experimental data (black histogram) with superimposed histograms resulting from the fit of the sum of the exponential function convoluted with the detector resolution, performed by means of the PALS Avalanche program25,40. The first component (yellow line) shows the contribution from p-Ps (mean lifetime: 0.125 ns), the second component (green line) originates from annihilations in the source (Kapton foil) (0.374 ns), the third component (light blue) shows the free annihilation lifetime (0.395 ns), and the fourth component (dark blue) illustrates the contribution from o-Ps. The sum of all contributions resulting from the fit is shown as a red curve.

Before each measurement, the system was calibrated using an empty chamber including only the 22Na source. To perform the measurement, spheroids were transferred from the centrifuge tube to the chamber by using a scalper, and 22Na radioactive sources were placed between the samples. Then, the spheroid sample was placed adjacent to the source, with no air between the source and spheroids. Then, the chamber was placed inside a holder that was connected to a heater to keep the cells at 37 °C. Finally, the holder was located between two BaF2 detectors, as indicated in Fig. 3.

The lifetime of positronium in spheroids at two different time points, the 4th and 8th days after seeding cells, was evaluated based on the 106 events collected for each studied case. In this study, only spheroids without medium, water, and any chemicals were evaluated.

The positronium lifetime was extracted from each recorded spectrum by fitting the sum of four exponential functions convoluted with the detector resolution function40,41. The exemplary spectrum with the result of the fit is shown in Fig. 4. The fitted components correspond to the decay of para-positronium (yellow curve), direct annihilation in the source/Kapton foil (green curve), and annihilation of ortho-positronium (dark blue). The full experiment was repeated three times.

The exemplary spheroids grown in the microplate are shown in Fig. 5. The WM266-4 cell line formed more spherical and concentrated spheroids than the WM115 cell line, as previously reported38. Spheroids in both cell lines showed an increase in their size and circularity over time. Generally, the rate of growth and size of spheroids depend on the size of cavities and the type of plates that have been seeded. The larger the microcavities are, the larger the formed spheroids⁠. Although the size of spheroids can be controlled by initial culture seeding, the scales of microwells can also affect the spheroid diameter42.

Comparison of spheroids from different melanoma cell lines. (A) Microscopic images of two different cell lines, WM266-4 and WM115, on days 4 and 8 of culturing. The density and circularity of spheroids increased during the culture time. (B) As a function of culturing time, WM266-4 spheroids showed a higher volume than WM115 spheroids. (C) The number of cells in the plate was increased during the culture time for both cell lines. The number of cells in one plate was calculated using the Luna II cell counter after spheroid harvesting.

Figure 5B,C shows that WM266-4 and WM115 spheroids grew over time. WM266-4 spheroids exhibited a faster proliferation rate than WM115 spheroids, which returned to their malignant characteristics.

Figure 5B shows that there was an increase in the number of cells from cell seeding until the 8th day after culturing. The division rate of WM266-4 cells in spheroids was higher than that of WM115 cells. WM266-4 resulted in 1.5- and 2.74-fold increases in cell number after 4 and 8 days, respectively. The number of cells in WM115 spheroids increased 1.4- and 1.7-fold after 4 and 8 days, respectively. A total of 36,000 spheroids were used in each positronium lifetime measurement.

The WM266-4 spheroids comprised cells with a diameter of 15.70 ± 0.10 μm at day 4 after cell seeding and 15.92 ± 0.08 μm at day 8, while the WM155 spheroid cell diameter was equal to 16.66 ± 0.20 μm at day 4 and 17.28 ± 0.25 μm on day 8. The mean size of stained cells in 2D cell culture was smaller than the size of cells from the 3D spheroids, 14.65 ± 0.09 μm and 16.27 ± 0.10 μm for WM266-4 and WM115, respectively.

The number of cells in both cell lines increased over time. The observed growth was faster for the WM266-4 cell line than for the WM115 cell line. WM115 has a longer doubling time than WM266-4, approximately 7.5 and 6 days, respectively, which means that cells of the WM115 line spend more time in their cell cycle than cells from the WM266-4 line. The doubling time (DT) of spheroids is usually characterized as the real tumor doubling time and is calculated by using spheroid volume via:

where V1 and V2 are spheroid volumes at times t1 and t2 = t1 + t after seeding, respectively43,44,45.

According to the fluorescence intensity analysis (Fig. 6), the region between 100 and 200 µm from the center of a spheroid was considered the outer layer, also called the proliferation rim, and the region between 50 and 100 µm of the inner part of the spheroid was considered as a necrotic core. However, the size of different layers of a spheroid depends on the size of the spheroid. Fluorescence intensity after 4 days in spheroids obtained from an initial number of 1000 cells at the edge of the WM266-4 spheroids was not significantly higher than that for WM115 (p = 0.20), while in the deeper region of spheroids, a higher intensity of glucose uptake was observed in WM266-4 spheroids than in WM115 spheroids (p = 0.00001). Larger spheroids formed from 2000 cells showed higher glucose uptake after 4 days, and the intensity was significantly higher in the proliferation rim (p = 0.04) and in the core area in WM266-4 spheroids vs. WM115 spheroids (p = 0.0018). On day 8, WM266-4 and WM115 spheroids showed a significant difference in the intensity of glucose uptake in the proliferation rim (p = 0.007) but not in the core (p = 0.064) (Fig. 6).

Determination of glucose distribution in melanoma spheroids. Glucose concentration is observed equally distributed in the proliferation rim in WM266-4 spheroids (A) and more dispersed WM115 cell lines (B) with significant differences in the larger and older spheroids. Plot of glucose distribution in WM266-4 (C) and WM115 (D) spheroids at two different culture time points as assessed with the 2-NBDG probe.

On day 4 of culture, the progression of hypoxia measured as fluorescence intensity in spheroids with an initial cell number of 1000 was not significantly different (p > 0.05), while on day 8, the degree of hypoxia in the center of WM266-4 spheroids was significantly higher than that of WM115 spheroids (p = 0.0001). In the larger spheroids during culture, the degree of hypoxia in WM266-4 spheroids was significantly higher (p < 0.001) than that in WM115 spheroids, while at day 8, WM266-4 and WM115 spheroids showed a significant difference in the progression of hypoxia in their center (p = 0.00048) but not in the outer part (p = 0.24) (Fig. 7).

Determination of hypoxia progression in melanoma spheroids. The hypoxic region is observed in the center of spheroids in (A) WM266-4 and (B) WM115 cell lines with significant differences in the larger and older spheroids. Plot of hypoxia distribution in WM266-4 (C) and WM115 (D) spheroids at two different culture time points as assessed with the Image-IT™ Green Hypoxia probe.

The mean lifetime of ortho-positronium atoms in the spheroids was established for two different cell lines characterized by different degrees of malignancy: WM115 and WM266-4. Both measurements were performed for spheroids at two different time points: on days 4 and 8 after seeding. The viability tests were evaluated before and after each PALS measurement and showed that spheroid viability remained constant above 85% until the 8th day. Thus, the spheroids were in good and stable condition during the experiment.

The measurement of the positronium lifetime was repeated three times. In each measurement, spheroids were harvested from the plate and placed inside the chamber.

Figure 8 demonstrates the results of the mean o-Ps lifetime and distribution in 3D melanoma spheroids with different malignancy levels at two different ages.

Comparison of the mean o-Ps lifetime and intensity for different melanoma cell lines. (Left) The lifetime of o-Ps in WM266-4 (black squares) and WM115 (red dots) spheroids in two different ages, 4 and 8 days after cell seeding. (Right) Intensity of o-Ps production in WM266-4 spheroids (black squares) and WM115 spheroids (red dots) as a function of time. To compare the mean o-Ps in relation to spheroid metabolism, please see supplementary file 1.

The obtained results given in Table 1 show that WM266-4 spheroids with more malignancy, proliferation rate, and high concentration of cells in spheroids have a shorter ortho-positronium mean lifetime than WM115 cells, therefore, smaller free intermolecular voids, than WM115 cells from the primary tumor with less division rate and lower concentration of cells over time. Although both cell lines show a gradual decrease in the lifetime of o-Ps during the time of culturing, o-Ps show a longer lifetime in WM115 than in WM266-4 spheroids. The intensity in WM266-4 spheroids remained almost constant while WM115 spheroids present a decrease in the intensity of < 3SD >, which can be considered as changes at the molecular level in WM115 spheroids cells.

Recently, positronium imaging was introduced14,15,19,46, and the first in vitro positronium images were demonstrated25,26, opening new possibilities for the improvement of cancer diagnosis by using positronium as a tissue pathology biomarker13,14.

The main aim of this study was to test the hypothesis that positronium may serve as a novel biomarker for assessing the different cancer activities and biological properties in cancer cells and to perform the first studies of positronium properties in 3D cell spheroids. To test this hypothesis, 3D spheroids were used because they can mimic the structure of real tumors under physiological conditions. The morphology and cell–cell interactions in 3D spheroids are completely different than those in monolayer cell culture. Moreover, they hold certain advantages over standard research methods, including low cost of use, high reproducibility, timesaving, and reduced need for laboratory animal models. The unique properties of 3D tumor spheroids make them invaluable for biological experiments and drug tests in a variety of experimental studies focused on chemotherapy and radiotherapy1,2,3,4,5,47,48.

The cell cycle time in 2D cells and 3D spheroids is significantly different because cells in monolayer cell culture have homogenous biological conditions; in a confluent monolayer, most of the cells are in the same cell cycle and have sufficient accessibility to nutrients and oxygen. In 3D cell culture with a multilayer structure, cells are in different cell cycles, and accessibility to necessary nutrients for survival depends on their distance from the surface49. During the cell cycle, cells in the G1 phase grow up to the end of the G2 phase when mitosis and division start after passing the checkpoints. In spheroids, space for proliferation is limited, and a necrotic core is formed surrounded by a quiescent layer where cells are arrested in the G1 phase of the cycle, and a proliferating rim at the outermost part of the spheroids where cells can stay in the S/G2/M phase50. A significant difference in the dynamics of spheroid formation and cell proliferation was observed in our study. We observed that the more malignant cell line (WM266-4) formed larger spheroids (both on days 4 and 8 in culture) with a higher number of cells (Fig. 5B,C), and the cell diameter in the spheroids was lower than that in the primary melanoma cell line (WM115).

Prior positronium lifetime research has been mostly performed on different types of tissues and monolayer cell cultures, such as skin cancer and colorectal cancer cells20,21,22,23,46,51,52. In this novel research, the lifetime of positronium was evaluated in 3D spheroids rather than 2D cell culture or tissue. Although much research has been performed on positronium in materials and biological samples, research using 3D cell aggregates is limited to the determination of the o-Ps lifetime in 3D colorectal cancer cell aggregates with a mixture of cells and collagen20. The o-Ps mean life time was shown to be a sensitive indicator of mean void volume in 3D cancer cultures stimulated with the growing factor affecting epithelial cell stimulation (TGF-β). Colonies treated with TGF-β showed growth with the mean o-Ps 2–3 weeks after stimulation, which decreased to the value under control conditions in the late phase of culture19. This o-Ps distribution indicated a possible relationship between cell dynamics and mean molecular voids. The diffusion behavior of small molecules (Å) has been investigated using PALS in different biological samples (hair, feather, cotton and silk), showing that o-Ps are sensitive to void investigation in biological polymers53,54. The lifetimes of o-Ps in the cell can be different from the lifetime measured in the case when the cell is in the medium containing water, collagen, or other chemical compounds because the cell surface, cell growth, and the shape of cells can be changed.

In our study, we focused on the early stage of spontaneous proliferation of two different cancer cell lines. These cell lines differ in cell size, with WM155 cells being larger than WM266-4 cells. Thus, during spheroid formation, WM266-4 cells compensate for their smaller size with a higher rate of proliferation to form larger spheroids. Their doubling time was shorter than that of WM115 cells, and WM266-4 cells spent less time in their cell cycle than WM115 cells. These important differences in the cell dynamics in the early stage of spheroid growth were indicated by the shorter o-Ps lifetime for the WM266-4 spheroids on day 4 of culture. In our previous study, we compared the percentage of viable cells as assessed by the FDA/PI fluorescent assay (fluorescein acetate and propidium iodide) in situ, and we did not observe any difference between cell viability in the early phase of spheroid formation (on day 4) for these two melanoma cell lines, but a larger necrotic core was recognized in the WM115 spheroids in comparison to WM266-4 (29.36 vs. 13.66% of dead cells)38. These results are in concordance with the glucose and hypoxia distribution assays, where the core of the spheroid was recognized as the region with the lowest glucose concentration and the highest hypoxic conditions. Decreased glucose concentration (possibly caused by reduced diffusion or increased anaerobic metabolism—glycolysis) was observed in older spheroids (day 8 in culture) and those formed from more cells (Fig. 6). Interestingly, the irregular shape of WM115 allowed the fluorescent probe to penetrate deeper into the core, which was observed as a peak in fluorescence intensity in the region between 50 and 75 µm from the center of the spheroid (Fig. 6D). This glucose availability may facilitate anaerobic glucose metabolism (glycolysis) in the WM115 cell line and may be related to the higher o-Ps mean life time in these conditions.

We also performed additional measurements to compare the o-Ps lifetime in 2D and 3D cell cultures. The obtained results (Fig. 9) demonstrate that the positronium lifetime is smaller in 3D cell culture than in 2D cell culture55. This graph shows how the mean o-Ps lifetime is different between 2 and 3D cell cultures, which is due to their diversity in biological properties and metabolism. The WM155 and WM266-4 cell lines were characterized by different gene expression, e.g., those responsible for amino acid transport, and many others that are related to the formation of cell junctions, migration and mobility12,56.

Comparison of mean positronium lifetime in 2D and 3D cell culture of melanoma. The lifetime of o-Ps in WM266-4 2D cell culture (black square) according to Ref.55, the lifetime of o-Ps in WM266-4 2D cell culture (black triangle), and the lifetime of o-Ps in WM266-4 3D cell culture (black dots) show the o-Ps lifetime in spheroids of two different ages, 4 and 8 days after cell seeding.

This study was conducted to evaluate the lifetime distribution of ortho-positronium in WM266-4 and WM115 melanoma cell spheroids characterized by different degrees of malignancy. The measurements were performed at two different times (4 and 8 days) after seeding. During this time, both spheroid cell lines showed a reduction in o-Ps lifetime, which backs to the gradual cell proliferation in spheroids. More malignant WM266-4 spheroids showed a shorter o-Ps lifetime than WM115 cell line spheroids.

In summary, the results show that the lifetime of o-Ps is a useful parameter to differentiate between cancer cells with different biological characteristics. The higher the degree of malignancy and the rate of proliferation of neoplastic cells, the shorter the lifetime of ortho-positronium. This research paves the way for the development of positronium biomarkers using a cell spheroid model.

The datasets generated and analysed during the current study containing numerical data, mathematical models, and calculation results, including fitting graphs, are available in the Jagiellonian University Repository (RUJ) at the link https://ruj.uj.edu.pl/xmlui/handle/item/29740557.

Stępień, E., Karimi, H., Leszczyński, B. & Szczepanek, M. Melanoma spheroids as a model for cancer imaging study. Acta Phys. Pol. B. 51, 159–164 (2020).

Article ADS Google Scholar

Kapałczyńska, M. & Kolenda, T. 2D and 3D cell cultures—A comparison of different types of cancer cell cultures. Arch. Med. Sci. 14, 910–919 (2018).

PubMed Google Scholar

Kunz-Schughart, L. A., Kreutz, M. & Knuechel, R. Multicellular spheroids: A three-dimensional in vitro culture system to study tumour biology. Int. J. Exp. Pathol. 79, 1–23 (1998).

Article CAS PubMed PubMed Central Google Scholar

Vörsmann, H. et al. Development of a human three-dimensional organotypic skin-melanoma spheroid model for in vitro drug testing. Cell Death Dis. 4, e719 (2013).

Article PubMed PubMed Central Google Scholar

Lagies, S. et al. Cells grown in three-dimensional spheroids mirror in vivo metabolic response of epithelial cells. Commun. Biol. 3(1), 246. https://doi.org/10.1038/s42003-020-0973-6 (2020).

Article CAS PubMed PubMed Central Google Scholar

https://www.wcrf.org/cancer-trends/skin-cancer-statistics.

Casalou, C., Moreiras, H., Mayatra, J. M., Fabre, A. & Tobin, D. J. Loss of "epidermal melanin unit" integrity in human skin during melanoma-genesis. Front. Oncol. 27(12), 878336. https://doi.org/10.3389/fonc.2022.878336 (2022).

Article Google Scholar

Wróbel, S., Przybyło, M. & Stępień, E. The clinical trial landscape for melanoma therapies. J. Clin. Med. 8, 368–382 (2019).

Article PubMed PubMed Central Google Scholar

Karimkhani, C. et al. The global burden of melanoma: Results from the global burden of disease study 2015. Br. J. Dermatol. 177, 134–140 (2017).

Article CAS PubMed PubMed Central Google Scholar

Arnold, M. et al. Global burden of cutaneous melanoma in 2020 and projections to 2040. JAMA Dermatol. 158(5), 495–503. https://doi.org/10.1001/jamadermatol.2022.0160 (2020).

Article Google Scholar

Meghnani, V., Vetter, S. & Leclerc, E. RAGE overexpression confers a metastatic phenotype to the WM115 human primary melanoma cell line. Biochim. Biophys. Acta Mol. Basis Dis. 1842, 1017–1027 (2014).

Article CAS Google Scholar

Szczepanek, M., Panek, D., Przybyło, M., Moskal, P. & Stępień, E. Transcriptomic data analysis of melanocytes and melanoma cell lines of LAT transporter genes for precise medicine. Bio-Algorithms Med-Syst. 18(1), 144–150. https://doi.org/10.2478/bioal-2022-0086 (2022).

Article Google Scholar

Moskal, P., Jasińska, B., Stępień, E. & Bass, S. Positronium in medicine and biology. Nat. Rev. Phys. 1, 527–529 (2019).

Article Google Scholar

Bass, S. D., Mariazzi, S., Moskal, P., Stepien, E. Positronium physics and biomedical applications. Rev. Mod. Phys. (2023) https://arxiv.org/abs/2302.09246.

Moskal, P. et al. Feasibility study of the positronium imaging with the J-PET tomograph. Phys. Med. Biol. 64, 055017 (2019).

Article CAS PubMed Google Scholar

Jasińska, B. et al. Determination of the 3 γ fraction from positron annihilation in mesoporous materials for symmetry violation experiment with J-PET scanner. Acta Phys. Pol. B. 47, 453–460 (2016).

Article ADS Google Scholar

Stepanov, P. S. et al. Interaction of positronium with dissolved oxygen in liquids. Phys. Chem. Chem. Phys. 22, 5123–5131 (2020).

Article CAS PubMed Google Scholar

Shibuya, K., Saito, H., Nishikido, F., Takahashi, M. & Yamaya, T. Oxygen sensing ability of positronium atom for tumor hypoxia imaging. Commun. Phys. 3, 1–8 (2020).

Article Google Scholar

Moskal, P. et al. Performance assessment of the 2 γ positronium imaging with the total-body PET scanners. EJNMMI Phys. 7, 1–6 (2020).

Article Google Scholar

Axpe, E. et al. Detection of atomic scale changes in the free volume void size of three-dimensional colorectal cancer cell culture using positron annihilation lifetime spectroscopy. PLoS ONE 9, e83838 (2014).

Article ADS PubMed PubMed Central Google Scholar

Jean, Y. & Ache, H. J. Positronium reactions in micellar systems. J. Am. Chem. Soc. 99, 7504–7509 (1977).

Article CAS Google Scholar

Liu, G., Chen, H., Chakka, L., Gadzia, J. E. & Jean, Y. C. Applications of positron annihilation to dermatology and skin cancer. Phys. Status Solidi C. 4, 3912–3915 (2007).

Article ADS CAS Google Scholar

Jean, Y. C. et al. Applications of slow positrons to cancer research: Search for selectivity of positron annihilation to skin cancer. Appl. Surf. Sci. 252, 3166–3171 (2006).

Article ADS CAS Google Scholar

Jasińska, B. et al. Human tissues investigation using PALS technique. Acta Phys. Pol. B. 48, 1737–1748 (2017).

Article ADS Google Scholar

Moskal, P. et al. Developing a novel positronium biomarker for cardiac myxoma imaging. EJNMMI Phys. 10(1), 22. https://doi.org/10.1186/s40658-023-00543-w (2023).

Article PubMed PubMed Central Google Scholar

Moskal, P. et al. Positronium imaging with the novel multiphoton PET scanner. Sci. Adv. 7, eabh4394 (2021).

Article ADS CAS PubMed Google Scholar

Moskal, P. & Stępień, E. Positronium as a biomarker of hypoxia. Bio-Algoritms Med-Syst. 17, 314–319 (2021).

Google Scholar

Moskal, P. et al. Testing CPT symmetry in ortho-positronium decays with positronium annihilation tomography. Nat. Commun. 12, 5658. https://doi.org/10.1038/s41467-021-25905-9 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Moskal, P. & Stępień, E. Ł. Prospects and clinical perspectives of total-body PET imaging using plastic scintillators. PET Clin. 15, 439–452 (2020).

Article PubMed Google Scholar

Matulewicz, T. Radioactive nuclei for β+ γ PET and theranostics: Selected candidates. Bio-Algoritms Med-Syst. 17, 235–239 (2021).

Google Scholar

Choiński, J. & Łyczko, M. Prospects for the production of radioisotopes and radiobioconjugates for theranostics. Bio-Algoritms Med-Syst. 17, 241–257 (2021).

Google Scholar

Badawi, R. D. et al. First human imaging studies with the EXPLORER total-body PET scanner. J. Nucl. Med. 60, 299–303 (2019).

Article CAS PubMed PubMed Central Google Scholar

Vera, D. B. et al. First-in-human total-body PET imaging of HIV with 89Zr-VRC01 on the EXPLORER. J. Nucl. Med. 61, 545 (2020).

Google Scholar

Vandenberghe, S., Moskal, P. & Karp, J. S. State of the art in total body PET. EJNMMI Phys. 7, 1–33 (2020).

Article Google Scholar

Karp, J. S. et al. PennPET explorer: Design and preliminary performance of a whole-body imager. J. Nucl. Med. 61, 136–143 (2020).

Article CAS PubMed PubMed Central Google Scholar

Alavi, A., Werner, T., Stępień, E. & Moskal, P. Unparalleled and revolutionary impact of PET imaging on research and day to day practice of medicine. Bio-Algorithms Med-Syst. 17(4), 203–212. https://doi.org/10.1515/bams-2021-0186 (2021).

Article Google Scholar

Moskal, P. & Stępień, E. Ł. Perspectives on translation of positronium imaging into clinics. Front. Phys. 10, 969806 (2022).

Article Google Scholar

Karimi, H. et al. X-ray microtomography as a new approach for imaging and analysis of tumor spheroids. Micron 137, 102917 (2020).

Article PubMed Google Scholar

Wassmer, C. H. et al. Engineering of primary pancreatic islet cell spheroids for three-dimensional culture or transplantation: A methodological comparative study. Cell Transplant. 29, 963689720937292. https://doi.org/10.1177/0963689720937292 (2020).

Article PubMed Google Scholar

Dulski, K. PALS avalanche: A new PAL spectra analysis software. Acta Phys. Pol. A. 137, 167–170. https://doi.org/10.12693/APhysPolA.137.167 (2020).

Article ADS CAS Google Scholar

Dulski, K. et al. The J-PET detector—A tool for precision studies of ortho-positronium decays. Nucl. Instrum. Methods. 1008, 165452 (2021).

Article CAS Google Scholar

Sakai, Y. & Nakazawa, K. Technique for the control of spheroid diameter using microfabricated chips. Acta Biomater. 3, 1033–1040. https://doi.org/10.1016/j.actbio.2007.06.004 (2007).

Article CAS PubMed Google Scholar

Song, Y. S. et al. Volume and mass doubling times of persistent pulmonary subsolid nodules detected in patients without known malignancy. Radiology 273, 276–284 (2014).

Article PubMed Google Scholar

Mehrara, E., Forssell-Aronsson, E., Ahlman, H. & Bernhardt, P. Specific growth rate versus doubling time for quantitative characterization of tumor growth rate. Cancer Res. 67, 3970–3975 (2007).

Article CAS PubMed Google Scholar

Schwartz, M. A biomathematical approach to clinical tumor growth. Cancer 14, 1272–1294 (1961).

3.0.CO;2-H" data-track-action="article reference" href="https://doi.org/10.1002%2F1097-0142%28196111%2F12%2914%3A6%3C1272%3A%3AAID-CNCR2820140618%3E3.0.CO%3B2-H" aria-label="Article reference 45" data-doi="10.1002/1097-0142(196111/12)14:63.0.CO;2-H">Article CAS PubMed Google Scholar

Moskal, P. Positronium imaging. In 2019 IEEE Nucl. Sci. Symp. Med. Imaging Conf. Proc. (NSS/MIC) 1–3 (2019).

Costa, E. et al. 3D tumor spheroids: An overview on the tools and techniques used for their analysis. Biotechnol. Adv. 34, 1427–1441 (2016).

Article PubMed Google Scholar

Pampaloni, F., Ansari, N. & Stelzer, E. High resolution deep imaging of live cellular spheroids with light-sheet-based fluorescence microscopy. Cell Tissue Res. 352, 161–177 (2013).

Article PubMed Google Scholar

Milotti, E., Vyshemirsky, V., Sega, M., Stella, S. & Chignola, R. Metabolic scaling in solid tumours. Sci. Rep. 3, 1938. https://doi.org/10.1038/srep01938 (2013).

Article ADS CAS PubMed PubMed Central Google Scholar

Beaumont, K. A., Anfosso, A., Ahmed, F., Weninger, W. & Haass, N. K. Imaging-and flow cytometry-based analysis of cell position and the cell cycle in 3D melanoma spheroids. J. Vis. Exp. JOVE. 106, e53486 (2015).

Google Scholar

Sane, P. et al. Temperature-induced structural transition in-situ in porcine lens—Changes observed in void size distribution. Biochim. Biophys. Acta Biomembr. 1798, 958–965 (2010).

Article CAS Google Scholar

Elias, M., Al-Mashhadani, A. & Al-Shiebani, Z. Temperature dependence of microstructure of biological tissues probed via the positronium method. Pure Sci. 28, 240–244 (2001).

CAS Google Scholar

Chandrashekara, M. N. & Ranganathaiah, C. Free volume size distribution in some natural polymers. AIP Conf. Proc. 1349, 228–229 (2011).

Article ADS CAS Google Scholar

Chandrashekara, M. N. & Ranganathaiah, C. Diffusion of permanent liquid dye molecules in human hair investigated by positron lifetime spectroscopy. Colloids Surf. B. 69, 129–134 (2009).

Article CAS Google Scholar

Kubicz, E. Biomedical applications of Positron Annihilation Lifetime Spectroscopy: Nano structural characterization of normal and cancer cells and tissues. PhD thesis, Jagiellonian University (2020).

Kubicz, E. & Stępień, E. Results transcriptomes sequencing dataset for human melanoma and melanocytes cell cultures. Jagiellonian University Repository. Raw data (2021) https://ruj.uj.edu.pl/xmlui/handle/item/269420

Stępień, E., Karimi, H. & Moskal, P. 3D melanoma spheroid model for the development of positronium biomarker. Jagiellonian University Repository. Raw data. https://ruj.uj.edu.pl/xmlui/handle/item/297405.

Download references

This work was supported by the Foundation for Polish Science (FNP) through the TEAM POIR.04.04.00-00-4204/17 program, by the National Science Centre of Poland through grants no. 2019/33/B/NZ3/01004, 2021/42/A/ST2/00423 and 2021/43/B/ST2/02150, the Jagiellonian University via projects CRP/0641.221.2020, SciMat and qLife Priority Research Areas budget under the program Excellence Initiative—Research University, and the DSC grant no. N17/MNS/000023 Mrs. Hanieh Karimi, for her PhD research.

Department of Medical Physics, M. Smoluchowski Institute of Physics, Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Łojasiewicza 11 Street, 30-348, Kraków, Poland

Hanieh Karimi & Ewa Ł. Stępień

Department of Biochemistry, University of Missouri, Columbia, USA

Hanieh Karimi

Department of Experimental Particle Physics and Applications, M. Smoluchowski Institute of Physics, Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Kraków, Poland

Paweł Moskal

Center for Theranostics, Jagiellonian University, Kraków, Poland

Paweł Moskal & Ewa Ł. Stępień

Faculty of Chemistry, Jagiellonian University, Kraków, Poland

Agata Żak

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H.K.: Methodology, investigation, writing—original draft, visualization, review & editing,P.M.: Conceptualization, methodology, software, formal analysis, resources, data curation, writing—original draft, review and editing, visualization, supervision, project administration, funding acquisition,A.Ż.: Methodology, investigation, E.Ł.S.: Conceptualization, methodology, investigation, resources, data curation, writing—review & editing, supervision, funding acquisition.

Correspondence to Ewa Ł. Stępień.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Karimi, H., Moskal, P., Żak, A. et al. 3D melanoma spheroid model for the development of positronium biomarkers. Sci Rep 13, 7648 (2023). https://doi.org/10.1038/s41598-023-34571-4

Download citation

Received: 19 June 2022

Accepted: 03 May 2023

Published: 11 May 2023

DOI: https://doi.org/10.1038/s41598-023-34571-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.