Abstract

Preeclampsia is associated with increased perinatal morbidity and mortality. There have been numerous efforts to determine preeclampsia biomarkers by means of biophysical, biochemical, and spectroscopic methods. In this study, the preeclampsia and control groups were compared via band component analysis and multivariate analysis using Raman spectroscopy as an alternative technique. The Raman spectra of serum samples were taken from nine preeclamptic, ten healthy pregnant women. The Band component analysis and principal component analysis-linear discriminant analysis were applied to all spectra after a sensitive preprocess step. Using linear discriminant analysis, it was found that Raman spectroscopy has a sensitivity of 78% and a specificity of 90% for the diagnosis of preeclampsia. Via the band component analysis, a significant difference in the spectra of preeclamptic patients was observed when compared to the control group. 19 Raman bands exhibited significant differences in intensity, while 11 of them decreased and eight of them increased. This difference seen in vibrational bands may be used in further studies to clarify the pathophysiology of preeclampsia.

1. Introduction

Preeclampsia is a pregnancy-specific disorder that affects up to 15% of all pregnancies. The specific conditions of the disease can basically be defined as hypertension (≥140/90 mmHg) and proteinuria. These conditions occur after the 20–22th week of pregnancy [1]. Preeclampsia brings risks to both mother and baby and causes maternal morbidity and fetal mortality [2].

One of the causes of preeclampsia is thought to be a lack of antioxidants [3]. Free radicals or reactive oxygen species can interact with lipids and proteins which may cause lipid peroxidation [4] and protein modification [5]. It is suggested that the alteration of lipid and protein levels is correlated with oxidative stress and vascular dysfunction in preeclampsia [4].

The aetiology and pathology of preeclampsia are still of great interest to researchers. There have been numerous efforts to determine preeclampsia biomarkers employing biophysical, biochemical, and spectroscopic methods. D’Aniello et al. [6] used a high-performance liquid chromatography (HPLC) fluorometric method to predict the preeclampsia serum biomarkers and reported reduced levels of free amino acids, especially arginine levels which dramatically decreased in preeclampsia patients. On the other hand, Glew et al. [7] found no significant difference between the amino acid levels of preeclamptic and control groups’ plasma samples. On the contrary, López-Quesada [8] reported an increase in most amino acid levels. Turner et al. [9] obtained higher tyrosine, histidine, and phenylalanine concentrations for preeclamptic plasma samples using nuclear magnetic resonance spectroscopy. In another study [10], they reported lower lipid concentrations for preeclamptic plasma samples. On the contrary, there are studies reporting elevated lipid levels in the serum of preeclamptic patients. Robinson et al. [11] showed altered lipid metabolism characterized by elevated circulating triglycerides and nonesterified free fatty acids in preeclampsia. Hubel et al. [12] found a decreased particle size of low-density lipoprotein and increased triglyceride, free fatty acids, total cholesterol, and LDL-cholesterol concentrations in preeclampsia.

A very recent study which used Fourier transform infrared spectroscopy [13] reported decreased intensities for protein bands from the spectra of lyophilized serum samples. Furthermore, they found a relationship between preeclampsia and lipid parameters, such as increased triglycerides and decreased HDL-cholesterol during preeclampsia.

A study of the available literature reveals conflicting views about determining the serum levels of these biomarkers in body fluids from healthy and preeclamptic women. Therefore, there is a need for complementary studies to explore the pathophysiology of preeclampsia.

Raman spectroscopy is a potent tool that gives information about molecular structures. Raman spectroscopy is based on vibrational transitions between one vibrational state and a virtual state as a result of the inelastic collision of light and molecules. The frequencies of transition lines and their intensities are characteristic for each molecular vibration which may provide information about the conformation of macromolecules such as proteins [14, 15], nucleic acids [16, 17], and lipids [18, 19]. Thus Raman spectroscopy can be a suitable tool to determine biomarkers which may play a role in the pathophysiological mechanism of diseases. Several diseases were investigated on human body fluids such as whole blood [20, 21] and blood serum [22, 23] by using Raman spectroscopy.

Multivariate methods with Raman spectroscopy in the determination of the concentration of some constituents in serum or whole blood were performed. The partial least squares method [2428] was used to predict the concentration of blood analytes including glucose, urea, cholesterol, triglycerides, total protein, albumin, and hemoglobin. Furthermore, a number of studies have elucidated the diagnostic potential of Raman spectroscopy combined with multivariate analysis techniques for several diseases [2934].

To our knowledge, this is the first Raman spectroscopic study on preeclampsia to exist in the literature. In the study, we measured the Raman spectra of serum samples taken from preeclamptic and healthy women. We analyzed the measured spectra with band component analysis and applied principal component analysis-linear discriminant Analysis (PCA-LDA) after a sensitive preprocess step to find the differences between preeclamptic and healthy women.

2. Materials and Methods

2.1. Sample Collection and Preparation

Informed consent for the study was obtained from all the women concerned, and approval was given by the Ethics Committee of the Faculty of Medicine, Istanbul University. The sample population consisted of ten normotensive women who had no underlying medical conditions and nine otherwise healthy preeclamptic women who were diagnosed at Cerrahpasa Medical Faculty, Istanbul University. The normal pregnancies were followed until their completion, and no preeclampsia occurred.

Blood samples taken in 10 mL gel tubes were centrifuged during a 10-minute period with 3000 rpm to obtain serum samples. Each serum sample was prepared in two aliquots which were measured on two different days to validate the stability of the experimental conditions. The samples were kept at −20°C until the observation day and were measured maximum in seven days. Frozen serum samples in eppendorf tubes were defrosted by hand warmth before the Raman measurements were taken. About 1 mL of serum samples were prepared in glass cuvettes which were cleaned properly after every measurement.

2.2. Experiment

The experimental setup consists of a 785 nm, 100 mW single mode diode laser (CrystaLaser) and a spectrograph ( 𝑓 = 3 0 3  mm, 𝑓 #4.3, Andor) equipped with a 600 lines/mm grating and with TE cooled CCD (−90°C, Andor iDus DU420A-OE). A Faraday isolator (EOTech, BB8-00011) was placed just in front of the diode laser to prevent back reflections which cause instabilities on single longitudinal mode operation and the power of the laser. A laser line filter was used to obtain a clear laser emission (Semrock, LL01-780-12.5). The laser beam was directed and then focused on the sample by suitable mirrors and lenses. Back scattered light from the sample was collected at 180° geometry by the same lens. The scattered light from the sample was filtered two times successively with Raman edge filters which transmit only the Raman scattering (inelastic scattering) beam while they reflect the Rayleigh (elastic scattering) beam. The Raman scattering beam was inserted into the spectrograph through 100 μm entrance slit by an achromat lens with a focal length of 50 mm.

The Raman spectra of 38 serum samples (nine preeclamptic and ten healthy; two from each patient) were taken with a 30 s exposure time, 20 successive scans over a 10-minute period. The laser power on the sample was about 50 mW. After cosmic-ray removal, the 20 scans were averaged for each sample. The wavenumber calibration of the spectrograph was performed with reference bands of toluene [35]. A toluene Raman spectrum was taken for one second after each Raman spectrum of serum sample to apply a developed iterative calibration technique which minimizes the effects of possible spectral shifts during the measurement day(s). These shifts may arise from some reasons such as temperature, pressure, humidity which may affect the diode laser, spectrograph, and all optical components. The mean Raman spectrum of all serum samples for control and preeclamptic groups and the Raman spectrum of water are shown in Figure 1(a). The Raman spectrum of water was taken at the end of each measurement day. The wide Raman band from water at around 1640 cm−1 and other wide bands from glass at around 1400 cm−1 and 1100 cm−1 are visible. However, the background correction by subtracting the water Raman spectrum from each serum spectrum eliminates these interferences.

Although the fluorescence could be significantly reduced with the excitation at 785 nm, the residual fluorescence which originates from macromolecules in serum [36], from glass and also stray light contributions, should be eliminated. Therefore, baseline corrections were applied by fitting a cubic spline curve for each spectrum and subtracting the curves from the spectra. The wavenumbers corresponding to the data points that constitute the curve for the baseline were chosen to be identical for each curve in order to assure objectivity among the processed spectra. Then all spectra were vector normalized. The baseline corrected, vector normalized mean Raman spectrum of control samples is given in Figure 1(b).

Second derivative spectra were used to determine the number and initial positions of the components under broad bands. By using this information all spectra were fitted with a Voigt profile using commercial software, Grams/AI 8.0, to determine band intensities and positions. Figure 1(c) is an example of a fit of the part of the serum Raman spectrum given in Figure 1(b).

Additionally, the Raman spectra of 15 amino acids, β-carotene, glucose, glutathione, adenine, urea, and serum albumin in water (0.1 M) were measured individually to determine the band positions and intensities. As an example, the Raman spectra of 0.1 M glycine and L-leucine in water are given in Figure 2 superposed with the mean Raman spectrum of serum samples of the control group to show the contribution of glycine and leucine bands. All of the chemicals used in this work were purchased from Merck and Sigma.

2.3. Statistical Analysis

All of the measured Raman spectra established a data set for PCA and LDA which are two of the most widely used data classification methods. In the PCA process, all spectra constitute a matrix whose rows are observations and columns are wavenumbers. The PCA software first applies a mean centering on this matrix and then determines the eigenvalues and eigenvectors of the covariance matrix of the data and thus maximizes the variance between variables. Therefore, defining eigenvector 𝑉 , transpose of eigenvector 𝑉 , scores 𝑇 and mean spectrum 𝐴 ; PCA score vectors are expressed as below: 𝑇 = 𝑉 𝐴 . ( 2 . 1 ) Having defined original spectrum 𝑋 , the expression below can be obtained by means of scores and loadings: 𝑋 = 𝐴 + 𝑇 𝑉 + 𝜀 , ( 2 . 2 ) where 𝜀 means residual matrix [48, 49]. In spectral analysis, the scores lead to a classification between groups. Since eigenvalues are sorted in descending order in the PCA process, when the first most important scores are plotted versus each other, spectral differences will occur as discriminations in score plot.

The LDA analysis was applied using Fisher’s criterion which is actually a Bayesian decision rule. Defining 𝑤 = Σ 1 𝜇 1 𝜇 2 , 𝑤 0 1 = 2 𝜇 1 + 𝜇 2 Σ 1 𝜇 1 𝜇 2 𝑇 𝜋 l o g 1 𝜋 2 , ( 2 . 3 ) Fisher’s criterion is given by 𝑤 𝑇 𝑥 + 𝑤 0 > 0 , ( 2 . 4 ) where 𝜇 1 and 𝜇 2 are group mean vectors, 𝜋 1 and 𝜋 2 are prior probabilities for each group, 𝐱 is a row vector which consists of two most important PC scores, and 𝚺 is pooled covariance matrix [48, 50].

PCA and LDA processes were applied to all preprocessed spectra. All PCA analyses were made on Grams PLS Plus/IQ software with multiplicative scatter correction preprocessing and leave-one-out cross-validation.

The LDA process was performed with software programmed on MATLAB 7.4.0.

3. Results

The preeclamptic and normotensive groups were well matched. There was no significant difference for age, parity, and gestational age at blood sampling between the groups ( 𝑃 > 0 . 0 5 ). The Raman spectra of 38 serum samples (18 preeclamptic and 20 healthy) were obtained. The vector normalized mean spectra of serum samples of the preeclamptic and control groups are shown in Figure 3. The figure demonstrates that there are differences in some band intensities between the two spectra indicated with arrows and corresponding wavenumbers. The differences of some Raman bands between normalized mean spectra in Figure 3 are not visible since many Raman bands consist of multiple band components (see Figure 1(c)). However, the changes in band intensities can better be seen in the difference spectrum.

All spectra were fitted with a Voigt profile as mentioned previously. The vibrational band intensities and positions in wavenumbers were determined for the control and preeclamptic groups according to fit results. The vibrational bands exhibiting differences in intensities ( 𝑃 < 0 . 0 5 ) are given in Table 1 with tentative band assignments. Additionally, the possible contributions of the constituents to vibrational bands on the Raman spectra of serum were determined with the measurements of the Raman spectra of amino acids and other ingredients mentioned above. Their contributions to the vibrational bands were also represented in Table 1.

The vibrational bands in Table 1 were divided into two groups according to alterations in band intensities. As can be seen from the table, the intensities of 11 Raman bands decreased and eight Raman bands increased in preeclamptic samples when compared to the control group. The seven Raman bands with reduced band intensities are all protein bands presented above with possible main contributions from amino acids glycine, glutamine, valine, leucine, phenylalanine, tyrosine, and histidine. Other significant decreases in band intensities for preeclamptic samples are observed at the bands at 958, 1342, 1448, and 1679 cm−1. Regarding these bands, 958 and 1448 cm−1  have contributions from both amino acids and lipids. The band at 1342 cm−1 can be ascribed to CH bending of many amino acids, and the band at 1679 cm−1 can be assigned to one of the coupled modes of the polypeptide backbone, the so-called amid I band. Lipids contribute to all of the bands of increasing intensities and amino acids also contribute to most of these bands.

A PCA was performed by using five factors which were determined by selecting an eigenvalue number that gives minimum slope of the scree plot and by taking into account the total variance explained values under the threshold after testing the accuracy of the analysis by leave-one-out cross-validation. PCA was applied to all spectra which build up a matrix as an input. This gives scores and loadings as outputs. The scatter plot of the scores PC2 versus PC3 and spectra of their loadings in spectral range of 600 cm−1–1760 cm−1 with the assigned Raman bands (given in Table 1) are shown together in Figure 4. Among the first five most important component scores, PC 2 and PC 3 were found to be significant in the 99% confidence level. Therefore, PC 2 and PC 3 loadings were used to provide chemical information.

To determine a quantitative discrimination, we applied our LDA code on the first two most important PCA scores. Fisher’s criterion was used to obtain a decision line to discriminate two groups from each other which is superposed in Figure 4(a). The graph shows that the two groups are clearly separated from each other. The duplicate numbers indicate the samples (1–10 control and 11–19 preeclampsia) taken from the same person measured on different days as mentioned previously. The decision line demonstrates that two healthy samples from one person and four diseased samples from two patients were misclassified. The results of PCA-LDA analysis applied to all serum samples summarized in Table 2 with 78% sensitivity and 90% specificity.

PC 2 and PC 3 loadings are presented in Figures 4(b) and 4(c) superposed with control mean Raman spectrum. Two sharp Raman bands at 1156 cm−1 and at 1522 cm−1 belong to betacarotene Raman bands appear in Figure 4(b) since they scattered strongly in intensity especially for the preeclamptic patients. They are absent in Table 1 because these bands showed a slight change in average intensity according to the band component analysis. Moreover, their intensities were also found to be insignificant ( 𝑃 > 0 . 0 5 ) according to the Mann-Whitney U test. PC 2 loading mainly shows negative variations for lipid bands which are at 1300, 1439, and 1741 cm−1, while CH bending mode at 1340 and amide I band at around 1650 cm−1 exhibit positive variation. Although PC 3 shows more noisy contributions than PC 2, it exhibits the same trend like in PC 2 for protein and lipid bands except amide I band given in Figure 4(c). Protein bands at 1642, 1667 (amide I), and 699 cm−1 show negative variations, and the ones at 900, 940, 1209, 1342, and 1448 cm−1 show positive variations. On the other hand, lipid bands at 717, 1086, 1269, and 1439 cm−1 exhibit negative variation. Some of these bands do not appear in Table 1 since the significance values of those bands’ intensities are higher than 0.05.

4. Discussion

By means of band component analysis, it is observable that the intensities of 11 Raman bands decreased and eight Raman bands increased in the preeclamptic samples. Although the decreasing bands in intensities are probably related to the bands of the individual amino acids (e.g., Figure 2) given in Table 1, it is difficult to conclude that a certain amino acid is responsible for the corresponding decrease because of the complexity of the spectrum. However, protein-related bands in the spectra tend to decrease in intensity for preeclamptic patients. It is well known that proteins are vulnerable to the attack from free radicals and related oxidants [4, 51]. When attacked by free radicals, amino acids are converted to their oxidative derivatives [52] which may be the reason for the decrease in the protein band intensities.

A significant increase in the vibrational band intensities of preeclamptic samples was found at 699, 717, 1269, 1303, 1427, 1439, 1460, and 1667 cm−1 which are all related to lipids as given in Table 1.

In some studies, [1113] elevated lipid concentrations are found in preeclampsia, which are also consistent with our results. Increased lipid concentrations are known to be a major risk factor for cardiovascular diseases which is a known cause of endothelial dysfunction also seen in preeclampsia [1].

It can be seen from Table 1 that although the differences in increasing or decreasing band intensities are almost within the standard deviations, they exhibit significant differences in the intensities between the two groups ( 𝑃 < 0 . 0 5 ). However, many of the bands are formed by multiple amino acids, lipids, or other constituents in the serum. Therefore, the individual contributions of amino acids to these differences in band intensities remain unknown. Moreover, some weak bands in Table 1 may not be considered as reliable components of spectra since their intensities are almost at the level of noise.

The use of Raman spectroscopy combined with multivariate statistical methods gives rise to the classification between samples and helps to determine which variables change in disease states when compared to healthy samples. The linear discriminant analysis calculation after PCA exhibits a clear separation between the preeclamptic and control groups which can be seen in Figure 4(a).

It is to be remarked that each pair of data points belonging to the same sample is closely located in the graph. Although the pairs were measured on different days, they gave similar results. This means that systematic errors due to experimental conditions do not have any noticeable effect on the results.

We predicted preeclampsia with 78% specificity and 90% sensitivity by using multivariate analysis with Raman spectroscopy. The results obtained may be interpreted as low for the diagnosis of a disease.

PC 2 and PC 3 loadings spectra are consistent with the results of band component analysis.

These results may be improved by applying further multivariate techniques with larger numbers of samples. Nevertheless, Raman spectroscopy may provide a noninvasive method for studies exploring the pathophysiology of preeclampsia.

5. Conclusion

This study is the first Raman spectroscopic study on preeclampsia. When the control and preeclamptic groups were compared, alterations at some Raman bands of serum samples were observed. The contributions of some biomarkers that may lead to these changes were determined by using Raman spectroscopic measurements. These alterations on lipid bands and also on protein bands assigned to main contributions from particular amino acids may be investigated in further studies in respect of lipid bands and protein bands that focus on explaining the pathophysiology of the disease. They may also be tested in studies trying to predict the existence of preeclampsia. Although studies that use more samples may be required in order to confirm our results, our preliminary findings may offer a pioneering basis for further studies.

Acknowledgment

S. Seninak and U. Parlatan are grateful to the support offered by the research fund of Istanbul University, Project no. 4328 and 4330, respectively.