Objective Dysregulated cellular metabolism is a distinct hallmark of human colorectal cancer (CRC). However, metabolic programme rewiring during tumour progression has yet to be fully understood.
Design We analysed altered gene signatures during colorectal tumour progression, and used a complex of molecular and metabolic assays to study the regulation of metabolism in CRC cell lines, human patient-derived xenograft mouse models and tumour organoid models.
Results We identified a novel RNA-binding protein, RALY (also known as hnRNPCL2), that is highly associated with colorectal tumour aggressiveness. RALY acts as a key regulatory component in the Drosha complex, and promotes the post-transcriptional processing of a specific subset of miRNAs (miR-483, miR-676 and miR-877). These miRNAs systematically downregulate the expression of the metabolism-associated genes (ATP5I, ATP5G1, ATP5G3 and CYC1) and thereby reprogramme mitochondrial metabolism in the cancer cell. Analysis of The Cancer Genome Atlas (TCGA) reveals that increased levels of RALY are associated with poor prognosis in the patients with CRC expressing low levels of mitochondrion-associated genes. Mechanistically, induced processing of these miRNAs is facilitated by their N6-methyladenosine switch under reactive oxygen species (ROS) stress. Inhibition of the m6A methylation abolishes the RALY recognition of the terminal loop of the pri-miRNAs. Knockdown of RALY inhibits colorectal tumour growth and progression in vivo and in organoid models.
Conclusions Collectively, our results reveal a critical metabolism-centric role of RALY in tumour progression, which may lead to cancer therapeutics targeting RALY for treating CRC.
- colorectal cancer
- molecular mechanisms
- gene expression
- molecular biology
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information. The authors are willing to share the data related to any of the work reported.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
Significance of this study
What is already known on this subject?
Dysregulation of metabolism in cancer cells is dynamically regulated at many cancer hallmarks.
Metabolic reprogramming in cancer cells could be impacted by post-translational events and epigenetic regulation of genes.
RALY was identified in the spliceosome complex and in the Drosha microprocessor complex.
The m6A modification is found to alter RNA structure to enhance HNRNPC binding, and HNRNPA2B1 could bind to the m6A marks in some pri-miRNA transcripts to promote their processing.
What are the new findings?
RALY is highly associated with metastasis of malignant colorectal cancer (CRC).
RALY acts as a key regulatory component in the Drosha complex, and promotes the post-transcriptional processing of a specific subset of miRNAs.
RALY-mediated miRNAs processing is modulated by N6-methyladenosine switch, and reprogrammes mitochondrial metabolism by downregulating the electron transport chain genes.
Inhibition of RALY suppresses colorectal tumour growth and progression in vivo and in the organoid models.
How might it impact on clinical practice in the foreseeable future?
Our report first reveals a critical metabolism-centric role of RALY in tumour progression. It provides a novel therapeutic strategy to target RALY for CRC treatment.
The acquired capability of unrestrained replicative potential has been a representative hallmark for cancer cells.1 To meet the requirement for hyperproliferation, cancer cells shift intracellular energy metabolism and mainly rely on glycolysis pathways instead of mitochondrial oxidative phosphorylation (OXPHOS) to support their energy demand and biomass synthesis, even when the oxygen level is sufficient.1 2 Dysregulation of metabolism in cancer cells is actively regulated at many cancer hallmarks, including proliferation, angiogenesis, invasion and metastasis.3 4 In a recent study, we reported that USP13 as a key regulator rewires energy metabolism through deubiquitinating ATP citrate synthase (ACLY) and 2-oxoglutarate dehydrogenase (OGDH) and promotes ovarian tumour progression.5 Emerging evidence has shown that metabolic reprogramming in cancer cells could be impacted by post-translational events and epigenetic modulation of gene expression.6–8 Identification of key regulators in mediating cancer cell metabolism could lead to a new therapeutic approach for cancer treatment.
As ‘energy factories’, mitochondria control the bioenergetics, biosynthetic and signalling organelles of cancer cells and serve as an attractive target for cancer therapy.9–11 Previous reports showed that inhibition of electron transport chain (ETC) genes could sensitise cancer cells to glucose depletion and reduce tumour progression.12 13 However, regulation of ETC gene expression in tumour progression remains largely unknown. Recent studies indicate that epigenetic regulation has an essential role in controlling cellular signalling pathways during tumourigenesis.14 15 Delineation of multifaceted epigenetic changes in mitochondrial metabolism during tumour progression would provide insights into the crucial mechanism by which cancer cells adapt to a changing microenvironment.
Recently, miRNAs have been found to act as master gene regulators, similar to transcription factors, to mediate gene expression patterns in various cellular activities, including cell metabolism.16 17 C-Myc suppression of miR-23a and miR-23b enhances mitochondrial glutaminase expression, leading to elevated glutamine catabolism in prostate cancer and B lymphoma.18 The RNA binding protein Lin28 inhibits let-7 biogenesis to promote the insulin-PI3K-mTOR pathway and increase glucose uptake.19 MiR-103, miR-107 and miR-143 are involved in regulating insulin-mediated glucose metabolism,20 21 while miR-33 cooperates with SREBPs to target ABCA1 and coordinate cholesterol biosynthesis and trafficking.22 MiRNAs have been known to provide a functional link between metabolic reprogramming and tumourigenesis. In this study, we identify a novel RNA-binding protein, RALY (HnRNP-Associated with Lethal Yellow, also called hnRNPCL2), as a crucial indicator for tumour progression in colorectal cancer (CRC). We present here a role of RALY in reprogramming mitochondrial metabolism via mediating miRNA processing in CRC, and highlight the potential therapeutic strategy by targeting RALY in CRC therapy.
Materials and methods
As previously described,23 RNA immunoprecipitation (RIP) assay was performed using the Magna RIP RNA-Binding Protein Immunoprecipitation Kit (17–700, Merck Millipore, USA) in accordance with manufacturer’s protocol. Stable RALY-overexpressed HEK293T cells were collected and total RNA was obtained by re-suspending and dispersing the cell pellet in complete RIP Lysis Buffer. Antibodies of RALY (ab170105, Abcam) or negative control normal rabbit IgG (PP64B) in magnetic beads were added to pull down specific RNAs. Enriched microRNAs were subject to high-throughput NGS sequencing and validated by quantitative RT-PCR.
ECAR and oxygen consumption rate (OCR) were monitored using a Seahorse Bioscience Analyzer (XF96) according to the manufacturer’s instructions. Cells were seeded in 96-well plates at a density of 5–8×103 cells per well with growth medium for at least 18 hours. Next, the medium was changed to base medium (DMEM, 143 mM NaCl, phenol red, pH 7.35). For ECAR analysis, cells were added with the medium containing 2 mM glutamine (pH 7.35) and monitored every 3 min following successive administration of 10 mM of glucose, and inhibitors (1 µM oligomycin and 50 mM 2-deoxyglucose). For OCR analysis, cells were added with Mito stress-test base medium (10 mM glucose, 2 mM glutamine, 1 mM pyruvate, pH 7.4) and monitored every 3 min following successive administration of inhibitors (2 µM oligomycin, 2 µM FCCP, or 0.5 µM rotenone/antimycin A).11
A detailed description of the materials and methods can be found in online supplemental materials and methods.
Upregulation of RALY is associated with aggressiveness in CRC
Specific focal gene amplifications, gene losses and altered gene-expression levels have been identified as signature events for CRC progression.24 We conducted intensive analysis of gene alterations in CRC from the TCGA database by using genomic set enrichment analysis (GSEA), and identified a specific group of genes related to RNA alternative splicing that is highly associated with tumour metastasis (figure 1A). Among those genes, a novel RNA binding protein, RALY was identified as a molecular signature associated with tumour aggressiveness from the integrative analysis of tumour stages, lymph node invasion, vascular invasion and distant metastasis status at the time of surgery (figure 1B). Amplification of the RALY gene occurs frequently in human cancers and locates in the proximity of E2F1 (figure 1C and D). Upregulation of RALY is found in 56% of CRC cases (figure 1E). We demonstrate that RALY is overexpressed in nine human CRC cell lines as compared with the normal human colon epithelial cell line (figure 1F–H). By using CRC patient tissue samples locally procured, we further validated the upregulation of RALY in the CRC tumours in comparison with their adjacent normal tissues (figure 1I and J and online supplemental figure 1) and (online supplemental table 1). Analysis of TCGA database and the local CRC database revealed that expression of RALY is tightly associated with its gene copy number (figure 1K) and high expression of RALY correlates with poor clinical outcome in CRC patients (figure 1L and online supplemental figure 2A). By correlation analysis with clinical data from TCGA (colon adenocarcinoma (COAD)+rectum adenocarcinoma (READ)), the expression of RALY is significantly correlated with several univariate clinical factors including tumour stage and distant metastasis, but not with the status of KRAS, BRAF and MSI, indicating that RALY may be used as an independent prognostic biomarker for CRC metastasis (online supplemental figure 2B) and (online supplemental table 2).
RALY interacts with the Drosha complex and regulates a subset of miRNA expression
Previously microRNA biogenesis has been shown to coordinate with alternative splicing under stress in the supraspliceosome,25 26 and RALY was identified in the spliceosome complex27 and in the Drosha microprocessor complex.28 However, the functional role of RALY in microRNA biogenesis remains unknown. We examined the physical interaction of RALY with the Drosha microprocessor (figure 2A). Endogenous RALY was validated in the Drosha immunoprecipitate, and conversely, Drosha and DGCR8 were detected in the RALY immunoprecipitate. We also verified the colocalisation of RALY and Drosha in the nucleus (online supplemental figure 2C). To examine the effect of RALY on miRNA biogenesis, we performed the analysis of global mRNA and miRNA expression by next-generation sequencing in the RALY-knockdown HCT116 cells, and by RALY RIP sequencing in the RALY-overexpressed 293 T cells (figure 2B and online supplemental table 3). Depletion of RALY significantly decreased the expression levels of a subset of 17 miRNAs (cutoff >2 fold), while overexpression of RALY significantly increased their level enrichment in the RALY immunoprecipitate complex (figure 2C and D). To determine whether RALY regulates the miRNA expression at transcriptional or post-transcriptional levels, we detected the levels of both pri-miRNAs and mature miRNAs in the RALY-knockdown or RALY-overexpression HCT116 cells. As a control, miR-21 was used in the assay as it does not depend on the RALY expression level. No notable differences were detected on the transcription of pri-miRNAs including pri-miR-96, pri-miR-483, pri-miR-618, pri-miR-676 and pri-miR-877 in the control cells and RALY-knockdown or RALY-overexpression cells. However, levels of mature forms of the RALY-dependent miRNA, but not control miR-21, were dramatically decreased when RALY was knocked down in HCT116 cells. Conversely, these RALY-dependent miRNAs were elevated in the RALY-overexpression HCT116 cells (figure 2E,F). To validate the universal RALY-dependency of these miRNAs in CRC, we further analysed the mature forms of these miRNAs in eight more CRC cell lines, including HCT8, HCT15, LOVO, RKO, HT29, SW620, SW480 and HT55 cell lines. Consistently, knockdown of RALY significantly suppressed the levels of RALY-dependent miRNAs, but not the level of the control miR-21, in these cell lines (figure 2G). These results indicate that RALY promotes the expression of specific miRNAs at the post-transcriptional level.
RALY reprogrammes mitochondrial metabolism in CRC through mediating miRNA expression
To investigate the functional role of RALY in CRC tumour progression, we establish stable RALY-knockdown models in HCT116, HCT15 and LOVO cells with two independent shRNAs (shRALY#1 and #2) (online supplemental figure 3A). To determine the biological function of RALY, we performed GSEA in the RALY-knockdown HCT116 cells, and found that the OXPHOS pathways including energy metabolism and mitochondrial respiratory ETC are significantly enriched in the RALY-depleted cells (figures 2B and 3A). Knockdown of RALY induced the expression of some mitochondrion-related ETC genes, including ATP5I, ATP5G1, ATP5G3, ATP5J2 and CYC1 at both RNA and protein levels (figure 3B and C). We reasoned that knockdown of RALY reprograms mitochondrial metabolism in CRC, and thereby limits the synthesis of intermediates in the TCA cycle from the OXPHOS pathways. Marked reduction of extracellular acidification rate (ECAR) levels was observed in the RALY-knockdown cells, indicating the influence of RALY on aerobic glycolysis, while increased basal OCR levels in the RALY-knockdown cells suggested a role of RALY in aerobic respiration. Conversely, overexpression of RALY elevated the ECAR levels and decreased the OCR levels (figure 3D and E and online supplemental figure 3B-F). Significant reduction of ATP levels and induction of ROS levels were also observed in the RALY-knockdown cells and the opposites were seen in the RALY-overexpressed cells (figure 3F and G and online supplemental figure 3GH). As a rapid turnover of NAD+ is required in cancer cells due to their elevated metabolic needs,29 depletion of RALY significant increased the NAD+levels and the NAD+/NADH ratios, indicating reduced reductive flux (figure 3H and online supplemental figure 3I,J). In line with our hypothesis, we find that knockdown of RALY significantly increased the expression of ETC genes, leading to reduced ATP production, induced ROS level and the NAD+/NADH ratio in CRC cells.
To validate the effects of RALY on mitochondrial metabolism are mainly mediated by mitochondrion-related ETC genes, we have constructed stable knockdown CRC cells (ATP5I-KD, ATP5G1-KD, ATP5G3-KD and CYC1-KD) by lentiviral shRNAs (figure 3I and online supplemental figure 4A,B). Knockdown of ATP5I, ATP5G1, ATP5G3 and CYC1 significantly reduced ROS levels in the RALY-depleted CRC cells (figure 3J and supplemental figure 4C). To determine how RALY regulates those mitochondrion-related ETC genes, we conducted miRNA-target prediction by using comprehensive tools including TargetScan, MiRanda and miRBase, and identified six specific binding sites at the 3' untranslated region (3’UTR) of those ETC genes for the RALY-dependent miRNAs including miR-483, miR-676 and miR-877 (figure 2B and supplemental figure 4D). To prove the direct effects of these miRNAs on ETC gene expression, we cloned the ETC genes’ 3’UTR sequences into a dual luciferase reporter construct and generated deletion mutant vectors (online supplemental figure 4DE). Overexpression of miR-483, miR-676 and miR-877 reduced the relative luciferase activity of the targeted ETC 3’UTRs while their specific inhibitors (antagomiRs) increased the luciferase activity signals. Deletion mutants on the miRNA-binding sites blocked the regulatory effects mediated by these miRNAs (online supplemental figure 4F). Consistently, overexpression of miR-483, miR-676 and miR-877 downregulated the expression of the target ETC genes at both RNA and protein levels (figure 3K and L and online supplemental figure 5A-D). Conversely, overexpression of miR-483, miR-676 and miR-877 rescued the effects on mitochondrial metabolism caused by RALY depletion in CRC cells (figure 3M and supplemental figure 5E-L). These results illustrate that RALY modulates mitochondrial metabolism through mediating miRNA expression.
RALY regulates CRC cell proliferation and apoptosis under ROS stress
Mitochondrial metabolism is a key role in cellular energy, biosynthesis and cell growth, and ROS is controlled by mitochondria. It has been shown to affect the tumourigenic and metastatic phenotypes.30 Since our results indicate that RALY regulates ROS production by mediating miRNA processing, we postulated that RALY may participate in mitochondrial ROS-involved stress signalling and thus affect tumour growth. Knockdown of RALY remarkably suppressed CRC cell proliferation as indicated by decreasing cell numbers (figure 4A and B and online supplemental figure 6AB), reducing EDU staining (figure 4E and F) and fewer cellular survival dots (figure 4G). To determine whether the mitochondrion-mediated ROS production by RALY is crucial for cell growth in CRC, we use a ROS scavenger, N-acetyl L-cysteine (NAC), to deplete the ROS level in cancer cells. Treatments of NAC partially rescued cancer cell growth in the RALY-knockdown cells, but the rescue effect was not dependent on NAC concentrations (figure 4C, E and G). Similarly, overexpression of ETC genes-related miRNAs (miR-483, miR-676 and miR-877) restored cell growth and proliferation in the RALY-knockdown cells (figure 4D, F and G). We next examined the function of RALY in cell survival by measuring the levels of Ki67 and caspase-3. RALY depletion induced significant apoptosis (indicated by cleaved caspase-3) by flow cytometry. Treatment of NAC or overexpression of the miRNAs related to ETC genes restored cell survival in the RALY-knockdown cells (figure 4H1), as validated by reduction of the ROS level in these rescue of cells compared with the RALY-knockdown cells (figure 4J). These results suggest that RALY plays a crucial role in controlling CRC cell growth and cell death by mediating mitochondrial metabolism.
RALY mediates miRNA processing through N6-methyladenosine switch
We next looked into the detailed molecular mechanism by which RALY efficiently recruits pri-miRNAs to the Drosha complex and promotes their processing. It has been previously known that the N6-methyladenosine (m6A) mark is associated with the processing of nuclear transcripts,31 and alternative splicing.32 The m6A modification is able to alter RNA structure to enhance HNRNPC binding,33 and HNRNPA2B1 could bind to m6A marks in some pri-miRNA transcripts to promote their processing.34 As a member of the HNRNPC subfamily, RALY shares an RRM motif and an acidic-rich auxiliary domain for formation of a strong and specific RNA interaction.35 We postulated that RALY may promote pri-miRNA processing via the m6A modification of the pri-miRNA. We identified the binding motif of RALY (HHCTG) in the pri-miRNAs from the RALY-RIP sequencing (figure 5A). By comparing the miRNA profiling in the METTL3-depleted cells (GSE70061) and RALY-depleted cells, 5 miRNAs were found in the overlapped regions, including miR-877 (figure 5B and online supplemental figure 7A). Knockdown of METTL3 downregulated the expression levels of miR-483, miR-676 and miR-877 in CRC cells (figure 5C and online supplemental figure 7B). The m6A signal is enriched in pri-miR-877, supporting the potential role of m6A modification in the pri-miRNA processing (figure 5D). As the RGAC motif is the recognition motif for METTL3-mediated m6A methylation,34 three potential m6A sites on miR-877 were identified and indicated as A1, A2 and A3 (figure 5E). We performed RNA pull-down assay with methylated oligonucleotides to validate the authentic m6A sites. The results showed that only m6A-methylation on A3 could significantly increase the pull-down amount of RALY (figure 5F and online supplemental figure 7C), and A3 locates at the terminal loop region of miR-877. Collectively, these results validate that m6A modification could alter its local RNA structure at the terminal loop region of pri-miRNA and enhance the accessibility of its nearby region to modulate RALY binding, thus promoting processing of miRNAs (figure 5G).
METTL3 regulates RALY-mediated mitochondrial metabolism in CRC
Since METTL3-mediated m6A modification is involved in the RALY-dependent miRNA processing, we next asked whether METTL3 plays a role in the RALY-mediated mitochondrial metabolism. By comparing METTL3-dependent m6A gene peaks (GSE94613) and RALY-dependent gene peaks, 454 altered genes were identified in the overlapped region. The GSEA analysis highlighted the overlapped genes enriched in metabolism-related pathways, including OXPHOS (figure 6A). Depletion of METTL3 by siRNAs remarkably reduced the expression of miR-483, miR-676 and miR-877 in the RALY-overexpressed HCT116 cells (figure 6B). We next performed RIP assays to detect whether the pri-miRNA-binding activity of RALY is affected by METTL3. Significantly decreased levels of pri-miR-483, pri-miR-676 and pri-miR-877 were bound by RALY after METTL3 depletion in HCT116 cells (figure 6C). We further generated stable METTL3-knockdown models in HCT116, HCT15 and LOVO cells with two independent shRNA sequences (shMETTL3#1 and #2). Knockdown of METTL3 significantly suppressed cell growth in CRC cells (online supplemental figure 8B). Accordingly, depletion of METTL3 induced the expression of ATP5I, ATP5G1, ATP5G3 and CYC1 at both protein and RNA levels (figure 6D and online supplemental figure 8A). Loss of METTL3 promoted mitochondrial respiration as indicated by increased OCR levels, and suppressed aerobic glycolysis as indicated by decreased ECAR levels in CRC cells. Overexpression of miR-877 rescued the effects on mitochondrial respiration and aerobic glycolysis caused by METTL3 depletion (figure 6E,F and online supplemental figure 8C,D). In line with it, knockdown of METTL3 significantly reduced ATP production and enhanced ROS levels, while overexpression of miR-877 diminished the effects caused by the METTL3 knockdown (figure 6G,H and online supplemental figure 8EF). A notable increase in the NAD+/NADH ratio was observed in the METTL3-depleted CRC cells, but overexpression of miR-877 abolished the induction (figure 6I and online supplemental figure 8G). These results suggest that METTL3-dependent m6A modification in pri-miRNAs is essential for the RALY-mediated mitochondrial metabolism in CRC cells.
RALY serves as a potential therapeutic target in CRC
We next examined the functional role of RALY in tumourigenesis in vivo. We generated inducible RALY-knockdown syngeneic mouse tumour models by subcutaneously injection of mouse MC38 cells into C57BL/6 mice (figure 7A). After initial tumour establishment, administration of doxycycline (1.0 mg/mL) in drinking water suppressed RALY expression and consequently inhibited tumour growth (figure 7B–E). Knockdown of RALY in the tumours induced the expression of ATP5I and decreased miR-676 expression, as indicated by Western blotting assay and RT-PCR, respectively (figure 7B,C). We then established an orthotopic mouse xenograft CRC model by injection of inducible RALY-knockdown LOVO cells into the cecal wall of NOD-SCID mice. The group treated with doxycycline showed fewer and smaller tumour nodules in the colon compared with the group treated with water (figure 7F and online supplemental figure 9A), as indicated by the H&E staining of the colon and tumour weight (figure 7G and H and online supplemental figure 9B). Using the organoid models generated from the MC38-derived tumours, we find that knockdown of RALY by treatment with doxycycline could significantly inhibit the growth of the organoids (figure 7I,J), and induced apoptosis in the organoid (figure 7K). To verify whether the inhibitory effect by RALY in vivo is mediated by RALY-dependent miRNAs, we overexpressed miR-676 in the RALY-knockdown HCT116 cells. Notably, overexpression of miR-676 in the RALY-knockdown HCT116 cells partly rescued xenograft tumour growth (figure 7L–N), and suppressed the expression of ATP5I in comparison with the RALY-knockdown xenograft tumours, as indicated by immunohistochemistry (IHC) and Western blotting assays (figure 7O–Q and online supplemental figure 9C).
Further, we generated CRC patient-derived xenograft (PDX) models and depleted RALY expression by using RNAi and nano-gold particle delivery tools (figure 8A). We employed two treatment strategies for the PDX models. The first was to administrate with the RNAi regimen on the second day after the tumours were embedded, while the second was to administrate with the RNAi regimen on the 20th day when the tumours regained growth to the size of 100 mm3. Our results show that both treatment strategies led to significantly reduced tumour growth (figure 8B–D) and induced the expression of ATP5I, as indicated by the H&E and IHC staining assays (figure 8E,F). We also established the patient-derived organoid (PDO) models from fresh CRC patient tumours. Treatment with siRALY dramatically reduced the organoid viability and made it vulnerable in shape and membrane permeability (figure 8G,H), and induced apoptosis phenotype in the model (figure 8I). To determine whether various levels of RALY would affect the therapeutic effects of siRNA treatments, we select 3 CRC patient tissues with various levels of RALY (online supplemental figure 10A). Targeting RALY by siRNAs in CRC tissue with high level of RALY achieved the maximal inhibitory effects on tumour growth in both PDO model and PDX model, while it also showed effective inhibition of tumour growth in CRC tissues with normal or low level of RALY (online supplemental figure 10B-E).
In analysis of clinical database (TCGA-COAD), the expression levels of RALY and METTL3 are significantly higher in tumours compared with normal tissues, while the levels of ATP5I, ATP5G1 and ATP5G3 are significantly lower in tumours (online supplemental figure 11A). Lower expression levels of these ETC genes indicate poorer clinical outcome for CRC patients (figure 8J). Further, we showed the positive correlation between RALY and its dependent miRNAs in CRC clinical samples from various CRC datasets (online supplemental figure 11BC) and (online supplemental table 4). Higher expression level of miR-877 and miR-483 shows significant poorer clinical outcome for CRC patients, which is consistent with higher expression level of RALY (online supplemental figure 11DE). Collectively, these results suggest that RALY promotes tumour growth in vivo and in the organoid models, highlighting its role as a potential therapeutic target in CRC.
Systematic analysis of genome-wide miRNA expression profiles indicates that aberrant miRNA expression contributes to cellular transformation and tumour progression.36 37 The stem-loop structure of miRNA was found to be critical for the Drosha-mediated cleavage.38 39 Current studies, including ours, have shown that RNA-binding proteins confer the specificity for recruiting and processing individual miRNAs.34 40–43 A growing list of RNA-binding proteins identified from the Drosha microprocessor complex is the members of the heterogeneous nuclear ribonucleoprotein family, including hnRNPA1, hnRNPA2B1, hnRNPK and hnRNPL.34 43–45 In this study, we identify a novel member of hnRNP, RALY, which is highly associated with CRC stages and aggressiveness, interacts with the Drosha complex and promotes the expression of a subset of miRNAs including miR-483, miR-676 and miR-877. Previously miR-483 was reported to be involved in initiation and invasion of several cancers.46–48 We show that these RALY-dependent miRNAs downregulate the ETC genes, ATP5I, ATP5G1, ATP5G3 and CYC1, and mediate reprogramming of mitochondrial metabolism in CRC cells, which provides the functional link between epigenetic regulation and mitochondrial metabolism. Our results demonstrate that inhibition of RALY abolishes ATP production, induces ROS level and increases the NAD+/NADH ratio in CRC cells, while overexpression of miR-483, miR-676 and miR-877 in the RALY-depleted cells could rescue the metabolic phenotype. We further knocked down the related ETC genes in RALY-depleted cells and validated that the metabolism-centric effects of RALY in OXPHOS are mainly mediated by these ETC genes (figure 9).
Recently, N6-methyladenosine modification has been shown as critical post-transcriptional regulation for gene expression in tumour development, and plays important roles in RNA splicing and degradation.49–52 Previous studies reported that hnRNP proteins promoted some miRNAs processing by recognising the m6A-modified pri-miRNAs transcripts.32–34 As a member of the hnRNPC subfamily, RALY is now found to recognise the terminal loop region of pri-miR-877 with m6A modification, and facilitate the binding of pri-miR-877 to the microprocessor complex, thus promoting its processing. Knockdown of METTL3, a primary m6A methylase, abolishes the upregulation of RALY-dependent miRNAs, and in part affects the RALY-mediated mitochondrial metabolism in CRC cells. Whether RALY acts as a direct m6A reader or modulates the miRNA processing through the m6A switch raised an intriguing question. Previously, hnRNPA2B1 was reported as a direct m6A reader to promote some miRNAs processing.34 However, subsequent structure analysis suggested that hnRNPA2B1 may promote pri-miRNAs processing through the m6A switch instead of acting as a direct m6A reader.53 Of note, hnRNPC was reported to modulate alternative splicing and binding activities of target RNAs through the m6A-switch-regulated mode to gain the abundance of gene expression.33 For RALY, we performed RNA pull-down assay with methylated oligonucleotides and found that the functional m6A site locates at the terminal loop region of pri-miR-877, adjacent to the RALY-binding motif, indicating that RALY may mediate pri-miRNA processing through the m6A switch mode.
Combinational therapy with BACH1 and ETC inhibitors showed favourable therapeutic effects in triple-negative breast cancer.11 Similarly, we demonstrate that RALY is required for the metastasis of aggressive CRC, and its gene expression is highly associated with poor clinical outcomes. Similar to BACH1, RALY functions as a negative regulator to regulate expression of the ETC genes and reprogramme mitochondrial metabolism, suggesting that it may be a promising target for cancer therapy. Inhibition of RALY sensitised resistant CRC cells to oxaliplatin.54 In terms of translational development, we developed a core-shell nano-gold platform for delivering in vivo optimised RALY siRNA to target RALY as our previous ‘nano-bomb’ platform.55 By using human PDX model and PDO model from clinical patient samples, we validated the significance of targeting the RALY-miRNAs-mediated mitochondrial metabolism axis with our nano-gold technology in CRC therapy. Taken together, our study not only reveals a crucial metabolism-centric role of RALY in reprogramming OXPHOS mediated by the m6A switch, but also highlights a potential therapeutic target in CRC.
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information. The authors are willing to share the data related to any of the work reported.
The authors thank Dr Bo Wei for help with CRC tissue isolation. They also acknowledge the support of the SYSU Animal Centre.
LS and AW are joint first authors.
LS and AW contributed equally.
Contributors Study concept and design: GW, LS and AW; Acquisition of data: LS, AW, ZZ, DC, HL, CL, SY, YN, ZL and SZ; Analysis and interpretation of data: GW, AX, XL, LS and AW; Statistical analysis: LS, AW, ZZ and DC; Drafting of the manuscript: GW and LS; Critical revision and final approval of the manuscript: GW, XL and AX; Administrative, technical support: SW, XB and WH; Study supervision: GW, XL and AX.
Funding This work was supported in part by grants from the National Natural Science Foundation of China (82073869, 31701114, 21672266, 81871994, 81701834, 91942301); National Key Research and Development Project (2019YFC1710104) from the Ministry of Science and Technology of China; Guangdong Basic and Applied Basic Research Foundation (2019A050510019, 2019B151502063); Guangdong Provincial Key Laboratory of Construction Foundation (2017B030314030, 2020B1212060034); Guangzhou Science and Technology Planning Programme (202002020051, 201902020018); National Engineering Research Centre for New Drug and Druggability Evaluation, Seed Program of Guangdong Province (2017B090903004).
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.