# Characterization of accessory genes in coronavirus genomes

Christian Jean Michel, Claudine Mayer, Olivier Poch, Julie D. Thompson

► To cite this version:

Christian Jean Michel, Claudine Mayer, Olivier Poch, Julie D. Thompson. Characterization of accessory genes in coronavirus genomes. Virology Journal, 2020, 17 (1), pp.131. 10.1186/s12985-020-01402-1. hal-03464220

HAL Id: hal-03464220

https://hal.science/hal-03464220

Submitted on 3 Dec 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH

Open Access

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

Check for
updates

# Characterization of accessory genes in coronavirus genomes

Christian Jean Michel¹, Claudine Mayer¹,²,³, Olivier Poch¹ and Julie Dawn Thompson¹*

**Abstract**

**Background:** The Covid19 infection is caused by the SARS-CoV-2 virus, a novel member of the coronavirus (CoV) family. CoV genomes code for a ORF1a / ORF1ab polyprotein and four structural proteins widely studied as major drug targets. The genomes also contain a variable number of open reading frames (ORFs) coding for accessory proteins that are not essential for virus replication, but appear to have a role in pathogenesis. The accessory proteins have been less well characterized and are difficult to predict by classical bioinformatics methods.

**Methods:** We propose a computational tool GOFIX to characterize potential ORFs in virus genomes. In particular, ORF coding potential is estimated by searching for enrichment in motifs of the X circular code, that is known to be over-represented in the reading frames of viral genes.

**Results:** We applied GOFIX to study the SARS-CoV-2 and related genomes including SARS-CoV and SARS-like viruses from bat, civet and pangolin hosts, focusing on the accessory proteins. Our analysis provides evidence supporting the presence of overlapping ORFs 7b, 9b and 9c in all the genomes and thus helps to resolve some differences in current genome annotations. In contrast, we predict that ORF3b is not functional in all genomes. Novel putative ORFs were also predicted, including a truncated form of the ORF10 previously identified in SARS-CoV-2 and a little known ORF overlapping the Spike protein in Civet-CoV and SARS-CoV.

**Conclusions:** Our findings contribute to characterizing sequence properties of accessory genes of SARS coronaviruses, and especially the newly acquired genes making use of overlapping reading frames.

**Keywords:** COVID-19, SARS-CoV-2, SARS-CoV, Coronavirus, Accessory genes, ORF prediction, Circular code motifs

**Background**

Coronaviruses (CoVs) cause respiratory and intestinal infections in animals and humans [1]. They were not considered to be highly pathogenic to humans until the last two decades, which have seen three outbreaks of highly transmissible and pathogenic coronaviruses, including SARS-CoV (severe acute respiratory syndrome coronavirus), MERS-CoV (Middle East respiratory syndrome coronavirus), and SARS-CoV-2 (which causes the disease COVID-19). Other human coronaviruses (such as HCoV-NL63, HCoV-229E, HCoV-OC43 or HKU1)

generally induce only mild upper respiratory diseases in immunocompetent hosts, although some may cause severe infections in infants, young children and elderly individuals [1].

Extensive studies of human coronaviruses have led to a better understanding of coronavirus biology. Coronaviruses belong to the family *Coronaviridae* in the order nidovirales. Whereas MERS-CoV is a member of the *Merbecovirus* subgenus, phylogenetic analyses indicated that SARS-CoV-2 clusters with SARS-CoV in the *Sarbecovirus* subgenus [2]. All human coronaviruses are considered to have animal origins. SARS-CoV, MERS-CoV and SARS-CoV-2 are assumed to have originated in bats [1]. It is widely believed that SARS-CoV and SARS-CoV-2 were transmitted directly to humans from market

<footnote>* Correspondence: thompson@unistra.fr
¹Laboratoire ICube, Department of Computer Science, CNRS, University of Strasbourg, F-67412 Strasbourg, France
Full list of author information is available at the end of the article</footnote>
civets and pangolin, respectively, based on the sequence
analyses of CoV isolated from these animals and from
infected patients.

All members of the coronavirus family are enveloped
viruses that possess long positive-sense, single-stranded
RNA genomes ranging in size from 27 to 33 kb. The cor-
onavirus genomes encode five major open reading
frames (ORFs), including a 5' frameshifted polyprotein
(ORF1a/ORF1ab) and four canonical 3' structural pro-
teins, namely the spike (S), envelope (E), membrane (M)
and nucleocapsid (N) proteins, which are common to all
coronaviruses [3]. In addition, a number of subgroup-
specific accessory genes are found interspersed among,
or even overlapping, the structural genes. Overlapping
genes originate by a mechanism of overprinting, in
which nucleotide substitutions in a pre-existing frame
induce the expression of a novel protein in an alternative
frame. The accessory proteins in coronaviruses vary in
number, location and size in the different viral sub-
groups, and are thought to contain additional functions
that are often not required for virus replication, but are
involved in pathogenicity in the natural host [4, 5].

In the face of the ongoing COVID-19 pandemic, ex-
tensive worldwide research efforts have focused on iden-
tifying coronavirus genetic variation and selection [6–8],
in order to understand the emergence of host/tissue
specificities and to help develop efficient prevention and
treatment strategies. These studies are complemented by
structural genomics [9–11], as well as transcriptomics
[12] and interactomics studies [13] of the structural and
putative accessory proteins.

However, there have been less studies of accessory
proteins, for two main reasons [14]. First, accessory pro-
teins are often not essential for viral replication or struc-
ture, but play a role in viral pathogenicity or spread by
modulating the host interferon signaling pathways for
example. This has led to some contradictory experimen-
tal results concerning the presence or functionality of
accessory proteins. For instance, in a recent experiment
[13] to characterize SARS-CoV-2 gene functions, 9 pre-
dicted accessory protein ORFs (3a, 3b, 6, 7a, 7b, 8, 9b,
9c, 10) were codon optimized and successfully expressed
in human cells, with the exception of ORF3b. However,
another recent study using DNA nanoball sequencing
[12] concluded that the SARS-CoV-2 expresses only five
canonical accessory ORFs (3a, 6, 7a, 7b, 8).

Second, bioinformatics approaches for the prediction
of accessory proteins are challenged by their complex
nature as short, overlapping ORFs. Such proteins are
known to have biased amino acid sequences compared
to non-overlapping proteins [15]. In addition, the
homology-based approaches widely used to predict
ORFs in genomes are less useful here, because many
accessory proteins are lineage- or subgroup-specific.

Thus, many state of the art viral genome annotation sys-
tems, such as Vgas [16], only predict overlapping pro-
teins if homology information is available. Other
methods have been developed dedicated specifically to
the ab initio prediction of overlapping genes, for ex-
ample based on multiple sequence alignments and statis-
tical estimates of the degree of variability at synonymous
sites [17] or sequence simulations and calculation of ex-
pected ORF lengths [18].

Here, we propose a computational tool GOFIX (Gene
prediction by Open reading Frame Identification using X
motifs) to predict potential ORFs in virus genomes.
Using a complete viral genome as input, GOFIX first lo-
cates all potential ORFs, defined as a region delineated
by start and stop codons. In order to predict functional
ORFs, GOFIX calculates the enrichment of the ORFs in
X motifs, i.e. motifs of the X circular code [19], a set of
20 codons that are over-represented in the reading
frames of genes from a wide range of organisms. For ex-
ample, in a study of 299,401 genes from 5217 viruses
[20] including double stranded and single stranded DNA
and RNA viruses, codons of the X circular code were
found to occur preferentially in the reading frame of the
genes. This is an important property of viral genes, since
it has been suggested that X motifs at different locations
in a gene may assist the ribosome to maintain and
synchronize the reading frame [21]. An initial evaluation
test of the GOFIX method on a large set of 80 virus ge-
nomes [15] showed that it achieves high sensitivity and
specificity for the prediction of experimentally verified
overlapping proteins (manuscript in preparation). A
major advantage of our approach is that it requires only
the sequence of the studied genome and does not rely
on any homology information. This allows us to detect
novel ORFs that are specific to a given lineage.

We applied GOFIX to study the SARS-CoV-2 genome
and related SARS genomes, with a main focus on the
accessory proteins. Using the extensive experimental
data concerning the SARS-CoV genome and the
expressed ORFs, we first show that the reading frames of
the SARS-CoV ORFs are enriched in X motifs, including
most of the overlapping accessory proteins. Exceptions
include SARS-CoV ORF3b and ORF8b which may not
be functional. Then, we use GOFIX to predict and com-
pare putative genes in related genomes of SARS-like vi-
ruses from bat, civet and pangolin hosts as well as
human SARS-CoV-2.

Methods

Genome sequences

Viral genome sequences were downloaded from the
Genbank database, as shown in Table 1. The Genbank
reference genomes were used as representative genomes
for SARS-CoV and SARS-CoV-2. For the Bat-CoV,
<table><caption>**Table 1** Genome sequences selected for the current study. Note that the SARS-CoV strain hTor02 is from humans infected during the middle and late phases of the SARS epidemic of 2013, and has a deletion of 29 nucleotides in the region of ORF8</caption><thead><tr><th></th><th>Description</th><th>Genbank accession number</th></tr></thead><tbody><tr><td>Bat-CoV</td><td>Bat SARS-like coronavirus isolate As6526</td><td>KY417142</td></tr><tr><td>Civet-CoV</td><td>Civet SARS coronavirus civet007</td><td>AY572034</td></tr><tr><td>SARS-CoV</td><td>Human severe acute respiratory syndrome-related coronavirus strain hTor02</td><td>NC_004718</td></tr><tr><td>Pangolin-CoV</td><td>Pangolin coronavirus isolate PCoV_GX-P2V</td><td>MT072864</td></tr><tr><td>SARS-CoV-2</td><td>Human severe acute respiratory syndrome coronavirus 2 isolate Wuhan-Hu-1</td><td>MT072688</td></tr></tbody></table>

Civet-CoV and Pangolin-CoV genomes, we selected well annotated Genbank entries having the highest number of annotated ORFs. All CDS annotations were extracted from the Genbank files, and ORF names were standardized according to the SARS-CoV-2 nomenclature (Table 2).

## Definition of X motif enrichment (XME) scores

The X circular code contains the following 20 codons $X = \{AAC, AAT, ACC, ATC, ATT, CAG, CTC, CTG, GAA, GAC, GAG, GAT, GCC, GGC, GGT, GTA, GTC, GTT, TAC, TTC\}$ (1) and has several strong mathematical properties [19]. In particular, it is self-complementary, i.e. 10 trinucleotides of X are complementary to the other 10 trinucleotides of X, and it is a circular code. A circular code is defined as a set of words such that any motif obtained from this set, allows to retrieve, maintain and synchronize the reading frame.

An X motif *m* is defined as a word containing only codons from the X circular code (1) with length $|m| \ge 3$ codons and cardinality (i.e. number of unique codons) $c \ge 2$ codons. The minimal length $|m| = 3$ codons was chosen based on a previous study showing that the probability of retrieving the reading frame with an X motif of at least 3 codons is 99.9% [22]. The class of X motifs with cardinality $c < 2$ are excluded here because they are mostly associated with the “pure” trinucleotide repeats often found in non-coding regions of genomes [23].

The total length $XL_f$ of all X motifs $m_f$ of nucleotide length $|m_f|$ in a frame *f* (the reading frame or one of the 2 shifted frames) of a nucleotide sequence *s* is defined as:

$$XL_f = \sum_{m_f \in s} |m_f|.$$

Then the X motif enrichment $XME_f$ in a frame *f* of a sequence *s* of nucleotide length *l* is defined as:

$$XME_f = \frac{100}{l_f} XL_f$$

where for non-overlapping ORFs: $l_f = l$, and for overlapping ORFs: $l_f = l - XL_g$ where $XL_g$ is the total length of all X motifs in the overlapped frame *g*.

Finally, for an ORF of length *l* and associated with a reading frame *f*, the X motif enrichment score $XME$ is defined as:

$$XME = XME_f.$$

## GOFIX method

The GOFIX method will be described in detail in a separate manuscript. Briefly, the method consists of two main steps:

(i.) Identification of all potential ORFs. Using the complete genome sequences as input, all potential ORFs in the positive sense are located, defined as a sequence region starting with a start codon (AUG) and ending with a stop codon (UAA, UAG, UGA). For a given region, if alternative start codons are found, the longest ORF is selected. In this study, we selected all ORFs having a minimum length of 120 nucleotides (40 amino acids).

(ii.) Calculation of X motif enrichment scores. For each potential ORF, all X motifs in the nucleotide sequence are identified in the three positive sense frames *f* using the computational method described in [24]. For each identified potential ORF, the X motif enrichment ($XME_f$ and $XME$) scores are calculated as defined above. We set the threshold for prediction of a functional ORF to be $XME \ge 5$ (i.e. XME score in the reading frame $\ge 5$), based on our benchmark studies (data not shown) of experimentally validated ORFs in a large set of 80 genomes [15] covering a wide range of viruses (including single-stranded and double-stranded DNA viruses and single-stranded and double-stranded RNA viruses). To avoid any statistical bias, coronavirus genome sequences were excluded from the benchmark.

## Results

### Initial study of SARS-CoV reference genome

We first analyzed the complete genome of the well-studied SARS-CoV and plotted the X motif enrichment ($XME_f$) scores calculated in a sliding window of 150
<table><caption>**Table 2** CDS annotations extracted from Genbank, with ORF names standardized according to the SARS-CoV-2 nomenclature</caption><thead><tr><th>Name</th><th>Start</th><th>Stop</th><th>Length</th></tr></thead><tbody><tr><td colspan="4"><strong>Bat-CoV</strong></td></tr><tr><td>ORF1a<sup>a</sup></td><td>265</td><td>13,398</td><td>13,134</td></tr><tr><td>ORF1b<sup>a</sup></td><td>13,398</td><td>21,485</td><td>8086</td></tr><tr><td>S</td><td>21,492</td><td>25,217</td><td>3726</td></tr><tr><td>ORF3a</td><td>25,227</td><td>26,051</td><td>825</td></tr><tr><td>ORF3b</td><td>25,648</td><td>25,992</td><td>345</td></tr><tr><td>E</td><td>26,076</td><td>26,306</td><td>231</td></tr><tr><td>M</td><td>26,357</td><td>27,022</td><td>666</td></tr><tr><td>ORF6</td><td>27,033</td><td>27,224</td><td>192</td></tr><tr><td>ORF7a</td><td>27,232</td><td>27,600</td><td>369</td></tr><tr><td>ORF7b</td><td>27,597</td><td>27,731</td><td>135</td></tr><tr><td>ORF8</td><td>27,738</td><td>28,103</td><td>366</td></tr><tr><td>N</td><td>28,118</td><td>29,386</td><td>1269</td></tr><tr><td colspan="4"><strong>Pangolin-CoV</strong></td></tr><tr><td>ORF1a</td><td>249</td><td>13,427</td><td>13,179</td></tr><tr><td>ORF1b</td><td>13,427</td><td>21,514</td><td>8086</td></tr><tr><td>S</td><td>21,522</td><td>25,331</td><td>3810</td></tr><tr><td>ORF3a</td><td>25,341</td><td>26,168</td><td>828</td></tr><tr><td>E</td><td>26,193</td><td>26,420</td><td>228</td></tr><tr><td>M</td><td>26,468</td><td>27,136</td><td>669</td></tr><tr><td>6</td><td>27,147</td><td>27,332</td><td>186</td></tr><tr><td>7a</td><td>27,339</td><td>27,704</td><td>366</td></tr><tr><td>7b</td><td>27,701</td><td>27,832</td><td>132</td></tr><tr><td>8</td><td>27,839</td><td>28,202</td><td>366</td></tr><tr><td>N</td><td>28,218</td><td>29,471</td><td>1254</td></tr><tr><td colspan="4"><strong>Civet-CoV</strong></td></tr><tr><td>ORF1a</td><td>239</td><td>13,366</td><td>13,128</td></tr><tr><td>ORF1b</td><td>13,366</td><td>21,459</td><td>8092</td></tr><tr><td>S</td><td>21,466</td><td>25,233</td><td>3768</td></tr><tr><td>ORF3a</td><td>25,242</td><td>26,066</td><td>825</td></tr><tr><td>ORF3b</td><td>25,663</td><td>26,127</td><td>465</td></tr><tr><td>E</td><td>26,091</td><td>26,321</td><td>231</td></tr><tr><td>M</td><td>26,372</td><td>27,037</td><td>666</td></tr><tr><td>ORF6</td><td>27,048</td><td>27,239</td><td>192</td></tr><tr><td>ORF7a</td><td>27,247</td><td>27,615</td><td>369</td></tr><tr><td>ORF7b</td><td>27,612</td><td>27,746</td><td>135</td></tr><tr><td>ORF8</td><td>27,753</td><td>28,121</td><td>369</td></tr><tr><td>N</td><td>28,123</td><td>29,391</td><td>1269</td></tr><tr><td>ORF9b</td><td>28,133</td><td>28,429</td><td>297</td></tr><tr><td>ORF9c</td><td>28,586</td><td>28,798</td><td>213</td></tr><tr><td colspan="4"><strong>SARS-CoV-2</strong></td></tr><tr><td>ORF1a</td><td>251</td><td>13,453</td><td>13,203</td></tr><tr><td>ORF1b</td><td>13,453</td><td>21,538</td><td>8086</td></tr><tr><td>S</td><td>21,521</td><td>25,369</td><td>3849</td></tr></tbody></table>
<table><caption>**Table 2** CDS annotations extracted from Genbank, with ORF names standardized according to the SARS-CoV-2 nomenclature
(Continued)</caption><thead><tr><th>Name</th><th>Start</th><th>Stop</th><th>Length</th></tr></thead><tbody><tr><td>ORF3a</td><td>25,378</td><td>26,205</td><td>828</td></tr><tr><td>ORF3b<sup>c</sup></td><td>25,509</td><td>25,680</td><td>172</td></tr><tr><td>E</td><td>26,230</td><td>26,457</td><td>228</td></tr><tr><td>M</td><td>26,508</td><td>27,176</td><td>669</td></tr><tr><td>ORF6</td><td>27,187</td><td>27,372</td><td>186</td></tr><tr><td>ORF7a</td><td>27,379</td><td>27,744</td><td>366</td></tr><tr><td>ORF7b<sup>c</sup></td><td>27,741</td><td>27,872</td><td>130</td></tr><tr><td>ORF8</td><td>27,879</td><td>28,244</td><td>366</td></tr><tr><td>ORFN</td><td>28,259</td><td>29,518</td><td>1260</td></tr><tr><td>ORF9b<sup>c</sup></td><td>28,269</td><td>28,562</td><td>294</td></tr><tr><td>ORF9c<sup>c</sup></td><td>28,719</td><td>28,940</td><td>222</td></tr><tr><td>ORF10</td><td>29,543</td><td>29,659</td><td>117</td></tr><tr><td>SARS-CoV</td><td></td><td></td><td></td></tr><tr><td>ORF1a</td><td>265</td><td>13,398</td><td>13,134</td></tr><tr><td>ORF1b</td><td>13,398</td><td>21,485</td><td>8086</td></tr><tr><td>S</td><td>21,492</td><td>25,259</td><td>3768</td></tr><tr><td>ORF3a</td><td>25,268</td><td>26,092</td><td>825</td></tr><tr><td>ORF3b</td><td>25,689</td><td>26,153</td><td>465</td></tr><tr><td>E</td><td>26,117</td><td>26,347</td><td>231</td></tr><tr><td>M</td><td>26,398</td><td>27,063</td><td>666</td></tr><tr><td>ORF6</td><td>27,074</td><td>27,265</td><td>192</td></tr><tr><td>ORF7a</td><td>27,273</td><td>27,641</td><td>369</td></tr><tr><td>ORF7b</td><td>27,638</td><td>27,772</td><td>135</td></tr><tr><td>ORF8a</td><td>27,779</td><td>27,898</td><td>120</td></tr><tr><td>ORF8b</td><td>27,864</td><td>28,118</td><td>255</td></tr><tr><td>N</td><td>28,120</td><td>29,388</td><td>1269</td></tr><tr><td>ORF9b</td><td>28,130</td><td>28,426</td><td>297</td></tr><tr><td>ORF9c<sup>b</sup></td><td>28,583</td><td>28,793</td><td>211</td></tr></tbody></table>

nucleotides for each of the three positive sense frames (Fig. 1). We then mapped the ORF1ab, the four structural proteins (S, E, M, N), and the nine generally accepted accessory genes (3a, 3b, 6, 7a, 7b, 8a, 8b, 9b, 9c) to the X enrichment plot.

We observe a tendency for the reading frames of the SARS-CoV ORFs to be enriched in X motifs. For example, ORF1ab is the longest ORF, encoding a polyprotein, which is translated by a – 1 programmed ribosomal frameshift at position 13,398. Sequences upstream and downstream of the frameshift are enriched in X motifs in the corresponding reading frame (green and yellow plots respectively in Fig. 1a). Other ORFs enriched in X motifs in the reading frame include the S protein (yellow plot in Fig. 1b) and

the E and M proteins (blue and green plots respectively in Fig. 1c). The S, E and M ORFs are conserved in all coronavirus genomes and code for structural proteins that together create the viral envelope.

The case of overlapping ORFs is more complex. For example, the last structural protein coded by the N ORF is overlapped by two accessory genes: ORF9b and ORF9c. The sequence regions containing the overlapping ORFs are characterized by an enrichment in X motifs in the 2 frames (green and blue plots in Fig. 1c).

## Characterization of known accessory genes in SARS-CoV

The SARS-CoV genome is known to contain four structural proteins and nine accessory proteins, namely ORFs

<footnote><sup>a</sup> For convenience, ORF1ab is split into 2 regions corresponding the ORF1ab gene regions upstream and downstream of the frameshift</footnote>

<footnote><sup>b</sup> SARS-CoV annotation for ORF9c was propagated from Genbank entry AY274119: SARS-CoV isolate Tor2, where it is annotated as ORF14</footnote>

<footnote><sup>c</sup> SARS-CoV-2 annotations for ORF3b, ORF7b, ORF9b and ORF9c were propagated from Genbank entry MN985325: Severe acute respiratory syndrome coronavirus 2 isolate 2019-nCoV/USA-WA1/2020</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** X motif enrichment (XME<sub>f</sub>) scores in the three frames *f* = 0, 1 and 2 (green, blue, yellow respectively) of the SARS-CoV genome, using a sliding window of length 150 nucleotides. Genomic organization of known ORFs is shown underneath the plots. **a** Polyprotein gene ORF1ab. **b** Spike protein. **c** C-terminal structural and accessory proteins. The colors used in the enrichment plot and in the boxes representing ORFs (green, blue, yellow) indicate the three frames 0,1 and 2 respectively</figcaption></figure>

3a, 3b, 6, 7a, 7b, 8a, 8b, 9b and 9c. To verify that our approach can predict the accessory genes in coronavirus genomes, we used GOFIX to identify all potential ORFs in the complete SARS-CoV genome and calculate their *X* enrichment. Figure 2 shows the *X* motif enrichment (XME<sub>f</sub>) scores calculated by GOFIX for the identified ORFs in the 3' terminal region of the SARS-CoV genome.

The overall performance of GOFIX is shown in Table 3. Initially, GOFIX found 25 potential ORFs (delineated by start and stop codons) in the 3' region (21,492–29,751) of SARS-CoV. Twelve of these 25 potential ORFs were predicted to be non-functional (see Methods), including 10 unknown ORFs mostly overlapping the S protein. Two previously annotated ORFs were also predicted to be non-

<figure><img src="image_4.png" /><figcaption>**Fig. 2** XME<sub>f</sub> scores calculated by GOFIX for potential ORFs in the 3' terminal region of the SARS-CoV genome, in the three frames *f* = 0, 1 and 2 (green, blue, yellow respectively). For clarity, only Genbank annotated ORFs or new ORFs predicted in this work are shown. The red line represents the threshold value XME = XME<sub>f</sub> = 5 (where *f* is the reading frame) for the prediction of a functional ORF. Known ORFs are indicated below the histogram using the color corresponding to the ORF reading frame. Known ORFs not predicted to be functional by GOFIX are outlined in red. Novel ORFs predicted by GOFIX are outlined in blue</figcaption></figure>
<table><caption>**Table 3** Prediction performance of the GOFIX method on the set of known ORFs in the SARS-CoV genome</caption><thead><tr><th></th><th>Predicted: YES</th><th>Predicted: NO</th><th>Total</th></tr></thead><tbody><tr><td>Known ORF</td><td>11</td><td>2</td><td>13</td></tr><tr><td>Unknown ORF</td><td>2</td><td>10</td><td>12</td></tr><tr><td>Total</td><td>13</td><td>12</td><td>25</td></tr><tr><td colspan="2">Sensitivity = 0.85</td><td colspan="2">Specificity = 0.83</td></tr></tbody></table>

functional, namely ORF3b (XME = 1.9) and ORF8b (XME = 0.0) that are discussed in detail below.

GOFIX predicts that 13 of the 25 potential ORFs are functional, with XME scores (in the reading frame) > 5. These include 11 previously annotated ORFs, namely S, 3a, E, M, 6, 7a, 7b, 8a, N, 9b, 9c. Two novel ORFs are also predicted by the GOFIX method: ORF10 (XME = 15.8) is located downstream of the N gene (29,415–29,496) and a new ORF we called ORFSa (XME = 7.6) that overlaps the S gene (22,732–22,928). These novel ORFs are discussed in more detail below. It should also be noted that some ORFs have larger XME scores in the shifted frames than in the reading frame. This is often linked to overlapping ORFs where proteins are coded in more than one frame, for example ORFs 3a and 3b, or ORFs 9b, 9c that overlap the N ORF.

## Comparative analyses of accessory proteins in coronavirus genomes

Having evaluated the GOFIX method on the well-studied SARS-CoV genome, we then used it to characterize and compare the accessory proteins in

representative strains of five coronavirus genera, including SARS-CoV, SARS-CoV-2 and three viruses from animal hosts with SARS-CoV-like infections. Bat is considered to be the most likely host origin of SARS-CoV and SARS-CoV-2. It is generally considered that transmission to humans occurred via an intermediate host. For SARS-CoV, civets probably acted as the intermediate host, while pangolin has been proposed as the intermediate host in SARS-CoV-2 animal-to-human transmission [25]. For each of the five genomes, we used GOFIX to predict all potential ORFs in the complete genomes and calculated the X motif enrichment (XME) scores for each ORF. Figure 3 gives an overview of the predicted ORFs in each genome, confirming for example that the structural proteins S, E, M and N, as well as the accessory proteins ORF6, ORF7a and ORF7b are conserved and have XME scores above the defined threshold XME = 5. However, important differences in XME scores are observed for the remaining accessory protein ORFs.

## ORF3b may not code for a functional protein in all CoVs

ORF3a codes for the largest accessory protein that comprises 274–275 amino acids (Fig. 4). In SARS-CoV, ORF3a is not required for virus replication, but contributes to pathogenesis by mediating trafficking of Spike (S protein) [4]. It is efficiently expressed on the cell surface, and was easily detected in a majority of SARS patients. The XME scores for ORF3a in all the genomes range from 13.8–19.3, i.e. almost 3 times greater than the defined threshold for functional ORFs.

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Prediction of ORFs in representative SARS-like coronavirus genomes. A schema is provided for each genome, showing the Genbank annotated ORFs and new ORFs predicted in this work. The numbers in the tables below each schema indicate the XME scores in the reading frame of each ORF. Genbank annotated ORFs that are not predicted to be functional by the GOFIX method are highlighted in red. Novel ORFs predicted by GOFIX are shown in blue. ORFs with conflicting annotations in Genbank, but predicted by GOFIX are shown in brown. Note that ORF3b in Civet-CoV and SARS-CoV is not homologous to ORF3b in Pangolin-CoV and SARS-CoV-2</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. 4** a Schematic view of genome organization of ORF3a, ORF3b and E gene. b Multiple alignment of ORF3a, ORF3b sequences, with *X* motifs in the reading frame of ORF3a shown in blue. The start and stop codons of the overlapping ORF3b sequences (in the + 1 reading frame of ORF3a) are indicated by purple and red boxes respectively. *X* motifs in the reading frame of ORF3b are shown in green</figcaption></figure>

The ORF3b coding sequence overlaps the + 1 reading frame of ORF3a and sometimes extends beyond the start codon of the E gene. In SARS-CoV, it is proposed to antagonize interferon (IFN) function by modulating the activity of IFN regulatory factor 3 (IRF3) [26]. However, immunohistochemical analyses of tissue biopsies and/or autopsies of SARS-CoV-infected patients have failed to demonstrate the presence of ORF3b in vivo, and the presence of ORF3b in SARS-CoV-infected Vero E6 cells is the only evidence for the expression of this protein [27]. Furthermore, when mice are infected with mutant SARS-CoV lacking ORF3b, the deletion viruses grow to levels similar to those of wild-type virus, which demonstrates that SARS-CoV is able to inhibit the host IFN response without the 3b gene [28].

Bat-Cov and Civet-CoV also present ORF3b overlapping the 3' region of ORF3a (start codon at nt 422), although the sequence of Bat-CoV ORF3b is shorter having a stop codon within the ORF3a sequence (nt 764). We observe a single *X* motif in the ORF3b reading frame of length 9 nucleotides (563–571), resulting in low XME scores of 2.6, 1.9 and 1.9 respectively for Bat-CoV, Civet-CoV and SARS-CoV ORF3b. We thus predict that ORF3b is not functional in these strains. This ORF is not predicted to be present in Pangolin-CoV or SARS-CoV-2 due to the introduction of a new stop codon

(indicated by *** in Fig. 4) and the loss of the *X* motif in the + 1 reading frame.

However, a completely different ORF is identified in the Pangolin-CoV and SARS-CoV-2 sequences, overlapping the 5' region of ORF3a (132–305). This ORF is not annotated in the SARS-CoV-2 reference genome (MT072688), but is annotated as ORF3b in the genome of another SARS-CoV-2 strain isolated from the first U.S. case of COVID-19 (MN985325). The Pangolin-CoV ORF3b sequence contains one *X* motif in the reading frame of length 9 nucleotides (183–191), with an XME score of 4.6. However, the *X* motif is lost in the SARS-CoV-2 genome, in agreement with recent ORF expression data [13].

## ORF8: a rapidly evolving region of SARS-CoV genomes

Previously shown to be a recombination hotspot, ORF8 is one of the most rapidly evolving regions among SARS-CoV genomes [29]. Furthermore, the evolution of ORF8 is supposed to play a significant role in adaptation to the human host following interspecies transmission and virus replicative efficiency [30].

In SARS-CoV isolated from bats and civets (as well as early human isolates of the SARS-CoV outbreak in 2003: data not shown), ORF8 encodes a single protein of length 122 amino acids (Fig. 5). However, in SARS-CoV
<figure><img src="image_7.png" /></figure>

isolated from humans during the peak of the epidemic, there is a 29-nt deletion in the middle of ORF8, resulting in the splitting of ORF8 into two smaller ORFs, namely ORF8a and ORF8b [31]. ORF8a and ORF8b encode a 39 amino acid and 84 amino acid polypeptide, respectively. The XME scores in these ORFs are in line with the known experimental evidence concerning their functions. ORF8a has an XME score of 15.3 in SARS-CoV and anti-p8a antibodies were identified in some patients with SARS [32]. In contrast, ORF8b has no X motifs in the reading frame, and protein 8b was not detected in SARS-CoV-infected Vero E6 cells [31].

It is interesting to note that although Civet-CoV has a full-length ORF8, it has a low XME score (XME = 4.9) compared to Bat-CoV (XME = 9.9). Thus, it is tempting to suggest that the loss of X motifs in transmission of the virus from bats to civets is somehow linked to the loss of ORF8 in the transmission from civets to humans. Both Pangolin-CoV and most SARS-CoV-2 strains contain the full length ORF8, with XME scores of 23.1 and 12.4 respectively. However, a 382-nt deletion has been reported recently covering almost the entire ORF8 of SARS-CoV-2 obtained from eight hospitalized patients in Singapore, that has been hypothesized to lead to an attenuated phenotype of SARS-CoV-2 [33].

## Characterization of ORFs overlapping the N gene

The annotation of functional ORFs overlapping the N gene is variable in the different genomes studied here. In SARS-CoV, only ORF9b has been observed to be translated, probably via a ribosomal leaky scanning mechanism and may have a function during virus assembly [30, 34]. ORF9b limits host cell interferon responses by targeting the mitochondrial-associated adaptor molecule (MAVS) signalosome. However, some SARS-CoV strains have an additional ORF9c, annotated as a hypothetical protein (e.g. Genbank:AY274119). For Bat-CoV and Pangolin-CoV, no overlapping genes are annotated in the corresponding Genbank entries. In contrast, the Civet-CoV genome is predicted to contain both overlapping genes, ORF9b and ORF9c. Similarly, the annotation of overlapping ORFs for SARS-CoV-2 is different depending on the strain: the reference strain has no overlapping ORFs of the N gene, while the U.S. strain has ORF9b and ORF9c (see Methods). ORF9c is described as a short polypeptide (70 amino acids) dispensable for viral replication, but there is no data yet providing evidence that the protein is expressed during SARS-CoV-2 infection.

Here, we predict that ORF9b and ORF9c are present in all genomes as overlapping ORFs within the N gene
(Fig. 6). Furthermore, Pangolin-CoV may also have an additional ORF, that we called ORF9d (XME = 12.7), in the 3' region of the N gene.

## Origin and evolution of ORF10

ORF10 is proposed as unique to SARS-CoV-2 [35] and codes for a peptide only 38 amino acids long. There is no data yet providing evidence that the protein is expressed during SARS-CoV-2 infection. Therefore, we wanted to investigate the potential origin of this protein. New proteins in viruses can originate from existing proteins acquired through horizontal gene transfer or through gene duplication for example, or can be generated de novo. To determine whether homologs of ORF10 are present in the other coronavirus genomes, we relaxed the GOFIX parameters used to predict functional ORFs, and set the minimum ORF length to 60

nucleotides. The predicted ORFs in the different genomes are shown in Fig. 7. The Pangolin-CoV genome contains a full-length ORF10 with XME = 10.4, compared to the SARS-CoV-2 ORF10 with XME = 20.2. A truncated version of ORF10 coding 26 amino acids is also detected in the Bat-CoV, Civet-CoV and SARS-CoV genomes, although this short ORF is probably not functional. We suggest that the ORF10 of SARS-CoV-2 thus evolved via the mutation of a stop codon (TAA) at nt 76 and the addition of a new X motif of length 15 nucleotides in the 3' region.

## Novel ORF overlapping the S gene

The GOFIX method predicts a novel ORF, that we called ORFSa, overlapping the RBD (Receptor Binding Domain) of the S (Spike) ORF in SARS-CoV (XME = 7.6) and Civet-CoV (XME = 7.6). ORFSa is found in the

<figure><img src="image_8.png" /><figcaption>**Fig. 6** a Schematic view of genome organization of ORF N, with overlapping genes ORF9b, 9c and the novel predicted 9d. **b** Multiple alignment of ORF N sequences, with X motifs in the reading frame of ORF N shown in blue, in ORF9b in green, in ORF9c in yellow. Start and stop codons of the overlapping genes are indicated by violet and red boxes, respectively. **c** The novel ORF9d predicted in Pangolin-CoV with X motifs in the reading frame shown in orange</figcaption></figure>
<table><tr><td>Bat-CoV</td><td>1 ATGGGCTAT</td><td>GTAAACGTTTTCGAATTCCGTTTACGATACATAGTCTACTCTTGTGCAGAATGAATTCTCGTAGCGTAAAC</td><td>80</td></tr><tr><td>Civet-CoV</td><td>1 ATGGGCTAT</td><td>GTAAACGTTTTCGAATTCCGTTTACGATACATAGTCTACTCTTGTGCAGAATGAATTCTCGTAAC</td><td>80</td></tr><tr><td>SARS-CoV</td><td>1 ATGGGCTAT</td><td>GTAAACGTTTTCGAATTCCGTTTACGATACATAGTCTACTCTTGTGCAGAATGAATTCTCGTAAC</td><td>80</td></tr><tr><td>Pangolin-CoV</td><td>1 ATGGGCTAT</td><td>GTAAACGTTTTCGAATTCCGTTTACGATACATAGTCTACTCTTGTGCAGAATGAATTCTCGTAGCTATAC</td><td>80</td></tr><tr><td>SARS-CoV-2</td><td>1 ATGGGCTAT</td><td>ATACGTTTTCGCTTTTCCGTTTACGATATATAGTCTACTCTTGTGCAGAATGAATTCTCGTAACTACAT</td><td>80</td></tr></table>

<table><caption>**Fig. 7** Multiple alignment of ORF10 sequences, with X motifs in the reading frame shown in blue. Stop codons are indicated by red boxes</caption><tr><td>Bat-CoV</td><td>81 AGCACAAGTAGGTTTAGTTAACTTTAATCTCACATAG</td><td>117</td></tr><tr><td>Civet-CoV</td><td>81 AGCACAAGTAGGTTTAGTTAACTTTAATCTCACATAG</td><td>117</td></tr><tr><td>SARS-CoV</td><td>81 AGCACAAGTAGGTTTAGTTAACTTTAATCTCACATAG</td><td>117</td></tr><tr><td>Pangolin-CoV</td><td>81 AGCACAAGTAGGTATAGTTAACTTTAATCTCACATAG</td><td>117</td></tr><tr><td>SARS-CoV-2</td><td>81 AGCACAACTAGATGTAGTTAACTTTAATCTCACATAG</td><td>117</td></tr></table>

+1 frame and codes for a protein with 64 amino acids, as shown in Fig. 8. As the ORFSa sequence was not present in the Bat-CoV reference genome, we also searched for the ORF in the genomes of other Bat-CoV strains, and found one occurrence (XME = 6.5) in the strain WIV16 (Genbank:KT444582) (Fig. 8), another bat coronavirus that is closely related to SARS-CoV [36].

To investigate whether the novel ORFSa might be a functional protein in SARS-CoV, we used BlastP to search the Genbank database for matches to viral proteins. A significant hit was obtained with a sequence identity of 100% to the protein AAR84376, described as “putative transmembrane protein 2d” from the genome of SARS coronavirus strain ZJ01 (AY28632). To further characterize this putative protein, the Phobius web site (phobius.sbc.su.se) was used to predict transmembrane

(TM) helices. Two potential TM helices of nearly twenty amino acids (residues 6–28 and 42–62) were predicted with a small inter-TM endodomain. Thus, this potential double-membrane spanning small protein might complement the set of already known SARS-CoV membrane proteins, namely the Spike (S), membrane (M) and envelope (E) proteins.

## Discussion

Coronaviruses are complex genomes with high plasticity in terms of gene content. This feature is thought to contribute to their ability to adapt to specific hosts and to facilitate host shifts [1]. It is therefore essential to characterize the coding potential of coronavirus genomes. Here, we used an ab initio approach to identify potential functional ORFs in the genomes of a set of

a.

<table><tr><td>Bat-CoV</td><td>1253</td><td>ACGACTTCACAGGCTGTGTTATTCCCTGGAATACT</td><td>GCTCAGCAAGATAAGGCTCAATAT</td><td>TATTATAGATCCTCTAGAAAGACAA</td><td>1336</td></tr><tr><td>Bat-CoV (WIV16)</td><td>1241</td><td>ATGATTTCACGGGTTGTGTCCTTGCTTGAATACTAGGAACATTGATGCTACTCAAACTCGTAATT</td><td>ATAATTATAAATATAGATCTCTCAGACATGGCA</td><td>1339</td></tr><tr><td>Civet-CoV</td><td>1241</td><td>ATGATTTCATGGGTTGTGTCCTTGCTTGAATACTAGGAACATTGATGCTACTTCAAACTCGTAATT</td><td>ATAATTATAAATATAGGTATCTTAGACATGGCA</td><td>1339</td></tr><tr><td>SARS-CoV</td><td>1241</td><td>ATGATTTCATGGGTTGTGTCCTTGCTTGAATACTAGGAACATTGATGCTACTTCAAACTCGTAATT</td><td>ATAATTATAAATATAGGTATCTTAGACATGGCA</td><td>1339</td></tr><tr><td>Pangolin-CoV</td><td>1280</td><td>ATGATTTTACTGGTTGTGTTATTGCTTGAACACTCAGTTAAGCAAGATGCTTGTACGTGGTGAATT</td><td>ATGCTGGTGAATTATGGTTATTGTATAGATTATTTAGAAAGTCTA</td><td>1378</td></tr><tr><td>SARS-CoV-2</td><td>1280</td><td>ATGATTTTACAGGCTGCGTTATAGCTTGAATTCTAACAATCTTGATTCTAAGCTTGGTCGAATT</td><td>AATTACCTGTATAGATTGTTTAGGAAGTCTA</td><td>1378</td></tr></table>

<table><tr><td>Bat-CoV</td><td>1337</td><td>AACTTAAACCTTTTGAGAGGACCTATCTTCA</td><td>GACGAAAATCGTGACGTACTCTTA</td><td>1393</td></tr><tr><td>Bat-CoV (WIV16)</td><td>1340</td><td>AGCTTAGCCTTTTGAGAGAGATATTCTAATGTGCCTTTCTCTCCTGATGGCAAACCTTGACCCACCT</td><td>GCTTTAATTCTTATTGGCCATTAA</td><td>1435</td></tr><tr><td>Civet-CoV</td><td>1340</td><td>AGCTTAGCCCTTTTGAGAGAGATATCTAATGTGCCTTTCTCTCTGATGGCAAACCTTGACCCACCT</td><td>GCTCCTAATTGTATTGGCCATTAA</td><td>1435</td></tr><tr><td>SARS-CoV</td><td>1340</td><td>AGCTTAGCCCTTTTGAGAGAGATATCTAATGTGCCTTTCTCTCCTGATGGCAAACCTTGACCCACCT</td><td>GCTCTTAATTGTATTGGCCATTAA</td><td>1435</td></tr><tr><td>Pangolin-CoV</td><td>1379</td><td>AGCTTAAACCTTTTGAGAGAGATATTCTCAACTGAATATACCAAGCCGGCAGCACACCCTGTAACGCTCAAGTTGGTCTAA</td><td>ATTGTTATTATCCTCTTTG</td><td>1477</td></tr><tr><td>SARS-CoV-2</td><td>1379</td><td>ATCTCAAACCTTTTGAGAGAGATATTTCAACTGAAATCTATCAGGCCGGTAGCACACCTTGTATCGTGTTAAGGTTTAATTGTTACTTTCCTTTAC</td><td>1477</td></tr></table>

b.

gat gat ttc atg ggt tgt gtc ctt gct tgg aat act agg aac att gat gct act tca act ggt aat
D D F M G C V L A W N T R N I D A T S T G N
M I S W V V S L L G I L G T L M L L Q L V I

tat aat tat aaa tat agg tat ctt aga cat ggc aag ctt agg ccc ttt gag aga gac ata tct aat
Y N Y K Y R Y L R H G K L R P F E R D I S N
I I I N I G I L D M A S L G P L R E T Y L M

gtg cct ttc tcc cct gat ggc aaa cct tgc acc cca cct gct ctt aat tgt tat tgg cca tta
V P F S P D G K P C T P P A L N C Y W P L
C L S P L M A N L A P H L L L I V I G H -

<caption>**Fig. 8** a Multiple alignment of ORFSa sequences, with X motifs in the reading frame of ORFS shown in blue and ORFSa in green. Start and stop codons of the overlapping genes are indicated by violet and red boxes, respectively. Bat-CoV (WIV16) sequence is from Genbank:KT444582. **b** Nucleotide and amino acid sequences of the novel ORF predicted to overlap the Spike protein in the genome of SARS-CoV. The nucleotide sequence segment (SARS-CoV-nt 22,732–22,926) encodes part (residues 414–478) of the RBD (residues 323–502) of the Spike protein (normal characters), while the reading frame +1 encodes a potential overlapping ORF (italics), which we named Sa</caption>
representative SARS or SARS-like coronaviruses. Our method allows comprehensive annotation of all ORFs. Surprisingly, the calculation of *X* motif enrichment is also accurate for the detection of overlapping genes, even though the codon usage and amino acid composition of overlapping genes is known to be significantly different from non-overlapping genes [15].

We showed that the predictions made by the GOFIX method have high sensitivity and specificity compared to the known functional ORFs in the well characterized SARS-CoV. For example, the annotated ORFs that have been described previously as non-functional or redundant, notably ORF3b and ORF8b, are not predicted to be functional by GOFIX. In contrast, we identified a putative small ORF overlapping the RBD of the Spike protein in SARS-CoV, that is conserved in Civet-CoV and Bat-CoV strain WIV16. Protein sequence analysis predicts that this novel ORF codes for a double-membrane spanning protein.

We then used the method GOFIX to compare all putative ORFs in representative genomes, and showed that most are conserved in all genomes, including the structural proteins (S, E, M and N) and accessory proteins 3a, 6, 7a, 7b, 9b and 9c. However, a number of ORFs were predicted to be non-functional, notably ORF8b in SARS-CoV and ORF3b in all genomes. We also identified potential new ORFs, including ORF9d in Pangolin-CoV and ORF10 in all genomes.

Concerning SARS-CoV-2, to date, the coding potential of SARS-CoV-2 remains partially unknown, and distinct studies have provided different genome annotations [37–39]. Overall, the genome of SARS-CoV-2 has 89% nucleotide identity with bat SARS-like-CoV (ZXC21) and 82% with that of human SARS-CoV [40]. A recent annotation [39] of the SARS-CoV and SARS-CoV-2 genomes identified 380 amino acid substitutions in 27 shared proteins, including the four structural proteins and eight accessory proteins, named 3a, 3b, p6, 7a, 7b, 8b, 9b and orf14 (corresponding to ORF9c here). Our analysis is in agreement with the previous studies showing that the genome organization is generally conserved. In particular, ORF9b and ORF9c are predicted to be expressed in SARS-CoV-2 genome. As expected, the structural proteins, S, E, M and N are conserved and have similar XME scores. ORF3a, ORF6 and ORF9b in SARS-CoV-2 also have similar XME scores to SARS-CoV.

Our ab initio analysis also allowed us to highlight some important specificities of the SARS-CoV-2 genome. Previously identified differences include some interferon antagonists and inflammasome activators encoded by SARS-CoV that are not conserved in SARS-CoV-2, in particular ORF8 in SARS-CoV-2 and ORF8a,b in SARS-CoV. Recent annotations of ORF3b are conflictual. For example, some authors have predicted that

SARS-CoV-2 ORF3b is homologous to SARS-CoV ORF3b [40], although the proposed SARS-CoV-2 protein is shorter with only 22 amino acids. In contrast, the SARS-CoV-2 ORF3b observed in [13] is not coded by the same region as SARS-CoV-2 ORF3b and the protein sequence is completely different. Here, we show that ORF3b has 0 *X* motifs in SARS-CoV-2, in agreement with the fact that little expression was observed in recent experiments aimed at characterizing the functions of SARS-CoV-2 proteins [13]. ORF10 is supposed to be unique to SARS-CoV-2, however it is also present in the Pangolin-CoV genome and its origin can be traced back to the Bat-CoV, where a truncated ORF of 26 amino acids, also present in the civet and human SARS-CoV genomes, can be found. Here, we observe that ORF7a, ORF7b and ORF9c have reduced XME scores in SARS-CoV-2. It remains to be seen whether these differences reflect functional divergences between SARS-CoV and SARS-CoV-2.

## Conclusions

In summary, we have developed a computational method GOFIX to characterize potential ORFs in virus genomes and applied the method to study the SARS-CoV-2 and related genomes. Our analysis of ORF coding potential helps to resolve some differences in current genome annotations. In addition, we suggest that some annotated ORFs may not be functional and predict novel putative ORFs in some genomes. Our findings contribute to characterizing sequence properties of accessory genes of SARS coronaviruses, and especially the newly acquired genes making use of overlapping reading frames.

## Abbreviations

*CoV*: Coronavirus; *GOFIX*: Gene prediction by Open reading Frame Identification using *X* motifs; *MERS*: Middle East Respiratory Syndrome; *nt*: Nucleotide; *ORF*: Open Reading Frame; *SARS*: Severe Acute Respiratory Syndrome; *XME*: *X* Motif Enrichment

## Acknowledgements

The authors would like to thank the BiGeSt Bioinformatics Platform for assistance.

## Authors' contributions

CJM, OP and JDT participated in the conceptualization of the work. CJM and JDT developed the methods and analyzed the data. CM and OP analyzed and interpreted the data. All authors contributed to writing the manuscript. All authors read and approved the final manuscript.

## Funding

This work was supported by the Institut Français de Bioinformatique (IFB), ANR-11-INBS-0013, the ANR projects Elixir-Excelerate: GA-676559 and RAin-RARE: ANR-18-RAR3-0006-02, and institute funds from the French Centre National de la Recherche Scientifique and the University of Strasbourg.

## Availability of data and materials

The datasets analysed during the current study are available in the Genbank viral genomes database, https://www.ncbi.nlm.nih.gov/genome/viruses/

## Ethics approval and consent to participate

Not applicable.
**Consent for publication**

Not applicable.

**Competing interests**

The authors declare that they have no competing interests.

**Author details**

¹Laboratoire ICube, Department of Computer Science, CNRS, University of Strasbourg, F-67412 Strasbourg, France. ²Unité de Microbiologie Structurale, Institut Pasteur, CNRS UMR 3528, 75724 Paris Cedex 15, France. ³Université Paris Diderot, Sorbonne Paris Cité, 75724 Paris Cedex 15, France.

Received: 3 June 2020 Accepted: 16 August 2020

Published online: 27 August 2020

**References**

1. Cui J, Li F, Shi Z. Origin and evolution of pathogenic coronaviruses. Nat Rev Microbiol. 2019;17:181-92.

2. Coronaviridae Study Group of the International Committee on Taxonomy of Viruses. The species severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol. 2020;5:536-44.

3. Ashour HM, Elkhatib WF, Rahman MM, Elshabrawy HA. Insights into the Recent 2019 Novel Coronavirus (SARS-CoV-2) in Light of Past Human Coronavirus Outbreaks. Pathogens. 2020;9:E186.

4. Schaecher SR, Pekosz A. SARS coronavirus accessory gene expression and function. Mol Biol SARS-Coronavirus. 2009;22:153-66.

5. Liu DX, Fung TS, Chong KK, Shukla A, Hilgenfeld R. Accessory proteins of SARS-CoV and other coronaviruses. Antivir Res. 2014;109:97-109.

6. Cagliani R, Forni D, Clerici M, Sironi M. Computational inference of selection underlying the evolution of the novel coronavirus, SARS-CoV-2. J Virol. 2020; 94:e00411-00420.

7. Khailany RA, Safdar M, Ozaslan M. Genomic characterization of a novel SARS-CoV-2. Gene Rep. 2020;19:100682.

8. Wang C, Liu Z, Chen Z, Huang X, Xu M, He T, Zhang Z. The establishment of reference sequence for SARS-CoV-2 and variation analysis. J Med Virol. 2020;92: 667-74.

9. Jin Z, Du X, Xu Y, et al. Structure of MPRO from COVID-19 virus and discovery of its inhibitors. Nature. 2020. https://doi.org/10.1038/s41586-020-2223-y.

10. Hussain M, Jabeen N, Raza F, Shabbir S, Baig AA, Amanullah A, Aziz B. Structural variations in human ACE2 may influence its binding with SARS-CoV-2 spike protein. J Med Virol. 2020. https://doi.org/10.1002/jmv.25832.

11. Srinivasan S, Cui H, Gao Z, Liu M, Lu S, Mkandawire W, Narykov O, Sun M, Korkin D. Structural genomics of SARS-CoV-2 indicates evolutionary conserved functional regions of viral proteins. Viruses. 2020;12:360.

12. Kim D, Lee JY, Yang JS, Kim JW, Kim VN, Chang H. The Architecture of SARS-CoV-2 Transcriptome. Cell. 2020;181:914-921.e10.

13. Gordon DE, et al. A SARS-CoV-2-human protein-protein interaction map reveals drug targets and potential drug repurposing. Nature. 2020. https://doi.org/10.1038/s41586-020-2286-9.

14. Yuen KS, Ye ZW, Fung SY, Chan CP, Jin DY. SARS-CoV-2 and COVID-19: the most important research questions. Cell Biosci. 2020;10:40.

15. Pavesi A, Vianelli A, Chirico N, Bao Y, Blinkova O, Belshaw R, Firth A, Karlin D. Overlapping genes and the proteins they encode differ significantly in their sequence composition from non-overlapping genes. PLoS One. 2018;13: e0202513.

16. Zhang KY, Gao YZ, Du MZ, Liu S, Dong C, Guo FB. Vgas: a viral genome annotation system. Front Microbiol. 2019;10:184.

17. Firth AE. Mapping overlapping functional elements embedded within the protein-coding regions of RNA viruses. Nucleic Acids Res. 2014;42:2420-39.

18. Schlub TE, Buchmann JP, Holmes EC. A simple method to detect candidate overlapping genes in viruses using single genome sequences. Mol Biol Evol. 2018;35:2572-81.

19. Arqués DG, Michel CJ. A complementary circular code in the protein coding genes. J Theor Biol. 1996;182:45-58.

20. Michel CJ. The Maximal C³ Self-Complementary Trinucleotide Circular Code X in Genes of Bacteria, Archaea, Eukaryotes, Plasmids and Viruses. Life (Basel). 2017;7:E20.

21. Dila G, Ripp R, Mayer C, Poch O, Michel CJ, Thompson JD. Circular code motifs in the ribosome: a missing link in the evolution of translation? RNA. 2019;25:1714-30.

22. Michel CJ. Circular code motifs in transfer and 16S ribosomal RNAs: a possible translation code in genes. Comput Biol Chem. 2012;37:24-37.

23. El Soufi K, Michel CJ. Unitary circular code motifs in genomes of eukaryotes. Biosystems. 2017;153:45-62.

24. El Soufi K, Michel CJ. Circular code motifs in genomes of eukaryotes. J Theor Biol. 2016;408:198-212.

25. Zhang T, Wu Q, Zhang Z. Probable pangolin origin of SARS-CoV-2 associated with the COVID-19 outbreak. Curr Biol. 2020;30:1578.

26. Kopecky-Bromberg SA, Martínez-Sobrido L, Frieman M, Baric RA, Palese P. Severe acute respiratory syndrome coronavirus open reading frame (ORF) 3b, ORF 6, and nucleocapsid proteins function as interferon antagonists. J Virol. 2007;81:548-57.

27. McBride R, Fielding BC. The role of severe acute respiratory syndrome (SARS)-coronavirus accessory proteins in virus pathogenesis. Viruses. 2012;4: 2902-23.

28. Yount B, Roberts RS, Sims AC, Deming D, Frieman MB, Sparks J, Denison MR, Davis N, Baric RS. Severe acute respiratory syndrome coronavirus group-specific open reading frames encode nonessential functions for replication in cell cultures and mice. J Virol. 2005;79:14909-22.

29. Ceraolo C, Giorgi FM. Genomic variance of the 2019-nCoV coronavirus. J Med Virol. 2020;92:522-8.

30. Xu K, Zheng B-J, Zeng R, Lu W, Lin Y-P. Severe acute respiratory syndrome coronavirus accessory protein 9b is a virion-associated protein. Virology. 2009;388:279-85.

31. Oostra M, de Haan CA, Rottier PJ. The 29-nucleotide deletion present in human but not in animal severe acute respiratory syndrome coronavirus disrupts the functional expression of open reading frame 8. J Virol. 2007;81: 13876-88.

32. Chen C-Y, Ping Y-H, Lee H-C, Chen K-H, Lee Y-M. Open reading frame 8a of the human severe acute respiratory syndrome coronavirus not only promotes viral replication but also induces apoptosis. J Infect Dis. 2007;196: 405-15.

33. Su YCF, Anderson DE, Young BE, Zhu F, Linster M, Kalimuddin S, Low JGH, Yan Z, Jayakumar J, Sun L, Yan GZ, Mendenhall IH, Leo Y-S, Lye DC, Wang L-F, Smith GJD. Discovery of a 382-nt deletion during the early evolution of SARS-CoV-2. bioRxiv. 2020;03(11):987222. https://doi.org/10.1101/2020.03.11. 987222.

34. Shukla A, Hilgenfeld R. Acquisition of new protein domains by coronaviruses: analysis of overlapping genes coding for proteins N and 9b in SARS coronavirus. Virus Genes. 2015;50:29-38.

35. Wu F, Zhao S, Yu B, et al. A new coronavirus associated with human respiratory disease in China. Nature. 2020;579:265-9.

36. Yang XL, Hu B, Wang B, Wang MN, Zhang Q, Zhang W, Wu LJ, Ge XY, Zhang YZ, Daszak P, Wang LF, Shi ZL. Isolation and characterization of a novel bat coronavirus closely related to the direct progenitor of severe acute respiratory syndrome coronavirus. J Virol. 2015;90:3253-6.

37. Zhou P, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579:270-3.

38. Chan JF, Kok KH, Zhu Z, Chu H, To KK, Yuan S, Yuen KY. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg Microbes Infect. 2020;9:221-36.

39. Wu A, Peng Y, Huang B, et al. Genome composition and divergence of the novel coronavirus (2019-nCoV) originating in China. Cell Host Microbe. 2020; 27:325-8.

40. Xu J, Zhao S, Teng T, Abdalla AE, Zhu W, Xie L, Wang Y, Guo X. Systematic comparison of two animal-to-human transmitted human coronaviruses: SARS-CoV-2 and SARS-CoV. Viruses. 2020;12:244.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.