Preprint
Article

Galactosylceramide Upregulates the Expression of the BCL2 Gene and Downregulates the Expression of TNFRSF1B and TNFRSF9 Genes, Acting as an Ani-apoptotic Molecule in Breast Cancer Cells

Altmetrics

Downloads

83

Views

48

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

29 November 2023

Posted:

05 December 2023

You are already at the latest version

Alerts
Abstract
Using a breast cancer (BC) cellular model: MDA-MB-231 cells with blocked synthesis of galacto-sylceramide (GalCer; MDA.Δ.UGT8.4) and MCF7 and T47D cells that over-produced GalCer (MCF7.UGT8 and T47D.UGT8), we found that GalCer increased resistance to doxorubicin, paclitaxel, and cisplatin serving as an anti-apoptotic molecule. These results were supported by a study using murine mammary carcinoma 4T1 cells over-producing GalCer, which were trans-planted into Balb/c mice and were more resistant to doxorubicin therapy than the control. Sub-sequent studies revealed that GalCer specifically downregulated the level of pro-apoptotic TNFRSF1B and TNFRSF9 genes and increased the level of the anti-apoptotic BCL2 gene expres-sion, suggesting that GalCer regulates their expression at the level of transcription. Consistent with this hypothesis, MDA-MB-231 and MCF7.UGT8 cells with high GalCer content showed lower activity of TNFRSF1B and TNFRSF9 promoters than cells lacking GalCer. In contrast, the activity of the BCL2 promoter in MCF7.UGT8 was higher than in MCF7 cells without GalCer. However, no difference in BCL2 promoter activity was observed between MDA-MB-231 cells with high GalCer level and GalCer-negative MDA.Δ.UGT8.4 cells. Instead, we found that GalCer increased the stability of Bcl-2 mRNA. As a candidate protein that simultaneously increases the expression of TNFRSF1B and TNFRSF9 genes and decreases the expression of BCL2 gene and stability of Bcl-2 transcripts is probably p53, its expression was analysed in BC cells containing various amounts of GalCer. In agreement with this hypothesis, BC cells with high GalCer con-tent are characterized by significantly lower expression of p53, and inhibition of p53 using siR-NA resulted in decreased expression of TNFRS1B and TNFRS9 genes and increased expression of BCL2 gene.
Keywords: 
Subject: Biology and Life Sciences  -   Biochemistry and Molecular Biology

1. Introduction

Breast cancer (BC) is one of the greatest oncological challenges in developed countries. High morbidity, together with a high mortality rate, makes BC exceptionally important medically. Drug resistance (DR) is one of the major concerns, particularly in patients with advanced cancer. Different mechanisms of DR are responsible for tumor treatment failures, including decreased intracellular drug concentrations, disturbances in cell cycle, resistance to apoptosis, DNA repair deficiencies, activation of intracellular pathways related to cancer progression, epigenetic alterations, and changes in drug targets [1]. One class of cellular compounds involved in the DR of cancer cells are glycosphingolipids (GSLs) [2]. The involvement of GSLs in DR was first described by Lavie et al. [3], who showed that DR cancer cells were characterized by the accumulation of glucosylceramide (GlcCer). Initially, it was proposed that increased synthesis of GlcCer highly decreases the intracellular content of pro-apoptotic ceramides, making cancer cells more resistant to drug-induced apoptosis [4]. However, further studies showed that overexpression of GCS and accumulation of GlcCer correlated with high expression of P-glycoprotein (Pgp, MDR1), suggesting that GlcCer positively regulates the expression of MDR1 gene [5]. A subsequent study revealed that increased levels of globo-series GSLs (globotriaosylceramide and globotetraosylceramide), but not GlcCer by itself, are responsible for up-regulation of Pgp in cancer cells [6,7,8]. In addition to globo-series GSLs, a high content of GM3 ganglioside was shown to protect murine lung carcinoma cells from cytotoxic effects of etoposide and doxorubicin [9], and small-cell lung cancer cells from cytotoxic effects of DOX and cisplatin [10]. However, in chronic myeloid leukemic K562 cells, increased GM3 content decreased resistance to etoposide and staurosporine [11]. In the case of GM3, increased or decreased resistance to chemotherapy has been linked to differences in the resistance of cancer cells to apoptosis induced by anti-cancer drugs.
Galactosylceramide (GalCer) is a GSL synthesised by galactosylceramide synthase (UGT8, C.E. 2.4.2.62) [12] encoded by the UGT8 gene [13], and is a major myelin-stabilizing component [14]. However, it was recently shown that BC cells expressing high level of UGT8 and accumulating GalCer, form tumors and metastatic colonies in the lungs more efficiently than BC cells with suppressed UGT8 expression [15]. In accordance with these results, tumor specimens with elevated GalCer content were characterized by higher Ki-67 proliferation index and a reduction in the number of apoptotic cells. Therefore, it was hypothesized that GalCer serves as an anti-apoptotic molecule, and its presence increases the ability of cancer cells to survive in the unfavourable tumor microenvironment [15,16]. In addition, GalCer was found to increase resistance of BC cells to doxorubicin, suggesting that GalCer can be an important factor in the BC cells resistance to chemotherapy DR. This hypothesis is endorsed by findings that GalCer content is increased in DR colon and ovarian cancer cells [17,18], and that Krabbe cells synthesizing large amounts of GalCer are more resistant to apoptosis induced by daunorubicin and cytosine arabinoside than Gaucher cells, which contain low levels of this GSL [19]. Furthermore, inhibition of GlCer in U937 and HL60 cells with such inhibitors as D,L-threo-1-phenyl-2-palmitoylamino-3-morpholino-1-propanol (PDMP) or 1-phenyl-2-palmitoylamino-3-morpholino-1-propanol (PPMP), which was accompanied by increased levels of GalCer, protects them from daunorubicin-induced apoptosis [19]. In addition, such cells with high content of exogenously added GalCer became more resistant to daunorubicin-induced apoptosis. Accumulation of GalCer was observed in Madin-Darby canine kidney cells exposed to hyperosmotic and heat stress [20,21]. Based on this data and generally accepted idea that ceramide is one of the important pro-apoptotic molecules [22,23], we proposed that the anti-apoptotic activity of GalCer can be liked to decreased levels of ceramide due to increased synthesis of GalCer [15]. However, it was found that BC cells with various amounts of GalCer contain essentially the same amounts of ceramide [16], suggesting that glycosylation of ceramide is not a major mechanism that increases resistance of tumor cells to drug-induced apoptosis. Consequently, a combination of molecular and pharmacological techniques were used in this research to examine the function of GalCer in BC cell resistance to anti-cancer agents and elucidate the mechanisms of its anti-apoptotic activity.

2. Materials and Methods

2.1. Cell lines and culture conditions

Human BC MDA-MB-231 and T47D cells were obtained from American Type Culture Collection (ATCC, VA, USA). The BC MCF7 and murine 4T1 mammary cancer cell lines were kindly provided by the Cell Line Collection of the Hirszfeld Institute of Immunology and Experimental Therapy (Wroclaw, Poland). MCF7 cell line were authenticated by the ATCC Cell Line Authentication Service. T47D.UGT8 cells overexpressing UGT8 and control T47D.C cells derived from T47D cells have previously been described by Suchanski et al. [16]. MDA-MB-231 and MCF7 cells were cultured in α-MEM supplemented with 10% fetal bovine serum (FBS, Cytogen, Poland), 2 mM L-glutamine, and antibiotics. T47D and 4T1 cells were cultured in RPMI-1640 Medium with 10% FBS (Cytogen), 2 mM L-glutamine and antibiotics. For T47D cells, the culture medium was supplemented with bovine insulin (0.2 U/mL) (Merck, Germany).

2.2. Vectors construction, virus production and transductions

Using PCR, human UGT8 cDNA was generated from the MDA-MB-231 cDNA library as a template for primers: forXhoI-UGT8 and revBamHI-UGT8 (Supplementary Table S1), and mouse galactosylceramide synthase (UGT8a) cDNA was amplified from the pCMV6-Entry vector containing UGT8a cDNA (MR208632, 1626 bp; OriGene, MD, USA) with forEcoRI-UGT8a and revMluI-UGT8a primers (Supplementary Table S1). To generate UGT8 and UGT8a-expressing vectors, the resulting inserts were cloned into a pLVX-Puro Vector (Clontech, CA, USA) as previously described [16] and named as pLVX-UGT8-Puro and pLVX-UGT8a-Puro, respectively. For knockout of the UGT8 gene, a UGT8 CRISPR/Cas9 vector encoding the Cas9 nuclease and UGT8-specific 20 nucleotide guide RNA (gRNA-TGTGATAGCTCATCTTTTAG) was purchased from Applied Biological Materials Inc. (Canada). For the generation of lentivirus, LentiX 293T cells (Clontech) were transfected at 50 – 60% confluence with Lenti-X Packaging Single Shots (Takara, Japan). The production of virus particles and cell transductions have been previously described [15]. UGT8 knockout cells were selected using puromycin (1 μg/mL) (Thermo Fisher Scientific, UK). Antibiotic-resistant cells were subcloned using the limiting dilution technique. The mutation introduced by CRISPR-Cas9 editing was identified by sequencing genomic DNA with the primers: F1-exon2/UGT8 and R1-exon2/UGT8 (Supplementary Table S1).

2.3. Quantitative PCR (qPCR) assay and apoptosis gene expression profiling with RT-qPCR array

RNA was purified from BC cells using an RNeasy Mini Kit (Qiagen, Germany). To synthesize cDNA, the SuperScipt RT Kit (Thermo Fisher Scientific) was used. The relative amounts of mRNAs were quantified by qPCR with iQ SYBR Green Supermix (Bio-Rad, CA, USA) using an iQ5 Optical System (Bio-Rad) and ACTB as the reference gene. The primers for qPCR are shown in Supplementary Table S1. The gene expression profiles of 46 apoptosis-related genes were determined by RT-qPCR (Supplementary Table S2) using the EvaGreen dye-based detection system (Biotium, CA, USA). Two reference genes, ACTB and GAPDH, were used to normalise RNA levels. PCR amplification consisted of an initial denaturation for 3 min at 95ºC, followed by 35 cycles of 20 s denaturation at 95ºC, annealing for 20 s at 5ºC and extension for 20 s at 72ºC. The relative changes in gene expression were calculated using the ΔΔCt (cycle threshold) method. Fold change values were calculated using the 2-ΔΔCt formula.

2.4. Western blotting

Western blotting was performed as previously described [16]. The following rabbit polyclonal antibodies were used to detect specific proteins: anti-UGT8 (Biorbyt, UK), anti-Bcl-2 (Abcam, UK), anti-TNFRSF1B (Merck), anti-TNFRSF9 (Biorbyt), anti-p53 (Cell Signaling, MA, USA), anti-CREB (Cell Signaling, MA, USA), anti NF-κB p65 (Abcam, UK) and mouse monoclonal anti-GAPDH antibody (Novus Biologicals, UK).

2.5. Purification of neutral GSLs

GSLs were extracted from 108–109 cells by the chloroform–methanol extraction method and neutral GSLs were purified as previously described .

2.6. Thin-layer chromatogram (TLC)-binding assay

High-performance-thin layer chromatography (HP-TLC) was performed on silica gel 60 HP-TLC plates (Merck). Neutral GSLs were separated using a 2-isopropanol:methyl acetate:15 M ammonium hydroxide:water solvent system (75:10:5:15, v/v/v/v). GSLs were visualised by spraying the plate with primuline reagent (0.05% primuline in acetone/water, 4:1 v/v) and heated for 1 min at 120°C. GalCer was detected by a TLC-binding assay as described previously [15], using rabbit polyclonal anti-GalCer antibodies (Merck).

2.7. Cell survival assay

The viability of the cells was assessed by MTT assay. Briefly, 5 × 103 cells/well were seeded in a 96-well plate (Nunc, NY USA) and incubated with increasing concentrations of doxorubicin hydrochloride (DOX, Pfizer, NY, USA), paclitaxel (PAX, Ebewe Pharma, Austria) or cisplatin cis-diamminedichloroplatinum (II) (CDDP, Merck) for 48 h. Then, the medium was discarded, and the cells were incubated with MTT solution (0.5 mg/mL, Merck) for 4 h at 37°C. After removing the MTT solution, 200 μL of DMSO (Chempur, Poland) was added for 15 min at 22°C and the optical density was measured at 570 nm.

2.8. Apoptotic assay

The cells (5 × 105) were seeded in multi-well plates (Nunc). The following day, the cells were treated with DOX with or without the selective inhibitor ABT-199 of Bcl-2 protein (Cayman Chemical, MI, USA) or TNFα (Thermo Fisher Scientific) for 48 h. The percentages of apoptotic cells were determined using the Dead Cell Apoptosis Kit (Thermo Fisher Scientific). Fluorescence was measured on a BD FACS Lyric (Becton-Dickinson, NJ, USA) and data were processed using the BD FACSuite™ Software (Becton-Dickinson). The cells in early apoptosis were stained with Annexin V APC , whereas cells in late apoptosis were stained with Annexin V APC and SYTOX Green.

2.9. In-vivo tumor growth assay

The animal experiments were approved by the Second Local Ethics Committee for Animal Experimentation (033/2019/P1,Wroclaw, Poland). Six-week old female BALB/c mice were obtained from the Mossakowski Medical Research Institute, Polish Academy of Sciences (Warsaw, Poland) and kept under specific pathogen-free conditions at room temperature. Suspensions of 4T1 cells (2 x 104 cells/100 µL PBS) were mixed with the equal volume of ice-cold BD Matrigel Matrix High Concentration (Becton-Dickinson) and injected subcutaneously (s.c.). Growth of tumors was monitored once a week by measuring tumour diameter using a calliper. The tumor volume (TV) was calculated as TV (mm3) = (d2x D)/2, where d is the shortest diameter and D is the longest diameter. For the therapeutic experiments, mice were injected intravenously (i.v.) twice a week with DOX (1.5 mg/kg body weight). Mice were euthanised by cervical dislocation following light anaesthesia by isoflurane (Forane, Abott Laboratories, IL, USA) inhalation. Tissue samples were collected in 10% buffered formalin and subjected to histological analyses.

2.10. Cloning of human TNFRSF1B, TNFRSF9 and BCL2 gene promoters and determination of their activity using luciferase reporter vectors

The Genomatix MatInspector software (http://www. genomatix.de) was applied to identify the nucleotide sequences of TNFRSF1B, TNFRSF9 and BCL2 promoters, which were generated by PCR using genomic DNA purified from MDA-MB-231 cells as a template and primers listed in Supplementary Table S1. They were amplified as follows: 3 min of initial denaturation at 95°C, followed by 30 cycles of 20 s denaturation at 95ºC, annealing at 55°C for 20 s, and 1 min of final extension at 72ºC. The resulting PCR products containing additional sequences corresponding to MluI and HindIII restriction sites were cloned into pGL3-Basic luciferase vector (Promega, USA). The obtained constructs were named pGL3-TNFRSF1B/LUC, pGL3-TNFRSF9/LUC, and pGL3-BCL2//LUC for the TNFRSF1B, TNFRSF9, and BCL2 promoters, respectively.
To determine the activity of the promoters, pGL3-TNFRSF1B/LUC, pGL3-TNFRSF9/LUC, and pGL3-BCL2/LUC constructs (2 μg) were co-transfected with control Renilla luciferase-expressing pRL-TK vector (2 μg) using polyethylenimine (Merck) into BC cells seeded into multi-well plates and grown to 70% confluence. After 48 h, to measure promoter activity, cells were analyzed using the Dual-Luciferase Reporter Assay System (Promega) following recommendations provided by the manufacturer.

2.10. mRNA stability assay

Cells (5 x 105) were treated with 2 μg/mL of actinomycin D (Thermo Fisher Scientific) in DMSO for 2, 4 and 6 h at 37°C. At the indicated time points, RNA was extracted using Fenozol reagent (A&A Biotechnology, Poland) and the level of Bcl-2 transcripts was analysed using RT-qPCR with iQ SYBR Green Supermix (Bio-Rad) as described above in qPCR assay section.

2.11. siRNA transfections

Transfection with mixture of p53 siRNAs (Santa Cruz Biotechnology, USA) was performed using Lipofectamine RNAiMAX transfection reagent (Thermo Fisher Scientific). Briefly, cells (2.5x105) were seeded in a 6-well plates in complete α-MEM and at 60% confluence incubated with siRNAs for 72 h. The transfection reagent was prepared according to RNAiMAX Transfection Procedure (Thermo Fisher Scientific). Cells were harvested by trypsinization and subjected to further experiments.

2.12. Statistical analysis

All statistical analyses were performed using Prism 5.0 (GraphPad, CA, USA). The two-way Anova Dunnetts post-test was used for statistical analysis of in vivo tumor growth experiments. The Mann–Whitney test was used to compare the groups of data that did not meet the assumptions of the parametric test. In all analyses, the results were considered statistically significant at p < 0.05.

3. Results

3.1. DR of BC cells depends on expression of UGT8 and synthesis of GalCer

We previously showed that GalCer is responsible for resistance of BC cells to DOX [15]. Based on these results, this study has been carried out to further assess the GalCer function in BC cell resistance to commonly used chemotherapeutics, using cells with different expression levels of this GSL. A model, representing a loss-of-function phenotype, was obtained by CRISPR-edited knockout of the UGT8 gene in MDA-MB-231 cells (Supplementary Figure S1). The two UGT8-negative, puromycin-resistant clones, named MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4, were chosen for further studies (Figure 1A). In addition, a control cell line MDA.Δ.C was created by transducing MDA-MB-231 cells with the CRISPR/Cas9 plasmid encoding a non-specific 20 nt gRNA. The model, representing the gain-of-function phenotype, was obtained overexpressing UGT8 in MCF7 cells (Figure 1B). The cells that accumulated GalCer were named MCF7.UGT8. The control MCF7.C cell line was generated by transducing MCF7 cells with an empty vector. A second gain-of-function cellular model used in this study was represented by T47D.UGT8 cells [16]. All cells used in these studies had the same morphology and proliferative potential as parental cells (Supplementary Figure S2).
To assess, whether the presence of GalCer affects the viability of BC cells treated with chemotherapeutic agents, cells were incubated with DOX, PAX, and CDDP at different concentrations. Using loss-of-function and gain-of-function cellular models, it was shown that all analysed GalCer-rich BC cells were more resistant to cytotoxic effects induced by DOX, PTX and CDDP than BC cells unable to synthesise this GSL (Figure 1C). Notably, the sensitivity of the control MDA.Δ.C, MCF7.C and T47D.C cells to anti-cancer drugs was the same as the respective parental wild-type BC cells.

3.2. Resistance of murine 4T1 mammary tumours to DOX treatment is dependent on the presence of GalCer

To further prove that GalCer is directly involved in increased resistance of BC cells to anticancer drugs, an in vivo study using murine 4T1.UGT8a mammary carcinoma cells overexpressing mouse UGT8a and GalCer was performed (Figure 2A). The 4T1.UGT8a cells were obtained by transducing 4T1 cells with a vector encoding mouse UGT8a cDNA. Similar to human BC cells, 4T1.UGT8a cells were more resistant to DOX as their viability was higher than that of control 4T1.C cells carrying the “empty” vector (Figure 2B). Initially, we assessed whether the presence of GalCer affected the tumorigenicity of murine mammary carcinoma cells. After s.c. transplantation into BALB/c mice, on the day 22, the mean tumor volumes (357±63 mm3) formed by 4T1.UGT8a cells were similar to tumor volumes (373 ±68 mm3) formed by control 4T1.C cells (Figure 2C). To further evaluate, whether the sensitivity of these tumors to DOX treatment can be affected by the accumulation of GalCer, mice were injected i.v. on days 7, 10, 13, and 16 after mammary carcinoma cell transplantation with DOX (Figure 2C).
Tumor regression was observed in mice transplanted with 4T1.UGT8a cells with high GalCer and GalCer-negative 4T1.C cells compared to mice treated with placebo. However, only in case of GalCer-negative 4T1.C tumors, the reduction in tumor volume was statistically significant in comparison to that in mice subjected to treatment with placebo. Therefore, these data confirm our proposal that GalCer is involved in DR of BC tumors.

3.3. GalCer is responsible for changes in the expression of specific apoptotic genes

To further elucidate the mechanism by which GalCer mediates cytoprotective effects on BC cells, we compared the expression of 46 apoptotic genes in BC cells with high and low GalCer content using qPCR (Supplementary Figures S3 and S4).
All analysed BC cells that had a high content of GalCer (MDA.Δ.C, MCF7.UGT8, and T47D.UGT8) were characterised by considerably increased levels of the anti-apoptotic BCL-2 gene and substantially decreased levels of pro-apoptotic TNFRSF1B and TNFRSF9 genes in comparison to BC cells with low content of GalCer (MDA.Δ.UGT8.1, MDA.Δ.UGT8.4, MCF7.C and T47D.C). Consistent differences between BC cells with high and low content were found for these three genes, and their expression in loss-of-function and gain-of-function cellular models was further validated at the mRNA and protein levels (Figure 3A,B). As expected, the expression of Bcl-2 protein was markedly increased in MDA.C and MCF7.UGT8 and T47D.UGT8 cells compared to MDA.Δ.UGT8.1, MDA.Δ.UGT8.4 and MCF7.C and T47D.C cells. In contrast, the levels of TNFRSF1B and TNFRSF9 were higher in MDA.Δ.UGT8.1, MDA.Δ.UGT8.4, MCF7.C and T47D.C cells than in MDA.Δ.C and MCF7.UGT8 and T47D.UGT8 cells.
To prove that the changes in Bcl-2 expression observed in BC cells with various contents of GalCer directly affects their sensitivity to anti-cancer drugs, MDA.C cells with high levels of GalCer and high Bcl-2 expression were treated with ABT-199, which is a selective Bcl-2 inhibitor [25], and subsequently incubated with DOX. Using flow cytometry, it was found that treatment of MDA.C cells with ABT-199 increased their sensitivity to DOX-induced apoptosis compared to MDA.C cells incubated with DOX only (Figure 4A, Supplementary Figure S5A). However, the sensitivity of the MDA.Δ.UGT8.4 cells with a lack of GalCer and low levels of Bcl-2 to DOX-induced apoptosis was not affected by treatment with the Bcl-2 inhibitor. These results were verified using a gain-of-function cellular model. When MCF7.UGT8 cells, also with high GalCer content and high Bcl-2 expression, were treated with ABT-199, their sensitivity to DOX-induced apoptosis was markedly increased compared to MCF7.UGT8 cells incubated with DOX only (Figure 4B, Supplementary Figure S5B). In contrast, DOX-induced apoptosis was not affected by the Bcl-2 inhibitor in GalCer-negative MCF7.C cells with low Bcl-2 levels.
The effects of TNFRSF1B and TNFRSF9 on BC cell apoptosis were assessed by treating cells with human recombinant TNFα and flow cytometry. MDA.Δ.UGT8.4 cells with UGT8 gene knockout increased the expression of TNFRSF1B and TNFRSF9 and were more prone to apoptosis triggered by TNFα than control MDA.Δ.C cells (Figure 4C, Supplementary Figure S6A). However, MCF7.UGT8 cells overexpressing GalCer with low levels of TNFRSF1B and TNFRSF9 expression, that were incubated with TNFα, the percentage of apoptotic cells was considerably lower than in control MCF7.C cells (Figure 4D, Supplementary Figure S6B).
In summary, we found that GalCer-driven changes in the expression of Bcl-2, TNFRSF1B, and TNFRSF9 were responsible for the sensitivity/resistance of BC cells to DOX-induced apoptosis.

3.4. GalCer regulates the expression of TNFRSF1B, TNFRSF9, and Bcl-2 on the level of transcription and mRNA stability

Changes in the expression of TNFRSF1B, TNFRSF9, and Bcl-2 mRNAs found in BC cells with different GalCer levels suggested that it may influence the transcription of these genes. Therefore, to address this hypothesis, the TNFRSF1B, TNFRSF9 and BCL2 promoters were cloned and sequenced (Supplementary Figure S7) and dual-luciferase reporter assay system was used to analyse their activities in BC cells.
The activities of the TNFRSF1B and TNFRSF9 promoters were found to be higher in MDA.Δ.UGT8.4 and MCF7.C cells lacking GalCer than in MDA.Δ.C and MCF7.UGT8 with high GalCer levels (Figure 5A,B). In contrast, in the case of BCL2, higher promoter activity was observed in cells with a high GalCer level (Figure 5C). However, no differences in BCL2 promoter activity among MDA.Δ.C and MDA.Δ.UGT.4 were observed. Based on this result, we hypothesized that in MDA-MB-231 cells, Bcl-2 expression is regulated at the post-transcriptional level as shown in HL60, Jurkat or HEK 293 cells [26]. One of these mechanisms entails the control of mRNA stability [27]. Therefore, we analyzed the Bcl-2 transcript levels in MDA.Δ.C and MDA.Δ.UGT.4 cells treated with actinomycin D and found that the stability of Bcl-2 mRNA was significantly higher than that of MDA.Δ.UGT.4 cells with blocked GalCer synthesis (Figure 5D).
Based on these results, an in silico analysis was performed to identify the binding sites for potential transcription factors affecting the activities of TNFRSF1B, TNFRSF9, and BCL2 promoters in BC cells. Using the Genomatix MatInspector software, numerous transcriptional elements were identified in each of the analysed promoter sequences. However, only the binding sites for the transcription factors: p53, NF-ĸB1 snd CREB1 were present in all three promoters (Supplementary Figure S7), suggesting that they may regulate TNFRSF1B, TNFRSF9 and BCL2 gene expression. Therefore, we analyzed the expression of CREB1, NF-ĸB1 and p53 in BC cells with different levels of GalCer by Western blotting. MDA.Δ.C cells with a high level of GalCer and high expression level of Bcl-2 expressed significantly less p53 than MDA.Δ.UGT.4 cells lacking GalCer, and with low expression levels of Bcl-2 (Figure 6A). In contrast, GalCer-negative MCF7.C cells with low levels of Bcl-2 are characterised by considerably higher expression of p53 than MCF7.UGT8 cells overexpressing GalCer. On the other hand, no differences in the expression of CREB1 and NF-ĸB1 were found between MDA-MB-231 and MCF7 cells with different expression of GalCer (Supplementary Figure S8). Furthermore, CREB1 and NF-ĸB1 upregulate the expression of TNFR family members as well as a Bcl-2 [28,29,30,31,32]. In contrast, p53 upregulates TNFRSF1B and TNFRSF9 genes, but downregulates the expression of BCL2 gene [33]. Therefore, only p53 meets the conditions for regulatory factor affecting simultaneously expression of TNFRSF1B, TNFRSF9 and BCL2 genes, but in opposite way. To prove that changes in the expression of TNFRSF1B, TNFRSF9 genes are directly related to changes in the expression of p53, MDA.Δ.UGT.4 and MCF7.C cells were subjected to treatment with p53 siRNAs.
Inhibition of p53 expression in MDA.Δ.UGT.4 as well as MCF7.C cells resulted in highly decreased expression of TNFRSF1B, TNFRSF9 genes (Figure 6B). On the other hand, treatment of MCF7.C cells with p53 siRNAs resulted in increased expression of BCL2 gene, and lack of differences in expression of BCL2 gene was found in p53 siRNA-treated MDA.Δ.UGT.4 cells. To further confirm that p53 directly regulates the expression of TNFRSF1B, TNFRSF9 and BCL2 genes at the transcriptional level, we analysed promoter activities in BC cells after inhibition of p53 expression using siRNA. In agreement with previous findings (see Figure 5), decreased p53 levels resulted in highly decreased activity of TNFRSF1B and TNFRSF9 promoters and increased activity of BCL2 promoter in MCF7.C cells and also decreased activity of TNFRSF1B and TNFRSF9 promoters in MDA.Δ.UGT.4 cells (Figure 6C). Also, in line with previous observations, inhibition of p53 expression increased the stability of the Bcl-2 transcript (Figure 6D).

4. Discussion

We previously showed that BC MDA-MB-231 cells, which synthesize large amounts of GalCer become more resistant to DOX [15]. In the present study, using three human BC cellular models with loss-of-function (MDA-MB-231 cells) or gain-of-function (MCF7 and T47D cells) phenotypes, we extended this finding and proved that BC cells are also more resistant to other anti-cancer drugs such as PTX and CDDP. These in vitro results are strongly supported by an in vivo study using mouse model. Interestingly, not only human, but also murine mammary carcinoma 4T1 cells overexpressing murine UGT8a and accumulating GalCer become more resistant to DOX. When such cells were injected into BALB/c mice, they showed increased resistance to DOX therapy than the control 4T1 cells. Taken together, these data strongly suggest that GalCer is a marker of DR in BC cells and a new target for breast cancer treatment. In agreement with the general view that the majority of anticancer drugs induce apoptosis in cancer cells [34], we also showed that GalCer serves as an anti-apoptotic molecule, increasing the resistance of BC cells to apoptosis induced by DOX [15]. Therefore, the present study aimed to further define the mechanisms by which GalCer increases the resistance of BC cells to drug-induced apoptosis. We found that GalCer specifically downregulated the expression of TNFRSF1B and TNFRSF9 genes and increased BCL-2 gene expression at the mRNA and protein levels, proving that changes in their expression directly affect the apoptotic properties of BC cells. Therefore, this results showed for the first time that the same species of GSL could produce simultaneous but opposite effects on the extrinsic and intrinsic (mitochondrial) apoptotic pathways. In contrast to GalCer, other GSLs involved in the regulation of apoptotic gene expression belong to the ganglio-series of GSLs and act mainly as pro-apoptotic molecules. One of these gangliosides, GD3, is involved in suppression of the BCL-2 gene expression at the transcriptional level via dephosphorylation of AKT and CREB [35]. The forced expression of Bcl-2 protected from GD3-induced apoptosis in T cell lymphoma CEM cells and mouse oligodendrocytes, suggesting that relevant GD3 targets are under Bcl-2 control [36], which we observed in case of GalCer, but with an opposite effect. However, GD3, in contrast to GalCer, activates the mitochondrial pathway at the level of the permeability transition pore complex causing cytochrome c release and caspase activation [37]. In addition, GD3 acts as an pro-apoptotic agents by preventing (through a currently unknown mechanism) the nuclear localisation of active NF-κB that suppresses the NF-κB-dependent survival pathway [38]. Sequestration of NF-κB in the cytoplasm protects the activation of genes dependent on NF-κB, sensitises hepatocytes against apoptotic stimuli such as TNFα, and sensitises human hepatoma cells to daunorubicin and/or ionising radiation [39]. In contrast, overexpression of ganglioside sialidase (Neu3) in colon carcinoma cells was associated with increased Bcl-2 and decreased caspase expression, which is another example of apoptosis suppression associated with GSLs metabolism [40]. Increased resistance to drug-induced apoptosis associated with Bcl-2 expression was observed in mouse 3LL Lewis lung cancer cells overexpressing GM3 ganglioside compared with GM3-negative parental cells [9]. In this case, upregulation of Bcl-2 was linked to the post-translational protein modifications.
To further prove that GalCer affects the expression of specific apoptotic genes at the transcriptional level, the promoters of TNFRSF1B, TNFRSF9, and BCL2 genes were cloned and their activities were determined using MDA-MB-231 and MCF7 cellular models. Consistent with our hypothesis, MDA.Δ.C and MCF7.UGT8 cells, both with high GalCer content, showed lower activity of TNFRSF1B and TNFRSF9 promoters than MDA.Δ.UGT8.4 and MCF7.C cells lacking GalCer. In MCF7.UGT8 cells, activity of the BCL2 promoter was higher than in MCF7.C cells. In silico analysis performed to identify potential transcription factors that bind to promoter sequences and affect the activities of TNFRSF1B, TNFRSF9 and BCL2 promoters revealed that the binding sites for CREB1, NF-ĸB1, and p53 were present in all three promoters, indicating that they can be simultaneously involved in the regulation of TNFRSF1B, TNFRSF9 and BCL2 gene expression. No data about transcription factors involved in regulation of TNFRSF1B and TNFRSF9 gene are available, however, CREB1, NF-ĸB1 and p53 play significant regulatory roles in case of other TNFR family members, upregulating their expression [28,29]. In contrast, much more is known about the regulation of BCL2 by these three transcription factors [30,31,32]. CREB1 and NF-ĸB1 increase [30,41], and p53 downregulates [33] the expression of BCL2 gene. With regards to MCF7.UGT8 and T47D.UGT8 cells, accumulation of GalCer leads to downregulation of TNFRSF1B and TNFRSF9 genes and upregulation of the BCL2 gene in comparison to MCF7.C and T47D.C cells that do not produce GalCer, therefore the involvement of p53 in their simultaneous but opposite regulation seems to be probable. This assumption is further supported by data obtained from MDA-MB-231 cells. In this cellular model, no difference in promoter activity of the BCL2 gene was observed between MDA.Δ.C and MDA.Δ.UGT8.4 cells, despite significant differences in mRNA and protein levels. To explain these apparent contradictions, we assumed, based on the available information that p53 can be engaged in post-transcriptional regulation of Bcl-2 expression by increasing the Bcl-2 transcript stability [42]. Therefore, an actinomycin D assay was performed on MDA.Δ.UGT8.4 and control MDA.C cells, which revealed that the stability of Bcl-2 mRNA was considerably lower in BC cells that did not express GalCer. These data also revealed for the first time that one species of GSL can regulate the expression of specific genes through different mechanisms. The involvement of p53 in the regulation of TNFRSF1B, TNFRSF9, and BCL2 genes was supported by the Western blot analysis, which revealed that the accumulation of GalCer in BC cells is associated with decreased level of p53. These correlational observations are strongly supported by the siRNA-based RNA interference study, in which p53 expression was silenced. Direct evidence that p53 is involved in the regulation of TNFRSF1B, TNFRSF9 and BCL2 expression came from the experiments with siRNA-mediated suppression of p53 in MCF7.C and MDA.Δ.UGT8.4 cells. First, it was found that inhibition of p53 strongly reduced the expression of the TNFRSF1B and TNFRSF9 genes. Importantly, the activities of the TNFRSF1B and TNFRSF9 promoters were also greatly reduced in these cells, confirming our proposal that the expression of these genes is regulated at the transcriptional level. Second, inhibition of p53 increased the expression of the BCL2 gene. In MCF7.C cells, this was due to increased promoter activity. However, p53 inhibition had no effect on BCL2 gene expression at the mRNA level and its promoter activity in MDA.Δ.UGT8.4 cells. Therefore, in light of our results showing that p53 can regulate Bcl-2 expression by increasing Bcl-2 transcript stability, siRNA-treated MDA.Δ.UGT8.4 cells were subjected to an actinomycin D assay, which showed that p53 is indeed involved in the post-transcriptional regulation of Bcl-2.
Since p53 functions primarily as a transcription factor [43], it is not surprising that in MCF7 cells carrying wild-type p53 [44], regulation of all three target genes takes place at the level of promoter activity. However, in the case of MDA-MB-231 cells carrying a mutant form of p53 [45], only TNFRSF1B and TNFRSF9 genes were regulated at the level of transcription, which was not observed in the case of BCl2 gene. These differences can be explained by the studies showing that different p53 variants have an altered promoter activity spectrum [46]. In our case, the R280K mutation did not affect p53 binding to the TNFRSF1B and TNFRSF9 promoters, but prevented it from binding to the BCL2 promoter. Instead, such mutant p53 acquires new properties and binds to the 5'UTR of Bcl2 mRNA as it was shown in case of genetic variant of p53 with the substitution of proline by arginine within the proline-rich domain expressed by mucous cells [42].
Studies are in progress to identify the specific mechanisms by which GalCer decreases the level of p53 and in this way subsequently affects TNFRSF1B, TNFRSF9, and BCL2 expression. At this stage of our research, we propose the following hypothesis: GalCer localised in lipid rafts interacts with a specific protein, most likely a protein kinase, to trigger signalling pathway leading to reduce expression of p53, which in turn causes decreased expression of TNFRSF1B and TNFRSF9 and increased expression of BCL2 genes (Figure 7).

5. Conclusions

Our data showed that DR is an intrinsic property of a subpopulation of BC cells overexpressing UGT8 and therefore accumulating GalCer [47], and not a result of chemotherapy as it was observed in the case of GCS and accumulation of GlcCer in BC cells selected by stepwise exposure to drugs [48]. Furthermore, we showed that GalCer is involved in the regulation of specific apoptotic genes (TNFRSF1B, TNFRSF9 and BCL2) and is the only known GSL that can simultaneously affect extrinsic and intrinsic (mitochondrial) apoptotic pathways. Downregulation of death receptor signalling, which contributes to cancer cell sensitivity to cellular stressors, is in agreement with our proposal that GalCer facilitates tumour cell survival in the hostile tumour microenvironment [15]. However, upregulation of the mitochondrial pathway is mainly responsible for the induction of resistance to cytotoxic drugs, as shown in the present study. Therefore, we propose that the inhibition of UGT8 by knockout of the UGT8 gene or using chemical inhibitors of the UGT8 enzyme would increase cancer cell sensitivity to conventional chemotherapy.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: Detection of CRISPR/Cas9-induced indel mutations in the human UGT8 gene.; Figure S2: Proliferative potential of breast cancer (BC) MDA-MB-231, MCF7, and T47D cells with different expression of GalCer; Figure S3: Venn diagram of apoptotic genes shared between various cellular models representing BC cells with different content of GalCer; Figure S4: cDNA transcriptional microarray analysis of apoptotic genes in human BC cell lines with knockout of UGT8 gene; Figure S5: Sensitivity of BC MDA-MB-231 and MCF7 cells with different GalCer levels to apoptosis induced by DOX; Figure S6: Sensitivity of BC MDA-MB-231 and MCF7 cells with different GalCer levels to apoptosis incubated with TNFα. Figure S7: Nucleotide sequences of BCL2, TNFRSF9, and TNFRSF1B promoters with potential sequences for transcriptional factors; Figure S8: Expression of CREB and NFĸB in BC MDA-MB-231 and MCF7 cells with different GalCer levels. Table S1: List of primers used in this study, related to the experimental procedures; Table S2: List of primers used in apoptosis gene expression profiling with a RT-qPCR array.

Author Contributions

Investigation, J.S., S.R., A.U. W.W. and B.K.; visualization, J.S. and S.R.; methodology, J.S., S.R., A.U., W.W. and B.K.; writing - review & editing, J.S. and S.R.; conceptualization, M.U. and J.S.; supervision, M.U.; funding acquisition, M.U.; project administration M.U. and A.U.; writing - original draft, M.U.

Funding

Grant support for this work was provided by the National Science Center, Poland (Grant no. 2019/35/B/NZ5/01392).

Institutional Review Board Statement

The animal study protocol was approved by the Second Local Ethics Committee for Animal Experimentation (033/2019/P1,Wroclaw, Poland).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Kruh, G.D. Introduction to resistance to anticancer agents. Oncogene 2003, 22, 7262–7264. [Google Scholar] [CrossRef]
  2. Giussani, P.; Tringali, C.; Riboni, L.; Viani, P.; Venerando, B. Sphingolipids: Key Regulators of Apoptosis and Pivotal Players in Cancer Drug Resistance. Int J Mol Sci 2014, 15, 4356–4392. [Google Scholar] [CrossRef] [PubMed]
  3. Lavie, Y.; Cao, H.T.; Bursten, S.L.; Giuliano, A.E.; Cabot, M.C. Accumulation of glucosylceramides in multidrug-resistant cancer cells. J Biol Chem 1996, 271, 19530–19536. [Google Scholar] [CrossRef] [PubMed]
  4. Ryland, L.K.; Fox, T.E.; Liu, X.; Loughran, T.P.; Kester, M. Dysregulation of sphingolipid metabolism in cancer. Cancer Biol Ther 2011, 11, 138–149. [Google Scholar] [CrossRef] [PubMed]
  5. Gouaze, V.; Yu, J.Y.; Bleicher, R.J.; Han, T.Y.; Liu, Y.Y.; Wang, H.T.; Gottesman, M.M.; Bitterman, A.; Giuliano, A.E.; Cabot, M.C. Overexpression of glucosylceramide synthase and P-glycoprotein in cancer cells selected for resistance to natural product chemotherapy. Mol Cancer Ther 2004, 3, 633–639. [Google Scholar] [CrossRef]
  6. Liu, Y.Y.; Gupta, V.; Patwardhan, G.A.; Bhinge, K.; Zhao, Y.F.; Bao, J.X.; Mehendale, H.; Cabot, M.C.; Li, Y.T.; Jazwinski, S.M. Glucosylceramide synthase upregulates MDR1 expression in the regulation of cancer drug resistance through cSrc and beta-catenin signaling. Mol Cancer 2010, 9, 145. [Google Scholar] [CrossRef]
  7. De Rosa, M.F.; Sillence, D.; Ackerley, C.; Lingwood, C. Role of multiple drug resistance protein 1 in neutral but not acidic glycosphingolipid biosynthesis. J Biol Chem 2004, 279, 7867–7876. [Google Scholar] [CrossRef] [PubMed]
  8. De Rosa, M.F.; Ackerley, C.; Wang, B.; Ito, S.; Clarke, D.M.; Lingwood, C. Inhibition of multidrug resistance by AdamantylGb(3), a globotriaosylceramide analog. J Biol Chem 2008, 283, 4501–4511. [Google Scholar] [CrossRef]
  9. Noguchi, M.; Kabayama, K.; Uemura, S.; Kang, B.W.; Saito, M.; Igarashi, Y.; Inokuchi, J. Endogenously produced ganglioside GM3 endows etoposide and doxorubicin resistance by up-regulating Bcl-2 expression in 3LL Lewis lung carcinoma cells. Glycobiology 2006, 16, 641–650. [Google Scholar] [CrossRef]
  10. Kiura, K.; Watarai, S.; Ueoka, H.; Tabata, M.; Gemba, K.I.; Aoe, K.; Yamane, H.; Yasuda, T.; Harada, M. An alteration of ganglioside composition in cisplatin-resistant lung cancer cell line. Anticancer Res 1998, 18, 2957–2960. [Google Scholar]
  11. Tringali, C.; Lupo, B.; Cirillo, F.; Papini, N.; Anastasia, L.; Lamorte, G.; Colombi, P.; Bresciani, R.; Monti, E.; Tettamanti, G.; et al. Silencing of membrane-associated sialidase Neu3 diminishes apoptosis resistance and triggers megakaryocytic differentiation of chronic myeloid leukemic cells K562 through the increase of ganglioside GM3. Cell Death Differ 2009, 16, 164–174. [Google Scholar] [CrossRef] [PubMed]
  12. Schulte, S.; Stoffel, W. Ceramide UDPgalactosyltransferase from Myelinating Rat-Brain - Purification, Cloning, and Expression. P Natl Acad Sci USA 1993, 90, 10265–10269. [Google Scholar] [CrossRef] [PubMed]
  13. Kapitonov, D.; Yu, R.K. Cloning, characterization, and expression of human ceramide galactosyltransferase cDNA. Biochem Bioph Res Co 1997, 232, 449–453. [Google Scholar] [CrossRef]
  14. Marcus, J.; Popko, B. Galactolipids are molecular determinants of myelin development and axo-glial organization. Bba-Gen Subjects 2002, 1573, 406–413. [Google Scholar] [CrossRef] [PubMed]
  15. Owczarek, T.B.; Suchanski, J.; Pula, B.; Kmiecik, A.M.; Chadalski, M.; Jethon, A.; Dziegiel, P.; Ugorski, M. Galactosylceramide Affects Tumorigenic and Metastatic Properties of Breast Cancer Cells as an Anti-Apoptotic Molecule. Plos One 2013, 8, e84191. [Google Scholar] [CrossRef] [PubMed]
  16. Suchanski, J.; Grzegrzolka, J.; Owczarek, T.; Pasikowski, P.; Piotrowska, A.; Kocbach, B.; Nowak, A.; Dziegiel, P.; Wojnar, A.; Ugorski, M. Sulfatide decreases the resistance to stress-induced apoptosis and increases P-selectin-mediated adhesion: a two-edged sword in breast cancer progression. Breast Cancer Res 2018, 20, 133. [Google Scholar] [CrossRef]
  17. Kok, J.W.; Veldman, R.J.; Klappe, K.; Koning, H.; Filipeanu, C.M.; Muller, M. Differential expression of sphingolipids in MRP1 overexpressing HT29 cells. Int J Cancer 2000, 87, 172–178. [Google Scholar] [CrossRef]
  18. Veldman, R.J.; Klappe, K.; Hinrichs, J.; Hummel, I.; van der Schaaf, G.; Sietsma, H.; Kok, J.W. Altered sphingolipid metabolism in multidrug-resistant ovarian cancer cells is due to uncoupling of glycolipid biosynthesis in the Golgi apparatus. Faseb J 2002, 16, 1111. [Google Scholar] [CrossRef]
  19. Grazide, S.; Terrisse, A.D.; Lerouge, S.; Laurent, G.; Jaffrezou, J.P. Cytoprotective effect of glucosylceramide synthase inhibition against daunorubicin-induced apoptosis in human leukemic cell lines. J Biol Chem 2004, 279, 18256–18261. [Google Scholar] [CrossRef]
  20. Niimura, Y.; Nagai, K.I. Metabolic responses of sulfatide and related glycolipids in Madin-Darby canine kidney (MDCK) cells under osmotic stresses. Comp Biochem Phys B 2008, 149, 161–167. [Google Scholar] [CrossRef]
  21. Niimura, Y.; Moue, T.; Takahashi, N.; Nagai, K.I. Modification of sphingoglycolipids and sulfolipids in kidney cell lines under heat stress: Activation of monohexosylceramide synthesis as a ceramide scavenger. Glycobiology 2010, 20, 710–717. [Google Scholar] [CrossRef]
  22. Bieberich, E. Integration of glycosphingolipid metabolism and cell-fate decisions in cancer and stem cells: Review and hypothesis. Glycoconjugate J 2004, 21, 315–327. [Google Scholar] [CrossRef]
  23. Ogretmen, B. Sphingolipid metabolism in cancer signalling and therapy. Nat Rev Cancer 2018, 18, 33–50. [Google Scholar] [CrossRef]
  24. Ugorski, M.; Pahlsson, P.; Dus, D.; Nilsson, B.; Radzikowski, C. Glycosphingolipids of Human Urothelial Cell-Lines with Different Grades of Transformation. Glycoconjugate J 1989, 6, 303–318. [Google Scholar] [CrossRef] [PubMed]
  25. Souers, A.J.; Leverson, J.D.; Boghaert, E.R.; Ackler, S.L.; Catron, N.D.; Chen, J.; Dayton, B.D.; Ding, H.; Enschede, S.H.; Fairbrother, W.J.; et al. ABT-199, a potent and selective BCL-2 inhibitor, achieves antitumor activity while sparing platelets. Nat Med 2013, 19, 202–208. [Google Scholar] [CrossRef] [PubMed]
  26. Lapucci, A.; Lulli, M.; Amedei, A.; Papucci, L.; Witort, E.; Di Gesualdo, F.; Bertolini, F.; Brewer, G.; Nicolin, A.; Bevilacqua, A.; et al. zeta-Crystallin is a bcl-2 mRNA binding protein involved in bcl-2 overexpression in T-cell acute lymphocytic leukemia. Faseb J 2010, 24, 1852–1865. [Google Scholar] [CrossRef]
  27. Ishimaru, D.; Zuraw, L.; Ramalingam, S.; Sengupta, T.K.; Bandyopadhyay, S.; Reuben, A.; Fernandes, D.J.; Spicer, E.K. Mechanism of Regulation of bcl-2 mRNA by Nucleolin and A plus U-rich Element-binding Factor 1 (AUF1). J Biol Chem 2010, 285, 27182–27191. [Google Scholar] [CrossRef] [PubMed]
  28. Puimege, L.; Libert, C.; Van Hauwermeiren, F. Regulation and dysregulation of tumor necrosis factor. receptor-1. Cytokine Growth F R 2014, 25, 285–300. [Google Scholar] [CrossRef] [PubMed]
  29. Liu, X.G.; Yue, P.; Khuri, F.R.; Sun, S.Y. p53 upregulates death receptor 4 expression through an intronic p53 binding site. Cancer Res 2004, 64, 5078–5083. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, P.J.; Wang, J.; Lang, H.Y.; Wang, W.X.; Liu, X.H.; Liu, H.Y.; Tan, C.C.; Lia, X.T.; Zhao, Y.M.; Wu, X.H. Knockdown of CREB1 promotes apoptosis and decreases estradiol synthesis in mouse granulosa cells. Biomed Pharmacother 2018, 105, 1141–1146. [Google Scholar] [CrossRef] [PubMed]
  31. Heckman, C.A.; Mehew, J.W.; Boxer, L.M. NF-kappa B activates Bcl-2 expression in t(14;18) lymphoma cells. Oncogene 2002, 21, 3898–3908. [Google Scholar] [CrossRef]
  32. Bourgarel-Rey, V.; Savry, A.; Hua, G.Q.; Carre, M.; Bressin, C.; Chacon, C.; Imbert, J.; Braguer, D.; Barra, Y. Transcriptional down-regulation of Bcl-2 by vinorelbine: Identification of a novel binding site of p53 on Bcl-2 promoter. Biochem Pharmacol 2009, 78, 1148–1156. [Google Scholar] [CrossRef]
  33. Wu, Y.L.; Mehew, J.W.; Heckman, C.A.; Arcinas, M.; Boxer, L.M. Negative regulation of bcl-2 expression by p53 in hematopoietic cells. Oncogene 2001, 20, 240–251. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, T.H.; Wang, H.S.; Soong, Y.K. Paclitaxel-induced cell death - Where the cell cycle and apoptosis come together. Cancer 2000, 88, 2619–2628. [Google Scholar] [CrossRef] [PubMed]
  35. Ha, K.T.; Lee, Y.C.; Kim, C.H. Overexpression of GD3 synthase induces apoptosis of ECV304 cells through downregulation of Bcl-2 (vol 568, pg 183, 2004). Febs Lett 2004, 576, 499–499. [Google Scholar] [CrossRef]
  36. Rippo, M.R.; Malisan, F.; Ravagnan, L.; Tomassini, B.; Condo, I.; Costantini, P.; Susin, S.A.; Rufini, A.; Todaro, M.; Kroemer, G.; et al. GD3 ganglioside directly targets mitochondria in a bcl-2-controlled fashion. Faseb J 2000, 14, 2047–2054. [Google Scholar] [CrossRef] [PubMed]
  37. Scorrano, L.; Petronilli, V.; Colonna, R.; Di Lisa, F.; Bernardi, P. Chloromethyltetramethylrosamine (Mitotracker Orange (TM)) induces the mitochondrial permeability transition and inhibits respiratory complex I - Implications for the mechanism of cytochrome c release. J Biol Chem 1999, 274, 24657–24663. [Google Scholar] [CrossRef] [PubMed]
  38. Collel, A.; Garcia-Ruiz, C.; Roman, J.; Ballesta, A.; Fernandez-Checa, J.C. Ganglioside GD3 enhances apoptosis by suppressing the nuclear factor-κB-dependent survival pathway. 2001. [Google Scholar] [CrossRef]
  39. Paris, R.; Morales, A.; Coll, O.; Sanchez-Reyes, A.; Garcia-Ruiz, C.; Fernandez-Checa, J.C. Ganglioside GD3 sensitizes human hepatoma cells to cancer therapy. J Biol Chem 2002, 277, 49870–49876. [Google Scholar] [CrossRef] [PubMed]
  40. Kakugawa, Y.; Wada, T.; Yamaguchi, K.; Yamanami, H.; Ouchi, K.; Sato, I.; Miyagi, T. Up-regulation of plasma membrane-associated ganglioside sialidase (Neu3) in human colon cancer and its involvement in apoptosis suppression. P Natl Acad Sci USA 2002, 99, 10718–10723. [Google Scholar] [CrossRef]
  41. Li, Z.Q.; Xuan, Z.F.; Chen, J.; Song, W.F.; Zhang, S.Y.; Jin, C.; Zhou, M.Q.; Zheng, S.S.; Song, P.H. Inhibiting the NF-kappa B pathway enhances the antitumor effect of cabazitaxel by downregulating Bcl-2 in pancreatic cancer. Int J Oncol 2020, 57, 161–170. [Google Scholar] [CrossRef]
  42. Chand, H.S.; Montano, G.; Huang, X.S.; Randell, S.H.; Mebratu, Y.; Petersen, H.; Tesfaigzi, Y. A genetic variant of p53 restricts the mucous secretory phenotype by regulating SPDEF and Bcl-2 expression. Nat Commun 2014, 5, 6567. [Google Scholar] [CrossRef]
  43. Clore, G.M.; Omichinski, J.G.; Sakaguchi, K.; Zambrano, N.; Sakamoto, H.; Appella, E.; Gronenborn, A.M. High-Resolution Structure of the Oligomerization Domain of P53 by Multidimensional Nmr. Science 1994, 265, 386–391. [Google Scholar] [CrossRef]
  44. Lu, X.H.; Errington, J.; Curtin, N.J.; Lunec, J.; Newell, D.R. The impact of status on cellular sensitivity to antifolate drugs. Clin Cancer Res 2001, 7, 2114–2123. [Google Scholar]
  45. Katayose, D.; Gudas, J.; Nguyen, H.; Srivastava, S.; Cowan, K.H.; Seth, P. Cytotoxic effects of adenovirus-mediated wild-type p53 protein expression in normal and tumor mammary epithelial cells. Clin Cancer Res 1995, 1, 889–897. [Google Scholar]
  46. Jordan, J.J.; Inga, A.; Conway, K.; Edmiston, S.; Carey, L.A.; Wu, L.; Resnick, M.A. Altered-Function p53 Missense Mutations Identified in Breast Cancers Can Have Subtle Effects on Transactivation. Mol Cancer Res 2010, 8, 701–716. [Google Scholar] [CrossRef] [PubMed]
  47. Dziegiel, P.; Owczarek, T.; Plazuk, E.; Gomulkiewicz, A.; Majchrzak, M.; Podhorska-Okolow, M.; Driouch, K.; Lidereau, R.; Ugorski, M. Ceramide galactosyltransferase (UGT8) is a molecular marker of breast cancer malignancy and lung metastases. Brit J Cancer 2010, 103, 524–531. [Google Scholar] [CrossRef] [PubMed]
  48. Gottesman, M.M.; Fojo, T.; Bates, S.E. Multidrug resistance in cancer: Role of ATP-dependent transporters. Nat Rev Cancer 2002, 2, 48–58. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Characteristics of BC cells with knock-out or overexpression of UGT8 gene. A(I) Level of UGT8 protein expression in parental MDA-MB-231, MDA-MB-231 cells transduced with “empty” vector (MDA.Δ.C), and MDA-MB-231 cell clones with knockout of UGT8 gene (MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4). A(II) Immunostaining of neutral GSLs from MDA-MB-231, MDA.Δ.C, MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4 cells. B(I) Level of UGT8 protein expression in parental MCF7 cells, MCF7 cells transduced with “empty” vector (MCF7.C), and MCF7 cells overexpressing UGT8 (MCF7.UGT8). B(II) Immunostaining of neutral GSLs from MCF7, MCF7.C, and MCF7.UGT8. The expression of UGT8 was analysed by Western blotting using rabbit polyclonal antibodies directed against UGT8. Cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control. For GalCer immunostaining, samples of neutral GSLs that had been purified from 1×107 cells were separated on HP-TLC plate. (C) Viability of BC MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4 cells with knockout of the UGT8 gene and MCF7.UGT8 and T47D.UGT8 cells overexpressing UGT8 and GalCer exposed to increased concentrations of anti-cancer drugs: doxorubicin (DOX) (I), paclitaxel (PTX) (II), and cisplatin cis-diamminedichloroplatinum (III) (CDDP) (III) for 48 h. Percentage of viable cells was determined by MTT assay. Data represent the mean ± standard deviation (SD) of 6 replicates from 2 independent measurements. Differences that are statistically significant (*p < 0.1, **p < 0.01).
Figure 1. Characteristics of BC cells with knock-out or overexpression of UGT8 gene. A(I) Level of UGT8 protein expression in parental MDA-MB-231, MDA-MB-231 cells transduced with “empty” vector (MDA.Δ.C), and MDA-MB-231 cell clones with knockout of UGT8 gene (MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4). A(II) Immunostaining of neutral GSLs from MDA-MB-231, MDA.Δ.C, MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4 cells. B(I) Level of UGT8 protein expression in parental MCF7 cells, MCF7 cells transduced with “empty” vector (MCF7.C), and MCF7 cells overexpressing UGT8 (MCF7.UGT8). B(II) Immunostaining of neutral GSLs from MCF7, MCF7.C, and MCF7.UGT8. The expression of UGT8 was analysed by Western blotting using rabbit polyclonal antibodies directed against UGT8. Cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control. For GalCer immunostaining, samples of neutral GSLs that had been purified from 1×107 cells were separated on HP-TLC plate. (C) Viability of BC MDA.Δ.UGT8.1 and MDA.Δ.UGT8.4 cells with knockout of the UGT8 gene and MCF7.UGT8 and T47D.UGT8 cells overexpressing UGT8 and GalCer exposed to increased concentrations of anti-cancer drugs: doxorubicin (DOX) (I), paclitaxel (PTX) (II), and cisplatin cis-diamminedichloroplatinum (III) (CDDP) (III) for 48 h. Percentage of viable cells was determined by MTT assay. Data represent the mean ± standard deviation (SD) of 6 replicates from 2 independent measurements. Differences that are statistically significant (*p < 0.1, **p < 0.01).
Preprints 91808 g001
Figure 2. Characteristics of murine mammary carcinoma 4T1 cells overexpressing murine UGT8a and GalCer. A(I) Expression of UGT8a protein in parental 4T1 cells, control 4T1 cells transduced with vector alone (4T1.C), and 4T1 cells overexpressing UGT8a (4T1.UGT8a); A(II) immunostaining of neutral GSLs from 4T1, 4T1.C and 4T1.UGTa cells. Western blotting with anti-human UGT8 antibodies was used to analyze the expression of UGT8a. Cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control For GalCer detection, samples of neutral GSLs purified from 1×107 cells was applied to an HP-TLC plate. (B) Viability of 4T1.C and 4T1.UGT8a cells overexpressing UGT8a and GalCer incubated with increased concentrations of DOX for 48 h. Percentage of viable cells was determined using MTT assay. Data represent the mean ±SD of 6 replicates from 2 independent measurements. Statistically significant differences (*p < 0.1). (C) Impact of i. v. administration of DOX on the growth of 4T1.C and 4T1.UGT8a tumors in BALB/c mice. 4T1.UGT8a – mice with 4T1.UGT8a tumors overexpressing UGT8a and GalCer treated with placebo (PBS); 4T1.UGT8 + DOX – mice with 4T1.UGT8a tumors treated with DOX; 4T1.C - mice with control 4T1.C tumors treated with placebo (PBS); 4T1.C + DOX - mice with 4T1.C tumours treated with DOX. Data was presented as mean tumor volume (±SD) for group of mice (n = 5) of the specified time intervals. Data was analysed using Prism 5.0 Two-way Anova Dunnetts post-test (p* < 0.05).
Figure 2. Characteristics of murine mammary carcinoma 4T1 cells overexpressing murine UGT8a and GalCer. A(I) Expression of UGT8a protein in parental 4T1 cells, control 4T1 cells transduced with vector alone (4T1.C), and 4T1 cells overexpressing UGT8a (4T1.UGT8a); A(II) immunostaining of neutral GSLs from 4T1, 4T1.C and 4T1.UGTa cells. Western blotting with anti-human UGT8 antibodies was used to analyze the expression of UGT8a. Cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control For GalCer detection, samples of neutral GSLs purified from 1×107 cells was applied to an HP-TLC plate. (B) Viability of 4T1.C and 4T1.UGT8a cells overexpressing UGT8a and GalCer incubated with increased concentrations of DOX for 48 h. Percentage of viable cells was determined using MTT assay. Data represent the mean ±SD of 6 replicates from 2 independent measurements. Statistically significant differences (*p < 0.1). (C) Impact of i. v. administration of DOX on the growth of 4T1.C and 4T1.UGT8a tumors in BALB/c mice. 4T1.UGT8a – mice with 4T1.UGT8a tumors overexpressing UGT8a and GalCer treated with placebo (PBS); 4T1.UGT8 + DOX – mice with 4T1.UGT8a tumors treated with DOX; 4T1.C - mice with control 4T1.C tumors treated with placebo (PBS); 4T1.C + DOX - mice with 4T1.C tumours treated with DOX. Data was presented as mean tumor volume (±SD) for group of mice (n = 5) of the specified time intervals. Data was analysed using Prism 5.0 Two-way Anova Dunnetts post-test (p* < 0.05).
Preprints 91808 g002
Figure 3. Levels of Bcl-2, TNFRS9, and TNFRS1B expression in BC MDA-MB-231, MCF7, and T47D cells with various contents of GalCer. (A) qPCR was used to determine mRNA levels, which were normalized against β-actin and cells with low expression levels of UGT8 mRNA served as calibrator sample. Data is expressed as mean ±SD. (B) Western blotting combined with densitometric analysis were used to quantify protein expression levels in BC cell lines. Specific cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control. Band intensities were analyzed with Image Lab software (Bio-Rad Laboratories) and GraphPad Prism version 5. (*p < 0.05). All values are means ±SD.
Figure 3. Levels of Bcl-2, TNFRS9, and TNFRS1B expression in BC MDA-MB-231, MCF7, and T47D cells with various contents of GalCer. (A) qPCR was used to determine mRNA levels, which were normalized against β-actin and cells with low expression levels of UGT8 mRNA served as calibrator sample. Data is expressed as mean ±SD. (B) Western blotting combined with densitometric analysis were used to quantify protein expression levels in BC cell lines. Specific cellular proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto nitrocellulose membrane. GAPDH was used as the internal control. Band intensities were analyzed with Image Lab software (Bio-Rad Laboratories) and GraphPad Prism version 5. (*p < 0.05). All values are means ±SD.
Preprints 91808 g003
Figure 4. Sensitivity to DOX-induced apoptosis of BC (A) MDA-MB-231 and (B) MCF7 cells with various GalCer levels and incubated with ABT-199 (specific inhibitor of Bcl-2). Cells were cultured with ABT-199 at a concentration of 0.125 µM or ABT-199 and DOX at concentrations of 0.125 µM and 1 µM, respectively, for 48 h. Sensitivity to TNFα-induced apoptosis of BC (C) MDA-MB-231 and (D) MCF7 cells with different GalCer content. Cells were cultured in the presence of TNFα at a concentration of 25 µM for 48 h. Apoptotic cells were determined by flow cytometry using Annexin V and SYTOX Green stain. MDA.Δ.C - MDA-MB-231 cells transduced with vector alone; MDA.Δ.UGT8.4 - MDA-MB-231 cell clones with knockout of the UGT8 gene; MCF7.C - MCF7 cells transduced with vector alone; MCF7.UGT8 – MCF7 cells overexpressing UGT8 and GalCer. Data represent the mean ±SD of 6 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Figure 4. Sensitivity to DOX-induced apoptosis of BC (A) MDA-MB-231 and (B) MCF7 cells with various GalCer levels and incubated with ABT-199 (specific inhibitor of Bcl-2). Cells were cultured with ABT-199 at a concentration of 0.125 µM or ABT-199 and DOX at concentrations of 0.125 µM and 1 µM, respectively, for 48 h. Sensitivity to TNFα-induced apoptosis of BC (C) MDA-MB-231 and (D) MCF7 cells with different GalCer content. Cells were cultured in the presence of TNFα at a concentration of 25 µM for 48 h. Apoptotic cells were determined by flow cytometry using Annexin V and SYTOX Green stain. MDA.Δ.C - MDA-MB-231 cells transduced with vector alone; MDA.Δ.UGT8.4 - MDA-MB-231 cell clones with knockout of the UGT8 gene; MCF7.C - MCF7 cells transduced with vector alone; MCF7.UGT8 – MCF7 cells overexpressing UGT8 and GalCer. Data represent the mean ±SD of 6 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Preprints 91808 g004
Figure 5. Promoter activities of (A) TNFRS1B, (B) TNFRS9, and (C) BCL2 genes cloned into pGL3 basic vector after transfection into BC cells with various contents of GalCer. MDA.Δ.C - MDA-MB-231 cells transduced with “empty” vector; MDA.Δ.UGT8.4 - MDA-MB-231 cell clones with knock-out of the UGT8 gene; MCF7.C – MCF7 cells transduced with “empty” vector; MCF7.UGT8 – MCF7 cells overexpressing UGT8 and GalCer. The promoter activities were measured using the dual-luciferase reporter assay system. The bars represent average luciferase activities compared with the control pGL3 vector. Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1). (D) GalCer increases the stability of Bcl-2 mRNA in BC cells. MDA.Δ.C and MDA.Δ.UGT8.4 cells were treated with actinomycin D for 2 - 6 h and after indicated time points, the Bcl-2 mRNA expression levels, which were normalised against β-actin were evaluated by qPCR. MDA MB 231.C and MDA MB 231.Δ4 cells not treated with actinomycin D were used as calibrator sample. The results are shown as the percentage of the control (0 h). Data represent the mean ±SD of 6 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Figure 5. Promoter activities of (A) TNFRS1B, (B) TNFRS9, and (C) BCL2 genes cloned into pGL3 basic vector after transfection into BC cells with various contents of GalCer. MDA.Δ.C - MDA-MB-231 cells transduced with “empty” vector; MDA.Δ.UGT8.4 - MDA-MB-231 cell clones with knock-out of the UGT8 gene; MCF7.C – MCF7 cells transduced with “empty” vector; MCF7.UGT8 – MCF7 cells overexpressing UGT8 and GalCer. The promoter activities were measured using the dual-luciferase reporter assay system. The bars represent average luciferase activities compared with the control pGL3 vector. Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1). (D) GalCer increases the stability of Bcl-2 mRNA in BC cells. MDA.Δ.C and MDA.Δ.UGT8.4 cells were treated with actinomycin D for 2 - 6 h and after indicated time points, the Bcl-2 mRNA expression levels, which were normalised against β-actin were evaluated by qPCR. MDA MB 231.C and MDA MB 231.Δ4 cells not treated with actinomycin D were used as calibrator sample. The results are shown as the percentage of the control (0 h). Data represent the mean ±SD of 6 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Preprints 91808 g005
Figure 6. (A) P53 expression in MDA-MB-231 cells transduced with “empty” vector (MDA.Δ.C), MDA-MB-231 cells with knockout of the UGT8 gene and lack of GalCer (MDA.Δ.UGT8.4), MCF7 cells transduced with vector alone (MCF7.C), and MCF7 cells overexpressing UGT8 and GalCer (MCF7.UGT8). Western blotting was used to analyse p53 expression in BC cell lines. Rabbit anti-p53 polyclonal antibodies were used to detect p53. Proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto a nitrocellulose membrane. GAPDH was used as an internal control. (B) TNFRS1B, TNFRS9, and BCL2-2 gene expression levels in MDA.Δ.UGT8.4 and MCF.7 cells transfected with siRNA directed against p53 mRNA. qPCR was used to determine mRNA levels, normalized against β-actin. Cells with low expression levels of UGT8 mRNA were used as calibrators. Data are expressed as mean ±SD. (C) TNFRS1B, TNFRS9 and BCL2 promoter activities in MDA.Δ.UGT8.4 and MCF.7 cells transfected with siRNA targeting p53 mRNA. The activities of the promoters cloned into the pGL3 basic vector were measured using the dual-luciferase reporter assay system. The bars represent the average luciferase activities compared to the control pGL3 vector. (D) The stability of Bcl-2 mRNA in MDA.Δ.UGT8.4 cells transfected with siRNA targeting p53 mRNA. Cells were treated with actinomycin D for 2 - 6 h and after indicated time points, the Bcl-2 mRNA expression levels, which were normalized against β-actin were evaluated by qPCR. MDA.Δ.UGT8.4 cells not treated with actinomycin D were used as calibrator samples. The results are shown as the percentage of the control (0 h). MDA.Δ.UGT8.4/siTP53 and MCF7.C/siTP53 – BC cells transfected with siRNA directed against p53 mRNA; MDA.Δ.UGT8.4 and MCF7.C – control BC cells transfected with scrambled siRNA. Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1). Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Figure 6. (A) P53 expression in MDA-MB-231 cells transduced with “empty” vector (MDA.Δ.C), MDA-MB-231 cells with knockout of the UGT8 gene and lack of GalCer (MDA.Δ.UGT8.4), MCF7 cells transduced with vector alone (MCF7.C), and MCF7 cells overexpressing UGT8 and GalCer (MCF7.UGT8). Western blotting was used to analyse p53 expression in BC cell lines. Rabbit anti-p53 polyclonal antibodies were used to detect p53. Proteins (40 µg) separated by SDS-PAGE on 12% gel were transferred onto a nitrocellulose membrane. GAPDH was used as an internal control. (B) TNFRS1B, TNFRS9, and BCL2-2 gene expression levels in MDA.Δ.UGT8.4 and MCF.7 cells transfected with siRNA directed against p53 mRNA. qPCR was used to determine mRNA levels, normalized against β-actin. Cells with low expression levels of UGT8 mRNA were used as calibrators. Data are expressed as mean ±SD. (C) TNFRS1B, TNFRS9 and BCL2 promoter activities in MDA.Δ.UGT8.4 and MCF.7 cells transfected with siRNA targeting p53 mRNA. The activities of the promoters cloned into the pGL3 basic vector were measured using the dual-luciferase reporter assay system. The bars represent the average luciferase activities compared to the control pGL3 vector. (D) The stability of Bcl-2 mRNA in MDA.Δ.UGT8.4 cells transfected with siRNA targeting p53 mRNA. Cells were treated with actinomycin D for 2 - 6 h and after indicated time points, the Bcl-2 mRNA expression levels, which were normalized against β-actin were evaluated by qPCR. MDA.Δ.UGT8.4 cells not treated with actinomycin D were used as calibrator samples. The results are shown as the percentage of the control (0 h). MDA.Δ.UGT8.4/siTP53 and MCF7.C/siTP53 – BC cells transfected with siRNA directed against p53 mRNA; MDA.Δ.UGT8.4 and MCF7.C – control BC cells transfected with scrambled siRNA. Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1). Data represent the mean ±SD of 3 replicates from 2 independent experiments. Statistically significant differences (*p < 0.1).
Preprints 91808 g006
Figure 7. Diagram showing the putative regulatory role of GalCer in regulation of TNFRSF1B, TNFRSF9 and BCL2 gene expression that results in resistance of BC cells to apoptosis (Created with BioRender.com). As hypothesized, GalCer modulates the activities of signaling pathway/pathways which leads to decreased expression of p53. Depending on the cellular context, p53, acting as a transcription factor, upregulates expression of the BCL2 gene and downregulates the expression of TNFRSF1B and TNFRSF9 genes at the transcriptional level, or p53 acts as mRNA-binding protein (RBP) and upregulates expression of Bcl-2 protein at the post-translational level. ↑ - upregulation, ↓ - downregulation.
Figure 7. Diagram showing the putative regulatory role of GalCer in regulation of TNFRSF1B, TNFRSF9 and BCL2 gene expression that results in resistance of BC cells to apoptosis (Created with BioRender.com). As hypothesized, GalCer modulates the activities of signaling pathway/pathways which leads to decreased expression of p53. Depending on the cellular context, p53, acting as a transcription factor, upregulates expression of the BCL2 gene and downregulates the expression of TNFRSF1B and TNFRSF9 genes at the transcriptional level, or p53 acts as mRNA-binding protein (RBP) and upregulates expression of Bcl-2 protein at the post-translational level. ↑ - upregulation, ↓ - downregulation.
Preprints 91808 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated