Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

This article is part of the supplement: Highlights from the Ninth International Society for Computational Biology (ISCB) Student Council Symposium 2013

Open Access Meeting abstract

Simulation of gene family histories

Maribel Hernandez-Rosales12*, Nicolas Wieseke2, Marc Hellmuth2 and Peter F Stadler12

Author Affiliations

1 Max-Planck-Institute for Mathematics in the Sciences, Inselstr. 22, D-04103 Leipzig, Germany

2 Department of Computer Science, Univ. Leipzig, Härtelstr. 16-18, D-04107 Leipzig, Germany

For all author emails, please log on.

BMC Bioinformatics 2014, 15(Suppl 3):A8  doi:10.1186/1471-2105-15-S3-A8

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2105/15/S3/A8


Published:11 February 2014

© 2014 Hernandez-Rosales et al; licensee 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.

Background

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.

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.

The user will introduce n number of genes (gene families), which will be placed at the root of the generated species tree T. T will then be traversed in a depth first order. For each visited edge a number of events is sampled from a stochastic Poisson Process <a onClick="popup('http://www.biomedcentral.com/1471-2105/15/S3/A8/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/15/S3/A8/mathml/M1">View MathML</a> where <a onClick="popup('http://www.biomedcentral.com/1471-2105/15/S3/A8/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2105/15/S3/A8/mathml/M2">View MathML</a> is the probability of the event to happen and l the branch length. The process may generate none, one or a series of these events: one gene gets duplicated (gene duplication), a group of genes gets duplicated (cluster duplication), the whole group of genes gets duplicated (genome duplication) and one gene of the species gets lost (gene loss). After each gene duplication, one of the copies will be lost with a user defined probability θ, based on the fact that when there is a gene duplication, one of the copies might be lost or become nonfunctional. In the case of a cluster or genome duplication, we apply this probability to every gene in the group, since it is known that in the wake of multiple gene duplications and in particular for genome duplications we have to expect that many duplicated genes are rapidly lost again through the formation of pseudogenes. A small example of a gene family history generated by our simulation is shown in Fig. 1. We also show the gene tree generated from the gene family history embedded in the species tree. Each leaf node represents a gene and each internal node represents an event (speciation or duplication). This tree is typically depicted as the reconciled tree as in Fig. 2.

thumbnailFigure 1. A one-gene family history: from a node parent to a node child, there could be duplications and losses of genes.

thumbnailFigure 2. The reconciled tree: the gene tree embedded in the species tree. Each internal node represents an event, either a speciation or a gene duplication.

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).

Conclusions

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.

References

  1. 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. Publisher Full Text OpenURL

  2. Page RD: GeneTree: comparing gene and species phylogenies using reconciled trees.

    Bioinformatics 1998, 14:819-820. PubMed Abstract | Publisher Full Text OpenURL

  3. 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. PubMed Abstract | Publisher Full Text OpenURL

  4. Vernot B, Stolzer M, Goldman A, Durand D: Reconciliation with non-binary species trees.

    Comput Syst Bioinformatics Conf 2007, 6:441-452. PubMed Abstract | Publisher Full Text OpenURL

  5. Oliver JC: AUGIST: inferring species trees while accommodating gene tree uncertainty.

    Bioinformatics 2008, 24:2932-2933. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. 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.3298v1

  7. Ohno S: Gene duplication and the uniqueness of vertebrate genomes circa 1970–1999.

    Seminars in Cell and Developmental Biology 1999, 10:517-522. PubMed Abstract | Publisher Full Text OpenURL

  8. 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. PubMed Abstract | Publisher Full Text OpenURL