OmniGenome: Aligning RNA Sequences with Secondary Structures in Genomic Foundation Models
Abstract: The structures of RNA sequences play a vital role in various cellular processes, while existing genomic foundation models (FMs) have struggled with precise sequence-structure alignment, due to the complexity of exponential combinations of nucleotide bases. In this study, we introduce OmniGenome, a foundation model that addresses this critical challenge of sequence-structure alignment in RNA FMs. OmniGenome bridges the sequences with secondary structures using structure-contextualized modeling, enabling hard in-silico genomic tasks that existing FMs cannot handle, e.g., RNA design tasks. The results on two comprehensive genomic benchmarks show that OmniGenome achieves state-of-the-art performance on complex RNA subtasks. For example, OmniGenome solved of complex puzzles, compared to SpliceBERT which solved only of the puzzles. Besides, OmniGenome solves most of the puzzles within hour, while the existing methods usually allocate hours for each puzzle. Overall, OmniGenome establishes wide genomic application cases and offers profound insights into biological mechanisms from the perspective of sequence-structure alignment.
Keywords: Genomic Foundation Model, RNA Design, Secondary Structure Prediction
1 Introduction
RNA, folded from diverse nucleotide sequences, is a critical type of molecule and enables the flow of genetic information from DNA to protein [1, 2]. Similar to the natural language texts, the nucleotide sequences are regarded as the sentences represented in the “language” of genomics [3]. Therefore, many studies aim to leverage pretrained genomic foundation models, a.k.a., language models (LMs), to address in-silico RNA modeling tasks, such as secondary structure prediction [4, 5, 6, 7], degradation rate prediction [8, 9], mRNA vaccine design[10, 11]. There is a vital impact of secondary structures111The “structure” term refers to “RNA secondary structure” in this work. of RNA molecules on various cellular processes [12], such as gene expression [13, 14], protein translation [15, 16]. However, the performance of existing FMs are unsatisfactory because they only model the RNA sequences without secondary structures. Consequently, leveraging the secondary structures in the RNA modeling is crucial for the FMs in solving complex genomic subtasks, e.g., RNA design [17]. One of the intuitive methods to utilize secondary structures is to build the bidirectional mapping between sequences and structures. We formulate this mapping as sequence-structure alignment, which can be divided into two parts: sequence-to-structure (Seq2Str) and structure-to-sequence (Str2Seq) prediction. Both two predictions are challenging and have not been settled in previous works because the sequences and structures are sophisticated due to exponential combinations of up to thousands of bases. In short, to align the sequence and structure in FM-based genomic modeling, we need to address the Seq2Str and Str2Seq predictions simultaneously.
Seq2Str prediction. Secondary structure prediction is a difficult task because the structures grow exponentially with increasing sequence lengths, and an identical sequence may fold into different sub-optimal structures because the folding principles of RNA sequences depend on a variety of biological factors [18], as indicated in Fig. 1. Moreover, biologically verified structures are unfortunately expensive, compared to raw sequences, to obtain and process. This problem causes an imbalanced data scale between RNA sequences and verified structures making it hard to train an FM for the Seq2Str prediction. This means we need a tremendous number of secondary structures to train the Seq2Str prediction capability. However, existing works on RNA modeling generally utilize the masked language modeling [19] objective adapted from natural language processing (NLP) to pretrain FMs and regard the secondary structure prediction as a fine-tuning downstream task [5, 6], which inevitably shadows the performance of FMs without the pretraining on Seq2Str prediction.
![Refer to caption](x1.png)
To address the data scarcity of secondary structures in Seq2Str, we leverage ViennaRNA [20] to compute the structures for billions of RNA sequences and introduce Seq2Str prediction into the pretraining phase. The pretraining on billions of structures will build an effective and robust mapping from sequences to structures. Our experiments in Section 3.3 show that the FM can achieve state-of-the-art performance even in zero-shot secondary structure prediction.
Str2Seq prediction. In addition to the data scarcity problem aforementioned in Seq2Str, Str2Seq prediction presents a substantial challenge because secondary structures can correspond to heterogeneous sequences as shown in Fig. 1, where the identical structure can be from different sequences. Therefore, it is impossible to reconstruct sequences according to only structures based on structure-to-sequence supervised training because there is no deterministic sequence for each structure. Consequently, the Str2Seq prediction has not been explored in the existing genomic FMs.
![Refer to caption](x2.png)
To overcome the challenge in Str2Seq prediction, we formulate the Str2Seq as a structure-contextualized sequence reconstruction task. Since RNA secondary structures can be represented in texts, we concatenate the sequences and structure pairs and then mask a small portion of the sequence or structure tokens before feeding them into the FM. In other words, we pretrain the FM to reconstruct the missing nucleotide bases given the context of structures, which will build the Str2Seq mapping ability for the FM. To perform this pretraining objective, we expand the vocabulary to represent the structure and deploy corresponding embeddings in the FM, because regular genomic FMs cannot recognize and embed the structures, as shown in the virtual example in Fig. 2.
Benchmark and Evaluation Results. We conduct two comprehensive in-silico genomic understanding benchmarks to evaluate the performance of our FM. The first benchmark compiled in the study is the RNA genomic benchmark (RGB) which contains diverse challenging genomic understanding tasks that benefit from the sequence-structure alignment, such as secondary structure prediction. The second benchmark is the plant genomic benchmark (PGB) which contains millions of DNA sequences to evaluate the DNA sequence understanding tasks. This benchmark is used to evaluate the generalizability of our FM among diversified species and genomes. The overall performance of our FM (up to M parameters) on both two benchmarks consistently outperforms existing genomic FMs with up to improvement even compared with agro-NT [14] that contains billion parameters.
We also conduct zero-shot Seq2Str and Str2Seq prediction experiments to verify the performance of sequence-structure alignment. As revealed in these experiments in Section 3.2 and Section 3.3 our FM achieves up to an macro-F1 score in zero-shot Seq2Str prediction, i.e., secondary structure prediction, outperforming some fine-tuned FMs and ViennaRNA. In terms of Str2Seq prediction performance, we leverage our FM to solve the complex RNA design task and solved of complex puzzles of the Eterna [17] V2 benchmark, while state-of-the-art FMs, e.g., SpliceBERT [21], only solved of simple (i.e., the target sequences are very short) puzzles. Besides, our FM takes at most hours to solve most of the puzzles, while the majority of RNA design methods can take up to hours to solve each puzzle.
Open-source Resources. The genomic sequence (e.g., RNA and DNA) modeling is still on its early stage, while the related open-source resources are extremely limited and undocumented. Besides, many of the early works on genomic FMs are not ready-to-use due to different reasons. This problem encourages us to release all the benchmarks, evaluation scripts and FM tutorials, etc. As a result, we have developed an open-source package 222https://github.com/yangheng95/OmniGenomeBench that includes step-by-step tutorials for FM pretraining and downstream tasks fine-tuning, etc. It provides prepared genomic benchmarks and automates the benchmarking process of FMs using the API with only a few lines of code. These resources avoid the tedious and time-consuming programming in future genomic FM research.
Finally, the FM proposed in this work is named OmniGenome, which may advance a deep mechanistic understanding of RNA sequences as well as molecular biology and accelerate our ability to design RNA molecules.
2 Methodology
This section delineates the implementation of OmniGenome step by step.
2.1 Data Preprocessing
![Refer to caption](x3.png)
Recent studies [21, 22] have shown that data diversity can enhance FM performance without significantly increasing model capacity. For the OmniGenome pretraining, we collected transcriptome data from the OneKP initiative [23], which compiles large-scale RNA-Seq data from plant species. Because the raw biological data are sparse, noisy, and characterized by long sequences, they are far from ready for effective FM training.
To address this problem, we developed a four-step data curation protocol to improve data quality. Raw RNA-Seq data are often excessively long with thousands of bases. We first sliced them into segments with a window size of bases to provide a sufficient context window for RNA sequence understanding. To enhance training efficiency and reduce bias, we removed all duplicate sequences. To tackle incomplete transcriptome data and other noises, we discard sequences shorter than bases. To facilitate the sequence-structure alignment training, we adopt ViennaRNA to obtain the secondary structures for the sequences.
2.2 Tokenization
In the realm of genomic sequence modeling, modeling resolution, and performance highly depend on the tokenization strategy [3, 21]. Because the traditional k-mers [24, 13] and BPE [19, 22] methods combine multiple bases into single tokens, they may compromise modeling resolution and thus fail to capture granular features. To achieve a single nucleotide (SN) resolution necessary for genomic tasks like RNA design and structure prediction, we employed single nucleotide tokenization (SNT) method [3, 21] (see an illustrative example in Fig. 3). Additionally, to simplify the implementation, we adopted a vocabulary ‘A’, ‘T’, ‘C’, ‘G’, ‘U’, ‘N’, ‘(’, ‘)’, ‘.’ to unify the tokenization of both nucleotide bases and RNA secondary structure information. Following BERT-style models [19], we incorporated special tokens, e.g., ‘<mask>’, to facilitate masked bases and structure prediction.
![Refer to caption](x4.png)
2.3 Pretraining Objectives
As discussed in Section 1, a key desideratum for complex genomic understanding and RNA sequence design is helping the model build an alignment between RNA sequences and their corresponding secondary structures. Bearing this in mind, we designed the following aggregated pretraining objective for OmniGenome:
(1) |
where is the regularization weight and represents the parameters of OmniGenome. The following paragraphs explain the design principles of each loss function used in equation (1).
-
•
is designed to enable OmniGenome to predict bases given partially masked sequences with the context of structures, so that OmniGenome can perceive RNA secondary structure by constructing Str2Seq prediction. To achieve this, we mask of the bases and structure tokens, encouraging the model to infer masked bases (i.e., ‘A’, ‘T’, ‘C’, ‘G’, ‘U’, ‘N’) and structure tokens (i.e., ‘(’, ‘)’, ‘.’). Specifically, is defined as the classic cross-entropy loss widely used in the masked language modeling:
(2) where is the number of masked nucleotide and structure tokens, and indicates the probability of predicting the masked nucleotide based on its context.
-
•
Complementary to , is designed to enable OmniGenome for Seq2Str prediction. Instead of directly feeding the secondary structure into OmniGenome during training, this objective employs the RNA secondary structures as supervised-training labels. Unlike Str2Seq, the Seq2Str prediction is generally deterministic because the optimal structures are unique. This task is implemented as a token-level classification, where the loss is defined in the following cross-entropy loss:
(3) where denotes the label of secondary structure at the -th position, and is the probability predicted by a linear classifier deployed on OmniGenome. is the length of an RNA sequence and denotes the number of the possible labels of structure, i.e., ‘(’, ‘)’, ‘.’.
-
•
The last one is similar to the conventional masked language modeling loss in NLP. It aims to improve the model’s understanding of latent patterns in RNA sequences by masking or randomly replacing of nucleotide bases. The definition of is similar to that of , but it only considers the prediction of masked nucleotide bases.
There will be a data leakage problem in objective while the structure is appended to the sequences. In practice, of RNA sequences are used for structure-contextualized pretraining for the objective, while the remaining are used for the other two objectives. This proportion is based on our experiment observation.
2.4 Model Architecture
OmniGenome adopts the classic Transformer encoder architecture with bidirectional multi-head attention. We designed two variants, dubbed OmniGenome and OmniGenome with and million parameters respectively. Some key model specifications are summarized in Table 1.
The training specifications are summarized as follows. The learning rate is set to and the weight decay is set to . We use AdamW as the optimizer with hyperparameters and . We use a linear decay strategy with a warm-up period of steps in the learning rate scheduler. The batch size is set to . No dropout is applied during pretraining, and we use the rotary position embeddings [25] to further enhance the model’s scalability to long RNA sequences. We built a distributed training environment with Nvidia RTX GPUs and pretrained for weeks, while its configuration is introduced in Appendix 6.
Models | # of | Embedding | Intermediate | # of | # of | Modeling | Position |
---|---|---|---|---|---|---|---|
Layers | dimension | dimension | heads | parameters | length | embedding | |
OmniGenome | M | Rotary | |||||
OmniGenome | M | Rotary |
2.5 Downstream Tasks
2.5.1 Predictive Tasks Fine-tuning
OmniGenome is designed as a general-purpose RNA FM that can be fine-tuned for a diverse set of downstream genomics predictive tasks. In this paper, we constructed a large-scale benchmark suite for RNA FMs. According to the category of genomes, we split the benchmark into two parts.
-
•
The first one is a collection of genomic understanding tasks aimed at the RNA genomic benchmark (RGB), as shown in Table 7. RGB contains SN-level tasks that are curated in this work or collected from published articles. The purpose of RGB is to benchmark genomic FMs in challenging SN-level modeling tasks like detection and repair of SN mutations, mRNA sequence degradation rates, and RNA secondary structure prediction. Due to the lack of a plant RNA benchmark dataset, the RGB considers the modeling of RNA sequences from a variety of species, e.g., plant and human. The sequence length in RGB ranges from to , which is enough for most RNA understanding tasks. In summary, these multi-species and SN-level tasks in RGB serve as the first comprehensive benchmark utilized to assess the RNA sequence modeling capabilities of OmniGenome and its baseline models. For detailed information on each dataset, such as their sources and sizes, please refer to Appendix 10.1.
-
•
The plant genomic benchmark [14] (PGB333https://huggingface.co/datasets/InstaDeepAI/plant-genomic-benchmark) shown in Table 9 provides a comprehensive suite of DNA genomic tasks designed to evaluate and improve the predictive capabilities of genomic models in plant biology. As a large-scale plant DNA FM benchmark, PGB involves types of DNA modeling subtasks, including a range of critical genomic tasks such as promoter strength prediction and gene expression regression. There are datasets in total with millions of DNA sequences to be evaluated in PGB, and the sequence lengths are up to , which is very long for most of the genomic FMs. Since the original evaluation protocol is not publicly available, we have re-implemented the auto-benchmark for all the subtasks from PGB in our package, including various aspects of plant genomics such as gene expression and chromatin accessibility. By integrating diverse genomic tasks, the PGB aims to facilitate advanced research and development in plant genomics, offering a robust platform for the assessment and enhancement of model performance across different plant species. Due to computational limitations, we randomly sample a maximum of k examples in all datasets in PGB to evaluate the FM’s performance.
2.5.2 RNA Sequence Design
One of the difficult practices addressed by OmniGenome is the RNA design task, which has been long-term neglected by existing FMs because of the incapability of Str2Seq prediction. We introduce a naive genetic algorithm based on OmniGenome’s Str2Seq capability and achieve state-of-the-art performance in RNA sequence design by precisely predicting bases according to the structure context. The main steps in the genetic algorithm and workflow visualization are available in Appendix 9 and Fig. 5, respectively. Please find the parameter setting details in Section 3.3, and the experiment script can be found in the supplemental materials.
3 Experiments
To comprehensively explore the potential of OmniGenome and assess its performance across genomic understanding and RNA design tasks, this section implements rich experimental evaluations on diverse genomic tasks. We first evaluate the sequence-structure alignment capability of OmniGenome. Next, we evaluate the overall performance of OmniGenome on two comprehensive genomic modeling benchmarks, i.e., the RNA Genomic Benchmark (RGB) and the Plant Genomic Benchmark (PGB), respectively. All the benchmark datasets are not from the pretraining database, i.e., OneKP.
3.1 Evaluation Methods
There are no direct counterparts to OmniGenome in plant genomic modeling. However, we can compare it with the following recent genomic FMs as potential baselines to help evaluate the performance of OmniGenome as shown in Table 2. Please refer to Appendix 8 for brief introductions of these FMs. We are aware that some FMs are also developed for RNA, such as RNA-FM [26], RNA-MSM [27], Uni-RNA [28], 5UTR-LM [16], etc. However, we cannot compare OmniGenome with them because their source codes are very hard to work with in our efforts or the models are not publicly available. In the secondary structure prediction tasks, we input the secondary structures computed by ViennaRNA that are different from the true structures into OmniGenome to avoid data leakage problems.
Model | Tokenization | # of Params | Pretraining Data | Species | Nucleic Acid |
---|---|---|---|---|---|
DNABERT-2 | BPE | M | B Tokens | Human + Species | DNA |
NT-V2-M | k-mers | M | B Tokens | Human + Species | DNA |
HyenaDNA-Large | SNT | M | B Tokens | Human | DNA |
Agro-NT-B | k-mers | M | B Tokens | Edible Plants | DNA |
SpliceBERT | SNT | M | M Sequences | Multi-Vertebrates | precursor-mRNA |
CDSBERT | SNT | M | M Sequences | 9 RNA Families | CDS |
3UTRBERT | k-mers | M | Sequences | Multi-Species | mRNA ’UTR |
OmniGenome | SNT | M | B Tokens | Plant Species | mRNA, CDS, UTR |
OmniGenome | M |
3.2 Zero-shot Structure to Sequence Prediction Evaluation
In this experiment, we demonstrate the Str2Seq prediction capability of OmniGenome based on RNA sequence design. We employed the Eterna [17] V2 benchmark, which consists of specified secondary structures. This task aims to design plausible RNA sequences based on the reference structures. It is intuitive that higher performance in RNA design tasks indicates better structure-to-sequence prediction capability. We develop a simple GA algorithm based on OmniGenome models to solve RNA design puzzles. The implementation details can be found in Fig. 5 in Appendix 9. Within our GA implementation, the population size is set at , with iterations, and the mutation rate for each base is . We set the maximum running time for each puzzle to hour because OmniGenome can easily solve most of the RNA design puzzles within hour, while each puzzle’s time budget is generally hours in existing RNA design methods. We utilize accuracy as our evaluation metric, meaning that only sequences with entirely correct structures are considered resolved. The experimental results are available in Table 3.
Model | Token. | EternaV2 |
Acc | ||
RNAInverse | — | |
SpliceBERT | SNT | |
CDSBERT | SNT | |
OmniGenome | SNT | |
OmniGenome | SNT | |
OmniGenome | SNT | |
OmniGenome | SNT |
It is observed from Table 3 that the FM-free method RNAInverse solved 30 of the RNA design puzzles, indicating a competent but not optimal capability in sequence design. We select the FMs employing SNT to conduct comparison experiments, as tokenizers like k-mers and BPE cannot handle the SN-level resolution modeling, i.e., masked nucleotide prediction. We cannot compare with HyenaDNA because it does not support masked nucleotide prediction. The models SpliceBERT and CDSBERT demonstrated trivial proficiency in the RNA sequence design task, solving and puzzles, respectively. This suggests these FMs cannot precisely predict the bases without any Str2Seq mapping ability. OmniGenome variants without structure context input have performance similar to the baselines, i.e., OmniGenome and OmniGenome solved and puzzles, respectively. In contrast, OmniGenome and OmniGenome significantly outperformed other models with and puzzles solved, respectively, underscoring the impact of Str2Seq prediction in sequence-structure alignment, which evidently bolsters the models’ abilities to predict bases according to the structure contexts. Besides, we expect an increase in performance with sufficient computational budgets and the findings provide crucial evidence of the significance of Str2Seq for RNA sequence design.
3.3 Zero-shot Sequence to Structure Prediction Evaluation
This subsection assesses OmniGenome in both Seq2Str and Str2Seq prediction aspects of sequence-structure alignment. The evaluation of Seq2Str is based on zero-shot secondary structure prediction. We use OmniGenome and OmniGenome without fine-tuning to predict the secondary structures of sequences from the test datasets and measure the macro-F1 score, where better structure prediction performance indicates a stronger capability for Seq2Str prediction. We evaluate Str2Seq mapping ability based on masked nucleotide base prediction. This task leverages the FM’s capacity to predict sequences with the context of structure. We randomly mask one base in the sequences from the test datasets and leverage the Str2Seq ability of OmniGenome to predict the masked bases in a zero-shot scenario. We assess the model’s accuracy in predicting masked bases both with and without structure context. The experimental results are available in Table 4.
Model | RNA Secondary Structure Prediction | Masked Nucleotide Prediction | ||||||
---|---|---|---|---|---|---|---|---|
Archive2 | bpRNA | RNAStralign | Rfam | Archive2 | bpRNA | RNAStralign | Rfam | |
F1 | F1 | F1 | F1 | Acc | Acc | Acc | Acc | |
ViennaRNA | — | — | — | — | ||||
SpliceBERT | — | — | — | — | ||||
OmniGenome | ||||||||
OmniGenome | ||||||||
OmniGenome | ||||||||
OmniGenome |
The results in Table 4 reveal OmniGenome’s effectiveness in both secondary structure prediction (i.e., Seq2Str) and masked nucleotide prediction (i.e., Str2Seq). More specifically, OmniGenome and OmniGenome show impressive performance for zero-shot sequence-structure alignment, while OmniGenome and OmniGenome outperform the base variants of OmniGenome. Even when ablating structure context, OmniGenome achieves performance comparable with ViennaRNA. Interestingly, if the structures calculated by ViennaRNA are fed into the model, OmniGenome models can predict the structure more precisely and outperform ViennaRNA on most of the secondary structure prediction datasets, e.g., bpRNA, RNAStralign, and Rfam. From the results of masked nucleotide prediction, it is observed that sequence reconstruction performance can be enhanced by the structure context. OmniGenome and OmniGenome significantly outperform baseline FMs like SpliceBERT.
Overall, the results in Table 4 provide a comprehensive evaluation of the FMs’ capabilities in sequence-structure alignment. This underscores OmniGenome’s performance and robustness in generalizing among different species without prior fine-tuning, highlighting their potential utility in SN-level in-silico genomic modeling.
3.4 RNA Genomic Benchmark
Model | mRNA | SNMD | SNMR | Archive2 | bpRNA | RNAStralign |
---|---|---|---|---|---|---|
RMSE | AUC | F1 | F1 | F1 | F1 | |
ViennaRNA | N.A. | N.A. | N.A. | 75.89 | 27.82 | 74.80 |
DNABERT2 | ||||||
HyenaDNA | ||||||
NT-V2 | ||||||
Agro-NT | ||||||
SpliceBERT | ||||||
3UTRBERT | ||||||
CDSBERT | ||||||
OmniGenome | ||||||
OmniGenome | ||||||
OmniGenome | ||||||
OmniGenome |
The results in Table 5 demonstrate the performance of OmniGenome and its generalizability across various fine-grained RNA genomic modeling tasks. It is observed that the base variants of OmniGenome models achieve better results than both RNA and DNA FM baselines, including Agro-NT and CDSBERT, which contain hundreds of millions of parameters. This is because the existing FMs usually adopt k-mers tokenization that cannot handle single nucleotide resolution tasks, e.g., single nucleotide mutation detection and repair. Because of the Seq2Str pretraining, OmniGenome and OmniGenome models exhibit strong results in secondary structure prediction, underscoring OmniGenome’s capabilities in SN-level RNA sequence understanding and manipulation.
3.5 Plant Genomic Benchmark
Model | PolyA | LncRNA | Chrom Acc | Prom Str | Term Str | Splice | Gene Exp | Enhancer |
---|---|---|---|---|---|---|---|---|
F1 | F1 | F1 | RMSE | RMSE | F1 | RMSE | F1 | |
DNABERT2 | ||||||||
HyenaDNA | ||||||||
NT-V2 | ||||||||
Agro-NT | ||||||||
SpliceBERT | ||||||||
3UTRBERT | ||||||||
CDSBERT | ||||||||
OmniGenome | ||||||||
OmniGenome | ||||||||
OmniGenome |
The PGB is a plant-oriented genomic benchmark. Although the benchmark datasets in PGB are DNA-based tasks, we can still evaluate the performance of OmniGenome and its generalizability on multi-modal (i.e., DNA and RNA) genomic tasks. The results in Table 6 reveal substantial variability in the performance of different FMs, where OmniGenome outperformed other baseline models across most tasks, particularly in tasks like Polyadenylation, Splice Site, and Enhancer Region classification, where they achieved the highest F1 scores. This suggests that OmniGenome’s architecture is particularly adept at handling complex genomic sequences. In comparison, existing FMs, e.g., CDSBERT and Agro-NT, showed lower performance with more parameters than OmniGenome. Besides, the performance of OmniGenome suggests that the structure context can further enhance the performance of genomic modeling.
Overall, OmniGenome models achieve state-of-the-art performance on both benchmarks, especially for OmniGenome variants. The results underscore the importance of sequence-structure alignment in achieving complex genomic modeling tasks.
4 Related Works
Biological sequence modeling, including DNA, RNA, and protein, has attracted attention in recent years. Protein modeling, e.g., AlphaFold [29, 30, 31] and ESM [32], has been studied for many years compared to DNA and RNA modeling. In the realm of genomic sequence modeling, several early works aimed at addressing diversified genome downstream subtasks. For instance, DNABERT [33] adapts the architecture of BERT [19] for genomic sequence modeling, showing preliminary performance for in-silico genomic tasks. DNABERT2 [22], a multi-species FM improved based on DNABERT, proposes replacing k-mers tokenization with BPE tokenization to improve model performance. To explore the performance of large-scale FMs, the nucleotide transformers (V1 & V2) [13], AgroNT [14] and SegmentNT [34] leveraged billions of parameters to boost genomic sequence modeling and achieved promising performance in understanding DNA genome, with model scales up to billion and billion parameters, respectively. Agro-NT [14] was pretrained on multi-species edible plant DNA sequences but failed to transfer effectively to RNA sequence modeling in our experiments. To address the modeling capacity problem caused by the remarkable lengths of genomes, there is growing focus on the necessity of long-range sequence modeling and the introduction of autoregressive FMs, namely, HyenaDNA [3] and Evo [2]. In terms of RNA genomic modeling, there have been some preliminary works, such as scBERT [35], RNABERT [36], RNA-FM [26], RNA-MSM [37], and RNAErnie [38], to name a few. However, these methods have only trained the FMs on a limited-scale database, as RNA sequences are generally expensive to obtain. Some FMs focus on specific types of RNA sequences, such as coding sequences (CDS) [39], 5’ untranslated regions (5’UTR) [40], 3’ untranslated regions (3’UTR) [24], or precursor mRNA sequences [21], thus limiting the models’ ability to capture the diversity of RNA sequences. Uni-RNA [28] has been reported to achieve good performance due to the large scale of the model and database, however, it is not open-sourced and cannot be compared in the experiments.
In short, the existing FMs usually neglect the sequence-structure alignment problem in RNA genomic modeling, while the 5UTR-LM [16] adopts the secondary structure prediction as a pre-training objective to achieve Seq2Str prediction in pretraining. However, these FMs are not available for Str2Seq mapping and suffer from limited model and data scales that fail to uncover the comprehensive efficacy of sequence-structure alignment on a wide set of genomic tasks. ERNIE-RNA [41] feeds the RNA structure along with the sequence into the model and improves the downstream tasks. However, it also ignores the significance of Str2Seq prediction capability. In a nutshell, existing FMs fail to achieve sequence-structure alignment without exception.
5 Conclusion
In this paper, we introduced OmniGenome, a foundation model (FM) that addresses the critical challenge of sequence-structure alignment in genomic modeling. Our approach leverages a structure-contextualized modeling to achieve accurate predictions for both Seq2Str and Str2Seq tasks. By utilizing the extensive OneKP dataset, we ensured the robustness and transferability of OmniGenome across various genomic tasks. Our experimental results demonstrate that OmniGenome significantly outperforms existing models in zero-shot secondary structure prediction, achieving up to macro-F1 score, and excels in solving complex RNA design puzzles with remarkable efficiency. These findings validate the efficacy of our approach and highlight the potential of OmniGenome in advancing RNA and broader genomic research.
Acknowledgment
This work was supported in part by the UKRI Future Leaders Fellowship under Grant MR/S017062/1 and MR/X011135/1; in part by NSFC under Grant 62376056 and 62076056; in part by the Royal Society under Grant IES/R2/212077; in part by the EPSRC under Grant 2404317; in part by the Kan Tong Po Fellowship (KTP\R1\231017); and in part by the Amazon Research Award and Alan Turing Fellowship.
References
- [1] N. Wang, J. Bian, Y. Li, X. Li, S. Mumtaz, L. Kong, and H. Xiong, “Multi-purpose rna language modelling with motif-aware pretraining and type-guided fine-tuning,” Nature Machine Intelligence, pp. 1–10, 2024.
- [2] E. Nguyen, M. Poli, M. G. Durrant, A. W. Thomas, B. Kang, J. Sullivan, M. Y. Ng, A. Lewis, A. Patel, A. Lou et al., “Sequence modeling and design from molecular to genome scale with evo,” bioRxiv, pp. 2024–02, 2024.
- [3] E. Nguyen, M. Poli, M. Faizi, A. W. Thomas, C. Birch-Sykes, M. Wornow, A. Patel, C. M. Rabideau, S. Massaroli, Y. Bengio, S. Ermon, S. A. Baccus, and C. Ré, “Hyenadna: Long-range genomic sequence modeling at single nucleotide resolution,” CoRR, vol. abs/2306.15794, 2023. [Online]. Available: https://doi.org/10.48550/arXiv.2306.15794
- [4] Z. Tan, Y. Fu, G. Sharma, and D. H. Mathews, “Turbofold ii: Rna structural alignment and secondary structure prediction informed by multiple homologs,” Nucleic acids research, vol. 45, no. 20, pp. 11 570–11 581, 2017.
- [5] P. Danaee, M. Rouches, M. Wiley, D. Deng, L. Huang, and D. Hendrix, “bprna: large-scale automated annotation and analysis of rna secondary structure,” Nucleic acids research, vol. 46, no. 11, pp. 5381–5394, 2018.
- [6] D. H. Mathews, “How to benchmark rna secondary structure prediction accuracy,” Methods, vol. 162, pp. 60–67, 2019.
- [7] I. Kalvari, E. P. Nawrocki, N. Ontiveros-Palacios, J. Argasinska, K. Lamkiewicz, M. Marz, S. Griffiths-Jones, C. Toffano-Nioche, D. Gautheret, Z. Weinberg et al., “Rfam 14: expanded coverage of metagenomic, viral and microrna families,” Nucleic Acids Research, vol. 49, no. D1, pp. D192–D200, 2021.
- [8] O. Yaish and Y. Orenstein, “Computational modeling of mrna degradation dynamics using deep neural networks,” Bioinformatics, vol. 38, no. 4, pp. 1087–1101, 2022.
- [9] H. K. Wayment-Steele, W. Kladwang, A. M. Watkins, D. S. Kim, B. Tunguz, W. Reade, M. Demkin, J. Romano, R. Wellington-Oguri, J. J. Nicol et al., “Deep learning models for predicting rna degradation via dual crowdsourcing,” Nature Machine Intelligence, vol. 4, no. 12, pp. 1174–1184, 2022.
- [10] K. S. Corbett, D. K. Edwards, S. R. Leist, O. M. Abiona, S. Boyoglu-Barnum, R. A. Gillespie, S. Himansu, A. Schäfer, C. T. Ziwawo, A. T. DiPiazza et al., “Sars-cov-2 mrna vaccine design enabled by prototype pathogen preparedness,” Nature, vol. 586, no. 7830, pp. 567–571, 2020.
- [11] F. Runge, J. K. Franke, D. Fertmann, R. Backofen, and F. Hutter, “Partial rna design,” bioRxiv, 2023. [Online]. Available: https://www.biorxiv.org/content/early/2023/12/30/2023.12.29.573656
- [12] Y. Wan, M. Kertesz, R. C. Spitale, E. Segal, and H. Y. Chang, “Understanding the transcriptome through rna structure,” Nature Reviews Genetics, vol. 12, no. 9, pp. 641–655, 2011.
- [13] H. Dalla-Torre, L. Gonzalez, J. Mendoza-Revilla, N. L. Carranza, A. H. Grzywaczewski, F. Oteri, C. Dallago, E. Trop, B. P. de Almeida, H. Sirelkhatim et al., “The nucleotide transformer: Building and evaluating robust foundation models for human genomics,” bioRxiv, pp. 2023–01, 2023.
- [14] J. Mendoza-Revilla, E. Trop, L. Gonzalez, M. Roller, H. Dalla-Torre, B. P. de Almeida, G. Richard, J. Caton, N. Lopez Carranza, M. Skwark et al., “A foundational large language model for edible plant genomes,” bioRxiv, pp. 2023–10, 2023.
- [15] B. Frey, A. Celaj, A. Gao, T. Lau, E. Holgersen, A. Lo, V. Lodaya, C. Cole, R. Denroche, C. Spickett et al., “An rna foundation model enables discovery of disease mechanisms and candidate therapeutics,” bioRxiv, pp. 2023–09, 2023.
- [16] Y. Chu, D. Yu, Y. Li, K. Huang, Y. Shen, L. Cong, J. Zhang, and M. Wang, “A 5’ utr language model for decoding untranslated regions of mrna and function predictions,” Nature Machine Intelligence, pp. 1–12, 2024.
- [17] J. Lee, W. Kladwang, M. Lee, D. Cantu, M. Azizyan, H. Kim, A. Limpaecher, S. Gaikwad, S. Yoon, A. Treuille et al., “Rna design rules from a massive open laboratory,” Proceedings of the National Academy of Sciences, vol. 111, no. 6, pp. 2122–2127, 2014.
- [18] I. Tinoco Jr and C. Bustamante, “How rna folds,” Journal of molecular biology, vol. 293, no. 2, pp. 271–281, 1999.
- [19] J. Devlin, M. Chang, K. Lee, and K. Toutanova, “BERT: pre-training of deep bidirectional transformers for language understanding,” in NAACL-HLT (1). Association for Computational Linguistics, 2019, pp. 4171–4186.
- [20] R. Lorenz, S. H. Bernhart, C. Höner zu Siederdissen, H. Tafer, C. Flamm, P. F. Stadler, and I. L. Hofacker, “Viennarna package 2.0,” Algorithms for molecular biology, vol. 6, pp. 1–14, 2011.
- [21] K. Chen, Y. Zhou, M. Ding, Y. Wang, Z. Ren, and Y. Yang, “Self-supervised learning on millions of pre-mrna sequences improves sequence-based rna splicing prediction,” bioRxiv, pp. 2023–01, 2023.
- [22] Z. Zhou, Y. Ji, W. Li, P. Dutta, R. V. Davuluri, and H. Liu, “DNABERT-2: efficient foundation model and benchmark for multi-species genome,” CoRR, vol. abs/2306.15006, 2023. [Online]. Available: https://doi.org/10.48550/arXiv.2306.15006
- [23] E. J. Carpenter, J. H. Leebens-Mack, and M. S. B. et al., “One thousand plant transcriptomes and the phylogenomics of green plants,” Nature, vol. 574, no. 7780, pp. 679–685, 2019.
- [24] Y. Yang, G. Li, K. Pang, W. Cao, X. Li, and Z. Zhang, “Deciphering 3’utr mediated gene regulation using interpretable deep representation learning,” bioRxiv, pp. 2023–09, 2023.
- [25] J. Su, M. H. M. Ahmed, Y. Lu, S. Pan, W. Bo, and Y. Liu, “Roformer: Enhanced transformer with rotary position embedding,” Neurocomputing, vol. 568, p. 127063, 2024.
- [26] J. Chen, Z. Hu, S. Sun, Q. Tan, Y. Wang, Q. Yu, L. Zong, L. Hong, J. Xiao, T. Shen et al., “Interpretable rna foundation model from unannotated data for highly accurate rna structure and function predictions,” bioRxiv, pp. 2022–08, 2022.
- [27] Y. Zhang, M. Lang, J. Jiang, Z. Gao, F. Xu, T. Litfin, K. Chen, J. Singh, X. Huang, G. Song et al., “Multiple sequence alignment-based rna language model and its application to structural inference,” Nucleic Acids Research, vol. 52, no. 1, pp. e3–e3, 2024.
- [28] X. Wang, R. Gu, Z. Chen, Y. Li, X. Ji, G. Ke, and H. Wen, “Uni-rna: universal pre-trained models revolutionize rna research,” bioRxiv, pp. 2023–07, 2023.
- [29] J. Jumper, R. Evans, A. Pritzel, T. Green, M. Figurnov, O. Ronneberger, K. Tunyasuvunakool, R. Bates, A. Žídek, A. Potapenko, A. Bridgland, C. Meyer, S. A. A. Kohl, A. J. Ballard, A. Cowie, B. Romera-Paredes, S. Nikolov, R. Jain, J. Adler, T. Back, S. Petersen, D. Reiman, E. Clancy, M. Zielinski, M. Steinegger, M. Pacholska, T. Berghammer, S. Bodenstein, D. Silver, O. Vinyals, A. W. Senior, K. Kavukcuoglu, P. Kohli, and D. Hassabis, “Highly accurate protein structure prediction with AlphaFold,” Nature, vol. 596, no. 7873, pp. 583–589, 2021.
- [30] R. Evans, M. O’Neill, A. Pritzel, N. Antropova, A. Senior, T. Green, A. Žídek, R. Bates, S. Blackwell, J. Yim, O. Ronneberger, S. Bodenstein, M. Zielinski, A. Bridgland, A. Potapenko, A. Cowie, K. Tunyasuvunakool, R. Jain, E. Clancy, P. Kohli, J. Jumper, and D. Hassabis, “Protein complex prediction with alphafold-multimer,” bioRxiv, 2021. [Online]. Available: https://www.biorxiv.org/content/early/2021/10/04/2021.10.04.463034
- [31] J. Abramson, J. Adler, J. Dunger, R. Evans, T. Green, A. Pritzel, O. Ronneberger, L. Willmore, A. J. Ballard, J. Bambrick et al., “Accurate structure prediction of biomolecular interactions with alphafold 3,” Nature, pp. 1–3, 2024.
- [32] Z. Lin, H. Akin, R. Rao, B. Hie, Z. Zhu, W. Lu, A. dos Santos Costa, M. Fazel-Zarandi, T. Sercu, S. Candido et al., “Language models of protein sequences at the scale of evolution enable accurate structure prediction,” BioRxiv, vol. 2022, p. 500902, 2022.
- [33] Y. Ji, Z. Zhou, H. Liu, and R. V. Davuluri, “DNABERT: pre-trained bidirectional encoder representations from transformers model for dna-language in genome,” Bioinform., vol. 37, no. 15, pp. 2112–2120, 2021.
- [34] B. P. de Almeida, H. Dalla-Torre, G. Richard, C. Blum, L. Hexemer, M. Gélard, J. Mendoza-Revilla, P. Pandey, S. Laurent, M. Lopez et al., “Segmentnt: annotating the genome at single-nucleotide resolution with dna foundation models,” bioRxiv, pp. 2024–03, 2024.
- [35] F. Yang, W. Wang, F. Wang, Y. Fang, D. Tang, J. Huang, H. Lu, and J. Yao, “scbert as a large-scale pretrained deep language model for cell type annotation of single-cell rna-seq data,” Nat. Mac. Intell., vol. 4, no. 10, pp. 852–866, 2022. [Online]. Available: https://doi.org/10.1038/s42256-022-00534-z
- [36] M. Akiyama and Y. Sakakibara, “Informative rna base embedding for rna structural alignment and clustering by deep representation learning,” NAR genomics and bioinformatics, vol. 4, no. 1, p. lqac012, 2022.
- [37] Y. Zhang, F. Ge, F. Li, X. Yang, J. Song, and D.-J. Yu, “Prediction of multiple types of rna modifications via biological language model,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, 2023.
- [38] N. Wang, J. Bian, Y. Li, X. Li, S. Mumtaz, L. Kong, and H. Xiong, “Multi-purpose rna language modelling with motif-aware pretraining and type-guided fine-tuning,” Nature Machine Intelligence, pp. 1–10, 2024.
- [39] L. Hallee, N. Rafailidis, and J. P. Gleghorn, “cdsbert-extending protein language models with codon awareness,” bioRxiv, 2023.
- [40] Y. Chu, D. Yu, Y. Li, K. Huang, Y. Shen, L. Cong, J. Zhang, and M. Wang, “A 5’ utr language model for decoding untranslated regions of mrna and function predictions,” Nature Machine Intelligence, pp. 1–12, 2024.
- [41] W. Yin, Z. Zhang, L. He, R. Jiang, S. Zhang, G. Liu, X. Zhang, T. Qin, and Z. Xie, “Ernie-rna: An rna language model with structure-enhanced representations,” bioRxiv, pp. 2024–03, 2024.
- [42] C. Peng, J. Shang, J. Guan, D. Wang, and Y. Sun, “Viralm: Empowering virus discovery through the genome foundation model,” bioRxiv, pp. 2024–01, 2024.
- [43] J. Kaplan, S. McCandlish, T. Henighan, T. B. Brown, B. Chess, R. Child, S. Gray, A. Radford, J. Wu, and D. Amodei, “Scaling laws for neural language models,” CoRR, vol. abs/2001.08361, 2020. [Online]. Available: https://arxiv.org/abs/2001.08361
- [44] J. Hoffmann, S. Borgeaud, A. Mensch, E. Buchatskaya, T. Cai, E. Rutherford, D. de Las Casas, L. A. Hendricks, J. Welbl, A. Clark, T. Hennigan, E. Noland, K. Millican, G. van den Driessche, B. Damoc, A. Guy, S. Osindero, K. Simonyan, E. Elsen, J. W. Rae, O. Vinyals, and L. Sifre, “Training compute-optimal large language models,” CoRR, vol. abs/2203.15556, 2022. [Online]. Available: https://doi.org/10.48550/arXiv.2203.15556
- [45] N. Muennighoff, A. M. Rush, B. Barak, T. L. Scao, A. Piktus, N. Tazi, S. Pyysalo, T. Wolf, and C. Raffel, “Scaling data-constrained language models,” CoRR, vol. abs/2305.16264, 2023. [Online]. Available: https://doi.org/10.48550/arXiv.2305.16264
6 Pretraining Environment
The pretraining of OmniGenome was conducted on a dedicated Linux computation node, equipped with Nvidia RTX GPUs. For distributed model training, we employed version of the Transformers library alongside version of the Accelerate library. Our implementation framework of choice for OmniGenome was PyTorch, specifically version . The ViennaRNA version is in our experiments. While some existing code was adapted for the modules within OmniGenome, the majority of the codebase, such as genomic sequences preprocessing, model pretraining, objective functions, and experiments, was meticulously crafted from scratch.
7 OmniGenome Package
Genomic modeling is still in its early stages, and resources are consequently very scarce. Typically, existing foundation models open-source only the model, without providing the training, fine-tuning, and benchmark evaluation codes, etc. To address this issue, we have developed a comprehensive open-source genomic modeling toolkit444https://github.com/yangheng95/OmniGenomeBench based on OmniGenome. This toolkit aims to provide extensive FM fine-tuning tutorials and a unified automated benchmark evaluation. The main features of the OmniGenome Package are as follows:
-
•
Fine-Tuning Tutorials: We provide tutorials for fine-tuning on all downstream genomic modeling tasks, including dataset processing, model implementation, and training processes. A fine-tuning example for secondary structure is included, covering both training and demonstration of secondary structure prediction. The fine-tuning tutorials are available at: https://github.com/yangheng95/OmniGenomeBench/blob/master/examples/secondary_structure_prediction_demo.ipynb
-
•
Automated Benchmark Evaluation: We offer an automated benchmark evaluation interface, which includes the built-in PGB and RGB benchmarks. By predefining the configurations for benchmark evaluation subtasks, such as hyperparameters, our tool supports the automated benchmark evaluation of future FMs and the addition of new benchmarks. The goal of automated benchmark evaluation is to ensure fairness and ease of use. We provide a tutorial on automated evaluation to guide users in benchmark evaluation. The automated benchmarking example is available at: https://github.com/yangheng95/OmniGenomeBench/blob/master/examples/benchmark/batch_rgb_benchmark.py
-
•
Genomic Repository Hub: We have created a hub for hosting and distributing open-source licensed datasets, model checkpoints, and benchmark evaluations. Additionally, we have designed flexible interfaces to support the sharing of datasets and models within the community. This approach helps mitigate the issue of resource scarcity. The hub will be available soon.
We are in the process of finalizing the necessary documentation and will officially release this tool in the near future.
8 Evaluation Methods
There are no direct counterparts to OmniGenome in plant genomic modeling. However, we can compare it with the following recent genomic FMs as potential baselines to help evaluate the performance of OmniGenome. The brief introductions of the FMs in Table 2 are as follows:
-
•
ViennaRNA [20]. ViennaRNA is a comprehensive genomic analysis tool that includes a diverse set of interfaces, such as RNAFold and RNAInverse design. ViennaRNA serves as the baseline for RNA structure prediction and RNA design in our experiments.
-
•
DNABERT2 [22]. DNABERT2 is one of the latest DNA FMs which improves the performance of DNABERT. The main modification of DNABERT2 is the tokenization method, which was changed to BPE from k-mers.
-
•
HyenaDNA [3]. HyenaDNA is an autoregressive FM optimized for long-range genomic data processing. HyenaDNA is based on the Hyena convolution architecture and capable of handling sequences up to M bases in length.
-
•
Nucleotide Transformer (NT) V2 [13]. The NT FMs were trained on DNA data, including the human reference genome and multi-species DNA sequences. They aim to capture the complex patterns within nucleotide sequences for various genomic applications.
-
•
Agricultural Nucleotide Transformer (Agro-NT) [14]. Agro-NT is a large-scale DNA FM (B parameters) akin to the Nucleotide Transformers but with a focus on plant DNA.
-
•
SpliceBERT [21]. It was trained on M precursor messenger RNA (pre-mRNA) and specialized in RNA splicing of pre-mRNA sequences.
-
•
3UTRBERT [24]. This model was trained on k 3’UTRs for 3’UTR-mediated gene regulation tasks. It uses k-mers tokenization instead of SNT.
-
•
CDSBERT [39]. It was trained on massive CDS regions within the genome. CDSBERT is adapted from ProteinBERT for predicting protein structures and functions and leverages the alignment between RNA and protein sequences.
-
•
OmniGenome. OmniGenome is the first genomic FM that highlights the importance of sequence-structure alignment and is available for RNA design tasks.
-
•
OmniGenome555Please find the finetuning example of OmniGenome in the supplemental materials.. OmniGenome is an enhanced variant of OmniGenome that feeds both sequences and structures into OmniGenome to aggregate the feature representations to improving modeling ability.
9 Genetic Algorithm for RNA Sequence Design
9.1 Genetic Algorithm
![Refer to caption](x5.png)
The working mechanism of our designed genetic algorithm based on OmniGenome is implemented as the following five-step process:
-
Step .
Given the target RNA secondary structure, we use OmniGenome to generate a set of candidate sequences .
-
Step .
If the termination criterion is not met, go to Step ; otherwise, output the current best sequence .
-
Step .
Based on , use single-point crossover and mutation to generate a population of offspring sequences .
-
Step .
Combine and to obtain , and use OmniGenome to predict the corresponding secondary structures of each sequence in . Evaluate the fitness values of sequences in .
-
Step .
Sort according to the fitness values and preserve the best sequences to constitute a new . Return to Step .
Note that the fitness value of a sequence , denoted as , is evaluated as the Hamming distance of the RNA secondary structure predicted by OmniGenome against the target structure. The above genetic algorithm is not terminated until the sequence for the target RNA secondary structure is identified or the allocated computational budget is exhausted.
![Refer to caption](x6.png)
9.2 RNA Design Case Study
As demonstrated in the zero-shot experiments in Table 3, OmniGenome models achieve state-of-the-art performance. We show several complex examples of puzzles from the EternaV2 benchmark. According to Fig. 6, puzzles #5 and #11 with approximately + bases are solved, while these puzzles are challenging to existing FMs. Even for puzzles that are not completely solved, e.g., puzzles #3 and #27, OmniGenome generates very similar structures, where the nucleotide base difference ratio between the designed structure and the target structure is only . This finding indicates the proficiency of OmniGenome models in solving challenging single-nucleotide resolution genomic tasks.
10 Benchmark Suites
10.1 RNA Genomic Benchmark
Details of the RGB can be found in Table 7. RGB contains SN-level tasks that are curated or collected from published articles. The purpose of RGB is to benchmark genomic FMs in challenging SN-level modeling tasks such as detection and repair of SN mutations, mRNA sequence degradation rates, and RNA secondary structure prediction. Due to the lack of a plant RNA benchmark dataset, RGB includes the modeling of RNA sequences from a variety of species, e.g., plant and human. The sequence length in RGB ranges from to , which is sufficient for most RNA understanding tasks. In summary, these multi-species and SN-level tasks in RGB serve as the first comprehensive benchmark utilized to assess the RNA sequence modeling capabilities of OmniGenome and its baseline models. The brief introduction of the datasets in RGB is as follows:
-
•
Single-Nucleotide Mutation Detection (SNMD): We developed a plant RNA dataset synthesizing the single-nucleotide mutations. Focused on identifying potential single nucleotide changes, this task is essential for detecting mutations linked to genetic disorders. The SNMD dataset introduces up to random mutations in the original sequences, regardless of variation ratios. Cross-entropy is utilized as the loss function for this binary token classification task.
-
•
Single-Nucleotide Mutation Repair (SNMR): This task challenges the model to suggest corrective actions at the single nucleotide level, aiding in gene therapy approaches. The SNMR dataset mirrors the SNMD dataset, with cross-entropy as the loss function, indicating a token 4-way (i.e., A, U, C, G) classification task.
-
•
mRNA Degrade Rate Prediction (mRNA): Estimating the decay rate of nucleotides in mRNA sequences, this task is vital for deciphering gene expression and regulation. The dataset originates from the Kaggle COVID-19 vaccine design competition666https://www.kaggle.com/competitions/stanford-covid-vaccine, focusing solely on sequence-based degradation rate prediction and excluding RNA structures. It’s a token regression task using MSE as the loss function, with the dataset resplit into training, validation, and testing sets for evaluation.
-
•
RNA Secondary Structure Prediction (bpRNA & Archive2 & RNAStralign & Rfam): Aiming to predict RNA folding into secondary structures, this task is fundamental to RNA functionality and interactions. We evaluated OmniGenome on four datasets, bpRNA [5] (TR0, VL0, TS0 sets), ArchiveII [6], RNAStralign [4] and Rfam [7]. Following existing works, we have excluded sequences over bases and complex structures, simplifying to three symbols: ‘(’, ‘.’, ‘)’Ṙesults may not directly compare with other studies due to these modifications. Cross-entropy serves as the loss function.
Task | Task Type | # of examples | # of classes | Metric | Sequence length | Source |
---|---|---|---|---|---|---|
SNMD | Token classification | AUC | This work | |||
SNMR | Token classification | F1 | This work | |||
mRNA | Token regression | — | RMSE | Kaggle | ||
bpRNA | Token classification | F1 | [5] | |||
AchiveII | Token classification | F1 | [6] | |||
RNAStrAlign | Token classification | F1 | [4] | |||
Rfam∗ | Token classification | F1 | [7] |
Please find the appendix for the input and output examples of each subtask in RGB. The detailed task descriptions for each nucleic acid and species, including the number of examples, classes, evaluation metric, and sequence length, are outlined in Table 7. Each task is carefully curated to reflect the complexity and variety inherent in genomic data, providing a robust framework for assessing the nuanced capabilities of state-of-the-art RNA FMs.
Table 8 show the virtual examples of different datasets in RGB. Please refer to our supplementary materials to find the datasets for more details.
Genome Type | Dataset | Examples | |
RNA | SNMD | Input Sequence | G A G T A T T G A G |
True Label | 0 0 1 0 0 0 0 1 0 0 | ||
Prediction | 0 0 0 0 0 0 0 1 0 0 | ||
SNMR | Input Sequence | T A C G A C T G A T | |
True Label | T A C A A G T A A T | ||
Prediction | T A C A A C T G A T | ||
mRNA | Input Sequence | G G A C | |
True Label | [0.1,0.3,0.2] [0.8,0.4,0.1][0.9,0.4,0.3] [0.5,0.2,0.6] | ||
Prediction | [0.1,0.3,0.2] [0.8,0.4,0.1][0.9,0.4,0.3] [0.5,0.2,0.6] | ||
bpRNA | Input Sequence | G G C G A C U U U U | |
True Label | ( ( ( ) ) ) | ||
Prediction | ( ( ( ( ) ) ) ) | ||
Archive2 | Input Sequence | A G U A G U U U G C U | |
True Label | ( ( ( ) ) ) | ||
Prediction | ( ( ( ) ) ) | ||
RNAStralign | Input Sequence | A G U A G U U U G C U | |
True Label | ( ( ( ) ) ) | ||
Prediction | ( ( ( ) ) ) | ||
Rfam | Input Sequence | A G U A G U U U G C U | |
True Label | ( ( ( ) ) ) | ||
Prediction | ( ( ( ) ) ) |
10.2 Plant Genomic Benchmark
The Plant Genomic Benchmark [14] (PGB) provides a comprehensive suite of datasets designed to evaluate and improve the predictive capabilities of genomic models in plant biology. This benchmark, as shown in Table 9, encompasses a range of critical genomic tasks777https://huggingface.co/InstaDeepAI/agro-nucleotide-transformer-1b, including binary classification, single and multi-variable regression, and multi-label classification, addressing various aspects of plant genomics such as RNA processing, gene expression, and chromatin accessibility. By integrating diverse genomic tasks, the PGB aims to facilitate advanced research and development in plant genomics, offering a robust platform for the assessment and enhancement of model performance across different plant species. To obtain a detailed description of PGB, please refer to Agro-NT [14].
Task | # of datasets | Task Type | Total # of examples | # of classes | Metric | Sequence length |
---|---|---|---|---|---|---|
Polyadenylation | Sequence classification | F1 | ||||
Splice site | Sequence classification | F1 | ||||
LncRNA | Sequence classification | F1 | ||||
Promoter strength | Sequence regression | — | RMSE | |||
Terminator strength | Sequence regression | — | RMSE | |||
Chromatin accessibility | Multi-label classification | F1 | ||||
Gene expression | Multi-variable regression | — | RMSE | |||
Enhancer region | Sequence classification | F1 |
11 OneKP Initiative
There has been a variety of FMs utilized in different species, e.g., humans [3, 13], bacteria [2], and viruses [42], which indicates the effectiveness of pretrained FMs on multi-species genomics. In this work, we aim to propose an FM for multi-species plant RNA sequence modeling. We leverage the OneKP initiative [23] to address the scarcity of plant RNA data, which contains species of plant transcriptomes. The scale of OneKP enables the development of a more robust and transferable RNA FM.
The Plant Transcriptomes Initiative (OneKP) was a comprehensive effort aimed at exploring genetic diversity across the green plant kingdom (Viridiplantae), sequencing the RNA from ( in other versions) samples that represent over species, encompassing all major taxa within Viridiplantae. This includes streptophyte and chlorophyte green algae, bryophytes, ferns, angiosperms, and gymnosperms. The initiative’s final or capstone publication presents three major analyses: inferring species trees, identifying whole genome duplications, and detecting gene family expansions. These findings are particularly valuable for plant and evolutionary scientists interested in specific gene families, whether their focus is across the entire green plant tree of life or within more narrowly defined lineages.
The sampling strategy for the 1KP was global and collaborative, with samples sourced from a wide range of environments including wild field collections, greenhouses, botanical gardens, laboratory specimens, and algal culture collections. The initiative prioritized the collection of live growing cells, such as young leaves, flowers, or shoots, to ensure a high abundance of expressed genes, though many samples also came from roots and other tissues. RNA extraction was performed using well-established protocols or commercial kits, facilitating the comprehensive analysis of transcribed RNA across this diverse set of species. This monumental effort not only sheds light on plant genetic diversity but also provides a rich data resource for ongoing and future research in plant science and evolutionary biology.
12 Limitations
The limitations of our work are rooted in resource constraints. First, according to the data scaling laws [43, 44, 45], we recognize that the scale of the proposed RNA foundation model remains small, and the parameterization may not be adequate to fully exploit the OneKP database’s potential. Due to resource limitations, we could not pretrain larger models. Moving forward, we aim to train larger-scale foundation models, anticipating that our contributions will accelerate the advancement of DNA and RNA foundation models. Secondly, while the modeling length of our foundation model satisfies most RNA and DNA sequence processing needs, there are tasks with extremely long sequences. Our future efforts will focus on enhancing the model’s capability to handle such lengthy sequences.
13 Ethics Statement
In this research, we utilized the open OneKP dataset, which does not contain human-related privacy concerns. We ensure that such data is not exploited without fair compensation and acknowledgment of the source communities. The pretraining sequences are plant-based genomic data that pose potential harm to ecological systems; we do not permit the use of our model beyond expectations, such as developing malicious bio-software or designing harmful RNA structures. The models and findings should support, not undermine, the conservation of plant species and their habitats. We adhere to principles of transparency and open science, using datasets that are publicly available and providing clear documentation of our methodologies and findings.
Overall, in conducting this research, we have committed to ethical scientific practices that respect biodiversity and aim to contribute positively to the field of genomic research. We encourage ongoing dialogue around the ethical use of plant RNA sequences and support initiatives that promote the sharing of benefits arising from such research with all stakeholders.