Wen Huang

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.

Identification and functional characterization of two executioner caspases in Crassostrea gigas.

Caspase-3 and caspase-7 are two key effector caspases that play important roles in apoptotic pathways that maintain normal tissue and organ development and homeostasis. However, little is known about the sequence, structure, activity, and function of effector caspases upon apoptosis in mollusks, especially marine bivalves. In this study, we investigated the possible roles of two executioner caspases in the regulation of apoptosis in the Pacific oyster Crassostrea gigas. A full-length capase-3–like gene named Cgcaspase-3 was cloned from C.gigas cDNA, encoding a predicted protein containing caspase family p20 and p10 domain profiles and a conserved caspase active site motif. Phylogenetic analysis demonstrated that both Cgcaspase-3 and Cgcaspase-1 may function as effector caspases clustered in the invertebrate branch. Although the sequence identities between the two caspases was low, both enzymes possessed executioner caspase activity and were capable of inducing cell death. These results suggested that Cgcaspase-3 and Cgcaspase-1 were two effector caspases in C. gigas. We also observed that nucleus-localized Cgcaspase-3, may function as a caspase-3–like protein and cytoplasm-localized Cgcaspase-1 may function as a caspase-7–like protein. Both Cgcaspase-3 and Cgcaspase-1 mRNA expression increased after larvae settled on the substratum, suggesting that both caspases acted in several tissues or organs that degenerated after oyster larvae settlement. The highest caspase expression levels were observed in the gills indicating that both effector caspases were likely involved in immune or metabolic processes in C. gigas.

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.

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.

A new L-type lectin ( Lv LTLC1) from the shrimp Litopenaeus vannamei facilitates the clearance of Vibrio harveyi

ABSTRACT According to structures and functions of lectins found in shrimp, they are classified into seven types, namely, L‐type, C‐type, P‐type, M‐type, galectins, fibrinogen‐like domain lectins, and calnexin/calreticulin. Until now, the researches of shrimp lectins are mainly focused on C‐type lectins. In this study, we identified a new L‐type lectin, designated as LvLTLC1, from the shrimp Litopenaeus vannamei. The cDNA of LvLTLC1 is 1184 bp with an open reading frame of 990 bp encoding a protein of 329 amino acids. The LvLTLC1 protein contained a putative signal peptide, an L‐type lectin‐like domain, and a transmembrane helix region. Phylogenetic analysis showed that LvLTLC1 belonged to VIP36‐like family. LvLTLC1 was expressed in all examined tissues but had higher expression level in gills and hepatopancreas than other tissues. LvLTLC1 expression was up‐regulated after immune challenge by Vibrio harveyi and lipopolysaccharide. The recombinant LvLTLC1 agglutinated Gram‐positive (Staphylococcus aureus) and Gram‐negative bacteria (V. harveyi, V. parahaemolyticus, V. alginolyticus, V. cholerae, V. vulnificus, Pseudomonas aeruginosa, P. fluorescens) in a calcium‐independent manner. Recombinant LvLTLC1 exerted the ability of enhancing the clearance of V. harveyi injected in shrimp. Our results indicated that LvLTLC1 functions in anti‐pathogen innate immunity of shrimp. HIGHLIGHTSA new L‐type lectin (LvLTLC1) from the shrimp Litopenaeus vannamei was identified.LvLTLC1 agglutinates various bacteria in vitro in a calcium‐independent manner.LvLTLC1 enhances the clearance of invading bacteria (V. harveyi) in vivo.

A fibrinogen-related protein, Lv FREP2, from Litopenaeus vannamei facilitates the clearance of Vibrio harveyi

ABSTRACT Fibrinogen‐related proteins (FREPs) play a crucial role in invertebrate immune response. In this study, we acquired a novel fibrinogen‐related protein gene in Litopenaeus vannamei coding for one kind of fibrinogen‐related protein, designated as LvFREP2. The complete cDNA sequence of LvFREP2 was 1903 bp long, containing an open reading frame of 1479 bp coding for LvFREP2. The LvFREP2 protein contained a putative signal peptide and a fibrinogen‐related protein domain. qRT‐PCRs indicated that LvFREP2 mRNA ubiquitously distributed in all examined tissues, and it was up‐regulated in gills after V. harveyi and LPS challenges. The recombinant LvFREP2 agglutinated Gram‐positive (Staphylococcus aureus) and Gram‐negative bacteria (Vibrio alginolyticus, V. cholerae, V. vulnificus, V. parahaemolyticus, V. harveyi, Pseudomonas aeruginosa, P. fluorescens) in a calcium‐dependent manner. LvFREP2 also facilitated the clearance of Vibrio harveyi in vivo. Therefore, our results suggested that lvFREP2 may have important roles in the anti‐bacterial immunity of L. vannamei.

Transcriptomic analyses on muscle tissues of Litopenaeus vannamei provide the first profile insight into the response to low temperature stress

The Pacific white shrimp (Litopenaeus vannamei) is an important cultured crustacean species worldwide. However, little is known about the molecular mechanism of this species involved in the response to cold stress. In this study, four separate RNA-Seq libraries of L. vannamei were generated from 13°C stress and control temperature. Total 29,662 of Unigenes and overall of 19,619 annotated genes were obtained. Three comparisons were carried out among the four libraries, in which 72 of the top 20% of differentially-expressed genes were obtained, 15 GO and 5 KEGG temperature-sensitive pathways were fished out. Catalytic activity (GO: 0003824) and Metabolic pathways (ko01100) were the most annotated GO and KEGG pathways in response to cold stress, respectively. In addition, Calcium, MAPK cascade, Transcription factor and Serine/threonine-protein kinase signal pathway were picked out and clustered. Serine/threonine-protein kinase signal pathway might play more important roles in cold adaptation, while other three signal pathway were not widely transcribed. Our results had summarized the differentially-expressed genes and suggested the major important signaling pathways and related genes. These findings provide the first profile insight into the molecular basis of L. vannamei response to cold stress.

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.