Cigarette smoke alters the transcriptome of non-involved lung tissue in lung adenocarcinoma patients

Characterization of lung adenocarcinoma based on immunophenotyping and constructing an immune scoring model to predict prognosis

Background: Lung cancer poses great threat to human health, and lung adenocarcinoma (LUAD) is the main subtype. Immunotherapy has become first line therapy for LUAD. However, the pathogenic mechanism of LUAD is still unclear. Methods: We scored immune-related pathways in LUAD patients using single sample gene set enrichment analysis (ssGSEA) algorithm, and further identified distinct immune-related subtypes through consistent clustering analysis. Next, immune signatures, Kaplan-Meier survival analysis, copy number variation (CNV) analysis, gene methylation analysis, mutational analysis were used to reveal differences between subtypes. pRRophetic method was used to predict the response to chemotherapeutic drugs (half maximal inhibitory concentration). Then, weighted gene co-expression network analysis (WGCNA) was performed to screen hub genes. Significantly, we built an immune score (IMscore) model to predict prognosis of LUAD. Results: Consensus clustering analysis identified three LUAD subtypes, namely immune-Enrich subtype (Immune-E), stromal-Enrich subtype (Stromal-E) and immune-Deprived subtype (Immune-D). Stromal-E subtype had a better prognosis, as shown by Kaplan-Meier survival analysis. Higher tumor purity and lower immune cell scores were found in the Immune-D subtype. CNV analysis showed that homologous recombination deficiency was lower in Stromal-E and higher in Immune-D. Likewise, mutational analysis found that the Stromal-E subtype had a lower mutation frequency in TP53 mutations. Difference in gene methylation (ZEB2, TWIST1, CDH2, CDH1 and CLDN1) among three subtypes was also observed. Moreover, Immune-E was more sensitive to traditional chemotherapy drugs Cisplatin, Sunitinib, Crizotinib, Dasatinib, Bortezomib, and Midostaurin in both the TCGA and GSE cohorts. Furthermore, a 6-gene signature was constructed to predicting prognosis, which performed better than TIDE score. The performance of IMscore model was successfully validated in three independent datasets and pan-cancer.

Genetic predisposition to lung adenocarcinoma outcome is a feature already present in patients' noninvolved lung tissue

Emerging evidence suggests that the prognosis of patients with lung adenocarcinoma can be determined from germline variants and transcript levels in nontumoral lung tissue. Gene expression data from noninvolved lung tissue of 483 lung adenocarcinoma patients were tested for correlation with overall survival using multivariable Cox proportional hazard and multivariate machine learning models. For genes whose transcript levels are associated with survival, we used genotype data from 414 patients to identify germline variants acting as cis-expression quantitative trait loci (eQTLs). Associations of eQTL variant genotypes with gene expression and survival were tested. Levels of four transcripts were inversely associated with survival by Cox analysis (CLCF1, hazard ratio [HR] = 1.53; CNTNAP1, HR = 2.17; DUSP14, HR = 1.78; and MT1F: HR = 1.40). Machine learning analysis identified a signature of transcripts associated with lung adenocarcinoma outcome that was largely overlapping with the transcripts identified by Cox analysis, including the three most significant genes (CLCF1, CNTNAP1, and DUSP14). Pathway analysis indicated that the signature is enriched for ECM components. We identified 32 cis-eQTLs for CNTNAP1, including 6 with an inverse correlation and 26 with a direct correlation between the number of minor alleles and transcript levels. Of these, all but one were prognostic: the six with an inverse correlation were associated with better prognosis (HR < 1) while the others were associated with worse prognosis. Our findings provide supportive evidence that genetic predisposition to lung adenocarcinoma outcome is a feature already present in patients' noninvolved lung tissue.

Cigarette Smoke Regulates the Expression of EYA4 via Alternation of DNA Methylation Status

Cigarette SMOKE (CS) considerably contributes to causing some diseases such as cancer, and it has a role in the alternation of gene expression through several mechanisms including epigenetics modification, particularly DNA methylation. EYA4 is one of the genes, that whose expression has been dysregulated in lung, colon, bladder, and breast cancer, leading to tumor progression. The alternation of DNA methylation levels has been implicated in regulating the expression of the EYA4 gene. Thus, in this study, we have shown the effect of CS on the DNA methylation level of the EYA4 promoter region as well as the methylation level on EYA4 expression. To determine the level of DNA methylation on the promoter region of the EYA4 gene, we have employed the bisulfite conversion treatment followed by the Sanger Sequence for 100 DNA samples taken from Saudi people (50 smokers and 50 nonsmokers). We found that 26% of DNA extracted from smoker samples is methylated, while there was no methylation identified in nonsmoker samples. Also, using the demethylating agents such as AZA on LoVo and Caco-2 cancer cell lines causes induction of transcription level of EYA4, implying the possible mechanism of DNA methylation in the upregulation of EYA4. These findings suggest the possible mechanism of CS in controlling the expression of EYA4 via changing the status of DNA methylation.

SPINK5 is a Tumor-Suppressor Gene Involved in the Progression of Nonsmall Cell Lung Carcinoma through Negatively Regulating PSIP1

This study aimed to elucidate how SPINK5 affects the malignant phenotypes of NSCLC and the molecular mechanism. NSCLC and adjacent normal tissues were collected to detect the differential level of SPINK5. The influence of SPINK5 on pathological indicators of NSCLC was analyzed. Cellular functions of NSCLC cells overexpressing SPINK5 were assessed by CCK-8, EdU, and transwell assay. By confirming the downstream target of SPINK5, its molecular mechanism on regulating NSCLC was finally explored through rescue experiments. SPINK5 was lowly expressed in NSCLC tissues, and it predicted tumor staging and lymphatic metastasis. In vitro overexpression of SPINK5 declined proliferative and migratory rates in NSCLC cells. PSIP1 was verified as the target gene binding SPINK5, and they displayed a negative correlation in NSCLC tissues. Overexpression of PSIP1 was able to reverse the inhibited proliferative and migratory potentials in NSCLC cells overexpressing SPINK5. SPINK5 level has a close relation to tumor staging and lymphatic metastasis in NSCLC. It serves as a tumor-suppressor gene that inhibits proliferation and migration of NSCLC through negatively regulating PSIP1.

Epigenetic Regulation of F2RL3 Associates With Myocardial Infarction and Platelet Function

BACKGROUND

DNA hypomethylation at the F2RL3 (F2R like thrombin or trypsin receptor 3) locus has been associated with both smoking and atherosclerotic cardiovascular disease; whether these smoking-related associations form a pathway to disease is unknown. F2RL3 encodes protease-activated receptor 4, a potent thrombin receptor expressed on platelets. Given the role of thrombin in platelet activation and the role of thrombus formation in myocardial infarction, alterations to this biological pathway could be important for ischemic cardiovascular disease.

METHODS

We conducted multiple independent experiments to assess whether DNA hypomethylation at F2RL3 in response to smoking is associated with risk of myocardial infarction via changes to platelet reactivity. Using cohort data (N=3205), we explored the relationship between smoking, DNA hypomethylation at F2RL3, and myocardial infarction. We compared platelet reactivity in individuals with low versus high DNA methylation at F2RL3 (N=41). We used an in vitro model to explore the biological response of F2RL3 to cigarette smoke extract. Finally, a series of reporter constructs were used to investigate how differential methylation could impact F2RL3 gene expression.

RESULTS

Observationally, DNA methylation at F2RL3 mediated an estimated 34% of the smoking effect on increased risk of myocardial infarction. An association between methylation group (low/high) and platelet reactivity was observed in response to PAR4 (protease-activated receptor 4) stimulation. In cells, cigarette smoke extract exposure was associated with a 4.9% to 9.3% reduction in DNA methylation at F2RL3 and a corresponding 1.7-(95% CI, 1.2-2.4, P=0.04) fold increase in F2RL3 mRNA. Results from reporter assays suggest the exon 2 region of F2RL3 may help control gene expression.

CONCLUSIONS

Smoking-induced epigenetic DNA hypomethylation at F2RL3 appears to increase PAR4 expression with potential downstream consequences for platelet reactivity. Combined evidence here not only identifies F2RL3 DNA methylation as a possible contributory pathway from smoking to cardiovascular disease risk but from any feature potentially influencing F2RL3 regulation in a similar manner.

Pregnancy exposure to phthalates and DNA methylation in male placenta - An epigenome-wide association study

BACKGROUND

Exposure to phthalates during pregnancy may alter DNA methylation in the placenta, a crucial organ for the growth and development of the fetus.

OBJECTIVES

We studied associations between urinary concentrations of phthalate biomarkers during pregnancy and placental DNA methylation.

METHODS

We measured concentrations of 11 phthalate metabolites in maternal spot urine samples collected between 22 and 29 gestational weeks in 202 pregnant women. We analyzed DNA methylation levels in placental tissue (fetal side) collected at delivery. We first investigated changes in global DNA methylation of repetitive elements Alu and LINE-1. We then performed an adjusted epigenome-wide association study using IlluminaHM450 BeadChips and identified differentially methylated regions (DMRs) associated with phthalate exposure.

RESULTS

Monobenzyl phthalate concentration was inversely associated with placental methylation of Alu repeats. Moreover, all phthalate biomarkers except for monocarboxy-iso-octyl phthalate and mono(2-ethyl-5-hydroxyhexyl) phthalate were associated with at least one DMR. All but three DMRs showed increased DNA methylation with increased phthalate exposure. The largest identified DMR (22 CpGs) was positively associated with monocarboxy-iso-nonyl phthalate and encompassed heat shock proteins (HSPA1A, HSPA1L). The remaining DMRs encompassed transcription factors and nucleotide exchange factors, among other genes.

CONCLUSIONS

This is the first description of genome-wide modifications of placental DNA methylation in association with pregnancy exposure to phthalates. Our results suggest epigenetic mechanisms by which exposure to these compounds could affect fetal development. Of interest, four identified DMRs had been previously associated with maternal smoking, which may suggest particular sensitivity of these genomic regions to the effect of environmental contaminants.

Bioinformatics Identification of Key Genes for the Development and Prognosis of Lung Adenocarcinoma

Objective: Lung adenocarcinoma (LUAD) is a common malignant tumor with a poor prognosis. The present study aimed to screen the key genes involved in LUAD development and prognosis. Methods: The transcriptome data for 515 LUAD and 347 normal samples were downloaded from The Cancer Genome Atlas and Genotype Tissue Expression databases. The weighted gene co-expression network and differentially expressed genes were used to identify the central regulatory genes for the development of LUAD. Univariate Cox, LASSO, and multivariate Cox regression analyses were utilized to identify prognosis-related genes. Results: The top 10 central regulatory genes of LUAD included IL6, PECAM1, CDH5, VWF, THBS1, CAV1, TEK, HGF, SPP1, and ENG. Genes that have an impact on survival included PECAM1, HGF, SPP1, and ENG. The favorable prognosis genes included KDF1, ZNF691, DNASE2B, and ELAPOR1, while unfavorable prognosis genes included RPL22, ENO1, PCSK9, SNX7, and LCE5A. The areas under the receiver operating characteristic curves of the risk score model in the training and testing datasets were .78 and .758, respectively. Conclusion: Bioinformatics methods were used to identify genes involved in the development and prognosis of LUAD, which will provide a basis for further research on the treatment and prognosis of LUAD.

Chronic E-Cigarette Aerosol Inhalation Alters the Immune State of the Lungs and Increases ACE2 Expression, Raising Concern for Altered Response and Susceptibility to SARS-CoV-2

Conventional smoking is known to both increase susceptibility to infection and drive inflammation within the lungs. Recently, smokers have been found to be at higher risk of developing severe forms of coronavirus disease 2019 (COVID-19). E-cigarette aerosol inhalation (vaping) has been associated with several inflammatory lung disorders, including the recent e-cigarette or vaping product use-associated lung injury (EVALI) epidemic, and recent studies have suggested that vaping alters host susceptibility to pathogens such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To assess the impact of vaping on lung inflammatory pathways, including the angiotensin-converting enzyme 2 (ACE2) receptor known to be involved in SARS-CoV-2 infection, mice were exposed to e-cigarette aerosols for 60 min daily for 1-6 months and underwent gene expression analysis. Hierarchical clustering revealed extensive gene expression changes occurred in the lungs of both inbred C57BL/6 mice and outbred CD1 mice, with 2,933 gene expression changes in C57BL/6 mice, and 2,818 gene expression changes in CD1 mice (>abs 1.25-fold change). Particularly, large reductions in IgA and CD4 were identified, indicating impairment of host responses to pathogens via reductions in immunoglobulins and CD4 T cells. CD177, facmr, tlr9, fcgr1, and ccr2 were also reduced, consistent with diminished host defenses via decreased neutrophils and/or monocytes in the lungs. Gene set enrichment (GSE) plots demonstrated upregulation of gene expression related to cell activation specifically in neutrophils. As neutrophils are a potential driver of acute lung injury in COVID-19, increased neutrophil activation in the lungs suggests that vapers are at higher risk of developing more severe forms of COVID-19. The receptor through which SARS-CoV-2 infects host cells, ACE2, was found to have moderate upregulation in mice exposed to unflavored vape pens, and further upregulation (six-fold) with JUUL mint aerosol exposure. No changes were found in mice exposed to unflavored Mod device-generated aerosols. These findings suggest that specific vaping devices and components of e-liquids have an effect on ACE2 expression, thus potentially increasing susceptibility to SARS-CoV-2. In addition, exposure to e-cigarette aerosols both with and without nicotine led to alterations in eicosanoid lipid profiles within the BAL. These data demonstrate that chronic, daily inhalation of e-cigarette aerosols fundamentally alters the inflammatory and immune state of the lungs. Thus, e-cigarette vapers may be at higher risk of developing infections and inflammatory disorders of the lungs.

Read-through transcripts in lung: germline genetic regulation and correlation with the expression of other genes

Transcripts originating from the transcriptional read through of two adjacent, similarly oriented genes have been identified in normal and neoplastic tissues, but their functional role and the mechanisms that regulate their expression are mostly unknown. Here, we investigated whether the expression of read-through transcripts previously identified in the non-involved lung tissue of lung adenocarcinoma patients was genetically regulated. Data on genome-wide single nucleotide variant genotypes and expression levels of 10 read-through transcripts in 201 samples of lung tissue were combined to identify expression quantitative trait loci (eQTLs). Then, to identify genes whose expression levels correlated with the 10 read-through transcripts, we used whole transcriptome profiles available for 154 patients. For 8 read-though transcripts, we identified 60 eQTLs (false discovery rate <0.05), including 17 cis-eQTLs and 43 trans-eQTLs. These eQTLs did not maintain their behavior on the 'parental' genes involved in the read-through transcriptional event. The expression levels of 7 read-through transcripts were found to correlate with the expression of other genes: CHIA-PIFO and CTSC-RAB38 correlated with CHIA and RAB38, respectively, while 5 other read-through transcripts correlated with 43 unique non-parental transcripts; thus offering indications about the molecular processes in which these chimeric transcripts may be involved. We confirmed 9 eQTLs (for 4 transcripts) in the non-involved lung tissue from an independent series of 188 lung adenocarcinoma patients. Therefore, this study indicates that the expression of four read-through transcripts in normal lung tissue is under germline genetic regulation, and that this regulation is independent of that of the genes involved in the read-through event.

Association between tobacco substance usage and a missense mutation in the tumor suppressor gene P53 in the Saudi Arabian population

The tumor suppressor gene TP53 and its downstream genes P21 and MDM2 play crucial roles in combating DNA damage at the G1/S cell cycle checkpoint. Polymorphisms in these genes can lead to the development of various diseases. This study was conducted to examine a potential association between tobacco substance usage (TSU) and single-nucleotide polymorphism (SNP) at the exon regions of the P53, P21, and MDM2 genes by comparing populations of smokers and non-smokers from Saudi Arabia. P53 rs1042522 (C/G), P21 rs1801270 (A/C), and MDM2 rs769412 (A/G) were investigated by genotyping 568 blood specimens: 283 from male/female smokers and 285 from male/female non-smokers. The results obtained from the smokers and their control non-smokers were compared according to age, sex, duration of smoking, and type of TSU. Heterozygous CG, homozygous GG, and CG+GG genotypes, as well as the G allele of rs1042522 were significantly associated with TSU in Saudi smokers compared with non-smokers. The C allele frequency of rs1801270 was also associated with TSU in smokers (OR = 1.33, p = 0.049) in comparison with non-smokers, in younger smokers (≤29 years) (OR = 1.556, p = 0.03280) in comparison with non-smokers of the same age, in smokers who had smoked cigarettes for seven years or less (OR = 1.596, p = 0.00882), and in smokers who had consumed shisha (OR = 1.608, p = 0.04104) in comparison with the controls. However, the genotypic and allelic frequencies for rs769412 did not show significant associations with TSU in Saudis. The selected SNP of P53 was strongly associated with TSU and may be linked to TSU-induced diseases in the Saudi Arabian population.

Prognostic Value of a Ten-Gene Signature in HNSCC Patients Based on Tumor-Associated Macrophages Expression Profiling

Tumor-associated macrophages (TAMs) are regarded as the most abundantly infiltrating immune cells around the tumor microenvironment (TME) in head and neck squamous cell carcinoma (HNSCC), which plays an essential role in immunosuppression and tumorigenesis. In the TCGA HNSCC cohort, 500 patients with clinical-pathological information and RNA sequence expression were randomly assigned to training for lasso regression and validation for verification, respectively. A TAM-based ten-gene signature (TBGs) was constructed, which divided the patients into high-risk and low-risk groups, could predict overall survival (OS) of HNSCC patients in the training dataset (p = 3.527e^-05) and validation dataset (p = 3.785e^-02). The result of Cox univariate and multivariate regression analyses showed that the risk score of TBGs could be an independent prognostic factor in HNSCC. ROC curve confirmed that the risk score of TBGs has good sensitivity and specificity for prognosis prediction (AUC = 0.659) and was also verified by the validation dataset (AUC = 0.621). We obtained key risk transcription factors (TFs)-EHF and SNAI2-by correlation analysis with TBGs. Moreover, we ran a gene set enrichment analysis (GSEA) to speculate that TBGs act on interstitial remodeling, tumor killing, metabolic reprogramming, and tumor immune-related pathways. Finally, we combined clinical-pathological features and risk score of TBGs to establish clinical nomograms, and calibration curves verified the accuracy of long-term clinical prognosis in the two datasets (C-index of 5-year OS = 0.721 and 0.716). In general, the TBGs we obtained may accurately predict the prognosis of HNSCC patients to provide personalized treatment.

Analysis of genomic and transcriptomic variations as prognostic signature for lung adenocarcinoma

Background

Lung cancer is the leading cause of the largest number of deaths worldwide and lung adenocarcinoma is the most common form of lung cancer. In order to understand the molecular basis of lung adenocarcinoma, integrative analysis have been performed by using genomics, transcriptomics, epigenomics, proteomics and clinical data. Besides, molecular prognostic signatures have been generated for lung adenocarcinoma by using gene expression levels in tumor samples. However, we need signatures including different types of molecular data, even cohort or patient-based biomarkers which are the candidates of molecular targeting.

Results

We built an R pipeline to carry out an integrated meta-analysis of the genomic alterations including single-nucleotide variations and the copy number variations, transcriptomics variations through RNA-seq and clinical data of patients with lung adenocarcinoma in The Cancer Genome Atlas project. We integrated significant genes including single-nucleotide variations or the copy number variations, differentially expressed genes and those in active subnetworks to construct a prognosis signature. Cox proportional hazards model with Lasso penalty and LOOCV was used to identify best gene signature among different gene categories.

We determined a 12-gene signature (BCHE, CCNA1, CYP24A1, DEPTOR, MASP2, MGLL, MYO1A, PODXL2, RAPGEF3, SGK2, TNNI2, ZBTB16) for prognostic risk prediction based on overall survival time of the patients with lung adenocarcinoma. The patients in both training and test data were clustered into high-risk and low-risk groups by using risk scores of the patients calculated based on selected gene signature. The overall survival probability of these risk groups was highly significantly different for both training and test datasets.

Conclusions

This 12-gene signature could predict the prognostic risk of the patients with lung adenocarcinoma in TCGA and they are potential predictors for the survival-based risk clustering of the patients with lung adenocarcinoma. These genes can be used to cluster patients based on molecular nature and the best candidates of drugs for the patient clusters can be proposed. These genes also have a high potential for targeted cancer therapy of patients with lung adenocarcinoma.

Serum levels of retinol-binding protein 4 and the risk of non-small cell lung cancer: A case-control study

Retinol binding protein 4 (RBP4), as an adipokine, has been identified to be associated with several types of cancer. However, no studies have assessed its effect on non-small cell lung cancer (NSCLC) risk. The objective of this study was to assess the association between serum RBP4 levels and the risk of NSCLC.A case-control study design was used to recruit 256 confirmed NSCLC cases and 256 age- and gender-matched healthy controls by frequency between August 2017 and January 2019. Serum RBP4 was measured using enzyme-linked immune absorbent assay before treatment. Unconditional logistic regression analysis was applied to estimate the odds ratio and 95% confidence interval (CI).Serum RBP4 level was significantly higher in NSCLC patients than those in the healthy control group (36.05 ± 8.28 vs 29.54 ± 7.71 μg/mL, P < .05). Higher serum RBP4 level was associated with increased risk of NSCLC (P trend = .001). Compare with those in the lowest tertile, the adjusted odds ratios were 1.85 (95% CIs 1.07-3.2) (P = .029) for the second tertile and 2.18 (95% CIs 1.37-3.45) (P = .001) for the highest tertile after adjusting for confounding variables. No interactions were observed after stratified analyses by body mass index and smoking status (P for interaction: .584 and .357).Our study indicated that serum RBP4 level was positively related to the risk of NSCLC. Additional studies with prospective design are required to confirm this finding.

Biomarker potential of repetitive-element transcriptome in lung cancer

Since repetitive elements (REs) account for nearly 53% of the human genome, profiling its transcription after an oncogenic change might help in the search for new biomarkers. Lung cancer was selected as target since it is the most frequent cause of cancer death. A bioinformatic workflow based on well-established bioinformatic tools (such as RepEnrich, RepBase, SAMTools, edgeR and DESeq2) has been developed to identify differentially expressed RNAs from REs. It was trained and tested with public RNA-seq data from matched sequencing of tumour and healthy lung tissues from the same patient to reveal differential expression within the RE transcriptome. Healthy lung tissues express a specific set of REs whose expression, after an oncogenic process, is strictly and specifically changed. Discrete sets of differentially expressed REs were found for lung adenocarcinoma, for small-cell lung cancer, and for both cancers. Differential expression affects more HERV-than LINE-derived REs and seems biased towards down-regulation in cancer cells. REs behaving consistently in all patients were tested in a different patient cohort to validate the proposed biomarkers. Down-regulation of AluYg6 and LTR18B was confirmed as potential lung cancer biomarkers, while up-regulation of HERVK11D-Int is specific for lung adenocarcinoma and up-regulation of UCON88 is specific for small cell lung cancer. Hence, the study of RE transcriptome might be considered another research target in cancer, making REs a promising source of lung cancer biomarkers.