1 May 2007 Sensitivity map of laser tweezers Raman spectroscopy for single-cell analysis of colorectal cancer
Author Affiliations +
Abstract
Raman spectroscopy on single, living epithelial cells captured in a laser trap is shown to have diagnostic power over colorectal cancer. This new single-cell technology comprises three major components: primary culture processing of human tissue samples to produce single-cell suspensions, Raman detection on singly trapped cells, and diagnoses of the cells by artificial neural network classifications. It is compared with DNA flow cytometry for similarities and differences. Its advantages over tissue Raman spectroscopy are also discussed. In the actual construction of a diagnostic model for colorectal cancer, real patient data were taken to generate a training set of 320 Raman spectra and a test set of 80. By incorporating outlier corrections to a conventional binary neural classifier, our network accomplished significantly better predictions than logistic regressions, with sensitivity improved from 77.5% to 86.3% and specificity improved from 81.3% to 86.3% for the training set and moderate improvements for the test set. Most important, the network approach enables a sensitivity map analysis to quantitate the relevance of each Raman band to the normal-to-cancer transform at the cell level. Our technique has direct clinic applications for diagnosing cancers and basic science potential in the study of cell dynamics of carcinogenesis.

1.

Introduction

Recently, laser tweezers and Raman spectroscopy (LTRS) were combined to diagnose single living human cells confined in a laser trap.1, 2 Both solid tumors surgically removed from human colorectal cancer sites1 and human lymphocytes2 were investigated. In the first case, primary culture protocols were followed to detach the cancer cells from their surroundings,1, 3 while in the second case, the healthy and transformed lymphocytes were already in single-cell form.

In contrast, DNA flow cytometry (FCM) is a currently adopted medical procedure to evaluate the malignancy of a tumor,4 in which single-cell suspensions are prepared, stained with a DNA-binding fluorophore, and analyzed by a flow cytometer. LTRS and DNA FCM carry several similarities. First, both are single-cell techniques, in which the cells are analyzed one by one; second, both share the same procedure for sample preparation (up to the staining); and third, in both cases, the cells are detected in an aqueous environment involving laser excitation and optical detection.

DNA FCM relies on ploidy information, such as irregularities of DNA contents, extra chromosomes, and increased percentage of S-phase cells undergoing DNA replication, to quantitate the aggressiveness of cancer. Thus, it cannot reveal the cell structural changes or the biochemical variations associated with the normal-to-cancer transition. Furthermore, a DNA FCM analysis typically requires 20,000 single cells in order to accumulate enough statistical significance5; therefore, such tests are unavailable to small amount samples from exfoliative cytology or biopsies on small tumors. Last, the accuracy of DNA FCM tests tends to fluctuate on several factors. A procedural guideline has been created to ensure quality control and reproducibility of the test results among different laboratories.5 According to the guideline, adequate neoplastic material in a specimen should be ensured (by other means) before a test can be invoked in order to justify the accuracy. For example, at least 20% tumor cells are necessary if tumor proliferation calculations are needed. Further, accidental mixing with normal tissues can deviate or may even invalidate the test results. The statistical nature of FCM (including conventional FCM) ensures that it is an ensemble averaging technique.

On the other hand, LTRS has several characteristics unavailable to DNA FCM. First, LTRS allows single living cells to be directly investigated without any staining, free of the artifacts introduced by the staining process. FCM is essentially a serial method with low efficiency of information, and difficulty in multicolor labeling prevents the simultaneous detection of multiple targets on the same cell, whereas LTRS collects the information about the whole biochemical composition of a cell and reports the contributions of vast chemical function groups in parallel. The abundance of information makes LTRS a true single-cell diagnosis technique, in which the decision is made on each cell, not the ensemble average over 20,000 cells. For the diagnosis of a patient, 20 cells are sufficient for an LTRS test.1 Therefore, small tumor biopsy, exfoliative cytology, and scrape and brush cytology can all be employed to collect materials enough for an LTRS run. Nevertheless, DNA FCM and LTRS are both achievements of laser technologies, and each has its own merits and should be used to supplement each other to assist doctors in the fight against cancers.

For two decades, Raman spectroscopy has demonstrated diagnostic power over many diseases, from various cancers, 6, 7, 8, 9 to atherosclerosis,10 to Alzheimer’s disease,11 etc. In the area of cancer studies, great efforts have been taken to obtain the Raman features of lesions from breasts,6, 7 bladders,8 skins,9 etc. Sophisticated diagnostic algorithms have been developed to differentiate, classify, and stage tumors. Over the years, progress in data collection schemes and spectral analysis methods has achieved high accuracy and pushed the detection limit to early stages.

As is well known, more than 85% of all cancers originate from the epithelia lining the internal surfaces of human organs. At early stages of cancer development, it is the surface tissue layer that carries diagnostic information, while the underlying bulk tissues are mostly irrelevant and even interfere with the detection. In addition, tissue components in the bulk are the major sources of fluorescence background that complicates the analysis of Raman spectra. In recent years, there is a trend to go from bulk-tissue spectroscopy to surface-tissue spectroscopy. Polarization gating has been exploited to selectively detect the photons singly backscattered from the surface layer within the depth of one scattering length.12 A natural step further would be to retrieve single cells from the epithelium by medical cytology means and study them in isolation from the bulks. This not only creates a disease-screening tool, but also allows investigation on the cell dynamics of carcinogenesis in a controlled environment. LTRS is a way of realizing these new ideas.1, 2 Unfortunately, current LTRS research mainly focuses on bacteria,13, 14 because the primary culture techniques are unfamiliar to the optical community and human cells are larger and harder to trap. Here we call for attention to LTRS research on cancers.

Successful Raman diagnosis also depends on spectral analysis. Visual inspection of Raman peaks for clues is subjective, correlations among multiple peaks are often hidden, and simple peak ratios inevitably overlook the overall spectral shape as well as fine spectral details. Sophisticated diagnostic models always resort to some kind of numerical algorithms. Principal component analysis (PCA) is commonly employed to reduce data dimensionality; linear regressions or logistic regressions are the prevailing algorithms for classification. The role of individual Raman peaks in the classification, i.e., which peak is important and which is not, is qualitatively argued (but not quantitatively evaluated) by visual inspection of the spectral patterns of the principal components (PCs) included in the regressions.1, 7 Recently, χ2 fittings of the tissue/cell spectra have been attempted, using the spectra of 10 major chemical constituents as the expansion basis.6, 15 However, this method seems oversimplifying, since a cell easily contains at least 8,000 different proteins and many other components. The incompleteness of the spectral basis set can impose a problem to the validity of the results. At the least, some fine structures of the Raman data could be missing in the base spectra and thus ignored in the fitting. In addition, the results can be prone to distortion due to the dominance of one or two strong Raman constituents.

Artificial neural network (ANN) is a nonlinear, multidimensional model capable of classifying spectral data.9, 16 The flexibility of ANN results in better discriminating power than any other regression models. But the nonparametric black-box nature of ANN has prevented its adoption compared with other explicit algorithms. To overcome this problem, Sigurdsson invoked sensitivity analysis in their ANN study of skin cancers,16 in which the contribution of each Raman frequency to the network output was numerically evaluated by a function called a sensitivity map.

In this paper, we report an LTRS investigation on colorectal cancers and the ANN classification and sensitivity map analysis of the data. Characteristic bands identified by the sensitivity map are linked to intracellular biochemical alterations. The diagnostic relevance of these bands is further studied by Student’s t-tests.

2.

Methods

2.1.

Sample Preparation

This study was approved by the Institutional Review Board of the Cancer Hospital of Fudan University. Written informed consent was obtained from each patient. A total of ten patients with sporadic colorectal adenocarcinomas were involved, including seven males and three females, all Asians, at ages from 33 to 79 with an average of 57.1. Tissue specimens were obtained by surgical resection from the patients. From each patient, the adenocarcinoma and a small amount of normal mucosa adjacent to the tumor site (with > 6cm separations) were removed. Histological assessments reported no tumor cell infiltration in all the normal mucosa sections. Upon removal, single-cell suspensions were prepared from small portions of the tissue specimens following primary culture protocols.3 Details of the processing procedure were described previously.1 Two suspensions were produced for each patient, one from the tumor and the other from the normal.

Histopathology analyses were conducted on the tissue specimens by two pathologists working independently. Only the cases in which consensus was reached are included in the following study. The pathology conclusions serve as the golden standard for the training and testing of ANN. The training set comprises seven moderately-differentiated cases and one poorly-differentiated case, and the test set comprises two well-differentiated cases.

2.2.

LTRS Data Collection on Single Living Cells

Data were acquired by using an LTRS system described previously.1, 17 Briefly, the excitation beam of a diode laser (wavelength 782.5nm ) was delivered into the back pupil of the objective ( 100×oil immersion) of a Nikon TE2000U differential interference contrast microscope. A laser trap was formed at the focal point. The excitation power was 11.5mW , measured near the focal point.

Single-cell suspensions were analyzed in the sample chamber mounted on the microscope stage. Single living epithelial cells were selected, captured, transferred to clear regions, and measured. The measurement of a cell comprised a 60-s signal integration and a 60-s background integration with the cell captured and released, respectively. Their difference was recorded as the raw data signal R(ν¯) , with ν¯ the Raman shift.

A total of 20 cancer cells and 20 normal cells were measured for each patient. The ten patients were divided into two groups according to the time sequence of their medical operations: data from the first eight were used to train the ANN, while data from the last two were used for testing. This generated a training set of 160 cancer and 160 normal cells and a test set of 40 cancer and 40 normal ones. Figure 1 presents the average spectra [cancer (curve a ); normal (curve b )] of the training set and their differences (curve c ).

Fig. 1

Average spectra of the cancerous (a) and normal (b) cells in the training set and their difference spectrum (c=ab) . The inset is a photo of a captured epithelial cell.

034002_1_028703jbo1.jpg

2.3.

Diagnostic Modeling by Neural Network Approach

2.3.1.

PCA reduction of spectral data

Principal component analysis is commonly used to reduce data dimensionality. The first step of PCA is a normalization of every spectrum in order to cancel out an overall factor caused by the fluctuation of excitation power. In our PCA, each spectrum was treated as an M -dimension vector, where M=1027 is the total number of pixels. Hence, a unit-vector normalization scheme was adopted,18 i.e.,

1

d(ν¯)=R(ν¯)[m=1MR2(ν¯m)]12.
The PCA performed on the 320 spectra in the training set revealed that the first 12 PCs account for 80% of the total variations. Each spectrum from Eq. 1 was then converted to the scores of corresponding PCs formulated as

2

Xi(n)=m=1M[d(n)(ν¯m)d(ν¯m)]PCi(ν¯m),i=1,,M,
where d(ν¯)(1320)n=1320d(n)(ν¯) , with index n denoting different cells of the set. In Eq. 2, the subscript m denotes the individual pixel of a spectrum, while the subscript i denotes the score. In PCA, the total number of scores equals the total number of pixels, but only the first tens of scores are useful. To facilitate network training, the scores of the training data were further rescaled to zero mean and unit variance, i.e.,

3

xi(n)=(Xi(n)Xi)σi,i=1,,M,
where Xi(1320)n=1320Xi(n) and σi[(1319)n=1320(Xi(n)Xi)2]12 .

2.3.2.

Neural network architecture

The ANN we employed is a two-layer feedforward network with backpropagation training. The network architecture is depicted in Fig. 2 and belongs to the perceptron category. The open circles labeled with h or y represent the neurons, while the solid circles labeled with x represent the inputs to the network. The two circles labeled with “1” are plotted in order to treat the biases in the same fashion as the weights. The network architecture is uniquely defined by the arrangement of the neurons and their interconnections. Each neuron can have multiple inputs but only one output. Each arrow into a neuron represents an input, and the arrow out represents the output. Each interconnection between two nodes is associated with a factor w , called weight. For example, the arrow linking xI to h1 carries a weighting factor w1I , which means that the contribution of xI to the input of neuron h1 is w1IxI [Eq. 6, shown later]. The x nodes stacked on the left serve as the inputs to the network, and their values are substituted with the scores from Eq. 3. The h neurons stacked in the middle form the hidden layer, and the y neuron on the right denotes the output layer. The numerical outputs of the h neurons and y neuron are given by Eqs. 6, 7 later, respectively. Here, the output layer is set to one single node, since a cell is classified as either normal or cancerous. So our network is essentially a binary neural classifier (with modifications described later). The number of nodes in the input column I and the hidden layer H were determined in the training session by trial and test. Best results were achieved at I=5 and H=8 .

Fig. 2

Structure of the binary feedforward neural network.

034002_1_028703jbo2.jpg

The extreme flexibility of multilayer perceptron networks to approximate arbitrary complex posterior probability functions calls for special procedures to avoid overfitting and detect outliers. Outliers are data points that have been assigned to a wrong category of patterns in the original data set. They are defects of the experimental data. For example, if a cell in the training set is truly normal, but instead determined as cancerous by pathology, then it is an outlier, or vice versa. The existence of outliers is an outcome of the imperfectness of data collection. In our experimental data, the pattern of a cell was set to the pathological classification of the tissue from which the cell was harvested: if the tissue was normal, the cell was normal, or vice versa. However, there was a small possibility that a cell inside a malignant tumor was actually normal, but erroneously labeled as cancerous based on the histopathology grading on the tumor tissue, resulting in an outlier.

Outliers in the training set can significantly shift the decision plane of the network if uncorrected and consequently affect the accuracy of network predictions. An outlier probability ϵ=[0,1] was introduced in Ref. 16 to correct the classification flipping due to the existence of outliers. In our binary case, this means that

4

{P(1x)P0(1x)(1ϵ)+ϵP0(0x),P(0x)P0(0x)(1ϵ)+ϵP0(1x),P0(0x)=1P0(1x),
where P0(1x) and P0(0x) are the posterior probability for cancer (class 1) and normal (class 0) at the absence of outliers, respectively and P(1x) and P(0x) are the actual posterior probability for cancer and normal, respectively. Here the term “posterior” means post-experiment, and the posterior probability is the conditional probability of the cell’s classification under the requirement that it has the characteristics to produce Raman spectrum x . Both P(1x) and P(0x) contain two terms, indicating that the cell data x carry 1ϵ probability of not being an outlier and ϵ probability of being an outlier. After straightforward calculations, Eq. 4 can be simplified to

5

{P(1x)=P0(1x)(12β)+β,P(0x)=1P(1x),
where βϵ . The connection of Eq. 5 to the ANN output (Fig. 2) is through P0(1x) .

In ANN, the response (the output) of a neuron to the stimuli (the inputs) is calculated through its activation function. In Fig. 2, the activation functions are hyperbolic tangent for the hidden layer, i.e.,

6

hj(x)=tanh(i=1Iwjixi+wj0),
and linear for the output layer, i.e.,

7

y(x)=j=1Hw̃jhj(x)+w̃0,
where wji and wj0 are the input-to-hidden weight and bias, and w̃j and w̃0 are the hidden-to-output weight and bias. The ANN estimate to the posterior probability without outlier corrections P0(1x) is predicted as

8

P̂0(1x)=11+exp[y(x)].
Consequently, the ANN estimates to Eq. 5 become

9

{P̂(1x)=12β1+exp[y(x)]+β,P̂(0x)=1β12β1+exp[y(x)].

2.3.3.

Input-output and training of the network

The training set comprises the cell data of eight patients, D={(x(n),t(n))n=1,2,,320} , with x(n) the rescaled scores in Eq. 3, and t(n) the target value defined as

10

t(n)={1,ifcelliscancerous;0,ifcellisnormal.
The procedures of Ref. 16 were executed to train the network, but with some modifications. The network optimization was done by minimizing a cost function redefined as

11

S(w)=n=1320e(n)(w,β)+αEW(w),
where

12

e(n)(w,β){t(n)P̂(1x(n))+(1t(n))[1P̂(1x(n))]},
and

13

EW(w)12i=1Wwi2.
The first term on the right side of Eq. 11 is the cross-entropy error function, while the second one is a weight-decay regularization to restrict network overfitting on the input noise. The weight vector w is a simplified notation for all the network weights and biases, i.e., w{wji,w̃ki=0,,I;j=1,,H;k=0,,H} and its dimension W=(I+2)H+1 . It is worth pointing out that Eq. 11 is different from the special case of c=2 in Ref. 16 for several reasons, where c is the number of output classes. The first is that our case employs only one instead of two sets of hidden-to-output weights, so the network is simpler and the dimension of w smaller; the second, P̂(1x(n)) in Eq. 12 is calculated using the logistic function, while in Ref. 16, the softmax function is used and the two calculated posterior probabilities are not independent, inducing problems in the evaluation of the inverse Hessian and causing a singularity of S(w)β at β=0 . Indeed, our numerical computations confirmed that Eqs. 11 to 13 produce better training results than the c=2 version of Ref. 16.

The unknown parameters are the network weights-biases w and the hyperparameters α [Eq. 11] and β [Eq. 9]. The training is carried out by iterations on two nested loops. In the inner loop, the weights are optimized using a Broyden-Fletcher-Goldfarb-Shanno (BFGS) quasi-Newton algorithm19 to minimize the cost function S(w) at fixed (α,β) ; in the outer loop, the hyperparameters (α,β) are adapted by maximizing the evidence20, 21 p(Dα,β) at fixed w . (For original descriptions, see Refs. 16 and 22). We rephrase the key points on adapting α and β here for the purpose of completeness. The Gauss-Newton approximation to the Hessian matrix is first computed as

14

A(w)=n=1320e(n)(w,β)we(n)(w,β)wT+αIW×W.
The W×W symmetric matrix A(w) is positive definite if α> 0 , so accordingly its determinant is always positive; as a result, the logarithmic in Eq. 16 is defined. The identity matrix in Eq. 14 comes from the second order derivatives of EW(w) in Eq. 11. The hyperparameter α is updated as

15

αnew=[WαTraceA1(w)]i=1Wwi2,
and the scaled outlier probability β is updated by minimizing the function

16

C(β)S(w)+12ln[detA(w)].

2.3.4.

Sensitivity maps

Despite the lack of physical interpretation of neural network parameters, sensitivity maps of ANN were suggested as an essential tool to establish explicit correlation between the inputs and the outputs.16 At the network inputs end (i.e., the scores), one could further map this correlation back to the spectral space through the relationship between the scores and PCs. In this way, the sensitivity map is capable of characterizing the contributions of individual vibrational bands to the classification of cancer versus normal.

Physically, this means a biochemical change correlated to a tumor progression is expressed as intensity variations at corresponding Raman bands, which in turn induce a perturbation in the network output. The higher the ratio between the perturbation to the input variation, the more relevant the corresponding chemical substance is to the course of disease transformation.

Mathematically, the sensitivity map is defined as the derivative of the estimated posterior probability P̂(1x) with respect to the Raman intensity. The absolute-value-average sensitivity is one type of such maps,16

17

s(ν¯m)=1320n=1320P̂(1x(n))d(ν¯m),m=1,,M,
where

18

P̂(1x)d(ν¯m)=(12β)exp[y(x)]{1+exp[y(x)]}2j=1Hw̃j[1hj2(x)]i=1IwjiσiPCi(ν¯m).
Finally, the sensitivity map is formulated as the unit-vector normalization of Eq. 17, i.e.,

19

s̃(ν¯)=s(ν¯)[m=1Ms2(ν¯m)]12.

3.

Results and Discussion

3.1.

Classification

The network was trained on the training set described in Sec. 2.2, comprising Raman spectra of 160 cancerous and 160 normal cells from eight patients. The network weights and biases were initialized with random numbers of Gaussian distribution at zero mean and σ2=1I . Best training results were accomplished at I=5 and H=8 by trial and test on various combinations of (I,H) ; the network performance started to deteriorate when I6 .

We therefore restrict the following discussions to the case of I=5 and H=8 . The cost function S(w) had 47% chance of hitting the global minimum after runs on different initializations. Results at local minima were also close but were discarded anyway. There were no observable differences among the network predictions once the global minimum was reached, although the settled network weights varied from run to run. The final values of the hyperparameters were α=0.7126 and βϵ=0.0374 . The average outlier probability,16 defined as the average percentage of the second term inside the two-term summation on the right side of Eq. 4, was found to be 0.0252.

After the values for the network weights and hyperparameters were settled, the estimated posterior P̂(1x) was calculated for each cell and compared with a threshold value pth to predict its class; it was cancerous if P̂(1x)> pth , or normal otherwise. A reasonable value pth=0.5 was chosen. In the final classifications of the training set, 138 out of 160 cancerous cells and 138 out of 160 normal cells were correctly identified, reaching a sensitivity of 86.3% and specificity of 86.3%. These were considerably better than the results given by logistic regressions1 (Table 1 ).

Table 1

Predictions of this model versus the logistic regression model.1 Major improvements are observed in the training set

Pathology analysisANNLogistic regressions
NormalCancerNormalCancer
Training setNormal1382213030
Cancer2213836124
Sensitivity86.3%77.5%
Specificity86.3%81.3%
Test set, patient 1Normal173173
Cancer317317
Sensitivity85%85%
Specificity85%85%
Test set, patient 2Normal200200
Cancer317416
Sensitivity85%80%
Specificity100%100%

A receiver operating characteristic (ROC) curve is commonly used by medical doctors to assess the performance of a diagnostic model. It displays the dynamic trend of the model’s predictions. In fact, it is possible to sacrifice specificity for higher sensitivity by reducing pth , or vice versa. For example, if we set pth=0 , all the cells will be identified as cancerous, corresponding to 100% sensitivity but 0% specificity. By sweeping pth from 0 to 1, one can map out the ROC curve. In a perfect model, the curve is a step function along the upper-left corner (Fig. 3 , the OAB line). On the other hand, a random-guess model possesses no predictive power and thus produces an ROC curve along the diagonal. The closer an ROC curve to the upper-left corner, the better the model’s performance. In Fig. 3, the ROC of the current model shows much improvement over that of the logistic regression model.1

Fig. 3

ROC curves for this model (solid line), logistic regression model1 (dashed line), and a random-guess model (dotted line). The OAB line along the upper-left corner represents a perfect model, while the random-guess model flips a coin to assign a cell’s class.

034002_1_028703jbo3.jpg

In a double-blinded test of the ANN, we fed the data from two new patients to evaluate its capability of predicting new results, the so-called generalization. The test set, described earlier in Sec. 2.2, comprised the Raman spectra of 40 cancerous and 40 normal cells. The network predicted correctly 17 of 20 cancer cells and 17 of 20 normal cells for the first patient and 17 of 20 cancer cells and 20 of 20 normal cells for the second patient. The overall sensitivity and specificity were 85% and 92.5% respectively, compared with the corresponding 82.5% and 92.5% of logistic regression. The limited improvement in the test may be due to the samples. In fact, the tissue specimens of the test set were graded as well-differentiated by pathology, and a higher degree of differentiation in pathology points to a smaller difference between the normal and abnormal. In pathological terms, well-differentiated tissues show a regular order similar to that of normal tissues, where there are clear structures. Also cells from well-differentiated tissues exhibit morphology similar to those in normal tissues, such as the cell shapes, nuclear sizes, and uniformity of cells inside the tissues. On the other hand, in poorly differentiated tissues, the cells are aligned in a chaotic way, indicating a higher degree of invasiveness. Here, our test results may provide evidence that morphological similarities and biochemical similarities between the well-differentiated cells and the normal cells are correlated.

Nevertheless, the generalization of our ANN is excellent. A common problem in ANN algorithms is that a well-trained network predicts very poorly on new data. This is absent in our model.

3.2.

Sensitivity Analysis

The sensitivity map is presented in Fig. 4 . The split-half resampling technique was applied to investigate the reproducibility of sensitivity maps.16 The Z scores of 100 split-half resamplings (corresponding to 100 pairs of independent networks trained separately on split-half data) were found to be above 2.3657, corresponding to a minimum 99.1% confidence interval. Therefore, Fig. 4 is reproducible at above 99.1% confidence level.

Fig. 4

Sensitivity map of the network and the training data. Major features are identified here and listed in Table 2. Do not confuse the sensitivity here with that in Fig. 3.

034002_1_028703jbo4.jpg

Table 2

Tentative assignments of the Raman bands identified in Fig. 4, mean intensity changes (increase + /decrease − ) of cancer relative to normal, and p-values of Student’s t-tests on band intensities. The spectral resolution of our system is 8cm−1 .

Bands (cm−1) Tentative assignmentMean changep-Value
623.985 CC twist Phe 0.000
645.438 CC twist Tyr 0.000
759.803Trp 0.324
841.568Glucose 0.001
940.988 CC BK str. α helix + 0.615
1004Sym. ring br. Phe + 0.333
1034.556 CH in-plane Phe 0.000
1096.433 PO2 str. DNA/RNA backbone + 0.395
1132.941 CC stretch protein 0.000
1156.272 CC & CN str. proteins, carotenoids 0.022
1209.015Phe, Trp + 0.145
1265.006 CH in-plane lipid + 0.000
1304.081 CH2 twist lipid/adenine, cytosine + 0.000
1341.605adenine, guanine + 0.188
1380.083thymine 0.000
1437.853 CH2 deformation lipid + 0.000
1520.245 CC carotenoid + 0.066
1655.304Amide I α helix, CC str. lipid + 0.000
1672.795Cholesterol + 0.001
1746.711 CO str. lipid 0.330

Prominent Raman bands are easily identifiable on Fig. 4. Tentative assignments23 are listed in Table 2 . Intensities at these peaks were read out from the data and subjected to Student’s t-test. We observed mean intensity increases from normal to cancer at 940.988, 1004, 1096.433, 1209.015, 1265.006, 1304.081, 1341.605, 1437.853, 1520.245, 1655.304, and 1672.795cm1 and decreases at the remaining peaks. The increases occur mostly at the proteins and nucleic acids bands, indicating that cancer cells have higher DNA/RNA and protein activities. It is a surprise that glucose (decreasing, p- value=0.001 ) and cholesterol (increasing, p- value=0.001 ) appear on the list. However, signal alterations at 759.803, 940.988, 1004, 1096.433, 1209.015, 1341.605, and 1746.711cm1 cannot be ruled out as statistical fluctuations (p- value> 0.05 ). Furthermore, the cancer cells systematically (p- value<0.05 ) have lower signals than the normal ones in the spectral region below 720cm1 (see Fig. 1, for example), which may require further investigation. This phenomenon is also reflected in Fig. 4 as increased sensitivity levels at the same region, and confirmed by the colon tissue spectra in Fig. 5B of Ref. 23. Finally, even though mean intensity changes with p-values less than 0.05 are considered as statistically significant by medical doctors, a p-value greater than 0.05 does not necessarily reject the correlation of a band to the disease transformation.

3.3.

Conclusions and Discussions

In this paper, we compared the technical characteristics of DNA FCM and LTRS, including working principles, sample preparations, hardware aspects, data collection schemes, operational procedures, and sample amount requirements, for single-cell investigation of cancers. LTRS can share the single-cell processing procedures already existing in hospitals for FCM, but removes the requirement of fluorescent dye staining. Unlike DNA FCM, which replies on the accumulative histogram of tens of thousands of cells, LTRS can decide on each individual cell, which makes LTRS very advantageous for diagnosing small tumors. In the meantime, by detecting the biochemical activities of cancers at the single-cell level, LTRS can provide deeper insight into the dynamics of carcinogenesis than tissue Raman spectroscopy, because cell functions are more fundamental than tissue functions.

Potentially, single-cell LTRS can achieve higher accuracy than tissue Raman spectroscopy, because the course of carcinogenesis starts within the epithelial cells, and the epithelial cells are the primary source of diagnostic information for early cancer detection. In literature, some tissue components like collagen are reported to show correlations to cancer developments. However, because genetical alterations are the foundation of structural changes, variations in these components are by-products of carcinogenic activities in the epithelial cells and thus serve only as the secondary source of diagnostic information. Furthermore, many other tissue components do not possess diagnostic information at all. Tissue Raman spectroscopy suffers from two sources of interference: one is the Raman signal from the irrelevant tissue components; the other is the strong fluorescence and its associated noise from the tissue bulk. Even at 785-nm excitation, tissue fluorescence is still considerably higher than tissue Raman signal, and mathematical preprocessing such as fifth-order polynomial fitting is typically required to subtract the slow-varying fluorescence background from the data. Unfortunately, the noises associated with the fluorescence are left in the processed data. Therefore, in Raman studies, it is desirable to avoid as many fluorescence sources as possible.

In a clinical study of colorectal adenocarcinoma, we measured the Raman spectra of 400 living epithelial cells from ten patients using LTRS. An ANN algorithm incorporating outlier probability corrections was developed to overcome the limitations existing in a previous study employing logistic regression.1 Remarkable improvements were achieved in the sensitivity and specificity of the predictions. Unavailable to logistic regression, the sensitivity map of ANN was evaluated to objectively quantify the contribution of each Raman frequency in the course of carcinogenesis. Important Raman bands were identified from the map, and Student’s t-tests were performed on their intensities. Unlike the subjective visual inspections commonly used in spectral analysis, the sensitivity map technique provides an objective, quantitative, and automatic way to discover important Raman peaks.

In our experiments, only living epithelial cells were measured. The following discussions are restricted to living cells only. We cannot rule out the possibility of potential damage to some membrane proteins on the cell surface during the process of cell isolation from the tissue, but the primary culture procedure we employed is well-established and should reduce such damage to a minimum. In the meantime, the intracellular substances should be intact. Over the past decades, there have been intense studies on cancer cell culture techniques to prevent damage to the cells. Today, these techniques are already mature and widely applied in various fields. One example is the cell lines on the market, which came from the primary cultures originally harvested from patients. The demand for membrane protein protection might originally come from conventional FCM, which requires the binding of fluorescence markers to specific targets on the cell surface. After years of development, different protocols for different tissue types are already available in standard textbooks. LTRS on single cells can directly enjoy the fruits of other fields to keep possible cell damage under control. Furthermore, the major signal of LTRS originates from the volume inside the cell, so we do not expect a detectable alteration of the cell’s Raman spectrum due to the isolation processing.

The average Raman spectra of single epithelial cells [Figs. 1a and 1b] share many common features with those of tissues (Fig. 5B in Ref. 23). They exhibit a general resemblance in the overall spectral shapes. However, differences in spectral details are obvious due to the extra components in tissue other than the epithelial cells. The difference spectra between the abnormal and the normal [Fig. 1c in this paper and Fig. 6B in Ref. 23] also show similarity in general trends, but differ in details.

To establish our diagnostic model, we employed the tissue histopathology assessments as the golden standard and involved only those cases with consensus reached between two independent pathologists. A cell’s classification was set to the pathology classification of the tissue where it originated. There was a small possibility that a cell from a malignant tumor might be actually normal, resulting in an outlier. Unfortunately, it is impractical to obtain a cell’s classification from a direct single-cell histochemical analysis. To the best of our knowledge, immuno-histochemical analysis on cells/tissues is not for diagnostic purpose and cannot determine the normal/malignant nature of samples. It is rather used to supply assistive information for histopathology assessments. Besides, immuno-histochemical information alone is not reliable for diagnosis because the targeted proteins are also expressed in normal cells/tissues. Our way to assign a cell’s classification should provide the best accuracy for model training. In both the training set and the test set, around 14% of all cancer cells were misidentified by our model. But in our opinion, the outliers are only a minor source for this error, since it is highly unlikely that our cancer tissue samples contained such a high percentage of normal cells. As a matter of fact, our diagnostic model was built purely on experimental data with no prior knowledge and could be affected by the noises in the signals as well as the statistical nature of the individual variations of the cells.

Unlike the training stage of our model, the application stage does not require prior histopathology assessment on a sample. LTRS diagnosis can be run independent of or without histopathology. In a clinical application, cells from suspected tumor sites can be retrieved through fine needle aspiration or exfoliative cytology methods like brushing and washing and then examined by LTRS to determine the nature of the cells. Frequently, doctors do not know the definitive characteristics of a suspected site without resorting to sophisticated diagnostic techniques. This is where histopathology, LTRS, DNA FCM, etc. come into play.

In many situations, doctors need to know only whether a suspected cancer transform has occurred in the first place, before further investigation can be triggered to find its exact location. Exfoliative cells in pleural fluid, ascites, and urine can be collected and analyzed by LTRS to help doctors make decisions.

The extension of our work to other epithelial cancers is straightforward. The single-cell diagnosis technique presented here can be directly applied to clinical tests. Thanks to FCM, preparation of single-cell suspensions is already routine in hospitals today. Further, exfoliative cytology can also provide samples to our technique. The low requirement for sample amounts makes our technique very attractive to cancer screening.

The current work deals only with the question “Normal or cancer?” Future studies will focus on the subclassification of different tumor stages.

Acknowledgments

We thank Dr. Sigurdsson for a private communication on artificial neural networks. This work was partially supported by Shanghai Optical-Tech Special Project Grant No. 036105016.

References

1. 

Chen K., Qin Y. J., Zheng F., Sun M. H., Shi D. R., “Diagnosis of colorectal cancer using Raman spectroscopy of laser-trapped single living epithelial cells,” Opt. Lett., 31 2015 –2017 (2006). 10.1364/OL.31.002015 0146-9592 Google Scholar

2. 

Chan J. W., Taylor D. S., Zwerdling T., Lane S. M., Ihara K., Huser T., “Micro-Raman spectroscopy detects individual neoplastic and normal hematopoietic cells,” Biophys. J., 90 648 –656 (2006). 10.1529/biophysj.105.066761 0006-3495 Google Scholar

3. 

MacLeod K. G., Langdon S. P., “Essential techniques of cancer cell culture,” Methods Mol. Med., 88 17 –29 (2004). 1543-1894 Google Scholar

4. 

Costa J., Cordon-Cardo C., Devita V. T., Hellman S., Rosenberg S. A., “Cancer diagnosis: molecular pathology,” 641 –657 (2001) Google Scholar

5. 

“Recommended guideline standards for DNA investigations,” Purdue Cytometry CD-ROM, Purdue University, West Lafayette, IN (1997). Google Scholar

6. 

Haka A. S., Shafer-Peltier K. E., Fitzmaurice M., Crowe J., Dasari R. R., Feld M. S., “Diagnosing breast cancer by using Raman spectroscopy,” Proc. Natl. Acad. Sci. U.S.A., 102 12371 –12376 (2005). 10.1073/pnas.0501390102 0027-8424 Google Scholar

7. 

Haka A. S., Shafer-Peltier K. E., Fitzmaurice M., Crowe J., Dasari R. R., Feld M. S., “Identifying microcalcifications in benign and malignant breast lesions by probing differences in their chemical composition using Raman spectroscopy,” Cancer Res., 62 5375 –5380 (2002). 0008-5472 Google Scholar

8. 

Crow P., Molckovsky A., Stone N., Uff J., Wilson B., Wongkeesong L. M., “Assessment of fiberoptic near-infrared Raman spectroscopy for diagnosis of bladder and prostate cancer,” Urology, 65 1126 –1130 (2005). 10.1016/j.urology.2004.12.058 0090-4295 Google Scholar

9. 

Gniadecka M., Philipsen P. A., Sigurdsson S., Wessel S., Nielsen O. F., Christensen D. H., Hercogova J., Rossen K., Thomsen H. K., Gniadecki R., Hansen L. K., Wulf H. C., “Melanoma diagnosis by Raman spectroscopy and neural networks: structure alterations in proteins and lipids in intact cancer tissue,” J. Invest. Dermatol., 122 443 –449 (2004). 10.1046/j.0022-202X.2004.22208.x 0022-202X Google Scholar

10. 

Motz J. T., Fitzmaurice M., Miller A., Gandhi S. J., Haka A. S., Galindo L. H., Dasari R. R., Kramer J. R., Feld M. S., “In vivo Raman spectral pathology of human atherosclerosis and vulnerable plaque,” J. Biomed. Opt., 11 021003 (2006). 10.1117/1.2190967 1083-3668 Google Scholar

11. 

Hanlon E. B., Manoharan R., Koo T. W., Shafer K. E., Motz J. T., Fitzmaurice M., Kramer J. R., Itzkan I., Dasari R. R., Feld M. S., “Prospects for in vivo Raman spectroscopy,” Phys. Med. Biol., 45 R1 –R59 (2000). 10.1088/0031-9155/45/2/201 0031-9155 Google Scholar

12. 

Smith Z. J., Berger A. J., “Surface-sensitive polarized Raman spectroscopy of biological tissue,” Opt. Lett., 30 1363 –1365 (2005). 10.1364/OL.30.003365 0146-9592 Google Scholar

13. 

Singh G. P., Volpe G., Creely C. M., Grotsch H., Geli I. M., Petrov D., “The lag phase and G(1) phase of a single yeast cell monitored by Raman microspectroscopy,” J. Raman Spectrosc., 37 858 –864 (2006). 10.1002/jrs.1520 0377-0486 Google Scholar

14. 

Xie C., Mace J., Dinno M. A., Li Y. Q., Tang W., Newton R. J., Gemperline P. J., “Identification of single bacterial cells in aqueous solution using conflocal laser tweezers Raman spectroscopy,” Anal. Chem., 77 4390 –4397 (2005). 10.1021/ac0504971 0003-2700 Google Scholar

15. 

Short K. W., Carpenter S., Freyer J. P., Mourant J. R., “Raman spectroscopy detects biochemical changes due to proliferation in mammalian cell cultures,” Biophys. J., 88 4274 –4288 (2005). 10.1529/biophysj.103.038604 0006-3495 Google Scholar

16. 

Sigurdsson S., Philipsen P. A., Hansen L. K., Larsen J., Gniadecka M., Wulf H. C., “Detection of skin cancer by classification of Raman spectra,” IEEE Trans. Biomed. Eng., 51 1784 –1793 (2004). 10.1109/TBME.2004.831538 0018-9294 Google Scholar

17. 

Deng J. L., Wei Q., Zhang M. H., Wang Y. Z., Li Y. Q., “Study of the effect of alcohol on single human red blood cells using near-infrared laser tweezers Raman spectroscopy,” J. Raman Spectrosc., 36 257 –261 (2005). 10.1002/jrs.1301 0377-0486 Google Scholar

18. 

Deinum G., Rodriguez D., Romer T. J., Fitzmaurice M., Kramer J. R., Feld M. S., “Histological classification of Raman spectra of human coronary artery atherosclerosis using principal component analysis,” Appl. Spectrosc., 53 938 –942 (1999). 10.1366/0003702991947829 0003-7028 Google Scholar

19. 

Nielsen H., “UCMINF—an algorithm for unconstrained nonlinear optimization,” (2000). Google Scholar

20. 

MacKay D. J. C., “A practical Bayesian framework for backpropagation networks,” Neural Comput., 4 448 –472 (1992). 0899-7667 Google Scholar

21. 

MacKay D. J. C., “The evidence framework applied to classification networks,” Neural Comput., 4 720 –736 (1992). 0899-7667 Google Scholar

22. 

MacKay D. J. C., “Comparison of approximate methods for handling hyperparameters,” Neural Comput., 11 1035 –1068 (1999). 10.1162/089976699300016331 0899-7667 Google Scholar

23. 

Stone N., Kendall C., Smith J., Crow P., Barr H., “Raman spectroscopy for identification of epithelial cancers,” Faraday Discuss., 126 141 –157 (2004). 10.1039/b304992b 0301-7249 Google Scholar
© (2007) Society of Photo-Optical Instrumentation Engineers (SPIE)
Feng Zheng, Feng Zheng, Kun Chen, Kun Chen, Yejun Qin, Yejun Qin, "Sensitivity map of laser tweezers Raman spectroscopy for single-cell analysis of colorectal cancer," Journal of Biomedical Optics 12(3), 034002 (1 May 2007). https://doi.org/10.1117/1.2748060 . Submission:
JOURNAL ARTICLE
9 PAGES


SHARE
Back to Top