- Research
- Open access
- Published:
The epigenetic downregulation of LncGHRLOS mediated by RNA m6A methylase ZCCHC4 promotes colorectal cancer tumorigenesis
Journal of Experimental & Clinical Cancer Research volume 43, Article number: 44 (2024)
Abstract
Background
m6A modification is currently recognized as a major driver of RNA function that maintains cancer cell homeostasis. Long non-coding (Lnc) RNAs control cell proliferation and play an important role in the occurrence and progression of colorectal cancer (CRC). ZCCHC4 is a newly discovered m6A methyltransferase whose role and mechanism in tumors have not yet been elucidated.
Methods
The EpiQuik m6A RNA methylation kit was used to detect the level of total RNA m6A in six types of digestive tract tumors. The Kaplan-Meier method and receiver operating characteristic curve were used to evaluate the prognostic and diagnostic value of the newly discovered m6A methyltransferase, ZCCHC4, in CRC. The effects on CRC growth in vitro and in vivo were studied using gain- and loss-of-function experiments. The epigenetic mechanisms underlying ZCCHC4 upregulation in CRC were studied using RIP, MeRIP-seq, RNA pull-down, and animal experiments.
Results
We reported that the ZCCHC4-LncRNAGHRLOS-KDM5D axis regulates the growth of CRC in vitro and in vivo. We found that ZCCHC4 was upregulated in primary CRC samples and could predict adverse clinical outcomes in patients with CRC. Mechanistically, ZCCHC4 downregulated LncRNAGHRLOS to promote CRC tumorigenesis. As a downstream molecule of LncRNAGHRLOS, KDM5D directly controls CRC cell proliferation, migration, and invasion.
Conclusion
This study suggests that the ZCCHC4 axis contributes to the tumorigenesis and progression of CRC and that ZCCHC4 may be a potential biomarker for this malignancy.
Background
Colorectal cancer (CRC) is the second most common malignant tumor of the digestive system worldwide, with high incidence and mortality rates. In China, approximately 376,000 new cases and 191,000 deaths occur annually, and the incidence and mortality rates keep rising [1, 2]. However, most patients show no apparent clinical symptoms in the early stages of the disease due to the lack of early detection indicators or characteristic molecular markers. Thus, most patients are usually diagnosed in the late stages and lose the opportunity for surgery or other treatments [3, 4].
In the 1980s, scientists found that when the nuclear DNA sequence did not change, gene function could undergo reversible and heritable changes [5, 6]. This field is referred to as epigenetics. DNA methylation is a well-known method for epigenetic modification. Similarly, there are many chemical modifications on eukaryotic RNA, including N6 methyladenine (m6A), N7 methylguanine (m7G), and N1 methyladenine (m1A), among which m6A is the most abundant chemical modification. Regarding m6A, three important components have been identified: methyltransferase, demethylase, and m6A binding protein. Methyltransferase is figuratively called “Writer,” and demethylase is called “Eraser.” The m6A binding protein recognizes m6A and is called the “Reader.” m6A is a dynamically regulated reversible chemical modification. Some studies have shown that once primary stimuli such as hypoxia occur, it leads to the uncontrolled expression of these genes, leading to angiogenesis and tumor growth progress [7,8,9,10].
In 2012, the m6A modification map of the entire human transcriptome was first revealed using methylated RNA immunoprecipitation combined with high-throughput sequencing (MeRIP-seq) technology [11, 12]. Since then, RNA m6A modification has attracted increasing attention from scientists. In addition to mRNA, m6A modification occurs in various non-coding RNA (ncRNAs), such as tRNA, rRNA, snRNA, snoRNA, and long non-coding RNA (lncRNA). For example, m6A modification occurs on LncRNA-MALAT1, which is closely related to dynamic atherosclerosis and is considered to be closely related to changes in the stem ring structure [13,14,15,16,17,18].
The role of m6A in various diseases, especially cancer, has attracted increasing attention. Many studies have confirmed that changes in the m6A modification intensity widely affect the progression of various diseases, and research in colorectal cancer is becoming increasingly in-depth [19,20,21,22,23]. The comprehensive extent of m6A modification is under dynamic regulation, with both its “writer” and “eraser” contributing to the regulatory process. Although the roles of methyltransferases in CRC have been widely reported, the newly discovered involvement of methyltransferases in CRC progression will continue to enhance our understanding of m6A. In this study, we identified the role of ZCCHC4, a key methyltransferase, in promoting the proliferation, migration, and metastasis of CRC cells. We conducted additional screening of essential downstream molecules and unveiled the mechanism by which ZCCHC4 modulates the behavior of CRC cells.
Materials and methods
CRC samples from patients and cell lines
Tumors and adjacent tissues of hepatocellular carcinoma (HCC), esophageal cancer (EC), gastric cancer (GC), gallbladder cancer, CRC, and pancreatic cancer (PC) were procured from patients at the Department of General Surgery, The First Affiliated Hospital of Anhui Medical University. The samples were promptly stored in liquid nitrogen for preservation, and informed consent was obtained from each participant. Human CRC cell lines were acquired from the National Collection of Authenticated Cell Cultures and cultured in Dulbecco’s Modified Eagle Medium (DMEM) containing 10% Fetal Bovine Serum (FBS; Gibco, Waltham, MA, USA) and 1% penicillin-streptomycin (Beyotime, Shanghai, China).
Constructs and transfections
Genepharm (Shanghai, China) designed and synthesized short hairpin RNA (shRNA) interfering with the expression of ZCCHC4, LncGHRLOS, and KDM5D. Overexpression of ZCCHC4, lncGHRLOS, and KDM5D was performed using expression plasmids, and empty vectors were used as negative controls. The cells were transfected with the virus in the presence of a polycoagulant. After 48 h, puromycin was added to the medium to identify stable clones. The efficacy of overexpression was assessed using qRT-PCR and western blot analyses. The transfection efficiency was evaluated by qRT-PCR and western blotting after 48 h.
The m6A RNA methylation assay
Total RNA was extracted according to the manufacturer’s instructions. The m6A RNA methylation quantitative kit (Abcam, Cambridge, UK) was used to determine the RNA m6A levels in total RNA. First, the RNA packets were incubated in wells for 90 min at 37 °C. Capture antibody, detection antibody, and booster solution were added. Finally, a color-developing solution was added to detect the signal and measure absorbance. The calorimetric m6A levels were quantified and calculated using a standard curve.
Cell cycle distribution
Cell cycle analysis was performed using the Muse Cell Cycle Kit and Muse Cell Analyzer (Millipore, Bedford, MA, USA). HCT-15 cells (5 × 105 cells/well) were treated with betulin (0–8 µM) for 24 h in 6-well plates. Propidium iodide (PI)-positive cells were stained and analyzed according to the manufacturer’s instructions.
Enzyme-linked immunosorbent assay (ELISA)
All blood samples were collected using anticoagulant tubes for further experiments. Each sample contained 5 ml of blood which was centrifuged at 3500 rpm for 10 min. The supernatant plasma was manually separated and stored at -80 ℃ for batch analysis. The concentration of KDM5D was determined at the Central Laboratory of the First Affiliated Hospital of Anhui Medical University. All samples were collected in batches using commercially available kits in accordance with the manufacturer’s instructions.
Immunohistochemistry (IHC)
Human CRC sections were obtained from the Department of Pathology of the Second Affiliated Hospital of Anhui Medical University. After dewaxing, the sections were repaired in citrate buffer at 95-100°C for 10 min, and incubated with antibodies against ZCCHC4(1:500) and KDM5D (1:400). A vertical microscope (Nikon, Tokyo, Japan) was used for image acquisition. The integrated optical density (IOD) of the stained slices was evaluated using the Image-Pro Plus software (Media Cybernetics, Rockville, MD, USA).
Western blot
PRO-PREP™ Protein Extraction Solution (iNtRon Biotech, Seoul, Korea) was used to extract total protein from both cells and tissues. The lysates were subsequently combined with a 5 × sample buffer volume and separated by gel electrophoresis to isolate total proteins. The specific target proteins were detected using antibodies obtained from Abcam. Secondary antibodies were used to bind to specific primary antibodies, and the FluorChem M System (ProteinSimple, San Jose, CA, USA) was used for visualization.
Animal models
Five-week-old male BALB/c nude mice aged 5 weeks were obtained from SPF Beijing Biotechnology Co. Ltd. (Beijing, China). After random assignment, HCT-116 cells suspended in 100 µl phosphate-buffered saline (PBS) were injected subcutaneously into nude mice. All nude mice were fed in an SPF environment and given sufficient water and food, and the tumor growth rate was detected. After 21 days, eight nude mice in the ZCCHC4 group were monitored for tumor progression using molecular imaging software (the other two were not monitored using live imaging as controls). Luciferase signal strength was maintained at the same scale. After in vivo imaging, the mice were killed, the tumors were photographed, weighed, and the tumor tissues were collected for immunohistochemical analysis. All procedures are carried out in accordance with the laboratory animal management code. LncGHRLOS and KDM5D group directly killed the mice, and the tumors were photographed and weighed and stained by IHC. In addition, 5-week nude mice in each group were used to monitor the survival conditions with tumors, and the death time of mice was recorded and the survival curve was drawn.
RNA extraction and quantitative real-time PCR (qRT-PCR)
The FastPure Cell/Tissue Total RNA Isolation Kit V2 (Vazyme Biotech Co., Ltd. Nanjing, China) was used, and reverse transcription and qRT-PCR for analyzing mRNA expression were conducted using the PrimeScriptTM RT reagent kit (TAKARA, Beijing, China) and TB GREEN SuperMix (TAKARA). The experiment was performed in triplicate.
RNA-seq, methylated RNA immunoprecipitation sequencing (MeRIP-seq)
Briefly, RNA was extracted and purified to deplete ribosomal RNA and prevent DNA contamination. After fragmentation and denaturation, RNA was sheared into approximately 100-nt fragments and then incubated with an anti-m6A antibody (Abcam) together with protein A/G magnetic beads (Thermo Fisher Scientific, Waltham, MA, USA) in an immunoprecipitation buffer at 4℃ overnight. An RNase inhibitor was added. The antibody combined with methylated RNA was eluted with m6A and purified for MeRIP sequencing using BioAcme (Wuhan, China).
RNA immunoprecipitation (RIP)
A Magna RIP RNA-binding protein immunoprecipitation kit (Millipore) was used. Briefly, protein A/G magnetic beads conjugated with rabbit immunoglobulin G (17–700, Millipore), ZCCHC4 (Proteintech, USA), or LncGHRLOS (Proteintech, Rosemont, IL, USA) antibodies were incubated with cell lysates supplemented with RNase inhibitor at 4 °C overnight. After washing 6 times, RNA–protein complexes were added to proteinase K buffer. Finally, the RNA was extracted using phenol–chloroform method. qPCR was performed to determine the relative interactions.
Migration and invasion assays
For the transwell assay, chambers were prepared in a 24-well culture table with 800 ul DMEM containing 10% FBS at the bottom. Then, 5 × 104 cells in 200 ul serum-free medium were added into the upper layer. After 24 h, the migrated cells in the lower layer were fixed, stained with crystal violet, and imaged. Matrigel (Corning, NY, USA) was added to the upper layer for the invasion assay.
Wound healing assay
CRC cells were plated in 24-well plates and incubated for 24–48 h until 100% growth was achieved. Then, 200 ul micropipette tips were used to generate wounds. The cells were washed twice with 500 ul PBS. Next, 500 uL serum-free medium was added to the well. Images were captured for 24 h and analyzed using ImageJ software (National Institute of Health, Bethesda, MD, USA).
RNA stability assay
The CRC cells were grown in 12-well plates. Actinomycin D (5 µg/mL, Cell Signaling Technology, Danvers, MA, USA) was used in this experiment. The cells were collected at constant time points for RNA extraction. The remaining LncRNAs were analyzed by qRT-PCR. The half-life of the mRNA was calculated using linear regression analysis and GAPDH was used for normalization.
RNA pull-down and mass spectrometry analysis
Biotin-labeled LncGHRLOS ssRNA probes were synthesized in vitro by Sangon Biotin (Shanghai, China), while the LncGHRLOS-CDS, CDS-mut1, and CDS-mut2 mRNA were first transcribed using the MEGAscript T7 Transcription Kit (Thermo Fisher Scientific), and then labelled with biotin using Pierce RNA 3′ End Desthiobiotinylation Kit (Thermo Fisher Scientific). Next, 20 pmol of biotinylated RNA and cell lysate were mixed with streptavidin agarose beads (Thermo Fisher Scientific) at 4 ℃ overnight. After six washes, streptavidin beads were collected for western blotting and mass spectrometry.
Statistical analysis
Statistical analyses were performed using GraphPad Prism 9 software (GraphPad Software, San Diego, CA, USA). Two-tailed unpaired Student’s t-test or one-way analysis of variance were used. Spearman’s correlation analysis evaluated the relationship between the expression of ZCCHC4, lncGHRLOS, and KDM5D. Survival curves of patients with CRC were analyzed using Kaplan–Meier analysis. Receiver operating characteristic (ROC) curves were generated using GraphPad Prism 9. Each experiment had a minimum of three replications.
Results
ZCCHC4 is upregulated in CRC tissues and predicted poor prognosis in CRC
To determine the extent of m6A expression across different digestive tract tumors, five pairs of clinical tumor tissues with adjacent tissue samples were randomly selected to assess the overall levels of m6A modification (Fig. 1A-F). Significantly higher levels of m6A were observed in HCC and CRC tissues than in the adjacent normal tissues. Considering the extensive research on m6A in HCC, we selected CRC for this study. Similarly, we detected some key m6A components in five CRC specimens. Finally, we confirmed significant differences in METTL3, FTO, and ZCCHC4 expression between cancerous and adjacent tissues (Fig. 1G). After the ZCCHC4 knockout, a decline in m6A expression was observed in most CRC cell lines (Fig. 1H). The role and mechanism of ZCCHC4 in CRC have become a potential research topic.
To investigate whether ZCCHC4 is abnormally expressed in CRC, IHC was used to detect the expression level of ZCCHC4 in 243 patients with CRC, and qPCR analysis was used to explore the expression of ZCCHC4 mRNA in 30 CRC tissue samples. IHC results showed that ZCCHC4 was highly expressed in the tumor tissues of 175 patients, while patients with low expression accounted for only 26.47% (Fig. 1I). The mRNA levels of ZCCHC4 were higher in most CRC tissues (Fig. 1J). Similar results were obtained after western blot analysis of the three pairs of CRC and adjacent cancer tissues (Fig. 1K).
Between January 2015 and July 2017, postoperative specimens and clinical and pathological information related to patients with CRC were fully preserved. The relationship between ZCCHC4 expression and clinicopathological variables was also explored, and its expression was found to be closely related to the TNM stage and differentiation degree (Table 1). Next, we studied the relationship between ZCCHC4 levels and the prognosis of patients with CRC. The results of Cox regression model analysis showed that patients with higher levels of ZCCHC4 had a shorter overall survival (OS) compared with those with low levels of ZCCHC4 (Table 2). All the patients had complete follow-up data. After 60–77 months of follow-up, the results showed that the prognosis of patients with CRC with high ZCCHC4 expression was significantly worse than that of the low expression group (P < 0.0001, Fig. 2A). These results indicated that high ZCCHC4 expression was associated with poor OS.
ZCCHC4 is associated with the biological behavior of CRC cells in vitro and in vivo
To determine the potential role of ZCCHC4 in CRC progression, we first assessed the expression of ZCCHC4 among CRC cell lines. We detected ZCCHC4 levels in several CRC cell lines and found that both the mRNA and protein levels of ZCCHC4 were elevated in most CRC cell lines when compared with the normal colorectal epithelial cell line NCM460 (Fig. 2B-C). Subsequently, stable overexpression and knockdown of ZCCHC4 were established in HCT-15 cells. Transwell cell experiments were conducted to investigate the effect of ZCCHC4 on the invasive and migratory abilities of CRC cells (Fig. 2D). The results showed that in the HCT-15 cell line, downregulation of ZCCHC4 expression significantly reduced the invasion and migration abilities of CRC cells compared to the control group. As expected, the CCK-8 assay revealed that the downregulation of ZCCHC4 markedly inhibited the proliferation of HCT-15 and HCT-116 cells (Fig. 2E-F). These results were confirmed by a wound healing assay experiment (Fig. 2G). By examining the cell cycle, it was observed that downregulation of ZCCHC4 expression led to a higher proportion of CRC cells in the G1 phase in both HCT-116 and HCT-15 cell lines (Fig. 2H-I). Given that ZCCHC4 is crucial for maintaining CRC cell proliferation, we validated the in vivo role of ZCCHC4 in CRC tumorigenesis using a nude mouse xenograft model. HCT-116 cells with stable and low ZCCHC4 expression mediated by lentiviruses were subcutaneously inoculated into nude mice to study the tumorigenicity of ZCCHC4 in CRC cells. It is worth noting that the high expression group of ZCCHC4 showed increased growth ability of CRC cells in nude mice, resulting in more volume, as confirmed by in vivo imaging (Fig. 3A-D), increased tumor weight (Fig. 3E) and poor prognosis (Fig. 3F). IHC results showed a decrease in the molecular marker Ki-67 for cell proliferation in ZCCHC4 deficient tumors compared to overexpression group, suggesting that ZCCHC4 may be a driving factor for CRC tumor development (Fig. 3G).
LncGHRLOS was identified as a downstream target of ZCCHC4
To explore the potential targets of ZCCHC4 in CRC, we conducted RNA analysis of control and CRC cell lines expressing ZCCHC4. RNA-seq results showed that most transcripts were downregulated after ZCCHC4 overexpression. A volcano map was used to visualize the overall distribution of differentially expressed genes. Two levels of difference, multiple (log2(fold change)) and significance, were evaluated (Fig. 4A). Gene expression cluster analysis was used to determine the clustering patterns of gene expression in different treatments (Fig. 4B). The protein sequences of all genes and identified differentially expressed genes (gene sequences) were compared with protein sequences from the UniProt database, and the comparison results were annotated by GO function according to the known protein GO annotation in the UniProt database and the pathway significance was tested using the hypergeometric distribution. The GO term screening condition for significant enrichment was a p-value of hypergeometric distribution test less than 0.05(Fig. 4C). The identified differentially expressed genes were enriched in the KEGG metabolic pathway using the “clusterpro digestion” package in R (Fig. 4D). In order to identify the main target genes that play a critical role downstream of ZCCHC4, based on MeRIP-seq results and combined with the genes that significantly changed after overexpression of ZCCHC4 in RNA seq results, we identified 11 common genes (Fig. 4E). After a literature search, we selected CPA4, Col5A3, lncGHRLOS, and CLEC18B as the potential downstream targets of ZCCHC4(Fig. 4F). After detecting the mRNA expression of these molecules in 10 pairs of CRC tissues and adjacent normal tissues, we found that only lncGHRLOS differed between the two groups, and we used lncGHRLOS as our candidate molecule for further verification (Fig. 4G-J). We first described the correlation between lncGHRLOS and ZCCHC4 using PCR results for the above 50 samples, and the results showed that the two presented a negative correlation (Fig. 4K-L). Moreover, after knocking down the expression of ZCCHC4, we detected significantly increased lncGHRLOS mRNA levels in the CRC cell lines (Fig. 4M). RNA pull-down and RIP experiments also showed binding between ZCCHC4 and lncGHRLOS transcripts (Fig. 4N-O). Notably, ZCCHC4 ablation slowed the degradation of lncGHRLOS in CRC cells, whereas ZCCHC4 overexpression accelerated the degradation of lncGHRLOS transcripts (Fig. 4P). Therefore, we speculate that lncGHRLOS may be a downstream effector of ZCCHC4.
LncGHRLOS inhibited the proliferation, invasion, and migration of CRC in vitro and in vivo
To verify the unique role of lncGHRLOS in CRC, we first analyzed lncGHRLOS in CRC cell lines (Fig. 5A). shRNA was used to specifically reduce the intracellular RNA levels of lncGHRLOS in HCT-15 cells and activate the expression of lncGHRLOS through a lentivirus-mediated system. The proliferative ability of CRC cells treated with lncGHRLOS was detected using a CCK-8 assay. The CCK-8 activity experiment showed that silencing lncGHRLOS effectively enhanced the proliferative ability of CRC cells, while enhanced expression of lncGHRLOS caused the opposite biological effect (Fig. 5B). A wound healing experiment was performed to determine the migration ability of the cells (Fig. 5C-D). Our results showed that after overexpression of lncGHRLOS, the migratory ability was significantly reduced compared to that of ordinary CRC cells, and there was a significant difference compared to the low expression group of lncGHRLOS. The Transwell assay yielded almost identical results (Fig. 5E).
Subsequently, we inoculated nude mice with CRC cells stably overexpressing lncGHRLOS to study the effect of lncGHRLOS on tumorigenicity. Notably, overexpression of lncGHRLOS reduced the growth of CRC cells in nude mice, resulting in a smaller tumor volume (Fig. 5F-G), reduced tumor weight (Fig. 5H), better prognosis (P = 0.045, Fig. 5I). Consistent with these observations, IHC results showed that the molecular marker Ki-67 for cell proliferation was downregulated in lncGHRLOS-overexpression tumors compared to the corresponding control group (Fig. 5J). Moreover, samples of patients with CRC with CRC used for PCR detection and ROC analysis showed that lncGHRLOS had significant differentiation in clinical CRC diagnosis (AUC = 0.7664, Fig. 5K). The survival curve of the patients was drawn based on the follow-up results. The results showed that patients in the high lncGHRLOS expression group had a longer survival period (P = 0.0300, Fig. 5L).
Correlation between ZCCHC4, LncGHRLOS and KDM5D
RNA pull-down experiments and mass spectrometry were performed to investigate the potential downstream targets of lncGHRLOS. Mass spectrometry shows that KDM5D may play a role in the m6A modification of LncGHRLOS (Fig. 6A). Using the median of PCR results, we divided patients with high and low expression of KDM5D in clinical specimens into groups and found that KDM5D was poorly expressed in CRC, and its low expression conveyed poor prognosis in patients with CRC (P = 0.009, Fig. 6B). ELISA results also showed that KDM5D has the potential to be a marker of CRC (AUC = 0.7638, Fig. 6C). CRC and adjacent tissues were selected for PCR experiments to confirm KDM5D mRNA levels (P = 0.0425, Fig. 6D). Furthermore, a significant positive correlation between LncGHRLOS and KDM5D mRNA expression was observed in CRC (P = 0.0017, Fig. 6E). Validation experiments in cell lines confirmed that the expression of KDM5D was significantly reduced after the knockdown of LncGHRLOS (Fig. 6F-H). Importantly, we evaluated the stability of lncGHRLOS mRNA after knockdown and overexpression of KDM5D. It was found that the mRNA stability of lncGHRLOS was significantly increased after overexpression of KDM5D under actinomycin D treatment but significantly decreased after knockdown of KDM5D (Fig. 6I). Significantly reduced expression of KDM5D mRNA was observed in many CRC cell lines (Fig. 6J). Next, we used a wound healing assay (Fig. 6K), CCK-8 assay (Fig. 6L), transwell assay (Fig. 6M) and flow cytometry (Fig. 6N) to verify the influence of KDM5D knockdown and overexpression on the biological behavior of CRC cells. We found that high expression of KDM5D reduced the proliferation and migration ability of HCT-116 cells. The cell cycle stagnates during the G0/G1 phase. Furthermore, we performed in vivo experiments to clarify the effects of KDM5D on CRC. Tumor formation experiments in nude mice showed that the low expression group of KDM5D showed increased growth ability of CRC cells in nude mice, resulting in more volume (Fig. 6O-P), increased tumor weight (Fig. 6Q) and Ki-67 index (Fig. 6S). However, unlike in previous experiments, there was no significant difference in survival time between the two groups of mice (P = 0.09, Fig. 6R). This may be due to a bias caused by the small number of experimental mice.
Discussion
Cell proliferation, migration, and metastasis are the fundamental differences between tumor and normal cells. Much has been published regarding RNA modifications and disease-related mechanisms in CRC. A deeper understanding of the role of RNA modifications in CRC progression is essential for the search for novel tumor markers and the evolution of therapeutic strategies [24–27]. m6A is the most prevalent type of RNA modification among mRNAs and ncRNAs. The “writer” and “eraser” form a dynamic regulatory process to control the level of m6A in the human body. Several studies suggest that m6A regulators play important roles in various physiological and pathological diseases by regulating the epigenetic transcriptional levels of genes. m6A modification is involved in various human diseases, such as cardiovascular, digestive, reproductive, nervous, and blood system diseases. Although many studies have investigated the role of key m6A components in CRC [28,29,30,31,32,33], as a newly discovered m6A methylase, the role of ZCCHC4 in disease progression remains largely unknown.
Through the quantitative analysis of m6A in different gastrointestinal tumor samples, we excluded HCC, which has been extensively studied in the past and focused on the role of m6A in CRC. A newly discovered methyltransferase was identified by continuing the analysis of the different components of m6A. In our study, we found that in CRC cell lines, high expression of ZCCHC4 promoted the proliferation, migration, and invasion of CRC cells, suggesting that ZCCHC4 is important for the survival and viability of CRC cells. Subcutaneous tumor-bearing experiments in mice also confirmed that the ZCCHC4 knockout inhibited CRC progression.
ZCCHC4, a member of the zinc finger protein family, is conserved in multicellular biological tissues but does not exist in yeast. ZCCHC4 contains a conserved catalytic domain, DPPF, which is a potential N6 methyladenosine methyltransferase domain [34,35,36]. In addition, ZCCHC4 contains CCHC Znf and Zf GRF domains [37]. Currently, there are very few reports on ZCCHC4. In January 2019, Chuan et al. found that ZCCHC4 is an m6A-modifying enzyme that mainly methylates 28 S rRNA and interacts with several mRNAs [38]. ZCCHC4 knockout reduced the m6A modification level of 28 S rRNA and the overall translation level and inhibited cell proliferation. Although many rRNA-related m6A modifying enzymes have been identified in bacteria, ZCCHC4 was the first rRNA-modified m6A methyltransferase identified in eukaryotes [39,40,41,42]. Currently, little is known about the function of ZCCHC4 in tumors. Our study first detected m6A levels in colorectal cancer tissues and further demonstrated that ZCCHC4 was a key factor in altering m6A levels, which made it possible to observe ZCCHC4, a newly identified methyltransferase from the perspective of m6A for the first time. Although ZCCHC4 can also participate in CRC development as an RBP, its specific mechanism of action has not been elucidated. Instead, m6A-seq continued to be used to search for downstream molecules of ZCCHC4 and finally confirmed that it could bind to lncGHRLOS to affect CRC progression [43,44,45,46].
However, we omitted delving into the bioinformatics aspects, primarily because of the contentious nature of ZCCHC4 in CRC. The novelty of ZCCHC4 as an m6A methylase raises concerns about the reliability of the results due to the limited available data. Our study conducted consistent clustering and principal component analysis (PCA) on the GSE 17,536 dataset, encompassing expression profiles and clinical data from 175 patients with CRC with a minimum overall survival of 30 days. The patients were categorized into two distinct groups: colon cancer 1 (N = 105) and colon cancer 2 (N = 70). Our findings indicated significant associations between RBM 15B, ZCCHC4, KIAA1429, IGF2BP2, and FTO and overall survival (OS) in patients with CRC. Using the median expression value of ZCCHC4 as a cut-off, patients were stratified into high- and low-expression groups. We established KM survival curves to evaluate the prognostic significance of these five genes in 175 patients with CRC from GSE 17,536. Surprisingly, ZCCHC4 did not demonstrate predictive value for OS in patients with CRC. Furthermore, analysis of the other four markers yielded negative results, highlighting the prognostic divergence between patients with negative and positive ZCCHC4 levels. Consequently, a comprehensive bioinformatics analysis is imperative for constructing a polygenic prognostic model.
ZCCHC4 was finally recognized because of its role as a member of the RBP family, which has been shown to play an important role in the development and treatment of cancer [47,48,49,50]. However, tumor-promoting RBP and their partners, which may serve as potential cancer therapeutic targets, need to be further identified. The results of this study suggest that ZCCHC4 is abnormally highly expressed in various digestive tract tumor tissues and may be a key regulator of CRC progression, and its knockdown can significantly inhibit tumor progression. Statistical analysis of the clinical features showed that low expression of ZCCHC4 was positively correlated with good prognosis and lymph node metastasis. We speculate that ZCCHC4 may have the opposite effect by regulating key target genes. In terms of mechanism, ZCCHC4 inhibits the apoptosis of CRC cells by interacting with a long-chain non-coding RNA (lncRNA) and promotes the proliferation, migration, and invasion of tumor cells. In conclusion, our study identified ZCCHC4 as a novel predictor of poor prognosis in CRC and a potential treatment target, providing mechanistic insights into its role in cancer progression.
Unlike epigenetic modifications such as DNA methylation and histone acetylation, which regulate gene expression at the transcriptional level, the anatomical mechanism of m6A-modified RNA has greatly improved our understanding of gene regulation at the post-transcriptional level in eukaryotic cells [51, 52]. The dynamic modification of m6A mRNA is necessary to maintain cell growth, metabolism, and differentiation. By regulating the abnormal modification of specific mRNA by m6A, tumor cells maintain self-repair and proliferation under different conditions [53]. In recent years, many studies have focused on the role of m6a methylated mRNA in cancer cells, and some lncRNAs have also been reported [54,55,56,57,58]. However, little is known about the regulatory role of m6A modifications of lncRNAs in CRC. Our study revealed that lncGHRLOS, which is not widely recognized, can combine with ZCCHC4 to control CRC by regulating KDM5D. Few studies have investigated its role in CRC, and no studies have explored its mechanism of action in CRC. Notably, we further confirmed that m6A methyltransferases can regulate the stability of lncRNAs in CRC cells, which expands our understanding of the involvement of lncRNAs in the pathogenesis of CRC.
Data availability
Data will be made available on request.
Change history
12 March 2024
A Correction to this paper has been published: https://doi.org/10.1186/s13046-024-03005-y
References
Deng X, Su R, Weng H, Huang H, Li Z, Chen J. RNA N (6)-methyladeno sine modification in cancers: current status and perspectives. Cell Res. 2018;28:507–17.
Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, et al. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methyla tion. Nat Chem Biol. 2014;10:93–5.
Ping XL, Sun BF, Wang L, Xiao W, Yang X, Wang WJ, et al. Mammalian WTAP is a regulatory subunit of the RNA N6-methyladenosine methyl transferase. Cell Res. 2014;24:177–89.
Li Z, Weng H, Su R, Weng X, Zuo Z, Li C, et al. FTO plays an oncogenic role in Acute myeloid leukemia as a N (6)-Methyladenosine RNA demethylase. Cancer Cell. 2017;31:127–41.
Zheng G, Dahl JA, Niu Y, Fedorcsak P, Huang CM, Li CJ, et al. ALKBH5 is a mammalian RNA demethylase that impacts RNA metabolism and mouse fertility. Mol Cell. 2013;49:18–29.
Wang X, Lu Z, Gomez A, Hon GC, Yue Y, Han D, et al. N6-methyladen osine-dependent regulation of messenger RNA stability. Nature. 2014;505:117–20.
Huang H, Weng H, Sun W, Qin X, Shi H, Wu H, et al. Recognition of RNA N (6)-methyladenosine by IGF2BP proteins enhances mRNA stability and translation. Nat Cell Biol. 2018;20:285–95.
Liu N, Dai Q, Zheng G, He C, Parisien M, Pan T. N (6)-methyladenosine dependent RNA structural switches regulate RNA-protein interactions. Nature. 2015;518:560–4.
Wang X, Zhao BS, Roundtree IA, Lu Z, Han D, Ma H, et al. N(6)-meth yladenosine modulates Messenger RNA translation E ciency. Cell. 2015;161:1388–99.
Geula S, Moshitch-Moshkovitz S, Dominissini D, Mansour AA, Kol N, Salmon Divon M, et al. Stem cells. m6A mRNA methylation facilitates resolution of naive pluripotency toward di erentiation. Science. 2015;347:1002–6.
Lin Z, Hsu PJ, Xing X, Fang J, Lu Z, Zou Q, et al. Mettl3– /Mettl14-mediated mRNA N(6)-methyladenosine modulates murine spermatogenesis. Cell Res. 2017;27:1216–30.
Li HB, Tong J, Zhu S, Batista PJ, Du y EE, Zhao J, et al. M (6) a mRNA methylation controls T cell homeostasis by targeting the IL-7/STAT5/SOCS pathways. Nature. 2017;548:338–42.
Zhang C, Chen Y, Sun B, Wang L, Yang Y, Ma D, et al. M(6) a modu lates haematopoietic stem and progenitor cell specification. Nature. 2017;549:273–6.
Li F, Yi Y, Miao Y, Long W, Long T, Chen S, et al. N (6)-Methyladenosine modulates nonsense-mediated mRNA decay in human glioblastoma. Cancer Res. 2019;79:5785–98.
Zhang C, Samanta D, Lu H, Bullen JW, Zhang H, Chen I, et al. Hypoxia induces the breast cancer stem cell phenotype by HIF-dependent and ALKBH5-mediated m6A-demethylation of NANOG mRNA. Proc Natl Acad Sci USA. 2016;113:E2047–56.
Yue B, Song C, Yang L, Cui R, Cheng X, Zhang Z, et al. METTL3-mediated N6-methyladenosine modification is critical for epithelial-mesenchymal transition and metastasis of gastric cancer. Mol Cancer. 2019;18:142.
Wu Y, Yang X, Chen Z, Tian L, Jiang G, Chen F, et al. M(6) A-induced lncRNA RP11 triggers the dissemination of colorectal cancer cells via upregulation of Zeb1. Mol Cancer. 2019;18:87.
Menon N, Patil V, Noronha V, Joshi A, Bhattacharjee A, Satam BJ, et al. Quality of life in patients with locally advanced head and neck cancer treated with concurrent chemoradiation with cisplatin and nimotu zumab versus cisplatin alone - additional data from a phase 3 trial. Oral Oncol. 2021;122:105517.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.
Borden EC. Interferons α and β in cancer: therapeutic opportunities from new insights. Nat Rev Drug Discov. 2019;18:219–34.
Ma H, Jin S, Yang W, Zhou G, Zhao M, Fang S, et al. Interferon-alpha enhances the antitumour activity of EGFR-targeted therapies by upregu lating RIG-I in head and neck squamous cell carcinoma. Br J Cancer. 2018;118:509–21.
Ma H, Chang H, Yang W, Lu Y, Hu J, Jin S. A novel IFNα-induced long noncoding RNA negatively regulates immunosuppression by interrupt ing H3K27 acetylation in head and neck squamous cell carcinoma. Mol Cancer. 2020;19:4.
Ma H, Yang W, Zhang L, Liu S, Zhao M, Zhou G, et al. Interferon-alpha promotes immunosuppression through IFNAR1/STAT1 signalling in head and neck squamous cell carcinoma. Br J Cancer. 2019;120:317–30.
Yang W, Jiang C, Xia W, Ju H, Jin S, Liu S, et al. Blocking autophagy flux promotes interferon-alpha-mediated apoptosis in head and neck squa mous cell carcinoma. Cancer Lett. 2019;451:34–47.
He L, Li H, Wu A, Peng Y, Shu G, Yin G. Functions of N6-methyladenosine and its role in cancer. Mol Cancer. 2019;18:176.
Chen Y, Zhao Y, Chen J, Peng C, Zhang Y, Tong R, et al. ALKBH5 suppresses malignancy of hepatocellular carcinoma via m(6) A-guided epigenetic inhibition of LYPD1. Mol Cancer. 2020;19:123.
Guo X, Li K, Jiang W, Hu Y, Xiao W, Huang Y, et al. RNA demethylase ALKBH5 prevents pancreatic cancer progression by posttranscriptional activation of PER1 in an m6A-YTHDF2-dependent manner. Mol Cancer. 2020;19:91.
Jin D, Guo J, Wu Y, Yang L, Wang X, Du J, et al. M (6) a demethylase ALKBH5 inhibits tumor growth and metastasis by reducing YTHDFs mediated YAP expression and inhibiting miR-107/LATS2-mediated YAP activity in NSCLC. Mol Cancer. 2020;19:40.
Tang B, Yang Y, Kang M, Wang Y, Wang Y, Bi Y, et al. M (6) a demethylase ALKBH5 inhibits pancreatic cancer tumorigenesis by decreasing WIF-1 RNA methylation and mediating wnt signaling. Mol Cancer. 2020;19:3.
Zhang S, Zhao BS, Zhou A, Lin K, Zheng S, Lu Z, et al. M (6) a demethylase ALKBH5 maintains tumorigenicity of Glioblastoma Stem-like cells by sustaining FOXM1 expression and cell proliferation program. Cancer Cell. 2017;31:591–606e596.
Wang J, Li Y, Wang P, Han G, Zhang T, Chang J, et al. Leukemogenic chromatin alterations promote AML Leukemia Stem cells via a KDM4C ALKBH5-AXL Signaling Axis. Cell Stem Cell. 2020;27:81–97. e88.
Zhao X, Cui L. Development and validation of a m(6) a RNA methylation regulators-based signature for predicting the prognosis of head and neck squamous cell carcinoma. Am J Cancer Res. 2019;9:2156–69.
Zhou X, Han J, Zhen X, Liu Y, Cui Z, Yue Z, et al. Analysis of genetic Altera tion signatures and prognostic values of m6A Regulatory genes in Head and Neck squamous cell carcinoma. Front Oncol. 2020;10:718.
Liu L, Wang J, Sun G, Wu Q, Ma J, Zhang X, et al. M (6) a mRNA methyla tion regulates CTNNB1 to promote the proliferation of hepatoblastoma. Mol Cancer. 2019;18:188.
Li T, Hu PS, Zuo Z, Lin JF, Li X, Wu QN, et al. METTL3 facilitates tumor progression via an m(6) A-IGF2BP2-dependent mechanism in colorectal carcinoma. Mol Cancer. 2019;18:112.
Gu C, Wang Z, Zhou N, Li G, Kou Y, Luo Y, et al. Mettl14 inhibits bladder TIC self-renewal and bladder tumorigenesis through N(6)-methyladenosine of Notch1. Mol Cancer. 2019;18:168.
Song H, Song J, Cheng M, Zheng M, Wang T, Tian S, et al. METTL3-medi ated m(6) a RNA methylation promotes the anti-tumour immunity of natural killer cells. Nat Commun. 2021;12:5522.
Lu M, Zhang Z, Xue M, Zhao BS, Harder O, Li A, et al. N (6)-methyladeno sine modification enables viral RNA to escape recognition by RNA sensor RIG-I. Nat Microbiol. 2020;5:584–98.
Lan Q, Liu PY, Haase J, Bell JL, Huttelmaier S, Liu T. The critical role of RNA m (6)a methylation in Cancer. Cancer Res. 2019;79(7):1285–92.
Meyer KD. Ja Rey SR. Rethinking m (6)a readers, writers, and Erasers. Annu Rev Cell Dev Biol. 2017;33:319–42.
Yang C, Hu Y, Zhou B, Bao Y, Li Z, Gong C, et al. The role of m (6)a modifica tion in physiology and disease. Cell Death Dis. 2020;11(11):960.
Wang J, Wang J, Gu Q, Ma Y, Yang Y, Zhu J, et al. The biological function of m6A demethylase ALKBH5 and its role in human disease. Cancer Cell Int. 2020;20:347.
Huang X, Lv D, Yang X, Li M, Zhang H. m6A RNA methylation regulators could contribute to the occurrence of chronic obstructive pulmonary disease. J Cell Mol Med. 2020;24(21):12706–15.
Chen M, Wong CM. The emerging roles of N6-methyladenosine (m6A) deregulation in liver carcinogenesis. Mol Cancer. 2020;19(1):44.
Shen C, Xuan B, Yan T, Ma Y, Xu P, Tian X, et al. M (6)A-dependent glycolysis enhances colorectal cancer progression. Mol Cancer. 2020;19(1):72.
Shen C, Sheng Y, Zhu AC, Robinson S, Jiang X, Dong L, et al. RNA demethylase ALKBH5 selectively promotes Tumorigenesis and Cancer Stem Cell Self Renewal in Acute myeloid leukemia. Cell Stem Cell. 2020;27(1):64–80e9.
Malacrida A, Rivara M, Di Domizio A, Cislaghi G, Miloso M, Zuliani V, et al. 3D proteome-wide scale screening and activity evaluation of a new ALKBH5 inhibi tor in U87 glioblastoma cell line. Bioorg Med Chem. 2020;28(4):115300.
Tang B, Yang Y, Kang M, Wang Y, Wang Y, Bi Y, et al. M (6)a demethylase ALKBH5 inhibits pancreatic cancer tumorigenesis by decreasing WIF-1 RNA methylation and mediating wnt signaling. Mol Cancer. 2020;19(1):3.
Zhu Z, Qian Q, Zhao X, Ma L, Chen P. N (6)-methyladenosine ALKBH5 pro motes non-small cell lung cancer progress by regulating TIMP3 stability. Gene. 2020;731:144348.
Zhang C, Zhi WI, Lu H, Samanta D, Chen I, Gabrielson E, et al. Hypoxia inducible factors regulate pluripotency factor expression by ZNF217- and ALKBH5-mediated modulation of RNA methylation in breast cancer cells. Oncotarget. 2016;7(40):64527–42.
Wang J, Li Y, Wang P, Han G, Zhang T, Chang J, et al. Leukemogenic chromatin alterations promote AML Leukemia Stem cells via a KDM4C ALKBH5-AXL Signaling Axis. Cell Stem Cell. 2020;27(1):81–97e8.
He Y, Hu H, Wang Y, Yuan H, Lu Z, Wu P, et al. ALKBH5 inhibits pancreatic Cancer motility by decreasing long non-coding RNA KCNK15-AS1 methylation. Cell Physiol Biochem. 2018;48(2):838–46.
Toledo CM, Ding Y, Hoellerbauer P, Davis RJ, Basom R, Girard EJ, et al. Genome-wide CRISPR-Cas9 screens reveal loss of redundancy between PKMYT1 and WEE1 in Glioblastoma Stem-like cells. Cell Rep. 2015;13(11):2425–39.
Xuan ZH, Wang HP, Zhang XN, Chen ZX, Zhang HY, Gu MM. PKMYT1 aggravates the progression of ovarian cancer by targeting SIRT3. Eur Rev Med Pharmacol Sci. 2020;24(10):5259–66.
Wang J, Wang L, Chen S, Peng H, Xiao L, Du E, et al. PKMYT1 is associated with prostate cancer malignancy and may serve as a therapeutic target. Gene. 2020;744:144608.
Zhang Q, Zhao X, Zhang C, Wang W, Li F, Liu D, et al. Overexpressed PKMYT1 promotes tumor progression and associates with poor survival in esophageal squamous cell carcinoma. Cancer Manag Res. 2019;11:7813–24.
Long HP, Liu JQ, Yu YY, Qiao Q, Li G. PKMYT1 as a potential target to improve the Radiosensitivity of Lung Adenocarcinoma. Front Genet. 2020;11:376.
Sun QS, Luo M, Zhao HM, Sun H. Overexpression of PKMYT1 indicates the poor prognosis and enhances proliferation and tumorigenesis in non-small cell lung cancer via activation of Notch signal pathway. Eur Rev Med Pharmacol Sci. 2019;23(10):4210–9.
Acknowledgements
We would like to thank Professor Wang Fengsong (School of Life Sciences, Anhui Medical University) for his comments on the revision.
Funding
Natural Science Foundation in Anhui Province of China (grant number: 2008085MH279);Key project of Anhui Translational Medicine Research Institute (grant number: 2022zhyx-B08); Clinical Medicine Discipline construction project of Anhui Medical University (grant number: 9301001807).
Author information
Authors and Affiliations
Contributions
Wenbin Wang and Zhengguang Wang provided grant support, conceived the project and supervised the work. Ke Chen, Jingcheng Zhang, Lei Meng and Lingshang Kong performed the experiments. Ke Chen performed all analysis. Ming Lu and Lingshang Kong both contributed to the manuscript and figures. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval
The research protocol was approved by the clinical research ethics committee of Anhui Medical University.
Animal studies
The animal studies protocol was approved by the animal research ethics committee of Anhui Medical University.
Conflict of interest
The funders had no role in study design, data collection and analysis, manuscript preparation, or submission for publication. All those named as authors have made a sufficient contribution to the final manuscript for submission. All authors declare no conflict of interest, and agree on consent to be listed as coauthors.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: errors were spotted in some of the entries in Table 1, Table2 and Figure 6
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Chen, K., Zhang, J., Meng, L. et al. The epigenetic downregulation of LncGHRLOS mediated by RNA m6A methylase ZCCHC4 promotes colorectal cancer tumorigenesis. J Exp Clin Cancer Res 43, 44 (2024). https://doi.org/10.1186/s13046-024-02965-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13046-024-02965-5