Morphological and Molecular Analysis Identified a Subspecies of Crassostrea ariakensis (Fujita, 1913) along the Coast of Asia.

Crassostrea ariakensis (Fujita, 1913) is one of the most important economic and ecological oysters that is naturally distributed along the coast of Asia, separated by the Yangtze River estuary. They are usually compared as different populations, while there is no consensus on whether C. ariakensis in northern and southern areas should be considered as two species or subspecies. Here, we analyzed morphological characteristics, COI, 16s rRNA, mitogenome sequences, and species delimitation analysis (ASAP and PTP) to resolve the intraspecific taxonomic status of the C. ariakensis. Phylogenetic and ASAP analysis highlight that C. ariakensis was divided into N-type and S-type. PTP was unable to differentiate between the two types of C. ariakensis. The divergence time of N-type and S-type C. ariakinsis is estimated to be 1.6 Mya, using the relaxed uncorrelated lognormal clock method. Additionally, significant morphological differences exist between the two groups in terms of the adductor muscle scar color. Despite these differences, the COI (0.6%) and 16S rRNA (0.6%) genetic distance differences between N-type and S-type C. ariakensis has not yet reached the interspecific level. These results suggest that N-type and S-type C. ariakensis should be treated as different subspecies and renamed as C. ariakensis ariakensis subsp. nov and C. ariakensis meridioyangtzensis subsp. nov.

Genetic divergence accessed with microsatellite markers reflects the time of Crassostrea gigas genetic breeding in Brazil.

The Pacific Oyster was introduced on Santa Catarina Island in 1987, experiencing processes of selection and genetic breeding since then. Such procedures may have led to the establishment of specific strains, given the saltier and warmer conditions of the Atlantic Ocean. This study employed microsatellite markers to compare allelic patterns of oysters cultivated in Santa Catarina, the USA, and Asia. Specific allelic patterns were revealed in the Santa Catarina samples, reflecting the time of selection/breeding of the oyster in this region. This result supports the effectiveness of the selection/breeding procedures and the demand for protection of this commercially important genetic resource.

Transient Receptor Potential (TRP) Channels in the Pacific Oyster (Crassostrea gigas): Genome-Wide Identification and Expression Profiling after Heat Stress between C. gigas and C. angulata.

Transmembrane proteins are involved in an array of stress responses, particularly in thermo-sensation and thermo-regulation. In this study, we performed a genome-wide identification and characterization of the Transient Receptor Potential (TRP) genes in the Pacific oyster (Crassostrea gigas) and investigated their expression profiles after heat stress to identify critical TRPs potentially associated with thermal regulation. A total of 66 TRP genes were identified in the C. gigas, which showed significant gene expansion and tandem duplication. Meta-analysis of the available RNA-Seq data generated from samples after acute heat stress revealed a set of heat-inducible TRPs. Further examination of their expression profiles under chronic heat stress, and comparison between C. gigas and C. angulata, two oyster species with different tolerance levels to heat stress, led to the identification of TRPC3.6, TRPC3.7, and TRPV4.7 as important TRPs involved in thermal regulation in oysters. This work provided valuable information for future studies on the molecular mechanism of TRP mediated thermal tolerance, and identification of diagnostic biomarker for thermal stress in the oysters.

Phylogenetic comparisons reveal mosaic histories of larval and adult shell matrix protein deployment in pteriomorph bivalves.

Molluscan shells are organo-mineral composites, in which the dominant calcium carbonate is intimately associated with an organic matrix comprised mainly of proteins and polysaccharides. However, whether the various shell matrix proteins (SMPs) date to the origin of hard skeletons in the Cambrian, or whether they represent later deployment through adaptive evolution, is still debated. In order to address this issue and to better understand the origins and evolution of biomineralization, phylogenetic analyses have been performed on the three SMP families, Von Willebrand factor type A (VWA) and chitin-binding domain-containing protein (VWA-CB dcp), chitobiase, and carbonic anhydrase (CA), which exist in both larval and adult shell proteomes in the bivalves, Crassostrea gigas and Pinctada fucata. In VWA-CB dcp and chitobiase, paralogs for larval and adult SMPs evolved before the divergence of these species. CA-SMPs have been taken as evidence for ancient origins of SMPs by their presumed indispensable function in biomineralization and ubiquitous distribution in molluscs. However, our results indicate gene duplications that gave rise to separate deployments as larval and adult CA-SMPs occurred independently in each lineage after their divergence, which is considerably more recent than hitherto assumed, supporting the "recent heritage and fast evolution" scenario for SMP evolution.

A novel tumor necrosis factor in the Pacific oyster Crassostrea gigas mediates the antibacterial response by triggering the synthesis of lysozyme and nitric oxide.

Tumor necrosis factors (TNFs) are a group of multifunctional inflammatory cytokines involved in various pathological and immune processes. Recently, a few primitive TNFs have been characterized from molluscs, which play important roles in modulating cell apoptosis, phagocytosis and production of immune-related enzymes. In the present study, a novel TNF (named as CgTNF-2) with the activity to mediate antibacterial response was identified from the Pacific oyster Crassostrea gigas. The open reading frame of CgTNF-2 was of 783 bp encoding a putative polypeptide of 261 amino acids with a typical TNF domain. The deduced amino acid sequence of CgTNF-2 shared high identity with that of TNFs previously identified from other molluscs, such as 96.1% identity with that in oyster C. hongkongensis, 33.7% identity with that in scallop Mizuhopecten yessoensis and 33.0% identity with CgTNF-1 in oyster C. gigas. There were two distinct TNF branches of vertebrate and invertebrate in the phylogenetic tree, and CgTNF-2 was firstly clustered with TNF-14 from C. hongkongensis, and then clustered with other molluscan TNFs. The mRNA transcripts of CgTNF-2 were widely expressed in various oyster tissues, with the highest expression level in hemocytes. The expression level of CgTNF-2 increased significantly at 6 h (2.45-fold and 6.20-fold, respectively, p < 0.05) after peptidoglycan and lipopolysaccharides treatments, and peaked at 12 h (31.86-fold and 7.90-fold, respectively, p < 0.05). The recombinant protein of CgTNF-2 (rCgTNF-2) inhibited the growth of human alveolar basal epithelial (A549) cells at a concentration of 800 ng/mL. After the oysters received an injection of rCgTNF-2, the serum from those oysters exhibited significantly higher antibacterial activity compared to that from control group, evidenced by inhibiting the growth of Vibrio splendidus. Moreover, the lysozyme activity as well as the contents of nitric oxide in the oyster serum also increased significantly. The above results collectively suggested that CgTNF-2 was a novel member of bivalve TNF-α family, which could prompt the antibacterial activity by inducing the lysozyme activity and the production of nitric oxide in the innate immune response of oyster.

Development of SNP Panels as a New Tool to Assess the Genetic Diversity, Population Structure, and Parentage Analysis of the Eastern Oyster (Crassostrea virginica).

Culture of the eastern oyster (Crassostrea virginica) is rapidly expanding. Combined with their continuing role as an environmental sentinel species and ecological model, this trend necessitates improved molecular tools for breeding and selection, as well as population assessment and genetic conservation. Here, we describe the development and validation of two panels of 58 single nucleotide polymorphism markers (SNPs) for the species. Population analyses revealed three distinct populations, based on FST values and STRUCTURE, among wild oysters sampled from Delaware Bay (1), northwest Florida (2), Alabama (2), Louisiana (2), and the Texas Gulf Coast (3), consistent with previous microsatellite and mtDNA analyses. In addition, utilizing the developed panels for parentage assignment in cultured oysters (Rutgers, New Jersey) resulted in a highly accurate identification of parent pairs (99.37%). The SNP markers could, furthermore, clearly discriminate between hatchery stocks and wild-sourced individuals. The developed SNP panels may serve as an important tool for more rapid and affordable genetic analyses in eastern oyster.

Conservation and divergence of mitochondrial apoptosis pathway in the Pacific oyster, Crassostrea gigas.

Apoptosis is considered a crucial part of the host defense system in oysters according to previous reports; however, the exact process by which this occurs remains unclear. Besides, mitochondrial apoptosis is the primary method of apoptosis in vertebrate cells, but has been poorly studied in invertebrates and is quite controversial. In this study, we investigated the molecular mechanism of mitochondrial apoptosis in the Pacific oyster Crassostrea gigas. Notably, we show that most key elements involved in the vertebrate mitochondrial apoptosis pathway - including mitochondrial outer membrane permeabilization, cytochrome c release, and caspase activation - are also present in C. gigas. In contrast, the lack of Bcl-2 homology 3-only subfamily members and apoptotic protease activating factor-1 (APAF-1) protein revealed evolutionary diversity from other phyla. Our results support that mitochondrial apoptosis in animals predates the emergence of vertebrates, but suggest that an unexpectedly diverse mitochondrial apoptosis pathway may exist in invertebrates. In addition, our work provided new clues for an improved understanding of how bivalve acclimate themselves to an inconstant environment.

Molecular identification and functional characterization of a tumor necrosis factor (TNF) gene in Crassostrea hongkongensis.

Tumor necrosis factor superfamily (TNFSF) represents a group of multifunctional inflammatory cytokines that have been shown to participate in a variety of pathological and immunological process. However, the functions of these proteins in oyster are still poorly understood. In the present study, an oyster TNF homolog (named ChTNF) was identified from a cDNA library of Crassostrea hongkongensis. The complete cDNA of ChTNF was 2457bp in length containing an open reading frame (ORF) of 1044bp, a 5'-untranslated region (UTR) of 381bp and a 3'-UTR of 1032bp. The deduced ChTNF protein consisted of 347 amino acids with a characteristic transmembrane domain and a typical TNF homology domain (THD). Quantitative real-time PCR analysis revealed that ChTNF was broadly expressed in various oyster tissues and different stages of embryonic development. In addition, transcriptional analysis indicated that ChTNF transcription levels in hemocytes were increased significantly in pathogen challenge groups (Vibrio alginolyticus and Staphylococcus haemolyticus) compared to that in the control. Moreover, in vitro PAMP (lipopolysaccharide and peptidoglycan) treatments showed a stimulatory effect on the expression of ChTNF in the primary cultured hemocytes of C. hongkongensis. Finally, dual-luciferase reporter assays revealed that ChTNF could activate the NF-κB-Luc reporter in a dose-dependent manner in HEK293T cells. Altogether, these findings may provide valuable information regarding oyster TNFs and its role in innate immunity.

First data on three bivalve species exposed to an intra-harbour polymetallic contamination (La Rochelle, France).

Evaluating diffuse sediment contamination in the environment is a major concern with the aim of reaching a good chemical and ecological state of the littoral zone. In this study the risks of chronic chemical contamination and consequences in the bivalves Crassostrea gigas, Mytilus sp. and Mimachlamys varia were evaluated in coastal environments. The objective here was to understand the anthropological phenomena that affect the functioning of the marina of La Rochelle (semi-closed environment). Harbours seeking ecomanagement accreditations (such as the international reference ISO 14001) constitute zones of interest to implement biomonitoring studies. The biological effects of chemical pollution in the Marina of La Rochelle were studied to develop a multi-biomarker biomonitoring approach on specific marine species of this site. Moreover, a genetic (DNA barcoding) approach was applied to validate the species identity of collected bivalves. Of the three species tested the scallop, M. varia, was the most sensitive to metal exposure.

A Carbonic Anhydrase Serves as an Important Acid-Base Regulator in Pacific Oyster Crassostrea gigas Exposed to Elevated CO2: Implication for Physiological Responses of Mollusk to Ocean Acidification.

Carbonic anhydrases (CAs) have been demonstrated to play an important role in acid-base regulation in vertebrates. However, the classification and modulatory function of CAs in marine invertebrates, especially their responses to ocean acidification remain largely unknown. Here, a cytosolic α-CA (designated as CgCAII-1) was characterized from Pacific oyster Crassostrea gigas and its molecular activities against CO2 exposure were investigated. CgCAII-1 possessed a conserved CA catalytic domain, with high similarity to invertebrate cytoplasmic or mitochondrial α-CAs. Recombinant CgCAII-1 could convert CO2 to HCO3- with calculated activity as 0.54 × 103 U/mg, which could be inhibited by acetazolamide (AZ). The mRNA transcripts of CgCAII-1 in muscle, mantle, hepatopancreas, gill, and hemocytes increased significantly after exposure to elevated CO2. CgCAII-1 could interact with the hemocyte membrane proteins and the distribution of CgCAII-1 protein became more concentrated and dense in gill and mantle under CO2 exposure. The intracellular pH (pHi) of hemocytes under CO2 exposure increased significantly (p < 0.05) and CA inhibition reduced the pHi value. Besides, there was no increase in CA activity in gill and mantle after CO2 exposure. The impact of CO2 exposure on CA activity coupled with the mRNA expression level and protein translocation of CgCAII-1 provided evidences that CgCAII-1 could respond to ocean acidification and participate in acid-base regulation. Such cytoplasmic CA-based physiological regulation mechanism might explain other physiological responses of marine organisms to OA.

Biochemical alterations in native and exotic oyster species in Brazil in response to increasing temperature.

The increase of temperature in marine coastal ecosystems due to atmospheric greenhouse gas emissions is becoming an increasing threat for biodiversity worldwide, and may affect organisms' biochemical performance, often resulting in biogeographical shifts of species distribution. At the same time, the introduction of non-native species into aquatic systems also threatens biodiversity and ecosystem functions. Oysters are among the most valuable socio economic group of bivalve species in global fishery landings, and also provide numerous ecosystem services. However, the introduction of non-native oyster species, namely Crassostrea gigas for aquaculture purposes may threaten native oyster species, mainly by out competing their native congeners. It is therefore of upmost importance to understand physiological and biochemical responses of native and introduced oyster species in a scenario of global temperature rise, in order to provide knowledge that may allow for better species management. Hence, we compared biochemical alterations of the introduced C. gigas and the native Crassostrea brasiliana, the most important oyster species in Brazil, in response to different thermal regimes for 28days (24, 28 and 32°C). For this, metabolism (ETS), energy content (GLY), antioxidant system (SOD, CAT and GSH/GSSG) and cellular damage (LPO) were assessed in adult and juvenile specimens of both species. Juvenile C. gigas were the most affected by increased temperatures, presenting higher mortality, more pronounced antioxidant response (SOD), whereas adults were more tolerant than juveniles, showing no mortality, no significant changes in antioxidant enzymes activity neither energy expenditure. Native C. brasiliana juveniles presented lower mortality and less pronounced biochemical alterations were noted at higher temperature comparing to non-native C. gigas juveniles. Adult C. brasiliana were the least responsive to tested temperatures. Results obtained in this study bring interesting new insights on different oyster species life stages' physiological and biochemical tolerance towards thermal stress. The native species C. brasiliana showed ability to maintain biochemical performance at higher temperatures, with less pronounced biochemical changes than the non-native species. The introduced (C. gigas) showed to be more sensitive, presenting biochemical alterations to cope with the increase of temperature. Despite the lower observed fitness of the introduced species to temperatures closer to those naturally experienced by the native species, the ability of C. gigas to cope with higher temperatures should still raise concerns towards the native species C. brasiliana management and protection.

Molecular Cloning and Characterization of Full-Length cDNA of Calmodulin Gene from Pacific Oyster Crassostrea gigas.

The shell of the pearl oyster (Pinctada fucata) mainly comprises aragonite whereas that of the Pacific oyster (Crassostrea gigas) is mainly calcite, thereby suggesting the different mechanisms of shell formation between above two mollusks. Calmodulin (CaM) is an important gene for regulating the uptake, transport, and secretion of calcium during the process of shell formation in pearl oyster. It is interesting to characterize the CaM in oysters, which could facilitate the understanding of the different shell formation mechanisms among mollusks. We cloned the full-length cDNA of Pacific oyster CaM (cgCaM) and found that the cgCaM ORF encoded a peptide of 113 amino acids containing three EF-hand calcium-binding domains, its expression level was highest in the mantle, hinting that the cgCaM gene is probably involved in shell formation of Pacific oyster, and the common ancestor of Gastropoda and Bivalvia may possess at least three CaM genes. We also found that the numbers of some EF hand family members in highly calcified species were higher than those in lowly calcified species and the numbers of these motifs in oyster genome were the highest among the mollusk species with whole genome sequence, further hinting the correlation between CaM and biomineralization.

Coverage recommendation for genotyping analysis of highly heterologous species using next-generation sequencing technology.

Next-generation sequencing (NGS) technology is being applied to an increasing number of non-model species and has been used as the primary approach for accurate genotyping in genetic and evolutionary studies. However, inferring genotypes from sequencing data is challenging, particularly for organisms with a high degree of heterozygosity. This is because genotype calls from sequencing data are often inaccurate due to low sequencing coverage, and if this is not accounted for, genotype uncertainty can lead to serious bias in downstream analyses, such as quantitative trait locus mapping and genome-wide association studies. Here, we used high-coverage reference data sets from Crassostrea gigas to simulate sequencing data with different coverage, and we evaluate the influence of genotype calling rate and accuracy as a function of coverage. Having initially identified the appropriate parameter settings for filtering to ensure genotype accuracy, we used two different single-nucleotide polymorphism (SNP) calling pipelines, single-sample and multi-sample. We found that a coverage of 15× was suitable for obtaining sufficient numbers of SNPs with high accuracy. Our work provides guidelines for the selection of sequence coverage when using NGS to investigate species with a high degree of heterozygosity and rapid decay of linkage disequilibrium.

High expression of new genes in trochophore enlightening the ontogeny and evolution of trochozoans.

Animals with trochophore larvae belong to Trochozoa, one of the main branches of Bilateria. In addition to exhibiting spiral cleavage and early cell fate determination, trochozoans typically undergo indirect development, which contributes to the most unique characteristics of their ontogeny. The indirect development of trochozoans has provoked discussion regarding the origin and evolution of marine larvae and is interesting from the perspective of phylogeny-ontogeny correspondence. While these phylo-onto correlations have an hourglass shape in Deuterostomia, Ecdysozoa, plants and even fungi, they have seldom been studied in Trochozoa, and even Lophotrochozoa. Here, we compared the ontogenetic transcriptomes of the Pacific oyster, Crassostrea gigas (Bivalvia, Mollusca), the Pacific abalone, Haliotis discus hannai (Gastropoda, Mollusca), and the sand worm Perinereis aibuhitensis (Polychaeta, Annelida) using several complementary phylotranscriptomic methods to examine their evolutionary trajectories. The results revealed the late trochophore stage as the phylotypic phase. However, this basic pattern is accompanied with increased use of new genes in the trochophore stages which marks specific adaptations of the larval body plans.

Cloning, characterization and comparative analysis of four death receptorTNFRs from the oyster Crassostrea hongkongensis.

Apoptosis plays an important role in homeostasis of the immune systems. The tumor necrosis factor receptors (TNFRs) play critical roles in the extrinsic apoptosis pathways and in determining cell fate. In this study, four death receptors (DR) named ChEDAR, ChTNFR27, ChTNFR5, and ChTNFR16 were identified from the oyster Crassostrea hongkongensis. These ChDRs proteins had 382, 396, 414 and 384 amino acids, respectively, with the typical domains of death receptors, such as the signal peptide (SP), transmembrane helix region (TM) and death domains. Phylogenetic analysis showed that the ChDR proteins clustered into three distinct groups, indicating that these subfamilies had common ancestors. mRNA expression of the ChDRs were detected in all 8 of the selected oyster tissues and at different stages of development. Furthermore, expression of all the genes was increased in the hemocytes of oysters challenged with pathogens or air stress. Fluorescence microscopy revealed that the full-length proteins of the ChDRs were located in the plasma membrane of HEK293T cells. Over-expression of the ChDRs activated the NF-κB-Luc reporter in HEK293T cells in a dose-dependent manner. These results indicate that the ChDRs may play important roles in the extrinsic apoptotic pathways in oysters.

Efecto de la temperatura y la salinidad en el ciclo reproductor de hembras y machos de Crassostrea virginica (Bivalvia: Ostreidae) / Effect of temperature and salinity on the reproductive cycle of males and females of Crassostrea virginica (Bivalvia: Ostreidae)

ResumenLa reproducción del ostión americano, Crassostrea virginica se encuentra ligada a la temperatura y la salinidad. Sin embargo el efecto de estos parámetros no ha sido evaluado por sexos y por tallas. En este trabajo se analiza el efecto de la temperatura y la salinidad sobre la reproducción del ostión americano, por sexos separados y por tallas (40-60mm y ≥ 61 mm). Se recolectaron 1 080 ostiones, 90 organismos por mes durante 2011 en la laguna de Tamiahua, Veracruz México. El desarrollo de madurez gonádica se clasificó en tres estadios: gametogénesis, maduros y liberación. La reproducción fue continua en ambos sexos; la gametogénesis presentó máximos entre enero-julio, la madurez entre junio y septiembre y la liberación en enero y febrero, y de septiembre a diciembre. Sólo se presentó correlación significativa entre la temperatura y las hembras de ambas tallas en madurez, pero no en machos. La disminución de la salinidad fue determinante en la liberación de gametos pero sólo para los machos ≥ 61 mm, en machos pequeños y hembras de ambas tallas, este parámetro no tuvo efecto. En este trabajo se puso en evidencia que la temperatura y la salinidad tuvieron efecto diferencial entre tallas y sexos del ostión americano.

AbstractThe American oyster, Crassostrea virginica, is one of the most important economic coastal fishery activities in Mexico. As for other bivalve molluscs, its sustainability is based on healthy populations, and additional information on current reproductive cycles, will reinforce local management strategies. The objective of this study was to evaluate the effect of temperature and salinity on the gonadal development stages in both sexes and in two size groups (40-60 mm and ≥ 61 mm) of C. virginica from Tamiahua Lagoon, Gulf of Mexico. Monthly surveys of 90 organisms were undertaken, for a yearly total of 1 080 oysters sampled from the lagoon in 2011. Both field and laboratory work used standard procedures. Our results showed that reproduction was continuous in both sexes, but some peaks of the reproductive cycle were observed along the year. Peaks registered from January to July were for gametogenesis; those of June and September were for mature individuals; and in January, February, and from September to December, for spawning/spent oysters. Maturity of females of both range sizes was positively correlated with temperature, but not for males. For spawning, spent individuals, males of sizes ≥ 61 mm were positively correlated with salinity, whereas smaller males (40-60 mm) and females of both sizes, did not show any correlation. In this study, we observed that temperature and salinity had a different effect in the reproduction of the American oyster of different sizes and sexes of Tamiahua lagoon. Rev. Biol. Trop. 64 (2): 449-459. Epub 2016 June 01.

The complete mitochondrial genome of Crassostrea gasar (Bivalvia: Ostreidae).

The complete mitochondrial genome of Crassostrea gasar was sequenced using the Ion Proton technology in combination with 454 Roche GS-FLX plataform data. We assembled a 17,686 bp complete circular mitochondrial genome, containing 13 protein-coding genes, a major non-coding region (MNR), two ribosomal RNA genes and 24 transfer RNA genes. Phylogenetic analysis of concatenated amino acid sequences from mitochondria showed monophyletic clades formed with high bootstrap values. This is the first complete mitochondrial sequence of an oyster from South America. Mitogenome sequence was deposited in GenBank under the accession number KR856227.

Characterization of the gene encoding a fibrinogen-related protein expressed in Crassostrea gigas hemocytes.

Four exons of the CgFrep1 gene (3333 bp long) encode a putative fibrinogen-related protein (324 aa) bearing a single C-terminal FBG domain. Transcripts of the gene obtained from hemocytes of different Pacific oysters show prominent individual variation based on SNP and indels of tandem repeats resulted in polymorphism of N-terminus of the putative CgFrep1 polypeptide. The polypeptide chain bears N-terminal coiled-coil region potentially acting as inter-subunit interface in the protein oligomerization. It is suggested that CgFrep1 gene encodes the oligomeric lectin composed of at least two subunits.

Intraspecific Variation in Mitogenomes of Five Crassostrea 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.

Expression Characterization of Stress Genes Under High and Low Temperature Stresses in the Pacific Oyster, Crassostrea gigas.

As a characteristic sessile inhabitant of the intertidal zone, the Pacific oyster Crassostrea gigas occupies one of the most physically stressful environments on earth. With high exposure to terrestrial conditions, oysters must tolerate broad fluctuations in temperature range. However, oysters' cellular and molecular responses to temperature stresses have not been fully characterized. Here, we analyzed oyster transcriptome data under high and low temperatures. We also identified over 30 key temperature stress-responsive candidate genes, which encoded stress proteins such as heat shock proteins and apoptosis-associated proteins. The expression characterization of these genes under short-term cold and hot environments (5 and 35 °C) and long-term cold environments (5 °C) was detected by quantitative real-time PCR. Most of these genes reached expression peaks during the recovery stage after 24 h of heat stress, and these genes were greatly induced around day 3 in long-term cold stress while responded little to short-term cold stress. In addition, in the second heat stress after 2 days of recovery, oysters showed milder expression in these genes and a lower mortality rate, which indicated the existence of plasticity in the oyster's response to heat stress. We confirmed that homeostatic flexibility and anti-apoptosis might be crucial centers of temperature stress responses in oysters. Furthermore, we analyzed stress gene families in 11 different species and found that the linage-specific expansion of stress genes might be implicated in adaptive evolution. These results indicated that both plasticity and evolution played an important role in the stress response adaptation of oysters.