HNRNPA2B1-mediated m6A modification of TLR4 mRNA promotes progression of multiple myeloma
Journal of Translational Medicine volume 20, Article number: 537 (2022)
Multiple myeloma (MM) is a malignancy of plasma cells that remains incurable. Toll-like receptor 4 (TLR4) acts as a stress-responsive signal, protecting mitochondria during proteasome inhibitor (PI) exposure, maintaining mitochondrial metabolism and increasing drug resistance in MM. However, the mechanism of TLR4 regulation remains elusive.
The purpose of this study was to investigate the methylation pattern of multiple myeloma and its effect on the expression of HNRNPA2B1 and downstream targets.
The methylation level in MM and normal bone marrow specimens was detected using a colorimetric assay. HNRNPA2B1 gene knockdown was achieved in RPMI 8226 MM cells via adenovirus transfection. CCK8 and flow cytometric assays were used to detect proliferation and apoptosis, respectively. Transcriptome sequencing and m6A methylation MeRIP sequencing were applied, and differentially expressed genes (DEGs) were detected. Three independent NCBI GEO datasets were applied to examine the effects of HNRNPA2B1 and TLR4 expression on MM patient survival.
HNRNPA2B1 promoted MM progression. Clinical data from database revealed that HNRNPA2B1 was adverse prognostic factor for survival among MM patients. Furthermore, transcriptome sequencing and methylation sequencing showed that HNRNPA2B1 recognized and was enriched at the m6A sites of TLR4 and TLR4 was down-regulated of both the m6A level and transcription level in HNRNPA2B1-knockdown MM cells. Moreover, TLR4 was an adverse survival prognostic factor based on database analysis.
Overall, our study implies that the RNA-binding protein HNRNPA2B1 increases cell proliferation and deregulates cell apoptosis in MM through TLR4 signaling. Our study suggests HNRNPA2B1 as a potential therapeutic target for MM.
Multiple myeloma (MM) is a clonal B-cell malignancy that accounts for approximately 10% of all hematologic malignancies and is characterized by abnormal plasma cell proliferation and abnormal globulin secretion in bone marrow . Although survival outcomes in MM patients have significantly improved over the past 20 years, the disease remains incurable [2, 3]. Additionally, chemotherapy resistance develops often in MM patients during treatment, stemming from a variety of factors including genetic abnormalities . Therefore, new therapeutic targets are urgently needed. Transcriptomic studies have provided important information on disease-related pathways and critical genes in addition to revealing new targets for the treatment of diseases [5,6,7].
N6-methyladenosine (m6A) methylation is the most common post-transcriptional modification found on eukaryotic mRNA so far and also widely exists in a variety of bacteria and RNA viruses . m6A modification participates in the regulation of mRNA translation, splicing processing, as well as nuclear transport and degradation, which determines the entire life process of mRNA . Over the past few decades, the biological functions regulated by m6A have been shown to be involved in several types of tumor progression . However, the critical role of m6A modification in MM is still elusive, and the regulatory mechanisms of MM in general are not fully understood.
Heterogeneous nuclear ribonucleoprotein A2B1 (HNRNPA2B1) is one of the main m6A readers [11, 12]. The reader is a methylated reading protein that recognizes RNA methylation modification and performs its functions by specifically binding to the m6A-modified region or altering the RNA secondary structure to facilitate binding to RNA . Recently, abnormal HNRNPA2B1 expression has been discovered in various disease states and linked closely to tumor progression. HNRNPA2B1 may play key roles in cancer development due to its ability to accelerate pre-mRNA processing through the function of RNA binding .
In the present study, we discovered the hypermethylated status of MM and the high expression of HNRNPA2B1 in MM. Through transcriptome and m6A sequencing, we identified Toll-like receptor 4 (TLR4) as a target of HNRNPA2B1-mediated m6A modification. HNRNPA2B1 knockdown decreased the TLR4 mRNA level and simultaneously abolished TLR4 mRNA m6A modification. Through analysis of gene expression omnibus (GEO) datasets, HNRNPA2B1 and TLR4 were found to be adverse prognostic factors for survival among MM patients. Our results reveal the biological role of HNRNPA2B1 in mediating m6A modification in MM, indicating that HNRNPA2B1 plays an important role in MM pathogenesis via epigenetic regulation a critical target TLR4. From these findings, we propose that HNRNPA2B1 may be a novel therapeutic target for evaluating MM progression under m6A-based post-transcriptional regulation through TLR4 pathway.
Materials and methods
MM specimens and cell lines
Bone marrow specimens were obtained from 15 newly diagnosed MM patients and 15 normal controls (NC) at Harbin Medical University Cancer Hospital, China, between January 2019 and December 2020. The National Comprehensive Cancer Network clinical practice guidelines were used to diagnose MM.
The MM cell line, RPMI 8226 was purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). The cells were maintained in RPMI 1640 medium (Gibco, Grand Island, NY, USA) supplemented with 10% (for RPMI 8226) fetal bovine serum (FBS, Hyclone, Logan, UT, USA) and 1% penicillin–streptomycin at 37 ℃ in a humidified incubator containing 5% CO2.
RNA extraction, reverse transcription (RT) and quantitative PCR (qPCR)
Total RNA was isolated using Trizol Reagent (Invitrogen, USA). cDNA was synthesized from 1 µg of total RNA using the Prime Script™ RT reagent kit (Takara Bio, Inc., Otsu, Japan). Relative gene expression was assessed by real-time RT-PCR using the SYBR Premix Ex Taq TM kit (Takara Bio, Inc., Otsu, Japan) on the ABI 7900HT Real-Time PCR system (Applied Biosystems Life Technologies, Foster City, CA, USA) according to the manufacturers’ instructions and suggested protocol. The primers used for qRT-PCR are listed in Additional file 1: Table S1. Each analysis was carried out in triplicate.
Cell proliferation and apoptosis assays
Aliquots of 3 × 104 cells were suspended in 100 µl RPMI 1640 medium and then seeded into 96-well plates. The CCK8 (Dojindo Molecular Technologies, Japan) assay was conducted according to the manufacturer’s instructions. Briefly, after cell culture, 10 µl CCK8 solution was added into each 96-well, and after incubation for 2 h, the absorbance at 450 nm wavelength was measured. Each analysis was carried out in triplicate.
Apoptosis was evaluated by Annexin V-APC/propidium iodide (PI) staining. Briefly, cells were resuspended in binding buffer, and then Annexin V-APC and PI were added to each sample for 15 min and incubation at room temperature in darkness. Flow cytometry CantoII (BD Biosciences, San Jose, CA, USA) was used to evaluate the numbers of apoptotic cells. Each analysis was carried out in triplicate.
Adenovirus vector plasmid for transfection of RPMI 8226 cells
Adenovirus expressing HNRNPA2B1-shRNAs was purchased from Hanheng Bioscience Co., Ltd. (Shanghai, China). The HNRNPA2B1 shRNA sequences were: sh1,5′-TCGAGGCCATGGCTGCAAGACCTCATTCAATTTCAAGAGAATGAATGAGGTCTTGCAGCCATGGTTTTTTA-3′; sh2, 5′-AGCTTAAAAAACCATGGCTGCAAGACCTCATTCAATTCTCTTGAAATTGAATGAGGTCTTGCAGCCATGGCC-3′; and sh3, 5′-TCGAGGAGGAACAGTTCCGTAAGCTCTTTATTTCAAGAGAATAAAGAGCTTACGGAACTGTTCCTTTTTTTA-3′. Aliquots of 5 × 106 cells were plated onto 25-cm2 plates for 48 h before transfection with Ad-eGFP or Ad-HNRNPA2B1 with RPMI 1640 medium for 4 h at 37 °C. The cells were then washed twice with RPMI 1640 medium and once with RPMI 1640 medium containing 10% FBS, and then incubated for 2 days at 37 °C. After the formation of a stable cell line, RT-qPCR and western blot analyses were conducted to screen the RPMI 8226 with HNRNPA2B1 knockdown.
RNA m6A quantification
After extraction, the quality and quantity of the total RNA were measured by the NanoDrop (Thermo Fisher Scientific, Waltham, MA, USA). The EpiQuik m6A RNA Methylation Quantification Kit was then used to examine the m6A modification level according to the manufacturer’s instructions. Briefly, RNA was coated onto the assay wells, which were then washed and incubated with capture antibody. After antibody addition, the m6A level was detected by measurement of the absorbance at 450 nm (optical density [OD]450). Each analysis was carried out in triplicate.
To assess the RNA degradation and contamination, 1% agarose gel electrophoresis was conducted, and the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA) was used to test the RNA purity. RNA integrity was assessed by using the Bioanalyzer 2100 System (Agilent Technologies, CA, USA).
For RNA-sequencing of 13 paired samples, mRNA Capture Beads with Oligo(dT) (VAHTS, Nanjing, China) were first used to capture mRNA. mRNA was then purified with binding and washing buffer. mRNA was then randomly fragmented and reverse transcribed into cDNA in the corresponding buffer. Random primers were used to synthesize the first cDNA strand, and dNTPs/DNA polymerase I (ABclonal, MA, USA) was used to synthesize the second cDNA strand. After dA—tailing of 3' DNA fragments, UMI (Novogene, Beijing, China) and sequence adapter were used to ligate DNA. After PCR, the product was purified and collected using DNA Clean Beads (Beckman Coulter, CA, USA) and nuclease-free H2O. Each analysis was carried out in triplicate. After quality assessment, the prepared library was sequenced on the Illumina platform (Illumina Novaseq, Inc., USA).
M6A-RNA immunoprecipitation (MeRIP) assay and m6A sequencing
The mRNA m6A of 15 paired samples was sequenced by MeRIP-seq at Novogene (Beijing, China). The degradation and contamination of the extracted RNA were tested as described above. For immunoprecipitation, fragmented mRNA (~ 100 nt) was incubated with anti-m6A antibody at 4 ℃ for 2 h (Synaptic Systems). m6A modified-mRNA was then enriched and detected through qRT-PCR or next generation sequencing using Illumina HiSeq 2000 (Illumina Inc.). Each analysis was carried out in triplicate. The NEBNext ultraRNA library prepare kit for Illumina (Illumina, Inc., USA) was used for library construction. Finally, the library preparation was sequenced by an Illumina Novaseq platform.
Dot blot for measuring m6A
Nitrocellulose membrane was fixed in plates. Then 10 µg RNA per well was added to the dot bolt apparatus. Then RNA was crosslinked to the membrane using an ultraviolet (UV) crosslinker. The membrane was then washed with washing buffer and blocked with blocking buffer at room temperature for 30 min. After that, the membrane was incubated with anti-m6A antibody at 4 °C. Then the membrane was washed, incubated with western blotting reagent, and exposed to autoradiography.
NCBI GEO microarray dataset analysis
GEO datasets GSE116294, GSE80608 and GSE141260 were used to compare the expression of HNRNPA2B1 at different stage of monoclonal gammopathy. GSE116294 consisted of 69 clinical samples, including 4 normal samples, 50 MM samples at diagnosis and 15 plasma cell leukemia (PCL) samples at diagnosis. GSE80608 contained 10 control, 10 monoclonal gammopathy of undetermined significance (MGUS) and 10 MM samples. GSE141260 contained 10 control and 10 MM samples. Additionally, GSE4202, GSE24080 and GSE136400 were used to estimate the relationship between the expression of HNRNPA2B1/TLR4 and the disease prognosis of MM.
Based on the Affymetrix-GPL570 platform, the data for gene expression and survival were obtained from the GEO datasets GSE4204 (newly diagnosed MM, n = 538) and GSE24080 (newly diagnosed MM, n = 559). In GSE4204, the samples were divided into two subgroups based on the median gene expression values, the survival curves were drawn to analyze overall survival (OS) between different expression subgroups. In GSE24080, the clinical endpoints were OS and event-free survival (EFS) at 24 months, and we compared the expression levels of hnRNPA2B1 and TLR4 between two subgroups with different survival outcomes.
In addition, GSE136400 dataset was used to analyze the relationship between TLR4 expression and survival of MM patients. The dataset consisted of 1424 samples from different times, including before treatment, post-treatment but pre-transplantation, post-transplantation, pre-consolidation, etc., and we selected all samples before treatment (n = 354) for further analysis. The gene expression data were analyzed using the GPL27143 platform Affymetrix Human Genome U133 Plus 2.0 Array. X-tile software (v3.6.1, Yale University, New Haven, CT, USA) was used to confirm the cut-off value of TLR4 expression, and the patients were divided into two subgroups: TLR4high and TLR4low. Survival curves were drawn to analyze OS and progression-free survival (PFS) between the two subgroups.
Western blot analysis
After collection of cells and tissues, radioimmunoprecipitation assay buffer (Beyotime Institute of Biotechnology, Shanghai, China) was added into the pellet to extract proteins. After extraction, concentration of total protein was measured using a BCA protein assay kit (Beyotime Institute of Biotechnology). Agarose gel electrophoresis was performed after the addition of 2X sodium dodecyl sulfate loading buffer and incubated at 100 °C for 5 min. After transfer onto polyvinylidene difluoride membranes and blocking using 5% milk, the membranes were incubated with primary antibodies against TLR4 (1:500; Affinity Biosciences, Jiangsu, China), HNRNPA2B1 (1:500; Affinity Biosciences) and GAPDH (1:500; Affinity Biosciences) overnight at 4 °C. Then the secondary antibodies were applied to the membranes. After a series of washing steps, the membranes were visualized using chemiluminescence (Beyotime Institute of Biotechnology). Each analysis was carried out in triplicate.
Statistical analyses were conducted by GraphPad Prism 8 software (GraphPad Software, Inc., San Diego, CA, USA). Survival curves were drawn by the Kaplan–Meier approach, and the log-rank test was used for comparisons between groups. The differences between two groups were analyzed by unpaired t tests, and differences between multiple groups were analyzed by one-way analysis of variance (ANOVA). In all statistical analyses, p-values < 0.05 were considered statistically significant.
M6A reader HNRNPA2B1 is highly expressed in MM
To investigate the role of m6A modification in MM, we first conducted a dot blot assay to examine m6A methylation levels in bone marrow specimens from MM patients and normal controls. The m6A methylation level in MM bone marrow was significantly higher than that in the normal control bone marrow (p < 0.0001), with overall m6A average values for 15 MM patients and 15 normal controls of 0.246 ng and 0.029 ng, respectively (Fig. 1A).
As m6A is regulated by the multiple “writer”, “eraser” and “reader” genes, including METTL3, METTL14, HNRNPA2B1, WTAP, FTO, ALKBH5, YTHDF1/2/3, IGF2BP1/2/3, eIF3 and others , we next explored dysregulation of these genes in MM. qRT-PCR was performed to explore the expression of these genes in the RPMI 8226 MM cell line. The methylation “reader” gene HNRNPA2B1 was highly expressed in RPMI 8226 cells (Fig. 1B). Furthermore, the expression level of HNRNPA2B1 in bone marrow was significantly higher in 13 newly diagnosed patients with MM than in 13 normal controls (Fig. 1C).
In addition to MM cell line RPMI 8226 and bone marrow specimen, the elevated expression of HNRNPA2B1 in MM was also authenticated in GEO dataset. In GSE116294, the expression of HNRNPA2B1 did not differ among normal control, MM and PCL (P = 0.0795). However, as shown in Fig. 1D, the expression level of HNRNPA2B1 in MM patients was significantly higher compared with that of control patients (P = 0.0009). In GSE80608 and GSE141260, the expression of HNRNPA2B1 appeared to be higher in the MM group than in the normal group, but the difference did not reach statistical significance (Fig. 1E and F). Together, these results indicate that HNRNPA2B1 expression and m6A methylation level were upregulated in MM.
HNRNPA2B1 Knockdown decreases proliferation and increases apoptosis among MM cells
The sharp increase in HNRNPA2B1 expression in MM indicates its potential role in myeloma process. Thus, we investigated its influence on cell proliferation and apoptosis in MM. First, RPMI 8226 myeloma cells were transfected with the constructed HNRNPA2B1 viral plasmid (sh-HNRNPA2B1-8226, Fig. 2A), and the transfection efficiency was verified by qPCR. The expression of HNRNPA2B1 in the transfected RPMI 8226 cells was reduced by more than 50%, and the protein level of HNRNPA2B1 in these cells also was reduced (Fig. 2B and C). CCK8 assay results indicated that MM cell proliferation was significantly reduced at 24, 28, and 72 h after HNRNPA2B1 knockdown in a time-dependent manner (Fig. 2D). Additionally, flow cytometric detection of apoptotic MM cells showed that apoptosis among these cells was increased significantly at 72 h after HNRNPA2B1 knockdown, with the greatest increase in late apoptosis (Fig. 2E and F).
HNRNPA2B1 plays multiple roles in myeloma cells
To analyze the role of HNRNPA2B1 in myeloma cells, we compared the transcriptome sequencing results for MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226. The analysis identified 346 differentially expressed genes (Additional file 2: S1), of which 148 were up-regulated and 198 were down-regulated (Fig. 3A). Gene Ontology (GO) enrichment showed that the effect of HNRNPA2B1 in myeloma cells is mainly concentrated in the two processes of biological process and cell composition, and that HNRNPA2B1 participates in the regulation of myeloma cells (Fig. 3B). FPKM based on gene read count was used to quantify the differences among the two groups (sh-HNRNPA2B1-8226 and the control sh-NC-8226) for transcriptomics analysis. In the overall FPKM hierarchical clustering diagram, the Log10 (FPKM + 1) value is normalized to scale number and clustered. Genes expressed at high levels are indicated by red, and genes expressed at low levels are indicated by blue. The color from red to blue indicates that log10 (FPKM + 1) is from high to low. As shown in the heatmap, HNRNPA2B1 regulated multiple gene clusters in MM (Fig. 3C). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed the differential gene enrichment was related to a variety of signaling pathways including proteins processing in endoplasmic reticulum (ER) stress and hypoxia induction, suggesting that HNRNPA2B1 plays an important role in ER stress in the regulation of myeloma (Fig. 3D).
To analyze how HNRNPA2B1 regulates m6A methylation in myeloma cells, m6A methylation immunoprecipitation sequencing (MeRIP) was performed in MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226. The results showed that 496 peaks were differentially expressed in MeRIP (Additional file 3: S2), of which 98 were up-regulated and 398 were down-regulated (Fig. 4A). The volcano plot indicated the distribution of RPM based on the difference peak between the two groups. The black line represents the median value in each group of experiments. The dotted line represents the average value of the two groups (Fig. 4B). The heatmap of differential peak signal distribution showed that HNRNPA2B1 regulates the RPM signal (Fig. 4C). GO enrichment indicated that the significantly different peaks were enriched in a variety of regulatory signal pathways of biological processes, cell components, and molecular functions (Fig. 4D). KEGG enrichment analysis showed that the pathways of transcription regulation, ER stress, oxidative phosphorylation, MAPK signaling pathway, EBV infection, and cell adhesion molecules were influenced. These findings suggest that HNRNPA2B1 participates in the regulation of m6A methylation levels of related genes in multiple signaling pathways (Fig. 4E).
TLR4 is a downstream target of HNRNPA2B1
To analyze how HNRNPA2B1 regulates m6A methylation in myeloma cells, m6A methylation immunoprecipitation sequencing (MeRIP) was performed in MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226. The m6A consensus sequence (RRACH) motif was shown in Fig. 5A. An enriched peak was detected around the stop codon region of TLR4 mRNA in control RPMI 8226 cells, and the peak was diminished upon HNRNPA2B1 knockdown (Fig. 5B). The percentage of peak distribution in different regions (5'UTR, CDS and 3'UTR) of sh-HNRNPA2B1-8226 and the control cells were shown in Fig. 5C. The distribution histogram of peak on four transcriptional functional regions (5'UTR, 3'UTR, CDs and intronic) is shown in Fig. 5D. To identify the key target of HNRNPA2B1 in the process of m6A methylation regulation, we analyzed correlations in the results of transcriptome sequencing and MeRIP sequencing. Of these, 23 genes were differentially expressed in both the transcriptome and methylation (Additional file 4: S3), including PLEKHO1, KCNH2, DDX18, NRCAM, RASSF2, STS, BMF, SYNGR2, SPP1, AHNAK, APOL2, TLR4, DUSP5, BMP2K, ARHGAP17, HSPD1, MX1, PBXIP1, RNF187, RAPH1, HEXIM1, HLA-DRA and HSPA1B (Fig. 5E). Importantly, TLR4 was identified as a direct target of m6A methylation in the result of our m6A-Seq analysis. In addition, studies have confirmed that TLR4 participates in the processes of activation of multiple inflammatory pathways and tumorigenesis . Moreover, TLR4 was shown to be associated with the myeloma tumor microenvironment and bortezomib resistance . Therefore, we selected TLR4 as a candidate target of HNRNPA2B1-mediated m6A modification for further investigation.
Enrichment of TLR4 in the GO differential gene pathways was analyzed. TLR4 was enriched in a variety of important signaling pathways, such as alkylation metabolism, cell metabolism regulation, cell negative regulation metabolism, cell cycle complex metabolism regulation, etc. KEGG enrichment analysis showed that the pathways of transcription regulation, ER stress, oxidative phosphorylation, MAPK signaling pathway, EBV infection, and cell adhesion molecules were influenced (Fig. 5F). These findings suggest that HNRNPA2B1 participates in the regulation of m6A methylation levels of related genes in multiple signaling pathways.
From the KEGG differential gene pathway analysis, TLR4 enrichment was also involved in a variety of signaling pathways, such as the PI3K-Akt pathway, influenza virus A pathway, hepatitis B signaling pathway, phagocytic signaling pathway, Toll-like receptor signaling pathway, Salmonella infection signaling pathway, cancer proteoglycan, hypoxia inducible factor signaling pathway, inflammatory bowel disease signaling pathway, rheumatoid arthritis signaling pathway, ER stress, nuclear factor (NF)-κB signaling pathway, and many other important signaling pathways related to tumors and inflammation.
Next, qRT-PCR and western blot analyses were conducted to verify the relationship between TLR4 and HNRNPA2B1. We observed that knockdown of HNRNPA2B1 suppressed the mRNA and protein expression levels of TLR4 (Fig. 5G and H). Together, our data suggested that HNRNPA2B1 promotes tumor progression of MM through the regulation of TLR4 expression.
HNRNPA2B1 and TLR4 expression is associated with the survival of MM patients in NCBI GEO microarray datasets
Further, we used datasets to further explore the role of HNRNPA2B1 and TLR4 expression in clinical specimens of MM. Three independent NCBI GEO datasets were applied to validate the effect of HNRNPA2B1 and TLR4 expression on the survival of MM patients. We found that the OS of the HNRNPA2B1high subgroup was significantly poorer than that of the HNRNPA2B1low subgroup in GSE4204 (Fig. 6A, P = 0.0086). In GSE24080 set, the subgroups with EFS > 24 months and OS > 24 months showed lower expression of HNRNPA2B1compared to the subgroups with EFS ≤ 24 months and OS ≤ 24 months, (Fig. 6B and C, P = 0.0007, 0.0123, respectively), which was consistent with the results from GSE4204 and indicates that MM patients with higher HNRNPA2B1 expression had relatively decreased survival times.
With regard to TLR4 expression in MM patients, we found that the TLR4high subgroup had significantly poorer survival, including both OS and PFS, compared to the TLR4low subgroup in two independent datasets of GSE136400 (Fig. 6D and E, P = 0.0095, 0.0114, respectively) and GSE4204 (Fig. 6F, P = 0.0397). In GSE24080, TLR4 expression was higher in the OS ≤ 24 months subgroup than in the OS > 24 months subgroup (Fig. 6G, P = 0.023). These results identify TLR4 acted as an adverse prognostic factor of survival among MM patients. Moreover, the consistency rate of the HNRNPA2B1/TLR4 expression was about 46% in the GSE4204 dataset, with 123 patients (22.8%) having both high HNRNPA2B1 and TLR4 expression, and 125 patients (23.2%) having both low HNRNPA2B1 and TLR4 expression. Altogether, the above results indicate a possible role of HNRNPA2B1 in the positive regulation of TLR4 in MM patient through m6A methylation.
M6A is the most common RNA chemical modification in human mRNA. Recently, with the advancement of high-throughput sequencing and antibody-based technology, researchers have been able to accurately map the site of m6A and further study its biological functions . Although studies have shown potential functions of m6A in regulating mRNA decay, translation, and processing, the pathologic significance of m6A in MM has remained unknown. In the present study, we found that m6A levels were significantly elevated in MM due to upregulation of the methyl reader HNRNPA2B1. Our observations in MM cells with HNRNPA2B1 knockdown indicated that HNRNPA2B1 plays a critical role in promoting MM proliferation and inhibiting MM apoptosis. In addition, we found that TLR4 is a downstream target of m6A modification mediated by HNRNPA2B1. In summary, our results suggest that HNRNPA2B1 and its associated m6A modification promote MM progression by regulating the expression of the key gene TLR4 at the post-transcriptional level. Therefore, our findings reveal a new epigenetic regulatory mechanism contributing to the progression of MM.
HNRNPA2B1 participates in multiple biological processes in many diseases, especially cancers. One study demonstrated that down-regulation of HNRNPA2B1 can inhibit the cell proliferation, tumor invasion, and cell cycle progression through the PI3K/AKT signaling pathway of cervical cancer cells, thus triggering cell apoptosis . HNRNPA2B1 also was identified as an oncogene in head and neck cancer, which promotes epithelial to mesenchymal transition through the AKT/PKB signaling pathway . Additionally, the protein level of HNRNPA2B1 can be modulated by its ubiquitination status through miR503HG, thereby regulating HCC metastasis and migration in hepatocellular carcinoma (HCC) . Recently, the role of HNRNPA2B1 in MM was also explored by Jiang and colleagues, who found that HNRNPA2B1 expression is a negative prognostic factor in MM. Mechanistically, interleukin enhancer binding factor 3 (ILF3) was identified as an m6A target site of HNRNPA2B1 . In the present study, we identified TLR4 as the HNRNPA2B1 target site. Notably, we did not observe ILF3 among the core genes identified by RNA-sequencing and MeRIP sequencing. These inconsistent results may be due to different experiment conditions, as well as slight differences in the sequencing technologies employed.
TLR4, as a pattern recognition receptor, mainly mediates endogenous immunity and immune presentation, regulating the inflammatory response and participating in the expression of inflammatory factors. Recent studies have suggested that TLR4 expression correlates with cancer progression and TLR4 plays a key role in disrupting tumor cell apoptosis regulation, decreasing tumor cell apoptosis, and promoting tumor cell proliferation. For example, TLR4 activation releases immunosuppressive exosomes, promotes tumor progression, and accelerates the metastatic process . In MM, TLR4 expression was increased in the bone marrow cells of MM patients compared to healthy volunteers. Flow cytometric analyses also indicated upregulation of TLR4 in MM cells . Mechanistically, TLR4 suppresses ER stress-related apoptosis and promotes MM cell proliferation and survival through the PERK-CHOP pathway . A recent study highlighted the overexpression of TLR4 according to disease stage and implied that TLR4 inhibition reduced mesenchymal stromal cell activity and decreased MM cell growth . Consistent with these previous studies, our results also imply that a high level of TLR4 expression correlates with poor prognosis in MM patients.
m6A methylation and its regulators links epigenetic transcriptomics to tumor development. Studies have shown that m6A modification plays critical and diverse roles in the progression of various cancers . m6A methylation and its enzyme METTL3 are upregulated in glioblastoma stem cells and HCC, and the high expression of METTL3 induces SOCS2 degradation by m6A methylation, thus promoting tumor cell growth and survival . In another study, fat mass and obesity-associated protein (FTO), an m6A eraser, showed altered expression in intrahepatic cholangiocarcinoma, and low FTO expression predicted poor prognosis . These studies imply oncogenic functions of m6A methylation. On the other hand, several studies have shown that m6A may work as a tumor suppressor in some cancers. For example, one study found that METTL4 down-regulation led to m6A modification dysregulation, which was related to tumor metastasis and poor prognosis in HCC . In our study, we suggested that m6A modification has an oncogenic function in MM. Dysregulation of the m6A level seems to act as a “double-edged sword” in cancer progression, and further studies are needed to elucidate the key role of m6A.
The findings of this study must be considered in light of some limitations. According to the results in Fig. 1B, HNRNPA2B1 and eukaryotic translation initiation factor 3 (eIF3, B, D, E, K) were both overexpressed genes. eIF3 is the most complex eukaryotic translation initiation factor and is crucial for tumorigenesis [30,31,32]. eIF3 regulates cap-dependent translation processes and also plays an important role in cap-dependent translation regulation by binding to putative internal ribosome entry sites [33, 34]. The role of eIF3 in MM is not fully understood, and this can be a future direction of research. eIF3 may also represent a potential mechanism for the regulation of TLR4 in MM. However, our study did not include experiments to exclude it. Further exploration of the role of eIF3 in MM is needed in the future.
In summary, the results of the present study indicate that HNRNPA2B1 acts as oncogene in MM development. Our experiments suggested the HNRNPA2B1 promotes MM proliferation and inhibits apoptosis through the epigenetic regulation of TLR4 by m6A modification. Moreover, HNRNPA2B1 expression was significantly increased in MM and shown to correlate with poor prognosis in MM patients. Altogether, our results suggest HNRNPA2B1 as a potential therapeutic target for MM.
The data analyzed during this study are included in the published article and supplementary materials. No more additional data is generated.
Differentially expressed genes
Fat mass and obesity-associated protein
Gene expression omnibus
Heterogeneous nuclear Ribonucleoprotein A2B1
Interleukin enhancer binding factor 3
Kyoto Encyclopedia of Genes and Genomes
Methylation immunoprecipitation sequencing
Nuclear factor -κB
Toll-like receptor 4
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin. 2020;70:7–30. https://doi.org/10.3322/caac.21590.
Xiao CL, Zhang ZP, Xiong S, Lu CH, Wei HP, Zeng HL, et al. Comparative proteomic analysis to discover potential therapeutic targets in human multiple myeloma. Proteomics Clin Appl. 2009;3:1348–60. https://doi.org/10.1002/prca.200900068.
Rajkumar SV. Multiple myeloma: every year a new standard? Hematol Oncol. 2019;37(Suppl 1):62–5. https://doi.org/10.1002/hon.2586.
Kazandjian D. Multiple myeloma epidemiology and survival: a unique malignancy. Semin Oncol. 2016;43:676–81. https://doi.org/10.1053/j.seminoncol.2016.11.004.
Cavill R, Jennen D, Kleinjans J, Briedé JJ. Transcriptomic and metabolomic data integration. Brief Bioinform. 2016;17:891–901. https://doi.org/10.1093/bib/bbv090.
Carrasco DR, Tonon G, Huang Y, Zhang Y, Sinha R, Feng B, et al. High-resolution genomic profiles define distinct clinico-pathogenetic subgroups of multiple myeloma patients. Cancer Cell. 2006;9:313–25. https://doi.org/10.1016/j.ccr.2006.03.019.
Broyl A, Hose D, Lokhorst H, de Knegt Y, Peeters J, Jauch A, et al. Gene expression profiling for molecular classification of multiple myeloma in newly diagnosed patients. Blood. 2010;116:2543–53. https://doi.org/10.1182/blood-2009-12-261032.
Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3’ UTRs and near stop codons. Cell. 2012;149:1635–46. https://doi.org/10.1016/j.cell.2012.05.003.
Zhao BS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Biol. 2017;18:31–42. https://doi.org/10.1038/nrm.2016.132.
Tuncel G, Kalkan R. Importance of m N(6)-methyladenosine (m(6)A) RNA modification in cancer. Med Oncol. 2019;36:36. https://doi.org/10.1007/s12032-019-1260-6.
Alarcón CR, Goodarzi H, Lee H, Liu X, Tavazoie S, Tavazoie SF. HNRNPA2B1 Is a mediator of m(6)A-dependent nuclear RNA processing events. Cell. 2015;162:1299–308. https://doi.org/10.1016/j.cell.2015.08.011.
Dreyfuss G, Matunis MJ, Piñol-Roma S, Burd CG. hnRNP proteins and the biogenesis of mRNA. Annu Rev Biochem. 1993;62:289–321. https://doi.org/10.1146/annurev.bi.62.070193.001445.
Torres IO, Fujimori DG. Functional coupling between writers, erasers and readers of histone and DNA methylation. Curr Opin Struct Biol. 2015;35:68–75. https://doi.org/10.1016/j.sbi.2015.09.007.
Dai S, Zhang J, Huang S, Lou B, Fang B, Ye T, et al. HNRNPA2B1 regulates the epithelial-mesenchymal transition in pancreatic cancer cells through the ERK/snail signalling pathway. Cancer Cell Int. 2017;17:12. https://doi.org/10.1186/s12935-016-0368-4.
Shetab Boushehri MA, Lamprecht A. TLR4-based immunotherapeutics in cancer: a review of the achievements and shortcomings. Mol Pharm. 2018;15:4777–800. https://doi.org/10.1021/acs.molpharmaceut.8b00691.
Giallongo C, Tibullo D, Puglisi F, Barbato A, Vicario N, Cambria D, et al. Inhibition of TLR4 signaling affects mitochondrial fitness and overcomes bortezomib resistance in myeloma plasma cells. Cancers. 2020. https://doi.org/10.3390/cancers12081999.
Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485:201–6. https://doi.org/10.1038/nature11112.
Shi X, Ran L, Liu Y, Zhong SH, Zhou PP, Liao MX, et al. Knockdown of hnRNP A2/B1 inhibits cell proliferation, invasion and cell cycle triggering apoptosis in cervical cancer via PI3K/AKT signaling pathway. Oncol Rep. 2018;39:939–50. https://doi.org/10.3892/or.2018.6195.
Gupta A, Yadav S, Pt A, Mishra J, Samaiya A, Panday RK, et al. The HNRNPA2B1-MST1R-Akt axis contributes to epithelial-to-mesenchymal transition in head and neck cancer. Lab Investig. 2020;100:1589–601. https://doi.org/10.1038/s41374-020-0466-8.
Wang H, Liang L, Dong Q, Huan L, He J, Li B, et al. Long noncoding RNA miR503HG, a prognostic indicator, inhibits tumor metastasis by regulating the HNRNPA2B1/NF-κB pathway in hepatocellular carcinoma. Theranostics. 2018;8:2814–29. https://doi.org/10.7150/thno.23012.
Jiang F, Tang X, Tang C, Hua Z, Ke M, Wang C, et al. HNRNPA2B1 promotes multiple myeloma progression by increasing AKT3 expression via m6A-dependent stabilization of ILF3 mRNA. J Hematol Oncol. 2021;14:54. https://doi.org/10.1186/s13045-021-01066-6.
Domenis R, Cifù A, Marinò D, Fabris M, Niazi KR, Soon-Shiong P, et al. Toll-like receptor-4 activation boosts the immunosuppressive properties of tumor cells-derived exosomes. Sci Rep. 2019;9:8457. https://doi.org/10.1038/s41598-019-44949-y.
Xu Y, Zhao Y, Huang H, Chen G, Wu X, Wang Y, et al. Expression and function of toll-like receptors in multiple myeloma patients: toll-like receptor ligands promote multiple myeloma cell growth and survival via activation of nuclear factor-kappaB. Br J Haematol. 2010;150:543–53. https://doi.org/10.1111/j.1365-2141.2010.08284.x.
Bagratuni T, Sklirou AD, Kastritis E, Liacos CI, Spilioti C, Eleutherakis-Papaiakovou E, et al. Toll-like receptor 4 activation promotes multiple myeloma cell growth and survival via suppression of the endoplasmic reticulum stress factor chop. Sci Rep. 2019;9:3245. https://doi.org/10.1038/s41598-019-39672-7.
Lemaitre L, Hamaidia M, Descamps JG, Do Souto Ferreira L, Joubert MV, Gadelorge M, et al. Toll-like receptor 4 selective inhibition in medullar microenvironment alters multiple myeloma cell growth. Blood Adv. 2021. https://doi.org/10.1182/bloodadvances.2020003704.
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. https://doi.org/10.1186/s12943-019-1109-9.
Chen M, Wei L, Law CT, Tsang FH, Shen J, Cheng CL, et al. RNA N6-methyladenosine methyltransferase-like 3 promotes liver cancer progression through YTHDF2-dependent posttranscriptional silencing of SOCS2. Hepatology. 2018;67:2254–70. https://doi.org/10.1002/hep.29683.
Rong ZX, Li Z, He JJ, Liu LY, Ren XX, Gao J, et al. Downregulation of fat mass and obesity associated (FTO) promotes the progression of intrahepatic cholangiocarcinoma. Front Oncol. 2019;9:369. https://doi.org/10.3389/fonc.2019.00369.
Ma JZ, Yang F, Zhou CC, Liu F, Yuan JH, Wang F, et al. METTL14 suppresses the metastatic potential of hepatocellular carcinoma by modulating N(6) -methyladenosine-dependent primary MicroRNA processing. Hepatology. 2017;65:529–43. https://doi.org/10.1002/hep.28885.
Unbehaun A, Borukhov SI, Hellen CU, et al. Release of initiation factors from 48S complexes during ribosomal subunit joining and the link between establishment of codon-anticodon base-pairing and hydrolysis of eIF2-bound GTP. Genes Dev. 2004;18(24):3078–93. https://doi.org/10.1101/gad.1255704.
Dong Z, Zhang JT. Initiation factor eIF3 and regulation of mRNA translation, cell growth, and cancer. Crit Rev Oncol Hematol. 2006;59(3):169–80. https://doi.org/10.1016/j.critrevonc.2006.03.005.
Yin JY, Dong Z, Liu ZQ, et al. Translational control gone awry: a new mechanism of tumorigenesis and noveltargets of cancer treatments. Biosci Rep. 2011;31(1):1–15. https://doi.org/10.1042/BSR20100077.
Wei ZY, Zhang P, Zhou ZC, et al. Crystal structure of human eIF3k, the first structure of eIF3 subunit. J Biol Chem. 2004;279(33):34983–90. https://doi.org/10.1074/jbc.M405158200.
Mayeur GL, Hershey JW. Malignant transformation by the eukaryotic translation initiation factor 3 subunit p48 (eIF3e). FEBS Lett. 2002;514(1):49–54. https://doi.org/10.1016/s0014-5793(02)02307-4.
We thank Novogene (Beijing, China) and GeneCreate Biological Engineering Co. Ltd. (Wuhan, China) for helping with RNA sequencing and MeRIP sequencing.
This study provides the first evidence of the mechanism by which m6A reader HNRNPA2B1 regulates m6A methylation in myeloma. RNA-sequencing and MeRIP sequencing were used to analyze the mechanism, and HNRNPA2B1 and TLR4 were revealed to be adverse prognostic factors for survival in myeloma patients.
This study was supported by the National Natural Science Foundation of China (81870151), the China National Cancer Center Pandeng Project (NCC201908A01), Beijing Science and Technology Innovation Medical Development Foundation (award no. KC2021-JF-0055-04), Harbin Medical University Cancer Hospital Haiyan Foundation (JJZD2023-07) and Haiyan Foundation (JJMS2023-11).
Ethics approval and consent to participate
The study was approved by the Scientific Ethics Committee of the Harbin Medical University Cancer Hospital, China. Informed consent (ICF) was obtained from each participant before the start of the study.
Consent for publication
All the authors have read and approved the final manuscript.
The authors declare no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original version of this article was revised: authors’ sequence and affiliations’ assignment were incorrectly published.
Primers used for qRT-PCR.
Transcriptome sequencing identified 346 differentially expressed genes for MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226.
MeRIP sequencing identified 496 differentially peaks for MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226.
The 23 genes differentially expressed in both the transcriptome and methylation sequencing for MM cells transfected with sh-HNRNPA2B1-8226 and the control sh-NC-8226.
About this article
Cite this article
Jia, C., Guo, Y., Chen, Y. et al. HNRNPA2B1-mediated m6A modification of TLR4 mRNA promotes progression of multiple myeloma. J Transl Med 20, 537 (2022). https://doi.org/10.1186/s12967-022-03750-8
- Multiple myeloma
- N6-methyladenosine methylation