賈新蕾，黃增朝，楊林狄，呂 靜，李妍萍，簡紀常，黃郁蔥

（1. 廣東海洋大學(xué) 水產(chǎn)學(xué)院/廣東省水產(chǎn)動物病害防控與健康養(yǎng)殖重點實驗室，廣東 湛江 524088;2. 南方海洋科學(xué)與工程廣東省實驗室（湛江），廣東 湛江 524006）

海豚鏈球菌（Streptococcus iniae）是一種重要的水生動物致病菌，1976 年P(guān)ier 等[1]從亞馬遜江豚（Inia geoffrensis）中首次分離獲得，之后陸續(xù)發(fā)現(xiàn)它可感染除南極洲以外的所有魚類[2]。我國也有海豚鏈球菌感染卵形鯧鲹（Trachinotus ovatus）[3 - 4]、尼羅羅非魚（Oreochromis niloticus）[5-6]、鱘魚（Acipenser sinensis）[7 - 8]、銀鼓魚（Selenotocamultifasciata）[9]、小黃魚（Larimichthys polyactis）[10]、黃鰭鯛（Acanthopagrus latus）[11]等發(fā)病的相關(guān)報道。海豚鏈球菌具有感染宿主范圍廣、傳染性強、導(dǎo)致死亡率高等特點，給養(yǎng)殖業(yè)造成巨大的經(jīng)濟損失，并受到了國內(nèi)外水產(chǎn)行業(yè)工作者的高度關(guān)注。

海豚鏈球菌病的暴發(fā)除了與病原菌本身的毒力有關(guān)，一些環(huán)境因子在疾病發(fā)生發(fā)展過程中亦發(fā)揮重要作用，如溫度、pH 值等[12]，尤其是養(yǎng)殖的環(huán)境溫度與該病的發(fā)生呈現(xiàn)比較明顯的相關(guān)性，發(fā)病時水溫通常在27 ℃以上[13]。祝璟琳等[14]開展的流行病學(xué)研究表明，在32 ℃以上的溫度下羅非魚更容易爆發(fā)鏈球菌病，羅非魚免疫力在高溫條件下受到明顯的抑制[15]，鏈球菌對宿主的粘附、定植、入侵和抵抗免疫清除能力明顯增強[16 - 17]。同時，較高溫度時海豚鏈球菌產(chǎn)生的胞外產(chǎn)物中可能含有較多致病因子，從而影響細菌毒力[18]，尼羅羅非魚在感染海豚鏈球菌后其肝臟、腎臟在高溫時病理損傷更嚴重[19]。以上研究均表明在較高的溫度下病原鏈球菌的毒力顯著增強，從而引起鏈球菌病的爆發(fā)。

轉(zhuǎn)錄組測序（RNA-Seq）技術(shù)是一項能夠檢測有機體整體轉(zhuǎn)錄水平的新興技術(shù)手段，目前已有學(xué)者利用轉(zhuǎn)錄組技術(shù)研究細菌的毒力調(diào)控。如：Li 等[20]利用RNA-seq 技術(shù)分析了提高溫度和營養(yǎng)成分對副溶血性弧菌毒力的影響，結(jié)果發(fā)現(xiàn)溫度和營養(yǎng)成分的提高可顯著提高副溶血性弧菌的毒力；王松艷[21]對不同溫度下培養(yǎng)的無乳鏈球菌進行轉(zhuǎn)錄組測序，篩選出7 個毒力相關(guān)基因；師若萍[22]運用高通量RNA 測序和非靶向代謝組學(xué)的研究方法，探究大腸桿菌在不同溫度下的生長及其在轉(zhuǎn)錄和代謝水平的變化，了解其對于不同溫度的生理調(diào)節(jié)機制，以上研究結(jié)果表明溫度可廣泛影響細菌毒力和基因的表達。但目前尚未見溫度對海豚鏈球菌致病性的影響以及其基因轉(zhuǎn)錄表達的研究報道。卵形鯧鲹是我國海水養(yǎng)殖魚類中產(chǎn)量位居第二的優(yōu)良品種，高溫期易暴發(fā)海豚鏈球菌病[3-4]。本研究以卵形鯧鲹源海豚鏈球菌Tozj-1為研究對象，選用RNA-Seq 技術(shù)對不同溫度（25 ℃和35 ℃）培養(yǎng)的海豚鏈球菌進行高通量測序分析，篩選差異表達基因，對差異表達基因進行GO 功能和KEGG 通路富集分析，探究溫度對海豚鏈球菌基因轉(zhuǎn)錄表達的影響及富集的功能和信號通路，以期為揭示溫度與海豚鏈球菌致病性積累基礎(chǔ)數(shù)據(jù)，為后續(xù)進一步研究海豚鏈球菌的致病機理提供理論依據(jù)。

1 材料和方法

1.1 材料本研究測序和試驗用菌株Tozj-1來源于廣東省水產(chǎn)動物病害防控與健康養(yǎng)殖重點實驗室。Tozj-1為實驗室于2020 年8 月從廣東省湛江市某深海網(wǎng)箱養(yǎng)殖的患病卵形鯧鲹（Trachinotus ovatus）分離獲得。分離的菌株經(jīng)16S rRNA 測序、生理生化鑒定為海豚鏈球菌（S.iniae），通過人工感染卵形鯧鲹確認其為病原菌，并保存于廣東省水產(chǎn)動物病害防控與健康養(yǎng)殖重點實驗室。

1.2 方法

1.2.1 生長曲線測定將Tozj-1菌株劃線接種于BHI 固體培養(yǎng)皿上，置于28 ℃恒溫培養(yǎng)箱中培養(yǎng)24 h。用無菌接種環(huán)挑取單個菌落接種于新鮮的無菌BHI 液體培養(yǎng)基，于28 ℃，120 r·min-1下振蕩過夜擴大培養(yǎng)。將上述菌液以1∶50（V/V）的比例接種到無菌的BHI 培養(yǎng)基，分別置于25 ℃和35 ℃恒溫搖床，150 r·min-1振蕩培養(yǎng)，每個溫度下設(shè)置3 個重復(fù)。每隔1 h 取樣，測定菌液的吸光值（OD600），并繪制不同溫度下的海豚鏈球菌生長曲線。

1.2.2 人工回歸感染試驗為了探究不同溫度對卵形鯧鲹源海豚鏈球菌毒力的影響，本研究使用卵形鯧鲹進行體內(nèi)攻毒試驗。實驗前，將28 ℃擴大培養(yǎng)的海豚鏈球菌以1∶50（V/V）的比例接種到新鮮BHI 培養(yǎng)基中，分別置于25、35 ℃恒溫搖床中培養(yǎng)至對數(shù)生長期（OD600=0.9～1.1），收集菌體，用無菌PBS 緩沖液將菌液稀釋至1.0×107CFU·mL-1備用。

將360 尾健康的卵形鯧鲹[體質(zhì)量（50±6）g·尾]隨機分為4 組（2 個實驗組和2 個對照組），每組設(shè)3 個重復(fù)，每個重復(fù)30 尾魚。25 ℃實驗組于水溫25 ℃飼養(yǎng)，每尾腹腔注射25 ℃體外培養(yǎng)的Tozj-1菌懸液，對照組注射無菌PBS；35 ℃實驗組于水溫35 ℃飼養(yǎng)，每尾腹腔注射35 ℃體外培養(yǎng)的Tozj-1菌懸液，對照組注射等量無菌PBS。觀察攻毒后14 d 內(nèi)試驗魚的發(fā)病情況，記錄并統(tǒng)計死亡魚數(shù)量，通過解剖瀕臨死亡魚并進行病原分離鑒定，以確定感染病原，鑒定方法同1.1。

1.2.3 總RNA 提取、cDNA 文庫構(gòu)建和高通量測序?qū)?8 ℃下培養(yǎng)的菌株Tozj-1以1∶50（V/V）的比例稀釋于新鮮BHI 肉湯中，分別于35 ℃和25 ℃條件下培養(yǎng)至對數(shù)生長期（OD600=0.9～1.1）。每個溫度設(shè)置3 個生物學(xué)重復(fù)。其中，以25 ℃培養(yǎng)的海豚鏈球菌作為對照組，標記為CK-1、CK-2、CK-3；以35 ℃培養(yǎng)的海豚鏈球菌作為實驗組，標記為Test-1、Test-2、Test-3。每個樣品取1 mL 菌液，12 000 r·min-1離心2 min，棄去上清，加入100 μL 溶菌酶，混勻，置于37 ℃裂解15 min，按照Trizol 試劑盒操作說明進行菌體的總RNA 抽提，使用Agilent 2100 生物分析儀評估RNA 質(zhì)量并通過RNase free 瓊脂糖凝膠電泳進行檢測，經(jīng)質(zhì)量檢測合格后純化總RNA，cDNA 文庫構(gòu)建和轉(zhuǎn)錄組測序工作委托廣州基迪奧生物科技有限公司完成。

1.2.4 原始數(shù)據(jù)質(zhì)控與差異表達基因分析原始數(shù)據(jù)的存儲采用FASTQ 格式，通過fastp[23]（https://github.com/OpenGene/fastp）對下機數(shù)據(jù)進行質(zhì)量檢測和數(shù)據(jù)過濾，去除原始數(shù)據(jù)中質(zhì)量低、接頭污染、全部都是A 堿基以及未知堿基N 含量過高的Reads 后最終獲得Clean reads。利用工具Bowtie 2[24]建立參考基因組索引，將質(zhì)控分析后的高質(zhì)量測序數(shù)據(jù)比對到海豚鏈球菌參考基因組SF1（GenBank 登錄號CP005941.1）上。采用FPKM 法計算每個注釋基因的表達水平，使用FDR 與log2FC 來篩選差異基因，篩選條件為FDR＜0.05 且|log2FC|＞1。

1.2.5 聚類分析根據(jù)每組篩選的差異表達基因，使用R 語言Pheatmap 軟件包對差異基因進行雙向聚類分析。首先，得到的差異表達基因在GO（The Gene Ontology）數(shù)據(jù)庫中的各個項目進行映射，然后通過GO 功能顯著性富集分析[25]，在宏觀上了解該物種的基因功能分布特征。然后，基于KEGG（Kyoto Encyclopedia of Genes and Genomes）數(shù)據(jù)庫，對每個KEGG pathway 不同層次上存在的差異表達基因數(shù)量進行統(tǒng)計，確定這些差異表達基因的主要代謝路徑和信號通路[26]，進一步了解基因的生物功能。最后，利用超幾何學(xué)方法，以Qvalue≤0.1 的Pathway 作為在差異表達基因中顯著富集的Pathway。

1.2.6 差異表達目標基因的qRT-PCR 驗證為了驗證海豚鏈球菌轉(zhuǎn)錄組數(shù)據(jù)的可靠性，將差異表達基因與毒力因子數(shù)據(jù)庫（Virulence Factors of Pathogenic Bacteria，VFDB）進行比對，篩選出22 個與毒力相關(guān)的差異表達基因，使用primer 5.0 軟件設(shè)計特異性定量引物（基因名稱和引物序列見表1），引物委托生工生物工程（上海）股份有限公司合成。根據(jù)cDNA 反轉(zhuǎn)錄試劑盒的使用說明書將質(zhì)檢合格的總RNA 反轉(zhuǎn)錄為cDNA，應(yīng)用實時熒光定量PCR（Real-time quantitative PCR,qPCR）技術(shù)對其進行轉(zhuǎn)錄水平的驗證，qPCR 反應(yīng)程序為95 ℃預(yù)變性1 min，95 ℃變性15 s，60 ℃退火20 s，72 ℃延伸30 s，共計40 個循環(huán)結(jié)束，60 ℃讀取熔解曲線。以16S rRNA 基因為內(nèi)參基因，采用2-ΔΔCt方法分析計算各毒力基因的相對表達量。

表1 檢測海豚鏈球菌毒力基因的qRT-PCR 特異性引物

2 結(jié)果與分析

2.1 不同溫度下海豚鏈球菌的生長曲線海豚鏈球菌在本試驗的2 個培養(yǎng)溫度下，35 ℃生長速度明顯快于25 ℃，能夠快速進入生長對數(shù)期，35 ℃在3 h 左右即達到對數(shù)生長期，25 ℃在9 h 左右才達到對數(shù)生長期。到達生長平臺期后，35 ℃的吸光值略高于25 ℃的吸光值。當(dāng)吸光值OD600在1.0 左右時，2 個培養(yǎng)溫度的海豚鏈球菌均處于對數(shù)生長期（圖1）。

然而目前職業(yè)學(xué)校，很多家長的文化程度本身較低，對孩子的教育重視程度不夠。而且很多家長忙于生計，工作時間較長，經(jīng)常會說“我不會教育”，“我管不了孩子”，“我沒時間”。對于這樣的家長，我覺得必須幫助他們認識到一個問題，那就是任何事業(yè)的成功，都不能彌補孩子教育的失敗。忙，不是借口；不會，不是理由?！盎畹嚼?，學(xué)到老”，父母也應(yīng)該學(xué)習(xí)在當(dāng)今時代如何與孩子溝通，如何陪伴孩子的成長。

圖1 海豚鏈球菌在不同培養(yǎng)溫度下的生長曲線

2.2 人工感染試驗結(jié)果不同溫度人工感染結(jié)果顯示（圖2），35 ℃時，在感染后1～3 d 卵形鯧鲹即出現(xiàn)大量死亡，之后死亡量逐漸趨于穩(wěn)定，第5 天停止死亡，患病卵形鯧鲹死亡前出現(xiàn)游泳異常、眼球突出、鰭條基部充血等癥狀，病原菌經(jīng)分離培養(yǎng)，16S rRNA 鑒定為海豚鏈球菌。25 ℃時，卵形鯧鲹在感染后3 天內(nèi)未出現(xiàn)死亡，第4 天開始出現(xiàn)少量死亡，死亡率較低，第7 天停止死亡。海豚鏈球菌在水溫35 ℃和25 ℃造成的卵形鯧鲹的累積死亡率分別為81.11%和12.23%，PBS 對照組均未出現(xiàn)死亡。

圖2 環(huán)境溫度對卵形鯧鲹海豚鏈球菌病死亡率的影響

2.3 總RNA 的提取提取培養(yǎng)至對數(shù)生長期（OD600=0.9～1.1）的海豚鏈球菌的總RNA，1%瓊脂糖凝膠電泳結(jié)果如圖3 所示，28 S 和18 S 的條帶清晰無彌散，RNA 完整性好，無降解。核酸蛋白檢測儀測定結(jié)果顯示樣品的A260/A280比值在2.0 左右，濃度為400～700 ng·μL-1，RNA質(zhì)量符合轉(zhuǎn)錄組測序的要求。

圖3 1%瓊脂糖凝膠電泳檢測海豚鏈球菌Tozj-1 株總RNA

2.4 樣品相關(guān)性分析本試驗中每個溫度設(shè)有3 個生物學(xué)重復(fù)，其中CK-1、CK-2、CK-3 分別表示25 ℃的3 個平行樣品，Test-1、Test-2、Test-3 分別表示35 ℃的3 個平行樣品。如圖4 所示，同一溫度的3 個平行樣品相關(guān)性系數(shù)在0.998 0以上，表明來自3 個重復(fù)性試驗的轉(zhuǎn)錄組數(shù)據(jù)具有高度可重復(fù)性，可用于后續(xù)分析。

圖4 樣品相關(guān)性分析熱圖

2.5 數(shù)據(jù)過濾與質(zhì)量評估Illumina 測序得到的數(shù)據(jù)經(jīng)過質(zhì)量檢測和數(shù)據(jù)過濾，最終獲得Clean reads，各樣本的Clean reads 與Raw reads 的Q20和Q30 分別達到了97%和93%以上，GC 含量在各樣本中較為一致，為40.19%～41.22%（表2）。以上數(shù)據(jù)均表明轉(zhuǎn)錄組測序的數(shù)據(jù)良好，可以用于后續(xù)進一步分析。

表2 過濾前后堿基信息統(tǒng)計

2.6 差異表達基因的分析使用Bowtie 2 軟件將6 個轉(zhuǎn)錄組數(shù)據(jù)映射到海豚鏈球菌SF1 基因組，結(jié)果顯示，所有樣品數(shù)據(jù)皆可映射到基因組上，并且這些樣本的映射率均在98%以上（表3）。說明該轉(zhuǎn)錄組數(shù)據(jù)可靠，可以用于后續(xù)的試驗分析。差異表達基因結(jié)果如圖5 所示，共篩選到927 個顯著差異表達基因，其中820 個上調(diào)表達，107 個下調(diào)表達。

圖5 差異基因統(tǒng)計火山圖

表3 原始數(shù)據(jù)比對到基因組

2.7 差異表達基因的GO 功能富集分析對獲得的差異表達基因進行GO 注釋分析，共得到1 611 個GO 功能注釋，其中生物學(xué)過程（biological process）959 個、細胞組分（cellular component）446 個和分子功能（molecular function）570 個。如圖6 所示，差異基因注釋在生物學(xué)過程有關(guān)的GO terms 數(shù)量最多，共計12 個，主要有代謝過程（metabolic process）、細胞過程（cellular process）和單一有機過程（single-organism process）等；其次是與分子功能有關(guān)的terms，共計8 個，主要有催化活性（catalytic activity）、結(jié)合（binding）等；與細胞組分有關(guān)的terms 最少，為7 個，主要有細胞（cell）、細胞成分（cell part）和膜（membrane）等。海豚鏈球菌代謝過程、細胞過程、單一有機過程、結(jié)合過程和催化活性這幾個亞類包含的差異基因比較多，它們均受溫度調(diào)控。

圖6 CK 與Test 差異基因GO 分類

2.8 差異表達基因的KEGG 富集分析將2 個樣本的差異基因進行KEGG 富集分析，結(jié)果顯示，820 個上調(diào)表達基因中有289 個基因可以富集到KEGG 的98 條pathways，其中顯著富集的pathways有5 個。這些顯著上調(diào)的基因主要與信號轉(zhuǎn)導(dǎo)、跨膜轉(zhuǎn)運和毒力因子表達有關(guān)，包括ABC 轉(zhuǎn)運、核糖體通路、群體感應(yīng)、脂肪酸代謝和合成通路（圖7-a）。除此之外，還包括少量與葉酸和淀粉蔗糖代謝的相關(guān)基因。值得注意的是，這些顯著上調(diào)的pathways 中有許多重復(fù)出現(xiàn)的基因（表4），特別是srtF基因（K710_0263）在ABC 轉(zhuǎn)運、群體感應(yīng)系統(tǒng)和雙組份調(diào)控通路中均上調(diào)表達。

圖7 差異表達基因的KEGG 富集分析（35 ℃相對于25 ℃）

表4 同時參與多條通路的上調(diào)表達基因

2.9 轉(zhuǎn)錄組數(shù)據(jù)的qRT-PCR 驗證為了評估轉(zhuǎn)錄組數(shù)據(jù)的可靠性，隨機篩選22 個差異表達的重要毒力基因，采用實時熒光定量PCR 的方法檢測毒力基因的轉(zhuǎn)錄水平，每個基因設(shè)置3 次生物學(xué)重復(fù)，以保證數(shù)據(jù)的準確性，結(jié)果如圖8 所示。2 種檢測方法經(jīng)SPSS 軟件的相關(guān)性分析，相關(guān)系數(shù)達到了0.824，相關(guān)極顯著（P＜0.001），說明qRTPCR 驗證結(jié)果與RNA-Seq 數(shù)據(jù)具有較好的一致性，并且基因表達變化的趨勢基本相同，進一步驗證了轉(zhuǎn)錄組數(shù)據(jù)的可靠性。

圖8 海豚鏈球菌在25 ℃和35 ℃條件下毒力相關(guān)基因的相對表達量

2.10 溫度與毒力因子基因表達的相關(guān)性根據(jù)文獻[27- 28]以及VFDB 數(shù)據(jù)庫的篩選，發(fā)現(xiàn)約58 個當(dāng)前已知或預(yù)測的海豚鏈球菌毒力因子。由表5 可知，上調(diào)表達的毒力因子有46 個，特別是溶血素（SagI）、黏附素（fap1）、轉(zhuǎn)錄調(diào)節(jié)因子（cysB、galR）和C5α 肽酶（scpB）等重要的毒力因子，在35 ℃培養(yǎng)條件下毒力因子基因的表達量與25 ℃培養(yǎng)條件下相比明顯升高（log2FC＞1.6）；下調(diào)表達的毒力因子基因相對較少，僅有12 個（表6），其中包含海豚鏈球菌毒力標志基因simA。

表5 上調(diào)表達的細菌毒力因子

表6 下調(diào)表達的細菌毒力因子

3 討 論

溫度是影響細菌毒力的重要環(huán)境因子。由病原菌引發(fā)的魚類疾病的暴發(fā)往往與水溫有著極為緊密的關(guān)聯(lián)性。郭富強等[29]的研究結(jié)果表明，隨著水溫的不斷升高羅非魚感染無乳鏈球菌的累積死亡率也隨之升高，在33 ℃死亡率最高；劉志剛等[17]發(fā)現(xiàn)無乳鏈球菌感染羅非魚在37 ℃死亡率最高；Zamri-Saad 等[30]認為鏈球菌水溫在32～37 ℃時易出現(xiàn)高發(fā)病率和高死亡率，而水溫26 ℃以下時較少發(fā)病。在本研究中，人工回歸感染試驗結(jié)果表明，卵形鯧鲹感染海豚鏈球菌后在35 ℃下的死亡率高達81.11%，而在25 ℃下的死亡率僅為12.23%，2 個溫度的PBS 對照組均未出現(xiàn)死亡，說明海豚鏈球菌在35 ℃時毒力更強，導(dǎo)致宿主死亡率更高，發(fā)病癥狀更為明顯。同時，本課題組在卵形鯧鲹鏈球菌病的流行病學(xué)調(diào)查研究中發(fā)現(xiàn)，卵形鯧鲹在25 ℃以下幾乎不感染發(fā)病，在30～35 ℃時死亡率較高。因此，本研究設(shè)置25 ℃培養(yǎng)的海豚鏈球菌作為對照組，35 ℃培養(yǎng)的海豚鏈球菌為實驗組，探究其在不同溫度下基因轉(zhuǎn)錄水平上的差異表達，揭示溫度對海豚鏈球菌毒力基因的影響。

在本研究中，當(dāng)設(shè)置FDR＜0.05 且|log2FC|＞1時，與25 ℃培養(yǎng)的海豚鏈球菌相比，35 ℃培養(yǎng)的海豚鏈球菌顯著上調(diào)的基因有820 個，而顯著下調(diào)的基因只有107 個。從整體來看，在較高的培養(yǎng)溫度下（35 ℃）菌株大多數(shù)基因處于上調(diào)狀態(tài)，菌株表現(xiàn)為更加“活躍”的狀態(tài)，而在較低的溫度下（25 ℃）菌株大多數(shù)基因則顯著下調(diào)，菌株顯得相對不那么“活躍”，具體表現(xiàn)為與糖類、碳類的代謝相關(guān)的基因顯著下調(diào)，而與細菌毒力相關(guān)的基因則顯著上調(diào)。Mereghetti[31]等在研究無乳鏈球菌全基因轉(zhuǎn)錄組測序時也有類似的發(fā)現(xiàn)，編碼毒力的相關(guān)基因在40 ℃時上調(diào)表達，而處于生長穩(wěn)定期的大量的新陳代謝相關(guān)基因在30 ℃時上調(diào)表達。

ABC 轉(zhuǎn)運蛋白是已知的最大蛋白家族之一，并在細菌中廣泛存在，其將ATP 水解與各種底物（例如蛋白質(zhì)、脂質(zhì)、肽、固醇、糖、離子及藥物）的主動運輸結(jié)合在一起[32-33]，這將導(dǎo)致細胞質(zhì)膜兩側(cè)的物質(zhì)頻繁交換。在本研究中，海豚鏈球菌在35 ℃時ABC 轉(zhuǎn)運通路相關(guān)的基因顯著上調(diào)，說明35 ℃海豚鏈球菌可能涉及各種物質(zhì)參與的能量交換，以及物質(zhì)的合成與代謝。群體感應(yīng)是一種協(xié)同調(diào)節(jié)系統(tǒng)，它以一種被稱為AIP（autoinduc ing peptide）的自誘導(dǎo)肽作為信號分子，通過感知菌群密度和環(huán)境因子的變化，并把周圍的環(huán)境信息傳遞到雙組份系統(tǒng)當(dāng)中，再由雙組份系統(tǒng)對相關(guān)基因進行表達調(diào)控[34]。一些AIP 的前體肽在核糖體中合成，但是AIP 無法自由地穿透細胞壁，因此必須通過ABC 運輸系統(tǒng)或其他的膜通道蛋白進入胞外發(fā)揮作用。因此部分上調(diào)表達的基因同時參與了ABC 轉(zhuǎn)運通路、雙組份信號系統(tǒng)通路、群體感應(yīng)通路、脂肪酸代謝和合成通路當(dāng)中的2～3 個通路，特別是顯著上調(diào)的srtF（K710_0263）基因同時參與了3 個通路的表達調(diào)控，具有重要的研究意義。研究發(fā)現(xiàn)，srtF屬于C 類sortase 家族，主要作用是在細菌表面形成菌毛并將其固定在細胞壁上，與細胞黏附、定植、生物膜形成等環(huán)境過程相互作用[35]。

在下調(diào)表達基因中比較典型的是雙組份系統(tǒng)中編碼檸檬酸合酶（Citrate lyase）相關(guān)的基因（K710_0369、K710_0374、K710_0375、K710_0379和K710_1345），由此表明，在35 ℃時，Tozj-1菌株的三羧酸循環(huán)被抑制，對葡萄糖的代謝能力明顯降低，而在25 ℃相對較低的環(huán)境溫度中，其對葡萄糖的代謝能力表現(xiàn)良好。這與胡文婷[16]研究的無乳鏈球菌和劉韜[36]研究的魯氏耶爾森菌結(jié)論一致，低溫更有利于病原菌葡萄糖的代謝。并且在35 ℃時，與磷酸轉(zhuǎn)移酶系統(tǒng)、果糖和蔗糖代謝、氨基糖和核苷酸糖代謝相關(guān)的基因顯著下調(diào)，這證明糖類、碳類的代謝更適合在較低溫度下進行。

病原菌的致病過程與毒力因子的共同介導(dǎo)作用密切相關(guān)。海豚鏈球菌幾種重要的毒力因子包括M 蛋白、磷酸葡萄糖苷酶（pgmA）、溶血素S、C5α 肽酶、脫乙?；浮AMP 因子。M 蛋白[37]是一種可以通過莢膜延伸到菌體表面形成菌毛的蛋白，它在細菌粘附、感染、抵御吞噬等方面起關(guān)鍵作用；磷酸葡萄糖苷酶（pgmA）[38]促進細胞壁和莢膜的生物合成，并能抵抗正電性抗菌肽；溶血素S[39- 40]能破壞宿主紅細胞、單核細胞和粒細胞，有可能對腦血管造成損傷；C5α 肽酶能水解中性粒細胞化學(xué)誘導(dǎo)物補體因子C5α，從而破壞宿主細胞的抗侵襲能力；脫乙酰基酶[41]能夠提高病原菌對溶菌酶的耐受性，逃避宿主免疫，黏附并入侵宿主的上皮細胞；CAMP 因子[42]通過與FC 區(qū)域的免疫球蛋白結(jié)合，從而有助于攔截抗體向補體呈遞抗原。這些毒力因子相互作用協(xié)助S.iniae侵入宿主，逃避宿主的免疫防御。在本研究中，不同溫度下培養(yǎng)的海豚鏈球菌，其毒力因子呈差異表達，在35 ℃時，除了溶血素、黏附素、轉(zhuǎn)錄調(diào)節(jié)因子等毒力因子轉(zhuǎn)錄上調(diào)之外，抗逆因子和一些酶類基因的轉(zhuǎn)錄表達量也大幅上升，只有少量的毒力因子下調(diào)表達。劉志剛等[17]、王松艷[21]和郭富強[29]對無乳鏈球菌的研究亦發(fā)現(xiàn)類似的結(jié)果。以上相關(guān)研究結(jié)果均表明，高溫有助于海豚鏈球菌毒力相關(guān)基因表達，導(dǎo)致其毒力增強，有助于其黏附、定植、入侵、胞內(nèi)存活和擴散等致病過程，引起宿主發(fā)病。

本研究結(jié)果初步證實溫度能夠廣泛影響海豚鏈球菌基因的轉(zhuǎn)錄水平，進而調(diào)控海豚鏈球菌的致病性，為進一步研究海豚鏈球菌的致病機理提供理論依據(jù)，但其調(diào)控機制仍需要進一步研究加以闡明。