# MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses

Sebastien A Choteau, Audrey Wagner, Philippe Pierre, Lionel Spinelli, Christine Brun

<figure><img src="image_1.png" /></figure>

To cite this version:

Sebastien A Choteau, Audrey Wagner, Philippe Pierre, Lionel Spinelli, Christine Brun. MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses. Database - The journal of Biological Databases and Curation, 2021, 2021, 10.1093/database/baab032. hal-03268485

HAL Id: hal-03268485

https://amu.hal.science/hal-03268485

Submitted on 23 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Original article

# MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses

Sebastien A. Choteau ID ¹,², Audrey Wagner¹, Philippe Pierre²,³,⁴, Lionel Spinelli¹,² and Christine Brun ID ¹,⁵,*

¹Aix-Marseille University, INSERM, TAGC, Turing Centre for Living Systems, 163 Avenue de Luminy, Marseille 13009, France, ²Aix-Marseille University, INSERM, CNRS, CIML, Turing Centre for Living Systems, 163 Avenue de Luminy, Marseille 13009, France, ³Department of Medical Sciences, Institute for Research in Biomedicine (iBiMED) and Ilidio Pinho Foundation, University of Aveiro, Aveiro 3810-193, Portugal, ⁴Shanghai Institute of Immunology, School of Medicine, Shanghai Jiao Tong University, Shanghai, China and ⁵CNRS, 31 Chemin Joseph Aiguier, Marseille 13009, France

*Corresponding author: Tel: +33 4 91 82 87 12; Fax: +33 4 91 82 87 01; Email: christine-g.brun@inserm.fr

Citation details: Choteau, S.A., Wagner, A., Pierre, P. et al. MetamORF: a repository of unique short open reading frames identified by both experimental and computational approaches for gene and metagene analyses. *Database* (2021) Vol. 2021: article ID baab032; doi:10.1093/database/baab032

Received 23 December 2020; Revised 8 April 2021; Accepted 17 May 2021

## Abstract

The development of high-throughput technologies revealed the existence of non-canonical short open reading frames (sORFs) on most eukaryotic ribonucleic acids. They are ubiquitous genetic elements conserved across species and suspected to be involved in numerous cellular processes. MetamORF (https://metamorf.hb.univ-amu.fr/) aims to provide a repository of unique sORFs identified in the human and mouse genomes with both experimental and computational approaches. By gathering publicly available sORF data, normalizing them and summarizing redundant information, we were able to identify a total of 1 162 675 unique sORFs. Despite the usual characterization of ORFs as short, upstream or downstream, there is currently no clear consensus regarding the definition of these categories. Thus, the data have been reprocessed using a normalized nomenclature. MetamORF enables new analyses at locus, gene, transcript and ORF levels, which should offer the possibility to address new questions regarding sORF functions in the future. The repository is available through an user-friendly web interface, allowing easy browsing, visualization, filtering over multiple criteria and export possibilities. sORFs can be searched starting from a gene, a transcript and an ORF ID, looking in a genome area or browsing the whole repository for a species. The database content has also been made available through track hubs at UCSC Genome Browser. Finally, we demonstrated
an enrichment of genes harboring upstream ORFs among genes expressed in response to reticular stress.

Database URL: https://metamorf.hb.univ-amu.fr/

## Introduction

Short open reading frames (sORFs) are usually defined as sequences delimited by a start codon and a stop codon and potentially translatable into proteins of <100 amino acids (1–8). They are present in all classes of transcripts [including presumptive long non-coding ribonucleic acids (lncRNAs)] and have been identified in most eukaryotic RNAs (2, 5, 8–15). In addition, their sequence often begins with a non-canonical start codon (8). Consequently, they have long been overlooked, and interest in their possible regulatory functions has only raised recently with the advent of the ribosome profiling method that strongly suggests their translation (1, 3, 5, 6, 16–22).

Several sORF categories have been defined according to their location on RNAs (Figure 1). For instance, upstream ORFs (uORFs) are located in the 5' untranslated regions (5' UTRs) of messenger RNAs (mRNAs) and have been defined as sORFs whose start codon precedes the main coding sequence (CDS; 6, 8, 17, 18, 23). They are conserved across species (5, 6, 11, 21, 24), but less conserved than canonical protein-coding ORFs (25). To date, uORFs have been essentially reported as gene-expression cis-regulatory elements that regulate the efficiency of translation initiation of the main CDS, notably alleviating the repression of translation during cellular stress (13, 17, 18, 20, 23, 26). Moreover, the discovery of uORF-encoded peptides, and more generally sORF-encoded peptides, led to the assumption that they may also play functional roles in trans (2–4, 7, 9, 10, 18, 24, 27–30), for instance as ligands of major histocompatibility complex molecules (12, 22, 23). Very interestingly, uORF-encoded peptides have also been shown to form protein complexes with the protein encoded by the main CDS of the same mRNA (31), and it has been suggested that polycistronic sequences may exist in eukaryotes (24, 31). Furthermore, given the increasing evidence on the regulatory functions of peptides encoded by sORFs located within mRNAs, introns of pre-mRNAs, lncRNAs and primary transcripts of microRNAs or ribosomal RNAs (2, 8–15, 26), there is an urgent need to study sORFs (i) individually and (ii) at the whole proteome scale. Indeed, the latter should reveal important features of sORFs, thus enabling the characterization and the identification of their functions. However, the fact that (i) the publicly available data are scattered across different databases and (ii) datasets are aligned on different genome builds, differently annotated and formatted, calls for an uniformed resource

where each sORF is individually described. With this in mind, we have built a resource database of publicly available sORFs identified in the human and mouse genomes, by gathering information from computational predictions and Ribo-seq and proteomic experiments. The curation of data, their homogenization in order to merge the redundant information into unique entries, the completion and computation of missing information (e.g. sequences and Kozak contexts) and the re-annotation of sORF classes represent the added value of this database. Notably, this enables the analyses at locus, gene, transcript and ORF levels. In this work, we propose (i) a pipeline to regularly update the content of the database in a reproducible manner, (ii) a database content that can be fully downloaded for custom computational analyses and (iii) an user-friendly web interface to ease data access to biologists.

## Material and methods

### MetamORF pipeline and database development

#### Inclusion criteria for publicly available sORF-related data

A total of 18 data sources, either *Homo sapiens* and *Mus musculus* original datasets or re-processed publicly available sORFs repositories, have been considered for inclusion in our database (Supplementary Table S1) (5, 7, 11, 12, 14, 15, 17–22, 32–37). These data sources provide results from computational predictions, Ribo-seq experiment analyses and mass spectrometry (proteomics/proteogenomics) analyses. The data sources not providing the absolute genomic coordinates of the ORF start and stop codons (5, 17, 20, 32–34) or fully included in another data source considered here (21) have been discarded. Databases that did not allow export of their content in a single file or automating the download of all the files from their website have also been discarded (19, 35). Despite their short size, it has been noticed that sORFs can be spliced. Theoretical lengths of the ORFs have been computed as the distance between the start and stop codons, eventually removing the intron length(s) when information about ORF splicing was provided. Due to splicing, the theoretical length and the one reported by the data source may be different. Data sources harboring this difference for >95% of their entries were discarded as this indicates the splicing information was missing (10). Finally, data sources for which we were not able to perform this assessment as they were not providing
<figure><img src="image_4.png" /><figcaption>**Figure 1.** MetamORF pipeline. This figure represents the workflow used to build MetamORF. First, the data from the sources selected have been inserted into the database, and the absolute genomic coordinates have been homogenized from their original annotation version to the most recent version (GRCh38 or GRCm38). Then, the redundant information, i.e. the entries describing the same ORFs (same start, stop and splicing), has been merged, allowing to get one single and unique entries for each ORF detected on the human and mouse genomes. The missing information (sequences and transcript biotypes) has been downloaded from Ensembl, and the ORF relative coordinates have been computed. Finally, the cell types and ORF classes have been normalized, and the Kozak contexts have been computed using the sequences flanking the start codons.</figcaption></figure>

information regarding (i) the splicing of the ORF and (ii) ORF length (15, 36) have not been included as well. Hence, the database has been made by collecting data from six distinct sources (Table 1), including either original datasets (Table 1 and Supplementary Table S2) (11, 12, 14, 18, 22) or reprocessed data (37), and discarding 12 of them (Supplementary Table S1). Notably, we have included data from sORFs.org (37), considered as

the main and most comprehensive repository of sORFs identified by genome-wide translation profiling (Ribo-seq), that currently integrates re-processed data from 73 original publications.

For each of these sources, a set of features essential to properly characterize the sORFs, related to their location, length, sequences, environmental signatures and cell types (i.e. cell lines, tissues or organs) in which they are expressed,
<table><caption>**Table 1.** Information about the data sources used to build MetamORF</caption><thead><tr><td>Publication</td><td>DOI</td></tr></thead><tbody><tr><td>Mackowiak et al., 2015, Genome Biol. (11)</td><td>10.1186/s13059-015-0742-x</td></tr><tr><td>Erhard et al., 2018, Nat. Methods (22)</td><td>10.1038/nmeth.4631</td></tr><tr><td>Johnstone et al., 2016, EMBO J. (18)</td><td>10.15252/embj.201592759</td></tr><tr><td>Laumont et al., 2016, Nat. Commun. (12)</td><td>10.1038/ncomms10238</td></tr><tr><td>Samandi et al., 2017, eLife (14)</td><td>10.7554/eLife.27860</td></tr><tr><td>Olexiouk et al., 2018, Nucleic Acids Res. (37)</td><td>10.1093/nar/gkx1130</td></tr></tbody></table>

See Supplementary Table S1 for more information about these data sources.

have been collected (see Table 2 for a full list of features considered for inclusion). When it was not provided by the source, the symbol of the gene related to the sORF was recovered using the transcript identifier (ID, if provided) or searching for the gene(s) or ncRNA(s) overlapping with the sORF coordinates in the original annotation version by querying Ensembl databases (38) in their appropriate versions (v74, 75, 76, 80, 90) with pyensembl (v1.8.5, https://github.com/openvax/pyensembl). In addition to these features, information regarding the transcript(s) harboring the ORFs have been collected from the data sources when available. This is of particular interest as some ORF features, such as the ORF class, may depend on the transcript they are located in (e.g. an ORF may be located in the 5' UTR of a transcript and be overlapping with the CDS of another transcript). Finally, 3379219 and 2066627 entries from these six data sources have been collected and inserted in MetamORF for *H. sapiens* and *M. musculus*, respectively (Table 3).

## Homogenization of genomic coordinates

As the data sources were providing genomic coordinates from different genome annotation versions (e.g. GRCh38 and GRCh37), all the genomic coordinates registered in our database have been lifted over the latest annotation version (GRCh38 for *H. sapiens* and GRCm38 for *M. musculus*) using pyliftover (v0.4, https://pypi.org/project/pyliftover). The liftover has been considered as failed for an entry if (i) at least one of the coordinates (i.e. start, stop or one of the start or end exon coordinates) was located on a strand different from all the others or (ii) the chromosome of the position changed during the liftover or (iii) the distance (in nucleotides) between the sORF start and stop codons has changed after the liftover. All the entries for which the liftover failed were removed from the database. Based on the previous assumptions, the liftover failed for 709 ORFs

(377 failed due to the last criteria) in *H. sapiens* and for none of the *M. musculus* entries (Table 3). The choice of such stringent criteria has been strengthened by the fact that these entries (i) only represent <0.05% of the entries for *H. sapiens* and (ii) are more susceptible to be unreliable entries.

## Merge of redundant information

As our database aims to provide a repository of unique identified sORFs of the human and mouse genomes, all the redundant entries describing the same sORFs have been merged. In a first step, we identified all the sORF entries for which all the identification features were provided (chromosome, strand, start position, stop position, splicing status and splicing coordinates). sORFs sharing the same feature values were merged. In a second step, we identified all the remaining entries with only partial identification features provided: the chromosome as well as either (i) both the strand and the start positions or (ii) both the strand and the stop positions or (iii) both the start and the stop positions. Those entries were merged to the best matching fully described entries identified in the first step. If no matching fully described entry was found, then the entries were removed. In order to keep track of the number of times a same sORF has been described in the original data sources, the initial number of entries merged together was registered for each sORF.

During this merging, information regarding the transcripts that harbor the sORFs has been registered too. Hence, when several sORFs were merged into one single entry in MetamORF, the resulting new sORF entry was registered as harbored by all the distinct transcripts related to the original entries. After this removal of redundant information, we were finally able to identify 664771 and 497904 unique sORFs for *H. sapiens* and *M. musculus*, respectively (Table 3).

It should be noticed that all unique sORF entries generated at this stage have been kept, including the ones describing ORFs longer than 100 amino acids. Entries describing such ORFs may be either coming from data sources that (i) did not remove the ORFs longer than 100 amino acids or (ii) used a higher threshold or (iii) described the ORF as unspliced while it is actually susceptible to be spliced (and thus has a shorter sequence on the transcript than the one expected).

## Completion of missing information and computation of relative coordinates

In the original data sources, the only information provided (when provided) on the transcripts was the transcript ID. Detailed information was retrieved from Ensembl
<table><caption>**Table 2.** Features allowing to characterize the sORFs</caption><thead><tr><th>Family</th><th>Feature</th><th>Details</th></tr></thead><tbody><tr><td rowspan="10">Location</td><td>Chromosome</td><td>The chromosome or scaffold on which the ORF is located</td></tr><tr><td>Strand</td><td>The strand of the sORF</td></tr><tr><td>ORF start</td><td>The absolute genomic coordinates of the start codon (position of the first nucleotide)</td></tr><tr><td>ORF stop</td><td>The absolute genomic coordinates of the stop codon (position of the third nucleotide)</td></tr><tr><td>Splicing status</td><td>Is the sORF spliced?</td></tr><tr><td>Splicing coordinates</td><td>The coordinates of the start and end of each exon constituting the sORF</td></tr><tr><td>Transcript</td><td>The name or ID of the transcript(s) related to the sORF (eventually with transcript strand, start and end positions and transcript biotype)</td></tr><tr><td>Gene</td><td>The name, symbol, alias or ID of the gene(s) related to the sORF (when not intergenic)</td></tr><tr><td>Lengths</td><td>Length</td><td>The length of the sORF (in nucleotides)</td></tr><tr><td>Category</td><td>Putative sPEP length</td><td>The length of the (putative) sORF-encoded peptide in amino acids</td></tr><tr><td>Sequence signature</td><td>Start codon sequence</td><td>The nucleic sequence of the sORF start codon</td></tr><tr><td>Environmental signature</td><td>Nucleic sequence</td><td>The nucleic sequence of the sORF</td></tr><tr><td>Conservation</td><td>Amino acid sequence</td><td>The amino acid sequence of the (putative) sORF-encoded peptide</td></tr><tr><td>Coding potential assessment</td><td>Kozak context</td><td>Does a Kozak context has been identified for the sORF start codon?</td></tr><tr><td>Biological context</td><td>PhyloCSF score</td><td>The PhyloCSF score computed for the sORF</td></tr><tr><td>Sequence signature</td><td>PhastCons score</td><td>The PhastCons score computed for the sORF</td></tr><tr><td>Environmental signature</td><td>FLOSS class and score</td><td>The FLOSS class and score computed for the sORF</td></tr><tr><td>Conservation</td><td>ORF score</td><td>The ORF score computed for the sORF</td></tr><tr><td>Conservation</td><td>Cell context</td><td>The cellular context in which the sORF has been identified or detected</td></tr></tbody></table>

<table><caption>**Table 3.** MetamORF most important statistics</caption><thead><tr><th>Feature</th><th></th><th>H. sapiens</th><th>M. musculus</th></tr></thead><tbody><tr><td rowspan="10">Original data sources</td><td>ORFs</td><td>1 344 978</td><td>1 249 176</td></tr><tr><td>Transcripts</td><td>101 597</td><td>85 653</td></tr><tr><td>Predicted ORFs for which the transcript is unknown</td><td>181 122</td><td>213 301</td></tr><tr><td>ORFs detected by Ribo-seq for which the transcript is unknown</td><td>79 422</td><td>8546</td></tr><tr><td>ORFs detected by MS for which the transcript is unknown</td><td>54</td><td>0</td></tr><tr><td>ORF to transcript associations</td><td>3 379 219</td><td>2 066 627</td></tr><tr><td>ORFs predicted</td><td>202 309</td><td>222 705</td></tr><tr><td>ORFs identified by ribosome profiling</td><td>1 142 669</td><td>1 026 471</td></tr><tr><td>ORFs identified by MS</td><td>166</td><td>0</td></tr><tr><td>ORFs for which the homogeneization of genomic coordinates failed</td><td>709</td><td>0</td></tr><tr><td rowspan="10">MetamORF database</td><td>ORFs</td><td>664 771</td><td>497 904</td></tr><tr><td>Transcripts</td><td>90 406</td><td>63 147</td></tr><tr><td>Predicted ORFs for which the transcript is unknown</td><td>13 440</td><td>14 327</td></tr><tr><td>ORFs detected by Ribo-seq for which the transcript is unknown</td><td>71 158</td><td>2</td></tr><tr><td>ORFs detected by MS for which the transcript is unknown</td><td>48</td><td>0</td></tr><tr><td>ORF for which the transcripts are unknown</td><td>83 403</td><td>14 329</td></tr><tr><td>ORF to transcript associations</td><td>729 793</td><td>696 785</td></tr><tr><td>ORFs predicted</td><td>17 027</td><td>14 500</td></tr><tr><td>ORFs identified by ribosome profiling</td><td>664 771</td><td>497 904</td></tr><tr><td>ORFs identified by MS</td><td>147</td><td>0</td></tr><tr><td>Genes harboring at least 1 sORF</td><td>23 767</td><td>15 869</td></tr><tr><td>ORFs having at least one class annotation (short, upstream)</td><td>630 953</td><td>497 904</td></tr></tbody></table>

MS: mass spectrometry.
databases (v90) through their REST API and inserted in our database: (i) the transcript biotype, (ii) the transcript start and end genomic coordinates, (iii) the codon of the canonical CDS (for protein-coding transcripts only) start and stop genomic coordinates and (iv) the full nucleic sequence. In addition, the sequence flanking the start codon (20) has been recovered. As the sORF nucleic and amino acid sequences were not systematically provided by the data sources, these were downloaded from the Ensembl databases using their genomic coordinates.

Moreover, when the transcript ID was available, sORF start and stop relative coordinates have been computed on each of their transcript using AnnotationHub (v2.18.0; 39) and ensembldb (v2.10.2, https://bioconductor.org/package s/release/bioc/html/ensembldb.html) R packages (R v3.6.0).

# Standardization of the cell types and ORF classes

## Cell types

Original data sources do not use a common thesaurus or ontology to name the cell types (e.g. 'HFF' and 'Human Foreskin Fibroblast') or use non-biological meaning names (e.g. sORFs.org (37) provides the name of the original publication as a cell type). In order to provide an uniform informative naming, we manually recovered the name of the cell line, tissue or organ used in these datasets and defined an unique name to be used in our database for each cell line, tissue or organ, trying to use the most commonly used nomenclature for cell lines (Supplementary Table S3). In addition, in order to ensure interoperability with other biological resources, we recovered the matching ontology terms from the following ontologies when feasible: the Cell Ontology (40), the Cell Line Ontology (41), the BRENDA Tissue Ontology (42), the Human Cell Atlas Ontology (43), the Foundational Model of the Anatomy Ontology (44), the Ontology for Biomedical Investigations (45), the NCI Thesaurus OBO Edition (46), the Experimental Factor Ontology (47), the BioAssay Ontology (48) and the Ontology for MIRNA Target (49), using the Ontology Lookup Service (EBI) (50) (Supplementary Table S4).

## ORF classification

Despite the use of a common nomenclature by the wide majority of the scientific community to annotate the open reading frames, based on their size and relative position on their transcript (e.g. short, upstream, downstream and overlapping), no clear consensus about the definitions of these categories nor their names has been defined so far (25). In order to homogenize this information in MetamORF, we created a new annotation of the ORFs using the ORF length, transcript biotype, relative positions and reading frame information when available (see Supplementary Methods). In this annotation, a threshold of 100 amino

acids has been used to define the 'short ORFs', as this value is the most commonly used for historical reasons (2, 4, 6, 8, 24).

## Computation of the Kozak contexts

The Kozak motif and context have been regarded as the optimal sequence context to initiate translation in all eukaryotes. We have thus assessed the Kozak context for each sORF, using the criteria defined by Hernández et al. (51). Briefly, for each ORF to transcript association, the Kozak context was computed looking for regular expression characterizing an optimal, strong, moderate or weak Kozak context (Supplementary Tables S5 and S6). Kozak-alike contexts were also computed for non-ATG initiated sORFs looking for the same patterns with flexibility regarding nucleotides at +1 to +4 positions.

## MetamORF software and languages

The pipeline used to build MetamORF has been developed using Python (v2.7) with SQLAlchemy ORM (sqlKElly.org, v1.3.5). The database has been handled using MySQL (mysql.com, v8.0.16). Docker (docker.com, v18.09.3) and Singularity (singularity.lbl.gov, v2.5.1) environments have been used in order to ensure reproducibility and to facilitate deployment on high-performance clusters.

The MetamORF web interface has been developed using the Laravel (laravel.com, v7.14.1) framework with PHP (v7.3.0), JavaScript 9, HTML 5 and CSS 3. The NGINX (v1.17.10) web server and PHP server (v7.3.0) were deployed with Docker (docker.com, v18.09.3) and Docker-compose (v1.24.0) to ensure stability.

## Enrichment analysis

### Gene lists

The list of genes harboring at least one uORF has been collected from MetamORF as a list of Ensembl identifiers using a SQL query.

The list of ATF4 and CHOP targets identified by ChIP-seq comes from Han et al. (52) (available as supplementary data on the editor's website). Genes congruently and translationally upregulated under endoplasmic reticular (ER) stress have been provided by Guan et al. (53) (upon request). As these lists of genes were provided as gene symbols, matching Ensembl IDs have been recovered using the g:Convert tool available on the gProfiler web interface (54).

The universe contains all protein-coding genes annotated at least once in Gene Ontology (55, 56) (downloaded from the g:Profiler web interface on 3 November 2020).

### Statistics

After discarding genes absent in the universe from the lists, the enrichment analysis was performed using
an hypergeometric test with R 3.6.0 (https://www project.org/). A Benjamini-Hochberg correction has been applied to allow for multiple comparisons, and a False Discovery Rate (FDR) threshold of 0.05 has been considered as significant.

## Database content, accessibility and web interface

### A new repository of short ORF-related data

MetamORF describes 664 771 and 497 904 unique ORFs in the human and mouse genomes, respectively, providing at least the information necessary to locate the ORF on the genome, its sequence and the gene it is located on (excepted for intergenic ORFs). Extensive information related to the transcripts is provided for 614 997 (~93%) and 497 904 (100%) sORFs for the human and mouse genomes, respectively. These features allowed us to classify 630 953 (~95%) human ORFs and 497 904 (100%) mouse ORFs in at least one class (Table 3, Figure 2, Supplementary Figure S1). Interestingly, it should be noticed that a large proportion (36% and 52% for *H. sapiens* and *M. musculus*, respectively) of ORFs are using an alternative frame to the main CDS. In addition, nearly 23% of the ORFs are located on non-coding RNAs for both species.

### User-friendly web interface and genome tracks

To provide users with a clear, fast and easy-to-use database, MetamORF can be queried through an user-friendly web interface at https://metamorf.hb.univ-amu.fr. A tutorial as well as a documentation page are available online. Briefly, the users may search for sORFs contained in the database starting with a gene symbol (symbol, alias, ID), a transcript ID (ID, name) and an ORF ID or screening a particular genomic area. The data are made accessible through four types of pages: (i) a 'gene' page (Figure 3) to allow visualizing information related to all transcripts and sORFs on a gene, (ii) a 'transcript' page to allow browsing information related to a transcript gene and all its sORFs, (iii) an 'ORF' page to allow fetching information related to all transcripts and gene that harbor the chosen ORF and finally (iv) a 'locus' page to allow getting information related to all sORFs located in a particular locus. In addition, the user may also browse across all sORFs related to a species or detected in a particular cell type. It is possible to navigate from one to another page easily to get extensive information about a sORF, a gene or a transcript (Supplementary Figure S2).

In each page, the results can be filtered on (i) the identification method (computational prediction, ribosome

profiling or mass spectrometry), (ii) the start codon, (iii) the Kozak context (as previously defined), (iv) the genomic length (defined as the sum of lengths of each exon constituting the ORF), (v) the transcript biotype (according to the Ensembl definitions), (vi) the ORF annotation (as previously defined) and (vii) the cell type (Supplementary Tables S3 and S4).

All results can be exported in an easily parsable format (comma-separated values file, CSV), as well as in FASTA or BED format.

On ORF, transcript and locus pages, a link allowing the user to easily visualize all the ORFs localized in a particular area on the UCSC Genome Browser (57) is proposed. We also implemented genome track hubs, to allow using UCSC Genome Browser advanced options, such as filtering on ORF categories, transcript biotypes, cell types and transcript IDs.

In addition to this user-friendly interface, it is possible to download from the website the content of the full MetamORF database in BED and FASTA formats.

### Using MetamORF to analyze the regulation of integrated stress response

Several studies have reported the role of uORFs in the regulation of the translation during the integrated stress response (ISR) (13, 23, 26, 28). Notably, the mechanism by which the repression of the translation is alleviated under an ER stress has been elucidated for the mammalian transcription factor ATF4, the targets of which are responsible for cell adaptation to stress. Briefly, ATF4 CDS is preceded by two functional uORFs (58), both highly expressed under normal growth and stress conditions. Under the ISR, the small ribosomal subunit is expected to remain bound to the mRNA, scan through the uORF2 and acquire the eIF2•GTP•Met-tRNAᵢMet and the large ribosomal subunit in time for initiation at the start codon of the CDS, a phenomenon known as 'leaky scanning'. In addition, it has been also suggested that the translation of the CDS under stress may result from the 're-initiation', a model in which the large ribosomal subunit and the initiation complex are recruited by the small subunit right after the termination of the translation of the uORF2, allowing thus the initiation at the CDS start codon. Both events are nevertheless technically difficult to distinguish and the exact process remains debated. Hence, assuming the presence of one uORF is sufficient to regulate the translation of the CDS (20), are targets of ATF4 and CHOP (another transcription factor activated upon stress) more likely to harbor uORFs than other genes? Are
A H. sapiens

<figure><img src="image_5.png" /></figure>

B M. musculus

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Count of ORFs in each class. The bar plots represent the count of ORFs annotated for each class for (A) *H. sapiens* and (B) *M. musculus*. The percentages displayed over the bars indicate the proportion of ORFs annotated in the class over the total number of ORFs registered in the database for the species. NMD: non-sense-mediated decay; NSD: non-stop decay.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** MetamORF gene-centric view. The page displays the transcripts and the ORFs related to *SGK3* gene. A filter has been applied to select exclusively the ORFs detected in HFF, Jurkat, RPE-1, HEK293 or HeLa cells. Other filters may be used and the results can be exported as CSV, FASTA or BED files.</figcaption></figure>

<table><caption>**Table 4.** Enrichment analysis</caption><thead><tr><th>Gene list<sup>a</sup></th><th>List size</th><th>Genes harboring uORFs</th><th>Intersection size</th><th>Universe<sup>a</sup> size</th><th>FDR</th><th>Odds ratio</th></tr></thead><tbody><tr><td>ATF4 targets</td><td>392</td><td>8863</td><td>256</td><td>19 985</td><td>5.52.10<sup>-17</sup></td><td>2.40</td></tr><tr><td>CHOP targets</td><td>256</td><td>8863</td><td>166</td><td>19 985</td><td>3.34.10<sup>-11</sup></td><td>2.34</td></tr><tr><td>Genes congruently upregulated</td><td>484</td><td>8863</td><td>268</td><td>19 985</td><td>5.41.10<sup>-7</sup></td><td>1.57</td></tr><tr><td>Genes transitionally upregulated</td><td>1068</td><td>8863</td><td>736</td><td>19 985</td><td>1.21.10<sup>-61</sup></td><td>2.94</td></tr></tbody></table>

genes translationally or congruently upregulated during an ER stress, enriched in genes harboring uORFs? To answer these questions, we performed enrichment analyses, getting the list of genes harboring uORFs by querying MetamORF, and using the published lists of target genes of ATF4 and CHOP identified by ChIP-seq (Supplementary Table S7). We demonstrated that ATF4 and CHOP targets as well as genes upregulated under an ER stress are more likely to harbor uORFs than expected by chance (OR<sub>ATF4</sub> = 2.40, pval = 2.76.10<sup>-17</sup> and OR<sub>CHOP</sub> = 2.34, pval = 2.50.10<sup>-11</sup>, respectively; Table 4). This suggests that the translation of these genes is likely to be under the control of uORFs, as it has been experimentally shown for PPP1R15A and PPP1R15B (23), two well-known targets of ATF4.

## Discussion and conclusion

MetamORF contains data about 1 162 675 unique sORFs for the human and mouse genomes identified by both experimental and computational approaches. While the Ribo-seq is considered by most as the ‘gold standard’ method to identify sORFs experimentally, the added value of predictive computational approaches, proteomics and peptidomics to characterize such biological sequences remains certain. Because these technologies are offering complementary information at genomic, transcriptomic and proteomic scales, we decided to include data from both experimental and computational experiments in our database. Nevertheless, data coming from distinct data sources may be difficult to compare, in particular because they are not necessarily using the same genome annotation and definitions of ORF

<footnote><sup>a</sup>See Supplementary Table S7 for more information about the gene lists.</footnote>
classes and Kozak contexts, for instance. By homogenizing
this information, MetamORF offers the possibility to com-
pare datasets coming from different sources. We noticed
that information regarding the Kozak context is missing
most of the time, and start flanking sequences are usually
not provided. Hence, MetamORF provides a new interest-
ing set of information. It is noteworthy that we discarded
12 of 18 datasets because they lack crucial information
regarding their integration into MetamORF. Although this
is a rather drastic method, this is performed for the sake of
data quality. In these conditions, the confidence in the data
and the reliability in the existence of the sORF of interest
can be assessed by the number of original experiments that
identify the sORF (column ‘EXP. COUNT’ in the tables
of the web interface). It is noteworthy that >97% of the
unique ORF entries registered in MetamORF have been
identified by at least one experimental method.

It should be noticed that a large amount (~80%) of
the sORFs contained in our database have been described
in the sORFs.org repository (37). Despite being the most
prominent sORF database and offering the community
data processed in a normalized way using their own work-
flow, sORFs.org does not provide metagene analyses (1). In
addition, such analysis is made difficult by the absence of
gene names and transcriptomic coordinates as well as the
high redundancy of information contained in the sORF.org
database (37), issues that we addressed with MetamORF. It
is noteworthy that another sORF resource, namely Open-
Prot (59), does not contain ORFs shorter than 30 amino
acids, whereas in MetamORF, sORFs of such size repre-
sent ~50% of the dataset. Of note, 54% of them have been
detected in at least two data sources, therefore reinforcing
their probability of existence. Hence, in comparison with
existing resources (Supplementary Table S8), MetamORF
is complementary and allows analyses at ORF, transcript,
gene and locus levels. In addition, it opens the possibility
of studying sORFs as a group, at a global scale.

The resource is accessible at https://metamorf.hb.univ-
amu.fr and provides an intuitive querying interface to
enable wet-laboratory researchers to easily question this
large set of information. The web interface comes with
advanced filters, notably on computed ORF classes, ORF
start codons, identification methods, Kozak contexts and
cell contexts. Such filters should help end-user biologists
without computational skills to identify and collect infor-
mation about the sORFs important for their topic of inter-
est. Moreover, the implementation of MetamORF content
in track hubs allows both quick and advanced visualiza-
tion of data through the UCSC Genome Browser. Finally,
the database content may be exported in various con-
venient formats widely used by the scientific community
(e.g. FASTA and BED).

We believe that MetamORF is of interest not only to
bioinformaticians working on short ORFs but also to a
wider community, including any biologist who may ben-
efit from knowledge regarding the sORFs located on their
gene, transcript or region of interest. As ribosome profiling
becomes more appreciated and proteomics starts allowing
accurate identification of short peptides, new data describ-
ing sORFs in various conditions will be published in the
next years, and our database is expected to grow accord-
ingly. In particular, the next release of MetamORF is
expected to include data describing the sORFs of other
organisms such as *Drosophila melanogaster*. As a conclu-
sion, we believe that MetamORF should help to address
new questions in the future, in particular regarding the
regulatory functions of the sORFs as well as the functions
of the short peptides they may encode.

**Supplementary data**

Supplementary data are available at *Database* online.

Acknowledgement

We thank Andreas Zanzoni for helpful discussions.

Funding

The project leading to this publication has received funding from the « Investissements d'Avenir » French Government program managed by the French National Research Agency (ANR-16-CONV-0001) and from Excellence Initiative of Aix-Marseille University - A*MIDEX. It has also received funding from the Excellence Initiative of Aix-Marseille University - A*Midex a French "Investissements d'Avenir programme"- Institute MarMaRa AMX-19-IET007.

*Conflict of interest.* The authors have no conflict of interest to declare.

Data availability

Data sources are available on the editor's website or
using the links provided in their original publications.
The source code used to create the database and the
full technical documentation (source code documentation,
manual, database structure and dockerfiles) are available
on GitHub (https://github.com/TAGC-NetworkBiology/
MetamORF). Full content of the database can be down-
loaded in BED and FASTA formats from the Meta-
mORF website, and up-to-date version of track hubs may
be downloaded and/or used with your favorite genome
browser from the link https://metamorf.hb.univ-amu.fr/
hubDirectory/hub.txt. The dump of the database is avail-
able on request.

**References**

1. Martinez,T.F., Chu,Q., Donaldson,C. et al. (2019) Accu-rate annotation of human protein-coding small open reading frames. *Nat. Chem. Biol.*, **16**, 458–468
2. Andrews,S.J. and Rothnagel,J.A. (2014) Emerging evidence for functional peptides encoded by short open reading frames. *Nat. Rev. Genet.*, **15**, 193-204.

3. Aspden,J.L., Eyre-Walker,Y.C., Phillips,R.J. et al. (2014) Extensive translation of small open reading frames revealed by Poly-Ribo-Seq. *eLife*, **3**, e03528.

4. Saghatelian,A. and Couso,J.P. (2015) Discovery and characterization of smORF-encoded bioactive polypeptides. *Nat. Chem. Biol.*, **11**, 909-916.

5. Chew,G.-L., Pauli,A. and Schier,A.F. (2016) Conservation of uORF repressiveness and sequence features in mouse, human and zebrafish. *Nat. Commun.*, **7**, 11663.

6. Ingolia,N.T., Brar,G.A., Stern-Ginossar,N. et al. (2014) Ribosome profiling reveals pervasive translation outside of annotated protein-coding genes. *Cell Rep.*, **8**, 1365-1379.

7. Hao,Y., Zhang,L., Niu,Y. et al. (2018) SmProt: a database of small proteins encoded by annotated coding and non-coding RNA loci. *Brief. Bioinform.*, **19**, 636-643.

8. Pueyo,J.I., Magny,E.G. and Couso,J.P. (2016) New peptides under the s(ORF)ace of the genome. *Trends Biochem. Sci.*, **41**, 665-678.

9. Raj,A., Wang,S.H., Shim,H. et al. (2016) Thousands of novel translated open reading frames in humans inferred by ribosome footprint profiling. *eLife*, **5**, e13328.

10. Zanet,J., Chanut-Delalande,H., Plaza,S. et al. (2016) Small peptides as newcomers in the control of Drosophila development. *Curr. Top. Dev. Biol.*, **117**, 199-219.

11. Mackowiak,S.D., Zauber,H., Bielow,C. et al. (2015) Extensive identification and analysis of conserved small ORFs in animals. *Genome Biol.*, **16**, 179.

12. Laumont,C.M., Daouda,T., Laverdure,J.-P. et al. (2016) Global proteogenomic analysis of human MHC class I-associated peptides derived from non-canonical reading frames. *Nat. Commun.*, **7**, 10238.

13. Couso,J.-P. and Patraquim,P. (2017) Classification and function of small open reading frames. *Nat. Rev. Mol.*, **18**, 575-589.

14. Samandi,S., Roy,A.V., Delcourt,V. et al. (2017) Deep transcriptome annotation enables the discovery and functional characterization of cryptic small proteins. *eLife*, **6**, e27860.

15. McGillivray,P., Ault,R., Pawashe,M. et al. (2018) A comprehensive catalog of predicted functional upstream open reading frames in humans. *Nucleic Acids Res.*, **46**, 3326-3338.

16. Olexiouk,V., Crappé,J., Verbruggen,S. et al. (2016) sORFs.org: a repository of small ORFs identified by ribosome profiling. *Nucleic Acids Res.*, **44**, D324-D329.

17. Wethmar,K., Barbosa-Silva,A., Andrade-Navarro,M.A. et al. (2014) uORFdb—a comprehensive literature database on eukaryotic uORF biology. *Nucleic Acids Res.*, **42**, D60-D67.

18. Johnstone,T.G., Bazzini,A.A. and Giraldez,A.J. (2016) Upstream ORFs are prevalent translational repressors in vertebrates. *EMBO J.*, **35**, 706-723.

19. Wang,H., Yang,L., Wang,Y. et al. (2019) RPFdb v2.0: an updated database for genome-wide information of translated mRNA generated from ribosome profiling. *Nucleic Acids Res.*, **47**, D230-D234.

20. Andreev,D.E., Arnold,M., Kiniry,S.J. et al. (2018) TASEP modelling provides a parsimonious explanation for the ability

of a single uORF to derepress translation during the integrated stress response. *eLife*, **7**, e32563.

21. Fields,A.P., Rodriguez,E.H., Jovanovic,M. et al. (2015) A regression-based analysis of ribosome-profiling data reveals a conserved complexity to mammalian translation. *Mol. Cell.*, **60**, 816-827.

22. Erhard,F., Halenius,A., Zimmermann,C. et al. (2018) Improved Ribo-seq enables identification of cryptic translation events. *Nat. Methods*, **15**, 363-366.

23. Starck,S.R., Tsai,J.C., Chen,K. et al. (2016) Translation from the 5' untranslated region shapes the integrated stress response. *Science*, **351**, aad3867.

24. Crappé,J., Van Criekinge,W. and Menschaert,G. (2014) Little things make big things happen: a summary of micropeptide encoding genes. *EuPA Open Proteom.*, **3**, 128-137.

25. Orr,M.W., Mao,Y., Storz,G. et al. (2019) Alternative ORFs and small ORFs: shedding light on the dark proteome. *Nucleic Acids Res.*, **48**, 1029-1042.

26. Plaza,S., Menschaert,G. and Payre,F. (2017) In search of lost small peptides. *Annu. Rev. Cell Dev. Biol.*, **33**, 391-416.

27. Hazarika,R.R., Sostaric,N., Sun,Y. et al. (2018) Large-scale docking predicts that sORF-encoded peptides may function through protein-peptide interactions in Arabidopsis thaliana. *PLoS One*, **13**, e0205179.

28. Andreev,D.E., O'Connor,P.,B.F., Fahey,C. et al. (2015) Translation of 5' leaders is pervasive in genes resistant to eIF2 repression. *eLife*, **4**, e03971.

29. Zanet,J., Benrabah,E., Li,T. et al. (2015) Pri sORF peptides induce selective proteasome-mediated protein processing. *Science*, **349**, 1356-1358.

30. Cabrera-Quio,L.E., Herberg,S. and Pauli,A. (2016) Decoding sORF translation – from small proteins to gene regulation. *RNA Biol.*, **13**, 1051-1059.

31. Chen,J., Brunner,A.-D., Cogan,J.Z. et al. (2020) Pervasive functional translation of noncanonical human open reading frames. *Science*, **367**, 1140-1146.

32. Rodriguez,C.M., Chun,S.Y., Mills,R.E. et al. (2019) Translation of upstream open reading frames in a model of neuronal differentiation. *BMC Genomics*, **20**, 391.

33. Sharipov,R.N., Yevshin,I.S., Kondrakhin,Y.V. et al. (2014) RiboSeqDB – a repository of selected human and mouse ribosome footprint and RNA-seq data. *Virtual Biol.*, **1**, 37-46.

34. Evans,V.C., Barker,G., Heesom,K.J. et al. (2012) De novo derivation of proteomes from transcriptomes for transcript and protein identification. *Nat. Methods*, **9**, 1207-1211.

35. Liu,W., Xiang,L., Zheng,T. et al. (2018) TranslatomeDB: a comprehensive database and cloud-based analysis platform for translatome sequencing data. *Nucleic Acids Res.*, **46**, D206-D212.

36. Lee,S., Liu,B., Lee,S. et al. (2012) Global mapping of translation initiation sites in mammalian cells at single-nucleotide resolution. *Proc. Natl. Acad. Sci. U.S.A.*, **109**, E2424-E2432.

37. Olexiouk,V., Van Criekinge,W. and Menschaert,G. (2018) An update on sORFs.org: a repository of small ORFs identified by ribosome profiling. *Nucleic Acids Res.*, **46**, D497-D502.
38. Yates, A.D., Achuthan, P., Akanni, W. *et al.* (2020) Ensembl 2020. *Nucleic Acids Res.*, **48**, D682–D688.

39. Hansen, K.D., Sabunciyan, S., Langmead, B. *et al.* (2014) AnnotationHub: large-scale hypomethylated blocks associated with Epstein-Barr virus-induced B-cell immortalization. *Genome Res.*, **24**, 177–184.

40. Diehl, A.D., Meehan, T.F., Bradford, Y.M. *et al.* (2016) The cell ontology 2016: enhanced content, modularization, and ontology interoperability. *J Biomed Semantics*, **7**, 44.

41. Sarntivijai, S., Lin, Y., Xiang, Z. *et al.* (2014) CLO: the cell line ontology. *J Biomed Semantics*, **5**, 37.

42. Chang, A., Schomburg, I., Placzek, S. *et al.* (2015) BRENDA in 2015: exciting developments in its 25th year of existence. *Nucleic Acids Res.*, **43**, D439–D446.

43. Welter, D., Osumi-Sutherland, D. and Jupp, S. (2018) Human cell atlas ontology. CEUR-WS.org, Vol. 2285, p. 1–2.

44. Golbreich, C., Grosjean, J. and Darmoni, S.J. (2013) The foundational model of anatomy in OWL 2 and its use. *Artif. Intell. Med.*, **57**, 119–132.

45. Bandrowski, A., Brinkman, R., Brochhausen, M. *et al.* (2016) The ontology for biomedical investigations. *PLoS One*, **11**, e0154556.

46. Sioutos, N., de Coronado, S., Haber, M.W. *et al.* (2007) NCI Thesaurus: a semantic model integrating cancer-related clinical and molecular information. *J Biomed Inform*, **40**, 30–43.

47. Malone, J., Holloway, E., Adamusiak, T. *et al.* (2010) Modeling sample variables with an experimental factor ontology. *Bioinformatics*, **26**, 1112–1118.

48. Abeyruwan, S., Vempati, U.D., Küçük-mcginty, H. *et al.* (2014) Evolving BioAssay Ontology (BAO): modularization, integration and applications. *J Biomed Semantics*, **5**, S5.

49. Huang, J., Gutierrez, F., Strachan, H.J. *et al.* (2016) OmniSearch: a semantic search system based on the Ontology for MicroRNA Target (OMIT) for microRNA-target gene interaction data. *J Biomed Semantics*, **7**, 25.

50. Jupp, S., Burdett, T., Malone, J. *et al.* (2015) A new ontology lookup service at EMBL-EBI. *SWAT4LS*, **2**, 118–119.

51. Hernández, G., Osnaya, V.G. and Pérez-Martínez, X. (2019) Conservation and variability of the AUG initiation codon context in eukaryotes. *Trends Biochem. Sci.*, **44**, 1009–1021.

52. Han, J., Back, S.H., Hur, J. *et al.* (2013) ER-stress-induced transcriptional regulation increases protein synthesis leading to cell death. *Nat. Cell Biol.*, **15**, 481–490.

53. Guan, B.-J., van Hoef, V., Jobava, R. *et al.* (2017) A unique ISR program determines cellular responses to chronic stress. *Mol. Cell*, **68**, 885–900.e6.

54. Raudvere, U., Kolberg, L., Kuzmin, I. *et al.* (2019) g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). *Nucleic Acids Res.*, **47**, W191–W198.

55. Ashburner, M., Ball, C.A., Blake, J.A. *et al.* (2000) Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium. *Nat. Genet.*, **25**, 25–29.

56. The Gene Ontology Consortium. (2019) The Gene Ontology Resource: 20 years and still GOing strong. *Nucleic Acids Res.*, **47**, D330–D338.

57. Kent, W.J., Sugnet, C.W., Furey, T.S. *et al.* (2002) The human genome browser at UCSC. *Genome Res.*, **12**, 996–1006.

58. Vattem, K.M. and Wek, R.C. (2004) Reinitiation involving upstream ORFs regulates ATF4 mRNA translation in mammalian cells. *Proc. Natl. Acad. Sci. U.S.A.*, **101**, 11269–11274.

59. Brunet, M.A., Brunelle, M., Lucier, J.-F. *et al.* (2019) OpenProt: a more comprehensive guide to explore eukaryotic coding potential and proteomes. *Nucleic Acids Res.*, **47**, D403–D410.