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.

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.

Evolutionary changes of multiple visual pigment genes in the complete genome of Pacific bluefin tuna.

Tunas are migratory fishes in offshore habitats and top predators with unique features. Despite their ecological importance and high market values, the open-ocean lifestyle of tuna, in which effective sensing systems such as color vision are required for capture of prey, has been poorly understood. To elucidate the genetic and evolutionary basis of optic adaptation of tuna, we determined the genome sequence of the Pacific bluefin tuna (Thunnus orientalis), using next-generation sequencing technology. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified five common fish visual pigment genes: red-sensitive (middle/long-wavelength sensitive; M/LWS), UV-sensitive (short-wavelength sensitive 1; SWS1), blue-sensitive (SWS2), rhodopsin (RH1), and green-sensitive (RH2) opsin genes. Sequence comparison revealed that tuna's RH1 gene has an amino acid substitution that causes a short-wave shift in the absorption spectrum (i.e., blue shift). Pacific bluefin tuna has at least five RH2 paralogs, the most among studied fishes; four of the proteins encoded may be tuned to blue light at the amino acid level. Moreover, phylogenetic analysis suggested that gene conversions have occurred in each of the SWS2 and RH2 loci in a short period. Thus, Pacific bluefin tuna has undergone evolutionary changes in three genes (RH1, RH2, and SWS2), which may have contributed to detecting blue-green contrast and measuring the distance to prey in the blue-pelagic ocean. These findings provide basic information on behavioral traits of predatory fish and, thereby, could help to improve the technology to culture such fish in captivity for resource management.

Complete genome sequence analysis of two Pseudomonas plecoglossicida phages, potential therapeutic agents.

Pseudomonas plecoglossicida is a lethal pathogen of ayu (Plecoglossus altivelis) in Japan and is responsible for substantial economic costs to ayu culture. Previously, we demonstrated the efficacy of phage therapy against P. plecoglossicida infection using two lytic phages (PPpW-3 and PPpW-4) (S. C. Park, I. Shimamura, M. Fukunaga, K. Mori, and T. Nakai, Appl Environ Microbiol 66:1416-1422, 2000, http://dx.doi.org/10.1128/AEM.66.4.1416-1422.2000; S. C. Park and T. Nakai, Dis Aquat Org 53:33-39, 2003, http://dx.doi.org/10.3354/dao053033). In the present study, the complete genome sequences of these therapeutic P. plecoglossicida phages were determined and analyzed for deleterious factors as therapeutic agents. The genome of PPpW-3 (myovirus) consisted of 43,564 bp with a GC content of 61.1% and 66 predicted open reading frames (ORFs). Approximately half of the genes were similar to the genes of the Escherichia coli phage vB_EcoM_ECO1230-10 (myovirus). The genome of PPpW-4 (podovirus) consisted of 41,386 bp with a GC content of 56.8% and 50 predicted ORFs. More than 70% of the genes were similar to the genes of Pseudomonas fluorescens phage ϕIBB-PF7A and Pseudomonas putida phage ϕ15 (podoviruses). The whole-genome analysis revealed that no known virulence genes were present in PPpW-3 and PPpW-4. An integrase gene was found in PPpW-3, but other factors used for lysogeny were not confirmed. The PCR detection of phage genes in phage-resistant variants provided no evidence of lysogenic activity in PPpW-3 and PPpW-4. We conclude that these two lytic phages qualify as therapeutic agents.

Full-genome sequence of a novel myovirus, GF-2, infecting Edwardsiella tarda: comparison with other Edwardsiella myoviral genomes.

Edwardsiellosis, which is caused by Edwardsiella tarda, a Gram-negative bacterium, is one of the most serious infectious diseases in both marine and freshwater fish farms worldwide. Previously, we reported the complete genome sequences of three E. tarda-lytic bacteriophages (two podoviruses and a myovirus), which were isolated from fish tissues and fish-rearing seawater. Further genomic information regarding E. tarda phages is important for understanding phage-host interactions as well as for applications of the phages for the control of disease. Here, we report the complete genome sequence of a novel E. tarda phage (GF-2) of myovirus morphology (family Myoviridae), isolated from tissue homogenates of a cultured Japanese flounder (Paralichthys olivaceus) that succumbed to edwardsiellosis in Japan. The size of the entire genome was 43,129 bp, with a GC content of 51.3 % and containing 82 open reading frames (ORFs). The GF-2 genome possesses lysogeny-related genes that have not been found in the reported Edwardsiella phage genomes. Comparative genomics of Edwardsiella myophages suggest that the C-terminal domains of the tail fiber proteins have relevance to their host specificity. Thus, GF-2 genome information provides a novel resource for our understanding of the molecular mechanisms involved in their host specificity and for detection of E. tarda in aquaculture environments.

Comparative transcriptomics of Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis.

BACKGROUND: Salmon species vary in susceptibility to infections with the salmon louse (Lepeophtheirus salmonis). Comparing mechanisms underlying responses in susceptible and resistant species is important for estimating impacts of infections on wild salmon, selective breeding of farmed salmon, and expanding our knowledge of fish immune responses to ectoparasites. Herein we report three L. salmonis experimental infection trials of co-habited Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha, profiling hematocrit, blood cortisol concentrations, and transcriptomic responses of the anterior kidney and skin to the infection. RESULTS: In all trials, infection densities (lice per host weight (g)) were consistently highest on chum salmon, followed by Atlantic salmon, and lowest in pink salmon. At 43 days post-exposure, all lice had developed to motile stages, and infection density was uniformly low among species. Hematocrit was reduced in infected Atlantic and chum salmon, and cortisol was elevated in infected chum salmon. Systemic transcriptomic responses were profiled in all species and large differences in response functions were identified between Atlantic and Pacific (chum and pink) salmon. Pink and chum salmon up-regulated acute phase response genes, including complement and coagulation components, and down-regulated antiviral immune genes. The pink salmon response involved the largest and most diverse iron sequestration and homeostasis mechanisms. Pattern recognition receptors were up-regulated in all species but the active components were often species-specific. C-type lectin domain family 4 member M and acidic mammalian chitinase were specifically up-regulated in the resistant pink salmon. CONCLUSIONS: Experimental exposures consistently indicated increased susceptibility in chum and Atlantic salmon, and resistance in pink salmon, with differences in infection density occurring within the first three days of infection. Transcriptomic analysis suggested candidate resistance functions including local inflammation with cytokines, specific innate pattern recognition receptors, and iron homeostasis. Suppressed antiviral immunity in both susceptible and resistant species indicates the importance of future work investigating co-infections of viral pathogens and lice.

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.

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.

Comparative genomics reveals that a fish pathogenic bacterium Edwardsiella tarda has acquired the locus of enterocyte effacement (LEE) through horizontal gene transfer.

BACKGROUND: Edwardsiella tarda is an enterobacterium which causes edwardsiellosis, a fatal disease of cultured fishes such as red sea bream, eel, and flounder. Preventing the occurrence of E. tarda infection has thus been an important issue in aquaculture. E. tarda has been isolated from other animals and from many environments; however, the relationship between the genotype and evolutionary process of this pathogen is not fully understood. To clarify this relationship, we sequenced and compared the genomes of pathogenic and non-pathogenic E. tarda strains isolated from fish, human, and eel pond using next-generation sequencing technology. RESULTS: Eight strains of E. tarda were sequenced with high accuracy (>99.9%) with coverages from 50- to 400-fold. The obtained reads were mapped to a public reference genome. By comparing single nucleotide and insertion/deletion polymorphisms, we found that an attenuated strain of E. tarda had a loss-of-function mutation in a gene related to the type III secretion system (T3SS), suggesting that this gene is involved in the virulence of E. tarda. A comprehensive gene comparison indicated that fish pathogenic strains possessed a type VI secretion system (T6SS) and pilus assembly genes in addition to the T3SS. Moreover, we found that an E. tarda strain isolated from red sea bream harbored two pathogenicity islands of T3SS and T6SS, which were absent in other strains. In particular, this T3SS was homologous to the locus of enterocyte effacement (LEE) in enteropathogenic and enterohemorrhagic Escherichia coli. Evolutionary analysis suggested that this locus, here named Et-LEE (E. tarda LEE), was introgressed into the E. tarda genome through horizontal transfer. CONCLUSIONS: We found significant differences in the presence/absence of virulence-related genes among E. tarda strains, reflecting their evolutionary relationship. In particular, a single genotype previously proposed for fish-pathogenic strains may be further divided into two subgroups. Furthermore, the current study demonstrated, for the first time, that a fish pathogenic bacterium carried a LEE-like pathogenicity island which was previously reported only in zoonotic pathogenic enterobacteria. These findings will contribute to the exploration of strain-specific drug targets against E. tarda in aquafarms, while also shedding light on the evolution of pathogenesis in enterobacteria.

Tetra-repeat microsatellite markers for the masu salmon (Oncorhynchus masou masou) and its application in cross-subspecies amplification.

We developed tetranucleotide-repeat microsatellite markers for the masu salmon (Oncorhynchus masou) complex. 454 pyrosequencing was used to discover repeat motifs, and seven polymorphic microsatellite-primer sets were identified. The number of alleles detected at each locus ranged from four to 24 and the expected heterozygosity varied from 0.57 to 0.92. Cross-subspecies amplification for O. m. masou, O. m. ishikawae and O. m. subsp. was successful. These microsatellites can be utilized in studies of genetic structure, genetic diversity, and intra- and inter-subspecific hybridization, making a contribution to conservation and management of the Oncorhynchus masou complex.

Transcriptomics of coping strategies in free-swimming Lepeophtheirus salmonis (Copepoda) larvae responding to abiotic stress.

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite of wild and farmed salmon in the Northern Hemisphere. Infections of farmed salmon are of economic and ecological concern. Nauplius and copepodid salmon lice larvae are free-swimming and disperse in the water column until they encounter a host. In this study, we characterized the sublethal stress responses of L. salmonis copepodid larvae by applying a 38K oligonucleotide microarray to profile transcriptomes following 24 h exposures to suboptimal salinity (30-10 parts per thousand (‰)) or temperature (16-4 °C) environments. Hyposalinity exposure resulted in large-scale gene expression changes relative to those elicited by a thermal gradient. Subsequently, transcriptome responses to a more finely resolved salinity gradient between 30 ‰ and 25 ‰ were profiled. Minimal changes occurred at 29 ‰ or 28 ‰, a threshold of response was identified at 27 ‰, and the largest response was at 25 ‰. Differentially expressed genes were clustered by pattern of expression, and clusters were characterized by functional enrichment analysis. Results indicate larval copepods adopt two distinct coping strategies in response to short-term hyposaline stress: a primary response using molecular chaperones and catabolic processes at 27 ‰; and a secondary response up-regulating ion pumps, transporters, a different suite of chaperones and apoptosis-related transcripts at 26 ‰ and 25 ‰. The results further our understanding of the tolerances of L. salmonis copepodids to salinity and temperature gradients and may assist in the development of salmon louse management strategies.

Regulation and expression of sexual differentiation factors in embryonic and extragonadal tissues of Atlantic salmon.

BACKGROUND: The products of cyp19, dax, foxl2, mis, sf1 and sox9 have each been associated with sex-determining processes among vertebrates. We provide evidence for expression of these regulators very early in salmonid development and in tissues outside of the hypothalamic-pituitary-adrenal/gonadal (HPAG) axis. Although the function of these factors in sexual differentiation have been defined, their roles in early development before sexual fate decisions and in tissues beyond the brain or gonad are essentially unknown. RESULTS: Bacterial artificial chromosomes containing salmon dax1 and dax2, foxl2b and mis were isolated and the regulatory regions that control their expression were characterized. Transposon integrations are implicated in the shaping of the dax and foxl2 loci. Splice variants for cyp19b1 and mis in both embryonic and adult tissues were detected and characterized. We found that cyp19b1 transcripts are generated that contain 5'-untranslated regions of different lengths due to cryptic splicing of the 3'-end of intron 1. We also demonstrate that salmon mis transcripts can encode prodomain products that present different C-termini and terminate before translation of the MIS hormone. Regulatory differences in the expression of two distinct aromatases cyp19a and cyp19b1 are exerted, despite transcription of their transactivators (ie; dax1, foxl2, sf1) occurring much earlier during embryonic development. CONCLUSIONS: We report the embryonic and extragonadal expression of dax, foxl2, mis and other differentiation factors that indicate that they have functions that are more general and not restricted to steroidogenesis and gonadogenesis. Spliced cyp19b1 and mis transcripts are generated that may provide regulatory controls for tissue- or development-specific activities. Selection of cyp19b1 transcripts may be regulated by DAX-1, FOXL2 and SF-1 complexes that bind motifs in intron 1, or by signals within exon 2 that recruit splicing factors, or both. The potential translation of proteins bearing only the N-terminal MIS prodomain may modulate the functions of other TGF ß family members in different tissues. The expression patterns of dax1 early in salmon embryogenesis implicate its role as a lineage determination factor. Other roles for these factors during embryogenesis and outside the HPAG axis are discussed.

Molecular cloning and expression analysis of interferon regulatory factor 10 (IRF10) in Japanese flounder, Paralichthys olivaceus.

Interferons (IFNs) are important for the defense against intracellular pathogens. Interferon regulatory factor 10 (IRF10) is a transcript factor involved in regulating IFN signaling induced by virus infections in birds, but little is know of its role in the immune response in non-avian vertebrates. Here, we identified and characterized the IRF10 gene and examined its expression pattern in a teleost fish, Japanese flounder, Paralichthys olivaceus. Japanese flounder IRF10 (PoIRF10) gene is a single copy gene, contains 8 exons and 7 introns and has 4918 nucleotides (nt) including an open reading frame (ORF) of 1212nt encoding 404 amino acids. The deduced PoIRF10 amino acid sequence contains a DNA-binding domain (DBD), nuclear localization signal (NLS) and IRF association domain (IAD). PoIRF10 is most closely related to chicken and zebrafish IRF10 by homology search and phylogenetic analysis. Putative binding sites for activator protein-1 (AP-1), CAAT-enhancer binding protein (C/EBP), C/EBPß, cAMP response element-binding protein (CRE-BP), IFN-stimulated response element (ISRE) and signal transducer and activator of transcription (STAT) were found in the 5' flanking region (2.0kb) of PoIRF10 gene. PoIRF10 mRNA was strongly expressed in gill, head kidney, heart, peripheral blood lymphocytes (PBLs), spleen and trunk kidney. PoIRF10 expression was up-regulated by Edwardsiella tarda, Streptococcus iniae and viral hemorrhagic septicemia virus (VHSV) infection in kidney. Lipopolysaccharide (LPS) and poly I:C also increased PoIRF10 expression level in PBLs. These results suggest that PoIRF10 is related to immune response in not only virus infection but also bacterial infection in teleost fish and should help to clarify the biological function of IRF10.

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.

Evolution of duplicated IgH loci in Atlantic salmon, Salmo salar.

BACKGROUND: The Atlantic salmon (Salmo salar) immunoglobulin heavy chain (IgH) locus possesses two parallel IgH isoloci (IGH-A and IGH-B), that are related to the genomic duplication event in the family Salmonidae. These duplicated IgH loci in Atlantic salmon provide a unique opportunity to examine the mechanisms of genome diversity and genome evolution of the IgH loci in vertebrates. In this study, we defined the structure of these loci in Atlantic salmon, and sequenced 24 bacterial artificial chromosome (BAC) clones that were assembled into the IGH-A (1.1 Mb) and IGH-B (0.9 Mb) loci. In addition, over 7,000 cDNA clones from the IgH variable (VH) region have been sequenced and analyzed. RESULTS: The present study shows that the genomic organization of the duplicated IgH loci in Atlantic salmon differs from that in other teleosts and other vertebrates. The loci possess multiple Cτ genes upstream of the Cµ region, with three of the Cτ genes being functional. Moreover, the duplicated loci possess over 300 VH segments which could be classified into 18 families. This is the largest number of VH families currently defined in any vertebrate. There were significant structural differences between the two loci, indicating that both IGH-A and -B loci have evolved independently in the short time after the recent genome duplication approximately 60 mya. CONCLUSIONS: Our results indicate that the duplication of the IgH loci in Atlantic salmon significantly contributes to the increased diversity of the antibody repertoire, as compared with the single IgH locus in other vertebrates.

Differential gene expression profiles in Japanese flounder (Paralichthys olivaceus) with different susceptibilities to edwardsiellosis.

Edwardsiellosis, which is caused by Edwardsiella tarda, a Gram-negative bacterium, is one of the most serious problems in Japanese flounder (Paralichthys olivaceus) culture. In this study, we examined the immune responses at the molecular level of genetic groups of Japanese flounder that are either susceptible or moderately susceptible to E. tarda infection using a cDNA microarray which was spotted with approximately 2000 different genes. Four different genetic groups of Japanese flounder (groups A, B, C and D) were infected with E. tarda by immersion. Mortality was 100% in groups A and C but only about 50% in groups B and D. Microarray analyses revealed 36 genes that were differentially expressed between the susceptible (A and C) and resistant (B and D) groups before E. tarda infection. Three days after the challenge, the resistant groups highly expressed MHC class I antigen processing and presentation-related genes, while the susceptible groups highly expressed genes involved in innate immune responses. The microarray results could be useful for selective breeding to enhance disease resistance of Japanese flounder against E. tarda and for studying strategies for controlling edwardsiellosis.

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.

Teleostean IL11b exhibits complementing function to IL11a and expansive involvement in antibacterial and antiviral responses.

Interleukin 11 is a class-1 helical cytokine, having the four-helix bundle structure, possessing pleiotropic characteristics involved in physiological processes including blood production, bone formation and placentation. The interleukin 11 paralogues (IL11a and IL11b) have been identified in fish with only IL11a from carp and trout have been characterized and analyzed for its expression thus far. Here, we cloned and studied the structure and expression of IL11b in Japanese flounder (Paralichthys olivaceus), and compared this with the existing information on fish IL11 paralogues. Japanese flounder IL11b (poIL11b) cDNA is composed of 1536 bp encoding for 201 aa residues with a 23 aa leader peptide, three cysteine residues (C12, C183 and C198) and four potential N-linked glycosylation sites. poIL11b does not show constitutive expression in tissues of adult fish except for the very slight expression in kidney and spleen, and the very high expression in peripheral blood leukocytes (PBLs). poIL11b is transiently up-regulated by bacterial lipopolysaccharide (LPS) and increasingly stimulated by the IFN inducer poly I:C in kidney, spleen and peripheral blood leukocytes of adult fish in vitro. It is likewise slightly stimulated by Edwardsiella tarda infection but is highly expressed after hirame rhabdovirus (HIRRV) infection in kidney of juvenile fish. The stimulation studies suggest that poIL11b, aside from its role in bacterial infection, is well involved in antiviral responses. Moreover, poIL11b structure and expression pattern appears to be slightly distinct and opposite to IL11a, respectively, suggesting a complementation of function of the duplicate fish IL11 genes.