key: cord-0816429-yqnaphm0
authors: Neches, Russell Y.; Kyrpides, Nikos C.; Ouzounis, Christos A.
title: Atypical Divergence of SARS-CoV-2 Orf8 from Orf7a within the Coronavirus Lineage Suggests Potential Stealthy Viral Strategies in Immune Evasion
date: 2021-01-19
journal: mBio
DOI: 10.1128/mbio.03014-20
sha: da7c6900a99a5644b511423895139306c4c186a5
doc_id: 816429
cord_uid: yqnaphm0

Orf8, one of the most puzzling genes in the SARS lineage of coronaviruses, marks a unique and striking difference in genome organization between SARS-CoV-2 and SARS-CoV-1. Here, using sequence comparisons, we unequivocally reveal the distant sequence similarities between SARS-CoV-2 Orf8 with its SARS-CoV-1 counterparts and the X4-like genes of coronaviruses, including its highly divergent “paralog” gene Orf7a, whose product is a potential immune antagonist of known structure. Supervised sequence space walks unravel identity levels that drop below 10% and yet exhibit subtle conservation patterns in this novel superfamily, characterized by an immunoglobulin-like beta sandwich topology. We document the high accuracy of the sequence space walk process in detail and characterize the subgroups of the superfamily in sequence space by systematic annotation of gene and taxon groups. While SARS-CoV-1 Orf7a and Orf8 genes are most similar to bat virus sequences, their SARS-CoV-2 counterparts are closer to pangolin virus homologs, reflecting the fine structure of conservation patterns within the SARS-CoV-2 genomes. The divergence between Orf7a and Orf8 is exceptionally idiosyncratic, since Orf7a is more constrained, whereas Orf8 is subject to rampant change, a peculiar feature that may be related to hitherto-unknown viral infection strategies. Despite their common origin, the Orf7a and Orf8 protein families exhibit different modes of evolutionary trajectories within the coronavirus lineage, which might be partly attributable to their complex interactions with the mammalian host cell, reflected by a multitude of functional associations of Orf8 in SARS-CoV-2 compared to a very small number of interactions discovered for Orf7a.

identity to the human/civet strain (35) . Further confirming the origin of Orf8 from bat species, it has also been argued that this region may play a role in tracing the origin of SARS strains in epidemic outbreaks, since the reported deletions do not affect the survival of the SARS-CoV-1 virus (20) . Importantly, SARS-CoV-2 has an intact Orf8 (non-split, Orf8ab) gene, which is known to be absent from the more lethal MERS strain (20) .

Motivation. The Orf7a protein has a known Ig-like structure and exhibits large genome-level variation and yet relative conservation as a single protein family, whereas Orf8, which is predicted to have a similar Ig-like structure, exhibits fewer genome-level variations apart from the split a/b pair or a full deletion and is quite variable at the protein sequence level. It should be noted that neither specific disordered regions have been detected for either of the two genes (36) nor any peculiar codon usage patterns (37) . However, this peculiar trajectory of Orf8/Orf7a in SARS-CoV-2 indicates a stealth viral strategy that might be key to control virulence and pathogenicity, with roles in host cell-virus interactions. Here, we investigate the origins and evolution of Orf8 across the coronavirus lineage by sequence and recombination analysis, ultimately merging the two Ig-like families by sensitive searches and identifying the shared conserved residues that define the common Ig fold.

Patterns of the multiple sequence alignment. The multiple alignment, as presented, reveals a turbulent evolutionary history across multiple coronavirus strains for this pair of SARS-CoV-2 genes and their homologs ( Fig. 1) . At the top of the alignment, there are nine members of the SARS Orf8a and Orf8b lineage, not always corresponding to a cognate pair, and the three truncated structures (block i, Fig. 1A ). The next set (block ii) is composed of 93 Orf7a homologs, excluding those 3 of known structure (2 for SARS-CoV-1 and 1 for SARS-CoV-2, intertwined within the first block, due to their truncated N termini, block i thus containing 12 members). The Orf7a group is followed by a quite heterogeneous set of 24 X4-like/Orf9 sequences from bats, Orf8 from pangolin and Orf9 from SARS-CoV-1 (block iii): despite different names and significant variation, they represent similar sequences of common origin. Finally, the bottom group (block iv) comprises 67 Orf8 sequences from SARS-CoV-2 intermixed with those from bat or pangolin hosts. The mutations V71L (V62L) and L96S (L84S) reported elsewhere (38) are clearly seen within the Orf8 family, among others (Fig. 1A , see DS04 at https:// doi.org/10.6084/m9.figshare.12678491.v1). Of 265 processed SARS-CoV-2 genomes, only two samples from Wuhan, China (WH19002/2019 and WH19004/2020), along with the USAWA-UW413/2020 case, exhibit no mutations compared to the original reference strain (Wuhan-Hu-1) (see DS.snp).

Sequence conservation across virus strains. The first 15 residues are considered to function as signal peptides for Orf7a and Orf8 proteins and were long thought to be the only element shared between them (13) . We now extend the signal peptide similarity throughout the superfamily and suggest a common origin for those genes, augmenting independently obtained results by parallel efforts (23) . The conservation pattern extends across the entire length of the superfamily, based on evidence from the database searches and the multiple alignment (Fig. 1A) . It is notable that the N-terminal region is aligned separately from the downstream protein sequence, suggesting a potential cleavage site. Interestingly, most residues at position 2 are lysines, with two exceptions shown for SARS-CoV-2: one substitution by glutamate (QJF76096) and another by arginine (QIU91254), similar to bat X4-like/Orf9 entries (in the middle of the alignment, block iii). Then, from position 15 onward, the proteins of known structure suggest the presence of the Ig-like beta-sandwich, and an insertion of ;20 residues in the bottom part of the group (block iv, X4-like and Orf8), mostly covered by the regular expression {W.[I,L,V][K,R][I,V,Y]}. The Orf8a/b pair of SARS-CoV-1 underlines the structural flexibility of this segment and the possibility of having two distinct protein products that associate to perform the role of Orf8 in this viral strain. A shorter insert of about 10 residues is also found in Orf8b and some of its Orf9 bat homologs, mostly covered by the regular expression {L[I,V].RC}. Note that these two segments are (52) . (b) Sequence-structure conservation for SARS-CoV-1 Orf7a (1xak), 10 conserved residues (apart from three residues found in the signal sequence but not available in the structure) are indicated by red; disulfide bridges are indicated by green, rendered by UCSF Chimera (59) . (c) Cladogram based on the final alignment, with SARS-CoV-2 sequences marked by red color, generated by FastTree (55) (DS.tree), visualized by IcyTree (57) . The block iv group containing Orf8 homologs exhibits the widest variation. Due to significant variation, a large number of gaps or a fragmented gene structure, a cladogram is selected instead of a more typical phylogenetic tree to depict the structure of the superfamily. matched onto members of the Orf8ab pair, underlining the likely origin of the intact SARS-CoV-2 Orf8 from those, as mentioned above.

Conserved positions and structural interpretation. We define 13 highly conserved positions in the full (final) alignment: M1, K2, and L7 (all in the predicted signal peptide) and H22, Q28, C30, P41, C42, Y47, P79, C95, C107, and V135 (residue coordinates follow Fig. 1A , see DS04 at https://doi.org/10.6084/m9.figshare.12678491.v1). It is known, for Orf7a, that two disulfide bonds are present in positions C30-C90 and C42-C107 (17) . Only a few sequences show mutations in positions C30F, C42S, and C107F. The 10 conserved positions are mapped onto the three-dimensional structure of SARS-CoV-1 Orf7a (1XAK) (17) , thus defining the invariant amino acid residues for the Ig-like fold of the SARS lineage (Fig. 1B) , presumed critical for structural and perhaps functional integrity for the entire protein superfamily. The dendrogram derived from the final alignment reveals the groupings for the superfamily and the distribution of SARS-CoV-2 sequences (DS.tree; Fig. 1C ).

Block annotation and consistency checking. An annotated version of the alignment was generated, manually trimming low-occupancy positions, and marking the block structure (108 residues long, see DS05 at https://doi.org/10.6084/m9.figshare .12678491.v1 see Fig. S2A ). From the annotated alignment (DS05), an automatically trimmed alignment (39) was also generated (54 residues long, see DS06 at https://doi .org/10.6084/m9.figshare.12678491.v1; see Fig. S2B ), with two of the conserved positions not included (Y47, P79), in addition to removing most of the variation from rapidly changing residues across family members. This alignment reveals a consistent picture further suggesting that the common origin of these genes is indeed determined by the positions that define the Orf7a Ig-like fold, with Orf8 being subjected to rapid evolutionary change.

Phylogenetic tree and sequence space. An orthogonal approach to tree inference is a multidimensional alignment embedding in sequence space (40) , which confirms the unique evolutionary history for the superfamily, not readily seen in a tree graph (Fig. 2) . The distances in the three-dimensional (3D) embedding represent variation seen across multiple groups: in particular, within-family distances from Orf7a groups are clearly less variable compared to the Orf8 groups, supporting the notion of Orf8 as a highly variable protein in coronaviruses. However, the distances between SARS-CoV-1 and SARS-CoV-2 strain groups are comparable (Fig. 2) . Predictably, when nonconserved residues are ignored, the within-family distances appear similar. SARS-CoV-1 Orf7a and Orf8 (also the concatenated ab) come closer than their SARS-CoV-2 counterparts, indicating that there is a faster evolution in the latter, also affected by more extensive sampling of the sequence space. This pattern is also reflected in a tree comparison tanglegram, with the major groups remaining stable, as the fast-evolving nodes alter the tree topology (see Fig. S3 ).

Phylogenetic profiling across the coronavirus pangenome. For both the Orf7a and Orf8 SARS-related families, the phylogenetic distribution of their members is distinctly restricted, supported by an ongoing coronavirus pangenome study (8) , which lists these genes separately. We have confirmed a single key exception (X4-like, block iii) in the alpha group (QCX35187.1, https://www.ncbi.nlm.nih.gov/nuccore/MK720946 .1/) and its X4-like homologs (e.g., QCX35176.1), as previously shown (23) . When mapped onto the tree of 89 representative coronavirus strains derived from DNAbased whole-genome alignment (41) , the restricted distribution pattern is clearly visible (see Fig. S4A ). In addition, no clear evidence for recombination for this region of the SARS-CoV-2 genome is detectable across the 89 strains, markedly challenging previous claims (35, 42) and further supporting the scenario of gene duplication and extreme divergence. No homologs of this superfamily have ever been detected outside the coronavirus group. Furthermore, we were not able to confirm any other similarities at the structural level (23); these remain particularly valuable predictions for future research.

From structure to function, as well as constraints from protein interactions. The discovered interactions of these proteins for SARS-CoV-2 with the mammalian host proteins reveal a contrasting picture (43) : while Orf7a is relatively more conserved, only two interactions have been detected for this protein (https://amp.pharm.mssm.edu/ covid19/genesets/20). On the other hand, 47 interactions have been detected for Orf8, which has a much faster evolutionary rate (https://amp.pharm.mssm.edu/covid19/ genesets/22) (43) . This could be an artifact related to the coverage of protein interactions between the protein complement of the virus and the host cell, and yet it may also point to a more complex interaction landscape for Orf8, underlining its importance in viral strategies to invade cells and propagate.

Here, we show for the first time that remote, nontrivial sequence similarities between the SARS-CoV-2 proteins Orf7a and Orf8 are detectable using supervised sequence space walks in database searches, aimed at precision and reproducibility (44) . The detection of similarity between Orf8 and Orf7a, a protein of known structure with an Ig-like fold, implies the importance of this gene duplication for virus biology and confirms previous results, independently derived by different approaches, i.e., profile-profile searches and protein modeling (23) . We assessed the extent at which sequence comparisons alone can establish unambiguously the homology between Orf8 and Orf7a family members within the coronavirus lineage and entire pangenome. The embedding of sequence conservation patterns in a multidimensional space reveals atypical divergence patterns, with "equidistant" groupings for Orf7a and Orf8 across SARS-CoV-1 and SARS-CoV-2 but significant divergence of Orf8 compared to within-family distances for Orf7a (Fig. 2) . Moreover, the conservation of Orf7a compared to Orf8 is clearly demonstrated not only by its similarity to bat viruses for SARS-CoV-1 and pangolin-Manis javanica viruses for SARS-CoV-2, respectively, but with the stark difference between Orf8 in these strains (;35% identity for SARS-CoV-1 Chinese bat-Rhinolophus sinicus virus homologs compared to ;88% identity for SARS-CoV-2 pangolin virus homologs), a puzzling pattern (Fig. 3) .

During the final stages of this study, the structure of Orf8 from SARS-CoV-2 has been announced in a preprint (PDB 7JTL), confirming the presence of an Ig-like fold (45) and further supporting the notion of a common origin between the Orf7a and Orf8 families. Of the 13 conserved alignment positions, 9 are available (4 are missing at the N terminus): Q28, C30 (strand 1), P41, C42 (loop between strands 2 and 3), Y47 (strand 3), P79 (loop before last three strands in both structures, at 16.141 Å distance of C-alpha atoms, due to the presence of an Orf8-specific region [45] , thus excluded from superposition measurements), and C95, C107, and V135 (the latter also excluded due to uncertainty) are shifted and interpreted differently between the automatically derived profile-driven alignment here and the structure-based alignment (45) . Remarkably, when the two cysteines are shifted, the root mean square deviation between 7 C-alpha atom pairs drops from 9.883 Å (from their sequence-based match) to 3.398 Å (to their structure-based match), thus confirming the rather different outcomes of the sequence-and structure-based alignments.

We provide strong evidence for the peculiar divergence of Orf8 from Orf7a, within an otherwise dense viral genome and delimit the phylogenetic distribution of these two genes across coronaviruses. Neither of these genes is found in the gamma or delta coronavirus groups, suggestive of a likely loss in those coronavirus sublineages (see Fig. S4BC ). It is quite perplexing that no member of this family is present in the MERS clade (23) . Given that Orf7a with an Ig-like structure is potentially an immune antagonist with a pivotal role in the viral infection strategy and the recent observation that Orf8 downregulates MHC-I (46), the Orf7a/Orf8 superfamily might be a key system for immune evasion, known for other analogous cases, including herpesviruses, poxviruses, and adenoviruses (47) . The detected protein interactions of SARS-CoV-2 exhibit such a sharp contrast between Orf7a and Orf8 (43) that, barring a technical artifact with respect to coverage, the hypothesis arises for Orf7a being used as a conserved template, to generate variants, such as Orf8, wreaking havoc through immune evasion in the host cell.

In summary, the pair Orf7a/Orf8 may be an under appreciated element in SARS-CoV-2 biology, given its peculiar and quite unusual patterns of divergence and functional properties that might be related to virulence and the pathogenicity of this strain which has caused the COVID-19 pandemic.

Database searches and documentation. Using the Orf8 protein sequence of SARS-CoV-2 as a primary query, we extensively searched NCBI's 'nr' protein collection (48), following an approach described elsewhere (44) , including masking by CAST (49) and permissive E value thresholds in supervised mode. This approach, named "sequence space walk," allows the use of permissive thresholds to detect weak sequence similarities, by inspecting all hits at each single iteration. We specifically searched the Virus subset of 'nr' using PSI-BLAST (50) and the following parameters: max target sequences, 500; Expect 10; word-size 2; BLOSUM62; gap costs, 11/1; compositional adjustment, none; Filter, none; and E value, 0.1. Against the full 'nr', the search exits early at ;200 hits with the same parameters (not shown). Fragments and sequences with unspecified residues were excluded (see below). Search statistics are provided (Table 1) , along with the hit tables and PSSMs (the PSSM is unavailable for step 1 only) (DS.runs).

Other software tools. Sequence alignments were performed with MAFFT and default parameters (51) and visualized with JalView 2 (52) . Family analysis was performed with AlignmentViewer (40) , and dimensionality reduction for sequence space exploration (2D and 3D) was aided by UMAP (53) , as implemented in AlignmentViewer. The color coding for sequence glyphs was selected according to mView (54) using a residue width of 2 pixels and a residue height of 1 pixel. Validation of sequence database searches, identifier tracking, and alignment quality checks were performed using various scripts and alignments were automatically trimmed with TrimAl (39), unless stated otherwise. Trees were inferred using FastTree (55) on the NGPhylogeny.fr servers with the LG/gamma options (56) and visualized with IcyTree (57). Tree annotation for FIG 3 Closest homologs of the SARS-CoV-1 and SARS-CoV-2 for Orf7a and Orf8. A restricted search against the virus subset of 'nr' that excludes SARS-1/-2 (taxid 694009 and 2697049, respectively) as targets reveals the similarity of superfamily members to coronavirus strains from other hosts. Orf7a of SARS-1 is most similar to Chinese bat (Rhinolophus sinicus) strains (;95% identity). The concatenated Orf8ab (from Orf8a/Orf8b) shows the closest-yet lower-similarity to Chinese bats (;35% identity), since the host species Rhinolophus ferrumequinum is excluded. In contrast, Orf7a of SARS-2 is most similar to pangolin (Manis javanica) strains (;97% identity), with Orf8 exhibiting high levels of similarity again to pangolin (;88% identity). Identifiers next to the host name are provided. The scale (sequence identity) is shown on the right. "NA" signifies that there are no similarities detected in this restricted sequence search. The diagram was generated by using ClustVis (66) . SARS-CoV-1 is referred to here as "SARS-1," and SARS-CoV-2 is referred to here as "SARS-2." Orf8 for "SARS-1" corresponds to a concatenated Orf8ab.

genes and clades was supported by MicroReact (58) . Protein structure analysis and annotation was performed with UCSF Chimera (59) . Protein interaction data were obtained from the Covid-19 drug/gene set library (60) . Variations were calculated using Virulign (61) over a representative sample of SARS-CoV-2 genomes. Recombination breakpoints were identified using BALi-Phy (62) to compute the pairwise homoplasy index across the genome alignment (63) . Regions without breakpoints were identified as nonrecombining regions (15) . Further details are given in Text S1 in the supplemental material.

Supervised sequence space walk. The initial query sequence Orf8 from SARS-CoV-2 detects its closest homologs within the Orf8 family (PF12093) (29) , followed at step 3 by distant homologs of the Orf7a family (PF08779) (17, 64) , variably called X4 or Orf10, from a range of viral hosts that include bats, pangolins and civets. The search converges at step 9 with high precision, i.e., no known false positives (Fig. 4a) . The detected region of homology spans the length of these proteins beyond the signal sequence, as observed independently (23) , and unifies the two Pfam entries into a single superfamily (65), with few invariant residues across its members (see below). All PSSMs are stored for future searches (see Text S1, DS.runs) and re-use by the community. Since the results are time sensitive at present (July 2020), the PSSM collection is provided to ensure reproducibility: very many and highly similar sequences will be deposited in the meanwhile, without substantially modifying the main conclusions or challenging the validity of the reported sequence search results.

Identifier processing. There are 24 sequence entries that are lost (as false negatives) in the iterative search, 20 of which reappear in subsequent iterations (Fig. 4b) . Of those, 1 entry is lost at step 2, 11 are lost at step 3, another 11 are lost at step 4, and just 1 is lost at step 6. All lost entries were eventually recovered, except for 4 (all of representing sequences shorter than 25 residues). A link with the 24 entries is provided for further inspection (identifiers are available in DS.runs).

Quality control and alignment. The iterative search returns 465 unique entries, many of which are partial sequences (of short length) or have multiple undefined residues (see DS01 at https://doi.org/10 .6084/m9.figshare.12678491.v1); just 1 entry (QJI07349.1) returns two hits in separate regions, due to 70 undefined middle positions (one duplicate ID). When partial/undefined sequences were removed, 288 entries remained, which were then aligned by MAFFT, revealing the conservation patterns of the superfamily (Fig. 4c, (31) .

Multiple alignment editing and statistics. To define a more accurate, reference alignment, all retained (n = 196) sequences were set to be $80 residues long (#100 for Orf7a and .100 for Orf8, plus other members of the superfamily to demonstrate variation), manually removing the C-terminal end for five low-occupancy positions. Exceptions were the Orf8a/b pair and three (truncated) structure entries from the PDB (1XAK 68/83 [17] , 1YO4 69/87 [64] , and 6W37 67/67). The reference alignment is provided (Fig. 4d, final alignment, 196 sequences, 141 residues long, see DS04 at https://doi.org/10.6084/m9.figshare.12678491.v1): at 80% identity threshold, we defined 13 conserved residues (see below), 9% (13/141) of the length of the aligned sequences.

Cross-family sequence similarities. The sequence space walk with Orf8 revealed its closest homologs at the first steps and at step 3 connected them with Orf7a and its homologs (Fig. 4a) , with minimum sequence identity dropping from 20 to 8% at convergence (Table 1) . This outcome is also confirmed by the resulting alignments, with an average sequence identity dropping from 40 to 10%. The alignment depiction as glyphs for both the initial (Fig. 4c ) and the final (Fig. 4d) alignments is also reflected in the matrix of pairwise sequence similarities for the former (see Fig. S1a ) and the latter (see Fig. S1b ), respectively; these cross-similarity patterns were generated by AlignmentViewer (40) . (Note that in this study, for data annotation purposes only, SARS-CoV-1 and SARS-CoV-2 are marked as "SARS-1" and "SARS-2," respectively.) 

Supplemental material is available online only. 

China Novel Coronavirus Investigating and Research Team. 2020. A novel coronavirus from patients with pneumonia in China

A new coronavirus associated with human respiratory disease in China

A pneumonia outbreak associated with a new coronavirus of probable bat origin

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

Origin and evolution of pathogenic coronaviruses

Bioinformatics and functional analyses of coronavirus nonstructural proteins involved in the formation of replicative organelles

Overview of lethal human coronaviruses

Functional pangenome analysis suggests inhibition of the protein E as a readily available therapy for COVID-2019

The architecture of SARS-CoV-2 transcriptome

Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus

Systematic comparison of two animal-to-human transmitted human coronaviruses: SARS-CoV-2 and SARS-CoV

On the origin and continuing evolution of SARS-CoV

The 29-nucleotide deletion present in human but not in animal severe acute respiratory syndrome coronaviruses disrupts the functional expression of open reading frame 8

Discovery of a 382-nt deletion during the early evolution of SARS-CoV-2

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

Are pangolins the intermediate host of the 2019 novel coronavirus (SARS-CoV-2)?

Structure and intracellular targeting of the SARS-coronavirus Orf7a accessory protein

Immunoglobulin domains in Escherichia coli and other enterobacteria: from pathogenesis to applications in antibody technologies

Cross-host evolution of severe acute respiratory syndrome coronavirus in palm civet and human

ORF8-related genetic evidence for Chinese horseshoe bats as the source of human severe acute respiratory syndrome coronavirus

Effects of a major deletion in the SARS-CoV-2 genome on the severity of infection and the inflammatory response: an observational cohort study

An 81-nucleotide deletion in SARS-CoV-2 ORF7a identified from sentinel surveillance in Arizona

Novel immunoglobulin domain proteins provide insights into evolution and pathogenesis of SARS-CoV-2-related viruses

Genomic surveillance of SARS-CoV-2 in Thailand reveals mixed imported populations, a local lineage expansion and a virus with truncated ORF7a

Identification of multiple large deletions in ORF7a resulting in inframe gene fusions in clinical SARS-CoV-2 isolates

The 8ab protein of SARS-CoV is a luminal ER membrane-associated protein and induces the activation of ATF6

SARS-coronavirus open reading frame-8b triggers intracellular stress pathways and activates NLRP3 inflammasomes

Accessory proteins 8b and 8ab of severe acute respiratory syndrome coronavirus suppress the interferon signaling pathway by mediating ubiquitin-dependent rapid degradation of interferon regulatory factor 3

The human severe acute respiratory syndrome coronavirus (SARS-CoV) 8b protein is distinct from its counterpart in animal SARS-CoV and down-regulates the expression of the envelope protein in infected cells

Molecular and biochemical characterization of the SARS-CoV accessory proteins ORF8a, ORF8b, and ORF8ab

Molecular biology of the SARS-coronavirus

The role of severe acute respiratory syndrome (SARS)-coronavirus accessory proteins in virus pathogenesis

Structural genomics and interactomics of 2019 Wuhan novel coronavirus, SARS-CoV-2, indicate evolutionary conserved functional regions of viral proteins

Emergence of SARS-CoV-2 through recombination and strong purifying selection

Comparative genomics suggests limited variability and similar evolutionary patterns between major clades of SARS-Cov-2

Severe acute respiratory syndrome (SARS) coronavirus ORF8 protein is acquired from SARS-related coronavirus from greater horseshoe bats through recombination

Dark proteome of newly emerged SARS-CoV-2 in comparison with human and bat coronaviruses

Sequence analysis of SARS-CoV-2 genome reveals features important for vaccine design

Virus against virus: a potential treatment for 2019-nCov (SARS-CoV-2) and other RNA viruses

trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses

AlignmentViewer: sequence analysis of large protein families (version 1)

A path towards SARS-CoV-2 attenuation: metabolic pressure on CTP synthesis rules the virus evolution

SARS-CoV-2 ORF8 and SARS-CoV ORF8ab: genomic divergence and functional convergence

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing

Sequence evidence for common ancestry of eukaryotic endomembrane coatomers

Structure of SARS-CoV-2 ORF8, a rapidly evolving coronavirus protein implicated in immune evasion

The ORF8 protein of SARS-CoV-2 mediates immune evasion through potently downregulating MHC-I

Immunoglobulin superfamily members encoded by viruses and their multiple roles in immune evasion

Database resources of the National Center for Biotechnology Information

CAST: an iterative algorithm for the complexity analysis of sequence tracts

Gapped BLAST and PSI-BLAST: a new generation of protein database search programs

MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization

Jalview Version 2-a multiple sequence alignment editor and analysis workbench

Dimensionality reduction for visualizing single-cell data using UMAP

MView: a web-compatible database search or multiple alignment viewer

FastTree 2-approximately maximumlikelihood trees for large alignments

NGPhylogeny.fr: new generation phylogenetic services for non-specialists

IcyTree: rapid browser-based visualization for phylogenetic trees and networks

Microreact: visualizing and sharing data for genomic epidemiology and phylogeography. Microb Genom 2:e000093

UCSF Chimera-a visualization system for exploratory research and analysis

The COVID-19 gene and drug set library

VIRULIGN: fast codoncorrect alignment and annotation of viral genomes

BAli-Phy: simultaneous Bayesian inference of alignment and phylogeny

Genome-wide inference of ancestral recombination graphs

Solution structure of the X4 protein coded by the SARS-related coronavirus reveals an immunoglobulin like fold and suggests a binding activity to integrin I domains

The Pfam protein families database in 2019

ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap

dendextend: an R package for visualizing, adjusting, and comparing trees of hierarchical clustering

We thank Konstantinos Kyritsis (School of Pharmacy, Aristotle University of Thessalonica) for help with the "dendextend" tree comparison. C.A.O. acknowledges support by the project Elixir-GR, implemented under the Action "Reinforcement of the Research and Innovation Infrastructure," funded by the Operational Program Competitiveness, Entrepreneurship, and Innovation (NSRF 2014-2020) and cofinanced by Greece and the European Union (European Regional