Identification of the Expression Patterns and Potential Prognostic Role of m6A-RNA Methylation Regulators in Wilms Tumor

doi:10.21203/rs.3.rs-2582463/v1

Download PDF

Research Article

Identification of the Expression Patterns and Potential Prognostic Role of m6A-RNA Methylation Regulators in Wilms Tumor

https://doi.org/10.21203/rs.3.rs-2582463/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 21 Sep, 2023

Read the published version in BMC Medical Genomics →

You are reading this latest preprint version

Background. To explore the potential role of m6A methylation modification in Wilms Tumor (WT) by m6A-RNA Methylation (m6A) regulators.

Methodology. The association of m6A modification patterns with immune and prognostic characteristics of tumors was systematically evaluated using 19 m6A regulators extracted from Wilms Tumor’s samples in public databases. A comprehensive model of "m6Ascore" was constructed using principal component analysis, and its prognostic value was evaluated.

Results. Almost all m6A regulators were differentially expressed between WT and normal tissues. Unsupervised clustering identified three distinct m6A clusters that differed in both immune cell infiltration and biological pathways. The m6Ascore was constructed to quantify m6A modifications in individual patients. Our analysis suggests that m6Ascore is an independent prognostic factor for WT and can be used as a novel predictor of WT prognosis.

Conclusions.This study comprehensively explored and systematically characterized m6A modifications in WT. m6A modification patterns play a critical role in the tumor immune microenvironment (TIME) and WT prognosis. m6Ascore provides a more comprehensive understanding of m6A modifications in WT and offers a practical tool for predicting WT prognosis. This study will help clinicians to identify valid indicators of WT to improve the poor prognosis of this disease.

m6A-RNA Methylation

Wilms Tumor

m6Ascore

Prognosis

Tumor immune microenvironment

Renal tumors are the fifth most common tumor in children, and Wilms Tumor (WT) is the most common renal malignancy in children¹. The long-term survival rate of WT has steadily increased over the past decades to 85% and has even reached 90% in developed countries^2–4. However, nearly 15% of patients still experience recurrence and associated complications^5–7. Genes such as WT1, WT2 and MYCN family have been suggested to be involved in the development of WT, and the detection of CTR9, DICER1, REST, TP53, TRIM28 and WT1 have been recommended as effective predictors of WT⁸. In addition, SPRY1, SPIN4, MAP7D3, C10orf71, and SPAG have also been analyzed for WT⁹. Recent studies suggest that cholesterol markers and methylation modifications may potentially influence WT^10,11. However, none of these studies could be included in the individualized assessment of specific patients. Newly discovered markers in WT are increasing, and these may eventually play a role in targeted therapies; unfortunately, only one genetic biomarker, LOH at chromosome 1p/16q, has been used in clinical therapy¹². Thus, there is an urgent need to find effective biomarkers to predict the prognosis of WT and to develop new targets for WT therapy.

RNA modifications are universal post-transcriptional modifications that play a crucial role in biological regulation^13,14. N6-methyladenosine (m6A) is the most abundant epigenetic-transcriptomic modification in eukaryotic mRNAs¹⁵. As a reversible RNA modification, m6A plays an important role in the regulation of biological processes such as RNA degradation and splicing¹⁶. There is increasing evidence that m6A modifications play an important role in tumorigenesis and tumor regulation. METTL14, a regulator of m6A, has been reported to affect WT progression and prognosis by regulating related gene expression and splicing patterns¹⁷, and three different multicenter case-control studies have suggested that m6A regulators ALKBH5, YTHDF1 and YTHDF2 all affect WT progression and prognosis to some extent^18–20. Recently, m6A alterations have been shown to be present and to influence tumorigenesis and prognosis in a variety of tumors, including colon, lung, pancreatic, cervical, ovarian, nasopharyngeal, and prostate cancers^21–27.

In the present study, we systematically evaluated m6A modification patterns and the tumor immune microenvironment (TIME) in WT patients. We identified three distinct m6A modification patterns in WT; these clusters differed significantly in prognosis, immune cell infiltration and biological pathways. Based on m6A regulators and related genes, we constructed a model (called "m6Ascore") to quantify the m6A modification patterns in individual patients. This study also suggests that m6Ascore may be a new practical tool to predict the prognosis of WT.

Data acquisition and processing

The RNA sequencing data of Wilms Tumor were obtained from the Cancer Genome Atlas database (TCGA, https://portal.gdc.cancer.gov/). Gene expression data (measured in fragments per kilobase of exon per million fragments mapped or FPKM) were converted to transcripts per kilobase million (TPM).

Analysis of m6A regulators in Wilms Tumor

Based on previous studies, we obtained 19 m6A regulators²⁸. These regulators included 7 “Writer” (METTL3, METTL14, WTAP, VIRMA, METTL16, RBM15, and RBM15B), 9 “Reader” (YTHDF1,YTHDF2,YTHDF3,YTHDC1,YTHDC2,IGF2BP1,IGF2BP2,IGF2BP3 and HNRNPA2B1), and 3 “Eraser” (FTO,ALKBH5 and ALKBH1). The expression profiles of these regulators were systematically extracted and analyzed in normal and tumor samples. Somatic mutations of WT were assessed using the R package "maftools". Tumor mutational burden (TMB) was calculated and the correlation of TMB with clinical features was evaluated. The prognostic value of m6A regulators was assessed using Kaplan-Meier (KM) curves and log-rank tests.

Clustering analysis based on m6A regulators

Based on the expression of m6A regulators, unsupervised clustering was performed using the "ConsensusClusterPlus" R package root to identify different m6A modification patterns in WT patients, and the stability of the clusters was ensured by 1,000 replicates. Survival analysis of different clusters was performed using the KM method. Differences in biological processes between clusters were investigated by gene set variation analysis (GSVA) using the R package "GSVA". The gene set "c2.cp.kegg. v7.4. symbol" was obtained from the Molecular Signatures Database (MSigDB). An adjusted p-value < 0.05 was considered statistically significant.

Comparison of tumor immune microenvironment between different m6A clusters

Single sample gene set enrichment analysis (ssGSEA) was used to quantify the relative infiltration levels of 29 immune cell types in WT samples. The proportion of immune stromal components in the tumor microenvironment (TME) was measured using the " estimate" R package. Wilcoxon rank sum test was used to analyze differences in TME between clusters. In addition, the "limma" R package was used to examine differences in the expression of molecules such as targeted immune checkpoints between different functional clusters.

Prognostic differential expression of genes between different m6A clusters

Principal component analysis (PCA) was used to investigate the different m6A modification patterns in WT. Differentially expressed genes (DEGs) between different m6A clusters were extracted using an empirical Bayesian approach. The significance criterion for DEGs was set to a p-value < 0.05. Gene Ontology (GO) biological process analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used to explore the rich functional annotation of DEGs, and univariate Cox regression analysis was performed to identify prognosis-related genes. Significance criteria were set at p-value < 0.05 and abs (logFC) > 1.

Construction of m6Ascore

A PCA-based scoring system (called "m6Ascore") was constructed to quantify the m6A modification pattern of individual WT patients. Principal components 1 and 2 and 3 were selected as signature scores. m6Ascore was defined using a method similar to the Genomic Grade Index (GGI)^29,30: m6Ascore=∑(PC1i+PC2i+PC3i), where i is the expression of overlapping genes with significant prognostic relevance for DEG in the m6A clusters. The samples were divided into high and low m6Ascore groups based on the scores. Correlation analysis was performed to investigate the relationship between m6Ascore and a number of relevant biological pathways, including survival analysis, immune correlation analysis, clinical correlation analysis, TMB and targeted immune checkpoint molecules.

Statistical analysis

All statistical analyses were performed with R software (version 4.1.2). Paired t-tests were used to compare m6A regulators expression levels in WT tissues. Patient survival was dichotomized for continuous variables using optimal cutoff values determined by the R package "survminer". For prognostic analysis, survival curves were constructed using the KM method, and log-rank tests were used to determine the significance of differences. Receiver operating characteristic (ROC) curves (R package "timeROC") and area under the curve (AUC) values were used to assess the prognostic value of the m6Ascore. Univariate and multivariate independent prognostic analyses were performed to evaluate whether the model was an independent prognostic factor for WT. All statistical p-values were two-sided, and p<0.05 was considered statistically significant.

Differential expression of m6A regulators in Wilms Tumor

We obtained expression data from the TCGA database for a total of 132 samples, including 6 normal tissues and 126 tumor tissues. The results showed significant differences in the expression of almost all m6A regulators between tumor and normal tissues, with the majority of m6A regulators being upregulated in WT tissues (P < 0.001), while ALKBH5 and YTHDC1 were downregulated in WT expression levels (Figure 1.A). In addition, our results showed that TMB in WT patients differed with different clinical characteristics, such as younger age group and diffusely anaplastic Wilms Tumor (DAWT) had higher TMB (Figure 1.B-C).

Univariate COX regression analysis and KM showed that m6A modulators were potential prognostic factors for WT patients, such as RBM15, WTAP and YTHDF2 showed high tumorigenicity (Figure 2.A-C). There was also a significant positive correlation between each m6A modulator (Figure 2.D). In conclusion, m6A regulators showed significant heterogeneity and differential expression in WT tissues compared with normal tissues, and m6A regulators may play a critical role in the development and progression of WT.

m6A modification patterns based on m6A regulators

Based on the expression of 19 m6A regulators, unsupervised cluster analysis was performed and three different m6A modification patterns (clusters 1-3) were identified, with 45 cases in cluster1, 25 cases in cluster2 and 62 cases in cluster3 (Figure 3.A). Survival analysis showed some differences in prognosis between the different modification patterns (Figure 3.B). Further analysis showed differences in regulator expression in the 3 different m6A alteration patterns (Figure 3.C). In addition, we grouped the three different modification patterns into two vehicles for GSVA analysis, and the results showed significant differences in biological behavior between the different modification patterns (Figure 3.D-F).

Immunological characteristics of different m6A modification patterns

The infiltration of 29 immune cells was analyzed by ssGSEA in different m6A modification patterns. The results showed that the immune infiltration was significantly different between the different clusters (Figure 4.A). The results obtained by ESTIMATE algorithm showed that the StromalScore of cluster3 was significantly higher than the other two clusters (Figure 4.B-D). In addition, we analyzed the expression levels of immune checkpoints, immune cell markers and WT metabolites in different clusters³¹. Among the immune checkpoints, cluster1 had higher expression of CD274 and PDCD1, while cluster3 had higher expression of LAG3; the star molecule in WT metabolites, GPC3, was expressed in the 3 clusters, and cluster2 was higher than the other two groups (Figure 5).

Generation of the m6Ascore model

PCA analysis showed that there were different m6A modification patterns in WT patients (Figure 6). GO enrichment analysis and KEGG pathway analysis of different m6A clusters using the R package "clusterProfiler" showed that DEGs were enriched in biological processes related to tumorigenesis and tumor progression, such as rRNA metabolic processes and cell cycle (Figure 7.A-D).

Univariate Cox regression analysis was used to understand the prognostic value of each DEG, and a total of 284 DEGs with significant effect on prognosis were screened to construct m6Ascore, and the scores were grouped into high and low scores based on the best critical value. The association between m6A modification pattern and m6Ascore was analyzed using analysis of variance, and the results showed that cluster1 had a lower score of m6Ascore (Figure 8.A). The KM method was used to show that the subgroup with low m6Ascore had a worse prognosis (Figure 8.B). We used univariate and multivariate Cox regression analysis including gender, age, tumor type, m6Ascore and tumor stage to confirm that m6Ascore was an independent prognostic factor for WT (Figure 8.C-D). The ROC curve showed that m6Ascore had some diagnostic value (Figure 8.E). We further selected a WT data set of m6A methylation modifications, GSE167054, for differential analysis, and the results showed that m6A modifications differed between normal and WT tissues (Figure 8.F). This validated the reliability of the m6Ascore model to some extent.

To further understand the potential biological mechanisms of m6Ascore, we analyzed the correlation between m6Ascore and some biological processes. m6Ascore was closely associated with the infiltration of immune cells, such as a significant positive correlation with dendritic cells and T cells, and a significant negative correlation with neutrophils (Figure 9.A). However, there was no significant correlation between TMB and m6Ascore (Figure 9.B). In addition, our results showed that patients with low m6Ascore had a higher percentage of death (Figure 9.C). Our results also showed significant differences in the expression of immune checkpoints, immune cell markers, and tumor metabolic markers in high and low m6Ascore groups, such as CCL2, CD8A, CD68, CTLA4, HAVCR2, and PDCD1LG2 were highly expressed in the high m6Ascore group, while two immune checkpoints, CD276 and PDCD1, were highly expressed in the low m6Ascore group (Figure 10). This suggests that m6Ascore may become a new evaluation index in immune-targeted therapy.

Currently, an increasing number of studies have shown that m6A modifications play an important role in tumorigenesis and progression^32–34. In general, m6A plays an important role in tumorigenesis and progression by regulating mRNA stability, expression and translation³⁵. For example, METTL3 can accelerate the metastasis and radiation resistance of Glioblastoma by enhancing the stability of SOX2. For example, METTL3 can accelerate the metastasis and radiation resistance of Glioblastoma by enhancing the stability of SOX2 leading to malignant events³⁶; decreased m6A methylation attenuates the expression of the AKT negative regulator PHLPP2, while increasing the expression of the AKT positive regulator mTORC2³⁷. The pathogenesis of Wilms Tumor is still not very clear, and recent several multicenter case studies have suggested that m6A modification patterns may contribute to the development and progression of WT^18–20.However, most of these studies were limited to single or intergenic studies and did not systematically analyze the impact of m6A regulator-related modification patterns on tumor progression and prognosis.

Our study demonstrates for the first time that m6A modifications play an important role in tumorigenesis and progression in WT, with potential prognostic value for WT and promise for treatment guidance. To investigate the role of m6A modifications in WT, we systematically analyzed the m6A modification patterns of WT samples extracted from public databases. Three different m6A modification patterns in WT were finally identified by an unsupervised clustering approach³⁸. The different modification patterns showed significant differences in immune cell infiltration, biological pathways, and prognosis. To further quantify m6A modifications in individual patients, we constructed a model (''m6Ascore") and demonstrated that this model is an independent prognostic factor for WT. Furthermore, our results suggest that m6A modifications are different in different patients and that related modulators may be a novel prognostic marker.

We found that almost all m6A regulators were expressed at significantly higher levels in WT tissues than in normal tissues, and a previous study suggested that the pathogenesis of WT may be related to changes in RNA methylation¹⁰. This suggests that altered m6A modifications are associated with the pathogenesis of WT. The three different modification patterns obtained by clustering analysis differed in terms of survival outcomes. m6A regulator expression was low but with poor survival gains in cluster2, and m6A regulator expression was higher and with better survival gains in cluster3. These results suggest that the expression levels of m6A regulators are closely related to tumor progression in WT. Considering the heterogeneity among m6A modifications, the m6Ascore model was constructed to quantify the m6A modification patterns of individual WT patients. Our results indicate that m6Ascore is an independent prognostic factor for WT patients. In addition, we found that m6Ascore was strongly associated with immune cell infiltration. m6Ascore was significantly positively correlated with dendritic cells and T cells and significantly negatively correlated with neutrophils.

A large number of existing studies have shown that the tumor microenvironment can not only influence tumor cell growth and metastasis, but also has great importance in therapy^39–41, in particular, the association between TME and methylation modifications and the influence of TME on tumor progression in WT has been reported^42,43. We quantified the infiltration of 29 different immune cell types in WT samples by ssGSEA analysis. The results showed a large difference in infiltration between different clusters. Among them, cluster2 had much lower immune infiltration than cluster1 and worse survival. Previous studies have also shown that a high degree of immune infiltration plays an anti-tumor role to some extent^44,45.

Immune checkpoint inhibitor (ICI) therapies targeting pathways such as PD-1 and PD-L1 have been widely used in the treatment of tumors⁴⁶.In particular, ICI therapies have revolutionized the way cancer is treated, offering new hope to a wide range of patients⁴⁷.Our study found that the expression levels of these molecules were different in different WT patients. A similar situation was also seen in the m6Ascore model.

The present study still has the following limitations. First, due to various limitations, we could only assess immune cell infiltration based on algorithms and lack of actual experimental data. Second, due to the lack of data, there is a lack of corresponding validation cohorts and an inability to directly examine the actual immunotherapeutic response in the high and low m6Ascore groups. Similarly, we do not have enough clinical cohorts to validate the prognostic value of m6Ascore in WT, and prospective studies of large cohorts are lacking.

In this study, we systematically analyzed the expression characteristics of m6A-related regulators in WT. m6A modification patterns play an important role in the mechanism and prognosis of WT. Our study provides practical tools for predicting the prognosis of WT, and this study may help clinical practitioners to identify valid indicators of Wilms Tumor for the poor prognosis of this disease.

Ethics approval and consent to participate

Not necessary. The data was from the public database.

Consent for publication

Not applicable.

Availability of data and materials

The raw data of this study are derived from the TCGA database (https://portal.gdc.cancer.gov/), the specific sample is TARGET-WT, which are provided on public databases.

Competing interests

The authors declare that they have no competing interests.

Funding

Supported by Grant 2019GSF108120 from the Key Research and Development Plan of Shandong Province, and by Grant ZR2020MH076 from Natural Science Foundation of Shandong Province.

Authors' contributions

CJ and HG contributed to the design of the study, and CJ drafting the paper. XL and MC contributed to acquisition of data and revision of the paper. WM contributed to analysis and interpretation of data, and revising the paper critically. FS contributed to design, acquisition, analysis, interpretation of data and revising the paper critically. All authors contributed to the article and approved the submitted version.

Acknowledgements

We appreciate the free use of TCGA database.

Jemal, A. et al. Cancer Statistics, 2006. CA. Cancer J. Clin. 56, 106–130 (2006).
Irtan, S., Ehrlich, P. F. & Pritchard-Jones, K. Wilms tumor: “State-of-the-art” update, 2016. Semin. Pediatr. Surg. 25, 250–256 (2016).
Aldrink, J. H. et al. Update on Wilms tumor. J. Pediatr. Surg. 54, 390–397 (2019).
Groenendijk, A. et al. Prognostic Factors for Wilms Tumor Recurrence: A Review of the Literature. Cancers 13, 3142 (2021).
Kalapurakal, J. A. et al. Management of Wilms’ tumour: current practice and future goals. Lancet Oncol. 5, 37–46 (2004).
Brok, J., Treger, T. D., Gooskens, S. L., van den Heuvel-Eibrink, M. M. & Pritchard-Jones, K. Biology and treatment of renal tumours in childhood. Eur. J. Cancer 68, 179–195 (2016).
Davidoff, A. M. Wilms Tumor. Adv. Pediatr. 59, 247–267 (2012).
Bhutani, N., Kajal, P. & Sharma, U. Many faces of Wilms Tumor: Recent advances and future directions. Ann. Med. Surg. 64, 102202 (2021).
Lin, X. et al. Identification of a five‐mRNA signature as a novel potential prognostic biomarker in pediatric Wilms tumor. Mol. Genet. Genomic Med. 8, (2020).
Jablonowski, C. M. et al. TERT Expression in Wilms Tumor Is Regulated by Promoter Mutation or Hypermethylation, WT1, and N-MYC. Cancers 14, 1655 (2022).
He, Y. et al. Using Elevated Cholesterol Synthesis as a Prognostic Marker in Wilms’ Tumor: A Bioinformatic Analysis. BioMed Res. Int. 2021, 1–9 (2021).
Theilen, T.-M. et al. Multidisciplinary Treatment Strategies for Wilms Tumor: Recent Advances, Technical Innovations and Future Directions. Front. Pediatr. 10, 852185 (2022).
Roundtree, I. A., Evans, M. E., Pan, T. & He, C. Dynamic RNA Modifications in Gene Expression Regulation. Cell 169, 1187–1200 (2017).
Boccaletto, P. et al. MODOMICS: a database of RNA modification pathways. 2017 update. Nucleic Acids Res. 46, D303–D307 (2018).
He, L. et al. Functions of N6-methyladenosine and its role in cancer. Mol. Cancer 18, 176 (2019).
Fazi, F. & Fatica, A. Interplay Between N6-Methyladenosine (m6A) and Non-coding RNAs in Cell Development and Cancer. Front. Cell Dev. Biol. 7, 116 (2019).
Zhuo, Z. et al. METTL14 gene polymorphisms decrease Wilms tumor susceptibility in Chinese children. BMC Cancer 21, 1294 (2021).
Liu, Y. et al. Impact of YTHDF1 gene polymorphisms on Wilms tumor susceptibility: A five‐center case‐control study. J. Clin. Lab. Anal. 35, (2021).
Hua, R. et al. ALKBH5 gene polymorphisms and Wilms tumor risk in Chinese children: A five‐center case‐control study. J. Clin. Lab. Anal. 34, (2020).
Wang, Z. et al. The contribution of YTHDF2 gene rs3738067 A>G to the Wilms tumor susceptibility. J. Cancer 12, 6165–6169 (2021).
Liu, Z. et al. Characterization of the m6A-Associated Tumor Immune Microenvironment in Prostate Cancer to Aid Immunotherapy. Front. Immunol. 12, 735170 (2021).
Liu, Z. et al. m6A Regulators Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration Characterization In Nasopharyngeal Carcinoma. Front. Immunol. 12, 762243 (2022).
Chang, L.-L. et al. Emerging role of m6A methylation modification in ovarian cancer. Cancer Cell Int. 21, 663 (2021).
Zhang, W. et al. m6A Regulator-Mediated Tumour Infiltration and Methylation Modification in Cervical Cancer Microenvironment. Front. Immunol. 13, 888650 (2022).
Yang, C., Yu, T., Li, Q., Xie, F. & Lin, Q. A signature based on m6A pattern and tumor microenvironment.
Sun, M. et al. m6A Methylation Modification Patterns and Tumor Microenvironment Infiltration Characterization in Pancreatic Cancer. Front. Immunol. 12, 739768 (2021).
Fan, Y. et al. m6A Regulator-Mediated Methylation Modification Patterns and Characterisation of Tumour Microenvironment Infiltration in Non-Small Cell Lung Cancer. J. Inflamm. Res. Volume 15, 1969–1989 (2022).
An, Y. & Duan, H. The role of m6A RNA methylation in cancer metabolism. Mol. Cancer 21, 14 (2022).
Zeng, D. et al. Tumor Microenvironment Characterization in Gastric Cancer Identifies Prognostic and Immunotherapeutically Relevant Gene Signatures. Cancer Immunol. Res. 7, 737–750 (2019).
Sotiriou, C. et al. Gene Expression Profiling in Breast Cancer: Understanding the Molecular Basis of Histologic Grade To Improve Prognosis. JNCI J. Natl. Cancer Inst. 98, 262–272 (2006).
Palmisani, F. et al. Systematic review of the immunological landscape of Wilms tumors. Mol. Ther. - Oncolytics 22, 454–467 (2021).
Chen, Y., Lin, Y., Shu, Y., He, J. & Gao, W. Interaction between N6-methyladenosine (m6A) modification and noncoding RNAs in cancer. Mol. Cancer 19, 94 (2020).
Yi, L., Wu, G., Guo, L., Zou, X. & Huang, P. Comprehensive Analysis of the PD-L1 and Immune Infiltrates of m6A RNA Methylation Regulators in Head and Neck Squamous Cell Carcinoma. Mol. Ther. - Nucleic Acids 21, 299–314 (2020).
Gao, Y. et al. m6A Modification Prevents Formation of Endogenous Double-Stranded RNAs and Deleterious Innate Immune Responses during Hematopoietic Development. Immunity 52, 1007-1021.e8 (2020).
Wang, T., Kong, S., Tao, M. & Ju, S. The potential role of RNA N6-methyladenosine in Cancer progression. Mol. Cancer 19, 88 (2020).
Visvanathan, A. et al. Essential role of METTL3-mediated m6A modification in glioma stem-like cells maintenance and radioresistance. Oncogene 37, 522–533 (2018).
Liu, J. et al. m6A mRNA methylation regulates AKT activity to promote the proliferation and tumorigenicity of endometrial cancer. Nat. Cell Biol. 20, 1074–1083 (2018).
Wilkerson, M. D. & Hayes, D. N. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics 26, 1572–1573 (2010).
Peña-Romero, A. C. & Orenes-Piñero, E. Dual Effect of Immune Cells within Tumour Microenvironment: Pro- and Anti-Tumour Effects and Their Triggers. Cancers 14, 1681 (2022).
Zeng, Z. et al. Immune and stromal scoring system associated with tumor microenvironment and prognosis: a gene-based multi-cancer analysis. J. Transl. Med. 19, 330 (2021).
Eulberg, D., Frömming, A., Lapid, K., Mangasarian, A. & Barak, A. The prospect of tumor microenvironment-modulating therapeutical strategies. Front. Oncol. 12, 1070243 (2022).
Zhang, F. et al. Crosstalk among m6A RNA methylation, hypoxia and metabolic reprogramming in TME: from immunosuppressive microenvironment to clinical application. J. Hematol. Oncol.J Hematol Oncol 15, 84 (2022).
Tian, X.-M. et al. Immune-related gene signature associates with immune landscape and predicts prognosis accurately in patients with Wilms tumour. Front. Immunol. 13, 920666 (2022).
Rossetti, R. A. M. et al. B lymphocytes can be activated to act as antigen presenting cells to promote anti-tumor responses. PLOS ONE 13, e0199034 (2018).
Eosinophils: Methods and Protocols. vol. 2241 (Springer US, 2021).
Callahan, M. K., Postow, M. A. & Wolchok, J. D. Targeting T Cell Co-receptors for Cancer Therapy. Immunity 44, 1069–1078 (2016).
Xie, D. et al. Immune checkpoint inhibitor plus tyrosine kinase inhibitor for unresectable hepatocellular carcinoma in the real world. Ann. Transl. Med.9, 652–652 (2021).

No competing interests reported.

Download PDF

Journal Publication

published 21 Sep, 2023

Read the published version in BMC Medical Genomics →

Editorial decision: Major revision
11 Jul, 2023
Reviews received at journal
07 Jul, 2023
Reviews received at journal
22 Jun, 2023
Reviewers agreed at journal
11 Jun, 2023
Reviewers agreed at journal
12 Apr, 2023
Reviewers agreed at journal
10 Apr, 2023
Reviewers invited by journal
25 Mar, 2023
Editor assigned by journal
25 Mar, 2023
Editor invited by journal
02 Mar, 2023
Submission checks completed at journal
02 Mar, 2023
First submitted to journal
13 Feb, 2023

You are reading this latest preprint version

Identification of the Expression Patterns and Potential Prognostic Role of m6A-RNA Methylation Regulators in Wilms Tumor

Status:

Journal Publication

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

Discussion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1