A program system for interactive measurements on digitized cell images.

Using a high precision image scanner and a PDP-8/F minicomputer, we have developed a program system for interactive measurements on microscopic images. By giving simple keyboard commands, the operator can run the image scanner and manipulate the digitized images. The interface between the operator and the microscope-computer system is a Tektronix 4010 graphic terminal. The system allows objects to be isolated and parameters to be calculated from each object, e.g., parameters characterizing shape of the object, irregularity in light transmission over the object, area, integrated light transmission, etc. Objects are isolated and parameters are calculated under complete operator control using interactive computer graphics technique. Calculated parameters may be stored in dedicated data records, which are stored in files for later statistical analysis. The system also includes a statistical evaluation part. Technically, the system consists of a command scanner, which translates commands into internal representation, a parser, which checks the syntax of the commands, and an interpreter, which executes the commands. The system is designed so that new commands can be added easily.     

The glycolyzer: Automated glycan annotation software for high performance mass spectrometry and its application to ovarian cancer glycan biomarker discovery

Human serum glycomics is a promising method for finding cancer biomarkers but often lacks the tools for streamlined data analysis. The Glycolyzer software incorporates a suite of analytic tools capable of identifying informative glycan peaks out of raw mass spectrometry data. As a demonstration of its utility, the program was used to identify putative biomarkers for epithelial ovarian cancer from a human serum sample set. A randomized, blocked, and blinded experimental design was used on a discovery set consisting of 46 cases and 48 controls. Retrosynthetic glycan libraries were used for data analysis and several significant candidate glycan biomarkers were discovered via hypothesis testing. The significant glycans were attributed to a glycan family based on glycan composition relationships and incorporated into a linear classifier motif test. The motif test was then applied to the discovery set to evaluate the disease state discrimination performance. The test provided strongly predictive results based on receiver operator characteristic curve analysis. The area under the receiver operator characteristic curve was 0.93. Using the Glycolyzer software, we were able to identify a set of glycan biomarkers that highly discriminate between cases and controls, and are ready to be formally validated in subsequent studies.     

Automated analysis of behavior: a computer-controlled system for drug screening and the investigation of learning

Efforts to understand cognition will be greatly facilitated by computerized systems that enable the automated analysis of animal behavior. A number of controversies in the invertebrate learning field have resulted from difficulties inherent in manual experiments. Driven by the necessity to overcome these problems during investigation of neural function in planarian flatworms and frog larvae, we designed and developed a prototype for an inexpensive, flexible system that enables automated control and analysis of behavior and learning. Applicable to a variety of small animals such as flatworms and zebrafish, this system allows automated analysis of innate behavior, as well as of learning and memory in a plethora of conditioning paradigms. We present here the schematics of a basic prototype, which overcomes experimenter effects and operator tedium, enabling a large number of animals to be analyzed with transparent on-line access to primary data. A scaled-up version of this technology represents an efficient methodology to screen pharmacological and genetic libraries for novel neuroactive reagents of basic and biomedical relevance.     

A computer-based videodensitometric system for studying banded human chromosomes illustrated by the analysis of the normal morphology of chromosome 18.

A computer-based, high-resolution, high-speed system is described for the digitization of television images of human G-banded chromosomes, on-line directly from a microscope. The digitized data are processed by a computer to determine the total length and centromere index of individual chromosomes and to obtain an integrated density profile representing the band pattern. The density profile is reduced by the computer to a series of dark and light bands, each with a defined position, width, and density. The results are displayed by the computer on a television monitor in a form that simulates routine laboratory preparations and that allows cytogenetic inspection. The system is easy to operate, and it includes a wide range of options for operator interaction to ensure reliability. The results of the computer analysis of the no. 18 chromosomes from five cells from each of 10 different persons demonstrated the usefulness of the system for detecting, measuring, and analyzing individual bands in human populations and its application to establishing band patterns of chromosomes in different states of contraction.     

A new semi-automated measurement technique based on X-ray pictures for ankle morphometry

Knowledge of human joint morphology is important in orthopaedic surgery and in prosthesis design. The literature on quantitative morphological analysis of the ankle joint is particularly scarce. A semi-automated radiographic measurement method was developed to collect morphological measures of the ankle joint. The method was based on standard lateral and antero-posterior X-ray pictures of the ankle joint. These were then scanned and analysed by means of specialized software designed for the purpose, which requires minimal operator contribution. Accuracy of the method was experimentally assessed by in vitro direct measurements. Intra- and inter-operator variability was also tested. Accuracy was assessed to be within 1mm for most measurements. Repeatability was not affected by operator skill and was within 2mm. The newly proposed method was applied successfully on 15 male adult subjects and relevant results are reported. The method allows ankle morphology to be analysed in a large number of subjects providing reliable data for anthropometric statistics.     

Quantitative analysis of electrophoresis data: novel curve fitting methodology and its application to the determination of a protein-DNA binding constant.

A computer program, GelExplorer, which uses a new methodology for obtaining quantitative information about electrophoresis has been developed. It provides a straightforward, easy-to-use graphical interface, and includes a number of features which offer significant advantages over existing methods for quantitative gel analysis. The method uses curve fitting with a nonlinear least-squares optimization to deconvolute overlapping bands. Unlike most curve fitting approaches, the data is treated in two dimensions, fitting all the data across the entire width of the lane. This allows for accurate determination of the intensities of individual, overlapping bands, and in particular allows imperfectly shaped bands to be accurately modeled. Experiments described in this paper demonstrate empirically that the Lorentzian lineshape reproduces the contours of an individual gel band and provides a better model than the Gaussian function for curve fitting of electrophoresis bands. Results from several fitting applications are presented and a discussion of the sources and magnitudes of uncertainties in the results is included. Finally, the method is applied to the quantitative analysis of a hydroxyl radical footprint titration experiment to obtain the free energy of binding of the lambda repressor protein to the OR1 operator DNA sequence.     

Automated analysis of behavior: A computer‐controlled system for drug screening and the investigation of learning

Efforts to understand cognition will be greatly facilitated by computerized systems that enable the automated analysis of animal behavior. A number of controversies in the invertebrate learning field have resulted from difficulties inherent in manual experiments. Driven by the necessity to overcome these problems during investigation of neural function in planarian flatworms and frog larvae, we designed and developed a prototype for an inexpensive, flexible system that enables automated control and analysis of behavior and learning. Applicable to a variety of small animals such as flatworms and zebrafish, this system allows automated analysis of innate behavior, as well as of learning and memory in a plethora of conditioning paradigms. We present here the schematics of a basic prototype, which overcomes experimenter effects and operator tedium, enabling a large number of animals to be analyzed with transparent on‐line access to primary data. A scaled‐up version of this technology represents an efficient methodology to screen pharmacological and genetic libraries for novel neuroactive reagents of basic and biomedical relevance. © 2006 Wiley Periodicals, Inc. J Neurobiol, 2006     

Robust functional principal component analysis via a functional pairwise spatial sign operator

Functional principal component analysis (FPCA) has been widely used to capture major modes of variation and reduce dimensions in functional data analysis. However, standard FPCA based on the sample covariance estimator does not work well if the data exhibits heavy-tailedness or outliers. To address this challenge, a new robust FPCA approach based on a functional pairwise spatial sign (PASS) operator, termed PASS FPCA, is introduced. We propose robust estimation procedures for eigenfunctions and eigenvalues. Theoretical properties of the PASS operator are established, showing that it adopts the same eigenfunctions as the standard covariance operator and also allows recovering ratios between eigenvalues. We also extend the proposed procedure to handle functional data measured with noise. Compared to existing robust FPCA approaches, the proposed PASS FPCA requires weaker distributional assumptions to conserve the eigenspace of the covariance function. Specifically, existing work are often built upon a class of functional elliptical distributions, which requires inherently symmetry. In contrast, we introduce a class of distributions called the weakly functional coordinate symmetry (weakly FCS), which allows for severe asymmetry and is much more flexible than the functional elliptical distribution family. The robustness of the PASS FPCA is demonstrated via extensive simulation studies, especially its advantages in scenarios with nonelliptical distributions. The proposed method was motivated by and applied to analysis of accelerometry data from the Objective Physical Activity and Cardiovascular Health Study, a large-scale epidemiological study to investigate the relationship between objectively measured physical activity and cardiovascular health among older women.     

Determining prevalence of bluetongue and epizootic hemorrhagic disease viruses in mule deer in Arizona (USA) using whole blood dried on paper strips compared to serum analyses

We investigated the feasibility of using whole blood dried on paper strips as a means to collect antibody prevalence data for the epizootic hemorrhagic disease viruses (EHDV) and bluetongue viruses (BTV) from hunter-harvested male mule deer (Odocoileus hemionus) in October 2002 from Arizona, USA. We compared antibody prevalence estimates in mule deer from paired paper strip and serum samples. Prevalence data obtained from elution of dried blood on paper strips proved to be consistent with results from serum in 94% of the samples tested. The paper strip method allows easy collection of blood from dead animals, with a smaller amount of blood being needed for analyses. Also, samples do not need to be refrigerated before analyses. We also used serum samples to determine hemorrhagic disease (HD) serotype exposure status of mule deer harvested from 4 distinct areas in Arizona. Antibodies to BTV and EHDV were identified in 3 of the 4 areas, with positive results to EHDV-1, EHDV-2, BTV-10, and BTV-11 being most common. Many animals did not have antibodies against the BTV serotypes. Exposure varied geographically and potentially with elevation. Hemorrhagic disease viruses commonly infect Arizona mule deer, except on the Kaibab Plateau in northern Arizona.     

THREE-TAXON STATEMENTS: A MORE PRECISE USE OF PARSIMONY?

Abstract— Binary characters can be represented in data matrices by the three-taxon statements they imply. Transforming characters into three-taxon statements may increase the sensitivity of parsimony to differences in the fit of data to alternative cladograms. Extrapolation of the technique to multistate features allows semi-additive characters to be coded accurately. In many cases, analysis of the transformed data produces fewer equally parsimonious solutions than does analysis of the raw data. In other cases, additional equally parsimonious solutions, or even different solutions, may be produced; in those cases, the results appear to accommodate information from a larger number of characters than do the results from raw data.     

An alpha-numeric code for representing N-linked glycan structures in secreted glycoproteins

Advances in high-throughput techniques have led to the creation of increasing amounts of glycome data. The storage and analysis of this data would benefit greatly from a compact notation for describing glycan structures that can be easily stored and interpreted by computers. Towards this end, we propose a fixed-length alpha-numeric code for representing N-linked glycan structures commonly found in secreted glycoproteins from mammalian cell cultures. This code, GlycoDigit, employs a pre-assigned alpha-numeric index to represent the monosaccharides attached in different branches to the core glycan structure. The present branch-centric representation allows us to visualize the structure while the numerical nature of the code makes it machine readable. In addition, a difference operator can be defined to quantitatively differentiate between glycan structures for further analysis. The usefulness and applicability of GlycoDigit were demonstrated by constructing and visualizing an N-linked glycosylation network.     

Bioinformatics support for high-throughput proteomics

In the "post-genome" era, mass spectrometry (MS) has become an important method for the analysis of proteome data. The rapid advancement of this technique in combination with other methods used in proteomics results in an increasing number of high-throughput projects. This leads to an increasing amount of data that needs to be archived and analyzed.To cope with the need for automated data conversion, storage, and analysis in the field of proteomics, the open source system ProDB was developed. The system handles data conversion from different mass spectrometer software, automates data analysis, and allows the annotation of MS spectra (e.g. assign gene names, store data on protein modifications). The system is based on an extensible relational database to store the mass spectra together with the experimental setup. It also provides a graphical user interface (GUI) for managing the experimental steps which led to the MS data. Furthermore, it allows the integration of genome and proteome data. Data from an ongoing experiment was used to compare manual and automated analysis. First tests showed that the automation resulted in a significant saving of time. Furthermore, the quality and interpretability of the results was improved in all cases.     

Low cost,microprocessor-based system for flow evaluation in a hydraulic model using the dilution method

An Automatic Dye Dilution Processor (ADP) based on a microprocessor is applied to a hydraulic model to process on-line dilution curves, using saline or thermodilution methods. A low cost, F8 microcomputer with minimum parts configuration was used in an Automatic Dilution Processor (ADP) system. Software allows an operator to select methods, amount of indicator to be injected, cell calibration constant, and temperature step. Care has been taken in the curve data processing, to avoid errors when recirculation occurs. The system is used, together with an open hydraulic model, for flow measurement demonstration, and the estimation of cell calibration constant.     

DNA signal splitting improves detection and analysis of tetraploid populations.

Detection of DNA tetraploid populations requires a high index of suspicion at the time of data acquisition and frequently requires subsequent off-line analysis for confirmation, including evaluation of the hypertetraploid region. To analyze these specimens, the flow cytometer operator must run all specimens with the G0G1 peak in low channels or rerun specimens in which tetraploidy is suspected with a lower photomultiplier tube (PMT) voltage or lower amplifier gain setting. Re-analysis may not be possible in specimens with few cells. A simple modification to the cytometer allows PMT signal splitting with simultaneous processing of the signal by two different amplifiers. This allows simultaneous acquisition of histograms optimized for both the hypotetraploid and hypertetraploid regions.     

The interaction of the recognition helix of lac repressor with lac operator.

We have constructed a system which allows systematic testing of repressor--operator interactions. The system consists of two plasmids. One of them carries a lac operon in which lac operator has been replaced by a unique restriction site into which synthetic operators can be cloned. The other plasmid carries the gene coding for the repressor, in our case a semisynthetic lacI gene of which parts can be exchanged in a cassette-like manner. A galE host allows us to select for mutants which express repressors with altered specificities. Here we report the change of specificity in the lac system by changing residues 1 and 2 of the recognition helix of lac repressor. The specificity changes are brought about cooperatively by the change of both residues. Exchanges of just one residue broaden the specificity. Our results hint that the recognition helix of lac repressor may possibly have the opposite orientation to those in Lambda cro protein or 434 CI repressor.     

A flow injection flow cytometry system for on-line monitoring of bioreactors

For direct and on-line study of the physiological states of cell cultures, a robust flow injection system has been designed and interfaced with flow cytometry (FI-FCM). The core of the flow injection system includes a microchamber designed for sample processing. The design of this microchamber allows not only an accurate on-line dilution but also on-line cell fixation, staining, and washing. The flow injection part of the system was tested by monitoring the optical density of a growing E.coli culture on-line using a spectrophotometer. The entire growth curve, from lag phase to stationary phase, was obtained with frequent sampling. The performance of the entire FI-FCM system is demonstrated in three applications. The first is the monitoring of green fluorescent protein fluorophore formation kinetics in E.coli by visualizing the fluorescence evolution after protein synthesis is inhibited. The data revealed a subpopulation of cells that do not become fluorescent. In addition, the data show that single-cell fluorescence is distributed over a wide range and that the fluorescent population contains cells that are capable of reaching significantly higher expression levels than that indicated by the population average. The second application is the detailed flow cytometric evaluation of the batch growth dynamics of E.coli expressing Gfp. The collected single-cell data visualize the batch growth phases and it is shown that a state of balanced growth is never reached by the culture. The third application is the determination of distribution of DNA content of a S. cerevisiae population by automatically staining cells using a DNA-specific stain. Reproducibility of the on-line staining reaction shows that the system is not restricted to measuring the native properties of cells; rather, a wider range of cellular components could be monitored after appropriate sample processing. The system is thus particularly useful because it operates automatically without direct operator supervision for extended time periods.     

Detection of serum proteins by native polyacrylamide gel electrophoresis using Blue Sepharose CL-6B-containing stacking gels

Analysis of serum proteins by native polyacrylamide gel electrophoresis is difficult because albumin is abundant in serum and interferes with the resolution of other proteins, especially alpha-antitrypsin which has mobility that is very similar to that of albumin. We present here a method in which serum proteins are separated by polyacrylamide gel electrophoresis using stacking gels containing Blue Sepharose CL-6B, which has a high affinity for albumin, lipoproteins, kinases, and pyridine-nucleotide-dependent oxidoreductases. During electrophoresis, proteins that bind to Blue Sepharose CL-6B stay in the stacking gel and do not migrate into the separating gel. As a consequence, certain proteins, including alpha(1)-antitrypsin, can be detected as clear bands. This method overcomes the requirement for fractionation of serum samples prior to electrophoresis to remove albumin and allows the simultaneous analysis of many samples.     

Rapid selection of tetracycline-controlled inducible cell lines using a green fluorescent-transactivator fusion protein

We describe a modification of the tetracycline-controlled expression system that facilitates the rapid identification of tetracycline-sensitive clones. The TetR/VP16 transactivator protein was tagged with the green fluorescent protein (GFP) at its N-terminus. This results in a functional transactivator which allows cells expressing high levels of the modified transactivator to be selected by fluorescent-activated cell sorting. After expansion, single cell clones that maintain a high level of GFP fluorescence can be tested for their ability to transactivate a luciferase gene under control of the Tet operator, leading to the rapid identification of clones with strong inducibility.     

Linear trend in multi-species time series

Investigation of permanent plots is the traditional approach to detect changes in species performance and floristic composition. When the time reserved for investigations is limited and statistically independent replicate samples for normal time series analysis do not exist, ordination of multi-species series is often applied. The approach is further developed here with time series data from wetland communities over six consecutive years. Random fluctuation and linear trend are the two mechanisms which can explain the observed changes. Trend analysis of species scores allows to smooth the data and hence the resulting ordination pattern. The expected scores are a conservative measure for trend, taking into account all the recorded time states of the system.