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.

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.

Triple jeopardy: The combined effects of viral, chemical, and thermal stress on kuruma prawn (Penaeus japonicus) juveniles.

Growing concerns have emerged over the combined effects of multiple stressors on ecosystems. Empirical evidence shows that the sensitivity of aquatic invertebrates to insecticides varies under thermally fluctuating conditions. Additionally, field surveys in estuarine areas of western Japan confirmed the presence of juvenile kuruma prawns (Penaeus japonicus) carrying the white spot syndrome virus (WSSV). Given the potential of co-exposure to multiple stressors, we performed a combined exposure experiment using a full-factorial design with three stressors: WSSV infection (presence or absence: initial 2 h exposure), fipronil (insecticide) exposure (0 or 0.1 µg/L: 14 d exposure), and temperature (20, 25, or 30 °C). We observed the highest mortality (75 %) in the WSSV + Fipronil treatment at 30 °C, with the associated specimens showing significant changes in the internal load of WSSV and concentrations of fipronil and its metabolite, fipronil sulfone. Severe perturbations of metabolites associated with increased energy expenditure and fatty acid utilization have been identified as potential factors underlying lethality in juvenile kuruma prawns. The results demonstrate that WSSV infection increases the susceptibility of thermally stressed juvenile kuruma prawns to fipronil. Therefore, further studies are required to determine the combined effects of multiple stressors in environmentally relevant scenarios on juvenile kuruma prawns as well as in estuarine ecosystems.

Dual action of serotonin on local excitatory and inhibitory neural circuits regulating the corticotropin-releasing factor neurons in the paraventricular nucleus of the hypothalamus.

Serotonergic neurons originating from the raphe nuclei have been proposed to regulate corticotropin-releasing factor (CRF) neurons in the paraventricular nucleus of the hypothalamus (PVH). Since glutamate- and γ-aminobutyric acid (GABA)-containing neurons, constituting the hypothalamic local circuits, innervate PVH CRF neurons, we examined whether they mediate the actions of serotonin (5-hydroxytryptamine [5-HT]) on CRF neurons. Spontaneous excitatory postsynaptic currents (sEPSCs) or spontaneous inhibitory postsynaptic currents (sIPSCs) were recorded in PVH CRF neurons, under whole cell patch-clamp, using the CRF-modified yellow fluorescent protein (Venus) ΔNeo mouse. Serotonin elicited an increase in the frequency of sEPSCs in 77% of the cells and a decrease in the frequency of sIPSCs in 71% of the cells, tested in normal medium. Neither the amplitude nor decay time of sEPSC and sIPSC was affected, thus the site(s) of action of serotonin may be presynaptic. In the presence of tetrodotoxin (TTX), serotonin had no significant effects on either parameter of sEPSC or sIPSC, indicating that the effects of serotonin are action potential-dependent, and that the presynaptic interneurons are largely intact within the slice; distant neurons may exist, though, since some 20%-30% of neurons did not respond to serotonin without TTX. We next examined through what receptor subtype(s) serotonin exerts its effects on presynaptic interneurons. DOI (5-HT2A/2C agonist) mimicked the action of serotonin on the sIPSCs, and the serotonin-induced decrease in sIPSC frequency was inhibited by a selective 5-HT2C antagonist RS102221. 8-OH-DPAT (5-HT1A/7 agonist) mimicked the action of serotonin on the sEPSCs, and the serotonin-induced increase in sEPSC frequency was inhibited by a selective 5-HT7 antagonist SB269970. Thus, serotonin showed a dual action on PVH CRF neurons, by upregulating glutamatergic- and downregulating GABAergic interneurons; the former may partly be mediated by 5-HT7 receptors, whereas the latter by 5-HT2C receptors. The CRF-Venus ΔNeo mouse was useful for the electrophysiological examination.

Occurrence of neonicotinoids and fipronil in estuaries and their potential risks to aquatic invertebrates.

This study aimed to evaluate and qualify field-based potential risks of seven neonicotinoid and phenylpyrazole (fipronil) insecticides on aquatic invertebrates, including estuary-resident marine crustaceans. One hundred and ninety-three estuarine water samples, with salinity ranging from 0.5 to 32.7, were collected from four estuarine sites in the Seto Inland Sea of Japan, in 2015-2018 and the insecticide levels were measured. Five neonicotinoid and fipronil insecticides were successfully identified, and their occurrence varied temporally. Marine crustaceans were simultaneously harvested every month from one of the estuarine water sampling sites in 2015-2017. Three predominant crustacean species, kuruma prawn (Penaeus japonicus), sand shrimp (Crangon uritai), and mysid (Neomysis awatschensis), were captured and their seasonal presence was species independent. A 96-h laboratory toxicity study with the insecticides using kuruma prawn, sand shrimp, and a surrogate mysid species (Americamysis bahia) indicated that fipronil exerted the highest toxicity to the three crustaceans. Using both toxicity data and insecticide occurrence in estuarine water (salinity ≥10, n = 169), the potential risks on the three marine crustaceans were quantified by calculating the proportion of mixture toxicity effects (Pmix). The Pmix of seven neonicotinoids on the crustaceans was less than 0.8%, which is likely to be too low to indicate adverse effects caused by the insecticides. However, short temporal detection of fipronil (exclusively in June and July) significantly affected the Pmix, which presented the maximal Pmix values of 21%, 3.4%, and 72% for kuruma prawn, sand shrimp, and mysid, respectively, indicating a significant effect on the organisms. As for estuarine water (salinity <10), some water samples contained imidacloprid and fipronil exceeding the freshwater benchmarks for aquatic invertebrates. The present study provides novel insights into the seasonally varying risks of insecticides to estuarine crustaceans and highlights the importance of considering whether ecological risk periods coincide with crustacean presence.

Variable region of betanodavirus RNA2 is sufficient to determine host specificity.

Betanodaviruses, the causative agents of viral nervous necrosis in marine fish, have bipartite positive-sense RNA genomes. The viruses have been classified into 4 distinct types based on nucleotide sequence similarities in the variable region (the so-called T4 region) of the smaller genomic segment RNA2 (1.4 kb). Betanodaviruses have marked host specificity, although the primary structures of the viral RNAs and encoded proteins are similar among the viruses. We have previously demonstrated, using reassortants between striped jack nervous necrosis virus (SJNNV) and redspotted grouper nervous necrosis virus (RGNNV), that RNA2, which encodes the coat protein, strictly controls host specificity. However, because RNA2 is large, we were unable to propose a mechanism underlying this RNA2-based host specificity. To identify the RNA2 region that controls host specificity, we constructed RNA2 chimeric viruses from SJNNV and RGNNV and tested their infectivity in the original host fish, striped jack Pseudocaranx dentex and sevenband grouper Epinephelus septemfasciatus. Among these chimeric viruses, SJNNV mutants containing the variable region of RGNNV RNA2 infected sevenband grouper larvae in a manner similar to RGNNV, while RGNNV mutants containing the variable region of SJNNV RNA2 infected striped jack larvae in a manner similar to SJNNV. Immunofluorescence microscopic studies using anti-SJNNV polyclonal antibodies revealed that these chimeric viruses multiplied in the brains, spinal cords and retinas of the infected fish, as in infections by the parental viruses. These results indicate that the variable region of RNA2 is sufficient to control host specificity in SJNNV and RGNNV.

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.

A functional genomics tool for the Pacific bluefin tuna: Development of a 44K oligonucleotide microarray from whole-genome sequencing data for global transcriptome analysis.

Bluefin tunas are one of the most important fishery resources worldwide. Because of high market values, bluefin tuna farming has been rapidly growing during recent years. At present, the most common form of the tuna farming is based on the stocking of wild-caught fish. Therefore, concerns have been raised about the negative impact of the tuna farming on wild stocks. Recently, the Pacific bluefin tuna (PBT), Thunnus orientalis, has succeeded in completing the reproduction cycle under aquaculture conditions, but production bottlenecks remain to be solved because of very little biological information on bluefin tunas. Functional genomics approaches promise to rapidly increase our knowledge on biological processes in the bluefin tuna. Here, we describe the development of the first 44K PBT oligonucleotide microarray (oligo-array), based on whole-genome shotgun (WGS) sequencing and large-scale expressed sequence tags (ESTs) data. In addition, we also introduce an initial 44K PBT oligo-array experiment using in vitro grown peripheral blood leukocytes (PBLs) stimulated with immunostimulants such as lipopolysaccharide (LPS: a cell wall component of Gram-negative bacteria) or polyinosinic:polycytidylic acid (poly I:C: a synthetic mimic of viral infection). This pilot 44K PBT oligo-array analysis successfully addressed distinct immune processes between LPS- and poly I:C- stimulated PBLs. Thus, we expect that this oligo-array will provide an excellent opportunity to analyze global gene expression profiles for a better understanding of diseases and stress, as well as for reproduction, development and influence of nutrition on tuna aquaculture production.

Visualization of corticotropin-releasing factor neurons by fluorescent proteins in the mouse brain and characterization of labeled neurons in the paraventricular nucleus of the hypothalamus.

Corticotropin-releasing factor (CRF) is the key regulator of the hypothalamic-pituitary-adrenal axis. CRF neurons cannot be distinguished morphologically from other neuroendocrine neurons in the paraventricular nucleus of the hypothalamus (PVH) without immunostaining. Thus, we generated a knock-in mouse that expresses modified yellow fluorescent protein (Venus) in CRF neurons (CRF-Venus), and yet its expression is driven by the CRF promoter and responds to changes in the interior milieu. In CRF-Venus, Venus-expressing neurons were distributed in brain regions harboring CRF neurons, including the PVH. The majority of Venus-expressing neurons overlapped with CRF-expressing neurons in the PVH, but many neurons expressed only Venus or CRF in a physiological glucocorticoid condition. After glucocorticoid deprivation, however, Venus expression intensified, and most Venus neurons coexpressed CRF. Conversely, Venus expression was suppressed by excess glucocorticoids. Expression of copeptin, a peptide encoded within the vasopressin gene, was induced in PVH-Venus neurons by glucocorticoid deprivation and suppressed by glucocorticoid administration. Thus, Venus neurons recapitulated glucocorticoid-dependent vasopressin expression in PVH-CRF neurons. Noradrenaline increased the frequency of glutamate-dependent excitatory postsynaptic currents recorded from Venus-expressing neurons in the voltage clamp mode. In addition, the CRF-iCre knock-in mouse was crossed with a CAG-CAT-EGFP reporter mouse to yield the Tg(CAG-CAT-EGFP/wt);CRF(iCre/wt) (EGFP/CRF-iCre) mouse, in which enhanced green fluorescent protein (EGFP) is driven by the CAG promoter. EGFP was expressed more constitutively in the PVH of EGFP/CRF-iCre mice. Thus, CRF-Venus may have an advantage for monitoring dynamic changes in CRF neurons and CRF networks in different glucocorticoid states.