莱姆病是1种由硬蜱传播的伯氏疏螺旋体（borrelia burgdorferi，Bb）引起的多阶段、多系统感染性疾病，常见于夏季，影响人类和动物。感染初期，患者可能出现圆形红色皮疹，此时抗生素治疗通常有效。然而，若未能及时诊断和治疗，Bb可进一步侵犯血液和神经系统，导致神经莱姆病（lyme neuro borreliosis，LNB），这是莱姆病的1种严重并发症，影响约10%的患者^[1]。LNB的临床表现多样，包括神经根神经炎、淋巴细胞性脑膜炎等，而脑炎和脊髓炎则较为罕见^[2]。

诊断LNB时，需检测患者的血清和脑脊液。感染后3周内，可检测到特异性IgM抗体；6周后，特异性IgG抗体出现。值得注意的是，即使在康复后，这些抗体也可能持续存在，因此血清抗体阳性并不能单独作为感染的确诊依据^[3]。

目前，LNB的治疗推荐包括静脉滴注头孢曲松或口服多西环素，但需注意孕妇禁用^[4−5]。尽管如此，某些情况下抗生素治疗可能效果不佳，或患者可能出现持久的神经损伤或认知障碍。鉴于传统诊断和治疗方法的局限性，本研究采用生物信息学分析，旨在探索新的诊断和治疗策略。

笔者通过Pubmed回顾了相关文献，发现色氨酸-天门冬氨酸重复序列36基因（WD40-repeat 36，WDR36）和WRN基因组织表达定位多处相同如：脑、脊髓、淋巴和脾脏等多个组织中表达显著^[6−8]。因此猜测WDR36也与该模型的炎症反应存在基因途径联系。这些发现促使笔者将研究放在WDR36和WRN基因上，以期揭示它们在神经莱姆病炎症反应中的作用。

色氨酸-天门冬氨酸重复序列36基因（WD40-repeat 36，WDR36）最早是由Monemi在研究2个大型高加索家系时发现，定位GLC1G，染色体5q22.1^[9]。WDR36大小为34.7 kB，有23个外显子，编码含950个氨基酸的蛋白，其编码产物含有4个保守域就包括鸟苷酸结合蛋白WD40重复域^[10]。其中WDR36所属的WD重复蛋白家族的成员参与多种细胞过程，如细胞周期、信号转导、细胞凋亡和基因调控^[11]，WD40重复（WDR）结构域是人类蛋白质组中最丰富的蛋白质相互作用结构域之一。迄今为止，已注释了360多个结构域，这些结构域是参与多种信号通路的多蛋白复合物的重要亚基，例如：DNA损伤传感和修复、基因表达和染色质组织的表观遗传调控、泛素信号传导和蛋白质降解、细胞周期和免疫相关途径^[12]。在功能方面上WDR36 基因编码的WDR36蛋白是一个多功能蛋白，其在加工rRNA，维持细胞核正常形态及脑、眼、消化道的正常发育方面发挥重要作用。另外，WDR36基因所编码的蛋白是T细胞激活蛋白，参与了T细胞的活化并且与白细胞介素（interleukin-2，IL-2）有高度的相互调节作用^[13]。据报道，IL-2主要是由CD4⁺T细胞产生，可以刺激自然杀伤细胞（natural killer cell，NK）的增殖分化，诱导细胞毒性淋巴细胞和淋巴因子活化杀伤细胞（lymphokine activated killer cells，LAK），促进B细胞增殖和分泌抗体和激活巨噬细胞^[14]。IL-2对机体的免疫应答和抗病毒感染有重要作用。因此，WDR36基因对机体的免疫起着必要的关系。

在 1992 年首先利用连锁分析技术将 Werner 基因（Werner syndrome gene，WRN）定位在8号染色体（8p12），后1996年应用传统的定位克隆方法确定了WRN基因，并对其进行了测序分析WRN基因跨度超过250 kB有35个外显子，其中34个为编码外显子，由1432个氨基酸残基组成、分子质量为160 kD的WRN蛋白（WRNp）编码^[15]。WRN有别于其它RecQ家族解旋酶的一点是其N端具有外切酶结构域WRN-exo，同时C端包含很长的低保守序列，据报道，这些保守性低的序列能够与多种DNA代谢相关的蛋白质有相互作用^[16-17]。WRN是DNA解旋酶RecQ家族的成员之一，参与多种生物学过程，包括DNA损伤修复、端粒维持、自噬、凋亡、基因组维护等^[18−19]。WRN基因功能障碍易导致过早衰老，使人类容易患上各种类型的肿瘤和癌症^[20]。

研究表明，WDR36基因在T细胞活化和IL-2产生中扮演关键角色，而WRN蛋白则参与调控机体的天然免疫反应^[14，17]。鉴于二者均与炎症反应相关，笔者推测它们可能在LNB的发病机制中发挥作用。LNB的主要病理过程涉及Bb激活神经胶质细胞，进而释放炎性因子，导致中枢和周围神经系统的损伤或功能紊乱。为了验证WDR36和WRN在Bb感染引起的LNB中的潜在作用，笔者进行实验验证。

1.   材料与方法

1.1   主要试剂

1×磷酸盐缓冲液（phosphate buffered saline，PBS）由Thermo Fisher Scientific（Gibico）公司供应、10×TBST由北京索莱宝公司提供、HMC3细胞购买自上海中乔新舟生物科技有限公司、MEM培养基由上海中乔新舟生物科技有限公司提供、PCR剂盒由宝生物工程（大连）有限公司提供、Borrelia burgdorferi 4680菌株从德国菌株保藏中心（DSMZ）购得、RIPA裂解液由北京索莱宝科技有限公司提供、RNAiso plus由宝生物工程（大连）有限公司提供、RNase-free water由宝生物工程（大连）有限公司提供、SYBR荧光定量试剂由宝生物工程（大连）有限公司提供、Tris base由生工生物工程（上海）股份有限公司提供、二甲基亚砜（DMSO）由德国Sigma公司提供、甘氨酸由北京鼎国生物技术有限公司提供、甲醇由生工生物工程（上海）股份有限公司提供、氯仿由天津市化学试剂一厂提供、琼脂粉由Thermo Fisher Scientific（OXOID）公司提供、四甲基乙二胺（TEMED）由生工生物工程（上海）股份有限公司提供、胎牛血清由Thermo Fisher Scientific（Gibico）公司提供、台盼蓝由生工生物工程（上海）股份有限公司提供、胰蛋白胨由赛默飞世尔科技（中国）有限公司提供、引物合成由北京擎科生物新业生物技术有限公司提供。

1.2   研究方法

1.2.1   转录组分析

本研究利用与人类基因组相似度高达93%以上的恒河猴作为实验模型，以模拟人类对LNB的反应^[21]。本研究通过昆明医科大学动物实验伦理审查委员会的审批（kmmu20211590）。实验方法基于恒河猴脑额叶组织与Bb共培养测序得到的转录组学数据。从转录组学数据中筛选了3个时间点的数据集，采用阈值|log Fc|>1和P < 0.05进行显著性分析，并利用Venn图得到了23个差异基因并直观展示不同时间点基因表达的差异性。

本研究通过STRING蛋白互作数据库构建了基因互作网络，并筛选出与免疫反应相关的差异表达基因（differentially expressed genes，DEGs）。

进一步利用Metascape工具，设定阈值为P < 0.01、最小计数≥3且富集因子 > 1.5的标准进行了GO功能富集分析。

为验证WDR36和WRN基因在Bb感染神经莱姆病的作用，发现小胶质细胞作为中枢神经系统的一种巨噬细胞，具有强大的吞噬功能和免疫调节作用。在静息状态下小胶质细胞缺乏吞噬功能，当中枢神经系统受到损伤或者炎症感染时可被激活，活化后使得炎症反应不断扩大。因此本研究以人脑小胶质细胞3（human brain microglia 3，HMC3）细胞株为模型，对实验组进行Bb 接种，分为MOI=1和MOI=10，分别在细胞铺板感染后6 h，12 h和24 h 3个时间点收集细胞悬液，用trizol法提取细胞悬液中RNA，通过qPCR法检测目的因子WDR36和WRN的mRNA相对表达情况。实验流程：首先利用细胞培养技术培养HMC3细胞株并对培养细胞进行计数、铺板、刺激；其次设立PBS对照组，实验组为Bb感染，分为MOI=1和MOI=10；另外分别在细胞铺板感染后6 h，12 h，24 h 3个时间点收集细胞的RNA悬液；最后用trizol法提取细胞总RNA，通过qPCR （SYBR Green） 法检测目的因子WDR36和WRN的mRNA相对表达情况。

1.2.2   转录组结果验证

（1）目的基因特异性序列查询：登录“NCBI”网站中的Gene网页（https://www.ncbi.nlm.nih.gov/gene/）查找目的基因特异序列。

（2）引物设计：登录引物设计网页“Primer3 version 4.0.0”（http://bioinfo.ut.ee/primer3-0.4.0/），将目的基因特异性序列输入，选择引物长度范围、物种等，得到需要的目的基因引物。

（3）引物特异性检验：将设计好的引物输入 “NCBI”网站“Primer-BLAST”数据库（https://www.ncbi.nlm.nih.gov/tools/primer-blast/），验证引物特异性。 挑选特异性良好的引物交由北京擎科生物新业生物技术有限公司进行引物序列合成。引物序列如：

WDR36 基因：上游引物（5'-3'）：ACTCGCCTTGGATGACTTCT、下游引物（5'-3'）：GCAATCCATCGCAGCACTTA。

WRN 基因：上游引物（5'-3'）：AGAAAGGCATGTGTTCGGAG、下游引物（5'-3'）：GCCACTCCATGTCAAATCCC。

1.2.3   qPCR 相对定量目的基因拷贝数

（1）使用TaKaRa公司SYBR荧光染料qPCR反应试剂盒进行试验，每个样本的目的基因（WDR36和WRN）设置2个复孔，在冰上完成反应体系配置和加样操作。反应体系：RNase-freed H₂O：8.5 μL；SYBR Premix ExTaqⅡ：12.5 μL；PCR Forward Primer：1 μL；PCR Reverse Primer：1 μL；cDNA模板：2 μL；总计25 μL。

目的基因（WDR36和WRN）qPCR反应条件：预变性：95 ℃，15 s；变性：95 ℃，5 s和退火：60 ℃，30 s共循环40次、Dissociation stage：65～95 ℃。

（2）最终结果判定标准：通过Bio-Rad CFX96TM软件得到qPCR（SYBR Green）检测的结果，熔解曲线平稳，1个熔解峰，扩增产物单一，判断引物特异良好。扩增曲线平滑，阈值循环数（Ct值）。目的基因（WDR36、WRN）Ct值为15～35次，2个复孔间Ct值相差小于0.1，此qPCR的检测结果可靠。

（3）反应结束后，读取Ct值，最终目的基因的mRNA相对表达量用2^-ΔΔCt值表示，本方法将对照组（PBS组）基因mRNA表达量设置为1，其他组mRNA表达量为对照组的倍数。计算公式：^ΔCt = Ct目的基因－Ct内参基因。

1.3   统计学处理

利用Graphpad Prism 6.0软件进行作图和统计分析，实验数据采用双因素方差分析法（Two-way ANOVA）进行统计学分析作图。P < 0.05表示差异有统计学意义。^ΔΔCt =^ΔCt实验组-^ΔCt对照组。

2.   结果

2.1   Venn图筛选出23个差异表达基因（DEGs）

以非人灵长动物（恒河猴）作模型，利用团队培养并测序的转录组学数据为背景，随机筛选了3个时间的数据集，设定阈值为|log Fc|>1，P value< 0.05，通过绘制Venn图筛选出23个差异表达基因（differentially expressed genes，DEGs），见图1。

图  1  23个差异表达基因（DEGs）的筛选

Figure  1.  Screening of 23 differentially expressed genes （DEGs）

下载: 全尺寸图片 幻灯片

2.2   免疫反应相关的DEGs的相互作用关系网络

以非人灵长动物（恒河猴）作模型，利用团队培养并测序的转录组学数据为背景，绘制与免疫反应相关的DEGs的相互作用关系网络图，其中13个显著的免疫差异基因被标红，包括CFH、OGN、BTAF1、WDR36、TMEM87A、ALCAM、MFN1、WRN、CENPC、SF3B1、CTNNAL1、ZBED5和FOLR2。WDR36和WRN的位置用红色方框表示，见图2。

图  2  免疫反应相关DEGs的相互作用关系网络

Figure  2.  The interaction network of immune response related DEGs

下载: 全尺寸图片 幻灯片

2.3   GO富集分析

参与的生物过程（biological process，BP）：细胞对外界刺激；基因的分子功能（molecular function，MF）：三磷酸核糖核苷磷酸酶活性，见图3。

图  3  GO富集分析含量

Figure  3.  GO enrichment analysis content

下载: 全尺寸图片 幻灯片

2.4   主要差异表达基因

结果显示，TMEM87A、WRN和FOLR2与细胞对外界刺激的反应（BP）相关；BTAF1、MFN1和WRN与三磷酸核糖核苷磷酸酶活性相关，见表1。其中，WRN基因出现了富集重合显示出了重要的研究潜力。

表  1  主要差异表达基因BP（TMEM87A、WRN、FOLR2）和MF（BTAF1、MFN1、WRN）

Table  1.  BP（TMEM87A、WRN、FOLR2）and MF（BTAF1、MFN1、WRN） selected genes

基因	GO富集BP	GO富集MF
CFH OGN BTAF1 WDR36 TMEM87A ALCAM MFN1 WRN CENPC SF3B1 CTNNAL1 ZBED5 FOLR2	未选择 未选择 未选择 未选择 选中 未选择 未选择 选中 未选择 未选择 未选择 未选择 选中	未选择 未选择 选中 未选择 未选择 未选择 选中 选中 未选择 未选择 未选择 未选择 未选择

下载: 导出CSV 

| 显示表格

2.5   WDR36和WRN基因表达的实验验证

目的基因WDR36在HMC3细胞被Bb感染后，6 h，12 h，24 h，PBS组（阴性对照组）和Bb感染组（实验组）的细胞做统计学分析，发现Bb MOI=1为适宜感染浓度，若Bb MOI=10，细胞在24 h耐受不佳，见图4。所以本研究采用Bb MOI=1感染组（实验组）和PBS组（阴性对照组），不同时间点（6 h，12 h，24 h），通过对目的基因WDR36的检测，2组实验对照发现WDR36的表达量在24 h与12 h对比有显著的上调（P < 0.01），且6 h的PBS组与试验组基因表达含量对比上升但不明显（P > 0.05），12 h、24 h上升明显（P < 0.01），见图5。

图  4  组内HMC3在不同时间点WDR36表达量的比较

^*P < 0.05，^**P < 0.01。

Figure  4.  Comparison of WDR36 secretion of HMC3 at different time points in the group

下载: 全尺寸图片 幻灯片

图  5  组间HMC3在不同时间点WDR36表达量的比较

^**P < 0.01。

Figure  5.  Comparison of WDR36 secretion of HMC3 at different time points between groups

下载: 全尺寸图片 幻灯片

目的基因WRN在HMC3细胞被Bb感染后，通过对6 h，12 h，24 h，PBS组（阴性对照组）和Bb感染组（实验组）的检测分析，表明Bb MOI=1为适宜感染浓度，若Bb MOI=10，细胞在24 h并未呈阶梯式增长，见图6。所以本研究采用Bb MOI=1感染组（实验组）和PBS组（阴性对照组），不同时间点（6 h，12 h，24 h），通过对目的基因WRN的检测，2组实验对照发现WRN的表达量在24 h与12 h对比有显著上调（P < 0.01），且12 h、24 h的实验组分别比其PBS组基因表达含量呈上升趋势（P < 0.01），见图7。

图  6  组内HMC3在不同时间点WRN表达量的比较

^**P < 0.01。

Figure  6.  Comparison of WRN secretion of HMC3 at different time points in the group

下载: 全尺寸图片 幻灯片

图  7  组间HMC3在不同时间点WRN表达量的比较

^**P < 0.01。

Figure  7.  Comparison of WRN secretion of HMC3 at different time points between groups

下载: 全尺寸图片 幻灯片

3.   讨论

在Bb感染的人脑胶质细胞HMC3细胞株中，观察到WDR36和WRN基因的转录水平在感染后显著上调。特别是在感染后的6 h和12 h，这2个基因的表达水平显著增加，而在24 h时，这种上调趋势有所减缓。这一现象表明，在Bb感染的早期阶段，WDR36和WRN基因可能在LNB的病理过程中发挥作用。

进一步的分析显示，在适宜的感染浓度（Bb MOI=1）下，24 h感染组的WDR36和WRN基因表达量明显高于12 h感染组，这表明随着感染时间的延长，这2个基因的表达水平可能会进一步增加。然而，笔者也注意到在高感染浓度（Bb MOI=10）下，长时间的培养（24 h）会导致细胞出现不耐受现象，这可能影响基因表达的稳定性。

在对比实验组与PBS对照组时，发现在6 h的感染条件下，WRN基因的表达上调幅度比WDR36更为显著。这表明WRN基因在Bb感染的早期反应中可能扮演更重要的角色，因此，WRN基因可能是1个更适合进行后续实验分析和研究的靶点。

LNB患者在接受抗生素漫长的治疗后会出现治疗后莱姆病治疗后综合征（post-treatment lyme disease syndrome，PTLDS），伴有疲劳、肌肉关节痛和认知障碍，当前未知此原因^[22]，Bb引起脑膜血管型伴脑血管梗死、萎缩型伴脑膜脑炎、皮质萎缩和胶质增生，可导致发展更为迅速的痴呆^[23]，中枢神经系统外的感染导致中枢神经系统内的免疫反应，可能与缓慢进展型痴呆症有关^[24]，因此推测WRN的基因稳定与患者抗生素治疗后出现过早衰老和痴呆、认知障碍、疲劳等症状相关。

《中国卫生健康统计年鉴2021》数据显示，2020年我国城市居民冠心病死亡率为126.91/10万，农村为135.88/10万^[25]，冠心病的治疗在我国仍是1项难题。据报道，在通脉养心丸治疗冠心病的研究中，通过筛选差异表达基因和检索蛋白相互作用时就有WDR36基因的存在^[26]。推测如遇到同时患有神经莱姆病和冠心病的患者使用通脉养心丸可能对2种疾病能起到治疗效果。

据报道，WDR36和WRN基因都与炎症反应有关系，在现阶段莱姆病炎症机制的研究指出伯氏疏螺旋体会致神经细胞凋亡^[27]，当Bb影响神经系统时，神经疏螺旋体病，免疫调节性炎症级联反应可触发神经胶质细胞凋亡，导致各种神经学表现^[28−30]。凋亡是种非炎症、程序性细胞死亡，而与凋亡密切相关的另1种炎症介导的程序性细胞死亡则是焦亡，细胞焦亡研究中它的经典途径是依赖于半胱氨酸蛋白酶-1（Caspase-1）与IL-1β和IL-18分泌，从而介导促炎效应^[31]在焦亡和凋亡细胞中都存一定作用机制^[32]，目前没有报道细胞焦亡与LNB的关系，所以怀疑神经莱姆病导致的炎症不仅只有凋亡因素，也存有细胞焦亡的可能性。

LNB的临床表现多样性受Bb的基因型和宿主细胞反应性的影响。本研究通过转录组学分析，结合人脑胶质细胞HMC3模型的实验验证，推测WDR36和WRN基因在Bb感染引起的神经莱姆病中表达上调，这可能与它们在神经炎症和细胞应激反应中的功能有关。此外，考虑到这两个基因在感染过程中的动态表达模式，它们可能作为潜在的治疗靶点，用于干预或缓解LNB的病程。未来的研究可以进一步探索这些基因在LNB发病机制中的具体作用，以及它们作为治疗靶点的潜力。

通过血清学检查确诊LNB患者的主要感染病原体为Bb时，WDR36和WRN基因的上调表达可能成为疾病进程中的生物标志物。这不仅为揭示LNB的分子机制提供了新的视角，也为针对Bb感染的新型治疗策略的开发提供了可能的分子靶点。未来研究将深入探讨这些基因在不同Bb基因型和宿主细胞背景下的功能，以及它们在LNB治疗中的潜在应用，以期为临床治疗提供更精准的指导。