Institute of Oceanology, Chinese Academy of Sciences
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
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
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 MaterialThe 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.
2.2 Genetic analysisIn 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):
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, Aσ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, Iσe2), where I is n×1 identity matrix, n is the number of individuals.
Model 2: Binary linear sire-dam model (BLSD):
where sirel+damm are the random additive genetic effect of sire l and dam m, with sirel or damm~(0, Aσ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):
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
Model 4: Probit threshold sire-dam model (PTSD):
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:
Genetic and phenotypic correlations between the two traits were calculated based on animal models (models 1 and 3) as follows:
The model solutions and genetic parameter estimations were processed in ASReml 4.0 (Gilmour et al., 2016).
2.3 Model comparisonTo 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 dataThe 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).
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 analysisVariance 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.
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).
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.
3.3 Comparisons among different modelsSpearman 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).
4 DISCUSSION 4.1 Genetic parameterUnderstanding 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 comparisonIn 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 CONCLUSIONSignificant 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 STATEMENTAll data generated and/or analyzed during this study are included in the manuscript.
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.
|