Reversible and permanent effects of tobacco smoke exposure on airway epithelial gene expression

Smoking-associated gene expression alterations in nasal epithelium reveal immune impairment linked to lung cancer risk

BACKGROUND

Lung cancer is the leading cause of cancer-related death in the world. In contrast to many other cancers, a direct connection to modifiable lifestyle risk in the form of tobacco smoke has long been established. More than 50% of all smoking-related lung cancers occur in former smokers, 40% of which occur more than 15 years after smoking cessation. Despite extensive research, the molecular processes for persistent lung cancer risk remain unclear. We thus set out to examine whether risk stratification in the clinic and in the general population can be improved upon by the addition of genetic data and to explore the mechanisms of the persisting risk in former smokers.

METHODS

We analysed transcriptomic data from accessible airway tissues of 487 subjects, including healthy volunteers and clinic patients of different smoking statuses. We developed a computational model to assess smoking-associated gene expression changes and their reversibility after smoking is stopped, comparing healthy subjects to clinic patients with and without lung cancer.

RESULTS

We find persistent smoking-associated immune alterations to be a hallmark of the clinic patients. Integrating previous GWAS data using a transcriptional network approach, we demonstrate that the same immune- and interferon-related pathways are strongly enriched for genes linked to known genetic risk factors, demonstrating a causal relationship between immune alteration and lung cancer risk. Finally, we used accessible airway transcriptomic data to derive a non-invasive lung cancer risk classifier.

CONCLUSIONS

Our results provide initial evidence for germline-mediated personalized smoke injury response and risk in the general population, with potential implications for managing long-term lung cancer incidence and mortality.

Impact of smoking cessation duration on lung cancer mortality: A systematic review and meta-analysis

BACKGROUND

Smoking history is a heterogeneous situation for different populations, and numerous studies suggest that smoking cessation is conducive to reduce the mortality of lung cancer. However, no quantitative meta-analysis regarding smoking cessation duration based on different populations has demonstrated it clearly.

METHODS

We systematically searched four electronic databases (PubMed, Embase, the Cochrane Central Register of Controlled Trials, and Scoups) till February 2023. Eligible studies reported the association between lung cancer survival and duration of smoking cessation. Additionally, we stratified the study population according to whether they had lung cancer at the time they quit smoking. Studies were pooled with the random-effects model.

RESULTS

Out of the 11,361 potential studies initially identified, we included 24 studies involving 969,560 individuals in our analysis. Lung cancer mortality varied across two groups: general quitters and peri-diagnosis quitters. For general quitters, those who had quit smoking for less than 10 years exhibited an RR of 0.64 (95% CI [0.55-0.76]), while those who quit for 10-20 years had an RR of 0.33 (0.25-0.43), over 20 years had an RR of 0.16 (0.11-0.24), and never-smokers had an RR at 0.11 (0.07-0.15). Among peri-diagnosis quitters, the 1-year Overall Survival (OS) showed an RR of 0.80 (0.67-0.96), the 2-year OS had an RR of 0.89 (0.80-0.98), the 3-year OS had an RR of 0.93 (0.84-1.03), and the 5-year OS had an RR of 0.85 (0.76-0.96).

CONCLUSIONS

Earlier and longer smoking cessation is associated with reduced lung cancer mortality, no matter in which cessation stage for two different populations.

Metabolomics profiles alterations in cigarette smokers and heated tobacco product users

BACKGROUND

Heated tobacco products (HTPs) have gained global popularity, but their health risks remain unclear. Therefore, the current study aimed to identify plasma metabolites associated with smoking and HTP use in a large Japanese population to improve health risk assessment.

METHODS

Metabolomics data from 9,922 baseline participants of the Tsuruoka Metabolomics Cohort Study (TMCS) were analyzed to determine the association between smoking habits and plasma metabolites. Moreover, alterations in smoking-related metabolites among HTP users were examined based on data obtained from 3,334 participants involved from April 2018 to June 2019 in a follow-up survey.

RESULTS

Our study revealed that cigarette smokers had metabolomics profiles distinct from never smokers, with 22 polar metabolites identified as candidate biomarkers for smoking. These biomarker profiles of HTP users were closer to those of cigarette smokers than those of never smokers. The concentration of glutamate was higher in cigarette smokers, and biomarkers involved in glutamate metabolism were also associated with cigarette smoking and HTP use. Network pathway analysis showed that smoking was associated with the glutamate pathway, which could lead to endothelial dysfunction and atherosclerosis of the vessels.

CONCLUSIONS

Our study showed that the glutamate pathway is affected by habitual smoking. These changes in the glutamate pathway may partly explain the mechanism by which cigarette smoking causes cardiovascular disease. HTP use was also associated with glutamate metabolism, indicating that HTP use may contribute to the development of cardiovascular disease through mechanisms similar to those in cigarette use.

Serial Sampling of the Small Airway Epithelium to Identify Persistent Smoking-dysregulated Genes

Rationale: The small airway epithelium (beyond the sixth generation), the initiation site of smoking-induced airway disorders, is highly sensitive to the stress of smoking. Because of variations over time in smoking habits, the small airway epithelium transcriptome is dynamic, fluctuating not only among smokers but also within each smoker. Objectives: To perform accurate assessment of the smoking-related dysregulation of the human small airway epithelium despite the variation of smoking within the same individual and of the effects of smoking cessation on the dysregulated transcriptome. Methods: We conducted serial sampling of the same smokers and nonsmoker control subjects over time to identify persistent smoking dysregulation of the biology of the small airway epithelium over 1 year. We conducted serial sampling of smokers who quit smoking, before and after smoking cessation, to assess the effect of smoking cessation on the smoking-dysregulated genes. Measurements and Main Results: Repeated measures ANOVA of the small airway epithelium transcriptome sampled four times in the same individuals over 1 year enabled the identification of 475 persistent smoking-dysregulated genes. Most genes were normalized after 12 months of smoking cessation; however, 53 (11%) genes, including CYP1B1, PIR, ME1, and TRIM16, remained persistently abnormally expressed. Dysregulated pathways enriched with the nonreversible genes included xenobiotic metabolism signaling, bupropion degradation, and nicotine degradation. Conclusions: Analysis of repetitive sampling of the same individuals identified persistent smoking-induced dysregulation of the small airway epithelium transcriptome and the effect of smoking cessation. These results help identify targets for the development of therapies that can be applicable to smoking-related airway diseases.

Smoking during pregnancy is associated with the placental proteome

Maternal smoking during pregnancy (MSDP) is a significant risk factor for the development of foetal, neonatal, and childhood morbidities. We hypothesized that infants exposed to MSDP have a distinct proteomic expression in their term placentas compared to infants without such an exposure. A total of 39 infants exposed (cord blood cotinine levels of >1 ng/mL) and 44 infants not exposed to MSDP were included in the study. Women with chronic disease, body mass index of > 30, or a history of uterine surgery were excluded. Total proteome abundance was analysed with quantitative mass spectrometry. For univariate analysis of differences in placental protein levels between groups, ANOVA with multiple testing corrections by the Benjamini-Hochberg method was used. For multivariate analysis, we used principal component analysis, partial least squares, lasso, random forest, and neural networks. The univariate analyses showed four differentially abundant proteins (PXDN, CYP1A1, GPR183, and KRT81) when heavy and moderate smoking groups were compared to non-smokers. With the help of machine learning, we found that an additional six proteins (SEPTIN3, CRAT, NAAA, CD248, CADM3, and ZNF648) were discriminants of MSDP. The placental abundance of these ten proteins together explained 74.1% of the variation in cord blood cotinine levels (p = 0.002). Infants exposed to MSDP showed differential abundance of proteins in term placentas. We report differential placental abundance of several proteins for the first time in the setting of MSDP. We believe that these findings supplement the current understanding of how MSDP affects the placental proteome.

Smoking-related dysregulation of plasma circulating microRNAs: the Rotterdam study

BACKGROUND

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression. Differential miRNA expression, which is widely shown to be associated with the pathogenesis of various diseases, can be influenced by lifestyle factors, including smoking. This study aimed to investigate the plasma miRNA signature of smoking habits, the potential effect of smoking cessation on miRNA levels, and relate the findings with lung cancer incidence.

RESULTS

A targeted RNA-sequencing approach measured plasma miRNA levels in 2686 participants from the population-based Rotterdam study cohort. The association between cigarette smoking (current versus never) and 591 well-expressed miRNAs was assessed via adjusted linear regression models, identifying 41 smoking-associated miRNAs that passed the Bonferroni-corrected threshold (P < 0.05/591 = 8.46 × 10^-5). Moreover, we found 42 miRNAs with a significant association (P < 8.46 × 10^-5) between current (reference group) and former smokers. Then, we used adjusted linear regression models to explore the effect of smoking cessation time on miRNA expression levels. The expression levels of two miRNAs were significantly different within 5 years of cessation (P < 0.05/41 = 1.22 × 10^-3) from current smokers, while for cessation time between 5 and 15 years we found 19 miRNAs to be significantly different from current smokers, and finally, 38 miRNAs were significantly different after more than 15 years of cessation time (P < 1.22 × 10^-3). These results imply the reversibility of the smoking effect on plasma levels of at least 38 out of the 41 smoking-miRNAs following smoking cessation. Next, we found 8 out of the 41 smoking-related miRNAs to be nominally associated (P < 0.05) with the incidence of lung cancer.

CONCLUSIONS

This study demonstrates smoking-related dysregulation of plasma miRNAs, which might have a potential for reversibility when comparing different smoking cessation groups. The identified miRNAs are involved in several cancer-related pathways and include 8 miRNAs associated with lung cancer incidence. Our results may lay the groundwork for further investigation of miRNAs as potential mechanism linking smoking, gene expression and cancer.

A high-throughput cigarette smoke-treated bronchosphere model for disease-relevant phenotypic compound screening

Cigarette smoking (CS) is the leading cause of COPD, and identifying the pathways that are driving pathogenesis in the airway due to CS exposure can aid in the discovery of novel therapies for COPD. An additional barrier to the identification of key pathways that are involved in the CS-induced pathogenesis is the difficulty in building relevant and high throughput models that can recapitulate the phenotypic and transcriptomic changes associated with CS exposure. To identify these drivers, we have developed a cigarette smoke extract (CSE)-treated bronchosphere assay in 384-well plate format that exhibits CSE-induced decreases in size and increase in luminal secretion of MUC5AC. Transcriptomic changes in CSE-treated bronchospheres resemble changes that occur in human smokers both with and without COPD compared to healthy groups, indicating that this model can capture human smoking signature. To identify new targets, we ran a small molecule compound deck screening with diversity in target mechanisms of action and identified hit compounds that attenuated CSE induced changes, either decreasing spheroid size or increasing secreted mucus. This work provides insight into the utility of this bronchopshere model to examine human respiratory disease impacted by CSE exposure and the ability to screen for therapeutics to reverse the pathogenic changes caused by CSE.

DNA methylation changes from primary cultures through senescence-bypass in Syrian hamster fetal cells initially exposed to benzo[a]pyrene

Current chemical testing strategies are limited in their ability to detect non-genotoxic carcinogens (NGTxC). Epigenetic anomalies develop during carcinogenesis regardless of whether the molecular initiating event is associated with genotoxic (GTxC) or NGTxC events; therefore, epigenetic markers may be harnessed to develop new approach methodologies that improve the detection of both types of carcinogens. This study used Syrian hamster fetal cells to establish the chronology of carcinogen-induced DNA methylation changes from primary cells until senescence-bypass as an essential carcinogenic step. Cells exposed to solvent control for 7 days were compared to naïve primary cultures, to cells exposed for 7 days to benzo[a]pyrene, and to cells at the subsequent transformation stages: normal colonies, morphologically transformed colonies, senescence, senescence-bypass, and sustained proliferation in vitro. DNA methylation changes identified by reduced representation bisulphite sequencing were minimal at day-7. Profound DNA methylation changes arose during cellular senescence and some of these early differentially methylated regions (DMRs) were preserved through the final sustained proliferation stage. A set of these DMRs (e.g., Pou4f1, Aifm3, B3galnt2, Bhlhe22, Gja8, Klf17, and L1l) were validated by pyrosequencing and their reproducibility was confirmed across multiple clones obtained from a different laboratory. These DNA methylation changes could serve as biomarkers to enhance objectivity and mechanistic understanding of cell transformation and could be used to predict senescence-bypass and chemical carcinogenicity.

The impact of stress and lifestyle factors on short-term sickness absence in a large Danish industrial company

Objectives: Work-related stress and certain lifestyles have been associated with sickness absence among employees. We analysed relations between stress and lifestyles (as risk factors) and short-term sickness absence among employees of a large Danish industrial company. Moreover, we analysed the impact of risk factors on short-term sickness absence. Methods: A self-administered questionnaire was completed by 7753 employees (67% response rate). Relations between risk factors (dyssomnia, low physical activity, alcohol, smoking, unhealthy food habits, overweight and stress) and short-term sickness absence were analysed by calculating etiologic fractions and binary logistic regression. Relations between the number of risk factors and short-term sickness absence were also analysed. Finally, the cost of short-term sickness absence from risk factors was estimated to illustrate the potential savings in avoided loss of productivity that could be gained from intervention programmes. Results: Stress, overweight, smoking and dyssomnia are significantly related to short-term sickness absence. Etiologic fractions revealed that these factors were associated with between 29.8% and 37.8% of short-term sickness absence. The number of risk factors was also related to the risk and length of sickness absence. Conclusions: This study identified risk factors that could be addressed by intervention programmes to reduce short-term sickness absence. Based on the results, focus on the risk factors that account for most short-term sickness absence and reduction of the number of risk factors could potentially reduce short-term sickness absence. A 30% reduction is equivalent to an avoidance of a loss of productivity of 4.5%, corresponding to €9.4 million per year.

Predictors of five-repetition sit-to-stand test performance in patients with lumbar degenerative disease

BACKGROUND

The five-repetition sit-to-stand test (5R-STS) has recently been validated as an objective measure of functional impairment in patients with lumbar degenerative disease (LDD). Knowledge of factors influencing 5R-STS performance is useful to correct for confounders, create personalized adjusted test times, and potentially identify prognostic subgroups. We evaluate factors predicting the 5R-STS performance in patients with LDD.

METHODS

Patients with LDD requiring surgery were included. Each participant performed the 5R-STS and completed a questionnaire that included their age, gender, weight, height, body mass index (BMI), smoking status, education level, employment type, ability to work, analgesic drug usage, history of previous spinal surgery, and EQ5D depression and anxiety domain. Surgical indication and index level of the spinal pathology were also recorded. Predictors of 5R-STS were identified through multivariable linear regression.

RESULTS

The cohort consisted of 240 patients, 47.9% being female (mean age, 47.7 ± 13.6 years). In the final multivariable model incorporating confounders, height (regression coefficient (RC), 0.08; 95% confidence interval (CI), 0.003/0.16, p = 0.042) and being an active smoker (RC, 2.44; 95%CI, 0.56/4.32, p = 0.012) were significant predictors of worse 5R-STS performance. Full ability to work (RC, - 2.39; 95%CI, - 4.39/ - 0.39, p = 0.020) was associated with a better 5R-STS performance. Age, height, surgical indication, index level of pathology, history of previous spine surgery, history of pain, analgesic drug use, employment type, and severity of anxiety and depression symptoms demonstrated confounding effect on the 5R-STS time.

CONCLUSIONS

Greater height, being an active smoker, and inability to work are significant predictors of worse 5R-STS performance in patients with LDD.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT03303300 and NCT03321357.

Smoking modulates different secretory subpopulations expressing SARS-CoV-2 entry genes in the nasal and bronchial airways

SARS-CoV-2 infection and disease severity are influenced by viral entry (VE) gene expression patterns in the airway epithelium. The similarities and differences of VE gene expression (ACE2, TMPRSS2, and CTSL) across nasal and bronchial compartments have not been fully characterized using matched samples from large cohorts. Gene expression data from 793 nasal and 1673 bronchial brushes obtained from individuals participating in lung cancer screening or diagnostic workup revealed that smoking status (current versus former) was the only clinical factor significantly and reproducibly associated with VE gene expression. The expression of ACE2 and TMPRSS2 was higher in smokers in the bronchus but not in the nose. scRNA-seq of nasal brushings indicated that ACE2 co-expressed genes were highly expressed in club and C15orf48+ secretory cells while TMPRSS2 co-expressed genes were highly expressed in keratinizing epithelial cells. In contrast, these ACE2 and TMPRSS2 modules were highly expressed in goblet cells in scRNA-seq from bronchial brushings. Cell-type deconvolution of the gene expression data confirmed that smoking increased the abundance of several secretory cell populations in the bronchus, but only goblet cells in the nose. The association of ACE2 and TMPRSS2 with smoking in the bronchus is due to their high expression in goblet cells which increase in abundance in current smoker airways. In contrast, in the nose, these genes are not predominantly expressed in cell populations modulated by smoking. In individuals with elevated lung cancer risk, smoking-induced VE gene expression changes in the nose likely have minimal impact on SARS-CoV-2 infection, but in the bronchus, smoking may lead to higher viral loads and more severe disease.

Early transcriptional responses of bronchial epithelial cells to whole cigarette smoke mirror those of in-vivo exposed human bronchial mucosa

BACKGROUND

Despite the well-known detrimental effects of cigarette smoke (CS), little is known about the complex gene expression dynamics in the early stages after exposure. This study aims to investigate early transcriptomic responses following CS exposure of airway epithelial cells in culture and compare these to those found in human CS exposure studies.

METHODS

Primary bronchial epithelial cells (PBEC) were differentiated at the air-liquid interface (ALI) and exposed to whole CS. Bulk RNA-sequencing was performed at 1 h, 4 h, and 24 h hereafter, followed by differential gene expression analysis. Results were additionally compared to data retrieved from human CS studies.

RESULTS

ALI-PBEC gene expression in response to CS was most significantly changed at 4 h after exposure. Early transcriptomic changes (1 h, 4 h post CS exposure) were related to oxidative stress, xenobiotic metabolism, higher expression of immediate early genes and pro-inflammatory pathways (i.e., Nrf2, AP-1, AhR). At 24 h, ferroptosis-associated genes were significantly increased, whereas PRKN, involved in removing dysfunctional mitochondria, was downregulated. Importantly, the transcriptome dynamics of the current study mirrored in-vivo human studies of acute CS exposure, chronic smokers, and inversely mirrored smoking cessation.

CONCLUSION

These findings show that early after CS exposure xenobiotic metabolism and pro-inflammatory pathways were activated, followed by activation of the ferroptosis-related cell death pathway. Moreover, significant overlap between these transcriptomic responses in the in-vitro model and human in-vivo studies was found, with an early response of ciliated cells. These results provide validation for the use of ALI-PBEC cultures to study the human lung epithelial response to inhaled toxicants.

Nicotine dose-dependent epigenomic-wide DNA methylation changes in the mice with long-term electronic cigarette exposure

Epigenomic-wide DNA methylation profiling holds the potential to reflect both electronic cigarette exposure-associated risks and individual poor health outcomes. However, a systemic study in animals or humans is still lacking. Using the Infinium Mouse Methylation BeadChip, we examined the DNA methylation status of white blood cells in male ApoE-/- mice after 14 weeks of electronic cigarette exposure with the InExpose system (2 hr/day, 5 days/week, 50% PG and 50% VG) with low (6 mg/ml) and high (36 mg/ml) nicotine concentrations. Our results indicate that electronic cigarette aerosol inhalation induces significant alteration of 8,985 CpGs in a dose-dependent manner (FDR<0.05); 7,389 (82.2%) of the CpG sites are annotated with known genes. Among the top 6 significant CpG sites (P-value<1e-8), 4 CpG sites are located in the known genes, and most (3/5) of these genes have been related to cigarette smoking. The other two CpGs are close to/associated with the Phc2 gene that was recently linked to smoking in a transcriptome-wide associations study. Furthermore, the gene set enrichment analysis highlights the activation of MAPK and 4 cardiomyocyte/cardiomyopathy-related signaling pathways (including adrenergic signaling in cardiomyocytes and arrhythmogenic right ventricular cardiomyopathy) following repeated electronic cigarette use. The MAPK pathway activation correlates well with our finding of increased cytokine mRNA expression after electronic cigarette exposure in the same batch of mice. Interestingly, two pathways related to mitochondrial activities, namely mitochondrial gene expression and mitochondrial translation, are also activated after electronic cigarette exposure. Elucidating the relationship between these pathways and the increased circulating mitochondrial DNA observed here will provide further insight into the cell-damaging effects of prolonged inhalation of e-cigarette aerosols.

Evaluating the Influences of Confounding Variables on Benchmark Dose using a Case Study in the Field of Ionizing Radiation

Purpose A vast amount of data regarding the effects of radiation stressors on transcriptional changes has been produced over the past few decades. These data have shown remarkable consistency across platforms and experimental design, enabling increased understanding of early molecular effects of radiation exposure. However, the value of transcriptomic data in the context of risk assessment is not clear and represents a gap that is worthy of further consideration. Recently, benchmark dose (BMD) modeling has shown promise in correlating a transcriptional point of departure (POD) to that derived using phenotypic outcomes relevant to human health risk assessment. Although frequently applied in chemical toxicity evaluation, our group has recently demonstrated application within the field of radiation research. This approach allows the possibility to quantitatively compare radiation-induced gene and pathway alterations across various datasets using BMD values and derive meaningful biological effects. However, before BMD modeling can confidently be used, an understanding of the impact of confounding variables on BMD outputs is needed. Methods: To this end, BMD modeling was applied to a publicly available microarray dataset (Gene Expression Omnibus #GSE23515) that used peripheral blood ex-vivo gamma-irradiated at 0.82 Gy/min, at doses of 0, 0.1, 0.5 or 2 Gy, and assessed 6 hours post-exposure. The dataset comprised six female smokers (F-S), six female non-smokers (F-NS), six male smokers (M-S), and six male non-smokers (M-NS). Results: A combined total of 412 genes were fit to models and the BMD distribution was noted to be bi-modal across the four groups. A total of 74, 41, 62 and 62 genes were unique to the F-NS, M-NS, F-S and M-S groups. Sixty-two BMD modeled genes and nine pathways were common across all four groups. There were no differential sensitivity of responses in the robust common genes and pathways. Conclusion: For radiation-responsive genes and pathways common across the study groups, the BMD distribution of transcriptional activity was unaltered by sex and smoking status. Although further validation of the data is needed, these initial findings suggest BMD values for radiation relevant genes and pathways are robust and could be explored further in future studies.

Nicotine-induced miR-21-3p promotes chemoresistance in lung cancer by negatively regulating FOXO3a

Lung cancer is the leading cause of cancer-related mortality worldwide and cigarette smoking is reported to contribute to the lung cancer-related mortality. The present study aimed to investigate the molecular mechanism underlying nicotine-induced chemoresistance in lung cancer. The expression of microRNA (miR)-21-3p and its predicted target FOXO3a in lung cancer cells was detected via reverse transcription-quantitative PCR, in the presence or absence of nicotine. The regulatory effect of miR-21-3p and FOXO3a on lung cancer cell proliferation and apoptosis induced by docetaxel or cisplatin treatment was evaluated by performing Cell Counting Kit-8 and Annexin V/PI staining assays, respectively. The interaction between miR-21-3p and FOXO3a was analyzed by performing luciferase reporter assays and western blotting. FOXO3a overexpression rescue experiments were conducted in vitro and in vivo using a xenograft mouse model to assess the function of miR-21-3p/FOXO3a in lung cancer. Nicotine induced miR-21-3p expression in lung cancer cells in a dose-dependent manner. miR-21-3p downregulated FOXO3a expression by directly binding to the 3′-untranslated region of FOXO3a. Moreover, miR-21-3p knockdown sensitized lung cancer cells to docetaxel or cisplatin treatment. Mechanistically, FOXO3a was predicted as a direct target of miR-21-3p. FOXO3a overexpression promoted the chemosensitivity of lung cancer cells to docetaxel or cisplatin treatment. Furthermore, FOXO3a overexpression antagonized the regulatory function of miR-21-3p on docetaxel- or cisplatin-treated lung cancer cells. In the docetaxel- or cisplatin-treated lung cancer xenograft mouse model, miR-21-3p promoted chemoresistance via negatively regulating FOXO3a. Therefore, the present study demonstrated that nicotine-induced miR-21-3p promoted chemoresistance to docetaxel or cisplatin treatment via negatively regulating FOXO3a, which may serve as a novel therapeutic strategy for the treatment of patients with chemoresistant lung cancer.

HIMA2: high-dimensional mediation analysis and its application in epigenome-wide DNA methylation data

Mediation analysis plays a major role in identifying significant mediators in the pathway between environmental exposures and health outcomes. With advanced data collection technology for large-scale studies, there has been growing research interest in developing methodology for high-dimensional mediation analysis. In this paper we present HIMA2, an extension of the HIMA method (Zhang in Bioinformatics 32:3150-3154, 2016). First, the proposed HIMA2 reduces the dimension of mediators to a manageable level based on the sure independence screening (SIS) method (Fan in J R Stat Soc Ser B 70:849-911, 2008). Second, a de-biased Lasso procedure is implemented for estimating regression parameters. Third, we use a multiple-testing procedure to accurately control the false discovery rate (FDR) when testing high-dimensional mediation hypotheses. We demonstrate its practical performance using Monte Carlo simulation studies and apply our method to identify DNA methylation markers which mediate the pathway from smoking to reduced lung function in the Coronary Artery Risk Development in Young Adults (CARDIA) Study.

Clinical Study of Aspirin and Zileuton on Biomarkers of Tobacco-Related Carcinogenesis in Current Smokers

The chemopreventive effect of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) on lung cancer risk is supported by epidemiologic and preclinical studies. Zileuton, a 5-lipoxygenase inhibitor, has additive activity with NSAIDs against tobacco carcinogenesis in preclinical models. We hypothesized that cyclooxygenase plus 5-lipoxygenase inhibition would be more effective than a placebo in modulating the nasal epithelium gene signatures of tobacco exposure and lung cancer. We conducted a randomized, double-blinded study of low-dose aspirin plus zileuton vs. double placebo in current smokers to compare the modulating effects on nasal gene expression and arachidonic acid metabolism. In total, 63 participants took aspirin 81 mg daily plus zileuton (Zyflo CR) 600 mg BID or the placebo for 12 weeks. Nasal brushes from the baseline, end-of-intervention, and one-week post intervention were profiled via microarray. Aspirin plus zilueton had minimal effects on the modulation of the nasal or bronchial gene expression signatures of smoking, lung cancer, and COPD but favorably modulated a bronchial gene expression signature of squamous dysplasia. Aspirin plus zileuton suppressed urinary leukotriene but not prostaglandin E2, suggesting shunting through the cyclooxygenase pathway when combined with 5-lipoxygenase inhibition. Continued investigation of leukotriene inhibitors is needed to confirm these findings, understand the long-term effects on the airway epithelium, and identify the safest, optimally dosed agents.

Five-Repetition Sit-to-Stand Test Performance in Healthy Individuals: Reference Values and Predictors From 2 Prospective Cohorts

Objective

The 5-repetition-sit-to-stand (5R-STS) test is an objective test of functional impairment- commonly used in various diseases, including lumbar degenerative disc diseases. It is used to measure the severity of disease and to monitor recovery. We aimed to evaluate reference values for the test, as well as factors predicting 5R-STS performance in healthy adults.

Methods

Healthy adults (> 18 years of age) were recruited, and their 5R-STS time was measured. Their age, sex, weight, height, body mass index (BMI), smoking status, education level, work situation and EuroQOL-5D Healthy & Anxiety category were recorded. Linear regression analysis was employed to identify predictors of 5R-STS performance.

Results

We included 172 individuals with mean age of 39.4±14.1 years and mean BMI of 24.0 ±4.0 kg/m². Females constituted 57%. Average 5R-STS time was 6.21 ±1.92 seconds, with an upper limit of normal of 12.39 seconds. In a multivariable model, age (regression coefficient [RC], 0.07; 95% confidence interval [CI], 0.05/0.09; p<0.001), male sex (RC, -0.87; 95% CI, -1.50 to -0.23; p=0.008), BMI (RC, 0.40; 95% CI, 0.10–0.71; p=0.010), height (RC, 0.13; 95% CI, 0.04–0.22; p=0.006), and houseworker status (RC, -1.62; 95% CI, -2.93 to -0.32; p=0.016) were significantly associated with 5R-STS time. Anxiety and depression did not influence performance significantly (RC, 0.82; 95% CI, -0.14 to 1.77; p=0.097).

Conclusion

The presented reference values can be applied as normative data for 5R-STS in healthy adults, and are necessary to judge what constitutes abnormal performance. We identified several significant factors associated with 5R-STS performance that may be used to calculate individualized expected test times.

Pathogenesis of COVID-19 described through the lens of an undersulfated and degraded epithelial and endothelial glycocalyx

The glycocalyx surrounds every eukaryotic cell and is a complex mesh of proteins and carbohydrates. It consists of proteoglycans with glycosaminoglycan side chains, which are highly sulfated under normal physiological conditions. The degree of sulfation and the position of the sulfate groups mainly determine biological function. The intact highly sulfated glycocalyx of the epithelium may repel severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) through electrostatic forces. However, if the glycocalyx is undersulfated and 3-O-sulfotransferase 3B (3OST-3B) is overexpressed, as is the case during chronic inflammatory conditions, SARS-CoV-2 entry may be facilitated by the glycocalyx. The degree of sulfation and position of the sulfate groups will also affect functions such as immune modulation, the inflammatory response, vascular permeability and tone, coagulation, mediation of sheer stress, and protection against oxidative stress. The rate-limiting factor to sulfation is the availability of inorganic sulfate. Various genetic and epigenetic factors will affect sulfur metabolism and inorganic sulfate availability, such as various dietary factors, and exposure to drugs, environmental toxins, and biotoxins, which will deplete inorganic sulfate. The role that undersulfation plays in the various comorbid conditions that predispose to coronavirus disease 2019 (COVID-19), is also considered. The undersulfated glycocalyx may not only increase susceptibility to SARS-CoV-2 infection, but would also result in a hyperinflammatory response, vascular permeability, and shedding of the glycocalyx components, giving rise to a procoagulant and antifibrinolytic state and eventual multiple organ failure. These symptoms relate to a diagnosis of systemic septic shock seen in almost all COVID-19 deaths. The focus of prevention and treatment protocols proposed is the preservation of epithelial and endothelial glycocalyx integrity.

Utilizing patient information to identify subtype heterogeneity of cancer driver genes

Identifying cancer driver genes is essential for understanding the mechanisms of carcinogenesis and designing therapeutic strategies. Although driver genes have been identified for many cancer types, it is still not clear whether the selection pressure of driver genes is homogeneous across cancer subtypes. We propose a statistical framework MutScot to improve the identification of driver genes and to investigate the heterogeneity of driver genes across cancer subtypes. Through simulation studies, we show that MutScot properly controls the type I error in detecting driver genes. In addition, we demonstrate that MutScot can identify subtype heterogeneity of driver genes. Applications to three studies in The Cancer Genome Atlas (TCGA) project showcase that MutScot has a desirable sensitivity for detecting driver genes and that MutScot identifies subtype heterogeneity of driver genes in breast cancer and lung cancer with regards to the status of hormone receptor and to the smoking status, respectively.

Impact of Adverse Childhood Experiences on Oral Health Among Women in the United States: Findings From the Behavioral Risk Factor Surveillance System

It is well established that adverse childhood experiences (ACEs) negatively impact physical and mental health. There is a paucity of research examining the impact of ACEs on oral health. Therefore, the aim of this study was to better understand how ACEs may impact oral health utilizing a nationally representative sample of females. Data from the 2010 Behavioral Risk Factor Surveillance System (BRFSS) were utilized. The sample comprised 36,249 females. Chi-square tests were used to examine associations between multiple forms of ACEs and oral health, mental health, smoking, and health-related quality of life. Logistic regression models were generated for the outcomes of interest: last dental cleaning, significant tooth loss, smoking, and health-related quality of life. Post-hoc analyses were conducted to examine whether smoking and dental cleaning moderated the relationship between ACEs and significant tooth loss. All ACE categories were significantly associated with poor health-related quality of life, being a current or former smoker, last dental cleaning >1 year, and significant tooth loss even after adjusting for demographic variables. Last dental cleaning and being a current or former smoker moderated the relationship between ACEs and significant tooth loss. This study affirms other findings that women with a history of ACEs were more likely to engage in health-risk behaviors, specifically cigarette smoking, which has been strongly associated with ACEs. These findings underscore the importance of including oral health in future studies given it has implications for disease development. Furthermore, dental visits may provide a window of opportunity for trauma-informed care.

Cigarette smoking-associated isoform switching and 3' UTR lengthening via alternative polyadenylation

Cigarette smoking induces a profound transcriptomic and systemic inflammatory response. Previous studies have focused on gene level differential expression of smoking, but the genome-wide effects of smoking on alternative isoform regulation have not yet been described. We conducted RNA sequencing in whole-blood samples of 454 current and 767 former smokers in the COPDGene Study, and we analyzed the effects of smoking on differential usage of isoforms and exons. At 10% FDR, we detected 3167 differentially expressed genes, 945 differentially used isoforms and 160 differentially used exons. Isoform switch analysis revealed widespread 3' UTR lengthening associated with cigarette smoking. The lengthening of these 3' UTRs was consistent with alternative usage of distal polyadenylation sites, and these extended 3' UTR regions were significantly enriched with functional sequence elements including microRNA and RNA-protein binding sites. These findings warrant further studies on alternative polyadenylation events as potential biomarkers and novel therapeutic targets for smoking-related diseases.

Validation of in vitro models for smoke exposure of primary human bronchial epithelial cells

RATIONALE

The bronchial epithelium is constantly challenged by inhalative insults including cigarette smoke (CS), a key risk factor for lung disease. In vitro exposure of bronchial epithelial cells using CS extract (CSE) is a widespread alternative to whole CS (wCS) exposure. However, CSE exposure protocols vary considerably between studies, precluding direct comparison of applied doses. Moreover, they are rarely validated in terms of physiological response in vivo and the relevance of the findings is often unclear.

METHODS

We tested six different exposure settings in primary human bronchial epithelial cells (phBECs), including five CSE protocols in comparison with wCS exposure. We quantified cell-delivered dose and directly compared all exposures using expression analysis of 10 well-established smoke-induced genes in bronchial epithelial cells. CSE exposure of phBECs was varied in terms of differentiation state, exposure route, duration of exposure, and dose. Gene expression was assessed by quantitative Real-Time PCR (qPCR) and Western Blot analysis. Cell type-specific expression of smoke-induced genes was analyzed by immunofluorescent analysis.

RESULTS

Three surprisingly dissimilar exposure types, namely chronic CSE treatment of differentiating phBECs, acute CSE treatment of submerged basal phBECs, and wCS exposure of differentiated phBECs performed best, resulting in significant upregulation of seven (chronic CSE) and six (acute wCS, acute submerged CSE exposure) out of 10 genes. Acute apical or basolateral exposure of differentiated phBECs with CSE was much less effective despite similar doses used.

CONCLUSIONS

Our findings provide guidance for the design of human in vitro CS exposure models in experimental and translational lung research.

Chronic E-Cigarette Exposure Alters Human Alveolar Macrophage Morphology and Gene Expression

INTRODUCTION

Alveolar macrophages (AMs) are lung-resident immune cells that phagocytose inhaled particles and pathogens, and help coordinate the lung's immune response to infection. Little is known about the impact of chronic e-cigarette use (i.e. vaping) on this important pulmonary cell type. Thus, we determined the effect of vaping on alveolar macrophage (AM) phenotype and gene expression.

METHODS

We recruited never-smokers, smokers, and e-cigarette users (vapers) and performed research bronchoscopies to isolate AMs from bronchoalveolar lavage fluid (BALF) samples and epithelial cells from bronchial brushings. We then performed morphological analyses and used the Nanostring platform to look for changes in gene expression.

RESULTS

AMs obtained from smokers and vapers were phenotypically distinct from those obtained from non-smokers, and from each other. Immunocytochemistry revealed that vapers AM's had significantly elevated inducible nitric oxide synthase (iNOS; M1) expression and significantly reduced CD301a (M2) expression compared to non-smokers or smokers. Vapers' AMs and bronchial epithelia exhibited unique changes in gene expression compared to non-smokers or smokers. Moreover, vapers' AMs were the most affected of all groups and had 124 genes uniquely downregulated. Gene ontology analysis revealed that vapers and smokers had opposing changes in biological processes.

CONCLUSIONS

These data indicate that vaping causes unique changes to AMs and bronchial epithelia compared to non-smokers and smokers which may impact pulmonary host defense.

IMPLICATIONS

These data indicate that normal "healthy" vapers have altered AMs and may be at risk of developing abnormal immune responses to inflammatory stimuli.

Smoking induces sex-specific changes in the small airway proteome

Introduction

Cigarette smoke triggers many cellular and signaling responses in the lung and the resulting inflammation plays a central role in smoke-related lung diseases, such as COPD. We explored the effects of smoking on the small airway proteome in samples obtained by collection of exhaled particles with the aim to identify specific proteins dysregulated by smoking.

Methods

Exhaled particles were obtained from 38 current smokers, 47 former smokers and 22 healthy controls with the PExA method. 120 ng of sample was collected from individual subjects and analyzed with the SOMAscan proteomics platform. General linear model-based statistics were performed.

Results

Two hundred and three proteins were detected in at least half of 107 total samples. Active smoking exerted a significant impact on the protein composition of respiratory tract lining fluid (RTLF), with 81 proteins altered in current smokers compared to never smokers (p < 0.05, q < 0.124). Among the proteins most clearly discriminating between current and never smokers were sRAGE, FSTL3, SPOCK2 and protein S, all of them being less abundant in current smokers. Analysis stratified for sex unveiled sex differences with more pronounced proteomic alterations due to active smoking in females than males. Proteins whose abundance was altered by active smoking in women were to a larger extent related to the complement system. The small airway protein profile of former smokers appeared to be more similar to that observed in never smokers.

Conclusions

The study shows that smoking has a strong impact on protein expression in the small airways, and that smoking affects men and women differently, suggesting PExA sampling combined with high sensitivity protein analysis offers a promising platform for early detection of COPD and identification of novel COPD drug targets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12931-021-01825-6.

Assessing the consequences of environmental exposures on the expression of the human receptor and proteases involved in SARS-CoV-2 cell-entry

The role of environmental condition on the infection by the novel pathogenic SARS-CoV-2 virus remains uncertain. In here, exploiting a large panel of publicly available genome-wide data, we investigated whether the human receptor ACE2 and human proteases TMPRSS2, FURIN and CATHEPSINs (B, L and V), which are involved in SARS-CoV-2 cell entry, are transcriptionally regulated by environmental cues. We report that more than 50 chemicals modulate the expression of ACE2 or human proteases important for SARS-CoV-2 cell entry. We further demonstrate that transcription factor AhR, which is commonly activated by pollutants, binds to the promoter of TMPRSS2 and enhancers and/or promoters of Cathepsin B, L and V encoding genes. Our exploratory study documents an influence of environmental exposures on the expression of genes involved in SARS-CoV-2 cell entry. These results could be conceptually and medically relevant to our understanding of the COVID-19 disease, and should be further explored in laboratory and epidemiologic studies.

Flow-Cytometric Analysis and Purification of Airway Epithelial-Cell Subsets

The human airway epithelium is essential in homeostasis, and epithelial dysfunction contributes to chronic airway disease. Development of flow-cytometric methods to characterize subsets of airway epithelial cells will enable further dissection of airway epithelial biology. Leveraging single-cell RNA-sequencing data in combination with known cell type-specific markers, we developed panels of antibodies to characterize and isolate the major airway epithelial subsets (basal, ciliated, and secretory cells) from human bronchial epithelial-cell cultures. We also identified molecularly distinct subpopulations of secretory cells and demonstrated cell subset-specific expression of low-abundance transcripts and microRNAs that are challenging to analyze with current single-cell RNA-sequencing methods. These new tools will be valuable for analyzing and separating airway epithelial subsets and interrogating airway epithelial biology.

Prospective study of plasma levels of coenzyme Q10 and lung cancer risk in a low-income population in the Southeastern United States

Background

Coenzyme Q10 (CoQ10) is a ubiquitous molecule in living organisms serving as a cofactor in energy production. Epidemiological studies have reported low CoQ10 levels being associated with an increased risk of various cancers. We conducted the first study to evaluate the association of CoQ10 concentrations with lung cancer risk.

Methods

A nested case‐control study including 201 lung cancer cases and 395 matched controls from the Southern Community Cohort Study was conducted. Plasma CoQ10 levels were measured using high‐performance liquid chromatography with photo‐diode array detection. Conditional logistic regression models were applied to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association between plasma CoQ10 levels and lung cancer risk.

Results

Plasma CoQ10 concentration was inversely associated with the risk of lung cancer. After adjusting for age, sex, race, and socioeconomic status, the OR (95% CI) comparing the third to first tertile was 0.57 (0.36–0.91, P for trend = 0.02). Further adjustments for smoking, alcohol, chronic obstructive pulmonary disease, and body mass index attenuated the point estimate slightly (OR = 0.60, 95% CI = 0.34–1.08, P for trend = 0.11), comparing third to first tertiles. Stratified analyses identified a significant inverse association between plasma CoQ10 levels and lung cancer risk in current smokers, but not in former/never smokers. The association was more evident in cases who were diagnosed within 1 year of blood draw than in cases diagnosed after 1 year.

Conclusions

Low plasma CoQ10 was significantly associated with increased lung cancer risk, particularly among current smokers. The stronger association seen shortly following the blood draw suggests that CoQ10 may be related to disease progression.

Association of Cigarette Smoking, COPD, and Lung Cancer With Expression of SARS-CoV-2 Entry Genes in Human Airway Epithelial Cells

SARS-CoV-2 enters into human airway epithelial cells via membrane fusion or endocytosis, and this process is dependent on ACE2, TMPRSS2, and cathepsin L. In this study, we examined the expression profiles of the three SARS-CoV-2 entry genes in primary human airway epithelial cells isolated from smokers, non-smokers, patients with chronic obstructive pulmonary disease or lung cancer. An exhaustive search of the GEO database was performed to identify eligible data on 1st June 2020. In total, 46 GEO datasets comprising transcriptomic data of 3,053 samples were identified as eligible data for further analysis. All meta-analysis were performed using RStudio. Standardized mean difference was utilized to assess the effect size of a factor on the expression of targeted genes and 95% confidence intervals (CIs) were calculated. This study revealed that (i) cigarette smoking is associated with an increased expression of ACE2 and TMPRSS2 and a decreased expression of cathepsin L; (ii) significant alternations in expression of ACE2, TMPRSS2, and cathepsin L were observed between current smokers and former smokers, but not between former smokers and never smokers; (iii) when compared with healthy controls with identical smoking status, patients with COPD or lung cancer showed negligible changes in expression of ACE2, TMPRSS2, and cathepsin L. Therefore, this study implicates cigarette smoking might contribute to the development of COVID-19 by affecting the expression of SARS-CoV-2 entry genes, while smoking cessation could be effective to reduce the potential risk.

Smoking and COVID-19: Similar bronchial ACE2 and TMPRSS2 expression and higher TMPRSS4 expression in current versus never smokers

Uncertainties remain concerning the pathophysiology, epidemiology, and potential therapeutics for COVID-19. Among unsettled controversies is whether tobacco smoking increases or protects from severe COVID-19. Several epidemiological studies reported reduced COVID-19 hospitalizations among smokers, while other studies reported the opposite trend. Some authors assumed that smokers have elevated airway expression of ACE2, the cell recognition site of the SARS-Cov-2 spike protein, but this suggestion remains unverified. We therefore performed data mining of two independent NCBI GEO genome-wide RNA expression files (GSE7894 and GSE994) and report that in both data sets, current smokers and never smokers have, on average, closely similar bronchial epithelial cell mRNA levels of ACE2, as well as TMPRSS2, coding for a serine protease priming SARS-Cov-2 for cell entry, and ADAM17, coding for a protease implicated in ACE2 membrane shedding. In contrast, the expression levels of TMPRSS4, coding for a protease that primes SARS-CoV-2 for cell entry similarly to TMPRSS2, were elevated in bronchial epithelial cells from current smokers compared with never smokers, suggesting that higher bronchial TMPRSS4 levels in smokers might put them at higher SARS-Cov-2 infection risk. The effects of smoking on COVID-19 severity need clarification with larger studies. Additionally, the postulated protective effects of nicotine and nitric oxide, which may presumably reduce the risk of a "cytokine storm" in infected individuals, deserve assessment by controlled clinical trials.

Invited review: human air-liquid-interface organotypic airway tissue models derived from primary tracheobronchial epithelial cells-overview and perspectives

The lung is an organ that is directly exposed to the external environment. Given the large surface area and extensive ventilation of the lung, it is prone to exposure to airborne substances, such as pathogens, allergens, chemicals, and particulate matter. Highly elaborate and effective mechanisms have evolved to protect and maintain homeostasis in the lung. Despite these sophisticated defense mechanisms, the respiratory system remains highly susceptible to environmental challenges. Because of the impact of respiratory exposure on human health and disease, there has been considerable interest in developing reliable and predictive in vitro model systems for respiratory toxicology and basic research. Human air-liquid-interface (ALI) organotypic airway tissue models derived from primary tracheobronchial epithelial cells have in vivo–like structure and functions when they are fully differentiated. The presence of the air-facing surface allows conducting in vitro exposures that mimic human respiratory exposures. Exposures can be conducted using particulates, aerosols, gases, vapors generated from volatile and semi-volatile substances, and respiratory pathogens. Toxicity data have been generated using nanomaterials, cigarette smoke, e-cigarette vapors, environmental airborne chemicals, drugs given by inhalation, and respiratory viruses and bacteria. Although toxicity evaluations using human airway ALI models require further standardization and validation, this approach shows promise in supplementing or replacing in vivo animal models for conducting research on respiratory toxicants and pathogens.

Improving lung cancer risk stratification leveraging whole transcriptome RNA sequencing and machine learning across multiple cohorts

BACKGROUND

Bronchoscopy for suspected lung cancer has low diagnostic sensitivity, rendering many inconclusive results. The Bronchial Genomic Classifier (BGC) was developed to help with patient management by identifying those with low risk of lung cancer when bronchoscopy is inconclusive. The BGC was trained and validated on patients in the Airway Epithelial Gene Expression in the Diagnosis of Lung Cancer (AEGIS) trials. A modern patient cohort, the BGC Registry, showed differences in key clinical factors from the AEGIS cohorts, with less smoking history, smaller nodules and older age. Additionally, we discovered interfering factors (inhaled medication and sample collection timing) that impacted gene expressions and potentially disguised genomic cancer signals.

METHODS

In this study, we leveraged multiple cohorts and next generation sequencing technology to develop a robust Genomic Sequencing Classifier (GSC). To address demographic composition shift and interfering factors, we synergized three algorithmic strategies: 1) ensemble of clinical dominant and genomic dominant models; 2) development of hierarchical regression models where the main effects from clinical variables were regressed out prior to the genomic impact being fitted in the model; and 3) targeted placement of genomic and clinical interaction terms to stabilize the effect of interfering factors. The final GSC model uses 1232 genes and four clinical covariates - age, pack-years, inhaled medication use, and specimen collection timing.

RESULTS

In the validation set (N = 412), the GSC down-classified low and intermediate pre-test risk subjects to very low and low post-test risk with a specificity of 45% (95% CI 37-53%) and a sensitivity of 91% (95%CI 81-97%), resulting in a negative predictive value of 95% (95% CI 89-98%). Twelve percent of intermediate pre-test risk subjects were up-classified to high post-test risk with a positive predictive value of 65% (95%CI 44-82%), and 27% of high pre-test risk subjects were up-classified to very high post-test risk with a positive predictive value of 91% (95% CI 78-97%).

CONCLUSIONS

The GSC overcame the impact of interfering factors and achieved consistent performance across multiple cohorts. It demonstrated diagnostic accuracy in both down- and up-classification of cancer risk, providing physicians actionable information for many patients with inconclusive bronchoscopy.

Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and non-smokers, we generate a comprehensive atlas of epithelial cell types and states, connect these into lineages, and define cell-specific responses to smoking. Our analysis infers multi-state lineages that develop into surface mucus secretory and ciliated cells and then contrasts these to the unique specification of submucosal gland (SMG) cells. Accompanying knockout studies reveal that tuft-like cells are the likely progenitor of both pulmonary neuroendocrine cells and CFTR-rich ionocytes. Our smoking analysis finds that all cell types, including protected stem and SMG populations, are affected by smoking through both pan-epithelial smoking response networks and hundreds of cell-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium.

Transcriptomic changes in the nasal epithelium associated with diesel engine exhaust exposure

BACKGROUND

Diesel engine exhaust (DEE) exposure causes lung cancer, but the molecular mechanisms by which this occurs are not well understood.

OBJECTIVES

To assess transcriptomic alterations in nasal epithelium of DEE-exposed factory workers to better understand the cellular and molecular effects of DEE.

METHODS

Nasal epithelial brushings were obtained from 41 diesel engine factory workers exposed to relatively high levels of DEE (17.2-105.4 μg/m3), and 38 unexposed workers from factories without DEE exposure. mRNA was profiled for gene expression using Affymetrix microarrays. Linear modeling was used to identify differentially expressed genes associated with DEE exposure and interaction effects with current smoking status. Pathway enrichment among differentially expressed genes was assessed using EnrichR. Gene Set Enrichment Analysis (GSEA) was used to compare gene expression patterns between datasets.

RESULTS

225 genes had expression associated with DEE exposure after adjusting for smoking status (FDR q < 0.25) and were enriched for genes in pathways related to oxidative stress response, cell cycle pathways such as MAPK/ERK, protein modification, and transmembrane transport. Genes up-regulated in DEE-exposed individuals were enriched among the genes most up-regulated by cigarette smoking in a previously reported bronchial airway smoking dataset. We also found that the DEE signature was enriched among the genes most altered in two previous studies of the effects of acute DEE on PBMC gene expression. An exposure-response relationship was demonstrated between air levels of elemental carbon and the first principal component of the DEE signature.

CONCLUSIONS

A gene expression signature was identified for workers occupationally exposed to DEE that was altered in an exposure-dependent manner and had some overlap with the effects of smoking and the effects of acute DEE exposure. This is the first study of gene expression in nasal epithelial cells of workers heavily exposed to DEE and provides new insights into the molecular alterations that occur with DEE exposure.

Ellagic acid ameliorates lung damage in rats via modulating antioxidant activities, inhibitory effects on inflammatory mediators and apoptosis-inducing activities

Phytochemicals is considered one of the most effective and safe alternative therapy against oxidative linked lung diseases. Ellagic acid (EA), an important component of fruits, nuts, and vegetables, are partly responsible for their beneficial health effects against oxidation-related diseases. In the present study, we investigated the ameliorative effect of EA on lung damage induced by carbon tetrachloride (CCl₄) in Wistar male albino rats. Thirty-six male rats (n = 36, 8-week old) were divided into 4 groups, each with 9 rats. The groups were: Control group: received standard diet; EA group: administered with EA (10 mg/kg body weight, intraperitoneal); CCl₄ group: administered with CCl₄ (1.5 mg/kg body weight, intraperitoneal); EA+CCl₄ group: administered with EA and CCl₄. . The rats were decapitated at the end of experimental period of 8 weeks and the lung tissues were examined. CCl₄-induced rats showed elevation in the expressions of inflammatory proteins, nuclear factor kappa b (NF-κB), cyclooxygenase-2 (COX-2), and pro-inflammatory cytokine, tumor necrosis factor alpha (TNF-α); and the indicator of lipid peroxidation, malondialdehyde (MDA). Intraperitoneal administration of EA significantly reduced the levels of these markers. EA administration increased the protein expression levels of nuclear factor erythroid 2-related factor 2 (Nrf-2) and enhanced the activity of glutathione (GSH) and catalase enzyme (CAT). In addition, EA administration increased the expression levels of the executioner protein of apoptosis, caspase-3, and decreasing pro-survival protein, B cell lymphoma-2 (Bcl-2). In conclusion, these results establishes the protective role of EA in the treatment of lung damage and that in the future, this may have the potential to be used as a medication for the prevention or attenuation of lung diseases. Graphical abstract.

Biomarkers of Potential Harm: Summary of an FDA-Sponsored Public Workshop

Introduction

Since 2009, the United States (US) Food and Drug Administration (FDA) Center for Tobacco Products (CTP) has had the authority to regulate the manufacture, distribution, and marketing of tobacco products in order to reduce the death and disease caused by tobacco use. Biomarkers could play an important role across a number of FDA regulatory activities, including assessing new and modified risk tobacco products and identifying and evaluating potential product standards.

Methods

On April 4-5, 2016, FDA/CTP hosted a public workshop focused on biomarkers of potential harm (BOPH) with participants from government, industry, academia, and other organizations. The workshop was divided into five sessions focused on: (1) overview of BOPH; (2) cardiovascular disease (CVD); (3) chronic obstructive pulmonary disease (COPD); (4) cancer; and (5) new areas of research.

Results and Conclusions

The deliberations from the workshop noted some promising BOPH but also highlighted the lack of systematic effort to identify BOPH that would have utility and validity for evaluating tobacco products. Research areas that could further strengthen the applicability of BOPH to tobacco regulatory science include the exploration of composite biomarkers as predictors of disease risk, "omics" biomarkers, and examining biomarkers using existing cohorts, surveys, and experimental studies.

Implications

This paper synthesizes the main findings from the 2016 FDA-sponsored workshop focused on BOPH and highlights research areas that could further strengthen the science around BOPH and their applicability to tobacco regulatory science.

Nickel-induced transcriptional changes persist post exposure through epigenetic reprogramming

Background

Nickel is an occupational and environmental toxicant associated with a number of diseases in humans including pulmonary fibrosis, bronchitis and lung and nasal cancers. Our earlier studies showed that the nickel-exposure-induced genome-wide transcriptional changes, which persist even after the termination of exposure may underlie nickel pathogenesis. However, the mechanisms that drive nickel-induced persistent changes to the transcriptome remain elusive.

Results

To elucidate the mechanisms that underlie nickel-induced long-term transcriptional changes, in this study, we examined the transcriptome and the epigenome of human lung epithelial cells during nickel exposure and after the termination of exposure. We identified two categories of persistently differentially expressed genes: (i) the genes that were differentially expressed during nickel exposure; and (ii) the genes that were differentially expressed only after the termination of exposure. Interestingly, > 85% of the nickel-induced gene expression changes occurred only after the termination of exposure. We also found extensive genome-wide alterations to the activating histone modification, H3K4me3, after the termination of nickel exposure, which coincided with the post-exposure gene expression changes. In addition, we found significant post-exposure alterations to the repressive histone modification, H3K27me3.

Conclusion

Our results suggest that while modest first wave of transcriptional changes occurred during nickel exposure, extensive transcriptional changes occurred during a second wave of transcription for which removal of nickel ions was essential. By uncovering a new category of transcriptional and epigenetic changes, which occur only after the termination of exposure, this study provides a novel understanding of the long-term deleterious consequences of nickel exposure on human health.

Characterizing smoking-induced transcriptional heterogeneity in the human bronchial epithelium at single-cell resolution

The human bronchial epithelium is composed of multiple distinct cell types that cooperate to defend against environmental insults. While studies have shown that smoking alters bronchial epithelial function and morphology, its precise effects on specific cell types and overall tissue composition are unclear. We used single-cell RNA sequencing to profile bronchial epithelial cells from six never and six current smokers. Unsupervised analyses led to the characterization of a set of toxin metabolism genes that localized to smoker ciliated cells, tissue remodeling associated with a loss of club cells and extensive goblet cell hyperplasia, and a previously unidentified peri-goblet epithelial subpopulation in smokers who expressed a marker of bronchial premalignant lesions. Our data demonstrate that smoke exposure drives a complex landscape of cellular alterations that may prime the human bronchial epithelium for disease.

Effect of Intermittent Versus Continuous Low-Dose Aspirin on Nasal Epithelium Gene Expression in Current Smokers: A Randomized, Double-Blinded Trial

A chemopreventive effect of aspirin (ASA) on lung cancer risk is supported by epidemiologic and preclinical studies. We conducted a randomized, double-blinded study in current heavy smokers to compare modulating effects of intermittent versus continuous low-dose ASA on nasal epithelium gene expression and arachidonic acid (ARA) metabolism. Fifty-four participants were randomized to intermittent (ASA 81 mg daily for one week/placebo for one week) or continuous (ASA 81 mg daily) for 12 weeks. Low-dose ASA suppressed urinary prostaglandin E2 metabolite (PGEM; change of -4.55 ± 11.52 from baseline 15.44 ± 13.79 ng/mg creatinine for arms combined, P = 0.02), a surrogate of COX-mediated ARA metabolism, but had minimal effects on nasal gene expression of nasal or bronchial gene-expression signatures associated with smoking, lung cancer, and chronic obstructive pulmonary disease. Suppression of urinary PGEM correlated with favorable changes in a smoking-associated gene signature (P < 0.01). Gene set enrichment analysis (GSEA) showed that ASA intervention led to 1,079 enriched gene sets from the Canonical Pathways within the Molecular Signatures Database. In conclusion, low-dose ASA had minimal effects on known carcinogenesis gene signatures in nasal epithelium of current smokers but results in wide-ranging genomic changes in the nasal epithelium, demonstrating utility of nasal brushings as a surrogate to measure gene-expression responses to chemoprevention. PGEM may serve as a marker for smoking-associated gene-expression changes and systemic inflammation. Future studies should focus on NSAIDs or agent combinations with broader inhibition of pro-inflammatory ARA metabolism to shift gene signatures in an anti-carcinogenic direction.

Gene network approach reveals co-expression patterns in nasal and bronchial epithelium

Nasal gene expression profiling is a new approach to investigate the airway epithelium as a biomarker to study the activity and treatment responses of obstructive pulmonary diseases. We investigated to what extent gene expression profiling of nasal brushings is similar to that of bronchial brushings. We performed genome wide gene expression profiling on matched nasal and bronchial epithelial brushes from 77 respiratory healthy individuals. To investigate differences and similarities among regulatory modules, network analysis was performed on correlated, differentially expressed and smoking-related genes using Gaussian Graphical Models. Between nasal and bronchial brushes, 619 genes were correlated and 1692 genes were differentially expressed (false discovery rate <0.05, |Fold-change|>2). Network analysis of correlated genes showed pro-inflammatory pathways to be similar between the two locations. Focusing on smoking-related genes, cytochrome-P450 pathway related genes were found to be similar, supporting the concept of a detoxifying response to tobacco exposure throughout the airways. In contrast, cilia-related pathways were decreased in nasal compared to bronchial brushes when focusing on differentially expressed genes. Collectively, while there are substantial differences in gene expression between nasal and bronchial brushes, we also found similarities, especially in the response to the external factors such as smoking.

Lifetime Smoking History and Risk of Lung Cancer: Results From the Framingham Heart Study

Background

The relative risk of lung cancer decreases with years since quitting (YSQ) smoking, but risk beyond 25 YSQ remains unclear. Current lung cancer screening guidelines, which exclude smokers with more than 15 YSQ, may not detect lung cancers in this population.

Methods

We analyzed data from Framingham Heart Study Original (n = 3905) and Offspring cohort (n = 5002) participants for lifetime smoking and lung cancer incidence from 1954 to 1958 (Exam 4) and 1971 to 1975 (Exam 1), respectively, through 2013. We used multivariable-adjusted Cox proportional hazards regression models to compare current, former, and never smokers and lung cancer risk. Smoking status and covariates were time-updated every two years (Original) or four years (Offspring). Primary analyses were restricted to heavy ever smokers with more than 21.3 pack-years; additional analyses included all ever smokers.

Results

On follow-up (median = 28.7 years), 284 lung cancers were detected: incidence rates/1000 person-years in current, former, and never smokers were 1.97 (95% confidence interval [CI] = 1.66 to 2.33), 1.61 (95% CI = 1.34 to 1.93), and 0.26 (95% CI = 0.17 to 0.39), respectively. Heavy former (vs never) smokers had elevated lung cancer risk at all YSQ (<5: hazard ratio [HR] = 12.12, 95% CI = 6.94 to 21.17; 5-9: HR = 11.77, 95% CI = 6.78 to 20.45; 10-14: HR = 7.81, 95% CI = 3.98 to 15.33; 15-24: HR = 5.88, 95% CI = 3.19-10.83; ≥25: HR = 3.85, 95% CI = 1.80 to 8.26). Heavy former (vs current) smokers had 39.1% lower lung cancer risk within five YSQ. Among all former smokers, 40.8% of lung cancers occurred after more than 15 YSQ.

Conclusions

Among heavy former smokers, lung cancer risk drops within five YSQ relative to continuing smokers, yet it remains more than threefold higher than never smokers after 25 YSQ. Four of ten lung cancers occurred in former smokers with more 15 YSQ, beyond the screening window of the current guideline.

Gene Expression Alterations in the Bronchial Epithelium of e-Cigarette Users

BACKGROUND

Although e-cigarette (ECIG) use has increased in the United States, their potential health effects remain uncertain. Understanding the effects of tobacco cigarette (TCIG) smoke on bronchial airway epithelial gene expression have previously provided insights into tobacco-related disease pathogenesis. Identifying the impact of ECIGs on airway gene expression could provide insights into their potential long-term health effects. We sought to compare the bronchial airway gene-expression profiles of former TCIG smokers now using ECIGs with the profiles of former and current TCIG smokers.

METHODS

We performed gene-expression profiling of bronchial epithelial cells collected from current TCIG smokers (n = 9), current ECIG users who are former TCIG smokers (n = 15), and former TCIG smokers (n = 21). We then compared our findings with previous studies of the effects of TCIG use on bronchial epithelium, as well an in vitro model of ECIG exposure.

RESULTS

Among 3,165 genes whose expression varied between the three study groups (q < 0.05), we identified 468 genes altered in ECIG users relative to former smokers (P < .05). Seventy-nine of these genes were up- or down-regulated concordantly among ECIG and TCIG users. We did not detect ECIG-associated gene-expression changes in known pathways associated with TCIG usage. Genes downregulated in ECIG users are enriched among the genes most downregulated by exposure of airway epithelium to ECIG vapor in vitro.

CONCLUSIONS

ECIGs induce both distinct and shared patterns of gene expression relative to TCIGs in the bronchial airway epithelium. The concordance of the genes altered in ECIG users and in the in vitro study suggests that genes altered in ECIG users are likely to be changed as the direct effect of ECIG exposure.

Tobacco-Related Alterations in Airway Gene Expression are Rapidly Reversed Within Weeks Following Smoking-Cessation

The physiologic response to tobacco smoke can be measured by gene-expression profiling of the airway epithelium. Temporal resolution of kinetics of gene-expression alterations upon smoking-cessation might delineate distinct biological processes that are activated during recovery from tobacco smoke exposure. Using whole genome gene-expression profiling of individuals initiating a smoking-cessation attempt, we sought to characterize the kinetics of gene-expression alterations in response to short-term smoking-cessation in the nasal epithelium. RNA was extracted from the nasal epithelial of active smokers at baseline and at 4, 8, 16, and 24-weeks after smoking-cessation and put onto Gene ST arrays. Gene-expression levels of 119 genes were associated with smoking-cessation (FDR < 0.05, FC ≥1.7) with a majority of the changes occurring by 8-weeks and a subset changing by 4-weeks. Genes down-regulated by 4- and 8-weeks post-smoking-cessation were involved in xenobiotic metabolism and anti-apoptotic functions respectively. These genes were enriched among genes previously found to be induced in smokers and following short-term in vitro exposure of airway epithelial cells to cigarette smoke (FDR < 0.05). Our findings suggest that the nasal epithelium can serve as a minimally-invasive tool to measure the reversible impact of smoking and broadly, may serve to assess the physiological impact of changes in smoking behavior.

Molecular subtyping reveals immune alterations associated with progression of bronchial premalignant lesions

Bronchial premalignant lesions (PMLs) are precursors of lung squamous cell carcinoma, but have variable outcome, and we lack tools to identify and treat PMLs at risk for progression to cancer. Here we report the identification of four molecular subtypes of PMLs with distinct differences in epithelial and immune processes based on RNA-Seq profiling of endobronchial biopsies from high-risk smokers. The Proliferative subtype is enriched with bronchial dysplasia and exhibits up-regulation of metabolic and cell cycle pathways. A Proliferative subtype-associated gene signature identifies subjects with Proliferative PMLs from normal-appearing uninvolved large airway brushings with high specificity. In progressive/persistent Proliferative lesions expression of interferon signaling and antigen processing/presentation pathways decrease and immunofluorescence indicates a depletion of innate and adaptive immune cells compared with regressive lesions. Molecular biomarkers measured in PMLs or the uninvolved airway can enhance histopathological grading and suggest immunoprevention strategies for intercepting the progression of PMLs to lung cancer.

Multiomic analysis of oral keratinocytes chronically exposed to shisha

BACKGROUND

Tobacco is smoked in different form including cigarettes and water pipes. One popular form of water pipe smoking especially in Middle Eastern countries is shisha smoking. Shisha has been associated with various diseases including oral cancer. However, genomic alterations and gene expression changes associated with chronic shisha exposure have not been previously investigated.

OBJECTIVES

Whole-exome sequencing and gene expression profiling of immortalized human oral keratinocytes (OKF6/TERT1) cells chronically treated with 0.5% shisha extract for a period of 8 months was undertaken to characterize molecular alterations associated with shisha exposure.

METHODS

Genomic DNA and RNA were extracted and preprocessed as per manufacturer's instruction and subjected to whole-exome and transcriptome sequencing using Illumina HiSeq2500 platform. Exome was analyzed using GATK pipeline whereas RNA-Seq data was analyzed using HiSat2 and HTSeq along with DESeq to elucidate differentially expressed genes.

RESULTS

Whole-exome sequence analysis led to identification of 521 somatic missense variants corresponding to 389 genes RNA-Seq data revealed 247 differentially expressed genes (≥2-fold, P-value<0.01) in shisha treated cells compared to parental cells. Pathway analysis of differentially expressed genes revealed that interferon-signaling pathway was significantly affected. We predict activation of MAPK1 pathway which is known to play a key role in oral cancer. We also observed allele specific expression of mutant LIMA1 based on RNA-Seq dataset.

CONCLUSION

Our findings provide insights into genomic alterations and gene expression pattern associated with oral keratinocytes chronically exposed to shisha.

Assessment of a Highly Multiplexed RNA Sequencing Platform and Comparison to Existing High-Throughput Gene Expression Profiling Techniques

The need to reduce per sample cost of RNA-seq profiling for scalable data generation has led to the emergence of highly multiplexed RNA-seq. These technologies utilize barcoding of cDNA sequences in order to combine multiple samples into a single sequencing lane to be separated during data processing. In this study, we report the performance of one such technique denoted as sparse full length sequencing (SFL), a ribosomal RNA depletion-based RNA sequencing approach that allows for the simultaneous sequencing of 96 samples and higher. We offer comparisons to well established single-sample techniques, including: full coverage Poly-A capture RNA-seq, microarrays, as well as another low-cost highly multiplexed technique known as 3' digital gene expression (3'DGE). Data was generated for a set of exposure experiments on immortalized human lung epithelial (AALE) cells in a two-by-two study design, in which samples received both genetic and chemical perturbations of known oncogenes/tumor suppressors and lung carcinogens. SFL demonstrated improved performance over 3'DGE in terms of coverage, power to detect differential gene expression, and biological recapitulation of patterns of differential gene expression from in vivo lung cancer mutation signatures.

Longitudinal analysis of epigenome-wide DNA methylation reveals novel smoking-related loci in African Americans

Changes in DNA methylation may be a potential mechanism that mediates the effects of smoking on physiological function and subsequent disease risk. Given the dynamic nature of the epigenome, longitudinal studies are indispensable for investigating smoking-induced methylation changes over time. Using blood samples collected approximately five years apart in 380 African Americans (mean age 60.7 years) from the Genetic Epidemiology Network of Arteriopathy (GENOA) study, we measured DNA methylation levels using Illumina HumanMethylation BeadChips. We evaluated the association between Phase 1 smoking status and rate of methylation change, using generalized estimating equation models. Among the 6958 CpG sites examined, smoking status was associated with methylation change for 22 CpGs (false discovery rate q < 0.1), with the majority (91%) becoming less methylated over time. Methylation change was greater in ever smokers than never smokers, and the absolute differences in rates of change ranged from 0.18 to 0.77 per decade in M value, equivalent to a β value change of 0.013 to 0.047 per decade. Significant enrichment was observed for CpG islands, enhancers, and DNAse hypersensitivity sites (p < 0.05). Although biological pathway analyses were not significant, most of the 22 CpGs were within genes known to be associated with cardiovascular disease, cancers, and aging. In conclusion, we identified epigenetic signatures for cigarette smoking that may have been missed in cross-sectional analyses, providing insight into the epigenetic effect of smoking and highlighting the importance of longitudinal analysis in understanding the dynamic human epigenome.

The roles of microRNAs in the pathogenesis of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is characterized by a progressive and irreversible airflow obstruction, with an abnormal lung function. The etiology of COPD correlates with complex interactions between environmental and genetic determinants. However, the exact pathogenesis of COPD is obscure although it involves multiple aspects including oxidative stress, imbalance between proteolytic and anti-proteolytic activity, immunity and inflammation, apoptosis, and repair and destruction in both airways and lungs. Many genes have been demonstrated to be involved in those pathogenic processes of this disease in patients exposed to harmful environmental factors. Previous reports have investigated promising microRNAs (miRNAs) to disclose the molecular mechanisms for COPD development induced by different environmental exposure and genetic predisposition encounter, and find some potential miRNA biomarkers for early diagnosis and treatment targets of COPD. In this review, we summarized the expression profiles of the reported miRNAs from studies of COPD associated with environmental risk factors including cigarette smoking and air pollution exposures, and provided an overview of roles of those miRNAs in the pathogenesis of the disease. We also highlighted the potential utility and limitations of miRNAs serving as diagnostic biomarkers and therapeutic targets for COPD.