Parameter estimation and identifiability analysis for a bivalent analyte model of monoclonal antibody-antigen binding.

Surface plasmon resonance (SPR) is an extensively used technique to characterize antigen-antibody interactions. Affinity measurements by SPR typically involve testing the binding of antigen in solution to monoclonal antibodies (mAbs) immobilized on a chip and fitting the kinetics data using 1:1 Langmuir binding model to derive rate constants. However, when it is necessary to immobilize antigens instead of the mAbs, a bivalent analyte (1:2) binding model is required for kinetics analysis. This model is lacking in data analysis packages associated with high throughput SPR instruments and the packages containing this model do not explore multiple local minima and parameter identifiability issues that are common in non-linear optimization. Therefore, we developed a method to use a system of ordinary differential equations for analyzing 1:2 binding kinetics data. Salient features of this method include a grid search on parameter initialization and a profile likelihood approach to determine parameter identifiability. Using this method we found a non-identifiable parameter in data set collected under the standard experimental design. A simulation-guided improved experimental design led to reliable estimation of all rate constants. The method and approach developed here for analyzing 1:2 binding kinetics data will be valuable for expeditious therapeutic antibody discovery research.

Focused goodness of fit tests for gene set analyses.

Gene set-based signal detection analyses are used to detect an association between a trait and a set of genes by accumulating signals across the genes in the gene set. Since signal detection is concerned with identifying whether any of the genes in the gene set are non-null, a goodness-of-fit (GOF) test can be used to compare whether the observed distribution of gene-level tests within the gene set agrees with the theoretical null distribution. Here, we present a flexible gene set-based signal detection framework based on tail-focused GOF statistics. We show that the power of the various statistics in this framework depends critically on two parameters: the proportion of genes within the gene set that are non-null and the degree of separation between the null and alternative distributions of the gene-level tests. We give guidance on which statistic to choose for a given situation and implement the methods in a fast and user-friendly R package, wHC (https://github.com/mqzhanglab/wHC). Finally, we apply these methods to a whole exome sequencing study of amyotrophic lateral sclerosis.

Bridging the gaps: using an NHP model to predict single dose radiation absorption in humans.

Purpose: Design and characterization of a radiation biodosimetry device are complicated by the fact that the requisite data are not available in the intended use population, namely humans exposed to a single, whole-body radiation dose. Instead, one must turn to model systems. We discuss our studies utilizing healthy, unexposed humans, human bone marrow transplant patients undergoing total body irradiation (TBI), non-human primates subjected to the same irradiation regimen received by the human TBI patients and NHPs given a single, whole-body dose of ionizing radiation.Materials and Methods: We use Bayesian linear mixed models to characterize the association between NHP and human expression patterns in radiation response genes when exposed to a common exposure regimen and across exposure regimens within the same species.Results: We show that population average differences in expression of our radiation response genes from one to another model system are comparable to typical differences between two randomly sampled members of a given model system and that these differences are smaller, on average, for linear combinations of the probe data and for the model-based combinations employed for dose prediction as part of a radiation biodosimetry device.Conclusions: Our analysis suggests that dose estimates based on our gene list will be accurate when applied to humans who have received a single, whole-body exposure to ionizing radiation.

Exploiting expression patterns across multiple tissues to map expression quantitative trait loci.

BACKGROUND: In order to better understand complex diseases, it is important to understand how genetic variation in the regulatory regions affects gene expression. Genetic variants found in these regulatory regions have been shown to activate transcription in a tissue-specific manner. Therefore, it is important to map the aforementioned expression quantitative trait loci (eQTL) using a statistically disciplined approach that jointly models all the tissues and makes use of all the information available to maximize the power of eQTL mapping. In this context, we are proposing a score test-based approach where we model tissue-specificity as a random effect and investigate an overall shift in the gene expression combined with tissue-specific effects due to genetic variants. RESULTS: Our approach has 1) a distinct computational edge, and 2) comparable performance in terms of statistical power over other currently existing joint modeling approaches such as MetaTissue eQTL and eQTL-BMA. Using simulations, we show that our method increases the power to detect eQTLs when compared to a tissue-by-tissue approach and can exceed the performance, in terms of computational speed, of MetaTissue eQTL and eQTL-BMA. We apply our method to two publicly available expression datasets from normal human brains, one comprised of four brain regions from 150 neuropathologically normal samples and another comprised of ten brain regions from 134 neuropathologically normal samples, and show that by using our method and jointly analyzing multiple brain regions, we identify eQTLs within more genes when compared to three often used existing methods. CONCLUSIONS: Since we employ a score test-based approach, there is no need for parameter estimation under the alternative hypothesis. As a result, model parameters only have to be estimated once per genome, significantly decreasing computation time. Our method also accommodates the analysis of next- generation sequencing data. As an example, by modeling gene transcripts in an analogous fashion to tissues in our current formulation one would be able to test for both a variant overall effect across all isoforms of a gene as well as transcript-specific effects. We implement our approach within the R package JAGUAR, which is now available at the Comprehensive R Archive Network repository.

Testing the effect of rare compound-heterozygous and recessive mutations in case--parent sequencing studies.

Compound heterozygous mutations are mutations that occur on different copies of genes and may completely "knock-out" gene function. Compound heterozygous mutations have been implicated in a large number of diseases, but there are few statistical methods for analyzing their role in disease, especially when such mutations are rare. A major barrier is that phase information is required to determine whether both gene copies are affected and phasing rare variants is difficult. Here, we propose a method to test compound heterozygous and recessive disease models in case-parent trios. We propose a simple algorithm for phasing and show via simulations that tests based on phased trios have almost the same power as tests using true phase information. A further complication in the study of compound heterozygous mutations is that only families where both parents carry mutations are informative. Thus, the informative sample size will be quite small even when the overall sample size is not, making asymptotic approximations of the null distribution of the test statistic inappropriate. To address this, we develop an exact test that will give appropriate P-values regardless of sample size. Using simulation, we show that our method is robust to population stratification and significantly outperforms other methods when the causal model is recessive.

Testing for risk and protective trends in genetic analyses of HIV acquisition.

Host genetics studies of HIV-1 acquisition are critically important for the identification of new targets for drug and vaccine development. Analyses of such studies typically focus on pairwise comparisons of three different groups: HIV-1 positive individuals, HIV-1 high-risk seronegative individuals, and population controls. Because there is a clear expectation of how gene frequencies of risk or protective alleles would be ordered in the three groups, we are able to construct a statistical framework that offers a consistent increase in power over a wide-range of the magnitude of risk/protective effects. In this paper, we develop tests that constrain the alternative hypothesis to appropriately reflect risk or protective trends jointly across the three groups and show that they lead to a substantial increase in power over the naive pairwise approach. We develop both likelihood-ratio and score statistics that test for genetic effects across the three groups while constraining the alternative hypothesis to reflect biologically motivated trends of risk or protection. The asymptotic distribution of both statistics (likelihood ratio and score) is derived. We investigate the performance of our approach via extensive simulation studies using a biologically motivated model of HIV-1 acquisition, and find that our proposed approach leads to an increase in power of roughly 10-28%. We illustrate our approach with an analysis of the effect of the CCR5Δ32 mutation on HIV acquisition.