SARS-CoV-2 genomes from Oklahoma, USA
=====================================

* Sai Narayanan
* John C Ritchey
* Girish Patil
* Teluguakula Narasaraju
* Sunil More
* Jerry Malayer
* Jeremiah Saliki
* Anil Kaul
* Akhilesh Ramachandran

## Abstract

Genomic sequencing has played a major role in understanding the pathogenicity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). With the current pandemic, it is essential that SARS-CoV-2 viruses are sequenced regularly to determine mutations and genomic modifications in different geographical locations. In this study we sequenced SARS-CoV-2 from 5 clinical samples obtained in Oklahoma, USA during different time points of pandemic presence in the state. One sample from the initial days of the pandemic in the state and 4 during the peak in Oklahoma were sequenced. Previously reported mutations including D614G in S gene, P4715L in ORF1ab, S194L, R203K and G204R in N gene were identified in the genomes sequenced in this study. Possible novel mutations were also detected such as G1167V in S gene, A6269S and P3371S in ORF1ab, T28I in ORF7b, G96R in ORF8. Phylogenetic analysis of the genomes showed similarity to viruses from across the globe. These novel mutations and phylogenetic analysis emphasize the contagious nature of the virus.

## Introduction

Towards the end of 2019, several individuals with signs of pneumonia reported to hospitals in Wuhan, the capital of Hubei province in Central China. The etiological agent was identified to be a novel coronavirus (SARS-CoV-2/ nCoV-19) on 7th January 2020(1). Human to human transmission was recorded around the same time(2). By the end of January 2020, WHO declared a ‘public health emergency of international concern’(3). As of September 5th, 2020, a total of 26,696,160 confirmed patients and 875,943 deaths have been reported worldwide and 63,187 registered cases and over 850 deaths in the state of Oklahoma, USA(4).

Numerous coronaviruses infecting different animal species including humans have been identified. SARS-CoV-2 is an enveloped positive sense single stranded RNA virus belonging to the genus *Betacoronaviru*s (subgenus *Sarbecovirus*) in the *Coronaviridae* family (order: *Nidovirales*). Based on genomic sequence analysis, it is reported to have originated from bats(5) and pangolins(6, 7). Other previously identified coronaviruses known to infect humans include SARS-CoV-1, MERS-CoV, HCoV-NL63, HCoV-229E, HCoV-OC43 and HCoV-HKU1(8–13).

Genomic sequence data is important in identifying, characterizing and understanding pathogens (14, 15). It can shed light on pathogenicity, virulence, drug/vaccine targets, mutation sites etc. and can also be critical in source attribution and determining microbial provenance (16). The first genome sequence of SARS-CoV-2 was made available on January 10th (GenBank ID: [MN908947.3](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=MN908947.3&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom)) (17). Multiple genomic sequences of SARS-CoV-2 from all over the world have since been deposited in public data bases such as GenBank and GISAID (Global Initiative on Sharing All Influenza Data; [www.epicov.org](http://www.epicov.org)) (18). This has facilitated extensive genomic studies leading to the identification of several mutations in the genome that can influence infectivity and virulence of the virus (16, 19, 20).

The genome of SARS-CoV-2 is similar to many other pathogenic coronaviruses and has multiple genes that code for different proteins such as S gene (Surface glycoprotein), N gene (nucleocapsid phosphoproteins), M gene (membrane glycoprotein), E gene (envelope protein), and open reading frames, such as ORF1a, ORF1b ORF3a,ORF6, ORF7a, ORF7b, ORF8 and ORF10 (16, 21). The leading sequence of the viral genome is the sequence for ORF1ab, which encodes for multiple proteins including replicase polyproteins, non-structural proteins, papain-like proteinase, RNA dependent RNA polymerase etc., which are essential for replication and survival in the host (21). Though the exact function of ORF3a is yet to be clearly understood, it is believed to play a major role in viral release after replication in SARS-CoV-1 (22). ORF7 and ORF8 code for accessory proteins, the functions of which are yet to be clearly understood (23, 24). Minimal roles of ORF7 and 8 in viral replication have been reported in SARS-CoV-1 alongside apoptosis stimulation of host cells(25).

The huge repository of sequence data in open-source databases such as the GenBank-NCBI and GISAID has facilitated the identification of numerous mutations and single nucleotide polymorphisms (SNPs) in the SARS-CoV-2 genome. SNPs in the genome that result in commonly reported non-synonymous mutations such as P4715L in ORF1ab, D614G in S gene, R203K and G204R in N gene are some of the commonly reported. P4715L in ORF1ab is believed to play a major role in interaction with other proteins that regulate RNA dependent RNA polymerase (26). D614G (Aspartate to Glycine) mutation in the S gene has been reported to result in increased transduction into human epithelial cells (20). N gene mutations R203K and G204R are believed to increase viral fitness, survival and adaptation to humans (27).

In this study, we sequenced the genome of SARS-CoV-2 from 5 human clinical samples received at the Oklahoma Animal Disease Diagnostic Laboratory (OADDL) at various times during the SARS-CoV-2 pandemic in Oklahoma. Multiple mutations were detected in the genomic sequences including those already reported as well as previously unreported mutations.

## Materials and methods

This study was approved by the Institutional Review Board (Application number: IRB-20-357) at Oklahoma State University, Stillwater OK 74078, USA.

### Clinical samples and processing

Five Nasopharyngeal swabs collected from human patients received at OADDL for COVID-19 testing between March 2020 – July 2020 were used in this study. Nucleic acid extraction was performed using MagMax Viral Pathogen Nucleic acid isolation kit (Thermofisher, MA, USA) as per manufacturer’s recommended protocols. Viral presence was detected by real-time PCR using TaqPath COVID-19 Multiplex Diagnostic Solutions (Thermofisher, MA, USA). All samples had a cycle threshold value between 18-22.

### Genomic sequencing

cDNA was synthesized from extracted RNA from five clinical samples that were positive for SARS-CoV-2. cDNA was then PCR amplified using ARTIC V3 primers ([https://github.com/artic-network/artic-ncov2019/tree/master/primer_schemes/nCoV-2019/V3](https://github.com/artic-network/artic-ncov2019/tree/master/primer_schemes/nCoV-2019/V3)) to obtain overlapping segments of the whole viral genome. DNA library repair (SQK-LSK-109, Oxford Nanopore Technologies, UK), Solid Phase Reversible Immobilization (SPRI) paramagnetic beads clean-up (AMPure XP, Beckmann Coulter, CA, USA) and adapter ligation and barcoding (NBD-001, Oxford Nanopore Technologies, UK) were done as per manufacturer recommendations. Libraries were then pooled and sequenced using MinION (Oxford Nanopore Technologies, UK) platform following manufacturer recommendations.

### Genome Assembly, Alignment and Phylogenetic analysis

Sequences obtained were assembled *de novo* using Canu (28). To further obtain a reliable consensus genome assembly, *de novo* assemblies and sequence output files were assembled to the SARS-CoV-2 reference genome Wuhan Hu-1 (GenBank ID: [MN908947.3](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=MN908947.3&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom)) with minimap2 (29) and Nanopolish (30).

To assess uniqueness of the genomes sequenced in this study, MAFFT (31, 32) was used to align whole sequences of the 5 genomes to SARS-CoV-2 reference genome (NC_045512.2).

Gene predictions on the consensus assemblies were made using Viral Genome ORF Reader 4 (33) (VIGOR4) using a curated library available in the Virus Pathogen Resource (ViPR) (34) database. Individual genes were aligned to SARS-CoV-2 Wuhan Hu-1 genome from NCBI (GenBank ID: [MN908947.3](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=MN908947.3&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom)) using MUSCLE aligner in MEGA-X (35) to identify SNPs and changes in the amino acid produced by the gene.

To assess similarity to previously reported genomes, a phylogenetic analysis was made using 9072 genomes of SARS-CoV-2 from the GenBank database and the 5 viruses sequenced in the study. A General Time Reversible (GTR) substitution model based Unweighted Pair Group Method with Arithmetic Mean (UPGMA) alignment was constructed using MAFFT and FastTree (36). Clade definitions for the sequences were identified using 9 marker variants reported for classification in the GISAID database (18).

## Results and Discussion

Five clinical rRT-PCR positive samples from different periods of the pandemic in Oklahoma, USA were chosen for the study. One of the samples (Oklahoma-ADDL-1) was received during the initial stages (April 2020) of the pandemic. An increased incidence rate of the disease was observed by the end of May 2020 (4). Three of the samples (Oklahoma-ADDL-2,3,4) sequenced were received during this period and the last sample (Oklahoma-ADDL-5) was received one month (June 2020) after this period. More than 4000X coverage was obtained at the end of sequencing for all the samples. Following *de novo* assembly with Canu, consensus genomes were obtained after reference genome assembly with nanopolish and minimap2.

The five viral genome assemblies were aligned to the reference genome (SARS-CoV-2 Wuhan-Hu-1 NC_045512.2) using MAFFT and genome similarities to the reference isolate were calculated (Table 1). Most of the genomes showed more than 99% similarity. For Oklahoma-ADDL-5, only a partial sequence could be generated and hence showed a lower 86.836% similarity. This could be due to reduced amplification during amplicon generation and also due to reference assembly. Other than Oklahoma-ADDL-5, Oklahoma-ADDL-1 sequenced from a clinical sample obtained during the beginning of the pandemic in Oklahoma (April 2020) showed a lower similarity to the reference genome when compared to the other four isolates.

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20195420/T1)

Table 1: 
Genome similarity of 5 sequenced genomes when compared to NCBI reference genome (NC_045512.2) along with total reads generated and genome coverage obtained

Genes were predicted using VIGOR4 and individual genes were aligned to their respective NCBI reference genes using MUSCLE aligner with UPGMA alignment in MEGA-X to assess mutations in the genome. Using MEGA-X visualization tool, various silent and missense mutations were detected. The missense/non-synonymous mutations detected in major genes are listed in Table2. While non-synonymous mutations were detected in ORF1ab, ORF1a, S, N, ORF3a, ORF7, ORF8, none were detected in Envelope (E) gene, Membrane glycoprotein (M) gene or ORF10 gene.

Non-synonymous mutation at amino acid location 81 (C→T; Arginine → Cysteine), which codes for nsp2 was present in all the genomes sequenced in the study. The potential implication of this mutation is unknown. nsp2 along with nsp3 are known to play major roles in pathogenesis (37) of SARS-CoV-2 in humans. A few other possibly novel and previously reported non-synonymous mutations in ORF1a and ORF1ab were also identified. A previously reported mutation in the ORF1ab gene P4715L(38) (Proline to Leucine) was recorded alongside novel mutations at various amino acid locations. P4715L mutation has been implicated to play a major role in interaction with other proteins that regulate RNA Dependent RNA polymerase activity. P3371S (Proline to Serine) mutation in ORF1ab and ORF1a was detected in Oklahoma-ADDL-5. Oklahoma-ADDL-4 carried mutation T4412A (Threonine to Alanine) in ORF1ab, while Oklahoma-ADDL-2 and 3 carried mutation A6269S (Alanine to Serine). ORF1ab has multiple functions including RNA dependent RNA polymerase activity, helicase activity, Fe-S cluster binding, Zn- binding activity, methyltransferase activity(39) etc. Functional implications of these mutations are still unknown and further studies to understand functional changes caused by these mutations may aid in better understanding the pathogenesis of the viral isolates found in Oklahoma.

Non-synonymous mutations were also found in S, ORF3a, ORF7b, ORF8 and N gene (Table 2). A previously reported mutation in S gene - D614G(40) (Aspartate to Glycine), was identified in all the genomes sequenced. The D614G mutation has been reported to cause a decrease in PCR cycle thresholds, suggestive of higher upper respiratory tract viral load (41, 42) in the host. A previously reported deleterious variation in the protein expressed from ORF3a, Q57H(40, 43), was recorded in Oklahoma-ADDL-1, the genome isolated at the beginning of the pandemic in Oklahoma. This mutation was not recorded in genomes isolated at other times in the state. The N gene also carried other previously reported mutations (Table 2). Oklahoma-ADDL-4, carried a non-synonymous mutation S194L(19) while Oklahoma-ADDL 2,3,5 carried R203K and G204R(40).

View this table:
[Table 2:](http://medrxiv.org/content/early/2020/09/18/2020.09.15.20195420/T2)

Table 2: 
Non-synonymous mutations detected and their respective amino acid changes when compared to NCBI reference genome (SARS-CoV-2, Wuhan Hu-1, NC 045512.2).

A few possible novel mutations were also identified in S, ORF3a, ORF7b, ORF8 and N genes (Table 2). A mutation in the S gene, G1167V (Glycine to Valine) was identified in Oklahoma-ADDL-4. T28I in ORF7b gene, a non-synonymous mutation (Threonine to Isoleucine) in the genome Oklahoma-ADDL-5 and in ORF8, G96R non-synonymous mutations resulting in Glycine to Arginine were noted in Oklahoma-ADDL - 2 and Oklahoma-ADDL-3. These mutations in the genome indicate presence of multiple variations of the virus in Oklahoma.

A phylogenetic tree was constructed using 9072 genomes available in the NCBI repository. The nearest neighbors to the genomes were found to be isolates reported from San Diego and Atlanta in USA and from Greece and Australia. While the Oklahoma-ADDL-1 and 5 were phylogenetically more related to isolates reported from Australia, Oklahoma-ADDL-2 and 3 were phylogenetically related to isolates reported from Greece and Atlanta, USA. Oklahoma-ADDL-4 was phylogenetically related to the isolate reported from San Diego USA. Oklahoma-ADDL-1 was identified to be in clade-GH and Oklahoma-ADDL-2,3,4,5 were identified to be in clade GR as per annotations of different marker variants in the GISAID database. The presence of multiple isolates within the Oklahoma sheds light on contagious nature of these viruses.

The genomes sequenced have been submitted to GISAID and GenBank (Table 1)

## Conclusion

We determined the sequence five SARS-COV-2 genomes, collected between March to July from Oklahoma USA, using MinION by Oxford Nanopore Technologies. Genome assembly and annotation studies identified several new mutations as well as previously reported mutations. Notably, presence of D614G mutation S gene was found in all the isolates, thus indicating presence of highly virulent strains in Oklahoma. Detection of multiple mutations in the viral genomes collected from a narrow geographic region within a few months of the Pandemic underscores the ability of SARS-CoV-2 to undergo rapid genomic alterations that can potentially influence host susceptibility, pathogenicity and virulence. Phylogenetic analysis of the virus genomes revealed high similarities in gene sequences with isolates reported from Australia, Greece, Atlanta and San Diego in USA, indicating origin of different isolates arising from multiple introductions in the state. Further mass sequencing of clinical isolates is needed to comprehensively identify genomic variations of SARS-CoV-2 in specific geographic locations.

## Data Availability

The data from the study are available in GISAID and GenBank database as mentioned in the manuscript

## Competing interests

None

## Acknowledgements

We thank all the laboratory personnel involved in COVID-19 testing at the Oklahoma Animal Disease Diagnostic Laboratory. This study was funded by FDA CVM VET-LIRN.

*   Received September 15, 2020.
*   Revision received September 15, 2020.
*   Accepted September 18, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Zheng J. 2020. SARS-CoV-2: an Emerging Coronavirus that Causes a Global Threat. Int J Biol Sci 16:1678–1685.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7150/ijbs.45053&link_type=DOI) 

2.  2.Nishiura H, Linton NM, Akhmetzhanov AR. 2020. Initial cluster of novel coronavirus (2019-nCoV) infections in Wuhan, China is consistent with substantial human-to-human transmission. Multidisciplinary Digital Publishing Institute.
    
    

3.  3.2020. Statement on the second meeting of the International Health Regulations (2005) Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV).
    
    

4.  4.2020. OSDH-COVID-19 Tracker.
    
    

5.  5.Hu B, Ge X, Wang L-F, Shi Z. 2015. Bat origin of human coronaviruses. Virol J 12:221.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12985-015-0422-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26689940&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

6.  6.Zhang T, Wu Q, Zhang Z. 2020. Pangolin homology associated with 2019-nCoV. bioRxiv 2020.02.19.950253.
    
    

7.  7.Zhang T, Wu Q, Zhang Z. 2020. Probable Pangolin Origin of SARS-CoV-2 Associated with the COVID-19 Outbreak. Curr Biol 30:1346-1351.e2.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2020.03.022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

8.  8.1.  Kielian, M, 
    2.  Mettenleiter, TC, 
    3.  Roossinck, MJBT-A
    
    Corman VM, Muth D, Niemeyer D, Drosten C. 2018. Chapter Eight - Hosts and Sources of Endemic Human Coronaviruses, p. 163–188. In Kielian, M, Mettenleiter, TC, Roossinck, MJBT-A in VR (eds.),. Academic Press.
    
    

9.  9.Annan A, Baldwin HJ, Corman VM, Klose SM, Owusu M, Nkrumah EE, Badu EK, Anti P, Agbenyega O, Meyer B. 2013. Human betacoronavirus 2c EMC/2012–related viruses in bats, Ghana and Europe. Emerg Infect Dis 19:456.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1903.121503&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23622767&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

10. 10.Graat JM, Schouten EG, Heijnen M-LA, Kok FJ, Pallast EGM, de Greeff SC, Dorigo-Zetsma JW. 2003. A prospective, community-based study on virologic assessment among elderly people with and without symptoms of acute respiratory infection. J Clin Epidemiol 56:1218–1223.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0895-4356(03)00171-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14680673&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000187582500014&link_type=ISI) 

11. 11.Mackay IM, Arden KE, Speicher DJ, O’Neil NT, McErlean PK, Greer RM, Nissen MD, Sloots TP. 2012. Co-circulation of four human coronaviruses (HCoVs) in Queensland children with acute respiratory tract illnesses in 2004. Viruses 4:637–653.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/v4040637&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22590689&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303261200010&link_type=ISI) 

12. 12.Owusu M, Annan A, Corman VM, Larbi R, Anti P, Drexler JF, Agbenyega O, Adu-Sarkodie Y, Drosten C. 2014. Human coronaviruses associated with upper respiratory tract infections in three rural areas of Ghana. PLoS One 9:e99782.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0099782&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25080241&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

13. 13.Van Elden LJR, Anton M AM, Van Alphen F, Hendriksen KAW, Hoepelman AIM, Van Kraaij MGJ, Oosterheert J-J, Schipper P, Schuurman R, Nijhuis M. 2004. Frequent detection of human coronaviruses in clinical specimens from patients with respiratory tract infection by use of a novel real-time reverse-transcriptase polymerase chain reaction. J Infect Dis 189:652–657.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/381207&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14767819&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188832900013&link_type=ISI) 

14. 14.Rambaut A, Pybus OG, Nelson MI, Viboud C, Taubenberger JK, Holmes EC. 2008. The genomic and epidemiological dynamics of human influenza A virus. Nature 453:615–619.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature06945&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18418375&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256185200035&link_type=ISI) 

15. 15.Salipante SJ, Roach DJ, Kitzman JO, Snyder MW, Stackhouse B, Butler-Wu SM, Lee C, Cookson BT, Shendure J. 2015. Large-scale genomic sequencing of extraintestinal pathogenic Escherichia coli strains. Genome Res 25:119–128.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI1LzEvMTE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMTgvMjAyMC4wOS4xNS4yMDE5NTQyMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

16. 16.van Dorp L, Acman M, Richard D, Shaw LP, Ford CE, Ormond L, Owen CJ, Pang J, Tan CCS, Boshier FAT, Ortiz AT, Balloux F. 2020. Emergence of genomic diversity and recurrent mutations in SARS-CoV-2. Infect Genet Evol 83:104351.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10/ggvz4h&link_type=DOI) 

17. 17.Wu F, Zhao S, Yu B, Chen Y-M, Wang W, Song Z-G, Hu Y, Tao Z-W, Tian J-H, Pei Y-Y. 2020. A new coronavirus associated with human respiratory disease in China. Nature 579:265–269.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2008-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

18. 18.Elbe S, Buckland-Merrett G. 2017. Data, disease and diplomacy: GISAID’s innovative contribution to global health. Glob challenges (Hoboken, NJ) 1:33–46.
    
    

19. 19.Banerjee A, Sarkar R, Mitra S, Lo M, Dutta S, Chawla-Sarkar M. 2020. The novel Coronavirus enigma: Phylogeny and mutation analyses of SARS-CoV-2 viruses circulating in India during early 2020. bioRxiv 2020.05.25.114199.
    
    

20. 20.Daniloski Z, Guo X, Sanjana NE. 2020. The D614G mutation in SARS-CoV-2 Spike increases transduction of multiple human cell types. bioRxiv 2020.06.14.151357.
    
    

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

22. 22.Lu W, Xu K, Sun B. 2010. SARS accessory proteins ORF3a and 9b and their functional analysis, p. 167–175. In Molecular Biology of the SARS-Coronavirus. Springer.
    
    

23. 23.alam I, Kamau A, Kulmanov M, Arold ST, Pain A, Gojobori T, Duarte CM. 2020. Functional pangenome analysis provides insights into the origin, function and pathways to therapy of SARS-CoV-2 coronavirus. bioRxiv 2020.02.17.952895.
    
    

24. 24.Zhang Y, Zhang J, Chen Y, Luo B, Yuan Y, Huang F, Yang T, Yu F, Liu J, Liu B, Song Z, Chen J, Pan T, Zhang X, Li Y, Li R, Huang W, Xiao F, Zhang H. 2020. The ORF8 Protein of SARS-CoV-2 Mediates Immune Evasion through Potently Downregulating MHC-I. bioRxiv 2020.05.24.111823.
    
    

25. 25.Liu DX, Fung TS, Chong KK-L, Shukla A, Hilgenfeld R. 2014. Accessory proteins of SARS-CoV and other coronaviruses. Antiviral Res 109:97–109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.antiviral.2014.06.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24995382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

26. 26.Pachetti M, Marini B, Benedetti F, Giudici F, Mauro E, Storici P, Masciovecchio C, Angeletti S, Ciccozzi M, Gallo RC. 2020. Emerging SARS-CoV-2 mutation hot spots include a novel RNA-dependent-RNA polymerase variant. J Transl Med 18:1–9.
    
    

27. 27.Leary S, Gaudieri S, Chopra A, Pakala S, Alves E, John M, Das S, Mallal S, Phillips E. 2020. Three adjacent nucleotide changes spanning two residues in SARS-CoV-2 nucleoprotein: possible homologous recombination from the transcription-regulating sequence. bioRxiv 2020.04.10.029454.
    
    

28. 28.Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. 2017. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res 27:722–736.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjI3LzUvNzIyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDkvMTgvMjAyMC4wOS4xNS4yMDE5NTQyMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

29. 29.Li H. 2018. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34:3094–3100.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty191&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29750242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

30. 30.Loman NJ, Quick J, Simpson JT. 2015. A complete bacterial genome assembled de novo using only nanopore sequencing data. Nat Methods 12:733–735.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.3444&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26076426&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

31. 31.Katoh K, Standley DM. 2013. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol Biol Evol 30:772–780.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/mst010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23329690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317002300004&link_type=ISI) 

32. 32.Katoh K, Asimenos G, Toh H. 2009. Multiple alignment of DNA sequences with MAFFT. Methods Mol Biol 537:39–64.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-1-59745-251-9_3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19378139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268335500003&link_type=ISI) 

33. 33.Wang S, Sundaram JP, Stockwell TB. 2012. VIGOR extended to annotate genomes for additional 12 different viruses. Nucleic Acids Res2012/06/04. 40:W186–W192.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gks528&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22669909&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306670900030&link_type=ISI) 

34. 34.Pickett BE, Sadat EL, Zhang Y, Noronha JM, Squires RB, Hunt V, Liu M, Kumar S, Zaremba S, Gu Z, Zhou L, Larson CN, Dietrich J, Klem EB, Scheuermann RH. 2012. ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Res 2011/10/17. 40:D593–D598.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkr859&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22006842&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298601300089&link_type=ISI) 

35. 35.Kumar S, Stecher G, Li M, Knyaz C, Tamura K. 2018. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol 35:1547–1549.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msy096&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29722887&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

36. 36.Price MN, Dehal PS, Arkin AP. 2010. FastTree 2 – Approximately Maximum-Likelihood Trees for Large Alignments. PLoS One 5:e9490.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0009490&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20224823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

37. 37.Angeletti S, Benvenuto D, Bianchi M, Giovanetti M, Pascarella S, Ciccozzi M. 2020. COVID-2019: The role of the nsp2 and nsp3 in its pathogenesis. J Med Virol 92:584–588.
    
    

38. 38.Banerjee S, Seal S, Dey R, Mondal KK, Bhattacharjee P. 2020. Mutational spectra of SARS-CoV-2 orf1ab polyprotein and Signature mutations in the United States of America. bioRxiv 2020.05.01.071654.
    
    

39. 39.Graham RL, Sparks JS, Eckerle LD, Sims AC, Denison MR. 2008. SARS coronavirus replicase proteins in pathogenesis. Virus Res2007/03/29. 133:88–100.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.virusres.2007.02.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17397959&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F18%2F2020.09.15.20195420.atom) 

40. 40.Laha S, Chakraborty J, Das S, Manna SK, Biswas S, Chatterjee R. 2020. Characterizations of SARS-CoV-2 mutational profile, spike protein stability and viral transmission. bioRxiv 2020.05.03.066266.
    
    

41. 41.Grubaugh ND, Hanage WP, Rasmussen AL. 2020. Making sense of mutation: what D614G means for the COVID-19 pandemic remains unclear. Cell.
    
    

42. 42.Korber B, Fischer WM, Gnanakaran S, Yoon H, Theiler J, Abfalterer W, Hengartner N, Giorgi EE, Bhattacharya T, Foley B, Hastie KM, Parker MD, Partridge DG, Evans CM, Freeman TM, de Silva TI, McDanal C, Perez LG, Tang H, Moon-Walker A, Whelan SP, LaBranche CC, Saphire EO, Montefiori DC, Group SC-19 G. 2020. Tracking changes in SARS-CoV-2 Spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell.
    
    

43. 43.Issa E, Merhi G, Panossian B, Salloum T, Tokajian S. 2020. SARS-CoV-2 and ORF3a: Non-Synonymous Mutations and Polyproline Regions. bioRxiv 2020.03.27.012013.