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.

Linking Pedigree Information to the Gene Expression Phenotype to Understand Differential Family Survival Mechanisms in Highly Fecund Fish: A Case Study in the Larviculture of Pacific Bluefin Tuna.

Mass spawning in fish culture often brings about a marked variance in family size, which can cause a reduction in effective population sizes in seed production for stock enhancement. This study reports an example of combined pedigree information and gene expression phenotypes to understand differential family survival mechanisms in early stages of Pacific bluefin tuna, Thunnus orientalis, in a mass culture tank. Initially, parentage was determined using the partial mitochondrial DNA control region sequence and 11 microsatellite loci at 1, 10, 15, and 40 days post-hatch (DPH). A dramatic proportional change in the families was observed at around 15 DPH; therefore, transcriptome analysis was conducted for the 15 DPH larvae using a previously developed oligonucleotide microarray. This analysis successfully addressed the family-specific gene expression phenotypes with 5739 differentially expressed genes and highlighted the importance of expression levels of gastric-function-related genes at the developmental stage for subsequent survival. This strategy demonstrated herein can be broadly applicable to species of interest in aquaculture to comprehend the molecular mechanism of parental effects on offspring survival, which will contribute to the optimization of breeding technologies.

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 Genome Sequence of a Novel Picornavirus Isolated from Japanese Eel (Anguilla japonica).

We report the draft genome sequence of a novel member of the order Picornavirales that was obtained from the gills of farmed Japanese eel (Anguilla japonica). A putative polyprotein encoded by the genome was similar to that of other picornaviruses and shared 31% amino acid identity with that of eel picornavirus 1.

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.

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.

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.

Antibody profiling using a recombinant protein-based multiplex ELISA array accelerates recombinant vaccine development: Case study on red sea bream iridovirus as a reverse vaccinology model.

Predicting antigens that would be protective is crucial for the development of recombinant vaccine using genome based vaccine development, also known as reverse vaccinology. High-throughput antigen screening is effective for identifying vaccine target genes, particularly for pathogens for which minimal antigenicity data exist. Using red sea bream iridovirus (RSIV) as a research model, we developed enzyme-linked immune sorbent assay (ELISA) based RSIV-derived 72 recombinant antigen array to profile antiviral antibody responses in convalescent Japanese amberjack (Seriola quinqueradiata). Two and three genes for which the products were unrecognized and recognized, respectively, by antibodies in convalescent serum were selected for recombinant vaccine preparation, and the protective effect was examined in infection tests using Japanese amberjack and greater amberjack (S. dumerili). No protection was provided by vaccines prepared from gene products unrecognized by convalescent serum antibodies. By contrast, two vaccines prepared from gene products recognized by serum antibodies induced protective immunity in both fish species. These results indicate that ELISA array screening is effective for identifying antigens that induce protective immune responses. As this method does not require culturing of pathogens, it is also suitable for identifying protective antigens to un-culturable etiologic agents.

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.

V-GAP: Viral genome assembly pipeline.

Next-generation sequencing technologies have allowed the rapid determination of the complete genomes of many organisms. Although shotgun sequences from large genome organisms are still difficult to reconstruct perfect contigs each of which represents a full chromosome, those from small genomes have been assembled successfully into a very small number of contigs. In this study, we show that shotgun reads from phage genomes can be reconstructed into a single contig by controlling the number of read sequences used in de novo assembly. We have developed a pipeline to assemble small viral genomes with good reliability using a resampling method from shotgun data. This pipeline, named V-GAP (Viral Genome Assembly Pipeline), will contribute to the rapid genome typing of viruses, which are highly divergent, and thus will meet the increasing need for viral genome comparisons in metagenomic studies.

Influences of diurnal sampling bias on fixed-point monitoring of plankton biodiversity determined using a massively parallel sequencing-based technique.

In this study, we investigated the influence of diurnal sampling bias on the community structure of plankton by comparing the biodiversity among seawater samples (n=9) obtained every 3h for 24h by using massively parallel sequencing (MPS)-based plankton monitoring at a fixed point conducted at Himedo seaport in Yatsushiro Sea, Japan. The number of raw operational taxonomy units (OTUs) and OTUs after re-sampling was 507-658 (558 ± 104, mean ± standard deviation) and 448-544 (467 ± 81), respectively, indicating high plankton biodiversity at the sampling location. The relative abundance of the top 20 OTUs in the samples from Himedo seaport was 48.8-67.7% (58.0 ± 5.8%), and the highest-ranked OTU was Pseudo-nitzschia species (Bacillariophyta) with a relative abundance of 17.3-39.2%, followed by Oithona sp. 1 and Oithona sp. 2 (Arthropoda). During seawater sampling, the semidiurnal tidal current having an amplitude of 0.3ms(-1) was dominant, and the westward residual current driven by the northeasterly wind was continuously observed during the 24-h monitoring. Therefore, the relative abundance of plankton species apparently fluctuated among the samples, but no significant difference was noted according to G-test (p>0.05). Significant differences were observed between the samples obtained from a different locality (Kusuura in Yatsushiro Sea) and at different dates, suggesting that the influence of diurnal sampling bias on plankton diversity, determined using the MPS-based survey, was not significant and acceptable.