Takifugu obscurus is a euryhaline fugu species very close to Takifugu rubripes and suitable for studying osmoregulation

Copper induced intestinal inflammation response through oxidative stress induced endoplasmic reticulum stress in Takifugu fasciatus

Copper (Cu) pollution in aquaculture water has seriously threatened the healthy and sustainable development of the aquaculture industry. Recently, many researchers have studied the toxic effects of Cu exposure on fish. However, the relationship between endoplasmic reticulum stress (ERS) and the inflammatory response, as well as its possible mechanisms, remain unclear. Particularly, information related to fish intestines must be expanded. Our study initially investigated the mechanisms underlying intestinal toxicity and inflammation resulting from Cu-induced ERS in vivo and in vitro in Takifugu fasciatus. In vivo study, T. fasciatus were treated with different concentrations (control, 20, and 100 µg/L) of Cu exposure for 28 days, causing intestinal oxidative stress, ERS, inflammatory responses, and histopathological and ultrastructural damage. Transcriptomic data further showed that Cu exposure caused ERS, as well as inflammatory responses, in the intestinal tracts of T. fasciatus. In vitro experiments on the intestinal cells of T. fasciatus showed that Cu exposure treatment (7.5 µg/mL) for 24 h induced ERS and increased mitochondrial numbers and inflammatory responses. In contrast, the addition of 4-phenylbutyric acid (4-PBA) alleviated ERS and inflammatory response in the Cu-exposed group. Furthermore, the reactive oxygen species (ROS) inhibitor, N-Acetyl-l-cysteine (NAC), effectively alleviated Cu-induced ERS. In conclusion, our in vivo and in vitro studies have confirmed that oxidative stress triggers the ERS pathway, which is involved in the intestinal inflammatory response. Our study provides new insights into the relationship among Cu-induced oxidative stress, ERS, and inflammatory responses in fish, as well as for the healthy culture of fish in aqueous environments.

Response of Intestinal Microbiota of Tiger Puffer (Takifugu rubripes) to the Fish Oil Finishing Strategy

The fish oil finishing (FOF) strategy, that is, re-feeding fish with fish oil (FO)-based diet after a certain period of feeding with alternative lipid source-based diets. On tiger puffer, the present study investigated the response of intestinal microbiota to FOF. Fish were fed four diets based on FO, soybean oil, palm oil and beef tallow as lipid sources, respectively, firstly for 50 days (growing-out period), and then fed the FO-based diet for 30 more days (FOF period). The results showed that dietary terrestrially sourced oils impaired the intestinal function in the growing-out period. However, the activities of amylase, trypsin and anti-oxidative enzymes (SOD, CAT, T-AOC), as well as gene expression of inflammatory cytokines (IL-1β, TNF-α, TGF-β) and tight junction protein (Claudin4, Claudin7, Claudin18, JAM, ZO-1) in the intestine were significantly recovered by FOF. The 16S rDNA sequencing analysis showed that FOF improved the similarity of bacterial community among the groups. The MetaStat analysis confirmed that FOF regulated the abundance of butyric acid-producing bacteria (Lachnospiraceae, Eubacterium, Butyricicoccus, Clostridium and Roseburia) and bacteria related to digestion and absorption (Sphingomonas, Romboutsia and Brevibacillus). In conclusion, FOF can recover the intestine function. The intestinal microbiota probably participated in and played a key role in the recovery process.

Seawater fish use an electrogenic boric acid transporter, Slc4a11A, for boric acid excretion by the kidney

Boric acid is a vital micronutrient in animals; however, excess amounts are toxic to them. Little is known about whole-body boric acid homeostasis in animals. Seawater (SW) contains 0.4 mM boric acid, and since marine fish drink SW, their urinary system was used here as a model of the boric acid excretion system. We determined that the bladder urine of a euryhaline pufferfish (river pufferfish, Takifugu obscurus) acclimated to fresh water and SW contained 0.020 and 19 mM of boric acid, respectively (a 950-fold difference), indicating the presence of a powerful excretory renal system for boric acid. Slc4a11 is a potential animal homolog of the plant boron transporter BOR1; however, mammalian Slc4a11 mediates H+ (OH-) conductance but does not transport boric acid. We found that renal expression of the pufferfish paralog of Slc4a11, Slc4a11A, was markedly induced after transfer from fresh water to SW, and Slc4a11A was localized to the apical membrane of kidney tubules. When pufferfish Slc4a11A was expressed in Xenopus oocytes, exposure to media containing boric acid and a voltage clamp elicited whole-cell outward currents, a marked increase in pHi, and increased boron content. In addition, the activity of Slc4a11A was independent of extracellular Na+. These results indicate that pufferfish Slc4a11A is an electrogenic boric acid transporter that functions as a B(OH)4- uniporter, B(OH)3-OH- cotransporter, or B(OH)3/H+ exchanger. These observations suggest that Slc4a11A is involved in the kidney tubular secretion of boric acid in SW fish, probably induced by the negative membrane potential and low pH of urine.

Effects of acute low-salinity stress on osmoregulation, antioxidant capacity, and growth of the black sea bream (Acanthopagrus schlegelii)

The black sea bream (Acanthopagrus schlegelii) is an important marine economic fish found on the southeast coast of China. Because of the frequent climate change, the salinity of the waters inhabited by A. schlegelii often decreases, which interferes with the fish's physiological homeostasis. The isotonic salinity of teleosts are usually lower than that of seawater, so maximum economic benefits cannot be obtained from conventional mariculture. This study was performed to preliminarily clarify the osmotic regulation and antioxidant mechanism of juvenile A. schlegelii and find an appropriate culture salinity value. We selected 5 psu, 10 psu, 15 psu, and 25 psu (control) to conduct physiological experiments for 96 h and growth experiments for 60 days. We found that the juvenile A. schlegelii could adjust their osmotic pressure within 12 h. The growth hormone and cortisol were found to be seawater-acclimating hormones, whereas prolactin was freshwater-acclimating hormone. The activity and mRNA expression of Na⁺/K⁺-ATPase showed a U-shaped trend with the decrease of in salinity at 12-96 h. Serum ion concentration and osmotic pressure remained at a relatively stable level after being actively adjusted from 6 to 12 h. At 96 h, the osmotic pressure of the serum isotonic point of juvenile A. schlegelii was approximately equal to that of water with 14.94 salinity. The number and volume of Cl^--secreting cells in the gills decreased. The glomeruli were more developed and structurally sound, with the renal tubules increasing in diameter and the medial brush border being more developed; this may indicate a decrease in salt secretion and an enhanced reabsorption function in the low salinity groups. The activities of superoxide dismutase and catalase and concentration of malondialdehyde were the lowest in the 15 psu group. In addition, the culture conditions of the 15 psu group improved the feed conversion rate without significant differences in weight gain when compared with the control group. Our results show that 15 psu salinity may be the best parameter for obtaining the maximum economic benefits.

Salinity and/or nanoparticles (Al2O3, TiO2) affect metal accumulation and ATPase activity in freshwater fish (Oreochromis niloticus)

The osmoregulation system of freshwater fish is sensitive to salinity increase in water. There is no satisfactory data to our knowledge on the accumulation of metal-oxide nanoparticles (NPs) in tissues of O. niloticus and their effects on ATPases (Na,K-ATPase, Mg-ATPase, Ca-ATPase) in differing salinities. Thus, this study investigated the effects of salinity (0 and 10 ppt) and Al₂O₃ and TiO₂ NPs (1 and 10 mg NPs/L) on the response of ATPases in acute (2 days) and chronic (20 days) durations. Data showed that nanoparticles accumulated in the tissues of fish, gill tissues having the highest levels of Al and Ti in both acute and chronic durations. Interestingly, the higher salinity significantly increased (P < 0.05) NP accumulations in the tissues in acute exposures, whereas it significantly decreased (P < 0.05) in chronic exposures. Salinity increase caused significant decreases (P < 0.05) in ATPase activities (up to 54 %) in control fish from both exposure protocols. Likewise, NP alone exposures (up to 80 %) and salt+NP (up to 83 %) exposures generally caused significant (P < 0.05) decreases in ATPase activities compared to their controls. Similarly, salt+NP exposures also decreased ATPase activities compared to NP exposures alone. The present data demonstrated that salinity and/or NP exposures decreased ATPase activities in the gill of freshwater fish, emphasizing the possible hazardous consequences of salt inputs and NP discharges into freshwater systems.

Membrane Transport Proteins Expressed in the Renal Tubular Epithelial Cells of Seawater and Freshwater Teleost Fishes

The kidney is an important organ that maintains body fluid homeostasis in seawater and freshwater teleost fishes. Seawater teleosts excrete sulfate and magnesium in small amounts of isotonic urine, and freshwater teleosts excrete water in large amounts of hypo-osmotic urine. The volume, osmolality, and ionic compositions of the urine are regulated mainly by membrane transport proteins expressed in the renal tubular epithelial cells. Gene expression, immunohistochemical, and functional analyses of the fish kidney identified membrane transport proteins involved in the secretion of sulfate and magnesium ions by the proximal tubules and reduction of urine volume by the collecting ducts in seawater teleosts, and excretion of water as hypotonic urine by the distal tubules and collecting ducts in freshwater teleosts. These studies promote an understanding of how the kidney contributes to the seawater and freshwater acclimation of teleosts at the molecular level.

An update on the evolutionary origin of aglomerular kidney with structural and ultrastructural descriptions of the kidney in three fish species

The kidney of fish contains numerous nephrons, each of which is divided into the renal corpuscle and renal tubules. This glomerular structure is the filtration unit of the nephron and is important for the kidney function, but it has been reported that the renal corpuscle was lost in at least four independent linages of fish (i.e., aglomerular kidney). In this study, the authors newly described renal structures for three species by histological and ultrastructural observations: two aglomerular kidneys from a seahorse Hippocampus barbouri and a toadfish Allenbatrachus grunniens and a glomerular kidney from a snake eel Pisodonophis boro. The renal development of H. barbouri was also described during 1-35 days after birth. In all species tested, the anterior kidney was comprised of haematopoietic tissues and a few renal tubules, whereas the posterior kidney contained more renal tubules. Although the glomerular structure was present in P. boro, light microscopic observations identified no glomeruli in the kidney of H. barbouri and A. grunniens. Ultrastructurally, abundant deep basal infoldings with mitochondria in the renal tubules were observed in A. grunniens compared to H. barbouri and P. boro, suggesting the possible role of basal infoldings in maintaining the osmotic balance. By integrating the results from the three species and comprehensive literature search, the authors further showed that 56 species have been reported to be aglomerular, and that the aglomerular kidney has evolved at least eight times in bony fishes.

Capacity for freshwater acclimation and differences in the transcription of ion transporter genes underlying different migratory life histories of Takifugu fish

The genus Takifugu is a group of approximately 20 species of puffer fishes living in a wide range of salinity environments around East Asian countries. This group presents a broad spectrum of evolutionary stages adapted to anadromy as a result of speciation that occurred a short time (2-5 million years) ago on an evolutionary timescale. This group thus can be considered as a model for studying the evolutionary mechanisms of anadromy. We firstly conducted a transfer experiment from seawater to low-salinity waters on five Takifugu species: two anadromous species T. obscurus and T. ocellatus, two euryhaline wanderer marine species T. rubripes and T. niphobles, and a strictly marine species T. snyderi, and confirmed that the capacity for acclimation to hypotonic environments was associated with their life history strategies. Next, transcriptomes of the gill and intestine of these species in hypotonic condition were compared to those under hypertonic condition for each species using RNA-Sequencing so as to determine possible candidate transporters playing an important role on freshwater adaptation. As this analysis suggested that cftr, encoding an important ion transporter for seawater acclimation in the gill, and ncc, encoding a transporter that is suggested to play important osmoregulatory roles in the intestine, are important candidates, their expression was validated by quantitative real-time PCR analysis. Expression of cftr was downregulated in the gills of the four euryhaline species under the hypotonic condition, but no change was detected in the gill of stenohaline T. snyderi, which may be one reason for the poor hypotonic acclimation capacity of T. snyderi. Expression of ncc was clearly upregulated in the intestines of the two anadromous species under the hypotonic condition, but not in other three species. Different ion transporter expression patterns between the five species indicate that the transcriptional regulation of cftr in the gill and ncc in the intestine may be important for the improvement of hypotonic acclimation capacity and evolution of anadromy in the Takifugu species.

Characterization of two myostatin genes in pufferfish Takifugu bimaculatus: sequence, genomic structure, and expression

Myostatin (MSTN) is a negative regulator of muscle growth, which restrains the proliferation and differentiation of myoblasts. To understand the role of two mstn genes of Takifugu bimaculatus, the full-length cDNAs of 1131 bp Tbmstn1 and 1,080 bp Tbmstn2 were obtained from the T. bimaculatus' genomic database, which encodes 376 and 359 amino acids, respectively. The results of qRT-PCR showed that Tbmstn1 was expressed in the eye, kidney, spleen, skeletal muscle, gill, and brain, and the expression level in the skeletal muscle was extremely significantly higher than in other examined tissues. Tbmstn2 was expressed in the skin, skeletal muscle, gill, and brain, and had the highest expression in the skeletal muscle, followed by expression in the brain. Meanwhile, in different stages of embryonic development, the expression of Tbmstn1 started from the gastrula stage. Its expression in the eye-pigment formation stage and hatching stage was significantly higher than that in other stages. The Tbmstn2 was expressed in all examined embryonic stages with different levels, and the highest expression was detected in the eye-pigment formation stage. These results suggested that Tbmstn1 and Tbmstn2 may involve in the development of skeletal muscle, and Tbmstn2 may be related to the formation of nervous system.

Transcriptome changes of Takifugu obscurus liver after acute exposure to the oxygenated-PAH 9,10-phenanthrenequione

Contamination with polycyclic aromatic hydrocarbons (PAHs) causes noticeable ecological problems in aquatic ecosystems. 9,10-Phenanthrenequione (9,10-PQ) is an oxidized PAH and is highly toxic to aquatic animals. However, the effects of 9,10-PQ on the molecular metabolism of fish remain largely unknown. In this study, Takifugu obscurus juveniles were acutely exposed to 44.30 µg/L 9,10-PQ for 3 days. The transcriptome profile changes in their livers were compared between the 9,10-PQ treatment group and the control using T. rubripes as the reference genome. The results identified 22,414 genes in our transcriptome. Among them, 767 genes were differentially expressed after exposure to 9,10-PQ, which enriched 16 KEGG pathways. Among them, the glycolysis, phagosome, and FOXO signaling pathways were significantly activated in 9,10-PQ treatment compared with the control. These data indicate that 9,10-PQ increased the glycolysis capacity to produce more energy for resistance and harmed immune function. Moreover, several genes related to tumorigenesis were significantly upregulated in response to 9,10-PQ, displaying the carcinogenic toxicity of 9,10-PQ to T. obscurus. Genes in steroid biosynthesis pathways were downregulated in the 9,10-PQ treatment group, suggesting interference with the endocrine system. Overall, these findings provide information to help evaluate the environmental risks that oxygenated-PAHs present to T. obscurus.

Immunohistochemical characterization and change in location of branchial ionocytes after transfer from freshwater to seawater in the euryhaline obscure puffer, Takifugu obscurus

The obscure puffer Takifugu obscurus is a euryhaline fish species suitable for studying the molecular mechanism of osmoregulation. The distributional changes of branchial ionocytes were detected following the transfer from freshwater (FW) to seawater (SW) based on two main ion transporters, Na⁺/K⁺-ATPase (NKA) and Na⁺/K⁺/ 2Cl^- cotransporter 1 (NKCC1). The mRNA and protein expression levels of NKA and NKCC1 in the gills all increased rapidly in the first four days after transfer to SW. Double immunofluorescence staining showed that NKCC1 and NKA were colocalized in the branchial ionocytes and the immunoreaction of NKCC1 was stronger after transfer. Moreover, following transfer to SW, the number of lamellar ionocytes in the gills is reduced and the number of filament ionocytes is increased significantly. Taken together, these findings indicated that SW transfer of obscure puffer promotes the changes of distribution, function and size of branchial ionocytes.

Effects of 9,10-phenanthrenequione on antioxidant indices and metabolite profiles in Takifugu obscurus plasma

Derived from polycyclic aromatic hydrocarbons (PAHs), oxygenated-PAHs (oxy-PAHs) may pose hazards to aquatic organisms, which remain largely unknown. Takifugu obscurus is an important anadromous fish species of high economic and ecological values. In the present study, T. obscurus was acutely exposed to 44.29 µg l-1 9,10-phenanthrenequione (9,10-PQ) for 96 h. Changes of antioxidant indices and metabolite profiles in plasma were compared between 9,10-PQ treatment and the control. The results showed that 9,10-PQ treatment significantly increased malondialdehyde (MDA) content during 6 to 96 h, increased superoxide dismutase (SOD) and catalase (CAT) activities at 6 h, but decreased them at 96 h. These results indicated that 9,10-PQ induced oxidative stress to fish. Ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) analysis revealed that four metabolic pathways were influenced in response to treatment with 9,10-PQ, including glycerophospholipid metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, purine metabolism and sulfur metabolism. These pathways are associated with antioxidant mechanisms, biosynthesis of neurotransmitters and innate immune functions. Thus, the as-obtained results confirmed that 9,10-PQ induced oxidative stress and raised concerns of neurotoxicity and immunotoxicity to fish. Overall, the present study posed a high environmental risk of oxy-PAHs to aquatic ecosystems.

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.

Histological, oxidative and immune changes in response to 9,10-phenanthrenequione, retene and phenanthrene in Takifugu obscurus liver

Polycyclic aromatic hydrocarbons (PAHs) are typical pollutants and may be alkylated and oxygenated to form alkyl-PAHs and oxygenated-PAHs (oxy-PAHs), respectively. Takifugu obscurus is an important anadromous fish species and displays a high risk of being exposed to PAHs-contaminated areas. In the present study, the effects of acute exposure to 44.29 µg L^-1 9,10-phenanthrenequione (9,10-PQ), retene and phenanthrene (Phe) on T. obscurus liver histology, antioxidant enzymes and immune indices were compared. After exposure to these three compounds, histological sections showed damages of hepatocyte, and the activities of alanine and aspartate aminotransferase increased in plasma, indicating direct hepatic toxicity. Hepatic malondialdehyde (MDA) content increased, but superoxide dismutase (SOD) and catalase (CAT) activities decreased in response to treatments with Phe, retene and 9,10-PQ. These results revealed peroxidative effects on T. obscurus hepatocytes. In addition, total immunoglobulin content and lysozyme activity in plasma increased in treatments with Phe, retene and 9,10-PQ, which might be resulted from the damaged liver cells and the subsequently hepatic inflammation. Besides, the changes were more severe in treatment with 9,10-PQ than those with Phe and retene, demonstrating higher toxicity of 9,10-PQ than the other two compounds. Overall, the present study posed a high environmental risk of PAH derivatives to aquatic ecosystems.

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.

Chromosomal-level assembly of Takifugu obscurus (Abe, 1949) genome using third-generation DNA sequencing and Hi-C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome-level reference genome for T. obscurus was constructed using long-read Pacific Biosciences (PacBio) Sequel sequencing and a Hi-C-based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single-molecule real-time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome-scale scaffolds using a Hi-C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one-to-one orthologous relationship was observed between the two species, supporting the chromosome-level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu-specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.

Evolutionary patterns of diadromy in fishes: more than a transitional state between marine and freshwater

Background

Across the tree of life there are numerous evolutionary transitions between different habitats (i.e., aquatic and terrestrial or marine and freshwater). Many of these dramatic evolutionary shifts parallel developmental shifts that require physiological, anatomical and behavioral changes for survival and reproduction. Diadromy (scheduled movement between marine and freshwater) has been characterized as a behavior that acts as an evolutionary intermediate state between marine and freshwater environments, implying that diadromous lineages are evolutionarily transient. This hypothesis comes with assumptions regarding the rates of evolutionary transitions in and out of diadromy as well as rates of speciation and extinction in diadromous fishes.

Results

Based on a published phylogeny of 7822 species of ray-finned fishes, state speciation and extinction models of evolutionary transition between marine, freshwater, and diadromous species suggest transition rates out of diadromy are 5–100 times higher that transition between marine and freshwater or into diadromy. Additionally, high speciation and low extinction rates separate diadromous fishes from marine and freshwater species. As a result, net diversification (net diversification = speciation – extinction) is about 7–40 times higher in diadromous fishes compared to freshwater and marine respectively. Together the transition, speciation, and extinction rates suggest diadromy is the least stable of the three states.

Conclusion

Evolutionary transitions to diadromy are rare in fishes. However, once established, diversification rates in diadromous lineages are high compared to both marine and freshwater species. Diadromous lineages tend to be more transient than marine or freshwater lineages and are found to give rise to marine and freshwater specialists in addition to diadromous descendants. Although diadromy is not a necessary evolutionary intermediate between marine and freshwater, these results support the interpretation of diadromy as an important, occasionally intermediate state, that contributes to biodiversity in fishes in all environments. This evolutionary instability of diadromous lineages is counteracted by their relatively high diversification rates. These findings highlight the importance of integrating the dynamics of diversification and major evolutionary transitions for understanding macroevolutionary patterns.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1492-2) contains supplementary material, which is available to authorized users.

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.

Individual and combined effects of salinity and lipopolysaccharides on the immune response of juvenile Takifugu fasciatus

Lipopolysaccharides (LPS) and salinity are important variables in aquatic environments. High concentration of LPS and large changes in salinity seriously threat the survival of a variety of organisms, including fish. To reveal the effects of salinity and LPS on a fish immune response, we measured the immune-related parameters (total leukocyte count, total serum protein, albumin and globulin concentrations, complement C3 concentration, and lysozyme activity) and genes (the expressions of TNF-α, IL-1β, and SOCS1-3 at the mRNA and protein levels) of juvenile Takifugu fasciatus exposed to phosphate buffered saline (PBS) or LPS (25 μg mL^-1) under different salinities (0, 15, and 30 ppt) for 24 h. Changes in key immunological indicators suggested that the LPS challenge induced considerable damage to T. fasciatus, whereas an increase in salinity mitigated the harmful effects. Moreover, although the immune responses in blood and other selected tissues (gill and kidney) were suppressed with an increase in salinity, the increased response in liver in saltwater enabled T. fasciatus to conquer large salinity variation during migration. The appropriate addition of salts appeared to be a sensible strategy to mitigate LPS-induced toxicity in the aquaculture of T. fasciatus.

iTRAQ-based quantitative proteomic analysis of Takifugu fasciatus liver in response to low-temperature stress

Low temperatures profoundly influence the physiological and behavioural processes of ectotherms, especially teleosts, which have made them the subjects of strong interest over time. However, the characteristics of fish cold-tolerance at the protein level remain unclear. Therefore, to shed further light on the molecular mechanisms of low temperature adaptation in fish, we conducted quantitative proteomics on the T. fasciatus liver using iTRAQ. Comparing the proteomic profiles of the T. fasciatus liver at 12 °C and 26 °C, a total of 3741 proteins were identified, and 160 were differentially abundant proteins (DAPs). Among the DAPs, the most significant changes were noted in proteins involved in oxidative stress (nine proteins), mitochondrial enzymes (eleven proteins) and signal transduction (thirteen proteins). The KEGG enrichment analysis indicated significant enhancement of D-arginine and D-ornithine metabolism, MAPK signalling, Wnt signalling and Gap junction pathway. Subsequently, three significantly up-regulated proteins (CIRB, HSP90 and GST) and two significantly down-regulated proteins (FLNB and A2ML1) were validated with parallel reaction monitoring (PRM) assays. Furthermore, the changes in abundance of proteins that are involved in oxidative stress, mitochondrial enzymes and signal transduction were validated at the transcriptional level with qPCR. These verification results show that the experimental data of iTRAQ are reliable. Our results not only deepen the understanding of the mechanisms underlying low-temperature tolerance in fish, but they also may contribute to the enhancement of cold tolerance during its breeding process. SIGNIFICANCE OF THE STUDY: The study focused on a comparative quantitative proteomics analysis of the T. fasciatus liver in response to low temperatures using iTRAQ, which has not yet been reported in the literatures. The results showed that the effect of low temperature on T. fasciatus is significant, including a detoxification of metabolic by-products and oxidative stress, an activation of the mitochondrial enzyme to strengthen energy metabolism, and a negative effect on signal transduction, which result in dysfunction or suboptimal performance. These low-temperature-related changes in the liver proteome of T. fasciatus can facilitate the understanding of the low temperature-related response that takes place in similar conditions in the liver and may contribute to the breeding of cold-resistant strains.

Effect of tributyltin chloride (TBT-Cl) exposure on expression of HSP90β1 in the river pufferfish (Takifugu obscurus): Evidences for its immunologic function involving in exploring process

HSP90β1 (known as glyco-protein 96, GP96) is a vital endoplasmic reticulum (ER) depended chaperonin among the HSPs (heat shock proteins) family. Furthermore, it always processes and presents antigen of the tumor and keeps balance for the intracellular environment. In the present study, we explored the effect of tributyltin chloride (TBT-Cl) exposure on HSP90β1 expression in river pufferfish, Takifugu obscurus. The full length of To-HSP90β1 was gained with 2775 bp in length, with an ORF (open reading frame) encoding an 803 aa polypeptide. A phylogenetic tree was constructed and showed the close relationship to other fish species. The HSP90β1 mRNA transcript was expressed in all tissues investigated with higher level in the gill and liver. After the acute and chronic exposure of TBT-Cl, the To-HSP90β1 mRNA transcript significantly was up-regulated in gills. Moreover, the histology study indicated the different injury degree of TBT-Cl in liver and gill. Immunohistochemistry (IHC) staining results implied the cytoplasm reorganization after TBT-Cl stress and the function of immunoregulation for To-HSP90β1 to TBT-Cl exposure. All the results indicated that HSP90β1 may be involved in the resistance to the invasion of TBT-Cl for keeping autoimmune homeostasis.

Molecular characterization and expression of suppressor of cytokine signaling (SOCS) 1, 2 and 3 under acute hypoxia and reoxygenation in pufferfish, Takifugu fasciatus

Hypoxia seriously affects the innate immune system of fish. However, the roles of suppressor of cytokine signaling (SOCS), pivotal anti-inflammatory genes, in response to hypoxia/reoxygenation remain largely unexplored. The primary objective of this study was to elucidate the function of SOCS genes under acute hypoxia and reoxygenation in pufferfish (Takifugu fasciatus). In the present study, SOCS1, 2 and 3 were identified in T. fasciatus referred to as TfSOCS1, 2 and 3. Then, qRT-PCR and western blot analysis were employed to assess their expressions at both the mRNA and protein levels. Tissue distribution demonstrated that the three SOCS genes were predominantly distributed in gill, brain and liver. Under hypoxia challenge (1.63 ± 0.2 mg/L DO for 2, 4, 6 and 8 h), the expressions of TfSOCS1 and 3 in brain and liver at the mRNA and protein levels were significantly decreased, while their expressions showed an opposite trend in gill. Different from the expressions of TfSOCS1 and 3, the expression of TfSOCS2 was inhibited in gill, along with its increased expression in brain and liver. After normoxic recovery (7.0 ± 0.3 mg/L of DO for 4 and 12 h), most of TfSOCS genes were significantly altered at R4 (reoxygenation for 4 h) and returned to the normal level at R12 (reoxygenation for 12 h). SOCS genes played vital roles in response to hypoxia/reoxygenation challenge. Our findings greatly strengthened the relation between innate immune and hypoxia stress in T. fasciatus.

Characterization and expression of Na+/K+-ATPase in gills and kidneys of the Teleost fish Oreochromis mossambicus, Oreochromis urolepis hornorum and their hybrids in response to salinity challenge

Tilapia (Oreochromis mossambicus, O. urolepis hornorum, their hybrids O. mossambicus♀ × O. hornorum♂ and O. hornorum♀ × O. mossambicus♂) were exposed to a high salinity environment to evaluate their osmoregulatory responses. The plasma osmolality of all the tilapia species were elevated with the salinity challenge. The activities of Na⁺/K⁺-ATPase (NKA) in both the gill and kidney showed a similar increased change tendency compared with the control. The distribution of NKA α1 mRNA in all the examined tissues suggested that NKA α1 has a possible housekeeping role for this isoform. The amount of NKA α1 mRNA in the gill and kidney was elevated in the four fishes with similar expression patterns after transfer from freshwater to seawater. The NKAα1 mRNA expression levels in the gill reached their peak level at 24 h after transfer (P < 0.01) compared to the freshwater group, following decreases in the pretreatment level at 48 h (P > 0.05). However, the NKAα1 mRNA expression levels in the kidney were not significantly affected with increasing environmental salinity (P > 0.05). The differences in the responses to saltwater challenge may be associated with differences in saltwater tolerance between the four tilapia. The drastic increase in the plasma osmolality, NKA activities and mRNA expression suggested that the hybrids (O. mossambicus♀ × O. hornorum♂) possess heterosis in salinity responsiveness compared to that of both the parents, indicating a maternal effect on the salinity tolerance of the tilapia hybrids. This study provides a theoretical basis to further study the mechanism of fish osmoregulation response to salinity challenge.

Physiological mechanism of osmoregulatory adaptation in anguillid eels

In recent years, the production of eel larvae has dramatic declines due to reductions in spawning stocks, overfishing, growth habitat destruction and access reductions, and pollution. Therefore, it is particularly important and urgent for artificial production of glass eels. However, the technique of artificial hatching and rearing larvae is still immature, which has long been regarded as an extremely difficult task. One of the huge gaps is artificial condition which is far from the natural condition to develop their capability of osmoregulation. Thus, understanding their osmoregulatory mechanisms will help to improve the breed and adapt to the changes in the environment. In this paper, we give a general review for a study progress of osmoregulatory mechanisms in eels from five aspects including tissues and organs, ion transporters, hormones, proteins, and high throughput sequencing methods.

Combined effects of cadmium and salinity on juvenile Takifugu obscurus: cadmium moderates salinity tolerance; salinity decreases the toxicity of cadmium

Obscure puffer Takifugu obscurus, a species of anadromous fish, experiences several salinity changes in its lifetime. Cadmium (Cd) is a toxic heavy metal that can potentially induce oxidative stress in fish. The present study aimed to detect the combined effects of Cd (0, 5, 10, 20 and 50 mg L(-1)) and salinity (0, 15 and 30 ppt) on juvenile T. obscurus. Results showed the juveniles could survive well under different salinities; however, with Cd exposure, the survival rates significantly decreased at 0 and 30 ppt. At 15 ppt, tolerance to Cd increased. Cd exposure clearly induced oxidative stress, and the responses among different tissues were qualitatively similar. Salinity acted as a protective factor which could reduce the reactive oxygen species and malondialdehyde levels. In addition, salinity could enhance the antioxidant defense system, including superoxide dismutase, catalase and glutathione. Na(+)/K(+)-ATPase activity significantly decreased under Cd exposure in gill, kidney and intestine. These findings indicated that Cd could moderate the adaptability of juvenile T. obscurus to high salinity and low salinity played a protective role upon Cd exposure. Thus, the role of salinity should be considered when evaluating the effect of heavy metals on anadromous and estuarine fishes.

Responses of Takifugu obscurus fertilized eggs and larvae to increased ammonia exposure

Ammonia is a common toxicant in aquatic systems; this substance has become a critical threat to fish, especially in early life stages. This study aimed to evaluate the effects of unionized ammonia (NH3-N: 0, 0.068, 0.138, 0.206, 0.275, 0.343, 0.412, and 0.481 mg L(-1)) on fertilized eggs and larvae of obscure puffer Takifugu obscurus, a fish species with potential economic value. Results showed that hatch time was significantly retarded and hatch rate was significantly decreased as NH3-N concentrations increased; newly hatched larvae exhibited high rate of abnormalities and low viability. The survival rate of larvae also decreased significantly as NH3-N concentrations increased; larvae could tolerate NH3-N to a less extent than embryos. NH3-N also caused a significant decrease in superoxide dismutase (SOD) and Na(+)/K(+) ATPase activities but not in malondialdehyde (MDA) levels of larvae. Two-way ANOVA indicated that there was a statistically significant interaction between NH3-N concentrations and exposure times on SOD activity but not on Na(+)/K(+) ATPase activity. Such responses indicated that an increase in ammonia concentration in surface water may negatively affect the early development of T. obscurus and thus likely impair population recruitment and persistence of this fish species.

Hypo-osmotic stress-induced physiological and ion-osmoregulatory responses in European sea bass (Dicentrarchus labrax) are modulated differentially by nutritional status

We investigated the impact of nutritional status on the physiological, metabolic and ion-osmoregulatory performance of European sea bass (Dicentrarchus labrax) when acclimated to seawater (32 ppt), brackish water (20 and 10 ppt) and hyposaline water (2.5 ppt) for 2 weeks. Following acclimation to different salinities, fish were either fed or fasted (unfed for 14 days). Plasma osmolality, [Na(+)], [Cl(-)] and muscle water content were severely altered in fasted fish acclimated to 10 and 2.5 ppt in comparison to normal seawater-acclimated fish, suggesting ion regulation and acid-base balance disturbances. In contrast to feed-deprived fish, fed fish were able to avoid osmotic perturbation more effectively. This was accompanied by an increase in Na(+)/K(+)-ATPase expression and activity, transitory activation of H(+)-ATPase (only at 2.5 ppt) and down-regulation of Na(+)/K(+)/2Cl(-) gene expression. Ammonia excretion rate was inhibited to a larger extent in fasted fish acclimated to low salinities while fed fish were able to excrete efficiently. Consequently, the build-up of ammonia in the plasma of fed fish was relatively lower. Energy stores, especially glycogen and lipid, dropped in the fasted fish at low salinities and progression towards the anaerobic metabolic pathway became evident by an increase in plasma lactate level. Overall, the results indicate no osmotic stress in both feeding treatments within the salinity range of 32 to 20 ppt. However, at lower salinities (10-2.5 ppt) feed deprivation tends to reduce physiological, metabolic, ion-osmo-regulatory and molecular compensatory mechanisms and thus limits the fish's abilities to adapt to a hypo-osmotic environment.

Identification and lateral membrane localization of cyclin M3, likely to be involved in renal Mg2+ handling in seawater fish

The kidney of marine teleosts is the major site of Mg(2+) excretion and produces urine with a high Mg(2+) concentration. However, the transporters involved in Mg(2+) excretion are poorly understood. The cyclin M (Cnnm; also known as ancient conserved domain protein) family comprises membrane proteins homologous to the bacterial Mg(2+) and Co(2+) efflux protein, CorC. To understand the molecular mechanism of Mg(2+) homeostasis in marine teleosts, we analyzed the expression of the Cnnm family genes in the seawater (SW) pufferfish, torafugu (Takifugu rubripes), and the closely related euryhaline species, mefugu (Takifugu obscurus). Database mining and phylogenetic analysis indicated that the Takifugu genome contains six members of the Cnnm family: two orthologs of Cnnm1, one of Cnnm2, one of Cnnm3, and two of Cnnm4. RT-PCR analyses indicated that Cnnm2, Cnnm3, and Cnnm4a are expressed in the kidney, whereas other members are mainly expressed in the brain. Renal expression of Cnnm3 was upregulated in SW mefugu, whereas renal expression of Cnnm2 was upregulated in freshwater (FW) mefugu. No significant difference was observed in renal expression of Cnnm4a between SW and FW mefugu. In situ hybridization and immunohistochemical analyses of the SW mefugu kidney revealed that Cnnm3 is expressed in the proximal tubule, and its product localizes to the lateral membrane. When Cnnm3 was expressed in Xenopus laevis oocytes, whole cellular Mg(2+) content and free intracellular Mg(2+) activity significantly decreased. These results suggest that Cnnm3 is involved in body fluid Mg(2+) homeostasis in marine teleosts.

Acclimation of brackish water pearl spot (Etroplus suratensis) to various salinities: relative changes in abundance of branchial Na(+)/K (+)-ATPase and Na (+)/K (+)/2Cl (-) co-transporter in relation to osmoregulatory parameters

The present study was conducted to elucidate the osmoregulatory ability of the fish pearl spot (Etroplus suratensis) to know the scope of this species for aquaculture under various salinities. Juvenile pearl spot were divided into three groups and acclimated to freshwater (FW), brackish water (BW) or seawater (SW) for 15 days. The fish exhibited effective salinity tolerance under osmotic challenges. Although the plasma osmolality and Na(+), K(+) and Cl(-) levels increased with the increasing salinities, the parameters remained within the physiological range. The muscle water contents were constant among FW-, BW- and SW-acclimated fish. Two Na+/K+-ATPase α-isoforms (NKA α) were expressed in gills during acclimation in FW, BW and SW. Abundance of one isoform was up-regulated in response to seawater acclimation, suggesting its role in ion secretion similar to NKA α1b, while expression of another isoform was simultaneously up-regulated in response to both FW and SW acclimation, suggesting the presence of isoforms switching phenomenon during acclimation to different salinities. Nevertheless, NKA enzyme activities in the gills of the SW and FW individuals were higher (p < 0.05) than in BW counterparts. Immunohistochemistry revealed that Na(+)/K(+)-ATPase immunoreactive (NKA-IR) cells were mainly distributed in the interlamellar region of the gill filaments in FW groups and in the apical portion of the filaments in BW and SW groups. The number of NKA-IR cells in the gills of the FW-acclimated fish was almost similar to that of SW individuals, which exceeded that of the BW individuals. The NKA-IR cells of BW and SW were bigger in size than their FW counterparts. Besides, the relative abundance of branchial Na(+)/K(+)/2Cl(-) co-transporter showed stronger evidence in favor of involvement of this protein in hypo-osmoregulation, requiring ion secretion by the chloride cells. To the best of our knowledge, this is the first study reporting the wide salinity tolerance of E. suratensis involving differential activation of ion transporters and thereby suggesting its potential as candidate for fish farming under different external salinities.

Salinity changes in the anadromous river pufferfish, Takifugu obscurus, mediate gene regulation

This study aimed to better understand the hydromineral regulatory response of the anadromous river pufferfish, Takifugu obscurus, to salinity changes through real-time RT-PCR. After abrupt transfer from 30 or 5 psu to 5 or 30 psu, respectively, we analyzed the mRNA expression of Na⁺/K⁺ ATPase, prolactin receptor, and aquaporin from osmoregulatory organs of the river pufferfish such as gills, kidney, and intestine. Na⁺/K⁺ ATPase showed notable changes in the gills and kidney when salinity was increased. In the gills, the expression level of Na⁺/K⁺ ATPase suddenly increased within a day after abrupt transfer from 5 to 30 psu and then slightly declined within 2 days after exposure. In the kidney, Na⁺/K⁺ ATPase has shown consistently high mRNA expression after the increase in salinity. Expression levels of the prolactin receptor gene increased when environmental salinity decreased. In the intestine, gene expression of the prolactin receptor remained high, even when salinity decreased. To the contrary, there was a steady increase or decrease in mRNA expression in the kidney in response to salinity decrease or increase, respectively. As for aquaporins, aquaporin 1 was mainly expressed in the intestine and kidney, and aquaporin 3 was mainly expressed in the gills and intestine. In the gills, increased expression of aquaporin 3 was found after transfer to lower salinity and in the intestine and kidney, a decrease in salinity followed by an abrupt decrease in aquaporin 1 and aquaporin 3. Contrastingly, the expression of these genes increased in the intestine after transfer to 30 psu. Osmoregulatory genes were expressed in diverse organs, apparently to overcome an influx or exhaust of water or ions. A superior adaptation ability of the river pufferfish to a wide range of salinities is most reasonably due to active osmoregulatory processes mediated by the genes monitored here.

Identification and proximal tubular localization of the Mg²⁺ transporter, Slc41a1, in a seawater fish

The second most abundant cation in seawater (SW), Mg²⁺, is present at concentrations of ~53 mM. Marine teleosts maintain plasma Mg²⁺ concentration at 1-2 mM by excreting Mg²⁺ into the urine. Urine Mg²⁺ concentrations of SW teleosts exceed 70 mM, most of which is secreted by the renal tubular epithelial cells. However, molecular mechanisms of the Mg²⁺ secretion have yet to be clarified. To identify transporters involved in Mg²⁺ secretion, we analyzed the expression of fish homologs of the Slc41 Mg²⁺ transporter family in various tissues of SW pufferfish torafugu (Takifugu rubripes) and its closely related euryhaline species mefugu (Takifugu obscurus). Takifugu genome contained five members of Slc41 genes, and only Slc41a1 was highly expressed in the kidney. Renal expression of Slc41a1 was markedly elevated when mefugu were transferred from fresh water (FW) to SW. In situ hybridization analysis and immunohistochemistry at the light and electron microscopic levels revealed that Slc41a1 is localized to vacuoles in the apical cytoplasm of the proximal tubules. These results suggest that pufferfish Slc41a1 is a Mg²⁺ transporter involved in renal tubular transepithelial Mg²⁺ secretion by mediating Mg²⁺ transport from the cytosol to the vacuolar lumen, and support the hypothesis that Mg²⁺ secretion is mediated by exocytosis of Mg²⁺-rich vacuoles to the lumen.

Cloning, expression, and characterization of TNFSF14 (LIGHT) gene in mefugu, Takifugu obscurus

LIGHT/TNFSF14 is a member of the tumor necrosis factor ligand superfamily, which plays important roles in inflammatory and immune responses. In this study, the cDNA of mefugu (Takifugu obscures) LIGHT (designated as fLIGHT) was amplified from spleen by reverse transcription polymerase chain reaction. The open reading frame of fLIGHT covers 765 bp and translates into a 254-aa protein. The predicted three-dimensional (3D) structure of the soluble LIGHT of mefugu (designated as fsLIGHT) monomer analyzed by comparative protein modeling revealed that it was very similar to its human counterpart. Real-time quantitative PCR analysis indicated that LIGHT is constitutively expressed in various tissues in mefugu. The soluble LIGHT binding of green fluorescent protein (GFP) (designated as GFP/fsLIGHT) had been cloned into the pET28a vector; SDS-PAGE and western blotting analysis confirmed that the recombinant protein pET28a-GFP/fsLIGHT was efficiently expressed in Escherichia coli BL21 (DE3). After purification, the GFP/fsLIGHT fusion protein obtained similar fluorescence spectrum to those of GFP. Laser scanning confocal microscopy analysis showed GFP/fsLIGHT could bind to its receptors. In view of our basic research, LIGHT may be a potential immunologic factor for enhancing immunological efficacy in fish, and our results might provide a platform for further study into the effects of LIGHT.

Expression of sodium/hydrogen exchanger 3 and cation-chloride cotransporters in the kidney of Japanese eel acclimated to a wide range of salinities

Reabsorption of monovalent ions in the kidney is essential for adaptation to freshwater and seawater in teleosts. To assess a possible role of Na(+)/H(+) exchanger 3 (NHE3) in renal osmoregulation, we first identified a partial sequence of cDNA encoding NHE3 from the Japanese eel kidney. For comparison, we also identified cDNAs encoding kidney specific Na(+)-K(+)-2Cl(-) cotransporter 2 (NKCC2α) and Na(+)-Cl(-) cotransporter (NCCα). In eels acclimated to a wide range of salinities from deionized freshwater to full-strength seawater, the expression of NHE3 in the kidney was the highest in eel acclimated to full-strength seawater. Meanwhile, the NCCα expression exhibited a tendency to increase as the environmental salinity decreased, whereas the NKCC2α expression was not significantly different among the experimental groups. Immunohistochemical studies showed that NHE3 was localized to the apical membrane of epithelial cells composing the second segments of the proximal renal tubule in seawater-acclimated eel. Meanwhile, the apical membranes of epithelial cells in the distal renal tubule and collecting duct showed more intense immunoreactions of NKCC2α and NCCα, respectively, in freshwater eel than in seawater eel. These findings suggest that renal monovalent-ion reabsorption is mainly mediated by NKCC2α and NCCα in freshwater eel and by NHE3 in seawater eel.

Identification and apical membrane localization of an electrogenic Na⁺/Ca²⁺ exchanger NCX2a likely to be involved in renal Ca²⁺ excretion by seawater fish

Seawater (SW) contains ∼10 mM Ca(2+), yet marine fish must drink seawater as their major water source. Thus marine teleosts fish need to excrete Ca(2+) to maintain whole body Ca(2+) homeostasis. In the intestine, seawater Ca(2+) interreacts with epithelial-secreted HCO(3)(-) by the intestinal epithelium, and the resulting CaCO(3) precipitates, which is rectally excreted. Recently the transporters involved in intestinal HCO(3)(-) secretion were identified. Ca(2+) is also excreted by the kidney, but the protein(s) involved in renal Ca(2+) excretion have not been identified. Here we identified a candidate transporter by using SW pufferfish torafugu (Takifugu rubripes) and its closely related euryhaline species mefugu (Takifugu obscurus), which are becoming useful animal models for studying molecular mechanisms of seawater adaptation. RT-PCR analyses of Na(+)/Ca(2+) exchanger (NCX) family members in various torafugu tissues demonstrated that only NCX2a is highly expressed in the kidney. Renal expression of NCX2a was markedly elevated when mefugu were transferred from freshwater to seawater. In situ hybridization and immunohistochemical analyses indicated that NCX2a is expressed in the proximal tubule at the apical membrane. NCX2a, expressed in Xenopus oocytes, conferred [Ca(2+)](out)- and Na(+)-dependent currents. These results suggest that NCX2a mediates renal Ca(2+) secretion at the apical membrane of renal proximal tubules and has an important role in whole body Ca(2+) homeostasis of marine teleosts.

Molecular structure, expression and bioactivity characterization of TNF13B (BAFF) gene in mefugu, Takifugu obscurus

B-cell activating factor of the TNF family (BAFF) induces B-cell survival, proliferation, immunoglobulin secretion and plays a role in enhancing immune responses. In the present study, a BAFF homolog has been identified in mefugu Takifugu obscurus, and its biological activities have been characterized. The mefugu BAFF (fBAFF) cDNA is 789 bp in length and translates into a 262-aa protein. The fBAFF genomic sequence consists of six exons and five introns, is approximately 1.8 kb in size. Amino acid sequence comparison indicated that fBAFF possessed the TNF signatures, a transmembrane domain, a putative furin protease cleavage site and three cysteine residues, which were the typical characteristics of TNF gene in mammals and birds. The predicted three-dimensional (3D) structure of the fsBAFF monomer analyzed by comparative protein modeling revealed that it was very similar to its human counterpart. Real-time quantitative PCR (qPCR) analysis revealed that fBAFF was predominantly expressed in mefugu lymphoid tissue spleen. The SUMO-fsBAFF and GFP/fsBAFF efficiently expressed in Escherichia coli Rosetta (DE3) were confirmed by SDS-PAGE and Western blotting analysis. After purification, the GFP/fsBAFF fusion protein obtained similar fluorescence spectrum to those of GFP. Laser scanning confocal microscopy analysis showed GFP/fsBAFF could bind to its receptors. In vitro, the MTT assays and flow cytometric analysis indicated that SUMO-fsBAFF could promote the proliferation of mefugu splenocytes or mouse splenic B cells together with/without anti-mouse IgM. These findings indicate that SUMO-fsBAFF plays an important role in proliferation of mefugu B cells and has functional cross-reactivity among mefugu and other mammalians. Therefore, BAFF may be a potential immunologic factor for enhancing immunological efficacy in fish.

Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals

The compact genome of fugu (Takifugu rubripes) has been used widely as a reference genome for understanding the evolution of vertebrate genomes. However, the fragmented nature of the fugu genome assembly has restricted its use for comparisons of genome architecture in vertebrates. To extend the contiguity of the assembly to the chromosomal level, we have generated a comprehensive genetic map of fugu and anchored the scaffolds of the assembly to the 22 chromosomes of fugu. The map consists of 1,220 microsatellite markers that provide anchor points to 697 scaffolds covering 86% of the genome assembly (http://www.fugu-sg.org/). The integrated genome map revealed a higher recombination rate in fugu compared with other vertebrates and a wide variation in the recombination rate between sexes and across chromosomes of fugu. We used the extended assembly to explore recent rearrangement events in the lineages of fugu, Tetraodon, and medaka and compared them with rearrangements in three mammalian (human, mouse, and opossum) lineages. Between the two pufferfishes, fugu has experienced fewer chromosomal rearrangements than Tetraodon. The gene order is more highly conserved in the three teleosts than in mammals largely due to a lower rate of interchromosomal rearrangements in the teleosts. These results provide new insights into the distinct patterns of genome evolution between teleosts and mammals. The consolidated genome map and the genetic map of fugu are valuable resources for comparative genomics of vertebrates and for elucidating the genetic basis of the phenotypic diversity of ~25 species of Takifugu that evolved within the last 5 My.

Multiple invasions into freshwater by pufferfishes (teleostei: tetraodontidae): a mitogenomic perspective

Pufferfishes of the Family Tetraodontidae are the most speciose group in the Order Tetraodontiformes and mainly inhabit coastal waters along continents. Although no members of other tetraodontiform families have fully discarded their marine lives, approximately 30 tetraodontid species spend their entire lives in freshwaters in disjunct tropical regions of South America, Central Africa, and Southeast Asia. To investigate the interrelationships of tetraodontid pufferfishes and thereby elucidate the evolutionary origins of their freshwater habitats, we performed phylogenetic analysis based on whole mitochondrial genome sequences from 50 tetraodontid species and closely related species (including 31 newly determined sequences). The resulting phylogenies reveal that the family is composed of four major lineages and that freshwater species from the different continents are independently nested in two of the four lineages. A monophyletic origin of the use of freshwater habitats was statistically rejected, and ancestral habitat reconstruction on the resulting tree demonstrates that tetraodontids independently entered freshwater habitats in different continents at least three times. Relaxed molecular-clock Bayesian divergence time estimation suggests that the timing of these invasions differs between continents, occurring at 0–10 million years ago (MA) in South America, 17–38 MA in Central Africa, and 48–78 MA in Southeast Asia. These timings are congruent with geological events that could facilitate adaptation to freshwater habitats in each continent.

Differential expression of Na+-Cl- cotransporter and Na+-K+-Cl- cotransporter 2 in the distal nephrons of euryhaline and seawater pufferfishes

The process of NaCl reabsorption in the distal nephron allows freshwater fishes to excrete hypotonic urine and seawater fishes to excrete urine containing high concentrations of divalent ions; the relevant transporters, however, have not yet been identified. In the mammalian distal nephron, NaCl absorption is mediated by Na(+)-K(+)-Cl(-) cotransporter 2 (NKCC2, Slc12a1) in the thick ascending limb, Na(+)-Cl(-) cotransporter (NCC, Slc12a3) in the distal convoluted tubule, and epithelial sodium channel (ENaC) in the collecting duct. In this study, we compared the expression profiles of these proteins in the kidneys of euryhaline and seawater pufferfishes. Mining the fugu genome identified one NKCC2 gene and one NCC gene, but no ENaC gene. RT-PCR and in situ hybridization analyses demonstrated that NKCC2 was highly expressed in the distal tubules and NCC was highly expressed in the collecting ducts of euryhaline pufferfish (mefugu, Takifugu obscurus). On the other hand, the kidney of seawater pufferfish (torafugu, Takifugu rubripes), which lacked distal tubules, expressed very low levels of NCC, and, in the collecting ducts, high levels of NKCC2. Acclimation of mefugu to seawater resulted in a 2.7× decrease in NCC expression, whereas NKCC2 expression was not markedly affected. Additionally, internalization of NCC from the apical surface of the collecting ducts was observed. These results suggest that NaCl reabsorption in the distal nephron of the fish kidney is mediated by NCC and NKCC2 in freshwater and by NKCC2 in seawater.

Effect of cadmium exposure on expression of antioxidant gene transcripts in the river pufferfish, Takifugu obscurus (Tetraodontiformes)

Cadmium (Cd) is a non-essential toxic heavy metal with the potential to induce oxidative stress. Cd toxicity and its capacity for accumulation in aquatic habitats have earned its recognition as a pollutant of immediate and widespread concern. To obtain a better understanding of oxidative stress-associated gene expression in different tissues, six antioxidant genes such as catalase (CAT), glutathione reductase (GR), glutathione peroxidase 1a (GPx1a), glutathione peroxidase 1b (GPx1b), Cu/Zn superoxide dismutase (Cu/Zn-SOD), and Mn superoxide dismutase (Mn-SOD) were cloned and fully sequenced in the river pufferfish, Takifugu obscurus. On tissue specific mRNA expression, the liver showed the highest expression when compared to other tissues, even though each antioxidant gene showed different modes of expression patterns in the examined tissues. Of the various antioxidant genes, GR was the most highly expressed in the liver, followed by CAT, GPx1, and Cu/Zn-SOD. For the time-course experiment, all the antioxidant genes were significantly induced over time except for Cu/Zn-SOD in the liver, and there was a 5-fold induction in hepatic GR, CAT, and Mn-SOD mRNA compared to the control. These findings indicate that the liver of T. obscurus has a robust antioxidant system. In addition, these results suggest that Cd exposure modulates the expression of antioxidant genes, and would indicate that the antioxidant genes would be a relevant biomarker of trace metal pollution such as Cd exposure in T. obscurus.

Co-ordination of osmotic stress responses through osmosensing and signal transduction events in fishes

This review centres upon the molecular regulation of osmotic stress responses in fishes, focusing on how osmosensing and signal transduction events co-ordinate changes in the activity and abundance of effector proteins during osmotic stress and how these events integrate into osmotic stress responses of varying magnitude. The concluding sections discuss the relevance of osmosensory signal transduction to the evolution of euryhalinity and present experimental approaches that may best stimulate future research. Iterating the importance of osmosensing and signal transduction during fish osmoregulation may be pertinent amidst the increased use of genomic technologies that typically focus solely on changes in the abundances of gene products, and may limit insight into critical upstream events that occur mainly through post-translational mechanisms.

Proteomic analysis on the alteration of protein expression in gills of ayu (Plecoglossus altivelis) associated with salinity change

Gill is the primary osmoregulatory organ for euryhaline fish to acclimate salinity change. The effect of salinity on gill proteome in ayu, Plecoglossus altivelis, was investigated by two-dimensional gel electrophoresis (2-DE) and matrix assisted laser desorption ionization time-of-flight tandem mass spectrometry (MALDI-TOF-MS/MS). Eight of eighteen altered proteins were successfully identified. They are involved in osmoregulation, cytoskeleton, energy metabolism, and stress response. Our results showed that vinculin, echinoderm microtubule-associated protein like protein 1, pyruvate kinase, betaine-homocysteine methyltransferase (BHMT), transaldolase, glyceraldehyde 3-phosphate dehydrogenase, and heat shock protein 70 (HSP70) were down-regulated, whereas cofilin was up-regulated when ayu transferred from fresh water (FW) to brackish water (BW). Partial cDNA sequences of BHMT, HSP70, Na(+)/K(+) ATPase (NKA) alpha-subunit and 18S rRNA genes were subsequently determined and used for 2-DE data verification by real-time PCR. Gill BHMT and HSP70 mRNAs decreased significantly in BW-transferred ayu, while NKA alpha-subunit mRNA had no significant change. It was suggested that cell volume-regulatory response, especially the protection by the BHMT/betaine system might play an important role in ayu acclimation to salinity change.