接下来,作者收集了患者的临床信息,以探讨肿瘤内微生物丰度与临床因素之间的关系。RDA/CAA环境因素分析显示,癌症抗原CA199对肿瘤组织细菌的影响最大,其次是白蛋白和白细胞(图4a),这些临床因素中,CA199和白蛋白是正相关的,而CA199与白细胞和高铁血红蛋白是负相关的。瘤内微生物分类和临床因素之间的相关性计算表明,Bacillus anthracis和P. azotoformans的数量与CA199呈正相关(图4b)。相比之下,Burkholderia tuberum和P. fungorum的数量与CA199呈负相关,并与ICC癌旁组织中该细菌的高丰度有关,由此作者推测这两种细菌可能具有抗肿瘤特性。
图4 环境因素分析显示细菌丰度与临床数据的相关性
5.功能预测表明细菌微生物群可能与氨基酸的生物合成有关
基于16S rRNA测序数据,结合eggNOG和KEGG数据库,作者预测了这些细菌的生物功能。结合eggNOG数据库,预测细菌的功能主要涉及代谢,特别是氨基酸代谢(图5a)。类似的,16S rRNA测序数据结合level 2的KEGG功能预测表明,细菌菌落功能主要与碳水化合物的氨基酸和核苷酸代谢途径有关(图5b)。结合level 3的KEGG数据信息,可以进一步确认细菌与氨基酸的生物合成、嘌呤代谢、嘧啶代谢和丙酮酸代谢的关联(图5c)。
图5 细菌微生物群的功能预测分析
6.差异分析显示癌旁组织中Burkholderia的丰度更高
考虑到癌和癌旁组织之间的差异,作者计算了它们不同菌落在物种水平上的特征。在门级水平上,线性判别分析(LDA)结果显示Proteobacteria在癌旁组织中的丰度较高,而Armatimonadetes、Verrucomicrobia和Fusobacteria在癌组织中的丰度较高(图6a)。使用具有高特异性的引物序列,通过RT-qPCR进一步验证了肿瘤组织中Verrucomicrobia DNA是异常丰富的(图6b)。通过qPCR进一步验证了肿瘤组织中Proteobacteria DNA的耗竭(图6d)。在科的水平上,大多数差异性生物在肿瘤组织中表现出较高的水平,包括Comamonadaceae、Staphylococcaceae、Verrucomicrobiaceae、Thiobacteriaceae和Peptococcaceae。在属的水平上,最丰富的细菌在肿瘤组织中表现出高水平,包括Acidovorax、Staphylococcus、Bdellovibrio、Roseimicrobium和Roseburia。在物种水平上,Ralstonia pickettii和Acinetobacter johnsonii的数量最多(图6c),然而,这两个物种有可能是环境中的污染性细菌。P. fungorum和P. azotoformans在癌和癌旁组织中有显著差异(图6c),二者在癌旁组织中的含量比在癌组织中的要高。作者由此猜想,癌旁组织中含量较高的细菌对肿瘤的发展有抑制作用,这与P. fungorum的荧光半定量分析结果相似(图6e)。
图6 癌和癌旁组织之间细菌含量的差异分析
7.P. fungorum对胆管癌细胞迁移和增殖的抑制作用
考虑到细菌液的上清液含有大分子和细菌的代谢物,作者以MTT(噻唑蓝)实验、Transwell实验和划痕实验来研究P. fungorum代谢物对癌细胞迁移的影响。Transwell实验的结果表明,加入P. fungorum上清液后,电镜下单位面积的细胞数明显减少,这代表细胞的增殖和迁移能力减弱(图7a)。划痕实验中,加入P. fungorum上清液后,QBC939细胞划痕愈合的百分比明显低于阴性对照样品(图7b),且划痕愈合率随着细菌上清液浓度的逐渐增加而降低。这些结果表明,细菌上清液抑制了不同类型的胆管癌细胞的迁移,这预示着细菌代谢物中可能存在抑制癌细胞迁移的因素。
接下来,作者在体内实验中研究了细菌对肿瘤表型的影响以及它们可能的功能。作者利用QBC939细胞系,在6周大的裸鼠中构建了一个皮下移植的肿瘤模型(图7c)。除了常规地将肿瘤细胞皮下移植到小鼠体内作为对照组外,还设计了另外三个实验组:肿瘤细胞加活菌、肿瘤细胞加细菌上清液、肿瘤细胞加活菌后在小鼠饮用水中加入抗生素(青霉素和链霉素的混合物)以消减小鼠体内的细菌。体内实验的结果显示,实验组小鼠的肿瘤体积相对于对照组有所减少(图7d)。与对照组相比,实验组中使用活菌的小鼠的肿瘤体积是最低的。在实验组中,从接种到3周,细菌上清液组和抗生素处理组在其体积变化上没有明显差异,不过抗生素组的最终平均肿瘤体积较大。作者还对经活菌处理的实验组小鼠的肿瘤切片进行了不同探针下的FISH实验,验证了小鼠肿瘤中P. fungorum的存在,这也在一定程度上说明了体内实验的准确性(图7e)。
图7 体外和体内试验显示P. fungorum对肿瘤生长的抑制作用
8.移植肿瘤的代谢组学和转录组学表明细菌可能影响氨基酸代谢途径
为了探索细菌对肿瘤的影响,作者收集了小鼠移植肿瘤的样本,并对其进行了基于液相色谱-质谱(LC-MS)的非靶向代谢组学研究。单变量统计分析结果显示,实验小鼠和对照小鼠移植肿瘤中共有380种代谢物存在显著差异,其中上调和下调的代谢物分别有187种和193种(图8a)。
主成分分析结果表明,两组样本数据可以很好地区分开来(图8b), 这表明两组样本之间存在不同的生物学特征。通过对比分析,用P. fungorum治疗移植瘤后,代谢产物2-丙基戊二酸、3-甲基二氧吲哚和二氢叶酸等呈上调趋势(图8c)。上调最显著的代谢物是2-丙基戊二酸和N-乙酰天冬氨酸谷氨酸。接下来,作者对差异代谢物进行了KEGG通路富集分析。结果表明,差异代谢物主要富集在丙氨酸、天冬氨酸和谷氨酸代谢途径中(图8d)。通过将功能预测分析的结果与之前的微生物组数据相结合,发现该通路中N-乙酰基-L-天冬氨酸的数量增加(图8e),这可能是由细菌作用引起的。
图8 小鼠皮下移植肿瘤的代谢组学分析