Skip to main content

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 16 February 2022
Sec. Signaling
Volume 9 - 2021 | https://doi.org/10.3389/fcell.2021.772380

The Dual Regulation Effects of ESR1/NEDD4L on SLC7A11 in Breast Cancer Under Ionizing Radiation

  • 1School of Public Health and Management, Wenzhou Medical University, Wenzhou, China
  • 2NHC Key Laboratory of Radiobiology, School of Public Health of Jilin University, Changchun, China
  • 3The Second Affiliated Hospital of Jilin University, Changchun, China
  • 4Key Laboratory of Watershed Science and Health of Zhejiang Province, Wenzhou Medical University, Wenzhou, China

Radiotherapy is one of the most important treatments for breast cancer. Ferroptosis is a recently recognized form of regulated cell death that is characterized by lipid peroxidation. However, whether ionizing radiation (IR) could induce ferroptosis in breast cancer and how it works remain unknown. Bioinformatics analysis were performed to screen ferroptosis-related genes differentially expressed in breast tumor tissue and normal tissue. Then, breast cancer cell lines with different estrogen receptor (ER) phenotypes were used for studies in vitro, including ER-positive (MCF-7 and ZR-75-1) and ER-negative (MDA-MB-231) cells. The dynamic changes of mRNA and protein levels were examined after x-ray of 8 Gy by qRT-PCR and Western blotting, respectively. Immunoprecipitation (IP) was used to explore the interaction between proteins. Luciferase assay was used to analyze the transcriptional regulation effect of ESR1 on SLC7A11. BODIPY C11 and trypan blue dyes were used to determine lipid peroxidation and cell death, respectively. The result showed that the ferroptosis-related gene SLC7A11 was higher in breast cancer tissues compared with normal tissues and associated with poor survival. A positive correlation exists between ESR1 and SLC7A11 expression. ESR1 promoted SLC7A11 expression at the early stage after IR. ESR1/SLC7A11 knockdown significantly enhanced IR-induced ferroptosis in ER-positive cells. At 12 h after IR, the IP data showed the interaction between E3 ubiquitin ligase NEDD4L and SLC7A11 increased, followed by the ubiquitylation and degradation of SLC7A11. Thus, SLC7A11 expression was regulated by both ESR1 and NEDD4L, in opposite ways. For the first time, we elucidated that ESR1 and NEDD4L functioned together after radiation treatment and finally induced ferroptosis in breast cancer cells, which provides novel insight into the guidance of clinical treatment of breast cancer.

Introduction

Breast cancer is a complex and clinically heterogeneous disease. Similar clinicopathologic characteristics and histologic types may have different clinical course and survival rates, and distinct molecular subtypes of tumors are inextricably bound up with clinical outcomes (Kumar Atri et al., 2017; Lian et al., 2017). According to the expression status of progesterone receptor (PR), estrogen receptor (ER), and human epidermal growth factor receptor 2 (HER2), immunohistochemistry is often used to characterize four subtypes: luminal A (ER+ and HER2), luminal B (HR+ and HER2+), triple-negative/basal-like (ER, PR, and HER2), and HER2-enriched (ER, HER2+, and PR) (Meyers et al., 2011). ERα and Erβ are two main forms of ER, encoded by separate genes (ESR1 and ESR2), respectively. Clinically, ESR1 is used to define the ER status of breast cancer (Lin et al., 2016). It has been reported that 70% of breast cancer are ER-positive at diagnosis (Khan et al., 2019). Under these circumstances, high expression of ER functions in tumor development and progression (Stoner et al., 2002; Ferreira Almeida et al., 2020). ER is a ligand-activated transcription factor involved in the modulation of cell division, which is vital to cell growth and survival (Cho et al., 2005; Carroll, 2016).

Ferroptosis is a novel type of programmed cell death distinct from apoptosis, necrosis, and autophagy, which is an iron-dependent form of regulated cell death involving accumulation of reactive oxygen species (ROS) and overwhelming lipid peroxidation (Stockwell et al., 2017). SLC7A11 is one of the critical regulators of ferroptosis, which promotes cystine uptake and glutathione biosynthetic and is used to generate the anti-oxidant glutathione (GSH), thus converting lipid hydroperoxides to lipid alcohol to protect cells from death (Liu et al., 2021). Furthermore, IR induces ferroptosis in cancer patients and improves the therapeutic effects and survival of radiation therapy (Wu et al., 2020a; Lei et al., 2020). Recently, studies have revealed that IR induces ferroptosis by the suppression of SLC7A11, which represents an important part of radiotherapy (Lang et al., 2019). It is also shown that IR induces the expression of SLC7A11, as an adaptive response (Lei et al., 2020). Many studies have indicated that ferroptosis also mediated cell death in breast cancer by SLC7A11 (Xu et al., 2020; Yang et al., 2021), which implied the important role of SLC7A11 in regulating the ferroptosis of breast cancer. However, the mechanisms for IR-induced ferroptosis and regulation in breast cancer cells remain largely unknown.

Ionizing radiation (IR) can directly interact with water in human bodies, and generate ROS, which could be amplified by Fenton reaction and generate an extremely reactive hydroxy radical (Ahsan et al., 2003; Yu, 2012). The hydroxy radical results in lipid peroxide and induces cell death, including ferroptosis (Kajarabille and Latunde-Dada, 2019). Radiotherapy is an essential mode for many cancers, including breast cancer, particularly those seeking breast-conserving surgery. Radiotherapy to breast conservation surgery can halve the rate of disease recurrence and reduce breast cancer mortality by about one-sixth (Early Breast Cancer Trialists’ Collaborative Group (EBCTCG) et al., 2011). According to NCCN clinical practice guidelines in oncology, the adjuvant therapy and overall treatment pathway should be adjusted according to the level of ER expression, although in the ER-positive group (Gradishar et al., 2020). In this study, we designed a systematic study to elucidate the underlying mechanisms of radio-resistance in ER-positive breast cells and managed to improve radiation killing effect through targeted intervention of IR-induced ferroptosis.

Materials and Methods

Data Collection

The gene expression profiles ( raw counts and FPKM-UQ) were downloaded from the TCGA database (https://gdc-portal.nci.nih.gov) by the TCGA biolinks R package (version 2.20.0) (Colaprico et al., 2016). Data from metastatic, male, or formalin-fixed paraffin-embedded tissue samples, as well as repeated data for the same samples, were excluded. After exclusion, 1,191 primary solid tumor (TP) and 112 solid tissue normal (NT) samples were retained for genomic data. The survival and clinical information of breast cancer were obtained from the Xena Functional Genomics Explorer of University of California Santa Cruz (https://xenabrowser.net/).

Reagents, Antibodies, and Plasmids

Monoclonal anti-FLAG antibody and anti-β-actin were obtained from Sigma (St. Louis, United States). Rabbit polyclonal anti-SLC7A11, and rabbit polyclonal anti-ESR1 were from Cell Signaling Technology (CST, United States). Protein G Plus-agarose was from Sangon Biotech (Shanghai, China). RNAiso Plus kits, reverse transcription kits, and real-time PCR kits were from Takara Biomedical Technology (Dalian, China). Ferrostatin-1 (Fer-1), Erastin, and carbobenzoxy-Leu-Leu-leucinal (MG132) were all purchased from MedChemExpress (Shanghai, China). The expression plasmids for MYC-Ub were from Dr. Haishan Huang’s lab (Wenzhou Medical University, China). SLC7A11 and ESR1 genes were obtained by PCR using MCF-7 cDNA as a template and then inserted into pcDNA3.1-FLAG vectors.

Cell Culture and siRNA, Plasmid Transfection

The human-source cell lines MCF-7, ZR-75-1, and MDA-MB-231 were obtained from the Cell Bank of Chinese Academy of Sciences and periodically tested for the absence of mycoplasma. The cell lines were grown in high-glucose DMEM medium (Gibco, United States) supplemented with 10% fetal bovine serum and 1% penicillin and streptomycin in a 37 C incubator maintained at 5% CO2. For transfection, when the cells had grown to 70%–80% confluence in six-well plates, the plasmid (1 μg) or siRNA (100 μM) was transfected using Lipofectamine 2000 (Invitrogen, United States) according to the manufacturer's instructions. The siRNAs were synthesized by GenePharma (Shanghai, China), and the sequences were listed as follows: siRNA for ESR1, CTA​CAG​GCC​AAA​TTC​AGA​TAA; NEDD4L, CGC​CTT​GAC​TTA​CCT​CCA​TAT.

Flow Cytometric Analysis

Cell death and lipid peroxidation were assessed by flow cytometry with trypan blue (Gibco, United States) and BODIPY™ 581/591 C11 (Thermo Fisher Scientific, United States) following the manufacturer’s instruction, respectively. In brief, cells were pretreated with 10 µM of Erastin or 10 µM of Fer-1 for 1 h and then exposed to IR. Cells were collected at 48 h for lipid peroxidation or 72 h for cell death testing after X-ray exposure. Cells were resuspended in 400 μl of serum-free medium that contained BODIPY™ 581/591 C11 (5 µM) at 37°C in the dark for 30 min, or in 200 μl of PBS containing 0.04% of trypan blue, respectively to detect lipid peroxidation and cell death. The data were analyzed using the NovoExpress software (ACEA, United States), and a minimum of 10,000 cells were analyzed per condition.

Western Blot and Immunoprecipitation

For Western blot, cells were lysed in ice-cold RIPA Lysis Buffer. Lysates were centrifuged at 12,000 rpm for 15 min at 4°C. The supernatant was collected, and the protein content was determined using the BCA Protein Assay Kit (Thermo-Fisher, United States) according to the manufacturer’s instructions. Equal amounts of total protein (20 μg) were separated by 10% or 12% SDS–PAGE and then transferred to an equilibrated PVDF membrane (Millipore, United States). β-actin was used as a loading control. For IP, MCF-7 cells were transfected with expression vectors for FLAG-SLC7A11, and were exposed to IR at 48 h post-transfection. At 24 h after IR, cells were lysed in 1% NP-40 lysis buffer. Whole cell lysates (WCL) were incubated with the specified antibodies overnight at 4°C with constant rotation. Then, protein G Plus-agarose was added for 2 h under rotation at 4°C. After incubation, immunocomplexes were washed five times with IP wash buffer (Sangon Biotech, China) and subsequently resolved by SDS-PAGE followed by Western blot analysis.

Ubiquitination Assay

Cells were transiently transfected with MYC-tagged ubiquitin (UB), and FLAG-tagged SLC7A11 with or without si-NEDD4L for 48 h, and were exposed to IR for 24 h. At 24 h after IR, cells were pretreated with MG132 (10 μM) for 4 h before harvest. The lysates were subjected to immunoprecipitation with FLAG antibodies overnight at 4°C. Ubiquitinated proteins were determined by immunoblotting with anti-MYC antibodies (Santa Cruz, United States).

Real-Time PCR for mRNA Quantification

Total RNA was isolated from cultured cells using RNAiso Plus kits according to the manufacturer’s instructions. mRNAs were first reverse transcribed into cDNA and quantitative reverse transcription-PCR (qPCR) was performed as the following thermal conditions: 95°C for 30 s followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. qPCR was performed on QuantStudio 3 (Applied Biosystems, United States), and the expression of target genes was normalized to that of GAPDH. The following primers were used: GAPDH forward CCA​TGG​GTG​GAA​TCA​TAT​TGG​A and reverse TCA​ACG​GAT​TTG​GTC​GTA​TTG​G; SLC7A11 forward GCG​TGG​GCA​TGT​CTC​TGA​C and reverse GCT​GGT​AAT​GGA​CCA​AAG​ACT​TC; and ESR1 forward GCT​TAC​TGA​CCA​ACC​TGG​CAG​A and reverse GGA​TCT​CTA​GCC​AGG​CAC​ATT​C.

Double Fluorescein–Reporter Enzyme

Double luciferase–reporter gene determination was conducted to determine whether SLC7A11 the downstream target gene of ESR1. Reporter plasmid of SLC7A11 promoter double luciferase was synthesized from Generalbiol (Anhui, China) and inserted into the pGL3-Basic vector. pGL3-Basic SLC7A11 promoter was co-transfected into cells with ESR1 overexpression vector or control vector by Lipofectamine 2000. After transfection for 48 h, luciferase activity was determined on a Centro LB 960 Luminometer (Berthold Technologies, Germany) and the activity of renilla luciferase was used as a standardized control.

Statistical Analysis

Statistical analysis was performed using R version 4.1.1 (http://www.r-project.org/) and GraphPrism 5.0 software (GraphPad Software, United States). The differentially expressed gene (DEG) analyses were done with R package edgeR (version 3.34.0) (Robinson et al., 2010). According to the suggested workflow, raw RNA-Seq read count data were normalized by Trimmed-mean M values (TMM) to balance the counts among the different samples (Maza, 2016). Genes were marked as DEGs if the absolute value of log2 transformed fold change (log2FC) ≥ 1.5 and the false discovery rate (FDR) < 10-10. The FPKM-UQ normalized gene expression data were log-transformed for further analyses. The R package survminer (version 0.4.9) was used to separate patients into higher and lower expression groups by gene expression. The relationship between SLC7A11 expression status and clinicopathologic characteristics of breast cancer was analyzed by Chi-square test. The association with prognosis of breast cancer patients was determined by Kaplan–Meier survival analysis with Log-rank test. A Cox proportional hazard regression model was conducted to perform a univariate and multivariate analysis. All experiments were performed in triplicate, and replicate data were used to calculate means ± standard deviation (SD). ANOVA plus subsequent Dunnett post hoc tests were performed for multiple comparisons. If not specified above, a p-value less than 0.05 was considered to be statistically significant, and p-values were two-tailed.

Result

Differential Expression of Ferroptosis-Related Genes in Breast Tumor Samples

Breast cancer data were collected from TCGA. Based on the screening principles, 1191 TP and 112 NT samples were selected into this study. DEGs were reflected in volcano plot (Figure 1A). According to the principle of statistics, 4 upregulated and 6 downregulated ferroptosis-related DEGs were associated with breast cancer. Among them, SLC7A11 was a critical regulator of ferroptosis (Koppula et al., 2020). The SLC7A11 mRNA expression was extremely upregulated in tumor tissues compared with normal tissues (Figure 1B). Excluding unavailable SLC7A11 expression value, 1,104 breast cancer samples including 217 high expressions (19.66%) and low 887 expressions (80.34%) were extracted for further study. Chi-Square analysis were applied to study the data and results revealed that there were many clinicopathologic features significantly correlated with SLC7A11 expression, including age, histological type, molecular subtype, ER, PR, and N classification (Supplementary Table S1). Differences in SLC7A11 expression are shown in Figure 1C according to ER subtype. SLC7A11 expression was lower in ER-positive tissues than ER-negative tissues.

FIGURE 1
www.frontiersin.org

FIGURE 1. Compared with normal tissues, ferroptosis-related gene SLC7A11 was upregulated in breast cancer tissues and contributed to poor survival. (A) The differentially expressed genes (DEGs) were identified between TP and NT with in breast cancer patients. Volcano plot shows the gene statistical significance with |log2FC| > 1.5 and FDR <10-10. The green points were ferroptosis-related genes. (B,C) Comparison of SLC7A11 expression in different subgroups. Normal and tumor (B); ER-positive and ER-negative (C). (D–F) Kaplan–Meier curves of overall survival (OS) in breast cancer according to SLC7A11 expression. Overall survival analysis (D); Subgroup analyses of OS in ER-positive (E) and ER-negative (F).

The High Expression of SLC7A11 Predicted Poor Prognosis in Breast Cancer

In univariate Cox regression for overall survival (OS) model, we included variables that correlated with SLC7A11. The following multivariate analyses confirmed SLC7A11 expression, age, and N classification as independent prognostic factors after adjusting other prognostic indicators (Supplementary Table S2). Analysis demonstrated that patients with higher SLC7A11 were associated with worse OS than those with lower SLC7A11 (Figure 1D). In addition, K–M curves were used to evaluate the prognostic value in differentiation ER subtypes, which revealed that SLC7A11 specifically associated with OS in ER-positive cancers (Figure 3E), but not in ER-negative cancers (Figure 1F).

Then, we further define the relationship between ferroptosis and ER-positive cancers, and eight ferroptosis-related genes were analyzed according to Wu et al. (2020b) and Tang et al. (2021). Results revealed that only low expression of SLC7A11 supported ferroptosis occurrence in ER-positive cancers (Supplementary Figure S1). The above data suggest that SLC7A11 may play an important role in the regulation of ER-positive cell ferroptosis.

IR Induced Ferroptosis by the Reduction of SLC7A 11 Expression

The IR-induced ferroptosis in different breast cancer cell lines were evaluated in our previous study (Ma et al., 2021). Fer-1 was able to significantly attenuate IR-induced cell death in ER-positive MCF-7 and ZR-75-1, but not in ER-negative MDA-MB-231 cells. To further verify the result, we analyzed the effects of Erastin and Fer-1 on lipid peroxidation accumulation and cell death in response to IR treatment. Importantly, Fer-1 significantly reduced IR-stimulated lipid peroxidation and cell death. Correspondingly, Erastin significantly increased lipid peroxidation accumulation and cell death, and showed synergistic effect with IR (Figures 2A–D). When we knockdown of SLC7A11 in ZR-75-1, IR-induced lipid peroxidation and cell death significantly increased, which could be rescued by Fer-1 (Figures 2E–J). These results suggested that ferroptosis was one pattern of IR-induced cell death in ER-positive cells, and SLC7A11 may play an important role in this process.

FIGURE 2
www.frontiersin.org

FIGURE 2. IR-induced ferroptosis in ER-positive breast cell lines. (A–D) MCF-7 cells were pretreated with Erastin (10 µM) or Fer-1 (10 µM) for 1 h and then exposed to IR (8 Gy). (E–J) Western blot and flow cytometry analysis of ferroptosis by SLC7A11 downregulation in ZR-75-1. Lipid peroxidation (A,B,G,H) and cell death (C,D,I,J) were assessed by flow cytometry using C11-BODIPY and trypan blue, respectively. Shown is mean ± SD, n ≥ 3, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

ESR1 Activates the Transcription of SLC7A11

Referring to the above results, we considered that the relationship between SLC7A11 expression and breast cancer survival was probably based on ESR1. Furthermore, a positive correlation of ESR1 with SLC7A11 was found in ER-positive cancers (Figure 3A) and the survival curve of ESR1 expression in ER-positive cancers is shown in Figure 3B. From the K–M curve, higher expression of ESR1 was also related to poor prognosis. From these results, we can conclude that SLC7A11 and ESR1 are both related to survival in ER-positive cancers. However, the precise relationship is unclear.

FIGURE 3
www.frontiersin.org

FIGURE 3. Correlation between ESR1 and SLC7A11 in ER-positive cancers. (A) Scatter plots for correlation analysis of ESR1 and SLC7A11. The dark gray line means a linear regression line and the gray region is 95% confidence region. (B) Kaplan–Meier curve of OS in ER-positive cancers according to ESR1 expression. (C) IP assay of MCF-7 cells transfected with a plasmid expressing FLAG-tagged SLC7A11. (D–G) Western blot and qRT-PCR analysis showed that si-ESR1 reduced SLC7A11 expression in MCF-7 (D,F,G) and ZR-75-1 (E). (H) Effect of ESR1 on SLC7A11 transcription by double fluorescein-reporter enzyme activity in MCF-7. Shown is mean ± SD, n ≥ 3, *p < 0.05; **p < 0.01.

Then, we designed siRNA targeting ESR1 to verify the function of ESR1. To determine whether ESR1 could interact with SLC7A11, an IP experiment was carried out. Our result suggested that there was no direct protein–protein interaction (Figure 3C). In order to determine the relationship between ESR1 and SLC7A11, MCF-7 and ZR-75-1 cells were transfected with siRNA-NC or siRNA-ESR1; 72 h later, cells were harvested, and the total mRNA and proteins were extracted and subjected to qPCR and Western blotting, respectively. As shown in Figures 3D–G, si-ESR1 significantly reduced ESR1 mRNA and protein levels, following the SLC7A11 levels decreased. Further dual-luciferase reporter gene assay was applied and the results showed that EGR1 was capable of activating the SLC7A11 promoter ESR1 was able to activate the transcription of SLC7A11 (Figure 3H). These findings support that ESR1 suppresses ferroptosis by promoting SLC7A11 expression in IR treatment.

In addition, the mRNA levels of SLC7A11 and ESR1 were evaluated by qRT-PCR based on time course. After exposure to X-ray, there was an increase in the mRNA levels of SLC7A11 and ESR1 (Figures 4A,B). The protein levels of SLC7A11 and ESR1 based on IR time course were evaluated by western blot (Figure 4C). The results suggested that IR-induced ferroptosis occurred at 48 h, but at 4–12 h SLC7A11 had a rising period to against cell death.

FIGURE 4
www.frontiersin.org

FIGURE 4. IR-induced SLC7A11 and ESR1 expression in ER-positive cell lines. Cells were collected at different time points after IR with 8 Gy x-rays. Western blot, and qRT-PCR were used to detect the expression of SLC7A11 and ESR1. (A,B) Relative mRNA expression of ESR1 (A) and SLC7A11 (B) were detected in MCF-7. Shown is mean ± SD. (C) SLC7A11 and ESR1 protein expression by IR treated in MCF-7 and ZR-75-1.

ESR1 Functions as a Cell Survival Response in Early Effects of IR

Lipid peroxidation and cell death were assayed by flow cytometry using the fluorescent probes C11-BODIPY and trypan blue, respectively. ESR1 knockdown rendered ER-positive cells susceptible for ferroptosis and lipid peroxidation accumulation induced by IR. When ESR1 was knockdown, the increase of SLC7A11 induced by IR was shown to significantly attenuate in MCF7 (Figure 5A) and ZR-75-1 (Figure 5F). The suppression of ESR1 by siRNA increased IR-induced cell death as well as rise of lipid peroxidation in MCF-7 (Figures 5B–E). Moreover, IR had similar effects on cell death and lipid peroxidation with ESR1 knockdown in ZR-75-1 (Figures 5G–J). Indeed, ESR1 knockdown showed a significant effect on ferroptosis induced by IR. The increased lipid peroxidation and cell death can be significantly rescued by Fer-1.

FIGURE 5
www.frontiersin.org

FIGURE 5. ESR1-SLC7A11 axis is involved in IR-induced ferroptosis. The results of Western blot assay show that ESR1 knockdown inhibited the increase of SLC7A11 induced by IR at 12 h in MCF-7 (A) and ZR-75-1(F). (B–E,G–J) Lipid peroxidation and cell death were assessed by flow cytometry using C11-BODIPY and trypan blue, respectively. Shown is mean ± SD, n ≥ 3, *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

On the basis of these findings, a rescue assay was conducted to confirm the role of ESR1 in SLC7A11-mediated IR-induced ferroptosis. As shown in Figures 6A–D, SLC7A11 can also rescue siESR1-increased lipid peroxidation and cell death. Over-expression of SLC7A11 rescue IR-induced ferroptosis in ESR1 knockdown cells.

FIGURE 6
www.frontiersin.org

FIGURE 6. Characterization of rescue cells and ESR1 overexpression. (A–D) MCF-7 cells were co-transfected with Si-ESR1 and pcDNA3.1-SLC7A11 for 48 h, then treated with IR (8 Gy). (E) Western blotting analysis of ESR1 overexpression in MDA-MB-231. (A–D, F–I) Lipid peroxidation and cell death were assessed by flow cytometry using C11-BODIPY and trypan blue, respectively. Shown is mean ± SD, n ≥ 3, **p < 0.01; ***p < 0.001.

For further verification of the ESR1 effect, we overexpressed ESR1 in ER-negative cells. The data showed that the upregulation of ESR1 could increase SLC7A11 (Figure 6E) and rescue IR-induced lipid peroxidation (Figures 6F,G) and cell death (Figures 6H,I). Our results indicate that ESR1 promotes SLC7A11 transcription to resistance ferroptosis and cytotoxic stress, thereby promoting tumor malignancy.

The NEDD4L-SLC7A11 Interaction Increased After Radiation Treatment

We next detected whether the reduction of SLC7A11 was regulated at the transcriptional level or posttranscriptional level. The mRNA and protein levels of SLC7A11 inconsistently changed in IR treated at 24 h (Figure 4B), indicating that IR regulated the expression of SLC7A11 post-transcriptionally. To test this assumption, MG132, which acts as a blocker in the ubiquitin–proteasome pathway, was used and we found that SLC7A11 expression was accumulated by MG132 (Figures 7A,B). These observations suggest that SLC7A11 is regulated by the ubiquitin–proteasome pathway. By using the UbiBrowser database, NEDD4L is the one with the highest confidence to interact with SLC7A11 (Supplementary Table S3) and serves as a potential E3 ligase (http://ubibrowser.ncpsb.org.cn/ubibrowser/) (Li et al., 2017). Our data showed that NEDD4L knockdown rescued IR-induced degradation of SLC7A11 (Figure 7C). We then confirmed the interaction between SLC7A11 and NEDD4L by IP experiments, and a stronger binding was detected after IR treatment (Figure 7D), suggesting that NEDD4L associates with SLC7A11 in ER-positive cells. Consistent with these observations, we found that knockdown of NEDD4L markedly reduced the ubiquitination of SLC7A11 in ER-positive cells (Figure 7E). Taken together, IR could increase the interaction between NEDD4L and SLC7A11.

FIGURE 7
www.frontiersin.org

FIGURE 7. SLC7A11 is degraded by ubiquitination after IR treatment. MCF-7 (A) and ZR-75-1 (B) cells were pretreated with MG132 (10 μM) or DMSO for 4 h before harvest. (C) MCF-7 cells were transfected with Si-NEDD4L. The protein level of SLC7A11 was analyzed by Western blot. (D) Western blot analysis of WCL and IP samples of IgG or anti-NEDD4L antibody were conducted to detect the interaction between SLC7A11 and NEDD4L. MCF-7 cells were transfected with FLAG-SLC7A11 for 48 h, then treated with IR. MG132 (10 μM) was added into the culture medium for 4 h before harvest. (E) Western blotting analysis of WCL and ubiquitin IP samples were performed with indicated antibody in MCF-7 cells. The cells were treated with MG132 (10 μM) for 4 h before harvest.

Discussion

Ferroptosis is a form of programmed cell death driven by cellular metabolism and iron-dependent lipid peroxidation, but its role in breast cancer remains unclear. Our study found that SLC7A11 was correlative to survival in breast cancer and high expression implied poor prognosis. Moreover, in multivariate analysis, high expression of SLC7A11 was an independent prognostic factor for OS of breast cancer. The prognostic analysis of SLC7A11 in subtypes showed that SLC7A11 is specifically associated with OS in ER-positive cancers. It is reported that SLC7A11 plays an important role in the protection of cells from oxidative stress and ferroptosis by promoting the uptake of cystine and the synthesis of GSH (Stockwell et al., 2017; Koppula et al., 2020; Sun et al., 2021). Our research suggested that ferroptosis occupies a very important position in IR-induced cell death in ER-positive cell lines. We also found that ESR1 expression was also related to prognosis in ER-positive cancers, and our finding was consistent with Ogłuszka’s group, which reported that lowered expression of ESR1 was more favorable for disease-free survival (Ogłuszka et al., 2019). ESR1 is a member of the nuclear receptor superfamily; after the activation by estrogen, it can translocate into the nucleus and function as a transcription factor (Nilsson and Gustafsson, 2011; Lai et al., 2015). Among total breast cancer cases, 70% of cases are ER-positive (Khan et al., 2019). In addition, the presence of elevated levels of ESR1 in benign breast epithelium appears to indicate an increased risk of breast cancer, and the downregulation of ESR1 results in induction of autophagy and apoptosis (Ali and Coombes, 2000; Zheng et al., 2015). A recent report showed the association of ESR1 with ferroptosis in breast cancer, especially in cells with low ESR1 expression (Yu et al., 2019). The unclear relationship between ESR1 and SLC7A11 is puzzling and in great need of clarification.

Ferroptosis is characterized by the overwhelming iron-dependent accumulation of lipid hydroperoxides. Many reports provide definitive evidence for IR-induced lipid peroxidation, and ferroptosis is as important as other forms of cell death (Lang et al., 2019; Lei et al., 2020; Ye et al., 2020). IR induces or reduces the expression of ferroptosis inhibitor SLC7A11 to regulate ferroptotic cell death in different cell lines (Lang et al., 2019; Lei et al., 2020). Our present study finds that SLC7A11 correlated to IR-induced ferroptosis, with SLC7A11 transiently increased as an adaptive response and reduced continuously thereafter in ER-positive cell lines. It is important to emphasize that the mechanism of IR-induced ferroptosis is complicated. In this study, we focused on IR-induced ferroptosis regulated by SLC7A11. ESR1 is an estrogen-activated transcription factor and is responsible for transcription of cell-critical genes such as those involved in regulating proliferation (Khorshid Shamshiri et al., 2020). We detected SLC7A11 changes following ESR1 knockdown and found a positive correlation between ESR1 and SLC7A11 expression. Dual-fluorescence reporter assay confirmed that ESR1 upregulated SLC7A11 by promoting SLC7A11 gene transcription. ESR1 knockdown suppressed the IR-induced SLC7A11 transient increase. In addition, ESR1 knockdown also had a significant positive effect on ferroptosis. Based on the above results, ESR1 a central role in upregulating anti-ferroptosis defense SLC7A11 for ER-positive cancer in IR treatment. SLC7A11 was transiently increased after IR and reduced continuously thereafter. IR was able to increase the expression of ESR1, which activated SLC7A11 expression. ESR1 played a significant role in the initial stage against IR-induced ferroptosis. However, IR increased SLC7A11 protein degradation stronger than the protein synthesis in the long term, and then ferroptosis occurred. However, the exact regulatory mechanism underlying the SLC7A11-mediated ferroptosis after IR in breast cancer remains unknown and needs further study.

In this study, we further investigated the mechanism of SLC7A11 degradation induced by IR. NEDD4L is an E3 ubiquitin ligase, containing four WW domains for substrate interaction (Yang and Kumar, 2010). To the best of our knowledge, we report for the first time that NEDD4L binds to SLC7A11 for SLC7A11 degradation in IR treatment. A recent evidence suggests that NEDD4L likely displays maximal activities when localized to membrane (Escobedo et al., 2014). NEDD4L activity is influenced by the conformational changes, with specific interaction between phosphatidylinositol 4,5-bisphosphate (Escobedo et al., 2014). These results provide available evidence for NEDD4L to catalyze ubiquitination of membrane proteins such as SLC7A11. In addition, the mechanism by which NEDD4L regulates SLC7A11 remains to be elucidated, and it is important to determine the motif residues involved in the interaction between SLC7A11 and Ub.

In conclusion, the downregulation of ESR1 leads to a decrease in SLC7A11 expression, resulting in acceleration of IR-induced ferroptosis. Our findings reveal a novel insight that IR-induced ferroptosis was via two-way regulation of ESR1/NEDD4L to SLC7A11, and ESR1 might be a critical factor responsible for conferring the resistance to radiotherapy and tumor malignancy. We propose that targeting ESR1 in combination with radiotherapy could enhance the radiotherapeutic efficacy of ER-positive cancers.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author Contributions

RL contributed to collection of data, statistical analysis, interpretation of results, and manuscript writing. LL, YB, SZ, and YW conducted parts of the experiments and writing of the final manuscript. HC and XJ contributed to data collection and analysis. GL, QC, CX, ML, and LL analyzed the data, and discussed the results and recommendations. XL and SM experimental design, conceived, and supervised the project. All authors provided revision and final approval of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (81773363, 31700738, 81673092, 81972969, and 81872558), Science and Technology Development Planning Project of Jilin Province (20180520127JH), and The National Key Research and Development Program of China (2017YFC0108602).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors kindly thank Haishan Huang and Honglei Jin for their technical help, and Yi Liu for careful reading of the manuscript and discussion.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.772380/full#supplementary-material

Abbreviations

CI, confidence interval; COX, Cox proportional-hazards model; DEGs, differentially expressed genes; DMSO, dimethyl sulfoxide; ER, estrogen receptor; ESR1, estrogen receptor 1; Fer-1, Ferrostatin-1; GSH, glutathione; HER2, human epidermal growth factor receptor 2; HR, hazard ratio; IP, immunoprecipitation; IR, ionizing radiation; MG132, carbobenzoxy-leu-leu-leucinal; NT, solid tissue normal; OS, overall survival; PR, progesterone receptor; ROS, reactive oxygen species; SD, standard deviation; SLC7A11, solute carrier family 7 member 11; TCGA, the cancer genome atlas; TNBC, triple-negative BRCA; TP, primary solid tumor; qRT-PCR, real-time quantitative reverse transcription PCR; UB, ubiquitin; WCL, whole cell lysate.

References

Ahsan, H., Ali, A., and Ali, R. (2003). Oxygen Free Radicals and Systemic Autoimmunity. Clin. Exp. Immunol. 131, 398–404. doi:10.1046/j.1365-2249.2003.02104.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, S., and Coombes, R. C. (2000). Estrogen Receptor Alpha in Human Breast Cancer: Occurrence and Significance. J. Mammary Gland Biol. Neoplasia 5, 271–281. doi:10.1023/a:1009594727358

CrossRef Full Text | Google Scholar

Carroll, J. S. (2016). EJE PRIZE 2016: Mechanisms of Oestrogen Receptor (ER) Gene Regulation in Breast Cancer. Eur. J. Endocrinol. 175, R41–R49. doi:10.1530/EJE-16-0124

CrossRef Full Text | Google Scholar

Cho, J., Kim, D., Lee, S., and Lee, Y. (2005). Cobalt Chloride-Induced Estrogen Receptor α Down-Regulation Involves Hypoxia-Inducible Factor-1α in MCF-7 Human Breast Cancer Cells. Mol. Endocrinol. 19, 1191–1199. doi:10.1210/me.2004-0162

PubMed Abstract | CrossRef Full Text | Google Scholar

Colaprico, A., Silva, T. C., Olsen, C., Garofano, L., Cava, C., Garolini, D., et al. (2016). TCGAbiolinks: an R/Bioconductor Package for Integrative Analysis of TCGA Data. Nucleic Acids Res. 44, e71. doi:10.1093/nar/gkv1507

PubMed Abstract | CrossRef Full Text | Google Scholar

Early Breast Cancer Trialists’ Collaborative Group (EBCTCG) Darby, S., Darby, S., McGale, P., Correa, C., Taylor, C., Arriagada, R., et al. (2011). Effect of Radiotherapy after Breast-Conserving Surgery on 10-year Recurrence and 15-year Breast Cancer Death: Meta-Analysis of Individual Patient Data for 10,801 Women in 17 Randomised Trials. Lancet 378, 1707–1716. doi:10.1016/S0140-6736(11)61629-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Escobedo, A., Gomes, T., Aragón, E., Martín-Malpartida, P., Ruiz, L., and Macias, M. J. (2014). Structural Basis of the Activation and Degradation Mechanisms of the E3 Ubiquitin Ligase Nedd4L. Structure 22, 1446–1457. doi:10.1016/j.str.2014.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira Almeida, C., Oliveira, A., João Ramos, M., Fernandes, P. A., Teixeira, N., and Amaral, C. (2020). Estrogen Receptor-Positive (ER+) Breast Cancer Treatment: Are Multi-Target Compounds the Next Promising Approach? Biochem. Pharmacol. 177, 113989. doi:10.1016/j.bcp.2020.113989

PubMed Abstract | CrossRef Full Text | Google Scholar

Gradishar, W. J., Anderson, B. O., Abraham, J., Aft, R., Agnese, D., Allison, K. H., et al. (2020). Breast Cancer, Version 3.2020, NCCN Clinical Practice Guidelines in Oncology. J. Natl. Compr. Canc. Netw. 18, 452–478. doi:10.6004/jnccn.2020.0016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kajarabille, N., and Latunde-Dada, G. O. (2019). Programmed Cell-Death by Ferroptosis: Antioxidants as Mitigators. Ijms 20, 4968. doi:10.3390/ijms20194968

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, A., Ashfaq-Ur-Rehman, , Junaid, M., Li, C.-D., Saleem, S., Humayun, F., et al. (2019). Dynamics Insights into the Gain of Flexibility by Helix-12 in ESR1 as a Mechanism of Resistance to Drugs in Breast Cancer Cell Lines. Front. Mol. Biosci. 6, 159. doi:10.3389/fmolb.2019.00159

PubMed Abstract | CrossRef Full Text | Google Scholar

Khorshid Shamshiri, A., Afzaljavan, F., Alidoust, M., Taherian, V., Vakili, F., Moezzi, A., et al. (2020). ESR1 Gene Variants, Haplotypes and Diplotypes May Influence the Risk of Breast Cancer and Mammographic Density. Mol. Biol. Rep. 47, 8367–8375. doi:10.1007/s11033-020-05823-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Koppula, P., Zhuang, L., and Gan, B. (2020). Cystine Transporter SLC7A11/xCT in Cancer: Ferroptosis, Nutrient Dependency, and Cancer Therapy. Protein Cell 12, 599–620. doi:10.1007/s13238-020-00789-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar Atri, S., Mohan Rana, V., Rana, M., and Ahmed, Z. (2017). Prevalence of Molecular Subtypes of Carcinoma of Breast- a Study from North India. jemds 6, 6204–6207. doi:10.14260/jemds/2017/1349

CrossRef Full Text | Google Scholar

Lai, A., Kahraman, M., Govek, S., Nagasawa, J., Bonnefous, C., Julien, J., et al. (2015). Identification of GDC-0810 (ARN-810), an Orally Bioavailable Selective Estrogen Receptor Degrader (SERD) that Demonstrates Robust Activity in Tamoxifen-Resistant Breast Cancer Xenografts. J. Med. Chem. 58, 4888–4904. doi:10.1021/acs.jmedchem.5b00054

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, X., Green, M. D., Wang, W., Yu, J., Choi, J. E., Jiang, L., et al. (2019). Radiotherapy and Immunotherapy Promote Tumoral Lipid Oxidation and Ferroptosis via Synergistic Repression of SLC7A11. Cancer Discov. 9, 1673–1685. doi:10.1158/2159-8290.CD-19-0338

PubMed Abstract | CrossRef Full Text | Google Scholar

Lei, G., Zhang, Y., Koppula, P., Liu, X., Zhang, J., Lin, S. H., et al. (2020). The Role of Ferroptosis in Ionizing Radiation-Induced Cell Death and Tumor Suppression. Cell Res. 30, 146–162. doi:10.1038/s41422-019-0263-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Xie, P., Lu, L., Wang, J., Diao, L., Liu, Z., et al. (2017). An Integrated Bioinformatics Platform for Investigating the Human E3 Ubiquitin Ligase-Substrate Interaction Network. Nat. Commun. 8, 347. doi:10.1038/s41467-017-00299-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Lian, W., Fu, F., Lin, Y., Lu, M., Chen, B., Yang, P., et al. (2017). The Impact of Young Age for Prognosis by Subtype in Women with Early Breast Cancer. Sci. Rep. 7, 11625. doi:10.1038/s41598-017-10414-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., Huang, J., Liao, X., Yuan, Z., Feng, S., Xie, Y., et al. (2016). Neo-tanshinlactone Selectively Inhibits the Proliferation of Estrogen Receptor Positive Breast Cancer Cells through Transcriptional Down-Regulation of Estrogen Receptor Alpha. Pharmacol. Res. 111, 849–858. doi:10.1016/j.phrs.2016.07.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Liu, R., Liu, Y., Li, G., Chen, Q., Liu, X., et al. (2021). Cystine‐glutamate Antiporter xCT as a Therapeutic Target for Cancer. Cell Biochem. Funct. 39, 174–179. doi:10.1002/cbf.3581

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, S., Fu, X., Liu, L., Liu, Y., Feng, H., Jiang, H., et al. (2021). Iron-Dependent Autophagic Cell Death Induced by Radiation in MDA-MB-231 Breast Cancer Cells. Front. Cell Dev. Biol. 9, 723801. doi:10.3389/fcell.2021.723801

PubMed Abstract | CrossRef Full Text | Google Scholar

Maza, E. (2016). In Papyro Comparison of TMM (edgeR), RLE (DESeq2), and MRN Normalization Methods for a Simple Two-Conditions-Without-Replicates RNA-Seq Experimental Design. Front. Genet. 7, 164. doi:10.3389/fgene.2016.00164

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyers, M. O., Klauber-Demore, N., Ollila, D. W., Amos, K. D., Moore, D. T., Drobish, A. A., et al. (2011). Impact of Breast Cancer Molecular Subtypes on Locoregional Recurrence in Patients Treated with Neoadjuvant Chemotherapy for Locally Advanced Breast Cancer. Ann. Surg. Oncol. 18, 2851–2857. doi:10.1245/s10434-011-1665-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nilsson, S., and Gustafsson, J. Å (2011). Estrogen Receptors: Therapies Targeted to Receptor Subtypes. Clin. Pharmacol. Ther. 89, 44–55. doi:10.1038/clpt.2010.226

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogłuszka, M., Orzechowska, M., Jędroszka, D., Witas, P., and Bednarek, A. K. (2019). Evaluate Cutpoints: Adaptable Continuous Data Distribution System for Determining Survival in Kaplan-Meier Estimator. Comp. Methods Programs Biomed. 177, 133–139. doi:10.1016/j.cmpb.2019.05.023

CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Stockwell, B. R., Friedmann Angeli, J. P., Bayir, H., Bush, A. I., Conrad, M., Dixon, S. J., et al. (2017). Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell 171, 273–285. doi:10.1016/j.cell.2017.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoner, M., Saville, B., Wormke, M., Dean, D., Burghardt, R., and Safe, S. (2002). Hypoxia Induces Proteasome-dependent Degradation of Estrogen Receptor α in ZR-75 Breast Cancer Cells. Mol. Endocrinol. 16, 2231–2242. doi:10.1210/me.2001-0347

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, L., Dong, H., Zhang, W., Wang, N., Ni, N., Bai, X., et al. (2021). Lipid Peroxidation, GSH Depletion, and SLC7A11 Inhibition Are Common Causes of EMT and Ferroptosis in A549 Cells, but Different in Specific Mechanisms. DNA Cel. Biol. 40, 172–183. doi:10.1089/dna.2020.5730

CrossRef Full Text | Google Scholar

Tang, D., Chen, X., Kang, R., and Kroemer, G. (2021). Ferroptosis: Molecular Mechanisms and Health Implications. Cel. Res. 31, 107–125. doi:10.1038/s41422-020-00441-1

CrossRef Full Text | Google Scholar

Wu, Y., Yu, C., Luo, M., Cen, C., Qiu, J., Zhang, S., et al. (2020a). Ferroptosis in Cancer Treatment: Another Way to Rome. Front. Oncol. 10, 571127. doi:10.3389/fonc.2020.571127

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Zhang, S., Gong, X., Tam, S., Xiao, D., Liu, S., et al. (2020b). The Epigenetic Regulators and Metabolic Changes in Ferroptosis-Associated Cancer Progression. Mol. Cancer 19, 39. doi:10.1186/s12943-020-01157-x

CrossRef Full Text | Google Scholar

Xu, X., Zhang, X., Wei, C., Zheng, D., Lu, X., Yang, Y., et al. (2020). Targeting SLC7A11 Specifically Suppresses the Progression of Colorectal Cancer Stem Cells via Inducing Ferroptosis. Eur. J. Pharm. Sci. 152, 105450. doi:10.1016/j.ejps.2020.105450

CrossRef Full Text | Google Scholar

Yang, B., and Kumar, S. (2010). Nedd4 and Nedd4-2: Closely Related Ubiquitin-Protein Ligases with Distinct Physiological Functions. Cell Death Differ. 17, 68–77. doi:10.1038/cdd.2009.84

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Zhou, Y., Xie, S., Wang, J., Li, Z., Chen, L., et al. (2021). Metformin Induces Ferroptosis by Inhibiting UFMylation of SLC7A11 in Breast Cancer. J. Exp. Clin. Cancer Res. 40, 206. doi:10.1186/s13046-021-02012-7

CrossRef Full Text | Google Scholar

Ye, L. F., Chaudhary, K. R., Zandkarimi, F., Harken, A. D., Kinslow, C. J., Upadhyayula, P. S., et al. (2020). Radiation-Induced Lipid Peroxidation Triggers Ferroptosis and Synergizes with Ferroptosis Inducers. ACS Chem. Biol. 15, 469–484. doi:10.1021/acschembio.9b00939

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H. (2012). Typical Cell Signaling Response to Ionizing Radiation: DNA Damage and Extranuclear Damage. Chin. J. Cancer Res. 24, 83–89. doi:10.1007/s11670-012-0083-1

CrossRef Full Text | Google Scholar

Yu, H., Yang, C., Jian, L., Guo, S., Chen, R., Li, K., et al. (2019). Sulfasalazine-induced F-erroptosis in B-reast C-ancer C-ells I-s R-educed by the I-nhibitory E-ffect of E-strogen R-eceptor on the T-ransferrin R-eceptor. Oncol. Rep. 42, 826–838. doi:10.3892/or.2019.7189

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, N., Zhang, P., Huang, H., Liu, W., Hayashi, T., Zang, L., et al. (2015). ERα Down-Regulation Plays a Key Role in Silibinin-Induced Autophagy and Apoptosis in Human Breast Cancer MCF-7 Cells. J. Pharmacol. Sci. 128, 97–107. doi:10.1016/j.jphs.2015.05.001

CrossRef Full Text | Google Scholar

Keywords: ESR1, SLC7A11, ionizing radiation, NEDD4L, ferroptosis, breast cancer

Citation: Liu R, Liu L, Bian Y, Zhang S, Wang Y, Chen H, Jiang X, Li G, Chen Q, Xue C, Li M, Liu L, Liu X and Ma S (2022) The Dual Regulation Effects of ESR1/NEDD4L on SLC7A11 in Breast Cancer Under Ionizing Radiation. Front. Cell Dev. Biol. 9:772380. doi: 10.3389/fcell.2021.772380

Received: 08 September 2021; Accepted: 20 December 2021;
Published: 16 February 2022.

Edited by:

Ruixue Huang, Central South University, China

Reviewed by:

Qiang Liu, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Lijian Shao, Nanchang University, China

Copyright © 2022 Liu, Liu, Bian, Zhang, Wang, Chen, Jiang, Li, Chen, Xue, Li, Liu, Liu and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaodong Liu, liuxd2014@126.com; Shumei Ma, shmm2001@126.com

Download