m6A-Atlas v2.0: updated resources for unraveling the N6-methyladenosine (m6A) epitranscriptome among multiple species.

N 6-Methyladenosine (m6A) is one of the most abundant internal chemical modifications on eukaryote mRNA and is involved in numerous essential molecular functions and biological processes. To facilitate the study of this important post-transcriptional modification, we present here m6A-Atlas v2.0, an updated version of m6A-Atlas. It was expanded to include a total of 797 091 reliable m6A sites from 13 high-resolution technologies and two single-cell m6A profiles. Additionally, three methods (exomePeaks2, MACS2 and TRESS) were used to identify >16 million m6A enrichment peaks from 2712 MeRIP-seq experiments covering 651 conditions in 42 species. Quality control results of MeRIP-seq samples were also provided to help users to select reliable peaks. We also estimated the condition-specific quantitative m6A profiles (i.e. differential methylation) under 172 experimental conditions for 19 species. Further, to provide insights into potential functional circuitry, the m6A epitranscriptomics were annotated with various genomic features, interactions with RNA-binding proteins and microRNA, potentially linked splicing events and single nucleotide polymorphisms. The collected m6A sites and their functional annotations can be freely queried and downloaded via a user-friendly graphical interface at: http://rnamd.org/m6a.

m7GHub V2.0: an updated database for decoding the N7-methylguanosine (m7G) epitranscriptome.

With recent progress in mapping N7-methylguanosine (m7G) RNA methylation sites, tens of thousands of experimentally validated m7G sites have been discovered in various species, shedding light on the significant role of m7G modification in regulating numerous biological processes including disease pathogenesis. An integrated resource that enables the sharing, annotation and customized analysis of m7G data will greatly facilitate m7G studies under various physiological contexts. We previously developed the m7GHub database to host mRNA m7G sites identified in the human transcriptome. Here, we present m7GHub v.2.0, an updated resource for a comprehensive collection of m7G modifications in various types of RNA across multiple species: an m7GDB database containing 430 898 putative m7G sites identified in 23 species, collected from both widely applied next-generation sequencing (NGS) and the emerging Oxford Nanopore direct RNA sequencing (ONT) techniques; an m7GDiseaseDB hosting 156 206 m7G-associated variants (involving addition or removal of an m7G site), including 3238 disease-relevant m7G-SNPs that may function through epitranscriptome disturbance; and two enhanced analysis modules to perform interactive analyses on the collections of m7G sites (m7GFinder) and functional variants (m7GSNPer). We expect that m7Ghub v.2.0 should serve as a valuable centralized resource for studying m7G modification. It is freely accessible at: www.rnamd.org/m7GHub2.

DirectRMDB: a database of post-transcriptional RNA modifications unveiled from direct RNA sequencing technology.

With advanced technologies to map RNA modifications, our understanding of them has been revolutionized, and they are seen to be far more widespread and important than previously thought. Current next-generation sequencing (NGS)-based modification profiling methods are blind to RNA modifications and thus require selective chemical treatment or antibody immunoprecipitation methods for particular modification types. They also face the problem of short read length, isoform ambiguities, biases and artifacts. Direct RNA sequencing (DRS) technologies, commercialized by Oxford Nanopore Technologies (ONT), enable the direct interrogation of any given modification present in individual transcripts and promise to address the limitations of previous NGS-based methods. Here, we present the first ONT-based database of quantitative RNA modification profiles, DirectRMDB, which includes 16 types of modification and a total of 904,712 modification sites in 25 species identified from 39 independent studies. In addition to standard functions adopted by existing databases, such as gene annotations and post-transcriptional association analysis, we provide a fresh view of RNA modifications, which enables exploration of the epitranscriptome in an isoform-specific manner. The DirectRMDB database is freely available at: http://www.rnamd.org/directRMDB/.

RMDisease V2.0: an updated database of genetic variants that affect RNA modifications with disease and trait implication.

Recent advances in epitranscriptomics have unveiled functional associations between RNA modifications (RMs) and multiple human diseases, but distinguishing the functional or disease-related single nucleotide variants (SNVs) from the majority of 'silent' variants remains a major challenge. We previously developed the RMDisease database for unveiling the association between genetic variants and RMs concerning human disease pathogenesis. In this work, we present RMDisease v2.0, an updated database with expanded coverage. Using deep learning models and from 873 819 experimentally validated RM sites, we identified a total of 1 366 252 RM-associated variants that may affect (add or remove an RM site) 16 different types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G, A-to-I, ac4C, Am, Cm, Um, Gm, hm5C, D and f5C) in 20 organisms (human, mouse, rat, zebrafish, maize, fruit fly, yeast, fission yeast, Arabidopsis, rice, chicken, goat, sheep, pig, cow, rhesus monkey, tomato, chimpanzee, green monkey and SARS-CoV-2). Among them, 14 749 disease- and 2441 trait-associated genetic variants may function via the perturbation of epitranscriptomic markers. RMDisease v2.0 should serve as a useful resource for studying the genetic drivers of phenotypes that lie within the epitranscriptome layer circuitry, and is freely accessible at: www.rnamd.org/rmdisease2.

m5C-Atlas: a comprehensive database for decoding and annotating the 5-methylcytosine (m5C) epitranscriptome.

5-Methylcytosine (m5C) is one of the most prevalent covalent modifications on RNA. It is known to regulate a broad variety of RNA functions, including nuclear export, RNA stability and translation. Here, we present m5C-Atlas, a database for comprehensive collection and annotation of RNA 5-methylcytosine. The database contains 166 540 m5C sites in 13 species identified from 5 base-resolution epitranscriptome profiling technologies. Moreover, condition-specific methylation levels are quantified from 351 RNA bisulfite sequencing samples gathered from 22 different studies via an integrative pipeline. The database also presents several novel features, such as the evolutionary conservation of a m5C locus, its association with SNPs, and any relevance to RNA secondary structure. All m5C-atlas data are accessible through a user-friendly interface, in which the m5C epitranscriptomes can be freely explored, shared, and annotated with putative post-transcriptional mechanisms (e.g. RBP intermolecular interaction with RNA, microRNA interaction and splicing sites). Together, these resources offer unprecedented opportunities for exploring m5C epitranscriptomes. The m5C-Atlas database is freely accessible at https://www.xjtlu.edu.cn/biologicalsciences/m5c-atlas.

Geographic encoding of transcripts enabled high-accuracy and isoform-aware deep learning of RNA methylation.

As the most pervasive epigenetic mark present on mRNA and lncRNA, N6-methyladenosine (m6A) RNA methylation regulates all stages of RNA life in various biological processes and disease mechanisms. Computational methods for deciphering RNA modification have achieved great success in recent years; nevertheless, their potential remains underexploited. One reason for this is that existing models usually consider only the sequence of transcripts, ignoring the various regions (or geography) of transcripts such as 3'UTR and intron, where the epigenetic mark forms and functions. Here, we developed three simple yet powerful encoding schemes for transcripts to capture the submolecular geographic information of RNA, which is largely independent from sequences. We show that m6A prediction models based on geographic information alone can achieve comparable performances to classic sequence-based methods. Importantly, geographic information substantially enhances the accuracy of sequence-based models, enables isoform- and tissue-specific prediction of m6A sites, and improves m6A signal detection from direct RNA sequencing data. The geographic encoding schemes we developed have exhibited strong interpretability, and are applicable to not only m6A but also N1-methyladenosine (m1A), and can serve as a general and effective complement to the widely used sequence encoding schemes in deep learning applications concerning RNA transcripts.

ConsRM: collection and large-scale prediction of the evolutionarily conserved RNA methylation sites, with implications for the functional epitranscriptome.

Motivation N6-methyladenosine (m6A) is the most prevalent RNA modification on mRNAs and lncRNAs. Evidence increasingly demonstrates its crucial importance in essential molecular mechanisms and various diseases. With recent advances in sequencing techniques, tens of thousands of m6A sites are identified in a typical high-throughput experiment, posing a key challenge to distinguish the functional m6A sites from the remaining 'passenger' (or 'silent') sites. Results: We performed a comparative conservation analysis of the human and mouse m6A epitranscriptomes at single site resolution. A novel scoring framework, ConsRM, was devised to quantitatively measure the degree of conservation of individual m6A sites. ConsRM integrates multiple information sources and a positive-unlabeled learning framework, which integrated genomic and sequence features to trace subtle hints of epitranscriptome layer conservation. With a series validation experiments in mouse, fly and zebrafish, we showed that ConsRM outperformed well-adopted conservation scores (phastCons and phyloP) in distinguishing the conserved and unconserved m6A sites. Additionally, the m6A sites with a higher ConsRM score are more likely to be functionally important. An online database was developed containing the conservation metrics of 177 998 distinct human m6A sites to support conservation analysis and functional prioritization of individual m6A sites. And it is freely accessible at: https://www.xjtlu.edu.cn/biologicalsciences/con.

Primary sequence-assisted prediction of m6A RNA methylation sites from Oxford nanopore direct RNA sequencing data.

Traditional epitranscriptome profiling approach relies on specific antibodies or chemical treatments to capture modified RNA molecules and then applies high throughput sequencing to identify their transcriptomic locations. However, due to the lack of suitable or high-quality antibodies, only a small proportion of the 170 documented RNA modifications were profiled with those approaches. Direct sequencing of native RNA molecules using Oxford Nanopore Technologies (ONT) enabled straight inspection of RNA modifications and offered a promising alternative solution. N6-methyladenosine (m6A) is known to cause characteristic changes and increased base call errors of ONT signals compared with non-modified adenosines, based on which, the m6A sites can be identified directly from ONT signals. Meanwhile, a number of studies have shown that it is possible to predict m6A sites from RNA primary sequences. Using the m6A sites revealed by Illumina technology as benchmark, we showed that, the accuracy of ONT-based m6A site prediction can be further increased by integrating additional information from the primary sequences of RNA (AUROC of 0.918), compared with using ONT signals only (AUROC 0.878 using Base call error features, and 0.804 using Tombo features), providing a new perspective for more reliable mining of the relatively noisy ONT signals.

m6AmPred: Identifying RNA N6, 2'-O-dimethyladenosine (m6Am) sites based on sequence-derived information.

N6,2'-O-dimethyladenosine (m6Am) is a reversible modification widely occurred on varied RNA molecules. The biological function of m6Am is yet to be known though recent studies have revealed its influences in cellular mRNA fate. Precise identification of m6Am sites on RNA is vital for the understanding of its biological functions. We present here m6AmPred, the first web server for in silico identification of m6Am sites from the primary sequences of RNA. Built upon the eXtreme Gradient Boosting with Dart algorithm (XgbDart) and EIIP-PseEIIP encoding scheme, m6AmPred achieved promising prediction performance with the AUCs greater than 0.954 when tested by 10-fold cross-validation and independent testing datasets. To critically test and validate the performance of m6AmPred, the experimentally verified m6Am sites from two data sources were cross-validated. The m6AmPred web server is freely accessible at: https://www.xjtlu.edu.cn/biologicalsciences/m6am, and it should make a useful tool for the researchers who are interested in N6,2'-O-dimethyladenosine RNA modification.

WHISTLE server: A high-accuracy genomic coordinate-based machine learning platform for RNA modification prediction.

The primary sequences of DNA, RNA and protein have been used as the dominant information source of existing machine learning tools, especially for contexts not fully explored by wet-experimental approaches. Since molecular markers are profoundly orchestrated in the living organisms, those markers that cannot be unambiguously recovered from the primary sequence often help to predict other biological events. To the best of our knowledge, there is no current tool to build and deploy machine learning models that consider genomic evidence. We therefore developed the WHISTLE server, the first machine learning platform based on genomic coordinates. It features convenient covariate extraction and model web deployment with 46 distinct genomic features integrated along with the conventional sequence features. We showed that, when predicting m6A sites from SRAMP project, the model integrating genomic features substantially outperformed those based on only sequence features. The WHISTLE server should be a useful tool for studying biological attributes specifically associated with genomic coordinates, and is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/whi2.

RMDisease: a database of genetic variants that affect RNA modifications, with implications for epitranscriptome pathogenesis.

Deciphering the biological impacts of millions of single nucleotide variants remains a major challenge. Recent studies suggest that RNA modifications play versatile roles in essential biological mechanisms, and are closely related to the progression of various diseases including multiple cancers. To comprehensively unveil the association between disease-associated variants and their epitranscriptome disturbance, we built RMDisease, a database of genetic variants that can affect RNA modifications. By integrating the prediction results of 18 different RNA modification prediction tools and also 303,426 experimentally-validated RNA modification sites, RMDisease identified a total of 202,307 human SNPs that may affect (add or remove) sites of eight types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G and Nm). These include 4,289 disease-associated variants that may imply disease pathogenesis functioning at the epitranscriptome layer. These SNPs were further annotated with essential information such as post-transcriptional regulations (sites for miRNA binding, interaction with RNA-binding proteins and alternative splicing) revealing putative regulatory circuits. A convenient graphical user interface was constructed to support the query, exploration and download of the relevant information. RMDisease should make a useful resource for studying the epitranscriptome impact of genetic variants via multiple RNA modifications with emphasis on their potential disease relevance. RMDisease is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/rmd.

m6A-Atlas: a comprehensive knowledgebase for unraveling the N6-methyladenosine (m6A) epitranscriptome.

N 6-Methyladenosine (m6A) is the most prevalent RNA modification on mRNAs and lncRNAs. It plays a pivotal role during various biological processes and disease pathogenesis. We present here a comprehensive knowledgebase, m6A-Atlas, for unraveling the m6A epitranscriptome. Compared to existing databases, m6A-Atlas features a high-confidence collection of 442 162 reliable m6A sites identified from seven base-resolution technologies and the quantitative (rather than binary) epitranscriptome profiles estimated from 1363 high-throughput sequencing samples. It also offers novel features, such as; the conservation of m6A sites among seven vertebrate species (including human, mouse and chimp), the m6A epitranscriptomes of 10 virus species (including HIV, KSHV and DENV), the putative biological functions of individual m6A sites predicted from epitranscriptome data, and the potential pathogenesis of m6A sites inferred from disease-associated genetic mutations that can directly destroy m6A directing sequence motifs. A user-friendly graphical user interface was constructed to support the query, visualization and sharing of the m6A epitranscriptomes annotated with sites specifying their interaction with post-transcriptional machinery (RBP-binding, microRNA interaction and splicing sites) and interactively display the landscape of multiple RNA modifications. These resources provide fresh opportunities for unraveling the m6A epitranscriptomes. m6A-Atlas is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/atlas.

MetaTX: deciphering the distribution of mRNA-related features in the presence of isoform ambiguity, with applications in epitranscriptome analysis.

MOTIVATION: The distribution of biological features strongly indicates their functional relevance. Compared to DNA-related features, deciphering the distribution of mRNA-related features is non-trivial due to the existence of isoform ambiguity and compositional diversity of mRNAs. RESULTS: We propose here a rigorous statistical framework, MetaTX, for deciphering the distribution of mRNA-related features. Through a standardized mRNA model, MetaTX firstly unifies various mRNA transcripts of diverse compositions, and then corrects the isoform ambiguity by incorporating the overall distribution pattern of the features through an EM algorithm. MetaTX was tested on both simulated and real data. Results suggested that MetaTX substantially outperformed existing direct methods on simulated datasets, and that a more informative distribution pattern was produced for all the three datasets tested, which contain N6-Methyladenosine sites generated by different technologies. MetaTX should make a useful tool for studying the distribution and functions of mRNA-related biological features, especially for mRNA modifications such as N6-Methyladenosine. AVAILABILITY AND IMPLEMENTATION: The MetaTX R package is freely available at GitHub: https://github.com/yue-wang-biomath/MetaTX.1.0. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

m7GHub: deciphering the location, regulation and pathogenesis of internal mRNA N7-methylguanosine (m7G) sites in human.

MOTIVATION: Recent progress in N7-methylguanosine (m7G) RNA methylation studies has focused on its internal (rather than capped) presence within mRNAs. Tens of thousands of internal mRNA m7G sites have been identified within mammalian transcriptomes, and a single resource to best share, annotate and analyze the massive m7G data generated recently are sorely needed. RESULTS: We report here m7GHub, a comprehensive online platform for deciphering the location, regulation and pathogenesis of internal mRNA m7G. The m7GHub consists of four main components, including: the first internal mRNA m7G database containing 44 058 experimentally validated internal mRNA m7G sites, a sequence-based high-accuracy predictor, the first web server for assessing the impact of mutations on m7G status, and the first database recording 1218 disease-associated genetic mutations that may function through regulation of m7G methylation. Together, m7GHub will serve as a useful resource for research on internal mRNA m7G modification. AVAILABILITY AND IMPLEMENTATION: m7GHub is freely accessible online at www.xjtlu.edu.cn/biologicalsciences/m7ghub. CONTACT: kunqi.chen@liverpool.ac.uk. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

WHISTLE: a high-accuracy map of the human N6-methyladenosine (m6A) epitranscriptome predicted using a machine learning approach.

N 6-methyladenosine (m6A) is the most prevalent post-transcriptional modification in eukaryotes, and plays a pivotal role in various biological processes, such as splicing, RNA degradation and RNA-protein interaction. We report here a prediction framework WHISTLE for transcriptome-wide m6A RNA-methylation site prediction. When tested on six independent datasets, our approach, which integrated 35 additional genomic features besides the conventional sequence features, achieved a major improvement in the accuracy of m6A site prediction (average AUC: 0.948 and 0.880 under the full transcript or mature messenger RNA models, respectively) compared to the state-of-the-art computational approaches MethyRNA (AUC: 0.790 and 0.732) and SRAMP (AUC: 0.761 and 0.706). It also out-performed the existing epitranscriptome databases MeT-DB (AUC: 0.798 and 0.744) and RMBase (AUC: 0.786 and 0.736), which were built upon hundreds of epitranscriptome high-throughput sequencing samples. To probe the putative biological processes impacted by changes in an individual m6A site, a network-based approach was implemented according to the 'guilt-by-association' principle by integrating RNA methylation profiles, gene expression profiles and protein-protein interaction data. Finally, the WHISTLE web server was built to facilitate the query of our high-accuracy map of the human m6A epitranscriptome, and the server is freely available at: www.xjtlu.edu.cn/biologicalsciences/whistle and http://whistle-epitranscriptome.com.

m6Acomet: large-scale functional prediction of individual m6A RNA methylation sites from an RNA co-methylation network.

BACKGROUND: Over one hundred different types of post-transcriptional RNA modifications have been identified in human. Researchers discovered that RNA modifications can regulate various biological processes, and RNA methylation, especially N6-methyladenosine, has become one of the most researched topics in epigenetics. RESULTS: To date, the study of epitranscriptome layer gene regulation is mostly focused on the function of mediator proteins of RNA methylation, i.e., the readers, writers and erasers. There is limited investigation of the functional relevance of individual m6A RNA methylation site. To address this, we annotated human m6A sites in large-scale based on the guilt-by-association principle from an RNA co-methylation network. It is constructed based on public human MeRIP-Seq datasets profiling the m6A epitranscriptome under 32 independent experimental conditions. By systematically examining the network characteristics obtained from the RNA methylation profiles, a total of 339,158 putative gene ontology functions associated with 1446 human m6A sites were identified. These are biological functions that may be regulated at epitranscriptome layer via reversible m6A RNA methylation. The results were further validated on a soft benchmark by comparing to a random predictor. CONCLUSIONS: An online web server m6Acomet was constructed to support direct query for the predicted biological functions of m6A sites as well as the sites exhibiting co-methylated patterns at the epitranscriptome layer. The m6Acomet web server is freely available at: www.xjtlu.edu.cn/biologicalsciences/m6acomet .

Construction and validation of a prognostic signature for mucinous colonic adenocarcinoma based on N7-methylguanosine-related long non-coding RNAs.

Background: Mucinous colonic adenocarcinoma remains a challenging disease due to its high propensity for metastasis and recurrence. N7-methylguanosine (m7G) and long non-coding RNA (lncRNA) are closely associated with the occurrence and progression of tumors. However, research on m7G-related lncRNA in mucinous colonic adenocarcinoma is lacking. Therefore, we sought to explore the prognostic impact of m7G-related lncRNAs in mucinous adenocarcinoma (MC) patients. Methods: In this study, Pearson analysis was used to identify m7G-related lncRNAs from transcriptome data in The Cancer Genome Atlas (TCGA). Univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression were used to further screen m7G-related lncRNAs and incorporate them into a prognostic signature. Based on the risk model, patients were divided into low- and high-risk groups and randomly assigned to the training set and test sets in a 6:4 ratio. Kaplan-Meier, receiver operating characteristic (ROC) curve, multivariate regression, and nomogram analyses were used to confirm the accuracy of the signature. The CIBERSORT algorithm was used to calculate the degree of immune cell infiltration (ICI). Finally, the correlation of the prognostic signature with tumor mutational burden (TMB) and immunophenotype score (IPS) was evaluated. Results: A total of 432 m7G-related lncRNAs were identified by Pearson analysis. Univariate Cox regression, LASSO regression and survival analysis were performed to further select six m7G-related lncRNAs (P<0.05): AC254629.1, LINC01133, LINC01134, MHENCR, SMIM2-AS1, and XACT. Based on the risk model, heat maps, Kaplan-Meier curves, and ROC curves were constructed, and the results showed that there were significant differences in expression levels and survival status between the two risk groups. The area under the ROC curve (AUC) values for 3-, 5-, and 10-year survival in the training set were 0.944, 0.957, and 1.000, respectively. And in the test set were 0.964, 1.000, and 1.000, respectively. Subsequently, univariate and multivariate regression analyses of clinical characteristics and risk score were performed. The results of risk score were [hazard ratio (HR): 6.458, 95% confidence interval (CI): 2.708-15.403, P<0.001; HR: 7.280, 95% CI: 2.500-21.203, P<0.001], respectively. Using the risk score as an independent prognostic factor, the AUC of it over 3, 5, and 10 years was 0.911, 0.955, and 0.961, respectively. Calibration plots for the nomogram show that the model calibration line is very close to the ideal calibration line, indicating good calibration. The level of ICI was significantly different in the different risk groups. Survival analysis showed that, regardless of TMB risk, patients with MC and a high-risk score consistently had a poor overall survival (OS). Conclusions: The m7G-related lncRNA prognostic signature has potential value for the prognosis of mucinous colonic adenocarcinoma.

Unveiling the Impact of ApoF Deficiency on Liver and Lipid Metabolism: Insights from Transcriptome-Wide m6A Methylome Analysis in Mice.

Lipid metabolism participates in various physiological processes and has been shown to be connected to the development and progression of multiple diseases, especially metabolic hepatopathy. Apolipoproteins (Apos) act as vectors that combine with lipids, such as cholesterol and triglycerides (TGs). Despite being involved in lipid transportation and metabolism, the critical role of Apos in the maintenance of lipid metabolism has still not been fully revealed. This study sought to clarify variations related to m6A methylome in ApoF gene knockout mice with disordered lipid metabolism based on the bioinformatics method of transcriptome-wide m6A methylome epitranscriptomics. High-throughput methylated RNA immunoprecipitation sequencing (MeRIP-seq) was conducted in both wild-type (WT) and ApoF knockout (KO) mice. As a result, the liver histopathology presented vacuolization and steatosis, and the serum biochemical assays reported abnormal lipid content in KO mice. The m6A-modified mRNAs were conformed consensus sequenced in eukaryotes, and the distribution was enriched within the coding sequences and 3' non-coding regions. In KO mice, the functional annotation terms of the differentially expressed genes (DEGs) included cholesterol, steroid and lipid metabolism, and lipid storage. In the differentially m6A-methylated mRNAs, the functional annotation terms included cholesterol, TG, and long-chain fatty acid metabolic processes; lipid transport; and liver development. The overlapping DEGs and differential m6A-modified mRNAs were also enriched in terms of lipid metabolism disorder. In conclusion, transcriptome-wide MeRIP sequencing in ApoF KO mice demonstrated the role of this crucial apolipoprotein in liver health and lipid metabolism.

Involvement of autophagy in mesaconitine-induced neurotoxicity in HT22 cells revealed through integrated transcriptomic, proteomic, and m6A epitranscriptomic profiling.

Background: Mesaconitine (MA), a diester-diterpenoid alkaloid extracted from the medicinal herb Aconitum carmichaelii, is commonly used to treat various diseases. Previous studies have indicated the potent toxicity of aconitum despite its pharmacological activities, with limited understanding of its effects on the nervous system and the underlying mechanisms. Methods: HT22 cells and zebrafish were used to investigate the neurotoxic effects of MA both in vitro and in vivo, employing multi-omics techniques to explore the potential mechanisms of toxicity. Results: Our results demonstrated that treatment with MA induces neurotoxicity in zebrafish and HT22 cells. Subsequent analysis revealed that MA induced oxidative stress, as well as structural and functional damage to mitochondria in HT22 cells, accompanied by an upregulation of mRNA and protein expression related to autophagic and lysosomal pathways. Furthermore, methylated RNA immunoprecipitation sequencing (MeRIP-seq) showed a correlation between the expression of autophagy-related genes and N6-methyladenosine (m6A) modification following MA treatment. In addition, we identified METTL14 as a potential regulator of m6A methylation in HT22 cells after exposure to MA. Conclusion: Our study has contributed to a thorough mechanistic elucidation of the neurotoxic effects caused by MA, and has provided valuable insights for optimizing the rational utilization of traditional Chinese medicine formulations containing aconitum in clinical practice.

Gelsenicine disrupted the intestinal barrier of Caenorhabditis elegans.

Gelsemium elegans Benth. (G. elegans) is a traditional medicinal herb that has anti-inflammatory, analgesic, sedative, and detumescence effects. However, it can also cause intestinal side effects such as abdominal pain and diarrhea. The toxicological mechanisms of gelsenicine are still unclear. The objective of this study was to assess enterotoxicity induced by gelsenicine in the nematodes Caenorhabditis elegans (C. elegans). The nematodes were treated with gelsenicine, and subsequently their growth, development, and locomotion behavior were evaluated. The targets of gelsenicine were predicted using PharmMapper. mRNA-seq was performed to verify the predicted targets. Intestinal permeability, ROS generation, and lipofuscin accumulation were measured. Additionally, the fluorescence intensities of GFP-labeled proteins involved in oxidative stress and unfolded protein response in endoplasmic reticulum (UPRER) were quantified. As a result, the treatment of gelsenicine resulted in the inhibition of nematode lifespan, as well as reductions in body length, width, and locomotion behavior. A total of 221 targets were predicted by PharmMapper, and 731 differentially expressed genes were screened out by mRNA-seq. GO and KEGG enrichment analysis revealed involvement in redox process and transmembrane transport. The permeability assay showed leakage of blue dye from the intestinal lumen into the body cavity. Abnormal mRNAs expression of gem-4, hmp-1, fil-2, and pho-1, which regulated intestinal development, absorption and catabolism, transmembrane transport, and apical junctions, was observed. Intestinal lipofuscin and ROS were increased, while sod-2 and isp-1 expressions were decreased. Multiple proteins in SKN-1/DAF-16 pathway were found to bind stably with gelsenicine in a predictive model. There was an up-regulation in the expression of SKN-1:GFP, while the nuclear translocation of DAF-16:GFP exhibited abnormality. The UPRER biomarker HSP-4:GFP was down-regulated. In conclusion, the treatment of gelsenicine resulted in the increase of nematode intestinal permeability. The toxicological mechanisms underlying this effect involved the disruption of intestinal barrier integrity, an imbalance between oxidative and antioxidant processes mediated by the SKN-1/DAF-16 pathway, and abnormal unfolded protein reaction.