Genetic Variation in Wild and Cultured Populations of the Sea Squirt Halocynthia roretzi Inferred from Microsatellite DNA Analysis
Genetic Variation in Wild and Cultured Populations of the Sea Squirt Halocynthia roretzi Inferred from Microsatellite DNA Analysis
Fisheries and aquatic sciences. 2012. Jun, 15(2): 151-155
Copyright ©2012, The Korean Society of Fisheries and Aquatic Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons. org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : February 06, 2012
  • Accepted : March 20, 2012
  • Published : June 30, 2012
Export by style
Cited by
About the Authors
Hyon-Sob, Han
West-Sea Fisheries Research Institute, National Fisheries Research and Development Institute (NFRDI), Incheon 400-420, Korea
Bo-Hye, Nam
Biotechnology Research Center, NFRDI, Busan 619-902, Korea
Jung-Ha, Kang
Biotechnology Research Center, NFRDI, Busan 619-902, Korea
Yi-Kyoung, Kim
Biotechnology Research Center, NFRDI, Busan 619-902, Korea
Young-Ju, Jee
Biotechnology Research Center, NFRDI, Busan 619-902, Korea
Young-Baek, Hur
Aquaculture Environmental Institute, NFRDI, Tongyeong 650-943, Korea
Moongeun, Yoon
Yangyang Salmon Station, Korea Fisheries Resources Agency, Yangyang 215-851, Korea

We used nine microsatellite DNA markers to estimate genetic variation among wild and cultured populations of the sea squirt Halocynthia roretzi . The loci were polymorphic, with 6-32 alleles, and allelic richness ranged from 6.0 to 26.1 in each population. The wild and the cultured populations had similar mean heterozygosities ( HO and HE ), allele numbers, and allelic richness. One cultured population with softness syndrome had a lower mean in the observed heterozygosity ( HO = 0.57) and higher mean inbreeding coefficient ( F IS = 0.261) than any other populations. This suggests that the loss of genetic variation in the diseased population might be due to increased inbreeding. A neighbor-joining tree and pairwise population estimates of F ST showed moderate genetic differentiation between the wild and the cultured populations. Additionally, the softness syndrome population was genetically divergent from wild populations, but it was genetically close to the cultured populations.
The natural distribution of the solitary ascidian Halocynthia roretzi (sea squirt) is restricted essentially to the waters around Korea and Japan, where it inhabits rocky littoral surfaces. H. roretzi has received considerable attention due to its high commercial importance in Korea. The Korean aquaculture industry produced 9,300 tonnes of H. roretzi in 2007, but this decreased rapidly to less than 1,000 tonnes in 2010 (Ministry of Maritime Affairs and Fisheries, 2010); this decrease has been the result of an unexplained syndrome, which has been termed “softness syndrome”. To date, the syndrome has caused the loss of more than $30 million to Korean producers.
The sea squirt are cultivated by the hanging culture method (Rho et al., 1993): seeds from a hatchery are attached to a vertical rope that is hung from a longline, much as oysters are cultivated (Kang et al., 1980). Hatcheries in Korea often use seeds that have been introduced from Iwate in Japan, but the details of the origins and breeding records of the introduced sea squirt are unknown. Also, large numbers of the cultured sea squirt in Korea have died from an unexplained syndrome over the last 20 years (Cho et al., 2008). Each dying Sea Squirt shows thinning of the tunica, and the meat becomes quite soft (thus the name “softness syndrome”). To manage and improve the sustainable production of this sea squirt resource in Korea, it is vital to identify the causative factors responsible for the
PPT Slide
Lager Image
Geographic locations of sampling sites. SAMW, the wild populations from Samcheok; TOYW, the wild populations from Tongyeong; TOYC, the cultured populations from Tongyeong; IWTC, cultured population from Japan; SOFT, the cultured population from Tongyeong with softness syndrome.
softness syndrome.
A major concern has been focused on genetic changes within cultured populations may threaten the integrity of wild populations (Allendorf and Ryman, 1987). Cultivation to enhance the stock has also been identified as one of the factors responsible for the depletion of wild populations, due to inbreeding, competition, and the spread of disease (Allendorf and Phelps, 1980). Furthermore, genetic variation is beneficial and important for the long-term survival of natural populations, because variation can provide the population with the ability to adapt to changing environmental conditions (Frankel and Soulé, 1981 ).
Genetic variation within and between populations provides a potential resource for future adaptations (Hurt and Hedrick, 2004). In recent years, molecular markers have been used to observe the genetic variation in sea squirt populations to study phylogeography and population genetics (e.g., Ben-Shlomo et al., 2001, 2006; Tarjuelo et al., 2004; López-Legentil and
Sampling locations, date of collection, the number of samples (N), observed (HO) and expected (HE) heterozygosities, inbreeding coefficient (FIS), and mean number of alleles (AM) used for Halocynthia roretzi msDNA analysis
PPT Slide
Lager Image
Sampling locations, date of collection, the number of samples (N), observed (HO) and expected (HE) heterozygosities, inbreeding coefficient (FIS), and mean number of alleles (AM) used for Halocynthia roretzi msDNA analysis
Turon, 2006). To date, no genetic study in H. roretzi by molecular markers has been reported.
In the present study, we investigated genetic variability between the wild and the cultured populations of H. roretzi collected in Korea and Japan by analyzing nine microsatellite loci.
Materials and Methods
- Samples and DNA extraction
In 2006 and 2007, sea squirt were collected from two wild populations at Samcheok (abbreviated as SAMW) and Tongyeong (TOYW), Korea, respectively, and three cultured populations: a population from Tongyeong (TOYC), a population from Tongyeong with softness syndrome (SOFT), and a population from Iwate, Japan (IWTC) ( Table 1 , Fig. 1 ). The TOYC population, which was introduced in 1989 from Iwate in Japan, originated from third generation spawners (body size > 10 cm) that had been reared continuously. The SOFT and IWTC populations had also been reared continuously, but no information is available about the origin of these populations. Samples of muscle from internal organs were collected from 426 individuals from the five populations. The muscle samples were preserved in 100% ethanol at room temperature until the DNA was extracted. Genomic DNA was extracted from the stored muscle tissues by an automated DNA extraction system (MagExtractor MFX-6100; Toyobo, Japan ) with a MagExtractor Genomic DNA Purification Kit (Toyobo). Extracted DNA was dissolved in 80 μL of Tris-EDTA buffer (TE; 10 mM Tris-HCl and 1 mM EDTA, pH 8.0) until genotyping.
- Microsatellite DNA analysis
Nine polymorphic microsatellite loci were examined by PCR amplification by primers previously developed for this species (Han et al., 2009). The forward primer of each pair was labeled with one of the following fluorescent dyes: 6FAM, NED or HEX. PCR amplification was performed with a PTC-200 thermocycler (MJ Research, Inc., Waltham, MA, USA ) and the following reaction conditions: initial denaturation for 5 min at 95℃, 40 cycles of 30 s at 95℃, 45 s at a primer-specific annealing temperature, 45 s at 72℃, and a final extension step at 72℃ for 30 min. The specific annealing temperature of each primer set is given in our previous paper (Han et al., 2009). The PCR products were analysed using an ABI 3130 XL genetic analyzer (Applied Biosystems, Foster City, CA, USA ) to screen for polymorphisms.
- Analysis of population genetics data
The inbreeding coefficient F IS was calculated from an analysis of variance following Weir and Cockerham (1984) by the GENEPOP program. Observed ( H O ) and expected ( H E ) heterozygosities per locus were calculated for each population using the programs FSTAT 2.9.3 (Goudet, 2001) and GENEPOP, version 3.4 (Rousset and Raymond, 1995). Pairwise values of F ST were used to estimate genetic differentiation between pairs of populations according to Slatkin (1995) by the program FSTAT 2.9.3 (Goudet, 2001). The significance of F ST values was tested via a permutation method (28,000 permutations).
A matrix of genetic distances, based on the set of allele frequencies in different populations, was constructed according to Nei (1972). A neighbor-joining (NJ) tree was constructed for each replicated matrix by bootstrapping 1,000 replications with NEIGHBOR in PHYLIP version 3.5 (Felsenstein, 1993). The consensus tree was generated following the 50% majority rule by CONSENSUS, also in the PHYLIP package (available at html ).
Results and Discussion
To our knowledge, this is the first reported study on the population genetics of H. roretzi in Korea to estimate current levels of genetic variation by nine microsatellite DNA markers. In total, 203 alleles were observed for all nine loci, ranging from 12 alleles at Haro 26 to 46 at Haro 15, in all populations. Allele sizes ranged from 125 to 232 bp across all microsatellite loci. As shown in Table 1 , the mean number of observed alleles ranged from 10.0 (SOFT) to 12.4 (TOYC), the mean H O per population ranged from 0.57 (SOFT) to 0.79 (SAMW), and the mean H E ranged from 0.76 (IWTC) to 0.79 (TOYW). These results indicated that the levels of genetic variation were similar in each of the five populations. No meaningful difference in heterozygosity was found between the wild and the cultured populations, except SOFT, suggesting that the cultured populations of H. roretzi showed similar levels of genetic variation to those of the wild populations, at least in terms of the observed levels of heterozygosity. The mean observed heterozygosity of SOFT was the lowest among the different populations. Additionally, SOFT showed a large disparity between H O and H E . However, we noted that the mean H O values were lower than the H E values for all populations, except SAMW, which might indicate a deficit in heterozygotes. The values of F IS and the evaluations of Hardy-Weinberg equilibrium are shown in Table 1 . The SAMW population had the lowest estimate of mean values of inbreeding (mean F IS = -0.021) of the five populations. TOYW, TOYC, and IWTC exhibited moderate mean F IS values, which ranged from 0.043 to 0.139. In contrast, the mean F IS was higher in SOFT than in the other populations.
Most of the sea suirt cultured in Tongyeong area have died during the last 20 years from the “softness syndrome” disease which has had a detrimental effect on the Korean aquaculture industry. Cho et al. (2008) discriminated normal and dying sea squirt by their differential patterns of gene expression. Low genetic variability in the cultured populations, which is caused by increased levels of inbreeding due to the use of a small effective number of parents, may generally lower resistance to various diseases. Although the mean H O and F IS values for the microsatellite loci did not reveal pronounced differences between the wild and the cultured populations, the results of our study confirmed that the diseased population did show less genetic variation. The population with “softness syndrome” had lower H O and higher F IS values at each locus compared with the healthy populations. Moreover, the disparity between the H O and H E in the diseased population was much higher than that in the healthy populations, suggesting that the diseased population was considerably inbred.
Pairwise population estimates of F ST are presented in Table 2 . The F ST values showed significant differences between the wild and the cultured populations, suggesting that genetic differentiation had occurred between these populations. Significant differentiation was inferred for all pairs of the cultured
PPT Slide
Lager Image
Neighbor-joining tree showing genetic distance according to Nei (1972) between five sea squirt populations. The inset shows the consensus phenogram with nodal bootstrap values of 1,000 replicates (%). SAMW, the wild populations from Samcheok; TOYW, the wild populations from Tongyeong; TOYC, the cultured populations from Tongyeong; IWTC, the cultured population from Japan; SOFT, the cultured population from Tongyeong with softness syndrome.
populations, except TOYC and IWTC, and for pairs of wild populations. The highest values of F ST were between SOFT and the wild populations. The sampled populations were clustered by the NJ method ( Fig. 2 ). All the wild populations were separated from the cluster of the cultured populations with 100% bootstrap support. These findings indicate that genetic differentiation has occurred between the wild and the cultured populations. The SOFT population was genetically distinct from the other populations, as suggested by the estimates of F ST . Topology of the NJ tree showed that TOYC (Korea) was genetically close to IWTC (Japan), by frequent import of the parent stocks.
Pairwise population estimates of F ST revealed highly significant genetic differentiation between the eastern (SAMW) and
FST values between samples (below diagonal) and probability of differentiation with P-value in FST estimate (above diagonal)SAMW, wild populations at Samcheok; TOYW, wild populations at Tongyeong; TOYC, cultured populations from Tongyeong; IWTC, cultured population from Japan; SOFT, cultured population from Tongyeong with softness syndrome.*Significant support forFSTvalues (P< 0.05).
PPT Slide
Lager Image
FST values between samples (below diagonal) and probability of differentiation with P-value in FST estimate (above diagonal) SAMW, wild populations at Samcheok; TOYW, wild populations at Tongyeong; TOYC, cultured populations from Tongyeong; IWTC, cultured population from Japan; SOFT, cultured population from Tongyeong with softness syndrome. *Significant support for FST values (P < 0.05).
southern coastal (TOYW) populations. The NJ tree also indicated that this wild population pair was genetically distinct. These results suggested low levels of gene flow between populations in different regions and long-term isolation between regional populations. Grosberg (1987) and Davis and Butler (1989) found that the larvae of colonial sea squirt were only dispersed over short distances. Thus, our findings suggest that genetic interactions between populations of H. roretzi only occur over short distances. This probably reflects the limited duration of the free-swimming stage before it settles into the habitat. Even though our findings indicate strongly restricted gene flow between populations, pairwise tests revealed that the Tongyeong cultured population is very closely related to the Iwate cultured population, possibly due to human-mediated transportation of stocks. As noted in the Introduction of this study, to produce a large number of offspring of H. roretzi at the Tongyeong aquaculture sites, many parent and seed have been introduced from Iwate during the last 20 years.
In conclusion, our results demonstrate that genotyping with microsatellite allele frequencies can be applied to the identification of stock and the estimation of genetic variability among and within wild and/or cultured populations. In particular, the results of the analysis suggested that genetic variation had decreased in a cultured population of the sea squirt H. roretzi with softness syndrome and that this population had differentiated genetically from the healthy populations. The study of genetic differences among the wild and the cultured populations will help in the development of a sound strategy for the management of the genetic stock programs and the subsequent ecological monitoring of sea squirt communities.
This research was supported by a grant from the NFRDI (RP-2008-BT-006)
Allendorf FW , Phelps SR. 1980 Loss of genetic variation in a hatchery stock of cutthroat trout. Trans Am Fish Soc 109 537 - 543
Allendorf FW , Ryman N. 1987 Genetic management of hatchery stocks. In: Population Genetics and Fishery Management. Ryman N and Utter F, eds. University of Washington Press Seattle, WA, US 141 - 159
Ben-Shlomo R , Douek J , Rinkevich B. 2001 Heterozygote deficiency and chimerism in remote populations of a colonial ascidian from New Zealand. Mar Ecol Prog Ser 209 109 - 117
Ben-Shlomo R , Paz G , Rinkevich B. 2006 Postglacial-period and recent invasions shape the population genetics of Botryllid ascidians along European Atlantic Coasts. Ecosystems 9 1118 - 1127
Cho HK , Nam BH , Kong HJ , Han HS , Hur YB , Choi TJ , Choi YH , Kim WJ , Cheong J. 2008 Identification of softness syndromeassociated candidate genes and DNA sequence variation in the sea squirt, Halocynthia roretzi. Mar Biotechnol 10 447 - 456
Davis AR , Butler AJ. 1989 Direct observations of larval dispersal in the colonial ascidian (Podoclavella moluccensis Sluiter): evidence for closed populations. J Exp Mar Biol Ecol 127 189 - 203
Felsenstein J. 1993 PHYLIP (Phylogeny Interference Package) version 3.5c. Department of Genetics Washington State University Seattle, WA, US.
Goudet J 2001 FSTAT, a program to estimate and test gene diversities and fixation indices (version 2. 9. 3) [Internet]. Universite de Lausanne Suisse, CH .
Grosberg RK. 1987 Limited dispersal and proximity-dependent mating success in the colonial ascidian Botryllus schlosseri. Evolution 41 372 - 384
Han HS , Nam BH , Kang JH , Choi TJ , Kim WJ , Hur YB , Yoon M. 2009 Isolation and characterization of polymorphic microsatellite DNA markers in the sea squirt, Halocynthia roretzi. Mol Ecol Resour 9 648 - 650
Hurt C , Hedrick P. 2004 Conservation genetics in aquatic species: general approaches and case studies in fishes and springsnails of arid lands. Aquat Sci 66 402 - 413
Kang PA , Kim Y , Yoon DS. 1980 Studies on the hanging culture of oyster, Crassostrea gigas, in the Korean coastal waters. 4. On the fouling organisms associated with culturing oysters at the oyster culture farms in Chungmu. Bull Fish Res Dev Agency Busan 25 29 - 34
López-Legentil S , Turon X. 2006 Population genetics, phylogeography and speciation of Cystodytes (Ascidiacea) in the western Mediterranean Sea. Biol J Linn Soc 88 203 - 214
Ministry of Maritime Affairs and Fisheries. 2010 Statistical Database for Fisheries Production [Internet]. Korea National Statistical Office Daejeon, KR
Nei M. 1972 Genetic distance between populations. Am Nat 106 283 - 292
Rho YG , Lee YH , Park MW. 1993 The environmental factors affecting mortality of cultured sea squirt, Halocynthia roretzi (Drasche). Bull Nat'l Fish Res Dev Agency (Korea) (47) 145 - 164
Rousset F , Raymond M. 1995 Testing heterozygote excess and deficiency. Genetics 140 1413 - 1419
Slatkin M. 1995 A measure of population subdivision based on microsatellite allele frequencies. Genetics 139 457 - 462
Tarjuelo I , Posada D , Crandall A , Pascual M , Turon X. 2004 Phylogeography and speciation of colour morphs in the colonial ascidian Pseudodistoma crucigaster. Mol Ecol 13 3125 - 3136
Weir BS , Cockerham CC. 1984 Estimating F-statistics for analysis of population structure. Evolution 38 1358 - 1370