A fast and powerful W-test for pairwise epistasis testing

Epistasis plays an essential role in the development of complex diseases. Interaction methods face common challenge of seeking a balance between persistent power, model complexity, computation efficiency, and validity of identified bio-markers. We introduce a novel W-test to identify pairwise epistasis effect, which measures the distributional difference between cases and controls through a combined log odds ratio. The test is model-free, fast, and inherits a Chi-squared distribution with data adaptive degrees of freedom. No permutation is needed to obtain the P-values. Simulation studies demonstrated that the W-test is more powerful in low frequency variants environment than alternative methods, which are the Chi-squared test, logistic regression and multifactor-dimensionality reduction (MDR). In two independent real bipolar disorder genome-wide associations (GWAS) datasets, the W-test identified significant interactions pairs that can be replicated, including SLIT3-CENPN, SLIT3-TMEM132D, CNTNAP2-NDST4 and CNTCAP2-RTN4R. The genes in the pairs play central roles in neurotransmission and synapse formation. A majority of the identified loci are undiscoverable by main effect and are low frequency variants. The proposed method offers a powerful alternative tool for mapping the genetic puzzle underlying complex disorders.     

A fast, powerful method for detecting identity by descent

We present a method, fastIBD, for finding tracts of identity by descent (IBD) between pairs of individuals. FastIBD can be applied to thousands of samples across genome-wide SNP data and is significantly more powerful for finding short tracts of IBD than existing methods for finding IBD tracts in such data. We show that fastIBD can detect facets of population structure that are not revealed by other methods. In the Wellcome Trust Case Control Consortium bipolar disorder case-control data, we find a genome-wide excess of IBD in case-case pairs of individuals compared to control-control pairs. We show that this excess can be explained by the geographical clustering of cases. We also show that it is possible to use fastIBD to generate highly accurate estimates of genome-wide IBD sharing between pairs of distant relatives. This is useful for estimation of relationship and for adjusting for relatedness in association studies. FastIBD is incorporated in the freely available Beagle software package.     

A testing framework for identifying susceptibility genes in the presence of epistasis

An efficient testing strategy called the "focused interaction testing framework" (FITF) was developed to identify susceptibility genes involved in epistatic interactions for case-control studies of candidate genes. In the FITF approach, likelihood-ratio tests are performed in stages that increase in the order of interaction considered. Joint tests of main effects and interactions are performed conditional on significant lower-order effects. A reduction in the number of tests performed is achieved by prescreening gene combinations with a goodness-of-fit chi2 statistic that depends on association among candidate genes in the pooled case-control group. Multiple testing is accounted for by controlling false-discovery rates. Simulation analysis demonstrated that the FITF approach is more powerful than marginal tests of candidate genes. FITF also outperformed multifactor dimensionality reduction when interactions involved additive, dominant, or recessive genes. In an application to asthma case-control data from the Children's Health Study, FITF identified a significant multilocus effect between the nicotinamide adenine dinucleotide (phosphate) reduced:quinone oxidoreductase gene (NQO1), myeloperoxidase gene (MPO), and catalase gene (CAT) (unadjusted P = .00026), three genes that are involved in the oxidative stress pathway. In an independent data set consisting primarily of African American and Asian American children, these three genes also showed a significant association with asthma status (P = .0008).     

Modeling the evolution of genetic architecture: A continuum of alleles model with pairwise AxA epistasis

A continuum of alleles model with pair-wise AxA epistasis is proposed and its transmission genetic, and variational properties are analysed. The basic idea is that genes control the values of underlying variables, which affect the genotypic value of phenotypic characters proportional to a "scaling factor". Epistasis is the influence of one gene on the average effect of another gene. In this model, epistasis is introduced as a mutational effect of one gene on the scaling factors of another gene. In accordance with empirical results, the model assumes that the average direct effect of mutations is zero, as is the average epistatic effect. The model predicts that, on average, a mutation at one locus increases the expected mutational variance of mutations at another interacting locus. The increase in mutational variance is predicted to be equal to the variance of the pair-wise epistatic effects. This result is consistent with the observation that mutant phenotypes tend to be more variable than the wildtype phenotype. Another generic result of this model is that the frequency of canalizing mutations can at most be equal to the frequency of de-canalizing mutations. Furthermore, it is predicted that the mutational variance of a character increases at least linearly with the size of the character; hence this model is scale variant. In the case of two characters it is shown that the dimensionality of the locus-specific mutational effect distribution is invariant, i.e. the rank of the mutational covariance matrix M is invariant. While in additive models the mutational covariance matrix is always and entirely invariant, the invariance in the case of epistatic models is unexpected. Epistatic interactions can change the magnitude of the mutational (co)variances at a locus and can thus influence the structure of the mutational covariance matrix. However, in the present model the dimensionality of the mutational effect distribution remains the same. A consequence of this result is that, in this model, the genetic architecture of a set of characters is always evolvable i.e. no hard constraints can evolve.     

An entropy-based approach for testing genetic epistasis underlying complex diseases

The genetic basis of complex diseases is expected to be highly heterogeneous, with complex interactions among multiple disease loci and environment factors. Due to the multi-dimensional property of interactions among large number of genetic loci, efficient statistical approach has not been well developed to handle the high-order epistatic complexity. In this article, we introduce a new approach for testing genetic epistasis in multiple loci using an entropy-based statistic for a case-only design. The entropy-based statistic asymptotically follows a χ² distribution. Computer simulations show that the entropy-based approach has better control of type I error and higher power compared to the standard χ² test. Motivated by a schizophrenia data set, we propose a method for measuring and testing the relative entropy of a clinical phenotype, through which one can test the contribution or interaction of multiple disease loci to a clinical phenotype. A sequential forward selection procedure is proposed to construct a genetic interaction network which is illustrated through a tree-based diagram. The network information clearly shows the relative importance of a set of genetic loci on a clinical phenotype. To show the utility of the new entropy-based approach, it is applied to analyze two real data sets, a schizophrenia data set and a published malaria data set. Our approach provides a fast and testable framework for genetic epistasis study in a case-only design.     

Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative

We revisit statistical tests for branches of evolutionary trees reconstructed upon molecular data. A new, fast, approximate likelihood-ratio test (aLRT) for branches is presented here as a competitive alternative to nonparametric bootstrap and Bayesian estimation of branch support. The aLRT is based on the idea of the conventional LRT, with the null hypothesis corresponding to the assumption that the inferred branch has length 0. We show that the LRT statistic is asymptotically distributed as a maximum of three random variables drawn from the chi(0)2 + chi(1)2 distribution. The new aLRT of interior branch uses this distribution for significance testing, but the test statistic is approximated in a slightly conservative but practical way as 2(l1- l2), i.e., double the difference between the maximum log-likelihood values corresponding to the best tree and the second best topological arrangement around the branch of interest. Such a test is fast because the log-likelihood value l2 is computed by optimizing only over the branch of interest and the four adjacent branches, whereas other parameters are fixed at their optimal values corresponding to the best ML tree. The performance of the new test was studied on simulated 4-, 12-, and 100-taxon data sets with sequences of different lengths. The aLRT is shown to be accurate, powerful, and robust to certain violations of model assumptions. The aLRT is implemented within the algorithm used by the recent fast maximum likelihood tree estimation program PHYML (Guindon and Gascuel, 2003).     

A pairwise relatedness estimator for polyploids

Studies in genetics and ecology often require estimates of relatedness coefficients based on genetic marker data. Many diploid estimators have been developed using either method‐of‐moments or maximum‐likelihood estimates. However, there are no relatedness estimators for polyploids. The development of a moment estimator for polyploids with polysomic inheritance, which simultaneously incorporates the two‐gene relatedness coefficient and various ‘higher‐order’ coefficients, is described here. The performance of the estimator is compared to other estimators under a variety of conditions. When using a small number of loci, the estimator is biased because of an increase in ill‐conditioned matrices. However, the estimator becomes asymptotically unbiased with large numbers of loci. The ambiguity of polyploid heterozygotes (when balanced heterozygotes cannot be distinguished from unbalanced heterozygotes) is also considered; as with low numbers of loci, genotype ambiguity leads to bias. A software, PolyRelatedness , implementing this method and supporting a maximum ploidy of 8 is provided.     

A powerful truncated tail strength method for testing multiple null hypotheses in one dataset

This article re-analyses a prey-predator model with a refuge introduced by one of the founders of population ecology Gause and his co-workers to explain discrepancies between their observations and predictions of the Lotka-Volterra prey-predator model. They replaced the linear functional response used by Lotka and Volterra by a saturating functional response with a discontinuity at a critical prey density. At concentrations below this critical density prey were effectively in a refuge while at a higher densities they were available to predators. Thus, their functional response was of the Holling type III. They analyzed this model and predicted existence of a limit cycle in predator-prey dynamics. In this article I show that their model is ill posed, because trajectories are not well defined. Using the Filippov method, I define and analyze solutions of the Gause model. I show that depending on parameter values, there are three possibilities: (1) trajectories converge to a limit cycle, as predicted by Gause, (2) trajectories converge to an equilibrium, or (3) the prey population escapes predator control and grows to infinity.     

A powerful strategy to account for multiple testing in the context of haplotype analysis

Haplotypes--that is, linear arrangements of alleles on the same chromosome that were inherited as a unit--are expected to carry important information in the context of association fine mapping of complex diseases. In consideration of a set of tightly linked markers, there is an enormous number of different marker combinations that can be analyzed. Therefore, a severe multiple-testing problem is introduced. One method to deal with this problem is Bonferroni correction by the number of combinations that are considered. Bonferroni correction is appropriate for independent tests but will result in a loss of power in the presence of linkage disequilibrium in the region. A second method is to perform simulations. It is unfortunate that most methods of haplotype analysis already require simulations to obtain an uncorrected P value for a specific marker combination. Thus, it seems that nested simulations are necessary to obtain P values that are corrected for multiple testing, which, apparently, limits the applicability of this approach because of computer running-time restrictions. Here, an algorithm is described that avoids such nested simulations. We check the validity of our approach under two disease models for haplotype analysis of family data. The true type I error rate of our algorithm corresponds to the nominal significance level. Furthermore, we observe a strong gain in power with our method to obtain the global P value, compared with the Bonferroni procedure to calculate the global P value. The method described here has been implemented in the latest update of our program FAMHAP.     

A new powerful non-parametric two-stage approach for testing multiple phenotypes in family-based association studies

We introduce a new powerful nonparametric testing strategy for family-based association studies in which multiple quantitative traits are recorded and the phenotype with the strongest genetic component is not known prior to the analysis. In the first stage, using a population-based test based on the generalized estimating equation approach, we test all recorded phenotypes for association with the marker locus without biasing the nominal significance level of the later family-based analysis. In the second stage the phenotype with the smallest p value is selected and tested by a family-based association test for association with the marker locus. This strategy is robust against population admixture and stratification and does not require any adjustment for multiple testing. We demonstrate the advantages of this testing strategy over standard methodology in a simulation study. The practical importance of our testing strategy is illustrated by applications to the Childhood Asthma Management Program asthma data sets.     

Connected populations for detecting quantitative trait loci and testing for epistasis: an application in maize

Quantitative trait loci (QTL) detection experiments have often been restricted to large biallelic populations. Use of connected multiparental crosses has been proposed to increase the genetic variability addressed and to test for epistatic interactions between QTL and the genetic background. We present here the results of a QTL detection performed on six connected F₂ populations of 150 F_2:3 families each, derived from four maize inbreds and evaluated for three traits of agronomic interest. The QTL detection was carried out by composite interval mapping on each population separately, then on the global design either by taking into account the connections between populations or not. Epistatic interactions between loci and with the genetic background were tested. Taking into account the connections between populations increased the number of QTL detected and the accuracy of QTL position estimates. We detected many epistatic interactions, particularly for grain yield QTL (R ² increase of 9.6%). Use of connections for the QTL detection also allowed a global ranking of alleles at each QTL. Allelic relationships and epistasis both contribute to the lack of consistency for QTL positions observed among populations, in addition to the limited power of the tests. The potential benefit of assembling favorable alleles by marker-assisted selection are discussed.     

tuple_plot: fast pairwise nucleotide sequence comparison with noise suppression

SUMMARY: The program tuple_plot identifies and visualizes local similarities between two genomic sequences, typically 100 kb or longer, by applying the well-known dotplot principle. A dictionary of sequence words built from the input sequences serves to construct a task-specific expectancy model that is used to attribute significance values to pairwise word hits. The dictionary-based approach allows fast computation, the computation time scaling to O(N log N), depending on the size of the input sequences. The proposed scoring scheme appreciably increases the signal-to-noise ratio and may help to improve other word-based sequence comparison approaches. AVAILABILITY: tuple_plot is available at http://genome.fli-leibniz.de/software.html and may be used under GNU public license.     

A maximum-likelihood estimation of pairwise relatedness for autopolyploids

Relatedness between individuals is central to ecological genetics. Multiple methods are available to quantify relatedness from molecular data, including method-of-moment and maximum-likelihood estimators. We describe a maximum-likelihood estimator for autopolyploids, and quantify its statistical performance under a range of biologically relevant conditions. The statistical performances of five additional polyploid estimators of relatedness were also quantified under identical conditions. When comparing truncated estimators, the maximum-likelihood estimator exhibited lower root mean square error under some conditions and was more biased for non-relatives, especially when the number of alleles per loci was low. However, even under these conditions, this bias was reduced to be statistically insignificant with more robust genetic sampling. We also considered ambiguity in polyploid heterozygote genotyping and developed a weighting methodology for candidate genotypes. The statistical performances of three polyploid estimators under both ideal and actual conditions (including inbreeding and double reduction) were compared. The software package POLYRELATEDNESS is available to perform this estimation and supports a maximum ploidy of eight.     

A network model for the correlation between epistasis and genomic complexity

The study of genetic interactions (epistasis) is central to the understanding of genome organization and evolution. A general correlation between epistasis and genomic complexity has been recently shown, such that in simpler genomes epistasis is antagonistic on average (mutational effects tend to cancel each other out), whereas a transition towards synergistic epistasis occurs in more complex genomes (mutational effects strengthen each other). Here, we use a simple network model to identify basic features explaining this correlation. We show that, in small networks with multifunctional nodes, lack of redundancy, and absence of alternative pathways, epistasis is antagonistic on average. In contrast, lack of multi-functionality, high connectivity, and redundancy favor synergistic epistasis. Moreover, we confirm the previous finding that epistasis is a covariate of mutational robustness: in less robust networks it tends to be antagonistic whereas in more robust networks it tends to be synergistic. We argue that network features associated with antagonistic epistasis are typically found in simple genomes, such as those of viruses and bacteria, whereas the features associated with synergistic epistasis are more extensively exploited by higher eukaryotes.     

Mathematical frameworks for phenotypical selection and epistasis

A mathematical approach to interactions between genotypes and phenotypes in a multilocus multiallele population is developed. No a priori information on a fitness function is required. In particular, some structural definitions of epistasis and the position effect are given in terms of a decomposition of phenotypical structures. On this base a distance to the additive non-epistasis is introduced and an explicit formula for it is obtained. A class of phenotypical structures including multilocus dominance is described in terms of directed graphs. The evolutionary equations are adjusted to a fitness function compatible with a phenotypical structure. Some results on the finiteness of the equilibria set are presented.     

Testing for pairwise independence

This article presents a method for testing the hypothesis of mutual pairwise independence of k events. The method, which is based on the weighted least squares approach (Grizzle, Starmer, and Koch, 1969, Biometrics 25, 489-504) can be generalized to two types of incomplete data: the multiple-recapture census, where one of the cells of the corresponding 2k contingency table cannot be observed, and situations allowing an "unknown" response to the question designed to determine whether an event has occurred.