Search

We describe a general approach to transforming molecular models between different levels of resolution, based on machine learning methods. The approach uses a matched set of models at both levels of resolution for training, but requires only the coordinates of their particles and no side information (e.g., templates for substructures, defined mappings, or molecular mechanics force fields). Once trained, the approach can transform further molecular models of the system between the two levels of resolution in either direction with equal facility.

Human mesenchymal stem cells (hMSCs) are widely represented in regenerative medicine clinical strategies due to their compatibility with autologous implantation. Effective bone regeneration involves crosstalk between macrophages and hMSCs, with macrophages playing a key role in the recruitment and differentiation of hMSCs. However, engineered biomaterials able to simultaneously direct hMSC fate and modulate macrophage phenotype have not yet been identified. A novel combinatorial chemistry-topography screening platform, the ChemoTopoChip, is used here to identify materials suitable for bone regeneration by screening 1008 combinations in each experiment for human immortalized mesenchymal stem cell (hiMSCs) and human macrophage response. The osteoinduction achieved in hiMSCs cultured on the "hit" materials in basal media is comparable to that seen when cells are cultured in osteogenic media, illustrating that these materials offer a materials-induced alternative to osteo-inductive supplements in bone-regeneration. Some of these same chemistry-microtopography combinations also exhibit immunomodulatory stimuli, polarizing macrophages towards a pro-healing phenotype. Maximum control of cell response is achieved when both chemistry and topography are recruited to instruct the required cell phenotype, combining synergistically. The large combinatorial library allows us for the first time to probe the relative cell-instructive roles of microtopography and material chemistry which we find to provide similar ranges of cell modulation for both cues. Machine learning is used to generate structure-activity relationships that identify key chemical and topographical features enhancing the response of both cell types, providing a basis for a better understanding of cell response to micro topographically patterned polymers., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

In a previous contribution, a new Riemannian shape space, named TPS space, was introduced to perform statistics on shape data. This space was endowed with a Riemannian metric and a flat connection, with torsion, compatible with the given metric. This connection allows the definition of a Parallel Transport of the deformation compatible with the three-fold decomposition in spherical, deviatoric, and non-affine components. Such a parallel transport also conserves the Γ-energy, strictly related to the total elastic strain energy stored by the body in the original deformation. A new approach is here presented in order to calculate the bending energy on the body alone (body bending energy) and to restrict it exclusively within physical boundaries of objects involved in the deformation analysis. The novelty of this new procedure resides in the fact that we propose a new metric to be preserved during the TPS direct transport. This allows transporting the shape change more coherently with the mechanical meaning of the deformation. The geometry of the TPS space is then further discussed in order to better represent the relationship between the Γ-energy, the strain energy, and the so-called bending energy densities., (© 2019 John Wiley & Sons, Ltd.)

Solid dispersions can be a successful way to enhance the bioavailability of poorly soluble drugs. Here 60 solid dispersion formulations were produced using ten chemically diverse, neutral, poorly soluble drugs, three commonly used polymers, and two manufacturing techniques, spray-drying and melt extrusion. Each formulation underwent a six-month stability study at accelerated conditions, 40 °C and 75% relative humidity (RH). Significant differences in times to crystallization (onset of crystallization) were observed between both the different polymers and the two processing methods. Stability from zero days to over one year was observed. The extensive experimental data set obtained from this stability study was used to build multiple linear regression models to correlate physicochemical properties of the active pharmaceutical ingredients (API) with the stability data. The purpose of these models is to indicate which combination of processing method and polymer carrier is most likely to give a stable solid dispersion. Six quantitative mathematical multiple linear regression-based models were produced based on selection of the most influential independent physical and chemical parameters from a set of 33 possible factors, one model for each combination of polymer and processing method, with good predictability of stability. Three general rules are proposed from these models for the formulation development of suitably stable solid dispersions. Namely, increased stability is correlated with increased glass transition temperature ( T g ) of solid dispersions, as well as decreased number of H-bond donors and increased molecular flexibility (such as rotatable bonds and ring count) of the drug molecule.

In landmarks-based Shape Analysis size is measured, in most cases, with Centroid Size. Changes in shape are decomposed in affine and non affine components. Furthermore the non affine component can be in turn decomposed in a series of local deformations (partial warps). If the extent of deformation between two shapes is small, the difference between Centroid Size and m-Volume increment is barely appreciable. In medical imaging applied to soft tissues bodies can undergo very large deformations, involving large changes in size. The cardiac example, analyzed in the present paper, shows changes in m-Volume that can reach the 60%. We show here that standard Geometric Morphometrics tools (landmarks, Thin Plate Spline, and related decomposition of the deformation) can be generalized to better describe the very large deformations of biological tissues, without losing a synthetic description. In particular, the classical decomposition of the space tangent to the shape space in affine and non affine components is enriched to include also the change in size, in order to give a complete description of the tangent space to the size-and-shape space. The proposed generalization is formulated by means of a new Riemannian metric describing the change in size as change in m-Volume rather than change in Centroid Size. This leads to a redefinition of some aspects of the Kendall's size-and-shape space without losing Kendall's original formulation. This new formulation is discussed by means of simulated examples using 2D and 3D platonic shapes as well as a real example from clinical 3D echocardiographic data. We demonstrate that our decomposition based approaches discriminate very effectively healthy subjects from patients affected by Hypertrophic Cardiomyopathy., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Most functional MRI (fMRI) studies map task-driven brain activity using a block or event-related paradigm. Sparse paradigm free mapping (SPFM) can detect the onset and spatial distribution of BOLD events in the brain without prior timing information, but relating the detected events to brain function remains a challenge. In this study, we developed a decoding method for SPFM using a coordinate-based meta-analysis method of activation likelihood estimation (ALE). We defined meta-maps of statistically significant ALE values that correspond to types of events and calculated a summation overlap between the normalized meta-maps and SPFM maps. As a proof of concept, this framework was applied to relate SPFM-detected events in the sensorimotor network (SMN) to six motor functions (left/right fingers, left/right toes, swallowing, and eye blinks). We validated the framework using simultaneous electromyography (EMG)-fMRI experiments and motor tasks with short and long duration, and random interstimulus interval. The decoding scores were considerably lower for eye movements relative to other movement types tested. The average successful rate for short and long motor events were 77 ± 13% and 74 ± 16%, respectively, excluding eye movements. We found good agreement between the decoding results and EMG for most events and subjects, with a range in sensitivity between 55% and 100%, excluding eye movements. The proposed method was then used to classify the movement types of spontaneous single-trial events in the SMN during resting state, which produced an average successful rate of 22 ± 12%. Finally, this article discusses methodological implications and improvements to increase the decoding performance. Hum Brain Mapp 38:5778-5794, 2017. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

Root architecture impacts water and nutrient uptake efficiency. Identifying exactly which root architectural properties influence these agronomic traits can prove challenging. In this paper, approximately 300 wheat (Triticum aestivum) plants were divided into four groups using two binary classifications, high versus low nitrogen uptake efficiency (NUpE), and high versus low nitrate in the growth medium. The root system architecture for each wheat plant was captured using 16 quantitative variables. The multivariate analysis tool, linear discriminant analysis, was used to construct composite variables, each a linear combination of the original variables, such that the score of the plants on the new variables showed the maximum between-group variability. The results show that the distribution of root-system architecture traits differs between low- and high-NUpE plants and, less strongly, between low-NUpE plants grown on low versus high nitrate media., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Experimental Biology.)

We aimed to identify novel molecular mechanisms for muscle growth during administration of anabolic agents. Growing pigs (Duroc/(Landrace/Large-White)) were administered Ractopamine (a beta-adrenergic agonist; BA; 20 ppm in feed) or Reporcin (recombinant growth hormone; GH; 10 mg/48 hours injected) and compared to a control cohort (feed only; no injections) over a 27-day time course (1, 3, 7, 13 or 27-days). Longissimus Dorsi muscle gene expression was analyzed using Agilent porcine transcriptome microarrays and clusters of genes displaying similar expression profiles were identified using a modified maSigPro clustering algorithm. Anabolic agents increased carcass (p = 0.002) and muscle weights (Vastus Lateralis: p < 0.001; Semitendinosus: p = 0.075). Skeletal muscle mRNA expression of serine/one-carbon/glycine biosynthesis pathway genes (Phgdh, Psat1 and Psph) and the gluconeogenic enzyme, phosphoenolpyruvate carboxykinase-M (Pck2/PEPCK-M), increased during treatment with BA, and to a lesser extent GH (p < 0.001, treatment x time interaction). Treatment with BA, but not GH, caused a 2-fold increase in phosphoglycerate dehydrogenase (PHGDH) protein expression at days 3 (p < 0.05) and 7 (p < 0.01), and a 2-fold increase in PEPCK-M protein expression at day 7 (p < 0.01). BA treated pigs exhibit a profound increase in expression of PHGDH and PEPCK-M in skeletal muscle, implicating a role for biosynthetic metabolic pathways in muscle growth.

We consider the problem of comparing sizes and shapes of objects when landmark data are prone to measurement error. We show that naive implementation of ordinary Procrustes analysis that ignores measurement error can compromise inference. To account for measurement error, we propose the conditional score method for matching configurations, which guarantees consistent inference under mild model assumptions. The effects of measurement error on inference from naive Procrustes analysis and the performance of the proposed method are illustrated via simulation and application in three real data examples. Supplementary materials for this article are available online.

This is a discussion of the paper: 'Overview of object oriented data analysis' by J. Steve Marron and Andrés M. Alonso., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

fMRI studies of brain activity at rest study slow (<0.1 Hz) intrinsic fluctuations in the blood-oxygenation-level-dependent (BOLD) signal that are observed in a temporal scale of several minutes. The origin of these fluctuations is not clear but has previously been associated with slow changes in rhythmic neuronal activity resulting from changes in cortical excitability or neuronal synchronization. In this work, we show that individual spontaneous BOLD events occur during rest, in addition to slow fluctuations. Individual spontaneous BOLD events were identified by deconvolving the hemodynamic impulse response function for each time point in the fMRI time series, thus requiring no information on timing or a-priori spatial information of events. The patterns of activation detected were related to the motor, visual, default-mode, and dorsal attention networks. The correspondence between spontaneous events and slow fluctuations in these networks was assessed using a sliding window, seed-correlation analysis, where seed regions were selected based on the individual spontaneous event BOLD activity maps. We showed that the correlation varied considerably over time, peaking at the time of spontaneous events in these networks. By regressing spontaneous events out of the fMRI signal, we showed that both the correlation strength and the power in spectral frequencies <0.1 Hz decreased, indicating that spontaneous activation events contribute to low-frequency fluctuations observed in resting state networks with fMRI. This work provides new insights into the origin of signals detected in fMRI studies of functional connectivity., (Copyright © 2011 Wiley Periodicals, Inc.)

The ability to detect single trial responses in functional magnetic resonance imaging (fMRI) studies is essential, particularly if investigating learning or adaptation processes or unpredictable events. We recently introduced paradigm free mapping (PFM), an analysis method that detects single trial blood oxygenation level dependent (BOLD) responses without specifying prior information on the timing of the events. PFM is based on the deconvolution of the fMRI signal using a linear hemodynamic convolution model. Our previous PFM method (Caballero-Gaudes et al., 2011: Hum Brain Mapp) used the ridge regression estimator for signal deconvolution and required a baseline signal period for statistical inference. In this work, we investigate the application of sparse regression techniques in PFM. In particular, a novel PFM approach is developed using the Dantzig selector estimator, solved via an efficient homotopy procedure, along with statistical model selection criteria. Simulation results demonstrated that, using the Bayesian information criterion to select the regularization parameter, this method obtains high detection rates of the BOLD responses, comparable with a model-based analysis, but requiring no information on the timing of the events and being robust against hemodynamic response function variability. The practical operation of this sparse PFM method was assessed with single-trial fMRI data acquired at 7T, where it automatically detected all task-related events, and was an improvement on our previous PFM method, as it does not require the definition of a baseline state and amplitude thresholding and does not compromise on specificity and sensitivity., (Copyright © 2011 Wiley Periodicals, Inc.)

The two key aims of this research were: (i) to conduct a probabilistic elicitation to quantify the variation in veterinarians' beliefs regarding the efficacy of systemic antibiotics when used as an adjunct to intra-mammary dry cow therapy and (ii) to investigate (in a Bayesian statistical framework) the strength of future research evidence required (in theory) to change the beliefs of practising veterinary surgeons regarding the efficacy of systemic antibiotics, given their current clinical beliefs. The beliefs of 24 veterinarians in 5 practices in England were quantified as probability density functions. Classic multidimensional scaling revealed major variations in beliefs both within and between veterinary practices which included: confident optimism, confident pessimism and considerable uncertainty. Of the 9 veterinarians interviewed holding further cattle qualifications, 6 shared a confidently pessimistic belief in the efficacy of systemic therapy and whilst 2 were more optimistic, they were also more uncertain. A Bayesian model based on a synthetic dataset from a randomised clinical trial (showing no benefit with systemic therapy) predicted how each of the 24 veterinarians' prior beliefs would alter as the size of the clinical trial increased, assuming that practitioners would update their beliefs rationally in accordance with Bayes' theorem. The study demonstrated the usefulness of probabilistic elicitation for evaluating the diversity and strength of practitioners' beliefs. The major variation in beliefs observed raises interest in the veterinary profession's approach to prescribing essential medicines. Results illustrate the importance of eliciting prior beliefs when designing clinical trials in order to increase the chance that trial data are of sufficient strength to alter the clinical beliefs of practitioners and do not merely serve to satisfy researchers., (Copyright © 2012 Elsevier B.V. All rights reserved.)

A general framework for a novel non-geodesic decomposition of high-dimensional spheres or high-dimensional shape spaces for planar landmarks is discussed. The decomposition, principal nested spheres, leads to a sequence of submanifolds with decreasing intrinsic dimensions, which can be interpreted as an analogue of principal component analysis. In a number of real datasets, an apparent one-dimensional mode of variation curving through more than one geodesic component is captured in the one-dimensional component of principal nested spheres. While analysis of principal nested spheres provides an intuitive and flexible decomposition of the high-dimensional sphere, an interesting special case of the analysis results in finding principal geodesics, similar to those from previous approaches to manifold principal component analysis. An adaptation of our method to Kendall's shape space is discussed, and a computational algorithm for fitting principal nested spheres is proposed. The result provides a coordinate system to visualize the data structure and an intuitive summary of principal modes of variation, as exemplified by several datasets.

The placenta is an essential component of the life-support system for the developing foetus, enabling nutrients and waste to be exchanged between the foetal and maternal circulations. Maternal blood flows between the densely packed branches of villous trees, within which are foetal vessels. Here, we explore some of the challenges in modelling maternal haemodynamic transport using homogenization approaches. We first show how two measures can be used to estimate the minimum distance over which the distribution of villous branches appears statistically homogeneous. We then analyse a simplified model problem (solute transport by a unidirectional flow past a distribution of point sinks) to assess the accuracy of homogenization approximations as a function of governing parameters (Péclet and Damköhler numbers) and the statistical properties of the sink distribution. The difference between the leading-order homogenization approximation and the exact solute distribution is characterized by large spatial gradients at the scale of individual villi and substantial fluctuations that can be correlated over lengthscales comparable to the whole domain. This study highlights the importance of quantifying errors owing to spatial disorder in multi-scale approximations of physiological systems.

This work presents a novel method of mapping the brain's response to single stimuli in space and time without prior knowledge of the paradigm timing: paradigm free mapping (PFM). This method is based on deconvolution of the hemodynamic response from the voxel time series assuming a linear response and using a ridge-regression algorithm. Statistical inference is performed by defining a spatio-temporal t-statistic and by controlling for multiple comparisons using the false discovery rate procedure. The methodology was validated on five subjects who performed self-paced and visually cued finger tapping at 7 Tesla, with moderate (TR = 2 s) and high (TR = 0.4 s) temporal resolution. The results demonstrate that detection of single-trial BOLD events is feasible without a priori information on the stimulus paradigm. The proposed method opens up the possibility of designing temporally unconstrained paradigms to study the cortical response to unpredictable mental events., (Copyright © 2010 Wiley-Liss, Inc.)

Some methods for the statistical analysis of surface shapes and asymmetry are introduced. We focus on a case study where magnetic resonance images of the brain are available from groups of 30 schizophrenia patients and 38 controls, and we investigate large-scale brain surface shape differences. Key aspects of shape analysis are to remove nuisance transformations by registration and to identify which parts of one object correspond with the parts of another object. We introduce maximum likelihood and Bayesian methods for registering brain images and providing large-scale correspondences of the brain surfaces. Brain surface size-and-shape analysis is considered using random field theory, and also dimension reduction is carried out using principal and independent components analysis. Some small but significant differences are observed between the the patient and control groups. We then investigate a particular type of asymmetry called torque. Differences in asymmetry are observed between the control and patient groups, which add strength to other observations in the literature. Further investigations of the midline plane location in the 2 groups and the fitting of nonplanar curved midlines are also considered.

Milk recording data collected in 2,128 dairy herds in England and Wales between 2004 and 2006 were used to predict the calving to conception intervals. The average cumulative milk production was 8,200 kg. Conception (or not) within 5 intervals measured in days (20 to 60 d, 61 to 81 d, 82 to 102 d, 103 to 123 d, 124 to 144 d) was modeled as a function of milk yields and milk constituents at the start of lactation using multilevel discrete-time survival models. Milk yield, weight and percentage of fat, protein, and lactose, and somatic cell counts on the first 2 test-days of lactation were corrected for either stage of lactation alone or stage of lactation and time of year. Five hundred and 1,628 herds, respectively, were used for parameter estimation and cross-validation. Covariates were retained in the final model if their coefficient was at least twice its standard error and their inclusion resulted in a decrease in the deviance. Overall, 73% of cows recalved. The percentage of cows that had conceived by d 20, 61, 82, 103, 124, and 145 were 0.5, 7.3, 17.9, 29.3, 38.7, and 46, respectively. The probability of conception before 145 d in milk increased with lower milk production on the second test-day, higher percentage of protein on the second test-day, and higher percentage of lactose on the first test-day. Positive associations were of a limited magnitude but nonetheless significant with the percentage of protein on the first test-day, the percentage of butterfat on the first test-day, and somatic cell count on both test-days. The model predicted the probability of conception in the cross validation data set very well. Despite the common use of fat to protein ratio as a measure of energy balance, this parameter exhibited wide variation with stage of lactation and time of the year and had a much-reduced ability to predict an early conception compared with other combinations of milk quantity and constituents., (Copyright © 2010 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Background: Microarrays are a powerful tool used for the determination of global RNA expression. There is an increasing requirement to focus on profiling gene expression in tissues where it is difficult to obtain large quantities of material, for example individual tissues within organs such as the root, or individual isolated cells. From such samples, it is difficult to produce the amount of RNA required for labelling and hybridisation in microarray experiments, thus a process of amplification is usually adopted. Despite the increasing use of two-cycle amplification for transcriptomic analyses on the Affymetrix ATH1 array, there has been no report investigating any potential bias in gene representation that may occur as a result., Results: Here we compare transcriptomic data generated using Affymetrix one-cycle (standard labelling protocol), two-cycle (small-sample protocol) and IVT-Express protocols with the Affymetrix ATH1 array using Arabidopsis root samples. Results obtained with each protocol are broadly similar. However, we show that there are 35 probe sets (of a total of 22810) that are misrepresented in the two-cycle data sets. Of these, 33 probe sets were classed as mis-amplified when comparisons of two independent publicly available data sets were undertaken., Conclusions: Given the unreliable nature of the highlighted probes, we caution against using data associated with the corresponding genes in analyses involving transcriptomic data generated with two-cycle amplification protocols. We have shown that the Affymetrix IVT-E labelling protocol produces data with less associated bias than the two-cycle protocol, and as such, would recommend this kit for new experiments that involve small samples.

Molecular-dynamics simulations of a normal DNA duplex show that breathing events typically occur on the microsecond time scale. This paper analyzes a 12 base pairs DNA duplex containing the "rogue" base difluorotoluene (F) in place of a thymine base (T), for which the breathing events occur on the nanosecond time scale. Starting from a nonlinear Klein-Gordon lattice model and adding noise and damping, we obtain a mesoscopic model of the DNA duplex close to that observed in experiments and all-atom molecular dynamics simulations. The mesoscopic model is calibrated to data from the all-atom molecular dynamics package AMBER for a variety of twist angles of the DNA duplex. Defects are considered in the interchain interactions as well as in the along-chain interactions. This paper also discusses the role of the fluctuation-dissipation relations in the derivation of reduced (mesoscopic) models, the differences between the potential of mean force and the potential energies used in Klein-Gordon lattices, and how breathing can be viewed as competition between the along-chain elastic energy and the interchain binding energy.

Familial correlations in the spinal shape of children are explored using statistical shape analysis. Measurements of the spine on the surface of the back were taken using an optical device, and recorded for several hundred children from the Leeds area of the U.K. A portion of the spine from the distal thoracic to proximal lumbar spine was used for the analysis, and measures of the shape of the spine line were determined using Procrustes analysis and principal components (PCs) analysis. Analysis was carried out on scans repeated in triplicate and scans taken six months apart, and good repeatability was demonstrated. Four groups of children were considered in the main study: monozygotic twins, dizygotic twins, same sex siblings and different sex siblings. Intra-familial correlations were calculated for the shape measures, and also for age and sex matched unrelated controls. The second PC score (which measures possible lordosis) had significant positive correlation in the family groups. It was observed that the familial correlations were higher for the genetically closer groups. Also, the same sex siblings had higher correlation than the different sex siblings. All sibling groups showed greater correlation of sagittal profile for the second PC score than unrelated controls. A significant correlation was observed in same sex pairings for the control data. This work suggests that some elements of spinal profile may be familial but also shows correlation with sex. Both of these observations may be important in the aetiology of idiopathic scoliosis., (Copyright (c) 2007 John Wiley & Sons, Ltd.)

We consider Bayesian methodology for comparing two or more unlabeled point sets. Application of the technique to a set of steroid molecules illustrates its potential utility involving the comparison of molecules in chemoinformatics and bioinformatics. We initially match a pair of molecules, where one molecule is regarded as random and the other fixed. A type of mixture model is proposed for the point set coordinates, and the parameters of the distribution are a labeling matrix (indicating which pairs of points match) and a concentration parameter. An important property of the likelihood is that it is invariant under rotations and translations of the data. Bayesian inference for the parameters is carried out using Markov chain Monte Carlo simulation, and it is demonstrated that the procedure works well on the steroid data. The posterior distribution is difficult to simulate from, due to multiple local modes, and we also use additional data (partial charges on atoms) to help with this task. An approximation is considered for speeding up the simulation algorithm, and the approximating fast algorithm leads to essentially identical inference to that under the exact method for our data. Extensions to multiple molecule alignment are also introduced, and an algorithm is described which also works well on the steroid data set. After all the steroid molecules have been matched, exploratory data analysis is carried out to examine which molecules are similar. Also, further Bayesian inference for the multiple alignment problem is considered.

We describe the application of statistical shape analysis to homologous landmarks on the cortical surface of the adult human brain. Statistical shape analysis has a sound theoretical basis. Landmarks are identified on the surface of a 3-D reconstruction of the segmented cortical surface from magnetic resonance image (MRI) data. Using publicly available software (morphologika) the location and size dependence of the landmarks are removed and the differences in landmark distribution across subjects are analysed using principal component analysis. These differences, representing shape differences between subjects, can be visually assessed using wireframe models and transformation grids. The MRI data of 58 adult brains (27 female and 15 left handed) were examined. Shape differences in the whole brain are described which concern the relative orientation of frontal lobe sulci. Analysis of all 116 hemispheres revealed a statistically significant difference (P < 0.001) between left and right hemispheres. This finding was significant for right- but not left-handed subjects alone. No other significant age, gender, handedness, or brain-size correlations with shape differences were found., (Copyright 2001 Academic Press.)

We propose modeling a nearly regular point pattern by a generalized Neyman-Scott process in which the offspring are Gaussian perturbations from a regular mean configuration. The mean configuration of interest is an equilateral grid, but our results can be used for any stationary regular grid. The case of uniformly distributed points is first studied as a benchmark. By considering the square of the interpoint distances, we can evaluate the first two moments of the K-function. These results can be used for parameter estimation, and simulations are used to both verify the theory and to assess the accuracy of the estimators. The methodology is applied to an investigation of regularity in plumes observed from swimming microorganisms.

In many disciplines, it is of great importance to match objects. Procrustes analysis is a popular method for comparing labeled point configurations based on a least squares criterion. We consider alternative procedures that are highly resistant to outlier points, and we describe an application in electrophoretic gel matching. We consider procedures based on S estimators, least median of squares, and least quartile difference estimators. Practical implementation issues are discussed, including random subset selection and intelligent subset selection (where subsets with small size or near collinear subsets are ignored). The relative performances of the resistant and Procrustes methods are examined in a simulation study.