Báo cáo sinh học: "Genotype by environment interaction for adult body weights of shrimp Penaeus vannamei when grown at low and high densities" ppt

Original articleGenotype by environment interaction for adult body weights of shrimp Penaeus vannamei when grown at low and high densities Ana M.. Phenotypic plasticity in the growth res

Original article

Genotype by environment interaction for adult body weights of shrimp

Penaeus vannamei when grown at low

and high densities Ana M IBARRA 1*, Thomas R FAMULA2

1

Aquaculture Genetics and Breeding Laboratory, Centro de Investigaciones Biolo´gicas del Noroeste S.C., Mar Bermejo 195, La Paz B.C.S 23090, Mexico

2

Department of Animal Science, University of California Davis, 1 Shields Avenue,

Davis, CA 95616, USA (Received 14 June 2007; accepted 1st April 2008)

Abstract – Shrimp is one of few marine species cultured worldwide for which several selective breeding programs are being conducted One environmental factor that can affect the response to selection in breeding programs is the density at which the shrimp are cultured (low-medium-high) Phenotypic plasticity in the growth response to different densities might be accompanied by a significant genotype by environment interaction, evidenced by a change in heritabilities between environments and by a genetic correlation less than one for a unique trait between environments Our goal was to understand whether different growth densities affect estimates of those genetic parameters for adult body weight (BW) in the Pacific white shrimp (Penaeus vannamei) BW heritabilities were significantly different between environments, with the largest at high density These differences resulted from both an increased additive genetic variance and a decreased environmental variance when grown at high density The genetic correlation between BWs at the two environmental conditions was significantly less than one Whereas these results might be suggestive for carrying out shrimp selective breeding for BW under high density conditions, further understanding of genetic correlations between growth and reproductive traits within a given environment is necessary, as there are indications of reduced reproductive fitness for shrimp grown at high densities.

Litopenaeus / heritability / body weight / density / genotype by environment interaction

1 INTRODUCTION

Production of the Pacific white shrimp, Penaeus (Litopenaeus) vannamei, is presently a worldwide activity This species, whose natural distribution is along the Pacific coast of the western American continent from Mexico to Peru, has been

*

Corresponding author: aibarra@cibnor.mx

 INRA, EDP Sciences, 2008

DOI: 10.1051/gse:2008020

www.gse-journal.org

Article published by EDP Sciences

introduced into the Atlantic coast from USA to Brazil and into different Asian countries [7] Between and within different countries, shrimp is cultured in a large variety of environmental and culture conditions and selective breeding programs are being developed Whereas slight changes in environmental conditions, partic-ularly salinity and temperature, have proven to result in no significant effects on inducing a change in performance rank among genotypes of P vannamei [8,14] and P japonicus [4], larger changes in grow-out temperature (24–30C) and uti-lization of different stocking densities for P japonicus shrimp culture have indi-cated that changes in family ranks can occur for growth traits and survival [5,10] Standard classification of stocking densities for shrimp culture has been rapidly changing, but reported densities range from the minimum stocking densities uti-lized for extensive pond culture (2–5 postlarvaeÆm2) to semi-intensive (8–25 m2), intensive (25–50 m2), and super-intensive (> 50 m2) [11,15] For super-intensive culture in raceways, densities up to 400 m2 have been reported [7] The challenge will come as we consider the possibility of a genotype by environment interaction (GEI) for production traits Commercial pro-duction systems typically rear animals in environments of less than super-intensive density, whereas many breeding programs, facing issues of biosafety, often turn to super-intensive densities The presence of GEI will therefore complicate selection programs Unrecognized changes in the performance rank of families (genotypes) across culture densities resulting from the occurrence of a GEI will lead to a reduced accuracy in breeding value estimation and less than expected genetic gain, although this could depend on the selection scheme being applied [12]

Falconer [6] stated that a unique phenotype ‘‘expressed in two environments may

be considered to be two genetically correlated character states’’ As such, not only

a significant change in heritability values between environments, but also an esti-mate of the genetic correlation less than one (i.e., rg < 1) between a specific pheno-type measured in two environments provide evidence of a significant GEI [17]

We present here an analysis of heritabilities of P vannamei shrimp BW at adult, mature age, after growing 102 replicated genotypes (families) at two grow-out density conditions, and an estimate of the genetic correlation for

BW between grow-out densities

2 MATERIALS AND METHODS

2.1 Parental generation

A two-generation pedigreed structured population of P vannamei shrimp was used for the analyses The base generation of the pedigreed population originated from a newly formed breeding line (produced in 2004, with 102 full-sib families)

was obtained from a commercial shrimp facility in Baja California Sur, Mexico In

2005, a second generation, also with a total of 102 families, was produced and progeny and parents from these 102 families were used for this study

The base generation was kept in replicate for security reasons in two grow-out sites: an earthen pond (origin A) and a plastic-liner-covered pond (origin B) In both of these origins, grow-out density at stocking (1–2 g) was 5–8 shrimpÆm2 Regardless of the same stocking densities, previous observations had indicated that growth in these two environments resulted in different growth rates, possi-bly associated with different natural productivity at each pond type Therefore to test for prior possible fixed effects of parent grow-out environments on growth, dams and sires utilized in the production of the progeny generation used in this study were identified as for their ‘grow-out origin’, their age, and BW obtained after 8 months of grow-out

2.2 Production of the families

All families were produced by artificial insemination using randomly chosen mature individuals Matings took place over a period of 67 days, from February

to April 2005 The only restriction for mating two individuals was if they were known to be related by the pedigree data A total of 102 families were produced, derived from 79 females and 92 males Amongst these 102 families, 8 were related

to one another through a common sire (i.e., as paternal half-sibs), 20 through a common dam, and 34 with both parents in common In addition there were 3 fam-ilies tied by a maternal half-sib trio and 1 family through a paternal half-sib trio Larvae production and juvenile rearing of the 102 families were done as described by Perez-Rostro and Ibarra [14] In short, each family was grown

in an individual 100-L tank from nauplii to postlarvae 15-day old (PL15), using

an initial stocking density of 100 naupliiÆL1, and decreasing it at postlarvae 1-day old (PL1) to 30 PL1ÆL1 Larvae feeding consisted of a microalgae diet (Chaetoceros calcitrans, Thalassiosira fluviatilis, and Tetraselmis suecica) from nauplii V to advanced zoea III, and the addition of Artemia sp from zoea III to PL15 The amount of microalgae and the proportion of Artemia sp added were adjusted daily for each larval stage After this period of approximately 28 days, each family was transferred and grown in 1 m2concrete tanks at an initial stock-ing density of 800 PL15, decreasstock-ing it to 400 when juveniles reached an average size of 0.5 g After 60 days, juveniles reached the size for tagging

2.3 Tagging

A total of 8160 juvenile (1–2 g) shrimp were tagged (80 per family) Families were identified by different color elastomers (NMT, Seattle, WA) injected at

different positions in their abdomen Half of the individuals tagged per family (40), or a total of 4080 individuals, were then transferred to each one of two grow-out densities in Baja California Sur, Mexico

2.4 Test environments

Half of the tagged shrimp (4080 individuals; 1–2 g) were stocked at a low density (5.9 m2) in a 690 m2 plastic-liner covered bottom pond-type tank at CIBNOR; 5–10% of the water volume was exchanged daily Supplemental aer-ation was provided to maintain oxygen levels above 4 mgÆL1 Density was not purposely adjusted during grow-out, but decreased by natural mortality to 1.8 adult shrimpÆm2 at the end of a 9-month grow-out period (June 2005– February 2006) The total number of tagged shrimp recovered was 1241, with

an estimated survival of 30% from 1–2 g to mature age ( 35 g) BWs of all shrimp were collected in February 2006 All families were represented among the survivors, with a mean number per family of 12.2 shrimp (minimum =

1 shrimp per family; maximum = 30 shrimp per family; mode = 11)

The remaining 4080 tagged individuals were stocked at high density (400 m2) at the commercial hatchery APBC, using a 308 m2rectangular race-way tank, and provided with a continuous water flow (15 LÆs1) To achieve this culture density in the raceway, unmarked same age shrimp were simultaneously grown with the tagged shrimp Supplemental aeration was also used, maintain-ing oxygen levels greater than 4 mgÆL1 Density was periodically decreased by natural mortalities, and after 3 months in the raceway, density was adjusted to

50 shrimpÆm2 by randomly transferring out unmarked shrimp At the end of the study after a grow-out period of 10 months (June 2005–March 2006), the density in the raceway was 38 adult shrimpÆm2 The total number of tagged shrimp recovered was 1064, with an estimated survival of 26% from 1–2 g to mature age ( 35 g) BWs of all shrimp were collected in March 2006 All fam-ilies were represented among the survivors, with a mean number per family

of 10.5 shrimp (minimum = 1 shrimp per family; maximum = 24 shrimp per family; mode = 10)

To test for differences in survival patterns of the families when grown at each of the two environments, a correlation analysis was done between the numbers of shrimp surviving per family at each environment The resulting correlation was positive and significant (r = 0.27, P = 0.007) Though signif-icantly different from zero, this correlation suggests a weak (though positive) relationship between survival patterns in the two densities, which is the cause for further concern about the impact of GEI in the construction of breeding programs

2.5 Feeding

Feeding in both cases was done on demand, using a commercial pellet (40% protein) and tray-feeding For the low density tank, food was added to the trays three times a day and for the high density raceway every 2 h

2.6 Statistical and genetic parameter analyses

2.6.1 Model evaluation

Means and standard deviation (SD) of BWs unadjusted for age for each grow-out density and sex were estimated using STATISTICA v.5 software

A linear model was used to evaluate the significance of fixed effects and covariates utilizing SOLAR (Sequential Oligogenic Linkage Analysis Routines; Version 4.0.5 for Linux, downloaded on September 8, 2006 from

covariates, and density, sex, origin, and their interactions were used for fixed effects Neither dam nor sire origin (see Sect 2 for origin A or B) nor their inter-action with density were significant (P > 0.05) and therefore they were excluded from all further analyses

2.6.2 Heritabilities

For the objective of estimating the heritability of BW across the two densities,

a single trait linear animal model was used This model included terms for the fixed effects; sex of shrimp, density class, a sex by density interaction and covar-iates for age and age-squared, and an animal genetic effect Algebraically, the model for the BW of the jth animal (yj) was:

yj ¼ l þ b0xjþ gj þ ej where l is a constant common to all animals, xj is the vector of covariates and fixed effects (i.e., age, sex, density, sex by density) for the jth individual,

b is the vector of corresponding unknown fixed effects, gj is modeled as the additive genetic, polygenic contribution to BW, and ej is the corresponding unknown residual The distribution of the unobservable additive genetic contribution and the residual were assumed to be multivariate normal, and independent of one another In addition, defining g as the vector of additive genetic values for all individuals in the pedigree, we assume that Var g½  ¼ A r2

g, where A is the matrix of numerator relationships and r2

g is the unknown additive genetic variance In a similar fashion, defining e as the vector of residuals, we assume that Var e½  ¼ I r2, where I is an identity

matrix and r2

e is the unknown residual variance It is important to note that, because of the mating system, g may contain effects beyond the additive genetic contribution to BW

2.6.3 Genetic correlations

Estimation of the genetic correlation across environments is a straightforward and commonly used technique for the identification of GEI (e.g., Via and Hawthorne [18]) Expansion of the above single trait model to a multiple trait model across environments is easily accomplished That is, we consider BW

in the low density environment to be a different trait from BW in the high den-sity environment, allowing for the establishment of a multiple trait model The model for BW across the two density classes (yjl or yih), the high density class (subscript ‘h’) and low density class (subscript ‘l’), between two shrimp i, nested

in the high density environment and shrimp j, nested in the low density environ-ment, was:

yjl ¼ ll þ b0lxjl þ gjl þ ejl

yih ¼ lhþ b0hxihþ gihþ eih where ll (lh) is a constant common to all animals in the low (high) density environment, xjl (xih) is the vector of covariates and fixed effects (i.e., age, sex) for the jth individual (ith individual) in the low (high) density environ-ment, bl (bh) is the vector of corresponding unknown fixed effects for the low (high) density environment, gjl(gih) is modeled as the additive genetic, polygenic contribution to BW in the low (high) density environment, and

ejl(eih) is the corresponding unknown residual for the low (high) density envi-ronment Moreover, defining the gl (gh) as vectors of additive genetic values for all individuals in the low (high) density environments,

Var

gl

gh

" #

2

l rlh

rlh r2 h

 A

where r2

l is the additive genetic variance for shrimp in the low density envi-ronment, r2

h is the additive genetic variance for shrimp in the high density environment, rlh is the additive genetic covariance across the low and high density environments, A is the numerator relationship matrix among all shrimp, and  is the symbol for a direct product In addition, define el and

eh as the vectors of residuals for the low and high density BWs, respectively, then

Var el

eh

 

2

el 0

0 r2 eh

 I

where r2

elis the residual variance for shrimp in the low density environment and

r2

ehis the residual variance for shrimp in the high density environment Note that the residual covariance is defined as zero Because individual shrimp are nested within a density class, no animal can be observed in both environmental classes (i.e., traits) Accordingly, the environmental correlation cannot be estimated from the present data Estimates of the unknown variances and covariances (along with the associated heritabilities and genetic correlation) were imple-mented through the public domain computer program MTDFREML [3]

3 RESULTS

Means of BWs by sex and density are presented in TableI

The fixed effects of dam and sire origin, and their interaction with density were not significant (P > 0.05) The covariate age (and age-squared) and the fixed effects of sex and density, as well as their interaction, were all significant (P < 0.001) and therefore included in all further analyses Specifically, the dif-ference between females and males was estimated as 2.91 g (± 0.23) at the high density and 4.39 g (± 0.24) at the low density, using estimates of the sex effects

in the mixed linear model

The heritability of shrimp BW at adult age in the combined environments using the reduced linear model that excluded only origin (i.e., leaving sex and

Table I Means for BW (and SD) of P vannamei male and female adult shrimp when grown at low vs high stocking densities.

Grown at

Sex

body weight (g) (SD)

Coefficient of variation (%)

Range in weights (g) Low density 1241

High density 1064

Coefficient of variation and ranges in weights are included.

age in the model) as a fixed effect was h2= 0.46 ± 0.065 (P < 0.001) TableII

presents estimates of the heritability for adult shrimp BW estimated separately for each density class In each case, the values of heritability were significantly different from zero, as indicated by their standard errors The different heritabil-ities across environments were caused by different additive genetic and environ-mental variances occurring at each environment for growth: a larger environmental variance and smaller additive genetic variance were observed for the low density environment

Our estimate of the genetic correlation between BWs recorded in the two density test environments (Tab.II) was 0.54 ± 0.12, providing strong evidence

of considerable GEI for growth in these two density classes Building crude 95% confidence intervals (i.e., twice the standard error) around the correlation esti-mate, we further get the indication that the genetic correlation is less than one

4 DISCUSSION

Previous studies have found indications of a GEI for growth of P japonicus when just a few families were grown at different stocking densities [5,10] and for growth, biomass, and survival when grown at extreme temperatures [4] Our results confirmed those previous findings when evaluated in P vannamei

It is known that a genetic correlation less than one for a specific trait when mea-sured in two environments is an indication of a significant GEI [18] The observed genetic correlation provided strong evidence that the genes involved

Table II (a) Heritabilities (on diagonal) and genetic correlation (above the diagonal) for BW across density class with standard error Variance (on diagonal) and covariance (above the diagonal) components of (b) genetic, and (c) environmental contribution to growth of P vannamei shrimp across two density classes.

Low density High density (a) Heritability and genetic correlation

(b) Genetic (co)variances

(c) Environmental (co)variances

in shrimp growth at a low density environment are either not the same suite of loci, interact differently, or are differently expressed when shrimp is grown in a high density environment

Further evidence for stating that different genes were involved in growth in the two environments came from the estimated additive genetic and environ-mental variances for each density in the present study, which resulted in different heritability estimates It is well known that a change in environmental conditions can affect the average expression of traits, the amount of genetic and environ-mental variance of the traits, and consequently the heritability of the traits [9] Whereas a change in heritabilities in stressful environments could be a function

of different environmental variances across environments, with the additive genetic variance being constant, that was not the case for the present study Not only was the environmental variance decreased when shrimp were grown

at a high density, but also the additive genetic variance was increased This was in agreement with previous findings in different species for which stressful conditions generally increase genetic variation [9]

Based on these results, it might be thought that selection for shrimp improve-ment should be done after growing it at high densities However, before selec-tion for BW is practiced in this type of environment, a further understanding of the effects on genetic correlations of BW with other traits within this type of environment is necessary In particular, the genetic correlation between BW and reproductive capacity (number of spawns, days to first spawn, fecundity, larval survival, etc.) has to be estimated It is known that genetic correlations among traits can be environment-dependent and can be changed when environ-mental conditions are modified [16] Early life environmental conditions appear

to have a significant effect on the maturation condition of sub-adult P vannamei, with poor conditions resulting in no development of vitellogenic oocytes [13], whereas when grown in normal conditions vitellogenic oocytes are already developing in sub-adults [2] The previous studies indicate that any existing (and beneficial) genetic correlation between BW and reproductive capacity might change at high densities, and that further research is required before pro-posing selective breeding at high densities This conclusion also derives from observations in this study: inasmuch as the present study was not aimed at eval-uating reproductive performance of shrimp grown at each density, data obtained for females selected from the top 10% BW distribution during spawning indicate

a diminished reproductive capacity (measured as egg vitellin concentration, number of spawns, and number of viable nauplii) of females when grown under high density condition when compared with those grown under low density con-dition (A.M Ibarra, unpublished results)

This research was supported by CIBNOR and by the commercial shrimp com-pany Acuacultores de la Peninsula de B.C., S.A de C.V through his General Director, Jaime J Malagamba-Ansotegui We thank Ximena Malagamba and Juan M Mackliz for technical support and data collection

REFERENCES

[1] Almasy L., Blangero J., Multipoint quantitative-trait linkage analysis in general pedigrees, Am J Hum Genet 62 (1998) 1198–1211.

[2] Arcos F.G., Racotta I.S., Palacios E., Ibarra A.M., Ovary development at the onset of gametogenesis is genetically determined and correlated with reproduc-tive traits at maturity in shrimp Litopenaeus (Penaeus) vannamei, Mar Biol 148 (2005) 339–346.

[3] Boldman K.G., Kriese L.A., van Vleck L.D., van Tassell C.P., Kachman S.D.,

A manual for use of MTDFREML A set of programs to obtain estimates of variances and covariances [Draft], US Department of Agriculture, Agricultural Research Service, 1995.

[4] Coman G.J., Crocos P.J., Preston N.P., Fielder D., The effects of temperature on the growth, survival and biomass of different families of juvenile Penaeus japonicus Bate, Aquaculture 214 (2002) 185–199.

[5] Coman G.J., Crocos P.J., Preston N.P., Fielder D., The effects of density on the growth and survival of different families of juvenile Penaeus japonicus Bate, Aquaculture 229 (2004) 215–223.

[6] Falconer D.S., The problem of environment and selection, Am Nat 86 (1952) 293–298.

[7] Funge-Smith S., Briggs M., The introduction of Penaeus vannamei and

P stylirostris into the Asia-Pacific region Case studies, in: Bartley D.M., Bhujel R.C., Funge-Smith S., Olin P.G., Phillips M.J (Eds.), International Mechanisms for the Control and Responsible use of Alien Species in Aquatic Ecosystems, Report of an Ad Hoc Expert Consultation, FAO, Rome, 2005, pp 149–167, http://library.enaca.org/NACA-Publications/alienspecies/a0113e.pdf [consulted: May 2007].

[8] Gitterle T., Rye M., Salte R., Cock J., Johansen H., Lozano C., Sua´rez J.A., Gjerde B., Genetic (co)variation in harvest body weight and survival in Penaeus (Litopenaeus) vannamei under standard commercial conditions, Aquaculture

243 (2005) 83–92.

[9] Hoffmann A.A., Parsons P.A., Evolutionary genetics and environmental stress, 1st edn., Oxford University Press, Oxford, 1991.

[10] Jerry D.R., Preston N.P., Crocos P.J., Keys S., Meadows J.R.S., Li Y., Application of DNA parentage analyses for determining relative growth rates of