Skip to main content
  • Research article
  • Open access
  • Published:

Effects of smoking and smoking cessation on human serum metabolite profile: results from the KORA cohort study



Metabolomics helps to identify links between environmental exposures and intermediate biomarkers of disturbed pathways. We previously reported variations in phosphatidylcholines in male smokers compared with non-smokers in a cross-sectional pilot study with a small sample size, but knowledge of the reversibility of smoking effects on metabolite profiles is limited. Here, we extend our metabolomics study with a large prospective study including female smokers and quitters.


Using targeted metabolomics approach, we quantified 140 metabolite concentrations for 1,241 fasting serum samples in the population-based Cooperative Health Research in the Region of Augsburg (KORA) human cohort at two time points: baseline survey conducted between 1999 and 2001 and follow-up after seven years. Metabolite profiles were compared among groups of current smokers, former smokers and never smokers, and were further assessed for their reversibility after smoking cessation. Changes in metabolite concentrations from baseline to the follow-up were investigated in a longitudinal analysis comparing current smokers, never smokers and smoking quitters, who were current smokers at baseline but former smokers by the time of follow-up. In addition, we constructed protein-metabolite networks with smoking-related genes and metabolites.


We identified 21 smoking-related metabolites in the baseline investigation (18 in men and six in women, with three overlaps) enriched in amino acid and lipid pathways, which were significantly different between current smokers and never smokers. Moreover, 19 out of the 21 metabolites were found to be reversible in former smokers. In the follow-up study, 13 reversible metabolites in men were measured, of which 10 were confirmed to be reversible in male quitters. Protein-metabolite networks are proposed to explain the consistent reversibility of smoking effects on metabolites.


We showed that smoking-related changes in human serum metabolites are reversible after smoking cessation, consistent with the known cardiovascular risk reduction. The metabolites identified may serve as potential biomarkers to evaluate the status of smoking cessation and characterize smoking-related diseases.

Peer Review reports


Smoking is responsible for 90% of all lung cancers, accounts for 25% of cancer deaths worldwide [13] and is a significant risk factor for cardiovascular disease (CVD) [47]. The benefits of smoking cessation are remarkable. Risk of CVD is reduced in former smokers (FS) compared with current smokers (CS) [810]; mortality and future cardiac events both decline in FS [11, 12]. Nevertheless, for cancers, especially for adenocarcinoma, the risk remains high in FS compared with never smokers (NS) [13, 14]. Studies have made attempts to find the molecular basis for the influence of smoking and smoking cessation on cardiovascular risks. For instance, smoking is associated with the increase of several CVD-related inflammatory markers, for example, c-reactive protein and fibrinogen [1517], and smoking cessation could largely reduce the level of these markers [18]. However, there is also evidence that other molecular changes associated with smoking are permanent, for example, loss of heterozygosity and hypermethylation in the promoter regions of cancer-related genes [1923].

The metabolomics approach provides a functional readout of activities located downstream of the gene expression level that are more closely related to the physiological status [24] and, thus, may be particularly useful for the study of environmental influences, namely the 'exposome' [25]. Studying a strong environmental factor, for example a lifestyle-related exposure to smoking, may be considered a very powerful approach for understanding the links between environmental exposure and the metabolome. In human lung epithelial cells, it has been shown that metabolite concentration changes in various pathways, for example, the urea cycle and polyamine metabolism and lipid metabolism under smoke exposure [26]. In a pilot study with 283 male participants from the Cooperative Research in the Region of Augsburg (KORA) F3 in Germany, we have shown that levels of diacyl-phosphatidylcholines (PCs) were higher in 28 CS compared with 101 NS, except for acyl-alkyl-PCs [1]. The reduced ratios of acyl-alkyl-to diacyl-PCs in CS may be regulated by the enzyme alkyl-dihydroxyacetone phosphate in both ether lipid and glycerophospholipid pathways [1]. However, little has been reported about the reversibility of the metabolite profile upon smoking cessation, which is important for comprehensive understanding of smoking effects. It is also known that metabolite profile is different between men and women [25], but whether lifestyle factors such as smoking may induce different metabolite patterns in men and women is still unknown.

In this study, we analyzed the association between smoking and the concentration of metabolites in 1,241 serum samples from the KORA baseline survey 4 (S4) and follow-up (F4) study, aiming to extend the knowledge of smoking-associated metabolites beyond our pilot study by including female CS at two time points over seven years, to investigate whether smoking-associated changes in metabolite profile are reversible after smoking cessation, and to provide insights into the pathophysiological consequences of smoking in protein-metabolite networks.


Ethics statement

Written informed consent was obtained from KORA S4 and F4 participants. The KORA study was approved by the ethics committee of the Bavarian Medical Association in Munich, Germany.

Study population

The KORA surveys are population-based studies conducted in the Region of Augsburg in Germany [27, 28]. Four surveys were conducted with 18,079 participants recruited from 1984 to 2001. The S4 consists of 4,261 individuals (25 to 74 years old) examined from 1999 to 2001. From 2006 to 2008, 3,080 participants (with an age range of 32 to 81 years) took part in the F4 survey. Each participant completed a lifestyle questionnaire providing information on a number of parameters including smoking status (current, former, never). Serum samples for metabolomics analysis were collected in parallel in the KORA S4 and F4 survey as described elsewhere [2931].

For metabolite profiles, serum samples from 1,614 people aged 55 to 74 years old were available [29]. Participants with non-fasting status (N = 216) or missing values (N = 22) were excluded from the analysis. We further excluded 145 people in KORA S4 and 116 people in the longitudinal data of KORA S4 → F4, whose spouses were CS, to rule out passive smoking effects. Furthermore, metabolite concentrations of serum samples from 1,036 participants were measured in both KORA S4 and F4.

Metabolite measurements

Liquid handling of serum samples (10 μl) was performed with Hamilton star robot (Hamilton Bonaduz AG, Bonaduz, Switzerland) and prepared for quantification using the AbsoluteIDQ P180 and P150 kits (BIOCRATES Life Science AG, Innsbruck, Austria) for the KORA S4 and F4 surveys, respectively. This allowed simultaneous quantification of 188 or 163 metabolites using liquid chromatography and flow injection analysis mass spectrometry as described previously [32, 33]. The complete analytical process was monitored by quality control steps, reference samples and the MetIQ software package, which is an integral part of the Absolute IDQ kit.

Because the two datasets were generated by different platforms, different quality control processes were introduced. The metabolite data quality control procedure for the KORA S4 samples was described in our recently published work [29]. There were 140 metabolites that passed the two quality controls: one hexose, 21 amino acids, eight biogenic amines, 21 acylcarnitines, 13 sphingomyelins (SMs), eight lysoPCs, 33 diacyl-PCs (PC aa Cx:y) and 35 acyl-alkyl-PCs (PC ae Cx:y). Lipid side chain composition is abbreviated as Cx:y, where × denotes the number of carbons in the side chain and y the number of double bonds. The precise position of the double bonds and the distribution of the carbon atoms in different fatty acid side chains cannot be determined with this technology. Concentrations of all analyzed metabolites are reported in μmol/L (μM). The data cleaning procedure for the KORA F4 samples has previously been described in detail [24, 30]. In total, 121 metabolites were measured in both S4 and F4, and used for the prospective study.

Statistical analysis

Differences in population characteristics (CS, FS and NS) were tested by a two-tailed student's t-test. The metabolite concentrations were log transformed for normalization. We tested cross-sectional association of each metabolite with smoking using logistic regression models adjusted for age, body mass index (BMI) and alcohol consumption (see Figure 1). To correct for multiple testing, false discovery rate (FDR) was calculated using the Benjamini-Hochberg method [34] and the cut-off for statistical significance was set at FDR <0.05.

Figure 1
figure 1

Flow diagram illustrating the analysis strategy. CS: current smokers; FS: former smokers; NS: never smokers.

Linear regression models were used to investigate whether smoking intensities measured in pack years and cessation time are associated with metabolite concentrations. In the case of CS, the years of smoking were calculated as the time period from starting smoking until the start of the survey. Pack year was calculated as the number of cigarettes per day multiplied by smoking duration and divided by 20 [35]. Cessation time (in years) was calculated according to the questionnaire. The models contained the log-transformed metabolite concentrations as the dependent variable and the smoking intensities as the explanatory variable, with age, BMI and alcohol consumption as covariates. Every unit change of one covariate corresponds to a relative change of the metabolite concentration by Δ (%):

Δ = ( exp ( β i ) - 1 ) × 100 %

where βi indicates the estimate of ith covariate in the model.

To assess the role of smoking cessation for the quitters, who were CS at S4 but FS at F4, we fitted the linear mixed models to the longitudinal data of KORA S4 → F4. The models contained the fixed effect of smoking status (CS, FS and NS), age, BMI and alcohol consumption with a random effect assigned to each participant. All calculations were performed in R (version 2.14.1).

Network and pathway analysis

We retrieved protein-protein interactions from the databases of the Search Tool for the Retrieval of Interacting Genes/Proteins [36] and the relationships between enzymes and metabolites from the Human Metabolome Database [37] to construct protein-metabolite networks containing links between metabolites, enzymes and smoking-related genes. Genes and metabolites were connected allowing for at most one intermediate enzyme by Dijkstra's algorithm [38], and optimized by eliminating edges with Search Tool for the Retrieval of Interacting Genes/Proteins scores less than 0.7. Each edge in the networks was manually checked. We have implemented this method in our previous studies [29, 39]. The analysis was performed using the R package igraph [40]. The network was visualized using Cytoscape [41]. Pathway analysis was performed by MetaboAnalyst [42].


Characteristics of participants of the cross-sectional KORA S4

Participants were divided into three groups according to their self-reported smoking status. Population characteristics are shown in Table 1. On average, CS were two to three years younger and had a lower BMI than FS and NS. Male CS showed higher alcohol consumption than male NS, but there was no significant difference observed in women. Furthermore, the statistics showed differences in lifestyle factors between men and women. Alcohol consumption was higher in men than women (P = 1.5e-11 (CS); P = 2.2e-18 (FS); P = 9.5e-17 (NS)), and smoking intensity (in pack years) was higher in male than in female CS (P = 6.0e-6).

Table 1 Characteristics of cross-sectional KORA S4.

Metabolomic differences between current, former and never smokers

We identified 18 metabolites in men and six in women that were significantly different (FDR <0.05) between CS and NS. Three metabolites (PC ae C34:3, PC aa C36:1 and glutamate) were identified in both men and women showing the same pattern of variation (higher or lower) (Table 2). Compared with FS and NS, in male CS the concentrations of four unsaturated diacyl-PCs (PC aa C34:1, PC aa C36:1, PC aa C38:3 and PC aa C40:4) and five amino acids (arginine, aspartate, glutamate, ornithine and serine) were higher, whereas three saturated diacyl-PCs, one lysoPC and four acyl-alkyl-PCs, as well as kynurenine showed lower concentrations. In female CS, we found higher levels of carnitine and PC aa C32:1, and a lower level of hydroxysphingomyeline (SM (OH)) C22:2.

Table 2 Smoking-related metabolites in KORA S4.

Among the 21 smoking-related metabolites (18 in men and six in women), 19 were found to be reversible (that is, significant difference between FS and CS but without significant difference between FS and NS; FDR <0.05). No irreversible metabolite was observed (that is, significant difference between FS and NS). Serine and PC aa C32:3 in men were not classified because their concentrations were not significantly different between CS and FS or between FS and NS (Table 2). A heat map representing the concentration profiles of the 21 identified metabolites in CS, FS and NS is shown in Figure 2, demonstrating the reversibility of metabolites after smoking cessation.

Figure 2
figure 2

Heat maps of smoking-related metabolites in (A) men and (B) women. The heat map shows mean residues of smoking-related metabolites in CS, FS and NS and the reversibility after smoking cessation. The color of each cell in the heat map represents the relative mean concentration of each metabolite in NS, FS or CS. The number of samples in each group is provided. The bar besides the metabolite names indicates the reversibility of these metabolites after smoking cessation. aa: diacyl-; ae: acyl-alkyl-; C0: carnitine; CS: current smokers; FS: former smokers; lysoPC: acyl-phosphatidylcholine; NS: never smokers; PC: phosphatidylcholine; SM (OH): hydroxysphingomyeline.

In women, SM (OH) C22:2 was significantly associated with cessation time (FDR <0.05); however, there was no such significant metabolite in men (Table S1 in Additional file 1), indicating a non-linear relationship between cessation time and the reversion of metabolite profile. In addition, we grouped the FS by stratified cessation years (0 to 10, 11 to 20, 21 to 30, 31 to 40, over 40 years). For some metabolites (for example, PC ae C38:0, PC aa C36:0 and ornithine), the greatest change of concentration occurred within the first 10 years of cessation compared with CS (Figure 3).

Figure 3
figure 3

Metabolite concentration variations in relation to smoking cessation time. Taking NS as baseline, figures show the mean residuals of metabolites in different groups of CS and FS, giving the trend of metabolite variation with cessation time. FS were grouped by stratified cessation time (≤10, 11 to 20, 21 to 30, 31 to 40, 41+). Residuals were calculated by linear regression model (regression of metabolite concentration on age, body mass index and alcohol consumption). aa: diacyl-; ae: acyl-alkyl-; CS: current smokers; FS: former smokers; NS: never smokers; PC: phosphatidylcholine.

Within CS, we found kynurenine and PC ae C34:3, PC ae C38:0 and PC ae C38:6 in men, and PC aa C36:1 in women showing significant association with pack years. In the linear regression model, pack years showed a negative relation (parameter estimation β <0) to these five metabolites (Table 3) (for example, one pack year increase will lead to a decrease of the kynurenine level in CS by 0.33%).

Table 3 Smoking intensity (pack years) related to metabolites

Prospective change of metabolite profiles (from KORA baseline S4 to follow-up F4)

The prospective dataset included 40 CS, 432 NS and 49 quitters (people who were CS in KORA S4 but FS in KORA F4) (Table 4). Among the 16 reversible metabolites in men, 13 (except kynurenine, glutamate and aspartate) were also measured in KORA F4 using a different kit (see Methods). We employed a linear mixed effect model to investigate the effects of smoking cessation on metabolite concentrations. Among these 13 metabolites, 10 metabolites showed a significant variation in quitters, with a period of smoking cessation from one to seven years, which indicated a reverting process. The arginine level decreased by 11.3% and ornithine by 14.8% in quitters compared with CS, whereas PC aa C36:0 increased by 18.5%. Figure 4 shows the prospective changes of the significant metabolites. For women, the same analysis was conducted. Because the number of female quitters was small (N = 10), five metabolites that were measured in both KORA S4 and F4 showed borderline significance (P <0.05). However, none of these metabolites was found to be significant considering FDR <0.05 (see Table 5).

Table 4 Characteristics of the prospective dataset (KORA S4 → F4).
Figure 4
figure 4

Changes of smoking-related metabolites in current, former and never smokers in KORA S4 → F4. Taking the NS as baseline, the concentration change of each metabolite is shown as the adjusted mean residue in KORA S4 and F4 in all three groups (CS, FS and NS). Only metabolites with significant prospective change in KORA S4 F4 are shown in the figure. Residuals were calculated from a linear regression model (regression of metabolite concentration on age, body mass index and alcohol consumption). aa: diacyl-; ae: acyl-alkyl-; CS: current smokers; FS: former smokers; lysoPC: acyl-phosphatidylcholine; NS: never smokers; PC: phosphatidylcholine.

Table 5 Association of reversible metabolites with smoking status change in the prospective dataset (KORA S4 → F4)

Smoking effects on metabolic network

Enrichment analysis of the 21 identified smoking-related metabolites on Kyoto Encyclopedia of Genes and Genomes pathways showed enrichment in a set of amino acid and lipid metabolism pathways (ether lipid, glycerophospholipid, arginine and proline metabolism). In addition, we analyzed the impact of the smoking-related metabolites in each pathway by measuring their structural importance (see Methods). These metabolites had high betweenness centrality and a strong impact on the enriched pathways (Figure 5 and Table S2 in Additional file 2).

Figure 5
figure 5

Pathway analyses of smoking-related metabolites. Figure shows enrichment and impact of smoking-related metabolites in Kyoto Encyclopedia of Genes and Genomes pathways. The enrichment scores are shown on y-axis, which was calculated as the negative logarithm of the P-value from an enrichment test. The x-axis indicates the structural impact with a score from 0 to 1 of the smoking-related metabolites in the enriched pathways.

To systematically investigate how the effects of smoking propagate over the metabolic networks, we evaluated the association between 175 smoking-related genes, previously reported [23], and the 21 smoking-related metabolites we found in this study by analyzing protein-metabolite networks (see Methods). In men, 15 metabolites (lysoPC a C18:2, PC aa C32:3,PC aa C34:1, PC aa C36:0, PC aa C36:1, PC aa C38:0, PC aa C38:3, PC aa C40:4, PC ae C34:3, PC ae C38:0, PC ae C38:6, PC ae C40:6, arginine, glutamate and serine) were found to be linked with 11 genes (ADH7, AKR1B1, DHRS3, FTL, GALE, GPC1, KRAS, S100A10, SLC7A11, SULF1, PLA2G10) by related enzymes. In women, four metabolites (PC aa C36:1, PC ae C34:3, PC aa C32:1 and glutamate) were closely linked with nine genes (ADH7, AKR1B1, DHRS3, FTL, GALE, GPC1, S100A10, SULF1, PLA2G10) (Figure 6A and Table S3 in Additional file 3). Similar to enrichment analysis, the network in men and in women could be generally divided into glycerophospholipids and tightly associated proteins as well as amino acids and the associated genes and enzymes. A description of the protein-metabolite and protein-protein interactions was listed in Table S3 in Additional file 3.

Figure 6
figure 6

Protein-metabolite networks and pathways of the smoking-related metabolites and genes. (A) Network linking metabolites and proteins encoded by smoking-related genes with maximum one intermediate. Node color indicates the reversibility after smoking cessation. (B, C) Effects of smoking on arginine and glutamate as well as on lipid metabolism. Metabolites are in regular font, protein coding genes are in italic, gender-specific gene (CPS1) is in bold italic font. aa: diacyl-; ae: acyl-alkyl-; APOA5: apolipoprotein A-V; BDH: 3-hydroxybutyrate dehydrogenase, type 1; cPLA2: cytosolic phospholipase A2; CS: current smokers; FS: former smokers; GIIC sPLA2: phospholipase A2, membrane associated; LRAT: lecithin retinol acyltransferase; LYPLA1: lysophospholipase I; lysoPC: acyl-phosphatidylcholine; NOS1: nitric oxide synthase 1; NS: never smokers; PC: phosphatidylcholine; PLA2G10: group 10 secretory phospholipase A2; SCGB1A1: uteroglobin; SDH: serine dehydratase; SLC3A2: solute carrier family 3 member 2

The smoking effects on the networks were reversible. With regards to gene expressions, with the exception of SULF1 and PLA2G10, all changes in the networks were reversible after smoking cessation [23]. All changes in metabolites in the network were also reversible, except serine.


In this study, we have used an 'omics' approach to investigate the association of metabolite concentrations with smoking, delineated the reversion of metabolite variations after smoking cessation and demonstrated the results using protein-metabolite networks. We identified strong associations of various metabolites with smoking, and confirmed part of the findings of our pilot study [1]. Among the 23 smoking-related metabolites identified in the pilot study, 11 metabolites were measured in this study, five of which (four unsaturated diacyl-PCs and one acyl-alkyl-PC) were validated in men, based on about five-fold larger CS samples. Consistent patterns of smoking effects on metabolite profile were observed in the current study. Among all the smoking-related metabolites, in CS we found higher unsaturated diacyl-PCs, but lower acyl-alkyl-PCs and saturated diacyl-PCs, which may indicate generally increased levels of unsaturated fatty acids in CS. Unsaturated fatty acids are more vulnerable to lipid peroxidation and influence the risk of different diseases [43, 44].

Smoking-related metabolites and cardiovascular disease

The study results implied the potential of metabolomics in revealing the role of an environmental factor, for example a smoking lifestyle, in the pathogenesis and prognosis of CVD.

One study on the peripheral blood metabolite profile showed an association of coronary artery disease and urea cycle-related metabolites, including arginine and glutamate [45], which were also identified in our study as smoking-related metabolites. By scrutinizing the smoking-related metabolites in metabolic pathways, we found further support for the pathophysiological relation between these metabolites and CVD. Previous findings indicated that the glutamate transporter in human lung epithelial cells, encoded by the SLC7A11 gene, is activated in CS [23, 46], which increases the transportation of glutamate and subsequently raises the levels of the downstream metabolites, arginine and ornithine (Figure 6B). The activation of the cysteine-glutamate transporter (encoded by SLC7A11) and the increased glutamate level as a response to oxidative stress is also of great importance to endothelial dysfunction involved at all stages of atherosclerotic plaque evolution, which leads to CVD [47, 48].

Ether lipid and glycerophospholipid metabolisms are associated with smoking [1, 49]. The decreased level of lysoPC a C18:2 reflects the inhibition of upstream synthesis and activation of downstream hydrolysis. As shown in Figure 6C, upregulation of S100A10 and GPC1 inhibits cytosolic phospholipase A2, which plays a role in the synthesis of lyso-PCs. The lysophospholipase I isoform, which hydrolyses lysoPC into glycerophosphocholine, is upregulated in CS [23]. Interestingly, one recent study showed that a disorder of phosphatidylcholine metabolism would promote CVD [50], which may establish a link between smoking-related phosphatidylcholine variation and cardiovascular events. For example, the phosphatidylcholine hydroperoxide will promote angiogenesis in endothelial cells that are associated with atherosclerotic development [51].

The reversibility of metabolite concentrations in a small time window may reveal a reduced risk of smoking-related diseases after stopping smoking. Concentrations of arginine and glutamate that are associated with both smoking and coronary artery diseases quickly returned to normal levels (within seven years) after smoking cessation, which is in line with epidemiological findings that the smoking effects on CVD are quickly and largely reduced after smoking cessation [8, 9, 52]. The reversed glutamate level indicates reduced oxidative stress after smoking cessation, and the reversion of arginine and ornithine reflects a reversion of functioning in the urea cycle. Our findings provide metabolic insight into the reduced risk of CVD after smoking cessation and provide support for the remarkable benefits people would gain by stopping smoking.

Concordance of reversibility in metabolic network

The protein-metabolite interaction network shows that the reversibility of metabolite concentrations also coincided with gene expression (Figure 6A). Arginine and glutamate were quickly reversed after smoking cessation, which was in line with the quick reversibility of SLC7A11 expression. Expression of enzyme coding genes for the hydrolysis of diacyl-PCs and acyl-alkyl-PCs, for instance lysophospholipase, cytosolic phospholipase A2 and S100 calcium binding protein A2, were quickly reversible and smoking-related diacyl-PCs and acyl-alkyl PCs shared the same reverse pattern.

Gender-specific effects of smoking

In this study, we found gender-specific effects of smoking on metabolite profiles (Table S1 in Additional file 1). This result supports the assumption that differences in smoking effects on men and women are not solely based on smoking intensity but are also gender-specific. Glutamate was higher in both male and female CS, however, the levels of arginine and ornithine were only higher in male CS. According to a previous study of the metabolomic and genetic biomarkers on sexual dimorphisms [30], the CPS1 gene, which regulates the formation of arginine, has a gender-specific manner in certain single nucleotide polymorphisms, with stronger effects in women than in men. The gender-specific genetic effect might cause a lower efficiency in women in regard to the transformation of extra glutamate to citrulline (Figure 6C).

Strengths and limitations

We used a systematic targeted metabolomics approach with 140 metabolites in a large population-based cohort. Analyzing the effects of smoking and smoking cessation in this prospective manner (follow-up of seven years) provides more power to investigate smoking effects by ruling out individual differences. However, our study is based on a limited range and number of metabolites and cannot fully represent the whole metabolome. Thus, an improved metabolomics technique measuring more metabolites is urgently needed for a comprehensive understanding of both reversible and permanent effects of smoking on human metabolism. It would be interesting for future studies to also include data on other environmental factors such as diet and lifestyle, which are known to have effects on the human metabolome [53, 54].


Our study shows the power of the metabolomics approach in investigating the molecular signature of lifestyle-related environmental exposures. We demonstrated that smoking is associated with concentration variations in amino acids, ether lipid and glycerophospholipid metabolism at an 'omics' level. The smoking-related changes in the human serum metabolite profile are reversible after stopping smoking. This indicates the remarkable benefits of smoking cessation and provides a link to CVD benefits. Furthermore, linking metabolomic knowledge to other 'omics' approaches, for example, transcriptomics, may have the potential to identify novel biomarkers as well as new risk assessment tools.







body mass index


current smokers


cardiovascular disease


false discovery rate


former smokers




never smokers





SM (OH):



  1. Wang-Sattler R, Yu Y, Mittelstrass K, Lattka E, Altmaier E, Gieger C, Ladwig KH, Dahmen N, Weinberger KM, Hao P, Liu L, Li Y, Wichmann H-E, Adamski J, Suhre K, Illig T: Metabolic profiling reveals distinct variations linked to nicotine consumption in humans--first results from the KORA study. PLoS One. 2008, 3: e3863-10.1371/journal.pone.0003863.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Xavier F, Henn L de A, Oliveira M, Orlandine L: Smoking and its relation to the histological type, survival, and prognosis among patients with primary lung cancer. Sao Paulo Med J. 1996, 114: 1298-1302.

    Article  CAS  PubMed  Google Scholar 

  3. Smoking-attributable mortality, years of potential life lost, and productivity losses--United States, 2000-2004. Morbidity and Mortality Weekly Report. 2008, 57: 1226-1228.

  4. Bøttcher M, Falk E: Pathology of the coronary arteries in smokers and non-smokers. J Cardiovasc Risk. 1999, 6: 299-302.

    Article  PubMed  Google Scholar 

  5. Howard G, Burke GL, Szklo M, Tell GS, Eckfeldt J, Evans G, Heiss G: Active and passive smoking are associated with increased carotid wall thickness. The Atherosclerosis Risk in Communities Study. Arch Intern Med. 1994, 154: 1277-1282. 10.1001/archinte.1994.00420110125014.

    Article  CAS  PubMed  Google Scholar 

  6. Libby P: Inflammation in atherosclerosis. Nature. 2002, 420: 868-874. 10.1038/nature01323.

    Article  CAS  PubMed  Google Scholar 

  7. Szmitko PE, Wang C-H, Weisel RD, Jeffries GA, Anderson TJ, Verma S: Biomarkers of vascular disease linking inflammation to endothelial activation: Part II. Circulation. 2003, 108: 2041-2048. 10.1161/01.CIR.0000089093.75585.98.

    Article  PubMed  Google Scholar 

  8. Danesh J, Whincup P, Walker M, Lennon L, Thomson A, Appleby P, Gallimore JR, Pepys MB: Low grade inflammation and coronary heart disease: prospective study and updated meta-analyses. BMJ. 2000, 321: 199-204. 10.1136/bmj.321.7255.199.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Danesh J, Muir J, Wong YK, Ward M, Gallimore JR, Pepys MB: Risk factors for coronary heart disease and acute-phase proteins. A population-based study. Eur Heart J. 1999, 20: 954-959. 10.1053/euhj.1998.1309.

    Article  CAS  PubMed  Google Scholar 

  10. Rigotti NA, Pasternak RC: Cigarette smoking and coronary heart disease: risks and management. Cardiol Clin. 1996, 14: 51-68. 10.1016/S0733-8651(05)70260-5.

    Article  CAS  PubMed  Google Scholar 

  11. Cook DG, Shaper AG, Pocock SJ, Kussick SJ: Giving up smoking and the risk of heart attacks. A report from The British Regional Heart Study. Lancet. 1986, 2: 1376-1380.

    Article  CAS  PubMed  Google Scholar 

  12. Aberg A, Bergstrand R, Johansson S, Ulvenstam G, Vedin A, Wedel H, Wilhelmsson C, Wilhelmsen L: Cessation of smoking after myocardial infarction. Effects on mortality after 10 years. Br Heart J. 1983, 49: 416-422. 10.1136/hrt.49.5.416.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Halpern MT, Gillespie BW, Warner KE: Patterns of absolute risk of lung cancer mortality in former smokers. J Natl Cancer Inst. 1993, 85: 457-464. 10.1093/jnci/85.6.457.

    Article  CAS  PubMed  Google Scholar 

  14. Ebbert JO, Yang P, Vachon CM, Vierkant RA, Cerhan JR, Folsom AR, Sellers TA: Lung cancer risk reduction after smoking cessation: observations from a prospective cohort of women. J Clin Oncol. 2003, 21: 921-926. 10.1200/JCO.2003.05.085.

    Article  CAS  PubMed  Google Scholar 

  15. Bakhru A, Erlinger TP: Smoking cessation and cardiovascular disease risk factors: results from the Third National Health and Nutrition Examination Survey. PLoS Med. 2005, 2: e160-10.1371/journal.pmed.0020160.

    Article  PubMed  PubMed Central  Google Scholar 

  16. De Maat MP, Pietersma A, Kofflard M, Sluiter W, Kluft C: Association of plasma fibrinogen levels with coronary artery disease, smoking and inflammatory markers. Atherosclerosis. 1996, 121: 185-191. 10.1016/0021-9150(95)05716-1.

    Article  CAS  PubMed  Google Scholar 

  17. Pradhan AD, Manson JE, Rossouw JE, Siscovick DS, Mouton CP, Rifai N, Wallace RB, Jackson RD, Pettinger MB, Ridker PM: Inflammatory biomarkers, hormone replacement therapy, and incident coronary heart disease: prospective analysis from the Women's Health Initiative observational study. JAMA. 2002, 288: 980-987. 10.1001/jama.288.8.980.

    Article  CAS  PubMed  Google Scholar 

  18. Yanbaeva DG, Dentener MA, Creutzberg EC, Wesseling G, Wouters EFM: Systemic effects of smoking. Chest. 2007, 131: 1557-1566. 10.1378/chest.06-2179.

    Article  CAS  PubMed  Google Scholar 

  19. Powell CA, Klares S, O'Connor G, Brody JS: Loss of heterozygosity in epithelial cells obtained by bronchial brushing: clinical utility in lung cancer. Clin Cancer Res. 1999, 5: 2025-2034.

    CAS  PubMed  Google Scholar 

  20. Wistuba II, Lam S, Behrens C, Virmani AK, Fong KM, LeRiche J, Samet JM, Srivastava S, Minna JD, Gazdar AF: Molecular damage in the bronchial epithelium of current and former smokers. J Natl Cancer Inst. 1997, 89: 1366-1373. 10.1093/jnci/89.18.1366.

    Article  CAS  PubMed  Google Scholar 

  21. Wistuba II, Mao L, Gazdar AF: Smoking molecular damage in bronchial epithelium. Oncogene. 2002, 21: 7298-7306. 10.1038/sj.onc.1205806.

    Article  CAS  PubMed  Google Scholar 

  22. Guo M, House MG, Hooker C, Han Y, Heath E, Gabrielson E, Yang SC, Baylin SB, Herman JG, Brock MV: Promoter hypermethylation of resected bronchial margins: a field defect of changes?. Clin Cancer Res. 2004, 10: 5131-5136. 10.1158/1078-0432.CCR-03-0763.

    Article  CAS  PubMed  Google Scholar 

  23. Beane J, Sebastiani P, Liu G, Brody JS, Lenburg ME, Spira A: Reversible and permanent effects of tobacco smoke exposure on airway epithelial gene expression. Genome Biol. 2007, 8: R201-10.1186/gb-2007-8-9-r201.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Yu Z, Zhai G, Singmann P, He Y, Xu T, Prehn C, Römisch-Margl W, Lattka E, Gieger C, Soranzo N, Heinrich J, Standl M, Thiering E, Mittelstraß K, Wichmann HE, Peters A, Suhre K, Li Y, Adamski J, Spector TD, Illig T, Wang-Sattler R: Human serum metabolic profiles are age dependent. Aging cell. 2012, 11: 960-967. 10.1111/j.1474-9726.2012.00865.x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ellis JK, Athersuch TJ, Thomas LD, Teichert F, Pérez-Trujillo M, Svendsen C, Spurgeon DJ, Singh R, Järup L, Bundy JG, Keun HC: Metabolic profiling detects early effects of environmental and lifestyle exposure to cadmium in a human population. BMC Med. 2012, 10: 61-10.1186/1741-7015-10-61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Vulimiri SV, Misra M, Hamm JT, Mitchell M, Berger A: Effects of mainstream cigarette smoke on the global metabolome of human lung epithelial cells. Chem Res Toxicol. 2009, 22: 492-503. 10.1021/tx8003246.

    Article  CAS  PubMed  Google Scholar 

  27. Holle R, Happich M, Löwel H, Wichmann HE: KORA--a research platform for population based health research. Gesundheitswesen. 2005, 67 (Suppl 1): S19-25.

    Article  PubMed  Google Scholar 

  28. Wichmann H-E, Gieger C, Illig T: KORA-gen--resource for population genetics, controls and a broad spectrum of disease phenotypes. Gesundheitswesen. 2005, 67 (Suppl 1): S26-30.

    Article  PubMed  Google Scholar 

  29. Wang-Sattler R, Yu Z, Herder C, Messias AC, Floegel A, He Y, Heim K, Campillos M, Holzapfel C, Thorand B, Grallert H, Xu T, Bader E, Huth C, Mittelstrass K, Döring A, Meisinger C, Gieger C, Prehn C, Roemisch-Margl W, Carstensen M, Xie L, Yamanaka-Okumura H, Xing G, Ceglarek U, Thiery J, Giani G, Lickert H, Lin X, Li Y, et al: Novel biomarkers for pre-diabetes identified by metabolomics. Mol Syst Biol. 2012, 8: 615.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Mittelstrass K, Ried JS, Yu Z, Krumsiek J, Gieger C, Prehn C, Roemisch-Margl W, Polonikov A, Peters A, Theis FJ, Meitinger T, Kronenberg F, Weidinger S, Wichmann HE, Suhre K, Wang-Sattler R, Adamski J, Illig T: Discovery of sexual dimorphisms in metabolic and genetic biomarkers. PLoS Genet. 2011, 7: e1002215-10.1371/journal.pgen.1002215.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Jourdan C, Petersen A-K, Gieger C, Döring A, Illig T, Wang-Sattler R, Meisinger C, Peters A, Adamski J, Prehn C, Suhre K, Altmaier E, Kastenmüller G, Römisch-Margl W, Theis FJ, Krumsiek J, Wichmann H-E, Linseisen J: Body fat free mass is associated with the serum metabolite profile in a population-based study. PLoS One. 2012, 7: e40009-10.1371/journal.pone.0040009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Illig T, Gieger C, Zhai G, Römisch-Margl W, Wang-Sattler R, Prehn C, Altmaier E, Kastenmüller G, Kato BS, Mewes H-W, Meitinger T, Angelis MH de, Kronenberg F, Soranzo N, Wichmann H-E, Spector TD, Adamski J, Suhre K: A genome-wide perspective of genetic variation in human metabolism. Nat Genet. 2010, 42: 137-141. 10.1038/ng.507.

    Article  CAS  PubMed  Google Scholar 

  33. Römisch-Margl W, Prehn C, Bogumil R, Röhring C, Suhre K, Adamski J: Procedure for tissue sample preparation and metabolite extraction for high-throughput targeted metabolomics. Metabolomics. 2012, 8: 133-142. 10.1007/s11306-011-0293-4.

    Article  Google Scholar 

  34. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.

    Google Scholar 

  35. NCI Dictionary of Cancer Terms. []

  36. Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, Jensen LJ, von Mering C: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011, 39: D561-568. 10.1093/nar/gkq973.

    Article  CAS  PubMed  Google Scholar 

  37. Wishart DS, Knox C, Guo AC, Eisner R, Young N, Gautam B, Hau DD, Psychogios N, Dong E, Bouatra S, Mandal R, Sinelnikov I, Xia J, Jia L, Cruz JA, Lim E, Sobsey CA, Shrivastava S, Huang P, Liu P, Fang L, Peng J, Fradette R, Cheng D, Tzur D, Clements M, Lewis A, De Souza A, Zuniga A, Dawe M, et al: HMDB: a knowledgebase for the human metabolome. Nucleic Acids Res. 2009, 37: D603-610. 10.1093/nar/gkn810.

    Article  CAS  PubMed  Google Scholar 

  38. Dijkstra EW: A note on two problems in connexion with graphs. Numerische Mathematik. 1959, 1: 269-271. 10.1007/BF01386390.

    Article  Google Scholar 

  39. He Y, Yu Z, Giegling I, Xie L, Hartmann AM, Prehn C, Adamski J, Kahn R, Li Y, Illig T, Wang-Sattler R, Rujescu D: Schizophrenia shows a unique metabolomics signature in plasma. Transl Psychiatry. 2012, 2: e149-10.1038/tp.2012.76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Csardi G, Nepusz T: The igraph software package for complex network research. InterJournal. 2006, Complex Systems: 1695.

    Google Scholar 

  41. Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27: 431-432. 10.1093/bioinformatics/btq675.

    Article  CAS  PubMed  Google Scholar 

  42. Xia J, Mandal R, Sinelnikov IV, Broadhurst D, Wishart DS: MetaboAnalyst 2.0--a comprehensive server for metabolomic data analysis. Nucleic Acids Res. 2012, 40: W127-133. 10.1093/nar/gks374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Kris-Etherton PM: AHA Science Advisory. Monounsaturated fatty acids and risk of cardiovascular disease. American Heart Association. Nutrition Committee. Circulation. 1999, 100: 1253-1258. 10.1161/01.CIR.100.11.1253.

    Article  CAS  PubMed  Google Scholar 

  44. Mozaffarian D, Rimm EB, Herrington DM: Dietary fats, carbohydrate, and progression of coronary atherosclerosis in postmenopausal women. Am J Clin Nutr. 2004, 80: 1175-1184.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Shah SH, Bain JR, Muehlbauer MJ, Stevens RD, Crosslin DR, Haynes C, Dungan J, Newby LK, Hauser ER, Ginsburg GS, Newgard CB, Kraus WE: Association of a peripheral blood metabolic profile with coronary artery disease and risk of subsequent cardiovascular events. Circ Cardiovasc Genet. 2010, 3: 207-214. 10.1161/CIRCGENETICS.109.852814.

    Article  CAS  PubMed  Google Scholar 

  46. Bridges CC, Kekuda R, Wang H, Prasad PD, Mehta P, Huang W, Smith SB, Ganapathy V: Structure, function, and regulation of human cystine/glutamate transporter in retinal pigment epithelial cells. Invest Ophthalmol Vis Sci. 2001, 42: 47-54.

    CAS  PubMed  Google Scholar 

  47. Harrison D, Griendling KK, Landmesser U, Hornig B, Drexler H: Role of oxidative stress in atherosclerosis. Am J Cardiol. 2003, 91: 7-11. 10.1016/S0002-9149(02)03144-2.

    Article  Google Scholar 

  48. Glass CK, Witztum JL: Atherosclerosis. The road ahead. Cell. 2001, 104: 503-516. 10.1016/S0092-8674(01)00238-0.

    Article  CAS  PubMed  Google Scholar 

  49. Palmerini CA, Saccardi C, Ferracci F, Arienti S: Lipid patterns in the saliva of smoking young adults. Hum Exp Toxicol. 2011, 30: 1482-1488. 10.1177/0960327111398672.

    Article  CAS  PubMed  Google Scholar 

  50. Wang Z, Klipfell E, Bennett BJ, Koeth R, Levison BS, Dugar B, Feldstein AE, Britt EB, Fu X, Chung Y-M, Wu Y, Schauer P, Smith JD, Allayee H, Tang WHW, DiDonato JA, Lusis AJ, Hazen SL: Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature. 2011, 472: 57-63. 10.1038/nature09922.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Nakagawa K, Shibata A, Saito T, Sookwong P, Kato S, Tsuduki T, Matsubara K, Miyazawa T: Phosphatidylcholine hydroperoxide promotes VEGF-induced angiogenesis in endothelial cells and rat aorta ring cultures. Biochim Biophys Acta. 2011, 1810: 1205-1211. 10.1016/j.bbagen.2011.08.018.

    Article  CAS  PubMed  Google Scholar 

  52. Danesh J, Collins R, Appleby P, Peto R: Association of fibrinogen, C-reactive protein, albumin, or leukocyte count with coronary heart disease: meta-analyses of prospective studies. JAMA. 1998, 279: 1477-1482. 10.1001/jama.279.18.1477.

    Article  CAS  PubMed  Google Scholar 

  53. Gibney MJ, Walsh M, Brennan L, Roche HM, German B, van Ommen B: Metabolomics in human nutrition: opportunities and challenges. Am J Clin Nutr. 2005, 82: 497-503.

    CAS  PubMed  Google Scholar 

  54. Pohjanen E, Thysell E, Jonsson P, Eklund C, Silfver A, Carlsson I-B, Lundgren K, Moritz T, Svensson MB, Antti H: A multivariate screening strategy for investigating metabolic effects of strenuous physical exercise in human serum. J Proteome Res. 2007, 6: 2113-2120. 10.1021/pr070007g.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references


We express our appreciation to all KORA study participants for donating their blood and time. We thank the field staff in Augsburg who conducted the KORA studies. The KORA research platform and the KORA Augsburg studies were initiated and financed by the Helmholtz Zentrum München, which is funded by the German Federal Ministry of Education, Science, Research and Technology and by the State of Bavaria. The KORA study group consists of A. Peters (speaker), J. Heinrich, R. Holle, R. Leidl, C. Meisinger, K. Strauch and their co-workers, who are responsible for the design and conduct of the KORA studies. We thank Julia Scarpa, Katharina Sckell and Arsin Sabunchi for metabolomics measurements performed at the Helmholtz Zentrum München, Genome Analysis Centre, Metabolomics

Core Facility. This study was supported in part by a grant from the German Federal Ministry of Education and Research (BMBF) to the German Centre for Diabetes Research (DZD e.V.) and this work was partly supported by the BMBF project 'Metabolomics of ageing' (FKZ: 01DO12030) and EU FP7 grant HEALTH-2009-2.2.1-3/242114 (Project OPTiMiSE).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Rui Wang-Sattler.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HEW, KS, JA, TI, AP and RWS initiated and designed the study. CP, WRM, WR, HEW, KHL and JA were involved in and performed the experiment. TX and ZY performed the data analysis. TX, CH, ZY and RWS wrote the manuscript, XD, EB, CP, KP, MJ, YL, HW, FT, JA and AP revised the manuscript. The manuscript has been approved by all authors.

Electronic supplementary material


Additional file 1: Table S1: Cessation time-related metabolites in FS. FDR was calculated by P-value adjusted for the number of smoking-related metabolites with Benjamini-Hochberg method. aa: diacyl-; ae: acyl-alkyl-; C0: carnitine; FS: former smokers; lysoPC: acyl-phosphatidylcholine; PC: phosphatidylcholine; SM (OH): hydroxysphingomyeline. (DOC 42 KB)


Additional file 2: Table S2: Enrichment and impact of smoking-related metabolites in Kyoto Encyclopedia of Genes and Genomes pathways. Table shows the enrichment and impact scores of smoking-related metabolites in Kyoto encyclopedia of Genes and Genomes pathways. The pathway analysis consists of enrichment and a structural impact analysis both based on Kyoto Encyclopedia of Genes and Genomes database. The -log (P) was considered as the enrichment score. Impact, scored between 0 and 1, indicated the pathway topological importance of the metabolites. In particular, the parameter Total is the total number of compounds in the pathway; the parameter Hits is the actual number of metabolites with significant variations in the pathway; the Raw P was the original P-value calculated from the enrichment analysis; the FDR was calculated as the P-value adjusted using Benjamini-Hochberg method. (DOC 56 KB)


Additional file 3: Table S3: Links between smoking-related metabolites, enzymes and genes. The table describes the links showed in Figure 6 of the main text. The smoking-related metabolites, enzymes and genes are listed in the first and second columns. The score of interaction is given according to the definition by the Search Tool for the Retrieval of Interacting Genes/Proteins [1]. A reference for each link and a short description is provided. The Column of reaction shows the possible biochemical reaction of the corresponding link or the type of protein interaction. The enzymes includes, phospholipase A2, membrane associated (GIIC sPLA2), cytosolic phospholipase A2 (cPLA2), group 10 secretory phospholipase A2 (PLA2G10), lysophospholipase I (LYPLA1), apolipoprotein A-V (APOA5), uteroglobin (SCGB1A1), lecithin retinol acyltransferase (LRAT), nitric oxide synthase 1 (NOS1), solute carrier family 3 member 2 (SLC3A2), serine dehydratase (SDH), 3-hydroxybutyrate dehydrogenase, type 1 (BDH). The smoking-related gene/protein includes, S100 calcium binding protein A10 (S100A10), glypican 1 (GPC1), sulfatase 1 (SULF1), alcohol dehydrogenase 7 (ADH7), dehydrogenase member 3 (DHRS3), aldose reductase (AKR1B1), acetoacetyl-CoA synthetase (AACS), V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS), solute carrier family 7 (SLC7A11) and three enzyme listed above, PLA2G10, LYPLA1, SCGB1A1. The links in the network for male and female CS are combined and listed together. Smoking-related genes are show in italic. aa: diacyl-; ae: acyl-alkyl-; C0: carnitine; lysoPC: acyl-phosphatidylcholine; PC: phosphatidylcholine; SM (OH): hydroxysphingomyeline. (DOC 130 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Xu, T., Holzapfel, C., Dong, X. et al. Effects of smoking and smoking cessation on human serum metabolite profile: results from the KORA cohort study. BMC Med 11, 60 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: