Acessibilidade / Reportar erro

Transcriptomics analysis identified ezrin as a potential druggable target in cervical and gastric cancer cells

Abstract

Objective:

Cancer genomics and transcriptomics studies have provided a large volume of data that enables to test of hypotheses based on real data from cancer patients. Ezrin (encoded by the EZR gene) is a highly expressed protein in cancer that contributes to linking the actin cytoskeleton to the cell membrane and signal transduction pathways involved in oncogenesis and disease progression. NSC305787 is a pharmacological ezrin inhibitor with potential antineoplastic effects. In the present study, the authors prospected EZR mRNA levels in a pan-cancer analysis and identified potential cancers that could benefit from anti-EZR therapies.

Methods:

This study analyzed TCGA data for 32 cancer types, emphasizing cervical squamous cell carcinoma and stomach adenocarcinoma. It investigated the impact of EZR transcript levels on clinical outcomes and identified differentially expressed genes. Cell lines were treated with NSC305787, and its effects were assessed through various cellular and molecular assays.

Results:

EZR mRNA levels are highly expressed, and their expression is associated with biologically relevant molecular processes in cervical squamous carcinoma and stomach adenocarcinoma. In cellular models of cervical and gastric cancer, NSC305787 reduces cell viability and clonal growth (p < 0.05). Molecular analyses indicate that the pharmacological inhibition of EZR induces molecular markers of cell death and DNA damage, in addition, to promoting the expression of genes associated with apoptosis and inhibiting the expression of genes related to survival and proliferation.

Conclusion:

The present findings provide promising evidence that ezrin may be a molecular target in the treatment of cervical and gastric carcinoma.

Keywords:
Ezrin; Transcriptomics; Cervical cancer; Gastric cancer; Pharmacology

HIGHLIGHTS

Ezrin (EZR) is involved in oncogenesis and disease progression.

EZR is highly expressed in cervical squamous cell carcinoma and stomach adenocarcinoma.

NSC305787, an EZR inhibitor, reduces cell viability and clonal growth in cervical and gastric cancer cells.

EZR may be a molecular target in the treatment of cervical and gastric carcinoma.

Introduction

In recent years, cancer genomics and transcriptomics studies have provided and publicly deposited a large volume of data and driven the development of analysis platforms for generating and testing hypotheses based on real data from cancer patients.11 Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov 2012;2(5):401–4.,22 Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res 2020;48(W1):W509–14. Among these studies, the initiative of The Cancer Genome Atlas (TCGA) stands out. These data allow evaluation of the impact of mutations or gene expression on clinical, biological, molecular, or survival outcomes.

Ezrin (encoded by EZR) is a member of the ERM (ezrin, radixin, moesin) protein family, and when phosphorylated, it is essential for linking the actin cytoskeleton to the cell membrane.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79. Ezrin binds to integral membrane proteins, adhesion molecules, multidrug resistance proteins, scaffold proteins, Rho-related proteins, tyrosine kinases, and PIP2.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79. In normal tissues, EZR is widely expressed, but its high levels are time-specific during human development.44 Xie JJ, Zhang FR, Tao LH, Lu Z, Xu XE, Jian S, et al. Expression of ezrin in human embryonic, fetal, and normal adult tissues. J Histochem Cytochem 2011;59(11):1001–8.,55 Xu WT, Shi LL, Xu J, Qian H, Zhou H, Wang LH. Ezrin expression in female reproductive tissues: A review of regulation and pathophysiological implications. Front Cell Dev Biol 2023;11:1125881. This protein also contributes to signal transduction pathways involved in cancer development and progression, including PI3K/AKT/mTOR and RhoGTPase signaling.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79. In cancer cells, EZR-mediated AKT activation modulates the BCL2 protein family, favoring survival and proliferation.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79. Thus, there is growing evidence that EZR may serve as a proto-oncogene, and several studies have demonstrated its high expression in tumors compared to their normal histological counterpart.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79.

Given the above, EZR has drawn up attention as a potential pharmacological target to compose the antineoplastic arsenal. There are two synthetic compounds described as selective ezrin inhibitors (NSC668394 and NSC305787) that have been better characterized in preclinical studies, NSC305787 being considered more potent and exhibiting a better pharmacokinetic profile.66 Bulut G, Hong SH, Chen K, Beauchamp EM, Rahim S, Kosturko GW, et al. Small molecule inhibitors of ezrin inhibit the invasive phenotype of osteosarcoma cells. Oncogene 2012;31(3):269–81.,77 Celik H, Bulut G, Han J, Graham GT, Minas TZ, Conn EJ, et al. Ezrin Inhibition Up-regulates Stress Response Gene Expression. J Biol Chem 2016;291(25):13257–70.

In the present study, the authors examined EZR mRNA levels across various cancers through a comprehensive pan-cancer analysis. These findings highlighted cancer types exhibiting notably elevated EZR expression, specifically cervical and gastric cancer. Consequently, the authors posited that these particular tumor types might offer promising targets for anti-EZR therapies. To explore this hypothesis further, the authors evaluated the efficacy of an EZR inhibitor, NSC305787, in representative cellular models of these diseases.

Material and methods

Bioinformatics and functional genomics analysis

Data from TCGA cohorts that included 32 types of cancer, of which 21 had data from adjacent normal tissue, were obtained from the TIMER 2.0 database (http://timer.cistrome.org/). Statistical significance was calculated using the Wilcoxon test.

Gene expression and clinical data from the TCGA cervical squamous cell carcinoma (n = 294) and stomach adenocarcinoma (n = 412) studies were obtained from cBioPortal (https://www.cbioportal.org/).88 TCGA. Integrated genomic and molecular characterization of cervical cancer. Nature 2017;543(7645):378–84.,99 TCGA. Comprehensive molecular characterization of gastric adenocarcinoma. Nature 2014;513(7517):202–9. The impact of EZR transcript levels on clinical outcomes (Overall Survival [OS], Disease-Free Survival [DFS], and Disease-Specific Survival [DSS]) was investigated. Dichotomization was realized according to the ROC curve and its respective Area Under the Curve (AUC) and the C-index using a maximization metric provided by the R package Cutpoint.1010 Heagerty PJ, Zheng Y. Survival model predictive accuracy and ROC curves. Biometrics 2005;61(1):92–105. The observational study follows the STROBE statement.

All transcripts from the RNAseq of the TCGA cervical squamous cell carcinoma and stomach adenocarcinoma cohorts were pre-ranked according to their differential expression by comparing samples with high and low EZR expression. Normalized quantile gene expression was used for classification using the limma-voom package at Galaxy (https://usegalaxy.org/). A heatmap was constructed using ClusterVis (https://biit.cs.ut.ee/clustvis/) to represent the top 50 differentially expressed genes between high and low EZR expression samples. Volcano plots correlating the adjusted -log10 p-value and log2-fold-change and Spearman correlation plot were constructed in GraphPad Prism 8.0 (GraphPad Software, Inc., San Diego, CA, USA). All differentially expressed genes obtained from the Galaxy tool (fold-change > 1.5 and p < 0.05) were used for Gene Ontology (GO) analysis using ShinyGo 0.77 database (http://bioinformatics.sdstate.edu/go/), and the top twenty upregulated and downregulated GO biological processes, GO molecular processes, and GO cellular components were illustrated.

Cell culture and reagent chemicals

A-431 (cervix squamous cell carcinoma) and HGC-27 (metastatic gastric cancer) cells were obtained from Banco de Células do Rio de Janeiro (BCRJ; Rio de Janeiro, Brazil). AGS (gastric adenocarcinoma) cells were kindly provided by Prof. Tiago Góss dos Santos (A. C. Camargo Cancer Center, São Paulo, Brazil). Cell culture conditions were performed according to the recommendations of the BCRJ and American Type Culture Collection (ATCC). NSC305787 was purchased from ChemScene (Monmouth Junction, NJ, USA) and prepared as a 10 mM stock solution in dimethyl sulfoxide (Me2SO4; DMSO).

Cell viability assay

Cell viability was determined by a Methylthiazoletetrazolium (MTT) assay. A total of 5 × 103 cells per well were plated in a 96-well plate and exposed to vehicle or increasing concentrations of NSC305787 (ranging from 0.4 to 50 µM) for 24h. Then, 10 μL MTT solution (5 mg/mL) was added, and the cells were incubated at 37°C in 5% CO2 for 4h. The reaction was stopped by adding 100 μL 0.1N HCl in anhydrous isopropanol. Cell viability was evaluated by measuring the absorbance at 570 nm. The IC50 values were calculated by performing a nonlinear regression analysis in GraphPad Prism 8 (GraphPad Software, Inc.).

Colony formation assay

Cancer cell lines (1 × 103 cells/35 mm2 plate) were incubated with vehicle or different concentrations of NSC305787. Colonies were detected after 10-15 days of culture by staining with 0.5% crystal violet (Sigma-Aldrich, St. Louis, MO, USA) in 10% ethanol. Images were acquired using the G:BOX Chemi XRQ (Syngene, Cambridge, UK) and analyzed using ImageJ software (US National Institutes of Health, Bethesda, MD, USA).

Morphology analysis by immunofluorescence

A total of cells at 1 × 104 per well were seeded in a 24-well plate with coverslips at the bottom. The plate was kept at 37°C in an incubator with 5% CO2 for 48 hours. Next, treatment with NSC305787 was carried out at concentrations of 1.6 and 2.4 μM for A-431 cells, 3.2 and 4.0 μM for AGS cells, and 4.0 and 4.8 μM for HGC-27 cells. After 24h of treatment, the culture medium was removed, the wells were washed 2 times with PBS and the cells were fixed with 3.7% formaldehyde for 15 minutes, permeabilized with 0.5% triton for 10 minutes, and incubated with phalloidin (1:400, Thermo Fisher Scientific, San Jose, CA, USA) for 1h. Then, a new wash with PBS was performed and the slides were prepared with a mounting medium for fluorescence (ProLongTM Gold antifade reagent with DAPI; Thermo Fisher Scientific). Images were performed using fluorescence microscopy (LionHeart FX automated microscope, Biotek, Winooski, VT, USA).

Cell death analysis by flow cytometry

A total of 1 × 105 cells per well were seeded in 12-well plates in the presence of either the vehicle or NSC305787 (1.6, 2.4, and 4.8 μM for A-431 cells, 3.2, 4, and 8 μM for AGS cells, and 4, 4.8 and 9.6 μM for HGC-27 cells) for 24h. Subsequently, the cells were washed twice with ice-cold PBS and resuspended in a binding buffer containing 1 μg/mL Propidium Iodide (PI) and 1 μg/mL FITC-labeled annexin V (BD Biosciences, San Jose, CA, USA). All samples were acquired by flow cytometry (FACS-Calibur; Becton Dickinson, Franklin Lakes, NJ, USA) after 15-minute incubation at room temperature in a light-protected area and analyzed using FlowJo software (Treestar, Inc., San Carlos, CA, USA).

Western blot analysis

A-431, AGS, and HGC-27 cells were exposed to vehicle or NSC305787 for 24h. Total protein was extracted using a buffer containing 100 mM Tris (pH 7.6), 1% Triton X-100, 2 mM PMSF, 10 mM Na3VO4, 100 mM NaF, 10 mM Na4P2O7, and 10 mM EDTA. Equal amounts of total protein were separated by SDS-PAGE and identified by western blotting with the relevant antibodies listed below. Protein signals were detected using the SuperSignalTM West Dura Extended Duration Substrate System (Thermo Fisher Scientific) and the G:BOX Chemi XX6 gel doc system (Syngene). Antibodies against PARP1 (#9542), γH2AX (#9718), and α-tubulin (#2144) were purchased from Cell Signaling Technology (Danvers, MA, USA).

Quantitative PCR (qPCR)

A-431, AGS, and HGC-27 were exposed to vehicle or NSC305787 for 24h. Total RNA was obtained using TRIzol reagent (Thermo Fisher Scientific). The cDNA was synthesized from 1 µg of RNA using a High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific). Quantitative PCR (qPCR) was performed using a QuantStudio 3 RealTime PCR System (Thermo Fisher Scientific) and a SybrGreen System and specific primers (Supplementary Table 1), using HPRT1 and ACTB as reference genes. Relative quantification values were calculated using the 2-ΔΔCT equation.1111 Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001;25(4):402–8. A negative ‘No Template Control’ was included for each primer pair.

Statistical analysis

Statistical analyses were performed using GraphPad Prism 8 (Graph-Pad Software, Inc.), SAS System for Windows 9.2 (SAS Institute, Inc., NC, USA), or Stata software (Stata Corp., College Station, TX, USA). The Cox proportional hazard model for regression analysis or a log-rank test (Mantel-Cox) was used to estimate OS, DFS, and DSS. OS was defined as the time from sample collection to the date of death; living patients were censored on the date of the last assessment. DFS was defined as the time interval between treatment and the date of disease progression, relapse, treatment failure, or death; patients alive and without disease progression were censored on the date of the last assessment. DSS was defined as the time from diagnosis or initiation of treatment to the time of death, specifically from the disease. For comparisons, the Student t-test or ANOVA and the Bonferroni post-test were used. A p-value <0.05 was considered statistically significant.

Results

The pan-cancer analysis identified cancers with high EZR expression

To begin the exploratory analysis, the authors took advantage of public databases to investigate the transcript levels of EZR in different cohorts of TCGA cancer patients. A significant increase in EZR expression was observed in patients with breast invasive carcinoma, cervical and endocervical cancer, cholangiocarcinoma, kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, stomach adenocarcinoma, and uterine carcinosarcoma (all p < 0.05). On the other hand, samples from patients with colon adenocarcinoma, kidney chromophobe, lung adenocarcinoma, lung squamous cell carcinoma, pheochromocytoma and paraganglioma, prostate adenocarcinoma, and thyroid carcinoma showed reduced levels of EZR when compared to non-tumor tissues (all p < 0.05; Fig. 1). The role of EZR as an oncogene has been widely reported in the literature for several types of tumors, but there are few studies on its potential as a pharmacological target.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79. Thus, due to limited data on EZR in cervical cancer and stomach adenocarcinoma, the authors decided to further the clinical and biological characterization of the role of EZR in these cohorts, as well as to investigate the potential of a pharmacological ezrin inhibitor in cellular models of these tumor types.

Fig. 1
EZR expression in a pan-cancer analysis from The Cancer Genome Atlas (TCGA) cohorts. Data from TCGA cohorts that included 32 types of cancer, 21 of which had data from adjacent normal tissue, were obtained from the TIMER2.0 database (http://timer.cistrome.org/). Statistical significance was calculated using the Wilcoxon test; *p < 0.05, **p < 0.01, ***p < 0.001. Abbreviations: ACC, adrenocortical carcinoma; BLCA, Bladder Urothelial Carcinoma; BRCA, Breast Invasive Carcinoma; CESC, Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma; CHOL, Cholangiocarcinoma; COAD, Colon Adenocarcinoma; DLBC, Diffuse Large B-Cell lymphoma; ESCA, Esophageal Carcinoma; GBM, Glioblastoma Multiforme; HNSC, Head and Neck Squamous Cell Carcinoma; KICH, Kidney Chromophobe; KIRC, Kidney Renal Clear Cell Carcinoma; KIRP, Kidney Renal Papillary Cell Carcinoma; LAML, Acute Myeloid Leukemia; LGG, Brain Lower Grade Glioma; LIHC, Liver Hepatocellular Carcinoma; LUAD, Lung Adenocarcinoma; LUSC, Lung Squamous Cell Carcinoma; MESO, Mesothelioma; OV, Ovarian Serous Cystadenocarcinoma; PAAD, Pancreatic Adenocarcinoma; PCPG, Pheochromocytoma and Paraganglioma; PRAD, Prostate Adenocarcinoma; READ, Rectum Adenocarcinoma; SARC, Sarcoma; SKCM, Skin Cutaneous Melanoma; STAD, Stomach Adenocarcinoma; TGCT, Testicular Germ Cell Tumors; THCA, Thyroid Carcinoma; THYM, Thymoma; UCEC, Uterine Corpus Endometrial Carcinoma; UCS, Uterine Carcinosarcoma; UVM, Uveal Melanoma.

EZR expression is associated with relevant biological and molecular processes in cervical squamous carcinoma and stomach adenocarcinoma

In patients with cervical squamous cell carcinoma, EZR levels are not associated with disease staging (Fig. 2A), but there is a tendency for increased EZR expression to negatively impact disease-free survival (Hazard Ratio [HR = 2.65], 95% Confidence Interval [95 %CI 1.16‒6.09], p = 0.06) and disease-specific survival (HR = 1.91; 95% CI 1.06‒3.46, p = 0.06; Fig. 2B). Functional transcriptomic analysis identified differences in expression between patients with high and low expression of EZR (Fig. 2C‒D), as well as the correlation of EZR with several other genes (Fig. 2E). Those differently expressed genes were associated with the downregulation or upregulation of relevant biological processes, cellular components, and molecular functions, of which the authors would like to highlight cell adhesion, cell-cell signaling, development, cell mobility, migration, epithelial differentiation, and the EGFR-mediated signaling pathway (Fig. 2F-G). Similarly, EZR expression was not associated with staging (Fig. 3A) or survival outcomes (Fig. 3B) in the stomach adenocarcinoma cohort. High and low EZR expression allowed the identification of distinct gene expression signatures in patients with stomach adenocarcinoma (Fig. 3C‒E). Among the cellular processes, cellular components, and associated molecular functions, the authors would like to highlight cell-cell signaling, cell motility, migration, development, differentiation, and signaling receptor binding (Fig. 3F‒G). The association of EZR expression with clinical and laboratory characteristics in patients with cervical squamous carcinoma and stomach adenocarcinoma are described in Tables 1 and 2.

Fig. 2
Association of EZR expression and clinical, prognosis, biological and molecular characteristics in cervical squamous cell carcinoma TCGA cohort. (A) Box plot displaying EZR expression distribution between the different grades of cervical squamous cell carcinoma patients. (B) Survival curves were estimated using the Kaplan-Meier method. The Hazard Ratio (HR), 95% Confidence Interval, and p-values are indicated for overall survival, disease-free survival, and disease-specific survival. (C) Heatmap constructed using ClusterVis (https://biit.cs.ut.ee/clustvis/) displaying the top 25 upregulated and 25 downregulated DEGs across high vs. low EZR expression. Color intensity represents the z-score within each row. (D) Volcano plot depicting the extent (x-axis) and significance (y-axis) of DEGs, comparing high vs. low EZR expression. (E) Graph displaying Spearman’s correlation of EZR expression with other transcripts’ expression. Representative ShinyGO plots reflect the main 20 biological processes, cellular components, and molecular functions associated with downregulated genes (F) or upregulated genes (G) in high vs. low EZR cervical squamous cell carcinoma patients.

Fig. 3
Clinical, biological, and molecular significance of EZR expression in stomach adenocarcinoma TCGA cohort. (A) Box plot displaying EZR expression distribution between the different grades of stomach adenocarcinoma patients. (B) Survival curves were estimated using the Kaplan-Meier method. The Hazard Ratio (HR), 95% Confidence Interval, and p-values are indicated for overall survival, disease-free survival, and disease-specific survival. (C) Heatmap constructed using ClusterVis (https://biit.cs.ut.ee/clustvis/) displaying the top 25 upregulated and 25 downregulated DEGs across high vs low EZR expression. Color intensity represents the z-score within each row. (D) Volcano plot depicting the extent (x-axis) and significance (y-axis) of DEGs, comparing high vs. low EZR expression. (E) Graph displaying Spearman’s correlation of EZR expression with other transcripts’ expression. Representative ShinyGO plots reflect the main 20 biological processes, cellular components, and molecular functions associated with downregulated genes (F) or upregulated genes (G) in high vs. low EZR stomach adenocarcinoma patients.

Table 1
Association of EZR expression with clinical and molecular factors in TCGA cervical squamous cell carcinoma cohort.
Table 2
Association of EZR expression with clinical and molecular factors in TCGA stomach adenocarcinoma cohort.

Pharmacological inhibition of ezrin reduces viability and clonal growth in cervical and stomach cancer cells

Next, the potential antineoplastic effects of the pharmacological ezrin inhibitor NSC305787 were investigated in cervical (A-431) and gastric (AGS and HGC-27) cancer cells. NSC305787 dose-dependently reduced the viability of all tested cell models, with IC50 values of 1.6, 2.2, and 3.2 µM for A-431, AGS, and HGC-27 cells, respectively (Fig. 4A). Notably, long-term exposure strongly reduced colony formation at concentrations starting at 1.6 µM for A-431 and AGS cells and 2.4 µM for HGC-27 cells (all p < 0.05; Fig. 4B). Morphological analysis indicates the presence of cells with pyknotic or fragmented nuclei, in addition to the presence of actin inclusion bodies, suggesting cytoskeletal disorganization and apoptosis (Fig. 4C). Flow cytometry analyses indicate that exposure to NSC305787 induces cell death in a concentrationdependent manner in A-431, AGS, and HCG-27 cells (Fig. 5).

Fig. 4
Pharmacological inhibition of ezrin reduces cell viability and clonal growth of cervical and gastric cancer cells. (A) Dose-response cytotoxicity was evaluated with the Methylthiazoletetrazolium (MTT) assay. A-431, AGS, and HGC-27 cells were treated with vehicle or different concentrations of NSC305787 for 24h. Values are expressed as the percentage of viable cells for each condition relative to vehicletreated cells. Results are shown as mean ± SD of at least 3 independent experiments. (B) Colony formation of the cells treated with vehicle or NSC305787 for 10-15 days. The bar graph represents the mean ± SD of the relative number of colonies (% of control). * p < 0.05, *** p < 0.001; ANOVA and Bonferroni post-test. (C) Immunofluorescence analysis of A-431, AGS, and HGC-27 treated with vehicle or NSC305787 for 24h, displaying phalloidin (red) and DAPI (blue) staining. The white arrows indicate pyknotic nuclei and the yellow arrows indicate actin inclusion bodies. Scale bar 100 or 200 µm.

Fig. 5
NSC305787 induces cell death in cervical and gastric cancer cells. (A) Apoptosis was detected by flow cytometry in A-431, AGS, and HGC-27 cells treated with vehicle or increasing concentrations of NSC305787 for 24 hours using annexin V/Propidium Iodide (PI) staining. Representative dot plots are shown for each condition; the upper and lower right quadrants (Q2 + Q3) cumulatively contain the cell death population (annexin V+ cells). (B) The bar graphs represent the mean ± SD from at least three independent experiments quantifying cell death. The p-values and cell lines are indicated in the graphs; *p < 0.05, **p < 0.01, ***p < 0.0001; ANOVA test and Bonferroni post-test.

Pharmacological inhibition of ezrin induces molecular markers of cell death

From a molecular point of view, NSC305787 induces apoptosis markers (cleaved PARP1) and DNA damage (γH2AX) in the evaluated cervical and gastric cancer cell models (Fig. 6A). To provide new insights into ezrin inhibitor action in the studied models, a panel of genes involved in cell cycle progression, apoptosis, DNA damage, and autophagy was investigated. In A-431 cells, drug treatment reduced MCL1 and CCNE1 expression (p < 0.05). In AGS cells, an increase in BNIP3L, MCL1, PMAIP1, and SQSTM1 upon exposure to NSC305787 was observed (p < 0.05). The expression levels of ATG5, GADD45A, MCL1, and SQSTM1 were upregulated in HGC-27 cells treated with the pharmacological inhibitor of ezrin (p < 0.05; Fig. 6B).

Fig. 6
NSC305787 induces apoptosis and DNA damage markers and modulates survival-related genes in cervical and gastric cancer cells. (A) Western blot analysis PARP1 and γH2AX in total cell extracts from A-431, AGS, and HGC-27 cells treated with vehicle or NSC305787 for 24h. The membranes were incubated with the indicated antibodies and developed with the SuperSignal™ West Dura Extended Duration Substrate and Gel Doc XR + system. (B) Heatmap for the gene expression analysis in A-431, AGS, and HGC-27 cells treated with vehicle or NSC305787 for 24h. The data represent the fold change of vehicle-treated cells, and downregulated and upregulated genes are shown in blue and red, respectively. Fold-Change (FC), Standard Deviation (SD), and p-values are indicated, Student t-test.

Discussion

Increasing attention has been given to the functional characterization of the ezrin protein and its relationship with different types of cancer. Studies have highlighted its role in important cellular processes and the correlation between high levels of EZR expression and tumor progression and metastasis, as well as worse prognosis.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79.,1212 Brambilla D, Fais S. The Janus-faced role of ezrin in "linking" cells to either normal or metastatic phenotype. Int J Cancer 2009;125(10):2239–45. In cervical carcinoma, the ezrin protein is highly expressed in all cervical cancer cell lines and tissues compared with normal cervical tissues and predicts poor prognosis in cervical cancer patients.1313 Kong J, Di C, Piao J, Sun J, Han L, Chen L, et al. Ezrin contributes to cervical cancer progression through induction of epithelial-mesenchymal transition. Oncotarget 2016;7(15):19631–42.,1414 Kong J, Li Y, Liu S, Jin H, Shang Y, Quan C, et al. High expression of ezrin predicts poor prognosis in uterine cervical cancer. BMC Cancer 2013;13:520. EZR expression is also associated with the severity of HPV-associated squamous intraepithelial lesions.1515 Auvinen E, Carpen O, Korpela T, Ronty M, Vaheri A, Tarkkanen J. Altered expression of ezrin, E-Cadherin and beta-Catenin in cervical neoplasia. Neoplasma 2013;60(1):56–61. EZR silencing reduces migration, invasion, and mesenchymal markers expression, and PI3K/AKT activation in HeLa and SiHa cells.1313 Kong J, Di C, Piao J, Sun J, Han L, Chen L, et al. Ezrin contributes to cervical cancer progression through induction of epithelial-mesenchymal transition. Oncotarget 2016;7(15):19631–42.,1414 Kong J, Li Y, Liu S, Jin H, Shang Y, Quan C, et al. High expression of ezrin predicts poor prognosis in uterine cervical cancer. BMC Cancer 2013;13:520. Similarly, EZR was found to be highly expressed in gastric cancer compared to non-tumor mucosae, which is associated with age, tumor size, location, differentiation stage, depth of invasion, vessel invasion, lymph node, and distant metastasis, and it is an independent prognostic factor in patients with gastric carcinoma.1616 Li L, Wang YY, Zhao ZS, Ma J. Ezrin is associated with gastric cancer progression and prognosis. Pathol Oncol Res 2011;17(4):909–15.,1717 Jin J, Jin T, Quan M, Piao Y, Lin Z. Ezrin overexpression predicts the poor prognosis of gastric adenocarcinoma. Diagn Pathol 2012;7(135). Furthermore, increased EZR expression is related to Helicobacter pylori infection, a risk factor for gastric cancer.1818 Bal N, Yildirim S, Nursal TZ, Bolat F, Kayaselcuk F. Association of ezrin expression in intestinal and diffuse gastric carcinoma with clinicopathological parameters and tumor type. World J Gastroenterol 2007;13(27):3726–9. The present findings corroborate these previous studies, indicating that the high expression of EZR may contribute to the malignant phenotype of these neoplasms.

Among the main biological processes and molecular functions, the authors highlight that a higher expression of EZR, both in cervical carcinoma and gastric adenocarcinoma, is related to the reduction of biological processes associated with signaling and cell-cell contact, negatively regulating cell junction, a factor that can facilitate the development of metastasis. Among the biological processes with increased function in patients with high EZR expression, there is cell motility, locomotion, migration, and epithelial differentiation, all factors strongly related to tumor formation. The increased molecular functions, on the other hand, highlight a strong association with growth factors, also related to tumor development and association with the binding of actin filaments, a function of ezrin already well described in the literature.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79.,1919 Barik GK, Sahay O, Paul D, Santra MK. Ezrin gone rogue in cancer progression and metastasis: An enticing therapeutic target. Biochim Biophys Acta Rev Cancer 2022;1877(4):188753.

Given the above, EZR proved to be a potential therapeutic target and inhibitor molecules began to be tested in preclinical studies for various tumors.33 Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79.,1919 Barik GK, Sahay O, Paul D, Santra MK. Ezrin gone rogue in cancer progression and metastasis: An enticing therapeutic target. Biochim Biophys Acta Rev Cancer 2022;1877(4):188753. Two pharmacological ezrin inhibitors have been better characterized in cancer models, NSC305787 and NSC668394, with NSC305787 being considered more potent and with a more favorable pharmacokinetic profile,66 Bulut G, Hong SH, Chen K, Beauchamp EM, Rahim S, Kosturko GW, et al. Small molecule inhibitors of ezrin inhibit the invasive phenotype of osteosarcoma cells. Oncogene 2012;31(3):269–81.,77 Celik H, Bulut G, Han J, Graham GT, Minas TZ, Conn EJ, et al. Ezrin Inhibition Up-regulates Stress Response Gene Expression. J Biol Chem 2016;291(25):13257–70. which led us to choose this compound for the studies in vitro in models of cervical and gastric cancer. NSC305787 reduced cell viability, and colony formation, and led to morphological changes compatible with apoptosis. From a molecular point of view, drug treatment induced markers of DNA damage and apoptosis. These findings are in line with those previously described in pancreatic cancer and acute leukemia.2020 Lipreri da Silva JC, Carvalho MFL, de Miranda LBL, de Almeida BO, Lima K, Machado-Neto JA. NSC305787, a pharmacological ezrin inhibitor, exhibits antineoplastic activity in pancreatic cancer cells. Invest New Drugs 2022;40(4):728–37.,2121 Lipreri da Silva JC, Coelho-Silva JL, Lima K, Vicari HP, Lazarini M, Costa-Lotufo LV, et al. Comprehensive analysis of cytoskeleton regulatory genes identifies ezrin as a prognostic marker and molecular target in acute myeloid leukemia. Cell Oncol (Dordr) 2021;44(5):1105–17.

In exploratory analysis of gene expression, pharmacological inhibition reduced MCL1 and CCNA2 expression in A-431 cells. MCL1 is a BCL2 family protein that can act in the control of apoptosis, interfering at an early stage in a cascade of events that lead to the release of cytochrome c from mitochondria.2222 Michels J, Johnson PW, Packham G. Mcl-1. Int J Biochem Cell Biol 2005;37(2):267–71. The relationship between MCL1 and cancer is well described in the literature, studies have shown that its overexpression contributes to cell survival and resistance to various chemotherapeutic agents.2323 Pratelli G, Carlisi D, Di Liberto D, Notaro A, Giuliano M, D’Anneo A, et al. MCL1 Inhibition Overcomes the Aggressiveness Features of Triple-Negative Breast Cancer MDA-MB-231 Cells. Int J Mol Sci 2023;24(13):11149.,2424 Wang L, Xi C, Liu R, Ye T, Xiang N, Deng J, et al. Dual targeting of Mcl-1 and Bcl-2 to overcome chemoresistance in cervical and colon cancer. Anticancer Drugs 2024;35(3):219–26. Research has shown that there is overexpression of MCL1 in cervical cancer, while adjacent normal tissue showed slight or no expression of this protein. Furthermore, increased levels of MCL1 are closely associated with the pathogenesis of cervical cancer, being associated with histological grade, tumor size, and lymph node involvement, and positively correlated with poor prognosis.2525 Zhang T, Zhao C, Luo L, Zhao H, Cheng J, Xu F. The expression of Mcl-1 in human cervical cancer and its clinical significance. Med Oncol 2012;29(3):1985–91. In contrast, the MCL1 gene was significantly more expressed in both gastric cancer cell lines, which could act as a drug resistance mechanism.2626 Li S, Guo W, Wu H. The role of post-translational modifications in the regulation of MCL1. Cell Signal 2021;81:109933. Cyclin A2 is associated with both CDK1 and CDK2 and has roles in both S phase and mitosis.2727 Yam CH, Fung TK, Poon RY. Cyclin A in cell cycle control and cancer. Cell Mol Life Sci. 2002;59(8):1317–26. Consistent with its role as a key cell cycle regulator, expression of this protein is elevated in a variety of tumors.2828 Zhang X, Yang L, Huang B, Yin J, Wei Y. Identification and validation of cyclin A2 and cyclin E2 as potential biomarkers in small cell lung cancer. Oncol Res Treat 2023;46(6):246–58.,2929 Walach MT, Nitschke K, Gross-Weege M, Grosshans J, Wildner L, Pause L, et al. Cyclin A2 Expression as Predictive Biomarker in Muscle-Invasive Upper Tract Urothelial Carcinoma. Urol Int 2024;108(2):128–36. In addition to potentially disrupting the cell cycle and checkpoints directly, cyclin A may also contribute to tumor formation through the phosphorylation of other oncogenic and tumor suppressor proteins.2727 Yam CH, Fung TK, Poon RY. Cyclin A in cell cycle control and cancer. Cell Mol Life Sci. 2002;59(8):1317–26. Research has shown that CCNA2 is significantly increased in cervical cancer tissues, suggesting a possible mechanism of anomalous mitosis in this condition.3030 Cheng J, Lu X, Wang J, Zhang H, Duan P, Li C. Interactome analysis of gene expression profiles of cervical cancer reveals dysregulated mitotic gene clusters. Am J Transl Res 2017;9(6):3048–59. The increase in this expression reflects a high rate of cell proliferation once the tumor has been developed.3131 Kanai M, Shiozawa T, Xin L, Nikaido T, Fujii S. Immunohistochemical detection of sex steroid receptors, cyclins, and cyclin-dependent kinases in the normal and neoplastic squamous epithelia of the uterine cervix. Cancer 1998;82(9):1709–19.

In the present study, gastric cancer cells showed increased expression of genes involved in the autophagy process, such as ATG5 and SQSTM1, after exposure to NSC305787. ATG5 encodes an E1-type activating enzyme essential for autophagy and for the transport of molecules from the cytoplasm to the vacuole,3232 Arakawa S, Honda S, Yamaguchi H, Shimizu S. Molecular mechanisms and physiological roles of Atg5/Atg7-independent alternative autophagy. Proc Jpn Acad Ser B Phys Biol Sci 2017;93(6):378–85. while SQSTM1 encodes the p62 protein, which acts as a receptor for the degradation of ubiquitinated proteins by autophagic or proteasome pathways.3333 Lamark T, Svenning S, Johansen T. Regulation of selective autophagy: the p62/SQSTM1 paradigm. Essays Biochem 2017;61(6):609–24. These findings indicate that EZR inhibition in gastric cancer cell lines may activate genes related to autophagy, which deserves further investigation. Furthermore, in the AGS cells, there was an increase in the expression of BNIP3L, which was initially recognized as coding for a pro-apoptotic protein of the BCL2 family. Specifically, BNIP3L interacts with the pro-apoptotic proteins of the BAX and BAK family to increase the permeability of the outer mitochondrial membrane. However, the most important activity of BNIP3L in the cell has been found to be the promotion of mitochondrial autophagy (mitophagy) through the recruitment of autophagosomes in response to cellular or environmental stress.3434 Li Y, Zheng W, Lu Y, Zheng Y, Pan L, Wu X, et al. BNIP3L/NIX-mediated mitophagy: molecular mechanisms and implications for human disease. Cell Death Dis 2021;13(1):14.,3535 Lu JJ, Wu PF, He JG, Li YK, Long LH, Yao XP, et al. BNIP3L/NIX-mediated mitophagy alleviates passive stress-coping behaviors induced by tumor necrosis factor-alpha. Mol Psychiatry 2023;28(12):5062–76. Also in this cell line, there was a decrease in the expression of the PMAIP1, which belongs to the pro-apoptotic subfamily within the BCL2 protein family and promotes the activation of caspases, changes in the mitochondrial membrane and efflux of apoptogenic proteins from the mitochondria.3636 Ishida M, Sunamura M, Furukawa T, Lefter LP, Morita R, Akada M, et al. The PMAIP1 gene on chromosome 18 is a candidate tumor suppressor gene in human pancreatic cancer. Dig Dis Sci 2008;53(9):2576–82.,3737 Zhao Y, Zhao C, Deng Y, Pan M, Mo G, Liao Z, et al. PMAIP1 promotes J subgroup avian leukosis virus replication by regulating mitochondrial function. Poult Sci 2024;103(6):103617.

Additionally, in HGC-27 cells, there was an increase in the expression of BAD and GADD45A. BAD encodes a protein member of the BCL2 family, which are regulators of programmed cell death. This protein positively regulates cell apoptosis by forming heterodimers with BCL2L1 and BCL2, reversing its death repressor activity.3838 Howells CC, Baumann WT, Samuels DC, Finkielstein CV. The Bcl-2-associated death promoter (BAD) lowers the threshold at which the Bcl-2-interacting domain death agonist (BID) triggers mitochondria disintegration. J Theor Biol 2011;271(1):114–23. GADD45A is a member of a group of genes whose transcription levels are increased after stressful conditions of growth arrest and DNA damage.3939 Zhan Q. Gadd45a, a p53- and BRCA1-regulated stress protein, in cellular response to DNA damage. Mutat Res 2005;569(1-2):133–43.,4040 Palomer X, Salvador JM, Grinan-Ferre C, Barroso E, Pallas M, Vazquez-Carrera M. GADD45A: With or without you. Med Res Rev 2024;44(4):1375–403.

In summary, this study reveals elevated EZR expression in cervical and gastric carcinoma, highlighting its potential as a molecular target for treating these cancers. The present preclinical findings pave the way for further research and support the need for clinical trials to explore the effectiveness of pharmacological EZR inhibitors in patients with cervical and gastric carcinoma.

Data available

All the data are available from the corresponding author upon request.

Acknowledgments

The authors would like to acknowledge all the research participants contributing to The Cancer Genome Atlas (TCGA) resource for providing high-quality data for analyses.

Supplementary materials

Supplementary material associated with this article can be found in the online version at doi:10.1016/j.clinsp.2024.100422.

  • Funding
    M.F.L.C., B.O.A., L.B.L.M., J.C.L.S., and K.L. received a fellowship from the São Paulo Research Foundation (FAPESP) (grants #2022/14505-6, #2022/11038-8, #2022/03316-8, #2020/12909-7, and #2020/12842-0). This study was supported by grants #2019/23864-7, #2021/11606-3, and #2023/12246-6 from FAPESP. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brasil (CAPES), Finance Code 001.
  • Ethical approval
    Not applicable.

References

  • 1
    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov 2012;2(5):401–4.
  • 2
    Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res 2020;48(W1):W509–14.
  • 3
    Lipreri Da Silva JC, Vicari HP, Machado-Neto JA. Perspectives for Targeting Ezrin in Cancer Development and Progression. Future pharmacology 2023;3(1):61–79.
  • 4
    Xie JJ, Zhang FR, Tao LH, Lu Z, Xu XE, Jian S, et al. Expression of ezrin in human embryonic, fetal, and normal adult tissues. J Histochem Cytochem 2011;59(11):1001–8.
  • 5
    Xu WT, Shi LL, Xu J, Qian H, Zhou H, Wang LH. Ezrin expression in female reproductive tissues: A review of regulation and pathophysiological implications. Front Cell Dev Biol 2023;11:1125881.
  • 6
    Bulut G, Hong SH, Chen K, Beauchamp EM, Rahim S, Kosturko GW, et al. Small molecule inhibitors of ezrin inhibit the invasive phenotype of osteosarcoma cells. Oncogene 2012;31(3):269–81.
  • 7
    Celik H, Bulut G, Han J, Graham GT, Minas TZ, Conn EJ, et al. Ezrin Inhibition Up-regulates Stress Response Gene Expression. J Biol Chem 2016;291(25):13257–70.
  • 8
    TCGA. Integrated genomic and molecular characterization of cervical cancer. Nature 2017;543(7645):378–84.
  • 9
    TCGA. Comprehensive molecular characterization of gastric adenocarcinoma. Nature 2014;513(7517):202–9.
  • 10
    Heagerty PJ, Zheng Y. Survival model predictive accuracy and ROC curves. Biometrics 2005;61(1):92–105.
  • 11
    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001;25(4):402–8.
  • 12
    Brambilla D, Fais S. The Janus-faced role of ezrin in "linking" cells to either normal or metastatic phenotype. Int J Cancer 2009;125(10):2239–45.
  • 13
    Kong J, Di C, Piao J, Sun J, Han L, Chen L, et al. Ezrin contributes to cervical cancer progression through induction of epithelial-mesenchymal transition. Oncotarget 2016;7(15):19631–42.
  • 14
    Kong J, Li Y, Liu S, Jin H, Shang Y, Quan C, et al. High expression of ezrin predicts poor prognosis in uterine cervical cancer. BMC Cancer 2013;13:520.
  • 15
    Auvinen E, Carpen O, Korpela T, Ronty M, Vaheri A, Tarkkanen J. Altered expression of ezrin, E-Cadherin and beta-Catenin in cervical neoplasia. Neoplasma 2013;60(1):56–61.
  • 16
    Li L, Wang YY, Zhao ZS, Ma J. Ezrin is associated with gastric cancer progression and prognosis. Pathol Oncol Res 2011;17(4):909–15.
  • 17
    Jin J, Jin T, Quan M, Piao Y, Lin Z. Ezrin overexpression predicts the poor prognosis of gastric adenocarcinoma. Diagn Pathol 2012;7(135).
  • 18
    Bal N, Yildirim S, Nursal TZ, Bolat F, Kayaselcuk F. Association of ezrin expression in intestinal and diffuse gastric carcinoma with clinicopathological parameters and tumor type. World J Gastroenterol 2007;13(27):3726–9.
  • 19
    Barik GK, Sahay O, Paul D, Santra MK. Ezrin gone rogue in cancer progression and metastasis: An enticing therapeutic target. Biochim Biophys Acta Rev Cancer 2022;1877(4):188753.
  • 20
    Lipreri da Silva JC, Carvalho MFL, de Miranda LBL, de Almeida BO, Lima K, Machado-Neto JA. NSC305787, a pharmacological ezrin inhibitor, exhibits antineoplastic activity in pancreatic cancer cells. Invest New Drugs 2022;40(4):728–37.
  • 21
    Lipreri da Silva JC, Coelho-Silva JL, Lima K, Vicari HP, Lazarini M, Costa-Lotufo LV, et al. Comprehensive analysis of cytoskeleton regulatory genes identifies ezrin as a prognostic marker and molecular target in acute myeloid leukemia. Cell Oncol (Dordr) 2021;44(5):1105–17.
  • 22
    Michels J, Johnson PW, Packham G. Mcl-1. Int J Biochem Cell Biol 2005;37(2):267–71.
  • 23
    Pratelli G, Carlisi D, Di Liberto D, Notaro A, Giuliano M, D’Anneo A, et al. MCL1 Inhibition Overcomes the Aggressiveness Features of Triple-Negative Breast Cancer MDA-MB-231 Cells. Int J Mol Sci 2023;24(13):11149.
  • 24
    Wang L, Xi C, Liu R, Ye T, Xiang N, Deng J, et al. Dual targeting of Mcl-1 and Bcl-2 to overcome chemoresistance in cervical and colon cancer. Anticancer Drugs 2024;35(3):219–26.
  • 25
    Zhang T, Zhao C, Luo L, Zhao H, Cheng J, Xu F. The expression of Mcl-1 in human cervical cancer and its clinical significance. Med Oncol 2012;29(3):1985–91.
  • 26
    Li S, Guo W, Wu H. The role of post-translational modifications in the regulation of MCL1. Cell Signal 2021;81:109933.
  • 27
    Yam CH, Fung TK, Poon RY. Cyclin A in cell cycle control and cancer. Cell Mol Life Sci. 2002;59(8):1317–26.
  • 28
    Zhang X, Yang L, Huang B, Yin J, Wei Y. Identification and validation of cyclin A2 and cyclin E2 as potential biomarkers in small cell lung cancer. Oncol Res Treat 2023;46(6):246–58.
  • 29
    Walach MT, Nitschke K, Gross-Weege M, Grosshans J, Wildner L, Pause L, et al. Cyclin A2 Expression as Predictive Biomarker in Muscle-Invasive Upper Tract Urothelial Carcinoma. Urol Int 2024;108(2):128–36.
  • 30
    Cheng J, Lu X, Wang J, Zhang H, Duan P, Li C. Interactome analysis of gene expression profiles of cervical cancer reveals dysregulated mitotic gene clusters. Am J Transl Res 2017;9(6):3048–59.
  • 31
    Kanai M, Shiozawa T, Xin L, Nikaido T, Fujii S. Immunohistochemical detection of sex steroid receptors, cyclins, and cyclin-dependent kinases in the normal and neoplastic squamous epithelia of the uterine cervix. Cancer 1998;82(9):1709–19.
  • 32
    Arakawa S, Honda S, Yamaguchi H, Shimizu S. Molecular mechanisms and physiological roles of Atg5/Atg7-independent alternative autophagy. Proc Jpn Acad Ser B Phys Biol Sci 2017;93(6):378–85.
  • 33
    Lamark T, Svenning S, Johansen T. Regulation of selective autophagy: the p62/SQSTM1 paradigm. Essays Biochem 2017;61(6):609–24.
  • 34
    Li Y, Zheng W, Lu Y, Zheng Y, Pan L, Wu X, et al. BNIP3L/NIX-mediated mitophagy: molecular mechanisms and implications for human disease. Cell Death Dis 2021;13(1):14.
  • 35
    Lu JJ, Wu PF, He JG, Li YK, Long LH, Yao XP, et al. BNIP3L/NIX-mediated mitophagy alleviates passive stress-coping behaviors induced by tumor necrosis factor-alpha. Mol Psychiatry 2023;28(12):5062–76.
  • 36
    Ishida M, Sunamura M, Furukawa T, Lefter LP, Morita R, Akada M, et al. The PMAIP1 gene on chromosome 18 is a candidate tumor suppressor gene in human pancreatic cancer. Dig Dis Sci 2008;53(9):2576–82.
  • 37
    Zhao Y, Zhao C, Deng Y, Pan M, Mo G, Liao Z, et al. PMAIP1 promotes J subgroup avian leukosis virus replication by regulating mitochondrial function. Poult Sci 2024;103(6):103617.
  • 38
    Howells CC, Baumann WT, Samuels DC, Finkielstein CV. The Bcl-2-associated death promoter (BAD) lowers the threshold at which the Bcl-2-interacting domain death agonist (BID) triggers mitochondria disintegration. J Theor Biol 2011;271(1):114–23.
  • 39
    Zhan Q. Gadd45a, a p53- and BRCA1-regulated stress protein, in cellular response to DNA damage. Mutat Res 2005;569(1-2):133–43.
  • 40
    Palomer X, Salvador JM, Grinan-Ferre C, Barroso E, Pallas M, Vazquez-Carrera M. GADD45A: With or without you. Med Res Rev 2024;44(4):1375–403.

Publication Dates

  • Publication in this collection
    16 Aug 2024
  • Date of issue
    2024

History

  • Received
    01 Feb 2024
  • Reviewed
    13 May 2024
  • Accepted
    10 June 2024
Faculdade de Medicina / USP Rua Dr Ovídio Pires de Campos, 225 - 6 and., 05403-010 São Paulo SP - Brazil, Tel.: (55 11) 2661-6235 - São Paulo - SP - Brazil
E-mail: clinics@hc.fm.usp.br