Next Article in Journal
Identification of Tentative Traceability Markers with Direct Implications in Polyphenol Fingerprinting of Red Wines: Application of LC-MS and Chemometrics Methods
Next Article in Special Issue
Chemical Profiling and Characterization of Different Cultivars of Cannabis sativa L. Inflorescences by SPME-GC-MS and UPLC-MS
Previous Article in Journal
What You Extract Is What You Get: Different Methods of Protein Extraction from Hemp Seeds
Previous Article in Special Issue
Gas Chromatography-Mass Spectrometry Analysis of Compounds Emitted by Pepper Yellow Leaf Curl Virus-Infected Chili Plants: A Preliminary Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Method for the Measurement of Human Scent Samples Using Comprehensive Two-Dimensional Gas Chromatography with Mass Detection

Department of Analytical Chemistry, Faculty of Chemical Engineering, University of Chemistry and Technology, Prague, Technická 5, 166 28 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Separations 2021, 8(12), 232; https://doi.org/10.3390/separations8120232
Submission received: 31 October 2021 / Revised: 19 November 2021 / Accepted: 1 December 2021 / Published: 3 December 2021

Abstract

:
Every human body is a source of a unique scent, which can be used for medical or forensic purposes. Human skin scent is a complex mixture of more or less volatile compounds with different chemical and physical properties, which often differ significantly in their concentrations. The most efficient technique for separating such complex samples is comprehensive two-dimensional gas chromatography (GC × GC). This work aimed to find the optimal arrangement of a two-dimensional chromatographic system and define a suitable chromatographic method for non-targeted analysis of human scent samples. Four different chromatographic columns (non-polar Rxi-5MS and TG-5HT, medium polar Rxi-17Sil MS and Rtx-200MS) and their different configurations were tested. The best system was the 30 m primary column Rtx-200MS (with the 2 m pre-column Rtx-200MS) and the 1 m secondary column TG-5HT in a reverse configuration. This system achieved the highest theoretical and conditional peak capacities, optimal resolution, and the lowest number of coelutions.

1. Introduction

The human scent and its molecular profile are one of the major topics of current forensic interest. The possibility of its chemical analysis and identification of the individual based on scent composition would mean a significant breakthrough in forensic practice, providing more convincing evidence than specially trained dogs [1,2,3,4,5]. In addition to forensic practice, monitoring changes in the molecular profile of an individual would be of great importance for the evaluation of physical or mental conditions in biometric and medical applications [6,7,8].
Concerning human scent, a number of publications [9,10] focus mainly on volatile organic compounds (VOCs), which are largely contained in human skin scent and which can be found in all body fluids and excrement: in breath, saliva, blood, milk, urine, feces, and skin secretions (sweat and follicular fluids). In 2014, de Lacy Costello et al. [6] published a review with several hundred VOCs of the human scent. However, some authors point to the importance of low volatile compounds [11,12,13]. In general, the human skin scent contains several thousand different compounds, e.g., acids, alcohols, aldehydes, amides, amines, esters, hydrocarbons, ketones, heterocyclic compounds, fatty acids, and fatty acid esters. These compounds, in addition, can vary in concentrations by several orders of magnitude [12]. The qualitative and quantitative variability of the human scent is influenced by many circumstances, e.g., by the sampling condition and sampled body parts, age, sex, heredity, diet, drugs, mood and emotions, cosmetics, and environmental factors [14,15].
For the analysis of volatile compounds in real complex samples, such as human scent, conventional gas chromatography is often an insufficient separation tool. Some analytes cannot be separated during a single analysis due to their similar properties, and either overlapping or insufficiently distributed peaks appear in the chromatogram, causing problems for both identification and quantification.
The aim of the experiments described in this work is to find the optimal arrangement of the 2D chromatographic system and a suitable chromatographic method to achieve the best separation condition for non-targeted analysis of human scent samples.

2. Materials and Methods

2.1. Materials

Hexane (quality for GC-MS) and all standards used (see Table S1) were obtained from Sigma-Aldrich (St. Louis, MO, USA). Ethanol (for UV–VIS spectroscopy, min. 99.8%) was purchased from Penta (CZ). The chromatographic columns Rxi-5MS, Rxi-17Sil MS, Rtx-200 MS (Restek, Bellefonte, PA, USA) and TG-5HT (Thermo Scientific, Waltham, MA, USA) had the same column diameter of 0.25 mm and a stationary phase thickness of 0.25 μm. The carrier gas was helium 5.5 (purity ≥ 99.9995) from Linde (CZ).

2.2. Preparation of a Mixture of Standards

To develop and optimize the chromatographic measurement method, a model mixture of human skin scent containing 98 compounds (see Supplementary Materials, Table S1) with a concentration of 100 mg·L−1 of each component in hexane was created. A solution with a concentration of 25 mg·L−1 of each component in hexane was used for the work. The compounds were chosen to cover the widest possible range (chemical family, polarity, molecular weight) of compounds found in human skin scent samples.

2.3. Preparation of Real Scent Samples

The developed chromatographic methods were also verified by measuring real scent samples. The scent samples were taken on glass beads of the palms of volunteers. Details on the purification of the sampling material and the scent sampling are given in Appendix A. The samples were subsequently extracted by double extraction into 1 mL of ethanol, each extraction step lasting 20 min, of which 10 min extraction was supported by shaking (1000 rpm, IKAVIBRAX VXR basic, Staufen im Breisgau, Germany) and 10 min by placing the samples in an ultrasonic bath (Fisherbrand FB15059, Fisher Scientific, Loughborough, UK). The individual extracts were gradually evaporated to dryness in a 1.5 mL vial in a Genevac EZ-2 concentrator under reduced pressure (approx. 100 mbar). The residue was then redissolved in 100 μL of a 1:1 mixture of ethanol and hexane. The samples prepared in this way were stored in a refrigerator at 7 °C until the analysis.

2.4. Instrumentation

The samples were analyzed on a 7890B GC gas chromatograph (Agilent, Santa Clara, CA, USA) with a Pegasus® 4D-C time-of-life mass detector from LECO Corp. (St. Joseph, MI, USA). A MultiPurpose Sampler MPS Robotic from Gerstel (Mülheim, Germany) was used to dose samples to the column. The sample measurement and data evaluation were performed in the ChromaTOF® program (version 4.72.0.0, LECO Corp., St. Joseph, MI, USA).
During the development of the method, five main column configurations were tested, namely:
  • Classic arrangement: non-polar Rxi-5MS (2 + 30 m) + medium polar Rxi-17Sil MS (1 m);
  • Classic arrangement: non-polar Rxi-5MS (2 + 30 m) + medium polar Rtx-200MS (1 m);
  • Reverse arrangement: medium polar Rxi-17Sil MS (2 + 30 m) + non-polar TG-5HT (1 m);
  • Reverse arrangement: medium polar Rtx-200 MS (2 + 30 m) + non-polar TG-5HT (1 m);
  • Reverse arrangement: medium polar Rtx-200 MS (2 + 15 m) + non-polar TG-5HT (1 m),
where the first bracket shows the length of the pre-column and the primary column, the pre-column being the same as the primary column, and the second bracket shows the total length of the secondary column, of which 10 cm is in the modulator, 21 cm in the GC-MS transfer line and 69 cm in the secondary oven.
The injection of samples (1 µL) took place at a temperature of 280 °C in the splitless mode with a period of 2 min. The carrier gas flow was set at a constant 1.5 mL/min. The optimized temperature program was: 40 °C—2 min—5 °C/min—320 °C—10 min, the total time of one analysis was 68 min (4080 s). The secondary column was heated to a temperature 5 °C higher than the primary column and the modulator to a temperature 15 °C higher than the secondary column. In the case of a 15 m primary column instead of 30 m, the temperature program was shortened by 10 min (the final holding of the final temperature was eliminated), and the secondary column was heated 10 °C higher than the primary column.
Cryogenic modulation was performed on a secondary column using dried air at −80 °C. The modulation periods and the corresponding hot and cold pulse times were different for different column configurations; see Table 1. The warm pulse length corresponds to 0.3 times the modulation period length and the cold pulse length to 0.2 times the modulation period length.
The transfer line was set to 280 °C. The mass detector was in electron ionization (EI) mode with an ionizing electron energy of 70 eV, and the temperature of the ion source was 250 °C. Data collection took place in the Total Ion Current (TIC) mode; in the range 29–800 m/z, the data collection speed was 200 spectra/s. The acquisition delay was set to 500 s for a 30 m primary column, and in the case of a 15 m column it was shortened to 300 s. The voltage on the detector was set 200 V higher than tuning to increase the detection sensitivity.

2.5. Choosing the Best System Configuration

The different systems were compared with each other in terms of orthogonality, the used separation space, the maximum theoretical peak capacity and the condition peak capacity of the system, the achieved resolution, and the number of coelutions.

2.5.1. Spearman’s Correlation Coefficient

Spearman’s correlation coefficient is used to measure monotonic nonlinear correlations between retention times and can be used to measure the orthogonality of systems [16]. First, the order of retention times in the first (rit1) and second dimensions (rit2) is determined for compound i. The Spearman’s correlation coefficient (ρ) is then calculated according to Equation (1),
ρ = 1 6 i = 1 n d i 2 n ( n 2 1 )
where di = rit1rit2 and n is the number of compounds.

2.5.2. Bin Counting Method

The bin counting method [17,18] serves both as an indicator of orthogonality and an indicator of the coverage of the separation space. In the bin counting method, retention times in both dimensions must first be normalized according to Equation (2),
t r , n o r m ( i ) = t r ( i ) t r ,   f i r s t t r ,   l a s t t r ,   f i r s t
where 1tr,norm(i) and 2tr,norm(i) are normalized retention times in the first and the second dimension, tr(i) is the retention time of compound i, and tr,first and tr,last are the retention times of the first and last observed peaks in the separation space, respectively.
Then, the separation space is divided into bins, whose number is the same or very close to the number of compounds in the sample. Subsequently, the number of bins containing a peak is counted and compared with the total number of bins according to Equation (3),
O = b i n s P m a x 0.63 P m a x
where Pmax is the total number of bins. Occupied space O takes values from 0 (completely non-orthogonal system) to 1 (completely orthogonal system). The coefficient of 0.63 is based on the observation that systems close to orthogonality have a ratio of filled and total bins (so-called fcoverage) close to 0.63.
For reverse systems, the separation space was divided into three parts, which corresponded to the individual modulation periods. The total value O of the given reverse system was calculated as a weighted average of partial O.

2.5.3. Peak Capacities

The theoretical peak capacities [19] of the system in both dimensions (1nc and 2nc) were calculated according to Equation (4).
n c = t r , l a s t t r ,   f i r s t w b
where tr is the retention time of the last of first observed peaks in the separation space, and wb is the measured average 4σ peak width. The maximum theoretical GC × GC capacity of the 2Dnc system was then determined as the product of the peak capacities in the first (1nc) and the second (2nc) dimensions (Equation (5)).
n 2 D c = n 1 c · n 2 c
However, orthogonality and surface coverage have a huge impact on the calculations of the conditional (practical) peak capacity 2Dnc,cond of a comprehensive 2D system [19]. It is given by Equation (6),
n 2 D c , c o n d = n 1 c · n 2 c · f c o v e r a g e · 1 β
where fcoverage is the ratio of filled and total bins (see Section 2.5.2), and β is an undersampling of the first dimension, and it is given by Equation (7),
β = 1 + 0.21 ( t s σ 1 ) 2
where ts is the first-dimension sampling time, and 1σ is the peak standard deviation in the first dimension.
For reverse systems, the separation space was divided into three parts, which corresponded to the individual modulation periods. The total value of theoretical peak capacity and the conditional peak capacity of the given reverse system were calculated as a sum of partial capacitances.

2.5.4. Resolutions

With knowledge of the maximum theoretical peak capacities (see Equation (4)), the average resolutions (Rs,av) in both dimensions can be calculated [20]. It is given by Equation (8),
R s ,   a v = n c n 1
where n is the number of compounds.
The resolutions of the two adjacent peaks in the individual dimensions (1Rs and 2Rs) were calculated from the retention times (tr1 and tr2) and the width of the peaks at the base (wb1 and wb2) according to Equation (9) [21],
R s = t r 2 t r 1 ( w b 1 + w b 2 ) / 2
and the total resolution in a 2D system of 2DRs was then calculated as the square root of the sum of the squares of the resolutions obtained in the first and second dimensions; see Equation (10) [21].
R 2 D s = R 1 s 2 + R 2 s 2

3. Results

3.1. Classic Arrangement: Non-Polar Rxi-5MS (2 + 30 m) + Medium Polar Rxi-17Sil MS (1 m)

The first column arrangement tested was the 30 m Rxi-5ms column with the 2 m Rxi-5ms precolumn in the first dimension and the 1.5 m Rxi-17Sil MS column in the second dimension. The Rxi-5ms column with a (5%-phenyl)-methylpolysiloxane phase has a low polarity, which is suitable for a wide range of analyses. The second Rxi-17Sil MS column with a (50%-phenyl)-methylpolysiloxane phase is characterized by medium polarity and higher selectivity to aromatic substances and is usually recommended for environmental analyses. After the initial analyses, however, it was observed that when using a column length of 1.5 m in the second dimension, the retention times in the second dimension are too high related to the modulation period and a large area of unused separation space was formed at the bottom of the chromatogram. For this reason, the secondary column was shortened to 1 m, and this column length was used as the standard for the other column arrangements.
The ramp of the temperature program was also tested under the above conditions (see Section 2.4). At a temperature rate of 3 °C/min, the temperature program was very long (105 min), at 8 °C/min, the temperature program was shortened to 47 min, however, the resolution of the peaks in the first dimension deteriorated significantly. The rate of 5 °C/min was chosen as the most suitable, as a compromise between the total length of the program and the achieved resolution in the first dimension.
Figure 1 shows the chromatogram obtained under the optimized conditions (see Section 2.4). The classification zones of homologues of various compounds are highlighted. n-alkanes are marked with black contour, with which amines partially coelute (lower dark yellow zone). For instance, heptadecane with N,N-dimethyltetradecanamide achieved an overall resolution of 2DRs = 0.65. Acids and esters (left red zone), aldehydes and ketones (yellow zone), alcohols (left white zone), and chloroalkanes (green zone) showed similar retention on the secondary column. For example, dodecan-1-ol with undecanoic acid achieved an overall resolution of 2DRs = 0.52 and dodecan-1-ol with 1-chlordodecane 2DRs = 0.56. Ethers (right white zone) and amides (upper dark yellow zone) are then clearly separated from these compounds. The right red zone then represents the area of higher fatty acid esters. Unlabeled compounds with higher retention on the secondary column represent compounds classified as “others”, which are often part of cosmetic products. The list of standards including retention indices calculated from absolute retention times according to [22] is given in Supplementary Material, Table S1.
Unfortunately, it was not possible to distinguish structural isomers such as hexadecyl ester of tetradecanoic acid and dodecyl ester of octadecanoic acid in the model mixture (2DRs = 0, both compounds have the general formula C30H60O2 and Mr = 452.80, their mixed mass spectrum see Figure 2). The total number of analyte coelutions (with 2DRs < 1.5) in the mixture was 38.
As described above, the Rxi-17Sil MS column is selective for aromatics, resulting in high retention of aromatic compounds in the second dimension. This was also shown in the subsequent analysis of the real sample, when the modulation period of 4 s was too short, and the so-called wrap-around effect occurred; see Figure 3.

3.2. Classic Arrangement: Non-Polar Rxi-5MS (2 + 30 m) + Medium Polar Rtx-200MS (1 m)

In another tested arrangement, the secondary column was replaced with the Rtx-200MS column with a (trifluoropropyl)-methylpolysiloxane phase, which, like Rxi-17Sil MS, belongs to the group of the mid-polar columns; however, the polarity is slightly higher. The great advantage of the Rtx-200MS column is its unique selectivity for electronegative substances, which changes the elution order and can ensure a resolution that cannot be achieved on the phenyl, cyano or polyethylene glycol phases.
The chromatogram of the model mixture measured on the Rxi-5ms and the Rtx-200MS columns under optimized conditions (see Section 2.4) is then shown in Figure 4. Several differences can be observed compared to the chromatogram in Figure 1. First, there was a clear separation of aldehydes and ketones in the second dimension (yellow classification zone) from acids and esters (first red zone). Second, there was a significant reduction in aromatic retention in the second dimension (see caffeine, which had 2tr = 3.18 s for the Rxi-17Sil MS column, 2tr = 1.66 s for the Rtx-200 MS column). However, the worst--separated places remained poorly differentiated. Amines (lower dark yellow zone) still coeluted with n-alkanes (black zone), e.g., the previously discussed heptadecane with N,N-dimethyltetradecanamine reached 2DRs = 0.57. Alcohols (left white zone) then still showed similar retention as acids (left red zone) and chloroalkanes (green zone). For example, dodecan-1-ol with undecanoic acid achieved an overall resolution of 2DRs = 0.65 and dodecan-1-ol with 1-chlordodecane 2DRs = 0.61. The structural isomers in the area of higher fatty acid esters (right red zone) were not distinguished at all (2DRs = 0). The total number of analyte coelutions (with 2DRs < 1.5) in the mixture was 34.
Due to the lower retention of aromatic compounds, in the case of a real sample, there was no longer a wrap-around effect during the modulation period of 4 s. See Figure 5 for the chromatogram of the real sample.

3.3. Reverse Arrangement: Medium Polar Rxi-17Sil MS (2 + 30 m) + Non-Polar TG-5HT (1 m)

Due to the failure of achieving satisfactory separation in the classical arrangement of the columns, the order of the tested columns was reversed. The first reverse system consisted of the 30 m primary column Rxi-17Sil MS with the 2 m precolumn, also Rxi-17Sil MS. A 1 m long TG-5HT was chosen as the secondary column, which is formed as the Rxi-5ms by a non-polar (5%-phenyl)-methylpolysiloxane phase; however, it is designed up to temperatures of 400 °C and the manufacturer guarantees up to 40% longer life compared to similar columns.
Figure 6 shows the measured chromatogram obtained under optimal conditions (see Section 2.4), from which it is clear, that the elution of different types of substances was reversed compared to the classical systems. The highest retention on the secondary column had n-alkanes (black classification zone), from which amines (upper dark yellow zone) were visibly separated, e.g., the mentioned problem site of heptadecane and N,N-dimethyltetradecanamine reached an overall resolution of 2DRs = 4.15. Acids and esters (left red zone) and aldehydes and ketones (yellow zone) were separated from each other within the first dimension, but not within the second dimension. Alcohols (left white zone), such as dodecan-1-ol, were further separated from the acids in the first dimension, and undecanoic acid reached a resolution of 2DRs = 1.15. Previously mentioned problematic places, namely, alcohols (white zone) and chloroalkanes (green zone), were successfully separated within the first and second dimensions (the distinction between dodecan-1-ol and 1-chlordodecane reached 2DRs = 3.51). The structural isomers in the area of higher fatty acid esters (right red zone) again remained indistinguishable (2DRs = 0). However, the distribution of amines and chloroalkanes deteriorated, e.g., for N,N-dimethyldodecanamine with 1-chlordodecane the resolution was 2DRs = 0.71 (in both classical column arrangements the resolutions were 2DRs~3.7). The total number of analyte coelutions (with 2DRs < 1.5) in the mixture was 24. The chromatogram of a real scent sample is shown in Figure 7.

3.4. Reverse Arrangement: Medium Polar Rtx-200 MS (2 + 30 m) + Non-Polar TG-5HT (1 m)

The second reverse system consisted of the 30 m primary column Rtx-200 MS with the 2 m precolumn, also Rtx-200 MS. The 1 m long TG-5HT remained as a secondary column. Figure 8 shows that all studied compounds’ groups have already been separated. In the second dimension, there was a more significant separation of amines (upper dark yellow classification zone) from chloroalkanes (green zone) within the second dimension, alcohols (left white zone) were more significantly separated from the area of acids and esters (left red zone), from whose zone also separated aldehydes and ketones (yellow zone). At the same time, however, there was a change in the retention of aromatic or polyunsaturated hydrocarbons (see squalene, cholesterol, caffeine). In this arrangement, squalene eluted near octacosane with a resolution of 2DRs = 1.34 and cholesterol near tetradecyl ester of hexadecanoic acid (2DRs = 2.15). It was again not possible to distinguish the structural isomers in the region of higher fatty acid esters (second red bubble, 2DRs = 0). The total number of analyte coelutions (with 2DRs < 1.5) in the mixture was 7. The chromatogram of a real scent sample is shown in Figure 9.

3.5. Reverse Arrangement 3: Medium Polar Rtx-200 MS (2 + 15 m) + Non-Polar TG-5HT (1 m)

To shorten the analysis time, an alternative to a 15 m primary column instead of a 30 m column was tested. While maintaining the original temperature gradient, the total analysis time was reduced by 10 min (for more details on setting up the analysis, see Section 2.4). In this configuration, however, the resolution of the compounds within the first dimension deteriorated (compared to the 30 m system, see Table 2) and the shape of the peaks deteriorated, in particular, their profile (see Figure 10). The total number of analyte coelutions (with 2DRs < 1.5) in the mixture was 21.

4. Discussion

Five different configurations of the 2D chromatographic system were tested to find its optimal arrangement and at the same time with the intention of finding a suitable chromatographic method for non-targeted analysis of human odor samples. In testing chromatographic methods, mainly the temperature ramp was optimized. The rate of 5 °C/min was chosen as the most suitable, as a compromise between the total length of the program and the achieved resolution in the first dimension. Individual chromatographic systems were evaluated in terms of orthogonality, use of separation space, peak capacities, resolutions, and finally the number of coelutions (see Table 2). The first of the tested systems in the classical arrangement (Rxi-5 ms/Rxi-17Sil MS) reached a high theoretical peak capacity; however, the use of separation space was the worst compared to the other tested systems. Thanks to this, the conditional peak capacity was considerably lower, and the number of coelutions of analytes of the model mixture was the highest (38). In the case of the analysis of a real scent sample, which contained, e.g., many aromatic heterocyclic compounds with high retention, the so-called wrap-around effect. The second system in the classical arrangement (Rxi-5ms/Rtx-200 MS) showed the highest orthogonality according to Spearman’s correlation coefficient, and the use of the separation space was considerably higher than in the first case. The theoretical and conditional peak capacities were generally lower than in the first case, as the peak capacity in the second dimension decreased due to the lower retention of aromatic compounds (e.g., caffeine). However, the average resolution of the analytes improved slightly, but multiple coelutions still occurred (34 in total). Due to the lower retention of aromatic compounds, there was no longer a wrap-around effect in the analysis of a real scent sample. Reverse systems (except Rxi-17Sil MS/TG-5HT system) according to Spearman’s correlation coefficient, showed lower orthogonality than systems in the classical arrangement; nevertheless, the use of separation space was considerably higher, thanks to which the maximum theoretical and conditional peak capacities were higher than in the case of classical systems. The highest conditional peak capacity was then achieved with the Rtx-200 MS/TG-5HT system. The average 2D resolution was lower than in the case of classical systems; however, the number of coelutions (7 in total) was the lowest in this case, mainly due to the much higher resolution in the second dimension than was the case in previous systems. This was also shown in the case of analyses of real samples when, in the case of Rxi-17 Sil MS/TG-5HT columns, the separation of the above-mentioned aromatic heterocyclic compounds was insufficient. In the case of the Rtx 200 MS/TG-5HT system, even these heterocyclic compounds were sufficiently separated.
The reverse system (Rtx-200 MS/TG-5HT) was evaluated as the optimal system for human scent analysis and an alternative of a 15 m primary column instead of a 30 m column was tested on this arrangement to shorten the analysis time. However, when using the 15 m column, the theoretical and conditional peak capacities decreased considerably, especially in the first dimension, which is due to the deterioration of the peak shape and, above all, peaks expansion.
The deterioration of the peak shape is mainly related to the long modulation periods (15, 18, and 20 s), which were needed to avoid the wrap-around effect. The high retention in the second dimension can be explained by the fact that, at the same temperature program, the analytes are eluted from the shorter primary column at a lower temperature, i.e., the separation on the secondary column also occurs at a lower temperature, which is associated with prolongation of retention. With a higher setting of the temperature offset between the primary and secondary ovens, there was a significant loss of resolution in the second dimension. Furthermore, it must be said that long modulation periods also correlate with long periods of cold and warm pulses (see Table 1 in Section 2.4). The problem is that short columns generate relatively narrow primary peaks, which require rapid and sharp reinjections into the secondary system to maintain properly performed modulation and cryofocusing, which was not followed in this case. The solution could be a further optimization of the modulation process, a change of the length of the secondary column, or change of the diameter and thickness of the stationary phase on the primary and secondary columns, respectively.
None of the systems was able to distinguish the structural isomers of higher fatty acid esters, which were represented in the model mixture by hexadecyl ester of tetradecanoic acid and dodecyl ester of octadecanoic acid (both compounds have the general formula C30H60O2 and Mr = 452.80). Due to their structure and the long aliphatic chains, it would probably be necessary to use a very non-polar column, such as Rxi-1ms, to distinguish between them. In addition to changing the stationary phase, it would also probably be necessary to increase the separation efficiency, either by increasing the length of the column, which would, however, lead to very long temperature programs, or by reducing the internal diameter of the column and the thickness of the stationary phase, which in turn could cause insufficient column capacity. However, although the isomers could not be resolved chromatographically, they could be distinguished based on the characteristic masses in the mass spectrum. These esters have a typical m/z, which corresponds to the loss of the alkyl chain of the alcohol associated with the migration of two hydrogens.

5. Conclusions

Using a model mixture containing 98 compounds commonly found in human skin scent, five different column configurations were tested for GC × GC analysis. The best results were obtained for the reverse system, consisting of the 2 m pre-column Rtx-200 MS, the 30 m primary column Rtx-200 MS and the 1 m secondary column TG-5HT. On this reverse system, using the chromatographic method with a temperature gradient of 5 °C/min, the highest utilization of the separation space and the highest theoretical and conditional peak capacities of the system were achieved, and the smallest number of coelutions occurred. Therefore, this system layout will continue to be used in analyses with real scent samples.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/separations8120232/s1, Table S1: The list of the standards in a model mixture of human skin scent and measured retention indices on systems with a 2 m pre-column, a 30 m primary column, and a 1 m secondary column.

Author Contributions

Conceptualization, P.P. and N.L.; methodology, P.P. and N.L.; formal analysis, P.P.; resources, Š.U.; data curation, P.P.; writing—original draft preparation, P.P.; writing—review and editing, N.L. and Š.U.; supervision, Š.U.; project administration, Š.U.; funding acquisition, Š.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of the Interior of the Czech Republic projects No.VI20172020075, VH20182021030, and VJ01010123.

Institutional Review Board Statement

This study was approved by the Institutional Review Board of University of Chemistry and Technology, Prague (approval number EK/8/2020) and complies with the Declaration of Helsinki for Medical Research involving Human Subjects.

Informed Consent Statement

The recruited volunteers were verbally instructed regarding the sampling method and analysis. They signed an agreement (written consent) regarding using his/her scent samples in this experiment and publishing the results obtained. The samples of the volunteers were used for this experiment only, and then these samples were destroyed.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Cleaning of the sampling material: Sodium–calcium glass beads with a diameter of 3.6 mm (Glass Sphere, CZ) were used for the scent collection. The glass beads were cleaned according to the methodology described here [12]. First, a chromium-sulfur mixture was prepared by gradually mixing 30 g of potassium dichromate in 0.5 L of concentrated sulfuric acid, in which the glass beads were placed for at least 24 h. The excess chromium–sulfur mixture was poured through a watch glass back into the storage bottle. The remainder of the chromium–sulfur mixture was then reduced with dilute hydrogen peroxide (approximately 5% solution). This was followed by three washes with super-pure water and solvents in the order of ethanol, hexane, and again ethanol. Then, the purified glass beads were placed in an oven (Memmert UF 30 Plus, Schwabach, Germany) and dried at 200 °C for three hours. The cleaned glass beads were stored in a desiccator until use, but for a maximum time of one week.
Human skin scent sampling: Before the scent sampling, volunteers were asked to refrain from using any cosmetic products (perfumes, antiperspirants, creams, makeup, etc.) 24 h prior to the sampling. The scent was sampled of the volunteers’ palms in the manner described herein [12]. First, the volunteers thoroughly lathered their hands and forearms with unscented Amadeus Neutral soap (Cormen, Bystřice nad Perštejnem, CZ), which they then rinsed thoroughly under warm running water. After drying their hands in the open air, the volunteers rubbed their palms together for 5 min. Then, cleaned glass beads (70 pcs, which corresponds to about 4 mL) were poured into their palms. The actual collection lasted 10 min and took place in such a way that the volunteers held and rubbed the cleaned glass beads in their hands. The sampled beads were then carefully transferred to a 10 mL headspace vial, in which a subsequent extraction took place (see Section 2.3).

References

  1. Cuzuel, V.; Cognon, G.; Rivals, I.; Sauleau, C.; Heulard, F.; Thiebaut, D.; Vial, J. Origin, Analytical Characterization, and Use of Human Odor in Forensics. J. Forensic Sci. 2017, 62, 330–350. [Google Scholar] [CrossRef] [PubMed]
  2. Cuzuel, V.; Sizun, A.; Cognon, G.; Rivals, I.; Heulard, F.; Thiebaut, D.; Vial, J. Human odor and forensics. Optimization of a comprehensive two-dimensional gas chromatography method based on orthogonality: How not to choose between criteria. J. Chromatogr. A 2018, 1536, 58–66. [Google Scholar] [CrossRef]
  3. Woidtke, L.; Dreßler, J.; Babian, C. Individual human scent as a forensic identifier using mantrailing. Forensic Sci. Int. 2018, 282, 111–121. [Google Scholar] [CrossRef]
  4. Ferry, B.; Ensminger, J.J.; Schoon, A.; Bobrovskij, Z.; Cant, D.; Gawkowski, M.; Hormila, I.; Kos, P.; Less, F.; Rodionova, E.; et al. Scent lineups compared across eleven countries: Looking for the future of a controversial forensic technique. Forensic Sci. Int. 2019, 302, 109895. [Google Scholar] [CrossRef] [PubMed]
  5. Cuzuel, V.; Leconte, R.; Cognon, G.; Thiebaut, D.; Vial, J.; Sauleau, C.; Rivals, I. Human odor and forensics: Towards Bayesian suspect identification using GCxGC-MS characterization of hand odor. J. Chromatogr. B 2018, 1092, 379–385. [Google Scholar] [CrossRef]
  6. de Lacy Costello, B.; Amann, A.; Al-Kateb, H.; Flynn, C.; Filipiak, W.; Khalid, T.; Osborne, D.; Ratcliffe, N.M. A review of the volatiles from the healthy human body. J. Breath Res. 2014, 8, 014001. [Google Scholar] [CrossRef]
  7. Buljubasic, F.; Buchbauer, G. The scent of human diseases: A review on specific volatile organic compounds as diagnostic biomarkers. Flavour Fragr. J. 2015, 30, 5–25. [Google Scholar] [CrossRef]
  8. de la Mata, A.P.; McQueen, R.H.; Nam, S.L.; Harynuk, J.J. Comprehensive two-dimensional gas chromatographic profiling and chemometric interpretation of the volatile profiles of sweat in knit fabrics. Anal. Bioanal. Chem. 2017, 409, 1905–1913. [Google Scholar] [CrossRef] [PubMed]
  9. Pandey, S.K.; Kim, K.-H. Human body-odor components and their determination. TrAC Trends Anal. Chem. 2011, 30, 784–796. [Google Scholar] [CrossRef]
  10. Prada, P.A.; Curran, A.M.; Furton, K.G. Comparison of extraction methods for the removal of volatile organic compounds (VOCs) present in sorbents used for human scent evidence collection. Anal. Methods 2010, 2, 470–478. [Google Scholar] [CrossRef]
  11. Doležal, P.; Kyjaková, P.; Valterová, I.; Urban, Š. Qualitative analyses of less-volatile organic molecules from female skin scents by comprehensive two dimensional gas chromatography–time of flight mass spectrometry. J. Chromatogr. A 2017, 1505, 77–86. [Google Scholar] [CrossRef] [PubMed]
  12. Pojmanová, P.; Ladislavová, N.; Škeříková, V.; Kania, P.; Urban, Š. Human scent samples for chemical analysis. Chem. Pap. 2020, 74, 1383–1393. [Google Scholar] [CrossRef] [Green Version]
  13. Doležal, P.; Furton, K.G.; Lněničková, J.; Kyjaková, P.; Škeříková, V.; Valterová, I.; Pinc, L.; Urban, Š. Multiplicity of human scent signature. Egypt. J. Forensic Sci. 2019, 9, 7. [Google Scholar] [CrossRef] [Green Version]
  14. Jha, S.K.; Hayashi, K. Molecular structural discrimination of chemical compounds in body odor using their GC–MS chromatogram and clustering methods. Int. J. Mass Spectrom. 2017, 423, 1–14. [Google Scholar] [CrossRef]
  15. Colon-Crespo, L.J.; Herrera-Hernandez, D.; Holness, H.; Furton, K.G. Determination of VOC marker combinations for the classification of individuals by gender and race/ethnicity. Forensic Sci. Int. 2017, 270, 193–199. [Google Scholar] [CrossRef] [PubMed]
  16. Bakain, R.; Rivals, I.; Sassiat, P.; Thiébaut, D.; Hennion, M.-C.; Euvrard, G.; Vial, J. Comparison of different statistical approaches to evaluate the orthogonality of chromatographic separations: Application to reverse phase systems. J. Chromatogr. A 2011, 1218, 2963–2975. [Google Scholar] [CrossRef]
  17. Gilar, M.; Olivova, P.; Daly, A.E.; Gebler, J.C. Orthogonality of separation in two-dimensional liquid chromatography. Anal. Chem. 2005, 77, 6426–6434. [Google Scholar] [CrossRef]
  18. Davis, J.M. Dependence of Effective Peak Capacity in Comprehensive Two-Dimensional Separations on the Distribution of Peak Capacity between the Two Dimensions. Anal. Chem. 2008, 80, 8122–8134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Stoll, D.R.; Wang, X.; Carr, P.W. Comparison of the practical resolving power of one- and two-dimensional high-performance liquid chromatography analysis of metabolomic samples. Anal. Chem. 2008, 80, 268–278. [Google Scholar] [CrossRef]
  20. Wang, X.; Stoll, D.R.; Schellinger, A.P.; Carr, P.W. Peak Capacity Optimization of Peptide Separations in Reversed-Phase Gradient Elution Chromatography:  Fixed Column Format. Anal. Chem. 2006, 78, 3406–3416. [Google Scholar] [CrossRef] [Green Version]
  21. Sanz, J. Chapter 3 Theoretical Considerations. In Comprehensive Analytical Chemistry; Ramos, L., Ed.; Elsevier: Amsterdam, The Netherlands, 2009; Volume 55, pp. 49–76. [Google Scholar]
  22. van Den Dool, H.; Kratz, P.D. A generalization of the retention index system including linear temperature programmed gas—liquid partition chromatography. J. Chromatogr. A 1963, 11, 463–471. [Google Scholar] [CrossRef]
Figure 1. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-17Sil MS (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Figure 1. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-17Sil MS (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Separations 08 00232 g001
Figure 2. The mixed mass spectrum of hexadecyl ester of tetradecanoic acid and dodecyl ester of octadecanoic acid. The arrows indicate their typical m/z, which correspond to the loss of the alkyl chain of the alcohol associated with the migration of two hydrogens.
Figure 2. The mixed mass spectrum of hexadecyl ester of tetradecanoic acid and dodecyl ester of octadecanoic acid. The arrows indicate their typical m/z, which correspond to the loss of the alkyl chain of the alcohol associated with the migration of two hydrogens.
Separations 08 00232 g002
Figure 3. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-17Sil MS (1 m) columns.
Figure 3. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-17Sil MS (1 m) columns.
Separations 08 00232 g003
Figure 4. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-200MS (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Figure 4. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rxi-200MS (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Separations 08 00232 g004
Figure 5. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rtx-200MS (1 m) columns.
Figure 5. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-5ms (2 + 30 m) and Rtx-200MS (1 m) columns.
Separations 08 00232 g005
Figure 6. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-17Sil MS (2 + 30 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Figure 6. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-17Sil MS (2 + 30 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Separations 08 00232 g006
Figure 7. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-17Sil MS (2 + 30 m) and TG-5HT (1 m) columns.
Figure 7. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-17Sil MS (2 + 30 m) and TG-5HT (1 m) columns.
Separations 08 00232 g007
Figure 8. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 30 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Figure 8. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 30 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Separations 08 00232 g008
Figure 9. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 30 m) and TG-5HT (1 m) columns.
Figure 9. The TIC 2D chromatogram of a real scent sample measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 30 m) and TG-5HT (1 m) columns.
Separations 08 00232 g009
Figure 10. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 15 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Figure 10. The TIC 2D chromatogram of the model mixture (composition see Table 1) measured under the conditions described in Section 2.4 on Rxi-200 MS (2 + 15 m) and TG-5HT (1 m) columns. The marked classification zones are the following: black (n-alkanes), lower dark yellow (amines), upper dark yellow (amides), green (chloroalkanes), left white (alcohols), right white (ethers), yellow (aldehydes and ketones), left red (acids and esters), right red (esters of higher fatty acids), and unlabeled (others).
Separations 08 00232 g010
Table 1. Modulation periods and durations of hot and cold pulses for different column configurations.
Table 1. Modulation periods and durations of hot and cold pulses for different column configurations.
Column ConfigurationTime [s]Modulation Period [s]Hot Pulse [s]Cold Pulse [s]
Rxi-5 ms/Rxi-17Sil MS500–408041.20.8
Rxi-5 ms/Rtx-200 MS500–408041.20.8
Rxi-17Sil MS/TG-5HT500–259241.20.8
2592–319861.81.2
3198–408082.41.6
Rtx-200 MS (30 m)/TG-5HT500–170461.81.2
1704–259282.41.6
2592–4080103.31.7
Rtx-200 MS (15 m)/TG-5HT300–1995154.53.0
1995–2499185.43.6
2499–3480206.04.0
Table 2. The evaluation of orthogonality, use of separation space, peak capacities, average resolutions of individual systems, and the number of coelutions.
Table 2. The evaluation of orthogonality, use of separation space, peak capacities, average resolutions of individual systems, and the number of coelutions.
IndicatorRxi-5 ms/Rxi-17Sil MSRxi-5 ms/Rtx-200 MSRxi-17Sil MS/TG-5HTRtx-200 MS (30 m)/TG-5HTRtx-200 MS (15 m)/TG-5HT
ρ0.535−0.1940.3560.6180.648
O0.5080.7300.7460.8170.885
2Dnc2236888189631311474
2Dnc,cond8054259141329589
1Rs,av2.6242.7032.5922.4171.334
2Rs,av0.0870.0340.2730.4320.324
2DRs,av2.6252.7032.6072.4551.373
2DRs < 1.5383424721
ρ is Spearman’s correlation coefficient, O is the occupied space, 2Dnc is the theoretical maximal peak capacity, 2Dnc,cond is the conditional peak capacity, Rs,av are resolutions in the first, second dimension or total 2D resolution, and 2DRs < 1.5 is the number of analyte pairs with 2DRs < 1.5.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pojmanová, P.; Ladislavová, N.; Urban, Š. Development of a Method for the Measurement of Human Scent Samples Using Comprehensive Two-Dimensional Gas Chromatography with Mass Detection. Separations 2021, 8, 232. https://doi.org/10.3390/separations8120232

AMA Style

Pojmanová P, Ladislavová N, Urban Š. Development of a Method for the Measurement of Human Scent Samples Using Comprehensive Two-Dimensional Gas Chromatography with Mass Detection. Separations. 2021; 8(12):232. https://doi.org/10.3390/separations8120232

Chicago/Turabian Style

Pojmanová, Petra, Nikola Ladislavová, and Štěpán Urban. 2021. "Development of a Method for the Measurement of Human Scent Samples Using Comprehensive Two-Dimensional Gas Chromatography with Mass Detection" Separations 8, no. 12: 232. https://doi.org/10.3390/separations8120232

APA Style

Pojmanová, P., Ladislavová, N., & Urban, Š. (2021). Development of a Method for the Measurement of Human Scent Samples Using Comprehensive Two-Dimensional Gas Chromatography with Mass Detection. Separations, 8(12), 232. https://doi.org/10.3390/separations8120232

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop