SNP Mining in Crassostrea gigas EST Data: Transferability to Four Other Crassostrea Species,Phylogenetic Inferences and Outlier SNPs under Selection

Oysters, with high levels of phenotypic plasticity and wide geographic distribution, are a challenging group for taxonomists and phylogenetics. Our study is intended to generate new EST-SNP markers and to evaluate their potential for cross-species utilization in phylogenetic study of the genus Crassostrea. In the study, 57 novel SNPs were developed from an EST database of C. gigas by the HRM (high-resolution melting) method. Transferability of 377 SNPs developed for C. gigas was examined on four other Crassostrea species: C. sikamea, C. angulata, C. hongkongensis and C. ariakensis. Among the 377 primer pairs tested, 311 (82.5%) primers showed amplification in C. sikamea, 353 (93.6%) in C. angulata, 254 (67.4%) in C. hongkongensis and 253 (67.1%) in C. ariakensis. A total of 214 SNPs were found to be transferable to all four species. Phylogenetic analyses showed that C. hongkongensis was a sister species of C. ariakensis and that this clade was sister to the clade containing C. sikamea, C. angulata and C. gigas. Within this clade, C. gigas and C. angulata had the closest relationship, with C. sikamea being the sister group. In addition, we detected eight SNPs as potentially being under selection by two outlier tests (fdist and hierarchical methods). The SNPs studied here should be useful for genetic diversity, comparative mapping and phylogenetic studies across species in Crassostrea and the candidate outlier SNPs are worth exploring in more detail regarding association genetics and functional studies.     

Interspecific hybridization in oysters: Restriction Enzyme Digestion Chromosome Banding confirms Crassostrea angulata × Crassostrea gigas F1 hybrids

The taxonomic status of the two commercially important cupped oysters, Crassostrea angulata, the Portuguese oyster (Lamarck, 1819) and Crassostrea gigas, the Japanese oyster (Thunberg, 1793) has long been in question. The recent observation of the hybridization between C. gigas and C. angulata and the production of fertile F1s led us to search for cytogenetic evidence of both parental genomes in the interspecific hybrids. The cytogenetic characterization of the hybrids was performed by the use of restriction endonuclease treatments. This technique has recently shown the potential for individual chromosome identification by banding in oysters. Chromosomes of C. gigas, C. angulata and their hybrids were treated with two different restriction enzymes (ApaI and HaeIII), stained with Giemsa, and examined for banding patterns. These chromosome markers allowed the parental haploid sets to be identified in the hybrids. The analysis of the banded karyotypes of the interspecific hybrids showed that for each chromosome pair, one of the homologues presented a banding pattern consistent with that of C. gigas and the other homologue presented a banding pattern consistent with that of C. angulata. These cytogenetic results substantiate the reported interspecific hybridization between C. gigas and C. angulata. In view of these results and taking into account the present expansion of C. gigas aquaculture in southern Europe, the question of the need for preservation of pure C. angulata stocks should be raised as only a few populations remain in the south of Spain and Portugal. Recently, changes in the genetic composition of populations in southern Portugal have indeed been observed, showing that human activities have created contact zones between the two taxa while no natural sympatric zones exist in Europe.     

Individual relationship between aneuploidy of gill cells and growth rate in the cupped oysters Crassostrea angulata, C. gigas and their reciprocal hybrids

The Portuguese oyster, Crassostrea angulata, is taxonomically close to the Pacific oyster, C. gigas, but there are clear genetic and phenotypic differences between these taxa. Among those differences, the faster growth of C. gigas compared with C. angulata has often been observed in the field. Crosses between C. angulata and C. gigas were performed to investigate the relationship between growth variation and somatic aneuploidy at the individual level in the two taxa and their reciprocal hybrids. The different progenies were reared in Ria Formosa (Portugal) under standard farming conditions. Growth rate and survival were significantly higher in C. gigas than in C. angulata, and the hybrids showed intermediate performances. Significant differences were also observed in the proportion of aneuploid cells (PAC) and of missing chromosomes (PMC) between the two taxa, C. angulata showing the highest values. Intermediate values of PAC and PMC were observed in the hybrids, supporting additive genetic bases of these parameters. Our results also confirm the negative correlation between somatic aneuploidy and growth rate at the individual level, as previously reported in C. gigas.     

Intraspecific Variation in Mitogenomes of Five <Emphasis Type="Italic">Crassostrea</Emphasis> Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China’s coast, which are shaped by unknown mechanisms in a north–south divide.     

Mitochondrial and nuclear DNA phylogeography of Crassostrea angulata, the Portuguese oyster endangered in Europe

The respective status of the Portuguese oyster, Crassostrea angulata, and the Pacific oyster, Crassostrea gigas, has long been a matter of controversy. Morphological and physiological similarities, homogeneity of allozyme allelic frequencies between populations of the two taxa and the demonstration of hybridization lead most authors to suggest that they should be regrouped within the same species. The risk of introgression and the present expansion of C. gigas aquaculture in Europe raises the question of the need for preservation of C. angulata in Europe, as only a few populations remain. We studied European and Asian populations of C. gigas and C. angulata using microsatellite and mitochondrial DNA markers to estimate their genetic diversity and differentiation. The analysis of genetic distances and the distribution of allelic and haplotype frequencies revealed significant genetic differences between taxa, showing two clusters: (1) C. gigas French and Japanese populations and (2) C. angulata Portuguese and Taiwanese populations. The Asian origin of the Crassostrea angulata taxa is therefore confirmed. Unlike previous studies based on allozymes, significant nuclear genome differences were noted between C. angulata and C. gigas. Despite the presumed history of the introduction of C. angulata into Southern Europe, these populations did not show any significant reduction of variability compared to Taiwanese populations. Any conservation plans for European C. angulata populations should take its non-native origin into account. They represent a valuable genetic resources for European breeding program.     

Environmental variability at the phosphoglucose isomerase locus in the genus Chlamys

Two geographically and reproductively isolated populations each of Chlamys opercularis and of Chlamys varia are compared with respect to polymorphism at the phosphoglucose isomerase (PGI) locus. One population of C. varia is further sub-divided into sub-littoral and inter-tidal samples. From the results it is inferred that PGI is a dimeric molecule, controlled by five alleles in C. varia and by at least four alleles in C. opercularis. acting at a single locus. The alleles recognised for C. opercularis have identical electrophoretic mobility as four of those seen in C. varia. The two populations of C. varia are remarkably similar, as were their habitats. However, the sub-littoral sample does show a lesser degree of homozygosity than the inter-tidal sample of the same population. The two populations of C. opercularis differ significantly and represent two separate and identifiable ‘races’. C. varia exhibits a much greater degree of genetic variability than the deep sub-littoral species, C. opercularis indicating that a more stable environment regulates a more monomorphic species.     

Characterization of 79 microsatellite DNA markers in the Pacific oyster Crassostrea gigas

We characterized 79 microsatellite DNA markers, which were obtained from genomic libraries enriched for CA, GA, ATG and TAGA motif repeats, in the Pacific oyster Crassostrea gigas. For eight F₁ grandparents or great‐grandparents of mapping families, the average heterozygosity, 0.705, and average number of alleles per locus, 5.7, did not vary among motif‐repeat or motif‐complexity categories. Non‐amplifying polymerase chain reaction null alleles, which were confirmed by segregation in the mapping families, were detected at 41 (51.9%) of the 79 loci. Cross‐species amplifications from C. angulata, C. sikamea, C. ariakensis and C. virginica showed a precipitous decline with distance from the focal species C. gigas.     

Iridoid glucosides in Griselinia, Aralidium and Toricellia

From three of five investigated species of Griselinia a new iridoid glucoside, griselinoside, was isolated. It was found to be present also in foliage of Aralidium pinnatifidum and Toricellia angulata, accompanied in the former by aralidioside another novel iridoid glucoside. The structures and absolute configurations of the two iridoids were elucidated by NMR spectroscopy and chemical conversions. From G. littoralis and T. angulata the glucosides magnolioside and syringoside respectively were isolated. ¹³C NMR spectra are given for thirteen iridoid derivatives.     

Complete mitochondrial DNA sequence of Crassostrea nippona: comparative and phylogenomic studies on seven commercial Crassostrea species

The complete mitochondrial genome of Crassostrea nippona was determined and compared with six other Crassostrea mitogenomes from GenBank in an attempt to shed light on the evolutionary relatedness within Crassostrea. The total length of the mitogenome was 20,030 bp for C. nippona, which was the largest among seven Crassostrea mitogenomes. Among six Asian oysters, the gene order of mitochondrial DNA was identical except for C. nippona with a transposition of trnG. While the American oyster C. virginica and Asian oysters showed broad differences in gene order with relocation of most tRNA genes and indels of duplicated tRNAs and rrnS, indicating the relatively distant relationships between the American oyster and Asian oysters. Different from other six Crassostrea oysters, C. nippona had two repeats of 66 bp in non-coding regions. Pairwise divergence among the seven Crassostrea oysters based on DNA sequences of 12 protein-coding genes ranged from 3.1 to 44.4% (Kimura two-parameter distance, K2P). The close relationship between C. nippona and C. hongkongensis was revealed by K2P of 18.9%. Phylogenetic analyses robustly revealed Crassostrea monophyly, with C. virginica at the basal position. The results of phylogenetic analyses strongly supported C. gigas and C. angulata had the closest relationship, with C. sikamea being the sister taxon. These findings presented here provide a better insight into the relationships within Crassostrea and will be useful for further evolution studies of oysters.     

Polymorphic microsatellite loci for two Atlantic oyster species: Crassostrea rhizophorae and C. gasar

Using a CA/CAA enriched library screening procedure, we isolated and characterised a total of seventeen polymorphic microsatellite loci for two species of Crassostrea with recognised economic importance. Eleven microsatellite loci were developed for C. rhizophorae, a Western Atlantic species for which no microsatellites were previously known. Another six loci were developed for C. gasar, a species that occurs on both sides of the South Atlantic, adding to the ten loci previously described for the species. The levels of polymorphism were estimated using 24 C. rhizophorae from Southeast Brazil (São Paulo) and 23 C. gasar individuals from North Brazil (Maranhão). The number of alleles per polymorphic locus varied from 3 to 27, and the observed and expected heterozygosities ranged between 0.174 and 0.958 and between 0.237 and 0.972 in C. rhizophorae and C. gasar, respectively. No linkage disequilibrium was found between any locus pair, and four of them exhibited deviations from Hardy–Weinberg expectations. Of the 17 loci developed, 8 cross-amplified in C. gigas and 13 in C. virginica. These markers are useful for evolution and population genetics studies of Crassostrea species and may provide fundamental data for the future cultivation of native oysters in Western Atlantic.     

Recently established <Emphasis Type="Italic">Crassostrea</Emphasis>-reefs versus native <Emphasis Type="Italic">Mytilus</Emphasis>-beds: differences in ecosystem engineering affects the macrofaunal communities (Wadden Sea of Lower Saxony,southern German Bight)

Since 1998 the non-indigenous Pacific oyster Crassostrea gigas (Thunberg 1793) has been invading the Wadden Sea of Lower Saxony, southern German Bight. C. gigas settles predominantly on intertidal Mytilus-beds (M. edulis L.) and subsequently create rigid reef-like structures. Both bivalve species are ecosystem engineers in sedimentary tidal flats. They provide hard substrate for sessile species, mobile organisms find refuge within the habitat matrix of dense suspension feeders, and biodeposits enrich the sediments with organic matter. The transformation of Mytilus-beds into Crassostrea-reefs gives rise to the question whether the invader may affect the native community. We investigated two parts of a changing bivalve bed in the backbarrier area of the island of Juist in March 2005. One part was still dominated by M. edulis whereas the other part was already densely colonized by C. gigas. Crassostrea-reefs compensate for the conceivable loss of Mytilus-beds in the intertidal of the Wadden Sea by replacing the ecological function of M. edulis. There was no indication of a suppression of indigenous species. This even applied to M. edulis, which persisted at the site invaded by C. gigas. The associated macrofaunal community showed increased species richness, abundance, biomass, and diversity in the Crassostrea-reef. The latter particularly favored sessile species like anthozoans, hydrozoans, and barnacles. Higher abundance and biomass for vagile epizoic species like the shore crab Carcinus maenas and the periwinkle Littorina littorea also occurred among oysters. Abundance of deposit feeding oligochaetes was enhanced by oysters as well. More opportunistic, facultative filter-feeding polychaetes occurred in the Crassostrea-reef.     

Crassostrea angulata Bindin Gene and the Divergence of Fucose-Binding Lectin Repeats Among Three Species of Crassostrea

Bindin is a major protein for species-specific recognition between sperm and congenetic egg in many free-spawning marine invertebrates. We cloned a novel bindin gene from the oyster Crassostrea angulata by 3′ and 5′ rapid amplification of cDNA ends. The full-length bindin cDNA was 1,049 bp with a 771-bp open reading frame encoding 257 amino acids. The deduced amino acid sequence contained a putative signal peptide of 24 amino acids. The length of the bindin genomic DNA was 8,508 bp containing four exons and three introns. Three haplotypes of F-lectin repeat were detected from seven sequences of F-lectin repeat of six male oysters. Both neighbor-joining and minimum-evolution phylogenetic trees show that haplotype an1 was close to Crassostrea gigas while an2 and an3 were close to Crassostrea sikamea. Intron-4 in the middle of F-lectin repeat is highly variable in both size and sequence. We classified intron-4 into three types according to their size and the F-lectin repeat they were located in. Intron-4 may play an important role in recombination. We compared the number of nonsynonymous substitutions (Dn) and synonymous substitutions (Ds) per nucleotide site among 19 F-lectin haplotypes of the three species. Dn/Ds ratios suggested that positive selection occurred between C. gigas and C. sikamea and between C. gigas and C. angulata. Nine positive selected positions (p > 90%) are identified among 19 haplotypes of three species. They are located on the F-lectin binding face around the three recognition motif residues. We assume that these nine clustered amino acids are related with species-specific recognition.     

Enzyme polymorphisms in the European oyster, Ostrea edulis L.

Extracts of the digestive diverticula of more than 300 individuals from four geographically separated populations of the European oyster, Ostrea edulis L. have been examined by electrophoresis for esterase and phosphoglucose isomerase polymorphisms. Three regions of esterase variability have been detected in the electrophoretograms. It is inferred that the variants in two of these regions are governed by two co-dominant alleles at each of two loci. Two phenotypes of phosphoglucose isomerase, which is inferred to be a dimeric enzyme encoded by two alleles at a single locus, have been observed. Allele frequencies and phenotype distributions are compared in the four populations. It is concluded that the populations differ genetically at the loci investigated, the magnitude of the differences being generally greatest between different 'physiological races'.     

Specific Detection of Pacific Oyster (Crassostrea gigas) Larvae in Plankton Samples Using Nested Polymerase Chain Reaction

Management of sustainable Pacific oyster fisheries would be assisted by an early, rapid, and accurate means of detecting their planktonic larvae. Reported here is an approach, based on polymerase chain reaction (PCR), for the detection of Pacific oyster larvae in plankton samples. Species-specific primers were designed by comparing partial mitochondrial cytochrome oxidase subunit I (COI) sequences from Crassostrea gigas, with other members of the family Ostreidae including those of Crassostrea angulata. Assay specificity was empirically validated through screening DNA samples obtained from several species of oysters. The assay was specific as only C. gigas samples returned PCR-positive results. A nested PCR approach could consistently detect 5 or more D-hinge-stage larvae spiked into a background of about 146 mg of plankton. The assay does not require prior sorting of larvae. We conclude that the assay could be used to screen environmental and ballast water samples, although further specificity testing against local bivalve species is recommended in new locations.     

Temperature,energy metabolism,and adaptive divergence in two oyster subspecies

Comparisons of related species that have diverse spatial distributions provide an efficient way to investigate adaptive evolution in face of increasing global warming. The oyster subjected to high environmental selections is a model species as sessile marine invertebrate. This study aimed to detect the adaptive divergence of energy metabolism in two oyster subspecies from the genus Crassostrea—C. gigas gigas and C. gigas angulata—which are broadly distributed along the northern and southern coasts of China, respectively. We examined the effects of acute thermal stress on energy metabolism in two oyster subspecies after being common gardened for one generation in identical conditions. Thermal responses were assessed by incorporating physiological, molecular, and genomic approaches. Southern oysters exhibited higher fluctuations in metabolic rate, activities of key energetic enzymes, and levels of thermally induced gene expression than northern oysters. For genes involved in energy metabolism, the former displayed higher basal levels of gene expression and a more pronounced downregulation of thermally induced expression, while the later exhibited lower basal levels and a less pronounced downregulation of gene expression. Contrary expression pattern was observed in oxidative stress gene. Besides, energy metabolic tradeoffs were detected in both subspecies. Furthermore, the genetic divergence of a nonsynonymous SNP (SOD‐132) and five synonymous SNPs in other genes was identified and validated in these two subspecies, which possibly affects downstream functions and explains the aforementioned phenotypic variations. Our study demonstrates that differentiations in energy metabolism underlie the plasticity of adaptive divergence in two oyster subspecies and suggest C. gigas angulata with moderate phenotypic plasticity has higher adaptive potential to cope with exacerbated global warming.     

A preliminary integrated genetic map distinguishes every chromosome pair and locates essential genes related to abiotic adaptation of <Emphasis Type="Italic">Crassostrea angulata</Emphasis>/<Emphasis Type="Italic">gigas</Emphasis>

Background

The re-sequencing of C. angulata has revealed many polymorphisms in candidate genes related to adaptation to abiotic stress that are not present in C. gigas; these genes, therefore, are probably related to the ability of this oyster to retain high concentrations of toxic heavy metals. There is, in addition, an unresolved controversy as to whether or not C. angulata and C. gigas are the same species or subspecies. Both oysters have 20 metacentric chromosomes of similar size that are morphologically indistinguishable. From a genomic perspective, as a result of the great variation and selection for heterozygotes in C. gigas, the assembly of its draft genome was difficult: it is fragmented in more than seven thousand scaffolds.

Results

In this work sixty BAC sequences of C. gigas downloaded from NCBI were assembled in BAC-contigs and assigned to BACs that were used as probes for mFISH in C. angulata and C. gigas. In addition, probes of H3, H4 histone, 18S and 5S rDNA genes were also used. Hence we obtained markers identifying 8 out the 10 chromosomes constituting the karyotype. Chromosomes 1 and 9 can be distinguished morphologically. The bioinformatic analysis carried out with the BAC-contigs annotated 88 genes. As a result, genes associated with abiotic adaptation, such as metallothioneins, have been positioned in the genome. The gene ontology analysis has also shown many molecular functions related to metal ion binding, a phenomenon associated with detoxification processes that are characteristic in oysters. Hence the provisional integrated map obtained in this study is a useful complementary tool for the study of oyster genomes.

Conclusions

In this study 8 out of 10 chromosome pairs of Crassostrea angulata/gigas were identified using BAC clones as probes. As a result all chromosomes can now be distinguished. Moreover, FISH showed that H3 and H4 co-localized in two pairs of chromosomes different that those previously escribed. 88 genes were annotated in the BAC-contigs most of them related with Molecular Functions of protein binding, related to the resistance of the species to abiotic stress. An integrated genetic map anchored to the genome has been obtained in which the BAC-contigs structure were not concordant with the gene structure of the C. gigas scaffolds displayed in the Genomicus database.

     

Phenotypes of Phosphoglucose isomerase,phosphoglucose mutase and general protein in some freshwater molluscs from Basrah,Iraq

Electrophoretic variants of phosphoglucose isomerase (EC.5.3.1.9) and phosphoglucose mutase (EC.2.7.5.1) have been studied in eight species of freshwater molluscs. Two phenotypes of phosphoglucose isomerase were observed in Melanopsis nodosa and one phenotype was observed in the rest of the species. One phenotype of phosphoglucose mutase was observed in all the species of molluscs studied. Phosphoglucose isomerase is inferred to be a dimer encoded at a single polymorphic locus in Melanoides nodosa. There are two alleles at this locus. Phosphoglucose mutase is inferred to be a monomer encoded at a single monomorphic locus in all species. The electrophoretic analysis revealed that phosphoglucose isomerase enzyme cannot be considered a good taxonomic criterion to differentiate the different members of the six families studied but, on the other hand, it is considered a good taxonomic criterion to differentiate Melanopsis nodosa and Theodoxus jordani. Phosphoglucose mutase is considered a good taxonomic criterion to differentiate the family Melanidae from the remaining five families studied. General protein can be considered a good taxonomic criterion to differentiate the family Corbicullidae from Melanidae, Viviparidae and Neritidae but, on the other hand, it seems to be a less useful taxonomic criterion to differentiate between the Viviparidae and Neritidae.     

EST‐SSR markers from the Pacific oyster,Crassostrea gigas

Ten polymorphic microsatellite repeat markers were identified from Crassostrea gigas, expressed sequence tags (EST) deposited in public sequence database. Number of alleles per locus ranged from three to 18, expected and observed heterozygosities ranged from 0.071 to 0.738 and from 0.306 to 0.913, respectively. Marker transferability was tested on other two Crassostrea species and polymorphic products were detected at nine loci. EST‐derived simple sequence repeats provide robust, informative and potentially transferable polymorphic markers suitable for population genetic, parentage, and mapping studies of C. gigas.     

Development of Novel Microsatellite DNA Markers from the Pacific Oyster Crassostrea gigas

We document the potential of novel microsatellites as a genetic tool in furthering our understanding of the Crassostrea gigas genetic structure. From the microsatellite-enriched libraries we constructed, 123 repeat regions that had sufficient sequence information to design polymerase chain reaction primer sets were isolated. From these, 9 primer pairs were screened in a C. gigas population of 67 individuals to evaluate the genetic variability. All but 1 of the 9 loci showed allelic variation (number of alleles, 2–20; observed heterozygosity, 0.119–0.925; unbiased expected heterozygosity, 0.139–0.914). Considerable discrepancy of genotypic proportions from the Hardy-Weinberg equilibrium was observed at 1 locus with an apparent heterozygote deficiency. Several loci were successfully amplified in 3 other related species with the appropriate allele size: 6 loci in C. sikamea, 4 loci in C. ariakensis, and 5 loci in C. nippona.     

板栗和锥栗同域居群的空间遗传结构

采用微卫星技术研究了板栗(Castanea mollissima)和锥栗(C. henryi)同域分布居群的空间遗传结构, 探讨了它们的遗传变异空间分布特征及其形成机制。结果表明,采用的7个微卫星在两物种中共扩增出173个等位基因,两物种都具有较高的遗传多样性且种间遗传分化水平较低(F_ST=0.051), 但空间自相关分析揭示板栗和锥栗在同域居群中具有不同的空间遗传结构,锥栗在100 m内有空间遗传结构,而板栗没有;同时,基于亲缘关系系数F_ij的Sp统计值也显示锥栗具有比板栗更强的空间遗传结构(板栗的Sp=0.002;锥栗的Sp=0.018)。居群遗传变异的空间结构是其传粉和种子散播及生境共同作用的结果,其中种子的近距离散播和居群密度可能是主要的因素。板栗和锥栗居群在小尺度上空间遗传结构的差异可能反映了它们种子的大小及扩散过程的差异。