Interaction of Cigarette Smoking and Polygenic Risk Score on Reduced Lung Function

The Gender Gap in the Relationship between Metabolic Syndrome and Restrictive Ventilatory Defects

BACKGROUND

Given the fundamental physiological differences between the sexes, this study aimed to investigate the effect of metabolic syndrome on ventilatory defects stratified by sex.

METHODS

We conducted a nationwide, pooled, cross-sectional study. Data from 45,788 participants (men, n = 15,859; women, n = 29,929) aged 30 years or more were obtained from the Taiwan Biobank. Age-sex-adjusted and multivariate logistic regression models were used to estimate the risk of developing impaired pulmonary function (restrictive or obstructive ventilatory defects) in individuals with or without metabolic syndromes. Separate models were also used to estimate the effect of metabolic syndrome scores and the effect of individual metabolic abnormalities on the risk of restrictive ventilatory defects.

RESULTS

The overall prevalence of metabolic syndrome was estimated to be 15.9% in Taiwan, much higher in men than in women (18.6% versus 14.4%). A significant association was observed between metabolic syndromes and the risk of restrictive ventilatory defects. The risk of developing a restrictive ventilator defect was 35% higher in participants with metabolic syndromes (odds ratio, 1.35; 95% confidence interval, 1.26-1.45) than in those without metabolic syndromes. Elevated blood pressure and a triglycerides abnormality were important predictors of restrictive ventilator defects. Sex-stratified subgroup analyses of the individual metabolic abnormalities indicated that men with abdominal obesity and women with dysglycemia were more likely to develop restrictive ventilatory defects.

CONCLUSIONS

Our study's evidence suggested that metabolic syndromes were important predictors of impaired pulmonary function and an increased risk of developing restrictive ventilatory defects, and its risk increased with increasing numbers of metabolic abnormalities.

Genetics of chronic respiratory disease

Chronic respiratory diseases, such as chronic obstructive pulmonary disease (COPD), asthma and interstitial lung diseases are frequently occurring disorders with a polygenic basis that account for a large global burden of morbidity and mortality. Recent large-scale genetic epidemiology studies have identified associations between genetic variation and individual respiratory diseases and linked specific genetic variants to quantitative traits related to lung function. These associations have improved our understanding of the genetic basis and mechanisms underlying common lung diseases. Moreover, examining the overlap between genetic associations of different respiratory conditions, along with evidence for gene-environment interactions, has yielded additional biological insights into affected molecular pathways. This genetic information could inform the assessment of respiratory disease risk and contribute to stratified treatment approaches.

circADAMTS6 via stabilizing CAMK2A is involved in smoking-induced emphysema through driving M2 macrophage polarization

Cigarette smoke (CS), an indoor environmental pollutant, is a prominent risk factor for emphysema, which is a pathological feature of chronic obstructive pulmonary disease (COPD). Emerging function of circRNAs in immune responses and disease progression shed new light to explore the pathogenesis of emphysema. In this research, we demonstrated, by single-cell RNA sequencing (scRNAseq), that the ratio of M2 macrophages were increased in lung tissues of humans and mice with smoking-related emphysema. Further, our data showed that circADAMTS6 was associated with cigarette smoke extract (CSE)-induced M2 macrophage polarization. Mechanistically, in macrophages, circADAMTS6 stabilized CAMK2A mRNA via forming a circADAMTS6/IGF2BP2/CAMK2A RNA-protein ternary complex to activate CREB, which drives M2 macrophage polarization and leads to emphysema. In addition, in macrophages of mouse lung tissues, downregulation of circADAMTS6 reversed M2 macrophage polarization, the proteinase/anti-proteinase imbalance, and the elastin degradation, which protecting against CS-induced emphysema. Moreover, for macrophages and in a model with co-cultured lung organoids, the target of circADAMTS6 restored the growth of lung organoids compared to CSE-treated macrophages. Our results also demonstrated that, for smokers and COPD smokers, elevation of circADAMTS6 negatively correlated with lung function. Overall, this study reveals a novel mechanism for circADAMTS6-driven M2 macrophage polarization in smoking-related emphysema and postulates that circADAMTS6 could serve as a diagnostic and therapeutic marker for smoking-related emphysema.

Association between n-3 PUFA and lung function: results from the NHANES 2007-2012 and Mendelian randomisation study

This study aimed to investigate the association between n-3 PUFA and lung function. First, a cross-sectional study was conducted based on the National Health and Nutrition Examination Survey (NHANES) 2007-2012 data. n-3 PUFA intake was obtained from 24-h dietary recalls. A multivariable linear regression model was used to assess the observational associations of n-3 PUFA intake with lung function. Subsequently, a two-sample Mendelian randomisation (MR) was performed to estimate the potential causal effect of n-3 PUFA on lung function. Genetic instrumental variables were extracted from published genome-wide association studies. Summary statistics about n-3 PUFA was from UK Biobank. Inverse variance weighted was the primary analysis approach. The observational study did not demonstrate a significant association between n-3 PUFA intake and most lung function measures; however, a notable exception was observed with significant findings in the highest quartile for forced vital capacity (FVC) and % predicted FVC. The MR results also showed no causal effect of circulating n-3 PUFA concentration on lung function (forced expiratory volume in one second (FEV1), β = 0·01301, se = 0·01932, P = 0·5006; FVC, β = -0·001894, se = 0·01704, P = 0·9115; FEV1:FVC, β = 0·03118, se = 0·01743, P = 0·07359). These findings indicate the need for further investigation into the impact of higher n-3 PUFA consumption on lung health.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function

BACKGROUND

Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function.

OBJECTIVES

We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions.

METHODS

We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings.

RESULTS

Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers.

CONCLUSION

Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Precision Approaches to Chronic Obstructive Pulmonary Disease Management

Chronic obstructive pulmonary disease (COPD) is a leading cause of morbidity and mortality worldwide. COPD heterogeneity has hampered progress in developing pharmacotherapies that affect disease progression. This issue can be addressed by precision medicine approaches, which focus on understanding an individual's disease risk, and tailoring management based on pathobiology, environmental exposures, and psychosocial issues. There is an urgent need to identify COPD patients at high risk for poor outcomes and to understand at a mechanistic level why certain individuals are at high risk. Genetics, omics, and network analytic techniques have started to dissect COPD heterogeneity and identify patients with specific pathobiology. Drug repurposing approaches based on biomarkers of specific inflammatory processes (i.e., type 2 inflammation) are promising. As larger data sets, additional omics, and new analytical approaches become available, there will be enormous opportunities to identify high-risk individuals and treat COPD patients based on their specific pathophysiological derangements. These approaches show great promise for risk stratification, early intervention, drug repurposing, and developing novel therapeutic approaches for COPD.

Association of Abnormal Lung Function and Its Subtypes With Arterial Stiffness: A Longitudinal Cohort Study

BACKGROUND

Prior studies have reported the cross-sectional relationship between lung function and arterial stiffness, while the longitudinal association remains unclear to date. This study aimed to investigate whether abnormal lung function and its subtypes at baseline are associated with increased arterial stiffness using a cohort.

METHODS AND RESULTS

This was a secondary analysis extracting 2461 participants from Beijing Health Management Cohort as baseline and annually followed for development of arterial stiffness. Abnormal lung function was defined by forced expiratory volume in 1s <80% of the predicted value, forced vital capacity of the predicted value, or forced expiratory volume in 1s/forced vital capacity ratio <70%. Increased arterial stiffness was determined by brachial-ankle pulse wave velocity ≥1400 cm/s. Cox proportional hazards model was used to calculate the hazard ratio and population attributable fraction. The mean age was 42.8±8.1 years, and 444 (18.0%) cases developed increased arterial stiffness during a median follow-up of 3.0 years. The adjusted hazard ratio (95% CI) of arterial stiffness was 1.47 (95% CI, 1.10-1.96) for abnormal lung function, with a population attributable fraction of 3.9% (95% CI, 0.8-7.1). Of subtypes, only obstructive ventilatory dysfunction was significantly associated with arterial stiffness (adjusted hazard ratio, 2.06 [95% CI, 1.27-3.36]), not restricted ventilatory dysfunction (adjusted hazard ratio, 0.95 [95% CI, 0.54-1.65]). Consistent results were observed on multiple sensitivity analyses.

CONCLUSIONS

Our study indicated a longitudinal association of abnormal lung function with increased arterial stiffness using a large cohort, especially for the obstructive ventilatory dysfunction.

Polygenic risk scores identify heterogeneity in asthma and chronic obstructive pulmonary disease

BACKGROUND

Asthma and chronic obstructive pulmonary disease (COPD) have distinct and overlapping genetic and clinical features.

OBJECTIVE

We sought to test the hypothesis that polygenic risk scores (PRSs) for asthma (PRSAsthma) and spirometry (FEV1 and FEV1/forced vital capacity; PRSspiro) would demonstrate differential associations with asthma, COPD, and asthma-COPD overlap (ACO).

METHODS

We developed and tested 2 asthma PRSs and applied the higher performing PRSAsthma and a previously published PRSspiro to research (Genetic Epidemiology of COPD study and Childhood Asthma Management Program, with spirometry) and electronic health record-based (Mass General Brigham Biobank and Genetic Epidemiology Research on Adult Health and Aging [GERA]) studies. We assessed the association of PRSs with COPD and asthma using modified random-effects and binary-effects meta-analyses, and ACO and asthma exacerbations in specific cohorts. Models were adjusted for confounders and genetic ancestry.

RESULTS

In meta-analyses of 102,477 participants, the PRSAsthma (odds ratio [OR] per SD, 1.16 [95% CI, 1.14-1.19]) and PRSspiro (OR per SD, 1.19 [95% CI, 1.17-1.22]) both predicted asthma, whereas the PRSspiro predicted COPD (OR per SD, 1.25 [95% CI, 1.21-1.30]). However, results differed by cohort. The PRSspiro was not associated with COPD in GERA and Mass General Brigham Biobank. In the Genetic Epidemiology of COPD study, the PRSAsthma (OR per SD: Whites, 1.3; African Americans, 1.2) and PRSspiro (OR per SD: Whites, 2.2; African Americans, 1.6) were both associated with ACO. In GERA, the PRSAsthma was associated with asthma exacerbations (OR, 1.18) in Whites; the PRSspiro was associated with asthma exacerbations in White, LatinX, and East Asian participants.

CONCLUSIONS

PRSs for asthma and spirometry are both associated with ACO and asthma exacerbations. Genetic prediction performance differs in research versus electronic health record-based cohorts.

Associations between polygenic risk scores and accelerated brain ageing in smokers

BACKGROUND

Smoking contributes to a variety of neurodegenerative diseases and neurobiological abnormalities, suggesting that smoking is associated with accelerated brain aging. However, the neurobiological mechanisms affected by smoking, and whether they are genetically influenced, remain to be investigated.

METHODS

Using structural magnetic resonance imaging data from the UK Biobank (n = 33 293), a brain age predictor was trained on non-smoking healthy groups and tested on smokers to obtain the BrainAge Gap (BAG). The cumulative effect of multiple common genetic variants associated with smoking was then calculated to acquire a polygenic risk score (PRS). The relationship between PRS, BAG, total gray matter volume (tGMV), and smoking parameters was explored and further genes included in the PRS were annotated to identify potential molecular mechanisms affected by smoking.

RESULTS

The BrainAge in smokers was predicted with very high accuracy (r = 0.725, MAE = 4.16). Smokers had a greater BAG (Cohen's d = 0.074, p < 0.0001) and higher PRS (Cohen's d = 0.63, p < 0.0001) than non-smokers. A higher PRS was associated with increased amount of smoking, mediated by BAG and tGMV. Several neurotransmitters and ion channel pathways were enriched in the group of smoking-related genes involved in addiction, brain synaptic plasticity, and some neurological disorders.

CONCLUSION

By using a simplified single indicator of the entire brain (BAG) in combination with the PRS, this study highlights the greater BAG in smokers and its linkage with genes and smoking behavior, providing insight into the neurobiological underpinnings and potential features of smoking-related aging.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Age-related loss of chromosome Y is associated with levels of sex hormone binding globulin and clonal hematopoiesis defined by TET2, TP53, and CBL mutations

Mosaic loss of the Y-chromosome (LOY) in peripheral blood leukocytes is the most common somatic alteration in men and linked to wide range of malignant and nonmalignant conditions. LOY is associated with age, smoking, and constitutional genetics. Here, we aimed to assess the relationships between LOY, serum biomarkers, and clonal hematopoiesis (CH). LOY in U.K. Biobank was strongly associated with levels of sex hormone binding globulin (SHBG), a key regulator of testosterone bioavailability. Mendelian randomization suggested a causal effect of SHBG on LOY but there was no evidence for an effect of LOY on SHBG. In contrast, age-related CH defined by somatic driver mutations was not associated with SHBG but was associated with LOY at clonal fractions above 30%. TET2, TP53, and CBL mutations were enriched in LOY cases, but JAK2 V617F was depleted. Our findings thus identify independent relationships between LOY, sex hormone levels, and CH.

Impaired pulmonary function and associated factors in the elderly with tuberculosis on admission: a preliminary report

BACKGROUND

Pulmonary tuberculosis (TB) can impair pulmonary function (PF), especially in the elderly. The risk factors associated with the severity of PF impairment in the elderly with pulmonary TB remain unclear. Hence, this retrospective study aimed to address this issue to help improve the management of TB in the elderly population.

METHODS

From January 2019 to February 2022, the elderly who were admitted to our hospital for pulmonary TB and underwent PF testing were included in this analysis. The forced expiratory volume in one second percent of predicted (FEV1% predicted) and clinical characteristics were collected and analyzed retrospectively. The extent of impaired PF was then categorized based on the FEV1% predicted and classified as grade 1-5. Logistic regression analysis was used to analyze the risk factors for impaired PF.

RESULTS

A total of 249 patients who met the enrollment criteria were included in this analysis. According to the results of FEV1% predicted, all patients were classified as grade 1 (n = 37), grade 2 (n = 46), grade 3 (n = 55), grade 4 (n = 56), or grade 5 (n = 55). Statistical analysis showed that albumin (adjusted odds ratio (aOR) = 0.928, P = 0.013), body mass index (BMI) < 18.5 kg/m² (aOR = 4.968, P = 0.046), lesion number ≥ 3 (aOR = 4.229, P < 0.001), male (aOR = 2.252, P = 0.009), respiratory disease (aOR = 1.669, P = 0.046), and cardiovascular disease (aOR = 2.489, P = 0.027) were related to the impairment of PF.

CONCLUSIONS

PF impairment is common in the elderly with pulmonary TB. The male sex, BMI < 18.5 kg/m², lesion number ≥ 3, hypoproteinemia, and respiratory and cardiovascular comorbidities were identified as risk factors for significant PF impairment. Our findings highlight the risk factors associated with PF impairment, which may be helpful to improve the current management of pulmonary TB in the elderly to save their lung function.

Inflammatory Bowel Diseases Before and After 1990

BACKGROUND AND AIMS

Inflammatory bowel disease (IBD) is caused by interaction of genetic and environmental risk factors. We evaluated potential determinants of the post-1990 increased incidence in North America.

METHODS

Using fitted generalized linear models, we assessed clinical features, smoking and genetic risk scores (GRS) for Crohn's disease (CD) and ulcerative colitis (UC) in the National Institutes of Diabetes, Digestion and Kidney Diseases IBD Genetics Consortium database, before and post 1990.

RESULTS

Among 2744 patients (55% CD, 42.2% UC), smoking status and GRS were the main determinants of diagnosis age. After 1990, smoking at diagnosis declined significantly in both UC and CD (34.1% vs 20.8%, P < .001, and 14.7% vs 8.7%, P = .06, respectively). In UC, ex-smoking increased (9% vs 15%, P < .001), and nonsmoking rates remained unchanged, whereas in CD, ex-smoking remained unchanged. CD-GRS and IBD-GRS were significantly associated with young diagnosis age, Jewish ethnicity, IBD family history, and surgery. CD-GRS showed a borderline significant decrease (P = .058) in multivariate analysis post 1990 but only when excluding surgery in the model; surgery significantly decreased post 1990 in both CD and UC. CD-GRS inversely correlated with smoking at diagnosis (P < .001) suggesting that, in the presence of smoking, CD may only require a low genetic risk to develop.

CONCLUSION

Significantly increase in ex-smoking correlates with UC incidence post 1990. Conversely, smoking risk decreased significantly post 1990 despite rising CD incidence. CD-GRS likewise trended to decrease post 1990 only when not accounting for a significant decrease in CD surgery. We therefore deduce that unaccounted risk factors (eg, dietary, obesity, antibiotic use, improved hygiene, etc.) or greater detection or presence of mild CD may underlie post-1990 increased CD incidence.

Polygenic scores for low lung function and the future risk of adverse health outcomes

Aims

Reduced lung function and adverse health outcomes are often observed. This study characterizes genetic susceptibility for reduced lung function and risk of developing a range of adverse health outcomes.

Methods

We studied 27,438 middle-aged adults from the Malmö Diet and Cancer study (MDCS), followed up to 28.8 years. Trait-specific Polygenic scores (PGS) for forced expiratory volume in 1 s (FEV₁) and forced vital capacity (FVC) were constructed for each participant using MDCS genetic data and summary statistics from the latest GWAS of lung function. Linear regression models and cox proportional hazards regression models were used to assess associations between adverse health outcomes and lung function-PGS.

Results

FEV₁-PGS and FVC-PGS were significantly associated with mean sBP at baseline after adjustments (FEV₁-PGS Q1 (highest PGS = highest lung function): 140.7mmHg vs. Q4: 141.5mmHg, p-value 0.008). A low FVC-PGS was significantly associated with the risk of future diabetic events after adjustments (Q4 vs. Q1 HR: 1.22 (CI 1.12–1.32), p-trend < 0.001) and had added value to risk prediction models for diabetes. Low FEV₁-PGS was significantly associated with future coronary events (Q4 vs. Q1 HR: 1.13 (CI: 1.04–1.22), p-trend 0.008). No significant association was found between PGS and sudden cardiac death, chronic kidney disease or all-cause mortality. Results remained largely unchanged in a subgroup of subjects when further adjusted for apolipoproteins.

Conclusion

Genetic susceptibility for reduced lung function is associated with higher sBP, increased risk of diabetes and to a lesser extent, future coronary events, suggesting etiological roles of lung function on these outcomes. Using PGS, high-risk groups could be early detected to implement early lifestyle changes to mitigate the risk.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12933-022-01661-y.

A robust and adaptive framework for interaction testing in quantitative traits between multiple genetic loci and exposure variables

The identification and understanding of gene-environment interactions can provide insights into the pathways and mechanisms underlying complex diseases. However, testing for gene-environment interaction remains a challenge since a.) statistical power is often limited and b.) modeling of environmental effects is nontrivial and such model misspecifications can lead to false positive interaction findings. To address the lack of statistical power, recent methods aim to identify interactions on an aggregated level using, for example, polygenic risk scores. While this strategy can increase the power to detect interactions, identifying contributing genes and pathways is difficult based on these relatively global results. Here, we propose RITSS (Robust Interaction Testing using Sample Splitting), a gene-environment interaction testing framework for quantitative traits that is based on sample splitting and robust test statistics. RITSS can incorporate sets of genetic variants and/or multiple environmental factors. Based on the user's choice of statistical/machine learning approaches, a screening step selects and combines potential interactions into scores with improved interpretability. In the testing step, the application of robust statistics minimizes the susceptibility to main effect misspecifications. Using extensive simulation studies, we demonstrate that RITSS controls the type 1 error rate in a wide range of scenarios, and we show how the screening strategy influences statistical power. In an application to lung function phenotypes and human height in the UK Biobank, RITSS identified highly significant interactions based on subcomponents of genetic risk scores. While the contributing single variant interaction signals are weak, our results indicate interaction patterns that result in strong aggregated effects, providing potential insights into underlying gene-environment interaction mechanisms.

Metabolic Differences between Ex-Smokers and Nonsmokers: A Metabolomic Analysis

The aim of this study was to compare changes in the metabolite levels of ex-smokers and nonsmokers using a metabolomics approach, accounting for the weight gain in ex-smokers. Volunteer ex-smokers and nonsmokers were recruited from two cohorts Shijingshan (174) and Xishan (78), respectively, at a 1 : 1 ratio for age and sex. Nontargeted metabolomics was performed on the volunteers' blood samples using liquid chromatography-mass spectrometry, and multivariate statistical analysis was performed using principal component analysis and orthogonal partial least squares discriminant analysis. Enrichment analysis was used to identify Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with differential metabolites and weighted gene co-expression network analysis and maximal correlation coefficient (MCC) algorithms were used to identify key metabolites. The results revealed no significant differences between the distribution of blood metabolite levels in the ex-smokers and nonsmokers. The biosynthesis of valine, leucine, and isoleucine was determined to be associated with differential metabolites, and five key metabolites were identified. Further analysis revealed differences in weight gain and regained metabolite levels in ex-smokers, and 10 differential metabolites were identified that may be associated with weight gain in ex-smokers. These findings suggest that quitting smoking restores metabolites to almost normal levels and results in weight gain. The identified key metabolites and metabolic pathways may also provide a basis for clinical studies.