廖書丹 譚艾娟 孫海峰 向潤 楊劍 張元鑫 宋旭琴 呂世明

摘要：【目的】探究環(huán)境脅迫因子對豬源多重耐藥大腸桿菌生存適應(yīng)性的影響，為創(chuàng)造多重耐藥菌株帶來高適應(yīng)代價的環(huán)境及加速降低種群耐藥水平提供理論依據(jù)?！痉椒ā恳苑蛛x自貴州省規(guī)模化豬場的多重耐藥大腸桿菌QL15為研究對象，以大腸桿菌質(zhì)控（敏感）菌株ATCC25922為對照，經(jīng)體外競爭試驗測定不同環(huán)境（培養(yǎng)基濃度、碳源、氮源、pH和溫度）脅迫下耐藥菌株QL15相對于敏感菌株ATCC25922的適應(yīng)性，并通過重測序變異位點和實時熒光定量PCR測定分析耐藥基因表型及外排泵和外膜孔道蛋白基因表達與細菌適應(yīng)性代價的關(guān)系。【結(jié)果】耐藥菌株QL15在體外獨立及混合培養(yǎng)時，除在稀釋LB培養(yǎng)基（1/4、1/8和1/16）中的長勢較敏感菌株ATCC25922低外，在低碳源（0.1%和0.5%）、低氮源（0.01%和0.05%）、酸堿性（pH=5.0、pH=6.0和pH=8.0）及低溫（20 ℃）和高溫（42 ℃）脅迫下的生長均優(yōu)于敏感菌株ATCC25922，且混合培養(yǎng)結(jié)果優(yōu)于獨立培養(yǎng)。將耐藥菌株QL15基因組序列與敏感菌株ATCC25922基因組序列進行比對，共檢測到509個InDel突變位點，其中Delete突變位點233個、Insert突變位點276個;共注釋到71個功能基因，以細胞壁與細胞膜相關(guān)基因突變最多（10個），占14.1%。與敏感菌株ATCC25922相比，耐藥菌株QL15在不同環(huán)境脅迫下其外排泵與外膜孔道蛋白相關(guān)基因emrB、acrB、tolC和ompC呈下調(diào)表達，而ompF基因呈上調(diào)表達。【結(jié)論】多重耐藥大腸桿菌耐藥性的獲得促使其具有更強的適應(yīng)性，但對可直接利用營養(yǎng)物質(zhì)的獲取弱于敏感菌株，即處于競爭劣勢。耐藥大腸桿菌的多重耐藥性主要是由多種外排泵和外膜孔道蛋白介導(dǎo)完成，其中，ompF基因在pH脅迫、全營養(yǎng)脅迫和高、低溫脅迫下均處于上調(diào)表達狀態(tài)，故推測ompF基因的高表達有利于大腸桿菌在上述脅迫條件下生長。

關(guān)鍵詞： 多重耐藥大腸桿菌;豬源;環(huán)境脅迫因子;相對適應(yīng)性;耐藥性

中圖分類號： S852.612? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2022）04-0946-12

Effects of different environmental stress factors on the growth of pig-origin drug-resistant Escherichia coli

LIAO Shu-dan1，TAN Ai-juan2，SUN Hai-feng2，XIANG Run2，YANG Jian1，

ZHANG Yuan-xin3，SONG Xu-qin1*，LYU Shi-ming1*

（1College of Animal Science，Guizhou University， Guiyang， Guizhou? 550025， China; 2College of Life Sciences，Guizhou University， Guiyang， Guizhou? 550025， China; 3 Department of Agriculture and Rural Affairs of

Guizhou Province， Guiyang， Guizhou? 550005， China）

Abstract：【Objective】To investigate the effects of environmental stress factors on the survival adaptability of pig-origin multi-drug resistant Escherichia coli， so as to provide a theoretical basis for creating an environment where multi-drug resistant strains could present a high cost of adaptability and for accelerating the reduction of resistance level in the population. 【Method】Multi-drug resistant E. coli QL15 isolated from a large-scale pig farm in Guizhou Province was taken as the study object， and E. coli quality control （sensitive） strain ATCC25922 was taken as the control. The adapta-bility of drug-resistant strain QL15 and that of sensitive strain ATCC25922 under different environmental stress （medium concentration， carbon source， nitrogen source， pH and temperature） was determined by in vitro competition assay. The relationship between drug-resistant gene phenotype and gene expression of efflux pump and outer membrane pore protein with the cost of bacterial adaptation was analyzed by resequencing variant loci and fluorescence quantitative PCR. 【Result】Except that the growth of drug-resistant strain QL15 was lower in dilute LB medium （1/4， 1/8 and 1/16） compared with the sensitive strain ATCC25922， the growth of drug-resistant strain QL15 was better than that of sensitive strain ATCC25922 in in vitro independent and mixed cultures under low carbon sources （0.1% and 0.5%）， low nitrogen sources （0.01% and 0.05%）， acid-base property （pH=5.0， pH=6.0 and pH=8.0）， low and high temperature （20 ℃ and 42 ℃）. Besides， the results of mixed culture were better than that of the independent culture. The genome sequence of drug-resistant strain QL15 was compared with that of sensitive strain ATCC25922， and a total of 509 InDel mutation sites were detected， including 233 delete mutation sites and 276 insert mutation sites; a total of 71 functional genes were annotated， most of which were mutations related to cell wall and cell membrane（10）， accounting for 14.1%. Compared with the sensitive strain ATCC25922， in the drug-resistant strain QL15， genes related to efflux pump and outer membrane pore protein， namely emrB， acrB， tolC and ompC showed down-regulated expression， and ompF gene showed up-regulated expression under different environmental stress conditions. 【Conclusion】The acquisition of drug resistance of multi-drug-resistant E. coli makes it more adaptive， while its acquisition of directly available nutrients is weaker than that of sensitive strains， that is， E. coli is at a competitive disadvantage. The multi-drug resistance of drug-resistant E. coli is mainly mediated by a variety of efflux pumps and outer membrane pore proteins. Among them， the ompF gene shows up-regulated expression under pH stress， allotrophic stress， high and low temperature stress， so it is presumed that the high expression of ompF gene is beneficial to the growth of E. coli under the above stress conditions.

Key words：multi-drug-resistant Escherichia coli; pig-origin; environmental stress factors; relative adaptability;drug resistance

Foundation items： National Natural Science Foundation of China （31760749）;Guizhou Science and Technology Project （QKHPTRC〔2017〕5788，QKHJC-ZK〔2022〕yiban129）;Guizhou University Natural Science Special Items （Special Post） Scientific Research Fund Project （GDTGHZ〔2020〕25）

0 引言

【研究意義】細菌性疾病嚴重威脅著畜牧養(yǎng)殖業(yè)的健康發(fā)展（劉菊枚等，2020;付霞麗等，2021），雖然抗菌藥物在動物疾病防治、促進動物生長、保障動物福利及提高養(yǎng)殖經(jīng)濟效益等方面發(fā)揮著重要作用（Bassetti et al.，2017;李瑋瑋等，2021），但因動物養(yǎng)殖過程中不合理使用及濫用抗菌藥物造成的細菌耐藥性問題日益嚴重，不僅直接影響動物性食品安全，還可通過食物鏈對人類健康構(gòu)成威脅，解決細菌耐藥問題已迫在眉睫（周瑤琴等，2018;劉昌孝，2019;劉躍華等，2019）。因此，在動物養(yǎng)殖“減抗、禁抗”背景下（Ghosh et al.，2019;顧曉曉等，2020），開展不同環(huán)境條件下耐藥菌的生長適應(yīng)性及相關(guān)耐藥基因表達調(diào)控研究，對尋找無抗菌藥物壓力下加速細菌敏感性恢復(fù)具有重要意義?！厩叭搜芯窟M展】減少抗菌藥物使用（減抗）是目前遏制細菌耐藥性發(fā)展的重要手段之一（袁萬哲等，2019），其主要理論依據(jù)是耐藥細菌在無藥環(huán)境中其適應(yīng)性代價高于敏感菌（Andersson and Levin，1999）。Gagneux等（2006）的研究結(jié)果顯示，耐利福平結(jié)核分枝桿菌產(chǎn)生了高額的適應(yīng)性代價;Vogwill等（2015）經(jīng)Meta分析發(fā)現(xiàn)，經(jīng)藥物壓力細菌耐藥性的產(chǎn)生付出了較大代價，其中染色體突變付出的適應(yīng)代價通常高于質(zhì)粒的獲得，同時表現(xiàn)出耐藥菌株生長趨勢弱于敏感菌株;de Sousa等（2017）對克隆出的耐利福平和鏈霉素大腸桿菌與敏感菌進行競爭試驗，發(fā)現(xiàn)在無藥物作用條件下克隆耐藥菌株競爭處于弱勢，但該現(xiàn)象可通過補償性進化得以改變。因此，細菌在抗菌藥物的選擇壓力下產(chǎn)生耐藥性的相關(guān)變化，會改變細菌體內(nèi)生理生化關(guān)鍵酶等結(jié)構(gòu)，而對細菌的生長造成額外負擔，進而影響細菌的部分功能，增加其適應(yīng)性代價（楊焰等，2019）。耐藥菌在停藥后的無藥環(huán)境中會因其適應(yīng)性代價高于敏感菌而處于競爭劣勢，逐漸被敏感菌群取代，從而降低種群的耐藥水平。此外，耐藥菌的適應(yīng)性代價受環(huán)境差異影響，各基因間存在復(fù)雜的協(xié)同與拮抗機制，同一細菌在不同環(huán)境中的適應(yīng)性代價也不同（方結(jié)紅，2018）。Petersen等（2009）研究表明，在環(huán)境脅迫下鏈霉素耐藥大腸桿菌適應(yīng)性代價顯著增加;商可心（2015）研究發(fā)現(xiàn)，鼠傷寒沙門氏菌rpsL突變介導(dǎo)的鏈霉素耐藥菌在營養(yǎng)豐富時的生長弱于敏感菌，但在碳源缺乏時耐藥菌因RpoS蛋白表達水平低，其生長優(yōu)于敏感菌;Maharjan和Ferenci（2017）研究表明，rpoB突變的利福平耐藥大腸桿菌在營養(yǎng)豐富時適應(yīng)性代價較高，而在營養(yǎng)缺乏時其適應(yīng)性更強;孫海峰等（2021）研究發(fā)現(xiàn)，在鹽脅迫下多重耐藥大腸桿菌的生長優(yōu)于敏感菌，其機制可能與外膜孔道蛋白基因表達上調(diào)有關(guān)?！颈狙芯壳腥朦c】研究耐藥菌的適應(yīng)性代價有助于認識干預(yù)耐藥軌跡、新藥研發(fā)及聯(lián)合用藥，提高抗生素敏感性和抑制耐藥性進化，延長抗生素使用年限。適應(yīng)性代價是細菌耐藥性進化的重要環(huán)節(jié)，而在無抗菌藥物壓力下環(huán)境脅迫是細菌在自然環(huán)境中常遇的生存條件，但至今有關(guān)適應(yīng)性代價在不同環(huán)境中對耐藥菌的影響及其作用機制鮮見研究報道?！緮M解決的關(guān)鍵問題】以分離自貴州省規(guī)?；i場的多重耐藥大腸桿菌QL15為研究對象，以大腸桿菌質(zhì)控（敏感）菌株ATCC25922為對照，經(jīng)體外競爭試驗測定不同環(huán)境脅迫下耐藥菌株QL15相對于敏感菌株ATCC25922的適應(yīng)性，并檢測其攜帶的耐藥基因，分析耐藥基因表型及外排泵和外膜孔道蛋白基因表達與細菌適應(yīng)性代價的關(guān)系，旨在探究環(huán)境脅迫因子對多重耐藥菌株生存適應(yīng)性的影響，為創(chuàng)造多重耐藥菌株帶來高適應(yīng)代價的環(huán)境及加速降低種群耐藥水平提供理論依據(jù)。

1 材料與方法

1. 1 試驗材料

大腸桿菌耐藥菌株QL15由貴州大學(xué)動物科學(xué)學(xué)院獸醫(yī)藥理實驗室保存提供，敏感菌株ATCC25922購自中國獸醫(yī)藥品監(jiān)察所。LB肉湯培養(yǎng)基、MH培養(yǎng)基及平板計數(shù)瓊脂培養(yǎng)基購自貴州為萊科技有限責任公司;瓊脂糖購自美國Sigma公司;細菌DNA提取試劑盒、RNA提取試劑盒、PCR相關(guān)試劑及RT-PCR相關(guān)試劑購自北京索萊寶科技有限公司。

1. 2 培養(yǎng)基制備

碳源和氮源脅迫培養(yǎng)基制備：參照M9基礎(chǔ)培養(yǎng)基配方（王希越等，2017）進行制備，其中，碳源脅迫培養(yǎng)基為0.1%和0.5%葡萄糖，氮源脅迫培養(yǎng)基為0.01%和0.05%氯化銨，其他組份含量不變，高壓滅菌后使用。酸堿脅迫培養(yǎng)基制備：常規(guī)LB培養(yǎng)基pH為7.0，分別以1 mol/L HCl和1 mol/L NaOH調(diào)節(jié)LB培養(yǎng)基pH至5.0、6.0、8.0和9.0，高壓滅菌后使用。全營養(yǎng)脅迫培養(yǎng)基：LB培養(yǎng)基與0.9%生理鹽水比例為1∶3、1∶7和1∶15，對應(yīng)的Luria-Bertani（LB）濃度分別為1/4、1/8和1/16，高壓滅菌后使用。溫度脅迫培養(yǎng)基：高溫和低溫脅迫條件選用LB培養(yǎng)基。

1. 3 不同環(huán)境脅迫下大腸桿菌生長曲線測定

采用比濁法測定敏感和耐藥大腸桿菌在不同環(huán)境因子脅迫下的生長曲線（姜興粲等，2020）。設(shè)活化敏感與耐藥大腸桿菌的OD600 nm=0.5，吸取50 μL菌液置于50 mL培養(yǎng)基中進行恒溫振蕩培養(yǎng)。碳源脅迫、氮源脅迫、酸堿脅迫和全營養(yǎng)脅迫的培養(yǎng)溫度均為37 ℃，轉(zhuǎn)速為150 r/min，溫度脅迫條件為42和20 ℃，每隔1 h測定其吸光值。

1. 4 不同環(huán)境脅迫下耐藥大腸桿菌相對適應(yīng)性測定

相對適應(yīng)性表現(xiàn)為不同環(huán)境脅迫下耐藥菌株和敏感菌株生存競爭能力的差異，主要通過體外混合培養(yǎng)競爭試驗測定。參考本課題組前期改良的方法（孫海峰等，2021）：分別取10 μL的敏感菌株ATCC25922和耐藥菌株QL15（活化后調(diào)節(jié)麥氏比濁度至0.5），置于5 mL不同脅迫條件的培養(yǎng)基中，在37 ℃、150 r/min的條件下進行培養(yǎng)（高溫脅迫培養(yǎng)溫度42 ℃，低溫脅迫培養(yǎng)溫度20 ℃）。取培養(yǎng)0和24 h的菌液，適當稀釋后涂布在無藥培養(yǎng)基和有藥（含5 μg/mL恩諾沙星）培養(yǎng)基上。無藥培養(yǎng)基的菌落數(shù)為敏感菌落和耐藥菌落總數(shù)，有藥培養(yǎng)基則為耐藥菌落數(shù)，二者的差值即為敏感菌落數(shù)?；旌吓囵B(yǎng)0 h的敏感菌株和耐藥菌株數(shù)量記錄為S0和R0，培養(yǎng)24 h的分別記錄為R24和S24，根據(jù)以下公式計算相對適應(yīng)度（Relative fitness，F(xiàn)）：

F=ln（R24/R0）/ln（S24/S0）

1. 5 菌株外排泵與外膜孔道蛋白基因表達量測定

采用本課題組設(shè)計的引物（孫海峰等，2021），委托生工生物工程（上海）股份有限公司合成。敏感菌株和耐藥菌株總RNA按照試劑盒操作說明進行提取。菌株外排泵與外膜蛋白基因表達量測定：將耐藥菌株QL15和敏感菌株ATCC25922在37 ℃、150 r/min的條件下活化培養(yǎng)12 h，取50 μL菌液加入到50 mL不同環(huán)境脅迫培養(yǎng)基中，待細菌生長至對數(shù)生長期時，提取細菌總RNA并反轉(zhuǎn)錄合成cDNA，采用實時熒光定量PCR測定耐藥菌株和敏感菌株外排泵與外膜孔道蛋白基因在不同環(huán)境脅迫下的相對表達量。實時熒光定量PCR反應(yīng)體系及擴增程序參照孫海峰等（2021）的方法。

2 結(jié)果與分析

2. 1 不同環(huán)境脅迫下大腸桿菌的生長曲線

2. 1. 1 碳源和氮源脅迫下大腸桿菌的生長曲線

如圖1所示，在不同濃度的碳源和氮源脅迫下，耐藥菌株QL15與敏感菌株ATCC25922相比其生長速率更快，能更快地進入對數(shù)生長期且峰濃度均高于敏感菌株ATCC25922，表明碳源和氮源的濃度變化對耐藥菌株QL15生長速率影響較小，其獨立生長適應(yīng)性更強。隨著碳源濃度的升高，耐藥菌株QL15和敏感菌株ATCC25922的生長速率均降低，而氮源濃度的升高能促進二者快速生長（圖2）。此外，由圖3可看出，培養(yǎng)10 h后不同碳源和氮源對敏感菌株ATCC25922的影響均達極顯著水平（P<0.01，下同），當碳源為0.1%和氮源為0.01%時對耐藥菌株QL15的影響極顯著，當?shù)礊?.05%時對耐藥菌株QL15的影響顯著（P<0.05，下同），但在碳源為0.5%時對耐藥菌株QL15生長無顯著影響（P>0.05，下同）。同時，碳源濃度為0.5%時耐藥菌株QL15的峰濃度是0.1%碳源濃度的2倍，繼續(xù)增加濃度至1.0%則對峰值影響不明顯，表明在一定范圍內(nèi)碳源對耐藥大腸桿菌種群的承載度影響明顯。

2. 1. 2 pH脅迫下大腸桿菌的生長曲線 不同pH對敏感菌株ATCC25922和耐藥菌株QL15生長的影響如圖4和圖5所示。與敏感菌株ATCC25922相比，耐藥菌株QL15在pH=5.0、pH=6.0和pH=8.0脅迫下均能較快地進入對數(shù)生長期，且其峰值濃度差異不顯著;在pH=9.0時，耐藥菌株QL15進入對數(shù)生長期的時間更快，但在對數(shù)生長期的生長速率慢于敏感菌株ATCC25922。pH=7.0是細菌生長的最適pH，pH的升高和降低均會抑制敏感菌株ATCC25922和耐藥菌株QL15的生長，尤其在pH=9.0脅迫下耐藥菌株QL15的生長受到明顯抑制。

2. 1. 3 全營養(yǎng)脅迫下大腸桿菌的生長曲線 不同LB濃度對敏感菌株ATCC25922和耐藥菌株QL15生長的影響如圖6和圖7所示。不同LB濃度脅迫下，敏感菌株ATCC25922和耐藥菌株QL15進入對數(shù)生長期的時間及對數(shù)生長期的生長速率差異均不顯著，但敏感菌株ATCC25922在穩(wěn)定期的峰值濃度均顯著高于耐藥菌株QL15，表明敏感菌株在不同營養(yǎng)物質(zhì)濃度下處于優(yōu)勢生長。此外，隨著LB培養(yǎng)基稀釋倍數(shù)的增加，敏感菌株ATCC25922和耐藥菌株QL15對數(shù)生長期的生長速率和穩(wěn)定期的峰值濃度均明顯降低，表明LB濃度脅迫對細菌的生長具有明顯影響。

2. 1. 4 不同溫度脅迫下大腸桿菌的生長曲線 不同溫度脅迫對敏感菌株ATCC25922和耐藥菌株QL15生長的影響如圖8所示。當溫度為20 ℃時，敏感菌株ATCC25922和耐藥菌株QL15的適應(yīng)期時間相近且二者在培養(yǎng)10 h內(nèi)均未到達穩(wěn)定期，但耐藥菌株QL15在對數(shù)生長期的生長速率更快，為優(yōu)勢菌株;當培養(yǎng)溫度為42 ℃時，敏感菌株ATCC25922和耐藥菌株QL15的生長曲線差異不明顯，且到達穩(wěn)定期的峰值濃度相近，但耐藥菌株QL15在對數(shù)生長期的生長速率仍高于敏感菌株ATCC25922。比較敏感菌株ATCC25922和耐藥菌株QL15在不同溫度脅迫下的生長曲線，結(jié)果（圖9）發(fā)現(xiàn)37 ℃下大腸桿菌的生長速率最高，42和20 ℃對細菌生長均有一定的抑制作用，且以20 ℃的抑制作用最大。

2. 2 耐藥大腸桿菌在不同環(huán)境脅迫下的相對適應(yīng)性

耐藥菌株QL15的相對適應(yīng)性如表1所示。在LB、碳源、氮源、pH=6.0～8.0、20 ℃和42 ℃脅迫下，耐藥菌株QL15的相對適應(yīng)度均大于1.00，表明其與敏感菌株ATCC25922在上述環(huán)境脅迫下競爭優(yōu)勢明顯;耐藥菌株QL15在pH=9.0及所有LB培養(yǎng)基稀釋條件下則表現(xiàn)出競爭劣勢（相對適應(yīng)度<1.00），其中1/4 LB濃度脅迫下的相對適應(yīng)度最低，僅為0.95。

2. 3 耐藥大腸桿菌基因突變位點檢測結(jié)果

將耐藥菌株QL15和敏感菌株ATCC25922的基因組序列進行比對，共檢測到509個InDel突變位點，其中Delete突變位點233個、Insert突變位點276個;共注釋到71個功能基因，以細胞壁與細胞膜相關(guān)基因突變最多（10個），占14.1%，其次為氨基酸的轉(zhuǎn)運與代謝（9個）占12.7%。突變涉及多種蛋白，包括磷酸轉(zhuǎn)移酶、羧肽酶、細胞分裂蛋白、青霉素結(jié)合蛋白、ABC型外排系統(tǒng)、細胞壁水解酶、甘油醛-3-磷酸脫氫酶及外膜蛋白TolC等均檢測到InDel位點。

2. 4 不同環(huán)境脅迫因子對大腸桿菌外排泵與外膜孔道蛋白基因表達量的影響

2. 4. 1 碳源和氮源脅迫下大腸桿菌外排泵與外膜孔道蛋白基因的相對表達量 在碳源和氮源脅迫下，大腸桿菌外排泵與外膜孔道蛋白基因的相對表達情況如圖10所示。與常規(guī)的LB培養(yǎng)條件相比，敏感菌株ATCC25922和耐藥菌株QL15中的外排泵與外膜孔道蛋白基因（acrB、ompF、emrB、tolC和ompC）相對表達量均呈下調(diào)趨勢，且耐藥菌株QL15的5種基因在0.5% C、0.05% N、0.01% N和1.0% C+0.1% N條件下的相對表達量較0.1% C脅迫的下調(diào)程度更大，敏感菌株ATCC25922則表現(xiàn)為除0.1% C和1.0% C+0.10% N條件外，其他條件下均呈顯著下調(diào)。

2. 4. 2 pH脅迫下大腸桿菌外排泵與外膜孔道蛋白基因的相對表達量 敏感菌株ATCC25922和耐藥菌株QL15在pH=5.0～9.0的條件下其外排泵與外膜孔道蛋白基因的相對表達量如圖11所示。與pH=7.0相比，敏感菌株ATCC25922和耐藥菌株QL15的emrB、acrB、tolC和ompC基因相對表達呈下調(diào)趨勢，而ompF基因表達呈上調(diào)趨勢。其中，ompC基因隨環(huán)境pH的變化未表現(xiàn)出明顯的變化趨勢（除耐藥菌株QL15在pH=9.0條件下的相對表達量外）。

2. 4. 3 全營養(yǎng)脅迫下大腸桿菌外排泵與外膜孔道蛋白基因的相對表達量 敏感菌株ATCC25922和耐藥菌株QL15的外排泵與外膜孔道蛋白基因表達在不同LB濃度脅迫下表現(xiàn)出不同變化規(guī)律。隨著LB稀釋倍數(shù)的增加，敏感菌株ATCC25922的emrB基因表達呈上調(diào)趨勢，而耐藥菌株QL15的emrB基因表達呈下調(diào)趨勢。2種大腸桿菌ompF基因在不同LB稀釋倍數(shù)下均呈上調(diào)表達，且耐藥菌株QL15的ompF基因隨稀釋倍數(shù)的增加其相對表達量呈遞增趨勢。此外，敏感菌株ATCC25922和耐藥菌QL15的acrB、tolC、ompC基因相對表達量在不同LB濃度脅迫下均呈下調(diào)趨勢（圖12）。

2. 4. 4 高溫和低溫脅迫下大腸桿菌外排泵與外膜孔道蛋白基因的相對表達量 敏感菌株ATCC25922和耐藥菌株QL15的外排泵與外膜孔道蛋白基因在高溫和低溫脅迫下的相對表達量如圖13所示。敏感菌株ATCC25922和耐藥菌QL15的emrB、acrB、tolC和ompC基因在20 ℃下的相對表達量較37 ℃下呈下調(diào)趨勢，在42 ℃下這4種基因的相對表達量則是上調(diào)和下調(diào)趨勢并存。2種大腸桿菌的ompF基因在20和42 ℃下均呈上調(diào)表達，且在42 ℃下的相對表達量高于20 ℃。

3 討論

據(jù)統(tǒng)計，我國每年約有9.7萬t抗菌藥物用于畜牧業(yè)，接近抗菌藥物年產(chǎn)總量的50%。由于抗菌藥物的大量使用及濫用，細菌耐藥性已成為全球廣泛關(guān)注的公共衛(wèi)生問題，意味著人類可能進入后抗生素時代（李瑞珍等，2011;Houston，2011）。目前，全球都在為降低細菌耐藥性而努力，主要措施包括從源頭控制抗菌藥物使用、阻斷耐藥病原傳播及開發(fā)新型藥物或疫苗等。通過阻斷耐藥性傳播和開發(fā)新型藥物等方法以降低細菌的耐藥性，其技術(shù)要求較高、研發(fā)較難、周期較長，且細菌耐藥速度遠大于新藥的產(chǎn)出速度，因此“減抗”被世界衛(wèi)生組織（WHO）列為主要措施。適應(yīng)性代價是降低細菌耐藥性的重要隱含機制，耐藥菌株在無藥物作用條件下較敏感菌株競爭劣勢，二者相互競爭生長后敏感菌株會成為優(yōu)勢菌群，其群體敏感性得以恢復(fù)。但在無藥物作用的環(huán)境中，耐藥基因的敏感性恢復(fù)速度極慢，因此，尋找無藥條件下加速耐藥菌株敏感性恢復(fù)的方法對降低或減緩細菌耐藥性具有重要意義。

本課題組前期對耐藥菌株QL15的耐藥基因進行檢測，結(jié)果發(fā)現(xiàn)存在15種耐藥基因（孫海峰等，2021），分別是aac（3）-IV、cmlA、ompC、ompF、acrA、tolC、emrB、emrD、tetA、blaTEM、acrB、sugE、mdtL、sul2和sul3。其中，與外排泵相關(guān)的基因有8種，外膜孔道蛋白基因有3種。AcrAB-TolC為RND外排家族，包括由融合蛋白（AcrA）、細胞質(zhì)膜轉(zhuǎn)運蛋白（AcrB）及外膜通道（TolC），是目前發(fā)現(xiàn)與大腸桿菌聯(lián)系最緊密的外排泵（Yasufuku et al.，2010），能外排氨基糖苷類、β-內(nèi)酰胺類、氯霉素類和大環(huán)內(nèi)酯類等多種抗菌藥物（黃嬌和穆青，2019）。在自然環(huán)境中，AcrAB-TolC的正常表達有助于大腸桿菌抵抗腸道膽鹽和脂肪酸等物質(zhì)的侵入，但在某些條件下AcrAB-TolC過量表達可使大腸桿菌表現(xiàn)出對多種化合藥物的耐受，甚至呈高水平耐受（甘露等，2020）?？梢姡馀疟煤屯饽た椎赖鞍讓Υ竽c桿菌耐藥性具有積極作用。

生長曲線（獨立生長）和相對適應(yīng)性（競爭生長）是測定細菌適應(yīng)性的主要方法（楊焰等，2019）。本研究中，耐藥菌株QL15在不同環(huán)境脅迫下獨立生長時，表現(xiàn)為在0.1% C、0.5% C、1.0% C+0.10% N、0.05% N、0.01% N、pH=5.0、pH=6.0、pH=8.0、20 ℃和42 ℃中的生長狀況均優(yōu)于敏感菌株ATCC25922，而在pH=9.0、1/4 LB、1/8 LB和1/16 LB中的生長狀況弱于敏感菌株ATCC25922。將2種大腸桿菌混合培養(yǎng)測定出的相對適應(yīng)性，表現(xiàn)為耐藥菌株QL15在0.1% C、0.5% C、1.0% C+0.10% N、0.05% N、0.01% N、pH=5.0、pH=6.0、pH=8.0、20 ℃和42 ℃中均處于競爭優(yōu)勢，而在pH=9.0、1/4 LB、1/8 LB和1/16 LB中處于劣勢生長。此外，本研究發(fā)現(xiàn)稀釋LB培養(yǎng)基不利耐藥菌株QL15生長，但隨著LB稀釋倍數(shù)的增加，其相對適應(yīng)性又逐漸增強，與林文芳等（2016）的研究結(jié)果相近，即貧營養(yǎng)和痕量抗生素對質(zhì)粒抗生素抗性適應(yīng)度代價有重要影響，較低的營養(yǎng)濃度和痕量抗生素均能增加耐藥菌株的適應(yīng)性代價。但某些耐藥菌株的基因突變在不同營養(yǎng)中對其生長影響可能不同，如耐利福平大腸桿菌的rpoB基因突變在營養(yǎng)豐富時能抑制其生長，而在營養(yǎng)豐富時（鐵元素、碳源和磷酸鹽）具有積極作用（Maharjan and Ferenci，2017）。本研究通過比對分析耐藥菌株QL15和敏感菌株ATCC25922的基因組序列，結(jié)果共檢測到509個InDel突變位點，且以細胞壁和細胞膜相關(guān)功能基因居多，占注釋到71個功能基因的14.1%。由此推測，耐藥菌株QL15在物質(zhì)轉(zhuǎn)運與代謝方面得到突變增強的同時，膜相關(guān)功能基因的突變可能促使其擁有更強的膜功能，包括通透性和物質(zhì)運輸能力等方面，致使耐藥菌株QL15在多數(shù)環(huán)境脅迫下較敏感菌株ATCC25922具有更強的適應(yīng)性。

腸道細菌的存活和生長取決于pH耐受性，pH可影響RND泵底物的積累和流出，細菌在酸性條件下的培養(yǎng)初期，acrA和acrB基因表達下調(diào)，培養(yǎng)18 h后檢測到基因表達顯著上調(diào)（Nové et al.，2020）。耐藥菌株QL15和敏感菌株ATCC25922在pH=5.0和pH=6.0脅迫下進入對數(shù)生長期的培養(yǎng)時間遠低于18 h，且二者的acrB基因相對表達量均呈下調(diào)趨勢，與Nové等（2020）的研究結(jié)果一致。此外，tolC和emrB基因的上調(diào)表達能促進大腸桿菌在極端酸性條件下存活（Deininger et al.，2011），但本研究中tolC和emrB基因的相對表達量并未上調(diào)，可能是由于pH=5.0尚未達到極端酸性條件，大腸桿菌或通過其他機制應(yīng)對酸脅迫。溫度是病原體在感染初期經(jīng)歷的關(guān)鍵環(huán)境信號。有研究表明，升高培養(yǎng)溫度能增加RND泵的表達量，如熒光假單胞菌RND泵（Adebusuyi and Foght，2011）和鼠傷寒沙門氏菌AcrAB泵（Hartog et al.，2008）在35和37 ℃下的表達量均高于25和28 ℃。本研究中，acrB和tolC基因在20 ℃脅迫下呈下調(diào)表達趨勢也進一步驗證了該結(jié)論。

細菌能通過抗菌藥物的非特異性通道蛋白OmpC和OmpF進行物質(zhì)交換，既能阻礙抗菌藥物進入細胞，又能促進營養(yǎng)吸收和有毒物質(zhì)的擴散（Knopp and Andersson，2015）。如大腸桿菌在葡萄糖中的生長速度取決于其通過中央碳代謝途徑的有效運輸及分解代謝，憑借外膜孔蛋白（OmpC、OmpF和LamB）將葡萄糖從細胞外環(huán)境運輸至周質(zhì)空間（Ferenci，2001;Masi et al.，2019）。根據(jù)細胞生長的環(huán)境條件差異，糖轉(zhuǎn)運系統(tǒng)的典型作用明顯不同。特定膜外蛋白（OMPs）在葡萄糖轉(zhuǎn)運中的作用因培養(yǎng)液中可用葡萄糖含量（痕量或過量）而異（Alva et al.，2020）。本研究中，耐藥菌株QL15和敏感菌株ATCC25922的ompC基因在不同LB濃度（1/4 LB、1/8 LB和1/16 LB）脅迫下均呈下調(diào)表達，ompF基因則呈上調(diào)表達，說明在此條件下糖類營養(yǎng)物質(zhì)可能主要通過OmpF蛋白進人細胞。但在不同碳源（0.1% C、0.5% C和1.0% C）脅迫下，2種大腸桿菌的ompC和ompF基因表達量均下調(diào)，可能是在此濃度脅迫下葡萄糖通過其他轉(zhuǎn)運體系進入細胞?？梢?，ompF基因?qū)毦鼓嫔L具有一定的積極作用。本研究中，ompF基因在pH脅迫、全營養(yǎng)脅迫和高、低溫脅迫下均呈上調(diào)表達，故推測ompF基因的高表達有利于大腸桿菌在上述脅迫條件下生長。

4 結(jié)論

多重耐藥大腸桿菌耐藥性的獲得促使其具有更強的適應(yīng)性，但對可直接利用營養(yǎng)物質(zhì)的獲取弱于敏感菌株，即處于競爭劣勢。耐藥大腸桿菌的多重耐藥性主要是由多種外排泵和外膜孔道蛋白介導(dǎo)完成，其中，ompF基因在pH脅迫、全營養(yǎng)脅迫和高、低溫脅迫下均處于上調(diào)表達狀態(tài)，故推測ompF基因的高表達有利于大腸桿菌在上述脅迫條件下生長。

參考文獻：

方結(jié)紅. 2018. 生豬屠宰廠大腸桿菌磺胺耐藥基因sul1、sul2和sul3的傳播與適應(yīng)性機制研究[D]. 杭州：浙江工商大學(xué). [Fang J H. 2018. Research on the dissemination and fitness mechanism of Escherichia coli sulfonamide-resistant genes sul1，sul2 and sul3 in a pig slaughterhouse[D]. Hangzhou：Zhejiang Gongshang University.] doi：10.27462/d.cnki. ghzhc.2018.000007.

付霞麗，鄭紫方，李洋，肖書奇，李爽. 2021. PRRSV感染繼發(fā)細菌性病原的致病性及耐藥性分析[J]. 中國獸醫(yī)學(xué)報，41（8）：1470-1475. [Fu X L，Zheng Z F，Li Y，Xiao S Q，Li S. 2021. Analysis of pathogenicity and drug resistance of bacterial pathogens after PRRSV infection[J]. Chinese Journal of Veterinary Science，41（8）：1470-1475.] doi：10. 16303/j.cnki.1005-4545.2021.08.04.

甘露，李耘，橋本重陽，張凡. 2020. 大腸埃希菌對碳青霉烯類抗生素的耐藥機制研究[J]. 中國臨床藥理學(xué)雜志，36（18）：2896-2900. [Gan L，Li Y，Hashimoto S，Zhang F. 2020. Resistance mechanism of Escherichia coli to carbapenems[J]. The Chinese Journal of Clinical Pharmaco-logy，36（18）：2896-2900.] doi：10.13699/j.cnki.1001-6821. 2020.18.043.

顧曉曉，鄔琴，陶喬孝慈，馬雪，李劼，周霞，韓猛立，黃新，吳桐忠，張星星，鐘發(fā)剛. 2020. 羊源大腸桿菌對氨基糖苷類藥物耐藥表型及耐藥基因的檢測[J]. 中國獸醫(yī)學(xué)報，40（8）：1566-1570. [Gu X X，Wu Q，Tao Q X C，Ma X，Li J，Zhou X，Han M L，Huang X，Wu T Z，Zhang X X，Zhong F G. 2020. Analysis of aminoglycoside antibiotics resistance phenotype and drug resistance gene of E. coli from sheep[J]. Chinese Journal of Veterinary Science，40（8）：1566-1570.] doi：10.16303/j.cnki.1005-4545.2020.08. 22.

黃嬌，穆青. 2019. 幾種菌的外排泵系統(tǒng)及其抑制劑[J]. 國外醫(yī)藥（抗生素分冊），40（1）：9-12. [Huang J，Mu Q. 2019. Efflux pump systems of several bacteria and their inhibitors[J]. World Notes on Antibiotics，40（1）：9-12.] doi：10.13461/j.cnki.wna.005186.

姜興粲，李冰，王瑋瑋，張繼瑜. 2020. 中國藥典中5種標準菌生長曲線的測定[J]. 浙江農(nóng)業(yè)科學(xué)，61（1）：40-42. [Jiang X C，Li B，Wang W W，Zhang J Y. 2020. Determination of growth curve of five standard bacteria in Chinese pharmacopoeia[J]. Journal of Zhejiang Agricultural Sciences，61（1）：40-42.] doi：10.16178/j.issn.0528-9017. 2020 0114.

李瑞珍，高沂，宋獻藝，李向陽，趙亮. 2011. 后抗生素時代生物飼料添加劑研發(fā)狀況與對策[J]. 山西農(nóng)業(yè)科學(xué)，39（2）：170-173. [Li R Z，Gao Y，Song X Y，Li X Y，Zhao L. 2011. Present situation and strategy for creative research and industry development of biofeed additives at post-antibiotic era in China[J]. Journal of Shanxi Agricultural Sciences，39（2）：170-173.] doi：10.3969/j.issn.1002-2481. 2011.02.23.

李瑋瑋，黃弈睿，吳澤純，薛雅文，冀霞. 2021. 豬源大腸桿菌對氨基糖苷類抗生素耐藥率的Meta分析[J]. 惠州學(xué)院學(xué)報，41（3）：35-44. [Li W W，Huang Y R，Wu Z C，Xue Y W，Ji X. 2021. Meta analysis of resistance rate of Escherichia coli from swine to aminoglycoside antibiotics[J]. Journal of Huizhou University，41（3）：35-44.] doi：10.16778/j.cnki.1671-5934. 2021.03.008.

林文芳，陳勝，萬堃，王春明，林惠榮，于鑫. 2016. 貧營養(yǎng)和痕量抗生素對質(zhì)?？股乜剐赃m應(yīng)度代價的影響[J]. 環(huán)境科學(xué)，37（8）：3224-3228. [Lin W F，Chen S，Wan K，Wang C M，Lin H R，Yu X. 2016. Effects of low-level nutrition and trace antibiotics on the fitness cost of plasmids bearing antibiotic resistance[J]. Environmental Science，37（8）：3224-3228.] doi：10.13227/j.hjkx.2016.08.051.

劉昌孝. 2019. 全球關(guān)注：重視抗生素發(fā)展與耐藥風險的對策[J]. 中國抗生素雜志，44（1）：1-8. [Liu C X. 2019. Global concern：Strategies for antibiotic development and risk of resistance[J]. Chinese Journal of Antibiotics，44（1）：1-8.] doi：10.3969/j.issn.1001-8689.2019.01.001.

劉菊枚，戚南山，廖申權(quán)，李秀，呂敏娜，吳彩艷，李娟，蔡海明，林栩慧，胡俊菁，于林增，肖文婉，張小慧，張健騑，孫銘飛，李國清. 2021. 豬場重要蟲媒及其傳播病原研究進展[J]. 中國畜牧獸醫(yī)，48（5）：1725-1734. [Liu J M，Qi N S，Liao S Q，Li X，Lü M N，Wu C Y，Li J，Cai H M，Lin X H，Hu J J，Yu L Z，Xiao W W，Zhang X H，Zhang J F，Sun M F，Li G Q. 2021. Research progress on the important lnsect vectors and their transmission pathogens in pig farms[J]. China Animal Husbandry & Veterinary Medicine，48（5）：1725-1734.] doi：10.16431/j.cnki.1671-7236.2021.05.024.

劉躍華，韓萌，冉素平，吳久鴻. 2019. 歐洲應(yīng)對抗生素耐藥問題的治理框架及行動方案[J]. 中國醫(yī)院藥學(xué)雜志，39（3）：219-223. [Liu Y H，Han M，Ran S P，Wu J H. 2019. Governance framework and action plan of European Union and member countries on antimicrobial resistance[J]. Chinese Journal of Hospital Pharmacy，39（3）：219-223.] doi：10.13286/j.cnki.chinhosppharmacyj.2019.03.01.

商可心. 2015. 大環(huán)內(nèi)酯類耐藥性對豬鏈球菌適應(yīng)性的影響[D]. 南京：南京農(nóng)業(yè)大學(xué). [Shang K X. 2015. Impact of macrolide resistance on the fitness in Stretococcus suis[D]. Nanjing：Nanjing Agricultural University.]

孫海峰，譚艾娟，呂世明，羅鈺嘉，紀琦琪，艾仁麗，路星星，馮丹. 2021. 高滲對多重耐藥大腸桿菌外排泵和外膜孔道蛋白基因表達及其生長的影響[J]. 中國畜牧獸醫(yī)，48（9）：3473-3482. [Sun H F，Tan A J，Lü S M，Luo Y J，Ji Q Q，Ai R L，Lu X X，F(xiàn)eng D. 2021. Effect of hypertonicity on gene expression and growth of efflux pump and outer membrane protein in multidrug resistant Escherichia coli and its growth[J]. China Animal Husbandry & Vete-rinary Medicine，48（9）：3473-3482.] doi：10.16431/j.cnki. 1671-7236.2021.09.040.

王希越，田媛媛，明明，連麗麗，婁大偉. 2017. 不同碳源條件下遺傳重組大腸桿菌的生長表型研究[J]. 黑龍江畜牧獸醫(yī)，（15）：195-197. [Wang X Y，Tian Y Y，Ming M，Lian L L，Lou D W. 2017. Growth phenotypes of genetically recombinant Escherichia coli under different carbon source conditions[J]. Heilongjiang Animal Science and Veterinary Medicine，（15）：195-197.] doi：10.13881/j.cnki.hljxmsy.2017.1358.

楊焰，游雪甫，李聰然. 2019. 耐藥細菌的適應(yīng)性代價及補償性進化[J]. 中國抗生素雜志，44（7）：784-792. [Yang Y，You X F，Li C R. 2019. The fitness cost and compensatory evolution of drug-resistant bacteria[J]. Chinese Journal of Antibiotics，44（7）：784-792.] doi：10.3969/j.issn.1001-8689.2019.07.004.

袁萬哲，趙興華，何欣，王建昌，張若曦，李麗敏. 2019. 禁抗、限抗、減抗背景下，獸藥創(chuàng)制的出路[J]. 中國動物保健，21（11）：6. [Yuan W Z，Zhao X H，He X，Wang J C，Zhang R X，Li L M. 2019. The way out for the creation of veterinary drugs under the background of prohibition，restriction and reduction of resistance[J]. China Animal Health，21（11）：6.] doi：10.3969/j.issn.1008-4754.2019.11. 002.

周瑤琴，游鳳，高彤，付雅慧，王豪舉，丁紅雷. 2018. 雞源致病性大腸桿菌的分離鑒定及耐藥性分析[J]. 西南大學(xué)學(xué)報（自然科學(xué)版），40（5）：23-29. [Zhou Y Q，You F，Gao T，F(xiàn)u Y H，Wang H J，Ding H L. 2018. Isolation，identification and antimicrobial resistance analysis of chicken-associated pathogenic Escherichia coli[J]. Journal of Southwest University（Natural Science），40（5）：23-29.] doi：10.13718/j.cnki.xdzk.2018.05.005.

Adebusuyi A A，F(xiàn)oght J M. 2011. An alternative physiological role for the EmhABC efflux pump in Pseudomonas fluorescens cLP6a[J]. BMC Microbiology，11（1）：252. doi：10.1186/1471-2180-11-252.

Alva A，Sabido-Ramos A，Escalante A，Bolívar F. 2020. New insights into transport capability of sugars and its impact on growth from novel mutants of Escherichia coli[J]. Applied Microbiology and Biotechnology，104（4）：1463-1479. doi：10.1007/s00253-019-10335-x.

Andersson D I，Levin B R. 1999. The biological cost of antibiotic resistance[J]. Current Opinion in Microbiology，2（5）：489-493. doi：10.1016/S1369-5274（99）00005-3.

Bassetti M，Poulakou G，Ruppe E，Bouza E，van Hal S J，Brink A. 2017. Antimicrobial resistance in the next 30 years，humankind，bugs and drugs：A visionary approach[J]. Intensive Care Medicine，43（1）：1464-1475. doi：10. 1007/s00134-017-4878-x.

de Sousa J M，Balbontin R，Dur?o P，Gordo I. 2017. Multidrug-resistant bacteria compensate for the epistasis between resistances[J]. PLoS Biology，15（4）：e2001741. doi：10. 1371/journal.pbio.2001741.

Deininger K N W，Horikawa A，Kitko R D，Tatsumi R，Rosner J L，Wachi M，Slonczewski J L. 2011. A requirement of TolC and MDR efflux pumps for acid adaptation and GadAB induction in Escherichia coli[J]. PLoS One，6（4）：e18960. doi：10.1371/journal.pone.0018960.

Ferenci T. 2001 . Hungry bacteria—Definition and properties of a nutritional state[J]. Environmental Microbiology，3（10）：605-611. doi：10.1046/j.1462-2920.2001.00238.x.

Gagneux S，Long C D，Small P M，Van T，Schoolnik G K，Bohannan B J M. 2006. The competitive cost of antibiotic resistance in Mycobacterium tuberculosis[J]. Science，312（5782）：1944-1946. doi：10.1126/science.1124410.

Ghosh C，Sarkar P，Issa R，Haldar J. 2019. Alternatives to conventional antibiotics in the era of antimicrobial resistance[J]. Trends in Microbiology，27（4）：323-338. doi：10.1016/j.tim.2018.12.010.

Hartog E，Ben-Shalom L，Shachar D，Matthews K R，Yaron S. 2008 . Regulation of marA，soxS，rob，acrAB and micF in Salmonella enterica serovar Typhimurium[J]. Microbio-logy and Immunology，52（12）：565-574. doi：10.1111/j. 1348-0421.2008.00075.x.

Houston D M. 2011. The use and misuse of antibiotics lea-ding to increasing bacterial resistance and “Antibiotic Choice Frustration”[J]. Advances in Small Animal Medicine & Surgery，24（7）：1-2. doi：10.1016/j.asams.2011. 06.001.

Knopp M，Andersson D I. 2015. Amelioration of the fitness costs of antibiotic resistance due to reduced outer membrane permeability by upregulation of alternative porins[J]. Molecular Biology and Evolution，32（12）：3252-63. doi：10.1093/molbev/msv195.

Maharjan R，F(xiàn)erenci T. 2017. The fitness costs and benefits of antibiotic resistance in drug-free microenvironments encountered in the human body[J]. Environmental Microbiology Reports，9（5）：635-641. doi：10.1111/1758-2229. 12564.

Masi M，Winterhalter M，Pagès J M. 2019. Outer membrane porins[J]. Bacterial Cell Walls and Membranes. Subcell Biochem，92：79-123. doi：10.1007/978-3-030-18768-2_4.

Nové M，Kincses A，Molnár J，Amaral L，Spengler G. 2020. The role of efflux pumps and environmental ph in bacterial multidrug resistance[J]. In Vivo，34（1）：65-71. doi：10.21873/invivo.11746.

Petersen A，Aarestrup F M，Olsen J E. 2009. The in vitro fitness cost of antimicrobial resistance in Escherichia coli varies with the growth conditions[J]. FEMS Microbio-logy Letters，299（1）：53-59. doi：10.1111/j.1574-6968.2009. 01734.x.

Vogwill T，MacLean R C. 2015. The genetic basis of the fitness costs of antimicrobial resistance：A meta-analysis approach[J]. Evolutionary Applications，8（3）：284-295. doi：10.1111/ eva.12202.

Yasufuku T，Shigemura K，Shirakawa T，Matsumoto M，Nakano Y，Tanaka K，Arakawa S，Kinoshita S，Kawabata M，F(xiàn)ujisawa M. 2010. Correlation of overexpression of efflux pump genes with antibiotic resistance in Escherichia coli strains clinically isolated from urinary tract infection patients[J]. Journal of Clinical Microbiology，49（1）：189-94. doi：10.1128/jcm.00827-10.

收稿日期：2021-12-28

基金項目：國家自然科學(xué)基金項目（31760749）;貴州省科技計劃項目（黔科合平臺人才〔2017〕5788號，黔科合基礎(chǔ)-ZK〔2022〕一般129號）;貴州大學(xué)自然科學(xué)專項（特崗）科研基金項目（貴大特崗合字〔2020〕25號）

通信作者：宋旭琴（1991-），https：//orcid.org/0000-0002-8671-3161，博士，主要從事獸醫(yī)藥理學(xué)與毒理學(xué)研究工作，E-mail：xqsong@gzu.edu.cn;呂世明（1963-），https：//orcid.org/0000-0001-5939-0574，博士，教授，主要從事獸醫(yī)藥理學(xué)與毒理學(xué)研究工作，E-mail：lvlvsm@163.com

第一作者：廖書丹（1994-），https：//orcid.org/0000-0001-8419-1702，研究方向為獸醫(yī)藥理毒理學(xué)，E-mail：871823310@qq.com