Variable microsatellites in the Pacific oyster Crassostrea gigas and other cupped oyster species

Mapping Genetic Loci for Quantitative Traits of Golden Shell Color, Mineral Element Contents, and Growth-Related Traits in Pacific Oyster (Crassostrea gigas)

Golden shell color and mineral content are important economic traits of Pacific oyster (Crassostrea gigas). In this study, we mapped a series of quantitative trait loci (QTLs) that control zinc (Zn) and magnesium (Mg) content, shell color and growth performance to two sex-averaged linkage maps from the FAM-A and FAM-B families. In total, ten QTLs were identified in seven linkage groups (LGs) in the FAM-B family, and seven QTLs were identified in four linkage groups in the FAM-A family. Two QTLs affecting the trait of golden shell color were identified in LG8 of the FAM-A and LG10 of the FAM-B families, which could explain 20.2 and 10.5% of the phenotypic variations, respectively. Two QTLs for Zn content were identified that could contribute to 17.9 and 34.44% of the phenotypic variations in FAM-A. Six QTLs for Zn and Mg contents were identified in four LGs (LG1, LG2, LG5, and LG9) in FAM-B, which explained 13.5-26.7% of the phenotypic variations. In addition, seven QTLs related to oyster growth were recognized in both FAM-A and FAM-B families accounting for 14.6-36.7% of the phenotypic variations. All of the DNA markers in QTL regions were blasted and 14 genes associated with above traits were identified. The mRNA expression of these genes was determined by quantitative RT-PCR. These QTLs and candidate genes could be used as potential targets for marker-assisted selection in C. gigas breeding.

Rapid expansion of the invasive oyster Crassostrea gigas at its northern distribution limit in Europe: Naturally dispersed or introduced?

The Pacific oyster, Crassostrea gigas, was introduced to Europe for aquaculture purposes, and has had a rapid and unforeseen northward expansion in northern Europe. The recent dramatic increase in number of C. gigas populations along the species’ northern distribution limit has questioned the efficiency of Skagerrak as a dispersal barrier for transport and survival of larvae. We investigated the genetic connectivity and possible spreading patterns between Pacific oyster populations on the southern Norwegian coast (4 localities) and Swedish and Danish populations by means of DNA microsatellite analysis of adult oysters, and by simulating larvae drift. In the simulations we used a 3D oceanographic model to explore the influence of recent climate change (1990–2010) on development, survival, and successful spreading of Danish and Swedish Pacific oyster larvae to Norwegian coastal waters. The simulations indicated adequate temperature conditions for development, survival, and settlement of larvae across the Skagerrak in warm years since 2000. However, microsatellite genotyping revealed genetic differences between the Norwegian populations, and between the Norwegian populations and the Swedish and Danish populations, the latter two populations being more similar. This patchwork pattern of genetic dissimilarity among the Norwegian populations points towards multiple local introduction routes rather than the commonly assumed unidirectional entry of larvae drifted from Denmark and Sweden. Alternative origins of introduction and implications for management, such as forecasting and possible mitigation actions, are discussed.

Additive transcriptomic variation associated with reproductive traits suggest local adaptation in a recently settled population of the Pacific oyster, Crassostrea gigas

Background

Originating from Northeast Asia, the Pacific oyster Crassostrea gigas has been introduced into a large number of countries for aquaculture purpose. Following introduction, the Pacific oyster has turned into an invasive species in an increasing number of coastal areas, notably recently in Northern Europe.

Methods

To explore potential adaptation of reproductive traits in populations with different histories, we set up a common garden experiment based on the comparison of progenies from two populations of Pacific oyster sampled in France and Denmark and their hybrids. Sex ratio, condition index and microarray gene expression in gonads, were analyzed in each progeny (n = 60).

Results

A female-biased sex-ratio and a higher condition index were observed in the Danish progeny, possibly reflecting an evolutionary reproductive strategy to increase the potential success of natural recruitment in recently settled population. Using multifarious statistical approaches and accounting for sex differences we identified several transcripts differentially expressed between the Danish and French progenies, for which additive genetic basis is suspected (showing intermediate expression levels in hybrids, and therefore additivity). Candidate transcripts included mRNA coding for sperm quality and insulin metabolism, known to be implicated in coordinated control and success of reproduction.

Conclusions

Observed differences suggest that adaptation of invasive populations might have occurred during expansion acting on reproductive traits, and in particular on a female-biased sex-ratio, gamete quality and fertility.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1972-8) contains supplementary material, which is available to authorized users.

Microsatellite analysis of genetic diversity and population structure of freshwater mussel (Lamprotula leai)

Lamprotula leai is one of the most commercially important freshwater pearl mussels in China, but there is limited data on its genetic diversity and population structure. In the present study, 119 individuals from four major geographical populations were investigated using 15 microsatellite loci identified via cross-species amplification. A total of 114 alleles were detected, with an average of 7.6 alleles per locus (range: 2 to 21). Among the four stocks, those from Hung-tse Lake and Poyang Lake had the lowest (0.412) and highest (0.455) observed heterozygosity respectively. The polymorphism information content (PIC) ranged from 0.374 to 0.927 (mean: 0.907). AMOVA showed that 12.56% and 44.68% genetic variances were among populations and within individuals, respectively. Pairwise Fst ranged from 0.073 to 0.146, indicating medium genetic differentiation among the populations. In aggregate, our results suggest that inbreeding is a crucial factor accounting for deviations from Hardy-Weinberg equilibrium at 12 loci. Moreover, the genetic distance among four stocks ranged from 0.192 to 0.890. Poyang Lake and Hung-tse Lake were clustered together, joined with Dongting Lake and Anqing Lake. Given that specimens from Hung-tse Lake showed the highest average allele richness, expected heterozygosity and PIC, this location may be the source of the highest quality germplasm resources and the stock from this area may be the best for future breeding efforts.

Population genetics of the black scar oyster, Crassostrea iredalei: repercussion of anthropogenic interference

Mitochondrial cytochrome oxidase subunit I (COI) gene was utilized to assess the population genetics of the commercially important black scar oyster, Crassostrea iredalei among 11 populations throughout the west and east coasts Peninsular Malaysia and Sabah (Malaysian Borneo). Overall, populations of C. iredalei demonstrated low nucleotide diversity π (0.000-0.004) and low-to-high haplotype diversity h (0.000-0.795) levels. Genetic structuring was detected between the Peninsular Malaysia and Sabah populations as revealed by the FST analysis. However, the COI gene analyses showed minimal and non-significant (p > 0.05) population differentiation within the east and west coasts Peninsular Malaysia and Sabah regions. This was attributed to both high larval dispersal along the east and west coasts and human-driven spat translocation between the two coastlines due to C. iredalei cultivation practices. Phylogeographic relationships inferences were also conducted to further support these hypotheses. The neutrality and mismatch distribution analyses suggested that C. iredalei had experienced a/several bottleneck event(s), followed by population expansion. The molecular information obtained from this study could be incorporated in a pragmatic aquaculture management strategy of wild broodstock and the hatchery lines of C. iredalei in Malaysia.

Population genomics shed light on the demographic and adaptive histories of European invasion in the Pacific oyster, Crassostrea gigas

Crassostrea gigas originated from the Pacific coast of Asia, but was introduced into several European countries in the early 1970s. Natural populations have now spread across the length of the western seaboard of Europe. To elucidate the demographic and selective processes at play during this rapid expansion, genome-scan analysis was performed on different populations. High diversities and low differentiation were observed overall, but significant genetic differentiation was found among newly established populations and between the newly established northern group and a nearly panmictic group composed of southern European populations and a population from Japan. Loss of genetic diversity was also seen in the north, likely caused by founder events during colonization. The few strongly supported outlier loci revealed a genetic structure uncorrelated with the north/south differentiation, but grouping two samples from the Danish fjords (northern group) and one from the Dutch Scheldt estuary (southern group) with the one from Japan. These findings might reflect the following: (i) parallel adaptation to similar environmental pressures (fjord-like environment) within each of the two groups or (ii) a footprint of a secondary introduction of an alternative genomic background maintained by multifarious isolation factors. Our results call for a closer examination of adaptive genetic structure in the area of origin.

Genetic evidence for the uncoupling of local aquaculture activities and a population of an invasive species--a case study of Pacific oysters (Crassostrea gigas)

Human-mediated introduction of nonnative species into coastal areas via aquaculture is one of the main pathways that can lead to biological invasions. To develop strategies to counteract invasions, it is critical to determine whether populations establishing in the wild are self-sustaining or based on repeated introductions. Invasions by the Pacific oyster (Crassostrea gigas) have been associated with the growing oyster aquaculture industry worldwide. In this study, temporal genetic variability of farmed and wild oysters from the largest enclosed bay in Ireland was assessed to reconstruct the recent biological history of the feral populations using 7 anonymous microsatellites and 7 microsatellites linked to expressed sequence tags (ESTs). There was no evidence of EST-linked markers showing footprints of selection. Allelic richness was higher in feral than in aquaculture samples (P = 0.003, paired t-test). Significant deviations from Hardy-Weinberg equilibrium due to heterozygote deficiencies were detected for almost all loci and samples, most likely explained by the presence of null alleles. Relatively high genetic differentiation was found between aquaculture and feral oysters (largest pairwise multilocus F(ST) 0.074, P < 0.01) and between year classes of oysters from aquaculture (largest pairwise multilocus F(ST) 0.073, P < 0.01), which was also confirmed by the strong separation of aquaculture and wild samples using Bayesian clustering approaches. A 10-fold higher effective population size (N(e)) and a high number of private alleles in wild oysters suggest an established self-sustaining feral population. The wild oyster population studied appears demographically independent from the current aquaculture activities in the estuary and alternative scenarios of introduction pathways are discussed.

Translocation of wild populations: conservation implications for the genetic diversity of the black-lipped pearl oyster Pinctada margaritifera

Translocation has been widely studied as a tool for conservation management to restore or enhance degraded populations. On the contrary, few studies have been conducted on translocation for commercial purposes. In this study, we evaluate the genetic consequences of translocation of wild individuals of Pinctada margaritifera on farmed and adjacent wild populations. We tested the hypotheses that translocations would induce high genetic heterogeneity in farmed populations and this heterogeneity would then leak into the adjacent wild populations. In fact, farmed samples exhibit high levels of heterogeneity and low pairwise relatedness compared to wild populations, highlighting the pooling of genetically divergent populations into farms. We also demonstrate that this heterogeneity is transmitted to adjacent wild populations as a result of interbreeding. Adjacent wild populations tend to have higher genetic diversity values and greater pairwise relatedness coefficient with farmed populations than wild populations. Overall, pearl culture in French Polynesia promotes the mixing of unrelated individuals in farmed locations and reduces genetic divergence among geographically distant populations as well as among farmed and wild populations of a same lagoon. We also studied for the first time a farmed population originating from spat collected in a lagoon where release of hatchery-produced larvae occurred 10 years ago and we were able to identify four distinct genetic groups. These groups contribute highly to reproduction and caused considerable genetic drift in the lagoon, suggesting that hatchery-produced larvae are neither sustainable method for pearl culture nor for conserving the diversity of P. margaritifera in French Polynesia.

Genetic mapping and QTL analysis of growth-related traits in the Pacific oyster

The Pacific oyster (Crassostrea gigas) is one of the most important oysters cultured worldwide. To analyze the oyster genome and dissect growth-related traits, we constructed a sex-averaged linkage map by combining 64 genomic simple sequence repeats, 42 expressed sequence tag-derived SSRs, and 320 amplified fragment length polymorphism markers in an F(1) full-sib family. A total of 426 markers were assigned to 11 linkage groups, spanning 558.2 cM with an average interval of 1.3 cM and 94.7% of genome coverage. Segregation distortion was significant for 18.8% of the markers (P < 0.05), and distorted markers tended to occur on some genetic regions or linkage groups. Most growth-related quantitative traits were highly significantly (P < 0.01) correlated, and principal component analysis obtained four principal components. Quantitative trait locus (QTL) analysis identified three significant QTLs for two principal components, which explained 0.6-13.9% of the phenotypic variation. One QTL for sex was detected on linkage group 6, and the inheritabilities of sex for parental alleles and maternal alleles on that locus C15 are 39.8% and 0.01%, respectively. The constructed linkage map and determined QTLs can provide a tool for further genetic analysis of the traits and be potential for marker-assisted selection in C. gigas breeding.

Quantitative trait locus analysis of stage-specific inbreeding depression in the Pacific oyster Crassostrea gigas

Inbreeding depression and genetic load have been widely observed, but their genetic basis and effects on fitness during the life cycle remain poorly understood, especially for marine animals with high fecundity and high, early mortality (type-III survivorship). A high load of recessive mutations was previously inferred for the Pacific oyster Crassostrea gigas, from massive distortions of zygotic, marker segregation ratios in F(2) families. However, the number, genomic location, and stage-specific onset of mutations affecting viability have not been thoroughly investigated. Here, we again report massive distortions of microsatellite-marker segregation ratios in two F(2) hybrid families, but we now locate the causative deleterious mutations, using a quantitative trait locus (QTL) interval-mapping model, and we characterize their mode of gene action. We find 14-15 viability QTL (vQTL) in the two families. Genotypic frequencies at vQTL generally suggest selection against recessive or partially recessive alleles, supporting the dominance theory of inbreeding depression. No epistasis was detected among vQTL, so unlinked vQTL presumably have independent effects on survival. For the first time, we track segregation ratios of vQTL-linked markers through the life cycle, to determine their stage-specific expression. Almost all vQTL are absent in the earliest life stages examined, confirming zygotic viability selection; vQTL are predominantly expressed before the juvenile stage (90%), mostly at metamorphosis (50%). We estimate that, altogether, selection on vQTL caused 96% mortality in these families, accounting for nearly all of the actual mortality. Thus, genetic load causes substantial mortality in inbred Pacific oysters, particularly during metamorphosis, a critical developmental transition warranting further investigation.

Variance in the reproductive success of flat oyster Ostrea edulis L. assessed by parentage analyses in natural and experimental conditions

In order to document further the phenomena of variance in reproductive success in natural populations of the European flat oyster Ostrea edulis, two complementary studies based on natural and experimental populations were conducted. The first part of this work was focused on paternity analyses using a set of four microsatellite markers for larvae collected from 13 brooding females sampled in Quiberon Bay (Brittany, France). The number of individuals contributing as the male parent to each progeny assay was highly variable, ranging from 2 to more than 40. Moreover, paternal contributions showed a much skewed distribution, with some males contributing to 50-100% of the progeny assay. The second part of this work consisted of the analysis of six successive cohorts experimentally produced from an acclimated broodstock (62 wild oysters sampled in the Quiberon Bay). Allelic richness was significantly higher in the adult population than in the temporal cohorts collected. Genetic differentiation (F(st) estimates) was computed for each pair of samples and all significant values ranged from 0.7 to 11.9%. A limited effective number of breeders (generally below 25) was estimated in the six temporal cohorts. The study gives first indications of the high variance in reproductive success as well as a reduced effective size, not only under experimental conditions but also in the wild. Surprisingly, the pool of the successive cohorts, based on the low number of loci used, appeared to depict a random and representative set of alleles of the progenitor population, indicating that the detection of patterns of temporal genetic differentiation at a local scale most likely depends on the sampling window.

QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (Crassostrea gigas)

Summer mortality is a phenomenon severely affecting the aquaculture production of the Pacific oyster (Crassostrea gigas). Although its causal factors are complex, resistance to mortality has been described as a highly heritable trait, and several pathogens including the virus Ostreid Herpes virus type 1 (OsHV-1) have been associated with this phenomenon. A QTL analysis for survival of summer mortality and OsHV-1 load, estimated using real-time PCR, was performed using five F(2) full-sib families resulting from a divergent selection experiment for resistance to summer mortality. A consensus linkage map was built using 29 SNPs and 51 microsatellite markers. Five significant QTL were identified and assigned to linkage groups V, VI, VII and IX. Analysis of single full-sib families revealed differential QTL segregation between families. QTL for the two-recorded traits presented very similar locations, highlighting the interest of further study of their respective genetic controls. These QTL show substantial genetic variation in resistance to summer mortality, and present new opportunities for selection for resistance to OsHV-1.

Genetic variation of wild and hatchery populations of the Pacific oyster Crassostrea gigas assessed by microsatellite markers

Microsatellite DNA technique was used to detect the genetic variation between five hatchery populations of the Pacific oyster from China and two wild populations from Japan. Seven microsatellite loci screened in this study showed high polymorphism in both hatchery and wild populations, as observed in an average number of allele per locus (19.1-29.9) and average expected heterozygosity (0.916-0.958). No significant difference in average allelic richness or expected heterozygosity was observed between Chinese hatchery populations and Japanese wild populations. Pairwise FST values and heterogeneity tests of allele frequencies showed significant genetic differentiation between all populations. According to the neighbor-joining tree constructed on the basis of the DC distance, the seven populations fell into three groups showing a clear division between hatchery and wild populations, and between the northern and southern hatchery populations. Assignment tests correctly assigned high percentages (97%-100%) of individuals to their original populations and demonstrated the feasibility of microsatellite analysis for discrimination between populations. The information obtained in this study is useful for designing suitable management guidelines and selective breeding programs for the Pacific oyster in China.

Development and characterization of EST-SSR markers in the eastern oyster Crassostrea virginica

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.

Genetic differentiation in relation to marine landscape in a broadcast-spawning bivalve mollusc (Placopecten magellanicus)

Marine bivalves are sessile or sedentary as adults but have planktonic larvae which can potentially disperse over large distances. Consequently larval transport is expected to play a prominent role in facilitating gene flow and determining population structure. The sea scallop (Placopecten magellanicus) is a dioecious species with high fecundity, broadcast spawning and a c. 30-day planktonic larval stage, yet it forms discrete populations or 'beds' which have significantly different dynamics and characteristics. We analysed variation at six microsatellite loci in 12 locations throughout the geographic range of the species from Newfoundland, Canada, to New Jersey, USA. Significant differentiation was present and the maximum pairwise theta value, between one of the Newfoundland samples in the north and a sample from the southern portion of the range, was high at 0.061. Other proximate pairs of samples had no detectable genetic differentiation. Mantel tests indicated a significant isolation by distance, but only when one of the populations was excluded. A landscape genetic approach was used to detect areas of low gene flow using a joint analysis of spatial and genetic information. The two major putative barriers inferred by Monmonier's algorithm were then used to define regions for an analysis of molecular variance (amova). That analysis showed a significant but low percentage (1.2%) of the variation to be partitioned among regions, negligible variation among populations within regions, and the majority of the variance distributed between individuals within populations. Prominent currents were concordant with the demarcation of the regions, while a novel approach of using particle tracking software to mimic scallop larval dispersal was employed to interpret within-region genetic patterns.

A complementary method for production of tetraploid Crassostrea gigas using crosses between diploids and tetraploids with cytochalasin b treatments

We present a new method to produce tetraploid Crassostrea gigas by cytochalasin B inhibition of polar body 2 expulsion in diploid females crossed with tetraploid males. This offers a means of direct introgression of genetic characters from selected diploid to tetraploid lines, avoiding a triploid step. Offspring larval ploidy shifted over time and depended on size, with tetraploids more frequent among the smaller larvae and triploids among the large. Viable tetraploids were found at 4 and 6 months, indicating the technique was successful. The possibility that gynogenesis occurred was tested by microsatellite analysis to confirm the presence of paternally inherited alleles. These were present in all animals of the 2n x 4n + CB (female first) cross. However, a 4n x 2n + CB cross produced triploids, including some gynogens. Our method illustrates for the first time that diploid C. gigas eggs, if selected for large size, can give viable tetraploid offspring.

Linkage maps of microsatellite DNA markers for the Pacific oyster Crassostrea gigas

We constructed male and female consensus linkage maps for the Pacific oyster Crassostrea gigas, using a total of 102 microsatellite DNA markers typed in 11-day-old larvae from three families. We identified 11 and 12 linkage groups in the male and female consensus maps, respectively. Alignment of these separate maps, however, suggests 10 linkage groups, which agrees with the haploid chromosome number. The male linkage map comprises 88 loci and spans 616.1 cM, while the female map comprises 86 loci and spans 770.5 cM. The male and the female maps share 74 loci; 2 markers remain unlinked. The estimated coverages for the consensus linkage maps are 79% for the male and 70-75% for the female, on the basis of two estimates of genome length. Ninety-five percent of the genome is expected to lie within 16 and 21 cM of markers on the male and female maps, respectively, while 95% of simulated minimum distances to the male and female maps are within 10.1 and 13.6 cM, respectively. Females have significantly more recombination than males, across 118 pairs of linked markers in common to the parents of the three families. Significant differences in recombination and orders of markers are also evident among same-sex parents of different families as well as sibling parents of opposite sex. These observations suggest that polymorphism for chromosomal rearrangements may exist in natural populations, which could have profound implications for interpreting the evolutionary genetics of the oyster. These are the first linkage maps for a bivalve mollusc that use microsatellite DNA markers, which should enable them to be transferred to other families and to be useful for further genetic analyses such as QTL mapping.

High genetic load in the Pacific oyster Crassostrea gigas

The causes of inbreeding depression and the converse phenomenon of heterosis or hybrid vigor remain poorly understood despite their scientific and agricultural importance. In bivalve molluscs, related phenomena, marker-associated heterosis and distortion of marker segregation ratios, have been widely reported over the past 25 years. A large load of deleterious recessive mutations could explain both phenomena, according to the dominance hypothesis of heterosis. Using inbred lines derived from a natural population of Pacific oysters and classical crossbreeding experiments, we compare the segregation ratios of microsatellite DNA markers at 6 hr and 2-3 months postfertilization in F(2) or F(3) hybrid families. We find evidence for strong and widespread selection against identical-by-descent marker homozygotes. The marker segregation data, when fit to models of selection against linked deleterious recessive mutations and extrapolated to the whole genome, suggest that the wild founders of inbred lines carried a minimum of 8-14 highly deleterious recessive mutations. This evidence for a high genetic load strongly supports the dominance theory of heterosis and inbreeding depression and establishes the oyster as an animal model for understanding the genetic and physiological causes of these economically important phenomena.