泪腺良性淋巴上皮病变基因网络分析及关键枢纽基因鉴定

张黎黎，刘文辉，刘小静，丁 岩，张 琪，刘 瑜，刘林嶓

(1.郑州大学第一附属医院整形外科，河南 郑州 450052；2.郑州大学第一附属医院眼科，河南 郑州 450052)

泪腺良性淋巴上皮病变(benign lymphoepithelial lesion，BLEL)是一种常表现为泪腺无痛性肿大的良性疾病，治疗方法以激素为主，还包括放疗、手术等[1]。但激素作用的个体差异明显，有些患者接受激素治疗的疗效并不理想。因此，进一步明确BLEL的发病机制是十分必要的。本研究基于BLEL基因芯片数据，构建了基因共表达网络，并鉴定其中的关键致枢纽病基因，有助于在分子水平上系统全面明确BLEL的发病机制，对BLEL的治疗也有一定的启示作用。

1 资料与方法

1.1 数据获取本研究选取了Gene Expression Omnibus(www.ncbi.nlm.nih.gov/geo/)数据库中编号为GSE76497的数据集，该数据集中包括18例BLEL和作为对照的18例眼眶海绵状血管瘤。该数据集也曾被Adzavon等[2]和Li等[3]用于研究BLEL的发病机制，但其研究方法与本研究并不相同。

1.2 差异基因及聚类分析为避免与BLEL发病关联较弱基因的干扰，本研究利用双尾Studentt检验求得了2组的差异表达基因，取检验水准α=0.05。在获得差异基因后，使用Mev软件[4]对2组进行聚类分析，排除异常样本的数据。

1.3 共表达网络构建基因共表达网络是基于加权基因共表达网络分析的方法构建的。WGCNA的方法在Web of Science数据库中被引已经超过5 000次[5]，其有效性得到了广泛的认可。不同于常规的分析方法，加权基因共表达网络分析将基因表达信息拟合为无尺度网络，再根据网络的拓扑学特性，将基因网络分为若干基因模块，并从中选择与表型相关性最高的关键枢纽基因模块，进一步明确模块中的关键枢纽基因[6]。……