Abstract

The recent outbreak of COVID-19 caused by a new zoonotic origin coronavirus (SARS-CoV-2 or 2019-nCoV) has sound the alarm for the potential spread of epidemic coronavirus crossing species. With the urgent needs to assist disease control and to provide invaluable scientific information, we developed the coronavirus database (CoVdb), an online genomic, proteomic and evolutionary analysis platform. CoVdb has brought together genomes of more than 5000 coronavirus strains, which were collected from 1941 to 2020, in more than 60 countries and in hosts belonging to more than 30 species, ranging from fish to human. CoVdb presents comprehensive genomic information, such as gene function, subcellular localization, topology and protein structure. To facilitate coronavirus research, CoVdb also provides flexible search approaches and online tools to view and analyze protein structure, to perform multiple alignments, to automatically build phylogenetic trees and to carry on evolutionary analyses. CoVdb can be accessed freely at http://covdb.popgenetics.net. Hopefully, it will accelerate the progress to develop medicines or vaccines to control the pandemic of COVID-19.

Introduction

Coronaviridae is a group of positive-sense, single-strand RNA viruses with a likely ancient origin, and human coronavirus repeatedly emerged during the past hundred years (1). Coronaviruses are classified into four distinct genera: alpha and beta coronavirus mainly infect mammals, whereas gamma and delta coronavirus circulate more often in avian hosts (2). As a potential dangerous zoonotic disease, the previous outbreaks of respiratory syndrome-related coronavirus (SARS-CoV) and Middle East respiratory syndrome-related coronavirus (MERS-CoV) have plagued the general public and researchers in the past years (3). Recently, a novel coronavirus, which may originated from wild animals, was first identified in Wuhan City, China. The virus is the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), also named 2019 novel coronavirus (2019-nCoV), causing coronavirus disease-2019 (COVID-19). Till now, it has resulted in more than 16 million confirmed infections worldwide (4, 5) with the number of infection cases still increasing. Although we have knowledge and experience in the virology, diagnosis, clinical characteristics and other aspects related to SARS-CoV and MERS-CoV, there are many unanswered questions about the new emerging SARS-CoV-2. The new outbreak coronavirus strongly reminds the continuous threat of zoonotic diseases caused by coronavirus to global health security. Sharing experience and knowledge across disciplines in historical and global scale should provide invaluable scientific knowledge to fight against the threat of coronavirus.

The aim of the construction of CoVdb is to provide coronavirus knowledge, to contribute to global coronavirus research, especially for the investigation of the emerging SARS-CoV-2. For previous works, ViPR (6) and ViralZone (7) are general data resources and are lack of analysis tools in population genomics and evolution. Different from those databases, CoVdb is specially designed for coronavirus. It combines, compares and annotates all published coronavirus genomes up to date (8–19). Compared to 2019nCoV (20), CoVdb provides more population genetic analysis information and contains several online sequence analysis tools. The new developed database provides the convenience for the identification of gene function and identity among Coronaviridae genomes. CoVdb provides information on subcellular location, function, protein topology, as well as population level through analyses. We will be dedicated to keep updating the genomic information and optimizing the database.

Materials and methods

Data processing

Coronavirus sequences and annotations are downloaded from the NCBI nucleotide database (21). We chose records with complete genomes. For newly sequenced strains without open reading frame (ORF) annotation, we did annotation through mapping known coronavirus proteins to the genome by GeneWise (22). We verified the quality of these proteins by known proteins, documented coronavirus proteins in NCBI and kept predicted proteins that have an identity >0.5 and a coverage >50%. We renamed all human isolates in the format of ‘Human_name_accession’ (‘name’ is 2019-nCoV (SARS-CoV-2), SARS, MERS or other human coronavirus strain names) and all nonhuman isolates in the format of ‘host_accession’ (‘host’ is bat, camel, cow or other coronavirus hosts). Accession is the strain’s GeneBank ID. According to previous methods to cluster homologous genes (23, 24), we grouped coronavirus proteins into 628 unified clusters by CD-HIT (25) with an identity >50% and a coverage >80%. We made classification for all documented coronavirus strains based on NCBI’s annotation on taxonomy. For an unclassified strain, we run BLAST with the strain’s genome against a ‘reference set’ of some lineage, such as Sarbecovirus, Setracovirus and so on. According to the BLAST output, we attribute the strain to a closest lineage. We list the sequenced strains documented in GISAID (26), and users are instructed to visit GISAID if they want to get the actual data.

We wrote Perl scripts to automatically BLAST coronavirus protein sequences against the UniProt database (27). We filtered out hits with an E-value <0.05 and only kept the one that has the highest score in alignments. Using matched UniProt accession numbers, we retrieved detailed proteomic information from UniProt. We looked for each gene’s possible protein 3D structure counterparts in the Protein Data Bank (28, 29) in the same way.

We did subcellular localization prediction for coronavirus gene clusters using an online tool MSLVP (30). We used TMHMM 2.0 (31) to predict the transmembrane helixes within protein sequences for all coronavirus genes and converted output images into PNG format by Magick (www.imagemagick.org).

We utilized CUDA ClustalW (32) to perform multiple alignments of all documented coronavirus genomes or proteins. The results are used to build phylogenetic trees. We aligned genomes by LASTZ (33), aligned CDS or protein sequences by MUSCLE (34, 35) and built phylogenetic trees by FastTree 2.1 (36). In the detection of selection signals, we performed sliding widow analyses with a window size 200 bp and a step size 50 bp. Population genetic tests were performed by VariScan 2.0 (37, 38) and SweepFineder2 (39). Due to the lack of software able to use multiple alignment to calculate the Fixation Index (Fst) (40), we wrote Perl scripts to calculate Fst according to reported algorithms (41).

A semi-automatic pipeline was developed to update the database, considering the fast increase of coronavirus’ sequences in NCBI. For specifics, the pipeline first refreshes the strain list according to NCBI’s coronavirus list (requiring complete genome), then updates corresponded annotation and finally make revision for some overall analysis results.

The development of engines, interfaces and tools

Similar with what we have done for SGID (42) and ASFVdb (43), we made the web interface of CoVdb based on SWAV (44). CoVdb also incorporates MSAViewer (45) to display multiple alignments and phylotree.js (46) to show phylogenetic trees. To improve the display of virus data, we made changes in these two softwares. We changed parameters to fit virus’s dense gene arrangements and added links within diagrams. The search engine is written by PHP integrated with SQL, BLAT (47) and NCBI BLAST (48). The protein 3D viewer in CoVdb uses the libraries of iCn3D (49). For ‘Pop Analyzer’ and ‘Aln Browser’, the background codes were written basing on C, Perl, VariScan 2.0, SweepFineder2, FastTree 2.1 and LASTZ. ‘Phylo Tree’ shows phylogenetic trees built by genomes or coronavirus’ major proteins, the Orf1ab polyprotein (Orf1), the spike glycoprotein (S), the envelope protein (E), the membrane protein (M) and the nucleocapsid protein (N).

Results and discussion

Data and information

CoVdb extensively collects published coronavirus data and have taken in genomes of 5709 strains after the update in 22 May 2020. The strains were collected from 32 organisms and in the years from 1941 to present, 2020 (Figure 1). A total of 3414 (59.8%) in CoVdb are human isolates and 217 (3.8%) are bat isolates, which are referred as the possible source of human coronavirus (5051). Porcine coronavirus also take a big percentage (945, 16.6%) and coronavirus used to make damages in the pig industry (52). The number of documented human isolates varied in years, and there are three peaks that reflect the outbreaks of SARS-CoV in 2003, MERS-CoV in 2014–2015 and SARS-CoV-2 in 2019–2020 separately. Using all documented coronavirus genomes in CoVdb, we generated a phylogenetic tree (Figure 2A), from which we observed that the nearest nonhuman isolate to 2019-nCoV is Bat_MN996532 (Bat-CoV-RaTG13), isolated from Rhinolophus affinis, a species of bat in the Rhinolophidae family. Strains isolated from pangolins are also in the vicinity of SARS-CoV-2. Pangolin was once considered as a potential intermediate host of SARS-CoV-2 (53). We developed search tools to enable users to search in the big phylogenetic tree (Figure 2B).

The distribution of documented coronavirus strains in CoVdb according to collection date (X-axis) and hosts (colored by different colors). Y is the number of coronavirus isolated from some organism. Red triangles points to peaks in the distribution of human coronavirus in years.
Figure 1.

The distribution of documented coronavirus strains in CoVdb according to collection date (X-axis) and hosts (colored by different colors). Y is the number of coronavirus isolated from some organism. Red triangles points to peaks in the distribution of human coronavirus in years.

(A) Partial display of the phylogenetic tree built by all coronavirus genomes documented in CoVdb. Red numbers are marginal likelihoods. (B) Snapshot showing that users can search a strain by name in a phylogenetic tree. Both A and B center on the split of Bat_MN996532 and 2019-nCoV (SARS-CoV-2).
Figure 2.

(A) Partial display of the phylogenetic tree built by all coronavirus genomes documented in CoVdb. Red numbers are marginal likelihoods. (B) Snapshot showing that users can search a strain by name in a phylogenetic tree. Both A and B center on the split of Bat_MN996532 and 2019-nCoV (SARS-CoV-2).

In average, there are 5–14 possible ORFs or genes in one coronavirus strain. We grouped homologous coronavirus genes (requiring identity >0.5 and coverage >0.8) into 628 clusters (for details, see Materials and Methods). This number indicates that the differentiation or diversity within coronavirus strains is not low. For these, we still performed a subcellular localization analysis for the 628 clusters to predict their roles in infection, although the structure of coronavirus is not complex. Based on prediction only, 21% (133 items) are predicted to be located in the host nucleus or host cytoplasm, while 40% (250 items) are predicted to be membrane proteins (Figure S1). CoVdb has included more than 50 000 function annotations and more than 300 000 GO records. Using WEGO (54), we found coronavirus genes enrich in the membrane (Figure S2). We searched for possible protein 3D structure for coronavirus genes in the Protein Data Bank (2829) and found more than 3 000 000 mappings with an E-value <0.05 and a coverage >50%.

For all coronavirus strains, using nine representative human coronavirus genomes as the reference, we did sliding window analyses on Pi (55), Tajima’s D (56), composite likelihood ratio (CLR) (5758) and Fst (40). For Pi, Tajima’s D and CLR, the target group are strains that belong to COVID-19, MERS, SARS or other human coronavirus diseases. We also did the same thing for human isolates, bat isolates and isolates of other hosts documented in CoVdb. Fst is between human coronavirus and one nonhuman coronavirus. All these data can be viewed in the genome browser.

Interface and analysis tools

The genome browser (GBrowser) in CoVdb follows a style with gene segments followed by analysis tracks (CLR, Pi, Tajima’s D and Fst). Users can view population genetic tracks of strains belonging to one or more hosts, such as avian, bat and so on. They can also view that of strains belonging to one or more specific diseases, such as COVID-19, SARS and 229E. Users can select by checkboxes. At the top of GBrowser, there are browsing tools, such as search by inputting a chromosome position, zoom in/out and position movements. It provides notes for whether a gene of ORF is already annotated in GeneBank or newly predicted in the top of the genome browser page. In addition to basic information, CoVdb shows a gene in function, subcellular localization, topology and protein structure. The search engine in CoVdb is powerful and supports fuzzy search, search with taxonomy (such as Alpha, Beta, Sarbecovirus, etc.), filtering, sorting, BLAT and BLAST. CoVdb also allows to search by cell location. For personalized analyses, CoVdb is able to provide gene links if inputting a list of chromosome positions or gene accessions.

CoVdb has tools to facilitate some specific use in coronavirus research, such as tracing origination, vaccine or drug design. In the tool ‘Protein’, the protein structure information is listed, where users can view the overlapped amino acids of a coronavirus protein in the 3D structure counterpart and do online protein structure analysis by an embedded application iCn3D (49) (Figure S3A, C). Users also can BLAST a protein sequence against CoVdb’s protein data and view the mapped region in a protein 3D view (Figure S3B, D). The tool ‘Aln Browser’, alignment browser, allows users to retrieve the multiple alignment of two or more strains at some position and build a phylogenetic tree using the alignment (Figure 3). With the tool ‘Pop Analyzer’, population genetics analyzer, users can do personalized online sliding window analyses. Users can choose the window size, the step size, the target region and one or more population genetic tests (Figure S4). Users can jump to the search list and select multiple strains to do analysis in ‘Aln Browser’ or in ‘Pop Analyzer’. For initial users, they can learn to use each analysis tool by clicking the button ‘Demo’. ‘Phylo Tree’ is a tool to view and search in phylogenetic trees made by genomic or proteomic sequences (Figure 2B). Users also can go to the GBrowser page by clicking the name of one strain.

A snapshot displaying the usage of ‘Aln Browser’, where users need to select the reference strain, the start position, the end position and the strains to be put in alignment. If clicking on the button ‘Retrieve Alignment’, a multiple alignment of selected strains will be shown below. If clicking on ‘Make Tree’, a phylogenetic tree will be built basing on the alignment and shown at the bottom.
Figure 3.

A snapshot displaying the usage of ‘Aln Browser’, where users need to select the reference strain, the start position, the end position and the strains to be put in alignment. If clicking on the button ‘Retrieve Alignment’, a multiple alignment of selected strains will be shown below. If clicking on ‘Make Tree’, a phylogenetic tree will be built basing on the alignment and shown at the bottom.

View population genetic tracks of the spike glycoprotein in GBrowser

The spike glycoprotein (S protein) plays a key role in the infection of COVID-19 while the receptor-binding domain (RBD) (59) is the region in S protein to interact with the human protein ACE2 (60). In CoVdb’s GBrowser (Figure S5), we observed that S protein is highly conserved within SARS-CoV-2/2019-nCoV strains (CLR, Pi and Tajima’s D are nearly fixed to zero for most points). In comparison, for SARS-CoV, there are CLR peaks and variations in Pi and Tajima’s D. For human coronavirus, the region near RBD is highly conserved (with a lowland of Pi and Tajima’s D, while Tajima’s D is −0.81 in median) compared to the region far away from RBD (with a plateau of Pi and Tajima’s D, while Tajima’ D is 1.39 in median). However, we did not observe a similar pattern for bat coronavirus. Human and bat coronavirus are of CLR peaks at different sites. We also observed Fst peaks between human and bat coronavirus. These indicated that the evolution is different not only between SARS-CoV-2/2019-nCoV and SARS-CoV but also between human and bat coronavirus.

Conclusion

Dedicated to assist researchers to combat the pandemic of COVID-19 and to provide a more specialized platform for coronavirus, we comprehensively gathered data and systematically constructed the coronavirus database, CoVdb. In the database, researchers can conveniently retrieve genomic or gene information of coronavirus and do online analyses in comparative genomics, protein structure and evolutionary biology. With the help of this database, we have successfully developed test strips able to detect SARS-CoV-2 (unpublished). With the increase of the number of sequenced coronavirus genomes, we will provide continuous update and maintenance of the database in the future. Hopefully, this database will play more important roles in fighting against the infection of coronavirus in the future.

Acknowledgements

We gratefully acknowledge the submitting and the originating laboratories where genetic sequence data were generated and shared via NCBI and the GISAID Initiative.

Supplementary data

Supplementary data are available at Database Online.

Funding

National Key Research and Development Program (2019YFC1604600), the National Natural Science Foundation of China (31200941), the Fundamental Research Funds for the Central Universities (106112016CDJXY290002) and the National Natural Science Foundation of HeBei Province (19226631D).

Conflict of interest

We declare no competing interests.

Data availability

All CoVdb data are publicly and freely accessible at http://covdb.popgenetics.net. Feedback on any aspect of the CoVdb and discussions of coronavirus gene annotations are welcome by email to zhuzl@cqu.edu.cn or mg@cau.edu.cn.

Author contributions

Z.Z. developed the web interface of the database, collected and compiled the data. K.M. performed the analyses. Z.Z. and G.L. are responsible for updating and maintenance. Z.Z. and G.M. conceived the idea, coordinated the project and wrote the manuscript.

References

1.

Forni
D.
,
Cagliani
R.
,
Clerici
M.
et al.  (
2017
)
Molecular evolution of human coronavirus genomes
.
Trends Microbiol.
,
25
,
35
48
.

2.

Wertheim
J.O.
,
Chu
D.K.
,
Peiris
J.S.
et al.  (
2013
)
A case for the ancient origin of coronaviruses
.
J. Virol.
,
87
,
7039
7045
.

3.

de Wit
E.
,
van Doremalen
N.
,
Falzarano
D.
et al.  (
2016
)
SARS and MERS: recent insights into emerging coronaviruses
.
Nat. Rev. Microbiol.
,
14
,
523
534
.

4.

Lu
H.
,
Stratton
C.W.
and
Tang
Y.W.
(
2020
)
Outbreak of pneumonia of unknown etiology in Wuhan China: the mystery and the miracle
.
J. Med. Virol.
,
92
,
401
402
.

5.

Hui
D.S.
,
E
I.A.
,
Madani
T.A.
et al.  (
2020
)
The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health - the latest 2019 novel coronavirus outbreak in Wuhan, China.
Int. J. Infect. Dis.
,
91
,
264
266
.

6.

Pickett
B.E.
,
Greer
D.S.
,
Zhang
Y.
et al.  (
2012
)
Virus pathogen database and analysis resource (ViPR): a comprehensive bioinformatics database and analysis resource for the coronavirus research community
.
Viruses
,
4
,
3209
3226
.

7.

Hulo
C.
,
de Castro
E.
,
Masson
P.
et al.  (
2011
)
ViralZone: a knowledge resource to understand virus diversity
.
Nucleic Acids Res.
,
39
,
D576
D582
.

8.

Boursnell
M.E.
,
Brown
T.D.
,
Foulds
I.J.
et al.  (
1987
)
Completion of the sequence of the genome of the coronavirus avian infectious bronchitis virus
.
J. Gen. Virol.
,
68
,
57
77
.

9.

Coley
S.E.
,
Lavi
E.
,
Sawicki
S.G.
et al.  (
2005
)
Recombinant mouse hepatitis virus strain A59 from cloned, full-length cDNA replicates to high titers in vitro and is fully pathogenic in vivo
.
J. Virol.
,
79
,
3097
3106
.

10.

St-Jean
J.R.
,
Jacomy
H.
,
Desforges
M.
et al.  (
2004
)
Human respiratory coronavirus OC43: genetic stability and neuroinvasion
.
J. Virol.
,
78
,
8824
8834
.

11.

Chouljenko
V.N.
,
Lin
X.Q.
,
Storz
J.
et al.  (
2001
)
Comparison of genomic and predicted amino acid sequences of respiratory and enteric bovine coronaviruses isolated from the same animal with fatal shipping pneumonia
.
J. Gen. Virol.
,
82
,
2927
2933
.

12.

van Boheemen
S.
,
de Graaf
M.
,
Lauber
C.
et al.  (
2012
)
Genomic characterization of a newly discovered coronavirus associated with acute respiratory distress syndrome in humans
.
mBio
,
3
,
e00473
12
.

13.

Vlasova
A.N.
,
Halpin
R.
,
Wang
S.
et al.  (
2011
)
Molecular characterization of a new species in the genus Alphacoronavirus associated with mink epizootic catarrhal gastroenteritis
.
J. Gen. Virol.
,
92
,
1369
1379
.

14.

Marra
M.A.
,
Jones
S.J.
,
Astell
C.R.
et al.  (
2003
)
The genome sequence of the SARS-associated coronavirus
.
Science
,
300
,
1399
1404
.

15.

Woo
P.C.
,
Lau
S.K.
,
Chu
C.M.
et al.  (
2005
)
Characterization and complete genome sequence of a novel coronavirus, coronavirus HKU1, from patients with pneumonia
.
J. Virol.
,
79
,
884
895
.

16.

Tang
X.C.
,
Zhang
J.X.
,
Zhang
S.Y.
et al.  (
2006
)
Prevalence and genetic diversity of coronaviruses in bats from China
.
J. Virol.
,
80
,
7481
7490
.

17.

Lau
S.K.
,
Woo
P.C.
,
Li
K.S.
et al.  (
2007
)
Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome
.
Virology
,
367
,
428
439
.

18.

Chu
D.K.
,
Peiris
J.S.
,
Chen
H.
et al.  (
2008
)
Genomic characterizations of bat coronaviruses (1A, 1B and HKU8) and evidence for co-infections in miniopterus bats
.
J. Gen. Virol.
,
89
,
1282
1287
.

19.

Woo
P.C.
,
Wang
M.
,
Lau
S.K.
et al.  (
2007
)
Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features
.
J. Virol.
,
81
,
1574
1585
.

20.

Zhao
W.M.
,
Song
S.H.
,
Chen
M.L.
et al.  (
2020
)
The 2019 novel coronavirus resource
.
Yi Chuan
,
42
,
212
221
.

21.

Sayers
E.W.
,
Beck
J.
,
Brister
J.R.
et al.  (
2020
)
Database resources of the National Center for Biotechnology Information
.
Nucleic Acids Res.
,
48
,
D9
D16
.

22.

Madeira
F.
,
Park
Y.M.
,
Lee
J.
et al.  (
2019
)
The EMBL-EBI search and sequence analysis tools APIs in 2019
.
Nucleic Acids Res.
,
47
,
W636
W641
.

23.

Yue
H.
,
Shu
D.
,
Wang
M.
et al.  (
2018
)
Genome-wide identification and expression analysis of the HD-zip gene family in wheat (Triticum aestivum L.)
.
Genes
,
9
,
70
.

24.

She
R.
,
Chu
J.S.
,
Wang
K.
et al.  (
2009
)
GenBlastA: enabling BLAST to identify homologous gene sequences
.
Genome Res.
,
19
,
143
149
.

25.

Li
W.
and
Godzik
A.
(
2006
)
Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences
.
Bioinformatics
,
22
,
1658
1659
.

26.

Shu
Y.
and
McCauley
J.
(
2017
)
GISAID: global initiative on sharing all influenza data - from vision to reality
.
Euro Surveill.
,
22
, 30494.

27.

Patient
S.
,
Wieser
D.
,
Kleen
M.
et al.  (
2008
)
UniProtJAPI: a remote API for accessing UniProt data
.
Bioinformatics
,
24
,
1321
1322
.

28.

Berman
H.M.
,
Westbrook
J.
,
Feng
Z.
et al.  (
2000
)
The Protein Data Bank
.
Nucleic Acids Res.
,
28
,
235
242
.

29.

Burley
S.K.
,
Berman
H.M.
,
Bhikadiya
C.
et al.  (
2019
)
RCSB Protein Data Bank: biological macromolecular structures enabling research and education in fundamental biology, biomedicine, biotechnology and energy
.
Nucleic Acids Res.
,
47
,
D464
D474
.

30.

Thakur
A.
,
Rajput
A.
and
Kumar
M.
(
2016
)
MSLVP: prediction of multiple subcellular localization of viral proteins using a support vector machine
.
Mol. Biosyst.
,
12
,
2572
2586
.

31.

Krogh
A.
,
Larsson
B.
,
von Heijne
G.
et al.  (
2001
)
Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes
.
J. Mol. Biol.
,
305
,
567
580
.

32.

Hung
C.L.
,
Lin
Y.S.
,
Lin
C.Y.
et al.  (
2015
)
CUDA ClustalW: an efficient parallel algorithm for progressive multiple sequence alignment on multi-GPUs
.
Comput. Biol. Chem.
,
58
,
62
68
.

33.

Harris
R.S.
(
2007
)
Improved pairwise alignment of genomic DNA
.
Ph.D. Thesis
.
Pennsylvania State University
.

34.

Edgar
R.C.
(
2004
)
MUSCLE: a multiple sequence alignment method with reduced time and space complexity
.
BMC Bioinform.
,
5
, 113.

35.

Edgar
R.C.
(
2004
)
MUSCLE: multiple sequence alignment with high accuracy and high throughput
.
Nucleic Acids Res.
,
32
,
1792
1797
.

36.

Price
M.N.
,
Dehal
P.S.
and
Arkin
A.P.
(
2010
)
FastTree 2—approximately maximum-likelihood trees for large alignments
.
PLoS One
,
5
, e9490.

37.

Hutter
S.
,
Vilella
A.J.
and
Rozas
J.
(
2006
)
Genome-wide DNA polymorphism analyses using VariScan
.
BMC Bioinform.
,
7
, 409.

38.

Vilella
A.J.
,
Blanco-Garcia
A.
,
Hutter
S.
et al.  (
2005
)
VariScan: Analysis of evolutionary patterns from large-scale DNA sequence polymorphism data
.
Bioinformatics
,
21
,
2791
2793
.

39.

DeGiorgio
M.
,
Huber
C.D.
,
Hubisz
M.J.
et al.  (
2016
)
SweepFinder2: increased sensitivity, robustness and flexibility
.
Bioinformatics
,
32
,
1895
1897
.

40.

Holsinger
K.E.
and
Weir
B.S.
(
2009
)
Genetics in geographically structured populations: defining, estimating and interpreting F(ST)
.
Nat. Rev. Genet.
,
10
,
639
650
.

41.

Fumagalli
M.
,
Vieira
F.G.
,
Korneliussen
T.S.
et al.  (
2013
)
Quantifying population genetic differentiation from next-generation sequencing data
.
Genetics
,
195
,
979
992
.

42.

Zhu
Z.
,
Guan
Z.
,
Liu
G.
et al.  (
2019
)
SGID: a comprehensive and interactive database of the silkworm
.
Database
,
2019
, baz134.

43.

Zhu
Z.
and
Meng
G.
(
2019
)
ASFVdb: An integrative resource for genomics and proteomics analyses of African swine fever
.
Database, 2019, baaa023
.

44.

Zhu
Z.
,
Wang
Y.
,
Zhou
X.
et al.  (
2020
)
SWAV: a web-based visualization browser for sliding window analysis
.
Sci. Rep.
,
10
, 149.

45.

Yachdav
G.
,
Wilzbach
S.
,
Rauscher
B.
et al.  (
2016
)
MSAViewer: interactive JavaScript visualization of multiple sequence alignments
.
Bioinformatics
,
32
,
3501
3503
.

46.

Shank
S.D.
,
Weaver
S.
and
Kosakovsky Pond
S.L.
(
2018
)
phylotree.js - a JavaScript library for application development and interactive data visualization in phylogenetics
.
BMC Bioinform.
,
19
, 276.

47.

Kent
W.J.
(
2002
)
BLAT—the BLAST-like alignment tool
.
Genome Res.
,
12
,
656
664
.

48.

Johnson
M.
,
Zaretskaya
I.
,
Raytselis
Y.
et al.  (
2008
)
NCBI BLAST: a better web interface
.
Nucleic Acids Res.
,
36
,
W5
W9
.

49.

Wang
J.
,
Youkharibache
P.
,
Zhang
D.
et al.  (
2020
)
iCn3D, a web-based 3D viewer for sharing 1D/2D/3D representations of biomolecular structures
.
Bioinformatics
,
36
,
131
135
.

50.

Lu
R.
,
Zhao
X.
,
Li
J.
et al.  (
2020
)
Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding
.
Lancet
,
395
,
565
574
.

51.

Hu
B.
,
Zeng
L.P.
,
Yang
X.L.
et al.  (
2017
)
Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus
.
PLoS Pathog.
,
13
, e1006698.

52.

Niederwerder
M.C.
and
Hesse
R.A.
(
2018
)
Swine enteric coronavirus disease: a review of 4 years with porcine epidemic diarrhoea virus and porcine deltacoronavirus in the United States and Canada
.
Transbound Emerg. Dis.
,
65
,
660
675
.

53.

Xiao
K.
,
Zhai
J.
,
Feng
Y.
et al.  (
2020
)
Isolation and characterization of 2019-nCoV-like coronavirus from Malayan Pangolins
. bioRxiv. doi:

54.

Ye
J.
,
Fang
L.
,
Zheng
H.
et al.  (
2006
)
WEGO: a web tool for plotting GO annotations
.
Nucleic Acids Res.
,
34
,
W293
W297
.

55.

Watterson
G.A.
(
1975
)
On the number of segregating sites in genetical models without recombination
.
Theor. Popul. Biol.
,
7
,
256
276
.

56.

Tajima
F.
(
1989
)
Statistical method for testing the neutral mutation hypothesis by DNA polymorphism
.
Genetics
,
123
,
585
595
.

57.

Nielsen
R.
,
Williamson
S.
,
Kim
Y.
et al.  (
2005
)
Genomic scans for selective sweeps using SNP data
.
Genome Res.
,
15
,
1566
1575
.

58.

Zhu
L.
and
Bustamante
C.D.
(
2005
)
A composite-likelihood approach for detecting directional selection from DNA sequence data
.
Genetics
,
170
,
1411
1421
.

59.

Chen
Y.
,
Guo
Y.
,
Pan
Y.
et al.  (
2020
)
Structure analysis of the receptor binding of 2019-nCoV
.
Biochem. Biophys. Res. Commun
,
525
,
135
140
.

60.

Cai
G.
(
2020
)
Bulk and single-cell transcriptomics identify tobacco-use disparity in lung gene expression of ACE2, the receptor of 2019-nCov
.
medRxiv
. doi:

Author notes

Citation details: Zhu,Z., Meng,K., Liu,G., et al. A database resource and online analysis tools for coronaviruses on a historical and global scale. Database (2020) Vol. 00: article ID baaa070; doi:10.1093/database/baaa070

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data