RNA modifications, especially methylation of the N6 position of adenosine （A）--m6A, rep- resent an emerging research frontier in RNA biology. With the rapid development of high-throughput sequencing technology, in-d...RNA modifications, especially methylation of the N6 position of adenosine （A）--m6A, rep- resent an emerging research frontier in RNA biology. With the rapid development of high-throughput sequencing technology, in-depth study of m6A distribution and function relevance becomes feasible. However, a robust method to effectively identify m6A-modified regions has not been available yet. Here, we present a novel high-efficiency and user-friendly analysis pipeline called MeRIP-PF for the signal identification of MeRIP-Seq data in reference to controls. MeRIP-PF provides a statistical P-value for each identified m6A region based on the difference of read distribution when compared to the con- trols and also calculates false discovery rate （FDR） as a cut offto differentiate reliable m6A regions from the background. Furthermore, MeRIP-PF also achieves gene annotation ofm6A signals or peaks and produce outputs in both XLS and graphical format, which are useful for further study. MeRIP-PF is implemented in Perl and is freely available at http：//software.big.ac.cn/MeRIP-PF.html.展开更多
Evidence showed that N6-methyladenosine(m^(6)A)modification plays a pivotal role in influencing RNA fate and is strongly associated with cell growth and developmental processes in many species.However,no information r...Evidence showed that N6-methyladenosine(m^(6)A)modification plays a pivotal role in influencing RNA fate and is strongly associated with cell growth and developmental processes in many species.However,no information regarding m^(6)A modification in Eimeria tenella is currently available.In the present study,we surveyed the transcriptome-wide prevalence of m^(6)A in sporulated oocysts and unsporulated oocysts of E.tenella.Methylated RNA immunoprecipitation sequencing(MeRIP-seq)analysis showed that m^(6)A modification was most abundant in the coding sequences,followed by stop codon.There were 3,903 hypermethylated and 3,178 hypomethylated mRNAs in sporulated oocysts compared with unsporulated oocysts.Further joint analysis suggested that m^(6)A modification of the majority of genes was positively correlated with mRNA expression.The mRNA relative expression and m^(6)A level of the selected genes were confirmed by quantitative reverse transcription PCR(RT-qPCR)and MeRIP-qPCR.GO and KEGG analysis indicated that differentially m^(6)A methylated genes(DMMGs)with significant differences in mRNA expression were closely related to processes such as regulation of gene expression,epigenetic,microtubule,autophagy-other and TOR signaling.Moreover,a total of 96 DMMGs without significant differences in mRNA expression showed significant differences at protein level.GO and pathway enrichment analysis of the 96 genes showed that RNA methylation may be involved in cell biosynthesis and metabolism of E.tenella.We firstly present a map of RNA m^(6)A modification in E.tenella,which provides significant insights into developmental biology of E.tenella.展开更多
Epigenetic changes in the spinal cord play a key role in the initiation and maintenance of nerve injury-induced neuro pathic pain.N6-methyladenosine(m6A)is one of the most abundant internal RNA modifications and plays...Epigenetic changes in the spinal cord play a key role in the initiation and maintenance of nerve injury-induced neuro pathic pain.N6-methyladenosine(m6A)is one of the most abundant internal RNA modifications and plays an essential function in gene regulation in many diseases.However,the global m6A modification status of mRNA in the spinal cord at different stages after neuropathic pain is unknown.In this study,we established a neuropathic pain model in mice by preserving the complete sural nerve and only damaging the common peroneal nerve.High-throughput methylated RNA immunoprecipitation sequencing res ults showed that after spared nerve injury,there were 55 m6A methylated and diffe rentially expressed genes in the spinal cord.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway results showed that m6A modification triggered inflammatory responses and apoptotic processes in the early stages after spared nerve injury.Over time,the diffe rential gene function at postoperative day 7 was enriched in "positive regulation of neurogenesis" and "positive regulation of neural precursor cell prolife ration." These functions suggested that altered synaptic morphological plasticity was a turning point in neuropathic pain formation and maintenance.Results at postoperative day 14 suggested that the persistence of neuropathic pain might be from lipid metabolic processes,such as "very-low-density lipoprotein particle clearance," "negative regulation of choleste rol transport" and "membrane lipid catabolic process." We detected the expression of m6A enzymes and found elevated mRNA expression of Ythdf2 and Ythdf3 after spared nerve injury modeling.We speculate that m6A reader enzymes also have an important role in neuropathic pain.These results provide a global landscape of mRNA m6A modifications in the spinal cord in the spared nerve injury model at diffe rent stages after injury.展开更多
BACKGROUND Severe acute pancreatitis(SAP)is a deadly inflammatory disease with complex pathogenesis and lack of effective therapeutic options.N6-methyladenosine(m6A)modification of circRNAs plays important roles in ph...BACKGROUND Severe acute pancreatitis(SAP)is a deadly inflammatory disease with complex pathogenesis and lack of effective therapeutic options.N6-methyladenosine(m6A)modification of circRNAs plays important roles in physiological and pathological processes.However,the roles of m6A circRNA in the pathological process of SAP remains unknown.AIM To identify transcriptome-wide map of m6A circRNAs and to determine their biological significance and potential mechanisms in SAP.METHODS The SAP in C57BL/6 mice was induced using 4%sodium taurocholate salt.The transcriptome-wide map of m6A circRNAs was identified by m6A-modified RNA immunoprecipitation sequencing.The biological significance of circRNAs with differentially expressed m6A peaks was evaluated through gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis.The underlying mechanism of m6A circRNAs in SAP was analyzed by constructing of m6A circRNAmicroRNA networks.The expression of demethylases was determined by quantitative polymerase chain reaction and western blot to deduce the possible mechanism of reversible m6A process in SAP.RESULTS Fifty-seven circRNAs with differentially expressed m6A peaks were identified by m6A-modified RNA immunoprecipitation sequencing,of which 32 were upregulated and 25 downregulated.Functional analysis of these m6A circRNAs in SAP found some important pathways involved in the pathogenesis of SAP,such as regulation of autophagy and protein digestion.In m6A circRNA–miRNA networks,several important miRNAs participated in the occurrence and progression of SAP were found to bind to these m6A circRNAs,such as miR-24-3p,miR-26a,miR-92b,miR-216b,miR-324-5p and miR-762.Notably,the total m6A level of circRNAs was reduced,while the demethylase alkylation repair homolog 5 was upregulated in SAP.CONCLUSION m6A modification of circRNAs may be involved in the pathogenesis of SAP.Our findings may provide novel insights to explore the possible pathogenetic mechanism of SAP and seek new potential therapeutic targets for SAP.展开更多
Background： The recently emerged technology of methylated RNA immunoprecipitation sequencing （MeRIP-seq） sheds light on the study of RNA epigenetics. This new bioinformatics question calls for effective and robust ...Background： The recently emerged technology of methylated RNA immunoprecipitation sequencing （MeRIP-seq） sheds light on the study of RNA epigenetics. This new bioinformatics question calls for effective and robust peaking calling algorithms to detect mRNA methylation sites from MeRIP-seq data. Methods： We propose a Bayesian hierarchical model to detect methylation sites from MeRIP-seq data. Our modeling approach includes several important characteristics. First, it models the zero-inflated and over-dispersed counts by deploying a zero-inflated negative binomial model. Second, it incorporates a hidden Markov model （HMM） to account for the spatial dependency of neighboring read enrichment. Third, our Bayesian inference allows the proposed model to borrow strength in parameter estimation, which greatly improves the model stability when dealing with MeRIP-seq data with a small number of replicates. We use Markov chain Monte Carlo （MCMC） algorithms to simultaneously infer the model parameters in a de novo fashion. The R Shiny demo is available at https：//qiwei. shinyapps.io/BaySeqPeak and the R/C ＋＋ code is available at https：//github.com/liqiwei2000/BaySeqPeak. Results： In simulation studies, the proposed method outperformed the competing methods exomePeak and MeTPeak, especially when an excess of zeros were present in the data. In real MeRIP-seq data analysis, the proposed method identified methylation sites that were more consistent with biological knowledge, and had better spatial resolution compared to the other methods. Conclusions： In this study, we develop a Bayesian hierarchical model to identify methylation peaks in MeRIP-seq data. The proposed method has a competitive edge over existing methods in terms of accuracy, robustness and spatial resolution.展开更多
基金supported by Grants from the Ministry of Science and Technology of China (Grant No. 2011CB944100) to JYthe Natural Science Foundation (Grant No. 30900831 and 31271372)Beijing Nova Program (Grant No. Z121105002512060)to SS
文摘RNA modifications, especially methylation of the N6 position of adenosine （A）--m6A, rep- resent an emerging research frontier in RNA biology. With the rapid development of high-throughput sequencing technology, in-depth study of m6A distribution and function relevance becomes feasible. However, a robust method to effectively identify m6A-modified regions has not been available yet. Here, we present a novel high-efficiency and user-friendly analysis pipeline called MeRIP-PF for the signal identification of MeRIP-Seq data in reference to controls. MeRIP-PF provides a statistical P-value for each identified m6A region based on the difference of read distribution when compared to the con- trols and also calculates false discovery rate （FDR） as a cut offto differentiate reliable m6A regions from the background. Furthermore, MeRIP-PF also achieves gene annotation ofm6A signals or peaks and produce outputs in both XLS and graphical format, which are useful for further study. MeRIP-PF is implemented in Perl and is freely available at http：//software.big.ac.cn/MeRIP-PF.html.
基金supported by the National Natural Science Foundation of China(31902298)the Shanxi Provincial Key Research and Development Program,China(2022ZDYF126)+2 种基金the Fund for Shanxi“1331 Project”,China(20211331-13)the Science and Technology Innovation Program of Shanxi Agricultural University,China(2017YJ10)the Special Research Fund of Shanxi Agricultural University for High-level Talents,China(2021XG001)。
文摘Evidence showed that N6-methyladenosine(m^(6)A)modification plays a pivotal role in influencing RNA fate and is strongly associated with cell growth and developmental processes in many species.However,no information regarding m^(6)A modification in Eimeria tenella is currently available.In the present study,we surveyed the transcriptome-wide prevalence of m^(6)A in sporulated oocysts and unsporulated oocysts of E.tenella.Methylated RNA immunoprecipitation sequencing(MeRIP-seq)analysis showed that m^(6)A modification was most abundant in the coding sequences,followed by stop codon.There were 3,903 hypermethylated and 3,178 hypomethylated mRNAs in sporulated oocysts compared with unsporulated oocysts.Further joint analysis suggested that m^(6)A modification of the majority of genes was positively correlated with mRNA expression.The mRNA relative expression and m^(6)A level of the selected genes were confirmed by quantitative reverse transcription PCR(RT-qPCR)and MeRIP-qPCR.GO and KEGG analysis indicated that differentially m^(6)A methylated genes(DMMGs)with significant differences in mRNA expression were closely related to processes such as regulation of gene expression,epigenetic,microtubule,autophagy-other and TOR signaling.Moreover,a total of 96 DMMGs without significant differences in mRNA expression showed significant differences at protein level.GO and pathway enrichment analysis of the 96 genes showed that RNA methylation may be involved in cell biosynthesis and metabolism of E.tenella.We firstly present a map of RNA m^(6)A modification in E.tenella,which provides significant insights into developmental biology of E.tenella.
基金National Natural Science Foundation of China,No.819 73305 (to ZQ)Science and Technology Planning Project of Guangzhou of China,No.20190401 0487 (to ZQ)+1 种基金Natural Science Foundation of Guangdong Province,China,No.2021A1515010897 (to TT)Discipline Construction Fund of Cen tral Peoples Hospital of Zhanjiang,Nos.2020A01 (to TT) and 2020A02 (to TT)。
文摘Epigenetic changes in the spinal cord play a key role in the initiation and maintenance of nerve injury-induced neuro pathic pain.N6-methyladenosine(m6A)is one of the most abundant internal RNA modifications and plays an essential function in gene regulation in many diseases.However,the global m6A modification status of mRNA in the spinal cord at different stages after neuropathic pain is unknown.In this study,we established a neuropathic pain model in mice by preserving the complete sural nerve and only damaging the common peroneal nerve.High-throughput methylated RNA immunoprecipitation sequencing res ults showed that after spared nerve injury,there were 55 m6A methylated and diffe rentially expressed genes in the spinal cord.Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway results showed that m6A modification triggered inflammatory responses and apoptotic processes in the early stages after spared nerve injury.Over time,the diffe rential gene function at postoperative day 7 was enriched in "positive regulation of neurogenesis" and "positive regulation of neural precursor cell prolife ration." These functions suggested that altered synaptic morphological plasticity was a turning point in neuropathic pain formation and maintenance.Results at postoperative day 14 suggested that the persistence of neuropathic pain might be from lipid metabolic processes,such as "very-low-density lipoprotein particle clearance," "negative regulation of choleste rol transport" and "membrane lipid catabolic process." We detected the expression of m6A enzymes and found elevated mRNA expression of Ythdf2 and Ythdf3 after spared nerve injury modeling.We speculate that m6A reader enzymes also have an important role in neuropathic pain.These results provide a global landscape of mRNA m6A modifications in the spinal cord in the spared nerve injury model at diffe rent stages after injury.
基金the National Natural Science Foundation of China,No.81772001the National Clinical Key Subject of China,No.41732113.
文摘BACKGROUND Severe acute pancreatitis(SAP)is a deadly inflammatory disease with complex pathogenesis and lack of effective therapeutic options.N6-methyladenosine(m6A)modification of circRNAs plays important roles in physiological and pathological processes.However,the roles of m6A circRNA in the pathological process of SAP remains unknown.AIM To identify transcriptome-wide map of m6A circRNAs and to determine their biological significance and potential mechanisms in SAP.METHODS The SAP in C57BL/6 mice was induced using 4%sodium taurocholate salt.The transcriptome-wide map of m6A circRNAs was identified by m6A-modified RNA immunoprecipitation sequencing.The biological significance of circRNAs with differentially expressed m6A peaks was evaluated through gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis.The underlying mechanism of m6A circRNAs in SAP was analyzed by constructing of m6A circRNAmicroRNA networks.The expression of demethylases was determined by quantitative polymerase chain reaction and western blot to deduce the possible mechanism of reversible m6A process in SAP.RESULTS Fifty-seven circRNAs with differentially expressed m6A peaks were identified by m6A-modified RNA immunoprecipitation sequencing,of which 32 were upregulated and 25 downregulated.Functional analysis of these m6A circRNAs in SAP found some important pathways involved in the pathogenesis of SAP,such as regulation of autophagy and protein digestion.In m6A circRNA–miRNA networks,several important miRNAs participated in the occurrence and progression of SAP were found to bind to these m6A circRNAs,such as miR-24-3p,miR-26a,miR-92b,miR-216b,miR-324-5p and miR-762.Notably,the total m6A level of circRNAs was reduced,while the demethylase alkylation repair homolog 5 was upregulated in SAP.CONCLUSION m6A modification of circRNAs may be involved in the pathogenesis of SAP.Our findings may provide novel insights to explore the possible pathogenetic mechanism of SAP and seek new potential therapeutic targets for SAP.
文摘Background： The recently emerged technology of methylated RNA immunoprecipitation sequencing （MeRIP-seq） sheds light on the study of RNA epigenetics. This new bioinformatics question calls for effective and robust peaking calling algorithms to detect mRNA methylation sites from MeRIP-seq data. Methods： We propose a Bayesian hierarchical model to detect methylation sites from MeRIP-seq data. Our modeling approach includes several important characteristics. First, it models the zero-inflated and over-dispersed counts by deploying a zero-inflated negative binomial model. Second, it incorporates a hidden Markov model （HMM） to account for the spatial dependency of neighboring read enrichment. Third, our Bayesian inference allows the proposed model to borrow strength in parameter estimation, which greatly improves the model stability when dealing with MeRIP-seq data with a small number of replicates. We use Markov chain Monte Carlo （MCMC） algorithms to simultaneously infer the model parameters in a de novo fashion. The R Shiny demo is available at https：//qiwei. shinyapps.io/BaySeqPeak and the R/C ＋＋ code is available at https：//github.com/liqiwei2000/BaySeqPeak. Results： In simulation studies, the proposed method outperformed the competing methods exomePeak and MeTPeak, especially when an excess of zeros were present in the data. In real MeRIP-seq data analysis, the proposed method identified methylation sites that were more consistent with biological knowledge, and had better spatial resolution compared to the other methods. Conclusions： In this study, we develop a Bayesian hierarchical model to identify methylation peaks in MeRIP-seq data. The proposed method has a competitive edge over existing methods in terms of accuracy, robustness and spatial resolution.
