水產(chǎn)生物技術(shù)發(fā)展戰(zhàn)略研究

陳松林，邵長偉,徐鵬

（1.中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所，山東青島266071；2.青島海洋科學(xué)與技術(shù)國家實驗室海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室，山東青島266273；3.中國水產(chǎn)科學(xué)研究院水產(chǎn)生物技術(shù)中心，北京100141）

水產(chǎn)生物技術(shù)發(fā)展戰(zhàn)略研究

陳松林1,2，邵長偉1,2,徐鵬3

（1.中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所，山東青島266071；2.青島海洋科學(xué)與技術(shù)國家實驗室海洋漁業(yè)科學(xué)與食物產(chǎn)出過程功能實驗室，山東青島266273；3.中國水產(chǎn)科學(xué)研究院水產(chǎn)生物技術(shù)中心，北京100141）

本文從水產(chǎn)養(yǎng)殖生物基因組測序、遺傳連鎖圖譜繪制、重要經(jīng)濟性狀相關(guān)分子標記/基因篩選、基因組編輯、基因組選擇以及細胞培養(yǎng)與種質(zhì)冷凍保存等方面綜合介紹了水產(chǎn)生物技術(shù)的發(fā)展現(xiàn)狀，并深入分析了水產(chǎn)生物技術(shù)研究中存在的主要問題，諸如基因功能分析平臺不完善，抗病與性控育種技術(shù)研究滯后，基因組編輯與全基因組選擇技術(shù)剛剛起步等。同時，圍繞上述主要問題，提出了水產(chǎn)生物技術(shù)亟待突破的關(guān)鍵技術(shù)，并建議“十三五”期間設(shè)立重點研究計劃專項，深入開展水產(chǎn)動物基因組資源開發(fā)與利用、重要經(jīng)濟性狀遺傳解析以及水產(chǎn)生物信息大數(shù)據(jù)平臺構(gòu)建等。

水產(chǎn)；生物技術(shù)；基因

DOI 10.15302/J-SSCAE-2016.03.008

一、前言

我國是世界第一大農(nóng)業(yè)生產(chǎn)與消費國，也是世界水產(chǎn)品生產(chǎn)大國，產(chǎn)量連續(xù)20多年居世界首位?！皣赞r(nóng)為本、農(nóng)以種為先”，種業(yè)作為現(xiàn)代農(nóng)業(yè)發(fā)展的第一產(chǎn)業(yè)要素，是確保主要農(nóng)產(chǎn)品有效供給的

重要物質(zhì)基礎(chǔ)，種業(yè)科技的不斷創(chuàng)新已成為提升養(yǎng)殖業(yè)核心競爭力的源泉。近年來，伴隨著國家一系列關(guān)于“現(xiàn)代種業(yè)”“藍色經(jīng)濟”等戰(zhàn)略性指導(dǎo)方針、規(guī)劃的提出，我國水產(chǎn)種業(yè)發(fā)展迎來了前所未有的重大機遇。

水產(chǎn)生物技術(shù)作為水產(chǎn)種業(yè)發(fā)展的核心技術(shù)，涉及到養(yǎng)殖學(xué)、遺傳學(xué)和育種學(xué)等多個學(xué)科，從群體、個體、細胞到基因等多種水平，從經(jīng)典技術(shù)到組學(xué)技術(shù)等多種技術(shù)，從魚類、蝦蟹類、貝類到藻類等各類生物，可以說是一個龐大的系統(tǒng)工程。特別是隨著基因組技術(shù)的發(fā)展，以水產(chǎn)基因組技術(shù)為核心的水產(chǎn)生物技術(shù)在水產(chǎn)良種培育、苗種繁育、種質(zhì)保存及病害防控等方面展現(xiàn)出巨大的應(yīng)用潛力和廣闊的應(yīng)用前景。現(xiàn)代水產(chǎn)生物技術(shù)的研究內(nèi)容主要包括水產(chǎn)生物全基因組解析和基因資源發(fā)掘技術(shù)，重要經(jīng)濟性狀基因鑒定及功能分析技術(shù)，基因組編輯、基因組選擇和分子設(shè)計技術(shù)，水產(chǎn)基因工程產(chǎn)品的研制與利用技術(shù)，種質(zhì)分子鑒定、評估與保存技術(shù)等。近年來，國內(nèi)外在水產(chǎn)生物全基因組測序和精細圖譜繪制上取得很大進展，獲得了海量的基因資源信息，但在重要經(jīng)濟性狀功能基因鑒定和應(yīng)用、基因組編輯和基因組選擇、水產(chǎn)動物胚胎批量化保存以及無脊椎動物細胞培養(yǎng)等技術(shù)上還存在一些問題，滿足不了水產(chǎn)種業(yè)和養(yǎng)殖業(yè)發(fā)展的需要。我國“十三五”及今后相當(dāng)長一段時間內(nèi)水產(chǎn)養(yǎng)殖業(yè)的可持續(xù)發(fā)展對水產(chǎn)生物技術(shù)具有重大需求，大力發(fā)展以水產(chǎn)基因組技術(shù)為核心的水產(chǎn)生物技術(shù)研發(fā)是我國“十三五”期間的重大研究課題。

二、 水產(chǎn)生物技術(shù)發(fā)展現(xiàn)狀

（一）水產(chǎn)生物全基因組測序和精細圖譜繪制

近5年來，國內(nèi)外科學(xué)家在水產(chǎn)養(yǎng)殖生物全基因組測序和精細圖譜繪制方面取得了一些重大進展。國外相繼完成了大西洋鱈[1]、三棘刺魚[2]、藍鰭金槍魚[3]、七鰓鰻[4]、腔棘魚[5]、虹鱒[6]和羅非魚[7]等水產(chǎn)動物的全基因組測序和精細圖譜繪制。

國內(nèi)，中國科學(xué)院海洋研究所牽頭完成了長牡蠣全基因組測序和精細圖譜繪制，并揭示了潮間帶逆境適應(yīng)的分子機制[8]；中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所牽頭完成了半滑舌鰨全基因組精細圖譜繪制，并揭示了ZW性染色體起源和進化以及適應(yīng)底棲生活的分子機制[9]，這也是世界上第一個完成全基因組精細圖譜繪制的比目魚類；中國水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所牽頭完成了鯉魚全基因組測序和遺傳多樣性分析[10]；浙江海洋大學(xué)牽頭完成了大黃魚全基因組測序和精細圖譜繪制[11]；中國科學(xué)院水生生物研究所牽頭完成了草魚全基因組測序和精細圖譜構(gòu)建[12]；深圳華大基因研究院牽頭完成了彈頭魚全基因組測序[13]；中國科學(xué)院海洋研究所完成了桔黃東方鲀?nèi)蚪M測序[14]。此外，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所完成了海帶全基因組圖譜的繪制[15]。上述水產(chǎn)生物全基因組結(jié)構(gòu)的解析為篩選重要性狀相關(guān)功能基因提供了基因資源。

（二）水產(chǎn)養(yǎng)殖生物重要性狀相關(guān)功能基因發(fā)掘

在水產(chǎn)養(yǎng)殖動物性別決定、抗病免疫、生殖發(fā)育等相關(guān)基因篩選與鑒定方面國內(nèi)外也取得了重要進展。近5年來，國外科學(xué)家在多種魚類上發(fā)現(xiàn)了性別決定基因，例如，青鳉的dmy基因、銀漢魚的amhy基因、河豚的amhr2基因、呂宋青鳉的gsdf基因、虹鱒的sdY 基因[16～20]。

國內(nèi)，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所克隆與表征了半滑舌鰨性別決定和生殖相關(guān)的10多個基因[21～23]，特別是發(fā)現(xiàn)dmrtl基因是半滑舌鰨Z染色體連鎖、精巢發(fā)育必不可少的雄性決定基因[9]。西南大學(xué)克隆了羅非魚雄性性腺發(fā)育必不可少的基因amhr[24]。在抗病免疫相關(guān)基因篩選方面，國內(nèi)多家單位都開展了不同水產(chǎn)養(yǎng)殖生物抗病免疫相關(guān)基因克隆的研究， 其中，代表性的單位包括：中國科學(xué)院水生生物研究所解析了鯽魚干擾素系統(tǒng)關(guān)鍵基因的抗病毒作用機理[25～28]，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所和中國科學(xué)院海洋研究所等單位克隆與表征了大菱鲆stat2、grim-l9、牙鲆akirinl、hepn、clq，半滑舌鰨ghClq、irfl等一批抗病免疫相關(guān)基因[29～34]。

（三）水產(chǎn)養(yǎng)殖生物高密度遺傳連鎖圖譜繪制

近5年來，水產(chǎn)養(yǎng)殖動物高密度遺傳連鎖圖譜的構(gòu)建主要采用單核苷酸多態(tài)性（SNP）標記進行，目前已在多種水產(chǎn)養(yǎng)殖動物構(gòu)建了高密度單核苷酸多態(tài)性遺傳圖譜。 國外，利用5 950個單核苷酸多態(tài)性標記構(gòu)建了大西洋鮭高密度遺傳圖譜[35]；利用

2 226個單核苷酸多態(tài)性標記構(gòu)建了虹鱒高密度遺傳圖譜[36]；利用2 235個單核苷酸多態(tài)性標記構(gòu)建了墨西哥脂鯉高密度遺傳圖譜，標記間平均距離為0.94 cM[37]；構(gòu)建了斑點叉尾高密度單核苷酸多態(tài)性遺傳連鎖圖譜[38]。

國內(nèi)，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所構(gòu)建了半滑舌鰨高密度單核苷酸多態(tài)性遺傳連鎖圖譜，定位單核苷酸多態(tài)性標記12 142個[9]，平均圖距為0.326 cM；構(gòu)建了牙鲆高密度單核苷酸多態(tài)性遺傳連鎖圖譜，平均圖距為0.47 cM[39]；構(gòu)建了大菱鲆高密度單核苷酸多態(tài)性遺傳圖譜，平均圖距為0.4 cM[40]。中國水產(chǎn)科學(xué)研究院黑龍江水產(chǎn)研究所構(gòu)建了鯉魚微衛(wèi)星遺傳連鎖圖譜，共包含1 025個標記[41]。中國科學(xué)院海洋研究所構(gòu)建了海灣扇貝高密度微衛(wèi)星遺傳連鎖圖譜[42]和南美白對蝦高密度單核苷酸多態(tài)性遺傳連鎖圖譜，標記間距達0.7 cM[43]。中國海洋大學(xué)構(gòu)建了扇貝高密度單核苷酸多態(tài)性遺傳連鎖圖譜，平均圖距為0.41 cM[44]。此外，中國海洋大學(xué)還構(gòu)建了海帶高密度單核苷酸多態(tài)性遺傳圖譜，平均圖距為0.362 cM[45]。

（四）水產(chǎn)養(yǎng)殖生物重要性狀相關(guān)分子標記篩選與應(yīng)用

國內(nèi)外近幾年來對水產(chǎn)養(yǎng)殖動物重要性狀相關(guān)分子標記的篩選非常重視，其中歐盟國家對魚類性別特異分子標記篩選投資很大，例如，法國主導(dǎo)啟動了20多種養(yǎng)殖魚類性別特異分子標記篩選的項目研究。美國發(fā)現(xiàn)了斑點叉尾雄性特異分子標記AUEST0678[46]。此外，法國科學(xué)家利用模式藻類水云性別決定區(qū)域的序列，開發(fā)了掌狀海帶、裙帶菜和巨藻的性別特異分子標記[47]。而在抗病相關(guān)分子標記篩選方面，初步篩選出與虹鱒抗細菌性冷水病 (CWD) 的12個相關(guān)標記以及與傳染性造血器官壞死病病毒(IHNV)相關(guān)的19個單核苷酸多態(tài)性標記[48]。發(fā)現(xiàn)大菱鲆與盾纖毛蟲病相關(guān)的微衛(wèi)星標記Sma-USC256[49]；發(fā)現(xiàn)斑節(jié)對蝦抗白斑?。╓SD）相關(guān)的RAPD-SCAR標記1個[50]。

國內(nèi)近幾年來在魚類性別特異分子標記篩選方面取得了重要進展，篩選出半滑舌鰨和黃顙魚等多種魚類性別特異擴增片段長度多態(tài)性（AFLP）標記和微衛(wèi)星標記[51]。其中，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所篩選出半滑舌鰨雌性特異擴增片段長度多態(tài)性標記和微衛(wèi)星標記，發(fā)明了鑒別雄魚、偽雄魚、雌魚和超雌魚遺傳性別的聚合酶鏈式反應(yīng)（PCR）技術(shù)，研制了半滑舌鰨高雌性苗種制種技術(shù)，將苗種的生理雌魚比例提高了20 %以上[9,52,53]。中國科學(xué)院水生生物研究所篩選出黃顙魚雄性特異擴增片段長度多態(tài)性標記，建立了黃顙魚遺傳性別鑒定的聚合酶鏈式反應(yīng)（PCR）方法[54]。浙江海洋大學(xué)篩選出條石鯛雄性特異分子標記，建立了條石鯛遺傳性別鑒定技術(shù)[55]。上海海洋大學(xué)利用熒光原位雜交篩選出海帶雌性特異標記[56]。在水產(chǎn)養(yǎng)殖動物抗病相關(guān)分子標記篩選方面報道較少，其中黃海水產(chǎn)研究所篩選出牙鲆抗鰻弧菌病相關(guān)微衛(wèi)星標記[57]。

（五）全基因組選擇技術(shù)

基因組選擇（GS）是挪威科學(xué)家Meuwissen在21世紀初期提出的方法。這一方法是以覆蓋整個基因組的單核苷酸多態(tài)性標記信息為基礎(chǔ)，能夠鑒定出導(dǎo)致目標性狀表型的所有遺傳變異，可以同時選擇抗病、高產(chǎn)、抗逆等多種經(jīng)濟性狀，達到多性狀選育的分子育種模式。這個方法特別適合于低遺傳力的數(shù)量性狀，例如抗病、抗逆等性狀。全基因組選擇技術(shù)在禽、牛和豬等養(yǎng)殖動物上得到應(yīng)用并產(chǎn)生了很好效果。而在水產(chǎn)養(yǎng)殖動物上，由于全基因組測序的滯后，迄今國際上有關(guān)魚類全基因組選擇育種大多停留在理論研究，投入實際應(yīng)用的還很少。其中，挪威科學(xué)家進行了大西洋鮭感染病原前后存活魚和死亡魚基因表達譜的比較研究，通過模擬計算，預(yù)測經(jīng)6～7代選擇后抗病力有可能提高1倍[58]。美國奧本大學(xué)開發(fā)了250 K的斑點叉尾鮰單核苷酸多態(tài)性標記的基因芯片，推進了全基因組選擇在鯰魚育種中的應(yīng)用[59]。挪威、美國等國家目前啟動了大西洋鮭、虹鱒魚等抗病全基因組選擇的項目研究。

“十二五”期間，我國啟動了半滑舌鰨、大黃魚、櫛孔扇貝和對蝦基于全基因組信息的遺傳選育的課題研究，其中中國海洋大學(xué)建立了扇貝全基因組選擇技術(shù)平臺；中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所完成了半滑舌鰨和牙鲆抗病參考群體的構(gòu)建及全基因組重測序，計算了全基因組單核苷酸多態(tài)性標記的遺傳效應(yīng)和基因組估計育種值（GEBV）；集美大學(xué)開展了大黃魚生長性狀的全基因組選擇技術(shù)

研究。上述進展為建立這些種類的全基因組選擇育種技術(shù)奠定了重要基礎(chǔ)。

（六）水產(chǎn)養(yǎng)殖生物基因組編輯技術(shù)

基因組編輯技術(shù)是近幾年來發(fā)展起來的對基因組進行精確修飾的一種先進技術(shù)，主要包括轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶（TALEN）和CRISPR/Cas9兩項技術(shù)。這些基因組編輯技術(shù)在模式魚類斑馬魚和青鳉上成功建立并得到應(yīng)用[60,61]。 有關(guān)水產(chǎn)養(yǎng)殖動物基因組編輯技術(shù)的成功報道目前還很少。Edvardsen 等[62]采用CRISPR/Cas9技術(shù)對大西洋鮭色素沉積相關(guān)的兩個基因tyrosinase和slc45a2進行了定點修飾，導(dǎo)致突變體在幼魚階段相比野生型魚減少了色素沉積。

西南大學(xué)通過轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶技術(shù)敲除了羅非魚性別相關(guān)的2個基因dmrtl和foxl2，揭示了這兩個基因在羅非魚性別分化過程中的功能[63]。蘇州大學(xué)建立了鯉魚轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶和CRISPR/Cas9基因組編輯技術(shù)，并對sp7、runx2、sppl和mstn基因進行了定點突變，研究了這些基因在鯉魚肌間刺形成中的作用[64]。黃海水產(chǎn)研究所成功建立了海水鲆鰈魚類半滑舌鰨受精卵顯微注射和轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶基因組編輯技術(shù)，采用該技術(shù)成功將雄性決定基因dmrtl進行了突變，并觀察到突變后的遺傳雄魚生長變快，接近正常雌魚[65]。隨著越來越多的水產(chǎn)動物的基因組被測序，它們的基因組功能研究會顯得日益重要?；蚪M編輯是基因功能研究和基因組改造的一個重要手段。在水產(chǎn)領(lǐng)域，基因組編輯技術(shù)還處于研究的初期階段，但其在基因敲除和基因組育種方面已經(jīng)表現(xiàn)出巨大潛力，將對水產(chǎn)動物基因功能研究和遺傳改良產(chǎn)生深遠影響。

（七）水產(chǎn)動物細胞培養(yǎng)

近幾年來，國外在水產(chǎn)養(yǎng)殖動物細胞系建立方面取得一些進展。如印度建立了熱帶觀賞魚Puntius fasciatus 和 Pristolepis fasciata的尾鰭細胞系[66]，建立了南亞野鯪魚鰓組織細胞系等[67]。特別是在水產(chǎn)無脊椎動物細胞培養(yǎng)上進行了一些有益探索，并取得可喜進展。Jayesh 等[68]研制出一種斑節(jié)對蝦細胞的施旺細胞培養(yǎng)液（SCCM），班節(jié)對蝦淋巴和卵巢細胞分別傳了2代。 Mercurio等進行了海膽卵巢組織的原代培養(yǎng)[69]。

近5年來，我國在水產(chǎn)養(yǎng)殖動物細胞系建立上取得一些進展，新建立了多種魚類組織細胞系，陳松林等主編出版了《魚類細胞培養(yǎng)理論與技術(shù)》專著[70]。其中黃海水產(chǎn)研究所建立了3個半滑舌鰨組織細胞系，包括卵巢細胞系[71]， 偽雄魚性腺細胞系[72]以及神經(jīng)膠質(zhì)細胞系[73]；建立了1個云紋石斑魚心臟細胞系[74]。

（八）水產(chǎn)養(yǎng)殖動物種質(zhì)冷凍保存

在水產(chǎn)養(yǎng)殖動物精子冷凍保存方面，國外已進入商業(yè)化應(yīng)用階段。其中精子冷凍保存技術(shù)在大西洋鮭種業(yè)發(fā)展和苗種推廣中起到重要作用，挪威的大西洋鮭苗種繁育實行了采用冷凍精子授精魚卵批量化生產(chǎn)魚苗的生產(chǎn)模式。 美國路易斯安那州立大學(xué)建立了魚類精子高效冷凍保存技術(shù)，3個人每天可以冷凍保存1 000管凍精的任務(wù)。 在魚類胚胎冷凍保存方面，伊朗科學(xué)家采用玻璃化冷凍方法，將波斯鱘魚胚胎在液氮中冷凍保存后復(fù)活，成功獲得冷凍復(fù)活的胚胎，最高復(fù)活率達69 %[75]， 從而證明玻璃化冷凍方法在魚類胚胎超低溫冷凍保存上的可行性和應(yīng)用潛力。

國內(nèi)近幾年來不僅建立了性逆轉(zhuǎn)七帶石斑魚、抗病牙鲆和抗病半滑舌鰨優(yōu)良家系的精子冷凍保存技術(shù)和精子庫，而且在石斑魚冷凍精子產(chǎn)業(yè)化應(yīng)用方面取得較大進展，實現(xiàn)了魚類冷凍精子的產(chǎn)業(yè)化應(yīng)用。在水產(chǎn)動物胚胎冷凍保存方面，自2005年我國首次突破海水魚類牙鲆胚胎玻璃化冷凍技術(shù)，并在–196℃液氮溫區(qū)獲得冷凍復(fù)活且孵化出魚苗的牙鲆胚胎以來，國內(nèi)一些學(xué)者又相繼在其他魚類獲得冷凍復(fù)活的胚胎，例如，采用玻璃化技術(shù)在七帶石斑魚上獲得了14粒冷凍復(fù)活胚胎，其中2個凍胚孵化出魚苗[76]。

三、水產(chǎn)生物技術(shù)研究中存在的主要問題

（一）基因組資源開發(fā)進展迅速，但水產(chǎn)養(yǎng)殖動物基因功能分析平臺尚未建立，滿足不了基因功能研究的需要

近年來我國學(xué)者相繼完成牡蠣、半滑舌鰨、鯉魚、草魚和大黃魚等的全基因組測序，成果陸續(xù)在《Nature》《Nature Genetics》等頂級期刊上發(fā)表。然而，

在后基因組時代，重點將是基因功能驗證和重要性狀的遺傳解析。但絕大多數(shù)水產(chǎn)養(yǎng)殖動物，我國尚未建立轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶和CRISPR/ Cas9等基因組編輯技術(shù)和平臺，嚴重影響了重要性狀關(guān)鍵基因的發(fā)掘及其在遺傳育種和水產(chǎn)養(yǎng)殖業(yè)中的應(yīng)用。

（二）水產(chǎn)養(yǎng)殖動物抗病基礎(chǔ)研究及抗病分子育種研究投入低、進展慢，滿足不了產(chǎn)業(yè)發(fā)展的需求

抗病良種培育是我國水產(chǎn)養(yǎng)殖業(yè)提質(zhì)增效的重要途徑。然而由于對水產(chǎn)養(yǎng)殖動物抗病基因的功能及其調(diào)控機制的認識不足，限制了水產(chǎn)養(yǎng)殖動物抗病分子育種的研究進程，滿足不了水產(chǎn)業(yè)對抗病高產(chǎn)良種的需求。

（三）性別特異分子標記和分子輔助性控技術(shù)研究滿足不了魚類養(yǎng)殖業(yè)發(fā)展的需求

水產(chǎn)動物中雌雄個體在產(chǎn)量性狀方面的性別差異非常普遍，因此性別控制育種和相關(guān)性別決定基礎(chǔ)研究是水產(chǎn)前沿技術(shù)研究的重要內(nèi)容。盡管目前已經(jīng)發(fā)現(xiàn)半滑舌鰨、黃顙魚、羅非魚、圓斑星鰈和條石鯛的性別特異分子標記，并進行了產(chǎn)業(yè)化應(yīng)用，但是，在許多其他水產(chǎn)養(yǎng)殖動物尚未開發(fā)出性別特異分子標記，分子性控育種研究也比較滯后。

（四）基因組選擇剛剛起步，離良種選育應(yīng)用尚有距離

基因組選擇將成為今后水產(chǎn)養(yǎng)殖動物育種的重要手段之一， 但目前只在少數(shù)水產(chǎn)養(yǎng)殖動物上進行了前期研究，距離建立全基因組選擇育種技術(shù)并在產(chǎn)業(yè)中進行應(yīng)用仍然有相當(dāng)距離，需要加大投入，重點突破。

四、水產(chǎn)生物技術(shù)發(fā)展戰(zhàn)略與關(guān)鍵技術(shù)

（一）發(fā)展戰(zhàn)略

以服務(wù)于水產(chǎn)養(yǎng)殖業(yè)創(chuàng)新驅(qū)動戰(zhàn)略，提高現(xiàn)代漁業(yè)科技競爭力為主要目標，加快水產(chǎn)生物技術(shù)的原始創(chuàng)新研究，圍繞水產(chǎn)重要經(jīng)濟物種基因組結(jié)構(gòu)和功能解析、經(jīng)濟性狀調(diào)控網(wǎng)絡(luò)、基因組育種技術(shù)體系以及種質(zhì)資源評價與保護，力爭在理論方面有重大發(fā)現(xiàn)，在關(guān)鍵共性技術(shù)方面有所突破，支撐和引領(lǐng)水產(chǎn)種業(yè)發(fā)展，為建設(shè)水產(chǎn)生物技術(shù)強國、引領(lǐng)世界水產(chǎn)養(yǎng)殖業(yè)發(fā)展、實現(xiàn)我國從水產(chǎn)養(yǎng)殖大國向水產(chǎn)養(yǎng)殖強國邁進做出貢獻。

（二）關(guān)鍵技術(shù)

1. 水產(chǎn)生物高復(fù)雜度基因組的組裝和分析技術(shù)

水產(chǎn)生物高復(fù)雜度基因組的組裝和分析技術(shù)主要包括高雜合度、高重復(fù)序列水產(chǎn)生物的全基因組或性染色體基因組拼接、組裝及生物信息學(xué)分析技術(shù)。

2. 重要水產(chǎn)養(yǎng)殖動物基因組編輯技術(shù)

重要水產(chǎn)養(yǎng)殖動物基因組編輯技術(shù)主要包括水產(chǎn)養(yǎng)殖動物（魚、蝦、蟹、貝等）基因組編輯（轉(zhuǎn)錄激活因子樣效應(yīng)物核酸酶和CRISPR/Cas9）技術(shù)以及基因組編輯平臺。

3. 水產(chǎn)動物重要經(jīng)濟性狀關(guān)鍵基因及分子標記篩選與鑒定技術(shù)

水產(chǎn)動物重要經(jīng)濟性狀關(guān)鍵基因及分子標記篩選與鑒定技術(shù)包括水產(chǎn)養(yǎng)殖動物生長、抗病、性別決定等重要性狀相關(guān)分子標記和關(guān)鍵基因的篩選與鑒定技術(shù)。

4. 水產(chǎn)動物基因組選擇和分子設(shè)計技術(shù)

水產(chǎn)動物基因組選擇和分子設(shè)計技術(shù)主要包括水產(chǎn)養(yǎng)殖動物抗病、抗逆和品質(zhì)性狀的全基因組選擇和分子設(shè)計技術(shù)。

五、水產(chǎn)生物技術(shù)發(fā)展建議

鑒于水產(chǎn)生物技術(shù)在我國水產(chǎn)種業(yè)和水產(chǎn)養(yǎng)殖業(yè)可持續(xù)發(fā)展中的重要意義和應(yīng)用潛力，考慮到國內(nèi)外目前的發(fā)展現(xiàn)狀和存在問題，建議我國圍繞如下幾個方面設(shè)立重點研究計劃專項。

（一）水產(chǎn)養(yǎng)殖生物全基因組結(jié)構(gòu)解析及基因資源的深度挖掘

進行魚、蝦、貝、蟹、藻等重要水產(chǎn)養(yǎng)殖生物全基因組測序和精細圖譜構(gòu)建，分析不同水產(chǎn)生物基因組的結(jié)構(gòu)特征和進化規(guī)律，闡釋不同水產(chǎn)養(yǎng)殖生物基因組的結(jié)構(gòu)和功能，進行基因資源的深度挖掘。

（二）水產(chǎn)養(yǎng)殖生物重要經(jīng)濟性狀的遺傳解析及分子設(shè)計的基礎(chǔ)研究

建立水產(chǎn)養(yǎng)殖生物基因組編輯公共平臺和技術(shù)

體系，篩選重要水產(chǎn)養(yǎng)殖生物生長、抗病、性別、環(huán)境適應(yīng)等性狀關(guān)鍵基因；開展抗病、抗逆和品質(zhì)性狀的全基因組關(guān)聯(lián)分析和基因組選擇及分子設(shè)計的基礎(chǔ)研究；構(gòu)建水產(chǎn)養(yǎng)殖生物重要基因高效重組表達技術(shù)體系，研制基因工程重組蛋白；探索重組表達產(chǎn)物提高水產(chǎn)養(yǎng)殖動物生長速率和抗病力的途徑。

（三）重要水產(chǎn)養(yǎng)殖生物經(jīng)濟性狀形成的表觀遺傳與環(huán)境調(diào)控機制

以水產(chǎn)養(yǎng)殖生物生長、生殖、發(fā)育和免疫等重要經(jīng)濟性狀為目標性狀，通過基因組甲基化分析、組蛋白修飾以及小核糖核酸（RNA）分析等技術(shù)手段解析環(huán)境因素介導(dǎo)的重要經(jīng)濟性狀形成的表觀遺傳調(diào)控機制；建立甲基化修飾因子在水產(chǎn)動物中的檢測技術(shù)體系，篩選與重要經(jīng)濟性狀相關(guān)的表觀標記。

（四）重要水產(chǎn)養(yǎng)殖生物基因信息大數(shù)據(jù)平臺構(gòu)建及應(yīng)用

構(gòu)建以種質(zhì)資源數(shù)據(jù)和基因組數(shù)據(jù)為核心的水產(chǎn)生物基因資源信息大數(shù)據(jù)平臺，收集、儲存和分析我國重要水產(chǎn)養(yǎng)殖物種豐富多樣的種質(zhì)資源和基因資源信息；通過數(shù)據(jù)的加工和大數(shù)據(jù)分析，為水產(chǎn)生物技術(shù)研發(fā)、遺傳育種、病害防治和資源保護等提供必要的數(shù)據(jù)支撐。

[1] Star B, Nederbragt A J, Jentoft S, et al. The genome sequence of Atlantic cod reveals a unique immune system [J]. Nature. 2011; 477(7363): 207–210.

[2] Jones F C, Grabherr M G, Chan Y F, et al. The genomic basis of adaptive evolution in threespine sticklebacks [J]. Nature. 2012; 484(7392): 55–61.

[3] Nakamura Y, Morib K, Saitoh K, et al. Evolutionary changes of multiple visual pigment genes in the complete genome of pacific bluefin tuna [J]. PNAS USA. 2013; 110(27): 11061–11066.

[4] Smith J J, Kuraku S, Holt C, et al. Sequencing of the Sea Lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution [J]. Nat Genet. 2013; 45(4): 415–421.

[5] Amemiya C T, Alfoldi J, Lee A P, et al. The African Coelacanth genome provides insights into tetrapod evolution [J]. Nature. 2013; 496(7445): 311–316.

[6] Berthelot C, Brunet F, Chalopin D, et al. The Rainbow Trout genome provides novel insights into evolution after whole-genome duplication in vertebrates [J]. Nat commun. 2014; 5: 3657.

[7] David B, Catherine E W, Yang I L, et al. The genomic substrate for adaptive radiation in African Cichlid Fish [J]. Nature. 2014; 513(7518): 375–381.

[8] Zhang G, Fang X, Guo X, et al. The Oyster genome reveals stress adaptation and complexity of shell formation [J]. Nature. 2012; 490(7418): 49–54.

[9] Chen S, Zhang G, Shao C, et al. Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle [J]. Nat Genet. 2014; 46(3): 253–260.

[10] Xu P, Zhang X, Wang X, et al. Genome sequence and genetic diversity of the common carp (Cyprinus carpio) [J]. Nat Genet. 2014; 46(11): 1212–1219.

[11] Wu C, Zhang D, Kan M, et al. The draft genome of the Large Yellow Croaker reveals well-developed innate immunity [J]. Nat Commun. 2014; 5: 5227.

[12] Wang Y, Lu Y, Zhang Y, et al. The draft genome of the Grass Carp (Ctenopharyngodon idellus) provides insights into its evolution and vegetarian adaptation [J]. Nat Genet. 2015; 47(8): 625–631.

[13] You X X, Bian C, Zan Q J, et al. Mudskipper genomes provide insights into the terrestrial adaptation of amphibious fishes [J]. Nat Commun. 2014; 5: 5594.

[14] Gao Y, Gao Q, Zhang H, et al. Draft sequencing and analysis of the genome of pufferfish takifugu flavidus [J]. DNA Res. 2014; 21(6): 627–637.

[15] Ye N, Zhang X, Miao M, et al. Saccharina genomes provide novel insight into kelp biology [J]. Nat Commun. 2015; 6: 6986.

[16] Matsuda M, Nagahama Y, Shinomiya A, et al. DMY is a Y-specific DM-domain gene required for male development in the medaka fish [J]. Nature. 2002; 417(6888): 559–563.

[17] Hattori R S, Murai Y, Oura M, et al. A Y-linked anti-mullerian hormone duplication takes over a critical role in sex determinat ion [J]. PNAS USA. 2012; 109(8): 2955–2959.

[18] Myosho T, Otake H, Masuyama H, et al. Tracing the emergence of a novel sex determining gene in medaka, Oryzias luzonensis [J]. Genetics. 2012; 191(1): 163–170.

[19] Kamiya T, Kai W, Tasumi S, et al. A trans-species missense SNP in Amhr2 is associated with sex determination in the tiger pufferfish, takifugu rubripes (Fugu) [J]. PLoS Genet. 2012; 8(7): e1002798.

[20] Yano A, Guyomard R, Nicol B, et al. An immune-related gene evolved into the master sex-determining gene in rainbow trout, Oncorhynchus mykiss [J]. Curr Biol. 2012; 22(15): 1423–1428.

[21] Meng L, Zhu Y, Zhang N, et al. Cloning and characterization of tesk1, a novel spermatogenesis-related gene, in half-smooth tongue sole (Cynoglossus semilaevis) [J]. PLoS One. 2014; 9(10): e107922.

[22] Hu Q, Zhu Y. Liu Y, et al. Cloning and characterization of wnt4a gene and evidence for positive selection in half-smooth tongue sole (Cynoglossus semilaevis) [J]. Sci Rep. 2014; 4: 7167.

[23] Shao C, Liu G, Liu S, et al. Characterization of the cyp19a1a gene from a BAC sequence in half-smooth tongue sole (Cynoglossus semilaevis) and analysis of its conservation among teleosts [J]. Acta Oceanol Sin. 2014; 32(8): 35–43.

[24] Li M H, Sun Y L, Zhao J Y, et al. A tandem duplicate of antimüllerian hormone with a missense SNP on the Y chromosome is essential for male sex determination in nile tilapia, Oreochromis niloticus [J]. PLoS Genet. 2015; 11(11): e1005678.

[25] Liu Y, Zhang Y B, Liu T K, et al. Lineage-specific expansion

of IFIT gene family: an insight into coevolution with IFN gene family [J]. PloS One. 2013; 8: e66859.

[26] Sun F, Zhang Y B, Jiang J, et al. Gig1, a novel antiviral effector involved in fish interferon response [J]. Virology. 2014; 448: 322–332.

[27] Wang B, Zhang Y B, Liu T K, et al. Fish viperin exerts a conserved antiviral function through RLR-triggered IFN signaling pathway [J]. Dev Comp Immunol. 2014; 47(1): 140–149.

[28] Zhang J, Zhang Y B, Wu M, et al. Fish MAVS is involved in RLR pathway-mediated IFN response [J]. Fish Shellfish Immunol. 2014; 41(2): 222–230.

[29] Wang N, Wang X, Yang C, et al. Molecular cloning, subcelluar location and expression profile of signal transducer and activator of transcription 2 (STAT2) from turbot, Scophthalmus maximus [J]. Fish Shellfish Immunol. 2014; 35(4): 1200–1208.

[30] Wang N, Wang X, Yang C, et al. Molecular cloning and multifunctional characterization of GRIM-19 (gene associated with retinoid-interferon-induced mortality 19) homologue from turbot (Scophthalmus maximus) [J]. Dev Comp Immunol. 2014; 43(1): 96–105.

[31] Chen S, Li W, Meng L, et al. Molecular cloning and expression analysis of a hepcidin antimicrobial peptide gene from turbot (Scophthalmus maximus) [J]. Fish Shellfish Immunol. 2014; 22(3): 172–181.

[32] Yang C, Wang X, Zhang B, et al. Screening and analysis of PoAkirin1 and two related genes in response to immunological stimulants in the Japanese flounder (Paralichthys olivaceus) [J]. BMC Mol Biol. 2013; 14: 10.

[33] Zeng Y, Xiang J, Lu Y, et al. SghC1q, a novel C1q family member from half-smooth tongue sole (Cynoglossus semilaevis): identification, expression and analysis of antibacterial and antiviral activities [J]. Dev Comp Immunol. 2015; 48(1): 151–163.

[34] Lu Y, Wang Q, Liu Y, et al. Gene cloning and expression analysis of IRF1 in half-smooth tongue sole (Cynoglossus semilaevis) [J]. Mol Biol Rep. 2014; 41(6): 4093–4101.

[35] Lien S, Gidskehaug L, Moen T, et al. A dense SNP-based linkage map for Atlantic Salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns [J]. BMC Genomics. 2011; 12: 615.

[36] Guyomard R, Boussaha M, Krieg F, et al. A synthetic Rainbow Trout linkage map provides new insights into the salmonid whole genome duplication and the conservation of synteny among teleosts [J]. BMC Genetics. 2012; 13: 15.

[37] Carlson B M, Onusko S W, Gross J B. A high-density linkage map for astyanax mexicanus using genotyping-by-sequencing technology [J]. Genes, Genomes, Genetics (Bethesda). 2014; 5(2): 241–251.

[38] Li Y, Liu S, Qin Z, et al. Construction of a high-density, highresolution genetic map and its integration with BAC-based physical map in channel catfish [J]. DNA Res. 2015; 22(1): 39–52.

[39] Shao C, Niu Y, Pasi R, et al. Genome-wide snp identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus) applied to QTL mapping of vibrio anguillarum disease resistance and comparative genomic analysis [J]. DNA Res. 2015; 22(2): 161–170.

[40] Wang W, Hu Y, Ma Y, et al. High-density genetic linkage mapping in turbot (Scophthalmus maximus) based on SNP markers and major sex-and growth-related regions detection [J]. PLoS One. 2015; 10(3): e0120410.

[41] Zhang X, Zhang Y, Zheng X, et al. A consensus linkage map provides insights on genome character and evolution in common carp (Cyprinus carpio) [J]. Mar Biotechnol. 2013; 15(3): 275–312.

[42] Li H, Liu X, Zhang G. A consensus microsatellite-based linkage map for the hermaphroditic bay scallop (Argopecten irradians) and its application in size-related QTL analysis [J]. PLoS One. 2012; 7(10): e46926.

[43] Yu Y, Zhang X J, Yuan J B, et al. Genome survey and high-density genetic map construction provide genomic and genetic resources for the pacific white shrimp litopenaeus vannamei [J]. Sci Rep. 2015; 5: 15612.

[44] Jiao W, Fu X, Dou J, et al. High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc [J]. DNA Res. 2014; 21(1): 85–101.

[45] Zhang N, Zhang L, Tao Y, et al. Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus [J]. BMC Genomics. 2015; 16(1): 189.

[46] Ninwichia P, Peatman E, Perera D, et al. Identification of a sexlinked marker for channel catfish [J]. Anim Genet. 2011; 43(4): 476–477.

[47] Lipinska A P, Ahmed S, Peters A F, et al. Development of PCR-based markers to determine the sex of kelps [J]. PLoS One. 2015; 10(10): e0140535.

[48] Campbell N R, LaPatra S E, Overturf K, et al. Association mapping of disease resistance traits in rainbow rout using restriction site associated DNA sequencing [J]. Genes, Genomes, Genetics (Bethesda). 2014; 4(12): 2473–2481.

[49] Rodriguez-Ramilo S T, Fernandez J, Toro M A, et al. Uncovering QTL for resistance and survival time to philasterides dicentrarchi in turbot (Scophthalmus maximus) [J]. Anim Genet. 2013; 44(2): 149–157.

[50] Dutta S, Biswas S, Mukherjee K, et al. Identification of RAPDSCAR marker linked to white spot syndrome virus resistance in populations of giant black tiger shrimp, penaeus monodon fabricius [J]. J Fish Dis. 2014; 37(5): 471–480.

[51] 陳松林. 魚類性別控制與細胞工程育種[M]. 北京: 科學(xué)出版社, 2013. Chen S L. Fish sex control and breeding by cell engineering [M]. Beijing: Science Press; 2013.

[52] Chen S L, Li J, Deng S P, et al. Isolation of female specific AFLP markers and identification of genetic sex in half-smooth tongue sole (Cynoglossus semilaevis) [J]. Mar Biotechnd (NY) 2007; 9(2): 272–280.

[53] Chen S L, Ji X S, Shao C W, et al. Induction of mitogynogenetic diploids and identification of WW super-female using sex-specific SSR markers in half-smooth tongue sole (Cynoglossus semilaevis) [J]. Mar Biotechnol. 2012; 14(1): 120–128.

[54] Dan C, Mei J, Wang D, et al. Genetic differentiation and efficient sex-specific marker development of a pair of Y- and X-linked markers in yellow catfish [J]. Int J Biol Sci. 2013; 9(10): 1043–1049.

[55] Xu D D, Lou B, Xu H X, et al. Isolation and characterization of male-specific DNA markers in the rock bream oplegnathus

Fasciatus [J]. Mar Biotechnol. 2013; 15(2): 221–229.

[56] Liu Y, Bi Y, Gu J, et al. Localization of a female-specific marker on the chromosomes of the brown seaweed saccharina japonica using fluorescence in situ hybridization [J]. PLoS One. 2012; 7(11): e48784.

[57] Wang L, Fan C, Liu Y, et al. A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder (Paralichthys olivaceus) by bulked segregant analysis [J]. Mar Biotechnol. 2014; 16(5): 513–521.

[58] Robinson N, Hayes B. Modelling the use of gene expression profiles with selective breeding for improved disease resistance in atlantic salmon (Salmo salar) [J]. Aquaculture. 2008; 285(1): 38–46.

[59] Liu S, Sun L, Li Y, et al. Development of the Catfish 250K SNP array for genome-wide association studies [J]. BMC Res Notes. 2014; 7: 135.

[60] Tatsumi Y, Takeda M, Matsuda M, et al. TALEN-mediated mutagenesis in zebrafish reveals a role for r-spondin 2 in fin ray and vertebral development [J]. FEBS Lett. 2014; 588(24): 4543–4550.

[61] Ansai S, Kinoshita M. Targeted mutagenesis using CRISPR/Cas system in medaka [J]. Biol Open. 2014; 3(5): 362–371.

[62] Edvardsen R B, Leininger S, Kleppe L, et al. Targeted mutagenesis in atlantic salmon (Salmo salar) using the CRISPR/Cas9 system induces complete knockout individuals in the F0 generation [J]. PLoS One. 2014; 9(9): e108622.

[63] Li M H, Yang H H, Li M R, et al. Antagonistic roles of dmrt1 and foxl2 in sex differentiation via estrogen production in tilapia as demonstrated by TALENs [J]. Endocrinology. 2013; 154: 4814–4825.

[64] Zhong Z M, Niu P F, Wang M Y, et al. Targeted disruption of sp7 and myostatin with CRISPR/Cas9 results in severe bone defects and more muscular cells in common carp [J]. Sci Rep. 2016; 6: 22953.

[65] 陳松林, 崔忠凱, 鄭漢其, 等. 一種基于基因組編輯的海水鲆鰈魚類種質(zhì)構(gòu)建方法及應(yīng)用[P]. 201610162019. 5. Chen S L, Cui Z K, Zheng H Q, et al. A genome editing-based breeding method in flatfish and its application [P]. 201610162019. 5.

[66] Swaminathan T R, Basheer V S, Gopalakrishnan A, et al. Establishment of caudal fin cell lines from tropical ornamental fishes Puntius fasciatus and Pristolepis fasciata endemic to the western ghats of India [J]. Acta Tropica. 2013; 28(3): 536–541.

[67] Abdul Majeed S, Nambi K S, Taju G, et al. Development, characterization and application of a new fibroblastic-like cell line from kidney of a freshwater air breathing fish channa striatus (Bloch, 1793) [J]. Acta Tropica. 2013; 127(1): 25–32.

[68] Jayesh P, Jose S, Philip R, et al. A novel medium for the development of in vitro cell culture system from penaeus monodon [J]. Cytotechnology. 2013; 65(3): 307–322.

[69] Mercurio S, Di Benedetto C, Sugni M, et al. Primary cell cultures from sea urchin ovaries: a new experimental tool [J]. In Vitro Cell Dev Biol Anim. 2014; 50(2): 139–145.

[70] 陳松林, 秦啟偉. 魚類細胞培養(yǎng)理論與技術(shù)[M]. 北京: 科學(xué)出版社, 2011. Chen S L, Qin Q W. Fish cell culture theory and technique [M]. Beijing: Science Press Ltd.; 2011.

[71] Sun A, Wang T Z, Wang N. et al. Establishment and characterization of an ovarian cell line from half-smooth tongue sole (Cynoglossus semilaevis) [J]. J Fish Biol. 2015; 86(1): 46–59.

[72] Sun A, Chen S, Gao F, et al. Establishment and characterization of a gonad cell line from half-smooth tongue sole (Cynoglossus semilaevis) pseudomale [J]. Fish Physiol Biochem. 2015; 41(3): 673–683.

[73] Wang T Z, Sun A, Wang N. et al. Establishment and characterization of an astroglial cell line derived from the brain of half-smooth tongue sole (Cynoglossus semilaevis) [J]. Zool Res. 2015; 36 (5): 305–310.

[74] 劉肖峰, 陳松林, 沙珍霞, 等. 云紋石斑魚心臟細胞系的建立與鑒定[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報, 2015, 23(10): 1394–1400. Liu X F, Chen S L, Sha Z X, et al. Establishment and characterization of a new cell line from heart of epinephelus moara [J]. J Agr Biotechnol. 2015; 23(10): 1394–1400.

[75] Keivanloo S, Sudagar M . Feasibility studies on vitrification of persian sturgeon (Acipenser persicus) embryos [J]. J Aquac Res Dev. 2013; 4: 172.

[76] Tian Y, Jiang J, Song L, et al. Effects of cryopreservation on the survival rate of the seven-band grouper (Epinephelus septemfasciatus) embryos [J]. Cryobiology. 2015; 71(3): 499–506.

Development Strategy for Fisheries Biotechnology

Chen Songlin1,2, Shao Changwei1,2, Xu Peng3
(1.Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, Shandong , China; 2. Function Laboratory for Marine Fisheries Science and Food Production Processes, National Laboratory for Marine Science and Technology, Qingdao 266273, Shandong, China; 3. Centre for Aquaculture Biotechnology, Chinese Academy of Fishery Sciences, Beijing 100141, China)

The recent developments in fisheries biotechnology are reviewed in this paper, including aquatic animal genome sequencing, genetic mapping construction, screening of molecular marker or genes for economically important traits, genome editing, genomic selection, and cell culture and germplasm conservation. The main problems addressed in the study of fisheries biotechnology are also analyzed, including the incomplete gene function analysis platform, the hysteresis of breeding techniques on sex control and disease resistance, the ongoing of genome editing and genomic selection and so on. The key technological breakthroughs that are needed to solve these problems are outlined. Meanwhile, this paper proposes a key major research project to be implemented, during the 13th Five Year Plan, in order to develop and utilize aquatic animal genome resources, analyze their commercially important genetic traits, and construct a megadata platform for aquatic organisms.

fisheries; biotechnology; gene

Q78

A

2016-04-27；

2016-05-15

陳松林，中國水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所，研究員，研究方向為水產(chǎn)生物技術(shù)和遺傳育種；E-mail: chensl@ysfri.ac.cn

中國工程院重點咨詢項目 “水產(chǎn)養(yǎng)殖業(yè)十三五規(guī)劃戰(zhàn)略研究”(2014-XZ-19-3)；國家自然科學(xué)基金重點項目(31130057)；山東省泰山學(xué)者攀登計劃專家項目

本刊網(wǎng)址：www.enginsci.cn