编辑推荐：
摘要
致病性立克次体（Rickettsia）为专性胞内细菌，必须在哺乳动物宿主或节肢动物载体细胞内存活。尽管这些细菌在感染过程中在不同胞内环境间转换，但其编码的假定转录因子极少，且其基因调控网络大多未知。由于其与真核细胞存在不可分割的联系，病原体的转录谱分析因
摘要
致病性立克次体（Rickettsia）为专性胞内细菌，必须在哺乳动物宿主或节肢动物载体细胞内存活。尽管这些细菌在感染过程中在不同胞内环境间转换，但其编码的假定转录因子极少，且其基因调控网络大多未知。由于其与真核细胞存在不可分割的联系，病原体的转录谱分析因宿主RNA的大量污染而复杂化，尤其在细菌载量本就很低的感染条件或阶段。本研究采用病原体杂交捕获技术（PatH-Cap），通过富集细菌转录本同时去除宿主和rRNA分子，以改进文库构建。利用PatH-Cap，研究人员探索了感染最初24小时内的转录变化，包括感染起始——这一感染阶段难以用标准文库构建方法进行分析。接着，研究人员基于基因的时间趋势进行聚类，揭示了在感染不同阶段表达上调或下调的基因群。研究还强调了在生长和致病性中具有已知作用的基因的多样时间表达趋势，包括翻译和细胞分裂基因、分泌效应蛋白及分泌系统组分。最后，研究人员鉴定出639个反义RNA分子，其中许多也显示出强烈的时间趋势。本工作表明，如PatH-Cap般灵敏的转录谱分析方法在剖析历史上构成重大技术挑战的胞内病原体感染驱动基因表达网络方面具有巨大前景。
一、研究背景与问题
斑点热群立克次体（Spotted fever group Rickettsia）是蜱传专性胞内细菌，其中许多可导致人类危及生命的血管疾病。由于与宿主和载体细胞存在不可分割的相互作用，立克次体在其复杂生命周期中演化出大量机制以操纵真核细胞并逃避免疫应答。然而，细菌如何在宿主体内感知环境信号并执行特定基因表达程序以转换这些生命周期阶段仍属未知。
立克次体经过还原进化，基因组极小（约130万碱基对），仅含约1200个基因，其中大量为功能未知的假定蛋白。此外，其编码的转录调控因子极少，这引发了关于基因表达程序如何在不同宿主物种间或整个感染周期中被介导的问题。转录谱分析可为理解新基因的表达模式提供关键见解，但胞内病原体的RNA测序因感染细胞总RNA提取物中存在大量宿主RNA而复杂化。此前研究通过多种技术组合提高了立克次体转录本的占比，但即使在多步富集过程后，映射到立克次体基因组的读数百分比通常仍较低，且高度依赖于初始感染量、感染持续时间和感染系统。因此，在生理相关感染条件下，尤其是在感染早期细菌载量低时，对动态基因表达变化进行系统水平分析一直存在技术缺口。
二、研究目的与意义
本研究旨在首次将病原体杂交捕获（PatH-Cap）方法应用于立克次体物种，以富集细菌转录本，实现对模型斑点热群立克次体——帕氏立克次体（Rickettsia parkeri）在哺乳动物宿主细胞内感染生命周期的时间性转录谱分析。通过克服宿主RNA污染的技术障碍，研究期望揭示立克次体在建立感染过程中的动态转录变化，鉴定具有时间特征的基因网络，预测操纵子结构，并识别反义RNA转录本，从而为理解立克次体的基因调控网络及其在感染过程中的作用提供重要基础。相关成果发表于《mSphere》期刊。
三、关键技术方法
研究采用以下关键技术：1. 设计并应用针对帕氏立克次体基因组的定制杂交捕获探针文库（包含10746个非重叠的120核苷酸探针），通过PatH-Cap技术从感染的人上皮细胞（A549）总RNA中富集细菌转录本，同时去除宿主RNA和rRNA。2. 在感染后多个时间点（0、1、4、8、12、16、24小时）收集RNA样本，构建链特异性RNA测序文库并进行测序。3. 利用生物信息学工具进行差异表达分析、加权基因共表达网络分析（WGCNA）聚类、基因本体（GO）富集分析、操纵子预测（使用Rockhopper软件）及反义RNA鉴定。
四、研究结果
1. 为帕氏立克次体建立PatH-Cap
研究人员成功设计了帕氏立克次体特异性杂交捕获文库。应用该技术后，感染后24小时样本中映射到细菌mRNA和小RNA（sRNA）的读段比例从捕获前的高水平宿主污染提升至88.77%，0小时样本（纯化细菌）达到96.60%。即使在每个宿主细胞低至一个细菌的感染复数下，PatH-Cap也能在感染早期（1、4小时）实现稳健的细菌转录本富集。
2. PatH-Cap实现精细时间进程转录谱分析
比较不同时间点的差异表达基因发现，24小时与48小时感染后仅有一个基因差异表达，而0小时与24小时比较则发现大量差异表达基因，表明帕氏立克次体在建立感染过程中经历了动态转录变化。主成分分析显示，1至24小时的样本按时间顺序聚类，证实了数据中存在时间特征。
3. WGCNA促进时间相关立克次体基因的聚类
通过加权基因共表达网络分析，基因被分为6个具有独特时间表达模式的簇。例如，簇2基因在感染早期（1-4小时）表达最高；簇3和簇4基因表达随时间下降；簇5和簇6基因表达随时间增加。簇1包含无明显共同时间特征的基因。利用CytoHubba分析潜在枢纽基因未发现注释的转录因子。仅在簇4基因上游100核苷酸区域发现一个富集基序“GSAATGACR”。
4. 基因表达与基因组位置无明显相关性
绘制基因在1小时和24小时的平均每百万转录本数（TPM）及簇归属图显示，高表达的核糖体和翻译相关基因在染色体上成簇分布，但并非靠近复制起点。各簇基因在染色体上呈分散分布，表明其时间网络不受基因组位置直接影响。
5. 簇包含具有不同功能分类的差异表达基因
对各簇进行基因本体富集分析发现，簇2（早期高表达）在跨膜转运相关基因上显著富集；簇6（表达随时间增加）在核糖体、翻译及RNA结合产物编码基因上显著富集；簇5（表达随时间逐渐增加）在脂质转运基因上富集程度最高。许多假定基因和假基因也显示出特定的时间表达模式。
6. 跨膜转运基因表达在建立感染期间达到峰值
簇2中跨膜转运相关基因的富集提示其在感染早期可能发挥重要作用。转座子诱变显示该簇及簇4中的某些转运蛋白基因突变会减弱帕氏立克次体在哺乳动物细胞感染中的噬斑大小，表明其潜在功能重要性。
7. 核糖体、翻译和脂质转运相关基因表达随感染进展而增加
簇6中核糖体和翻译相关基因的表达在≥8小时达到峰值，与细菌退出滞后期进入指数生长期的时间一致，以满足蛋白质合成增加的需求。簇5中脂质转运基因（如脂多糖转运蛋白编码基因）的表达逐渐增加，可能反映了复制期间细胞被膜组装和维持所需脂质转运基因的时间性转录控制。
8. 促进宿主相互作用的基因产物呈现多样的转录趋势
分析已知在感染特定阶段起作用的细菌基因发现，尽管肌动蛋白成尾相关蛋白RickA和Sca2在功能上存在时间差异，但其编码基因转录趋势相似，均属于簇6，表达在0-4小时增加后保持稳定。已鉴定的分泌效应蛋白编码基因表现出不同的时间模式：如Pat1、Risk1、Sca4属于簇6；而RLip和RARP2属于簇2，表达在1小时达峰后下降。一组新发现的分泌效应蛋白编码基因大多属于簇6，SrfF属于簇4。
9. T4SS操纵子分裂产生不同的时间特征
对编码IV型分泌系统（T4SS）的基因进行分析发现，与农杆菌（Agrobacterium）简单的T4SS操纵子不同，立克次体的T4SS基因在基因组上分散分布，并预测存在多个操纵子。表达分析显示，即使位于同一基因组区域或预测操纵子内的T4SS基因，也呈现出多样的时间表达趋势，反映了其复杂的转录调控。
10. 反义RNA显示多样的表达模式及其与对应mRNA的关系
研究人员鉴定出639个基因（占所有编码基因的42.77%）存在潜在的反义RNA转录本。这些基因在催化活性（如ATP结合和水解）相关的GO术语上显著富集。对反义RNA与对应mRNA表达水平的皮尔逊相关性分析显示，35.99%的基因对呈正相关，17.68%呈强负相关，其余无明显相关。以细胞分裂基因（zapE, ftsI）、DNA聚合酶I基因（polA）、DNA修复基因（recG）和赖氨酸甲基转移酶基因（pkmt2）为例，展示了其mRNA与反义RNA之间或正或负的强相关性。
五、讨论与结论
本研究通过应用PatH-Cap技术，首次在低细菌载量下实现了帕氏立克次体感染早期和晚期时间性转录谱的全面、稳健分析。研究揭示了在感染过程中具有特定时间表达模式的基因网络，改进了操纵子预测，并鉴定了大量反义RNA转录本。
研究结论可总结为：1. PatH-Cap能有效富集帕氏立克次体转录本，即使在高宿主RNA污染和低细菌载量下也能实现高读段深度，从而能够分析感染早期阶段的转录变化。2. 时间进程分析揭示，帕氏立克次体在感染最初24小时内经历动态的转录重编程，基因可聚类为具有不同时间表达特征的组别，这些组别在功能上富集于跨膜转运、翻译、核糖体生物发生、脂质转运等过程。3. 许多与宿主相互作用相关的基因（包括效应蛋白和T4SS组分）表现出多样的时间表达模式，与其在感染生命周期中的阶段特异性功能相符。4. 在帕氏立克次体中鉴定出大量反义RNA，其中许多与其对应mRNA的表达呈现强时间相关性，提示其可能参与转录后调控。
本研究建立的数据集为未来探索特定基因的转录谱或在遗传或环境扰动下的全局变化提供了资源。PatH-Cap方法可应用于更多样的感染条件，以进一步阐明立克次体应对不同环境和引发疾病的调控网络。
