Using standard positions and image fusion to create proteome maps from collections of two-dimensional gel electrophoresis images

Databases for two-dimensional protein gels pose new challenges in extracting meaningful information from large numbers of experiments. In order to create expression profiles, positions of corresponding protein spots across all gel images have to be established. In larger gel sets errors may accumulate rapidly during this spot matching process, effectively limiting the number of samples available for data mining. Here we present a novel approach for organizing spot data based on the concept of a standard position for a protein species. Standard positions are meaningful average positions that are determined using all occurrences of a protein species. They can be extended to spots that are not annotated via interpolation. The standard position of a spot can serve as a unifying index across all gels in a database, thus allowing creation and analysis of expression profiles that span the whole collection. The standard position gives a much more accurate estimation of a spot's position on a gel than can be obtained using theoretical isoelectric point and molecular weight. Positional indexing is a complement to a priori identifications (e.g. by mass spectrometry or Edman degradation). Moreover it can be used in advance to select spots that are worth identifying because they show relevant expression profiles. Furthermore, we show how to combine all spots that occur on any of the gels into one synthetic but nevertheless realistic-looking image. This composite image is produced such that all spots have their standard positions. It can serve as a proteome reference map for an organism. As an application, we have computed a reference map from 23 gel images of Bacillus subtilis, using an enhanced prerelease version of the gel analysis software Delta2D (DECODON, Greifswald, Germany).     

Primary skin fibroblasts as human model system for proteome analysis

Elucidation of cellular processes and their changes at the level of protein expression and post-translational modification patterns may allow identification of novel proteins and thereby mechanisms involved in the pathogenesis of multigenic diseases. The aim of this study was to test cultured, nontransformed primary fibroblasts derived from human skin biopsies as a suitable model system for proteome analysis. Therefore soluble protein fractions were separated on several overlapping ultrazoom gels covering the pH range from 3.5-9. Correlation analysis of gel-pairs revealed a highly reproducible protein expression pattern within (intra-assay) and between (inter-assay) independent experiments of a single fibroblast cell line (intra-cell line comparison). Spot intensity variations were less than a factor of two for more than 80% of identical spots. In addition, inter-cell line comparison exhibits no significant variations in spot intensities. To achieve further improvements in reproducibility we generated master gels for each pH range by combining averaged spot information derived from two different cell lines each analysed by two independent experiments using the raw master gel algorithm of the Z3 image analysis software. The resulting reference images of primary human fibroblasts provided a basis for investigating regulation by extracellular stimuli and drugs as well as their alterations in patients with different diseases.     

Hierarchical analysis of large-scale two-dimensional gel electrophoresis experiments

Large-scale two-dimensional gel experiments have the potential to identify proteins that play an important role in elucidating cell mechanisms and in various stages of drug discovery. Such experiments, typically including hundreds or even thousands of related gels, are notoriously difficult to perform, and analysis of the gel images has until recently been virtually impossible. In this paper we describe a scalable computational model that permits the organization and analysis of a large gel collection. The model is implemented in Compugen's Z4000 system. Gels are organized in a hierarchical, multidimensional data structure that allow the user to view a large-scale experiment as a tree of numerous simpler experiments, and carry out the analysis one step at a time. Analyzed sets of gels form processing units that can be combined into higher level units in an iterative framework. The different conditions at the core of the experiment design, termed the dimensions of the experiment, are transformed from a multidimensional structure to a single hierarchy. The higher level comparison is performed with the aid of a synthetic "adaptor" gel image, called a Raw Master Gel (RMG). The RMG allows the inclusion of data from an entire set of gels to be presented as a gel image, thereby enabling the iterative process. Our model includes a flexible experimental design approach that allows the researcher to choose the condition to be analyzed a posteriori. It also enables data reuse, the performing of several different analysis designs on the same experimental data. The stability and reproducibility of a protein can be analyzed by tracking it up or down the hierarchical dimensions of the experiment.     

Identification and purification of overlapping cosmid clones of the region Xq24-qter of human X chromosome using HPLC.

The assembly of a large physical map of genomes requires simultaneous analysis of many cosmid clones for overlapping regions. The search for overlapping regions may be achieved by various means. High-performance liquid chromatography (HPLC) provides an alternative to gel electrophoresis since microgram amounts of each DNA fragment may be collected into individual test tubes for further analysis. HPLC has been used to identify overlapping cosmid clones from a pool of cosmid DNA containing the terminal portion of the long arm of the human X chromosome (Xq24-qter). Among 400 cosmids analyzed, 3 were shown to overlap.     

Clones from an 840-kb fragment containing the 5' region of the DMD locus enriched by pulsed field gel electrophoresis

Detailed analysis of a large region of genomic DNA is facilitated by generating overlapping clones covering the entire region. These clones are usually obtained by bidirectional "walking" using either bacteriophage lambda or cosmid cloning vectors. This is a slow procedure when starting from a single start site. Multiple start sites are an advantage, and here we describe a method of generating clones from an extensive region of the Duchenne muscular dystrophy locus by preparative pulsed field gel electrophoresis using the chromosome of interest isolated in a cell hybrid. We have generated 12 clones mapping to an 840-kb SfiI fragment of DNA from the Xp2.1 region of the X chromosome, where the DMD gene has been localized. Further localization of these clones to the four subregions of the 840-kb fragment indicates that the clones are distributed throughout the fragment. The feasibility of using this approach to generate probes close to other loci is discussed.     

Challenges related to analysis of protein spot volumes from two-dimensional gel electrophoresis as revealed by replicate gels

Assumptions that need to be considered prior to statistical analysis of protein spot volumes from two-dimensional gel electrophoresis (2-DE) data are studied using replicate gels of the same sample. The most important observation is that the data tables of protein spot volumes from 2-DE images contain a large number of missing values, which are not consistent with the presence or absence of the proteins. This implies both loss of information and problems for the subsequent statistical analysis. Challenges with 2-DE protein spot volumes are viewed in light of multiple gel comparisons and multivariate data analysis.     

Combinatorial use of mRNA and two-dimensional electrophoresis expression data to choose relevant features for mass spectrometric identification

It is only recently that quantitative studies of differential proteome analysis (DPA) have become possible. In this paper the issues involved in quantitative DPA are discussed and novel tools to select features for identification by mass spectrometry (MS) are described. The problem of comparing two sets of gels on a global level is explored as well as how to find specific protein features that differentiate two sets of two-dimensional electrophoresis gels. The concept of a 'virtual' gel, derived from gene expression data, is introduced. The virtual gel enables the co-analysis of data from gene and protein expression. We discuss the value of such an approach, and consider what new information can be gained by using gene and protein expression together. These tools are illustrated by analysis of data from tandem gene and protein expression experiments. Features that are highlighted by the above methods are putative candidates for MS identification. Tools are described that integrate the process of feature selection, cutting, and MS analysis.     

Pixel-based analysis of multiple images for the identification of changes: a novel approach applied to unravel proteome patterns [corrected] of 2-D electrophoresis gel images

A novel approach for revealing patterns of proteome variation among series of 2-DE gel images is presented. The approach utilises image alignment to ensure that each pixel represents the same information across all gels. Gel images are normalised, and background corrected, followed by unfolding of the images to 1-D pixel vectors and analysing pixel vectors by multivariate data modelling. Information resulting from the data analysis is refolded back to the image domain for visualisation and interpretation. The method is rapid and suitable for automatic routines applied after the gel alignment. The approach is compared with spot volume analysis to illustrate how this approach can solve persistent problems like mismatch of protein spots, erroneous missing values and failure to detect variation in overlapping proteins. The method may also detect variation in the border area of saturated proteins. The approach is given the name pixel-based analysis of multiple images for the identification of changes (PMC). The method can be used for multiple images in general. Effects of pretreatment of the images are discussed.     

An optical comparator for measuring two-dimensional polyacrylamide gel electrophoresis records using an on-line microcomputer

A comparator which makes it possible to compare two wet gels or photographic negatives or autoradiograms through a flickering light system has been built. The system consists of two special-purpose projectors which combine the images on a digitizing platform. When the lights are switched On and off out of phase, the positions of the common components remain unchanged, whereas those that are spatially displaced appear to jump from side to side and those present in one image but not the other switch on and off. This produces a flickering image in which differences are readily seen. Commercial camera lenses were used to construct the projectors and the overall specifications for the system are given. The coordinates of both the displaced components, as well as the selected standards from the two images, are digitized and entered automatically into an on-line microcomputer. By using an iterative procedure for collecting records from several superimposable records of the gel, it is possible to compensate for the lack of total reproducibility over the whole gels. These coordinates are then normalized and superimposed on a master map through a television display using a curser to adjust the coordinates. The whole procedure can be repeated for many gels using a common reference gel in the comparator, and the result is a set of normalized coordinates which can be plotted on a single map to provide a final record of the experiments.     

Rapid screening of a YAC library by pulsed-field gel Southern blot analysis of pooled YAC clones

A new method for screening of YAC libraries is described. Individual YACs were pooled into groups of 384 clones and prepared as samples suitable for pulsed-field gel electrophoresis. A five hit human YAC library (Brownstein et al., 1989) containing approximately 60,000 clones was condensed into 150 such pools and chromosomal DNAs in each sample were separated on three pulsed field gels containing 50 samples each. Southern blots prepared from these gels were hybridized with probes of interest to identify pools containing homologous YACs. Further purification was performed using standard colony hybridization procedures. Twenty-one probes used thus far have identified 47 positive pools and corresponding YACs have been purified from 28 of these. Some significant advantages of this method include avoidance of DNA sequence analysis and primer generation prior to YAC screening and the ability to handle the entire library on three filters. The screening approach described here permits rapid isolation of YACs corresponding to unsequenced loci and will accelerate establishment of YAC contigs for large chromosomal segments.     

High-resolution restriction map for a 240-kilobase region spanning 91 to 96 minutes on the Salmonella typhimurium LT2 chromosome.

A hierarchical approach allows the completion of contiguous sets of overlapping clones for small regions of a genome, one at a time rather than tackling the whole genome at once. On the basis of the BlnI restriction map for Salmonella typhimurium LT2, we dissected the chromosome into 21 different fragments by using a Tn5 transposon carrying a BlnI site. Dissected chromosomal fragments were purified by pulsed-field gel electrophoresis and used as probes for sorting a lambda DASHII genomic library of 2,304 primary clones. A total of 129 clones identified as spanning the region from 91 min to 98 min were partly ordered on the basis of the intensity of hybridization with mitomycin-induced Mud-P22 phage DNAs from insertions with pac sites in opposite orientations at 93 min used as probes. Decreased signal intensity with the Mud-P22 probes corresponded to the increased distance of the clone from the site of Mud-P22 insertion and allowed the clones to be placed in two groups from 91 min to 93 min and from 93 min to 98 min and into four intensity categories within the two groups. A member of each category was used to generate a riboprobe from the T3 promoter flanking the insert. This probe identified overlapping clones among the 129 clones. This subchromosomal library was then screened again with riboprobes from nonoverlapping clones. After four cycles of this strategy, a minimal contiguous sequence of 19 partly overlapping clones was selected for restriction mapping. A detailed map of 378 sites for eight restriction enzymes is presented for a region of about 240 kb. Working clockwise, the following genes were placed on this physical map on the basis of their restriction maps: malFEK, lamB, malM, lexA, qor, dnaB, alr, uvrA, proP, pmrB, pmrA, melA, melB, phoN, amiB, mutL, and miaA.     

A bacterial artificial chromosome based physical map of the Ustilago maydis genome.

Ustilago maydis, a basidiomycete, is a model organism among phytopathogenic fungi. A physical map of U. maydis strain 521 was developed from bacterial artificial chromosome (BAC) clones. BAC fingerprints used polyacrylamide gel electrophoresis to separate restriction fragments. Fragments were labeled at the HindIII site and co-digested with HaeIII to reduce fragments to 50-750 bp. Contiguous overlapping sets of clones (contigs) were assembled at nine stringencies (from P < or = 1 x 10(-6) to 1 x 10(-24)). Each assembly nucleated contigs with different percentages of bands overlapping between clones (from 20% to 97%). The number of clones per contig decreased linearly from 41 to 12 from P < or = 1 x 10(-7) to 1 x 10 (-12). The number of separate contigs increased from 56 to 150 over the same range. A hybridization-based physical map of the same BAC clones was compared with the fingerprint contigs built at P < or = 1 x 10(-7). The two methods provided consistent physical maps that were largely validated by genome sequence. The combined hybridization and fingerprint physical map provided a minimum tile path composed of 258 BAC clones (18-20 Mbp) distributed among 28 merged contigs. The genome of U. maydis was estimated to be 20.5 Mbp by pulsed-field gel electrophoresis and 24 Mbp by BAC fingerprints. There were 23 separate chromosomes inferred by both pulsed-field gel electrophoresis and fingerprint contigs. Only 11 of the tile path BAC clones contained recognizable centromere, telomere, and subtelomere repeats (high-copy DNA), suggesting that repeats caused some false merges. There were 247 tile path BAC clones that encompassed about 17.5 Mbp of low-copy DNA sequence. BAC clones are available for repeat and unique gene cluster analysis including tDNA-mediated transformation. Program FingerPrint Contigs maps aligned with each chromosome can be viewed at http://www.siu.edu/~meksem/ustilago_maydis/.     

Elsie 4: quantitative computer analysis of sets of two-dimensional gel electrophoretograms

We have developed and refined a system for quantitative computer analysis of two-dimensional polyacrylamide gel electrophoretograms. The system, named Elsie 4, is based on one described by Vo et al. (Anal. Biochem. 112, 258 (1981]. It is highly automated. Elsie 4 can find, and measure the intensity of, almost any spot resolvable on two-dimensional gels, including spots visible only as shoulders off larger spots and spots so close together that there is no "valley" between them. It can automatically match the spot patterns of different gels, potentially without the need for a user to provide landmark matches. The matches between paired gels let us follow the synthesis of any spot through a set of gels. Information about a group of matched spots can be obtained by referring to any spot in the group. There is generally no need for a standard or reference gel. Data for two experiments can be combined and compared by matching any gel in one experiment with any gel in the other. There are ways to automatically find possible mismatches in sets of gels. Scans and the results of the analysis can be shown on an image displayer. The programs use function libraries; this helps ensure consistency and increase portability. The programs and functions can be linked together in many ways; this lets users build custom programs for analysis of specific experiments.     

The generation of ordered sets of cosmid DNA clones from human chromosome region 11p.

We describe progress in a continuing project aimed at the generation of an overlapping cosmid DNA clone map of the short arm of human chromosome 11. The automated procedures used to prepare DNA samples and the computerized data collection and recording systems are described. We also demonstrate the use of the clones as reagents for the rapid isolation of genomic DNAs containing smaller probed regions. We have isolated approximately 4700 human cosmid DNA clones from mouse/human hybrid cell lines that contain predominantly human chromosomal region 11p. Of the DNA in the cell lines, 60% is derived from this chromosomal region, and the remaining 40% is derived from regions of chromosomes 3, 19, and 20. A total of 4159 clones have been fingerprinted to identify potential overlaps, and we have developed 535 sets ("contigs"). Using random modeling, it is estimated that 65% of 11p must be contained in the analyzed cosmids. The database of clones has been used to identify single or overlapping clones from noncosmid DNA probes. Examples are presented. It is proposed that cosmid reference filters be distributed to requesting laboratories.     

Examination of 2-DE in the Human Proteome Organisation Brain Proteome Project pilot studies with the new RAIN gel matching technique

The Human Proteome Organisation (HUPO) Brain Proteome Project (BPP) pilot studies have generated over 200 2-D gels from eight participating laboratories. This data includes 67 single-channel and 60 DIGE gels comparing 30 whole frozen C57/BL6 female mouse brains, ten each at embryonic day 16, postnatal day 7 (juvenile) and postnatal day 54-56 (adult); and ten single-channel and three DIGE gels comparing human epilepsy surgery of the temporal front lobe with a corresponding post-mortem specimen. The samples were generated centrally and distributed to the participating laboratories, but otherwise no restrictions were placed on sample preparation, running and staining protocols, nor on the 2-D gel analysis packages used. Spots were characterised by MS and the annotated gel images published on a ProteinScape web server. In order to examine the resultant differential expression and protein identifications, we have reprocessed a large subset of the gels using the newly developed RAIN (Robust Automated Image Normalisation) 2-D gel matching algorithm. Traditional approaches use symbolic representation of spots at the very early stages of the analysis, which introduces persistent errors due to inaccuracies in spot modelling and matching. With RAIN, image intensity distributions, rather than selected features, are used, where smooth geometric deformation and expression bias are modelled using multi-resolution image registration and bias-field correction. The method includes a new approach of volume-invariant warping which ensures the volume of protein expression under transformation is preserved. An image-based statistical expression analysis phase is then proposed, where small insignificant expression changes over one gel pair can be revealed when reinforced by the same consistent changes in others. Results of the proposed method as applied to the HUPO BPP data show significant intra-laboratory improvements in matching accuracy over a previous state-of-the-art technique, Multi-resolution Image Registration (MIR), and the commercial Progenesis PG240 package.     

Towards a physical map of the Drosophila melanogaster genome: mapping of cosmid clones within defined genomic divisions.

A physical map of the D. melanogaster genome is being constructed, in the form of overlapping cosmid clones that are assigned to specific polytene chromosome sites. A master library of ca. 20,000 cosmids is screened with probes that correspond to numbered chromosomal divisions (ca. 1% of the genome); these probes are prepared by microdissection and PCR-amplification of individual chromosomes. The 120 to 250 cosmids selected by each probe are fingerprinted by Hinfl digestion and gel electrophoresis, and overlaps are detected by computer analysis of the fingerprints, permitting us to assemble sets of contiguous clones (contigs). Selected cosmids, both from contigs and unattached, are then localized by in situ hybridization to polytene chromosomes. Crosshybridization analysis using end probes links some contigs, and hybridization to previously cloned genes relates the physical to the genetic map. This approach has been used to construct a physical map of the 3.8 megabase DNA in the three distal divisions of the x chromosome. The map is represented by 181 canonical cosmids, of which 108 clones in contigs and 32 unattached clones have been mapped individually by in situ hybridization to chromosomes. Our current database of in situ hybridization results also includes the beginning of a physical map for the rest of the genome: 162 cosmids have been assigned by in situ hybridization to 129 chromosomal subdivisions elsewhere in the genome, representing 5 to 6 megabases of additional mapped DNA.     

Toward a high resolution 2-DE profile of the normal human liver proteome using ultra-zoom gels

The human liver is the largest organ in the body and has many important physiological functions. A global analysis of human liver proteins is essential for a better understanding of the molecular basis of the normal functions of the liver and of its diseases. As part of the Human Liver Proteome Project (HLPP), the goal of the present study was to visualize and detect as many proteins as possible in normal human livers using two-dimensional gel electrophoresis (2-DE). We have constructed a reference map of the proteins of human normal liver that can be used for the comprehensive analysis of the human liver proteome and other related research. To improve the resolution and enhance the detection of low abundance proteins, we developed and optimized narrow pH range ultra-zoom 2-DE gels. High resolution patterns of human liver in pH gradients 4.5–5.5, 5–6, 5.5–6.7, 6–9 and 6–11 are presented. To improve the poor resolution in the alkaline pH range of 2-DE gels, we optimized the isoelectric focusing protocol by including sample application using cup loading at the anode and incorporating 1.2% hydroxyethyl disulfide, 15% 2-propanol and 5% glycerol in the rehydration buffer. Using the optimized protocol, we obtained reproducibly better resolution in both analytical and preparative 2-DE gels. Compared with the 2386 and 1878 protein spots resolved in the wide range 3–10 and 4–7 pH gradients respectively, we obtained 5481 protein spots from the multiple (overlapping) narrow pH range ultra-zoom gels in the range of pH 4.5–9. The visualized reference map of normal human liver proteins presented in this paper will be valuable for comparative proteomic research of the liver proteome.     

The QUEST system for quantitative analysis of two-dimensional gels

The strategies and methods used by the QUEST system for two-dimensional gel analysis are described, and the performance of the system is evaluated. Radiolabeled proteins, resolved on two-dimensional gels and detected using calibrated exposures to film, are quantified in units of disintegrations per minute or as a fraction of the total protein radioactivity applied to the gel. Spot quantitation and resolution of overlapping spots is performed by two-dimensional gaussian fitting. Pattern matching is carried out for groups of gels called matchsets, and within each matchset every gel is matched to every other gel. During the matching process, spots are automatically added to each pattern at positions where unmatched spots were detected in other patterns. This results in enhanced accuracy for both spot detection and for matching. The spot fitting procedure is repeated after matching. Tests show that up to 97% of spots in each pattern can be matched and that fewer than 1% of the spots are matched inconsistently. Approximately 2000 proteins are detected from typical gels. Of these 1600 are high quality spots. Tests to measure the coefficient of variation of spot quantitation versus spot quality show that the average coefficient of variation for high quality spots is 21%. The intensities of the detected proteins range from 4 to 20,000 ppm of total protein synthesis. The QUEST analysis system has been used to build a quantitative database for the proteins of normal and transformed REF52 cells, as presented in the accompanying reports (Garrels, J., and Franza, B. R., Jr. (1989) J. Biol. Chem. 264, 5283-5298, 5299-5312).     

Rapid pulsed field separation of DNA molecules up to 250 kb.

Pulsed field gel electrophoresis (PFGE) is capable of resolving a wide size range of DNA molecules which would all co-migrate in conventional agarose gels. We describe pulsed field gel conditions which permit DNA fragments of up to 250 kilobases (kb) to be separated in only 3.5 h. The separations, which employ commercially available gel boxes, are achieved using conditions which deviate significantly from traditional pulsed field conditions. PFGE separations have been thought to require reorientation angles greater than 90 degrees to be effective. However, reorientation angles of 90 degrees and even less will resolve DNA fragments a few hundred kb and smaller approximately 5 x faster than with standard pulsed field conditions. The mobility of DNA fragments separated with 90 degrees reorientation angles is switch time-dependent, as is seen for DNA run with the commonly used reorientation angle of 120 degrees. With DNA fragments of several hundred kb and smaller, higher field strengths may be used, resulting in still greater increases in separation speed. The conditions described allow DNA from large insert bacterial clones, such as those using cosmid, Fosmid, P1, bacterial artificial chromosome (BAC), or P1-derived artificial chromosome (PAC) vectors, to be prepared, digested and analyzed on gels within a single working day.     

Elastic Moduli of Collagen Gels Can Be Predicted from Two-Dimensional Confocal Microscopy

We quantitatively compare data obtained from imaging two-dimensional slices of three-dimensional unlabeled and fluorescently labeled collagen gels with confocal reflectance microscopy (CRM) and/or confocal fluorescence microscopy (CFM). Different network structures are obtained by assembling the gels over a range of concentrations at various temperatures. Comparison between CRM and CFM shows that the techniques are not equally sensitive to details of network structure, with CFM displaying higher fidelity in imaging fibers parallel to the optical axis. Comparison of CRM of plain and labeled collagen gels shows that labeling itself induces changes in gel structure, chiefly through inhibition of fibril bundling. Despite these differences, image analyses carried out on two-dimensional CFM and CRM slices of collagen gels reveal identical trends in structural parameters as a function of collagen concentration and gelation temperature. Fibril diameter approximated from either CRM or CFM is in good accord with that determined via electron microscopy. Two-dimensional CRM images are used to show that semiflexible polymer theory can relate network structural properties to elastic modulus successfully. For networks containing bundled fibrils, it is shown that average structural diameter, rather than fibril diameter, is the length scale that sets the magnitude of the gel elastic modulus.