Toxo: a library for calculating penetrance tables of high-order epistasis models.

BACKGROUND: Epistasis is defined as the interaction between different genes when expressing a specific phenotype. The most common way to characterize an epistatic relationship is using a penetrance table, which contains the probability of expressing the phenotype under study given a particular allele combination. Available simulators can only create penetrance tables for well-known epistasis models involving a small number of genes and under a large number of limitations. RESULTS: Toxo is a MATLAB library designed to calculate penetrance tables of epistasis models of any interaction order which resemble real data more closely. The user specifies the desired heritability (or prevalence) and the program maximizes the table's prevalence (or heritability) according to the input epistatic model boundaries. CONCLUSIONS: Toxo extends the capabilities of existing simulators that define epistasis using penetrance tables. These tables can be directly used as input for software simulators such as GAMETES so that they are able to generate data samples with larger interactions and more realistic prevalences/heritabilities.

Myriads: P-value-based multiple testing correction.

Motivation: There are many multiple testing correction methods. Some of them are robust to various dependencies in the data while others are not. Some of the implementations have problems for managing high dimensional list of P-values as currently demanded by microarray and other omic technologies. Results: The program Myriads, formerly SGoF+, provides some of the most important P-value-based correction methods jointly with a test of dependency and a P-value simulator. Myriads easily manage hundreds of thousands of P-values. Availability and implementation: http://myriads.webs.uvigo.es. Contact: myriads@uvigo.es. Supplementary information: Supplementary data are available at Bioinformatics online.

Geographic and Temporal Variation of Distinct Intracellular Endosymbiont Strains of Wolbachia sp. in the Grasshopper Chorthippus parallelus: a Frequency-Dependent Mechanism?

Wolbachia is an intracellular endosymbiont that can produce a range of effects on host fitness, but the temporal dynamics of Wolbachia strains have rarely been experimentally evaluated. We compare interannual strain frequencies along a geographical region for understanding the forces that shape Wolbachia strain frequency in natural populations of its host, Chorthippus parallelus (Orthoptera, Acrididae). General linear models show that strain frequency changes significantly across geographical and temporal scales. Computer simulation allows to reject the compatibility of the observed patterns with either genetic drift or sampling errors. We use consecutive years to estimate total Wolbachia strain fitness. Our estimation of Wolbachia fitness is significant in most cases, within locality and between consecutive years, following a negatively frequency-dependent trend. Wolbachia spp. B and F strains show a temporal pattern of variation that is compatible with a negative frequency-dependent natural selection mechanism. Our results suggest that such a mechanism should be at least considered in future experimental and theoretical research strategies that attempt to understand Wolbachia biodiversity.

Impact of deep coalescence and recombination on the estimation of phylogenetic relationships among species using AFLP markers.

Deep coalescence and the nongenealogical pattern of descent caused by recombination have emerged as a common problem for phylogenetic inference at the species level. Here we use computer simulations to assess whether AFLP-based phylogenies are robust to the uncertainties introduced by these factors. Our results indicate that phylogenetic signal can prevail even in the face of extensive deep coalescence allowing recovering the correct species tree topology. The impact of recombination on tree accuracy was related to total tree depth and species effective population size. The correct tree topology could be recovered upon many simulation settings due to a trade-off between the conflicting signals resulting from intra-locus recombination and the benefits of the joint consideration of unlinked loci that better matched overall the true species tree. Errors in tree topology were not only determined by deep coalescence, but also by the timing of divergence and the tree-building errors arising from an insufficient number of characters. DNA sequences generally outperformed AFLPs upon any simulated scenario, but this difference in performance was nearly negligible when a sufficient number of AFLP characters were sampled. Our simulations suggest that the impact of deep coalescence and intra-locus recombination on the reliability of AFLP trees could be minimal for effective population sizes equal to or lower than 10,000 (typical of many vertebrates and tree plants) given tree depths above 0.02 substitutions per site.

Multiple hypothesis testing in proteomics: a strategy for experimental work.

In quantitative proteomics work, the differences in expression of many separate proteins are routinely examined to test for significant differences between treatments. This leads to the multiple hypothesis testing problem: when many separate tests are performed many will be significant by chance and be false positive results. Statistical methods such as the false discovery rate method that deal with this problem have been disseminated for more than one decade. However a survey of proteomics journals shows that such tests are not widely implemented in one commonly used technique, quantitative proteomics using two-dimensional electrophoresis. We outline a selection of multiple hypothesis testing methods, including some that are well known and some lesser known, and present a simple strategy for their use by the experimental scientist in quantitative proteomics work generally. The strategy focuses on the desirability of simultaneous use of several different methods, the choice and emphasis dependent on research priorities and the results in hand. This approach is demonstrated using case scenarios with experimental and simulated model data.

1LocusSim a mobile-friendly simulator for teaching population genetics.

Motivation: Biology students often struggle with the fundamental concepts of evolutionary genetics, including genetic drift, mutation and selection. To address this problem, 1LocusSim was developed to simulate the interaction of different factors, such as population size, mutation, selection and dominance, to study their effect on allelic frequency during evolution. With 1LocusSim, students can compare theoretical results with simulation outputs and solve and analyze different problems of population genetics. The 1LocusSim web has a responsive design which means that it has been specifically designed to be used on smartphones. Results: To demonstrate its use, I review the classical overdominance model of population genetics and highlight a characteristic that is often not explicitly stated. Specifically, it is emphasized that the equilibrium of the model does not depend on the homozygous selection coefficients but rather on the ratio of the selection coefficients. This is already clear from the classical formula but maybe not so much for students. Also it implies that the equilibrium can be expressed solely in terms of the dominance coefficient h. To verify this theoretical prediction, I utilize the simulator and calculate the equilibrium for the well-known case of sickle cell anemia. By utilizing this tool, students can learn at their own pace and convenience, anywhere and anytime. Availability and implementation: 1LocusSim if freely available at https://1LocusSim-biosdev.pythonanywhere.com/. Website implemented under the Bottle micro web-framework for Python, with all major browsers supported. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Multiple hypothesis testing in proteomics: A strategy for experimental work.

In quantitative proteomics work, the differences in expression of many separate proteins are routinely examined to test for significant differences between treatments. This leads to the multiple hypothesis testing problem: when many separate tests are performed many will be significant by chance and be false positive results. Statistical methods such as the false discovery rate (FDR) method that deal with this problem have been disseminated for more than one decade. However a survey of proteomics journals shows that such tests are not widely implemented in one commonly used technique, quantitative proteomics using two-dimensional electrophoresis (2-DE). We outline a selection of multiple hypothesis testing methods, including some that are well known and some lesser known, and present a simple strategy for their use by the experimental scientist in quantitative proteomics work generally. The strategy focuses on the desirability of simultaneous use of several different methods, the choice and emphasis dependent on research priorities and the results in hand. This approach is demonstrated using case scenarios with experimental and simulated model data.

Detecting Local Adaptation between North and South European Atlantic Salmon Populations.

Pollution and other anthropogenic effects have driven a decrease in Atlantic salmon (Salmo salar) in the Iberian Peninsula. The restocking effort carried out in the 1980s, with salmon from northern latitudes with the aim of mitigating the decline of native populations, failed, probably due to the deficiency in adaptation of foreign salmon from northern Europe to the warm waters of the Iberian Peninsula. This result would imply that the Iberian populations of Atlantic salmon have experienced local adaptation in their past evolutionary history, as has been described for other populations of this species and other salmonids. Local adaptation can occur by divergent selections between environments, favoring the fixation of alleles that increase the fitness of a population in the environment it inhabits relative to other alleles favored in another population. In this work, we compared the genomes of different populations from the Iberian Peninsula (Atlantic and Cantabric basins) and Scotland in order to provide tentative evidence of candidate SNPs responsible for the adaptive differences between populations, which may explain the failures of restocking carried out during the 1980s. For this purpose, the samples were genotyped with a 220,000 high-density SNP array (Affymetrix) specific to Atlantic salmon. Our results revealed potential evidence of local adaptation for North Spanish and Scottish populations. As expected, most differences concerned the comparison of the Iberian Peninsula with Scotland, although there were also differences between Atlantic and Cantabric populations. A high proportion of the genes identified are related to development and cellular metabolism, DNA transcription and anatomical structure. A particular SNP was identified within the NADP-dependent malic enzyme-2 (mMEP-2*), previously reported by independent studies as a candidate for local adaptation in salmon from the Iberian Peninsula. Interestingly, the corresponding SNP within the mMEP-2* region was consistent with a genomic pattern of divergent selection.

Evaluation of Existing Methods for High-Order Epistasis Detection.

Finding epistatic interactions among loci when expressing a phenotype is a widely employed strategy to understand the genetic architecture of complex traits in GWAS. The abundance of methods dedicated to the same purpose, however, makes it increasingly difficult for scientists to decide which method is more suitable for their studies. This work compares the different epistasis detection methods published during the last decade in terms of runtime, detection power and type I error rate, with a special emphasis on high-order interactions. Results show that in terms of detection power, the only methods that perform well across all experiments are the exhaustive methods, although their computational cost may be prohibitive in large-scale studies. Regarding non-exhaustive methods, not one could consistently find epistasis interactions when marginal effects are absent. If marginal effects are present, there are methods that perform well for high-order interactions, such as BADTrees, FDHE-IW, SingleMI or SNPHarvester. As for false-positive control, only SNPHarvester, FDHE-IW and DCHE show good results. The study concludes that there is no single epistasis detection method to recommend in all scenarios. Authors should prioritize exhaustive methods when sufficient computational resources are available considering the data set size, and resort to non-exhaustive methods when the analysis time is prohibitive.

Data and models from multi-model inference of non-random mating from an information theoretic approach.

This is a co-submission with Multi-model inference of non-random mating from an information theoretic approach [1]. These data corresponds to the complete simulated data set jointly with the set of models defined for analysing the data. The simulated data set was obtained using the program MateSim [2]. The simulated cases correspond to one-sex competition and mate choice models. For each simulation run, the population frequencies (premating individuals) and the sample of 500 mating pairs were generated randomly for a hypothetical trait with two classes at each sex. Some datasets represent larger population size species (n = 10 000) and the mating process was represented as a sampling with replacement, and the population frequencies were constant over the mating season. The minimum phenotype frequency (MPF) allowed was 0.1. Five different model cases were simulated, namely random mating, female competition with mate choice (with independent or compound parameters) and male competition with mate choice (with independent or compound parameters). Each case was simulated 1000 times. Other datasets represent monogamous species (with large or small population size) and the mating process was without replacement (from the point of view of the available phenotypes). These data sets were used to test the performance of the multi-model inference methodology proposed in [1]. The data may be useful for testing any new/old statistics for measuring sexual selection or assortative mating patterns.

A new multitest correction (SGoF) that increases its statistical power when increasing the number of tests.

BACKGROUND: The detection of true significant cases under multiple testing is becoming a fundamental issue when analyzing high-dimensional biological data. Unfortunately, known multitest adjustments reduce their statistical power as the number of tests increase. We propose a new multitest adjustment, based on a sequential goodness of fit metatest (SGoF), which increases its statistical power with the number of tests. The method is compared with Bonferroni and FDR-based alternatives by simulating a multitest context via two different kinds of tests: 1) one-sample t-test, and 2) homogeneity G-test. RESULTS: It is shown that SGoF behaves especially well with small sample sizes when 1) the alternative hypothesis is weakly to moderately deviated from the null model, 2) there are widespread effects through the family of tests, and 3) the number of tests is large. CONCLUSION: Therefore, SGoF should become an important tool for multitest adjustment when working with high-dimensional biological data.

GENOMEPOP: a program to simulate genomes in populations.

BACKGROUND: There are several situations in population biology research where simulating DNA sequences is useful. Simulation of biological populations under different evolutionary genetic models can be undertaken using backward or forward strategies. Backward simulations, also called coalescent-based simulations, are computationally efficient. The reason is that they are based on the history of lineages with surviving offspring in the current population. On the contrary, forward simulations are less efficient because the entire population is simulated from past to present. However, the coalescent framework imposes some limitations that forward simulation does not. Hence, there is an increasing interest in forward population genetic simulation and efficient new tools have been developed recently. Software tools that allow efficient simulation of large DNA fragments under complex evolutionary models will be very helpful when trying to better understand the trace left on the DNA by the different interacting evolutionary forces. Here I will introduce GenomePop, a forward simulation program that fulfills the above requirements. The use of the program is demonstrated by studying the impact of intracodon recombination on global and site-specific dN/dS estimation. RESULTS: I have developed algorithms and written software to efficiently simulate, forward in time, different Markovian nucleotide or codon models of DNA mutation. Such models can be combined with recombination, at inter and intra codon levels, fitness-based selection and complex demographic scenarios. CONCLUSION: GenomePop has many interesting characteristics for simulating SNPs or DNA sequences under complex evolutionary and demographic models. These features make it unique with respect to other simulation tools. Namely, the possibility of forward simulation under General Time Reversible (GTR) mutation or GTRxMG94 codon models with intra-codon recombination, arbitrary, user-defined, migration patterns, diploid or haploid models, constant or variable population sizes, etc. It also allows simulation of fitness-based selection under different distributions of mutational effects. Under the 2-allele model it allows the simulation of recombination hot-spots, the definition of different frequencies in different populations, etc. GenomePop can also manage large DNA fragments. In addition, it has a scaling option to save computation time when simulating large sequences and population sizes under complex demographic and evolutionary situations. These and many other features are detailed in its web page [1].

Disease progression and evolution of the HIV-1 env gene in 24 infected infants.

We have studied the relationship between disease progression and HIV-1 evolution in 24 infants classified as rapid or non-rapid progressors, during nearly the entire disease progression cycle from infection to AIDS. Specifically, we examined the temporal relationship between clinical status and changes in genetic diversity, divergence, selection and recombination at the C2V3C3 region of the env gene during a period of 3 years. Statistical analyses were performed using linear mixed models that are particularly well-suited for longitudinal studies in which repeated measures are taken from the same patients. We did not observe significant differences in genetic diversity or overall substitution rates between clinical categories. However, the nonsynonymous substitution rate per nonsynonymous site (dN) evolved differently between groups. Changes in dN explained the evolutionary slowdown of the dN/dS ratio in the rapid progressors, while in non-rapid progressors the dN/dS ratio continuously increased through time. The number of positively selected sites had limited power for predicting disease progression. Recombination rate estimates were different among groups, although not significantly in the linear mixed models analysis. They showed some power predicting clinical categories and, interestingly, they were significantly correlated with the frequency of positively selected sites. Overall, the results obtained confirm that viral adaptation in the C2V3C3 region of the env gene is related to disease progression, although the statistical characterization of such pattern seems rather difficult.

Simulation of genomes: a review.

There is an increasing role of population genetics in human genetic research linking empirical observations with hypotheses about sequence variation due to historical and evolutionary causes. In addition, the data sets are increasing in size, with genome-wide data becoming a common place in many empirical studies. As far as more information is available, it becomes clear that simplest hypotheses are not consistent with data. Simulations will provide the key tool to contrast complex hypotheses on real data by generating simulated data under the hypothetical historical and evolutionary conditions that we want to contrast. Undoubtedly, developing tools for simulating large sequences that at the same time allow simulate natural selection, recombination and complex demography patterns will be of great interest in order to better understanding the trace left on the DNA by different interacting evolutionary forces. Simulation tools will be also essential to evaluate the sampling properties of any statistics used on genome-wide association studies and to compare performance of methods applied at genome-wide scales. Several recent simulation tools have been developed. Here, we review some of the currently existing simulators which allow for efficient simulation of large sequences on complex evolutionary scenarios. In addition, we will point out future directions in this field which are already a key part of the current research in evolutionary biology and it seems that it will be a primary tool in the future research of genome and post-genomic biology.

Recombination favors the evolution of drug resistance in HIV-1 during antiretroviral therapy.

We studied the relationship between recombination and the fixation time of multiple drug resistance mutations after HIV-1 drug therapy, under a set of different realistic scenarios. We have generalized a previous model by Bretscher et al. [Bretscher, Althaus, Muller, Bonhoeffer, 2004. Recombination in HIV and the evolution of drug resistance: for better or for worse? Bioessays 26(2), 180-188] in order to explore different implementations of phenotypic mixing and more realistic demographic and selective regimes. Using computer simulations we show that the effect of recombination on the evolution of drug resistance depends strongly on the intensity of selection, as well as on the viral population size. Under the high selection pressure expected during antiretroviral therapy, the strength of the Hill-Robertson effect increases and recombination favors the evolution of resistance under a wide range of population sizes, independently of the sign of the epistatic interaction. Our results suggest that recombination plays an important role in the evolution of drug resistance in HIV-1 under various realistic scenarios. These findings could be taken into account in order to develop optimal HIV-1 drug treatments.

Evidence of recombination within human alpha-papillomavirus.

BACKGROUND: Human papillomavirus (HPV) has a causal role in cervical cancer with almost half a million new cases occurring each year. Presence of the carcinogenic HPV is necessary for the development of the invasive carcinoma of the genital tract. Therefore, persistent infection with carcinogenic HPV causes virtually all cervical cancers. Some aspects of the molecular evolution of this virus, as the putative importance of recombination in its evolutionary history, are an opened current question. In addition, recombination could also be a significant issue nowadays since the frequency of co-infection with more than one HPV type is not a rare event and, thus, new recombinant types could be currently being generated. RESULTS: We have used human alpha-PV sequences from the public database at Los Alamos National Laboratory to report evidence that recombination may exist in this virus. A model-based population genetic approach was used to infer the recombination signal from the HPV DNA sequences grouped attending to phylogenetic and epidemiological information, as well as to clinical manifestations. Our results agree with recently published ones that use a different methodology to detect recombination associated to the gene L2. In addition, we have detected significant recombination signal in the genes E6, E7, L2 and L1 at different groups, and importantly within the high-risk type HPV16. The method used has recently been shown to be one of the most powerful and reliable procedures to detect the recombination signal. CONCLUSION: We provide new support to the recent evidence of recombination in HPV. Additionally, we performed the recombination estimation assuming the best-fit model of nucleotide substitution and rate variation among sites, of the HPV DNA sequence sets. We found that the gene with recombination in most of the groups is L2 but the highest values were detected in L1 and E6. Gene E7 was recombinant only within the HPV16 type. The topic deserves further study because recombination is an important evolutionary mechanism that could have high impact both in pharmacogenomics (i.e. on the influence of genetic variation on the response to drugs) and for vaccine development.

HacDivSel: Two new methods (haplotype-based and outlier-based) for the detection of divergent selection in pairs of populations.

The detection of genomic regions involved in local adaptation is an important topic in current population genetics. There are several detection strategies available depending on the kind of genetic and demographic information at hand. A common drawback is the high risk of false positives. In this study we introduce two complementary methods for the detection of divergent selection from populations connected by migration. Both methods have been developed with the aim of being robust to false positives. The first method combines haplotype information with inter-population differentiation (FST). Evidence of divergent selection is concluded only when both the haplotype pattern and the FST value support it. The second method is developed for independently segregating markers i.e. there is no haplotype information. In this case, the power to detect selection is attained by developing a new outlier test based on detecting a bimodal distribution. The test computes the FST outliers and then assumes that those of interest would have a different mode. We demonstrate the utility of the two methods through simulations and the analysis of real data. The simulation results showed power ranging from 60-95% in several of the scenarios whilst the false positive rate was controlled below the nominal level. The analysis of real samples consisted of phased data from the HapMap project and unphased data from intertidal marine snail ecotypes. The results illustrate that the proposed methods could be useful for detecting locally adapted polymorphisms. The software HacDivSel implements the methods explained in this manuscript.

A novel method for estimating the strength of positive mating preference by similarity in the wild.

Mating preference can be a driver of sexual selection and assortative mating and is, therefore, a key element in evolutionary dynamics. Positive mating preference by similarity is the tendency for the choosy individual to select a mate which possesses a similar variant of a trait. Such preference can be modelled using Gaussian-like mathematical functions that describe the strength of preference, but such functions cannot be applied to empirical data collected from the field. As a result, traditionally, mating preference is indirectly estimated by the degree of assortative mating (using Pearson's correlation coefficient, r) in wild captured mating pairs. Unfortunately, r and similar coefficients are often biased due to the fact that different variants of a given trait are nonrandomly distributed in the wild, and pooling of mating pairs from such heterogeneous samples may lead to "false-positive" results, termed "the scale-of-choice effect" (SCE). Here we provide two new estimators of mating preference (Crough and Cscaled) derived from Gaussian-like functions which can be applied to empirical data. Computer simulations demonstrated that r coefficient showed robust estimations properties of mating preference but it was severely affected by SCE, Crough showed reasonable estimation properties and it was little affected by SCE, while Cscaled showed the best properties at infinite sample sizes and it was not affected by SCE but failed at biological sample sizes. We recommend using Crough combined with the r coefficient to infer mating preference in future empirical studies.

JMATING: a software for the analysis of sexual selection and sexual isolation effects from mating frequency data.

BACKGROUND: Many different sexual isolation and sexual selection statistics have been proposed in the past. However, there is no available software that implements all these statistical estimators and their corresponding tests for the study of mating behaviour. RESULTS: JMATING is an easy-to-use program developed in Java for the analysis of mating frequency data to study sexual selection and sexual isolation effects from laboratory experiments as well as descriptive studies accomplished in the wild. The software allows the re-organization of the data previous to the analysis, the estimation of the most important estimators, and a battery of complementary statistical tests. CONCLUSION: JMATING is the first complete and versatile software for the analyses of mating frequency data. It is available at http://www.uvigo.es/webs/c03/webc03/XENETICA/XB2/JMsoft.htm and requires the Java runtime environment.

The scale-of-choice effect and how estimates of assortative mating in the wild can be biased due to heterogeneous samples.

The mode in which sexual organisms choose mates is a key evolutionary process, as it can have a profound impact on fitness and speciation. One way to study mate choice in the wild is by measuring trait correlation between mates. Positive assortative mating is inferred when individuals of a mating pair display traits that are more similar than those expected under random mating while negative assortative mating is the opposite. A recent review of 1134 trait correlations found that positive estimates of assortative mating were more frequent and larger in magnitude than negative estimates. Here, we describe the scale-of-choice effect (SCE), which occurs when mate choice exists at a smaller scale than that of the investigator's sampling, while simultaneously the trait is heterogeneously distributed at the true scale-of-choice. We demonstrate the SCE by Monte Carlo simulations and estimate it in two organisms showing positive (Littorina saxatilis) and negative (L. fabalis) assortative mating. Our results show that both positive and negative estimates are biased by the SCE by different magnitudes, typically toward positive values. Therefore, the low frequency of negative assortative mating observed in the literature may be due to the SCE's impact on correlation estimates, which demands new experimental evaluation.