Sui Qihai, Liang Jiaqi, Hu Zhengyang, Chen Zhencong, Bi Guoshu, Huang Yiwei, Li Ming, Zhan Cheng, Lin Zongwu, Wang Qun
Department of Thoracic Surgery, Zhongshan Hospital, Fudan University, Shanghai, China.
Eight-Year Program Clinical Medicine, Grade of 2016, Shanghai Medical College, Fudan University, Shanghai, China.
Transl Lung Cancer Res. 2020 Aug;9(4):1407-1421. doi: 10.21037/tlcr-20-276.
Non-smoking-related lung adenocarcinoma (LUAD) has its own characteristics. Genetic and microenvironmental differences in smoking and non-smoking LUAD patients were analyzed to elucidate the oncogenesis of non-smoking-related LUAD, which will improve our understanding of the underlying molecular mechanism and be of clinical use in the future.
The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) databases were used for clinical and genomic information. Various bioinformatics tools were used to analyze differences in somatic mutations, RNA and microRNA (miRNA) expression, immune infiltration, and stemness indices. GO, KEGG, and GSVA analyses were performed with R. A merged protein-protein interaction (PPI) network was constructed and analyzed. A miRNA-differentially expressed gene network was constructed with miRNet. qRT-PCR was used for validation of 4 most significantly differently expressed genes and 2 miRNAs in tumor samples obtained from 20 pairs of non-smoking and smoking patients.
Five hundred and one patients with LUAD were obtained, including 210 in the non-smoking group and 292 in the smoking group. A total of 174 significantly altered somatic mutations were detected, including mutations in tumor protein p53 and epidermal growth factor receptor, which were downregulated in non-smoking-related LUAD. At the RNA level, 231 significantly differentially expressed genes were obtained; 124 were upregulated and 107 downregulated in the non-smoking group. GSVA analysis revealed 42 significant pathways. Other functional and enrichment analyses of somatic mutations and RNA expression levels revealed that these genes were significantly enriched in receptor activity regulation and receptor binding. Differences in microenvironments including immune infiltration (e.g., CD8 T cells and resting mast cells) and stemness indices were also found between groups. A 79-pair interaction was found between differentially expressed genes and miRNAs, of which miR-335-5p and miR-34a-5p were located in the center. Twenty-one genes, including vitronectin, neurotensin, and neuronatin, were differentially expressed in both non-smoking LUAD patients and DMSO-treated A549 cells. And the different expression of neurotensin, neuronatin, trefoil factor family2, regenerating family member 4, miR-377-5p, miR-34a were verified with the same tendency in our own samples.
Non-smoking LUAD patients, compared to smokers, have different characteristics in terms of somatic mutation, gene, and miRNA expression and the microenvironment, indicating a diverse mechanism of oncogenesis.
非吸烟相关肺腺癌(LUAD)有其自身特点。分析吸烟和非吸烟LUAD患者的基因及微环境差异,以阐明非吸烟相关LUAD的肿瘤发生机制,这将增进我们对潜在分子机制的理解,并在未来具有临床应用价值。
利用癌症基因组图谱(TCGA)、基因表达综合数据库(GEO)获取临床和基因组信息。使用多种生物信息学工具分析体细胞突变、RNA和微小RNA(miRNA)表达、免疫浸润及干性指数的差异。用R软件进行基因本体(GO)、京都基因与基因组百科全书(KEGG)及基因集变异分析(GSVA)。构建并分析合并的蛋白质-蛋白质相互作用(PPI)网络。用miRNet构建miRNA-差异表达基因网络。采用qRT-PCR对20对非吸烟和吸烟患者的肿瘤样本中4个差异最显著的表达基因和2个miRNA进行验证。
共纳入501例LUAD患者,其中非吸烟组210例,吸烟组292例。共检测到174个显著改变的体细胞突变,包括肿瘤蛋白p53和表皮生长因子受体的突变,这些突变在非吸烟相关LUAD中下调。在RNA水平,获得231个显著差异表达基因;非吸烟组中124个上调,107个下调。GSVA分析揭示了42条显著通路。对体细胞突变和RNA表达水平的其他功能及富集分析表明,这些基因在受体活性调节和受体结合方面显著富集。两组之间在包括免疫浸润(如CD8 T细胞和静息肥大细胞)和干性指数在内的微环境方面也存在差异。在差异表达基因和miRNA之间发现了79对相互作用,其中miR-335-5p和miR-34a-5p位于中心。包括玻连蛋白、神经降压素和神经调素在内的21个基因在非吸烟LUAD患者和二甲基亚砜处理的A549细胞中均有差异表达。并且在我们自己的样本中验证了神经降压素、神经调素、三叶因子家族2、再生家族成员4、miR-377-5p、miR-34a的差异表达具有相同趋势。
与吸烟者相比,非吸烟LUAD患者在体细胞突变、基因和miRNA表达以及微环境方面具有不同特征,表明肿瘤发生机制多样。