Draft genome sequence and tissue expression panel of Pacific saury (Cololabis saira).

Pacific saury (Cololabis saira) is an important fish in several countries. Notably, the catch of this fish has markedly decreased recently, which might be due to environmental changes, including feeding habitat changes. However, no clear correlation has been observed. Therefore, the physiological basis of Pacific saury in relation to possible environmental factors must be understood. We sequenced the genome of Pacific saury and extracted RNA from nine tissues (brain, eye, gill, anterior/posterior guts, kidney, liver, muscle, and ovary). In 1.09 Gb assembled genome sequences, a total of 26,775 protein-coding genes were predicted, of which 26,241 genes were similar to known genes in a public database. Transcriptome analysis revealed that 24,254 genes were expressed in at least one of the nine tissues, and 7,495 were highly expressed in specific tissues. Based on the similarity of the expression profiles to those of model organisms, the transcriptome obtained was validated to reflect the characteristics of each tissue. Thus, the present genomic and transcriptomic data serve as useful resources for molecular studies on Pacific saury. In particular, we emphasize that the gene expression data, which serve as the tissue expression panel of this species, can be employed in comparative transcriptomics on marine environmental responses.

Transcriptome characterization of gonadal sex differentiation in Pacific bluefin tuna, Thunnus orientalis (Temminck et Schlegel).

Tunas (genus Thunnus) are one of the most ecologically and commercially important fish worldwide. To establish a biological basis for reproduction in this globally essential species, we have recently studied crucial reproductive aspects of the Pacific bluefin tuna (T. orientalis; PBT), as a model of tuna species, based on our closed-cycle aquaculture technology. In this study, we clarified the global expression profile of the genes regulating gonadal sex differentiation in PBT, as this developmental process is vital to sexual reproduction. Based on the results of our comparative (RNA-sequencing) and temporal (qRT-PCR) transcriptome analyses using the updated genome dataset, we propose the molecular mechanisms of gonadal sex differentiation in PBT. In female gonads, foxl2 and cyp19a1a (coding aromatase) are expressed at the onset of sex differentiation. Active aromatase-mediated estrogen biosynthesis, which includes positive regulation of cyp19a1a expression by Foxl2, induces ovarian differentiation. By contrast, dmrt1 and gsdf are upregulated in differentiating male gonads lacking active estrogen synthesis. Dmrt1 and Gsdf would mainly promote testicular differentiation. Furthermore, androgen biosynthesis is upregulated in differentiating male gonad. Endogenous androgens may also be vital to testicular differentiation. This study provides the first comprehensive data clarifying the molecular basis for gonadal sex differentiation in tunas.

Genome analysis of Parmales, the sister group of diatoms, reveals the evolutionary specialization of diatoms from phago-mixotrophs to photoautotrophs.

The order Parmales (class Bolidophyceae) is a minor group of pico-sized eukaryotic marine phytoplankton that contains species with cells surrounded by silica plates. Previous studies revealed that Parmales is a member of ochrophytes and sister to diatoms (phylum Bacillariophyta), the most successful phytoplankton group in the modern ocean. Therefore, parmalean genomes can serve as a reference to elucidate both the evolutionary events that differentiated these two lineages and the genomic basis for the ecological success of diatoms vs. the more cryptic lifestyle of parmaleans. Here, we compare the genomes of eight parmaleans and five diatoms to explore their physiological and evolutionary differences. Parmaleans are predicted to be phago-mixotrophs. By contrast, diatoms have lost genes related to phagocytosis, indicating the ecological specialization from phago-mixotrophy to photoautotrophy in their early evolution. Furthermore, diatoms show significant enrichment in gene sets involved in nutrient uptake and metabolism, including iron and silica, in comparison with parmaleans. Overall, our results suggest a strong evolutionary link between the loss of phago-mixotrophy and specialization to a silicified photoautotrophic life stage early in diatom evolution after diverging from the Parmales lineage.

Prediction of the Sex-Associated Genomic Region in Tunas (Thunnus Fishes).

Fish species have a variety of sex determination systems. Tunas (genus Thunnus) have an XY genetic sex determination system. However, the Y chromosome or responsible locus has not yet been identified in males. In a previous study, a female genome of Pacific bluefin tuna (T. orientalis) was sequenced, and candidates for sex-associated DNA polymorphisms were identified by a genome-wide association study using resequencing data. In the present study, we sequenced a male genome of Pacific bluefin tuna by long-read and linked-read sequencing technologies and explored male-specific loci through a comparison with the female genome. As a result, we found a unique region carrying the male-specific haplotype, where a homolog of estrogen sulfotransferase gene was predicted to be encoded. The genome-wide mapping of previously resequenced data indicated that, among the functionally annotated genes, only this gene, named sult1st6y, was paternally inherited in the males of Pacific bluefin tuna. We reviewed the RNA-seq data of southern bluefin tuna (T. maccoyii) in the public database and found that sult1st6y of southern bluefin tuna was expressed in all male testes, but absent or suppressed in the female ovary. Since estrogen sulfotransferase is responsible for the inactivation of estrogens, it is reasonable to assume that the expression of sult1st6y in gonad cells may inhibit female development, thereby inducing the individuals to become males. Thus, our results raise a promising hypothesis that sult1st6y is the sex determination gene in Thunnus fishes or at least functions at a crucial point in the sex-differentiation cascade.

Estimation of tuna population by the improved analytical pipeline of unique molecular identifier-assisted HaCeD-Seq (haplotype count from eDNA).

Many studies have investigated the ability to identify species from environmental DNA (eDNA). However, even when individual species are identified, the accurate estimation of their abundances by traditional eDNA analyses has been still difficult. We previously developed a novel analytical method called HaCeD-Seq (Haplotype Count from eDNA), which focuses on the mitochondrial D-loop sequence. The D-loop is a rapidly evolving sequence and has been used to estimate the abundance of eel species in breeding water. In the current study, we have further improved this method by applying unique molecular identifier (UMI) tags, which eliminate the PCR and sequencing errors and extend the detection range by an order of magnitude. Based on this improved HaCeD-Seq pipeline, we computed the abundance of Pacific bluefin tuna (Thunnus orientalis) in aquarium tanks at the Tokyo Sea Life Park (Kasai, Tokyo, Japan). This tuna species is commercially important but is at high risk of resource depletion. With the developed UMI tag method, 90 out of 96 haplotypes (94%) were successfully detected from Pacific bluefin tuna eDNA. By contrast, only 29 out of 96 haplotypes (30%) were detected when UMI tags were not used. Our findings indicate the potential for conducting non-invasive fish stock surveys by sampling eDNA.

Draft Whole-Genome Sequence of Triparma laevis f. inornata (Parmales, Bolidophyceae), Isolated from the Oyashio Region, Western North Pacific Ocean.

We present the first draft whole-genome sequence for the Parmales (Bolidophyceae, Heterokonta), a picoplanktonic sister group of diatoms, using a Triparma laevis f. inornata strain that was isolated from the Oyashio region in the western North Pacific Ocean.

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.

A novel jumbo Tenacibaculum maritimum lytic phage with head-fiber-like appendages.

A novel jumbo bacteriophage (myovirus) is described. The lytic phage of Tenacibaculum maritimum, which is the etiological agent of tenacibaculosis in a variety of farmed marine fish worldwide, was plaque-isolated from seawater around a fish aquaculture field in Japan. The phage had an isometric head 110-120 nm in diameter, from which several 50- to 100-nm-long flexible fiber-like appendages emanate, and a 150-nm-long rigid contractile tail. The full genomes of the two representative phages (PTm1 and PTm5) were 224,680 and 226,876 bp long, respectively, both with 29.7% GC content, and the number of predicted open reading frames (ORFs) was 308 and 306, respectively. The average nucleotide sequence identity between PTm1 and PTm5 was 99.95%, indicating they are quite similar to each other. A genetic relationship was found in 15.0-16.6% of the predicted ORFs among the T. maritimum phages PTm1 and PTm5, the Tenacibaculum spp. phage pT24, and the Sphingomonas paucimobilis phage PAU. Phylogenetic analysis based on the terminase large subunit genes revealed that these four phages (PTm1, PTm5, pT24 and PAU) are more closely related than the other 10 jumbo myoviruses that have similar genome sizes. Transmission electron microscopy observations suggest that the head fibers of the T. maritimum phage function as tentacles to search and recognize the host cell surface to facilitate infection.

Monitoring of the toxic dinoflagellate Alexandrium catenella in Osaka Bay, Japan using a massively parallel sequencing (MPS)-based technique.

Since 2002, blooms of Alexandrium catenella sensu Fraga et al. (2015) and paralytic shellfish toxicity events have occurred almost yearly in Osaka Bay, Japan. To better understand the triggers for reoccurring A. catenella blooms in Osaka Bay, phytoplankton community was monitored during the spring seasons of 2012-2015. Monitoring was performed using massively parallel sequencing (MPS)-based technique on amplicon sequences of the 18S rRNA gene. Dense blooms of A. catenella occurred every year except in 2012, however, there was no significant correlation with the environmental parameters investigated. Plankton community diversity decreased before and middle of the A. catenella blooms, suggesting that the decline in diversity could be an indicator for the bloom occurrence. The yearly abundance pattern of A. catenella cells obtained by morphology-based counting coincided with the relative sequence abundances, which supports the effectiveness of MPS-based phytoplankton monitoring.

Molecular diet analysis of Anguilliformes leptocephalus larvae collected in the western North Pacific.

Natural diets of leptocephalus larvae have been enigmatic. In this study, we collected DNA samples from the gut contents and body surface of leptocephali belonging to the five Anguilliform families (Anguillidae, Chlopsidae, Congridae, Muraenidae, and Serrivomeridae) from the northwest Pacific and performed next-generation 18S rDNA sequencing. Wide variety of eukaryotes was detected in both samples, from which eight eukaryotic groups (jellyfish, conoid parasite, tunicate, copepod, krill, segmented worm, fungi, and dinoflagellate) were selected on the basis of abundance. All groups except conoid parasites were common in both the samples. Cnidarian 18S rDNA reads were the most abundant in both the samples; however, the number of samples having cnidarian reads and the read counts were significantly higher in the body surface scraping samples than in the gut content samples, regardless of careful rinsing of the body surface. These results indicate that the cnidarian DNAs are most likely found because of cross contamination from the body surface and/or environment. 18S rDNA read counts of copepod and tunicate in the gut contents were greater than or comparable with those in the body surface scraping samples, which may correspond to the previous observations of fecal pellets and larvacean houses in the leptocephali gut. Thus, the present study supports previous implications that leptocephali utilize detritus materials, so called marine snow.

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.

Prediction of Horizontally and Widely Transferred Genes in Prokaryotes.

Horizontal gene transfer (HGT) is the process whereby an organism acquires exogenous genes (horizontally transferred genes or HT genes) that are not inherited from the parent, but are derived from another organism. In prokaryotes, HGT has been considered as one of the important driving forces of evolution. Previously, genome-wide analyses have been conducted for estimating the proportion of HT genes in prokaryotic genomes, but the number of species examined at the time was limited, and gene annotation was relatively poor. Currently, tens of thousands of prokaryotic genomes have been published and gene annotation resources have improved. In the present study, HT gene prediction method was modified so that the estimate was robust to gene length, conducting a comprehensive search using 3017 representative prokaryotic genomes belonging to 1348 species. The result showed that an average of 13% (ranging from 0% to 30% across species) of protein-coding genes was predicted as being of horizontal origin. The proportion of the predicted HT genes per species was associated with the species' habitat, while a positive correlation between the proportion and genomic nucleotide frequency was also observed. Moreover, the functions of the predicted HT genes were inferred and compared according to two popular databases, the Clusters of Orthologous Groups and the Kyoto Encyclopedia of Genes and Genomes. As a result, both databases indicated that many of the widely transferred genes were involved in mobile genetic elements (transposons, phages, and plasmids) as expected. Notably, the present study predicted that six as-yet-uncharacterized genes were widely distributed HT genes, and therefore, will be interesting targets for evolutionary studies. Thus, this study demonstrates that a data-driven approach using massive sequence data may contribute to a broader understanding of HGT in prokaryotes.

The yellowtail (Seriola quinqueradiata) genome and transcriptome atlas of the digestive tract.

Seriola quinqueradiata (yellowtail) is the most widely farmed and economically important fish in aquaculture in Japan. In this study, we used the genome of haploid yellowtail fish larvae for de novo assembly of whole-genome sequences, and built a high-quality draft genome for the yellowtail. The total length of the assembled sequences was 627.3 Mb, consisting of 1,394 scaffold sequences (>2 kb) with an N50 length of 1.43 Mb. A total of 27,693 protein-coding genes were predicted for the draft genome, and among these, 25,832 predicted genes (93.3%) were functionally annotated. Given our lack of knowledge of the yellowtail digestive system, and using the annotated draft genome as a reference, we conducted an RNA-Seq analysis of its three digestive organs (stomach, intestine and rectum). The RNA-Seq results highlighted the importance of certain genes in encoding proteolytic enzymes necessary for digestion and absorption in the yellowtail gastrointestinal tract, and this finding will accelerate development of formulated feeds for this species. Since this study offers comprehensive annotation of predicted protein-coding genes, it has potential broad application to our understanding of yellowtail biology and aquaculture.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication.

BACKGROUND: Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods. RESULTS: In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions. CONCLUSIONS: Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

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.

Analysis of the complete genome sequence of Nocardia seriolae UTF1, the causative agent of fish nocardiosis: The first reference genome sequence of the fish pathogenic Nocardia species.

Nocardiosis caused by Nocardia seriolae is one of the major threats in the aquaculture of Seriola species (yellowtail; S. quinqueradiata, amberjack; S. dumerili and kingfish; S. lalandi) in Japan. Here, we report the complete nucleotide genome sequence of N. seriolae UTF1, isolated from a cultured yellowtail. The genome is a circular chromosome of 8,121,733 bp with a G+C content of 68.1% that encodes 7,697 predicted proteins. In the N. seriolae UTF1 predicted genes, we found orthologs of virulence factors of pathogenic mycobacteria and human clinical Nocardia isolates involved in host cell invasion, modulation of phagocyte function and survival inside the macrophages. The virulence factor candidates provide an essential basis for understanding their pathogenic mechanisms at the molecular level by the fish nocardiosis research community in future studies. We also found many potential antibiotic resistance genes on the N. seriolae UTF1 chromosome. Comparative analysis with the four existing complete genomes, N. farcinica IFM 10152, N. brasiliensis HUJEG-1 and N. cyriacigeorgica GUH-2 and N. nova SH22a, revealed that 2,745 orthologous genes were present in all five Nocardia genomes (core genes) and 1,982 genes were unique to N. seriolae UTF1. In particular, the N. seriolae UTF1 genome contains a greater number of mobile elements and genes of unknown function that comprise the differences in structure and gene content from the other Nocardia genomes. In addition, a lot of the N. seriolae UTF1-specific genes were assigned to the ABC transport system. Because of limited resources in ocean environments, these N. seriolae UTF1 specific ABC transporters might facilitate adaptation strategies essential for marine environment survival. Thus, the availability of the complete N. seriolae UTF1 genome sequence will provide a valuable resource for comparative genomic studies of N. seriolae isolates, as well as provide new insights into the ecological and functional diversity of the genus Nocardia.

Complete Genome Sequence of Ichthyobacterium seriolicida JBKA-6T, Isolated from Yellowtail (Seriola quinqueradiata) Affected by Bacterial Hemolytic Jaundice.

Ichthyobacterium seriolicida is a fish bacterial pathogen that causes hemolytic jaundice in farmed yellowtail in Japan. To understand more about the characteristics of this bacterium, we determined its complete genome sequence. Two hemolysin genes which may be important for its pathogenicity were identified in the I. seriolicida genome.

Full-Genome Sequencing and Confirmation of the Causative Agent of Erythrocytic Inclusion Body Syndrome in Coho Salmon Identifies a New Type of Piscine Orthoreovirus.

Erythrocytic inclusion body syndrome (EIBS) causes mass mortality in farmed salmonid fish, including the coho salmon, Onchorhynchus kisutchi, and chinook salmon, O. tshawytscha. The causative agent of the disease is a virus with an icosahedral virion structure, but this virus has not been characterized at the molecular level. In this study, we sequenced the genome of a virus purified from EIBS-affected coho salmon. The virus has 10 dsRNA genomic segments (L1, L2, L3, M1, M2, M3, S1, S2, S3, and S4), which closely resembles the genomic organization of piscine orthoreovirus (PRV), the causative agent of heart and skeletal inflammation (HSMI) in Atlantic salmon and HSMI-like disease in coho salmon. The genomic segments of the novel virus contain at least 10 open reading frames (ORFs): lambda 1 (λ1), λ2, λ3, mu 1 (µ1), µ2, µNS, sigma 1 (σ1), σ2, σ3, and σNS. An additional ORF encoding a 12.6-kDa protein (homologue of PRV p13) occurs in the same genomic segment as σ3. Phylogenetic analyses based on S1 and λ3 suggest that this novel virus is closely related to PRV, but distinctly different. Therefore, we designated the new virus 'piscine orthoreovirus 2' (PRV-2). Reverse transcription-quantitative real-time PCR revealed a significant increase in PRV-2 RNA in fish blood after the artificial infection of EIBS-naïve fish but not in that of fish that had recovered from EIBS. The degree of anemia in each fish increased as the PRV-2 RNA increased during an epizootic season of EIBS on an inland coho salmon farm. These results indicate that PRV-2 is the probable causative agent of EIBS in coho salmon, and that the host acquires immunity to reinfection with this virus. Further research is required to determine the host range of PRV species and the relationship between EIBS and HSMI in salmonid fish.

Complete Genome Sequence of Nonagglutinating Lactococcus garvieae Strain 122061 Isolated from Yellowtail in Japan.

Nonagglutinating Lactococcus garvieae has been isolated from diseased farmed yellowtail in Japan since 2012. In this study, the complete genome and plasmid sequence of nonagglutinating L. garvieae strain 122061 was determined, to our knowledge, for the first time.

CHARACTERIZATION OF 17 NEW MICROSATELLITE MARKERS FOR THE DINOFLAGELLATE ALEXANDRIUM FUNDYENSE (DINOPHYCEAE), A HARMFUL ALGAL BLOOM SPECIES.

Alexandrium fundyense is the toxic marine dinoflagellate responsible for "red tide" events in temperate and sub-arctic waters worldwide. In the Gulf of Maine (GOM) and Bay of Fundy in the Northwest Atlantic, blooms of A. fundyense recur annually, and are associated with major health and ecosystem impacts. In this region, microsatellite markers have been used to investigate genetic structure and gene flow; however, the loci currently available for this species were isolated from populations from Japan and the North Sea, and only a subset are suitable for the analysis of A. fundyense populations in the Northwest Atlantic. To facilitate future studies of A. fundyense blooms, both in this region and globally, we isolated and characterized 17 polymorphic microsatellite loci from 31 isolates collected from the GOM and from the Nauset Marsh System, an estuary on Cape Cod, MA, USA. These loci yielded between two and 15 alleles per locus, with an average of 7.1. Gene diversities ranged from 0.297 to 0.952. We then analyzed these same 31 isolates using previously published markers for comparison. We determined the new markers are sufficiently variable and better suited for the investigation of genetic structure, bloom dynamics, and diversity in the Northwest Atlantic.