IF31.7｜类器官结合TraDIS揭示志贺菌肠道定植全基因组图谱

用肠类器官+转座子插入测序+贝叶斯模型，首次系统绘制志贺菌感染人类肠上皮所需基因图谱。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

志贺菌（Shigella flexneri）是引起细菌性痢疾的主要人类适应性病原体，其分子致病机制研究长期受限于缺乏能模拟人类生理的小动物模型。此前已知志贺菌通过III型分泌系统（T3SS）及其效应蛋白入侵肠上皮细胞，但大多数研究基于肿瘤来源的细胞系，这些细胞系在形态、代谢、信号通路和细胞死亡等方面与天然上皮存在显著差异。近年来，来源于干细胞的肠道类器官（enteroids/colonoids）为宿主-病原互作研究提供了更接近生理状态的三维模型。然而，将全基因组功能筛选技术（如TraDIS）应用于类器官感染模型仍面临感染瓶颈、富集偏倚等挑战，导致志贺菌在生理排列的肠上皮中定植所需的完整基因集尚未被系统定义。

科学问题

卡在哪里： 此前志贺菌毒力研究主要依赖肿瘤细胞系，缺乏在生理相关肠上皮中的全基因组功能筛选。感染模型中的群体瓶颈和细菌富集步骤的选择性生长效应会扭曲突变体表型，导致假阳性或假阴性，限制了TraDIS等技术的应用。

本文要回答： 本文旨在建立可扩展的肠类器官感染模型，结合TraDIS和贝叶斯统计建模，系统绘制志贺菌感染人类肠上皮所需的全基因组图谱，并解析关键毒力因子的作用机制。

技术路线

作者主要用了：肠类器官感染、TraDIS、RNA-seq、贝叶斯ZINB模型、条形码竞争感染、蛋白质组。

作者首先建立基底朝外（BO）极性的人肠类器官悬浮培养感染模型，通过优化感染条件和富集温度，减少选择偏倚。随后构建包含约13万个突变体的志贺菌Tn5转座子文库，分成43个子文库分别感染类器官，对输入和输出群体进行TraDIS测序。针对感染瓶颈导致的随机突变体丢失，开发了基于零膨胀负二项分布（ZINB）的贝叶斯模型，利用假基因估计技术参数，从而准确提取基因水平的定植适应度。最后通过条形码竞争感染实验验证关键基因，并结合蛋白质组和翻译分析揭示调控机制。

图 1　Fig. 1。图1展示了类器官感染模型和瓶颈评估。a,b显示野生型志贺菌在BO肠类器官中的时间序列成像，3小时出现感染灶，9小时扩散明显。c显示野生型与ΔmxiD共感染6小时后CFU差异超过60倍。d,e通过模拟和实验估计感染瓶颈约2000个细菌。f-h显示在37°C富集时，约1200个基因的突变体表现出适应性差异，而30°C时仅约130个，结合RNA-seq数据量化了毒力基因表达的成本。i,j通过条形码竞争实验验证ΔvirF和ΔvirB突变体在37°C富集时显著富集。

核心亮点

亮点 1｜建立可扩展的类器官感染模型

作者开发了高通量庆大霉素保护试验，使用基底朝外（BO）极性的人肠类器官悬浮培养。感染野生型志贺菌后，3小时出现感染灶，9小时内细胞间扩散明显。通过混合野生型和ΔmxiD（非侵袭）菌株感染，6小时后回收约10万CFU野生型菌，而ΔmxiD减少60倍以上，提供了宽动态范围。通过模拟和实验估计感染瓶颈约为2000个细菌，指导了后续子文库设计。

图 2　Fig. 2。图2展示了主TraDIS筛选的设计和结果。a为优化后的实验流程示意图。b显示43个子文库的转座子插入位点分布，合并后覆盖整个染色体和毒力质粒。c显示各子文库的测序读数，子文库11因读数过低被剔除。d为ZINB贝叶斯模型示意图。e,f分别显示染色体和毒力质粒上基因的log2FC热图，鉴定出105个染色体基因和38个毒力质粒基因对定植有显著贡献。

亮点 2｜全基因组定植必需基因图谱

通过43个子文库的TraDIS筛选，覆盖86%的注释基因，鉴定出105个染色体基因和38个毒力质粒基因对肠类器官定植有贡献。其中超过70个染色体基因此前未与志贺菌上皮定植关联，包括应激调控因子RpoS、转运蛋白、tRNA修饰酶MnmE/G等。毒力质粒上的必需基因包括virF、virB、T3SS结构组分和icsA等，验证了筛选的饱和度和可靠性。

图 3　Fig. 3。图3展示了T3SS效应蛋白在不同上皮环境中的需求。a为TraDIS筛选中效应基因的log2FC热图。b-e显示野生型志贺菌在BO/AO肠类器官和结肠类器官中的共聚焦成像。f-i显示野生型与ΔmxiD在不同类器官中的CFU比较，AO极性下野生型定植效率降低。j-m通过条形码竞争感染显示IpgB1、IpgB2、IpaA和VirA是主要效应蛋白，且IpaA在顶端侵袭中尤为重要。

亮点 3｜T3SS效应蛋白的上下文依赖性

通过条形码竞争感染在BO/AO极性的肠类器官和结肠类器官中验证，发现IpgB1、IpgB2、IpaA和VirA是早期定植的主要效应蛋白。IpaA在顶端侵袭中尤为重要，其缺失在AO类器官中导致约5倍定植减弱，而在BO类器官中影响较小。IpgB1/B2双缺失在AO结肠类器官中定植减弱超过10倍，表明效应蛋白需求具有上皮极性和肠段特异性。

图 4　Fig. 4。图4展示了MnmE/G tRNA修饰酶调控毒力表达的机制。a为MnmE/G依赖的tRNA修饰示意图。b显示ΔmnmE和ΔmnmG突变体在条形码竞争感染中定植指数显著降低。c,d为蛋白质组火山图，显示ΔmnmE和ΔmnmG突变体中25个毒力蛋白下调。e,f显示毒力蛋白编码基因中AGA和GGA密码子使用频率显著高于染色体基因。g-j显示IpaA蛋白在ΔmnmE和ΔmnmG突变体中表达降低，而mRNA水平不变，互补实验可恢复。k-n显示互补mnmE/G可恢复突变体的定植缺陷。

亮点 4｜tRNA修饰酶全局调控毒力表达

TraDIS筛选发现mnmE和mnmG突变体定植能力降低6-10倍。蛋白质组分析显示，ΔmnmE和ΔmnmG突变体中154个蛋白显著下调，其中25个为毒力质粒编码的T3SS相关蛋白。这些毒力蛋白的编码基因富含AGA和GGA密码子，而MnmE/G负责修饰识别这些密码子的tRNA。通过诱导表达IpaA证明MnmE/G直接调控其翻译，互补实验可恢复定植缺陷。这揭示了一种基于密码子偏好和tRNA修饰的毒力表达全局调控机制。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 肠类器官感染、TraDIS、RNA-seq、贝叶斯ZINB模型、条形码竞争感染、蛋白质组
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事病原微生物功能基因组学、宿主-病原互作或类器官感染模型研究的科研人员。

带走一句

将类器官感染模型与TraDIS和贝叶斯建模结合，可克服感染瓶颈和富集偏倚，实现病原体-宿主互作的全基因组功能筛选。