A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster.

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture.

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

Evolution of Epidermal Growth Factor (EGF)-like and Zona Pellucida Domains Containing Shell Matrix Proteins in Mollusks.

Several types of shell matrix proteins (SMPs) have been identified in molluskan shells. Their diversity is the consequence of various molecular processes, including domain shuffling and gene duplication. However, the evolutionary origin of most SMPs remains unclear. In this study, we investigated the evolutionary process EGF-like and zona pellucida (ZP) domains containing SMPs. Two types of the proteins (EGF-like protein (EGFL) and EGF-like and ZP domains containing protein (EGFZP)) were found in the pearl oyster, Pinctada fucata. In contrast, only EGFZP was identified in the gastropods. Phylogenetic analysis and genomic arrangement studies showed that EGFL and EGFZP formed a clade in bivalves, and their encoding genes were localized in tandem repeats on the same scaffold. In P. fucata, EGFL genes were expressed in the outer part of mantle epithelial cells are related to the calcitic shell formation. However, in both P. fucata and the limpet Nipponacmea fuscoviridis, EGFZP genes were expressed in the inner part of the mantle epithelial cells are related to aragonitic shell formation. Furthermore, our analysis showed that in P. fucata, the ZP domain interacts with eight SMPs that have various functions in the nacreous shell mineralization. The data suggest that the ZP domain can interact with other SMPs, and EGFL evolution in pterimorph bivalves represents an example of neo-functionalization that involves the acquisition of a novel protein through gene duplication.

Whole-Genome Sequencing Highlights Conservative Genomic Strategies of a Stress-Tolerant, Long-Lived Scleractinian Coral, Porites australiensis Vaughan, 1918.

Massive corals of the genus Porites, common, keystone reef builders in the Indo-Pacific Ocean, are distinguished by their relative stress tolerance and longevity. In order to identify genetic bases of these attributes, we sequenced the complete genome of a massive coral, Porites australiensis. We developed a genome assembly and gene models of comparable quality to those of other coral genomes. Proteome analysis identified 60 Porites skeletal matrix protein genes, all of which show significant similarities to genes from other corals and even to those from a sea anemone, which has no skeleton. Nonetheless, 30% of its skeletal matrix proteins were unique to Porites and were not present in the skeletons of other corals. Comparative genomic analyses showed that genes widely conserved among other organisms are selectively expanded in Porites. Specifically, comparisons of transcriptomic responses of P. australiensis and Acropora digitifera, a stress-sensitive coral, reveal significant differences in regard to genes that respond to increased water temperature, and some of the genes expanded exclusively in Porites may account for the different thermal tolerances of these corals. Taken together, widely shared genes may have given rise to unique biological characteristics of Porites, massive skeletons and stress tolerance.

The 'Shellome' of the Crocus Clam Tridacna crocea Emphasizes Essential Components of Mollusk Shell Biomineralization.

Molluscan shells are among the most fascinating research objects because of their diverse morphologies and textures. The formation of these delicate biomineralized structures is a matrix-mediated process. A question that arises is what are the essential components required to build these exoskeletons. In order to understand the molecular mechanisms of molluscan shell formation, it is crucial to identify organic macromolecules in different shells from diverse taxa. In the case of bivalves, however, taxon sampling in previous shell proteomics studies are focused predominantly on representatives of the class Pteriomorphia such as pearl oysters, edible oysters and mussels. In this study, we have characterized the shell organic matrix from the crocus clam, Tridacna crocea, (Heterodonta) using various biochemical techniques, including SDS-PAGE, FT-IR, monosaccharide analysis, and enzyme-linked lectin assay (ELLA). Furthermore, we have identified a number of shell matrix proteins (SMPs) using a comprehensive proteomics approach combined to RNA-seq. The biochemical studies confirmed the presence of proteins, polysaccharides, and sulfates in the T. crocea shell organic matrix. Proteomics analysis revealed that the majority of the T. crocea SMPs are novel and dissimilar to known SMPs identified from the other bivalve species. Meanwhile, the SMP repertoire of the crocus clam also includes proteins with conserved functional domains such as chitin-binding domain, VWA domain, and protease inhibitor domain. We also identified BMSP (Blue Mussel Shell Protein, originally reported from Mytilus), which is widely distributed among molluscan shell matrix proteins. Tridacna SMPs also include low-complexity regions (LCRs) that are absent in the other molluscan genomes, indicating that these genes may have evolved in specific lineage. These results highlight the diversity of the organic molecules - in particular proteins - that are essential for molluscan shell formation.

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.

Functional shell matrix proteins tentatively identified by asymmetric snail shell morphology.

Molluscan shell matrix proteins (SMPs) are essential in biomineralization. Here, we identify potentially important SMPs by exploiting the asymmetric shell growth in snail, Lymnaea stagnalis. Asymmetric shells require bilaterally asymmetric expression of SMP genes. We examined expression levels of 35,951 transcripts expressed in the left and right sides of mantle tissue of the pond snail, Lymnaea stagnalis. This transcriptome dataset was used to identify 207 SMPs by LC-MS/MS. 32 of the 207 SMP genes show asymmetric expression patterns, which were further verified for 4 of the 32 SMPs using quantitative PCR analysis. Among asymmetrically expressed SMPs in dextral snails, those that are more highly expressed on the left side than the right side are 3 times more abundant than those that are more highly expressed on the right than the left, suggesting potentially inhibitory roles of SMPs in shell formation. The 32 SMPs thus identified have distinctive features, such as conserved domains and low complexity regions, which may be essential in biomineralization.

Divergent northern and southern populations and demographic history of the pearl oyster in the western Pacific revealed with genomic SNPs.

In the open ocean without terrain boundaries, marine invertebrates with pelagic larvae can migrate long distances using ocean currents, suggesting reduced genetic diversification. Contrary to this assumption, however, genetic differentiation is often observed in marine invertebrates. In the present study, we sought to explain how population structure is established in the western Pacific Ocean, where the strong Kuroshio Current maintains high levels of gene flow from south to north, presumably promoting genetic homogeneity. We determined the population structure of the pearl oyster, Pinctada fucata, in the Indo-Pacific Ocean using genome-wide genotyping data from multiple sampling localities. Cluster analysis showed that the western Pacific population is distinct from that of the Indian Ocean, and that it is divided into northern (Japanese mainland) and southern (Nansei Islands, China, and Cambodia) populations. Genetic differentiation of P. fucata can be explained by geographic barriers in the Indian Ocean and a local lagoon, and by environmental gradients of sea surface temperature (SST) and oxygen concentration in the western Pacific. A genome scan showed evidence of adaptive evolution in genomic loci, possibly associated with changes in environmental factors, including SST and oxygen concentration. Furthermore, Bayesian simulation demonstrated that the past population expansion and division are congruent with ocean warming after the last glacial period. It is highly likely that the environmental gradient forms a genetic barrier that diversifies P. fucata populations in the western Pacific. This hypothesis helps to explain genetic differentiation and possible speciation of marine invertebrates.

Dual Gene Repertoires for Larval and Adult Shells Reveal Molecules Essential for Molluscan Shell Formation.

Molluscan shells, mainly composed of calcium carbonate, also contain organic components such as proteins and polysaccharides. Shell organic matrices construct frameworks of shell structures and regulate crystallization processes during shell formation. To date, a number of shell matrix proteins (SMPs) have been identified, and their functions in shell formation have been studied. However, previous studies focused only on SMPs extracted from adult shells, secreted after metamorphosis. Using proteomic analyses combined with genomic and transcriptomic analyses, we have identified 31 SMPs from larval shells of the pearl oyster, Pinctada fucata, and 111 from the Pacific oyster, Crassostrea gigas. Larval SMPs are almost entirely different from those of adults in both species. RNA-seq data also confirm that gene expression profiles for larval and adult shell formation are nearly completely different. Therefore, bivalves have two repertoires of SMP genes to construct larval and adult shells. Despite considerable differences in larval and adult SMPs, some functional domains are shared by both SMP repertoires. Conserved domains include von Willebrand factor type A (VWA), chitin-binding (CB), carbonic anhydrase (CA), and acidic domains. These conserved domains are thought to play crucial roles in shell formation. Furthermore, a comprehensive survey of animal genomes revealed that the CA and VWA-CB domain-containing protein families expanded in molluscs after their separation from other Lophotrochozoan linages such as the Brachiopoda. After gene expansion, some family members were co-opted for molluscan SMPs that may have triggered to develop mineralized shells from ancestral, nonmineralized chitinous exoskeletons.

Structural and functional analyses of calcium ion response factors in the mantle of Pinctada fucata.

The pearl oyster, Pinctada fucata, is cultured for pearl production in Japan. The shell of the pearl oyster consists of calcium carbonate and a small amount of organic matrix. Despite many studies of the shell matrix proteins, the mechanism by which calcium elements are transported from the mantle to the shell remains unclear. Investigating the molecular mechanism of calcium transportation, we prepared artificial seawater with a high concentration of calcium ions (10ASW) to induce calcification in the pearl oyster. When pearl oysters were cultured in 10ASW, unusual nanoparticles were precipitated on the surface of the nacreous layer. SDS-PAGE and 2D-PAGE analyses revealed that some calcium-sensing proteins (Sarcoplasmic Ca-binding Protein (Pf-SCP) and Pf-filamin A) might be related to the synthesis of these nanoparticles. The recombinant proteins of Pf-SCP can bind to calcium ions and accumulate nanoparticles of calcium carbonate crystals. However, transcriptomic analysis of the pearl oysters grown in 10ASW showed that the matrix protein genes in the shell did not differ before and after treatment with 10ASW. These results suggest that, despite increasing calcium transportation to the shell, treatment with a high concentration of calcium ions does not induce formation of the organic framework in the shell microstructure. These findings offer meaningful insights into the transportation of calcium elements from the mantle to the shell.

The habu genome reveals accelerated evolution of venom protein genes.

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Biochemical characterization of the skeletal matrix of the massive coral, Porites australiensis - The saccharide moieties and their localization.

To construct calcium carbonate skeletons of sophisticated architecture, scleractinian corals secrete an extracellular skeletal organic matrix (SOM) from aboral ectodermal cells. The SOM, which is composed of proteins, saccharides, and lipids, performs functions critical for skeleton formation. Even though polysaccharides constitute the major component of the SOM, its contribution to coral skeleton formation is poorly understood. To this end, we analyzed the SOM of the massive colonial coral, Porites australiensis, the skeleton of which has drawn great research interest because it records environmental conditions throughout the life of the colony. The coral skeleton was extensively cleaned, decalcified with acetic acid, and organic fractions were separated based on solubility. These fractions were analyzed using various techniques, including SDS-PAGE, FT-IR, in vitro crystallization, CHNS analysis, chromatography analysis of monosaccharide and enzyme-linked lectin assay (ELLA). We confirmed the acidic nature of SOM and the presence of sulphate, which is thought to initiate CaCO3 crystallization. In order to analyze glycan structures, we performed ELLA on the soluble SOM for the first time and found that it exhibits strong specificity to Datura stramonium lectin (DSL). Furthermore, using biotinylated DSL with anti-biotin antibody conjugated to nanogold, in situ localization of DSL-binding polysaccharides in the P. australiensis skeleton was performed. Signals were distributed on the surfaces of fiber-like crystals of the skeleton, suggesting that polysaccharides may modulate crystal shape. Our study emphasizes the importance of sugar moieties in biomineralization of scleractinian corals.

Two divergent Symbiodinium genomes reveal conservation of a gene cluster for sunscreen biosynthesis and recently lost genes.

BACKGROUND: The marine dinoflagellate, Symbiodinium, is a well-known photosynthetic partner for coral and other diverse, non-photosynthetic hosts in subtropical and tropical shallows, where it comprises an essential component of marine ecosystems. Using molecular phylogenetics, the genus Symbiodinium has been classified into nine major clades, A-I, and one of the reported differences among phenotypes is their capacity to synthesize mycosporine-like amino acids (MAAs), which absorb UV radiation. However, the genetic basis for this difference in synthetic capacity is unknown. To understand genetics underlying Symbiodinium diversity, we report two draft genomes, one from clade A, presumed to have been the earliest branching clade, and the other from clade C, in the terminal branch. RESULTS: The nuclear genome of Symbiodinium clade A (SymA) has more gene families than that of clade C, with larger numbers of organelle-related genes, including mitochondrial transcription terminal factor (mTERF) and Rubisco. While clade C (SymC) has fewer gene families, it displays specific expansions of repeat domain-containing genes, such as leucine-rich repeats (LRRs) and retrovirus-related dUTPases. Interestingly, the SymA genome encodes a gene cluster for MAA biosynthesis, potentially transferred from an endosymbiotic red alga (probably of bacterial origin), while SymC has completely lost these genes. CONCLUSIONS: Our analysis demonstrates that SymC appears to have evolved by losing gene families, such as the MAA biosynthesis gene cluster. In contrast to the conservation of genes related to photosynthetic ability, the terminal clade has suffered more gene family losses than other clades, suggesting a possible adaptation to symbiosis. Overall, this study implies that Symbiodinium ecology drives acquisition and loss of gene families.

A Spirochaete is suggested as the causative agent of Akoya oyster disease by metagenomic analysis.

Mass mortality that is acompanied by reddish browning of the soft tissues has been occurring in cultured pearl oyster, Pinctada fucata martensii. The disease is called Akoya oyster disease (AOD). Although spreading pattern of the disease and transmission experiments suggest that the disease is infectious, the causative agent has not yet been identified. We used shotgun and 16S rRNA-based metagenomic analysis to identify genes that are present specifically in affected oysters. The genes found only in diseased oysters were mostly bacterial origin, suggesting that the causative agent was a bacterial pathogen. This hypothesis was supported by the inhibition of AOD development in naïve oysters injected with the hemolymph of diseased animals followed immediately with penicillin bath-administration. Further analyses of the hemolymph and mantle specifically and universally detected genes of bacteria that belong to phylum Spirochaetes in diseased pearl oysters but not in healthy oysters. By in situ hybridization or immunostaining, a Brachyspira-like bacterium was observed in the smears of hemolymph from affected oysters, but not from healthy oysters. Phylogenetic analysis using 16S rRNA sequences showed that the presumptive causative bacterium was outside of but most closely related to family Brachyspiraceae. We propose 'Candidatus Maribrachyspira akoyae' gen. nov, sp nov., for this bacterium.

Metabolic and physiological interdependencies in the Bathymodiolus azoricus symbiosis.

The hydrothermal vent mussel Bathymodiolus azoricus lives in an intimate symbiosis with two types of chemosynthetic Gammaproteobacteria in its gills: a sulfur oxidizer and a methane oxidizer. Despite numerous investigations over the last decades, the degree of interdependence between the three symbiotic partners, their individual metabolic contributions, as well as the mechanism of carbon transfer from the symbionts to the host are poorly understood. We used a combination of proteomics and genomics to investigate the physiology and metabolism of the individual symbiotic partners. Our study revealed that key metabolic functions are most likely accomplished jointly by B. azoricus and its symbionts: (1) CO2 is pre-concentrated by the host for carbon fixation by the sulfur-oxidizing symbiont, and (2) the host replenishes essential biosynthetic TCA cycle intermediates for the sulfur-oxidizing symbiont. In return (3), the sulfur oxidizer may compensate for the host's putative deficiency in amino acid and cofactor biosynthesis. We also identified numerous 'symbiosis-specific' host proteins by comparing symbiont-containing and symbiont-free host tissues and symbiont fractions. These proteins included a large complement of host digestive enzymes in the gill that are likely involved in symbiont digestion and carbon transfer from the symbionts to the host.

Stepwise Evolution of Coral Biomineralization Revealed with Genome-Wide Proteomics and Transcriptomics.

Despite the importance of stony corals in many research fields related to global issues, such as marine ecology, climate change, paleoclimatogy, and metazoan evolution, very little is known about the evolutionary origin of coral skeleton formation. In order to investigate the evolution of coral biomineralization, we have identified skeletal organic matrix proteins (SOMPs) in the skeletal proteome of the scleractinian coral, Acropora digitifera, for which large genomic and transcriptomic datasets are available. Scrupulous gene annotation was conducted based on comparisons of functional domain structures among metazoans. We found that SOMPs include not only coral-specific proteins, but also protein families that are widely conserved among cnidarians and other metazoans. We also identified several conserved transmembrane proteins in the skeletal proteome. Gene expression analysis revealed that expression of these conserved genes continues throughout development. Therefore, these genes are involved not only skeleton formation, but also in basic cellular functions, such as cell-cell interaction and signaling. On the other hand, genes encoding coral-specific proteins, including extracellular matrix domain-containing proteins, galaxins, and acidic proteins, were prominently expressed in post-settlement stages, indicating their role in skeleton formation. Taken together, the process of coral skeleton formation is hypothesized as: 1) formation of initial extracellular matrix between epithelial cells and substrate, employing pre-existing transmembrane proteins; 2) additional extracellular matrix formation using novel proteins that have emerged by domain shuffling and rapid molecular evolution and; 3) calcification controlled by coral-specific SOMPs.

Organic matrices in metazoan calcium carbonate skeletons: Composition, functions, evolution.

Calcium carbonate skeletal tissues in metazoans comprise a small quantity of occluded organic macromolecules, mostly proteins and polysaccharides that constitute the skeletal matrix. Because its functions in modulating the biomineralization process are well known, the skeletal matrix has been extensively studied, successively via classical biochemical approaches, via molecular biology and, in recent years, via transcriptomics and proteomics. The optimistic view that the deposition of calcium carbonate minerals requires a limited number of macromolecules has been challenged, in the last decade, by high-throughput approaches. Such approaches have made possible the rapid identification of large sets of mineral-associated proteins, i.e., 'skeletal repertoires' or 'skeletomes', in several calcifying animal models, ranging from sponges to echinoderms. One of the consequences of this expanding set of data is that a simple definition of the skeletal matrix is no longer possible. This increase in available data, however, makes it easier to compare skeletal repertoires, shedding light on the fundamental evolutionary mechanisms affecting matrix components.

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

INTRODUCTION: Bivalve 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. RESULTS: We 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. CONCLUSIONS: Both 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.

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 Lingula's unique biomineralization. This Lingula genome provides resources for further studies of lophotrochozoan evolution.

Massive gene transfer and extensive RNA editing of a symbiotic dinoflagellate plastid genome.

Genome sequencing of Symbiodinium minutum revealed that 95 of 109 plastid-associated genes have been transferred to the nuclear genome and subsequently expanded by gene duplication. Only 14 genes remain in plastids and occur as DNA minicircles. Each minicircle (1.8-3.3 kb) contains one gene and a conserved noncoding region containing putative promoters and RNA-binding sites. Nine types of RNA editing, including a novel G/U type, were discovered in minicircle transcripts but not in genes transferred to the nucleus. In contrast to DNA editing sites in dinoflagellate mitochondria, which tend to be highly conserved across all taxa, editing sites employed in DNA minicircles are highly variable from species to species. Editing is crucial for core photosystem protein function. It restores evolutionarily conserved amino acids and increases peptidyl hydropathy. It also increases protein plasticity necessary to initiate photosystem complex assembly.