Development, characterization, and inheritance of 113 novel EST-SSR markers in the Pacific oyster (Crassostrea gigas)

A total of 113 novel EST-derived simple sequence repeat (EST-SSR) markers were developed in the Pacific oyster (Crassotrea gigas). Polymorphisms of these markers were evaluated in a wild population of 30 individuals. The number of alleles per locus ranged from 2 to 27 with an average of 6.3, and the observed and expected heterozygosities varied from 0 to 0.9667 and from 0.0333 to 0.9701, respectively. Mendelian segregations were tested for 24 of the markers that were polymorphic in one family produced by single-pair mating. Null alleles were discovered at four loci. Nine tests of segregation ratios revealed significant departures from expected Mendelian ratios. As a useful addition to the collection of the microsatellites that are now available for C. gigas, these EST-SSR markers will help the advance in investigation of QTL mapping and genetic diversity in this species.     

Genetic and epigenetic variation in mass selection populations of Pacific oyster Crassostrea gigas

Selective breeding often produces an improvement in phenotype. Much of the phenotypic change within a species is a consequence of genetic variation. However, there is growing evidence for phenotypic change even in the absence of DNA sequence polymorphisms, termed epigenetic variation. This study’s goal was to investigate the genetic and epigenetic variation in the mass selection populations of the Pacific oyster (Crassostrea gigas), determine if any correlation exists between the genetic and epigenetic variations. This can serve as a first step in investigating the potential role epigenetic variations have in selective breeding. Amplified fragment length polymorphism analysis and methylation-sensitive amplified polymorphism methodology were used to monitor genetic and epigenetic variation in two populations (the base stock and the third selected generation) from a mass selection line in the Pacific oyster. The correlation between genetic and epigenetic variation was evaluated by Co-Inertia Analysis. The genetic difference was mainly found in the gene frequency shift revealed by the F _ST value (0.0151, P < 0.01) and no significant reduction in genetic diversity was detected. The percentage of methylation in C. gigas was 26.4 %. No significant difference was observed on the average state of methylation, but a few bands showed different frequencies between the two populations. Co-Inertia Analysis revealed a significant association between the genetic and epigenetic profiles (P < 0.01).     

Microsatellite analysis of gynogenetic families in the Pacific oyster, Crassostrea gigas

Five families of gynogenetic diploid Pacific oyster (Crassostrea gigas) were induced by inhibiting the second polar body in meiotic cell division of eggs fertilized with UV-irradiated sperm. Segregation patterns of eight microsatellite loci were investigated in the gynogenetic diploid offspring; the proportion of heterozygous progeny was used to estimate microsatellite-centromere (M-C) distances. Mendelian inheritance was confirmed for the eight loci by examining the genotypic segregation in the control crosses. Three of the eight microsatellite loci showed the existence of null alleles in four control crosses. All gynogenetic offspring only possessed the alleles of the mother, indicating 100% success level for the five families. The M-C recombination frequency estimates ranged from 0.62 to 0.77 (0.72 mean), comparable to those in the oyster based on allozyme markers and suggesting that meiotic gynogenesis does not appear to be a very efficient inbreeding method in the oyster. Recombination frequencies observed were often higher than the theoretical maximum of 0.67, indicating the existence of positive interference after a single chiasma formation in some chromosomes. Information on the positions of centromeres in relation to the microsatellite loci will represent a contribution toward assembly of genetic maps in C. gigas.     

Development and Characterization of EST-SSR Markers in the Eastern Oyster <Emphasis Type="Italic">Crassostrea virginica</Emphasis>

Simple sequence repeat (SSR) markers were developed from expressed sequence tags (ESTs) in the eastern oyster (Crassostrea virginica). ESTs of the eastern oyster were downloaded from GenBank and screened for SSRs with at least eight units of dinucleotide or five units of tri-, tetra-, penta-, and hexa-nucleotide repeats. The screening of 9101 ESTs identified 127 (1.4%) SSR-containing sequences. Primers were designed for 88 SSR-containing ESTs with good and sufficient flanking sequences. Polymerase chain reaction (PCR) amplification was successful for 71 primer pairs, including 19 (27%) pairs that amplified fragments longer than expected sizes, probably due to introns. Sixty-six pairs that produced fragments shorter than 800 bp were screened for polymorphism in five oysters from three populations via polyacrylamide gels, and 53 of them (80%) were polymorphic. Fifty-three polymorphic SSRs were labeled and genotyped in 30 oysters from three populations via an automated sequencer. Five of the SSRs amplified more than two fragments per oyster, suggesting locus duplication. The remaining 48 SSRs had 2 alleles per individual, including 11 with null alleles. In the 30 oysters analyzed, the SSRs had an average of 9.3 alleles per locus, ranging from 2 to 24. Forty-three loci segregated in a family with 100 progeny, with nine showing significant deviation from Mendelian ratios (three after Bonferroni correction). Seventy percent of the loci were successfully amplified in C. rhizophorae and 34% in C. gigas. This study demonstrates that ESTs are valuable resources for the development of SSR markers in the eastern oyster, and EST-derived SSRs are more transferable across species than genomic SSRs.     

Estimates of Heritability for Growth and Shell Color Traits and Their Genetic Correlations in the Black Shell Strain of Pacific Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis>

The Pacific oyster Crassostrea gigas has been introduced widely and massively and became an economically important aquaculture species on a global scale. We estimated heritabilities of growth and shell color traits and their genetic correlations in black shell strain of C. gigas. Analyses were performed on 22 full-sib families in a nested mating design including 410 individuals at harvest (24 months of age). The parentage assignment was inferred based on four panels of multiplex PCR markers including 10 microsatellite loci and 94.9% of the offspring were unambiguously assigned to single parent pairs. The Spearman correlation test (r = ? 0.992, P < 0.001) demonstrated the high consistency of the shell pigmentation (SP) and L* and their same efficacy in shell color measurements. The narrow-sense heritability estimated under the animal model analysis was 0.18 ± 0.12 for shell height, 0.25 ± 0.16 for shell length, 0.10 ± 0.09 for shell width, 0.42 ± 0.20 for total weight, 0.32 ± 0.18 for shell weight, and 0.68 ± 0.16 for L*, 0.69 ± 0.16 for shell pigmentation, respectively. The considerable additive genetic variation in growth and shell color traits will make it feasible to produce genetic improvements for these traits in selective breeding program. High genetic and phenotypic correlations were found among growth traits and among shell color traits. To optimize a selection strategy for both fast growth and pure dark shell strain of C. gigas, it is proposed to take both total weight and black shell as joint objective traits in selective breeding program. Our study offers an important reference in the process of selective breeding in black shell color stain of C. gigas and will facilitate to develop favorable breeding strategies of genetic improvements for this economically important strain.     

Development of gene-targeted SNP markers for genomic mapping in Pacific abalone Haliotis discus hannai Ino

Useful and novel DNA markers are needed for aquaculture genetics and breeding. In this study, we report the discovery and development of gene-targeted single nucleotide polymorphisms (SNPs) for genomic mapping in the Pacific abalone Haliotis discus hannai Ino. Single EST or EST-contigs from 66 genes that had positive BLASTx matches (E-value ≤ 1e-8) were used for polymerase chain reaction (PCR) amplification. PCR products from the two parents of one mapping family were directly sequenced, and 83 SNP loci were found from 17 genes. Allele-specific PCR (AS-PCR) was developed and optimized for genotyping of 11 SNP loci in 120 progeny of the mapping family. Nine of the loci conformed to the expected Mendelian ratio of 1:1 based on the χ² test (P > 0.05) and could potentially be used for linkage map construction. Our data also indicate that the sequencing of two parents may be a practical strategy for the discovery of informative SNPs for linkage mapping in a particular mapping population.     

Comparative genetic diversity of wild and hatchery-produced Pacific oyster (Crassostrea gigas) populations in Korea using multiplex PCR assays with nine polymorphic microsatellite markers

The Pacific oyster, Crassostrea gigas, is the most important and valuable commercial fishery species in Korea. Its farming started 20 years ago and is still rapid expansion in Korea. In this study, to maintain the genetic diversity of this valuable marine resource, possible genetic similarity and differences between the wild population and hatchery population in Tongyeong, Korea were accessed using multiplex assays with nine highly polymorphic microsatellite loci. A total of 250 different alleles were found over all loci. Despite a long history of hatchery practices, very high levels of polymorphism (mean alleles = 22.89 and mean heterozygosity = 0.92) were detected between the two populations. No statistically significant reductions were found in heterozygosity or allelic diversity in the hatchery population compared with the wild population. However, significant genetic heterogeneity was found between two populations. These results provide no evidence to show that hatchery practice of Pacific oyster in Korea has significantly affected the genetic variability of the hatchery stock. Although further studies are needed for comprehensive determinations of the hatchery and wild populations with increased number of Pacific oyster sample collections, information on the genetic variation and differentiation obtained in this study can be applied for genetic monitoring of aquaculture stocks, genetic improvement by selective breeding and designing of more efficient conservation management guidelines for these valuable genetic materials.     

Genetic Variation and Breeding Signature in Mass Selection Lines of the Pacific Oyster (Crassostrea gigas) Assessed by SNP Markers

In breeding industries, a challenging problem is how to keep genetic diversity over generations. To investigate genetic variation and identify breeding signatures in mass selected lines of Pacific oyster (Crassostrea gigas), three sixth-generation selected lines and four wild populations were assessed using 103 single nucleotide polymorphism (SNP) markers. The genetic diversity data indicated that the selected lines exhibited a significant reduction in the observed heterozygosity and observed number of alleles per locus compared with the wild populations (P≤0.05), indicating the selected lines tended to lose genetic diversity contrasted with the wild populations. The unweighted pair-group method with arithmetic mean (UPGMA) analysis showed that the wild populations and selected lines were not separated into two groups. Using four outlier tests, a total of 17 loci were found under selection at two levels. The global outlier detection suggested that 4 common outlier loci were subject to selection using both the hierarchical island model and Bayesian likelihood approaches. At regional level, 3 SNPs were detected as outlier using at least two outlier tests and one outlier SNP (CgSNP309) was overlapped in the two wild-selected population comparisons. The candidate outlier SNPs provide valuable resources for future association studies in C. gigas.     

Association between polymorphism in the insulin receptor-related receptor gene and growth traits in the Pacific oyster Crassostrea gigas

Insulin receptor-related receptor (IRR) is an orphan receptor tyrosine kinase of the insulin receptor family, and involved in the growth and reproduction processes of the Pacific oyster Crassostrea gigas. Polymorphisms of the IRR gene were evaluated for associations with growth performance of 336 individuals in five families, and further confirmed in 206 individuals from three selectively bred strains for fast growth. Two of the six identified synonymous mutations (C.1996G > A and C.2110C > T) were significantly associated with growth performance in the families and strains. Five diplotypes were constructed based on the two growth-related SNPs, and diplotypes analysis revealed that D₃ (GGTT) might be the most advantageous diplotype for growth traits. The results suggest that two SNPs (C.1996G > A and C.2110C > T) in IRR gene are potentially associated with growth performance of C. gigas, and could serve as genetic markers for fast growth in oyster breeding.     

AFLP-Based Genetic Linkage Maps of the Pacific Oyster <Emphasis Type="Italic">Crassostrea gigas</Emphasis> Thunberg

Amplified fragment length polymorphisms (AFLPs) were used for genome mapping in the Pacific oyster Crassostrea gigas Thunberg. Seventeen selected primer combinations produced 1106 peaks, of which 384 (34.7%) were polymorphic in a backcross family. Among the polymorphic markers, 349 were segregating through either the female or the male parent. Chi-square analysis indicated that 255 (73.1%) of the markers segregated in a Mendelian ratio, and 94 (26.9%) showed significant (P < 0.05) segregation distortion. Separate genetic linkage maps were constructed for the female and male parents. The female framework map consisted of 119 markers in 11 linkage groups, spanning 1030.7 cM, with an average interval of 9.5 cM per marker. The male map contained 96 markers in 10 linkage groups, covering 758.4 cM, with 8.8 cM per marker. The estimated genome length of the Pacific oyster was 1258 cM for the female and 933 cM for the male, and the observed coverage was 82.0% for the female map and 81.3% for the male map. Most distorted markers were deficient for homozygotes and closely linked to each other on the genetic map, suggesting the presence of major recessive deleterious genes in the Pacific oyster.     

Context- and density-dependent effects of introduced oysters on biodiversity

Pacific oysters, Crassostrea gigas, have been introduced throughout much of the world, become invasive in many locations and can alter native assemblage structure, biodiversity and the distribution and abundance of other species. It is not known, however, to what extent their effects on biodiversity change as their cover increases, and how these effects may differ depending on the environmental context. Experimental plots with increasing cover of oysters were established within two estuaries in two different habitats commonly inhabited by C. gigas, (mussel-beds and mud-flats) and were sampled after 4 and 15 months. Within mud-flat habitats, macroscopic species living on or in the substratum increased in richness, Shannon–Wiener diversity and number of individuals with oyster cover. In mussel-bed habitats, however, these indices were unaffected by the cover of oysters except at one estuary after 15 months when species richness was significantly lower in plots with the greatest cover of oysters. Assemblage structure differed with oyster cover in mud-flats but not in mussel-beds, except at 100 % cover in one location and at one time. Within mud-flats at one location and time (of four total tests), assemblages became more homogenous with increasing cover of oysters leading to a significant decrease in β-diversity. These responses were primarily underpinned by the facilitation of several taxa including a grazing gastropod (Littorina littorea), an invasive barnacle (Austrominius modestus) and a primary producer (Fucus vesiculosus) with increasing cover of oysters. Although there were consistent positive effects of C. gigas on mud-flat biodiversity, effects were weak or negative at higher cover on mussel-beds. This highlights the need for the impacts of invasive species to be investigated at a range of invader abundances within different environmental contexts.     

Exploiting EST databases for the development and characterization of EST-SSRs in the Pacific oyster (Crassostrea gigas)

A total of 147 microsatellite-containing expressed sequence tags (ESTs) (3.63%) were detected from 4053 ESTs of the Pacific oyster (Crassostrea gigas) in GenBank. The average density of simple sequence repeats (SSRs) was 1 per 8.25 kb of EST after redundancy elimination. Dinucleotide repeat motifs appeared to be the most abundant type. Sixteen new polymorphic EST-SSRs were developed. The number of alleles per locus varied from 3 to 12, with an average of 5.9 alleles per locus. Marker transferability was tested on 2 other Crassostrea species, and 14 loci gave successful amplifications in both species. Twenty EST-SSRs were tested on 3 families of C. gigas for examination of inheritance mode of EST-SSRs. Thirty-five tests of segregation ratios revealed 5 significant departures from expected Mendelian ratios, 4 of which confirmed Mendelian expectations when accounting for the presence of null alleles. Null alleles were detected at 3 loci (15.0%) of the 20 loci, and the frequency of null alleles at EST-SSRs was lower than that in genomic SSRs in C. gigas. The results obtained in this study suggest that C. gigas EST-SSRs will complement the currently available genomic SSR markers and may be useful for comparative mapping, marker-assisted selection, and evolutionary studies.     

Fosmid library construction and end sequences analysis of the Pacific oyster,Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is globally distributed and is one of the most commercially and ecologically important marine organisms. However, little is known about the genome of this species. In this study, a C. gigas fosmid library was constructed that contains 459,936 clones with an average insert size of approximately 40 kb, representing 22.34-fold haploid genome equivalents. End sequencing generated 90,240 fosmid end sequences (FESs) with an average length of 384.27 base pairs (bp), covering approximately 2.58% of the Pacific oyster genome. The FESs were subsequently assembled and annotated, resulting in 6332 sequences with predicted open reading frames≥300 and 1,189,100 bp repeats. Furthermore, a total of 3200 microsatellite repeats were identified, and dinucleotide repeats were found to occur most abundantly, with AG and AAT being the most abundant repeat class of dinucleotides and trinucleotides. We also found that the repeat number was generally negatively proportional to the repeat element length. Microsatellites composition between the transcribed sequences and genomic sequences was shown to be different. Point mutations of microsatellite were non-random and underwent strong selection stress. Overall, a comprehensive sequence resource for the Pacific oyster was created, including annotated transposable elements, tandem repeats, protein coding sequences and microsatellites. These initial findings will serve as resources for further in-depth studies of physical mapping, gene discovery, microsatellite marker developing and evolution studies.     

Identification of candidate AFLP markers for shell color of the Pacific oyster (Crassostrea gigas) under artificial selection

The shell color of the Pacific oyster (Crassostrea gigas) is a desirable trait, but only a few genetic studies on shell color have been documented. Through successive selective breeding, four shell color variants of white (W), gold (G), black (B) and purple (P) C. gigas have been developed. The amplified fragment length polymorphism (AFLP) technique was used to scan the genomes of the four variants with different shell colors and one wild population (C) to identify candidate markers for shell polymorphism. Fifteen AFLP primer combinations were used, 1079 loci were scored as polymorphic loci, and the percentage of polymorphic bands was 95.5%. In the gold, white, black, purple and wild populations, the percentages of polymorphic loci were estimated to be 90.5% (G), 90.0% (W), 91.1% (B), 95.3% (P) and 93.2% (C); the expected heterozygosity values were 0.3115 (G), 0.3044 (W), 0.3102 (B), 0.3285 (P) and 0.3105 (C). The white shell variant was observed to have slightly lower genetic diversity than others, with a F_ST value of 0.1483. These results indicated that the four different shell color variants had high genetic diversity and that the genetic differentiation of populations mostly results from genetic diversity of individuals within populations. Furthermore, 11 outlier loci were considered candidate markers for shell color. This work provides new insights on relationships among color variants of C. gigas.     

Development of SNP‐genotyping arrays in two shellfish species

Use of SNPs has been favoured due to their abundance in plant and animal genomes, accompanied by the falling cost and rising throughput capacity for detection and genotyping. Here, we present in vitro (obtained from targeted sequencing) and in silico discovery of SNPs, and the design of medium‐throughput genotyping arrays for two oyster species, the Pacific oyster, Crassostrea gigas, and European flat oyster, Ostrea edulis. Two sets of 384 SNP markers were designed for two Illumina GoldenGate arrays and genotyped on more than 1000 samples for each species. In each case, oyster samples were obtained from wild and selected populations and from three‐generation families segregating for traits of interest in aquaculture. The rate of successfully genotyped polymorphic SNPs was about 60% for each species. Effects of SNP origin and quality on genotyping success (Illumina functionality Score) were analysed and compared with other model and nonmodel species. Furthermore, a simulation was made based on a subset of the C. gigas SNP array with a minor allele frequency of 0.3 and typical crosses used in shellfish hatcheries. This simulation indicated that at least 150 markers were needed to perform an accurate parental assignment. Such panels might provide valuable tools to improve our understanding of the connectivity between wild (and selected) populations and could contribute to future selective breeding programmes.     

Eleven polymorphic simple sequence repeat markers from expressed sequence tags of Pacific oyster Crassostrea gigas EST database

Oyster (Crassostrea gigas) is widely distributed in coastal areas of world. We developed and evaluated simple sequence repeat (SSR) markers from expressed sequence tags (ESTs) of C. gigas and to amplify EST-SSR in C. gigas. Characteristics of 11 EST-SSR loci were investigated using 45 C. gigas individuals. The number of alleles per locus ranged from two to thirteen. The observed heterozygosity (H _o) ranged from 0.0889 to 0.7333 and the expected heterozygosity (H _e) ranged from 0.0859 to 0.8981. Because of their high level of polymorphism, our 11 single-locus EST-SSR markers will be valuable tools for research on mating system, population genetics and systemic evolution of oyster in the future.     

Arsenic and Arsenic Species in Cultured Oyster (Crassostrea gigas and C. corteziensis) from Coastal Lagoons of the SE Gulf of California, Mexico

The aim of this study was to evaluate the bioavailability of arsenic (As) through cultured oyster Crassostrea gigas and Crassostrea corteziensis from four coastal lagoons (SE Gulf of California). Organisms were collected in two seasons (rainy and dry season), and they were analyzed for total arsenic and chemical speciation of this element. The concentrations of As in oyster soft tissue fluctuated between 5.44 and 9.56 μg/g for rainy season and 6.46 and 8.33 μg/g for dry season (dry weight) in C. gigas. In C. corteziensis, the As concentrations were <5 μg/g for both seasons (dry weight). Arsenic speciation indicated arsenobetaine as the major arseno-compound accounting for 43.2–76.3 % of total content of As. Lower contributions were obtained for non-extractable As (11.3–17.5 %) and other molecules such as arsenocholine and methyl-arsonate (<5 %). Inorganic arsenic was detectable in only two samples, at concentrations lower than <0.1 μg/g. These As data are the first generated for these mollusks in NW Mexico and indicate that C. gigas and C. corteziensis farmed in this area are safe for human consumption in terms of arseno-compounds.     

EST-SSR markers derived from Laminaria digitata and its transferable application in Saccharina japonica

The public availability of numerous expressed sequence tag (EST) enables EST-based SSR (simple sequence repeat) markers to be widely used for genetics and breeding studies. In the present study, EST-SSR markers were developed from ESTs of Laminaria digitata and were transferred to the non-congeneric species Saccharina japonica. Among the 2,668 non-redundant ESTs, 83 (3.1%) ESTs containing SSR were identified totally, with an average of one SSR per 13.6 kb. Analysis of SSR motifs revealed that the trinucleotide and tetranucleotide were major motifs, accounted for 44.58% and 16.87%, respectively. Based on the 83 ESTs containing SSR, we designed 45 pairs of primers in the flanking regions of the SSR, of which 13 pairs showed polymorphism in a wild S. japonica population, and the mean alleles per locus was 3.6 (ranging from 2 to 6). The observed (Ho) and expected (He) heterozygosities of these EST-SSRs were 0.234–0.632 and 0.260–0.635, respectively. All loci were in Hardy–Weinberg equilibrium in the wild population and no linkage disequilibrium was detected among loci. The obtained EST-SSR markers can facilitate and promote related research such as ecological investigation, genetic diversity assessment and breeding practice of S. japonica as well.