LOSITAN: A workbench to detect molecular adaptation based on a F st -outlier method
© Antao et al; licensee BioMed Central Ltd. 2008
Received: 15 February 2008
Accepted: 28 July 2008
Published: 28 July 2008
Testing for selection is becoming one of the most important steps in the analysis of multilocus population genetics data sets. Existing applications are difficult to use, leaving many non-trivial, error-prone tasks to the user.
Here we present LOSITAN, a selection detection workbench based on a well evaluated F st -outlier detection method. LOSITAN greatly facilitates correct approximation of model parameters (e.g., genome-wide average, neutral F st ), provides data import and export functions, iterative contour smoothing and generation of graphics in a easy to use graphical user interface. LOSITAN is able to use modern multi-core processor architectures by locally parallelizing fdist, reducing computation time by half in current dual core machines and with almost linear performance gains in machines with more cores.
LOSITAN makes selection detection feasible to a much wider range of users, even for large population genomic datasets, by both providing an easy to use interface and essential functionality to complete the whole selection detection process.
Understanding the contribution of selection and molecular adaptation in shaping genome wide variation is among the most exciting and widely researched problems with many applications ranging from human health to conservation of endangered species. Among the many selection detection strategies , F st outlier approaches are becoming widely used [2, 3] because they are important not only for studying the genetic basis of adaptation but also for eliminating non-neutral outlier loci from data sets before computing most population genetic parameters (e.g., F st , N m , N e ), that require neutral loci . This is particularly important in a time where production of data sets with information from hundreds of loci is becoming fairly common.
One such F st method is described in [2, 5] (but see also  and ) and is implemented in the fdist program and can be used for any codominant genetic molecular markers including microsatellites, Single Nucleotide Polymorphisms (SNPs) and allozymes. This method evaluates the relationship between F st and H e (expected heterozygosity) in an island model , describing the expected distribution of Wright's inbreeding coefficient F st vs. H e under an island model of migration with neutral markers. This distribution is used to identify outlier loci that have excessively high or low F st compared to neutral expectations. Such outlier loci are candidates for being subject to selection.
Using fdist can be a challenging task for those not familiarized with command-line applications and requires a specific data format not used by other applications . Furthermore, several independent runs are usually needed to tune parameters (e.g., determine the appropriate average F st ) before a final execution is made in a process that is prone to human introduced mistakes. Fdist, not being one of the most computationally intensive programs available, can still take up to one hour for a single run (especially if smooth contours for confidence intervals are required), and, in most cases, multiple runs are needed for parameter tuning. Large population genomic datasets can take even longer. In this context, fdist requires experienced computer users, and its usage is error prone (e.g., by incorrectly converting data files or not approximating average F st appropriately).
We designed LOSITAN (LOoking for Selection In a TANgled dataset), a selection detection workbench constructed around fdist. LOSITAN is a Java Web Start application coded mostly in Jython with a small part in Java, allowing direct execution from the web. LOSITAN provides the following features:
2. Data import in Genepop  format.
3. Generation of graphics in several formats (PNG, SVG and PDF).
4. Data export in a format suitable for import into statistical packages like R  or commonly used spreadsheet software.
In case the user desires to further analyze the data or have total flexibility in generating graphics, LOSITAN makes available both the confidence intervals computed and the F st s and heterozygosities for each locus.
A simple R script is supplied in order to facilitate loading the data into R. Loading in spreadsheet software is done simply by importing as a tab delimited file.
5. Choice of which populations and/or loci are studied.
6. Approximating mean neutral F st (in the real dataset) by removing potential selected loci.
The initial mean dataset F st is often not neutral in the sense that (initially unknown) selected loci are often included in the computation. LOSITAN can optionally be run once to determine a first candidate subset of selected loci in order to remove them from the computation of the neutral F st . This value will be, in most cases, a better approximation of the neutral F st . The procedure works as follows: LOSITAN is run a first time, using all loci to estimate the mean neutral F st . After the first run, all loci that are outside the desired confidence intervals (e.g. 99% CIs) are removed and the mean neutral F st is computed again using only putative neutral loci that were not removed. A second and final run of LOSITAN, using all loci, is then conducted using the last computed mean. This procedure lowers the bias on the estimation of the mean neutral F st by removing the most extreme loci from the estimation. Naturally all loci will be present in the last run will have their estimated selection status reported.
7. Approximating average simulated F st to the average value found in the real dataset even when the experimental conditions are far from the ones where the theoretical formula holds (e.g. low number of demes or the usage of the stepwise mutation model, common in microsatellite markers).
To be able to (optionally) approximate the average F st in conditions far from the theoretical optimum, LOSITAN starts by running fdist for 10,000 realizations using the theoretical value, calculating the average simulated F st , if the value is too far from the real average F st , LOSITAN uses a bisection approximation algorithm running 10,000 realizations for every tentative bisection point. The algorithm works by iteratively slicing the interval of possible F st values (i.e., between 0 and 1) in half at each iteration and choosing the mean of the bounds on each iteration (with the exception of the first iteration where one of the extremes is chosen). An example is provided to make the approach clearer:
In a certain demographic scenario we want to simulate a neutral F st of 0.08. The algorithm starts by trying 0.08. If the result is higher than desired then 0.0 will be tried (creating an absolute lower bound limit), after that 0.04 (0.0 + 0.08)/2 will be tried, if the result is too low, 0.06 will be used next (i.e. (0.04 + 0.08)/2), the process repeats until the error margin is acceptable.
In practical terms the method was able to converge to the desired value in all cases tested (a completely trivial bisection approach is not possible as the method for computing F st is stochastic and results might vary for the same input conditions).
8. Iterative smoothing of confidence interval contours.
Contour smoothing is achieved by running fdist an extra 5,000 realizations. The user can request smoothing an unlimited number of times until the result is deemed satisfactory.
9. Ability to use multiple CPU cores and processors when running fdist.
To be able to use multiple cores, LOSITAN divides the number of desired simulation repeats among all available cores (although the application detects the number of existing cores, the user is able to change the number of simultaneous concurrent processes), this is possible because fdist simulation runs are independent, thus making parallelization a simple task. Tests show a near linear relationship between the number of cores used and performance gains, an existing 5–10% penalty is due mainly to joining the partial results together. LOSITAN, although being directly executable from the web is a client-side application and all computational intensive operations occur on the user computer and not on the server.
10. Automatic and transparent download of the latest version of fdist.
We maintain the latest version of the fdist application on the server, which is downloaded transparently by the client application whenever there is a new version. At the time of this writing the supported version is fdist2.
The interface includes tips for all the less obvious parameters and enforces constraints for all the user inputs which the system can infer are not correct.
Results and discussion
In a beta test release to users the feedback was generally very positive stressing essentially that the application is easy to use, allows to easily input and output data and deal with non-trivial parameter determination like calculating neutral F st . Most importantly it made users aware of issues in data analysis that they were not aware of. For example, users were not aware of how to estimate the genome wide average neutral F st from their empirical data set by removing one or a few strong outlier loci, and the recomputing the average F st . Although LOSITAN helps avoid many pitfalls involved with using F st -outlier approaches in general, it is not able to solve fundamental issues regarding these approaches, for instance the non-linear behavior of when F st approaches zero can make it difficult to detect low F st -outliers especially when selection is not strong. As such an easy to use application should not be seen by users as a excuse to avoid critical reasoning around the the whole selection detection process. Feedback from users also allows to chart possible future work, like supporting dominant markers or supporting other selection detection approaches like .
Our solution to use all the available computing power on new multi-core hardware is an example of an "embarrassingly simple parallel" computation approach. We contend that having a simple approach is a good principle: The point in this application is to make all computational power available to the users and not to develop new concurrent algorithms. A simple, highly efficient, elegant and less bug-prone approach is what responds to the users needs, as the objective of this work is not to develop new algorithms, but to use them.
LOSITAN is built along the principles exposed in , namely that intuitiveness and user empowerment should be fundamental guidelines for software construction targeting biologists. This is done, not only by supplying an easy to use web interface for an, otherwise, hard to use application, but also allowing the use of widely utilized population genetic data formats, automating the tuning of nuisance parameters and lowering the computational costs on modern hardware. In addition, strong emphasis is put on trying to avoid errors on the usage of the software either by both enforcing constraints and giving suggestions on less obvious features. This will lower the barriers to usage of the underlying application, allowing for a wider user base which will be able to concentrate more on the biological problems and less on unnecessary application complexity.
We are in the dawn of the era of multi-core computing. The vast majority of existing software cannot make use of the extra computational power made available on new machines. Our approach, based on partitioning a computational intensive task into smaller ones, can be used to leverage the extra computational power even without changing existing code on applications which can be broken into smaller independent running units. This partitioning approach can be performed in some cases by users on existing software or by programmers in new applications that take advantage of multiple cores. With the current trend of supplying many more cores with new computers, strategies like the one presented here will be mandatory in order to take full advantage of all the existing processing power. LOSITAN is one of the first of many applications to explore the multi-core programming paradigm.
Future planned developments will include addition of other F-outlier methods and simulation facilities for explore the effects of different demographic scenarios on F st variance and the detection of outliers. All the code to handle GenePop and fdist file formats and applications was also donated to the Biopython project and is publicly available starting from version 1.44.
Availability and requirements
Project name LOSITAN
Operating systems Platform independent
Programming language Java and Jython
Other requirements Browser with JavaWebStart to run over the internet (software can be run locally).
Windows: At least Windows 2000 and Java 1.6.
Mac OS X: 10.4 (Tiger) and Java 1.5 (Most current 10.4 installations will require a freely available Java update).
Linux: Java 1.6 and the free GNU C compiler.
License GNU GPL
Any restrictions to use by non-academics None
This work was partially supported by the Bill & Melinda Gates Foundation (grant #39777).
TA was supported by research grant SFRH/BD/30834/2006, RJL by SFRH/BPD/14953/2004 and AB-P by SFRH/BPD/17822/2004 and this work was supported by POCI/CVT/567558/2004 all from Fundacao para a Ciencia e Tecnologia (FCT), Portugal. GL was supported by the Luso-American Foundation, UP, CIBIO and research grant PTDC/BIA-BDE/65625/2006 from FCT.
- Nielsen R: Molecular signatures of natural selection. Annual Reviews in Genetics 2005, 39: 197–218. 10.1146/annurev.genet.39.073003.112420View ArticleGoogle Scholar
- Beaumont MA: Adaptation and speciation: what can Fst tell us? Trends Ecol Evol 2005, 20(8):435–440. 10.1016/j.tree.2005.05.017View ArticlePubMedGoogle Scholar
- Vitalis R, Dawson K, Boursot P: Interpretation of variation across marker loci as evidence of selection. Genetics 2001, 158(4):1811–1823.PubMed CentralPubMedGoogle Scholar
- Luikart G, England PR, Tallmon D, Jordan S, Taberlet P: The power and promise of population genomics: from genotyping to genome typing. Nat Rev Genet 2003, 4(12):981–994. 10.1038/nrg1226View ArticlePubMedGoogle Scholar
- Beaumont MA, Nichols RA: Evaluating loci for use in the genetic analysis of population structure. Proceedings of the Royal Society B 1996, 363: 1619–1626. 10.1098/rspb.1996.0237View ArticleGoogle Scholar
- Cavalli-Sforza LL: Population Structure and Human Evolution. Proc R Soc Lond B Biol Sci 1966, 164: 362–379.View ArticlePubMedGoogle Scholar
- Lewontin RC, Krakauer J: Letters to the editors: Testing the heterogeneity of F values. Genetics 1975, 80(2):397–398.PubMed CentralPubMedGoogle Scholar
- Wright S: Evolution in Mendelian Populations. Genetics 1931, 16(2):97–159.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Heckel G: Computer programs for population genetics data analysis: a survival guide. Nat Rev Genet 2006, 7(10):745–758. 10.1038/nrg1904View ArticlePubMedGoogle Scholar
- Raymond M, Rousset F: GENEPOP: population genetics software for exact tests and ecumenicism. Journal of Heredity 1995, 86: 248–249.Google Scholar
- R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2007.Google Scholar
- Kumar S, Dudley J: Bioinformatics software for biologists in the genomics era. Bioinformatics 2007, 23(14):1713–1717. 10.1093/bioinformatics/btm239View ArticlePubMedGoogle 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.