BRIEF RESEARCH REPORT article

Front. Genet., 21 August 2023

Sec. Computational Genomics

Volume 14 - 2023 | https://doi.org/10.3389/fgene.2023.1254827

Comparative evaluation and analysis of DNA N4-methylcytosine methylation sites using deep learning

  • 1. Heilongjiang Agricultural Engineering Vocational College, Harbin, China

  • 2. Engineering Research Center of Integration and Application of Digital Learning Technology, Ministry of Education, Hangzhou, China

  • 3. Beidahuang Industry Group General Hospital, Harbin, China

  • 4. Department of Neurosurgical Laboratory, The First Affiliated Hospital of Harbin Medical University, Harbin, China

Abstract

DNA N4-methylcytosine (4mC) is significantly involved in biological processes, such as DNA expression, repair, and replication. Therefore, accurate prediction methods are urgently needed. Deep learning methods have transformed applications that previously require sequencing expertise into engineering challenges that do not require expertise to solve. Here, we compare a variety of state-of-the-art deep learning models on six benchmark datasets to evaluate their performance in 4mC methylation site detection. We visualize the statistical analysis of the datasets and the performance of different deep-learning models. We conclude that deep learning can greatly expand the potential of methylation site prediction.

Introduction

The rapid progress in genome sequencing technologies has facilitated the investigation of the functional effects of DNA chemical modifications with unprecedented precision (; ; ). DNA methylation, as a vital epigenetic modification, plays a crucial role in normal organism development and essential biological processes (). In the genomes of both prokaryotic and eukaryotic organisms, the most prevalent kinds of DNA methylation include N6-methyladenine (6mA) (; ; ), C5-methylcytosine (5mC) (), and N4-methylcytosine (4mC) (; ; ; ; ). The distribution of 4mC sites in the genome is highly significant as they play a crucial role in regulating gene expression and maintaining genome stability. Accurate identification and analysis of 4mC sites allow for a deeper understanding of the role of DNA methylation in gene regulation and disease mechanisms. This has important implications for the study of epigenetics, cancer etiology, biological evolution, and potential therapeutic strategies. Therefore, the development of efficient and accurate methods for detecting and identifying 4mC sites is of great importance for understanding biological processes and disease research (; ).

Several experimental techniques have been utilized to identify epigenetic 4mC sites. These methodologies include methylation-specific PCR, mass spectrometry, 4mC-Tet-assisted bisulfite-sequencing (4mCTABseq), whole-genome bisulfite sequencing, nanopore sequencing, and single-molecule real-time (SMRT) sequencing (; ; ; ; ). These experiment-based methods suffer from limitations such as low throughput, high cost, and restricted detection sensitivity. Nowadays, machine learning has been widely utilized and are successful technology in bioinformatics for extracting knowledge from huge data (; ; ; ; ; ; ; ; ) and numerous computer techniques have been created to anticipate DNA 4mC sites. Both standard machine learning techniques and more current deep learning algorithms have been used to provide a strong result. In the field of 4mC site prediction, researchers have made significant strides by leveraging machine learning algorithms. These approaches utilize computational models to identify and classify 4mC sites within DNA sequences. Various machine learning techniques have been explored, including support vector machine (SVM) (), random forest (RF), Markov model (MM), and ensemble methods. Additionally, advanced techniques such as extreme gradient boosting (XGBoost) and Laplacian Regularized Sparse Representation have also been employed in this context (; ; ; ; ; ; ). However, traditional machine learning algorithms rely significantly on data representations known as features for appropriate performance, and it's tough to figure out which features are best for a certain task. Deep learning overcomes the limitations of traditional methods by offering adaptivity, fault tolerance, nonlinearity, and improved input-to-output mapping. Deep learning methods, such as convolutional neural networks (CNNs) and recurrent neural networks (RNNs), have been developed for the detection of 4mC sites, leveraging their ability to capture sequence patterns and dependencies, thereby contributing to accurate identification of these sites and enhancing our understanding of DNA methylation in gene regulation and epigenetics (; ). Yet there are still many deep learning methods that have not been applied, which have achieved great success in various application scenarios, including computer vision, speech recognition, biomarker identification (; ) and drug discovery (; ; ; ; ; ).

Choosing an appropriate deep learning model for bioinformatics problems poses a significant challenge for biologists. Understanding and comparing the performance of different models on specific datasets is of paramount importance for guiding practical applications. Therefore, our research focuses on evaluating the performance of multiple deep learning models on the 4mC datasets, aiming to assist biologists in making informed decisions when selecting suitable models.

We selected several common deep learning models, including RNN (Recurrent Nerual Network) (), long short-term memory (LSTM) (), bi-directional long short-term memory (Bi-LSTM) (; ), text convolutional neural network (Text-CNN) (), and bidirectional encoder representations from transformers (BERT) (; ), and compared their performances on the 4mC datasets through optimization of model hyperparameters. Our research findings provide strong evidence-based support for biologists, aiding them in making informed choices when addressing bioinformatics problems on the 4mC datasets. By comparing the performance of multiple models, we can offer recommendations tailored to different problems and datasets, enabling biologists to better understand and leverage the advantages of deep learning models.

Materials and methods

The implementation of our experiments relies on the DeepBIO () platform, which provides a wide selection of deep learning models and a visual comparison of multiple models. Figure 1 illustrates the overall framework of our works. We selected four deep learning models (RNN, LSTM, Bi-LSTM, Text-CNN) and pre-trained BERT models from the DeepBIO platform, and BERT is used as our main method to compare with other methods.

FIGURE 1

Datasets

The first step in creating a strong and trustworthy classification model is creating high-quality benchmark datasets. In this study, six benchmark datasets were utilized (). Table 1 provides a statistical summary of the datasets. The positive samples consisted of sequences that were 41 base pairs (bp) in length and contained a 4mC (4-methylcytosine) site located in the middle. These datasets have undergone rigorous preprocessing and quality control measures to ensure data accuracy and consistency (). By training and evaluating the model on data from multiple species, including humans, animals, and plants, we ensure its broad applicability and provide valuable insights for biologists in selecting deep learning models.

TABLE 1

SpeciesPositive sampleNegative sampleTotal
C. elegans1,5541,5543,108
D. melanogaster1,7691,7693,538
A. thaliana1,9781,9783,956
E. coli388388776
G. subterraneus9069061,812
G. pickeringii5695691,138

Statistical summary of benchmark datasets.

Input feature matrix

Deep learning algorithms possess the capability to autonomously extract valuable features from data, distinguishing them from conventional machine learning methods that necessitate manual feature engineering. Nonetheless, when dealing with a string of nucleotide letters (A, C, G, and T), a conversion into a matrix format is required prior to feeding it into a neural network layer. Unlike prior methods that used several feature encodings schemes to represent the sequence as the input to train the model, this method uses a single feature encoding scheme. We took the dictionary encoding approaches for representing DNA sequences. To represent DNA sequences, we utilized a dictionary encoding method where each nucleotide (A, C, G, and T) is assigned a numeric value. Specifically, A is represented by 1, C by 2, G by 3, and T by 4. This encoding scheme allows us to convert the sequence into an N-dimensional vector, facilitating its input into the neural network for further analysis.

Model construction and parameters

We have selected deep-learning models that have received a lot of attention in recent years as follows: RNN, LSTM, Bi-LSTM, Text-CNN, and BERT. The first four deep learning models we used are the models provided by the DeepBIO platform with parameters already set and the BERT model we used is pre-trained DNABERT (; ), which achieves the best performance on several DNA sequence classification tasks.

RNN is a type of neural network where the output of the previous neuron is fed back as input to the current neuron, creating temporal memory and enabling the processing of dynamic input sequences. RNNs find wide applications in various domains, including voice recognition, time series analysis, DNA sequences, and sequential data processing. One notable variant of RNNs that addresses the issue of capturing long-term dependencies is Long Short-Term Memory (LSTM). LSTM introduces a cell state that serves as a memory component, allowing the network to retain relevant information over extended periods. The forget gate in LSTM controls which information should be discarded and retained by using a sigmoid activation function. Additionally, Bidirectional LSTM (BiLSTM) processes input data in both forward and backward directions, effectively incorporating information from both past and future states. This bidirectional approach enables BiLSTM to capture intricate sequential relationships between words and sentences, making it particularly advantageous for Natural Language Processing (NLP) tasks that require contextual information from both preceding and succeeding elements in the input sequence. The RNN, LSTM, and Bi-LSTM architectures consist of stacked RNN cells, LSTM cells, and bidirectional LSTM cells, respectively. All these architectures share a similar structure, featuring 128 hidden neurons and a single layer for optimal performance. To prevent overfitting and promote generalization, a dropout rate of 0.2 was applied, and the output layer utilized sigmoid activation with a single neuron.

Text-CNN, a powerful deep learning approach for language classification tasks, such as sentiment analysis and question categorization, is a convolutional neural network tailored for text processing. The core structure comprises four layers: an embedding layer, a convolution layer, a pooling layer, and a fully connected layer. In our implementation, we set four convolutional kernel sizes (1, 2, 4, 8), and the number of convolutional kernels is uniformly set to 128. The embeddings undergo convolutional operations with a sliding kernel, producing convolutions that are subsequently downsampled through a Max Pooling layer to manage complexity and computational requirements. The scalar pooling outputs are then concatenated to form a vector representation of the input sequence. To mitigate overfitting, regularization methods, including a dropout layer with a rate of 0.2 and ReLU activation, are employed in the penultimate layer, preventing overfitting of the hidden layer.

BERT, an abbreviation for Bidirectional Encoder Representations from Transformers, originates from the Transformer architecture. In the Transformer model, every output element is intricately connected to every input element, with dynamically calculated weightings based on their connections. BERT is a pre-trained model that benefits from its ability to learn rich contextualized representations by considering the entire input sequence during training. Our study employs the pre-trained DNABert model, which has demonstrated superior performance in several DNA sequence classification tasks. We specifically fine-tune the 6mer-BERT variant on the 4mC methylation site benchmark dataset. Fine-tuning a pre-trained model on a task-specific dataset allows us to transfer the knowledge acquired during pre-training, enabling the model to achieve state-of-the-art performance in predicting DNA 4mC methylation sites. The incorporation of BERT’s pre-trained knowledge provides significant advantages, as the model has already learned from vast amounts of data and captures intricate sequence patterns and dependencies. By leveraging pre-trained models like BERT, we achieve robust and accurate predictions, even in scenarios with limited training data.

Evaluation metrics

In order to compare with previous related work, we selected the commonly used evaluation indicators comprised of accuracy (ACC), sensitivity (SN), specificity (SP), Matthews’ coefficient correlation (MCC), and area under the receiver operating characteristic curve (AUC). These indicators are calculated by the following formula:where TP represents true positives, which is the number of correctly predicted positive samples; TN represents true negatives, which is the number of correctly predicted negative samples; FP represents false positives, which is the number of negative samples wrongly predicted as positive; and FN represents false negatives, which is the number of positive samples wrongly predicted as negative.

Experimental setup

In our experimental design, we adopted the default settings of DeepBIO for other hyperparameters. For instance, when performing data set deduplication, we limited the duplication rate to 0.8 using the CDHIT algorithm integrated in the DeepBIO platform. Furthermore, we conducted a grid search on hyperparameters such as learning rate and batch size for each model. Grid search is a method for hyperparameter tuning, where different combinations of hyperparameters are tried to determine the optimal model configuration. Such experimental settings ensure that the models achieve their maximum potential performance while maintaining the reliability, fairness, and accuracy of the experiments.

Result

In this section, we evaluate the performance of the different models and analyze the features extracted by the different models. In addition, we also compare the features learnt from deep-learning models with the traditional manual feature extraction methods applied in other studies to further demonstrate the superiority of deep learning in solving the 4mC methylation site detection problem. To ensure a balanced representation, the samples were randomly divided into training and test datasets for each species. The division was done in a ratio of 9:1, with 90% of the samples allocated to the training dataset and the remaining 10% assigned to the test dataset.

Performance evaluation of multiple models

We conducted a comprehensive performance evaluation of four different models on six datasets to assess their performance in various data environments. The evaluation process involved the use of common binary classification metrics, such as accuracy (ACC), sensitivity, specificity, area under the curve (AUC), and Matthews correlation coefficient (MCC), to provide a comprehensive understanding of the models’ classification capabilities and highlight their performance differences. In addition to these metrics, we also employed receiver operating characteristic (ROC) curves and precision-recall curves (PRC) to further analyze the models’ performance.

Throughout our evaluation, we observed variations in performance across different datasets. While certain models demonstrated superior predictive performance on most datasets, their performance might vary on specific datasets. As shown in Figures 2A, B, the RNN and TextCNN models exhibited promising performance on the G. pickeringii dataset, while DNABERT outperformed others on the G. subterraneus dataset. Overall, DNABERT consistently showcased superior performance across the evaluated datasets.

FIGURE 2

Furthermore, let’s consider the results obtained on the E. coli dataset. The density distribution of prediction confidence by different deep learning models (Figure 2C) provides insights into the prediction preferences of each model. In the case of LSTM and Text-CNN, their density distribution shows a preference towards the center part of the X-axis, around 0.5. This indicates their poor binary classification ability and confusion in distinguishing between positive and negative instances. On the other hand, the density distribution for DNABERT is skewed towards the right side of the X-axis, indicating a better classification performance. This suggests that DNABERT exhibits a stronger ability to differentiate between positive and negative instances. And this is consistent with the conclusions drawn from the performance comparison in Figure 2A.

We also performed statistics on the overlap of predictions between different models for the same dataset. Take the results obtained on the G. subterraneus dataset as an example, the distribution of sets classified as negative classes by different models in the test set is shown in Figure 2D. In the VN diagram on the left, 41.4% of the test set is judged as negative by all models (negative classes account for 50% of the test set in total). The difference in quantity is shown more clearly in the right figure, and we can find that DNABERT may be one of the less effective models for classification under this dataset, as it predicts more negative cases individually. However, given that most of the model predictions converge on the same, we can conclude that most of the models are consistent in their classification results.

Deep learning model feature analysis

We conduct a comparative study on the features learned by deep learning from biological information. This includes comparisons between different deep learning models as well as comparisons between deep learning features and manually designed features. By conducting feature comparisons, we aim to further validate the superiority of deep learning methods and enhance the interpretability of deep learning models. We select ANF, binary, CKSNAP, and DNC approaches to extracting features and using SVM for unsupervised classification to compare with our deep learning models. Figure 3A presents the ROC and PR curves for all models on the G. pickeringii dataset. We only display the two best-performing traditional manual feature methods for comparison. It is evident that most of the deep learning methods outperform the traditional approaches in terms of classification performance.

FIGURE 3

To visualize the results of deep learning features, we utilized UMAP (Uniform Manifold Approximation and Projection) and SHAP (Shapley Additive Explanations) plots for display (Figure 3B). The UMAP plot reduces the dimensionality of the features while preserving the underlying data structure. It enables data clustering and categorization by mapping high-dimensional features into a lower-dimensional space, allowing for an analysis of feature similarity between positive and negative instances. The SHAP plot facilitates the understanding of feature importance and contribution to model predictions, providing interpretability to the model and enabling comparison of feature impacts. It helps to comprehend the significance of features in model predictions, enhancing interpretability and facilitating comparison among different features. In the feature visualization figure, each row corresponds to a specific feature, and the x-axis represents the snap value, providing a clearer understanding of the feature. The color gradient indicates the feature value, with higher values represented by redder colors and lower values represented by bluer colors. Each line represents a feature, and the horizontal position represents the SHAP value assigned to that feature in a particular sample. Each point represents a sample. The intensity of the color reflects the impact of the feature, with redder colors indicating a larger impact and bluer colors indicating a smaller impact. The scattered distribution of points indicates a greater influence of the feature.

Conclusion

In this study, we use several currently popular deep learning models on the problem of 4mC methylation detection of DNA. We first present the current status of DNA 4mC methylation site detection, followed by the design of deep learning model workflows on six benchmark datasets, and finally, we evaluate the output of all models and conclude that deep learning has great potential for methylation detection, leading the way to future sequencing technologies along with newer bio-experimental methods. In fact, deep learning methods consistently outperformed traditional machine learning methods on all datasets. Furthermore, it was observed that pre-trained deep learning models with a higher number of parameters exhibited even better performance. We believe this may be because deep learning models with more parameters capture more features and analyze the features acquired by each model. By attempting to explain the model’s internal workings and shed light on its internal representations, we aim to define its “black box” behavior.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Author contributions

HJ: Data curation, Validation, Writing–original draft, Writing–review and editing. JB: Writing–review and editing. JJ: Data curation, Writing–review and editing. YC: Data curation, Writing–review and editing. XC: Writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research work was supported by the Innovation Fund of the Ministry of Education’s Engineering Research Center for the Integration and Application of Digital Learning Technologies, under project grant number 1221001.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AoC.JiaoS.WangY.YuL.ZouQ. (2022a). Biological sequence classification: A review on data and general methods. Research2022, 0011. 10.34133/research.0011

  • 2

    AoC.ZouQ.YuL. (2022b). NmRF: Identification of multispecies RNA 2’-O-methylation modification sites from RNA sequences. Briefings Bioinforma.23, bbab480. 10.1093/bib/bbab480

  • 3

    BuryanovY. I.ShevchukT. (2005). DNA methyltransferases and structural-functional specificity of eukaryotic DNA modification. Biochem. Mosc.70, 730742. 10.1007/s10541-005-0178-0

  • 4

    CaiL.WangL.FuX.ZengX. (2021). Active semisupervised model for improving the identification of anticancer peptides. ACS Omega6, 2399824008. 10.1021/acsomega.1c03132

  • 5

    CaoC.WangJ.KwokD.CuiF.ZhangZ.ZhaoD.et al (2022). webTWAS: a resource for disease candidate susceptibility genes identified by transcriptome-wide association study. Nucleic Acids Res.50, D1123D1130. 10.1093/nar/gkab957

  • 6

    ChenJ.ZouQ.LiJ. (2022). DeepM6ASeq-EL: Prediction of human N6-methyladenosine (m6A) sites with LSTM and ensemble learning. Front. Comput. Sci.16, 162302. 10.1007/s11704-020-0180-0

  • 7

    ChenK.ZhaoB. S.HeC. (2016). Nucleic acid modifications in regulation of gene expression. Cell Chem. Biol.23, 7485. 10.1016/j.chembiol.2015.11.007

  • 8

    ChenW.YangH.FengP.DingH.LinH. (2017). iDNA4mC: identifying DNA N4-methylcytosine sites based on nucleotide chemical properties. Bioinformatics33, 35183523. 10.1093/bioinformatics/btx479

  • 9

    ChenY.YangX.WangJ.SongB.ZengX. (2021). Muffin: Multi-scale feature fusion for drug–drug interaction prediction. Bioinformatics37, 26512658. 10.1093/bioinformatics/btab169

  • 10

    DongJ.ZhaoM.LiuY.SuY.ZengX. (2022). Deep learning in retrosynthesis planning: Datasets, models and tools. Briefings Bioinforma.23, bbab391. 10.1093/bib/bbab391

  • 11

    DwyerD. B.FalkaiP.KoutsoulerisN. (2018). Machine learning approaches for clinical psychology and psychiatry. Annu. Rev. Clin. Psychol.14, 91118. 10.1146/annurev-clinpsy-032816-045037

  • 12

    GravesA. (2012). “Long Short-Term Memory,” in Supervised Sequence Labelling with Recurrent Neural Networks. Studies in Computational Intelligence (Berlin, Heidelberg: Springer) 385, 37–45.

  • 13

    GravesA.SchmidhuberJ. (2005). Framewise phoneme classification with bidirectional LSTM and other neural network architectures. Neural Netw.18, 602610. 10.1016/j.neunet.2005.06.042

  • 14

    HamdyR.MaghrabyF. A.OmarY. M. K. (2022). ConvChrome: Predicting gene expression based on histone modifications using deep learning techniques. Curr. Bioinforma.17, 273283. 10.2174/1574893616666211214110625

  • 15

    HasanM. M.ManavalanB.ShoombuatongW.KhatunM. S.KurataH. (2020). i4mC-Mouse: Improved identification of DNA N4-methylcytosine sites in the mouse genome using multiple encoding schemes. Comput. Struct. Biotechnol. J.18, 906912. 10.1016/j.csbj.2020.04.001

  • 16

    HeW.JiaC.ZouQ. (2019). 4mCPred: Machine learning methods for DNA N4-methylcytosine sites prediction. Bioinformatics35, 593601. 10.1093/bioinformatics/bty668

  • 17

    HuY.SunJ. Y.ZhangY.ZhangH.GaoS.WangT.et al (2021). rs1990622 variant associates with Alzheimer's disease and regulates TMEM106B expression in human brain tissues. BMC Med.19, 11. 10.1186/s12916-020-01883-5

  • 18

    HuY.ZhangH.LiuB.GaoS.WangT.HanZ.et al (2020). rs34331204 regulates TSPAN13 expression and contributes to Alzheimer's disease with sex differences. Brain143, e95. 10.1093/brain/awaa302

  • 19

    HuY.ZhangY.ZhangH.GaoS.WangL.WangT.et al (2022a). Cognitive performance protects against Alzheimer's disease independently of educational attainment and intelligence. Mol. Psychiatry27, 42974306. 10.1038/s41380-022-01695-4

  • 20

    HuY.ZhangY.ZhangH.GaoS.WangL.WangT.et al (2022b). Mendelian randomization highlights causal association between genetically increased C-reactive protein levels and reduced Alzheimer's disease risk. Alzheimers Dement.18, 20032006. 10.1002/alz.12687

  • 21

    HuangQ. F.ZhangJ.WeiL. Y.GuoF.ZouQ. (2020). 6mA-RicePred: A method for identifying DNA N (6)-methyladenine sites in the rice genome based on feature fusion. Front. Plant Sci.11, 4. 10.3389/fpls.2020.00004

  • 22

    JiY.ZhouZ.LiuH.DavuluriR. V. (2021). Dnabert: Pre-trained bidirectional encoder representations from transformers model for DNA-language in genome. Bioinformatics37, 21122120. 10.1093/bioinformatics/btab083

  • 23

    JiaoY.DuP. (2016). Performance measures in evaluating machine learning based bioinformatics predictors for classifications. Quant. Biol.4, 320330. 10.1007/s40484-016-0081-2

  • 24

    JinJ.YuY.WangR.ZengX.PangC.JiangY.et al (2022). iDNA-ABF: multi-scale deep biological language learning model for the interpretable prediction of DNA methylations. Genome Biol.23, 219223. 10.1186/s13059-022-02780-1

  • 25

    KimY. (2014). Convolutional neural network for sentence classification[J]. Waterloo, ON: University of Waterloo. arXiv preprint arXiv:1408.5882.

  • 26

    KulisM.EstellerM. (2010). DNA methylation and cancer. Adv. Genet.70, 2756. 10.1016/B978-0-12-380866-0.60002-2

  • 27

    LairdP. W. (2010). Principles and challenges of genome-wide DNA methylation analysis. Nat. Rev. Genet.11, 191203. 10.1038/nrg2732

  • 28

    LarranagaP.CalvoB.SantanaR.BielzaC.GaldianoJ.InzaI.et al (2006). Machine learning in bioinformatics. Briefings Bioinforma.7, 86112. 10.1093/bib/bbk007

  • 29

    LiJ.HeS. D.GuoF.ZouQ. (2021). HSM6AP: A high-precision predictor for the Homo sapiens N6-methyladenosine (m^6 A) based on multiple weights and feature stitching. Rna Biol.18, 18821892. 10.1080/15476286.2021.1875180

  • 30

    LiZ.ZhuS.ShaoB.ZengX.WangT.LiuT. Y. (2023). DSN-DDI: An accurate and generalized framework for drug–drug interaction prediction by dual-view representation learning. Briefings Bioinforma.24, bbac597. 10.1093/bib/bbac597

  • 31

    LiuC.SongJ.OgataH.AkutsuT. (2022). MSNet-4mC: Learning effective multi-scale representations for identifying DNA N4-methylcytosine sites. Bioinformatics38, 51605167. 10.1093/bioinformatics/btac671

  • 32

    LvH.DaoF. Y.ZhangD.YangH.LinH. (2021). Advances in mapping the epigenetic modifications of 5‐methylcytosine (5mC), N6‐methyladenine (6mA), and N4‐methylcytosine (4mC). Biotechnol. Bioeng.118, 42044216. 10.1002/bit.27911

  • 33

    ManavalanB.BasithS.ShinT. H.LeeD. Y.WeiL.LeeG. (2019). 4mCpred-EL: An ensemble learning framework for identification of DNA N4-methylcytosine sites in the mouse genome. Cells8, 1332. 10.3390/cells8111332

  • 34

    MooreL. D.LeT.FanG. (2013). DNA methylation and its basic function. Neuropsychopharmacology38, 2338. 10.1038/npp.2012.112

  • 35

    NiP.HuangN.ZhangZ.WangD. P.LiangF.MiaoY.et al (2019). DeepSignal: Detecting DNA methylation state from nanopore sequencing reads using deep-learning. Bioinformatics35, 45864595. 10.1093/bioinformatics/btz276

  • 36

    PanX.LinX.CaoD.ZengX.YuP. S.HeL.et al (2022). Deep learning for drug repurposing: Methods, databases, and applications. Wiley Interdiscip. Rev. Comput. Mol. Sci.12, e1597. 10.1002/wcms.1597

  • 37

    PlongthongkumN.DiepD. H.ZhangK. (2014). Advances in the profiling of DNA modifications: Cytosine methylation and beyond. Nat. Rev. Genet.15, 647661. 10.1038/nrg3772

  • 38

    RazinA.CedarH. (1991). DNA methylation and gene expression. Microbiol. Rev.55, 451458. 10.1128/mr.55.3.451-458.1991

  • 39

    RenS. J.YuL.GaoL. (2022). Multidrug representation learning based on pretraining model and molecular graph for drug interaction and combination prediction. Bioinformatics38, 43874394. 10.1093/bioinformatics/btac538

  • 40

    RumelhartD. E.HintonG. E.WilliamsR. J. (1986). Learning representations by back-propagating errors. nature323, 533536. 10.1038/323533a0

  • 41

    SharmaA. K.SrivastavaR. (2021). Protein secondary structure prediction using character bi-gram embedding and Bi-LSTM. Curr. Bioinforma.16, 333338. 10.2174/1574893615999200601122840

  • 42

    SongB.LuoX.LuoX.LiuY.NiuZ.ZengX. (2022). Learning spatial structures of proteins improves protein–protein interaction prediction. Briefings Bioinforma.23, bbab558. 10.1093/bib/bbab558

  • 43

    TranH. V.NguyenQ. H. (2022). iAnt: Combination of convolutional neural network and random forest models using PSSM and BERT features to identify antioxidant proteins. Curr. Bioinforma.17, 184195. 10.2174/1574893616666210820095144

  • 44

    WangR.JiangY.JinJ.YinC.YuH.WangF.et al (2022). DeepBIO is an automated and interpretable deep-learning platform for biological sequence prediction, functional annotation, and visualization analysis, 2022.2009.2029.509859. bioRxiv. 10.1101/2022.09.29.509859

  • 45

    XiaoZ. C.WangL. Z.DingY. J.YuL. A. (2022). iEnhancer-MRBF: Identifying enhancers and their strength with a multiple Laplacian-regularized radial basis function network. Methods208, 18. 10.1016/j.ymeth.2022.10.001

  • 46

    XuH.JiaP.ZhaoZ. (2021). Deep4mC: Systematic assessment and computational prediction for DNA N4-methylcytosine sites by deep learning. Briefings Bioinforma.22, bbaa099. 10.1093/bib/bbaa099

  • 47

    XuJ.MengY.LuC.CaiL.ZengX.NussinovR.et al (2023). Graph embedding and Gaussian mixture variational autoencoder network for end-to-end analysis of single-cell RNA sequencing data. Cell Rep. Methods3, 100382. 10.1016/j.crmeth.2022.100382

  • 48

    YuY.HeW.JinJ.XiaoG.CuiL.ZengR.et al (2021). iDNA-ABT: advanced deep learning model for detecting DNA methylation with adaptive features and transductive information maximization. Bioinformatics37, 46034610. 10.1093/bioinformatics/btab677

  • 49

    ZengX.WangF.LuoY.KangS. G.TangJ.LightstoneF. C.et al (2022a). Deep generative molecular design reshapes drug discovery. Cell Rep. Med.4, 100794. 10.1016/j.xcrm.2022.100794

  • 50

    ZengX.XiangH.YuL.WangJ.LiK.NussinovR.et al (2022b). Accurate prediction of molecular properties and drug targets using a self-supervised image representation learning framework. Nat. Mach. Intell.4, 10041016. 10.1038/s42256-022-00557-6

  • 51

    ZengX.ZhuS.LuW.LiuZ.HuangJ.ZhouY.et al (2020). Target identification among known drugs by deep learning from heterogeneous networks. Chem. Sci.11, 17751797. 10.1039/c9sc04336e

  • 52

    ZhangY.LinJ.ZhaoL.ZengX.LiuX. (2021). A novel antibacterial peptide recognition algorithm based on BERT. Briefings Bioinforma.22, bbab200. 10.1093/bib/bbab200

  • 53

    ZhaoZ.ZhangX.ChenF.FangL.LiJ. (2020). Accurate prediction of DNA N4-methylcytosine sites via boost-learning various types of sequence features. BMC genomics21, 627. 10.1186/s12864-020-07033-8

  • 54

    ZulfiqarH.HuangQ. L.LvH.SunZ. J.DaoF. Y.LinH. (2022b). Deep-4mCGP: A deep learning approach to predict 4mC sites in geobacter pickeringii by using correlation-based feature selection technique. Int. J. Mol. Sci.23, 1251. 10.3390/ijms23031251

  • 55

    ZulfiqarH.SunZ. J.HuangQ. L.YuanS. S.LvH.DaoF. Y.et al (2022a). Deep-4mCW2V: A sequence-based predictor to identify N4-methylcytosine sites in Escherichia coli. Methods203, 558563. 10.1016/j.ymeth.2021.07.011

Summary

Keywords

4mC DNA methylation, deep learning, classification, feature, visualization, interpretable ability

Citation

Ju H, Bai J, Jiang J, Che Y and Chen X (2023) Comparative evaluation and analysis of DNA N4-methylcytosine methylation sites using deep learning. Front. Genet. 14:1254827. doi: 10.3389/fgene.2023.1254827

Received

07 July 2023

Accepted

31 July 2023

Published

21 August 2023

Volume

14 - 2023

Edited by

Lei Chen, Shanghai Maritime University, China

Reviewed by

Lesong Wei, King Abdullah University of Science and Technology, Saudi Arabia

Hao Lin, University of Electronic Science and Technology of China, China

Updates

Copyright

*Correspondence: Xin Chen, ; Jie Bai,

† These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics