Josep A. Calduch-Giner1
Jaume Pérez-Sánchez- 1Nutrigenomics and Fish Growth Endocrinology Group, Biology, Culture and Pathology of Marine Species, Institute of Aquaculture Torre de la Sal (IATS-CSIC), Castellón, Spain
- 2Fish Pathology Group, Biology, Culture and Pathology of Marine Species, Institute of Aquaculture Torre de la Sal (IATS-CSIC), Castellón, Spain
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.
Introduction
Teleost fish constitute the most abundant vertebrate group, exhibiting an amazing level of biodiversity affecting their morphology, ecology, and behavior as well as many other aspects of their biology. This huge variability makes fish extremely attractive for the study of many biological questions. They show a high assortment of morphology, anatomy and histology of their gastrointestinal (GI) tract in association with their numerous specialized functions (Wilson and Castro, 2010), as it is involved not only in digestion and feed absorption but also in water and electrolyte balance, nutrient sensing, and immunity (Cain and Swan, 2010). This diversity is now starting to be elucidated and molecular approaches are helping to understand the many vital functions conducted by the GI tract in the context of integrative and comparative physiology. Thus, fish GI tract is an important immunological site acting as a physical and chemical barrier against invading organisms, and the cells within the intestine produce a range of chemical substances to enhance barrier function, innate immunity, and humoral immune responses (Rombout et al., 2011; Xia et al., 2013). This is particularly important at the posterior intestine, where the uptake of macromolecules and foreign antigens mainly occurs (Fuglem et al., 2010; Løkka et al., 2014). Nutrient availability is also essential for the generation and maintenance of a protective effector immune system in humans and other model species (Fox et al., 2005; Michalek and Rathmell, 2010). For instance, a deficiency in glucose uptake negatively affects T cell function with impairment of both proliferation and cytokine expression. Similarly, deficiencies in amino acids, such as tryptophan, arginine, glutamine and cysteine, reduce immune cell activation. Furthermore, short-chain fatty acids provide one of the clearest examples of how host diet and nutrient processing by the gut microbiota shape immune responses (Kau et al., 2011). This is extensive to fish, as dietary butyrate was able to revert most of the changes induced in the intestinal transcriptome of gilthead sea bream (Sparus aurata) by the replacement of fish meal and fish oil with plant ingredients (Estensoro et al., 2014). However, experimental evidence in Atlantic halibut (Hippoglossus hippoglossus) (Murray et al., 2010), Atlantic salmon (Salmo salar) (Skugor et al., 2011), and Atlantic cod (Gadus morhua) (Morais et al., 2012a; De Santis et al., 2015) suggests that the nutritionally mediated changes in the intestinal transcriptome and function are highly fish species-specific.
Changes in the intestinal transcriptome also reflect the effect of nutritional background on the disease outcome of gilthead sea bream infected with the intestinal parasite Enteromyxum leei (Estensoro et al., 2011; Calduch-Giner et al., 2012). Furthermore, the high expression level of an intestinal mucin specific to the fish lineage has been proven a suitable marker of diagnostic and prognostic value of this parasitic enteritis (Pérez-Sánchez et al., 2013). In this context, the functional outline of the intestinal transcriptome of Mediterranean, Perciform, carnivorous, and gastric fish, such as gilthead sea bream and European sea bass (Dicentrarchus labrax), will allow comparative studies with other earlier teleostean groups, such as Salmoniformes (salmon, trout) and Cypriniformes (zebrafish, carp), which have different feeding behaviors or are stomachless (agastric fish). It would also help to understand the impact of dietary restriction, intermittent feeding, or compensatory growth (Inness and Metcalfe, 2008). In addition, this transcriptomic profiling can be a valuable tool to define the nutritional value of novel fish feeds and their effect on fish immune system.
In European sea bass, massive gene expression analysis have undertaken for gills, liver, brain, eggs, and larvae (Geay et al., 2011; Magnanou et al., 2014; Nuñez Ortiz et al., 2014; Kaitetzidou et al., 2015). A high-quality chromosome-scale genome assembly is also available (Tine et al., 2014) for this species, but in-depth analysis of the intestinal transcriptome remains to be fully addressed. Thus, the main goal of this work was to compare three different intestinal segments of European sea bass in order to generate a precise functional map of the intestine of this valuable fish species. For such purpose, transcriptomic data were generated through the construction of next-generation sequencing libraries. These sequences, combined with others from public and private repositories, were used to construct a reliable assembly nucleotide database (http://www.nutrigroup-iats.org/seabassdb) enriched on intestinal sequences, being used the annotated assembled sequences to design a custom oligo-microarray.
Materials and Methods
Animals and Tissue Sampling
European sea bass of Mediterranean origin were reared under standard conditions at densities lower than 15 kg/m3 at indoor experimental facilities of the Institute of Aquaculture Torre de la Sal (IATS-CSIC). Photoperiod and temperature followed the natural changes at our latitude (40°5′N, 0°10′E), and fish were fed to visual satiety on 5 days per week with a commercial diet (Biomar, Palencia, Spain). At midsummer (fast growing fish), overnight fasted fish of 200 g average weight were anesthetized with 100 mg/L of 3-aminobenzoic acid ethyl ester (MS-222, Sigma-Aldrich, Madrid, Spain), bled, and samples of anterior (AI), middle (MI), and posterior (PI) intestinal segments were rapidly excised, frozen in liquid nitrogen and stored at −80°C until RNA extraction. AI was taken as a segment immediately after the pyloric caeca. MI was taken as a representative portion at the middle of the intestine. PI was taken as a segment including the rectum. All procedures were approved by the Ethics and Animal Welfare Committee of the Institute of Aquaculture Torre de la Sal according to national (Royal Decree RD53/2013) and EU legislation (2010/63/EU) on the handling of experimental animals.
RNA Extraction
Total RNA was isolated by means of the Ambion MagMax-96 for Microarray kit (Applied Biosystems, Foster City, CA, USA) after tissue homogenization in TRI reagent at a concentration of 100 mg/ml following the manufacturers' instructions. RNA quantity and purity was determined by Nanodrop (Thermo Scientific, Wilgminton, DE, USA) and Agilent 2100 bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) measurements with absorbance ratios at 260 nm/280 nm above 1.9 and RNA integrity numbers between 9.2 and 10, which are indicative of clean and intact RNA.
Transcriptome Libraries Sequencing and Asssembly
Tissue samples from four fish were used for the construction of next-generation sequencing libraries of AI and PI segments. cDNA synthesis was performed with 700 ng total RNA by means of the MINT kit (Evrogen, Heidelberg, Germany). To increase the recovery rate of rare and unique transcripts, amplified cDNAs were normalized by duplex-specific nuclease with the Trimmer kit (Evrogen) following the manufacturer's instructions (Zhulidov et al., 2004). Normalized cDNA samples were quantified with the Quant-iT PicoGreen dsDNA quantification Kit (Life Technologies, Carlsbad, CA, USA) using a VersaFluorTM Fluorometer system (Bio-Rad, Hercules, CA, USA). Normalized cDNA (500 ng) was sheared into small fragments (250–600 nt) by nebulization with compressed nitrogen. Then, sequencing adapters were ligated to the blunt ends of the fragments, and an emulsion PCR (emPCR) was performed. After emPCR, the enriched beads were loaded onto the 454 microtiter plate and the amplicons were sequenced with a Titanium GS FLX 454 platform (Roche, Mannheim, Germany). The sff files containing all reads for each library have been deposited to NCBI Short Read Archive under accession SRR1560808.
The quality of the reads was assessed with PERL scripts developed at Lifesequencing S.L. (Valencia, Spain) for trimming of adaptors and validation of high-quality sequences. All sequences were edited to remove vector and adaptor sequences, and cleaned and filtered before assembly. Cleaning involved masking of poor-quality bases and low-complexity sequences, such as polyA tails. Filtering removed contaminating sequences (bacteria, yeast) and only high-quality sequences (q > 25) of more than 100 bases in length were retained. Assembly of high-quality reads was performed by means of Newbler version 2.8 program (454 Life Science-Roche) using the default parameters and the cDNA option.
European Sea Bass Transcriptome Database Construction
According to the scheme pipeline shown in Figure 1, assembled contigs and singletons from the de novo assembly (85,350 sequences) were combined with European sea bass annotated nucleotide sequences extracted from public (European sea bass complete mRNAs from Genbank, 492 sequences) and private (AQUAFIRST EU project at http://www.sigenae.org/aquafirst, 6953 sequences) databases. Redundant sequences were removed retaining the longest ones, and filtered sequences were annotated by similarity searches using protein reference databases, such as UniProt/Swissprot, RefSeq protein, TrEMBL, and Pfam. The e-value threshold to determine similarities was set to 1E-5, and the Uniprot/Swissprot entry to which they received the highest similarity was usually assigned as the gene identity. After annotation, the frameshift edition algorithm previously developed for sequencing corrections at homopolymer regions in gilthead sea bream transcriptome assemblies was applied (Calduch-Giner et al., 2013). Functional characterization by means of gene ontology analysis was made from the most representative contigs/singletons for each gene identity using the Blast2GO software (Conesa et al., 2005) with a threshold cutoff set at 1E-3. Multilevel GO term analysis was performed under default settings (minimum cutoff 1E-04 BlastX) and results were filtered by annotation score.
Figure 1. Schematic representation of the data processing pipeline for the de novo assembly and annotation of European sea bass intestine transcriptome. Distribution of biological process multilevel GO annotation terms in European sea bass sequences. The number of sequences for each GO term is represented between parentheses.
Microarray Construction, Hybridization, and Data Analysis
A custom high-density oligo-microarray (8 × 15 K) from the assembled nucleotide European sea bass sequences was designed and printed using the eArray web tool (Agilent). The array comprised 60-oligomer probes for 14,147 different European sea bass annotated sequences. Total RNA (150 ng) from individuals (n = 8 for each intestine segment: AI, MI, and PI) were labeled with cyanine 3-CTP (Low Input Quick Amp Labeling Kit, Agilent), and 600 ng of each labeled cRNA were hybridized to microarray slides that were analyzed with an Agilent G2565C Microarray Scanner according to the manufacturer's protocol. Data were extracted using the Agilent Feature Extraction Software 11.5.1.1. Data analysis of differentially expressed (DE) genes was carried out with the Genespring GX 13.0 software (Agilent). Functional pathway analysis of DE genes was performed with the Ingenuity Pathway Analysis (IPA) software (http://www.ingenuity.com). For each gene, the Uniprot accession of the annotation equivalent for one of the three higher vertebrates model species in IPA (human, rat, or mouse) was assigned.
Real-Time qPCR Validation
Up to 29 genes that were considered representative of DE genes were validated on individual AI and PI samples (n = 6 for each intestinal segment) by real-time qPCR, using an iCycler IQ Real-time Detection System (Bio-Rad). cDNA synthesis was performed using the High-Capacity cDNA Archive Kit (Applied Biosystems) with random decamers. For this purpose, 500 ng total RNA were reverse transcribed in a final volume of 100 μl. RT reactions were incubated for 10 min at 25°C and 2 h at 37°C. Negative control reactions were incubated in the absence of reverse transcriptase. Diluted RT reactions were conveniently used for PCR reactions in a 25 μl volume. Each PCR well contained a SYBR Green Master Mix (Bio-Rad) and specific primers (Table S1) were used at a final concentration of 0.9 μM. DNA Polymerase was activated and cDNA denatured by preincubation for 3 min at 95°C; the template was amplified for 40 cycles of denaturation at 95°C for 15 s, and annealing/extension at 60°C for 60 s. β-Actin was used as the housekeeping gene and the efficiency of PCR reactions for target and reference genes varied between 90 and 98%. The amount of product in a particular sample was determined by interpolation of the cycle threshold (Ct) value. The specificity of the reaction was verified by analysis of melting curves and by electrophoresis and sequencing of PCR-amplified products. Fluorescence data acquired during the extension phase were ultimately normalized to β-actin by the ΔΔCt method (Livak and Schmittgen, 2001). For each selected gene, fold-change variations were calculated between posterior and anterior regions.
Statistical Analysis
Microarray results for the different intestinal segments were compared in pairs by means of a t-test (corrected P < 0.05, Benjamini-Hochberg multiple testing correction).
Results
Transcriptome Assembly and Annotation
As stated in Table 1, sequencing of AI and PI segments yielded similar results in terms of the number of high-quality reads and average read length. The combined de novo assembly of all high-quality reads resulted in a total of 85,350 unique sequences (20,554 contigs and 64,796 singletons) with an average length of 788 nucleotides. After combination with annotated assembled reads from Aquafirst and GenBank complete mRNAs, filtering resulted in 24,906 non-redundant annotated sequences with an average length of 1050 nucleotides (N50 = 1323 nucleotides) that were uploaded to our nucleotide database (http://www.nutrigroup-iats.org/seabassdb). The assembled transcriptome yielded a total of 15,367 different gene descriptions for annotated sequences. Blast2GO analysis showed a high abundance among others of GO terms related to organ development (GO:0048513, 1921 different genes), cell communication, cell surface receptor signaling and signal transduction (GOs:0010646, 0007166, 0009966; 1683, 1663, and 1534 genes, respectively), cellular assembly (GO:0022607, 1464 genes), protein transport (GO:0015031, 1349 genes), immune system (GO:0002376, 1285 genes), cell development (GO:0048468, 1215 genes), and lipid metabolism (GO:0006629, 1195 genes) (Figure 1).
Table 1. Statistics for 454 pyrosequencing libraries.
Most of the non-redundant assembled sequences had a reliable annotation against the UniProt/Swissprot database (17,181 sequences, 69%). Sequences with no significant equivalent in this database were annotated by similarity to RefSeq-Prot (4004 sequences), Trembl (3160 sequences), or Pfam (561 sequences) (Figure 2A). In all cases, reliable annotations were obtained with e-values lower than 1E-30 for more than 60% of sequences (Figure 2B).
Figure 2. Annotation statistics for the European sea bass transcriptome. (A) Top annotation hits by reference databases. (B) Annotation best hit e-value distribution.
Functional Analysis of Differentially Expressed (DE) Genes
Microarray gene expression profiling of individual fish showed a gradient expression pattern along the intestinal tract with a close association of AI and MI by principal component analysis, with the first two components accounting for 41% of the total variance (Figure 3). Hence, t-tests (Benjamini-Hochberg, P < 0.05) were not able to detect DE genes between AI and MI samples, and these two intestinal segments were combined into one category (AI-MI) for further analysis comparing AI-MI vs. PI. This strategy yielded 5770 DE genes (Benjamini-Hochberg t-test, P < 0.05), reduced to 1906 (960 upregulated in AI-MI; 946 upregulated in PI) after filtering with a fold-change cutoff of 2.0. The complete list of DE genes with data on fold-changes is shown as Table S2.
Figure 3. Principal component analysis of the gene expression profile of intestine sections of European sea bass. Anterior intestine samples are marked as AI (red), middle intestine as MI (blue) and posterior intestine as PI (green).
Functional pathway analysis of DE genes revealed a different pattern of molecular and cellular functions throughout the intestine. In AI-MI, there was an overrepresentation of genes related to digestion and transport of nutrients (lipids, vitamins, minerals, carbohydrates, amino acids) (Figure 4). This was consistent with the observation that highly expressed genes in the AI-MI segments in comparison with PI included a vast array of proteins (pepsin A, chymotrypsin-like elastase family member 2A, elastase-1), lipid (phosphodiesterase family member 7, neutral ceramidase, bile salt-activated lipase), and carbohydrate (lactase-phlorizin hydrolase) digestion-related genes, as well as genes involved in lipoprotein metabolism (14 kDa apolipoprotein, apolipoprotein C-II) and intestinal motility (neuromedin-B-like) (Table 2). In PI, the most relevant molecular and cellular functions (17 out of 20) were coincident with those of AI-MI, but they were differently ranked, and the most significant ones were related to the preservation of tissue integrity and cell-to-cell communication (cellular movement, cellular function and maintenance, cell-to-cell signaling and interaction, cell morphology, cellular development, cellular growth and proliferation, cell death and survival) (Figure 5). At the same time, the list of the highest expressed genes in PI in comparison with AI-MI (Table 3) included markers of cell positional identity and epithelial barrier function (sorting nexin-6-like, homeobox protein Hox-D4b, sorting nexin 10). Other genes with higher expression in the PI were those related to absorption of vitamin B12 (protein amnionless, cubilin, gastric intrinsic factor-like, transcobalamin-2) and bile acids (ileal sodium/bile acid cotransporter, gastrotropin). Genes related to immune or antimicrobial response were also represented in the top list of genes with a specific or upregulated expression in the PI (immunoglobulin-like and fibronectin type III domain-containing protein 1, transmembrane and immunoglobulin domain-containing protein 1, histidine ammonia-lyase, B-L beta chain class II histocompatibility antigen, and beta defensin) (Table 3). Other immune-related DE genes in the PI, not included in the top list, were lysozymes, mieloperoxidases, and some TLRs and cytokines (Table S2). This was supported by the effector network analysis conducted by means of IPA software, which identified the immune-related effector network (24 genes connecting chemotaxis, inflammatory responses and leukocyte activation and proliferation) as one of the most distinctive gene signatures of the distal intestine (Figure 6). Likewise, cell signaling comprising a vast array of G-protein coupled receptors (GPCRs) was identified as a highly significant network (Figure 7). Interestingly, many of the receptors that were mainly expressed in AI-MI segments were related to intestinal secretion and motility (GPR39; GPR112; alpha-2C adrenergic receptor, ADRA2C; neuropeptide B/W receptor types 1 and 2, NPBWR1, NPBWR2), whereas those GPCRs with intestinal cell proliferation and inflammatory function (GPR18; GPR63; GPR84; hydroxycarboxylic acid receptor 2-like, HCAR2; frizzled class receptor 10, FZD10) usually had a higher expression level in the PI segment.
Figure 4. Molecular and cellular functions of differentially expressed genes in the anterior-middle intestine segment (AI-MI). For a given molecular function, blue bars indicate the −log of the P-value, and the number of differentially expressed genes is represented.
Table 2. Filtering by fold-change of differentially expressed genes through the intestine segments.
Figure 5. Molecular and cellular functions of differentially expressed genes in the posterior intestine segment (PI). For a given molecular function, blue bars indicate the −log of the P-value, and the number of differentially expressed genes is represented.
Table 3. Filtering by fold-change of differentially expressed genes through the intestine segment.
Figure 6. Regulator effector network of genes differentially expressed in the posterior intestine of European sea bass. It integrates the upstream regulators (upper line), differentially expressed genes (middle line), and downstream effects results (lower line).
Figure 7. Gene network generated by IPA analysis. Red boxes indicate GPCRs differentially and highly expressed at posterior intestine (fold-change, expression ratio between posterior intestine and anterior-middle intestine, >1). Green boxes indicate GPCRs differentially and highly expressed at anterior-middle intestine (fold-change < 1).
Real-Time qPCR Validation
Selected genes for qPCR validation of microarray results were representative of the main processes along the intestinal tract: (i) amino acid transport and proteolysis [B(0,+)-type amino acid transporter 1, excitatory amino acid transporter 3, cathepsin L], (ii) lipid metabolism and digestion (apolipoprotein C-II, bile salt-activated lipase, phosphoethanolamine-N-methyltransferase 3), (iii) integrity of the epithelial barrier (sorting nexin 6), (iv) antimicrobial action (acidic mammalian chitinase, beta-defensin, liver-expressed antimicrobial peptide 2, unconventional myosin-Vb), (v) inflammatory and immune responses (histidine ammonia-lyase, immunoglobulin-like and fibronectin type III domain-containing protein 1, lysozyme C, neuromedin-B, platelet glycoprotein 4, transmembrane, and immunoglobulin domain-containing protein 1), (vi) intestinal chemosensation (GPR18, GPR39, GPR63, GPR84, GPR112, NPBWR1), and (vii) metabolism of bile acids (gastrotropin, ileal sodium/bile acid cotransporter) and vitamin B12 (cubilin, protein amnionless, transcobalamin-2). A scatter plot of microarray and qPCR fold-change values for all the above-mentioned genes highlighted a strong linear correlation (r = 0.97) near to equality (Figure 8).
Figure 8. qPCR validation of microarray results. Correlation plot of fold-change (posterior/anterior) values for the selected genes analyzed by qPCR (X-axis) and microarray (Y-axis).
Discussion
In the present work, the assembly and annotation of 454 normalized libraries resulted in a transcriptome that, in terms of non-redundant contigs (85,350) or different annotated sequences (15,367), is similar to those previously reported for other tissues in this fish species. The resulting information is accessible through an easy-to-use web interface hosted at http://www.nutrigroup-iats.org/seabassdb that contains data from assembled, annotated and semi-automatically curated singletons/contigs. This database has been built as previously reported for gilthead sea bream (http://www.nutrigroup-iats.org/seabreamdb) (Calduch-Giner et al., 2013), and offers the possibility for data requisition by several Blast options or direct word searches for annotated names or gene ontology terms. These genomic resources are useful tools for functional phenotyping of Mediterranean farmed fish and have been used to develop several pathway-focused PCR arrays for the simultaneous and semi-automated profiling of selected markers of mitochondria function and biogenesis (Bermejo-Nogales et al., 2014, 2015), lipid metabolism (Benedito-Palos et al., 2013, 2014), immune response (Pérez-Cordón et al., 2014; Pérez-Sánchez et al., 2015), and muscle growth (Benedito-Palos et al., 2016) in gilthead sea bream, and lipid metabolism in European sea bass (Rimoldi et al., 2016). In the current study, this knowledge has been used to design and validate a high density oligo-microarray, enriched in actively transcribed genes of the intestinal tract.
The advent of transcriptome analysis by RNA-seq has revolutionized the field of fish transcriptomics with the discovery of new transcripts, interrogation of post-transcriptional modifications and analysis of single nucleotide polymorphism (Louro et al., 2014; Qian et al., 2014). However, DNA microarray technology still constitutes a solid and widely used approach because of its competitiveness in terms of cost, sample preparation, and interpretation of data, despite of their methodological limitations including the dependence on the existing knowledge on expressed sequence tags or genomic sequences (Sedighi and Li, 2014; Martin et al., 2016). Indeed, the current custom European sea bass microarray was able to evidence a pronounced transcriptional divergence of PI in relation to AI-MI segments, which was validated by qPCR of 29 selected genes representative of the main biological process and functions across the intestine. The magnitude of the transcriptional difference (5770 genes, P < 0.05; 1906 with a fold-change cutoff = 2) between the intestinal segments found in the current study is higher than that reported in most of the studies dealing with changes induced in fish intestine by different experimental approaches. Table 4 shows the number of DE intestinal genes in fish with different feeding habits, exposed to environmental stress, pathogens, or diet interventions. Although comparisons across studies in different species is always difficult, transcriptome differences found between the two intestinal segments are of the same order than those found when comparing different organs and tissues, either in fish or other animal models. To put it in a proper perspective, the intestinal segment differences found in European sea bass are comparable to those found among heart and red and white skeletal muscle in gilthead sea bream, with more than 3000 DE genes (Calduch-Giner et al., 2014). Also, the level of DE probes among two different tissues in gene expression atlas of human and other model organisms is in the range of 50–90% (Son et al., 2005; Briggs et al., 2011; Freeman et al., 2012; Kapushesky et al., 2012).
Table 4. Fish transcriptomic profiling studies for the characterization of gut responses to environmental stress, pathogens, or diet interventions.
At the morphological level, differences across the intestinal length have been described in teleosts, with a progressive decrease in the diameter, the density of mucosal foldings and goblet cells (main cell producers of mucus), whereas rodlet cells closely linked to osmoregulation and the immune system become more abundant in the posterior segment (Reite, 2005; Ballester-Lozano et al., 2015). In European sea bass, a different profile of digestive enzyme activities between the anterior and the posterior segments has been shown, being the activity of amylase, alkaline phosphatase, total alkaline proteases, and trypsin higher at the AI segment (Castro et al., 2016). Functional and transcriptional compartmentalization of the digestive tract has also been described in zebrafish (Danio rerio), a stomachless teleost, with 2558 DE genes across the seven segments considered in that study (Wang et al., 2010). Thus, the fish intestine has evolved as a complex tissue with a gradual change in nutrient absorptive capacity from AI to PI segments, whereas the posterior segment appears as part of the first line of defense against pathogens (Løkka et al., 2014). The results reported herein indicate that this change is more gradual than initially envisaged, as some immune-relevant genes are much more expressed in the AI-MI segments than in PI, such as natterin-3-like and acidic mammalian chitinase. Natterin-like proteins contain a mannose-binding lectin-like domain that might act as a pathogen recognition protein with an important role in the acute-phase response of fish (Magnadóttir, 2006). Chitinases are primarily associated with the stomach, but they have also been found in the intestine with a role in food digestion and immune defense in several vertebrate species, including fish (Gutowska et al., 2004; Tran et al., 2011). Moreover, administration of recombinant chitinases boosts several serum immune parameters in orange-spotted grouper (Epinephelus coioides) (Zhang et al., 2012). All these results point out the importance of antimicrobial function also in this intestinal portion. On the other hand, the PI also plays a key role in the reabsorption of bile salts and some micronutrients. This dualism is, thereby, of relevance when considering the regulation of the intestine as whole, as further discussed below for some particular genes or biological processes.
Fish intestine has been reported to progressively tighten from the anterior to the posterior part, and intestinal integrity and intercellular permeability is maintained by proteins, such as integrins, claudins, cadherins, and gap junction proteins (Tsukita et al., 2001; Lu et al., 2013). Genes encoding these proteins were largely represented in our microarray, although no clear distinctive differential expression pattern was detected with the established cutoff value. However, the expression profile of these intestinal markers was quite similar to that reported in gilthead sea bream (Pérez-Sánchez et al., 2015). This finding suggests that the structure and function of tight junctions is maintained through the evolution of modern fish species, although some species differences can occur even when comparing close-related fish.
Functional network effector analysis clearly evidenced that active immune surveillance is a key role of the distal intestine of European sea bass, in comparison with the role of mammalian large intestine. This is probably due to fact that fish are continuously and directly exposed to a microbial-rich environment; e.g., 1 liter of marine water contains 10,000 million virus and 9000 million of bacteria (Fenical, 1993; Fuhrman, 1999). Thus, compared to terrestrial animals, they have to cope with high microbial loads. The immune cell repertoire of teleost intestine is governed by the gut-associated lymphoid tissue, and regional differences have already been reported (reviewed in Rombout et al., 2011; Salinas and Parra, 2015). In European sea bass, the predominance of T cells in the PI has been long established (Abelli et al., 1997), and the intestinal regionalization of the expression of T cell-markers was suggested for the first time in a teleost in this marine fish (Picchietti et al., 2011). Our present results confirm and expand these data with over-representation of several genes related to cell-mediated intestinal immunity in the PI segment. Similarly, a marked increase in neutrophil density was observed in intestinal segment 6 of zebrafish, in coincidence with an increase of the neutrophil-expressed genes and several cytokines implicated in intestinal immunity (Wang et al., 2010). Cathepsin L-1 also had the highest expression in segment 6, as cathepsin L-like in the PI of European sea bass. Cathepsin L family has multifunctional roles in many biochemical pathways of vertebrates, including intracellular protein degradation, antigen presentation, and cellular development (Zhou et al., 2015). Nevertheless, their function in fish is just starting to be elucidated and further studies are needed to reveal if they have an immunological role in the intestine of fish. In addition, infection models in a wide-range of fish species have shown the importance of the local intestinal immune response (Estensoro et al., 2012; Li et al., 2012; Pérez-Cordón et al., 2014; Dezfuli et al., 2015). In our study, the transcriptome of MI did not differ from that of AI, however, in a quantitative PCR analysis of 27 antimicrobial genes (Oehlers et al., 2011), the MI exhibited elevated expression of dual oxidase, the defensin beta-like, and peptidoglycan recognition protein families, and also presented the highest numbers of leukocytes and endocytic cells, supporting a specialized immunological role. In agreement with this, beta-defensin was among the top DE genes in the PI of European sea bass.
As expected, the highest expression level of genes related to digestion was observed in the AI-MI segments, since an absorption gradient is present along the length of the intestine of all fish. Water-soluble nutrients are mostly absorbed in the pyloric caeca and the midgut (Sundell and Rønnestad, 2011) and the majority of lipids are absorbed in these first segments with the concurrence of secreted bile salts to hydrolyse triglycerides to free fatty acids and glycerol. In coincidence with data from zebrafish (Wang et al., 2010), several well-known molecular markers of mammalian small intestine, fabp2, apoa1, and vil1l, had a higher expression at the AI-MI segments than at PI in European sea bass (Table S2). Furthermore, we also found that the gene markers (fabp2 and transcription factor gata5) found to be stable across the five most proximal segments of adult zebrafish (Oehlers et al., 2011), were higher in AI-MI segments than in the PI one. Typically, bile acids continue along the whole intestine and are absorbed by enterocytes of the PI segment and are transported to the liver as a part of the enterohepatic circulation. This process ensures the rapid turnover of bile acids by the participation of gastrotropin (also termed ileal lipid-binding protein or fatty acid-binding protein 6, fabp6), which is exclusively expressed in the intestine of mammals (Gong et al., 1994; Fujita et al., 1995; Besnard et al., 2002). The higher expression of gastrotropin in the PI segment of European sea bass is consistent with previous results in zebrafish (Alves-Costa et al., 2008; Oehlers et al., 2011) and gilthead sea bream (Pérez-Sánchez et al., 2015), in which the detection of gastrotropin transcripts was exclusive to the distal portions of intestine. The enterohepatic circulation also serves to efficiently use the water-soluble vitamin, vitamin B12. This vitamin plays a key role in the normal functioning of the brain and nervous system, and it is synthesized exclusively by bacteria in the gut (Raux et al., 2000). Deficiency of this vitamin results in decreased growth, anemia and abnormal fish behavior (NRC, 2011), and the dietary replacement of fish meal by plant proteins is accompanied by an increased risk of vitamin B12 deficiency in Atlantic cod (Hansen et al., 2007). The underlying regulatory mechanisms remain to be established in fish, but the observed high expression of vitamin B12 binders (gastric intrinsic factor-like protein, transcobalamin-2) and binder-vitamin complex receptors (protein amnionless, cubilin) in the PI of European sea bass ensures an efficient intestinal uptake and absorption of this essential micronutrient.
Our results also indicate a different expression pattern of GPCRs throughout the AI-MI and PI. GPCRs are the largest family of signaling receptors in vertebrates and most of them are ubiquitously expressed in the enteroendocrine cells of the GI tract. The presence of these receptors along the gut constitutes a key gastrointestinal chemosensory system for the regulation of appetite, nutrient digestion, intestinal motility, and mucosal defense mechanisms, contributing to the integration of a vast array of pathways linking the gut with the brain and metabolically active tissues (Reimann et al., 2012). In this study, the GPCR with the highest intestine spatial regulation (upregulated in the AI-MI) was GPR39. This is an orphan receptor that belongs to the ghrelin/motilin receptor subfamily, and it is emerging as an important regulator of gastrointestinal motility and secretion (Depoortere, 2012). Two forms of GPR39 have been characterized in black sea bream (Acanthopagrus schlegeli), and intestinal expression of sbGPR39-1a was found to decrease significantly during food deprivation (Zhang et al., 2008). Another orphan receptor, GPR112, was highly DE in the AI-MI of European sea bass. This GPCR was also highly expressed in the intestinal mucosal layer of zebrafish (Harty et al., 2015) and rodents (Ito et al., 2009; Badiali et al., 2012). NPBWR1 and NPBWR2 were also differentially upregulated in the AI-MI of European sea bass. Neuropeptide B and neuropeptide W are both endogenous peptide ligands for NPBWR1 and NPBWR2, which are key players in the regulation of feeding and energy metabolism (Mondal et al., 2003; Tanaka et al., 2003). These neuropeptides are mostly expressed in the brain, although a high expression of neuropeptide B has also been found in the intestine of medaka (Oryzias latipes) (Hiraki et al., 2014). Another GPCR with a high level of expression in the AI-MI of European sea bass was ADRA2C, a known ligand for norepinephrine and epinephrine (Lefkowitz and Caron, 1988). The regulatory role of alpha adrenergic receptors in the proliferation of intestinal epithelial crypt cells is well known in mammals (Schaak et al., 2000) and their tissue expression patterns parallel those of their mammalian orthologs in zebrafish (Ruuskanen et al., 2005) and perhaps, European sea bass.
Regarding the upregulated GPCRs in the PI segment, many of them were involved in cell proliferation and immune response. The exception was the orphan receptor GPR97, with a lowered expression in the PI of European sea bass despite its proven role in the remodeling of lymphatic cells of the mouse intestine (Valtcheva et al., 2013). This contrasted with the observation that FZD10, a receptor for molecules in the Wnt pathway that is a good marker of colorectal tumors in humans (Nagayama et al., 2009), was the most upregulated GPCR in the PI of European sea bass. Another highly and differentially expressed GPCR at PI was the hydrocarboxylic acid receptor 2, also termed GPR109. This receptor has anti-inflammatory properties acting as the metabolite sensor of butyrate, the end product of the microbial fermentation of dietary fiber (Singh et al., 2014). This short-chain fatty acid has been considered a promising feed additive in aquaculture, improving the immunological status and intestinal condition in carp (Cyprinus carpio) (Liu et al., 2014). Recent data also support the promising effects of butyrate supplementation in European sea bass, and Rimoldi and co-workers have found reported improved growth and increased gene expression of the oligopeptide transporter 1 (PEPT1) in the hindgut of fish fed low-fish meal diets (unpublished). Other GPCRs that showed increased expression levels in the PI of European sea bass were GPR84, GPR63, and GPR18. GPR84 is a putative receptor of medium-chain fatty acids that mediates proinflammatory responses in human myeloid cells (Suzuki et al., 2013), whereas GPR63 binds phosphatidic acid and is also involved in the inflammatory response (Kostenis, 2004). Lastly, GPR18, also known as N-arachidonyl glycine receptor, is ubiquitously expressed in channel catfish (Ictalurus punctatus) tissues, and its immunostimulatory action in response to infection with Aeromonas hydrophila has been proven (Pridgeon and Klesius, 2013).
In summary, the functional phenotyping at the molecular level of the intestinal tract of European sea bass has been assessed with a specific oligo-microarray. Similar molecular signatures were found for AI and MI segments, that were clearly different to PI one. The consistency of this finding was highly supported by qPCR results for the genes selected as representative markers of several functions or processes, including intestinal permeability, macronutrient and micronutrient digestion and absorption, immune response and intestinal chemosensing. With the available transcriptomic data, we can match the AI-MI segments of European sea bass to the small intestine of mammals and the first segments of zebrafish (segments 1–5 as defined in Wang et al., 2010, 1–4 as defined in Oehlers et al., 2011), whereas the PI of European sea bass would correspond to the large intestine of mammals and last segments of zebrafish (segments 6–7 as in Wang et al., 2010, or segment 6 as in Oehlers et al., 2011). In any case, fish gut transcriptome should not be considered a static feature, and significant changes and evolution could be envisaged not only in a spatial basis, but also related to changes in season and feeding regimes. In this sense, further work is necessary, not only in European sea bass, but in other fish species of interest, to establish the precise tuning of intestine plasticity from a transcriptional and functional perspective.
Author Contributions
AS and JP conceived and designed the experiment. JC sampled animals and performed RNA extractions. JC and JP analyzed microarray data. JC, AS, and JP wrote the manuscript.
Funding
This work was funded by the EU seventh Framework Programme by the ARRAINA (Advanced Research Initiatives for Nutrition and Aquaculture; KBBE-2011-288925) project. It does not necessarily reflect the views of the EU and in no way anticipates the Commission's future policy in this area. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. Additional funding was obtained from the Spanish Ministerio de Economía y Competitividad through the MI2-FISH project (Unraveling Metabolic, Intestinal and Immunopathological Fish Status; AGL2013-48560) and Generalitat Valenciana (PROMETEO FASE II-2014/085).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphys.2016.00359
References
Abelli, L., Picchietti, S., Romano, N., Mastrolia, L., and Scapigliati, G. (1997). Immunohistochemistry of gut-associated lymphoid tissue of the sea bass Dicentrarchus labrax (L.). Fish Shellfish Immunol. 7, 235–245. doi: 10.1006/fsim.1996.0079
Alves-Costa, F. A., Denovan-Wright, E. M., Thisse, C., Thisse, B., and Wright, J. M. (2008). Spatio-temporal distribution of fatty acid-binding protein 6 (fabp6) gene transcripts in the developing and adult zebrafish (Danio rerio). FEBS J. 275, 3325–3334. doi: 10.1111/j.1742-4658.2008.06480.x
Badiali, L., Cedernaes, J., Olszewski, P. K., Nylander, O., Vergoni, A. V., and Schiöth, H. B. (2012). Adhesion GPCRs are widely expressed throughout the subsections of the gastrointestinal tract. BMC Gastroenterol. 12:134. doi: 10.1186/1471-230X-12-134
Ballester-Lozano, G. F., Benedito-Palos, L., Estensoro, I., Sitjà-Bobadilla, A., Kaushik, S., and Pérez-Sánchez, J. (2015). Comprehensive biometric, biochemical and histopathological assessment of nutrient deficiencies in gilthead sea bream fed semi-purified diets. Br. J. Nutr. 114, 713–726. doi: 10.1017/S0007114515002354
Benedito-Palos, L., Ballester-Lozano, G. F., and Pérez-Sánchez, J. (2014). Wide-gene expression analysis of lipid-relevant genes in nutritionally challenged gilthead sea bream (Sparus aurata). Gene 547, 34–42. doi: 10.1016/j.gene.2014.05.073
Benedito-Palos, L., Ballester-Lozano, G. F., Simó, P., Karalazos, V., Ortiz, A., Calduch-Giner, J., et al. (2016). Lasting effects of butyrate and low FM/FO diets on growth performance, blood haematology/biochemistry and molecular growth-related markers in gilthead sea bream (Sparus aurata). Aquaculture 454, 8–18. doi: 10.1016/j.aquaculture.2015.12.008
Benedito-Palos, L., Calduch-Giner, J. A., Ballester-Lozano, G. F., and Pérez-Sánchez, J. (2013). Effect of ration size on fillet fatty acid composition, phospholipid allostasis and mRNA expression patterns of lipid regulatory genes in gilthead sea bream (Sparus aurata). Br. J. Nutr. 109, 1175–1187. doi: 10.1017/S000711451200311X
Bermejo-Nogales, A., Calduch-Giner, J. A., and Pérez-Sánchez, J. (2015). Unraveling the molecular signatures of oxidative phosphorylation to cope with the nutritionally changing metabolic capabilities of liver and muscle tissues in farmed fish. PLoS ONE 10:e0122889. doi: 10.1371/journal.pone.0122889
Bermejo-Nogales, A., Nederlof, M., Benedito-Palos, L., Ballester-Lozano, G. F., Folkedal, O., Olsen, R. E., et al. (2014). Metabolic and transcriptional responses of gilthead sea bream (Sparus aurata L.) to environmental stress: new insights in fish mitochondrial phenotyping. Gen. Comp. Endocrinol. 205, 305–315. doi: 10.1016/j.ygcen.2014.04.016
Besnard, P., Niot, I., Poirier, H., Clément, L., and Bernard, A. (2002). New insights into the fatty acid-binding protein (FABP) family in the small intestine. Mol. Cell. Biochem. 239, 139–147. doi: 10.1023/A:1020505512364
Briggs, J., Paoloni, M., Chen, Q. R., Wen, X., Khan, J., and Khanna, C. (2011). A compendium of canine normal tissue gene expression. PLoS ONE 6:e17107. doi: 10.1371/journal.pone.0017107
Cain, K., and Swan, C. (2010). “Barrier function and immunology,” in Fish Physiology Vol. 30. The Multifunctional Gut of Fish, eds M. Grosell, A. P. Farrell, and C. Brauner (San Diego, CA: Academic Press), 111–134.
Calduch-Giner, J. A., Bermejo-Nogales, A., Benedito-Palos, L., Estensoro, I., Ballester-Lozano, G., Sitjà-Bobadilla, A., et al. (2013). Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts. BMC Genomics 14:178. doi: 10.1186/1471-2164-14-178
Calduch-Giner, J. A., Echasseriau, Y., Crespo, D., Baron, D., Planas, J. V., Prunet, P., et al. (2014). Transcriptional assessment by microarray analysis and large-scale meta-analysis of the metabolic capacity of cardiac and skeletal muscle tissues to cope with reduced nutrient availability in gilthead sea bream (Sparus aurata L.). Mar. Biotechnol. 16, 423–435. doi: 10.1007/s10126-014-9562-3
Calduch-Giner, J. A., Sitjà-Bobadilla, A., Davey, G. C., Cairns, M. T., Kaushik, S., and Pérez-Sánchez, J. (2012). Dietary vegetable oils do not alter the intestine transcriptome of gilthead sea bream (Sparus aurata), but modulate the transcriptomic response to infection with Enteromyxum leei. BMC Genomics 13:470. doi: 10.1186/1471-2164-13-470
Castro, C., Couto, A., Pérez-Jiménez, A., Serra, C. R., Díaz-Rosales, P., Fernandes, R., et al. (2016). Effects of fish oil replacement by vegetable oil blend on digestive enzymes and tissue histomorphology of European sea bass (Dicentrarchus labrax) juveniles. Fish Physiol. Biochem. 42, 203–217. doi: 10.1007/s10695-015-0130-1
Conesa, A., Götz, S., García-Gómez, J. M., Perol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Davey, G., Calduch-Giner, J. A., Houeix, B., Talbot, A., Sitjà-Bobdailla, A., Prunet, P., et al. (2011). Molecular profiling of the gilthead sea bream (Sparus aurata) response to chronic exposure to the myxosporean parasite Enteromyxum leei. Mol. Immunol. 48, 2102–2112. doi: 10.1016/j.molimm.2011.07.003
De Santis, C., Bartie, K. L., Olsen, R. E., Taggart, J. B., and Tocher, D. R. (2015). Nutrigenomic profiling of transcriptional processes affected in liver and distal intestine in response to a soybean meal-induced stress in Atlantic salmon (Salmo salar). Comp. Biochem. Physiol. D Genomics Proteomics 15, 1–11. doi: 10.1016/j.cbd.2015.04.001
Depoortere, I. (2012). GI functions of GPR39: novel biology. Curr. Opin. Pharmacol. 12, 647–625. doi: 10.1016/j.coph.2012.07.019
Dezfuli, B. S., Bo, T., Lorenzoni, M., Shinn, A. P., and Giari, L. (2015). Fine structure and cellular responses at the host-parasite interface in a range of fish-helminth systems. Vet. Parasitol. 208, 272–279. doi: 10.1016/j.vetpar.2015.01.002
Doñate, C., Balasch, J. C., Callol, A., Bobe, J., Tort, L., and MacKenzie, S. (2010). The effects of immunostimulation through dietary manipulation in the rainbow trout; evaluation of mucosal immunity. Mar. Biotechnol. 12, 88–99. doi: 10.1007/s10126-009-9203-4
Estensoro, I., Ballester-Lozano, G., Benedito-Palos, L., Karalazos, V., Mallo, J. J., Sitjà-Bobadilla, A., et al. (2014). Histopathological and Transcriptional Scoring of Intestinal Traits in Gilthead Sea Bream (Sparus aurata) Fed Low Fish Meal and Fish Oil Diets with Butyrate Supplementation. San Sebastián: EAS Aquaculture.
Estensoro, I., Benedito-Palos, L., Palenzuela, O., Kaushik, S., Sitjà-Bobadilla, A., and Pérez-Sánchez, J. (2011). The nutritional background of the host alters the disease course in a fish-myxosporean system. Vet. Parasitol. 175, 141–150. doi: 10.1016/j.vetpar.2010.09.015
Estensoro, I., Calduch-Giner, J. A., Kaushik, S., Pérez-Sánchez, J., and Sitjà-Bobadilla, A. (2012). Modulation of the IgM gene expression and IgM immunoreactive cell distribution by the nutritional background in gilthead sea bream (Sparus aurata) challenged with Enteromyxum leei (Myxozoa). Fish Shellfish Immunol. 33, 401–410. doi: 10.1016/j.fsi.2012.05.029
Fenical, W. (1993). Chemical Studies of marine bacteria: developing a new resource. Chem. Rev. 93, 1673–1683. doi: 10.1021/cr00021a001
Fox, C. J., Hammerman, P. S., and Thompson, C. B. (2005). Fuel feeds function: energy metabolism and the T-cell response. Nat. Rev. Immunol. 5, 844–852. doi: 10.1038/nri1710
Freeman, T. C., Ivens, A., Baillie, J. K., Beraldi, D., Barnett, M. W., Dorward, D., et al. (2012). A gene expression atlas of the domestic pig. BMC Biol. 10:90. doi: 10.1186/1741-7007-10-90
Fuglem, B., Jirillo, E., Bjerkås, I., Kiyono, H., Nochi, T., Yuki, Y., et al. (2010). Antigen-sampling cells in the salmonid intestinal epithelium. Dev. Comp. Immunol. 34, 768–774. doi: 10.1016/j.dci.2010.02.007
Fuhrman, J. A. (1999). Marine viruses and their biogeochemical and ecological effects. Nature 399, 541–548. doi: 10.1038/21119
Fujita, M., Fujii, H., Kanda, T., Sato, E., Hatakeyama, K., and Ono, T. (1995). Molecular cloning, expression and characterization of a human intestinal 15-kDa protein. Eur. J. Biochem. 233, 406–413. doi: 10.1111/j.1432-1033.1995.406_2.x
Geay, F., Ferraresso, S., Zambonino-Infante, J. L., Bargelloni, L., Quentel, C., Vandeputte, M., et al. (2011). Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (Dicentrarchus labrax) half-subfamilies showing different growth rates with the plant-based diet. BMC Genomics 12:522. doi: 10.1186/1471-2164-12-522
Gong, Y. Z., Everett, E. T., Schwartz, D. A., Norris, J. S., and Wilson, F. A. (1994). Molecular cloning, tissue distribution, and expression of a 14-kDa bile acid-binding protein from rat ileal cytosol. Proc. Natl. Acad. Sci. U.S.A. 91, 4741–4745. doi: 10.1073/pnas.91.11.4741
Gutowska, M. A., Drazen, J. C., and Robison, B. H. (2004). Digestive chitinolytic activity in marine fishes of Monterey Bay, California. Comp. Biochem. Physiol. A 139, 351–358. doi: 10.1016/j.cbpb.2004.09.020
Hansen, A. C., Rosenlund, G., Karlsen, Ø., Koppe, W., and Hemre, G. I. (2007). Total replacement of fish meal with plant proteins in diets for Atlantic cod (Gadus morhua L.) I-Effects on growth and protein retention. Aquaculture 272, 599–611. doi: 10.1016/j.aquaculture.2007.08.034
Harty, B. L., Krishnan, A., Sanchez, N. E., Schiöth, H. B., and Monk, K. R. (2015). Defining the gene repertoire and spatiotemporal expression profiles of adhesion G protein-coupled receptors in zebrafish. BMC Genomics 16:62. doi: 10.1186/s12864-015-1296-8
Hiraki, T., Nakasone, K., Hosono, K., Kawabata, Y., Nagahama, Y., and Okubo, K. (2014). Neuropeptide B is female-specifically expressed in the telencephalic and preoptic nuclei of the medaka brain. Endocrinology 155, 1021–1032. doi: 10.1210/en.2013-1806
Inness, C. L. W., and Metcalfe, N. B. (2008). The impact of dietary restriction, intermittent feeding and compensatory growth on reproductive investment and lifespan in a short-lived fish. Proc. Biol. Sci. 275, 1703–1708. doi: 10.1098/rspb.2008.0357
Ito, J., Ito, M., Nambu, H., Fujikawa, T., Tanaka, K., Iwaasa, H., et al. (2009). Anatomical and histological profiling of orphan G-protein-coupled receptor expression in gastrointestinal tract of C57BL/6J mice. Cell Tissue Res. 338, 257–269. doi: 10.1007/s00441-009-0859-x
Kaitetzidou, E., Xiang, J., Antonopoulou, E., Tsigenopoulos, C. S., and Sarropoulou, E. (2015). Dynamics of gene expression patterns during early development of the European sea bass (Dicentrarchus labrax). Physiol. Genomics 47, 158–169. doi: 10.1152/physiolgenomics.00001.2015
Kapushesky, M., Adamusiak, T., Burdett, T., Culhane, A., Farne, A., Filippov, A., et al. (2012). Gene Expression Atlas update: a value-added database of microarray and sequencing-based functional genomics experiments. Nucleic Acids Res. 40, D1077–D1081. doi: 10.1093/nar/gkr913
Kau, A. L., Ahern, P. P., Griffin, N. W., Goodman, A. L., and Gordon, J. I. (2011). Human nutrition, the gut microbiome and the immune system. Nature 474, 327–336. doi: 10.1038/nature10213
Kostenis, E. (2004). Novel clusters of receptors for sphingosine-1-phosphate, sphingo-sylphosphorylcholine and (lyso)-phosphatidic acid: new receptors for “old” ligands. J. Cell. Biochem. 92, 923–936. doi: 10.1002/jcb.20092
Król, E., Douglas, A., Tocher, D. R., Crampton, V. O., Speakman, J. R., Secombes, C. J., et al. (2016). Differential responses of the gut transcriptome to plant protein diets in farmed Atlantic salmon. BMC Genomics 17:156. doi: 10.1186/s12864-016-2473-0
Lefkowitz, R. J., and Caron, M. G. (1988). Adrenergic receptors. Models for the study of receptors coupled to guanine nucleotide regulatory proteins. J. Biol. Chem. 263, 4993–4996.
Li, C., Zhang, Y., Wang, R., Lu, J., Nandi, S., Mohanty, S., et al. (2012). RNA-seq analysis of mucosal immune response reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus. Fish Shellfish Immunol. 32, 816–827. doi: 10.1016/j.fsi.2012.02.004
Li, L., Liang, X. F., He, S., Sun, J., Wen, Z. Y., He, Y. H., et al. (2015). Transcriptome analysis of grass carp (Ctenopharyngodon idella) fed with animal and plant diets. Gene 574, 371–379. doi: 10.1016/j.gene.2015.08.030
Liu, W., Yang, Y., Zhang, J., Gatlin, D. M., Ringo, E., and Zhou, Z. (2014). Effects of microencapsulated sodium butyrate on growth, intestinal mucosal morphology, immune response and adhesive bacteria in juvenile common carp (Cyprinus carpio) pre-fed with or without oxidized oil. Br. J. Nutr. 112, 15–29. doi: 10.1017/S0007114514000610
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Løkka, G., Austbø, L., Falk, K., Bromage, E., Fjelldal, P. G., Hansen, T., et al. (2014). Immune parameters in the intestine of wild and reared unvaccinated and vaccinated Atlantic salmon (Salmo salar L.). Dev. Comp. Immunol. 47, 6–16. doi: 10.1016/j.dci.2014.06.009
Louro, B., Power, D. M., and Canario, A. V. (2014). Advances in European sea bass genomics and future perspectives. Mar. Genomics 18, 71–75. doi: 10.1016/j.margen.2014.06.010
Lu, Z., Ding, L., Lu, Q., and Chen, Y. H. (2013). Claudins in intestines: distribution and functional significance in health and diseases. Tissue Barriers 1:e24978. doi: 10.4161/tisb.24978
Magnadóttir, B. (2006). Innate immunity of fish (overview). Fish Shellfish Immunol. 20, 137–151. doi: 10.1016/j.fsi.2004.09.006
Magnanou, E., Klopp, C., Noirot, C., Besseau, L., and Falcón, J. (2014). Generation and characterization of the sea bass Dicentrarchus labrax brain and liver transcriptomes. Gene 544, 56–66. doi: 10.1016/j.gene.2014.04.032
Martin, S. A. M., Dehler, C. E., and Król, E. (2016). Transcriptomic responses in the fish intestine. Dev. Comp. Immunol. 64, 103–117. doi: 10.1016/j.dci.2016.03.014
Michalek, R. D., and Rathmell, J. C. (2010). The metabolic life and times of a T-cell. Immunol. Rev. 236, 190–202. doi: 10.1111/j.1600-065X.2010.00911.x
Mondal, M. S., Yamaguchi, H., Date, Y., Shimbara, T., Toshinai, K., Shimomura, Y., et al. (2003). A role for neuropeptide W in the regulation of feeding behavior. Endocrinology 144, 4729–4733. doi: 10.1210/en.2003-0536
Morais, S., Edvardsen, R. B., Tocher, D. R., and Bell, J. G. (2012a). Transcriptomic analyses of intestinal gene expression of juvenile Atlantic cod (Gadus morhua) fed diets with Camelina oil as replacement for fish oil. Comp. Biochem. Physiol. B 161, 282–293. doi: 10.1016/j.cbpb.2011.12.004
Morais, S., Silva, T., Cordeiro, O., Rodrigues, P., Guy, D. R., Bron, J. E., et al. (2012b). Effects of genotype and dietary fish oil replacement with vegetable oil on the intestinal transcriptome and proteome of Atlantic salmon (Salmo salar). BMC Genomics 13:448. doi: 10.1186/1471-2164-13-448
Murray, H. M., Lall, S. P., Rajaselvam, R., Boutilier, L. A., Blanchard, B., Flight, R. M., et al. (2010). A nutrigenomic analysis of intestinal response to partial soybean meal replacement in diets for juvenile Atlantic halibut, Hippoglossus hippoglossus, L. Aquaculture 298, 282–293. doi: 10.1016/j.aquaculture.2009.11.001
National Research Council (NRC) (2011). Nutrient Requirement of Fish and Shellfish. Washington, DC: National Academy Press.
Nagayama, S., Yamada, E., Kohno, Y., Aoyama, T., Fukukawa, C., Kubo, H., et al. (2009). Inverse correlation of the up-regulation of FZD10 expression and the activation of ß-catenin in synchronous colorectal tumors. Cancer Sci. 100, 405–412. doi: 10.1111/j.1349-7006.2008.01052.x
Nguyen, T. V., Jung, H., Nguyen, T. M., Hurwood, D., and Mather, P. (2016). Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach. Mar. Genomics 25, 75–88. doi: 10.1016/j.margen.2015.11.010
Nuñez Ortiz, N., Gerdol, M., Stocchi, V., Marozzi, C., Randelli, E., Bernini, C., et al. (2014). T cell transcripts and T cell activities in the gills of the teleost fish sea bass (Dicentrarchus labrax). Dev. Comp. Immunol. 47, 309–318. doi: 10.1016/j.dci.2014.07.015
Oehlers, S. H., Flores, M. V., Chen, T., Hall, C. J., Crosier, K. E., and Crosier, P. S. (2011). Topological distribution of antimicrobial genes in the zebrafish intestine. Fish Shellfish Immunol. 35, 385–391. doi: 10.1016/j.dci.2010.11.008
Pérez-Cordón, G., Estensoro, E., Benedito-Palos, L., Calduch-Giner, J. A., Sitjà-Bobadilla, A., and Pérez-Sánchez, J. (2014). Interleukin gene expression is strongly modulated at the local level in a fish-parasite model. Fish Shellfish Immunol. 37, 201–208. doi: 10.1016/j.fsi.2014.01.022
Pérez-Sánchez, J., Benedito-Palos, L., Estensoro, I., Petropoulos, Y., Calduch-Giner, J. A., Browdy, C. L., et al. (2015). Effects of dietary NEXT ENHANCE 150 on growth performance and expression of immune and 2 intestinal integrity related genes in gilthead sea bream (Sparus aurata L.). Fish Shellfish Immunol. 44, 117–128. doi: 10.1016/j.fsi.2015.01.039
Pérez-Sánchez, J., Estensoro, I., Redondo, M. J., Calduch-Giner, J. A., Kaushik, S., and Sitjà-Bobadilla, A. (2013). Mucins as diagnostic and prognostic biomarkers in a fish-parasite model: transcriptional and functional analysis. PLoS ONE 8:e65457. doi: 10.1371/journal.pone.0065457
Picchietti, S., Guerra, L., Bertoni, F., Randelli, E., Belardinelli, M. C., Buonocore, F., et al. (2011). Intestinal T cells of Dicentrarchus labrax (L.): gene expression and functional studies. Fish Shellfish Immunol. 30, 609–617. doi: 10.1016/j.fsi.2010.12.006
Pridgeon, J. W., and Klesius, P. H. (2013). G-protein coupled receptor 18 (GPR18) in channel catfish: expression analysis and efficacy as immunostimulant against Aeromonas hydrophila infection. Fish Shellfish Immunol. 35, 1070–1078. doi: 10.1016/j.fsi.2013.07.017
Qian, X., Ba, Y., Zhuang, Q., and Zhong, G. (2014). RNA-Seq technology and its application in fish transcriptomics. Omics 18, 98–110. doi: 10.1089/omi.2013.0110
Raux, E., Schubert, H. L., and Warren, M. J. (2000). Biosynthesis of cobalamin (vitamin B12): a bacterial conundrum. Cell. Mol. Life Sci. 57, 1880–1893. doi: 10.1007/PL00000670
Reimann, F., Tolhurst, G., and Gribble, F. M. (2012). G-protein coupled receptors in intestinal chemosensation. Cell Metab. 15, 421–431. doi: 10.1016/j.cmet.2011.12.019
Reite, O. B. (2005). The rodlet cells of teleostean fish: their potential role in host defence in relation to the role of mast cells/eosinophilic granule cells. Fish Shellfish Immunol. 19, 253–267. doi: 10.1016/j.fsi.2005.01.002
Rimoldi, S., Benedito-Palos, L., Terova, G., and Pérez-Sánchez, J. (2016). Wide-targeted gene expression approach to infer the tissue-specific molecular signatures of lipid metabolism in fed and fasted juvenile European sea bass (Dicentrarchus labrax). Rev. Fish Biol. Fish. 26, 93–108. doi: 10.1007/s11160-015-9408-8
Robledo, D., Ronza, P., Harrison, P. W., Losada, A. P., Bermúdez, R., Pardo, B. J., et al. (2014). RNA-seq analysis reveals significant transcriptome changes in turbot (Scophthalmus maximus) suffering severe enteromyxosis. BMC Genomics 15:1149. doi: 10.1186/1471-2164-15-1149
Rombout, J. H. W. M., Abelli, L., Picchietti, S., Scapigliati, G., and Kiron, V. (2011). Teleost intestinal immunology. Fish Shellfish Immunol. 31, 616–626. doi: 10.1016/j.fsi.2010.09.001
Ronkin, D., Seroussi, E., Nitzan, T., Doron-Faigenboim, A., and Cnaani, A. (2015). Intestinal transcriptome analysis revealed differential salinity adaptation between two tilapiine species. Comp. Biochem. Physiol. D 13, 35–43. doi: 10.1016/j.cbd.2015.01.003
Rurangwa, E., Sipkema, D., Kals, J., Ter Veld, M., Forlenza, M., Bacanu, G. M., et al. (2015). Impact of a novel protein meal on the gastrointestinal microbiota and the host transcriptome of larval zebrafish Danio rerio. Front. Physiol. 6:133. doi: 10.3389/fphys.2015.00133
Ruuskanen, J. O., Peitsaro, N., Kaslin, J. V., Panula, P., and Scheinin, M. (2005). Expression and function of alpha-adrenoceptors in zebrafish: drug effects, mRNA and receptor distributions. J. Neurochem. 94, 1559–1569. doi: 10.1111/j.1471-4159.2005.03305.x
Sahlmann, C., Sutherland, B. J., Kortner, T. M., Koop, B. F., Krogdahl, Å., and Bakke, A. M. (2013). Early response of gene expression in the distal intestine of Atlantic salmon (Salmo salar L.) during the development of soybean meal induced enteritis. Fish Shellfish Immunol. 34, 599–609. doi: 10.1016/j.fsi.2012.11.031
Salinas, I., and Parra, D. (2015). “Fish mucosal immunity: intestine,” in Mucosal Health in Aquaculture, eds B. H. Beck and E. Peatman (London: Academic Press), 135–171.
Sarropoulou, E., Galindo-Villegas, J., Garcia-Alcazar, A., Kasapidis, P., and Mulero, V. (2012). Characterization of European sea bass transcripts by RNA SEQ after oral vaccine against V. anguillarum. Mar. Biotechnol. 14, 634–642. doi: 10.1007/s10126-012-9466-z
Schaak, S., Cussac, D., Cayla, C., Devedjian, J. C., Guyot, R., Paris, H., et al. (2000). Alpha(2) adrenoceptors regulate proliferation of human intestinal epithelial cells. Gut 47, 242–250. doi: 10.1136/gut.47.2.242
Sedighi, A., and Li, P. C. H. (2014). “Challenges and future trends in DNA microarray analysis,” in Fundamentals of Advanced Omics Technologies: From Genes to Metabolites, eds C. Simó, A. Cifuentes, and V. García-Cañas (Philadelphia, PA: Elsevier), 25–46.
Shi, M., Huang, R., Du, F., Pei, Y., Liao, L., Zhu, Z., et al. (2014). RNA-seq profiles from grass carp tissues after reovirus (GCRV) infection based on singular and modular enrichment analyses. Mol. Immunol. 61, 44–53. doi: 10.1016/j.molimm.2014.05.004
Singh, N., Gurav, A., Sivaprakasam, S., Brady, E., Padia, R., Shi, H., et al. (2014). Activation of GPR109a, receptor for niacin and the commensal metabolite butyrate, suppresses colonic inflammation and carcinogenesis. Immunity 40, 128–139. doi: 10.1016/j.immuni.2013.12.007
Skugor, S., Grisdale-Helland, B., Refstie, S., Afanasyev, S., Vielma, J., and Krasnov, A. (2011). Gene expression responses to restricted feeding and extracted soybean meal in Atlantic salmon (Salmo salar L.). Aquacult. Nutr. 17, 505–517. doi: 10.1111/j.1365-2095.2010.00832.x
Son, C. G., Bilke, S., Davis, S., Greer, B. T., Wei, J. S., Whiteford, C. C., et al. (2005). Database of mRNA gene expression profiles of multiple human organs. Genome Res. 15, 443–450. doi: 10.1101/gr.3124505
Sundell, K. S., and Rønnestad, I. (2011). “Intestinal absorption,” in Encyclopedia of Fish Physiology, ed P. F. Anthony (San Diego, CA: Academic Press), 1311–1321.
Suzuki, M., Takaishi, S., Nagasiki, M., Onozawa, Y., Iino, I., Maeda, H., et al. (2013). Medium-chain fatty-acid sensing receptor, GPR84, is a proinflammatory receptor. J. Biol. Chem. 288, 10684–10691. doi: 10.1074/jbc.M112.420042
Tanaka, H., Yoshida, T., Miyamoto, N., Motoike, T., Kurosu, H., Shibata, K., et al. (2003). Characterization of a family of endogenous neuropeptide ligands for the G protein-coupled receptors GPR7 and GPR8. Proc. Natl. Acad. Sci. U.S.A. 100, 6251–6256. doi: 10.1073/pnas.0837789100
Tine, M., Kuhl, H., Gagnaire, P.-A., Louro, B., Desmarais, E., Martins, R. S. T., et al. (2014). European sea bass genome and its variation provide insights into adaptation to euryhalinity and speciation. Nat. Commun. 5:5770. doi: 10.1038/ncomms6770
Tran, H. T., Barnich, N., and Mizoguchi, E. (2011). Potential role of chitinases and chitin-binding proteins in host-microbial interactions during the development of intestinal inflammation. Histol. Histopathol. 26, 1453–1464. doi: 10.14670/HH-26.1453
Tsukita, S., Furuse, M., and Itoh, M. (2001). Multifunctional strands in tight junctions. Nat. Rev. Mol. Cell Biol. 2, 285–293. doi: 10.1038/35067088
Valtcheva, N., Primorac, A., Jurisic, G., Hollmén, M., and Detmar, M. (2013). The orphan adhesion G-protein coupled receptor GPR97 regulates migration of lymphatic endothelial cells via the small GTPases RhoA and Cdc42. J. Biol. Chem. 288, 35736–35748. doi: 10.1074/jbc.M113.512954
Wang, Z., Du, J., Lam, S. H., Mathavan, S., Matsudaira, P., and Gong, Z. (2010). Morphological and molecular evidence for functional organization along the rostrocaudal axis of the adult zebrafish intestine. BMC Genomics 11:192. doi: 10.1186/1471-2164-11-392
Wilson, J. M., and Castro, L. F. C. (2010). “Morphological diversity of the gastrointestinal tract in fishes,” in Fish Physiology Vol. 30. The Multifunctional Gut of Fish, eds M. Grosell, A. P. Farrell, and C. Brauner (San Diego, CA: Academic Press), 1–55.
Wong, M. K., Ozaki, H., Suzuki, Y., Iwasaki, W., and Takei, Y. (2014). Discovery of osmotic sensitive transcription factors in fish intestine via a transcriptomic approach. BMC Genomics 15:1134. doi: 10.1186/1471-2164-15-1134
Xia, J. H., Liu, P., Liu, F., Lin, G., Sun, F., Tu, R. J., et al. (2013). Analysis of stress-responsive transcriptome in the intestine of Asian seabass (Lates calcarifer) using RNA-Seq. DNA Res. 20, 449–460. doi: 10.1093/dnares/dst022
Zhang, Y., Feng, S., Chen, J., Qin, C., Lin, H., and Li, W. (2012). Stimulatory effects of chitinase on growth and immune defense of orange-spotted grouper (Epinephelus coioides). Fish Shellfish Immunol. 32, 844–854. doi: 10.1016/j.fsi.2012.02.009
Zhang, Y., Liu, Y., Huang, X., Liu, X., Jiao, B., Meng, Z., et al. (2008). Two alternatively spliced GPR39 transcripts in seabream: molecular cloning, genomic organization, and regulation of gene expression by metabolic signals. J. Endocrinol. 199, 457–470. doi: 10.1677/JOE-07-0608
Zhou, J., Zhang, Y. Y., Li, Q. Y., and Cai, Z. H. (2015). Evolutionary history of cathepsin L (L-like) family genes in vertebrates. Int. J. Biol. Sci. 11, 1016–1025. doi: 10.7150/ijbs.11751
Keywords: European sea bass, intestine, transcriptome, next-generation sequencing, microarray, immune response
Citation: Calduch-Giner JA, Sitjà-Bobadilla A and Pérez-Sánchez J (2016) Gene Expression Profiling Reveals Functional Specialization along the Intestinal Tract of a Carnivorous Teleostean Fish (Dicentrarchus labrax). Front. Physiol. 7:359. doi: 10.3389/fphys.2016.00359
Received: 14 June 2016; Accepted: 05 August 2016;
Published: 25 August 2016.
Edited by:
Rubina Sirri, University of Bologna, ItalyCopyright © 2016 Calduch-Giner, Sitjà-Bobadilla and Pérez-Sánchez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jaume Pérez-Sánchez, amFpbWUucGVyZXouc2FuY2hlekBjc2ljLmVz