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.

Author Correction: Targeted mutagenesis of the ryanodine receptor by Platinum TALENs causes slow swimming behaviour in Pacific bluefin tuna (Thunnus orientalis).

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Targeted mutagenesis of the ryanodine receptor by Platinum TALENs causes slow swimming behaviour in Pacific bluefin tuna (Thunnus orientalis).

In bluefin tuna aquaculture, high mortalities of hatchery-reared juveniles occur in sea cages owing to wall collisions that are caused by high-speed swimming in panic due to changes in illuminance. Here, we report that targeted gene mutagenesis of the ryanodine receptor (RyR1b), which allows the sarcoplasmic reticulum to release Ca2+ in fast skeletal muscle, using highly active Platinum TALENs caused slow swimming behaviour in response to external stimuli in Pacific bluefin tuna (PBT) larvae. This characteristic would be a useful trait to prevent wall collisions in aquaculture production. A pair of Platinum TALENs targeting exons 2 and 43 of the PBT ryr1b gene induced deletions in each TALEN target site of the injected embryos with extremely high efficiency. In addition, ryr1b expression was significantly decreased in the mutated G0 larvae at 7 days after hatching (DAH). A touch-evoked escape behaviour assay revealed that the ryr1b-mutated PBT larvae swam away much less efficiently in response to mechanosensory stimulation at 7 DAH than did the wild-type larvae. Our results demonstrate that genome editing technologies are effective tools for determining the functional characterization of genes in a comparatively short period, and create avenues for facilitating genetic studies and breeding of bluefin tuna species.

A ddRAD-based genetic map and its integration with the genome assembly of Japanese eel (Anguilla japonica) provides insights into genome evolution after the teleost-specific genome duplication.

BACKGROUND: Recent advancements in next-generation sequencing technology have enabled cost-effective sequencing of whole or partial genomes, permitting the discovery and characterization of molecular polymorphisms. Double-digest restriction-site associated DNA sequencing (ddRAD-seq) is a powerful and inexpensive approach to developing numerous single nucleotide polymorphism (SNP) markers and constructing a high-density genetic map. To enrich genomic resources for Japanese eel (Anguilla japonica), we constructed a ddRAD-based genetic map using an Ion Torrent Personal Genome Machine and anchored scaffolds of the current genome assembly to 19 linkage groups of the Japanese eel. Furthermore, we compared the Japanese eel genome with genomes of model fishes to infer the history of genome evolution after the teleost-specific genome duplication. RESULTS: We generated the ddRAD-based linkage map of the Japanese eel, where the maps for female and male spanned 1748.8 cM and 1294.5 cM, respectively, and were arranged into 19 linkage groups. A total of 2,672 SNP markers and 115 Simple Sequence Repeat markers provide anchor points to 1,252 scaffolds covering 151 Mb (13%) of the current genome assembly of the Japanese eel. Comparisons among the Japanese eel, medaka, zebrafish and spotted gar genomes showed highly conserved synteny among teleosts and revealed part of the eight major chromosomal rearrangement events that occurred soon after the teleost-specific genome duplication. CONCLUSIONS: The ddRAD-seq approach combined with the Ion Torrent Personal Genome Machine sequencing allowed us to conduct efficient and flexible SNP genotyping. The integration of the genetic map and the assembled sequence provides a valuable resource for fine mapping and positional cloning of quantitative trait loci associated with economically important traits and for investigating comparative genomics of the Japanese eel.

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.

Reproduction-related genes in the pearl oyster genome.

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

Homozygous and heterozygous GH transgenesis alters fatty acid composition and content in the liver of Amago salmon (Oncorhynchus masou ishikawae).

Growth hormone (GH) transgenic Amago (Oncorhynchus masou ishikawae), containing the sockeye GH1 gene fused with metallothionein-B promoter from the same species, were generated and the physiological condition through lipid metabolism compared among homozygous (Tg/Tg) and heterozygous GH transgenic (Tg/+) Amago and the wild type control (+/+). Previously, we have reported that the adipose tissue was generally smaller in GH transgenic fish compared to the control, and that the Δ-6 fatty acyl desaturase gene was down-regulated in the Tg/+ fish. However, fatty acid (FA) compositions have not been measured previously in these fish. In this study we compared the FAs composition and content in the liver using gas chromatography. Eleven kinds of FA were detected. The composition of saturated and monounsaturated fatty acids (SFA and MUFA) such as myristic acid (14:0), palmitoleic acid (16:1n-7), and cis-vaccenic acid (cis-18:1n-7) was significantly (P<0.05) decreased in GH transgenic Amago. On the other hand, the composition of polyunsaturated fatty acids (PUFAs) such as linoleic acid (18:2n-6), arachidonic acid (20:4n-6), and docosapentaenoic acid (22:5n-3) was significantly (P<0.05) increased. Levels of serum glucose and triacylglycerol were significantly (P<0.05) decreased in the GH transgenics compared with +/+ fish. Furthermore, 3'-tag digital gene expression profiling was performed using liver tissues from Tg/Tg and +/+ fish, and showed that Mid1 interacting protein 1 (Mid1ip1), which is an important factor to activate Acetyl-CoA carboxylase (ACC), was down-regulated in Tg/Tg fish, while genes involved in FA catabolism were up-regulated, including long-chain-fatty-acid-CoA ligase 1 (ACSL1) and acyl-coenzyme A oxidase 3 (ACOX3). These data suggest that liver tissue from GH transgenic Amago showed starvation by alteration in glucose and lipid metabolism due to GH overexpression. The decrease of serum glucose suppressed Mid1ip1, and caused a decrease of de novo FA synthesis, resulting in a decrease of SFA and MUFA. This induced expression of ACSL1 and ACOX3 to produce energy through ß-oxidation in the GH transgenic Amago.

Identification of the hybrid between Oryzias latipes and Oryzias curvinotus using nuclear genes and mitochondrial gene region.

Commercialized transgenic fish should be sterilized to prevent their genetic effects on native wild fish in the event of their accidental release to the environment. Diploid and allotriploid hybrids between Oryzias latipes and Oryzias curvinotus are sterile, which contributes to the sterilization of transgenic O. latipes or O. curvinotus strains, for example, fluorescent aquarium fish. However, it is very difficult to distinguish these hybrids from their parental species by morphological measurements. Thus, we confirmed whether our previously developed species identification method for O. latipes and O. curvinotus by restriction fragment length polymorphism (RFLP) analysis of polymerase chain reaction (PCR) products of nuclear DNA gene regions, namely, PCR-RFLP using Hinf I or Hsp92 II for the aromatase gene, Rsa I for the calmodulin gene, and Hae III for the caspase-6 gene, is useful for distinguishing the diploid and allotriploid hybrids from their parental Oryzias species. We found that the hybrids have maternal mitochondrial DNA by PCR-RFLP analysis using Mse I of the 16S rRNA gene.

Genetic structure and polymorphisms of the N16 gene in Pinctada fucata.

The molluscan shell is a composite of inorganic crystals comprising calcium carbonate and a minute amount of organic matrix. The organic matrix (OM) is intimately involved in every step of shell formation and has consequently received much attention in recent years. However, most of the deposited information has resulted from cDNA analysis, with little analysis of the genome, including the presence and effects of polymorphic genes encoding OM proteins. The current study aimed to clarify the genome structure of the N16 gene from the Japanese pearl oyster, Pinctada fucata, with particular reference to polymorphisms. The N16 gene was analyzed using PCR and DNA sequencing. 23 polymorphic variants were identified from 28 individuals. The variants were analyzed for their relationship to shell formation. All the variations detected by genomic PCR appeared in the cDNAs, implying that all the polymorphisms were transcribed and translated into N16 proteins. Additional genome analysis revealed at least two N16 genes, which were sequentially positioned, each of them comprising four exons and three introns. Further analyses of the transcriptional regulation and function of the N16 genes may provide new insights into their role in molluscan biomineralization.

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.

Effects of growth hormone on the salmon pituitary proteome.

Growth hormone 1 (GH1), a pituitary hormone, plays a key role in the regulation of growth. Both excess GH1 treatment and overexpression of a GH1 transgene promote growth of salmon, but these animals exhibit physiological abnormalities in viability, fertility and metabolism, which might be related to pituitary function. However, the molecular dynamics induced in the pituitary by excess GH1 remain unknown. In this study, we performed iTRAQ proteome analysis of the amago salmon pituitary, with and without excess GH1 treatment, and found that the expression levels of proteins related to endocrine systems, metabolism, cell growth and proliferation were altered in the GH1-treated pituitary. Specifically, pituitary hormone prolactin (2.29 fold), and somatolactin α (0.14 fold) changed significantly. This result was confirmed by proteome and transcriptome analyses of pituitary from the GH1-transgenic (GH1-Tg) amago salmon. The dynamics of protein and gene expression in the pituitary of GH1-Tg amago salmon were similar to those of pituitary treated with excess GH1. Our findings suggest that not only excess GH1 hormone, but also the quantitative changes in other pituitary hormones, might be essential for the abnormal growth of amago salmon. These data will be useful in future attempts to increase the productivity of fish farming.

A second generation genetic linkage map of Japanese flounder (Paralichthys olivaceus).

BACKGROUND: Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine species in Northeast Asia. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. Commercial production of Japanese flounder could be increased by developing disease-resistant fish and improving commercially important traits. Previous maps have been constructed with AFLP markers and a limited number of microsatellite markers. In this study, improved genetic linkage maps are presented. In contrast with previous studies, these maps were built mainly with a large number of codominant markers so they can potentially be used to analyze different families and populations. RESULTS: Sex-specific genetic linkage maps were constructed for the Japanese flounder including a total of 1,375 markers [1,268 microsatellites, 105 single nucleotide polymorphisms (SNPs) and two genes]; 1,167 markers are linked to the male map and 1,067 markers are linked to the female map. The lengths of the male and female maps are 1,147.7 cM and 833.8 cM, respectively. Based on estimations of map lengths, the female and male maps covered 79 and 82% of the genome, respectively. Recombination ratio in the new maps revealed F:M of 1:0.7. All linkage groups in the maps presented large differences in the location of sex-specific recombination hot-spots. CONCLUSIONS: The improved genetic linkage maps are very useful for QTL analyses and marker-assisted selection (MAS) breeding programs for economically important traits in Japanese flounder. In addition, SNP flanking sequences were blasted against Tetraodon nigroviridis (puffer fish) and Danio rerio (zebrafish), and synteny analysis has been carried out. The ability to detect synteny among species or genera based on homology analysis of SNP flanking sequences may provide opportunities to complement initial QTL experiments with candidate gene approaches from homologous chromosomal locations identified in related model organisms.

Characterization of Japanese flounder karyotype by chromosome bandings and fluorescence in situ hybridization with DNA markers.

The chromosomes of Japanese flounder, Paralichthys olivaceus, were examined by conventional differential staining methods including G-, Q-, C-, silver (Ag)-, fluorochrome, and replication R-bandings and by fluorescence in situ hybridization (FISH) with 5S and 18S rDNAs and telomeric DNA as probes. Replication R-banding substantially made it possible to identify 24 homologous pairs by their RBG-banding pattern and relative length. Both rDNA loci were mapped to chromosome 1, where 5S and 18S rDNA loci were located at the centromeric region and secondary constriction, respectively. C-banding revealed that both rDNA loci were heterochromatic, and 18S rDNA loci were positive for chromomycin A(3) but negative for 4',6-diamidino-2-phenylindole (DAPI) staining. Telomeric FISH signals were observed at all chromosome ends and at the interstitial region of some chromosomes. The observed results were discussed in relation to the karyotype evolution in the order Pleuronectiformes.