Metabonomics analysis of ovaries of Coilia nasus in seawater and freshwater based on liquid chromatography-mass spectrometry
-
摘要: 为探究海、淡水养殖环境对刀鲚 (Coilia nasus) 卵巢发育的影响,采用非靶向代谢组学的方法,测定了海、淡水环境下刀鲚卵巢中代谢产物的差异情况,并与基因组百科全书 (KEGG) 数据库进行比对,找出相对应的代谢通路并分析其原因。结果表明,海水组和淡水组样品共鉴定出47种差异代谢物 (P<0.05、FC>1、VIP>1),与海水组相比,淡水组表达差异倍数最明显的为碳环血氧烷A2 (Carbocyclic thromboxane A2)、半乳糖神经酰胺(Galactosyl ceramide),差异倍数分别为10.40、2.78倍;与海水组相比,淡水组卵巢组织内皮质醇升高了1.61倍;对47种差异代谢物进行KEGG分析发现,变化显著的通路有氨酰-tRNA的生物合成和嘧啶代谢通路 (P<0.05),皮质醇、氨酰-tRNA的生物合成通路、嘧啶代谢通路和鞘磷脂代谢通路可能与刀鲚生殖洄游过程中卵巢发育有关。Abstract: In order to clarify the effects of seawater and freshwater on the ovary development of Coilia nasus, we analyzed their differences by using non-targeted metabolomics, and compared with database of KEGG directly to find out the corresponding metabolic pathways, then analyzed its causes. The results show that a total of 47 metabolites had significant difference between the two groups (P<0.05, FC>1, VIP>1). Compared with the seawater group, the most significant differences in the expression were carbocyclic thromboxane A2, Galactosylceramide, and their differences were 10.40 and 2.78 times, respectively. The cortisol in the ovarian tissue of the freshwater group increased by 1.61 times. According to the analysis on KEGG metabolic pathways of 47 different metabolites, the biosynthesis of aminoacyl-tRNA and pyridine metabolic pathways changed in the seawater and freshwater environments significantly (P<0.05). The biosynthesis pathway of cortisol, aminoyl-tRNA, pyrimidine metabolism pathway and sphingo-lipid metabolism pathway may be related to the ovarian development during the reproductive migration of C. nasus.
-
Keywords:
- Coilia nasus /
- LC-MS /
- Salinity /
- Ovary /
- Metabolomics /
- Metabolic pathways
-
哈维弧菌 (Vibrio harveyi) 呈短杆状,为革兰氏阴性菌,能够感染无脊椎动物[1]、脊椎动物以及人类[2]。它是海水养殖动物的常见致病菌,曾引起点带石斑鱼 (Epinephelus coioides)[3]、中国花鲈 (Lateolabrax maculatus)[4]等多种鱼类大量死亡,造成巨大经济损失[5]。哈维弧菌黏附在宿主表面后会产生生物膜,在多种毒力因子参与下入侵到宿主体内,引起发病[6]。目前发现,哈维弧菌的毒力因子包括细胞外产物、膜蛋白、生物膜、黏附素及多种分泌系统等[7]。
溶血素 (Hemolysin) 是细菌的一种重要毒力因子,其能溶解宿主红细胞膜,导致红细胞内容物溢出,引起溶血现象,在感染过程中发挥重要作用[8]。弧菌溶血素主要包括热稳定性直接溶血素 (Tdh)、热稳定性溶血素 (δ-Vph)、热不稳定性溶血素 (Tlh) 和E1 Tor溶血素 (HlyA)[9]。哈维弧菌溶血素 (Vhh) 是Tlh家族的成员[10],它是一种单体蛋白,分子量约为44.3 kD[11],具有溶血磷脂酶活性[12]。Vhh的氨基酸序列与副溶血弧菌 (V. parahemolyticus)Tlh溶血素[13]、创伤弧菌 (V. vulnificus) 磷脂酶VpL (Genbank AF291424)、拟态弧菌 (V. mimicus) 卵磷脂酶(pHL)[14]和霍乱弧菌 (V. cholerae) 卵磷脂酶(Lec)[15]的序列相似性分别为85.6%、75.6%、65.3%和64.3%。
国内外学者对Vhh做了相关的功能研究,Bai等[11]纯化了Vhh蛋白,研究其对牙鲆 (Paralichthys olivaceus) 鳃细胞和红细胞的影响;Zhong等[12]通过大肠杆菌 (Escherichia coli) 表达Vhh,研究了不同温度和pH对其溶血活性、磷脂酶活性、细胞毒性及对鱼体致病性的影响。上述研究集中在体外活性检测,仍缺乏对vhh基因在菌体内作用机制的认识。本研究通过同源重组方法构建哈维弧菌vhh基因缺失突变株,并分析缺失株相关的生物学特性,为进一步了解vhh的基因功能提供基础资料。
1. 材料与方法
1.1 材料
1.1.1 菌株、质粒和引物
哈维弧菌V. harveyi 345、自杀质粒pSW7848[16]、大肠杆菌E. coli Π3813[17]以及大肠杆菌E. coli GEB883[18]保存于本实验室;本研究所用菌株和质粒见表1,引物见表2。
表 1 本研究所用菌株和质粒Table 1. Strains and plasmids used in this study菌株、质粒
Strain, plasmid相关特征
Relevant characteristics来源
Source菌株 Strain V. harveyi 345 分离自中国南部沿海患病石斑鱼肾脏 NCBI database CP025537.1-CP025540.1 V. harveyi 345-Δvhh V. harveyi 345缺失vhh基因 本文构建 E. coli Π3813 lacIQ,thi1,supE44,endA1,recA1,hsdR17,gyrA462,zei298::Tn10[Tc],ΔthyA::erm-pir116;自杀质粒pSW7848的中间宿主; 文献[17] E. coli GEB883 大肠杆菌野生株K12 ΔdapA::ermpir RP4-2 ΔrecAgyrA462,zei298::Tn10;接合作用供体菌 文献[18] 质粒 Plasmid pSW7848 氯霉素抗性;自杀质粒,R6K起点,复制需要Pir蛋白,具ccdB毒性基因 文献[16] pSW7848-Δvhh 氯霉素抗性;包括vhh的上下游同源臂 本文构建 表 2 本研究所用引物Table 2. Primers used in this study引物
Primer序列 (5'–3')
Sequence用途
ApplicationpSW7848-F GTCTGATTCGTTACCAATTATGACAAC 扩增pSW7848片段 pSW7848-R GAATTCGATATCAAGCTTATCGATAC vhh-UP-F aagcttgatatcgaattcTACCGTGAGCGTTGTTATTAC 扩增vhh上游片段 vhh-UP-R ctctggatGATTCACTCTCTGTGATGATTATC vhh-DOWN-F gagtgaatcATCCAGAGCAACACGGTTTTG 扩增vhh下游片段 vhh-DOWN-R ttggtaacgaatcagacAAGATTGATGTTTTCTATTGAGACACAAAAG Del-check-pSW7848-F TCACTGTCCCTTATTCGCACC 验证pSW7848重组质粒是否构建成功 Del-check-pSW7848-R CTGCTTTTGAGCACTACCCG Δvhh-check-F TTGAAAGTTGAACGCGGTAC 检验vhh基因是否敲除成功 Δvhh-check-R GAGGTATGTCTTGAGAGGTCAA gyrA-354-F AGGTGTTCGCGGTATGAAGC 荧光定量管家基因gyrA扩增 gyrA-354-R CGTTGGGTATTCCGCTAGTTC uvrA-354-F ATCGTGGTGACATTCGCGTT 荧光定量管家基因uvrA扩增 uvrA-354-R TCGCAAGCTCTTCACGTACC recA-354-F TGCACGTTCTGGTGCTATTG 荧光定量管家基因recA扩增 recA-354-R AGCTTACGCATCGCTTGAGA k14645-354-F ATAACGCCTTTGATGCCTTCC 荧光定量vhh上游基因k14645扩增 k14645-354-R AACGATCTGAGTGGTAACCGAGT vhh-F CACTTATGTCCGCTGCTGGTAC 荧光定量vhh基因扩增 vhh-R CTTGAAAGAAACCGAACTCCAC na-F CAGAAATCACTAAGCGAGCGACAT 荧光定量vhh下游基因na扩增 na-R GTGCGTTTAACTCATCCAACAAGG 1.1.2 主要试剂
大肠杆菌E.coli Π3813生长所需的2'-Deoxythymidine (Thy) 以及用于铁离子吸收测试的2-2'-联吡啶 (2,2'-Dipyridyl,DIP),购自生工生物工程 (上海) 股份有限公司;反转录试剂盒PrimeScriptTM RT reagent Kit with gDNA Eraser (Perfect Real Time) 和荧光定量反应试剂盒TB GreenTM Premix Ex TaqTM Π (Tli RNaseH Plus),购自TaKaRa Bio Ine (日本);抗生素纸片,购自杭州天和微生物试剂有限公司;PBS缓冲液,购自赛默飞世尔 (苏州) 仪器有限公司;脱纤维棉羊红细胞,购自北京索莱宝科技有限公司。
1.1.3 主要培养基
LB (Luria-bertani) 培养基购自广东环凯生物科技有限公司。1 L LB含胰蛋白胨10 g,酵母提取物5 g,氯化钠 (NaCl) 10 g。1 L固体培养基含琼脂15 g。在LB培养基基础上将NaCl质量浓度提高到30 g·L−1即配制得到LBS (Luria-bertani-salt) 培养基。
1.2 方法
1.2.1 vhh基因序列分析
将本实验室V. harveyi 345基因组测序获得的vhh基因 (CU052_24955)ORF序列[19]进行BLAST分析,选取其在不同细菌中的同源序列进行多序列比对,并用MEGA 7.0软件采用NJ邻接法构建系统进化树。
1.2.2 vhh缺失株的构建
构建vhh缺失株参考Deng等[20]的方法并进行了相关修改,vhh的序列在NCBI数据库下载 (登陆号CU052_24955),扩增线性化自杀质粒pSW7848使用引物对pSW7848-F/R,扩增上下游同源臂片段使用引物对vhh-UP-F/R和vhh-DOWN-F/R,等温组装线性化质粒片段和上下游同源臂片段,对重组自杀质粒进行PCR和测序检测使用引物对Del-check-pSW7848-F/R,组装成功的质粒大小约为5.3 kb,将其命名为pSW7848-Δvhh。然后重组自杀质粒pSW7848-Δvhh通过接合作用将供体菌大肠杆菌GEB883-pSW7848-Δvhh送入到受体菌哈维弧菌野生株V. harveyi 345。基于自杀质粒pSW7848在哈维弧菌中不能复制,具有氯霉素抗性,毒性基因ccdB能够在受阿拉伯糖诱导后表达的特点,使用氯霉素筛选后,重组质粒和哈维弧菌基因组上的vhh上下游同源臂的相同片段发生第一次同源重组;使用阿拉伯糖进行第二次筛选,其发生第二次同源重组,最终获得vhh缺失株,将其命名为V. harveyi 345-Δvhh。重组自杀质粒构建和同源重组筛选缺失突变株流程见图1。
1.2.3 vhh基因及其上下游基因转录水平分析
挑取野生株V. harveyi 345和突变株V. harveyi 345-Δvhh于28 ℃、200 r·min–1振摇培养16 h过夜,稀释至OD600为0.01,培养至对数早期 (OD600=0.5)。RNAiso Plus裂解细菌,提取总RNA,并将其进行反转录。利用SYBR Premix Ex TaqTM Π进行荧光定量PCR扩增,特异性引物 (表2) 分别检测vhh上游丝氨酸蛋白酶基因K14645、溶血素基因vhh和vhh下游未知功能基因na的表达水平,3个管家基因 (recA、uvrA和gyrA) 作为内参。对每个单菌落进行3次技术重复,通过2−ΔΔCT方法[21]分析各基因的相对表达量。
1.2.4 菌株生长测试
将V. harveyi 345和缺失突变株V. harveyi 345-Δvhh的单克隆分别接种至LBS液体培养基,于28 ℃、200 r·min–1振摇培养16 h过夜,稀释1 000倍于无菌锥形瓶中,在相同条件下培养,利用分光光度计 (BIO-RADSmart-SpecTM Plus) 测定不同时间点下细菌的OD600,绘制生长曲线。
1.2.5 对绵羊红细胞的溶血活性测定
同1.2.4过夜培养菌液,1 500 r·min–1离心5 min,弃上清。加PBS洗涤,离心后弃上清,重复两次。用PBS配置成4%绵羊红细胞悬液。参照Takahiro等[22]方法稍作改变,将过夜培养的V. harveyi 345和V. harveyi 345-Δvhh用PBS洗涤并调到吸光度OD600为3.0。取50 μL菌液加入到500 μL 4%脱纤维棉羊红细胞中,混匀后放置28 ℃孵育5 h。然后12 000 r·min−1离心1 min,取200 μL上清液测OD600,分别加入PBS和ddH2O作为阴性和阳性对照,溶血活性以样品相对于阳性对照OD600比值表示。
1.2.6 游动性和涌动性测试
同1.2.4过夜培养菌液,用LBS培养基将其调整到吸光度OD600为3.0,取3 μL点样到含0.3%琼脂和1.5%琼脂的LBS平板上测试其游动性及涌动性,进行3次技术重复,放置28 ℃分别培养16 h和24 h,记录菌斑直径。
1.2.7 胞外蛋白酶分泌测试
同1.2.4过夜培养菌液,同1.2.6调整培养液吸光度,取3 μL点样到LBS平板 (含1.5%脱脂奶粉) 上,进行3次技术重复,放置28 ℃培养24 h,测量菌落直径以及蛋白质透明分解圈直径,蛋白质透明分解圈直径与菌落直径的比值表示蛋白酶的分泌能力。
1.2.8 对过氧化氢(H2O2)、铜离子(Cu2+)抗性及对铁离子的吸收测试
参照Deng等[23-24]的方法,稍作改变,同1.2.4过夜培养菌液,同1.2.6调整培养液吸光度,用LBS培养基将其稀释到最高倍数108,取3 μL点样到3种LBS平板 [分别含有0.000 2% H2O2、0.2 mmol·L–1 硫酸铜(CuSO4)和120 μmol·L–1 DIP] 上,不含添加剂的LBS平板作为对照,放置28 ℃培养24 h。查看各平板上的菌株生长情况。
1.2.9 抗生素的敏感性测试
同1.2.4过夜培养菌液,取600 μL于10 mL 0.9%生理盐水中,全部倾入平板中润湿培养基表面,弃去菌液,静置10 min,将各抗生素纸片置于晾干的平板上,放置28 ℃培养24 h,记录抑菌圈直径。
1.2.10 生物膜测试
同1.2.4过夜培养菌液,同1.2.6调整培养液吸光度,生物膜测试参考顾丹[25]的方法,取50 μL菌液接种到装有5 mL的LBS液体培养基的无菌玻璃试管中,放置28 ℃培养48 h后,加入200 μL结晶紫染色5 min,将菌液倒掉并用水反复冲洗几次,加入200 μL 33%冰乙酸溶解结晶紫,充分溶解后,吸取200 μL至96孔板,利用酶标仪 (SUNRISE吸光酶标仪) 在600 nm处测定其吸光值。
1.2.11 数据处理与分析
利用软件IBM SPSS Statistics 19统计分析实验数据,显著性差异用P<0.05表示。
2. 结果
2.1 vhh系统进化
哈维弧菌V. harveyi 345溶血素基因vhh片段大小为1 257 bp,编码419个氨基酸;将vhh核苷酸序列与其他弧菌溶血素核苷酸序列用MEGA7.0构建NJ系统进化树 (图2)。结果显示哈维弧菌V. harveyi 345 (vhh) 和哈维弧菌V. harveyi GDH 11388 (vhhA)、哈维弧菌V. harveyi ATCC 35084 (vhhB) 同源性最高,聚为一支,与创伤弧菌V. vulnificus (vpl)、溶藻弧菌V. alginolyticus R4 (tlh) 和霍乱弧菌V. cholerae (lec) 均有较高的同源性,但与霍乱弧菌V. cholerae (hlyA) 同源性较低。
2.2 构建vhh缺失突变株
扩增得到的线性化质粒片段大小约3.3 kb,上下游同源臂片段大小均约为1 kb (图3-a、图3-b),等温组装质粒线性化片段和上下游同源臂片段后,PCR扩增出2 307 bp的片段 (图3-c),重组质粒pSW7848-Δvhh构建成功。重组自杀质粒上的vhh上下游同源臂片段通过接合作用与野生株V. harveyi 345的vhh上下游同源臂片段发生两次同源交换,这可能使染色体上的vhh缺失,利用引物对Δvhh-check-F/R分别对野生株和vhh缺失株进行PCR检测,电泳结果分别出现1 658 bp和628 bp两条带 (图3-d),测序发现vhh序列成功缺失,获得vhh缺失突变株。
图 3 pSW7848 PCR线性片段 (a)、vhh上下游同源臂片段 (b)、重组自杀质粒的PCR鉴定扩增片段 (c) 及vhh缺失株的PCR鉴定 (d)M1. DNA Marker DL5 000;1. pSW7848线性片段;M2. DNA Marker DL2 000;2/3. vhh上/下游片段;M3. DNA Marker DL5 000;4. 重组子 pSW7848-Δvhh检测片段;M4. DNA Marker DL2 000;5. 野生株V. harveyi 345;6. vhh候选突变株Figure 3. Linearized segment of pSW7848 (a), upstream and downstream homologous segments of vhh (b), segment of identification of recombinant suicide vector by PCR (c) and identification of vhh deletion mutant by PCR (d)M1. DNA Marker DL5 000; Lane 1. Linearized segment of pSW7848; Lane M2. DNA Marker DL2 000; Lane 2/3. Up/down segments of vhh; Lane M3. DNA Marker DL5 000; Lane 4. Segment of identification of the recombinant suicide vector by PCR; Lane M4. DNA Marker DL2 000; Lane 5. Wild type V. harveyi 345; Lane 6. Deletion candidate of vhh2.3 vhh基因缺失对邻近基因的表达影响
为研究vhh基因的缺失对邻近基因的mRNA相对表达量的变化,检测上游丝氨酸蛋白酶基因k14645、溶血素基因vhh和下游未知功能基因na。突变株V. harveyi 345-Δvhh的vhh基因未表达,与野生株V. harveyi 345相比,表达差异极显著 (P<0.01),k14645和na基因的表达与野生株相比,差异不显著 (P>0.05,图4)。结果显示,vhh基因缺失不影响邻近基因的表达。
2.4 vhh对哈维弧菌在LBS培养基中生长的影响
在LBS培养基中,vhh缺失株与野生株生长曲线几乎重合,缺失株的生长相对野生型无显著差异 (P>0.05,图5)。
2.5 vhh基因缺失对哈维弧菌溶血活性的影响
阳性对照对绵羊红细胞的溶血百分比为100%,阴性对照几乎没有发生溶血,均与预期结果一致。通过比较哈维弧菌菌株对绵羊红细胞的溶血百分比,发现野生株V. harveyi 345和突变株V. harveyi 345-Δvhh均无溶血活性,并且突变株相对野生株对绵羊红细胞的溶血百分比差异不显著 (P>0.05,图6)。
2.6 vhh对哈维弧菌运动性的调控
野生株V. harveyi 345在0.3%琼脂LBS平板上没有游动,在1.5%琼脂LBS平板上也无涌动,但是突变株V. harveyi 345-Δvhh在这两种平板上均表现出明显的游动及涌动,其游动与涌动菌落直径明显变大,相对野生株差异极显著 (P<0.01,图7-a、图7-b)。
2.7 vhh对哈维弧菌胞外蛋白酶分泌的影响
vhh缺失后,在1.5%脱脂奶粉LBS平板上,胞外蛋白酶透明分解圈与野生株一致,缺失株胞外蛋白酶的分泌相对野生株无显著差异 (P>0.05),说明两者具有相同的胞外蛋白酶分泌能力 (图8-a、图8-b)。
图 8 哈维弧菌V. harveyi 345和V. harveyi 345-Δvhh在LBS+1.5%脱脂奶粉平板上胞外蛋白酶的分泌 (a) 和统计分析 (b)内箭头指向菌斑边缘,外箭头指向蛋白酶分泌圈边缘Figure 8. Extracellular protease activities of V. harveyi 345 and V. harveyi 345-Δvhh mutant on LBS agar plates containing 1.5% skimmed milk (a) and statistical analysis (b)The inner arrow points to the edge of the plaque, and the outer arrow points to the edge of the proteolytic ring.2.8 对H2O2、Cu2+抗性以及对铁离子的吸收测试
在添加有H2O2、CuSO4以及DIP的3种LBS平板上,随着稀释倍数增加,突变株和野生株菌落数均减少,到108稀释倍数时,菌落数明显变少;两者对H2O2、CuSO4抗性以及铁离子的获取无显著差异 (图9)。
2.9 vhh对哈维弧菌抗生素敏感性的影响
野生株和突变株均对利福平、多西环素、万古霉素等5种抗生素耐药,对呋喃唑酮、阿莫西林、氯霉素等6种抗生素敏感。突变株对氟苯尼考敏感,而野生株则表现中介敏感,但两者的抑菌圈直径差异不显著 (P>0.05);野生株对妥布霉素、庆大霉素和红霉素均敏感,而突变株都表现为中介敏感,两者的抑菌圈直径差异不显著 (P>0.05,表3)。
表 3 哈维弧菌V. harveyi 345和V. harveyi 345-Δvhh对各抗生素抗性Table 3. Antibiotics resistance of V. harveyi 345 and V. harveyi 345-Δvhh抗生素
Antibiotic浓度
Concentration/(μg·片−1)抑菌圈直径 Size of inhibition zone/mm V. harveyi 345 V. harveyi 345-Δvhh 利福平 Rifampin 5 7.89±0.60/R 7.89±0.78/R 多西环素 Doxycycline 30 7.11±0.33/R 7.00±0.00/R 万古霉素 Vancomycin 30 0.00/R 0.00/R 复方新诺明 Compound sulfamethoxazole 1.25 0.00/R 0.00R 四环素 Tetracycline 30 0.00/R 0.00/R 妥布霉素 Tobramycin 10 14.67 ±1.32/S 13.22±0.97/I 氟苯尼考 Florfenicol 30 17.22±0.67/I 16.78±2.05/S 庆大霉素 Gentamicin 10 15.67±1.32/S 14.89±1.69/I 呋喃唑酮 Furazolidone 100 19.78±1.39/S 18.00±2.40/S 阿莫西林 Amoxicillin 20 22.56±1.67/S 21.00±1.32/S 氯霉素 Chloramphenicol 30 22.22±0.97/S 21.89±1.05/S 红霉素 Ergomycin 15 23.44±0.53/S 21.22±1.79/I 麦迪霉素 Medicamycin 30 19.11±0.33/S 19.00±1.12/S 环丙沙星 Ciprofloxacin 5 32.33±1.32/S 32.22±0.83/S 诺氟沙星 Norfloxacin 10 28.00±1.12/S 29.33±0.50/S 注:S. 敏感;I. 中敏;R. 耐药 Note: S. Susceptible; I. Intermediate; R. Resistance 2.10 vhh对哈维弧菌生物膜形成的影响
野生株V. harveyi 345和突变株V. harveyi 345-Δvhh在玻璃试管中静置培养,发现两者管壁上都形成了少量生物膜,并且结晶紫染色后的管壁上出现宽度几乎一样的染色圈;统计分析结果显示,突变株V. harveyi 345-Δvhh的生物膜形成量略微高于野生株V. harveyi 345,但两者的生物膜形成能力并无显著差异 (P>0.05,图10-a、图10-b)。
图 10 哈维弧菌V. harveyi 345和V. harveyi 345-Δvhh的生物膜形成能力.a-1、a-2分别为V. harveyi 345、V. harveyi 345-Δvhh培养液;a-3、a-4分别为结晶紫染色;b. 定量结晶紫染色后聚苯乙烯微量滴定板,OD600测定生物膜形成量Figure 10. Biofilm formation ability of V. harveyi 345 and V. harveyi 345-Δvhha-1 and a-2 are V. harveyi 345 and V. harveyi 345-Δvhh culture, respectively; a-3 and a-4 are crystal violet stainings; b. Quantitative crystal violet stained polystyrene microtiter plate; OD600 determination of biofilm formation3. 讨论
多种细菌可分泌溶血素,该物质能破坏宿主红细胞膜,造成细胞内容物溢出,导致溶血现象[26]。在作用机理上,可把溶血素分为成孔蛋白和磷脂酶[27]。其中以Tdh为代表的成孔溶血素常以可溶形式分泌,通过组装成跨膜孔发挥其催化作用;而磷脂酶具有两亲性质,作用类似去垢剂分子,直接作用细胞膜,使磷脂分子水解,导致细胞膜损伤。Bai等[11]证实,哈维弧菌溶血素Vhh是溶血磷脂酶,属于脂肪酶家族。本研究发现,哈维弧菌野生株及vhh缺失突变株对绵羊红细胞均没有溶血活性,这与Rattanama等[28]的研究类似,其在哈维弧菌分离株中均检测到vhh基因,但部分菌株并未表现出溶血活性。类似的,Fan等[29]发现部分具有完整hlyA基因的霍乱弧菌对绵羊红细胞也未表现出溶血活性,向非溶血性霍乱弧菌中转入hlyA过表达质粒,这些菌株则出现溶血活性,暗示HlyA分泌正常。这些结果说明,在非溶血性菌株中,hlyA基因不表达或表达量太低。本研究的vhh基因转录水平较高,但翻译水平未知,对绵羊红细胞也不表现溶血活性,推测可能存在Vhh蛋白未翻译或者翻译但分泌不正常,具体机制仍需进一步阐明。
胞外蛋白酶主要是金属蛋白酶或丝氨酸蛋白酶,他们能催化蛋白质水解成氨基酸等多种有机物质用于病原体的生长和繁殖,加速其扩散到宿主组织中,在感染过程中发挥重要作用[30]。本研究溶血素基因vhh的缺失并未导致哈维弧菌胞外蛋白酶活性的改变,暗示vhh基因可能不影响胞外蛋白酶的分泌。
应激反应是细菌对外界环境的响应,它能帮助细菌个体感知种群中其他细胞的状态并做出反应,增强细菌防御不利环境的能力。氧化应激是细菌在海洋环境中可能遇到的情况,细菌本身代谢和侵染宿主过程中会产生大量活性氧物质 (ROS),这些ROS分子可能会损害病原体本身,致病菌为了能生长繁殖和成功定植宿主,必须要有一套有效的氧化应激系统来清除这些有害分子,比如产生抗氧化剂过氧化氢酶,该酶将H2O2转化为H2O和O2来对其进行解毒[31]。Chen等[32]报道溶藻弧菌 (V. alginolyticus) 超氧化物歧化酶B基因sodB缺失突变体对H2O2更加敏感,这说明sodB基因的缺失降低了溶藻弧菌对外界环境中氧化自由基的防御能力,即sodB正调控抗氧化剂的分泌。本研究中vhh基因的缺失并没有改变哈维弧菌对H2O2的耐受力。过量的铜离子催化大量羟基自由基的形成,对细胞产生毒性,摄铁系统也是哈维弧菌重要致病因子之一[25]。本研究发现vhh基因未改变对铜离子的耐受能力和对铁离子的吸收利用能力。细菌为适应复杂多变的环境而形成生物膜,生物膜由蛋白质、核酸、脂质和胞外多糖组成,可有效帮助细菌抵抗环境中的不利因素[33]。当致病菌定植在宿主表面后,为使其不受宿主免疫系统和抗生素的影响而形成大量生物膜,提高病原体感染宿主的能力[34]。鞭毛介导的运动帮助细菌附着在宿主细胞表面,促进其形成生物膜[35]。本实验发现vhh基因的缺失提高了哈维弧菌的运动性,但没有改变其生物膜的形成能力,暗示哈维弧菌运动性可能与其生物膜的形成无关。
抗生素是抑制或杀死细菌生长繁殖的化学物质,大多数是由细菌合成的天然化学武器,可杀死周围环境中的其他微生物并保持细菌体内平衡[36]。大多数抗生素通过抑制细胞壁合成、破坏细胞膜完整性、防止DNA、RNA和蛋白质合成以及阻断细胞基本代谢途径而发挥抑菌作用。滥用抗生素会导致细菌产生耐药性[37],这是细菌对抗菌药物的自然反应,是适应环境的一种表现。Das等[36]发现某些rRNA甲基转移酶的缺乏会引起细菌产生抗生素耐药性,并且rulC基因缺失后能显著提高霍乱弧菌和其他肠道病原体对克林霉素、利奈唑胺和泰妙菌素的耐药性。本实验发现vhh基因不会显著影响哈维弧菌对所测试抗生素的耐药性,表明vhh不参与调控哈维弧菌对这些抗生素的耐药过程。
细菌能够独立响应环境,通过包括鞭毛和菌毛在内的物理附属物游动。细菌运动的机制主要取决于它的鞭毛数量及其在细胞表面的分布[38]。哈维弧菌具有两套鞭毛系统[39],一套主要由一根能够在液体环境中泳动的极性鞭毛组成;另一套由多根侧生鞭毛组成,使细菌能在固体表面上进行群集运动。鞭毛结构由大量严格受调控的基因编码,上游基因编码调节运动的转录激活因子,中游基因编码鞭毛输出系统和钩基体结构成分,下游基因编码鞭毛细丝、动力形成器和趋化信号转导蛋白[40]。Yang等[41]证实哈维弧菌群体感应 (Quorum sensing,QS) 不仅能调节鞭毛中下游基因,还可调节鞭毛上游转录激活因子flaK,说明哈维弧菌QS能通过控制鞭毛转录因子的表达来调控运动。本研究vhh基因缺失株运动能力增强,暗示vhh可能通过作用于QS系统,降低了鞭毛合成相关基因的表达,从而负调控哈维弧菌的运动性。细菌胞外聚合物 (EPS) 由多糖、蛋白质和脂质等组成,这种基质可以帮助细胞抵抗干旱、抗生素、重金属和紫外线辐射等不利环境因素,据报道运动性与EPS的分泌能力呈负相关[42]。本研究vhh基因缺失引起菌株运动性增强,暗示vhh基因可能正调控EPS的分泌,下一步将分析哈维弧菌野生株和vhh缺失突变株鞭毛合成及QS系统相关基因的表达,并量化EPS的分泌,以期为探究哈维弧菌溶血素基因vhh与运动性之间的联系提供新证据。
-
表 1 差异代谢物信息表
Table 1 Differential metabolites information sheet
差异代谢物
Differential metabolite分子式
Molecular formula变量投影重要度
VIPFC
Fold changeP 变化趋势
Variation trend磷脂酰肌醇 PI [20:4(5Z,8Z,11Z,14Z)/0:0] C29H49O12P 1.58 0.93 0.017 下降 谷氨酰胺色氨酸 Glutaminyltryptophan C16H20N4O4 1.05 0.96 0.018 下降 L-异亮氨酸 L-Isoleucine C6H13NO2 1.33 0.96 0.005 下降 木麻黄6-α-D-葡萄糖苷 Casuarine 6-alpha-D-glucoside C14H25NO10 1.13 0.94 0.043 下降 甘油一脂 MG(10:0/0:0/0:0) C13H26O4 3.59 0.71 0.025 下降 皮质醇 Cortisol C21H30O5 3.28 1.61 0.008 上升 碳环血氧烷A2 Carbocyclic thromboxane A2 C22H36O3 2.88 10.40 0.007 上升 (±)9-十八碳二烯酸 (±)9-HPODE C18H32O4 2.90 0.58 0.020 下降 N-棕榈酰蛋氨酸 N-palmitoyl methionine C21H41NO3S 1.86 0.75 0.022 下降 13-羟基十八酸 13-hydroxyoctadecanoic acid C18H36O3 1.95 0.86 0.005 下降 咖啡酰环戊醇 Caffeoylcycloartenol C39H56O4 2.07 0.81 0.036 下降 11,13-二十碳二烯酸 15-OxoEDE C20H34O3 2.51 0.70 0.032 下降 2,3-二氢苯并呋喃 2,3-dihydrobenzofuran C8H8O 1.19 1.08 0.033 上升 4-甲酰基吲哚 4-formyl indole C9H7NO 1.34 0.83 0.040 下降 红花素C Safflomin C C30H30O14 2.04 0.85 0.003 下降 6-[(2-羧基乙酰基)氧]-3,4,5-三羟基氧烷-2-羧酸6-[(2-carboxyacetyl)oxy]-3,4,5-trihydroxyoxane-2-carboxylic acid C9H12O10 2.07 1.12 0.001 上升 异戊二烯 Isoputreanine C7H16N2O2 2.33 0.84 0.001 下降 脯氨酸 L-Proline C5H9NO2 1.80 0.91 0.003 下降 羟脯氨酸 Hydroxyprolyl-hydroxyproline C10H16N2O5 3.13 0.72 0.007 下降 5-O-阿魏酰黑精 5-O-Feruloylnigrumin C21H25NO10 1.0 0.95 0.037 下降 壬二酸 Azelaic acid C9H16O4 1.12 0.96 0.012 下降 尿苷 Uridine C9H12N2O6 1.35 0.95 0.001 下降 视黄酯 Retinyl ester C20H30O2 1.24 0.96 0.004 下降 顺-9,10-环氧硬脂酸 cis-9,10-Epoxystearic acid C18H34O3 1.66 0.92 0.025 下降 异柠檬酸盐 Isocitrate C6H8O7 1.55 1.12 0.047 上升 戊二酸 Glutaric acid C5H8O4 2.53 1.38 0.015 上升 6-脱氧噬菌体胺 (6-脱氧花青) 6-Deoxyfagomine C6H13NO2 1.36 0.94 0.012 下降 乙酰-L-酪氨酸 Acetyl-L-tyrosine C11H13NO4 1.38 0.87 0.010 下降 γ-谷氨酰鸟氨酸 Gamma glutamyl ornithine C10H19N3O5 1.87 1.58 0.023 上升 对茴香酸异戊酯 (异戊基异茴香酸酯) Isoamyl p-anisate C13H18O3 1.76 0.67 0.041 下降 3-氧十二酸 3-Oxododecanoic acid C12H22O3 1.01 0.92 0.043 下降 C-2神经酰胺 C-2 Ceramide C20H39NO3 2.21 0.75 0.028 下降 磷脂酰乙醇胺 PE(15:0/16:1(9Z)) C36H70NO8P 1.42 1.10 0.030 上升 半乳糖神经酰胺 Galactosyl ceramide (d18:1/14:0) C38H73NO8 1.97 2.78 0.045 上升 N-棕榈酰甘氨酸 N-Palmitoyl glycine C18H35NO3 1.97 0.86 0.000 下降 二十碳五烯酸 Eicosapentaenoic acid C20H30O2 1.01 0.98 0.004 下降 溶血磷脂酰乙醇胺 LysoPE(0:0/20:2(11Z,14Z)) C25H48NO7P 1.33 1.12 0.037 上升 9-羟基癸酸 9-Hydroxydecanoic acid C10H20O3 1.21 0.92 0.028 下降 甘油一脂 MG(a-13:0/0:0/0:0)[rac] C16H32O4 1.92 0.79 0.040 下降 反式-2-十二碳烯二酸 Traumatic acid C12H20O4 1.37 0.86 0.025 下降 9-氧壬酸 9-Oxo-nonanoic acid C9H16O3 1.15 0.93 0.019 下降 8-羟基-5,6-辛二烯酸 8-Hydroxy-5,6-octadienoic acid C8H12O3 1.04 0.94 0.016 下降 (S)-3-磺酸盐 (S)-3-Sulfonatolactate C3H6O6S 2.57 3.05 0.046 上升 γ-谷氨酰缬氨酸 gamma-Glutamylvaline C10H18N2O5 1.80 0.87 0.008 下降 穗花牡荆苷 Agnuside C22H26O11 1.41 0.91 0.005 下降 假尿苷 Pseudouridine C9H12N2O6 1.19 0.93 0.003 下降 雌三醇7-(6-反式-对-香豆酰基葡萄糖苷)
Eriodictyol 7-(6-trans-p-coumaroylglucoside)C30H28O13 1.44 0.92 0.020 下降 注:FC表示某差异代谢物在淡水组相对于海水组的表达倍数变化,FC>1表示该代谢物上调,FC<1表示该代谢物下调;P<0.05表示差异显著,P<0.01表示差异极显著。 Note: FC value indicates the change of expression multiple of a differential metabolite in FOV group compared with SOV group. FC>1 indicates that the metabolite is upregulated; FC<1 indicates that the metabolite is down-regulated. P<0.05 indicates significant difference, and P<0.01 indicates extremely significant difference. 表 2 差异代谢产物的KEGG通路富集表
Table 2 KEGG pathway enrichment table of differential metabolites
富集的差异代谢产物名称
Enriched differential metabolite name富集个数
Numble富集通路
ID Pathway通路描述
Pathway desciptionP 半乳糖神经酰胺 Galactosyl ceramide (d18:1/14:0) 1 map00600 鞘脂代谢 0.059 脯氨酸 L-Proline 1 map00330 精氨酸和脯氨酸代谢 0.176 异柠檬酸盐 Isocitrate 1 map00020 柠檬酸循环 0.057 戊二酸 Glutaric acid 1 map00310 赖氨酸降解 0.130 皮质醇 Cortisol 1 map00140 类固醇激素生物合成 0.207 尿苷 Uridine;假尿苷 Pseudo uridine 2 map00240 嘧啶代谢 0.011 异柠檬酸盐 Isocitrate 1 map00630 乙醛酸和二羧酸代谢 0.137 L-异亮氨酸 L-Isoleucine;脯氨酸 L-Proline 2 map00970 氨酰tRNA生物合成 0.008 L-异亮氨酸 L-Isoleucine 1 map00290 缬氨酸、亮氨酸和异亮氨酸生物合成 0.064 二十碳五烯酸 Eicosapentaenoic acid 1 map01040 不饱和脂肪酸生物合成 0.097 L-异亮氨酸 L-Isoleucine 1 map00280 缬氨酸、亮氨酸和异亮氨酸降解 0.107 (S)-3-磺酸盐 (S)-3-Sulfonatolactate 1 map00270 半膀氨酸和蛋氨酸代谢 0.141 -
[1] SU M, DUAN Z, SHI H, et al. The effects of salinity on reproductive development and egg and larvae survival in the spotted scat Scatophagus argus under controlled conditions[J]. Aquac Res, 2019, 50(7): 1782-1794. doi: 10.1111/are.14056
[2] MOHANTY B, GUPTA K, BABU G, et al. Exposure to salinity stress cause ovarian disruption in a stenohaline freshwater teleost, Heteropneustes fossilis (Bloch, 1794)[J]. Aquac Res, 2020, 51(5): 1964-1972. doi: 10.1111/are.14548
[3] PHAM H Q, KJØRSVIK E, NGUYEN A T, et al. Reproductive cycle in female Waigieu seaperch (Psammoperca waigiensis) reared under different salinity levels and the effects of dopamine antagonist on steroid hormone levels[J]. J Exp Mar Biol Ecol, 2010, 383(2): 137-145. doi: 10.1016/j.jembe.2009.12.010
[4] PHAM H Q, NGUYEN A T, KJØRSVIK E, et al. Seasonal reproductive cycle of Waigieu seaperch (Psammoperca waigiensis)[J]. Aquac Res, 2012, 43(6): 815-830. doi: 10.1111/j.1365-2109.2011.02894.x
[5] 牛景彦, 刘占才. 影响鱼类性腺发育的生态因素研究[J]. 农业与技术, 2016, 36(16): 109. [6] 杜学芳. 盐度对凡纳滨对虾繁殖及家系生长、存活的影响[D]. 上海: 上海海洋大学, 2013: 9-17. [7] 吴旭干, 赵亚婷, 何杰, 等. 低盐度海水和淡水对中华绒螯蟹性腺发育及交配行为的影响[J]. 动物学杂志, 2013, 48(4): 555-561. [8] 吴建辉, 王家启, 戴小杰, 等. 基于概率模型的长江口鱼类空间共现模式分析[J]. 南方水产科学, 2019, 15(1): 1-9. doi: 10.12131/20180112 [9] ZHU G L, WANG L J, TANG W Q, et al. Identification of olfactory receptor genes in the Japanese grenadier anchovy Coilia nasus[J]. Genes Genomics, 2017, 39(5): 521-532. doi: 10.1007/s13258-017-0517-8
[10] DUAN J R, ZHOU Y F, XU D P, et al. Ovary transcriptome profiling of Coilia nasus during spawning migration stages by Illumina sequencing[J]. Mar Genomics, 2015, 21: 17-19.
[11] LI W X, ZOU H, WU S G, et al. Richness and diversity of helminth communities in the JapaneseI grenadier anchovy, Coilia nasus, during its anadromous migration in the Yangtze River, China[J]. J Parasitol, 2012, 98(3): 449-452. doi: 10.1645/GE-2983.1
[12] LI W X, SONG R, WU S G, et al. Seasonal occurrence of helminths in the anadromous fish Coilia nasus (Engraulidae): parasite indecators of fish migratory movements[J]. J Parasitol, 2011, 97(2): 192-196. doi: 10.1645/GE-2621.1
[13] ZHANG H, WU G, XIE P, et al. Role of body size and temporal hydrology in the dietary shifts of shortjaw tapertail anchovy Coi-lia brachygnathus (Actinopterygii, Engraulidae) in a large floodplain lake[J]. Hydrobiologia, 2013, 703(1): 247-256. doi: 10.1007/s10750-012-1370-z
[14] 代培, 严燕, 朱孝彦, 等. 长江刀鲚国家级水产种质资源保护区(安庆段)刀鲚资源现状[J]. 中国水产科学, 2020, 27(11): 1267-1276. [15] SHEN H, GU R, XU G, et al. In-depth transcriptome analysis of Coilia ectenes, an important fish resource in the Yangtze River: de novo assembly, gene annotation[J]. Mar Genomics, 2015, 23: 15-17. doi: 10.1016/j.margen.2015.03.002
[16] 鲜博, 高建操, 徐钢春, 等. 盐度对刀鲚生长、抗氧化应激和渗透压调节能力的影响[J]. 海洋湖沼通报, 2020(2): 152-159. [17] 王武. 鱼类增养殖学[M]. 北京: 中国农业出版社, 2000: 193. [18] XU G, DU F, LI Y, et al. Integrated application of transcripto-mics and metabolomics yields insights into population-asynchronous ovary development in Coilia nasus[J]. Sci Rep, 2016, 6(1): 1-11. doi: 10.1038/s41598-016-0001-8
[19] YIN D, LIN D, YING C, et al. Metabolic mechanisms of Coilia nasus in the natural food intake state during migration[J]. Genomics, 2020, 112(5): 3294-3305. doi: 10.1016/j.ygeno.2020.05.027
[20] ZHAO H J, XU J K, YAN Z H, et al. Microplastics enhance the developmental toxicity of synthetic phenolic antioxidants by disturbing the thyroid function and metabolism in developing zebrafish[J]. Environ Int, 2020, 140: 105750. doi: 10.1016/j.envint.2020.105750
[21] DING J, HUANG Y, LIU S, et al. Toxicological effects of nano-and micro-polystyrene plastics on red tilapia: are larger plastic particles more harmless?[J]. J Hazard Mater, 2020, 396: 122693. doi: 10.1016/j.jhazmat.2020.122693
[22] ZHANG H, LIU Y, ZHOU L, et al. Metabonomic insights into the sperm activation mechanisms in ricefield eel (Monopterus albus)[J]. Genes, 2020, 11(11): 1259. doi: 10.3390/genes11111259
[23] HUANG Y, ZHANG Y, ZHENG J, et al. Metabolic profiles of fish nodavirus infection in vitro: RGNNV induced and exploited cellular fatty acid synthesis for virus infection[J]. Cell Microbiol, 2020, 22(9): e13216.
[24] OZEN G, NOREL X. Prostanoids in the pathophysiology of human coronary artery[J]. Prostag Oth Lipid M, 2017, 133: 20-28.
[25] YAN H, ZHANG M Z, WONG G, et al. Mechanisms of U46619-induced contraction in mouse intrarenal artery[J]. Clin Exp Pharmacol Physiol, 2019, 46(7): 643-651. doi: 10.1111/1440-1681.13087
[26] SMYTH E M. Thromboxane and the thromboxane receptor in cardiovascular disease[J]. Clin Lipidol, 2010, 5(2): 209-219. doi: 10.2217/clp.10.11
[27] GU W, MADRID D M D, YANG G, et al. Unaltered influenza disease outcomes in swine prophylactically treated with α-galactosylceramide[J]. Dev Comp Immunol, 2021, 114: 103843. doi: 10.1016/j.dci.2020.103843
[28] WANG T, LIU F, TIAN G, et al. Lineage/species-specific expansion of the Mx gene family in teleosts: differential expression and modulation of nine Mx genes in rainbow trout Oncorhynchus mykiss[J]. Fish Shellfish Immunol, 2019, 90: 413-430. doi: 10.1016/j.fsi.2019.04.303
[29] HOEHN K, MARIEB E N. Human anatomy & physiology[M]. San Francisco: Benjamin Cummings, 2010: 1064-1094.
[30] SOPINKA N M, CAPELLE P M, SEMENIUK C A D, et al. Glucocorticoids in fish eggs: variation, interactions with the environment, and the potential to shape offspring fitness[J]. Physiol Biochem Zool, 2017, 90(1): 15-33. doi: 10.1086/689994
[31] MOUSA M A, IBRAHIM M G, KORA M F, et al. Experimental studies on the reproduction of the thin-lipped mullet, Liza ramada[J]. Egypt J Aquat Biol Fish, 2018, 22(3): 125-138. doi: 10.21608/ejabf.2018.9455
[32] CUI W, MA A, WANG X, et al. Myo-inositol enhances the low-salinity tolerance of turbot (Scophthalmus maximus) by modula-ting cortisol synthesis[J]. Biochem Biophys Res Commun, 2020, 526(4): 913-919. doi: 10.1016/j.bbrc.2020.04.004
[33] 胡静, 叶乐, 吴开畅, 等. 急性盐度胁迫对克氏双锯鱼幼鱼血浆皮质醇浓度和 Na+-K+-ATP 酶活性的影响[J]. 南方水产科学, 2016, 12(2): 116-120. doi: 10.3969/j.issn.2095-0780.2016.02.017 [34] SARAVANAN M, RAMESH M, PETKAM R, et al. Influence of environmental salinity and cortisol pretreatment on gill Na+/K+-ATPase activity and survival and growth rates in Cyprinus carpio[J]. Aquac Rep, 2018, 11: 1-7. doi: 10.1016/j.aqrep.2018.04.002
[35] KWON N H, FOX P L, KIM S. Aminoacyl-tRNA synthetases as therapeutic targets[J]. Nat Rev Drug Discov, 2019, 18(8): 629-650. doi: 10.1038/s41573-019-0026-3
[36] GOMEZ M A R, IBBA M. Aminoacyl-tRNA synthetases[J]. RNA, 2020, 26(8): 910-936. doi: 10.1261/rna.071720.119
[37] GARAVITO M F, NARVÁEZ-ORTIZ H Y, ZIMMERMANN B H. Pyrimidine metabolism: dynamic and versatile pathways in pathogens and cellular development[J]. J Genet Genomics, 2015, 42(5): 195-205. doi: 10.1016/j.jgg.2015.04.004
[38] LAVIEU G, SCARLATTI F, SALA G, et al. Regulation of autophagy by sphingosine kinase 1 and its role in cell survival during nutrient starvation[J]. J Biol Chem, 2006, 281(13): 8518-8527. doi: 10.1074/jbc.M506182200
[39] HANNUN Y A, OBEID L M. The ceramide-centric universe of lipid-mediated cell regulation: stress encounters of the lipid kind[J]. J Biol Chem, 2002, 277(29): 25847-25850. doi: 10.1074/jbc.R200008200
-
期刊类型引用(4)
1. 王美凤,王秀华,李晨,吕若萱,杨冰. 水生动物环境DNA研究概况及其在水生动物疫病传播研究中的应用前景. 中国动物检疫. 2024(01): 87-94 . 百度学术
2. 刘洋,应方,杨俊,王焕英,洪文杰,王庭璋. 鱼类环境DNA保存条件的比较研究. 生物学杂志. 2024(03): 114-120 . 百度学术
3. 郑若丹,陈炳耀,张帅,麦俊晓,蒋佩文,王文欣,李敏. 基于环境DNA的涠洲岛周围海域布氏鲸种群分布初探. 南方水产科学. 2024(03): 36-46 . 本站查看
4. 陈炳耀,信誉,路方婷,孙婧,刘杉,李敏,吴斌,王崇瑞,张衡,周玉香,白玲. 长江江豚监测现状及展望. 中国环境监测. 2023(02): 1-10 . 百度学术
其他类型引用(2)