Case-control genome-wide association studies (CC-GWAS) might provide valuable clues to the underlying pathophysiologic mechanisms of complex diseases, such as neurodegenerative disease and cancer. A commonly overlooked complication is that multiple distinct disease states might present with the same set of symptoms and hence share a clinical diagnosis. These disease states can only be distinguished based on a biomarker evaluation that might not be feasible in the whole set of cases in the large number of samples that are typically needed for CC-GWAS. Instead, the biomarkers are measured on a subset of cases. Or an external reliability study estimates the frequencies of the disease states of interest within the clinically diagnosed set of cases. These frequencies often vary by the genetic and/or nongenetic variables. We derive a simple approximation that relates the genetic effect estimates obtained in a traditional logistic regression model with the clinical diagnosis as the outcome variable to the genetic effect estimates in the relationship to the true disease state of interest. We performed simulation studies to assess the accuracy of the approximation that we have derived. We next applied the derived approximation to the analysis of the genetic basis of the innate immune system of Alzheimer's disease.

Case-control Genome-Wide Association Studies (GWAS) provide a rich resource for studying the genetic architecture of complex diseases. A key is to elucidate how the genetic effects vary by the environment, what is traditionally defined by Gene-Environment interactions (GxE). The overlooked complication is that multiple, distinct pathophysiologic mechanisms may lead to the same clinical diagnosis and often these mechanisms have distinct genetic bases. In this paper, we first show that using the clinically diagnosed status can lead to severely biased estimates of GxE interactions in situations when the frequency of the pathologic diagnosis of interest, as compared to other diagnoses, depends on the environment. We then propose a pseudo-likelihood solution to correct the bias. Finally, we demonstrate our method in extensive simulations and in a GWAS of Alzheimer's disease.

One of the most important research areas in case-control Genome-Wide Association Studies is to determine how the effect of a genotype varies across the environment or to measure the gene-environment interaction (G × E). We consider the scenario when some of the "healthy" controls actually have the disease and when the frequency of these latent cases varies by the environmental variable of interest. In this scenario, performing logistic regression with the clinically diagnosed disease status as an outcome variable and will result in biased estimates of G × E interaction. Here, we derive a general theoretical approximation to the bias in the estimates of the G × E interaction and show, through extensive simulation, that this approximation is accurate in finite samples. Moreover, we apply this approximation to evaluate the bias in the effect estimates of the genetic variants related to mitochondrial proteins a large-scale prostate cancer study.

Background and objectives

Biological systems with intertwined feedback loops pose a challenge to mathematical modeling efforts. Moreover, rare events, such as mutation and extinction, complicate system dynamics. Stochastic simulation algorithms are useful in generating time-evolution trajectories for these systems because they can adequately capture the influence of random fluctuations and quantify rare events. We present a simple and flexible package, BioSimulator.jl, for implementing the Gillespie algorithm, τ-leaping, and related stochastic simulation algorithms. The objective of this work is to provide scientists across domains with fast, user-friendly simulation tools.

Methods

We used the high-performance programming language Julia because of its emphasis on scientific computing. Our software package implements a suite of stochastic simulation algorithms based on Markov chain theory. We provide the ability to (a) diagram Petri Nets describing interactions, (b) plot average trajectories and attached standard deviations of each participating species over time, and (c) generate frequency distributions of each species at a specified time.

Results

BioSimulator.jl's interface allows users to build models programmatically within Julia. A model is then passed to the simulate routine to generate simulation data. The built-in tools allow one to visualize results and compute summary statistics. Our examples highlight the broad applicability of our software to systems of varying complexity from ecology, systems biology, chemistry, and genetics.

Conclusion

The user-friendly nature of BioSimulator.jl encourages the use of stochastic simulation, minimizes tedious programming efforts, and reduces errors during model specification.

Background

Many important applications in bioinformatics, including sequence alignment and protein family profiling, employ sequence weighting schemes to mitigate the effects of non-independence of homologous sequences and under- or over-representation of certain taxa in a dataset. These schemes aim to assign high weights to sequences that are 'novel' compared to the others in the same dataset, and low weights to sequences that are over-represented.

Results

We formalise this principle by rigorously defining the evolutionary 'novelty' of a sequence within an alignment. This results in new sequence weights that we call 'phylogenetic novelty scores'. These scores have various desirable properties, and we showcase their use by considering, as an example application, the inference of character frequencies at an alignment column-important, for example, in protein family profiling. We give computationally efficient algorithms for calculating our scores and, using simulations, show that they are versatile and can improve the accuracy of character frequency estimation compared to existing sequence weighting schemes.

Conclusions

Our phylogenetic novelty scores can be useful when an evolutionarily meaningful system for adjusting for uneven taxon sampling is desired. They have numerous possible applications, including estimation of evolutionary conservation scores and sequence logos, identification of targets in conservation biology, and improving and measuring sequence alignment accuracy.

Abstract Background Terabyte-scale collections of string-encoded data are expected from consortia efforts such as the Human Microbiome Project http://nihroadmap.nih.gov/hmp. Intra- and inter-project data similarity searches are enabled by rapid k-mer matching strategies. Software applications for sequence database partitioning, guide tree estimation, molecular classification and alignment acceleration have benefited from embedded k-mer searches as sub-routines. However, a rapid, general-purpose, open-source, flexible, stand-alone k-mer tool has not been available. Results Here we present a stand-alone utility, Simrank, which allows users to rapidly identify database strings the most similar to query strings. Performance testing of Simrank and related tools against DNA, RNA, protein and human-languages found Simrank 10X to 928X faster depending on the dataset. Conclusions Simrank provides molecular ecologists with a high-throughput, open source choice for comparing large sequence sets to find similarity.

Intestinal microbiota play a significant role in nutrient metabolism, modulation of the immune system, obesity, and possibly in carcinogenesis, although the underlying mechanisms resulting in disease or impacts on longevity caused by different intestinal microbiota are mostly unknown. Herein we use isogenic Atm-deficient and wild type mice as models to interrogate changes in the metabolic profiles of urine and feces of these mice, which are differing in their intestinal microbiota. Using high resolution mass spectrometry approach we show that the composition of intestinal microbiota modulates specific metabolic perturbations resulting in a possible alleviation of a glycolytic phenotype. Metabolites including 3-methylbutyrolactone, kyneurenic acid and 3-methyladenine known to be onco-protective are elevated in Atm-deficient and wild type mice with restricted intestinal microbiota. Thus our approach has broad applicability to study the direct influence of gut microbiome on host metabolism and resultant phenotype. These results for the first time suggest a possible correlation of metabolic alterations and carcinogenesis, modulated by intestinal microbiota in A-T mice.