Downregulation of GLYAT facilities tumor growth and metastasis and poor clinical outcomes through the PI3K/AKT/Snail pathway in human breast cancer

doi:10.21203/rs.3.rs-125337/v1

Download PDF

Research

Downregulation of GLYAT facilities tumor growth and metastasis and poor clinical outcomes through the PI3K/AKT/Snail pathway in human breast cancer

https://doi.org/10.21203/rs.3.rs-125337/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Background

The Glycine N-acyltransferase (GLYAT) gene encodes a protein that catalyzes the transfer of acyl groups from acyl CoA to glycine, resulting in acyl glycine and coenzyme A. Aberrant GLYAT expression is associated with several malignant tumors, but its clinical importance in malignant tumors, especially human breast cancer (BC), has yet to be fully addressed. This study aims to evaluate the clinical function of GLYAT in BC patients.

Methods

GLYAT expression was determined by immune blot and immunohistochemistry in three BC cell lines and primary cancer tissues. The MDA-MB 231 cell line was used for GLYAT gene knockdown experiments while the MCF7 cell line for overexpression experiments. Colony formation experiments, soft agar experiments, and transwell assays were utilized for further inspection of cell migration and proliferation capabilities. Immunofluorescence and western blot were used to detect markers of the epithelial-mesenchymal transition (EMT) and changes in the PI3K/AKT/Snail pathway. The role of GLYAT in tumor growth and metastasis was also assessed in nude mice in vivo. Also, a correlation analysis was performed between clinicopathological features and GLYAT expression in BC patients.

Results

GLYAT was decreased in human BC tissues and cell lines. Functional analysis showed that knockdown of GLYAT augmented BC cell proliferation in vitro and in vivo. However, this phenomenon was reversed when GLYAT was overexpressed in the transfected cells. Moreover, GLYAT inhibited the migratory properties of BC cells, likely through the activation of PI3K/ATK/Snail signaling, which subsequently induced the EMT. IHC analysis indicated that GLYAT was decreased in human BC tissues and lower GLYAT expression was correlated with histological grade, tumor TNM stage, Ki-67 status, and poorer survival in BC patients. Furthermore, lower GLYAT expression seemed as an independent risk factor related to poor prognosis in BC patients based on Cox regression analyses.

Conclusions

Our findings demonstrate that downregulation of GLYAT expression in human breast cancer is correlated with EMT via the PI3K/ATK/Snail pathway and is also associated with histological grade, tumor TNM stage, Ki-67 status, and poor survival in breast cancer patients.

Clinical Pharmacology

Cancer Biology

Oncology

GLYAT

breast cancer

EMT

PI3K/ATK

clinicopathological features

prognosis

Breast cancer (BC) amongst women has been on the rise and it is fast becoming the most common malignancy in this population [1]. Leading cancer statistics show that the incidence and mortality of breast cancer remain high worldwide [2, 3]. Although surgery, chemoradiotherapy, endocrine therapy, and molecular targeted therapy have greatly improved the survival of BC patients, those with advanced BC are often left with very limited therapeutic options[4, 5]. Many factors influence the prognosis and probability of responding to systemic therapies (e.g., TNM stage, ER, PR, and HER2 status), but clinically, we have observed that even patients with the same TNM stage, molecular typing, and other diagnostic markers still experience varied prognoses despite being given the same treatment modalities. While early diagnosis and initiation of treatment are crucial in improving survival rates of this debilitating disease, there is an urgent need for a deeper understanding of its molecular biology in order to aid in the development of more personalized and targeted treatment.

The EMT involves the dedifferentiation of epithelial cells into mesenchymal cells [6, 7]. Cells which have undergone EMT demonstrate more malignant features that promote metastasis, treatment resistance, as well as relapse [8–10]. The EMT participates in tumor invasion and metastasis through multiple pathways. Many studies have indicated the EMT is strongly linked to breast cancer development [11–13]. The EMT is regulated by multiple factors, among which the primary signaling pathway is the PI3K/AKT pathway [14].

The glycine N-acyltransferase (GLYAT) gene was first discovered in the mitochondria of the bovine liver in 1953 and isolated in the human liver and kidneys in 1976 [15]. It contains more than 23,000 base pairs over six exons and is situated chromosome 11 at position 11q12 [16]. GLYAT catalyzes acyl group transfer from acyl CoA to glycine, resulting in acyl glycine and coenzyme A (acyl-CoA) [17]. Several catabolic and anabolic reactions are catalyzed by acyl-CoA esters [18]. Almost all catabolic reactions produce acyl-CoA, the product of which is an important source of oxidative phosphorylation and lipogenesis. Studies have shown that GLYAT expression is suppressed in human hepatocellular carcinomas and may be a critical molecule in the transition between the differentiation and carcinogenesis of liver cells [19]. Nevertheless, literature is scarce surrounding GLYAT expression and its impact on human breast cancer.

The current investigation uncovered markedly suppressed GLYAT expression in breast cancer cells and tissues, which correlated with poorer prognosis and highly malignant clinicopathologic features in individuals with breast cancer. Interestingly, the EMT pathway appeared to be augmented in the presence of GLYAT, working to decrease in vivo and in vitro breast cancer cell migration through alteration of the PI3K/ATK/Snail signaling pathway. In conclusion, our study is the first of its kind to implicate GLYAT to be a breast cancer anti-oncogene. These findings lay the foundation for future research on the biological behaviour and targeted therapy of BC.

Datasets and GLYAT expression analysis

The public data used in this study were obtained from UALCAN dataset (http://ualcan.path.uab.edu/index.html) [20], GEPIA (http://gepia.cancer-pku.cn/index.html) [21], and the Human Protein Atlas dataset (http://www.proteinatlas.org) [22]. For UALCAN data, we analyzed the heat map of GLYAT expression between normal and breast cancer samples; For GEPIA data, we scrutinized differently GLYAT expression level and survival data in breast cancer. The relationship between patient survival and GLYAT levels were also analyzed using information available from the Human Protein Atlas data.

Cell lines

Human BC cell lines MDA-MB-231, MCF-7, and SKRB-3 were supplied by The Stem Cell Bank, Chinese Academy of Sciences (Shanghai, China). Cells were maintained at 37°C in a humidified cubicle containing 5% CO₂and 10% fetal bovine serum (FBS) in DMEM (all from Gibco, Carlsbad, CA, USA).

Plasmidformation of GLYAT knockdown and overexpression in breast cancer cells

The pGPU6/mCherry/Puro-shRNA-GLYAT plasmid was purchased from GenePharma Company (Shanghai, China). The pOGP-T2A-CKNeo-GLYAT overexpression plasmid was

purchased from JTS Scientific Company (Wuhan, China). The target sequences of GLYAT were indicated in (Additional file 1: Table S1). Lipofectamine 3000 (Thermo Fisher Scientific, US) was used to transfect plasmids into cells as reported by the manufacturer’s instructions. Cells transfected with GPU6/mCherry/Puro-shNC or pOGP-T2A-CKNeo-NC were used as controls. Transfected cells were screened using 3 ug/mL puromycin or 400 ug/mL G418 for 2 weeks. Stable cells were cultured in complete medium with 0.25 ug/mL puromycin or 100 ug/mL G418 (Beyotime, Nanjing, China). Positive clones were then selected and amplified for further analyses.

Cell proliferation assay

For colony forming experiments, six-well plates were used to house stable GLYAT-shRNA and NC cells. The cells were maintained at a concentration of 500 cells/well for 12 days in DMEM supplemented with 10% FBS, with media changed once every three days. Cells were then methanol-fixed and treated with crystal violet (Sigma, St Louis, MO, USA) prior to manual counting and photographing of the visible colonies. For soft agar experiments, Cells are harvested and pipetted well to become single-cell suspension in complete culture media in 1x 106 /ml. A mixture of 0.9 ml 4% soft-agar (Sigma) with 4.1 ml pre-warmed 10% FBS DMEM was added into a 60-mm culture dish to make the bottom layer. The top layer contained 3 x 104 cells in 3 ml of 10% FBS DMEM and 0.36% agar. The soft-agar colony dish was marked and placed at a 37 °C incubator for 3 weeks. cells were harvested and fully mixed to form a single-cell suspension in media of 1 x 10⁶/ml. The upper layer included 3 x 10⁴ cells, 3 ml of 10% fetal bovine serum, DMEM, combined with 0.36% agar (Sigma), while the bottom layer contained 0.9 ml 4% soft-agar, 4.1 ml pre-heated 10% fetal bovine serum, and DMEM in a 60-mm culture dish. Labeled colony dishes were placed in an incubator at 37°C for 3 weeks.

Transwell migration assay

Costar chambers containing 8 μm pore inserts were used (Millipore-Sigma, Danvers, MA, USA). The top chamber was used to house suspended cells transfected with GLYAT KD, GLYAT OE, or NC in 200 μl of serum-free medium, while media with 20% FBS was applied to the bottom chamber. The upper chamber was removed after 36 hours of incubation, with the bottom cartridge fixed with methanol and DAPI-stained to determine the migrated number of cells using a microscope to visualize the cells across five randomly selected fields (Olympus, Tokyo, Japan).

Western blot analysis

The proteins in BC cells were first extracted and separated. Skimmed milk (5%) was used to block endogenous reactions. Membranes were left to incubate with antibodies of GLYAT（1:1000，Abcam, Cambridge, MA, USA）, E-cadherin, vimentin, N-cadherin, fibronectin (1:1000 dilution respectively, Cell Signaling Technology, Boston, MA, USA), p-AKT Ser473, AKT (1:1000 and 1:5000 dilution, Signalway Antibody, Baltimore, MD, USA), PI3K p58 and SNAI1 (1:10000 and 1:1000 dilution, Proteintech, Rosemont, IL, USA) in agreement with the protocols set by the manufacturer. Anti-β-actin (1:10000 dilution, Proteintech, Rosemont, IL, USA) was used as the internal control. After three rinses with TBS-T for 5 minutes, the membranes were incubated with the corresponding secondary antibody and observed with an ECL Plus kit.

Immunofluorescence staining

Cells were subjected into 8 well chamber slides (Millipore-Sigma, Danvers, MA, USA) and for 20 hours incubation then rinsed with PBS containing 10% FBS. Pre-cooled methanol was used to fix the cells prior to further incubation with 0.2% Triton. BSA (5%) was used to block cells. The cells were next incubated overnight with primary antibodies E-cadherin, vimentin (1:200 dilution respectively) before being washed and re-incubated with fluorescein-conjugated secondary antibodies (1:500 ZSGB-Bio, Beijing China) for one hour. Then, DAPI and a fluorescence microscope (Olympus, Tokyo, Japan) were used for visualization and intensity analysis. For dewaxed sections, the process was similar like above, and primary antibodies dilution of E-cadherin, vimentin and p-AKT Ser473 were 1:100 respectively.

In vivo metastasis assay

The HFK BIOSCIENCE CO., LTD (Beijing, China) provided 5 weeks old BALB/c (nu/nu) female nude mice (18–20 g). These animals were divided into four cohorts with five mice each. Food and drinking water were provided every day. All mice received subcutaneous injections into their left breast fat pads with 100 microliters of PBS containing 2.0 × 10⁶ MCF7 cells with or without GLYAT OE and MDA-MB-231 cells with or without GLYAT KD. At time of MCF7 cells injection, E2 pellets (60-day release, 1.5 mg/pellet; Innovative Research of America,Sarasota, Florida USA) was implanted subcutaneously at the mammary fat pad. Four days measurements and volume calculation ((mm³) =length × width²/2) of the tumors were performed. The mice were sacrificed and underwent tumor removal after 24 days. The tumors were processed with 4% paraformaldehyde solution for immunofluorescence staining analysis.

Patients and samples

Two independent breast cancer patient cohorts from the Shengjing Hospital of CMU were collected for this investigation. The first cohort consisted of 21 breast cancer patients (seven for Luminal A, seven for Luminal B, and seven for TNBC) and 20 normal controls. Samples from the cancer patients and normal controls were analyzed for GLYAT expression. The second cohort comprised of 310 chemoradiotherapy-naïve breast cancer patients who had distinctive pathological diagnoses and possessed complete follow-up data from 2006 to 2008. Additional follow-up of these patients was carried out until September 2012.

Immunohistochemistry

The dewaxed sections were retrieved with Tris/EDTA (pH9.0) and quenched with 0.3% H₂O₂. Sections then underwent an overnight incubation with GLYAT primary antibody (1:100，Abcam, Cambridge, MA, USA). PBS was used to rinse the samples prior to repeat incubation with secondary biotinylated antibodies at room temperature for 45 minutes. Coloration was visualized with 3,3 '-diaminobenzidine tetrachoric acid (Sigma-Aldrich), Vector hematoxylin QS (Vector Laboratories) was reverse-stained, Zeiss Mirax MidiSlide scanner was analyzed, and image collection was performed with 3 CCD color cameras and Panoramic Audience (3DHISTECH, Budapest Hungary). H-score was calculated as the percentage of positive tumor cells multiplied by the intensity of staining (0, no staining; 1, weak; 2, mild to moderate; 3, strong staining) and the overall H-score ranged from 0 to 300. An H-score of more than 100 was set as cut-off value to define high/low level for GLYAT expression.

Statistical analysis

Data were analyzed using SPSS version 25.0 (SPSS, Inc., Chicago, IL, USA) and expressed as means ± SD. The Kaplan-Meier method was used to estimate survival curves. Cox regression model was applied to carry out univariate and multivariate statistical analysis. Comparison between two groups was assessed by Chi-square. P value of less than 0.05 was considered statistically significant.

Characterization of GLYAT expression profiling in human BC

The heat map analysis of the UALCAN dataset showed lower expression of GLYAT protein and mRNA in BC tissues in contrast to healthy samples (Fig. 1a and 1c). However, the expression of GLYAT was not statistically significant change in other cancer types, besides BRCA, LIHC and CHOL (Additional file 2: Figure S1), which indicated that GLYAT may have specific function in breast cancer. Additionally, analysis of data published in another study in Nature [23] which compared mRNA expression levels of GLYAT in normal breast tissues with invasive ductal and invasive lobular carcinoma revealed that two breast carcinoma subtypes had lower GLYAT mRNA levels (Fig. 1b). A Kaplan-Meier analysis from GEPIA dataset and a scatter plot from the Human Protein Atlas dataset suggested that patients with low GLYAT expression in breast cancer tissues had worse overall survival (OS) (Fig. 1d-e), further highlighting the suppressive function of GLYAT in BC progression.

GLYAT suppresses BC cell proliferation and metastasis

Three typical breast cancer cell lines were scrutinized for GLYAT expression. These cell lines included SK-BR3 (Her-2 positive subtype), MDA-MB-231 (triple negative subtype) and MCF-7 (luminal A subtype). Western blotting showed that the MDA-MB-231 cell line had the highest GLYAT expression, while the MCF-7 cell line had the lowest amount (Fig. 2a). Therefore, these two cell lines above were selected for further experiments.

In vitro assays were carried out in MDA-MB-231 cells with GLYAT knockdown (KD) and overexpression (OE) in MCF-7 cells to determine the impact of GLYAT expression on breast cancer behavior (Fig. 2b-c). The tablet colony formation and soft agar colony formation experiments revealed that GLYAT KD in MDA-MB-231 cells significantly increased colony numbers in contrast to control cell lines (Fig. 2d-e). Transwell migration assay revealed artificially altered MDA-MB-231 cells with GLYAT KD also had stronger migratory abilities in contrast to the control cells (Fig. 2f). In line with these results, the converse was seen in MCF-7 cells with OE GLYAT (Fig. 2g). These results indicate that GLYAT may hinder the proliferative and metastatic abilities of breast cancer cells.

GLYAT suppresses the EMT phenotype in BC cells

Given the prominent role of EMT in tumor spread, we postulated that GLYAT may be involved in the EMT process. In vitro assays were performed in MDA-MB-231 with GLYAT KD and in MCF with GLYAT OE. Both vimentin and E-cadherin were analysed for their nucleic location via immunofluorescence assays. In GLYAT KD cells, we found an increase in vimentin but a decrease in E-cadherin expression (Fig. 3a-b). Western blotting of EMT-related proteins was consistent with the immunofluorescence results, showing a reduction in E-cadherin, but increase in levels of vimentin, N-cadherin, and fibronectin in MDA-MB-231 GLYAT KD cells (Fig. 3c-g). On the contrary, MCF-7 GLYAT OE cells demonstrated raised E-cadherin protein expression and lowered vimentin levels (Fig. 3h-j). These results suggest that EMT formation could be suppressed by GLYAT.

GLYAT suppresses the EMT of BC cells via the PI3K/Akt/Snail pathway

EMT process is highly regulated by the PI3k/Akt/Snail pathway. We therefore hypothesized that GLYAT might suppress breast cancer metastasis by modulation of this pathway. Western blot assays indicated that p-AKT, AKT, PI3K p58, and SNAI1 were significantly activated in cells with GLYAT KD (Fig. 4a-e). OE cells on the other hand, demonstrated markedly suppressed levels of p-AKT, AKT, PI3K p58, and SNAI1 in contrast to control cells (Fig. 4f-j). Our results indicate the involvement of the PI3K/Akt/Snail signaling pathway in GLYAT-mediated EMT suppression in BC cells.

GLYAT suppresses in vivoBC EMT, proliferation and metastasis

In vivo experiments were performed on nude mice that were received subcutaneous injection with either stable MDA-MB-231 GLYAT KD or MCF-7 GLYAT OE cells (Fig. 5a). The mice injected with stable GLYAT KD cells developed bigger and heavier tumors (Fig. 5b-d). In contrast, the mice injected with stable GLYAT OE MCF-7 cells had markedly smaller and lighter tumors (Fig. 5e-g). Then immunofluorescence assays for E-cadherin, vimentin, and p-AKT were performed using serial sections of mouse tumor tissues. These assays revealed that GLYAT silencing in MDA-MB-231 cells significantly decreased E-cadherin levels, but increased vimentin and p-AKT levels (Fig. 5h). In contrast, in GLYAT OE MCF-7 cells, E-cadherin was increased, whereas vimentin and p-AKT were decreased (Fig. 5i).

Lower GLYAT expression is correlated with poorer prognosis and malignant clinicopathologicalfeatures in human breast cancer tissues

Immunohistochemical experiments found that breast cancer tissues possessed markedly decreased GLYAT expression in contrast to healthy breast tissue (Fig. 6a-b). Additionally, GLYAT expression in BC tissues correlated positively with the degree of pathological differentiation. (Fig.6c). Based on these results, we conclude GLYAT may act as an anti-oncogene in BC.

We then studied the correlation between clinicopathological features and GLYAT expression in 310 patients who were recruited for this investigation. Table 1 demonstrates their characteristics. All individuals were women between 29 and 78 years old. No one had a previous history of malignancy and were chemoradiotherapy-naïve. The patients were grouped based on their GLYAT levels-low (n=143) or high (n=167) expression groups (tumor tissue/adjacent, normal samples). We found GLYAT expression correlated with TNM stage, histological grade, and Ki-67 status (Table 1; all P<0.05). Additionally, we did not uncover any relationship between GLYAT level and menopausal status, age, receptor status (HER, ER, PR), tumor size, node status, pathologic or molecular subtype (Table 1; all P>0.05). These findings suggest GLYAT is associated with malignant clinicopathological characters of BC.

Table 1 Correlation between GLYAT expression level and clinicopathological features in patients with breast cancer

Clinical characteristics	No. (n=310)	GLYAT		P value
Clinical characteristics	No. (n=310)	Low (143)	High (167)
Age (years)
<50	139	59	80	0.241
≥50	171	84	87
Menopausal status
Premenopausal	165	73	92	0.477
Postmenopausal	145	70	75
Tumor Size (cm)
≤2	193	93	100	0.351
>2	117	50	67
Nodes status
Positive	102	41	61	0.142
Negative	208	102	106
TNM stage
Ⅰ	128	40	68	0.033
Ⅱ	133	64	69
Ⅲ	49	39	30
Pathologic type
Invasive ductal carcinoma	202	97	105	0.361
Others	108	46	62
Histological grade
Ⅰ	51	16	35	0.040
Ⅱ	201	95	106
Ⅲ	58	32	26
ER status
Positive	223	102	121	0.826
Negative	87	41	46
PR status
Positive	217	103	114	0.471
Negative	93	40	53
Her-2 status
Positive	45	19	26	0.570
Negative	265	124	141
Ki-67 status
High	91	50	41	0.045
Low	219	93	126
Molecular subtype
Luminal A	156	71	85	0.971
Luminal B	71	32	39
HER2 Positive	30	14	16
Thriple Negative	53	26	27
GLYAT, Glycine N-acyltransferase; TNM, tumor-node-metastasis.

BC patients were analyzed for their GLYAT status and prognosis via Kaplan-Meier survival analysis as well as log-rank tests after being followed up for an average of 61.75 months (range, 9–77 months). We found that those with decreased GLYAT levels experienced poorer disease free survival (DFS) (Fig. 6d; P=0.012). GLYAT expression was a significant indicator of DFS rates in breast cancer patients based on univariate Cox regression analyse, demonstrating a hazard ratio (HR) of 1.565 [P<0.05; 95% confidence interval (CI) 1.098-2.232] (Table 2). Further multivariate analysis revealed that GLYAT expression was significantly related to DFS (HR, 1.475; 95% CI, 0.991-2.197; P<0.05) (Table 2) and represented an independent risk factor of prognosis for BC.

Table 2 Univariate and multivariate analyses of clinicopathological risk factors for disease-free survival among breast cancer patients

	DFS
Variables	Univariate analysis		Multivariate analysis
Variables	HR (95% CI)	P value	HR (95% CI)	P value
Age	1.109 (0.778-1.580)	0.567	1.779 (0.995-3.182)	0.052
Menopausal status	0.835 (0.586-1.191)	0.321	0.562 (0.313-1.009)	0.053
T stage	1.139 (0.784-1.653)	0.044	0.457 (0.245-0.851)	0.014
N stage	0.951 (0.655-1.381)	0.023	0.466 (0.235-0.927)	0.029
TNM stage	0.630 (0.375-1.058)	0.043	0.432 (0.225-0.741)	0.004
Pathological type	0.914 (0.633-1.319)	0.630	0.937 (0.635-1.383)	0.744
Histological grade	1.035 (0.580-1.846)	0.908	0.878 (0.403-1.914)	0.743
ER status	0.726 (0.480-1.097)	0.128	0.510 (0.215-1.210)	0.127
PR status	0.737 (0.496-1.095)	0.131	0.788 (0.422-1.469)	0.453
HER-2 status	1.181 (0.618-2.257)	0.615	1.986 (0.689-5.726)	0.204
Ki-67 status	0.864 (0.574-1.300)	0.483	1.408 (0.621-3.192)	0.412
Molecular subtype	1.677 (0.925-3.041)	0.089	2.910 (0.696-4.177)	0.143
GLYAT status	1.565 (1.098-2.232)	0.013	1.475 (0.991-2.197)	0.046

Patients were then divided into five subgroups based on their molecular subtype, Ki67, ER, PR and HER2 status in order to investigate the effect of GLYAT on prognosis.

In the subgroup analyses by molecular subtype, patients were grouped based on triple negative or HER2, Luminal A, or Luminal B positivity. We found that GLYAT had a negative impact on DFS in all groups, especially in the Luminal A subgroup, where there was a statistically significant difference (P=0.014, Fig. 6e-h).

In the subgroup analyses based on Ki67, those with low GLYAT levels were found to have shorter DFS in comparison to those with high GLYAT, regardless of whether they had high or low Ki67. This difference was statistically significant in those with low Ki67 (P=0.014, Fig. 6i). Similarly, those with low GLYAT levels also experienced shorter DFS in contrast to those with high GLYAT, irrespective of ER+ or ER- status. There was a significantly statistical difference in those of ER+ status (P=0.011, Fig. 6j). The same pattern of findings also applied to subgroup analyses based on PR. We also observed a statistically significant difference in the subgroup of PR+ patients (P=0.011, Fig. 6k). Likewise, HER2 status did not change the impact of low GLYAT on shorter DFS, with statistically significant changes observed in both HER+ or HER- status groups (P=0.046 and P=0.045 , Fig. 6l).

Advanced breast cancer is a condition that carries high morbidity and mortality, despite the existence of several diagnostic markers and treatment strategies [24]. Therefore, we sought to determine novel molecular components involved in the biology of BC that may be useful in the management of this condition.

A potential candidate is the GLYAT protein, which is involved in glycine conjugation of xenobiotics such as benzoic acid, and plays a role in many anabolic and catabolic reactions [25]. In cancerous cells, GLYAT functions as a key metabolite that inhibits glycine uptake or biosynthesis, impairing growth of these cells likely through inhibition of nucleic acid synthesis [26]. Previous reports have highlighted the role of this molecule in musculoskeletal growth and development as well as in hepatocellular carcinoma progression [19]. Nevertheless, there has been no study on its oncological effect in breast cancer.

Our series of experiments found that human breast cancer cells and tissues contained remarkably suppressed levels of GLYAT. We further uncovered that lower GLYAT level was associated with more malignant clinicopathological characteristics, including higher TNM stage, histological grade, and Ki-67 status. DFS of those with lower GLYAT levels were also shorter. Our findings are consistent with information available in public databases. As breast cancer is related to hormone which has different molecular types, we then investigated the relation between GLYAT and patient prognosis based on their molecular subtypes, namely their HER, PR, ER, and Ki67 statuses. Lower GLYAT levels resulted in poorer DFS irrespective of molecular profile, as demonstrated by statistically significant differences in those of the luminal A subgroup, the Ki67 low status subgroup, the ER + subgroup, the PR + subgroup, and both the HER + and HER- subgroups (all P < 0.05). We postulate that the lack of statistical significance for other molecular subtype subgroups, such as the Ki67 high status subgroup, the ER- subgroup, and the PR- subgroup, can be attributed to small sample sizes or an inadequate follow-up time. However, perhaps more importantly, ER and PR status are intricately linked to breast cancer initiation, development and prognosis. GLYAT may act in combination with ER or PR status to impact breast cancer prognosis. Additional experiments are needed to clarify this question in further research.

We further discovered GLYAT to be an independent prognostic factor for BC. These findings indicated that GLYAT acts as antioncogene and is associated with malignant clinicopathological features that may enhance breast cancer metastasis and progression. Our findings based on DFS alone are suggestive of the potential of this molecule to prognosticate breast cancer. Our findings are consistent with bioinformatics analyses and the data reported by Liu et al., who indicated that GLYAT was suppressed in human hepatocellular carcinomas [19]. Given that our study is the first to analyze the correlation between GLYAT and breast cancer, future investigations and follow up investigations regarding OS information (not available in our study due to short follow-up time) are still required.

We further demonstrated that GLYAT regulates breast cancer migration and invasion via EMT modulation. EMT is reduced via alteration of the PI3K/ATK/Snail signaling pathway both in vitro and in vivo. The EMT process participate in cancer metastasis [27]. Numerous studies reported that EMT is regulated by several signaling pathways, for example Notch-, Wnt-, PI3K/Akt- and NF-κB-dependent pathways [28–30]. Previous studies have confirmed that the EMT enhances tumor cell mobility and invasiveness, thus, contributing to the development of chemotherapy resistance and metastasis [31, 32]. Additionally, breast cancer recurrence and prognosis have been reported to be affected strongly by EMT [33]. We believe there may be other mechanisms in addition to the PI3K/AKT/Snail pathway by which GLYAT promotes metastasis that need to be further explored. One example by Ren et al. using a Drosophila model demonstrated GLYAT to be a critical modulator of the c-Jun N-terminal kinase (JNK) signaling pathway [34]. GLYAT downregulation suppresses JNK-dependent ROS activation, which is a key facilitator of several important biological functions.

Besides the impact on metastasis, the proliferation ability of BC cells was increased after inhibiting GLYAT both in vitro and in vivo, further highlighting the tumor suppressed role of GLYAT. Our research is the only study to illustrate the role of GLYAT in BC patients, revealing that down-regulated GLYAT induces EMT and tumor metastasis via PI3K/ATK/Snail signaling.

Collectively, our data revealed that GLYAT downregulated in human BC cells and tissues, and lower GLYAT expression was related to poor clinical outcomes. GLYAT suppresses BC cell proliferation and migration through EMT induction via the PI3K/ATK/Snail pathway in vitro and in vivo. Our study puts forth GLYAT as a novel biomarker in BC and may also represent a therapeutic target for BC treatment.

BC:Breast cancer; NC: Normal control; ER: Estrogen receptor; PR: Progesterone receptor;

HER2: Human epidermal growth factor receptor 2; TNBC: Triple negative breast cancer;

BRCA: Breast invasive carcinoma; IDBC: Invasive ductal breast carcinoma; ILBC: Invasive

lobular breast carcinoma; LIHC: Liver hepatocellular carcinoma; CHOL: Cholangio carcinoma

ACC: Adrenocortical carcinomaACC Adrenocortical carcinoma; BLCA Bladder Urothelial CarcinomaBLCA: Bladder Urothelial Carcinoma; CESC: Cervical

squamous cell carcinoma and endocervical adenocarcinoma; COAD: Colon adenocarcinoma;

DLBC: Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA: Esophageal carcinoma;

GBM: Glioblastoma multiforme; HNSC: Head and Neck squamous cell carcinoma; KIRC:

Kidney renal clear cell carcinoma; LAML: Acute Myeloid Leukemia; LGG: Brain Lower Grade

Glioma; LUAD: Lung adenocarcinoma; LUSC: Lung squamous cell carcinoma; OV: Ovarian

serous cystadenocarcinoma; PAAD: Pancreatic adenocarcinoma; PRAD: Prostate

adenocarcinoma; READ: Rectum adenocarcinoma; STAD: Stomach adenocarcinoma; THCA:

Thyroid carcinoma; THYM: Thymoma;

Acknowledgements

Not applicable.

Authors’ contributions

LCG and GS conceived and designed the experiments as well as contributed to the writing of the manuscript. TX, WLN, JM, ZZY, MJN and WR performed the experiments. TX performed clinical analysis and helped with drawing the figures. LCG and GS revised the paper. All authors read and approved the final manuscript.

Funding

The study was funded by the PhD research startup foundation of Liaoning Province (No. 201805400242), and 345 Talent Project, Shengjing Hospital of China Medical University (No. M0397).

Availability of data and materials

Please contact author for data request.

Ethics approval and consent to participate

All participants signed informed written consent and this study has been passed by the ethics committee of Shengjing Hospital of China Medical University (NO.2018PS234K).

Consent for publication

All authors consent this manuscript to be published.

Competing interests

The authors declare that they have no competing interests.

Author details

¹Department of Oncology, Shengjing Hospital of China Medical University, Shenyang, 110004, P.R. China. ²Key Laboratory of Shengjing Hospital, China Medical University, Shenyang, 110004, P.R. China.

[1] Yang H, Pawitan Y, He W, Eriksson L, Holowko N, Hall P, et al. Disease trajectories and mortality among women diagnosed with breast cancer. BREAST CANCER RES. 2019;21:95.

[2] Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394-424.

[3] Sancho-Garnier H, Colonna M. [Breast cancer epidemiology]. PRESSE MED. 2019;48:1076-84.

[4] Scimeca M, Urbano N, Bonfiglio R, Duggento A, Toschi N, Schillaci O, et al. Novel insights into breast cancer progression and metastasis: A multidisciplinary opportunity to transition from biology to clinical oncology. Biochim Biophys Acta Rev Cancer. 2019;1872:138-48.

[5] Claessens A, Ibragimova K, Geurts S, Bos M, Erdkamp F, Tjan-Heijnen V. The role of chemotherapy in treatment of advanced breast cancer: an overview for clinical practice. Crit Rev Oncol Hematol. 2020;153:102988.

[6] Boulding T, McCuaig RD, Tan A, Hardy K, Wu F, Dunn J, et al. LSD1 activation promotes inducible EMT programs and modulates the tumour microenvironment in breast cancer. Sci Rep. 2018;8:73.

[7] Lourenco AR, Ban Y, Crowley MJ, Lee SB, Ramchandani D, Du W, et al. Differential Contributions of Pre- and Post-EMT Tumor Cells in Breast Cancer Metastasis. CANCER RES. 2020;80:163-9.

[8] Krstic M, Kolendowski B, Cecchini MJ, Postenka CO, Hassan HM, Andrews J, et al. TBX3 promotes progression of pre-invasive breast cancer cells by inducing EMT and directly up-regulating SLUG. J PATHOL. 2019;248:191-203.

[9] Bill R, Christofori G. The relevance of EMT in breast cancer metastasis: Correlation or causality? FEBS LETT. 2015;589:1577-87.

[10] Kar R, Jha NK, Jha SK, Sharma A, Dholpuria S, Asthana N, et al. A "NOTCH" Deeper into the Epithelial-To-Mesenchymal Transition (EMT) Program in Breast Cancer. Genes (Basel). 2019;10.

[11] Kim BG, Sung JS, Jang Y, Cha YJ, Kang S, Han HH, et al. Compression-induced expression of glycolysis genes in CAFs correlates with EMT and angiogenesis gene expression in breast cancer. Commun Biol. 2019;2:313.

[12] Sun X, Chang X, Wang Y, Xu B, Cao X. Oroxylin A Suppresses the Cell Proliferation, Migration, and EMT via NF-kappaB Signaling Pathway in Human Breast Cancer Cells. BIOMED RES INT. 2019;2019:9241769.

[13] Ren Z, Yang T, Zhang P, Liu K, Liu W, Wang P. SKA2 mediates invasion and metastasis in human breast cancer via EMT. MOL MED REP. 2019;19:515-23.

[14] Jiang H, Zhou Z, Jin S, Xu K, Zhang H, Xu J, et al. PRMT9 promotes hepatocellular carcinoma invasion and metastasis via activating PI3K/Akt/GSK-3beta/Snail signaling. CANCER SCI. 2018;109:1414-27.

[15] van der Sluis R, Badenhorst CP, Erasmus E, van Dyk E, van der Westhuizen FH, van Dijk AA. Conservation of the coding regions of the glycine N-acyltransferase gene further suggests that glycine conjugation is an essential detoxification pathway. GENE. 2015;571:126-34.

[16] Guo YF, Zhang LS, Liu YJ, Hu HG, Li J, Tian Q, et al. Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study. HUM GENET. 2013;132:189-99.

[17] Dempsey DR, Bond JD, Carpenter AM, Rodriguez OS, Merkler DJ. Expression, purification, and characterization of mouse glycine N-acyltransferase in Escherichia coli. Protein Expr Purif. 2014;97:23-8.

[18] van der Sluis R, Badenhorst CP, van der Westhuizen FH, van Dijk AA. Characterisation of the influence of genetic variations on the enzyme activity of a recombinant human glycine N-acyltransferase. GENE. 2013;515:447-53.

[19] Matsuo M, Terai K, Kameda N, Matsumoto A, Kurokawa Y, Funase Y, et al. Designation of enzyme activity of glycine-N-acyltransferase family genes and depression of glycine-N-acyltransferase in human hepatocellular carcinoma. Biochem Biophys Res Commun. 2012;420:901-6.

[20] Chandrashekar DS, Bashel B, Balasubramanya S, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. NEOPLASIA. 2017;19:649-58.

[21] Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. NUCLEIC ACIDS RES. 2017;45:W98-102.

[22] Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait BH, et al. A subcellular map of the human proteome. SCIENCE. 2017;356.

[23] Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. NATURE. 2012;486:346-52.

[24] Gu X, Xue J, Ai L, Sun L, Zhu X, Wang Y, et al. SND1 expression in breast cancer tumors is associated with poor prognosis. Ann N Y Acad Sci. 2018;1433:53-60.

[25] Badenhorst CP, van der Sluis R, Erasmus E, van Dijk AA. Glycine conjugation: importance in metabolism, the role of glycine N-acyltransferase, and factors that influence interindividual variation. Expert Opin Drug Metab Toxicol. 2013;9:1139-53.

[26] Lino CC, Bourgine J, Cauffiez C, Allorge D, Lo-Guidice JM, Broly F, et al. Genetic polymorphisms of glycine N-acyltransferase (GLYAT) in a French Caucasian population. XENOBIOTICA. 2010;40:853-61.

[27] Matysiak M, Kapka-Skrzypczak L, Jodlowska-Jedrych B, Kruszewski M. EMT promoting transcription factors as prognostic markers in human breast cancer. ARCH GYNECOL OBSTET. 2017;295:817-25.

[28] De Francesco EM, Maggiolini M, Musti AM. Crosstalk between Notch, HIF-1alpha and GPER in Breast Cancer EMT. INT J MOL SCI. 2018;19.

[29] Yu J, Luo Y, Wen Q. Nalbuphine suppresses breast cancer stem-like properties and epithelial-mesenchymal transition via the AKT-NFkappaB signaling pathway. J Exp Clin Cancer Res. 2019;38:197.

[30] Baek SH, Ko JH, Lee JH, Kim C, Lee H, Nam D, et al. Ginkgolic Acid Inhibits Invasion and Migration and TGF-beta-Induced EMT of Lung Cancer Cells Through PI3K/Akt/mTOR Inactivation. J CELL PHYSIOL. 2017;232:346-54.

[31] Chaffer CL, San JB, Lim E, Weinberg RA. EMT, cell plasticity and metastasis. Cancer Metastasis Rev. 2016;35:645-54.

[32] Huang J, Li H, Ren G. Epithelial-mesenchymal transition and drug resistance in breast cancer (Review). INT J ONCOL. 2015;47:840-8.

[33] Ren Z, Yang T, Zhang P, Liu K, Liu W, Wang P. SKA2 mediates invasion and metastasis in human breast cancer via EMT. MOL MED REP. 2019;19:515-23.

[34] Ren P, Li W, Xue L. GLYAT regulates JNK-mediated cell death in Drosophila. Sci Rep. 2017;7:5183.

Additionalfile1TableS1.docx
The sequence of pGPU6/mCherry/Puro-shRNA-GLYAT and pOGP-T2A-CKNeo-GLYAT-over-expression.
Additionalfile1TableS1.docx
The sequence of pGPU6/mCherry/Puro-shRNA-GLYAT and pOGP-T2A-CKNeo-GLYAT-over-expression.
Additionalfile2FigureS1.pdf
Box plot analysis showed the expression level of GLYAT in different tumor tissues. (A) The GLYAT mRNA levels was decreased in BRCA, LIHC, CHOL. (B) There was no difference of the GLYAT mRNA levels in ACC, BLCA, CESC, COAD, DLBC and normal tissues. (C) The GLYAT mRNA levels was decreased in KIRC, but there was no difference in ESCA, GBM, HNSC, LAML and normal tissues. (D) The GLYAT mRNA levels was decreased in LGG and PAAD, but there was no difference of the GLYAT mRNA levels in LUAD, LUSC, OV and normal tissues. (E) There was no difference of the GLYAT mRNA levels in PRAD, READ, STAD, THCA, THYM and normal tissues.*p<0.05, **p<0.01.
Additionalfile2FigureS1.pdf
Box plot analysis showed the expression level of GLYAT in different tumor tissues. (A) The GLYAT mRNA levels was decreased in BRCA, LIHC, CHOL. (B) There was no difference of the GLYAT mRNA levels in ACC, BLCA, CESC, COAD, DLBC and normal tissues. (C) The GLYAT mRNA levels was decreased in KIRC, but there was no difference in ESCA, GBM, HNSC, LAML and normal tissues. (D) The GLYAT mRNA levels was decreased in LGG and PAAD, but there was no difference of the GLYAT mRNA levels in LUAD, LUSC, OV and normal tissues. (E) There was no difference of the GLYAT mRNA levels in PRAD, READ, STAD, THCA, THYM and normal tissues.*p<0.05, **p<0.01.

Download PDF

Version 1

posted

You are reading this latest preprint version

Downregulation of GLYAT facilities tumor growth and metastasis and poor clinical outcomes through the PI3K/AKT/Snail pathway in human breast cancer

Status:

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

Background

Materials And Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Supplementary Files

Status:

Version 1