Simulation of gene family histories
© Hernandez-Rosales et al; licensee BioMed Central Ltd. 2014
Published: 11 February 2014
The way gene families and genomes evolve can be understood in detail only when the location of gene duplication episodes in the tree of life can be deciphered. Since most genes belong to larger gene families, the analysis of the gene family histories thus plays an important role in the study of genome evolution. Empirically, one frequently observes that the tree that describes the evolution of species, the species tree, is inconsistent with the tree that is obtained from a group of genes of a gene family (the gene tree). Goodman et al. deduced that this inconsistency might be the result of mistaking paralogs for orthologs. Orthologous genes refer to copies of genes that reveal the phylogeny of species, while paralogous genes have been created by duplication events. Phylogeny reconstruction can help to understand how gene families evolved and to identify the chronology of duplications within a gene family of a single species. Several software tools, including GeneTree, DupTree, NOTUNG, and AUGIST have been developed for this task. There is, however, lack of both test data and evaluation procedures to test, compare, and benchmark their performance and results.
We present here a simulation environment designed to generate large gene families with complex duplication histories on which reconstruction algorithms can be tested and software tools can be benchmarked.
The simulation of gene family histories starts with the generation of species trees. Within these rooted bifurcating trees the nodes represent species and edges their relation. Specifically, internal nodes represent ancient species whereas leaf nodes represent extant species. Given a number of species N, we generate a random tree T under the Age Model described in Keller-Schmidt et al. This model starts with a rooted tree with two leaves. In an iterative process one of the leaves is selected and two new leaves are attached to it until the tree has N leaves. This model makes use of the idea that the longer a leaf has not been involved in a speciation, the less likely it will be in the future.
Finally, the algorithm will generate one gene tree for each species, i.e. the pruned reconciled tree containing only genes of a certain species. Furthermore, for each gene family the orthology and homology matrices are computed. To generate the orthology matrix, we say that two genes are orthologous if their lowest common ancestor (LCA) in the reconciled tree represents a speciation event. To generate the homology matrix, a gene a from species i is homologous to gene b from species j if for every gene c from species i and every gene d from species j the LCA(a, b) ≤ LCA(c, b) and LCA(a, b) ≤ LCA(a, d).
We propose an algorithm that simulates gene family histories akin to real data. This will allow reconstruction algorithms to measure their accuracy and performance. Given a certain reconstruction method one might ask if the orthology matrix could be deduced from the inferred reconciled tree or if the homology relation between the genes was predicted correctly. Furthermore it could be analysed if the method was able to infer the gene duplications and losses. A method that is able to detect large scale duplications will then identify the cluster and genome duplications generated by our algorithm.
- Goodman M, Czelusniak J, Moore GW, Romero-Herrera AE, Matsuda G: Fitting the Gene Lineage into its Species Lineage, a Parsimony Strategy Illustrated by Cladograms Constructed from Globin Sequences. Systematic Zoology. 1979, 28: 132-163. 10.2307/2412519.View ArticleGoogle Scholar
- Page RD: GeneTree: comparing gene and species phylogenies using reconciled trees. Bioinformatics. 1998, 14: 819-820. 10.1093/bioinformatics/14.9.819.View ArticlePubMedGoogle Scholar
- Wehe A, Bansal MS, Burleigh JG, Eulenstein O: DupTree: a program for large-scale phylogenetic analyses using gene tree parsimony. Bioinformatics. 2008, 24: 1540-1541. 10.1093/bioinformatics/btn230.View ArticlePubMedGoogle Scholar
- Vernot B, Stolzer M, Goldman A, Durand D: Reconciliation with non-binary species trees. Comput Syst Bioinformatics Conf. 2007, 6: 441-452.View ArticlePubMedGoogle Scholar
- Oliver JC: AUGIST: inferring species trees while accommodating gene tree uncertainty. Bioinformatics. 2008, 24: 2932-2933. 10.1093/bioinformatics/btn556.PubMed CentralView ArticlePubMedGoogle Scholar
- Keller-Schmidt S, Tugrul M, Eguiluz VM, Hernandez-Garcia E, Klemm K: An Age Dependent Branching Model for Macroevolution. (Submitted on 15 Dec 2010). arXiv:1012.3298v1Google Scholar
- Ohno S: Gene duplication and the uniqueness of vertebrate genomes circa 1970–1999. Seminars in Cell and Developmental Biology. 1999, 10: 517-522. 10.1006/scdb.1999.0332.View ArticlePubMedGoogle Scholar
- Sonja Prohaska, Claudia Fried, Christoph Flamm, Günter P Wagner, Peter F Stadler: Surveying Phylogenetic Footprints in Large Gene Clusters: Applications to Hox Cluster Duplications. Mol Phylogenet Evol. 2004, 31: 581-604. 10.1016/j.ympev.2003.08.009.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.