Maternal prednisone exposure during pregnancy elevates susceptibility to osteoporosis in female offspring: The role of mitophagy/FNDC5 alteration in skeletal muscle

Maternal exposure to glucocorticoids has been associated with adverse outcomes in offspring. However, the consequences and mechanisms of gestational exposure to prednisone on susceptibility to osteoporosis in the offspring remain unclear. Here, we found that gestational prednisone exposure enhanced susceptibility to osteoporosis in adult mouse offspring. In a further exploration of myogenic mechanisms, results showed that gestational prednisone exposure down-regulated FNDC5/irisin protein expression and activation of OPTN-dependent mitophagy in skeletal muscle of adult offspring. Additional experiments elucidated that activated mitophagy significantly inhibited the expression of FNDC5/irisin in skeletal muscle cells. Likewise, we observed delayed fetal bone development, downregulated FNDC5/irisin expression, and activated mitophagy in fetal skeletal muscle upon gestational prednisone exposure. In addition, an elevated total m6A level was observed in fetal skeletal muscle after gestational prednisone exposure. Finally, gestational supplementation with S-adenosylhomocysteine (SAH), an inhibitor of m6A activity, attenuated mitophagy and restored FNDC5/irisin expression in fetal skeletal muscle, which in turn reversed fetal bone development. Overall, these data indicate that gestational prednisone exposure increases m6A modification, activates mitophagy, and decreases FNDC5/irisin expression in skeletal muscle, thus elevating osteoporosis susceptibility in adult offspring. Our results provide a new perspective on the earlier prevention and treatment of fetal-derived osteoporosis.

Co-transcriptional R-loops-mediated epigenetic regulation drives growth retardation and docetaxel chemosensitivity enhancement in advanced prostate cancer

R-loops are prevalent three-stranded nucleic acid structures, comprising a DNA-RNA hybrid and a displaced single-stranded DNA, that frequently form during transcription and may be attributed to genomic stability and gene expression regulation. It was recently discovered that RNA modification contributes to maintain the stability of R-loops such as N6-methyladenosine (m6A). Yet, m6A-modified R-loops in regulating gene transcription remains poorly understood. Here, we demonstrated that insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs) recognize R-loops in an m6A-dependent way. Consequently, IGF2BPs overexpression leads to increased overall R-loop levels, cell migration inhibition, and cell growth retardation in prostate cancer (PCa) via precluding the binding of DNA methyltransferase 1(DNMT1) to semaphorin 3 F (SEMA3F) promoters. Moreover, the K homology (KH) domains of IGF2BPs are required for their recognition of m6A-containing R-loops and are required for tumor suppressor functions. Overexpression of SEMA3F markedly enhanced docetaxel chemosensitivity in prostate cancer via regulating Hippo pathway. Our findings point to a distinct R-loop resolution pathway mediated by IGF2BPs, emphasizing the functional importance of IGF2BPs as epigenetic R-loop readers in transcriptional genetic regulation and cancer biology.The manuscript summarizes the new role of N6-methyladenosine in epigenetic regulation, we introduce the distinct R-loop resolution mediated by IGF2BP proteins in an m6A-dependent way, which probably lead to the growth retardation and docetaxel chemotherapy resistance in prostate cancer. Moreover, our findings first emphasized the functional importance of IGF2BPs as epigenetic R-loop readers in transcriptional genetic regulation and cancer biology. In addition, our research provides a novel RBM15/IGF2BPs/DNMT1 trans-omics regulation m6A axis, indicating the new crosstalk between RNA m6A methylation and DNA methylation in prostate cancer.

ALKBH5 SUMOylation-mediated FBXW7 m6A modification regulates alveolar cells senescence during 1-nitropyrene-induced pulmonary fibrosis

Our previous study revealed that 1-nitropyrene (1-NP) exposure evoked pulmonary fibrosis in mice. However, the exact mechanism remained elusive. We found that 1-NP induced telomere damage and cellular senescence in mice lungs, and two alveolar epithelial cells lines. 1-NP downregulated telomere repeat binding factor 2 (TRF2), and upregulated FBXW7. Mechanistically, 1-NP-caused TRF2 ubiquitination and proteasomal degradation depended on E3 ubiquitin ligase activity of FBXW7. Moreover, 1-NP upregulated FBXW7 m6A modification via an ALKBH5-YTHDF1-dependent manner. Further analysis suggested 1-NP promoted ALKBH5 SUMOylation and subsequent proteasomal degradation. Additionally, 1-NP evoked mitochondrial reactive oxygen species (mtROS) overproduction. Mito-TEMPO, a mitochondrial-targeted antioxidant, mitigated 1-NP-caused mtROS overproduction, ALKBH5 SUMOylation, FBXW7 m6A modification, TRF2 degradation, cellular senescence, and pulmonary fibrosis. Taken together, mtROS-initiated ALKBH5 SUMOylation and subsequent FBXW7 m6A modification is indispensable for TRF2 degradation and cellular senescence in alveolar epithelial cells during 1-NP-induced pulmonary fibrosis. Our study provides target intervention measures towards 1-NP-evoked pulmonary fibrosis.

N6-methyladenosine-modified circ_0000337 sustains bortezomib resistance in multiple myeloma by regulating DNA repair

Studies have shown that bortezomib resistance in multiple myeloma (MM) is mediated by the abnormalities of various molecules and microenvironments. Exploring these resistance mechanisms will improve the therapeutic efficacy of bortezomib. In this study, bone marrow tissues from three patients with MM, both sensitive and resistant to bortezomib, were collected for circRNA high-throughput sequencing analysis. The relationship between circ_0000337, miR-98-5p, and target gene DNA2 was analyzed by luciferase detection and verified by RT-qPCR. We first found that circ_0000337 was significantly upregulated in bortezomib-resistant MM tissues and cells, and overexpression of circ_0000337 could promote bortezomib resistance in MM cells. circ_0000337 may act as a miR-98-5p sponge to upregulate DNA2 expression, regulate DNA damage repair, and induce bortezomib resistance. Furthermore, it was determined that the increased circ_0000337 level in bortezomib-resistant cells was due to an increased N6-methyladenosine (m6A) level, resulting in enhanced RNA stability. In conclusion, the m6A level of circ_0000337 and its regulation may be a new and potential therapeutic target for overcoming bortezomib resistance in MM.

Oxidative stress-mediated activation of FTO exacerbates impairment of the epithelial barrier by up-regulating IKBKB via N6-methyladenosine-dependent mRNA stability in asthmatic mice exposed to PM2.5

In order to comprehend the underlying mechanisms contributing to the development and exacerbation of asthma resulting from exposure to fine particulate matter (PM2.5), we established an asthmatic model in fat mass and obesity-associated gene knockdown mice subjected to PM2.5 exposure. Histological analyses using hematoxylin-eosin (HE) and Periodic Acid-Schiff (PAS) staining revealed that the down-regulation of the fat mass and obesity-associated gene (Fto) expression significantly ameliorated the pathophysiological alterations observed in asthmatic mice exposed to PM2.5. Furthermore, the down-regulation of Fto gene expression effectively attenuated damage to the airway epithelial barrier. Additionally, employing in vivo and in vitro models, we elucidated that PM2.5 modulated FTO expression by inducing oxidative stress. Asthmatic mice exposed to PM2.5 exhibited elevated Fto expression, which correlated with increased levels of reactive oxygen species. Similarly, when cells were exposed to PM2.5, FTO expression was up-regulated in a ROS-dependent manner. Notably, the administration of N-acetyl cysteine successfully reversed the PM2.5-induced elevation in FTO expression. Concurrently, we performed transcriptome-wide Methylated RNA immunoprecipitation Sequencing (MeRIP-seq) analysis subsequent to PM2.5 exposure. Through the implementation of Gene Set Enrichment Analysis and m6A-IP-qPCR, we successfully identified inhibitor of nuclear factor kappa B kinase subunit beta (IKBKB) as a target gene regulated by FTO. Interestingly, exposure to PM2.5 led to increased expression of IKBKB, while m6A modification on IKBKB mRNA was reduced. Furthermore, our investigation revealed that PM2.5 also regulated IKBKB through oxidative stress. Significantly, the down-regulation of IKBKB effectively mitigated epithelial barrier damage in cells exposed to PM2.5 by modulating nuclear factor-kappa B (NF-κB) signaling. Importantly, we discovered that decreased m6A modification on IKBKB mRNA facilitated by FTO enhanced its stability, consequently resulting in up-regulation of IKBKB expression. Collectively, our findings propose a novel role for FTO in the regulation of IKBKB through m6A-dependent mRNA stability in the context of PM2.5-induced oxidative stress. Therefore, it is conceivable that the utilization of antioxidants or inhibition of FTO could represent potential therapeutic strategies for the management of asthma exacerbated by PM2.5 exposure.

Revealing PPP1R12B and COL1A1 as piRNA pathway genes contributing to abdominal aortic aneurysm through integrated analysis and experimental validation

BACKGROUND

Abdominal aortic aneurysm (AAA) is a permanent dilation of the abdominal aorta, with a high mortality rate when rupturing. Although lots of piRNA pathway genes (piRPGs) have recently been linked to both neoplastic and non-neoplastic illnesses, their role in AAA is still unknown. Utilizing integrative bioinformatics methods, this research discovered piRPGs as biomarkers for AAA and explore possible molecular mechanisms.

METHODS

The datasets were obtained from the Gene Expression Omnibus and piRPGs were identified from the Genecards database. The "limma" and "clusterProfiler" R-packages were used to discover differentially expressed genes and perform enrichment analysis, respectively. Hub piRPGs were further filtered using least absolute shrinkage and selection operator regression, random forests, as well as receiver operating characteristic curve. Additionally, multi-factor logistic regression (MLR), extreme gradient boosting (XGboost), and artificial neural network (ANN) were employed to construct prediction models. The relationship between hub piRPGs and immune infiltrating cells and sgGSEA were further studied. The expression of hub piRPGs was verified by qRT-PCR, immunohistochemistry, and western blotting in AAA and normal vascular tissues and analyzed by scRNA-seq in mouse AAA model. SRAMP and cMAP database were utilized for the prediction of N6-methyladenosine (m6A) targets therapeutic drug.

RESULTS

34 differentially expressed piRPGs were identified in AAA and enriched in pathways of immune regulation and gene silence. Three piRPGs (PPP1R12B, LRP10, and COL1A1) were further screened as diagnostic genes and used to construct prediction model. Compared with MLR and ANN, Xgboost showed better predictive ability, and PPP1R12B might have the ability to distinguish small and large AAA. Furthermore, the expression levels of PPP1R12B and COL1A1 were consistent with the results of bioinformatics analysis, and PPP1R12B showed a downward trend that may be related to m6A.

CONCLUSION

The results suggest that piRPGs might serve a significant role in AAA. PPP1R12B, COL1A1, and LRP10 had potential as diagnostic-specific biomarkers for AAA and performed better in XGboost model. The expression and localization of PPP1R12B and COL1A1 were experimentally verified. Besides, downregulation of PPP1R12B caused by m6A might contribute to the formation of AAA.

m6A-related lncRNAs as potential biomarkers and the lncRNA ELFN1-AS1/miR-182-5p/BCL-2 regulatory axis in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoid subtype. However, unsatisfactory survival outcomes remain a major challenge, and the underlying mechanisms are poorly understood. N6-methyladenosine (m6A), the most common internal modification of eukaryotic mRNA, participates in cancer pathogenesis. In this study, m6A-associated long non-coding RNAs (lncRNA) were retrieved from publicly available databases. Univariate, LASSO, and multivariate Cox regression analyses were performed to establish an m6A-associated lncRNA model specific to DLBCL. Kaplan-Meier curves, principal component analysis, functional enrichment analyses and nomographs were used to study the risk model. The underlying clinicopathological characteristics and drug sensitivity predictions against the model were identified. Risk modelling based on the three m6A-associated lncRNAs was an independent prognostic factor. By regrouping patients using our model-based method, we could differentiate patients more accurately for their response to immunotherapy. In addition, prospective compounds that can target DLBCL subtypes have been identified. The m6A-associated lncRNA risk-scoring model developed herein holds implications for DLBCL prognosis and clinical response prediction to immunotherapy. In addition, we used bioinformatic tools to identify and verify the ceRNA of the m6A-associated lncRNA ELFN1-AS1/miR-182-5p/BCL-2 regulatory axis. ELFN1-AS1 was highly expressed in DLBCL and DLBCL cell lines. ELFN1-AS1 inhibition significantly reduced the proliferation of DLBCL cells and promoted apoptosis. ABT-263 inhibits proliferation and promotes apoptosis in DLBCL cells. In vitro and in vivo studies have shown that ABT-263 combined with si-ELFN1-AS1 can inhibit DLBCL progression.

The role of epigenetic and epitranscriptomic modifications in plants exposed to non-essential metals

Contamination of the soil with non-essential metals and metalloids is a serious problem in many regions of the world. These non-essential metals and metalloids are toxic to all organisms impacting crop yields and human health. Crop plants exposed to high concentrations of these metals leads to perturbed mineral homeostasis, decreased photosynthesis efficiency, inhibited cell division, oxidative stress, genotoxic effects and subsequently hampered growth. Plants can activate epigenetic and epitranscriptomic mechanisms to maintain cellular and organism homeostasis. Epigenetic modifications include changes in the patterns of cytosine and adenine DNA base modifications, changes in cellular non-coding RNAs, and remodeling histone variants and covalent histone tail modifications. Some of these epigenetic changes have been shown to be long-lasting and may therefore contribute to stress memory and modulated stress tolerance in the progeny. In the emerging field of epitranscriptomics, defined as chemical, covalent modifications of ribonucleotides in cellular transcripts, epitranscriptomic modifications are postulated as more rapid modulators of gene expression. Although significant progress has been made in understanding the plant's epigenetic changes in response to biotic and abiotic stresses, a comprehensive review of the plant's epigenetic responses to metals is lacking. While the role of epitranscriptomics during plant developmental processes and stress responses are emerging, epitranscriptomic modifications in response to metals has not been reviewed. This article describes the impact of non-essential metals and metalloids (Cd, Pb, Hg, Al and As) on global and site-specific DNA methylation, histone tail modifications and epitranscriptomic modifications in plants.

AL360181.1 promotes proliferation and invasion in colon cancer and is one of ten m6A-related lncRNAs that predict overall survival

Background

N6-methyladenosine (m6A) exerted a pivotal role in colon cancer. Nevertheless, the long non-coding RNAs (lncRNAs) associated with this process have yet to be elucidated.

Methods

The open-access data used for analysis was downloaded from The Cancer Genome Atlas (TCGA) database for analysis, employing the R software for computational evaluations. The RNA level of specific molecules was assessed using the quantitative real-time PCR. CCK8, colony formation and transwell assay were used to evaluate the proliferation, invasion and migration ability of colon cancer cells.

Results

Here, we identified the m6A regulators from TCGA data and subsequently pinpointed lncRNAs with a -Cor- > 0.3 and P < 0.05, categorizing them as m6A-associated lncRNAs. Moreover, we formulated a prognosis signature rooted in ten m6A-related lncRNAs, consisting of AL360181.1, PCAT6, SNHG26, AC016876.1, AC104667.2, AL114730.3, LINC02257, AC147067.1, AP006621.3 and AC009237.14. This signature exhibited notable predictive accuracy in gauging patient survival. Immune-related evaluations revealed varied immune cell infiltration patterns across different risk groups, with our findings suggesting superior immunotherapy response in low-risk patients. Biological enrichment analysis indicated that the high-risk patients had a higher activity of multiple carcinogenic pathways, including glycolysis. The previously unreported lncRNA, AL360181.1, displayed a connection to glycolytic activity and diminished survival rates, warranting further investigation. The result indicated that AL360181.1 was correlated with more aggressive clinical characteristics. Immune infiltration assessments found AL360181.1 to have a positive correlation with Tcm infiltration, but an inverse relationship with entities like Th2 cells, T cells, neutrophils and macrophages. Biological enrichment analysis indicated that the pathways of WNT/β-catenin, pancreas beta cells, hedgehog signaling and some metabolism pathways were upregulated in high AL360181.1 patients. In vitro experiments showed that AL360181.1 was upregulated in the colon cancer cells. Moreover, AL360181.1 significantly promotes the proliferation, invasion and migration of colon cancer cells.

Conclusions

Our results can provide direction for future studies on m6A-related lncRNA in colon cancer.

m6A regulates breast cancer proliferation and migration through stage-dependent changes in Epithelial to Mesenchymal Transition gene expression

While many factors have been implicated in breast cancer progression, effective treatments are still lacking. In recent years, it has become clear that posttranscriptional regulation plays a key role in the aberrant gene expression underlying malignancy and metastasis. For example, the mRNA modification N6-methyladenosine (m6A) is involved in numerous post-transcriptional regulation processes and has been implicated in many cancer types, including breast cancer. Despite intense study, even within a single type of cancer, there is little consensus, and often conflicting results, as to the role of m6A, suggesting other factors must influence the process. The goal of this study was to determine if the effects of m6A manipulation on proliferation and migration differed based on the stage of disease progression. Using the MCF10 model of breast cancer, we reduced m6A levels by targeting METTL3, the main cellular m6A RNA methyltransferase. Knocking down Mettl3 at different stages of breast cancer progression indeed shows unique effects at each stage. The early-stage breast cancer line showed a more proliferative phenotype with the knockdown of Mettl3 while the transformed breast cancer line showed a more migratory phenotype. Interestingly, the metastasized breast cancer cell line showed almost no effect on phenotype with the knockdown of Mettl3. Furthermore, transcriptome wide analysis revealed EMT as the probable pathway influencing the phenotypic changes. The results of this study may begin to address the controversy of m6A's role in cancer and suggest that m6A may have a dynamic role in cancer that depends on the stage of progression.

N6-methyladenosine methylation in kidney injury

Multiple mechanisms are involved in kidney damage, among which the role of epigenetic modifications in the occurrence and development of kidney diseases is constantly being revealed. However, N6-methyladenosine (M6A), a well-known post-transcriptional modification, has been regarded as the most prevalent epigenetic modifications in higher eukaryotic, which is involved in various biological processes of cells such as maintaining the stability of mRNA. The role of M6A modification in the mechanism of kidney damage has attracted widespread attention. In this review, we mainly summarize the role of M6A modification in the progression of kidney diseases from the following aspects: the regulatory pattern of N6-methyladenosine, the critical roles of N6-methyladenosine in chronic kidney disease, acute kidney injury and renal cell carcinoma, and then reveal its potential significance in the diagnosis and treatment of various kidney diseases. A better understanding of this field will be helpful for future research and clinical treatment of kidney diseases.

Multi-omics analysis of MRPL-13 as a tumor-promoting marker from pan-cancer to lung adenocarcinoma

BACKGROUND

As a member of the mitochondrial ribosomal protein family, mitochondrial ribosomal protein L13 (MRPL13) is responsible for synthesizing mitochondrial proteins in cells. Several studies have indicated that MRPL13 is associated with the proliferation cycle, migration ability, apoptosis and autophagy of cancer cells. However, a thorough examination of MRPL13 across cancers remains uncertain. Therefore, we tried to clarify the relationship between MRPL13 and pan-cancer, and verified it in lung adenocarcinoma by various methods. Finally, our research is expected to reveal new targets for pan-cancer treatment and improve the prognosis of cancer patients.

METHODS

Using bioinformatics tools, we quantified the differential expression of MRPL13 between cancer tissues and corresponding or noncorresponding normal tissues across cancers. We also analyzed the relationships between MRPL13 expression levels and several factors, including diagnosis, prognosis, mutation, functional signaling pathways, immune infiltration, RNA modification, and the relationship with cuproptosis-related genes. Furthermore, we studied the relationship between the expression level of MRPL13 across cancers and the change in cancer functional status through single-cell data. In addition, quantitative experiments (PCR and Western blot) proved that the expression of MRPL13 was significantly different between LUAD and control samples. Finally, the effect of knocking out MRPL13 on cancer cells was compared by gene silencing experiments. In summary, we used a combination of bioinformatics and experimental applications to study the potential roles of MRPL13 in cancer.

RESULTS

After conducting a multidimensional analysis, we found that the application of MRPL13 multigroup analysis can effectively improve the diagnostic efficiency of various cancers and predict the prognosis of cancer. Moreover, MRPL13 in pan-cancer is related to the cancer immune infiltration pattern, methylation level and cuproptosis-related genes. Furthermore, single-cell data analysis showed that the modules of metastasis, EMT, cell cycle, DNA repair, invasion, DNA damage and proliferation were positively correlated with the expression of MRPL13 in LUAD (Lung adenocarcinoma), while the modules of hypoxia and inflammation were negatively correlated. Moreover, through quantitative experiments, we observed higher expression of MRPL13 in cancer tissues at the RNA or protein level. Knockdown of MRPL13 in LUAD led to decreased cancer cell survival, delayed tumor division and migration, reduced invasion, and increased cancer cell apoptosis.

CONCLUSIONS

Our study demonstrates the potential of using MRPL13 as a molecular biomarker for diagnosing and suggesting the prognosis of certain malignant tumors. Furthermore, our research shows that MRPL13 may be an effective therapeutic target for lung adenocarcinoma.

Characteristics of m6A-related LncRNAs in breast cancer as prognostic biomarkers and immunotherapy

N6-methyladenosine (m6A) is a common RNA modification in coding and non-coding RNAs and plays an important role in the occurrence and development of breast cancer (BC). However, the role of m6A-related lncRNAs in breast cancer prognosis is unclear. This study aimed to help verify the biological function of m6A-related lncRNAs in breast cancer prognosis through bio-informatics techniques. First, we screened 18 m6A-related lncRNAs from the TCGA database: AL137847.1, AC137932.2, OTUD6B-AS1, MORF4L2-AS1, AC078846.1, AC012442.1, AL118556.1, AL138955.1, AC009754.1, AC024257.4, AL391095.1, AC024270.3, AC087392.1, LINC02649, AC090948.2, AL158212.1, ITGA6-AS1, AL133243.2 and constructed a risk-prognosis model based on this. Based on the model's median risk score, BC patients were divided into high-risk and low-risk groups. Then, the predictive value of the model was verified by Cox regression, Lasso regression, Kaplan-Meier curve and ROC curve analysis, and biological differences between the two groups were verified by GO enrichment analysis, tumor mutation burden, immune indications and in vitro tests. Importantly, the risk score of this prognostic model is an excellent independent prognostic factor, and m6A regulators are differentially expressed in patients with different risks. In addition, based on patients' different sensitivities to drugs, some drug candidates for different risk populations are screened to provide targets for breast cancer treatment. The difference in immune function between high-risk and low-risk patients also affected the sensitivity to immunotherapy. In the validation of clinical samples, we analyzed the expression of relevant lncRNAs in different risk groups and speculated the possible impact on the prognosis of breast cancer patients. The risk assessment tool built based on the full analysis of these m6A-related genes and m6A-related lncRNA libraries, as well as the m6A-related lncRNAs, has a high prognostic prediction ability, which may provide a supplementary screening method for accurately judging the prognosis of BC and a new perspective for personalized treatment of breast cancer patients.

M6A transcriptome-wide map of circRNAs identified in the testis of normal and AZ-treated Xenopus laevis

BACKGROUND

Evidence showed that N6-methyladenosine (m6A) is strongly associated with male germline development. However, the role of m6A methylation on circRNAs in amphibians remains unknown. In this study, we conducted m6A sequencing analysis to explore the m6A transcriptome-wide profile of circRNAs in testis tissues of Xenopus laevis (X. laevis) with and without treatment with 100 µg/L atrazine (AZ).

RESULTS

The analysis showed that m6A modification of circRNAs enriched in sense overlapping in testes of X. laevis. We identified the differential m6A modification sites within circRNAs in testes of AZ-exposed X. laevis and compared that with animals from control group. The results showed that a total of 1507 methylated m6A sites was induced by AZ (760 up-methylated and 747 down-methylated). The cross-analysis exhibited a negative correlation of differentially methylated m6A peaks and circRNAs expression level. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that 20 key pathways may be involved in the mechanism of testis damage of AZ-exposed X. laevis.

CONCLUSIONS

These findings indicated that differentially m6A-methylated circRNAs may play important roles in abnormal testis development of AZ-exposed X. laevis. This study is the first report about a map of m6A modification of circRNAs in male X. laevis and provides a basis for further studying on the function and mechanism of m6A methylation of circRNAs in the testis development of amphibian.

N6-methyladenosine modification: an important player in the tumor immune microenvironment

The decoration of RNA with N6-methyladenosine (m6A) is a reversible post-transcriptional modification that plays an important regulatory role in all eukaryotic life activities. The m6A modification of RNA regulates the development and progression of tumors, including bladder cancer, melanoma, Lewis lung carcinoma, and hepatocellular carcinoma. The tumor immune microenvironment (TIME) includes immune cells, cytokines, and cell surface molecules, which interact with each other and ultimately determine the flow of tumor immunity. The onset of cancer implies that the TIME has been reshaped into a pro-tumor state. The key to cancer treatment lies in reshaping the TIME to reset the anti-tumor immune response. Here, we have reviewed how RNA m6A modification affects the TIME, and discussed the merits of using m6A regulator inhibitors as an individual treatment strategy as well as in combination with immune checkpoint blockade therapy.

Prognostic and tumor microenvironmental feature of clear cell renal cell carcinoma revealed by m6A and lactylation modification-related genes

Background

Both lactylation and m6A modification have important implications for the development of clear cell renal cell carcinoma (ccRCC), and we aimed to use crosstalk genes of both to reveal the prognostic and immunological features of ccRCC.

Methods

Our first step was to look for lactylation-related genes that differed between normal and tumor tissues, and then by correlation analysis, we found the genes associated with M6A. Following that, ccRCC subtypes will be identified and risk models will be constructed to compare the prognosis and tumor microenvironment among different subgroups. A nomogram was constructed to predict the prognosis of ccRCC, and in vitro, experiments were conducted to validate the expression and function of key genes.

Results

We screened 100 crosstalk genes and identified 2 ccRCC subtypes. A total of 11 prognostic genes were screened for building a risk model. we observed higher immune scores, elevated tumor mutational burden, and microsatellite instability scores in the high-risk group. Therefore, individuals classified as high-risk would derive greater benefits from immunotherapy. The nomogram's ability to predict overall survival with a 1-year AUC of 0.863 demonstrates its significant practical utility. In addition, HIBCH was identified as a potential therapeutic target and its expression and function were verified by in vitro experiments.

Conclusion

In addition to developing a precise prognostic nomogram for patients with ccRCC, our study also discovered the potential of HIBCH as a biomarker for the disease.

Expression, prognostic value and mechanism of SP100 family in pancreatic adenocarcinoma

BACKGROUND

Pancreatic adenocarcinoma (PAAD) is one of the most aggressive malignancies with a very poor prognosis. Exploring more therapeutic targets and prognostic biomarkers is of great significance to improve the prognosis of PAAD patients. Increasing evidence supports that the speckled protein (SP) 100 family is associated with human cancer and immune disorders. However, the function of the SP100 family members in PAAD is still unclear.

METHODS

R, Cytoscape, cBioPortal, and other software and online databases were used to comprehensively analyze the expression characteristics, prognostic value, and oncogenic mechanism of the SP100 family in PAAD.

RESULTS

The high expression of SP100 family members in PAAD was significantly correlated with poor clinicopathological features and poor prognosis of PAAD patients. Mechanistically, TP53 mutations were significantly associated with the expression levels of the SP100 family members, which were significantly coexpressed with M6A methylation regulators and were activated in multiple oncogenic pathways, including the EMT pathways. Moreover, we found that their expression levels were significantly correlated with the sensitivity of multiple traditional chemotherapeutic drugs.

CONCLUSION

The SP100 family is closely related to the occurrence and development of PAAD and can be used as a new biomarker and therapeutic target for patients with PAAD.

Comprehensive analysis of differences in N6-methyladenosine RNA methylomes in Helicobacter pylori infection

Background: Helicobacter pylori (H.pylori) infection is an important factor in the occurrence of human gastric diseases, but its pathogenic mechanism is not clear. N6-methyladenosine (m6A) is the most prevalent reversible methylation modification in mammalian RNA and it plays a crucial role in controlling many biological processes. However, there are no studies reported that whether H. pylori infection impacts the m6A methylation of stomach. In this study, we measured the overall level changes of m6A methylation of RNA under H. pylori infection through in vitro and in vivo experiment. Methods: The total quantity of m6A was quantified in gastric tissues of clinical patients and C57 mice with H. pylori infection, as well as acute infection model [H. pylori and GES-1 cells were cocultured for 48 h at a multiplicity of infection (MOI) from of 10:1 to 50:1]. Furthermore, we performed m6A methylation sequencing and RNA-sequencing on the cell model and RNA-sequencing on animal model. Results: Quantitative detection of RNA methylation showed that H. pylori infection group had higher m6A modification level. M6A methylation sequencing identified 2,107 significantly changed m6A methylation peaks, including 1,565 upregulated peaks and 542 downregulated peaks. A total of 2,487 mRNA was upregulated and 1,029 mRNA was downregulated. According to the comprehensive analysis of MeRIP-seq and RNA-seq, we identified 200 hypermethylation and upregulation, 129 hypermethylation but downregulation, 19 hypomethylation and downregulation and 106 hypomethylation but upregulation genes. The GO and KEGG pathway analysis of these differential methylation and regulatory genes revealed a wide range of biological functions. Moreover, combining with mice RNA-seq results, qRT- PCR showed that m6A regulators, METTL3, WTAP, FTO and ALKBH5, has significant difference; Two key genes, PTPN14 and ADAMTS1, had significant difference by qRT- PCR. Conclusion: These findings provide a basis for further investigation of the role of m6A methylation modification in H. pylori-associated gastritis.

Identification of m6a-related signature genes in esophageal squamous cell carcinoma by machine learning method

Background: We aimed to construct and validate the esophageal squamous cell carcinoma (ESCC)-related m6A regulators by means of machine leaning. Methods: We used ESCC RNA-seq data of 66 pairs of ESCC from West China Hospital of Sichuan University and the transcriptome data extracted from The Cancer Genome Atlas (TCGA)-ESCA database to find out the ESCC-related m6A regulators, during which, two machine learning approaches: RF (Random Forest) and SVM (Support Vector Machine) were employed to construct the model of ESCC-related m6A regulators. Calibration curves, clinical decision curves, and clinical impact curves (CIC) were used to evaluate the predictive ability and best-effort ability of the model. Finally, western blot and immunohistochemistry staining were used to assess the expression of prognostic ESCC-related m6A regulators. Results: 2 m6A regulators (YTHDF1 and HNRNPC) were found to be significantly increased in ESCC tissues after screening out through RF machine learning methods from our RNA-seq data and TCGA-ESCA database, respectively, and overlapping the results of the two clusters. A prognostic signature, consisting of YTHDF1 and HNRNPC, was constructed based on our RNA-seq data and validated on TCGA-ESCA database, which can serve as an independent prognostic predictor. Experimental validation including the western and immunohistochemistry staining were further successfully confirmed the results of bioinformatics analysis. Conclusion: We constructed prognostic ESCC-related m6A regulators and validated the model in clinical ESCC cohort as well as in ESCC tissues, which provides reasonable evidence and valuable resources for prognostic stratification and the study of potential targets for ESCC.

Identification and experimental validation of key m6A modification regulators as potential biomarkers of osteoporosis

Osteoporosis (OP) is a severe systemic bone metabolic disease that occurs worldwide. During the coronavirus pandemic, prioritization of urgent services and delay of elective care attenuated routine screening and monitoring of OP patients. There is an urgent need for novel and effective screening diagnostic biomarkers that require minimal technical and time investments. Several studies have indicated that N6-methyladenosine (m6A) regulators play essential roles in metabolic diseases, including OP. The aim of this study was to identify key m6A regulators as biomarkers of OP through gene expression data analysis and experimental verification. GSE56815 dataset was served as the training dataset for 40 women with high bone mineral density (BMD) and 40 women with low BMD. The expression levels of 14 major m6A regulators were analyzed to screen for differentially expressed m6A regulators in the two groups. The impact of m6A modification on bone metabolism microenvironment characteristics was explored, including osteoblast-related and osteoclast-related gene sets. Most m6A regulators and bone metabolism-related gene sets were dysregulated in the low-BMD samples, and their relationship was also tightly linked. In addition, consensus cluster analysis was performed, and two distinct m6A modification patterns were identified in the low-BMD samples. Subsequently, by univariate and multivariate logistic regression analyses, we identified four key m6A regulators, namely, METTL16, CBLL1, FTO, and YTHDF2. We built a diagnostic model based on the four m6A regulators. CBLL1 and YTHDF2 were protective factors, whereas METTL16 and FTO were risk factors, and the ROC curve and test dataset validated that this model had moderate accuracy in distinguishing high- and low-BMD samples. Furthermore, a regulatory network was constructed of the four hub m6A regulators and 26 m6A target bone metabolism-related genes, which enhanced our understanding of the regulatory mechanisms of m6A modification in OP. Finally, the expression of the four key m6A regulators was validated in vivo and in vitro, which is consistent with the bioinformatic analysis results. Our findings identified four key m6A regulators that are essential for bone metabolism and have specific diagnostic value in OP. These modules could be used as biomarkers of OP in the future.

Identification of an m6A-related lncRNA prognostic signature in colorectal cancer

Data sets of colorectal cancer (CRC) were obtained from The Cancer Genome Atlas (TCGA), three N6-methyladenosine (m6A) subtypes were identified using 21 m6A-related long noncoding RNAs (lncRNAs) and differential m6A subtypes of different CRC tumors were determined in this study to evaluate the m6A expression and the prognosis of patients with CRC. Subsequently, eight key lncRNAs were identified based on co-expression with 21 m6A-related genes in CRC tumors using the single-factor Cox and least absolute shrinkage and selection operator. Finally, an m6A-related lncRNA risk score model of CRC tumor was established using multifactor Cox regression based on the eight important lncRNAs and found to have a better performance in evaluating the prognosis of patients in the TCGA-CRC data set. TCGA-CRC tumor samples were divided based on the risk scores: high and low. Then, the clinical characteristics, tumor mutation load, and tumor immune cell infiltration difference between the high- and low-risk-score groups were explored, and the predictive ability of the risk score was assessed for immunotherapeutic benefits. We found that the risk score model can determine the overall survival, be a relatively independent prognostic indicator, and better evaluate the immunotherapeutic benefits for patients with CRC. This study provides data support for accurate immunotherapy in CRC.

Weighted gene co-expression network analysis and CIBERSORT screening of key genes related to m6A methylation in Hirschsprung's disease

Hirschsprung's disease (HSCR) is a neural crest disease that results from the failure of enteric neural crest cells (ENCCs) to migrate to the corresponding intestinal segment. The RET gene, which regulates enteric neural crest cell proliferation and migration, is considered one of the main risk factors for HSCR and is commonly used to construct HSCR mouse models. The epigenetic mechanism of m6A modification is involved in HSCR. In this study, we analyzed the GEO database (GSE103070) for differentially expressed genes (DEGs) and focused on m6A-related genes. Comparing the RNA-seq data of Wide Type and RET Null, a total of 326 DEGs were identified, of which 245 genes were associated with m6A. According to the CIBERSORT analysis, the proportion of Memory B-cell in RET Null was significantly higher than that of Wide Type. Venn diagram analysis was used to identify key genes in the selected memory B-cell modules and DEGs associated with m6A. Enrichment analysis showed that seven genes were mainly involved in focal adhesion, HIV infection, actin cytoskeleton organization and regulation of binding. These findings could provide a theoretical basis for molecular mechanism studies of HSCR.

Identification of disease-related genes and construction of a gene co-expression database in non-alcoholic fatty liver disease

Background: The mechanism of NAFLD progression remains incompletely understood. Current gene-centric analysis methods lack reproducibility in transcriptomic studies. Methods: A compendium of NAFLD tissue transcriptome datasets was analyzed. Gene co-expression modules were identified in the RNA-seq dataset GSE135251. Module genes were analyzed in the R gProfiler package for functional annotation. Module stability was assessed by sampling. Module reproducibility was analyzed by the ModulePreservation function in the WGCNA package. Analysis of variance (ANOVA) and Student's t-test was used to identify differential modules. The receiver operating characteristic (ROC) curve was used to illustrate the classification performance of modules. Connectivity Map was used to mine potential drugs for NAFLD treatment. Results: Sixteen gene co-expression modules were identified in NAFLD. These modules were associated with multiple functions such as nucleus, translation, transcription factors, vesicle, immune response, mitochondrion, collagen, and sterol biosynthesis. These modules were stable and reproducible in the other 10 datasets. Two modules were positively associated with steatosis and fibrosis and were differentially expressed between non-alcoholic steatohepatitis (NASH) and non-alcoholic fatty liver (NAFL). Three modules can efficiently separate control and NAFL. Four modules can separate NAFL and NASH. Two endoplasmic reticulum related modules were both upregulated in NAFL and NASH compared to normal control. Proportions of fibroblasts and M1 macrophages are positively correlated with fibrosis. Two hub genes Aebp1 and Fdft1 may play important roles in fibrosis and steatosis. m6A genes were strongly correlated with the expression of modules. Eight candidate drugs for NAFLD treatment were proposed. Finally, an easy-to-use NAFLD gene co-expression database was developed (available at https://nafld.shinyapps.io/shiny/). Conclusion: Two gene modules show good performance in stratifying NAFLD patients. The modules and hub genes may provide targets for disease treatment.

N6-methyladenosine mRNA methylation is important for the light response in soybean

N6-methyladenosine (m⁶A) modification of messenger RNA (mRNA) is the most prevalent and abundant modification in eukaryotic mRNA and posttranscriptionally modulates the transcriptome at almost all stages of mRNA metabolism. In plants, m⁶A is crucial for embryonic-phase growth, flowering time control, microspore generation and fruit maturation. However, the role of m⁶A in plant responses to light, the most important environmental stimulus, remains unexplored. Here, we profile the m⁶A transcriptome of Williams 82, a soybean cultivar, and reveal that m⁶A is highly conserved and plays an important role in the response to light stimuli in soybean. Similar to the case in Arabidopsis, m⁶A in soybean is enriched not only around the stop codon and within the 3'UTR but also around the start codon. Moreover, genes with methylation occurring in the 3'UTR have higher expression levels and are more prone to alternative splicing. The core genes in the light signaling pathway, GmSPA1a, GmPRR5e and GmBIC2b, undergo changes in methylation modification and transcription levels in response to light. KEGG pathway analysis revealed that differentially expressed genes with differential m⁶A peaks were involved in the "photosynthesis" and "circadian rhythm" pathways. Our results highlight the important role played by epitranscriptomic mRNA methylation in the light response in soybean and provide a solid basis for determining the functional role of light on RNA m⁶A modification in this plant.

Constructing and validating of m6a-related genes prognostic signature for stomach adenocarcinoma and immune infiltration: Potential biomarkers for predicting the overall survival

Background

Stomach adenocarcinoma (STAD) arises from the mutations of stomach cells and has poor overall survival. Chemotherapy is commonly indicated for patients with stomach cancer following surgical resection. The most prevalent alteration that affects cancer growth is N6-methyladenosine methylation (m6A), although the possible function of m6A in STAD prognosis is not recognized.

Method

The research measured predictive FRGs in BLCA samples from the TCGA and GEO datasets. Data on the stemness indices (mRNAsi), gene mutations, copy number variations (CNV), tumor mutation burden (TMB), and corresponding clinical characteristics were obtained from TCGA and GEO. STAD from TCGA and GEO at 24 m6A was investigated. Lasso regression was used to construct the prediction model to assess the m6A prognostic signals in STAD. In addition, the correlation between m6a and immune infiltration in STAD patients was discussed using GSVA and ssGSEA analysis. Based on these genes, GO and KEGG analyses were performed to identify key biological functions and key pathways.

Result

A significant relationship was discovered between numerous m6A clusters and the tumor immune microenvironment, as well as three m6A alteration patterns with different clinical outcomes. Furthermore, GSVA and ssGSEA showed that m6A clusters were significantly associated with immune infiltration in the STAD. The low-m6Ascore group had a lower immunotherapeutic response than the high-m6Ascore group. ICIs therapy was more effective in the group with a higher m6Ascore. Three writers (VIRMA, ZC3H13, and METTL3) showed significantly lower expression, whereas five authors (METTL14, METTL16, WTAP, RBM15, and RBM15B) showed considerably higher expression. Three readers (YTHDC2, YTHDF2, and LRPPRC) had higher levels of expression, whereas eleven readers (YTHDC1, YTHDF1, YTHDF3, HNRNPC, FMR1, HNRNPA2B1, IGFBP1, IGFBP2, IGFBP3, and RBMX) had lower levels. As can be observed, the various types of m6 encoders have varied ramifications for STAD control.

Conclusion

STAD occurrence and progression are linked to m6A-genes. Corresponding prognostic models help forecast the prognosis of STAD patients. m6A-genes and associated immune cell infiltration in the tumor microenvironment (TME) may serve as potential therapeutic targets in STAD, which requires further trials. In addition, the m6a-related gene signature offers a viable alternative to predict bladder cancer, and these m6A-genes show a prospective research area for STAD targeted treatment in the future.

Prognostic value, immune signature and molecular mechanisms of the SUMO family in pancreatic adenocarcinoma

Background: Pancreatic adenocarcinoma (PAAD) has a high degree of malignancy and a very poor prognosis, and the 5-year overall survival rate of patients is approximately 7%. To improve the prognosis of patients with PAAD, a more comprehensive and in-depth study of the pathogenesis of PAAD and the identification of new diagnostic markers and treatment targets are urgently needed. Increasing evidence supports that the small ubiquitin-like modifier (SUMO) family is closely related to the occurrence and development of a variety of cancers. However, the function of the SUMO family in PAAD is not clear, and related research is very scarce. Methods: R, Cytoscape, cBioPortal, and other software and online databases were used to comprehensively analyze the expression characteristics, prognostic value, and oncogenic mechanism of the SUMO family in PAAD. Results: SUMO family members are highly expressed in PAAD, and high expression of SUMO family members is significantly associated with poor clinicopathological features and poor prognosis in PAAD patients. In addition, SUMO family members are significantly coexpressed with M6A methylation regulators and various oncogenes and play an activating role in various oncogenic pathways, including EMT. Furthermore, it is worth noting that the close association between SUMO family members and TP53 mutation status and the negative regulatory effect of SUMO1/2 on PAAD immunity may represent the potential mechanism by which SUMO family members promote the development of PAAD. Moreover, the coexpression characteristics of SUMO family members and a variety of cancer-promoting immune checkpoint genes, as well as the positive correlation between SUMO4 expression level and the sensitivity of various targeted or chemotherapeutic drugs, including gemcitabine, paclitaxel, and doxorubicin, suggest future clinical directions of this study. Conclusion: The SUMO family is closely related to the occurrence and development of PAAD and can be used as a new biomarker and therapeutic target for patients with PAAD.

Analysis of M6A associated lncRNAs in prognosis and immune response of NSCLC patients

Distinguishing between N6-methyladenosine (m6A)-associated long noncoding RNAs (lncRNAs) is crucial in non-small-cell lung cancer (NSCLC) patients. In this research, the prognosis and immunotherapeutic response of lncRNAs and m6A in NSCLC were examined. lncRNAs related to m6A were identified using co-expression analyses, and their prognostic impact on patients with NSCLC was assessed using univariate Cox regression analysis. Sixty-three m6A-associated lncRNAs were determined as prognostic lncRNAs, and on this basis, 25 m6A-associated lncRNAs were screened by least absolute shrinkage and selection operator (lasso) Cox regression. Multivariable Cox analysis obtained 14 m6A-associated lncRNAs for the construction of risk model. The NSCLC patients were grouped into different risk subgroups in accordance with the median of the risk fraction in each data, and we evaluated the differences of potential immunotherapeutic characteristics and drug sensitivity prediction between the two subgroups. By using this model to recombine patients, they can be effectively distinguished in terms of the immunotherapy response. Furthermore, candidate compounds for the differentiation of NSCLC subtypes were identified. The model based on 14 m6A-associated lncRNAs is a promising prognostic biomarker, which may help to predict the efficacy of immunotherapy in NSCLC patients and provide a theoretical basis for improving the outcome of patients.

Corrigendum: m6A regulatory gene-mediated methylation modification in glioma survival prediction

[This corrects the article DOI: 10.3389/fgene.2022.873764.].

Expression patterns of platinum resistance-related genes in lung adenocarcinoma and related clinical value models

The purpose of this study was to explore platinum resistance-related biomarkers and mechanisms in lung adenocarcinoma. Through the analysis of gene expression data of lung adenocarcinoma patients and normal patients from The Cancer Genome Atlas, Gene Expression Omnibus database, and A database of genes related to platinum resistance, platinum resistance genes in lung adenocarcinoma and platinum resistance-related differentially expressed genes were obtained. After screening by a statistical significance threshold, a total of 252 genes were defined as platinum resistance genes with significant differential expression, of which 161 were up-regulated and 91 were down-regulated. The enrichment results of up-regulated gene Gene Ontology (GO) showed that TOP3 entries related to biological processes (BP) were double-strand break repair, DNA recombination, DNA replication, the down-regulated gene GO enriches the TOP3 items about biological processes (BP) as a response to lipopolysaccharide, muscle cell proliferation, response to molecule of bacterial origin. Gene Set Enrichment Analysis showed that the top three were e2f targets, g2m checkpoint, and rgf beta signaling. A prognostic model based on non-negative matrix factorization classification showed the characteristics of high- and low-risk groups. The prognostic model established by least absolute shrinkage and selection operator regression and risk factor analysis showed that genes such as HOXB7, NT5E, and KRT18 were positively correlated with risk score. By analyzing the differences in m6A regulatory factors between high- and low-risk groups, it was found that FTO, GPM6A, METTL3, and YTHDC2 were higher in the low-risk group, while HNRNPA2B1, HNRNPC, TGF2BP1, IGF2BP2, IGF2BP3, and RBM15B were higher in the high-risk group. Immune infiltration and drug sensitivity analysis also showed the gene characteristics of the platinum-resistant population in lung adenocarcinoma. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 were lower in the tumor expression group, and that the survival of the low expression group was worse than that of the high expression group. In conclusion, the results of this study show that platinum resistance-related differentially expressed genes in lung adenocarcinoma are mainly concentrated in biological processes such as DNA recombination and response to lipopolysaccharide. The validation set proved that the high-risk group of our prognostic model had poor survival. M6A regulatory factor analysis, immune infiltration, and drug sensitivity analysis all showed differences between high and low-risk groups. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 could be protective factors. Further exploration of the potential impact of these genes on the risk and prognosis of drug-resistant patients with lung adenocarcinoma would provide theoretical support for future research.

LDHA is a prognostic biomarker on the immune response in pancreatic adenocarcinoma and associated with m6A modification

PURPOSE

N6-methyladenosine (m⁶A) is tightly associated with the progression of pancreatic ductal adenocarcinoma (PDAC). Another prominent feature of PDAC is metabolic reprogramming, which provides sufficient nutrients to support rapid cell growth via the tumor microenvironment. However, the underlying influences of m⁶A-associated metabolic genes on the PDAC microenvironment remain poorly understood. Therefore, we sought to construct a survival prediction model using m⁶A-related genes to clarify the molecular characteristics of PDAC.

METHODS

In the present study, m⁶A-related metabolic genes were obtained from The Cancer Genome Atlas (TCGA) pancreatic adenocarcinoma dataset and subjected to coexpression analysis. Consensus clustering recognized two distinct subgroups with different immune cell infiltration patterns according to the expression of m⁶A-associated metabolic genes. Multivariate Cox regression analyses and least absolute shrinkage and selection operator (LASSO) analysis were adopted to create an m⁶A-related metabolism model. A nomogram including clinical features and the risk score based on the expression of m⁶A-related metabolism regulators was constructed.

RESULTS

A four-gene signature comprising ATP8B2, GMPS, LDHA and SDR39U1 was built to predict the overall survival (OS) of PDAC patients. This signature also robustly predicted survival in two independent validation cohorts from the International Cancer Genome Consortium (ICGC) and ArrayExpress (E-MTAB-6134). The four-gene signature divided patients into high- and low-risk groups with distinct OS values as verified by the log-rank test. Among the four genes, LDHA was upregulated in both PDAC tissues and cell lines.

CONCLUSIONS

Collectively, we analyzed the immune microenvironment, predicted drug sensitivity and assessed the implications of the mutation landscape based on the crosstalk between m⁶A regulators and metabolic rewiring, and we also constructed a novel signature based on m⁶A-associated metabolic genes to predict PDAC prognosis.

Comprehensive analysis of N6-methyladenosine-related lncRNAs reveals distinct hepatocellular carcinoma subtypes with immunotherapeutic implications

Accumulating studies have demonstrated critical roles of N6-methyladenosine (m6A) modification and long noncoding RNAs (lncRNAs) in the biological processes leading to occurrence, development and chemoresistance of cancers. However, the specific identities and functional roles of lncRNAs associated with m6A modification in hepatocellular carcinoma (HCC) remain elusive. In this study, eighty-two prognostic m6A-related lncRNAs (m6A-LncRNAs) were identified in HCC datasets. Patients with HCC were classified into three subtypes (C1, C2 and C3) based on the expression of the m6A-LncRNAs. The three subtypes showed significant differences in clinical features, immune and stromal infiltration signatures, and immunotherapy sensitivity. Subclass C1 was notable for high immune and stromal cell infiltration and active immune responses, low serum α-fetoprotein (AFP) levels and high sensitivity to immune checkpoint inhibitors (ICIs). Subclass C2 showed high metabolic activities and absence of immune infiltration with favorable prognosis. Subclass C3 was associated with an exhausted immune environment, high serum AFP and poor prognosis. Notably, subclass C3 displayed high expression of immune checkpoints but failed to respond to ICIs. Finally, 12 m6A-LncRNA signatures were identified for HCC classification and validated in an external dataset. This integrated analysis indicated that the interactions between m6A methylation and lncRNAs are involved in immune and stromal cell infiltration in HCC, and may provide novel insights into precision diagnostics as well as therapeutics for HCC patients.

Prognostic value of comprehensive typing based on m6A and gene cluster in TNBC

BACKGROUND

Triple-negative breast cancer (TNBC) is resistant to targeted therapy with HER2 monoclonal antibodies and endocrine therapy, because it lacks the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2). TNBC is a subtype of breast cancer with the worst prognosis and the highest mortality rate compared to other subtypes. N6-methyladenosine (m6A) modification is significant in cancer and metastasis, because it can alter gene expression and function at numerous levels, such as RNA splicing, stability, translocation, and translation. There are limited investigations into the connection between TNBC and m6A.

MATERIALS AND METHODS

Breast cancer-related data were retrieved from the Cancer Genome Atlas (TCGA) database, and 116 triple-negative breast cancer cases were identified from the data. The GSE31519 data set, which included 68 cases of TNBC, was obtained from the Gene Expression Omnibus (GEO) database. Survival analysis was used to determine the prognosis of distinct m6A types based on their m6A group, gene group, and m6A score. To investigate the potential mechanism, GO and KEGG analyses were performed on the differentially expressed genes.

RESULTS

The expression of m6A-related genes and their impact on prognosis in TNBC patients were studied. According to the findings, m6A was crucial in determining the prognosis of TNBC patients, and the major m6A-linked genes in this process were YTHDF2, RBM15B, IGFBP3, and WTAP. YTHDF2, RBM15B and IGFBP3 are associated with poor prognosis, while WTAP is associated with good prognosis. By cluster analysis, the gene cluster and the m6A cluster were beneficial in predicting the prognosis of TNBC patients. The m6A score based on m6A and gene clusters was more effective in predicting the prognosis of TNBC patients. Furthermore, the tumor microenvironment may play an important role in the process of m6A, influencing TNBC prognosis.

CONCLUSIONS

N6-adenylic acid methylation (m6A) was important in altering the prognosis of TNBC patients, and the key m6A-associated genes in this process were YTHDF2, RBM15B, IGFBP3, and WTAP. Furthermore, the comprehensive typing based on m6A and gene clusters was useful in predicting TNBC patients' prognosis, showing potential as valuable evaluating tools for TNBC.

ALKBH5/MAP3K8 axis regulates PD-L1+ macrophage infiltration and promotes hepatocellular carcinoma progression

Hepatocellular carcinoma is one of the most common malignant tumors.M6A is a novel epigenetic modification that have been emerged as vital regulators for the progression of HCC. However, the regulatory role, clinical significance and the details of the modification, such as the impact on the local tumor environment, remain largely unclear. Our study showed that ALKBH5 was highly expressed in HCC and high ALKBH5 expression predicted a worse prognosis of HCC patients. Prediction of ALKBH5 function by tissue samples and single cell sequencing Gene Set Variation Analysis. Primary CD3 + T lymphocytes and bone marrow-derived macrophages were used to evaluate the effect of ALKBH5 on immune microenvironment. The results indicated that ALKBH5 promote HCC cell proliferation, metastasis and PD-L1+macrophage recruitment. Mechanistically the results showed that ALKBH5 regulates MAP3K8 expression in a m6A dependent manner which mediates the proliferation and metastasis of HCC cells. ALKBH5 also promotes the activation of JNK and ERK pathways through upregulating MAP3K8, thus regulating the expression of IL-8 and promoting macrophage recruitment. Taken together, these data show that ALKBH5 promotes HCC growth, metastasis and macrophage recruitment through ALKBH5/MAP3K8 axis and it may serve as a potential diagnostic marker and target for treatment of HCC patients.

A Novel YTHDF3-Based Model to Predict Prognosis and Therapeutic Response in Breast Cancer

Background: Due to high tumor heterogeneity, breast cancer (BC) patients still suffer poor survival outcomes. YTHDF3 plays a critical role in the prognosis of BC patients. Hence, we aimed to construct a YTHDF3-based model for the prediction of the overall survival (OS) and the sensitivity of therapeutic agents in BC patients.

Methods: Based on The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/) database, we obtained BC patients’ data (n = 999) with YTHDF3 expression profiles. The association between YTHDF3 expression and 5-year OS was determined via Cox proportional hazards regression (CPHR) analysis. By integrating the variables, we established a prognostic nomogram. The model was estimated via discrimination, calibration ability, and decision curve analysis (DCA). The performance of the model was compared with the TNM stage system through receiver operating characteristic (ROC) curves and DCA. By means of the Genomics of Drug Sensitivity in Cancer (GDSC) database (https://www.cancerrxgene.org/), the therapeutic agents’ response was estimated. Gene set enrichment analysis (GSEA) demonstrated possible biological mechanisms related to YTHDF3. TIMER and CIBERSORTx were employed to analyze the association between YTHDF3 and tumor-infiltrating immune cells.

Results: The high YTHDF3 expression was significantly correlated with poor 5-year OS in BC patients. Through multivariate CPHR, four independent prognostic variables (age, TNM stage, YTHDF3 expression, and molecular subtype) were determined. On the basis of the four factors, a YTHDF3-based nomogram was built. The area under the curve (AUC) of the ROC curve for the model surpassed that of the TNM stage system (0.72 vs. 0.63, p = 0.00028). The model predictions showed close consistency with the actual observations via the calibration plot. Therapeutic response prediction was conducted in high- and low-risk groups and compared with each other. The BC patients with higher risk scores showed more therapeutic resistance than those with a lower risk score.

Conclusion: YTHDF3 was verified as a prognostic biomarker of BC, and a novel YTHDF3-based model was constructed to predict the 5-year OS of BC patients. Our model could be applied to effectively predict the therapeutic response of commonly used agents for BC patients.

Establishment of Prognostic Signatures of N6-Methyladenosine-Related lncRNAs and Their Potential Functions in Hepatocellular Carcinoma Patients

N6-methyladenosine (m6a)-related mRNAs and lncRNAs have been explored for their functions in several cancers. The present study aimed to identify potential signatures of m6a-related lncRNAs in hepatocellular carcinoma (HCC). We downloaded the expression and clinical data from The Cancer Genome Atlas (TCGA) database. The interacted mRNAs and lncRNAs, prognosis-related lncRNAs, potential metabolic pathways of lncRNAs, immune infiltration of various cells, and CD274 (PD-L1) -related lncRNAs were analyzed. Then, in vitro experiments explored the role of AC012073.1 (LOC105377626) in HCC cell lines. We found that candidate 14 lncRNA signatures play functions in HCC maybe by affecting immune infiltration, cell cycle, Notch signaling pathway, etc. LncRNA AC012073.1 (LOC105377626) functions as oncogenic roles in affecting HCC prognosis.

Meclofenamic Acid Restores Gefinitib Sensitivity by Downregulating Breast Cancer Resistance Protein and Multidrug Resistance Protein 7 via FTO/m6A-Demethylation/c-Myc in Non-Small Cell Lung Cancer

Background and Objective

Gefitinib (GE) is a first-line epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI) for patients with advanced non-small cell lung cancer (NSCLC) carrying EGFR activating mutations. However, drug resistance limits the clinical efficacy of gefitinib and ultimately leads to extremely poor clinical benefit. Meclofenamic acid (MA) is a non-steroidal anti-inflammatory drug (NSAID) that relieves moderate and severe pain. In the present study, we aim to determine the MA sensibilization of GE in NSCLC.

Methods

MTT assay was conducted to determine the synergistic effect of MA with GE in GE-sensitive and -resistant cell lines based on the Chou–Talalay method. The Annexin V-PI flow cytometry analysis was conducted to evaluate apoptosis. Western blot assay was used to detect alterations of EGFR downstream molecules. Tritium-labeled GE accumulation analysis was used to determine the efflux activity of GE. Dot blot assays were conducted to determine m6A levels after the MA and GE co-administration. Western blot evaluated the expression of FTO, c-Myc, MRP7, BCRP, and apoptotic proteins.

Results

MA showed a significant synergistic effect with GE in GE-resistant NSCLC cells; co-administration of MA with GE induced caspase-related apoptosis in resistant NSCLC cells. Moreover, EGFR downstream molecules, including Akt and MAPKs pathways, were significantly inhibited by the MA-GE combination. Short-term incubation of MA did not alter the efflux of GE; however, after incubation for 24 h, the accumulation of tritium-labeled GE significantly increased. A mechanism study showed that co-administration of MA and GE significantly downregulated BCRP and MRP7 expression in GE-resistant cells; increased N6-methylation was also observed after co-administration. The FTO/c-Myc was determined as target pathways on MA and GE co-administration mechanisms.

Conclusion

Our findings provide novel therapeutic approaches for GE-resistant NSCLC by combination use with MA through FTO-mediated N6-demethylation.

m6A Regulatory Gene-Mediated Methylation Modification in Glioma Survival Prediction

The median survival of patients with gliomas is relatively short. To investigate the epigenetic mechanisms associated with poor survival, we analyzed publicly available datasets from patients with glioma. This analysis revealed 12 prognosis-related m6A regulatory genes that may be responsible for poor prognosis. These genes may be involved in genomic changes inherent to oxidative phosphorylation, adipogenesis, hedgehog signaling, and Myc signaling. We reconstructed a risk model with univariate and multivariate Cox analyses and identified older age and the m6A risk score as independent risk factors for predicting the prognosis of glioma patients, which is associated with glioma immune infiltration. In conclusion, m6A regulatory genes may serve as both reliable biomarkers and potential targets to increase the chance of survival of patients with glioma.

Diagnostic, Therapeutic, and Prognostic Value of the m6A Writer Complex in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) has poor prognosis and is usually diagnosed only at an advanced stage. Identification of novel biomarkers is critical to early diagnosis and better prognosis for HCC patients. N6-methyladenosine (m⁶A) RNA methylation regulators play important roles in the development of many tumors. However, the m⁶A writer complex, a key executor of m⁶A methylation modification, has not been independently investigated, and its specific bioinformatics analysis has not yet been performed in HCC. In this study, we used multiple public databases to evaluate the diagnostic, therapeutic, and prognostic value of the m⁶A writers in HCC. The results showed that expression levels of METTL3, VIRMA and CBLL1 were significantly increased, while expression levels of METTL14 and ZC3H13 were significantly decreased in HCC, which was closely related to clinicopathological factors, such as tumor stage and prognosis. Bioinformatics further explored the possible underlying mechanisms by which the m⁶A writer complex are involved in activation of tumor-promoting pathways and/or inhibition of tumor-suppressing pathways, including apoptosis, cell cycle, DNA damage response and EMT. Furthermore, we showed that the m⁶A writer complex is correlated with immune cell infiltration and immunoregulator expression in HCC. In conclusion, the m⁶A writer complex may represent a promising biomarker and target that can guide targeted therapy or immunotherapy for HCC patients.

Analysis of m6A Methylation Modification Patterns and Tumor Immune Microenvironment in Breast Cancer

Increasing evidence indicates that the abnormal expression of N6-methyladenosine (m6A) modification is closely related to the epigenetic regulation of immune response in breast cancer (BC). However, the potential roles of m6A modification in the tumor microenvironment (TME) of BC remain unclear. For addressing this issue, we comprehensively analyzed the m6A modification patterns in 983 samples and correlated these modification patterns with TME immune cell infiltration, based on 23 kinds of m6A regulators. Principal component analysis (PCA) was used to construct the m6A scoring system to quantify the modification pattern of m6A of BC individuals. Consequently, three different m6A modification patterns were identified, and the infiltrating characteristics of TME cells were consistent with the three immune phenotypes, including immune rejection, immune inflammation, and immune desert. Besides, our analysis showed that the prognosis of patients could be predicted by evaluating the m6A modification pattern in the tumor. The low m6Ascore corresponded to increased mutation burden and immune activation, while stroma activation and lack of immune infiltration were observed in high m6Ascore subtypes. In addition, a low m6Ascore was associated with enhanced response to anti-PD-1/PD-L1 immunotherapy. In conclusion, the m6A modification pattern was closely related to the BC immune landscape. This well-validated score model of the m6A modification patterns will provide a valuable tool to depict the tumor immune state and guide effective tumor immunotherapy for combating BC.

Isovitexin protects against acute liver injury by targeting PTEN, PI3K and BiP via modification of m6A

Isovitexin (IVT) has been shown to have a potential therapeutic effect on acute liver injury (ALI), but its underlying mechanisms especially the targets remain unclear, which was investigated in the present study. Briefly, the targets of IVT were predicted by bioinformatics and then were verified by multiple examinations using molecular docking, cellular thermal shift assay (CETSA), and Lipopolysaccharide/D-Galactosamine (LPS/D-GalN)-induced ALI animal model. The bioinformatic analysis predicted that the target genes of IVT against ALI were enriched into the PI3K/Akt and ERS-related pathways, in which, molecular docking and CETSA examination verified that the binding sites of IVT likely were PTEN, PI3K and BiP. Furthermore, the possible targets were also verified by animal experiments. The results revealed that IVT significantly ameliorated the hepatic injury, as evidenced by the attenuation of histopathological changes and the reduction in serum aminotransferase and total bilirubin activities. In addition, IVT treatment led to the reduction of PTEN, BiP and ERS-related targets expressions, as well as the elevation of PI3K, Akt and mTOR expressions. Notably, IVT significantly decreased total hepatic m6A level and m6A enrichment of PTEN and BiP, suggesting IVT regulated PTEN and BiP by modulating m6A modification. To sum up, the results indicate that IVT significantly ameliorates ALI, which is attributed to its ability to regulate the PI3K/Akt pathway and ERS by targeting PTEN, PI3K and BiP via modification of m6A. Our finding demonstrates that IVT may be a promising natural medicine for the treatment of ALI.

Deciphering Obesity-Related Gene Clusters Unearths SOCS3 Immune Infiltrates and 5mC/m6A Modifiers in Ossification of Ligamentum Flavum Pathogenesis

BACKGROUND

Ossification of ligamentum flavum (OLF) is an insidious and debilitating heterotopic ossifying disease with etiological heterogeneity and undefined pathogenesis. Obese individuals predispose to OLF, whereas the underlying connections between obesity phenotype and OLF pathomechanism are not fully understood. Therefore, this study aims to explore distinct obesity-related genes and their functional signatures in OLF.

METHODS

The transcriptome sequencing data related to OLF were downloaded from the GSE106253 in the Gene Expression Omnibus (GEO) database. The obesity-related differentially expressed genes (ORDEGs) in OLF were screened, and functional and pathway enrichment analysis were applied for these genes. Furthermore, protein-protein interactions (PPI), module analysis, transcription factor enrichment analysis (TFEA), and experiment validation were used to identify hub ORDEGs. The immune infiltration landscape in OLF was depicted, and correlation analysis between core gene SOCS3 and OLF-related infiltrating immune cells (OIICs) as well as 5mC/m6A modifiers in OLF was constructed.

RESULTS

Ninety-nine ORDEGs were preliminarily identified, and functional annotations showed these genes were mainly involved in metabolism, inflammation, and immune-related biological functions and pathways. Integrative bioinformatic algorithms determined a crucial gene cluster associated with inflammatory/immune responses, such as TNF signaling pathway, JAK-STAT signaling pathway, and regulation of interferon-gamma-mediated signaling. Eight hub ORDEGs were validated, including 6 down-regulated genes (SOCS3, PPARG, ICAM-1, CCL2, MYC, and NT5E) and 2 up-regulated genes (PTGS2 and VEGFA). Furthermore, 14 differential OIICs were identified by ssGSEA and xCell, and SOCS3 was overlapped to be the core gene, which was associated with multiple immune infiltrates (dendritic cells, macrophage, and T cells) and six m6A modifiers as well as four 5mC regulators in OLF. Reduced SOCS3 and FTO expression and up-regulated DNMT1 level in OLF were validated by Western blotting.

CONCLUSION

This study deciphered immune/inflammatory signatures of obesity-related gene clusters for the first time, and defined SOCS3 as one core gene. The crosstalk between 5mC/m6A methylation may be a key mediator of SOCS3 expression and immune infiltration. These findings will provide more insights into molecular mechanisms and therapeutic targets of obesity-related OLF.

N6-methyladenosine-mediated downregulation of miR-374c-5p promotes cadmium-induced cell proliferation and metastasis by targeting GRM3 in breast cancer cells

Cadmium (Cd) is a toxic heavy metal that can facilitate the development and progression of breast cancer (BC). Emerging evidence has indicated that the progression of Cd-exposed BC is related to the dysregulation of microRNAs (miRNAs). The purpose of our study was to investigate the expression pattern and underlying mechanisms of miR-374c-5p in Cd-mediated BC progression. In this study, T-47D cells and MCF-7 cells were treated with different concentrations of Cd (0.1, 1 and 10 μM) for 72 h. MiR-374c-5p expression was downregulated, and transfection of miR-374c-5p mimics significantly decreased BC cell proliferation, migration and invasion induced by 10 μM Cd. Importantly, we used the Cytoscape software plugin cytoHubba to analyse the intersected genes between our RNA-Seq results and the mirDIP database, and six hub genes (CNR1, CXCR4, GRM3, RTN1, SLC1A6 and ZEB1) were identified as potential direct targets of miR-374c-5p in our model; however, luciferase reporter assays indicated that miR-374c-5p only repressed GRM3 by directly binding to its 3'-untranslated region (UTR). Of note, we verified that suppression of N6-methyladenosine (m6A) modification led to miR-374c-5p downregulation by decreasing its RNA transcript stability. Together, these findings demonstrated that m6A modification of pri-miRNA-374c blocks miRNA-374c-5p maturation and then activates GRM3 expression, which drives BC cell metastasis after Cd exposure.

The role of RNA m6A methylation in lipid metabolism

The m6A methylation is the most numerous modification of mRNA in mammals, coordinated by RNA m6A methyltransferases, RNA m6A demethylases, and RNA m6A binding proteins. They change the RNA m6A methylation level in their specific manner. RNA m6A modification has a significant impact on lipid metabolic regulation. The "writer" METTL3/METTL14 and the "eraser" FTO can promote the accumulation of lipids in various cells by affecting the decomposition and synthesis of lipids. The "reader" YTHDF recognizes m6A methylation sites of RNA and regulates the target genes' translation. Due to this function that regulates lipid metabolism, RNA m6A methylation plays a pivotal role in metabolic diseases and makes it a great potential target for therapy.

Role of Main RNA Methylation in Hepatocellular Carcinoma: N6-Methyladenosine, 5-Methylcytosine, and N1-Methyladenosine

RNA methylation is considered a significant epigenetic modification, a process that does not alter gene sequence but may play a necessary role in multiple biological processes, such as gene expression, genome editing, and cellular differentiation. With advances in RNA detection, various forms of RNA methylation can be found, including N6-methyladenosine (m6A), N1-methyladenosine (m1A), and 5-methylcytosine (m5C). Emerging reports confirm that dysregulation of RNA methylation gives rise to a variety of human diseases, particularly hepatocellular carcinoma. We will summarize essential regulators of RNA methylation and biological functions of these modifications in coding and noncoding RNAs. In conclusion, we highlight complex molecular mechanisms of m6A, m5C, and m1A associated with hepatocellular carcinoma and hope this review might provide therapeutic potent of RNA methylation to clinical research.

Identification of m6A-Related lncRNAs Associated With Prognoses and Immune Responses in Acute Myeloid Leukemia

Background: Acute myeloid leukemia (AML) remains the most common type of hematopoietic malignancy in adults and has an unfavorable outcome. Herein, we aimed to construct an N6-methylandenosine (m6A)-related long noncoding RNAs (lncRNAs) signature to accurately predict the prognosis of patients with AML using the data downloaded from The Cancer Genome Atlas (TCGA) database.

Methods: The RNA-seq and clinical data were obtained from the TCGA AML cohort. First, Pearson correlation analysis was performed to identify the m6A-related lncRNAs. Next, univariate Cox regression analysis was used to determine the candidate lncRNAs with prognostic value. Then, feature selection was carried out by Least absolute shrinkage and selection operator (LASSO) analysis, and seven eligible m6A-related lncRNAs were included to construct the prognostic risk signature. Kaplan–Meier and receiver operating characteristic (ROC) curve analyses were performed to evaluate the predictive capacity of the risk signature both in the training and testing datasets. A nomogram was used to predict 1-year, 2-year, and 3-year overall survival (OS) of AML patients. Next, the expression levels of lncRNAs in the signature were validated in AML samples by qRT-PCR. Functional enrichment analyses were carried out to identify probable biological processes and cellular pathways. The ceRNA network was developed to explore the downstream targets and mechanisms of m6A-related lncRNAs in AML.

Results: Seven m6A-related lncRNAs were identified as a prognostic signature. The low-risk group hold significantly prolonged OS. The nomogram showed excellent accuracy of the signature for predicting 1-year, 2-year and 3-year OS (AUC = 0.769, 0.820, and 0.800, respectively). Moreover, the risk scores were significantly correlated with enrichment in cancer hallmark- and malignancy-related pathways and immunotherapy response in AML patients.

Conclusion: We developed and validated a novel risk signature with m6A-related lncRNAs which could predict prognosis accurately and reflect the immunotherapy response in AML patients.

Demethyltransferase FTO alpha-ketoglutarate dependent dioxygenase (FTO) regulates the proliferation, migration, invasion and tumor growth of prostate cancer by modulating the expression of melanocortin 4 receptor (MC4R)

N6-Methyladenosine (m6A) is the most abundant modifications in human messenger RNAs (mRNAs). This study aimed at investigating the function and mechanism of demethyltransferase fat mass and obesity-associated protein (FTO) in prostate cancer(PCa). The expression level of FTO in PCa was detected by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and western blot. Besides, the impacts of FTO on the proliferation, migration and invasion of PCa cells were also detected by cell counting kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU) and transwell assays. Furthermore, we also explored the potential mechanism of FTO in PCa. The results showed that FTO expression was decreased in PCa, and the low expression of FTO showed an obvious relevance to the clinical characteristics. Downregulation of FTO facilitated the proliferation, migration, invasion and tumor growth of PCa cells. Besides, MC4R displayed a remarkably high expression in PCa tissues, whose expression and m6A level were regulated by FTO. Meanwhile, the in vitro experiments revealed that highly expressed FTO partially reversed the facilitating effect of highly expressed MC4R on the malignant phenotype of PCa cells. Overall, FTO was downregulated in PCa and its expression level showed a relevance to the prognosis of PCa patients. Additionally, FTO could regulate the proliferation, migration and invasion of PCa via regulating the expression level of MC4R.

Comprehensive analysis of differences of N6-methyladenosine of lncRNAs between atrazine-induced and normal Xenopus laevis testis

Background

Increasing evidence suggested N⁶-methyladenosine (m⁶A) modification is crucial for male germline development. However, m⁶A modification of lncRNAs gains a little attention in amphibians in recent years. Xenopus laevis (X. laevis) was chosen to be an ideal model organism for testing environmental endocrine disrupting chemicals (EDCs) exposure and resultant effects. Atrazine (AZ) as an endocrine disrupt can effect development of testis in amphibians. Our previous study revealed that m⁶A is a highly conserved modification across the species.

Results

The results of m⁶A sequences showed that m⁶A-methylated lncRNAs enriched in intergenic region in testes of X. laevis. We further examined the differential expression of lncRNAs m⁶A sites in testes of AZ-exposed and compared with that in animals from control group. The results indicated that up to 198 differentially methylated m⁶A sites were detected within 188 lncRNAs, in which 89 significantly up-methylated sites and 109 significantly down-methylated sites. Data from KEGG pathway analysis indicated that AZ-affected lncRNAs m⁶A sites were mainly involved in 10 pathways in which 3 mutual pathways were found in the result of differentially m⁶A-methylated mRNAs.

Conclusions

These findings suggested that differentially m⁶A-methylated lncRNAs and these 3 pathways may act on regulatory roles in abnormal testis development of AZ-exposed X. laevis. This study for the first time provides insights into the profile of lncRNAs m⁶A modifications in amphibian species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41021-021-00223-0.

Deregulation of N6-Methyladenosine RNA Modification and Its Erasers FTO/ALKBH5 among the Main Renal Cell Tumor Subtypes

(1) Background: Methylation of N⁶-adenosine (m⁶A) is the most abundant messenger RNA (mRNA) modification in eukaryotes. We assessed the expression profiles of m⁶A regulatory proteins in renal cell carcinoma (RCC) and their clinical relevance, namely, as potential biomarkers. (2) Methods: In silico analysis of The Cancer Genome Atlas (TCGA) dataset was use for evaluating the expression of the m⁶A regulatory proteins among RCC subtypes and select the most promising candidates for further validation. ALKBH5 and FTO transcript and protein expression were evaluated in a series of primary RCC (n = 120) and 40 oncocytomas selected at IPO Porto. (3) Results: In silico analysis of TCGA dataset disclosed altered expression of the major m⁶A demethylases among RCC subtypes, particularly FTO and ALKBH5. Furthermore, decreased FTO mRNA levels associated with poor prognosis in ccRCC and pRCC. In IPO Porto’s cohort, FTO and ALKBH5 transcript levels discriminated ccRCC from oncocytomas. Furthermore, FTO and ALKBH5 immunoexpression differed among RCC subtypes, with higher expression levels found in ccRCC comparatively to the other RCC subtypes and oncocytomas. (4) Conclusion: We conclude that altered expression of m⁶A RNA demethylases is common in RCC and seems to be subtype specific. Specifically, FTO and ALKBH5 might constitute new candidate biomarkers for RCC patient management, aiding in differential diagnosis of renal masses and prognostication.