Clinical Characteristics of Varicella-zoster Virus Infection in 38 Immunocompromised Patients
-
摘要:
目的 探讨免疫低下人群感染水痘-带状疱疹病毒(varicella-zoster virus,VZV)的临床特征。 方法 收集昆明市第三人民医院2020年1月至2021年12月间,收治的免疫低下人群感染VZV38例为观察组,40例普通人群感染VZV为对照组,对其临床资料进行回顾性的分析,观察2组的心肌酶、肝功能、肾功能等评估脏器功能的指标;以及IL-6、IL-8、CD+4 T淋巴细胞等评估机体免疫的指标。 结果 将2组的肝功能、心肌酶进行对比,观察组GGT/AST/ALT、LDH/CK 均高于对照组,差异有统计学意义(P < 0.05)。将2组的CREA进行对比: 2组人群CREA,差异无统计学意义( P > 0.05)。2组的细胞因子IL-6、IL-8对比;观察组IL-6、IL-8均高于对照组,差异有统计学意义( P < 0.05)。2组CD4 +T淋巴细胞数对比,观察组CD4+T淋巴细胞数低于对照组,差异有统计学意义(P < 0.05)。 结论 免疫低下人群感染VZV,由于免疫被抑制,更容易诱发细胞因子风暴,病情更重,常易合并器官功能障碍甚至衰竭,病情进展快,若能早期识别,积极救治,能改善其预后。 Abstract:Objective To investigate the clinical characteristics of varicella-zoster virus (VZV) infection in immunocompromised population. Methods From January 2020 to December 2021 in the Third People’s Hospital of Kunming, 38 cases of VZV infection in immunosuppressed population were selected as the observation group, and 40 cases of VZV infection in the general population were selected as the control group. The clinical data were retrospectively analyzed, and the indexes of cardiac enzymes, liver function, kidney function and other organs function evaluation in the two groups were observed. And IL-6, IL-8, CD+4 T lymphocytes and other indicators’s data were analyzed to evaluate the body’s immunity. Results The GGT/AST/ALT and LDH/CK of the observation group were significantly higher than those of the control group (P < 0.05). There was no statistical significance in CREA between the two groups ( P > 0.05). IL-6 and IL-8 in the observation group were significantly higher than those in the control group ( P < 0.05). The number of CD4+T lymphocytes in the observation group was significantly lower than that in the control group ( P < 0.05). Conclusion Immunocompromised population with VZV infection are more likely to induce cytokine storm, and are more often complicated with organ dysfunction or even organ failure. The disease progresses faster and is more severe. Early identification and treatment can improve the prognosis. -
Key words:
- Immunocompromised population /
- Varicella-zoster virus /
- Clinical signs /
- Analysis
-
垂体腺瘤是最常见的神经内分泌肿瘤之一,发病率约17%,占中枢神经系统肿瘤的10%~25%[1]。垂体腺瘤按生物学特性分为非侵袭性垂体腺瘤和侵袭性垂体腺瘤。非侵袭性垂体腺瘤呈膨胀性生长,一般都具有清晰的边界。通常将侵袭性垂体腺瘤的生物学特性及其产生的影像学表现相结合作为诊断标准,表现为腺瘤侵入海绵窦、硬脑膜和颅骨等周围脑组织并包绕颈内动脉等重要结构[2]。侵袭性垂体腺瘤手术切除极其困难,术后容易残留和复发,给临床治疗带来极大困难[3]。因此,研究侵袭性垂体腺瘤的侵袭机制,对提高临床诊治效果具有重要意义。
基因组学发现,不到2%的人类基因被转录并翻译成蛋白质,而超过90%的基因被转录成非编码RNA (ncRNA)[4]。LncRNA在许多生命活动中起着重要的作用,如表观遗传调控、细胞周期调控、细胞分化调控等[5]。大量研究表明,lncRNA与肿瘤发生密切相关[6]。近年来,研究者开始研究lncRNA在侵袭性垂体腺瘤中的作用,郭辉等[7]发现垂体腺瘤细胞中IFNGAS1过表达可能通过下调上皮剪接调控蛋白2 (ESRP2)而调控垂体腺瘤的侵袭性作用;郭强等[8]研究发现,lncRNA1 C5orf66-AS1可抑制垂体功能腺瘤的发生和侵袭。虽然有侵袭性垂体腺瘤相关的lncRNA的报道,但lncRNA在侵袭性垂体瘤发病机制中的研究较少,尚无明确的lncRNA作为侵袭性垂体腺瘤的特异性标志物。
笔者使用RNA-seq分析来识别差异表达的lncRNA,以寻找可能作为分子标志物的lncRNA。鉴定差异表达的mRNA并建立lncRNA-mRNA网络来预测lncRNA与mRNA在侵袭性垂体瘤中可能的相互作用。GO和KEGG分析预测侵袭性垂体腺瘤相关的信号通路。lncRNA在垂体腺瘤中的表达有显著差异,异常表达的lncRNA可作为生物学标志物和特异性治疗靶点。这些发现可能为侵袭性垂体瘤今后的早期诊断和预后预测提供潜在的生物标志物。
1. 材料与方法
1.1 样本收集与测序准备
4例侵袭性垂体腺瘤和3例非侵袭性垂体腺瘤均来自昆明医科大学第一附属医院神经外科。所有病例均经影像学、手术及病理证实。本研究经昆明医科大学第一附属医院伦理委员会批准,所有参与者均已签署知情同意书。将手术切除的垂体瘤组织用生理盐水中反复漂洗,放入无菌、氯化的2.0 mL冷冻管。然后将样本浸泡在Trizol试剂中,放入−80 ℃冰箱保存,准备进行组织中的总RNA提取和后续分析。
1.2 总RNA提取、文库制备和测序
每个样本提取3 μg RNA作为RNA-seq准备材料。NanoPhotometer® spectrophotometer (IMPLEN,CA,USA)测量RNA纯度。Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies,CA,USA)测定RNA浓度。通过安捷伦生物分析仪2100系统(Agilent Technologies,CA,USA)的RNA Nano 6000检测试剂盒评估RNA完整性。用Epicentre Ribo-zero™ rRNA删除工具(Epicentre,USA)去除样本中的rRNA,并用乙醇沉淀法去除rRNA残留。按照制造商的说明使用rRNA-depleted RNA by NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB,USA)生成测序库。采用AMPure XP系统对产品进行纯化,并在安捷伦生物分析仪2100系统上对文库质量进行评价。根据制造商的说明,使用TruSeq PE Cluster Kit v3-cBot-HS (Illumia)在cBot集群生成系统上对索引编码的样本进行聚类。集群生成之后,在Illumina Hiseq 4000平台上进行测序。
1.3 RNA定量及差异表达
RNA表达水平由TPM(转录本/百万)进行表达量归一化处理,通过以下标准进行估计:标准化表达式 = Mapped readcount/总读取数×100000011[9]。使用 R语言DEGseq包对样本进行差异表达分析,差异lncRNA默认筛选条件为:P < 0.05,|log2(foldchange)| > 1。
1.4 GO和KEGG通路
为了研究侵袭性垂体瘤样本与非侵袭性垂体瘤样本失调的lncRNA的功能,笔者对其共表达的mRNA进行了功能分析,包括基因本体论(GO)分析和《京都基因与基因组百科全书》[10](KEGG)通路分析。用GOseq[11]对与lncRNA共表达的mRNA作GO富集分析,包括细胞组分(CC)、分子功能(BP)和生物学过程(MF)。KOBAS软件[12]检测KEGG通路中mRNA的富集分析情况。
1.5 lncRNA-mRNA共表达网络的构建
皮尔逊相关系数(PCC)被广泛用于测量线性相关[13]。一般来说,PCC值为1表示完全正线性相关,PCC值为−1表示完全负线性相关,PCC值为0表示无线性相关。使用差异表达的lncRNA和差异表达的mRNA来检测所有可能的lncRNA- mRNA相关关系。对于每个RNA,笔者使用其FPKM标准化表达值来计算相关性。对于每个lncRNA-mRNA,计算PCC来评估lncRNA和mRNA的表达是否在所有4个侵袭性垂体腺瘤或3个非侵袭性垂体腺瘤样本中一致线性相关。如果| PCC | > 0.95和P < 0.05,笔者认为lncRNA与mRNA呈线性相关,即具有共表达关系。使用Cytoscape ( http://www.cytoscape.org)将相关的lncRNA-mRNA对可视化为共表达网络。分子网络的拓扑性质包括度、介数和紧密度。因此,对这些常见的拓扑测量进行分析,以揭示网络的特性。Cytohubba是Cytoscape 插件,其计算的结果包括11种算法,选择最常用的算法度(Degree)筛选网络中排名前五的lncRNA[14]。
2. 结果
2.1 差异表达的lncRNA和mRNA
共检测到1997个lncRNA,其中表达具有显著差异的共35个,其中上调21个,下调14个,见表1。为了探讨lncRNA与mRNA之间可能的相互作用。对mRNA进行RNA-Seq分析,检测到463个差异表达基因的mRNA,其中229个表达上调,234个表达下调,差异最显著的前20个mRNA,见表2。差异表达基因的聚类热图,见图1。
表 1 差异表达的lncRNATable 1. Differentially expressed lncRNAs基因名 Log2(Foldchange) P 异常表达情况 RP11-252P19.2 7.110230855 0.025594461 上调 LINC00473 6.892052835 0.000316802 上调 STXBP5-AS1 5.640589233 0.007878966 上调 RP11-71H17.7 5.221628067 0.017764884 上调 TPTEP1 4.89206949 0.003097546 上调 RP11-262H14.3 3.952787834 0.03606433 上调 RP11-12O16.1 2.807343125 0.011324588 上调 RP11-111M22.3 2.28836899 0.002792906 上调 RP11-250B2.3 1.906609127 0.035774432 上调 FAM66A 1.848133787 0.002209197 上调 RP11-262H14.4 1.813382724 0.009868508 上调 CTD-3193K9.4 1.721455664 0.001307362 上调 MAFG-AS1 1.536359193 0.044840975 上调 SNHG8 1.461730762 0.049804872 上调 AC139100.4 1.448409941 0.045905713 上调 FAM27C 1.415494365 0.029287703 上调 AC005062.2 1.349936987 0.033475711 上调 RP11-646I6.5 1.272730756 0.046225249 上调 RN7SL674P 1.260917811 0.008860604 上调 RP11-662M24.2 1.188403274 0.037409712 上调 LRRC75A-AS1 1.042079322 0.002182732 上调 ZFAS1 −1.030295183 0.02204746 下调 TP53TG1 −1.093552638 0.027463784 下调 CTB-129O4.1 −1.167587296 0.041920971 下调 RP11-1020A11.1 −1.227419109 0.025518005 下调 AP006621.5 −1.265098387 0.004293786 下调 C5orf66 −1.323699343 0.024003471 下调 RP11-1114A5.4 −1.346165585 0.036015113 下调 AC017116.11 −1.355620383 0.018311217 下调 AC019181.2 −1.556768999 0.031834266 下调 UMAD1 −2.049223764 0.026775377 下调 AC005618.6 −2.103056949 0.012754031 下调 SOCS2-AS1 −2.450904018 0.010320293 下调 RP11-246K15.1 −3.394395932 0.006322703 下调 ARMCX5-GPRASP2 −3.509890448 0.042606052 下调 表 2 排名前20的差异表达的mRNATable 2. The top 20 differentially expressed mRNAmRNAs基因名 Log2(Foldchange) P 异常表达情况 SLC4A3 3.130642688 0.000067700 上调 ZNF444 1.555389418 0.000111993 上调 HDAC4 2.871048277 0.000146683 上调 MRO −3.344017691 0.000209461 下调 ECHDC1 10.92744203 0.000217674 上调 KIAA1958 −1.739496332 0.000224985 下调 TBC1D5 22.81214317 0.000254824 上调 NCALD 2.281098096 0.000325798 上调 LAMTOR5 −4.230692295 0.000366157 下调 TMPRSS6 7.265115275 0.000399492 上调 RPL26 −4.965389762 0.000462772 下调 RPSA 2.332207622 0.000464748 上调 CTSC 1.469599612 0.000464889 上调 MYO10 4.166108193 0.000472723 上调 ATP6AP2 3.197467642 0.000514505 上调 EPHX2 2.326018659 0.000616419 上调 PARP11 −1.579039067 0.000625344 下调 AGTRAP 1.045079801 0.000845595 上调 MXRA7 −1.391168544 0.000967649 下调 SLC29A2 −9.596231979 0.001050552 下调 2.2 lncRNA-mRNA共表达网络的构建与分析
用Pearson相关系数分析lncRNA与mRNA的相关性。35个差异表达的lncRNA与463个差异表达的mRNA共进行了35×463个相关分析。最后筛选出449对相关系数不同的碱基对,包括34个lncRNA和274个mRNA,见图2。列出Degree排名最高的5个lncRNA,见表3。
表 3 Degree排名前五的lncRNATable 3. The top 5 lncRNAs in Degree基因名 Degree FAM66A 52 TPTEP1 42 AC017116.11 33 RP11-1114A5.4 29 RP11-246K15.1 28 2.3 共表达网络中mRNA的GO和KEGG通路
侵袭性垂体瘤样本中失调的lncRNA与非侵袭性垂体瘤样本相比的功能,对与其共表达的mRNA进行 GO和KEGG通路分析。GO分析中显著富集的生物学过程结果显示lncRNA可能与细胞增殖、细胞周期和凋亡过程有关,图3。KEEG通路分析的结果表明差异表达的lncRNA可能通过各种代谢途径参与了垂体瘤的侵袭过程,见图4。
3. 讨论
近年来,越来越多的研究表明,lncRNA在包括癌症在内的许多疾病中发挥着重要作用,已受到越来越多的关注。随着高通量测序技术的发展,大量的lncRNA被发现,其作用机制的研究也取得了一定的进展。然而,lncRNA在侵袭性垂体腺瘤中的全基因组表达谱尚未见报道,笔者的研究填补了这一空白。
在本研究中,笔者使用RNA-seq来鉴定lncRNA和mRNA在4例侵袭性垂体腺瘤和3例非侵袭性垂体腺瘤中的表达谱。发现35个lncRNA的表达存在显著差异,其中21个表达上调,14个表达下调。另外,笔者的研究中发现229个mRNA表达上调,234个mRNA表达下调。因为lncRNA的大部分功能是未知的,特别是在侵袭性垂体腺瘤中,因此,笔者建立了共表达网络。基于Pearson相关系数筛选了34个lncRNA和274个mRNA,构建了一个lncRNA-mRNA网络。在lncRNA- mRNA网络中,通过Cytoscape插件计算得到了排名前五lncRNA。最后,对共表达网络中的mRNA进行了GO和KEGG分析,预测侵袭性垂体腺瘤可能的发生机制。结果显示差异基因主要与肿瘤的侵袭和代谢有关,如代谢通路、肿瘤内通路、内吞作用等。
在共表达网络中寻找到了FAM66A、TPTEP1、AC017116.11、RP11-1114A5.4和RP11-246K15.1共五个关键的lncRNA。它们在网络中的Degree排名位列前五,数值分别52、42、33、29和28。LncRNA TPTEP1已被报道通过影响IL-6/STAT3信号通路抑制肝癌细胞的进展[15]。由于对lncRNA的研究还处于早期阶段,大多数lncRNA的功能还不清楚。在本研究中FAM66A、AC017116.11、RP11-1114A5.4和RP11-246K15.1尚未见有报道,但是它们在共表达网络中处于核心位置,具有作为垂体瘤分子标志物潜在的研究价值。
本研究存在一定的局限性。首先,样本数量有限,笔者只用了4个侵袭性垂体腺瘤样本和3个非侵袭性垂体腺瘤样本来获取转录组数据。其次,笔者没有使用功能实验来验证本假设。但这项研究也有明显的优势,笔者的研究是第一次分析侵袭性垂体腺瘤完整的lncRNAs表达谱。
总之,本研究首次提供了lncRNA在侵袭性和非侵袭性垂体组织中的表达。生物信息学分析丰富了与病理机制相关的信号通路。通过lncRNA-RNA共表达网络,筛选出5种可能与侵袭性垂体腺瘤病理机制密切相关的关键lncRNA。笔者的发现将有助于了解侵袭性垂体腺瘤的发病机制,并可为侵袭性垂体腺瘤提供早期诊断生物标志物和潜在的治疗靶点。
-
表 1 观察组与对照组出疹持续时间及发热持续时间对比( $ \bar x \pm s $)
Table 1. Comparison of eruption duration and fever duration between observation group and control group ( $ \bar x \pm s $)
分组 出疹持续时间(d) 发热持续时间(d) 对照组 3.0 ± 0.6 3.1 ± 0.7 观察组 6.00 ± 1.7 6.000 ± 1.7 t 9.439 -9.240 P 0.000* 0.000* *P < 0.05。 表 2 观察组与对照组肝功能、肾功能、心肌酶对比( $ \bar x \pm s $)
Table 2. Comparison of liver function,renal function and myocardial enzymes between the observation group and the control group ( $ \bar x \pm s $)
分组 肝功能 肾功能 心肌酶 GGT(U/L) AST(U/L) ALT(U/L) CREA(μmol/L) LDH
(U/L)CK
(U/L)观察组 100.3 ± 158.91 84.11 ± 135.97 80.171 ± 130.11 80.80 ± 146.89 573.26 ± 952.70 128.20 ± 148.34 对照组 36.79 ± 31.79 29.76 ± 14.78 24.82 ± 3.14 57.68 ± 18.11 232.85 ± 68.30 74.38 ± 32.83 t 2.286 2.317 2.456 2.078 2.066 P 0.025* 0.024* 0.017* 0.366 0.042* 0.043* *P < 0.05。 表 3 观察组与对照组的IL-6、IL-8、CD4+T淋巴细胞数对比( $ \bar x \pm s $)
Table 3. IL-6,IL-8 and CD4+T lymphocytes counts was compared between the observation group and the control group ( $ \bar x \pm s $)
分组 IL-6(pg/mL) IL-8(pg/mL) CD4+T淋巴细胞数(个/μL) 对照组 15.24 ± 13.49 31.56 ± 20.69 786.75 ± 359.14 观察组 88.66 ± 213.02 93.88 ± 100.58 445.66 ± 354.26 t 2.005 3.540 −3.971 P 0.049* 0.001* 0.000* *P < 0.05。 -
[1] 王宇明, 李梦东. 实用传染病学[M]. 北京: 人民卫生出版社, 2016: 483-484. [2] World Health Organization. Varicella and herpes zoster vaccines: WHO position paper, June 2014[J]. Wkly Epidemiol Rec, 2014, 89 (25): 265–87. [3] Kennedy P G E,Mogensen T H. Determinants of neurological syndromes caused by Varicella-Zoster virus (VZV)[J]. J NeuroVirol,2020,26(4):482-495. doi: 10.1007/s13365-020-00857-w [4] 王春又,王欢,宋志强,等. 成人重症水痘合并多器官衰竭二例并文献复习[J]. 中国麻风皮肤病杂志,2021,37(2):96-99. [5] 李兰娟, 任红. 传染病学[M]. 9 版. 北京: 人民卫生出版社, 2018: 85-88. [6] 杨栋,甘云,伍晓梅,等. 重症水痘合并弥散性血管内凝血的临床分析[J]. 临床医学研究与实践,2020,5(22):84-85,88. doi: 10.19347/j.cnki.2096-1413.202022031 [7] 韩洋,费文敏,李可可,等. 免疫抑制剂的概念和类型[J]. 皮肤科学通报,2021,38(4):291-297. [8] 龚非力. 医学免疫学[M]. 北京: 科学出版社, 2014: 47-48, 117. [9] Cakir F B,Timur C,Yoruk A,et al. Seroconversion status after single dose and double doses of varicella vaccination in children with leukemia[J]. Pediatr Hematol Oncol,2012,29(2):191-194. doi: 10.3109/08880018.2011.593093 [10] 秦新华,王芳,赵保玲,等. 丙种球蛋白用于小儿重症出血性水痘治疗中的效果分析[J]. 医药论坛杂志,2021,42(2):119-122. [11] 陈莹莹. 益肾活血法干预小儿原发性肾病综合征对 NGAL、 KIM-1、IL-18 的影响[D]. 福州: 福建中医药大学硕士学位论文, 2017. [12] 李珍,梅文龙,柴艳婷,等. IL-6 基因 rs1800796 多态性与小 儿热性惊厥易感性的 Meta 分析[J]. 重庆医学,2017,46(15):2089-2093. doi: 10.3969/j.issn.1671-8348.2017.15.025 [13] 李由. 对普通水痘患儿与突破性水痘患儿T淋巴细胞亚群指标的对比分析[J].当代医药论丛, 2016, 14(22): 90-91. [14] 杨美蓉,蒋亦明. 水痘患者细胞免疫功能的检测及意义[J]. 内蒙古中医药,2008,27(11):107-108. doi: 10.3969/j.issn.1006-0979.2008.11.083 [15] 陈娓,刘军连. 带状疱疹发病机制的研究进展[J]. 中国医学文摘:皮肤科学,2017,34(1):33-38. [16] Monneret D,Mestari F,Djiavoudine S,et al. Wide-range CRP versus high-sensitivity CRP on Roche analyzers:Focus on lowgrade inflammation ranges and high-sensitivity cardiac troponin T levels[J]. Scand J Clin Lab Invest,2018,78(5):346-351. doi: 10.1080/00365513.2018.1471618 期刊类型引用(9)
1. 李早慧,都晓伟. 滇重楼药材的研究进展及其质量标志物(Q-Marker)预测分析. 中草药. 2023(09): 3032-3048 . 百度学术
2. 王晓菲,李鹏,杨玥,徐旭,林娟,杨洪军,李楚源. 重楼皂苷抗炎、抗氧化及抗肿瘤作用机制研究进展. 药物评价研究. 2023(12): 2699-2708 . 百度学术
3. 刘平安,张国民. 七叶一枝花药理机制研究进展. 湖南中医杂志. 2022(06): 202-206 . 百度学术
4. 张毅,郭亮,谢辉,田源,苏文莉. 中药巴布剂研究进展. 实用中医药杂志. 2022(07): 1269-1271 . 百度学术
5. 王宇飞,江媛,杨成金,王婧,徐志超,刘颖琳,段宝忠. 滇重楼化学成分、药理作用和临床应用研究进展. 中草药. 2022(23): 7633-7648 . 百度学术
6. Guoqing ZHAO,Lamei ZHOU. Research Progress of External Treatment of Acute Gouty Arthritis with Traditional Chinese Medicine. Medicinal Plant. 2021(03): 7-9 . 必应学术
7. 李梅,王燕,邝丽娟. 中医外治法干预急性痛风性关节炎白细胞介素的研究进展. 风湿病与关节炎. 2020(02): 73-77 . 百度学术
8. 杨建宇,刘冠军,刘白云,刘桂香,张举昌,仪忠宝,祝之友,唐祖宣. 中华中医药道地药材系列汇讲(8)道地药材滇重楼的研究近况. 现代医学与健康研究电子杂志. 2020(08): 117-120 . 百度学术
9. 周娟,修光辉,熊伟,凌斌,孙洁,陈献忠. 重楼治疗脓毒症研究进展. 中医学报. 2019(08): 1625-1629 . 百度学术
其他类型引用(7)
-