A simple regression method for mapping quantitative trait loci in line crosses using flanking markers

The use of flanking marker methods has proved to be a powerful tool for the mapping of quantitative trait loci (QTL) in the segregating generations derived from crosses between inbred lines. Methods to analyse these data, based on maximum-likelihood, have been developed and provide good estimates of QTL effects in some situations. Maximum-likelihood methods are, however, relatively complex and can be computationally slow. In this paper we develop methods for mapping QTL based on multiple regression which can be applied using any general statistical package. We use the example of mapping in an F(2) population and show that these regression methods produce very similar results to those obtained using maximum likelihood. The relative simplicity of the regression methods means that models with more than a single QTL can be explored and we give examples of two lined loci and of two interacting loci. Other models, for example with more than two QTL, with environmental fixed effects, with between family variance or for threshold traits, could be fitted in a similar way. The ease, speed of application and generality of regression methods for flanking marker analysis, and the good estimates they obtain, suggest that they should provide the method of choice for the analysis of QTL mapping data from inbred line crosses.     

Unravelling the functional interaction structure of a cellular network from temporal slope information of experimental data

Due to the unavoidable nonbiological variations accompanying many experiments, it is imperative to consider a way of unravelling the functional interaction structure of a cellular network (e.g. signalling cascades or gene networks) by using the qualitative information of time-series experimental data instead of computation through the measured absolute values. In this spirit, we propose a very simple but effective method of identifying the functional interaction structure of a cellular network based on temporal ascending or descending slope information from given time-series measurements. From this method, we can gain insight into the acceptable measurement error ranges in order to estimate the correct functional interaction structure and we can also find guidance for a new experimental design to complement the insufficient information of a given experimental dataset. We developed experimental sign equations, making use of the temporal slope sign information from time-series experimental data, without a specific assumption on parameter perturbations for each network node. Based on these equations, we further describe the available specific information from each part of experimental data in detail and show the functional interaction structure obtained by integrating such information. In this procedure, we use only simple algebra on sign changes without complicated computations on the measured absolute values of the experimental data. The result is, however, verified through rigorous mathematical definitions and proofs. The present method provides us with information about the acceptable measurement error ranges for correct estimation of the functional interaction structure and it further leads to a new experimental design to complement the given experimental data by informing us about additional specific sampling points to be chosen for further required information.     

A SNP streak model for the identification of genetic regions identical-by-descent

The availability of very dense genetic maps is changing in a fundamental way the methods used to identify the genetic basis of both rare and common inherited traits. The ability to directly compare the genomes of two related individuals and quickly identify those regions that are inherited identical-by-descent (IBD) from a recent common ancestor would be of utility in a wide range of genetic mapping methods. Here, we describe a simple method for using dense SNP maps to identify regions of the genome likely to be inherited IBD by family members. This method is based on identifying obligate recombination events and examining the pattern of distribution of such events along the genetic map. Specifically, we use the length of a consecutive set of biallelic markers that have a high probability of having avoided such obligate recombination events. This ;;SNP streak" is derived from subsets of samples within a pedigree and allows us to make statistical inferences about the ancestry of the region(s) containing stretches of markers with these properties. We show that the use of subsets of more than two samples has the advantage of identifying shorter shared subsegments as significant. This mitigates the effects of errors in SNP calls. We provide specific examples of microarray-based SNP data, using a family with a complex pedigree and with a rare form of inherited kidney disease, to illustrate this approach.     

Analysing ecological networks of species interactions

Network approaches to ecological questions have been increasingly used, particularly in recent decades. The abstraction of ecological systems – such as communities – through networks of interactions between their components indeed provides a way to summarize this information with single objects. The methodological framework derived from graph theory also provides numerous approaches and measures to analyze these objects and can offer new perspectives on established ecological theories as well as tools to address new challenges. However, prior to using these methods to test ecological hypotheses, it is necessary that we understand, adapt, and use them in ways that both allow us to deliver their full potential and account for their limitations. Here, we attempt to increase the accessibility of network approaches by providing a review of the tools that have been developed so far, with – what we believe to be – their appropriate uses and potential limitations. This is not an exhaustive review of all methods and metrics, but rather, an overview of tools that are robust, informative, and ecologically sound. After providing a brief presentation of species interaction networks and how to build them in order to summarize ecological information of different types, we then classify methods and metrics by the types of ecological questions that they can be used to answer from global to local scales, including methods for hypothesis testing and future perspectives. Specifically, we show how the organization of species interactions in a community yields different network structures (e.g., more or less dense, modular or nested), how different measures can be used to describe and quantify these emerging structures, and how to compare communities based on these differences in structures. Within networks, we illustrate metrics that can be used to describe and compare the functional and dynamic roles of species based on their position in the network and the organization of their interactions as well as associated new methods to test the significance of these results. Lastly, we describe potential fruitful avenues for new methodological developments to address novel ecological questions.     

On Genetic Map Functions

Various genetic map functions have been proposed to infer the unobservable genetic distance between two loci from the observable recombination fraction between them. Some map functions were found to fit data better than others. When there are more than three markers, multilocus recombination probabilities cannot be uniquely determined by the defining property of map functions, and different methods have been proposed to permit the use of map functions to analyze multilocus data. If for a given map function, there is a probability model for recombination that can give rise to it, then joint recombination probabilities can be deduced from this model. This provides another way to use map functions in multilocus analysis. In this paper we show that stationary renewal processes give rise to most of the map functions in the literature. Furthermore, we show that the interevent distributions of these renewal processes can all be approximated quite well by gamma distributions.     

Two Statistical Methods for Resolving Healthy Individuals and Those with Congestive Heart Failure Based on Extended Self-similarity and a Recursive Method

In this paper we introduce two methods for measuring irregularities in human heartbeat time series (HHTS). First we consider the multi-fractal structure of HHTS to distinguish healthy individuals and from those with congestive heart failure. In this way we modify the Extended Self-Similarity (ESS) method and apply it to HHTS. Our second approach is based on the recursive method, which we use to predict the duration of the next heartbeat by considering a few previous ones. We use standard physiological data and show that these approaches lead to very satisfactory methods to resolve the healthy and CHF individuals. These methods can be used potentially in portable electronic heart alarm systems.     

Allele frequency estimation from ambiguous data: using resampling schema in validating frequency estimates and in selective neutrality testing

The development of molecular typing techniques applied to the study of population genetic diversity originates data with increasing precision but at the cost of some ambiguities. As distinct techniques may produce distinct kinds of ambiguities, a crucial issue is to assess the differences between frequency distributions estimated from data produced by alternative techniques for the same sample. To that aim, we developed a resampling scheme that allows evaluating, by statistical means, the significance of the difference between two frequency distributions. The same approach is then shown to be applicable to test selective neutrality when only sample frequencies are known. The use of these original methods is presented here through an application to the genetic study of a Munda human population sample, where three different HLA loci were typed using two different molecular methods (reverse PCR-SSO typing on microbeads arrays based on Luminex technology and PCR-SSP typing), as described in details in the companion article by Riccio et al. [The Austroasiatic Munda population from India and its enigmatic origin: An HLA diversity study. Hum. Biol. 38:405-435 (2011)]. The differences between the frequency estimates of the two typing techniques were found to be smaller than those resulting from sampling. Overall, we show that using a resampling scheme in validating frequency estimates is effective when alternative frequency estimates are available. Moreover, resampling appears to be the unique way to test selective neutrality when only frequency data are available to describe the genetic structure of populations.     

Modelling the interactions of protein side-chains

The study of small model molecules containing the relevant functional groups can help us to understand the interactions between side-chains in proteins.Ab initio quantum chemical techniques allow the interactions between the model molecules to be studied with much greater accuracy than is possible for an entire protein, where the use of simple empirical potentials is the norm. In particular, the use ofab initio methods on model molecules permits us to incorporate the atom-atom anisotropic directionality of these interactions. We survey various methods of obtaining the components of theab initio interaction energy. These are then applied to three systems of biological interest. The first of these is the arginine/aspartate pair found in salt bridges, which involves hydrogen bonding between two charged species. Secondly, we look at the arginine/phosphotyrosine interaction found in complexes between SH2 domains and peptide ligands: here we find that the arginine/phosphate part of the interaction is energetically far more important than the arginine/aromatic part. Finally, we describe a detailed study of amino/aromatic interactions in proteins: unconventional hydrogen bonds are found to be remarkably uncommon relative to stacked geometries, and the reasons for this are examined.     

A Method for Using Player Tracking Data in Basketball to Learn Player Skills and Predict Team Performance

Player tracking data represents a revolutionary new data source for basketball analysis, in which essentially every aspect of a player’s performance is tracked and can be analyzed numerically. We suggest a way by which this data set, when coupled with a network-style model of the offense that relates players’ skills to the team’s success at running different plays, can be used to automatically learn players’ skills and predict the performance of untested 5-man lineups in a way that accounts for the interaction between players’ respective skill sets. After developing a general analysis procedure, we present as an example a specific implementation of our method using a simplified network model. While player tracking data is not yet available in the public domain, we evaluate our model using simulated data and show that player skills can be accurately inferred by a simple statistical inference scheme. Finally, we use the model to analyze games from the 2011 playoff series between the Memphis Grizzlies and the Oklahoma City Thunder and we show that, even with a very limited data set, the model can consistently describe a player’s interactions with a given lineup based only on his performance with a different lineup.     

Rich probabilistic models for gene expression

Clustering is commonly used for analyzing gene expression data. Despite their successes, clustering methods suffer from a number of limitations. First, these methods reveal similarities that exist over all of the measurements, while obscuring relationships that exist over only a subset of the data. Second, clustering methods cannot readily incorporate additional types of information, such as clinical data or known attributes of genes. To circumvent these shortcomings, we propose the use of a single coherent probabilistic model, that encompasses much of the rich structure in the genomic expression data, while incorporating additional information such as experiment type, putative binding sites, or functional information. We show how this model can be learned from the data, allowing us to discover patterns in the data and dependencies between the gene expression patterns and additional attributes. The learned model reveals context-specific relationships, that exist only over a subset of the experiments in the dataset. We demonstrate the power of our approach on synthetic data and on two real-world gene expression data sets for yeast. For example, we demonstrate a novel functionality that falls naturally out of our framework: predicting the "cluster" of the array resulting from a gene mutation based only on the gene's expression pattern in the context of other mutations.     

The effect of multiple binding modes on empirical modeling of ligand docking to proteins.

A popular approach to the computational modeling of ligand/receptor interactions is to use an empirical free energy like model with adjustable parameters. Parameters are learned from one set of complexes, then used to predict another set. To improve these empirical methods requires an independent way to study their inherent errors. We introduce a toy model of ligand/receptor binding as a workbench for testing such errors. We study the errors incurred from the two state binding assumption--the assumption that a ligand is either bound in one orientation, or unbound. We find that the two state assumption can cause large errors in free energy predictions, but it does not affect rank order predictions significantly. We show that fitting parameters using data from high affinity ligands can reduce two state errors; so can using more physical models that do not use the two state assumption. We also find that when using two state models to predict free energies, errors are more severe on high affinity ligands than low affinity ligands. And we show that two state errors can be diagnosed by systematically adding new binding modes when predicting free energies: if predictions worsen as the modes are added, then the two state assumption in the fitting step may be at fault.     

A unified framework for unraveling the functional interaction structure of a biomolecular network based on stimulus-response experimental data

We propose a unified framework for the identification of functional interaction structures of biomolecular networks in a way that leads to a new experimental design procedure. In developing our approach, we have built upon previous work. Thus we begin by pointing out some of the restrictions associated with existing structure identification methods and point out how these restrictions may be eased. In particular, existing methods use specific forms of experimental algebraic equations with which to identify the functional interaction structure of a biomolecular network. In our work, we employ an extended form of these experimental algebraic equations which, while retaining their merits, also overcome some of their disadvantages. Experimental data are required in order to estimate the coefficients of the experimental algebraic equation set associated with the structure identification task. However, experimentalists are rarely provided with guidance on which parameters to perturb, and to what extent, to perturb them. When a model of network dynamics is required then there is also the vexed question of sample rate and sample time selection to be resolved. Supplying some answers to these questions is the main motivation of this paper. The approach is based on stationary and/or temporal data obtained from parameter perturbations, and unifies the previous approaches of Kholodenko et al. (PNAS 99 (2002) 12841-12846) and Sontag et al. (Bioinformatics 20 (2004) 1877-1886). By way of demonstration, we apply our unified approach to a network model which cannot be properly identified by existing methods. Finally, we propose an experiment design methodology, which is not limited by the amount of parameter perturbations, and illustrate its use with an in numero example.     

Incorporating spatial autocorrelation in rarefaction methods: Implications for ecologists and conservation biologists

Recently, methods for constructing Spatially Explicit Rarefaction (SER) curves have been introduced in the scientific literature to describe the relation between the recorded species richness and sampling effort and taking into account for the spatial autocorrelation in the data. Despite these methodological advances, the use of SERs has not become routine and ecologists continue to use rarefaction methods that are not spatially explicit. Using two study cases from Italian vegetation surveys, we demonstrate that classic rarefaction methods that do not account for spatial structure can produce inaccurate results. Furthermore, our goal in this paper is to demonstrate how SERs can overcome the problem of spatial autocorrelation in the analysis of plant or animal communities. Our analyses demonstrate that using a spatially-explicit method for constructing rarefaction curves can substantially alter estimates of relative species richness. For both analyzed data sets, we found that the rank ordering of standardized species richness estimates was reversed between the two methods. We strongly advise the use of Spatially Explicit Rarefaction methods when analyzing biodiversity: the inclusion of spatial autocorrelation into rarefaction analyses can substantially alter conclusions and change the way we might prioritize or manage nature reserves.     

The incidence of calorie labeling on fast food choices: A comparison between stated preferences and actual choices

In order to test the effect of calorie information on fast food choices, we conducted a questionnaire employing two types of stated preferences methods (the best-worst-scaling and intentional questions) and a follow-up randomized field experiment in a sample of 119 participants. This combined approach allowed us to test the internal validity of preferences for fast food meals across elicitation scenarios. The results showed that calorie information reduces the probability of selecting high calorie meals only in the questionnaire, while it did not have any significant impact on actual purchasing behavior in the field experiment. Thus, the findings show that there is a clear difference between the role of calorie information on immediate stated preference choices, and the relatively low level of responsiveness in real choices in a restaurant. We believe that the current results are quite suggestive, indicating the limits of predicting actual fast food behavior, and may open the way to using data sources that combine stated methods with field experiments.     

A Novel Scoring Approach for Protein Co-Purification Data Reveals High Interaction Specificity

Large-scale protein interaction networks (PINs) have typically been discerned using affinity purification followed by mass spectrometry (AP/MS) and yeast two-hybrid (Y2H) techniques. It is generally recognized that Y2H screens detect direct binary interactions while the AP/MS method captures co-complex associations; however, the latter technique is known to yield prevalent false positives arising from a number of effects, including abundance. We describe a novel approach to compute the propensity for two proteins to co-purify in an AP/MS data set, thereby allowing us to assess the detected level of interaction specificity by analyzing the corresponding distribution of interaction scores. We find that two recent AP/MS data sets of yeast contain enrichments of specific, or high-scoring, associations as compared to commensurate random profiles, and that curated, direct physical interactions in two prominent data bases have consistently high scores. Our scored interaction data sets are generally more comprehensive than those of previous studies when compared against four diverse, high-quality reference sets. Furthermore, we find that our scored data sets are more enriched with curated, direct physical associations than Y2H sets. A high-confidence protein interaction network (PIN) derived from the AP/MS data is revealed to be highly modular, and we show that this topology is not the result of misrepresenting indirect associations as direct interactions. In fact, we propose that the modularity in Y2H data sets may be underrepresented, as they contain indirect associations that are significantly enriched with false negatives. The AP/MS PIN is also found to contain significant assortative mixing; however, in line with a previous study we confirm that Y2H interaction data show weak disassortativeness, thus revealing more clearly the distinctive natures of the interaction detection methods. We expect that our scored yeast data sets are ideal for further biological discovery and that our scoring system will prove useful for other AP/MS data sets.     

RNA recognition by the joint action of two nucleolin RNA-binding domains: genetic analysis and structural modeling.

The interaction of nucleolin with a short stem-loop structure (NRE) requires two contiguous RNA-binding domains (RBD 1+2). The structural basis for RNA recognition by these RBDs was studied using a genetic system in Escherichia coli. Within each of the two domains, we identified several mutations that severely impair interaction with the RNA target. Mutations that alter RNA-binding specificity were also isolated, suggesting the identity of specific contacts between RBD 1+2 amino acids and nucleotides within the NRE stem-loop. Our data indicate that both RBDs participate in a joint interaction with the NRE and that each domain uses a different surface to contact the RNA. The constraints provided by these genetic data and previous mutational studies have enabled us to propose a three-dimensional model of nucleolin RBD 1+2 bound to the NRE stem-loop.     

Combining parametric and nonparametric models to estimate treatment effects in observational studies

Performing causal inference in observational studies requires we assume confounding variables are correctly adjusted for. In settings with few discrete-valued confounders, standard models can be employed. However, as the number of confounders increases these models become less feasible as there are fewer observations available for each unique combination of confounding variables. In this paper, we propose a new model for estimating treatment effects in observational studies that incorporates both parametric and nonparametric outcome models. By conceptually splitting the data, we can combine these models while maintaining a conjugate framework, allowing us to avoid the use of Markov chain Monte Carlo (MCMC) methods. Approximations using the central limit theorem and random sampling allow our method to be scaled to high-dimensional confounders. Through simulation studies we show our method can be competitive with benchmark models while maintaining efficient computation, and illustrate the method on a large epidemiological health survey.     

Computational studies of the extraction of visual spatial information from binocular and motion cues

This paper reviews some of the contributions that work in computational vision has made to the study of biological vision systems. We concentrate on two areas where there has been strong interaction between computational and experimental studies: the use of binocular stereo to recover the distances to surfaces in space, and the recovery of the three-dimensional shape of objects from relative motion in the image. With regard to stereo, we consider models proposed for solving the stereo correspondence problem, focussing on the way in which physical properties of the world constrain possible methods of solution. We also show how critical observations regarding human stereo vision have helped to shape these models. With regard to the recovery of structure from motion, we focus on how the constraint of object rigidity has been used in computational models of this process.