中国汉族人类风湿关节炎外周血基因表达谱研究

类风湿关节炎（rheumatoid arthritis, RA）是以持续性滑膜炎、全身性炎症及自身抗体的出现为特征的一组疾病，这类自身抗体以类风湿因子（rheumatoid factor, RF）和抗环瓜氨酸多肽（CCP）抗体（anti-cyclic citrullinated peptide antibody, anti-CCPantibody）为代表。RA是最常见的风湿性疾病之一，在我国，RA患病率为0.2-0.4%。该病病程迁延，致残率高，严重影响患者的生活质量，同时高额的医疗费用给患者和国家造成沉重的经济负担。因此，提高RA发病机制的研究水平，从而推进临床诊断和治疗水平是降低RA致残率，改善患者生活质量和保护社会劳动力的重要方法，是国家亟待解决的重大问题。 虽然RA发病的发病机制迄今未明，但是公认为RA的发病是由遗传和环境因素（如：吸烟）共同发挥作用而致病。RA是在环境因素，如吸烟和感染的影响下发生免疫紊乱导致活化的CD4+T细胞和MHC-Ⅱ型抗原递呈细胞浸润滑膜，并在细胞因子的参与下发病。遗传学因素可以解释RA大约50%的易感性。迄今为止，通过基因组学的研究，特别是全基因组关联分析研究（GWAS）的引入后，已发现40余个RA易感基因，为疾病的发病机制研究和风险预测提供了依据。同时，对这些发现的基因进行进一步挖掘，从群体遗传学中易感位点过渡到实际运用于诊疗的临床分子标记，完成真正的从实验室成果向临床诊疗手段的转化成了当今RA研究的当务之急。作为基因水平和蛋白水平的桥梁，转录水平的研究显得非常重要。 目前国际上在转录组水平研究RA的文章数量不菲，但是有以下局限性：一，目前尚无中国汉族人转录组水平的研究。绝大多数研究是基于欧裔人群的数据，而既往的遗传学研究证实类风湿关节炎是具有种族特异性的疾病。二，大多数研究使用的方法为芯片杂交技术。基因芯片技术需要根据已知基因来设计探针，只限用于已知序列，因为无法检测新的mRNA。基因芯片技术的灵敏度也非常有限,难以检测低丰度的目标和重复序列，这就需要更多的样本来进行研究。同时，它也很难检测出异常转录产物。RNA测序（RNA sequencing, RNA-Seq）技术为解决第二个问题提供了良好的平台。相对于芯片杂交技术，RNA测序技术能够提供更精确的数字化信号，更高的检测通量以及更广泛的检测范围。 本研究将运用RNA测序技术，在mRNA水平探索RA免疫遗传学发病机制：对已发现的易感基因进行进一步挖掘，力图从群体遗传学中易感位点过渡到实际运用于诊疗的临床分子标记，完成真正的从实验室成果向临床诊疗手段的转化。 第一部分基于RNA-Seq技术类风湿关节炎外周血基因表达谱 目的：基于RNA-seq技术检测中国汉族人类风湿、关节炎患者与正常对照外周血差异表达的基因，探讨汉族人RA的发病机制。 方法：提取5例处于活动期的抗CCP抗体阳性女性类风湿关节炎患者与性别匹配5例正常对照（Normal Control, NC）的外周血单个核细胞mRNA，使用Truseq RNAsample preparation kit（Illumina）构建测序文库，应用HiSeq2000系统进行测序分析，使用Cufflinks软件进行RNA表达量的计算，Cuffdiff软件默认参数进行显著性差异表达分析（P值<0.05并且表达差异2倍以上）。表达差异显著的基因进行聚类分析、GO分析和KEGG Pathway分析。通过测序结果与既往类风湿关节炎转录组报道文献比对、与既往报道的全基因关联分析研究（GWAS）基因比对确定下一步进行验证的候选基因。 结果：32106个转录本在RA组或NC组中至少有一个样本表达，P值<0.05且表达差异2倍以上的转录本274个（138个转录本表达上调；136个转录本表达下调）。聚类分析显示显著差异表达的转录本能够将测序样本明确分为2类；274个差异表达的转录本进行GO分析，富集通路为：神经系统相关通路：轴突、神经元的发育，突触的装配、结构和激活；免疫系统相关通路：多肽、多聚糖抗原的MHC-II分子处理和递呈，T细胞的分化、激活和增殖，单核细胞的增殖，细胞对细菌成分（如内毒素）的反应；蛋白泛素化等等。根据比对既往文献及GO分析的结果选取43个基因作为验证试验的候选基因。 结论：测序实验差异基因的GO分析提示固有免疫和适应性免疫共同参与类风湿关节炎的发病发展；环境因素细菌感染特别是内毒素在类风湿关节炎发病发展中的作用；蛋白泛素化在发病机制中的研究值得深入探讨。 第二部分大样本验证类风湿关节炎相关基因 目的：在大样本人群中验证测序结果和既往报道类风湿关节炎相关基因，寻找可应用于临床诊治的分子标记。 方法：收集类风湿关节炎患者外周血101例，正常对照55例，并登记相关临床资料。通过Trizol法提取总RNA，逆转录为cDNA。使用实时荧光定量PCR技术首先在小样本（15例RA患者和15例健康对照）中进行验证中检测43个基因表达情况；最终筛选9个基因在大样本（86例RA和40例NC）中进行基因表达情况的检测。根据临床资料对患者进行分层及关联分析，统计基因表达的差异。 结果：小样本验证实验中证实以下基因在病例组和正常对照组中差异表达显著：GPRIN3(р<0.01)、CXCR2(р<0.01)、FCGBP(р<0.05)、CD14(р<0.05)、RCAN1(р<0.05)、IKZF3(р<0.05)、BIRC3(р<0.05)、AFF3(р<0.05)、TSC1(р<0.05)、TYK2(р<0.05)、PADI4(р<0.05)。选取其中感兴趣的9个基因在大样本实验中进行验证：TSC1(р<0.01)、CD14(р<0.01)、AFF3(р<0.01)、BIRC3(р<0.01)、GPRIN3(р<0.05)、FCGBP(р<0.05)、RCAN1(р<0.05)差异表达显著，IKZF3、CXCR2表达无差异。其中，抗环瓜氨酸抗体阴性患者TSC1表达下降较抗环瓜氨酸抗体阳性患者明显(р<0.05)。中晚期患者CD14表达下降较早期患者下降明显(р<0.05)。 结论：TSC1和BIRC3是首次报道的与汉族人RA相关的基因；除CD14、AFF3基因介导的免疫相关性通路，BIRC3、TSC1基因参与的蛋白异寡聚化及BIRC3介导的凋亡与抗凋亡在类风湿关节炎病理过程中可能具有重要作用。