Kazuyoshi Endo

Draft Genome of the Pearl Oyster Pinctada fucata: A Platform for Understanding Bivalve Biology

The study of the pearl oyster Pinctada fucata is key to increasing our understanding of the molecular mechanisms involved in pearl biosynthesis and biology of bivalve molluscs. We sequenced ∼1150-Mb genome at ∼40-fold coverage using the Roche 454 GS-FLX and Illumina GAIIx sequencers. The sequences were assembled into contigs with N50 = 1.6 kb (total contig assembly reached to 1024 Mb) and scaffolds with N50 = 14.5 kb. The pearl oyster genome is AT-rich, with a GC content of 34%. DNA transposons, retrotransposons, and tandem repeat elements occupied 0.4, 1.5, and 7.9% of the genome, respectively (a total of 9.8%). Version 1.0 of the P. fucata draft genome contains 23 257 complete gene models, 70% of which are supported by the corresponding expressed sequence tags. The genes include those reported to have an association with bio-mineralization. Genes encoding transcription factors and signal transduction molecules are present in numbers comparable with genomes of other metazoans. Genome-wide molecular phylogeny suggests that the lophotrochozoan represents a distinct clade from ecdysozoans. Our draft genome of the pearl oyster thus provides a platform for the identification of selection markers and genes for calcification, knowledge of which will be important in the pearl industry.

Biphasic and Dually Coordinated Expression of the Genes Encoding Major Shell Matrix Proteins in the Pearl Oyster Pinctada fucata

Regional expression patterns of shell matrix protein genes of Pinctada fucata were investigated using real-time quantitative polymerase chain reaction (PCR) and in situ hybridization. Six shell matrix proteins examined in this study indicated a distinct biphasic pattern of expression, falling into one of the following three groups: (1) expressed only in the more dorsal region of the mantle (MSI60 and N16); (2) expressed only in the more ventral region (MSI31, Prismalin-14, and Aspein); and (3) expressed in both regions (nacrein). The ubiquity of the last protein probably reflects its general role as a carbonate-producing enzyme, while the other groups are interpreted as corresponding to the distinction between the two varieties of shell layers, the aragonitic nacreous layer and the calcitic prismatic layer. In addition, the constituent genes of each of these two groups indicated similar levels of relative expression among different sites even among different individuals, suggesting that the genes of each group share a single upstream regulatory factor, respectively, and that these genes are expressed in a dually coordinated fashion.

The complete primary structure of molluscan shell protein 1 (MSP-1), an acidic glycoprotein in the shell matrix of the scallop Patinopecten yessoensis.

Abstract: The complete primary structure of MSP-1, a major water-soluble glycoprotein in the foliated calcite shell layer of the scallop Patinopecten yessoensis, is reported. The full-length complementary DNA for MSP-1 isolated by polymerase chain reaction contained a sequence for a signal peptide of 20 amino acids followed by a polypeptide of 820 amino acids with calculated molecular mass of 74.5 kDa. The deduced amino acid sequence of MSP-1 includes a high proportion of Ser (32%), Gly (25%), and Asp (20%), and the predicted isoelectric point is 3.2; in these respects, MSP-1 is a typical acidic glycoprotein of mineralized tissues. A repeated modular structure characterizes MSP-1, with a sequence unit between 158 and 177 amino acids in length being repeated 4 times in tandem in the middle part of the protein. The repeated unit comprises 3 modules (SG, D, and K domains), each having a distinct amino acid composition and sequence. The SG domain is almost exclusively composed of Ser and Gly residues. The D domain is rich in Asp residues, potential N-glycosylation and phosphorylation sites. The K domain is rich in Gly residues and has a core of basic residues. The Asp residues are arranged more or less regularly in the D domains, exhibiting some repeated motifs such as Asp-Gly-Ser-Asp and Asp-Ser-Asp. Further, the 4 D domains indicate remarkable overall sequence similarities to each other. These observations suggest that the regular arrangements of COO− groups in the D domain side chains may be important for specific control of crystal growth.

In vitro regulation of CaCO3 crystal polymorphism by the highly acidic molluscan shell protein Aspein

Biominerals, especially molluscan shells, generally contain unusually acidic proteins. These proteins are believed to function in crystal nucleation and inhibition. We previously identified an unusually acidic protein Aspein from the pearl oyster Pinctada fucata. Here we show that Aspein can control the CaCO3 polymorph (calcite/aragonite) in vitro. While aragonite is preferentially formed in Mg2+‐rich solutions imitating the extrapallial fluids of marine molluscs, Aspein exclusively induced calcite precipitation. Our results suggest that Aspein is involved in the specific calcite formation in the prismatic layer. Experiments using truncated Aspein demonstrated that the aspartic acid rich domain is crucial for the calcite precipitation.

The Lingula genome provides insights into brachiopod evolution and the origin of phosphate biomineralization.

The evolutionary origins of lingulid brachiopods and their calcium phosphate shells have been obscure. Here we decode the 425-Mb genome of Lingula anatina to gain insights into brachiopod evolution. Comprehensive phylogenomic analyses place Lingula close to molluscs, but distant from annelids. The Lingula gene number has increased to ∼34,000 by extensive expansion of gene families. Although Lingula and vertebrates have superficially similar hard tissue components, our genomic, transcriptomic and proteomic analyses show that Lingula lacks genes involved in bone formation, indicating an independent origin of their phosphate biominerals. Several genes involved in Lingula shell formation are shared by molluscs. However, Lingula has independently undergone domain combinations to produce shell matrix collagens with EGF domains and carries lineage-specific shell matrix proteins. Gene family expansion, domain shuffling and co-option of genes appear to be the genomic background of Lingulas unique biomineralization. This Lingula genome provides resources for further studies of lophotrochozoan evolution.

The Diversity of Shell Matrix Proteins: Genome-Wide Investigation of the Pearl Oyster, Pinctada fucata

In molluscs, shell matrix proteins are associated with biomineralization, a biologically controlled process that involves nucleation and growth of calcium carbonate crystals. Identification and characterization of shell matrix proteins are important for better understanding of the adaptive radiation of a large variety of molluscs. We searched the draft genome sequence of the pearl oyster Pinctada fucata and annotated 30 different kinds of shell matrix proteins. Of these, we could identified Perlucin, ependymin-related protein and SPARC as common genes shared by bivalves and gastropods; however, most gastropod shell matrix proteins were not found in the P. fucata genome. Glycinerich proteins were conserved in the genus Pinctada. Another important finding with regard to these annotated genes was that numerous shell matrix proteins are encoded by more than one gene; e.g., three ACCBP-like proteins, three CaLPs, five chitin synthase-like proteins, two N16 proteins (pearlins), 10 N19 proteins, two nacreins, four Pifs, nine shematrins, two prismalin-14 proteins, and 21 tyrosinases. This diversity of shell matrix proteins may be implicated in the morphological diversity of mollusc shells. The annotated genes reported here can be searched in P. fucata gene models version 1.1 and genome assembly version 1.0 (http://marinegenomics.oist.jp/pinctada_fucata). These genes should provide a useful resource for studies of the genetic basis of biomineralization and evaluation of the role of shell matrix proteins as an evolutionary toolkit among the molluscs.

Molecular evolution and functionally important structures of molluscan Dermatopontin: implications for the origins of molluscan shell matrix proteins.

A major shell matrix protein originally obtained from a freshwater snail is a molluscan homologue of Dermatopontins, a group of Metazoan proteins also called TRAMP (tyrosine-rich acidic matrix protein). We sequenced and identified 14 molluscan homologues of Dermatopontin from eight snail species belonging to the order Basommatophora and Stylommatophora. The bassommatophoran Dermatopontins fell into three types, one is suggested to be a shell matrix protein and the others are proteins having more general functions based on gene expression analyses. N-glycosylation is inferred to be important for the function involved in shell calcification, because potential N-glycosylation sites were found exclusively in the Dermatopontins considered as shell matrix proteins. The stylommatophoran Dermatopontins fell into two types, also suggested to comprise a shell matrix protein and a protein having a more general function. Phylogenetic analyses using maximum likelihood and Bayesian methods revealed that gene duplication events occurred independently in both basommatophoran and stylommatophoran lineages. These results suggest that the dermatopontin genes were co-opted for molluscan calcification at least twice independently after the divergence of basommatophoran and stylommatophoran lineages, or more recently than we have expected.

Bivalve-specific gene expansion in the pearl oyster genome: implications of adaptation to a sessile lifestyle

IntroductionBivalve molluscs have flourished in marine environments, and many species constitute important aquatic resources. Recently, whole genome sequences from two bivalves, the pearl oyster, Pinctada fucata, and the Pacific oyster, Crassostrea gigas, have been decoded, making it possible to compare genomic sequences among molluscs, and to explore general and lineage-specific genetic features and trends in bivalves. In order to improve the quality of sequence data for these purposes, we have updated the entire P. fucata genome assembly.ResultsWe present a new genome assembly of the pearl oyster, Pinctada fucata (version 2.0). To update the assembly, we conducted additional sequencing, obtaining accumulated sequence data amounting to 193× the P. fucata genome. Sequence redundancy in contigs that was caused by heterozygosity was removed in silico, which significantly improved subsequent scaffolding. Gene model version 2.0 was generated with the aid of manual gene annotations supplied by the P. fucata research community. Comparison of mollusc and other bilaterian genomes shows that gene arrangements of Hox, ParaHox, and Wnt clusters in the P. fucata genome are similar to those of other molluscs. Like the Pacific oyster, P. fucata possesses many genes involved in environmental responses and in immune defense. Phylogenetic analyses of heat shock protein70 and C1q domain-containing protein families indicate that extensive expansion of genes occurred independently in each lineage. Several gene duplication events prior to the split between the pearl oyster and the Pacific oyster are also evident. In addition, a number of tandem duplications of genes that encode shell matrix proteins are also well characterized in the P. fucata genome.ConclusionsBoth the Pinctada and Crassostrea lineages have expanded specific gene families in a lineage-specific manner. Frequent duplication of genes responsible for shell formation in the P. fucata genome explains the diversity of mollusc shell structures. These duplications reveal dynamic genome evolution to forge the complex physiology that enables bivalves to employ a sessile lifestyle in the intertidal zone.

Expression patterns of engrailed and dpp in the gastropod Lymnaea stagnalis.

We isolated the full-length cDNAs of engrailed and dpp-BMP2/4 orthologues from the pond snail Lymnaea stagnalis and examined their expression patterns during development by the whole mount in situ hybridization. At the gastrula and trochophore stages, engrailed is expressed in the peripheral ectoderm of the presumptive and invaginating shell gland, corroborating its role in the shell formation that is widely conserved among molluscs. At the same stages, dpp-BMP2/4 is expressed in the right-hand side ectoderm of the shell gland and in the invaginating stomodaeum. Unlike in the gastropod Patella vulgata, our results suggested that dpp-BMP2/4 has a role in the shell formation, rather than in the regional specification and that it could be involved in the specification pathway of the left–right asymmetry of the developing shell in L. stagnalis.

Novel Repetitive Structures, Deviant Protein-Encoding Sequences and Unidentified ORFs in the Mitochondrial Genome of the Brachiopod Lingula anatina

Complete sequence determination of the brachiopod Lingula anatina mtDNA (28,818 bp) revealed an organization that is remarkably atypical for an animal mt-genome. In addition to the usual set of 37 animal mitochondrial genes, which make up only 57% (16,555 bp) of the entire sequence, the genome contains lengthy unassigned sequences. All the genes are encoded in the same DNA strand, generally in a compact way, whereas the overall gene order is highly divergent in comparison with known animal mtDNA. Individual genes are generally longer and deviate considerably in sequence from their homologues in other animals. The genome contains two major repeat regions, in which 11 units of unassigned sequences and six genes (atp8, trnM, trnQ, trnV, and part of cox2 and nad2) are found in repetition, in the form of nested direct repeats of unparalleled complexity. One of the repeat regions contains unassigned repeat units dispersed among several unique sequences, novel repetitive structure for animal mtDNAs. Each of those unique sequences contains an open reading frame for a polypeptide between 80 and 357 amino acids long, potentially encoding a functional molecule, but none of them has been identified with known proteins. In both repeat regions, tRNA genes or tRNA gene-like sequences flank major repeated units, supporting the view that those structures play a role in the mitochondrial gene rearrangements. Although the intricate repeated organization of this genome can be explained by recurrent tandem duplications and subsequent deletions mediated by replication errors, other mechanisms, such as nonhomologous recombinations, appear to explain certain structures more easily.