High-Resolution Melting Analysis of COI Sequences Distinguishes Pufferfish Species (Takifugu spp.) in China

Pufferfish is a traditional, delicious dish in Asia. However, eating wild or improperly processed pufferfish causes serious poisoning. This study aimed to exploit the high-resolution melting (HRM) method for authenticating four species of Takifugu pufferfish (Takifugu xanthopterus, T. fasciatus, T. flavidus, and T. rubripes). Candidate DNA barcodes, including the cytochrome c oxidase subunit I (COI), cytochrome oxidase b (Cytb), and the control region (D-loop), were analyzed, with COI selected as the optimal DNA barcode. An HRM method was developed to identify 57 commercial fish samples in China, including 33 commercial pufferfish products and 24 unlabeled fish products. The findings revealed that the pufferfish products were T. rubripes or T. fasciatus, and four T. xanthopterus samples were detected in unlabeled fish products. These results showed that DNA barcode coupled with HRM analysis was a rapid and efficient tool to identify pufferfish, which might aid in the prevention of consumer fraud or mislabeling of fish products.

Expression dynamics of the genes for the hypothalamo-pituitary-gonadal axis in tiger puffer (Takifugu rubripes) at different reproductive stages

Tiger puffer, Takifugu rubripes, a commercially important long-distance migratory fish, return to specific spawning grounds for reproduction. To clarify reproductive neuroendocrine system of the tiger puffer, the changes in the expression levels of the genes encoding three gonadotropin-releasing hormones (GnRHs), gonadotropin-inhibitory hormone (GnIH), GnIH receptor (GnIH-R), kisspeptin and kisspeptin receptor in the brain and gonadotropin (GTH) subunits, growth hormone (GH) and prolactin (PRL) in the pituitary were examined in the tiger puffer captured in the wild at different reproductive stages, namely immature and mature fish of both sexes, and post-ovulatory females that were obtained by hormonal treatment. The amounts of three gnrh mRNAs, gnih, gnih-r, fshb and lhb were substantially increased in the mature fish compared to the immature fish, especially in the females, and these augmented expressions were drastically decreased in the post-ovulatory females. gh expression showed a slight increase in the mature males. In contrast, kiss2, kiss2r and prl did not show significant changes in the males but significantly decreased in the post-ovulatory females. The present results demonstrate the expression dynamics of the hypothalamo-pituitary-gonadal axis genes associated with the reproductive conditions and the possible involvement of the GnRH/GnIH/GTH system in the regulation of the sexual maturation and spawning in the wild tiger puffer.

Effects of the Thyroid Endocrine System on Gonadal Sex Ratios and Sex-Related Gene Expression in the Pufferfish Takifugu rubripes

To examine the effect and mechanism of thyroid hormone on gonadal sex differentiation, Takifugu rubripes larvae were treated with goitrogen (methimazole, MET, 1000 g/g), and thyroxine (T4, 2nM) from 25 to 80 days after hatching (dah). Gonadal histology and sex ratios of fish were then determined at 80 dah. MET treatment induced masculinization, but T4 treatment did not induce feminization in T. rubripes larvae. Transcriptomic analysis of gonads at 80 dah was then conducted. Among the large number of differentially expressed genes between the groups, the expression of foxl2, cyp19a1a, and dmrt1 was altered. The expression of foxl2, cyp19a1a, dmrt1 and gsdf at 25, 40, 55 days after treatment (dat) was further analyzed by qPCR. MET treatment suppressed the expression of foxl2 and cyp19a1a, and induced the expression of dmrt1 in genetic females (p < 0.05). Additionally, T4 treatment induced an increase in the expression of cyp19a1a in genetic XY gonads only at 25 dat. However, the increase in cyp19a1a expression did not continue to 40 and 55 dat. This may explain why feminization of larvae was not found in the T4-treated group. Thus, the present study provides the first evidence that MET treatment causes masculinization in teleost fish. The effects of MET-induced masculinization in T. rubripes may act primarily via suppression of the expression of foxl2 and cyp19a1a, and stimulation of the expression of dmrt1. Moreover, the effects of higher concentrations of T4 or different concentrations of T3, on sex differentiation require further testing.

Genomic selection for heterobothriosis resistance concurrent with body size in the tiger pufferfish, Takifugu rubripes

Parasite resistance traits in aquaculture species often have moderate heritability, indicating the potential for genetic improvements by selective breeding. However, parasite resistance is often synonymous with an undesirable negative correlation with body size. In this study, we first tested the feasibility of genomic selection (GS) on resistance to heterobothriosis, caused by the monogenean parasite Heterobothrium okamotoi, which leads to huge economic losses in aquaculture of the tiger pufferfish Takifugu rubripes. Then, using a simulation study, we tested the possibility of simultaneous improvement of parasite resistance, assessed by parasite counts on host fish (HC), and standard length (SL). Each trait showed moderate heritability (square-root transformed HC: h2 = 0.308 ± 0.123, S.E.; SL: h2 = 0.405 ± 0.131). The predictive abilities of genomic prediction among 12 models, including genomic Best Linear Unbiased Predictor (GBLUP), Bayesian regressions, and machine learning procedures, were also moderate for both transformed HC (0.248‒0.344) and SL (0.340‒0.481). These results confirmed the feasibility of GS for this trait. Although an undesirable genetic correlation was suggested between transformed HC and SL (rg = 0.228), the simulation study suggested the desired gains index can help achieve simultaneous genetic improvements in both traits.

Diel Variation in CC Chemokine Gene Expression in the Japanese Pufferfish Takifugu rubripes

CC chemokines are key molecules in the regulation of leukocyte trafficking to the site of injury, infection, or inflammation. In recent years, some mammalian chemokines have been shown to exhibit rhythmic expression, regulated by clock genes. However, the rhythmic expression of chemokines in teleost fish remains unknown. In the present study, the diel variation of teleost CC chemokine genes was investigated using the model fish, Fugu (Takifugu rubripes). Diel variation analysis revealed that clock (bmal1, clock1, per2, rorα, and rev-erbβ) and CC chemokine (ccl18l, ccl19, and ccl25l) genes show diel expression under 12:12 light-dark cycle (LD12:12) conditions. CC chemokine genes, which exhibit diel expression, contain RORE (ccl18l, ccl19, ccl25l) and/or E-box (ccl25l) motifs in their transcription regulatory region. Moreover, in vitro head kidney stimulation with lipopolysaccharide (LPS) at different zeitgeber times (ZT) under LD12:12 conditions affected the degree of ccl18l, ccl19, and ccl25l expression; high and low responsiveness to LPS stimulation at ZT12 and ZT0 (ccl25l), and ZT16 and ZT4 (ccl18l and ccl19), respectively, were observed. These results suggest that the expression of some fish CC chemokines is affected by the diel variation regulated by clock proteins, and that responsiveness against bacterial infection depends on the time zone.

Classification of Takifugu rubripes, T. chinensis and T. pseudommus by genotyping-by-sequencing

Takifugu rubripes is more expensive than other species of the genus because of its high protein content and special flavor. However, it is easily confused with imported T. chinensis and T. pseudommus because they have similar morphological characteristics. We identified single nucleotide polymorphism (SNP) markers of T. rubripes by genotyping-by-sequencing (GBS) and evaluated their ability to distinguish among T. rubripes, T. chinensis, and T. pseudommus. In all, 18 polymorphic SNPs were subjected to phylogenetic analyses of the three Takifugu species. Additionally, we subjected a second set of samples to Sanger sequencing to verify that the polymorphic SNPs could be used to evaluate the genetic variation among the three Takifugu species. A phylogenetic tree that included the analyzed sequence of set A, which is referred to as the reference sequence, and a validation sequence of set B with 18 SNPs were produced. Based on this phylogenetic tree and STRUCTURE analyses, T. rubripes, T. chinensis and T. pseudommus have low genetic variation and should be considered the same gene pool. Our findings suggest that further studies are needed to estimate the genetic association of the three Takifugu species.

Interleukin-17A/F1 from Japanese pufferfish (Takifugu rubripes) stimulates the immune response in head kidney and intestinal cells

In mammals, interleukin (IL)-17A and IL-17F, mainly produced by Th17 cells, are hallmark inflammatory cytokines that play important roles in the intestinal mucosal immune response. In contrast, three mammalian IL-17A and IL-17F counterparts (IL-17A/F1-3) have been identified in teleosts, and most of their functions have been described in the lymphoid organs. However, their function in the intestinal mucosal immune response is poorly understood. In this study, a recombinant (r) tiger puffer fish fugu (Takifugu rubripes) IL-17A/F1 was produced and purified using a mammalian expression system, and was used to stimulate cells isolated from fugu head kidney and intestines. The gene expression levels of TNF-α, IL-1β, IL-6, and β-defensin-like protein-1 (BD-1) genes were evaluated at 0, 3, 6 and 12 h post-stimulation (hps). Phagocytic activity and superoxide anion production were evaluated at the same time points using an NBT assay. The rIL-17A/F1 protein was shown to induce the expression of pro-inflammatory cytokines and antimicrobial peptides in both head kidney and intestinal cells. Expression levels for IL-1β, TNF-α, and IL-6 were all up-regulated between 3 and 12 hps. In addition, stimulation with rIL-17A/F1 enhanced phagocytic activity at 24 hps. Superoxide anion production was increased at 48 hps in the head kidney cells and moderately increased at 48 hps in intestinal cells. This study suggests that fugu IL-17A/F1 plays an important role in promoting the innate immune response and may act as a bridge between innate and adaptive immunity in the head kidney and intestine.

Functional characterization of an immunity-related GTPase gene in immune defense from obscure puffer (Takifugu obscurus)

Immunity-related GTPases (IRGs) are a family of large interferon-inducible GTPases that function in effective host defense against invading pathogens. IRGs have been extensively studied in mammals for their roles in the elimination of intracellular pathogens; however, their homologs in lower vertebrates are not well known. In this study, an IRG from obscure puffer (Takifugu obscurus), ToIRG, was identified and further characterized for its functional activity. The ToIRG gene encodes a protein of 396 amino acids containing a typical N-terminal GTPase domain with three conserved motifs. Phylogenetic analysis revealed that it has a closer evolutionary relationship with mammalian GKS IRGs. Gene expression profile analysis revealed that ToIRG was ubiquitously expressed in all tested healthy tissues of obscure puffer and upregulated in response to Aeromonas hydrophila or Edwardsiella tarda challenge. The subcellular localization of ToIRG is characterized as condensed forms around the nucleus. Importantly, an antimicrobial assay in vitro suggested that ToIRG enhanced the ability of host cells to resist both intracellular (E. tarda) and extracellular pathogens (A. hydrophila). Taken together, these results provide the functional characterization of obscure puffer IRGs in immune defense, which is the first study to reveal the function of IRGs in bony fish and will provide important insights into the evolutionary divergence of IRGs.

Establishment of a Japanese medaka (Oryzias latipes) transgenic line expressing Takifugu rubripes pufferfish saxitoxin and tetrodotoxin binding protein 1, and evaluation of tributyltin toxicity via in ovo nanoinjection

Pufferfish saxitoxin and tetrodotoxin binding proteins (PSTBPs) play an important role in the toxification of certain species of pufferfish. Recombinant Takifugu rubripes PSTBP1 (rTrub.PSTBP1) is reported to bind to tributyltin, and so it has been suggested that rTrub.PSTBP1 may reduce the toxicity of tributyltin. However, the role of PSTBP1 in vivo remains to be elucidated. Here, we established a transgenic medaka line showing whole-body Renilla reniformis green fluorescent protein and Trub.PSTBP1 expression, as confirmed by real-time polymerase chain reaction and mRNA-Seq analysis. mRNA-Seq analysis also showed that cytochrome P450 superfamily genes and the gene encoding ATP-binding cassette sub-family G member 2 were highly expressed in the transgenic medaka. Using embryos of the transgenic medaka line, we conducted an in ovo nanoinjection test to examine the effect of Trub.PSTBP1 in vivo, and obtained data suggesting that Trub.PSTBP1 expression may have reduced the toxicity of tributyltin in our transgenic medaka line. Our findings will be useful for future functional analyses of Trub.PSTBP1.

Tissue distribution of transcription for 29 lipid metabolism-related genes in Takifugu rubripes, a marine teleost storing lipid predominantly in liver

The tissue distribution pattern of lipid is highly diverse among different fish species. Tiger puffer has a special lipid storage pattern, storing lipid predominantly in liver. In order to better understand the lipid physiology in fish storing lipid in liver, the present study preliminarily investigated the tissue distribution of transcription for 29 lipid metabolism-related genes in tiger puffer, which are involved in lipogenesis, fatty acid oxidation, biosynthesis and hydrolysis of glycerides, lipid transport, and relevant transcription regulation. Samples of eight tissues, brain, eye, heart, spleen, liver, intestine, skin, and muscle, from fifteen juvenile tiger puffer were used in the qRT-PCR analysis. The intestine and brain had high transcription of lipogenic genes, whereas the liver and muscle had low expression levels. The intestine also had the highest transcription level of most apolipoproteins and lipid metabolism-related transcription factors. The transcription of fatty acid β-oxidation-related genes was low in the muscle. The peroxisomal fatty acid oxidation may dominate over mitochondrial β-oxidation in the liver and intestine of tiger puffer, and the MAG pathway probably predominates over the G3P pathway in re-acylation of absorbed lipids in the intestine. The intracellular glyceridases were highly transcribed in the brain, eye, and heart. In conclusion, in tiger puffer, the intestine could be a center of lipid metabolism whereas the liver is more likely a pure storage organ for lipid. The lipid metabolism in the muscle could also be inactive, possibly due to the very low level of intramuscular lipid.

Transcriptome analysis of Takifugu obscurus liver in response to acute retene exposure

Retene (1-methyl-7-isopropyl-phenanthrene, RET) is an alkyl polycyclic aromatic hydrocarbon (PAH) with environmental risk to aquatic animals. Takifugu obscurus is a migratory fish species with high economic and ecological value. To assess the toxic effects of RET on molecular metabolism, juvenile T. obscurus in this study were acutely exposed to 44.30 µg/L of RET for four days. The transcriptome profiles of livers were compared between RET treatment group and the control, and the results revealed that 1,897 genes were significantly differentially expressed (DEGs) after exposure to RET, which enriched 17 KEGG pathways. Among these, glycerolipid metabolism, glycerophospholipid metabolism, insulin signaling pathway, and FOXO signaling pathways were significantly activated. Further exploration indicated that RET exposure disrupted glucose metabolism, stimulated insulin metabolism, and activated cell proliferation genes. Overall, these findings help explain the molecular mechanisms underlying RET toxicity, and may offer evidence to support T. obscurus protection.

Annual patterns of ocular melatonin level in the female grass puffer, Takifugu alboplumbeus: possible involvement in seasonal reproductive response

The aim of this study was to investigate the expression patterns of ocular melatonin in the annual reproductive cycle of the female grass puffer. Spawning season of the female grass puffer is from June to July in Jeju, South Korea. Time-resolved fluoroimmunoassay revealed that levels of ocular melatonin, which show an annual change, peaked in May (spawning season). Additionally, expression of reproductive-related genes also showed annual patterns: GnRH1 peaked in August, GnRH2 peaked in February, GnRH3, Kiss2, and LPXRFa peaked in November. These results suggest that ocular melatonin may be related to the annual reproductive cycle in the grass puffer. To better understand the photic regulation of AANAT1a mRNA in the retina, we observed the nocturnal pattern of ocular melatonin levels daily, which shows a nocturnal pattern in both short photoperiod (SD) and long photoperiod (LD) conditions. In the brain, AANAT2 mRNA also shows a nocturnal pattern in both SD and LD; however, the time of peak expression of AANAT2 mRNA was unchanged in both conditions. Following intraperitoneal injection of melatonin for 2 weeks, expression of GnRH2 and LPXRFa mRNA in the brain significantly increased, while that of Kiss2 mRNA was decreased, suggesting that melatonin has a reproduction-related effect. Furthermore, under SD and LD conditions for 14 weeks, the gonadosomatic index more increased and the maturity of the ovary progressed under LD compared with those under SD, suggesting that the SD photoperiodic signal inactivated ovarian development. These results indicate that the ocular melatonin may have a possible role in the reproductive endocrinology of the grass puffer.

Toll-like receptors (TLRs) respond to tributyltin chloride (TBT-Cl) exposure in the river pufferfish (Takifugu obscurus): Evidences for its toxic injury function

Tributyltin chloride (TBT-Cl) residual in water body had become a noticeable ecological problem for aquatic ecosystems. Toll-like receptors (TLRs) are an ancient family of pattern recognition receptors that play key roles in detecting nonself antigens and immune system activation. In this study, we explored the effect of TBT-Cl exposure on four TLRs expression in river pufferfish, Takifugu obscurus. The four T. obscurus Toll-like receptors (To-TLRs) contained different types of domains such as leucine-rich repeats (LRRs), leucine-rich repeats, typical subfamily (LRR_TYP) and other special domains. The To-TLRs mRNA transcripts expressed in all tissues, also To-TLR2 was investigated with higher level in kidney, as well as To-TLR3 in kidney, while To-TLR18 in liver and To-TLR22 in intestine. After the acute and chronic exposure of TBT-Cl, To-TLR2 and To-TLR3 mRNA transcripts were significantly down-regulated in gill. However, To-TLR18 and To-TLR22 were significantly up-regulated in gill and liver. Moreover, the histology and immunohistochemistry (IHC) results showed the different injury degrees of TBT-Cl in liver and gill and implied the cytoplasm reorganization after TBT-Cl stress and the function of immunoregulation for To-TLRs to TBT-Cl exposure. All the results indicated that To-TLRs might involve in sensing and mediating innate immune responses caused by TBT-Cl for keeping detoxification homeostasis.

Construction of a High-Density Genetic Linkage Map and QTL Mapping for Growth-Related Traits in Takifugu bimaculatus

Takifugu bimaculatus is a euryhaline species, distributed ranging from the southern Yellow Sea to the South China Sea. Their tolerance to a wide range of salinity and temperature, coupled with a desirable firm texture, makes T. bimaculatus a strong candidate for Takifugu aquaculture in subtropics areas. Due to the increasing demand in markets and emerging of the Takifugu aquaculture industry, close attention has been paid to improvement on the T. bimaculatus production. In aquaculture, the great effort has been put into marker-assisted selective breeding, and efficient improvement was realized. However, few genetic resources on T. bimaculatus are provided so far. Aiming at understanding the genetic basis underlying important economic growth traits, facilitating genetic improvement and enriching the genetic resource in T. bimaculatus, we constructed the first genetic linkage map for T. bimaculatus via double digestion restriction-site association DNA sequencing and conducted quantitative traits locus (QTL) mapping for growth-related traits. The map comprised 1976 single nucleotide polymorphism markers distributed on 22 linkage groups (LG), with a total genetic distance of 2039.74 cM. Based on the linkage map, a chromosome-level assembly was constructed whereby we carried out comparative genomics analysis, verifying the high accuracy on contigs ordering of the linkage map. On the other hand, 18 QTLs associated with growth traits were detected on LG6, LG7, LG8, LG10, LG20, and LG21 with phenotypical variance ranging from 15.1 to 56.4%. Candidate genes participating in cartilage development, fat accumulation, and other growth-related regulation activities were identified from these QTLs, including col11a1, foxa2, and thrap3. The linkage map provided a solid foundation for chromosomes assembly and refinement. QTLs reported here unraveled the genomic architecture of some growth traits, which will advance the investigation of aquaculture breeding efforts in T. bimaculatus.

Biochemical and molecular identification of a novel hepcidin type 2-like antimicrobial peptide in the skin mucus of the pufferfish Takifugu pardalis

Fish skin mucus is considered to act as the first line of defense against waterborne pathogens and to be potential source of novel antimicrobial components. Here we report the purification and characterization of a novel hepcidin type 2-like antimicrobial peptide (TpHAMP2) from the skin mucus of the pufferfish Takifugu pardalis. The purified TpHAMP2 comprised of 23 amino acids (AAs) with eight Cys residues that form four intramolecular disulfide bonds. The TpHAMP2 gene shared overall structural characteristics with all known hepcidins, which have a tripartite exon-intron gene organization and three structural signatures in the precursor protein. Phylogenetically, TpHAMP2 was classified as HAMP2 class in acanthopterygian fish. Interestingly, the AA sequence of TpHAMP2 did not contain a proprotein cleavage site (RXXR motif) that conserved in most hepcidins and showed a highly positive charged (RKR-) short N-terminus and Val¹⁸ and Gly²² residues, which are distinctive structures compared to other known active hepcidins. Recombinant TpHAMP2 identical to the native form exhibited a broad spectrum and potent antimicrobial activity against tested gram-positive and -negative bacteria. Expression of TpHAMP2 mRNA was predominant in the liver and was upregulated in the liver, the spleen, the intestine, and the skin of T. pardalis post immune challenge. Thus, our findings suggests that TpHAMP2 might be of importance in the framework of discovering the fish hepcidins, especially type 2s, and provide noteworthy insight into its gene structure and expression and in the innate immunity as well as the mucosal immunity in regard to hepcidins' evolutionary history in fish species.

Effects of high and low temperature on expression of GnIH, GnIH receptor, GH and PRL genes in the male grass puffer during breeding season

Gonadotropin-inhibitory hormone (GnIH) is a multifunctional hypophysiotropic neurohormone and has a stimulatory role in the control of reproduction in the grass puffer. To clarify the neuroendocrine mechanisms underlying the effect of changes in water temperature on reproduction in fish, we previously revealed that, in parallel to gonadal regression, both low and high temperature significantly decreased the expressions of the genes encoding kisspeptin (kiss2), kisspeptin receptor (kiss2r), gonadotropin-releasing hormone 1 (gnrh1) in the brain and gonadotropin (GTH) subunits (fshb and lhb) in the pituitary of sexually mature male grass puffer. In this study, we examined the changes in expression of gnih and GnIH receptor gene (gnihr) in the brain and pituitary along with the genes for growth hormone (gh) and prolactin (prl) in the pituitary of male grass puffer exposed to low temperature (14 °C), normal temperature (21 °C, as initial control) and high temperature (28 °C) conditions for 7 days. The levels of gnih and gnihr mRNAs were significantly decreased in both low and high temperature conditions compared to normal temperature in the brain and pituitary. Similarly, the gh mRNA levels were significantly decreased in both low and high temperature conditions. The prl mRNAs showed no significant changes at high temperature, whereas drastically decreased at low temperature possibly by dysfunctional cold stress. Taken together, the present results suggest that, in addition to the inhibitory effect of temperature changes on the Kiss2/GnRH1/GTH system, the suppression of GnIH/GH system may also be involved in the termination of reproduction by high temperature at the end of breeding season.

The sequence and de novo assembly of Takifugu bimaculatus genome using PacBio and Hi-C technologies

Takifugu bimaculatus is a native teleost species of the southeast coast of China where it has been cultivated as an important edible fish in the last decade. Genetic breeding programs, which have been recently initiated for improving the aquaculture performance of T. bimaculatus, urgently require a high-quality reference genome to facilitate genome selection and related genetic studies. To address this need, we produced a chromosome-level reference genome of T. bimaculatus using the PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The genome was assembled into 2,193 contigs with a total length of 404.21 Mb and a contig N50 length of 1.31 Mb. After chromosome-level scaffolding, 22 chromosomes with a total length of 371.68 Mb were constructed. Moreover, a total of 21,117 protein-coding genes and 3,471 ncRNAs were annotated in the reference genome. The highly accurate, chromosome-level reference genome of T. bimaculatus provides an essential genome resource for not only the genome-scale selective breeding of T. bimaculatus but also the exploration of the evolutionary basis of the speciation and local adaptation of the Takifugu genus.

Alternations in oxidative stress, apoptosis, and innate-immune gene expression at mRNA levels in subadult tiger puffer (Takifugu rubripes) under two different rearing systems

Tiger puffer (Takifugu rubripes) is one of the major aquaculture fish species in China due to its high economic value. In this study, the transcriptions of hepatic antioxidant enzyme, stress, apoptosis, and immune-related genes of sub-adult tiger puffers (Takifugu rubripes) were evaluated under two different rearing systems [offshore sea cage aquaculture system (OSCS) and recirculating aquaculture system (RAS)]. Results showed that the mRNA expression levels of the antioxidant enzyme (mn-sod, cu/zn-sod, gpx, and gr) and stress-related (hsp70 and hsp90) genes of male tiger puffers reared in the OSCS were significantly higher than female fish reared in the OSCS and fish reared in the RAS. The anti-apoptotic gene bcl2 exhibited the similar results. By contrast, the mRNAs of the pro-apoptotic genes (p53, caspase8, caspase9, and caspase3) of male tiger puffers reared in the OSCS were significantly lower than female fish reared in the OSCS and fish reared in the RAS. Male tiger puffers reared in the OSCS displayed significantly higher complement components (c3) and inflammatory cytokine (il-6) mRNAs, whereas B-cell activating factor (baf) and tumor necrosis factor α (tnf-α) mRNAs remained unchanged. Meanwhile, the mRNA levels of pro-apoptotic (bax, caspase8) and immunity-related (c3, il-6 and il-7) genes of female tiger puffers reared in the OSCS were significantly lower and higher than female fish reared in the RAS, respectively. In conclusion, the hepatic antioxidant, anti-apoptosis, and innate immunity of tiger puffers reared in the OSCS were better than fish in the RAS, male tiger puffer obtained the best values. These results expand the knowledge on the combined RAS and OSCS alternative aquaculture model for tiger puffers and aid in their management in captive.

Changes in DNA methylation during epigenetic-associated sex reversal under low temperature in Takifugu rubripes

DNA methylation has frequently been implicated in sex determination and differentiation in teleost species. In order to detect the DNA methylation patterns established during sexual differentiation in tiger pufferfish T. rubripes, we performed comprehensive whole genome methylation sequencing and analyses of the gonads of male, female, and pseudo male. We obtained a total of 33.12, 32.44, and 31.60 Gb clean data for male, female, and pseudo male, with a sequencing depth of 66.44×, 60.47× and 54.86×, respectively. The methylation level of cytosine (C) residues in the genomic DNA from gonads was 11.016%, 10.428%, and 11.083% in male, female, and pseudo male, respectively. More than 65% of C methylation was at CpG sites, and less than 1% was at CHG and CHH sites. In each regulatory element, there were low methylation levels on both sides of the transcription start site, and higher methylation levels in exons, introns, and downstream of genes. The highest mCpG was on chromosome 8 and the lowest mCpG was on chromosome 5. Comparisons of whole-genome DNA methylation between pairs of samples revealed that there were 3,173 differentially methylated regions (DMRs) between female and male, and 3,037 DMRs between male and pseudo male, corresponding to 0.232% and 0.223% of the length of the genome, respectively. There were only 1,635 DMRs between female and pseudo male, representing 0.127% of the length of the genome. A number of differentially methylated genes (DMGs) related to sex determination and differentiation were selected, such as amhr2 and pfcyp19a. After Bisulfite Sequencing PCR (BSP) verification, amhr2 was exhibited low methylation level in normal males and pseudo male, and high methylation level in normal females but pfcyp19a showed low methylation level in normal females and high methylation level in normal males and pseudo males. These results provide information about the molecular epigenetic mechanisms of DNA methylation during low-temperature induced masculinization of tiger pufferfish, and increase our understanding of the mechanisms of sex determination and differentiation in this important aquaculture fish species.

The key differentially expressed genes and proteins related to immune response in the spleen of pufferfish (Takifugu obscurus) infected by Aeromonas hydrophila

The immune mechanism elicited in pufferfish (Takifugu obscurus) against the invasion of Aeromonas hydrophila is still poorly understood. We examined the spleen of pufferfish at the transcriptome and proteome levels by using Illumina-seq and TMT coupled mass spectrometry after 12 h infection by A. hydrophila, respectively. A total of 2,339 genes (1,512 up-regulated and 827 down-regulated) and 537 (237 up-regulated and 300 down-regulated) proteins were identified. GO and KEGG analyses revealed that the responses to stimulus were the main biological processes, intestinal immune network for IgT production and calcium signaling pathway. Fourteen genes (8 up-regulated and 6 down-regulated) and proteins (5 up-regulated and 9 down-regulated) involved immune responses or signal transduction were validated by qRT-PCR and parallel reaction monitoring to confirm the reliability of the transcriptomic and proteomic analyses, respectively. Moreover, qRT-PCR and flow cytometry were used to detect dynamics of the genes in calcium signaling pathway and changes of concentration of cytoplasm Ca²⁺ in spleen cells within a 72 h challenge. This study provides the findings regarding immune response, especially intestinal immune network for IgT production pathway and calcium signaling pathway at the molecular, protein and cellular in pufferfish after infection by A. hydrophila. These results would provide a new insight and molecular targets into the response to pathogenic infection in pufferfish.

Integrated application of multi-omics provides insights into cold stress responses in pufferfish Takifugu fasciatus

BACKGROUND

T. fasciatus (Takifugu fasciatus) faces the same problem as most warm water fish: the water temperature falls far below the optimal growth temperature in winter, causing a massive death of T. fasciatus and large economic losses. Understanding of the cold-tolerance mechanisms of this species is still limited. Integrated application of multi-omics research can provide a wealth of information to help us improve our understanding of low-temperature tolerance in fish.

RESULTS

To gain a comprehensive and unbiased molecular understanding of cold-tolerance in T. fasciatus, we characterized mRNA-seq and metabolomics of T. fasciatus livers using Illumina HiSeq 2500 and UHPLC-Q-TOF MS. We identified 2544 up-regulated and 2622 down-regulated genes in the liver of T. fasciatus. A total of 40 differential metabolites were identified, including 9 down-regulated and 31 up-regulated metabolites. In combination with previous studies on proteomics, we have established an mRNA-protein-metabolite interaction network. There are 17 DEMs (differentially-expressed metabolites) and 14 DEGs-DEPs (differentially co-expressed genes and proteins) in the interaction network that are mainly involved in fatty acids metabolism, membrane transport, signal transduction, and DNA damage and defense. We then validated a number of genes in the interaction network by qRT-PCR. Additionally, a number of SNPs (single nucleotide polymorphisms) were revealed through the transcriptome data. These results provide key information for further understanding of the molecular mechanisms of T. fasciatus under cold stress.

CONCLUSION

The data generated by integrated application of multi-omics can facilitate our understanding of the molecular mechanisms of fish response to low temperature stress. We have not only identified potential genes and SNPs involved in cold tolerance, but also show that some nutrient metabolites may be added to the diet to help the overwintering of T. fasciatus.

Comparative transcriptomic analysis reveals the gene expression profiles in the liver and spleen of Japanese pufferfish (Takifugu rubripes) in response to Vibrio harveyi infection

Japanese pufferfish (Takifugu rubripes) is one of the main marine aquatic fish species cultured in Asia due to its high nutritional value. In recent years, disease caused by Vibrio harveyi infections have led to serious mortality in Japanese pufferfish industry. To understand the complex molecular mechanisms between V. harveyi and Japanese pufferfish, we performed a transcriptome analysis of liver and spleen samples from Japanese pufferfish at 1 and 2 day post-infection. Between-group comparisons revealed 922 genes that were significantly differentially expressed. The altered genes emphasized the function in several immune related pathways including MAPK signaling pathway, JAK-STAT signaling pathway, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction and lysosomal pathway. The data generated in this study provided insight into the responses of Japanese pufferfish against V. harveyi at the transcriptome level, promoting our comprehensive understanding of immune responses for aquatic animal against V. harveyi.

Molecular characterization and expression analysis of galectins in Japanese pufferfish (Takifugu rubripes) in response to Vibrio harveyi infection

Galectins are a family of proteins with conserved carbohydrate recognition domains (CRDs) that bind to specific glycans, including the glycans on the surface of pathogens, and therefore play a role in cytokine secretion, cell activation, migration, adhesion and apoptosis. Currently, galectins have been extensively studied in mammalian species but rarely studied in teleost fish species. In this study, a total of 12 galectin genes were characterized to understand the molecular mechanisms of galectin function in Japanese pufferfish (Takifugu rubripes). Phylogenetic analyses and syntenic analyses confirmed their correct annotation and suggested the strongest relationships to tetraodon. Furthermore, expression analyses were conducted in healthy tissues of Japanese pufferfish and after infection with Vibrio harveyi in the intestine, liver and spleen. The results showed that galectin genes were widely expressed in all examined tissues; however, most of the galectin genes were highly expressed in mucosal tissues (skin, gill and intestine). Moreover, majority of the galectin genes were significantly regulated after V. harveyi infection in the intestine, liver and spleen, suggesting that galectins were involved in the immune response to V. harveyi infection in Japanese pufferfish. This study established the foundation for future studies of galectin gene functions.

Comparative transcriptomics analysis of the river pufferfish (Takifugu obscurus) by tributyltin exposure: Clues for revealing its toxic injury mechanism

TBT residual in water had become a noticeable ecological problem for aquatic ecosystems. The river pufferfish (Takifugu obscurus) is a kind of an anadromous fish species and widely distributed in the East China Sea and the Yellow Sea. Because of the water contamination, the pufferfish wild resource had a sudden decline in recent years. Therefore, the study on the response of pufferfish to the TBT exposure may contribute to reveal toxic injury mechanism of T. obscurus under TBT exposure. In this study, the transcriptional library of T. obscurus liver and gill was constructed and sequenced by an improved Illumina HiseqX10 high-throughput sequencing platform under different concentrations of TBT acute stress. The blood cell numbers distinctly decreased after TBT exposure, showing the adverse effects of TBT invasion and self-adjusting ability of the pufferfish. The production of reactive oxygen species increased, demonstrating the oxidation resistance of T. obscurus when exposed to TBT. The obtained data were compared with the genome data of Takifugu rubripes and transcriptional resource database. On this basis, gene function annotation, analysis and classification were carried out by bioinformatics method, and differential genes related to toxic injury function were screened out. Meanwhile, new toxic related genes and related signal pathways were sought to provide new theoretical guidance for the pathogenesis of T. obscurus exposed to TBT. This study not only enriched the transcriptome data of T. obscurus under environmental stress, but also provided a new research method for the response mechanism of T. obscurus under the stimulation of environmental factors.

Genome editing of pufferfish saxitoxin- and tetrodotoxin-binding protein type 2 in Takifugu rubripes

The pufferfish saxitoxin- and tetrodotoxin-binding protein 2 (PSTBP2), which is involved in toxin accumulation, was knocked out in Takifugu rubripes embryos by using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 genome-editing technology. Treating the embryos with one of two single-guide RNA (sgRNA) resulted in mutation rates of 57.1% and 62.5%, respectively, as estimated using a heteroduplex mobility assay at 3 days postfertilization. Both sgRNAs might induced frameshift mutations that knocked out the T. rubripes PSTBP2.