- Research
- Open access
- Published:
M6A associated TSUC7 inhibition contributed to Erlotinib resistance in lung adenocarcinoma through a notch signaling activation dependent way
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 325 (2021)
Abstract
Background
The small tyrosine kinase inhibitors (TKIs) subversively altered the lung cancer treatments, but patients will inevitably face the therapy resistance and disease recurrence. We aim to explore the potential roles of non-coding RNAs in sensitizing the TKIs effects. Methods: Multiple cellular and molecular detections were applied to confirm the mechanistic regulations and intracellular connections.
Results
We explored the specific gene features of candidates in association with resistance, and found that m6A controlled the stemness of EMT features through METTL3 and YTHDF2. The miR-146a/Notch signaling was sustained highly activated in a m6A dependent manner, and the m6A regulator of YTHDF2 suppressed TUSC7, both of which contributed to the resistant features. Functionally, the sponge type of TUSC7 regulation of miR-146a inhibited Notch signaling functions, and affected the cancer progression and stem cells’ renewal in Erlotinib resistant PC9 cells (PC9ER) and Erlotinib resistant HCC827 cells (HCC827ER) cells. The Notch signaling functions manipulated the cMYC and DICER inner cytoplasm, and the absence of either cMYC or DICER1 lead to TUSC7 and miR-146a decreasing respectively, formed the closed circle to maintain the balance.
Conclusion
PC9ER and HCC827ER cells harbored much more stem-like cells, and the resistance could be reversed by Notch signaling inactivation. The intrinsic miR-146 and TUSC7 levels are monitored by m6A effectors, the alternation of either miR-146 or TUSC7 expression could lead to the circling loop to sustain the new homeostasis. Further in clinics, the combined delivery of TKIs and Notch specific inhibitory non-coding RNAs will pave the way for yielding the susceptibility to targeted therapy in lung cancer.
Background
Lung cancer is the most dangerous cancer type worldwide, topping the cancer related mortality [1]. The lung cancer incidence tends to be stable, and even decreased slightly in western world [1, 2], but situations on cancer prevention were severe in developing China [3,4,5]. Exploring more sensitive screening strategy, improving the radical operation methods, or developing more effective adjuvant therapeutic agents was so urgent than ever [6, 7]. Lung adenocarcinoma consists of lots of therapy targets, and the EGFR Tyrosin Kinase related Inhibitors (TKIs) have been widely and effectively applied in clinical treatments for a decade, shortening the suffering process [8, 9]. However, inevitably, the targeted therapies must face the recurrence, and receive the second or third line of chemo-radiotherapy [10,11,12], and to identify more novel and effective therapeutic molecules and agents will be helpful and promising.
Traditional non-coding genes were always referred to the miRNAs when researching the post-translational regulations, and the lncRNAs (long non-coding RNAs) and circRNAs (circular RNAs) were later revealed for supplementing the affection of RNA family of non-coding members [13,14,15,16,17]. Individually, they could act as crucial modulator toward to certain downstream genes in many ways [18, 19]. The N-6-methyladenosine (m6A) regulation of RNAs attracted lots attention due to its specific and strong modification ability of epigenetic functions [20,21,22]. The m6A controller of writers, readers, and erasers could modulate the mRNA stability and translation, to mediate downstream effects [23, 24]. However, the roles of m6A in controlling the non-coding RNAs biogenesis were still not clear. To further identify the candidates to improve the TKIs treatments sensitization, we tentatively explored the supporting role of TUSC7 in cancer suppression, trying to establish the m6A correlated lncRNA functions in modulating the TKIs therapies resistance.
Materials and methods
Cell lines and cell culture
The human lung adenocarcinoma cells lines (PC9 and HCC827), and human embryonic cell line (HEK-293 T) were purchased from ATCC (American Type Culture Collection, VA, USA) or the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). The Erlotinib-resistant cell lines (PC9/ER and HCC827/ER) were established by chronic exposure to increasing concentrations of drugs. The ectopic expression of TUSC or miR-146 in cells were constructed and restored as pervious study elucidated [25]. All cell lines were cultured in DMEM medium supplemented with 10% FBS (Gibco), 1% penicillin and 1% streptomycin, and incubated at 37 °C in a humidified atmosphere with 5% CO2.
Materials and agents
The Notch inhibitors were chosen as FLI-06 (inhibitor-1, ab120633, Abcam), and γ-Secretase inhibitor (inhibitor-2, ab146170, Abcam). Erlotinib was purchased from Sigma-Aldrich (SML156-50MG, Merck KGaA, Darmstadt, Germany).
Quantitative real-time PCR and western blot
Total RNA was extracted from cells using Trizol, according to the manufacturer’s protocol, and reverse-transcribed into cDNA by using SYBR RT-PCR kit (Takara, JAPAN). Real-time quantitative PCR (RT-qPCR) was performed with SYBR Premix ExTaqTM II Kit (Takara, JAPAN). The sequences of the primers for PCR were synthesized by Sangon Company (Shanghai, China) and were listed in Supplemental Materials and Methods. The relative expression of mRNA and miRNA were calculated by using the formula:2-ΔΔCΤ. For western blot analysis, the total protein from cell extracts was harvested using RIPA buffer contained protease inhibitors. The protein extracts were fractionated by 10% SDS-PAGE, transferred onto a nitrocellulose membrane, and then incubated with primary antibodies at 4 °C overnights, followed by HRP-conjugated secondary antibody (1:5000, #7074, Cell Signaling Technology) and visualized by using ECL Blotting Detection Reagents (Merck Millipore). The primary antibodies were as follows: anti-Notch1 (1:500; Rabbit mAb, #3608, Cell Signaling Technology), anti-Notch2 (1:1000, ab8926, Abcam), anti-NCSTN (1:1000, ab189125, Abcam), anti-Vinculin (1:8000, #4650, Cell Signaling Technology), anti-NUMB (1:1500, #2761, Cell Signaling Technology), anti-EGFR (1:2000, #4267, Cell Signaling Technology), anti-DICER1 (1:500, #5362, Cell Signaling Technology), anti-CMYC (1:1500, (9E10): sc-40, Santa Cruz), anti- Snail (1:1000, ab31787, Abcam), anti-METTL3 (1:2000, ab240595, Abcam), anti-YTHDF2 (1:1000, EPR20318, ab220163, Abcam).
Sphere formation assay
Single-cell suspensions (1000 cells per well) of different groups were plated on six-well ultralow attachment plates (Corning Incorporated) in serum-free DMEM/F12 Medium supplemented with 20 ng/mL EGF (Invitrogen, Carlsbad, CA), b-FGF (Invitrogen, Carlsbad, CA) and 4 μg/mL insulin (Sigma-Aldrich, St. Louis, MO). After 1 week culture, the spheres of>50 μm was quantified by using an inverted microscope.
Dual luciferase report assay
The putative sequences or mutant sequences of miR-146a target sites for TUSC7–3’UTR was synthesized and cloned into the pGL3 reporter vector (Promega). These constructed reporters were named pGL3-TUSC7-WT, pGL3-TUSC7-MUT. For luciferase assay, the cells were seeded onto 24-well plates and co-transfected with 200 ng of pGL3-TUSC7/EGFR-WT or pGL3-TUSC7/EGFR -MUT, 20 ng of pRL-TK plasmid as normalization control, together with miR-146a-5p mimic or miR-146a-5p control (GenePharma, Shanghai, China). After 48 h of transfection, Luciferase assays were performed by using the Dual Luciferase Reporter Assay System (Promega, WI, USA).
ALDEFLUOR assay and fluorescence-activated cell sorting (FACS)
Aldehyde dehydrogenase (ALDH) enzyme activity in lung cells was determined by ALDEFLUOR assay kit (Stem Cell Technologies, Grenoble, France) according to the manufacturer’s instructions. Briefly, 1 × 106/ml cells mixed with Aldefluor® assay buffer containing 1.5μΜ bodipy-aminoacetaldehyde (BAAA, an ALDH substrate). Then, the cell/substrate mixture was incubated for 1 h at 37 °C. Diethylaminobenzaldehyde (DEAB), a specific ALDH1 enzyme inhibitor, was used as negative control. The ALDH + population was detected in the green fluorescence channel (520–540 nm) of FACSAria (Becton Dickinson). Stained cells analyzed and sorted by utilizing FACSDiva (BD Biosciences) and Flow-Jo software (Treestar, Ashland, OR). Nonviable cells were excluded using 1 μg propidium iodide (PI; Sigma Aldrich, Vienna, Austria).
RNA-immunoprecipitation
The cells (5 × 106)/ml were harvested and washed twice with ice-cold 1 × PBS buffer. Collected cell pellet was lysed for 15 min on ice by RIP buffer for 30 min, and pretreated with a 1:10 dilution in NT2 buffer. The cell lysate was further centrifuged at 15,000×g for 15 min at 4 °C, followed by treated with magnetic beads conjugated to human anti-ATG3 antibody (1:50) or the control IgG for 18 h at 4 °C and further washed twice with cold NT2 Buffer. Magnetic beads subsequently mixed with the diluted lysates (10 μl beads/ml lysate). The mixture was re-suspended in 100 μl NT2 Buffer containing 30 μg proteinase K to digest the protein. Co-purified RNA was extracted by the TRIzol reagent and used in subsequent qRT-PCR assay.
M6A methylation quantification assay
The m6A methylation status of cells detected using the m6A RNA Methylation Quantification Kit (Epigentek, Cat#P-9005-113) according to the manufacturer’s instructions. In brief, 200 ng of total RNA was used as an input respectively. Then RNA samples were captured and detected by spectrophotometer (Bio Tek Instruments, Inc. US) at 450 nm. The level of m6A methylation was calculated according to the manufacturer’s instructions.
As for LC–MS/MS assay, mRNA was purified from the total RNA using via oligo dT magnetic beads. Then 200 ng mRNAs were incubated with 0.5 U nuclease P1 in reaction system at 42 °C for about 1 h. After that, mRNAs were incubated with 3 μL of 1 M NH4HCO3 and 1 μL of 1 U/μL alkaline phosphatase at 37 °C for 2 h. And mRNAs were diluted and filtered. After which, a C18 column were used to separate mRNAs [26,27,28,29]. Then mRNAs were analyzed by an Agilent (6410 QQQ) triple-quadrupole LC mass spectrometer. Calibration curves were used to calculate the Ratio of m6A to A.
Nude mouse xenograft model
A total of 21 4-week-old female BALB/cA-nu nude mice were purchased from Beijing Huafukang Biosciences (Beijing, China), then we maintained them in specific pathogen-free conditions. Control vector, TUSC7 knockout, FLI-06 treated H1975 cells (1*107) cells were suspended in 100 μL of serum-free DMEM medium (Hyclone, USA), mixed with matrix gel (Corning, USA), and then were injected subcutaneously. The changes in the tumor size were recorded every 3 or 5 days. We calculated the tumor volume using the following formula: V = 1/2 × l × w2 (l is the longer axis, 2 is the shorter axis). All mice were sacrificed 28 days after the injection of cells. The dissected tumor samples were immersed in 4% paraformaldehyde (BioSharp, China) and embedded in paraffin.
Statistical analysis
Statistical analysis was carried out by using Graph Pad Prism 6 and SPSS 20.0 software (SPSS Inc., Chicago, IL, USA). All numerical data were expressed as mean ± standard deviation (SD). Experiments were carried out with three or more replicates. Two or more groups were assessed by using Student’s t test or ANOVA individually. P < 0.05 was considered to be statistically significant.
Results
Notch signaling activation in lung adenocarcinoma pointed to poor survival expectance
The important members of Notch signaling [30, 31] were screened for expression patterns with using Pan-Cancer Atlas of TCGA data base. The total samples of 507 patients were collected, and the heat maps indicated the universal overexpression of notch signaling participants (Fig. 1A-B), and EGFR was correlated with aberrant Notch expression. Deep analysis from TCGA (Nat Genet 2016 data base) indicated the grouped enrichment of Notch signaling factors (Fig. 1C), and the changing was consistency in groups harboring most irregulated Notch functions (Fig. 1D).
The overall survival (OS) and Disease/Progression-free survival (DFS/PFS) data were acquired from TCGA data at the CBIOPORTAL FOR CANCER GENOMIC (http://www.cbioportal.org/) [32, 33]. Notch signaling activation decreased the survival time (Fig. 1E), and patients tended to bear relapse or resistance in shorter follow-up periods (Fig. 1F). The overall survival and first progression estimates were analyzed by applying Kaplan-Meier analysis (http://kmplot.com/analysis/) [34,35,36]. The highly expressed TUSC7 indicated better progression-free estimates in adenocarcinoma, comparing to the lower expressed groups (Fig. 1G). The functional EGFR signaling could transcriptionally activate multiple downstream pathways, and the positive relationship between miR-146a and EGFR was found (Fig. 1H). MiR-146 also correlated with Notch signaling factors expressions (Fig. 1I-J).
Erlotinib resistant cells harboring notch signaling activation and TUSC7 inhibition
Lung adenocarcinoma cell lines of PC9 and HCC827 were selected for their characteristics of specific EGFR mutant status, and the signatures of PC9 and HCC827 cells referring to erlotinib treatment were shown in Figure S1. The concentration of 0.2 μM Erlotinib was chosen as the function candidate, which decreased the Notch1 (Figure S2A) and Notch2 (Figure S2B) mRNA levels in PC9 and HCC827 cells. The Erlotinib treatment also increased the TUSC7 expression level significantly (Figure S2C), together with EGFR level slightly decreased (Figure S2D). Blotting results showed that the Erlotinib treatment decreased the Notch signaling factors in PC9 and HCC827 cells (Figure S2E).
We constructed Erlotinib resistant lung cancer cells with advancing concentration gradient (Figure S3), and the differentially expressed lncRNAs between Erlotinib sensitive and resistant cells were primarily detected by gene panel selection (Figure S2F). The differentially expressed lncRNAs were pasted for GO analysis of functional identification (Figure S2G). The detailed information showed that TUSC7 of both PC9ER (Figure S2H) and HCC827ER (Figure S2I) decreased greatly, and did not react to Erlotinib treatment, comparing to the Erlotinib sensitive cells. Last but not least Erlotinib failed to inhibit the Notch signaling functions in PC9ER and HCC827ER cells (Figure S2H-I).
Notch signaling inhibition was required for TUSC7 alleviating of the Erlotinib resistance
The stem-like cells were accused for treatment resistance, and to study the stem cells’ renewal ability in contributing to Erlotinib resistance, the ALDH1 phenotype and spheres formation nature were applied. The Erlotinib resistant PC9ER and HCC827ER cells consisted of more stem-like cells (Fig. 2A-B), and Erlotinib treatment failed to decrease the stem cells number (Fig. 2C). Higher stem cells ratio indicated resistant status, and resulted in group resistance. The Notch inhibitors were chosen as FLI-06 (inhibitor-1, ab120633, Abcam), and γ-Secretase inhibitor (inhibitor-2, ab146170, Abcam). Both PC9 and HCC827 cells responded to Notch signaling inhibitors, with stem-like cells ratios decreasing significantly (Fig. 2D-E). The lowered concentration of Notch signaling inhibitors sensitized resistant cells to Erlotinib treatment (Fig. 2F, Figure S4A-B), indicating the synergistic effects of TKI agents and Notch signaling inhibitors.
To study the TUSC7 functions, Lentiviral based TUSC7 expression vesical was introduced into PC9ER and HCC827ER cells. Erlotinib alone inhibited the Notch signaling slightly, but TUSC7 effectively enhanced the Erlotinib functions in PC9ER (Fig. 2G) and HCC827ER cells (Fig. 2H), and stimulated the suppressive functions of Erlotinib in both PC9ER (Figure S4C) and HCC827ER cells (Figure S4D). Moreover, the addition of lowered concentration of Notch signaling inhibitor strengthened TUSC7 functions (Fig. 2G-H).
Next, we found that combined TUSC7 and Erlotinib decreased the stem cells ratio greatly in both PC9ER (Fig. 2I) and HCC827ER cells (Fig. 2J). The stem cells’ renewal suppression evaluation did not show significant differences between TUSC7 alone and the combination of TUSC7 and FLI-06 (Fig. 2K-L, Figure S4E-G). The combined TUSC7 and Erlotinib decreased the stem cells associated signatures, decreasing EMT markers in PC9ER cells (Figure S5A) and HCC827ER cells (Figure S5B).
TUSC7 sensitization of Erlotinib through miR-146a/notch signaling inhibition was dependent on NUMB restoration
To reveal the transduction mechanisms, bioinformatic screening of the possible connections between TUSC7 and its binding partners was conducted. We noticed that miR-146a shared common sequences with the untranslated regions of TUSC7 (Fig. 3A), and miR-146a mimics decreased the Luc-activity of TUSC7 significantly (Fig. 3B) in 293 T cells, which was also confirmed in PC9ER (Fig. 3C) cells and HCC827ER cells (Fig. 3D). To testify the binding probabilities between TUSC7 and its downstream partner, biotin labeled sense and anti-sense RNAs of TUSC7 were used for RNA pull down detection, and the connection between TUSC7 and NUMB in resistant lung adenocarcinoma cells was identified (Fig. 3E). Further, RNA immunoprecipitation revealed that TUSC7 was enriched with NUMB expression in PC9ER cells (Fig. 3F, left) and HCC827ER cells (Fig. 3F right). Informatic screening of the potential miRNAs’ targets suggested that miR-146a may bind to NUMB (Fig. 3G), and NUMB decreased greatly in cancer group (Fig. 3H). The alignment of miR-146a and the 3’UTR of NUMB was constructed through using the enhanced green fluorescent protein (EGFP) reporter assay. The wild-type 3’UTR sequence and the mutant 3’UTR sequence of NUMB were cloned downstream from the EGFP-coding sequence respectively, to construct the reporter plasmid and the mutant vector. The co-transfection of miR-146a mimics and the wide-type reporter plasmid strongly reduced the EGFP intensity (Fig. 3I, left), but not happened in mutant-type reporter plasmid (Fig. 3I, right). On the contrary, TUSC7 did not reduce the EGFP activity of the NUMB (Fig. 3J). MiR-146a decreased the NUMB expression, which could be rescued by TUSC7, and the TUSC7 inhibition (TUSC7-in) also decreased the NUMB at the protein level (Fig. 3K).
m6A in resistant cells contributed to TUSC7 inhibition and miR-146a overexpression
The pluripotency status of resistant PC9ER and HCC827ER cells contributed to specific miR-146a and TUSC7 patterns, and to characterize the roles of m6A in therapy resistance, we investigated the variations of m6A levels, and identified that the m6A levels of total RNAs from resistant cells were statistically more abundant than sensitive original cells by using LC/MS (Fig. 4A). To further characterize the roles of m6A in generating the resistance, we used siRNAs to tentatively test the m6A related processers in controlling of TUSC7 and miR-146a. METTL3 affected the miR-146a level (Fig. 4B), and YTHDF2 affected TUSC7 level (Fig. 4C),and results were all confirmed by using the lentiviral based METTL3/YTHDF2 knock-down systems (Fig. 4D-E). Dysregulated METTL3 (Fig. 4F) and YTHDF2 (Fig. 4G) affected the EMT and pluripotency features through controlling m6A status. Moreover, decreasing either METTL3 or YTHDF2 significantly altered Snail and EMT factors expressions. m6A at Snai1 mRNA was greatly increased in resistant cells [26, 37], and we found the METTL3 inhibition decreased m6A at Snai1(Fig. 4H-I), which later failed to activate the miR-146a promoter (Fig. 4J) [38]. In the other hand, the m6A at TUSC7 level increased in resistant cells (Fig. 4K), and the recognition of TUSC7 m6A peak by YTHDF2 degraded TUSC7 [39, 40]. The Me-RIP (methylated RNA immunoprecipitation) assay confirmed that the high abundance of m6A modification in cells with YTHDF2 inhibition (Fig. 4L). All the results revealed that m6A determined the expression levels of miR-146a and TUSC7 in resistant cells, sustaining the balanced status.
TUSC7 formed feedback loop with miR-146/notch signaling and lead the Erlotinib re-sensitization
The unique m6A at different RNA sequences determined intrinsic RNA expressions, cellular homeostasis may be disturbed when manipulating certain lncRNA candidate. CMYC was proved previously to promote the translation of lncRNAs [41,42,43], and DICER1 was critical for miRNAs maturation [44,45,46]. Results showed that in PC9ER and HCC827ER cells, Notch signaling inhibition lead to cMYC promoter activity decreasing, and TUSC7 inhibition strongly increased the activity (Fig. 5A). Inhibition of TUSC7 increased miR-146a expression, and then stimulated the DICER1 activity (Fig. 5B).
The scheme image represented the cMYC promoter detection modes (Fig. 5C), and the mutation blocked the EGFR/Notch signaling feedback loop on lncRNA transcription induction. Western blotting further confirmed that in PC9ER and HCC827ER cells, miR-146 and TUSC7 acted the opposite way to promote the DICER1/CCND1 expression (Fig. 5D). TUSC7 and miR-146a was sustained by cMYC (Fig. 5E, left) and DICER1 (Fig. 5E, right) respectively. The regulative signaling pathways were drafted and illustrated for detailed exhibition (Fig. 5F), the dysregulated TUSC7 and miR-146a formed the feedback loop with their downstream effectors to sustain the new homeostasis.
In vivo study confirmed the effective inhibition of TUSC7 exhibited on tumor growth, and the Notch signaling inactivation by using FLI-06 also suppressed the in vivo tumor expansion (Fig. 5G), further proved its suppressive functions.
Conclusion
Lung cancer treatments have been refined greatly, with constantly and novelly emerging components, and the small molecular compounds improved the life quality and expected response to therapies in patients carrying sensitive EGFR mutants [47,48,49]. Surgery manners have been evolved continuously with little operation wound and duration, however, the improvement has brought little progress on overall survival [50,51,52]. The 1st generation of small tyrosine kinase inhibitors subversively altered the lung cancer treatments, and the “small step” has made the “major progress”, increasing the PFS with little suffering.
Reality is cruel. No matter how ideal results the TKIs treatment has brought, patients with advanced lung cancer will inevitably face the therapy resistance and disease recurrence [53,54,55]. To improve the TKIs sensitization and to overcome resistance will be the headline goal in the current emergency [56]. Many factors were thought to be correlated with 1st generation of TKI functions’ achieving, and certain genes and compounds sensitized the treatments by oncogenic signaling repression, but detailed mechanistic regulations have not been illustrated clearly ever [10, 57].
We have been focusing on noncoding RNAs related cancer emergence, progression, stemness features and therapy response for decades, and have revealed some important factors that may contribute to [44, 45, 58]. Previously, we identified the miR-146 functions in lung adenocarcinoma with suppressive affections on stem cells’ renewal, and the tentative exploration of TKIs treatment improvements will facilitate its practical application.
In this study, we identified the inhibitive roles of TUSC7 in lung cancer progression, and after acquiring the Erlotinib resistant cells, gene panel was used for massive assessing of the dysregulated non-coding RNAs. We first explored the features of candidate genes in association with resistance, and the activated miR-146a/Notch signaling was sustained in resistant stem cells in a m6A dependent manner of and METTL3/Snail cascade. M6A related YTHDF2 mediation of suppressed TUSC7 also contributed to resistant features. Functionally, m6A controlled the stemness of EMT features through METTL3 and YTHDF2 in resistant PC9ER and HCC827ER cells.
In detail, TUSC7 sensitized the Erlotinib effects, and decreased the stem cells ratio through Notch signaling inhibition. Bench study showed that, TUSC7 sponged to miR-146 and then released NUMB to control Notch signaling, the latter of which was critical for maintaining cancer stem cells (CSCs) pool. DICER1 and cMYC activity was critical for sustaining the non-coding RNAs maturation [41, 43, 46, 59, 60], the sponge style of TUSC7 regulation toward to miR-146 released the controlling of NUMB expression in PC9ER and HCC827ER cells, which manipulated the DICER and cMYC cascade inner cytoplasm [43, 46, 61, 62]. The absence of either cMYC or DICER1 will lead to TUSC7 and miR-146 decreasing respectively, formed the closed circle to maintain the new balance.
To be concluded, the lncRNA of TUSC7 affected the cancer progression and stem cells renewal, and TUSC7 suppression of Notch signaling determines the Erlotinib treatment response. PC9ER and HCC827ER cells harbors much more stem-like groups, which dominated in therapy response, and their resistance could be reversed by Notch signaling inactivation. Interestingly and importantly, the intrinsic miR-146 and TUSC7 levels are monitored and sustained by m6A effectors, and disturbing the miR-146 and TUSC7 expression patterns will push themselves to form the circling loop to sustain the new homeostasis. Further in clinics, the combined using of TKIs and Notch specific inhibitory non-coding RNAs will pave the way for yielding the susceptibility to targeted therapy in lung cancer.
Availability of data and materials
The datasets during and/or analyzed during the current study available from the corresponding author on reasonable request.
Change history
24 July 2023
A Correction to this paper has been published: https://doi.org/10.1186/s13046-023-02760-8
Abbreviations
- TKIs:
-
Tyrosin Kinase related Inhibitors
- lncRNAs:
-
long non-coding RNAs
- circRNAs:
-
circular RNAs
- m6A:
-
N-6-methyladenosine
- RT-qPCR:
-
Real-time quantitative PCR
- ALDH:
-
Aldehyde dehydrogenase
- BAAA:
-
Bodipy-aminoacetaldehyde
- DEAB:
-
Diethylaminobenzaldehyde
- EGFP:
-
Green fluorescent protein
- Me-RIP:
-
Methylated RNA immunoprecipitation
References
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(6):394–424.
Bogos K, Kiss Z, Gálffy G, Tamási L, Ostoros G, Müller V, et al. Revising Incidence and Mortality of Lung Cancer in Central Europe: An Epidemiology Review From Hungary. Front Oncol. 2019;9:1051.
Zheng R, Zeng H, Zuo T, Zhang S, Qiao Y, Zhou Q, et al. Lung cancer incidence and mortality in China, 2011. Thorac Cancer. 2016;7(1):94–9.
Chen WQ, Zuo TT, Zheng RS, Zeng HM, Zhang SW, He J. Lung cancer incidence and mortality in China in 2013. Zhonghua Zhong Liu Za Zhi. 2017;39(10):795–800.
Raaschou-Nielsen O, Andersen ZJ, Beelen R, Samoli E, Stafoggia M, Weinmayr G, et al. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European study of cohorts for air pollution effects (ESCAPE). Lancet Oncol. 2013;14(9):813–22.
Cai W-Y, Wei T-Z, Luo Q-C, Wu Q-W, Liu Q-F, Yang M, et al. The Wnt–β-catenin pathway represses let-7 microRNA expression through transactivation of Lin28 to augment breast cancer stem cell expansion. J Cell Sci. 2013;126(13):2877–89.
Chang W, Wang H, Kim W, Liu Y, Deng H, Liu H, et al. Hormonal Suppression of Stem Cells Inhibits Symmetric Cell Division and Gastric Tumorigenesis. Cell Stem Cell. 2020;26(5):739–54.e8.
Sini C, Tuzi A, Rossi G, Russo A, Pezzuto A. Acquired resistance in oncogene-addicted non-small-cell lung cancer. Future Oncol. 2018;14(13s):29–40.
Jeong H, Phan ANH, Choi JW. Anti-cancer Effects of Polyphenolic Compounds in Epidermal Growth Factor Receptor Tyrosine Kinase Inhibitor-resistant Non-small Cell Lung Cancer. Pharmacogn Mag. 2017;13(52):595–9 PMCID: 5701397.
Wang H, Lu B, Ren S, Wu F, Wang X, Yan C, et al. Long noncoding RNA LINC01116 contributes to Gefitinib resistance in non-small cell lung Cancer through regulating IFI44. Mol Ther Nucleic Acids. 2020;19:218–27.
Shu D, Xu Y, Chen W. Knockdown of lncRNA BLACAT1 reverses the resistance of afatinib to non-small cell lung cancer via modulating STAT3 signalling. J Drug Target. 2020;28(3):300–6.
He J, Jin S, Zhang W, Wu D, Li J, Xu J, et al. Long non-coding RNA LOC554202 promotes acquired gefitinib resistance in non-small cell lung cancer through upregulating miR-31 expression. J Cancer. 2019;10(24):6003–13 PMCID: 6856583.
Wang M, Yu F, Wu W, Zhang Y, Chang W, Ponnusamy M, et al. Circular RNAs: a novel type of non-coding RNA and their potential implications in antiviral immunity. Int J Biol Sci. 2017;13(12):1497–506.
Liu Y, Cheng Z, Pang Y, Cui L, Qian T, Quan L, et al. Role of microRNAs, circRNAs and long noncoding RNAs in acute myeloid leukemia. J Hematol Oncol. 2019;12(1):51.
Sulaiman SA, Muhsin NIA, Jamal R. Regulatory Non-coding RNAs Network in Non-alcoholic Fatty Liver Disease. Front Physiol. 2019;10:279.
Pang Y, Mao C, Liu S. Encoding activities of non-coding RNAs. Theranostics. 2018;8(9):2496–507.
Li X, Wang M, Du N, Liang T, Xiao G-D, Li K, et al. Matrine inhibitory effect on self-renewal and re-sensitization of 5-FU resistant NSCLC stem cells were through let-7b dependent Downregulation of CCND1. Cell Cycle. 2020;19(23):3249–59.
Shabaninejad Z, Vafadar A, Movahedpour A, Ghasemi Y, Namdar A, Fathizadeh H, et al. Circular RNAs in cancer: new insights into functions and implications in ovarian cancer. J Ovarian Res. 2019;12(1):84.
Zhou R, Chen KK, Zhang J, Xiao B, Huang Z, Ju C, et al. The decade of exosomal long RNA species: an emerging cancer antagonist. Mol Cancer. 2018;17(1):75.
Xie F, Huang C, Liu F, Zhang H, Xiao X, Sun J, et al. CircPTPRA blocks the recognition of RNA N6-methyladenosine through interacting with IGF2BP1 to suppress bladder cancer progression. Mol Cancer. 2021;20(1):68.
Rickert D, Bartl J, Picard D, Bernardi F, Qin N, Lovino M, et al. Circular RNA profiling distinguishes medulloblastoma groups and shows aberrant RMST overexpression in WNT medulloblastoma. Acta Neuropathol. 2021;141(6):975–8.
Huang H, Weng H, Chen J. m(6)A Modification in Coding and Non-coding RNAs: Roles and Therapeutic Implications in Cancer. Cancer Cell. 2020;37(3):270–88 PMCID: 7141420.
Zhang H, Xiao X, Wei W, Huang C, Wang M, Wang L, et al. CircLIFR synergizes with MSH2 to attenuate chemoresistance via MutSα/ATM-p73 axis in bladder cancer. Mol Cancer. 2021;20(1):70.
Si W, Li Y, Ye S, Li Z, Liu Y, Kuang W, et al. Methyltransferase 3 Mediated miRNA m6A Methylation Promotes Stress Granule Formation in the Early Stage of Acute Ischemic Stroke. Front Mol Neurosci. 2020;13:103 PMCID: 7289951.
Huang G, Wang M, Li X, Wu J, Chen S, Du N, et al. TUSC7 suppression of notch activation through sponging MiR-146 recapitulated the asymmetric cell division in lung adenocarcinoma stem cells. Life Sci. 2019;232:116630.
Lin X, Chai G, Wu Y, Li J, Chen F, Liu J, et al. RNA m6A methylation regulates the epithelial mesenchymal transition of cancer cells and translation of Snail. Nat Commun. 2019;10(1):2065.
Cho S, Lee G, Pickering BF, Jang C, Park JH, He L, et al. mTORC1 promotes cell growth via m6A-dependent mRNA degradation. Mol Cell. 2021;81(10):2064–75.e8.
Xu F, Li J, Ni M, Cheng J, Zhao H, Wang S, et al. FBW7 suppresses ovarian cancer development by targeting the N6-methyladenosine binding protein YTHDF2. Mol Cancer. 2021;20(1):45.
Wang M, Liu J, Zhao Y, He R, Xu X, Guo X, et al. Upregulation of METTL14 mediates the elevation of PERP mRNA N6 adenosine methylation promoting the growth and metastasis of pancreatic cancer. Mol Cancer. 2020;19(1):130.
Ito T, Kudoh S, Ichimura T, Fujino K, Hassan WA, Udaka N. Small cell lung cancer, an epithelial to mesenchymal transition (EMT)-like cancer: significance of inactive notch signaling and expression of achaete-scute complex homologue 1. Hum Cell. 2017;30(1):1–10.
Nowell CS, Radtke F. Notch as a tumour suppressor. Nat Rev Cancer. 2017;17(3):145–59.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):pl1 PMCID: 4160307.
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.
Lánczky A, Nagy Á, Bottai G, Munkácsy G, Szabó A, Santarpia L, et al. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.
Szasz AM, Lanczky A, Nagy A, Forster S, Hark K, Green JE, et al. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget. 2016;7(31):49322–33 PMCID: 5226511.
Győrffy B, Bottai G, Lehmann-Che J, Kéri G, Őrfi L, Iwamoto T, et al. TP53 mutation-correlated genes predict the risk of tumor relapse and identify MPS1 as a potential therapeutic kinase in TP53-mutated breast cancers. Mol Oncol. 2014;8(3):508–19.
Pan Z, Cai J, Lin J, Zhou H, Peng J, Liang J, et al. A novel protein encoded by circFNDC3B inhibits tumor progression and EMT through regulating Snail in colon cancer. Mol Cancer. 2020;19(1):71.
Hwang WL, Jiang JK, Yang SH, Huang TS, Lan HY, Teng HW, et al. MicroRNA-146a directs the symmetric division of snail-dominant colorectal cancer stem cells. Nat Cell Biol. 2014;16(3):268–80.
Liu J, Zhang X, Chen K, Cheng Y, Liu S, Xia M, et al. CCR7 Chemokine Receptor-Inducible lnc-Dpf3 Restrains Dendritic Cell Migration by Inhibiting HIF-1α-Mediated Glycolysis. Immunity. 2019;50(3):600–15.e15.
Wang X, Zhao Boxuan S, Roundtree Ian A, Lu Z, Han D, Ma H, et al. N6-methyladenosine modulates messenger RNA translation efficiency. Cell. 2015;161(6):1388–99.
Arman K, Möröy T. Crosstalk between MYC and lncRNAs in hematological malignancies. Front Oncol. 2020;10:1–15.
Wang C, Yang Y, Zhang G, Li J, Wu X, Ma X, et al. Long noncoding RNA EMS connects c-Myc to cell cycle control and tumorigenesis. Proc Natl Acad Sci U S A. 2019;116(29):14620–9 PMCID: 6642410.
Tokgun PE, Tokgun O, Kurt S, Tomatir AG, Akca H. MYC-driven regulation of long non-coding RNA profiles in breast cancer cells. Gene. 2019;714:143955.
Wang M, Li Y, Xiao GD, Zheng XQ, Wang JC, Xu CW, et al. H19 regulation of oestrogen induction of symmetric division is achieved by antagonizing let-7c in breast cancer stem-like cells. Cell Prolif. 2019;52(1):e12534.
Sun X, Xu C, Tang SC, Wang J, Wang H, Wang P, et al. Let-7c blocks estrogen-activated Wnt signaling in induction of self-renewal of breast cancer stem cells. Cancer Gene Ther. 2016;23(4):83–9.
Sun X, Tang S-C, Xu C, Wang C, Qin S, Du N, et al. DICER1 regulated let-7 expression levels in p53-induced cancer repression requires cyclin D1. J Cell Mol Med. 2015;19(6):1357–65.
Liu Y. Small cell lung cancer transformation from EGFR-mutated lung adenocarcinoma: A case report and literatures review. Cancer Biol Ther. 2018;19(6):445–9 PMCID: 5927699.
Herbst RS, Morgensztern D, Boshoff C. The biology and management of non-small cell lung cancer. Nature. 2018;553(7689):446–54.
Soria J-C, Ohe Y, Vansteenkiste J, Reungwetwattana T, Chewaskulyong B, Lee KH, et al. Osimertinib in untreated EGFR-mutated advanced non–small-cell lung Cancer. N Engl J Med. 2017;378(2):113–25.
Song S-W, Lee H-S, Kim MS, Lee JM, Kim JH, Nam B-H, et al. Preoperative serum fibrinogen level predicts postoperative pulmonary complications after lung Cancer resection. Ann Thorac Surg. 2006;81(6):1974–81.
Wao H, Mhaskar R, Kumar A, Miladinovic B, Djulbegovic B. Survival of patients with non-small cell lung cancer without treatment: a systematic review and meta-analysis. Syst Rev. 2013;2:10.
Oshima Y, Tanimoto T, Yuji K, Tojo A. EGFR-TKI-Associated Interstitial Pneumonitis in Nivolumab-Treated Patients With Non-Small Cell Lung Cancer. JAMA Oncol. 2018;4(8):1112–5 PMCID: 5885195.
Li H, Hu H, Wang R, Pan Y, Wang L, Li Y, et al. Primary concomitant EGFR T790M mutation predicted worse prognosis in non-small cell lung cancer patients. Onco Targets Ther. 2014;7:513–24.
Rotow J, Bivona TG. Understanding and targeting resistance mechanisms in NSCLC. Nat Rev Cancer. 2017;17(11):637–58.
Wu SG, Shih JY. Management of acquired resistance to EGFR TKI-targeted therapy in advanced non-small cell lung cancer. Mol Cancer. 2018;17(1):38 PMCID: 5817870.
Sun X, Li K, Zhao R, Sun Y, Xu J, Peng Z-Y, et al. Lung cancer pathogenesis and poor response to therapy were dependent on driver oncogenic mutations. Life Sci. 2021;265:118797.
Sousa AC, Silveira C, Janeiro A, Malveiro S, Oliveira AR, Felizardo M, et al. Detection of rare and novel EGFR mutations in NSCLC patients: implications for treatment-decision. Lung Cancer. 2020;139:35–40.
Xu C, Xiao G, Zhang B, Wang M, Wang J, Liu D, et al. CCAT1 stimulation of the symmetric division of NSCLC stem cells through activation of the Wnt signalling cascade. Gene Ther. 2018;25(1):4–12.
Wang C, Yang Y, Zhang G, Li J, Wu X, Ma X, et al. Long noncoding RNA EMS connects c-Myc to cell cycle control and tumorigenesis. (1091–6490 (Electronic)).
Xiao G, Li X, Li G, Zhang B, Xu C, Qin S, et al. MiR-129 blocks estrogen induction of NOTCH signaling activity in breast cancer stem-like cells. Oncotarget. 2017;8(61):103261–73 PMCID: 5732725.
Lü J, Zhao Q, Ding X, Guo Y, Li Y, Xu Z, et al. Cyclin D1 promotes secretion of pro-oncogenic immuno-miRNAs and piRNAs. Clin Sci (Lond). 2020;134(7):791–805.
Yu Z, Wang L, Wang C, Ju X, Wang M, Chen K, et al. Cyclin D1 induction of Dicer governs microRNA processing and expression in breast cancer. Nat Commun. 2013;4(1):2812.
Acknowledgments
The authors acknowledge assistants in the Center for Translational Medicine, Cancer Center of the First Affiliated Hospital of Xi’an Jiaotong University, and authors appreciated all the staff of the Thoracic Department, for their technical assistance. We also appreciated very much for Dr. Mallory Ellis and Dr. Leslie Robinson for her careful revision on words editing and language polishing.
Funding
The research leading to these results received funding from National Science Foundation for Young Scientists of China, grant No. 81602597 (Referred to Xin Sun), National Science Foundation for Young Scientists of China, Grant No. 82003140 (Referred to Guo-Dong Xiao), Foundation Research Project of Shaanxi Province, 2021SF-117 (Referred to Xin Sun), The Natural Science Basic Research Program of Shaanxi Province, No. 2018JM7017 (Referred to Xin Sun). This work was also partly supported by the Scientific and Technological Planning project of Shaanxi Province, No. 2020JM-403 (Referred to Jing Zhang).
Author information
Authors and Affiliations
Contributions
KL: Bioscientific experiments; Cells culturing, Experimental tests, draft revision; RNA/Protein tests. ZP: Bioscientific experiments; Study designation. QW: Experimental tests; Cells culturing. XL: Bioscientific experiments; Study designation. GX: Paper drafting; RNA/Protein tests; Statistical analysis; Data base screening; Figures preparation. JZ: RNA/Protein tests. HR: Study designation; Statistical analysis. ST: Paper drafting; Study designation; Statistical analysis; Images quality control, draft revision. XS: Paper drafting; Study designation; Statistical analysis; Data base screening; Figures preparation; References cross checking, draft revision. The authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
All procedures performed in studies involving animals were in accordance with the ethical standards of the institution or practice at which the studies were conducted. We stated that the protocol for the research project has been approved by the Ethics Committee of the First Affiliated Hospital of Xi’an Jiaotong university, and that it conforms to the provisions of the Declaration of Helsinki.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests” in this section.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: an error was identified Figure 2 and 4. The original alternative image was mistakenly used in Figure 2H, while in the editing of Figure 4G, one Western blot band was not carefully discerned and was duplicated.
Supplementary Information
Additional file 1: Supplemental Figure 1.
The sensitivity of Erlotinib in treating Lung adenocarcinoma harboring mutant EGFR. The sensitivity analysis was carried with using shared data of Genomics of Drug Sensitivity in Cancer at the SANGER site. A. The Erlotinib sensitivity referring to lung cancer samples of PAN data were rankly exhibited. B. The Erlotinib sensitivity referring to lung adenocarcinoma samples were rankly exhibited. Both PC9 (C) and HCC827 (D) were sensitive to Erlotinib treatment with concentration much lower than IC50. E. EGFR mutant lung cancer cells are very sensitive to Erlotinib treatment, comparing to that of lung cancer cells with wild type EGFR.
Additional file 2: Supplemental Figure 2.
Signatures of the Erlotinib resistant lung adenocarcinoma cells. A. Erlotinib decreased the Notch1 mRNA level in both PC9 and HCC827 cells. B. Erlotinib decreased the Notch2 mRNA level in both PC9 and HCC827 cells. C. Erlotinib increased the TUSC7 expression level significantly in both PC9 and HCC827 cells. D. Addition of Erlotinib did not change the EGFR expression level in both PC9 and HCC827 cells. E. Erlotinib decreased the Notch signaling factors in PC9 and HCC827 cells. F. The newly established PC9ER and HCC827ER were analyzed for lncRNAs expression patterns, and the primarily results were showed in Heatmap. G. The differentially expressed LncRNAs between Erlotinib resistant cells and sensitive cells were categorized by using GO analysis, and TUSC7 was supposed to be inhibited in PC9ER and HCC827ER cells. Real-time PCR further confirmed the suppression of TUSC7 in PC9ER (H) and HCC827ER (I) cells.
Additional file 3: Supplemental Figure 3.
Establishing the Erlotinib resistant cells. A. Osimertinib was dissolved in dimethyl sulfoxide (DMSO), and a total of 1 × 106 cells/ml of H1975 cells were seeded in a 6-wells plate and incubated in RPMI-1640 medium containing Osimertinib. The initial concentrations of Osimertinib were started with a concentration equal to the half-maximal inhibitory concentration of H1975 cells. After a cycle of Osimertinib treatment, only a small percentage of cells remained. Once cells had resumed normal growth and returned to 80% confluence under the light microscope, the next cycle began. The drug concentration was gradually increased for the next cycle until cells could survive with 10 μM Osimertinib. After 6 months, the H1975OR cells were successfully established and were then harvested for RNAs analysis. B. The images of gene probes detection were exhibited for illustration.
Additional file 4: Supplemental Figure 4.
TUSC7 re-sensitized the resistant PC9ER cells and HCC827ER cells through Notch signaling inhibition. A-B. All lung cancer cells responded to Notch signaling inhibitors greatly, and the much-lowered concentration of Notch signaling inhibitor sensitized both PC9ER and HCC827ER cells to Erlotinib treatment. TUSC7 stimulated the suppressive functions of Erlotinib in both PC9ER (C) and HCC827ER cells (D). E. Differences of the proliferation inhibition ratios between TUSC7 alone and the combination of TUSC7 and FLI-06 was insignificant. F-G. Supplemented images.
Additional file 5: Supplemental Figure 5.
Blotting results to prove that TUSC7 functioned through Notch signaling inhibition in resistant PC9ER cells and HCC827ER cells. The combined TUSC7 and Notch signaling inhibitor decreased the stem cells associated signatures, and decreased the EMT markers in PC9ER cells (A) and HCC827ER cells (B), but the inhibitory result was similar to that of using either TUSC7 or Notch signaling inhibitor alone.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, K., Peng, ZY., Gao, S. et al. M6A associated TSUC7 inhibition contributed to Erlotinib resistance in lung adenocarcinoma through a notch signaling activation dependent way. J Exp Clin Cancer Res 40, 325 (2021). https://doi.org/10.1186/s13046-021-02137-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13046-021-02137-9