REVIEW article

Front. Med., 08 April 2025

Sec. Precision Medicine

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1503229

DNA sequence analysis landscape: a comprehensive review of DNA sequence analysis task types, databases, datasets, word embedding methods, and language models

  • 1. German Research Center for Artificial Intelligence GmbH, Kaiserslautern, Germany

  • 2. Intelligentx GmbH (intelligentx.com), Kaiserslautern, Germany

  • 3. Department of Computer Science, Technical University of Kaiserslautern, Kaiserslautern, Germany

Abstract

Deoxyribonucleic acid (DNA) serves as fundamental genetic blueprint that governs development, functioning, growth, and reproduction of all living organisms. DNA can be altered through germline and somatic mutations. Germline mutations underlie hereditary conditions, while somatic mutations can be induced by various factors including environmental influences, chemicals, lifestyle choices, and errors in DNA replication and repair mechanisms which can lead to cancer. DNA sequence analysis plays a pivotal role in uncovering the intricate information embedded within an organism's genetic blueprint and understanding the factors that can modify it. This analysis helps in early detection of genetic diseases and the design of targeted therapies. Traditional wet-lab experimental DNA sequence analysis through traditional wet-lab experimental methods is costly, time-consuming, and prone to errors. To accelerate large-scale DNA sequence analysis, researchers are developing AI applications that complement wet-lab experimental methods. These AI approaches can help generate hypotheses, prioritize experiments, and interpret results by identifying patterns in large genomic datasets. Effective integration of AI methods with experimental validation requires scientists to understand both fields. Considering the need of a comprehensive literature that bridges the gap between both fields, contributions of this paper are manifold: It presents diverse range of DNA sequence analysis tasks and AI methodologies. It equips AI researchers with essential biological knowledge of 44 distinct DNA sequence analysis tasks and aligns these tasks with 3 distinct AI-paradigms, namely, classification, regression, and clustering. It streamlines the integration of AI into DNA sequence analysis tasks by consolidating information of 36 diverse biological databases that can be used to develop benchmark datasets for 44 different DNA sequence analysis tasks. To ensure performance comparisons between new and existing AI predictors, it provides insights into 140 benchmark datasets related to 44 distinct DNA sequence analysis tasks. It presents word embeddings and language models applications across 44 distinct DNA sequence analysis tasks. It streamlines the development of new predictors by providing a comprehensive survey of 39 word embeddings and 67 language models based predictive pipeline performance values as well as top performing traditional sequence encoding-based predictors and their performances across 44 DNA sequence analysis tasks.

1 Introduction

Deoxyribonucleic acid (DNA) functions as the blueprint of life as it contains essential instructions for the development, operation, growth, and reproduction of all living organisms (). Organisms utilize cell division process to grow from fertilized egg to a multicellular adult. Throughout an organism's lifespan, the health of tissues and organs is maintained through a continuous cycle of cell replacement. In this cycle, worn-out or damaged cells are systematically replaced with new, healthy cells. When a cell divides, each new cell requires an exact copy of the DNA to function correctly (). DNA replication and repair processes ensure that each daughter cell receives the same genetic information as the parent cell, which is essential for the survival and proper functioning of all living organisms (). DNA sequence changes occur through two fundamental mechanisms: germline mutations inherited from parents and somatic mutations acquired during an individual's lifetime (). Germline mutations are present in all cells and can be passed to offspring, underlying hereditary conditions. Somatic mutations occur post-conception and can be caused by various factors including internal factors such as cellular metabolites, replication errors, and spontaneous chemical changes and external factors such as ionizing radiation, chemical mutagens, environmental pollutants, and lifestyle factors (, ). Understanding these distinct mutation types is crucial as they require different analytical approaches. Germline mutation analysis typically involves comparing an individual's sequence to population databases, while somatic mutation analysis often requires comparing affected tissue to unaffected tissue from the same individual. Regardless of type, mutations in genetic information can lead to complex diseases and disorders such as cancer (). To detect susceptibility, initiation, and progression of such diseases at early stages, scientists perform large-scale DNA sequence analysis (). Through DNA sequence analysis, scientists can decode the intricate genetic data by uncovering the origins of genetic mutations and disorders (). In addition, this analysis is crucial for the development of targeted therapies and the advancement of personalized medicine ().

DNA sequence analysis through traditional wet-lab experiments is expensive and time-consuming (, ). This is because wet-lab experiments require specialized equipment, e.g., PCR machines, and costly reagents (e.g., enzymes and chemicals). Detailed experiments on multiple patient samples may take weeks or even months. Moreover, experimentation requires careful execution and validation to prevent incorrect interpretations of genetic mutations due to errors or inconsistencies. The influx of next-generation sequencing and high-throughput approaches has given rise to huge sequences data. This abundance of genomic information has created both opportunities and challenges for comprehensive analysis. To expedite genomics sequence analysis, researchers are analyzing publicly available sequences data by harnessing the capabilities of Artificial Intelligence (AI) methods. It is important to mention that AI approaches serve to augment rather than replace experimental methods in DNA sequence analysis. For example, in precision medicine, AI models trained on large genomic databases can help to interpret patient-specific data by identifying relevant patterns and potential functional impacts. However, patient-specific experimental data remain essential, particularly for understanding unique aspects of individual cases such as tumor mutations. Thus, AI methods provide a valuable tool for generating hypotheses and guiding experimental design while working in concert with traditional molecular biology approaches.

While DNA sequence analysis encompasses a broad range of computational approaches in bioinformatics, from genome assembly and variant detection to evolutionary analysis and microbiome studies, this review focuses specifically on DNA sequence analysis tasks that involve pattern recognition and prediction, where artificial intelligence approaches can be effectively applied. These tasks include predicting functional elements, identifying regulatory regions, and classifying sequence types applications where AI can learn complex sequence patterns that may not be apparent through traditional computational methods.

Most of the AI-based genomics sequence analysis methods fall under the hood of regression and classification paradigms (). Figure 1 illustrates a unified workflow of AI-based predictive pipelines for genomics sequence analysis tasks. It is evident in the Figure that, overall, AI predictive pipelines can be divided into 4 different stages (). First stage emphasizes on the collection and development of quality benchmark datasets using public databases (). Second stage focuses on the characterization of raw DNA sequences in terms of statistical vectors using different kinds of sequence encoders (). This is primarily done to address the inherent dependency of AI predictive pipelines on statistical vectors (). In entire predictive pipeline, this stage is the most crucial one because highly informative and discriminative statistical vectors help the predictors to learn comprehensive useful patterns for accurate prediction (). It is widely accepted that with quality statistical vectors, even simple machine learning predictors can produce promising performance. On contrary, with less informative and discriminative statistical representations, even sophisticated deep learning predictors fail to produce decent performance ().

Figure 1

There is a marathon of developing powerful sequence encoders for generating highly informative and discriminative statistical vectors of raw sequences. To date, hundreds of sequence encoding methods have been developed () that can be broadly classified into four categories: Physico-chemical properties based methods, statistical methods (, ), neural word embedding methods (), and language models (22). While physico-chemical properties based methods generate statistical vectors of raw sequences using pre-computed physical and chemical values of nucleotides, statistical methods rely on occurrence frequencies of individual or group of nucleotides with DNA sequences (). Physico-chemical properties based and statistical methods capture the intrinsic characteristics of biological sequences, such as nucleotide composition and distributional information. However, these methods lack to capture complex relationships of nucleotides such as long range interactions of nucleotides in the sequences (, 23). In addition, these methods may not fully capture the semantic and functional similarities between sequences (, 23). Neural word embedding methods learn distributed representations of nucleotides in the continuous vector space. These methods capture the syntactic and semantic similarities of nucleotides by mapping them to vectors in a high-dimensional space. This enables the representation of residues with similar contexts to be closer together in the vector space. Neural word embeddings methods efficiently capture semantic and contextual information of nucleotides. However, these methods lack to efficiently handle different contexts of same nucleotides (). Language models also learn representation of individual nucleotides or groups of nucleotides (k-mers) in an unsupervised fashion by predicting masked nucleotides based on the context of surrounding nucleotides. Language models based methods capture complex nucleotide relations; however, these methods require large amount of sequence data for training and hyperparameter optimization (22).

Third stage includes predictors that make best use of statistical vectors produced by second stage to extract informative patterns for creating decision boundaries. Overall, these predictors can be classified into two categories: machine learning and deep learning (). Machine learning predictors require less data and computational power for training. However, these predictors lack to capture comprehensive complex relationships of nucleotide (), whereas deep learning predictors (24) are capable to learn highly complex relationships of nucleotide. However, these predictors require a huge amount of training data and computational power (). In fourth stage, comprehensive evaluation of predictors using different experimental settings and evaluation measures is performed (24).

AI researchers have been endeavoring to complement wet-lab-based DNA sequence analysis methods by incorporating more innovative sequence encoders at second stage and predictors at third stage of predictive pipeline. However, there is still ample room for the development of more powerful predictive pipelines. Different fields such as Natural Language Processing (NLP), Energy, and Computer Vision have seen substantial progress in the development of diverse predictive pipelines. Whereas, the DNA sequence analysis field is known for its wide range of tasks, still the progress of AI applications in this area is hindered mainly due to the lack of integration between molecular biologist and AI experts. For instance, the field of NLP has made strides with multi-task learning predictors. However, the DNA sequence analysis field lags behind due to AI experts limited understanding of the diverse range of DNA analysis tasks that could support the development of multi-task learning predictors. Furthermore, the efficacy of AI applications hinges on the availability of benchmark datasets. Although developing datasets in DNA sequence analysis is relatively straightforward due to abundance of public databases which contain raw biological sequences along with associated labels, there is a tendency among researchers to overlook existing benchmark datasets, develop new benchmark datasets, and neglect comprehensive performance comparisons with existing predictors. This oversight often complicates the determination of the most effective predictors for specific tasks. For example, up to date, according to our best of knowledge, approximately 127 predictive models have been developed and published in 59 different conferences and journals for widely studied 44 different DNA sequence analysis tasks. To enhance the performance of predictive models developed for diverse DNA sequence analysis tasks, researchers need to conduct a comprehensive examination of existing literature to find most effective algorithms for different stages of new predictive pipelines. With an aim to expedite progress in the development of fair and robust AI applications for DNA sequence analysis, numerous review articles have emerged. However, these reviews typically focus on isolated tasks rather than providing a holistic overview. Considering the need and significance of a comprehensive study that bridges the gap between AI specialists and biologists, this paper makes manifold contributions:

  • It bridges the gap between DNA sequence analysis and artificial intelligence fields by presenting a diverse range of DNA analysis tasks and AI methodologies.

  • It empowers AI researchers by equipping them with essential biological knowledge related to 44 distinct DNA sequence analysis tasks. It categorizes 44 different DNA sequence analysis tasks into 8 different categories on the basis of sequence analysis goals. This categorization provides a structured overview to biologists and AI researchers in navigating the complex landscape of genomics studies more efficiently.

  • It streamlines the integration of AI into DNA sequence analysis by consolidating information of 36 diverse biological databases being used to develop benchmark datasets for 44 different DNA sequence analysis tasks.

  • It sheds light on the nature of 44 different DNA sequence analysis tasks and categorizes them into three primary categories: regression, classification, and clustering, and three secondary categories: binary classification, multi-class classification, and multi-label classification. This categorization assists computer scientists in efficient selection of most suitable algorithms for each task category, development of more effective and specialized computational frameworks, and to significantly accelerate advancements in AI-driven genomic research.

  • It provides insights of 140 benchmark datasets related to 44 distinct DNA sequence analysis tasks to ensure performance comparisons between new and existing AI predictors.

  • It presents word embeddings and language models applications for 44 distinct DNA sequence analysis tasks.

  • It streamlines the development of new predictors by providing a comprehensive survey of current top predictors, their performances across 44 DNA sequence analysis tasks, and their public accessibility. This comprehensive overview serves as a valuable resource for researchers developing and validating predictive pipelines in computational genomics.

It is important to note that our categorization of 44 DNA sequence analysis tasks emerges from the AI and computational biology literature rather than representing a definitive biological taxonomy. We have organized these tasks into biologically relevant groupings based on their functional and analytical similarities, while recognizing that many tasks span multiple biological domains. This organization aims to bridge the gap between computational methodologies and biological applications, although we acknowledge that future refinements with deeper domain expert input would further enhance this framework.

2 Research methodology

This section provides a detailed overview of the research methodology used to identify articles focused on word embeddings and large language models applications in DNA sequence analysis landscape (, ). Figure 2 illustrates two stage processes for article identification and selection.

Figure 2

2.1 Article searching

To identify a wide range of relevant scholarly articles, initial stage involves formulation of quality search queries using different keywords. In Figure 2, article identification module contains keywords cell of three different categories, namely, DNA tasks, word embedding methods, and Language models. To formulate quality search queries, keywords within same category are combined using OR ∨ operator, while keywords of different categories are combined using AND ∧ operator. For instance, few sample search queries include DNA Replication Origins Identification using BERT language model, DNA Replication Origins Identification using DeepWalk word embedding method, etc. To acquire relevant papers, formulated search queries are executed on academic search engines such as Google Scholar,1 ACM Digital Library,2 Elsevier,3 IEEEXplore,4 Wiley Online Library,5 Springer,6 and ScienceDirect.7 In addition, snowballing method is employed to explore sources referenced in extracted papers to identify more research articles. This technique is particularly useful in research contexts where access to resources is limited, such as niche topics or hard-to-reach communities, as it expands the pool of resources for a study. Execution of queries across multiple academic databases acquired approximately 238 research articles which are screened and filtered in second stage.

2.2 Article screening and filtering

Second stage selects most relevant articles in two steps. In the first step, titles and abstracts of 113 word embeddings and 125 large language models related studies were reviewed. This review analysis identified 80 word embeddings and 104 language models related relevant articles. Second step involves full-text assessment of articles selected in first step, resulting in 39 word embeddings and 67 language models related articles.

Our selection criteria focused on DNA sequence analysis tasks where (1) raw DNA sequence data serve as the primary input, (2) AI methods extract patterns from these sequences, and (3) the analysis predicts specific biological properties or functions. This allowed us to examine AI's impact on genomic sequence interpretation while acknowledging that bioinformatics encompasses many other types of analyses not covered here.

3 Biological foundations of DNA sequence analysis goals and tasks

With an aim to find molecular basis of diseases initiation and progression, their effective detection at early stages, and development of potent drugs, researchers are trying to understand DNA sequence language by performing a variety of sequence analysis tasks. Every unique DNA sequence analysis task aims to enhance the understanding of one specific aspect of DNA, and a bunch of tasks can enhance the understanding of specific major biological goal. To summarize the biological background of 44 distinct DNA sequence analysis tasks, we have categorized them into 8 major biological goals. Figure 3 depicts the biological categorization of 44 unique DNA sequences analysis tasks into 8 different goals, namely, genome structure and stability, gene expression regulation, gene analysis, gene network analysis, DNA modification prediction, DNA functional analysis, environmental and microbial genomics, and disease analysis. This biologically informed organization was developed by analyzing both the computational biology literature and aligning with biological processes in genomics research. While computational researchers often approach these tasks through the lens of AI methodologies, we have endeavored to categorize them according to their biological relevance and function. Our categorization into 8 major biological goals represents an attempt to bridge computational approaches with biological understanding. Although we recognize the inherent complexity and interconnectedness of biological systems which indicates that many tasks could reasonably be classified in multiple categories, thus, this categorization represents one of several possible ways to organize these tasks. This categorization reflects the diverse biological applications where AI-based sequence analysis has made significant contributions. However, we recognize that DNA sequence analysis in bioinformatics extends beyond these pattern recognition tasks to include other critical applications such as genome assembly, variant detection, and population genetics studies. We specifically examine how modern AI approaches are transforming our ability to extract meaningful biological insights from sequence data through pattern-based prediction tasks.

Figure 3

In living organisms, DNA is packaged at multiple levels to condense vast genetic information into a well-organized structure within the cell nucleus (). At the first level, DNA is wrapped around histone octamers also known as nucleosomes. These nucleosomes further assemble into chromatin, which then folds and condenses into an even more compact structure known as the genome (). The exploration of genome structure and stability is pivotal in understanding the biological intricacies and potential therapeutic avenues. Genome structure can affect how genes are accessed and used. Disruptions in this structure, such as missing or misplaced DNA sections, or changes in how tightly DNA is wrapped around histone octamers, or irregularities in nucleosomes positions can lead to genes being turned on or off at the wrong times or in the wrong amounts (). This can cause various diseases and biological disorders. DNA is an instruction manual that controls biological functioning within living organisms. If genome gets unstable, the manual gets messed up such as typos and missing sections. It can lead to uncontrolled growth of the cells (cancer) and improper working of the genes (many diseases) (). In a nutshell, a stable genome possesses clear, complete instruction manual, essential for keeping biological functions working smooth. To better understand genome structure and stability, it is essential to explore various tasks such as DNA Replication Origins Prediction (25, 26), Genome Structure Analysis (27, 28), Nucleosome Position Detection (29, 30), Chromatin Accessibility Prediction (3133), Chromatin Feature Prediction (31, 34, 35), Long-range Chromatin Interaction Prediction (36, 37), and YY1-Mediated Chromatin Loops Prediction (38, 39). These tasks are crucial for comprehending the intricate mechanisms governing genetic information processing and regulation within cells (40).

DNA replication origin prediction is fundamental as accurate replication of the genome is vital for maintaining genomic stability (25). The prediction of replication origins involves calculating DNA structural properties to identify sites crucial for initiating DNA replication (25). Understanding where these sites are located and how they are specified is essential for comprehending DNA replication and ensuring genome integrity (41). Genome structure analysis plays a pivotal role in deciphering the organization and arrangement of genetic material within the cell (27). By analyzing the structural features of the genome, researchers can gain insights into the functional and spatial organization of chromosomes, aiding in the identification of genomic elements involved in gene regulation and phenotypic variations (27, 42). Furthermore, nucleosome position detection is essential for understanding how nucleosomes, the basic units of genome, are arranged along the DNA strand (29, 43). This information is crucial for elucidating gene regulation mechanisms and chromatin dynamics within the cell (29, 43). Chromatin accessibility prediction is a key task that involves determining the regions of chromatin that are accessible for transcription factors and other regulatory proteins to bind (3133). Prediction of chromatin accessibility across different cellular contexts provides valuable insights into gene regulation and chromatin dynamics (3133). Chromatin feature prediction complements accessibility prediction by identifying specific chromatin features and epigenetic markers that influence gene expression and regulatory processes (31, 34, 35, 44). These features include transcription factor (TF) binding sites, DNase I-hypersensitive sites (DHS), and histone marks (HM). By understanding these features, researchers can unravel the mechanisms underlying chromatin regulation and gene expression (34). Long-range chromatin interactions make bridges between distant enhancers and promoters. These interactions enable interactions between enhancers and promoters by bringing them closer to each other (36, 37). YY1-mediated chromatin loop prediction provides comprehensive understanding about gene regulation (38, 39, 45). YY1 is a protein that makes loop between enhancers and promoters. These loops are essential for gene regulation, and by predicting these loops, we can see which genes can be controlled through YY1 protein (38, 39, 45). This knowledge is valuable for understanding diseases where gene regulation goes wrong. To sum up, only through multi-dimensional exploration of genome structure and stability, researchers can discriminate healthy cellular processes from malfunctioned processes, find the root causes of diseases, and develop potent therapies.

Another major goal of molecular biologists behind is gene expression regulation. Gene expression regulation provides fundamental insights into how genes are activated or repressed in response to various cellular cues (46). Specifically, researchers are trying to unravel the intricate mechanisms that control when and up to what extent specific genes are turned on or off in different cells and tissues (46). This knowledge forms the basis for understanding the functional behavior of genes in different biological contexts and sets the stage for further analyses. Hence, it holds immense promise for scientists and pharmaceutical industries. This helps scientists to detect irregularities in normal gene expression regulation, the way diseases develop at the molecular level, and identify potential drug targets (46). Furthermore, this understanding can assist pharmaceutical industries to develop improved diagnostic tools, innovative personalized therapies, and targeted interventions, which will ultimately contribute to advancements in personalized healthcare (46). In addition, it can provide a deeper understanding of biological systems which can lead to breakthroughs in biotechnology (46). For better understanding of gene expression regulation, researchers are performing nine different DNA sequence analysis tasks, namely, enhancer identification (47), promoter identification (48), enhancer-promoter interactions prediction (49), transcription site prediction (50), transcription factor binding site prediction (51), transcription factor binding affinity prediction (52), protein-DNA binding site prediction (53), splice site prediction (53), and translation initiation site prediction (54). Enhancers (47, 5575) and promoters identification (48, 7681), along with their interactions (8286) prediction are important to decipher a complex control panel for gene expression (4749). Enhancers are known as distant switches of genes, while promoters are the landing sites where gene activation starts. Identification of these elements and predicting how they loop together provide a comprehensive understanding of gene regulation, including which genes are activated or repressed, the intensity of their expression, and the specific cell types involved (87, 88). This knowledge reveals the intricate regulatory code that governs gene expression and offers valuable insights into the mechanisms underlying normal cellular function as well as the dysregulation that may contribute to various diseases.

Furthermore, prediction of different genomic sites including transcription sites (50), transcription factor binding sites (8993), transcription factor binding site affinity (52), protein-DNA binding site (53, 9496), splice site (93, 97100), and translation initiation site (50, 101) provide deep insights into gene expression regulation. A transcription site refers to the specific location on the DNA where the process of transcription takes place. Transcription is the synthesis of RNA from a DNA template, and the transcription site represents the region where the RNA polymerase enzyme binds and initiates the transcription process, whereas transcription factor binding sites are specific DNA sequences where transcription factors (proteins), that regulate gene expression, bind. These binding sites are typically located near the transcription start site and are recognized by transcription factors to control the initiation or repression of transcription. In contrast, transcription factor binding site affinity refers to the strength or affinity with which a transcription factor binds to its specific binding site on DNA. It represents the likelihood of a transcription factor binding to its target site and influencing gene expression. A protein-DNA binding site refers to any region on the DNA where a protein binds. This can include transcription factors, as mentioned earlier, as well as other proteins involved in various cellular processes such as DNA replication, repair, and chromatin remodeling. Splice sites are specific sequences within a gene's DNA that mark the boundaries of introns and exons. During the process of RNA splicing, introns are removed from the pre-mRNA molecule, and exons are joined together to form the mature mRNA. Splice sites are essential for the accurate and precise splicing of RNA. Translation initiation site (TIS) is the specific location on the mRNA molecule where the process of translation begins. TIS prediction seems like a RNA sequence analysis task; however, in molecular biology research, to study gene expression, researchers are synthesizing complementary DNA (cDNA) data from messenger RNA (mRNA) template through a process called reverse transcription. In the context of cDNA data, the translation initiation site (TIS) represents the position where the ribosome, the cellular machinery responsible for protein synthesis, binds to the mRNA to initiate translation. The TIS is typically identified by the presence of specific start codons, such as AUG, which serve as signals for the ribosome to start protein synthesis.

To better understand gene functions and their roles in disease initiation, researchers are exploring various aspects such as gene expression prediction (102, 103), identification of essential (104109) and disease-specific genes (110), gene function prediction (111, 112), pseudo-gene function prediction (111), target gene classification (113), and candidate gene prioritization (114). Overall together, these tasks provide a comprehensive platform for disease diagnosis and development of treatment strategies by uncovering disease mechanisms, identifying potential therapeutic targets, and organizing genes into functional categories. Specifically, gene expression prediction provides useful information about the level of gene activity in different cells or tissues (115). This task is vital for understanding the molecular mechanisms underlying complex diseases such as cancer and identifying potential therapeutic targets. Essential gene identification is another critical task in gene analysis that helps researchers pinpoint genes that are crucial for an organism's survival and development (116, 117). This task is particularly important in understanding gene function and the genetic basis of various disorders. Gene function prediction elucidates the roles of genes in different pathways and biological processes and provides valuable insights into disease mechanisms and potential therapeutic interventions.

Apart from gene function prediction, pseudo-gene function prediction has gained a lot of attention as a critical task in gene analysis (111). Pseudogenes were once thought to be useless DNA because they cannot code for proteins due to mutations that happened over time. However, recent studies have shown that pseudogenes actually play important roles in controlling genes, especially in cancer. For instance, the pseudogene PTENP1 helps to regulate the tumor suppressor gene PTEN in various cancer conditions, showing that pseudogenes can have important functions. Pseudogene function prediction offers numerous advantages, including better understanding of gene regulation, disease mechanisms, evolutionary biology, and the potential for new biomarkers and drug targets. In addition, disease gene prediction is a pivotal task in gene analysis that focuses on identifying genes associated with specific diseases or disorders (118). By pinpointing disease-related genes, researchers can unravel the genetic basis of diseases, discover novel biomarkers for diagnosis and prognosis, and develop targeted therapies. This task is instrumental in precision medicine approaches, where understanding the genetic underpinnings of diseases is crucial for personalized treatment strategies. Target gene classification involves categorizing genes based on their functions, interactions, or regulatory mechanisms (119). By classifying target genes, researchers can better understand gene networks, signaling pathways, and biological processes. This task is essential for deciphering the complex relationships between genes and their roles in health and disease. Candidate gene prioritization and selection are critical tasks in gene analysis that aim to identify genes with the highest likelihood of being involved in a particular biological process or disease (120). By prioritizing candidate genes, researchers can focus their efforts on studying genes that are most likely to have significant effects, accelerating the discovery of novel gene functions and disease mechanisms. This task is crucial for efficiently allocating research resources and maximizing the impact of genetic studies. Aforementioned seven DNA sequence analysis tasks are essential for advancing our understanding of genes and their roles in health and disease. By leveraging these tasks, researchers can unravel the complexities of the genome, uncover novel gene functions, and pave the way for innovative diagnostic and therapeutic strategies in various fields of biology and medicine.

Furthermore, gene network analysis is a promising goal that seeks to comprehend the intricate interactions and relationships between genes within a biological system. Two primary tasks within Gene Network Analysis are Gene Taxonomy Classification and Gene Network Reconstruction. Gene Taxonomy Classification (121123) involves categorizing genes based on their evolutionary relationships and functional similarities, providing a structured framework for organizing genetic information. Gene Taxonomy Classification plays a crucial role in gene network analysis by offering a foundational structure for understanding the evolutionary history and functional relationships between genes. By classifying genes into taxonomic groups based on shared characteristics and evolutionary relatedness, researchers can infer valuable insights into the origins and evolutionary trajectories of genes within a network (124). This classification allows for the identification of core genes that have remained conserved throughout evolution, providing a basis for inferring phylogenetic relationships and understanding the fundamental building blocks of gene networks. Moreover, Gene Taxonomy Classification enables researchers to utilize existing knowledge about gene functions and evolutionary relationships to guide Gene Network Reconstruction. By categorizing genes into taxonomic groups, researchers can pinpoint gene clusters with similar functions or evolutionary origins, facilitating the identification of modules within gene networks that exhibit coordinated activity (125). This classification serves as a roadmap for exploring the functional roles of genes within a network and understanding how these roles have evolved over time. On the other hand, Gene Network Reconstruction (126128) involves creating a detailed map of the interactions and regulatory relationships between genes within a cell or an organism. The primary input for gene network reconstruction is gene expression data obtained through high-throughput techniques such as RNA sequencing (RNA-seq) or microarrays. This task is pivotal for understanding how genes work together to control various biological functions and processes (129). By reconstructing gene networks, researchers can uncover key regulatory hubs involving highly connected genes, clusters of closely interacting genes, pathways, and interactions that steer cellular functions and responses to external stimuli (130).

DNA modification prediction is also a crucial goal where researchers aim is to decipher how tiny tweaks to the DNA code can lead to big changes in cellular functions (131133). In DNA modifications, distinct chemical groups are added to specific locations on the DNA molecule. These additions do not change the actual sequence of nucleotides (A, C, G, T) but can alter the physical properties of DNA sequence. Understanding these modifications, such as 4-Methylcytosine (4mc) (134143), Methyladenine (6ma) (144151), 5-methylcytosine (5mc) (152, 153), 5-hydroxymethylcytosine (5hmc) (154157), and methylation modifications (146, 154159), are essential for advancing our comprehension of epigenetic regulation (160162). Specifically, methylation modifications that occur due to the addition of methyl groups to DNA molecules play a pivotal role in regulating gene expression and maintaining genomic integrity. Similarly, methyladenine modifications, such as DNA N6-methyladenine (6mA), occur due to the addition of a methyl group to the adenine base of DNA. DNA 6mA modifications dynamically influence DNA thermal stability, curvature, and transcription factor interactions, impacting gene expression in a heritable manner. Understanding the prediction of 6mA sites is pivotal for both basic and clinical research as it aids in the identification of gene expression patterns and potential epigenetic changes induced by environmental factors. These predictions enhance our ability to study the role of 6mA modifications in diseases and could lead to improved therapeutic strategies, highlighting the relevance of accurate prediction methods in unraveling the complexities of DNA modifications. Moreover, 5-methylcytosine (5mc) modification occurs due to the addition of a methyl group to the cytosine base of DNA, whereas 5-hydroxymethylcytosine (5hmc) modification is an oxidized derivative of 5mc, where an additional hydroxyl group (-OH) is added to the methyl group of 5mc. Prediction of 5-methylcytosine (5mc) and 5-hydroxymethylcytosine (5hmc) modifications is essential for decoding their roles in gene regulation, developmental processes, and disease states. These critical epigenetic modifications are dynamically regulated by enzymes and influence gene expression crucial for neuronal differentiation and cellular proliferation. Abnormal levels of these modifications have been linked to diseases such as cancer. Precise prediction of 5mc and 5hmc sites is useful for the development of targeted therapies and improved prognostic assessments.

Functional genomics is also a critical goal that encompasses multiple sub-tasks including species classification (44), conserved non-coding element (NCE) classification (163), functional prioritization of non-coding variants (34), prediction of context specific functional impact of genetic variants (36), exon and intron region classification (164), and recombination spots identification (165). Each of these tasks plays a vital role in unraveling the complexities of genetic regulation and molecular mechanisms within the genome. In biomedical research, understanding the genetic similarities and differences between humans and other species is crucial for modeling diseases and studying genetic disorders. Majority of the genome is conserved across different species which makes it difficult to distinguish humans and non-human species. Despite very high genetic similarity across species (< 10% sequence divergence), small differences are extremely valuable and they have significant biological implications. Species classification determines the source species of genetic sequences based on such differences and pave way for better modeling diseases and studying genetic disorders (44). Conserved non-coding element classification is another critical task in functional genomics that focuses on identifying and understanding non-coding regions of the genome that are evolutionarily conserved across different species (163). It is essential for advancing our understanding of gene regulation, evolutionary biology, and the genetic basis of diseases. By elucidating the functions of these non-coding regions, researchers can gain insights into the intricate regulatory networks that govern gene expression and cellular processes and contribute to the development of targeted therapies.

Functional prioritization of non-coding variants (34) is another crucial task for making sense of the vast amount of genetic data generated by modern sequencing technologies. By identifying which variants have significant biological impacts, researchers can gain a deeper understanding of the genetic architecture of complex diseases, uncover novel therapeutic targets, and advance the field of precision medicine. This prioritization is essential for translating genomic research into practical health benefits and ultimately improving patient outcomes and advancing our knowledge of human biology (34). As functional prioritization of non-coding variants task involves identifying which non-coding variants among millions are likely to have functional consequences, it does not account for the specific context in which these variants might exert their effects, whereas prediction of context-specific functional impact of genetic variants aims to provide a detailed understanding of how specific variants influence gene function in different contexts (e.g., specific tissue) (36). This is particularly important for genetic studies that seek to uncover the mechanisms by which variants contribute to disease phenotypes. Unlike functional prioritization of non-coding variants task which only filters the variants that are most likely to have functional significance. Prediction of context-specific functional impact of genetic variants provides a finer level of detail by predicting the actual effect of a variant on gene expression or other functional outcomes in specific tissues. This granularity is essential for precisely understanding the specific biological mechanisms and for developing targeted therapies (36).

Exon and intron region classification is crucial for understanding gene structure and function within the genome. Exons are coding regions that are translated into proteins, while introns are non-coding regions that are spliced out during mRNA processing. By classifying exons and introns, researchers can describe gene boundaries, identify functional elements, and elucidate the mechanisms of gene expression regulation (166). This task is essential for deciphering the genetic code and unraveling the complexities of gene regulation in health and disease. Recombination spots identification is a pivotal task in functional genomics that focuses on mapping regions of the genome where genetic recombination events occur. Genetic recombination is a natural process where DNA segments are exchanged between two chromosomes during cell division. Recombination plays a vital role in generating genetic diversity, ensuring proper chromosome segregation, and driving evolution (167). By identifying recombination hot spots, researchers can gain insights into the mechanisms underlying genetic diversity and genome evolution, shedding light on the processes that shape genetic variation and adaptation in populations. In conclusion, the tasks related to functional genomics, including species classification, conserved non-coding element classification, functional prioritization of non-coding variant, prediction of context-specific functional impact of genetic variants, exon and intron region classification, and recombination spots identification, are essential for advancing our understanding of genetic regulation, molecular mechanisms, and disease pathogenesis. By delving into these tasks, researchers can unravel the complexities of the genome, decipher the genetic basis of diseases, and pave the way for precision medicine and personalized healthcare interventions tailored to an individual's genetic profile.

Another goal of researchers is to study overlap between two distinct fields namely environmental science and microbial genomics (27). This interdisciplinary study enables researchers to explore how environmental factors such as pollution, climate change, and agricultural practices affect on function and diversity of microbial communities (27). A key area of focus in this field is the nitrogen cycle prediction. By examining the genomes of microbes involved in nitrogen fixation, nitrification, and denitrification, scientists can predict how these processes might respond to environmental changes (168). This prediction provides understanding about potential impacts of environmental shifts on ecosystem health (169) and nitrogen availability, which are essential for plant growth and overall biogeochemical cycles (170).

From all eight different biological goals, disease analysis goal has received huge attention in scientific community as it aims to understand, diagnose, and treat various illnesses. Within this field, several tasks play a vital role in enhancing our comprehension of diseases. One such task is Pathogen Signatures Identification (171), which involves identifying specific markers or characteristics of pathogens that can aid in their detection and classification (172). By pinpointing these signatures, researchers can develop targeted diagnostic tools and therapies, ultimately improving disease management and control. Mutation Susceptibility Analysis (173) is another essential task in disease analysis. This task focuses on investigating the genetic variations that make individuals more prone to developing certain diseases (174). Understanding mutation susceptibility can aid in personalized medicine approaches, where individuals at higher risk can be identified early for preventive interventions or closer monitoring. Phage-Host Interactions Prediction (175177) is a task that delves into the relationships between bacteriophages and their host bacteria (178). By predicting these interactions, researchers can gain insights into how phages influence bacterial populations, which is crucial for developing phage-based therapies to combat bacterial infections and antibiotic resistance. Disease Risks Estimation (90) is a fundamental aspect of disease analysis that involves assessing the likelihood of an individual developing a particular condition based on various factors such as genetics, lifestyle, and environmental exposures (179). Accurately estimating disease risks enables healthcare providers to offer targeted interventions and counseling to high-risk individuals, potentially preventing the onset or progression of diseases. Tumor Type Prediction (180) is a significant task in disease analysis that focuses on identifying the specific type of tumor a patient may have based on various characteristics such as genetic markers, imaging features, and histopathological findings (181). Predicting tumor types is essential for determining the most effective treatment strategies and prognostic outcomes for patients with cancer. Pathogenicity Potential Assessment (27) is a critical task that involves evaluating the ability of pathogens to cause disease in a host. By assessing the pathogenicity potential of different microorganisms, researchers can prioritize the development of interventions against the most virulent pathogens, thereby improving disease prevention and control strategies. Phylogenetic Analysis () is a key component of disease analysis that involves studying the evolutionary relationships between different strains of pathogens or tumor cells. Phylogenetic analysis provides insights into the origins, spread, and diversification of diseases, aiding in the development of targeted interventions and understanding disease transmission dynamics.

4 A look on DNA sequence analysis tasks from the perspective of computer scientists

While Section 3 presents a biologically motivated categorization of DNA sequence analysis tasks, this section reframes these same tasks from a computational perspective. This dual categorization approach (biological and computational) aims to facilitate interdisciplinary understanding between life scientists and AI researchers. With the influx of biological data and rise of AI, researchers are increasingly applying AI in diverse areas of molecular biology. Development of large scale AI applications requires a good understanding of variety of sequence analysis tasks. However, there exist a huge domain gap between computer scientists and molecular biologists. Molecular biologists know the need, biological importance, and pharmaceutical worth of different sequence analysis tasks. However, they do not know which machine or deep learning models are most appropriate to use to either replace or complement experimental work. Similarly, computer scientists know which Artificial Intelligence predictive pipeline can potentially perform better with specific type of data; however, they do not know the nature of biological sequence analysis tasks. For instance, DNA sequence analysis tasks such as gene function prediction, gene network reconstruction, gene expression prediction, and disease risk estimation can be challenging for computer scientists to grasp. However, a comprehensive literature review that explains the basics of these tasks can bridge this gap. For example, gene function prediction is a multi-label classification tasks, gene expression prediction is a regression task, while gene network reconstruction and disease risk estimation are binary classification tasks. With this foundational understanding, computer scientists can more easily develop predictive pipelines for these binary, multi-label classification, and regression tasks. To empower AI experts, we have presented 44 DNA sequence analysis tasks in computer scientist language in Figure 4. A simple look on Figure 4 reveals that nature of DNA sequence analysis tasks can be categorized into three primary types: regression, clustering, and classification where classification can be further divided into three secondary types: binary classification, multi-class classification, and multi-label classification. Let us mathematically formulate the possible natures of DNA sequence analysis tasks.

Figure 4

In binary classification, researchers aim to predict the outcome of a binary variable (0 or 1). Given a dataset with features X ∈ ℝnxd, binary labels y ∈ 0, 1, and training dataset (x1, y1), (x2, y2), …, our goal is to learn a decision function f:XY that maps inputs to binary outputs 0, 1 on the basis of hypothesis function h(x) learned from the training data.

In multi-class classification, researchers aim to predict the outcome from more than two classes. Specifically, given a dataset having sequences X ∈ ℝnxd, labels y ∈ 1, 2, …, K where K is the number of classes, and training dataset (x1, y1), (x2, y2), …, (xn, yn) where xiX and yiY, our goal is to learn a decision function f:XY that assigns inputs to one of the classes.

where hk(x) is the hypothesis function for class k learned from the training data. On the other hand, in multi-label classification, each input can be assigned to multiple classes simultaneously. Given a dataset with features X ∈ ℝnxd, labels y ∈ 1, 2, …, K where K is the number of classes, and training dataset (x1, y1, y2, ..), (x2, y1, y4, …), …, (xn, y5, yn, ….) where xiX and yiY, our goal is to learn a decision function f:X → 0, 1K that assigns inputs to multiple classes simultaneously using hypothesis function hk(x) for class k learned from the training data.

Furthermore, in regression, researchers goal is to predict a continuous outcome variable. Given a dataset with sequences X ∈ ℝnxd, labels y ∈ ℝ, and training dataset (x1, y1), (x2, y2), …, (xn, yn) where xiX and yiY, our goal is to learn a function f:X → ℝ that predicts continuous outputs using hypothesis function h(x) learned from the training data.

In clustering, the goal is to group similar data points into same clusters. Given a dataset with data points X = x1, x2, …, xn, where each , our goal is to find a partition of the data into clusters C = C1, C2, …, CK. This is done on the basis of a distance metric d(x, μc) between data point x and the centroid μc of cluster c.

5 DNA sequence analysis databases

This section provides a comprehensive overview of various databases employed to develop benchmark datasets for development of AI-based applications for 44 distinct DNA sequence analysis tasks. A total of 45 DNA sequence databases have been identified from 127 existing studies. Among these, 36 databases are publicly accessible, while the remaining 9 databases are either inaccessible or no longer exist. To ease the lives of researchers and practitioners, Table 1 summarizes accessible databases in terms of their release year, types of inherent genetic data (DNA, RNA, protein), details of species and organisms, statistics of raw sequences, and supported data formats.

Table 1

Database nameRelease dateType of data (DNA/RNA/ Protein)Organism nameSpecies nameSequences statisticsData format
Descartes2020DNA, RNAAnimalsMus musculus, Homo sapiensHuman Gene Expression During Development: 4M Cells, 121 Tissues, 15 Organs; Human Chromatin Accessibility During Development: 720K Cells, 53 Tissues, 15 Organs; Mouse: ~2M Cells, 61 Embryos.RDS
PPD2020DNA, RNABacteria, Archaea63 species129,148 Promoter Sequences, 63 Species with 74 strains.csv
EnhancerAtlas 2.02019DNAAnimals, BacteriaHomo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Danio rerio, Rattus norvegicus, Gallus gallus, Sus scrofa, Saccharomyces cerevisiae13,494,603 Enhancers, 586 tissue.csv, BED, R Object
DREAM Base2018DNA, RNA, ProteinAnimalsHomo sapiensscRNA-Seq Data: 93,3704 cells; RNA-Seq Data: ~18,196; ChiP-Seq Data: ~10,000; RNA Modification Data: ~500; ribo-Seq Data: 1,570; DNase-Seq Data: 599; CLIP-Seq Data: 568.excel, .txt
EmExplorer database2018DNA, ProteinAnimalsBos taurus, Homo sapiens, Mus musculus, Rattus norvegicus, Sus scrofa158,000 items that contain more than 32,000 development-related Genes under 306 related pathways.txt
COSMIC2018DNAAnimalsHomo sapiensTotal Genomic variants = 24,599,940; Genomic non-coding variants = 16,748,366,406; Genomic mutations within Exons = 768; Genomic mutations within Intronic and other intragenic regions = 9,217,664; Samples = 1,531,613; Fusions = 19,428; Gene expression variants = 9,215,470; Differentially Methylated CpGs = 7,930,489.FASTA, .tsv
DisGeNet2015DNAAnimalsHomo sapiens1,134,942 GDAs between 21,671 Genes, 30,170 diseases and traits; 369,554 VDAs between 194,515 variants and 14,155 diseases and traits.txt, RDF, SQL Dump
genomAD2014DNAAnimalsHomo sapiens730,947 Exomes, 76,215 whole GenomesVCF, Hail Table
ClinVar2013DNA, RNA, ProteinAnimalsHomo sapiensRecords = 4391341, Total Genes = 92225.xml, VCF, .tsv
HOCOMOCO Human v11 database2013DNAAnimalsMus musculus, Homo sapiens1,443 TF binding models including secondary motif subtypes for 949 human TFs and 720 Mouse orthologsPWM, PFM, PCM, Flat text files
DeOri2012DNAAnimalsHomo sapiens, Mus musculus, Arabidopsis thaliana, Kluyveromyces lactis, Schizosaccharomyces pombe, Drosophila melanogaster189,743 entries.FASTA
BioLip2012DNA, RNA, ProteinAnimalsHomo sapiens873,925 Entries, 448,816 regular ligands, 191,485 mental ligands, 37,492 Peptide ligands, 43,448 DNA ligands, 152,684 RNA ligands, 873,925 binding affinity data, 451,485 Protein receptors.FASTA
DeOri6.02011DNAAnimals, Plants, Fungi17 species189,740 eukaryotic replication origins.FASTA
GWAS2008DNAAnimalsHomo sapiens146,394 TASs.tsv, OWL/RD
Broad DepMap2008DNAAnimalsHomo sapiens2,000 Human cancer cell lines.csv, .txt
CCLE2008DNA, RNA, ProteinAnimalsHomo sapiens1,019 RNA cell lines, 954 microRNA expression profiles, 899 Protein lines, 897 Genome-wide histone modifications, 843 DNA methylation, 329 whole Genome Sequencing, 326 whole exome Sequencing.csv
GENCODE2006DNAAnimalsHomo sapiens, Mus musculusHomo sapiens: Total genes = 63,086, Total transcripts = 254,070, Total distinct Translations = 65,650; Mus musculus: Total Genes = 57,132, Total Transcripts = 149,138, Total distinct Translations = 44,819.txt
Consensus Coding Sequence Database2005DNAAnimalsHomo sapiens, Mus musculus35,608 CCDS IDs that correspond to 19,107 Genes, with 48,062 Protein Sequences.FASTA
MSigDB2005DNAAnimalsHomo sapiens, Mus musculus8,380 Gene set.gct, .res, .pcl, .txt, .cls, .gmx, .gmt, .grp, .xml, .chip, .rnk
Gene Ontology2004DNAAnimals, Bacteria, Fungi, PlantsEscherichia coli, Homo sapiens, Oryza sativa, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Mus musculus, many moreAnnotated Gene products = 1,536,921; Annotated Species = 5,409; Annotated Species with over 1,000 annotations = 183OBO, OWL, GAF, GPAD, GPI, JSON
JASPAR2004DNA, RNA, ProteinFungi, Insects, Nematoda, Plants, Urochordata, Vertebrata34 species4279 Profiles.txt
Database of Essential Genes2004DNA, RNA, ProteinBacteria, Archaea, Eukaryotes51 species53,885 essential Genes, 786 essential non-coding Sequences.csv, DAT
ENCODE2003DNA, RNA, ProteinAnimalsHomo sapiens, Mus musculus_.FASTA, BAM, BigWig, BED, VCF
DataBase of Transcriptional Start Sites2002DNAAnimalsHomo sapiens, Mus musculus491M TSS tag Sequences from a total of 20 tissues and 7 cell cultures.FASTA, .csv, .xlsx, .txt
MGC2001DNAAnimalsMus musculus, Homo sapiens, Rat, BovineTotal MGC Full ORF Clones: Homo sapiens: 29,818; Mus musculus: 27,285; Rat: 6,763; Bovine: 9,104; Non-redundant Genes: Homo sapiens: 17,592; Mus musculus: 17,701; Rat: 6,486; Bovine: 8,724_
GEO2000DNA, RNA, ProteinAnimals20 speciesSamples = 7,209,691SOFT, MINiML, .txt
Exon-Intron Database1999DNA, RNAAnimals, PlantsHomo sapiens, Mus musculus, Drosophila melanogaster, C. Elegans, S. Pombe42,460 Genes (243,589 Exons).FASTA
Ensembl1999DNA, RNA, ProteinAnimalsHomo sapiens, Mus musculus, Danio rerio, Sus scrofaGenomes = 44,048, Ensembl Fungi = 1,014 Genomes, Ensembl Metazoa = 78 Genomes (Invertebrate species), Genomes for vertebrate Species = 236, Ensembl Plants = 67 Genomes, Ensembl Protists = 237 Genomes.FASTA, GTF, GFF, MySQL Dump
RegulonDB1998DNA, RNA, ProteinBacteriaEscherichia coli4,748 Genes, 2,590 Operons, 287 Regulons, 3,718 Transcription Unit, 4050 Promoters.tsv, .csv
EPD21998DNAAnimals, Plant, Fungi, Protists139 species4,806 Promoters.FASTA, EMBL
dbSNP1998DNAHumanHomo sapiensNearly 2 billion submissions representing more than 675 million distinct variants; 23.7 million refSNP entries (14.5 million validated)ASN.1, FASTA, XML
KEGG1995DNA, RNA, ProteinAnimals, Plants, Fungi, Protists, Bacteria, ArchaeaEuryarchaeota Candidatus, Thermoplasmatota, Thermoproteota, Chordata, Echinodermata, Hemichordata, Ascomycota, Basidiomycota, Atribacterota Candidatus, SaccharibacteriaGenes = 53,674,741, Addendum Proteins = 4,181, Viral Genes=688,823, Viral mature Peptides = 377KGML, .FASTA, .txt
NCBI1988DNA, RNA, ProteinMultiple organismsMultiple speciesHub of databases.FASTA, XML
Eukaryotic Promoter Database1986DNAAnimals, Plants, Fungi, InvertebratesHomo sapiens, Macaca mulatta, Mus musculus, Rattus norvegicus, Gallus gallus, Canis familiaris, Drosophila melanogaster, Apis mellifera, Danio rerio, Caenorhabditis elegans, Arabidopsis thaliana; Zea mays, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Plasmodium falciparum192,586 Promoters, 163,676 Genes.FASTA, EMBL
GenBank1982DNAAnimals, Archaea, Bacteria, Fungi, Plants, Virus557000 species3,213,818,003,787 Bases, 250,803,006 Sequences.gb
OMIM1960DNAAnimalsHomo sapiens17,290 Gene descriptions, 18 Gene and Phenotypes, 8,361 Phenotype description, 1,736 Phenotypes suspected mendelian.txt

Summary of publicly accessible biological databases, their inherent data types, species diversity, and statistics of raw sequences related to different genomic and proteomic data.

A holistic view of the Table 1 reveals that 12 databases provide RNA and protein sequences as well in addition to providing DNA sequences. As word embeddings methods and large language models are trained in unsupervised fashion and when they are trained on large sequence data usually, they produce better representations. To efficiently train word embedding methods and large language models, raw data can be acquired from these databases. To facilitate researchers, we have categorized 36 databases into three different categories on the basis of volume of raw sequences: low sequence facilitators, medium sequence facilitators, and high sequence facilitators. Specifically, 13 low sequence facilitators, namely, HOCOMOCO Human v11 database (182), Consensus Coding Sequence Database (183), MSigDB (184), Broad DepMap (185), JASPAR (186), Database of Essential Genes (187), ENCODE (188), MGC (189), Exon-Intron Database (190), Ensembl (191), RegulonDB (192), EPD2 (193), offer up to 100,000 DNA sequences each, while 9 medium sequence facilitators, namely, PPD (194), DREAM (195), EmExplorer database (196), GenomAD (197), DeOri (198), BioLip (199), DeOri6.0 (198), GWAS (200), Eukaryotic Promoter Database (193), provide up to 1 million DNA sequences. In contrast, 13 high sequence facilitators such as Descartes (201), EnhancerAtlas 2.0 (202), COSMIC (203), DisGeNet (204), ClinVar (205), CCLE (206), GENCODE (207), Gene Ontology (208), DataBase of Transcriptional Start Sites (209), GEO (210), KEGG (211), NCBI (212), GenBank (213), and dbSNP (214, 215) offer more than 1 million DNA sequences each. These databases predominantly house DNA sequences from a diverse array of species, including humans, mice, plants, bacteria, and fungi. A comprehensive analysis reveals that approximately 22 databases, namely, Descartes (201), DREAM (195), EmExplorer database (196), COSMIC (203), DisGeNet (204), GenomAD (197), ClinVar (205), HOCOMOCO Human v11 database (182), DeOri (198), BioLip (199), GWAS (200), Broad DepMap (185), CCLE (206), GENCODE (207), Consensus Coding Sequence Database (183), MSigDB (184), ENCODE (188), DataBase of Transcriptional Start Sites (209), MGC (189), GEO (210), Ensembl (191), and OMIM (216), focus on animal DNA sequences, 4 databases including PPD (194), Database of Essential Genes (187), RegulonDB and (192) on bacterial sequences, and JASPAR (186) on plant DNA sequences. EnhancerAtlas 2.0 (202) is the only database that facilitates with both animal and bacterial DNA sequences, while 4 databases namely DeOri6.0 (198), Exon-Intron Database (190), EPD2 (193), and Eukaryotic Promoter Database (193) focus on animal and plant DNA sequences, whereas Gene Ontology (208), KEGG (211), and GenBank (213) provide DNA sequences for animal, plant, and bacteria. In addition, sequences from other organisms such as eukaryotes, invertebrates, fungi, and various microorganisms are also well-represented. Some databases encompass a broad spectrum of species. For instance, the EDP2 (193) database includes genomics data for 139 species, GenBank (213) houses sequences for 557,000 species, and PPD (194) has genomics data of 63 species.

Moreover, Table 1 includes data formats utilized by various databases to manage and provide access to DNA sequences. TXT and FASTA format are universally accepted by almost all DNA sequence analysis programs. Each entry in both format types contains at least two lines: First line or header includes accession number, species name, or identification details, while next line contains nucleotide sequences. CSV and TSV are text-based formats in which values in rows are separated by commas or tabs, respectively. In both file formats, first row specifies headers which defines names of columns (“SeqID”, “SeqName”, “Type”, “Function”) and subsequent rows represent data. In VCF format, first row specifies headers which defines names of columns, but this format is specifically used to store genetic variation data including single nucleotide polymorphisms (SNPs), insertions, deletions, and structural variants. In addition, XLSX formats represent complex datasets that contain information computed with various formulas across multiple columns, whereas EMBL format includes structured sections for sequence data, feature annotations (genes and other biological features), organism information, references, and other details. An extensive analysis of Table 1 reveals that most widely used data formats are FASTA, TXT, CSV, XLSX, and EMBL in DNA sequence analysis.

A rigorous analysis of Table 1 reveals that out of 36 publicly accessible databases, several key categories of data emerge. Four databases, namely, Broad DepMap (185), genomAD, COSMIC, and MGC, provide data for DNA functional analysis tasks such as prediction of context-specific functional impact of genetic variants and conserved non-coding element classification. Seven databases, namely, BioLip, HOCOMOCO Human v11, GWAS, EnhancerAtlas 2.0, DataBase of Transcriptional Start Sites, Exon-Intron Database, and Eukaryotic Promoter Database, offer data on gene expression regulation. Three databases, namely, PPD, CCLE, and EmExplorer, focus on DNA modification data including methylcytosine and methyladenine modifications. In addition, DeOri, Descartes, DeOri6.0, and JASPAR provide information on gene structure and stability, including chromatin accessibility prediction, YY1-mediated chromatin loop identification, and DNA replication origins identification. GENCODE, Consensus Coding Sequence Database, MSigDB, Gene Ontology, DisGeNet, Database of Essential Genes, KEGG, and NCBI offer comprehensive gene analysis data. Furthermore, eight other databases, namely, EPD, ENCODE, RegulonDB, GEO, Ensembl, ClinVar, GenBank, and OMIM, provide a range of data on gene expression regulation, DNA modification prediction, genome structure and stability, DNA functional analysis, disease information, and gene analysis.

6 DNA sequence analysis benchmark datasets

The quality and quantity of datasets utilized in AI-driven DNA sequence analysis applications are vital determinants of their effectiveness and functionality. This section aims to provide a comprehensive overview of datasets relevant to 44 distinct DNA sequence analysis tasks. Overall, these datasets fall into two primary categories: publicly available datasets and in-house datasets. This categorization serves to illuminate the significance of dataset accessibility and its implications for the advancement of AI-driven DNA sequence analysis. Specifically, publicly available datasets are accessible to the wider research community and are commonly employed in the development of AI-based predictive models. They serve as foundational resources that facilitate the advancement of AI-driven DNA sequence analysis pipelines by ensuring accessibility, reusability, and transparency in research endeavors. Furthermore, the utilization of publicly available datasets fosters collaboration and knowledge exchange within the scientific community, thereby contributing to the overall progress of the field. In contrast, in-house datasets are proprietary in nature and are developed within specific research laboratories or institutions. These datasets often contain sensitive data tailored to particular research objectives. As in-house datasets cannot be shared publicly, their proprietary nature may limit broader access, reproducibility, and applicability of findings.

Rigorous assessment of 127 existing studies reveals that a total of 242 benchmark datasets related to 44 distinct DNA sequence analysis tasks are constructed or acquired from existing literature. Specifically, among these 242 benchmark datasets, 199 are publicly available and 43 are in-house datasets. Table 2 provides the distribution of public and in-house datasets for 44 distinct DNA sequence analysis tasks. It provides information about which of these datasets are used by word embeddings, large language models, nucleotide composition, and positional information-based predictive pipelines.

Table 2

Task nameTask typeDatasets used in language modelsDatasets used in word embeddingsDatasets used in other methods
PublicIn-housePublicIn-housePublicIn-house
DNA replication origins identificationBinary classification_Gao et al. (A. thaliana) (25)_Wu et al. datasets (S. cerevisiae Dataset, S. pombe Dataset, K. lactis Dataset, P. pastoris Dataset) (329)__
Nucleosome position detectionBinary classificationGangi et al. datasets (CE, DM, YS, HM, DM-5U, DM-PM, DM-LC, HM-5U, HM-LC, HM-PM, YS-PM) (330)_____
Chromatin accessibility predictionBinary classificationDeepSEA datasets (TF, DHS) (32), DNase-Seq experiment data (31)_____
YY1-Mediated chromatin loops predictionBinary classification__Dao et al. DeepYY1 datasets (HCT116, K562) (39)_Zhang et al. DeepYY1 datasets (HCT116, K562) (38)_
Genome structure analysisMulti-class classificationNCycDB dataset (27)_____
Chromatin feature predictionMulti-label classificationLogo919 (34), Logo2002 (34), Logo3357 (34)_____
Long-range chromatin interaction predictionInteractionChip-seq dataset (36)_____
Enhancers identificationBinary classificationLiu et al. dataset (57), Liao et al. datasets (HEK293, NHEK, K652, GM12878, HMEC, HSMM, NHLF, HUVEC) (58)_Liu et al. dataset (57)_DiseaseEnhancer (55), EnDisease (55), CancerEnD (55)_
Promoter identificationBinary classificationYang et al. dataset (34), Ji et al. dataset (90), Xiao et al. dataset (331)_Wang et al. dataset [H. Sapiens-I (TATA-containing), H. Sapiens-II (TATA-less), R. Norvegicus-I (TATA-containing), R. Norvegicus-II (TATA-less), D. melanogaster-I (TATA-containing), D. melanogaster-II (TATA-less), Z. mays-I (TATA-containing), Z. mays-II (TATA-less)] (236), Zhang et al. dataset (K562, GM12878, HeLa-S3, HUVEC) (76), Xiao et al. dataset (331)_Yang et al. dataset (34), Xiao et al. (331)_
Enhancer-promoter interactions predictionInteraction/binary classificationYang et al. datasets (FoeT, Mon, nCD4, tB, tCD4, tCD8) (249)_Whalen et al. dataset (GM12878, HUVEC, HeLa-S3, IMR90, K562, NHEK) (329)_Whalen et al. dataset (GM12878, HUVEC, HeLa-S3, IMR90, K562, NHEK) (329)Zhang et al. (GM12878 cell line, HeLa cell line) (332)
Transcription sites predictionBinary classificationClauwaert et al. dataset (50)_____
Transcription factor binding sites predictionBinary classificationChIP-Seq dataset (94), TSSs dataset (91), 497 TF ChIP-Seq dataset (90), 690 ChIP-Seq (51)_Shen et al. datasets (A549 dataset, MCF-7 Dataset, H1-HESCDataset, HUVEC dataset) (92)___
Transcription factor binding affinity predictionMulti-class classificationWeirauch et al. dataset (PBM Dataset), Jolma et al. dataset (HT-SELEX Dataset) (52)_____
Protein-DNA binding sites predictionInteraction/binary classification690 ChIP-Seq dataset (96), Patiyal et al. dataset (53), Xia et al. (dataset 2) (53), Liu and Tian (dataset 1, dataset 2) (95)_____
Splice sites predictionBinary classificationWang et al. dataset (100), Ji et al. dataset (93)___Splice-junction gene sequence dataset (97), Degroeve et al. (98), Liu et al. Datasets [O. sativa (Acceptor, Donor), A. Thaliana (Acceptor, Donor), H. sapiens (Acceptor, Donor)] (99)_
Translation initiation sitesBinary classificationClauwaert et al. dataset (50)___Kalkatawi et al. TIS (101)_
Essential genes identificationBinary classificationMa et al. datasets (S. cerevisiae, E. coli, H. sapiens, D. melanogaster) (109)_Ma et al. datasets (S. cerevisiae, E. coli, H. sapiens, D. melanogaster) (109), Campos et al. datasets (D. melanogaster, M. maripaludis, H. sapiens, C. elegans) (333), Zhang et al. (106), Xiao et al. dataset (107)_Campos et al. datasets (D. melanogaster, M. maripaludis, H. sapiens, C. elegans) (333), Sharma et al. (104)_
Disease genes predictionBinary classificationNunes et al. dataset (110)_____
Pseudogene function predictionInteraction/binary classification__Fan et al. dataset (CC, MF, BP) (111)___
Target gene classificationMulti-class classification__Arango et al. dataset (113)___
Candidate gene prioritization & SelectionMulti-class classification_Toufiq et al. dataset (114)____
Gene functions predictionMulti-label classificationHu et al. (Human gene set from gene ontology) (112)_GTEx (CC, MF, BP) (111)___
Gene expression predictionRegressionReddy et al. dataset (Jurkat, K-562, THP-1) (334)_Al Taweraqi et al., dataset (103)___
Gene taxonomy classificationMulti-class classificationMock et al. dataset (123)_Verma et al. dataset (121)CAMI2 Airway dataset (122)__
Gene network reconstructionMulti-class classification__SynTReN dataset (128) DREAM5 dataset (128)__Pio et al. dataset (126), Schaffter et al. datasets (DREAM4 10, DREAM4 100) (127), Jozefczuk et al. datasets (E.coli cold, E.coli head, E.coli oxidative) (127)
4mc-Methyl-cytosine modification predictionBinary classificationXu et al. datasets (C. elegans, D. Malenogaster, A. Thaliana, E.coli, G. subterraneus, G. pickeringii) (141), Chen et al., datasets (C. elegans, D. Malenogaster, A. Thaliana, E.coli, G. subterraneus, G. pickeringii) (135)_Khanal et al. datasets (F. vesca, R. chinensis) (138), Zulfiqar et al. dataset (136), Zeng et al. dataset (137)_Khanal et al. (138), Chen et al., Datasets (C. elegans, D. Malenogaster, A. Thaliana, E.coli, G. subterraneus, G. pickeringii) (135)Manavalan et al. (139)
5mc-Methyl-Cytosine modification predictionBinary classificationWang et al. dataset (282), Stanojevic et al. dataset (GM24385, NA12878, NA19240, H1ESc, K562, HX1) (152)_Wang et al. dataset (282), Hyb_2021 (C. elegans, D. Malenogaster, A. Thaliana, E.coli, G. subterraneus, G. pickeringii) (142)_Xu et al. (141), Rao Zeng et al. (140), Saha et al. (135)Nguyen-Vo et al. (139)
5hmc-hydroxy-methylcytosine modification predictionBinary classification__Lv et al. dataset (M. musculus, H. sapiens) (156)___
6mA-methyladenine modification predictionBinary classificationAbbas et al. dataset (A. thaliana, H. sapiens, M. musculus, S. cerevisiae) (335), DNA 6 mA dataset (281)_Lv et al. dataset (156)_Zhou et al. dataset (A. thaliana, C. elegans, C. equisetifolia, D. melanogaster, F. vesca, H. sapiens, R. chinensis, S. cerevisiae, T. thermophile, Xos. BLS256) (148), Fan et al. dataset (1, D.malenogaster, 3, 4, 5) (149)_
Methylation modification predictionBinary classificationLv et al. dataset 6mA (T. thermophile, A. thaliana, H. sapiens, Xos. BLS256, D. melanogaster, C. elegans, C. equisetifolia, S. cerevisiae, Tolypocladium, F. vesca, R. chinensis) 5hmC (M. musculus, H. sapiens) 4mC (F. vesca, Tolypcladium, S. cerevisiae, C. equisetifolia) (156)_____
Conserved non-coding element classificationBinary classification___Polychronopoulos et al. dataset (163)__
Functional prioritization of non-coding variantsMulti-class classificationLogo919 (34) Logo2002 (34) Logo3357 (34)_____
Exon & Intron region classificationBinary classificationAkalin et al. dataset (164)_____
Recombination spots identificationBinary classification__Liu et al. dataset (165)___
Species classificationMulti-class classificationMouse enhancers (44) Coding vs. intergenomic (44) human vs. worm (44) Human enhancers cohn (44) human enhancers ensembl (44) Human Regulatory (44) human nontata promoter (44) human OCR ensembl ()_____
Prediction of context-specific functional impact of genetic variantsMulti-class classificationeQTLs dataset (36)_____
Nitrogen cycle predictionMulti-class classification_NCycDB (27)____
Pathogen signature identificationBinary classification___DS500 dataset (171)__
Phage-host interactions predictionInteraction/binary classification__ESKAPE dataset (177), Wang et al. dataset (176)_Qiu et al. dataset (Kingdom, Phylum, Class, Order, Family, Genus) (175)_
Mutation susceptibility analysisMulti-class classification___Yilmaz et al. dataset (Human, Mouse) (173)__
Tumor type predictionMulti-class classificationTCGA pan-cancer dataset (180)_____
Pathogenicity potential assessmentMulti-class classification_E-K12 (27), CARD-A (27), CARD-D (27), CARD-R (27), VFDB (27), ENZYME (27), PATRIC (27), NCycDB (27)____
Phylogenetic analysisClustering___Ren et al. datasets (111)__
Disease risks estimationBinary classificationHSCR-RET (90), HSCR-RET-Long (90)_____

Overview of 199 public and 43 in-house datasets used across 44 different DNA sequence analysis tasks.

For each DNA sequences analysis task, public and in-house datasets are distributed as DNA Replication Origins Identification (0, 5), Nucleosome Position Detection (11, 0), Chromatin Accessibility Prediction (2, 0), YY1-Mediated Chromatin Loop Prediction (4, 0), Genome structure analysis (0, 1), Chromatin Feature Prediction (3, 0), Long-range chromatin interaction prediction (1, 0), Enhancers Identification (12, 0), Promoter Identification (15, 0), Enhancer-Promoter Interactions Prediction (18, 2), Transcription Site Prediction (1, 0), Transcription Factor Binding Site Prediction (4, 4), Transcription Factor Binding Affinity Prediction (2, 0), Protein-DNA Binding Site Prediction (5, 0), Splice Site Prediction (10, 0), Translation Initiation Sites (1, 1), Essential Gene Identification (6, 5), Disease Gene Prediction (1, 0), Pseudogene Function Prediction (3, 0), Target Gene Classification (1, 0), Candidate Gene Prioritization/ Identification (0, 1), Gene Function Prediction (4, 0), Gene Expression Prediction (4, 0), Gene Taxonomy Classification (2, 1), Gene Network Reconstruction (2, 6), 4mc-Methylcytosine Site Prediction (16, 0), 6mA-Methyladenine Site Prediction (5, 0), 5mc-Methylcytosine Site Prediction (24, 1), 5hmc-Methylcytosine Site Prediction (2, 0), Methylation Site Prediction (17, 0), Conserved Non-Coding Elements Classification (0, 1), Functional Priorizitation of non-coding variants (3, 0), Exon and Intron Region Classification (0, 1), Recombination Spots Identification (1, 0), Species Classification (8, 0), Prediction of context-specific functional impact of genetic variant (1, 0), Nitrogen Cycle Prediction (0, 1), Pathogen Signatures Identification (0, 1), Phage-Host Interactions Prediction (8, 0), Mutation Susceptibility Analysis (0, 2), Tumor Type Prediction (1, 0), Pathogenicity Potential Assessment (0, 8), Phylogenetic Analysis (0, 1), and Disease Risks Estimation (2, 0). First entry in brackets refers to count of public datasets, and second entry indicates total number of in-house datasets for a particular task. For example, in “Essential Gene Identification (6, 5)” task, 6 refers to public datasets while 5 represents in-house datasets.

A holistic view of Table 2 reveals 110 public and 18 in-house datasets are employed to develop both word embeddings and language models based predictive pipelines for 12 DNA sequence analysis tasks, namely, DNA replication origins identification, enhancers identification, promoters identification, enhancer-promoter interaction prediction, transcription factor binding site prediction, essential gene identification, gene function prediction, gene expression prediction, gene taxonomy classification, 4mC-methyl cytosine modification prediction, 5mC-methl cytosine modification, and 6mA-methyl modification prediction. Notably, both types of predictive pipelines have utilized 1 common dataset to evaluate the performance of predictive models developed for three tasks, namely, enhancer identification, essential gene identification, and 5mC-methyl cytosine modification prediction.

Furthermore, 112 public and 15 in-house datasets are used to develop both word embedding and nucleotide compositional and positional information-based predictive pipelines for 11 DNA sequence analysis tasks including essential gene identification, gene network reconstruction, 4mC-methyl cytosine modification prediction, 5mC-modification prediction, 6mA-methyl adenine modification prediction, and phage-host interaction prediction. However, both predictive pipelines have used 9 common pubic dataset for only three tasks. Specifically, six public datasets for enhancer-promoter interactions prediction, one public data for essential gene identification, and two public datasets for 4mC-Methyl cytosine modification prediction are commonly employed by both predictive pipelines.

Moreover, Table 2 highlights that 107 public and 9 in-house datasets are utilized by 9 DNA sequence analysis tasks, namely, enhancers identification, promoters identification, enhancer-promoter interaction prediction, splice site prediction, translation initiation sites identification, essential gene identification, 4mC-methyl cytosine modification prediction, 5mC-methl cytosine modification, and 6mA-methyl modification prediction for developing both language models and nucleotide compositional and positional information-based predictive pipelines. Merely, 7 public datasets are used commonly by both predictive pipelines for two tasks: one for promoter identification and six for 4mC-methyl cytosine modification prediction.

Although all three different types of representation learning-based predictive pipelines are employed across six different DNA sequence analysis tasks, namely, enhancers identification, promoters identification, enhancer-promoter interaction prediction 4mC-methyl cytosine modification prediction, 5mC-methl cytosine modification, and 6mA-methyl modification prediction, only on one task, namely, promoters identification, all three kinds of predictive pipelines are evaluated on one common dataset. These statistics reveal that researchers have focused on creating new datasets for each kind of predictive pipelines instead of using existing datasets. Consequently, this domain lacks a fair performance comparison between different kinds of predictive pipelines.

7 A brief look on representation learning methods and predictors used in DNA sequence analysis predictive pipelines

This section dives into 12 most commonly used word embedding approaches, 8 large language models, 9 machine learning, 8 deep learning, and 3 statistical algorithms that are used in development of predictive pipelines for 44 different DNA sequence analysis tasks.

7.1 DNA sequence representation learning using word embeddings

In the domain of natural language processing (NLP), the introduction of word embedding techniques represented a significant advancement by enabling the development of more accurate machine and deep learning predictive models. These approaches assign statistical vectors to words by capturing contextual representations of words within extensive, unlabelled corpora (217, 218). The primary objective is to assign comparable vectors to semantically similar words and distinct vectors to dissimilar words (217, 218). Leveraging transfer learning strategies, these contextual word representations have empowered data-hungry deep learning models to achieve exceptional performance, even with limited training data. Following the success of word embeddings in various NLP tasks (217220), researchers have adopted these approaches for genomic and proteomic sequence analysis tasks, which share similarities with NLP tasks. This section offers a comprehensive overview of 12 distinct word embedding approaches that are utilized in DNA sequence analysis predictive pipelines. Figure 5 visually illustrates the utilization of various word embedding methods in conjunction with different machine and deep learning algorithms.

Figure 5

These word embedding approaches leveraged for DNA sequence analysis tasks can be categorized into two types: (1) non-graph-based methods and (2) graph-based methods. Non-graph-based methods segregate DNA sequences into overlapping or non-overlapping k-mers. Specifically, overlapping k-mers are generated by sliding a fixed-size window over sequence with a smaller stride as compared to window size. For instance, for ACGTG sequence with a window size of 4 and a stride of 1, the k-mers generated are ACGT and CGTA. Alternatively, in non-overlapping k-mers generation, window and stride size must be equal in size. For same sequence used in the overlapping case, this non-overlapping approach generates only one k-mer, such as ACGT. The length of the k-mer is determined by the window size. Researchers often create pre-trained embeddings with different k-mer sizes and then select the size which yields best performance in downstream tasks. Once k-mers are generated, these k-mers sequences are passed to traditional word embedding models (Word2vec, FastText, GloVe) to generate representation.

A high level overview of Figure 5 indicates that various studies have explored the potential of the Word2Vec embedding method in combination with 13 different machine and deep learning algorithms as well as 2 statistical algorithms. The era of word embedding approaches begin in 2013 with introduction of Word2Vec (221). Word2vec has two different embeddings generation paradigm: (1) SkipGram and (2) Continuous Bag of Words (CBoW). SkipGram learns representations of k-mers by predicting surrounding k-mers for every k-mers of corpus. The number of surrounding k-mers is a hyper-parameter that can be adjusted according to available data. Contrarily, CBoW model learns k-mers representations by predicting single k-mer based on the context of its surrounding k-mers. Similar to SkipGram model, here context of surrounding k-mers is a hyper-parameter. Word2vec architecture is comprised of input layer, hidden layer, and an output layer. At input layer, a random d-dimensional vector is initialized for each k-mer, while the hidden layer extracts relationships between k-mers. These relationships are further passed to output layer, which predicts probabilities of output k-mers based on the context of input k-mers. The predicted probabilities are passed to loss function which computes loss value. To facilitate readers, Equation 10 embodies mathematical expressions for computing loss values of both variants.

In above expression, N refer to number of k-mers, wi indicates target k-mers, wj is one of k-mers within contextual window, and Ws refers to set of k-mers in contextual windows of k-mers wi. After computing loss, weights are updated during back propagation which eventually helps in generating similar vectors for similar k-mers and distinct vectors for dissimilar k-mers.

Pennington et al. (222) proposed another k-mers embedding approach named Global Vectors (GloVe) which generates k-mers vectors by capturing both global and local contextual information of k-mer within corpora. It can be seen in Figure 5, in the context of DNA sequence analysis, the potential of Glove k-mers embedding method is explored with two distinct deep learning methods. Primarily, this embedding generation method computes local and global contextual information by incorporating occurrence frequencies of k-mer pairs into an objective function shown in Equation 7.

In above expression, wi and wj are k-mers within a pair, bi and bj are corresponding biases, and f(Cij) is a weighted function to normalize co-occurrence matrix values and eradicate biases and their impact of noise on k-mers embeddings.

Figure 5 shows that Word2Vec is the most commonly explored word embedding method, followed by FastText. Mikolov et al. (223) proposed FastText approach by extending the working paradigm of Word2Vec model. Primarily, this approach handles out-of-vocabulary (OOV) k-mers by discretizing k-mers into sub k-mers. After generating sub k-mers. it takes average of sub k-mers vectors to generate k-mers vectors and passed them to word2vec model. During back propagation, it updates vectors of both k-mers and sub k-mers. Through this strategy, vectors are generated for both k-mers and sub k-mers.

Furthermore, in NLP domain, with an aim to generate more comprehensive vectors of k-mers by capturing k-mers informative patterns from textual corpora, researchers have proposed different graph-based methods. These approaches include DeepWalk (224), Node2Vec (225), Graph2Vec (226), SDNE (227), SocDim (228), GraRep (229), Laplacian Eigenmaps (230), Locally Linear Embedding (231), and OPA2Vec (232). Figure 5 highlights that within the context of DNA sequence analysis, the potential of the 7 graph-based methods is less explored compared to the two foundational word embedding methods, Word2vec and FastText. In addition, among the graph-based methods, Node2Vec (225) has been investigated more extensively than DeepWalk (224), Graph2Vec (226), SDNE (227), SocDim (228), GraRep (229), Laplacian Eigenmaps (230), Locally Linear Embedding (231), and OPA2Vec (232). Similar to non-graph-based methods, graph-based methods segregate sequences into k-mers and generate k-mers pairs by sliding a 2 size window over k-mers sequences. By using k-mers paris, a graph is formed where nodes represent k-mers, and edges represent relationships between the k-mers. For example, to generate a graph from the DNA sequence ACTGCA with k = 3, first, overlapping k-mers (ACT, CTG, TGC, GCA) are generated. By sliding a window of size 2 over these k-mers sequence, k-mers pairs [(ACT, CTG), (CTG, TGC), and (TGC, GCA)] are created. These pairs form edges of graph, with k-mers serving as nodes. Perrozi et al. (224) proposed DeepWalk approach that utilizes graphical space to generate new sequences by capturing extensive relationships between k-mers. After generating new sequences, it makes use of Word2Vec model for generation of k-mers vectors. In contrast, Grover et al. (225) proposed Node2Vec approach that utilizes a distinct strategy for generation of new sequences. Primarily, Node2Vec employs second order random walk sampling strategy which reaps the benefits of breath first search (BFS) and depth first search (DFS) algorithms. This strategy computes probability of visiting next node depending on the previously visited nodes rather than just randomly selecting one of neighboring nodes. Naeayanan et al., (226) introduced another embedding generation approach namely Graph2Vec. It extracts root node, its sub-graph, and degree of intended sub-graph to generate a sorted list of nodes which is then passed to SkipGram with negative sampling (SGNS) model.

Matrix factorization embedding approaches extend graph-based embedding approaches by using adjacency matrix rather than generating new sequences directly from graph. Adjacency matrix encodes the relationships between nodes within the graph which is then decomposed using matrix factorization methods namely SVD and NMF. These approaches also decompose adjacency matrix of graph into lower-dimensional matrices which represents node embeddings. These embeddings extract nodes latent features and relationships between them. Mainly, matrix factorization methods aim to minimize reconstruction error between original adjacency matrix and reconstructed matrix from node embeddings. These methods include Laplacian Eigenmaps (230), Locally Linear Embedding (231), SDNE (227), SocDim (228), GraRep (229), and OPA2Vec (232). A closer view of Figure 5 indicates that 6 matrix factorization embedding approaches method are least explored as compared to foundational word embedding methods (Word2vec, FastText) and graph-based methods.

Laplacian Eigenmaps (230) approach derives degree matrix from adjacency matrix and computes graph Laplacian matrix by computing the difference between degree matrix and adjacency matrix. Next, it computes eigen values and constructs eigenvectors corresponding to smallest non-zero eigenvalues which results in generating lower-dimensional k-mer embeddings and preserving local k-mers relationships. Another matrix representation approach graph representations (GraRep) (229) make use of adjacency (Ai, j) and degree (Di, j) matrices driven from nodes and edges of graph. Equation 8 depicts mathematical expression for computing proximity matrix from Ai, j and Di, j matrices.

where V represents total number of nodes in graphs. Calculated proximity matrix is further passed to singular value decomposition approach for generating k-mer embeddings. Moreover, this approach focuses on extracting similarities between nodes by using k-step information of neighbors where levels of neighbors can be represented through k-steps. Similar to GraRep approach, SocDim (228) generates k-mer representations by incorporating social dimensions, namely, attributes and network structures. Specifically in SocDim, adjacency and degree matrices are used to compute modularity matrix by using Equation 9.

where m represents edges, and D represents degree matrix. Similar to GraRep, modularity matrix is passed to SVD for k-mers embeddings generation.

Moreover, structural deep network embedding (SDNE) (227) leverages deep auto-encoders to generate k-mer embeddings by determining first and second order proximities to ensure connected k-mers have similar embeddings. SDNE model architecture is trained to optimize combined loss function that incorporates both proximities and finally generates low-dimensional representations by capturing non-linear relationships between nodes and encoding structural information into embeddings. Afterward, structural embedding aims to address limitations of k-mer embeddings approaches in capturing structural and semantic information of nodes and edges in heterogeneous networks. Among structural embedding approaches, Opa2Vec (232) makes use of individual entities containing structured knowledge or characterized classes axioms and unstructured information or metadata, i.e., textual annotations and passes them semantic reasoner tool (Elk/HermiT) for generating ontology sequence which is then passed to Word2Vec model for generating representations. Locally linear embedding (LLE) (231) method identifies neighboring k-mers for each k-mer in the sequence and determines weights by employing graph Laplacian concept which linearly reconstructs each k-mer from its neighbors. Afterward, it computes sum of edges between close k-mers by using heat-kernel method which ensures weights of connected k-mers as 1 and unconnected k-mers as 0, ultimately maintaining the reconstruction relationship. These weights extract both semantic and syntactic information and maintain the reconstruction relationship. By optimizing reconstruction error and computing eigenvectors, LLE generates embeddings for each k-mer in the sequence. These embeddings represent the k-mers in a reduced-dimensional space, where similar k-mers in context are closer together.

Specifically, for DNA sequence analysis tasks, word embeddings methods are being utilized to generate pre-trained embeddings in 2 different ways: In one way, sequences are segregated into k-mers and embeddings of k-mers are generated. In second way, embeddings are generated for whole DNA sequence. Moreover, most of the DNA sequence analysis predictors follow first way to generate embeddings (, 39, 45, 4749, 58, 59, 6466, 76, 78, 79, 82, 92, 105, 106, 113, 121, 122, 136138, 151, 153, 163, 165, 171, 173, 233235), but second way is utilized by only few tasks including gene-disease association prediction (110), pseudogene function prediction (111), promoter identification (236), essential gene prediction (107, 108), gene network reconstruction (128), and gene expression prediction (103). In this section, we have defined methods from first way perspective. A comprehensive detail about second way is available in following articles (103, 107, 108, 110, 111, 128, 236).

In a nutshell, word embedding approaches have significantly propelled 44 distinct DNA sequence analysis tasks, enriching the research community with the development of robust and precise models. Notably, conventional word embedding techniques such as Word2Vec, GloVe, and FastText excel in capturing k-mers context and sub k-mers information effectively. In contrast, innovative techniques such as Graph2Vec, Node2Vec, DeepWalk, and GraRep harness graph-based methodologies to enhance embeddings based on connectivity and proximities. In addition, SocDim and OPA2Vec offer distinctive perspectives by integrating social and ontological elements, while SDNE combines local and global structural insights through deep autoencoders. Locally linear embedding (LLE) and Laplacian eigenmaps are dedicated to preserving local geometric properties. Ultimately, each approach makes a distinctive contribution to driving significant progress in DNA sequence analysis.

7.2 DNA sequence representation learning using language models

In the evolving landscape of natural language processing (NLP), the inception of the Transformer model has announced a new era of advancements, setting the precedent for subsequent developments in language models (237, 238). The Transformer and distinct language models, including BERT, GPT-3, and ELECTRA, have significantly contributed to pushing the boundaries of what machines can understand and generate in terms of human language (237, 238). The importance of these models lies not only in their ability to comprehend and produce text but also in their application across different domains including genomics and proteomics sequence analysis (239). These models have found multifarious applications in genomics and proteomics sequence analysis tasks by generating highly effective representations of biological sequences (239). To facilitate DNA sequence analysis researchers, here we briefly delve into the key features, advantages, and disadvantages of commonly used eight modern sophisticated language models, namely, Transformer (102), Transformer-XL (50), XLNet (156), ULMFIT, BERT (156), ALBERT (156), ELECTRA (156), and GPT-3 (240). Table 3 presents 8 distinct language models and their variants, categorized into 4 different groups based on their architectures. These architectures include trivial LSTM-based language model, encoder-decoder architecture, encoder-only architecture, and decoder-only architecture. Table 3 also provides information about language model architecture and outlines number of layers as well as count of encoders or decoders and their respective layers.

Table 3

Architecture typeLanguage model, release yearLanguage model variantsNumber of layers in encodersNumber of layers in decoders
Trivial LSTM based language modelULMFiT, (243), 2018AWD-LSTM language model1_
Encoder-decoderTransformer_XL (336), 2019Transformer_XL Large (WikiText-103)2424
24L Transformer_XL (text8)2424
12L Transformer (enwik8)1212
18L Transformer (enwik8)1818
24L Transformer (enwik8)2424
Transformer_XL base (Billion Word)1212
Transformer_XL large (Billion Word)2424
Encoder-decoderTransformer, (241), 2017Base66
Big66
ELECTRA, (246), 2020Small12_
Base12_
Large24_
ALBERT, (245), 2020Base12_
Large24_
xLarge24_
xxLarge12_
BERT, (244), 2019Base12_
Large24_
XL-Net, (242), 2019Base12_
Large24_
Decoder-onlyGPT, 2018GPT-1 (337)_12
GPT-2 small (338)_12
GPT-2 medium (338)_24
GPT-2 Large (338)_36
GPT-3 (247)_96
GPT-4 (339)_120

Summary of 8 contemporary language models used in DNA sequence analysis.

The Transformer model, introduced in 2017 by Vaswani et al., (241) marks a significant departure from previous models that relied on recurrent or convolutional neural networks for processing sequential data. This model utilizes a unique architecture that focuses on attention mechanisms which allows to handle long-range dependencies and understand the context and semantics of sequences more effectively (102, 241). Key innovations of the Transformer include positional encoding and self-attention mechanisms (102, 241). Positional encoding assigns a unique number to each individual k-mer or group of k-mers and helps in grasping k-mers order and sequence context. The self-attention mechanism allows the model to weigh the importance of each k-mer in relation to others, enhancing its ability to process and predict scientific language patterns (102, 241). The main advantage of the Transformer is its efficiency in training and inference due to parallel processing of sequences (102, 241). However, it requires substantial computational resources, which can be a limiting factor in resource-constrained environments. Despite this, its flexibility and scalability in handling diverse genomics tasks make it a preferred choice in many advanced AI applications (102).

Transformer-XL extends the Transformer architecture to address the limitation of fixed-length context by incorporating mechanisms that capture long-range dependencies more effectively (50). This model enhances the ability to maintain context over longer sequences than standard Transformer models, which significantly improves performance in various genomics and proteomics sequence analysis tasks (50). The core innovations of Transformer-XL include the introduction of a segment-level recurrence mechanism and a novel relative positional encoding (50). These features allow the model to reuse past information and thereby extend the context window across different segments. This design enables Transformer-XL to handle longer biological sequences efficiently and provides a substantial improvement over traditional models where each segment is processed in isolation (50). One of the main advantages of Transformer-XL is its capability to learn dependencies that are significantly longer than those captured by traditional models, leading to improvements in both short and long sequence analysis tasks (50). However, the model demands more memory due to its recurrence mechanism and larger context handling, which could be a limitation in resource-constrained environments.

XLNet extends the Transformer-XL model using an autoregressive method (242). This approach allows XLNet to learn bidirectional contexts by maximizing the expected likelihood over all permutations of the input sequence order which significantly enhances its scientific language understanding capabilities (242). Primary innovation of XLNet is its permutation language modeling (PLM), which enables the model to predict the likelihood of a sequence by considering all permutations of the k-mers within it (156, 242). This method allows XLNet to capture a comprehensive bidirectional context, unlike traditional autoregressive models only consider a single direction. In addition, XLNet incorporates a two-stream self-attention mechanism which enhances its ability to manage the context more effectively during the prediction process (156, 242). One of the main advantages of XLNet is its robustness in modeling bidirectional contexts, which significantly outperforms previous models such as BERT in numerous genomics sequence analysis tasks (156, 242). However, the complexity of its training process, which involves permutation of input sequences and a two-stream attention mechanism, may pose challenges in terms of computational resources and time (156, 242).

Universal Language Model Fine-tuning (ULMFiT) has revolutionized natural language processing by introducing effective transfer learning techniques for various NLP tasks. It is developed by Jeremy Howard and Sebastian Ruder in 2018 (243), and it typically leverages a pre-trained language model which is fine-tuned on specific DNA sequence analysis tasks having minimal sequences (243). ULMFiT utilizes Average Stochastic Gradient Descent - Long Short-Term Memory (AWD-LSTM) architecture to learn the distribution and contextual relationships of k-mers in DNA sequences (57). It employs self-supervised learning that predicts the next k-mer based on the previous known k-mers and enables the model to capture the semantics and discriminative potential of the sequences (57). The core innovation of ULMFiT lies in its ability to fine-tune pre-trained language models using techniques such as discriminative fine-tuning and the slanted triangular learning rates policy. Discriminative fine-tuning considers that different layers of neural network capture different kind of information; hence, it tunes every layer with distinct learning rates (243), whereas slanted triangular learning rate describes a unique learning rate scheduler that initially increases the learning rate and afterward drops it in a linear fashion (243). The short increase stage enables the model to quickly converge to a parameter space suitable for the task, while the extended decay period allows for more effective fine-tuning (243). By adjusting the learning rate for different layers, it prevents catastrophic forgetting and stabilizes the training process across various tasks (243). ULMFiT incorporates dropout techniques to regularize learnable parameters and prevent overfitting which ensures model's generalization ability (57). Another advantage of ULMFiT is its ability to achieve high performance with significantly less data compared to traditional models. However, the complexity of fine-tuning and the need for careful calibration of learning rates can be challenging, requiring a nuanced understanding of model behavior across different layers (57).

Bidirectional Encoder Representations from Transformers (BERT) is developed by Google in 2018 (244). It is pretrained on a large corpus of text data, such as Wikipedia and books (244). It has revolutionized NLP tasks by employing a transformer-based architecture that enables the model to consider the context of k-mers from both directions simultaneously, rather than a single direction at a time (244). BERT is distinctive for its deep bidirectional nature, achieved through the application of the transformer model, specifically using mechanisms such as Masked Language Modeling (MLM) and Next Sentence Prediction (NSP) (244). This approach allows BERT to understand the context of a k-mer based on all other k-mers in a sequence, rather than just those preceding it. Specifically, it learns to capture the semantics and contextual information of the input text exceptionally well through self-supervised learning tasks such as MLP and NSP (244). In the case of DNA sequence analysis, BERT is used to transform DNA sequences into statistical feature space and then fine-tuned on specific downstream tasks, such as enhancer identification and strength prediction (156). BERT captures the semantics of DNA sequences by dynamically learning their representations through a multihead self-attention mechanism. BERT leverages transfer learning by pre-training on a large corpus and then fine-tuning on specific DNA sequence analysis task, allowing it to adapt to different applications (156). BERT uses MLM and NSP tasks during pre-training to learn the contextual relationships between k-mers in DNA sequences (156).

The primary advantages of BERT include its high accuracy and efficiency across various DNA sequence analysis tasks, due to its robust handling of context and bidirectional training (156). BERT captures both discriminative and semantical relationships of k-mers, making it effective in characterizing DNA sequences (156). BERT-based models have shown improved performance compared to traditional approaches in different DNA sequence analysis tasks such as enhancer identification and strength prediction (156). In addition, BERT can be adapted to specific application scenarios by pre-training on domain-specific custom corpora (156). BERT is a large model that requires significant computational resources for training and inference on extensive datasets. BERT performs best when trained on large and diverse datasets, which may not always be available for specific DNA sequence analysis tasks. In addition, while BERT provides state-of-the-art results in many scenarios, it requires fine-tuning for specific tasks, which can be resource-intensive. BERT performance can degrade with longer texts and the complex architecture of BERT makes it challenging to interpret the learned representations and understand the underlying biological mechanisms (156).

ALBERT, introduced by Google researchers, is a streamlined version of BERT designed to provide state-of-the-art results in NLP with significantly fewer parameters (245). This model enhances the efficiency and scalability of BERT by incorporating innovative techniques such as factorized embedding parameterization, cross-layer parameter sharing, and sentence order prediction (156, 245). Factorized embedding parameterization technique reduces the size of the embedding matrix by separating the vocabulary and hidden layer sizes, which decreases the number of parameters significantly (156, 245). In cross-layer parameter sharing, parameters are shared across all layers of the model, reducing the total parameter count and improving training efficiency (156, 245). It replaces the next sequence prediction with sequence order prediction to enhance the model's ability to understand sequence coherence without requiring task prediction, making it more effective for downstream tasks (156, 245). The primary advantage of ALBERT is its reduced parameter size, which allows for faster training times and less memory usage compared to BERT, without a significant loss in performance. However, the extensive parameter sharing might lead to a slight decrease in model flexibility, potentially affecting task-specific fine-tuning (156, 245). Efficiently Learning an Encoder that Classifies Token Replacements Accurately (ELECTRA) (156, 246) has introduced a novel pre-training method for language models. ELECTRA operates on a replaced token detection (RTD) mechanism, where it differs from traditional masked language models such as BERT (156, 246). Instead of masking k-mers, ELECTRA corrupts the input by replacing tokens or k-mers with outputs from a generator model, challenging the discriminator to identify changes (156, 246). This approach allows the model to learn from the entire input sequence, enhancing training efficiency. The primary advantage of ELECTRA lies in its efficiency, requiring less computational power and time to reach or exceed the benchmarks set by larger models (156, 246). However, the complexity of its dual-model architecture, involving both a generator and a discriminator, might pose challenges in training stability and hyperparameter tuning.8

GPT-3 is one of the most advanced AI language models developed by OpenAI (247). It is recognized for its ability to generate text that closely mimics human writing, making it a pivotal development in natural language processing. GPT-3 builds upon the transformer architecture, which utilizes self-attention mechanisms to process input data (247). Unlike GPT-2, which had 1.5 billion parameters, GPT-3 boasts a staggering 175 billion parameters. This exponential increase in parameters enhances its ability to generate coherent and contextually relevant text (247). GPT-3 differs from models such as BERT and XLNet by maintaining an autoregressive nature. From scientific perspective, this implies that it predicts the next k-mer in a sequence based on the previous k-mers, while BERT uses bidirectional context (240, 247). One of the innovative aspects of GPT-3 is its use of alternating dense and locally banded sparse attention patterns. Dense attention considers all input k-mers simultaneously, while sparse attention focuses on a subset of k-mers, making the model more efficient and scalable. This combination enables GPT-3 to handle long-range dependencies and maintain computational efficiency (240, 247). One of GPT-3's standout capabilities is its performance in few-shot settings. Unlike fine-tuned models that require large amounts of task-specific data, GPT-3 can perform well on new tasks with minimal sequences. This flexibility is a significant advantage over models such as BERT, which typically require extensive fine-tuning for each specific task. GPT-3 demonstrates strong performance across various tasks, often matching or exceeding that of fine-tuned models. This capability makes it a versatile tool for a wide range of applications (240, 247).

For instance, in context of cell biology, scientific researchers have used GPT-3 to learn gene and cell embeddings effectively (240). Scientific researchers have utilized the text summaries of genes from the NCBI database, which contain curated information about gene functionalities and properties. The gene text summaries are passed through the GPT-3 language model, which generates gene embeddings that capture the underlying biology described in the gene summaries (240). The gene embeddings are averaged, weighted by the expression levels of each gene in the cell. These averaged embeddings are then normalized to a unit l2 norm to generate single-cell embeddings (240). In another strategy, each cell is represented by a natural language sentence constructed based on the ranked gene expressions. The gene names are ordered by descending normalized expression levels, and this sentence representation is passed through the GPT-3 model to obtain the cell embeddings (240). Extrinsic performance analysis of GPT-3 embeddings on tasks such as classifying gene properties or cell types has shown supreme effectiveness (240). While GPT-3's capabilities are groundbreaking, it faces challenges such as potential biases in training data and high computational demands. Moreover, its “black box” nature makes it difficult to discern how decisions are made, posing ethical and operational concerns. GPT-3's massive size requires significant computational resources for both training and inference. This makes it less accessible for smaller organizations or researchers without high-end hardware.

7.3 Machine and deep learning predictors

Machine and deep learning algorithms need statistical vectors to extract useful patterns for specific sequence analysis task. A comprehensive literature review of 127 studies reveals that 12 word embedding and 8 large language models have been used to generate statistical vectors of raw sequences to feed 28 different algorithm available within predictive pipelines of 44 DNA sequence analysis tasks. Based on working paradigms, these algorithms are categorized into 3 different categories, namely, statistical algorithms, machine learning algorithms, and deep learning algorithms. From 28 algorithms, 3 algorithms, namely, conditional random fields (CRF) (248), k-means clustering algorithm (), and cosine similarity algorithm (173), belong to statistical algorithms. Machine learning algorithms involves 8 algorithms, namely, support vector machine (SVM) (352), Naive Bayes (NB) (95), multilayer perceptron (MLP) (77), predictive clustering tree (PCT) (128), random forest (RF) (103), deep forest (DF) (61), XGBoost (352), and CatBoost (143). Furthermore, deep learning algorithms include convolutional neural network (CNN) (91), graph neural network (GNN) (104), temporal convolutional network (TCN) (235), graph convolutional network (GCN) (55), graph attention network (GAT) (108), long short-term memory (LSTM) (97), bidirectional long short-term memory (BiLSTM) (58), and bidirectional gated recurrent unit (BiGRU) (82). Similarly, 8 algorithms, namely, ELECTRA (89), ALBERT (249), Transformer-XL (50), XL-Net (156), Transformer (250), ULMFit (146), GPT-3 (114), and BERT (251), belong to language modeling algorithms, five algorithms, namely, LSTM + CNN, CNN + BiLSTM, CNN + BiLSTM + BiGRU, RF + CNN, and CNN + BiGRU, belong to hybrid algorithms, whereas 1 meta-predictor reaps benefits of both machine and deep learning algorithms, namely, KNN, RF, SVM, MLP, and CNN.

Statistical algorithms provide a framework for understanding DNA sequence distribution and characteristics. They offer valuable advantages in terms of interpretability by facilitating researchers to assess statistical significance of genomic features. Among three statistical algorithms, Conditional Random Fields (CRF) (248) calculate the conditional probability of class labels of sequences by using neighboring k-mers. By capturing dependencies between adjacent labels, CRF allows for more accurate predictions of sequence features by taking into account both local sequence context and broader genomic patterns. K-means clustering algorithm () groups sequences into k distinct clusters based on similarity. It starts by initializing k centroids and assigns each sequence to the nearest cluster by calculating Euclidean distance between sequence and centroid. These centroids are updated iteratively by averaging the sequences in each cluster until they stabilize. Cosine similarity can be advantageous in DNA sequence analysis for tasks such as similarity comparison and clustering, where measuring the similarity between sequences is essential (252). Cosine similarity can handle high-dimensional data efficiently and is suitable for tasks requiring similarity-based analysis (252). These models also have limitations, including potential difficulties in managing complex and high-dimensional data. In addition, they rely on strong assumptions about underlying data distribution, which may not always align with real-world DNA sequence analysis scenarios. Despite these challenges, statistical models remain indispensable tools in DNA sequence analysis and provide valuable insights.

From 8 different machine learning algorithms, support vector machine (SVM) operates by finding the optimal hyperplane that best separates data points into different classes. SVMs are known for their ability to handle high-dimensional data and work well in cases where the data are not linearly separable, as they can use kernel functions to transform the data into higher dimensions where separation is possible (62). However, SVMs can have limitations in terms of training time, especially with large datasets, as they need to solve a complex optimization problem to find the best hyperplane that separates the classes. Naive Bayes (NB) is a probabilistic algorithm based on Bayes' theorem with the assumption of independence between features. NB is efficient, simple to implement, and works well with high-dimensional data, making it suitable for tasks where feature independence assumptions hold (253). However, it may not always hold true in practice, especially in complex biological datasets where features are correlated.

In addition to SVM, tree-based algorithms are fundamentally built upon decision tree algorithm. Decision tree algorithm uses independent variables to construct a tree-like structure, where data are split at decision nodes into branches connected to leaf nodes to make predictions. This foundational algorithm is extended into more advanced algorithms, namely, Random Forest (RF), Deep Forest (DF), XgBoost, CatBoost, and Predictive Clustering Tree (PCT) (128). All of these advanced algorithms enhance basic decision tree by incorporating techniques such as ensembling, and boosting for improved accuracy and generalization. Random Forest (RF) algorithm is an ensemble learning method that constructs a multitude of decision trees during training and outputs the mode of the classes as the prediction. RF is known for its robustness to overfitting, feature importance estimation, and ability to handle high-dimensional data with ease (254). However, RF may not perform as well when dealing with imbalanced datasets or when there are many irrelevant features present in the data. Deep forest (DF) algorithm is another ensemble learning method that utilizes a cascade structure of multiple random forests to make predictions. DF can be advantageous in DNA sequence analysis for tasks such as clustering and species classification based on DNA barcodes (255). DFs are capable of learning hierarchical representations of data and can capture complex patterns in high-dimensional spaces effectively (255). Nonetheless, the main drawback of DF lies in its computational complexity and the need for substantial computational resources, which can limit its practicality in large-scale DNA sequence analysis projects. XGBoost combines multiple weak learners to create a strong predictive model. XGBoost can handle large datasets with high dimensionality and is known for its efficiency in boosting the performance of weak learners (256). However, XgBoost may require fine-tuning of hyperparameters to achieve optimal performance, and it could be sensitive to noisy data. CatBoost is another ensemble learning method designed to handle categorical features efficiently. CatBoost can automatically handle categorical features and is known for its robustness to overfitting and efficiency in training models with categorical data (256). Nevertheless, CatBoost's training time might be longer compared to other algorithms, especially when dealing with large genetic datasets.

Predictive clustering tree (PCT) (128) is a versatile predictor that integrates elements of both clustering and supervised learning. Unlike traditional decision trees, random forests, or support vector machines, PCTs are designed to handle hierarchical multi-label classification tasks, making them particularly effective for complex, high-dimensional data (128). PCTs operate by viewing a decision tree as a hierarchy of clusters. The root node represents a single cluster containing all training examples, which is recursively partitioned into smaller clusters as one moves down the tree. This approach allows PCTs to simultaneously perform clustering and classification, leveraging the hierarchical structure to predict multiple labels for each instance (128). One of the key strengths of PCTs is their ability to manage complex data with multiple interrelated labels. They can identify relevant features across different levels of the hierarchy, providing interpretable results that are valuable for domain experts. In addition, PCTs are capable of handling large datasets efficiently, making them suitable for various real-world applications (128). Despite their strengths, they can be computationally intensive, especially for large and deep hierarchies, and may require careful parameter tuning to avoid overfitting. In addition, while PCTs offer interpretability, the complexity of the hierarchical structure can sometimes make the results harder to interpret compared to simpler models (128). Apart from this, researchers have also designed customized meta-predictors which utilize the powers of five or more than five distinct algorithms, namely, kNN, RF, SVM, MLP, and CNN (105).

Multilayer perceptron (MLP) is composed of multiple layers of nodes that can learn complex patterns in data. MLPs are powerful algorithms for feature extraction and predictive modeling in DNA sequence analysis, capable of capturing intricate relationships in the data (252). MLPs excel in tasks requiring non-linear decision boundaries and can handle large amounts of data effectively. However, training MLPs can be computationally expensive, especially with large datasets, and they are prone to overfitting if not properly regularized.

Among all categories, deep learning algorithms are most extensively used for efficient DNA sequence analysis. A total of eight deep learning algorithms are most commonly used by scientific community for DNA sequence analysis. Convolutional neural network (CNN) is a deep learning algorithm designed to process structured grid-like data, such as images. In DNA sequence analysis, CNNs can be applied to DNA sequence analysis tasks to capture spatial dependencies in data. They are effective for tasks that require feature hierarchies and translation invariance (257). However, CNNs may struggle with capturing long-range dependencies in sequences, which can be crucial in DNA analysis where distant k-mers may interact. Graph neural network (GNN) is a type of neural network designed to operate on graph-structured data. GNNs are suitable for tasks involving relational data, such as molecular structures, making them applicable to DNA sequence analysis for tasks such as clustering (258). GNNs can effectively capture dependencies between nodes in a graph and are capable of learning representations that incorporate both local and global information (258). However, GNNs may encounter challenges in efficiently scaling to large graphs, and interpreting the learned representations in GNNs can be complex, limiting their interpretability. Temporal convolutional network (TCN) is a type of neural network designed to process sequential data efficiently. TCNs are suitable for tasks involving temporal dependencies, making them applicable to DNA sequence analysis for tasks like predicting DNA binding sites for transcription factors. TCNs can capture long-range dependencies in sequential data and are known for their parallel processing capabilities, enabling faster training times (259). However, TCNs may struggle with modeling complex temporal dynamics compared to recurrent models such as LSTMs. Graph convolutional network (GCN) is a type of neural network designed to operate on graph-structured data. GCNs can leverage graph structures to learn representations of nodes and edges, enabling tasks such as node classification and link prediction in DNA sequences (256). However, GCNs may require meticulous graph construction and preprocessing, and they can be computationally intensive, especially for large graphs, which can hinder their scalability.

Graph attention network (GAT) is a type of neural network that incorporates attention mechanisms to learn the importance of neighboring nodes in a graph. GATs are suitable for tasks involving relational data, as they can adaptively weigh the contributions of neighboring nodes, enabling more flexible and accurate learning on graph-structured data (260). However, GATs may be sensitive to noisy or sparse graphs, and designing optimal attention mechanisms can be challenging, impacting their performance in certain scenarios. Long short-term memory (LSTM) is a type of recurrent neural network designed to capture long-term dependencies in sequential data. LSTMs are effective in DNA sequence analysis for tasks such as hypersensitive DNA sequence classification (261). LSTMs can retain information over long sequences and are suitable for tasks requiring memory of past events, making them ideal for tasks such as classification of DNA sequences (261). However, LSTMs may encounter vanishing or exploding gradient problems during training, which can affect their ability to capture long-term dependencies accurately. Bidirectional long short-term memory (BiLSTM) is an extension of LSTM that processes sequences in both forward and backward directions. BiLSTMs are advantageous in DNA sequence analysis for tasks where contextual information from both past and future is essential (262). BiLSTMs can capture dependencies in both directions and are effective in tasks requiring bidirectional context understanding (262). However, BiLSTMs may be computationally intensive due to processing sequences in two directions, which can impact their training and inference speed. Bidirectional gated recurrent unit (BiGRU) is another type of recurrent neural network that combines the advantages of bidirectionality and gating mechanisms. BiGRUs can capture bidirectional dependencies efficiently and are known for their simpler architecture compared to LSTMs, making them computationally more efficient (256). However, BiGRUs may struggle with capturing very long-term dependencies compared to LSTMs, which can limit their effectiveness in tasks requiring extensive memory retention.

For different DNA sequence analysis tasks, eight contemporary language models, namely, ELECTRA, ALBERT, Transformer-XL, XLnet, Transformer, ULMFit, GPT-3, and BERT have been used in two different settings. In first setting, the addition of classification layers to these language models adapts the general-purpose language models to specific classification tasks by learning to map the rich contextual embeddings to the desired output classes. In second setting, rich contextual embeddings of these 8 language models are passed to standalone machine learning algorithms, deep learning algorithms, and ensemble or hybrid algorithms for accurate classification of DNA sequences.

A total of five hybrid algorithms combine different types of models to leverage the strengths of each component. LSTM + CNN, CNN + BiLSTM, CNN + BiLSTM + BiGRU, RF + CNN (76), and CNN + BiGRU are some of the examples of hybrid algorithms that integrate deep learning and traditional machine learning techniques to enhance predictive performance (263) for different DNA sequence analysis tasks. These hybrid models aim to capitalize on the complementary advantages of different algorithms to achieve superior results in various tasks.

8 Uncovering evaluation measures for DNA sequence analysis predictive pipelines

AI-driven DNA sequence analysis predictive pipelines are evaluated using two different experimental settings: (1) k-fold cross-validation (48, 78) and (2) Train-test split (108, 110). In k-fold cross-validation, dataset is splitted into k folds, where k−1 folds are used for training and one fold is used for testing. In next iterations, from k-folds, another fold is reserved for testing whereas remaining k−1 folds are used for training. In this way, pipelines are trained and tested k times on whole data. This method offers more precise assessment of model generalization capability. Specifically for deep learning models (236), an additional set, known as validation set, is created from training set which typically uses 10% of training data. This validation set is used to optimize the model's hyperparameters. In train-test split experimental setting, data are divided into two distinct sets: (a) train set and (b) test set. Train set comprises majority of data (usually 70%–80%), while test set contains remaining 20%–30%. Similar to k-fold cross-validation, validation set is also created from train set for deep learning models.

Among 127 DNA sequence analysis studies, 67 studies have utilized 5-fold cross-validation-based experimental setting. Thirty five studies have used 10-fold cross-validation-based setting and 17 studies have used train test split-based setting. Eight studies have used both k-fold cross validation and train test split-based setting. Performance and effectiveness of trained predictive pipelines highly depends on ability to handles new and unseen data. To assess effectiveness and performance of predictive pipelines from different perspectives, various evaluation measures have been proposed. Based on task type, these measures are categorized into four classes: binary (92, 153)/multi-class classification (27, 173), multi-label classification (111, 112), regression (102, 103), and clustering (). Following subsections summarize details of all four types of evaluation measures.

8.1 Binary or multi-class classification evaluation criteria

Most commonly used evaluation measures in this category are accuracy (152, 153, 264), precision (137, 152), recall (152, 153), specificity (137, 153), F1 Score (137, 152), and MCC (137, 153). These measures are typically calculated using confusion matrix, which consists of four entities: true positives (TP), false positives (FP), true negatives (TN), and false negatives (FN) (265). Figure 6 makes use of aforementioned four entities to compute distinct evaluation measures.

Figure 6

It can be seen in Figure 6 that TP and TN indicate correct positive and negative predictions, while FP and FN signify incorrect positive and negative predictions. Equation 10 embodies mathematical expressions to compute aforementioned measures.

An in-depth assessment of existing DNA sequence analysis predictive pipelines reveals that most widely used evaluation measures for balanced datasets are F1-score, precision, accuracy, recall, specificity, and Matthews correlation coefficient (MCC). However, for imbalanced datasets, micro, macro, and weighted versions of these measures are used. To address class imbalance issue, weighted precision (266) considers both precision and relative weight of each class. Precision of a class is ratio of true positives to total number of positives for that class, while relative weight is proportion of samples of that class relative to total number of samples. Similarly, weighted recall (266) and weighted F1-score (27) are calculated by determining weights, recall, and F1-score for each class. Macro precision (146) calculates precision for each class independently and then averages these values. Macro recall (146) and macro F1-score (146) average recall and F1-score across all classes by considering each class equally regardless of size. In contrast, micro precision (146) calculates precision globally by considering all true positives and false positives across all classes together. Micro recall (146) and micro F1-score (146) aggregate TP, FP, and FN across all classes and provides a fair and balanced evaluation of predictor performance. Equation 11 provides mathematical expressions for computing these measures.

where , , and refer to true positives, false positive, and false negatives in class i, respectively. Pri, Ri, and F1−scorei signify precision, recall, and F1-score of class i. wi indicates relative weight of class i, and i refers to ith class among n classes.

8.2 Multi-label classification evaluation measures

Performance evaluation of multi-label classification predictive pipelines is more challenging compared to binary and multi-class classification predictive pipelines (267). In binary or multi-class classification, each sample is assigned to only one class at a time, so predicted class label will be either true or false. Contrarily, in multi-label classification, a sample belongs to two or more labels simultaneously and predictive pipelines predicts multiple labels (267). Among predicted labels, some labels can be correct, some labels can be incorrect, or all predicted labels can be correct or incorrect. This partial correctness introduces complexity. To address this problem, researchers have developed various evaluation measures, namely, accuracy (268), precision (268), recall (268), and hamming loss (269).

Equation 12 illustrates the mathematical expressions for these evaluation measures.

In these equations, N represents total number of samples, ni denotes ith sample out of N samples, Ai is actual class label, and Pi is the predicted label for ni sample. l represents sample length, j denotes class index, ∨ signifies logical OR operator, and ∧ represents logical AND operator. Similar to evaluation measures in binary or multi-class classification and includes micro precision, micro recall, micro F1-score, macro precision, macro recall, macro F1-score, weighted precision, weighted recall, and weighted F1-score. A rigorous analysis of existing literature on multi-label classification tasks in DNA sequence analysis reveals that most widely used evaluation metrics are accuracy, precision, recall, F1-score, MCC, sensitivity, and specificity.

8.3 Regression evaluation criteria

Regression tasks differ fundamentally from classification tasks where model predicts continuous numerical values rather than discrete class labels. Regression task-related predictive pipelines are evaluated using distinct evaluation measures including mean squared error (MSE) (270), mean absolute error (MAE) (271), mean bias error (MBE) (272), mean absolute percentage error (MAPE) (273), root mean square error (RMSE) (271), R2 (274), relative mean absolute error (rMAE) (275), relative mean square error (rMSE) (275), relative root mean square error (rRMSE) (275), and relative mean bias error (rMBE) (275).

MAE assesses predictor performance by measuring absolute difference between predicted and actual values (271). MSE quantifies deviation by averaging squared differences between actual and predicted values (270). Similarly, RMSE calculates standard deviation of prediction errors and demonstrates how tightly data points cluster around regression line (271). MBE assesses predictor performance in terms of under and overfitting by enumerating average difference between predicted and actual value (272). MAPE calculates percentage variation between predicted and actual values (273). Smaller the values of MAE, MBE, MSE, and MAPE, better will be predictor performance. Higher value of R2 score signifies promising predictor performance as it measures proportion of variance in predicted dependent variable explained by independent variable to determine strength of relationship.

MAE, MSE, RMSE, and MAPE measures compute average error value for N number of data points. Relative performance evaluation can improve quality of performance evaluation by reducing the noise from data. For relative performance evaluation, the percentage error of each metric is computed relative to the average of actual values (275). It facilitates in controlling factors that influence predictor performance by relatively calculating ratio of particular error with average of actual values (275). Since data continuously vary and produce varying predicted values at different time intervals, an overall percentage error is computed to obtain relative error of all data points. Equation 13 embodies mathematical expressions for aforementioned evaluation metrics.

In Equation 13, M denotes total number of samples, Ai represents actual value, and Pi is predicted value where i denotes the sample number and Ā is the average of total actual values.

8.4 Clustering evaluation measures

In contrast to first three categories explained, clustering tasks aim to group similar samples based on their features without predefined class labels. In this task, prime objective is to use clustering algorithms and identify inherent patterns or structures within data. In these tasks, clusters of data samples with similar features are created, and predictors assign new data points to appropriate clusters (276). A higher similarity to a cluster indicates that data sample belongs to that cluster (276). To assess clustering predictive pipeline performance, researchers have introduced various evaluation measures including accuracy (277), normalized mutual information (NMI) (277), silhouette score (SS) (278), dunn index (DI) (279), and Davies-Bouldin index (DBI) (280).

Accuracy (277) is the proportion of correctly predicted samples to total number of samples. NMI (277) quantifies quality of predictor by measuring mutual information between predicted clusters and actual clusters. Mutual information refers to computed joint probability between predicted clusters and actual clusters. Silhouette score (278) measures how similar data samples are within a cluster compared to other clusters. BDI (280) evaluates average similarity ratio of each cluster with its most similar cluster. DI (279) computes ratio of minimum inter-cluster distance to maximum intra-cluster distance. Equation 14 embodies mathematical expressions for these evaluations measures.

In Equation 14, yi refers to predicted cluster, ci and cj indicate ith and jth clusters among n clusters. Moreover, I(yi, ci) signifies mutual information, H(yi) and H(ci) show entropy of predicted and actual clusters. d(yi) is the average distance from yi to all points in other clusters, and a(yi) is the average distance of yi to all points in that clusters. d(ci, cj) represents inter-cluster distance between cluster i and cluster j, represents mean distance from cluster mean for all observations in cluster i, while denotes mean distance from cluster median for all observations in cluster j. An extensive analysis of existing literature reveals that most commonly used evaluation measures are accuracy and normalized mutual information.

9 Open-source DNA sequence analysis predictive models

The public availability of predictor source codes, pretrained language models, and word embeddings significantly benefits researchers by preventing the need to reinvent the wheel. These resources enable researchers to build on existing work, utilizing pre-trained models and complete predictive pipelines to develop new, enhanced applications. By integrating new strategies into these established pipelines, they can create more powerful predictors. In addition, open-source access to these codes allows for the reproduction of predictor performance, fostering transparency and reliability in research. To expedite the establishment of more precise, robust, reliable, and efficient AI models for DNA sequence analysis and ultimately accelerate advancements in genomics and bioinformatics research, this section provides a summary of open source predictive pipelines developed using two representation learning approaches: word embeddings and large language models for 44 different DNA sequence analysis tasks.

Our analysis reveals that, out of 39 existing word embedding based DNA sequence analysis studies, only 25 studies have made the source codes of their predictive pipelines publicly accessible. In addition, source code of only 38 studies is publicly available out of 67 existing DNA sequence analysis studies based on large language models. Tables 4, 5 offer details on open-source codes for DNA sequence analysis predictive pipelines based on word embeddings and large language models, respectively. They also provide a summary of the representation learning methods and machine/deep learning predictors employed, along with links to the corresponding source codes.

Table 4

ReferencesTaskEmbedding approachClassifierSource Code
Ratajczak et al. (340)Disease genes identificationNode2VecMLPhttps://github.com/fratajcz/speos
Pan et al. (177)Phage-host interactions predictionSDNE, Word2VecMLPhttps://github.com/NWUJiePan/Code
Chen et al. (341)Chromatin accessibility predictionGraph2VecNAhttps://github.com/pinellolab/simba
Han et al. (342)Nucleosome position detectionWord2VecCNN + BiLSTM + BiGRUhttps://github.com/lliqi-echo/Nucleosome-positioning-based-on
-DNA-sequence-word-vector-and-deep-learning
Liao et al. (58)Enhancer identificationWord2VecCNN + BiLSTMhttps://github.com/WamesM/iEnhancer-DCLA
Inayat et al. (59)Enhancer identificationFastTextMLPhttps://github.com/salman-khan-mrd/IEnhancer-DFH
Zhang et al. (76)Promoter identificationWord2VecRF + CNNhttps://github.com/HaoWuLab-Bioinformatics/iPro-WAEL
Le et al. (78)Promoter identificationFastTextCNNhttps://github.com/khanhlee/deepPromoter
Hong et al. (82)Enhancer-promoter interactions predictionWord2VecCNN + BiGRUhttps://github.com/hzy95/EPIVAN
Min et al. (49)Enhancer-promoter interactions predictionWord2VecCNN + BiGRU + Matching Heauristichttps://github.com/Xzenglab/EPI-DLMH
Dao et al. (45)YY1-mediated chromatin loops predictionWord2VecCNNhttp://lin-group.cn/server/DeepYY1
Tran et al. (153)Methylcytosine sites predictionFastTextXGBoosthttps://github.com/khucnam/5mC_Pred
Zulfiqar et al. (136)Methylcytosine sites predictionWord2VecCNNhttps://github.com/linDing-groups/Deep-4mCW2V
Fang et al. (137)Methylcytosine sites predictionWord2VecCNNhttps://github.com/mat310/W2VC
Khanal et al. (138)Methylcytosine sites predictionWord2VecCNNhttp://nsclbio.jbnu.ac.kr/tools/4mC-w2vec/
Huang et al. (151)Methyladenine sites predictionWord2VecBiLSTMhttp://39.100.246.211:5004/6mA_Pred/
Le et al. (105)Essential genes identificationFastTextEnsemble (kNN + RF + SVM + MLP + CNN)https://github.com/khanhlee/eDNN-EG
Zhang et al. (106)Essential genes identificationNode2VecMLPhttps://github.com/xzhang2016/DeepHE
Nunes et al. (110)Disease genes predictionOPA2VecRFhttps://github.com/liseda-lab/KGE_Predictions_GD
Fan et al. (111)Pseudogene function predictionNode2VecGCNhttps://github.com/yanzhanglab/Pseudo2GO
Yilmaz et al. (173)Mutation susceptibility analysisWord2VecCosine similarityhttps://github.com/alperyilmaz/dna2vec_snp
Arango et al. (113)Target gene classificationFastTextMLPhttps://bitbucket.org/gaarangoa/metamlp/src/master
Shi et al. (122)Gene taxonomy classificationLSH + FastTextMLPhttps://github.com/Lizhen0909/LSHVec
Li et al. (23)DNA-binding proteins binding sites identificationWord2Vec, FastText, GloVeCNN, DCNN, CNN + BiLSTM, DCNN + BiLSTMhttp://bliulab.net/BioSeq-BLM/
Do et al. (165)Recombination spots identificationFastTextSVMhttps://github.com/khanhlee/fastspot

Summary of open-source word embedding based models in existing studies.

Table 5

ReferencesTask nameLanguage modelClassifierPre-train/Self-trainSource code
Zhang et al. (31)Chromatin accessibility predictionTransformerTransformerSelf-trainhttps://github.com/ykzhang0126/semantic-CAP
Nguyen et al. (44)Species classification, chromatin accessibility predictionTransformerTransformerSelf-trainhttps://github.com/HazyResearch/hyena-dna
Luo et al. (96)Protein-DNA binding sites predictionBERTBERTSelf-trainhttps://github.com/lhy0322/TFBert
Liu et al. (53)Protein-DNA binding sites predictionBERTCNNPretrainedhttps://github.com/YAndrewL/clape
Gao et al. (36)Long-range chromatin interaction prediction, prediction of context-specific functional impact of genetic variantsTransformerTransformerSelf-trainhttps://github.com/ZjGaothu/EpiGePT
Ni et al. (83)Enhancer-promoter interaction predictionTransformerTransformerSelf-trainhttps://github.com/NWAFUniyu/EPI-Mind
Reddy et al. (102)Gene expression predictionTransformerCNNSelf-trainhttps://github.com/anikethjr/promoter-models
Fishman et al. (32)Gene expression predictionTransformerTransformerSelf-trainhttps://github.com/AIRI-Institute/GENA_LM
Osseni et al. (180)Tumor type predictionTransformerTransformerSelf-trainhttps://github.com/dizam92/multiomic-predictions
Le et al. (281)Methyladenine sites predictionBERTCNNSelf-trainhttps://github.com/khanhlee/bert-dna
Tsukiyama et al. (144)Methyladenine sites predictionBERTCNN + BiLSTMSelf-trainhttps://github.com/kuratahiroyuki/BERT-6mA.git
Wang et al. (158)Methylation sites predictionTransformerTransformerSelf-trainhttps://github.com/sb111169/tf-5mc
Zhou et al. (250)Methylation sites predictionTransformerTransformerSelf-trainhttps://github.com/LieberInstitute/INTERACT
Jeong et al. (159)Methylation sites predictionBERTBERTSelf-trainhttps://github.com/CompEpigen/methyl-seq_simulation
Huson et al. (248)Methylation sites predictionBERTCRFSelf-trainhttps://github.com/husonlab/MR-DNA
Yu et al. (154)Methylation sites predictionBERTBERTPretrainedhttps://github.com/YUYING07/iDNA_ABT
Zhou et al. (155)Methylation sites predictionBERTBERTPretrainedhttps://github.com/wrab12/StableDNAm
Jin et al. (157)Methylation sites predictionBERTFGMPretrainedhttps://github.com/FakeEnd/iDNA_ABF
Stanojevic et al. (152)Methylcytosine sites predictionTransformerTransformerPretrainedhttps://github.com/lbcb-sci/rockfish
Wang et al. (282)Methylcytosine site predictionBERTBERTSelf-trainhttps://zenodo.org/records/10143217
Yang et al. (143)Methylcytosine sites predictionBERTCatBoostPretrainedhttps://github.com/abcair/4mCBERT
Yang et al. (249)Conserved non-coding element classificationTransformer, ALBERTTransformer, ALBERTSelf-trainhttps://github.com/melobio/LOGO
Fazeel et al. (29)Nucleosome positioning predictionBERTLSTMSelf-trainhttps://github.com/FAhtisham/Nucleosome-position-prediction
Dalla et al. (100)Promoter identification, enhancers identification, splice site identification, chromatin accessibility predictionTransformerTransformerSelf-trainhttps://github.com/instadeepai/nucleotide-transformer
Ji et al. (93)Promoter prediction, splice sites prediction, transcription factor binding sites predictionBERTBERTSelf-trainhttps://github.com/jerryji1993/DNABERT
Xu et al. (52)Transcription factor binding affinity predictionBERTBERTSelf-trainhttps://github.com/ericcombiolab/TRAFICA
Kabir et al. (94)Transcription factor binding affinity predictionBERTCNNPretrainedhttps://github.com/lanl/EPBD-BERT
Chen et al. (343)Transcription factor binding affinity predictionGPTRFPretrainedhttps://github.com/yiqunchen/GenePT
Clauwaert et al. (50)Transcription sites prediction, translation initiation sites, methylation sites predictionTransformer-XLCNNSelf-trainhttps://github.com/jdcla/DNA-transformer
Zhang et al. (344)Conserved non-coding element classificationGPTGPTSelf-trainTencentAILabHealthcare/DNAGPT(github.com)
Zhang et al. (344)Translation initiation sites identificationGPTGPTSelf-trainhttps://github.com/TencentAILabHealth-care/DNAGPT
Wang et al. (61)DNA replication origins predictionBERTBERTPretrainedhttps://github.com/CongWang3/PLANNER
Le et al. (60)Enhancer identificationBERTCNNPretrainedhttps://github.com/khanhlee/bert-enhancer
Wang et al. (61)Enhancer identificationBERTDFPretrainedhttps://github.com/no-banana/SMFM-master
Clauwaert et al. (91)Gene functions predictionTransformerCNNPretrainedhttps://github.com/jdcla/DNA-transformer
Martinek et al. (283)Enhancers identification, promoter identificationBERTBERTPretrainedhttps://github.com/ML-Bioinfo-CEITEC/
Li et al. (284)Protein-DNA interface hotspots predictionBERTRFPretrainedhttps://github.com/lixiangli01/PDH-EH
Toufiq et al. (114)Candidate gene prioritization and selectionGPTGPTPretrainedhttps://github.com/Drinchai/A37_LLM

Summary of open-source language model-based models in existing studies.

A close look at Table 4 reveals that 25 AI-driven predictive pipelines are developed for 16 unique DNA sequence analysis tasks. These tasks include disease gene identification, phage-host interactions prediction, nucleosome position detection, enhancer identification, promoter identification, enhancer-promoter interactions prediction, YY1-mediated chromatin loop prediction, methylcytosine site prediction, methyladenine site prediction, essential gene identification, disease gene prediction, pseudogene function prediction, mutation susceptibility analysis, target gene classification, gene taxonomy classification, protein-DNA binding sites identification, and recombination spots identification. In addition, a high-level overview of Table 4 illustrates that a total of 2 node2vec and OPA2Vec word embedding approaches along with MLP and RF classifiers have made their source code publicly available for disease gene identification. In addition, source code of 4 Word2vec and FastText word embedding approach-based predictive pipelines is publicly available for promoter and enhancer identification tasks. Furthermore, two open-source FastText and node2vec word embedding approach-based predictive pipelines are developed for essential gene identification. Moreover, four Word2vec and one FastText word embeddings based predictive pipelines are developed for DNA methylation modification predictive pipelines.

Overall, Table 4 encompasses source codes of 7 unique word embedding approaches (Word2Vec, FastText, GloVe, Node2Vec, OPA2Vec, Graph2Vec, SDNE). Furthermore, a total of 3 machine learning classifiers, namely, RF, SVM, and XGBoost, 4 standalone deep learning classifiers, namely, MLP, CNN, GCN, and BiLSTM, and 6 hybrid deep learning models are used for the development of 26 predictive pipelines for 19 distinct DNA sequence analysis tasks.

Analysis of Table 5 demonstrates that 38 predictive pipelines are developed using 4 unique large language models, namely, BERT, ALBERT, GPT, and Transformer, and 9 unique classifiers, namely, RF, CatBoost, DF, MLP, CRF, CNN, LSTM, FGM, and hybrid (CNN, BiLSTM). Overall, these 38 large language models based predictive models are evaluated across 24 unique DNA sequence analysis tasks. These 24 tasks include chromatin accessibility prediction, species classification, protein-DNA binding site prediction, long-range chromatin interaction prediction, prediction of context-specific functional impact of genetic variants, enhancer-promoter interaction prediction, gene expression prediction, tumor type prediction, methyladenine modification prediction, methylation modification prediction, methylcytosine modification prediction, conserved non-coding element classification, nucleosome position prediction, promoter identification, splice site prediction, transcription factor binding site prediction, transcription factor binding affinity prediction, transcription site prediction, translation initiation site prediction, dna replication origin prediction, enhancer identification, gene function prediction, protein-dna interface hotspots prediction, and candidate gene prioritization and selection. A high level overview of Table 5 reveals that a total of two open-source chromatin accessibility predictive pipelines and two open-source gene expression prediction pipelines use transformers. In contrast, two open-source protein-DNA binding site identification pipelines and two open-source transcription factor binding site identification pipelines use BERT. In addition, three open-source transcription factor binding affinity prediction pipelines use GPT and BERT language models, whereas 8 methyl-adenine and 4 methyl-cytosine modification prediction pipelines use BERT and transformers.

Predictive pipelines can use language models in two different ways: (1) training a language model from scratch (self-training) on a large corpus and (2) leveraging a pre-trained open-source language model and fine-tuning it for specific downstream tasks. Overall, a critical analysis of existing studies reveals that source codes of 20 BERT, 13 Transformer, 4 GPT, and 1 Transformer-XL based predictive pipelines are publicly available. A holistic view of Table 5 reveals that 23 open-source predictive pipelines perform self-training of different language models from scratch for 20 tasks, whereas 15 open-source predictive pipelines have used pre-trained language models for 11 different tasks.

Specifically in 20 BERT-based predictive pipelines, 9 BERT models are self-trained from scratch for nine different tasks, namely, protein-DNA binding site prediction (96), 6mA-methyl adenine modification prediction (144, 281), DNA methylation modification prediction (159, 248), 5mC-methyl cytosine modification prediction (282), nucleosome positioning prediction (29), promoter prediction (93), splice site prediction (93), transcription factor binding site prediction (93), and transcription factor binding affinity prediction (52). In contrast, 11 pre-trained BERT models are utilized to perform 7 downstream tasks, namely, protein-DNA binding site prediction (53), DNA methylation modification (154, 155, 157), 4mC-methyl cytosine modification prediction (143), transcription factor binding affinity prediction (94), DNA replication origin prediction (25), enhancer identification (60, 61, 283), and protein-DNA interface hotspots prediction (284). To facilitate readers, we have summarized uniquely pre-trained language models along with pre-training data for DNA sequence analysis tasks in Table 6.

Table 6

Language modelPre-trained dataLanguage modelPre-trained dataLanguage modelPre-trained data
Zhang et al., Transformer (31)Human reference Genome Sequences GRCh37 dataStanojevic et al., Transformer (152)893k sequences from ONT GM24385 DatasetElnaggar et al. BERT (345)UniRef100 and BFD-100 datasets
Nguyen et al., Transformer (44)Human reference Genome Sequences GRCh37 dataClauwaert et al., Transformer-XL (50)2.7M Genome Data for TSSs, TISs, MethSMRT for the 4mC-MethylationsDevlin et al., BERT (244)BooksCorpus (800M words), English Wikipedia (2,500M words)
Gao et al., Transformer (36)EpiGenomic dataLuo et al., BERT (96)690 ChIP-Seq Dataset (20,464,149 Samples)Lai et al., BERT (346)PubMed abstracts
Ni et al., Transformer (83)Human reference Genome SequencesLe at al., BERT (281)Cased Text in the top 104 Languages with the Largest CorpusYang et al., Transformer, ALBERT (249)Human reference Genome Sequences hg19 data
Reddy et al., Transformer (102)MPRA dataTsukiyama et al., BERT (144)R. chinensisZhang et al., GPT (344)Genomes from 9 species: (Arabidopsis thaliana, Caenorhabditis elegans, Bos taurus, Danio rerio, Drosophila melanogaster, Escherichia coli gca 001721525, Homo sapiens, Mus musculus, Saccharomyces cerevisiae)
Fishman et al., Transformer (32)Human T2T v2 GenomeJeong et al., BERT (159)HG19, MM10Zhang et al., GPT (344)(a) approx.10B bps (b) approx. 200B bps
Osseni et al., Transformer (180)Omics datasetHuson et al., BERT (248)DNA Methylation and taxonomy DataCui et al., GPT (347)NCBI text descriptions of individual genes
Wang et al., Transformer (158)WGBS datasetWang et al., BERT (282)1,825,095 Promoter SequencesToufiq et al., GPT 4 (114)Co–expression gene set (M9.2) from the BloodGen3 repertoire associated with circulating erythroid cells
Zhou et al., Transformer (250)WGBS datasetFazeel et al., BERT (29)Human reference Genome Sequences with the length of sequences between 5 and 510 with 3-merToufiq et al., Claud (114)Co-expression gene set (M9.2) from the BloodGen3 repertoire associated with circulating erythroid cells
Dalla et al., Transformer (100)A total of 850 species, whose Genomes add up to 174B nucleotidesJi et al., BERT (93)Human reference Genome Sequences__
Clauwaert et al., Transformer (91)9 283 204 full genome sequenceXu et al., BERT (52)ATAC-Seq dataset with over 13M nucleotide sequences__

Summary of uniquely pre-trained language models along with pre-training data for DNA sequence analysis tasks.

10 DNA sequence analysis predictive pipeline performance analysis

This section facilitates AI researchers by providing details of performance figures achieved over diverse benchmark datasets for all three kinds of predictive pipelines, namely, word embedding, language models, and nucleotide compositional and positional information across 44 distinct DNA sequence analysis tasks. To assist researchers for developing novel predictive pipelines, we have thoroughly analyzed literature and identified the current state-of-the-art predictor for each task. Section 3 provides categorization of 44 DNA sequence analysis tasks into 8 different categories. In this section, we have summarized the performance values of the predictive pipelines for these 44 tasks into 7 different Tables. Each Table corresponds to the predictive pipelines of tasks within one category, except for 1 Table that includes predictive pipelines related to tasks from 3 different categories. Within these Tables, highlighted predictors represent state-of-the-art performance values on public datasets across each task. Furthermore, this section also facilitates crucial information that which of the tasks of every goal offer more room for improvement through the development of more robust and effective predictive pipelines.

Table 7 summarizes the crucial details of seven DNA sequence analysis tasks classified under the hood of genome structure and stability. Overall, for genome structure and stability goal, four unique representation learning methods, namely, BERT, Transformer, Word2vec, and multi-scale convolution, in conjunction with bi-directional gated recurrent unit methods are used across seven different tasks. Similarly, six unique classifiers, namely, BERT, LSTM, Transformer, CNN+LSTM, CapsNet, and CNN, are used in seven different task predictive pipelines. Most commonly used representation learning scheme for this goal is BERT followed by Transformer. BERT is most commonly used with a self-classifier for three different tasks and used with LSTM classifier for one task. Transformer is used with only self-classifier for three different tasks. Word2vec potential is explored with CNN-based classifiers for two different tasks and multi-scale convolution in conjunction with bi-directional gated recurrent unit method is only explored with CNN classifier for one task. Overall, among all predictive pipelines, BERT with self-classifier or LSTM classifier manages to achieve top performance figures as compared to transformer-based predictive pipelines. Among all 7 tasks, genome structure analysis and long range chromatin interaction prediction tasks provide a lot of room for improvement as the performance of their predictive models fall below 70%. BERT or Transformer with CapsNet classifier-based predictive pipeline can potentially enhance the performance on either or both of these tasks.

Table 7

Task typeTask nameReferencesDatasetRepresentation learningClassifierPerformance evaluation
Binary classificationDNA replication origins prediction(25)Gao et al. Dataset (A. thaliana)BERT_A. thaliana: AUROC = 0.9811
(329)Wu et al. Datasets (S. cerevisiae Dataset, S. pombe Dataset, K. lactis Dataset, P. pastoris Dataset)Word2VecCNNAccuracy [S. cerevisiae (S1): 0.975, S. pombe (S2): 0.765, K. lactis (S3): 0.885, P. pastoris (S4): 0.967]; MCC [S. cerevisiae (S1): 0.940, S. pombe (S2): 0.530, K. lactis (S3): 0.771, P. pastoris (S4): 0.934]; AUC [S. cerevisiae (S1): 0.975, S. pombe (S2): 0.800, K. lactis (S3): 0.888, P. pastoris (S4): 0.981]
Nucleosome position detection(29)Gangi et al. Dataset 1 [a. C. elegans (CE), b. D. melanogester (DM), c. S. cerevisiae (YS), d. H. sapiens(HM)], Gangi et al. Dataset 2 (DM-5U, DM-PM, DM-LC, HM-5U, HM-PM, HM-LC, YS-PM, YS-WG)BERTLSTMDataset 1 a. CE: Acc = 90.5, Sn = 91.8, Sp = 92.1, Precision = 91.8, MCC = 80.5, AUROC = 95.8 b. DM: Acc = 85.1, Sn = 84.8, Sp = 85.6, Precision = 85.3, MCC = 70.5, AUROC = 92.4 c. YS: Acc = 100, Sn = 100, Sp = 99.8, Precision = 99.8, MCC = 99.82, AUROC = 100 d. HM: Acc = 88.3, Sn = 88.3, Sp = 88.4, Precision = 88.5, MCC = 76.8, AUROC = 94.4 Dataset 2 a. DM-5U: Acc = 69.5, Sn = 41.1, Sp = 85.8, Precision = 63.8, MCC = 30.8, AUROC = 68.3 b. DM-PM: Acc = 73.6, Sn = 40.1, Sp = 93.6, Precision = 80.4, MCC = 42.0, AUROC = 73.7 c. DM-LC: Acc = 71.3, Sn = 43.1, Sp = 90.0, Precision = 75.2, MCC = 38.7, AUROC = 72.0 d. HM-5U: Acc = 81.8, Sn = 51.6, Sp = 94.3, Precision = 80.0, MCC = 53.4, AUROC = 80.2 e. HM-LC: Acc = 91.1, Sn = 83.7, Sp = 96.1, Precision = 93.7, MCC = 81.7, AUROC = 95.1 f. HM-PM: Acc = 85.1, Sn = 75.8, Sp = 92.4, Precision = 89.1, MCC = 70.1, AUROC = 90.4 g. YS-PM: Acc = 92.4, Sn = 63.1, Sp = 97.2, Precision = 79.6, MCC = 66.3, AUROC = 93.5 h. YS-WG: Acc = 94.3, Sn = 60.3, Sp = 98.4, Precision = 94.3, MCC = 67.2, AUROC = 94.5
Chromatin accessibility prediction(32)DeepSEA dataset (TF & DHS)BERT_TF: AUROC = 96.81±0.1 DHS: AUROC = 92.8±0.03
(31)DNase-Seq experiment DataTransformer_Average AUROC = 0.8977 Average AUPRC = 0.8983
Binary classificationYY1-mediated chromatin loops prediction(39)1. HCT116 2. K562Word2VecCNN+LSTMK562: AUROC = 98.2, Acc = 92.9 HCT116: AUROC = 95.7, Acc = 88.5
(38)1. HCT116 2. K562MSC+BiGRUCapsNet10-fold cross-validation HCT116: Acc = 0.9544, AUROC = 0.9886 K562: Acc = 0.9680, AUROC = 0.9924 Independent test setting HCT116: Acc = 0.9622, AUROC = 0.9913, AUPRC = 0.9917, F1-score = 0.9564 K562: Acc = 0.9560, AUROC = 0.9912, AUPRC = 0.9920, F1-score = 0.9583
(45)1. HCT116 2. K562Word2VecCNN1. AUROC = 0.93 2. AUROC = 0.93
Multi-class classificationGenome structure analysis(27)NCycDBBERT_Macro F1-score = 61.8, Weighted F1-score = 65.4
Chromatin feature prediction(34)Yang et al. Datasets 1. LOGO-919, 2. LOGO-2002 (TF & DHS), LOGO-3357 (TF & DHS)Transformer_LOGO-919: AUROC = 0.703 LOGO-2002: TF: AUROC = 0.954, DHSs: AUROC = 0.913, HM: AUROC = 0.883 LOGO-3357: TF: AUROC = 0.926, DHSs: AUROC = 0.928, HM: AUROC = 0.883
InteractionLong range chromatin interaction prediction(36)ChIP-seq DatasetTransformer_AUPRC = 0.647

Genome structure and stability related 10 distinct DNA sequence analysis task predictive pipeline performance.

Bold values highlight top performers on unique datasets related to distinct tasks.

Furthermore, Table 8 provides a high-level overview of the performance achieved by 48 predictors for 9 DNA sequence analysis tasks classified under the hood of gene expression regulation. Overall, for gene expression regulation goal, 10 unique representation learning methods, namely, ULMFiT, BERT, One-hot encoding, Word2vec, FastText, C2+NCP, Node2vec+SocDim+Grarep, ELECTRA, ALBERT, and Transformer, are used across 9 different tasks. Overall, 14 unique classifiers are employed in these predictive pipelines, namely, CNN, MLP, DF, CNN+BiLSTM, CNN+LSTM, SVM, Siamese network, DenseNet, CatBoost, TCN, XGBoost, RF+CNN, LSTM, BiGRU, and CNN+BiGRU.

Table 8

Task typeTask nameReferencesDatasetRepresentation learningClassifierPerformance evaluation
Binary classificationEnhancer identification(57)Liu et al. dataset (1, 2)ULMFiTCNN1. Enhancer/ Non-Enhancer Cross-Validation: Acc = 0.946, Sn = 0.946, Sp = 0.949, MCC = 0.892 Independent: Acc = 0.843, Sn = 0.842, Sp = 0.87, MCC = 0.686 2. Weak/ Strong-Enhancer Cross-Validation: Acc = 0.90, Sn = 0.90, Sp = 0.896, MCC = 0.8 Independent: Acc = 0.875, Sn = 0.873, Sp = 0.75, MCC = 0.774
(56)Liu et al. dataset 1BERTCNNSn = 1.00, Sp = 1.00, Acc = 1.00, MCC = 1.00, AUROC = 1.00
(55)1. DiseaseEnhancer 2. EnDisease 3. CancerEnDOne-hot EncodingMLP1. AUROC = 0.9645±0.0057, AUPRC = 0.9647±0.0043, Acc = 0.8986±0.0169, Precision = 0.8765±0.0266, Recall = 0.9290±0.0132, F1-score = 0.9018±0.0169 2. AUROC = 0.9546±0.0036, AUPRC = 0.9474±0.0118, Acc = 0.8959±0.0141, Precision = 0.8583±0.0261, Recall = 0.9469±0.0103, F1-score = 0.9003±0.0182 3. AUROC = 0.9755±0.0026, AUPRC = 0.9673±0.0047, Acc = 0.9306±0.0053, Precision = 0.9373±0.0051, Recall = 0.9261±0.0085, F1-score = 0.9317±0.0053
(67)Liu et al. dataset 1, Basith et al. dataset (HEK293, NHEK, K652, GN12878, HMEC, HSMM, NHLF, HUVEC)BERT_(Liu et al. dataset 1) Acc = 0.8300, Sn = 0.8000, Sp = 0.8600, MCC = 0.6612, AUROC = 0.8560 (Basith's dataset) HEK293: Acc = 0.8732, Sn = 0.8666, Sp = 0.8798, MCC = 0.7283, AUROC = 0.9443 NHEK: Acc = 0.7766, Sn = 0.7229, Sp = 0.8303, MCC = 0.5453, AUROC = 0.8716 K652: Acc = 0.7974, Sn = 0.8180, Sp = 0.7767, MCC = 0.5679, AUROC = 0.8712 GM12878: Acc = 0.8222, Sn = 0.7564, Sp = 0.8879, MCC = 0.6475, AUROC = 0.9179 HMEC: Acc = 0.7645, Sn = 0.7638, Sp = 0.7652, MCC = 0.5068, AUROC = 0.8631 HSMM: Acc = 0.7193, Sn = 0.7191, Sp = 0.7194, MCC = 0.4179, AUROC = 0.7948 NHLF: Acc = 0.7884, Sn = 0.8236, Sp = 0.7532, MCC = 0.5479, AUROC = 0.8623 HUVEC: Acc = 0.7334, Sn = 0.7691, Sp = 0.6977, MCC = 0.4417, AUROC = 0.8045
(61)Liu et al. datasetBERTDFAUROC = 0.808, Acc = 0.822, MCC = 0.655, Sn = 0.834, Sp = 0.810
(58)Liu et al. dataset (1, 2)Word2VecCNN + BiLSTM1: Acc = 83.32, Sn = 84.18, Sp = 82.45, MCC = 0.6668 2: Acc = 83.30, Sn = 89.27, Sp = 77.33, MCC = 0.6736
(65)Liu et al. dataset (1, 2)FastTextLSTM + CNNEnhancer Prediction: Acc = 0.7525, MCC = 0.5051; Enhancer type Prediction: Acc = 0.6972, MCC = 0.3954;
(59)Liu et al. dataset (1, 2)FastTextMLP1: Sn = 85.81, Sp = 86.35, Acc = 86.07, MCC = 0.722; 2: Sn = 69.90, Sp = 69.32, Acc = 69.59, MCC = 0.392
(66)Liu et al. dataset (1, 2)Word2VecCNN1: Acc = 0.784, Sn = 0.811, Sp = 0.758, MCC = 0.567; 2: Acc = 0.749, Sn = 0.961, Sp = 0.537, MCC = 0.505
(60)Liu et al. datasetBERTCNNSn = 80, Sp = 71.2, Acc = 75.6, MCC = 0.514
(47)Liu et al. dataset (1, 2)Word2VecCNN1: Sn = 75.88, Sp = 88.88, Acc = 80.63, MCC = 0.6929, AUROC = 0.8957 2: Sn = 73.64, Sp = 76.80, Acc = 76.43, MCC = 0.4505, AUROC = 0.8109
(64)Liu et al. dataset (1, 2)FastTextSVMCross-validation 1: Sn = 81.1, Sp = 83.5, Acc = 82.3, MCC = 0.65 2: Sn = 75.3, Sp = 60.8, Acc = 68.1, MCC = 0.37 independent 1: Sn = 82, Sp = 76, Acc = 79, MCC = 0.58 2: Sn = 74, Sp = 53, Acc = 63.5, MCC = 0.28
Promoter identification(77)Yang et al. datasetOne-hot encodingSiamese networkAcc = 96.80, Sn = 95.08, Sp = 98.56, MCC = 0.9367
(251)Xiao et al. dataset 1BERT_MCC = 0.81, AUPRC = 0.98
(331)Xiao et al. dataset (1, 2)C2+NCPDenseNet1: Sn = 0.9389±0.0106, Sp = 0.9471±0.0117, Acc = 0.9429±0.0038, MCC = 0.8865±0.0076, AUROC = 0.9774±0.0017, F1-score = 0.9432±0.0040 2: Sn = 0.8711±0.0168, Sp = 0.9211±0.0088, Acc = 0.8967±0.0086, MCC = 0.7947±0.0165, AUROC = 0.9353±0.0042, F1-score = 0.8880±0.0117
(236)Wang et al., Datasets (1-8)Node2Vec+SocDim+GraRepCatBoost1. H. Sapiens-I: Sn = 0.9565, Sp = 0.9782, Acc = 0.9673, MCC = 0.9350, Precision = 0.9772, F1-score = 0.9670, AUROC = 0.9952 2. H. Sapiens-II: Sn = 0.8646, Sp = 0.9849, Acc = 0.9248, MCC = 0.8558, Precision = 0.9829, F1-score = 0.9199, AUROC = 0.9844 3. R. Norvegicus-I: Sn = 0.9113, Sp = 0.9757, Acc = 0.9424, MCC = 0.8908, Precision = 0.9761, F1-score = 0.9425, AUROC = 0.9975 4. R. Norvegicus-II: Sn = 0.8905, Sp = 0.9849, Acc = 0.9377, MCC = 0.8793, Precision = 0.9833, F1-score = 0.9346, AUROC = 0.9832 5. D. melanogaster-I: Sn = 0.9350, Sp = 0.9670, Acc = 0.9560, MCC = 0.9123, Precision = 0.9662, F1-score = 0.9555, AUROC = 0.9923 6. D. melanogaster-II: Sn = 0.9425, Sp = 0.7819, Acc = 0.8596, MCC = 0.7281, Precision = 0.8112, F1-score = 0.8697, AUROC = 0.9448 7. Z. mays-I: Sn = 0.934, Sp = 0.9450, Acc = 0.9395, MCC = 0.8791, Precision = 0.9444, F1-score = 0.9392, AUROC = 0.9841 8. Z. mays-II: Sn = 0.9516, Sp = 0.7806, Acc = 0.8661, MCC = 0.7433, Precision = 0.8127, F1-score = 0.8767, AUROC = 0.9485
(235)Xiao et al. dataset 1Word2VecTCNAcc = 91.86, Sn = 92.74, Sp = 91
(348)Ji et al. DatasetBERT_Acc = 0.8613, AUROC = 0.9354, MCC = 0.7226, Precision = 0.8569, Recall = 0.8624
(349)Liu et al. dataset (1, 2)BERTXGBoost1. Promoter Identification: Sn = 84.34, Sp = 86.56, Acc = 85.45 2. Promoter strength classification: Sn = 70.85, Sp = 81.63, Acc = 76.92
(89)Human datasetELECTRA_Acc = 0.862, AUROC = 0.935, F1-score = 0.862, MCC = 0.725, Precision = 0.863, Recall = 0.862
(76)1. K562 2. GM12878 3. HeLa-S3 4. HUVECWord2VecRF + CNN1. AUROC = 0.9809, Acc = 0.9415, MCC = 0.8831 2. AUROC = 0.9783, Acc = 0.9334, MCC = 0.8668 3. AUROC = 0.9824, Acc = 0.9374, MCC = 0.8749 4. AUROC = 0.9847, Acc = 0.9481, MCC = 0.8963
(48)Xiao et al. datasetWord2VecLSTMAcc = 90.59, MCC = 0.8114, Sn = 90.28, Sp = 90.94
(249)Yang et al. datasetALBERT_AUROC = 0.743
(90)Ji et al. datasetBERTCNNPrecision = 0.805, Recall = 0.803, Acc = 0.894
(90)Ji et al. datasetBERT_Precision = 0.805, Recall = 0.803, Acc = 0.894
(34)Yang et al. datasetTransformer_Recall = 0.921, Precision = 0.940, F1-score = 0.933
(93)Human datasetBERT_AUROC = 0.981, AUPRC = 0.982
(79)Xiao et al. dataset (1, 2)Word2VecCNN1: Acc = 91.42 2: Acc = 82.42
(78)Xiao et al. dataset (1, 2)FastTextCNN1: Acc = 85.41 2: Acc = 73.1
Transcription sites prediction(50)Clauwaert et al. DataesetTransformer-XLCNNAUROC = 0.977
Transcription factor binding sites prediction(94)ChIP-Seq DatasetBERTCNNAUROC = 0.949, AUPRC = 0.326
(89)690 ChIP-Seq datasetELECTRA_Acc = 0.856, AUROC = 0.935, F1-score = 0.851, MCC = 0.727, Precision = 0.859, Recall = 0.856
(90)TF ChIP-Seq datasetBERT_Precision = 0.937, Recall = 0.935, Acc = 0.989
(91)TSSs datasetTransformerCNNAUROC = 0.981. AUPRC = 0.141
(93)497 TF ChIP-Seq datasetBERT_Mean Acc = 0.903, Mean AUROC = 0.954, Mean F1-score = 0.901, Mean MCC = 0.807, Mean Precision = 0.898, Mean Recall = 0.909
(92)1. A549 dataset 2. MCF-7 dataset 3. H1-HESC dataset 4. HUVEC datasetWord2VecBiGRUHUVEC: Average Precision = 0.9618, Average AUROC = 0.9608 MCF7: Average Precision = 0.9653, Average AUROC = 0.9643 A549: Average Precision = 0.9608, Average AUROC = 0.9593 H1-HESC: Average Precision = 0.9528, Average AUROC = 0.9524
Splice sites prediction(97)Splice-junction Gene sequence dataset_LSTMAcc = 0.98, Precision = 0.98, Recall = 0.99, F1-score = 0.98
(98)Degroeve et al. Dataset_CNNF1-score = 0.9187±0.0070, MCC = 0.9028±0.0085
(99)Liu et al. DatasetOne-hot EncodingMLPDonor: Acc = 96.57 Acceptor: Acc = 95.82
(100)Wang et al. DatasetTransformer_AUPRC = 0.984
(93)Ji et al. DatasetBERT_Acc = 0.939, AUROC = 0.992, F1-score = 0.937, MCC = 0.903, Precision = 0.961, Recall = 0.919
Translation initiation site prediction(101)Kalkatawi et al. TIS datasetk-mer EmbeddingBi-GRUAcc = 96.40, AUROC = 96.40, AUPRC = 94.87
(50)TIS datasetTransformer-XLCNNAUROC = 0.998
InteractionEnhancer-promoter interactions prediction(350)Whalen et al. datasets (GM12878, HUVEC, HeLa-S3, IMR90, K562, NHEK)GCN, K-merGCNF1-score: GM12878 = 0.8679, HUVEC = 0.8954, HeLa-S3 = 0.9175, IMR90 = 0.7949, NHEK = 0.6085
(233)Whalen et al. datasets (HMEC, IMR90, K562, NHEK)TransformerMulti-Scale CNNHMEC: AUROC = 0.9344, AUPRC = 0.6852, IMR90 AUROC = 0.8936, AUPRC = 0.5878, K562 AUROC = 0.8513, AUPRC = 0.2101, NHEK AUROC = 0.8243, AUPRC = 0.4760
(332)Zheng et al. Datasets (GM12878, HeLa)DeeptoolsRF1. Sn = 0.578, Sp = 0.964, Precision = 0.799, Acc = 0.887, AUROC = 0.919, AUPRC = 0.773, 2. Sn = 0.363, Sp = 0.953, Precision = 0.660, Acc = 0.836, AUROC = 0.831, AUPRC = 0.601
(177)ESKAPE datasetSDNE, Word2VecMLPAcc = 86.65 ± 1.55, Sn = 88.40 ± 1.81, Sp = 84.91 ± 1.96, Precision = 85.43 ± 1.74, F1-score = 86.88 ± 1.53, AUROC = 0.9208 ± 0.0119
(249)Yang et al. Datasets (FoeT, Mon, nCD4, tB, tCD4, tCD8)TransformerMLPFoeT: AUPRC = 0.9447, Mon: AUPRC = 0.9414, nCD4: AUPRC = 0.9457, tB: AUPRC = 0.9474, tCD4: AUPRC = 0.9475, tCD8: AUPRC = 0.9387
(82)Whalen et al. Datasets (GM12878, HUVEC, HeLa-S3, IMR90, K562, NHEK)Word2VecCNN + BiGRUGM12878: AUROC = 0.965, AUPRC = 0.819, HUVEC: AUROC = 0.950, AUPRC = 0.773, HeLa-S3: AUROC = 0.960, AUPRC = 0.820, IMR90: AUROC = 0.962, AUPRC = 0.801, K562: AUROC = 0.959, AUPRC = 0.814, NHEK: AUROC = 0.985, AUPRC = 0.899
(49)Whalen et al. Datasets (GM12878, HUVEC, HeLa-S3, IMR90, K562, NHEK)Word2VecCNN+BiGRUGM12878: AUROC = 0.949, AUPRC = 0.819, F1-score = 0.766, HUVEC: AUROC = 0.948, AUPRC = 0.720, F1-score = 0.649, HeLa-S3: AUROC = 0.952, AUPRC = 0.824, F1-score = 0.78, IMR90: AUROC = 0.948, AUPRC = 0.818, F1-score = 0.778, K562: AUROC = 0.955, AUPRC = 0.826, F1-score = 0.795, NHEK: AUROC = 0.977, AUPRC = 0.893, F1-score = 0.861
Protein-DNA binding sites prediction(53)Patiyal et al. dataset (dataset 1) Xia et al. dataset (dataset 2)BERTCNNdataset 1: TE46 Sp = 0.835, Recall = 0.747, Precision = 0.306, F1-score = 0.434, MCC = 0.401, AUROC = 0.871 dataset 2: TE129 Sp = 0.955, Recall = 0.464, Precision = 0.396, F1-score = 0.427, MCC = 0.389, AUROC = 0.881
(96)690 ChIP-Seq DatasetBERT_AUROC = 0.947±0.041, Acc = 0.880±0.062, Precision = 0.882±0.061, Recall = 0.880±0.062, F1-score = 0.880±0.062, MCC = 0.762±0.122
(95)Liu et al., Dataset, Tian et al., DatasetBERTRFdataset 1: Sp = 0.529, Precision = 0.106, Recall = 0.574, F1-score = 0.179, AUROC = 0.551, MCC = 0.025 dataset 2: Sp = 0.724, Precision = 0.119, Recall = 0.536, F1-score = 0.194, AUROC = 0.630, MCC = 0.067
RegressionTranscription factor binding affinity prediction(52)1. Weirauch et al. dataset, 2. Jolma et al. datasetBERT_1: Average PCC = 0.649 2: Average R2 = 0.977

Gene expression regulation related 48 distinct DNA sequence analysis task predictive pipeline performance.

Bold values highlight top performers on unique datasets related to distinct tasks.

For this goal, most commonly used representation learning approach is BERT. BERT is used with five different classifiers across six unique tasks. BERT with self-classifier is evaluated across all six unique tasks, whereas BERT with other four classifiers is evaluated on some of these six tasks. Specifically, BERT with CNN classifier is evaluated on four common tasks, namely, enhancer identification, promoter identification, protein-DNA binding site prediction, and transcription factor binding site prediction. BERT with DF, RF, and XGBoost is evaluated on one common task each including enhancer identification, protein-DNA binding site prediction, and promoter identification. Among all BERT-based predictive pipelines, BERT achieves the best performance with CNN classifier on two tasks, namely, transcription factor binding site prediction and protein-DNA binding site prediction. Second most common representation learning approach for this goal is Word2vec that is explored with seven different classifiers for four different tasks. Specifically, Word2Vec with CNN and CNN+BiLSTM is used for one task, with LSTM, TCN, and RF+CNN for one task, with BiGRU for one task, and with CNN+BiGRU for one task. Among all Word2vec-based predictive pipelines, Word2Vec achieves the best performance with CNN+BiGRU classifier on enhancer-promoter interaction prediction task. From two most common approaches, BERT with CNN and self-classifiers manages to yield top performance values as compared to Word2vec-based predictive pipelines. Beyond BERT and Word2vec, transformer-XL is used with CNN for two different tasks, ALBERT and ELECTRA are used with self-classifiers for a single task, and ULMFiT is used with CNN for a single task. In addition, potential of transformer is explored with CNN for one task and with a self-classifier for two tasks. FastText-based representation learning is used with MLP and SVM classifiers for a single task and potential of three graph embedding, namely, Node2Vec, SocDim, and GraRep, is explored with CatBoost classifier for a single task. Overall, among all approaches, ULMFiT manages to achieve best performance with CNN classifier on enhancer identification task. From all nine tasks, protein-DNA binding site prediction and protein-DNA binding affinity prediction have some room for improvement. Considering the promising performance trends for this goal, Word2vec potential can be explored with CNN+BiGRU classifier and ULMFiT potential can be explored with standalone CNN as well as ensemble CNN+BiGRU classifier to enhance the performance of under-performing tasks.

In addition, Table 9 summarizes predictive models developed for seven different DNA sequence analysis tasks classified under the hood of gene analysis. For gene analysis goal, 12 unique representation learning methods are used that include Gapped K-mer Encoding, ESM-2, Flux Sampling, Node2Vec, FastText, OPA2Vec, Transformer, Laplacian eigenmaps + Locally linear Embedding + DeepWalk + Node2Vec, and GPT. Overall, eight unique classifiers, namely, GCNN, GNN, GAT, kNN, RF, SVM, CNN, and MLP, are used in different predictive pipelines. Most commonly used representation learning scheme for this goal is Node2vec followed by FastText. Node2vec is used with GCNN classifier for three different tasks and with GAT and MLP classifiers for two different tasks, whereas FastText is used with an ensemble and MLP classifiers for two different tasks. From most commonly used approaches, Node2Vec performs better in the majority of tasks and achieve top performance with GAT classifier as compared to FastText approach. Apart from Node2Vec and FastText, potential of ESM-2 along with a self-classifier and flux sampling with GNN is explored for a single task, transformer with CNN is evaluated for two different tasks, and GPT is used with a self-classifier for two tasks. In addition, OPA2Vec is used with RF for one task, Gapped K-mer encoding with GCN classifier is used for one task, and potential of Laplacian eigenmaps, Locally linear Embedding, DeepWalk, and Node2Vec is explored with RF for one task. In the holistic view, among all approaches, Gapped K-mer Encoding approach with GCNN classifier achieves the best performance for essential gene identification task. Among all seven distinct tasks, multi-label classification task, namely, gene function prediction provides a lot of room for improvement as the performance of its respective predictive pipeline based on Node2vec and GCNN classifier fall below 60%. Considering the promising performance of Gapped K-mer Encoding method, Gapped K-mer Encoding method and GCNN-based predictive pipeline can prove fruitful for various low performance tasks such as gene function prediction.

Table 9

Task typeTask nameReferencesdatasetRepresentation learningClassifierPerformance evaluation
Binary classificationEssential genes identification(333)Hu et al. dataset (D. melanogaster, M. maripaludis, H. sapiens, C. elegans)Gapped K-mer encodingGCNN1. Sn = 0.8333, Sp = 0.9939, Acc = 0.9847, MCC = 0.8545, AUROC = 0.8283 2. Sn = 0.9052, Sp = 0.9304, Acc = 0.9221, MCC = 0.8265, AUROC = 0.8422 3. Sn = 0.9048, Sp = 0.9566, Acc = 0.9501, MCC = 0.7961, AUROC = 0.8655 4. Sn = 0.8362, Sp = 0.9368, Acc = 0.9242, MCC = 0.6983, AUROC = 0.7834
(109)Ma et al. dataset (S.cerevisiae, E.coli, H.sapiens, D.melanogaster)ESM-2__
(104)Campos et al. datasetFlux samplingGNNAcc = 0.871±0.012, Precision = 0.769±0.030, Recall = 0.718±0.037, F1-score = 0.743±0.023
(108)Ma et al. dataset (S. cerevisiae, E. coli, H. sapiens, D. melanogaster)Node2VecGATAUROC S. cerevisiae: (DIP: 76.57±0.74, BioGrid: 87.66±2.58, STRING: 90.13±1.08); E.coli: (DIP: 79.96±2.20, BioGrid: 92.35±1.15, STRING: 97.02±0.50); H.sapiens: (DIP: 75.61±0.90, BioGrid: 88.39±0.52, STRING: 90.95±0.54); D.melanogaster: (DIP: 32.90±5.34, BioGrid: 78.78±4.24, STRING: 77.82±1.95)
(105)Hu et al. datasetFastTextkNN + RF + SVM + CNNSn = 60.2, Sp = 84.6, Acc = 76.3, MCC = 0.449, AUROC = 0.814
(106)Zhang et al. datasetNode2VecMLPAUROC = 94.15, Sp = 94.75, Average Precision = 90.64, Acc = 90.88
(107)Xiao et al. datasetNode2VecMLPAverage AUROC = 95.17, AUPRC = 92.21, Acc = 91.59, F1-score = 78.71
Binary classificationTarget gene classification(113)Argoty et al. datasetFastTextMLPPrecision = 0.99, Recall = 0.99
Binary classificationDisease genes prediction(110)Nunes et al. datasetOPA2VecRFMedian WAF score = 0.768
Multi-label classificationGene function prediction(112)(1). 1,000 human gene set from gene ontology, (2). 100 omic gene setGPT_Semantic similarity = 0.50, Gene covered = 30
Multi-class classificationGene expression prediction(102)Reddy et al. dataset (Jurkat, K-562, THP-1)TransformerCNN(1). Jurkat PCC = 0.6389±0.0036, SRCC = 0.5996±0.0093 (2). K-562 PCC = 0.6152±0.0082, SRCC = 0.6043±0.0045 (3). THP-1 PCC = 0.5672±0.0131, SRCC = 0.4742±0.0136
(103)Al Taweraqi et al. datasetLaplacian eigenmaps + Locally linear Embedding + DeepWalk + Node2VecRF_
Multi-class classificationPseudo-gene function prediction(111)Fan et al. dataset (CC, MF, BP)Node2VecGCNCC: (AUPRC = 0.587 ± 0.02, F1-score = 0.380 ± 0.01) MF: (AUPRC = 0.463 ± 0.02, F1-score = 0.319 ± 0.01) BP: (AUPRC = 0.362 ± 0.01, F1-score = 0.193 ± 0.01)
Multi-class classificationCandidate gene prioritization & Selection(114)Toufiq et al. datasetGPT__

Gene analysis related 15 distinct DNA sequence analysis task predictive pipeline performance.

Bold values highlight top performers on unique datasets related to distinct tasks.

Moreover, Table 10 summarizes eight DNA sequence analysis tasks classified across three unique biological goals, namely, DNA modification prediction, environmental and microbial genomics, and gene network analysis. For DNA modification prediction goal, across 4 DNA sequence analysis tasks, 27 predictors are developed. In the predictive pipelines, overall 10 unique representation learning methods are used which include PSeKNC, BERT, nucleotide physico-chemical properties and occurrence frequency based encoder, Transformer-XL, Word2Vec, One-hot encoding, FastText, ULMFIT, and BERT+ ALBERT+XLNet+ ELECTRA. Similarly, nine unique classifiers, namely, Structural Sparse Regularized Random Vector Functional Link Network, CatBoost, KNN, CNN, BiLSTM, CNN+BiLSTM, SVM, XGBoost, and FGM, are employed in different predictive pipelines. For DNA modification prediction goal, most commonly used representation learning approach is BERT. BERT is used with five different classifiers for all four tasks. Specifically, BERT with a self-classifier is evaluated for three tasks, namely, 4mC-methyl cytosine, 5mC-methyl cytosine, and DNA methylation modification prediction. BERT is also used with two other classifiers, namely, CatBoost, and FGM, for two common tasks, namely, 4mC-methyl cytosine modification prediction, and DNA methylation modification prediction. In addition, BERT is used with CNN and CNN+BiLSTM classifier for one task, namely, 6mA-methyl adenine modification prediction. Among all BERT-based predictive pipelines, BERT with a self-classifier showed top performance values for two tasks, namely, 5mC-methyl cytosine modification prediction, and DNA methylation modification prediction.

Table 10

Task typeTask nameReferencesDatasetRepresentation learningClassifierPerformance evaluation
Goal: DNA modification prediction
Binary classification4-Methylcytosine (4mc) modification prediction(134)Chen et al. datasets (1-6)PseKNCStructural Sparse Regularized Random Vector Functional Link Network1. Acc = 0.8761, Sn = 0.8630, Sp = 0.8893, MCC = 0.7530 2. Acc = 0.8753, Sn = 0.8739, Sp = 0.8768, MCC = 0.7512 3. Acc = 0.8278, Sn = 0.8256, Sp = 0.8301, MCC = 0.6566 4. Acc = 0.9601, Sn = 0.8641, Sp = 0.9562, MCC = 0.9210 5. Acc = 0.9011, Sn = 0.8895, Sp = 0.9127, MCC = 0.8031 6. Acc = 0.9139, Sn = 0.9087, Sp = 0.9190, MCC = 0.8289
(143)Yang et al. dataset (A. thaliana, C.elegans, D. melanogaster, E. coli, G. pickeringii, G. subterraneous)BERTCatBoostA. thaliana: MCC = 0.6954, F1-score = 0.8521 C. elegans: MCC = 0.8326, F1-score = 0.9183 D. melanogaster: MCC = 0.7924, F1-score = 0.8986 E. coli: MCC = 0.9356, F1-score = 0.9679 G. pickeringii: MCC = 0.8904, F1-score = 0.9431 G. subterraneous: MCC = 0.8796, F1-score = 0.9363
(351)Chen et al. datasetBERT_AUROC = 0.897, AUPRC = 0.907
(135)Khanal et al. dataset (C. elegans, D. melanogaster, A. thaliana, E. coli, G. subterraneus, G. pickeringi, F. vesca, R. chinensis)Number of codons, Number of occurrences of each codon, Proportion of each codon, Number of Nucleotides, Average number of Nucleotides per codon, Percentage of GC, Percentage of purines AG, Percentage of pyrimidines CT, Percentage of AT, Molecular weight of the sequence, Melting temperature, Proportion of Nucleotide DNA sequence, Protein sequence from DNA sequence, Number of amino acids, Percentage of amino acids, Aromaticity, Instability index, Isoelectric point, Molecular weight of Portion, GravyKNNC.elegans: Acc = 92.20, AUROC = 91.99, Precision = 89.47, Recall = 95.67D. melanogaster: Acc = 92.79, AUROC = 92.80, Precision = 88.54, Recall = 98.30A. thaliana: Acc = 90.27, AUROC = 90.28, Precision = 87.52, Recall = 93.93E. coli: Acc = 91.02, AUROC = 91.03, Precision = 86.36, Recall = 97.43G. subterraneus: Acc = 93.09, AUROC = 93.09, Precision = 91.48, Recall = 95.02G. pickeringi: Acc = 90.78, AUROC = 90.79, Precision = 87.20, Recall = 95.61 F. vesca: Acc = 90.67, AUROC = 90.68, Precision = 85.31, Recall = 98.26R. chinensis: Acc = 91.87, AUROC = 91.88, Precision = 87.35, Recall = 97.93
(136)Zulifiqar et al. datasetWord2VecCNNAcc = 0.946, Sn = 0.938, Sp = 0.881, MCC = 0.778, AUROC = 0.989
(50)Clauwaert et al. datasetTransformer-XLCNNAUROC = 0.985
(138)Khanal et al. dataset (F. vesca, R. chinensis)Word2VecCNNF. vesca: Sn = 0.8976, Sp = 0.8417, Acc = 0.8697, MCC = 0.7407, AUROC = 0.9400 R. chinensis: Sn = 0.8219, Sp = 0.8854, Acc = 0.8541, MCC = 0.7093, AUROC = 0.9370
(137)Zeng et al. datasetWord2VecCNNAcc = 0.9321, MCC = 0.8559, Sn = 0.9508, Sp = 0.9161, AUROC = 0.9712
Methyladenine (6ma) modification Prediction(149)1. A. thaliana2. D. melanogaster3. 6mA-rice-Chen 4. 6mA-rice-Lv 5. RosaceaeOne-hot EncodingBiLSTM1. Sn = 0.896, Sp = 0.935, Acc = 0.915, MCC = 0.831, AUROC = 0.967 2. Sn = 0.903, Sp = 0.952, Acc = 0.927, MCC = 0.855, AUROC = 0.963 3. Sn = 0.850, Sp = 0.917, Acc = 0.882, MCC = 0.763, AUROC = 0.947 4. Sn = 0.947, Sp = 0.930, Acc = 0.938, MCC = 0.877, AUROC = 0.976 5. Sn = 0.962, Sp = 0.961, Acc = 0.962, MCC = 0.924, AUROC = 0.990
(150)1. Homo sapiens dataset (Train, Independent) 2. Mus musculus datasetTransformerCNNHomo sapiens (Train): Acc = 96.5 Homo sapiens (Independent): Acc = 93.75 Mus musculus: Acc = 96.86
(144)Lv et al. dataset (A. thaliana, C. elegans, C. equisetispolia, D. melanogaster, F. vesva, H. sapiens, R. chinensis, S. cerevisiae, T. thermophilus, Ts. SUP5-1, Xoc. BLS256)BERTCNN + BiLSTMA. thaliana: AUROC = 0.927 C. elegans: AUROC = 0.962 C. equisetifpolia: AUROC = 0.800 D. melanogaster: AUROC = 0.967 F. vesca: AUROC = 0.976 H. sapiens: AUROC = 0.963 R. chinensis: AUROC = 0.876 S. cerevisiae: AUROC = 0.892 T. thermophilus: AUROC = 0.938 Ts. SUP5-1: AUROC = 0.829 Xoc. BLS256: AUROC = 0.937
(281)DNA 6 mA datasetBERTCNNCross-Validation Sn = 86.4, Sp = 68.8, Acc = 77.6, MCC = 0.651 Independent Sn = 84.3, Sp = 73.1, Acc = 79.3, MCC = 0.580
(145)A. thaliana datasetTransformer_Acc = 0.9633, Sn = 0.9655, Sp = 0.9611, MCC = 0.9266
(151)1. Rice dataset 2. Mus musculus datasetWord2VecBiLSTMRice dataset: Sn = 95.66, Sp = 92.38, Acc = 94.02, MCC = 0.88, AUROC = 0.981 Mus musculus: Sn = 93.28, Sp = 100, Acc = 96.73, MCC = 0.93
(234)Chen et al. datasetFastTextSVMSn = 86.48, Sp = 89.09, Acc = 87.78, MCC = 0.756
5-methylcytosine (5mc) modification prediction(282)Wang et al. datasetBERT_Acc = 0.932, MCC = 0.653, AUROC = 0.966
(153)Wang et al. datasetFastTextXGBoostAcc = 91.8, Sp = 92.0, Sn = 89.9, MCC = 0.626, AUROC = 0.962
(152)Stanojevic et al. datasets (GM24385, NA12878, NA19240, HIESc, K562, HX1)Transformer_GM24385: Acc = 0.9988, Precision = 0.9990, Recall = 0.9988, FPR = 0.0011, F1-score = 0.9989 NA12878: Acc = 0.9936, Precision = 0.9921, Recall = 0.9953, FPR = 0.0081, F1-score = 0.9937 NA19240: Acc = 0.9872, Precision = 0.9678, Recall = 0.9765, FPR = 0.0096, F1-score = 0.9721 H1ESc: Acc = 0.9938, Precision = 0.9995, Recall = 0.9935, FPR = 0.0040, F1-score = 0.9965 K562: Acc = 0.9972, Precision = 0.9512, Recall = 0.9964, FPR = 0.0028, F1-score = 0.9733 HX1: Acc = 0.9950, Precision = 0.9993, Recall = 0.9951, FPR = 0.0057, F1-score = 0.9972
Methylation modification prediction(159)Jeong et al. datasetBERT_Precision = 0.98
(351)Yu et al. datasetBERT_AUROC = 0.897, AUPRC = 0.907
(155)Lv et al. datasets (5hmc: M. musculus, H. sapiens, 4mc: C. equisetifolia, F. vesca, S. cerevisiae, Tolypocladium, 6mA: A. thaliana, C. elegans, C. equisetifolia, D. melanogaster, F. vesca, H. sapiens, R. chinensis, S. cerevisiae, T. thermophile, Tolypocladium, XocBLS256)BERT_5hmC_M. sapiens: Acc = 0.949, AUROC = 0.967, MCC = 0.900 5hmC_M. musculus: Acc = 0.968, AUROC = 0.981, MCC = 0.936 4mC_C. equisetifolia: Acc = 0.853, AUROC = 0.896, MCC = 0.706 4mC_F. vesca: Acc = 0.853, AUROC = 0.928, MCC = 0.706 4mC_S. cerevisiae: Acc = 0.710, AUROC = 0.776, MCC = 0.423 4mC_Tolypocladium: Acc = 0.743, AUROC = 0.819, MCC = 0.487 6mA_A. thaliana: Acc = 0.861, AUROC = 0.934, MCC = 0.722 6mA_C. elegans: Acc = 0.909, AUROC = 0.966, MCC = 0.818 6mA_C. equisetifolia: Acc = 0.745, AUROC = 0.816, MCC = 0.494 6mA_D. melanogaster: Acc = 0.923, AUROC = 0.971, MCC = 0.846 6mA_F. vesca: Acc = 0.939, AUROC = 0.981, MCC = 0.878 6mA_H. sapiens: Acc = 0.907, AUROC = 0.969, MCC = 0.815 6mA_R. chinensis: Acc = 0.818, AUROC = 0.881, MCC = 0.635 6mA_S. cerevisiae: Acc = 0.827, AUROC = 0.905, MCC = 0.654 6mA_T. thermophile: Acc = 0.882, AUROC = 0.944, MCC = 0.772 6mA_Tolypocladium: Acc = 0.768, AUROC = 0.845, MCC = 0.538 6mA_XocBLS256: Acc = 0.877, AUROC = 0.949, MCC = 0.756
(158)CCLE datasetTransformer_Sn = 0.831, Sp = 0.991, Acc = 0.978, MCC = 0.871, AUROC = 0.989
(156)Lv et al. datasets (5hmc:M. musculus, H. sapiens, 4mc:C. equisetifolia, F. vesca, S. cerevisiae, Tolypocladium, 6mA:A. thaliana, C. elegans, C. equisetifolia, D. melanogaster, F. vesca, H. sapiens, R. chinensis, S. cerevisiae,T. thermophile, Tolypocladium, XocBLS256BERT + ALBERT + XLNet + ELECTRA_(6mA)T. thermophile: AUROC = 0.9467, Acc = 0.8840, F1-score = 0.8923, Recall = 0.9611, AUPRC = 0.9321 A. thaliana: AUROC = 0.9378, Acc = 0.8649, F1-score = 0.8615, Recall = 0.8401, AUPRC = 0.9423H. sapiens: AUROC = 0.9687, Acc = 0.9077, F1-score = 0.9068, Recall = 0.8975, AUPRC = 0.9721 Xoc. BLS256: AUROC = 0.9446, Acc = 0.8742, F1-score = 0.8712, Recall = 0.8511, AUPRC = 0.9421D. melanogaster: AUROC = 0.9730, Acc = 0.9276, F1-score = 0.9275, Recall = 0.9258, AUPRC = 0.9761 C. elegans: AUROC = 0.9684, Acc = 0.9131, F1-score = 0.9138, Recall = 0.9219, AUPRC = 0.9674 C. equisetifolia: AUROC = 0.8350, Acc = 0.7590, F1-score = 0.7481, Recall = 0.7158, AUPRC = 0.8492 S. cerevisiae: AUROC = 0.9082, Acc = 0.8325, F1-score = 0.8233, Recall = 0.7802, AUPRC = 0.9198 Tolypocladium: AUROC = 0.8669, Acc = 0.7895, F1-score = 0.7824, Recall = 0.7567, AUPRC = 0.8730 F. vesca: AUROC = 0.9821, Acc = 0.9407, F1-score = 0.9403, Recall = 0.9336, AUPRC = 0.9831 R. chinensis: AUROC = 0.9654, Acc = 0.9164, F1-score = 0.9167, Recall = 0.9197, AUPRC = 0.9691 (4mC) C. equisetifolia: AUROC = 0.9108, Acc = 0.8333, F1-score = 0.8272, Recall = 0.7978, AUPRC = 0.9221 F. vesca: AUROC = 0.9256, Acc = 0.8522, F1-score = 0.8554, Recall = 0.8739, AUPRC = 0.9144 S. cerevisiae: AUROC = 0.8064, Acc = 0.7376, F1-score = 0.7253, Recall = 0.6926, AUPRC = 0.8215 Tolypocladium: AUROC = 0.8149, Acc = 0.7380, F1-score = 0.7285, Recall = 0.7031, AUPRC = 0.80889 (5hmC)M. musculus: AUROC = 0.9817, Acc = 0.9649, F1-score = 0.9651, Recall = 0.9685, AUPRC = 0.9782H. sapiens: AUROC = 0.9680, Acc = 0.9484, F1-score = 0.9500, Recall = 0.9787, AUPRC = 0.9485
(157)Lv et al. datasets (5hmc: M. musculus, H. sapiens, 4mc: C. equisetifolia, F. vesca, S. cerevisiae, Tolypocladium, 6mA: A. thaliana, C. elegans, C. equisetifolia, D. melanogaster, F. vesca, H. sapiens, R. chinensis, S. cerevisiae, T. thermophile, Tolypocladium, XocBLS256)BERTFGM5hmC_H. sapiens: Acc = 0.9501, Sn = 0.9838, Sp = 0.9164, AUROC = 0.9501, MCC = 0.9022 5hmC_M. musculus: Acc = 0.9679, Sn = 0.969, Sp = 0.9668, AUROC = 0.9679, MCC = 0.9358 4mC_C. equisetifolia: Acc = 0.8579, Sn = 0.8743, Sp = 0.8415, AUROC = 0.8579, MCC = 0.7162 4mC_F. vesca: Acc = 0.8524, Sn = 0.8535, Sp = 0.8512, AUROC = 0.8524, MCC = 0.7047 4mC_S. cerevisiae: Acc = 0.723, Sn = 0.6876, Sp = 0.7583, AUROC = 0.723, MCC = 0.447 4mC_Tolypocladium: Acc = 0.7434, Sn = 0.7385, Sp = 0.7483, AUROC = 0.7434, MCC = 0.4868 6mA_A. thaliana: Acc = 0.8603, Sn = 0.8264, Sp = 0.8942, AUROC = 0.8603, MCC = 0.7223 6mA_C. elegans: Acc = 0.9138, Sn = 0.9256, Sp = 0.902, AUROC = 0.9138, MCC = 0.8279 6mA_C. equisetifolia: Acc = 0.7399, Sn = 0.6713, Sp = 0.8084, AUROC = 0.7399, MCC = 0.4843 6mA_D. melanogaster: Acc = 0.9228, Sn = 0.9301, Sp = 0.9155, AUROC = 0.9228, MCC = 0.8457 6mA_F. vesca: Acc = 0.9413, Sn = 0.9452, Sp = 0.9375, AUROC = 0.9413, MCC = 0.8827 6mA_R. chinensis: Acc = 0.8629, Sn = 0.8328, Sp = 0.893, AUROC = 0.8629, MCC = 0.7271 6mA_S. cerevisiae: Acc = 0.8278, Sn = 0.7966, Sp = 0.859, AUROC = 0.8278, MCC = 0.6569 6mA_T. thermophile: Acc = 0.8804, Sn = 0.9442, Sp = 0.8167, AUROC = 0.8804, MCC = 0.7671 6mA_Tolypocladium: Acc = 0.7771, Sn = 0.7649, Sp = 0.7892, AUROC = 0.7771, MCC = 0.5543 6mA_Xoc BLS256: Acc = 0.8817, Sn = 0.8808, Sp = 0.8827, AUROC = 0.8817, MCC = 0.7634
(250)DNAm dataset (Brain, Blood, Buccal, Saliva)Transformer_SRCC: (Brain: Mean = 0.82, SD = 0.004) (Blood: Mean = 0.78, SD = 0.005) (Buccal: Mean = 0.79, SD = 0.007) (Saliva: Mean = 0.79, SD = 0.010) Mean squared error: (Brain: Mean = 0.030, SD = 0.0026) (Blood: Mean = 0.043, SD = 0.0023) (Buccal: Mean = 0.049, SD = 0.0080) (Saliva: Mean = 0.040, SD = 0.055)
(146)1. Maize 5mC dataset 2. Nipponbare 5mC dataset 3. 6mA Chen dataset 4. 6mA Lv datasetULMFiT_1. Maize 5mC Acc = 0.9524, AUNP = 0.97, AUNU = 0.95, Macro precision = 0.9378, Min precision = 0.9524, Macro Sn = 0.9188, Micro Sn = 0.9524, Macro F1-score = 0.9269, Min F1-score = 0.9524 2. Nipponbare 5mC Acc = 0.8106, AUNP = 0.88, AUNU = 0.88, Macro precision = 80.74, Min precision = 81.93, Macro Sn = 81.93, Micro Sn = 81.93, Macro F1-score = 80.32, Min F1-score = 81.93 3. 6mA Chen Sn = 0.9303, Sp = 0.9225, Acc = 0.9265, MCC = 0.85, AUROC = 0.93 4. 6mA Lv Sn = 0.9605, Sp = 0.9248, Acc = 0.9426, MCC = 0.89, AUROC = 0.94
(154)Lv et al. datasets (5hmc: M. musculus, H. sapiens, 4mc: C. equisetifolia, F. vesca, S. cerevisiae, Tolypocladium, 6mA: A. thaliana, C. elegans, C. equisetifolia, D. melanogaster, F. vesca, H. sapiens, R. chinensis, S. cerevisiae, T. thermophile, Tolypocladium, XocBLS256)BERT_5hmC_H. sapiens: Acc = 94.92, Sn = 98.63, Sp = 91.21, MCC = 90.09, AUROC = 95.53, F1-score = 95.1 5hmC_M. musculus: Acc = 96.85, Sn = 97.06, Sp = 96.63, MCC = 93.69, AUROC = 97.57, F1-score = 96.85 4mC_C. equisetifolia: Acc = 82.51, Sn = 79.23, Sp = 85.79, MCC = 65.17, AUROC = 85.55, F1-score = 81.92 4mC_F. vesca: Acc = 84.2, Sn = 85.2, Sp = 83.21, MCC = 68.42, AUROC = 90.7, F1-score = 84.36 4mC_S. cerevisiae: Acc = 70.27, Sn = 66.94, Sp = 73.61, MCC = 40.64, AUROC = 75.37, F1-score = 69.25 4mC_Tolypocladium: Acc = 73.83, Sn = 72.16, Sp = 75.49, MCC = 47.68, AUROC = 80.57, F1-score = 73.39 6mA_A. thaliana: Acc = 85.38, Sn = 82.33, Sp = 88.42, MCC = 70.88, AUROC = 91.84, F1-score = 84.92 6mA_C. elegans: Acc = 89.03, Sn = 88.17, Sp = 89.9, MCC = 78.08, AUROC = 94.33, F1-score = 88.94 6mA_C. equisetifolia: Acc = 73.28, Sn = 68.91, Sp = 77.65, MCC = 46.73, AUROC = 79.02, F1-score = 72.06 6mA_D. melanogaster: Acc = 91.22, Sn = 90.38, Sp = 92.05, MCC = 82.44, AUROC = 95.44, F1-score = 91.14 6mA_F. vesca: Acc = 92.68, Sn = 92.33, Sp = 93.04, MCC = 82.44, AUROC = 95.44, F1-score = 92.66 6mA_H. sapiens: Acc = 89.8, Sn = 89.4, Sp = 90.2, MCC = 79.6, AUROC = 95.1, F1-score = 89.76 6mA_R. chinensis: Acc = 82.61, Sn = 80.94, Sp = 84.28, MCC = 65.25, AUROC = 87.89, F1-score = 82.31 6mA_S. cerevisiae: Acc = 80.11, Sn = 72. 37, Sp = 87.85, MCC = 60.96, AUROC = 87.09, F1-score = 78.44 6mA_T. thermophile: Acc = 87.4, Sn = 93.34, Sp = 81.54, MCC = 75.4, AUROC = 93.1, F1-score = 88.14 6mA_Tolypocladium: Acc = 77.38, Sn = 71.76, Sp = 83.01, MCC = 55.12, AUROC = 83.61, F1-score = 76.04 6mA_Xoc BLS256: Acc = 86.94, Sn = 88.9, Sp = 84.92, MCC = 73.94, AUROC = 92.61, F1-score = 87.2
Goal: environmental and microbial genomics tasks
Multi-class classificationNitrogen cycle prediction(27)NCycDBBERT_Macro F1-score = 99.5, Weighted F1-score = 99.2
Goal: gene network analysis
Multi-class classificationGene taxonomy classification(121)Verma et al. datasetFastTextMLPMacro F1-score = 0.92 ± 0.0054
(123)Mock et al. dataset (Superkingdom, Phylum)BERT_Superkingdom: Acc = 94.78 Phylum: Acc = 85.55
(122)ActinoMock dataset (LSH, Decimal, FNV)LSH+FastTextMLPLSH: Acc = 1.00, Precision = 1.00, Recall = 1.00, F1-score = 0.99 Decimal: Acc = 0.93, Precision = 0.94, Recall = 0.93, F1-score = 0.91 FNV: Acc = 0.937, Precision = 0.94, Recall = 0.94, F1-score = 0.92
Binary classificationGene network reconstruction(127)DREAM4 10, DREAM4 100, E. coli cold_LightGBMDREAM4 10: AUROC = 0.956, AUPRC = 0.891 DREAM4 100: AUROC = 0.909, AUPRC = 0.445 E. coli cold: AUROC = 0.602, AUPRC = 0.030
(126)Pio et al. datasetMetabolic feature encodingClustering_
(128)Ceci et al. datasetsNode2VecPCT_

Distinct predictive pipeline performance related to DNA modification, environmental and microbial genomics tasks, and gene network analysis.

Bold values highlight top performers on unique datasets related to distinct tasks.

Second most commonly used representation learning approach in this goal is Word2vec which is explored with two unique classifiers for two different tasks. Specifically, Word2vec with CNN classifier is used for one task, namely, 4mC-cytosine modification prediction, and with BiLSTM for one task, namely, 6mA-methyl adenine modification prediction. From most common approaches, BERT manages to achieve best performance with a self-classifier as compared to Word2vec-based predictive pipeline. Apart from BERT and Word2vec, Transformer with a self-classifier is used for three tasks, namely, 5mC-methycyctosine modification prediction, 6mA-methyl adenine modification prediction, and DNA methylation modification prediction. Transformer is also used with CNN classifier on one of the common task, namely, 6mA-methyl adenine modification prediction. Furthermore, potential of transformer-XL is explored with CNN classifier for one task, ULMFiT and hybrid encoding scheme (BERT+ALBERT+XLNet+ELECTRA) with self-classifier for one task, FastText with SVM classifier for one task and FastText with XGBoost classifier for one task. Overall, among all approaches, PseKNC encoding approach with structural sparse regularized random vector functional link network classifier manages to achieve best predictive performance on 4mC-methylcytosine modification prediction. Among all four tasks, DNA methylation modification prediction and 5mC-methyl cytosine modification prediction have some room for improvement. Building on the performance trends of predictive pipelines developed for different tasks of this goal, potential of BERT or PseKNC representation learning approach with structural sparse regularized random vector functional link network classifier can be explored to enhance the performance of under-performing tasks.

For environmental and microbial genomics goal, only potential of BERT representation learning is explored with a self-classifier. However, the potential of neural word embeddings and domain specific encoders based predictive pipelines remains unexplored.

For gene network analysis goal, across two different tasks, four unique representation learning approaches, namely, FastText, Node2vec, BERT, and metabolic encoding, along with four unique classifiers, namely, MLP, LightGBM, Clustering algorithm, and PCT, are used by six predictors. FastText representation learning is most commonly used among all approaches. Specifically, FastText along with MLP classifier is used for gene taxonomy classification task. Second most common representation learning approach is BERT that is used with a self-classifier for same gene taxonomy classification task. Apart from this, Node2Vec representation is explored with PCT classifier and metabolic feature encoding is explored with clustering algorithm for gene network reconstruction task. Among all approaches, FastText and MLP classifier-based predictive pipeline manages to achieve best performance for gene taxonomy classification task. Among all tasks of this goal, gene taxonomy classification offers some room for improvement. Building on promising performance achieved by contemporary language models for different sequence analysis tasks, hierarchical graph transformer and sophisticated machine or deep learning-based ensemble classifier can further raise the predictive performance on gene taxonomy classification task. Furthermore, advanced graph-based representation learning methods, GraRep, HOPE, and LINE with deep classifiers can also potentially raise the predictive performance on gene taxonomy classification task.

In addition, Table 11 provides an overview of six DNA sequence analysis tasks classified under the goal of DNA functional analysis. For this goal, four unique representation learning methods, namely, Word2Vec, Transformer, BERT, and FastText, are used in conjunction with three different predictors, namely, LogR, SVM, and cosine similarity. Among all four representation learning methods, Transformer is most commonly used followed by Word2vec and BERT. Transformer is used in three different tasks with self-classifier, Word2vec, and BERT are used with LogR and self-classifier in two different tasks. In addition, potential of FastText is explored with SVM for 1 task. Among all representation learning methods, Transformer with self-classifier manages to achieve top performance for tumor type prediction. Among all six tasks, disease risk estimation task offers a room for improvement as its respective BERT and self-classifier-based predictive pipeline accuracy falls approximately 56%. Hybrid approaches combining the powers of Transformer, BERT, and Word2vec with sophisticated machine learning classifier such as deep forest or deep learning classifiers such as CNN and CNN+BiGRU can potentially enhance the performance on under-performing tasks. Furthermore, except for two tasks, namely, species classification and functional prioritization, of non-coding variants, all other four tasks are evaluated on a single benchmark dataset. Considering deep learning models require huge amount of data to achieve promising performance, development, and utilization of more datasets in model building and validation can also prove fruitful for enhancing the predictive performance.

Table 11

Task typeTask nameReferencesDatasetRepresentation learningClassifierPerformance evaluation
Binary classificationConserved non-coding elements classification(163)Polychronopoulos et al. dataset (D1, D2, D3)Word2VecLogRD1: F1-score = 83.0, D2: F1-score = 85.5, D3: F1-score = 77.4
Multi-class classificationFunctional prioritization of non-coding variant(34)Yang et al. dataset (LOGO-919, LOGO-2002, LOGO-3357)Transformer__
Binary classificationExon and intron region classification(164)Akalin et al. datasetBERT_Precision = 100, Sn = 75, Sp = 100, Acc = 88.88, F1-score = 85.71
Binary classificationRecombination spots identification(165)Liu et al. datasetFastTextSVMSn = 90, Sp = 94.76, Acc = 92.6, MCC = 0.851
Multi-class classificationSpecies classification(44)1. Mouse enhancers dataset 2. Coding vs. intergenomic dataset 3. Human vs. worm dataset 4. Human enhancers cohn dataset 5. Human enhancers ensembl dataset 6. Human regulatory dataset 7. Human nontata promoter dataset 8. Human OCR ensembl datasetTransformer_Mouse enhancers: Acc = 85.1 Coding vs. Intergenomic: Acc = 91.3 human vs. worm: Acc = 96.6 human enhancers cohn: Acc = 74.2 human enhancers ensembl: Acc = 89.2 human regulatory: Acc = 93.8 human nontata promoter: Acc = 96.6 human OCR ensembl: Acc = 80.9
InteractionPrediction of context-specific functional impact of genetic variants(36)eQTLs datasetTransformer_AUPRC = 0.922

DNA functional analysis task predictive pipeline performance.

Bold values highlight top performers on unique datasets related to distinct tasks.

Finally, Table 12 summarizes predictive models developed for seven unique DNA sequence analysis tasks categorized under the goal of disease analysis. For this goal, seven unique representation learning methods, namely, Node2Vec, Graph2Vec, BERT, Graph Embedding, SDNE, Word2Vec, and Transformer, and five predictors, namely, MLP, RF, cosine similarity, clustering, and BERT self-classifier, are used in different tasks. Among all representation learning approaches, BERT and Word2vec are most commonly used. BERT is used with self-classifier on two different tasks, and word2vec is used with cosine similarity for a multi-class classification task, namely, mutation susceptibility analysis and with clustering algorithm for an only clustering task, namely, phylogenetic analysis. Apart from BERT and Word2vec, other representation learning methods Node2vec+Graph2vec, Graph Embedding, SDNE+Word2vec, and Transformer are used on one classification task each with MLP and self-classifiers. Overall, among all approaches, Transformer with self-classifier-based predictive pipelines manages to achieve best performance on tumor type prediction task. Among all seven tasks, disease risk estimation task offers a lot of room for improvement as its respective BERT with self-classifier-based predictive pipeline performance falls approximately 56%. Taking the transformer performance trends into account, latest sophisticated language models such as hierarchical graph transformer, ELECTRA, and GPT-4 along with ensemble machine or deep learning predictors can achieve significance performance rise in under-performing classification and clustering tasks.

Table 12

Task typeTask nameReferencesDatasetRepresentation learningClassifierPerformance evaluation
Binary classificationPathogen signatures identification(171)DS500 dataset, DS5000 datasetNode2Vec+Graph2VecMLP1. Acc = 73.49, 2. Acc = 89.7
Disease risks estimation(90)HSCR-RET, HSCR-RET-LongBERT_HSCR-RET: Precision = 0.770, Recall = 0.519, Acc = 0.562; HSCR-RET-Long: Precision = 0.768, Recall = 0.513, Acc = 0.541
InteractionPhage-host interactions prediction(175)Qiu et al. dataset_RFAcc = 0.801, Recall = 0.801, Sp = 0.801, Precision = 0.803, F1-score = 0.801, AUROC = 0.801
(176)Wang et al. datasetGraph EmbeddingMLPAUROC = 0.88317
(177)ESKAPE datasetSDNE+Word2VecMLPAcc = 86.65 ± 1.55, Sn = 88.40 ± 1.81, Sp = 84.91 ± 1.96, Precision = 85.43 ± 1.74, F1-score = 86.88 ± 1.53, AUROC = 0.9208 ± 0.0119
Multi-class classificationMutation susceptibility analysis(173)Yilmaz et al. dataset (Human, Mouse)Word2VecCosine similarityHuman data: Acc = 0.7974, mouse data: Acc = 0.8322
Tumor type prediction(180)TCGA pan-cancer datasetTransformer_Acc = 98.4, Precision = 98.50, Recall = 98.4, F1-score = 98.37
Pathogenicity potential assessment(27)E-K12, CARD-A, CARD-D, CARD-R, VFDB, ENZYME, PATRIC, NCycDBBERT_E-K12: Macro F1-score = 61.8, Weighted F1-score = 65.4; CARD-A AMR: Macro F1-score = 78.6, Weighted F1-score = 90.1; CARD-D: Macro F1-score = 57.4, Weighted F1-score = 85.2; CARD-R: Macro F1-score = 69.4, Weighted F1-score = 91.4; VFDB: Macro F1-score = 75.7, Weighted F1-score = 90.2; ENZYME: Macro F1-score = 99.1, Weighted F1-score = 98.8; PATRIC: Macro F1-score = 99.3, Weighted F1-score = 99.0; NCycDB: Macro F1-score = 99.5, Weighted F1-score = 99.2
ClusteringPhylogenetic analysis()Ren et al. datasetWord2VecClusteringAcc = 0.84

Disease analysis task predictive pipeline performance.

Bold values highlight top performers on unique datasets related to distinct tasks.

In a nutshell, a comprehensive analysis of state-of-the-art predictive pipelines developed using word embeddings, language models, and nucleotide compositional and positional information-based encoders reveals interesting trends. From 44 DNA sequence analysis tasks classified under the hood of 8 major biological goals, 24 tasks belong to binary classification, 4 belong to interaction prediction, 11 belong to multi-class classification, only 3 belong to multi-label classification, 1 belong to regression, and 1 belong to clustering. Overall, 25 unique representation learning methods and 28 predictors are explored for developing robust predictive pipelines for 44 DNA sequence analysis tasks classified under the hood of 8 major biological goals. Across all eight goals, language model-based representation learning approaches and deep learning classifiers are achieving better performance across majority of the tasks. Researchers can explore the performance potential of latest transformer-based language models such as Hierarchical graph transformer, GPT-4, and hybrid representation learning methods along with sophisticated ensemble machine learning or deep learning predictors for different classification, regression, and clustering tasks.

11 Publisher and journal-wise distribution of research articles

This section provides an overview of 44 distinct DNA sequence analysis task-related articles distribution across conferences, journals, and publishers. Before paper submission, identification of relevant journals for a study publication in the interdisciplinary field of AI applications in DNA sequence analysis is an important task. There are three types of journals in AI and DNA sequence analysis fields: (1) Journals focusing on core AI algorithms, (2) Journals dedicated to core biological findings, (3) Hybrid journals that publish research integrating both AI algorithms and biological data. Researchers often face desk rejections when submitting to core AI or biology journals. Instead, they should target hybrid journals. While many tools exist to find suitable journals, this comprehensive guide provides detailed information to help researchers to identify journals where applications using word embeddings and large language models for DNA sequence analysis are published.

Figure 7 graphically depicts distribution of 127 studies across 53 journals, 1 transactions, 3 conferences, and 2 pre-print repositories. Among all journals, more studies are published in Briefings in Bioinformatics followed by Bioinformatics, Computational Biology, and Chemistry, and International Journal of Molecular Sciences. Similarly, among all conferences, more studies are published in the International Conference on Bioinformatics and Biomedicine (BIBM) followed by the 11th Hellenic Conference on Artificial Intelligence, Proceedings of the 12th and 13th ACM International Conference on Bioinformatics, Computational Biology, and Health Informatics. Moreover, 5 studies are published in ACM transaction of computational biology. In the light of rapid development in research findings, researchers have also published 24 studies in bioRxiv, ArXiV, and MedXiv platforms. However, researchers generally prefer journal publications for their sustained impact.

Figure 7

Furthermore, Figure 8 illustrate that 127 DNA sequence analysis studies are published by 17 different publishers, namely, Springer,9 Elsevier (see text footnote 3), Oxford University Press,10 Cold Spring Harbor Laboratory,11 IEEE,12 Ozer UYGUN,13 ACS Publication,14 Frontier Media SA,15 Gazi University,16 Marry Ann Liebert,17 MDPI,18 National Acad Sciences,19 Nature Publishing Group UK London,20 PeerJ Inc.,21 Public Library of science,22 ACM (see text footnote 2), and pre-prints.23 Notably, approximately 60 out of 127 DNA sequence analysis studies are published by Oxford University Press, Elsevier, and Cold Spring Harbor Laboratory. In addition, IEEE, Springer, and MDPI have contributed 30 relevant papers. Furthermore, 32 DNA sequence analysis research articles are published by ACS Publications, Frontiers Media SA, Mary Ann Liebert, Inc., National Acad Sciences, Nature Publishing Group UK London, Public Library of Science, PeerJ Inc, and others. Collectively, 96 are journal publications, 6 are conference papers, 1 is transaction articles, and 24 are pre-prints out of 127 DNA sequence analysis studies published by 21 different publishers. This comprehensive analysis across various journals, conferences, transactions, and pre-print repositories highlights diverse and extensive research landscape in DNA sequence analysis.

Figure 8

12 Discussion

We acknowledge that “DNA sequence analysis” encompasses a much broader range of bioinformatics applications than covered in this review. Our focus is specifically on AI-based approaches that analyze raw DNA sequence data to predict biological functions and features. Other important areas of bioinformatics such as genome assembly, comprehensive variant analysis, phylogenomics, and many aspects of population genetics utilize different computational approaches and would benefit from separate dedicated reviews. A comprehensive review of existing literature on AI-driven DNA sequence analysis tasks reveals a significant inconsistency in the evaluation of predictive pipelines across similar datasets. Researchers have developed numerous datasets tailored to specific tasks, and most of the researchers have evaluated their proposed predictors solely on their own datasets.

The creation of new datasets is essential because public databases are frequently updated with new sequence information. These new datasets can incorporate the most recent sequence data alongside existing information. Moreover, existing datasets tend to be smaller, while deep learning models perform better with larger datasets. To address performance comparison inconsistency, there is an urgent need to standardize dataset utilization. One potential solution is to benchmark existing predictors on newly developed datasets and compare the performance of proposed predictors against these benchmarks. This approach would provide a more objective evaluation of proposed predictor performance. However, in this context, a significant challenge is the limited availability of source codes for existing predictors. Many studies make their source code private, hindering the reproducibility of results and direct comparison with other methods. To streamline the integration of innovative methods and ensure methodological advancement, it is crucial to analyze task-specific datasets and create standardized datasets with detailed descriptions. By benchmarking the performance of existing predictors on these standardized datasets, researchers can establish a common ground for comparison and facilitate more accurate evaluations of models. This approach would enhance transparency and reproducibility in DNA sequence analysis studies.

The development of AI-driven predictive pipelines for DNA sequence analysis relies heavily on effective sequence representation learning methods and appropriate machine or deep learning models. Machine and deep learning models inherently depend on statistical vectors and cannot process raw DNA sequences directly. Therefore, the role of representation learning methods in these pipelines is crucial. These methods are responsible for transforming raw DNA sequences into statistical vectors by capturing and encoding the most informative nucleotide patterns.

In the current landscape of AI-driven DNA sequence analysis, researchers have employed a variety of representation learning methods, including 12 distinct word embedding techniques and 8 language models. However, when it comes to other genetic molecules such as RNA and proteins, researchers have explored an additional set of 17 word embedding methods and 13 language models that have not yet been applied to DNA sequence analysis. These unexplored word embedding methods include DANE (285), ELMo (286288), GATNE (289), GEMSEC (290), MetaGraph2Vec (291), HAKE (292), HIN2Vec (293), HOPE (294, 295), LINE (296298), Mashup (299, 300), Random Watcher-Walker (RW2) (301), RotatE (292, 302, 303), RWR (304), Struc2Vec (305, 306), SVD (307, 308), Topo2Vec (309), and TransE (310), while the unexplored language models include AlphaFold (311315), AlphaFold2 (316, 317), BigBird (318), ESM-1 (315, 319, 320), ESM-2 (109, 286, 316, 320), Graph Transformer Network (321), Heterogeneous Graph Transformer (322), IgFold (323), LongFormer (318), RoBERTa (324, 325), T5 (320, 326328), and Vision Transformer (288). Integrating these advanced word embedding techniques and large language models into AI-driven DNA sequence analysis pipelines could potentially enhance their performance and robustness.

Within 127 AI-driven DNA sequence analysis predictive pipelines, researchers have utilized 18 machine and deep learning algorithms at the predictor level. In some cases, they have developed meta-predictors by combining multiple machine learning and deep learning algorithms to enhance predictive performance. However, similar to the representation learning stage, there are 24 distinct methods at the predictor level that have not yet been explored, representing untapped potential for improving the accuracy and robustness of these AI-driven pipelines.

Our categorization of 44 distinct tasks into 8 biological goals provides a structured framework that serves as a valuable starting taxonomy for both computer scientists and life scientists. This organization is informed by both computational and biological literature and creates a common reference point that bridges these disciplines while facilitating interdisciplinary communication. We recognize the inherent complexity of biological systems and the interconnected nature of many of these tasks. For example, enhancer identification categorized under gene expression regulation shares biological connections with chromatin accessibility prediction categorized under genome structure and stability. Nevertheless, this framework offers a practical organizing principle that will naturally evolve and be refined over time. The taxonomy presented here lays groundwork that future collaborative efforts between AI researchers and domain specialists in genomics can build upon. We anticipate gradual development into a more nuanced framework that maintains practical utility while better reflecting biological realities. Similar to many scientific classification systems, we expect this taxonomy to mature through iterative refinement as the field advances.

13 Conclusion

This review serves as a comprehensive resource for researchers working at the intersection of AI and DNA sequence analysis. It provides a structured foundation for future innovations in the rapidly evolving field of computational genomics. It bridges the critical gap between molecular biology and artificial intelligence by systematically analyzing 44 different DNA sequence analysis tasks, their associated databases, datasets, and AI methodologies. It identifies 36 biological databases and 140 benchmark datasets that provide a robust foundation for developing and evaluating AI predictors. Furthermore, our examination of existing predictive pipelines demonstrates the successful application of 39 word embeddings and 67 language models across various DNA sequence analysis tasks. Our analysis reveals that while significant progress has been made in developing AI-driven predictive pipelines for DNA sequence analysis, several challenges and opportunities remain unexplored. Several promising directions emerge for the advancement of this field. First, the integration of 17 unexplored word embedding methods and 13 language models (currently utilized only for RNA and protein analysis) could significantly enhance DNA sequence analysis capabilities. Second, the development of standardized benchmark datasets and evaluation protocols would facilitate fair comparisons between different predictive models and accelerate progress in the field. Third, the adoption of 24 untapped machine learning and deep learning algorithms at the predictor level presents an opportunity to improve prediction accuracy and robustness.

Future research should focus on developing multi-task learning frameworks that can simultaneously handle multiple DNA sequence analysis tasks, thereby improving computational efficiency and leveraging shared biological features. Furthermore, ensuring public accessibility of source codes and detailed documentation of predictive pipelines would foster reproducibility and collaborative advancement in the field. The establishment of standardized performance metrics and evaluation protocols across different DNA sequence analysis tasks would enable more meaningful comparisons between various approaches and guide future developments. As DNA sequence data continue to grow exponentially, the integration of more sophisticated AI architectures, particularly those capable of handling large-scale genomic data efficiently, will become increasingly important. Our categorization of 44 tasks reflects common AI applications in DNA sequence analysis literature and provides a starting point for interdisciplinary discourse. We recognize that deeper collaboration between AI researchers and life scientists would further strengthen the biological relevance of this framework. A greater amount of input from geneticists and bioinformaticians would be essential to develop a more comprehensive and biologically relevant tasks taxonomy.

Statements

Author contributions

MA: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. MI: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing. AZ: Data curation, Formal analysis, Writing – original draft. AD: Supervision, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

MA, MI, AZ, and AD were employed by German Research Center for Artificial Intelligence GmbH. MA and AD are shareholders of company Intelligentx GmbH.

Generative AI statement

The author(s) declare Gen AI was used in the creation of this manuscript. During the preparation of this study, the authors used Grammarly tool to fix language and grammar issues and ChatGpt for outlining, better understanding different studies, and expansion of concepts. After using these tools, the authors reviewed and edited the content as needed and take full responsibility for the content of the publication.

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.

    AllisonLA. Fundamental Molecular Biology. New York: John Wiley & Sons. (2021).

  • 2.

    GartnerAEngebrechtJ. DNA repair, recombination, and damage signaling. Genetics. (2022) 220:iyab178. 10.1093/genetics/iyab178

  • 3.

    ChenJPotlapalliRQuanHChenLXieYPouriyehSet al. Exploring DNA damage and repair mechanisms: a review with computational insights. BioTech. (2024) 13:3. 10.3390/biotech13010003

  • 4.

    AseervathamGSBSivasudhaTJeyadeviRArul AnanthD. Environmental factors and unhealthy lifestyle influence oxidative stress in humans–an overview. Environ Sci Pollut Res. (2013) 20:435669. 10.1007/s11356-013-1748-0

  • 5.

    LiaoXZhuWZhouJLiHXuXZhangBet al. Repetitive DNA sequence detection and its role in the human genome. Commun Biol. (2023) 6:954. 10.1038/s42003-023-05322-y

  • 6.

    JoiretMLeclercqMLambrechtsGRapinoFClosePLouppeGet al. Cracking the genetic code with neural networks. Front Artif Intell. (2023) 6:1128153. 10.3389/frai.2023.1128153

  • 7.

    LaubVDevrajKEliasLSchulteD. Bioinformatics for wet-lab scientists: practical application in sequencing analysis. BMC Genomics. (2023) 24:382. 10.1186/s12864-023-09454-7

  • 8.

    ElshafeiAAl-ToubatMFeibusAHKoulKJazayeriSBLelaniNet al. Genetic mutations in smoking-associated prostate cancer. Prostate. (2023) 83:122937. 10.1002/pros.24554

  • 9.

    GaoYLiuM. Application of machine learning based genome sequence analysis in pathogen identification. Front Microbiol. (2024) 15:1474078. 10.3389/fmicb.2024.1474078

  • 10.

    BiZDipSAHajialigolDKommuSLiuHLuMet al. AI for biomedicine in the era of large language models. arXiv preprint arXiv:240315673. (2024).

  • 11.

    ZhangQDingKLyvTWangXYinQZhangYet al. Scientific large language models: a survey on biological & chemical domains. arXiv preprint arXiv:240114656. (2024).

  • 12.

    AsimMN. An efficient automated machine learning framework for genomics and proteomics sequence analysis. Rheinland-Pfälzische Technische Universität Kaiserslautern-Landau. (2023).

  • 13.

    O'LearyNACoxEHolmesJBAndersonWRFalkRHemVet al. Exploring and retrieving sequence and metadata for species across the tree of life with NCBI datasets. Sci Data. (2024) 11:732. 10.1038/s41597-024-03571-y

  • 14.

    AbbasiAFAsimMNAhmedSDengelA. Long extrachromosomal circular DNA identification by fusing sequence-derived features of physicochemical properties and nucleotide distribution patterns. Sci Rep. (2024) 14:9466. 10.1038/s41598-024-57457-5

  • 15.

    JinYTTanYGanZHHaoYDWangTYLinHet al. Identification of DNase I hypersensitive sites in the human genome by multiple sequence descriptors. Methods. (2024) 229:12532. 10.1016/j.ymeth.2024.06.012

  • 16.

    LiYWeiXYangQXiongALiXZouQet al. msBERT-Promoter: a multi-scale ensemble predictor based on BERT pre-trained model for the two-stage prediction of DNA promoters and their strengths. BMC Biol. (2024) 22:126. 10.1186/s12915-024-01923-z

  • 17.

    YangBGuneriDYuHWrightEPChenWWallerZAet al. Prediction of DNA i-motifs via machine learning. Nucleic Acids Res. (2024) 52:218897. 10.1093/nar/gkae092

  • 18.

    WangXQiaoLQuPYangQ. TBCA: prediction of transcription factor binding sites using a deep neural network with lightweight attention mechanism. IEEE J Biomed Health Inf. (2024) 28:23972407. 10.1109/JBHI.2024.3355758

  • 19.

    WangXDuQWangR. Mus4mCPred: accurate identification of DNA N4-methylcytosine sites in mouse genome using multi-view feature learning and deep hybrid network. Processes. (2024) 12:1129. 10.3390/pr12061129

  • 20.

    ChenZZhaoPLiCLiFXiangDChenYZet al. iLearnPlus: a comprehensive and automated machine-learning platform for nucleic acid and protein sequence analysis, prediction and visualization. Nucleic Acids Res. (2021) 49:e60. 10.1093/nar/gkab122

  • 21.

    RenRYinCS-T YauS. kmer2vec: a novel method for comparing DNA sequences by word2vec embedding. J Comput Biol. (2022) 29:100121. 10.1089/cmb.2021.0536

  • 22.

    ZhouZWuWHoHWangJShiLDavuluriRVet al. Dnabert-s: Learning species-aware DNA embedding with genome foundation models. arXiv:2402.08777. (2024).

  • 23.

    LiHLPangYHLiuB. BioSeq-BLM: a platform for analyzing DNA, RNA and protein sequences based on biological language models. Nucleic Acids Res. (2021) 49:e129e129. 10.1093/nar/gkab829

  • 24.

    AirlanggaG. Comparative analysis of deep learning architectures for DNA sequence classification: performance evaluation and model insights. J Comput Syst Inf. (2024) 5:70918.

  • 25.

    WangCHeZJiaRPanSCoinLJSongJet al. PLANNER: a multi-scale deep language model for the origins of replication site prediction. IEEE J Biomed Health Inf. (2024) 28:24452454. 10.1109/JBHI.2024.3349584

  • 26.

    YinZNLaiFLGaoF. Unveiling human origins of replication using deep learning: accurate prediction and comprehensive analysis. Brief Bioinf. (2024) 25:bbad432. 10.1093/bib/bbad432

  • 27.

    DuanCZangZXuYHeHLiuZSongZet al. FGBERT: function-driven pre-trained gene language model for metagenomics. arXiv preprint arXiv:240216901. (2024).

  • 28.

    LiJChiuTPRohsR. Predicting DNA structure using a deep learning method. Nat Commun. (2024) 15:1243. 10.1038/s41467-024-45191-5

  • 29.

    FazeelAAghaADengelAAhmedS. NP-BERT: a two-staged BERT based nucleosome positioning prediction architecture for multiple species. In: BIOINFORMATICS. (2023). p. 175187. 10.5220/0011679200003414

  • 30.

    Han GS LiQLiY. Comparative analysis and prediction of nucleosome positioning using integrative feature representation and machine learning algorithms. BMC Bioinform. (2021) 22:129. 10.1186/s12859-021-04006-w

  • 31.

    ZhangYChuXJiangYWuHQuanL. SemanticCAP: chromatin accessibility prediction enhanced by features learning from a language model. Genes. (2022) 13:568. 10.3390/genes13040568

  • 32.

    FishmanVKuratovYPetrovMShmelevAShepelinDChekanovNet al. GENA-LM: a family of open-source foundational DNA language models for long sequences. bioRxiv. (2023). p. 202306. 10.1101/2023.06.12.544594

  • 33.

    GuoYZhouDNieRRuanXLiW. DeepANF: A deep attentive neural framework with distributed representation for chromatin accessibility prediction. Neurocomputing. (2020) 379:30518. 10.1016/j.neucom.2019.10.091

  • 34.

    YangMHuangHHuangLZhangNWuJYangHet al. LOGO, a contextualized pre-trained language model of human genome flexibly adapts to various downstream tasks by fine-tuning. (2021). 10.21203/rs.3.rs-448927/v1

  • 35.

    TaoHLiHXuKHongHJiangSDuGet al. Computational methods for the prediction of chromatin interaction and organization using sequence and epigenomic profiles. Brief Bioinform. (2021) 22:bbaa405. 10.1093/bib/bbaa405

  • 36.

    GaoZLiuQZengWJiangRWongWH. EpiGePT: a pretrained transformer model for epigenomics. bioRxiv. (2023). 10.1101/2023.07.15.549134

  • 37.

    VanhaerenTDivinaFGarcía-TorresMGómez-VelaFVanhoofWMartínez-GarcíaPM. A comparative study of supervised machine learning algorithms for the prediction of long-range chromatin interactions. Genes. (2020) 11:985. 10.3390/genes11090985

  • 38.

    ZhangZLiFZhaoJZhengC. CapsNetYY1: identifying YY1-mediated chromatin loops based on a capsule network architecture. BMC Genomics. (2023) 24:448. 10.1186/s12864-023-09217-4

  • 39.

    AbbasiAFAsimMNTryggJDengelAAhmedS. Deep learning architectures for the prediction of YY1-mediated chromatin loops. In: International Symposium on Bioinformatics Research and Applications. Springer (2023). p. 7284. 10.1007/978-981-99-7074-2_6

  • 40.

    ChenWFengPLinH. Prediction of replication origins by calculating DNA structural properties. FEBS Lett. (2012) 586:9348. 10.1016/j.febslet.2012.02.034

  • 41.

    NieduszynskiCHiragaSAkPBenhamCDonaldsonA. Oridb: a DNA replication origin database. Nucleic Acids Res. (2006) 35:D406. 10.1093/nar/gkl758

  • 42.

    JordanKHeFSotoMAkhunovaAAkhunovE. Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes. Genome Biol. (2020) 21:176. 10.1186/s13059-020-02093-1

  • 43.

    MinnoyeLMarinovGKrausgruberTPanLMarandASecchiaSet al. Chromatin accessibility profiling methods. Nat Rev Methods Primers. (2021) 1:10. 10.1038/s43586-020-00008-9

  • 44.

    NguyenEPoliMFaiziMThomasAWornowMBirch-SykesCet al. HyenaDNA: long-range genomic sequence modeling at single nucleotide resolution. In: Advances in Neural Information Processing Systems. (2024). p. 36.

  • 45.

    DaoFYLvHZhangDZhangZMLiuLLinH. DeepYY1: a deep learning approach to identify YY1-mediated chromatin loops. Brief Bioinform. (2021) 22:bbaa356. 10.1093/bib/bbaa356

  • 46.

    LiZGaoEZhouJHanWXuXGaoX. Applications of deep learning in understanding gene regulation. Cell Rep Methods. (2023) 3:100384. 10.1016/j.crmeth.2022.100384

  • 47.

    KhanalJTayaraHChongKT. Identifying enhancers and their strength by the integration of word embedding and convolution neural network. IEEE Access. (2020) 8:5836976. 10.1109/ACCESS.2020.2982666

  • 48.

    LiQZhangLXuLZouQWuJLiQ. Identification and classification of promoters using the attention mechanism based on long short-term memory. Front Comput Sci. (2022) 16:164348. 10.1007/s11704-021-0548-9

  • 49.

    MinXYeCLiuXZengX. Predicting enhancer-promoter interactions by deep learning and matching heuristic. Brief Bioinform. (2021) 22:bbaa254. 10.1093/bib/bbaa254

  • 50.

    ClauwaertJWaegemanW. Novel transformer networks for improved sequence labeling in genomics. IEEE/ACM Trans Computat Biol Bioinform. (2020) 19:97106. 10.1109/TCBB.2020.3035021

  • 51.

    GhoshNSantoniDSahaIFeliciG. Predicting transcription factor binding sites with deep learning. Int J Mol Sci. (2024) 25:4990. 10.3390/ijms25094990

  • 52.

    XuYWangCXuKDingYLyuAZhangL. TRAFICA: improving transcription factor binding affinity prediction using deep language model on ATAC-seq data. bioRxiv. (2023). 10.1101/2023.11.02.565416

  • 53.

    LiuYTianB. Protein-DNA binding sites prediction based on pre-trained protein language model and contrastive learning. Brief Bioinform. (2024) 25:bbad488. 10.1093/bib/bbad488

  • 54.

    ZhangSHuHJiangTZhangLZengJ. TITER predicting translation initiation sites by deep learning. Bioinformatics. (2017) 33:i23442. 10.1093/bioinformatics/btx247

  • 55.

    XuJSunWLiKZhangWZhangWZengYet al. MNESEDA: a prior-guided subgraph representation learning framework for predicting disease-related enhancers. Knowl-Based Syst. (2024) 294:111734. 10.1016/j.knosys.2024.111734

  • 56.

    WangY. EnhancerBD identifing sequence feature. bioRxiv. (2024). 10.1101/2024.03.05.583459

  • 57.

    MehmoodFArshadSShoaibM. ADH-Enhancer: an attention-based deep hybrid framework for enhancer identification and strength prediction. Brief Bioinform. (2024) 25:bbae030. 10.1093/bib/bbae030

  • 58.

    LiaoMZhaoJPTianJZhengCH. iEnhancer-DCLA: using the original sequence to identify enhancers and their strength based on a deep learning framework. BMC Bioinform. (2022) 23:480. 10.1186/s12859-022-05033-x

  • 59.

    InayatNKhanMIqbalNKhanSRazaMKhanDMet al. iEnhancer-DHF: identification of enhancers and their strengths using optimize deep neural network with multiple features extraction methods. IEEE Access. (2021) 9:4078396. 10.1109/ACCESS.2021.3062291

  • 60.

    LeNQKHoQTNguyenTTDOuYY. A transformer architecture based on BERT and 2D convolutional neural network to identify DNA enhancers from sequence information. Brief Bioinform. (2021) 22:bbab005. 10.1093/bib/bbab005

  • 61.

    WangYHouZYangYWongKCLiX. Genome-wide identification and characterization of DNA enhancers with a stacked multivariate fusion framework. PLoS Comput Biol. (2022) 18:e1010779. 10.1371/journal.pcbi.1010779

  • 62.

    LeeDKarchinRBeerMA. Discriminative prediction of mammalian enhancers from DNA sequence. Genome Res. (2011) 21:216780. 10.1101/gr.121905.111

  • 63.

    MinXZengWChenSChenNChenTJiangR. Predicting enhancers with deep convolutional neural networks. BMC Bioinform. (2017) 18:478. 10.1186/s12859-017-1878-3

  • 64.

    LeNQKYappEKYHoQTNagasundaramNOuYYYehHY. iEnhancer-5Step: identifying enhancers using hidden information of DNA sequences via Chou's 5-step rule and word embedding. Anal Biochem. (2019) 571:5361. 10.1016/j.ab.2019.02.017

  • 65.

    GengQYangRZhangLA. deep learning framework for enhancer prediction using word embedding and sequence generation. Biophys Chem. (2022) 286:106822. 10.1016/j.bpc.2022.106822

  • 66.

    YangRWuFZhangCZhangL. iEnhancer-GAN: a deep learning framework in combination with word embedding and sequence generative adversarial net to identify enhancers and their strength. Int J Mol Sci. (2021) 22:3589. 10.3390/ijms22073589

  • 67.

    LiJWuZLinWLuoJZhangJChenQet al. iEnhancer-ELM: improve enhancer identification by extracting position-related multiscale contextual information based on enhancer language models. Bioinform Adv. (2023) 3:vbad043. 10.1093/bioadv/vbad043

  • 68.

    OmarNWong YS LiXChongYLAbdullahMTLeeNK. Enhancer prediction in proboscis monkey genome: a comparative study. J Telecommun Electr Comput Eng. (2017) 9:1759.

  • 69.

    LiuBFangLLongRLanXChouKC. iEnhancer-2L: A Two-Layer Predictor for Identifying Enhancers and Their Strength by Pseudo k-Tuple Nucleotide Composition. Oxford: Oxford University Press. (2016). 10.1093/bioinformatics/btv604

  • 70.

    JiaCHeW. EnhancerPred: a predictor for discovering enhancers based on the combination and selection of multiple features. Sci Rep. (2016) 6:38741. 10.1038/srep38741

  • 71.

    LiuBLiKHuangDSChouKC. iEnhancer-EL: identifying enhancers and their strength with ensemble learning approach. Bioinformatics. (2018) 34:383542. 10.1093/bioinformatics/bty458

  • 72.

    LiuB. iEnhancer-PsedeKNC: Identification of enhancers and @articlebgroups based on Pseudo degenerate kMER nucleotide composition. Neurocomputing. (2016) 217:4652. 10.1016/j.neucom.2015.12.138

  • 73.

    AsimMNIbrahimMAMalikMIDengelAAhmedS. Enhancer-DSNet: a supervisedly prepared enriched sequence representation for the identification of enhancers and their strength. In: International Conference on Neural Information Processing. Springer (2020). p. 3848. 10.1007/978-3-030-63836-8_4

  • 74.

    HeWJiaC. EnhancerPred2.0: predicting enhancers and their strength based on position-specific trinucleotide propensity and electron-ion interaction potential feature selection. Molec Biosyst. (2017) 13:76774. 10.1039/C7MB00054E

  • 75.

    LiuQChenPWangBZhangJLiJ. Hot spot prediction in protein-protein interactions by an ensemble system. BMC Syst Biol. (2018) 12:8999. 10.1186/s12918-018-0665-8

  • 76.

    ZhangPZhangHWuH. Ipro-wael: a comprehensive and robust framework for identifying promoters in multiple species. Nucleic Acids Res. (2022) 50:1027889. 10.1093/nar/gkac824

  • 77.

    YangGLiJHuJShiJY. Recognition of cyanobacteria promoters via Siamese network-based contrastive learning under novel non-promoter generation. Brief Bioinform. (2024) 25:bbae193. 10.1093/bib/bbae193

  • 78.

    LeNQKYappEKYNagasundaramNYehHY. Classifying promoters by interpreting the hidden information of DNA sequences via deep learning and combination of continuous fasttext N-grams. Front Bioeng Biotechnol. (2019) 7:305. 10.3389/fbioe.2019.00305

  • 79.

    TahirMHayatMGulSChongKT. An intelligent computational model for prediction of promoters and their strength via natural language processing. Chemometr Intell Lab Syst. (2020) 202:104034. 10.1016/j.chemolab.2020.104034

  • 80.

    LiFChenJGeZWenYYueYHayashidaMet al. Computational prediction and interpretation of both general and specific types of promoters in Escherichia coli by exploiting a stacked ensemble-learning framework. Brief Bioinform. (2021) 22:212640. 10.1093/bib/bbaa049

  • 81.

    LiuBYangFHuangDSChouKC. iPromoter-2L: a two-layer predictor for identifying promoters and their types by multi-window-based PseKNC. Bioinformatics. (2018) 34:3340. 10.1093/bioinformatics/btx579

  • 82.

    HongZZengXWeiLLiuX. Identifying enhancer-promoter interactions with neural network based on pre-trained DNA vectors and attention mechanism. Bioinformatics. (2020) 36:103743. 10.1093/bioinformatics/btz694

  • 83.

    NiYFanLWangMZhangNZuoYLiaoM. EPI-mind: identifying enhancer-promoter interactions based on transformer mechanism. Interdiscipl Sci. (2022) 14:78694. 10.1007/s12539-022-00525-z

  • 84.

    TalukderASaadatSLiXHuH. EPIP a novel approach for condition-specific enhancer-promoter interaction prediction. Bioinformatics. (2019) 35:387783. 10.1093/bioinformatics/btz641

  • 85.

    ZhuangZShenXPanW. A simple convolutional neural network for prediction of enhancer-promoter interactions with DNA sequence data. Bioinformatics. (2019) 35:2899906. 10.1093/bioinformatics/bty1050

  • 86.

    SinghSYangYPóczosBMaJ. Predicting enhancer-promoter interaction from genomic sequence with deep neural networks. Quant Biol. (2019) 7:12237. 10.1007/s40484-019-0154-0

  • 87.

    LiuSXuXYangZZhaoXLiuSZhangW. EPIHC: improving enhancer-promoter interaction prediction by using hybrid features and communicative learning. In: IEEE/ACM Transactions on Computational Biology and Bioinformatics. (2021). p. 11.

  • 88.

    ZhangMHuYZhuM. Epishilbert: prediction of enhancer-promoter interactions via hilbert curve encoding and transfer learning. Genes. (2021) 12:1385. 10.3390/genes12091385

  • 89.

    AnWGuoYBianYMaHYangJLiCet al. MoDNA: motif-oriented pre-training for DNA language model. In: Proceedings of the 13th ACM International Conference on Bioinformatics, Computational Biology and Health Informatics. (2022). p. 15. 10.1145/3535508.3545512

  • 90.

    MoSFuXHongCChenYZhengYTangXet al. Multi-modal self-supervised pre-training for regulatory genome across cell types. arXiv preprint arXiv:211005231. (2021).

  • 91.

    ClauwaertJMenschaertGWaegemanW. Explainability in transformer models for functional genomics. Brief Bioinform. (2021) 22:bbab060. 10.1093/bib/bbab060

  • 92.

    ShenZBaoWHuangDS. Recurrent neural network for predicting transcription factor binding sites. Sci Rep. (2018) 8:15270. 10.1038/s41598-018-33321-1

  • 93.

    JiYZhouZLiuHDavuluriRV. DNABERT pre-trained bidirectional encoder representations from transformers model for DNA-language in genome. Bioinformatics. (2021) 37:211220. 10.1093/bioinformatics/btab083

  • 94.

    KabirABhattaraiMRasmussenShehuABishopARAlexandrovBet al. Advancing transcription factor binding site prediction using DNA breathing dynamics and sequence transformers via cross attention. bioRxiv. (2024). p. 202401. 10.1101/2024.01.16.575935

  • 95.

    MuradTAliSChourasiaPPattersonM. Advancing protein-DNA binding site prediction: integrating sequence models and machine learning classifiers. bioRxiv. (2023). p. 202308. 10.1101/2023.08.23.554389

  • 96.

    LuoHShanWChenCDingPLuoL. Improving language model of human genome for DNA-protein binding prediction based on task-specific pre-training. Interdiscipl Sci. (2023) 15:3243. 10.1007/s12539-022-00537-9

  • 97.

    ShaMRahamathullaMP. Splice site recognition-deciphering Exon-Intron transitions for genetic insights using Enhanced integrated Block-Level gated LSTM model. Gene. (2024) 915:148429. 10.1016/j.gene.2024.148429

  • 98.

    KabangaEYunSVan MessemADe NeveW. Impact of U2-type introns on splice site prediction in Arabidopsis thaliana using deep learning. bioRxiv. (2024). p. 202405. 10.1101/2024.05.13.593811

  • 99.

    LiuXZhangHZengYZhuXZhuLFuJ. DRANetSplicer: a splice site prediction model based on deep residual attention networks. Genes. (2024) 15:404. 10.3390/genes15040404

  • 100.

    Dalla-TorreHGonzalezLMendoza-RevillaJCarranzaNLGrzywaczewskiAHOteriFet al. The nucleotide transformer: Building and evaluating robust foundation models for human genomics. bioRxiv. (2023). p. 202301. 10.1101/2023.01.11.523679

  • 101.

    LiWGuoYWangBYangB. Learning spatiotemporal embedding with gated convolutional recurrent networks for translation initiation site prediction. Pattern Recognit. (2023) 136:109234. 10.1016/j.patcog.2022.109234

  • 102.

    ReddyAJHerschlMHKolliSLuAXGengXKumarAet al. Strategies for effectively modelling promoter-driven gene expression using transfer learning. bioRxiv. (2023). 10.1101/2023.02.24.529941

  • 103.

    Al TaweraqiNKingRD. Improved prediction of gene expression through integrating cell signalling models with machine learning. BMC Bioinform. (2022) 23:323. 10.1186/s12859-022-04787-8

  • 104.

    SharmaKMarucciLAbdallahZS. FluxGAT: integrating flux sampling with graph neural networks for unbiased gene essentiality classification. arXiv preprint arXiv:240318666. (2024).

  • 105.

    LeNQKDoDTHungTNKLamLHTHuynhTTNguyenNTKet al. computational framework based on ensemble deep neural networks for essential genes identification. Int J Mol Sci. (2020) 21:9070. 10.3390/ijms21239070

  • 106.

    ZhangXXiaoWXiaoW. DeepHE: Accurately predicting human essential genes based on deep learning. PLoS Comput Biol. (2020) 16:e1008229. 10.1371/journal.pcbi.1008229

  • 107.

    XiaoWZhangXXiaoW. A Deep Learning Framework for Predicting Human Essential Genes by Integrating Sequence and Functional data. bioRxiv. (2020). p. 202008.

  • 108.

    SchapkeJTavaresARecamonde-MendozaM. EPGAT: gene essentiality prediction with graph attention networks. IEEE/ACM Trans Comput Biol Bioinform. (2021) 19:161526. 10.1109/TCBB.2021.3054738

  • 109.

    MaJSongJYoungNDChangBCKorhonenPKCamposTLet al. ‘Bingo'–a large language model-and graph neural network-based workflow for the prediction of essential genes from protein data. Brief Bioinform. (2024) 25:bbad472. 10.1093/bib/bbad472

  • 110.

    NunesSSousaRTPesquitaC. Multi-domain knowledge graph embeddings for gene-disease association prediction. J Biomed Semantics. (2023) 14:11. 10.1186/s13326-023-00291-x

  • 111.

    FanKZhangY. Pseudo2GO: a graph-based deep learning method for pseudogene function prediction by borrowing information from coding genes. Front Genet. (2020) 11:538028. 10.3389/fgene.2020.00807

  • 112.

    HuMAlkhairySLeeIPillichRTFongDSmithKet al. Evaluation of large language models for discovery of gene set function. Nat Methods. (2023) 22:8291. 10.1038/s41592-024-02525-x

  • 113.

    Arango-ArgotyGAHeathLSPrudenAVikeslandPJZhangL. MetaMLP: a fast word embedding based classifier to profile target gene databases in metagenomic samples. J Comput Biol. (2021) 28:106374. 10.1089/cmb.2021.0273

  • 114.

    ToufiqMRinchaiDBettacchioliEKabeerBSAKhanTSubbaBet al. Harnessing large language models (LLMs) for candidate gene prioritization and selection. J Transl Med. (2023) 21:728. 10.1186/s12967-023-04576-8

  • 115.

    AiH. Gsea-sdbe: a gene selection method for breast cancer classification based on gsea and analyzing differences in performance metrics. PLoS ONE. (2022) 17:e0263171. 10.1371/journal.pone.0263171

  • 116.

    HulsenTHuynenMVliegJGroenenP. Untitled. Genome Biol. (2006) 7:R31. 10.1186/gb-2006-7-4-r31

  • 117.

    AbbasiAFAsimMNAhmedSVollmerSDengelA. Survival prediction landscape: an in-depth systematic literature review on activities, methods, tools, diseases, and databases. medRxiv. (2024). p. 202401. 10.1101/2024.01.05.24300889

  • 118.

    DengFYangKZhengG. Period family of clock genes as novel predictors of survival in human cancer: a systematic review and meta-analysis. Dis Markers. (2020) 2020:19. 10.1155/2020/6486238

  • 119.

    ZengWRaoNZhengJWanYWangGLiZet al. Identification of candidate targeted genes in molecular subtypes of gastric cancer. J Biomed Sci Eng. (2017) 10:4553. 10.4236/jbise.2017.105B005

  • 120.

    DağliyanOÜney YüksektepeFKavaklíTürkayM. Optimization based tumor classification from microarray gene expression data. PLoS ONE. (2011) 6:e14579. 10.1371/journal.pone.0014579

  • 121.

    VermaBParkinsonJ. HiTaxon: a hierarchical ensemble framework for taxonomic classification of short reads. Bioinform Adv. (2024) 4:vbae016. 10.1093/bioadv/vbae016

  • 122.

    ShiLChenB. LSHvec: a vector representation of DNA sequences using locality sensitive hashing and FastText word embeddings. In: Proceedings of the 12th ACM Conference on Bioinformatics, Computational Biology, and Health Informatics. (2021). p. 110. 10.1145/3459930.3469521

  • 123.

    MockFKretschmerFKrieseABöckerSMarzM. Taxonomic classification of DNA sequences beyond sequence similarity using deep neural networks. Proc Natl Acad Sci. (2022) 119:e2122636119. 10.1073/pnas.2122636119

  • 124.

    CoenyeTGeversDPeerYVandammePSwingsJ. Towards a prokaryotic genomic taxonomy. FEMS Microbiol Rev. (2005) 29:14767. 10.1016/j.femsre.2004.11.004

  • 125.

    CommichauxSLuanTMuralidharanHSPopM. Database size positively correlates with the loss of species-level taxonomic resolution for the 16S rRNA and other prokaryotic marker genes. bioRxiv. (2023) 14:571439. 10.1101/2023.12.13.571439

  • 126.

    PioGMignonePMagazzúGZampieriGCeciMAngioneC. Integrating genome-scale metabolic modelling and transfer learning for human gene regulatory network reconstruction. Bioinformatics. (2022) 38:48793. 10.1093/bioinformatics/btab647

  • 127.

    DuZZhongXWangFUverskyVN. Inference of gene regulatory networks based on the light gradient boosting machine. Comput Biol Chem. (2022) 101:107769. 10.1016/j.compbiolchem.2022.107769

  • 128.

    PioGCeciMPrisciandaroFMalerbaD. Exploiting causality in gene network reconstruction based on graph embedding. Mach Learn. (2020) 109:123179. 10.1007/s10994-019-05861-8

  • 129.

    LangfelderPHorvathS. WGCNA: an r package for weighted correlation network analysis. BMC Bioinform. (2008) 9:559. 10.1186/1471-2105-9-559

  • 130.

    StanleyDWatson-HaighNCowledCMooreR. Genetic architecture of gene expression in the chicken. BMC Gen. (2013) 14:13. 10.1186/1471-2164-14-13

  • 131.

    Nabeel AsimMAli IbrahimMFazeelADengelAAhmedS. DNA-MP: a generalized DNA modifications predictor for multiple species based on powerful sequence encoding method. Brief Bioinform. (2023) 24:bbac546. 10.1093/bib/bbac546

  • 132.

    AsimMNMalikMIDengelAAhmedS. K-MER neural embedding performance analysis using amino acid codons. In: 2020 International Joint Conference on Neural Networks (IJCNN). IEEE (2020). p. 18. 10.1109/IJCNN48605.2020.9206892

  • 133.

    AsimMNIbrahimMAMalikMIRazzakIDengelAAhmedS. Histone-net: a multi-paradigm computational framework for histone occupancy and modification prediction. Complex Intell Syst. (2023) 9:399419. 10.1007/s40747-022-00802-w

  • 134.

    XieHDingYQianYTiwariPGuoF. Structured sparse regularization based random vector functional link networks for DNA N4-methylcytosine sites prediction. Expert Syst Appl. (2024) 235:121157. 10.1016/j.eswa.2023.121157

  • 135.

    SahaSHalderRKUddinMN. Particle swarm optimization-assisted multilayer ensemble model to predict DNA 4mC sites. Inform Med Unlocked. (2023) 42:101374. 10.1016/j.imu.2023.101374

  • 136.

    ZulfiqarHSunZJHuangQLYuanSSLvHDaoFYet al. Deep-4mCW2V: A sequence-based predictor to identify N4-methylcytosine sites in Escherichia coli. Methods. (2022) 203:55863. 10.1016/j.ymeth.2021.07.011

  • 137.

    FangGZengFLiXYaoL. Word2vec based deep learning network for DNA N4-methylcytosine sites identification. Procedia Comput Sci. (2021) 187:2707. 10.1016/j.procs.2021.04.062

  • 138.

    KhanalJTayaraHZouQChongKT. Identifying DNA N4-methylcytosine sites in the rosaceae genome with a deep learning model relying on distributed feature representation. Comput Struct Biotechnol J. (2021) 19:16129. 10.1016/j.csbj.2021.03.015

  • 139.

    Nguyen-VoTHTrinhQHNguyenLNguyen-HoangPURahardjaSNguyenBP. i4mC-GRU: Identifying DNA N4-Methylcytosine sites in mouse genomes using bidirectional gated recurrent unit and sequence-embedded features. Comput Struct Biotechnol J. (2023) 21:304553. 10.1016/j.csbj.2023.05.014

  • 140.

    ZengRLiaoM. Developing a multi-layer deep learning based predictive model to identify DNA N4-methylcytosine modifications. Front Bioeng Biotechnol. (2020) 8:274. 10.3389/fbioe.2020.00274

  • 141.

    XuHJiaPZhaoZ. Deep4mC: systematic assessment and computational prediction for DNA N4-methylcytosine sites by deep learning. Brief Bioinform. (2021) 22:bbaa099. 10.1093/bib/bbaa099

  • 142.

    LiangYWuYZhangZLiuNPengJTangJ. Hyb4mC: a hybrid DNA2vec-based model for DNA N4-methylcytosine sites prediction. BMC Bioinform. (2022) 23:258. 10.1186/s12859-022-04789-6

  • 143.

    YangSYangZYangJ. 4mCBERT: A computing tool for the identification of DNA N4-methylcytosine sites by sequence-and chemical-derived information based on ensemble learning strategies. Int J Biol Macromol. (2023) 231:123180. 10.1016/j.ijbiomac.2023.123180

  • 144.

    TsukiyamaSHasanMMDengHWKurataH. BERT6mA: prediction of DNA N6-methyladenine site using deep learning-based approaches. Brief Bioinform. (2022) 23:bbac053. 10.1093/bib/bbac053

  • 145.

    ShiHLiSSuX. Plant6mA: a predictor for predicting N6-methyladenine sites with lightweight structure in plant genomes. Methods. (2022) 204:12631. 10.1016/j.ymeth.2022.02.009

  • 146.

    NavarezAMRoxasR. (2022). An evaluation of multitask transfer learning methods in identifying 6mA and 5mC methylation sites of rice and maize. 10.2139/ssrn.4178244

  • 147.

    LiZJiangHKongLChenYLangKFanXet al. Deep6ma: a deep learning framework for exploring similar patterns in DNA n6-methyladenine sites across different species. PLoS Comput Biol. (2021) 17:e1008767. 10.1371/journal.pcbi.1008767

  • 148.

    ZhouZXiaoCYinJSheJDuanHLiuCet al. PSAC-6mA: 6mA site identifier using self-attention capsule network based on sequence-positioning. Comput Biol Med. (2024) 171:108129. 10.1016/j.compbiomed.2024.108129

  • 149.

    FanXQLinBHuJGuoZY. I-DNAN6mA: accurate identification of DNA N6-methyladenine sites using the base-pairing map and deep learning. J Chem Inf Model. (2023) 63:107686. 10.1021/acs.jcim.2c01465

  • 150.

    AbbasZRehmanMUChongKT. TC-6mA-Pred: prediction of DNA N6-methyladenine sites using CNN with transformer. In: 2022 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). IEEE (2022). p. 25062510. 10.1109/BIBM55620.2022.9995083

  • 151.

    HuangQZhouWGuoFXuLZhangL. 6mA-Pred: identifying DNA N6-methyladenine sites based on deep learning. PeerJ. (2021) 9:e10813. 10.7717/peerj.10813

  • 152.

    StanojevićDLiZFooRŠikićM. Rockfish: a transformer-based model for accurate 5-methylcytosine prediction from nanopore sequencing. bioRxiv. (2022). p. 202211. 10.1101/2022.11.11.513492

  • 153.

    TranTAPhamDMOuYYet al. An extensive examination of discovering 5-Methylcytosine Sites in Genome-Wide DNA Promoters using machine learning based approaches. IEEE/ACM Trans Comput Biol Bioinform. (2021) 19:8794. 10.1109/TCBB.2021.3082184

  • 154.

    YuYHeWJinJXiaoGCuiLZengRet al. iDNA-ABT: advanced deep learning model for detecting DNA methylation with adaptive features and transductive information maximization. Bioinformatics. (2021) 37:460310. 10.1093/bioinformatics/btab677

  • 155.

    ZhuoLWangRFuXYaoX. StableDNAm: towards a stable and efficient model for predicting DNA methylation based on adaptive feature correction learning. BMC Gen. (2023) 24:742. 10.1186/s12864-023-09802-7

  • 156.

    ZengWGautamAHusonDH. MuLan-Methyl-multiple transformer-based language models for accurate DNA methylation prediction. GigaScience. (2023) 12:giad054. 10.1093/gigascience/giad054

  • 157.

    JinJYuYWangRZengXPangCJiangYet al. iDNA-ABF: multi-scale deep biological language learning model for the interpretable prediction of DNA methylations. Genome Biol. (2022) 23:219. 10.1186/s13059-022-02780-1

  • 158.

    WangZXiangSZhouCXuQ. DeepMethylation: a deep learning based framework with GloVe and Transformer encoder for DNA methylation prediction. PeerJ. (2023) 11:e16125. 10.7717/peerj.16125

  • 159.

    JeongYGerhäuserCSauterGSchlommTRohrKLutsikP. MethylBERT: a Transformer-based model for read-level DNA methylation pattern identification and tumour deconvolution. bioRxiv. (2023). p. 202310. 10.1101/2023.10.29.564590

  • 160.

    SongCDiaoJBrungerAQuakeS. Simultaneous single-molecule epigenetic imaging of DNA methylation and hydroxymethylation. Proc Nat Acad Sci. (2016) 113:433843. 10.1073/pnas.1600223113

  • 161.

    MaX. Deep5hmc: predicting genome-wide 5-hydroxymethylcytosine landscape via a multimodal deep learning model. Bioinformatics. (2024) 40:btae528. 10.1093/bioinformatics/btae528

  • 162.

    KohliRZhangY. Tet enzymes, TDG and the dynamics of DNA demethylation. Nature. (2013) 502:4729. 10.1038/nature12750

  • 163.

    GialitsisNGiannakopoulosGAthanasouliM. Evaluation of distributed DNA representations on the classification of conserved non-coding elements. In: 11th Hellenic Conference on Artificial Intelligence. (2020). p. 4147. 10.1145/3411408.3411463

  • 164.

    AkalınFYumuşakN. Classification of exon and intron regions on DNA sequences with hybrid use of sbert and anfis approaches. Politeknik Dergisi. (2023) 27:10431053. 10.2339/politeknik.1187808

  • 165.

    DoDTLeNQK. A sequence-based approach for identifying recombination spots in Saccharomyces cerevisiae by using hyper-parameter optimization in FastText and support vector machine. Chemometr Intell Labor Syst. (2019) 194:103855. 10.1016/j.chemolab.2019.103855

  • 166.

    SoldnerFJaenischR. Dissecting risk haplotypes in sporadic Alzheimer's disease. Cell Stem Cell. (2015) 16:3412. 10.1016/j.stem.2015.03.010

  • 167.

    PaigenKPetkovP. Mammalian recombination hot spots: properties, control and evolution. Nat Rev Genet. (2010) 11:22133. 10.1038/nrg2712

  • 168.

    LiZTianDWangBWangJWangSChenHet al. Microbes drive global soil nitrogen mineralization and availability. Glob Chang Biol. (2019) 25:107888. 10.1111/gcb.14557

  • 169.

    KnopsJBradleyKWedinD. Mechanisms of plant species impacts on ecosystem nitrogen cycling. Ecol Lett. (2002) 5:45466. 10.1046/j.1461-0248.2002.00332.x

  • 170.

    SantosPFangZMasonSSetúbalJDixonR. Distribution of nitrogen fixation and nitrogenase-like sequences amongst microbial genomes. BMC Gen. (2012) 13:162. 10.1186/1471-2164-13-162

  • 171.

    NarayananSRamachandranAAakurSNBagavathiA. Genome sequence classification for animal diagnostics with graph representations and deep neural networks. arXiv preprint arXiv:200712791. (2020).

  • 172.

    ZhouFGanRZhangFRenCYuLSiYet al. PHISDetector: a tool to detect diverse in silico phage-host interaction signals for virome studies. Gen Proteom Bioinform. (2022) 20:50823. 10.1016/j.gpb.2022.02.003

  • 173.

    YılmazA. Assessment of mutation susceptibility in DNA sequences with word vectors. J Intel Syst. (2020) 3:16. 10.38016/jista.674910

  • 174.

    ChenZLiXYangMZhangHXuX. Optimization of deep learning models for the prediction of gene mutations using unsupervised clustering. J Pathol Clin Res. (2022) 9:317. 10.1002/cjp2.302

  • 175.

    QiuJNieWDingHDaiJWeiYLiDet al. PB-LKS: a python package for predicting phage-bacteria interaction through local K-MER strategy. Brief Bioinform. (2024) 25:bbae010. 10.1093/bib/bbae010

  • 176.

    WangYSunHWangHLiDZhaoWJiangXet al. An effective model for predicting phage-host interactions via graph embedding representation learning with multi-head attention mechanism. IEEE J Biomed Health Inf. (2023) 27:30613071. 10.1109/JBHI.2023.3261319

  • 177.

    PanJYouWLuXWangSYouZSunYet al. A novel deep learning model for predicting phage-host interactions via multiple biological information. Comput Struct Biotechnol J. (2023) 21:340413. 10.1016/j.csbj.2023.06.014

  • 178.

    EdwardsRMcNairKFaustKRaesJDutilhB. Computational approaches to predict bacteriophage-host relationships. FEMS Microbiol Rev. (2015) 40:25872. 10.1093/femsre/fuv048

  • 179.

    HarbeckNThomssenC. A new look at node-negative breast cancer. Oncologist. (2010) 15:2938. 10.1634/theoncologist.2010-S5-29

  • 180.

    OsseniMATossouPLavioletteFCorbeilJ. MOT: a Multi-Omics Transformer for multiclass classification tumour types predictions. BioRxiv. (2022). p. 202211. 10.1101/2022.11.14.516459

  • 181.

    GhoshASinghTSinglaVBaggaRSrinivasanRKhandelwalN. DTI histogram parameters correlate with the extent of myoinvasion and tumor type in endometrial carcinoma: a preliminary analysis. Acta Radiol. (2019) 61:67584. 10.1177/0284185119875019

  • 182.

    KulakovskiyIVVorontsovIEYevshinISSharipovRNFedorovaADRumynskiyEIet al. HOCOMOCO: towards a complete collection of transcription factor binding models for human and mouse via large-scale ChIP-Seq analysis. Nucleic Acids Res. (2018) 46:D2529. 10.1093/nar/gkx1106

  • 183.

    PujarSO'LearyNAFarrellCMLovelandJEMudgeJMWallinCet al. Consensus coding sequence (CCDS) database: a standardized set of human and mouse protein-coding regions supported by expert curation. Nucleic Acids Res. (2018) 46:D2218. 10.1093/nar/gkx1031

  • 184.

    LiberzonASubramanianAPinchbackRThorvaldsdóttirHTamayoPMesirovJP. Molecular signatures database (MSigDB) 30. Bioinformatics. (2011) 27:173940. 10.1093/bioinformatics/btr260

  • 185.

    ShimadaKMuhlichJLMitchisonTJ. A tool for browsing the Cancer Dependency Map reveals functional connections between genes and helps predict the efficacy and selectivity of candidate cancer drugs. bioRxiv. (2019). p. 201912. 10.1101/2019.12.13.874776

  • 186.

    FornesOCastro-MondragonJAKhanAVan der LeeRZhangXRichmondPAet al. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. (2020) 48:D8792. 10.1093/nar/gkz1001

  • 187.

    ZhangROuHYZhangCT. DEG: a database of essential genes. Nucl Acids Res. (2004) 32:D271D272. 10.1093/nar/gkh024

  • 188.

    deSouza N. The ENCODE project. Nat Methods. (2012) 9:10461046. 10.1038/nmeth.2238

  • 189.

    El AllaliARoseJR. MGC: Gene calling in metagenomic sequences. In: 8 TH International Symposium On Bioinformatics Research and Applications.

  • 190.

    SaxonovSDaizadehIFedorovAGilbertW. EID the Exon-Intron Database–an exhaustive database of protein-coding intron-containing genes. Nucleic Acids Res. (2000) 28:18590. 10.1093/nar/28.1.185

  • 191.

    HubbardTBarkerDBirneyECameronGChenYClarkLet al. The Ensembl genome database project. Nucleic Acids Res. (2002) 30:3841. 10.1093/nar/30.1.38

  • 192.

    SalgadoHSantosAGarza-RamosUvan HeldenJDíazECollado-VidesJ. RegulonDB (version 20): a database on transcriptional regulation in Escherichia coli. Nucl Acids Res. (1999) 27:5960. 10.1093/nar/27.1.59

  • 193.

    PérierRCPrazVJunierTBonnardCBucherP. The eukaryotic promoter database (EPD). Nucleic Acids Res. (2000) 28:302303. 10.1093/nar/28.1.302

  • 194.

    SuWLiuMLYangYHWang JS LiSHLvHet al. PPD: a manually curated database for experimentally verified prokaryotic promoters. J Mol Biol. (2021) 433:166860. 10.1016/j.jmb.2021.166860

  • 195.

    LiSLiLMengXSunPLiuYSongYet al. DREAM: a database of experimentally supported protein-coding RNAs and drug associations in human cancer. Mol Cancer. (2021) 20:16. 10.1186/s12943-021-01436-1

  • 196.

    HuBZhengLLongCSongMLiTYangLet al. EmExplorer: a database for exploring time activation of gene expression in mammalian embryos. Open Biol. (2019) 9:190054. 10.1098/rsob.190054

  • 197.

    KarczewskiKFrancioliL. The genome aggregation database (gnomAD). MacArthur Lab. (2017). p. 110.

  • 198.

    GaoFLuoHZhangCT. DeOri: a database of eukaryotic DNA replication origins. Bioinformatics. (2012) 28:15512. 10.1093/bioinformatics/bts151

  • 199.

    YangJRoyAZhangY. BioLiP: a semi-manually curated database for biologically relevant ligand-protein interactions. Nucleic Acids Res. (2012) 41:D1096103. 10.1093/nar/gks966

  • 200.

    BunielloAMacArthurJALCerezoMHarrisLWHayhurstJMalangoneCet al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. (2019) 47:D100512. 10.1093/nar/gky1120

  • 201.

    CalvinC. Descartes: A new generation system for neutronic calculations. In: International Topical Meeting on Mathematics and Computation, Supercomputing, Reactor Physics and Nuclear and Biological Applications (M and C 2005). Avignon (2005).

  • 202.

    GaoTQianJ. EnhancerAtlas 2.0: an updated resource with enhancer annotation in 586 tissue/cell types across nine species. Nucleic Acids Res. (2020) 48:D5864. 10.1093/nar/gkz980

  • 203.

    TateJGBamfordSJubbHCSondkaZBeareDMBindalNet al. COSMIC: the catalogue of somatic mutations in cancer. Nucleic Acids Res. (2019) 47:D9417. 10.1093/nar/gky1015

  • 204.

    Pi neroJQueralt-RosinachNBravoADeu-PonsJBauer-MehrenABaronMet al. DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes. Database. (2015) 2015:bav028. 10.1093/database/bav028

  • 205.

    LandrumMJLeeJMBensonMBrownGChaoCChitipirallaSet al. ClinVar: public archive of interpretations of clinically relevant variants. Nucleic Acids Res. (2016) 44:D8628. 10.1093/nar/gkv1222

  • 206.

    GhandiMHuangFWJané-ValbuenaJKryukovGVLoCCMcDonald IIIERet al. Next-generation characterization of the cancer cell line encyclopedia. Nature. (2019) 569:5038. 10.1038/s41586-019-1186-3

  • 207.

    FrankishADiekhansMJungreisILagardeJLovelandJEMudgeJMet al. GENCODE 2021. Nucleic Acids Res. (2021) 49:D91623. 10.1093/nar/gkaa1087

  • 208.

    ConsortiumGO. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. (2004) 32:D258D261. 10.1093/nar/gkh036

  • 209.

    SuzukiYYamashitaRNakaiKSuganoSDBTSS. DataBase of human transcriptional start sites and full-length cDNAs. Nucleic Acids Res. (2002) 30:32831. 10.1093/nar/30.1.328

  • 210.

    CloughEBarrettT. The gene expression omnibus database. In:MathéE.DavisS., editor. Statistical Genomics. Methods in Molecular Biology. New York, NY: Humana Press.

  • 211.

    KanehisaM. The KEGG database. In: ‘In silico' simulation of biological processes: Novartis Foundation Symposium. Wiley Online Library (2002). p. 91103. 10.1002/0470857897.ch8

  • 212.

    GeerLYMarchler-BauerAGeerRCHanLHeJHeSet al. The NCBI biosystems database. Nucleic Acids Res. (2010) 38:D492D496. 10.1093/nar/gkp858

  • 213.

    BensonDAKarsch-MizrachiILipmanDJOstellJRappBAWheelerDL. GenBank. Nucleic Acids Res. (2000) 28:158. 10.1093/nar/28.1.15

  • 214.

    SherrySTWardMSirotkinK. dbSNP–database for single nucleotide polymorphisms and other classes of minor genetic variation. Genome Res. (1999) 9:6779. 10.1101/gr.9.8.677

  • 215.

    SherrySTWardMHKholodovMBakerJPhanLSmigielskiEMet al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. (2001) 29:30811. 10.1093/nar/29.1.308

  • 216.

    HamoshAScottAFAmbergerJValleDMcKusickVA. Online Mendelian inheritance in man (OMIM). Hum Mutat. (2000) 15:5761. 10.1002/(SICI)1098-1004(200001)15:1<57::AID-HUMU12>3.0.CO;2-G

  • 217.

    LillebergJZhuYZhangY. Support vector machines and word2vec for text classification with semantic features. In: 2015 IEEE 14th International Conference on Cognitive Informatics & Cognitive Computing (ICCI* CC). IEEE (2015). p. 136140. 10.1109/ICCI-CC.2015.7259377

  • 218.

    ZouWYSocherRCerDManningCD. Bilingual word embeddings for phrase-based machine translation. In: Proceedings of the 2013 Conference on Empirical Methods in Natural Language Processing. (2013). p. 13931398. 10.18653/v1/D13-1141

  • 219.

    Al-AminMIslamMSUzzalSD. Sentiment analysis of Bengali comments with Word2Vec and sentiment information of words. In: 2017 International Conference on Electrical, Computer and Communication Engineering (ECCE). IEEE (2017). p. 186190. 10.1109/ECACE.2017.7912903

  • 220.

    VermaPKAgrawalPAmorimIProdanR. WELFake: Word embedding over linguistic features for fake news detection. IEEE Trans Comput Soc Syst. (2021) 8:88193. 10.1109/TCSS.2021.3068519

  • 221.

    MikolovTChenKCorradoGDeanJ. Efficient estimation of word representations in vector space. arXiv preprint arXiv:13013781. (2013).

  • 222.

    PenningtonJSocherRManningCD. Glove: Global vectors for word representation. In: Proceedings of the 2014 Conference on Empirical Methods in Natural Language Processing (EMNLP). (2014). p. 15321543. 10.3115/v1/D14-1162

  • 223.

    MikolovTGraveEBojanowskiPPuhrschCJoulinA. Advances in pre-training distributed word representations. arXiv preprint arXiv:171209405. (2017).

  • 224.

    PerozziBAl-RfouRSkienaS. Deepwalk: online learning of social representations. In: Proceedings of the 20th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. (2014). p. 701710. 10.1145/2623330.2623732

  • 225.

    GroverALeskovecJ. node2vec: Scalable feature learning for networks. In: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. (2016). p. 855864. 10.1145/2939672.2939754

  • 226.

    NarayananAChandramohanMVenkatesanRChenLLiuYJaiswalS. graph2vec: Learning distributed representations of graphs. arXiv preprint arXiv:170705005. (2017).

  • 227.

    WangDCuiPZhuW. Structural deep network embedding. In: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (2016). p. 12251234. 10.1145/2939672.2939753

  • 228.

    TangLLiuH. Relational learning via latent social dimensions. In: Proceedings of the 15th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. (2009). p. 817826. 10.1145/1557019.1557109

  • 229.

    CaoSLuWXuQ. Grarep: learning graph representations with global structural information. In: Proceedings of the 24th ACM International on Conference on Information and Knowledge Management. (2015). p. 891900. 10.1145/2806416.2806512

  • 230.

    BelkinMNiyogiP. Laplacian eigenmaps and spectral techniques for embedding and clustering. In: Advances in Neural Information Processing Systems. (2001). p. 14. 10.7551/mitpress/1120.003.0080

  • 231.

    RoweisSTSaulLK. Nonlinear dimensionality reduction by locally linear embedding. Science. (2000) 290:23236. 10.1126/science.290.5500.2323

  • 232.

    SmailiFZGaoXHoehndorfR. OPA2Vec: combining formal and informal content of biomedical ontologies to improve similarity-based prediction. Bioinformatics. (2019) 35:213340. 10.1093/bioinformatics/bty933

  • 233.

    SongTSongHPanZGaoYYangQWangX. DeepDualEPI: predicting promoter-enhancer interactions based on DNA sequence and genomic signals. In: 2023 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). IEEE (2023). p. 28892895. 10.1109/BIBM58861.2023.10385972

  • 234.

    LeNQK. iN6-methylat (5-step): identifying DNA N 6-methyladenine sites in rice genome using continuous bag of nucleobases via Chou's 5-step rule. Molec Gen Genom. (2019) 294:117382. 10.1007/s00438-019-01570-y

  • 235.

    RazaATahirMAlamW. iPro-TCN: Prediction of DNA promoters recognition and their strength using temporal convolutional network. IEEE Access. (2023). 10.1109/ACCESS.2023.3285197

  • 236.

    WangYTaiSZhangSShengNXieX. PromGER: promoter prediction based on graph embedding and ensemble learning for eukaryotic sequence. Genes. (2023) 14:1441. 10.3390/genes14071441

  • 237.

    MinBRossHSulemEVeysehAPBNguyenTHSainzOet al. Recent advances in natural language processing via large pre-trained language models: a survey. ACM Comput Surv. (2023) 56:140. 10.1145/3605943

  • 238.

    ChangYWangXWangJWuYYangLZhuKet al. A survey on evaluation of large language models. ACM Trans Intell Syst Technol. (2024) 15:145. 10.1145/3641289

  • 239.

    YueTWangYZhangLGuCXueHWangWet al. Deep learning for genomics: from early neural nets to modern large language models. Int J Mol Sci. (2023) 24:15858. 10.3390/ijms242115858

  • 240.

    ChenYZouJ. GenePT: a simple but effective foundation model for genes and cells built from ChatGPT. bioRxiv. (2023). 10.1101/2023.10.16.562533

  • 241.

    VaswaniAShazeerNParmarNUszkoreitJJonesLGomezANet al. Attention is all you need. In: Advances in Neural Information Processing Systems. (2017). p. 30.

  • 242.

    YangZDaiZYangYCarbonellJSalakhutdinovRLeQV. XLNet: generalized autoregressive pretraining for language understanding. arXiv preprint arXiv:190608237. (2019).

  • 243.

    HowardJRuderS. Universal language model fine-tuning for text classification. arXiv preprint arXiv:180106146. (2018).

  • 244.

    DevlinJChangMWLeeKToutanovaK. Bert: Pre-training of deep bidirectional transformers for language understanding. arXiv preprint arXiv:181004805. (2018).

  • 245.

    LanZChenMGoodmanSGimpelKSharmaPSoricutR. Albert: a lite bert for self-supervised learning of language representations. arXiv preprint arXiv:190911942. (2019).

  • 246.

    ClarkKLuongMTLeQVManningCD. Electra: Pre-training text encoders as discriminators rather than generators. arXiv preprint arXiv:200310555. (2020).

  • 247.

    BrownTMannBRyderNSubbiahMKaplanJDDhariwalPet al. Language models are few-shot learners. Adv Neural Inf Process Syst. (2020) 33:1877901. 10.48550/arXiv.2005.14165

  • 248.

    HusonDZengW. MR-DNA: flexible 5mC-methylation-site recognition in DNA sequences using token classification. bioRxiv. (2023).

  • 249.

    YangMHuangLHuangHTangHZhangNYangHet al. Integrating convolution and self-attention improves language model of human genome for interpreting non-coding regions at base-resolution. Nucleic Acids Res. (2022) 50:e81e81. 10.1093/nar/gkac326

  • 250.

    ZhouJChenQBraunPRPerzel MandellKAJaffeAETanHYet al. Deep learning predicts DNA methylation regulatory variants in the human brain and elucidates the genetics of psychiatric disorders. Proc Nat Acad Sci. (2022) 119:e2206069119. 10.1073/pnas.2206069119

  • 251.

    LigetiBSzepesi-NagyIBodnárBLigeti-NagyNJuhászJ. ProkBERT family: genomic language models for microbiome applications. Front Microbiol. (2024) 14:1331233. 10.3389/fmicb.2023.1331233

  • 252.

    CarvalhoVNepomucenoTPoletoTTuretJCostaA. Mining public opinions on COVID-19 vaccination: a temporal analysis to support combating misinformation. Trop Med Infect Dis. (2022) 7:256. 10.3390/tropicalmed7100256

  • 253.

    LestariN. A comparison of artificial neural network and naive Bayes classification using unbalanced data handling. Barekeng Jurnal Ilmu Matematika Dan Terapan. (2023) 17:158594. 10.30598/barekengvol17iss3pp1585-1594

  • 254.

    AyyadevaraV. Random forest. (2018). p. 105116. 10.1007/978-1-4842-3564-5_5

  • 255.

    WeitschekEFisconGFeliciG. Supervised DNA barcodes species classification: analysis, comparisons and results. Biodata Min. (2014) 7:4. 10.1186/1756-0381-7-4

  • 256.

    GataWBayhaqyA. Analysis sentiment about islamophobia when Christchurch attack on social media. Telkomnika. (2020) 18:1819. 10.12928/telkomnika.v18i4.14179

  • 257.

    TamrakarSBalBThapaR. Aspect based sentiment analysis of Nepali text using support vector machine and naive Bayes. Technical Journal. (2020) 2:229. 10.3126/tj.v2i1.32824

  • 258.

    MuludiKAkbarMShofianaDSyarifA. Sentiment analysis of energy independence tweets using simple recurrent neural network. IJCCS. (2021) 15:339. 10.22146/ijccs.66016

  • 259.

    ArdiantoRRivanieTAlkhalifiYNugrahaFGataW. Sentiment analysis on e-sports for education curriculum using naive Bayes and support vector machine. Jurnal Ilmu Komputer Dan Informasi. (2020) 13:10922. 10.21609/jiki.v13i2.885

  • 260.

    ArilyaRAzharYChandranegaraD. Sentiment analysis on work from home policy using naïve bayes method and particle swarm optimization. Jurnal Ilmiah Teknik Elektro Komputer Dan Informatika. (2021) 7:433. 10.26555/jiteki.v7i3.22080

  • 261.

    KamathUShehuAJongKA. two-stage evolutionary approach for effective classification of hypersensitive DNA sequences. J Bioinform Comput Biol. (2011) 09:399413. 10.1142/S0219720011005586

  • 262.

    JannahN. Comparison of naïve Bayes and SVM in sentiment analysis of product reviews on marketplaces. Sinkron. (2024) 8:72733. 10.33395/sinkron.v8i2.13559

  • 263.

    GunasekaranHRamalakshmiKArokiarajAKanmaniSDhasC. Analysis of DNA sequence classification using CNN and hybrid models. Comput Math Methods Med. (2021) 2021:112. 10.1155/2021/1835056

  • 264.

    AsimMNMalikMIZeheCTryggJDengelAAhmedSet al. robust and precise convnet for small non-coding RNA classification (RPC-snrc). IEEE Access. (2020) 9:1937990. 10.1109/ACCESS.2020.3037642

  • 265.

    AsimMNIbrahimMAImran MalikMDengelAAhmedS. Circ-LocNet: a computational framework for circular RNA sub-cellular localization prediction. Int J Mol Sci. (2022) 23:8221. 10.3390/ijms23158221

  • 266.

    ChongSPeytonPJ. A meta-analysis of the accuracy and precision of the ultrasonic cardiac output monitor (USCOM). Anaesthesia. (2012) 67:126671. 10.1111/j.1365-2044.2012.07311.x

  • 267.

    AsimMNIbrahimMAMalikMIZeheCCloarecOTryggJet al. EL-RMLocNet: An explainable LSTM network for RNA-associated multi-compartment localization prediction. Comput Struct Biotechnol J. (2022) 20:39864002. 10.1016/j.csbj.2022.07.031

  • 268.

    SaloomRHKhafajiHK. Mutation types and pathogenicity classification using multi-label multi-class deep networks. In: AIP Conference Proceedings. AIP Publishing (2024). 10.1063/5.0213291

  • 269.

    ButuceaCNdaoudMStepanovaNATsybakovAB. Variable selection with Hamming loss. Ann Statist. (2018) 46:18371875. 10.1214/17-AOS1572

  • 270.

    AllenDM. Mean square error of prediction as a criterion for selecting variables. Technometrics. (1971) 13:46975. 10.1080/00401706.1971.10488811

  • 271.

    HodsonTO. Root mean square error (RMSE) or mean absolute error (MAE): when to use them or not. Geosci Model Dev Disc. (2022) 2022:110. 10.5194/gmd-2022-64

  • 272.

    KatoT. Prediction of photovoltaic power generation output and network operation. In: Integration of Distributed Energy Resources in Power Systems. Elsevier (2016). p. 77108. 10.1016/B978-0-12-803212-1.00004-0

  • 273.

    KimSKimHA. new metric of absolute percentage error for intermittent demand forecasts. Int J Forecast. (2016) 32:66979. 10.1016/j.ijforecast.2015.12.003

  • 274.

    CameronACWindmeijerFA. An R-squared measure of goodness of fit for some common nonlinear regression models. J Econom. (1997) 77:32942. 10.1016/S0304-4076(96)01818-0

  • 275.

    AzharMBlancPAsimMImranSHayatNShahidHet al. The evaluation of reanalysis and analysis products of solar radiation for Sindh province, Pakistan. Renewable Energy. (2020) 145:34762. 10.1016/j.renene.2019.04.107

  • 276.

    DudoitSGentlemanR. Cluster analysis in DNA microarray experiments. In: Bioconductor Short Course Winter. (2002). p. 506.

  • 277.

    RahmanianMMansooriEG. An unsupervised gene selection method based on multivariate normalized mutual information of genes. Chemometr Intell Lab Syst. (2022) 222:104512. 10.1016/j.chemolab.2022.104512

  • 278.

    LordEDialloABMakarenkovV. Classification of bioinformatics workflows using weighted versions of partitioning and hierarchical clustering algorithms. BMC Bioinform. (2015) 16:119. 10.1186/s12859-015-0508-1

  • 279.

    HossenMBAuwulMR. Comparative study of K-means, partitioning around medoids, agglomerative hierarchical, and DIANA clustering algorithms by using cancer datasets. Biomed Stat Inform. (2020) 5:20. 10.11648/j.bsi.20200501.14

  • 280.

    SinghAKMittalSMalhotraPSrivastavaYV. Clustering evaluation by davies-bouldin index (DBI) in cereal data using k-means. In: 2020 Fourth international conference on computing methodologies and communication (ICCMC). IEEE (2020). p. 306310. 10.1109/ICCMC48092.2020.ICCMC-00057

  • 281.

    LeNQKHoQT. Deep transformers and convolutional neural network in identifying DNA N6-methyladenine sites in cross-species genomes. Methods. (2022) 204:199206. 10.1016/j.ymeth.2021.12.004

  • 282.

    WangSLiuYLiuYZhangYZhuX. BERT-5mC: an interpretable model for predicting 5-methylcytosine sites of DNA based on BERT. PeerJ. (2023) 11:e16600. 10.7717/peerj.16600

  • 283.

    MartinekVCechakDGresovaKAlexiouPSimecekP. Fine-tuning transformers for genomic tasks. bioRxiv. (2022). p. 202202. 10.1101/2022.02.07.479412

  • 284.

    LiXWangGAWeiZWangHZhuX. Protein-DNA interface hotspots prediction based on fusion features of embeddings of protein language model and handcrafted features. Comput Biol Chem. (2023) 107:107970. 10.1016/j.compbiolchem.2023.107970

  • 285.

    LiWLiuWGuoYWangBQingH. Deep contextual representation learning for identifying essential proteins via integrating multisource protein features. Chinese J Electr. (2023) 32:86881. 10.23919/cje.2022.00.053

  • 286.

    YanYLiWWangSHuangT. Seq-RBPPred: predicting RNA-binding proteins from sequence. ACS Omega. (2024) 9:1273442. 10.1021/acsomega.3c08381

  • 287.

    ZhangXGuoHZhangFWangXWuKQiuSet al. HNetGO: protein function prediction via heterogeneous network transformer. Brief Bioinform. (2023) 24:bbab556. 10.1093/bib/bbab556

  • 288.

    JhaKSahaSKarmakarS. Prediction of protein-protein interactions using vision transformer and language model. IEEE/ACM Trans Comput Biol Bioinform. (2023) 20:321525. 10.1109/TCBB.2023.3248797

  • 289.

    Yu DL YuZGHan GS LiJAnhV. Heterogeneous types of miRNA-disease associations stratified by multi-layer network embedding and prediction. Biomedicines. (2021) 9:1152. 10.3390/biomedicines9091152

  • 290.

    CaiKZhuY. A method for identifying essential proteins based on deep convolutional neural network architecture with particle swarm optimization. In: 2022 Asia Conference on Advanced Robotics, Automation, and Control Engineering (ARACE). IEEE (2022). p. 712. 10.1109/ARACE56528.2022.00010

  • 291.

    DuanTKuangZWangJMaZ. GBDTLRL2D predicts LncRNA-disease associations using MetaGraph2Vec and K-means based on heterogeneous network. Front Cell Dev Biol. (2021) 9:753027. 10.3389/fcell.2021.753027

  • 292.

    WangHZhengHChenDZTANGO. A GO-term embedding based method for protein semantic similarity prediction. IEEE/ACM Trans Comput Biol Bioinform. (2022) 20:694706. 10.1109/TCBB.2022.3143480

  • 293.

    Wei MM YuCQLiLPYouZHRenZHGuanYJet al. LPIH2V: LncRNA-protein interactions prediction using HIN2Vec based on heterogeneous networks model. Front Genet. (2023) 14:1122909. 10.3389/fgene.2023.1122909

  • 294.

    ZhouJRYouZHChengLJiBY. Prediction of lncRNA-disease associations via an embedding learning HOPE in heterogeneous information networks. Molec Ther-Nucleic Acids. (2021) 23:27785. 10.1016/j.omtn.2020.10.040

  • 295.

    WangZGuYZhengSYangLLiJ. MGREL: a multi-graph representation learning-based ensemble learning method for gene-disease association prediction. Comput Biol Med. (2023) 155:106642. 10.1016/j.compbiomed.2023.106642

  • 296.

    WangJZhangJCaiYDengL. Deepmir2go: inferring functions of human micrornas using a deep multi-label classification model. Int J Mol Sci. (2019) 20:6046. 10.3390/ijms20236046

  • 297.

    SuXRHuLYouZHHuPWZhaoBW. Multi-view heterogeneous molecular network representation learning for protein-protein interaction prediction. BMC Bioinform. (2022) 23:234. 10.1186/s12859-022-04766-z

  • 298.

    WangLWuMWuYZhangXLiSHeMet al. Prediction of the disease causal genes based on heterogeneous network and multi-feature combination method. Comput Biol Chem. (2022) 97:107639. 10.1016/j.compbiolchem.2022.107639

  • 299.

    WuKZhouDSlonimDHuXCowenL. MELISSA: semi-supervised embedding for protein function prediction across multiple networks. bioRxiv. (2023). p. 202308.

  • 300.

    WanCCozzettoDFaRJonesDT. Using deep maxout neural networks to improve the accuracy of function prediction from protein interaction networks. PLoS ONE. (2019) 14:e0209958. 10.1371/journal.pone.0209958

  • 301.

    MadedduLStiloGVelardiP. Network-based methods for disease-gene prediction. arXiv preprint arXiv:190210117. (2019). 10.1504/IJDMB.2020.109502

  • 302.

    CaoWChenYYangJYXueFYYuZHFengJet al. Metapath-aggregated multilevel graph embedding for miRNA-disease association prediction. In: 2023 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). IEEE (2023). p. 468473. 10.1109/BIBM58861.2023.10385762

  • 303.

    PrabhakarVLiuK. Unsupervised co-optimization of a graph neural network and a knowledge graph embedding model to prioritize causal genes for Alzheimer's Disease. medRxiv. (2022). p. 202210. 10.1101/2022.10.03.22280657

  • 304.

    TianZHanCXuLTengZSongW. MGCNSS: miRNA-disease association prediction with multi-layer graph convolution and distance-based negative sample selection strategy. Brief Bioinform. (2024) 25:bbae168. 10.1093/bib/bbae168

  • 305.

    WangXFYuCQYouZHQiaoYLiZWHuangWZet al. KS-CMI: A circRNA-miRNA interaction prediction method based on the signed graph neural network and denoising autoencoder. Iscience. (2023) 26:107478. 10.1016/j.isci.2023.107478

  • 306.

    ChuXGuanBDaiLLiuJxLiFShangJ. Network embedding framework for driver gene discovery by combining functional and structural information. BMC Genom. (2023) 24:426. 10.1186/s12864-023-09515-x

  • 307.

    ZhouLPengXZengLPengL. Finding potential lncRNA-disease associations using a boosting-based ensemble learning model. Front Genet. (2024) 15:1356205. 10.3389/fgene.2024.1356205

  • 308.

    LiJLiJKongMWangDFuKShiJet al. predicting lncRNA-disease associations by Singular Value Decomposition and node2vec. BMC Bioinform. (2021) 22:118. 10.1186/s12859-021-04457-1

  • 309.

    MallickKBandyopadhyaySChakrabortySChoudhuriRBoseS. Topo2vec: a novel node embedding generation based on network topology for link prediction. IEEE Trans Comput Soc Syst. (2019) 6:130617. 10.1109/TCSS.2019.2950589

  • 310.

    VilelaJAsifMMarquesARSantosJXRasgaCVicenteAet al. Biomedical knowledge graph embeddings for personalized medicine: Predicting disease-gene associations. Expert Systems. (2023) 40:e13181. 10.1111/exsy.13181

  • 311.

    FrankeJKRungeFKoeksalRBackofenRHutterF. RNAformer: a simple yet effective deep learning model for rna secondary structure prediction. bioRxiv. (2024) p. 202402. 10.1101/2024.02.12.579881

  • 312.

    HuMYuanFYangKJuFSuJWangHet al. Exploring evolution-aware &-free protein language models as protein function predictors. Adv Neural Inf Process Syst. (2022) 35:3887384. 10.48550/arXiv.2206.06583

  • 313.

    ShahSMAOuYY. Disto-TRP: an approach for identifying transient receptor potential (TRP) channels using structural information generated by AlphaFold. Gene. (2023) 871:147435. 10.1016/j.gene.2023.147435

  • 314.

    WangJChenSYuanQChenJLiDWangLet al. Predicting the effects of mutations on protein solubility using graph convolution network and protein language model representation. J Comput Chem. (2024) 45:43645. 10.1002/jcc.27249

  • 315.

    HouXWangYBuDWangYSunS. EMNGly: predicting N-linked glycosylation sites using the language models for feature extraction. Bioinformatics. (2023) 39:btad650. 10.1093/bioinformatics/btad650

  • 316.

    RocheRMoussadBShuvoMHTarafderSBhattacharyaD. EquiPNAS: improved protein-nucleic acid binding site prediction using protein-language-model-informed equivariant deep graph neural networks. Nucleic Acids Res. (2024) 52:e27e27. 10.1093/nar/gkae039

  • 317.

    MaJZhaoZLiTLiuYMaJZhangR. GraphsformerCPI: graph transformer for compound-protein interaction prediction. Interdiscip Sci. (2024) 16:361377. 10.1007/s12539-024-00609-y

  • 318.

    DaiZDengF. LncPNdeep: a long non-coding RNA classifier based on Large Language Model with peptide and nucleotide embedding. bioRxiv. (2023). p. 202311. 10.1101/2023.11.29.569323

  • 319.

    MeierJRaoRVerkuilRLiuJSercuTRivesA. Language models enable zero-shot prediction of the effects of mutations on protein function. Adv Neural Inf Process Syst. (2021) 34:29287303. 10.1101/2021.07.09.450648

  • 320.

    XuSOnodaA. Accurate and fast prediction of intrinsically disordered protein by multiple protein language models and ensemble learning. J Chem Inf Model. (2023) 64:290111. 10.1021/acs.jcim.3c01202

  • 321.

    LiYGuoZWangKGaoXWangG. End-to-end interpretable disease-gene association prediction. Brief Bioinform. (2023) 24:bbad118. 10.1093/bib/bbad118

  • 322.

    ZouHJiBZhangMLiuFXieXPengS. MHGTMDA: molecular heterogeneous graph transformer based on biological entity graph for miRNA-disease associations prediction. Molec Ther-Nucleic Acids. (2024) 35:102139. 10.1016/j.omtn.2024.102139

  • 323.

    MelnykIChenthamarakshanVChenPYDasPDhurandharAPadhiIet al. Reprogramming pretrained language models for antibody sequence infilling. In: International Conference on Machine Learning. PMLR (2023). p. 2439824419.

  • 324.

    SaadatMBehjatiAZare-MirakabadFGharaghaniS. Drug-target binding affinity prediction using transformers. (2021). 10.1101/2021.09.30.462610

  • 325.

    LennoxMRobertsonNDevereuxB. Modelling drug-target binding affinity using a BERT based graph neural network. In: 2021 43rd Annual International Conference of the IEEE Engineering in Medicine & Biology Society (EMBC). IEEE (2021). p. 43484353. 10.1109/EMBC46164.2021.9629695

  • 326.

    ElnaggarAEssamHSalah-EldinWMoustafaWElkerdawyMRochereauCet al. Ankh: optimized protein language model unlocks general-purpose modelling. arXiv preprint arXiv:230106568. (2023).

  • 327.

    HaselbeckFJohnMZhangYPirnayJFuenzalida-WernerJPCostaRDet al. Superior protein thermophilicity prediction with protein language model embeddings. NAR Gen Bioinform. (2023) 5:lqad087. 10.1093/nargab/lqad087

  • 328.

    YuanQTianCSongYOuPZhuMZhaoHet al. GPSFun: geometry-aware protein sequence function predictions with language models. Nucleic Acids Res. (2024) 52:gkae381. 10.1093/nar/gkae381

  • 329.

    WuFYangRZhangCZhangL. A deep learning framework combined with word embedding to identify DNA replication origins. Sci Rep. (2021) 11:844. 10.1038/s41598-020-80670-x

  • 330.

    ZhouYWuTJiangYLiYLiKQuanLet al. DeepNup: prediction of nucleosome positioning from DNA sequences using deep neural network. Genes. (2022) 13:1983. 10.3390/genes13111983

  • 331.

    LeiRJiaJQinLWeiX. iPro2L-DG: Hybrid network based on improved densenet and global attention mechanism for identifying promoter sequences. Heliyon. (2024) 10:e27364. 10.1016/j.heliyon.2024.e27364

  • 332.

    ZhengLLiuLZhuWDingYWuF. Predicting enhancer-promoter interaction based on epigenomic signals. Front Genet. (2023) 14:1133775. 10.3389/fgene.2023.1133775

  • 333.

    HuWLiMXiaoHGuanL. Essential genes identification model based on sequence feature map and graph convolutional neural network. BMC Genomics. (2024) 25:47. 10.1186/s12864-024-09958-w

  • 334.

    ReddyAJHerschlMHKolliSLuAXGengXKumarAet al. Pretraining strategies for effective promoter-driven gene expression prediction. bioRxiv. (2023).

  • 335.

    LeNQKXuL. Optimizing hyperparameter tuning in machine learning to improve the predictive performance of cross-species N6-methyladenosine sites. ACS Omega. (2023) 8:394206. 10.1021/acsomega.3c05074

  • 336.

    DaiZYangZYangYCarbonellJLeQVSalakhutdinovR. Transformer-xl: Attentive language models beyond a fixed-length context. arXiv preprint arXiv:190102860. (2019).

  • 337.

    RadfordANarasimhanKSalimansTSutskeverI. Improving language understanding by generative pre-training. bioRxiv. (2018).

  • 338.

    RadfordAWuJChildRLuanDAmodeiDSutskeverIet al. Language models are unsupervised multitask learners. OpenAI blog. (2019) 1:9.

  • 339.

    AchiamJAdlerSAgarwalSAhmadLAkkayaIAlemanFLet al. Gpt-4 technical report. arXiv preprint arXiv:230308774. (2023).

  • 340.

    RatajczakFJoblinMHildebrandtMRingsquandlMFalter-BraunPHeinigM. Speos: an ensemble graph representation learning framework to predict core gene candidates for complex diseases. Nat Commun. (2023) 14:7206. 10.1038/s41467-023-42975-z

  • 341.

    ChenHRyuJVinyardMELererAPinelloL. SIMBA: single-cell embedding along with features. Nat Methods. (2023) 21:10031013. 10.1038/s41592-023-01899-8

  • 342.

    Han GS LiQLiY. Nucleosome positioning based on DNA sequence embedding and deep learning. BMC Genom. (2022) 23:301. 10.1186/s12864-022-08508-6

  • 343.

    ChenYZouJ. GenePT: a simple but effective foundation model for genes and cells built from ChatGPT. bioRxiv. (2024). 10.1101/2023.10.16.562533

  • 344.

    ZhangDZhangWHeBZhangJQinCYaoJ. Dnagpt: A generalized pretrained tool for multiple DNA sequence analysis tasks. bioRxiv. (2023). p. 202307. 10.1101/2023.07.11.548628

  • 345.

    ElnaggarAHeinzingerMDallagoCRehawiGWangYJonesLet al. Prottrans: Toward understanding the language of life through self-supervised learning. IEEE Trans Pattern Anal Mach Intell. (2021) 44:711227. 10.1109/TPAMI.2021.3095381

  • 346.

    LaiPTLuZ. BERT-GT cross-sentence n-ary relation extraction with BERT and graph transformer. Bioinformatics. (2020) 36:567885. 10.1093/bioinformatics/btaa1087

  • 347.

    CuiHWangCMaanHPangKLuoFDuanNet al. scGPT: toward building a foundation model for single-cell multi-omics using generative AI. Nat Methods. (2024) 21:14701480. 10.1038/s41592-024-02201-0

  • 348.

    LiZJinJLongWWeiL. PLPMpro: enhancing promoter sequence prediction with prompt-learning based pre-trained language model. Comput Biol Med. (2023) 164:107260. 10.1016/j.compbiomed.2023.107260

  • 349.

    LeNQKHoQTNguyenVNChangJS. BERT-Promoter: An improved sequence-based predictor of DNA promoter using BERT pre-trained model and SHAP feature selection. Comput Biol Chem. (2022) 99:107732. 10.1016/j.compbiolchem.2022.107732

  • 350.

    TenekeciSTekirS. Identifying promoter and enhancer sequences by graph convolutional networks. Comput Biol Chem. (2024) 110:108040. 10.1016/j.compbiolchem.2024.108040

  • 351.

    JuHBaiJJiangJCheYChenX. Comparative evaluation and analysis of DNA N4-methylcytosine methylation sites using deep learning. Front Genet. (2023) 14:1254827. 10.3389/fgene.2023.1254827

  • 352.

    HearstMADumaisSTOsunaEPlattJScholkopfB. (1998). Support vector machines. IEEE Intell. Syst.13, 1828.

Summary

Keywords

computational biology, computational genomics, DNA sequence analysis, artificial intelligence, deep learning

Citation

Asim MN, Ibrahim MA, Zaib A and Dengel A (2025) DNA sequence analysis landscape: a comprehensive review of DNA sequence analysis task types, databases, datasets, word embedding methods, and language models. Front. Med. 12:1503229. doi: 10.3389/fmed.2025.1503229

Received

28 September 2024

Accepted

10 March 2025

Published

08 April 2025

Volume

12 - 2025

Edited by

Alice Chen, Consultant, Potomac, MD, United States

Reviewed by

Matt Field, James Cook University, Australia

Ankan Bhattacharya, Hooghly Engineering and Technology College, India

Updates

Copyright

*Correspondence: Muhammad Nabeel Asim

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