Next Article in Journal
Improved Diazo-Transfer Reaction for DNA-Encoded Chemistry and Its Potential Application for Macrocyclic DEL-Libraries
Next Article in Special Issue
Searching for Potential Markers of Glomerulopathy in Urine by HS-SPME-GC×GC TOFMS
Previous Article in Journal
Solution and Solid State Studies of Urea Derivatives of DITIPIRAM Acting as Powerful Anion Receptors
Previous Article in Special Issue
The Influence of Smoking Status on Exhaled Breath Profiles in Asthma and COPD Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Needle Trap Device-GC-MS for Characterization of Lung Diseases Based on Breath VOC Profiles

by
Fernanda Monedeiro
1,
Maciej Monedeiro-Milanowski
1,
Ileana-Andreea Ratiu
1,2,3,
Beata Brożek
4,
Tomasz Ligor
1,3,* and
Bogusław Buszewski
1,3
1
Interdisciplinary Centre of Modern Technologies, Nicolaus Copernicus University in Toruń, 4 Wileńska St., 87-100 Toruń, Poland
2
“Raluca Ripan” Institute for Research in Chemistry, Babeş-Bolyai University, 30 Fântânele St., RO-400294 Cluj-Napoca, Romania
3
Department of Environmental Chemistry and Bioanalytics, Faculty of Chemistry, Nicolaus Copernicus University in Toruń, 7 Gagarina St., 87-100 Toruń, Poland
4
Department of Lung Diseases, Provincial Polyclinic Hospital in Toruń, 4 Krasińskiego St., 87-100 Toruń, Poland
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(6), 1789; https://doi.org/10.3390/molecules26061789
Submission received: 24 February 2021 / Revised: 16 March 2021 / Accepted: 19 March 2021 / Published: 22 March 2021

Abstract

:
Volatile organic compounds (VOCs) have been assessed in breath samples as possible indicators of diseases. The present study aimed to quantify 29 VOCs (previously reported as potential biomarkers of lung diseases) in breath samples collected from controls and individuals with lung cancer, chronic obstructive pulmonary disease and asthma. Besides that, global VOC profiles were investigated. A needle trap device (NTD) was used as pre-concentration technique, associated to gas chromatography-mass spectrometry (GC-MS) analysis. Univariate and multivariate approaches were applied to assess VOC distributions according to the studied diseases. Limits of quantitation ranged from 0.003 to 6.21 ppbv and calculated relative standard deviations did not exceed 10%. At least 15 of the quantified targets presented themselves as discriminating features. A random forest (RF) method was performed in order to classify enrolled conditions according to VOCs’ latent patterns, considering VOCs responses in global profiles. The developed model was based on 12 discriminating features and provided overall balanced accuracy of 85.7%. Ultimately, multinomial logistic regression (MLR) analysis was conducted using the concentration of the nine most discriminative targets (2-propanol, 3-methylpentane, (E)-ocimene, limonene, m-cymene, benzonitrile, undecane, terpineol, phenol) as input and provided an average overall accuracy of 95.5% for multiclass prediction.

1. Introduction

Respiratory diseases are conditions which affect the airways and other structures of the lungs and they are represented by lung cancer, asthma, tuberculosis, chronic obstructive pulmonary disease (COPD) and pneumonia, being the leading causes of mortality and morbidity globally. Smoking or exposure to secondhand smoke is the main risk factor associated to most of respiratory diseases, with current smokers 11 times more likely to develop lung cancer compared to non-smokers [1]. Globally, respiratory diseases affect 1 billion people and account for 7% of all deaths worldwide. Nevertheless, even considering that lung cancer is one of the leading causes of death worldwide, COPD and asthma are predominant lung diseases that represent a burden on society in terms of health care costs [2]. The diagnosis of asthma or COPD is usually made by non-invasive techniques based on spirometry, however lung cancer is often diagnosed in late stages, due to the lack of noticeable clinical manifestations, or because these can be easily associated with other symptoms. This fact may reduce the chance of applying a timely and effective treatment. Currently used diagnostic methods for respiratory diseases includes physical examination followed by a set of chemical, imaging, endoscopic and immunological procedures [3]. Because different lung diseases are characterized by inflammation and other correspondent symptoms, direct assessment of airways may be applied, by using invasive procedures such as: computer tomography, bronchoscopy, bronchoalveolar lavage or biopsy. These are costly, time consuming and/or invasive procedures [2]. Consequently, a simple, reliable, low-cost and non-invasive test, able to achieve the diagnosis in real time (minutes up to hours), using a mere sample of exhaled breath in highly desired.
Therefore, fast detection and characterization of volatile organic compounds (VOCs) emitted from different biological matrices (breath, sweat, saliva, plasma, tissues, exudates, urine, etc.) as a tool for diagnosis was approached [4,5,6,7,8,9,10,11]. Breath tests are minimally invasive procedures, which are more easily accepted by the patients. An exhaled breath sample consists of VOCs and the breath aerosol [12]. Breath consists of almost 3000 compounds which are present in different combinations and quantities. Consequently, not only specific biomarkers, but the global VOC profile can be potentially associated to a characteristic fingerprint for each disease [2]. Exhaled breath is largely composed of nitrogen, oxygen, carbon dioxide, water, and inert gases. Trace components—volatile substances that are generated in the body or absorbed from the environment—present in the nmol/L–pmol/L (ppb volume—ppt volume) range make up the rest of the breath. The endogenous VOCs are generated by the cellular biochemical processes of the body, hence VOCs existent in human breath can reflect endogenous metabolic processes which occur in the tissues. VOCs-patterns in exhaled breath have been associated with various respiratory diseases such as cancer, asthma, COPD, cystic fibrosis, tuberculosis, etc. [13,14]. Breath samples are probably the most adequate to reach the rapid diagnosis of respiratory conditions, once substances from surrounding blood vessels and tissue can be exchanged in the alveoli and be available in the exhaled air. A large number of VOCs has been reported in scientific literature as markers of various diseases, as well as bacterial infections. These compounds can be divided into different chemical groups [15,16,17]: saturated hydrocarbons (stable end products of lipid peroxidation) and unsaturated hydrocarbons (e.g., from mevalonic pathway of cholesterol synthesis) [6,16], alcohols (which can be addressed as oxidized products of hydrocarbons and their precursors) [16], aldehydes (associated with inflammatory processes, resulting from lipid peroxidation) [5,18], ketones (products of fatty acid decarboxylation processes in the liver, associated to a diet rich in proteins and fat) [16], aromatic VOC–typically related to exogenous sources such as tobacco smoke and pollution [19], sulfur-containing compounds generated by incomplete metabolism of methionine in the transamination pathway and also associated with bacterial activity [20,21]), and nitrogen-containing compounds (such as ammonia, dimethylamine and trimethylamine, derived e.g., when conversion to urea is limited due to an impairment of liver function) [17].
Nowadays, gas chromatography–mass spectrometry (GC-MS) is considered a gold standard for VOC analysis [22]. Solid phase microextraction (SPME) or sorption tubes followed by thermal desorption are the most frequently used pre-concentration techniques in breath analyses. A prominent sampling tool is the needle trap device (NTD), which consists of a sorbent material packed inside a needle, working as an extraction trap [23]. This solventless technique provides exhaustive extraction and has potential for laboratory automation [24,25]. In the present work, NTD was used as extraction technique, followed by GC-MS analysis. VOCs were analyzed in breath samples belonging to healthy controls and patients with lung cancer, asthma and COPD, in an attempt to develop a classification model able to discriminate between these lung diseases, which have in common inflammatory processes in the lungs. In this sense, besides the assessment of global VOC profiles, 29 target compounds previously reported as potential biomarkers of the referred respiratory diseases were also investigated and quantified in breath samples.
The present study describes the non-invasive assessment of asthma, COPD and lung cancer, based on breath analysis of VOCs. Once all of these are lung diseases involving inflammatory mechanisms, the applied design of data analysis intended to find specific VOC patterns able to provide discrimination between these illnesses. The comparison between self-annotated discriminating features and compounds reported by literature as indicators of lung diseases represents an original approach for the validation of candidate biomarkers. The outline of the work presents the application of NTD for the determination of VOCs in breath. The found results aim to support the implementation of breath analysis to the clinical practice, as an accurate and reliable diagnostic tool.

2. Results and Discussion

2.1. Calibration Method and Quantitation of Analytes

Table S1 presents information regarding calibration method, while Table 1 displays data concerning the quantitation of analytes in breath samples. Obtained limits of quantitation (LOQs) ranged from 0.003 (3-methylpentane, 2-butanone, toluene, isododecane, 1,2,4-trimethylbenzene, (E)-ocimene, limonene, m-cymene and benzonitrile) to 6.21 ppbv (tridecane). Higher limits were obtained for heavier and more polar analytes, which also displayed wider linearity ranges. Lower limits were associated to compounds with higher volatility, a factor that seemed to contribute for their more efficient recovery, besides their expected greater stability in samples. Relative standard deviation (RSD%) did not exceed 10%, demonstrating that the proposed method provided adequate reproducibility. In general, suitability of NTD for preconcentration of analytes in gas mixtures could be inferred. Among the targets, isoprene and 1-propanol were found in each breath sample. Styrene, decane and phenol were observed in lowest frequency of appearance. Ethanol, isoprene and acetoin were the targets which occurred in higher concentrations in all sample’s cohorts. Carry-over effect was not observed, indicating that there is no influence of previously analyzed samples on the current ones.

2.2. VOCs Detected in Breath

Regarding the obtained VOC global profiles, a total number of 112 different VOCs were detected. The VOCs most frequently observed in the samples were hydrocarbons, alcohols, aldehydes and ketones. A graph displaying the distribution of VOCs according to the functional groups in profiles belonging to the different studied groups is presented in Figure 1a. In general, the number of compounds belonging to each of the chemical classes seems to be proportional when evaluating the different studied conditions, however, some particularities of the qualitative composition of each group of profiles can be evidenced. Lung cancer and COPD profiles appear to be associated to a greater variety of compounds (103 and 95 detected VOCs, respectively), while asthma profiles are composed by smaller number of compounds (84 detected VOCs). An increased number of hydrocarbons is observed in the VOC composition in breath of lung cancer patients. Moreover, samples from patients with lung cancer and COPD appear related to a greater variety of aldehydes (12 and 11, respectively, against 9 found in healthy). This observation can be due to the fact that hydrocarbons and aldehydes are frequently reported as the most characteristic products of oxidative stress induced by inflammatory process [26,27,28].
A matrix displaying number and percentage of overlapping VOCs in the acquired profiles is presented in Figure 1b. By the content of coincident compounds, the level of similarity regarding the qualitative composition of breath samples can be inferred. In this sense, lung cancer and COPD profiles, display the greater similarity between each other, followed by the VOC profile of lung cancer and healthy individuals, while asthma breath samples present to be the most distinct in terms of composition.

2.3. Differential Distribution of VOCs

Principal component analysis (PCA) was performed intending to identify relationships and existing patterns within datasets. Peak area data regarding the global VOC profiles was used as input for generation of the score plot depicted in Figure 2a, in which 78.04% of variance was represented by the two first principal components. When using as input the calculated concentration values of the 29 preselected analytes in samples, the plot presented in Figure 2b is produced. In this case, 79.72% of total variance was described by the components 1 and 2. In both cases, around 80% of the total variance can be assigned to the observed distribution. Although both score plots indicate a discrimination between control cases and remaining samples, a clearer grouping can be observed when considering the global profile, once in Figure 1b control samples appear confined to an isolated cluster. Still, in both situations, the lack of a distinct grouping according to each of the investigated conditions demonstrates that other factors play a relevant role in the observed pattern of distribution of VOCs. This can be mainly related to the variability in the nature and extension of the involved pathophysiological mechanisms, inherent to the different lung diseases. Therefore, the usage of supervised approaches is essential to achieve the classification of samples in agreement with the related diagnosis.
A volcano plot was built in order to present found discriminating features when considering obtained global VOCs profiles. In Figure 3a the overall trend of the detected VOCs (variables) is graphically represented. The variables located above the dashed line refer to the compounds which displayed greater statistically relevant changes in their responses when compared to the control group. The variables located along the y-axis correspond to VOCs absent in the healthy group and detected solely in positive samples. In the left part of the plot are displayed compounds with decreased responses in the positive samples, while in the right side of the plot are displayed VOCs presenting an increased response in samples of diseased. The VOCs located towards the top of the graph expressed the greatest statistical significance. The names of the most discriminative components are exhibited in the plot.
Figure 3b presents a bar graph showing the distribution of all compounds classified as discriminant features, considering as criteria p ≤ 0.05. Most of the compounds which displayed significant alteration in their responses when compared to those presented in the healthy group belong to the class of hydrocarbons, followed by alcohols and aldehydes. In lung cancer profiles, a greater number of discriminating VOCs was verified (41 compounds). For asthma and COPD samples, 26 and 24 altered VOCs were indicated, respectively. As presented in Figure 1b, around 92, 88 and 74% of the compounds observed in lung cancer, COPD and asthma samples, respectively, were shown to be conserved in the healthy group profiles. This indicates that the differential abundance of VOCs in samples is determinant to discriminate between samples’ group, once the similarity between the qualitative profiles belonging to the four studied groups is not so divergent. Such observation highlights the importance of validated quantitative assays’ application regarding breath samples for diagnosis purposes.
Few compounds presented a more expressive incidence within the group of active smokers’ individuals, thus possibly being ascribed as products of cigarette smoke. 1,3-Cyclopentadiene was identified solely in this group, in 40% of the samples; 2,5-Dimethylfuran was detected in 80% of samples from active smokers, which represented around 73% of its total incidence across samples. Other substances commonly related to tobacco smoke composition, such as benzene and toluene [29], did not present a specific distribution within samples of smokers, probably because these can be originated from other various sources.
With respect to the VOCs found altered, acetonitrile is typically present in cigarette smoke, although also present in automobile exhaust and other anthropogenic emissions [30]. Considering that most of the enrolled subjects were not smokers, differentiated levels of this substance would not be expected. However, together with the decreasing trend observed for p-xylene, the reduction in the abundances of such compounds in positive group can be an indicative of diminished ability of elimination of exogenous through exhaled air, or a consequence of the augmented activity of cytochrome P450 isoforms documented in lung cancer [31], which could be responsible for the rapid metabolization of inhaled compounds in the lungs.
The two main lung cancer types are small-cell lung carcinoma (SCLC) and non-small-cell lung carcinoma (NSCLC). Two hypothesis involve SCLC histogenesis: the first assumes that SCLC derives from cells of the diffuse endocrine system, i.e., the amine precursor uptake decarboxylation (APUD)-system, the second suggests this type of lung cancer originates from the endodermbronchial lining [32,33]. Adenocarcinoma (NSCLC subtype) arises from glandular cells of bronchial mucosa, whereas squamous lung cancer origins from the modified bronchial epithelial cells and adenosquamous carcinoma contains two types of cells: squamous cells (thin, flat cells that line certain organs) and gland-like cells. Finally, large cell (undifferentiated) carcinoma originates from epithelial cells of the lung [32]. The origin and nature of the malignant cells is crucial for different treatment strategies. Tumor tissue releases different protein biomarkers according to subtype of cancer. The same concerns different types or amounts of certain VOCs secreted by various malignant part of cell. The oxidation of fatty acids present in the cell membranes is pointed out as the source of VOCs associated to oxidative stress condition. The mentioned process is initiated by the reactive oxygen species (ROS) which are found in increased levels in inflamed tissues [17,28].
Due to the ROS activity, mechanism chain reactions occur, with radicals tending to be stabilized through alpha and beta scissions [34], leading to the formation of a variety of shorter chain fatty acids, alkanes, alkenes, alcohols and aldehydes. In addition, formed compounds can be subjected to other reactions, aiming their transformation into smaller and more polar molecules [35]. Cancer cells are characterized by their enhanced metabolism and altered functions in several biochemical pathways [36]. Therefore, metabolite profile consisting of a greater variety of compounds may be expected. Hexane can be possibly formed during the oxidation of oleic acid [34], while can be addressed as an exogenous substance as well. Hexane showed decreased abundance in cancer and COPD samples. This fact can be explained by three hypotheses: impaired excretion through exhalation [37], enhanced conversion of the specie into oxidized forms [38] and favoring of alternative mechanism, which gives rise to different products, during lipid oxidation associated to oxidative stress particular to the referred conditions.
1-Pentanol can be interpreted as a pentane oxidation product, caused by cytochrome P450, and recognized as a metabolite of reactive oxygen species reactions with omega-6 fatty acids [26]. Methyl ketones such as 2-dodecanone can be formed by the decarboxylation of β-keto acids during the metabolism of fatty acids [39]. Nonanal can be also formed by different mechanisms during ROS attack on oleic acid from cell membranes [34]. Medium-chain branched alkanes, such as 2-methyldecane and 4-methyloctane, were pointed out by previous works as oxidative stress indicators [40,41]. However, their generation by human organism due to the oxidation of lipids is questionable, as cell membranes contain only linear chain lipids [26].
Branched alkanes can be originated from microbial lipids, mostly produced in the fatty acid pathway of bacteria, by using amino acids as precursor molecules which are submitted to elongation in this biochemical path [42,43]. Considering this, the occurrence of methylated branched alkanes in breath could be connected with bacterial activity. Alternatively, these could be products of transformation/degradation of prenyl molecules in organism, a mechanism that also remains undescribed. Aromatic species, such as p-xylene (decreased in COPD) and 1,2,4-trimethylbenzene (increased in COPD), are frequently addressed as pollutants, although also possibly formed by bacterial shikimate and related pathways [44].
Regarding the 29 compounds belonging the set of selected targets, 15 of them presented themselves as discriminating features (p < 0.05) when assessing solely controls against positive samples, all of them displaying increased concentration in the positive group. They were 2-propanol, 2-methylpentane, 3-methylpentane, 1-propanol, 2-butanone, styrene, isododecane, 1,2,4-trimethylbenzene, (E)-ocimene, m-cymene, phenol, undecane, dodecane, terpineol and tridecane. However, as demonstrated in the next section, compounds other than these displayed usefulness in the characterization of studied groups, presenting themselves as discriminating variables related to disease type. A combination of mechanisms involved in carcinogenesis, inflammatory processes and microbiota activity–which develop important role in pathogenesis of several diseases, may play a part in the alterations observed for certain compounds in breath samples.
The propionic acid formed during microbial fermentation and the propionyl-CoA generated during amino acids degradation enters in the propanoate metabolism, which takes place in the mitochondria and comprehend a series of reactions coupled with other pathways related to cell energetics. In microorganisms, 1-propanol is a product of propanoyl-CoA transformation [45], while 2-propanol can be formed by the reduction of acetone produced during the synthesis of ketone bodies [46]. 2-Butanone is a secondary ketone, therefore its origin can be associated to the β-oxidation of fatty acids. The acetyl-CoA units generated in this process fuel the citric acid cycle, supplying energy generation [47]. Terpenoids are very diverse natural products synthetized by plants, but also by bacteria. These metabolites are associated to the mevalonate and deoxyxylulose phosphate pathways [48,49]. Although their biosynthesis in human so far remain unknown, studies have reported terpenoid derivatives as potential cancer indicators. Considering this, increased concentration of compounds such as (E)-ocimene, m-cymene and terpineol can either be a consequence of deficient metabolic function impairing proper elimination of these substances coming from diet [50], an indicative of specific bacterial activity, or even a product of transformation of isoprenoids derivatives due to the dysregulated mevalonate pathway in human during carcinogenesis [51].
Isododecane is known as a synthetic chemical with several applications in the industry [52], without any identified biosynthetic pathway so far. Styrene is a constituent of polymers, nevertheless, there is evidence that some microorganisms can produce styrene using phenylalanine as precursor molecule [53]. On the other hand, phenol is often reported as product of bacterial catabolism of aromatic amino acid species previously documented as elevated in gastroesophageal neoplasms [54].
Their formation of the n-alkanes undecane, dodecane and tridecane can be related to the oxidation of lipids, more precisely, a formed alkoxyl radical undergoes scission, generating an alkyl radical which abstracts a hydrogen atom, turning into a stable alkane [17,26]. 2-Methylpentane and 3-methylpentane are other branched species possibly derive from the oxidation of branched chain fatty acids generated by bacteria.

2.4. Diagnosis Prediction–Global Profiles

Most of the studies comprising the detection of diseases based on VOC analysis in biological samples compare paired data from healthy and diseased groups. Many of the compounds addressed as candidate biomarkers by literature are explained as produced by oxidative stress–a process promoted by typical inflammatory immune responses and thus non-specific. In this sense, illnesses sharing common etiological and pathological processes may play a part as confounding factors when a specific diagnosis is intended. For this reason, the present and following sections of the manuscript were dedicated to the development of statistical models able to identify and discriminate specific VOC patterns, allowing simultaneous differentiation of the studied lung diseases.
A random forest (RF) analysis was conducted on global profiles data, aiming to classify obtained VOC fingerprints into the four investigated categories. Variance importance was assessed based on the mean decrease Gini when one of the questioned variables is removed from a preliminary RF model. Gini impurity can be interpreted as the chance of a case sampled randomly to be incorrectly classified in relation to a given class, thus being related to the purity of cases within a tree node [55]. Therefore, greater decreases in this measurement indicate greater importance of a given variable. The resulting plot is presented in Figure 4a, the compounds are ranked from the most essential to those less relevant for the obtaining of homogenous classes. The 12 most important variables were assigned to compose the RF final model, the selected compounds appear depicted as the gray diamonds, in the upper part of the graph. The intention was to obtain the greater model overall accuracy as possible, including a minimum number of features.
Predict probabilities of a case of the validation set to belong to a class were provided by RF modeling. The receiver operating characteristic (ROC) curves presenting the ability of the model to predict a certain condition are showed in Figure 4b, information on parameters regarding classification performance are presented in Table 2. It can be observed that class recognition was performed with at least 93% of sensitivity and 87.5% of specificity for lung cancer, asthma and healthy groups. Regarding the later mentioned groups, prediction with accuracy above 87% was achieved. The lower prediction capability obtained in case of COPD (67%). An exemplary decision tree, from the 1000 generated during modeling, is presented in Figure 4c.

2.5. Diagnosis Prediction–Target Analysis

In this section, in accordance with the criteria described in the Material and Methods section and empiric observations drawn from multinomial logistic regression (MLR) performance using different set of variables, 2-propanol, 3-methylpentane, (E)-ocimene, limonene, m-cymene, benzonitrile, undecane, terpineol, phenol were the compounds selected to build the MLR final model. A clearer depiction of variables distribution according to their importance can be observed in Figure 5. Table 3 presents information regarding the developed model, which, when applied to the train and test datasets provided 100% and 90.5% of accuracy, respectively (average overall accuracy = 95.3%).
In MLR, coefficients can be multiplied by the quantitative inputs for the calculation of probabilities of a case to belong to a specific condition. Equation (1) presents the model regression equation, where ln[P/(1 − P)] represents the log-odds pertinent to a specific disease, β0 is the intercept and β1…k are the coefficients provided by the MLR model, referring to the variables X (in the case, the concentration of the selected targets). A case for which the calculated probabilities are greater than 50%, can be assigned as belonging to that class.
ln[P/(1 − P)] = β 0 + β 1 · X 1 + … + β k · X k
The numerical coefficients provided by MLR can be interpreted as weights, or the contribution of these variables to the designated classes. Positive coefficients are related to compounds with increased response when comparing to the reference class (“Healthy”), while negative coefficients are associated to targets which were present in lower concentrations in positive samples. In a closer interpretation, the coefficients express multinomial log-odds. For example, assuming that all other variables remain constant, an increase of one unit in the concentration of 2-propanol multiplies the odds of a sample belonging to the asthma group instead of healthy group by 0.56. On the other hand, an increase in one unit of (E)-ocimene concentration in breath implies the log-odds of COPD to decrease by 64.27, in an assumption that the remaining variables are kept constant. Considering this, increased levels of limonene and m-cymene are characteristic from samples of asthma patients, while increased level of undecane and decreased concentrations of benzonitrile are observed for breath of individual with lung cancer; Moreover, greater concentrations of phenol and lower concentrations of m-cymene are particularly observed in samples from COPD patients. Values fitted for the train set and predictions performed by MLR method solely for the validation set were used as input to build ROC curves (Figure 6a–d and Figure 6e, respectively). Values of area under the curve (AUC) presented in Figure 6a–d represent the probability of samples belonging to a given group to be classified as the state condition. For each class specified in the model, AUC was 1.0, meaning that 100% of sensitivity and specificity was obtained. On the other hand, cases not assigned as the state variable provided AUC ≤ 0.5 (curves below random guessing line). When considering the performance of the model on the test data, an overall accuracy of 91% was obtained, resulting in an average accuracy of 95.5% when both evaluated sets are considered. Detailed information regarding MLR performance is presented in Table 4.

3. Materials and Methods

3.1. Apparatus and Standards

The analyses were conducted on a model 6890 A gas chromatograph coupled with a 5975 Inert XL MSD (Agilent Technologies, Waldbronn, Germany). Inlet temperature was kept at 260 °C and carrier gas (helium 6.0) flow was set at 2.2 mL min−1. A DB-624 capillary column (Agilent) 60 m × 0.32 mm × 1.8 µm was used. The oven temperature program was as follows: initial temperature was 35 °C (held for 3 min), ramped to 50 °C, then 75 °C, 200 °C and finally 240 °C, at rates of 3 °C min−1, 5 °C min−1, 15 °C min−1 and 10 °C min−1, respectively. The last temperature was kept for 15 min, resulting in a run time of 41.33 min. Full scan spectra were acquired within a range of 30–300 m/z, at electron ionization (EI) of 70 eV. The ion source and transfer line were set to 250 °C. Chromatographic data acquisition was performed using MSD ChemStation E.02.00.493 software (Agilent). Compounds identification was performed using NIST05 mass spectra library. Each peak was searched manually, including baseline subtraction and averaging over a peak. Forward match quality of at least 750/1000 was applied as the lower match threshold.
Needle trap device 700-60d-PXC (PDMS + Carbopack X + Carboxen 1000) was purchased from PAS Technology (Magdala, Germany). The air pump flow was conducted by a sampling case model SC-B (PAS Technology), designed for controlled air sampling. Prior first use, NTDs were conditioned in a heated conditioner (PAS Technology) at 300 °C under helium flow (1 bar), for 30 min, in order to remove VOC’s contaminations from sorbent. One liter-Tedlar bags were obtained from SKC (Eighty Four, PA, USA).
Chemicals used as standards (2-methylbutane, pentane, ethanol, isoprene, 2-propanol, 2-methylpentane, 3-methylpentane, 1-propanol, methylcyclopentane, 2-butanone, benzene, acetoin, toluene, ethylbenzene, p-xylene, styrene, decane, 6-methyl-2-heptanone, isododecane, 1,2,4-trimethylbenzene, ocimene, D-limonene, m-cymene, benzonitrile, phenol, undecane, dodecane, terpineol and tridecane) were purchased from Sigma-Aldrich (St. Louis, MO, USA), all with purity > 98%.

3.2. Breath Collection

The study was approved by the local Ethics Committee of Collegium Medicum in Bydgoszcz (No. KB 621/2016–25.10.2016). Individuals aged over 18, with positive clinical diagnosis for lung cancer (non-small cell lung cancer, subtype: adenocarcinoma) (n = 16), chronic obstructive pulmonary disease (n = 12) and asthma (n = 8) were recruited at the Department of Lung Diseases of the Provincial Polyclinic Hospital in Toruń. Samples from enrolled cancer patients were obtained before any medical intervention (such as neoadjuvant therapies or surgery).
Individuals were refrained to eat, drink or smoke 2 h prior sample collection. No special dietary regimes were applied. All individuals gave informed consent to participation in the study. The patients completed a questionnaire describing their age, gender and current smoking status (active smokers, non-smokers). Samples of mixed alveolar breath gas (alveolar and dead space gas) were collected in Tedlar bags with parallel collection of ambient air at the same room. Breath samples were obtained after approximately after 10 min rest in the same ambient. Each subject provided breath samples using a disposable plastic straw connected to the Tedlar bag.
Control samples (n = 20) were collected from healthy individuals aged over 18 years, without any history of positive diagnosis for cancer or respiratory diseases, who were not suffering from any other inflammatory disease. All samples were analyzed within 2–3 h after collection–this timeframe was considered adequate to avoid the interference of gas composition losses [56]. In the total, 56 breath samples were collected. Information regarding enrolled volunteers is summarized in Table 5 (details regarding presented significance probabilities are described in the section “Data analysis and chemometrics approaches”).
Tedlar bags involved in sample collection and calibration experiments were daily treated with several cycles of cleaning, each consisting of consecutively filling and evacuating argon 5.0 from the bag. Afterwards, the bags filled with argon were kept in an oven at 65 °C. The content bag was tested before breath sampling, by means of GC-MS, in order to verify the effectiveness of cleaning procedure.

3.3. Selection of Targets

The compounds selected as targets were VOCs previously reported as potential breath biomarkers of lung cancer, COPD and asthma, in accordance with previous studies on this theme. A literature search was performed in the electronic database Web of Science Core Collection (from Clarivate Analytics; Philadelphia, PA, USA), as well as Google Scholar. The searched terms were: “volatile organic compounds”, “gas chromatography”, “biomarker”, “lung cancer”, “COPD” and “asthma”, considering a time span from 1999 to 2016. The indexed literature is presented in the Supplementary Material (Table S2) [57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80].

3.4. Calibration Procedure

Gas mixtures of the analytes were prepared by injection of 1 µL of liquid standards into 1 L glass bulb (Supelco, Bellefonte, PA, USA) previously evacuated. Methanol HPLC was used for the preparation of 50:50 (v/v) dilution of acetoin, phenol and terpineol, which are solids at room temperature. After the complete vaporization of the liquids, the interior of the bulb was equilibrated with nitrogen, generating a gas mixture of the compounds of interest. Using a gas-tight syringe, different volumes of the stock gas solution were transferred to Tedlar bags filled with 1 L of nitrogen, in order to obtain the desired concentrations.
The concentrations were calculated in terms of part per billion per volume of analyte (ppbv), taking in consideration their molar volume. Six calibration levels were used in the construction of calibration curves, all analyzed in triplicates. The limit of detection (LOD) was defined as the lowest detectable concentration of analyte, considering a signal-to-noise (S/N) ratio of at least 3. LOQ was considered as the lowest concentration of analyte with imprecision of at least 15%, considering a minimum S/N value equal to 10. Calibration was conducted by linear regression analysis, using the obtained experimental data. Linearity was evaluated by the method of least squares and reported as the coefficient of determination (R2). Linearity was confirmed for values of R² above 0.99. Inter-assay imprecision was assessed from the evaluation of assays in triplicate, these were expressed in terms of relative standard deviation (RSD%). Reported RSD% values are the average of imprecision calculated for lower (LOQs), medium (5.17–17.25 ppbv) and high calibration levels (9.52–3452.0 ppbv)—which concentrations varied depending on the linearity range displayed by the analyte.

3.5. Sample Extraction

Prior to sample extraction process, NTDs were conditioned for 10 min, at 300 °C, under helium 6.0 flow (Air Products, Warsaw, Poland). Samples in Tedlar bags were drawn through the air pump, at a flow rate of 30 mL min−1. The fixed volume of 50 mL was sampled from each bag. Once extraction was complete, the loaded NTD was desorbed into GC inlet port for 2 min.

3.6. Data Analysis and Chemometrics Approaches

For the building of main dataset, area of peaks belonging to ambient air samples were subtracted from respective samples obtained from patients. Evaluation of normality of distributions was conducted using Kolmogorov-Smirnov test. Differences between volunteers’ ages was assessed by t-test. Principal component analysis was performed in order to evaluate data patterns regarding sample’s group. Mann-Whitney test was used to indicate VOCs which presented statistically relevant differences in their responses in the studied groups, p ≤ 0.05 was considered as the relevance criteria. For the above cited tests, IBM SPSS Statistics v. 24 software (IBM Corp., Armonk, NY, USA) was used. The following approaches were executed in R environment, using RStudio console v. 1.1.463 (RStudio, PBC, Boston, MA, USA). Significant differences between the proportions of volunteers assigned to a certain group were assessed by the test of equal or given proportions, employing the R function “prop.test”. For chemometrics approaches, the packages “gplots”, “RandomForest”, “caret”, “ROCR” were employed. Random forest is a machine learning method based on recognition of latent patterns within global VOC profiles. RF was conducted in order to develop a multiclass model, able to distinguish between studied conditions. RF input consisted of peak table data converted into binary entries–once this algorithm was dedicated to non-quantified data, this format of dataset was considered as more appropriate than to express RF outcome in terms of peak area. Variable importance plots were assessed for selection of variables to be included in the model. Half of the data was randomly selected to compose the training set (bootstrap sampling method) and the remaining data was applied in the validation process. Receiver operating characteristic curves were built based on calculated probabilities output from RF modeling. Ultimately, the development of a classificatory model based solely on target compounds was aimed, for that, variables (targets) were selected according to their discriminative potential between all four investigated conditions. The criteria comprised most unique targets which presented higher discriminative relevance when considering a given condition against all others. MLR was performed using the package “nnet”, employing the data regarding quantitation of the selected targets in analyzed samples. This multiclass categorical method performs a linear combination of features, allowing prediction through the calculated probabilities of an input (set of features’ values) to belong to each class specified in the model. Sixty percent of the data regarding quantitation of targets in the samples was randomly addressed as the training set, while the remaining data was addressed to a testing set. “Healthy” group was defined as the reference class. ROC curves were prepared based on the predictions computed by developed MLR model for fitted values and test data.

4. Conclusions

The developed NTD-GC-MS method was demonstrated to be suitable for the determination of target VOCs in breath samples, providing considerably low limits of detection and quantitation, as well as appropriate reproducibility. From the 29 targets selected from literature, more than half of them presented significant differentiated responses among control and positive groups – found discriminating features were 2-propanol, 2-methylpentane, 3-methylpentane, 1-propanol, 2-butanone, styrene, isododecane, 1,2,4-trimethylbenzene, (E)-ocimene, m-cymene, phenol, undecane, dodecane, terpineol and tridecane, limonene and benzonitrile (which proved to serve for further differentiation between diseases). Built statistical models (using both self-annotated discriminating variables and quantified targets) aimed to simultaneously classify VOC profiles into lung cancer, COPD or asthma cases. Both classification models (RF and MLR), provided an overall accuracy above 80%. The distinction between VOC profiles related to clinical conditions involving concomitant molecular mechanisms is extremely relevant in order to assess cofounding aspects of breath analysis diagnosis. In this sense, machine learning tools and other mathematical models can be useful to identify disease-specific latent patterns. Cross-validated studies, comparing candidate biomarkers found by different research groups by means of different techniques, are essential for a future implementation of breath screening tests in a clinical setting. Such an approach can also enable a focused investigation of the pathways involved in the modulation of these potential biomarkers, as well as it can contribute to the establishment of optimized analysis protocols.

Supplementary Materials

The following are available online, Table S1: Data regarding calibration method of gas mixtures (LOD = limit of detection, LOQ = limit of quantitation, ppbv = part per billion per volume, R² = determination coefficient, RSD = relative standard deviation). Table S2: References which reported the targets selected in the present study as potential biomarker of lung diseases in breath samples, where: LC–lung cancer; COPD–chronic obstructive pulmonary disease.

Author Contributions

Conceptualization, F.M., T.L. and I.-A.R.; methodology, F.M. and T.L.; software, F.M.; validation, F.M. and M.M.-M.; formal analysis, F.M.; investigation, F.M. and M.M.-M.; resources, T.L. and Beata Brożek (B.B.); data curation, F.M., M.M.-M., I.-A.R. and Beata Brożek (B.B.); writing—original draft preparation, F.M. and I.-A.R.; writing—review and editing, F.M., M.M.-M. and T.L.; visualization, F.M.; supervision, T.L. and Bogusław Buszewski (B.B.); project administration, T.L. and Bogusław Buszewski (B.B.); funding acquisition, T.L. and Bogusław Buszewski (B.B.). All authors have read and agreed to the published version of the manuscript.

Funding

This research was financed by The National Centre for Research and Development (Warsaw, Poland) in frame of Polish-Turkish bilateral project “A comparative study of volatile organic compound biomarkers in breath and urine samples collected from Polish and Turkish communities for monitoring of various respiratory diseases” (POLTUR2/4/2018).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the local Ethics Committee of Collegium Medicum in Bydgoszcz (No. KB 621/2016–25.10.2016).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data is contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Varghese, C.; Troisi, G.; Schotte, K.; Prasad, V.M.; St Claire, S.M. World No Tobacco Day 2019 puts the spotlight on lung health. J. Thorac. Dis. 2019, 11, 2639–2642. [Google Scholar] [CrossRef] [PubMed]
  2. Van de Kant, K.D.; Van der Sande, L.J.; Jöbsis, Q.; Van Schayck, O.C.; Dompeling, E. Clinical use of exhaled volatile organic compounds in pulmonary diseases: A systematic review. Respir. Res. 2012, 13, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Hashoul, D.; Haick, H. Sensors for detecting pulmonary diseases from exhaled breath. Eur. Respir. Rev. 2019, 28, 190011. [Google Scholar] [CrossRef] [Green Version]
  4. Sun, X.; Shao, K.; Wang, T. Detection of volatile organic compounds (VOCs) from exhaled breath as noninvasive methods for cancer diagnosis. Anal. Bioanal. Chem. 2016, 408, 2759–2780. [Google Scholar] [CrossRef] [PubMed]
  5. Monedeiro, F.; Dos Reis, R.B.; Peria, F.M.; Sares, C.T.G.; De Martinis, B.S. Investigation of sweat VOC profiles in assessment of cancer biomarkers using HS-GC-MS. J. Breath Res. 2020, 14, 026009. [Google Scholar] [CrossRef]
  6. Monedeiro, F.; Milanowski, M.; Ratiu, I.-A.; Zmysłowski, H.; Ligor, T.; Buszewski, B. VOC Profiles of Saliva in Assessment of Halitosis and Submandibular Abscesses Using HS-SPME-GC/MS Technique. Molecules 2019, 24, 2977. [Google Scholar] [CrossRef] [Green Version]
  7. Ratiu, I.-A.; Ligor, T.; Monedeiro, F.; Al-Suod, H.; Bocos-Bintintan, V.; Szeliga, J.; Jackowski, M.; Buszewski, B. “Features of infected versus uninfected chemical profiles released from human exudates ”. Stud. Univ. Babeș-Bolyai Chem. 2019, 64, 207–216. [Google Scholar] [CrossRef]
  8. Silva, C.L.; Passos, M.; Câmara, J.S. Investigation of urinary volatile organic metabolites as potential cancer biomarkers by solid-phase microextraction in combination with gas chromatography-mass spectrometry. Br. J. Cancer 2011, 105, 1894–1904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Buszewski, B.; Kęsy, M.; Ligor, T.; Amann, A. Human exhaled air analytics: Biomarkers of diseases. Biomed. Chromatogr. 2007, 21, 553–566. [Google Scholar] [CrossRef]
  10. Filipiak, W.; Mochalski, P.; Filipiak, A.; Ager, C.; Cumeras, R.; E. Davis, C.; Agapiou, A.; Unterkofler, K.; Troppmair, J. A Compendium of Volatile Organic Compounds (VOCs) Released By Human Cell Lines. Curr. Med. Chem. 2016, 23, 2112–2131. [Google Scholar] [CrossRef] [Green Version]
  11. Amann, A.; Miekisch, W.; Schubert, J.; Buszewski, B.; Ligor, T.; Jezierski, T.; Pleil, J.; Risby, T. Analysis of Exhaled Breath for Disease Detection. Annu. Rev. Anal. Chem. 2014, 7, 455–482. [Google Scholar] [CrossRef]
  12. Ratiu, I.-A.; Bocos-Bintintan, V.; Monedeiro, F.; Milanowski, M.; Ligor, T.; Buszewski, B. An Optimistic Vision of Future: Diagnosis of Bacterial Infections by Sensing Their Associated Volatile Organic Compounds. Crit. Rev. Anal. Chem. 2019, 1–12. [Google Scholar] [CrossRef]
  13. Schmidt, K.; Podmore, I. Current Challenges in Volatile Organic Compounds Analysis as Potential Biomarkers of Cancer. J. Biomarkers 2015, 2015, 1–16. [Google Scholar] [CrossRef] [Green Version]
  14. Mazzone, P.J. Analysis of Volatile Organic Compounds in the Exhaled Breath for the Diagnosis of Lung Cancer. J. Thorac. Oncol. 2008, 3, 774–780. [Google Scholar] [CrossRef] [Green Version]
  15. Ratiu, I.-A.; Ligor, T.; Bocos-Bintintan, V.; Buszewski, B. Mass spectrometric techniques for the analysis of volatile organic compounds emitted from bacteria. Bioanalysis 2017, 9, 1069–1092. [Google Scholar] [CrossRef] [PubMed]
  16. Nardi-Agmon, I.; Peled, N. Exhaled breath analysis for the early detection of lung cancer: Recent developments and future prospects. Lung Cancer Targets Ther. 2017, 8, 31–38. [Google Scholar] [CrossRef] [Green Version]
  17. Miekisch, W.; Schubert, J.K.; Noeldge-Schomburg, G.F. Diagnostic potential of breath analysis—Focus on volatile organic compounds. Clin. Chim. Acta 2004, 347, 25–39. [Google Scholar] [CrossRef] [PubMed]
  18. Milanowski, M.; Monedeiro, F.; Złoch, M.; Ratiu, I.-A.; Pomastowski, P.; Ligor, T.; De Martinis, B.S.; Buszewski, B. Profiling of VOCs released from different salivary bacteria treated with non-lethal concentrations of silver nitrate. Anal. Biochem. 2019, 578, 36–44. [Google Scholar] [CrossRef] [PubMed]
  19. Vishinkin, R.; Haick, H. Nanoscale Sensor Technologies for Disease Detection via Volatolomics. Small 2015, 11, 6142–6164. [Google Scholar] [CrossRef]
  20. Ratiu, I.A.; Bocos-Bintintan, V.; Patrut, A.; Moll, V.H.; Turner, M.; Thomas, C.L.P. Discrimination of bacteria by rapid sensing their metabolic volatiles using an aspiration-type ion mobility spectrometer (a-IMS) and gas chromatography-mass spectrometry GC-MS. Anal. Chim. Acta 2017, 982, 209–217. [Google Scholar] [CrossRef]
  21. Ratiu, I.-A.; Ligor, T.; Bocos-Bintintan, V.; Al-Suod, H.; Kowalkowski, T.; Rafińska, K.; Buszewski, B. The effect of growth medium on an Escherichia coli pathway mirrored into GC/MS profiles. J. Breath Res. 2017, 11, 036012. [Google Scholar] [CrossRef] [Green Version]
  22. Mametov, R.; Ratiu, I.-A.; Monedeiro, F.; Ligor, T.; Buszewski, B. Evolution and Evaluation of GC Columns. Crit. Rev. Anal. Chem. 2019, 1–24. [Google Scholar] [CrossRef]
  23. Dobrzyńska, E.; Buszewski, B. Needle-trap device for the sampling and determination of chlorinated volatile compounds. J. Sep. Sci. 2013, 36, 3372–3378. [Google Scholar] [CrossRef] [PubMed]
  24. Lord, H.L.; Zhan, W.; Pawliszyn, J. Fundamentals and applications of needle trap devices. Anal. Chim. Acta 2010, 677, 3–18. [Google Scholar] [CrossRef] [PubMed]
  25. Filipiak, W.; Filipiak, A.; Ager, C.; Wiesenhofer, H.; Amann, A. Optimization of sampling parameters for collection and preconcentration of alveolar air by needle traps. J. Breath Res. 2012, 6, 027107. [Google Scholar] [CrossRef]
  26. Frank Kneepkens, C.M.; Lepage, G.; Roy, C.C. The potential of the hydrocarbon breath test as a measure of lipid peroxidation. Free Radic. Biol. Med. 1994, 17, 127–160. [Google Scholar] [CrossRef]
  27. Fuchs, P.; Loeseken, C.; Schubert, J.K.; Miekisch, W. Breath gas aldehydes as biomarkers of lung cancer. Int. J. Cancer 2010, 126, 2663–2670. [Google Scholar] [CrossRef]
  28. Hakim, M.; Broza, Y.Y.; Barash, O.; Peled, N.; Phillips, M.; Amann, A.; Haick, H. Volatile Organic Compounds of Lung Cancer and Possible Biochemical Pathways. Chem. Rev. 2012, 112, 5949–5966. [Google Scholar] [CrossRef] [PubMed]
  29. Charles, S.M.; Batterman, S.A.; Jia, C. Composition and emissions of VOCs in main- and side-stream smoke of research cigarettes. Atmos. Environ. 2007, 41, 5371–5384. [Google Scholar] [CrossRef]
  30. Kischkel, S.; Miekisch, W.; Sawacki, A.; Straker, E.M.; Trefz, P.; Amann, A.; Schubert, J.K. Breath biomarkers for lung cancer detection and assessment of smoking related effects—Confounding variables, influence of normalization and statistical algorithms. Clin. Chim. Acta 2010, 411, 1637–1644. [Google Scholar] [CrossRef]
  31. Elfaki, I.; Mir, R.; Almutairi, F.M.; Abu Duhier, F.M. Cytochrome P450: Polymorphisms and roles in cancer, diabetes and atherosclerosis. Asian Pacific J. Cancer Prev. 2018, 19, 2057–2070. [Google Scholar] [CrossRef]
  32. Zamay, T.; Zamay, G.; Kolovskaya, O.; Zukov, R.; Petrova, M.; Gargaun, A.; Berezovski, M.; Kichkailo, A. Current and Prospective Protein Biomarkers of Lung Cancer. Cancers 2017, 9, 155. [Google Scholar] [CrossRef] [Green Version]
  33. Karachaliou, N.; Pilotto, S.; Lazzari, C.; Bria, E.; De Marinis, F.; Rosell, R. Cellular and molecular biology of small cell lung cancer: An overview. Transl. Lung Cancer Res. 2016, 5, 2–15. [Google Scholar] [CrossRef] [PubMed]
  34. Schaich, K.M.; Shahidi, F.; Zhong, Y.; Eskin, N.A.M. Lipid Oxidation. In Biochemistry of Foods; Eskin, N.A., Shahidi, F., Eds.; Elsevier: Cambridge, UK, 2013; pp. 419–478. ISBN 9780122423529. [Google Scholar]
  35. Altomare, D.F.; Di Lena, M.; Porcelli, F.; Trizio, L.; Travaglio, E.; Tutino, M.; Dragonieri, S.; Memeo, V.; De Gennaro, G. Exhaled volatile organic compounds identify patients with colorectal cancer. Br. J. Surg. 2013, 100, 144–150. [Google Scholar] [CrossRef] [PubMed]
  36. Hanahan, D.; Weinberg, R.A. Hallmarks of Cancer: The Next Generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Phillips, M.; Herrera, J.; Krishnan, S.; Zain, M.; Greenberg, J.; Cataneo, R.N. Variation in volatile organic compounds in the breath of normal humans. J. Chromatogr. B Biomed. Sci. Appl. 1999, 729, 75–88. [Google Scholar] [CrossRef]
  38. Mutti, A.; Falzoi, M.; Lucertini, S.; Arfini, G.; Zignani, M.; Lombardi, S.; Franchini, I. n-Hexane metabolism in occupationally exposed workers. Occup. Environ. Med. 1984, 41, 533–538. [Google Scholar] [CrossRef] [Green Version]
  39. Forney, F.W.; Markovetz, A.J. The biology of methyl ketones. J. Lipid Res. 1971, 12, 383–395. [Google Scholar] [CrossRef]
  40. Abaffy, T.; Duncan, R.; Riemer, D.D.; Tietje, O.; Elgart, G.; Milikowski, C.; DeFazio, R.A. Differential Volatile Signatures from Skin, Naevi and Melanoma: A Novel Approach to Detect a Pathological Process. PLoS ONE 2010, 5, e13813. [Google Scholar] [CrossRef] [Green Version]
  41. Phillips, M.; Cataneo, R.N.; Cummin, A.R.C.; Gagliardi, A.J.; Gleeson, K.; Greenberg, J.; Maxfield, R.A.; Rom, W.N. Detection of Lung Cancer With Volatile Markers in the Breatha. Chest 2003, 123, 2115–2123. [Google Scholar] [CrossRef] [Green Version]
  42. Bos, L.D.J.; Sterk, P.J.; Schultz, M.J. Volatile Metabolites of Pathogens: A Systematic Review. PLoS Pathog. 2013, 9, e1003311. [Google Scholar] [CrossRef] [Green Version]
  43. Egge, H.; Murawski, U.; Ryhage, R.; György, P.; Chatranon, W.; Zilliken, F. Minor constitutents of human milk IV: Analysis of the branched chain fatty acids. Chem. Phys. Lipids 1972, 8, 42–55. [Google Scholar] [CrossRef]
  44. Gosset, G. Production of aromatic compounds in bacteria. Curr. Opin. Biotechnol. 2009, 20, 651–658. [Google Scholar] [CrossRef] [PubMed]
  45. Xue, J.; Murrieta, C.M.; Rule, D.C.; Miller, K.W. Exogenous or l-Rhamnose-Derived 1,2-Propanediol Is Metabolized via a pduD-Dependent Pathway in Listeria innocua. Appl. Environ. Microbiol. 2008, 74, 7073–7079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ramachandriya, K.D.; Wilkins, M.R.; Delorme, M.J.M.; Zhu, X.; Kundiyana, D.K.; Atiyeh, H.K.; Huhnke, R.L. Reduction of acetone to isopropanol using producer gas fermenting microbes. Biotechnol. Bioeng. 2011, 108, 2330–2338. [Google Scholar] [CrossRef] [PubMed]
  47. Houten, S.M.; Violante, S.; Ventura, F.V.; Wanders, R.J.A. The Biochemistry and Physiology of Mitochondrial Fatty Acid β-Oxidation and Its Genetic Disorders. Annu. Rev. Physiol. 2016, 78, 23–44. [Google Scholar] [CrossRef] [Green Version]
  48. Ludwiczuk, A.; Skalicka-Woźniak, K.; Georgiev, M.I. Terpenoids. In Pharmacognosy; Badal, S., Delgoda, R., Eds.; Elsevier: London, UK, 2017; pp. 233–266. ISBN 9780128020999. [Google Scholar]
  49. Audrain, B.; Farag, M.A.; Ryu, C.M.; Ghigo, J.M. Role of bacterial volatile compounds in bacterial biology. FEMS Microbiol. Rev. 2015, 39, 222–233. [Google Scholar] [CrossRef] [Green Version]
  50. Friedman, M.I.; Preti, G.; Deems, R.O.; Friedman, L.S.; Munoz, S.J.; Maddrey, W.C. Limonene in expired lung air of patients with liver disease. Dig. Dis. Sci. 1994, 39, 1672–1676. [Google Scholar] [CrossRef]
  51. Mullen, P.J.; Yu, R.; Longo, J.; Archer, M.C.; Penn, L.Z. The interplay between cell signalling and the mevalonate pathway in cancer. Nat. Rev. Cancer 2016, 16, 718–731. [Google Scholar] [CrossRef] [PubMed]
  52. Johnson, W.; Bergfeld, W.F.; Belsito, D.V.; Hill, R.A.; Klaassen, C.D.; Liebler, D.; Marks, J.G.; Shank, R.C.; Slaga, T.J.; Snyder, P.W.; et al. Safety Assessment of Isoparaffins as Used in Cosmetics. Int. J. Toxicol. 2012, 31, 269S–295S. [Google Scholar] [CrossRef]
  53. Fujiwara, R.; Noda, S.; Tanaka, T.; Kondo, A. Styrene production from a biomass-derived carbon source using a coculture system of phenylalanine ammonia lyase and phenylacrylic acid decarboxylase-expressing Streptomyces lividans transformants. J. Biosci. Bioeng. 2016, 122, 730–735. [Google Scholar] [CrossRef] [PubMed]
  54. Wiggins, T.; Kumar, S.; Markar, S.R.; Antonowicz, S.; Hanna, G.B. Tyrosine, Phenylalanine, and Tryptophan in Gastroesophageal Malignancy: A Systematic Review. Cancer Epidemiol. Biomarkers Prev. 2015, 24, 32–38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pal, M. Random forest classifier for remote sensing classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  56. Steeghs, M.M.L.; Cristescu, S.M.; Harren, F.J.M. The suitability of Tedlar bags for breath sampling in medical diagnostic research. Physiol. Meas. 2007, 28, 73–84. [Google Scholar] [CrossRef] [Green Version]
  57. Bajtarevic, A.; Ager, C.; Pienz, M.; Klieber, M.; Schwarz, K.; Ligor, M.; Ligor, T.; Filipiak, W.; Denz, H.; Fiegl, M.; et al. Noninvasive detection of lung cancer by analysis of exhaled breath. BMC Cancer 2009, 9, 348. [Google Scholar] [CrossRef] [Green Version]
  58. Poli, D.; Carbognani, P.; Corradi, M.; Goldoni, M.; Acampa, O.; Balbi, B.; Bianchi, L.; Rusca, M.; Mutti, A. Exhaled volatile organic compounds in patients with non-small cell lung cancer: Cross sectional and nested short-term follow-up study. Respir. Res. 2005, 6, 71. [Google Scholar] [CrossRef] [Green Version]
  59. Fens, N.; Zwinderman, A.H.; Van der Schee, M.P.; De Nijs, S.B.; Dijkers, E.; Roldaan, A.C.; Cheung, D.; Bel, E.H.; Sterk, P.J. Exhaled Breath Profiling Enables Discrimination of Chronic Obstructive Pulmonary Disease and Asthma. Am. J. Respir. Crit. Care Med. 2009, 180, 1076–1082. [Google Scholar] [CrossRef]
  60. Ligor, M.; Ligor, T.; Bajtarevic, A.; Ager, C.; Pienz, M.; Klieber, M.; Denz, H.; Fiegl, M.; Hilbe, W.; Weiss, W.; et al. Determination of volatile organic compounds in exhaled breath of patients with lung cancer using solid phase microextraction and gas chromatography mass spectrometry. Clin. Chem. Lab. Med. 2009, 47, 550–560. [Google Scholar] [CrossRef] [PubMed]
  61. Ulanowska, A.; Kowalkowski, T.; Trawińska, E.; Buszewski, B. The application of statistical methods using VOCs to identify patients with lung cancer. J. Breath Res. 2011, 5, 046008. [Google Scholar] [CrossRef]
  62. Phillips, C.O.; Syed, Y.; Mac Parthaláin, N.; Zwiggelaar, R.; Claypole, T.C.; Lewis, K.E. Machine learning methods on exhaled volatile organic compounds for distinguishing COPD patients from healthy controls. J. Breath Res. 2012, 6, 036003. [Google Scholar] [CrossRef] [PubMed]
  63. Rudnicka, J.; Kowalkowski, T.; Ligor, T.; Buszewski, B. Determination of volatile organic compounds as biomarkers of lung cancer by SPME–GC–TOF/MS and chemometrics. J. Chromatogr. B 2011, 879, 3360–3366. [Google Scholar] [CrossRef]
  64. Fens, N.; Roldaan, A.C.; Van der Schee, M.P.; Boksem, R.J.; Zwinderman, A.H.; Bel, E.H.; Sterk, P.J. External validation of exhaled breath profiling using an electronic nose in the discrimination of asthma with fixed airways obstruction and chronic obstructive pulmonary disease. Clin. Exp. Allergy 2011, 41, 1371–1378. [Google Scholar] [CrossRef]
  65. Smolinska, A.; Klaassen, E.M.M.; Dallinga, J.W.; Van de Kant, K.D.G.; Jobsis, Q.; Moonen, E.J.C.; Van Schayck, O.C.P.; Dompeling, E.; van Schooten, F.J. Profiling of Volatile Organic Compounds in Exhaled Breath As a Strategy to Find Early Predictive Signatures of Asthma in Children. PLoS ONE 2014, 9, e95668. [Google Scholar] [CrossRef] [Green Version]
  66. Phillips, M.; Gleeson, K.; Hughes, J.M.B.; Greenberg, J.; Cataneo, R.N.; Baker, L.; McVay, W.P. Volatile organic compounds in breath as markers of lung cancer: A cross-sectional study. Lancet 1999, 353, 1930–1933. [Google Scholar] [CrossRef]
  67. Phillips, M.; Cataneo, R.N.; Ditkoff, B.A.; Fisher, P.; Greenberg, J.; Gunawardena, R.; Kwon, C.S.; Tietje, O.; Wong, C. Prediction of breast cancer using volatile biomarkers in the breath. Breast Cancer Res. Treat. 2006, 99, 19–21. [Google Scholar] [CrossRef]
  68. Gaspar, E.M.; Lucena, A.F.; Duro da Costa, J.; Chaves das Neves, H. Organic metabolites in exhaled human breath—A multivariate approach for identification of biomarkers in lung disorders. J. Chromatogr. A 2009, 1216, 2749–2756. [Google Scholar] [CrossRef] [PubMed]
  69. Gaida, A.; Holz, O.; Nell, C.; Schuchardt, S.; Lavae-Mokhtari, B.; Kruse, L.; Boas, U.; Langejuergen, J.; Allers, M.; Zimmermann, S.; et al. A dual center study to compare breath volatile organic compounds from smokers and non-smokers with and without COPD. J. Breath Res. 2016, 10, 026006. [Google Scholar] [CrossRef]
  70. Ibrahim, B.; Basanta, M.; Cadden, P.; Singh, D.; Douce, D.; Woodcock, A.; Fowler, S.J. Non-invasive phenotyping using exhaled volatile organic compounds in asthma. Thorax 2011, 66, 804–809. [Google Scholar] [CrossRef] [Green Version]
  71. Martinez-Lozano Sinues, P.; Meier, L.; Berchtold, C.; Ivanov, M.; Sievi, N.; Camen, G.; Kohler, M.; Zenobi, R. Breath Analysis in Real Time by Mass Spectrometry in Chronic Obstructive Pulmonary Disease. Respiration 2014, 87, 301–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Song, G.; Qin, T.; Liu, H.; Xu, G.-B.; Pan, Y.-Y.; Xiong, F.-X.; Gu, K.-S.; Sun, G.-P.; Chen, Z.-D. Quantitative breath analysis of volatile organic compounds of lung cancer patients. Lung Cancer 2010, 67, 227–231. [Google Scholar] [CrossRef] [PubMed]
  73. Dragonieri, S.; Schot, R.; Mertens, B.J.A.; Le Cessie, S.; Gauw, S.A.; Spanevello, A.; Resta, O.; Willard, N.P.; Vink, T.J.; Rabe, K.F.; et al. An electronic nose in the discrimination of patients with asthma and controls. J. Allergy Clin. Immunol. 2007, 120, 856–862. [Google Scholar] [CrossRef] [Green Version]
  74. Oguma, T.; Nagaoka, T.; Kurahashi, M.; Kobayashi, N.; Yamamori, S.; Tsuji, C.; Takiguchi, H.; Niimi, K.; Tomomatsu, H.; Tomomatsu, K.; et al. Clinical contributions of exhaled volatile organic compounds in the diagnosis of lung cancer. PLoS ONE 2017, 12, e0174802. [Google Scholar] [CrossRef]
  75. Chen, X.; Xu, F.; Wang, Y.; Pan, Y.; Lu, D.; Wang, P.; Ying, K.; Chen, E.; Zhang, W. A study of the volatile organic compounds exhaled by lung cancer cells in vitro for breath diagnosis. Cancer 2007, 110, 835–844. [Google Scholar] [CrossRef]
  76. Cazzola, M.; Segreti, A.; Capuano, R.; Bergamini, A.; Martinelli, E.; Calzetta, L.; Rogliani, P.; Ciaprini, C.; Ora, J.; Paolesse, R.; et al. Analysis of exhaled breath fingerprints and volatile organic compounds in COPD. COPD Res. Pract. 2015, 1, 7. [Google Scholar] [CrossRef] [Green Version]
  77. Van de Kant, K.D.G.; Van Berkel, J.J.B.N.; Jobsis, Q.; Lima Passos, V.; Klaassen, E.M.M.; Van der Sande, L.; Van Schayck, O.C.P.; De Jongste, J.C.; Van Schooten, F.J.; Derks, E.; et al. Exhaled breath profiling in diagnosing wheezy preschool children. Eur. Respir. J. 2013, 41, 183–188. [Google Scholar] [CrossRef] [PubMed]
  78. Van Berkel, J.J.B.N.; Dallinga, J.W.; Möller, G.M.; Godschalk, R.W.L.; Moonen, E.J.; Wouters, E.F.M.; Van Schooten, F.J. A profile of volatile organic compounds in breath discriminates COPD patients from controls. Respir. Med. 2010, 104, 557–563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Van Vliet, D.; Smolinska, A.; Jöbsis, Q.; Rosias, P.P.R.; Muris, J.W.M.; Dallinga, J.W.; Van Schooten, F.J.; Dompeling, E. Association between exhaled inflammatory markers and asthma control in children. J. Breath Res. 2016, 10, 016014. [Google Scholar] [CrossRef] [PubMed]
  80. Peng, G.; Hakim, M.; Broza, Y.Y.; Billan, S.; Abdah-Bortnyak, R.; Kuten, A.; Tisch, U.; Haick, H. Detection of lung, breast, colorectal, and prostate cancers from exhaled breath using a single array of nanosensors. Br. J. Cancer 2010, 103, 542–551. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) VOCs distribution according to main chemical classes, in profiles belonging to the different studied groups, the contoured box displays the total number of compounds found in each group; (b) Similarity matrix displaying number and percentage of overlapping VOCs in the acquired profiles.
Figure 1. (a) VOCs distribution according to main chemical classes, in profiles belonging to the different studied groups, the contoured box displays the total number of compounds found in each group; (b) Similarity matrix displaying number and percentage of overlapping VOCs in the acquired profiles.
Molecules 26 01789 g001
Figure 2. PCA plots using as input (a) VOCs’ responses in global profile analysis, (b) responses of the targets quantified in the samples (triangles = control samples; squares = positive samples).
Figure 2. PCA plots using as input (a) VOCs’ responses in global profile analysis, (b) responses of the targets quantified in the samples (triangles = control samples; squares = positive samples).
Molecules 26 01789 g002
Figure 3. (a) Volcano plot displaying the most discriminating features, in terms of fold change (x-axis) and statistical relevance (y-axis), in which the dashed line represents the point of y-axis in which p = 0.01; (b) Graph of distribution of number of all compounds assigned as discriminating features, according to disease and chemical class (significance criteria: p ≤ 0.05).
Figure 3. (a) Volcano plot displaying the most discriminating features, in terms of fold change (x-axis) and statistical relevance (y-axis), in which the dashed line represents the point of y-axis in which p = 0.01; (b) Graph of distribution of number of all compounds assigned as discriminating features, according to disease and chemical class (significance criteria: p ≤ 0.05).
Molecules 26 01789 g003
Figure 4. (a) Variable importance plot in terms of mean decrease Gini (node purity), obtained in the first training of RF model. Diamonds refer to VOCs selected for generation of the final classificatory model; (b) ROC curves based on RF’s final model output regarding the test set, using a panel of 12 VOCs; (c) Example decision tree produced by RF analysis, in which obtained accuracy was 81% (AS = asthma, CA = lung cancer, COPD = chronic obstructive pulmonary disease, H = healthy).
Figure 4. (a) Variable importance plot in terms of mean decrease Gini (node purity), obtained in the first training of RF model. Diamonds refer to VOCs selected for generation of the final classificatory model; (b) ROC curves based on RF’s final model output regarding the test set, using a panel of 12 VOCs; (c) Example decision tree produced by RF analysis, in which obtained accuracy was 81% (AS = asthma, CA = lung cancer, COPD = chronic obstructive pulmonary disease, H = healthy).
Molecules 26 01789 g004
Figure 5. Plot of −(log 10) of p values when applying Mann-Whitney test for specific classes: (a) lung cancer, (b) COPD or (c) asthma, against all other conditions. Dashed line represents where p ≤ 0.05. Variables represented by triangle shape icon were those included in MLR final model.
Figure 5. Plot of −(log 10) of p values when applying Mann-Whitney test for specific classes: (a) lung cancer, (b) COPD or (c) asthma, against all other conditions. Dashed line represents where p ≤ 0.05. Variables represented by triangle shape icon were those included in MLR final model.
Molecules 26 01789 g005
Figure 6. ROC curves generated from fitted values (train set) created by MLR model, labelling (a) healthy, (b) lung cancer, (c) COPD and (d) asthma as the state variables; (e) ROC curves generated from predictions computed by the MLR model for the test set. Colored numerals refer to values of AUC obtained for each depicted curve.
Figure 6. ROC curves generated from fitted values (train set) created by MLR model, labelling (a) healthy, (b) lung cancer, (c) COPD and (d) asthma as the state variables; (e) ROC curves generated from predictions computed by the MLR model for the test set. Colored numerals refer to values of AUC obtained for each depicted curve.
Molecules 26 01789 g006
Table 1. Data regarding quantitation of the targets in breath samples (H = healthy, CA = lung cancer, COPD = chronic obstructive pulmonary disease, AS = asthma, SD = standard deviation, nd = not detected, (−) = SD not calculated because analyte was quantified just in a single sample, nd = not detected).
Table 1. Data regarding quantitation of the targets in breath samples (H = healthy, CA = lung cancer, COPD = chronic obstructive pulmonary disease, AS = asthma, SD = standard deviation, nd = not detected, (−) = SD not calculated because analyte was quantified just in a single sample, nd = not detected).
AnalyteAverage Concentration (ppbv)Frequency in Samples (%)
H (SD)CA (SD)COPD (SD)AS (SD)TotalHCACOPDAS
2-Methylbutane1.52(1.32)3.73(6.15)1.63(0.50)1.72(1.85)37.525.050.033.350.0
Pentane1.66(0.67)2.21(1.09)1.87(0.53)2.11(2.26)51.845.062.541.750.0
Ethanol70.60(95.14)179.08(132.87)218.64(216.1)100.89(108.96)98.295.093.8100.0100.0
Isoprene32.85(34.01)34.19(30.08)34.61(20.9)48.70(16.95)100.095.0100.0100.0100.0
2-Propanol10.55(9.30)230.66(190.62)258.37(255.01)123.42(67.37)85.755.0100.0100.0100.0
2-Methylpentane1.24(0.30)3.44(2.41)2.61(2.07)4.59(5.57)55.425.075.075.050.0
3-Methylpentane0.24(0.12)0.93(0.72)1.27(0.49)1.07(1.25)35.710.068.833.325.0
1-Propanol14.59(14.63)34.10(37.73)28.15(38.54)9.94(5.77)100.095.0100.0100.0100.0
Methylcyclopentane1.80(0.53)2.49(1.11)2.20(0.47)2.20(0.27)87.575.093.883.3100.0
2-Butanone1.74(1.15)1.93(1.30)1.45(1.00)1.26(0.80)80.455.0100.083.387.5
Benzene1.13(0.83)0.29(−)0.57(−)0.60(0.09)16.125.06.38.325.0
Acetoin44.02(19.8)60.39(51.63)55.22(28.95)41.72(17.93)53.645.056.350.075.0
Toluene6.23(8.38)0.98(1.40)0.63(0.42)0.89(0.60)55.440.075.058.350.0
Ethylbenzene0.650(0.65)2.73(2.50)0.34(0.36)1.41(−)17.95.025.033.312.5
p-Xylene1.15(0.92)1.62(1.86)1.40(1.11)1.97(0.95)41.125.050.058.337.5
Styrene0.27(0.26)3.78(6.26)1.61(1.30)0.73(0.59)53.65.081.375.087.5
Decanend(−)nd(−)0.23(−)nd(−)1.80.00.08.30.0
6-Methyl-2-heptanone1.65(−)4.42(2.85)1.72(−)6.46(−)12.55.025.08.312.5
Isododecane0.69(0.49)1.59(1.57)0.98(0.48)0.52(0.35)76.845.093.883.3100.0
1,2,4-Trimethylbenzene0.83(0.61)2.55(1.94)2.60(2.46)1.42(0.76)82.150.093.8100.0100.0
(E)-Ocimene1.16(0.80)4.64(4.03)2.95(1.99)2.98(1.54)82.150.0100.091.7100.0
Limonene1.57(1.20)1.87(1.80)1.71(1.75)5.05(2.15)89.375.093.891.7100.0
m-Cymene0.61(0.21)0.41(0.35)0.38(0.23)0.32(0.21)46.410.062.550.087.5
Benzonitrile1.44(1.23)3.57(3.95)4.57(2.64)2.12(1.67)78.650.093.891.787.5
Phenolnd(−)52.78(47.13)75.02(72.16)nd(−)16.10.037.525.00.0
Undecane0.80(0.11)3.83(3.09)2.44(1.38)1.78(0.30)41.120.075.041.725.0
Dodecane5.18(0.72)10.58(8.66)9.51(7.18)6.27(3.19)73.245.087.591.787.5
Terpineol3.57(0.30)17.36(21.72)26.53(36.34)6.87(6.38)71.415.0100.0100.0100.0
Tridecane3.43(1.69)42.16(38.89)28.36(21.59)8.59(7.87)51.810.075.075.075.0
Table 2. RF model performance (AUC = area under the curve, CI = confidence interval).
Table 2. RF model performance (AUC = area under the curve, CI = confidence interval).
Statistics by ClassSensitivitySpecificityAUCBalanced Accuracy
Asthma75.0%100%0.87287.5%
Lung cancer93.8%87.5%0.95690.6%
COPD67.0%97.7%0.93582.2%
Healthy95.0%94.5%0.99494.7%
RF overall accuracy (95% CI)85.7% (73.7–93.6)
Table 3. Description of MLR model (AS = asthma, CA = lung cancer, COPD = chronic obstructive pulmonary disease, SE = standard error).
Table 3. Description of MLR model (AS = asthma, CA = lung cancer, COPD = chronic obstructive pulmonary disease, SE = standard error).
ConditionCoefficients
Intercept2-Propanol3-Methyl-pentane(E)-OcimeneLimonenem-CymeneBenzonitrileUndecaneTerpineolPhenol
AS−557.790.5617.0350.39100.40212.50−22.69−88.9432.33−12.95
(SE) (42.62)(0.37)(2.55)(2.48)(10.34)(11.52)(1.32)(3.64)(1.31)(5.31)
CA−127.120.61−129.3232.9232.3331.11−96.0726.5032.81−1.89
(SE)(29.25)(0.38)(17.43)(8.11)(9.90)(6.47)(16.91)(11.17)(8.40)(0.86)
COPD−23.490.49−269.07−64.2711.09−72.04−7.72−111.5237.911.73
(SE)(3.35)(0.83)(6.92)(8.57)(3.39)(0.66)(1.29)(4.22)(8.18)(0.65)
Table 4. MLR model performance (AUC = area under the curve, CI = confidence interval).
Table 4. MLR model performance (AUC = area under the curve, CI = confidence interval).
Statistics by ClassSensitivitySpecificityAUCBalanced Accuracy
Asthma68.0%100%0.70083.4%
Lung cancer83.4%93.4%0.95088.4%
COPD100%94.1%0.97197.1%
Healthy100%100%1.000100%
MLR overall accuracy (95% CI)91.0% (70.0–99.0)
Table 5. Main information regarding volunteers (SD = standard deviation, CA = lung cancer, COPD = chronic obstructive pulmonary disease, AS = asthma).
Table 5. Main information regarding volunteers (SD = standard deviation, CA = lung cancer, COPD = chronic obstructive pulmonary disease, AS = asthma).
GroupControlPositivep
n% n%
Total 20 36 0.367
GenderMale1365.0%2775.0%0.325
Female735.0%925.0%0.437
Age (SD) 41.2 (10.1)66.8 (8.22)0.078
Smoking statusActive smoker210.0%5 (2 COPD, 3 CA)13.9%0.287
Ex-smoker210.0%22 (12 CA, 10 COPD)61.1%0.083
Non-smoker1680.0%9 (1 CA, 8 AS)25.0%0.640
ConditionLung cancer1644.4%
COPD1233.3%
Asthma822.2%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Monedeiro, F.; Monedeiro-Milanowski, M.; Ratiu, I.-A.; Brożek, B.; Ligor, T.; Buszewski, B. Needle Trap Device-GC-MS for Characterization of Lung Diseases Based on Breath VOC Profiles. Molecules 2021, 26, 1789. https://doi.org/10.3390/molecules26061789

AMA Style

Monedeiro F, Monedeiro-Milanowski M, Ratiu I-A, Brożek B, Ligor T, Buszewski B. Needle Trap Device-GC-MS for Characterization of Lung Diseases Based on Breath VOC Profiles. Molecules. 2021; 26(6):1789. https://doi.org/10.3390/molecules26061789

Chicago/Turabian Style

Monedeiro, Fernanda, Maciej Monedeiro-Milanowski, Ileana-Andreea Ratiu, Beata Brożek, Tomasz Ligor, and Bogusław Buszewski. 2021. "Needle Trap Device-GC-MS for Characterization of Lung Diseases Based on Breath VOC Profiles" Molecules 26, no. 6: 1789. https://doi.org/10.3390/molecules26061789

APA Style

Monedeiro, F., Monedeiro-Milanowski, M., Ratiu, I. -A., Brożek, B., Ligor, T., & Buszewski, B. (2021). Needle Trap Device-GC-MS for Characterization of Lung Diseases Based on Breath VOC Profiles. Molecules, 26(6), 1789. https://doi.org/10.3390/molecules26061789

Article Metrics

Back to TopTop