Longitudinal hierarchical Bayesian models of covariate effects on airway and alveolar nitric oxide.

Biomarkers such as exhaled nitric oxide (FeNO), a marker of airway inflammation, have applications in the study of chronic respiratory disease where longitudinal studies of within-participant changes in the biomarker are particularly relevant. A cutting-edge approach to assessing FeNO, called multiple flow FeNO, repeatedly assesses FeNO across a range of expiratory flow rates at a single visit and combines these data with a deterministic model of lower respiratory tract NO to estimate parameters quantifying airway wall and alveolar NO sources. Previous methodological work for multiple flow FeNO has focused on methods for data from a single participant or from cross-sectional studies.

Functional human genes typically exhibit epigenetic conservation.

Recent DepMap CRISPR-Cas9 single gene disruptions have identified genes more essential to proliferation in tissue culture. It would be valuable to translate these finding with measurements more practical for human tissues. Here we show that DepMap essential genes and other literature curated functional genes exhibit cell-specific preferential epigenetic conservation when DNA methylation measurements are compared between replicate cell lines and between intestinal crypts from the same individual.

Hierarchical Bayesian estimation of covariate effects on airway and alveolar nitric oxide.

Exhaled breath biomarkers are an important emerging field. The fractional concentration of exhaled nitric oxide (FeNO) is a marker of airway inflammation with clinical and epidemiological applications (e.g.

Bayesian parameter estimation for automatic annotation of gene functions using observational data and phylogenetic trees.

Gene function annotation is important for a variety of downstream analyses of genetic data. But experimental characterization of function remains costly and slow, making computational prediction an important endeavor. Phylogenetic approaches to prediction have been developed, but implementation of a practical Bayesian framework for parameter estimation remains an outstanding challenge.

Genomic diversity and genome-wide association analysis related to yield and fatty acid composition of wild American oil palm.

Existing Elaeis guineensis cultivars lack sufficient genetic diversity due to extensive breeding. Harnessing variation in wild crop relatives is necessary to expand the breadth of agronomically valuable traits. Using RAD sequencing, we examine the natural diversity of wild American oil palm populations (Elaeis oleifera), a sister species of the cultivated Elaeis guineensis oil palm.

: a web application to identify the most similar mutational signature using shiny.

There are two frameworks for characterizing mutational signatures which are commonly used to describe the nucleotide patterns that arise from mutational processes. Estimated mutational signatures from fitting these two methods in human cancer can be found online, in the Catalogue Of Somatic Mutations In Cancer (COSMIC) website or a GitHub repository. The two frameworks make differing assumptions regarding independence of base pairs and for that reason may produce different results.

Epigenetic Conservation Is a Beacon of Function: An Analysis Using Methcon5 Software for Studying Gene Methylation.

Purpose: Different epigenetic configurations allow one genome to develop into multiple cell types. Although the rules governing what epigenetic features confer gene expression are increasingly being understood, much remains uncertain. Here, we used a novel software package, Methcon5, to explore whether the principle of biologic conservation can be used to identify expressed genes.

Impact of different fixed flow sampling protocols on flow-independent exhaled nitric oxide parameter estimates using the Bayesian dynamic two-compartment model.

Exhaled nitric oxide (FeNO) is an established respiratory biomarker with clinical applications in the diagnosis and management of asthma. Because FeNO depends strongly on the flow (exhalation) rate, early protocols specified that measurements should be taken when subjects exhaled at a fixed rate of 50 ml/s. Subsequently, multiple flow (or "extended") protocols were introduced which measure FeNO across a range of fixed flow rates, allowing estimation of parameters including C NO and C NO which partition the physiological sources of NO into proximal airway wall tissue and distal alveolar regions (respectively).

Mutational signatures in colon cancer.

Objective: Recently, many tumor sequencing studies have inferred and reported on mutational signatures, short nucleotide patterns at which particular somatic base substitutions appear more often. A number of signatures reflect biological processes in the patient and factors associated with cancer risk. Our goal is to infer mutational signatures appearing in colon cancer, a cancer for which environmental risk factors vary by cancer subtype, and compare the signatures to those in adult stem cells from normal colon.

HiLDA: a statistical approach to investigate differences in mutational signatures.

We propose a hierarchical latent Dirichlet allocation model (HiLDA) for characterizing somatic mutation data in cancer. The method allows us to infer mutational patterns and their relative frequencies in a set of tumor mutational catalogs and to compare the estimated frequencies between tumor sets. We apply our method to two datasets, one containing somatic mutations in colon cancer by the time of occurrence, before or after tumor initiation, and the second containing somatic mutations in esophageal cancer by sex, age, smoking status, and tumor site.

Bayesian model selection for the Drosophila gap gene network.

Background: The gap gene system controls the early cascade of the segmentation pathway in Drosophila melanogaster as well as other insects. Owing to its tractability and key role in embryo patterning, this system has been the focus for both computational modelers and experimentalists. The gap gene expression dynamics can be considered strictly as a one-dimensional process and modeled as a system of reaction-diffusion equations.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants.

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis.

Threshold response to stochasticity in morphogenesis.

During development of biological organisms, multiple complex structures are formed. In many instances, these structures need to exhibit a high degree of order to be functional, although many of their constituents are intrinsically stochastic. Hence, it has been suggested that biological robustness ultimately must rely on complex gene regulatory networks and clean-up mechanisms.

Conservation of social effects (Ψ between two species of despite reversal of sexual dimorphism.

Indirect genetic effects (IGEs) describe the effect of the genes of social partners on the phenotype of a focal individual. Here, we measure indirect genetic effects using the "coefficient of interaction" (Ψ) to test whether Ψ evolved between and . We compare Ψ for locomotion between ethanol and nonethanol environments in both species, but only utilizes ethanol ecologically.

Translating natural genetic variation to gene expression in a computational model of the Drosophila gap gene regulatory network.

Annotating the genotype-phenotype relationship, and developing a proper quantitative description of the relationship, requires understanding the impact of natural genomic variation on gene expression. We apply a sequence-level model of gap gene expression in the early development of Drosophila to analyze single nucleotide polymorphisms (SNPs) in a panel of natural sequenced D. melanogaster lines.

Basic reversal-learning capacity in flies suggests rudiments of complex cognition.

The most basic models of learning are reinforcement learning models (for instance, classical and operant conditioning) that posit a constant learning rate; however many animals change their learning rates with experience. This process is sometimes studied by reversing an existing association between cues and rewards, and measuring the rate of relearning. Augmented reversal-learning, where learning rates increase with practice, can be an important component of behavioral flexibility; and may provide insight into higher cognition.

Sure enough: efficient Bayesian learning and choice.

Probabilistic decision-making is a general phenomenon in animal behavior, and has often been interpreted to reflect the relative certainty of animals' beliefs. Extensive neurological and behavioral results increasingly suggest that animal beliefs may be represented as probability distributions, with explicit accounting of uncertainty. Accordingly, we develop a model that describes decision-making in a manner consistent with this understanding of neuronal function in learning and conditioning.

Social effects for locomotion vary between environments in Drosophila melanogaster females.

Despite strong purifying or directional selection, variation is ubiquitous in populations. One mechanism for the maintenance of variation is indirect genetic effects (IGEs), as the fitness of a given genotype will depend somewhat on the genes of its social partners. IGEs describe the effect of genes in social partners on the expression of the phenotype of a focal individual.

Early mutation bursts in colorectal tumors.

Tumor growth is an evolutionary process involving accumulation of mutations, copy number alterations, and cancer stem cell (CSC) division and differentiation. As direct observation of this process is impossible, inference regarding when mutations occur and how stem cells divide is difficult. However, this ancestral information is encoded within the tumor itself, in the form of intratumoral heterogeneity of the tumor cell genomes.

GWASeq: targeted re-sequencing follow up to GWAS.

Background: For the last decade the conceptual framework of the Genome-Wide Association Study (GWAS) has dominated the investigation of human disease and other complex traits. While GWAS have been successful in identifying a large number of variants associated with various phenotypes, the overall amount of heritability explained by these variants remains small. This raises the question of how best to follow up on a GWAS, localize causal variants accounting for GWAS hits, and as a consequence explain more of the so-called "missing" heritability.

Exact likelihood-free Markov chain Monte Carlo for elliptically contoured distributions.

Recent results in Markov chain Monte Carlo (MCMC) show that a chain based on an unbiased estimator of the likelihood can have a stationary distribution identical to that of a chain based on exact likelihood calculations. In this paper we develop such an estimator for elliptically contoured distributions, a large family of distributions that includes and generalizes the multivariate normal. We then show how this estimator, combined with pseudorandom realizations of an elliptically contoured distribution, can be used to run MCMC in a way that replicates the stationary distribution of a likelihood based chain, but does not require explicit likelihood calculations.

al3c: high-performance software for parameter inference using Approximate Bayesian Computation.

Motivation: The development of Approximate Bayesian Computation (ABC) algorithms for parameter inference which are both computationally efficient and scalable in parallel computing environments is an important area of research. Monte Carlo rejection sampling, a fundamental component of ABC algorithms, is trivial to distribute over multiple processors but is inherently inefficient. While development of algorithms such as ABC Sequential Monte Carlo (ABC-SMC) help address the inherent inefficiencies of rejection sampling, such approaches are not as easily scaled on multiple processors.

Simulation-based Bayesian Analysis of Complex Data.

Our ability to collect large datasets is growing rapidly. Such richness of data offers great promise in terms of addressing detailed scientific questions in great depth. However, this benefit is not without scientific difficulty: many traditional analysis methods become computationally intractable for very large datasets.