Genetic Characterization of Six Stocks of <i>Litopenaeus vannamei</i> Used in Cuba for Aquaculture by Means of Microsatellite Loci

Pérez-Beloborodova, Anna; Artiles-Valor, Adriana; Pérez-Jar, Lourdes; Hernández-Martínez, Damir; Guerra-Aznay, Missael; Espinosa-López, Georgina

doi:https://doi.org/10.1155/2012/352165

International Journal of Zoology

On this page

Abstract Introduction Materials and Methods Results Discussion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2012 | Article ID 352165 | https://doi.org/10.1155/2012/352165

Genetic Characterization of Six Stocks of Litopenaeus vannamei Used in Cuba for Aquaculture by Means of Microsatellite Loci

Anna Pérez-Beloborodova,¹Adriana Artiles-Valor,²Lourdes Pérez-Jar,²Damir Hernández-Martínez,¹Missael Guerra-Aznay,²and Georgina Espinosa-López³

Academic Editor: Pung-Pung Hwang

Received30 Aug 2011

Revised01 Nov 2011

Accepted11 Nov 2011

Published22 Feb 2012

Abstract

Four microsatellite loci were used to achieve genetic characterization of six stocks from Litopenaeus vannamei used for aquaculture in Cuba: second generation from first introduction (S2-1), first generation from the second one (S1-2), from the third one (S1-3), and the fourth one (S1-4) and the crossings from two parental population: first generation from the first with first generation from the third (S1-1 × S1-3) and first generation from the second with first generation from the third (S1-2 × S1-3). 66% (16/24) of genetic systems in total loci were in genetic disequilibrium. The four microsatellite loci were polymorphic for all six stocks. Major quantities of allelic variants correspond to locus Pvan 1758, which is at the same time that one where there are private alleles from first generation of the third. All Fst comparisons were significant. This indicates big differences between stocks. The highest values are those in which there is presence of the second introduction. This introduction and its descendants are also more consanguineous.

1. Introduction

In Latin America shrimp-producer countries, the Pacific white shrimp, Litopenaeus vannamei, is the most representative species, with about 90% of production. Native from East Pacific, and from the tropical American continent, this species has shown an excellent culture adaptation and has been more resistant to salinity, oxygen, and temperature fluctuations. That is why, in the last twenty years, Litopenaeus vannamei, has been introduced in many culture programs, and nowadays it is the second culture after Penaeus monodon.

In 2003 the first introduction of two stocks of White Pacific Shrimp, Litopenaeus vannamei, [1] was achieved in Cuba, imported from USA, Shrimp Improving System and so handling, nutrition, and health techniques are well established, as well as the assessment of genetic variation in farms that had before cultured the indigenous species Litopenaeus schmitti. In total, five stocks have been introduced, and all of them have been characterized using microsatellite techniques [2–4].

Genetic studies have a capital importance in shrimp industry, in order to determine genetic variability level either in natural or cultured populations, but mainly to know when the latter could be enriched with new specimens [5]. It is moreover important to know the structure of natural population from which those specimens will be taken [6] and also to have good markers that allow population and family studies. It is clear that the priority should be given to the domestication and handling of broodstocks through the application of genetic techniques [7–9]. The aim of the present work is to characterize different cultured stocks of Litopenaeus vannamei, used in Cuba for aquaculture, by means of microsatellite markers.

2. Materials and Methods

2.1. Samples

Samples were taken from pleopods of 30 shrimps, specifically from the fourth pair, between exo-and endopodite. The same male and female quantity was taken randomly in shrimps from the Postlarvae Production Hatchery “YAGUACAM,” Cienfuegos, Cuba.

Individuals from first or second generation of introduced shrimps were taken for this genetic characterization as described below. Two first parental generations from first and second introduced stocks were randomly crossed with descendants of the first generation of the third introduced stock resulting in (S1-2 S1-3) and (S1-1 × S1-3) as it is called in this work.

Thus, characterized lots were second generation from first introduced stock (S2-1), first generation from the second one (S1-2), from the third one (S1-3), and from the fourth one (S1-4) and the crossings: first generation from the first with first generation from the third (S1-1 × S1-3) and first generation from the second with first generation from the third (S1-2 × S1-3).

2.2. Microsatellite Genotyping

DNA isolation, amplification programs, and electrophoresis procedures were carried out as described in [2]. Used loci were M-1, isolated from L. vannamei by [10] and Pvan 0040, Pvan 1758, and Pvan 1815 obtained from the same species by [11]. As weight allelic controls, samples previously genotyped by [2], as well as PGEM, were used.

2.3. Statistics

Allele number by locus, allele frequency, and observed and expected heterozygosity for each locus, as well as, the Hardy-Weinberg equilibrium, were determined by GeneAlEx 6.1 program [12]. A locus was considered as polymorphic if it presented at least two alleles and when the most common frequency of an allele did not exceed 95%.

The FSTAT statistical package version 2.93 [13] was used to calculate the linkage disequilibrium, and also the inbreeding coefficient within populations, values [14], and pairwise and values between populations, after 1500 permutations. Significance levels were assessed through the Markov chains using 10000 Dememorisation, 100 batches, and 5000 iterations per batch using the Genepop Version 4.1.0 [15].

Assignment tests were used to clarify the belonging of individuals to the stocks. The Bayesian assignment test was performed [16] among shrimp stocks from Cuba using GeneClass 2 [17]. Simulations were run to determine the probability of assignment, using a probability of rejection set at .

Relatedness coefficient () according to [18] between individuals’ pairs was calculated by GeneAlEx 6.1 program [12]. The equation for this coefficient taking into consideration codominant marker is where represents individuals, all loci, allelic positions (2 for diploids, 1 for haploids), the frequency of the individual for the locus and allelic position l, the frequency of the allele in the group or individual with it is compared , and the total frequency in the population.

Relatedness coefficient must be to unrelated individuals, to half sibs, and to full sibs [18].

Relatedness means coefficients were compared with the Kruskal Wallis test (nonparametric ANOVA) by means of graphPad InStat version 3.00 [19]

3. Results

3.1. General Parameters of Genetic Variation of Six Stocks Descendants from Original Introductions: Hardy-Weinberg Deviation of Equilibrium

As it has been previously reported, L. vannamei shrimps introduced from SIS are characterized [2–4].

In the present work, the genetic variation analysis was performed to six relative stocks of those introductions by means of four microsatellite loci: M1, isolated from L. vannamei, and Pvan 0040, Pvan 1758, and Pvan 1815 also obtained from these species. Linkage disequilibrium was analyzed, and as it is known [2, 4, 11, 20], those loci for L. vannamei are not linked.

In Table 1 are shown the main genetic estimation parameters, such as number of alleles (Na) and number of effective alleles (Ne), observed and expected heterozigocities (Ho and He), and so, the Hardy-Weinberg equilibrium deviations (F_IS values). As several stocks are compared, private alleles are also provided (Np). Note that first generation from second introduction (S1-2) is the one of less genetic variability, within first generation stocks, and even lesser than the only second generation stock analyzed. Highest variation stock belongs to crossing of first generation from the first with first generation from the third (S1-1 × S1-3). The lowest relative variation of any generations and stocks which included second introduction is remarkable and stocks which included second introduction.

Table 1

Genetic variability parameters and the Hardy-Weinberg deviation of equilibrium () and its probability (Pfis) of four microsatellites loci isolated from L. vannamei, M-1 [10], Pvan 0040, Pvan 1758 and Pvan 1815 [11], in six different cultured stocks of L. vannamei used for aquaculture in Cuba. n: number of samples. Na: number of alleles, Ne: number of effective alleles, Np: number of private alleles, He: expected heterozygosity, Ho: observed heterozygosity. The associated probability was estimated using the Markov chains (10000 Dememorisation, 100 batches, and 5000 iterations per batch using Genepop program Version 4.1.0. Numbers in parenthesis are standard mean error.

Normally, the Hardy-Weinberg equilibrium conditions are not accomplished in cultures. In Table 1 it is observed that 66% (16/24) of genetic systems in total loci are significant in genetic disequilibrium by homozygotes excess.

3.2. Differentiation between Stocks

In Figure 1 frequencies of alleles of each microsatellite for all stocks are shown. Allele sizes are on range reported by other authors [2–4, 10, 20–22]. All loci were polymorphic and those with more allelic variants were Pvan 1758 (12 variants) and Pvan 1815 (9 variants), which is coincident with previous studies. Private alleles found in this paper corresponded to those loci. In the case of M1, distribution of frequencies is bimodal, with two predominant sizes, which includes the majority of stocks: 216, which contains all of them, and 202 in which it only lacks S1-4. In total there are eight allelic variants, and there are not private alleles for this microsatellite region.

(a)

(b)

(c)

(d)

Pvan 0040 locus is that one which has less allelic variants, with only five. Its frequency distribution is unimodal and in 141 size all stocks are represented. Two other loci, Pvan 1758 and Pvan 1815, are unimodal and have private alleles that contribute to the highest quantity of allelic variants. Private alleles correspond to the first generations of the third and fourth introductions. Sizes and frequencies of appeared alleles (within parenthesis) are, for locus 1758: 182 (0.025), 184 (0.250), and 193 (0.025), which appeared in first generation of the third and 210 (0.023) from first generation of the 4th. To 1815 there is only a private allele from first generation of the 4th: 111 (0.091).

The assignment test results are represented in Table 2. It revealed that 93.3% (125/134) from the shrimp were correctly assigned. S1-3 stock was the most failed assignation with only 80% of right assignation.

calculations are shown on Table 3, comparing all stocks among them. It is observed that all comparisons are statistically significant. The biggest values are those that have the presence of the second introduction. On the other hand, less values are observed in comparisons between S1-3, S1-4, and S2-1.

3.3. Relatedness Values

Relatedness coefficients distributions between individual pairs of six stocks are shown in Figure 2, calculated according to [18]. All stocks, except S1-2 and crossing that contains it (S1-2 × S1-3), present unimodal distributions. Three of them present the mode in negative relatedness values (S1-3, S1-4, S1-1 × S1-3), and the only one from a second generation S2–1 presents the mode in zero. In the case of exceptions, for S1-2, one of the modes is near zero, while the second one is near positive values but lower than 0.5. For crossing (S1-2 × S1-3), the first mode is in negative values and the other one near zero. Mean relatedness coefficients values (Figure 2, right inside) have not significant differences and have standard deviation values. Those are for each stock, S2-1: −0.028; S1-2: 0.044; S1-3: −0.063; S1-4: −0.023; S1-2 × S1-3: −0.019; S1-1 × S1-3: −0.044.

4. Discussion

4.1. Genetic Diversity within L. vannamei Stocks

The four microsatellite loci were polymorphic for all six stocks. The majority of expected and observed heterozigocities are on range according to many authors using microsatellites for penaeid shrimps [2–4, 20–26]. However, they are low and above range for one of the markers, Pvan 0040, which is the less variable loci of all stocks in which the second introduction is involved: in the first generation of the second introduction Ho = 0.24 and He = 0.21 and in crossing of first generation of second introduction with first generation of the first. Previously, [2–4] reported a great contribution of this locus to the decrease of total heterozigocities. Pvan 0040 was monomorphic for the original fourth and fifth introductions [3, 4], and so it could be considered that this locus is sufficiently sensitive to changes in different stocks.

Equilibrium deviations are significant in most loci by stocks, which in accord with other authors’ results [3, 4, 20, 22].

Most of allelic frequency graphics are unimodals, with exception of M1 locus. The biggest quantities of allelic variants correspond to locus Pvan 1758, which is at the same time that one in which there are private alleles from first generation of the third. This is in agreement with genetic variability reported for that locus in previous works [20, 22, 25].

4.2. Differentiation between Stocks

Indeed all comparisons were significant. They indicate big differences between stocks. The biggest values are those in which there is presence of the second introduction. In a similar way, other authors also obtained significant differences in natural populations as well as in cultured stocks [20, 22]. Using allozymes, the authors in [27, 28] also obtained a genetic differentiation in wild population of L. shmitti in Cuba and in wild and cultured L. vannamei stocks in Norwest Mexico, respectively. Those differences are in agreement with the assigning exercise, in which only first generation of the third stock had low percentages of assigning. It could be because there are components of those stocks in others obtained by meanings of crossing with this third introduction.

Although studies with L. vannamei in Cuba indicate a decrease in genetic variability, existing stocks at the moment of this genetic characterization, even with little variations, seem, however, significantly different between them. The reason could be that in their origins they constitute genetic lines well diverse in their origins in respect to the main objective for production (e.g., illness resistance, increasing growth, etc.). However, once the selection pressure is increased, genetic variation is decreased as it could be seen in different genetic variability reports for several introductions [4, 29].

4.3. Relatedness Coefficients

Borrell et al. [2] pointed out the risk that could implicate using individuals from the second introduced stock for crossing with others without any methodology such as the following, of relatedness coefficient, that offers information about inbreeding levels. For the second stock, they obtained a mean relatedness coefficient indicative of more related individuals than those of the first one, similar to the results of this work, in which the most consanguineous stock is descendants of this second introduced stock. The highest value, indicative of a great relatedness, was obtained for the fifth introduced stock of L. vannamei in Cuba [4], and that is why the use of this molecular tool is of crucial importance for designing and following crossings. The rest of the stocks show means and distributions in agreement with individuals presumed not related, similar to that obtained by [30] for turbot, near zero. The second introduction and its descendants, as well as crossings that contain it, the same as the fifth [4], should be followed with this methodology, and also estimation of heritability for interested characters for culture would be profitable.

4.4. General Remarks

A progressive decrease in genetic variability is observed in successive introductions of L. vannamei into Cuba since 2003 and henceforth. It means that the first broodstock should be maintained and combined with others or new stock. Relatedness coefficients support the use of this stock and highlight the fact that the second stock should be carefully used in crossings because of their tendency to positive values.

However, the sustainability of the culture has been achieved, as well as maintenance of no viral diseases that have caused considerable losses in Latin America [31]. In that way, if surveillance program [32], starts with the premise of introducing SPF animals (Shrimp Pathogens Free) for health maintenance in cultures, a strategy of genetic management with microsatellite markers will warranty an adequate productive performance that does not drive to an irreversible endogamy. In agreement with [2], the maintenance of a common origin, besides of health warranty ever crossing stocks genetically different, avoids phenomena as outbreeding depression [33–35] that could damage productive yields in a short time.

Acknowledgments

Many thanks are due to workers and technicians from the Postlarvae Production Hatchery “YAGUACAM,” in Cienfuegos, Cuba, for technical support in maintenance of lines, especially to its director, Ms. Angela Moreno. Special thanks are due to Dr. Vicente Berovides for his critical revision of the text and also to MSc. Román Machado for his comments. The authors also thank the anonymous referees for improving the document with their comments and remarks. They also want to express their gratitude to Ms. Mercedes Escobar for editing and improving English from the original manuscript.

References

R. Tizol, B. Jaime, R. Laria et al., “Introduction of Pacific White Shrimp L. vannamei in Cuba. Quarantine I step (In Spanish). Ocean Docs,” http://hdl.handle.net/1834/3588, 2004.
View at: Google Scholar
Y. J. Borrell, G. Espinosa, E. Vazquez, J. A. Sánchez, and G. Blanco, “Genetic variability of microsatellite loci in the two first stocks of Litopenaeus vannamei introduced to Cuba for aquaculture,” Revista de Investigaciones Marinas, vol. 27, no. 3, pp. 237–244, 2006 (Spanish).
View at: Google Scholar
R. Machado, Assessment of genetic variability in two lots of white shrimp, Litopenaeus vannamei (Boone, 1931) introduced to Cuba, M.S. thesis, International Fisheries Management. Department of Aquatic Biosciences. Norwegian College of Fishery Science. University of Tronso, Norway, Oslo, 2006.
A. Artiles, I. Rodríguez, A. Pérez, L. Pérez, and G. Espinosa, “Low genetic variability in the fifth introduction of Litopenaeus vannamei in Cuba, as estimated with microsatellite markers,” Biotecnología Aplicada, vol. 28, pp. 142–146, 2011.
View at: Google Scholar
D. K. García, M. A. Faggart, L. Rhoades et al., “Genetic diversity of cultured Penaeus vannamei shrimp using three molecular genetic techniques,” Molecular Marine Biology and Biotechnology, vol. 3, no. 5, pp. 270–280, 1994.
View at: Google Scholar
S. Sunden and K. Davis, “Evaluation of genetic variation in a domestic population of Penaeus vannamei (Boone): a comparison with three natural populations,” Aquaculture, vol. 97, no. 2-3, pp. 131–142, 1991.
View at: Google Scholar
R. A. Dunham, K. Majumdar, E. Hallerman et al., “Review of the status of aquaculture genetics,” in Aquaculture in the Third Millennium. Technical Proceedings of the Conference on Aquaculture in the Third Millennium (February 2000, Bangkok, Thailand), R. P. Subasinghe, P. Bueno, M. J. Phillips, C. Hough, S. E. McGladdery, and J. R. Arthur, Eds., pp. 137–166, NACA/FAO, 2001.
View at: Google Scholar
G. Hulata, “Genetic manipulations in aquaculture: a review of stock improvement by classical and modern technologies,” Genetica, vol. 111, no. 1-3, pp. 155–173, 2001.
View at: Publisher Site | Google Scholar
Z. J. Liu and J. F. Cordes, “DNA marker technologies and their applications in aquaculture genetics,” Aquaculture, vol. 238, no. 1–4, pp. 1–37, 2004.
View at: Publisher Site | Google Scholar
G. M. Wolfus, G. K. García, and A. Alcivar-Warren, “Application of the microsatellite technique for analyzing genetic diversity in shrimp breeding programs,” Aquaculture, vol. 152, no. 1–4, pp. 35–47, 1997.
View at: Publisher Site | Google Scholar
P. Cruz, C. H. Mejía-Ruiz, R. Pérez-Enriquez, and A. M. Ibarra, “Isolation and characterization of microsatellites in Pacific white shrimp Penaeus (Litopenaeus) vannamei,” Molecular Ecology Notes, vol. 2, no. 3, pp. 239–241, 2002.
View at: Publisher Site | Google Scholar
R. Peakall and P. E. Smouse, “GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research,” Molecular Ecology Notes, vol. 6, no. 1, pp. 288–295, 2006.
View at: Publisher Site | Google Scholar
J. Goudet, “FSTAT, a program to estimate and test gene diversities and fixation indexes. (version 2.9.3),” Journal of Heredity, vol. 86, pp. 485–486, 2002.
View at: Google Scholar
S. Wright, “The interpretation of population structure by F-statistics with special regards to systems of mating,” Evolution, vol. 19, pp. 395–420, 1965.
View at: Google Scholar
F. Rousset, Genepop 4.1 for Windows/Linux/Mac OS X, 2008.
B. Rannala and J. L. Mountain, “Detecting immigration by using multilocus genotypes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 94, no. 17, pp. 9197–9201, 1997.
View at: Publisher Site | Google Scholar
S. Piry, A. Alapetite, J. M. Cornuet, D. Paetkau, L. Baudouin, and A. Estoup, “GeneClass 2: a software for genetic assignment and first-generation migrant detection,” Journal of Heredity, vol. 95, no. 6, pp. 536–539, 2004.
View at: Publisher Site | Google Scholar
D. C. Quelle and K. F. Goodnight, “Estimating relatedness using genetic markers,” Evolution, vol. 43, no. 2, pp. 258–275, 1989.
View at: Google Scholar
GraphPad InStat version 5.04 for Windows, GraphPad Software, La Jolla California USA, http://www.graphpad.com.
R. Valles-Jiménez, P. Cruz, and R. Pérez-Enriquez, “Population genetic structure of Pacific white shrimp (Litopenaeus vannamei) from Mexico to Panama: microsatellite DNA variation,” Marine Biotechnology, vol. 6, no. 5, pp. 475–484, 2005.
View at: Google Scholar
P. Cruz, A. M. Ibarra, H. Mejia-Ruiz, P. M. Gaffney, and R. Pérez-Enríquez, “Genetic variability assessed by microsatellites in a breeding program of pacific white shrimp (Litopenaeus vannamei),” Marine Biotechnology, vol. 6, no. 2, pp. 157–164, 2004.
View at: Google Scholar
R. Pérez-Enríquez, F. Hernández-Martínez, and P. Cruz, “Genetic diversity status of White shrimp Penaeus (Litopenaeus) vannamei broodstock in Mexico,” Aquaculture, vol. 297, no. 1–4, pp. 44–50, 2009.
View at: Publisher Site | Google Scholar
N. Bierne, I. Beuzart, V. Vonau, F. Bonhomme, and E. Bedier, “Microsatellite—associated heterosis in hatchery—propagated stokcs of the srimp Penaeus stylirostris,” Aquaculture, vol. 184, pp. 203–219, 2000.
View at: Google Scholar
Z. Xu, J. H. Primavera, L. D. de la Pena et al., “Genetic diversity of white and cultured black tiguer shrimp (Penaeus monodon) in the Philippines using microsatellites,” Aquaculture, vol. 199, pp. 13–40, 2001.
View at: Google Scholar
E. Luvesuto, P. Dominguez de Freitas, and P. M. Galetti Junior, “Genetic variation in a closed line of the white shrimp Litopenaeus vannamei (Penaeidae),” Genetics and Molecular Biology, vol. 30, no. 4, pp. 1156–1160, 2007.
View at: Google Scholar
J. A. H. Benzie, “Population genetic structure in penaeid prawns,” Aquaculture Research, vol. 31, no. 1, pp. 95–119, 2000.
View at: Publisher Site | Google Scholar
G. Espinosa López, R. Díaz Fernández, U. Becker Zúñiga et al., “Population analysis of Cuban white shrimp Litopenaeus schmitti using allozymes as genetic markers,” Revista de Investigaciones Marinas, vol. 24, no. 1, pp. 11–16, 2003 (Spanish).
View at: Google Scholar
J. Soto-Hernández and J. M. Grijalva-Chon, “Genetic differentiation in hatchery strains and wild white shrimp Penaeus (Litopenaeus) vannamei (Boone, 1931) from northwest Mexico,” Aquaculture International, vol. 12, no. 6, pp. 593–601, 2004.
View at: Publisher Site | Google Scholar
V. Sbordoni, E. de Matthaeis, M. Cobolli Sbordoni, G. La Rosa, and M. Mattoccia, “Bottleneck effects and the depression of genetic variability in hatchery stocks of Penaeus japonicus (Crustacea, Decapoda),” Aquaculture, vol. 57, no. 1–4, pp. 239–251, 1986.
View at: Google Scholar
Y. J. Borrell, J. Álvarez, E. Vázquez et al., “Applying microsatellites to the management of farmed turbot stocks (Scophthalmus maximus L.) in hatcheries,” Aquaculture, vol. 241, no. 1–4, pp. 133–150, 2004.
View at: Publisher Site | Google Scholar
A. Artiles, M. Rubio, E. Gonzalez, R. Laria, and R. Silveira, “Crustacean virus of obligatory declaration by the OIE performance in cultured Litopenaeus vannamei in Cuba from 2003 to 2009,” Ocean Docs Digital Repository, vol. 28, no. 1, pp. 12–18, 2011 (Spanish).
View at: Google Scholar
R. Silveira-Coffigny, “Aquaculture health in Cuba,” in Proceedings of the 1st Meeting of Interamerican OIE Comité for Aquatic Animals, Panama, Panama City, January 2006.
View at: Google Scholar
S. Edmands, “Heterosis and outbreeding depression in interpopulation crosses spanning a wide range of divergence,” Evolution, vol. 53, no. 6, pp. 1757–1768, 1999.
View at: Google Scholar
M. Keller, J. Kollmann, and P. J. Edwards, “Genetic introgression from distant provenances reduces fitness in local weed populations,” Journal of Applied Ecology, vol. 37, no. 4, pp. 647–659, 2000.
View at: Publisher Site | Google Scholar
S. Granier, C. Audet, and L. Bernatchez, “Heterosis and outbreeding depression between strains of young-of the-year brook trout (Salvelinus fontinalis),” Canadian Journal of Zoology, vol. 89, pp. 190–198, 2011.
View at: Google Scholar

Copyright

Copyright © 2012 Anna Pérez-Beloborodova et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1740

Downloads

1381

Citations