RNA-seq Explores the Mechanism of Oxygen-Boosted Sonodynamic Therapy Based on All-in-One Nanobubbles to Enhance Ferroptosis for the Treatment of HCC

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

Download PDF

Research

RNA-seq Explores the Mechanism of Oxygen-Boosted Sonodynamic Therapy Based on All-in-One Nanobubbles to Enhance Ferroptosis for the Treatment of HCC

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

This work is licensed under a CC BY 4.0 License

Journal Publication

published 01 Jan, 2022

Read the published version in International Journal of Nanomedicine →

Version 1

posted

You are reading this latest preprint version

Background

The combination of sonodynamic therapy and oxygenation strategy is widely used in cancer treatment. However, due to the complexity, heterogeneity and irreversible hypoxic environment produced by Hepatocellular carcinoma (HCC) tissues, oxygen-enhancing sonodynamic therapy (SDT) have failed to achieve the desired results. With the emergence of ferroptosis with reactive oxygen species (ROS) cytotoxicity, this novel cell death method has attracted widespread attention.

Methods

In this study, nanobubbles (NBs) were connected with the sonosensitizer Indocyanine green (ICG) to construct a 2-in-1 nanoplatform loaded with RAS-selective lethal (RSL3, ferroptosis promoter) (RSL3@O₂-ICG NBs), combined with oxygen-enhanced SDT and potent ferroptosis. In addition, nanobubbles (NBs) combined with low-frequency ultrasound (LFUS) are called ultrasound-targeted nanobubble destruction (UTND) to ensure specific drug release and improve safety.

Results

MDA/GSH and other related experimental results show that RSL3@O₂-ICG NBs can enhance SDT and ferroptosis. Through RNA sequencing (RNA-seq), the differential expression of LncRNA and mRNA before and after synergistic treatment was identified, and then GO and KEGG pathways were used to enrich and analyze target genes and pathways related ferroptosis sensitivity. Then we searched for the expression of differentially expressed genes in the TCGA Hepatocellular carcinoma cohort. At the same time, we evaluated the proportion of immune cell infiltration and the identification of co-expression network modules and related prognostic analysis. We found that it was significantly related to the tumor microenvironment of hepatocellular carcinoma. The prognostic risk genes "SLC37A2" and "ITGB7" may represent new hepatocellular carcinoma ferroptosis-inducing markers and have guiding significance for treating hepatocellular carcinoma.

Conclusions

The therapeutic effect of the in vitro synergistic treatment has been proven to be significant, revealing the prospect of 2-in-1 nanobubbles combined with SDT and ferroptosis in treating HCC.

Cancer Biology

Oncology

Ferroptosis

Sonodynamic therapy

HCC

RSL3

Synergistic therapy.

Hepatocellular carcinoma, also known as primary liver cancer, as the fourth cancer incidence, has been the fastest-rising cause of cancer-related deaths worldwide (1). There are approximately 840,000 new cases of HCC each year, and at least 780,000 people die from HCC (2). The main treatment options for liver cancer are surgical resection, liver transplantation, local radiotherapy, chemotherapy (3). However, due to the invisibility of the onset of HCC and the lack of specific early markers, most patients are usually diagnosed as advanced HCC (4), the above treatment methods cannot get a thorough treatment, and sorafenib has been the mainstream targeted therapy for advanced HCC. Many nanomedicines have been invented for HCC and rodent hepatoma model treatment(5–7). In addition, some nanocapsules could also promote chemotherapy and induce an immune response, which synergistically performed antitumor efficacy(8–10).

However, drug resistance gradually develops in the later stages of treatment(11, 12). According to reports, sorafenib is mainly causing cell ferroptosis to achieve anti-tumor, and the role of apoptosis in sorafenib resistance in HCC has been demonstrated (13). Induce ferroptosis and increase the sensitivity of ferroptosis, avoid apoptosis, effectively improve the treatment of hepatocellular carcinoma, and reduce the occurrence of drug resistance (14).

Ferroptosis was first proposed by Dixon as a novel cell death in 2012 (15, 16). Scientists found a new way of cell death related to increased ROS levels in cells, which can be prevented by iron chelation or inhibition of iron uptake, and they suggested using the name ferroptosis to describe this phenotype. Unlike other apoptotic and non-apoptotic death forms, iron and lipid peroxides are the two main features of ferroptosis. The redox balance in normal cells mainly depends on glutathione peroxidase 4 (GPX4), a repair enzyme for cell lipid oxidation. Therefore, directly or indirectly inhibiting the activity of GPX4, leading to the accumulation of LPO, can lead to subsequent ferroptosis. (17) RSL3, an enzyme inhibitor of GPX4, a small molecule that causes ferroptosis, was selectively lethal to oncogenic RAS mutant cell lines. The RAS family of small GTPases (HRAS, NRAS, and KRAS) are commonly mutated in HCC. Finding compounds that are selectively lethal to RAS mutant tumor cells is also a high priority.

Ultrasound-triggered SDT has attracted much attention as a non-invasive way to treat tumors (18, 19). SDT combines a sonosensitizer to generate ROS to kill cancer cells. As enough ROS are produced, a series of biological events can be activated, including DNA fragmentation, cytoskeleton contraction, and chromatin condensation, leading to apoptosis (20).In summary, the possible potential mechanisms of SDT include ultrasonic cavitation, generation of free radicals, apoptosis, or any combination of these (21, 22). SDT with LFUS can avoid photodynamic therapy (PDT) limitations with low tissue penetration depth (23). PDT could generate immunogenic cell death (ICD), accompanied by the release of high-mobility group box 1 protein (HMGB1) and adenosine triphosphate (ATP) and the exposure of calreticulin (CRT), sending the “eat me” signal and promoting the antigen presentation and maturation of dendritic cells (DCs)(24). Also, Photothermal therapy (PTT) is a developing field that uses the strategy of employing near-infrared (NIR) irradiation to trigger drug delivery(25). These emerging non-invasive techniques have been widely used in anti-tumor research(26–29). However, due to the intrinsic or obtained resistance of cancer cells on apoptosis, the effectiveness of SDT in tumor-induced induction is limited. Therefore, the other non-apoptotic cell death model has opened up new treatment pathways for eliminating cancer cells and reducing the resistance brought about apoptosis. ICG, A near-infrared fluorescent dye to diagnose and treat hepatocellular carcinoma, which is widely used in liver surgery and fluorescence imaging system (30–33). Many research shows that ICG can produce a large amount of ROS after being triggered by LFUS, which is a superior sonosensitizer (34–38). At the same time, ICG has high HCC specificity. However, the further application of ICG is limited for several reasons, such as its tendency to aggregate, non-specific binding to proteins, and rapid clearance from the body (39, 40). Based on this, scientists continue to explore various carrier materials to overcome these limitations and maintain the stable optical performance of ICG (41, 42).

Nanobubbles have been widely used as vehicles with LFUS called UTND compared to other nanocarriers. Under ultrasound stimulation, NBs may expand, compress, and destroy, generating shock waves to increase cell membranes' permeability further, thereby promoting the delivery of intracellular drugs and sonosensitizer, due to the destruction of NBs, local temperature increases to promote ROS generation. Several studies have shown that lipid chitosan core NBs can successfully reverse tumor hypoxia by down-regulating the expression of HIF-1α (43, 44). Thus, NBs are a good candidate for research.

Next-generation high-throughput RNA sequencing (RNA-Seq) is the first sequencing-based technology to perform high-throughput and quantitative measurements of the entire transcriptome. RNA-Seq is a more accurate method of presenting transcript and isotype levels than any other sequencing technology. In order to further explore the mechanism of combination therapy to increase the sensitivity of ferroptosis and find related gene changes, combining RNA-Seq technology in our experiments can obtain new insights into tumor treatment mechanisms.

Herein, RSL3@O₂-ICG NBs, which have the ability to oxygen-boosted SDT and increase sensitization ferroptosis, was used in the combined treatment of HCC. Then RNA-Seq further explores its mechanism and finds related gene changes, which provide a promising strategy and therapeutic target for HCC therapy.

2.1 chemicals

1,2-Distearoyl-sn-glycerol-3-phosphocholine(DSPC);1,2-stearoyl-sn-glycerol-3-phosphoethanolamine-N-(methoxy[polyethylene glycol]-2000) (DSPE-PEG-2000) were got from Avanti Polar Lipids (Alabaster, AL). 20,70-Dichlorofluorescin diacetate (DCFH-DA),MDA and GSH were acquired from Beyotime Biotechnology (Shanghai, China). ICG was received from Shanghai yuanye Bio-Technology Co, Ltd (Shanghai, China). Annexin V-FITC Apoptosis Detection Kit and Cell counter kit-8 (CCK-8) were obtained from BD Pharmingen. RSL3 was gained from Shanghai MCE Chemicals Technology (Shanghai, China).

2.2 Preparation of RSL3@O₂-ICG NBs

As reported in previous studies, we have prepared RSL3@O₂-ICG NBs using the thin-film hydration ultrasound method (45). First, mix DSPC and DSPE-PEG-2000 at a mass ratio of 9:1, and then add 10ul of 5 mM RSL3 in methylene chloride and methanol (2:1, vol/vol). After entirely dissolving and mixing, transfer the above solution to a beaker and evaporate to form a lipid film. The dried film was hydrated into a lipid suspension with 5ml of 100ug/ml ICG in PBS. Next, the suspension was extruded 20 times through a 200nm membrane through a micro-extruder (Avanti Polar Lipids, Alabaster, AL). Then, the squeezed suspension was transferred to a sealed vial, and the syringe was evacuated and replaced with pure oxygen. The mixture was then mechanically shaken in a dental mixer (YJT Medical Apparatus and Instruments, Shanghai, China) for 60 seconds and resuspended in 2 ml of PBS solution for storage at 4°C. All processes are carried out in a dark environment. Use the same method to prepare nanobubbles without RSL3.

2.3 Characterization

Observe the distribution and morphology of RSL3@O₂-ICG NBs under a scanning electron microscope (SEM. Hitachi, Japan), and observe the element enrichment. Dynamic light scattering (DLS) analyzer (ZEN3600, Malvern Instruments) can analyze its particle size distribution and zeta potential. An ultraviolet-visible spectrophotometer (Thermo Fisher) was used to obtain the absorption spectrum of the particles.

2.4 cell culture

HepG2 and Huh7 derived from the human hepatocellular carcinoma cell line were included in this experiment. All cells were cultured in DMEM medium containing 10% FBS in a humidified incubator at 37°C and 5% CO₂.

2.5 Cell viability of HepG2 and Huh7 cells

HepG2 and Huh7 cells were seeded into 96-well plates at a density of 8000 cells per well. After the application of different ICG concentrations (40, 80, 100, 150, 200, 250, 300µg/ml), they were exposed to LFUS for 30 seconds (1.0 W/cm²) and then assessed by CCK8 assay to kill the cells. Then different RSL3 concentrations (4, 6, 8, 10, 12, 15, 20 µM) and specific concentrations of ICG under the same experimental environment.

2.6 ROS detection

In order to study the ability of O₂-ICG-NB to enhance SDT, RSL3 was not carried in the nanobubbles in this experiment. HepG2 and Huh7 cells were respectively seeded in 6-well plates, cultured overnight to remove the medium, and the cells were treated with PBS, ICG, ICG-NBs, and O₂-ICG-NBs for another 4 hours. Subsequently, DCFH-DA was added to the cells and incubated at 37°C for 30 minutes, protected from light. After removing the medium, the cells were irradiated with ultrasound for 30 seconds (1.0 W/cm²) and washed twice with DMEM, then 1ml of ice PBS was added, and the cells were quickly observed under a fluorescence microscope. Flow cytometry is performed by digesting cells from a 6-well plate, and the operation is processed as before.

2.7 Analysis of GSH

After various treatments (PBS, RSL3, O₂-ICG-NBs + LFUS, RSL3@O₂-ICG NBs + LFUS, RSL3@O₂-ICG NBs + LFUS + Fer-1), the cells were washed 3 times with PBS and trypsinized Collect the cells and centrifuge at 1000 rpm at 4°C for 5 min twice. According to the weight of the cell pellet, resuspend each milligram in 3 µl of protein removal solution M. After vortexing, place it in liquid nitrogen and in a 37 ℃ water bath for two quick freezing and thawing times, and then place it at 4 ℃. Centrifuge for 10 minutes at 10,000 rpm for 5 minutes. The supernatant is used to determine the content of GSH in the sample. Prepare the appropriate amount of detection working solution and supernatant sample according to the instructions, mix them and incubate at room temperature for 5 minutes, and then add NADPH to the system to trigger the reaction. After 60 minutes, measure the TNB absorbance at 412 nm and calculate the GSH level according to the instructions.

2.8 MDA assay

Collect the protein in HepG2 and Huh7 through various treatments (PBS, RSL3, O₂-ICG-NBs + LFUS, RSL3@O₂-ICGNBs + LFUS, RSL3@O₂-ICG NBs + LFUS + Fer-1), and after lysis 40min After centrifugation, the supernatant was mixed with the TBA detection working solution, heated in a boiling water bath for 15 minutes, cooled to room temperature, 200µl was transferred to a 96-well plate, the absorbance was measured at 532nm, compare the standard curve to calculate the MDA content.

2.9 Measurement of Mitochondrial Membrane Potential

Mitochondrial membrane potential was detected by JC-1 (Beyotime, C2005) detection kit. After treatment in the 4 treatment groups (PBS, RSL3, O₂-ICG-NBs + LFUS, RSL3@O₂-ICG NBs + LFUS), HepG2 and Huh7 cells were incubated with JC-1 staining working solution for 20 minutes, and then JC-1 Wash twice with buffer solution. After adding ice PBS solution, the stained cells were observed under a fluorescence microscope.

2.10 Cell Cycle Analysis

After various treatments (PBS, RSL3, O₂-ICG-NBs + LFUS, RSL3@O₂-ICG NBs + LFUS), HepG2 and Huh7 cells were digested and collected with trypsin, and then fixed with 70% ice-cold ethanol at 4°C overnight. Subsequently, add 500µl of the prepared staining working solution (keygentec. Nanjing, China) for staining. The DNA content of each group of cells was determined by flow cytometry.

2.11. Whole transcriptome library construction, sequencing, and analysis

The total RNA of HepG2/Huh7 cells before and after RSL3@O₂-ICG NBs + LFUS treatment was taken, respectively, named A, B, C, D group, and quality control was performed. The preparation and deep sequencing of the whole transcriptome library were performed by Novogene Bioinformatics (Beijing, China).

2.11.1 Whole transcriptome library construction and sequencing

A total amount of 3 µg RNA per sample was used as input material for the RNA sample preparations. According to the manufacturer's instructions, the clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia). After cluster generation, the libraries were sequenced on an Illumina HiSeq 2500 platform, and 125 bp paired-end reads were generated. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE178573 (https://www.ncbi.nlm.nih.gov/geo/query /acc.cgi?acc = GSE14 0729).

2.11.2 RNA-Seq data analysis

Quality control and comparative analysis

The generation of sequencing data goes through multiple steps such as RNA extraction, library construction, and sequencing. These steps will produce some low-quality or invalid data. For example, there will be a deviation in the length of the library during the library construction stage, and a sequencing error will occur during the sequencing stage. Therefore, it is necessary to perform quality control on the raw data obtained from the machine to ensure the accuracy of subsequent analysis. The raw data obtained by sequencing contains a small number of reads with sequencing adapters or low sequencing quality. In order to ensure the quality and reliability of data analysis, it is necessary to filter the original data. The filtering content is as follows: (1) Remove the adapter (2) Remove the reads with N (N means that the base information cannot be determined) greater than 0.002; (3) When the number of low-quality bases contained in a single-ended read exceeds read When the length is 50%, you need to remove the paired reading. After raw data filtering, sequencing error rate check, and GC content distribution check, clean reads for subsequent analysis are obtained. Use Hisat2 to compare the valuable sequencing data (clean reads) to the reference genome.

Quantitative expression level of transcripts

Stringtie software can accurately and accurately splice genes and realize gene quantification. Therefore, we use Stringtie software to assemble the reads into genes and quantify gene expression based on the results of the alignment to the genome.

In this article, the gene expression level of RNA-seq uses FPKM (Fragments Per Kilobase of transcript sequence per Millions of base pairs sequenced) calculation. FPKM is the number of fragments from a gene/transcript per kilobase length per million fragments. It also corrects the influence of sequencing depth and gene length on the count of fragments. It is currently the most commonly used method for estimating gene expression levels.

2.11.3 Differential expression analysis

After the quantitative analysis is completed, the gene expression matrix of all samples is obtained, and the expression level of the gene level can be analyzed for the significance of the difference to find the functional genes related to the treatment group. In this work, we use edgeR software to analyze the significance of gene expression differences. Use p-value or corrected p-value (padj) to determine the significance level. Padj is the value obtained by using the BH method to perform multiple test corrections on the p-value when the false positive rate is high. The smaller the corrected P-value value is, the more significant it is. Use padj < 0.05 as the difference significance criterion (if padj is < 0.05 to screen for too few differences, use p-value < 0.05 for differential screening, and finally get a list of genes with significant differences for subsequent analysis.

2.11.4 GO and KEGG enrichment analysis

In organisms, different genes coordinate with each other to perform their biological functions. By exploring the pathways that are significantly enriched in functional genes, we will further explore the important functional genes involved in important biochemical metabolic pathways and signal transduction pathways among differentially expressed genes. This article uses GOseq software for GO enrichment analysis; KOBAS (2.0) for pathway enrichment analysis. GO is a comprehensive database describing gene functions, divided into molecular function, biological process, and cellular component. GO enrichment takes padj < 0.05 (default) or p-value < 0.05 as significant enrichment. In organisms, different genes coordinate with each other to perform their biological functions, and the most important biochemical metabolic pathways and signal transduction pathways involved in candidate target genes can be determined through Pathway significant enrichment. KEGG (Kyoto Encyclopedia of Genes and Genomes) is a comprehensive database that integrates genome, chemistry and system function information. Pathway significant enrichment analysis takes KEGG Pathway as the unit and applies a hypergeometric test to find pathways that are significantly enriched in candidate target genes compared with the background of the entire genome. Similarly, KEGG pathway enrichment takes p-value or padj < 0.05 as significant enrichment.

2.11.5 Expression and immune infiltration in TCGA.

In order to further explore whether the genomic changes caused by RSL3@O₂-ICG NBs + LFUS treatment will affect the immunity of hepatocellular carcinoma,we obtained the RNA-seq data of hepatocellular carcinoma from the TCGA database (https://portal.gdc.cancer.gov/). Including gene expression information from 373 samples of hepatocellular carcinoma patients༌we find some of the differential genes that overlap with this research༌and use overlapping differentially expressed genes to perform heat map cluster analysis on TCGA hepatocellular carcinoma patients. At the same time, we used Xcell database (https://xcell.ucsf.edu/) to evaluate the proportion of immune cell infiltration of all liver cancer samples from TCGA. Explore the difference in the proportion of immune cell infiltration between different clusters.

2.11.6. Prognostic analysis of TCGA cohort

In order to further confirm that the genomic changes caused by RSL3@O₂-ICG NBs + LFUS treatment can regulate the tumor microenvironment and further affect the survival status of hepatocellular carcinoma, we used the tumor microenvironment score obtained from the Xcell database in the previous step. We used the pearson method to screen out the results. Among the differentially expressed genes, the genes whose correlation with the tumor microenvironment score is more significant than 0.4 and the significance P-value is less than 0.05 will be analyzed later. 373 TCGA hepatocellular carcinoma patients were randomly divided into training set and test set according to the ratio of 7:3, and single-factor cox risk regression and multi-factor risk regression were used to screen the prognostic risk genes that were significantly related to the prognosis (p < 0.05). Calculate each sample's ferroptosis risk score by multiplying the risk gene's expression value by the sum of the risk coefficient of the risk regression and using the median value of the risk score to divide hepatocellular carcinoma samples into high-risk and low-risk groups. Use log-rank test to explore the prognostic difference between high and low risk groups.

2.12 Statistical analysis

All experiments were conducted in triplicate. The biological data are calculated as an average of three replicates ± standard deviation and presented in percentages. All sample values were related to the control (untreated cells) value, which was considered 100% and compared using Student's t-test by GraphPad Prism 8.0 software. The differences between control and samples values were considered to be significant if P < 0.05, highly significant if P < 0.01 and extremely significant if P < 0.001.

3.1 Preparation and characterization of RSL3@O₂-ICG NBs.

It can be seen from Fig. 1A, the prepared RSL3@O₂-ICG NBs had a green suspension compared with blank NBs and were well dispersed in the PBS solution. At the same time, the electron microscopy energy spectrum measurement was carried out, and it was observed that the representative Chlorine element of RSL3 and Sodium Sulfur elements of ICG were enriched on the phospholipid nanobubbles (Fig. 1B-F). The mean diameter of RSL3@O₂-ICG NBs measured by DLS was 280.8 ± 16.8 nm (Fig. 1G). At the same time, the zeta potential of the complex is − 3.86 ± 0.35 mV ( Table S1). As shown in Fig. 1H, the composite has absorption peaks at about 770nm and 240nm, indicating that ICG and RSL3 are successfully loaded into RSL3@O₂-ICG NBs, which is consistent with the previous SEM imaging results.

3.2 Detection of intracellular ROS of SDT.

CCK8 experiment verified the cell viability of different concentrations of ICG on Hepg2/Huh7 cells in Supplementary Figure S1A, B. We choose 100ug/ml of ICG used for the next experiment. The fluorescence microscope shows that HepG2/Huh7 cells presented the most amount of green fluorescence in the O₂-ICG-NBs + LFUS group (Fig. 2A). Flow cytometry also showed that the number of fluorescence intensities in the O₂-ICG-NBs + LFUS group was much higher, which was similarly displayed in fluorescence microscopy imaging. Therefore, the role of O₂-ICG-NBs in enhancing SDT has been proven. (Fig. 2B and Supplementary Figure S1C, D).

3.3 The enhancement effect of ferroptosis.

First of all, CCK8 experiment found that different concentrations of RSL3 combined with SDT, the HepG2/Huh7 cells viability decreased significantly compared to RSL3 treatment alone (Supplementary Figure S2A, B). Mitochondrial damage is a significant feature of iron death. We found that in the combined treatment group RSL3@O₂-ICG NBs + LFUS, the red fluorescence in HepG2/Huh7 cells decreased significantly, while the green fluorescence increased, indicating the decrease in mitochondrial membrane potential and damage (Fig. 3A). Then we conducted GSH and MDA experiments. GSH level is usually used to evaluate the intensity of ferroptosis (46, 47). As seen from Fig. 3C and Fig. 3E, Compared with the PBS treatment group, the GSH level of the combination treatment group was significantly reduced. In order to prove that the GSH level reduction in the combination treatment group was mainly caused by ferroptosis, Fer-1 (1µM), a ferroptosis inhibitor, the drugs were added to our combined treatment group, and it was found that GSH did not decrease significantly. At the same time, MDA is an advanced fat oxidation end product. The results in Fig. 3B and Fig. 3D show that the expression level increases throughout the process, and there is a negative correlation between GSH.

3.4 Influence on Cell Cycle Progression In Vitro.

After various treatments, more than half of HepG2/Huh7 cells stagnated in G0/G1 phase in the combination therapy group compared to PBS groups (Fig. 4). These data suggested that RSL3@O₂-ICG NBs with LFUS could inhibit cell cycle progression, thus affecting HepG2/ Huh7 cells growth.

3.5 Differential expression analysis based on RNA-seq.

Estimate the expression of mRNA and lncRNA by expression level (FPKM). Show the difference between differentially expressed mRNA and lncRNA through the volcano graph (P < 0.05). As a result, a total of 3391 mRNAs in group A and group B (AvsB) were differentially expressed significantly, of which 2307 were up-regulated (Fig. 5A). Similarly, compared with group B, 2281 lncRNAs showed significant differences in expression, including 1516 up-regulated genes and 765 down-regulated genes (Fig. 5B). A total of 2582 mRNAs expression differences between C group and D group (CvsD), of which 1335 were down-regulated (Fig. 5C). 2227 lncRNA expressions are significantly different, including 1328 up-regulated genes and 902 down-regulated genes (Fig. 5D). A total of 886 mRNAs and 570 lncRNAs overlapped between the differentially expressed genes between group A and group B and the differentially expressed genes between group C and D (Supplementary Figure S3A, B).

3.6 Differentially expressed analysis.

In the previous step of the differentially expressed gene we obtained, we found that A vs. B, C vs. D are also significantly enriched in cell functions such as regulation of protein phosphorylation, leukocyte chemotaxis, MAPK cascade, and response to alcohol, regulation of cell proliferation (Table 1). At the same time, the differentially expressed mRNAs of the two groups A vs. B, C vs. D of cell lines were enriched by the KEGG pathway, respectively. The bubble chart (Fig. 6A, B) shows the top 20 main enriched pathways. After comparison, we found that after RSL3@O₂-ICG NBs + LFUS treatment, the common enrichment pathway is ECM-receptor interaction. We predict the target genes of all differentially expressed lncRNAs (Table S2,3 ), that is, predict the target genes of lncRNAs through the co-location of lncRNAs and protein-coding genes. GO and KEGG enrichment analyses were performed on target genes to predict the function of lncRNA. Functional enrichment of the predicted target genes found that the target genes of AvsB and CvsD were also significantly enriched in protein phosphorylation, ubiquitin protein ligase binding, regulation of Wnt receptor signaling pathway, cell cycle process, autophagy, and other pathways (Table 2). The result of KEGG pathway enrichment shows that the differentially expressed LncRNA target genes are also significantly enriched in the Alcoholism pathway (Fig. 6C, D ).

Table 1

The representive results of GO Enrichment Analysis with differential mRNA
GO accession	Description	Term type	Ajusted P value (AvsB)	Ajusted P value (CvsD)
GO:0006468	protein phosphorylation	biological_process	0.00028462	0.031899
GO:0010562	positive regulation of phosphorus metabolic process	biological_process	0.0027964	0.0069588
GO:0001932	regulation of protein phosphorylation	biological_process	0.0026933	0.024195
GO:0001934	positive regulation of protein phosphorylation	biological_process	0.0018498	0.0051221
GO:0045937	positive regulation of phosphate metabolic process	biological_process	0.0027964	0.0069588
GO:0042325	regulation of phosphorylation	biological_process	0.014019	0.029001
GO:0042327	positive regulation of phosphorylation	biological_process	0.0017908	0.0082935
GO:0071621	granulocyte chemotaxis	biological_process	0.012215	0.030822
GO:0030595	leukocyte chemotaxis	biological_process	0.012618	0.023316
GO:0000165	MAPK cascade	biological_process	6.14E-06	0.026917
GO:0043410	positive regulation of MAPK cascade	biological_process	0.0029229	0.027081
GO:0097305	response to alcohol	biological_process	0.021317	0.00094936
GO:0008283	cell proliferation	biological_process	8.95E-05	0.019719
GO:0042127	regulation of cell proliferation	biological_process	1.70E-05	0.026033

Table 2

The representive results of GO Enrichment Analysis with differential LncRNA targets
GO accession	Description	Term type	P value (AvsB)	P value (CvsD)
GO:0006468	protein phosphorylation	biological_process	3.85E-06	0.00014357
GO:0006793	phosphorus metabolic process	biological_process	1.83E-06	1.33E-07
GO:0006796	phosphate-containing compound metabolic process	biological_process	1.02E-06	1.88E-07
GO:0019220	regulation of phosphate metabolic process	biological_process	0.017031	0.0019835
GO:0019902	phosphatase binding	molecular_function	0.00096843	0.0041649
GO:0042325	regulation of phosphorylation	biological_process	0.026018	0.0036226
GO:0046434	organophosphate catabolic process	biological_process	0.030242	0.00059334
GO:0051174	regulation of phosphorus metabolic process	biological_process	0.018107	0.0018457
GO:0007049	cell cycle	biological_process	6.61E-05	0.011614
GO:0022402	cell cycle process	biological_process	0.0020298	0.020299
GO:0031625	ubiquitin protein ligase binding	molecular_function	0.003529	0.022639
GO:0006914	autophagy	biological_process	0.038656	0.022267
GO:0030111	regulation of Wnt receptor signaling pathway	biological_process	0.0020955	0.022372
GO:0030178	negative regulation of Wnt receptor signaling pathway	biological_process	0.023016	0.031455

3.7 The expression in the TCGA and the immune infiltration status.

We obtained the RNA-seq data of hepatocellular carcinoma from the TCGA database (https://portal.gdc.cancer.gov/), including the gene expression information of 373 samples of hepatocellular carcinoma patients. The combination therapy group and the control group included a total of 886 significantly different mRNA overlaps between group A and group B, group C and group D. Disturbances in the expression of these genes may indicate different clinical outcomes in patients with hepatocellular carcinoma. We used 886 overlapping differentially expressed genes (DEGs) to perform heat map cluster analysis of TCGA hepatocellular carcinoma patients, and we found that the expression of these differential genes can divide 373 hepatocellular carcinoma patients into two distinct categories. We define these two categories as ferroptosis-sensitive class 1 and ferroptosis insensitive class 2 (Fig. 7A, B). At the same time, we evaluated the immune cell infiltration ratio of all hepatocellular carcinoma samples from TCGA using the Xcell database. For class1 and class2, the scores of immune cell abundance, stromal cell abundance, and tumor microenvironment scores all show significant differences. It indicates that the treatment of the RSL3@O₂-ICG NBs + LFUS treatment group can adjust the proportion of immune cells in the hepatocellular carcinoma tumor microenvironment to induce the sensitivity of ferroptosis further (Fig. 7C).

3.8. Prognostic analysis of the TCGA cohort

We used the tumor microenvironment score obtained from the Xcell database in the previous step and used pearson to screen out 122 differential genes related to the tumor microenvironment score for subsequent analysis. First, in the training set, two prognostic risk genes were obtained, namely "SLC37A2" and "ITGB7" genes, and then these two genes were used to calculate the risk score of ferroptosis sensitivity for each sample, and the training set samples were divided high-risk group and low-risk group (Fig. 7D, E). The log-rank test shows that the prognosis of the high-risk group and the low-risk group is significantly different. Low-risk group means a worse prognosis. In the test set, the log-rank test also found that the low-risk group means a worse prognosis. The above results indicate that the prognostic risk genes "SLC37A2" and "ITGB7" which are significantly related to the tumor microenvironment of hepatocellular carcinoma, may represent new hepatocellular carcinoma ferroptosis inducing markers and have guiding significance for the treatment of hepatocellular carcinoma. The expression of the two prognostic risk genes between the class1 and class2 groups is also significantly different, and both are highly expressed in the class1 group relative to the class2 group (Fig. 8A, B). At the same time, the scatter plot shows that the Pearson correlations between the prognostic risk genes "SLC37A2" and "ITGB7" and the hepatocellular carcinoma tumor microenvironment score are 0.46 and 0.41, respectively (Fig. 8C).

3.9 Weighted co-expression network analysis in the TCGA cohort.

In order to explore the correlation pattern between the gene expression of hepatocellular carcinoma in the class1 and class2 groups in the TCGA data, we use the WGCNA method to hope to find gene co-expression modules related to ferroptosis sensitivity. WGCNA first divides the gene co-expression network of 886 overlapping genes in hepatocellular carcinoma into several highly correlated characteristic modules, representing several groups of highly coordinated gene sets. When the soft threshold is set to 3, we found that four different gene co-expression feature modules can be obtained, namely blue, brown, turquoise, and yellow modules, containing 157, 116, 164, and 60 genes, respectively. Use the gene co-expression similarity matrix to draw a heat map and display it in combination with the cluster tree of the characterization module. It can be seen that the higher the expression similarity between genes, the darker the color. Genes in the same module have a high co-expression pattern, but there are apparent differences between different modules. In order to further determine which module is most related to the iron death-sensitive phenotype, we calculated the pearson correlation coefficients between the feature vector ME of the four modules and the phenotype (Table S4 ). We found that the blue module is most relevant to the ferroptosis-sensitive phenotype. On this basis, we calculated the overall pearson correlation coefficient between the blue module and the gene to be 0.58. In the end, we screened out 6 genes with a correlation greater than 0.4 in the blue module. We believe that these 6 co-expressed genes are the essential genes that affect the sensitivity to ferroptosis of hepatocellular carcinoma, including "CPNE5", "GLIPR2", "FAIM2", "ANKDD1A", "P2RX1", "LYPD5". GO and KEGG analysis further show that they are enriched in membrane-related and glycosylation pathways (Fig. 8D-G).

In recent years, combination therapy has shown good advantages in the treatment of HCC. Combination therapy can overcome the shortcomings of monotherapy and reduce the resistance caused by a single medication(48, 49). Sonodynamic therapy, as a non-invasive ultrasound-induced treatment method, due to the hypoxia and high GSH of HCC TME, the effect of single sonodynamic therapy is not good(49). Ferroptosis is a type of programmed cell death characterized by iron overload and lipid peroxidation. Since ROS of ferroptosis has great similarity with SDT's cytotoxicity, homologous ROS generation therapies bypass the apoptosis pathway and avoid the occurrence of multi-drug resistance (MDR) (50), alleviating the inherent Resistance. Indocyanine green was selected for the construction of nanobubbles as a sonosensitizer. ICG can produce abundant ROS under the stimulation of ultrasound. We constructed RSL3@O₂-ICG NBs that took advantage of the amphiphilic nature of phospholipid molecules, loaded hydrophilic ICG on the surface of the hydrophilic head of the phospholipid, and packaged the hydrophobic RSL3. With the core of oxygen encapsulated in lipid NBs, NBs in this study can also be used as oxygen carrier combined with US to relieve hypoxia in the tumor microenvironment and achieve UTND simultaneously. UTND can selectively deliver the sonosensitizer ICG and the molecularly targeted drug RSL3 to the tumor site through NBs to increase the drug concentration in the tumor site and increase the permeability of cell membranes through the cavitation effect and promote the entry of drugs into cells. Furthermore, it has ultrasound imaging function. In the next experiment, we will add RSL3@O₂-ICG NBs combined with LFUS to achieve the combined treatment of SDT and ferroptosis for HCC.

Firstly, as shown in Fig. 1A, Spherical nanoparticles with a diameter of about 300 nm can be seen by SEM. It was observed that the representative Chlorine element of RSL3 and Sodium Sulfur elements of ICG were enriched on the phospholipid nanobubbles (Fig. 1B-F). The Zeta potential of the nanoplatform was shown in Table 1, and it showed a negative charge at about − 3.86 mV. The size of 300 nm diameter, negative potential and good solubility of the NBs all indicated that might have EPR (Enhanced Permeability and Retention) effect. As shown in (Fig. 1H), NBs had specific absorption of ICG and RSL3 at approximately 770nm and 240nm, indicating the successful incorporation of ICG and RSL3 into RSL3@O₂-ICG NBs, which is consistent with the previous SEM imaging results.

Secondly, the role of O₂-ICG-NBs in enhancing SDT has been proven(Fig. 2B), which proved that the ROS generated by LFUS increased and supported us to carry out the following experiment. Since the intracellular GSH will be reduced due to the cytotoxic •OH produced by the Fe Fenton reaction in ferroptosis, the GSH level is generally used to evaluate the intensity of ferroptosis (51). It can be seen from Fig. 3C and Fig. 3E that compared with the PBS treatment group, the GSH level of the combined treatment group decreased significantly, while the decreased GSH levels of the other two groups treated with SDT or RSL3 alone were considered to be negligible. In order to prove more convincingly that most of the decrease in GSH level in the combination treatment group was caused by ferroptosis, we chose to add Fer-1 (1µM). The results showed that the GSH level was similar to the blank group, indicating that GSH depletion was mainly due to ferroptosis. At the same time, kinds of literature indicate that LPO in cells is an important biomarker of ferroptosis (52, 53), and its accumulation contributes to the occurrence of ferroptosis. As an end product of fat oxidation, MDA has increased expression during the whole process, which can also indicate the occurrence of ferroptosis process (16). The results in Fig. 3B and Fig. 3D show that the MDA content of the combined treatment group increased. It can be seen that there is a negative correlation between GSH and GSH. In contrast, MMP analysis showed similar results; the combination therapy effectively enhanced the green fluorescence of the cells. After determining the method of leading induced cells, we will further explore the effects of cell division through cell cycle experiments; after various treatments, more than half of HepG2/Huh7 cells stagnated in G0/G1 phase in combination therapy group compared to PBS groups (Fig. 4). These data suggested that RSL3@O₂-ICG NBs with LFUS could inhibit cell cycle progression, thus affecting HepG2/ Huh7 cells growth.

Thirdly, RNA-seq analysis further clarify the primary mechanism. Moreover, the target genes of the differential lncRNAs were predicted, and then GO and KEGG function analyses were performed on the differential mRNAs and the predicted target gene of the differential lncRNAs. We found that different mRNAs before and after treatment in different cell lines are collectively enriched in “protein phosphorylation,” “positive regulation of phosphorus metabolic process,” “regulation of protein phosphorylation,” “positive regulation of protein phosphorylation,” “positive regulation of phosphate metabolic process,” “regulation of phosphorylation,” “positive regulation of Phosphorylation” and other protein phosphorylation-related functional pathways (Table 1). It suggests that targeted protein phosphorylation may help increase the sensitivity of HCC to ferroptosis. In addition, mRNA was also enriched in “granulocyte chemotaxis,” “leukocyte chemotaxis,” and other pathways (Table 1). It showed that the treatment caused the change of tumor cell TME and caused the change of related immune cells. To further confirm it, the immune cell infiltration ratio of all hepatocellular carcinoma samples from TCGA was evaluated, which found that the scores of stromal cell abundance and tumor microenvironment scores showed significant differences between the ferroptosis sensitive and insensitive groups that we considered. It shows that combined treatment can adjust the proportion of immune cells in the hepatocellular carcinoma tumor microenvironment to improve ferroptosis sensitivity further. Moreover, the WGCNA weighted co-expression network found six key modular genes related to ferroptosis sensitivity and found them enriched in membrane-related pathways and glycosylation-related pathways. It will become a new target for HCC.

Finally, we analyze the 122 differential genes related to the tumor microenvironment score for subsequent prognostic analysis, then discovered that the prognostic risk genes "SLC37A2" and "ITGB7" which are significantly related to the tumor microenvironment of hepatocellular carcinoma, which may represent new ferroptosis inducing markers for hepatocellular carcinoma, this has guiding significance for the clinical treatment of hepatocellular carcinoma.

However, this study still has some shortcomings and limitations; the potential therapeutic targets obtained by RNA sequencing have not been verified by QPCR, and animal experiments have not been carried out.

In this study, RSL3@O₂-ICG NBs provide a two-in-one treatment platform for SDT and sensitization ferroptosis treatment of liver cancer, showing high anti-tumor efficacy and safety. We found that it was significantly related to the tumor microenvironment of hepatocellular carcinoma. The prognostic risk genes "SLC37A2" and "ITGB7" may represent new hepatocellular carcinoma ferroptosis-inducing markers and have guiding significance for treating hepatocellular carcinoma.

HCC: hepatocellular carcinoma; SDT: sonodynamic therapy; PDT: photodynamic therapy; ROS: reactive oxygen species; LFUS: low-frequency ultrasound; RNA-seq: RNA sequencing; MMP: mitochondrial membrane potential; RSL3: RAS-selective lethal 3; UTND: ultrasound-targeted nanobubble destruction; GPX4: glutathione peroxidase 4.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets used and analyzed during the current study were uploaded and deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE178573 (https://www.ncbi.nlm.nih.gov/geo/query /acc.cgi?acc=GSE14 0729).

Competing interests

The authors have declared that no competing interest exists.

Funding

This research was supported by the National Natural Science Foundation of China (grant number 81873900), Natural Science Foundation of Heilongjiang Province (grant number YQ2020H037), and HAI YAN Science Foundation of Harbin Medical University Cancer Hospital (grant number JJMS2021-30).

Author Contributions

W.C., B.W. conceived this project. B.W., Y.C. designed and supervised the project and commented on the project. H.S., Y.C. and C.W. structured and characterized the nanobubble complexes, Y.C. detected the curative effect in vitro, Z.J. and S.Z. performed the RNA-seq, and Y.C. analyzed the data and wrote the paper. All the authors contributed to the discussion during the whole project.

Acknowledgments

Not applicable.

McGlynn KA, Petrick JL, El-Serag HB. Epidemiology of Hepatocellular Carcinoma. Hepatology (Baltimore, Md). 2021;73 Suppl 1(Suppl 1):4-13.
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: a cancer journal for clinicians. 2018;68(6):394-424.
Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet (London, England). 2018;391(10127):1301-14.
Hollebecque A, Malka D, Ferté C, Ducreux M, Boige V. Systemic treatment of advanced hepatocellular carcinoma: from disillusions to new horizons. European journal of cancer (Oxford, England : 1990). 2015;51(3):327-39.
Wang Q, Zhang P, Li Z, Feng X, Lv C, Zhang H, et al. Evaluation of Polymer Nanoformulations in Hepatoma Therapy by Established Rodent Models. Theranostics. 2019;9(5):1426-52.
Ding J, Chen J, Gao L, Jiang Z, Zhang Y, Li M, et al. Engineered nanomedicines with enhanced tumor penetration. Nano Today. 2019;29.
Zheng P, Liu Y, Chen J, Xu W, Li G, Ding J. Targeted pH-responsive polyion complex micelle for controlled intracellular drug delivery. Chinese Chemical Letters. 2020;31(5):1178-82.
Zhao D, Zhu T, Li J, Cui L, Zhang Z, Zhuang X, et al. Poly(lactic-co-glycolic acid)-based composite bone-substitute materials. Bioact Mater. 2021;6(2):346-60.
Feng X, Xu W, Liu J, Li D, Li G, Ding J, et al. Polypeptide nanoformulation-induced immunogenic cell death and remission of immunosuppression for enhanced chemoimmunotherapy. Science Bulletin. 2021;66(4):362-73.
Wei L, Chen J, Ding J. Sequentially stimuli-responsive anticancer nanomedicines. Nanomedicine (London, England). 2021;16(4):261-4.
Kuczynski EA, Lee CR, Man S, Chen E, Kerbel RS. Effects of Sorafenib Dose on Acquired Reversible Resistance and Toxicity in Hepatocellular Carcinoma. Cancer research. 2015;75(12):2510-9.
Wang H, Lu Z, Zhao X. Tumorigenesis, diagnosis, and therapeutic potential of exosomes in liver cancer. Journal of hematology & oncology. 2019;12(1):133.
Méndez-Blanco C, Fondevila F, García-Palomo A, González-Gallego J, Mauriz JL. Sorafenib resistance in hepatocarcinoma: role of hypoxia-inducible factors. Experimental & molecular medicine. 2018;50(10):1-9.
Shen Z, Song J, Yung BC, Zhou Z, Wu A, Chen X. Emerging Strategies of Cancer Therapy Based on Ferroptosis. Advanced materials (Deerfield Beach, Fla). 2018;30(12):e1704007.
Zhou Z, Song J, Tian R, Yang Z, Yu G, Lin L, et al. Activatable Singlet Oxygen Generation from Lipid Hydroperoxide Nanoparticles for Cancer Therapy. Angewandte Chemie (International ed in English). 2017;56(23):6492-6.
Conrad M, Pratt DA. The chemical basis of ferroptosis. Nature chemical biology. 2019;15(12):1137-47.
Xu T, Ma Y, Yuan Q, Hu H, Hu X, Qian Z, et al. Enhanced Ferroptosis by Oxygen-Boosted Phototherapy Based on a 2-in-1 Nanoplatform of Ferrous Hemoglobin for Tumor Synergistic Therapy. ACS nano. 2020;14(3):3414-25.
Yumita N, Sasaki K, Umemura S, Yukawa A, Nishigaki R. Sonodynamically induced antitumor effect of gallium-porphyrin complex by focused ultrasound on experimental kidney tumor. Cancer letters. 1997;112(1):79-86.
Zhang R, Yan F, Chen Y. Exogenous Physical Irradiation on Titania Semiconductors: Materials Chemistry and Tumor-Specific Nanomedicine. Advanced science (Weinheim, Baden-Wurttemberg, Germany). 2018;5(12):1801175.
Pan X, Wang H, Wang S, Sun X, Wang L, Wang W, et al. Sonodynamic therapy (SDT): a novel strategy for cancer nanotheranostics. Science China Life sciences. 2018;61(4):415-26.
Liu Y, Bai L, Guo K, Jia Y, Zhang K, Liu Q, et al. Focused ultrasound-augmented targeting delivery of nanosonosensitizers from homogenous exosomes for enhanced sonodynamic cancer therapy. Theranostics. 2019;9(18):5261-81.
Qian X, Han X, Chen Y. Insights into the unique functionality of inorganic micro/nanoparticles for versatile ultrasound theranostics. Biomaterials. 2017;142:13-30.
Huang J, Liu F, Han X, Zhang L, Hu Z, Jiang Q, et al. Nanosonosensitizers for Highly Efficient Sonodynamic Cancer Theranostics. Theranostics. 2018;8(22):6178-94.
Li Z, Zhu L, Sun H, Shen Y, Hu D, Wu W, et al. Fluorine assembly nanocluster breaks the shackles of immunosuppression to turn the cold tumor hot. Proc Natl Acad Sci U S A. 2020;117(52):32962-9.
Yang X, Li H, Qian C, Guo Y, Li C, Gao F, et al. Near-infrared light-activated IR780-loaded liposomes for anti-tumor angiogenesis and Photothermal therapy. Nanomedicine : nanotechnology, biology, and medicine. 2018;14(7):2283-94.
Li Z, Wang Y, Shen Y, Qian C, Oupicky D, Sun M. Targeting pulmonary tumor microenvironment with CXCR4-inhibiting nanocomplex to enhance anti-PD-L1 immunotherapy. Science advances. 2020;6(20):eaaz9240.
Li Z, Chen G, Ding L, Wang Y, Zhu C, Wang K, et al. Increased Survival by Pulmonary Treatment of Established Lung Metastases with Dual STAT3/CXCR4 Inhibition by siRNA Nanoemulsions. Mol Ther. 2019;27(12):2100-10.
Li H, Yang X, Gao F, Qian C, Li C, Oupicky D, et al. Bioreduction-ruptured nanogel for switch on/off release of Bcl2 siRNA in breast tumor therapy. Journal of controlled release : official journal of the Controlled Release Society. 2018;292:78-90.
Chen G, Wang K, Wang Y, Wu P, Sun M, Oupicky D. Fluorination Enhances Serum Stability of Bioreducible Poly(amido amine) Polyplexes and Enables Efficient Intravenous siRNA Delivery. Adv Healthc Mater. 2018;7(5).
Lieto E, Galizia G, Cardella F, Mabilia A, Basile N, Castellano P, et al. Indocyanine Green Fluorescence Imaging-Guided Surgery in Primary and Metastatic Liver Tumors. Surgical innovation. 2018;25(1):62-8.
Souzaki R, Kawakubo N, Matsuura T, Yoshimaru K, Koga Y, Takemoto J, et al. Navigation surgery using indocyanine green fluorescent imaging for hepatoblastoma patients. Pediatric surgery international. 2019;35(5):551-7.
Yang J, Tao HS, Cai W, Zhu W, Zhao D, Hu HY, et al. Accuracy of actual resected liver volume in anatomical liver resections guided by 3-dimensional parenchymal staining using fusion indocyanine green fluorescence imaging. Journal of surgical oncology. 2018;118(7):1081-7.
Peyrat P, Blanc E, Guillermet S, Chen Y, Ferlay C, Perol D, et al. HEPATOFLUO: A prospective monocentric study assessing the benefits of indocyanine green (ICG) fluorescence for hepatic surgery. Journal of surgical oncology. 2018;117(5):922-7.
Nguyen Cao TG, Kang JH, You JY, Kang HC, Rhee WJ, Ko YT, et al. Safe and Targeted Sonodynamic Cancer Therapy Using Biocompatible Exosome-Based Nanosonosensitizers. ACS applied materials & interfaces. 2021;13(22):25575-88.
Ma R, Wu Q, Si T, Chang S, Xu RX. Oxygen and Indocyanine Green loaded microparticles for dual-mode imaging and sonodynamic treatment of cancer cells. Ultrasonics sonochemistry. 2017;39:197-207.
Liu Z, Li J, Chen W, Liu L, Yu F. Light and sound to trigger the Pandora's box against breast cancer: A combination strategy of sonodynamic, photodynamic and photothermal therapies. Biomaterials. 2020;232:119685.
Wu P, Sun Y, Dong W, Zhou H, Guo S, Zhang L, et al. Enhanced anti-tumor efficacy of hyaluronic acid modified nanocomposites combined with sonochemotherapy against subcutaneous and metastatic breast tumors. Nanoscale. 2019;11(24):11470-83.
Xie W, Zhu S, Yang B, Chen C, Chen S, Liu Y, et al. The Destruction Of Laser-Induced Phase-Transition Nanoparticles Triggered By Low-Intensity Ultrasound: An Innovative Modality To Enhance The Immunological Treatment Of Ovarian Cancer Cells. International journal of nanomedicine. 2019;14:9377-93.
Chen Q, Liang C, Wang C, Liu Z. An imagable and photothermal "Abraxane-like" nanodrug for combination cancer therapy to treat subcutaneous and metastatic breast tumors. Advanced materials (Deerfield Beach, Fla). 2015;27(5):903-10.
Proulx ST, Luciani P, Derzsi S, Rinderknecht M, Mumprecht V, Leroux JC, et al. Quantitative imaging of lymphatic function with liposomal indocyanine green. Cancer research. 2010;70(18):7053-62.
Liu Y, Xu M, Chen Q, Guan G, Hu W, Zhao X, et al. Gold nanorods/mesoporous silica-based nanocomposite as theranostic agents for targeting near-infrared imaging and photothermal therapy induced with laser. International journal of nanomedicine. 2015;10:4747-61.
Ma Y, Sadoqi M, Shao J. Biodistribution of indocyanine green-loaded nanoparticles with surface modifications of PEG and folic acid. International journal of pharmaceutics. 2012;436(1-2):25-31.
Khan MS, Hwang J, Seo Y, Shin K, Lee K, Park C, et al. Engineering oxygen nanobubbles for the effective reversal of hypoxia. Artificial cells, nanomedicine, and biotechnology. 2018;46(sup3):S318-s27.
Tan H, Tian Y, Yang H, Liu Z, Liang X, Li B, et al. Oxygen-sufficient lipid nanobubbles combined with UTMD for enhanced sonodynamic therapy of Hep-G2 cells. J Biomed Mater Res B Appl Biomater. 2021.
Wu B, Yuan Y, Liu J, Shang H, Dong J, Liang X, et al. Single-cell RNA sequencing reveals the mechanism of sonodynamic therapy combined with a RAS inhibitor in the setting of hepatocellular carcinoma. Journal of nanobiotechnology. 2021;19(1):177.
Deng F, Sharma I, Dai Y, Yang M, Kanwar YS. Myo-inositol oxygenase expression profile modulates pathogenic ferroptosis in the renal proximal tubule. The Journal of clinical investigation. 2019;129(11):5033-49.
Doll S, Freitas FP, Shah R, Aldrovandi M, da Silva MC, Ingold I, et al. FSP1 is a glutathione-independent ferroptosis suppressor. Nature. 2019;575(7784):693-8.
Xu M, Zhou L, Zheng L, Zhou Q, Liu K, Mao Y, et al. Sonodynamic therapy-derived multimodal synergistic cancer therapy. Cancer letters. 2021;497:229-42.
Liang S, Deng X, Ma P, Cheng Z, Lin J. Recent Advances in Nanomaterial-Assisted Combinational Sonodynamic Cancer Therapy. Advanced materials (Deerfield Beach, Fla). 2020;32(47):e2003214.
Tang D, Kroemer G. Ferroptosis. Current biology : CB. 2020;30(21):R1292-r7.
Stockwell BR, Jiang X, Gu W. Emerging Mechanisms and Disease Relevance of Ferroptosis. Trends in cell biology. 2020;30(6):478-90.
Xie Y, Hou W, Song X, Yu Y, Huang J, Sun X, et al. Ferroptosis: process and function. Cell death and differentiation. 2016;23(3):369-79.
Stockwell BR, Friedmann Angeli JP, Bayir H, Bush AI, Conrad M, Dixon SJ, et al. Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell. 2017;171(2):273-85.

Additionalfile1FigureS1.tif
Additional file 1: Figure S1. Analysis of ROS.
Additionalfile2FigureS2.tif
Additional file 2: Figure S2. The function of RSL3@O2-ICG NBs in ferroptosis.
Additionalfile3FigureS3.tif
Additional file 3: Figure S3. The Venn diagram of differential expressions.
Additionalfile4TableS1.docx
Additional file 4: Table S1. Zeta potential of RSL3@O2-ICG NBs;
Additionalfile5TableS2.xls
Additional file 5: Table S2. Target gene prediction of differentially expressed LncRNA in A vs. B.
Additionalfile6TableS3.xls
Additional file 6: Table S3. Target gene prediction of differentially expressed LncRNA in C vs. D.
Additionalfile7TableS4.docx
Additional file 7: Table S4. The pearson coefficient between the four co-expression modules and the phenotype.
Additionalfile8Supplementaryfigurelegends.docx
Additional file 8: Supplementary Figure Legends;

Download PDF

Journal Publication

published 01 Jan, 2022

Read the published version in International Journal of Nanomedicine →

Version 1

posted

You are reading this latest preprint version

RNA-seq Explores the Mechanism of Oxygen-Boosted Sonodynamic Therapy Based on All-in-One Nanobubbles to Enhance Ferroptosis for the Treatment of HCC

Status:

Journal Publication

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

1. Background

2. Methods

2.1 chemicals

2.2 Preparation of RSL3@O2-ICG NBs

2.3 Characterization

2.4 cell culture

2.5 Cell viability of HepG2 and Huh7 cells

2.6 ROS detection

2.7 Analysis of GSH

2.8 MDA assay

2.9 Measurement of Mitochondrial Membrane Potential

2.10 Cell Cycle Analysis

2.11. Whole transcriptome library construction, sequencing, and analysis

2.11.1 Whole transcriptome library construction and sequencing

2.11.2 RNA-Seq data analysis

2.11.3 Differential expression analysis

2.11.4 GO and KEGG enrichment analysis

2.11.5 Expression and immune infiltration in TCGA.

2.11.6. Prognostic analysis of TCGA cohort

2.12 Statistical analysis

3. Results

3.1 Preparation and characterization of RSL3@O2-ICG NBs.

3.2 Detection of intracellular ROS of SDT.

3.3 The enhancement effect of ferroptosis.

3.4 Influence on Cell Cycle Progression In Vitro.

3.5 Differential expression analysis based on RNA-seq.

3.6 Differentially expressed analysis.

3.7 The expression in the TCGA and the immune infiltration status.

3.8. Prognostic analysis of the TCGA cohort

3.9 Weighted co-expression network analysis in the TCGA cohort.

4. Discussion

5. Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1

2.2 Preparation of RSL3@O₂-ICG NBs

3.1 Preparation and characterization of RSL3@O₂-ICG NBs.