广东汉族鼻咽癌患者HLA等位基因易感关联分析及中医易感体质初探

【摘要】鼻咽癌(nasopharyngealcarcinoma,NPC)为我国南方常见恶性肿瘤,尤其是在广东地区,据统计广东地区的世界年龄标化发病率可达20/10万以上。鼻咽癌一般位置较深,周围邻近大量重要器官及组织,手术难以彻底切除干净,故手术并非鼻咽癌治疗的首选。鼻咽癌一般对放射线极为敏感,其解剖学及生物学特点决定了其主要治疗方式为放射治疗。早期鼻咽癌单纯放射治疗5年生存率可达90%以上。然而,放疗本身也难以避免带来一系列不良反应。更遗憾的是,由于鼻咽癌早期多无特异性的症状,使得大多数患者就诊时即已为晚期。局部晚期鼻咽癌患者单纯放疗效果差,多需联合化疗,目前一般主张同步放化疗。放疗与化疗的联合应用,毒副反应也叠加,部分患者难以耐受,使得不能完成全程的同步放化疗,故生存率的提高更无从谈起。可见,鼻咽癌是严重影响我国南方尤其是广东地区人民身体健康的疾病。鼻咽癌确切发病原因尚不明确。目前普遍认为与理化环境因素、EBV感染和HLA基因等遗传因素有关。HLA(人类白细胞抗原)基因具有高度多态性,其编码的HLA分子在免疫识别“非我”抗原、呈递内外源性抗原,对免疫应答中起着决定性的作用,同时HLA基因具有的地区和种族差异也可以作为一种遗传标记,用于有遗传倾向性疾病的关联研究,连锁分析提示HLA具有NPC遗传易感性。文献报道,中国南方人群中HLA-A*02.B*58和DRB1*03可能是NPC易感基因,B*40可能是NPC保护基因,但是关于广东汉族鼻咽癌患者与]HLLA-A、B、DRB1高分辨等位基因易感关联分析却鲜有报道。中医体质学是以中医理论为指导,研究人体体质特征、体质类型的生理病理特点,分析疾病反应状态、病变性质及发展趋势的一门学科。而HLA系统的特性与中医学“体质”学说之间存在某些内在联系,反映了体质的先天遗传特征和稳定性。HLA系统为数众多的等位基因具有的多态性与中医体质类型的多态性相似。对广东汉族鼻咽癌患者进行体质分析及关联研究,将有助于解释不同的中医体质对NPC发生发展的影响,进一步认识中医体质与疾病的相关性。基于此,本课题通过研究广东汉族鼻咽癌患者和广东汉族健康志愿者进行的HLA等位基因,探讨广东汉族鼻咽癌患者的HLA易感关联。同时对上述两组研究对象中年龄介乎18-65周岁且同意进行中医体质调查的NPC患者和健康志愿者进行中医体质调查,探讨广东汉族鼻咽癌患者的中医易感体质,试图为鼻咽癌的早防早治提供中西医两方面的思考。目的：为探寻广东汉族鼻咽癌患者的HLA基因易感性及中医易感体质,试图为鼻咽癌的早防早治提供一些思路和参考。方法：2012年1月-2012年10月期间在广东省内两所医院(均为大学附属医院)随机选取广东汉族鼻咽癌患者114例为NPC患者组和121例广东汉族健康志愿者为对照组,运用SBT基因测序技术对该114例广东汉族鼻咽癌患者和该121例广东汉族健康志愿者进行HLA-A、B、DRB1高分辨等位基因分析,探讨广东汉族鼻咽癌的HLA基因易感关联。在上述114例广东汉族NPC患者中随机选取年龄介乎18-65周岁且同意进行中医体质调查的NPC患者(共收集到100例NPC患者)；在上述121例健康志愿者体检中随机选取年龄介乎18-65周岁且同意进行中医体质调查的健康志愿者(共收到100例健康志愿者)。对该100例NPC患者和该100例健康志愿者运用中华中医药学会2009年4月颁布的《中医体质分类和评定标准》中的中医体质量表进行中医体质问卷调查。在中医体质调查前,告知NPC组被调查对象根据本次患病一个月前,一年以内的实际情况完成问卷；告知健康志愿者组被调查对象,根据近一年以内的实际情况完成问卷。广东汉族鼻咽癌患者HLA等位基因易感关联分析中采用以下统计学方法：应用Arlequin软件3.11对健康对照组进行Hardy-Weinberg平衡吻合度检验(P>0.05,说明所该组达到遗传平衡)。分别计算NPC组和健康对照组中HLA等位基因的表现型频率以及基因频率,并比较两组间HLA等位基因频率分布的差异,采用Pearson卡方检验对患者组与对照组的HLA-A、B、DRB1抗原频率逐一进行关联分析(纯合子计作一个),必要时采用Yates检验方法与Fisher确切概率法计算确切P值,同时按Woolf公式计算相OR值及95%的可信区间(CI)。计量资料用均数±标准差表示,计数资料用率表示,均数比较用两样本t检验,率的比较用卡方检验,均为双侧检验,检验水准a=0.05。采用Excel表格和SPSS13.0统计软件进行统计分析。广东汉族鼻咽癌患者中医易感体质初探中采用以下统计学方法：两组率的比较用卡方检验,双侧检验,检验水准a=0.05,其中气虚质(qideficiencyconstitution)、气郁质(QidepressionconstitutionO的比较使用了Bonferroni校正检验水准a=0.025。结果：对NPC组进行HLA-A、B、DRB1基因多态性分析结果表明,114例NPC患者共表达75个HLA等位基因,占目前发现的总HLA等位基因的1.0267%,其中HLA-A、B、DRB1分别为16、34和25个,分别占目前发现的总HLA等位基因的0.620%,1.035%和1.388%(IMGT/HLAdatabase2014年01月17日公布Release3.15.0,HLA-A、B、DRB1等位基因数为7305个,其中HLA-A位点2579个,B位点3285个,DRB1位点1441个)。NPC患者组与健康对照组共表达21个HLA-A等位基因,其中NPC组共表达16个,占]MGT/HLAdatabase已公布的0.6226%,对照组表达20个,占HLA-A等位基因0.7782%,其中HLA-A*24:10仅在患者组表达,A*01:02、A*11:04、A*24:07、A*31:01和A*66：01仅在对照组中表达。NPC组中,基因频率>0.1的等位基因为A*02:07(0.194)、A*11:01(0.257)、A*24:02(0.157)和A*33：03(0.107),健康对照组基因频率>0.1的等位基因为A*02：07(0.142)、A*11：01(0.357)和A*24：02(0.167)。分别对两组间HLA-A等位基因频率进行比较,其中NPC组中A*11：01的基因频率与正常对照组相比存在显著的统计学差异(A*11:01,0.257vs0.357,OR=0.570,95%CI=0.340-0.956,X2=4.57,P=0.033)。两组间其余HLA-A等位基因频率的分布均无显著统计学差异。NPC患者组与健康对照组共表达44个HLA-B等位基因,其中NPC组共表达34个,占IMGT/HLAdatabase已公布HLA-B等位基因的0.6226%,对照组表达36个,占全部HLA-B等位基因的1.035%,其中B*15:05、B*15:17、B*15:25、B*27:04、B*27:05、B*27:06、B*37:01、B*52:01、B*54:02、B*54:02和B*55：02仅在患者组表达,B*08：01、B*14：02、B*15:18、B*44:02、B*50:01、B*55:01、B*56:01和B*57：01仅在对照组中表达。NPC患者组频率>0.1的基因为B*40：01(0.111)、B*46：01(0.233)和B*58：01(0.121),健康对照组频率>0.1的基因为B*40：01(0.157)、B*46：01(0.162)和B*58：01(0.0662)；分别对两组间HLA-B等位基因频率进行比较,其中NPC组中B*58：01基因频率与正常对照组存在显著的统计学差异(0.121vs0.064,OR=2.0879,95%CI=1.041-4.186,X2=4.17,P=0.036)。两组间其余HLA-B等位基因频率的分布均无显著统计学差异。NPC患者组与健康对照组共表达个27个HLA-DRB1等位基因,其中NPC组共表达25个,占MGT/HLAdatabase已公布HLA-DRB1等位基因的1.7349%,对照组表达25个,占HLA-DRB1等位基因的1.7349%,其中DRB1*04:07、DRB1*08:02、DRB1*08:09和DRB1*14:07仅在患者组表达,DRB1*01:02、DRB1*04:01、DRB1*11:03和DRB1*13:01仅在对照组中表达。NPC患者组基因频率>0.1的基因为DRB1*03:01(0.116)、DRB1*09:01(0.141)和DRB1*15:01(0.121),健康对照组基因频率>0.1的基因为DRB1*09:01(0.167)、DRB1*12:02(0.100)和DRB1*15:01(0.128),分别比较两组间HLA-DRB1*等位基因频率,其中NPC组中DRB1*03:01和DRB1*04:03基因频率与正常对照组存在显著的统计学差异(DRB1*03:01:0.116vs0.047,OR=2.088,95%CI=1.311-6.020,X2=7459,P=0.006；DRB1*04:03:0.022vs0.0064,OR=0.324,95%CI=0.114-0.923,χ2=4.57,P=0.028)。两组间其余HLA-DRB1等位基因频率的分布均无显著统计学差异。114例NPC患者中,EBV-EA-IgG阳性94例,占82.456%,阴性20例,占17.544%。对EBV-EA-IgG阳性及阴性的NPC患者的HLA-A、B、DRB1等位基因频率分布进行比较,其中HLA-A*33:03、B*38:01和DRB1*01:01在两组间分布差异存在统计学意义(A*33:03:OR=1.282,95%CI:1.149-1.429,X2=6.130,P=0.013；B*38:01,OR=6.222,95%CI:4.075-9.501,P=0.029；DRB1*01:01,OR=6.222,95%CI：4.075-9.501,P=0.029),其余HLA等位基因频率均无显著统计学差异。对114例NPC患者使用荧光定量PCR技术(real-timePCR进行EBV-DNA检测,其中NPC患者组92例为阳性,占79.8246%,22例为阴性,占20.1754%。对EBV-DNA阳性及阴性的NPC患者的HLA-A、B、DRB1等位基因频率分布进行比较,其中对HLA-B*15:19和B*46：01在两组间的分布差异存在统计学意义(B*15：19,OR=5.600,95%CI:3.764-8.331,P=0.036；B*46:01:OR=6.222,95％CI:1.240-12.575,χ2=5.976,P=0.015),其余HLA等位基因频率均无显著统计学差异。中医体质调查发现：广东汉族鼻咽癌患者此次患病一个月以前、近一年以内,其中医体质以偏颇体质为主,其中以气虚质最多见,其次为气郁质。健康志愿者组近一年以内的中医体质以平和质为主,共占全部健康志愿者的55%。结论：从等位基因层次来说,在广东汉族人群中HLA-B*58:01和DRB1*03:01可能是NPC易感关联的致病基因,而HLA-A*11:01和DRB1*04:03可能是NPC易感关联的保护基因。研究显示,NPC患者中携带HLA-A*33:03等位基因者EBV-EA-IgG抗体产生呈正相关,携带B*38：01和DRB1*01：01等位基因者可能EBV-EA-IgG抗体产生呈负相关,而NPC患者中携带HLA-B*46:01等位基因者与EBVDNA阳性呈正相关,携带B*15：19等位基因者可能与EBVDNA阳性呈负相关。因此,机体对在EBV的固有免疫的免疫应答过程中,EBV-EA-IgG抗体的出现与否,可能与HLA多态性有关,然而,携带保护性HLA等位基因的机体,可以更有效提呈EBV,其CTL(细胞毒性T淋巴细胞)对EBV感染细胞的杀伤能力更强,患NPC的风险也相应减低。对广东汉族鼻咽癌患者中医易感体质调查,结果提示气虚质、气郁质可能是广东汉族鼻咽癌患者的中医易感体质；平和质可能是广东汉族人群鼻咽癌的保护性体质
【作者】高卓维；
【导师】周迎春；
【作者基本信息】南方医科大学，中西医结合临床（专业学位），2014，博士
【关键词】鼻咽癌；人类白细胞抗原；中医体质；

【参考文献】
[1]孙振.一种面向时间不确定性问题的故障诊断方法研究[D].郑州大学,计算机软件与理论,2013,硕士.
[2]王坤.在线投票系统的设计与实现[D].华中科技大学,软件工程,2013,硕士.
[3]刘喜华.循环经济视角的酒店企业绿色管理模式研究[D].中南林业科技大学,生态学,2012,硕士.
[4]曾垂省.可控多孔生物陶瓷的制备及表征[D].武汉理工大学,生物医学工程,2004,硕士.
[5]龙泓君.HS公司薪酬体系再设计研究[D].大连理工大学,工商管理（专业学位）,2012,硕士.
[6]朱少英.管理诊断方案生成过程的主要逻辑环节分析研究[D].广西大学,政治经济学,2004,硕士.
[7]路达.长春市老年人体育健身现状的调查与分析[D].辽宁师范大学,体育教育训练学,2011,硕士.
[8]申屠胜男,徐龙稳,孟彬,阮健.双向比例电磁铁的设计及仿真分析[J].机械设计与研究,2015,01:27-30.
[9]王崇蕊.电解黑色磷化膜制备及性能研究[D].沈阳理工大学,应用化学,2012,硕士.
[10]艾恕.全球光通信市场最新动向[J].世界电信.2004(06)
[11]彭叶苹.基于遗传算法的测试数据自动生成方法研究[D].广东工业大学,计算机应用技术,2013,硕士.
[12]赵渊,黄旭,谢开贵.计及异步化同步发电机的电网可靠性评估模型[J].电力系统自动化,2011,19:6-11.
[13]谢坚.存量房批量评估系统的研究与实现[D].杭州电子科技大学,计算机技术,2013,硕士.
[14]宋雪.Hermite矩阵及四元数体上矩阵的一些不等式的问题研究[D].重庆大学,应用数学,2014,硕士.
[15]苏东涛,胡亚澜.大银鱼移植增殖技术[J].农民科技培训.2002(09)
[16]张锐兵.后危机时代我国金融消费者权益保护的法律思考[D].华中科技大学,法律,2013,硕士.
[17]孙会勤.青岛市市北区现代服务业发展中的问题研究[D].中国海洋大学,公共管理（专业学位）,2014,硕士.
[18]梁山,石飞光章,石为人,鲜晓东.离散时间MIMO系统极限零点的稳定性[J].自动化学报,2007,04:439-442.
[19]李勃燊.江苏联通移动新闻采编产品的营销策略研究[D].华东理工大学,高级工商管理（专业学位）,2014,硕士.
[20]魏轲.南极酵母重金属胁迫的蛋白组学及MT基因的克隆和分析[D].哈尔滨工业大学,微生物学,2014,硕士.
[21]张晓东,雍虎,熊建,陈显宝,刘波,田学军,郝建彬.如何看待中国快递市场对外开放[J].中国物流与采购,2014,22:32-34.
[22]吴海滨.闭环供应链混合回收渠道模型及协调研究[D].重庆大学,管理科学与工程,2014,硕士.
[23]邱阳.电子签章认证的相关法律问题研究[D].山东大学,法律（专业学位）,2013,硕士.
[24]吕一.基于现代优化算法的曲线识别方法[D].中北大学,应用数学,2013,硕士.
[25]贺铿.阿托伐他汀钙序贯治疗延缓冠状动脉病变进展的临床研究[D].南方医科大学,内科学,2012,硕士.
[26]武思敏.二语在三语习得中的迁移[D].吉林大学,英语语言文学,2004,硕士.
[27]史娅霖.高中历史教科书中科技人物詹天佑的编写研究[D].华中师范大学,学科教学,2013,硕士.
[28]周典恩.福建新教教会医院之研究（1842—1949）[D].福建师范大学,专门史,2004,硕士.
[29]邱祥龙.基于定性比较分析的政府新药政策研究[D].哈尔滨工业大学,行政管理,2013,硕士.
[30]陈长园.基于项目的学习促进学生有效合作的个案研究[D].华中师范大学,现代教育技术,2014,硕士.
[31]汪犇.microRNA-217对人皮肤成纤维细胞衰老的作用及机制研究[D].中南大学,临床医学,2014,博士.
[32]张保稳,何华灿,冯红伟.发现模糊状态演化模式[J].计算机科学,2002,03:11-13.
[33]温倩.住宅小区物业管理系统的设计与实现[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[34]白方兵.螺旋桨/机翼气动干扰的数值模拟研究[D].南京航空航天大学,飞行器设计,2014,硕士.
[35]张雅丽.酮症酸中毒对1型糖尿病患儿认知功能的影响[D].吉林大学,临床医学,2013,硕士.
[36]王永旺.室内运动目标姿态测量系统的研究与实现[D].中国科学院研究生院（西安光学精密机械研究所）,光学工程,2013,硕士.
[37]李志娟.我国承接欧美服务外包影响因素研究[D].河北经贸大学,国际贸易学,2012,硕士.
[38]叶灿.中国民营企并购国有企研究[D].华中科技大学,国际商务,2013,硕士.
[39]邵俊俊.基于工作过程系统化高职课程实践研究[D].宁夏大学,教育（专业学位）,2014,硕士.
[40]姜奉祥.面向现代制造的CAD/CAM应用研究[D].浙江大学,计算机应用,2004,硕士.
[41]沈显璞,张杰,吴宝榕,雷廷权.双相钢中马氏体相变热力学分析[J].金属学报.1989(02)
[42]季海波,陈志福,奚宏生,王俊.随机Marino-Tomei系统的输出反馈鲁棒自适应跟踪[J].自动化学报,2004,05:796-800.
[43]张新香.综合数据仓库的新型决策支持系统[J].现代计算机(专业版),2004,02:65-67+72.
[44]魏伟,王大伟,杨全红.锂空气电池中的碳基材料:优势与挑战[J].新型炭材料,2014,04:265-271.
[45]何觅然.操纵证券市场罪相关问题研究[D].吉林大学,刑法学,2013,硕士.
[46]邓毅凡.盈余管理水平与企业投资效率的关系研究[D].广东外语外贸大学,世界经济,2013,硕士.
[47]汪琦.评价和销量信息线索影响消费者在线购物行为的机理[D].浙江大学,管理科学与工程,2012,硕士.
[48]刘涵.指甲花及其制品的质量控制研究[D].北京中医药大学,中药学,2014,硕士.
[49]赵晖.针对电力机车牵引系统的单相整流有源电力滤波器[D].华中科技大学,电力电子与电气工程,2013,硕士.
[50]代忠美.高速车辆转向架稳定性机理研究[D].西南交通大学,载运工具运用工程,2014,硕士.