銅綠假單胞菌對鱟素耐藥前后的差異表達基因及SNP變化研究

洪軍 衛夏怡 吉冰潔 葉延欣 程天賜

（河南城建學院生命科學與工程學院，平頂山 467036）

銅綠假單胞菌（Pseudomonas aeruginosa，PA）作為一種人畜共患條件性致病菌，通常感染肺部、尿道、燒傷等部位，已成為醫院和養殖廠等感染的三大致病菌之一。由于對抗生素的多重耐藥性，在臨床抗感染治療中有著特殊重要的地位。該菌極易形成生物膜，單獨使用抗生素的治療效果并不理想。目前對銅綠假單胞菌感染的防治是醫學界迫切需要解決的問題，對其耐藥性產生的研究為預防和治療由銅綠假單胞菌引起的感染具有重要的指導意義。

鱟素（Tachyplesin I）是1988年由Nakamura等［1］從中國鱟血細胞的酸性提取物中分離出來的一種具有17個氨基酸殘基的典型環狀β-折疊結構的陽離子抗菌肽。已研究證明鱟素具有廣譜抗菌活性，可開發成抗菌藥物、抗腫瘤藥物、免疫調節劑和飼料添加劑等。2019年12月份國家農業農村部發布194號公告指出：自2020年元旦起，我國飼料中全面禁止添加抗生素，減少濫用抗生素造成的危害，維護動物源食品安全和公共衛生安全。因此尋找一些新型抗菌劑替代抗生素已成為新時代綠色環保養殖行業的重中之重。作為抗生素替代品的侯選開發抗菌肽之一，除了作為抗菌藥物使用外，也可以作為新型功能性添加劑應用到飼料中。然而，已有研究報道表明一些細菌（尤其人類病原菌）已對部分抗菌肽產生了耐藥性［2-3］。本課題組在實驗室條件下通過長期連續增高濃度鱟素誘導后，銅綠假單胞菌CGMCC 1.2620和ATCC 27853菌株均能對鱟素產生不同程度的耐藥性［4-5］，但是銅綠假單胞菌耐藥前后能引起哪些基因的差異表達以及是否引起SNP和sRNA的變化，以及sRNA是否能調控這些差異表達基因尚不清楚。

隨著高通量測序技術的迅速發展，細菌轉錄組學的研究可幫助人們探究細菌耐藥前后發生差異表達的基因以及篩選出具有調控作用的非編碼RNA。sRNA（small RNA）是一類廣泛存在于原核生物基因組中能被轉錄但不編碼蛋白質的一類RNA 分子，長度為50-300 nt的核苷酸，位于基因組2%-12%的基因間區域。有研究表明，細菌為適應各種不同的外界環境，通過較為復雜的基因調控網絡來對其做出反應，而其中sRNA發揮了重要的作用［6-7］，并揭示sRNA與細菌耐藥性密切相關［8］。通過轉錄組學研究手段篩選細菌中存在的sRNA并研究sRNA所發揮的功能，對揭示細菌耐藥基因轉錄后調控具有重要意義。Miller等［9］表明銅綠假單胞菌中存在著一系列sRNA，其中Rsm W的過度表達及Rsm Z的表達水平降低均有助于增強生物膜的形成，而Rsm Z的過度表達則會抑制生物膜的形成。劉翠翠［10］利用RNA-Seq技術對福氏志賀菌的耐藥株中差異表達的sRNA進行了GO富集分析和KEGG Pathway分析，為揭示福氏志賀菌sRNA與其耐藥性的關系奠定基礎。Zhang等［11］表明在環境適應過程中，細菌利用sRNAs來抑制或激活其蛋白質組的大部分表達，并且已經鑒定出了兩個sRNAs（Sr0161和ErsA）主要是與負責碳青霉烯類抗生素攝取的主要孔蛋白OprD mRNA相互作用。這兩個sRNAs與oprD的5'UTR堿基配對導致細菌對美羅培南的耐藥性增加。Sr0161在沒有Hfq的情況下，正調控PagL的表達，降低其促炎特性并導致多粘菌素抗性。通過RNASeq技術等轉錄組學研究手段可篩選出對細菌外排泵系統、生物被膜形成等具有調控作用的sRNA及靶基因，從RNA水平揭示了sRNA對細菌耐藥性形成的影響。

SNP（single nucleotide polymorphism）主要是指基因組在核苷酸水平上發生變異引起的DNA序列多態性，包括單個堿基轉換、顛換以及單個堿基插入和缺失等。在基因組DNA中，任何堿基都可能發生突變，因此SNP可存在于基因的任意序列中，可能位于非編碼序列中，也可能位于基因序列中。張棟等［12］采用直接測序法從P0代尼羅羅非魚的β2m基因組序列中篩選出30個SNP，獲取其位置信息，并利用Popgen32和picc -calc軟件分析其遺傳參數，結果表明24個SNPs的基因型和等位基因頻率與無乳鏈球菌的耐感/易感態顯著相關。

本研究通過對銅綠假單胞菌野生株CGMCC1.2620（PA1.2620）和 高 抗 鱟 素 突 變 株CGMCC1.2620-99（PA-99）進行轉錄組測序，在轉錄組水平上通過對差異表達基因以及差異表達的sRNA的靶基因進行功能富集和SNP分析，找到與抗藥相關的差異表達的sRNA基因及其相關通路信息以及SNP變化，為進一步探討銅綠假單胞菌對鱟素的耐藥性機制提供理論依據。

1 材料與方法

1.1 材料

實驗菌株：銅綠假單胞菌CGMCC1.2620（PA1.2620）由中國微生物菌種保藏中心（CGMCC）提供。耐鱟素突變株銅綠假單胞菌PA-99是由PA1.2620菌株在本實驗室長期連續暴露的增高鱟素濃度選擇壓力下轉接99代獲得的突變株。鱟素對PA1.2620和PA-99菌株的最小抑制濃度值分別為10 μg/mL和150 μg/mL。

抗菌肽：鱟素（tachyplesin I，含有兩個二硫鍵的環行肽）由吉爾生化（上海）有限公司合成，純度為95%以上。

1.2 方法

1.2.1 測序樣品制備及轉錄組測序 PA1.2620原始菌株和PA-99突變株在30℃恒溫搖床上180 r/min培養12 h后，按OD600=0.2的濃度轉接到新鮮的MHB培養基中，并生長至OD600=1.0時開始收集。用細菌總RNA純化試劑盒提取和純化總RNA（Sangon Biotech Co.，Ltd.，Shanghai，China）。為 了 確 保總RNA能滿足轉錄組測序分析的質量，分別采用Nanodrop、Qubit 2.0、Agilent 2100、電泳方法檢測RNA樣品的純度、濃度、完整性和是否有基因組DNA污染等。每個樣品3個重復。總RNA樣本檢測合格后，委托百邁客生物科技有限公司對樣品進行cDNA文庫構建和有參轉錄組測序分析；測序平臺Illumina HiSeq2500。

1.2.2 基因表達量及差異表達基因功能分析 利用DESeq軟件進行差異表達分析，基因表達量的計算使用RPKM（Reads Per Kilobase of transcript per Million fragments mapped）法。將Fold Change≥2且FDR＜0.01作為篩選差異表達基因標準，即log2FC≥1 是上調表達，log2FC≤1下調表達［13］。得到差異表達基因后，對差異表達基因進行功能注釋、KEGG Pathway和GO功能分析［14-16］。利用KEGG和Gene Ontology數據庫找出與整個基因組背景相比在顯著差異表達基因中q-Value≤0.05的顯著性富集的Pathway及GO條目。

1.2.3 RT-qPCR驗證 用RT-qPCR來驗證轉錄組測序檢測結果的準確性。分別從PA-99突變株和原始菌株PA1.2620中選擇5個靶基因進行RT-qPCR分析。每個基因設3個重復。此外，使用KAPA SYBR FAST qPCR試劑盒在ABI7900實時熒光定量PCR儀上進行RT-qPCR測定。利用Primer Premier 5軟件根據銅綠假單胞菌基因組序列設計特異性引物。反應條件：95℃ 2 min，94℃ 20 s，60℃ 20 s，72℃ 30 s，40個循環。以16S rRNA為內參，基因的相對表達水平用2-ΔΔCt方法定量［17］。

1.2.4 差異表達sRNA及靶基因的分析 利用miRDeep2 軟件將測序Read比對到基因組上后，看那段區域是否能形成發夾結構，且Read位于臂上，得到一個分值，篩選后得到預測結果，用的軟件是miRDeep2。篩選出PA1.2620與PA-99差異表達的sRNA；對篩選出來的差異sRNA靶基因數據進行處理，找出差異表達的基因，并對其進行GO功能和KEGG通路富集分析。

1.2.5 差異表達基因和sRNA靶基因的SNP分析 基于各樣品Reads與參考基因組序列的Bowtie比對結果，使用GATK軟件識別測序樣品與參考基因組間的單堿基錯配，查找基因區潛在的單核苷酸多態性SNP位點。

在有參轉錄組測序的結果上進行比對，對實驗中所用對照組與處理組的各個樣品組內發生的SNP進行比對統計，并對其類型和數目統計結果進行分析。通過百邁客云平臺的基于基因ID過濾GFF信息小工具，在銅綠假單胞菌的gff文件中提取得到基因對應的位置信息。將SNP的位置信息與基因對應的位置信息比對，可以得到差異基因中產生SNP的數目和類型。所得數據主要通過百邁客云平臺分析。

2 結果

2.1 測序數據及表達量分析

對原始測序數據進行過濾和質量控制，之后對兩個樣品的轉錄組測序數據進行評估，包括與參考基因組的比對（比對效率都在87.65%及以上）、rRNA統計、測序隨機性評價以及Reads在參考基因組上的分布等，這些分析結果表明該測序結果較為理想。在銅綠假單胞菌6 318個基因中，PA-99與PA1.2620菌株的表達差異的基因有1 459個，其中上調表達的基因為672個，下調表達的基因為787個。其中log2FC值上調表達倍數為3.0以上有35個基因，下調表達倍數為3.0以上基因有58個，且下調倍數log2FC大于11為28個基因。此結果表明，銅綠假單胞菌對鱟素產生耐藥性導致許多基因發生差異表達。

2.2 PA1.2620與PA-99菌株差異表達基因的功能富集分析

2.2.1 差異表達基因的GO功能富集分析 由GO富集的前20個GO條目圖1可知，在分子功能方面發生富集的亞類主要有：核苷酸結合（nucleotide binding，18）、甲酸脫氫酶（NAD+）活性（formate dehydrogenase（NAD+）activity，6）、錯 配DNA結合（mismatched DNA binding，4）、受體結合（receptor binding，3）等等；在生物過程發生主要富集的GO條目有：細胞呼吸（cellular respiration，4）、錯配修 復（mismatch repair，4）、細 胞 溶 解（cytolysis，4）等；在細胞組分主要發生富集的GO 條目是細胞膜的組成部分（integral component of membrane，158）。最為富集的前3條GO Term依次為：integral component of membrane、nucleotide binding、formate dehydrogenase（NAD+）activity。在GO功能研究中，我們重點關注GO富集的前20個GO條目。

圖 1 PA1.2620 與 PA-99 菌株差異基因的 GO 富集圖Fig. 1 GO enrichment analysis of differentially expressed genes between PA1.2620 vs. PA-99 strains.

2.2.2 KEGG通路富集分析 根據PA1.2620與PA-99菌株差異表達基因的KEGG結果可知，此次分析將395個基因注釋到KEGG的98條通路上，涉及到各種氨基酸、脂質、能量、糖類等代謝（266個基因，最多）、ABC轉運（36個基因）、雙組分系統（43個基因）、細菌分泌系統（16個基因）、細菌趨化性（16個基因）、氨酰生物合成（14個基因）等通路。由圖2富集的前20個通路可知，無顯著富集的通路。其中細菌趨化性（bacterial chemotaxis）、氨酰生物合成（aminoacyl-tRNA biosynthesis）通路中所涉及的基因通常與抗藥性密切相關。但Toluene degradation、Chlorocycchexane and chiorobenzene degradgfion、Phosphonate and phosphinate melaboisn、Novobiocin biosynthesis、Aminobenzoate degradation這些通路需要我們重點關注。

圖2 PA1.2620與PA-99菌株差異基因的KEGG富集圖Fig.2 KEGG pathway enrichment analysis of differentially expressed genes between PA1.2620 vs. PA-99 strains

2.3 RT-qPCR驗證

為了驗證轉錄組測序及后續SNP和sRNA靶基因等分析的結果，根據文獻報道挑選了疑似耐藥性的ParR、phoB、PA3573和opr86差異表達基因以及上調表達量最高的PA0912基因進行RT-qPCR測定。如圖3所示，所有選擇的差異表達基因在RNA-Seq和RT-qPCR結果之間顯示一致的表達模式，但是整體而言RT-qPCR數值較RNA-Seq表達量數值偏高。

圖3 RT-qPCR驗證結果Fig. 3 RT-qPCR validation results

2.4 SNP分析

2.4.1 SNP突變類型統計 對PA-99突變株各樣品對比PA1.2620樣品篩選出SNP位點數目、轉換類型比例、顛換類型比例以及雜合型SNP位點比例進行統計，如表1所示。以及對比分析得到各個處理樣品組的SNP類型進行數據統計如圖4所示。根據圖和表可知：相對于對照基因組來講，各個樣品組間發生的SNP類型頻率相對一致，在其中發生轉換類型的頻率遠高于顛換類型，在轉換之中數量最多的為A/G，其次是C/T。

表1 SNP位點統計表Table 1 Statistical table of SNP sites

圖4 PA-99突變株樣品的SNP突變類型分布圖Fig. 4 Distribution of SNP mutation types in PA-99 mutant samples

2.4.2 差異表達基因上產生的SNP分析 以PA1.2620菌株作為對照組，將PA-99組菌株與其對比結果匯總在一起查找差異表達基因的SNP。由PA1.2620與PA-99菌株的差異表達基因上產生的SNP可知，在1 459個差異表達基因中，共有SNP發生的基因數目為43個（25個上調和18個下調基因），其中新基因是10個如Novel_244、Novel_137、Novel_28等，差異表達值最大的基因是未知功能基因（ID為Novel_28基因）。對這些發生SNP基因過濾掉僅在一個生物學重復中發生的堿基突變，在這43個差異表達基因中SNP為PA-99菌株3個重復組共同發生的基因數目為22個（12個表達上調和10個表達下調基因）如表2所示，此結果也表明不僅不同處理組間存在SNP變化，而且重復組內也存在SNP變化的差別。這22個發生SNP的基因中包括5個新基因，7個假定蛋白基因，10個已知功能基因，SNP主要發生在基因表達量log2FC值在1-3倍之間。Novel_745基因的差異表達值最大（上調表達log2FC=2.58倍以上），是一種未知功能蛋白，但僅發生了一個堿基突變。這些發生SNP突變的已知基因，如pagL 基因是脂質A脫酰基酶，是參與脂質A代謝過程功能，下調表達log2FC=2.3倍以下；opr86是一個編碼外膜蛋白基因；PmrB雙組分調節系統信號傳導激酶，與脂多糖的修飾有關。

表2 PA1.2620與PA-99菌株差異表達基因的SNP信息Table 2 SNP analysis of differentially expressed genes in PA1.2620 vs. PA-99 strains

2.5 PA1.2620與PA-99菌株差異表達sRNA分析

2.5.1 差異表達sRNA及其靶基因數目分析 由表3和圖5可知：按照篩選差異表達基因的標準，篩選到PA1.2620與PA-99菌株的差異表達sRNA共有11個，其中1個sRNA上調表達，10個下調表達，且均是未知功能sRNA。通過對其靶基因進行篩選，篩選到的靶基因總數為5 840個，其中差異表達靶基因為1 247個。每個sRNA所對應的靶基因個數差別較大，其中sRNA Novel_583有968個靶基因包括223差異表達靶基因，而Novel_33有25個靶基因，僅存在4個差異表達靶基因。這些結果表明，sRNA所對應的靶基因越多，調控的基因功能越多，且存在不同sRNA基因調控同一靶基因的現象。因其存在不同sRNA有相同的靶基因，過濾掉重復后其靶基因總數為4 027個，其中差異表達靶基因為863個。隨后重點關注這些sRNA的差異表達靶基因。

圖5 PA1.2620與PA-99菌株差異表達的sRNA靶基因數目統計圖Fig.5 Statistical diagram of the number of sRNA target genes of differentially expressed genes in PA1.2620 vs. PA-99 strains

表3 PA1.2620與PA-99菌株差異表達sRNATable 3 sRNA of differentially expressed genes in PA1.2620 vs. PA-99 strains

2.5.2 sRNA差異表達靶基因GO功能富集分析 根據PA1.2620與PA-99的 差 異 表 達sRNA靶 基 因的GO富集前20個GO Term可知（圖6），在分子功能大類下主要發生富集的GO Term有：高絲氨酸激酶活性（homoserine kinase activity，3）、5-羧甲基-2-羥基黏液酸δ-異構酶活性（5-carboxymethyl-2-hydroxymuconate delta-isomerase activity，3）、轉錄因子結合（transcription factor binding，8）等等；在生物學過程大類下主要發生富集的GO Term有：組氨酸生物合成（histidine biosynthetic process，7）、蘇氨酸氨基甲酰胺腺苷生物合成過程（threonylcarbamoyladenosine biosynthetic process，3）、細胞蛋白質代謝過程（cellular protein metabolic process，13）等。最富集的前5個GO Term依次為：組氨酸生物合成（histidine biosynthetic process）、高絲氨酸激酶活性（homoserine kinase activity）、5-carboxymethyl-2-hydroxymuconate delta-isomerase activity（3）、threonylcarbamoyladenosine biosynthetic process、transcription factor binding，均為極顯著富集水平。

圖6 PA1.2620與PA-99菌株表達差異sRNA靶基因的GO富集圖Fig.6 GO enrichment of sRNA target genes of DEGs in PA1.2620 vs. PA-99 strains

2.5.3 差異sRNA靶基因的KEGG富集分析 通過KEGG富集圖7可知，沒有顯著富集的通路。但是Biosynthesis of unsaturated fatty acids、Phosphotransferase system（PTS）、Degradation of aromatic compounds、Carbon metabolism等通路所涉及的差異表達基因富集較可靠。

圖7 PA1.2620與PA-99菌株表達差異sRNA靶基因的KEGG富集圖Fig.7 KEGG enrichment of differentially expressed sRNA target genes in PA1.2620 vs. PA-99 strains

2.5.4 差異表達sRNA靶基因的SNP分析 對PA1.2620與PA-99菌株的差異表達的11個sRNA發生的SNP的數據為0個。對其差異表達sRNA的863個差異表達靶基因上發生的SNP的數據進行分析，共有22個靶基因發生SNP變化，其中在PA-99菌株的3個生物學重復中兩者或兩者以上共同發生的有10個靶基因（表4）。sRNA靶基因發生SNP變化的較少，且多發生在較低差異表達靶基因，僅有pvdL基因有2個位點突變，其余9個基因均發生一個位點突變。這些SNP的變化主要涉及的sRNA所作用的靶基因除了3個靶點新基因外，還涉及到編碼外膜蛋白基因opr86；雙組分調節系統信號傳導激酶pmrB，該雙組分調控系統調控PA3552-PA3559 脂多糖修飾操縱子，該操縱子是陽離子抗菌肽和多粘菌素B陽離子抗菌肽耐藥所必需的；以及編碼肽合成酶基因pvdL，主要參與銅綠假單胞菌鐵載體合成和脂類生物學合成過程，參與催化活性的分子功能。sRNA如何參與這些靶基因的功能有待進一步研究及驗證。

表4 sRNA相關靶基因的SNP變化分析Table 4 SNP analysis of sRNA target genes

3 討論

通過對耐藥菌株中差異表達基因的分析，有助于研究者們從眾多差異表達的基因中篩選出與細菌耐藥性密切相關的基因及其可能參與的代謝途徑，可為揭示細菌耐藥性機制提供新的思路，同時為新型抗耐藥菌藥物的研究奠定理論基礎。sRNA是一類新發現的基因表達調控因子，通過與靶mRNA或靶蛋白配對，從而調控細胞的生理功能以應對各種環境變化，尤其在細菌耐藥性的產生方面至關重要。這些差異表達的sRNA主要是通過作用于相應的靶mRNA，進而發揮其調控作用。在基因組水平上，由單個核苷酸變異引起的DNA序列多態性即稱為SNP，SNP可進行復雜性狀的關聯分析，從而揭示相關性狀的遺傳機理。基因突變產生的耐藥性是一個長期積累的過程，不同于轉錄組，基因組較穩定，不易發生突變，在經過特殊的誘導處理后可產生定向的突變以應對環境中的選擇性壓力。這些突變可能是菌株生存環境變化后為適應新環境做出的一些應對調整，可能與耐藥性、毒力等相關基因的特異性突變無關。對于差異表達基因、sRNA所對應靶基因與相應SNP進行分析，可從轉錄組和基因組層面來預測銅綠假單胞菌對鱟素耐藥性產生的原因，從而對于耐藥性機制的預測作為相互補充。因此，本研究基于轉錄組測序的基礎上，對差異表達基因的SNP和sRNA的靶基因進行了預測分析，以期為銅綠假單胞菌對鱟素的耐藥性的產生提供一定的理論支持。

本研究通過對高抗鱟素銅綠假單胞菌PA-99菌株與PA1.2620敏感菌株進行轉錄組測序，在轉錄組水平上通過對差異表達基因及sRNA靶基因進行富集分析，發現這些基因大多與氨基酸、糖類等物質的代謝通路、DNA轉錄、金屬離子的轉運以及代謝、細菌趨化性、雙組分調控系統等通路有關以及涉及到細胞膜組成部分、核苷酸結合及酶的活性等功能。對于差異表達基因和sRNA靶基因的SNP變化，主要涉及到編碼脂質A脫酰基酶、外膜蛋白、冷休克蛋白、鐵代謝、以及與脂多糖修飾相關等已知基因和一些編碼的假定蛋白、部分編碼的新基因有關，這些靶基因的SNP變化是否與抗藥性相關有待進一步驗證。這些結果表明sRNA在調控銅綠假單胞菌對鱟素的抗藥性方面也發揮了重要作用。

典型的雙組分調控系統由組氨酸激酶和相關的應答調控蛋白組成，組氨酸激酶感受外界環境變化后調節應答蛋白的磷酸化水平，而后對一系列應答基因進行調控。雙組分調控系統在銅綠假單胞菌的致病性、毒力、生物膜形成和耐藥性等方面起著重要作用［18］。PmrA-PmrB雙組分調控系統在革蘭氏陰性菌對多黏菌素B耐藥過程中起重要作用［19］。本研究發現耐鱟素銅綠假單胞菌PA-99菌株中雙組分調控系統有43個差異表達基因，包括一些受調控的抗性基因，如pmrA-pmrB、反應調節因子PhoP（基因ID：Novel_194）。PhoA編碼的堿性磷酸酶（基因ID：Novel_498，499）和PA-99突變體中的兩個成分反應調節因子PhoB；并預測到sRNA靶基因中有29個差異表達基因，且僅編碼修飾脂多糖PmrB基因的SNP發生堿基突變。此結果表明，雙組分調控系統在銅綠假單胞菌對鱟素的耐藥中同樣發揮作用，且sRNA也有可能參與調節雙組分系統表達。

氨酰-tRNA生物合成（aminoacyl-tRNA biosynthesis）的功能是將氨基酸與含有相應反密碼子的tRNAs精確匹配。除了在蛋白質合成和某些細菌革蘭氏陽性菌肽聚糖交聯或膜磷脂修飾中發揮重要作用外，氨酰基-tRNA生物合成還直接參與抗生素生物合成和耐藥性［20］。例如，MprF蛋白通過催化細胞內膜脂的氨基酰化來調控細胞壁對陽離子抗菌肽的通透性，從而介導抗菌肽的耐藥性［21］。膜脂的氨酰化通過調控細胞膜雙分子層的凈負電荷，其方式依賴于氨酰化 tRNA底物。在本研究中，氨酰基-tRNA生物合成通路在高抗鱟素銅綠假單胞菌PA-99菌株中有14個差異表達基因以及預測的sRNA靶基因中有10個差異表達基因在此通路中，但這些基因無SNP發生變化。這個結果表明銅綠假單胞菌對鱟素抗藥性也涉及到此通路。

生物膜形成的耐藥機制在銅綠假單胞菌中起著關鍵作用。銅綠單胞菌生物膜的形成途徑包括群體感應、細菌分泌系統和化學感覺系統等。在轉錄組數據中，我們發現有許多與生物形成相關的差異表達基因及部分sRNA靶基因，如鞭毛組裝（ko02040）、細菌分泌系統（ko03070）和細菌趨化（ko02030）通路，這些通路涉及基因之前被認為與生物膜形成途徑有關。在前期研究中發現，PA-99突變體比原始菌株PA1.2620更容易形成生物膜，胞外多糖含量也高于原始菌株［5］。已有學者探討了sRNA在某些細菌的鐵代謝、糖代謝以及氨基酸代謝中所發揮的調控作用［22］。針對鐵代謝，許多細菌都會進化出自身特殊的鐵調控系統來限制鐵的攝取。Fur作為一種較為重要的鐵調控蛋白參與到sRNA對鐵的調節過程中。在本研究中，sRNA的靶基因vreR是與鐵代謝相關基因，也發生了SNP堿基突變。葡萄糖、氨基酸等是細菌的能源物質，在細菌耐藥過程中也是受到了嚴格的調控。有研究指出細菌可以通過sRNA調節群體感應系統，進而調控毒力產生、抗菌藥物耐藥等過程。但是，相對于大腸桿菌、沙門氏菌等細菌，關于銅綠假單胞菌的sRNA調控耐藥機制的相關報道相對來說較少，這其中的機制尚不清楚。

4 結論

銅綠假單胞菌對鱟素的耐藥性導致許多基因發生差異表達，且有些差異表達基因受到sRNA的表達調控，并有極個別相對低差異表達基因的SNP發生堿基突變，這些突變基因主要與膜蛋白修飾、鐵代謝和編碼的假定蛋白等有關。盡管本研究未涉及這些SNP突變基因的功能驗證，但在文獻中已有報道這些已知基因是與抗藥性相關，且已選擇5個代表性差異表達基因通過熒光定量證實差異表達基因的表達趨勢與轉錄組測序表達量基本一致。