Intestinal transcriptome analysis revealed differential salinity adaptation between two tilapiine species

Biochemical indices, gene expression, and SNPs associated with salinity adaptation in juvenile chum salmon (Oncorhynchus keta) as determined by comparative transcriptome analysis

Chum salmon (Oncorhynchus keta) migrate from freshwater to saltwater, and incur developmental, physiological and molecular adaptations as the salinity changes. The molecular regulation for salinity adaptation in chum salmon is currently not well defined. In this study, 1-g salmon were cultured under 0 (control group, D0), 8‰ (D8), 16‰ (D16), and 24‰ (D24) salinity conditions for 42 days. Na⁺/K⁺-ATPase and Ca²⁺/Mg²⁺-ATPase activities in the gill first increased and then decreased in response to higher salinity environments where D8 exhibited the highest Na⁺/K⁺ATPase and Ca²⁺/Mg²⁺-ATPase activity and D24 exhibited the lowest. Alkaline phosphatase (AKP) activity was elevated in all salinity treatment groups relative to controls, while no significant difference in acid phosphatase (ACP) activity was observed across treatment groups. De novo transcriptome sequencing in the D0 and D24 groups using RNA-Seq analysis identified 187,836 unigenes, of which 2,143 were differentially expressed in response to environmental salinity (71 up-regulated and 2,072 down-regulated). A total of 56,020 putative single nucleotide polymorphisms (SNPs) were also identified. The growth, development, osmoregulation and maturation factors of N-methyl-D-aspartate receptors (nmdas) expressed in memory formation, as well as insulin-like growth factor 1 (igf-1) and igf-binding proteins (igfbps) were further investigated using targeted qRT-PCR. The lowest expression of all these genes occurred in the low salinity environments (D8 or D16), while their highest expression occurred in the high salinity environments (D24). These results provide preliminary insight into salinity adaptation in chum salmon and a foundation for the development of marker-assisted breeding for this species.

Salinity significantly affects intestinal microbiota and gene expression in striped catfish juveniles

In the present study, juvenile striped catfish (Pangasianodon hypophthalmus), a freshwater fish species, have been chronically exposed to a salinity gradient from freshwater to 20 psu (practical salinity unit) and were sampled at the beginning (D20) and the end (D34) of exposure. The results revealed that the intestinal microbial profile of striped catfish reared in freshwater conditions were dominated by the phyla Bacteroidetes, Firmicutes, Proteobacteria, and Verrucomicrobia. Alpha diversity measures (observed OTUs (operational taxonomic units), Shannon and Faith's PD (phylogenetic diversity)) showed a decreasing pattern as the salinities increased, except for the phylogenetic diversity at D34, which was showing an opposite trend. Furthermore, the beta diversity between groups was significantly different. Vibrio and Akkermansia genera were affected differentially with increasing salinity, the former being increased while the latter was decreased. The genus Sulfurospirillium was found predominantly in fish submitted to salinity treatments. Regarding the host response, the fish intestine likely contributed to osmoregulation by modifying the expression of osmoregulatory genes such as nka1a, nka1b, slc12a1, slc12a2, cftr, and aqp1, especially in fish exposed to 15 and 20 psu. The expression of heat shock proteins (hsp) hsp60, hsp70, and hsp90 was significantly increased in fish reared in 15 and 20 psu. On the other hand, the expression of pattern recognition receptors (PRRs) were inhibited in fish exposed to 20 psu at D20. In conclusion, the fish intestinal microbiota was significantly disrupted in salinities higher than 10 psu and these effects were proportional to the exposure time. In addition, the modifications of intestinal gene expression related to ion exchange and stressful responses may help the fish to adapt hyperosmotic environment. KEY POINTS: • It is the first study to provide detailed information on the gut microbiota of fish using the amplicon sequencing method. • Salinity environment significantly modified the intestinal microbiota of striped catfish. • Intestinal responses may help the fish adapt to hyperosmotic environment.

Salinity affects growth performance, physiology, immune responses and temperature resistance in striped catfish (Pangasianodon hypophthalmus) during its early life stages

In this study, striped catfish larvae were gradually exposed to the increase of different salinities, and then they reached the levels of 0, 5, 10, 15, and 20 psu after 10 days, followed by heat shock at 39 °C to determine stress tolerance. After the 10-day experiment, the survival rate of fish exposed to the 20 psu treatment was only 28.6 ± 4%, significantly lower than that of the other treatments. The results showed that the osmolality of the whole-body (WB) homogenate was gradually and significantly increased with salinity elevation, except in fish exposed to freshwater and 5 psu treatments, while there were no significant changes in WB Na⁺/K⁺-ATPase activity. Digestive enzymatic activities, i.e., pepsin, α-amylase, alkaline phosphatase, and leucine alanine peptidase (leu-ala) generally increased with salinity, but not aminopeptidase and trypsin. Lysozyme and peroxidase activities increased in fish larvae exposed to 15 and 20 psu. These increases proportionally improved growth performance, with the lowest and the highest final weights observed in fish reared at 0 psu (0.08 ± 0.03 g/larvae) and 20 psu (0.11 ± 0.02 g/larvae), respectively, although the average growth recorded at 20 psu could be biased by the high mortality in this group. Occurrence of skeleton deformities, such as in caudal vertebrae and branchiostegal rays, was significantly higher in fish exposed to the higher osmotic conditions (15.0 ± 1.2% and 10.3 ± 2.1% respectively at 0 psu vs. 31.0 ± 2.9% and 49.0 ± 5.6%, respectively at 15 psu). After the 12.5-h heat shock, survival rates significantly differed between treatments with the highest survival observed in fish submitted to 5 psu (68.9%), followed by those exposed to 0 (27%) and 10 (20%) while all fish died at 15 psu. These findings suggest that the striped catfish larvae could be reared in salinity up to 5 to 10 psu with a higher survival and tolerance to thermal stress when compared to fish maintained in freshwater.

Different transcriptomic architecture of the gill epithelia in Nile and Mozambique tilapia after salinity challenge

Tilapiine fishes of the genus Oreochromis vary in their euryhaline capabilities, therefore inhabiting aquatic environments of different salinities across the African continent. We analyzed the differential gene expression in the gills before and after 6 weeks salinity challenge between the highly tolerant Mozambique tilapia (Oreochromis mossambicus) and the less tolerant Nile tilapia (O. niloticus). The pathways triggered by salinity in both tilapia species reveal immune and cell stress responses as well as turnover of ionocytes. Nevertheless, the actual differential expressed genes vary between these two species, pointing at differential transcriptomic architecture, which likely contribute to the species osmoregulation capabilities in elevated salinities.

The time course of molecular acclimation to seawater in a euryhaline fish

The Arabian pupfish, Aphanius dispar, is a euryhaline fish inhabiting both inland nearly-freshwater desert ponds and highly saline Red Sea coastal lagoons of the Arabian Peninsula. Desert ponds and coastal lagoons, located respectively upstream and at the mouths of dry riverbeds (“wadies”), have been found to potentially become connected during periods of intense rainfall, which could allow the fish to migrate between these different habitats. Flash floods would therefore flush Arabian pupfish out to sea, requiring a rapid acclimation to a greater than 40 ppt change in salinity. To investigate the molecular pathways of salinity acclimation during such events, a Red Sea coastal lagoon and a desert pond population were sampled, with the latter exposed to a rapid increase in water salinity. Changes in branchial gene expression were investigated via genome-wide transcriptome measurements over time from 6 h to 21 days. The two natural populations displayed basal differences in genes related to ion transport, osmoregulation and immune system functions. These mechanisms were also differentially regulated in seawater transferred fish, revealing their crucial role in long-term adaptation. Other processes were only transiently activated shortly after the salinity exposure, including cellular stress response mechanisms, such as molecular chaperone synthesis and apoptosis. Tissue remodelling processes were also identified as transient, but took place later in the timeline, suggesting their importance to long-term acclimation as they likely equip the fish with lasting adaptations to their new environment. The alterations in branchial functional pathways displayed by Arabian pupfish in response to salinity increases are diverse. These reveal a large toolkit of molecular processes important for adaptation to hyperosmolarity that allow for successful colonization to a wide variety of different habitats.

Nonlinear effects of environmental salinity on the gill transcriptome versus proteome of Oreochromis niloticus modulate epithelial cell turnover

A data-independent acquisition (DIA) assay library for targeted quantitation of thousands of Oreochromis niloticus gill proteins using a label- and gel-free workflow was generated and used to compare protein and mRNA abundances. This approach generated complimentary rather than redundant data for 1899 unique genes in gills of tilapia acclimated to freshwater and brackish water. Functional enrichment analyses identified mitochondrial energy metabolism, serine protease and immunity-related functions, and cytoskeleton/ extracellular matrix organization as major processes controlled by salinity in O. niloticus gills. Non-linearity in salinity-dependent transcriptome versus proteome regulation was revealed for specific functional groups of genes. The relationship was more linear for other molecular functions/ cellular processes, suggesting that the salinity-dependent regulation of O. niloticus gill function relies on post-transcriptional mechanisms for some functions/ processes more than others. This integrative systems biology approach can be adopted for other tissues and organisms to study cellular dynamics for many changing ecological contexts.

Transcriptomic responses of saline-adapted Nile tilapia (Oreochromis niloticus) to rearing in both saline and freshwater

The Nile tilapia, Oreochromis niloticus, is an important species for global aquaculture. Recently, a single genetic line of Nile tilapia was developed using estimated breeding values (EBVs) for body weight under moderately saline water that showed significant improvement in growth performance. To explore the molecular mechanisms underlying this enhanced growth capacity, RNA-Seq was used to profile differences in gene expression in the liver and pituitary gland of high- and low-growth performance families of male Nile tilapia progeny, reared in either saline or freshwater environments. Comparisons of tissues from high- and low-EBV families, and also between fish reared in either saline or freshwater, revealed 142 and 2208 differentially expressed genes (DEGs), respectively. DEGs identified between the EBV groups comprised a number of genes involved in the regulation of growth and reproduction. We found an overexpression of hormone genes involved in growth-inhibition in the pituitary of Low-EBV tilapia including 2 somatostatin genes (GHIH), corticoliberin (CRH) and tachykinin-3-like protein. Furthermore, several genes associated with the cAMP pathway were underexpressed in low-EBV tilapia pituitary together with several early response genes. This study provides insight into the transcriptomic factors associated with growth performance in saline-adapted Nile tilapia reared in environments with high and low salinity levels and provides valuable knowledge for the future development of selection strategies to improve growth performance in this species.

Liver and spleen transcriptome reveals that Oreochromis aureus under long-term salinity stress may cause excessive energy consumption and immune response

To investigate the physiological responses of Oreochromis aureus to salinity fluctuations at the molecular level. We used RNA-seq to explore the differentially expressed genes (DEGs) in the liver and spleen of O. aureus at 0, 3, 7 and 11 ppt (parts per thousand) salinity levels. Herein, De novo assembly generated 71,009 O. aureus unigenes, of which 34,607 were successfully mapped to the four major databases. A total of 120 shared DEGs were identified in liver and spleen transcripts, of which 83 were up-regulated and 37 were down-regulated. GO and KEGG analysis found a total of 26 significant pathways, mainly including energy metabolism, immune response, ion transporters and signal transduction. The trend module category of DEGs showed that the genes (e.g., FASN, ODC1, CD22, MRC, TRAV and SLC7 family) involved in the change-stable-change (1) and the constant-change categories (2) were highly sensitive to salinity fluctuations, which were of great value for further study. Based on these results, it would help provide basic data for fish salinity acclimation, and provide new insights into evolutionary response of fish to various aquatic environments in the long-term stress adaptation mechanism.

Growth hormone regulates intestinal gene expression of nutrient transporters in tilapia (Oreochromis mossambicus)

Among the various ways that growth hormone (GH) underlies the growth physiology of teleost fishes, GH stimulates transport pathways that facilitate the absorption of nutrients across intestinal epithelia. The current study investigated the effects of GH on the gene expression of nutrient transporters in an omnivorous teleost, the Mozambique tilapia (Oreochromis mossambicus). We employed pituitary gland removal (hypophysectomy) and hormone replacement to assess whether GH directs the gene expression of the GH receptor (ghr2), the peptide transporters, pept1a, pept1b and pept2, the amino acid transporter, slc7a9, the Na+/glucose cotransporter, sglt1, the glucose transporter, glut2, and the myo-inositol transporter, smit2, in anterior, middle, and posterior intestine. ghr2 was predominantly expressed in posterior intestine, while pept1a, pept1b, slc7a9, sglt1, glut2, and smit2 exhibited the highest mRNA levels in anterior and/or middle intestine. While hypophysectomized tilapia exhibited diminished expression of ghr2, pept1a, pept1b, slc7a9, and glut2 compared with intact and sham-operated controls, only ghr2, pept1a, pept1b and glut2 levels were restored by GH replacement. Our findings indicate that GH supports growth, at least in part, by stimulating the gene expression of its cognate receptor and key nutrient transporters in the intestine.

The round goby genome provides insights into mechanisms that may facilitate biological invasions

Background

The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success.

Results

We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby’s capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions.

Conclusions

The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0731-8) contains supplementary material, which is available to authorized users.

Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation

Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.

Transcriptome profiling of the low-salinity stress responses in the gills of the juvenile Pseudopleuronectes yokohamae

Salinity is an important environmental factor that affects the life cycle of fish, including their growth, development and reproduction. The marbled flounder, Pseudopleuronectes yokohamae, is an important economic resource and serves as a good model to investigate osmoregulation, as it can adapt to a wide range of salinity levels. However, the lack of genomic resources for this species has hampered the understanding of the mechanisms underlying its salinity tolerance. In this study, RNA-Seq analysis was conducted to identify genes related to salinity adaptation and osmotic regulation in the gill tissue of marbled flounder exposed to different concentrations of environmental salinity (6 and 30 ppt). After de novo assembly, 19,265 genes were annotated by the Nr database. A comparison of expression between the two salinity groups revealed 673 differentially expressed genes, of which 180 were upregulated and 493 were downregulated in the low salinity group relative to the high salinity group. The related molecular biological processes were explored from several important perspectives, and potential functions were determined by enrichment analyses, including those of metabolites in ion transportation, energy metabolism and protein synthesis, and immune responses. This study is the first transcriptomic study conducted on marbled flounder, and it revealed many novel sequences for further biological analyses. In addition, the candidate genes identified in the gene expression analysis provided insights into responses to salinity change and molecular mechanisms underlying osmoregulation in the gills of marbled flounder.

Peptide Transporters in the Primary Gastrointestinal Tract of Pre-Feeding Mozambique Tilapia Larva

Fish larvae differ greatly from the adult form in their morphology and organ functionality. The functionality of the gastrointestinal tract depends on the expression of various pumps, transporters, and channels responsible for feed digestion and nutrients absorption. During the larval period, the gastrointestinal tract develops from a simple closed tube, into its complex form with differentiated segments, crypts and villi, as found in the adult. In this study, we characterized the expression of three peptide transporters (PepT1a, PepT1b, and PepT2) in the gastrointestinal tract of Mozambique tilapia (Oreochromis mossambicus) larvae along 12 days of development, from pre-hatching to the completion of yolk sac absorption. Gene expression analysis revealed differential and complimentary time-dependent expression of the PepT1 variants and PepT2 along the larval development period. Immunofluorescence analysis showed differential protein localization of the three peptide transporters (PepTs) along the gastrointestinal tract, in a similar pattern to the adult. In addition, PepT1a was localized in mucosal cells in the larvae esophagus, in much higher abundance than in the adults. The results of this study demonstrate specialization of intestinal sections and absorbance potential of the enterocytes prior to the onset of active exogenous feeding, thus pointing to an uncharacterized function and role of the gastrointestinal tract and its transporters during the larval period.

Transcriptome Analysis Reveals Common and Differential Response to Low Temperature Exposure Between Tolerant and Sensitive Blue Tilapia (Oreochromis aureus)

Tilapias are very important to the world's aquaculture. As befitting fish of their tropical origin, their distribution, and culture practices are highly affected by low temperatures. In this study, we used genetic and genomic methodologies to reveal pathways involved in the response and tolerance of blue tilapia (Oreochromis aureus) to low temperature stress. Cold tolerance was characterized in 66 families of blue tilapia. Fish from cold-tolerant and cold-sensitive families were sampled at 24 and 12°C, and the transcriptional responses to low-temperature exposure were measured in the gills and liver by high-throughput mRNA sequencing. Four hundred and ninety four genes displayed a similar temperature-dependent expression in both tolerant and sensitive fish and in the two tissues, representing the core molecular response to low temperature exposure. KEGG pathway analysis of these genes revealed down-regulation of focal-adhesion and other cell-extracellular matrix (ECM) interactions, and up-regulation of proteasome and various intra-cellular proteolytic activities. Differential responses between cold-tolerant and cold-sensitive fish were found with genes and pathways that were up-regulated in one group and down-regulated in the other. This reverse response was characterized by genes involved in metabolic pathways such as glycolysis/gluconeogenesis in the gills and biosynthesis of amino-acids in the liver, with low temperature down-regulation in tolerant fish and up-regulation in sensitive fish.

Liver transcriptome analysis reveals extensive transcriptional plasticity during acclimation to low salinity in Cynoglossus semilaevis

BACKGROUND

Salinity is an important abiotic stress that influences the physiological and metabolic activity, reproduction, growth and development of marine fish. It has been suggested that half-smooth tongue sole (Cynoglossus semilaevis), a euryhaline fish species, uses a large amount of energy to maintain osmotic pressure balance when exposed to fluctuations in salinity. To delineate the molecular response of C. semilaevis to different levels of salinity, we performed RNA-seq analysis of the liver to identify the genes and molecular and biological processes involved in responding to salinity changes.

RESULTS

The present study yielded 330.4 million clean reads, of which 83.9% were successfully mapped to the reference genome of C. semilaevis. One hundred twenty-eight differentially expressed genes (DEGs), including 43 up-regulated genes and 85 down-regulated genes, were identified. These DEGs were highly represented in metabolic pathways, steroid biosynthesis, terpenoid backbone biosynthesis, butanoate metabolism, glycerolipid metabolism and the 2-oxocarboxylic acid metabolism pathway. In addition, genes involved in metabolism, osmoregulation and ion transport, signal transduction, immune response and stress response, and cytoskeleton remodeling were affected during acclimation to low salinity. Genes acat2, fdps, hmgcr, hmgcs1, mvk, pmvk, ebp, lss, dhcr7, and dhcr24 were up-regulated and abat, ddc, acy1 were down-regulated in metabolic pathways. Genes aqp10 and slc6a6 were down-regulated in osmoregulation and ion transport. Genes abat, fdps, hmgcs1, mvk, pmvk and dhcr7 were first reported to be associated with salinity adaptation in teleosts.

CONCLUSIONS

Our results revealed that metabolic pathways, especially lipid metabolism were important for salinity adaptation. The candidate genes identified from this study provide a basis for further studies to investigate the molecular mechanism of salinity adaptation and transcriptional plasticity in marine fish.

Osmolality/salinity-responsive enhancers (OSREs) control induction of osmoprotective genes in euryhaline fish

Fish respond to salinity stress by transcriptional induction of many genes, but the mechanism of their osmotic regulation is unknown. We developed a reporter assay using cells derived from the brain of the tilapia Oreochromis mossambicus (OmB cells) to identify osmolality/salinity-responsive enhancers (OSREs) in the genes of Omossambicus Genomic DNA comprising the regulatory regions of two strongly salinity-induced genes, inositol monophosphatase 1 (IMPA1.1) and myo-inositol phosphate synthase (MIPS), was isolated and analyzed with dual luciferase enhancer trap reporter assays. We identified five sequences (two in IMPA1.1 and three in MIPS) that share a common consensus element (DDKGGAAWWDWWYDNRB), which we named "OSRE1." Additional OSREs that were less effective in conferring salinity-induced trans-activation and do not match the OSRE1 consensus also were identified in both MIPS and IMPA1.1 Although OSRE1 shares homology with the mammalian osmotic-response element/tonicity-responsive enhancer (ORE/TonE) enhancer, the latter is insufficient to confer osmotic induction in fish. Like other enhancers, OSRE1 trans-activates genes independent of orientation. We conclude that OSRE1 is a cis-regulatory element (CRE) that enhances the hyperosmotic induction of osmoregulated genes in fish. Our study also shows that tailored reporter assays developed for OmB cells facilitate the identification of CREs in fish genomes. Knowledge of the OSRE1 motif allows affinity-purification of the corresponding transcription factor and computational approaches for enhancer screening of fish genomes. Moreover, our study enables targeted inactivation of OSRE1 enhancers, a method superior to gene knockout for functional characterization because it confines impairment of gene function to a specific context (salinity stress) and eliminates pitfalls of constitutive gene knockouts (embryonic lethality, developmental compensation).

RNA-Seq analysis of salinity stress-responsive transcriptome in the liver of spotted sea bass (Lateolabrax maculatus)

Salinity is one of the most prominent abiotic factors, which greatly influence reproduction, development, growth, physiological and metabolic activities of fishes. Spotted sea bass (Lateolabrax maculatus), as a euryhaline marine teleost, has extraordinary ability to deal with a wide range of salinity changes. However, this species is devoid of genomic resources, and no study has been conducted at the transcriptomic level to determine genes responsible for salinity regulation, which impedes the understanding of the fundamental mechanism conferring tolerance to salinity fluctuations. Liver, as the major metabolic organ, is the key source supplying energy for iono- and osmoregulation in fish, however, little attention has been paid to its salinity-related functions but which should not be ignored. In this study, we perform RNA-Seq analysis to identify genes involved in salinity adaptation and osmoregulation in liver of spotted sea bass, generating from the fishes exposed to low and high salinity water (5 vs 30ppt). After de novo assembly, annotation and differential gene expression analysis, a total of 455 genes were differentially expressed, including 184 up-regulated and 271 down-regulated transcripts in low salinity-acclimated fish group compared with that in high salinity-acclimated group. A number of genes with a potential role in salinity adaptation for spotted sea bass were classified into five functional categories based on the gene ontology (GO) and enrichment analysis, which include genes involved in metabolites and ion transporters, energy metabolism, signal transduction, immune response and structure reorganization. The candidate genes identified in L. maculates liver provide valuable information to explore new pathways related to fish salinity and osmotic regulation. Besides, the transcriptomic sequencing data supplies significant resources for identification of novel genes and further studying biological questions in spotted sea bass.

Salinity-Dependent Shift in the Localization of Three Peptide Transporters along the Intestine of the Mozambique Tilapia (Oreochromis mossambicus)

The peptide transporter (PepT) systems are well-known for their importance to protein absorption in all vertebrate species. These symporters use H⁺ gradient at the apical membrane of the intestinal epithelial cells to mediate the absorption of small peptides. In fish, the intestine is a multifunctional organ, involved in osmoregulation, acid-base regulation, and nutrient absorption. Therefore, we expected environmental stimuli to affect peptide absorption. We examined the effect of three environmental factors; salinity, pH and feeding, on the expression, activity and localization of three PepT transporters (PepT1a, PepT1b, PepT2) along the intestine of the Mozambique tilapia (Oreochromis mossambicus). Quantitative real time PCR (qPCR) analysis demonstrated that the two PepT1 variants are typical to the proximal intestinal section while PepT2 is typical to the distal intestinal sections. Immunofluorescence analysis with custom-made antibodies supported the qPCR results, localized both transporters on the apical membrane of enterocytes and provided the first evidence for the participation of PepT2 in nutrient absorption. This first description of segment-specific expression and localization points to a complementary role of the different peptide transporters, corresponding to the changes in nutrient availability along the intestine. Both gene expression and absorption activity assays showed that an increase in water salinity shifted the localization of the PepT genes transcription and activity down along the intestinal tract. Additionally, an unexpected pH effect was found on the absorption of small peptides, with increased activity at higher pH levels. This work emphasizes the relationships between different functions of the fish intestine and how they are affected by environmental conditions.

Di- and tripeptide transport in vertebrates: the contribution of teleost fish models

Solute Carrier 15 (SLC15) family, alias H⁺-coupled oligopeptide cotransporter family, is a group of membrane transporters known for their role in the cellular uptake of di- and tripeptides (di/tripeptides) and peptide-like molecules. Of its members, SLC15A1 (PEPT1) chiefly mediates intestinal absorption of luminal di/tripeptides from dietary protein digestion, while SLC15A2 (PEPT2) mainly allows renal tubular reabsorption of di/tripeptides from ultrafiltration, SLC15A3 (PHT2) and SLC15A4 (PHT1) possibly interact with di/tripeptides and histidine in certain immune cells, and SLC15A5 has unknown function. Our understanding of this family in vertebrates has steadily increased, also due to the surge of genomic-to-functional information from 'non-conventional' animal models, livestock, poultry, and aquaculture fish species. Here, we review the literature on the SLC15 transporters in teleost fish with emphasis on SLC15A1 (PEPT1), one of the solute carriers better studied amongst teleost fish because of its relevance in animal nutrition. We report on the operativity of the transporter, the molecular diversity, and multiplicity of structural-functional solutions of the teleost fish orthologs with respect to higher vertebrates, its relevance at the intersection of the alimentary and osmoregulative functions of the gut, its response under various physiological states and dietary solicitations, and its possible involvement in examples of total body plasticity, such as growth and compensatory growth. By a comparative approach, we also review the few studies in teleost fish on SLC15A2 (PEPT2), SLC15A4 (PHT1), and SLC15A3 (PHT2). By representing the contribution of teleost fish to the knowledge of the physiology of di/tripeptide transport and transporters, we aim to fill the gap between higher and lower vertebrates.

Gene Expression Profiling Reveals Functional Specialization along the Intestinal Tract of a Carnivorous Teleostean Fish (Dicentrarchus labrax)

High-quality sequencing reads from the intestine of European sea bass were assembled, annotated by similarity against protein reference databases and combined with nucleotide sequences from public and private databases. After redundancy filtering, 24,906 non-redundant annotated sequences encoding 15,367 different gene descriptions were obtained. These annotated sequences were used to design a custom, high-density oligo-microarray (8 × 15 K) for the transcriptomic profiling of anterior (AI), middle (MI), and posterior (PI) intestinal segments. Similar molecular signatures were found for AI and MI segments, which were combined in a single group (AI-MI) whereas the PI outstood separately, with more than 1900 differentially expressed genes with a fold-change cutoff of 2. Functional analysis revealed that molecular and cellular functions related to feed digestion and nutrient absorption and transport were over-represented in AI-MI segments. By contrast, the initiation and establishment of immune defense mechanisms became especially relevant in PI, although the microarray expression profiling validated by qPCR indicated that these functional changes are gradual from anterior to posterior intestinal segments. This functional divergence occurred in association with spatial transcriptional changes in nutrient transporters and the mucosal chemosensing system via G protein-coupled receptors. These findings contribute to identify key indicators of gut functions and to compare different fish feeding strategies and immune defense mechanisms acquired along the evolution of teleosts.

Developmental transcriptome analysis and identification of genes involved in formation of intestinal air-breathing function of Dojo loach, Misgurnus anguillicaudatus

Dojo loach, Misgurnus anguillicaudatus is a freshwater fish species of the loach family Cobitidae, using its posterior intestine as an accessory air-breathing organ. Little is known about the molecular regulatory mechanisms in the formation of intestinal air-breathing function of M. anguillicaudatus. Here high-throughput sequencing of mRNAs was performed from six developmental stages of posterior intestine of M. anguillicaudatus: 4-Dph (days post hatch) group, 8-Dph group, 12-Dph group, 20-Dph group, 40-Dph group and Oyd (one-year-old) group. These six libraries were assembled into 81300 unigenes. Totally 40757 unigenes were annotated. Subsequently, 35291 differentially expressed genes (DEGs) were scanned among different developmental stages and clustered into 20 gene expression profiles. Finally, 15 key pathways and 25 key genes were mined, providing potential targets for candidate gene selection involved in formation of intestinal air-breathing function in M. anguillicaudatus. This is the first report of developmental transcriptome of posterior intestine in M. anguillicaudatus, offering a substantial contribution to the sequence resources for this species and providing a deep insight into the formation mechanism of its intestinal air-breathing function. This report demonstrates that M. anguillicaudatus is a good model for studies to identify and characterize the molecular basis of accessory air-breathing organ development in fish.

Flexible selection of diversified Na(+)/K(+)-ATPase α-subunit isoforms for osmoregulation in teleosts

BACKGROUND AND METHODS

Multiple Na+/K+-ATPase (NKA) α-subunit isoforms express differentially in response to salinity transfer in teleosts but we observed that the isoform nomenclature is inconsistent with the phylogenetic relationship of NKA α-genes. We cloned the catalytic NKA α-subunit isoforms in eels and medaka, analyzed the time course of their expressions in osmoregulatory tissues after transfer from freshwater (FW) to seawater (SW), and performed phylogenetic analyses to deduce an evolutionary scenario that illustrates how various duplication events have led to the current genomic arrangement of NKA α-genes in teleosts.

RESULTS AND DISCUSSION

Five and six α-subunits were cloned in eels and medaka respectively. In eels, the commonly-reported α1a and α1b isoforms were absent while the α1c isoform was diversified instead (α1c-1, α1c-2, α1c-3, α2, and α3 in eels). Phylogenetic estimation indicated that the α1a and α1b isoforms from salmon, tilapia, and medaka were generated by independent duplication events and thus they are paralogous isoforms. Re-examination of expression changes of known isoforms after salinity challenge revealed that the isoforms selected as predominant SW-types varied among teleost lineages. Diversification of α1 isoforms occurred by various types of gene duplication, or by alternative transcription among tandem genes to form chimeric transcripts, but there is no trend for more α1 copies in euryhaline species. Our data suggest that the isoform switching between FW (α1a predominates) and SW (α1b predominates) that occurs in salmonids is not universal in teleosts. Instead, in eels, α1c-1 was the major α-subunit upregulated gill, intestine, and kidney in SW. Localization of both NKA mRNA and protein showed consistent upregulation in gill and intestine in SW eels, but not in renal distal and collecting tubules, where low transcript expression levels were accompanied by high protein levels, suggesting a tissue-specific translational regulation that determines and fine-tunes the NKA expression. In medaka, α1b was upregulated in SW in anterior intestine while most other α-subunit isoforms were less responsive to salinity changes.

CONCLUSION

By integrating gene expression and phylogenetic results, we propose that the major NKA α-subunits for SW acclimation were not ancestrally selected, but rather were flexibly determined in lineage-specific fashion in teleosts.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.