Showing posts with label software. Show all posts
Showing posts with label software. Show all posts

Tuesday, 31 March 2015

ClonalFrameML: accounting for recombination in bacterial phylogenies

Horizontal gene transfer in bacteria, mediated by transformation, transduction or conjugation, can result in gain, loss and replacement of genes. The replacement of horizontally transferred genes or gene fragments in a process known as homologous recombination has far-reaching effects on bacterial phylogenetics - the study of relatedness between bacteria. A new method published by Xavier Didelot and me last month in PLoS Computational Biology corrects for these distorting effects of homologous recombination on bacterial phylogenies.

Two forms of phylogenetic distortion are caused by recombination. The first affects the shape of the tree topology. Although this is a potentially serious difficulty, Jessica Hedge and I recently showed that phylogenies estimated from whole bacterial genomes are surprisingly robust to this problem. The second affects the lengths of the branches. When genetic material is replaced by a homologous but distantly related sequence, it gives the appearance of a cluster of substitutions in the genome, and this can exaggerate branch lengths. ClonalFrameML detects these clusters of substitutions, identifies them as recombination events, and corrects the branch lengths of the tree.

Correcting for recombination is important in a variety of settings. In transmission studies, recent transmission between patients can be detected by comparing the genomes of the infecting bacteria. As we show in the paper, ClonalFrameML improves detection of transmission events by accounting for the tendency of recombination to elevate the evolutionary distance between genomes. We also report the discovery of a remarkably large chromosomal replacement event spanning 310 kilobases that may have led to the evolution of the ST582 strain of Staphylococcus aureus, underlining the importance of recombination over short and long timescales.

ClonalFrameML is a much faster implementation of the popular ClonalFrame method by Xavier and Daniel Falush. It is based on the same underlying assumptions and the same explicit evolutionary model, so it provides interpretable estimates of rates of recombination, the length of DNA imported by recombination, and the relative impact of recombination versus mutation. However, it can now analyse thousands of whole bacterial genomes in a matter of hours, representing a substantial improvement over the earlier method.

Thursday, 16 August 2012

gammaMap available for download

The software gammaMap - which implements the analyses developed in Wilson, Hernandez, Andolfatto and Przeworski (2011) PLoS Genetics 7: e1002395 - is available for download. It is provided as part of a flexible program called GCAT (general computational analysis tool) which is designed to rapidly facilitate novel variations on the standard analyses. GCAT has its own google code page, http://code.google.com/p/gcat-project. GCAT resembles BEAST and BUGs in that a statistical model is specified (using XML) and parameters are then estimated using MCMC or maximum likelihood. Future extensions to GCAT are planned that implement new fast approximations to gammaMap and omegaMap, and parallel processing, allowing the analyses to be scaled more readily to whole genomes.