|
1.IntroductionNoninvasive diagnostics are one of the most important directions for the development of modern medicine. Recently, the interest has been focused on patients’ exhaled air research as a noninvasive diagnostic method for bronchopulmonary, cardiovascular, gastrointestinal, and other diseases.1,2 The basis of similar methods is related to the variation of concentrations of volatile organic compounds (VOCs) in exhaled air according to specific diseases.3 For example, it is ascertained that bronchial asthma exacerbation is characterized by ammonia () concentration increasing in exhaled air by 250 to 300 times.3,4 A high level of propane () in the exhaled air was identified for different clinical forms of pulmonary tuberculosis.5,6 The issue is significant because this approach, due to the sparse sampling, precludes pain and physical and emotional discomfort of the patient, the possibility of transmission of the blood-borne infections, and provides safety for the diagnostic studies. On one hand, noninvasive diagnostic methods can be used on an outpatient basis—this provides their widespread application; on the other hand, for patients at resuscitation departments, as the severity of the patient’s state is not a contraindication for their application. In this paper, we discuss the abilities of the methods of infrared (IR) laser spectroscopy and the principal component analysis (PCA) technique for a totally noninvasive express diagnostic of chronic obstructive pulmonary disease (COPD) on the basis of absorption spectra analysis of the patient’s exhaled air. 2.Techniques and MethodsThe method of laser photoacoustic spectroscopy (LPAS) is convenient for measuring the concentration of VOCs in the exhaled air because of the simplicity of its practical implementation, safety, cost-effectiveness, and extremely high sensitivity (minimal measured concentration for some chemicals at atmospheric conditions is about 1 ppb). LPAS has advantages in the detection of the gases which have overlapping of the absorption lines with the lasing lines (e.g., CO or -laser). LPAS is based on the generation of acoustic waves in a gas excited by the modulated laser beam at the wavelength corresponding to the absorption line of the VOCs in gaseous samples and on the measurement of the parameters of these acoustic waves using sensitive microphones. 2.1.LPAS Combined with CO2-Laser-lasers are one of the most suitable radiation sources for photoacoustic detection because of their narrow lasing line and commercial availability. Waveguide laser exited by radio frequency generator (144 MHz, 1 kW) with a wavelength tuning from 9.2 to was developed by Special Technologies Ltd. Selection of the necessary lasing lines was provided by diffraction grating. We realized gas analyzers based on the LPAS method related with the mentioned above -laser with intracavity (ILPA) and extracavity (LGA-2) photoacoustic detector location in cooperation with the Institute of Laser Physics SB RAS and Institute of Atmospheric Optics SB RAS. Both devices have high sensitivity (over 1 ppb level) and spectral resolutions (over ) at pulsed mode. Main technical parameters of the developed gas analyzers are shown in the Table 1. 2.2.LPAS Combined with Optical Parametric OscillatorFrequency conversion using optical parametric oscillator (OPO) is one of the effective ways to generate widely tunable coherent light in a spectral range from visible to mid-IR ranges. These laser sources play a particularly important role in the IR range, where VOCs have their fundamental absorption lines. This is because of the OPOs’ ability to provide continuous wavelength tuning over a wide spectral range7 that the PAS method combined with OPO allows concentration determination of a number of different gases. We developed the gas analyzer LaserBreeze based on the LPAS method and OPO with a tuning range from 2.5 to . The main technical characteristics of the LaserBreeze gas analyzer are shown in Table 2. We used two types of nonlinear elements in the optical scheme: periodically poled lithium niobate structure (PPLN) and mercury thiogallate crystal (HGS). Special cavities were designed for each element and an Nd:YLF laser (10 ns, 0.5 to 1.5 kHz, 1.5 mJ) was used as a pump source. The linewidth of the developed OPOs was 3 to . The average power of OPO based on PPLN structure was 20 mW (1700 Hz). The average power of OPO based on HGS crystal was 9 mW (900 Hz). The double-channel resonant photoacoustic cell was used for recording the absorption spectra of gaseous samples. The LaserBreeze gas analyzer is described in detail in Ref. 8. Table 1Main technical parameters of the developed gas analyzers.
Table 2Main technical characteristics of the LaserBreeze gas analyzer.
The main VOCs that can be detected by gas analyzers based on the LPAS method with a -laser and OPO are shown in Table 3. The sensitivity for measured gases was not lower than 1 ppb. Table 3Main VOCs.
Given in Table 3 wavelengths were chosen so that they are out of the absorption lines of other gases in the mixture or, if it is impossible, absorption of other gases is minimized. Moreover, for some gases, it is reasonable to carry out concentration measurements in two so-called spectral measurement channels:9 the wavelength of one spectral channel is close to or coincides with the center of one of the absorption lines of the gas (usually the most intense). The spectral channel is located on the edge of the absorption line. Due to the overlapping spectra of the individual gas components, the task of selecting the spectral measurement channels becomes complex and sometimes ambiguous.9 In this case, the use of special computational algorithms capable of selecting a set of spectral measurement channels where the errors of recovery of the gas concentrations would be minimal or close to the minimum. We used this approach for selecting the wavelengths at which the concentrations were measured. 2.3.Principal Component Analysis and Data PreprocessingThe measuring of VOCs’ concentrations in exhaled air is a promising tool for diagnostics in the future, but it should be pointed out that a significant part of the VOCs is not highly specific. For example, asthma causes the essential growth of exhaled NO and moderate growth of CO, COPD causes a small NO growth and essential growth of CO.3 Additionally, taking into account the individual variability in metabolism, it is obvious that for the diagnostics it is more expedient to use the “profile” of the set of VOCs or to directly profile the absorption spectrum of a breath sample as a “fingerprint” of the patient’s medical condition without component analysis of the sample. In this situation, various methods of data mining promise to be effective to analyze data.10,11 One of them is the PCA.12,13 The basic idea of PCA is to find the minimum number of new features that are enough for the recovery of the basic features by linear transformation, possibly with insignificant errors. PCA projects correlated variables into a lower number of uncorrelated variables called principal components (PCs). A specific feature of PCA is that the hidden connections and patterns that are typical for the investigated data set can be revealed. The mathematical background of PCA consists of the decomposition of initial experimental data two-dimensional (2-D) matrix into the form of a matrix product11 where is the quantity of samples of experimental data, is the quantity of the features of investigated objects, is the score’s matrix, is the loading’s matrix, is the residual’s matrix, and is the quantity of PCs. In our case, these features of the state under investigation are the set of absorption coefficients of the exhaled air sample in the laser source frequency detuning branch of the used gas analyzer.The loading’s matrix contains weight coefficients which characterize the contribution of features to a specific PC. The score’s matrix contains coordinates of the samples in the space of PCs. PCA is useful if . In this case, the method allows, first of all, to separate the most informative features of the state, or in other words, to reduce the dimensions of the feature space and to decrease noise, and second, to estimate the relative position of the studied objects in the reduced space of PCs. 3.Results and DiscussionThe experimental research was carried out according to the principles of Good Clinical Practices. The protocol of the research was approved by the Ethic Committee of the Siberian State Medical University (Tomsk, Russia), Ref. Number 2882 on 24 November 2011. All participants were informed about details of the research and signed “Informed agreement” for the actions carried out. The study involved 11 healthy nonsmoking volunteers (control group) and seven patients with COPD (target group). The COPD patients were males with verified diagnoses who passed treatment at the Pulmonological Division of the Regional State Autonomous Institution of Public Health “Municipal Clinical Hospital No. 3” (Tomsk, Russia). The average age of this group was 59.6 years. We did not included COPD patients with an unverified diagnosis, the presence of pneumonia, asthma, and other respiratory pathologies. The control group consisted of conventionally healthy nonsmoking male volunteers. Inclusion criteria were the absence of acute illness within two weeks prior to sample collection, without chronic pathologies of bronchopulmonary, cardiovascular, digestive, urinary and reproductive systems, and the absence of the factor “smoking” in anamnesis vitae. The average age in this group was 21.1 years. The procedure of exhaled air sampling was as follows. All samples were taken before eating or 2 h thereafter. The air was collected in standard test tubes. Prior to sampling, participants rinsed their mouth with running water. The study does not imply special cleaning of the oral cavity. Then the participant did some calm breaths through a sterile plastic tube into the test tube, which was then sealed with a sterile cotton wad. All exhaled air samples were analyzed using the LGA-2 and the LaserBreeze LPAS gas analyzers. Five scans of the absorption spectrum of each sample were recorded and averaged to reduce random errors. Most informative subranges of the measured profiles of the absorption spectra were determined by PCA. The criterion was the best spatial separation of the target group from the control group in the space of the PC. The results below are focused on the most informative subranges. The number of PC (in other words, the dimensions of the above-mentioned space) is usually chosen to describe at least 70% of the variation of initial data [it is the so-called explained variance (EV)].11 Here, initial data involved in the absorption spectra of all exhaled air samples. The value of EV in the used spectral subranges that is dependent on the quantity of the PC is presented in Table 4. According to the data from Table 4, the 2-D space of PC is enough to analyze profiles of the absorption spectra of the exhaled air samples. Table 4The dependence of the explained variance (EV) on the quantity of the used principal components (PCs).
To select the most informative set of absorption coefficients, we apply the method that is similar to the well-known “method of broken sticks” to the loading matrix.14 The results are shown in Table 5. Here, the initial quantity is the quantity of absorption coefficients which were contained in the definite spectral subrange before PCA application. Table 5The quantity of informative absorption coefficients.
According to the PCA, every sample is represented by the point in the space of PC. We used the freeware “ViDaExpert”15 to estimate the spatial distribution of the exhaled air samples in the space of the PC. The results of point estimates of the exhaled air absorption spectra profiles of COPD patients and healthy volunteers in the 9.2 to subrange are shown in Fig. 1. The distance between the point estimates on the plane of the PC characterized the difference in the absorption spectra profiles of participants. This is caused by variations in metabolism and, hence, is the difference in the VOCs profile of the samples. The similar results of point estimates of the measured spectra of the exhaled air of COPD patients and healthy volunteers from the control group using gas analyzer LaserBreeze in the most informative subranges from 2.59 to are shown in Figs. 2Fig. 3–4. Figures 2–4 show that the methods of IR LPAS and PCA allow separating patients with COPD and healthy nonsmoking volunteers. Taking into account that the VOCs absorption bands correspond to the chosen spectral ranges, we can assume that spatial separation of the target and control groups is probably caused by the difference of hydrocarbon content in the COPD patients’ breath and the healthy nonsmoking volunteers’ breath. This matches to the results obtained by means of other methods for the analysis of the exhaled air.2 For further analysis, the absorption spectrum profile of a breath sample is suggested to be recognized as a “fingerprint” of the medical state of a patient. In order to estimate the possibility of diagnostics based on such “fingerprints,” the algorithm of soft*independent modeling of class analogy (SIMCA) was applied. SIMCA classification includes two stages. The training stage. Each class of objects from the training set is independently modeled using PCA. In this result, the initial data are presented in the cloud form in the space of PCs. The coordinate origin is placed at the center of gravity of the cloud. Each object can be represented as the sum of two vectors: one lying in the cloud (projection) and another perpendicular to the first (residues). The average value (range) and deviation of the lengths of these vectors are the indicators belonging to this class of objects. The testing stage. The classification procedure is as follows. Each new object is projected onto the built space (cloud). The obtained range and deviation are compared with the critical levels specified in the training stage. In Table 6, we presented the examples of SIMCA classification using the profiles of the absorption spectra of breath samples in the range of 2.59 to for various sets of samples for the training and testing stages. Table 6The classification of the absorption spectra of exhaled air of patients with COPD and healthy volunteers by SIMCA in the range of 2.59 to 2.817 μm.
The results in Table 6 show that analysis of the profile of exhaled air absorption spectra in the IR region allows us to separate COPD patients from the control group with a high enough accuracy. 4.ConclusionWe described two types of laser photoacoustic gas analyzers which were developed by Special Technologies Ltd. for medical applications. Laboratory research of the exhaled air of patients with COPD and healthy nonsmoking volunteers was carried out at Siberian State Medical University (Russia) and in the Tomsk State University (Russia). The PCA method was used to select the most informative ranges of the absorption spectra of patients’ exhaled air in terms of the separation of the studied groups. It is shown that analysis of the profile of the exhaled air absorption spectrum allows identifying COPD patients in comparison to the control group. The most informative ranges of the absorption spectra of the COPD patients’ exhaled air and healthy nonsmoking volunteers’ exhaled air are 9.2 to 9.8, 2.59 to 2.817, and 3.272 to . The presented results are the base for the future construction of the classification rules for the noninvasive express diagnostics methods. There are two ways for classification rules construction. The first is to use the profile of the absorption spectrum of a breath sample as a “fingerprint” of the patient’s medical state. Another one consists of two steps: first, to carry out component analysis of breath samples for various groups, then to define the profile of the set of informative VOCs as a “fingerprint” of the patient’s medical state. AcknowledgmentsThe work was carried out with partial financial support of the FCPIR Contract No. 14.578.21.0082 (ID RFMEFI57814X0082). ReferencesD. Smith and A. Amann, Breath Analysis for Clinical Diagnosis and Therapeutic Monitoring, World Scientific, Singapore
(2005). Google Scholar
D. Smith and A. Amann, Volatile Biomarkers: Non-Invasive Diagnosis in Physiology and Medicine, 1st ed.Elsevier, UK
(2013). Google Scholar
S. A. Kharitonov and P. J. Barnes,
“Exhaled markers of pulmonary disease,”
Am. J. Respir. Crit. Care Med., 163
(7), 1693
–1722
(2001). http://dx.doi.org/10.1164/ajrccm.163.7.2009041 AJCMED 1073-449X Google Scholar
M. Yamara,
“Exhaled carbon monoxide levels during treatment of acute asthma,”
Eur. Respir. J., 13 757
–760
(1999). http://dx.doi.org/10.1034/j.1399-3003.1999.13d10.x ERJOEI 0903-1936 Google Scholar
S. Kwiatkowska,
“Elevated exhalation of hydrogen peroxide and circulating IL-18 in patients with pulmonary tuberculosis,”
Respir. Med., 101
(3), 574
–580
(2007). http://dx.doi.org/10.1016/j.rmed.2006.06.015 RMEDEY 0954-6111 Google Scholar
O. B. Pikas,
“Effects of alcoholic beverages on the fatty acid spectrum of the expired air condensate lipids in patients with tuberculosis of the respiratory organs,”
Lik. Sprava, 7–8 30
–33
(2000). LISPEC Google Scholar
F. K. Tittel, D. Richter and A. Fried,
“Mid-infrared laser applications in spectroscopy,”
Solid-State Mid-Infrared Laser Sources, 445
–516 Springer-Verlag, Berlin, Heidelberg
(2003). http://dx.doi.org/10.1007/3-540-36491-9_11 Google Scholar
A. A. Karapuzikov et al.,
“LaserBreeze gas analyzer for noninvasive diagnostics of air exhaled by patients,”
Phys. Wave Phenom., 22
(3), 189
–196
(2014). http://dx.doi.org/10.3103/S1541308X14030054 1541-308X Google Scholar
V. I. Kozintsev et al., Laser Photo-Acoustic Analysis of Multicomponent Gaseous Mixture, Science and Education, Moscow
(2003). http://dx.doi.org/10.7463/0612.0368798 Google Scholar
M. Phillips et al.,
“Volatile organic compounds in breath as markers of lung cancer: a cross-sectional study,”
Lancet, 353
(9168), 1930
–1933
(1999). http://dx.doi.org/10.1016/S0140-6736(98)07552-7 LANCAO 0140-6736 Google Scholar
D. Poli et al.,
“Exhaled volatile organic compounds in patients with non-small cell lung cancer: cross sectional and nested short-term follow-up study,”
Respir. Res., 6
(1), 71
(2005). http://dx.doi.org/10.1186/1465-9921-6-71 RREEBZ 1465-9921 Google Scholar
A. L. Pomerantsev and O. Ye Rodionova,
“Concept and role of extreme objects in PCA/SIMCA,”
J. Chemom., 28
(5), 429
–438
(2014). http://dx.doi.org/10.1002/cem.v28.5 JOCHEU 0886-9383 Google Scholar
A. D. Wilson and M. Baietto,
“Advances in electronic-nose technologies developed for biomedical applications,”
Sensors, 11 1105
–1176
(2011). http://dx.doi.org/10.3390/s110101105 SNSRES 0746-9462 Google Scholar
R. Cangelosi and A. Goriely,
“Component retention in principal component analysis with application to cDNA microarray data,”
Biol. Direct, 2 21
(2007). http://dx.doi.org/10.1186/1745-6150-2-2 BDIIBV 1745-6150 Google Scholar
A. N. Gorban and A. Y. Zinovyev,
“Principal graphs and manifolds,”
Handbook of Research on Machine Learning Applications and Trends: Algorithms, Methods and Techniques, 28
–59 Information Science Reference, IGI Global, Hershey, Pennsylvania
(2009). http://dx.doi.org/10.4018/978-1-60566-766-9 Google Scholar
Biography |