周佳崢, 郭永軍, 全 峰, 梁 健, 李永仁, 黃亞冬, 梁 爽

不同水交換量下皺紋盤鮑轉(zhuǎn)錄組測(cè)序與分析

周佳崢1, 郭永軍1, 全 峰1, 梁 健1, 李永仁1, 黃亞冬2, 梁 爽1

(1. 天津農(nóng)學(xué)院水產(chǎn)學(xué)院, 天津 300380; 2. 天津市海升水產(chǎn)養(yǎng)殖有限公司, 天津 300270)

養(yǎng)殖過程的循環(huán)水交換量能夠顯著影響皺紋盤鮑()的生理狀態(tài), 進(jìn)而影響其養(yǎng)殖存活率和生長(zhǎng)率, 然而不同循環(huán)水量對(duì)于皺紋盤鮑的生理狀態(tài)影響的分子機(jī)制尚不明確。本實(shí)驗(yàn)以皺紋盤鮑()為研究對(duì)象, 將鮑分別置于循環(huán)水交換頻率為1次/d、2次/d、4次/d的條件下(編號(hào)為C0、T2和T4組)暫養(yǎng)3 d, 然后利用Illumina Hiseq 2500 技術(shù)測(cè)定不同處理組鮑魚的肌肉組織中轉(zhuǎn)錄組表達(dá)情況。結(jié)果表明, 單次換水組(C0)和多次換水組(T2和T4)間差異表達(dá)的基因個(gè)數(shù)較多, 分別為4 095和3 617個(gè), 而T2與T4組之間的差異表達(dá)基因個(gè)數(shù)則較少, 僅為232個(gè)。對(duì)差異表達(dá)基因進(jìn)行GO和KEGG功能分析顯示, 單次換水組(C0)和多次換水組(T2和T4)間差異表達(dá)的基因主要集中在能量代謝和免疫反應(yīng)相關(guān)通路或基因, 隨機(jī)選取10個(gè)差異基因進(jìn)行熒光定量PCR驗(yàn)證, 定量PCR結(jié)果與轉(zhuǎn)錄譜結(jié)果一致, 證實(shí)了轉(zhuǎn)錄組測(cè)序結(jié)果的可靠性。本研究說明了更高的水交換量下, 鮑的代謝和免疫相關(guān)通路被激活, 可能有助于鮑魚維持較好的生理和免疫狀態(tài), 從而表現(xiàn)出更好的生長(zhǎng)、存活等表型。

皺紋盤鮑; 水交換量; 轉(zhuǎn)錄組; 能量代謝; 免疫

在水產(chǎn)生物養(yǎng)殖過程中, 水交換量對(duì)養(yǎng)殖生物的生理狀態(tài)有顯著影響。通常情況下, 水交換量的增加能夠改善生物的生長(zhǎng)和存活等狀態(tài)[1-2]。例如, 在皺紋盤鮑幼鮑的養(yǎng)殖過程中, 當(dāng)水交換量率由流水、24 h/次、48 h/次和72 h/次逐漸降低時(shí), 幼鮑的存活率和增重率也隨之下降[3]。提高水交換量被認(rèn)為能夠顯著減少養(yǎng)殖生物在攝食、排泄等過程中產(chǎn)生的氨氮、亞硝酸鹽等自污染因子的濃度, 同時(shí)提高溶解氧、透明度等指標(biāo), 從而使水生生物維持一個(gè)較好的生理狀態(tài), 最終改善其生長(zhǎng)、存活情況等[4-5]。例如, 研究換水頻率對(duì)刺參養(yǎng)殖水質(zhì)的影響時(shí)發(fā)現(xiàn), 日換水率由0、10%、20%、30%、100%逐步提高時(shí), 刺參養(yǎng)殖水體中的氨氮和亞硝酸鹽的濃度逐漸降低[2]; 當(dāng)換水頻率為10 d/月時(shí), 則能維持優(yōu)良的刺參池塘養(yǎng)殖環(huán)境[6]。

轉(zhuǎn)錄組分析能夠了解特定條件下生物體內(nèi)不同基因的表達(dá)情況, 從而解析生物在生長(zhǎng)發(fā)育過程中或外界環(huán)境因子影響下, 其體內(nèi)的生長(zhǎng)發(fā)育、抗逆應(yīng)激等相關(guān)基因的反應(yīng)情況[7-9], 目前已被廣泛用于研究水產(chǎn)動(dòng)物的生理功能(如水產(chǎn)動(dòng)物毒理學(xué)、發(fā)育生物學(xué)、系統(tǒng)進(jìn)化生物學(xué)等)及免疫功能(包括病毒病、細(xì)菌病、寄生蟲病免疫生物學(xué))等研究[10-11]。轉(zhuǎn)錄組技術(shù)在研究水體自污染物質(zhì)(氨氮、亞硝酸鹽等)對(duì)水產(chǎn)動(dòng)物的脅迫效應(yīng)時(shí)也有應(yīng)用[12]。例如, 氨氮脅迫下的日本沼蝦轉(zhuǎn)錄組結(jié)果顯示, 脅迫前后有887個(gè)基因呈差異性表達(dá), 且主要富集在16條免疫相關(guān)通路中, 如補(bǔ)體和凝血級(jí)聯(lián)、T細(xì)胞受體信號(hào)通路、Toll樣受體信號(hào)通路與B細(xì)胞受體等通路[13]。

本課題組過去的研究發(fā)現(xiàn), 水交換量的變化能夠顯著影響皺紋盤鮑的成活率、攝食率和耗氧率情況, 具體表現(xiàn)為: 隨著水交換量由1次/d提升到4次/d, 鮑的存活率和攝食量均逐漸提高; 且隨著水交換量的提高, 鮑的耗氧率逐漸降低, 1次換水組的鮑耗氧率(0.47±0.05) mg/(g·h)顯著高于4 次換水組(0.37±0.05) mg/(g·h); 相應(yīng)地, 1次換水組的水體溶解氧量則小于4次換水組。此外, 水交換量提升降低了鮑養(yǎng)殖水體中的氨氮、亞硝酸氮等自污染因子的濃度[14]。然而, 水交換量變化影響鮑生理狀態(tài)的分子機(jī)制尚不明晰。本文擬采用轉(zhuǎn)錄組測(cè)序的手段, 研究皺紋盤鮑在不同水交換量的條件下, 其肌肉組織內(nèi)轉(zhuǎn)錄本表達(dá)變化, 以闡明不同水交換量對(duì)于皺紋盤鮑的生理狀態(tài)影響的分子機(jī)制, 為皺紋盤鮑養(yǎng)殖條件的優(yōu)化提供理論支撐, 為貝類健康生態(tài)養(yǎng)殖的理論奠定基礎(chǔ)。

1 材料與方法

1.1 材料

實(shí)驗(yàn)所用的皺紋盤鮑, 引自福建晉江某養(yǎng)殖公司。殼長(zhǎng)(2.24±0.30) cm, 殼寬(1.44±0.21) cm, 殼高(1.32±0.19) cm。低溫空運(yùn)至天津。實(shí)驗(yàn)前進(jìn)行7 d的暫養(yǎng), 暫養(yǎng)條件為鹽度33、pH 7.8、溫度(24±1) ℃。暫養(yǎng)期間按皺紋盤鮑體質(zhì)量的1%投喂鮑配合飼料。暫養(yǎng)期間每天早晚全量換水, 并及時(shí)清除死鮑以及殘餌糞便, 實(shí)驗(yàn)前24 h暫停投喂餌料。

1.2 方法

1.2.1 不同水交換量處理

試由上述馴養(yǎng)的群體中, 隨機(jī)取180只皺紋盤鮑, 平均分成3組, 平均分配于9個(gè)100 L的PVC養(yǎng)殖箱中, 按換水頻率分為1次每天組、2次每天組、4次每天組(每天全量換水1次、2次、4次, 每組3個(gè)平行), 每個(gè)換水組實(shí)驗(yàn)持續(xù)3 d, 實(shí)驗(yàn)期間無鮑魚死亡。

1.2.2 樣品取樣

換水實(shí)驗(yàn)結(jié)束后, 從每個(gè)換水處理組中的PVC整理養(yǎng)殖箱中各取9只鮑(每個(gè)平行箱3只), 對(duì)其進(jìn)行肌肉組織的取樣: 使用已滅菌并去除RNA酶的剪刀、鑷子剪取100 mg左右的肌肉組織, 放入干凈的凍存管中, 然后迅速放入液氮中速凍, 之后轉(zhuǎn)移至?80 ℃冰箱保存。

1.2.3 RNA提取

對(duì)保存的組織進(jìn)行液氮研磨, 充分研磨后使用Trizol法[15]提取各換水組鮑魚肌肉組織的總RNA。對(duì)每個(gè)換水組的9個(gè)RNA樣品進(jìn)行三三等量混合, 每組建立3個(gè)平行RNA混合樣品。各混樣總RNA利用Aligent 2000、Nanodrop、瓊脂糖凝膠電泳3種方法進(jìn)行濃度及純度檢測(cè)。1次每天組、2次每天組、4次每天組的混合樣本分別標(biāo)記為C0、T2和T4。

1.2.4 轉(zhuǎn)錄組測(cè)序及分析

轉(zhuǎn)錄組的建庫及測(cè)序工作委托北京百邁克生物科技有限公司完成, 每個(gè)換水組的鮑建立3個(gè)平行的轉(zhuǎn)錄組測(cè)序文庫(共建庫9個(gè))并進(jìn)行測(cè)序Illumina 2000雙端建庫測(cè)序。獲得原始測(cè)序序列(Sequenced reads)后, 進(jìn)行以下生物信息學(xué)分析: Clean Data的分析; Unigene的拼接和功能注釋; Unigene功能注釋, 包括與NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam數(shù)據(jù)庫的比對(duì)進(jìn)行; Unigene庫基因結(jié)構(gòu)分析主要包括SSR分析、CDS預(yù)測(cè)和SNP分析; 不同樣品中的基因進(jìn)行基于FPKM的表達(dá)量分析; 基于基因在不同樣品中的表達(dá)量, DEseq識(shí)別鑒定差異表達(dá)基因; 對(duì)差異表達(dá)基因進(jìn)行模式聚類、GO功能注釋以及KEGG富集分析。

1.2.5 熒光定量PCR驗(yàn)證

從轉(zhuǎn)錄組結(jié)果里隨機(jī)選取10個(gè)差異基因進(jìn)行qRT-PCR驗(yàn)證, 根據(jù)所挑選的差異表達(dá)基因的序列, 遵循引物設(shè)計(jì)的原則, 使用引物設(shè)計(jì)軟件Beacon designer 7設(shè)計(jì)兩對(duì)定量引物。經(jīng)標(biāo)準(zhǔn)曲線法測(cè)定引物的擴(kuò)增效率, 保證各定量引物的擴(kuò)增效率類似, 且在95%～105%范圍內(nèi)。以上述反轉(zhuǎn)錄合成的cDNA為模板, β-actin基因?yàn)閮?nèi)參[16], 在熒光定量PCR儀上進(jìn)行定量PCR反應(yīng)。反應(yīng)體系包括: 1 μL的cDNA模板(10倍稀釋), 各1 μL的上下游引物(10pmol L-1), 10 μL的2×DNAmo ColorFlash Master Mix(Thermo, USA), 以及8 μL的滅菌超純水。反應(yīng)程序設(shè)置為: 95 ℃ 7 min, (95 ℃ 20 s, 60 ℃ 60 s)×35, 72 ℃ 10 min, 每個(gè)反應(yīng)后均檢查熔解曲線以確保產(chǎn)物單一性, 每個(gè)反應(yīng)均設(shè)置3個(gè)重復(fù)。

2 結(jié)果與分析

對(duì)9個(gè)樣品進(jìn)行Illumina 2000雙端建庫測(cè)序, 共得到207 435 516條原始測(cè)序序列, 測(cè)序正確率高達(dá)99.9%。對(duì)原始測(cè)序序列進(jìn)行質(zhì)量篩選, 即截除Reads中的測(cè)序接頭、引物序列和低質(zhì)量值數(shù)據(jù)等, 共獲得61.81 Gb Clean data, 其中各樣品GC含量在46.60%～ 47.76%之間, Q30堿基百分比不小于94.39%, 表明測(cè)序數(shù)據(jù)完整性好, 質(zhì)量高(表1)。

表1 樣品測(cè)序數(shù)據(jù)質(zhì)量評(píng)估

將得到的所有Unigene序列分別與七個(gè)數(shù)據(jù)庫對(duì)比, 比對(duì)結(jié)果如表2所示, 可知在COG、GO、KEGG、KOG、Pfam、Swissprot、eggNOG七個(gè)數(shù)據(jù)庫中成功注釋的Unigene數(shù)量分別為6 518、7 750、9 116、12 139、15 772、11 014、16 863, 最終總共獲得20 480個(gè)至少在一個(gè)數(shù)據(jù)庫中有注釋信息的Unigene。其中長(zhǎng)度大于等于1 000的Unigene數(shù)為16 025個(gè), 長(zhǎng)度介于300至1 000之間的Unigene數(shù)為4 455個(gè)。

表2 Unigene注釋統(tǒng)計(jì)

利用皮爾遜相關(guān)系數(shù)(Pearson’s correlation coefficient)對(duì)樣品間相關(guān)性的指標(biāo)評(píng)估。如圖1所示, C0、T2、T4各組的生物學(xué)重復(fù)文庫之間的相關(guān)性均較高(2均大于0.98), 而C0和T2, C0和T4的樣品文庫之間的相關(guān)性則較低(2均小于0.25)。

圖1 樣品的相關(guān)性熱圖

如表3所示, C0和T2、C0和T4、T4和T2的差異表達(dá)基因數(shù)目分別為4 095、3 617、232, 且1次換水組與多次換水組間(C0和T2和C0和T4)的差異表達(dá)基因中, 上調(diào)基因的數(shù)目明顯小于下調(diào)基因的數(shù)目。

表3 差異表達(dá)基因數(shù)目

分別將單次換水組(C0)和多次換水組(T2和T4)的差異表達(dá)基因進(jìn)行GO功能注釋, 按照GO三大類型再次進(jìn)行分類, 共分為52個(gè)功能亞類, 如圖2所示, C0和T2和C0和T4的差異表達(dá)基因所注釋到的GO類群比較相似: 在生物學(xué)過程分類下, 主要富集在細(xì)胞學(xué)過程(cellular process)、定位(localization)、代謝過程(metabolic process)和單生物進(jìn)程(single-organism process)等功能; 在細(xì)胞組分分類下, 主要富集在細(xì)胞(cell)、細(xì)胞部分(cell part)、細(xì)胞膜部分(membrane part)和細(xì)胞膜部分(membrane)等功能; 在分子功能中, 主要富集在催化活性(catalytic activi-ty)、附著(bin-ding)和運(yùn)輸活性(transporter activity)等功能。

分別將單次換水組(C0)和多次換水組(T2和T4)的差異表達(dá)基因進(jìn)行KEGG功能注釋, 以此獲知不同換水組之間的差異基因可能參與的生理生化代謝通路。結(jié)果如圖3所示, C0和T2的差異表達(dá)基因主要顯著富集于淀粉和蔗糖代謝(starch and sucrose meta-bolism), 溶酶體(lysosome), 糖胺聚糖降解(glycosa-minoglycan degradation), 戊糖和葡萄糖醛酸酯的相互轉(zhuǎn)化(pentose and glucuronate interconversions), 氨基酸和核苷酸糖的代謝(amino sugar and nucleotide sugar metabolism)等通路, C0和T4的差異表達(dá)基因主要顯著富集于淀粉和蔗糖代謝(starch and sucrose metabolism), 溶酶體(lysosome), 糖胺聚糖降解(gly-cosa-mi-no-glycan degradation), 半乳糖代謝(galactose meta-bolism), 氨基酸和核苷酸糖代謝(amino sugar and nucleotide sugar metabolism)等通路。

根據(jù)FPKM(表達(dá)量)值, 選取不同水交換量處理組之間差異表達(dá)最顯著的基因, 對(duì)其表達(dá)量及NR注釋信息進(jìn)行分析、篩選和統(tǒng)計(jì), 結(jié)果如表4和表5所示, 在C0和T2差異表達(dá)最顯著基因中, 與生物消化代謝(胰蛋白酶trypsin等)與組織構(gòu)成(透明質(zhì)蛋白perlucin, 黏液相關(guān)蛋白mucus-associated protein, 肌球蛋白myosin heavy chain等)等功能相關(guān)的基因顯著在多次換水組T2和T4中上調(diào)表達(dá), 而在1次換水組C0中表達(dá)量較低; 而透明質(zhì)素樣蛋白(perlucin-like protein)海藻酸裂解酶(alginate lyase)內(nèi)切葡聚糖酶4樣亞型(endoglucanase 4-like isoform X1)等基因則在多次換水組中表達(dá)量較低(表4)。C0和T4差異表達(dá)最顯著基因中, 與組織構(gòu)成(睪蛋白testisin等), 生物消化代謝(胰蛋白酶trypsin等), 信號(hào)轉(zhuǎn)導(dǎo)(wnt信號(hào)通路受體frizzled 1, RyR受體ryanodine receptor-like等)在多次換水組T2和T4中上調(diào)表達(dá), 而在1次換水組C0中表達(dá)量較低; 而透明質(zhì)蛋白(perlucin), 纖維素酶(cellulase), 黏蛋白(mucin-5AC-like isoform X3), 內(nèi)切葡聚糖酶(endo-glucanase 14-like)等基因則在多次換水組中表達(dá)量較低(表5)。

從轉(zhuǎn)錄組數(shù)據(jù)庫隨機(jī)選取10個(gè)差異基因進(jìn)行熒光定量PCR驗(yàn)證。運(yùn)用2–ΔΔCt方法對(duì)不同水交換量組中的差異表達(dá)基因進(jìn)行相對(duì)表達(dá)水平定量分析。結(jié)果表示, 選取的10個(gè)差異表達(dá)基因(Unigene編號(hào)c100088, c100118, c100155, c100277, c100299, c100496, c100761, c100829, c100882)在三個(gè)換水組鮑中的轉(zhuǎn)錄本表達(dá)量(以FPKM的對(duì)數(shù)表示)和熒光定量PCR表達(dá)量(以倍數(shù)變化FC的對(duì)數(shù)表示)一致(圖4), 表明轉(zhuǎn)錄組測(cè)序結(jié)果是可靠的。

3 討論

隨著新一代高通分子測(cè)序技術(shù)的不斷完善和發(fā)展與成熟, 轉(zhuǎn)錄組研究在水產(chǎn)生物生理過程機(jī)制的研究中愈發(fā)常見, 已有大量研究將轉(zhuǎn)錄組學(xué)技術(shù)應(yīng)用于水產(chǎn)生物疾病與免疫學(xué)研究[17-18], 繁殖和發(fā)育研究[19-20], 毒性和壓力過程研究[21-22], 生長(zhǎng)和營養(yǎng)[23-24]等方面。本研究利用高通量測(cè)序平臺(tái), 對(duì)在三種不同水交換量條件下養(yǎng)殖的皺紋盤鮑的肌肉進(jìn)行轉(zhuǎn)錄組測(cè)序, 得到了質(zhì)量良好、結(jié)果可信的轉(zhuǎn)錄組測(cè)序結(jié)果數(shù)據(jù)。比較C0和T2、C0和T4與T2和T4數(shù)據(jù), 分別檢測(cè)到4 095、3 617、232個(gè)差異表達(dá)基因, 說明隨著水交換量的增加, 尤其是由1次/d的水交換量提升至2次/d或4次/d時(shí), 顯著改變鮑中相關(guān)基因的表達(dá)水平, 說明換水量的增加對(duì)鮑的生理生化過程產(chǎn)生了比較明顯的影響, 也印證了使用轉(zhuǎn)錄組的手段研究水交換量對(duì)于鮑的生理狀態(tài)的影響機(jī)制是可行的。

圖3 差異表達(dá)基因京都基因與基因組數(shù)據(jù)庫(KEGG)通路富集散點(diǎn)圖

Fig. 3 Scatter diagram of the differentially expressed genes enriched in the KEGG pathways

表4 1次換水組和2次換水組差異表達(dá)最顯著基因功能和表達(dá)量

表5 1次換水組和4次換水組差異表達(dá)最顯著基因功能和表達(dá)量

圖4 熒光定量聚合酶鏈?zhǔn)椒磻?yīng)(PCR)驗(yàn)證結(jié)果

注: FC: 定量PCR差異倍數(shù); FPKM: 轉(zhuǎn)錄組表達(dá)量

本研究中, 通過對(duì)C0組和T2、T4組的差異表達(dá)基因進(jìn)行GO和KEGG富集分析, 發(fā)現(xiàn)差異表達(dá)基因主要富集到了與消化和代謝相關(guān)的基因通路上, 如淀粉和蔗糖代謝, 糖胺聚糖降解, 戊糖和葡萄糖醛酸酯的相互轉(zhuǎn)化, 氨基酸和核苷酸糖的代謝, 半乳糖代謝等通路。在對(duì)不同換水組之間的差異表達(dá)最顯著的基因進(jìn)行功能分析時(shí), 同樣發(fā)現(xiàn)了許多基因的功能屬于消化、代謝類, 如胰蛋白酶、內(nèi)切葡聚糖酶、纖維素酶等。結(jié)合我們過去的研究, 在投喂顆粒狀飼料的鮑養(yǎng)殖水體中, 1次換水組、2次換水組、3次換水組、4次換水組的水體中的氨氮含量分別為(1.46±0.42) mg/L、(0.93±0.59) mg/L、(0.46± 0.11) mg/L、(0.39±0.26) mg/L[14]。說明了換水量的提高很有可能通過改善水質(zhì), 提升了鮑魚機(jī)體對(duì)于營養(yǎng)物質(zhì)的消化和吸收效率。類似地, 日本沼蝦在應(yīng)對(duì)不同濃度的氨氮脅迫時(shí), 其糖酵解相關(guān)通路、蛋白質(zhì)消化和吸收相關(guān)通路也會(huì)發(fā)生相應(yīng)的變化[13], 說明水生生物在水質(zhì)變化時(shí), 其消化代謝相關(guān)生理過程的活性會(huì)發(fā)生變化, 從而從能量供給的角度幫助機(jī)體應(yīng)對(duì)由水質(zhì)變化帶來的不利影響。

此外, 我們發(fā)現(xiàn)溶酶體相關(guān)通路和基因在不同換水組鮑中有顯著的變化。溶酶體是無脊椎動(dòng)物免疫系統(tǒng)中一個(gè)關(guān)鍵性的酶復(fù)合體, 由多種抗體、受體和酶成分組成, 包括溶酶體酶、凝集素、抗菌肽等體液免疫因子, 它以殺菌、促進(jìn)吞噬等方式參與貝類的免疫防御[25-26]。有研究表明, 皺紋盤鮑體內(nèi)的溶酶體凝集素能夠參與對(duì)外源病原的識(shí)別和消滅過程[27]。雜色鮑的溶酶體酚氧化酶活性在注射副溶血弧菌菌液24 h后顯著增加, 且顯著受到其他外部環(huán)境因子如溫度和溶解氧的影響[28-29]。結(jié)合本實(shí)驗(yàn)的結(jié)果, 說明了不同水交換量不僅能夠影響皺紋盤鮑的消化代謝功能, 而且可能通過影響其溶酶體相關(guān)的免疫系統(tǒng)功能, 從而影響皺紋盤鮑整體的免疫生理狀態(tài), 最終使得鮑表現(xiàn)出更好的生理狀態(tài)。

4 結(jié)論

本文研究了不同循環(huán)水交換量對(duì)皺紋盤鮑()轉(zhuǎn)錄組表達(dá)的影響, 研究結(jié)果證實(shí)了不同水交換頻率下, 鮑魚的肌肉轉(zhuǎn)錄本呈現(xiàn)顯著的差異性表達(dá)。這些差異表達(dá)的基因主要功能為能量代謝或免疫反應(yīng)相關(guān), 預(yù)示更高的水交換量可能對(duì)皺紋盤鮑代謝和免疫相關(guān)通路的激活有利, 未來應(yīng)進(jìn)一步聚焦代表性差異表達(dá)基因的具體生物學(xué)功能。

[1] 馬貴范, 郭文學(xué), 王玲玲, 等. 養(yǎng)殖密度和換水量及頻率對(duì)凡納濱對(duì)蝦生長(zhǎng)的影響[J]. 漁業(yè)現(xiàn)代化, 2018, 45(4): 21-26.

MA Guifan, GUO Wenxue, WANG Lingling, et al. The effects of aquaculture density and water exchange rate and frequency on the growth of[J]. Fishery Modernization, 2018, 45(4): 21-26.

[2] 王光玉, 姜佳惠, 祝楠, 等. 換水率和密度對(duì)刺參生長(zhǎng)和水質(zhì)的影響[J]. 漁業(yè)現(xiàn)代化, 2017, 44(6): 55-61.

WANG Guangyu, JIANG Jiahui, ZHU Nan, et al. Effects of water exchange rate and density on the growth and water quality of sea cucumbers[J]. Fishery Moder-nization, 2017, 44(6): 55-61.

[3] 楊鳳, 雷衍之, 王仁波, 等. 皺紋盤鮑自污染及其對(duì)幼鮑生長(zhǎng)及成活率的影響[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2003, 18(1): 1-6.

YANG Feng, LEI Yanzhi, WANG Renbo, et al. Self- contamination ofand its effect on growth and survival rate of juvenile abalone[J]. Journal of Dalian Fisheries University, 2003, 18(1): 1-6.

[4] 唐黎標(biāo). 水質(zhì)對(duì)水產(chǎn)養(yǎng)殖的影響和對(duì)策[J]. 漁業(yè)致富指南, 2017, 18: 36-37.

TANG Libiao. The impact of water quality on aquaculture and its countermeasures[J]. Guide to Getting Rich in Fisheries, 2017, 18: 36-37.

[5] 張偉. 水產(chǎn)養(yǎng)殖系統(tǒng)的水質(zhì)預(yù)測(cè)預(yù)警的應(yīng)用研究[D]. 大連: 大連海洋大學(xué), 2019.

ZHANG Wei. Application research on water quality prediction and early warning of aquaculture system[D]. Dalian: Dalian Ocean University, 2019.

[6] 盧超超, 孫永軍, 高勤峰, 等. 不同換水頻率對(duì)刺參養(yǎng)殖池塘底泥環(huán)境的影響[J]. 海洋湖沼通報(bào), 2014(4): 59-66.

LU Chaochao, SUN Yongjun, GAO Qinfeng, et al. The impact of different water exchange frequencies on the sediment environment of sea cucumber culture ponds[J]. Transactions of Oceanology and Limnology, 2014(4): 59-66.

[7] CINGOLANI P. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome ofstrain w1118; iso-2; iso-3[J]. Fly, 2012, 6(2): 1-13.

[8] HUSSING C, KAMPMANN M L, MOGENSEN H S, et al. Comparison of techniques for quantification of next- generation sequencing libraries[J]. Forensic ence Inter-national Genetics Supplement Series. 2015, 5(4): e276- e278.

[9] 王旭. 基于比較基因組學(xué)和轉(zhuǎn)錄組學(xué)技術(shù)揭示單增李斯特菌毒力因子的研究[D]. 上海: 上海海洋大學(xué), 2016.

WANG Xu. Research on the virulence factors ofmonocytogenes revealed based on comparative genomics and transcriptomics techniques[D]. Shanghai: Shanghai Ocean University, 2016.

[10] 趙婷. 青蛤()轉(zhuǎn)錄組文庫測(cè)序及免疫相關(guān)序列分析[D]. 天津: 天津師范大學(xué), 2014.

ZHAO Ting. Sequencing and immune-related sequence analysis oftranscriptome library[D]. Tianjin: Tianjin Normal University, 2014.

[11] 黃玉婷. 基于轉(zhuǎn)錄組數(shù)據(jù)的長(zhǎng)牡蠣未折疊蛋白反應(yīng)(UPR)熱脅迫響應(yīng)機(jī)制初步研究[D]. 大連: 大連海洋大學(xué), 2018.

HUANG Yuting. Preliminary study on the heat stress response mechanism of the unfolded protein response (UPR) of Pacific oyster based on transcriptome data[D]. Dalian: Dalian Ocean University, 2018.

[12] 章瓊, 蔣高中, 李冰. 水產(chǎn)動(dòng)物對(duì)氨氮脅迫響應(yīng)的轉(zhuǎn)錄組分析研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué), 2015, 43(3): 227-230.

ZHANG Qiong, JIANG Gaozhong, LI Bing. Transcri-ptome analysis research progress of aquatic animals’ response to ammonia nitrogen stress[J]. Jiangsu Agricultural Sciences, 2015, 43(3): 227-230.

[13] 孫健. 日本沼蝦抗氨氮脅迫轉(zhuǎn)錄組與代謝組分析[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2017.

SUN Jian. Analysis of transcriptome and metabolome ofresistance to ammonia nitrogen stress[D]. Taian: Shandong Agricultural University, 2017.

[14] 陳俊男, 郭永軍, 梁健, 等. 不同養(yǎng)殖方式對(duì)皺紋盤鮑養(yǎng)殖效果及養(yǎng)殖水質(zhì)的影響[J]. 海洋科學(xué), 2020, 44(3): 104-112.

CHEN Junnan, GUO Yongjun, LIANG Jian, et al. The effects of different culture methods on the aquaculture effect and the quality of aquaculture water of[J]. Marine Sciences, 2020, 44(3): 104-112.

[15] HUMMON A B, LIM S R, DIFILIPPANTONIO M J, et al. Isolation and solubilization of proteins after TRIzol extraction of RNA and DNA from patient material following prolonged storage[J]. Biotechniques, 2007, 42(4): 472.

[16] LI J, SUN H, LIU X. Acclimation-dependent expression of heat shock protein 70 in Pacific abalone (Ino) and its acute response to thermal exposure[J]. Chinese Journal of Oceanology & Limnology, 2012, 30(1): 146-151.

[17] HUANG L X, LIU W J, JIANG Q L, et al. Integration of transcriptomic and proteomic approaches reveals the temperature-dependent virulence of[J]. Frontiers in Cellular and Infection Microbiology, 2018, 8: 207.

[18] YE H, LIN Q, LUO H. Applications of transcriptomics and proteomics in understanding fish immunity[J]. Fish & Shellfish Immunology, 2018, 77: 319-327.

[19] CHAPMAN R W, READING B J, SULLIVAN C V, et al. Ovary Transcriptome profiling via artificial intelligence reveals a transcriptomic fingerprint predicting egg quality in striped bass,[J]. PLoS One, 2014, 9(5): e96818.

[20] AVARRE J C, DUGUé R, ALONSO P, et al. Analysis of the black-chinned tilapiareproducing under a wide range of salinities: from RNA-seq to candidate genes.[J]. Molecular Ecology Resources, 2013, 14(1): 139-149.

[21] COLLI-DULA R C, FANG X, MORAGA-AMADOR D, et al. Transcriptome analysis reveals novel insights into the response of low-dose benzo(a)pyrene exposure in male tilapia[J]. Aquatic Toxicology, 2018, 201(3/4): 162-173.

[22] CHUEH T C, HSU L S, KAO C M, et al. Transcriptome analysis of zebrafish embryos exposed to deltamethrin[J]. Environmental Toxicology, 2016, 32(5): 1548-1557.

[23] SUN S, XUAN F, GE X, et al. Identification of differentially expressed genes in hepatopancreas of oriental river prawn,exposed to environmental hypoxia[J]. Gene, 2014, 534(2): 298- 306.

[24] DIEGO R, CARLOS F, MIGUEL H, et al. Integrative transcriptome, genome and quantitative trait loci resources identify single nucleotide polymorphisms in candidate genes for growth traits in turbot[J]. International journal of molecular sciences, 2016, 17(2): 243.

[25] CHANDHINI S, REJISH KUMAR V J. Transcripto-mics in aquaculture: current status and applications[J]. Reviews in Aquaculture, 2019, 11(4): 1379-1397.

[26] 杜麗, 張巍, 陸逵, 等. 貝類免疫機(jī)制研究進(jìn)展[J]. 動(dòng)物醫(yī)學(xué)進(jìn)展, 2008(3): 77-81.

DU Li, ZHANG Wei, LU Kui, et al. Research progress in shellfish immune mechanism[J]. Advances in Veterinary Medicine, 2008(3): 77-81.

[27] WANG N, WHANG I, LEE J. A novel C-type lectin from abalone,, agglutinates[J]. Developmental & Comparative Immunology, 2008, 32(9): 1034-1040.

[28] CHENG W, HSIAO I S, HSU C H, et al. Change in water temperature on the immune response of Taiwan abalonesupertexta and its susceptibility to[J]. Fish & Shellfish Immunology, 2004, 77: 319-327.

[29] CHENG W, LI C H, CHEN J C. Effect of dissolved oxygen on the immune response ofsupertexta and its susceptibility to[J]. Aquaculture, 2004, 232(14): 103-115.

Sequencing and analysis of thetranscriptome under different water exchange rates

ZHOU Jia-zheng1, GUO Yong-jun1, QUAN Feng1, LIANG Jian1, LI Yong-ren1, HUANG Ya-dong2, LIANG Shuang1

(1. Fisheries College, Tianjin University of Agriculture, Tianjin 300380, China; 2. Tianjin Haisheng Aqua-cu-l-ture Limited Company, Tianjin 300270, China)

The rate of circulating water exchange can remarkably affect the physiology of abalone, including survival and growth. However, the molecular mechanism of this effect is unclear. In this study, abalones were placed under a circulating water exchange frequency of 1, 2, or 4 times/day (called C0, T2, and T4, respectively) for 3 days. Illumina Hiseq 2500 technology was used to investigate the transcriptome expression in the muscle tissues of abalone from the different water exchange groups. The results showed that many differentially expressed genes were detected between the single water exchange group (C0) and the multiple water exchange groups (T2 and T4), such as 4, 095 and 3, 617, respectively, while the number of differentially expressed genes between the T2 and T4 groups was relatively small at 232. The GO and KEGG functional analyses showed that the differentially expressed genes between the single water exchange group (C0) and the multiple water exchange group (T2 and T4) were mainly classified into energy metabolism and immune response-related pathways or genes. Ten differentially expressed genes were selected for verification by fluorescence quantitative polymerase chain reaction, and the results were consistent with the transcription profile results, confirming the reliability of the transcriptome sequencing results. This study shows that metabolism and the immune-related pathways of abalone were activated under the higher water exchange rate, which may help abalone maintain better physiological and immune status, thus leading to better growth, survival, and other phenotypes.

; water exchange rate; transcriptome; energy metabolism; immunity

Oct. 3, 2020

S966.9

A

1000-3096(2022)12-0159-10

10.11759/hykx20201003001

2020-10-03;

2020-11-09

天津市種業(yè)科技重大專項(xiàng)(19ZXZYSN00080); 天津市自然科學(xué)基金(18JCQNJC84500); 廣東省水產(chǎn)健康安全養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室開放基金(GDKLHSA0806); 2020年天津市大學(xué)生創(chuàng)新創(chuàng)業(yè)項(xiàng)目(202010061010); 2020年現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)資金(CARS49)

[Tianjin Seed Industry Science and Technology Major Project No. 19ZXZYSN00080; Tianjin Natural Science Foundation, No. 18JCQNJC84500; Guangdong Key Laboratory of Aquatic Health and Safety Aquaculture Open Fund, No. GDKLHSA0806; Tianjin University Student Innovation and Entrepreneurship Project in 2020, No. 202010061010; Modern Special Fund for Construction of Agricultural Industrial Technology System, No. CARS49]

周佳崢(1999—), 男, 本科生, 水產(chǎn)養(yǎng)殖學(xué)專業(yè), E-mail: 3055912470@qq.com; 梁爽(1989—), 男,通信作者, 博士, 講師, 主要從事貝類養(yǎng)殖生理生態(tài)研究, E-mail: shuangliang@tjau.edu.cn

(本文編輯: 楊 悅)