大肠杆菌是人类和多种动物中最常见的条件性致病菌,可利用许多毒力因子和效应分子破坏宿主细胞的功能介导致病,是引起肠道和肠外疾病的主要病原体,引发胃肠炎、尿路感染和血液感染
[1-2]。大肠杆菌可能通过种间和种内水平基因转移(Horizontal gene transfer,HGT)获得耐药基因
[3],对多种抗生素产生高度耐药性
[4],导致多重耐药性(Multi-drug resistance,MDR)大肠杆菌数量在世界范围内迅速增加
[5]。在临床感染治疗中,同样面临着原有抗菌药物有效性下降的问题,迫切需要开发与传统抗生素结构和作用方式不同的新型抗菌化合物对抗相关细菌引起的感染。
小分子抑制剂是一类能够与靶蛋白结合降低靶蛋白生物活性的小分子化合物,某些信号小分子通过非杀菌机制参与调控细菌生长过程,不易产生耐药性问题,可以作为抗生素发展的互补策略
[6]。基础研究中用于抗菌抑制剂的靶标主要针对细菌因子,如针对大肠杆菌脂蛋白运输途径、Ⅲ型分泌系统ATPase EscN蛋白、AraC家族转录因子Rns蛋白等,这类抑制剂一般作用于单一靶标,药物靶标的微小变化将影响药物作用效果
[7-8]。邻菲罗啉(1,10-Phenanthroline monohydrate,1,10-PHEN)由3个六元芳环组成,是一种金属螯合剂,可防止链脲佐菌素诱导染色体畸变,也是基质金属蛋白酶(Matrix metalloproteinases,MMP)的抑制剂。多项研究表明,邻菲罗啉及其衍生化合物均表现出广谱抗菌和抗真菌活性,可在金属螯合作用和生物功能之间建立紧密联系
[9-10]。邻菲罗啉可干扰白色念珠菌酶促反应中关键金属的获取和利用,强烈抑制其金属肽酶及锌依赖性金属酶水解活性,有选择性地隔离细胞代谢必需金属离子,干扰真菌的细胞稳态,导致细胞死亡
[11]。邻菲罗啉类药物与细菌结合非单一结合靶点,脱靶可能性小,已在临床上作为局部抗菌剂应用
[12]。
实验室前期研究通过虚拟筛选技术获得候选小分子抑制剂邻菲罗啉,初步研究发现邻菲罗啉对大肠杆菌具有抑菌活性,且能抑制大肠杆菌生物被膜的形成。为进一步了解邻菲罗啉对大肠杆菌细菌转录调控的影响和抗菌作用机制,本研究借助RNA-seq技术分析邻菲罗啉影响的细菌基因表达谱,为邻菲罗啉作用机制的阐释提供理论支撑。
1 材料与方法
1.1 菌株及试剂
大肠杆菌ATCC 25922由本实验室于-80 ℃条件保存。
邻菲罗啉购自上海陶术生物科技有限公司;细菌RNA快速提取试剂盒和TRUEscript SYBR Green qPCR Kit购自北京艾德莱生物科技有限公司;二甲基亚砜(DMSO)、LB琼脂培养基和LB液体培养基购自北京索莱宝科技有限公司。
1.2 试验方法
1.2.1 大肠杆菌ATCC 25922菌液制备
取出冻存的ATCC 25922菌种甘油管,划线接种于LB培养基,37 ℃恒温培养约16 h。挑取单菌落接种于LB液体培养基,37 ℃摇床200 r/min培养至OD
630为0.5,对菌液进行微量接种菌落计数
[13],调整菌液浓度为1×10
6 CFU/mL,置于4 ℃保存备用。
1.2.2 邻菲罗啉溶液的配制
准确称取邻菲罗啉25 mg,用DMSO溶解配置成190 mmoL/mL母液,使用时用LB液体培养基稀释为所需浓度。
1.2.3 邻菲罗啉的最小抑菌浓度(MIC)测定
采用微量稀释法
[14]和OD值测定法,分别测定不同浓度的邻菲罗啉对大肠杆菌的最小抑菌浓度(MIC)。使用含邻菲罗啉的LB培养基(终浓度为6.25、12.5、25、50、100、200、400 和800 µmol/L)将大肠杆菌培养物浓度调整为10
6 CFU/mL,大肠杆菌菌液为阳性对照,LB培养基为空白对照,稀释后菌液接种于96孔培养板中,37 ℃孵育24 h后使用酶标仪测量在630 nm处的OD值。
1.2.4 邻菲罗啉的最小生物被膜抑制浓度测定(MBIC)
使用含邻菲罗啉的LB培养基(终浓度为 1/32 MIC、1/16 MIC、1/8 MIC、1/4 MIC、1/2 MIC、1 MIC、2 MIC)将大肠杆菌培养物浓度调整为106 CFU/mL,接种于96孔培养板中,37 ℃静置培养48 h后取出培养物,用结晶紫染色法对生物被膜进行定量分析,使用酶标仪测量在570 nm处的OD值。
1.2.5 扫描电镜(SEM)观察邻菲罗啉处理后大肠杆菌的形态
将6 mm的细胞爬片放入12孔板底部,邻菲罗啉组加入300 µL含浓度为1/2 MIC邻菲罗啉的LB培养基及300 µL菌液,EC组(不添加邻菲罗啉),37 ℃孵育24 h。孵育后用PBS洗涤,爬片用2.5%戊二醛4 ℃固定过夜,乙醇分级脱水。样品临界点干燥后溅射涂金,在扫描电镜下观察。
1.2.6 邻菲罗啉胁迫下大肠杆菌转录组学分析
大肠杆菌在LB培养基中过夜培养,设置2个组,未处理对照组、邻菲罗啉组(添加1/2 MIC邻菲罗啉),37 ℃培养至OD630为0.6。使用细菌RNA快速提取试剂盒分离各样品中的总RNA,在1%琼脂糖凝胶上检测RNA是否降解和污染,使用Agilent 2100生物分析仪(Agilent Technologies, CA, USA)完成RNA定量,使用NanoDrop分光光度计评估质量和完整性,送至北京奥维森基因科技有限公司完成转录组测序。
1.2.7 差异表达基因的筛选、功能注释和富集分析
根据转录组测序结果,以大肠杆菌K-12 MG1655基因组为参考序列,使用Trimmomatic软件、Bowtie 2(Bowtie 2-2.3.5.1)比对软件及DESeqR包(1.12.0)
[15],对未处理大肠杆菌EC组、邻菲罗啉处理组(1/2 MIC邻菲罗啉)进行差异基因表达分析,差异基因进行筛选的阈值为:|log
2(Fold Change)|>1且
q value<0.005。对差异表达基因进行GO和KEGG富集分析。
1.2.8 邻菲罗啉对大肠杆菌差异表达基因的qPCR验证
为验证RNA-seq结果,选取大肠杆菌抗生素生物合成及生物被膜形成相关的4个显著上调和5个显著下调的基因作为目标扩增基因,根据NCBI上大肠杆菌序列设计引物。根据TRUEscript SYBR Green qPCR Kit说明书,以16S rRNA为内参基因,用qPCR对所选基因进行相对定量检测。通过2-∆∆Ct法计算目的基因mRNA的转录水平。
1.3 统计分析
应用Graphpad Prism 8.0 对所得数据进行分析,并制作结果图。获得数据以平均值±标准差(±SD)表示。*P<0.05,与未处理大肠杆菌组相比差异显著;**P<0.01,与未处理大肠杆菌组相比差异极显著。
2 结果与分析
2.1 邻菲罗啉对大肠杆菌的抗菌活性分析
最小抑菌浓度(MIC)、最小生物被膜抑制浓度(MBIC)用于评估各种抗菌化合物的抗菌功效。不同浓度邻菲罗啉作用大肠杆菌24 h后,测定其MIC及MBIC,结果如
图1所示。邻菲罗啉对大肠杆菌的MIC为100 µmol/L(
图1(a)),MBIC为50 µmol/L(图
1(b)和
1(c)),由此可知,邻菲罗啉对大肠杆菌具有良好的抑菌活性,对大肠杆菌生物被膜有明显抑制作用且呈浓度依赖性。
2.2 邻菲罗啉对大肠杆菌细胞形态的影响
扫描电子显微镜可以清晰显示细菌的表面形态和内部超微结构。邻菲罗啉作用大肠杆菌48 h后,通过电子显微镜观察细胞形态的变化,结果如
图2所示。未处理的大肠杆菌细胞形态完整,成短杆状,菌体饱满,表面较为光滑,细胞壁完整(
图2(a));经邻菲罗啉处理后的大肠杆菌细胞出现破损,菌体皱缩,表面粗糙且有多处凹陷,细胞边界模糊,细胞膜及细胞壁破坏严重(
图2(b))。结果表明,邻菲罗啉处理大肠杆菌后,可显著破坏细菌细胞的正常形态,导致细胞受到破坏。
2.3 邻菲罗啉对大肠杆菌转录组的影响
2.3.1 差异基因表达分析
利用Illumina二代高通量测序平台,采用PE150测序策略,利用Bowtie 2和 Rockhopper软件完成比对和转录本拼接分析,对所有基因进行定量分析,以未处理大肠杆菌组与邻菲罗啉处理组作为差异组合,分别获得16 379 242和16 461 344个序列,与大肠杆菌K-12的基因组比对,筛选标准为基因表达差异倍数大于2且FDR组比对,筛选显著差异的基因,结果如
图3所示,共获得显著上调基因786个,下调基因850个。
2.3.2 差异基因的GO和KEGG富集分析
邻菲罗啉处理大肠杆菌后,对差异表达基因进行GO分析,如
图4(a)所示,邻菲罗啉处理后前30个差异富集最显著的基因通路包括:在生物学过程方面主要参与有机氮化合物合成与代谢、氧化还原反应、小分子合成与代谢、蛋白质代谢等过程;在细胞组分方面主要参与细胞内、细胞内细胞器和非膜界细胞器等;在分子功能方面主要富集于催化活性、离子结合、氧化还原酶活性、结构分子活性等,这些功能分类包括上调基因和下调基因。其中,与生物素合成、青霉素结合相关基因均表达上调,与跨膜运输、对抗生素反应、转运酶活性、生殖调节、应激反应调节相关基因均表达下调。结果表明,邻菲罗啉可影响大肠杆菌催化活性和转运蛋白活性等分子功能,导致大肠杆菌的耐药和毒力表型发生变化。
对差异表达基因进行KEGG分析,结果如
图4(b)所示,大肠杆菌经邻菲罗啉处理后共有96条通路发生差异富集,其中前20条富集最显著的通路涉及多个代谢途径,包括精氨酸和脯氨酸等氨基酸合成与代谢、淀粉和蔗糖合成与代谢、柠檬酸循环等;次生代谢产物的生物合成、抗生素的生物合成、细菌趋化性、阳离子抗菌肽耐药性、β-内酰胺抗性、鞭毛集合体、ABC转运、生物被膜的形成、群体感应等多个途径,这些途径可能直接或间接与大肠杆菌的致病性及耐药性相关。
2.3.3 邻菲罗啉对大肠杆菌代谢的调节作用
邻菲罗啉可造成大肠杆菌能量代谢的紊乱。转录组分析结果显示,共有595个差异表达基因涉及多个代谢途径,如
图5所示,包括次生代谢物的生物合成、碳代谢、糖酵解与糖异生及氨基酸的生物合成等。其中涉及的基因包括
fdhDEF(编码甲酸脱氢酶)、
hycABCDE(参与
hyc和
hyp操纵子的转录抑制)、
galS(编码半乳糖超诱导因子)、
garD(编码半乳糖酸脱水酶)、
fucO(编码乙醇脱氢酶)、
hybABCDEO(参与编码氢化酶)、
lysA(参与赖氨酸合成)、
pepEN(编码丙氨酸氨基肽酶)、
argI(编码精氨酸酶)等显著下调,
dcuABC、
dctA、
glpA、
aspA、
frdABCD、
fum和
citDEF等显著表达上调,其中
fum基因编码延胡索酸水解酶可催化延胡索酸与苹果酸的相互转化,
citD、
citE和
citF基因产物分别为柠檬酸裂解酶的γ、β和α亚基,
citT基因编码产物为柠檬酸盐/琥珀酸盐抗转运剂,这些基因的差异表达表明邻菲罗啉可以影响大肠杆菌电子传递链、三羧酸循环过程。此外,
nar基因可编码硝酸盐还原酶复合物允许细菌在厌氧生长过程中使用硝酸盐作为电子受体在邻菲罗啉的刺激下发生上调,
hmp基因在有氧过程中参与NO解毒,也发生了显著下调。表明在邻菲罗啉的作用下大肠杆菌存在呼吸抑制现象。
2.3.4 邻菲罗啉对大肠杆菌环境应答的分子调控
邻菲罗啉对大肠杆菌双组份系统、ABC转运系统、生物被膜形成、群体感应、细菌趋化性等多个环境应答调控机制具有显著调节作用。转录组分析结果显示,共有91个双组份系统及85个ABC转运蛋白调控基因发生差异表达(
图6),其中
citT、
atoABE和
hyaC等47个基因显著下调,共34个基因发生了显著上调如
torACD和
narGHZ。共有26个群体感应及28个生物被膜形成相关基因发生差异表达,包括
sdiA、
flu和
bcsA等21个上调基因,
gcvA、
dgcC、
glgCP、
ygiS和
ompW等30个下调基因。表明邻菲罗啉可影响大肠杆菌的ATP形成和能量供应,降低大肠杆菌对环境的适应力。
2.3.5 邻菲罗啉对大肠杆菌铁载体的调控
邻菲罗啉对大肠杆菌铁载体调控基因的表达也产生了影响,共有6个与铁载体生物合成直接相关基因发生表达上调(
图7),
ent基因参与细菌铁载体肠杆菌素的生物合成,
suf基因可在铁限制下促进细胞外铁螯合剂中铁的摄取,异氢肟酸型铁载体(
fhu)参与络合铁的吸收,
efe参与外源性血红素铁的回收,这些基因均发生显著上调。其余间接相关如铁硫簇蛋白形成相关的基因
ftnAB、ytfE和
nuoIGH,与调节金属结合相关的基因
ydjY、
yiaY和
mocA均发生显著下调。表明邻菲罗啉可影响环境及大肠杆菌内铁离子的浓度,减少铁离子摄入及利用。
2.4 邻菲罗啉对大肠杆菌差异表达基因的qPCR验证
为验证RNA-seq分析差异表达基因的准确性,根据基因功能,选择与生物被膜形成及在耐药性产生方面具有潜在作用的9个差异表达功能基因进行了qPCR验证。结果如
图8所示,在邻菲罗啉处理组中,
degP、ldtC、
ompW、
hybA和
frdC共5个基因表达下调,
flgB、
marA、
flu和
entF共4个基因表达上调,qPCR验证结果与RNA-seq分析结果一致。由此说明,qPCR和RNA-seq分析之间的相关性及基因表达的真实性。
3 讨论
大肠杆菌可引起胃肠道感染、食物中毒、伤口细菌性感染、褥疮以及阑尾炎或憩室炎,在治疗过程中易产生耐药性,耐药谱广,耐药机制复杂,解决细菌耐药性的关键在于新型抗菌药物的筛选和研发。邻菲罗啉是一种金属螯合剂,通过干扰酶促反应中关键金属的获取和利用从而抑制胞外金属型肽酶活性。本研究表明,邻菲罗啉对大肠杆菌具有显著的抗菌活性,当浓度达到100 µmol/L可完全抑制细菌生长,50 µmol/L可显著抑制细菌生物被膜形成。使用邻菲罗啉处理大肠杆菌后进行RNA-seq,表明邻菲罗啉对大肠杆菌具有调控大肠杆菌代谢及耐药基因表达水平的作用。
细菌关键代谢途径三羧酸循环和电子传递链等途径产生多种代谢产物和能量供细菌生长所需
[16]。转录组结果显示,在邻菲罗啉的作用下大肠杆菌代谢过程中富马酸代谢(如
cit、
frd和
nar)及细菌厌氧生长过程中的相关基因(如
hmp)表达上调,因此推测邻菲罗啉使大肠杆菌的三羧酸循环过程及电子传递链受到影响,进而发生呼吸抑制。富马酸作为氧化剂在生物合成和氧化还原平衡中具有关键作用,是肠杆菌在厌氧生长过程中的关键底物,富马酸代谢过程相关基因如
citT基因主要负责柠檬酸盐的摄取与琥珀酸盐、富马酸盐或酒石酸盐的外排交换,可影响细菌三羧酸循环过程及ATP形成和能量供应
[17-19];富马酸还原酶操纵子(FRD)可以有效地催化还原反应,并允许富马酸在厌氧生长中作为电子受体发挥作用
[20]。此外,大肠杆菌
nar基因可编码硝酸盐还原酶复合物允许细菌在厌氧生长过程中使用硝酸盐作为电子受体,
hmp基因在有氧过程中参与NO解毒。这些基因表达上调均可说明邻菲罗啉影响大肠杆菌对氧的摄取,导致细菌生长所需能量供应不足,同时在有氧呼吸过程中产生有毒自由基,抑制其生长。
氨基酸的代谢是逆转耐药性的一种策略,如丙氨酸氨基肽酶PepN是大肠杆菌唯一不依赖ATP降解胞质蛋白的肽酶,可能是抵抗细菌感染的有效靶点
[21-22]。转录组结果显示,邻菲罗啉使大肠杆菌氨基酸代谢相关的基因显著上调,
fabI、
frdD、
marR和
acrR等基因与β内酰胺酶、多药外排泵系统和膜通透性的转录调节相关基因发生显著下调
[23]。表明邻菲罗啉可能促进大肠杆菌对β内酰胺类抗生素、氨基糖苷类抗生素、其他小分子抑制剂转运和敏感性,具有协同抗菌作用。本研究发现,邻菲罗啉可抑制大肠杆菌生物被膜的形成,且邻菲罗啉处理组可在扫描电镜下观察到细胞膜及细胞壁破坏严重,转录组结果也显示,与生物被膜形成相关基因发生表达下调,如可以抑制细菌的运动以减少生物被膜形成的
gcvA基因,在革兰氏阴性菌生物被膜形成、黏附、应激耐受方面发挥了重要作用的
ompW基因
[24-25]。同时,大肠杆菌在邻菲罗啉的刺激下,其双组份系统通路的正调控基因也均发生了下调。生物被膜可影响细菌对抗生素的耐受性和对机体的免疫抗性,但生物被膜的形成依赖于细菌的群体感应通信系统,破坏生物被膜可阻断细菌通讯信号
[26]。双组份系统是细胞内重要的环境信号感受器,通过感受多种调节细胞的生命活动,抗生素抗性决定簇可作为质粒和整合素在病原体之间传播,导致耐药性分布于细菌种群中
[27]。以上结果均可说明,邻菲罗啉可降低大肠杆菌对环境的适应力及对抗生素的耐药性。
通过转录组分析可知,在邻菲罗啉的作用下使大肠杆菌铁硫簇蛋白形成相关基因与调节金属结合相关的基因均发生显著下调,参与细菌铁载体肠杆菌素合成的基因发生显著上调,表明大肠杆菌在邻菲罗啉刺激下发生铁限制并增强了细胞外铁螯合剂、络合铁中铁的摄取及外源性血红素铁的回收。铁是细菌生存所必需的营养素,细菌分泌铁载体,并在外膜上表达特异性转运蛋白,该蛋白负责摄取铁载体复合物,产生特定环境信号(如局部铁浓度)的毒力因子(如铁载体、大肠杆菌素和铁载体-微菌素),使大肠杆菌能够适应肠道和/或血液的竞争环境
[28]。其中肠杆菌素可被用于传递分子的载体,增强革兰氏阴性菌外膜通透性,运送小分子物质到大肠杆菌细胞质,影响细菌毒力,也可杀灭无法分泌特异性肠杆菌素免疫蛋白的近缘细菌,有调节菌群数量的作用
[29-30]。邻菲罗啉有较强的金属螯合作用,可与金属离子形成稳定配合物。因此,推测邻菲罗啉通过降低环境及大肠杆菌内铁离子的浓度,减少细菌铁离子摄入及利用,造成细菌铁死亡,同时促进大肠杆菌铁载体肠杆菌素的生物合成,提高细菌对抗生素的运输效率,增强大肠杆菌对邻菲罗啉的敏感性。
4 结论
邻菲罗啉对大肠杆菌具有显著的抑制作用,可能通过螯合金属离子并形成稳定的螯合物,影响环境及大肠杆菌内铁离子的摄入及利用而抑制大肠杆菌生长;同时,邻菲罗啉可能刺激大肠杆菌在有氧呼吸过程中产生有毒自由基,对细菌产生毒性从而发挥抗菌作用。本研究为寻找抗生素替代药物提供新思路。
2023年北京农学院学位与研究生教育改革与发展项目(5056516030/003)
北京市优秀人才培养资助项目(2017000026833ZK18)
现代农业产业技术体系北京市创新团队BAIC08-2023
北京农学院青年科学基金项目(QNKJ202107)