This is a limited proof of concept to search for research data, not a production system.

Search the MIT Libraries

Title: roblanf/sarscov2phylo: 18-08-20

Type Software roblanf (2020): roblanf/sarscov2phylo: 18-08-20. Zenodo. Software. https://zenodo.org/record/3995110

Author: roblanf ;

Links

Summary

Citation and reuse

Please cite this release as:

Lanfear, Rob (2020). A global phylogeny of SARS-CoV-2 from GISAID data, including sequences deposited up to 18-August-2020. Zenodo DOI: 10.5281/zenodo.3958883

You can visit that DOI here:

Please note - you cannot publish papers that use this tree without following the GISAID data sharing and attribution rules. These rules are important - they protect the data uploaders, and create trust in a global system of data sharing with potentially vast public health benefits. By building and maintaining trust we ensure that people keep sharing their data, and that the public health benefits keep flowing. I do not want the existence of this tree to be some kind of attribution laundering service (e.g. where people feel free to use the tree without following the GISAID data sharing rules), so please don't use it in that way. For example, if you are going to interpret other people's data from GISAID and publish the results, including by using this tree, you should get in touch with the people that submitted the data. The code in this repo is covered by the GNU license, and you can use that however you like.

Details

The trees in this release were generated with the following command line:

bash global_tree_gisaid_start_tree.sh -i gisaid_hcov-19_2020_08_21_04.fasta -o global.fa -s ft_SH_16-08-20.tree -t 54

The raw sequence file contains all available SARS-CoV-2 genomes in GISAID available on the 18th of August 2020, determined by the 'submission date' filter on GISAID.

The ZIP file contains the code necessary to reproduce the trees themselves, and the README in the zip file also describes the methods used in detail. I also include the tree itself here so that it can be easily downloaded without downloading the entire repo. The file ' ft_SH_16-8-20.tree' is the 'ft_SH.tree' file from the 16-8-20 release.

The lnL of the final tree is: -513889.945

Filtering statistics sequences downloaded from GISAID 54551 // alignment stats of global alignment Alignment number: 1 Format: aligned FASTA Number of sequences: 53041 Alignment length: 29903 Total # residues: 1582325418 Smallest: 29030 Largest: 29903 Average length: 29832.1 Average identity: 100% // alignment stats of global alignment after masking sites Alignment number: 1 Format: aligned FASTA Number of sequences: 53041 Alignment length: 29903 Total # residues: 1573527101 Smallest: 28961 Largest: 29675 Average length: 29666.2 Average identity: 100% // alignment stats after filtering out short/ambiguous sequences Alignment number: 1 Format: aligned FASTA Number of sequences: 52758 Alignment length: 29903 Total # residues: 1565142267 Smallest: 29054 Largest: 29675 Average length: 29666.4 Average identity: 100% // alignment stats of global alignment after trimming sites that are >50% gaps Alignment number: 1 Format: aligned FASTA Number of sequences: 52758 Alignment length: 29661 Total # residues: 1561355620 Smallest: 28437 Largest: 29661 Average length: 29594.7 Average identity: 100% // After filtering sequences with TreeShrink Type: Phylogram #nodes: 94593 #leaves: 52747 #dichotomies: 39994 #leaf labels: 52747 #inner labels: 37785 Number of new sequences added this iteration 363 alignment_names_new.txt Notable changes to the scripts in this release None Notable aspects of the trees The likelihood of this tree is lower than of the previous release. This is quite interesting, and despite this tree having ~300 new sequences in it. It suggests the SPR moves found a change that makes substantial improvements to the likelihood of the tree, potentially (though this would need confirmation) underlining the benefit of an iterative approach.

More information

  • DOI: 10.5281/zenodo.3995110

Dates

  • Publication date: 2020
  • Issued: August 21, 2020

Rights

  • info:eu-repo/semantics/openAccess Open Access

Much of the data past this point we don't have good examples of yet. Please share in #rdi slack if you have good examples for anything that appears below. Thanks!

Format

electronic resource

Relateditems

DescriptionItem typeRelationshipUri
IsSupplementTohttps://github.com/roblanf/sarscov2phylo/tree/18-08-20
IsVersionOfhttps://doi.org/10.5281/zenodo.3958883
IsPartOfhttps://zenodo.org/communities/covid-19
IsPartOfhttps://zenodo.org/communities/zenodo