基于高通量靶向测序技术的肿瘤原发部位差异研究.pdf

西北大学学报（自然科学版）2023年8 月，第53卷第4期，Aug.，2 0 2 3,V o l.53,No.4Journal of Northwest University(Natural Science Edition)JNWU肿瘤学基于高通量靶向测序技术的肿瘤原发部位差异研究罗祥1,鲁四海,戴鹏高,杨栓盈（1.西安交通大学第二附属医院呼吸科,陕西西安7 10 0 0 4;2.西北大学生命科学学院，陕西西安7 10 0 6 9）摘要探究肺内不同病灶部位非小细胞肺癌患者肿瘤突变特征及导致临床表现差异的潜在分子机制。对铜川市人民医院收治的6 2 例肺癌患者的新鲜肿瘤组织标本进行高通量测序，并收集TCGA数据库中50 2 份非小细胞肺癌测序数据，从不同肺叶、分布侧等角度分析基因的突变特征与功能学影响。结果显示,99.1%（559/56 4）的样本中均检出突变，各标本中双突变比例最高（6 3.0 1%8 2.93%），各组间均检出特有突变，但其与年龄、性别等无相关性。功能分析表明，叶分组中特有突变参与癌症中心碳代谢、药物抵抗等过程，而左、右侧分组中主要与PD-L1表达调控、细胞增殖、错配修复及凋亡等有关。该研究揭示了原发部位不同的肺癌患者组织突变特征及生物学功能差异，对非小细胞肺癌致病机制研究及临床辅助诊断原发灶具有重要意义。关键词原发部位；突变；基因功能；高通量测序；非小细胞肺癌中图分类号：Q786Study of differences in tumor primary sites based on high-throughputLUO Xiang2,LU Sihai?,DAI Penggao*,YANG Shuanying(1.Department of Respiratory Medicine,Second Hospital of Xian Jiaotong University,Xian 710004,China;2.The College of Life Sciences,Northwest University,Xi an 710069,China)Abstract To explore the tumor mutation characteristics of patients with non-small cell lung cancer at differentlesion locations in the lung and the underlying molecular mechanisms leading to differences in clinical manifes-tations.High-throughput sequencing was performed on fresh tumor tissue samples from 62 lung cancer patientsadmitted to our hospital,and the sequencing data of 502 non-small cell lung cancers in the TCGA databasewere collected,and the characteristics and functional characteristics of gene mutation were analyzed from theperspectives of dfferent lung lobes and distribution sides.Mutations were detected in 99.1%(559/564)ofthe samples.The proportion of samples with two mutations detected in different groups is the highest(63.01%82.93%).There were significant differences in the types and frequencies of mutated genes and loci amongthe groups,but there was no correlation between them and clinical characteristics(age,gender).The func-tional prediction results show that the non-shared genes between different lung lobe groups are mainly related tothe regulation of carbon metabolism,drug resistance in the cancer center,while the non-shared genes betweenthe left and right lung cancer groups are mainly related to the regulation of PD-L1 expression,cell prolifera-收稿日期：2 0 2 2-11-15基金项目：陕西省重点研发计划项目（2 0 2 2 ZDXM-SF-04）第一作者：罗祥，男，副主任医师，从事肺癌诊治相关研究，136 191957 18 16 。通信作者：杨栓盈，男，教授，从事肺癌发病机制、耐药和综合防治研究，。D0I:10.16152/ki.xdxbzr.2023-04-013targeted sequencing technology 620tion,mismatch repair and apoptosis signaling pathway,etc.There are significant differences in tissue mutationcharacteristics and biological function changes in lung cancer patients with different primary sites,which is ofgreat significance for the study of the pathogenesis of non-small cell lung cancer and the auxiliary clinical diag-nosis of the primary tumor.Keywords primary site;mutations;gene function;high-throughput sequencing;non-small cell lung cancer非小细胞肺癌是我国发病率、致死率较高的组）。收集2 0 18 年9 月至2 0 2 1年12 月在铜川市恶性肿瘤之一 ,其病因与致病机制尚不清楚。人民医院就诊、未经治疗、经病理学明确诊断且无研究显示,非小细胞肺癌患者的临床治疗响应及合并其他肿瘤的6 2 例原发性肺癌患者的新鲜肿预后等差异化临床表现可能与肺内原发灶分布位瘤组织样本。入组患者在被采集组织样本前均已置密切相关 2 4。目前,研究发现不同生理部位的被告知并签署知情同意书，同意收集其临床信息。结构与功能、局部血液供给等均被证实存在显著1.3DNA提取与纯化性差异 5-7 。同时，也探讨了不同病灶部位中分子采集符合纳入标准的肺癌新鲜组织样本,置遗传学水平如基因突变类型、变异频率及表达水于组织保存液中低温暂存。采集后2 小时内，依平的差异化特征，以及运用基因突变谱分析以解据PureLinkTM基因组DNA小量提取试剂盒进行决同步性非小细胞肺癌与肺内转移扩散辨别诊断基因组DNA提取与纯化,Nanodrop测定吸光度问题 8-9,但是尚无一致性结论，且国内缺乏多病值，A260/A280接近1.8 视为合格。灶部位、大样本量的临床研究。1.4靶基因高通量测序随着高通量测序技术的不断发展，靶向重测经检测质量合格的DNA样本送往陕西佰美序技术已经成为挖掘鉴定肿瘤驱动基因变异及药基因股份有限公司进行文库构建，并采用Ilumina物靶点基因的有效工具 10 。本研究通过肿瘤组测序平台针对52 个肺癌热点基因（AKT1、A LK、织靶向测序技术与肿瘤大数据挖掘相结合，旨在APC、A RA F、A T M、BC L2 L11、BRA F、BRC A 1、探究不同肺叶部位以及左侧肺、右侧肺等肿瘤起BRCA2、C D K N2 A、C H EK 2、C T NNB1、D D R2、EG-源差异的肺癌基因组突变谱特征及功能学影响，FR、EPCA M、ERBB2、ERCC1、FBXW 7、FG FR1、FG-进一步揭示导致原发部位不同的肺癌患者临床表FR2、FG FR3、I D H 1、I G F1R、JA K 2、K D R、K I T、现差异产生的潜在分子机制，以期对临床制定个KRAS、M A P2 K 1、M A P2 K 2、M ET、M LH 1、M SH 2、体化诊疗策略、改善患者临床结局提供更多的实MSH6、M T O R、NF1、NRA S、NT RK 1、PI K 3C A、验数据与理论支撑。PIK3R1、PM S2、PO LE、PT EN、RET、RO S1、SMARCA4、SM A RCB1、ST K 11、T ERT、T P53、T SC1、1材料与方法TSC2、T W I ST 1）的靶向深度测序，上述基因均为FDA批准的靶向药物或可用于评估术前术后化1.1TCGA数据获取与筛选疗药效及毒副作用的相关基因,各样本测得数据本研究中不同叶分组测序数据及临床信息来量约为 1 GB。源于 TCGA 数据库（the cancer genome atlas pro-1.5生物信息学分析gram,TCGA）,此数据库储存了大量的肿瘤突变数基于Linux操作系统对所收集肿瘤组织样本据集以及临床资料。下载获取TCGA数据库中50 2的测序数据进行了数据质量控制、突变信息挖掘例非小细胞肺腺癌患者的突变数据与临床信息数及功能注释和分析。Fastx toolkit(http：/h a n n o n-据,利用python软件对突变数据进行获取并将其lab.cshl.edu/fastx_toolkit/index.html）用于全外与临床信息数据进行匹配与筛选。依据临床诊断显子测序数据的质量验证和测序数据预处理。信息，将其分为上、中、下叶肺癌组（叶分组）。GATK（h t t p s：/g a t k.b r o a d in s t it u t e.o r g）软件用于1.2入组患者筛选各样本中突变检测及注释，突变报告标准为突变基于不同肺叶、不同分布侧原发性肺癌基因丰度0.1%。R脚本用于高频突变基因（突变频突变特征分析的研究思路，临床样本用以补充率排位前10）、位点以及非公有突变基因筛选。TCGA数据库缺失的左、右侧肺癌突变数据（侧分基因变异通常会导致蛋白表达或功能上的变化，西北大学学报（自然科学版）第53卷第4期从而对细胞内部的生物学过程产生影响,DAVID（h t t p s：/d a v i d.n c i f c r f.g o v/）数据库被用来进行突变特征基因的通路分析及基因功能注释。1.6数据统计分析各组样本中突变检出率与年龄关联性利用线性回归分析等方法实现。不同分组间样本量、突变频数分布、列联表等数据间差异显著性检验采用MannWhitney方法、Fisher精确检验或。单因素方差分析用于多组数据间差异显著性分析。各组间突变基因及位点差异聚类分析、韦恩图、信号通路等分析及绘图由ggplot2、h e a t ma p 2 和Venn等R（h t t p s：/w w w.r-p r o j e c t.o r g/）软件包完成。显著性判定水平为p65岁155female170male132Missense mutation588Nonsense mutation97Deletions69Insertions24注：单因素方差分析（*）或卡方检验（#)用于各组间差异的统计学检验,P0.05。2.2肿瘤突变基因及位点分布差异分析2.2.1叶分组叶分组中各组突变基因的检出种类与频率分布差异分析结果显示【见图1（a)】，上叶肺癌组突变基因种类覆盖率为10 0%（52/52），中叶、下叶肺癌组突变基因覆盖率分别为8 8.5%（46/52）和44.2%（2 3/52）。突变基因特征分析结果发现,相较于下叶肺癌组，中叶肺癌组未能检出AKT1、A RA F、BCL2 L11、CD K N2 A、CH EK 2、D D R2、EPCAM、ERBB2、ERC C 1、FG FR2、FG FR3、JA K 2、KIT、M A P2 K 1、M A P2 K 2、M LH 1、M SH 2、M SH 6、罗祥，等：基于高通量靶向测序技术的肿瘤原发部位差异研究Lower lobe2117900672149499612833227383324428 621部样本中均检出突变（50 2/50 2）。如表1所示，从各组收录样本来看，上叶肺癌组占比最高（6 0.1%），各组样本量均随年龄增高呈增多趋势,但各组样本在年龄分布上不存在显著性差异（p=0.2 8 48）,且在性别分类上各组间也不存在显著性差异(p=0.4541)。在突变导致氨基酸变化分析中发现，叶分组各组均以错义突变为主,但不同分组间未发现显著性差异（p=0.1629）。2.1.2侧分组左、右侧肺癌组共纳入样本6 2 例,9 1.9%(57/62)样本均可检出突变,其中5例无突变样本中2 例来自左侧肺癌组,其余3例为右侧肺癌组。收录样本量变化趋势与叶分组类似,但左、右侧肺癌组中样本的性别组成有明显差异（p=0.0183）。突变引起的氨基酸改变分析结果显示左、右侧组中也均以错义突变为主,且各分组间未发现显著差异（双侧组：p=0.1081）（见表1）。pvalueLeft side260120.284 8*14180.454 1*834250.162 9*0MTOR、NRA S、NT RK 1、PM S2、PT EN、RET、SMARCB1、T ERT、T SC1、T SC2、T W I ST 1等2 3个突变基因，约占总体的44.2%。相较于上叶肺癌组，下叶肺癌组未能于EPCAM、ERCC1、M A P2 K 2、NRAS、NT RK 1、T ERT 等6 个基因中检出突变。对突变位点水平的分析结果显示【见图1（c）】,共同检出的突变位点远低于非共有突变位点，其中上、中、下叶肺癌组分别检出7 2 7 种（7 2 7/8 0 0，90.9%）、32 种（32/46,6 9.6%）、2 6 9 种（2 6 9/341,78.9%）特有突变位点。在各样本可检出突变位点数量中,叶分组均表现为双突变类型检出Right side36012241422467155p value0.428 6#0.018 3#0.108 1#622频率最高,其次为单一突变类型或3种复合突变类型。此外,本研究分析了不同分组中各样本突变位点检出频率与年龄、性别等临床特征之间的关联性，仅发现在数据库来源的数据集中突变频率与年龄分布存在显著负相关(p=0.0196)。2.2.2侧分组侧分组中各组突变基因的检出种类与频率分布差异分析结果显示【见图1（b）】，侧分组中左、右侧肺癌组检出突变基因覆盖率仅为2 3.1%（12/52）和34.6%（18/52），非共有突变基因检出率分别为7.7%（BCL2L11、K D R、M SH 6、PI K 3R1,(a)Upper lobegroup610146西北大学学报（自然科学版）4/52）、19.2%（A LK、A PC、BRCA 2、CH EK 2、ERBB2、FG FR2、NF1、PM S2、PT EN、RET,10/52）,检出的突变基因种类明显低于叶分组来源的样本组。对突变位点水平的分析结果显示【见图1（c）,左、右侧肺癌组中分别检出34种（34/42，80.9%）、52 种（52/6 0,8 6.7%）。在各样本中可检出突变位点数量中，各组均表现为双突变类型检出频率最高,其次为单一突变类型或3种复合突变类型,基本与叶分组各组表现相近。但侧分组中未发现各组样本突变频率与年龄分布、性别存在关联性。Middle lobegroup(b)Lowerlobegroup027第53卷Middle lobegroup32Upper lobegroup07155012Lowerlobegroup4043226803Left sidegroupRight sidegroup(c)234568123456789101112145Middlelobegroup812-080O0(a)各组间突变基因一致性分析；(b)各组间突变位点一致性分析；（c)各组间单突变、双突变、多突变样本分布分析；（d)叶分组中各样本突变种类数与年龄的关联性分析图1不同原发部位肿瘤突变基因分布差异及突变数量与临床特征相关性分析Fig.1 Differences in the distribution of mutated genes in tumors at different primary sites and the correlation analysis be-tween mutations frequency and clinical features2.3高频突变基因及位点差异分析2.3.1叶分组为了进一步分析导致不同原发部位肺癌临床表现差异的潜在分子机制，本研究还分析了高频突变基因及位点分布差异。由图2(a)可知,叶分组中上、下叶肺癌组总体突变率最高的基因均为TP53(154/778,19.79%vs.66/358,18.44%)。但中叶肺癌组中以KRAS（10/44,2 2.7 3%）占比00335(d)Rightside groupTCGALLUAD870111213415Upper lobegroup2323121212_Left sidegroup26207TCGALUAD1545Lower lobe group45333667NGS_sample45Left side group45Right side group45915666最高，TP53（7/44，15.91%）次之。各组间高频突变基因突变率分布差异显著(p0.0001）。对高频突变位点分析结果显示【见图2（b）】，叶分组KRASp.G12C突变频率最高，其在上、中、下叶肺癌组中检出频率分别为8.9 4%（2 7/30 2）、28.57%(6/21)和11.7 3%(2 1/17 9)。2.3.2侧分组左、右侧肺癌组中总体突变率最高的基因分10P=0.0196500R2=0.01252040Age al diagnosis6080100第4期别为TP53（17/41，41.46%）和EGFR（2 4/7 3，32.88%）。各分组间高频突变基因突变率分布差异不显著（p=0.5357）,且各组间差异均与各组样本量无关（p=0.70）【见图2（a)】。两组位点水平分析结果显示【见图2（b)】,左、右侧肺癌组中以EGFR p.L858R和ECFR_p.E746_A750del突变频率居首，分别为30.7 7%（8/2 6）和3ERBB22BCL2L111CDKN2A0PTEN-1EGFR-2TP53-3KRASTCGALLUADUpperlobe罗祥，等：基于高通量靶向测序技术的肿瘤原发部位差异研究STK11KDRNF1SMARCA4BRCA2BRCA1ATMBRAFPOLEALKAPCPIK3CAIGF1RIDH1Middle lobeLowerlobe 623 13.89%（5/36）。侧分组中ECFR的突变率（左侧：2 9.2 7%，右侧：32.8 8%）显著低于叶分组（上叶：5.6 6%，中叶：2.2 7%，下叶：7.2 6%）（p=0.0015）。此外，侧分组与叶分组中突变位点检出频率差异明显，仅ECFRp.L858R和PIK3CAp.E545K等2 个基因同时分布于叶分组、侧分组高频突变基因中。4BRCA2_p.G2063RBRCA1_p.P1749SBRAF_p.G596SATM_p.1923*APC_p.G1312RALK_p.S439SALK_p.F607FALK_p.L80MTP53_p.R248WEGFR_CNVCDKN2A_p.D153NBRCA1_p.R841WCDKN2A_p.A118VNGS_sampleLeft sideKRAS_p.G12CKRAS_p.G12VEGFR_p.L858REGFR_p.E746_A750delEGFR_p.L747_P753delinsSPIK3CA_p.E542KTP53_p.Q192*EGFR_p.T79OMALK_EML4-ALKPIK3CA_p.E545KPIK3CA_CNVBCL2L11_BIM intron2deletionEGFR_CNVKRAS_p.G12DKRAS_p.G12ATP53_p.T125TBRAF_p.G466VERBB2_p.Y772_A775dupTP53_p.R273LBRAF_p.G469VBRAF_p.V600ECTNNB1_P.S37FEGFR_p.E709_T710delinsDALK_p.P666PRight sideTCGA_LUAD202-4UpperlobeMiddlelobeLowerlobeNGS_sampleLeft sideRight side(a)高频突变基因热图分析图2 各组中检出样本量高频突变基因及位点分布差异分析Fig.2 Analysis of the top 10 mutated genes and loci detected samples in tumors at different primary sites2.4非共有突变基因分子功能与通路富集分析2.4.1叶分组功能预测及通路富集结果显示【见图3（a）、图3（c），上、中、下叶肺癌组之间EPCAM、ER-CC1、M A P2 K 2、NRA S、NT RK 1、T ERT 等6 种非共有突变基因,可能与中心碳代谢（central carbonmetabolism in cancer）、EC FR激酶抑制剂抵抗(EGFR tyrosine kinase inhibitor resistance)、PI3K-Akt信号通路（PI3K-Akt signalingpathway）、MAPK信号通路（MAPKsignalingpathway）、内分泌抵抗（endocrine resistance）及药物抵抗等信号(b)高频突变位点热图分析通路调控相关,生物学功能富集的5个生物学过程均与小RNA表达调控相关。2.4.2侧分组左、右侧肺癌组中共检出 BCL2L11、K D R、MSH6、PI K 3R1、A LK、A PC、BRCA 2、CH EK 2、ERBB2、FG FR2、NF1、PM S2、PT EN、RET 等14个非共有突变基因,其通路分析结果【见图3（b)】与上、中、下叶肺癌组特有基因富集结果高度相似,均表明中心碳代谢（central carbon metabolismincancer）、EG FR激酶抑制剂抵抗（EGFRtyrosinekinase inhibitorresistance）等多种通路调控与之相 624关,但左、右侧肺癌发生过程还可能与PD-L1表达调控通路（PD-L1 expression and PD-1 check-point pathway in cancer）、钙离子信号通路（calci-um signaling pathway）、错配修复（mismatch re-pair)等通路调控有关。此外,其生物学功能富集(a)Central carbon metabolism in cancerEGFRtyrosinekinase inhibitorresistanceBreastcancerNon-small cell lung cancerEndometrialcancerProstatecancerMelanomaGliomaPancreatic cancerskemredPI3K-Akt signalingpathwayEndocrineresistanceHepatocellularcarcinomaBladdercancerMicroRNAs in cancerGastriccancerMAPK signaling pathwayPlatinum drugresistanceFoxO signaling pathwaySignaling pathways regulating pluripotency of stem cellsColorectalcancer()regulationof production of smallpositive regulation of stem cellagingdsR processingproductionof small RinvolvedR polymerasecomplexnuclearchromosomc,telomericrcgiontelomerasecatalyticcore complexTERT-RMRPcomplexnucleotide-excision repairfactor1complexchromosome,telomericregionR-directedRpolymerasecomplerERCC4-ERCC1complex0Gene_Number（a）、（c)叶分组非共有突变基因通路及功能富集结果；(b）、（d)侧分组非共有突变基因通路及功能富集结果图3不同原发部位非共有突变基因的功能及通路富集分析Fig.3Function and pathway enrichment analysis of non-shared mutated genes in different primary sites3综合分析及结论目前,非小细胞肺癌仍是全球临床发病率与致死率较高的恶性肿瘤之一,对于肿瘤源性突变特征研究,将有利于进一步深入理解其在发生发展以及治疗响应等临床表现差异的潜在分子机西北大学学报（自然科学版）结果说明【见图3（d)】,左、右侧肺癌组中存在特征基因与凋亡信号通路（intrinsic apoptotic signa-ling pathway、e x t r i n s i c a p o p t o t i c s i g n a l i n g p a t h w a y)、蛋白酶B信号通路、PI3K介导的信号通路、上皮细胞增殖调控等。(b)EGFRtyrosinckinascinhibitorresistance-Central carbon metabolism in cancerEndometrialcancerBrcastcanccrNon-small cell lung cancerColorectalcancerPI3K-Akt signaling pathwayCountPrustatecancer6Gastriccancer7MicroRNAsincancer8Platinumdrug resistance9Pancreaticcancer10Focal adhesion-1112PD-L1cxpression and PD-1 chcckpoint pathway in canccr-Ras signaling pathwayP.adjustCalcium signaling pathway-1e-07Mismatch repair-2e-07MAPK signaling pathway-3e-07FoxO signalingpathway-4e-o07 Sigaling pathways regulating pluripotency f stem clls681012GeneNumber23第53卷Count234567P.adjust-0.002-0.0050.007-0.01234567GeneNumber(d)intrinsic apoptotice signaling pathway inproteinkinase Bsignalingphosphatidylinositol3-kinase signalingphosphatidylinositol-mediated signalinginositol lipid-mediated signalingregulation ofcpithelialcellproliferation-epithelial cellproliferationcxtrinsicapoptotic signalingpathwaycell-matrix adhesionproteinautophosphorylationregulation of protcinkinascBnegative regulation ofD-dependent Dlimbic systemdevelopmentregulation of cell-matrix adhesionintrinsicapoptoticsignaling pathwayurogenital systemdevclopmentnegative regulation of Dreplicationnegativeregulation ofcell-matrixadhesioncell-substrateadhesionpeptidyl-tyrosine phosphorylationpeptidyl-tyrosincmodificationnegativeregulation ofepithelial cell-mismatch repaircomplertransmcmbranereceptorprotcin tyrosinekinase-transmembrane receptor proteinkinase activityprotein tyrosinekinase activitysingle base insertion or deletionErbB-3class receptorbindingbiological processD insertion ordeletion hindingcellularcomponentmismatch repair compler bindingmolecularfunctionmismatchedD bindinggrowthfactorbinding制。研究报道,新辅助治疗后右肺切除术的肺癌患者在手术期死亡率显著高于左侧切除术患者 ，提示肺癌患者全肺切除术后生存期与肺内原发灶起源位置密切相关。此外，原发性肺癌的位置也与特定的淋巴结转移有关。研究发现,与肺中叶、肺下叶非小细胞肺癌患者相比,肺上叶癌患者更易发生纵隔淋巴结转移，且其转移多局限0246Gene_Number第4期于上纵隔区域淋巴结，而肺中叶、肺下叶癌则易出现上、下纵隔跳跃式淋巴结转移 2-3。左、右侧原发性肺癌患者术后治疗响应及预后等也存在显著性差异，原发灶起源部位也能被用于评估肺切除术后辅助治疗响应的预测指标。有研究指出,左侧肺切除术患者癌症特异性生存率与总生存率均明显优于右侧肺切除术患者，且辅助治疗显著延长左侧肺癌切除术后患者生存期，但与右侧肺癌切除术生存期无相关性 4。在ALK 重排的肺癌样本中,右侧肺癌患者样本中检出更高水平的EGFRp.L858R突变。然而，仍有部分研究指出I-IIA期左侧和右侧非小细胞肺癌的预后差异不显著，原发肿瘤分布侧不应作为一种预后因素 12 。不同原发部位肿瘤临床表现差异研究尚无一致性结论。此外,部分学者认为可能由于不同肺叶或左、右侧肺部生理结构及生理功能、血液供给以及肿瘤侵袭途径、基因变异等导致了位置特异性的肿瘤生物学特征变化 2 6.131,但具体机制尚不清楚。因此，本研究通过靶向测序技术对6 2例非小细胞肺癌患者新鲜组织样本进行分析,并通过大数据挖掘50 2 例TCGA非小细胞肺腺癌与之整合分析，从而进一步揭示不同肺叶及分布侧原发性非小细胞肺癌肿瘤源性突变特征及生物学功能影响,为探究非小细胞肺癌致病机制提供新的思路。非小细胞肺癌的发生伴有多种遗传和表观遗传学改变,且其变化特征可能与生理结构或组织学类型上的差异有关。本研究纳人的样本中仅有5例未检出突变（0.8 8%,5/56 4）。不同叶分组与左、右侧肺癌组样本均表现为共有突变位点远低于非共有突变位点，其中共有突变主要来源于TP53、EG FR、K RA S 等基因,与其他学者研究结果相似 14-15。进一步分析显示不同原发部位样本间存在多种特有突变基因（EPCAM、ERCC1等），且高频突变基因如TP53、K RA S等在各组中的突变频率存在差异，但EGFR在TCGA来源的数据集中突变频率（2.2 7%7.2 6%）显著低于新纳入临床标本水平（2 9.2 7%32.8 8%），这可能与数据集样本来源多为非亚裔人群有关，以往研究中发现亚裔肺癌患者中具有更高水平的突变频率 1-17 。KRAS突变的频率因不同种族而异,在以往的报道中非小细胞肺癌患者的检出率为19%30%不等 18 。本研究中KRAS的突变频率仅为10.8%（54/50 2），这与前期研究报道如罗祥，等：基于高通量靶向测序技术的肿瘤原发部位差异研究 625Yang等人也发现患者组（53/10 33,5.13%）和腺癌患者组（45/7 59人，5.93%）的KRAS突变率较低的结论基本一致 7 ,即亚洲人表现出较低的KRAS突变频率。此外,在Liu等人利用二代测序技术对17 113个非小细胞肺癌样本的测序分析中发现2 7%的样本携带KRAS突变,其中最常见突变为p.G12C突变（40%），其次是p.G12V（19%）和p.G12D（15%）19。此结果在本研究中也有相似发现，且深入分析还发现不同叶分组中KRASp.G12C检出率仍为最高，但相较之下，中叶肺癌组中检出率（2 8.57%）明显高于其他两组（8.9 4%和11.7 3%）。然而,左、右侧肺癌组中则分别以ECFRp.L858R、EG FRp.E7 46 _A750del(19号外显子缺失突变）为最多检出的突变类型,这两种突变的检出情况与Yang等人的研究发现携带EGFRp.L858R和19 号外显子缺失突变更容易分布在右侧的结果有所不同 7 。为了评估不同肺叶非小细胞肺癌患者特征突变基因的潜在功能。本研究通过分别筛选出上、中、下叶肺癌组间以及左、右侧肺癌组间的特征突变基因，并对特征突变基因可能参与的信号调控、生物学过程等功能富集分析，研究发现叶分组特征基因主要参与了癌症中心碳代谢以及药物抵抗等通路，而侧分组特征基因主要参与了错配修复及凋亡调控等通路。此外，与以往的研究结果不同 7.2 0 ,本研究中肺内不同原发部位总突变检出量及各基因突变频率分布与不同临床特征（年龄、性别)之间无相关性，提示纳人的患者临床特征差异对本研究分组影响较小。因此，上述结果既能说明病灶起源及研究分组不同均会对结果产生显著影响，又提示原发部位不同的肺癌可能具有不同的基因突变特征。但由于上述特征基因为低频突变基因（伴随突变），其通过哪些调控机制促进不同原发部位肺癌发生发展，目前相关研究报道较少，仍有许多值得研究和思考的问题。综上，本研究采用组织靶基因高通量测序技术与肿瘤大数据挖掘相结合的方法，探讨了56 4例肺内原发部位不同的肺癌患者突变特征及功能学影响,筛选并分析了位置特异性突变候选基因及其可能参与的信号通路及生物学过程,这将对非