论文部分内容阅读
第一部分基于基因表达谱的儿童噬血性淋巴组织细胞增生症候选基因筛选及相关实验研究目的应用生物信息学方法,筛选新的儿童噬血性淋巴组织细胞增生症(He-mophagocytic lymphohistiocytosis, HLH)候选疾病基因并进行实验验证,探讨儿童HLH发病和预后的分子机制。方法(1)从基因表达综合数据库(Gene expression omnibus, GEO)中下载GSE26050儿童HLH基因芯片数据集,然后在R语言软件包中进行差异表达基因(Differentially expressed genes, DEGs)分析,将获得的DEGs定义为“检测基因集”;采用Genecards和Fable文献挖掘已知儿童HLH疾病基因,将其定义为“训练基因集”;最后,利用Toppgene在线分析工具筛选出儿童HLH候选疾病基因。(2)以100例HLH患儿为病例组,146例健康儿童为对照组,采集外周静脉血,采用多重聚合酶链式反应(Polymerase chain reaction, PCR)和SNaPshot技术对部分候选基因单核苷酸多态性(Single nucleotide polymorphism, SNP)进行检测。(3)用Haploview v.4.2软件进行各SNP位点哈-温平衡(Hardy-Weinberg eq-uilibrium, HWE)和相互间的连锁不平衡(Linkage disequilibrium, LD)检验。用SPSS17.0软件分析各位点基因型频率、等位基因频率分布差异与儿童HLH发病的相关性。风险度采用比值比(Odds ratio, OR)及其95%可信区间(Confidence interval,CI)表示。用PHASE2.1软件进行单体型的构建和分析。P小于0.05视为差异有统计学意义。(4)对入选HLH患儿进行病历资料收集并随访。以确诊HLH为观察起点,死亡日期作为观察终点,失访、末次随访时仍存活(包括缓解、疾病活动、疾病复发的患儿)的病例则确定为截尾病例。以有意义的SNPs不同基因型作为研究因素,采用Kaplan-Meier曲线进行生存分析,应用Log-rank检验比较生存曲线。以P<0.05为差异有统计学意义。结果(1)获得一个含86个基因的“检测基因集”和一个含73个基因的“训练基因集”;采用Toppgene共获得45个儿童HLH候选疾病基因,其中表达上调的基因42个,表达下调的基因3个,为TNFRSF17、CX3CR1和CCR2。(2)TNFRSF17基因多态性分析结果:病例组rs3743591GG基因型频率高于对照组,差异有统计学意义(10.0% vs.1.4%,P=0.002),提示rs3743591GG基因型是儿童HLH发病的易感因素(OR=8.000,95% CI=1.714-37.349);病例组rs3743591G等位基因频率高于对照组,差异有统计学意义(20.5% vs.10.6%,P=0.002),提示rs3743591G等位基因是儿童HLH发病的易感因素(OR=2.171,95% CI=1.308-3.603)。病例组rs2017662AA基因型频率高于对照组,差异有统计学意义(11.0% vs.1.4%,P=0.001),提示rs2017662AA基因型是儿童HLH发病的易感因素(OR=8.899,95% CI=1.928-41.082);病例组rs2017662A等位基因频率高于对照组,差异有统计学意义(21.0% vs.10.3%,P=0.001),提示rs2017662A等位基因是儿童HLH发病的易感因素(OR=2.322,95% 0=1.396-3.860)。(3) CX3CR1基因多态性分析结果:病例组rs2853712TT基因型频率高于对照组,差异有统计学意义(72.0% vs.50.0%,P=0.001),提示rs2853712TT基因型是儿童HLH发病的易感因素(OR=2.571,95% CI=1.493-4.430)。(4)CCR2基因多态性分析结果:病例组rs1799865TT基因型频率高于对照组,差异有统计学意义(39.0% vs.20.5%,P=0.002),提示rs1799865TT基因型是儿童HLH发病的易感因素(OR=2.472,95% CI=1.401-4.363);病例组rS 1799865T等位基因频率高于对照组,差异有统计学意义(53.0% vs.37.0%,P=0.000),提示rs1799865T等位基因是儿童HLH发病的易感因素(OR=1.921,95% CI=1.333-2.769)。(5)单体型分析结果:单体型A-T-G-G-G-C-T-T-T-T-T-A-G-A-C (SNPs顺序:rs3743591-rs11570151-rs2017662-rs2071336-rs2669850-rs17793056-rs13088991-rs13062158-rs2853712-rs2669841-rs2853711-rs3762823-rs309296 3-rs3092962-rsl 799865)在病例组出现的频率高于对照组,差异有统计学意义(5.7% vs.1.3%,P=0.003),提示单体型A-T-G-G-G-C-T-T-T-T-T-A-G-A-C是儿童HLH发病的易感因素(OR=5.025,95% CI=1.530-16.505).单体型G-T-A-G-G-T-C-C-T-C-G-A-G-A-C (SNPsJl顺序:rs3743591-rs11570151-rs2017662-rs2071336-rs2669850-rs17793056-rs13088991-rs13062158-rs28537 12-rs2669841-rs2853711-rs3762823-rs3092963-rs3092962-rs 1799865)在病例组出现的频率高于对照组,差异有统计学意义(5.2%VS.0%,P=4.43×10-5),提示单体型G-T-A-G-G-T-C-C-T-C-G-A-G-A-C是儿童HLH发病的易感因素。(6)儿童HLH病例中,携带rs3743591GG基因型者中位生存时间小于携带rs3743591AA/AG基因型者,差异有统计学意义(86.0个月vs.118.0个月,P=0.009);携带rs2017662AA基因型者中位生存时间小于携带rs2017662GG/GA基因型者,差异有统计学意义(88.0个月vs.118.0个月,P=0.004);携带rs2853712TT基因型者平均生存时间与携带rs2853712CC/CT基因型者比较差别无统计学意义(96.5个月vS.108.6个月,P=0.124);携带rs1799865TT基因型者中位生存时间小于携带rs1799865CC/CT基因型者,差异有统计学意义(92.0个月vs.107.0个月,P=0.021)。结论(1)本研究应用生物信息学方法筛选出45个新的儿童HLH候选疾病基因,包括42个上调表达基因和3个下调表达基因。(2) TNFRSF17、CX3CR1和CCR2基因在儿童HLH中表达下调,提示其可能参与儿童HLH的发病机制。(3)TNFRSF17基因rs3743591、rs2017662, CX3CR1基因rs2853712, CCR2基因rs1799865多态位点与儿童HLH发病有关。(4) TNFRSF17基因、CX3CR1基因和CCR2基因多态位点构成的单体型A-T-G-G-G-C-T-T-T-T-T-A-G-A-C和G-T-A-G-G-T-C-C-T-C-G-A-G-A-C与儿童HLH发病有关。(5) TNFRSF17基因rs3743591、rs2017662, CCR2基因rs1799865多态位点有可能作为儿童HLH预后判断的分子标记。A-T-G-G-G-C-T-T-T-T-T-A-G-A-C和G-T-A-G-G-T-C-C-T-C-G-A-G-A-C与儿童HLH发病有关。第二部分儿童噬血性淋巴组织细胞增生症SH2D1A、XIAP基因变异和临床意义目的了解SH2区蛋白1A (SH2 domain protein 1A, SH2D1A)基因和X连锁凋亡抑制蛋白(X-linked inhibitor of apoptosis protein, XIAP)基因突变和序列变异在HLH患儿中的发生情况,探讨SH2D1A、XIAP基因变异与儿童HLH的关系。方法病例组为100例HLH患儿,对照组为100例健康儿童。采集外周静脉血,提取基因组DNA,采用PCR法扩增SH2D1A基因1-4外显子和XIAP基因1-6外显子编码区及其侧翼内含子区片段,PCR产物经ABI PRISM3730全自动DNA测序仪直接测序,然后与NCBI数据库中SH2D1A和XIAP基因的序列进行比对,获得基因变异信息,利用生物信息学方法预测SH2D1A和XIAP基因突变对蛋白质功能的影响。对基因突变的病例回顾性分析其临床资料,将其临床表型与基因型进行对比分析。结果(1)在HLH病例组中发现XIAP基因1个新的错义突变,30544 T>A,位于第1外显子上,导致所编码的氨基酸发生改变,Phe27Ile。该位点突变在对照组中未发现,突变发生率为1%(1/100)。(2)用PolyPhen2和SIFT软件预测XIAP基因错义突变功能,显示XIAP基因Phe27Ile突变有可能影响蛋白质功能。(3)在病例组和对照组中共发现5个SNPs位点,包括SH2D1A基因的rs72610640、rs182286559多态和XIAP基因的rs5956583、rs140997240、rs28382740多态;在病例组和对照组中比较了这5个SNPs位点的基因型、等位基因的频率分布,差异均无统计学意义,提示这些SNPs位点可能与儿童HLH发病无关。(4)临床特征分析:发生Phe27Ile突变的1例患儿无明显HLH阳性家族史和特殊既往史,13岁起病,患儿临床表现为HLH,按照HLH-2004方案接受化疗曾经得到缓解,但停药后复发,治疗效果差,预后不良,于16岁时死于多器官功能衰竭(Multiple system organ failure, MSOF)。结论本研究发现了1个新的XIAP基因错义突变,Phe27Ile,该突变可能与HLH患儿的发生和不良预后有关。