PhyBin: tree binning by topology

Just submitted our first paper to PeerJ - the awesome new, open access journal aimed at tipping the entire publishing establishment on its head.  I'm looking forward to a smooth review process -- hopefully as sleek and helpful as the submission process. <UPDATE: our paper was accepted.  I will post a link to the paper once it's online!>

Our manuscript presents PhyBin, a computer program aimed at binning precomputed sets of trees in Newick format, a file format produced by the majority of tree building software.  As we assert in the manuscript, PhyBin is a utility rather than a complete solution; it can serve as a component in many genomics pipelines, and provides a useful addition to the landscape of tools for dissecting and visualizing large numbers of trees.  After the user applies their chosen ortholog prediction and tree-building algorithms, PhyBin offers a quick way to visualize and browse the different evolutionary histories, either binned by topology and sorted by bin size, or in the form of a full hierarchical clustering based on Robinson-Foulds (RF) distance: i.e. a tree of trees.

In the manuscript, we explore to functionalities in PhyBin: 1) the ability to bin trees with identical topologies and 2) the ability to cluster similar trees by RF distance.  Lots of folks interested in the "landscape" of topologies produced by orthologous genes across a genome use RF distance as a measure of topological similarity.  What is RF distance?  It is essentially the number of different steps you'd have to take to create one tree out of another -- it's the edit distance between two topologies.  So, according to the original Robinson-Foulds publication, for example, the trees below (trees 1 and 2) are edit distance 2 apart because in order to convert one to the other, you must collapse a node and then reform it.

PhyBin does some pretty neat pre-processing of trees to facilitate comparissons. For example, you can set a branch length threshold to collapse branches that are essentially noise in your dataset (say, from very closely related taxa).  It also checks your dataset for number of taxa and is quite robust to file formatting.  Then what PhyBin does is calculate the edit distance for a large group of trees (a distance matrix) and then also displays these distances as a tree of trees - as a dendrogram that links each tree in the dataset to each other based on the edit distances between them (how you'd get from one to another).  In our manuscript, we used a set of orthologs generated from 10 published Wolbachia genomes.  Here's what the dendrogram looks like for those 508 trees without (A) and with (B) clustering by RF distance. In the figure, you can see that many of the trees in the Wolbachia ortholog set are similar, they cluster into 9 large clusters, many of which support the monophyly of Wolbachia supergrops but others (in fact, a good fraction!) which do not.  


Want to download it?  Check it out here





Comments

Post a Comment

Popular posts from this blog

The Wolbachia "holy grail"

Image
The once-obscure alpha-proteobacterium Wolbachia pipientis, was catapulted into medical relevance by the discovery that it inhibits the replication of RNA viruses .  However, it was first discovered for its odd quirk of manipulating insect reproduction . That's right, Wolbachia were called "reproductive parasites" and for a long time, no fitness benefit was attributed to them.  These intracellular masters do things like kill male offspring, feminize male offspring, induce parthenogenesis, and the most common phenotype, sperm-egg incompatibility (also known by its more complicated name: cytoplasmic incompatibility or CI) (Figure 1).   Figure 1. Modified from Werren et al., 2008 .  Wolbachia cause four distinct reproductive phenotypes in a range of arthropod orders (top).  The mechanisms behind CI - meaning how does  Wolbachia  induce it -  have been a mystery for some time.   When you make a cross between infected males and uninfected...
Read more

What's in a name? Would a Wolbachia by any other name, block RNA viruses as sweetly?

Image
I've ranted about bacterial nomenclature in the past (see h ere ), but this post is specifically about Wolbachia pipientis, and how we should name strain-level variants within a particular host-associated clade (such as w Mel, for example).  Currently, there seem to be multiple ways to write these names and to present them in databases or the literature.  For example: For variants within w Mel, we find examples of the following types of strain names: w MelPop-<insert specific substrain here after the dash> The stuff that comes after the strain name ( w MelPop) sometimes refers to the host insect from which the strain was isolated, or where it currently resides (such as in w MelPop-PYGP) or sometimes refers to how the strain was altered, through sequential passage (such as in w MelPop-CLA).  This variability in the naming schemas used in Wolbachia could lead to confusion down the road. We are currently adding strain level information to some of the stocks th...
Read more

ASM's MRA Journal Supports Undergraduate Research

Image
Disclaimer: Dr. Irene Newton is EiC of ASM's Microbiology Resource Announcements and Catherine Putonti is a Senior Editor. Undergraduate research changed my life. Truly. Along with many thousands of other undergraduate students, I started my career really thinking I was going to be an MD. That career path appealed for many reasons - my immigrant family understood it and approved, it was prestigious, and it meant my life was planned out for the next decade. But one semester of research with my undergraduate mentor was enough to get me hooked. I was privileged to have the opportunity to do one-on-one research at a small liberal arts institution but not all of our students have that opportunity. I now teach at one of the "Big 10" - with over 45,000 undergrads on campus there is no way we can provide all of them the research experience I had. But we can and are trying to integrate true scientific research into our undergraduate classes. Some good examples of this are the FI...
Read more