Open Access
1 September 2008 Serum protein profile study of normal and cervical cancer subjects by high performance liquid chromatography with laser-induced fluorescence
Sujatha N. Raja, Lavanya Rai M.D., Prathap Kumar, Krishna K. Mahato, V. B. Kartha, Santhosh Chidangil
Author Affiliations +
Abstract
High performance liquid chromatography with high sensitivity laser-induced fluorescence detection is used to study the protein profiles of serum samples from healthy volunteers and cervical cancer subjects. The protein profiles are subjected to principal component analysis (PCA). PCA shows that the large number of chromatograms of a given class of serum samples—say normal/malignant—can be expressed in terms of a small number of factors (principal components). Three parameters—scores of the factors, squared residuals, and Mahalanobis distance—are derived from PCA. The parameters are observed to have a narrow range for protein profiles of standard calibration sets formed from groups of clinically confirmed normal/malignant classes. Limit tests using match/no match of the parameters of any test sample with parameters derived for the standard calibration sets give very good discrimination between malignant and normal samples with high sensitivity (~100%) and specificity (~94%).

1.

Introduction

Cervical cancer is one of the most common female cancers, second only to breast cancer,1 with a combined worldwide incidence of almost half a million new cases annually. It is the most common cancer of the female genital tract in India, with approximately 100 000 new cases occurring each year. This accounts for about 20% of all new cases diagnosed worldwide annually.2 Squamous cell carcinoma is the most common malignant cervical tumor, but the incidence of adenocarcinomas has been rising during the past few decades.3 Cervical cancer progresses in stages from dysplasia through CIN I, CIN II, CIN III, and finally to invasive cervical carcinoma.4, 5, 6, 7 At present, the Pap smear is the standard screening technique, and this is said to be responsible for a 70% decrease of cervical cancer deaths in advanced countries. However, false negative results are reported to be high (about 15 to 25%) for early detection.8 Over the last several years, Raman spectroscopy9, 10, 11, 12 and fluorescence spectroscopy 13, 14, 15, 16, 17, 18, 19, 20, 21, 22 of tissues have been developed as useful techniques for the discrimination of normal/malignant conditions and early detection of cervical malignancy. These optical methods are fast, highly objective, can be applied repeatedly in vitro/in vivo, and give direct information on the biochemical changes as the cells become malignant.

Many molecular and functional changes precede and continue to take place at the cellular level in neoplasia before morphologic or other alterations happen.4 Many of the biochemical changes, such as immune responses and production of degradative enzymes/fetal antigens, will be reflected in blood circulation.23 It is also known that in many epithelial cancers, even in the early stages, some of the abnormal cells escape into the blood stream and migrate to different locations. Clinically useful cancer signature materials (biomarkers) should thus be observable in an easily accessible body fluid such as serum. It is therefore highly likely that the protein profile of serum samples can possibly provide a diagnostic tool for early detection of cervical cancer. Proteomics is increasingly recognized as a very efficient area for possible early diagnosis of many diseases, like cancer and hypertension, which remain clinically silent for long periods.24, 25

Though laser-induced fluorescence and Raman spectroscopy of tissue also have been developed for discrimination between normal and malignant conditions, protein profiling has some advantages over these spectroscopic methods. The spectroscopic methods are applied to tissue samples, which are always heterogeneous, and mistakes can happen by placing the probe at the wrong site. Biopsy is done by physical examination (colposcopy) and errors like “past pointing” are possible so that you may not be getting the sample from the correct location. Protein profiling is done on homogeneous samples (serum), sampling is easy, and the sample is always representative of the subject irrespective of from where it is obtained. Though laser-induced fluorescence (LIF) and Raman are quite fast (a few minutes per subject) compared to high performance liquid chromatography (HPLC) protein profiling (30to40min) , a major advantage of the HPLC-protein profile method is that it gives much more information on changes that take place during the various stages of the disease. LIF gives broad structureless spectra, and only limited information can be derived on biochemical changes that take place during the progression of the disease. Though Raman spectra is better in this context, we still cannot get detailed information, since all biopolymers of a given class (say, proteins) will have very similar Raman spectra and it is very difficult to judge whether changes have taken place in protein composition, etc.

Protein profiling and proteomics have been suggested for early detection of diseases that remain clinically silent for a long period of time.26, 27, 28, 29 Proteomics is the study of functional genomics at the protein level. Several techniques like 2-D gel electrophoresis, SELDI, MALDI-TOF, and protein chips are available in proteomics.24, 25, 29 However, all these methods have several disadvantages compared to the HPLC-based technique.30

High performance liquid chromatography combined with laser-induced fluorescence detection (HPLC-LIF) is a highly sensitive technique for separation and detection of complex mixtures of proteins or other biomolecules,26, 27, 28, 30 even at subfemto mole levels. We have, in the present work, employed HPLC-LIF for protein profile analysis of serum samples from normal healthy volunteers and cervical cancer patients. The results show that there are noticeable differences between the two types of serum samples, even at a visual inspection level. We have applied, to the best of our knowledge for the first time, the method of principal component analysis (PCA) to the protein profiles of serum samples. This study has shown that serum protein profiles provide an objective method for the diagnosis of cervical cancer. The results of our studies are presented and discussed.

2.

Materials and Methods

2.1.

High Performance Liquid Chromatography Laser-Induced Fluorescence System

An HP 1100 gradient HPLC system with G1322A degasser, G1311A pump, and a manual injector (model number 7725, Rheodyne, Perkin Elmer, Massachusettes, USA) coupled to a Vydac 219TP52 biphenyl reversed phase narrow bore column (diphenyl, 2.1×250mm , 5μm , 300Å ), California, USA, was used for the separation of proteins. The effluent from the column was sent into a capillary flow cell fabricated in our laboratory. The cell is made of quartz capillary ( 75-μm I.D., 300-μm O.D. Hewlett Packard, G1600-64311) connected to the column with appropriate sleeves (Upchurch Scientific, F-130X, Washington, USA). It is mounted on a precision mount for accurate positioning and alignment for excitation and collection of fluorescence. Fluorescence excitation of proteins in the sample was achieved by illumination with the 257-nm laser emission (15mW) from a frequency doubled Ar+ laser (Innova 90C FreD, Coherent, California, USA), focused onto the capillary cell. The fluorescence was collected and focused onto the entrance slit of a monochromator (Jobin Yvon DH10 SPEX, New Jersey, USA), set at 340nm . The fluorescence was detected by a photomultiplier (Hamamatsu R 453, New Jersey, USA) operated at 850V , coupled through a preamplifier (EG&G model 5113, Maryland, USA) to a lock-in amplifier (EG&G model 7265). The fluorescence was chopped with an EG&G model 651 chopper at the entrance slit at 20Hz , for lock-in detection.

2.2.

Sampling and Sample Preparation

Normal samples of serum were collected from volunteers who were judged to be clinically normal (with respect to cervical cancer), and age matched, as far as possible. All normal subjects did not have any diagnosed disease. The age of the normal subjects varied from 22 to 77. Malignant samples were collected from subjects diagnosed with cervical cancer at the Department of Obstetrics and Gynecology, Kasturba Hospital, Manipal, India. Samples from 25 normal and 33 malignant subjects were analyzed. All samples were collected after informed consent. Ethical clearance of the Institute Ethical Committee was obtained for this work.

Blood samples were transported to the laboratory immediately after collection. The samples were stored at room temperature in an upright position for about 30min . The separated liquid portion was centrifuged at 3000rpm for 5min . The serum thus obtained was subjected to HPLC immediately. If storage was necessary, serum samples were stored at 80°C in the deep freezer. They were passively thawed to room temperature just before use. It was seen that the chromatograms of normal samples did not show any changes, even when recorded after several weeks of storage in the deep freeze.28 Table 1 gives details of the samples.

Table 1

Samples taken for analysis.

SamplenumberAgeClinicalcondition
1 to 2522 to 77Normal
26 to 3737 to 62Stage 2B
38 to 5437 to 72Stage 3B
5555Stage 2A
56,5745,50Stage 1
5860Stage 4A

The samples were diluted 100500 times with HPLC grade water. Fiftyμl of the diluted sample was injected into the HPLC system, which was fitted with a 20-μl loop. The sample was then eluted under a gradient run with [water, 0.1% v/v trifluoro acetic acid (TFA)] 70 to 40%, [acetonitrile (ACN), 0.1% v/v TFA], 30 to 60% in 0to60min , followed by 60 to 100% in 1min . The 100% ACN run was continued for another 14min , taking the total run time to 75min . After each run, the column was regenerated with HPLC grade water 0.1% TFA for 15min .

2.3.

Data Analysis

To facilitate intercomparison of protein profiles recorded over several months, a rigorous protocol was followed for data analysis. All the chromatograms were subjected to a background correction with a polynomial fit to remove the background fluorescence. Background in the chromatographic runs comes from several sources like laser light scattered from the walls of the capillary, Rayleigh scattering from the eluent, PMT dark count, fluorescence of the acetonitrile-TFA complex on excitation with 257nm , fluorescence of biomolecules that are not effectively separated by the biphenyl column, etc. All these contributions remain more or less the same from run to run. The background thus has the same shape from run to run, starting with about 0.25V (30% acetonitrile-TFA) in our case and increasing to about 0.4V at the end of the run (60% acetonitrile-TFA). The background variation is a continuous function (no sudden discontinuities) and so can be expressed with a polynomial of the form, background=a0+a1t+a2t2+a3t3+higher terms (negligible), where t is the time of observation (retention time). We can give four or more points on the chromatogram, at places free from sample contribution. These points can then be fitted to the polynomial to derive the constant coefficients a0 , a1 , etc. The background can then be calculated at any point of time and subtracted from the total signal at that time to give the background-free chromatogram.

Minor shifts in peak positions from run to run were corrected by calibration of all the chromatograms along the time scale, using mean values of protein peaks of species like transferrin, human serum albumin (HSA), creatin kinase, etc., common to all samples. Presence of these species in all samples provided convenient internal standards for calibration of all chromatograms to the same time scale across the entire range. The chromatograms were then normalized with respect to an HSA peak at 1666s . Because of the relatively constant intensity of this peak in all samples, it has been used for normalization of all chromatograms for comparison purposes.

To see whether information could be obtained on the changes in serum profiles, in a first step, analysis of the HPLC peaks were attempted using curve resolution techniques. For this, band shapes were decided using single peaks fit with different types of functions (Gaussian, Lorentzian, Voigt, etc.), and the function that gave the best fit for the single peaks (in the present case, Gaussian) was used for the other regions also. Difference chromatograms were also computed to identify any significant differences in the protein profiles of different classes of samples.

For precise classification of serum protein profiles by an objective mathematical model, the chromatograms of a sufficiently large number of samples (15 to 30) were subjected to PCA (PLS PLUS/IQ software, Galactic Corporation, Salem, New Hampshire). In our method of PCA, the mean of all samples in the dataset is first formed. The differences of each sample from this mean are calculated to give the variations of each sample from the mean. With n samples each having p data points, we thus get an n×p matrix of these variations. Because all the samples contain more or less the same components (HSA, transferrin, immunoglobulins, etc.), the large amount of data can be represented by a much smaller set of component peaks, their contributions to the chromatogram varying from sample to sample, depending on their concentrations. In matrix language, this implies that the n×p matrix of variations discussed previously is highly redundant. It will have only a few nonzero eigenvectors (principal components), and the eigenvalues will rapidly come down to almost zero after the first few. Solving the eigenvalue eigenvector problem gives us the principal components (factors), percent variance (contribution of the factors to the variations in the dataset), and scores of factors for each sample. The scores for a given sample correspond to the contribution of each principal component to the variation of that sample from the mean. It is therefore possible to simulate the original chromatogram of any sample by multiplying the eigenvectors with their respective scores for that sample and adding these products to the mean of the dataset.

Runs were made with 20, 15, 12, 9, and 7 factors. Seven factors were found to contribute to more than 95% of variance for the combined data. An analysis of the combined data, though it gives very good discrimination between normal and malignant cases (see below), has several drawbacks. To overcome these draw backs, PCA was finally done with protein profiles of standard sets for routine diagnostic applications, discussed in detail later.

3.

Results and Discussion

Figure 1 shows typical protein profiles of serum samples from subjects judged as clinically “normal.” A visual inspection of all normal samples shows more or less the same pattern, irrespective of age, physiological state (menopause, pregnancy, married, with/without children, etc.), life style, or social status.

Fig. 1

Protein profiles of the serum samples of healthy volunteers. Inset: top is 50years without background correction; below is after background correction. Thick gray line is 50years ; thin black line is 23years ; and dotted line is 45years .

054062_1_061805jbo1.jpg

Figure 2 shows the results for typical cases of malignant subjects. The malignant samples show noticeable changes from normal, in the peaks around 1200to1700-s region. The main differences are 1. a shoulder is observed for 1484-s peak in some of the cases, and 2. new peaks are observed in the 1200to1300-s region in some cases. These differences between normal and other classes of chromatograms can be seen better by looking at difference spectra, discussed later. From the prior discussion it is seen that even a routine protein profiling of serum may give an indication of possible cervical malignancy.

Fig. 2

Protein profiles of the serum samples of cervical cancer patients: thick gray line is stage 2, 60years ; thin black line is stage 3, 55years ; and dotted line is stage 4, 50years .

054062_1_061805jbo2.jpg

Many groups31, 32, 33 have used PCA scores of Raman/fluorescence spectra of normal and malignant tissue samples combined for discrimination between normal and malignant conditions. We have shown9, 11, 18 that a better approach for classification of tissue spectra as belonging to normal, inflammatory, premalignant, or malignant is to use the technique of matching several parameters from PCA for test samples to corresponding parameters of certified calibration sets from each class of samples. To decide whether such a situation exists for protein profiles by HPLC also, and to optimize conditions of PCA for discrimination, PCA was run first with all the samples combined, irrespective of whether they belong to the normal or malignant group. The PCA analysis was performed by taking 12 factors, as mentioned earlier. Figure 3 shows the variations in eigenvalues and percent variance with increasing factor numbers. From Fig. 3 it is clear that seven factors are adequate to represent the data, since they contribute up to more than 95% of the total variance in the dataset.

Fig. 3

Eigenvalues and total percentage variance for serum protein profiles.

054062_1_061805jbo3.jpg

The 58 samples shown in Table 1 were used in this PCA. The sample number versus score plot is shown in Fig. 4 for scores of factor 1, which is contributing 63% to the variance. It is seen that the normal and malignant groups form clusters falling in different ranges of factor 1 scores except for very few exceptions. Scores of 22 out of 25 normal samples lie on the positive side of the plot and three samples are lying on the negative score value. Similarly, 70% of the malignant samples have their scores on the negative side of the plot. So it is clear from the plot that the score values can discriminate to a very good extent between normal and malignant samples.

Fig. 4

Scores of factor 1 verses sample number.

054062_1_061805jbo4.jpg

It may be noted that PCA determines the factors necessary to express the data to a desired level of accuracy. It is therefore possible that factors other than the first or a combination of two/three factors may give better discrimination. However, in the present case this has not happened. (Data are not shown.) We have done cluster analysis of serum samples also.34 The plot of scores of PCA is similar to cluster analysis. But in our method of PCA, we can simulate a chromatogram of any sample (e.g., stage 1) with factors of any calibration set (e.g., stage 2), thus accounting for all stage 2 components in the stage 1 sample. If we now subtract the simulated chromatogram of the sample from its actual chromatogram, we can get more correct information, even on small differences between the two, since the simulation will account only for amounts of components that can be generated by factors of the stage 2 calibration set. These will be canceled on subtraction, while all peaks not produced by simulation will show up.

The simple approach of using scores for discrimination of sample types has several disadvantages. To begin with, it may not discriminate very well between normal and premalignant types (or between premalignant and malignant classes), because of possible closer similarities between chromatograms of such two classes of samples. Second, in a PCA with all types of samples, the factors have to account for all chromatograms, and so even factor number 1 may contain contributions from all types: normal, premalignant, and malignant. The discrimination between types will be thus diluted considerably. Third, in a PCA with all types combined, the results can become weighted more strongly toward that class that has more representation in it, and so may not give a correct picture for the other types of samples.

In view of this, we have developed a more reliable approach for diagnosis, combining all the information that is available from a PCA. The third property mentioned before tells us that if PCA is done with only one class of samples, say normal, then the factors will be weighted highly for that class and will correctly account for only protein profiles of that class. Profiles from any other class of samples will be rejected as not belonging to that class with a high degree of probability.

A much more important aspect of such an approach is that, like in any analytical technique where standards with calibration curves are used for routine analysis, profiles of a set of clinically/pathologically diagnosed samples can be used as a standard calibration set. This standard calibration set can be subjected to PCA to derive parameters that will be highly characteristic for samples of that type. Any test sample can then be added to the set and the corresponding parameters for the test sample can be compared to the mean parameters for the set to decide whether the test sample belongs to that set, and if so, with what statistical probability.

A third advantage of this method is that sets of standard calibration profiles can be prepared for each class—normal, premalignant (CIN I, CIN II, CIN III, etc.), and malignant (different stages). Any test sample can be compared against each of these calibration sets and a decision can be made more accurately, depending on with which set it matches best.

It should also be mentioned here that calibration sets can be prepared in a main hospital, where qualified clinicians/pathologists are available and a large number of subjects may be examined in a reasonable period, say, in a few months to a year. These calibration sets can be supplied to any small hospital/clinic, and clinicians there can carry out objective diagnosis on their own with a blood sample routinely collected.

The possible errors that may arise from visual examination (fatigue factor in Pap smear and inexperience in colposcopy) are reduced in optical methods entirely based on recording of an optical signal and analysis of the data with standard mathematical procedures by a computer. No visual decision making is involved, and the system (HPLC instrument system plus a computer) is completely blind as to what sample is given for analysis.

Our diagnostic approach thus consists of the following steps. Record protein profile (chromatogram) from a reasonably large number of subjects who are clinically/pathologically diagnosed as normal, premalignant, or malignant. Take all the profiles from samples, clinically/pathologically diagnosed as belonging to one class, and do a PCA. Determine the number of significant factors required to represent that class. Check and remove any “outliers.” (In practice, a very small number, one or two, in a given class may stand out as not belonging to that class, for reasons such as instrument malfunction, ill-defined disease condition, etc.). From the set, randomly select 15 to 20 profiles to form the standard calibration set. Do PCA with this standard set to determine the scores of significant factors and other parameters. Test every member of the set by rotating them out one at a time for membership of the set within a desired range of standard deviation of the chosen parameters. Reject any that do not match to form the final calibration set. Do this for every class, where sufficient numbers of clinically/pathologically certified samples are available. It is helpful to remember that, since usually a maximum of 4 to 6 factors are quite sufficient to determine more than 95% of the variations from the mean for a given class, if a minimum of 10 to 12 chromatograms are available in a class, they are enough to form a calibration set. Once calibration sets are prepared for the different classes of samples, any test sample profile can be matched against each of the calibration sets. The sample is diagnosed as belonging to that class with which it matches best.

In addition to the scores and spectral residual, a statistical measure known as the Mahalanobis distance is also used for the match-no match test described above. This is calculated as the distance of the test sample point as measured from the mean of all the remaining points in the class. The distance is scaled in units of standard deviation for the range of variation in the class in all dimensions, and then used to assign a probability weight to the sample in terms of standard deviation.

Any sample that lies outside a desired range of standard deviation from the mean can be considered to be out of the group. The range can be decided by the clinician who can place a cut-off range. In Tables 2 and 3 a sample is considered “no match” if the Mahalanobis distance M is greater than 2.5. Since M is in units of standard deviation, a value of > 2.5 gives a probability of less than 1% of that sample belonging to the set with which it is matched. If the clinician wants a less rigorous cut-off, a higher value (e.g., 3) can be given that gives a probability of less than 0.1% for that sample to be out.

The Mahalanobis distance matrix equation is

M=(SS(n1)),
where M is an f×f Mahalanobis matrix, S is an n×f matrix of training sample PCA scores, n is the number of samples, and f is the number of PCA factors. The Mahalanobis distance D for a test sample is then given by
D2=(Sunk)M1(Sunk),
where D2 is the square of the Mahalanobis distance in terms of standard deviations of the set, represented by M .

There is one problem regarding the discrimination purely based on PCA scores. It is that any “impurities” that are in the unknown spectra, but were not present in the training samples, will not appear in the score calculations. To overcome this problem and get more specific diagnoses, residuals can be included in the vectors. The purpose of the protein profiling technique thus is to classify any test sample as belonging to only one of the well-defined sets, without or with limited prior knowledge. For this, standard calibration models are built from profiles of sets of samples, clinically/pathologically certified as normal, premalignant, and malignant-different stages. The parameters defined earlier—Mahalanobis distance, scores of factors, and squared residuals—are then used to discriminate between the sample types. This gives a very sensitive discrimination and shows how well an unknown test sample matches or does not match with any given calibration set.

It is quite reasonable to expect (and this has been verified experimentally) that all normal samples will more or less give the same profile. On the other hand, it is possible that malignant samples in different stages (stages 2, 3, etc.) may differ slightly from each other depending on the stage of the disease. However, it is highly probable that all such samples will be quite different from normal. Also, it is very likely that they may have some common features in their profiles because of conditions common in malignancy. Therefore, if all malignant samples are clubbed together as a class separate from normal, it is possible that in PCA some factors will contribute to variations common to all of them, while a few other factors may take care of intragroup variations. As mentioned earlier, it is also to be noted that in a mixed sample lot, the factors will contribute to common variations with weights corresponding to number of members in the subgroups. Hence if we have a sufficiently large number of samples of different stages, the factors will adjust to represent all members of the set. To test this hypothesis in the present work, we first formed standard model calibration sets of normal and malignant groups by randomly selecting 15 samples respectively from the normal and malignant groups, disregarding the stage of malignancy.

PCA was then done first with the normal calibration set. All normal and malignant samples were then matched with this calibration set. The 15 normal samples forming the standard set were tested for match/no match by rotating them out one by one (retrospectively), and the remaining ten normal samples were tested against the standard set as unknown samples (prospectively). It is seen that all the normal samples matched with the normal standard set. This corresponds to specificity [true negative/(true negative+false positive)] of 100%. When the 33 malignant samples were matched against the normal standard set, only two malignant samples (29 and 33) showed a match. The result is thus still better than that seen from the score versus sample number plot (Fig. 4), where nine malignant samples were in the range of the normal cluster. The match/no match test was then performed with the malignant standard set. It was seen that all normal samples did not match with the malignant set. Samples 29 and 33, which matched with normal, matched here also. Only two malignant samples, 44 and 53, did not match with the malignant set. The sensitivity (true positive/true positive+false negative) is thus about 94%. These results are shown in Table 2 .

Table 2

Discrimination analysis with normal and malignant standard set.

Discrimination analysis with normal standard set
SampleMatchM distanceLimit testsSpecresidual
1 to 25Yes0.63 to 2.3Pass1.83 to 4.03
26 to 28No5.97 to 8.07Fail7.70 to 10.26
29Yes1.51Pass2.98
30 to 32No3.96 to 5.34Fail6.02 to 7.52
33Yes2.12Pass3.63
34 to 58No2.58 to 14.36Fail3.61 to 19.45
Discrimination analysis with malignant standard set
SampleMatchM distanceLimit testsSpecresidual
1 to 25No2.86 to 11.95Fail4.00 to 18.16
26 to 43Yes0.62 to 2.27Pass0.33 to 5.12
44No3.07Fail5.89
45 to 52Yes0.74 to 1.48Pass0.88 to 4.01
53No10.82Fail17.34
54 to 58Yes1.04 to 2.13Pass1.12 to 5.12

3.1.

Difference Chromatograms

As mentioned earlier, to observe the changes from one stage to another, the difference chromatograms were calculated for serum samples from different stages. The changes from normal serum for the different types of samples are shown in Fig. 5 . The differences here are given by mean of calibration set of normal samples—chromatogram of mean of any given class simulated with the normal calibration set. In Fig. 5a, we have the mean of all normal samples as the test sample and the differences are practically zero (less than ±0.000 000 03), as it should be. The difference chromatograms for the different stages, on the other hand, show that there are noticeable differences of fairly large magnitude from the mean normal for all the malignant samples. The important point to note is that in addition to the large differences from the mean normal, the mean profiles for different stages of the disease showed appreciable differences between them also, as seen from the peaks marked with their retention time.

Fig. 5

Difference chromatograms: (a) mean of calibration set of normal samples—chromatogram of mean of a normal simulated calibration set; (b) mean of calibration set of normal samples—chromatogram of mean of stage 1 simulated calibration set; (c) mean of calibration set of normal samples—chromatogram of mean of stage 2 simulated calibration set; (d) mean of calibration set of normal samples—chromatogram of mean of stage 3 simulated calibration set; and (e) mean of calibration set of normal samples—chromatogram of a stage 4 simulated calibration set.

054062_1_061805jbo5.jpg

3.2.

Discrimination between the Stages of Samples

To see whether discrimination of different stages of samples is possible, we carried out match/no match PCA with calibration sets of stages 2 and 3 samples for which a reasonable number (13 and 17, respectively) of pathologically certified samples were available. PCA analysis was carried out using two samples of stage 1, 13 of stage 2, 17 of stage 3, and one of stage 4. With ten samples of stage 2 forming a calibration set, all stage 2 samples and also the stage 4 sample showed a match. All stage 3 and stage 1 samples showed no match. Out of 17 stage 3 samples, 12 samples were taken to form a calibration set. It is found that, except for two samples (29 and 33), all 11 samples from stage 2 did not match with the calibration set. Three samples (42, 44, and 48) from the stage 3 group were also not matching with the calibration set of stage 3. This preliminary study shows that it may be possible to discriminate between different stages of the disease by protein profile analysis. Table 3 shows the match/no match results. Figure 6 shows the spectral residual versus M-distance plot in this case. It can be seen that even with two parameters, the residual and Mahalanobis distance, there is a clear separation of stage 2 and stage 3 samples when either of these calibration sets is used to classify samples. The samples of the calibration set all cluster together, while samples from the other stages are widely separated from the calibration set.

Fig. 6

Discriminate analysis with (a) standard calibration set of stage 2 and (b) stage 3 cervical cancer serum samples.

054062_1_061805jbo6.jpg

Table 3

Discrimination analysis with stage 2 and stage 3 malignant standard set.

Discrimination with stage 2 malignant standard set
SampleM distanceLimit testsSpec residual
26 to 370.72 to 1.11Pass0.02 to 1.03
38 to 543.87 to 37.91Fail1.76 to 13.56
551.08Pass0.25
56, 578.80 to 14.61Fail4.0 to 5.9
580.89Pass0.51
Discrimination with stage 3 malignant standard set
SampleM distanceLimit testsSpec residual
26 to 282.95 to 16.34Fail3.98 to 17.22
291.56Pass4.61
30 to 322.96 to 4.21Fail1.74 to 5.10
330.79Pass5.34
34 to 372.92 to 10.61Fail0.67 to 11.67
38 to 410.79 to 1.17Pass0.62 to 10.50
429.49Fail1.27
430.79Pass3.91
442.29Fail1.52
45 to 470.83 to 1.18Pass0.61 to 5.81
484.46Fail1.93
49 to 540.75 to 1.65Pass1.07 to 19.06
55 to 582.77 to 10.61Fail3.97 to 11.67

3.3.

Curve Resolution Studies

To see whether additional information can be obtained on the different stages of cervical cancer serum, we carried out curve resolution studies of the regions involved. The peaks were resolved using the Gaussian function. The two unresolved peaks from 1300 to 1800 regions were resolved by this method. In stage 4 samples, additional peaks were observed. These curve fit results are shown in Fig. 7a for the normal sample, and the results for the samples in different stages are given in Table 4 . The data in Table 4 illustrate the very good reproducibility in the different chromatograms for similar peaks, as indicated by the peak positions, intensities (heights), and half widths for many of the peaks. It may be noted that as a result of the calibration procedure, the variation in peak positions from run to run is minimized. For example, for the 2440-s region peaks, the peak positions for the four different types of samples vary only about ± 1s , less than 0.04%. The area under the curve corresponding to the third peak has decreased continuously in going from normal to stage 4 (Table 4).

Fig. 7

Curve fitting for the protein profile of the serum samples: (a) region 1300to1800s and (b) region 2300to2600s . Original chromatogram (in red) and fitted curve (in green). (Color online only)

054062_1_061805jbo7.jpg

Table 4

Curve fitting for values for the protein profile of the serum samples. Peak numbers 1 to 4 (region 1300to1800s ) and 5, 6, and 7 (region 2300to2600s ).

Type ofsamplePeaknumberCenterHeightWidthArea
Normal11345.91.2425.834.2
21484.51.927.257.9
31530.16.327.9187.9
41667.20.8626.724.3
52350.60.4345.520.9
62438.30.7363.149.7
72566.10.2966.220.6
Stage 111343.91.0127.129.3
21482.11.933.0170.2
31529.55.529.7172.7
41666.10.928.526.4
52348.80.629.919.1
62439.10.738.530.53
72571.20.339.3512.92
Stage 211345.90.9925.524.9
21484.82.0428.561.6
31530.25.0628.3152.2
41667.30.8625.623.0
52349.50.542.3222.13
62438.80.6459.4341.01
72565.70.3268.3723.43
Stage 311346.60.9626.227.02
21485.02.0227.358.4
31530.44.6730.1150.0
41667.70.8727.625.5
52350.40.4545.4721.9
62438.30.7265.049.9
72567.40.3367.023.5
Stage 411348.21.1223.8328.5
1a1454.30.1931.926.43
21485.21.1625.8932.04
31529.43.2525.9289.7
3a1611.90.0778.346.26
41668.90.925.9824.9
52348.40.7534.928.01
62438.70.5543.3725.6
72570.80.3242.2814.5

Similar curve fit results for the other regions (2300to2600s) are shown in Fig. 7b for the normal sample, and the results for the different stages are given in Table 4. These peaks also showed noticeable differences between the different stages. From Table 4, for the peaks 5, 6 and 7, it is clearly seen that the half-width and peak intensities change for many component proteins, even when their retention times are relatively unchanged. These variations provide the basis for discrimination in PCA. In our method, we are doing the decision making based on principal component analysis (PCA) and match/no match of the parameters derived from PCA. The elution time, peak height, and width are built into the factors derived from the PCA. The parameters used include scores of factors, and residuals=(observed chromatogram-chromatogram simulated using PCA factors and scores).2 The elution time, peak height, and width automatically come into the scores of factors and simulated chromatograms. It is obvious that all these data are thus used in the analysis and decision making.

In view of such differences between normal and different classes of serum, it may be tempting to use the corresponding proteins as tumor markers for diagnosis, prognosis, and follow up in therapy. However, it is well recognized4, 35 that many tumor markers are useful only as indicators, and do not provide unambiguous diagnosis because of the possibility that they may also be present under other conditions, like pregnancy. On the other hand, the present method, which makes use of the entire protein profile, will be much more reliable since it uses all observed variations in the sample.

The minor discrepancies observed between the pathological and HPLC-LIF protein profiling results can arise from several factors. These include: (1) experimental variations (errors) in pathology and HPLC-LIF, (2) data processing errors, and (3) inherent variations from sample to sample.

The experimental variations in pathology can come from sampling errors like “past pointing” in colposcopy,36 fatigue factor in examining large numbers of slides, or the inexperience of pathologist. But in the present work, both calibration standards and test samples of malignant conditions consist of samples diagnosed by biopsy and pathology. Since the pathology results are based on morphological changes, when they are diagnosed as malignant, it is highly unlikely they can be in error. Variations in HPLC-LIF can come from blood collection and serum preparation procedures, storage effects, contamination, or errors in HPLC runs. However, the fact that all the normal samples matched with one another very well, even though the samples were collected and run over a period of more than a year, and different samples were stored for periods ranging from a few hours to weeks, suggest that experimental variations in HPLC-LIF may not be responsible for any discrepancies.

In many developing countries, where advanced medical facilities are restricted to major towns/hospitals, what is required is a system usable as a screening method for susceptible populations. The system should be able to perform as a diagnostic tool for suspect cases, using readily available samples, the quality of which will not depend on the personal experience of the examining physician. It should also be able to do follow-up without the need for repeated biopsy and pathology. The system has to be fairly rugged, relatively fast (reasonable number of samples per day), usable by a clinical technician (to be used in a single clinician/surgeon’s office or in a small hospital in a village), and should have specificity and sensitivity comparable to currently available techniques.

We feel that the present technique meets all these requirements and can be used in a routine manner in small clinics and hospitals without the need for qualified pathologists to examine cytological smear and biopsy tissue samples. Also the technique, highly objective and operator-independent, will serve as a very good complimentary method for conventional screening techniques, reducing errors due to “fatigue factor,” sampling errors, any lack of experience of the pathologist, and consequent subjective diagnosis. The protein profile study of the cervical and normal serum samples contains much information on the structural or conformational changes occurring within the proteins. We have seen that significant changes occur between the protein profiles of normal and malignant samples. These changes may be due to changes in either structure or conformation produced during the process of carcinogenesis. The analysis of 58 normal and malignant samples using PCA was found to give very good results, with 100% specificity and 94% sensitivity.

4.

Conclusions

Analysis of protein profiles of serum by the HPLC-LIF technique provides a method for screening for early detection of cancer. The method is minimally invasive (only routine blood sample required), relatively fast ( <1h /sample), highly objective (no visual observation or personal judgment), and can be used even by a technician in a small clinic/hospital. It can be used for initial screening, follow-up in therapy, and early detection of any recurrence, and can be applied repeatedly, since no biopsy is involved.

Acknowledgment

The work was done under the project “Study of the protein profile of the clinical samples for the early diagnosis of female cancers,” Department of Science and Technology, Government of India, Project number SR/S2/LOP/05/2003.

References

1. 

E. L. Franco, N. F. Schlecht, and D. Saslow, “The epidemiology of cervical cancer,” Cancer J., 9 (5), 348 –359 (2003). 1528-9117 Google Scholar

2. 

V. Shanta, S. Krishnamurthi, C. K. Gajalakshmi, R. Swaminathan, and K. Ravichandran, “Epidemiology of cancer of the cervix: global and national perspective,” J. Indian Med. Assoc., 98 (2), 49 –52 (2000). 0019-5847 Google Scholar

3. 

S. G. Silverberg and O. B. Ioffe, “Pathology of cervical cancer,” Cancer J., 9 (5), 335 –347 (2003). 1528-9117 Google Scholar

4. 

V. Kumar, R. S. Cotran, and S. L. Robbins, “Neoplasia,” Basic Pathology, Prism Books, India (1992). Google Scholar

5. 

A. Giuliano, “Cervical carcinogenesis: The role of co-factors and generation of reactive oxygen species,” Salud Publica Mex, 45 354 –360 (2003). 0036-3634 Google Scholar

6. 

B. C. Das, V. Gopalkrishna, S. Hedau, and S. Katiyar, “Cancer of the uterine cervix and human papilloma virus infection,” Curr. Sci., 78 52 –63 (2000). 0011-3891 Google Scholar

7. 

D. C. Walker, B. H. Brown, A. D. Blackett, J. Tidy, and R. H. Smallwood, “A study of the morphological parameters of cervical squamous epithelium,” Physiol. Meas, 24 (1), 121 –135 (2003). 0967-3334 Google Scholar

8. 

J. Monsonego, “Global challenge of cervical cancer screening and prevention,” Eur. J. Gynaecol. Oncol., 21 (6), 533 –539 (2000). 0392-2936 Google Scholar

9. 

C. M. Krishna, N. B. Prathima, B. M. Vadhiraja, R. Malini, D. J. Fernandes, M. S. Vidyasagar, and V. B. Kartha, “Raman spectroscopy studies for diagnosis of cancer in human uterine cervix,” Vib. Spectrosc., 41 136 –141 (2006). https://doi.org/10.1016/j.vibspec.2006.01.011 0924-2031 Google Scholar

10. 

K. Venkatakrishna, J. Kurien, K. M. Pai, C. M. Krishna, G. Ullas, and V. B. Kartha, “Optical pathology of oral tissue—A Raman spectroscopic diagnostic method,” Curr. Sci., 80 101 –105 (2001). 0011-3891 Google Scholar

11. 

M. V. P. Chowdary, K. K. Kumar, J. Kurian, S. Mathew, and C. M. Krishna, “Discrimination of normal, benign, malignant breast tissue by Raman spectroscopy,” Biopolymers, 83 (5), 556 –569 (2006). https://doi.org/10.1002/bip.20586 0006-3525 Google Scholar

12. 

M. D. Keller, E. M. Kanter, and A. Mahadevan-Jansen, “Raman spectroscopy for cancer diagnosis,” Spectrosc., 21 33 –41 (2006). Google Scholar

13. 

D. L. Heintzelman, U. Utzinger, H. Fuchs, A. Zuluaga, K. Gossage, “Optimal excitation wavelengths for in-vivo detection of oral neoplasia using fluorescence spectroscopy,” Photochem. Photobiol., 72 103 –113 (2000). https://doi.org/10.1562/0031-8655(2000)072<0103:OEWFIV>2.0.CO;2 0031-8655 Google Scholar

14. 

D. De Veld, G. M. Witjes, H. Sterenborg, and J. Roodenburg, “The status of in vivo autofluorescence spectroscopy and imaging for oral oncology,” Oral Oncol., 41 117 –131 (2005). https://doi.org/10.1016/j.oraloncology.2004.07.007 0964-1955 Google Scholar

15. 

S. Madhuri, N. Vengadesan, P. Aruna, D. Koteeswaran, P. Venkatesan, S. Ganesan, “Native fluorescence spectroscopy of blood plasma in the characterization of oral malignancy,” Photochem. Photobiol., 78 197 –204 (2003). https://doi.org/10.1562/0031-8655(2003)078<0197:NFSOBP>2.0.CO;2 0031-8655 Google Scholar

16. 

B. K. Manjunath, K. Jacob, L. Rao, C. M. Krishna, M. S. Chidananda, K. Venkatakrishna, and V. B. Kartha, “Autofluorescence of oral tissue for optical pathology in oral malignancy,” J. Photochem. Photobiol., B, 73 49 –58 (2004). https://doi.org/10.1016/j.jphotobiol.2003.09.004 1011-1344 Google Scholar

17. 

G. S. Nayak, D. S. Kamath, M. K. Pai, A. Sarkar, R. Satadru, J. Kurien, L. D’Almeida, B. R. Krishnanand, C. Santhosh, V. B. Kartha, and K. K. Mahato, “Principal component analysis (PCA) and artificial neural network (ANN) analysis of oral tissue fluorescence spectra. Classification of normal premalignant and malignant pathological conditions,” Biopolymers, 82 152 –166 (2006). https://doi.org/10.1002/bip.20473 0006-3525 Google Scholar

18. 

M. S. Chidananda, K. Satyamoorthy, R. Lavanya, A. P. Manjunath, and V. B. Kartha, “Optical diagnosis of cervical cancer by fluorescence spectroscopy technique,” Int. J. Cancer, 119 139 –145 (2006). https://doi.org/10.1002/ijc.21825 0020-7136 Google Scholar

19. 

C. Brookner, U. Utzinger, “Effects of biographical variables on cervical fluorescence emission spectra,” J. Biomed. Opt., 8 (3), 479 –483 (2003). https://doi.org/10.1117/1.1578642 1083-3668 Google Scholar

20. 

G. M. Palmer, P. J. Keely, T. M. Breslin, and N. Ramanujam, “Autofluorescence spectroscopy of normal and malignant human breast cell lines,” Photochem. Photobiol., 78 (5), 462 –469 (2003). https://doi.org/10.1562/0031-8655(2003)078<0462:ASONAM>2.0.CO;2 0031-8655 Google Scholar

21. 

N. Ramanujam, “Fluorescence spectroscopy of neoplastic and non-neoplastic tissues,” Neoplasia, 2 89 –117 (2000). https://doi.org/10.1038/sj.neo.7900077 1522-8002 Google Scholar

22. 

N. Ramanujam, M. F. Mitchell, A. Mahadevan, S. Thomsen, E. Silva, and R. Richards-Kortum, “Fluorescence spectroscopy: a diagnostic tool for cervical intraepithelial neoplasia (CIN),” Gynecol. Oncol., 52 31 –38 (1994). https://doi.org/10.1006/gyno.1994.1007 0090-8258 Google Scholar

23. 

J. Mendelson, “Principle of neoplasia,” Harrison’s Principles of Internal Medicine, 11th Ed.Mc Graw Hill, New York (1987). Google Scholar

24. 

C. P. Paweletz, “Rapid protein display profiling of cancer progression directly from human tissue using a protein biochip,” Drug Dev. Res., 49 34 –42 (2000). 0272-4391 Google Scholar

25. 

P. R. Shrinivas, S. Srivastava, S. Hanesh, and G. L. Wright, “Proteomics in early detection of cancer,” Clin. Chem., 47 1901 –1911 (2001). 0009-9147 Google Scholar

26. 

Sujatha, L. Rai, B. R. Krishnanand, K. K. Mahato, V. B. Kartha, and C. Santhosh, “Protein profile study of the cervical cancer using HPLC-LIF,” Proc. SPIE, 6092 60920V (2006). 0277-786X Google Scholar

27. 

S. K. Singh, R. L. Martis, Sujatha, R. A. Bhat, P. Kushtagi, L. Rai, V. B. Kartha, and C. Santhosh, “Protein profile study of clinical sample of ovarian cancer using high performance liquid chromatography-laser induced fluorescence (HPLC-LIF),” Proc. SPIE, 6092 60920N (2006). 0277-786X Google Scholar

28. 

K. Venkatakrishna, V. B. Kartha, M. K. Pai, C. M. Krishna, O. Ravikiran, J. Kurian, M. Alexander, and G. Ullas, “HPLC-LIF for early detection of oral cancer,” Curr. Sci., 84 551 –557 (2003). 0011-3891 Google Scholar

29. 

E. P. Diamandis, “Proteomic patterns in biological fluids: Do they represent the future of cancer diagnostics,” Clin. Chem., 49 (8), 1272 –1278 (2003). 0009-9147 Google Scholar

30. 

S. Menon, Sujatha, K. K. Kumar, J. Kurien, B. R. Krishnanand, K. K. Mahato, S. D. George, V. B. Kartha, and C. Santhosh, “Protein profile study of breast cancer tissues using HPLC-LIF—A pilot study,” Proc. SPIE, 6430 64300W1 (2007). 0277-786X Google Scholar

31. 

A. M. Jansen, M. F. Mitchell, N. Ramanujam, A. Malpica, S. Thomsen, U. Utzinger, and R. Richards-Kortum, “Near-infrared Raman spectroscopy for in vitro detection of cervical pre-cancers,” Photochem. Photobiol., 68 (1), 123 –132 (1998). https://doi.org/10.1562/0031-8655(1998)068<0123:NIRSFV>2.3.CO;2 0031-8655 Google Scholar

32. 

K. Chang, M. Y. Dawood, G. Staerkel, U. Utzinger, E. N. Atkinson, R. R. Richards-Kortum, and M. Follen, “Fluorescence spectroscopy for cervical precancer detection: Is there variance across the menstrual cycle,” J. Biomed. Opt., 7 (4), 595 –602 (2002). https://doi.org/10.1117/1.1509753 1083-3668 Google Scholar

33. 

C. G. Diana, D. Veld, M. Skurichina, M. J. H. Witjes, R. P. W. Duin, D. J. C. M. Sterenborg, W. M. Star, and J. L. N. Roodenberg, “Autofluorescence characteristics of healthy oral mucosa at different anatomical sites,” Lasers Surg. Med., 32 367 –376 (2003). https://doi.org/10.1002/lsm.10185 0196-8092 Google Scholar

34. 

C. Santhosh, “Serum protein profile study of clinical samples using high performance liquid chromatography laser induced fluorescence techniques,” SPIE Photonics West-2008. Technical Program, 117 (2008) Google Scholar

35. 

C. R. A. Booth and E. P. Diamandis, “Tumor marker,” Teitz Textbook of Clinical Chemistry, 745 –786 Wiley, New York (2004). Google Scholar

36. 

M. F. Mitchell, “Accuracy of colposcopy,” Clin. Consult. Obstetrics Gynecol., 6 70 –73 (1994). 1043-0660 Google Scholar
©(2008) Society of Photo-Optical Instrumentation Engineers (SPIE)
Sujatha N. Raja, Lavanya Rai M.D., Prathap Kumar, Krishna K. Mahato, V. B. Kartha, and Santhosh Chidangil "Serum protein profile study of normal and cervical cancer subjects by high performance liquid chromatography with laser-induced fluorescence," Journal of Biomedical Optics 13(5), 054062 (1 September 2008). https://doi.org/10.1117/1.2992166
Published: 1 September 2008
Lens.org Logo
CITATIONS
Cited by 16 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Proteins

Calibration

Principal component analysis

Cervical cancer

Mahalanobis distance

Statistical analysis

Laser induced fluorescence

Back to Top