Huayong Que

The oyster genome reveals stress adaptation and complexity of shell formation

The Pacific oyster Crassostrea gigas belongs to one of the most species-rich but genomically poorly explored phyla, the Mollusca. Here we report the sequencing and assembly of the oyster genome using short reads and a fosmid-pooling strategy, along with transcriptomes of development and stress response and the proteome of the shell. The oyster genome is highly polymorphic and rich in repetitive sequences, with some transposable elements still actively shaping variation. Transcriptome studies reveal an extensive set of genes responding to environmental stress. The expansion of genes coding for heat shock protein 70 and inhibitors of apoptosis is probably central to the oyster’s adaptation to sessile life in the highly stressful intertidal zone. Our analyses also show that shell formation in molluscs is more complex than currently understood and involves extensive participation of cells and their exosomes. The oyster genome sequence fills a void in our understanding of the Lophotrochozoa.

Oyster Shell Proteins Originate from Multiple Organs and Their Probable Transport Pathway to the Shell Formation Front

Mollusk shell is one kind of potential biomaterial, but its vague mineralization mechanism hinders its further application. Mollusk shell matrix proteins are important functional components that are embedded in the shell, which play important roles in shell formation. The proteome of the oyster shell had been determined based on the oyster genome sequence by our group and gives the chance for further deep study in this area. The classical model of shell formation posits that the shell proteins are mantle-secreted. But, in this study, we further analyzed the shell proteome data in combination with organ transcriptome data and we found that the shell proteins may be produced by multiple organs though the mantle is still the most important organ for shell formation. To identify the transport pathways of these shell proteins not in classical model of shell formation, we conducted a shell damage experiment and we determined the shell-related gene set to identify the possible transport pathways from multiple organs to the shell formation front. We also found that there may exist a remodeling mechanism in the process of shell formation. Based on these results along with some published results, we proposed a new immature model, which will help us think about the mechanism of shell formation in a different way.

Discovery and validation of genic single nucleotide polymorphisms in the Pacific oyster Crassostrea gigas

The economic and ecological importance of the oyster necessitates further research on the molecular mechanisms, which both regulate the commercially important traits of the oyster and help it to survive in the variable marine environment. Single nucleotide polymorphisms (SNPs) have been widely used to assess genetic variation and identify genes underlying target traits. In addition, high‐resolution melting (HRM) analysis is a potentially powerful method for validating candidate SNPs. In this study, we adopted a rapid and efficient pipeline for the screening and validation of SNPs in the genic region of Crassostrea gigas based on transcriptome sequencing and HRM analysis. Transcriptomes of three wild oyster populations were sequenced using Illumina sequencing technology. In total, 50–60 million short reads, corresponding to 4.5–5.4 Gbp, from each population were aligned to the oyster genome, and 5.8 × 105 SNPs were putatively identified, resulting in a predicted SNP every 47 nucleotides on average. The putative SNPs were unevenly distributed in the genome and high‐density (≥2%), nonsynonymous coding SNPs were enriched in genes related to apoptosis and responses to biotic stimuli. Subsequently, 1,671 loci were detected by HRM analysis, accounting for 64.7% of the total selected candidate primers, and finally, 1,301 polymorphic SNP markers were developed based on HRM analysis. All of the validated SNPs were distributed into 897 genes and located in 672 scaffolds, and 275 of these genes were stress inducible under unfavourable salinity, temperature, and exposure to air and heavy metals. The validated SNPs in this study provide valuable molecular markers for genetic mapping and characterization of important traits in oysters.

Identification of Conserved and Novel MicroRNAs in the Pacific Oyster Crassostrea gigas by Deep Sequencing

MicroRNAs (miRNAs) play important roles in regulatory processes in various organisms. To date many studies have been performed in the investigation of miRNAs of numerous bilaterians, but limited numbers of miRNAs have been identified in the few species belonging to the clade Lophotrochozoa. In the current study, deep sequencing was conducted to identify the miRNAs of Crassostrea gigas (Lophotrochozoa) at a genomic scale, using 21 libraries that included different developmental stages and adult organs. A total of 100 hairpin precursor loci were predicted to encode miRNAs. Of these, 19 precursors (pre-miRNA) were novel in the oyster. As many as 53 (53%) miRNAs were distributed in clusters and 49 (49%) precursors were intragenic, which suggests two important biogenetic sources of miRNAs. Different developmental stages were characterized with specific miRNA expression patterns that highlighted regulatory variation along a temporal axis. Conserved miRNAs were expressed universally throughout different stages and organs, whereas novel miRNAs tended to be more specific and may be related to the determination of the novel body plan. Furthermore, we developed an index named the miRNA profile age index (miRPAI) to integrate the evolutionary age and expression levels of miRNAs during a particular developmental stage. We found that the swimming stages were characterized by the youngest miRPAIs. Indeed, the large-scale expression of novel miRNAs indicated the importance of these stages during development, particularly from organogenetic and evolutionary perspectives. Some potentially important miRNAs were identified for further study through significant changes between expression patterns in different developmental events, such as metamorphosis. This study broadened the knowledge of miRNAs in animals and indicated the presence of sophisticated miRNA regulatory networks related to the biological processes in lophotrochozoans.

Cloning and characterization of a novel caspase-8-like gene in Crassostrea gigas

Cysteine-dependent aspartate-directed proteases, or caspases, play key roles in apoptosis and immune defense. In this study, we cloned the first caspase-8-like gene (CgCaspase8-2) identified in the pacific oyster, Crassostrea gigas. The 2572-bp cDNA encodes a putative protein of 714 amino acids that contains two tandem death effector domains (DEDs) at the N-terminal, and P20 and P10 domains at the C-terminal. The conserved pentapeptide motif QACQG was also identified in the deduced CgCaspase8-2 protein. Phylogenetic analysis indicated that CgCaspase8-2 was clustered with initiator caspases in the invertebrate subgroup, but the similarity between CgCaspase8-2 and other invertebrate caspase-8s was low. CgCaspase8-2 protein was localized in the cytoplasm, and over-expression of CgCaspase8-2 in HEK293T cells induced cell death, suggesting a role in apoptosis. Quantitative real-time PCR results demonstrated that CgCaspase8-2 was widely expressed in various tissues and developmental stages, with the highest CgCaspase8-2 expression levels detected in hemolymph and the blastula stage. Furthermore, CgCaspase8-2 transcripts showed no change in response to a bacterial challenge but exhibited notable up-regulation post-poly (I:C) challenge, suggesting that CgCaspase8-2 is specifically involved in immune responses against viruses. In summary, CgCaspase8-2 is involved in both apoptotic and immune function.

Identification of Thyroid Hormones and Functional Characterization of Thyroid Hormone Receptor in the Pacific Oyster Crassostrea gigas Provide Insight into Evolution of the Thyroid Hormone System

Thyroid hormones (THs) play important roles in development, metamorphosis, and metabolism in vertebrates. During the past century, TH functions were regarded as a synapomorphy of vertebrates. More recently, accumulating evidence has gradually convinced us that TH functions also occur in invertebrate chordates. To date, however, TH-related studies in non-chordate invertebrates have been limited. In this study, THs were qualitatively detected by two reliable methods (HPLC and LC/MS) in a well-studied molluscan species, the Pacific oyster Crassostrea gigas. Quantitative measurement of THs during the development of C. gigas showed high TH contents during embryogenesis and that oyster embryos may synthesize THs endogenously. As a first step in elucidating the TH signaling cascade, an ortholog of vertebrate TH receptor (TR), the most critical gene mediating TH effects, was cloned in C. gigas. The sequence of CgTR has conserved DNA-binding and ligand-binding domains that normally characterize these receptors. Experimental results demonstrated that CgTR can repress gene expression through binding to promoters of target genes and can interact with oyster retinoid X receptor. Moreover, CgTR mRNA expression was activated by T4 and the transcriptional activity of CgTR promoter was repressed by unliganded CgTR protein. An atypical thyroid hormone response element (CgDR5) was found in the promoter of CgTR, which was verified by electrophoretic mobility shift assay (EMSA). These results indicated that some of the CgTR function is conserved. However, the EMSA assay showed that DNA binding specificity of CgTR was different from that of the vertebrate TR and experiments with two dual-luciferase reporter systems indicated that l-thyroxine, 3,3′,5-triiodothyronine, and triiodothyroacetic acid failed to activate the transcriptional activity of CgTR. This is the first study to functionally characterize TR in mollusks. The presence of THs and the functions of CgTR in mollusks contribute to better understanding of the evolution of the TH system.

Candidate Gene Polymorphisms and their Association with Glycogen Content in the Pacific Oyster Crassostrea gigas.

Background The Pacific oyster Crassostrea gigas is an important cultivated shellfish that is rich in nutrients. It contains high levels of glycogen, which is of high nutritional value. To investigate the genetic basis of this high glycogen content and its variation, we conducted a candidate gene association analysis using a wild population, and confirmed our results using an independent population, via targeted gene resequencing and mRNA expression analysis. Results We validated 295 SNPs in the 90 candidate genes surveyed for association with glycogen content, 86 of were ultimately genotyped in all 144 experimental individuals from Jiaonan (JN). In addition, 732 SNPs were genotyped via targeted gene resequencing. Two SNPs (Cg_SNP_TY202 and Cg_SNP_3021) in Cg_GD1 (glycogen debranching enzyme) and one SNP (Cg_SNP_4) in Cg_GP1 (glycogen phosphorylase) were identified as being associated with glycogen content. The glycogen content of individuals with genotypes TT and TC in Cg_SNP_TY202 was higher than that of individuals with genotype CC. The transcript abundance of both glycogen-associated genes was differentially expressed in high glycogen content and low glycogen content individuals. Conclusions This study identified three polymorphisms in two genes associated with oyster glycogen content, via candidate gene association analysis. The transcript abundance differences in Cg_GD1 and Cg_GP1 between low- and the high-glycogen content individuals suggests that it is possible that transcript regulation is mediated by variations of Cg_SNP_TY202, Cg_SNP_3021, and Cg_SNP_4. These findings will not only provide insights into the genetic basis of oyster quality, but also promote research into the molecular breeding of oysters.

Iodothyronine deiodinase gene analysis of the Pacific oyster Crassostrea gigas reveals possible conservation of thyroid hormone feedback regulation mechanism in mollusks

Iodothyronine deiodinase catalyzes the initiation and termination of thyroid hormones (THs) effects, and plays a central role in the regulation of thyroid hormone level in vertebrates. In non-chordate invertebrates, only one deiodinase has been identified in the scallop Chlamys farreri. Here, two deiodinases were cloned in the Pacific oyster Crassostrea gigas (CgDx and CgDy). The characteristic in-frame TGA codons and selenocysteine insertion sequence elements in the oyster deiodinase cDNAs supported the activity of them. Furthermore, seven orthologs of deiodinases were found by a tblastn search in the mollusk Lottia gigantea and the annelid Capitella teleta. A phylogenetic analysis revealed that the deiodinase gene originated from an common ancestor and a clade-specific gene duplication occurred independently during the differentiation of the mollusk, annelid, and vertebrate lineages. The distinct spatiotemporal expression patterns implied functional divergence of the two deiodinases. The expression of CgDx and CgDy was influenced by L-thyroxine T4, and putative thyroid hormone responsive elements were found in their promoters, which suggested that the oyster deiodinases were feedback regulated by TH. Epinephrine stimulated the expression level of CgDx and CgDy, suggesting an interaction effect between different hormones. This study provides the first evidence for the existence of a conserved TH feedback regulation mechanism in mollusks, providing insights into TH evolution.

Phylogeny of forkhead genes in three spiralians and their expression in Pacific oyster Crassostrea gigas

The Fox genes encode a group of transcription factors that contain a forkhead domain, which forms a structure known as a winged helix. These transcription factors play a crucial role in several key biological processes, including development. High-degree identity in the canonical forkhead domain has been used to divide Fox proteins into 23 families (FoxA to FoxS). We surveyed the genome of three spiralians, the oyster Crassostrea gigas, the limpet Lottia gigantea, and the annelid Capitella teleta. We identified 25 C. gigas fox genes, 21 L. gigantea fox genes, and 25 C. teleta fox genes. The C. gigas fox and L. gigantea fox genes represented 19 of the 23 families, whereas FoxI, Q1, R, and S were missing. The majority of the Fox families were observed within the C. teleta fox genes, with the exception of FoxR and S. In addition, the foxAB-like gene, foxY-like gene, and foxH gene were also present in the three genomes. The conserved FoxC-FoxL1 cluster, observed in mammals, was also found in C. gigas. The diversity of temporal expression patterns observed across the developmental process implies the C. gigas fox genes exert a wide range of functions. Further functional studies are required to gain insight into the evolution of Fox genes in bilaterians.

Evolution of a novel nuclear receptor subfamily with emphasis on the member from the Pacific oyster Crassostrea gigas

Nuclear receptors (NRs) belong to the transcription factor superfamily that regulates development, homeostasis, differentiation, and reproduction in metazoans via control of gene expression. Recently, rapid advances in genome projects on various metazoans have provided new opportunities for studying the evolution and function of NRs. Typically structured NRs are divided into six subfamilies. Here, the gene for a typically structured NR (CgNR8A1) was cloned from the Pacific oyster Crassostrea gigas. However, this novel receptor could not be assigned to a known NR subfamily. By data mining, nine other CgNR8A1 gene homologs were identified in metazoans such as cnidarians, mollusks, annelids, echinoderms, hemichordates, and cephalochordates. Phylogenetic analysis showed that these receptors belonged to a novel NR subfamily, hereafter designated as NR8. Evolutionary analysis revealed that the NR8 subfamily was phylogenetically the third-oldest NR subfamily, and it originated from a common ancestor of Eumetazoa; several gene loss events occurred independently in ancestors of vertebrates, ecdysozoans, and platyhelminths, which do not have NR8 members. Furthermore, the function of CgNR8A1 was investigated to provide an insight into the functions of this novel NR subfamily. A nuclear localization signal peptide, GKHRNKKPRLD, was identified in CgNR8A1, and a recombinant full-length protein of CgNR8A1 was localized in the nuclei of HeLa cells. The mRNA expression profile of CgNR8A1 suggested that it might be involved in the embryogenesis of C. gigas. The electrophoretic mobility shift assay showed that CgNR8A1 binds strongly to conserved DNA core motifs DR0, DR2, and DR4 and weakly to DR1, DR3, DR5, Half, and Pal0. In summary, the novel NR8 subfamily identified in this study improves our understanding of NR evolution, and the functional analysis of CgNR8A1 provided further insights into the functions of NR8A1s.