A model of immunological distances in systematics

While immunological distances among taxa have had wide use in systematics, there has been some doubt about their utility because of the observed non-metricity of such distance matrices. A model is presented here relating observed immunological distance to the actual number of antigenic site differences between taxa. This model accounts for the observed departures of these distances from the metric condition of reciprocity and triangle inequality. Based upon the model, two procedures are suggested for the transformation of immunological distances to metric distances appropriate for phylogenetic analysis. The model implies that the usual scaling adjustments applied to the immunological distance matrix are inappropriate; however, the same transformation applied instead to an initial similarity matrix will solve a scaling problem. Non-reciprocity of the distances is shown to remain a problem independent of this initial scaling problem. It is suggested that further transformation of these re-scaled distances may be obtained through an extension of the ADCLUS procedure developed in psychology. This approach suggests a general strategy for a transformation to metric distances, given a particular model of non-metricity for the data.     

Simplification and shift in cognition of political difference: applying the geometric modeling to the analysis of semantic similarity judgment

Perceiving differences by means of spatial analogies is intrinsic to human cognition. Multi-dimensional scaling (MDS) analysis based on Minkowski geometry has been used primarily on data on sensory similarity judgments, leaving judgments on abstractive differences unanalyzed. Indeed, analysts have failed to find appropriate experimental or real-life data in this regard. Our MDS analysis used survey data on political scientists' judgments of the similarities and differences between political positions expressed in terms of distance. Both distance smoothing and majorization techniques were applied to a three-way dataset of similarity judgments provided by at least seven experts on at least five parties' positions on at least seven policies (i.e., originally yielding 245 dimensions) to substantially reduce the risk of local minima. The analysis found two dimensions, which were sufficient for mapping differences, and fit the city-block dimensions better than the Euclidean metric in all datasets obtained from 13 countries. Most city-block dimensions were highly correlated with the simplified criterion (i.e., the left-right ideology) for differences that are actually used in real politics. The isometry of the city-block and dominance metrics in two-dimensional space carries further implications. More specifically, individuals may pay attention to two dimensions (if represented in the city-block metric) or focus on a single dimension (if represented in the dominance metric) when judging differences between the same objects. Switching between metrics may be expected to occur during cognitive processing as frequently as the apparent discontinuities and shifts in human attention that may underlie changing judgments in real situations occur. Consequently, the result has extended strong support for the validity of the geometric models to represent an important social cognition, i.e., the one of political differences, which is deeply rooted in human nature.     

Extended dissimilarity: a method of robust estimation of ecological distances from high beta diversity data

It is widely accepted that reliable ordination of ecological data requires a strong linear or ordinal relationship between the dissimilarity of sites, based on species composition, and the ecological distance between them. Certain dissimilarity measures, having the property that they take a fixed maximum value when sites have no species in common, have been shown to be strongly correlated with ecological distance. For ecological gradients of moderate length (moderate beta diversity), such measures, in conjunction with non-metric multidimensional scaling, will reliably yield successful ordinations. However, as beta diversity increases, more sites have no species in common, and such measures invariably under-estimate ecological distance for such sites. Thus ordinations of data with high species turnover (high beta diversity) may fail.Extended dissimilarities are defined using an iterative adaptation of flexible shortest path adjustment applied to the matrix of dissimilarities with fixed maximum values. By means of theoretical argument and simulations, this is shown to lead to far stronger correlations between the adjusted site dissimilarity and ecological distance for ecological gradients of greater length than previously considered. Hence ordinations of extended dissimilarities, by means of either metric or non-metric scaling techniques, are shown to outperform corresponding ordinations of unadjusted dissimilarities, with the difference increasing with increasing beta diversity.     

Sensitivity in metric scaling and analysis of distance

Assessing the sensitivity or sampling variability of multivariate ordination methods is essential if inferences are to be drawn from the analysis, but such assessment has to date been notably absent in many applications of multidimensional scaling (MDS). The only available technique seems to be the one by DeLeeuw and Meulman who proposed a special jackknife in a general MDS setting, but this method does not appear to have been widely used to date. A possible reason for this is that it is perceived to be computationally daunting. However, if attention is focused on classical metric scaling (principal coordinate analysis) then known analytical results can be used and the apparent computational complexity disappears. The purpose of this article is to set out these results, to indicate their use in more general analysis of distance, and to illustrate the methodology on some biometric examples.     

Approximate geodesic distances reveal biologically relevant structures in microarray data

MOTIVATION: Genome-wide gene expression measurements, as currently determined by the microarray technology, can be represented mathematically as points in a high-dimensional gene expression space. Genes interact with each other in regulatory networks, restricting the cellular gene expression profiles to a certain manifold, or surface, in gene expression space. To obtain knowledge about this manifold, various dimensionality reduction methods and distance metrics are used. For data points distributed on curved manifolds, a sensible distance measure would be the geodesic distance along the manifold. In this work, we examine whether an approximate geodesic distance measure captures biological similarities better than the traditionally used Euclidean distance. RESULTS: We computed approximate geodesic distances, determined by the Isomap algorithm, for one set of lymphoma and one set of lung cancer microarray samples. Compared with the ordinary Euclidean distance metric, this distance measure produced more instructive, biologically relevant, visualizations when applying multidimensional scaling. This suggests the Isomap algorithm as a promising tool for the interpretation of microarray data. Furthermore, the results demonstrate the benefit and importance of taking nonlinearities in gene expression data into account.     

Optimal Representation of Supplementary Variables in Biplots from Principal Component Analysis and Correspondence Analysis

This paper treats the topic of representing supplementary variables in biplots obtained by principal component analysis (PCA) and correspondence analysis (CA). We follow a geometrical approach where we minimize errors that are obtained when the scores of the PCA or CA solution are projected onto a vector that represents a supplementary variable. This paper shows that optimal directions for supplementary variables can be found by solving a regression problem, and justifies that earlier formulae from Gabriel are optimal in the least squares sense. We derive new results regarding the geometrical properties, goodness of fit statistics and the interpretation of supplementary variables. It is shown that supplementary variables can be represented by plotting their correlation coefficients with the axes of the biplot only when the proper type of scaling is used. We discuss supplementary variables in an ecological context and give illustrations with data from an environmental monitoring survey.     

Trajectory analysis of NMR structure calculations

Summary NMR as well as X-ray crystallography are used to determine the three-dimensional structures of macromolecules at atomic resolution. Structure calculation generates coordinates that are compatible with NMR data from randomly generated initial structures. We analyzed the trajectory taken by structures during NMR structure calculation in conformational space, assuming that the distance between two structures in conformational space is the root-mean-square deviation between the two structures. The coordinates of a structure in conformational space were obtained by applying the metric multidimensional scaling method. As an example, we used a 22-residue peptide, -Conotoxin GIIIA, and a simulated annealing protocol of XPLOR. We found that the three-dimensional solution of the multidimensional scaling analysis is sufficient to describe the overall configuration of the trajectories in conformational space. By comparing the trajectories of the entire calculation with those of the converged calculation, random sampling of conformational space is readily discernible. Trajectory analysis can also be used for optimization of protocols of NMR structure calculation, by examining individual trajectories.Abbreviations MD molecular dynamics - MDS multidimensional scaling - rmsd root-mean-square deviation - armsd angular rmsd - R multiple correlation coefficient - YASAP yet another simulated annealing protocol - PCA principal component analysis     

Multivariate data analysis in empirical research. A look on the bright side

The interpretive benefits of employing multivariate analysis methods on experimental data with more than one dependent variable are described heuristically and illustrated on a set of data from a simply designed experiment in physiological psychology. Multivariate analysis of variance (MANOVA) is performed on the 9 dependent variables contained in the sample data and on the four composites derived from a principal components analysis (PCA) of the variability of the nine. A linear discriminant analysis (LDA) is conducted following both MANOVA results, and 5 methods of determining the "important" dependent variables in the experimental-control group difference are presented and discussed in terms of the data at hand.     

Scaling properties of pyramidal neurons in mice neocortex

Dendritic morphology is the structural correlate for receiving and processing inputs to a neuron. An interesting question then is what the design principles and the functional consequences of enlarged or shrinked dendritic trees might be. As yet, only a few studies have examined the effects of neuron size changes. Two theoretical scaling modes have been analyzed, conservative (isoelectrotonic) scaling (preserves the passive and active response properties) and isometric scaling (steps up low pass-filtering of inputs). It has been suggested that both scaling modes were verified in neuroanatomical studies. To overcome obvious limitations of these studies like small size of analyzed samples and restricted validity of utilized scaling measures, we considered the scaling problem of neurons on the basis of large sample data and by employing a more general method of scaling analysis. This method consists in computing the morphoelectrotonic transform (MET) of neurons. The MET maps the neuron from anatomical space into electrotonic space using the logarithm of voltage attenuation as the distance metric. The theory underlying this approach is described and then applied to two samples of morphologically reconstructed pyramidal neurons (cells from neocortex of wildtype and synRas transgenic mice) using the NEURON simulator. In a previous study, we could verify a striking increase of dendritic tree size in synRas pyramidal neurons. Surprisingly, in this study the statistical analysis of the sample MET dendrograms revealed that the electrotonic architecture of these neurons scaled roughly in a MET-conserving mode. In conclusion, our results suggest only a minor impact of the Ras protein on dendritic electroanatomy, with non-significant changes of most regions of the corresponding METs.     

Oscillatory neural network for pattern recognition: trajectory based classification and supervised learning

Computer algorithms that match human performance in recognizing written text or spoken conversation remain elusive. The reasons why the human brain far exceeds any existing recognition scheme to date in the ability to generalize and to extract invariant characteristics relevant to category matching are not clear. However, it has been postulated that the dynamic distribution of brain activity (spatiotemporal activation patterns) is the mechanism by which stimuli are encoded and matched to categories. This research focuses on supervised learning using a trajectory based distance metric for category discrimination in an oscillatory neural network model. Classification is accomplished using a trajectory based distance metric. Since the distance metric is differentiable, a supervised learning algorithm based on gradient descent is demonstrated. Classification of spatiotemporal frequency transitions and their relation to a priori assessed categories is shown along with the improved classification results after supervised training. The results indicate that this spatiotemporal representation of stimuli and the associated distance metric is useful for simple pattern recognition tasks and that supervised learning improves classification results.     

Interpretation of ANOVA models for microarray data using PCA

MOTIVATION: ANOVA is a technique, which is frequently used in the analysis of microarray data, e.g. to assess the significance of treatment effects, and to select interesting genes based on P-values. However, it does not give information about what exactly is causing the effect. Our purpose is to improve the interpretation of the results from ANOVA on large microarray datasets, by applying PCA on the individual variance components. Interaction effects can be visualized by biplots, showing genes and variables in one plot, providing insight in the effect of e.g. treatment or time on gene expression. Because ANOVA has removed uninteresting sources of variance, the results are much more interpretable than without ANOVA. Moreover, the combination of ANOVA and PCA provides a simple way to select genes, based on the interactions of interest. RESULTS: It is shown that the components from an ANOVA model can be summarized and visualized with PCA, which improves the interpretability of the models. The method is applied to a real time-course gene expression dataset of mesenchymal stem cells. The dataset was designed to investigate the effect of different treatments on osteogenesis. The biplots generated with the algorithm give specific information about the effects of specific treatments on genes over time. These results are in agreement with the literature. The biological validation with GO annotation from the genes present in the selections shows that biologically relevant groups of genes are selected. AVAILABILITY: R code with the implementation of the method for this dataset is available from http://www.cac.science.ru.nl under the heading "Software".     

Studying spike trains using a van Rossum metric with a synapse-like filter

Spike trains are unreliable. For example, in the primary sensory areas, spike patterns and precise spike times will vary between responses to the same stimulus. Nonetheless, information about sensory inputs is communicated in the form of spike trains. A challenge in understanding spike trains is to assess the significance of individual spikes in encoding information. One approach is to define a spike train metric, allowing a distance to be calculated between pairs of spike trains. In a good metric, this distance will depend on the information the spike trains encode. This method has been used previously to calculate the timescale over which the precision of spike times is significant. Here, a new metric is constructed based on a simple model of synaptic conductances which includes binding site depletion. Including binding site depletion in the metric means that a given individual spike has a smaller effect on the distance if it occurs soon after other spikes. The metric proves effective at classifying neuronal responses by stimuli in the sample data set of electro-physiological recordings from the primary auditory area of the zebra finch fore-brain. This shows that this is an effective metric for these spike trains suggesting that in these spike trains the significance of a spike is modulated by its proximity to previous spikes. This modulation is a putative information-coding property of spike trains.     

Spatial Distribution of Algal Assemblages in a Temperate Lowland Peat Bog

Samples of phytobenthos were collected during three different seasons in 2005 along a linear transect of a lowland peat bog at various spatial scales (10 cm, 1 m, 10 m) to investigate the seasonal dynamics, diversity, and factors influencing the spatial patterns of microalgal communities. Non‐metric multidimensional scaling (NMDS), similarity percentage (SIMPER) analyses, ANOSIM, Mantel tests and diversity indices were used to analyze the data. Seasonal dynamics were exhibited by an increase in diversity, and a decrease in dominance from May to October, with significant differences in species composition. Mantel tests showed the significant influence of distance, microhabitat type, and conductivity on maintaining the similarity of species composition on scales of 1 m and 10 m. The small‐scale processes (colonization and niche differentiation), microhabitat type, geographic distance and conductivity were found to be the main factors influencing the distribution of algal assemblages. We conclude that these factors are related to winter disturbance, and the consequent colonization and subsequent niche differentiation. (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Compositional dissimilarity as a robust measure of ecological distance

The robustness of quantitative measures of compositional dissimilarity between sites was evaluated using extensive computer simulations of species' abundance patterns over one and two dimensional configurations of sample sites in ecological space. Robustness was equated with the strength over a range of models, of the linear and monotonic (rank-order) relationship between the compositional dissimilarities and the corresponding Euclidean distances between sites measured in the ecological space. The range of models reflected different assumptions about species' response curve shape, sampling pattern of sites, noise level of the data, species' interactions, trends in total site abundance, and beta diversity of gradients.The Kulczynski, Bray-Curtis and Relativized Manhattan measures were found to have not only a robust monotonic relationship with ecological distance, but also a robust linear (proportional) relationship until ecological distances became large. Less robust measures included Chord distance, Kendall's coefficient, Chisquared distance, Manhattan distance, and Euclidean distance.A new ordination method, hybrid multidimensional scaling (HMDS), is introduced that combines metric and nonmetric criteria, and so takes advantage of the particular properties of robust dissimilarity measures such as the Kulczynski measure.We thank M. P. Austin for encouraging this study, and I. C. Prentice, E. Van der Maarel, and an anonymous reviewer for helpful comments. E. M. Adomeit provided technical assistance.     

Limited sampling of conformational space by the distance geometry algorithm: implications for structures generated from NMR data

Calculations with a metric matrix distance geometry algorithm were performed that show that the standard implementation of the algorithm generally samples a very limited region of conformational space. This problem is most severe when only a small amount of distance information is used as input for the algorithm. Control calculations were performed on linear peptides, disulfide-linked peptides, and a double-stranded DNA decamer where only distances defining the covalent structures of the molecules (as well as the hydrogen bonds for the base pairs in the DNA) were included as input. Since the distance geometry algorithm is commonly used to generate structures of biopolymers from distance data obtained from NMR experiments, simulations were performed on the small globular protein basic pancreatic trypsin inhibitor (BPTI) that mimic calculations performed with actual NMR data. The results on BPTI and on the control peptides indicate that the standard implementation of the algorithm has two main problems: first, that it generates extended structures; second, that it has a tendency to consistently produce similar structures instead of sampling all structures consistent with the input distance information. These results also show that use of a simple root-mean-square deviation for evaluating the quality of the structures generated from NMR data may not be generally appropriate. The main sources of these problems are identified, and our results indicate that the problems are not a fundamental property of the distance geometry algorithm but arise from the implementations presently used to generate structures from NMR data. Several possible methods for alleviating these problems are discussed.     

A Scaling Rule for Landscape Patches and How It Applies to Conserving Soil Resources in Savannas

Scaling issues are complex, yet understanding issues such as scale dependencies in ecological patterns and processes is usually critical if we are to make sense of ecological data and if we want to predict how land management options, for example, are constrained by scale. In this article, we develop the beginnings of a way to approach the complexity of scaling issues. Our approach is rooted in scaling functions, which integrate the scale dependency of patterns and processes in landscapes with the ways that organisms scale their responses to these patterns and processes. We propose that such functions may have sufficient generality that we can develop scaling rules—statements that link scale with consequences for certain phenomena in certain systems. As an example, we propose that in savanna ecosystems, there is a consistent relationship between the size of vegetation patches in the landscape and the degree to which critical resources, such as soil nutrients or water, become concentrated in these patches. In this case, the features of the scaling functions that underlie this rule have to do with physical processes, such as surface water flow and material redistribution, and the ways that patches of plants physically “capture” such runoff and convert it into plant biomass, thereby concentrating resources and increasing patch size. To be operationally useful, such scaling rules must be expressed in ways that can generate predictions. We developed a scaling equation that can be used to evaluate the potential impacts of different disturbances on vegetation patches and on how soils and their nutrients are conserved within Australian savanna landscapes. We illustrate that for a 10-km² paddock, given an equivalent area of impact, the thinning of large tree islands potentially can cause a far greater loss of soil nitrogen (21 metric tons) than grazing out small grass clumps (2 metric tons). Although our example is hypothetical, we believe that addressing scaling problems by first conceptualizing scaling functions, then proposing scaling rules, and then deriving scaling equations is a useful approach. Scaling equations can be used in simulation models, or (as we have done) in simple hypothetical scenarios, to collapse the complexity of scaling issues into a manageable framework. Received 8 December 1998; accepted 17 August 1999.     

Multivariate analysis of benthic invertebrate communities: the implication of choosing particular data standardizations,measures of association,and ordination methods

Benthic invertebrate data from thirty-nine lakes in south-central Ontario were analyzed to determine the effect of choosing particular data standardizations, resemblance measures, and ordination methods on the resultant multivariate summaries. Logarithmic-transformed, 0–1 scaled, and ranked data were used as standardized variables with resemblance measures of Bray-Curtis, Euclidean distance, cosine distance, correlation, covariance and chi-squared distance. Combinations of these measures and standardizations were used in principal components analysis, principal coordinates analysis, non-metric multidimensional scaling, correspondence analysis, and detrended correspondence analysis. Correspondence analysis and principal components analysis using a correlation coefficient provided the most consistent results irrespective of the choice in data standardization. Other approaches using detrended correspondence analysis, principal components analysis, principal coordinates analysis, and non-metric multidimensional scaling provided less consistent results. These latter three methods produced similar results when the abundance data were replaced with ranks or standardized to a 0–1 range. The log-transformed data produced the least consistent results, whereas ranked data were most consistent. Resemblance measures such as the Bray-Curtis and correlation coefficient provided more consistent solutions than measures such as Euclidean distance or the covariance matrix when different data standardizations were used. The cosine distance based on standardized data provided results comparable to the CA and DCA solutions. Overall, CA proved most robust as it demonstrated high consistency irrespective of the data standardizations. The strong influence of data standardization on the other ordination methods emphasizes the importance of this frequently neglected stage of data analysis.     

Projecting human pharmacokinetics of therapeutic antibodies from nonclinical data

The pharmacokinetics (PK) of therapeutic antibodies is determined by target and non-target mediated mechanisms. These antibody-specific factors need to be considered during prediction of human PK based upon preclinical information. Principles of allometric scaling established for small molecules using data from multiple animal species cannot be directly applied to antibodies. Here, different methods for projecting human clearance (CL) from animal PK data for 13 therapeutic monoclonal antibodies (mAbs) exhibiting linear PK over the tested dose ranges were examined: simple allometric scaling (CL versus body weight), allometric scaling with correction factors, allometric scaling based on rule of exponent and scaling from only cynomolgus monkey PK data. A better correlation was obtained between the observed human CL and the estimated human CL based on cynomolgus monkey PK data and an allometric scaling exponent of 0.85 for CL than other scaling approaches. Human concentration-time profiles were also reasonably predicted from the cynomolgus monkey data using species-invariant time method with a fixed exponent of 0.85 for CL and 1.0 for volume of distribution. In conclusion, we expanded our previous work and others and further confirmed that PK from cynomolgus monkey alone can be successfully scaled to project human PK profiles within linear range using simplify allometry and Dedrick plots with fixed exponent.