喉鳞状细胞癌的转录组数据分析及生物标志物筛选

671
《癌症进展》2021年4月第19卷第7期
*论著*
ONCOLOGY PROGRESS,Apr2021V ol.19,No.7喉鳞状细胞癌的转录组数据分析及生物标志物筛选
王娜娜1，杨川2,3#
1延安大学附属医院耳鼻咽喉科，陕西延安716000
2河北医科大学第二医院耳鼻咽喉科，石家庄050000
3平昌县人民医院耳鼻咽喉科，四川巴中636040
摘要
摘要：：目的通过数据挖掘分析长链非编码RNA（lncRNA）、微小RNA（miRNA）和信使RNA（mRNA）在喉鳞状细胞癌中的相互作用机制。

方法从肿瘤基因图谱（TCGA）中获取喉鳞状细胞癌的lncRNA、miRNA、mRNA表达谱数据及相关临床数据，从TCGA和GEO数据库中获取喉鳞状细胞癌的DNA甲基化数据。

进行蛋白相互作用（PPI）、基因文本（GO）、京都基因和基因组百科全书（KEGG）信号通路分析和Kaplan-Meier生存分析，构建竞争性内源RNA（ceRNA）网络和转录因子调控网络。

结果LINC00278、LINC00689、MYHA S、miRNA-99a和miRNA-301a属于低风险基因（HR﹤1）；MNX1-AS1、LINC02575、HOXB-AS4、LSAMP-AS1、LINC02086、IG-FL2-AS1、LINC02253、CASC20、miRNA-383、miRNA-196a-2、miRNA-196a-1、miRNA-100和miRNA-4652属于高风险基因（HR﹥1）。

ceRNA网络中有11个lncRNA（LINC02576、LINC02086、AC020659.1、LINC00528、LINC00689、HOXB-AS4、MNX1-AS1、LINC00278、AC010624.1、AC016773.1、MYHAS）、5个miRNA（has-miRNA-206、has-miR-NA-573、has-miRNA-3662、has-miRNA-133b、has-miRNA-449a）和12个mRNA（STC2、P AX3、NETO2、EIF5A2、AD-PRHL1、SYNM、ACTA1、GPR156、CLIC5、BMP3、HNF4A、FOXL2）。

MYHAS的共表达mRNA（cor-mRNA）主要富集在钙信号通路、环磷酸鸟苷（cGMP）-cGMP依赖性蛋白激酶（PKG）信号通路和催产素信号通路等。

转录因子网络中有9个转录因子（MYOD、RSRFC4、TAL1BETAITF2、YY1、P300、PAX3、PAX5、ZID和OLF1）。

喉鳞状细胞癌DNA甲基化分析得到75个甲基化位点，对应高甲基化基因56个，低甲基化基因15个。

结论ceRNA网络并不能完全阐释lncRNA、miRNA和mRNA在喉鳞状细胞癌发生中的作用机制，lncRNA、miRNA和mRNA在喉鳞状细胞癌发生过程中既相互联系又相互独立。

筛选出的lncRNA、miRNA和mRNA为试验验证缩小了范围并提供了理论基础。

MYHAS在喉鳞状细胞癌中发挥重要作用，是一个潜在的生物学靶点。

：喉鳞状细胞癌；长链非编码RNA；竞争性内源RNA；转录因子网络；DNA甲基化
关键词：
关键词
：A doi：10.11877/j.issn.1672-1535.2021.19.07.06
文献标志码：
739..65文献标志码
中图分类号：
：R739
中图分类号
Transcriptome data analyses and biomarker identification in laryngeal squamous cell carcinoma
WANG Nana1,YANG Chuan2,3#
1Department of Otorhinolaryngology,Yan’an University Affiliated Hospital,Yan’an716000,Shaanxi,China
2Department of Otorhinolaryngology,the Second Hospital of Hebei Medical University,Shijiazhuang050000,Hebei,China
3Department of Otorhinolaryngology,Pingchang County People’s Hospital,Bazhong636040,Sichuan,China
Abstract:Objective To analyze the interaction mechanism of long noncoding RNA(lncRNA),microRNA(miR-NA)and messenger RNA(mRNA)in laryngeal squamous cell carcinoma(LSCC)by using data mining.Method The transcriptome data of lncRNA,miRNA and mRNA in LSCC from TCGA database and related clinical data of were ob-tained;data of DNA methylation in LSCC was obtained from the TCGA and GEO database;the signaling pathway analy-sis of protein-protein interaction(PPI),gene ontology(GO),Kyoto Encyclopedia of Genes and Genomes(KEGG),and Kaplan-Meier survival analysis were performed to construct the competitive endogenous RNA(ceRNA)and transcription factor regulation network.Result LINC00278,LINC00689,MYHAS,miRNA-99a and miRNA-301a were low risk genes(HR<1);MNX1-AS1,LINC02575,HOXB-AS4,LSAMP-AS1,LINC02086,IGFL2-AS1,LINC02253,CASC20,miR-NA-383,miRNA-196a-2,miRNA-196a-1,miRNA-100and miRNA-4652were high risk genes(HR>1).The ceRNA net-work included11lncRNA(LINC02576,LINC02086,AC020659.1,LINC00528,LINC00689,HOXB-AS4,MNX1-AS1, LINC00278,AC010624.1,AC016773.1and MYHAS),5miRNA(has-miRNA-206,has-miRNA-573,has-miRNA-3662, has-miRNA-133b and has-miRNA-449a)and12mRNA(STC2,P AX3,NETO2,EIF5A2,ADPRHL1,SYNM,ACTA1, GPR156,CLIC5,BMP3,HNF4A,FOXL2);the co-expression mRNA(cor-mRNA)of MYHAS was chiefly enriched in the calcium signaling pathway,cyclic guanosine monophosphate(cGMP)-cGMP-dependent protein kinase(PKG)signal-
《癌症进展》2021年4月第19卷第7期
ing pathway and the oxytocin signaling pathway;there were9transcription factors in the transcription factor network (MYOD,RSRFC4,TAL1BETAITF2,YY1,P300,PAX3,PAX5,ZID and OLF1);a total of75methylation loci were ob-tained for the DNA methylation in LSCC,corresponding to56high methylation genes and16low methylation genes. Conclusion The ceRNA network cannot adequately explain the interaction mechanism of lncRNA,miRNA and mRNA
in LSCC.The lncRNA,miRNA and mRNA are both interconnected and independent.The lncRNA,miRNA and mRNA screened decrease the scope of testing identification and provide a theory basis.The MYHAS plays important role in the LSCC,which may be a potential biomarker.
Key words:laryngeal squamous cell carcinoma;long noncoding RNA;competitive endogenous RNA;transcription factor network;DNA methylation
Oncol Prog,2021,19(7)
喉癌是耳鼻咽喉-头颈外科第二高发肿瘤，其中喉鳞状细胞癌占全部病理分型的95%[1-2]。

喉鳞状细胞癌的病因复杂，近年来其发病率不断增加，患者的晚期生存率仍较低[3]。

随着机器人辅助手术在头颈部肿瘤治疗中的开展及食管、气管一期重建技术的不断成熟，中国临床肿瘤学会（Chinese Society of Clinical Oncology，CSCO）头颈部肿瘤诊疗指南也强调手术治疗对提高生存率的重要性，但仍面临术后功能障碍、生活质量下降等严峻挑战[4]。

另一方面，头颈部肿瘤手术仍被认为是外科“皇冠”，手术水平直接关系到患者的生存预后。

因此，进一步阐明喉鳞状细胞癌的分子机制，寻找潜在的生物学靶点，对喉鳞状细胞癌进行早期诊断和治疗仍十分必要。

信使RNA（messenger RNA，mRNA）、微小RNA（microRNA，miRNA）、长链非编码RNA（long noncoding RNA，lncRNA）在喉鳞状细胞癌的发生、发展中发挥重要作用。

本研究通过肿瘤基因图谱（The Cancer Genome Atlas，TCGA）转录组数据挖掘，分别阐明lncRNA、miR-NA和mRNA在喉鳞状细胞癌中的作用机制及其潜在联系，现报道如下。

1资料与方法
1.1数据获取及差异分析
在TCGA中下载喉鳞状细胞癌转录组数据和相关临床信息，其中肿瘤样本111例，正常样本12例。

采用“edgR”包进行差异分析，设|logFC|﹥2，错误发现率（false discovery rate，FDR）﹤0.01，分别得到差异lncRNA、差异mRNA和差异miRNA，采用R软件绘制可视化图形并进行统计分析。

1.2蛋白相互作用网络分析
采用STRING11.0（https:///cgi/in-put.pl）在线工具将喉鳞状细胞癌Top500差异mRNA构建蛋白相互作用（protein-protein interac-tion，PPI）网络，相互作用值=0.7；Cytoscape插件CentiScape2.2、cytoHubba0.1用于寻找关键基因，并进行可视化展示。

1.3功能富集和通路富集分析进行细胞组分（cellular component，CC）、生物学过程（biological process，BP）、分子功能（molecular function，MF）的基因文本（gene ontology，GO）富集及京都基因和基因组百科全书（Kyoto Encyclope-dia of Genes and Genomes，KEGG）信号通路富集。

1.4MYHAS的共表达mRNA（cor-mRNA）
采用Pearson方法计算cor-mRNA，并取相关系数cor≥0.4。

1.5构建竞争性内源RNA（endoge
competitive endoge--nous RNA，ceRNA）网络和转录因子网络
ceRNA网络：采用本地版miRanda和Tar-getscan分别预测lncRNA结合的miRNA，并与差异miRNA取交集。

再通过miRTarBase、Targetscan、miRBD分别预测miRNA结合的mRNA，并与差异mRNA取交集，建立lncRNA-miRNA-mRNA网络联系，并通过Cytoscape3.7.2软件进行可视化展示。

转录因子网络：将钙信号通路相关的mRNA 于DA VID6.8在线网站（https:///）中进行转录因子富集，建立lncRNA相关mRNA的转录因子联系，并通过Cytoscape3.7.2软件进行可视化展示。

2结果
2.1一般资料
肿瘤患者的年龄、肿瘤T分期、病理分级、饮酒、吸烟情况对肿瘤患者总体生存率的影响不明显。

远处转移对患者总生存率的影响显著，临床转移患者的5年生存率为0%，明显低于无临床转移患者的50.33%，差异有统计学意义（χ2=12.230，P=0.0005）；病理转移患者的5年生存率为0%，明显低于无病理转移患者的44.83%，差异有统计学意义（χ2=18.260，P=0.0001）。

女性患者的5年生存率为30.00%，低于男性患者的50.11%，差异有统计学意义（χ2=4.454，P=0.0348）。

发生病理淋巴结转移患者的5年生存率为48.08%，低于未发生病理淋巴结转移患者的70.14%，差异有统计学意义（χ2= 4.539，P=0.0331），但临床上发生淋巴结转移与未发生淋巴结转移患者的5年生存率比较，差异无统2
672
ONCOLOGY PROGRESS,Apr 2021V ol.19No.7
2.2lncRNA 、miRNA 、mRNA 差异表达情况
lncRNA 、miRNA 和mRNA 各自的主成分分析结果显示，肿瘤组与对照组间表达谱数据相互独立，组内重复性较好，数据可信度高（图1A 、1C 和1E ）。

当|logFC|﹥2，FDR ﹤0.01时，得到的差异ln-cRNA 共687个，其中上调561个，下调126个（图1B ）；得到的差异miRNA 共54个，其中上调28个，下调26个（图1D ）；得到的差异mRNA 共1818个，其中上调950个，下调868个（图1F ）。

2.3差异lncRNA 和miRNA 表达的喉鳞状细胞癌预后的单因素分析
根据RNA-Seq 临床数据计算了差异lncRNA 表达的喉鳞状细胞癌预后的单因素分析结果，取P ﹤0.05，共有64个lncRNA 满足条件，根据是否能在NCBI 中查到相关lncRNA 信息，总生存率是否有统计学意义，筛选了11个lncRNA ，其中LINC00278、LINC00689、MYHAS 属于低风险基因（HR ﹤1），MNX1-AS1、LINC02575、HOXB-AS4、LSAMP-AS1、LINC02086、IGFL2-AS1、LINC02253、CASC20属于高风险基因（HR ﹥1）。

同理，计算出了差异miRNA 表达的喉鳞状细胞癌预后的单因素383、miRNA-196a-2、miRNA-196a-1、miRNA-100、
miRNA-4652是高风险基因（HR ﹥1），miRNA-99a 、miRNA-301a 是低风险基因（HR ﹤1）。

（表1）
2.4差异mRNA 功能分析和关键基因筛选
将1818个差异mRNA 用Clusterprofiler 进行GO 分析和KEGG 信号通路分析，取P 值或P .adjust 值﹤0.05。

分子功能（MF ）以通道活性和激活跨膜转运活性为主，细胞组分（CC ）以细胞外基质（ex-tracelluar matrix ，ECM ）包含的胶原纤维和转运复合体为主，生物学过程（BP ）以细胞外基质形成和调节离子跨膜转运为主（图2A~2C ）。

KEGG 信号通路中以细胞因子-细胞因子受体间相互作用和钙信号通路为主要作用通路，其中也包括环磷酸腺苷（cyclic adenosine monophosphate ，cAMP ）信号通路、环磷酸鸟苷（cyclic guanosine monophosphate ，cGMP ）-cGMP 依赖性蛋白激酶（cGMP-dependent protein kinase ，PKG ）信号通路和ECM 受体相互作用等肿瘤发生、发展的经典信号通路（图2D ）。

将差异倍数前500的mRNA 进行String 在线分析，构建PPI 网络，并通过Centiscape 插件进行关键基因筛选，根据Degree 进行筛选，绿色﹤10，青色≥10，紫色≥15，红色≥20，其中ACTN2、ACTN3、MYL1、MYL2、MYL3、MYH6、TTN 网络连接最多，属于关键基因（图2E ）。

2.5MYHAS 的cor-mRNA 、关键基因筛选及功能富集
MYHAS 的共表达mRNA （cor ≥0.4）共253个，采用插件cytoHubba 分别计算Closeness 、Degree 、EPC 、MCC 、MNC 和Radiality 6种算法下前10的基因并取交集，得到关键基因：TNN 、MYL1、MYL3、
B
A 图1lncRNA 、miRNA 、mRNA 的主成分分析的主成分分析、、火山图
注：A 、C 和E 分别代表lncRNA 、miRNA 和mRNA 的主成分分析，红点代表肿瘤样本，蓝点代表对照样本；B 、D 和F 分别代表lncRNA 、miRNA 和mRNA 的火山图，横坐标表示FDR 值，纵坐标表示差异倍数，红色表示上调的基因，绿色表示下调的基因
D
C F
E 673
《癌症进展》2021年4月第19卷第7
期
C
B
A
D
E
注：A~C.GO 分析中分子功能（MF ）、细胞组分（CC ）和生物学过程（BP ）的富集情况，横坐标表示富集基因占所有基因的比例，纵坐标表示不同的富集条目名称，圆圈大小表示富集基因的多少，圆圈越大，富集的基因越多，圆圈越小，富集的基因越少，颜色由绿到红表示P .adjust 值由大到小；D.富集的KEGG 信号通路，横坐标表示富集基因占所有基因的比例，纵坐标表示不同的富集条目名称，圆圈大小表示富集基因的多少，圆圈越大，富集的基因越多，圆圈越小，富集的基因越少，颜色由绿到红表示P 值由大到小；E.PPI 网络，连线越多说明该基因在喉鳞状细胞癌发生过程中越关键，绿色<10，青色≥10，紫色≥15，红色≥20
674
ONCOLOGY PROGRESS,Apr2021V ol.19No.7
调节金属离子转运、跨膜转运复合体和肌动蛋白结合等（图3A~3C）。

信号通路分析显示，cor-mRNA主要富集在钙信号通路、cGMP-PKG信号通路和催产素信号通路等（图3D）。

其中钙离子信号通路中富集的基因分别是PRKACA、CAMK2A、SLC8A3、CAMK2B、RYR1、PLN、CASQ2、PHKG1、HRC、MYLK3、TNNC1、CACNA1S、MYLK2、CASQ1、ATP2A1、TNNC2。

2.6ceRNA网络和转录因子网络
将单因素分析与生存分析中P﹤0.05的ln-cRNA取交集，得到21个lncRNA，并分别进行mi-Randa和Targetscan靶基因预测并取交集，得到2363个与lncRNA可以结合的miRNA，然后将预测的miRNA与差异miRNA的成熟体取交集，得到66个miRNA。

对66个miRNA分别进行miRTarBase、Targetscan和MiRDB靶基因预测后与差异mRNA 取交集，并构建ceRNA网络（图4A）。

ceRNA网络中有11个lncRNA（青色菱形），分别为LINC02576、LINC02086、AC020659.1、LINC00528、LINC00689、HOXB-AS4、MNX1-AS1、LINC00278、AC010624.1、AC016773.1、MYHAS；5个miRNA （绿色三角形），分别为has-miRNA-206、has-miR-NA-573、has-miRNA-3662、has-miRNA-133b、has-miRNA-449a；12个mRNA（其中7个红色6边型为上调mRNA，5个蓝色圆形为下调mRNA），分别为STC2、P AX3、NETO2、EIF5A2、ADPRHL1、SYNM、ACTA1、GPR156、CLIC5、BMP3、HNF4A、FOXL2。

将cor-mRNA钙信号通路中的16个基因进行转录因子富集，并绘制转录因子调控网络，网络中绿色椭圆形代表mRNA，蓝色菱形代表转录因子，这些转录因子分别是MYOD、RSRFC4、TAL1BETAITF2、YY1、P300、PAX3、PAX5、ZID和OLF1（图4B）。

2.7DNA甲基化分析
为了进一步阐明关键基因在喉鳞状细胞癌中的调控作用，进一步分析了TCGA和GEO中的喉鳞状细胞癌DNA甲基化数据。

设Δβ绝对值≥0.4，P﹤0.05，取二者交集，共得到75个甲基化位点；其中cg15700197和cg17892556分别为甲基化差异最大的位点，cg15700197在正常组织中高甲基化，在喉鳞状细胞癌组织中低甲基化，cg17892556在正常组织中低甲基化，在喉鳞状细胞癌组织中高甲基化。

这75个差异甲基化位点主要位于1stExon、TSS1500和TSS200，其次为5'非翻译区（untranslat-ed region，UTR）和基因体区，3'非翻译区和基因间隔区（intergenic region，IGR）分布较少。

上述位点对应高甲基化基因56个，低甲基化基因15个，ZNF625是Δβ值最大的高甲基化基因，OR10J1是Δβ
值最大的低甲基化基因。

B
A
D
C
注：A~C.GO富集结果，D.KEGG信号通路富集结果，横坐标表示富集基因占所有基因的比例，纵坐标表示不同的富集条目名称，圆圈大小表示富集基因的多少，圆圈越大，富集的基因越多，圆圈越小，富集的基因越少，颜色由绿到红表示P值或P.adjust值由大到小
675
《癌症进展》2021年4月第19卷第7期
3讨论
喉鳞状细胞癌是头颈部常见的肿瘤之一，病因较多，影响因素复杂。

本研究针对喉鳞状细胞癌患者的一般情况进行分析，结果显示，患者年龄、T 分期、病理分级、饮酒、吸烟情况对总生存率的影响不明显，而远处转移、淋巴结转移会降低患者的总生存率。

关于肿瘤分期（T 分期）和病理分级对肿瘤预后无影响，可能与目前分期系统的局限性有关；因此针对喉鳞状细胞癌的分期，根据肿瘤萌芽活性（tumor budding activity ，BA ）和细胞巢大小（cell nest size ，CNS ）进行分期对预后似乎更具指导意义[5]。

本研究结果中淋巴结浸润对喉鳞状细胞癌预后的影响与Ho 等[6]研究结果一致，尤其是病理组织上的淋巴结转移对指导临床治疗意义重大，当然这也可能与从临床到病理诊断的过程增加了“开始治疗时间（time to treatment initiation ，TTI ）”，从而降低了预后有关[7]。

近年来有关lncRNA 在喉鳞状细胞癌中作用的研究较多，尤其是miRNA “海绵”靶基因是喉癌发生、发展的重要机制之一[8]。

作为ceRNA ，lncRNA 通过靶向结合miRNA 间接调控靶基因，是较好的生物标志物和潜在的治疗靶点[9]。

生存分析显示，miRNA-210、miRNA-455和miRNA-4326的表达量与喉鳞状细胞癌预后相关（数据未展示），其中有关miRNA-210在肺癌、胶质瘤、口腔癌、骨肉瘤等恶性肿瘤中的作用已广泛报道，并且miRNA-210在组织、血清和外泌体中均可存在[10]。

miRNA-455也可通过自噬、p21活化蛋白激酶2（p21-activated kinase 2，PAK2）轴和Janus 激酶1（Janus kinase 1，JAK1）轴等作用方式抑制肿瘤生长，并且miRNA-455上述功能的发挥有赖于不同lncRNA 的作用，这
些lncRNA 作为ceRNA 竞争性结合miRNA-455[11-13]
；当然，也有较多miRNA-455单独发挥作用的报道[14]。

[15]基因表达促进肺癌细胞增殖。

但遗憾的是，目前
尚未见上述3个miRNA 在喉鳞状细胞癌中被报道，有待于进一步验证。

在具有抑癌作用的lncRNA 中，LINC00689、LINC00278已分别在胶质瘤和食管鳞状细胞癌中报道，尤其是LINC00278编码的多肽YY1BM ，可抑制YY1与雄激素受体结合从而降低真核生物延伸因子2激酶（eukaryotic elongation factor 2kinase ，eEF2K ）的表达，促进肿瘤细胞凋亡[16-17]。

喉结是第二性征标志，雄激素受体高表达，由此可推断LINC00278在喉鳞状细胞癌中是一个较好的潜在靶点[18]。

在具有促癌作用的lncRNA 中，IGFL2-AS1、LSAMP-AS1和MNX1-AS1已分别在胃癌、前列腺癌和食管鳞状细胞癌中报道，为下一步喉鳞状细胞癌的研究提供了较好的参考意义[19-21]。

基于目前报道的有关lncRNA-miRNA-mRNA 的相互作用关系，本研究利用差异lncRNA 、差异miRNA 和差异mRNA 构建了ceRNA 网络，通过该网络发现，网络中miRNA 既无生存意义，又无单因素意义；网络中mRNA 亦不属于Centiscape 中的关键基因；同时，与miRTarBase 疾病搜索中喉鳞状细胞癌的miRNA 无交集。

因此认为lncRNA 并不能完全主导miRNA 和mRNA 发挥作用，lncRNA 与miRNA 在喉鳞状细胞癌发病过程中各自存在独立的作用机制。

为了进一步说明lncRNA 的功能，本研究选择了MYHAS 作为代表，通过分析MYHAS 共表达基因cor-mRNA ，发现cor-mRNA 得到的5个关键基因全部与差异mRNA 得到的关键基因吻合；并且cor-mRNA 的通路分析中主要信号通路钙信号通路、cGMP-PKG 信号通路和cAMP 信号通路亦是差异mRNA 的主要信号通路。

有关cAMP 、cGMP-PKG 和钙信号通路在不同肿瘤中的调控作用已有很多报道[22-24]，有研究报道了cGMP-PKG 在
B
A 图4ceRNA 网络和转录因子网络
注：A.lncRNA-miRNA-mRNA 构成的ceRNA 网络，青色菱形代表lncRNA ，绿色三角形代表miRNA ，红色六边型为上调mRNA ，蓝色圆形为下调mRNA ；B.钙离子信号通路中cor-mRNA 富集的转录因子网络，蓝色菱形是转录因子，绿色椭圆形是cor-mRNA
676
ONCOLOGY PROGRESS,Apr2021V ol.19No.7
头颈部鳞状细胞癌中是较好的治疗靶点，并且全球第一个鸟苷酸环化酶激动剂利那洛肽已于2019年初在中国上市[25-26]。

此外，Park和Juhnn[27]亦报道cAMP信号激活可通过增加去乙酰化酶8的表达促进顺铂对肺癌细胞的凋亡作用。

而关于钙信号通路在肿瘤中的作用，Raynal等[28]研究表明，钙信号通路中一系列激酶相互作用后，引起钙/钙调素依赖蛋白激酶Ⅱ（calcium/calmodulin-dependent pro-tein kinaseⅡ，CAMKⅡ）激活，使甲基化结合域蛋白甲基化CPG结合蛋白2（methyl-CPG binding pro-tein2，MeCP2）出核，导致基因组发生去甲基化，使已经甲基化的抑癌基因重新激活，以发挥抑癌作用。

MYHAS属于抑癌基因，其通过相关基因cor-mRNA在喉鳞状细胞癌组织中参与钙信号通路的失活，从而引起下游靶基因高甲基化。

基于上述结果，本研究又分析了数据库中喉鳞状细胞癌的DNA甲基化数据，共筛选出了71个基因；本实验室已经对ZNF667基因进行了试验验证，其在喉鳞状细胞癌组织和细胞中均呈高甲基化、低表达，发挥抑制喉鳞状细胞癌细胞增殖、侵袭的作用[29]。

由此认为MYHAS可以较好地代表喉鳞状细胞癌中的差异lncRNA，在喉鳞状细胞癌中是一个较好的生物标志物和治疗靶点，具有潜在研究价值。

同时，本研究针对cor-mRNA中参与钙信号通路的16个基因构建了转录因子调控网络，拟进一步说明非编码RNA在喉鳞状细胞癌信号通路中的上下游调控关系。

这些转录因子与MYHAS相互作用，直接在转录水平调控靶基因表达。

网络中PAX5是B淋巴细胞特异性转录因子，在B细胞淋巴瘤的发生、发展中具有重要作用，但也见于在其他恶性肿瘤中的报道[30]，并且Brazao等[31]报道了PAX5对不同阶段B细胞中lncRNA的调控作用。

而转录因子PAX3的功能较为复杂，这是因为PAX3具有7个可变剪切体，不同剪切体在肿瘤中表现的功能不同[32]。

P300被认为是转录共激活因子，是一种组蛋白乙酰转移酶，通常被招募到转录增强子中，并通过乙酰化染色质来调节基因表达。

Liang等[33]研究报道，P300可在食管鳞状细胞癌中调控LINC00460的表达，这亦为喉鳞状细胞癌的研究提供了可行性。

转录因子YY1是研究较多的肿瘤相关转录因子之一，其可调控lncRNA的表达，引起直肠癌、恶性胶质瘤、乳腺癌和肝癌等肿瘤的发生、发展，其中有lncRNA NPCCAT1通过上调YY1促进鼻咽癌发生发展的报道，这对研究YY1在喉鳞状细胞癌中的作用具有重要参考意义[34]。

MYOD属于肌源性转录因子，在横纹肌肉瘤中与同类转录因子MYF5共同发挥重要作用[35]；此外Dodd等[36]亦报道了MYOD通过直接结合miRNA-182的启动子增加miRNA-182的表达，从而发挥促肿瘤作用。

因此，在喉鳞状细胞癌中转录因子MYOD与lncRNA的作用还是值得期待。

综上所述，公共数据分析显示喉鳞状细胞癌中存在众多差异lncRNA和miRNA，其中lncRNA MYHAS通过抑制钙信号通路引起下游靶基因甲基化可能是喉鳞状细胞癌发生、发展的主要机制之一，lncRNA MYHAS是一个较好的生物标志物和治疗靶点，具有潜在研究价值。

参考文献
[1]Bray F,Ferlay J,Soerjomataram I,et al.Global cancer sta-
tistics2018:GLOBOCAN estimates of incidence and mor-tality worldwide for36cancers in185countries[J].CA Cancer J Clin,2018,68(6):394-424.
[2]Gale N,Poljak M,Zidar N.Update from the4th edition of
the World Health Organization classification of head and neck tumours:what is new in the2017WHO blue book for tumours of the hypopharynx,larynx,trachea and parapha-ryngeal space[J].Head Neck Pathol,2017,11(1):23-32. [3]Steuer CE,El-Deiry M,Parks JR,et al.An update on lar-
ynx cancer[J].CA Cancer J Clin,2017,67(1):31-50.
[4]Chinese Society of Clinical Oncology(CSCO)diagnosis
and treatment guidelines for head and neck cancer2018 (English version)[J].Chin J Cancer Res,2019,31(1):84-
98.
[5]Boxberg M,Kuhn PH,Reiser M,et al.Tumor budding and
cell nest size are highly prognostic in laryngeal and hypo-pharyngeal squamous cell carcinoma:further evidence for a
Pathol,2019,43(3):303-313.
[6]Ho AS,Kim S,Tighiouart M,et al.Association of quantita-
tive metastatic lymph node burden with survival in hypo-pharyngeal and laryngeal cancer[J].JAMA Oncol,2018,4
(7):985-989.
[7]Xiao R,Ward MC,Yang K,et al.Increased pathologic up-
staging with rising time to treatment initiation for head and neck cancer:a mechanism for increased mortality[J].Can-cer,2018,124(7):1400-1414.
[8]Cossu AM,Mosca L,Zappavigna S,et al.Long non-coding
RNAs as important biomarkers in laryngeal cancer and oth-er head and neck tumours[J].Int J Mol Sci,2019,20(14): 3444.
[9]Anastasiadou E,Jacob LS,Slack FJ.Non-coding RNA net-
works in cancer[J].Nat Rev Cancer,2018,18(1):5-18. [10]Feng S,He A,Wang D,et al.Diagnostic significance of
miR-210as a potential tumor biomarker of human cancer detection:an updated pooled analysis of30articles[J].On-
677
《癌症进展》2021年4月第19卷第7期
profiles analysis revealed lncRNA OXCT1-AS1promoted bladder cancer cell aggressiveness via miR-455-5p/JAK1 signaling[J].J Cell Physiol,2019,234(8):13592-13601. [12]Ni X,Ding Y,Yuan H,et al.Long non-coding RNA ZEB1-
AS1promotes colon adenocarcinoma malignant progres-sion via miR-455-3p/PAK2axis[J].Cell Prolif,2020,53
(1):e12723.
[13]Zhang H,Luo C,Zhang G.LncRNA MCM3AP-AS1regu-
lates epidermal growth factor receptor and autophagy to promote hepatocellular carcinoma metastasis by interact-ing with miR-455[J].DNA Cell Biol,2019,38(8):857-864.
[14]Zhan T,Zhu Q,Han Z,et al.miR-455-3p functions as a tu-
mor suppressor by restraining Wnt/beta-Catenin signaling via TAZ in pancreatic cancer[J].Cancer Manag Res,2020, 12:1483-1492.
[15]Xu G,Zhang Z,Zhang L,et al.miR-4326promotes lung
cancer cell proliferation through targeting tumor suppres-sor APC2[J].Mol Cell Biochem,2018,443(1-2):151-157.
[16]Liu X,Zhu Q,Guo Y,et al.LncRNA LINC00689pro-
motes the growth,metastasis and glycolysis of glioma cells by targeting miR-338-3p/PKM2axis[J].Biomed Pharmacother,2019,117:109069.
[17]Wu S,Zhang L,Deng J,et al.A novel micropeptide encod-
ed by Y-linked LINC00278links cigarette smoking and AR signaling in male esophageal squamous cell carcinoma [J].Cancer Res,2020,80(13):2790-2803.
[18]Fei M,Zhang J,Zhou J,et al.Sex-related hormone recep-
tor in laryngeal squamous cell carcinoma:correlation with androgen estrogen-a and prolactin receptor expression and influence of prognosis[J].Acta Otolaryngol,2018,138(1): 66-72.
[19]Ma Y,Liu Y,Pu YS,et al.LncRNA IGFL2-AS1functions
as a ceRNA in regulating ARPP19through competitive binding to miR-802in gastric cancer[J].Mol Carcinog, 2020,59(3):311-322.
[20]Hua X,Liu Z,Zhou M,et al.LSAMP-AS1binds to mi-
croRNA-183-5p to suppress the progression of prostate cancer by up-regulating the tumor suppressor DCN[J].
EBioMedicine,2019,50:178-190.
[21]Chu J,Li H,Xing Y,et al.LncRNA MNX1-AS1promotes
progression of esophageal squamous cell carcinoma by regulating miR-34a/SIRT1axis[J].Biomed Pharmacother, 2019,116:109029.
[22]Kumar N,Prasad P,Jash E,et al.Insights into exchange
factor directly activated by cAMP(EPAC)as potential tar-get for cancer treatment[J].Mol Cell Biochem,2018,447
(1-2):77-92.
[23]Peak TC,Richman A,Gur S,et al.The role of PDE5inhib-
itors and the NO/cGMP pathway in cancer[J].Sex Med Rev,2016,4(1):74-84.
[24]Cui C,Merritt R,Fu L,et al.Targeting calcium signaling
in cancer therapy[J].Acta Pharm Sin B,2017,7(1):3-17.
[25]Tuttle TR,Mierzwa ML,Wells SI,et al.The cyclic GMP/
protein kinase G pathway as a therapeutic target in head and neck squamous cell carcinoma[J].Cancer Lett,2016, 370(2):279-285.
[26]程能能.新型、口服肠上皮细胞鸟苷酸环化酶C激动
剂——利那洛肽[J].中国临床药学杂志,2019,28(4): 291-295.
[27]Park JY,Juhnn YS.cAMP signaling increases histone
deacetylase8expression via the Epac2-Rap1A-Akt path-way in H1299lung cancer cells[J].Exp Mol Med,2017, 49(2):e297.
[28]Raynal NJ,Lee JT,Wang Y,et al.Targeting calcium sig-
naling induces epigenetic reactivation of tumor suppressor genes in cancer[J].Cancer Res,2016,76(6):1494-1505. [29]Meng W,Cui W,Zhao L,et al.Aberrant methylation and
downregulation of ZNF667-AS1and ZNF667promote the malignant progression of laryngeal squamous cell carcino-ma[J].J Biomed Sci,2019,26(1):13.
[30]Dong BW,Zhang WB,Qi SM,et al.Transactivation of
PTGS2by PAX5signaling potentiates cisplatin resistance in muscle-invasive bladder cancer cells[J].Biochem Bio-phys Res Commun,2018,503(4):2293-2300.
[31]Brazao TF,Johnson JS,Müller J,et al.Long noncoding
RNAs in B-cell development and activation[J].Blood, 2016,128(7):e10-e19.
[32]Wang Q,Kumar S,Slevin M,et al.Functional analysis of
alternative isoforms of the transcription factor PAX3in melanocytes in vitro[J].Cancer Res,2006,66(17):8574-8580.
[33]Liang Y,Wu Y,Chen X,et al.A novel long noncoding
RNA linc00460up-regulated by CBP/P300promotes car-cinogenesis in esophageal squamous cell carcinoma[J].
Biosci Rep,2017,37(5):BSR20171019.
[34]Su H,Liu L,Zhang Y,et al.Long noncoding RNA NPC-
CAT1promotes nasopharyngeal carcinoma progression via upregulating YY1[J].Biochimie,2019,157:184-194.
[35]Tenente IM,Hayes MN,Ignatius MS,et al.Myogenic reg-
ulatory transcription factors regulate growth in rhabdo-myosarcoma[J].Elife,2017,6:e19214.
[36]Dodd RD,Sachdeva M,Mito JK,et al.Myogenic tran-
scription factors regulate pro-metastatic miR-182[J].On-cogene,2016,35(14):1868-1875.
（收稿日期：2020-08-31）
678。