BPG is committed to discovery and dissemination of knowledge
Basic Study Open Access
Copyright: ©Author(s) 2026. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution-NonCommercial (CC BY-NC 4.0) license. No commercial re-use. See permissions. Published by Baishideng Publishing Group Inc.
World J Gastroenterol. Nov 7, 2026; 32(41): 121890
Published online Nov 7, 2026. doi: 10.3748/wjg.121890
Epitranscriptomic regulation of the COX-2/PGE2 axis by hnRNPA2B1 promotes esophageal squamous cell carcinoma progression
Jiao-Qian Lu, Xue-Fen Fang, Xiao-Xiong Guo, Can-Mei Zhong, Qiu-Yan Lin, Xiao-Yun Huang, Xiao-Zhong Wang, Feng-Lin Chen, Bi-Yun Zheng, Department of Gastroenterology, Fujian Medical University Union Hospital, Fuzhou 350001, Fujian Province, China
Jiao-Qian Lu, Xue-Fen Fang, Xiao-Xiong Guo, Hua-Long Zheng, Can-Mei Zhong, Qiu-Yan Lin, Xiao-Yun Huang, Xiao-Zhong Wang, Feng-Lin Chen, Bi-Yun Zheng, Department of Gastroenterology, Fujian Clinical Research Center for Digestive System Tumors and Upper Gastrointestinal Disease, Fuzhou 350001, Fujian Province, China
Miao Liu, Department of Gastrointestinal Endoscopy Nursing, Fujian Medical University Union Hospital, Fuzhou 350001, Fujian Province, China
Ji-Hong Lin, Department of Thoracic Surgery, Fujian Medical University Union Hospital, Fuzhou 350001, Fujian Province, China
Hua-Long Zheng, Department of Gastric Surgery, Fujian Medical University Union Hospital, Fuzhou 350001, Fujian Province, China
Hao Tian, Eye Center, The Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310009, Zhejiang Province, China
Jian-Wen Que, Department of Medicine, Columbia University Irving Medical Center, New York, NY 10032, United States
ORCID number: Bi-Yun Zheng (0000-0001-7680-8661).
Co-first authors: Jiao-Qian Lu and Xue-Fen Fang.
Co-corresponding authors: Feng-Lin Chen and Bi-Yun Zheng.
Author contributions: Lu JQ and Fang XF contributed equally as co-first authors; Zheng BY and Chen FL conceptualized and supervised the study as co-corresponding authors; Lu JQ, Fang XF, Guo XX, Liu M, Zheng HL, Zhong CM, Lin QY, Huang XY, Wang XZ, Tian H, Que JW, Chen FL, and Zheng BY developed the methodology; Lu JQ, Fang XF, Guo XX, Liu M, Lin JH, Zheng HL, and Zhong CM performed the investigation; Lu JQ, Fang XF, Guo XX, Liu M, and Zheng BY conducted formal analysis and curated the data; Chen FL, Que JW, and Zheng BY provided resources; Lu JQ and Zheng BY drafted the manuscript; Que JW, Chen FL, and Zheng BY reviewed and edited the manuscript; all authors approved the final version of the article.
AI contribution statement: AI tools were not used to generate scientific data, perform data analysis, or draw scientific conclusions. The authors take full responsibility for the integrity, accuracy, and scientific content of the manuscript.
Supported by National Natural Science Foundation of China Young Scientists Fund, No. 82203308; National Natural Science Foundation of China General Program, No. 82570660; Fujian Provincial Natural Science Foundation of China, No. 2026J010043; Fujian Provincial Natural Science Foundation of China, No. 2025J01116; Fujian Natural Science Foundation, No. 2023J01651; Young and Middle-Aged Talents Training Project of Fujian Provincial Health Commission, No. 2023GGA018; and National Key Clinical Specialty Construction Project of Fujian Province, China, No. 2023-1594.
Institutional review board statement: The study was reviewed and approved by the Institutional Ethics Committee of Fujian Medical University Union Hospital, No. 2022(85).
Institutional animal care and use committee statement: All animal experiments were reviewed and approved by the Institutional Animal Care and Use Committee of Fujian Medical University, No. 2022-NSFC-0403.
Conflict-of-interest statement: All the authors report no relevant conflicts of interest for this article.
ARRIVE guidelines statement: The authors have read the ARRIVE guidelines, and the manuscript was prepared and revised according to the ARRIVE guidelines.
Data sharing statement: Bulk RNA-Seq data are available on reasonable request. All other data relevant to the study are included in the article or uploaded as online supplemental information.
Corresponding author: Bi-Yun Zheng, MD, PhD, Associate Chief Physician, Associate Professor, Department of Gastroenterology, Fujian Medical University Union Hospital, No. 29 Xinquan Road, Fuzhou 350001, Fujian Province, China. drzhengby@163.com
Received: April 7, 2026
Revised: June 30, 2026
Accepted: July 20, 2026
Published online: November 7, 2026
Processing time: 167 Days and 15.7 Hours

Abstract
BACKGROUND

Esophageal squamous cell carcinoma (ESCC) remains a highly lethal malignancy with limited therapeutic options. Cyclooxygenase-2 (COX-2) is frequently overexpressed in ESCC, yet the mechanisms leading to its dysregulation remain poorly understood. Heterogeneous nuclear ribonucleoprotein A2/B1 (hnRNPA2B1), an oncogenic N6-methyladenosine (m6A) reader, has been implicated in various cancers, but its role in regulating COX-2 expression and promoting ESCC progression has not been defined.

AIM

To investigate how hnRNPA2B1 promotes ESCC progression through COX-2 regulation.

METHODS

Expression of hnRNPA2B1 and COX-2 was analyzed in ESCC tissues and cell lines by reverse transcription-quantitative polymerase chain reaction (qPCR), western blotting, and immunohistochemistry. Stable hnRNPA2B1 knockdown and overexpression models were established to assess cellular proliferation, migration, invasion, and xenograft tumor growth. Bulk RNA-seq, m6A RNA immunoprecipitation-qPCR and RNA immunoprecipitation-qPCR, and actinomycin D assays were used to delineate post-transcriptional regulation of COX-2. Functional and therapeutic relevance of the hnRNPA2B1-COX-2/prostaglandin E2 (PGE2) axis was evaluated using PGE2 supplementation and celecoxib treatment in vitro and in vivo.

RESULTS

HnRNPA2B1 was significantly upregulated in ESCC and early intraepithelial neoplasia, correlating with advanced stage and poor survival. Functionally, hnRNPA2B1 promoted proliferation, colony formation, migration, invasion, and xenograft growth, whereas its silencing had opposite effects. Transcriptomic and mechanistic analyses identified COX-2 as a downstream effector. HnRNPA2B1 bound to m6A-modified COX-2 transcripts, stabilized its mRNA, and enhanced COX-2/PGE2 axis. Exogenous PGE2 supplementation rescued the inhibitory effects of hnRNPA2B1 knockdown, while celecoxib treatment suppressed hnRNPA2B1-driven malignant phenotypes both in vitro and in vivo.

CONCLUSION

HnRNPA2B1 drives ESCC initiation and progression by epitranscriptomically activating the COX-2/PGE2 axis. Targeting COX-2 or hnRNPA2B1-mediated m6A stabilization may offer a therapeutic strategy for ESCC.

Key Words: Esophageal squamous cell carcinoma; Heterogeneous nuclear ribonucleoprotein A2/B1; Cyclooxygenase-2; Prostaglandin E2; N6-methyladenosine

Core Tip: Heterogeneous nuclear ribonucleoprotein A2/B1 (HnRNPA2B1) is upregulated in esophageal squamous cell carcinoma and early intraepithelial neoplasia and promotes tumor progression. Mechanistically, hnRNPA2B1 stabilizes cyclooxygenase-2 mRNA in an m6A-related manner, leading to activation of the cyclooxygenase-2/prostaglandin E2 axis. Prostaglandin E2 supplementation rescues the inhibitory effects of hnRNPA2B1 silencing, whereas celecoxib suppresses hnRNPA2B1-driven malignant phenotypes. These findings reveal a novel epitranscriptomic mechanism of esophageal squamous cell carcinoma progression and highlight hnRNPA2B1 as a potential therapeutic target.



INTRODUCTION

Esophageal squamous cell carcinoma (ESCC) is one of the most prevalent and lethal malignancies worldwide, particularly in East Asia. Despite advances in surgery, radiotherapy, chemotherapy, and immunotherapy, the overall survival of patients with ESCC remains dismal, primarily due to late diagnosis, high recurrence rates, and resistance to therapy[1,2]. A deeper understanding of the molecular mechanisms driving ESCC progression is therefore essential to identify effective biomarkers and therapeutic targets.

N6-methyladenosine (m6A) is the most abundant internal modification of eukaryotic mRNA, influencing RNA splicing, stability, translation, and degradation[3,4]. Dysregulation of m6A regulators has been increasingly linked to the initiation and progression of various cancers[5]. m6A modification is dynamically regulated by methyltransferase “writers” (METTL3/METTL14), demethylase “erasers” (FTO, ALKBH5), and binding “readers”, including YTHDF proteins, IGF2BPs, and heterogeneous nuclear ribonucleoprotein A2/B1 (hnRNPA2B1)[6,7].

HnRNPA2B1, originally identified as an RNA-binding protein involved in alternative splicing and microRNA processing, has recently been recognized as an m6A reader that modulates transcript stability and oncogenic signaling[8,9]. Increasing evidence indicates that hnRNPA2B1 plays a pivotal role in multiple tumor types[10-12]. Pan-cancer analyses have reported elevated hnRNPA2B1 expression in ESCC, correlating with advanced stage and poor prognosis[13,14]. Functional studies suggest that hnRNPA2B1 promotes malignant progression by regulating lipid metabolism and the miR-17-92 cluster[15,16]. However, the downstream effectors and pathways through which hnRNPA2B1 drives ESCC progression remain largely unclear.

In this study, we found that the levels of hnRNPA2B1 were increased in ESCC. We further demonstrated that hnRNPA2B1 regulated the mRNA stability of cyclooxygenase-2 (COX-2), a rate-limiting enzyme in prostaglandin biosynthesis[17] in an m6A-related manner. Prostaglandin E2 (PGE2) is one of the major downstream prostaglandins produced through COX-2 activity and has been implicated in tumor cell proliferation, invasion, inflammation, and immune modulation. Celecoxib, a selective COX-2 inhibitor, is widely used to suppress COX-2/PGE2 axis and was therefore used in this study to evaluate the functional relevance of this axis. Our in vitro and in vivo studies revealed that manipulation of the hnRNPA2B1-COX-2/PGE2 axis inhibited ESCC progression.

MATERIALS AND METHODS
Clinical specimens and tissue microarrays

ESCC tissues, matched adjacent tissues, and distal normal tissues were obtained from patients undergoing surgical resection at Fujian Medical University Union Hospital. High-grade intraepithelial neoplasia (HGIN) and low-grade intraepithelial neoplasia (LGIN) samples were collected from patients who underwent endoscopic submucosal dissection (ESD). None of the patients had received preoperative chemotherapy or radiotherapy. All tissue samples were reviewed and confirmed by two experienced pathologists. A tissue microarray containing 30 ESCC cases with matched adjacent and distant normal tissues was purchased from Outdo Biotech Co., Ltd. (Shanghai, China). Written informed consent was obtained from all patients, and the study protocol was reviewed and approved by the Institutional Ethics Committee of Fujian Medical University Union Hospital, No. 2022(85). All procedures were conducted in accordance with the Declaration of Helsinki.

Bioinformatic analysis

RNA-seq data and clinical information for esophageal carcinoma (ESCA) patients were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/projects/TCGA) and Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/) databases. Expression differences of hnRNPA2B1 and prostaglandin-endoperoxide synthase 2 (PTGS2) between tumor and normal tissues, correlations with clinicopathological parameters, and associations with overall survival were analyzed using Gene Expression Profiling Interactive Analysis (http://gepia.cancer-pku.cn/) and Kaplan-Meier Plotter (http://kmplot.com).

Bulk RNA-seq and differential expression analysis

Total RNA was isolated from Eca109 cells with stable hnRNPA2B1 knockdown and corresponding control cells using TRIzol reagent (Invitrogen, Thermo Fisher Scientific, Waltham, MA, United States). RNA purity and integrity were assessed using a NanoDrop spectrophotometer and an Agilent 2100 Bioanalyzer, and samples with RNA integrity number ≥ 7.0 were selected for library construction. Polyadenylated mRNA was enriched from 1 μg of total RNA, fragmented, and used for strand-specific complementary DNA library preparation with the TruSeq Stranded mRNA Library Prep Kit (Illumina, San Diego, CA, United States) according to the manufacturer’s instructions. Indexed libraries were quantified, pooled, and sequenced on an Illumina NovaSeq 6000 platform to generate 150-bp paired-end reads.

Adaptor sequences, low-quality reads, and reads containing more than 10% ambiguous bases were removed during quality control. Clean reads were aligned to the human reference genome (GRCh38) using HISAT2, and gene-level read counts were generated with featureCounts. Differential expression analysis between knockdown and control groups was performed using DESeq2 with default normalization settings, and genes with an adjusted P value < 0.05 were defined as differentially expressed. Hierarchical clustering and volcano plots were generated based on normalized gene expression values. Functional annotation of differentially expressed genes was performed using clusterProfiler, including Gene Ontology enrichment analysis for biological processes.

Cell culture and lentiviral transduction

Human ESCC cell lines (Eca109 and KYSE150) and the immortalized esophageal epithelial cell line Het-1A were cultured in RPMI-1640 medium (Gibco, Thermo Fisher Scientific, Waltham, MA, United States) supplemented with 10% fetal bovine serum (Gibco, Thermo Fisher Scientific, Waltham, MA, United States), 100 U/mL penicillin, and 100 μg/mL streptomycin at 37 °C in a humidified incubator with 50 mL/L CO2. All cell lines were routinely tested and confirmed to be free of mycoplasma contamination. Cell identities were verified by short tandem repeat profiling. Cells were used at low passage numbers for all experiments.

Lentiviruses carrying hnRNPA2B1 short hairpin RNAs (sh-hnRNPA2B1-56-1, sh-hnRNPA2B1-57-12, and sh-hnRNPA2B1-58-1; for clarity, these are referred to as sh56-1, sh57-12, and sh58-1 throughout the manuscript) and the corresponding control (shCON), as well as lentiviruses for hnRNPA2B1 overexpression (OE) and its negative control (NC), were purchased from GeneChem Co., Ltd. (Shanghai, China). Cells were transduced at a multiplicity of infection of 20 and selected with puromycin (2 μg/mL) for 7 days. Infection efficiency was assessed by fluorescence microscopy, and knockdown or overexpression efficiency was validated by western blotting and reverse transcription-quantitative polymerase chain reaction (RT-qPCR).

Western blotting

Total protein was extracted using RIPA lysis buffer (Beyotime Biotechnology, Shanghai, China) containing protease and phosphatase inhibitors. Protein concentration was determined by bicinchoninic acid assay, and equal amounts were separated by sodium-dodecyl sulfate gel electrophoresis and transferred to polyvinylidene fluoride membranes. Membranes were blocked with 5% nonfat milk and incubated overnight at 4 °C with the following primary antibodies: Anti-hnRNPA2B1 (1:1000; Abcam, Cambridge, United Kingdom; Cat# ab259894), anti-COX-2 (1:1000; Cell Signaling Technology, Danvers, MA, United States; Cat# 12282S), anti-GAPDH (1:5000; Proteintech, Wuhan, China; Cat# 10494-1-AP), and anti-β-actin (1:5000; Proteintech, Wuhan, China; Cat# 20536-1-AP). After washing, membranes were incubated with HRP-conjugated secondary antibodies and visualized using enhanced chemiluminescence substrate (Clarity™ Western ECL Substrate; Bio-Rad, Cat#1705061). Band intensity was quantified by ImageJ software (version 1.53, National Institutes of Health, Bethesda, MD, United States). The full-length western blots are provided in Supplementary Figure 1.

RT-qPCR

Total RNA was extracted using TRIzol reagent, reverse-transcribed with PrimeScript RT kit (Takara Bio Inc., Shiga, Japan), and amplified with SYBR Green Master Mix (Takara Bio Inc., Shiga, Japan). Relative expression was calculated by the 2-ΔΔCt method using GAPDH as the internal control. Primer sequences are listed in Supplementary Table 1.

Cell proliferation and colony formation assays

Cell proliferation was assessed using the cell counting kit-8 (Biosharp, Hefei, China). Briefly, 2 × 103 cells were seeded in 96-well plates and incubated for up to 5 days. Absorbance at 450 nm was measured daily. For colony formation assays, 1000 cells were seeded in 6-well plates and cultured for 14 days. Colonies were fixed with methanol, stained with 0.1% crystal violet, photographed, and counted.

For all in vitro functional assays, cells were treated with PGE2 (1 μmol/L; Selleck Chemicals, Houston, TX, United States; Cat# S3003)[18,19] or celecoxib (10 μmol/L; Selleck Chemicals, Houston, TX, United States; Cat# S1261)[20,21] as indicated, with 0.1% dimethyl sulfoxide (DMSO) serving as the vehicle control. Drug treatments were initiated 12 hours after cell seeding to allow cell attachment and were maintained throughout the experimental period unless otherwise specified.

Migration and invasion assays

For wound-healing assays, confluent cells were scratched with a sterile pipette tip, washed to remove debris, and incubated in serum-free medium. Wound closure was imaged at 0, 12, and 24 hours under a microscope to assess wound closure. Cells were treated with PGE2 (1 μM) or celecoxib (10 μM) as indicated, with 0.1% DMSO serving as the vehicle control.

For Transwell migration and invasion assays, 5 × 104 cells were seeded in the upper chambers of Transwell inserts (8 μm pore size; Labselect, Beijing, China; Cat# 14347) with or without Matrigel coating (for invasion). Medium containing 10% fetal bovine serum was added to the lower chamber as chemoattractant. After 24 hours, migrated or invaded cells were fixed, stained with crystal violet, and counted in five random fields. Cells were treated with PGE2 (1 μM) or celecoxib (10 μM) as indicated, and 0.1% DMSO served as the vehicle control.

ELISA

PGE2 concentrations in culture supernatants were measured using an ELISA kit (Upingbio, Shanghai, China; Cat# YPJ1896) according to the manufacturer’s instructions. Absorbance was recorded at 450 nm and normalized to cell number.

m6A RNA immunoprecipitation-qPCR and RNA immunoprecipitation-qPCR

For m6A RNA immunoprecipitation (MeRIP)-qPCR, total RNA was extracted, fragmented to approximately 100 nt, and incubated with an anti-m6A antibody (BersinBio, Guangzhou, China; Cat# Bes5203) or control IgG for immunoprecipitation. A portion of fragmented RNA was reserved as input. The immunoprecipitated RNA was then purified and subjected to RT-qPCR analysis using COX-2-specific primers.

For RIP-qPCR, cell lysates were prepared and incubated with magnetic beads conjugated with an anti-hnRNPA2B1 antibody or control IgG using the Magna RIP Kit (MilliporeSigma, St. Louis, MO, United States; Cat# 17-700) according to the manufacturer’s instructions. A portion of the lysate was saved as input. The co-immunoprecipitated RNAs were purified and analyzed by RT-qPCR to determine the enrichment of COX-2 mRNA.

Actinomycin D assay

Eca109 cells were treated with actinomycin D (5 μg/mL; Selleck Chemicals, Houston, TX, United States; Cat# S8964) to inhibit transcription. Total RNA was extracted at the indicated time points (0, 2, 4, and 6 hours) using TRIzol reagent (Invitrogen, United States). COX-2 mRNA levels were measured by RT-qPCR and normalized to GAPDH. Relative COX-2 mRNA abundance at each time point was calculated relative to that at 0 hour to assess mRNA decay.

In vivo xenograft assay

Male BALB/c nude mice (6 weeks old) were obtained from Guangdong Yaokang Biotechnology Co., Ltd. (Foshan, Guangdong, China). Six mice were included in each group. The sample size was determined based on previous xenograft studies, experimental feasibility, and ethical considerations to minimize animal use. All animal procedures were designed to minimize pain and discomfort. Mice were housed under specific pathogen-free conditions at 22 ± 2 °C with 50% ± 10% humidity under a 12 hours light/dark cycle, with ad libitum access to food and water, and were acclimatized to the animal facility for 1 week before experimentation. A total of 5 × 106 Eca109 cells stably transduced with shCON, sh-hnRNPA2B1-56-1, NC, or OE were injected subcutaneously. Tumor volumes were measured with calipers (V = 0.5 × L × W2). At the endpoint, mice were euthanized by CO2 inhalation, and tumors were excised, weighed, and processed for hematoxylin and eosin and immunohistochemistry (IHC) staining. Tumor inhibition rate was calculated based on endpoint tumor weight using the following formula: Tumor inhibition rate (%) = (1 - mean tumor weight of treatment group/mean tumor weight of corresponding control group) × 100.

For pharmacological studies, tumor-bearing mice were monitored for 21 days in total. Treatment with celecoxib (20 mg/kg, intraperitoneally every other day) or vehicle control began 7 days after tumor formation and continued for 14 days[22,23]. Tumor size and weight were recorded. All animal experiments were reviewed and approved by the Institutional Animal Care and Use Committee of Fujian Medical University, No. 2022-NSFC-0403 and were performed in accordance with institutional guidelines.

IHC

Paraffin-embedded tissue sections were deparaffinized in xylene and rehydrated through graded ethanol solutions (100%, 95%, 80%, and 70%). Antigen retrieval was performed by heating sections in citrate buffer (pH 6.0; Servicebio, Wuhan, China; Cat# G1202) for 15 minutes. Endogenous peroxidase activity was blocked using 3% hydrogen peroxide (Servicebio, Wuhan, China; Cat# G0115) for 10 minutes, followed by blocking with 5% goat serum (Servicebio, Wuhan, China; Cat# G5001) for 30 minutes at room temperature.

Sections were incubated overnight at 4 °C with the following primary antibodies: Anti-hnRNPA2B1 (1:250; Abcam, Cambridge, United Kingdom; Cat# ab259894), anti-COX-2 (1:250; Cell Signaling Technology, Danvers, MA, United States; Cat# 12282S), and anti-Ki-67 (1:250; Abcam, Cambridge, United Kingdom; Cat# ab16667). After washing with phosphate-buffered saline, sections were incubated with an HRP-conjugated secondary antibody (Servicebio, Wuhan, China; Cat# G1215) for 50 minutes at room temperature. Signal was developed using a DAB substrate kit (ZSGB-Bio, Beijing, China; Cat# ZLI-9018), and nuclei were counterstained with hematoxylin (Servicebio, Wuhan, China; Cat# G1004). Finally, sections were dehydrated, cleared, and mounted with neutral resin (Servicebio, Wuhan, China; Cat# G8590) for microscopic observation. IHC-positive cells were quantified using identical image acquisition and analysis criteria across groups, and the quantification was reviewed by investigators blinded to group information.

Statistical analysis

All in vitro experiments were independently performed at least three times as biological replicates, with technical replicates included where applicable. Data are presented as the mean ± SD. Comparisons between two groups were analyzed using the Student’s t-test. For comparisons among more than two groups, one-way ANOVA followed by appropriate post hoc tests was performed. When experiments involved two independent factors, two-way ANOVA was applied. Categorical variables were evaluated using the χ2 test or Fisher’s precision probability test. Correlations between continuous variables were assessed using Pearson’s correlation analysis. Overall survival was analyzed using the Kaplan-Meier method and compared with the log-rank test. A two-sided P < 0.05 was considered statistically significant. All statistical analyses were conducted using GraphPad Prism 10.0 (GraphPad Software, La Jolla, CA, United States) and SPSS Statistics 27.0 (IBM, Armonk, NY, United States).

RESULTS
HnRNPA2B1 is aberrantly upregulated in ESCC and is associated with poor prognosis

To investigate the clinical significance of hnRNPA2B1 in ESCC, we first analyzed transcriptomic data from TCGA and Gene Expression Omnibus cohorts. HnRNPA2B1 was markedly upregulated in ESCA tissues compared with normal counterparts (Figure 1A and B). Expression levels increased with advancing clinical stage, with higher HnRNPA2B1 expression observed in more advanced tumors (Figure 1C). Kaplan-Meier survival analysis further demonstrated that patients with high hnRNPA2B1 expression had significantly worse overall survival than those with low expression (Figure 1D).

Figure 1
Figure 1 HnRNPA2B1 upregulation is an early event in esophageal squamous cell carcinoma progression and is associated with poor prognosis. A: Pan-cancer analysis of hnRNPA2B1 expression in normal (left) and tumor (right) tissues from The Cancer Genome Atlas, with the esophageal cohort highlighted; B: HnRNPA2B1 expression in esophageal carcinoma (ESCA) and normal esophageal tissues in four Gene Expression Omnibus datasets (GSE20347, GSE77861, GSE161533, and GSE111011); C: HnRNPA2B1 expression in normal tissues and ESCA at different pathological stages in The Cancer Genome Atlas; D: Kaplan-Meier overall survival analysis of ESCA patients with high or low HnRNPA2B1 expression; E: Representative immunohistochemistry staining of hnRNPA2B1 in normal esophageal epithelium, LGIN, and HGIN from ESD specimens; F: Quantification of hnRNPA2B1-positive cells in normal, LGIN, and HGIN tissues; G: Representative immunohistochemistry staining of hnRNPA2B1 in esophageal squamous cell carcinoma (ESCC) and matched adjacent and distal tissues; H: Quantification of hnRNPA2B1-positive cells in ESCC, adjacent, and distal tissues; I: Western blotting analysis of hnRNPA2B1 in paired adjacent (A) and tumor (T) ESCC tissues from three representative cases, with GAPDH as the loading control; J: Relative hnRNPA2B1 mRNA expression in paired adjacent and ESCC tissues determined by reverse transcription-quantitative polymerase chain reaction; K: Western blotting analysis of hnRNPA2B1 in Het-1A, Eca109, and KYSE150 cells, with β-actin as the loading control; L: Relative hnRNPA2B1 mRNA expression in Het-1A, Eca109, and KYSE150 cells determined by reverse transcription-quantitative polymerase chain reaction. Data are presented as mean ± SD. Statistical significance was determined by Student’s t-test for two-group comparisons or one-way ANOVA for multiple-group comparisons, as appropriate. Scale bars = 100 μm. aP < 0.05, bP < 0.01, and cP < 0.001. ESCA: Esophageal carcinoma; HGIN: High-grade intraepithelial neoplasia; LGIN: Low-grade intraepithelial neoplasia; ESCC: Esophageal squamous cell carcinoma.

To assess hnRNPA2B1 expression in early lesions, we performed IHC staining on samples from ESD. HnRNPA2B1 expression was elevated in LGIN and further increased in HGIN (Figure 1E and F). These findings are consistent with recent single-cell spatial atlases revealing early molecular alterations during esophageal carcinogenesis[24]. In line with these observations, tissue microarray analysis revealed stronger hnRNPA2B1 staining in ESCC compared with adjacent and distant normal tissues (Figure 1G and H). Elevated hnRNPA2B1 expression correlated positively with pathological grade (Table 1). In matched ESCC patient samples, both hnRNPA2B1 protein and mRNA levels were significantly higher in tumors than in adjacent normal tissues (Figure 1I and J). Moreover, elevated hnRNPA2B1 expression was also detected in ESCC cell lines (Eca109, KYSE150) compared with the normal esophageal epithelial cell line Het-1A (Figure 1K and L). Collectively, these results demonstrate that hnRNPA2B1 is aberrantly upregulated in ESCC, detectable even at the stage of LGIN, and positively associated with tumor stage, grade, and poor prognosis, suggesting a role in early malignant transformation.

Table 1 The correlation between expression of hnRNPA2B1 and clinicopathologic characteristics of 30 esophageal squamous cell carcinoma patients.
VariablesHnRNPA2B1 expression
P1 value
High
Low
Age (years)
< 601171.000
≥ 6084
Gender
Female431.000
Male158
Grade
I690.021
II-III132
Tumor size (cm)
< 51660.104
≥ 535
Lymph node metastasis
N01150.412
N+44
COX-2 expression
High1850.004
Low16
HnRNPA2B1 promotes ESCC cell proliferation, migration, and invasion in vitro

To investigate the functional role of hnRNPA2B1, we established stable knockdown in Eca109 and KYSE150 cells using three independent short hairpin RNA plasmids (sh56-1, sh57-12, and sh58-1) (Supplementary Figure 2A and B). Knockdown efficiency was confirmed by green fluorescent protein fluorescence imaging, western blotting, and RT-qPCR (Supplementary Figure 2C; Figure 2A and B). Among them, sh56-1 achieved the most robust suppression and was selected for subsequent experiments. In parallel, stable hnRNPA2B1 overexpression was generated in Eca109, KYSE150, and the non-malignant esophageal epithelial cell line Het-1A (Supplementary Figure 2D), with overexpression validated by green fluorescent protein fluorescence, protein and mRNA analysis (Supplementary Figure 2E; Figure 2C and D). These models provided a reliable system for functional characterization of hnRNPA2B1.

Figure 2
Figure 2 HnRNPA2B1 promotes esophageal squamous cell carcinoma cell proliferation, migration, and invasion in vitro. A: Western blotting analysis of hnRNPA2B1 knockdown efficiency in Eca109 and KYSE150 cells transduced with control short hairpin RNA or three independent hnRNPA2B1-targeting short hairpin RNAs (sh56-1, sh57-12, and sh58-1), with β-actin as the loading control; B: Relative hnRNPA2B1 mRNA expression in Eca109 and KYSE150 cells after hnRNPA2B1 knockdown with sh56-1, as determined by reverse transcription-quantitative polymerase chain reaction; C: Western blotting analysis of hnRNPA2B1 overexpression in Eca109, KYSE150, and Het-1A cells transduced with negative control or hnRNPA2B1 vector (overexpression), with β-actin as the loading control; D: Relative hnRNPA2B1 mRNA expression in Eca109, KYSE150, and Het-1A cells after hnRNPA2B1 overexpression, as determined by reverse transcription-quantitative polymerase chain reaction; E: Cell counting kit-8 assays showing the effects of hnRNPA2B1 knockdown or overexpression on the proliferation of Eca109, KYSE150, and Het-1A cells; F: Representative images of colony formation assays in Eca109 cells after hnRNPA2B1 knockdown or overexpression; G: Quantification of colony numbers in (F); H: Representative images of Transwell migration and invasion assays in Eca109 cells after hnRNPA2B1 knockdown or overexpression; I: Quantification of migrated and invaded cells in (H); J: Representative images of wound-healing assays in Eca109 cells at 0, 12, and 24 hours after hnRNPA2B1 knockdown or overexpression; K: Quantification of wound closure in (J). Data are presented as mean ± SD. Statistical significance was determined by Student’s t-test for two-group comparisons or one-way ANOVA for multiple-group comparisons, as appropriate. aP < 0.05, bP < 0.01, and cP < 0.001. NC: Negative control; OE: Overexpression; shCON: Control short hairpin RNA; shA2B1: HnRNPA2B1 short hairpin RNA.

Using these models, we next evaluated the impact of hnRNPA2B1 on cell growth. Cell counting kit-8 assays showed that hnRNPA2B1 knockdown markedly inhibited proliferation of Eca109 and KYSE150 cells. By contrast, hnRNPA2B1 overexpression enhanced proliferation in both ESCC cells and Het-1A cells (Figure 2E). Consistently, colony formation assays in Eca109 cells showed that hnRNPA2B1 knockdown reduced clonogenic capacity, whereas hnRNPA2B1 overexpression promoted colony formation (Figure 2F and G).

We further assessed the effects of hnRNPA2B1 on cell motility. Transwell assays demonstrated that hnRNPA2B1 silencing significantly suppressed both migration and invasion of Eca109 cells, whereas hnRNPA2B1 overexpression markedly promoted these phenotypes (Figure 2H and I). Wound-healing assays confirmed that hnRNPA2B1 knockdown impaired, while overexpression accelerated migratory capacity (Figure 2J and K). These findings demonstrate that hnRNPA2B1 enhances ESCC cell proliferation, migration, and invasion in vitro.

RNA-seq analyses identify PTGS2/COX-2 as a downstream effector of hnRNPA2B1

To identify potential downstream effectors of hnRNPA2B1, we performed bulk RNA-seq in Eca109 cells with or without hnRNPA2B1 knockdown. PTGS2, which encodes COX-2, was among the significantly downregulated genes after hnRNPA2B1 silencing (Figure 3A and B). Gene Ontology enrichment analysis revealed that hnRNPA2B1-regulated genes were mainly enriched in metabolic and redox-related pathways, including oxidation-reduction processes, organic acid metabolism, and oxidoreductase activity (Figure 3C). These findings suggest that hnRNPA2B1 may regulate multiple biological processes in ESCC.

Figure 3
Figure 3 Transcriptomic analyses identify prostaglandin-endoperoxide synthase 2/cyclooxygenase-2 as a downstream target of hnRNPA2B1 in esophageal carcinoma. A: Volcano plot of differentially expressed genes in Eca109 cells transduced with control short hairpin RNA or hnRNPA2B1 short hairpin RNA, based on bulk RNA-seq; B: Heatmap showing the relative expression of prostaglandin-endoperoxide synthase 2 (PTGS2) and HnRNPA2B1 in control short hairpin RNA and hnRNPA2B1 short hairpin RNA samples; C: Gene Ontology enrichment analysis of HNRNPA2B1-regulated differentially expressed genes, presented as a bubble plot; D: Pan-cancer analysis of PTGS2 expression in normal (left) and tumor (right) tissues from The Cancer Genome Atlas, with the esophageal cohort highlighted; E: Distribution of PTGS2 and HnRNPA2B1 expression in normal, primary tumor, and metastatic tissues based on The Cancer Genome Atlas data; F: Correlation analysis between HnRNPA2B1 and PTGS2 expression in esophageal carcinoma samples. Statistical significance was assessed as indicated in each panel. Pearson’s correlation analysis was used in (F). aP < 0.05. shCON: Control short hairpin RNA; shA2B1: HnRNPA2B1 short hairpin RNA; PTGS2: Prostaglandin-endoperoxide synthase 2.

To further support the association between hnRNPA2B1 and COX-2 in patient datasets, TCGA pan-cancer profiling showed that PTGS2 expression was elevated across multiple tumor types, including ESCA (Figure 3D). Comparative analysis of normal, primary, and metastatic tissues demonstrated increased expression of both hnRNPA2B1 and PTGS2 during tumor progression (Figure 3E). Correlation analysis in ESCA samples further confirmed a positive association between hnRNPA2B1 and PTGS2 expression (Figure 3F).

HnRNPA2B1 promotes COX-2 expression and stabilizes its mRNA in an m6A-related manner

To further validate the association between hnRNPA2B1 and COX-2, we first examined COX-2 expression in ESCC cell lines. Western blotting and RT-qPCR revealed that COX-2 levels were markedly higher in tumor tissues and ESCC cell lines compared with normal controls (Figure 4A-D), consistent with previous findings that COX-2 is significantly overexpressed in ESCC or hepatocellular carcinoma compared with adjacent non-tumorous tissues[25,26]. Silencing hnRNPA2B1 decreased COX-2 expression. In contrast, hnRNPA2B1 overexpression led to increased levels of COX-2 (Figure 4E-H). Similar to hnRNPA2B1 expression, immunohistochemical analysis of ESD specimens revealed that COX-2 expression was elevated in LGIN and further upregulated in HGIN (Figure 4I and J). Consistently, COX-2 staining was markedly stronger in ESCC tissues compared with adjacent or distant normal tissue (Figure 4K and L).

Figure 4
Figure 4 HnRNPA2B1 stabilizes m6A-modified cyclooxygenase-2 mRNA and promotes its expression in esophageal squamous cell carcinoma. A: Western blotting analysis of cyclooxygenase-2 (COX-2) expression in paired adjacent (A) and tumor (T) esophageal squamous cell carcinoma (ESCC) tissues from three representative patients; B: Relative COX-2 mRNA expression in paired adjacent and ESCC tissues determined by reverse transcription-quantitative polymerase chain reaction (RT-qPCR); C: Western blotting analysis of COX-2 expression in Het-1A cells and ESCC cell lines (Eca109 and KYSE150); D: Relative COX-2 mRNA expression in Het-1A, Eca109, and KYSE150 cells determined by RT-qPCR; E: Western blotting analysis of COX-2 expression in Eca109 and KYSE150 cells transduced with control short hairpin RNA or hnRNPA2B1 short hairpin RNA (sh56-1); F: Relative COX-2 mRNA expression in Eca109 and KYSE150 cells after hnRNPA2B1 knockdown, as determined by RT-qPCR; G: Western blotting analysis of COX-2 expression in Eca109 and KYSE150 cells transduced with negative control or hnRNPA2B1 overexpression vector; H: Relative COX-2 mRNA expression in Eca109 and KYSE150 cells after hnRNPA2B1 overexpression, as determined by RT-qPCR; I: Representative immunohistochemistry staining of COX-2 in normal esophageal epithelium, low-grade intraepithelial neoplasia, and high-grade intraepithelial neoplasia from endoscopic submucosal dissection specimens; J: Quantification of COX-2-positive cells in normal, low-grade intraepithelial neoplasia, and high-grade intraepithelial neoplasia tissues; K: Representative immunohistochemistry staining of COX-2 in ESCC tissues and matched adjacent and distal tissues; L: Quantification of COX-2-positive cells in ESCC, adjacent, and distal tissues; M: Methylated RNA immunoprecipitation-qPCR analysis in Eca109 cells showing enrichment of COX-2 mRNA in the m6A-immunoprecipitated fraction relative to the IgG control; N: RNA immunoprecipitation-qPCR analysis in Eca109 cells showing enrichment of COX-2 mRNA in hnRNPA2B1 immunoprecipitates, with reduced enrichment after hnRNPA2B1 knockdown and increased enrichment after hnRNPA2B1 overexpression; O: Actinomycin D assay in Eca109 cells showing COX-2 mRNA decay after hnRNPA2B1 knockdown or overexpression. (Panels A, C, E, and G were derived from the same membranes as Figure 1I and K, Figure 2A and C, respectively, so that hnRNPA2B1 and COX-2 expression could be assessed in the same samples. The corresponding loading controls are shown in Figure 1I and K, Figure 2A and C, respectively.) Data are presented as mean ± SD. Statistical significance was determined by Student’s t-test for two-group comparisons or one-way ANOVA for multiple-group comparisons, as appropriate. Scale bars = 100 μm. aP < 0.05, bP < 0.01, and cP < 0.001. COX-2: Cyclooxygenase-2; ESCC: Esophageal squamous cell carcinoma; shCON: Control short hairpin RNA; NC: Negative control; OE: Overexpression; shA2B1: HnRNPA2B1 short hairpin RNA; LGIN: Low-grade intraepithelial neoplasia; HGIN: High-grade intraepithelial neoplasia; MeRIP: Methylated RNA immunoprecipitation; RIP: RNA immunoprecipitation.

We then performed MeRIP-qPCR to determine whether COX-2 mRNA harbors m6A modification, and observed significant enrichment in the m6A immunoprecipitated fraction (Figure 4M). RIP-qPCR further confirmed the association between hnRNPA2B1 and COX-2 transcripts, which was reduced by hnRNPA2B1 knockdown and enhanced by hnRNPA2B1 overexpression (Figure 4N). Actinomycin D assays showed that hnRNPA2B1 knockdown accelerated, whereas hnRNPA2B1 overexpression delayed, COX-2 mRNA decay (Figure 4O). Collectively, these data indicate that hnRNPA2B1 promotes COX-2 expression by stabilizing its mRNA through an m6A-related mechanism, and that COX-2 expression follows a pattern similar to that of hnRNPA2B1 during ESCC progression.

Targeting the hnRNPA2B1-COX-2/PGE2 axis suppresses ESCC cell growth and motility in vitro

To determine the functional significance of the hnRNPA2B1-COX-2/PGE2 axis, we first measured PGE2 production. ELISA revealed that PGE2 secretion was reduced after hnRNPA2B1 knockdown and elevated upon hnRNPA2B1 overexpression (Figure 5A). We next investigated whether PGE2 supplementation or COX-2 inhibition modulates hnRNPA2B1-driven phenotypes. Exogenous PGE2 partially rescued the proliferation defects induced by hnRNPA2B1 silencing, while celecoxib treatment markedly suppressed proliferation and colony formation in hnRNPA2B1-overexpressing cells (Figure 5B-D). Consistently, functional assays showed that PGE2 supplementation restored, whereas celecoxib significantly impaired, the migration and invasion capacities of ESCC cells (Figure 5E-I). Together, these results indicate that hnRNPA2B1 promotes ESCC cell growth and motility through the COX-2/PGE2 axis, and that pharmacological inhibition of COX-2 effectively counteracts hnRNPA2B1-driven malignant phenotypes in vitro.

Figure 5
Figure 5 Targeting the hnRNPA2B1-cyclooxygenase-2/prostaglandin E2 axis suppresses esophageal squamous cell carcinoma cell growth and motility in vitro. A: Relative prostaglandin E2 (PGE2) production measured by ELISA in esophageal squamous cell carcinoma cells after hnRNPA2B1 knockdown or overexpression; B: Cell counting kit-8 assays showing the effects of PGE2 supplementation (1 μM) on hnRNPA2B1-silenced cells and celecoxib treatment (10 μM) on hnRNPA2B1-overexpressing cells, with 0.1% DMSO as the vehicle control; C: Representative images of colony formation assays under the indicated conditions; D: Quantification of colony numbers in (C); E: Representative images of Transwell migration and invasion assays under the indicated conditions; F: Quantification of migrated and invaded cells in (E); G: Representative images of wound-healing assays at 0, 12, and 24 hours under the indicated conditions; H: Quantification of wound closure in the PGE2 rescue experiments; I: Quantification of wound closure in the celecoxib inhibition experiments. Data are presented as mean ± SD. Statistical significance was determined by Student’s t-test for two-group comparisons or one-way ANOVA for multiple-group comparisons, as appropriate. aP < 0.05, bP < 0.01, and cP < 0.001. shCON: Control short hairpin RNA; NC: Negative control; OE: Overexpression; shA2B1: HnRNPA2B1 short hairpin RNA; PGE2: Prostaglandin E2; DMSO: Dimethyl sulfoxide.
Targeting the hnRNPA2B1-COX-2/PGE2 axis suppresses ESCC tumor growth in vivo

To assess whether the in vitro observations translated into tumorigenic potential in vivo, we established xenograft models using Eca109 cells with stable hnRNPA2B1 knockdown or overexpression. Tumors derived from hnRNPA2B1-knockdown cells exhibited markedly slower growth and reduced final volume and weight, whereas hnRNPA2B1 overexpression significantly accelerated tumor growth (Figure 6A and B). Based on endpoint tumor weight, the tumor inhibition rate of hnRNPA2B1 knockdown was 43.94% compared with the shCON group. Histological examination confirmed tumor formation, and IHC demonstrated decreased hnRNPA2B1, COX-2 and Ki-67 staining in knockdown tumors, with the opposite pattern observed in hnRNPA2B1-overexpressing xenografts (Figure 6C and D). These results indicate that hnRNPA2B1 promotes ESCC tumor growth in vivo.

Figure 6
Figure 6 Targeting the hnRNPA2B1-cyclooxygenase-2/prostaglandin E2 axis suppresses esophageal squamous cell carcinoma tumor growth in vivo. A: Representative images of subcutaneous xenograft tumors derived from Eca109 cells with hnRNPA2B1 knockdown (hnRNPA2B1 short hairpin RNA) or corresponding control (control short hairpin RNA), and hnRNPA2B1 overexpression or corresponding control (negative control); B: Quantification of tumor volume and tumor weight in the indicated groups; C: Representative hematoxylin and eosin and immunohistochemistry staining of hnRNPA2B1, cyclooxygenase-2 (COX-2), and Ki-67 in xenograft tumors from the indicated groups; D: Quantification of hnRNPA2B1-, COX-2-, and Ki-67-positive cells in (C); E: Schematic diagram of the xenograft study design and celecoxib treatment regimen; F: Representative images of xenograft tumors from mice bearing negative control or overexpression Eca109 xenografts with or without celecoxib treatment, and parental Eca109 xenografts with or without celecoxib treatment; G: Quantification of tumor volume and tumor weight in the groups shown in (F); H: Representative hematoxylin and eosin and immunohistochemistry staining of hnRNPA2B1, COX-2, and Ki-67 in xenograft tumors with or without celecoxib treatment, as indicated; I: Quantification of hnRNPA2B1-, COX-2-, and Ki-67-positive cells in (H). Data are presented as mean ± SD. Statistical significance was determined by Student’s t-test for two-group comparisons or one-way ANOVA for multiple-group comparisons, as appropriate. Scale bars = 1 cm in (A) and (F), and 100 μm in (C) and (H). aP < 0.05, bP < 0.01, and cP < 0.001. shCON: Control short hairpin RNA; NC: Negative control; OE: Overexpression; shA2B1: HnRNPA2B1 short hairpin RNA; COX-2: Cyclooxygenase-2; H&E: Hematoxylin and eosin; i.p.: Intraperitoneal; s.c.: Subcutaneous.

To further evaluate the therapeutic relevance of the hnRNPA2B1-COX-2/PGE2 axis, we administered celecoxib in xenograft-bearing mice (Figure 6E). Celecoxib treatment significantly suppressed tumor growth, as evidenced by reduced tumor volume and weight (Figure 6F and G). Based on endpoint tumor weight, the tumor inhibition rates were 67.31% for celecoxib treatment in hnRNPA2B1-overexpressing xenografts compared with the OE + vehicle group and 60.50% for celecoxib treatment in Eca109 xenografts compared with the Eca109 + vehicle group. IHC analysis further revealed decreased COX-2 and Ki-67 expression in treated tumors (Figure 6H and I). Together, these findings validate that the oncogenic effects of hnRNPA2B1 in ESCC are mediated, at least in part, by the COX-2/PGE2 axis, and that pharmacological inhibition of COX-2 attenuates hnRNPA2B1-driven tumor growth in vivo.

DISCUSSION

In this study, we identified hnRNPA2B1 as an oncogenic m6A reader that promotes the progression of ESCC through stabilization of COX-2 mRNA and activation of the COX-2/PGE2 axis. HnRNPA2B1 was markedly upregulated in ESCC and early intraepithelial neoplasia, correlating with pathological grade and poor survival. Functional experiments demonstrated that hnRNPA2B1 enhances tumor cell proliferation, migration, invasion, and xenograft growth. Mechanistically, hnRNPA2B1 was associated with m6A -modified COX-2 transcripts, stabilized their mRNA, and augmented COX-2/PGE2 axis (Figure 7). Importantly, exogenous PGE2 supplementation restored the suppressed phenotypes caused by hnRNPA2B1 silencing, whereas celecoxib treatment abrogated hnRNPA2B1-induced tumorigenic effects both in vitro and in vivo.

Figure 7
Figure 7 Schematic model illustrating the hnRNPA2B1-cyclooxygenase-2/prostaglandin E2 axis in esophageal squamous cell carcinoma progression. (Created in BioRender. https://BioRender.com/5wlddw2). LGIN: Low-grade intraepithelial neoplasia; HGIN: High-grade intraepithelial neoplasia; ESCC: Esophageal squamous cell carcinoma; PGE2: Prostaglandin E2; COX-2: Cyclooxygenase-2; m6A: N6-methyladenosine.

Previous studies have demonstrated diverse oncogenic roles of hnRNPA2B1 across multiple tumor types, including lung, breast, and liver cancers[27,28]. In ESCC, Guo et al[29] reported that hnRNPA2B1 upregulates ACLY and ACC1 to enhance de novo fatty acid synthesis, while Li et al[16] demonstrated that hnRNPA2B1 modulates the oncogenic miR-17-92 cluster to influence tumor progression. HnRNPA2B1 has been reported to dynamically shuttle between the nucleus and cytoplasm, where it participates in mRNA export and stability[27-30]. This dynamic localization may partly explain its potential dual function in transcriptional and post-transcriptional regulation of COX-2, as suggested in previous studies. Unlike these metabolic or microRNA-mediated mechanisms, our findings identify a distinct m6A-related post-transcriptional mechanism by which hnRNPA2B1 stabilizes COX-2 mRNA, thereby contributing to inflammation-associated tumor progression.

COX-2, a rate-limiting enzyme in prostaglandin biosynthesis[17], is frequently overexpressed in ESCC and contributes to proliferation, invasion, angiogenesis, and immune evasion[31,32]. COX-2 regulation by hnRNPA2B1 appears to be multilayered. Xuan et al[33] reported transcriptional activation of COX-2 via hnRNPA2B1-p300 interaction in lung cancer, whereas our work demonstrates post-transcriptional mRNA stabilization in ESCC. Beyond RNA methylation, broader epigenetic alterations also contribute to ESCC initiation; for example, CBX family upregulation correlates with TP53 mutations and immune infiltration[34]. Consistent with this epigenetic landscape, another study recently established a genetically engineered esophageal organoid model showing that Trp53 knockout or KrasG12D expression is sufficient to induce squamous neoplastic transformation, recapitulating early transcriptional and chromatin alterations characteristic of ESCC initiation[35]. Together, these findings suggest that early genetic alterations, particularly TP53 loss, create an epigenetically permissive environment that cooperates with inflammation-associated pathways such as hnRNPA2B1-mediated COX-2 activation to drive ESCC development. The identification of an hnRNPA2B1-COX-2/PGE2 axis carries important translational implications. COX-2/PGE2 axis is a well-established driver of tumorigenesis[36,37], and our data provide a mechanistic rationale for therapeutic targeting of this pathway in ESCC. In this study, PGE2 was selected as the major functional downstream readout of COX-2 activity, whereas other COX-2-derived prostanoids were not examined and may be worthy of further investigation in future studies. Although our study did not directly assess immune modulation, the known role of COX-2/PGE2 axis in recruiting myeloid-derived suppressor cells and regulatory T cells and suppressing cytotoxic T-cell responses[38,39] suggests that hnRNPA2B1-driven activation of this pathway may also contribute to an immunosuppressive tumor microenvironment in ESCC. While the clinical benefit of celecoxib monotherapy has been limited[40], this may partly reflect the lack of biomarker-based patient selection. Given that hnRNPA2B1 upregulation was associated with enhanced COX-2 expression and activation of the COX-2/PGE2 axis in our study, ESCC patients with high hnRNPA2B1 expression may represent a subgroup more likely to benefit from COX-2 inhibition. Future clinical studies incorporating hnRNPA2B1 expression and COX-2/PGE2 axis activity may help identify patients who could respond better to celecoxib or other COX-2-targeted strategies. Although celecoxib has been extensively characterized and widely used as a selective COX-2 inhibitor, potential COX-2-independent effects cannot be completely excluded. Future studies using genetic PTGS2/COX-2 inhibition, such as short hairpin RNA or siRNA knockdown, together with rescue experiments, would further validate the specificity of the hnRNPA2B1-COX-2/PGE2 axis.

Despite these advances, several limitations should be acknowledged. Although our findings were validated in cell lines, tissue microarrays, and xenograft models, larger clinical cohorts with complete follow-up and clinicopathological information are needed to confirm the prognostic value of hnRNPA2B1 and its correlation with COX-2 in ESCC. In addition, our survival analysis was based primarily on public TCGA-ESCA data; although this dataset provides useful exploratory evidence, the number of available cases and clinicopathological variables suitable for comprehensive Cox regression analysis were limited. Moreover, while COX-2 was identified as a key downstream effector, hnRNPA2B1 may also regulate additional oncogenic transcripts and pathways, including metabolic and redox-related processes suggested by our RNA-seq analysis, which warrant further exploration[41-43]. Although our in vivo experiments using subcutaneous xenograft models support a role of hnRNPA2B1 in promoting ESCC tumor growth, orthotopic and metastatic ESCC models were not included in the present study. Therefore, the effects of the hnRNPA2B1-COX-2/PGE2 axis on local invasion, tumor-microenvironment interactions, and metastatic dissemination require further investigation. In addition, although xenograft tumor tissues were analyzed for hnRNPA2B1, COX-2, and Ki-67 expression by IHC, serum and intratumoral PGE2 levels and additional pathway activity markers were not directly examined and may be evaluated in future studies. Although our data support an m6A-related post-transcriptional mechanism, the precise m6A sites on COX-2 mRNA and their binding affinity with hnRNPA2B1 remain to be mapped using MeRIP-seq, CLIP-seq, or site-directed mutagenesis.

CONCLUSION

In conclusion, our study identifies hnRNPA2B1 as an oncogenic m6A reader that stabilizes COX-2 mRNA and promotes ESCC progression through the COX-2/PGE2 axis. These findings provide new mechanistic insight into COX-2 dysregulation and highlight hnRNPA2B1 as a potential therapeutic target and biomarker for precision treatment of ESCC.

References
1.  Yang H, Wang F, Hallemeier CL, Lerut T, Fu J. Oesophageal cancer. Lancet. 2024;404:1991-2005.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 188]  [Cited by in RCA: 220]  [Article Influence: 110.0]  [Reference Citation Analysis (0)]
2.  Morgan E, Soerjomataram I, Rumgay H, Coleman HG, Thrift AP, Vignat J, Laversanne M, Ferlay J, Arnold M. The Global Landscape of Esophageal Squamous Cell Carcinoma and Esophageal Adenocarcinoma Incidence and Mortality in 2020 and Projections to 2040: New Estimates From GLOBOCAN 2020. Gastroenterology. 2022;163:649-658.e2.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 968]  [Cited by in RCA: 896]  [Article Influence: 224.0]  [Reference Citation Analysis (9)]
3.  Meyer KD, Jaffrey SR. Rethinking m(6)A Readers, Writers, and Erasers. Annu Rev Cell Dev Biol. 2017;33:319-342.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 689]  [Cited by in RCA: 936]  [Article Influence: 104.0]  [Reference Citation Analysis (4)]
4.  Patil DP, Pickering BF, Jaffrey SR. Reading m(6)A in the Transcriptome: m(6)A-Binding Proteins. Trends Cell Biol. 2018;28:113-127.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 300]  [Cited by in RCA: 524]  [Article Influence: 58.2]  [Reference Citation Analysis (3)]
5.  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.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 359]  [Cited by in RCA: 1078]  [Article Influence: 154.0]  [Reference Citation Analysis (6)]
6.  Zaccara S, Ries RJ, Jaffrey SR. Reading, writing and erasing mRNA methylation. Nat Rev Mol Cell Biol. 2019;20:608-624.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 2141]  [Cited by in RCA: 1950]  [Article Influence: 278.6]  [Reference Citation Analysis (17)]
7.  Liu J, Yue Y, Han D, Wang X, Fu Y, Zhang L, Jia G, Yu M, Lu Z, Deng X, Dai Q, Chen W, He C. A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylation. Nat Chem Biol. 2014;10:93-95.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 2974]  [Cited by in RCA: 2816]  [Article Influence: 234.7]  [Reference Citation Analysis (23)]
8.  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-1308.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 1336]  [Cited by in RCA: 1302]  [Article Influence: 118.4]  [Reference Citation Analysis (4)]
9.  Wang J, Zhang J, Liu H, Meng L, Gao X, Zhao Y, Wang C, Gao X, Fan A, Cao T, Fan D, Zhao X, Lu Y. N6-methyladenosine reader hnRNPA2B1 recognizes and stabilizes NEAT1 to confer chemoresistance in gastric cancer. Cancer Commun (Lond). 2024;44:469-490.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 62]  [Cited by in RCA: 62]  [Article Influence: 31.0]  [Reference Citation Analysis (0)]
10.  Tang J, Chen Z, Wang Q, Hao W, Gao WQ, Xu H. hnRNPA2B1 Promotes Colon Cancer Progression via the MAPK Pathway. Front Genet. 2021;12:666451.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 6]  [Cited by in RCA: 39]  [Article Influence: 7.8]  [Reference Citation Analysis (0)]
11.  Hao W, Chen Z, Tang J, Yang R, Gao WQ, Xu H. hnRNPA2B1 promotes the occurrence and progression of hepatocellular carcinoma by downregulating PCK1 mRNA via a m6A RNA methylation manner. J Transl Med. 2023;21:861.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in RCA: 21]  [Reference Citation Analysis (0)]
12.  Zhu F, Yang T, Yao M, Shen T, Fang C. HNRNPA2B1, as a m(6)A Reader, Promotes Tumorigenesis and Metastasis of Oral Squamous Cell Carcinoma. Front Oncol. 2021;11:716921.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 7]  [Cited by in RCA: 34]  [Article Influence: 6.8]  [Reference Citation Analysis (1)]
13.  Huang T, Zhu G, Chen F. The Potential Impact of HNRNPA2B1 on Human Cancers Prognosis and Immune Microenvironment. J Immunol Res. 2024;2024:5515307.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in RCA: 6]  [Reference Citation Analysis (0)]
14.  Chen C, Huang L, Sun Q, Yu Z, Wang X, Bu L. HNRNPA2B1 Demonstrates Diagnostic and Prognostic Values Based on Pan-Cancer Analyses. Comput Math Methods Med. 2022;2022:9867660.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in RCA: 5]  [Reference Citation Analysis (0)]
15.  Liu WW, Zheng SQ, Li T, Fei YF, Wang C, Zhang S, Wang F, Jiang GM, Wang H. RNA modifications in cellular metabolism: implications for metabolism-targeted therapy and immunotherapy. Signal Transduct Target Ther. 2024;9:70.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 52]  [Cited by in RCA: 118]  [Article Influence: 59.0]  [Reference Citation Analysis (18)]
16.  Li K, Chen J, Lou X, Li Y, Qian B, Xu D, Wu Y, Ma S, Zhang D, Cui W. HNRNPA2B1 Affects the Prognosis of Esophageal Cancer by Regulating the miR-17-92 Cluster. Front Cell Dev Biol. 2021;9:658642.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 15]  [Cited by in RCA: 27]  [Article Influence: 5.4]  [Reference Citation Analysis (0)]
17.  Smith WL, DeWitt DL, Garavito RM. Cyclooxygenases: structural, cellular, and molecular biology. Annu Rev Biochem. 2000;69:145-182.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 2237]  [Cited by in RCA: 2105]  [Article Influence: 81.0]  [Reference Citation Analysis (3)]
18.  Tanaka MN, Diaz BL, de Souza W, Morgado-Diaz JA. Prostaglandin E2-EP1 and EP2 receptor signaling promotes apical junctional complex disassembly of Caco-2 human colorectal cancer cells. BMC Cell Biol. 2008;9:63.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 24]  [Cited by in RCA: 25]  [Article Influence: 1.4]  [Reference Citation Analysis (0)]
19.  Wong CT, Ahmad E, Li H, Crawford DA. Prostaglandin E2 alters Wnt-dependent migration and proliferation in neuroectodermal stem cells: implications for autism spectrum disorders. Cell Commun Signal. 2014;12:19.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 29]  [Cited by in RCA: 43]  [Article Influence: 3.6]  [Reference Citation Analysis (0)]
20.  Valverde A, Peñarando J, Cañas A, López-Sánchez LM, Conde F, Hernández V, Peralbo E, López-Pedrera C, de la Haba-Rodríguez J, Aranda E, Rodríguez-Ariza A. Simultaneous inhibition of EGFR/VEGFR and cyclooxygenase-2 targets stemness-related pathways in colorectal cancer cells. PLoS One. 2015;10:e0131363.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 25]  [Cited by in RCA: 34]  [Article Influence: 3.1]  [Reference Citation Analysis (0)]
21.  Wang S, Liu Q, Zhang Y, Liu K, Yu P, Liu K, Luan J, Duan H, Lu Z, Wang F, Wu E, Yagasaki K, Zhang G. Suppression of growth, migration and invasion of highly-metastatic human breast cancer cells by berbamine and its molecular mechanisms of action. Mol Cancer. 2009;8:81.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 95]  [Cited by in RCA: 108]  [Article Influence: 6.4]  [Reference Citation Analysis (0)]
22.  Li YX, Wang JL, Gao M, Tang H, Gui R, Fu YF. Celecoxib-erlotinib combination delays growth and inhibits angiogenesis in EGFR-mutated lung cancer. Am J Cancer Res. 2016;6:1494-1510.  [PubMed]  [DOI]
23.  Dilly AK, Honick BD, Lee YJ, Bartlett DL, Choudry HA. Synergistic apoptosis following endoplasmic reticulum stress aggravation in mucinous colon cancer. Orphanet J Rare Dis. 2020;15:211.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 6]  [Cited by in RCA: 13]  [Article Influence: 2.2]  [Reference Citation Analysis (1)]
24.  Chang J, Lu J, Liu Q, Xiang T, Zhang S, Yi Y, Li D, Liu T, Liu Z, Chen X, Dong Z, Li C, Yi H, Yu S, Huang L, Qu F, Wang M, Wang D, Dong H, Cheng G, Zhu L, Li J, Li C, Wu P, Xie X, Teschendorff AE, Lin D, Wang X, Wu C. Single-cell multi-stage spatial evolutional map of esophageal carcinogenesis. Cancer Cell. 2025;43:380-397.e7.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in RCA: 56]  [Reference Citation Analysis (0)]
25.  Xie L, Li R, Zheng B, Xie Z, Fang X, Wang Y, Cuny GD, Li Z, Lin B, Chen X, Hu M. Development of Rofecoxib-Based Fluorescent Probes and Investigations on Their Solvatochromism, AIE Activity, Mechanochromism, and COX-2-Targeted Bioimaging. Anal Chem. 2021;93:11991-12000.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 5]  [Cited by in RCA: 13]  [Article Influence: 2.6]  [Reference Citation Analysis (0)]
26.  Xie L, Li R, Zheng B, Xie Z, Fang X, Dai T, Wang X, Li L, Wang L, Cuny GD, Eriksen J, Tu D, Chen Z, Wang X, Chen X, Hu M. One-Step Transformation from Rofecoxib to a COX-2 NIR Probe for Human Cancer Tissue/Organoid Targeted Bioimaging. ACS Appl Bio Mater. 2021;4:2723-2731.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 6]  [Cited by in RCA: 19]  [Article Influence: 3.8]  [Reference Citation Analysis (0)]
27.  Jin T, Yang L, Chang C, Luo H, Wang R, Gan Y, Sun Y, Guo Y, Tang R, Chen S, Meng D, Dai P, Liu M. HnRNPA2B1 ISGylation Regulates m6A-Tagged mRNA Selective Export via ALYREF/NXF1 Complex to Foster Breast Cancer Development. Adv Sci (Weinh). 2024;11:e2307639.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 27]  [Cited by in RCA: 33]  [Article Influence: 16.5]  [Reference Citation Analysis (1)]
28.  Lu Y, Zou R, Gu Q, Wang X, Zhang J, Ma R, Wang T, Wu J, Feng J, Zhang Y. CRNDE mediated hnRNPA2B1 stability facilitates nuclear export and translation of KRAS in colorectal cancer. Cell Death Dis. 2023;14:611.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in RCA: 14]  [Reference Citation Analysis (0)]
29.  Guo H, Wang B, Xu K, Nie L, Fu Y, Wang Z, Wang Q, Wang S, Zou X. m(6)A Reader HNRNPA2B1 Promotes Esophageal Cancer Progression via Up-Regulation of ACLY and ACC1. Front Oncol. 2020;10:553045.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 84]  [Cited by in RCA: 95]  [Article Influence: 15.8]  [Reference Citation Analysis (4)]
30.  Wang M, Chen L, He J, Xia W, Ye Z, She J. Structural insights into IL-6 signaling inhibition by therapeutic antibodies. Cell Rep. 2024;43:113819.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 2]  [Cited by in RCA: 20]  [Article Influence: 10.0]  [Reference Citation Analysis (0)]
31.  Huang JX, Xiao W, Chen WC, Lin MS, Song ZX, Chen P, Zhang YL, Li FY, Qian RY, Salminen E. Relationship between COX-2 and cell cycle-regulatory proteins in patients with esophageal squamous cell carcinoma. World J Gastroenterol. 2010;16:5975-5981.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in RCA: 10]  [Reference Citation Analysis (0)]
32.  Jin K, Qian C, Lin J, Liu B. Cyclooxygenase-2-Prostaglandin E2 pathway: A key player in tumor-associated immune cells. Front Oncol. 2023;13:1099811.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in RCA: 141]  [Reference Citation Analysis (0)]
33.  Xuan Y, Wang J, Ban L, Lu JJ, Yi C, Li Z, Yu W, Li M, Xu T, Yang W, Tang Z, Tang R, Xiao X, Meng S, Chen Y, Liu Q, Huang W, Guo W, Cui X, Deng W. hnRNPA2/B1 activates cyclooxygenase-2 and promotes tumor growth in human lung cancers. Mol Oncol. 2016;10:610-624.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 36]  [Cited by in RCA: 40]  [Article Influence: 4.0]  [Reference Citation Analysis (0)]
34.  Fang X, Wang J, Chen J, Zhuang M, Huang T, Chen Z, Huang Y, Zheng B, Wang X. Identification and Validation of Chromobox Family Members as Potential Prognostic Biomarkers and Therapeutic Targets for Human Esophageal Cancer. Front Genet. 2022;13:851390.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 2]  [Cited by in RCA: 7]  [Article Influence: 1.8]  [Reference Citation Analysis (0)]
35.  Zheng B, Ko KP, Fang X, Wang X, Zhang J, Jun S, Kim BJ, Luo W, Kim MJ, Jung YS, Cervantes CL, Park JI. A new murine esophageal organoid culture method and organoid-based model of esophageal squamous cell neoplasia. iScience. 2021;24:103440.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 2]  [Cited by in RCA: 29]  [Article Influence: 5.8]  [Reference Citation Analysis (0)]
36.  Sinha P, Clements VK, Fulton AM, Ostrand-Rosenberg S. Prostaglandin E2 promotes tumor progression by inducing myeloid-derived suppressor cells. Cancer Res. 2007;67:4507-4513.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 544]  [Cited by in RCA: 590]  [Article Influence: 31.1]  [Reference Citation Analysis (4)]
37.  Greenhough A, Smartt HJ, Moore AE, Roberts HR, Williams AC, Paraskeva C, Kaidi A. The COX-2/PGE2 pathway: key roles in the hallmarks of cancer and adaptation to the tumour microenvironment. Carcinogenesis. 2009;30:377-386.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 848]  [Cited by in RCA: 945]  [Article Influence: 55.6]  [Reference Citation Analysis (2)]
38.  Wang D, Dubois RN. Eicosanoids and cancer. Nat Rev Cancer. 2010;10:181-193.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 1559]  [Cited by in RCA: 1455]  [Article Influence: 90.9]  [Reference Citation Analysis (5)]
39.  Zelenay S, van der Veen AG, Böttcher JP, Snelgrove KJ, Rogers N, Acton SE, Chakravarty P, Girotti MR, Marais R, Quezada SA, Sahai E, Reis e Sousa C. Cyclooxygenase-Dependent Tumor Growth through Evasion of Immunity. Cell. 2015;162:1257-1270.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 1002]  [Cited by in RCA: 950]  [Article Influence: 86.4]  [Reference Citation Analysis (5)]
40.  Limburg PJ, Wei W, Ahnen DJ, Qiao Y, Hawk ET, Wang G, Giffen CA, Wang G, Roth MJ, Lu N, Korn EL, Ma Y, Caldwell KL, Dong Z, Taylor PR, Dawsey SM. Randomized, placebo-controlled, esophageal squamous cell cancer chemoprevention trial of selenomethionine and celecoxib. Gastroenterology. 2005;129:863-873.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in Crossref: 78]  [Cited by in RCA: 75]  [Article Influence: 3.6]  [Reference Citation Analysis (4)]
41.  Petri BJ, Piell KM, Wilt AE, Howser AD, Winkler L, Whitworth MR, Valdes BL, Lehman NL, Clem BF, Klinge CM. MicroRNA regulation of the serine synthesis pathway in endocrine-resistant breast cancer cells. Endocr Relat Cancer. 2023;30:e230148.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 9]  [Cited by in RCA: 12]  [Article Influence: 4.0]  [Reference Citation Analysis (0)]
42.  Li Y, Li K, Lou X, Wu Y, Seery S, Xu D, Pei Y, Qian B, Wu Y, Liang S, Wu K, Cui W. HNRNPA2B1-Mediated MicroRNA-92a Upregulation and Section Acts as a Promising Noninvasive Diagnostic Biomarker in Colorectal Cancer. Cancers (Basel). 2023;15:1367.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Cited by in RCA: 6]  [Reference Citation Analysis (0)]
43.  Peng WZ, Zhao J, Liu X, Li CF, Si S, Ma R. hnRNPA2B1 regulates the alternative splicing of BIRC5 to promote gastric cancer progression. Cancer Cell Int. 2021;21:281.  [RCA]  [PubMed]  [DOI]  [Full Text]  [Full Text (PDF)]  [Cited by in Crossref: 2]  [Cited by in RCA: 33]  [Article Influence: 6.6]  [Reference Citation Analysis (1)]
Footnotes

Peer review: Externally peer reviewed.

Peer-review model: Single blind

Specialty type: Gastroenterology and hepatology

Country of origin: China

Peer-review report’s classification

Scientific quality: Grade A, Grade A, Grade A

Novelty: Grade A, Grade A, Grade A

Creativity or innovation: Grade A, Grade A, Grade A

Scientific significance: Grade A, Grade A, Grade B

P-Reviewer: Mo S, Associate Professor, PhD, China; Sathish S, Head, Professor, India S-Editor: Wu S L-Editor: A P-Editor: Zhang YL

Write to the Help Desk