Quantitative analysis of gel electrophoretograms by image analysis and least squares modeling

Summit: Automated Analysis of Arrayed Single-Cell Gel Electrophoresis

New pipelines are required to automate the quantitation of emerging high-throughput electrophoretic (EP) assessment of DNA damage, or proteoform expression in single cells. EP cytometry consists of thousands of Western blots performed on a microscope slide-sized gel microwell array for single cells. Thus, EP cytometry images pose an analysis challenge that blends requirements for accurate and reproducible analysis encountered for both standard Western blots and protein microarrays. Here, we introduce the Summit algorithm to automate array segmentation, peak background subtraction, and Gaussian fitting for EP cytometry. The data structure storage of parameters allows users to perform quality control on identically processed data, yielding a ~6.5% difference in coefficient of quartile variation (CQV) of protein peak area under the curve (AUC) distributions measured by four users. Further, inspired by investigations of background subtraction methods to reduce technical variation in protein microarray measurements, we aimed to understand the trade-offs between EP cytometry analysis throughput and variation. We found an 11%-50% increase in protein peaks that passed quality control with a subtraction method similar to microarray "average on-boundary" versus an axial subtraction method. The background subtraction method only mildly influences AUC CQV, which varies between 1% and 4.5%. Finally, we determined that the narrow confidence interval for peak location and peak width parameters from Gaussian fitting yield minimal uncertainty in protein sizing. The AUC CQV differed by only ~1%-2% when summed over the peak width bounds versus the 95% peak width confidence interval. We expect Summit to be broadly applicable to other arrayed EP separations, or traditional Western blot analysis.

Quantitative computerized western blotting in detail

The analysis of antibody reactivity against multiple antigens separated according to their molecular weights is facilitated by western blotting. The distinction between immune dominant and recessive antigens is often difficult and carried out by qualitative or empirical means. Quantitative computerized western blotting (QCWB) analyzes reactivity to specific antigens by providing a statistically measurable value for each band allowing differentiation between immunodominant and immunorecessive determinants. QCWB is useful for both single time point analysis and longitudinal studies where multiple time points are evaluated and the relativities against individual bands compared. This technique can be employed to study humoral responses to complex antigenic mixtures such as allergens and infectious agents, or identify serologic markers for early diagnosis of cancer, autoimmune or infectious diseases, or to monitor patient's clinical status.

Spot quantification in two dimensional gel electrophoresis image analysis: comparison of different approaches and presentation of a novel compound fitting algorithm

Background

Various computer-based methods exist for the detection and quantification of protein spots in two dimensional gel electrophoresis images. Area-based methods are commonly used for spot quantification: an area is assigned to each spot and the sum of the pixel intensities in that area, the so-called volume, is used a measure for spot signal. Other methods use the optical density, i.e. the intensity of the most intense pixel of a spot, or calculate the volume from the parameters of a fitted function.

Results

In this study we compare the performance of different spot quantification methods using synthetic and real data. We propose a ready-to-use algorithm for spot detection and quantification that uses fitting of two dimensional Gaussian function curves for the extraction of data from two dimensional gel electrophoresis (2-DE) images. The algorithm implements fitting using logical compounds and is computationally efficient. The applicability of the compound fitting algorithm was evaluated for various simulated data and compared with other quantification approaches. We provide evidence that even if an incorrect bell-shaped function is used, the fitting method is superior to other approaches, especially when spots overlap. Finally, we validated the method with experimental data of urea-based 2-DE of Aβ peptides andre-analyzed published data sets. Our methods showed higher precision and accuracy than other approaches when applied to exposure time series and standard gels.

Conclusion

Compound fitting as a quantification method for 2-DE spots shows several advantages over other approaches and could be combined with various spot detection methods.

The algorithm was scripted in MATLAB (Mathworks) and is available as a supplemental file.

Automatic analysis of 2D polyacrylamide gels in the diagnosis of DNA polymorphisms

Introduction

The analysis of polyacrylamide gels is currently carried out manually or automatically. In the automatic method, there are limitations related to the acceptable degree of distortion of lane and band continuity. The available software cannot deal satisfactorily with this type of situations. Therefore, the paper presents an original image analysis method devoid of the aforementioned drawbacks.

Material

This paper examines polyacrylamide gel images from Li-Cor DNA Sequencer 4300S resulting from the use of the electrophoretic separation of DNA fragments. The acquired images have a resolution dependent on the length of the analysed DNA fragments and typically it is M_G×N_G=3806×1027 pixels. The images are saved in TIFF format with a grayscale resolution of 16 bits/pixel. The presented image analysis method was performed on gel images resulting from the analysis of DNA methylome profiling in plants exposed to drought stress, carried out with the MSAP (Methylation Sensitive Amplification Polymorphism) technique.

Results

The results of DNA polymorphism analysis were obtained in less than one second for the Intel Core™ 2 Quad CPU Q9300@2.5GHz, 8GB RAM. In comparison with other known methods, specificity was 0.95, sensitivity = 0.94 and AUC (Area Under Curve) = 0.98.

Conclusions

It is possible to carry out this method of DNA polymorphism analysis on distorted images of polyacrylamide gels. The method is fully automatic and does not require any operator intervention. Compared with other methods, it produces the best results and the resulting image is easy to interpret. The presented method of measurement is used in the practical analysis of polyacrylamide gels in the Department of Genetics at the University of Silesia in Katowice, Poland.

A computerized methodology for improved virus typing by PCR-RFLP gel electrophoresis

The analysis of digitized images from polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP)gel electrophoresis examinations is a popular method for virus typing, i.e., for identifying the virus type(s) that have infected an investigated biological sample. However, being mostly manual, the conventional virus typing protocol remains laborious, time consuming, and error prone. In order to overcome these shortcomings,we propose a computerized methodology for improving virus typing via PCR-RFLP gel electrophoresis. A novel realistic observation model of the viral DNA motion on the gel matrix is employed to assist in exploiting additional virus-related information in comparison to the conventional approaches. The extracted rich information is fed to a novel typing algorithm, resulting in faster and more accurate decisions. The proposed methodology is evaluated for the case of the human papillomavirus typing on a dataset of 80 real and 1500 simulated samples, producing very satisfactory results.Ind

GelBandFitter--a computer program for analysis of closely spaced electrophoretic and immunoblotted bands

GelBandFitter is a computer program that uses non-linear regression techniques to fit mathematical functions to densitometry profiles of protein gels. This allows for improved quantification of gels with partially overlapping and potentially asymmetric protein bands. The program can also be used to analyze immunoblots with closely spaced bands. GelBandFitter was developed in Matlab and the source code and/or a Windows executable file can be downloaded at no cost to academic users from http://www.gelbandfitter.org.

Quantitative computerized western blotting

Western blotting allows analysis of antibody reactivity against multiple antigens separated according to their molecular weights. The distinction between immune dominant and recessive antigens is often difficult and carried out by qualitative or empirical means. Quantitative computerized western blotting (QCWB) addresses this difficulty by analyzing reactivity to specific antigens and providing a statistically measurable value for each band. This allows differentiation between immunodominant and immunorecessive determinants. QCWB is appropriate for either single time point analysis or longitudinal studies where multiple time points are evaluated and the reactivities against individual bands compared. This technique can be used to study humoral responses to complex antigenic mixtures such as allergens and infectious agents, or to identify serologic markers for early diagnosis of cancer, autoimmune, or infectious diseases, or to monitor patient's clinical status.

Kinetics and diagnostic and prognostic potential of quantitative Western blot analysis and antigen-specific enzyme-linked immunosorbent assay in experimental canine leishmaniasis

Quantitative computerized Western blot analysis of antibody responses during experimental canine Leishmania infantum infection distinguished between immunodominant and nonimmunodominant protein bands. Six infected beagles, positive by both PCR and parasite culture, were monitored over 75 weeks postinfection and during a 12-week allopurinol treatment course. All dogs were symptomatic at the time of treatment. Of 12 antigenic bands examined, the immunodominant bands (12, 14, 24, 29, 48, and 68 kDa) showed significantly increased intensities (P<0.01) and higher frequencies of recognition than the nonimmunodominant bands at all time points. Detection of the former bands at 6 weeks postinfection preceded seroconversion by enzyme-linked immunosorbent assay (ELISA) both on crude Leishmania antigen or the recombinant proteins rK39 and HSP70. Reactivity with the 14-, 48-, and 68-kDa bands signified early infection, whereas increased reactivity with the 14-, 24-, and 29-kDa bands was associated with posttreatment parasite persistence and potential unfavorable prognosis. Total lane intensity (TLI) emerged as a sensitive marker for early infection and increased as early as 4 weeks postinfection. TLI had a significantly higher (P<0.01) relative increase rate than crude Leishmania antigen or HSP70 or rK39 ELISA at all time points. These immunodominant antigens and TLI, as determined by quantitative Western blotting, will be valuable for early detection and treatment evaluation of canine leishmaniasis.

Quantitative microfluidic separation of DNA in self-assembled magnetic matrixes

We present an experimental study of the microfluidic electrophoresis of long DNA in self-assembling matrixes of magnetic bead columns. Results are presented for the rapid separation of lambda-phage, 2lambda-DNA, and bacteriophage T4 DNA, where separation resolutions greater than 2 between lambda and T4 are achieved in times as short as 150 s. The use of a computer-piloted flow control system and injection results in high reproducibility between separations. We compare the experimentally measured mobility and dispersion with an exactly solvable lattice Monte Carlo model. The theory predicts that the mean velocity scales linearly with the field, the band broadening scales with the inverse of the field, and the resolution is independent of the field for intermediate fields-all of which are in accord with the experimental results. Moreover, reasonable quantitative agreement is achieved for band broadening for longer DNA (2lambda and T4) when the average postengagement time is measured experimentally. This work demonstrates the possibility of achieving fast microfluidic separation of large DNA on a routine basis.

An alternative method for electrophoretic gel image analysis in the GelMaster software

A novel methodology of electrophoretic gel image analysis has been proposed that is based on two-dimensional image processing methods instead of previously used one-dimensional Gaussian deconvolution. The crucial problem of the analysis of imperfect gels, that consists in band detection, is solved using the algorithms of band boundary detection and intensity homogeneity indication. The template approach represents the core element of the developed algorithms. The GelMaster software system has been developed in which the novel algorithms are implemented. It involves two-stage interaction with the user: detection of the true bands and deleting the false band detections. The main features of the GelMaster system and the most important algorithms are described.

Chlamydomonas sensory rhodopsins A and B: cellular content and role in photophobic responses

Two retinylidene proteins, CSRA and CSRB, have recently been shown by photoelectrophysiological analysis of RNAi-transformants to mediate phototaxis signaling in Chlamydomonas reinhardtii. Here we report immunoblot detection of CSRA and CSRB apoproteins in C. reinhardtii cells enabling assessment of the cellular content of the receptors. We obtain 9 x 10(4) CSRA and 1.5 x 10(4) CSRB apoprotein molecules per cell in vegetative cells of the wild-type strain 495, a higher value than that for functional receptor cellular content estimated previously from photosensitivity measurements and retinal extraction yields. Exploiting our ability to control the CSRA/CSRB ratio by transformation with receptor gene-directed RNAi, we report analysis of the CSRA and CSRB roles in the photophobic response of the organism by action spectroscopy with automated cell tracking/motion analysis. The results show that CSRA and CSRB each mediate the photophobic swimming response, a second known retinal-dependent photomotility behavior in C. reinhardtii. Due to the different light saturation and spectral properties of the two receptors, CSRA is dominantly responsible for photophobic responses, which appear at high light intensity.

A software system for data analysis in automated DNA sequencing

Software for gel image analysis and base-calling in fluorescence-based sequencing consisting of two primary programs, BaseFinder and GelImager, is described. BaseFinder is a framework for trace processing, analysis, and base-calling. BaseFinder is highly extensible, allowing the addition of trace analysis and processing modules without recompilation. Powerful scripting capabilities combined with modularity and multilane handling allow the user to customize BaseFinder to virtually any type of trace processing. We have developed an extensive set of data processing and analysis modules for use with the program in fluorescence-based sequencing. GelImager is a framework for gel image manipulation. It can be used for gel visualization, lane retracking, and as a front end to the Washington University Getlanes program. The programs were designed using a cross-platform development environment, currently allowing them to run in Windows NT, Windows 95, Openstep/Mach, and Rhapsody. Work is ongoing to deploy the software on additional platforms, including Solaris, Linux, and MacOS. This software has been thoroughly tested and debugged in the analysis of >2 million bp of raw sequence data from human chromosome 19 region q13. Overall sequencing accuracy was measured using a significant subset of these data, consisting of approximately 600 sequences, by comparing the individual shotgun sequences against the final assembled contigs. Also, results are reported from experiments that analyzed the accuracy of the software and two other well-known base-calling programs for sequencing the M13mp18 vector sequence. [The sequence data described in this paper have been submitted to the GenBank data library under accession no. AF025422]

Basic strategies for valid cytometry using image analysis

The present review provides a starting point for setting up an image analysis system for quantitative densitometry and absorbance or fluorescence measurements in cell preparations, tissue sections or gels. Guidelines for instrumental settings that are essential for the valid application of image analysis in cytophotometry and cytofluorometry are described. The general principles of the working mechanism of CCD cameras in combination with general methods to improve the behaviour of the cameras are presented. Optimization of illumination of microscopical and macroscopical objects receives special attention because of its importance for valid cytometry. Sources of errors in quantitative measurements are listed and step-by-step charts for tuning the CCD camera, frame grabber and illumination for the optimal use of the systems are described. Suggestions are given for improvement of image arithmetic in difficult imaging situations, such as low fluorescence signals and high absorbance signals.

Plasmid pGA1 from Corynebacterium glutamicum codes for a gene product that positively influences plasmid copy number

The complete nucleotide sequence (4,826 bp) of the cryptic plasmid pGA1 from Corynebacterium glutamicum was determined. DNA sequence analysis revealed four putative coding regions (open reading frame A [ORFA], ORFA2, ORFB, and ORFC). ORFC was identified as a rep gene coding for an initiator of plasmid replication (Rep) according to the high level of homology of its deduced amino acid sequence with the Rep proteins of plasmids pSR1 (from C. glutamicum) and pNG2 (from Corynebacterium diphtheriae). This function was confirmed by deletion mapping of the minimal replicon of pGA1 (1.7 kb) which contains only ORFC. Deletion derivatives of pGA1 devoid of ORFA exhibited significant decreases in the copy number in C. glutamicum cells and displayed segregational instability. Introduction of ORFA in trans into the cells harboring these deletion plasmids dramatically increased their copy number and segregational stability. The ORFA gene product thus positively influences plasmid copy number. This is the first report on such activity associated with a nonintegrating bacterial plasmid. The related plasmids pGA1, pSR1, and pNG2 lacking significant homology with any other plasmid seem to be representatives of a new group of plasmids replicating in the rolling-circle mode.

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.

Protein profiles of Streptomyces aureofaciens producing tetracyclines: reappraisal of the effect of benzyl thiocyanate

Cell protein profiles of submerged cultures of Streptomyces aureofaciens cultivated in the absence or presence of 12 microM benzyl thiocyanate (BT) were analyzed by one-dimensional SDS polyacrylamide gel electrophoresis. Substantial increase in the intensity of the 13, 35, 37, 60, and 100 kDa protein bands was observed in cultures treated with BT. Similar increase in the 35, 37, and 60 kDa bands was found in a mutant blocked in the last chlortetracycline biosynthesis step. Effect of BT on the solid medium-grown cultures was also observed, with a more intensive substrate mycelium pigmentation and alteration in the spore size and shape as the most characteristic features. Earlier studies of BT effect involving those on the stimulation of chlortetracycline biosynthesis are summarized and a possible signal-transducing mechanism is discussed from the point of view of adaptation of S. aureofaciens to the uncoupling of oxidative phosphorylation.

Two-dimensional electrophoresis of proteins: an updated protocol and implications for a functional analysis of the genome

The two-dimensional electrophoresis (2-DE) technique developed by Klose in 1975 (Humangenetik 1975, 26, 211-234), independently of the technique developed by O'Farrell (J. Biol. Chem. 1975, 250, 4007-4021), has been revised in our laboratory and an updated protocol is presented. This protocol is the result of our experience in using this method since its introduction. Many modifications and suggestions found in the literature were also tested and then integrated into our original method if advantageous. Gel and buffer composition, size of gels, use of stacking gels or not, necessity of isoelectric focusing (IEF) gel incubation, freezing of IEF gels or immediate use, carrier ampholytes versus Immobilines, regulation of electric current, conditions for staining and drying the gels - these and other problems were the subject of our concern. Among the technical details and special equipment which constitute our 2-DE method presented here, a few features are of particular significance: (i) sample loading onto the acid side of the IEF gel with the result that both acidic and basic proteins are well resolved in the same gel; (ii) use of large (46 x 30 cm) gels to achieve high resolution, but without the need of unusually large, flat gel equipment; (iii) preparation of ready-made gel solutions which can be stored frozen, a prerequisite, among others, for high reproducibility. Using the 2-DE method described we demonstrate that protein patterns revealing more than 10 000 polypeptide spots can be obtained from mouse tissues. This is by far the highest resolution so far reported in the literature for 2-DE of complex protein mixtures. The 2-DE patterns were of high quality with regard to spot shape and background. The reproducibility of the protein patterns is demonstrated and shown to be thoroughly satisfactory. An example is given to show how effectively 2-DE of high resolution and reproducibility can be used to study the genetic variability of proteins in an interspecific mouse backcross (Mus musculus x Mus spretus) established by the European Backcross Collaborative Group for mapping the mouse genome. We outline our opinion that the structural analysis of the human genome, currently pursued most intensively on a worldwide scale, should be accompanied by a functional analysis of the genome that starts from the proteins of the organism.

The use of protein synthesis elongation factor EF-Tu as internal calibration standard in two-dimensional electrophoretic studies of differentiation in Streptomyces

Protein synthesis elongation factor EF-Tu is presented as an internal calibration standard for quantitative analysis of two-dimensional (2-D) protein electrophoresis gels. EF-Tu was selected on the basis of concentration measurements in cell-free extracts from Streptomyces aureofaciens, grown under conditions leading to production of tetracyclines, and separated on one-dimensional (1-D) and 2-D electrophoresis gels. The results demonstrated that the amount of EF-Tu synthesized in S. aureofaciens under conditions of slow growth during production of tetracyclines is constant in proportion to all other de novo synthesized proteins regardless of their total number. This makes EF-Tu an ideal internal protein standard for quantitation of protein spots on 2-D electrophoresis gels. For such quantitative analysis we developed a computer-aided image analysis system which provides preparation of a gel image for further analysis including calibration, background subtraction and cleaning for streaking in both directions. The system can locate any resolvable spot in the gel and measure the integrated density of the spot, even in the case of irregular spot shape in crowded and overlapping spot regions.