Molecular Characterization of the Oncogene BTF3 and Its Targets in Colorectal Cancer

Introduction

Colorectal cancer (CRC) is one of the most commonly diagnosed and leading causes of cancer mortality worldwide (Bray et al., 2018). Although extensive progress has been made in chemotherapy, radiotherapy, immunotherapy, and surgery, the prognosis of patients with CRC is still unsatisfactory (Miller et al., 2016). Genetic and epigenetic molecular alterations at the transcriptional and post-transcriptional level are related to CRC occurrence, progression, treatment resistance, and patient outcomes (Keum and Giovannucci, 2019; Jung et al., 2020; Nguyen et al., 2020). A deeper understanding of the molecular mechanisms of CRC will provide novel insights into the pathogenesis of the disease and identify new treatment options.

Basic transcription factor 3 (BTF3), also known as nascent-polypeptide-associated complex (NAC) beta (NACB), was reported to serve as an oncogene and convey worse prognosis in gastric cancer (Liu et al., 2013; Zhang et al., 2017), pancreatic cancer (Kusumawidjaja et al., 2007), osteosarcoma (Liu et al., 2019a), cervical cancer (Wu et al., 2020), hypopharyngeal squamous cell carcinoma (Zhang et al., 2019), prostate cancer (Hu et al., 2019), breast cancer (Ding et al., 2019), and CRC (Liu et al., 2019b). BTF3 initiates transcription by binding to promoter elements, such as TATA box and CAAT box sequences, in the promoter region (Cavallini et al., 1988; Kanno et al., 1992). Due to alternative splicing, BTF3 is present in two different isoforms, BTF3a and BTF3b. BTF3a is the transcriptionally active form of BTF3, while the isoform lacking the first 44 amino acids of the BTF3a N-terminus, BTF3b, is transcriptionally inactive (Parvin et al., 1994). Interestingly, two distinct roles of BTF3 in cancer have been reported. On the one hand, it serves as a transcriptional regulator by targeting the estrogen receptor (ER; Ding et al., 2019) and tumor associated genes (Kusumawidjaja et al., 2007). On the other hand, it is associated with NAC subunit alpha (NACA). The binding of BTF3 with NACA prevents inappropriate targeting of non-secretory nascent polypeptides produced by ribosomes to the endoplasmic reticulum (Koplin et al., 2010). The NAC complex is linked to protein ubiquitination (Panasenko et al., 2009) and proteolysis (Kirstein-Miles et al., 2013) and has been reported to have a related role with BTF3, such as stabilizing BMI1 (Hu et al., 2019) and binding to the N-terminal domain of ER (Green et al., 2007). In our previous study, we reported the oncogenic function of BTF3 in CRC cell lines (Li et al., 2019). According to Liu et al. (2019b), reduced MAD2L2, MCM3, and PLK1 activity may be involved. However, the direct targets of BTF3 as a transcription factor and part of the NAC complex, and the reason for its elevated expression in CRC, remain obscure.

In this study, we first identified the transcriptional targets of BTF3 by applying combined RNA-Seq and chromatin immunoprecipitation (ChIP)-Seq analysis. Thereafter, we conducted immunoprecipitation (IP)-MS and E3 ubiquitin ligase analysis to identify potential interacting targets of BTF3-NACA. Finally, miRNAs targeting BTF3 were predicted and validated. This study provides a comprehensive basis for understanding the oncogenic mechanism of BTF3 in CRC.

Materials and Methods

Cell Culture

The colon cancer HT-29 cell line was purchased from the American Type Culture Collection (Manassas, VA, United States). All cells were grown in complete Dulbecco’s Modified Eagle’s Medium (DMEM; HyClone; GE Healthcare Life Sciences, Logan, UT, United States) supplemented with 10% fetal bovine serum (HyClone; GE Healthcare Life Sciences), and 100 U/ml penicillin and streptomycin (Gibco; Thermo Fisher Scientific, Waltham, MA, United States) and incubated at 37°C in a 5% CO₂ atmosphere.

Reverse Transcription Quantitative Polymerase Chain Reaction (RT qPCR)

Total RNA was extracted from BTF3-knockdown cell lines (those transfected with BTF3-targeted shRNA), the empty vector control cell line (vector), and non-transfected cells using an RNeasy mini kit (Qiagen China Co., Ltd., Shanghai, China) according to the manufacturer’s protocol. cDNA was generated by reverse transcription of 1-μg aliquots of RNA using the Takara PrimeScript RT Reagent kit (Takara Biotechnology Co., Ltd., Dalian, China) according to the manufacturer’s protocol. cDNA was used for qPCR using the SYBR Premix Ex Taq kit (Takara Biotechnology Co., Ltd.) on a CFX96 instrument qPCR system (Bio Rad Laboratories, Inc., Hercules, CA, United States). PCR was performed according to the manufacturer’s instructions. Initial denaturation was at 95°C for 10 min followed by 30 cycles at 95°C for 1 min, annealing at 53°C for 1 min, extension at 72°C for 1 min, and final extension at 72°C for 5 min. All expression data were normalized to β-actin levels using the 2-ΔΔCq method. Primer sequences were as follows: β-actin, 5′-CGAGCGCGGCTACAGCT-3′ (forward) and 5′-TCCTTAATGTCACGCACGATTT-3′ (reverse); and BTF3 5′-AGCTTGGTGCGGATAGTCTGA-3′ (forward) and 5′-GTGCTTTTCCATCCACAGATTG-3′ (reverse).

Cell Counting Kit-8 Assay

Cellular viability was determined by the Cell Counting Kit-8 (CCK8; Beyotime Biotechnology). Briefly, 100 μl of H1299 lung cancer cells per well was plated into 96-well plates. After treatment with corresponding reagents, 10 μl of CCK8 solution was incubated with cell medium for 2 h at 37°C. The absorbance of each well was detected at 450 nm by Multiskan FC Microplate spectrophotometer.

Transwell Assay

Cell migration assay was performed using a modified Boyden chamber plate with 8-μm pore size polycarbonate membrane filters (Corning Incorporated, Corning, NY, United States). ShBTF3- and empty vector-transfected HT-29 cells, along with non-transfected HT-29 cells, were incubated at 37°C in DMEM for 6 h. Subsequently, 1 × 10⁴ cells from each group and cell type were added to the upper part of the Boyden chamber, and the bottom chamber was filled with DMEM containing 20% serum. Cells were allowed to migrate to the underside of the membrane during incubation for 48 h at 37°C. Next, cells on the membrane filter were fixed with 4% paraformaldehyde and stained with 0.05% Giemsa (Sigma-Aldrich; Merck KGaA). The migration index was defined as the number of cells that migrated to the membrane filter by cell counting in at least three random fields per filter using a light microscope (magnification, ×200).

Wound Healing Assay

Six-well plates were chosen for the wound healing assay. When HT-29 cells were cultured for 24 h after transfection, a plastic pipette tip was used to scratch a line across the cell surface in each plate. The remaining cells were washed three times with PBS to remove the floating cells and debris. Images of the healing process were digitally captured 0 and 24 h after wounding.

Nude Mouse Xenograft Model

Four-week-old nude female mice were obtained from Slaccas Company. A total of 2 × 10⁶ HT-29 cells of each group were subcutaneously injected 12 h after transfection into mice. For the same cell line, the BTF3 knockdown and normal control groups were implanted into the left posterior flank with the knockdown group in the right of the same mouse. Tumor size was measured every week by calipers. Three weeks later, mice were euthanized, and tumors were statistically analyzed. This study was approved by the Ethics Committees of Changhai Hospital (approval no. 20171001089).

RNA Sequencing Analysis

We utilized a quartile normalization algorithm to subtract and correct the background. Then, we used the R software Limma package to screen differentially expressed genes (DEGs) by filtering p value of Student’s t-test and fold change (FC; Diboun et al., 2006). Finally, with a threshold of p value <0.05 and absolute value of FC >2, a volcano plot was created using the R software ggplot2 package to identify DEGs with statistical significance between two groups. Hierarchical clustering and combined analyses were performed to identify DEGs. Gene Ontology (GO) enrichment analyses of differentially expressed mRNAs were implemented with EnrichR (Kuleshov et al., 2016). GO terms, including Molecular Function, Biological Processes, and Cellular Components, with p values less than 0.05 were considered significantly enriched by DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for understanding high-level functions and effects of biological systems.¹ Cytoscape (version 3.40) was used to visualize the protein–protein interaction (PPI) relationships conducted by String database.²

Gene Expression and Tumor Infiltrate Correlation Analysis

Expression of BTF3 and CHD1L in The Cancer Genome Atlas (TCGA), colon adenocarcinoma (COAD), and rectum adenocarcinoma (READ) datasets were conducted using Gene Expression Profiling Interactive Analysis (GEPIA)³ (Tang et al., 2017). Expression of miR-497-5p in Gene Expression Omnibus (GEO) datasets of GSE128446, GSE81581 (Sayagues et al., 2016), and GSE35982 (Fu et al., 2012) were analyzed with the R software Limma package. Gene expression correlation analysis of tumor infiltrates was conducted using TIMER⁴ (Li T. et al., 2017).

ChIP Sequencing Analysis

HT29 cells were treated with 1% formaldehyde and then quenched with glycine for 5 min at room temperature. ChIP assays were performed using a chromatin IP kit (Cell Signaling Technology, Danvers, MA, United States) according to the manufacturer’s instructions. To perform ChIP analysis of BTF3 binding to the CHD1L promoter, the transcriptional start site was identified using the UCSC Genome browser.⁵ The analysis was conducted and is shown with IGV (Integrative Genomics Viewer).

Immunoprecipitation

Immunoprecipitation assays were performed using a Pierce Direct Magnetic IP/Co-IP kit (Thermo Scientific, Rockford, IL, United States) according to the manufacturer’s instructions. Briefly, cells were lysed with IP lysis buffer, and protein levels were quantified. Twenty-five microliters of beads coupled with 5 μg of BTF3 antibody were added into 500-μl lysis solutions (1–2 mg/ml) and incubated overnight at 4°C. Eluted proteins were subjected to mass spectrometry analysis using LC–MS/MS (PTM Biolabs, China). IgG (Cell Signaling Technology, Danvers, MA, United States) was used as a negative control.

Dual Luciferase Reporter Assay

Luciferase reporters were generated based on the psiCHECK2 vector (Promega). The complete 3′UTR of BTF3 mRNA, which included the predicted miR-497-5p binding sites, was PCR amplified and cloned into the psiCHECK2 vector. According to the manufacturer’s guidelines, luciferase reporter genes were co-transfected with miR-497-5p mimics and miR-NC into HT29 cells using Lipofectamine 2000. Dual-Luciferase Reporter Assay System (Promega) and Infinite M200 PRO microplate reader (Tecan) were used to measure the relative luciferase activity.

Statistical Analysis

All results are expressed as the mean ± standard error of the mean (SEM) and were analyzed using the statistical software SPSS 16.0 (SPSS Inc., Chicago, IL, United States). Chi-square test was performed to assess associations between BTF3 expression and CRC clinicopathological parameters. The Kaplan–Meier plot was performed for survival analysis, and the log-rank test was used to estimate differences in survival. A p value <0.05 was considered statistically significant.

Results

The Oncogenic Role and Potential Pathways of BTF3 in CRC

The function of BTF3 in CRC was investigated in vivo using a nude mouse xenograft model. As shown in Figure 1A, BTF3 knockdown HT29 cells exhibited reduced tumor size (42.3%, p < 0.05, N = 6), compared to the control group. Tumor volume was analyzed and is shown in Figure 1B. Consistent with in vitro results in our previous research, the in vivo phenotype supported the conclusion that BTF3 acts as an oncogene in CRC. Thereafter, BTF3 knockdown and vector-transfected HT29 cells were applied for RNA sequencing and bioinformatics analysis. In total, we identified 292 significantly DEGs with the criteria of expression FC >2 or <0.5 and p value <0.05 (Supplementary Table 1). A volcano plot (Figure 1C) and hierarchical clustered heatmap (Figure 1D) were utilized to portray the DEGs. Last, the function of DEGs was analyzed, and the top 10 enriched KEGG pathways are shown in Figure 1E. Neuroactive ligand–receptor interaction, ErbB2, and PPAR were the most significant items. Furthermore, PPI analysis enriched a network composed of AGT, CXCL5, CXCL11, C3AR1, CD36, and so on (Figure 1F).

FIGURE 1

Figure 1. Oncogenic role and potential pathways of BTF3 in colorectal cancer. (A) BTF3 knockdown attenuated colorectal cancer proliferation in a xenograft nude mouse model. (B) Statistical analysis showing tumor volume in BTF3 knockdown and control groups. (C) Volcano plot of DEGs in BTF3 knockdown HT29 cells. (D) Heatmap of DEG expression in BTF3 knockdown and control groups. (E) Functional enrichment of the top 10 pathways in the BTF3 knockdown group. (F) Protein–protein interaction network of DEGs based on the String database. Increased DEGs are shown in red nodes, and decreased DEGs are shown in blue nodes. The larger the node, the higher centrality the gene in the network.

BTF3 Activates CHD1L Transcription in CRC

As a transcription factor, BTF3 has been proven to regulate gene expression by binding DNA elements in cancer. Here, to identify potential targets of BTF3, we conducted ChIP sequencing for BTF3 knockdown and control HT29 cells. In total, ChIP-Seq analysis revealed 149 genes with significant differential peaks between the two groups (Supplementary Table 2). As shown in Figure 2A, combined analysis of RNA-Seq and ChIP-Seq results identified two genes with both differential expression and BTF3 binding capability, Hemicentin 1 (HMCN1) and chromodomain helicase DNA binding protein 1 Like (CHD1L). HMCN1 and CHD1L showed increased and decreased expression, respectively, in the BTF3 knockdown group. As a decreased binding peak of BTF3 was located in the promoter region of CHD1L (Figure 2B) and CHD1L has been reported as an oncogene in CRC (Ji et al., 2013) and other cancer types (Chen et al., 2009, 2010, 2011; Li et al., 2013; Liu et al., 2016), we chose CHD1L as a target candidate. qRT-PCR results further revealed that BTF3 knockdown attenuated CHD1L expression (Figure 2C). Next, we examined the function of BTF3 and CHD1L in CRC cell proliferation, invasion, and migration. As shown in Figure 2D, CCK8 assay shows that BTF3 overexpression promoted, while CHD1L knockdown inhibited, HT29 proliferation. In addition, a rescue assay indicated that CHD1L knockdown reduced proliferation capability in BTF3-overexpressing HT29 cells. Similarly, wound scratch assay and Transwell assay results indicated that BTF3 and CHD1L promoted cell migration and invasion (Figures 2E–H), and CHD1L is a target of BTF3. Significantly higher expression of BTF3 and CHD1L was observed in TCGA CRC datasets (Figures 3A,B), and BTF3 showed an evident positive correlation with CHD1L expression. Additionally, correlation analysis of BTF3 and CHD1L expression with clinical pathological parameters was conducted (Figure 3C). As shown in Table 1, BTF3 demonstrated significant correlations with lymphatic invasion and pathologic stage, while CHD1L was correlated with lymph node count and lymphatic invasion. Finally, expression of BTF3 and CHD1L in tumor-infiltrating immune cells was further analyzed, and both genes showed significant correlation with CD8+ T cells (Figures 3D,E), suggesting their potential role in CD8+ T cell regulation. In summary, these results imply that CHD1L is a transcriptional target of BTF3.

FIGURE 2

Figure 2. CHD1L is a potential transcriptional target of BTF3 in colorectal cancer. (A) Venn diagram showing two common genes, HMCN1 and CHD1L, between significant DEGs with expression and BTF3 binding peak differences. (B) ChIP peaks in the CHD1L gene promoter region in BTF3 control, BTF3 knockdown, and input groups. (C) Expression of CHD1L mRNA expression in BTF3 knockdown groups. (D) CCK8 showing the proliferation effect of BTF3 and CHD1L on HT29 cells. (E,F) Would healing assay showing the migration effect of BTF3 and CHD1L on HT29 cells. (G,H) Transwell assay showing the invasion effect of BTF3 and CHD1L on HT29 cells. *P < 0.05.

FIGURE 3

Figure 3. Clinical significance analysis of BTF3 and CHD1L in CRC. (A) BTF3 demonstrated significantly higher expression in TCGA COAD and READ datasets. (B) CHD1L showed significantly increased expression in TCGA COAD and READ datasets. (C) BTF3 was positively correlated with CHD1L expression. (D) Correlation analysis of BTF3 with tumor purity and tumor-infiltrating immune cells. (E) Correlation analysis of CHD1L with tumor purity and tumor-infiltrating immune cells.

TABLE 1

Table 1. Clinical significance analysis of CHD1L, BTF3, and miR-497-5p expression with clinicopathological parameters in CRC.

BTF3 May Regulate p53 Expression Through the E3 Ligase HERC2

In addition to the transcription factor role of BTF3, how BTF3 acts as a subunit of NAC was also studied. To identify potential BTF3 binding proteins, BTF3 IP and MS were applied (Figure 4A). Using IgG as a negative control, 542 proteins exhibited BTF3 specific expression (Figure 4A and Supplementary Table 3). Indeed, NACA was identified as binding to BTF3 with different scores. Next, GO annotation and enrichment were conducted. As expected, protein targeting to the ER was the most significant biological process (BP; Figure 4B). For cellular component (CC) analysis, proteins in cytosolic ribosome also supported NACA-related function (Figure 4C). Molecular function (MF) analysis enriched RNA binding as the most evident factor (Figure 4D). Of note, the CC analysis enriched 73 nuclear proteins (e.g., TOP2A, DDX3X, DDX47, CHD7, and ADAR) and MF analysis identified 64 proteins with DNA binding potential (e.g., TOP2A, DDX3X, THRB, MCM7, CHD4, and PHF6), which also supports the transcriptional regulatory role of BTF3. Then, as BTF3 was associated with ubiquitination-mediated proteolysis, we analyzed E3 ligases using two databases, Ubibrowser (Li Y. et al., 2017)⁶ and iUUCD 2.0 (Gao et al., 2013; Zhou et al., 2018). The analyses identified an E3 ubiquitin ligase called HECT and an RLD domain containing E3 ubiquitin protein ligase 2 (HERC2) in all three. Receptor for activated C kinase 1 (RACK1) was another protein in the Ubibrowser and BTF3 interacting protein list. Another eight E3 ligases (TRIM28, TRIM25, FUS, RANBP2, RBBP4, GNB2, RBBP7, and WDR5) were identified by the iUUCD 2.0 database (Figure 4E). Next, the substrates of HERC2 were predicted (Figure 4F), and p53 was the only one predicted in the top 20 protein list (Figure 4G). The interacting domains of HERC2 with p53 are shown in Figure 4H. Taken together, we propose that BTF3 may also play an oncogenic role in CRC via HERC2-mediated p53 ubiquitination and degradation.

FIGURE 4

Figure 4. BTF3 may regulate HERC2-mediated p53 ubiquitination in CRC. (A) Immunoprecipitation-MS was conducted and revealed 542 BTF3 specific proteins. The biological processes (B), cellular component (C), and molecular function (D) of 542 proteins were analyzed, and the top 10 items are shown. (E) E3 ligases in Ubibrowser and iUUCD databases were obtained, and HERC2 was enriched in all three lists. (F) Substrates of HERC2 were predicted with Ubibrowser. (G) Venn diagram of the top 20 predicted HERC2 substrates with 542 BTF3 binding proteins. (H) Predicted binding domain of BERC2 with p53.

MiR-497 Decreases BTF3 Expression in CRC

A previous study indicated that BTF3 was downregulated by miR-802 in ovarian cancer (Wu et al., 2020). BTF3 expression is significantly increased in CRC, while the role of miRNAs remains obscure. Therefore, we first predicted highly conserved miRNAs using miRwalk, Targetscan, and miRDB (Supplementary Table 4), and as shown in Figures 5A,B, we acquired six miRNA candidates (miR-15a-5p, miR-219a-2-3p, miR-361-5p, miR-497-5p, miR-503-5p, and miR-6838-5p). Then, their expression was analyzed in CRC samples of GEO datasets (GSE128446, GSE81581, and GSE35982). Only one miRNA, miR-497-5p, showed consistent significantly decreased expression in all three datasets (Figures 5C–E). In the next step, we transfected miR-497-5p mimics in HT29 cells, and BTF3 showed significantly decreased expression (Figure 5F). In addition, correlation analysis showed a significantly negative coefficient between miR-497-5p and BTF3 (Figure 5G). Next, the predicted binding sequence of BTF3 3′UTR was mutated (Figure 5H), and a dual luciferase assay was conducted. Indeed, luciferase activity in the wild-type (wt) group was decreased, while that in the mutation group was not significant (Figure 5I). Furthermore, the correlation analysis revealed a significant correlation of miR-497-5p with lymphatic invasion, pathologic M, pathologic N, and pathological stage (Table 1). In summary, these data suggest that elevated expression of BTF3 is caused by decreased expression of miR-497-5p in CRC.

FIGURE 5

Figure 5. MiR-497-5p is a negative regulator of BTF3 in CRC. (A) MiRNAs were predicted with miRwalk, miRDB, and Targetscan with conservation analysis, and six miRNAs were acquired. (B) The binding sequence was highly conserved in multiple species. MiR-497-5p was significantly downregulated in CRC samples of GEO datasets, GSE128446 (C), GSE81581 (D), and GSE35982 (E). (F) qRT-PCR results of BTF3 expression in miR-497-5p mimics-transfected HT29 cells. (G) Significantly negative correlation between miR-497-5p and BTF3 expression in TCGA COAD and READ datasets. (H) Predicted binding sequence of miR-497-5p was mutated. (I) Luciferase activity in the wild-type group was decreased, while that in the mutation group was not significant. *P < 0.05; N.S., not significant.

Discussion

Transcriptional and post-transcriptional gene and pathway dysregulation are important cancer signatures in CRC (Porcellini et al., 2018). BTF3, namely, basic transcription factor 3, has been proven to be an oncogene and hazardous prognosticator in CRC (Liu et al., 2019b). As a transcription factor, BTF3 is believed to regulate gene expression by binding target promoters, and this mechanism has been reported in breast cancer (Porcellini et al., 2018). Here, we refer to this transcription factor characteristic of BTF3 as the transcriptional mechanism, a molecular function primarily occurring in the nucleus. In CRC, the function of BTF3 was first reported to be linked to MAD2L2, MCM3, and PLK1 (Liu et al., 2019b). In addition, the alias name of BTF3 is nascent-polypeptide-associated complex beta (NACB). NAC is highly conserved from yeast to human and is believed to bind nascent peptides from ribosomes. According to S. Rospert et al., multiple hypotheses of NAC molecular functions have been proposed, including transcriptional coactivators and regulators of peptide translocation to ER and mitochondria (Rospert et al., 2002). Recent studies demonstrated that NAC complexes are involved in protein ubiquitination (Panasenko et al., 2009) and proteolysis (Kirstein-Miles et al., 2013), a molecular function primarily occurring in the cytoplasm. As a subunit of NAC, BTF3 has been reported to inhibit ubiquitin-mediated BMI1 degradation (Hu et al., 2019). Based on the BTF3 IP result, cytoplasmic and nuclear proteins were both identified (Figure 4C and Table 1), indicating that both mechanisms may exist for BTF3 in CRC, so we conducted a study on transcriptional and NAC-related mechanism.

To investigate the transcriptional mechanism of BTF3, RNA sequencing and ChIP sequencing were applied, and we identified CHD1L as a target of BTF3. CHD1L is an oncogene in several cancer types, particularly in hepatocellular carcinoma (Chen et al., 2009, 2010, 2011; Li et al., 2013; Liu et al., 2016). CHD1L is also known as amplified in liver cancer protein 1, a chromatin remodeling protein induced in response to DNA damage and a transcriptional regulator. In CRC, CHD1L was reported to promote tumor progression and to predict poor survival (Liu et al., 2019b). Here, we found that CHD1L itself is transcriptionally activated by BTF3. In addition, we observed that, in BTF3 binding proteins, other transcription factors and epigenetic regulators, such as TCF7L2, CHD4, and HDAC2, were identified. Based on the results of this study, we could not rule out the possibility that BTF3 may regulate target gene expression as a transcriptional coactivator by interacting with other transcription factors and epigenetic regulators. In the future, we will explore this possibility in additional studies.

Another finding in this study is the possible mechanism of BTF3 as a regulator of HERC2-mediated p53 ubiquitination and proteolysis. P53 is a well-known tumor suppressor that is regulated by the ubiquitin E3 ligase MDM2 in CRC (Liu et al., 2017; Choi et al., 2018; Zeng et al., 2018). HERC2 is an E3 ligase involved in tumor progression through degradation of cancer-specific substrates (Wu et al., 2010; Izawa et al., 2011; Zhu et al., 2014). Previous studies demonstrated that HERC2 is involved in regulating p53 oligomerization and the MDM2–p53 pathway (Kaustov et al., 2007; Cubillos-Rojas et al., 2014; Garcia-Cano et al., 2020). However, there are still limited studies concerning the role of HERC2 in CRC (Yoo et al., 2011). Here, we found that BTF3 binds HERC2 and p53, and considering the function of BTF3 as an oncogene and BTF3 as a tumor suppressor, we speculated that BTF3 might promote p53 ubiquitination and degradation by recruiting HERC2. To further explore this mechanism, more evidence, including functional analysis of HERC2 and p53 with BTF3, the impact of BTF3 and HERC2 p53 ubiquitination and expression level, co-IP of HERC2 and p53 with BTF3, transcriptome RNA sequencing and ChIP sequencing comparison of p53 and BTF3, etc., would be conducted in a subsequent study.

MicroRNAs (miRNAs) are endogenous single-stranded RNAs of 18–25 nt that are well recognized as post-transcriptional gene expression regulators in CRC (Xuan et al., 2015; Moridikia et al., 2018). In this study, we first predicted six highly conserved miRNAs for BTF3: miR-15a-5p, miR-219a-2-3p, miR-361-5p, miR-497-5p, miR-503-5p, and miR-6838-5p. MiR-15a-5p showed significantly increased expression and was a hazardous prognostic marker in CRC linked to recurrence (Kontos et al., 2017; Al Mustanjid et al., 2020). For miR-219a, miR-219a-5p has been reported to inhibit CRC progression by repressing oncogenic targets (Cheng et al., 2015; Xiong et al., 2015; Huang et al., 2017; Wang et al., 2017), while the function of miR-219a-2-3p in CRC is unclear. Similarly, miR-497-5p also exhibits a tumor-suppressing role in CRC (Wang et al., 2019; Bai et al., 2020; Gharib et al., 2020; Yu and Zhang, 2020). In this study, by expression validation and luciferase assay, we concluded that elevated expression of BTF3 in CRC was mediated by miR-497-5p, whose expression is significantly decreased in CRC.

Collectively, we conclude that BTF3 is an oncogene in CRC, which may operate using two distinct functional mechanisms. Acting as a transcription factor, BTF3 promotes the transcription and expression of CHD1L. In addition, BTF3 may also inhibit E3 ubiquitin ligase HERC2-mediated p53 degradation. Decreased miR-497-5p expression is responsible for higher levels of BTF3 post-transcriptionally (Figure 6).

FIGURE 6

Figure 6. Schematic diagram of BTF3 functional and expression mechanism in CRC. BTF3 is an oncogene in CRC, which may operate using two distinct functional mechanisms. Acting as a transcription factor, BTF3 promotes the transcription and expression of CHD1L, and on the other hand, BTF3 may also inhibit E3 ubiquitin ligase HERC2-mediated p53 degradation. Decreased miR-497-5p expression is responsible for higher levels of BTF3 post-transcriptionally.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Ethics Statement

The animal study was reviewed and approved by Ethics Committees of Changhai Hospital.

Author Contributions

HW and JX conceived and designed the experiments. WW performed visualization. GL and HH analyzed the data. MS, YL, and HC prepared figures and/or tables. XL drafted the manuscript. All authors read and approved the final manuscript.

Funding

This study was supported by grants from the National Natural Science Foundation of China (Nos. 81700717 and 81702297).

Conflict of Interest

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

Supplementary Material

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

Footnotes

1. ^ http://www.genome.jp/kegg/
2. ^ https://string-db.org/
3. ^ http://gepia.cancer-pku.cn/
4. ^ https://cistrome.shinyapps.io/timer/
5. ^ http://genome.ucsc.edu
6. ^ http://ubibrowser.nb.org

References

Al Mustanjid, M., Mahmud, S. M. H., Royel, M. R. I., Rahman, M. H., Islam, T., Rahman, M. R., et al. (2020). Detection of molecular signatures and pathways shared in inflammatory bowel disease and colorectal cancer: a bioinformatics and systems biology approach. Genomics 112, 3416–3426. doi: 10.1016/j.ygeno.2020.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Bai, J., Xu, J., Zhao, J., and Zhang, R. (2020). lncRNA SNHG1 cooperated with miR-497/miR-195-5p to modify epithelial-mesenchymal transition underlying colorectal cancer exacerbation. J. Cell. Physiol. 235, 1453–1468. doi: 10.1002/jcp.29065

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 68, 394–424. doi: 10.3322/caac.21492

PubMed Abstract | CrossRef Full Text | Google Scholar

Cavallini, B., Huet, J., Plassat, J. L., Sentenac, A., Egly, J. M., and Chambon, P. (1988). A yeast activity can substitute for the HeLa cell TATA box factor. Nature 334, 77–80. doi: 10.1038/334077a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Chan, T. H., Yuan, Y. F., Hu, L., Huang, J., Ma, S., et al. (2010). CHD1L promotes hepatocellular carcinoma progression and metastasis in mice and is associated with these processes in human patients. J. Clin. Invest. 120, 1178–1191. doi: 10.1172/jci40665

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Hu, L., Chan, T. H., Tsao, G. S., Xie, D., Huo, K. K., et al. (2009). Chromodomain helicase/adenosine triphosphatase DNA binding protein 1-like (CHD1l) gene suppresses the nucleus-to-mitochondria translocation of nur77 to sustain hepatocellular carcinoma cell survival. Hepatology 50, 122–129. doi: 10.1002/hep.22933

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Yuan, Y. F., Li, Y., Chan, T. H., Zheng, B. J., Huang, J., et al. (2011). Clinical significance of CHD1L in hepatocellular carcinoma and therapeutic potentials of virus-mediated CHD1L depletion. Gut 60, 534–543. doi: 10.1136/gut.2010.224071

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, J., Deng, R., Zhang, P., Wu, C., Wu, K., Shi, L., et al. (2015). miR-219-5p plays a tumor suppressive role in colon cancer by targeting oncogene Sall4. Oncol. Rep. 34, 1923–1932. doi: 10.3892/or.2015.4168

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, T. G., Nguyen, M. N., Kim, J., Jo, Y. H., Jang, M., Nguyen, N. N. Y., et al. (2018). Cyclophilin B induces chemoresistance by degrading wild-type p53 via interaction with MDM2 in colorectal cancer. J. Pathol. 246, 115–126. doi: 10.1002/path.5107

PubMed Abstract | CrossRef Full Text | Google Scholar

Cubillos-Rojas, M., Amair-Pinedo, F., Peiró-Jordán, R., Bartrons, R., Ventura, F., and Rosa, J. L. (2014). The E3 ubiquitin protein ligase HERC2 modulates the activity of tumor protein p53 by regulating its oligomerization. J. Biol. Chem. 289, 14782–14795. doi: 10.1074/jbc.m113.527978

PubMed Abstract | CrossRef Full Text | Google Scholar

Diboun, I., Wernisch, L., Orengo, C. A., and Koltzenburg, M. (2006). Microarray analysis after RNA amplification can detect pronounced differences in gene expression using limma. BMC Genomics 7:252. doi: 10.1186/1471-2164-7-252

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Wang, X., Zhang, Y., Sang, X., Yi, J., Liu, C., et al. (2019). Inhibition of BTF3 sensitizes luminal breast cancer cells to PI3Kalpha inhibition through the transcriptional regulation of ERalpha. Cancer Lett. 440-441, 54–63. doi: 10.1016/j.canlet.2018.09.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, J., Tang, W., Du, P., Wang, G., Chen, W., Li, J., et al. (2012). Identifying microRNA-mRNA regulatory network in colorectal cancer by a combination of expression profile and bioinformatics analysis. BMC Syst. Biol. 6:68. doi: 10.1186/1752-0509-6-68

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, T., Liu, Z., Wang, Y., Cheng, H., Yang, Q., Guo, A., et al. (2013). UUCD: a family-based database of ubiquitin and ubiquitin-like conjugation. Nucleic Acids Res. 41, D445–D451.

Google Scholar

Garcia-Cano, J., Sánchez-Tena, S., Sala-Gaston, J., Figueras, A., Viñals, F., Bartrons, R., et al. (2020). Regulation of the MDM2-p53 pathway by the ubiquitin ligase HERC2. Mol. Oncol. 14, 69–86. doi: 10.1002/1878-0261.12592

PubMed Abstract | CrossRef Full Text | Google Scholar

Gharib, E., Nasri Nasrabadi, P., and Reza Zali, M. (2020). miR-497-5p mediates starvation-induced death in colon cancer cells by targeting acyl-CoA synthetase-5 and modulation of lipid metabolism. J. Cell. Physiol. 235, 5570–5589. doi: 10.1002/jcp.29488

PubMed Abstract | CrossRef Full Text | Google Scholar

Green, C. D., Thompson, P. D., Johnston, P. G., and El-Tanani, M. K. (2007). Interaction between transcription factor, basal transcription factor 3, and the NH2-terminal domain of human estrogen receptor alpha. Mol. Cancer Res. 5, 1191–1200. doi: 10.1158/1541-7786.mcr-07-0123

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Sun, F., Chen, W., Zhang, J., Zhang, T., Qi, M., et al. (2019). BTF3 sustains cancer stem-like phenotype of prostate cancer via stabilization of BMI1. J. Exp. Clin. Cancer Res. 38:227.

Google Scholar

Huang, L. X., Hu, C. Y., Jing, L., Wang, M. C., Xu, M., Wang, J., et al. (2017). microRNA-219-5p inhibits epithelial-mesenchymal transition and metastasis of colorectal cancer by targeting lymphoid enhancer-binding factor 1. Cancer Sci. 108, 1985–1995. doi: 10.1111/cas.13338

PubMed Abstract | CrossRef Full Text | Google Scholar

Izawa, N., Wu, W., Sato, K., Nishikawa, H., Kato, A., Boku, N., et al. (2011). HERC2 Interacts with Claspin and regulates DNA origin firing and replication fork progression. Cancer Res. 71, 5621–5625. doi: 10.1158/0008-5472.can-11-0385

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, X., Li, J., Zhu, L., Cai, J., Zhang, J., Qu, Y., et al. (2013). CHD1L promotes tumor progression and predicts survival in colorectal carcinoma. J. Surg. Res. 185, 84–91. doi: 10.1016/j.jss.2013.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, G., Hernández-Illán, E., Moreira, L., Balaguer, F., and Goel, A. (2020). Epigenetics of colorectal cancer: biomarker and therapeutic potential. Nat. Rev. Gastroenterol. Hepatol. 17, 111–130. doi: 10.1038/s41575-019-0230-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanno, M., Chalut, C., and Egly, J. M. (1992). Genomic structure of the putative BTF3 transcription factor. Gene 117, 219–228. doi: 10.1016/0378-1119(92)90732-5

CrossRef Full Text | Google Scholar

Kaustov, L., Lukin, J., Lemak, A., Duan, S., Ho, M., Doherty, R., et al. (2007). The conserved CPH domains of Cul7 and PARC are protein-protein interaction modules that bind the tetramerization domain of p53. J. Biol. Chem. 282, 11300–11307. doi: 10.1074/jbc.m611297200

PubMed Abstract | CrossRef Full Text | Google Scholar

Keum, N., and Giovannucci, E. (2019). Global burden of colorectal cancer: emerging trends, risk factors and prevention strategies. Nat. Rev. Gastroenterol. Hepatol. 16, 713–732. doi: 10.1038/s41575-019-0189-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirstein-Miles, J., Scior, A., Deuerling, E., and Morimoto, R. I. (2013). The nascent polypeptide-associated complex is a key regulator of proteostasis. EMBO J. 32, 1451–1468. doi: 10.1038/emboj.2013.87

PubMed Abstract | CrossRef Full Text | Google Scholar

Kontos, C. K., Tsiakanikas, P., Avgeris, M., Papadopoulos, I. N., and Scorilas, A. (2017). miR-15a-5p, A Novel Prognostic Biomarker, Predicting Recurrent Colorectal Adenocarcinoma. Mol. Diagn. Ther. 21, 453–464. doi: 10.1007/s40291-017-0270-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Koplin, A., Preissler, S., Ilina, Y., Koch, M., Scior, A., Erhardt, M., et al. (2010). A dual function for chaperones SSB-RAC and the NAC nascent polypeptide-associated complex on ribosomes. J. Cell Biol. 189, 57–68. doi: 10.1083/jcb.200910074

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, W90–W97.

Google Scholar

Kusumawidjaja, G., Kayed, H., Giese, N., Bauer, A., Erkan, M., Giese, T., et al. (2007). Basic transcription factor 3 (BTF3) regulates transcription of tumor-associated genes in pancreatic cancer cells. Cancer Biol. Ther. 6, 367–376. doi: 10.4161/cbt.6.3.3704

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77, e108–e110.

Google Scholar

Li, X., Sui, J., Xing, J., Cao, F., Wang, H., Fu, C., et al. (2019). Basic transcription factor 3 expression silencing attenuates colon cancer cell proliferation and migration in vitro. Oncol. Lett. 17, 113–118.

Google Scholar

Li, Y., Chen, L., Chan, T. H., Liu, M., Kong, K. L., Qiu, J. L., et al. (2013). SPOCK1 is regulated by CHD1L and blocks apoptosis and promotes HCC cell invasiveness and metastasis in mice. Gastroenterology 144, 179.e4–191.e4.

Google Scholar

Li, Y., Xie, P., Lu, L., Wang, J., Diao, L., Liu, Z., et al. (2017). An integrated bioinformatics platform for investigating the human E3 ubiquitin ligase-substrate interaction network. Nat. Commun. 8:347.

Google Scholar

Liu, M., Chen, L., Ma, N. F., Chow, R. K., Li, Y., Song, Y., et al. (2016). CHD1L promotes lineage reversion of hepatocellular carcinoma through opening chromatin for key developmental transcription factors. Hepatology 63, 1544–1559. doi: 10.1002/hep.28437

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Jiang, L., Wang, W., and Jiang, T. (2019a). BTF3 silencing inhibits the proliferation of osteosarcoma cells. J. Cancer 10, 1855–1861. doi: 10.7150/jca.28476

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Q., Wu, J., Lu, T., Fang, Z., Huang, Z., Lu, S., et al. (2019b). Positive expression of basic transcription factor 3 predicts poor survival of colorectal cancer patients: possible mechanisms involved. Cell Death Dis. 10:509.

Google Scholar

Liu, Q., Zhou, J. P., Li, B., Huang, Z. C., Dong, H. Y., Li, G. Y., et al. (2013). Basic transcription factor 3 is involved in gastric cancer development and progression. World J. Gastroenterol. 19, 4495–4503. doi: 10.3748/wjg.v19.i28.4495

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Tackmann, N. R., Yang, J., and Zhang, Y. (2017). Disruption of the RP-MDM2-p53 pathway accelerates APC loss-induced colorectal tumorigenesis. Oncogene 36, 1374–1383. doi: 10.1038/onc.2016.301

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, K. D., Siegel, R. L., Lin, C. C., Mariotto, A. B., Kramer, J. L., Rowland, J. H., et al. (2016). Cancer treatment and survivorship statistics, 2016. CA Cancer J. Clin. 66, 271–289.

Google Scholar

Moridikia, A., Mirzaei, H., Sahebkar, A., and Salimian, J. (2018). MicroRNAs: potential candidates for diagnosis and treatment of colorectal cancer. J. Cell. Physiol. 233, 901–913. doi: 10.1002/jcp.25801

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, L. H., Goel, A., and Chung, D. C. (2020). Pathways of colorectal carcinogenesis. Gastroenterology 158, 291–302. doi: 10.1053/j.gastro.2019.08.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Panasenko, O. O., David, F. P., and Collart, M. A. (2009). Ribosome association and stability of the nascent polypeptide-associated complex is dependent upon its own ubiquitination. Genetics 181, 447–460. doi: 10.1534/genetics.108.095422

PubMed Abstract | CrossRef Full Text | Google Scholar

Parvin, J. D., Shykind, B. M., Meyers, R. E., Kim, J., and Sharp, P. A. (1994). Multiple sets of basal factors initiate transcription by RNA polymerase II. J. Biol. Chem. 269, 18414–18421.

Google Scholar

Porcellini, E., Laprovitera, N., Riefolo, M., Ravaioli, M., Garajova, I., and Ferracin, M. (2018). Epigenetic and epitranscriptomic changes in colorectal cancer: diagnostic, prognostic, and treatment implications. Cancer Lett. 419, 84–95. doi: 10.1016/j.canlet.2018.01.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Rospert, S., Dubaquie, Y., and Gautschi, M. (2002). Nascent-polypeptide-associated complex. Cell Mol. Life. Sci. 59, 1632–1639. doi: 10.1007/pl00012490

PubMed Abstract | CrossRef Full Text | Google Scholar

Sayagues, J. M., Corchete, L. A., Gutiérrez, M. L., Sarasquete, M. E., Del Mar Abad, M., Bengoechea, O., et al. (2016). Genomic characterization of liver metastases from colorectal cancer patients. Oncotarget 7, 72908–72922.

Google Scholar

Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45, W98–W102.

Google Scholar

Wang, H., Yu, M., Hu, W., Chen, X., Luo, Y., Lin, X., et al. (2019). Linc00662 promotes tumorigenesis and progression by regulating miR-497-5p/AVL9 axis in colorectal cancer. Front Genet. 10:1385. doi: 10.3389/fgene.2019.01385

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Zhu, L., Jiang, Y., Xu, J., Wang, F., and He, Z. (2017). miR-219-5p suppresses the proliferation and invasion of colorectal cancer cells by targeting calcyphosin. Oncol. Lett. 13, 1319–1324. doi: 10.3892/ol.2017.5570

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W., Sato, K., Koike, A., Nishikawa, H., Koizumi, H., Venkitaraman, A. R., et al. (2010). HERC2 is an E3 ligase that targets BRCA1 for degradation. Cancer Res. 70, 6384–6392. doi: 10.1158/0008-5472.can-10-1304

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, X., Liu, L., and Zhang, H. (2020). miR802 inhibits the epithelialmesenchymal transition, migration and invasion of cervical cancer by regulating BTF3. Mol. Med. Rep. 22, 1883–1891. doi: 10.3892/mmr.2020.11267

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, G. B., Zhang, G. N., Xiao, Y., Niu, B. Z., Qiu, H. Z., Wu, B., et al. (2015). MicroRNA-219-5p functions as a tumor suppressor partially by targeting platelet-derived growth factor receptor alpha in colorectal cancer. Neoplasma 62, 855–863. doi: 10.4149/neo_2015_104

CrossRef Full Text | Google Scholar

Xuan, Y., Yang, H., Zhao, L., Lau, W. B., Lau, B., Ren, N., et al. (2015). MicroRNAs in colorectal cancer: small molecules with big functions. Cancer Lett. 360, 89–105. doi: 10.1016/j.canlet.2014.11.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoo, N. J., Park, S. W., and Lee, S. H. (2011). Frameshift mutations of ubiquitination-related genes HERC2, HERC3, TRIP12, UBE2Q1 and UBE4B in gastric and colorectal carcinomas with microsatellite instability. Pathology 43, 753–755. doi: 10.1097/pat.0b013e32834c7e78

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, C., and Zhang, F. (2020). LncRNA AC009022.1 enhances colorectal cancer cells proliferation, migration, and invasion by promoting ACTR3B expression via suppressing miR-497-5p. J. Cell Biochem. 121, 1934–1944. doi: 10.1002/jcb.29428

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, K., Chen, X., Hu, X., Liu, X., Xu, T., Sun, H., et al. (2018). LACTB, a novel epigenetic silenced tumor suppressor, inhibits colorectal cancer progression by attenuating MDM2-mediated p53 ubiquitination and degradation. Oncogene 37, 5534–5551. doi: 10.1038/s41388-018-0352-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D. Z., Chen, B. H., Zhang, L. F., Cheng, M. K., Fang, X. J., and Wu, X. J. (2017). Basic transcription factor 3 is required for proliferation and epithelial-mesenchymal transition via regulation of FOXM1 and JAK2/STAT3 signaling in gastric cancer. Oncol. Res. 25, 1453–1462. doi: 10.3727/096504017x14886494526344

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Gross, N., Li, Z., Yin, G., Zhong, Q., Liu, C., et al. (2019). Upregulation of BTF3 affects the proliferation, apoptosis, and cell cycle regulation in hypopharyngeal squamous cell carcinoma. Biomed. Pharmacother. 118:109211. doi: 10.1016/j.biopha.2019.109211

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., Xu, Y., Lin, S., Guo, Y., Deng, W., Zhang, Y., et al. (2018). iUUCD 2.0: an update with rich annotations for ubiquitin and ubiquitin-like conjugations. Nucleic Acids Res. 46, D447–D453.

Google Scholar

Zhu, M., Zhao, H., Liao, J., and Xu, X. (2014). HERC2/USP20 coordinates CHK1 activation by modulating CLASPIN stability. Nucleic Acids Res. 42, 13074–13081. doi: 10.1093/nar/gku978

PubMed Abstract | CrossRef Full Text | Google Scholar