《Journal of Eukaryotic Microbiology》：Different Transcriptomic Responses After Moderate and Severe Heat Shock in Trypanosoma cruzi Epimastigotes
编辑推荐：
克氏锥虫（Trypanosoma cruzi）在其生活史中会经历热应激。本研究对不同克氏锥虫株的上鞭毛体（epimastigotes）进行了热休克响应的转录组学分析。在克隆CL Brener和Y株中，37°C下观察到的转录组统计学显著变化极少甚至没有，而在40
克氏锥虫（Trypanosoma cruzi）在其生活史中会经历热应激。本研究对不同克氏锥虫株的上鞭毛体（epimastigotes）进行了热休克响应的转录组学分析。在克隆CL Brener和Y株中，37°C下观察到的转录组统计学显著变化极少甚至没有，而在40°C下则观察到差异表达基因，其中CL Brener中有1018个上调转录本和734个下调转录本，Y株中有1426个上调转录本和1223个下调转录本。上调转录本富集于分解代谢通路，而下调转录本富集于分解代谢和合成代谢通路。此外，多个热休克蛋白（heat shock protein，HSP）基因家族不同成员的表达模式存在差异。表面蛋白基因家族成员（如转唾液酸酶（trans-sialidases）、GP63、黏蛋白（mucins）和黏蛋白相关蛋白（mucin-associated surface proteins，MASP））的表达模式也有所变化，尽管这些变化的生理意义尚不明确。研究得出结论：在克氏锥虫上鞭毛体中，37°C下的转录组调控对于8°C的温度升高而言不如预期显著，而在40°C下，诱导和抑制的转录本均出现功能富集，表征了典型的严重热休克响应。
克氏锥虫热休克转录组响应研究解读
一、研究背景与意义
克氏锥虫（Trypanosoma cruzi）是一种具有鞭毛的原生动物寄生虫，属于动质体目（Kinetoplastida），是恰加斯病（Chagas disease）的病原体，全球有数百万人感染，目前尚无有效治疗药物或疫苗。该寄生虫具有复杂的生活史，在昆虫媒介和哺乳动物宿主之间传播，经历多个发育阶段的转化。上鞭毛体（epimastigotes）在昆虫体内分化形成具有感染性的锥鞭毛体（metacyclic trypomastigotes），进入哺乳动物宿主细胞后分化为无鞭毛体（amastigotes）并增殖，最终释放为血流锥鞭毛体（bloodstream trypomastigotes）。
在生活史过程中，寄生虫会经历多种应激，其中热应激是重要类型之一。热应激可诱导分子伴侣和蛋白酶的产生，而营养应激和渗透应激则分别涉及P小体/应激颗粒形成和离子/渗透物质外排等不同响应机制。在锥虫科（trypanosomatids）生物中，基因表达调控具有特殊机制：蛋白编码基因由RNA聚合酶II转录为多顺反子转录本，通过多聚腺苷酸化和反式剪接（trans-splicing）加工为成熟mRNA，基因调控主要发生在转录后水平，主要通过mRNA稳定性和翻译调控实现。
热休克响应通常由仅几度的温度升高触发，足以引起蛋白质去折叠、非特异性聚集和蛋白质稳态的普遍丧失。经典的热休克响应涉及分子伴侣的诱导、蛋白水解系统组分的表达、RNA和DNA修饰酶的调控、代谢酶的重组、基因表达调控蛋白的参与以及细胞骨架和膜功能相关蛋白的恢复。
克氏锥虫不同株之间存在广泛的基因型和表型异质性，对应激的响应也存在差异。因此，本研究选取克隆CL Brener（DTU VI）和Y株（DTU II）两个株系的上鞭毛体，在29°C（标准培养条件，反映无脊椎动物媒介的典型环境温度）、37°C（模拟中度热应激，接近哺乳动物体温）和40°C（严重热休克，高于生理范围，用于触发经典热休克蛋白响应）三个温度条件下进行转录组分析，以探究寄生虫对热休克的基本生理响应以及不同株系间响应的潜在差异。该研究发表在《Journal of Eukaryotic Microbiology》。
二、主要技术方法
研究人员采用RNA测序（RNA-seq）技术对克氏锥虫两个株系（克隆CL Brener和Y株）的上鞭毛体在29°C、37°C和40°C处理3小时后的转录组进行分析。实验设置三个生物学重复。RNA样品经poly(A)+富集后构建测序文库，在Ion Torrent个人基因组测序仪上完成测序。测序数据经质量控制后比对到克氏锥虫CL Brener Esmeraldo-like参考基因组（版本v46）。差异基因表达分析采用CLC Genomics Workbench的RNA-seq模块以及edgeR和EnhancedVolcano R包，以倍数变化>1.5且FDR校正p值<0.05为显著性阈值。基因本体（Gene Ontology，GO）富集和通路富集分析使用EuPathDB在线工具进行。
三、研究结果
3.1 测序与参考基因组比对
研究人员共获得44,451,412条reads，平均读长范围为62.7至115.3 bp，至少80%的reads成功比对到参考基因组。在CL Brener上鞭毛体中，检测到10,831个注释蛋白编码基因中的9207个（85%）在所有温度条件下均有转录本检出，另有72、56和110个基因分别仅在29°C、37°C和40°C条件下检出。在Y株中，8908个基因（82%）在所有温度下均有检出，82、16和262个基因分别仅在29°C、37°C和40°C条件下检出。
3.2 差异基因表达
差异基因表达分析显示，两个株系的热休克转录组响应模式相似。令人意外的是，在37°C条件下，CL Brener中未观察到统计学显著差异表达基因，Y株中也仅有极少数基因出现显著变化，尽管这一温度变化对应8°C的升温幅度。相比之下，在40°C条件下，两个株系均观察到大量差异表达基因：CL Brener中鉴定出1752个差异表达基因（1018个上调，734个下调），Y株中鉴定出2649个差异表达基因（1426个上调，1223个下调）。
3.3 GO和通路富集分析
由于两个株系表现出高度相似的模式，研究人员仅展示了CL Brener的结果。40°C下上调的677个基因主要富集于有机物和含氮化合物代谢、大分子代谢和蛋白质代谢等GO生物学过程，对应类固醇降解、血红蛋白降解、乙醛酸降解、膦酸降解以及含硫化合物代谢（如半胱氨酸生物合成）等通路。下调的458个基因除富集于相同的GO生物学过程外，还富集于多种大分子和含氮化合物的生物合成过程，对应氨基酸降解、乙酰辅酶A生物合成、糖酵解、糖异生、磷酸戊糖途径、核苷酸代谢、三羧酸循环、脂肪酸降解以及胆固醇和其他类固醇降解等通路。此外，多个脂质代谢相关通路基因也受到调控，提示膜脂的稳态黏度适应（homeoviscous adaptation）机制。
3.4 热休克蛋白基因表达
热休克蛋白（HSP）基因家族不同成员在热休克后的表达模式存在差异。HSP110、HSP20和HSP10基因的转录本在40°C下大多上调；部分HSP60家族成员上调，这与HSP10和HSP60蛋白形成单一功能单元一致。相比之下，ClpB/HSP104和HSP33基因以及部分HSP70基因的转录本水平保持不变，这一表达模式与正常条件下存在的组成型HSP（cognate HSP）功能一致。两个株系的表达模式基本相似，但HSP60和HSP10家族成员存在株系间差异。值得注意的是，部分HSP40、HSP60、HSP70和HSP85/90基因家族成员的转录本在40°C下下调幅度可达4倍，其中下调成员的比例从不足20%到超过60%不等。HSP40家族作为锥虫科中最大的HSP基因家族，大多数成员转录本水平保持不变，约四分之一的成员分别上调或下调。HSP70家族成员在热休克后表现为不变或下调，未观察到上调，这与既往研究一致。两个HSP70基因（TcCLB.511211.170_1和TcCLB.511211.160_1）的转录本丰度比其他家族成员高10至20倍以上，线粒体HSP70基因的转录本丰度则高至少5倍。
3.5 表面蛋白基因表达
表面蛋白基因家族成员在40°C热休克后的表达模式同样存在差异。转唾液酸酶（trans-sialidases）家族526个检出成员中约70%上调、25%不变、5%下调。GP63家族172个检出成员中，CL Brener约65%上调、15%不变、20%下调，而Y株仅约50%上调。分散基因家族1（dispersed gene family 1，DGF-1）174个检出成员中，CL Brener约60%上调、30%不变、10%下调，而Y株仅约30%上调。黏蛋白（mucins）家族328个检出成员中约30%上调、60%不变、10%下调，尤其在CL Brener中大多保持不变。黏蛋白相关表面蛋白（MASP）家族531个检出成员中约70%上调、20%不变、10%下调。逆转录转座子热点（retrotransposon hot spot，RHS）蛋白家族233个检出成员中约5%上调、45%不变、50%下调。
四、讨论与结论
研究人员指出，37°C下基因表达缺乏统计学显著变化与既往观点一致，即寄生虫已为感染哺乳动物宿主做好准备，克氏锥虫及其他锥虫科生物在非热休克条件下存在HSP可能参与分化过程。研究结果提示这种准备状态可能延伸至整个转录组乃至蛋白质组。相比之下，40°C下观察到的差异表达基因模式与严重热休克的经典响应一致，表现为代谢通路和细胞过程的重组。GO和通路富集分析揭示的变化可能反映了高温下能量供应的重组，与在其他真核生物中的描述一致。脂质代谢相关基因的调控可能构成膜脂的稳态黏度适应。蛋白质稳态通过分子伴侣和蛋白酶（包括HSP）的诱导得以维持，多个组成型HSP在生理温度下即存在以协助蛋白质折叠。部分HSP基因家族（特别是HSP40和HSP85）成员表现出明显不同的转录本变化模式，研究人员推测不同HSP40成员可能赋予HSP70底物特异性，寄生虫HSP40基因在热休克期间的差异调控可能反映相应底物的需求。部分HSP基因家族成员在40°C下转录本下调的现象在乳腺癌细胞中也有描述，被视为一种失调，但其在克氏锥虫中的生理意义尚不明确。表面蛋白基因家族成员表达模式的差异可能对应不同的功能亚类（如细胞黏附与细胞保护），但需进一步研究加以阐明。本研究的局限性在于未评估热休克期间的蛋白质组，尽管蛋白质水平变化通常大致跟随mRNA水平变化，但翻译调控可能显著影响某些基因的表达水平。
研究结论：本研究表明，克隆CL Brener和Y株上鞭毛体在37°C热休克后转录组基本保持不变，提示寄生虫已为感染哺乳动物宿主做好准备。相比之下，两个株系在40°C下表现出典型的严重热休克转录响应。GO和通路富集分析显示40°C下多条通路受到调控，可能反映了寄生虫能量供应的重组以及膜脂的稳态黏度适应。此外，HSP和表面蛋白基因家族成员在40°C下表现出上调与下调转录本的不同模式，其中若干变化的生理意义尚不清楚。
