Journal of Oceanology and Limnology   2022, Vol. 40 issue(4): 1662-1668     PDF       
http://dx.doi.org/10.1007/s00343-021-0316-6
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

LÜ Ding, ZHANG Chengsong, YU Yang, XIANG Jianhai, LI Fuhua
Heritability estimations of ammonia tolerance and survival of ridgetail white prawn Exopalaemon carinicauda
Journal of Oceanology and Limnology, 40(4): 1662-1668
http://dx.doi.org/10.1007/s00343-021-0316-6

Article History

Received Sep. 10, 2020
accepted in principle Mar. 29, 2021
accepted for publication Aug. 27, 2021
Heritability estimations of ammonia tolerance and survival of ridgetail white prawn Exopalaemon carinicauda
Ding LÜ1,2,3#, Chengsong ZHANG1,2,3#, Yang YU1,2,3, Jianhai XIANG1,2,3, Fuhua LI1,2,3,4     
1 Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2 Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266237, China;
3 Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China;
4 The Innovation of Seed Design, Chinese Academy of Sciences, Wuhan 430072, China
Abstract: Ridgetail white prawn Exopalaemon carinicauda is an important economic marine species in China. To improve the germplasm quality of cultured E. carinicauda, for the first time, we estimated the heritability of ammonia tolerance and survival of E. carinicauda, covering 263 full-sib families in multiple generations. Analysis models including binary linear animal model (BLA), binary linear sire-dam model (BLSD), probit threshold animal model (PTA), and probit threshold sire-dam model (PTSD) were used. The estimated heritability of tolerance to ammonia in different models ranged from 0.09 to 0.22, and that of survival rate ranged from 0.06 to 0.16. Concerning the tolerance to ammonia, both the PTSD and BLSD models showed comparatively better prediction accuracies (0.48 and 0.47, respectively) than PTA and BLA models (0.40 and 0.26, respectively), while to the survival rate, all four models presented relatively low predictive abilities (0.16-0.26). This study provided important information for the future breeding program of E. carinicauda.
Keywords: heritability    Exopalaemon carinicauda    tolerance to ammonia    survival rate    
1 INTRODUCTION

The ridgetail white prawn, Exopalaemon carinicauda (Crustacean: Decapoda: Palaemonidae), is an economically important marine prawn, distributed naturally along China's coast and the west coast of the Korean Peninsula (Li et al., 2003). With fast growth, strong reproduction ability, wide environmental adaptability, and delicious meat, it is one of the most important species in both aquaculture and fishery activities in China. However, farming on this species mainly relies on the capture of post-larvae or parent prawns from coastal areas, and this primitive method depletes wild resources (Bashirullah et al., 1989) and introduces the risk of pathogen dissemination, especially virulent viruses and bacteria (Lightner et al., 1998), which has led to serious economic loss in shrimp industry around the world (Flegel, 2009; Thitamadee et al., 2016). Besides, culturing wild prawns without artificial breeding greatly restricted the production efficiency. To continue developing aquaculture industry of E. carinicauda, genetic improvement on the main economic traits is particularly important.

Ammonia, especially nonionic ammonia, is the most common toxic environmental factor in aquaculture, due to direct discharge of cultivated animals and ammonization of residual feed and excrement (Armstrong et al., 1978; Yu and Hirayama, 1986). It was reported that high concentrations of ammonia can lead to harmful effects on cultured prawns, such as reducing growth (Chen and Kou, 1992; Chen and Lin, 1992a, b), weakening immune system (Liu and Chen, 2004; Wang and Chen, 2005) and increasing mortality (Chen et al., 1990; Lin et al., 1993). Therefore, tolerance of prawns to ammonia stress is an important trait, which enables them to be adaptive to relatively high concentration of ammonia thus to raise the production. Additionally, the output of E. carinicauda depends greatly on the survival rate during cultivation.

Selective breeding is a basic approach for genetic improvement in aquatic species, which offers an opportunity of continuous genetic gain, and can be multiplied and expressed in a huge number of progenies in aquatic species (Gjedrem et al., 2012). In the past decade, numerous shrimp selective breeding programs in Asia and America were initiated and have made considerable achievements in a broad range of traits, including growth, survival, and resistance to a specific pathogen (Hetzel et al., 2000; Argue et al., 2002; Gitterle et al., 2007; Luan et al., 2012, 2014; Sui et al., 2016). However, studies and breeding programs of E. carinicauda are scarce.

To improve the germplasm quality of cultured E. carinicauda, a breeding project aimed at increasing the tolerance to ammonia and survival rate was conducted in 2012–2016. Although complete pedigree of the prawn of five generations and their phenotypic data were available, information on the genetic parameters remained limited. In this study, the heritability of ammonia tolerance and survival of E. carinicauda was estimated in four statistical models. The predictive performance of these models was compared to determine the most optimal fitting model. The results obtained in this study provide valuable information for subsequent selection breeding.

2 MATERIAL AND METHOD 2.1 Material

The breeding program was performed at Nantong Juhao Company, Jiansu Province, China. The base generation (G0), including 57 full-sib families, was bred by collecting females carrying eggs from Laoba port near Hai'an County, Jiangsu Province in 2012. When the offspring were cultured to the stage of post-larvae, 500 individuals from each family were selected randomly and reared in 200-L white plastic tank separately based on a family-culturing unit. To ensure that the rearing conditions were kept as consistently as possible, all the plastic tanks were put in the same workshop. Eighty percent of the water was changed every 2 d, using disinfected seawater from the same reservoir. Water temperature was varied with air temperature (22–32 ℃), at the salinity of 20±2, pH of 8.2–8.5. The light source was natural light. After 50 d, the survival rate of each family was calculated. Meanwhile, another 40 individuals from each family were selected randomly to test their tolerance to ammonia exposure. The testing ammonia concentrations were determined by preliminary experiments to ensure that about half of the prawns died within 96 h after ammonia exposure. The total ammonia nitrogen concentrations used for formal testing were 50 mg/L in G0 and G1; 55 mg/L in G2– G4. After the tolerance to ammonia nitrogen was tested for each family, the top 20 families for each trait were selected to produce G1 by mating among families in 2013. The subsequent generations (G2– G4) were established using similar approaches as G1. A total of 263 full-sib families in 5 generations were used for data analysis. In total, 10 220 prawns were involved in the ammonia tolerance test and 59 250 prawns in the survival test (Table 1). Phenotypes of alive and dead individuals after tolerance to ammonia nitrogen test and survival test were recorded as 1 and 0, respectively.

Table 1 Numbers of full-sib families and prawns tested
2.2 Genetic analysis

In this study, genetic parameters for two traits were analyzed using both linear mixed model and threshold generalized linear mixed model. Additionally, each type of model contained both animal model and sire-dam model. Four models were applied as follows:

Model 1: Binary linear animal model (BLA):

    (1)

where, yij is the survival status of individual i (0 for dead, 1 for alive); μ is the overall mean; yearj is the fixed effect of year j (five levels); c is the regression coefficient for the covariate for day-age of ith offspring bi; ai is the random additive genetic effect of animal i, with a~(0, a2), where A is the additive genetic relationship matrix among all individuals, and it was constructed using pedigree information in the current study; eij is the random error term, with eij ~(0, e2), where I is n×1 identity matrix, n is the number of individuals.

Model 2: Binary linear sire-dam model (BLSD):

    (2)

where sirel+damm are the random additive genetic effect of sire l and dam m, with sirel or damm~(0, sd2)(σsd2 = σsire2 = σdam2), where A is the additive genetic relationship matrix among all individuals; and the other parameters are as described in model 1.

Model 3: Probit threshold animal model (PTA):

    (3)

where λij is the underlying liability of yij, which is assumed to be a cumulative standard normal distribution and the residual variance of λij(σe2) is assumed to be 1; the other parameters are as described in model 1. The probit link function η=Φ-1(pi) is used, with inverse link , where Φ is the cumulative normal density function, and pi is the probability of alive for individual i. The Bernoulli distribution for a binary trait for an individual prawn with yi=1 (alive) and yi=0 (dead) is the probability (yi|pi)=(pi)yi(1–pi)i1–y.

Model 4: Probit threshold sire-dam model (PTSD):

    (4)

Model 4 is a probit threshold sire-dam model. The parameters in model 4 are well described in model 2 and 3 and the probit link function was also applied. In particular, σe2=1 was also used in model 4.

When using animal models (models 1 and 3), the heritability of both traits was estimated from a single trait model. Phenotypic and genetic correlations were obtained using bivariate analyses, using the same statistical model. Phenotypic variance was the sum of additive genetic variance and residual variance as: . Heritability was calculated as the ratio between additive genetic variance and phenotypic variance: . Under sire-dam model (models 2 and 4), Heritability was calculated as: . Besides, the heritability would be overestimated in probit threshold model, so the estimates in models 3 and 4 were adjusted according to Robertson and Lerner (1949).

Genetic and phenotypic correlations between the two traits were calculated based on animal models (models 1 and 3) as follows: where C is the genetic or phenotypic covariance, σX2 and σY2 are the additive genetic or phenotypic variance of trait X and trait Y, respectively. The Z-score, computed as the ratio between estimate of each parameter and its standard error (SE), was used for testing whether each parameter estimate was significantly different from zero.

The model solutions and genetic parameter estimations were processed in ASReml 4.0 (Gilmour et al., 2016).

2.3 Model comparison

To compare predictive abilities of four models in this study, we used phenotype dataset containing the first four generations (2012–2015) to calculate family estimated breeding value (EBVs) in 2016 from each model. Then, Spearman correlation coefficients between family EBVs and actual family average phenotype (survival rate) in 2016 were calculated. Since we paid more attention to the ranking difference of two sets of data, the Spearman correlation coefficient was adopted. The EBVs for families were estimated using ASReml and Spearman correlation coefficients were calculated using R (R Core Team, 2018).

3 RESULT 3.1 Descriptive statistics for phenotype data

The coefficient of variation of the survival rate in G0 was 27.21%, being higher than those of subsequent generations, ranging from 5.40% to 10.67%, and the data of cultured survival rate in 2013 were missing (Fig. 1).

Fig.1 Distribution of family survival rate after being raised for 50 d under normal conditions (a) and after 96-h acute ammonia stress (b) in each generation The black dots are outliers.

The coefficients of variation in family survival rate after 96-h acute ammonia stress test are relatively larger, ranging from 20.50% to 59.08% (Fig. 1). And the overall survival rate after tolerance to ammonia test in 2012–2016 was 47.72%, 49.35%, 50.76%, 48.71%, and 44.22%, respectively.

3.2 Genetic analysis

Variance components, heritability for ammonia tolerance and survival calculated using different models are shown in Table 2. The heritability estimates of tolerance to ammonia ranged from 0.09±0.01 to 0.22±0.03, and those of survival ranged from 0.06±0.01 to 0.16±0.02. Heritability for tolerance to ammonia was overall larger than that for survival. For both traits, estimates from sire-dam models (models 2 and 4) were close and they were higher than those from the other models. Estimated heritability using different models was highly significantly different from zero (P < 0.01) through Z-score.

Table 2 Variance components, heritability of tolerance to ammonia and survival in four statistical models

Genetic and phenotypic correlations between ammonia tolerance and survival calculated using animal models are shown in Table 3. Genetic correlations were 0.20±0.10 and 0.20±0.04 based on model 1 and model 3, respectively. Phenotypic correlations were 0.29±0.13 and 0.09±0.02, respectively. All of them were low but significantly different from zero (P < 0.05).

Table 3 Genetic correlation (γg) and phenotypic correlation (γp) between tolerance to ammonia and survival in models 1 and 3

The family EBVs of two traits were also obtained from different models and their relationships are presented in Table 4. The Pearson correlation coefficients between EBVs calculated using different models were all greater than 0.9.

Table 4 Pearson correlation coefficients among EBVs from different models
3.3 Comparisons among different models

Spearman correlation coefficients between family EBVs based on data of 2012–2015 and true family average phenotype (survival rate) in 2016 are presented in Table 5. For tolerance to ammonia, Spearman correlation coefficients under sire-dam models (models 2 and 4) were larger than those obtained from animal models (models 1 and 3), which reached a moderate correlation level. However, for survival, all models showed a low predictive ability (Table 5).

Table 5 Spearman correlation coefficients between family EBVs within G4 based on dataset consisting of G0–G3 and actual family average phenotype (survival rate) within G4 in different models
4 DISCUSSION 4.1 Genetic parameter

Understanding the additive genetic variance existing in the target population was the foundation of initiating a selection-breeding program. One of the main objectives of this study was to estimate the heritability of tolerance to ammonia and survival in E. carinicauda. The results showed that heritability for tolerance to ammonia or survival based on four models was all highly significantly different from zero (P < 0.01), which indicated that improvement on these two traits in E. carinicauda was theoretically feasible. Furthermore, heritability for tolerance to ammonia estimated in this study belonged to the medium, based on the following categorization: low (0.05–0.15), medium (0.20–0.40), high (0.45–0.60), and very high (> 0.65) (Cardellino and Rovira, 1987), and that for survival was lower than tolerance to ammonia. Until present, knowledge about the heritability of ammonia tolerance or survival in crustaceans is still limited. The heritability for tolerance to ammonia in Litopenaeus vannamei juveniles was estimated to be 0.148–0.264 (Li et al., 2016; Lu et al., 2017), and it was reported as 0.13±0.06 in Marsupenaeus japonicus (Jiang et al., 2017). For the survival trait, the heritability was estimated to be 0.04–0.10 in L. vannamei (Gitterle et al., 2005, 2007) and 0.02–0.14 in Macrobrachium rosenbergii (Luan et al., 2014; Vu et al., 2017). Estimates in the present and previous studies were similar for both traits. What called for special attention was that we did our best to maintain the environmental conditions consistent among the families in the trait-testing process, but subtle differences were inevitable. Besides, the absence of common environmental effects may lead to the bias of heritability estimates (Gallardo et al., 2010).

Analysis of the genetic correlation and phenotypic correlation between ammonia tolerance and survival showed a low but significant positive genetic correlation (0.20±0.10 and 0.20±0.04, respectively). This result suggested that selection on one trait would not cause a negative response to another one.

4.2 Model comparison

In previous quantitative genetics studies, the predictive abilities of statistic models were evaluated through model cross-validation (Gitterle et al., 2006; Saatchi et al., 2011; Rezende et al., 2012; Lever et al., 2016). In this study, we used a correlation between the actual family average phenotype in G4 and their family EBVs calculated from each model as the criteria for model predictive abilities to make the statistical model closer to practical application. Although the family EBVs calculated from different statistical models showed high correlations (Table 4), the predictive abilities were varied considerably among models (Table 5). For the tolerance to ammonia, the predictive abilities of threshold sire-dam model (model 4) were slightly higher than the linear sire-dam model (model 2), and clearly higher than animal models (models 1 and 3), suggesting that sire-dam models (both linear and threshold) are the optimal models. Previous studies indicated that animal threshold model was not suitable, because it tended to give a biased estimation of genetic parameters (Hoeschele and Tier, 1995; Luo et al., 2001; Stock et al., 2007). Therefore, threshold sire-dam model should be a preferred model for genetic evaluation on threshold traits in future selection breeding.

The predictive abilities of different models for survival traits were all unsatisfactory. There was no or weak correlation between predicted family EBVs and actual phenotype values because survival is complicated and influenced by many factors. During aquaculture, the survival trait contains various traits, such as resistance to different pathogens, tolerance to different physical and chemical stresses and other lethal factors. Considering that the low predictive accuracies were obtained in survival trait, it is not an ideal breeding target. In fact, selective breeding aimed at increasing resistance to specific diseases, which is measured as survival under challenge testing with a specific pathogen, seemed to be a more targeted approach to improve the survival rate of aquaculture (Ødegård et al., 2011). Besides, as the accumulation of data and pedigree depth, selection based on EBVs instead of phenotypes is becoming increasingly urgent and effective to accelerate the progress of selective breeding.

5 CONCLUSION

Significant additive genetic variations for tolerance to ammonia and survival existed in E. carinicauda. Sire-dam models presented relatively high prediction accuracies for tolerance to ammonia of the prawn, which should be a preferred statistical model for genetic evaluation.

6 DATA AVAILABILITY STATEMENT

All data generated and/or analyzed during this study are included in the manuscript.

References
Argue B J, Arce S M, Lotz J M, Moss S M. 2002. Selective breeding of Pacific white shrimp (Litopenaeus vannamei) for growth and resistance to Taura syndrome virus. Aquaculture, 204(3-4): 447-460. DOI:10.1016/S0044-8486(01)00830-4
Armstrong D A, Chippendale D, Knight A W, Colt J E. 1978. Interaction of ionized and un-ionized ammonia on short-term survival and growth of prawn larvae, Macrobrachium rosenbergh. The Biological Bulletin, 154(1): 15-31. DOI:10.2307/1540771
Bashirullah A K M, Mahmood N, Matin A K M A. 1989. Aquaculture and coastal zone management in Bangladesh. Coastal Management, 17(2): 119-127. DOI:10.1080/08920758909362080
Cardellino R, Rovira J. 1987. Mejoramiento Genético Animal. Hemisferio Sur, Buenos Aires. 253p.
Chen J C, Kou Y Z. 1992. Effects of ammonia on growth and molting of Penaeus japonicus juveniles. Aquaculture, 104(3-4): 249-260. DOI:10.1016/0044-8486(92)90207-2
Chen J C, Lin C Y. 1992a. Effect of ammonia on growth of Penaeus penicillatus juveniles. Comparative Biochemistry and Physiology Part C: Comparative Pharmacology, 101(3): 443-447. DOI:10.1016/0742-8413(92)90067-H
Chen J C, Lin C Y. 1992b. Effects of ammonia on growth and molting of Penaeus monodon juveniles. Comparative Biochemistry and Physiology Part C: Comparative Pharmacology, 101(3): 449-452. DOI:10.1016/0742-8413(92)90068-I
Chen J C, Ting Y Y, Lin J N, Lin M N. 1990. Lethal effects of ammonia and nitrite on Penaeus chinensis juveniles. Marine Biology, 107(3): 427-431. DOI:10.1007/BF01313424
Flegel T W. 2009. Current status of viral diseases in Asian shrimp aquaculture. Israeli Journal of Aquaculture-Bamidgeh, 61(3): 229-239.
Gallardo J A, Lhorente J P, Neira R. 2010. The consequences of including non-additive effects on the genetic evaluation of harvest body weight in Coho salmon (Oncorhynchus kisutch). Genetics Selection Evolution, 42(1): 19. DOI:10.1186/1297-9686-42-19
Gilmour A R, Gogel B J, Cullis B R, Thompson R. 2016. ASReml User Guide Release 4.0. VSN International Ltd., Hemel Hempstead.
Gitterle T, Johansen H, Erazo C, Lozano C, Cock J, Salazar M, Rye M. 2007. Response to multi-trait selection for harvest weight, overall survival, and resistance to white spot syndrome virus (WSSV) in Penaeus (Litopenaeus) vannamei. Aquaculture, 272(S1): S262.
Gitterle T, Ødegård J, Gjerde B, Rye M, Salte R. 2006. Genetic parameters and accuracy of selection for resistance to White Spot Syndrome Virus (WSSV) in Penaeus (Litopenaeus) vannamei using different statistical models. Aquaculture, 251(2-4): 210-218. DOI:10.1016/j.aquaculture.2005.06.008
Gitterle T, Rye M, Salte R, Cock J, Johansen H, Lozano C, Suárez J A, Gjerde B. 2005. Genetic (co)variation in harvest body weight and survival in Penaeus (Litopenaeus) vannamei under standard commercial conditions. Aquaculture, 243(1-4): 83-92. DOI:10.1016/j.aquaculture.2004.10.015
Gjedrem T, Robinson N, Rye M. 2012. The importance of selective breeding in aquaculture to meet future demands for animal protein: a review. Aquaculture, 350-353: 117-129. DOI:10.1016/j.aquaculture.2012.04.008
Hetzel D J S, Crocos P J, Davis G P, Moore S S, Preston N C. 2000. Response to selection and heritability for growth in the Kuruma prawn, Penaeus japonicus. Aquaculture, 181(3-4): 215-223. DOI:10.1016/S0044-8486(99)00237-9
Hoeschele I, Tier B. 1995. Estimation of variance components of threshold characters by marginal posterior modes and means via Gibbs sampling. Genetics Selection Evolution, 27(6): 519-540. DOI:10.1186/1297-9686-27-6-519
Jiang X, Zheng J J, Xie M, Liu Y K, Zeng F X, Liu J Y. 2017. Estimation of genetic parameters for high ammonia nitrogen tolerance and growth in kuruma prawn Marsupenaeus japonicus. Fisheries Science, 36(6): 700-706. (in Chinese with English abstract)
Lever J, Krzywinski M, Altman N. 2016. Model selection and overfitting. Nature Methods, 13(9): 703-704. DOI:10.1038/nmeth.3968
Li W J, Lu X, Luan S, Luo K, Sui J, Kong J. 2016. Heritability of body weight and resistance to ammonia in the Pacific white shrimp Litopenaeus vannamei juveniles. Chinese Journal of Oceanology and Limnology, 34(5): 1025-1033. DOI:10.1007/s00343-016-5034-0
Li X Z, Liu R Y, Liang X Q. 2003. The zoogeography of Chinese Palaemonoidea fauna. Biodiversity Science, 11(5): 393-406. (in Chinese with English abstract) DOI:10.17520/biods.2003047
Lightner D V, Hasson K W, White B L, Redman R M. 1998. Experimental infection of western hemisphere penaeid shrimp with Asian white spot syndrome virus and Asian yellow head virus. Journal of Aquatic Animal Health, 10(3): 271-281. DOI:10.1577/1548-8667(1998)010<0271:EIOWHP>2.0.CO;2
Lin H P, Thuet P, Trilles J P, Mounet-Guillaume R, Charmantier G. 1993. Effects of ammonia on survival and osmoregulation of various development stages of the shrimp Penaeus japonicus. Marine Biology, 117(4): 591-598. DOI:10.1007/BF00349770
Liu C H, Chen J C. 2004. Effect of ammonia on the immune response of white shrimp Litopenaeus vannamei and its susceptibility to Vibrio alginolyticus. Fish & Shellfish Immunology, 16(3): 321-334.
Lu X, Luan S, Cao B X, Meng X H, Sui J, Dai P, Luo K, Shi X L, Hao D C, Han G M, Kong J. 2017. Estimation of genetic parameters and genotype-by-environment interactions related to acute ammonia stress in Pacific white shrimp (Litopenaeus vannamei) juveniles at two different salinity levels. PLoS One, 12(3): e0173835. DOI:10.1371/journal.pone.0173835
Luan S, Yang G L, Wang J Y, Luo K, Chen X F, Gao Q, Hu H L, Kong J. 2014. Selection responses in survival of Macrobrachium rosenbergii after performing five generations of multi-trait selection for growth and survival. Aquaculture International, 22(3): 993-1007. DOI:10.1007/s10499-013-9722-x
Luan S, Yang G L, Wang J Y, Luo K, Zhang Y F, Gao Q, Hu H L, Kong J. 2012. Genetic parameters and response to selection for harvest body weight of the giant freshwater prawn Macrobrachium rosenbergii. Aquaculture, 362-363: 88-96. DOI:10.1016/j.aquaculture.2012.05.011
Luo M F, Boettcher P J, Schaeffer L R, Dekkers J C M. 2001. Bayesian inference for categorical traits with an application to variance component estimation. Journal of Dairy Science, 84(3): 694-704. DOI:10.3168/jds.S0022-0302(01)74524-9
Ødegård J, Gitterle T, Madsen P, Meuwissen T H E, Yazdi M H, Gjerde B, Pulgarin C, Rye M. 2011. Quantitative genetics of Taura syndrome resistance in pacific white shrimp (Penaeus vannamei): a cure model approach. Genetics Selection Evolution, 43(1): 14. DOI:10.1186/1297-9686-43-14
R Core Team. 2018. R: a Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna.
Rezende F M, Ferraz J B S, Eler J P, Silva R C G, Mattos E C, Ibáñez-Escriche N. 2012. Study of using marker assisted selection on a beef cattle breeding program by model comparison. Livestock Science, 147(1-3): 40-48. DOI:10.1016/j.livsci.2012.03.017
Robertson A, Lerner I M. 1949. The heritability of all-or-none traits: viability of poultry. Genetics, 34(4): 395-411. DOI:10.1093/genetics/34.4.395
Saatchi M, McClure M C, McKay S D, Rolf M M, Kim J, Decker J E, Taxis T M, Chapple R H, Ramey H R, Northcutt S L, Bauck S, Woodward B, Dekkers J C, Fernando R L, Schnabel R D, Garrick D J, Taylor J F. 2011. Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation. Genetics Selection Evolution, 43(1): 40. DOI:10.1186/1297-9686-43-40
Stock K F, Distl O, Hoeschele I. 2007. Influence of priors in Bayesian estimation of genetic parameters for multivariate threshold models using Gibbs sampling. Genetics Selection Evolution, 39(2): 123-137. DOI:10.1186/1297-9686-39-2-123
Sui J, Luan S, Luo K, Meng X H, Cao B X, Liu N, Li W J, Lu X, Chai Z, Kong J. 2016. Genetic parameters and response to selection of harvest body weight of the Chinese shrimp Fenneropenaeus chinensis after five generations of multi-trait selection. Aquaculture, 452: 134-141.
Thitamadee S, Prachumwat A, Srisala J, Jaroenlak P, Salachan P V, Sritunyalucksana K, Flegel T W, Itsathitphaisarn O. 2016. Review of current disease threats for cultivated penaeid shrimp in Asia. Aquaculture, 452: 69-87.
Vu N T, Trong T Q, Nguyen N H. 2017. Effects of selection for fast growth on survival rate during grow-out phase in giant freshwater prawn (Macrobrachium rosenbergii). BMC Genetics, 18(1): 56.
Wang L U, Chen J C. 2005. The immune response of white shrimp Litopenaeus vannamei and its susceptibility to Vibrio alginolyticus at different salinity levels. Fish & Shellfish Immunology, 18(4): 269-278.
Yu J P, Hirayama K. 1986. The effect of un-ionized ammonia on the population growth of the rotifer in mass culture. Nippon Suisan Gakkaishi, 52(9): 1509-1513.