王明珠 孟憲紅 孔 杰 曹家旺王 軍 馮亞萍 強光峰

(1. 農業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室 中國水產科學研究院黃海水產研究所 青島 266071；2. 青島海洋科學與技術國家實驗室海洋漁業(yè)科學與食物產出過程功能實驗室 青島 266071；3. 上海海洋大學水產與生命學院 上海 201306)

中國明對蝦(Fenneropenaeus chinensis)主要分布于我國黃渤海區(qū)，包括山東、河北、遼寧、天津及江蘇近海，在朝鮮半島西海岸和南海岸也有批量生產，是我國重要的漁業(yè)資源。中國明對蝦對水溫的耐受范圍為 4℃～38℃，生長的適溫范圍是 18℃～30℃，以25℃為最宜(王清印, 2013)。中國明對蝦的適宜生長溫度限制了對蝦的養(yǎng)殖季節(jié)和地域，更重要的是對低溫的不耐受使得中國明對蝦親蝦越冬的安全性降低，越冬成本也大幅上升，影響了養(yǎng)殖的產量和效益。因此，有必要對中國明對蝦進行遺傳改良，選育出生長快，同時對低溫耐受度較高的優(yōu)良新品種。

對耐低溫性狀進行遺傳參數評估可為選擇育種提供理論依據和基礎。水產動物中溫度耐受性遺傳參數評估的研究在國內外已開展多年，如王俊等(2011)對哲羅鮭(Hucho taimen)幼魚在不同溫度下體重(Body weight, BW)和體長(Body length, BL)的遺傳力進行了估計，其低溫組(9±1.5)℃和高溫組(13±1.5)℃體重的父系全同胞遺傳力分別為 0.413～0.653和0.396～0.545，體長的父系全同胞遺傳力分別為0.297～0.701 和 0.194～0.608。Charo-Karisa等(2005)對尼羅羅非魚(Oreochromis niloticus)耐低溫性狀進行研究，其低溫累計存活小時數(Cooling degree hours, CDH)和死亡時溫度(Temperature at death, TAD)的遺傳力分別為(0.08±0.19)和(0.09±0.17)，均為低等遺傳力；同時，發(fā)現(xiàn)體重對低溫耐受力存在顯著影響，個體較小的耐受性較差。Cnaani等(2000)對羅非魚雜交品種的耐寒性進行了研究，認為體長與低溫累計存活天數(Coolingdegree days, CDD)相關性很小。目前，對蝦的相關研究較少，僅Li等(2015、2016)對中國明對蝦成蝦和凡納濱對蝦(Litopenaeus vannamei)耐低溫性狀遺傳參數進行了初步研究。對于耐低溫性狀遺傳參數的準確估計以及評估其與生長性狀的遺傳相關，有重要的理論指導意義，能夠直接影響到我們育種策略的制定。

本研究通過室內人工低溫脅迫的方式，對中國明對蝦成蝦在低溫下生長性狀和耐低溫性狀進行遺傳評估，分析了中國明對蝦成蝦在低溫下的生長性狀和耐低溫性狀的遺傳力以及相關性，以期為中國明對蝦耐低溫新品種的選育提供理論依據。

1 材料與方法

1.1 實驗材料

實驗于中國水產科學研究院黃海水產研究所海水養(yǎng)殖遺傳育種中心即墨市鰲山衛(wèi)實驗基地進行。實驗材料為 2015年構建的中國明對蝦“黃海 2號”G10代家系。家系平均體重達到 3 g后，對其進行熒光(Visible implant elastomer, VIE)標記，以區(qū)分家系。將所有對蝦置于 2個 100 m3的室外圓形水泥池中進行混養(yǎng)，混養(yǎng)期間，每日換水量為30%，溫度為自然水溫。每天根據其營養(yǎng)需求投喂2次混合飼料和2次菲律賓蛤仔(Ruditapes philippinarum)肉。本研究共使用69個全同胞家系(包含4個半同胞家系)，1327尾成蝦，每個家系的個體數在10～41尾之間。

1.2 成蝦低溫脅迫處理

由于中國明對蝦的存活耐受低溫為4℃，所以本實驗設置4℃為低溫脅迫水溫。當外界自然水溫降至4℃時，將所有對蝦移入室內6 m×6 m×1 m的水泥池中，并將水溫維持在(4±1)℃。實驗所用池塘為3個。每天吸底、換水1次，投喂蝦料4次。從第1尾死蝦出現(xiàn)開始，每2 h查看1次水池，撈取死亡個體。死亡定義為對蝦失去平衡，倒向一側，并對外界刺激失去反應。記錄死亡個體的死亡時間、家系信息、性別、體重、體長等信息。待全部對蝦死亡達到50%時，停止低溫脅迫。記錄所有存活對蝦的家系信息、性別、體重、體長。低溫脅迫共持續(xù)7 d。

1.3 數據處理

實驗數據用Excel軟件進行整理，描述性統(tǒng)計用SPSS 19.0軟件進行分析。

1.3.1 中國明對蝦成蝦低溫下體重、體長的遺傳力估計

建立線性混合模型，利用平均信息約束極大似然法(Average information restricted maximum likelihood,AIREML)估計中國明對蝦成蝦體重和體長的方差組分。生長數據通過ASReml軟件進行分析。池塘和性別作為固定效應，經檢驗所有的固定效應均具有顯著意義(P＜0.01)。采用2種單性狀個體動物模型對體重、體長的遺傳方差組分進行估計，如下：

式中，模型1比模型2中多加了共同環(huán)境效應，Yijkl和 Yijl表示第 l尾蝦的體重或體長的觀測值，μ是總體平均值，Tanki是第i個池塘的固定效應，Sexj是第j個性別的固定效應，al是第l尾蝦的體重或體長的加性遺傳效應，Ck為第k個全同胞家系單獨養(yǎng)殖時產生的共同環(huán)境效應， eijlk和 eijl是第l尾蝦的隨機殘差。表型方差是所有方差組分的和，計算公式：

遺傳力是加性遺傳方差和表型方差的比值，計算公式：

式中，為表型方差，是加性方差，是共同環(huán)境方差，是殘差組分。

1.3.2 中國明對蝦成蝦耐低溫性狀的遺傳力估計

耐低溫性狀的遺傳參數選用半致死時的存活狀態(tài)(Survival status at half lethal time, SS50)作為指標進行估計。SS50是一個非連續(xù)性狀，因此，用廣義線性混合模型(Generalized linear mixed model, GLMM)方法(連接函數 Probit)估計其方差組分。由于加入共同環(huán)境效應后，模型不收斂，所以，模型中沒有加入共同環(huán)境效應，育種分析模型為父母本模型：

式中，ijkly 是第l尾蝦的存活狀態(tài)，存活的個體記錄為1，死亡個體記錄為0；ijklλ是ijkly的潛在變量；μ是總體平均值；Sexj是第j個性別的固定效應；Sirei和Damk則分別是第i個父本和第 k個母本的加性遺傳效應；ijkle 是第 l尾蝦的隨機殘差。表型方差組分的計算公式為：

遺傳力的計算公式則為：

式中，是表型方差，是父本加性方差，是母本加性方差，是殘差組分。

1.3.3 不同動物模型間的比較 用似然比檢驗(Likelihood ratio test, LRT)對1.3.1中含共同環(huán)境效應和不含共同環(huán)境效應的2種動物模型進行比較。檢驗統(tǒng)計量為:

式中，LR為似然比值，ln()rL 和ln()fL 分別是2個不同模型下的最大似然函數值，模型 r是模型 f的子模型。LR值服從卡方分布，自由度為2個模型中估計的參數個數之差，本研究中即為1。

1.3.4 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀間的相關分析 不同性狀間表型和遺傳相關均利用兩性狀動物模型進行評估。體重和體長之間的表型和遺傳相關估計所用的兩性狀動物模型是模型(2)的擴展，模型中每一項的含義同模型(2)描述一致。體重、體長和SS50之間的表型和遺傳相關分析中，SS50被設置為第1個變量，其模型同模型(5)；體重和體長分別被設置為第2個變量，模型同模型(2)。模型中各項的描述同前。

1.3.5 顯著性檢驗分析 運用 Z-score來檢驗上述各性狀遺傳力和相關評估值是否顯著(Nguyen et al,2007)，公式如下：

式中， xi和 xj是指各性狀的遺傳力或兩性狀間的相關系數值，是指相應遺傳力或相關系數的標準誤。當檢驗遺傳力是否與0有顯著差異時， xj和σj則均定義為0，當檢驗相關系數與1是否差異顯著時，xj和σj分別定義為1和0。

2 結果

2.1 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀的描述性統(tǒng)計

中國明對蝦成蝦低溫下生長和耐低溫性狀的樣本量、平均值、最大值、最小值、標準差以及變異系數如表1所示。從表1可以看出，體重、體長和半致死時各家系存活率(Survival rate for each family at half lethal time, SR50)平均值分別為17.81 g、114.98 mm、48.45%，變異系數分別為26.41%、8.62%和30.01%，其中，SR50的變異系數值最大(30.01%)，表明不同家系間的半致死存活率存在較大變異；體長的變異系數值最小(8.62%)，表明體長在家系間的變異相對較小。

2.2 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀的遺傳參數評估及2種模型間的比較

中國明對蝦成蝦低溫下生長和耐低溫性狀的遺傳參數以及生長性狀 2種模型間的比較結果見表2。由模型 1獲得的體重和體長的遺傳力分別為(0.206±0.177)和(0.187±0.179)，二者均與0無顯著差異，并且二者之間也無顯著差異(Z=0.075, P＞0.05)。由模型 2獲得的體重和體長的遺傳力分別為(0.317±0.065)和(0.298±0.063)，且二者均與0存在極顯著差異(P＜0.01)。經LRT檢驗體重、體長兩性狀在模型1和模型2之間的似然比值分別為1.640和1.764，均不存在顯著差異(P＞0.05)，表明模型1和模型2之間不存在顯著差異，模型2為較優(yōu)遺傳模型。以SS50作為耐低溫指標，遺傳力為0.169±0.078，經Z檢驗與0存在顯著差異(P＜0.05)。

表1 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀的描述性統(tǒng)計Tab.1 Description statistics of growth traits under low temperature and cold tolerance traits in F. chinensis

2.3 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀的相關分析

中國明對蝦成蝦的生長性狀與耐低溫性狀間的兩兩相關性分析見表3。從表3可以看出，體重和體長二者間的表型和遺傳相關系數分別為(0.823±0.010)和(0.969±0.018)，為高度正相關，且二者均與1無顯著差異(P＞0.05)。體重、體長和SS50表型相關系數分別為(0.187±0.030)、(0.218±0.030)，遺傳相關系數分別為(0.517±0.205)、(0.538±0.203)。由以上可知，中國明對蝦的體重、體長和 SS50之間存在較高的遺傳正相關，而表型相關則相對較弱。

3 討論

3.1 中國明對蝦成蝦低溫下生長性狀遺傳力及 2種模型間的比較

在本研究中，通過2種模型獲得的中國明對蝦成蝦在低溫下體重和體長的遺傳力均為中低等遺傳力。由模型 1獲得的體重、體長的遺傳力分別為(0.206±0.177)和(0.187±0.179)，模型 2獲得的體重、體長的遺傳力分別為(0.317±0.065)和(0.298±0.063)，且模型1獲得的遺傳力均與0無顯著差異(P＞0.05)，而模型2獲得的遺傳力均與 0存在顯著差異(P＜0.01)。2種生長性狀不同模型間的似然比值分別為1.640和1.764，均無顯著差異(P＞0.05)，故認為共同環(huán)境效應對模型的影響不大。當模型不加共同環(huán)境效應時，兩性狀的遺傳力較高一些，但加入共同環(huán)境效應后，兩性狀遺傳力都較低，且標準誤較大，這可能是由于本研究中69個全同胞家系中僅包含了 4個半同胞家系，家系間缺乏遺傳聯(lián)系，加性遺傳效應和共同環(huán)境效應存在一定程度的混淆，使得加性遺傳效應剖分不夠準確(Castillo-Juárez et al, 2007; Luan et al, 2015)。此前，Li等(2016)也用2種動物模型對中國明對蝦成蝦低溫脅迫下體重的遺傳力進行了估計，在不加共同環(huán)境效應時遺傳力為(0.3271±0.0447)，但加入共同環(huán)境效應后則降為(0.1320±0.0269)，二者均顯著差異于 0(P＜0.01)，與本研究結果基本一致。

徐利永等(2014)對大菱鲆(Scophthalmus maximus)低溫條件下生長性狀進行遺傳參數評估時發(fā)現(xiàn)，其低溫下體重和體長的遺傳力分別為(0.32±0.04)和(0.47±0.06)，為中高等遺傳力，略高于本研究中獲得的遺傳力。此外，本研究結果與 He等(2011)采用全同胞組內相關法估計的 150日齡中國明對蝦體重的遺傳力(0.04～ 0.29)結果大體一致，但體長的遺傳力估計要略低于其結果(0.36～0.51)。而與田燚等(2008)對 145日齡中國明對蝦生長性狀的遺傳力的研究結果基本一致，其體重、體長遺傳力分別為(0.14±0.16)和(0.10±0.06)。黃付友等(2008)采用全同胞組內相關法對“黃海1號”中國明對蝦3月齡和4月齡時的體長遺傳力進行評估，其估計值在0.4～0.5之間，屬于高等遺傳力。Zhang等(2011)對中國明對蝦成蝦的體重遺傳力分別用8種模型進行了分析，結果均為高等遺傳力，估計值在0.44～0.74之間。不同研究中獲得的遺傳力不同，分析主要有以下原因：由于估計遺傳力時所采用的方法不同，分析的家系數量也不同，分析過程中環(huán)境效應和不可預測的遺傳效應會得到不同的遺傳力估計值(He et al, 2011)。不同品種或相同品種的不同生長時期，群體遺傳背景和結構不同，會影響到遺傳力估計值(蔣湘等, 2013)。群體所處生長階段不同，遺傳參數中各方差組分所占比例不同，隨著養(yǎng)殖時間的逐漸增加，顯性效應和母體效應逐漸減少，因此，獲得的遺傳力可能較低(蔣湘等, 2013; 鄭靜靜等,2016)。而本研究中，生長性狀為中低等遺傳力可能是受其生長周期影響，較其他研究中，本研究用蝦日齡較大(平均日齡為252 d)，隨著生長周期增長，母體效應減小，因此，獲得的遺傳力較低。

表2 中國明對蝦成蝦低溫下生長和耐低溫性狀的遺傳參數評估Tab.2 Variance components, heritability for growth traits under low temperature and cold tolerance traits in F. chinensis

表3 中國明對蝦成蝦低溫下生長和耐低溫性狀間的相關分析Tab.3 Correlation analysis among growth traits under low temperature and cold tolerance traits in F. chinensis

3.2 中國明對蝦成蝦耐低溫性狀遺傳力

水產動物的低溫耐受能力通常會選 TAD、CDH和 SS50等指標來描述(Li et al, 2016; 盧其西等,2011)。但本研究基于選育留種的需求，僅選用了SS50指標對中國明對蝦的低溫耐受性進行評估。本研究中獲得的SS50的遺傳力為(0.169±0.078)，經Z檢驗與0存在顯著差異(P＜0.05)，為低等遺傳力。Li等(2016)在中國明對蝦成蝦和凡納濱對蝦成蝦耐低溫性狀的研究中，分別獲得的SS50和CDH的遺傳力也均為低等遺傳力[分別為(0.0192±0.0235)和(0.0258±0.0205)]。在其他水產動物溫度耐受性的研究中，劉寶鎖等(2011)利用2種動物模型對大菱鲆耐高溫性狀進行遺傳力的估計，所獲得耐熱性(Upper thermal tolerance,UTT)的遺傳力分別為(0.026±0.034)和(0.026±0.053)，也均為低等遺傳力。Charo-Karisa等(2005)對于尼羅羅非魚幼魚耐低溫性狀的遺傳力進行了研究，其CDH的遺傳力為(0.08±0.19)，TAD 的遺傳力為(0.09±0.17)。Elderkin等(2004)估計的斑馬貽貝(Driessena polymorpha)面盤幼蟲耐高溫性狀的遺傳力為(–0.125±0.095)。以上研究結果均表明，水產動物中溫度耐受性性狀為低等遺傳力性狀。同時，在水產動物的其他抗逆性的研究中也發(fā)現(xiàn)了抗逆性狀多為低等遺傳力。袁瑞鵬等(2016)對7周和14周凡納濱對蝦高氨氮耐受性的遺傳力進行估計，其高氨氮耐受存活率遺傳力分別為(0.13±0.11)和(0.17±0.08)；張嘉晨等(2016)進行了凡納濱對蝦耐低溶氧性狀的遺傳力估計，也為低等遺傳力(0.07±0.03)。由此可見，水產動物抗逆性狀的遺傳力一般較低。為了給中國明對蝦耐低溫品種的選育工作提供更精確的數據支持，仍需要對中國明對蝦耐低溫性狀的遺傳參數進行進一步的研究。

3.3 中國明對蝦成蝦低溫下生長性狀和耐低溫性狀間的相關分析

在本研究的遺傳參數相關性分析中，體重、體長與 SS50的遺傳相關系數分別為(0.517±0.205)和(0.538±0.203)，表型相關系數分別為(0.187±0.030)和(0.218±0.030)。此研究結果顯示，中國明對蝦成蝦耐低溫性狀與生長性狀間的遺傳相關較強，而表型相關較弱。Charo-Karisa等(2005)對于尼羅羅非魚幼魚耐低溫性狀研究中，體重與耐低溫性狀CDH和TAD的表型相關分別為(0.33±0.17)和(–0.26±0.18)，遺傳相關分別為(0.72±0.81)和(–0.68±0.84)，遺傳相關強于表型相關與本研究結果一致。但Li等(2015、2016)對中國明對蝦成蝦耐低溫性狀研究中，體重與SS50的遺傳相關系數低于表型相關系數；而其對凡納濱對蝦耐低溫的研究中，體重與CDH和SS50的遺傳相關值分別為(–0.7702±0.4583)和(–0.8253±0.4553)，均為負相關，且均與0無顯著差異。

劉寶鎖等(2011)對 40個大菱鲆家系耐熱性的研究表明，體重與 UTT的表型相關和遺傳相關系數分別為 0.04和–1.00，而在 Perry等(2005)對虹鱒魚(Oncorhynchus mykiss)耐熱性的研究表明，體重與UTT之間的遺傳相關幾乎為0 (ra=–0.03±0.18)。在選擇育種中，一些性狀通過直接選育往往得不到理想的效果，由于生物體存在基因連鎖效應和基因多效性，使得不同的性狀間具有一定的相關程度，因此，可以通過選育性狀相關程度較高的性狀進行間接選育(Moore et al, 1999)。本研究表明，中國明對蝦耐低溫性狀SS50與生長性狀間具有較高的遺傳相關系數，因此，中國明對蝦耐低溫新品種的選育可以結合生長性狀選育同時進行。

Castillo-Juárez H, Casares JCQ, Campos-Montes G, et al.Heritability for body weight at harvest size in the Pacific white shrimp, Penaeus (Litopenaeus) vannamei, from a multi-environment experiment using univariate and multivariate animal models. Aquaculture, 2007, 273(1):42–49

Charo-Karisa H, Rezk MA, Bovenhuis H, et al. Heritability of cold tolerance in Nile tilapia, Oreochromis niloticus,juveniles. Aquaculture, 2005, 249(1–4): 115–123

Cnaani A, Gall GAE, Hulata G. Cold tolerance of tilapia species and hybrids. Aquaculture International, 2000, 8(4): 289–298

Elderkin CL, Stoeckel JA, Klerks PL, et al. Heritability of heat tolerance in zebra mussel veligers. Journal of Great Lakes Research, 2004, 30(3): 360–366

He YY, Wang QY, Tan LY, et al. Estimates of heritability and genetic correlations for growth traits in Chinese shrimp Fenneropenaeus chinensis. Agricultural Science and Technology, 2011, 12(4): 613–616

Huang FY, He YY, Li J, et al. Estimates for the heritability of body length in Shrimp Fenneropenaeus chinensis Named“Huanghai No.1”. Periodical of Ocean University of China, 2008, 38(2): 269–274 [黃付友, 何玉英, 李健, 等.“黃海 1號”中國對蝦體長遺傳力的估計. 中國海洋大學學報, 2008, 38(2): 269–274]

Jiang X, Liu JY, Lai ZF. Genetic parameter estimation of growth traits in small abalone, Haliotis diversicolor supertexta.Journal of Fishery Sciences of China, 2013, 20(3): 544–550[蔣湘, 劉建勇, 賴志服. 九孔鮑幼鮑生長性狀的遺傳參數估計. 中國水產科學, 2013, 20(3): 544–550]

Li WJ, Luan S, Luo K, et al. Genetic parameters and genotype by environment interaction for cold tolerance, body weight and survival of the Pacific white shrimp Penaeus vannamei at different temperatures. Aquaculture, 2015, 441: 8–15

Li WJ, Luan S, Luo K, et al. Genetic parameters for cold tolerance and body weight of Chinese fleshy prawn,Fenneropenaeus chinensis. Acta Oceanologica Sinica, 2016,35(6): 12–18

Liu BS, Zhang TS, Kong J, et al. Estimation of genetic parameters for growth and upper thermal tolerance traits in turbot Scophthalmus maximus. Journal of Fisheries of China,2011, 35(11): 1601–1606 [劉寶鎖, 張?zhí)鞎r, 孔杰, 等. 大菱鲆生長和耐高溫性狀的遺傳參數估計. 水產學報, 2011,35(11): 1601–1606]

Lu QX, Lin Y, Bin SY, et al. Analysis of cold tolerance measure at low-temperature in six tilapia families. Journal of Guangxi Normal University (Natural Science Edition), 2011;29(2): 104–109 [盧其西, 林勇, 賓石玉, 等. 羅非魚 6 個家系的低溫耐寒測定分析. 廣西師范大學學報(自然科學版), 2011, 29(2): 104–109]

Luan S, Wang J, Yang G, et al. Genetic parameters of survival for six generations in the giant freshwater prawn Macrobrachium rosenbergii. Aquaculture Research, 2015, 46(6): 1345–1355

Moore SS, Whan V, Davis GP, et al. The development and application of genetic markers for the kuruma prawn Penaeus japonicus. Aquaculture, 1999, 173 (1–4): 19–32

Nguyen NH, Khaw HL, Ponzoni RW, et al. Can sexual dimorphism and body shape be altered in Nile tilapia(Oreochromis niloticus) by genetic means? Aquaculture,2007, 272(S1): S38–S46

Perry GML, Martyniuk CM, Ferguson MM, et al. Genetic parameters for upper thermal tolerance and growth-related traits in rainbow trout (Oncorhynchus mykiss). Aquaculture,2005, 250(1–2): 120–128

Tian Y, Kong J, Luan S, et al. Estimation of genetic parameters for growth traits of Chinese shrimp Fenneropenaeus chinensis. Marine Fisheries Research, 2008, 29(3): 1–6 [田燚, 孔杰, 欒生, 等. 中國對蝦生長性狀遺傳參數的估計.海洋水產研究, 2008, 29(3): 1–6]

Wang J, Kuang YY, Tong GX, et al. Genetic parameters of growth traits in Hucho taimen at different temperature.Journal of Fishery Sciences of China , 2011, 18(1): 75?82[王俊, 匡友誼, 佟廣香, 等. 不同溫度下哲羅鮭幼魚生長性狀的遺傳參數估計. 中國水產科學, 2011, 18(1):75?82]

Wang QY. Principles and practices of breeding in aquatic organisms. Beijing: Science Press, 2013, 42–77 [王清印.水產生物育種理論與實踐. 北京: 科學出版社, 2013,42–77]

Xu LY. Research on growth performance and estimates of heritabilty for growth traits of turbot (Scophthalmus maximus L.) under low temperature. Master′s Thesis of Nanjing Agricultural University, 2014, 38–42 [徐利永. 大菱鲆低溫條件下生長性能研究及生長性狀遺傳力的估計.南京農業(yè)大學碩士研究生學位論文, 2014, 38–42]

Yuan RP. The selective breeding research on the growth,reproduction and high ammonia nitrogen tolerance in Litopenaeus vannamei. Master′s Thesis of Guangdong Ocean University, 2016, 19–33 [袁瑞鵬. 凡納濱對蝦生長、繁殖及高氨氮耐受性的選擇育種研究. 廣東海洋大學碩士研究生學位論文, 2016, 19–33]

Zhang JC, Cao FJ, Liu JY, et al. Estimation on genetic parameters and genetic gain in growth and hypoxic tolerance traits of Litopenaeus vannamei. Oceanologia et Limnologia Sinica, 2016, 47(4): 869–875 [張嘉晨, 曹伏君, 劉建勇,等. 凡納濱對蝦(Litopenaeus vannamei)生長和耐低溶氧性狀的遺傳參數估計和遺傳獲得評估. 海洋與湖沼, 2016,47(4): 869–875]

Zhang TS, Kong J, Luan S, et al. Estimation of genetic parameters and breeding values in shrimp Fenneropenaeus chinensis using the REML/BLUP procedure. Acta Oceanologica Sinica, 2011, 30(1): 78–86

Zheng JJ, Liu JY, Liu JH, et al. Genetic parameterization and breeding value on body weight of Litopenaeus vannamei in different growth stages. Oceanologia et Limnologia Sinica,2016, 47(5): 1005–1012 [鄭靜靜, 劉建勇, 劉加惠, 等. 凡納濱對蝦(Litopenaeus vannamei)不同生長階段體重的遺傳參數和育種值估計. 海洋與湖沼, 2016, 47(5):1005–1012]