Transcriptome -Wide m6A Methylome and Potential Drug Targets Prediction Analysis of Rheumatoid Arthritis
Ontology highlight
ABSTRACT: This research aims at assessing the m6A methylome in the RA peripheral blood mononuclear cells (PBMCs) and perform potential drug targets prediction analysis. Five RA samples and ten control samples were obtained from China-Japan Friendship Hospital. The various expression of m6A methylation and genes in RA and control group were compared with methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq). Bioinformatics explorations were also made to explore the enriched biological roles and paths of the differentially expressed m6A methylation and genes. Potential drug targets prediction and validation were also constructed to provide potential therapeutic targets for RA. In total, 583 dysregulated m6A peaks, of which 295 were greatly upregulated and 288 were greatly downregulated, were identified. Similarly, 1570 differentially genes were acquired by RNA-seq, including 539 upregulated and 1031 downregulated. According to deeper joint exploration, the m6A methylation and mRNA expression degrees of 35 genes varied greatly, of which including 13 hyper-methylated as well as upregulated genes (hyper-up), 12 hyper-methylated and downregulated genes (hyper-down), 5 hypo-methylated as well as upregulated genes (hypo-up), and 5 hypo-methylated and downregulated genes (hypo-down) GO and KEGG explorations showed the enrichment of these distinct genes in “protein transport” “protein binding” and “lysosome”. In addition, the mRNA levels of TUBB2A, IGF2BP3, DYNC1I1 and FOSL1 were increased and consistent with the trend of our sequencing results. While after the treatment of TP and MTX, their mRNA levels were decreased. This research set up the transcriptional map of m6A in RA PBMCs and displayed the hidden association between RNA methylation alternation and associated genes in RA. The outcomes highlight the importance of m6A modification as a gene regulatory system in RA. Furthermore, TUBB2A, IGF2BP3, DYNC1I1 and FOSL1 might be potential therapeutic targets of TP and MTX during RA treatment.
Project description:Background: Methotrexate (MTX) is the first line treatment of rheumatoid arthritis (RA), and methylation changes in bulk T cells have been reported after treatment with MTX. We have investigated cell-type specific DNA methylation changes across the genome in naïve and memory CD4+ T cells before and after MTX treatment of RA patients. DNA methylation profiles of newly diagnosed RA patients (N=9) were assessed by reduced representation bisulfite sequencing. Results: We found that MTX treatment significantly influenced DNA methylation levels at multiple CpG sites in both cell populations. Interestingly, we identified differentially methylated sites annotated to two genes; TRIM15 and SORC2, previously reported to predict treatment outcome in RA patients when measured in bulk T cells. Furthermore, several of the genes, including STAT3, annotated to the significant CpG sites are relevant for RA susceptibility or the action of MTX. Conclusion: We detected CpG sites that were associated with MTX treatment in CD4+ naïve and memory T cells isolated from RA patients. Several of these sites overlap genetic regions previously associated with RA risk and MTX treatment outcome.
Project description:We performed the transcriptome-wide m6A profiling (MeRIP) with poly(A) RNA samples isolated from shoots and roots of wild type (Nipp), FTO homogeneous transgenic (FTO) and demethylation activity dead mutant variant of FTO transgenic (FTOmut) rice plants. The m6A motif and the distribution of the detected m6A sites along transcripts were consistent with reported results from previous m6A sequencing studies in plants. We then analyzed the differentially methylated m6A sites in poly(A) RNA species including those transcribed from gene loci and from repetitive elements (repeat RNAs). Compared to Nipp plants, both shoots and roots of FTO transgenic plants had showed more hypomethylated m6A peaks (hypo-m6A, 222 in shoots and 3313 in roots) than hypermethylated m6A peaks (hyper-m6A, 63 in shoots and 127 in roots) in mRNAs, while FTOmut rice had showed more hyper-m6A in shoots (349 for hyper vs 31 for hypo) but more hypo-m6A in roots (291 for hyper vs 380 for hypo) than the WT. The hypomethylated m6A peaks in FTO transgenic plants were highly enriched within 3'UTR. Gene Ontology (GO) analysis of these identified hypomethylated m6A-containing coding genes in both FTO transgenic rice shoots and roots revealed an enrichment for functional annotations related to “cellular homeostatic process”, “one-carbon and small molecule metabolic process”, and “gene expression”. We also found the extent of m6A hypomethylation was more pronounced for repeat RNAs than for mRNA species.The quantitative RNA-seq analysis that included ERCC spike-in controls confirmed that FTO plants do indeed accumulate higher overall levels of poly(A) RNA than do WT shoots and roots. Above all, FTO mediated m6A demethylation in plant mRNA and repeats RNA, leading to the activation of transcription.
Project description:Purpose: To systematically analyze the overall m6A modification pattern in hyperplastic scars (HS). Methods: The m6A modification patterns in HS and normal skin (NS) tissues were described by m6A sequencing and RNA sequencing, and subsequently bioinformatics analysis was performed. The m6A-related RNA was immunoprecipitated and verified by real-time quantitative PCR. Results: The appearance of 14,791 new m6A peaks in the HS sample was accompanied by the disappearance of 7,835 peaks. The unique m6A-related genes in HS were thus associated with fibrosis-related pathways. We identified the differentially expressed mRNA transcripts in HS samples with hyper-methylated or hypo-methylated m6A peaks. Conclusion: This study is the first to map the m6A transcriptome of human HS, which may help clarify the possible mechanism of m6A-mediated gene expression regulation.
Project description:Purpose: To systematically analyze the overall m6A modification pattern in hyperplastic scars (HS). Methods: The m6A modification patterns in HS and normal skin (NS) tissues were described by m6A sequencing and RNA sequencing, and subsequently bioinformatics analysis was performed. The m6A-related RNA was immunoprecipitated and verified by real-time quantitative PCR. Results: The appearance of 14,791 new m6A peaks in the HS sample was accompanied by the disappearance of 7,835 peaks. The unique m6A-related genes in HS were thus associated with fibrosis-related pathways. We identified the differentially expressed mRNA transcripts in HS samples with hyper-methylated or hypo-methylated m6A peaks. Conclusion: This study is the first to map the m6A transcriptome of human HS, which may help clarify the possible mechanism of m6A-mediated gene expression regulation.
Project description:Glutathione peroxidase 8 (GPX8) is a key regulator of redox homeostasis. Whether its antioxidant activity participates in the regulation of m6A modification is a crucial issue, which has important application value in cancer treatment. MeRIP-seq was used to explore the characteristics of transcriptome-wide m6A modification in GPX8-deficient oral cancer cells in this study. Oxidative stress caused by lack of GPX8 resulted in 1,279 hyper- and 2,287 hypo-methylated m6A peaks and 2,036 differentially expressed genes in GPX8-KO cells. Twenty-eight differentially expressed genes were related to the cell response to oxidative stress, and half of them changed their m6A modification. In GPX8-KO cells, m6A regulators IGF2BP2 and IGF2BP3 were upregulated, while FTO, RBM15, VIRMA, ZC3H13 and YTHDC2 were downregulated. After H2O2 treatment, the expression changes of RBM15, IGF2BP2 and IGF2BP3 were further enhanced. These data indicated that GPX8-mediated redox homeostasis regulated m6A modification, thereby affecting the expression and function of downstream genes. This study highlights the possible significance of GPX8 and the corresponding m6A regulatory or regulated genes as novel targets for antioxidant intervention in cancer therapy
Project description:The aim of the study is to identify differentially methylated positions (DMPs) and regions (DMRs) that predict response to Methotrexate (MTX) in early rheumatoid arthritis (RA) patients. DNA from baseline peripheral blood mononuclear cells was extracted from treatment naive RA patients. DNA methylation, quantified using the Infinium MethylationEPIC, was assessed in relation to response to MTX (combination) therapy (deltaDAS28) over the first 3 months in 69 RA patients.
Project description:N6-methyladenosine (m6A) modification plays a central role in both RNA and tumor biology. Yet comprehensive m6A epitranscriptome profiling of primary tumors remains uncharted. To fill this gap, we profiled the m6A epitranscriptome of 10 non-neoplastic lung (NL) tissues and 51 lung adenocarcinoma (LUAD) tumors, integrating them with corresponding transcriptome, proteome and extensive clinical annotation. A total of 8,030 genes are marked with m6A modification, with variable abundances across tumors. Distinct clusters and genes were found to be linked exclusively to disease progression through m6A, rather than mRNA or protein levels. In comparison with NL tissues, 430 transcripts were hypo-methylated in tumors, and 222 were hyper-methylated. Fully 69% of these had no influence on RNA abundances. Amongst these genes, EML4 emerged as a novel metastatic driver, displaying significant hyper-methylation in tumors. m6A modification promoted EML4 translation and led to widespread EML4 overexpression in primary tumors. Functionally, EML4 modulated cytoskeleton dynamics and cell morphology, triggering cellular protrusions and enhancing cellular motility, local invasion, and metastasis. Clinically, high EML4 protein abundance correlated with features of metastasis. A METTL3 small molecule inhibitor markedly diminished both EML4 m6A and protein abundance, and efficiently suppressed lung metastases in vivo. Our findings unveil a dynamic and functional epitranscriptomic landscape in lung adenocarcinoma, pinpointing the hyper-methylation of EML4 as a key driver of tumor metastasis.
Project description:N6-methyladenosine (m6A) modification plays a central role in both RNA and tumor biology. Yet comprehensive m6A epitranscriptome profiling of primary tumors remains uncharted. To fill this gap, we profiled the m6A epitranscriptome of 10 non-neoplastic lung (NL) tissues and 51 lung adenocarcinoma (LUAD) tumors, integrating them with corresponding transcriptome, proteome and extensive clinical annotation. A total of 8,030 genes are marked with m6A modification, with variable abundances across tumors. Distinct clusters and genes were found to be linked exclusively to disease progression through m6A, rather than mRNA or protein levels. In comparison with NL tissues, 430 transcripts were hypo-methylated in tumors, and 222 were hyper-methylated. Fully 69% of these had no influence on RNA abundances. Amongst these genes, EML4 emerged as a novel metastatic driver, displaying significant hyper-methylation in tumors. m6A modification promoted EML4 translation and led to widespread EML4 overexpression in primary tumors. Functionally, EML4 modulated cytoskeleton dynamics and cell morphology, triggering cellular protrusions and enhancing cellular motility, local invasion, and metastasis. Clinically, high EML4 protein abundance correlated with features of metastasis. A METTL3 small molecule inhibitor markedly diminished both EML4 m6A and protein abundance, and efficiently suppressed lung metastases in vivo. Our findings unveil a dynamic and functional epitranscriptomic landscape in lung adenocarcinoma, pinpointing the hyper-methylation of EML4 as a key driver of tumor metastasis.
Project description:During development and differentiation, enhancers, and not promoters are most dynamic in their DNA methylation status. However, the causal relationship between enhancer activity and methylation is not clear. Here, we describe that during early zebrafish development, enhancer activity has little influence on DNA methylation, and that hypo-methylation is a unique feature of primed enhancers, whereas active enhancers are hyper-methylated. Hypo-methylated enhancers (hypo-enhancers) are enriched close to important transcription factors (TFs) that act later in development, are specifically de-methylated before the maternal-zygotic transition (MZT), and reside in a unique epigenetic environment. Finally, we demonstrate that hypo-enhancers are functionally active and that they are physically associated with the transcriptional start site (TSS) of target-genes, irrespective of target-gene activity. In conclusion, we demonstrate that early development in zebrafish represents a time-window characterized by non-canonical DNA methylation-enhancer relationships, including global DNA hypo-methylation of inactive enhancers and DNA hyper-methylation of active enhancers. DNA methylation profiles of 4 time points during zebrafish early development.
Project description:N6-Methyladenosine (m6A) is the most abundant post-transcriptional modification in eukaryotes, the imbalance of which is reported to be associated with various pathological processes, including drug resistance. In this study, we analyzed the methylated RNA immunoprecipitation combined with next-generation sequencing (MeRIP-seq) data of AML cell line HL60 and its adriamycin-resistant cell line HL60/ADR. We found a total of 40550 peaks, representing 15640 genes in HL60, and a total of 38834 peaks, representing 15285 genes in HL60/ADR. A total of 4437 differentially methylated m6A peaks within 3461 genes have been found between HL60 and HL60/ADR. Among them, 3587 differentially m6A peaks within 2790 genes were hyper-methylated, and 850 m6A peaks within 671 genes were hypo-methylated. KEGG pathway analysis showed that pathways were enriched in tumor and drug-resistant related signaling pathway. Results of MeRIP-seq showed that fold enrichment of global m6A peaks was higher in HL60/ADR compared to HL60. This study provides a framework for the application of comprehensive mRNA m6A profiling towards acute myeloid leukemia cell line (HL60) and its adriamycin-resistant acute myeloid leukemia cell line (HL60/ADR).