{"id":344,"date":"2019-01-17T08:03:41","date_gmt":"2019-01-17T07:03:41","guid":{"rendered":"http:\/\/wp.unil.ch\/phylo\/?page_id=344"},"modified":"2019-01-17T08:03:41","modified_gmt":"2019-01-17T07:03:41","slug":"coev","status":"publish","type":"page","link":"https:\/\/wp.unil.ch\/phylo\/tools\/coev\/","title":{"rendered":"Coev"},"content":{"rendered":"
\n

Coev<\/span>\u00a0is a software that predicts\u00a0coevolving<\/span>\u00a0positions and their evolutionary profile based on the aligned sequences and a phylogenetic tree.<\/p>\n<\/div>\n

\n

Overview:<\/h3>\n

The software is the implementation of the model described in Dib et al. 2013, where we present a new substitution model that describes the\u00a0coevolving<\/span>\u00a0process along the tree and reconstruct the ancestral states of co-evolving pairs of positions. We propose a new probabilistic Markov model that not only identifies\u00a0coevolving<\/span>\u00a0positions but also estimates the associated\u00a0coevolving<\/span>\u00a0profile. The necessity to estimate profiles emerges when considering the nucleotide alphabet that induces 192 different\u00a0coevolving<\/span>\u00a0profiles and amino acid alphabet that induces a bigger number of profiles. The Coev Markov model is based on 16 states instantaneous matrix where each state represented the a combination. The instantaneous rate matrix Q contains 4 continuous parameters and a discrete parameter. The implementation of the model is written in C (Dennis Ritchie, 1969) with a user-friendly interface. Given an alignment file, the phylogenetic tree and a pair of positions the user will be able to assess the score of coevolution and estimate the profile with either the maximum likelihood or Bayesian framework. The user can also change within this interface the default values of the priors.<\/p>\n

Download:<\/h3>\n

To download, please click on the following\u00a0link<\/a>. A manual can be found\u00a0here<\/a>, which also explain how to use the Coev-web platform that will be released soon.<\/p>\n

System requirements: The software Coev uses several external tools that should be installed. We optimized the implementation using Lapack (Linear Algebra PACKage) and nlopt (library for nonlinear optimization) downloadable at\u00a0https:\/\/ab-initio.mit.edu\/wiki\/index.php\/NLopt<\/a>\u00a0and tested it on Linux and mac operating systems. Additionally, for the GUI interface we used Tcl\/Tk and wish softwares.<\/p>\n

Instructions for running Coev with a GUI:<\/h3>\n

Unzip Coev.zip in a directory of your choice.<\/p>\n

Once all external software installed, to compile run:make<\/p>\n

Then to execute run: wish Coev.tcl<\/p>\n

Instructions for running Coev with a command line:<\/h3>\n

Unzip Coev.zip in a directory of your choice.\u00a0Once all external software installed, to compile run:make<\/p>\n

Then to execute run:<\/p>\n

.\/coev [OPTIONS]\\r\\n  -method s:  Use method s (either ''bayes'' or ''ml''). Default is ''ml: maximum likelihood''.\r\n\u00a0-IT n:         Run for n iterations. Default is n=10000.\r\n-sfreq n:     Write every n''th iteration to file. Default is n=1000.\r\n\u00a0-pfreq n:     Print every n''th iteration on the screen. Default is n=1000.\r\n-burnin n:    Number of burn-in iterations. Default is n=0.\r\n-s v:          ??. Default is v=1.\r\n\u00a0-d v:          ??. Default is v=1.\r\n\u00a0-w1 v:        ??. Default is v=1.\r\n\u00a0-w2 v:        ??. Default is v=1.\r\n\u00a0-tree s:       The name of the file containing the input tree in Newick format. Default is s=treeInput.txt.\r\n\u00a0-align s:      The name of the file containing the sequence alignment in FASTA format. Default is s=alignment.txt.\r\n\u00a0-out s:        The name of the log file to write the results to. Default is s=output.log.\r\n\u00a0-cols n1 n2:  The columns in the alignment to use for the analysis. Default is n1=1 and n2=2.\r\n\u00a0-h:             Print this help screen and exit.\r\n<\/pre>\n
Reading output:<\/h3>\n
To read the output, the user can either look at the log file or at the output console. In the log file, each line corresponds to a profile.<\/p>\n
In the case of a Maximum Likelihood framework:<\/h4>\n
Each line can be truncated in 9 values:<\/p>\n
    \n
  • Value 1 corresponds to the profile set: it is a binary vector of length 16 where 1 corresponds\u00a0to a combination that belong to a profile and 0 when it doesn”t. The binary vector follows the following order AA,AC,AG,AT,CA,CC,CG,CT,GA,GC,GG,GT,TA,TC,TG,TT. For example when the value is\u00a01000010000000000 then it means that the profile is {AA,CC}.<\/li>\n
  • Value 2 corresponds\u00a0to the estimated value of w1 in the Null model: it\u00a0\u00a0is a double, ex:6.928542e-02<\/li>\n
  • Value 3 corresponds\u00a0to the estimated value of w2 in the\u00a0Null model:\u00a0it\u00a0\u00a0is a\u00a0double, ex:2.014202e-01<\/li>\n
  • Value 4 corresponds\u00a0to the estimated value of Log Likelihood in the Null model:\u00a0it\u00a0is a\u00a0doube,\u00a0ex:-7.405995e+01<\/li>\n
  • Value 5 corresponds\u00a0to the estimated value of s in the model:\u00a0it\u00a0is a\u00a0double, ex:2.735160e-02<\/li>\n
  • Value 6 corresponds\u00a0to the estimated value of d in the\u00a0Coev model:\u00a0it is a\u00a0double, ex:-4.970079e-01<\/li>\n
  • Value 7 corresponds\u00a0to the estimated value of r1 in the\u00a0Coev model:\u00a0it is a\u00a0double, ex:3.103843e-02<\/li>\n
  • Value 8 corresponds\u00a0to the estimated value of r2 in the\u00a0Coev model:\u00a0it is a\u00a0double, ex:2.922180e-01<\/li>\n
  • Value 9 corresponds\u00a0to the estimated value of\u00a0Log\u00a0Likelihood\u00a0in the\u00a0of the Coev model: it is a\u00a0double, ex:-6.814893e+01<\/li>\n<\/ul>\n
    In the case of a Bayesian framework:<\/h4>\n
    \n
    The log file is conceived to be read by tracer. The first line is a header. The other lines contains values where<\/div>\n<\/div>\n
      \n
    • Value\u00a01 corresponds to the iteration identifier: it is an interger<\/li>\n
    • Value 2 corresponds to the posterior\u00a0probabilityr: it is a double<\/li>\n
    • Value 3 corresponds to the\u00a0LogLikelihood value of Coev: it is a double<\/li>\n
    • Value 4 corresponds to the prior value:\u00a0 it is a double<\/li>\n
    • Value 5 corresponds to s:\u00a0it is a double<\/li>\n
    • Value 6 corresponds to d:\u00a0it is a double<\/li>\n
    • Value 7 corresponds to r1:\u00a0it is a double<\/li>\n
    • Value 8 corresponds to the r2:\u00a0it is a double<\/li>\n
    • The rest of the values correspond\u00a0to the \u00a0profile vector\u00a0\u00a0where 1 corresponds\u00a0to a combination that belong to a profile and 0 when it doesn”t.\u00a0The binary vector follows the following order AA,AC,AG,AT,CA,CC,CG,CT,GA,GC,GG,GT,TA,TC,TG,TT. For example when the value is\u00a01000010000000000 then it means that the profile is {AA,CC}.<\/li>\n<\/ul>\n
      Contacts:<\/h3>\n
      For questions, please send an email to wwwphylo [at]\u00a0unil.ch<\/p>\n
      Reference:<\/h3>\n
      When using this software, please cite:<\/p>\n
      L. Dib, D. Silvestro, N.Salamin. 2014. Evolutionary footprint of coevolving positions in genes. Bioinformatics. 30(9):1241-9.<\/address>\n<\/div>\n","protected":false},"excerpt":{"rendered":"
      Coev\u00a0is a software that predicts\u00a0coevolving\u00a0positions and their evolutionary profile based on the aligned sequences and a phylogenetic tree. Overview: The software is the implementation of the model described in Dib et al. 2013, where we present a new substitution model<\/p>\n","protected":false},"author":1124,"featured_media":0,"parent":132,"menu_order":1,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_seopress_robots_primary_cat":"","_seopress_titles_title":"","_seopress_titles_desc":"","_seopress_robots_index":"","footnotes":""},"class_list":["post-344","page","type-page","status-publish"],"_links":{"self":[{"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/pages\/344","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/users\/1124"}],"replies":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/comments?post=344"}],"version-history":[{"count":0,"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/pages\/344\/revisions"}],"up":[{"embeddable":true,"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/pages\/132"}],"wp:attachment":[{"href":"https:\/\/wp.unil.ch\/phylo\/wp-json\/wp\/v2\/media?parent=344"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}