{"id":203,"date":"2014-01-23T11:49:41","date_gmt":"2014-01-23T11:49:41","guid":{"rendered":"http:\/\/pop-gen.eu\/wordpress\/?page_id=203"},"modified":"2024-07-01T20:06:38","modified_gmt":"2024-07-01T20:06:38","slug":"comus-coalescent-of-multiple-species","status":"publish","type":"page","link":"https:\/\/pop-gen.eu\/wordpress\/software\/comus-coalescent-of-multiple-species","title":{"rendered":"CoMuS: Coalescent of Multiple Species and CoMuStats"},"content":{"rendered":"<h2><a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2012\/04\/image_t6.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\" wp-image-34 alignright\" src=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2012\/04\/image_t6-201x300.jpg\" alt=\"image_t6.jpg\" width=\"188\" height=\"281\" srcset=\"https:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2012\/04\/image_t6-201x300.jpg 201w, https:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2012\/04\/image_t6.jpg 590w\" sizes=\"auto, (max-width: 188px) 100vw, 188px\" \/><\/a>Introduction<\/h2>\n<p>Here, you can download the source code of CoMuS, CoMuStats, as well as R scripts that can be used for different analyses. The R scripts implement specific examples that have been used in CoMuS&#8217; manuscript or, in general, I consider them interesting. You can download them, modify them according to your needs.<\/p>\n<p>You can download the most recent version from the github:<br \/>\n<a href=\"https:\/\/github.com\/idaios\/comus\">https:\/\/github.com\/idaios\/comus<\/a><\/p>\n<p>or type in terminal: git clone https:\/\/github.com\/idaios\/comus.git<\/p>\n<h3><strong>Instructions:<\/strong><\/h3>\n<ol>\n<li>download the code from the previous link<\/li>\n<li>tar xvfz comus.tar.gz<\/li>\n<li>cd comus<\/li>\n<li>to compile: make -f Makefile.gcc (you may need to remove the *.o files)<br \/>\ni.e. rm *.o and then make -f Makefile.gcc<\/li>\n<li>now you should have the comus executable<\/li>\n<li>IMPORTANT: there is a pre-compiled executable in the comus.tar.gz. However, sometimes, depending on the system you will not be able to execute it.Instead, you will see the error message:<em>comus: \/lib64\/libc.so.6: version `GLIBC_2.14&#8242; not found (required by comus). <\/em>Please first clean: <strong>make clean -f Makefile.gcc<\/strong><br \/>\nand then recompile the code: <strong>make -f Makefile.gcc\u00a0<\/strong><\/li>\n<\/ol>\n<h3>Scripts and demonstration<\/h3>\n<p>In CoMuS manuscript we have used several scripts, simulations and inference examples to demonstrate CoMuS usage. These scripts are provided here either for demonstration or testing purposes.<\/p>\n<ul>\n<li><strong><strong><a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/ancestral_sampling_code_results_noOutput.tar.gz\">ancestral sampling<\/a><br \/>\n<\/strong><\/strong><\/p>\n<p align=\"justify\"><span style=\"font-family: Arial,sans-serif;\">CoMuS allows the simultaneous simulation of both modern and ancestral samples. To facilitate the simultaneous coalescent simulation of modern and ancestral samples, CoMuS implementation starts at the present-day with the whole dataset (modern and ancestral). However, all events that involve ancestral samples or their population are forbidden until sampling (time proceeds backwards). After sampling (backwards in time), evolutionary processes (e.g. recombination, mutations, coalescent, migration etc) take place as usually according to the model parameters.<br \/>\n<\/span><\/p>\n<p><span style=\"font-family: Arial,sans-serif;\">Here, we demonstrate the usage of CoMuS to infer potential ancestral gene flow between an extant population and an extinct sample (fossil). The simulation scenario is as follows: we assume a sample of 10 sequences from species A sampled at present, and a sample of 10 sequences from an extinct species B sampled at time 0.2 (phylogenetic time units). The time of the MRCA has been set to 0.5 (phylogenetic time units). We assume no gene flow after speciation between the species A and B. We illustrate the dendrogram for this scenario in Figure 1 (below)<\/span><span style=\"font-family: Arial,sans-serif;\">. Assuming that the above scenario represents the true evolutionary history for extant species A and extinct species B, our goal is to infer: (i) whether gene flow between A and B is absent or present, and (ii<\/span><span style=\"color: #0000ff;\"><span style=\"font-family: Arial,sans-serif;\">)<\/span><\/span><span style=\"font-family: Arial,sans-serif;\"> the time of sampling for species B. The time of the MRCA (= 0.5 phylogenetic units) as well as \u03b8 (= 100) value is assumed to be known. The length of the simulated region is 1kb and we assumed a mutation model with equal mutation rates between each pair of bases.<\/span><\/p>\n<p>Download the script <a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/ancestral_sampling_code_results_noOutput.tar.gz\">here<\/a><br \/>\n<a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/tree_ancestral.png\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-medium wp-image-398\" src=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/tree_ancestral-200x300.png\" alt=\"tree_ancestral\" width=\"200\" height=\"300\" srcset=\"https:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/tree_ancestral-200x300.png 200w, https:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/tree_ancestral-683x1024.png 683w, https:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/tree_ancestral.png 1600w\" sizes=\"auto, (max-width: 200px) 100vw, 200px\" \/><\/a><\/li>\n<\/ul>\n<p>Figure 1: a coalescent tree example comprises both present-day and ancestral sampling<\/p>\n<p><strong>Testing species delimitation software<\/strong><\/p>\n<ul>\n<li><strong>testing species delimitation with gene-flow<\/strong><\/li>\n<\/ul>\n<p>An apparent usage of CoMuS is to test species delimitation software. More specifically various parameters such as<\/p>\n<ul>\n<li>changes of population size<\/li>\n<li>population subdivision<\/li>\n<li>ghost populations<\/li>\n<li>etc<\/li>\n<\/ul>\n<p>can be examined and the performance of species delimitation software can be assessed. Here, we test the <a href=\"http:\/\/species.h-its.org\/\">PTP<\/a> software developed by Jiajie Zhang et al. (including myself) in the group of Alexandros Stamatakis in Heidelberg. The manuscript is available from <a href=\"http:\/\/bioinformatics.oxfordjournals.org\/content\/29\/22\/2869\">here<\/a>.<\/p>\n<p>The scripts used for this demonstration can be found <a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/subpopulations-givenTree.tar.gz\">here<\/a>. (NOTE: to run the full set of commands in the scripts you need to install raxml. see<a href=\"https:\/\/github.com\/stamatak\/standard-RAxML\"> raxml-github<\/a>)<\/p>\n<ul>\n<li>Testing species delimitation with various birth rates<\/li>\n<\/ul>\n<p>Scripts can be found <a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/birth_rate.tar.gz\">here<\/a>. The ideas are similar as those presented above.<\/p>\n<p><strong>Inferring parameters values<\/strong><\/p>\n<p>CoMuS can be used to infer parameter using the ABC framework. We have used two scenarios: (i) 2 species, each of 10 sequences sampled, inference of the birth rate <em>b<\/em> (ii) 10 species, each of 10 sequences, inference of the birth rate <em>b <\/em>and the time of the most recent common ancestor, i.e. the time that the two populations find common ancestor.<\/p>\n<ul>\n<li>Scripts for scenario (i) can be found <a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/figure_inference_2x10.tar.gz\">here<\/a><\/li>\n<li>Scripts for scenario (ii) can be found <a href=\"http:\/\/pop-gen.eu\/wordpress\/wp-content\/uploads\/2015\/06\/figure_inference_10x10.tar.gz\">here<\/a><\/li>\n<\/ul>\n<h3><strong>Examples and Manual<\/strong><\/h3>\n<ul>\n<li>Inside the comus directory you will find a directory called &#8216;manual&#8217; which contains the manual, i.e. manual.pdf<\/li>\n<li>Also there is a directory called examples. There is a run.sh file that contains commands as well as useful notes that explain most of the results. Please consult them first.<\/li>\n<\/ul>\n<h2>CoMuStats<\/h2>\n<p>This software can be used to calculate summary statistics from single or multi-FASTA alignment files. Multi-alignment fasta files should be separated by \/\/.<\/p>\n<p>For example:<\/p>\n<blockquote><p>&gt;seq1<br \/>\nACGTG<br \/>\n&gt;seq2<br \/>\nACGGG<br \/>\n\/\/<br \/>\n&gt;seq1<br \/>\nACCTC<br \/>\n&gt;seq2<br \/>\nACCCC<\/p><\/blockquote>\n<h3><strong>Bugs and Previous versions<\/strong><\/h3>\n<ol>\n<li>CoMuS v 1.0 : -oFormat unrecognized<\/li>\n<li>CoMuS v2.0 March 2016 (CoMuStats was not running properly the sliding window code).<\/li>\n<\/ol>\n<h3><strong>Log Report<\/strong><\/h3>\n<ol>\n<li>CoMuStats 1.0.1 accepts an outgroup. Run CoMuStats without arguments for further details.<\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Introduction Here, you can download the source code of CoMuS, CoMuStats, as well as R scripts that can be used for different analyses. The R scripts implement specific examples that have been used in CoMuS&#8217; manuscript or, in general, I consider them interesting. You can download them, modify them according to your needs. You can &hellip; <\/p>\n<p class=\"link-more\"><a href=\"https:\/\/pop-gen.eu\/wordpress\/software\/comus-coalescent-of-multiple-species\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;CoMuS: Coalescent of Multiple Species and CoMuStats&#8221;<\/span><\/a><\/p>\n","protected":false},"author":4,"featured_media":0,"parent":21,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"class_list":["post-203","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/pages\/203","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/users\/4"}],"replies":[{"embeddable":true,"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/comments?post=203"}],"version-history":[{"count":43,"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/pages\/203\/revisions"}],"predecessor-version":[{"id":794,"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/pages\/203\/revisions\/794"}],"up":[{"embeddable":true,"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/pages\/21"}],"wp:attachment":[{"href":"https:\/\/pop-gen.eu\/wordpress\/wp-json\/wp\/v2\/media?parent=203"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}