-
摘要:目的
探讨蛋白激酶Cι（PKCι）、Yes相关蛋白1（YAP1）与高危型人乳头瘤病毒（HPV）感染在宫颈癌局部免疫微环境中的相关性。
方法选择宫颈正常组织（NCT）、宫颈低级别鳞状上皮内病变（LSIL）、宫颈高级别鳞状上皮内病变（HSIL）以及宫颈鳞状细胞癌（SCC）各80例。收集4组蜡块组织，实时荧光PCR法测定4组组织高危型HPV感染率；免疫组织化学法检测4组组织中PKCι、YAP1、CD4及CD8表达水平并进行相关性以及临床病理分析。
结果NCT、LSIL、HSIL以及SCC各组间高危型HPV感染率、PKCι、YAP1、CD4、CD8阳性率差异均有统计学意义（P < 0.05）；宫颈病变级别与高危型HPV感染、PKCι、YAP1及CD8阳性表达均呈正相关（P < 0.05），与CD4阳性表达呈负相关（P < 0.05）；在SCC中HPV感染、PKCι、YAP1以及CD8阳性表达相互之间的均呈正相关，而与CD4阳性表达均呈负相关（P < 0.05）；HPV感染、PKCι、YAP1及CD8阳性表达在SCC不同分化程度及有无脉管瘤栓之间差异均有统计学意义（P < 0.05）。
结论宫颈病变患者常伴有高危型HPV感染与PKCι、YAP1、CD4及CD8表达异常，相互之间可能存在协同作用，造成SCC的局部免疫抑制微环境，为SCC的发病机制研究及其诊治提供可能的策略。
Abstract:ObjectiveTo investigate the role of PKCι, YAP1 and high-risk HPV infection in the local immune microenvironment of cervical cancer.
MethodsWe chose 80 cases of normal tissue of the cervix (NCT), cervical low-grade squamous intraepithelial lesion (LSIL), cervical high-grade squamous intraepithelial lesion (HSIL) and early cervical squamous cell carcinoma (SCC) each. Four groups were collected.The infection rate of high-risk HPV in four groups was determined by real-time fluorescence PCR method. The expression levels of PKCι, YAP1, CD4 and CD8 in four groups were measured and correlated by IHC and clinicopathologic features were also analyzed.
ResultsThe differences of high-risk HPV infection rate and PKCι, YAP1, CD4, CD8 positive rate among groups of NCT, LSIL, HSIL and SCC had statistical significance (P < 0.05). The level of cervical lesions was positively associated with high-risk HPV infection and positive PKCι, YAP1, CD8 expression (P < 0.05), while negatively associated with positive CD4 expression (P < 0.05). HPV infection and positive PKCι, YAP1, CD8 expression were positively correlated with each other in SCC, while were all negatively correlated with positive CD4 expression(P < 0.05). The differences of HPV infection, PKCι, YAP1 and CD8 positive expression were significant in different levels of differentiation and vascular invasion of SCC (P < 0.05).
ConclusionThe patients with cervical lesions are often accompanied by high-risk HPV infection and abnormal expression of PKCι, YAP1, CD4 and CD8, which may have synergistic effects on each other, causing the local immunosuppression microenvironment of SCC. It provides a possible strategy for the study of pathogenesis, diagnosis and treatment of cervical cancer.
-
Key words:
- Cervical carcinoma /
- HPV /
- Protein kinase Cι /
- Yes-associated protein 1 /
- Immune microenviroment
-
0 引言
骨肉瘤（osteosarcoma）是骨科最常见的恶性肿瘤[1]，其发生、发展是一个复杂的、涉及多基因调控网络失衡导致的细胞异常代谢的结果。基因表达受转录水平及翻译水平的调控，转录因子（transcription factors, TFs）及microRNAs（miRNAs）在调控基因的表达方面扮演了重要的角色。TFs和miRNAs之间存在广泛的相互作用和合作调控，共同组成复杂的调控网络，对于了解细胞的生理过程、生物学功能、疾病的发生机制等都起到了重要作用。目前面临的首要难题是骨肉瘤往往不能做到早期诊断，诊断不及时会耽误治疗时机，严重影响患者预后[2-3]。因此，建立新的、有效的骨肉瘤防治措施势在必行，还需寻找并发现新的功能基因靶向治疗骨肉瘤，根据治疗效果建立分子水平预后指标，开发安全的不同于传统治疗骨肉瘤的放化疗药物，是当前科研及临床研究中的热门话题[4-5]。本研究将利用高通量基因芯片技术分析骨肉瘤中差异表达的基因和miRNAs，结合生物信息学软件预测与调控骨肉瘤相关的TFs、miRNAs及其靶基因，分析调控网络中的关键基因及关键miRNAs的功能，构建骨肉瘤发生发展的“TFs-miRNAs”共调控网络。
1 材料与方法
1.1 miRNA数据来源
收集美国国家生物技术信息中心（National Center of Biotechnology Information, NCBI）维护的GEO数据库中的两组数据库（GPL8227、GPL19631）中所有骨肉瘤miRNA数据，收集数据时间截至2016年12月31日。
1.2 数据库资源
用基因本体论数据库（gene ontology, GO）对所分析得到的差异表达基因进行功能注释；使用京都基因与基因组百科全书（kyoto encyclopedia of genes and genomes, KEGG）数据库对差异表达基因进行代谢通路富集分析；使用转录作用元件数据库（transcriptional regulatory element database, TRED）对差异表达基因进行转录因子筛选。
1.3 方法
1.3.1 差异表达基因筛选及注释
为探究骨肉瘤细胞发生发展的分子机制，用R语言下载GEO数据库中GSE28425、GSE65071数据集的原始数据。使用MAS 5.0算法对数据集中的原始数据进行标准化预处理。使用Limma程序包中所提供的统计分析函数，计算与健康者细胞系相比，骨肉瘤患者中表达水平发生变化的基因。选取差异表达倍数≥2.0（Fold Change≥2.0）为筛选阈值，筛选获得差异表达的miRNA。通过差异表达分析获得两组平台共有的骨肉瘤细胞系表达水平发生显著性改变的miRNA，利用DIANA TOOLS数据库中收录的人类基因信息，找到具有显著差异的miRNA靶基因。利用R语言，将筛选得到的差异表达基因与GEO数据库提供的GPL平台注释信息进行匹配，对差异表达基因进行信息注释。
1.3.2 差异表达基因KEGG代谢通路富集分析
为进一步探究差异表达基因参与生物代谢的代谢通路，对筛选得到的差异表达基因进行KEGG代谢通路富集分析。提取差异表达基因GPL平台注释信息中的“Entrez_Gene_ID”信息，应用R语言KEGG程序包与KEGG数据库信息，对每一个差异表达基因向已知的代谢通路进行富集映射，探究差异表达基因集中映射的细胞代谢通路，初步推测骨肉瘤所影响的主要代谢通路，并筛选出与肿瘤通路关系密切的基因，进行GO功能注释分析。
1.3.3 差异表达基因GO功能注释分析
为探究差异表达且与肿瘤通路关系密切的基因所参与的生物功能，对筛选得到的与肿瘤通路关系密切的差异表达基因进行GO功能注释分析。提取差异表达基因GPL平台注释信息中的“Entrez_Gene_ID”信息，使用R语言GO程序包与GO数据库信息，对每一个差异表达基因进行功能注释，探究骨肉瘤组织与正常组织差别的基因功能，初步推测骨肉瘤作用的分子机制。
1.3.4 分子共调控网络构建
为进一步挖掘骨肉瘤的分子调控模式，探究差异表达基因之间的相互调控作用关系，提取TRED数据库中所收录的所有人类转录因子数据。通过Pearson相关性分析，提取基因间分子共调控作用关系。使用Cytoscape version 3.4.0软件对得到的分子共调控作用关系进行可视化展示，构建出最终的分子共调控网络。
2 结果
2.1 差异表达基因的筛选
本研究分别从GPL8227、GPL19631两组数据库中筛选出差异表达miRNA 232个（上调129个、下调103个）和265个（上调205个、下调60个）。共筛选得到两组均包含的差异表达miRNA52个，其中31个miRNA上调表达，21个miRNA下调表达见表 1。
表 1 差异表达的miRNAs信息Table 1 Information of differentially-expressed miRNAs
2.2 差异表达基因KEGG代谢通路富集分析
将筛选出的52个差异表达miRNA结合DIANA TOOLS数据库中收录的人类基因信息进行分析。其中KEGG代谢通路富集分析结果见图 1，差异表达基因主要参与肿瘤相关的诸多代谢通路。差异表达的52个miRNA的整体调控模式及每个miRNA所调控的代谢通路情况见图 2。肿瘤代谢通路的分子调控关系网络见图 3。
![]()
图 1 差异表达基因KEGG代谢通路富集分析Figure 1 KEGG metabolic pathway enrichment analysis of differentially-expressed genes![]()
图 2 差异表达52个miRNA聚类分析热图Figure 2 Heat map of clustering analysis of 52 differentially-expressed miRNAs![]()
图 3 差异表达miRNA靶基因肿瘤代谢通路Figure 3 Tumor metabolic pathway of differentially-expressed miRNA target genes2.3 差异表达基因GO功能注释分析
将上述筛选出与肿瘤通路关系最紧密的5个miRNA及其靶基因进行GO功能注释分析。注释分析结果显示这些基因主要参与维生素E的代谢路径（response to vitamin E）、RNA聚合酶Ⅱ启动子的转录调控反应（regulation of transcription from RNA polymerase Ⅱ promoter in response etc）、调控神经元的成熟（regulation of neuron maturation）、受体激动剂活性（receptor agonist activity）、正调控血管内皮生长因子受体活性（positive regulation of vascular endothelial growth factor receptor activity）、通过MAPK级联非经典Wnt信号途径（non-canonical wnt signaling pathway via MAPK cascade）、间充质上皮细胞信号转导（mesenchymal-epithelial cell signaling）、Ⅳ型胶原（collagen type Ⅳ）、细胞对生长激素刺激的反应（cellular response to growth hormone stimulus）、蛋白激酶活性的激活（activation of protein kinase A activity）等生物功能，见图 4。
![]()
图 4 差异表达基因GO功能注释分析Figure 4 GO functional annotation analysis of differentially-expressed genes2.4 分子共调控网络构建
最后基于与肿瘤通路密切相关的差异表达基因及TRED数据库中所收录的人类转录因子信息，对差异表达基因进行分子偶联，构建了分子共调控网络，见图 5。
3 讨论
目前miRNA被广泛认为是血液学和肿瘤组织的生物标志物。miRNA有可能通过扭曲、突变、成纤维细胞生长因子（FGF）和骨形态发生蛋白（BMP）信号参与骨的形成。研究发现miRNA与多种骨科疾病有关系，比如有研究认为miRNA-140与骨性关节炎有关，miRNA-503、miRNA-3960与骨质疏松有关，miRNA-34、miRNA-146等与类风湿关节炎有关[6]。研究表明众多的miRNA在骨肉瘤发病的机制中异常表达，被miRNA决定的许多基因被上调或下调。比如Karbasy等[7]通过研究发现miR-125b和miR-300的下调与骨肉瘤的进展相关且两者都起到抑癌基因的作用。Goudarzi等[8]通过研究发现miR-185上调miR-218下调与软骨肉瘤的进展相关联，并且两者均可以作为软骨肉瘤的抑癌基因。Bahador等[9]通过研究发现在青少年、儿童以及青年人骨肉瘤患者中，如果miR-29b和miR-422a的表达下调往往表明患者的预后较差。Taheriazam等[10]通过研究发现下调microRNA-26a和微小RNA-27A的上调有助于判断人类骨肉瘤的攻击性进展。高表达的miR-130b水平和miR-218的低表达与差的临床病理特征有关。此外，有研究表明miR-130b可对预测骨肉瘤的进展起关键作用[10-11]。另有研究表明miR-145、miR-128、miR-34a在骨肉瘤中的表达均发生了改变。miR-335、miR-128、miR-340、miR-214等均与骨肉瘤的侵袭性有重要关系，表达增高后预示骨肉瘤的恶性程度较高，可以较早发生转移，患者的生存期缩短[12]。并且近期有研究表明血浆中miR-34b表达下降可能与骨肉瘤已经发生转移有密切的联系[13]，如果能明确这种作用关系会对骨肉瘤的早期诊断起到非凡的作用。就目前研究现状来看骨肉瘤发生发展的分子作用机制尚不十分明确，本研究利用高通量基因芯片技术，通过生物信息学手段，系统性预测构建了一张与肿瘤通路关系密切的TF、miRNA及其靶基因共调控网络蓝图。该分子共调控网络从整体的角度鸟瞰与肿瘤通路关系密切的miRNA及其靶基因的相互作用关系，展示了差异表达基因主要参与肿瘤相关的诸多代谢通路及与肿瘤通路有关的基因主要参与的生物功能。为骨肉瘤诊断标志物簇的筛选提供可靠的选择方向；为临床诊断及预后分析分子靶点的筛选提供实用工具，同时也为其他癌症的同类研究提供新思路。
Competing interests: The authors declare that they have no competing interests.作者贡献陈胜民：资料收集、实验操作、数据统计及论文撰写王国平：资料收集、实验操作及数据统计林志仁：实验操作及结果判读袁峰：病例采集、实验操作、结果判读及论文撰写 -
![]()
图 1 PKCι(A)、YAP1(B)、CD4(C)及CD8(D)在SCC中的表达
Figure 1 Expression of PKCι(A), YAP1(B), CD4(C) and CD8(D) in SCC tissues
表 1 各组高危型HPV感染率以及PKCι、YAP1、CD4、CD8阳性率比较(n(%))
Table 1 Comparison of HPV infection rate, PKCι, YAP1, CD4 and CD8 positive rates among groups (n(%))
下载: 导出CSV
表 2 SCC患者HPV感染、PKCι、YAP1、CD4及CD8表达相关性
Table 2 Correlation among HPV infection, PKCι, YAP1, CD4 and CD8 expression in SCC patients
下载: 导出CSV
表 3 HPV感染以及各蛋白表达与早期SCC临床病理特征的关系(n)
Table 3 Relation between HPV infection, proteins expression and clinicopathological features of early SCC (n)
下载: 导出CSV
-
[1] 赵婧, 钱延玲. 宫颈局部免疫微环境与宫颈癌[J]. 世界最新医学信息文摘, 2018, 18(43): 61-63. https://www.cnki.com.cn/Article/CJFDTOTAL-WMIA201843024.htm Zhao J, Qian YL. The Relationship between Immune Microenvironment and Cervical Cancer[J]. Shi Jie Zui Xin Yi Xue Xin Xi Wen Zhai, 2018, 18(43): 61-63. https://www.cnki.com.cn/Article/CJFDTOTAL-WMIA201843024.htm
[2] Macartney K, Phillips A, Patel C, et al. Authors'reply: Safety of Human Papillomavirus Vaccines[J]. Drug Saf, 2018: 41(5): 541-543. doi: 10.1007/s40264-018-0655-1
[3] Sharmistha S, Christopher AB, Prasenjit D, et al. PRKCI promotes immune suppression in ovarian cancer[J]. Genes Dev, 2017, 31(11): 1109-1121. doi: 10.1101/gad.296640.117
[4] 陈胜民, 袁峰. 蛋白激酶Cι/Yes相关蛋白1与人乳头瘤病毒在肿瘤微环境中的作用研究进展[J]. 广西医学, 2020, 42(24): 3224-3227. https://www.cnki.com.cn/Article/CJFDTOTAL-GYYX202024026.htm Chen SM, Yuan F. Progress in the role of PKCι/YAP1 and HPV in the immune microenvironment of tumor[J]. Guangxi Yi Xue, 2020, 42(24): 3224-3227. https://www.cnki.com.cn/Article/CJFDTOTAL-GYYX202024026.htm
[5] 陈胜民. PTEN、PI3K和Akt蛋白表达与卵巢浆液性癌临床病理学的相关性研究[D]. 石河子: 石河子大学, 2011: 4-5. Chen SM. Study on correlation of clinical pathology between PTEN and PI3K/Akt protein expression in ovarian serous tumors[D]. Shihezi: Shihezi University, 2011: 4-5.
[6] Sichero L, Rollison DE, Amorrortu RP, et al. Beta Human Papillomavirus and Associated Diseases[J]. Acta Cytol, 2019, 63(2): 100-108. doi: 10.1159/000492659
[7] Gilkey MB, Margolis MA. Human Papillomavirus Vaccination Coverage Gaps in Young Adolescents[J]. J Infect Dis, 2019, 220(5): 727-729. doi: 10.1093/infdis/jiy683
[8] Hanbazazh M, Gyure KA. Ocular Human Papillomavirus Infections[J]. Arch Pathol Lab Med, 2018, 142(6): 706-710. doi: 10.5858/arpa.2017-0571-RA
[9] Harden ME, Munger K. Human papillomavirus molecular biology[J]. Mutat Res Rev Mutat Res, 2017, 772: 3-12. doi: 10.1016/j.mrrev.2016.07.002
[10] Schiffman M, Doorbar J, Wentzensen N, et al. Carcinogenic human papillomavirus infection[J]. Nat Rev Dis Primers, 2016, 2: 16086. doi: 10.1038/nrdp.2016.86
[11] Westricha JA, Warrena CJ, Pyeona D. Evasion of host immune defenses by human papillomavirus[J]. Virus Res, 2017, 231: 21-33. doi: 10.1016/j.virusres.2016.11.023
[12] Pearsall EA, Connor KM. Atypical Protein Kinase C: Breaking Down Barriers in Ocular Disease?[J]. Am J Pathol, 2018, 188(10): 2142-2146. doi: 10.1016/j.ajpath.2018.07.006
[13] Parker PJ, Justilien V, Riou P, et al. Atypical Protein Kinase Cι as a human oncogene and therapeutic target[J]. Biochem Pharmacol, 2014, 88(1): 1-11. doi: 10.1016/j.bcp.2013.10.023
[14] 葛彩云, 蔡红兵, 宋紫烨, 等. HPV16-E6上调MIF对宫颈癌细胞增殖和凋亡的影响[J]. 肿瘤防治研究, 2018, 45(7): 468-474. doi: 10.3971/j.issn.1000-8578.2018.17.1406 Ge CY, Cai HB, Song ZY, et al. HPV16-E6 Influences Growth and Apoptosis of Cervical Cancer Cells Through Upregulation of MIF[J]. Zhong Liu Fang Zhi Yan Jiu, 2018, 45(7): 468-474. doi: 10.3971/j.issn.1000-8578.2018.17.1406
-
期刊类型引用(1)
1. 冯晓慧，朱正秋. NLR对晚期食管鳞癌一线免疫治疗疗效的预测价值. 医学研究杂志. 2023(07): 157-160+123 . 
百度学术
其他类型引用(2)
计量
- 文章访问数: 1249
- HTML全文浏览量: 390
- PDF下载量: 349
- 被引次数: 3
