The Ontogeny and Brain Distribution Dynamics of the Appetite Regulators NPY, CART and pOX in Larval Atlantic Cod (Gadus morhua L.)

Abstract

Similar to many marine teleost species, Atlantic cod undergo remarkable physiological changes during the early life stages with concurrent and profound changes in feeding biology and ecology. In contrast to the digestive system, very little is known about the ontogeny and the localization of the centers that control appetite and feed ingestion in the developing brain of fish. We examined the expression patterns of three appetite regulating factors (orexigenic: neuropeptide Y, NPY; prepro-orexin, pOX and anorexigenic: cocaine- and amphetamine-regulated transcript, CART) in discrete brain regions of developing Atlantic cod using chromogenic and double fluorescent in situ hybridization. Differential temporal and spatial expression patterns for each appetite regulator were found from first feeding (4 days post hatch; dph) to juvenile stage (76 dph). Neurons expressing NPY mRNA were detected in the telencephalon (highest expression), diencephalon, and optic tectum from 4 dph onward. CART mRNA expression had a wider distribution along the anterior-posterior brain axis, including both telencephalon and diencephalon from 4 dph. From 46 dph, CART transcripts were also detected in the olfactory bulb, region of the nucleus of medial longitudinal fascicle, optic tectum and midbrain tegmentum. At 4 and 20 dph, pOX mRNA expression was exclusively found in the preoptic region, but extended to the hypothalamus at 46 and 76 dph. Co-expression of both CART and pOX genes were also observed in several hypothalamic neurons throughout larval development. Our results show that both orexigenic and anorexigenic factors are present in the telencephalon, diencephalon and mesencephalon in cod larvae. The telencephalon mostly contains key factors of hunger control (NPY), while the diencephalon, and particularly the hypothalamus may have a more complex role in modulating the multifunctional control of appetite in this species. As the larvae develop, the overall progression in temporal and spatial complexity of NPY, CART and pOX mRNAs expression might be correlated to the maturation of appetite control regulation. These observations suggest that teleost larvae continue to develop the regulatory networks underlying appetite control after onset of exogenous feeding.

About the Authors:

Hoang T. M. D. Le

Current address: Institute of Aquaculture, Nha Trang University, Nha Trang, Khanh Hoa, Vietnam

Affiliation: Department of Biology, University of Bergen, Bergen, Norway

Anna Rita Angotzi

Affiliation: Department of Biology, University of Bergen, Bergen, Norway

Lars O. E. Ebbesson

Affiliation: Uni Research Environment, Uni Research AS, Bergen, Norway

Ørjan Karlsen

Affiliation: Institute of Marine Research, Storebø, Norway

Ivar Rønnestad

* E-mail: [email protected]

Affiliation: Department of Biology, University of Bergen, Bergen, Norway

Introduction

Atlantic cod (Gadus morhua) has a key role in North Atlantic fisheries and several major attempts to introduce cod into aquaculture have occurred over the last decades [1, 2]. Following the recent sequencing of the Atlantic cod genome [3], this species is also being explored as a model for understanding the genetic control of development and physiology in relation to environmental and nutritional factors [4]. During the early life stages, Atlantic cod undergo significant changes in feeding biology. At the onset of first feeding (4 dph), the larvae are 4–5 mm long and possess a simple and straight gut lumen. The intestine develops both in length and complexity including the development of elaborate folding of the mucosa layer. The stomach starts to differentiate from the foregut at about 25 dph [5, 6] although it takes a long time (>50 dph) before it becomes fully functional [7]. Larvae have been observed to continue to feed despite a full gut, especially just after onset of exogenous feeding [8]. These observations raise the question whether appetite control (particularly satiety) and digestive processes are fully functional in early larval stages [8]. The neuronal network underlying appetite control in vertebrates includes both orexigenic and anorexigenic signals, mainly neuropeptides produced in the brain, where the hypothalamus plays a pivotal role in the regulation of feeding [9–11]. Neuropeptide Y (NPY) appears to stimulate feeding behavior in several studies on teleost fish [11]. In Atlantic cod, RT-PCR analysis on whole larvae demonstrated NPY mRNA expression as early as 3 dph and a decrease during further development with the final assessment in larvae aged 60 dph [12]. Cocaine- and amphetamine-regulated transcript (CART) is a well-known anorexigenic factor that mediates feeding behavior and body-weight control in both mammals [13] and fish [11]. In Atlantic cod, the first detection of CART mRNA expression was found at 3 dph [12], i.e. at onset of first exogenous feeding. The Orexin A (hypocretin-1) and Orexin B (hypocretin-2) are two orexigenic neuropeptides produced from a single gene precursor known as prepro-orexin (pOX) [14]. In the teleost fish investigated so far, the onset of pOX expression occurs at very early embryonic stages [15]. Consistently, RT-PCR studies in Atlantic cod revealed pOX expression as early as cleavage stage [16].

Although recent studies have described the distribution of the neuropeptides NPY, CART, and pOX in the brain of adult Atlantic cod and other teleost fish [17–21], the spatial and temporal mRNA distribution of these genes during the critical first feeding stages and subsequent development are unknown. In the present study, we describe the mRNA-expression patterns of NPY, CART and pOX in the brain of Atlantic cod from first feeding (4dph) to 76 dph, with emphasis on homologous regions to mammalian appetite control regions located in both forebrain and midbrain. We further investigated coexpression of orexigenic and anorexigenic factors (NPY/CART, and pOX/CART) in these brain regions.

Materials and Methods

Fish material and sampling

The study was performed at Austevoll Research station (outside Bergen, Norway). The larvae originated from five-year old Atlantic cod broodstock (3–5 kg), grand-offspring of wild fish raised in net-pens at ambient temperature and salinity in Parisvannet, Norway (60°37’N, 4°48’E). The broodstocks were kept at 1:1 sex ratio and fed ad libitum on artificial feed three times a week. Eggs stripped from females were fertilized and incubated at 6°C under a 12L:12D light regime. The larvae were transferred to start feeding units and maintained according to best practice protocols at Institute of Marine Research (IMR) [22] at 8°C with a 16L:8D photoperiod and fed three meals of zooplankton per day [22]. The zooplanktons were filtered from a nearby pond, and from 4 to 34 days post hatching (dph) consisted of nauplii and copepodites (5–34 ind./mL), and from 35 to 46 dph mainly copepodities (3–7 ind./mL). From 47 to 55 dph, the larvae were fed with copepodities (3–7 ind./mL) and started weaning by AgloNorse^® 300–500 μm (Tromsø Fiskeindustri AS, Tromsø, Norway). From 56 dph, larvae were fed on only AgloNorse^® 400–600 μm. Experiments were carried out in triplicate tanks.

The mean standard length of the larvae sampled at 4, 20, 46 and 76 dph were 4.5, 5.9, 14.1 and 37.6 mm, indicating that they were in the Stages 1, 2, 4 and juvenile according to Sæle et al. (Ø. Sæle, NIFES, Bergen, Norway, personal communication). Cod larvae in Stage 1 are characterized an open mouth, pigmented eyes and the utilization about ¾ of the yolk sac has been used. Stage 2 is prior to metamorphosis, the ossification in the head is dominated by the jaws, and some ossification is also observed in the hyoid arch, the gill arches and some parts of the neurocranium. In Stage 4, after metamorphosis, the fins are developed but unpigmented, flexion is present, and most of the bones in the head are being ossified. The stomach is differentiated and the larval fin fold has disappeared. In the juvenile fish the pigmentation is complete.

Identification and cloning of NPY, CART and pOX genes

NPY, CART, and pOX sequences were retrieved from GeneBank sequence database (NCBI) with accession number AY822596, DQ167209, and DQ486137, respectively. cDNA sequence fragments of NPY (346 bp), CART (300 bp), and pOX (453bp) were synthesized from total RNA extracted from cod brain, using oligo(dT)12-18 primer and Superscript III (Invitrogen, Carlsbad, CA, USA), by standard procedures. The primers used were: ggaactctgaccgagggat NPY(F) and gccctctgatgacaaatca NPY(R) for NPY; gagtgtggaccagagccttg CART(F) and gcagtcacacatcttcccaat CART(R) for CART; aatgaagtggtcctccacagtgt Orex(F) and tcaagcggtgaagtcttgctgc Orex(R) for pOX. PCR amplification was performed with an initial step of 94°C for 5 min, 30 cycles of 94°C for 25s, 55°C for 30s, 68°C for 90s and extension at 68°C for 7 min. PCR products were purified using QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany) and cloned into TOPO cloning vector (Invitrogen). The inserts were sequenced at the University of Bergen Sequencing Facility, and used to synthesize antisense and sense (control) probes with digoxigenin (DIG RNA labeling mix, Roche Diagnostic) and fluorescein (Fluorescein RNA Labeling Mix, Roche Diagnostic) haptens for ISH methods.

Tissue preparation for cryosectioning

Cod larvae (4, 20, 46, 76 dph; correlating to Stages 1, 2, 4 and juvenile, respectively; Sæle et al., unpublished) were killed by an anesthetic overdose of MS– 222 and fixed in 4% paraformaldehyde (PF) in 0.1 M Sørensen´s phosphate buffer (28 mM NaH₂PO₄, 72 mM Na₂HPO₄, pH 7.2) at 4°C for 30 to 48 hours depending on the larval size. Samples were rinsed in 1X phosphate buffered saline (PBS, pH 7.4) and treated overnight in 25% sucrose in 1X PBS/25% Tissue-Tek (O.C.T) Compound (Sakura Finetek Europe B.V, Netherland) at 4°C. Fixed larvae were then embedded in 100% O.C.T compound and maintained at– 80°C until sectioning. Transverse sections were cut at thickness of 10–12 μm using Leica CM 1850 cryostat, stretched on glass, dried at 65°C and stored at -80°C until used.

In situ hybridization

Chromogenic in situ hybridization (CISH) was used to detect the location of NPY, CART, and pOX expression separately. Dual fluorescent in situ hybridization (FISH) was carried out on gene pairs (NPY/CART and pOX/CART) in the same sections. Hybridizations with NPY, CART and pOX digoxigenin or fluorescein labeled antisense riboprobes were revealed with either colorimetric NBT/BCIP staining (Roche Diagnostic), or fluorescent TSA and Fast Red staining (TSA PLUS Fluorescein Kit, PerkinElmer and Fast Red Roche Diagnostic, respectively). ISH was performed as described by Sandbakken M et. al [23], with some modifications as described below.

Pre-hybridization treatment.

Transverse brain sections were air-dried, rehydrated and washed in 2X saline-sodium citrate (2X SSC), followed by proteinase K treatment (10 μg/mL in 100 mM Tris-HCl pH 8.0 and 50 mM ethylenediamine tetracetic acid (EDTA) pH 8.0) and post-fixation with 4% PF in 1X PBS (2 min each). Samples were then treated with 0.25% acetic anhydride in 0.1 M TEA (10 min) prior washing in 2X SSC (1 min), before incubation with riboprobes.

Hybridization.

For CISH, 100 μl of hybridization solution (10% dextran sulphate, 10 mM Tris-HCl pH 7.5, 300 mM NaCl, 20 mM EDTA pH 8.0, 50% deionized formamide (Sigma-Aldrich), 0.2% tween-20, 1% blocking solution) with 200 ng digoxigenin or fluorescent labeled probes was applied to each slide and incubated overnight at 65°C in humid chamber with 2X SSC.

Post-hybridization.

Following overnight incubation, sections were washed in 2X SSC (2x30 min) and treated with 50% deionized formamide in 2X SSC (65°C, 30 min), prior RNase treatment (RNase A 0.02 mg/ml in 10 mM Tris-HCl pH 7.5, 500 mM NaCl, 1 mM EDTA) (37°C, 30 min).

Color detection.

For CISH, slides were washed with 2% blocking solution in 2X SSC- 0.05% triton x-100 (2 hours) and in 1X maleate buffer (100 mM maleic acid, 150 mM NaCl, pH 7.5; 2x5 min), before overnight incubation with 1:2000 anti-digoxigenin (FAB fragments) conjugated to alkaline phosphatase antibody (Roche Diagnostic) in antibody buffer (100mM maleic acid, 1% blocking solution, 0.3% triton x-100). Redundant antibody was removed by washes in 1X maleate buffer (2x10 min). Sections were washed in visualization buffer (100nM Tris-HCl pH 9.5, 5 mM MgCl₂, 100 mM NaCl) for 10 minutes before being incubated with chromogen substrates Nitroblue Tetrazolium (NBT) and Bromo-4-chloro-3-indoylphosphate (BCIP; Roche Diagnostic) to develop the staining and then washed in stop buffer (10 mM Tris-HCl pH 7.5, 1 mM EDTA pH 8.0, 150 mM NaCl, pH 7.5) prior mounting in 60% glycerol in stop buffer, and stored at 4°C until visualization by Leica DM6000B microscope (4 and 20 dph) and Leica M420 microscope (46 and 76 dph) equipped with digital cameras Leica DFC 420 and Cool Snap-Pro color, respectively.

For dual florescence detection, sections were incubated overnight at 4°C with anti-fluorescein FAB fragments conjugated to horseradish peroxidase (POD) antibody (Roche Diagnostic) in antibody buffer (1:250). After washing in 1X maleate buffer, sections were rinsed in TNT washing buffer (100 mM Tris-HCl pH 7.5, 150 mM NaCl, 0.05% tween20) and incubated overnight at 4°C in darkness with fluorescein amplification reagent tyramide (TSA PLUS Fluorescein Kit, PerkinElmer) in TNT buffer (1:200). The reaction was then stopped in TNT buffer and the sections washed in immuno-buffer and 1X maleate buffer prior overnight incubation (4°C) with anti-digoxigenin (FAB fragments) conjugated to alkaline phosphatase antibody (Roche Diagnostic) in antibody buffer (1:2000). Fast red staining was performed by incubation with 0.05 mg naphtol substrate, 0.2 mg Fast Red chromogen, and 0.04 mg levamisole in 100 mM Tris-HCl pH 8.2 and 0.1% Tween-20 overnight at 4°C in darkness, according to manufacturer’s protocol (Fast Red tablets—Roche Diagnostic). Sections were then washed in 1X PBS, mounted with VECTASHIELD Mounting Medium (Vector laboratories) and analyzed with Leica TCS laser scanning confocal microscope using 10 and 20X oil immersion objectives. Image stacks of 1 μm optical sections were acquired with Leica PowerScan 2.5 software followed by brightness and contrast adjustments using Adobe Photoshop CS6 (San Jose, CA) before being mounted on panels.

Ethics Statement

Austevoll Research station has the following permission for catch and maintenance of Atlantic cod: H-AV 77, H-AV 78 and H-AV 79, all issued by the Norwegian Directorate of Fisheries. Austevoll Research station has a general permission to run as a Research Animal facility and to conduct experiments involving all developmental stages of fish (code 93) provided by the national Animal Care and Use Committee (ACUC, the Norwegian Animal Research Authority, FDU, www.fdu.no). Although the sampling of stages of larval cod used in the current study do not require a specific permit according to Norwegian law of Animals in Research, Regulation of the 15^th of January 1996; the method of euthanization of fish larvae with the use of MS222 (metacaine) follows § 16 Euthanasia of Laboratory Animals stating that the selected euthanasia should involve no signs of unnecessary suffering of the animal.

This study is not categorized as a field study. The pond was merely used to filter different stages of mainly calanoid copepods which were used as feed for the fish larvae. The owner of the land has given permission to conduct such filtering. Calanoid copepods are not endangered or protected species and are not recognized as research animals in the Regulation on Animal Experimentation of 15. January 1996, last changed 6. August 2010, and no specific permission is therefore needed in order to use copepods as feed. The specific experiment described in this study was approved by the local representative of the Norwegian Research Animal Committee (Forsøksdyrutvalget, FDU, www.fdu.no) with the acceptance reference FOTS id 5448.

Results

To increase our understanding of the gene distribution and putative roles of neuropeptides NPY, CART and pOX in selected brain regions of Atlantic cod larvae, we examined their mRNA expression patterns in serial brain sections by using both chromogenic and dual fluorescent in situ hybridization methods.

Identification of brain regions and nuclei were mostly based on Atlas of Early Zebrafish Brain Development [24] and available literatures on similar studies in other teleost species [17, 25–29]. NPY, CART and pOX mRNA expression patterns in the forebrain and midbrain are summarized in Table 1 and Figs 1–4.

[Figure omitted. See PDF.]

Table 1. Distribution of NPY, CART, and pOX mRNA brain expression by ISH methods in Atlantic cod at different larval developmental stages.

https://doi.org/10.1371/journal.pone.0153743.t001

[Figure omitted. See PDF.]

Fig 1. Schematic representation of rostro-caudal series of transverse sections of Atlantic cod brain at 4 days post hatch (dph); Stage 1.

mRNA expression of NPY (purple dots), CART (green dots), pOX (blue dots), mRNA co-expression of NPY and CART (black dots), and pOX and CART (orange dots). Larval age (days post hatch; dph) correlates to Stage (based on Ø. Sæle, NIFES, Bergen, Norway, personal communication). Drawings in upper left corners show relative position of brain transverse sections.

https://doi.org/10.1371/journal.pone.0153743.g001

[Figure omitted. See PDF.]

Fig 2. Schematic representation of rostro-caudal series of transverse sections of Atlantic cod brain at 20 days post hatch (dph); Stage 2.

mRNA expression of NPY (purple dots), CART (green dots), pOX (blue dots), mRNA co-expression of NPY and CART (black dots), and pOX and CART (orange dots). Other captions as in Fig 1.

https://doi.org/10.1371/journal.pone.0153743.g002

[Figure omitted. See PDF.]

Fig 3. Schematic representation of rostro-caudal series of transverse sections of Atlantic cod brain at 46 days post hatch (dph); Stage 4.

mRNA expression of NPY (purple dots), CART (green dots), pOX (blue dots), mRNA co-expression of NPY and CART (black dots), and pOX and CART (orange dots). Other captions as in Fig 1.

https://doi.org/10.1371/journal.pone.0153743.g003

[Figure omitted. See PDF.]

Fig 4. Schematic representation of rostro-caudal series of transverse sections of Atlantic cod brain at 76 days post hatch (dph); Stage Juvenile.

mRNA expression of NPY (purple dots), CART (green dots), pOX (blue dots), mRNA co-expression of NPY and CART (black dots), and pOX and CART (orange dots). Other captions as in Fig 1.

https://doi.org/10.1371/journal.pone.0153743.g004

NPY

NPY transcripts were detected in the telencephalon, diencephalon and optic tectum from 4 dph and its distribution pattern varied between all stages.

NPY mRNA in the telencephalon.

At 4 dph larvae, most of the NPY mRNA expressing neurons were localized within the ventral telencephalon (VTe, Fig 5A), whereas at 20, 46 and 76 dph, two distinct symmetric clusters developed in the lateral part of the dorsal telencephalon (DlTe, Fig 5D, 5G and 5K). The telencephalon was the brain region showing the highest NPY mRNA- expression at all stages investigated.

[Figure omitted. See PDF.]

Fig 5. Localization of NPY gene expression in serial sections of brain in Atlantic cod larvae.

At 4dph, (A) NPY mRNA expression in the ventral telencephalon (VTe), (B) several cells contained NPY gene expression in the ventral layer of the optic tectum (OT) and in the preoptic region (Po), (C) scattered cells expressing NPY mRNAs in the OT. At 20dph, (D) NPY mRNA expression in the lateral of dorsal telencephalon (DlTe) and in the VTe, (E) symmetric clusters of NPY expressing cells in the dorsal thalamus (DTh) and in the ventral layer of the OT, (F) hybridization signal in the ventral layer of the OT and in the hypothalamus (Hy). At 46dph, (G) symmetric clusters of NPY- expressing cells in the DlTe, (H) NPY in the periventricular gray zone (PGZ) and in the lateral tuberal nucleus of hypothalamus (NLT), (I) NPY in the PGZ of the caudal OT. At 76dph, (K) hybridization signal in the DlTe and VTe, (L) NPY- symmetric expression in the preoptic region (Po), (M) NPY symmetric expression in the NLT, (N) NPY symmetric expression in the PGZ of the OT, in the periventricular pretectal nucleus (PPN) and in the NLT.

https://doi.org/10.1371/journal.pone.0153743.g005

NPY mRNA in the diencephalon and mesencephalon.

In the preoptic region, NPY transcripts were found only in 4 and 76 dph old larvae (Po, Fig 5B and 5L). At 20, 46 and 76 dph, NPY mRNA containing neurons were also detected dorsally and laterally to the pituitary gland, in the presumptive lateral tuberal nucleus of the hypothalamus (Hy, Fig 5F; NLT, Fig 5H, 5M and 5N). NPY expression in the dorsal thalamus was only observed at the 20 dph stage (DTh, Fig 5E and 5F). At 76 dph, two symmetric clusters of NPY expressing cells were also observed within the presumptive periventricular pretectal nucleus (PPN, Fig 5N). The optic tectum exhibited NPY mRNA expression throughout its rostro-caudal brain extent, in its ventral layer at 4 and 20 dph (OT, Fig 5B, 5C, 5E and 5F) and in the periventricular gray zone at older stages (PGZ, Fig 5H, 5I and 5N).

CART

CART mRNA showed a wide temporal distribution in the brain of the larvae, being expressed in both the telencephalon and diencephalon from 4 dph onwards. At 46 and 76 dph, new populations of CART positive neurons were also identified in the olfactory bulb, in the region of the nucleus of medial longitudinal fascicle (N), the periventricular gray zone of the optic tectum (PGZ) and in the midbrain tegmentum.

CART mRNA in the telencephalon.

In the olfactory bulbs, two distinct populations of CART mRNA cells were visible at 46 dph larvae (OB, Fig 6H). The number of CART mRNA labeled cells greatly increased in two larger and rounded domains at 76 dph stage (OB, Fig 6M). Several scattered cells expressing CART were detected in the presumptive lateral part of the dorsal telencephalon (DlTe), at all stages investigated (Fig 6A, 6E, 6I and 6N).

[Figure omitted. See PDF.]

Fig 6. Localization of CART gene expression in serial sections of brain in Atlantic cod larvae.

At 4dph, (A) CART positive cells in the lateral parts of dorsal telencephalon (DlTe), (B) CART mRNAs in the dorsal thalamus (DTh) and the preoptic region (Po), (C) symmetrical clusters of CART-containing neurons in the posterior tuberculum (PTc), (D) hybridization signal of CART in PTc and in the hypothalamus (Hy). At 20dph, (E) CART mRNAs in the DlTe, (F) CART expression in the rostral hypothalamus (RHy), (G) CART in the region of the nucleus of medulla longitudinal fascicle (N) and in the PTc. At 46dph, (H) two distinct clusters of CART mRNA- expressing cells in the olfactory bulbs (OB), (I) clusters of labeled neurons in the DlTe and in the Po, K: CART gene expression in the Hy, PTc and in the periventricular gray zone of the optic tectum (PGZ), (L) clusters of CART expressed cells in the region of the nucleus of medulla longitudinal fascicle (N). At 76dph, (M, N, O) show: CART mRNA- cells in the OB, in the DlTe and in the Po, respectively. (P) CART mRNA in the dorsal thalamus (DTh), Hy and in the diffuse nucleus of the inferior hypothalamic lobe (DIL). (Q) CART containing cells distribute in the PGZ, N and new clusters of labeled neurons in the putative lateral tuberal nucleus (NLT). (R) CART in the midbrain tegmentum (MTeg).

https://doi.org/10.1371/journal.pone.0153743.g006

CART mRNA in the diencephalon and mesencephalon.

CART mRNAs had a wide distribution in the diencephalon, particularly in the preoptic region at all stages investigated (Po, Fig 6B, 6I and 6O for 4, 46, 76 dph, respectively and Fig 7A₃ for 20 dph) and in the caudal and rostral hypothalamus at 4 and 20 dph, respectively (Hy, Fig 6D; RHy, Fig 6F). At 76 dph, CART- expressing neurons were found in the presumptive lateral tuberal nucleus (NLT, Fig 6Q) and in the diffuse nucleus of the inferior hypothalamic lobe (DIL, Fig 6P and 6Q). CART mRNAs- were also widely spread within the thalamus and in the posterior tuberculum. In the 4 dph larvae, a weak expression was seen in the dorsal thalamus (DTh, Fig 6B) while distinct symmetric clusters of CART positive neurons were detected in the presumptive lateral parts of the ventral posterior tuberculum (PTc, Fig 6C and 6D). From 20 dph onward, new populations of CART mRNA-cells were found in the presumptive nucleus of the medial longitudinal fascicle (N, Fig 6G, 6L and 6Q). At both 46 and 76 dph CART mRNAs were found in the presumptive periventricular gray zone of the optic tectum (PGZ, Fig 6K, 6L and 6Q) and in the midbrain tegmentum (MTeg, Fig 6R).

[Figure omitted. See PDF.]

Fig 7. Dual FISH photomicrographs showing localization of pOX and CART gene expression in brain of Atlantic cod larvae at 20 and 76 dph.

The left, middle, and right columns show fluorescent signals corresponding to red, green and overlap channels as reported above. (A_1-3), pOX and CART mRNAs in the preoptic region (Po) at 20 dph; (B_1-3), and (C_1-3), in the rostral (RHy) and intermediate hypothalamus respectively, at 76 dph.

https://doi.org/10.1371/journal.pone.0153743.g007

pOX

Prepro-orexin (pOX) had a restricted distribution in the cod brain throughout the four larval stages studied. At 4 and 20 dph, pOX mRNA were distributed in two clearly defined symmetric neuron clusters in the preoptic region (Po, Fig 8A and 8B). Later, at both 46 and 76 dph, these cell populations had a more spread distribution from rostral to intermediate regions of the hypothalamus (RHy and Hy, Fig 8C–8F).

[Figure omitted. See PDF.]

Fig 8. Localization of pOX gene expression in serial sections of brain in Atlantic cod larvae.

ISH signals were detected in the preoptic region (Po) at 4 dph (A) and 20 dph (B) and in the rostral (RHy) and intermediate (Hy) hypothalamus at 46 dph (C, D), and 76 dph (E, F).

https://doi.org/10.1371/journal.pone.0153743.g008

Co-localization of CART/pOX and CART/NPY in the brain of Atlantic cod larvae

On the basis of single-gene detection studies presented here, dual fluorescent ISH protocol was additionally used to detect potential mRNA co-localization of CART/pOX and CART/NPY in same brain sections and at the level of single neurons. Dual fluorescent ISH of CART/ pOX genes identified adjacent CART and pOX expressing cells in the preoptic region at 20 dph (Po, Fig 7A_1-3). Interestingly, at 76 dph, CART and pOX mRNA-expressing cells shared cellular domains in the intermediate hypothalamus where some cells clearly co-expressed both genes (Fig 7B_1-3 and 7C_1-3). At 46 dph, both CART and NPY expressing cells were detected in the lateral parts of the dorsal telencephalon (DlTe, Fig 9A_1-3) and in the periventricular gray zone of the optic tectum (PGZ, Fig 9B_1-3). Particularly, mRNA co-localization was seen in few neurons of the DlTe (Fig 9A₃).

[Figure omitted. See PDF.]

Fig 9. Dual FISH photomicrographs showing localization of NPY and CART gene expression in brain of Atlantic cod larvae at 46 dph.

Fluorescence signals were recorded in the Fast Red (NPY, left column; detection of wavelengths greater than 560 nm) and TSA (CART, middle column; detection of wavelengths from 500 nm to 540 nm) detection channels. Right column shows both channels overlap to visualize co-localization of NPY and CART gene expression (yellow). (A_1-3), NPY and CART mRNAs in the lateral parts of the dorsal telencephalon (DlTe); (B_1-3), in the periventricular gray zone of the optic tectum (PGZ).

https://doi.org/10.1371/journal.pone.0153743.g009

Discussion

The present study describes for the first time the development expression patterns of NPY, CART and pOX genes in discrete brain regions of the marine teleost Atlantic cod, from start of exogenous feeding until juvenile stage. The mRNA distribution indicates that all three brain regions investigated (telencephalon, diencephalon and mesencephalon) contain both orexigenic and anorexigenic key elements in appetite control. The spatial and temporal expression patterns indicate a progressive development of the brain regulatory networks that control appetite during larval ontogeny.

NPY

The NPY mRNA expression changed spatially during larval development. At the first feeding stage (4 dph), high NPY mRNA expression levels were found in the ventral telencephalon, diencephalon, and optic tectum. From Stage 2 (20 dph) to the juvenile stage (76 dph), two distinct and growing symmetric clusters of labeled cells appeared in the lateral part of ventral telencephalon. From 20 dph onwards NPY- gene expressing cells were also seen in the ventrolateral parts of hypothalamus. In the mesencephalon, NPY mRNA was found throughout the rostro-caudal extent of the optic tectum in the presumptive periventricular gray zone at 4 and 20 dph and in the same developed region at older stages.

NPY RT-PCR expression during early developmental stages of fish has also been previously reported; for instance at 3 dph in Atlantic cod larvae [12], at blastula stage in orange spotted grouper [19] and at zygote stage in blunt snout bream (Megalobrama amblycephala) [30, 31]. Moreover, immunoreactivity-studies in zebrafish have shown that NPY-ir cell bodies appeared for the first time at 60 h post fertilization (hpf) in the hypothalamus [26]. NPY-ir neurons were also found in the olfactory placode at 16 dpf in chum salmon (Oncorhynchus keta) [31], and in 40 mm-staged embryos of cloudy dogfish (Scyliorhinus torazame) [31]. Thus, it is likely that ontogeny of the NPY system, both at mRNA and protein levels, is required for a proper brain development at very early embryonic stages.

The telencephalon showed the highest NPY expression at all stages investigated, in accordance with previous studies in juvenile Atlantic cod [32] and other fish species at adult stages (e.g. sea bass Dicentrarchus labrax [17], goldfish (Carassius auratus) [18, 28, 33, 34], winter skate (Raja ocellata) [35], killifish (Fundulus heteroclitus) [36], African lungfish (Protopterus annectens) [37, 38], and white sturgeon (Acipenser transmontanus) [39, 40]. Due to the proposed role of NPY as an appetite stimulator in fish [11, 41, 42], this could indicate that the telencephalon may possess some components in hunger regulation in Atlantic cod and that this system is already present from onset of exogenous feeding.

The expression of NPY mRNAs in the diencephalon is in agreement with previous studies that report NPY mRNA in the preoptic region from the early larval stages, juvenile and adult fish (e.g. goldfish [18, 34], sea bass [17], killifish [36], African lungfish [37], and gilthead seabream (Sparus aurata) [43]). ISH signals were mainly observed in the ventro-lateral parts of hypothalamus, referred to as the putative NLT at the two older stages investigated. In previous studies, NPY-ir and/or mRNA transcript detection in NLT were observed in sea bass [17], dogfish (Scyliorhinus canicula) [44], white sturgeon [39], killifish [36], and pejerrey (Odontesthes bonariensis) [45].

In goldfish, NPY mRNA is widely distributed in the thalamus [18]. In contrast, we found NPY expression only in the dorsal thalamus at 20 dph stage. Given the organization of the pretectum in cichlid fish [27, 46], the two clusters of NPY mRNA expressing cells identified in the dorsal midbrain of cod larvae (Fig 5N) may correspond to the periventricular pretectal nucleus (PPN), particularly to the ventral periventricular pretectal nucleus (PPv). Similar observations have also been reported in adult sea bass [17].

In agreement with our results in Atlantic cod larvae, NPY mRNA expression has been observed within the optic tectum of goldfish [18, 33], sea bass [17] and winter skate [47], and numerous studies have reported that NPY- ir cells and/or fibers appear in the optic tectum in goldfish [28, 34], in Atlantic salmon [48, 49], in pejerrey [45], in gilthead seabream [43] and in lamprey Petromyzon marinus [50].

CART

In the present study, CART mRNA expression was detected in the telencephalon, preoptic region and diencephalon from 4 dph onwards. After metamorphosis, at 46 and 76 dph, CART was also found in olfactory bulb and optic tectum. The early mRNA expression found in the telencephalon and its late appearance in the olfactory bulbs is in agreement with a previous study on zebrafish where CART immunoreactivity (ir) cells in the telencephalon were detected as early as 24 hpf, while in the olfactory bulbs they first appeared at 46 hpf [25]. We found the highest level of CART mRNA in the olfactory bulbs at the latest larval stage investigated, also in agreement with previous observations in goldfish [21]. Taken together, these results suggest that the CART peptide may be associated with the regulation of feeding in the early larval stages and might have additional roles in sensory (olfactory and visual) information/transmission at later stages.

The expression of CART mRNA in the diencephalon showed interesting spatio-temporal patterns during cod larval development. CART mRNAs were initially found in the preoptic region as early as 20 dph, with expression gradually increasing at later stages [25]. The preoptic region is known to play an important role in both reproduction and feeding, in fish [51, 52] as in mammals [53]. In the sexually inmature cod larvae the early onset of CART mRNA expression in the preoptic region may therefore suggest that this region aquire an involvement in feeding behavior from the early larval stages.

CART mRNA-expressing cells were detected in rostral and intermediate hypothalamus from 4 dph onwards. The hypothalamic domains expressing CART mRNA at 76 dph were identified the presumptive NLT and the DIL. Similarly, in zebrafish CART(ir) fibers in the hypothalamus were first detected at 4 dpf in rostral, intermedial, and central zones of the periventricular hypothalamus, DIL, and lateral hypothalamus [25].

The high expression of CART, combined with the absence of the orexigenic modulator NPY in the DIL, suggest that this nucleus may serves as a “satiety center” in Atlantic cod. DIL may be the functional homologue to the mammalian ventromedial nucleus (VMN) situated in the intermediate hypothalamus and associated with satiety control [54, 55].

We also observed a wide distribution of CART mRNAs in the thalamus and in the posterior tuberculum (PTc) especially from 46 dph, and new groups of labeled cells appeared in the putative (N) from 46 dph onwards. The presence of CART mRNA expression in the thalamus of cod larvae is in agreement with previous studies on CART-(ir) cell distribution in other teleost fish [25, 56]. In fish, thalamic nuclei are involved in relaying sensory inputs to the telencephalon [57], and the wide distribution of CART mRNAs we found in the same brain regions reinforces the idea that this peptide may be related to sensory information processing in cod larvae as well.

CART mRNA-expressing cells were also detected in the optic tectum and in the midbrain tegmentum, from 46 dph onwards. The teleostean optic tectum is the primary target of the optic nerve fibers [57]. Its involvement in processing of visual information has been shown well documented in teleost species [29, 58]. The detection of CART mRNAs in the ganglion cell layer (GCL) in the present study (data not shown) could indicate that the CART peptide is involved in visual processing.

In summary, CART mRNA expression was apparent at the start of feeding and the occurrence of a wide and dense CART expression in thalamus and hypothalamus suggests that neuronal cell populations identifying key components controlling satiety may be located in these areas. Moreover, the presence of CART mRNA in the olfactory bulb and optic tectum suggests that CART functions outside that of appetite control, perhaps modulating sensory systems.

pOX

The pOX mRNA expression was exclusively located in the preoptic region at 4 and 20 dph, and extended to the hypothalamus from 46 dph onwards. The hypothalamic pOX expression in our study is consistent with other ISH studies in zebrafish at both embryonic to larval stages [59] and in adults [20, 60, 61]. However using RT-PCR methods in adult Atlantic cod, pOX mRNAs were detected in several other brain regions including the olfactory bulb, telencephalon, hypothalamus and optic tectum [16]. Similar RT-PCR expression profiles were also found in mature specimens of winter skate [62], winter flounder [63], orange-spotted grouper [15] and cichlids [20]. In adult fish, the pOX peptide has been demonstrated using immunohistochemistry in the ventral telencephalon, thalamus, and hypothalamus of goldfish [64], extensively throughout the brain in medaka [61], and in the preoptic region and hypothalamus in cichlids [20]. Together, these results indicate that in adult fish pOX transcripts are detected in several brain regions while at early embryonic and larval stages pOX expression/functions are primarily exerted through neurons located in the preoptic area and hypothalamus. As environmental challenges and related feeding behaviors become more complex in adult fish, it is reasonable that the mechanisms underlying these brain functions may require more sophisticated neurochemical and anatomical brain networks.

Co-localization of CART/NPY and CART/pOX in the brain of Atlantic cod larvae

CISH for single-gene detection indicated some potential sites of co-localization of mRNAs for CART and NPY or pOX genes. DFISH detection protocol was implemented using combination pairs of the anorexigenic factor CART with either NPY or pOX orexigenic genes in the same brain sections. The distribution of both NPY and CART mRNAs by DFISH in the telencephalon and periventricular gray zone of the optic tectum were consistent with results recorded by CISH. Numerous NPY and CART mRNA positive neurons were detected in same brain regions, e.g. telencephalon and optic tectum, however only few neurons co-expressing these two genes were found in the DlTe.

Both NPY and CART mRNA transcripts were found in the lateral tuberal nucleus (NLT) at 76 dph juveniles. In mammals, the hypothalamic nucleus that releases both orexigenic and anorexigenic peptides is the arcuate nucleus [65], suggesting that the NLT might be the teleostean homolog of the mammalian arcuate nucleus (ARC) [66, 67].

Interestingly, we also found clear co-localization of CART and pOX mRNAs in hypothalamic neurons of the caudal part of pOX cell population, which overlaps with the rostral part of CART cell population. The co-expression of CART and NPY in the same neurons supports the idea that this hypothalamic area exerts a multifunctional control of appetite.

Overall, CISH and DFISH data show a wide distribution of NPY, pOX, and CART expression in the hypothalamus of cod larvae substantiating its central role in control of food intake as previously showed in other teleost fish [17,45] and mammals [9–11].

Moreover, the temporal and spatial differences in NPY, CART, and pOX mRNA expression from first feeding stage (4 dph), through metamorphosis and until juvenile (76 dph) stages indicate a dynamic distribution consistent with progression of brain development and function. These observations support the statement that the neurotransmitter systems underlying appetite control continues to develop during the larval stages of teleosts as well as after the onset of exogenous feeding [68, 69].

Acknowledgments

We thank Valentina Tronci at Uni Research for support in in situ hybridization method, Tom Ole Nilsen for designing in situ hybridization primers for cod NPY and CART. We also thank Ann-Elise Olderbakk Jordal for equipment support, Jon Vidar Helvik and Martina Raasholm for microscope use guide.

Author Contributions

Conceived and designed the experiments: IR ARA LOEE. Performed the experiments: ØK. Analyzed the data: ARA HTMDL LOEE. Wrote the paper: HTMDL ARA LOEE ØK IR.

Citation: Le HTMD, Angotzi AR, Ebbesson LOE, Karlsen Ø, Rønnestad I (2016) The Ontogeny and Brain Distribution Dynamics of the Appetite Regulators NPY, CART and pOX in Larval Atlantic Cod (Gadus morhua L.). PLoS ONE 11(4): e0153743. https://doi.org/10.1371/journal.pone.0153743

References

1. Øiestad V, Kvenseth PG, Folkvord A. Mass production of Atlantic cod juveniles Gadus morhua in a Norwegian saltwater pond. Trans Am Fish Soc. 1985 1985/07/01;114(4):590–5.

2. FAO. Food and Agriculture Organization of the United Nations Cultured aquatic species information programme Gadus morhua. Cultured aquatic species information programme. Text by Håkon Otterå. In: FAO Fisheries and Aquaculture Department [online]. Rome. Updated 1 January 2004. 2004 [5 February 2014].

3. Johansen SD, Coucheron DH, Andreassen M, Karlsen BO, Furmanek T, Jørgensen TE, et al. Large-scale sequence analyses of Atlantic cod. New Biotechnol. 2009 06/;25(5):263–71.

4. Rønnestad I. Norway Research Council invest NOK 21 million into research on cod larvae. Int Aquafeed Mag. 2013;May June:3.

5. Auditore PJ, Lough RG, Broughton EA. A review of the comparative development of Atlantic cod (Gadus morhua L.) and haddock (Melanogrammus aeglefinus L.) based on an illustrated series of larvae and juveniles from Georges Bank. NAFO Science Council Studies, pp. 7–18., Northwest Atlantic Fisheries Organisation. 1994.

6. Kamisaka Y, Rønnestad I. Reconstructed 3D models of digestive organs of developing Atlantic cod (Gadus morhua) larvae. Mar Biol. 2011 2011/01/01;158(1):233–43. pmid:24489408

7. Kvåle A, Mangor-Jensen A, Moren M, Espe M, Hamre K. Development and characterisation of some intestinal enzymes in Atlantic cod (Gadus morhua L.) and Atlantic halibut (Hippoglossus hippoglossus L.) larvae. Aquaculture. 2007 4/6/;264(1–4):457–68.

8. Rønnestad I, Yúfera M, Ueberschär B, Ribeiro L, Sæle Ø, Boglione C. Feeding behaviour and digestive physiology in larval fish: current knowledge, and gaps and bottlenecks in research. Reviews in Aquaculture. 2013;5:S59–S98.

9. Dietrich MO, Horvath TL. Feeding signals and brain circuitry. Eur J Neurosci. 2009;30(9):1688–96. pmid:19878280

10. Williams KW, Elmquist JK. Lighting up the hypothalamus: Coordinated control of feeding behavior. Nat Neurosci. 2011 Mar;14(3):277–8. pmid:21346745

11. Volkoff H, Canosa LF, Unniappan S, Cerdá-Reverter JM, Bernier NJ, Kelly SP, et al. Neuropeptides and the control of food intake in fish. Gen Comp Endocrinol. 2005;142(1–2):3–19. pmid:15862543

12. Kortner TM, Overrein I, Øie G, Kjørsvik E, Bardal T, Wold P-A, et al. Molecular ontogenesis of digestive capability and associated endocrine control in Atlantic cod (Gadus morhua) larvae. Comp Biochem Physiol A Physiol. 2011 10//;160(2):190–9.

13. Hunter RG, Philpot K, Vicentic A, Dominguez G, Hubert GW, Kuhar MJ. CART in feeding and obesity. Trends Endocrinol Metab. 2004 11//;15(9):454–9. pmid:15519893

14. Sakurai T. Roles of orexins in the regulation of feeding and arousal. Sleep Med. 2002 12//;3, Supplement 2(0):S3–S9.

15. Yan A, Zhang L, Tang Z, Zhang Y, Qin C, Li B, et al. Orange-spotted grouper (Epinephelus coioides) orexin: Molecular cloning, tissue expression, ontogeny, daily rhythm and regulation of NPY gene expression. Peptides. 2011 7//;32(7):1363–70. pmid:21600944

16. Xu M, Volkoff H. Molecular characterization of prepro-orexin in Atlantic cod (Gadus morhua): Cloning, localization, developmental profile and role in food intake regulation. Mol Cell Endocrinol. 2007 6/15/;271(1–2):28–37. pmid:17434256

17. Cerdá-Reverter JM, Anglade I, Martínez-Rodríguez G, Mazurais D, Muñoz-Cueto JA, Carrillo M, et al. Characterization of neuropeptide Y expression in the brain of a perciform fish, the sea bass (Dicentrarchus labrax). J Chem Neuroanat. 2000 9//;19(4):197–210. pmid:11036237

18. Peng C, Gallin W, Peter RE, Blomqvist AG, Larhammar D. Neuropeptide-Y gene expression in the goldfish brain: Distribution and regulation by ovarian steroids. Endocrinology. 1994 Mar;134(3):1095–103. pmid:8119148

19. Chen R, Li W, Lin H. cDNA cloning and mRNA expression of neuropeptide Y in orange spotted grouper, Epinephelus coioides. Comp Biochem Physiol B Biochem Mol Biol. 2005 9//;142(1):79–89. pmid:16005251

20. Pérez Sirkin DI, Suzuki H, Cánepa MM, Vissio PG. Orexin and neuropeptide Y: Tissue specific expression and immunoreactivity in the hypothalamus and preoptic area of the cichlid fish Cichlasoma dimerus. Tissue Cell. 2013 12//;45(6):452–9. pmid:24138942

21. Volkoff H, Peter RE. Characterization of two forms of cocaine- and amphetamine-regulated transcript (CART) peptide precursors in goldfish: Molecular cloning and distribution, modulation of expression by nutritional status, and interactions with leptin. Endocrinology. 2001 Dec;142(12):5076–88. pmid:11713200

22. Karlsen Ø, van der Meeren T, Rønnestad I, Mangor-Jensen A, Galloway TF, Kjørsvik E, et al. Copepods enhance nutritional status, growth, and development in Atlantic cod (Gadus morhua L.) larvae—can we identify the underlying factors? PeerJ. 2015;3:e902. pmid:26038712

23. Sandbakken M, Ebbesson L, Stefansson S, Helvik JV. Isolation and characterization of melanopsin photoreceptors of atlantic salmon (Salmo salar). J Comp Neuro. 2012;520(16):3727–44.

24. Mueller T, Wullimann MF. Atlas of early zebrafish brain development: A tool for molecular neurogenetics: Elsevier; 2005.

25. Mukherjee A, Subhedar NK, Ghose A. Ontogeny of the cocaine- and amphetamine-regulated transcript (CART) neuropeptide system in the brain of zebrafish, Danio rerio. J Comp Neuro. 2012;520(4):770–97.

26. Mathieu M, Tagliafierro G, Bruzzone F, Vallarino M. Neuropeptide tyrosine-like immunoreactive system in the brain, olfactory organ and retina of the zebrafish, Danio rerio, during development. Dev Brain Res. 2002 12/15/;139(2):255–65.

27. Huffman LS, O’Connell LA, Kenkel CD, Kline RJ, Khan IA, Hofmann HA. Distribution of nonapeptide systems in the forebrain of an African cichlid fish, Astatotilapia burtoni. J Chem Neuroanat. 2012 7//;44(2):86–97. pmid:22668656

28. Kah O, Pontet A, Danger JM, Dubourg P, Pelletier G, Vaudry H, et al. Characterization, cerebral distribution and gonadotropin release activity of neuropeptide Y (NPY) in the goldfish. Fish Physiol Biochem. 1989 1989/06/01;7(1–6):69–76. pmid:24221756

29. Castro A, Becerra M, Manso MJ, Anadón R. Development of immunoreactivity to neuropeptide Y in the brain of brown trout (Salmo trutta fario). J Comp Neuro. 1999;414(1):13–32.

30. Ping HC, Feng K, Zhang GR, Wei KJ, Zou GW, Wang WM. Ontogeny expression of ghrelin, neuropeptide Y and cholecystokinin in blunt snout bream, Megalobrama amblycephala. J Anim Physiol Anim Nutr (Berl). 2013 May 21.

31. Chiba A, Oka S, Honma Y. Ontogenetic changes in neuropeptide Y-like-immunoreactivity in the terminal nerve of the chum salmon and the cloudy dogfish, with special reference to colocalization with gonadotropin-releasing hormone-immunoreactivity. Neurosci Lett. 1996 07/;213(1):49–52. pmid:8844710

32. Kehoe AS, Volkoff H. Cloning and characterization of neuropeptide Y (NPY) and cocaine and amphetamine regulated transcript (CART) in Atlantic cod (Gadus morhua). Comp Biochem Physiol A Physiol. 2007 3//;146(3):451–61.

33. Vecino E, Perez MT, Ekstrom P. In situ hybridization of neuropeptide Y (NPY) mRNA in the goldfish brain. Neuroreport. 1994 Dec 30;6(1):127–31. pmid:7703400

34. Pontet A, Danger JM, Dubourg P, Pelletier G, Vaudry H, Calas A, et al. Distribution and characterization of neuropeptide Y-like immunoreactivity in the brain and pituitary of the goldfish. Cell Tissue Res. 1989 1989/03/01;255(3):529–38. pmid:23885381

35. MacDonald E, Volkoff H. Neuropeptide Y (NPY), cocaine- and amphetamine-regulated transcript (CART) and cholecystokinin (CCK) in winter skate (Raja ocellata): cDNA cloning, tissue distribution and mRNA expression responses to fasting. Gen Comp Endocrinol. 2009 4//;161(2):252–61. pmid:19523382

36. Subhedar N, Cerda J, Wallace RA. Neuropeptide Y in the forebrain and retina of the killifish, Fundulus heteroclitus. Cell Tissue Res. 1996 Feb;283(2):313–23. pmid:8593660

37. Vallarino M, Tranchand-Bunel D, Thoumas JL, Masini MA, Conlon JM, Fournier A, et al. Neuropeptide tyrosine in the brain of the African lungfish, Protopterus annectens: Immunohistochemical localization and biochemical characterization. J Comp Neuro. 1995 Jun 12;356(4):537–51.

38. Vallarino M, Angela Masini M, Trabucchi M, Mathieu M, Vaudry H. Autoradiographic distribution of neuropeptide tyrosine binding sites in the brain of the African lungfish, Protopterus annectens. Neurosci Lett. 1998 9/18/;254(1):5–8. pmid:9780078

39. Chiba A, Honma Y. Neuropeptide Y-immunoreactive structures in the telencephalon and diencephalon of the white sturgeon, Acipenser transmontanus, with special regard to the hypothalamo-hypophyseal system. Arch Histol Cytol. 1994 Mar;57(1):77–86. pmid:8198837

40. Peter RE. The brain and feeding behavior. In: Hoar W.S R DJ, Brett JR, editors. Fish Physiology: Academic Press; 1979. p. 121–59.

41. Lin X, Volkoff H, Narnaware Y, Bernier NJ, Peyon P, Peter RE. Brain regulation of feeding behavior and food intake in fish. Comp Biochem Physiol A Physiol. 2000 8/1/;126(4):415–34.

42. Volkoff H. The role of neuropeptide Y, orexins, cocaine and amphetamine-related transcript, cholecystokinin, amylin and leptin in the regulation of feeding in fish. Comp Biochem Physiol A Physiol. 2006 Jul;144(3):325–31.

43. Pirone A, Lenzi C, Marroni P, Betti L, Mascia G, Giannaccini G, et al. Neuropeptide Y in the brain and retina of the adult teleost gilthead seabream (Sparus aurata L.). Anat Histol Embryol. 2008;37(3):231–40. pmid:18205886

44. Vallarino M, Danger JM, Fasolo A, Pelletier G, Saint-Pierre S, Vaudry H. Distribution and characterization of neuropeptide Y in the brain of an elasmobranch fish. Brain Res. 1988 5/10/;448(1):67–76. pmid:3390717

45. Traverso JM, Ravaglia MA, Vissio PG, Maggese MC, Paz DA. Localization of neuropeptide y-like immunoreactive structures in the brain of the pejerrey, Odontesthes bonariensis (Teleostei, Atheriniformes). Anat Histol Embryol. 2003;32(1):29–35. pmid:12733270

46. Fernald RD, Shelton LC. The organization of the diencephalon and the pretectum in the cichlid fish, Haplochromis burtoni. J Comp Neuro. 1985;238(2):202–17.

47. MacDonald E, Volkoff H. Cloning, distribution and effects of season and nutritional status on the expression of neuropeptide Y (NPY), cocaine and amphetamine regulated transcript (CART) and cholecystokinin (CCK) in winter flounder (Pseudopleuronectes americanus). Horm Behav. 2009;56(1):58–65. pmid:19303880

48. Vecino E, Ekström P. Distribution of met-enkephalin, leu-enkephalin, substance P, neuropeptide Y, FMRFamide, and serotonin immunoreactivities in the optic tectum of the Atlantic salmon (Salmo salar L.). J Comp Neuro. 1990;299(2):229–41.

49. Silverstein JT, Breininger J, Baskin DG, Plisetskaya EM. Neuropeptide Y-like gene expression in the salmon brain increases with fasting. Gen Comp Endocrinol. 1998 5//;110(2):157–65. pmid:9570936

50. Barreiro-Iglesias A, Anadón R, Rodicio MC. New insights on the neuropeptide Y system in the larval lamprey brain: Neuropeptide Y immunoreactive neurons, descending spinal projections and comparison with tyrosine hydroxylase and GABA immunoreactivities. Neuroscience. 2010;167(2):396–413. pmid:20167263

51. Bruin JC. Telencephalon and behavior in teleost fish. In: Ebbesson SE, editor. Comparative Neurology of the Telencephalon: Springer US; 1980. p. 175–201.

52. Kyle AL, Sorensen PW, Stacey NE, Dulka JG. Medial olfactory tract pathways controlling sexual reflexes and behavior in teleosts. Ann N Y Acad Sci. 1987;519(1):97–107.

53. Roos J, Roos M, Schaeffer C, Aron C. Sexual differences in the development of accessory olfactory bulbs in the rat. J Comp Neuro. 1988;270(1):121–31.

54. King BM. The rise, fall, and resurrection of the ventromedial hypothalamus in the regulation of feeding behavior and body weight. Physiol Behav. 2006 2/28/;87(2):221–44. pmid:16412483

55. Williams G, Bing C, Cai XJ, Harrold JA, King PJ, Liu XH. The hypothalamus and the control of energy homeostasis: Different circuits, different purposes. Physiol Behav. 2001 11/12/;74(4–5):683–701. pmid:11790431

56. Singru PS, Mazumdar M, Sakharkar AJ, Lechan RM, Thim L, Clausen JT, et al. Immunohistochemical localization of cocaine- and amphetamine-regulated transcript peptide in the brain of the catfish, Clarias batrachus (Linn.). J Comp Neuro. 2007;502(2):215–35.

57. Meek J, Nieuwenhuys R. Holosteans and teleosts. In: Nieuwenhuys R, Donkelaar HJ, Nicholson C, editors. The Central Nervous System of Vertebrates: With Posters. Verlag, Berlin: Springer; 1998. p. 759–937.

58. Schwassmann HO, Kruger L. Organization of the visual projection upon the optic tectum of some freshwater fish. J Comp Neuro. 1965;124(1):113–26.

59. Faraco JH, Appelbaum L, Marin W, Gaus SE, Mourrain P, Mignot E. Regulation of hypocretin (orexin) expression in embryonic zebrafish. J Biol Chem. 2006 Oct 6;281(40):29753–61. pmid:16867991

60. Kaslin J, Nystedt JM, Ostergard M, Peitsaro N, Panula P. The orexin/hypocretin system in zebrafish is connected to the aminergic and cholinergic systems. J Neurosci. 2004 Mar 17;24(11):2678–89. pmid:15028760

61. Amiya N, Amano M, Oka Y, Iigo M, Takahashi A, Yamamori K. Immunohistochemical localization of orexin/hypocretin-like immunoreactive peptides and melanin-concentrating hormone in the brain and pituitary of medaka. Neurosci Lett. 2007 10/29/;427(1):16–21. pmid:17935885

62. MacDonald EE, Volkoff H. Molecular cloning and characterization of preproorexin in winter skate (Leucoraja ocellata). Gen Comp Endocrinol. 2010 12/1/;169(3):192–6. pmid:20875823

63. Buckley C, MacDonald EE, Tuziak SM, Volkoff H. Molecular cloning and characterization of two putative appetite regulators in winter flounder (Pleuronectes americanus): Preprothyrotropin-releasing hormone (TRH) and preproorexin (OX). Peptides. 2010 9//;31(9):1737–47. pmid:20685285

64. Huesa G, van den Pol AN, Finger TE. Differential distribution of hypocretin (orexin) and melanin-concentrating hormone in the goldfish brain. J Comp Neuro. 2005;488(4):476–91.

65. Kristensen P, Judge ME, Thim L, Ribel U, Christjansen KN, Wulff BS, et al. Hypothalamic CART is a new anorectic peptide regulated by leptin. Nature. 1998 05/07/print;393(6680):72–6. pmid:9590691

66. Cerdá-Reverter JM, Peter RE. Endogenous melanocortin antagonist in fish: Structure, brain mapping, and regulation by fasting of the goldfish Agouti-related protein gene. Endocrinology. 2003;144(10):4552–61. pmid:12960082

67. Cerdá-Reverter JM, Ringholm A, Schiöth HB, Peter RE. Molecular cloning, pharmacological characterization, and brain mapping of the melanocortin 4 receptor in the goldfish: Involvement in the control of food intake. Endocrinology. 2003;144(6):2336–49. pmid:12746294

68. Gomes AS, Jordal A-EO, Olsen K, Harboe T, Power DM, Rønnestad I. Neuroendocrine control of appetite in Atlantic halibut (Hippoglossus hippoglossus): Changes during metamorphosis and effects of feeding. Comp Biochem Physiol A Physiol. 2015 5//;183(0):116–25.

69. Gomes A, Kamisaka Y, Harboe T, Power D, Rønnestad I. Functional modifications associated with gastrointestinal tract organogenesis during metamorphosis in Atlantic halibut (Hippoglossus hippoglossus). BMC Dev Biol. 2014;14(1):11.