趙婷 潘寶平

摘要 利用病原微生物革蘭氏陽(yáng)性藤黃微球菌（Micrococcus luteus）和革蘭氏陰性鰻弧菌（Vibrio anguillarum）分別侵染活體青蛤，采用第二代測(cè)序MiSeq技術(shù)測(cè)序，構(gòu)建了青蛤血淋巴中應(yīng)答2種病原物的轉(zhuǎn)錄組文庫(kù)。根據(jù)KEGG等數(shù)據(jù)庫(kù)注釋信息對(duì)Unigene進(jìn)行生物學(xué)通路注釋，預(yù)測(cè)出青蛤免疫信號(hào)通路的差異性表達(dá)基因及相關(guān)注釋信息。結(jié)果發(fā)現(xiàn)，共注釋并篩選出2 605個(gè)差異表達(dá)基因，其中893個(gè)基因注釋到15個(gè)免疫相關(guān)的通路中。該研究可為探索青蛤的免疫識(shí)別方式和信號(hào)傳導(dǎo)路徑提供依據(jù)。

關(guān)鍵詞 青蛤；轉(zhuǎn)錄組文庫(kù)；鰻弧菌；藤黃微球菌；免疫相關(guān)基因

中圖分類號(hào) S917.4 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2018）27-0084-04

Comparison of Immunerelated Genes in Transcriptome Library of Cyclina sinensis Hemolymph under the Stress of Different Pathogenic Bacteria

ZHAO Ting1， PAN Baoping2

（1.Tianjin Medical College， Tianjin 300222；2.Tianjin Key Laboratory of Animal and Plant Resistance， College of Life Sciences， Tianjin Normal University， Tianjin 300387）

Abstract Cyclina sinensis were injected with Micrococcus luteus and Vibrio anguillarum respectively.

Using the Illumina Miseq system of secondgeneration sequencing technology， two transcriptome library of hemolymph was constructed. Unigene could be annotated to biological pathway by KEGG database to predict differential expression network of immune signaling pathway in C. sinensis. The results showed that 2 605 differentially expressed genes were found and a total of 893 genes annotated to 15 immune related pathways. This research could provide basis for exploring immune recognition and signal transduction pathway in C. sinensis.

Key words Cyclina sinensis；Transcriptome library；Vibrio anguillarum；Micrococcus luteus；Immunerelated genes

基金項(xiàng)目 天津市科委應(yīng)用基礎(chǔ)與前沿技術(shù)重點(diǎn)項(xiàng)目（12JCZDJC22800）。

作者簡(jiǎn)介 趙婷（1988—），女，山東煙臺(tái)人，助教，碩士，從事生物化學(xué)與分子生物學(xué)研究。*通訊作者，教授，博士，碩士生導(dǎo)師，從事軟體動(dòng)物分子系統(tǒng)學(xué)、海產(chǎn)貝類養(yǎng)殖與病害防治工作。

收稿日期 2018-06-05

青蛤養(yǎng)殖是我國(guó)傳統(tǒng)的海濱經(jīng)濟(jì)產(chǎn)業(yè)，蘊(yùn)含著巨大的發(fā)展前景[1]。近年來(lái)，由于青蛤養(yǎng)殖產(chǎn)業(yè)中出現(xiàn)的高密度投放、種質(zhì)退化及養(yǎng)殖環(huán)境污染等問(wèn)題，屢次發(fā)生病害及大面積死亡現(xiàn)象，給青蛤養(yǎng)殖業(yè)造成了巨大的經(jīng)濟(jì)損失[2-4]，有關(guān)貝類病害問(wèn)題已成為制約貝類養(yǎng)殖業(yè)持續(xù)發(fā)展的重要因素之一。青蛤的病害防治應(yīng)從外因和內(nèi)因2個(gè)方面入手，外因方面要求改善貝類的養(yǎng)殖環(huán)境，內(nèi)因方面則要求提高貝類自身的抗病能力，即通過(guò)分子免疫學(xué)技術(shù)尋找物種的免疫相關(guān)基因和探索免疫應(yīng)答機(jī)制，最終揭示貝類的免疫調(diào)控分子機(jī)制并為抗病品系選育提供有價(jià)值的數(shù)據(jù)。

目前，有關(guān)貝類轉(zhuǎn)錄組的研究已成為貝類免疫學(xué)研究的熱點(diǎn)，亦成為研究者獲得大量免疫抗病基因的有效途徑之一。Hou等[5]通過(guò)454 高通量測(cè)序方法對(duì)蝦夷扇貝成體各組織和各發(fā)育階段的轉(zhuǎn)錄組進(jìn)行分析，一些與免疫、抗逆性、生長(zhǎng)和繁殖相關(guān)的基因被發(fā)現(xiàn)，為蝦夷扇貝的基因組研究和遺傳提供重要的平臺(tái)。Huan等[6]將文蛤不同發(fā)育階段的幼蟲(chóng)進(jìn)行轉(zhuǎn)錄組測(cè)序，分析了文蛤 4 個(gè)幼蟲(chóng)發(fā)育階段的表達(dá)差異，通過(guò)對(duì)序列拼接及基因功能注釋，發(fā)現(xiàn)與生長(zhǎng)、發(fā)育及免疫相關(guān)的基因。結(jié)合文獻(xiàn)檢索結(jié)果和GenBank中注冊(cè)青蛤EST的情況可知，國(guó)內(nèi)外對(duì)青蛤免疫功能基因的研究剛剛起步[7]，因此構(gòu)建青蛤轉(zhuǎn)錄組文庫(kù)大規(guī)模進(jìn)行測(cè)序是十分必要的。筆者采用第二代測(cè)序技術(shù)對(duì)在不同病原物刺激下的青蛤血淋巴轉(zhuǎn)錄組進(jìn)行了比較，以期為揭示青蛤的免疫應(yīng)答與免疫調(diào)控機(jī)制提供有價(jià)值的數(shù)據(jù)，同時(shí)也為青蛤養(yǎng)殖實(shí)踐中的病害防治策略提供新的研究途徑。

1 材料與方法

1.1 試驗(yàn)材料

從天津大港海灘沿岸采集的青蛤作為試驗(yàn)材料，隨機(jī)選取沒(méi)有損傷、形態(tài)大小一致的青蛤樣品，暫養(yǎng)于通氧氣的海水中，海水密度約1.030 g/cm3，溫度保持在22 ℃左右，海水pH 7.0，投喂5‰的小球藻，7 d后做建庫(kù)試驗(yàn)準(zhǔn)備。

1.2 試驗(yàn)方法

1.2.1 侵染試驗(yàn)。

鰻弧菌在28 ℃條件下用2216E大量培養(yǎng)菌株，藤黃微球菌在37 ℃條件下用LB大量培養(yǎng)菌株，分別進(jìn)行菌液制備。24 h后菌液渾濁，經(jīng)轉(zhuǎn)速3 000 r/min離心10 min，將2種菌體分別用過(guò)濾滅菌的海水清洗3次，去除多余殘留培養(yǎng)基，用同樣海水重新懸浮菌體，OD600控制在04。隨機(jī)抽樣12只健康青蛤，備用于制備轉(zhuǎn)錄組測(cè)序樣品。然后，采用隨機(jī)分組的方式分成鰻弧菌刺激組和藤黃微球菌刺激組。在閉殼肌內(nèi)分別注射等量的菌懸液，分別在注射后6和12 h后取青蛤血淋巴。血淋巴在4 ℃條件下，8 000 r/min離心10 min收集血細(xì)胞，加入1 mL Trizol，置于-80 ℃超低溫冰箱中備用。

1.2.2 青蛤轉(zhuǎn)錄組文庫(kù)的構(gòu)建。

采用Trizol法提取青蛤血淋巴總RNA，產(chǎn)物使用QIAGEN公司的Oligotex mRNA Kits進(jìn)行分離與純化。純化后加入裂解緩沖液將mRNA打斷成短的片段，采用隨機(jī)引物逆轉(zhuǎn)錄合成cDNA。純化的雙鏈cDNA再經(jīng)末端修復(fù)、加A尾、加測(cè)序接頭、連接產(chǎn)物純化、PCR富集接頭連接DNA等步驟完成整個(gè)文庫(kù)的制備。利用Life Tech（Invitrogen）Qubit 2.0熒光計(jì)及生化分析儀（Agilent 2100高靈敏度芯片）檢測(cè)文庫(kù)質(zhì)量，檢測(cè)流程參照其說(shuō)明書(shū)。經(jīng)質(zhì)量檢測(cè)合格的文庫(kù)使用Illumina MiSeq測(cè)序儀上進(jìn)行pair end雙端模式（PE250）測(cè)序。

1.2.3 原始數(shù)據(jù)的預(yù)處理及拼接組裝。

測(cè)序后得到原始圖像數(shù)據(jù)轉(zhuǎn)化成序列數(shù)據(jù)（Raw data 或 Raw reads），原始測(cè)序數(shù)據(jù)中包含接頭信息、低質(zhì)量堿基和未測(cè)出的堿基，通過(guò)深度質(zhì)量控制（QC）過(guò)濾篩選，最終得到的數(shù)據(jù)即為有效數(shù)據(jù)（Clean data或Clean reads）。統(tǒng)計(jì)原始測(cè)序數(shù)據(jù)量以及經(jīng)過(guò)堿基質(zhì)量控制處理后reads平均長(zhǎng)度、reads總數(shù)目和總數(shù)據(jù)量產(chǎn)出等信息。根據(jù)樣本經(jīng)過(guò)質(zhì)量控制分析所保留下來(lái)的測(cè)序reads數(shù)據(jù)進(jìn)行樣本轉(zhuǎn)錄組拼接組裝。利用Trinity生物分析軟件將短的讀長(zhǎng)拼接，得到組裝片段（Contig），將Contig互相拼接，最后得到具有一定長(zhǎng)度的非冗余Unigene。

1.2.4 Unigene Pathway生物通路注釋。

KEGG是系統(tǒng)研究基因的生物代謝途徑和相關(guān)基因功能的大型數(shù)據(jù)庫(kù)，利用KEGG數(shù)據(jù)庫(kù)對(duì)合并拼接的Unigene進(jìn)行Pathway生物學(xué)通路的注釋和預(yù)測(cè)。通過(guò)KEGG pathway注釋，對(duì)免疫相關(guān)基因數(shù)目進(jìn)行統(tǒng)計(jì)。

1.2.5

差異基因的篩選及表達(dá)量分析。

將拼接序列與基因組序列比對(duì)，利用編碼RNA表達(dá)計(jì)算，這是通過(guò)統(tǒng)計(jì)方法估計(jì)的表達(dá)量，也用于后續(xù)樣本差異表達(dá)量比較分析。RNA表達(dá)量的計(jì)算采用FPKM計(jì)算度量指標(biāo)。表達(dá)量的計(jì)算有多種標(biāo)準(zhǔn)化方法，采用參考基因區(qū)域內(nèi)匹配read標(biāo)準(zhǔn)化模式（cufflink軟件）計(jì)算表達(dá)量與差異基因表達(dá)。

針對(duì)isoform（即Unigene）以及gene（即components）分別進(jìn)行表達(dá)量計(jì)算。表達(dá)量用PFKM表示，并給出unigene對(duì)應(yīng)NR與Uniprot蛋白數(shù)據(jù)庫(kù)的最佳匹配結(jié)果作為其注釋。表達(dá)差異分析就是組間差異表達(dá)譜分析。差異分析可分別對(duì)基因（gene FPKM）和基因異構(gòu)體（isoform FPKM，即components）層面分別進(jìn)行考慮計(jì)算。差異分析采用edgeR軟件，在R軟件環(huán)境下實(shí)現(xiàn)。表達(dá)量差異倍數(shù)大于2（FPKM≥2）且FDR<10-5作為篩選標(biāo)準(zhǔn)。

1.2.6 Unigene差異表達(dá)基因GO功能分析。

分別對(duì)Unigene進(jìn)行g(shù)ene以及isoform差異結(jié)果的GO富集度分析。GO功能聚類常用于對(duì)目標(biāo)組基因進(jìn)行功能注釋與分類。全部或差異mRNA/基因功能聚類（富集度）分析是采用GO注釋系統(tǒng)，功能注釋針對(duì)全部差異mRNA結(jié)果展開(kāi)。富集度分析采用超幾何分布算法和并采用多重假設(shè)檢驗(yàn)校驗(yàn)超幾何分布統(tǒng)計(jì)量的P值，獲得校正后P值，即q值。

2 結(jié)果與分析

2.1 測(cè)序數(shù)據(jù)量統(tǒng)計(jì)

對(duì)青蛤被藤黃微球菌刺激的樣本（樣本T）進(jìn)行測(cè)序，共獲得12 916 040條Raw reads，去除低質(zhì)量的和載體序列等，共10 966 322條Clean reads（表1）。青蛤被鰻弧菌刺激的樣本（樣本M）測(cè)序，共獲得15 930 210條Raw reads，去除低質(zhì)量的和載體序列等，共13 389 388條Clean reads（表2）。樣本T的測(cè)序原始數(shù)據(jù)量為3.22 G，平均Q20為88.62 %，GC含量平均值為43.10%；樣本M的測(cè)序原始數(shù)據(jù)量為3.98 G，平均Q20為87.72 %，GC含量平均值為42.67%。轉(zhuǎn)錄組測(cè)序質(zhì)量高，建庫(kù)成功。

2.2 轉(zhuǎn)錄組測(cè)序數(shù)據(jù)的組裝拼接效果統(tǒng)計(jì)

M樣本轉(zhuǎn)錄組序列組裝拼接完成的總轉(zhuǎn)錄本數(shù)量為148 238條，拼接完成的總基因數(shù)量為62 152條，轉(zhuǎn)錄本長(zhǎng)度均值（average contig）約618.50 bp。T樣本轉(zhuǎn)錄組序列組裝拼接完成的總轉(zhuǎn)錄本數(shù)量為137 940條，拼接完成的總基因數(shù)量為57 877條，轉(zhuǎn)錄本長(zhǎng)度平均值（average contig）約644.32 bp。

2.3 免疫相關(guān)基因數(shù)目統(tǒng)計(jì)

通過(guò)KEGG pathway注釋（E值≤10-5作為篩選標(biāo)準(zhǔn)），共893個(gè)基因注釋到15個(gè)免疫相關(guān)pathway，T細(xì)胞受體信號(hào)通路的免疫防御基因居多，其次是趨化因子信號(hào)通路，白細(xì)胞跨內(nèi)皮遷移，F(xiàn)cγR介導(dǎo)的吞噬作用及Toll樣受體信號(hào)通路中也包含大量免疫防御相關(guān)基因（表3）。經(jīng)序列同源性比較、功能基因注釋和代謝途徑分析，此次轉(zhuǎn)錄組測(cè)序中發(fā)掘青蛤免疫防御系統(tǒng)的相關(guān)抗病因子，主要有熱冷休克蛋白家族、親環(huán)素、非消化性的酶和酶抑制劑家族、抗菌肽家族、Toll樣受體家族及其下游信號(hào)因子等。

2.4 Unigene差異表達(dá)分析

由圖1可知，青蛤經(jīng)過(guò)革蘭氏陽(yáng)性菌藤黃微球菌和革蘭氏陰性菌鰻弧菌的分別刺激，共得到2 605個(gè)差異表達(dá)基因，其中1 091個(gè)基因樣品M中上調(diào)表達(dá)（樣品T中下調(diào)）和1 514個(gè)基因在樣品T中上調(diào)表達(dá)（樣品M中下調(diào)）。

2.5 差異表達(dá)基因Pathway功能富集分布

由圖2可知，共1 942個(gè)差異基因被注釋到KEGG并顯著富集到110個(gè)pathway。差異基因主要集中在以下通路：亨廷頓氏病、帕金森病、非酒精性脂肪性肝病、抗原加工提呈、類風(fēng)濕性關(guān)節(jié)炎。差異基因被富集到110個(gè)pathway中，有11個(gè)pathway與免疫相關(guān)（該研究中有15個(gè)與免疫相關(guān)pathway被注釋到）。它們分別為造血細(xì)胞譜系、補(bǔ)體系統(tǒng)、Toll樣受體信號(hào)通路、RIG-I樣受體信號(hào)通路、自然殺傷細(xì)胞介導(dǎo)細(xì)胞毒性、抗原加工提呈、B細(xì)胞受體信號(hào)通路、FcεRI信號(hào)通路、FcγR-介導(dǎo)吞噬、白細(xì)胞跨內(nèi)皮遷移、趨化因子信號(hào)通路。

3 討論與結(jié)論

青蛤是具有高經(jīng)濟(jì)價(jià)值的海洋經(jīng)濟(jì)動(dòng)物，近年來(lái)研究表明，引發(fā)青蛤養(yǎng)殖大規(guī)模病害和死亡的最主要因素是病原菌

的感染[8]。目前研究表明微球菌可以引發(fā)水產(chǎn)動(dòng)物發(fā)生疾病，如變異微球菌（Micrococcus varians）是牙鲆出血病的病原菌，患病牙鲆會(huì)表現(xiàn)出上下頜出血發(fā)紅，體表彌散性出血。與水產(chǎn)動(dòng)物病害相關(guān)的革蘭氏陽(yáng)性菌廣泛分布并造成大面積病害，如藤黃微球菌（Micrococcus luteus）分散在水、土壤、空氣、植物和動(dòng)物表面，一般會(huì)引起傷口等局部組織的感染，也可能會(huì)引起嚴(yán)重感染[9-10]。彭彬等[11]從患出血性疾病的黃鱔體內(nèi)首次分離到藤黃微球菌，這也是首次在水生動(dòng)物體內(nèi)分離到病原菌藤黃微球菌。革蘭氏陽(yáng)性菌生長(zhǎng)過(guò)程中細(xì)胞會(huì)分泌外毒素，會(huì)引起機(jī)體的免疫反應(yīng)，一定時(shí)間內(nèi)大量免疫活性物質(zhì)活性將在一定程度上升。鰻弧菌（Vibrio anguillarum）是青蛤疾病的重要病原微生物，鰻弧菌的致病性與其自身的一些有毒物質(zhì)有關(guān)[12]，在脂多糖、胞外蛋白酶、溶血素等許多毒力因子的共同作用下，引發(fā)青蛤大規(guī)模病害甚至死亡的現(xiàn)象[13]?；【闹虏∵^(guò)程首先通過(guò)吸附宿主細(xì)胞，侵入細(xì)胞中，在體內(nèi)大量進(jìn)行增殖，發(fā)出有毒的物質(zhì)進(jìn)而導(dǎo)致宿主死亡[14]。筆者所在課題組前期研究表明一定濃度的鰻弧菌脅迫青蛤后，其會(huì)產(chǎn)生大量免疫相關(guān)的活性物質(zhì)，相關(guān)免疫基因的表達(dá)也會(huì)發(fā)生一定程度的上調(diào)[15-17]。

為深入探索導(dǎo)致貝類病害和死亡的主要病原細(xì)菌及動(dòng)物的抗逆機(jī)制，國(guó)內(nèi)外對(duì)于貝類轉(zhuǎn)錄組的研究已有報(bào)道[18]，包括對(duì)多種經(jīng)濟(jì)貝類進(jìn)行了高通量測(cè)序，均獲得了大量的EST 序列。貝類免疫系統(tǒng)中的免疫防御相關(guān)的抗病因子，包括免疫識(shí)別信號(hào)Toll樣受體因子及其免疫信號(hào)通路以及下游的熱休克蛋白、親環(huán)素、非消化性的酶和酶抑制劑、抗菌肽等，為探索貝類免疫應(yīng)答方式、免疫識(shí)別防御系統(tǒng)及其調(diào)控網(wǎng)絡(luò)提供了大量有價(jià)值的實(shí)驗(yàn)數(shù)據(jù)。該研究利用第二代高通量測(cè)序技術(shù)對(duì)青蛤的血淋巴細(xì)胞的轉(zhuǎn)錄組文庫(kù)進(jìn)行建庫(kù)及注釋分析，在確保獲得的轉(zhuǎn)錄組拼接結(jié)果良好的條件下，共篩選到893個(gè)免疫相關(guān)基因，并將其注釋到15個(gè)免疫相關(guān)的信號(hào)通路中。綜合序列同源性比較和功能基因注釋信息，發(fā)現(xiàn)青蛤多種免疫防御相關(guān)抗病基因，特別是在Toll樣受體信號(hào)通路中，MyD88樣轉(zhuǎn)接蛋白（Mal）、Trif 相關(guān)轉(zhuǎn)接分子（TRAM）、泛素連接酶（TRAF6），轉(zhuǎn)錄因子 NF-κB及抑制物 IκB等免疫相關(guān)基因在此次青蛤轉(zhuǎn)錄組測(cè)序中首次被拼接和注釋到，為后續(xù)研究青蛤Toll信號(hào)通路信號(hào)傳導(dǎo)機(jī)制奠定重要基礎(chǔ)[19]。

在水產(chǎn)動(dòng)物的免疫研究中，采用外界病原生物刺激的方法獲得大量差異表達(dá)基因，這為后期綜合分析其免疫應(yīng)答方式和免疫防御提供數(shù)據(jù)支持。Astrofsky等[20]將白斑病毒感染后的藍(lán)蝦與健康藍(lán)蝦進(jìn)行基因表達(dá)差異比較，獲得32個(gè)差異表達(dá)基因，其中某些差異表達(dá)基因與提高藍(lán)蝦抗病能力密切相關(guān)，該研究為選育抗病品種奠定堅(jiān)實(shí)基礎(chǔ)。Huvet等[21]利用亮弧菌分別刺激抗性牡礪和易感性牡礪，通過(guò)基因表達(dá)的比較分析，共獲得150個(gè)差異表達(dá)基因，經(jīng)分析抗性牡礪和易感性牡礪存在免疫防御功能方面存在差異。該研究利用革蘭氏陽(yáng)性菌和陰性菌刺激青蛤后，從構(gòu)建的轉(zhuǎn)錄組文庫(kù)中拼接獲得了2 605個(gè)差異表達(dá)基因。在上述差異表達(dá)基因中，有1 091個(gè)基因在樣品M中上調(diào)表達(dá)（相對(duì)于樣品T），而1 514個(gè)基因在樣品T中上調(diào)表達(dá)（相對(duì)于樣品M）。此外，通過(guò)對(duì)差異基因Pathway功能富集，其中共有1 942個(gè)差異基因顯著富集到110個(gè)通路中，其中有11個(gè)pathway跟免疫相關(guān)。上述試驗(yàn)結(jié)果表明青蛤?yàn)閼?yīng)答不同病原物的入侵，不同免疫因子參與免疫應(yīng)答反應(yīng)，而特定基因的差異性表達(dá)證明其對(duì)于革蘭氏陽(yáng)性菌和陰性菌的識(shí)別和清除作用有一定的差異。綜上所述，該研究結(jié)果為后期研究青蛤的先天性免疫的應(yīng)答機(jī)制以及病害防治提供新途徑。

參考文獻(xiàn)

[1] 王斌，王翔，王莉明.我國(guó)主要養(yǎng)殖貝類微生物性疾病研究進(jìn)展[J].海洋環(huán)境科學(xué)，2002，21（3）：76-80.

[2] 孫國(guó)銘，萬(wàn)夕和，劉培庭，等.通州海區(qū)灘涂青蛤死亡原因的初步分析[J].水產(chǎn)養(yǎng)殖，2004，25（2）：26-27.

[3] 曹華.沿海灘涂青蛤死亡原因初探及對(duì)策[J].科學(xué)養(yǎng)魚(yú)，2004（4）：47-48.

[4] 王蘭萍，耿榮慶，陳燕，等.灘涂養(yǎng)殖青蛤體內(nèi)細(xì)菌菌群的初步分析[J].江蘇農(nóng)業(yè)科學(xué)，2007（2）：152-153.

[5] HOU R，BAO Z M，WANG S，et al.Transcriptome sequencing and de novo analysis for Yesso scallop（Patinopecten yessoensis）using 454 GS FLX[J].PLoS One，2001，6（6）：1-8.

[6] HUAN P，WANG H X，LIU B Z.Transcriptomic analysis of the clam Meretrix meretrix on different larval stages[J].Marine biotechnology，2012，14（1）：69-78.

[7] 高瑋瑋，潘寶平.青蛤（Cyclina sinensis）cDNA文庫(kù)構(gòu)建及免疫相關(guān)因子基因的篩選[J].海洋與湖沼，2013，44（1）：90-94.

[8] 于瑞海，王如才，田傳遠(yuǎn)，等.櫛孔扇貝大面積死亡原因分析及預(yù)防的探討[J].海洋湖沼通報(bào)，1998（3）：69-72.

[9] 東秀珠，蔡妙英.常見(jiàn)細(xì)菌系統(tǒng)鑒定手冊(cè)[M].北京：科學(xué)出版社，2001：117-118.

[10] 郭青，錢(qián)鳴亮，俞開(kāi)康.牙鲆出血性敗血病病原菌的初步研究[J].湛江海洋大學(xué)學(xué)報(bào)，2000，20（1）：9-14.

[11] 彭彬，楊光友，陳曉利，等.黃鱔藤黃微球菌的分離鑒定及藥敏試驗(yàn)[J].上海海洋大學(xué)學(xué)報(bào)，2011，20（3）：406-411.

[12] 陳吉祥，于德華，李筠.魚(yú)類病原鰻弧菌致病相關(guān)因子及分子生物學(xué)[J].海洋科學(xué)，2003，27（8）：11-15.

[13] 葛端陽(yáng)，趙婷，潘寶平.鰻弧菌（Vibrio anguillarum）對(duì)青蛤（Cyclina sinensis）的毒性及半致死濃度研究[J].海洋與湖沼，2012，43（6）：1192-1195.

[14] BOWER S M，MCGLADDERY S E，PRICE I M.Synopsis of infections diseases and parasites of commercially exploited shellfish[J].Annual review of fish diseases，1994，4：1-199.

[15] 潘寶平，宋欣，羅凱婭，等.青蛤（Cyclina sinensis）溶菌酶基因在鰻弧菌（Vibrio anguillarum）刺激下的表達(dá)[J].海洋與湖沼，2010，41（6）：901-906.

[16] 任毅鵬，高晶，潘寶平，等.青蛤（Cyclina sinensis）TLR2 基因的克隆與表達(dá)分析[J].海洋與湖沼，2014，45（5）：1037-1043.

[17] 丁丹，潘寶平，王玉梅，等.青蛤（Cyclina sinensis）AP-1基因的克隆及在鰻弧菌（Vibrio anguillarum）侵染下的表達(dá)分析[J].海洋與湖沼，2018，49（1）：192-197.

[18] 劉敏，葉晟，楊鳳，等.高通量測(cè)序技術(shù)及其在貝類轉(zhuǎn)錄組研究中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué)，2017，45（7）：130-133，192.

[19] BELL J K，MULLEN G E D，LEIFER C A，et al.Leucinerich repeats and pathogen recognition in Tolllike receptors[J].Trends in immunology，2003，24（10）：528-533.

[20] ASTROFSKY K M，ROUX M M，KLIMPEL K R，et al.Isolation of differentially expressed genes from white spot virus（WSV）infected Pacific blue shrimp（Penaeus stylirostris）[J].Archives of virology，2002，147（9）：1799-1812.

[21] HUVET A，HERPIN A，DGREMONT L，et al.The identification of genes from the oyster Crassostrea gigas that are differentially expressed in progeny exhibiting opposed susceptibility to summer mortality[J].Gene，2004，343（1）：211-220.