《Inflammation Research》：Temporal transcriptomic profiling identifies core regulators in cytokine-induced gut barrier disruption in differentiated Caco-2 cells
编辑推荐：
摘要：肠上皮屏障维持肠道稳态，其破坏可导致炎症性肠病（IBD）等炎症性疾病。虽然屏障修复是一种有前景的治疗策略，但目前大多数方法靶向免疫反应而非上皮功能。由于屏障损伤与修复经历不同阶段，本研究旨在明确时间分辨的上皮转录程序，并鉴定可用于屏障导向治疗干预的候选调
摘要：肠上皮屏障维持肠道稳态，其破坏可导致炎症性肠病（IBD）等炎症性疾病。虽然屏障修复是一种有前景的治疗策略，但目前大多数方法靶向免疫反应而非上皮功能。由于屏障损伤与修复经历不同阶段，本研究旨在明确时间分辨的上皮转录程序，并鉴定可用于屏障导向治疗干预的候选调控节点。研究人员使用分化的Caco-2上皮单层模型细胞因子诱导的屏障破坏，并通过时间分辨bulk RNA-seq（批量RNA测序）进行分析。利用成对差异表达分析、似然比检验（LRT）、通路富集和上游调控因子推断来鉴定转录动力学和候选调控因子。通过药物扰动测试候选节点，并通过qPCR（定量聚合酶链反应）和组织学在人IBD结肠组织中评估疾病相关性。时间分辨分析揭示涉及Rho GTP酶循环、sirtuin信号、非经典WNT信号以及连接/EMT（上皮-间质转化）相关通路的动态上皮程序。上游调控因子推断突出JAK2、TNF、gasdermin D和β-estradiol，这些药物扰动减轻了细胞因子诱导的屏障丧失。HNF4A（肝细胞核因子4α）作为调节包括SATB2和SLC26A3等靶点的肠上皮转录枢纽出现。值得注意的是，SLC26A3在IBD结肠组织中显著降低。研究人员得出结论，该研究提供了一个时间分辨的上皮框架，并强调了屏障导向干预的治疗机会。
论文解读：时间分辨转录组解析细胞因子诱导肠上皮屏障破坏的核心调控机制
研究背景方面，肠上皮屏障是维持肠道与全身稳态的核心结构，可选择性允许营养素吸收并限制腔面抗原、病原体及有害大分子通过。
其破坏参与多种病理过程，尤以炎症性肠病（IBD，包括溃疡性结肠炎UC和克罗恩病CD）最为典型，这类慢性复发?缓解性疾病以免疫失调和肠上皮屏障完整性受损为特征。
当前IBD治疗主要以免疫抑制为主，如靶向细胞因子信号通路和白细胞招募，Janus激酶（JAK）抑制剂Upadacitinib等在UC中显示疗效，但专门针对肠上皮屏障功能恢复的药物仍属空白。
已有证据表明，肠屏障高通透性可在CD患者健康一级亲属中出现，并预示CD复发风险，说明屏障缺陷可能是疾病发生发展的早期因素。
上皮屏障完整性主要由顶端紧密连接（TJ）复合体调控，关键跨膜蛋白如claudins、occludin及胞质支架蛋白如zonula occludens?1（ZO?1）构成动态信号枢纽；炎症细胞因子如TNF?α、IL?13可上调孔隙形成蛋白claudin?2而增加旁细胞阳离子通透性，促炎信号还通过myosin light chain kinase（MLCK）磷酸化myosin light chain（MLC）致周连接肌动蛋白环收缩而破坏TJ；缺氧诱导因子（HIF）通路经由脯氨酰羟化酶域蛋白（PHD）抑制稳定HIF1α可增强屏障保护基因表达；Rho GTP酶家族控制TJ组装与细胞骨架动态，其失调与肠通透性及黏膜炎症相关。
尽管分化Caco?2人结肠腺癌细胞单层因表达TJ、刷状缘酶、离子转运体等而广泛用于屏障机制与化合物筛选，但既往基因组尺度分析多局限于单时间点或微阵列共培养，缺乏上皮单独的、多时间点RNA?seq时间进程数据。
为此，研究人员通过建立细胞因子（TNF?α＋IFN?γ）诱导的分化Caco?2单层屏障破坏模型，进行时间分辨bulk RNA?seq（批量RNA测序），整合成对差异表达与似然比检验（LRT）解析动态转录程序，利用Ingenuity Pathway Analysis（IPA， ingenuity通路分析）做通路富集与上游调控因子推断，进而通过药物扰动、qPCR（定量聚合酶链反应）、3D微流控OrganoPlate模型及人IBD组织cDNA阵列与组织染色验证，旨在绘制上皮内在的时序分子架构，鉴定可干预的核心调控节点，为屏障导向治疗提供时间与机制依据。论文发表于《Inflammation Research》。
主要关键技术方法如下：研究人员以Cas9表达的分化Caco?2细胞在Transwell聚酯膜上培养12天形成极性单层，施加TNF?α＋IFN?γ双细胞因子处理，设0?h（对照）及6、9、18、30、54?h时间点，每组生物学重复4～8例；通过transepithelial electrical resistance（TEER，跨上皮电阻）与FITC?葡聚糖通透性评估屏障功能，CellTiter?Glo检测活力；bulk RNA?seq（批量RNA测序）建库用TruSeq Stranded mRNA试剂盒、NovaSeqX测序，比对Human.B38基因组、GENCODE v33注释得基因计数；成对差异表达用DESeq2（FDR＜0.05、|FC|＞1.5），低表达（FPKM＜1）非蛋白编码基因剔除，union DEGs做层次聚类；时间动态用DESeq2似然比检验（LRT，full model含Time＋Treatment＋Treatment:Time，reduced model去交互项，adj?p＜1e?10），rlog变换后DEGreport degPatterns聚类；通路与上游调控因子分析用IPA；候选上游节点用相应抑制剂/抗体/化合物共处理做TEER与通透性验证，并在3D微流控OrganoPlate（OrganoReady Colon Caco?2）中复测；人样本方面采用商用Crohn’s/Colitis cDNA阵列（非IBD?13例、CD?35例、UC?43例）qPCR检测SLC26A3，以及正常结肠9例、CD?6例、UC?9例石蜡标本做SLC26A3免疫组化定量；统计用ANOVA＋Sidak、Mann?Whitney?U、Kruskal?Wallis＋Dunn等，GraphPad Prism分析。
结果部分，研究人员首先建立模型并验证表型。通过细胞因子诱导肠屏障破坏模型表征与转录组概览：研究人员用TNF?α＋IFN?γ处理分化Caco?2单层，发现54?h内TEER显著降低约25%，FITC?葡聚糖通透性增高两倍以上，细胞活力未受影响；qPCR显示IFN?刺激基因MX1显著上调，TNF?诱导基因IL?8持续升高，屏障相关基因与细胞黏附基因如occludin、claudin?1等多数在54?h下调，分化标志物ALPI、SI、VIL1、DPP4在对照中稳定、细胞因子下显著下调；bulk RNA?seq的PCA显示PC1（80%方差）对应细胞因子主效应，PC2（11%）捕获早（6～9?h）晚（18～54?h）时间趋势，说明存在明确时序转录信号。
成对差异表达分析与通路富集：研究人员对每时间点细胞因子组vs对照做DESeq2成对差异表达，共获得5352个DEGs（差异表达基因），层次聚类划分为8个簇（1～3持续下调，4早期下调，5～7持续上调，8早期上调）；IPA富集显示簇1～3共同核心通路为I相化合物功能化（药物/异源代谢）、CLEAR（溶酶体协调表达与调控）信号、鞘脂代谢、离子通道转运；簇5～7共同核心含40条通路，选出Rho GTP酶循环与sirtuin信号网络做qPCR验证；Rho GTP酶循环相关基因（如NCF2、NOX1、IQGAP3、BAIAP2L2、RND1、JAG1等）及sirtuin信号相关基因（SLC2A1、NAMPT、NOS2、MYC、SOD2、DUSP6等）在细胞因子处理后显著上调，与RNA?seq趋势一致，说明Rho GTP酶循环介导细胞骨架重塑与ROS产生、sirtuin信号参与代谢?氧化?炎症调控是屏障破坏中动态程序。
时间进程LRT分析与通路：研究人员用似然比检验（LRT）捕捉Time×Treatment交互的时序基因（adj?p＜1e?10），18初始簇合并为10主簇（A～J，I、J因DEG少弃用），其中A～C为全程上调，D为早期瞬时上调，E～G为全程下调，H为早期瞬时下调；IPA显示A～C重叠核心通路为Rho GTP酶循环与中性粒细胞脱颗粒，E～G重叠核心为突触生成信号；鉴于紧密连接信号在B、C共现、非经典WNT信号在A、B共现且与屏障相关，研究人员将紧密连接、细胞连接组织、EMT相关归为上皮屏障重塑网络，非经典WNT（WNT/Ca2+、PCP）选代表基因验证；qPCR证实CLDN1、F11R（JAM?A）、LAMC2、TNFRSF11B持续上调（可能代偿性连接重塑/修复），MYO10、CTNNB1（β?catenin）呈双相（早期上调后期下调反映黏附协调崩溃）；非经典WNT相关PSME1、PIK3R3、WNT4、RAC1、SMURF1上调（STK31变化不显著），说明非经典WNT与细胞骨架?极性?可塑性调控被动态激活。
上游调控因子推断与功能验证：研究人员对成对DE簇和LRT簇分别做IPA上游调控因子分析，筛选出显著分子类别（转录调控因子、细胞因子、化学药物等），关键上游含JAK2、TNF、gasdermin D（GSDMD）、β?estradiol、HNF4A等；功能验证表明，JAK1/2抑制剂tofacitinib、TNF中和抗体adalimumab、gasdermin D抑制剂disulfiram共处理能显著减弱TEER下降与通透性增加；β?estradiol也增强屏障抗性；在3D微流控OrganoPlate中，tofacitinib（100?μM）、disulfiram（10?μM、1?μM）同样恢复TEER至近对照水平；HNF4A激动剂N?反式?caffe酰酪胺（NCT）在模型中未逆转屏障破坏且增加通透性，可能与炎症强度高、HNF4A本身下调、SATB2等染色质环境改变及Caco?2转化背景有关。
HNF4A?SATB2?SLC26A3轴与IBD组织关联：研究人员发现HNF4A、SATB2、SLC26A3在细胞因子54?h均显著下调，SLC26A3自6?h即持续抑制；商用IBD cDNA阵列显示SLC26A3 mRNA在CD（35例）和UC（43例）病变结肠显著低于非IBD（13例）；组织染色定量证实SLC26A3蛋白面积占比在CD、UC结肠上皮明显低于正常结肠（9例），与转录趋势一致，支持HNF4A为上皮转录枢纽、其下游SLC26A3（氯/碳酸氢根交换体，维持离子/液体稳态与黏液扩张）在炎症与IBD中一致下调。
讨论部分总结：研究人员指出TNF?α＋IFN?γ协同破坏上皮屏障已有报道（IFN?γ先上调TNFR2而后TNF?α触发 dysfunction），本模型中54?h内屏障受损但无细胞毒性，适合解析上皮内在程序。时间分辨转录谱区分早（6～9?h）晚（18～54?h）响应：早期以快速信号与重塑（上皮屏障重塑网络、非经典WNT、部分Rho GTP酶、sirtuin相关早期成分）为主，晚期伴持续屏障基因抑制、稳态基因下调及某些持续炎症?代谢程序（Rho GTP酶循环整体持续上调含NCF2、NOX1促ROS、IQGAP3促肌动蛋白聚合、RND1弱化细胞黏附、JAG1提示Notch可能代偿；sirtuin通路中SLC2A1上调反映SIRT1/6失抑制、NAMPT强诱导但SIRT6依赖抑制可能有保护、NOS2促硝基应激、MYC促增殖/炎症、SOD2上调或为抗氧化代偿）。Rho GTP酶上下文依赖：RhoA维护屏障但过度ROCK有害，本研究未直接阻ROCK，文献提示ROCK抑制剂可减轻类似破坏。Sirtuin结果解释为炎症中sirtuin轴基因改变可能兼顾保护性代谢?抗氧化代偿与选择性炎症贡献。上皮屏障重塑网络内CLDN1、JAM?A、LAMC2、TNFRSF11B持续上调或为修复尝试，MYO10、CTNNB1双相反映早期骨架适配到后期失协调；非经典WNT（WNT4、RAC1、PIK3R3、SMURF1、PSME1）动态活化指向细胞极性?可塑性与蛋白水解调控参与。上游推断验证证明JAK2、TNF、GSDMD（disulfiram阻断非致死质膜孔而非仅细胞死亡）、β?estradiol可减轻屏障丧失，其中JAK抑制剂、TNF抗体已是临床IBD药，在纯上皮系统亦有效，说明上皮自身这些节点足以介导部分屏障破坏，模型可优先挖掘上皮可成药点；GSDMD非经典质膜穿孔提示非裂解膜损伤参与屏障缺陷。HNF4A为关键上皮转录枢纽，其靶基因SATB2（染色质组织者限制HNF4A结合）、SLC26A3在细胞因子模型与IBD组织一致下调，但HNF4A激动剂NCT未能挽救，因本模型协同炎症强、HNF4A蛋白低、可能需SATB2等协同、且Caco?2转化背景与原代上皮有别。研究局限含Caco?2肿瘤来源缺正常上皮全复杂性、单剂量细胞因子不模拟体内动态波动（TNF?α 100?ng/mL偏高但常用、培养基未更换致后期浓度可能降但仍维持表型）、时间起点6?h错过即刻早期（＜2?h）转录爆发、RNA?seq某些时间点重复偏少侧重通路而非单基因严格显著、未做HNF4A?SLC26A3因果功能验证（需后续增益/缺失）、上皮单独缺免疫?基质?微生物互作。方法上学成对DE加LRT互补：两者共凸现Rho GTP酶循环为核心，成对DE更抓各时间点特异通路（sirtuin等），LRT抓交互驱动时序模式（非经典WNT、上皮屏障重塑等）；上游推断结合药物扰动把转录信号转为功能保护证据。结论为：本研究提供了分化Caco?2细胞因子诱导屏障破坏的时间分辨上皮转录框架，确定Rho GTP酶循环、sirtuin信号、非经典WNT、上皮屏障重塑/EMT等为动态核心通路，鉴定上游可干预节点JAK2、TNF、GSDMD、β?estradiol及转录枢纽HNF4A?SATB2?SLC26A3轴（SLC26A3在IBD人类组织一致下调），该上皮优先、时间意识模型可绘制屏障失效时序架构、提名阶段特异性可成药点，为屏障导向治疗提供机制假设与优先化平台。
