Analysis of CRISPR/Cas9 screening in HCC Lenvatinib-resistance cells
We established HCC Lenvatinib-resistance cell lines and performed CRISPR/Cas9 Screen on the model so that we could identify the key regulatory factors regulating HCC Lenvatinib-resistance. We confirmed resistance to Lenvatinib of HCC cell lines by analysis of IC50(Figure S1C). Using CRISPR/Cas9 high-throughput screening technology for positive or negative screening in tumor cells, candidate gene-directed RNA (gRNA) can be enriched to search for target genes of interest. The complete process of this technology includes the following steps: lentiviral library construction, virus infection of cells within the library, experimental screening of cells, genomic DNA extraction and second-generation sequencing library construction, second generation sequencing and bioinformatics analysis (Figure 1A). For PCR amplification products without any sequencing adapter sequence, after passing the quality inspection, take 100~500ng DNA fragments and go through End Repair, 5'Phosphorylation and dA-Tailing, Adapter Ligation, Clean Up, PCR Enrichment and Clean Up and other steps to construct and form a high-throughput sequencing library on the Illumina platform (Figure 1B). After obtaining raw data with sufficient sequencing depth, bioinformatic analysis was performed, which consists of the following three main stages. The detailed analysis process was shown in Figure 1C. We used volcano plots to show genes with significant differences in positive and negative screening results (Figure 1D-1E). GO analysis revealed that the positive group was mainly enriched in ERBB, cell development etc., while the negative group was mainly enriched in fatty acid, immune response etc. (Figure 1F, Figure S1A). KEGG analysis indicated that the positive group was mainly enriched in signaling by RANBP2, KRAS etc., while the negative group was mainly enriched in NF-κB ,TP53 etc. (Figure 1G, Figure S1B ). These results suggested that CRISPR/Cas9 screening technology indeed provides a large number of reliable differential genes for promoting or inhibiting HCC Lenvatinib-resistance, which can be targeted for subsequent studies.
Whole transcriptome results of HCC Lenvatinib-resistance cells
In this study, we further performed whole-transcriptome sequencing including coding and non-coding RNAs of HCCLM3-LR and HEP3B-LR cell lines. For mRNAs, heat map and volcano map showed differential genes in susceptible and resistant groups, respectively (Figure 2A-2B). Venn diagrams compared the differential genes common to and unique to different comparison combinations (Figure 2C). GO analysis revealed that the different genes in HCCLM3-LR/LS was mainly enriched in Developmental process, Chromosome, Molecular function etc., while the different genes in HEP3B-LR/LS was mainly enriched in Metabolic process, Intracellular, Protein binding etc. (Figure 2D-2E). KEGG analysis indicated that the different genes in HCCLM3-LR/LS was mainly enriched in Apoptosis, Alcoholism, Toll like receptor signaling pathway etc., while the different genes in HEP3B-LR/LS was mainly enriched in Systemic lupus erythematosus, Alcoholism, Necroptosis etc. (Figure 2F-2G).
For lncRNAs, heat map and volcano map showed differential genes in susceptible and resistant groups, respectively(Figure S2A-2B).Venn diagrams compared the differential genes common to and unique to different comparison combinations (Figure S2C).GO analysis revealed that the different genes in HCCLM3-LR/LS was mainly enriched in Phosphorus metabolic process, Mitochondrial membrane, NADH dehydrogenase(ubiquinone) activity etc., while the different genes in HEP3B-LR/LS was mainly enriched in Regulation of metabolic process, Non-membrane bounded organelle, Molecular function etc (Figure S2D-S2E). KEGG analysis indicated that the different genes in HCCLM3-LR/LS was mainly enriched in Thermogenesis, Oxidative phosphorylation, Chemokine signaling pathway etc., while the different genes in HEP3B-LR/LS was mainly enriched in Necroptosis, Focal adhesion, TNF signaling pathway etc. (Figure S2F-S2G).
For circRNAs, volcano map showed differential genes in susceptible and resistant groups, respectively (Figure S3A). GO analysis revealed that the different genes in HCCLM3-LR/LS was mainly enriched in Cellular metabolic process, Intracellular, Protein binding etc, while the different genes in HEP3B-LR/LS was mainly enriched in Metabolic process, Organelle, Molecular function etc(Figure S3B). KEGG analysis indicated that the different genes in HCCLM3-LR/LS was mainly enriched in Ubiquitin mediated proteolysis, Lysine degradation, Autophagy etc., while the different genes in HEP3B-LR/LS was mainly enriched in Protein process in endoplasmic reticulum, Cellular senescence, Cell cycle etc (Figure S3C). According to the prediction results of miRNA binding sites, the regulatory network of miRNA-circRNA was constructed (Figure S3D).
For miRNAs, we demonstrated the first base bias of known miRNAs ranging in length from 18 to 30nt in four cells (HCCLM3-LR/LS, HEP3B-LR/LS), respectively (Figure 3A). Heat map and volcano map showed differential genes in susceptible and resistant groups, respectively(Figure 3B-3C). GO analysis revealed that the different genes in HCCLM3-LR/LS was mainly enriched in Biological process, Intracellular, Molecular function etc., while the different genes in HEP3B-LR/LS was mainly enriched in Cellular process, Intracellular, Molecular function etc (Figure 3D,3F). KEGG analysis indicated that the different genes in HCCLM3-LR/LS was mainly enriched in VEGF signaling pathway, Mitophagy, Apoptosis, while the different genes in HEP3B-LR/LS was mainly enriched in VEGF, MAPK pathway,etc. (Figure 3E, 3G).
We also performed a joint analysis of the above RNAs and showed the KEGG analysis diagram. miRNA-mRNA joint showed that the different genes in HCCLM3-LR/LS was mainly enriched in Alcoholism, Tyrosine metabolism, Metobolic pathways etc., while the different genes in HEP3B-LR/LS was mainly enriched in MAPK signaling pathway, MicroRNAs in cancer, P53 signaling pathway etc (Figure S4A). circRNA-mRNA joint showed that the different genes in HCCLM3-LR/LS was mainly enriched in TGF-βsignaling pathway, Hippo signaling pathway, Cell cycle etc., while the different genes in HEP3B-LR/LS was mainly enriched in Serotonergic synapse, RAS signaling pathway, Hippo signaling pathway etc (Figure 4). LncRNA-mRNA joint showed that the different genes in HCCLM3-LR/LS was mainly enriched in TNF signaling pathway, P53 signaling pathway, Ferroptosis etc., while the different genes in HEP3B-LR/LS was mainly enriched in TNF signaling pathway, Focal adhesion, Autophagy etc (Figure S4C). LncRNA-miRNA-mRNA joint showed that the different genes in HCCLM3-LR/LS was mainly enriched in MicroRNAs in cancer, MAPK signaling pathway, P53 signaling pathway etc, while the different genes in HEP3B-LR/LS was mainly enriched in Metabolic pathways, Histidine metabolism, Serotonergic synapse etc., (Figure S4D). CircRNA-miRNA-mRNA joint showed that the different genes in HCCLM3-LR/LS was mainly enriched in MicroRNAs in cancer, MAPK signaling pathway, P53 signaling pathway etc., while the different genes in HEP3B-LR/LS was mainly enriched in Metabolic pathways, Histidine metabolism, Serotonergic synapse etc. (Figure S4E). Whole transcriptome sequencing systematically studies various RNA molecules and their interactions, providing a rich database for HCC Lenvatinib-resistance.
ASB2 facilitated the progressioninHCC Lenvatinib-resistance cells
Based on the cross-analysis of HCCLM3-LR and HEP3B-LR cell lines, we found that ASB2 was involved in the resistance of both cell lines, so it has the potential to be further studied. Three shRNAs against ASB2 (shASB2) were developed to silence ASB2 in HCCLM3-LR and HEP3B-LR cells, which was confirmed by qRT-PCR and western blotting and sh3-ASB2 was selected for further experiment (Figure 4A-4B). shASB2 inhibited proliferation in HCCLM3-LR and HEP3B-LR cells, as per the results of EdU assays (Figure 4C). Transwell and wound healing assays showed that shASB2 significantly hindered the invasion and migration functions of cancer cells (Figure 4D-4E). These results demonstrated that ASB2 promoted the progression of HCC Levatinib-resistant cells in vitro.
ASB2activated NF-κB pathway by promoting IκBα ubiquitination to mediate HCC Lenvatinib-resistance
The above results raised our curiosity about the reason why ASB2 promotes HCC Lenvatinib-resistance. We enriched the CRISPR/Cas9 results and surprisingly found differential promoting genes were enriched NF-κB pathway (Figure 1G ). Continued activation of NF-κB is thought to be associated with cancer resistance, including Sorafenib resistance in HCC [9-10]. So we wanted to explore how does ASB2 activate NF-κB and promote HCC Lenvatinib-resistance. STRING analysis predicted that ASB2 protein was associated with NF-κB pathway-related proteins(Figure S5). Based on the fact that ASB2 is a ubiquitinase, we hypothesized that ASB2 activates the NF-κB pathway by promoting the ubiquitination and degradation of IκBα. After transfection with ASB2 plasmid, the ubiquitination level of IκBα increased, and the expression of IκBα decreased in 293T cells (Figure 5A). We further examined the effect of ASB2 knockdown on IκBα ubiquitination in HCCLM3-LR and HEP3B-LR cells and found that the ubiquitination of IκBα was decreased by ASB2 knockdown and the expression of IκBα was increased (Figure 5B).To verify the interaction between ASB2 and IκBα, we performed co-immunoprecipitation experiments. As shown in Figure 5C, ASB2 and IκBα interacted with each other in 293 T cells. We observed co-localization of ASB2 and IκBα through immunofluorescence experiments in 293 T cells(Figure 5C-5D). In addition ,western blot confirmed that the expression of P-P65 was down-regulated in shASB2 group but P65 was not differentially expressed, suggesting that the NF-κB pathway was inhibited (Figure 5E).These results verified ASB2 promoted the ubiquitination of IκBα, activating NF-κB pathway in HCC Lenvatinib-resistance cells.
ASB2 inhibitedferroptosis by promoting P53 ubiquitination to mediate HCC Lenvatinib-resistance
CRISPR/Cas9 enrichment results also revealed differential gene enrichment in the fatty acid pathway (Figure 1F). Ferroptosis is a common type of fatty acid metabolism [11-12]. We found that ASB2 production was significantly negatively correlated with GPX4 in HCC via GEPIA database (Figure S6). Therefore, we tried to investigate that whether ASB2 mediates ferroptosis to promote drug resistance. To prove our hypothesis, we first examined the ROS content in shASB2 and shNC groups. The ROS probe suggested that the ROS content in shASB2 group was significantly higher than that of the control group (Figure 6A). Electron micrographs showed that the mitochondria of shASB2 group were pyknotic, smaller in size, higher in electron density, and expanded in cristae, consistent with the characteristics of ferroptosis (Figure 6B). So how does ASB2 affect ferroptosis deserves further study. Ubibrowser predicted that P53 may be degraded by ASB2 ubiquitination (Figure S7). Recent studies have found that P53 acts as a positive regulator of ferroptosis by promoting ROS production. P53 directly regulates the metabolic versatility of cells by favoring mitochondrial respiration, leading to ROS-mediated ferroptosis[13-14]. After transfection with ASB2 plasmid, the ubiquitination level of P53 increased in 293T cells (Figure 6C). We detected the effect of ASB2 knockdown on P53 ubiquitination and results showed that the ubiquitination of P53 was decreased by ASB2 knockdown and the expression of P53 was increased in HCCLM3-LR and HEP3B-LR cells (Figure 6D). Western blot confirmed that ferroptosis-related proteins including GPX4 and SLC7A11 proteins were all down-regulated in shASB2 group, indicating ferroptosis increased and iron inducers-erastin was used as a control (Figure 6E). These results confirmed ASB2 promoted the ubiquitination of P53, inhibiting ferroptosis pathway in HCC Lenvatinib-resistance cells.
NOTCH1 transcriptionally regulated the expression of ASB2 inHCC Lenvatinib-resistance cells
The above studies confirmed the downstream pathway of ASB2, so who regulates the expression of ASB2? Previous study indicated that NOTCH1 could induced ASB2 expression to promote protein ubiquitination by forming non-canonical E3 ligase complexes [15]. Therefore, we hypothesized that NOTCH1 activates the expression of ASB2 through transcriptional regulation and then participates in NF-κB and ferroptosis pathways leading to persistent drug resistance in HCC Lenvatinib resistance cells. Results showed that the addition of NOTCH1 inhibitor (LY3039478, iNOTCH1) decreased the proliferation, invasion and migration in HCCLM3-LR and HEP3B-LR cells (Figure S8). Both ASB2 and ferroptosis genes including GPX4 and SLC7A11 in mRNA level were significantly reduced after the addition of NOTCH1 inhibitor (Figure 7A). Western blot results indicated that in iNOTCH1 group, ferroptosis-related proteins were all down-regulated and the coexistence of iNOTCH1 and shASB2 further enhanced ferroptosis (Figure 7B). The expression of P-P65 was down-regulated in iNOTCH1 group but P65 was not differentially expressed, and the expression of IκBα was increased (Figure 7B).The ROS probe suggested that the ROS content in iNOTCH1 group was significantly higher than that of the control group (Figure 7C). Electron microscopic observation showed that the mitochondria of iNOTCH1 group were pyknotic, small in size, and expanded in cristum, which were consistent with the characteristics of ferroptosis (Figure 7D). The ChIP data from GSE39263 showed that NOTCH1 had a significant peak in the upstream of ASB2, suggesting that there was an obvious binding site between NOTCH1 and ASB2 (Figure 7E).
Knocked down of ASB2 inhibited the growth of HCC Lenvatinib-resistance tumors in vivo
To explore the relevance of ASB2 to HCC growth in vivo, the shASB2 and shNC HCCLM3-LR and HEP3B-LR cells were injected into the axilla of nude mice, respectively. Results revealed that the down-regulated expression of ASB2 significantly inhibited the growth of HCCLM3-LR and HEP3B-LR compared to the shNC group (Figure 8A-8C). Immunofluorescence results showed that shASB2 decreased the expression of SLC7A11 and GPX4 in tumor tissues, while the group of shNC showed increased the expression of SLC7A11 and GPX4 (Figure 8D). These evidences indicated that knocked down of ASB2 promoted ferroptosis and inhibited the growth of HCC Lenvatinib-resistance tumors in vivo.
Venetoclax acting as ASB2 inhibitor decreased the progression of HCC Lenvatinib-resistance tumors in vivo
As there is no readily available drug to inhibit ASB2, we used virtual screening to predict small molecule compounds that might inhibit ASB2. A total of 100 small molecule compounds with potential binding to ASB2 were screened, and the first 20 were shown (Figure S89). After evaluating it, we found that Venetoclax can bind to ASB2 and has a potential inhibitory function using virtual screening (Figure 8E). According to the results of Edu, transwell and wound healing assays, Venetoclax significantly inhibited the proliferation, invasion and migration of HCCLM3-LR and HEP3B-LR cells (Figure S10). In vivo experiments confirmed that Venetoclax inhibited the growth of HCCLM3-LR and HEP3B-LR compared to the PBS group. The combination of Venetoclax and Lenvatinib significantly inhibited the progression of HCC, and the efficacy was better than Lenvatinib alone (Figure 8F-8H). These results demonstrate that Venetoclax as an ASB2 inhibitor can reduce the progression of HCC Lenvatinib-resistant tumors.