EWS-FLI1 regulates and cooperates with core regulatory circuitry in Ewing sarcoma.

Core regulatory circuitry (CRC)-dependent transcriptional network is critical for developmental tumors in children and adolescents carrying few gene mutations. However, whether and how CRC contributes to transcription regulation in Ewing sarcoma is unknown. Here, we identify and functionally validate a CRC 'trio' constituted by three transcription factors (TFs): KLF15, TCF4 and NKX2-2, in Ewing sarcoma cells. Epigenomic analyses demonstrate that EWS-FLI1, the primary fusion driver for this cancer, directly establishes super-enhancers of each of these three TFs to activate their transcription. In turn, KLF15, TCF4 and NKX2-2 co-bind to their own and each other's super-enhancers and promoters, forming an inter-connected auto-regulatory loop. Functionally, CRC factors contribute significantly to cell proliferation of Ewing sarcoma both in vitro and in vivo. Mechanistically, CRC factors exhibit prominent capacity of co-regulating the epigenome in cooperation with EWS-FLI1, occupying 77.2% of promoters and 55.6% of enhancers genome-wide. Downstream, CRC TFs coordinately regulate gene expression networks in Ewing sarcoma, controlling important signaling pathways for cancer, such as lipid metabolism pathway, PI3K/AKT and MAPK signaling pathways. Together, molecular characterization of the oncogenic CRC model advances our understanding of the biology of Ewing sarcoma. Moreover, CRC-downstream genes and signaling pathways may contain potential therapeutic targets for this malignancy.

dbInDel: a database of enhancer-associated insertion and deletion variants by analysis of H3K27ac ChIP-Seq.

SUMMARY: Cancer hallmarks rely on its specific transcriptional programs, which are dysregulated by multiple mechanisms, including genomic aberrations in the DNA regulatory regions. Genome-wide association studies have shown many variants are found within putative enhancer elements. To provide insights into the regulatory role of enhancer-associated non-coding variants in cancer epigenome, and to facilitate the identification of functional non-coding mutations, we present dbInDel, a database where we have comprehensively analyzed enhancer-associated insertion and deletion variants for both human and murine samples using ChIP-Seq data. Moreover, we provide the identification and visualization of upstream TF binding motifs in InDel-containing enhancers. Downstream target genes are also predicted and analyzed in the context of cancer biology. The dbInDel database promotes the investigation of functional contributions of non-coding variants in cancer epigenome. AVAILABILITY AND IMPLEMENTATION: The database, dbInDel, can be accessed from http://enhancer-indel.cam-su.org/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Super-enhancer-associated MEIS1 promotes transcriptional dysregulation in Ewing sarcoma in co-operation with EWS-FLI1.

As the second most common malignant bone tumor in children and adolescents, Ewing sarcoma is initiated and exacerbated by a chimeric oncoprotein, most commonly, EWS-FLI1. In this study, we apply epigenomic analysis to characterize the transcription dysregulation in this cancer, focusing on the investigation of super-enhancer and its associated transcriptional regulatory mechanisms. We demonstrate that super-enhancer-associated transcripts are significantly enriched in EWS-FLI1 target genes, contribute to the aberrant transcriptional network of the disease, and mediate the exceptional sensitivity of Ewing sarcoma to transcriptional inhibition. Through integrative analysis, we identify MEIS1 as a super-enhancer-driven oncogene, which co-operates with EWS-FLI1 in transcriptional regulation, and plays a key pro-survival role in Ewing sarcoma. Moreover, APCDD1, another super-enhancer-associated gene, acting as a downstream target of both MEIS1 and EWS-FLI1, is also characterized as a novel tumor-promoting factor in this malignancy. These data delineate super-enhancer-mediated transcriptional deregulation in Ewing sarcoma, and uncover numerous candidate oncogenes which can be exploited for further understanding of the molecular pathogenesis for this disease.

Targetable BET proteins- and E2F1-dependent transcriptional program maintains the malignancy of glioblastoma.

Competitive BET bromodomain inhibitors (BBIs) targeting BET proteins (BRD2, BRD3, BRD4, and BRDT) show promising preclinical activities against brain cancers. However, the BET protein-dependent glioblastoma (GBM)-promoting transcriptional network remains elusive. Here, with mechanistic exploration of a next-generation chemical degrader of BET proteins (dBET6), we reveal a profound and consistent impact of BET proteins on E2F1- dependent transcriptional program in both differentiated GBM cells and brain tumor-initiating cells. dBET6 treatment drastically reduces BET protein genomic occupancy, RNA-Pol2 activity, and permissive chromatin marks. Subsequently, dBET6 represses the proliferation, self-renewal, and tumorigenic ability of GBM cells. Moreover, dBET6-induced degradation of BET proteins exerts superior antiproliferation effects compared to conventional BBIs and overcomes both intrinsic and acquired resistance to BBIs in GBM cells. Our study reveals crucial functions of BET proteins and provides the rationale and therapeutic merits of targeted degradation of BET proteins in GBM.

Master transcription factors form interconnected circuitry and orchestrate transcriptional networks in oesophageal adenocarcinoma.

OBJECTIVE: While oesophageal squamous cell carcinoma remains infrequent in Western populations, the incidence of oesophageal adenocarcinoma (EAC) has increased sixfold to eightfold over the past four decades. We aimed to characterise oesophageal cancer-specific and subtypes-specific gene regulation patterns and their upstream transcription factors (TFs). DESIGN: To identify regulatory elements, we profiled fresh-frozen oesophageal normal samples, tumours and cell lines with chromatin immunoprecipitation sequencing (ChIP-Seq). Mathematical modelling was performed to establish (super)-enhancers landscapes and interconnected transcriptional circuitry formed by master TFs. Coregulation and cooperation between master TFs were investigated by ChIP-Seq, circularised chromosome conformation capture sequencing and luciferase assay. Biological functions of candidate factors were evaluated both in vitro and in vivo. RESULTS: We found widespread and pervasive alterations of the (super)-enhancer reservoir in both subtypes of oesophageal cancer, leading to transcriptional activation of a myriad of novel oncogenes and signalling pathways, some of which may be exploited pharmacologically (eg, leukemia inhibitory factor (LIF) pathway). Focusing on EAC, we bioinformatically reconstructed and functionally validated an interconnected circuitry formed by four master TFs-ELF3, KLF5, GATA6 and EHF-which promoted each other's expression by interacting with each super-enhancer. Downstream, these master TFs occupied almost all EAC super-enhancers and cooperatively orchestrated EAC transcriptome. Each TF within the transcriptional circuitry was highly and specifically expressed in EAC and functionally promoted EAC cell proliferation and survival. CONCLUSIONS: By establishing cancer-specific and subtype-specific features of the EAC epigenome, our findings promise to transform understanding of the transcriptional dysregulation and addiction of EAC, while providing molecular clues to develop novel therapeutic modalities against this malignancy.

dbCoRC: a database of core transcriptional regulatory circuitries modeled by H3K27ac ChIP-seq signals.

Core transcription regulatory circuitry (CRC) is comprised of a small group of self-regulated transcription factors (TFs) and their interconnected regulatory loops. Studies from embryonic stem cells and other cellular models have revealed the elementary roles of CRCs in transcriptional control of cell identity and cellular fate. Systematic identification and subsequent archiving of CRCs across diverse cell types and tissues are needed to explore both cell/tissue type-specific and disease-associated transcriptional networks. Here, we present a comprehensive and interactive database (dbCoRC, http://dbcorc.cam-su.org) of CRC models which are computationally inferred from mapping of super-enhancer and prediction of TF binding sites. The current version of dbCoRC contains CRC models for 188 human and 50 murine cell lines/tissue samples. In companion with CRC models, this database also provides: (i) super enhancer, typical enhancer, and H3K27ac landscape for individual samples, (ii) putative binding sites of each core TF across the super-enhancer regions within CRC and (iii) expression of each core TF in normal or cancer cells/tissues. The dbCoRC will serve as a valuable resource for the scientific community to explore transcriptional control and regulatory circuitries in biological processes related to, but not limited to lineage specification, tissue homeostasis and tumorigenesis.

Potent Activity of the Bromodomain Inhibitor OTX015 in Multiple Myeloma.

Several studies demonstrate that the bromodomain inhibitor OTX015 has an antitumor activity in cancers. However, translation of these data to molecules suitable for clinical development has yet to be accomplished in multiple myeloma (MM). Here, we identified genes and biologic processes that substantiated the antimyeloma activity of OTX015 with global transcriptomics. OTX015 exerted a strong antiproliferative effect and induced cell cycle arrest in vitro. Gene expression profiling uncovered that OTX015 targeted NF-κB, EGFR, cell cycle regulation, and the cancer proliferation signaling pathway. Gene expression signatures displaying various levels of sensitivity to OTX015 were also identified. The data also showed that oral administration of OTX015 displayed significant antitumor activity in the mice model of disseminated human myeloma. In addition, our study provided the first evidence and rationale that OTX015 could promote osteoblast differentiation of mesenchymal stem cells (MSCs) and inhibited osteoclast formation and resorption in vivo experiments. Herein our results expanded the understanding of the mechanism for BET inhibitors OTX015 in MM. Our study provided an impressive basis for the clinical application of the novel antimyeloma agent OTX015 and uncovered signaling pathways that may play key roles in myeloma cell proliferation.

Upregulation of lncRNA MEG3 Promotes Osteogenic Differentiation of Mesenchymal Stem Cells From Multiple Myeloma Patients By Targeting BMP4 Transcription.

Multiple myeloma (MM) is characterized by the impaired osteogenic differentiation of mesenchymal stromal cells (MSCs). However, the underlying molecular mechanisms are still poorly understood. Long noncoding RNAs (lncRNAs) are emerging as important regulatory molecules in tumor-suppressor and oncogenic pathways. Here we showed that MSCs from MM expressed less lncRNA MEG3 relative to those from normal donors during osteogenic differentiation. To evaluate the effect of MEG3 on osteogenesis, bone marrow MSCs with enhanced or reduced MEG3 were prepared. We observed that MEG3 knockdown significantly reduced the expression of key osteogenic markers, including Runt-related transcription factor 2, osterix, and osteocalcin, while overexpression of MEG3 enhanced their expression. Additionally, MEG3 knockdown decreased BMP4 transcription. Here we showed that MEG3 was critical for SOX2 transcriptional repression of the BMP4. MEG3, which is located near the BMP4 gene, could dissociate the transcription factor SOX2 from the BMP4 promoter. A stable complex containing the MEG3, SOX2, and the SOX2 consensus site of BMP4 suggested that MEG3 activated transcriptional activity by directly influencing SOX2 activity. By using assays such as luciferase, chromatin immunoprecipitation, and RNA immunoprecipitation, we showed that MEG3 had a critical function in a mechanism of promoter-specific transcriptional activation. These results suggested that MEG3 played an essential role in osteogenic differentiation in bone marrow MSCs, partly by activating BMP4 transcription. Our data provided novel evidence for the biological and clinical significance of lncRNA MEG3 expression as a potential biomarker for identifying patients with MM and as a potential therapeutic target in MM.

Dual roles of FBXL3 in the mammalian circadian feedback loops are important for period determination and robustness of the clock.

The mammalian circadian clock is composed of interlocking feedback loops. Cryptochrome is a central component in the core negative feedback loop, whereas Rev-Erbα, a member of the nuclear receptor family, is an essential component of the interlocking loop. To understand the roles of different clock genes, we conducted a genetic interaction screen by generating single- and double-mutant mice. We found that the deletion of Rev-erbα in F-box/leucine rich-repeat protein (Fbxl3)-deficient mice rescued its long-circadian period phenotype, and our results further revealed that FBXL3 regulates Rev-Erb/retinoic acid receptor-related orphan receptor-binding element (RRE)-mediated transcription by inactivating the Rev-Erbα:histone deacetylase 3 corepressor complex. By analyzing the Fbxl3 and Cryptochrome 1 double-mutant mice, we found that FBXL3 also regulates the amplitudes of E-box-driven gene expression. These two separate roles of FBXL3 in circadian feedback loops provide a mechanism that contributes to the period determination and robustness of the clock.

Transcriptome analysis of BmN cells following over-expression of BmSTAT.

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway are involved in immune response, cell proliferation, differentiation, cell migration and apoptosis. In order to better understand the role of the JAK/STAT pathway in insects we chose Bombyx mori as an experimental model system. Over-expression of BmSTAT in a BmN cell line increased the number of cells in the G2 phase of the cell cycle. Genome-wide transcriptome analysis was performed to identify genes that were differentially expressed following BmSTAT overexpression. Transcriptome data showed that 10,853 and 10,129 expressed genes were obtained from the normal BmN cells and transformed cells, respectively. A total of 800 differentially expressed genes (DEGs) were detected, of which 787 were up-regulated and 13 were down-regulated with T test. In case of FC-test, 252 DEGs were detected, and 123 were expressed in the transformed cells and remaining were in the normal cells. Gene ontology (GO) annotation predicted a functional role for DEGs in catalytic activity, binding, transport, biological regulation, cellular and metabolic processes and pigmentation, while Kyoto encyclopedia of genes and genomes (KEGG) analysis revealed the affected genes to be involved in a multitude of cell signaling pathways. Our findings implicate JAK/STAT signaling in regulating the cell cycle in Bombyx mori, probably in combination with other pathways. These findings justify further investigation into the functional role of the BmSTAT gene.

The gene expression profile of resistant and susceptible Bombyx mori strains reveals cypovirus-associated variations in host gene transcript levels.

High-throughput paired-end RNA sequencing (RNA-Seq) was performed to investigate the gene expression profile of a susceptible Bombyx mori strain, Lan5, and a resistant B. mori strain, Ou17, which were both orally infected with B. mori cypovirus (BmCPV) in the midgut. There were 330 and 218 up-regulated genes, while there were 147 and 260 down-regulated genes in the Lan5 and Ou17 strains, respectively. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment for differentially expressed genes (DEGs) were carried out. Moreover, gene interaction network (STRING) analyses were performed to analyze the relationships among the shared DEGs. Some of these genes were related and formed a large network, in which the genes for B. mori cuticular protein RR-2 motif 123 (BmCPR123) and the gene for B. mori DNA replication licensing factor Mcm2-like (BmMCM2) were key genes among the common up-regulated DEGs, whereas the gene for B. mori heat shock protein 20.1 (Bmhsp20.1) was the central gene among the shared down-regulated DEGs between Lan5 vs Lan5-CPV and Ou17 vs Ou17-CPV. These findings established a comprehensive database of genes that are differentially expressed in response to BmCPV infection between silkworm strains that differed in resistance to BmCPV and implied that these DEGs might be involved in B. mori immune responses against BmCPV infection.

In Situ Spectroscopic Elucidation of the Electrochemical Potential Drop at Polyelectrolytes/Au Interfaces.

Polyelectrolytes have been widely applied in electrochemical devices. Understanding the polyelectrolyte/electrode interfaces is pivotal for polyelectrolyte-based applications. Here, we measured the electrochemical potential drop and the local activity of the mobile ion of H+ or OH- at the polyelectrolytes/Au interfaces by in situ electrochemical surface-enhanced Raman spectroscopy and voltammetry in three-electrode cells. We found that the potential dependences of the electrochemical potential drop in polyelectrolytes were smaller than that in conventional electrolyte solutions. The interfacial activity of H+ or OH- was much lower than that of bulk polyelectrolytes. The potential-dependent molecular dynamics simulations showed that the mobility of ionomers of polyelectrolytes in an electrostatic field was limited by a polymer matrix. These results suggested a characteristically thicker compact layer in the electrical double layer of a polyelectrolyte/electrode interface due to the accumulation of mobile H+ or OH- with a thicker hydration layer and immobile ionomers.

The deubiquitinase USP2a promotes tumor immunosuppression by stabilizing immune checkpoint B7-H4 in lung adenocarcinoma harboring EGFR-activating mutants.

B7-H4 is an immune checkpoint crucial for inhibiting CD8+ T-cell activity. A clinical trial is underway to investigate B7-H4 as a potential immunotherapeutic agent. However, the regulatory mechanism of B7-H4 degradation via the ubiquitin-proteasome pathway (UPP) remains poorly understood. In this study, we discovered that proteasome inhibitors effectively increased B7-H4 expression, while EGFR-activating mutants promoted B7-H4 expression through the UPP. We screened B7-H4 binding proteins by co-immunoprecipitation and mass spectrometry and found that USP2a acted as a deubiquitinase of B7-H4 by removing K48- and K63-linked ubiquitin chains from B7-H4, leading to a reduction in B7-H4 degradation. EGFR mutants enhanced B7-H4 stability by upregulating USP2a expression. We further investigated the role of USP2a in tumor growth in vivo. Depletion of USP2a in L858R/LLC cells inhibited tumor cell proliferation, consequently suppressing tumor growth in immune-deficient nude mice by destabilizing downstream molecules such as Cyclin D1. In an immune-competent C57BL/6 mouse tumor model, USP2a abrogation facilitated infiltration of CD95+CD8+ effector T cells and hindered infiltration of Tim-3+CD8+ and LAG-3+CD8+ exhausted T cells by destabilizing B7-H4. Clinical lung adenocarcinoma samples showed a significant correlation between B7-H4 abundance and USP2a expression, indicating the contribution of the EGFR/USP2a/B7-H4 axis to tumor immunosuppression. In summary, this study elucidates the dual effects of USP2a in tumor growth by stabilizing Cyclin D1, promoting tumor cell proliferation, and stabilizing B7-H4, contributing to tumor immunosuppression. Therefore, USP2a represents a potential target for tumor therapy.

Advances in sequencing-based studies of microDNA and ecDNA: Databases, identification methods, and integration with single-cell analysis.

Extrachromosomal circular DNA (eccDNA) is a class of circular DNA molecules that originate from genomic DNA but are separate from chromosomes. They are common in various organisms, with sizes ranging from a few hundred to millions of base pairs. A special type of large extrachromosomal DNA (ecDNA) is prevalent in cancer cells. Research on ecDNA has significantly contributed to our comprehension of cancer development, progression, evolution, and drug resistance. The use of next-generation (NGS) and third-generation sequencing (TGS) techniques to identify eccDNAs throughout the genome has become a trend in current research. Here, we briefly review current advances in the biological mechanisms and applications of two distinct types of eccDNAs: microDNA and ecDNA. In addition to presenting available identification tools based on sequencing data, we summarize the most recent efforts to integrate ecDNA with single-cell analysis and put forth suggestions to promote the process.

Circlehunter: a tool to identify extrachromosomal circular DNA from ATAC-Seq data.

In cancer, extrachromosomal circular DNA (ecDNA), or megabase-pair amplified circular DNA, plays an essential role in intercellular heterogeneity and tumor cell revolution because of its non-Mendelian inheritance. We developed circlehunter ( https://github.com/suda-huanglab/circlehunter ), a tool for identifying ecDNA from ATAC-Seq data using the enhanced chromatin accessibility of ecDNA. Using simulated data, we showed that circlehunter has an F1 score of 0.93 at 30× local depth and read lengths as short as 35 bp. Based on 1312 ecDNAs predicted from 94 publicly available datasets of ATAC-Seq assays, we found 37 oncogenes contained in these ecDNAs with amplification characteristics. In small cell lung cancer cell lines, ecDNA containing MYC leads to amplification of MYC and cis-regulates the expression of NEUROD1, resulting in an expression pattern consistent with the NEUROD1 high expression subtype and sensitive to Aurora kinase inhibitors. This showcases that circlehunter could serve as a valuable pipeline for the investigation of tumorigenesis.

Identification of a novel eighteen-gene signature of recurrent metastasis neuroblastoma.

Neuroblastoma is the most common malignant tumor in childhood, and metastases occur in more than 30% patients. Recurrent metastasis is the main cause of poor prognosis and high mortality in neuroblastoma. In this regard, there is still a lack of sufficient biomarkers and effective therapies. Therefore, we performed a multi-omics analysis of neuroblastoma patients from Therapeutically Applicable Research To Generate Effective Treatments (TARGET). With clinical relapse site information, tumor samples derived from the primary site were divided into recurrent metastasis and primary tumor groups. The initial gene signature was obtained by comparing RNA-Seq and copy number variation differences. Survival data was used to further filter prognosis-related genes. This 18-gene signature consists of three clusters: tumor suppression, cell proliferation, and immunity. A super enhancer is involved in the enhanced expression of NCAPG in cluster2 together with IRF3. Based on the gene signature expression in primary neuroblastoma, it is possible to predict tumor metastasis before it occurs. According to the anticancer drug dataset of Genomics of Drug Sensitivity in Cancer (GDSC), vinorelbine and docetaxel were predicted to have high sensitivity against recurrent metastatic neuroblastoma. In conclusion, our study offers a novel metastasis biomarker and helps understand the mechanisms of tumor recurrent metastasis. KEY MESSAGES: We identified a novel eighteen-gene signature of recurrent metastasis neuroblastoma and build risk and classification models. We dissected the regulatory role of NCAPG in signatures. We found immune exhaustion and immunosuppression in recurrent metastasis neuroblastoma. Vinorelbine and docetaxel were predicted to have high sensitivity against recurrent metastatic neuroblastoma.

Rhythmic Component Analysis Tool (RCAT): A Precise, Efficient and User-Friendly Tool for Circadian Clock Genes Analysis.

High-throughput next-generation sequencing (NGS) technologies and real-time circadian dynamics reporting systems produce large amounts of experimental data on RNA and protein levels in the field of circadian rhythm and therefore require statistical knowledge and computational skills for quantitative analysis. Although there are many software applications that can process these data, they are often difficult to use and computationally inefficient. Hence, a convenient, user-friendly tool that can accurately acquire rhythmic components (period, amplitude, and phase) of circadian clock genes is necessary. Here, we develop a new analysis tool named rhythmic component analysis tool (RCAT), which has an easily understood interface featuring a one-button operation, that presents all results as tables and images and automatically saves them as CSV files. We use the relative amplitude error (RAE), widely-adopted criteria on the circadian research field to estimate the quality of results. To illustrate the analytical ability of the RCAT under different situations, we generate four groups of time-series data by CircaInSilico (a web server for generating synthetic genome biology data to benchmark statistical methods for studying biological rhythms) with different collection intervals and amplitude ranges and use RCAT to analyze them. To demonstrate the effectiveness of RCAT, we analyze two sets of case studies with time-series data: one set uses microarray and RNA-Seq data from the gene expression omnibus (GEO) repository to identify core clock genes (CCGs) with significant periodicity in the liver, and the other set uses real-time fluorescence reporting data collected by Lumicycle® (a commonly-used luminometer) to calculate the precise period, amplitude and phase. In these examples, most cycling samples are successfully detected by the RCAT within a short collection time, and accurate rhythmic components are also successfully computed. These results indicate that RCAT improves flexibility and convenience in periodic oscillation data analysis. RCAT, is freely available at: https://github.com/lzbbest/Rhythmic-Component-Analysis-Tool/releases . It, as a cross-platform software, can be run not only on Linux, but also on Win10, Win8 and Win7.

B7-H4 is increased in lung adenocarcinoma harboring EGFR-activating mutations and contributes to immunosuppression.

PD-1/PD-L1 inhibitors have shown clinical benefit in lung adenocarcinoma (LUAD). However, the immunotherapy strategy is less effective in patients with EGFR-activating mutations (EGFR MT). Studies showed that besides low expression of PD-L1, the absence of TILs and distinct expression profile of immune checkpoint molecules might be associated with low response of the patient subset. In this study, we first compared CD8A, GZMB and PRF1 mRNA levels in different LUAD subtypes harboring different driver mutations by dataset analyses and investigated the association between 15 well-defined B7-CD28 family members and driver mutations. The results showed that the decreases in the density and function of CD8+ TILs, CD274 (PD-L1 gene), and CD86 and increases in VTCN1 (B7-H4 gene) and HHLA2 were associated with LUAD with EGFR-activating mutations. Immunohistochemical staining studies further supported that PD-L1 was downregulated and B7-H4 was upregulated in the subtype. Furthermore, PD-L1 expression was positively associated with levels of CD8A and granzyme B, while B7-H4 expression was negatively associated with granzyme B levels. In lung cancer cell lines, EGFR-activating mutations effectively upregulated B7-H4 and downregulated PD-L1. MEK/ERK-pathway activation upregulated B7-H4, and PI3K/Akt activation upregulated PD-L1. EGFR 19Del mutation was associated with inhibition of CD8+ T-cell function, while knocking down B7-H4 could reverse the inhibition, and further showed tumor-growth inhibition and longer survival in vivo. Taken together, this study shed light on that B7-H4 might be an alternative immune-checkpoint molecule and a potential therapeutic target for LUAD with EGFR MT.

Isovalerylspiramycin I suppresses non-small cell lung carcinoma growth through ROS-mediated inhibition of PI3K/AKT signaling pathway.

Novel drugs are required for non-small cell lung cancer (NSCLC) treatment urgently. Repurposing old drugs as new treatments is a practicable approach with time and cost savings. Some studies have shown that carrimycin, a Chinese Food and Drug Administration (CFDA)-approved macrolide antibiotic, possesses potent anti-tumor effects against oral squamous cell carcinoma. However, its detailed component and underlying mechanisms in anti-NSCLC remain unknown. In our study, isovalerylspiramycin I (ISP-I) was isolated from carrimycin and demonstrated a remarkable anti-NSCLC efficacy in vitro and in vivo with a favorable safety profile. It has been proven that in NSCLC cell lines H460 and A549, ISP-I could induce G2/M arrest and apoptosis, which was mainly attributed to ROS accumulation and subsequently PI3K/AKT signaling pathway inhibition. Numerous downstream genes including mTOR and FOXOs were also changed correspondingly. An observation of NAC-induced reverse effect on ISP-I-leading cell death and PI3K/AKT pathway inhibition, emphasized the necessity of ROS signaling in this event. Moreover, we identified ROS accumulation and PI3K/AKT pathway inhibition in tumor xenograft models in vivo as well. Taken together, our study firstly reveals that ISP-I is a novel ROS inducer and may act as a promising candidate with multi-target and low biological toxicity for anti-NSCLC treatment.