Chinese Journal of Lung Cancer
中国肺癌杂志编辑部
image
利用高通量转录组测序对肺鳞癌基因型的
Volume: 20, Issue: 11
DOI 10.3779/j.issn.1009-3419.2017.11.01
  • PDF   
  • XML   
  •       

Table of Contents

Highlights

Notes

Abstract

背景与目的目前尚无有效的靶向药物用于肺鳞癌的临床治疗,在肺腺癌中有效的靶向药物却无法让鳞癌患者获益,而研究者对肺鳞癌的靶点知之甚少。本研究重点筛选并鉴定肺鳞癌的特异关键基因,为肺鳞癌的治疗提供新的靶点。方法转录组测序检测5例肺鳞癌患者的配对肺癌组织及正常肺组织样本,使用生物信息学分析筛选两组之间的差异编码基因,进而应用实时定量PCR验证关键差异基因在肺癌细胞系的表达情况。结果转录组测序结果显示,与正常肺组织相比,肺鳞癌组织上调的差异基因数目为534条。位于前十位肺癌组织高表达的基因包括GAGE12JSPRR3PRAMESPRR1ASPRR2EMAGEA3SPRR1BIL36GTMPRSS11DSPRR2D。进而于不同生物特征的肺癌细胞系H520、GLC82、A549、H1299及PC9验证SPRR家族表达状态,我们发现在淋巴结转移细胞系H1299处于高表达水平。结论高通量转录组测序筛选到肺鳞癌异常高表达基因SPRR家族,此家族与肺癌淋巴结转移相关,为肺癌的靶向治疗提供了新的思路。

Keywords
利用高通量转录组测序对肺鳞癌基因型的

目前肺癌的发病率和死亡率在全球范围内位于前列,对公众健康造成极大危害。肺癌的死亡人数在过去的三十年上升了465%,约增长了5倍,已经成为我们国家恶性肿瘤死亡原因中第一位。所以现在人们把肺癌称为我们国家癌症领域的“第一杀手”。其中非小细胞肺癌(non-small cell lung cancer, NSCLC)的发病率占到了85%,主要包括腺癌、鳞癌以及大细胞癌,其中肺鳞癌是最常见的病理类型之一,占NSCLC的30%-40%[1, 2]。肺鳞癌又称肺鳞状上皮细胞癌,以中央型肺癌多见,并有胸管腔内生长的倾向,早期常引发支气管狭窄或阻塞性肺炎[2]

目前分子靶向治疗的出现彻底改变了癌症治疗格局。在肺腺癌中,分子靶向治疗获得了较为广泛的应用,但是在肺鳞状细胞癌患者的治疗中还没有实现。高通量二代测序作为一种新的实验技术,能够为我们提供详细、特异的信息,最重要的是它的正确率和高特异性[3, 4]。所以从现实意义的角度出发,本研究的目的在于利用高通量二代测序技术,通过查找正常组织与肿瘤组织在基因表达上的差异,为实现肺鳞癌的分子靶向治疗提供基础。因此,我们针对5对肺癌鳞癌患者的肿瘤组织和正常肺组织进行高通量转录组测序,发现差异表达基因,并且对其差异基因SPRR进行肺癌细胞系的验证,有望为肺鳞癌靶向治疗提供新的靶点。

1. 资料与方法

1.1. 临床资料

收集5例经病理诊断为肺鳞癌患者的肿瘤组织及正常肺组织的新鲜标本,年龄45岁-65岁。标本收集后分立即置于液氮中速冻,保存于-80 ℃冰箱中,用于实验检测。肺癌细胞系H520、GLC82、A549、H1299及PC9使用1640培养基、37 ℃的CO2温箱培养。

1.2. 总RNA提取

取等量肺癌癌组织及正常肺组织,液氮中磨碎;加入TRIZOL Reagent,室温孵育5 min;加入200 μL氯仿剧烈振荡15 s,室温孵育10 min;4 ℃,12, 000 rpm,离心15 min。吸取上层含有RNA的水相入新管,加500 μL异丙醇沉淀RNA,室温孵育10 min;4 ℃,12, 000 rpm,离心15 min。弃上清,加1 mL的75%乙醇洗涤RNA沉淀,4 ℃,7, 500 rpm,离心5 min,用20 pL DEPC水溶解RNA。用紫外分光光度计测定RNA的含量,-80 ℃冰箱保存。

1.3. 转录组文库制备及上机测序

利用Oligo dT微珠富集纯化mRNA并进行片段化处理,在3末端加碱基A、连接测序接头,割胶,纯化并回收200 bp-500 bp之间的cDNA片段。PCR扩增后完成测序文库的制备。对测序文库进行质量控制,构件好的文库用Illumina Hiseq2500进行测序,以fastq格式输出。

1.4. 生物信息分析

将转录组测序样本,分别用TopHat回贴人基因组,其样本比对率需达到90%以上。用Cufflinks分别重新构建转录本,对转录组进行鉴定,合并转录本,得到一个统一的转录本集合。对转录本进行鉴定质量的评估后,最后采用cuffdiff鉴定正常组织样本组和相对应的肺癌组织样本组之间的差异编码基因。

1.5. 实时定量PCR

使用TRIzol提取肺癌细胞系总RNA,2 μg总RNA逆转成cDNA,使用SYBR Green PCR试剂盒,罗氏480仪器检测差异基因的表达。反应条件,95 ℃ 5 min;95 ℃ 15 s,60 ℃ 30 s,40个循环;72 ℃ 5 min。β-actin作为对照,计算公式:2-ΔCt,ΔCt=Ct目的基因-Ct对照

2. 结果

2.1. 肺鳞癌转录组概况

利用高通量转录组测序技术,检测5对肺癌组织及正常肺组织的标本,利用生物信息学分析筛选差异基因。用cuffdiff鉴定正常组织样本组和相对应的肺癌组织样本组之间的差异基因,发现按照cuffdiff的标准(上下调2倍;p-value≤0.05;q-value≤0.05),鉴定出来40个差异基因。不考虑q-value之后,差异基因数目明显上升,为534个(q-value是基于p-value的再统计,p-value是样本的一个基本检验几率)。

2.2. 筛选肺鳞癌组织高表达的前10个基因

基于转录组测序检测5对配对肺癌组织及正常肺组织的标本,重点关注于肺癌肿瘤组织高表达的前10个关键基因(表 1)。正常肺组织测序的FPKM平均值为0.222, 4,肺癌肿瘤组织FPKM平均值为233.4,则肿瘤组织前10位差异基因的表达量是正常组织的1, 049倍。

1
肺鳞癌组织中高表达的前10个基因The top ten genes highly expressed in the lung squamous carcinoma
GenesNormal tissues FPKMTumor tissues FPKMLog2 (Fold change)P
GAGE12J0.018, 025, 354.51911.562, 50.006, 067, 27
SPRR30.198, 758307.33610.594, 69.24E-06
PRAME0.008, 184, 3711.32910.434, 90.039, 201, 7
SPRR1A0.515, 89575.83810.124, 40.010, 568, 3
SPRR2E0.125, 185132.15610.0445.65E-07
MAGEA30.006, 615, 166.476, 359.935, 190.008, 148, 43
SPRR1B0.511, 568493.1499.912, 888.55E-07
IL36G0.035, 212, 133.301, 59.885, 30.025, 304, 8
TMPRSS11D0.005, 171, 984.680, 649.821, 780.000, 461, 89
SPRR2D0.799, 756715.0349.804, 246.16E-06

在这10个基因中,SPRR1ASPRR1BSPRR2DSPRR2ESPRR3为同一家族成员,其FPKM值在正常肺组织与肿瘤组织中分别为0.430, 2和444.72,二者相差1, 033倍。其中又以SPRR3在肺肿瘤组织与正常组织之间的差异最大,为1, 546倍(表 2)。

2
SPRRs家族基因在正常肺组织与肺癌肿瘤组织中的表达Expression of SPRRs family in the normal tissues and lung cancer tissues
GenesNormal tissues
FPKM
Tumor tissues
FPKM
Fold change
SPRR30.198, 758307.3361, 546
SPRR1A0.515, 89575.8381, 116
SPRR2E0.125, 185132.1561, 055
SPRR1B0.511, 568493.149963
SPRR2D0.799, 756715.034894

2.3. 肺癌细胞系中SPRR家族基因的表达情况

采用荧光定量PCR的方法检测了前10位差异基因在肺癌细胞系中的表达情况,见图 1。H520为肺鳞癌细胞,GLC82、A549、H1299、PC9为肺腺癌细胞。其中H1299为淋巴结转移的肺癌细胞系,PC9是表皮生长因子受体酪氨酸激酶抑制剂(epithelial growth factor receptor tyrosine kinase inhibitors, EGFR-TKIs)敏感的肺癌细胞系。从图 1中可以看出SPRR1A、SPRR1B、SPRR2D、SPRR2E和SPRR3在肺癌细胞中广泛表达,此5类亚型于淋巴结转移的H1299细胞中均有较高表达,均高于在未转移的肺腺癌细胞系A549的表达水平。

柱状图显示SPRRs家族基因包括SPRR1A、SPRR1B、SPRR2D、SPRR2E和SPRR3在肺癌细胞系中的表达情况(2-ΔCt value)。
1
柱状图显示SPRRs家族基因包括SPRR1ASPRR1BSPRR2DSPRR2ESPRR3在肺癌细胞系中的表达情况(2-ΔCt value)。The histogram shows expression of SPRRs family genes including SPRR1A, SPRR1B, SPRR2D, SPRR2E and SPRR3 in the lung cancer cells lines (2-ΔCt value).

3. 讨论

近十年来,基于NSCLC驱动基因的个体化靶向诊治研究拉开了帷幕,并取得了革命性进展而应用于临床,例如EGFR突变基础上的TKIs,基于EML4-ALK融合基因的抑制剂等主要应用于肺腺癌[5]。然而,针对肺鳞癌的靶向诊治研究十分有限。此外,随着新一代高通量测序技术的快速发展,转录组测序已成为基因表达和转录组分析新的重要手段。转录组测序能够快速并全面的获取组织中几乎所有转录本RNA的序列信息,是目前常用的高通量分析手段,已成为肿瘤研究的重要技术支撑。

本研究正是利用转录组测序技术检测了5对肺癌组织及正常肺组织的基因型的改变,并利用生物信息学分析筛选出与肿瘤密切相关的排在前10位的基因,包括SPRR1ASPRR1BSPRR2DSPRR2ESPRR3PRAMEMAGEA3IL36GGAGE12JTMPRSS11D。其中,位于差异基因前5位的SPRR1ASPRR1BSPRR2DSPRR2ESPRR3基因同属于SPRRs家族,编码一类富含脯氨酸的蛋白[6],而有研究表明这些基因亚型在表达上有一定的偏向性。SPRR1ASPRR1B基因表达与口腔上皮细胞的屏障作用相关[7]SPRR2基因亚型在消化道疾病以及皮肤病的炎症反应中表达水平有所增加[8, 9]

IL36G基因在组织发生损伤与修复过程中表达增高,可见与炎症和免疫反应密切相关[10, 11]。也有报道显示IL36G基因与皮肤病的炎性反应相关[12]。据报道炎症反应能够促进肿瘤的发生,并与肿瘤发生转移密切相关。而本研究结果显示SPRRs家族及IL36G在癌组织高表达的情况可能参与了炎症反应且具有促进肺癌发生的作用。GAGE12J基因在胎儿和肿瘤组织中表达,是肿瘤/睾丸抗原家族成员之一。有报道[13, 14]发现其在脑膜瘤和神经鞘瘤中可与其他基因相互作用。TMPRSS11D属于Ⅱ型跨膜丝氨酸蛋白酶亚家族成员,在宫颈鳞癌和食管鳞癌的形成过程中缺失[15]。PRAME是肿瘤/睾丸抗原的一种,可作为潜在的免疫治疗的靶点,其表达状态与上皮性卵巢癌的启动子低甲基化有关[16]。MAGEA3同样也属于肿瘤/睾丸抗原的一种,ORIS与其启动子结合而调控肺癌的转录活性[17]

我们进一步在肺癌细胞系中检测差异基因SPRRs家族的表达状态发现,在淋巴结转移的细胞系H1299中,SPRR1ASPRR1BSPRR2DSPRR2ESPRR3均呈高表达状态。前期的研究发现SPRRs家族基因在胆管上皮细胞中能够增加对损伤的抵抗,并与EMT密切相关[18]。另有研究[19]表明长链非编码RNA MALAT-1通过调控SPRRs家族基因而影响了舌鳞状细胞癌的远端转移。而我们的研究亦发现此家族基因在局部转移的细胞中高表达,提示有可能在肺癌细胞发生转移过程中发挥着重要作用。

目前随着临床研究的深入,证明了个体化的肿瘤治疗相比传统的方法更加安全有效。而本研究利用转录组测序技术,通过查找正常肺组织与肺癌肿瘤组织在基因表达上的差异,为实现肺鳞癌的分子靶向治疗以及寻找更多的分子治疗靶点提供了可行性。

References

1 

    Kobayashi-Watanabe N, Sato A, Watanabe T, . Functional analysis of discoidin domain receptor 2 mutation and expression in squamous cell lung cancer. Lung Cancer 2017. 110: , pp.35-41, doi: 10.1016/j.lungcan.2017.05.017

2 

    Perez-Moreno P, Brambilla E, Thomas R, . Squamous cell carcinoma of the lung: molecular subtypes and therapeutic opportunities. Clin Cancer Res 2012. 18: 9, pp.2443-2451, doi: 10.1158/1078-0432.CCR-11-2370

3 

4 

    D'Haene N, Le Mercier M, De Neve N, . Clinical validation of targeted next generation sequencing for colon and lung cancers. PLoS One 2015. 10: 9, pp.e0138245, doi: 10.1371/journal.pone.0138245

5 

    Hu W, Liu Y, Chen J. Concurrent gene alterations with EGFR mutation and treatment efficacy of EGFR-TKIs in Chinese patients with non-small cell lung cancer. Oncotarget 2017. 8: 15, pp.25046-25054http://www.ncbi.nlm.nih.gov/pubmed/28212572

6 

    Gibbs S, Fijneman R, Wiegant J, . Molecular characterization and evolution of the SPRR family of keratinocyte differentiation markers encoding small proline-rich proteins. Genomics 1993. 16: 3, pp.630-637, doi: 10.1006/geno.1993.1240

7 

    Lee CH, Marekov LN, Kim S, . Small proline-rich protein 1 is the major component of the cell envelope of normal human oral keratinocytes. FEBS Lett 2000. 477: 3, pp.268-272, doi: 10.1016/S0014-5793(00)01806-8

8 

    Steinert PM, Candi E, Tarcsa E, . Transglutaminase crosslinking and structural studies of the human small proline rich 3 protein. Cell Death Differ 1999. 6: 9, pp.916-930, doi: 10.1038/sj.cdd.4400568

9 

    Tarcsa E, Candi E, Kartasova T, . Structural and transglutaminase substrate properties of the small proline-rich 2 family of cornified cell envelope proteins. J Biol Chem 1998. 273: 36, pp.23297-23303, doi: 10.1074/jbc.273.36.23297

10 

    Birse KD, Romas LM, Guthrie BL, . Genital injury signatures and microbiome alterations associated with depot medroxyprogesterone acetate usage and intravaginal drying practices. J Infect Dis 2017. 215: 4, pp.590-598http://www.ncbi.nlm.nih.gov/pubmed/28011908

11 

    Keermann M, Koks S, Reimann E, . Transcriptional landscape of psoriasis identifies the involvement of IL36 and IL36RN. BMC Genomics 2015. 16: , pp.322, doi: 10.1186/s12864-015-1508-2

12 

    Li N, Yamasaki K, Saito R, . Alarmin function of cathelicidin antimicrobial peptide LL37 through IL-36gamma induction in human epidermal keratinocytes. J Immunol 2014. 193: 10, pp.5140-5148, doi: 10.4049/jimmunol.1302574

13 

    Gjerstorff MF, Ditzel HJ. An overview of the GAGE cancer/testis antigen family with the inclusion of newly identified members. Tissue Antigens 2008. 71: 3, pp.187-192, doi: 10.1111/tan.2008.71.issue-3

14 

    Torres-Martin M, Lassaletta L, Isla A, . Global expression profile in low grade meningiomas and schwannomas shows upregulation of PDGFD, CDH1 and SLIT2 compared to their healthy tissue. Oncol Rep 2014. 32: 6, pp.2327-2334, doi: 10.3892/or.2014.3526

15 

    Duhaime MJ, Page KO, Varela FA, . Cell surface human airway trypsin-like protease is lost during squamous cell carcinogenesis. J Cell Physiol 2016. 231: 7, pp.1476-1483, doi: 10.1002/jcp.25173

16 

17 

    Bhan S, Negi SS, Shao C, . BORIS binding to the promoters of cancer testis antigens, MAGEA2, MAGEA3, and MAGEA4, is associated with their transcriptional activation in lung cancer. Clin Cancer Res 2011. 17: 13, pp.4267-4276, doi: 10.1158/1078-0432.CCR-11-0653

18 

    Demetris AJ, Specht S, Nozaki I, . Small proline-rich proteins (SPRR) function as SH3 domain ligands, increase resistance to injury and are associated with epithelial-mesenchymal transition (EMT) in cholangiocytes. J Hepatol 2008. 48: 2, pp.276-288, doi: 10.1016/j.jhep.2007.09.019

19 

    Fang Z, Zhang S, Wang Y, . Long non-coding RNA MALAT-1 modulates metastatic potential of tongue squamous cell carcinomas partially through the regulation of small proline rich proteins. BMC Cancer 2016. 16: , pp.706, doi: 10.1186/s12885-016-2735-x
https://www.researchpad.co/tools/openurl?pubtype=article&doi=10.3779/j.issn.1009-3419.2017.11.01&title=利用高通量转录组测序对肺鳞癌基因型的&author=&keyword=肺肿瘤,高通量测序,SPRR家族,Lung neoplasms,Transcriptome sequence,SPRR family,&subject=临床研究,Clinical Research,