m6A Regulatory Genes Research Articles

The m6A and immune regulatory gene signature predicts the prognosis and correlates with immune infiltration of head and neck squamous cell carcinoma

Recent investigations have underscored the epigenetic modulation of the immune response; however, the interplay between RNA N6-methyladenosine (m6A) modification and immunomodulation in head and neck squamous cell carcinoma (HNSC) remains relatively unexplored. To bridge this knowledge gap, we undertook an extensive examination of the potential contributions of m6A modification and immunomodulation in HNSC. We amalgamated and deduplicated 27 m6A -related genes (m6AGs) and 1342 immune regulation-related genes (IMRGs), resulting in a comprehensive dataset encompassing 1358 genes. This dataset was scrutinized for m6A modification and immunomodulatory patterns within HNSC specimens. Employing Cox regression analysis and the Least Absolute Shrinkage and Selection Operator (LASSO) technique, we developed a prognostic risk model for m6A regulator-mediated methylation modification and immunomodulation-related differentially expressed genes (m6A&IMRDEGs). Our differential expression analysis delineated 29 m6A&IMRDEGs, and Weighted Gene Co-expression Network Analysis (WGCNA) elucidated two module genes (IL11 and MMP13) subjected to correlation analysis. The prognostic prediction models revealed that the clinical predictive efficacy peaked for 1-year forecasts, followed sequentially by 3-year and 5-year predictions. The risk scores derived from the model adeptly categorized HNSC patients into high- and low-risk cohorts, with the high-risk group exhibiting a more unfavorable prognosis. Protein-Protein Interaction (PPI) analysis identified 7 hub genes implicated in m6A and immune regulation, namely BPIFB1, BPIFB2, GP2, MUC5B, MUC7, PIP, and SCGB3A1. Furthermore, we noted marked disparities in the expression profiles of 18 immune cell types between the high- and low-risk groups. Our results substantiate that the clustering subpopulations and risk models associated with m6A and immune regulatory genes portend a poor prognosis in HNSC. The risk score emerges as a potent prognostic biomarker and predictive metric for HNSC patients. A thorough assessment of m6A and immune regulatory genes in HNSC will augment our comprehension of the tumor immune microenvironment and facilitate the advancement of HNSC therapeutics.

Transcriptome-Wide Identification of m6A Writers, Erasers and Readers and Their Expression Profiles under Various Biotic and Abiotic Stresses in Pinus massoniana Lamb.

N6-methyladenosine (m6A) RNA modification is the most prevalent form of RNA methylation and plays a crucial role in plant development. However, our understanding of m6A modification in Masson pine (Pinus massoniana Lamb.) remains limited. In this study, a complete analysis of m6A writers, erasers, and readers in Masson pine was performed, and 22 m6A regulatory genes were identified in total, including 7 m6A writers, 7 m6A erases, and 8 readers. Phylogenetic analysis revealed that all m6A regulators involved in Masson pine could be classified into three distinct groups based on their domains and motifs. The tissue expression analysis revealed that the m6A regulatory gene may exert a significant influence on the development of reproductive organs and leaves in Masson pine. Moreover, the results from stress and hormone expression analysis indicated that the m6A regulatory gene in Masson pine might be involved in drought stress response, ABA-signaling-pathway activation, as well as resistance to Monochamus alternatus. This study provided valuable and anticipated insights into the regulatory genes of m6A modification and their potential epigenetic regulatory mechanisms in Masson pine.

Identifying subtypes and developing prognostic models based on N6-methyladenosine and immune microenvironment related genes in breast cancer

Breast cancer (BC) is the most prevalent cancer in women globally. The tumor microenvironment (TME), comprising epithelial tumor cells and stromal elements, is vital for breast tumor development. N6-methyladenosine (m6A) modification plays a key role in RNA metabolism, influencing its various aspects such as stability and translation. There is a notable link between m6A methylation and immune cells in the TME, although this relationship is complex and not fully deciphered. In this research, BC expression and clinicopathological data from TCGA were scrutinized to assess expression profiles, mutations, and CNVs of 31 m6A genes and immune microenvironment-related genes, examining their correlations, functions, and prognostic impacts. Lasso and Cox regression identified prognostic genes for constructing a nomogram. Single-cell analyses mapped the distribution and patterns of these genes in BC cell development. We investigated associations between gene-derived risk scores and factors like immune infiltration, TME, checkpoints, TMB, CSC indices, and drug response. As a complement to computational analyses, in vitro experiments were conducted to confirm these expression patterns. We included 31 m6A regulatory genes and discovered a correlation between these genes and the extent of immune cell infiltration. Subsequently, a 7-gene risk score was generated, encompassing HSPA2, TAP1, ULBP2, CXCL1, RBP1, STC2, and FLT3. It was observed that the low-risk group exhibited better overall survival (OS) in BC, with higher immune scores but lower tumor mutational burden (TMB) and cancer stem cell (CSC) indices, as well as lower IC50 values for commonly used drugs. To enhance clinical applicability, age and stage were incorporated into the risk score, and a more comprehensive nomogram was constructed to predict OS. This nomogram was validated and demonstrated good predictive performance, with area under the curve (AUC) values for 1-year, 3-year, and 5-year OS being 0.848, 0.807, and 0.759, respectively. Our findings highlight the profound impact of prognostic-related genes on BC immune response and prognostic outcomes, suggesting that modulation of the m6A-immune pathway could offer new avenues for personalized BC treatment and potentially improve clinical outcomes.

Antiretroviral Therapy Suppresses RNA N6-Methyladenosine Modification in Peripheral Blood Mononuclear Cells from HIV-1-Infected Individuals.

RNA N6-methyladenosine (m6A) modification is important for regulating gene expression and innate immune responses to viral infection. HIV-1 in vitro infection induces a significant increase in m6A modification of cellular RNA; however, whether m6A levels of cellular RNA are affected by HIV-1 replication or by antiretroviral therapy (ART) in infected individuals remains unknown. Using dot blot or enzyme-linked immunosorbent assay, we measured RNA m6A levels of peripheral blood mononuclear cells (PBMCs) from healthy donors or HIV-1-infected individuals with or without ART. Using a reverse transcription-quantitative polymerase chain reaction array, we quantified expression levels of 84 type-I interferon (IFN-I)-responsive genes in PBMCs from some individuals of these three groups. RNA m6A levels in PBMCs from HIV-1 viremic patients (n = 10) were significantly higher (p ≤ .0001) compared with ART-treated individuals (n = 22) or 1.5-fold higher compared with healthy donors (n = 14). However, the increase in RNA m6A levels did not correlate with changes in the expression of 10 m6A-regulatory genes. We found significant upregulation and downregulation in the expression of several IFN-I-responsive genes from HIV-1 viremic patients (n = 4) and ART-treated patients (n = 6) compared with healthy donors (n = 5), respectively. Our results suggest that post-transcriptional m6A modification may contribute to the regulation of IFN-I-responsive gene expression during HIV-1 infection and ART.

Integrated analysis of m6A regulator-mediated RNA methylation modification patterns and immune characteristics in Sjögren's syndrome

The epigenetic modifier N6-methyladenosine (m6A), recognized as the most prevalent internal modification in messenger RNA (mRNA), has recently emerged as a pivotal player in immune regulation. Its dysregulation has been implicated in the pathogenesis of various autoimmune conditions. However, the implications of m6A modification within the immune microenvironment of Sjögren's syndrome (SS), a chronic autoimmune disorder characterized by exocrine gland dysfunction, remain unexplored. Herein, we leverage an integrative analysis combining public database resources and novel sequencing data to investigate the expression profiles of m6A regulatory genes in SS. Our cohort comprised 220 patients diagnosed with SS and 62 healthy individuals, enabling a comprehensive evaluation of peripheral blood at the transcriptomic level. We report a significant association between SS and altered expression of key m6A regulators, with these changes closely tied to the activation of CD4+ T cells. Employing a random forest (RF) algorithm, we identified crucial genes contributing to the disease phenotype, which facilitated the development of a robust diagnostic model via multivariate logistic regression analysis. Further, unsupervised clustering revealed two distinct m6A modification patterns, which were significantly associated with variations in immunocyte infiltration, immune response activity, and biological function enrichment in SS. Subsequently, we proceeded with a screening process aimed at identifying genes that were differentially expressed (DEGs) between the two groups distinguished by m6A modification. Leveraging these DEGs, we employed weight gene co-expression network analysis (WGCNA) to uncover sets of genes that exhibited strong co-variance and hub genes that were closely linked to m6A modification. Through rigorous analysis, we identified three critical m6A regulators - METTL3, ALKBH5, and YTHDF1 - alongside two m6A-related hub genes, COMMD8 and SRP9. These elements collectively underscore a complex but discernible pattern of m6A modification that appears to be integrally linked with SS's pathogenesis. Our findings not only illuminate the significant correlation between m6A modification and the immune microenvironment in SS but also lay the groundwork for a deeper understanding of m6A regulatory mechanisms. More importantly, the identification of these key regulators and hub genes opens new avenues for the diagnosis and treatment of SS, presenting potential targets for therapeutic intervention.

Genome-wide identification of the N6-methyladenosine regulatory genes reveals NtFIP37B increases drought resistance of tobacco (Nicotiana tabacum L.)

BackgroundN6-methyladenosine (m6A) is one of the common internal RNA modifications found in eukaryotes. The m6A modification can regulate various biological processes in organisms through the modulation of alternative splicing, alternative polyadenylation, folding, translation, localization, transport, and decay of multiple types of RNA, without altering the nucleotide sequence. The three components involved in m6A modification, namely writer, eraser, and reader, mediate the abundance of RNA m6A modification through complex collaborative actions. Currently, research on m6A regulatory genes in plants is still in its infancy.ResultsIn this study, we identified 52 candidate m6A regulatory genes in common tobacco (Nicotiana tabacum L.). Gene structure, conserved domains, and motif analysis showed structural and functional diversity among different subgroups of tobacco m6A regulatory genes. The amplification of m6A regulatory genes were mainly driven by polyploidization and dispersed duplication, and duplicated genes evolved through purified selection. Based on the potential regulatory network and expression pattern analysis of m6A regulatory genes, a significant number of m6A regulatory genes might play important roles in growth, development, and stress response processes. Furthermore, we have confirmed the critical role of NtFIP37B, an m6A writer gene in tobacco, in enhancing drought resistance.ConclusionsThis study provides useful information for better understanding the evolution of m6A regulatory genes and the role of m6A modification in tobacco stress response, and lays the foundation for further elucidating the function of m6A regulatory genes in tobacco.

Identification of m6A-Related Biomarkers in Systemic Lupus Erythematosus: A Bioinformation-Based Analysis.

Systemic Lupus Erythematosus (SLE), a prototypical autoimmune disorder, presents a challenge due to the absence of reliable biomarkers for discerning organ-specific damage within SLE. A growing body of evidence underscores the pivotal involvement of N6-methyladenosine (m6A) in the etiology of autoimmune conditions. The datasets, which primarily encompassed the expression profiles of m6A regulatory genes, were retrieved from the Gene Expression Omnibus (GEO) repository. The optimal model, selected from either Random Forest (RF) or Support Vector Machine (SVM), was employed for the development of a predictive nomogram model. To identify pivotal genes associated with SLE, a comprehensive screening process was conducted utilizing LASSO, SVM-RFE, and RF techniques. Within the realm of SLE susceptibility, Weighted Gene Co-expression Network Analysis (WGCNA) was harnessed to delineate relevant modules and hub genes. Additionally, MeRIP-qPCR assays were performed to elucidate key genes correlated with m6A targets. Furthermore, a Mendelian randomization study was conducted based on genome-wide association studies to assess the causative influence of MMP9 on ischemic stroke (IS), which is not only a severe cerebrovascular event but also a common complication of SLE. Twelve m6A regulatory genes was identified, demonstrating statistical significance (p < 0.05) and utilized for constructing a nomogram model using the RF algorithm. EPSTI1, USP18, HP, and MMP9, as the hub genes, were identified. MMP9 uniquely correlates with m6A modification and was causally linked to an increased risk of IS, as indicated by our inverse variance weighting analysis showing an odds ratio of 1.0134 (95% CI=1.0004-1.0266, p = 0.0440). Our study identified twelve m6A regulators, shedding light on the molecular mechanisms underlying SLE risk genes. Importantly, our analysis established a causal relationship between MMP9, a key m6A-related gene, and ischemic stroke, a common complication of SLE, thereby providing critical insights for presymptomatic diagnostic approaches.

Deciphering the Divergent Gene Expression Landscapes of m6A/m5C/m1A Methylation Regulators in Hepatocellular Carcinoma Through Single-Cell and Bulk RNA Transcriptomic Analysis.

RNA modifications mediated by the m6A, m1A, and m5C regulatory genes are crucial for the progression of malignancy. This study aimed to explore the expression of regulator genes for m6A/m5C/m1A methylation at the single-cell level and to validate their expression in cancerous and adjacent para-cancerous liver tissues of adult patients with HCC who underwent tumor resection. The bulk sequencing from The Cancer Genome Atlas (TCGA) database and the single-cell RNA sequencing (scRNA-seq) data obtained from the Gene Expression Omnibus (GEO) database were used to identify the dysregulated m6A/m5C/m1A genes for hepatocellular carcinoma (HCC). A real-time polymerase chain reaction (real-time PCR) was used to measure the expression of dysregulated m6A/m5C/m1A genes in collected human HCC tissues and compared with adjacent para-cancerous liver tissues. Immune cell infiltration with these significantly expressed methylation-related genes was evaluated using Timer2.0. A discrepancy in m6A/m5C/m1A gene expression was observed between bulk sequencing and scRNA-seq. The clustered heatmap of the scRNA-seq-identified dysregulated m6A/m5C/m1A genes in TCGA cohort revealed heterogeneous expression of these methylation regulators within the cancer, whereas their expression in the adjacent liver tissues was more homogeneous. The real-time PCR validated the significant overexpression of DNMT1, NSUN5, TRMT6, IGF2BP1, and IGFBP3, which were identified using scRNA-seq, and IGFBP2, which was identified using bulk sequencing. These dysregulated methylation genes are mainly correlated with the infiltration of natural killer cells. This study suggests that cellular diversity inside tumors contributes to the discrepancy in the expression of methylation regulator genes between traditional bulk sequencing and scRNA-seq. This study identified five regulatory genes that will be the focus of further studies regarding the function of m6A/m5C/m1A in HCC.

HNRNPC mediated m6A methylation of 5-methyltetrahydrofolate-homocysteine methyltransferase and involved in the occurrence of RSA

N6‐methyladenosine methylated modification has been shown to play roles in recurrent spontaneous abortion. We aimed to explore role of heterogeneous nuclear ribonucleoprotein C in the occurrence of recurrent spontaneous abortion. We collected embryonic villous tissues from 3 patients with recurrent spontaneous abortion (RSA group) and 3 normal control pregnancy patients. Methylated RNA immunoprecipitation sequencing, RNA sequencing, methylated RNA immunoprecipitation quantitative PCR were conducted to detect the differentially expressed m6A methylation modification gene and regulatory gene in patients with recurrent spontaneous abortion. Methylated RNA immunoprecipitation sequencing and RNA sequencing results showed that the mRNA expression level of heterogeneous nuclear ribonucleoprotein C significantly decreased in RSA group and mRNA expression level of 5-methyltetrahydrofolate-homocysteine methyltransferase increased. Real-time quantitative PCR confirmed the differential expression of heterogeneous nuclear ribonucleoprotein C and 5-methyltetrahydrofolate-homocysteine methyltransferase. Methylated RNA immunoprecipitation quantitative PCR result showed that mRNA m6A modification level of 5-methyltetrahydrofolate-homocysteine methyltransferase decreased in RSA group. The results of western blotting, real-time quantitative PCR, immunofluorescence, matrigel invasion and wound healing assays indicated that heterogeneous nuclear ribonucleoprotein C might regulate the expression of 5-methyltetrahydrofolate-homocysteine methyltransferase by mediating m6A modification, thereby reducing the proliferation and migration of trophoblast cell line, ultimately leading to the occurrence of recurrent spontaneous abortion.

M6A RNA MODIfiCATION AND A DRIVER OF CELLULAR PLASTICITY AND INTRA-TUMOURAL HETEROGENEITY IN GLIOBLASTOMA

Abstract AIMS The immense cellular plasticity of glioblastoma is a barrier to successful therapy. N6-methyladenosine (m6A) methylation of RNA is increasingly recognised as a key driver of cellular differentiation. Crucially, it is a reversible and dynamic modification that alters cellular phenotype via the regulation of mRNA metabolism. We hypothesise that m6A methylation mediates glioblastoma plasticity, regulating transcriptional networks and thus influencing cellular adaptation under different microenvironmental conditions. Our aims are to assess whether total m6A levels and m6A-regulatory gene expression are responsive to change in oxygenation levels and cellular confluency. Furthermore, we aim to characterise the spatial distribution of m6A in 3D cell culture models. METHOD We have investigated m6A levels in patient-derived cell lines from the invasive region of adult glioblastoma. Mass spectrometry and anti-m6A immunofluorescence confocal microscopy was used to quantify total m6A levels. We have compared total m6A levels in normoxic versus hypoxic culture conditions and in low versus high confluency 2D cell culture, including cell migration during a scratch assay. We have measured expression levels of key m6A-regulatory genes using rtPCR. RESULTS Total m6A levels and regulatory gene expression were not significantly different in normoxic versus hypoxic conditions. However, m6A levels were decreased in high confluency versus low confluency conditions. We will further characterise the spatial distribution of m6A in 3D cell culture models using immunofluorescence confocal microscopy. CONCLUSIONS We demonstrate a degree of variability in m6A levels in response to micro-environmental conditions, such as cellular confluency. These findings indicate a role for m6A in tumoral heterogeneity and cellular adaptation.

Analyses of m6A regulatory genes and subtype classification in atrial fibrillation.

To explore the role of m6A regulatory genes in atrial fibrillation (AF), we classified atrial fibrillation patients into subtypes by two genotyping methods associated with m6A regulatory genes and explored their clinical significance. We downloaded datasets from the Gene Expression Omnibus (GEO) database. The m6A regulatory gene expression levels were extracted. We constructed and compared random forest (RF) and support vector machine (SVM) models. Feature genes were selected to develop a nomogram model with the superior model. We identified m6A subtypes based on significantly differentially expressed m6A regulatory genes and identified m6A gene subtypes based on m6A-related differentially expressed genes (DEGs). Comprehensive evaluation of the two m6A modification patterns was performed. The data of 107 samples from three datasets, GSE115574, GSE14975 and GSE41177, were acquired from the GEO database for training models, comprising 65 AF samples and 42 sinus rhythm (SR) samples. The data of 26 samples from dataset GSE79768 comprising 14 AF samples and 12 SR samples were acquired from the GEO database for external validation. The expression levels of 23 regulatory genes of m6A were extracted. There were correlations among the m6A readers, erasers, and writers. Five feature m6A regulatory genes, ZC3H13, YTHDF1, HNRNPA2B1, IGFBP2, and IGFBP3, were determined (p < 0.05) to establish a nomogram model that can predict the incidence of atrial fibrillation with the RF model. We identified two m6A subtypes based on the five significant m6A regulatory genes (p < 0.05). Cluster B had a lower immune infiltration of immature dendritic cells than cluster A (p < 0.05). On the basis of six m6A-related DEGs between m6A subtypes (p < 0.05), two m6A gene subtypes were identified. Both cluster A and gene cluster A scored higher than the other clusters in terms of m6A score computed by principal component analysis (PCA) algorithms (p < 0.05). The m6A subtypes and m6A gene subtypes were highly consistent. The m6A regulatory genes play non-negligible roles in atrial fibrillation. A nomogram model developed by five feature m6A regulatory genes could be used to predict the incidence of atrial fibrillation. Two m6A modification patterns were identified and evaluated comprehensively, which may provide insights into the classification of atrial fibrillation patients and guide treatment.

Genome-Wide Identification and Expression Analysis of m6A Writers, Erasers, and Readers in Litchi (Litchi chinensis Sonn.).

N6-methyladenosine (m6A) RNA modification is the most prevalent type of RNA methylation and plays a pivotal role in the development of plants. However, knowledge of the m6A modification in litchi remains limited. In this study, a complete analysis of m6A writers, erasers, and readers in litchi was performed and 31 litchi m6A regulatory genes were identified in total, including 7 m6A writers, 12 m6A erases, and 12 readers. Phylogeny analysis showed that all three of the kinds of litchi m6A regulatory proteins could be divided into three groups; domains and motifs exhibited similar patterns in the same group. MiRNA target site prediction showed that 77 miRNA target sites were located in 25 (80.6%) litchi m6A regulatory genes. Cis-elements analysis exhibited that litchi m6A regulatory genes were mainly responsive to light and plant hormones, followed by environmental stress and plant development. Expression analysis revealed litchi m6A regulatory genes might play an important role during the peel coloration and fruit abscission of litchi. This study provided valuable and expectable information of litchi m6A regulatory genes and their potential epigenetic regulation mechanism in litchi.

YTHDF1 amplification is correlated with worse outcome and lower immune cell infiltrations in breast cancer.

N6-methyladenosine (m6A) is a common RNA modification on eukaryotic mRNA and some of the m6A regulatory proteins play a crucial role in breast cancer. However, the copy number variations for m6A regulatory proteins and their role in clinicopathological characteristics and survival in breast cancer remain unclear. In this study, we screened the m6A related genes alterations in breast cancer by analyzing the Molecular Taxonomy of Breast Cancer International Consortium and The Cancer Genome Atlas database, and further analyzed the clinical prognostic value of YTHDF1 amplification. The YTH domain family (YTHDF3 and YTHDF1) amplification exhibited higher alteration rates among 10 m6A regulatory genes. YTHDF1 and YTHDF3 amplification resulted in higher mRNA expression (P< 0.0001). Protein expression of YTHDF1 and YTHDF3 were higher in breast cancer (P< 0.0001). YTHDF1 amplification presented a high correlation with worse clinicopathological characteristics and overall survival in patients with breast cancer. Cox regression analysis showed that YTHDF1 amplification was an independent risk factor for 10-year overall survival in breast cancer (Hazard ratio: 1.663; 95% confidence interval: 1.298-2.131; P< 0.001). Gene set enrichment analysis revealed that the downstream target of YTHDF1 may be related to MYC signaling regulation and T cell differentiation. Moreover, YTHDF1 amplification and high expression resulted in lower immune cell infiltration. YTHDF1 knockdown retrained proliferation, migration and invasion in breast cancer cells in vitro. We found significant worse clinical characteristics and lower immune infiltrates in patients with YTHDF1 amplification. The findings indicate that YTHDF1 amplification may be a potential target for the treatment of breast cancer.

The functions of N6-methyladenosine (m6A) RNA modifications in colorectal cancer.

Colorectal cancers (CRC), which includes colon cancer (CC) and rectal cancer (RC), are some of the most common malignant tumors that are prone to distant metastasis. Its high incidence rate and high mortality rate have attracted much attention. In recent years, epigenetics has attracted increasing attention and has been the focus of many research studies. N6-methyladenosine(m6A) RNA modifications can modify eukaryotic mRNA to impact metabolism. The changes in the m6A regulatory genes are related to the occurrence and development of CRC and play an important role in the pathogenesis of CRC. The effect of m6A RNA modification is regulated by its related regulatory factors ("writer", "eraser", "reader"). In this review, we comprehensively analyzed the effect of m6A methylation on CRC and the relationship between the expression of related regulatory factors and the development and occurrence of CRC. Then, we summarized the roles of m6A and its regulatory factors in CRC and its potential clinical value, which provides a basis for further research on the mechanism of m6A methylation in CRC.

Transcriptome-wide m6A methylome analysis uncovered the changes of m6A modification in oral pre-malignant cells compared with normal oral epithelial cells.

As the most common post-transcriptional RNA modification, m6A methylation extensively regulates the structure and function of RNA. The dynamic and reversible modification of m6A is coordinated by m6A writers and erasers. m6A reader proteins recognize m6A modification on RNA, mediating different downstream biological functions. mRNA m6A modification and its corresponding regulators play an important role in cancers, but its characteristics in the precancerous stage are still unclear. In this study, we used oral precancerous DOK cells as a model to explore the characteristics of transcriptome-wide m6A modification and major m6A regulator expression in the precancerous stage compared with normal oral epithelial cell HOEC and oral cancer cell SCC-9 through MeRIP-seq and RT-PCR. Compared with HOEC cells, we found 1180 hyper-methylated and 1606 hypo-methylated m6A peaks and 354 differentially expressed mRNAs with differential m6A peaks in DOK cells. Although the change of m6A modification in DOK cells was less than that in SCC-9 cells, mRNAs with differential m6A in both cell lines were enriched into many identical GO terms and KEGG pathways. Among the 20 known m6A regulatory genes, FTO, ALKBH5, METTL3 and VIRMA were upregulated or downregulated in DOK cells, and the expression levels of 10 genes such as METTL14/16, FTO and IGF2BP2/3 were significantly changed in SCC-9 cells. Our data suggest that precancerous cells showed, to some extent, changes of m6A modification. Identifying some key m6A targets and corresponding regulators in precancerous stage may provide potential intervention targets for the prevention of cancer development through epigenetic modification in the future.

Cluster classification and clinical prognostic modeling based on m6A RNA methylation regulators in liver cancer

Objective: Cluster classification based on m6A methylation regulators and construct prognostic evaluation model. Methods: Utilizing consensus cluster to classify the liver cancer samples form TCGA based on the expression of 13 m6A methylation regulators, and verify the function and prognostic significance of the clustered subtypes. Marker genes were further screened to construct a risk prediction model for evaluating the prognosis of liver cancer patients. Results: The two clustered subtypes based on m6A methylation regulators showed significant differences in the prognosis value of liver cancer patients (P=0.048), and 38 prognostic markers related to m6A methylation in liver cancer were screened from the subgroup with poor prognosis. Two m6A regulatory genes, YTHDF1 and YTHDF2, are proved with adverse prognosis by univariate cox analysis (P<0.05, Hazard ratio>1). We used Lasso regression method to build risk assessment model and effectively predicted the prognosis status of liver cancer patients within 4 years (4-year AUC=0.685, 3-year AUC=0.669). Moreover, the assessment model was validated in another dataset of Asia liver cancer patients. Conclusion: The study provided ideas for studying m6A methylation in liver cancer, and the risk prediction model can be used to evaluate the short-term prognosis of liver cancer patients.

The diagnostic significance of integrating m6A modification and immune microenvironment features based on bioinformatic investigation in aortic dissection.

PurposeThe aim of this study was to investigate the role of m6A modification and the immune microenvironment (IME) features in aortic dissection (AD) and establish a clinical diagnostic model for AD based on m6A and IME factors.MethodsGSE52093, GSE98770, GSE147026, GSE153434, and GSE107844 datasets were downloaded from the GEO database. The expression of 21 m6A genes including m6A writers, erasers, readers, and immune cell infiltrates was analyzed in AD and healthy samples by differential analysis and ssGSEA method, respectively. Both correlation analyses between m6A genes and immune cells were conducted by Pearson and Spearman analysis. XGboost was used to dissect the major m6A genes with significant influences on AD. AD samples were classified into two subgroups via consensus cluster and principal component analysis (PCA) analysis, respectively. Among each subgroup, paramount IME features were evaluated. Random forest (RF) was used to figure out key genes from AD and healthy shared differentially expressed genes (DEGs) and two AD subgroups after gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, we constructed an AD diagnostic model combining important m6A regulatory genes and assessed its efficacy.ResultsAmong 21 m6A genes, WTAP, HNRNPC, and FTO were upregulated in AD samples, while IGF2BP1 was downregulated compared with healthy samples. Immune cell infiltrating analysis revealed that YTHDF1 was positively correlated with γδT cell level, while FTO was negatively correlated with activated CD4+ T cell abundance. FTO and IGF2BP1 were identified to be crucial genes that facilitate AD development according to the XGboost algorithm. Notably, patients with AD could be classified into two subgroups among which 21 m6A gene expression profiles and IME features differ from each other via consensus cluster analysis. The RF identified SYNC and MAPK1IP1L as the crucial genes from common 657 shared common genes in 1,141 DEGs between high and low m6A scores of AD groups. Interestingly, the AD diagnostic model coordinating SYNC and MAPK1IP1L with FTO and IGF2BP1 performed well in distinguishing AD samples.ConclusionThis study indicated that FTO and IGF2BP1 were involved in the IME of AD. Integrating FTO and IGF2BP1 and MAPK1IP1L key genes in AD with a high m6A level context would provide clues for forthcoming AD diagnosis and therapy.

Bioinformatics Analysis Revealing the Correlation between NF-κB Signaling Pathway and Immune Infiltration in Gastric Cancer.

Although the emerging of immunotherapy conferred a new landscape of gastric cancer (GC) treatment, its response rate was of significant individual differences. Insight into GC immune microenviroment may contribute to breaking the dilemma. To this end, the enrichment score of NF-κB signaling pathway was calculated in each GC sample from The Cancer Genome Atlas (TCGA) via ssGSEA algorithm, and its association with immune infiltration was estimated. Based on NF-κB-related genes, a risk score was established and its involvement in immune infiltration, tumor mutational burden (TMB), and N6-methyladenosine (M6A) modification was analyzed in GC. The results showed that NF-κB signaling pathway promoted the infiltration of immune cells in GC. In addition, GC samples were divided into low- and high-risk groups according to a seven-gene (CARD11, CCL21, GADD45B, LBP, RELB, TRAF1, and VCAM1) risk score. Although the high-risk group displayed high immune infiltration and high expression of M6A regulatory genes, it remains in an immunosuppressive microenviroment and whereby suffers a poorer outcome. Of note, most of hub genes were related to immune infiltration and could serve as an independent prognostic biomarker. Conclusively, our study emphasized the crucial role of NF-κB signaling pathway in GC immune microenviroment and provided several candidate genes that may participate in immune infiltration.

Comprehensive Analysis of N6-Methyladenosine Regulatory Genes from Citrus grandis and Expression Profilings in the Fruits of “Huajuhong” (C. grandis “Tomentosa”) during Various Development Stages

Citrus grandis “Tomentosa” (“Huajuhong”) is a famous traditional Chinese medicine. The aim of the present study is to provide a comprehensive characterization of the m6A regulatory genes from C. grandis, and examine their expression patterns in fruits of C. grandis “Tomentosa” during various developmental stages. A total of 26 N6-methyladenosine (m6A) regulatory proteins were identified from the genome of C. grandis, which were distributed across nine chromosomes in C. grandis. Phylogenetic relationships revealed that all m6A regulatory genes were divided into groups of m6A writers, erasers, and readers. The m6A writer groups included CgMTA, CgMTB, and CgMTC three MTs (methyltransferases), one CgVIR (virilizer), one CgHAKAI (E3 ubiquitin ligase HAKAI), and one CgFIP37 (FKBP interacting protein 37). Moreover, 10 CgALKBH (α-ketoglutarate-dependent dioxygenase homolog) members (numbered from CgALKBH1 to CgALKBH10) and 10 CgECT (C-terminal region) members (numbered from CgECT1 to CgECT10) in C. grandis were identified as m6A erasers and readers, respectively. The domain structures and motif architectures among the groups of m6A writers, erasers, and readers were diverse. Cis-acting elements in the promoters of the 26 m6A regulatory genes predicted that the abscisic acid-responsive (ABA) element (ABRE) was present on the promoters of 19 genes. In addition, the expression profiles of all m6A regulatory genes were examined in the fruits of two varieties of C. grandis “Tomentosa” during different growth stages to give basic hints for further investigation of the function of the N6-methyladenosine regulatory genes in C. grandis “Tomentosa”.

Identification of epitranscriptomic methylation marker genes in Arabidopsis and their expression profiling in response to developmental, anatomical, and environmental modulations

RNA methylation has already conquered attention in cancer research and anticipated a similar contribution in plant stress modulation mechanism. Among all kinds of mRNA methylation, N6-methyladenosine (m6A) modifications are ubiquitous and conserved, making m6A regulatory genes worthwhile to investigate. This present study aims to identify m6A regulatory genes in Arabidopsis and their correlation with stress conditions. A total of 29 m6A regulatory genes (6 writers, 6 erasers, and 17 readers) were identified in Arabidopsis that encode 55 proteins (8 writers, 14 erasers, and 33 readers). Among these 29 identified genes, one eraser gene pair and two reader gene pairs are found to have segmental duplication. Most of the m6A regulatory genes are highly expressed in callus tissue. In senescence, writers and readers are shown higher expression than erasers. In wounding stress, all six writers such as AtMTA, AtMTB, AtMTC, AtHAKAI, AtVIRILIZER, and AtFIP37 were constantly downregulated; two readers AtCPSF30-L3, and AtECT10 were constantly up-regulated and three erasers such as AtALKBH9A, AtALKBH10A, and AtALKBH8B were down-regulated in both early and late stages. Differential expression pattern of m6A regulatory genes have substantiated their significant correlation in the stress modulation pathway and demands further investigation. This study will encourage researchers to functionally validate and elucidate the molecular mechanism of plants utilizing RNA methylation in stress tolerance, thus, will facilitate crop engineering.