PIK3CD correlates with prognosis, epithelial–mesenchymal transition and tumor immune infiltration in breast carcinoma

Background Breast carcinoma (BRCA) is one of the most common, fatal, and aggressive cancers, with increasing morbidity that has a major impact on human health. PIK3CD appears to have important roles in the beginning and advancement of various forms of human cancer, according to mounting data. However,the particular role and mechanism of PIK3CD in BRCA remains not fully identified. Methodology The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/), Genotype-Tissue Expression (GTEx) data and the UCSC Xena browser (https://xenabrowser.net) data were used in this study’s initial pan-cancer analysis of PIK3CD expression and prognosis. Circular RNAs (circRNAs) that regulated the expression of PIK3CD were subsequently found using a combination of in silico investigations of expression, correlation, and survival. Measurements of PIK3CD expression and an analysis of the in vitro function of PIK3CD in BRCA cells were performed using real-time RT-PCR, Western blotting and Transwell assays. Results In BRCA GLI2, RAB32, LAMB1, MGAT2, ITGA8, CHF, COL6A3 and PRRX1-miR-30b-5p axis was identified as the most likely upstream CircRNA-related route of PIK3CD. PIK3CD was correlated with the expression of EMT markers. The PIK3CD cDNA improved the capacity for invasion and migration. The expression of PIK3CD was linked to some of the m1A/m5C/m6A regulators. Additionally, it was discovered that the expression of PIK3CD was found to be highly connected to the expression of immunological checkpoints, immune cell biomarkers, and tumor immune cell invasion. Conclusions Our findings reveal that PIK3CD expression is associated with prognosis, EMT, and tumor immune infiltration in BRCA patients.


Introduction
The most frequent primary cancer and the second greatest cause of deaths due to cancer in women in the US is breast cancer (BRCA).Breast cancer is still difficult to prevent globally because it is a multi-step procedure that involves many types of cells.One of the best ways to avoid breast cancer is through early diagnosis [1].Several risk factors, including the immune system, have been connected to the beginning and development of BRCA.The prognosis for people with BRCA remains poor despite the enormous advances made in the fields of diagnosis, treatment, and prognosis.Thus, it is crucial to find viable prognostic biomarkers or create efficient treatment targets for BRCA.
The phosphoinositide 3-kinase (PI3K) is one of the most prevalent regulatory mechanisms in a range of human malignancies.Based on their substrate specificity and mammalian structural characteristics, PI3K is split into three types.The PI3K subfamily of class I appears to be the one most associated with human malignancies.According to their adaptors, Class I PI3K is further separated into subclasses IA and IB.Class IA catalytic isoforms p110α, p110β and p110δ are encoded three distinct genes termed PIK3CA, PIK3CB and PIK3CD.The functions and fundamental mechanisms involving PIK3CD in cancer, however, are poorly understood [2].
The epithelial-mesenchymal transition (EMT) can promote metastasis in a variety of cancer types.EMT was linked to the exclusion of immune cells important in the immune reaction to cancer as well as the expression of several immunosuppressive cytokines.There may be a link between EMT and tumor immune infiltration in cancer [3].
In our study, we first examined the expression and survival of PIK3CD in diverse human cancers.Following that, the CircRNA-related regulation of PIK3CD in BRCA was investigated.Subsequently, we examined the connection between PIK3CD expression and influx of immunological cells, biological markers, or immunological checkpoints in BRCA.These results suggest that PIK3CD is associated with tumor immune infiltration and prognosis in BRCA individuals.

Analysis of the StarBase database
StarBase is a database for miRNA research (http:// starb ase.sysu.edu.cn/).In BRCA, starBase was used to conduct an expression analysis of correlation for miRNA-PIK3CD, circRNA-miRNA, or circRNA-PIK3CD.StarBase was also used to examine miRNA-30b-5p expression in BRCA and normal controls.Furthermore, starBase was employed to anticipate candidate circRNAs that might bind to miRNA-30b-5p.

Prediction of candidate miRNAs
Several target gene prediction programs, including PITA, RNA22, miRmap, microT, miRanda, PicTar and Tar-getScan, anticipated upstream binding miRNAs of PIK3CD.The following analyses included only anticipated miRNAs that showed often in more than two programs, as previously noted.These expected miRNAs were considered candidate miRNAs of PIK3CD.

Cell culture
The Shanghai Institute of Cell Biology's and the Chinese Academy of Sciences' Cell Banks provided the breast cancer MCF-7 cells and MDA-MB-231 cells for this study.The cell lines were grown in Roswell Park Memorial Institute (RPMI)-1640 medium (Invitrogen) containing 10% fetal bovine serum (Invitrogen) and antibiotics (100 g/mL streptomycin and 100 U/mL penicillin) at 37 °C and 5% CO2 in a humidified tissue culture chamber.Cells from the period of logarithmic growth were employed in further tests.

Lentivirus production and transduction
The NC (negative control) mRNA precursor sequences and the PIK3CD sequences (cDNA) were amplified and inserted into the pLV3 vector, both from GeneChem in Shanghai, China.For virus production, HEK293T cells were transfected 1 3 with the packaging vectors pGag/Pol, pRev and pVSV-G in combination with the recombinant lentiviral plasmids.After 48 h of transfection, the lentivirus-containing medium was recovered and filtered.The viral supernatant containing polybrene was then introduced to the breast cancer MDA-MB-231 cells and MCF-7 cells to cause infection.The transduced cells were then chosen for 14 days of puromycin to identify those that were consistently producing the constructs.It was determined how much PIK3CD was present in these stable cells using quantitative real-time PCR [2].

Transwell assay
The migration experiment was conducted by seeding cells in the upper chambers without matrigel using inserts (Costar, Cambridge, MA).In invasion assay (BD Bioscience), the inserts were completely covered with matrigel.Cells were resuspended in FBS-free culture and introduced to the upper Transwell chamber, including Matrigel, while the lower chamber received 800 µl of culture medium with 10% FBS.
After the cells were incubated at 37 ℃ for 24 h, the inserts were fixed with methanol and stained with crystal violet at room temperature for 20 min.The cells were taken away from the upper side of the filter membrane, and the invading cells were stained for 10-15 min at room temperature with a 0.1% crystal violet solution.An inverted microscope was used to take the pictures, and ten random fields' worth of cells were counted.Displayed were the statistical analyses and illustrative pictures.Each group's experiments were carried out in triplicate.

TIMER database analysis
TIMER (https:// cistr ome.shiny apps.io/ timer/) is a web server for extensive research of tumor-infiltrating immune cells.TIMER was employed to look into the link between PIK3CD expression and infiltration of immune cells or immune-related checkpoint expression in BRCA.A p value of 0.05 was regarded as statistically significant.

Statistical analysis
The aforementioned internet database generated the statistical analysis for this study automatically.Data from three or more studies were analyzed using R 4.2.1 and given as mean ± SD.The statistical difference was ascertained using either the Student's t test or a one-way ANOVA.The threshold for statistical significance was set at a p value of 0.05 or logrank p value of 0.05.For the results, *p < 0.05,**p < 0.01 and ***p < 0.001 were regarded as statistically significant.
There was no statistical relevance of PIK3CD in predicting patient outcome in other kinds of cancer (Fig. 2E-H).In patients with BRCA, PIK3CD could be used as a predictive biomarker.

Prediction and evaluation of PIK3CD's upstream miRNAs
MiRNAs are widely known to be responsible for gene expression regulation.To determine whether PIK3CD was controlled by miRNAs, we predicted upstream miRNAs that could potentially bind to PIK3CD and eventually discovered 6 miRNAs.Based on the competing endogenous RNA (ceRNA) hypothesis, miRNA regulates target gene expression and there should be an adverse correlation between miRNA and PIK3CD.As a result, the expression correlation analysis was carried out.In BRCA, PIK3CD was significantly negatively linked with hsa-miR-30b-5p(|R|= 0.169) and favorably correlated with hsa-miR-30e-5p and hsa-miR-7-5p, as shown in Table 1.PIK3CD and hsa-miR-7-5p have no statistically significant expression correlations.Finally, hsa-miR-30b-5p expression and prognostic significance in BRCA were determined.As shown in Fig. 3A and B, hsa-miR-30b-5p was significantly downregulated in BRCA patients, and its overexpression was associated with a better prognosis.All of these findings point to miRNA-30b-5p as the most likely regulating miRNA of PIK3CD in BRCA.

Prediction and analysis of upstream circRNAs of hsa-miR-30b-5p
Circular RNAs (circRNAs) have been linked to the progression of cancer [4].We used bioinformatics to analyze and identify the circRNAs that might probably mediate PIK3CD expression in BRCA.The upstream circRNAs of hsa-miR-30b-5p were then predicted utilizing the starBase database.We investigated the top ten circRNAs.The levels of expression of these circRNAs were subsequently shown to be negatively related to the levels of hsa-miR-30b-5p, which were found to be negatively related to the expression of PIK3CD in BRCA.All of these circRNAs' expression was linked to PIK3CD (Tables 2 and 3).The starBase database was then used to analyze the amounts of expression of these circRNAs in BRCA.Only CD99 was not significantly regulated in BRCA as compared to normal controls, as seen in Fig. 4. The prognostic values of the ten circRNAs in BRCA were then assessed.BRCA patients with higher levels of CFH, GLI2, RAB32, LAMB1, MGAT2, CD99 and ITGA8 expression had a better prognosis,while those of COL6A3 and PRRX1 had a worse one, as illustrated in Fig. 5.

PIK3CD correlates with EMT, migration and invasion in BRCA
The link between PIK3CD and EMT variables was then investigated.In the majority of the 5 cancer types, PIK3CD expression was linked favorably with the expression of mesenchymal markers such as VIM (Vimentin), TWIST1, SNAI1 (SNAIL), SNAI2 (SLUG), and CDH2 (Fig. 6A).Furthermore, PIK3CD was found to be strongly inversely related to the epithelial marker CDH1 (E-Cadherin).The correlation heatmaps revealed a strong relationship between PIK3CD and the EMT variables.These findings showed that PIK3CD may be involved in EMT.
We studied the impact of PIK3CD on cell migration, invasion, and EMT of BRCA cells after verifying its expression and prognostic value for BRCA.MCF-7 BRCA cell lines were transfected with PIK3CD cDNA.To evaluate the effect of PIK3CD on cell migration and invasion, Transwell Assay was utilized.BRCA cells that received PIK3CD cDNA displayed higher migration and invasion abilities (Fig. 6B-D).After PIK3CD cDNA transfection, the expression of EMT-related markers was evaluated using Western blot.E-cadherin expression was downregulated in the PIK3CD cDNA group, while VIM and SNAI1 expression were elevated.
These findings demonstrated that PIK3CD cDNA increased cell EMT, migration, and invasion, indicating that PIK3CD plays an active role in BRCA cell activities.

PIK3CD correlates with influx of immunological cells in BRCA
The TIMER database was used to demonstrate the link between PIK3CD expression and cell infiltration levels in BRCA cells.The amount of immune cell infiltration significantly changed with different PIK3CD copy counts in BRCA, as illustrated in Fig. 7A.Somatic copy number alterations (SCNA) of PIK3CD were closely correlated with the abundance of immune infiltration.Consequently, the relationship between PIK3CD expression and immune cell infiltration was assessed.According to Fig. 7B, the expression of PIK3CD was highly positively linked with every Clinical treatment sensitivity and disease diagnosis are significantly impacted by the microenvironment, which includes immune cells, extracellular matrix, inflammatory chemicals, and other growth factors.Using the CIBER-SORT technique, the fraction of 22 immune cells was estimated using a histogram in this study.The results showed that the tumor group had considerably larger proportions of most immune cells than the control groups, with the exception of some immune cells such as CD56 dim natural killer cells in BRCA,LUAD, and STAD (Fig. 7C).

Correlation of PIK3CD expression and immune cell indicators in BRCA
The GEPIA database was used to examine the expression connection of PIK3CD with immune cell indicators in BRCA, able to better investigate the function of PIK3CD in tumor immunity.According to the findings, PIK3CD expression was strongly positively linked with immunity-associated gene marker expression, such as biomarkers of B cells including CD19 and CD79A, and biomarkers of dendritic cells in BRCA (Table 4, *p value < 0.05; **p value < 0.01; ***p value < 0.001).These findings revealed that increased PIK3CD levels were related to increase immunological infiltration as well as immune exhaustion, suggesting that PIK3CD could be a potential immunotherapeutic target for BRCA treatment.These results help to some extent establish the favorable relationship between PIK3CD and immune cell infiltration.

Link between PIK3CD and immunological checkpoints in BRCA
Immune checkpoints like PD1/PD-L1 and CTLA-4 are essential for tumor immunity escape.Given the putative oncogenic involvement of PIK3CD in BRCA, the association of PIK3CD with PD1, PD-L1,or CTLA-4 was investigated.As demonstrated in Fig. 8A, PIK3CD expression was strongly positively linked with PD1, PD-L1 and CTLA-4 in BRCA,which was modified by purity.We also found that PD1, PD-L1, CTLA-4 and PIK3CD in BRCA had a highly significant relationship, similar to what we found in the analysis of the TIMER data.Our findings show that PIK3CD-mediated tumorigenesis of BRCA may entail tumor immune evasion.

Relationships between PIK3CD expression, methylation level and the amount of lymphocytes that infiltrate tumors
We discovered that PIK3CD expression had a positive relationship with the abundance of 28 different types of tumorinfiltrating lymphocytes in the TISIDB database (Fig. 8B, Table 5).We investigated the association between PIK3CD   methylation and tumor-infiltrating lymphocyte and discovered that it was strongly inversely linked with tumor-infiltrating lymphocyte abundance (Fig. 8B, Table 5).The results indicated that PIK3CD level was favorably related to the high number of lymphocytes infiltrating tumors, but PIK3CD methylation was negatively associated, and it might be used as one of the immunotherapeutic targets for BRCA therapy.We also discovered links between PIK3CD methylation and clinicopathological characteristics in BRCA In TNBC and HER2-Positive breast cancer tissues, the promoter methylation level of PIK3CD was lower than in Luminal (Fig. 8B).Thus, reduced PIK3CD methylation could be an underlying indication reflecting clinical BRCA features.

Identification of the m1A/m5C/m6A regulators relevant to PIK3CD in BRCA
The pan-cancer correlations between PIK3CD and the expression levels of 40 kinds of the m1A/m5C/m6A regulators in 5 cancer types were shown in Fig. 9A.There was a strong association (P < 0.05) between the expression levels of the 41 different m1A/m5C/m6A regulators and PIK3CD in 5 different cancer types.Furthermore,in BRCA and LUAD, significant co-expression of PIK3CD with most regulator genes was detected.In LIHC, expression of PIK3CD was negatively correlated with some molecules such as ZC3H13, YTHDF3 and TRMT61A, though for some of them not to a significant degree.It is also interesting that expression of PIK3CD was negatively correlated with TRMT61B, FMR1, FTO and METTL14 in LIHC,but the levels of PIK3CD were positively associated with them in BRCA (Fig. 9A).

Discussion
BRCA is still known for having a terrible prognosis.The discovery of useful prognostic biomarkers or the creation of useful therapy targets may result from an understanding of the molecular process of BRCA.Numerous human cancers' genesis and development have been shown to be significantly influenced by PIK3CD, among other things.However, our understanding of PIK3CD in BRCA remains limited and has to be expanded.
In this investigation, we first used data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) to do a pan-cancer analysis of PIK3CD expression.PIK3CD was significantly downregulated in 10 tumor types, including BRCA, COAD, LIHC, LUAD and LUSC (Fig. 1).As an intrinsic cancer-causing driver, PIK3CD was found in primary fibroblasts but was scarcely detectable in breast cancer cell lines [5].According to the survival analysis for the cancer types of interest, patients with high PIK3CD expression had a favourable outlook for survival (OS) in LUAD and BRCA and a poor outcome in PRAD and LUSC.In individuals with laryngeal squamous cell carcinoma, PIK3CD expression was elevated and could predict a bad prognosis [6].Fibroblast-expressed PIK3CD was a distinct predictive factor for overall and disease-free survival, indicating it as a therapeutic target for TNBC [5].
We used PicTar, RNA22, PITA, miRanda, microT, miRmap, and TargetScan to identify candidate miRNAs that may potentially bind to PIK3CD in order to examine the upstream regulatory miRNAs of PIK3CD.Six miRNAs were discovered in the end (Table 1).In laryngeal squamous cell carcinoma, the miR-30c-5p/PIK3CD axis is related to malignant properties and predicts a prognosis [6].The bulk of these miRNAs have been identified as tumor-suppressive miRNAs in BRCA patients.Following correlation, expression, and survival analyses, miR-30b-5p was identified as the most likely upstream tumor suppressive miRNA of PIK3CD.It is reported that miR-30b-5p inhibits lung cancer cell growth and increases cisplatin sensitivity [7].
It has been extensively established that miRNAs and circRNAs contribute to gene expression regulation by communicating with one another via the ceRNA mechanism [8,9].Based on the ceRNA hypothesis, the putative circRNAs should be highly positively associated with the level of PIK3CD in BRCA.
PIK3CD expression correlated favorably with the expression of mesenchymal markers and was strongly inversely related to the epithelial marker.The correlation heatmaps revealed a strong relationship between PIK3CD and the EMT variables.PIK3CD cDNA increased cell EMT, migration, and invasion, indicating that PIK3CD played an active role in BRCA cell activities (Fig. 6).These findings showed that PIK3CD may be involved in EMT.EMT was associated with immune infiltration [18] and EMT-related genes were related to prognosis.PIK3CD significantly boosted colorectal cancer cell motility, invasion, and proliferation vitro, as well as tumor development in vivo [2].
To our interest, PIK3CD enhanced the migration, invasion, and EMT of BRCA cells,but PIK3CD might be used as a favorable prognostic biomarker in BRCA patients.PIK3CD was shown to be related to the expression of SNAI1.Luminal A tumors with high SNAI1 expression had a higher OS marker, while luminal B, Her2-positive, and basal-like malignancies had a poorer OS marker [19].As a result, we hypothesized that PIK3CD was associated with a favorable prognostic biomarker for OS in the luminal A subtype.The most common subtype of breast cancer is luminal/estrogen receptor-positive (ER +) breast cancer, which includes luminal A and luminal B subtypes) [20].
Several studies have shown that immune cell infiltration in tumors can affect the efficacy of chemotherapy [21], radiotherapy [22], or immunotherapy [23], as well as cancer patients' prognosis [24].PIK3CD correlated with immune cell infiltration in BRCA (Figs. 7-9, Table 4).In BRCA a substantial positive relationship was seen between PIK3CD and B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells (Fig. 7B, C, Table 4).Furthermore, PIK3CD was found to be significantly related to biomarkers of these infiltrated immune cells.These results suggested that tumor immune infiltration might contribute to PIK3CD-mediated carcinogenic functions in BRCA.
The effectiveness of immunotherapy is also influenced by the expression of immunological checkpoints and the number of immune cells entering the tumor microenvironment [25].As a result, we investigated the link between PIK3CD and immunological checkpoints.The findings showed that high PIK3CD expression was highly connected to PD1, PD-L1, or CTLA-4 in BRCA, implying that highly expressed PIK3CD may improve the efficacy of immunotherapy in BRCA (Fig. 8A).
We discovered that PIK3CD expression had a positive relationship with the abundance of 28 different types of tumor-infiltrating lymphocytes and that PIK3CD methylation was strongly inversely linked with tumor-infiltrating lymphocyte abundance in the TISIDB database.We also discovered links between PIK3CD methylation and clinicopathological characteristics in BRCA patients (Fig. 8B, Table 5).Reduced PIK3CD methylation could be an underlying indication reflecting clinical BRCA features.Tumor infiltrating lymphocytes (TIL) are critical in influencing treatment response and enhancing clinical outcomes in all subtypes of breast cancer [26].Both HER2-positive and TNBCs have a better prognosis and treatment response when there are many of tumor-infiltrating lymphocytes (TILs) [27].Tumors with low T cell infiltration are more resistant to immunogenic chemotherapies and checkpoint inhibition than cancers with high T cell infiltration [28].Recently, DNA methylation profiles were discovered to be an indication of the tumor immune microenvironment [29].PIK3CD expression was found to be favorably related to the infiltration of immune cells and gene-marker levels, which could be linked to methylation-related changes.
RNA modifications such as m6A, m5C and m1A have been shown in recent years to influence the immune response and cancer [30].Correlation analyses found a high correlation between PIK3CD and the m1A/m5C/m6A regulators in 5 cancer types (Fig. 9).In hepatocellular carcinoma, m6A/m5C/m1A regulated genes were linked to a poor prognosis and an immunological microenvironment [31].
In the present study, we systematically investigated PIK3CD's function in BRCA.We identified that PIK3CD had the potential to operate as a detection and prognosis biomarker for BRCA by examining its underlying mechanism.We found an upstream regulatory mechanism for PIK3CD in BRCA, which included GLI2, RAB32, LAMB1, MGAT2, ITGA8, CHF, COL6A3 and PRRX1/miR-30b-5p/PIK3CD. Furthermore, our current findings suggest that PIK3CD may have an important function by boosting EMT, immune checkpoint expression and tumor immune cell infiltration.However, these results require confirmation in the future through much more fundamental research and substantial clinical trials.To date, there have been relatively few studies on PIK3CD in BRCA.Our comprehensive examination of PIK3CD and its conceivable mechanisms may serve as the basis for further investigation.

Fig. 6
Fig. 6 PIK3CD correlates with EMT in BRCA.A Correlation of PIK3CD expression with expression of EMT biomarkers across 5 cancer types.B The levels of PIK3CD mRNA in BRCA cell lines were detected by using qRT-PCR.a MDA-MB-231 cells.b MCF-7 cells.C Upregulation of PIK3CD increased MDA-MB-231 BRCA cell migration and invasion.Transwell assay of cell migration ability following PIK3CD cDNA treatment.Transwell assessment of cell invasion capabilities following PIK3CD cDNA transfection (a and b).(c) Western blot examination of EMTrelated marker expression alterations following PIK3CD cDNA transfection.The statistical analysis was displayed.D Upregulation of PIK3CD enhanced MCF-7 BRCA cell motility and invasion.Transwell experiment of cell migration ability following PIK3CD cDNA treatment.Transwell test of cell invasion capabilities following PIK3CD cDNA transfection (a and b).(c) Western blot study of EMT-related marker expression following PIK3CD cDNA transfection.The statistical analysis was displayed

Fig. 7
Fig. 7 The association between immune cell infiltration and PIK3CD levels in BRCA.A The level of invasion of various immune cells in BRCA with varying copy levels of PIK3CD from the TIMER database.B The level of PIK3CD expression in BRCA is correlated with B cell (b), CD8 + T cell (c), CD4 + T cell (d), macrophage (e), neutrophil (f ), or dendriticcell (g) infiltration.C Correlation of PIK3CD with immune cells 5 cancer types.The CIBERSORT technique was used in this work to estimate the fraction of 22 immune cells

Fig. 8
Fig. 8 Relationship between and immune checkpoints/promoter methylation level of PIK3CD in BRCA.A Correlation of PIK3CD expression in BRCA with PD-1, PD-L1, and CTLA-4 expression.(a) TIMER-adjusted Spearman association of PIK3CD with PD-1 expression in BRCA.(b) The GEPIA database determined the expression connection of PIK3CD with PD1 in BRCA.(c) PIK3CD-PD-L1 TIMER-adjusted Spearman correlation in BRCA.(d) The GEPIA database determined the expression connection of PIK3CD with PD-L1 in BRCA.(e) TIMER Spearman association of PIK3CD with CTLA-4 expression in BRCA corrected for purity.(f ) The GEPIA database determined the expression connection of PIK3CD with CTLA-4 in BRCA.B The correlations between PIK3CD expression, methylation level, and Tumor-infiltrating lymphocyte (TIL) quantity.(a) Relationships between 28 different types of tumor-infiltrating lymphocytes (TILs) and PIK3CD expression in several tumors.(b) Correlations between 28 different types of TILs and PIK3CD methylation in several tumors.(c) Correlations between the promoter methylation level of PIK3CD and the clinicopathological traits of BRCA patients.UALCAN database was used to examine the relative promoter methylation level of PIK3CD based on major subclasses

Fig. 9
Fig. 9 The m1A/m5C/m6A regulators relevant to PIK3CD in BRCA.A Relations between the levels of PIK3CD and the expressions of 41 kinds of the m1A/m5C/m6A regulators in various cancers.B The m1A/m5C/m6A regulator expression variations between tumor and normal samples.The vioplot shows the top 10 genes (|R|) that are significantly differentially expressed in breast tumors and normal tissues.C The overall survival (OS) analysis for the top 10 m1A/m5C/m6A regulators significantly differentially expressed (|R|) in both breast tumor and normal samples determined by Kaplan-Meier Plotter

Table 1
The starBase database